Next: Longitudinal and Torsional Waves
Up: Timoshenko's Beam Equations
Previous: Simulation: Timoshenko's System for
Improved MDKC for Timoshenko's System Via Balancing
In the preceding simulation of the steel beam of rectangular crosssection and linearlyvarying thickness, we have, from (5.24),
; the time step is thus restricted to be quite small. We will now show how the balancing or preconditioning approach applied to the (1+1)D transmission line problem in §3.12 can be used to drastically increase the maximum allowable time step for a given grid spacing.
Suppose we scale the dependent variables according to
and allow the scaling parameters
to be arbitrary smooth functions of . Timoshenko's system (5.17)(5.18) can then be rewritten as

(5.27a) 
where primes above the ,
indicate differentiation.
If we choose:

(5.28) 
Then system (5.26) becomes

(5.29a) 
and the constantproportional terms appear antisymmetrically (note, from (5.27), that
). In the MDKC shown in Figure 5.15, these terms are all interpreted as gyrator couplings, where the gyrator coefficients are spatially varying. It is easy to check that this system is still symmetric hyperbolic according to definition (3.1).
Figure 5.15:
Balanced MDKC for Timoshenko's system.

The MDWD network (not shown) implied by the MDKC will be slightly more difficult to program, because of the additional reflectionfree ports which will necessarily be introduced, but it has the same memory requirements, and the operation count is slightly larger (due chiefly to the postscaling of the MDKC currents which must now be performed in order to obtain the physical dependent variables). We now have, however, that
where
is the maximum group velocity given in (5.20). is now optimal (in the CFL sense), for a constant grid spacing. Referring to the simulation of §5.2.5, it is easy to see that that due to the quartic dependence of the moment of inertia on (for a beam of linearly varying thickness), the maximum time step allowed by the previous approach (also the maximum time step allowed in [131]) will be severely constrained. Using a balanced formulation and MDKC, we now have
. Thus for a given grid spacing, the maximum time step is now 9 times larger. From a practical standpoint, this is a huge computational advantage.
We repeat that balancing is unnecessary if there is no spatial variation in the problem parameters, and that in a region of the material for which the parameters do not vary, we may simply drop the additional gyrator couplings and from the network entirely. We also note that it is possible to incorporate the scaling of the dependent variables into the MDKC itself by introducing transformers with turn ratios
in all the circuit loops; while useful for showing the MDpassivity of the system under scaling, there is no practical reason for doing so.
Next: Longitudinal and Torsional Waves
Up: Timoshenko's Beam Equations
Previous: Simulation: Timoshenko's System for
Stefan Bilbao
20020122