The Matrix Market

Home Search Browse Resources

CYLSHELL: Finite element analysis of cylindrical shells

from the Independent Sets and Generators

Source: Reijo Kouhia, Helsinki University of Technology,
Discipline: Structural mechanics
Accession: Summer 1997

These matrices result from finite element discretization of an octant of a cylindrical shell. The ends of the cylinder are free. A summary of the problem is given below; further details, including plots of the finite element meshes can be found in a separate description.

Notation: T indicates the Radius to thickness ratio, R/t = 10T; MO indicates the model type: DK = Discrete-Kirchhoff, RM = Reissner-Mindlin; EL indicates element type: T3 = 3-node triangular, Q4 = 4-node quadrilateral, etc.; MX indicates the mesh type: M1 = 30x30 uniform, M2 = 100x150 uniform, M3 = graded mesh of 1666 triangles (grading near one boundary and one vertex). The naming convention, cTMOELMX, is based on these.

The finite element is a facet type shell element (3/4 nodes) with drilling rotation incorporated by the Hughes-Brezzi technique and using the penalty parameter value of G/1000 (regularization parameter) , where G is the shear modulus. In order to improve the coarse mesh accuracy the membrane interpolation is amended by the Allman type quadratic modes linked to the drilling rotation. The bending formulation utilizes the stabilized MITC technique with the stabilization parameter equal to 0.4.

For iterative conjugate gradient type solvers, the problem gets harder when the radius to thickness ratio R/t increases. For quadrilateral meshes using the IC(0) preconditioner, iteration counts are about 100 (for R/t=10) and 180 (for R/t=1000) in reaching the relative residual norm of 10-9 depending slightly on the right-hand side vector. For matrices corresponding to triangular meshes the number of iterations doubles in comparison to quadrilateral ones.

The matrices are all symmetric positive definite.

The following estimates of largest and smallest eigenvalues and spectral condition numbers were provided by the author. These were computed using the Lanczos algorithm from QMRPACK (subroutine DSLAL) with 100 steps.

MatrixMax eigenvalueMin eigenvalueCondition No.

The following matrices have been withdrawn from this set: C1RMQ4M1, C2RMT3M1, C3RMQ4M1, C1RMT3M1, C3DKQ4M2, C3RMT3M1, C2RMQ4M1, C3DKT3M2, C3RMT3M3. These matrices are identical to the current matrices with names starting with S, except that the current matrices are represented using 17 digits of precision. Some of the previous matrices were indefinite in the 6 digits in which they were represented.

Matrices in this set:

The Matrix Market is a service of the Mathematical and Computational Sciences Division / Information Technology Laboratory / National Institute of Standards and Technology.

[ Home ] [ Search ] [ Browse ] [ Resources ]

Last change in this page: Wed Sep 22 13:34:26 US/Eastern 2004 [Comments: ]