Shell obstacle course: Scordelis Lo roof
The linear static Scordelis-Lo roof test considers a section of cylindrical shell fixed by rigid diaphragms at it's axial ends subject to a pseudo-gravity distributed load of a magnitude 90. Isotropic material properties are as per the figure below. Due to symmetry, only a quarter of the shell is modeled.
Problem definition [2]
The key result is the vertical displacement of the lateral side at the midpoint, denoted by "u_z" in the diagram above, for which the reference value is "u_z" = 0.3024 [2].
The following Z-displacement contour of the Kratos thin quad element (mesh = 192 elements) is provided for context.
Scordelis-Lo roof results: Z-displacement contour of Kratos thin quad element
The results of the test for the thin and thick triangle Kratos shell elements are presented below.
Scordelis-Lo roof results: triangle elements
The results of the test for the thin and thick quadrilateral Kratos shell elements are presented below.
Scordelis-Lo roof results: quadrilateral elements
Both graphs above indicate all Kratos triangular and quadrilateral shell elements agree with the reference solution.
- Ted Belytschko et al. “Stress projection for membrane and shear locking in shell finite elements”. In: Computer Methods in Applied Mechanics and Engineering 51.1-3 (1985), pp. 221–258.
- Robin Bouclier, Thomas Elguedj, and Alain Combescure. “Efficient isogeometric NURBS-based solid-shell elements: Mixed formulation and method”. In: Computer Methods in Applied Mechanics and Engineering 267 (2013), pp. 86 –110.