I got problem with EQ analysis

[SunnyPark]

hi all.

i have some problems with my code

i want to do EQ analyzis with SAC ground motion of LA21

but there is errors....................

i'm an king begginer.....with EQ analysis....somebody plz help me..

there are my codes (this codes are work properly with push-over analysis)

------------------------------
# As Built
wipe;
model BasicBuilder -ndm 2 -ndf 3;

# nodal coordinates:
node 1 0 0;
node 2 0 400;
node 3 0 1400;

# Single point constraints -- Boundary Conditions
fix 1 1 1 1;


# Define ELEMENTS & SECTIONS -------------------------------------------------------------
set ColSecTag 1001;
set coverCol 40;
set barAreaCol 133;

set IDconcCore 1;
set IDconcCover 2;
set IDreinf 3;
set TorMatID 4;
set IDShear 5;

# MATERIAL
set dam "energy"
set Fy 300;

uniaxialMaterial Pinching4 $IDconcCover 8 0.003 8.3 0.01 8.5 0.015 8 0.18 0.5 0.25 0.05 1.0 1.3 0.3 0.2 0.9 0.5 0.5 2.0 1.5 0.0 1.0 1.9 1.5 0.0 0.0 10 $dam
uniaxialMaterial Pinching4 $IDconcCore 8 0.003 8.3 0.01 8.5 0.015 8 0.18 0.5 0.25 0.05 1.0 1.3 0.3 0.2 0.9 0.5 0.5 2.0 1.5 0.0 1.0 1.9 1.5 0.0 0.0 10 $dam
uniaxialMaterial Hysteretic $IDreinf $Fy 0.003 350 0.02 345 0.3 -$Fy -0.003 -361 -0.02 -345 -0.3 1 0.9 0.01 0.0;

# SECTION
section Fiber 101 {
patch circ $IDconcCore 32 8 0.0 0.0 0.0 160.0 0.0 360.0
patch circ $IDconcCover 32 2 0.0 0.0 160.0 200.0 0.0 360.0
layer circ $IDreinf 16 $barAreaCol 0.0 0.0 160.0 0.0 360.0
}


set ColTransfTag 1;
set ColTransfType Corotational;
geomTransf $ColTransfType $ColTransfTag ;


# element connectivity:
set numIntgrPts 5;
element nonlinearBeamColumn 1 1 2 $numIntgrPts 101 $ColTransfTag;
element nonlinearBeamColumn 2 2 3 $numIntgrPts 101 $ColTransfTag;

puts "column1 assigned"


# define GRAVITY -------------------------------------------------------------
pattern Plain 1 Constant {

load 3 0.0 -1500.0 0.0

}

integrator LoadControl 0 1 0 0

mass 3 0 -0.5 0;

set scale 0.90

set dt 0.02
set g 0.3936
set outFile LA21.g3
set accelSeries "Path -filePath $outFile -dt $dt -factor [expr $scale*$g]"

source ReadSMDFile.tcl
ReadSMDFile LA21.AT2 $outFile dt

pattern UniformExcitation 2 1 -accel $accelSeries
rayleigh 0.0 0.0 0.0 [expr 2*0.045/pow([eigen 1],0.5)]
#rayleigh 0.0 0.0 0.0 [expr 2*0.025/pow(47.60,0.5)]
system BandGeneral
constraints Plain
test EnergyIncr 1.0e-3 60
#test NormDispIncr 1.0e-2 30
algorithm Newton
integrator Newmark 0.5 0.25
numberer RCM
analysis Transient

#Result script
#Record

set lamda [eigen 1]
set omega [expr pow($lamda, 0.5)]
set Tperiod [expr 2*3.14159/$omega]

puts $Tperiod

print node

# Define RECORDERS -------------------------------------------------------------
recorder Node -file MA00ResultGM.out -time -node 3 -dof 1 disp
recorder Element -file MA00Stress01GM.out -time -ele 1 section 1 fiber 0.0 0.0 stressStrain
recorder Element -file MA00Stress02GM.out -time -ele 1 section 1 fiber 0.0 200.0 stressStrain
recorder Element -file MA00Stress03GM.out -time -ele 1 section 1 fiber 0.0 100.0 stressStrain

analyze 3000 0.02

[mklondyn1]

I had also same error. Plz help me out if some one knows how to fix it.

[JaneParker]

I have the same problem! How to fix it?

[SunnyPark]

I fixed it
It was mass problem.
loadConst -time 0.0

I changed my codes like below (after analyzing gravity load).

===========================

puts "Gravity analysis complete"

mass 2 [expr $m/1400*350] 1.e-9 0.0;
mass 3 [expr $m/1400*350] 1.e-9 0.0;
mass 4 [expr $m/1400*350] 1.e-9 0.0;
mass 5 [expr $m/1400*150+$mg] 1.e-9 0.0;

set record_length 15000;
set dt 0.01;
set outFile gm1.g3;
set factor $g;

source ReadSMDFile.tcl
ReadSMDFile gm1.AT2 $outFile dt
#set accelSeries "Path -filePath $outFile -dt $dt -factor 1"
set accelSeries "Series -dt $dt -filePath $outFile -factor [expr $factor]";
pattern UniformExcitation 2 1 -accel $accelSeries

wipeAnalysis

set damp 0.0413
set ww [eigen 2]
set wi [expr sqrt([lindex $ww 0])]
set wj [expr sqrt([lindex $ww 1])]
set alpha [expr $damp*(2*$wi*$wj)/($wi+$wj)]
set beta [expr $damp*(2)/($wi+$wj)]
rayleigh $alpha 0.0 $beta 0.0
#rayleigh 0.0 0.0 0.0 [expr 2*0.05/pow([eigen 1],0.5)]