Results.
| Example | Reference | Records | Largest deviation |
|---|---|---|---|
| Scordelis-Lo roof | Published source | 4 | 4.68% |
| Pinched cylinder with end diaphragms | Published source | 3 | 25.80% |
| Hemisphere with an 18° hole | Published source | 3 | 5.23% |
Problem. The curved shell tests measure whether the element carries membrane and bending action together on curved surfaces. The Scordelis-Lo roof is solved under its own weight, and the pinched cylinder and the hemisphere with a hole under opposite point loads.
Reference document. The reference values are those of MacNeal and Harder (1985) for the Scordelis-Lo roof and the hemisphere, and of Belytschko et al. (1985) for the cylinder.
Numerical model. The curved surfaces are represented piecewise with flat shell elements. The models use symmetry: a quarter of the roof and the hemisphere and an eighth of the cylinder. On the symmetry edges the corresponding translations and rotations are held, and at the ends of the roof and the cylinder the translations in the plane of the section. Each problem is solved on meshes refined from four to thirty-two elements per side.
Scordelis-Lo roof
Published source. MacNeal and Harder (1985), A proposed standard set of problems to test finite element accuracy
| Quantity | Condition | Reference | Computed | Deviation |
|---|---|---|---|---|
| Vertical deflection at the middle of the free edge | mesh = 4×4 | 0.3024 | 0.316549 | 4.68% |
| Vertical deflection at the middle of the free edge | mesh = 8×8 | 0.3024 | 0.304034 | 0.540% |
| Vertical deflection at the middle of the free edge | mesh = 16×16 | 0.3024 | 0.301696 | 0.233% |
| Vertical deflection at the middle of the free edge | mesh = 32×32 | 0.3024 | 0.301377 | 0.338% |
Note. The reference 0.3024 is the value used by MacNeal and Harder. Krysl and Chen (2023) propose 0.302 for shear-flexible (Mindlin) shells, against which the deviation is 0.3%.
Pinched cylinder with end diaphragms
Published source. Belytschko et al. (1985), standard problems for shell elements
| Quantity | Condition | Reference | Computed | Deviation |
|---|---|---|---|---|
| Radial displacement at the load point | mesh = 8×8 | 0.000018248 | 0.0000135354 | 25.80% |
| Radial displacement at the load point | mesh = 16×16 | 0.000018248 | 0.0000169547 | 7.09% |
| Radial displacement at the load point | mesh = 32×32 | 0.000018248 | 0.0000180415 | 1.13% |
Note. The reference values come from thin shell theory. Under a point load the exact solution of shear-flexible shell theory is not finite, so on very fine meshes the result can keep rising slowly; the comparison is meaningful for the standard meshes.
Hemisphere with an 18° hole
Published source. MacNeal and Harder (1985), A proposed standard set of problems to test finite element accuracy
| Quantity | Condition | Reference | Computed | Deviation |
|---|---|---|---|---|
| Radial displacement at the load point | mesh = 8×8 | 0.094 | 0.0890824 | 5.23% |
| Radial displacement at the load point | mesh = 16×16 | 0.094 | 0.0930722 | 0.987% |
| Radial displacement at the load point | mesh = 32×32 | 0.094 | 0.0934263 | 0.610% |
Note. The reference values come from thin shell theory. Under a point load the exact solution of shear-flexible shell theory is not finite, so on very fine meshes the result can keep rising slowly; the comparison is meaningful for the standard meshes.