StructuralMind

Standard tests, curved shells

Updated:

On this page

Results.

ExampleReferenceRecordsLargest deviation
Scordelis-Lo roofPublished source44.68%
Pinched cylinder with end diaphragmsPublished source325.80%
Hemisphere with an 18° holePublished source35.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

QuantityConditionReferenceComputedDeviation
Vertical deflection at the middle of the free edgemesh = 4×40.30240.3165494.68%
Vertical deflection at the middle of the free edgemesh = 8×80.30240.3040340.540%
Vertical deflection at the middle of the free edgemesh = 16×160.30240.3016960.233%
Vertical deflection at the middle of the free edgemesh = 32×320.30240.3013770.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

QuantityConditionReferenceComputedDeviation
Radial displacement at the load pointmesh = 8×80.0000182480.000013535425.80%
Radial displacement at the load pointmesh = 16×160.0000182480.00001695477.09%
Radial displacement at the load pointmesh = 32×320.0000182480.00001804151.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

QuantityConditionReferenceComputedDeviation
Radial displacement at the load pointmesh = 8×80.0940.08908245.23%
Radial displacement at the load pointmesh = 16×160.0940.09307220.987%
Radial displacement at the load pointmesh = 32×320.0940.09342630.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.