Group Cylshell

Group Description
Finite element analysis of cylindrical shells, R. Kouhia, Helsinki.

Set CYLSHELL

Source:	    Reijo Kouhia, Helsinki University of Technology, reijo.kouhia at hut.fi

Discipline:	Structural mechanics
Accession (in Matrix Market collection):  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.

Matrix	Max eigenvalue	Min eigenvalue	Condition No.
S1RMQ4M1    6.8743215039276076E+05  3.7969627252353000E-01  1.810479E+06
S2RMQ4M1    6.8743398166879822E+04  3.8746327712024738E-04  1.774191E+08
S3RMQ4M1    6.8745248885112433E+03  3.8936812501755643E-07  1.765559E+10
S1RMT3M1    9.6684220885668334E+05  3.7976650205204426E-01  2.545886E+06
S2RMT3M1    9.6681138121311058E+04  3.8747124280659988E-04  2.495182E+08
S3RMT3M1    9.6688590706791765E+03  3.8956266582265278E-07  2.481978E+10
S3DKQ4M2    4.6016534362482462E+03  2.4268618877322254E-08  1.896133E+11
S3DKT3M2    8.7984363691286471E+03  2.4269393992599454E-08  3.625322E+11
S3RMT3M3    9.5986080894852857E+03  3.9983547883056739E-07  2.400639E+10
Displaying all 7 collection matrices
Id Name Group Rows Cols Nonzeros Kind Date Download File
1605 s1rmq4m1 Cylshell 5,489 5,489 262,411 Structural Problem 1997 MATLAB Rutherford Boeing Matrix Market
1608 s1rmt3m1 Cylshell 5,489 5,489 217,651 Structural Problem 1997 MATLAB Rutherford Boeing Matrix Market
1606 s2rmq4m1 Cylshell 5,489 5,489 263,351 Structural Problem 1997 MATLAB Rutherford Boeing Matrix Market
1609 s2rmt3m1 Cylshell 5,489 5,489 217,681 Structural Problem 1997 MATLAB Rutherford Boeing Matrix Market
1607 s3rmq4m1 Cylshell 5,489 5,489 262,943 Structural Problem 1997 MATLAB Rutherford Boeing Matrix Market
1610 s3rmt3m1 Cylshell 5,489 5,489 217,669 Structural Problem 1997 MATLAB Rutherford Boeing Matrix Market
1611 s3rmt3m3 Cylshell 5,357 5,357 207,123 Structural Problem 1997 MATLAB Rutherford Boeing Matrix Market