The Naghdi-Cooper System II Formulation

Another reason for choosing this particular shell model is that it can be simply written as Mindlin's plate system (in coordinates and instead of and ) coupled with the membrane shell. We can write this system in the form of (3.1), where the dependent variable is defined by

Here, the first 8 variables are precisely those that appear in the Mindlin theory (see §5.4), but in cylindrical coordinates; we have changed subscripts to , and written the instead of for the radial transverse velocity. The other five variables appear in subsystem (5.47) of the membrane shell theory (see §5.5.1). The system matrices are

where , , , , and are as defined in §5.4, and appears in §5.5.1. The coupling matrix is

Note that this coupling disappears in the limit as the shell radius becomes large (effectively leaving us with with the Mindlin system). is indeed anti-symmetric, so we are guaranteed a lossless MDKC network; in fact, this network can be directly constructed from the two networks shown in Figures 5.16 and 5.23, by attaching terminals and in the former to and in the latter. The scaling parameters , and can be chosen optimally according to (5.37), (5.40) and (5.41), and and can be set as in (5.48) and (5.49), giving a bound for passivity on the combined network,

where is the bounding space step/time step ratio for the Mindlin network, from (5.43) (in cylindrical coordinates), and is the same quantity for the membrane shell system, from (5.50).

Because this system can be constructed entirely by connecting subnetworks that we have already examined in detail, it seems unnecessary to show the discrete MDWD network or the alternate MDKC and its discrete form suitable for DWN implementation. The MDWD network will be exactly the combination of the ``Mindlin'' system, shown at bottom in Figure 5.16, and the MDWD network corresponding to the MDKC for the membrane shell, shown in Figure 5.23; recall that the MDKC for the membrane shell system (with a free port with terminals and ) is identical in form to that of the uncoupled five-variable Mindlin subsystem, and thus its MDWD counterpart will be of the same form as well.

For the transformed network to be used to generate a DWN, a few comments are in order. For the Mindlin system, we first applied network theoretic rules in order to arrive at a modified form, shown at top in Figure 5.17. In this case, the transverse velocity (renamed in this section) and the bending moments , and (renamed , and ) have been interpreted as voltages instead of currents. This transformed network can then be connected (via terminals and in Figure 5.17) to a transformed form of the membrane shell system, shown in Figure 5.23; for the shell subsystem, , and will be treated as voltages, and and as currents.

After connecting these transformed subnetworks, and applying the usual alternative discretization rules, we end up with a DWN that will operate on an interleaved grid as shown in Figure 5.24. Grids (a) and (b) are precisely the Mindlin grid shown in Figure 5.18, and are coupled instantaneously to a third grid (c), which adds the effect of curvature to the system.