Thick Shell Theory (Degenerate Solid)
Thick Shell Theory (Degenerate Solid): Theoretical Foundations
What is a Degenerated Solid Shell?
Professor, what is a "degenerated solid"? Is it a hybrid of shell and solid elements?
Exactly. The Degenerated Solid Shell Element is created by degenerating (reducing) the degrees of freedom in the thickness direction of a 3D solid element into a shell element. It was proposed by Ahmad-Irons-Zienkiewicz (1970).
The concept is simple:
1. Start with a 3D solid element having nodes on its top and bottom surfaces.
2. Use the displacement of the mid-surface and the rotation angle about the normal direction as degrees of freedom.
3. Assume a linear displacement distribution in the thickness direction (Mindlin assumption).
Starting from 3D and adding assumptions to make it 2D. That's the opposite approach of conventional shell theory.
Yes. Conventional shell theory starts from 2D equations, but the degenerated solid approach starts from 3D and "degenerates the unused degrees of freedom." The result converges to the same Mindlin shell, but the implementation is simpler because it's based on 3D.
Handling Thick Shells
What is meant by a "thick shell"?
Shells with moderate thickness where $R/t$ is around 10 to 30. It's the intermediate region, neither thin ($R/t > 30$) nor thick ($R/t < 10$, essentially solid).
For thick shells:
- Shear deformation cannot be ignored.
- The through-thickness stress $\sigma_z$ is not exactly zero.
- Membrane-bending coupling is strong.
Can Mindlin shell elements handle this?
They can handle shear deformation, but $\sigma_z \neq 0$ cannot be treated by conventional shell elements. To handle this, you need solid shell elements (shell elements with displacement degrees of freedom in the thickness direction) or solid elements.
Solid Shell Elements
What kind of elements are "solid shells"?
They look like solid elements (HEX8 or HEX20), but their internal formulation is optimized for shells.
| Element | Solver | Features |
|---|---|---|
| SC8R | Abaqus | 8-node solid shell. Reduced integration + anti-locking measures |
| SOLSH190 | Ansys | Solid shell. Bending representation with one element in thickness direction |
| CHEXA(solid-shell) | LS-DYNA | LSDYNA implementation of solid shell |
What are the advantages of solid shells?
Summary
Let me organize the thick shell theory.
Key points:
- Degenerated solid — Creating a shell by degenerating the thickness direction from a 3D solid
- Intermediate region $R/t = 10 \sim 30$ — Neither thin nor solid
- Solid shell elements — Look like solids, behave like shells. Handle contact surfaces and $\sigma_z$
- SC8R (Abaqus), SOLSH190 (Ansys) — Representative elements
- Bending representation with one element in thickness direction — Efficient
So we choose: thin → Mindlin shell, intermediate → solid shell, thick → solid element.
Judgment based on $R/t$ is fundamental. If in doubt, solve with both and compare the results.
Mindlin-Reissner Thick Shell Theory
The Mindlin-Reissner theory, which forms the basis for thick shells, was formulated independently by Raymond Mindlin and Eric Reissner between 1945 and 1951. Unlike the Kirchhoff assumption, it allows "normals to tilt due to shear deformation," explicitly treating transverse shear strains εxz and εyz as degrees of freedom. This extends the applicable range to plates with thickness/span ratios up to about 1/5, making it an essential theory for analyzing interlaminar shear in composite laminates.
Computational Methods for Thick Shell Theory (Degenerate Solid)
Implementation of Solid Shells
Please tell me about implementation considerations for solid shells.
Solid shells are "thin" solid elements, so the aspect ratio in the thickness direction becomes very large. For regular solid elements, accuracy degrades at aspect ratio > 5, but solid shells internally correct for this.
Locking Countermeasures
Locking issues in solid shells:
1. Shear locking — Occurs in thin plate bending. Countered by ANS method.
2. Volumetric locking — Occurs with incompressible materials. Countered by EAS method or B-bar method.
3. Trapezoidal locking — When the element tapers in the thickness direction. Unique to solid shells.
4. Curvature thickness locking — When elements become trapezoidal in the thickness direction on curved surfaces.
Is trapezoidal locking unique to solid shells?
Yes. When there is taper (trapezoidal shape) in the thickness direction, regular solid elements cannot correctly represent bending. Solid shell elements eliminate trapezoidal locking using EAS (Enhanced Assumed Strain).
Usage Tips
Points to note when using solid shell elements:
- One element in the thickness direction — Two or more are unnecessary (design philosophy of solid shells).
- Correctly specify the element's "thickness direction" — Specify stack direction in Abaqus's *SOLID SECTION.
- Meshing curved surfaces — Mesh the top and bottom CAD surfaces separately, then connect them in the thickness direction.
It's efficient that one element in the thickness direction is enough. With regular solid HEX8, you needed at least 4 elements in the thickness direction.
That's precisely the biggest advantage of solid shells. They achieve bending accuracy equivalent to regular shell elements with just one HEX8-equivalent element in the thickness direction. The DOF count is similar to shell elements, but they offer benefits in contact and thickness variation.
Summary
Let me organize the numerical methods for solid shells.
Key points:
- Represent bending with one element in the thickness direction — Efficient.
- Locking countermeasures are essential — Combination of ANS + EAS.
- Trapezoidal locking is unique to solid shells — Countered by EAS method.
- Specifying stack direction is important — Correctly define the thickness direction.
- Optimal for problems with contact surfaces on both sides — An advantage not found in regular shell elements.
MITC Element's Shear Locking Countermeasure
The MITC (Mixed Interpolation of Tensorial Components) method is a thick shell locking countermeasure technique developed by Bathe and Dvorkin at MIT in 1986. By independently interpolating shear strains, it ensures uniform accuracy from thin to thick plates. MITC4 corresponds to 4-node shells, MITC9 to 9-node shells, demonstrating high performance with displacement errors within 5% even for thickness/span ratios of 1/1000.
Experience the theory firsthand with the interactive simulator for this field
All Simulators