Error in Analysis

[MBCA]

Hi everyone!
I am getting this error message when I run the analysis. The eigen Value analysis gets completed but after that it throws this error:
!!!Pfect():diagonal 3 is Zero!!!
In SymSparsein Solver error in factorization.

Please can someone have a look at the code and suggest some possible solution. Thank you!

##It is one bay frame with 2 horizontal beams
wipe all;
model BasicBuilder -ndm 3 -ndf 6
##Information about the isolator##
set FyExt 18300; #Yield Strength
set alphaExt 0.1; #Post Yield Stiffness Ratio
set GExt 0.4; #Shear Modulus of Elastomeric Bearing
set KExt 1200; #Bulk Modulus of Rubber
set D1Ext 50; #Internal Diameter
set D2Ext 480; #Outer Diameter excluding the cover thickness
set tsExt 4; #Single Steel Shim layer Thickness
set trExt 5; #Single Rubber Layer Thickness
set nExt 25; #Number of Rubber Layers
##NODE DEFINITION##
node 1 0 0 0; fix 1 1 1 1 1 1 1 ;
node 2 6000 0 0 ; fix 2 1 1 1 1 1 1 ;
node 3 0 200 0 -mass 10 10 0 0 0 0;
node 4 6000 200 0 -mass 10 10 0 0 0 0;
node 5 0 4000 0 -mass 10 10 0 0 0 0;
node 6 6000 4000 0 -mass 10 10 0 0 0 0;
##MATERIAL DEFINITION##
set Fy 250.0;
set E 200000.0;
set alpha 0.01;
set numIntg 4;
uniaxialMaterial Steel02 1 $Fy $E $alpha 18.0 0.95 0.4;
##SECTION d w b t
section WFSection2d 2 1 525.0 8.9 165.0 11.4 4 2;
geomTransf Linear 1 0 0 1 ;
geomTransf Linear 2 0 0 1 ;
##BEAMS##
element dispBeamColumn 3 3 4 $numIntg 2 1;
element dispBeamColumn 4 5 6 $numIntg 2 1;

##COLUMNS##
element dispBeamColumn 1 3 5 $numIntg 2 2;
element dispBeamColumn 2 4 6 $numIntg 2 2;

##ISOLATORS##
element LeadRubberX 5 1 3 $FyExt $alphaExt $GExt $KExt $D1Ext $D2Ext $tsExt $trExt $nExt 0 1 0 1 0 0 10 0.5 1 0.5 0 0 0 1.12e-8 130e6 50 14.1 0 0 0 0 0
element LeadRubberX 6 2 4 $FyExt $alphaExt $GExt $KExt $D1Ext $D2Ext $tsExt $trExt $nExt 0 1 0 1 0 0 10 0.5 1 0.5 0 0 0 1.12e-8 130e6 50 14.1 0 0 0 0 0

timeSeries Path 2 -dt 0.02 -filePath elcentro.txt -factor 9810

##RECORDERS##
file mkdir output3
recorder Node -file output3/acceleration(Absolute).txt -time -timeSeries 2 -node 1 3 4 2 -dof 1 accel;

##EIGEN VALUE ANALYSIS##
puts "Eigen Value Analysis Begins"
set eigenvalues [eigen 2];
set lambda1 [lindex $eigenvalues 0]
set lambda2 [lindex $eigenvalues 1]
set w1 [expr pow ($lambda1, 0.5)]
set w2 [expr pow ($lambda2, 0.5)]
set T1 [expr 2*3.14/$w1];
set T2 [expr 2*3.14/$w2];
puts "T1 is $T1";
puts "T2 is $T2";

set dampratio 0.05;
rayleigh 0 0 [expr 2*$dampratio/$w1] 0;

##EQ ANALYSIS##
pattern UniformExcitation 3 1 -accel 2;
system SparseSYM;
constraints Plain;
numberer RCM;
algorithm Linear;
#algorithm NewtonLineSearch;
integrator Newmark 0.5 0.25;
analysis Transient ;
analyze 1561 0.02;
puts "Dynamic Analysis is Complete";

[fmk]

if you look at the output you will see that although it computed the eigenvalue it was not happy with the matrices .. the SparseSYM message confirms this .. as the solver fails in the initial step the problem is obviously in the model definition.

[MBCA]

Thank You! You were right. The Model was not restrained out of plane i.e why it was displaying high value of Time Period.