FEM, cylindrical shell, 150x100 quad. mesh, R/t=1000
Name s3dkq4m2
Group GHS_psdef
Matrix ID 1275
Num Rows 90,449
Num Cols 90,449
Nonzeros 4,427,725
Pattern Entries 4,820,891
Kind Structural Problem
Symmetric Yes
Date 1997
Author R. Kouhia
Editor R. Boisvert, R. Pozo, K. Remington, B. Miller, R. Lipman, R. Barrett, J. Dongarra
Structural Rank 90,449
Structural Rank Full true
Num Dmperm Blocks 1
Strongly Connect Components 1
Num Explicit Zeros 393,166
Pattern Symmetry 100%
Numeric Symmetry 100%
Cholesky Candidate yes
Positive Definite yes
Type real
Download MATLAB Rutherford Boeing Matrix Market
%FILE  s3dkq4m2.mtx                                                            
%TITLE Cyl shell R/t=1000 unif 150x100 quad mesh DKQ elem with drill rot       
%KEY   s3dkq4m2                                                                
%CONTRIBUTOR Reijo Kouhia (                                
%REFERENCE   M. Benzi, R. Kouhia, M.Tuma: An assesment of some                 
%            preconditioning techniques in shell problems                      
%            Technical Report LA-UR-97-3892, Los Alamos National Laboratory    
%BEGIN DESCRIPTION                                                             
% Matrix from a static analysis of a cylindrical shell                         
% Radius to thickness ratio R/t = 1000                                         
% Length to radius ratio    R/L = 1                                            
% One octant discretized with uniform 150 x 100 quadrilateral mesh             
% element:                                                                     
% facet-type shell element where the bending part is formulated                
% using the discrete Kirchhoff quadrilateral (Lyons-Crisfield version),        
% the membrane part includes drilling rotations using                          
% the Hughes-Brezzi formulation with (regularizing parameter = G/1000,         
% where G is the shear modulus)                                                
% full 2x2 Gauss-Legendre integration                                          
% --------------------------------------------------------------------------   
% Note:                                                                        
% The sparsity pattern of the matrix is determined from the element            
% connectivity data assuming that the element matrix is full.                  
% Since this case the  material model is linear isotropically elastic          
% and the FE mesh is  uniform there exist some zeros.                          
% Since the removal of those zero elements is trivial                          
% but the reconstruction of the current sparsity                               
% pattern is impossible from the sparsified structure without any further      
% knowledge of the element connectivity, the zeros are retained in this file.  
% ---------------------------------------------------------------------------  
%END DESCRIPTION