Next: Boundary Conditions
Up: Elastic Solids
Previous: Phase and Group Velocities
Scattering Networks for the Navier System
Nitsche has represented this system as an MDKC in [131]. Choosing current-like variables
he derived an MDKC for the Navier system, which we have reproduced (with some minor changes) at top in Figure 5.25, where coordinates defined by the transformation matrix (3.24) have been used. Notice that again, because
is not diagonal, the time derivatives of the components
,
and
are coupled. In the circuit representation, this is represented by a three-port coupled inductance of the form of (2.48), with a matrix inductance
defined by
where the inductances
,
,
and
appear in the figure.
The positivity condition on the inductances
,
,
,
,
and
, as well as a positive definiteness condition on the matrix
gives, as optimal choices of the parameters
and
,
and a bound on
,
When the material parameters are constant,
reduces to
A modified MDKC is shown at top in Figure 5.26, where velocities are treated as currents, and stresses as voltages; as such, the coupled inductance in Figure 5.25 has become a coupled capacitance. This network may be discretized in a way very similar to the network for Maxwell's Equations, as described in §4.10.6. Under the spectral mappings defined by (4.112) (with step sizes of
,
), the connecting LSI two-ports decompose into series/parallel connections of multidimensional unit elements. The one-port inductances and capacitances, as well as the coupled capacitance are discretized using the trapezoid rule, with a step size of
. The resulting multidimensional DWN is shown at bottom in Figure 5.26, and the interleaved computational grid in Figure 5.27 (which is very similar to the grid for the DWN for Maxwell's Equations, as shown in Figure 4.55).
Next: Boundary Conditions
Up: Elastic Solids
Previous: Phase and Group Velocities
Stefan Bilbao
2002-01-22