Next: Elastic Solids Up: Cylindrical Shells Previous: The Membrane Shell

The Naghdi-Cooper System II Formulation

As we mentioned, the membrane model of the cylindrical shell neglects certain important effects, in particular the crucial transverse shear effects. Many so-called higher-order shell theories have appeared in the literature; for a good survey of these theories, we refer to [77,81]. We have decided to focus on the Naghdi-Cooper system II shell model [31,77,128] because it can be simply interpreted as a passive circuit. We note in passing that not all shell theories have this property--Mirsky-Herrmann theory [77], for example, does not, and Naghdi and Cooper's system II, for example, is simplified from their proposed system I which also does not. The problem, more specifically, is that these systems can apparently not be written in the special symmetric hyperbolic form of (3.1), for which the matrices and will be independent of and . This property is essential here in that such dependence considerably complicates the inter-loop coupling in an MDKC (notice that the port-resistances of the Jaumann two-ports which realize this coupling been constant for every system we have looked at so far).

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

 (5.53)

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.

Next: Elastic Solids Up: Cylindrical Shells Previous: The Membrane Shell
Stefan Bilbao 2002-01-22