The OpenSees Community

hi
i have decided to modeling and calculating euler buckling load factor for column with one end restrained and other end free
how to write this code and how to record euler buckling load factor?
thank you

Dear Frank, I create code but dont work and more error

model BasicBuilder -ndm 2 -ndf 3
file mkdir Kbrace
foreach L {3.0 3.5 4.0 4.5 5.0 5.5 6.0} H {3.0 3.5 4.0 4.5 5.0} l {[expr 0.1*L] [expr 0.2*L] [expr 0.3*L] [expr 0.4*L] [expr 0.5*L]} h {[expr 0.1*H] [expr 0.2*H] [expr 0.3*H] [expr 0.4*H] [expr 0.5*H]} P {[expr 0.5*p] [expr 0.75*p] [expr 1.0*p] [expr 1.25*p] [expr 1.5*p]} {
wipe
node 1 0. 0.
node 2 $L 0.
node 3 0. $H
node 4 $L $H
node 5 0. 0.
node 6 [expr ($L-($l*0.5))] [expr ($H-($h*0.5))]
node 7 [expr ($L-$l)] $H
node 8 $L [expr ($H-$h)]
node 9 [expr ($L-($l*0.5))] [expr ($H-($h*0.5))]
node 10 [expr ($L-$l)] $H
node 11 $L [expr ($H-$h)]
equalDOF 1 5 1 2
equalDOF 6 9 1 2
equalDOF 7 10 1 2
equalDOF 8 11 1 2
equalDOF 3 4 2 3
fix 1 1 1 1
fix 2 1 1 1
set M [expr 1000]
mass 3 [expr $M*L*0.5] 1e-100 1e-100
mass 4 [expr $M*L*0.5] 1e-100 1e-100
set p [expr $M*L*0.5*0.16]
geomTransf Linear 1
set Abc 1.6e-3
set Abrace 1.6e-3
set E 2.1e10
set I 4.92e-6
set Ibrace 4.92e-6
element elasticBeamColumn 1 1 3 $Abc $E $I 1
element elasticBeamColumn 2 2 4 $Abc $E $I 1
element elasticBeamColumn 3 3 4 $Abc $E $I 1
element elasticBeamColumn 4 5 6 $Abrace $E $Ibrace 1
element elasticBeamColumn 5 7 8 $Abrace $E $Ibrace 1
pattern Plain 1 Linear {
load 3 $P 0. 0.
}
pattern Plain 2 Linear {
load 4 $P 0. 0.
}
set pi [expr 2.0*asin(1.0)]
set nEigenI 1
set nEigenJ 2
set lambdaN [eigen [expr $nEigenJ]]
set lambdaI [lindex $lambdaN [expr 0]]
set lambdaJ [lindex $lambdaN [expr $nEigenJ-1]]
set W1 [expr pow($lambdaI,0.5)]
set W2 [expr pow($lambdaJ,0.5)]
set T1 [expr 2.0*$pi/$W1]
set T2 [expr 2.0*$pi/$W2]
puts "T1 = $T1 s"
puts "T2 = $T2 s"
pattern Plain 3 Linear {
load 3 10. 0. 0.
}
recorder display "Kbrace" 500 0 500 500 -wipe
prp 0 0 1
vup 0 1 0
vpn 0 0 1
viewWindow -5 5 -5 5
display 1 5 1
constraints Plain
numberer Plain
system BandGeneral
test EnergyIncr 1.e-10 1000
algorithm KrylovNewton
integrator DisplacementControl 4 1 0.0001
analysis Static
analyze 1000
loadConst -time 0.0
}

help me please

thank you dear frank
bests for you

hi friend
im modeling 2D one story one bay steel brace frame in opensees.
L is length of beam
H is Height of column
L values are 3m,3.5m,4m,4.5m,5m,5.5m,6m
H values are 2m,2.5m,3m,3.5m,4m,4.5m,5m
section of beam,columns and brace are 10 sections with different values for each elements
pushover load apply to model with 10 factor. 0.25p,0.5p,0.75p,1p,1.25p,1.5p,1.75p,2p,2.25p,2.5p
i want to do analysis one model with all above parameter values successive. for example i create one model in opensees with one parameter values.after finishing analysis opensees automatic start another analysis with the next parameter and after finishing this start another analysis with the next parameter and this process perform automatic until all parameter used in analysis.
analysis number1. L=3m H=2m , p=0.25
analysis number2. L=3.5m H=2.5m , p=0.25
analysis number3. L=3.5m H=2m , p=0.5
analysis number4. L=4m H=3m , p=1
.
.
.
.
.
.
.
.
analysis number99. L=6m H=4.5m , p=2.5
analysis number100. L=6m H=5m , p=2.5
can opensees do this analysis with one model?
how i assigne multiple values for length and height and so all?
thank you

hi
i have 2D-2story model with 2 different damping (C) for each story- [C1 & C2]
i assign C1 and C2 to each story with region command.(C=Zeta*2*m*w)
i decided write active algorithm for change damping during analysis between C1 and C2. (LQR or PID algorithm)
for example: when t=0.1s damping for story 1 is C1 and damping for story 2 is C2 and when t=0.5s damping for story 1 is C2 and damping for story 2 is C1. and this process repeat continuously for each time during the time history analysis.
opensees can do this active algorithm? or i force to use matlab for this purpose?
please help me
thank you

I want to apply 2 different damping ratio for my model by use region commmand. I wrote a code but it did'nt work correctly , please tell me my mistake?

wipe

model BasicBuilder -ndm 2 -ndf 3

file mkdir nonclasic

set L [expr 5]
set H [expr 5]

node 1 0 0
node 2 $L 0
node 3 0 $H
node 4 $L $H
node 5 0 [expr 2*$H]
node 6 $L [expr 2*$H]

fix 1 1 1 1
fix 2 1 1 1

equalDOF 3 4 1
equalDOF 5 6 1

set M1 [expr 1000]
set M2 [expr 1000]

mass 3 $M1 0. 1e-10
mass 4 $M1 0. 1e-10
mass 5 $M2 0. 1e-10
mass 6 $M2 0. 1e-10

geomTransf Linear 1

set E1 [expr 2.1e10]
set E2 [expr 2.1e10]

set Fy1 [expr 2.4e7]
set Fy2 [expr 2.4e7]

uniaxialMaterial Steel01 1 $Fy1 $E1 0.02

section fiberSec 1 {
patch quad 1 2 4 -0.1 -0.1 0.1 -0.1 0.1 -0.09 -0.1 -0.09
patch quad 1 2 8 -0.005 -0.09 0.005 -0.09 0.005 0.09 -0.005 0.09
patch quad 1 2 4 -0.1 0.09 0.1 0.09 0.1 0.1 -0.1 0.1
}

uniaxialMaterial Steel01 2 $Fy2 $E2 0.02

section fiberSec 2 {
patch quad 2 2 4 -0.1 -0.1 0.1 -0.1 0.1 -0.09 -0.1 -0.09
patch quad 2 2 8 -0.005 -0.09 0.005 -0.09 0.005 0.09 -0.005 0.09
patch quad 2 2 4 -0.1 0.09 0.1 0.09 0.1 0.1 -0.1 0.1
}

#set I1 [expr 4.1e-5]
#set I2 [expr 4.1e-5]
#set K1 [expr (2*12*$E1*$I1)/pow($H,3)]
#set K2 [expr (2*12*$E2*$I2)/pow($H,3)]
#set X1 [expr 0.15]
#set X2 [expr 0.1]
#set K [expr $K1+$K2]
#set M [expr $M1+$M2]
#set W [expr sqrt($K/$M)]
#set C1 [expr $X1*(2*$M*$W)]
#set C2 [expr $X2*(2*$M*$W)]

element nonlinearBeamColumn 1 1 3 10 1 1
element nonlinearBeamColumn 2 2 4 10 1 1
element nonlinearBeamColumn 3 3 4 10 1 1
element nonlinearBeamColumn 4 3 5 10 2 1
element nonlinearBeamColumn 5 4 6 10 2 1
element nonlinearBeamColumn 6 5 6 10 2 1

set accel "Series -dt 0.01 -filePath ElCentroArray9180(dt=0.01)PGA=0.4g.txt -factor [expr 9.81]"

pattern UniformExcitation 3 1 -accel $accel

set pi [expr 2.0*asin(1.0)]
set nEigenI 1
set nEigenJ 2
set lambdaN [eigen [expr $nEigenJ]]
set lambdaI [lindex $lambdaN [expr 0]]
set lambdaJ [lindex $lambdaN [expr $nEigenJ-1]]
set w1 [expr pow($lambdaI,0.5)]
set w2 [expr pow($lambdaJ,0.5)]
set T1 [expr 2.0*$pi/$w1]
set T2 [expr 2.0*$pi/$w2]
puts "T1 = $T1 s"
puts "T2 = $T2 s"

set n1 3

set zeta1 10.01
set a01 [expr $zeta1*2.0*$w1*$w2/($w1 + $w2)]
set a11 [expr $zeta1*2.0/($w1 + $w2)]
set a11_mod [expr $a11*(1.0+$n1)/$n1]

region 1 -eleRange 1 3 -rayleigh 0.0 0.0 $a11_mod 0.0
region 2 -node 1 2 3 4 -rayleigh $a01 0.0 0.0 0.0

set n2 3

set zeta2 0.01
set a02 [expr $zeta2*2.0*$w1*$w2/($w1 + $w2)]
set a12 [expr $zeta2*2.0/($w1 + $w2)]
set a12_mod [expr $a12*(1.0+$n2)/$n2]

region 3 -eleRange 4 6 -rayleigh 0.0 0.0 $a12_mod 0.0
region 4 -node 5 6 -rayleigh $a02 0.0 0.0 0.0

recorder Node -file nonclasic/node4disp.out -time -node 4 -dof 1 disp
recorder Node -file nonclasic/node6disp.out -time -node 6 -dof 1 disp
recorder Node -file nonclasic/RBase1.out -time -node 1 -dof 1 2 3 reaction
recorder Node -file nonclasic/RBase2.out -time -node 2 -dof 1 2 3 reaction
recorder Drift -file nonclasic/Drift1.out -time -iNode 2 -jNode 4 -dof 1 -perpDirn 2
recorder Drift -file nonclasic/Drift2.out -time -iNode 4 -jNode 6 -dof 1 -perpDirn 2
recorder Element -file nonclasic/FCol1.out -time -ele 1 globalForce;
recorder Element -file nonclasic/FCol2.out -time -ele 2 globalForce;
recorder Element -file nonclasic/FCol4.out -time -ele 4 globalForce;
recorder Element -file nonclasic/FCol5.out -time -ele 5 globalForce;
recorder Element -file nonclasic/FBeam3.out -time -ele 3 globalForce;
recorder Element -file nonclasic/FBeam6.out -time -ele 6 globalForce;

recorder display "nonclasic" 500 0 500 500 -wipe
prp 0 0 1
vup 0 1 0
vpn 0 0 1
viewWindow -10 10 -10 10
display 1 5 5

constraints Plain
numberer Plain
system BandGeneral
test NormDispIncr 1.0e-10 10
algorithm Newton
integrator Newmark 0.5 0.25
analyze 5000 0.01

set a [eigen 2]
set eigen1 [lindex $a 0]
set eigen2 [lindex $a 1]
puts "eigens are $eigen1,$eigen2"
set Period1 [expr 2*3.1416/sqrt($eigen1)]
set Period2 [expr 2*3.1416/sqrt($eigen2)]
puts "Periods are $Period1,$Period2"

set numModes 2
for { set k 1 } { $k <= $numModes } { incr k } {
recorder Node -file nonclasic/mode$k.out -nodeRange 1 6 -dof 1 2 3 "eigen $k"
}

set lambda [eigen $numModes];
set omega {}
set f {}
set T {}
set pi 3.141593

foreach lam $lambda {
lappend omega [expr sqrt($lam)]
lappend f [expr sqrt($lam)/(2*$pi)]
lappend T [expr (2*$pi)/sqrt($lam)]
}

set period "Periods.txt"
set Periods [open $period "w"]
foreach t $T {
puts $Periods " $t"
}
close $Periods


#Thank You

thank you frank

if nod4 disp > 0.01
assign C1=0.01 for region 1 (region1 -ele 1 2 3 4)
if node6 disp > 0.05
assign C2=0.2 for region 2 (region2 -ele 5 6)
if Drift1<0.005
assign Copt=0.1 for all

how to write and apply this algoritm to my opensees code?

Thank you frank
i am decided to apply 3 damping ratio to model.(active control problem-TTSOONG).for specifiec time apply specified damping ratio to model. damping ratio are: Cmin,Copt,Cmax. in other words 3 damping ratio can switching together. for this purpose i apply some criteria to model for switching damping together like a displacement node 4 more than certain content. i write this code but nnot correct properly

for {set i 1} {$i <= 5000} {incr i 1} {
analyze 5000 0.01
if {[expr $Dispnode4 1] > 0.01} {
puts "[getTime] [nodeDisp 4 1]"
set x3Damp 1;
set M3propSwitch 1.0;
set K3currSwitch 0.0;
set K3commSwitch 1.0;
set K3initSwitch 0.0;
set nEigenI 1;
set nEigenJ 2;
set lambdaN [eigen [expr $nEigenJ]];
set lambdaI [lindex $lambdaN [expr $nEigenI-1]];
set lambdaJ [lindex $lambdaN [expr $nEigenJ-1]];
set omegaI [expr pow($lambdaI,0.5)];
set omegaJ [expr pow($lambdaJ,0.5)];
set alphaM3 [expr $M3propSwitch*$x3Damp*(2*$omegaI*$omegaJ)/($omegaI+$omegaJ)];
set betaK3curr [expr $K3currSwitch*2.*$x3Damp/($omegaI+$omegaJ)];
set betaK3comm [expr $K3commSwitch*2.*$x3Damp/($omegaI+$omegaJ)];
set betaK3init [expr $K3initSwitch*2.*$x3Damp/($omegaI+$omegaJ)];
rayleigh $alphaM3 $betaK3curr $betaK3init $betaK3comm;
region 1 -eleRange 1 3 rayleigh $alphaM3 $betaK3curr $betaK3init $betaK3comm;
region 3 -nodeRange 1 4 rayleigh $alphaM3 $betaK3curr $betaK3init $betaK3comm;
} if {[expr abs([nodeDisp 6 1])]>0.001} {
puts "[getTime] [nodeDisp 6 1]"
set x4Damp 2;
set M2propSwitch 1.0;
set K4currSwitch 0.0;
set K4commSwitch 1.0;
set K4initSwitch 0.0;
set nEigenI 1;
set nEigenJ 2;
set lambdaN [eigen [expr $nEigenJ]];
set lambdaI [lindex $lambdaN [expr $nEigenI-1]];
set lambdaJ [lindex $lambdaN [expr $nEigenJ-1]];
set omegaI [expr pow($lambdaI,0.5)];
set omegaJ [expr pow($lambdaJ,0.5)];
set alphaM4 [expr $M4propSwitch*$x4Damp*(2*$omegaI*$omegaJ)/($omegaI+$omegaJ)];
set betaK4curr [expr $K4currSwitch*2.*$x4Damp/($omegaI+$omegaJ)];
set betaK4comm [expr $K4commSwitch*2.*$x4Damp/($omegaI+$omegaJ)];
set betaK4init [expr $K4initSwitch*2.*$x4Damp/($omegaI+$omegaJ)];
rayleigh $alphaM2 $betaK4curr $betaK4init $betaK4comm;
region 2 -eleRange 4 6 rayleigh $alphaM4 $betaK4curr $betaK4init $betaK4comm;
region 4 -nodeRange 5 6 rayleigh $alphaM4 $betaK4curr $betaK4init $betaK4comm;
} else {
constraints Plain
numberer Plain
system BandGeneral
test NormDispIncr 1.0e-10 10
algorithm Newton
integrator Newmark 0.5 0.25
analysis Transient
analyze 1 0.01
set i [expr $i+1]
}
}

thank you

which region command suitable for my target? (for assign 2 different damping ratio to beams and columns in each story "story 1 and story2". mass assign to each joint and dont assign mass to elements )
region node or region element?

hi friends
(first): I decided to modeling 2D , 2story frame with 2 different damping ratio (each story have a uniqe damping ratio) . i use region command in tcl but unfortunately dont work correct. i guess , i write uncorrect region command.
(second): after modeling 2story frame with 2 different damping ratio i will write command for this model in order that these damping ration can change together in each time. (C1 and C2 can switching together depending with time)
can any body help me?
my tcl will attach




wipe

model BasicBuilder -ndm 2 -ndf 3

file mkdir nonclasic

set L [expr 5]
set H [expr 5]

node 1 0 0
node 2 $L 0
node 3 0 $H
node 4 $L $H
node 5 0 [expr 2*$H]
node 6 $L [expr 2*$H]

fix 1 1 1 1
fix 2 1 1 1

equalDOF 3 4 1
equalDOF 5 6 1

set M1 [expr 1000]
set M2 [expr 1000]

mass 3 $M1 0. 1e-10
mass 4 $M1 0. 1e-10
mass 5 $M2 0. 1e-10
mass 6 $M2 0. 1e-10

geomTransf Linear 1

set E1 [expr 2.e11]
set E2 [expr 2.e11]

set Fy1 [expr 2.e8]
set Fy2 [expr 2.e8]

uniaxialMaterial Steel01 1 $Fy1 $E1 0.02

section fiberSec 1 {
patch quad 1 2 4 -0.1 -0.1 0.1 -0.1 0.1 -0.09 -0.1 -0.09
patch quad 1 2 8 -0.005 -0.09 0.005 -0.09 0.005 0.09 -0.005 0.09
patch quad 1 2 4 -0.1 0.09 0.1 0.09 0.1 0.1 -0.1 0.1
}

uniaxialMaterial Steel01 2 $Fy2 $E2 0.02

section fiberSec 2 {
patch quad 2 2 4 -0.1 -0.1 0.1 -0.1 0.1 -0.09 -0.1 -0.09
patch quad 2 2 8 -0.005 -0.09 0.005 -0.09 0.005 0.09 -0.005 0.09
patch quad 2 2 4 -0.1 0.09 0.1 0.09 0.1 0.1 -0.1 0.1
}

#set I1 [expr 4.1e-5]
#set I2 [expr 4.1e-5]
#set K1 [expr (2*12*$E1*$I1)/pow($H,3)]
#set K2 [expr (2*12*$E2*$I2)/pow($H,3)]
#set X1 [expr 0.15]
#set X2 [expr 0.1]
#set K [expr $K1+$K2]
#set M [expr $M1+$M2]
#set W [expr sqrt($K/$M)]
#set C1 [expr $X1*(2*$M*$W)]
#set C2 [expr $X2*(2*$M*$W)]

element nonlinearBeamColumn 1 1 3 10 1 1
element nonlinearBeamColumn 2 2 4 10 1 1
element nonlinearBeamColumn 3 3 4 10 1 1
element nonlinearBeamColumn 4 3 5 10 2 1
element nonlinearBeamColumn 5 4 6 10 2 1
element nonlinearBeamColumn 6 5 6 10 2 1

set accel "Series -dt 0.01 -filePath ElCentroArray9180(dt=0.01)PGA=0.4g.txt -factor [expr 10*9.81]"

pattern UniformExcitation 3 1 -accel $accel

set x1Damp 0.1;
set M1propSwitch 1.0;
set K1currSwitch 0.0;
set K1commSwitch 1.0;
set K1initSwitch 0.0;
set nEigenI 1;
set nEigenJ 2;
set lambdaN [eigen [expr $nEigenJ]];
set lambdaI [lindex $lambdaN [expr $nEigenI-1]];
set lambdaJ [lindex $lambdaN [expr $nEigenJ-1]];
set omegaI [expr pow($lambdaI,0.5)];
set omegaJ [expr pow($lambdaJ,0.5)];
set alphaM1 [expr $M1propSwitch*$x1Damp*(2*$omegaI*$omegaJ)/($omegaI+$omegaJ)];
set betaK1curr [expr $K1currSwitch*2.*$x1Damp/($omegaI+$omegaJ)];
set betaK1comm [expr $K1commSwitch*2.*$x1Damp/($omegaI+$omegaJ)];
set betaK1init [expr $K1initSwitch*2.*$x1Damp/($omegaI+$omegaJ)];

region 1 -eleRange 1 3 rayleigh $alphaM1 $betaK1curr $betaK1init $betaK1comm;
region 3 -nodeRange 1 4 rayleigh $alphaM1 $betaK1curr $betaK1init $betaK1comm;

rayleigh $alphaM1 $betaK1curr $betaK1init $betaK1comm;

set x2Damp 0.01;
set M2propSwitch 1.0;
set K2currSwitch 0.0;
set K2commSwitch 1.0;
set K2initSwitch 0.0;
set nEigenI 1;
set nEigenJ 2;
set lambdaN [eigen [expr $nEigenJ]];
set lambdaI [lindex $lambdaN [expr $nEigenI-1]];
set lambdaJ [lindex $lambdaN [expr $nEigenJ-1]];
set omegaI [expr pow($lambdaI,0.5)];
set omegaJ [expr pow($lambdaJ,0.5)];
set alphaM2 [expr $M2propSwitch*$x2Damp*(2*$omegaI*$omegaJ)/($omegaI+$omegaJ)];
set betaK2curr [expr $K2currSwitch*2.*$x2Damp/($omegaI+$omegaJ)];
set betaK2comm [expr $K2commSwitch*2.*$x2Damp/($omegaI+$omegaJ)];
set betaK2init [expr $K2initSwitch*2.*$x2Damp/($omegaI+$omegaJ)];

region 2 -eleRange 4 6 rayleigh $alphaM2 $betaK2curr $betaK2init $betaK2comm;
region 4 -nodeRange 5 6 rayleigh $alphaM2 $betaK2curr $betaK2init $betaK2comm;

rayleigh $alphaM2 $betaK2curr $betaK2init $betaK2comm;

recorder Node -file nonclasic/node4disp.out -time -node 4 -dof 1 disp
recorder Node -file nonclasic/node6disp.out -time -node 6 -dof 1 disp
recorder Node -file nonclasic/RBase1.out -time -node 1 -dof 1 2 3 reaction
recorder Node -file nonclasic/RBase2.out -time -node 2 -dof 1 2 3 reaction
recorder Drift -file nonclasic/Drift1.out -time -iNode 2 -jNode 4 -dof 1 -perpDirn 2
recorder Drift -file nonclasic/Drift2.out -time -iNode 4 -jNode 6 -dof 1 -perpDirn 2
recorder Element -file nonclasic/FCol1.out -time -ele 1 globalForce;
recorder Element -file nonclasic/FCol2.out -time -ele 2 globalForce;
recorder Element -file nonclasic/FCol4.out -time -ele 4 globalForce;
recorder Element -file nonclasic/FCol5.out -time -ele 5 globalForce;
recorder Element -file nonclasic/FBeam3.out -time -ele 3 globalForce;
recorder Element -file nonclasic/FBeam6.out -time -ele 6 globalForce;

recorder display "nonclasic" 500 0 500 500 -wipe
prp 0 0 1
vup 0 1 0
vpn 0 0 1
viewWindow -10 10 -10 10
display 1 5 5

constraints Plain
numberer Plain
system BandGeneral
test NormDispIncr 1.0e-10 10
algorithm Newton
integrator Newmark 0.5 0.25
analyze 5000 0.01

#for {set i 1} {$i <= 5000} {incr i 1} {
#analyze 5000 0.01
if {[expr $Dispnode4 1] > 0.01} {
puts "[getTime] [nodeDisp 4 1]"
set x3Damp 1;
set M3propSwitch 1.0;
set K3currSwitch 0.0;
set K3commSwitch 1.0;
set K3initSwitch 0.0;
set nEigenI 1;
set nEigenJ 2;
set lambdaN [eigen [expr $nEigenJ]];
set lambdaI [lindex $lambdaN [expr $nEigenI-1]];
set lambdaJ [lindex $lambdaN [expr $nEigenJ-1]];
set omegaI [expr pow($lambdaI,0.5)];
set omegaJ [expr pow($lambdaJ,0.5)];
set alphaM3 [expr $M3propSwitch*$x3Damp*(2*$omegaI*$omegaJ)/($omegaI+$omegaJ)];
set betaK3curr [expr $K3currSwitch*2.*$x3Damp/($omegaI+$omegaJ)];
set betaK3comm [expr $K3commSwitch*2.*$x3Damp/($omegaI+$omegaJ)];
set betaK3init [expr $K3initSwitch*2.*$x3Damp/($omegaI+$omegaJ)];
rayleigh $alphaM3 $betaK3curr $betaK3init $betaK3comm;
region 1 -eleRange 1 3 rayleigh $alphaM3 $betaK3curr $betaK3init $betaK3comm;
region 3 -nodeRange 1 4 rayleigh $alphaM3 $betaK3curr $betaK3init $betaK3comm;
} if {[expr abs([nodeDisp 6 1])]>0.001} {
puts "[getTime] [nodeDisp 6 1]"
set x4Damp 2;
set M2propSwitch 1.0;
set K4currSwitch 0.0;
set K4commSwitch 1.0;
set K4initSwitch 0.0;
set nEigenI 1;
set nEigenJ 2;
set lambdaN [eigen [expr $nEigenJ]];
set lambdaI [lindex $lambdaN [expr $nEigenI-1]];
set lambdaJ [lindex $lambdaN [expr $nEigenJ-1]];
set omegaI [expr pow($lambdaI,0.5)];
set omegaJ [expr pow($lambdaJ,0.5)];
set alphaM4 [expr $M4propSwitch*$x4Damp*(2*$omegaI*$omegaJ)/($omegaI+$omegaJ)];
set betaK4curr [expr $K4currSwitch*2.*$x4Damp/($omegaI+$omegaJ)];
set betaK4comm [expr $K4commSwitch*2.*$x4Damp/($omegaI+$omegaJ)];
set betaK4init [expr $K4initSwitch*2.*$x4Damp/($omegaI+$omegaJ)];
rayleigh $alphaM2 $betaK4curr $betaK4init $betaK4comm;
region 2 -eleRange 4 6 rayleigh $alphaM4 $betaK4curr $betaK4init $betaK4comm;
region 4 -nodeRange 5 6 rayleigh $alphaM4 $betaK4curr $betaK4init $betaK4comm;
} else {
constraints Plain
numberer Plain
system BandGeneral
test NormDispIncr 1.0e-10 10
algorithm Newton
integrator Newmark 0.5 0.25
analysis Transient
analyze 1 0.01
set i [expr $i+1]
}
}

set a [eigen 2]
set eigen1 [lindex $a 0]
set eigen2 [lindex $a 1]
puts "eigens are $eigen1,$eigen2"
set Period1 [expr 2*3.1416/sqrt($eigen1)]
set Period2 [expr 2*3.1416/sqrt($eigen2)]
puts "Periods are $Period1,$Period2"

set numModes 2
for { set k 1 } { $k <= $numModes } { incr k } {
recorder Node -file nonclasic/mode$k.out -nodeRange 1 6 -dof 1 2 3 "eigen $k"
}

set lambda [eigen $numModes];
set omega {}
set f {}
set T {}
set pi 3.141593

foreach lam $lambda {
lappend omega [expr sqrt($lam)]
lappend f [expr sqrt($lam)/(2*$pi)]
lappend T [expr (2*$pi)/sqrt($lam)]
}

set period "Periods.txt"
set Periods [open $period "w"]
foreach t $T {
puts $Periods " $t"
}
close $Periods



thank you

setare u must use this code


geomTransf Corotational 1 0 0 -1 ; #colmn
geomTransf Corotational 2 0 0 1 ; #in plane beam
geomTransf Corotational 3 1 0 0 ; #out plane beam

be successfull

i thank you for reply.
i'm model friction damper according to what Dr. Imad H. Mualla used his papers. and i want to change this code as from friction damper to rotation friction damper.please help me to write this code.
for represented friction damper i used element zeroLength frome node 14 to 13. "element zeroLength 12 14 13 -mat 1 -dir 6"
best regards

Dear Vesna
please help me to correct this error, i'm so confuse to correct < corotCrdTrans3d: compute elemt length and orien: 0 length
< force beam column 3d: set domain .. error initializing coordinate transformation >

wipe;
model BasicBuilder -ndm 3 -ndf 6;
file mkdir payam;
set dataDir payam;
set H 3.5
set LB 5
set LG 5
node 1 0 0 0
node 2 [expr {1*$LB}] 0 0
node 3 [expr {2*$LB}] 0 0
node 4 [expr {3*$LB}] 0 0
node 5 [expr {4*$LB}] 0 0
node 6 [expr {5*$LB}] 0 0
node 7 [expr {6*$LB}] 0 0
node 8 0 0 [expr {1*$LG}]
node 9 [expr {1*$LB}] 0 [expr {1*$LG}]
node 10 [expr {2*$LB}] 0 [expr {1*$LG}]
node 11 [expr {3*$LB}] 0 [expr {1*$LG}]
node 12 [expr {4*$LB}] 0 [expr {1*$LG}]
node 13 [expr {5*$LB}] 0 [expr {1*$LG}]
node 14 [expr {6*$LB}] 0 [expr {1*$LG}]
node 15 0 0 [expr {2*$LG}]
node 16 [expr {1*$LB}] 0 [expr {2*$LG}]
node 17 [expr {2*$LB}] 0 [expr {2*$LG}]
node 18 [expr {3*$LB}] 0 [expr {2*$LG}]
node 19 [expr {4*$LB}] 0 [expr {2*$LG}]
node 20 [expr {5*$LB}] 0 [expr {2*$LG}]
node 21 [expr {6*$LB}] 0 [expr {2*$LG}]
node 22 0 0 [expr {3*$LG}]
node 23 [expr {1*$LB}] 0 [expr {3*$LG}]
node 24 [expr {2*$LB}] 0 [expr {3*$LG}]
node 25 [expr {3*$LB}] 0 [expr {3*$LG}]
node 26 [expr {4*$LB}] 0 [expr {3*$LG}]
node 27 [expr {5*$LB}] 0 [expr {3*$LG}]
node 28 [expr {6*$LB}] 0 [expr {3*$LG}]
node 29 0 0 [expr {4*$LG}]
node 30 [expr {1*$LB}] 0 [expr {4*$LG}]
node 31 [expr {2*$LB}] 0 [expr {4*$LG}]
node 32 [expr {3*$LB}] 0 [expr {4*$LG}]
node 33 0 0 [expr {5*$LG}]
node 34 [expr {1*$LB}] 0 [expr {5*$LG}]
node 35 [expr {2*$LB}] 0 [expr {5*$LG}]
node 36 [expr {3*$LB}] 0 [expr {5*$LG}]
node 37 0 0 [expr {6*$LG}]
node 38 [expr {1*$LB}] 0 [expr {6*$LG}]
node 39 [expr {2*$LB}] 0 [expr {6*$LG}]
node 40 [expr {3*$LB}] 0 [expr {6*$LG}]
node 41 0 [expr {1*$H}] 0
node 42 [expr {1*$LB}] [expr {1*$H}] 0
node 43 [expr {2*$LB}] [expr {1*$H}] 0
node 44 [expr {3*$LB}] [expr {1*$H}] 0
node 45 [expr {4*$LB}] [expr {1*$H}] 0
node 46 [expr {5*$LB}] [expr {1*$H}] 0
node 47 [expr {6*$LB}] [expr {1*$H}] 0
node 48 0 [expr {1*$H}] [expr {1*$LG}]
node 49 [expr {1*$LB}] [expr {1*$H}] [expr {1*$LG}]
node 50 [expr {2*$LB}] [expr {1*$H}] [expr {1*$LG}]
node 51 [expr {3*$LB}] [expr {1*$H}] [expr {1*$LG}]
node 52 [expr {4*$LB}] [expr {1*$H}] [expr {1*$LG}]
node 53 [expr {5*$LB}] [expr {1*$H}] [expr {1*$LG}]
node 54 [expr {6*$LB}] [expr {1*$H}] [expr {1*$LG}]
node 55 0 [expr {1*$H}] [expr {2*$LG}]
node 56 [expr {1*$LB}] [expr {1*$H}] [expr {2*$LG}]
node 57 [expr {2*$LB}] [expr {1*$H}] [expr {2*$LG}]
node 58 [expr {3*$LB}] [expr {1*$H}] [expr {2*$LG}]
node 59 [expr {4*$LB}] [expr {1*$H}] [expr {2*$LG}]
node 60 [expr {5*$LB}] [expr {1*$H}] [expr {2*$LG}]
node 61 [expr {6*$LB}] [expr {1*$H}] [expr {2*$LG}]
node 62 0 [expr {1*$H}] [expr {3*$LG}]
node 63 [expr {1*$LB}] [expr {1*$H}] [expr {3*$LG}]
node 64 [expr {2*$LB}] [expr {1*$H}] [expr {3*$LG}]
node 65 [expr {3*$LB}] [expr {1*$H}] [expr {3*$LG}]
node 66 [expr {4*$LB}] [expr {1*$H}] [expr {3*$LG}]
node 67 [expr {5*$LB}] [expr {1*$H}] [expr {3*$LG}]
node 68 [expr {6*$LB}] [expr {1*$H}] [expr {3*$LG}]
node 69 0 [expr {1*$H}] [expr {4*$LG}]
node 70 [expr {1*$LB}] [expr {1*$H}] [expr {4*$LG}]
node 71 [expr {2*$LB}] [expr {1*$H}] [expr {4*$LG}]
node 72 [expr {3*$LB}] [expr {1*$H}] [expr {4*$LG}]
node 73 0 [expr {1*$H}] [expr {5*$LG}]
node 74 [expr {1*$LB}] [expr {1*$H}] [expr {5*$LG}]
node 75 [expr {2*$LB}] [expr {1*$H}] [expr {5*$LG}]
node 76 [expr {3*$LB}] [expr {1*$H}] [expr {5*$LG}]
node 77 0 [expr {1*$H}] [expr {6*$LG}]
node 78 [expr {1*$LB}] [expr {1*$H}] [expr {6*$LG}]
node 79 [expr {2*$LB}] [expr {1*$H}] [expr {6*$LG}]
node 80 [expr {3*$LB}] [expr {1*$H}] [expr {6*$LG}]
node 81 0 [expr {2*$H}] 0
node 82 [expr {1*$LB}] [expr {2*$H}] 0
node 83 [expr {2*$LB}] [expr {2*$H}] 0
node 84 [expr {3*$LB}] [expr {2*$H}] 0
node 85 [expr {4*$LB}] [expr {2*$H}] 0
node 86 [expr {5*$LB}] [expr {2*$H}] 0
node 87 [expr {6*$LB}] [expr {2*$H}] 0
node 88 0 [expr {2*$H}] [expr {1*$LG}]
node 89 [expr {1*$LB}] [expr {2*$H}] [expr {1*$LG}]
node 90 [expr {2*$LB}] [expr {2*$H}] [expr {1*$LG}]
node 91 [expr {3*$LB}] [expr {2*$H}] [expr {1*$LG}]
node 92 [expr {4*$LB}] [expr {2*$H}] [expr {1*$LG}]
node 93 [expr {5*$LB}] [expr {2*$H}] [expr {1*$LG}]
node 94 [expr {6*$LB}] [expr {2*$H}] [expr {1*$LG}]
node 95 0 [expr {2*$H}] [expr {2*$LG}]
node 96 [expr {1*$LB}] [expr {2*$H}] [expr {2*$LG}]
node 97 [expr {2*$LB}] [expr {2*$H}] [expr {2*$LG}]
node 98 [expr {3*$LB}] [expr {2*$H}] [expr {2*$LG}]
node 99 [expr {4*$LB}] [expr {2*$H}] [expr {2*$LG}]
node 100 [expr {5*$LB}] [expr {2*$H}] [expr {2*$LG}]
node 101 [expr {6*$LB}] [expr {2*$H}] [expr {2*$LG}]
node 102 0 [expr {2*$H}] [expr {3*$LG}]
node 103 [expr {1*$LB}] [expr {2*$H}] [expr {3*$LG}]
node 104 [expr {2*$LB}] [expr {2*$H}] [expr {3*$LG}]
node 105 [expr {3*$LB}] [expr {2*$H}] [expr {3*$LG}]
node 106 [expr {4*$LB}] [expr {2*$H}] [expr {3*$LG}]
node 107 [expr {5*$LB}] [expr {2*$H}] [expr {3*$LG}]
node 108 [expr {6*$LB}] [expr {2*$H}] [expr {3*$LG}]
node 109 0 [expr {2*$H}] [expr {4*$LG}]
node 110 [expr {1*$LB}] [expr {2*$H}] [expr {4*$LG}]
node 111 [expr {2*$LB}] [expr {2*$H}] [expr {4*$LG}]
node 112 [expr {3*$LB}] [expr {2*$H}] [expr {4*$LG}]
node 113 0 [expr {2*$H}] [expr {5*$LG}]
node 114 [expr {1*$LB}] [expr {2*$H}] [expr {5*$LG}]
node 115 [expr {2*$LB}] [expr {2*$H}] [expr {5*$LG}]
node 116 [expr {3*$LB}] [expr {2*$H}] [expr {5*$LG}]
node 117 0 [expr {2*$H}] [expr {6*$LG}]
node 118 [expr {1*$LB}] [expr {2*$H}] [expr {6*$LG}]
node 119 [expr {2*$LB}] [expr {2*$H}] [expr {6*$LG}]
node 120 [expr {3*$LB}] [expr {2*$H}] [expr {6*$LG}]
node 121 0 [expr {3*$H}] 0
node 122 [expr {1*$LB}] [expr {3*$H}] 0
node 123 [expr {2*$LB}] [expr {3*$H}] 0
node 124 [expr {3*$LB}] [expr {3*$H}] 0
node 125 [expr {4*$LB}] [expr {3*$H}] 0
node 126 [expr {5*$LB}] [expr {3*$H}] 0
node 127 [expr {6*$LB}] [expr {3*$H}] 0
node 128 0 [expr {3*$H}] [expr {1*$LG}]
node 129 [expr {1*$LB}] [expr {3*$H}] [expr {1*$LG}]
node 130 [expr {2*$LB}] [expr {3*$H}] [expr {1*$LG}]
node 131 [expr {3*$LB}] [expr {3*$H}] [expr {1*$LG}]
node 132 [expr {4*$LB}] [expr {3*$H}] [expr {1*$LG}]
node 133 [expr {5*$LB}] [expr {3*$H}] [expr {1*$LG}]
node 134 [expr {6*$LB}] [expr {3*$H}] [expr {1*$LG}]
node 135 0 [expr {3*$H}] [expr {2*$LG}]
node 136 [expr {1*$LB}] [expr {3*$H}] [expr {2*$LG}]
node 137 [expr {2*$LB}] [expr {3*$H}] [expr {2*$LG}]
node 138 [expr {3*$LB}] [expr {3*$H}] [expr {2*$LG}]
node 139 [expr {4*$LB}] [expr {3*$H}] [expr {2*$LG}]
node 140 [expr {5*$LB}] [expr {3*$H}] [expr {2*$LG}]
node 141 [expr {6*$LB}] [expr {3*$H}] [expr {2*$LG}]
node 142 0 [expr {3*$H}] [expr {3*$LG}]
node 143 [expr {1*$LB}] [expr {3*$H}] [expr {3*$LG}]
node 144 [expr {2*$LB}] [expr {3*$H}] [expr {3*$LG}]
node 145 [expr {3*$LB}] [expr {3*$H}] [expr {3*$LG}]
node 146 [expr {4*$LB}] [expr {3*$H}] [expr {3*$LG}]
node 147 [expr {5*$LB}] [expr {3*$H}] [expr {3*$LG}]
node 148 [expr {6*$LB}] [expr {3*$H}] [expr {3*$LG}]
node 149 0 [expr {3*$H}] [expr {4*$LG}]
node 150 [expr {1*$LB}] [expr {3*$H}] [expr {4*$LG}]
node 151 [expr {2*$LB}] [expr {3*$H}] [expr {4*$LG}]
node 152 [expr {3*$LB}] [expr {3*$H}] [expr {4*$LG}]
node 153 0 [expr {3*$H}] [expr {5*$LG}]
node 154 [expr {1*$LB}] [expr {3*$H}] [expr {5*$LG}]
node 155 [expr {2*$LB}] [expr {3*$H}] [expr {5*$LG}]
node 156 [expr {3*$LB}] [expr {3*$H}] [expr {5*$LG}]
node 157 0 [expr {3*$H}] [expr {6*$LG}]
node 158 [expr {1*$LB}] [expr {3*$H}] [expr {6*$LG}]
node 159 [expr {2*$LB}] [expr {3*$H}] [expr {6*$LG}]
node 160 [expr {3*$LB}] [expr {3*$H}] [expr {6*$LG}]
node 161 0 [expr {4*$H}] 0
node 162 [expr {1*$LB}] [expr {4*$H}] 0
node 163 [expr {2*$LB}] [expr {4*$H}] 0
node 164 [expr {3*$LB}] [expr {4*$H}] 0
node 165 [expr {4*$LB}] [expr {4*$H}] 0
node 166 [expr {5*$LB}] [expr {4*$H}] 0
node 167 [expr {6*$LB}] [expr {4*$H}] 0
node 168 0 [expr {4*$H}] [expr {1*$LG}]
node 169 [expr {1*$LB}] [expr {4*$H}] [expr {1*$LG}]
node 170 [expr {2*$LB}] [expr {4*$H}] [expr {1*$LG}]
node 171 [expr {3*$LB}] [expr {4*$H}] [expr {1*$LG}]
node 172 [expr {4*$LB}] [expr {4*$H}] [expr {1*$LG}]
node 173 [expr {5*$LB}] [expr {4*$H}] [expr {1*$LG}]
node 174 [expr {6*$LB}] [expr {4*$H}] [expr {1*$LG}]
node 175 0 [expr {4*$H}] [expr {2*$LG}]
node 176 [expr {1*$LB}] [expr {4*$H}] [expr {2*$LG}]
node 177 [expr {2*$LB}] [expr {4*$H}] [expr {2*$LG}]
node 178 [expr {3*$LB}] [expr {4*$H}] [expr {2*$LG}]
node 179 [expr {4*$LB}] [expr {4*$H}] [expr {2*$LG}]
node 180 [expr {5*$LB}] [expr {4*$H}] [expr {2*$LG}]
node 181 [expr {6*$LB}] [expr {4*$H}] [expr {2*$LG}]
node 182 0 [expr {4*$H}] [expr {3*$LG}]
node 183 [expr {1*$LB}] [expr {4*$H}] [expr {3*$LG}]
node 184 [expr {2*$LB}] [expr {4*$H}] [expr {3*$LG}]
node 185 [expr {3*$LB}] [expr {4*$H}] [expr {3*$LG}]
node 186 [expr {4*$LB}] [expr {4*$H}] [expr {3*$LG}]
node 187 [expr {5*$LB}] [expr {4*$H}] [expr {3*$LG}]
node 188 [expr {6*$LB}] [expr {4*$H}] [expr {3*$LG}]
node 189 0 [expr {4*$H}] [expr {4*$LG}]
node 190 [expr {1*$LB}] [expr {4*$H}] [expr {4*$LG}]
node 191 [expr {2*$LB}] [expr {4*$H}] [expr {4*$LG}]
node 192 [expr {3*$LB}] [expr {4*$H}] [expr {4*$LG}]
node 193 0 [expr {4*$H}] [expr {5*$LG}]
node 194 [expr {1*$LB}] [expr {4*$H}] [expr {5*$LG}]
node 195 [expr {2*$LB}] [expr {4*$H}] [expr {5*$LG}]
node 196 [expr {3*$LB}] [expr {4*$H}] [expr {5*$LG}]
node 197 0 [expr {4*$H}] [expr {6*$LG}]
node 198 [expr {1*$LB}] [expr {4*$H}] [expr {6*$LG}]
node 199 [expr {2*$LB}] [expr {4*$H}] [expr {6*$LG}]
node 200 [expr {3*$LB}] [expr {4*$H}] [expr {6*$LG}]
node 201 [expr {2.5*$LB}] [expr {1*$H}] [expr {2.5*$LG}]
node 202 [expr {2.5*$LB}] [expr {2*$H}] [expr {2.5*$LG}]
node 203 [expr {2.5*$LB}] [expr {3*$H}] [expr {2.5*$LG}]
node 204 [expr {2.5*$LB}] [expr {4*$H}] [expr {2.5*$LG}]
fix 201 0 1 0 1 0 1
fix 202 0 1 0 1 0 1
fix 203 0 1 0 1 0 1
fix 204 0 1 0 1 0 1
fixY 0 1 1 1 1 1 1
rigidDiaphragm 2 201 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80
rigidDiaphragm 2 202 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120
rigidDiaphragm 2 203 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160
rigidDiaphragm 2 204 161 162 163 164 165 166 127 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 158 159 200
mass 41 1000 1e-10 1000 0 0 0
mass 42 1000 1e-10 1000 0 0 0
mass 43 1000 1e-10 1000 0 0 0
mass 44 1000 1e-10 1000 0 0 0
mass 45 1000 1e-10 1000 0 0 0
mass 46 1000 1e-10 1000 0 0 0
mass 47 1000 1e-10 1000 0 0 0
mass 48 1000 1e-10 1000 0 0 0
mass 49 1000 1e-10 1000 0 0 0
mass 50 1000 1e-10 1000 0 0 0
mass 51 1000 1e-10 1000 0 0 0
mass 52 1000 1e-10 1000 0 0 0
mass 53 1000 1e-10 1000 0 0 0
mass 54 1000 1e-10 1000 0 0 0
mass 55 1000 1e-10 1000 0 0 0
mass 56 1000 1e-10 1000 0 0 0
mass 57 1000 1e-10 1000 0 0 0
mass 58 1000 1e-10 1000 0 0 0
mass 59 1000 1e-10 1000 0 0 0
mass 60 1000 1e-10 1000 0 0 0
mass 61 1000 1e-10 1000 0 0 0
mass 62 1000 1e-10 1000 0 0 0
mass 63 1000 1e-10 1000 0 0 0
mass 64 1000 1e-10 1000 0 0 0
mass 65 1000 1e-10 1000 0 0 0
mass 66 1000 1e-10 1000 0 0 0
mass 67 1000 1e-10 1000 0 0 0
mass 68 1000 1e-10 1000 0 0 0
mass 69 1000 1e-10 1000 0 0 0
mass 70 1000 1e-10 1000 0 0 0
mass 71 1000 1e-10 1000 0 0 0
mass 72 1000 1e-10 1000 0 0 0
mass 73 1000 1e-10 1000 0 0 0
mass 74 1000 1e-10 1000 0 0 0
mass 75 1000 1e-10 1000 0 0 0
mass 76 1000 1e-10 1000 0 0 0
mass 77 1000 1e-10 1000 0 0 0
mass 78 1000 1e-10 1000 0 0 0
mass 79 1000 1e-10 1000 0 0 0
mass 80 1000 1e-10 1000 0 0 0
mass 81 1000 1e-10 1000 0 0 0
mass 82 1000 1e-10 1000 0 0 0
mass 83 1000 1e-10 1000 0 0 0
mass 84 1000 1e-10 1000 0 0 0
mass 85 1000 1e-10 1000 0 0 0
mass 86 1000 1e-10 1000 0 0 0
mass 87 1000 1e-10 1000 0 0 0
mass 88 1000 1e-10 1000 0 0 0
mass 89 1000 1e-10 1000 0 0 0
mass 90 1000 1e-10 1000 0 0 0
mass 91 1000 1e-10 1000 0 0 0
mass 92 1000 1e-10 1000 0 0 0
mass 93 1000 1e-10 1000 0 0 0
mass 94 1000 1e-10 1000 0 0 0
mass 95 1000 1e-10 1000 0 0 0
mass 96 1000 1e-10 1000 0 0 0
mass 97 1000 1e-10 1000 0 0 0
mass 98 1000 1e-10 1000 0 0 0
mass 99 1000 1e-10 1000 0 0 0
mass 100 1000 1e-10 1000 0 0 0
mass 101 1000 1e-10 1000 0 0 0
mass 102 1000 1e-10 1000 0 0 0
mass 103 1000 1e-10 1000 0 0 0
mass 104 1000 1e-10 1000 0 0 0
mass 105 1000 1e-10 1000 0 0 0
mass 106 1000 1e-10 1000 0 0 0
mass 107 1000 1e-10 1000 0 0 0
mass 108 1000 1e-10 1000 0 0 0
mass 109 1000 1e-10 1000 0 0 0
mass 110 1000 1e-10 1000 0 0 0
mass 111 1000 1e-10 1000 0 0 0
mass 112 1000 1e-10 1000 0 0 0
mass 113 1000 1e-10 1000 0 0 0
mass 114 1000 1e-10 1000 0 0 0
mass 115 1000 1e-10 1000 0 0 0
mass 116 1000 1e-10 1000 0 0 0
mass 117 1000 1e-10 1000 0 0 0
mass 118 1000 1e-10 1000 0 0 0
mass 119 1000 1e-10 1000 0 0 0
mass 120 1000 1e-10 1000 0 0 0
mass 121 1000 1e-10 1000 0 0 0
mass 122 1000 1e-10 1000 0 0 0
mass 123 1000 1e-10 1000 0 0 0
mass 124 1000 1e-10 1000 0 0 0
mass 125 1000 1e-10 1000 0 0 0
mass 126 1000 1e-10 1000 0 0 0
mass 127 1000 1e-10 1000 0 0 0
mass 128 1000 1e-10 1000 0 0 0
mass 129 1000 1e-10 1000 0 0 0
mass 130 1000 1e-10 1000 0 0 0
mass 131 1000 1e-10 1000 0 0 0
mass 132 1000 1e-10 1000 0 0 0
mass 133 1000 1e-10 1000 0 0 0
mass 134 1000 1e-10 1000 0 0 0
mass 135 1000 1e-10 1000 0 0 0
mass 136 1000 1e-10 1000 0 0 0
mass 137 1000 1e-10 1000 0 0 0
mass 138 1000 1e-10 1000 0 0 0
mass 139 1000 1e-10 1000 0 0 0
mass 140 1000 1e-10 1000 0 0 0
mass 141 1000 1e-10 1000 0 0 0
mass 142 1000 1e-10 1000 0 0 0
mass 143 1000 1e-10 1000 0 0 0
mass 144 1000 1e-10 1000 0 0 0
mass 145 1000 1e-10 1000 0 0 0
mass 146 1000 1e-10 1000 0 0 0
mass 147 1000 1e-10 1000 0 0 0
mass 148 1000 1e-10 1000 0 0 0
mass 149 1000 1e-10 1000 0 0 0
mass 150 1000 1e-10 1000 0 0 0
mass 151 1000 1e-10 1000 0 0 0
mass 152 1000 1e-10 1000 0 0 0
mass 153 1000 1e-10 1000 0 0 0
mass 154 1000 1e-10 1000 0 0 0
mass 155 1000 1e-10 1000 0 0 0
mass 156 1000 1e-10 1000 0 0 0
mass 157 1000 1e-10 1000 0 0 0
mass 158 1000 1e-10 1000 0 0 0
mass 159 1000 1e-10 1000 0 0 0
mass 160 1000 1e-10 1000 0 0 0
mass 161 1000 1e-10 1000 0 0 0
mass 162 1000 1e-10 1000 0 0 0
mass 163 1000 1e-10 1000 0 0 0
mass 164 1000 1e-10 1000 0 0 0
mass 165 1000 1e-10 1000 0 0 0
mass 166 1000 1e-10 1000 0 0 0
mass 167 1000 1e-10 1000 0 0 0
mass 168 1000 1e-10 1000 0 0 0
mass 169 1000 1e-10 1000 0 0 0
mass 170 1000 1e-10 1000 0 0 0
mass 171 1000 1e-10 1000 0 0 0
mass 172 1000 1e-10 1000 0 0 0
mass 173 1000 1e-10 1000 0 0 0
mass 174 1000 1e-10 1000 0 0 0
mass 175 1000 1e-10 1000 0 0 0
mass 176 1000 1e-10 1000 0 0 0
mass 177 1000 1e-10 1000 0 0 0
mass 178 1000 1e-10 1000 0 0 0
mass 179 1000 1e-10 1000 0 0 0
mass 180 1000 1e-10 1000 0 0 0
mass 181 1000 1e-10 1000 0 0 0
mass 182 1000 1e-10 1000 0 0 0
mass 183 1000 1e-10 1000 0 0 0
mass 184 1000 1e-10 1000 0 0 0
mass 185 1000 1e-10 1000 0 0 0
mass 186 1000 1e-10 1000 0 0 0
mass 187 1000 1e-10 1000 0 0 0
mass 188 1000 1e-10 1000 0 0 0
mass 189 1000 1e-10 1000 0 0 0
mass 190 1000 1e-10 1000 0 0 0
mass 191 1000 1e-10 1000 0 0 0
mass 192 1000 1e-10 1000 0 0 0
mass 193 1000 1e-10 1000 0 0 0
mass 194 1000 1e-10 1000 0 0 0
mass 195 1000 1e-10 1000 0 0 0
mass 196 1000 1e-10 1000 0 0 0
mass 197 1000 1e-10 1000 0 0 0
mass 198 1000 1e-10 1000 0 0 0
mass 199 1000 1e-10 1000 0 0 0
mass 200 1000 1e-10 1000 0 0 0
set IDconcCore 11;
set IDconcCover 12;
set IDreinf 13;
set fc [expr -2.1e6];
set Ec [expr 2.15e9];
set Kfc 1;
set fc1C [expr $Kfc*$fc];
set eps1C [expr 2.*$fc1C/$Ec];
set fc2C [expr 0.2*$fc1C];
set eps2C [expr 5*$eps1C];
set fc1U $fc;
set eps1U -0.003;
set fc2U [expr 0.2*$fc1U];
set eps2U -0.01;
set lambda 0.1;
set ftC [expr -0.14*$fc1C];
set ftU [expr -0.14*$fc1U];
set Ets [expr $ftU/0.002];
set Fy [expr 2.4e7];
set Es [expr 2.05e10];
set Bs 0.01;
set R0 18;
set cR1 0.925;
set cR2 0.15;
set HSec [expr 0.3];
set BSec [expr 0.3];
set coverH [expr 0.05];
set coverB [expr 0.05];
set numBarsTop 3;
set numBarsBot 3;
set numBarsIntTot 1;
set barAreaTop [expr 3.141592654];
set barAreaBot [expr 3.141592654];
set barAreaInt [expr 3.141592654];
set SecTag 1;
set coverY [expr $HSec/2.0];
set coverZ [expr $BSec/2.0];
set coreY [expr $coverY-$coverH];
set coreZ [expr $coverZ-$coverB];
set nfY 10;
set nfZ 2;
set numBarsInt [expr $numBarsIntTot/2];
uniaxialMaterial Concrete02 $IDconcCore $fc1C $eps1C $fc2C $eps2C $lambda $ftC $Ets;
uniaxialMaterial Concrete02 $IDconcCover $fc1U $eps1U $fc2U $eps2U $lambda $ftU $Ets;
uniaxialMaterial Steel02 $IDreinf $Fy $Es $Bs $R0 $cR1 $cR2;
geomTransf Corotational 1 0 0 -1
geomTransf Corotational 2 0 0 1
geomTransf Corotational 3 1 0 0
section fiberSec $SecTag {;
patch quadr $IDconcCore $nfZ $nfY -$coreY $coreZ -$coreY -$coreZ $coreY -$coreZ $coreY $coreZ;
patch quadr $IDconcCover 1 $nfY -$coverY $coverZ -$coreY $coreZ $coreY $coreZ $coverY $coverZ;
patch quadr $IDconcCover 1 $nfY -$coreY -$coreZ -$coverY -$coverZ $coverY -$coverZ $coreY -$coreZ
patch quadr $IDconcCover $nfZ 1 -$coverY $coverZ -$coverY -$coverZ -$coreY -$coreZ -$coreY $coreZ
patch quadr $IDconcCover $nfZ 1 $coreY $coreZ $coreY -$coreZ $coverY -$coverZ $coverY $coverZ
layer straight $IDreinf $numBarsInt $barAreaInt -$coreY $coreZ $coreY $coreZ;
layer straight $IDreinf $numBarsInt $barAreaInt -$coreY -$coreZ $coreY -$coreZ;
layer straight $IDreinf $numBarsTop $barAreaTop $coreY $coreZ $coreY -$coreZ;
layer straight $IDreinf $numBarsBot $barAreaBot -$coreY $coreZ -$coreY -$coreZ;
};
set SecTagTorsion 3;
set SecTag3D 2;
uniaxialMaterial Elastic $SecTagTorsion 1e10;
section Aggregator $SecTag3D $SecTagTorsion T -section $SecTag;
element nonlinearBeamColumn 1 41 42 5 2 2
element nonlinearBeamColumn 2 42 43 5 2 2
element nonlinearBeamColumn 3 43 44 5 2 2
element nonlinearBeamColumn 4 44 45 5 2 2
element nonlinearBeamColumn 5 45 46 5 2 2
element nonlinearBeamColumn 6 46 47 5 2 2
element nonlinearBeamColumn 7 48 49 5 2 2
element nonlinearBeamColumn 8 49 50 5 2 2
element nonlinearBeamColumn 9 50 51 5 2 2
element nonlinearBeamColumn 10 51 52 5 2 2
element nonlinearBeamColumn 11 52 53 5 2 2
element nonlinearBeamColumn 12 53 54 5 2 2
element nonlinearBeamColumn 13 55 56 5 2 2
element nonlinearBeamColumn 14 56 57 5 2 2
element nonlinearBeamColumn 15 57 58 5 2 2
element nonlinearBeamColumn 16 58 59 5 2 2
element nonlinearBeamColumn 17 59 60 5 2 2
element nonlinearBeamColumn 18 60 61 5 2 2
element nonlinearBeamColumn 19 62 63 5 2 2
element nonlinearBeamColumn 20 63 64 5 2 2
element nonlinearBeamColumn 21 64 65 5 2 2
element nonlinearBeamColumn 22 65 66 5 2 2
element nonlinearBeamColumn 23 66 67 5 2 2
element nonlinearBeamColumn 24 67 68 5 2 2
element nonlinearBeamColumn 25 69 70 5 2 2
element nonlinearBeamColumn 26 70 71 5 2 2
element nonlinearBeamColumn 27 71 72 5 2 2
element nonlinearBeamColumn 28 73 74 5 2 2
element nonlinearBeamColumn 29 74 75 5 2 2
element nonlinearBeamColumn 30 75 76 5 2 2
element nonlinearBeamColumn 31 77 78 5 2 2
element nonlinearBeamColumn 32 78 79 5 2 2
element nonlinearBeamColumn 33 79 80 5 2 2
element nonlinearBeamColumn 34 81 82 5 2 2
element nonlinearBeamColumn 35 82 83 5 2 2
element nonlinearBeamColumn 36 83 84 5 2 2
element nonlinearBeamColumn 37 84 85 5 2 2
element nonlinearBeamColumn 38 85 86 5 2 2
element nonlinearBeamColumn 39 86 87 5 2 2
element nonlinearBeamColumn 40 88 89 5 2 2
element nonlinearBeamColumn 41 89 90 5 2 2
element nonlinearBeamColumn 42 90 91 5 2 2
element nonlinearBeamColumn 43 91 92 5 2 2
element nonlinearBeamColumn 44 92 93 5 2 2
element nonlinearBeamColumn 45 93 94 5 2 2
element nonlinearBeamColumn 46 95 96 5 2 2
element nonlinearBeamColumn 47 96 97 5 2 2
element nonlinearBeamColumn 48 97 98 5 2 2
element nonlinearBeamColumn 49 98 99 5 2 2
element nonlinearBeamColumn 50 99 100 5 2 2
element nonlinearBeamColumn 51 100 101 5 2 2
element nonlinearBeamColumn 52 102 103 5 2 2
element nonlinearBeamColumn 53 103 104 5 2 2
element nonlinearBeamColumn 54 104 105 5 2 2
element nonlinearBeamColumn 55 105 106 5 2 2
element nonlinearBeamColumn 56 106 107 5 2 2
element nonlinearBeamColumn 57 107 108 5 2 2
element nonlinearBeamColumn 58 109 110 5 2 2
element nonlinearBeamColumn 59 110 111 5 2 2
element nonlinearBeamColumn 60 111 112 5 2 2
element nonlinearBeamColumn 61 113 114 5 2 2
element nonlinearBeamColumn 62 114 115 5 2 2
element nonlinearBeamColumn 63 115 116 5 2 2
element nonlinearBeamColumn 64 117 118 5 2 2
element nonlinearBeamColumn 65 118 119 5 2 2
element nonlinearBeamColumn 66 119 120 5 2 2
element nonlinearBeamColumn 67 121 122 5 2 2
element nonlinearBeamColumn 68 122 123 5 2 2
element nonlinearBeamColumn 69 123 124 5 2 2
element nonlinearBeamColumn 70 124 125 5 2 2
element nonlinearBeamColumn 71 125 126 5 2 2
element nonlinearBeamColumn 72 126 127 5 2 2
element nonlinearBeamColumn 73 128 129 5 2 2
element nonlinearBeamColumn 74 129 130 5 2 2
element nonlinearBeamColumn 75 130 131 5 2 2
element nonlinearBeamColumn 76 131 132 5 2 2
element nonlinearBeamColumn 77 132 133 5 2 2
element nonlinearBeamColumn 78 133 134 5 2 2
element nonlinearBeamColumn 79 135 136 5 2 2
element nonlinearBeamColumn 80 136 137 5 2 2
element nonlinearBeamColumn 81 137 138 5 2 2
element nonlinearBeamColumn 82 138 139 5 2 2
element nonlinearBeamColumn 83 139 140 5 2 2
element nonlinearBeamColumn 84 140 141 5 2 2
element nonlinearBeamColumn 85 142 143 5 2 2
element nonlinearBeamColumn 86 143 144 5 2 2
element nonlinearBeamColumn 87 144 145 5 2 2
element nonlinearBeamColumn 88 145 146 5 2 2
element nonlinearBeamColumn 89 146 147 5 2 2
element nonlinearBeamColumn 90 147 148 5 2 2
element nonlinearBeamColumn 91 149 150 5 2 2
element nonlinearBeamColumn 92 150 151 5 2 2
element nonlinearBeamColumn 93 151 152 5 2 2
element nonlinearBeamColumn 94 153 154 5 2 2
element nonlinearBeamColumn 95 154 155 5 2 2
element nonlinearBeamColumn 96 155 156 5 2 2
element nonlinearBeamColumn 97 157 158 5 2 2
element nonlinearBeamColumn 98 158 159 5 2 2
element nonlinearBeamColumn 99 159 160 5 2 2
element nonlinearBeamColumn 100 161 162 5 2 2
element nonlinearBeamColumn 101 122 163 5 2 2
element nonlinearBeamColumn 102 163 124 5 2 2
element nonlinearBeamColumn 103 164 165 5 2 2
element nonlinearBeamColumn 104 165 166 5 2 2
element nonlinearBeamColumn 105 166 167 5 2 2
element nonlinearBeamColumn 106 168 169 5 2 2
element nonlinearBeamColumn 107 169 170 5 2 2
element nonlinearBeamColumn 108 170 171 5 2 2
element nonlinearBeamColumn 109 171 172 5 2 2
element nonlinearBeamColumn 110 172 173 5 2 2
element nonlinearBeamColumn 111 173 174 5 2 2
element nonlinearBeamColumn 112 175 176 5 2 2
element nonlinearBeamColumn 113 176 137 5 2 2
element nonlinearBeamColumn 114 177 178 5 2 2
element nonlinearBeamColumn 115 178 179 5 2 2
element nonlinearBeamColumn 116 179 180 5 2 2
element nonlinearBeamColumn 117 180 181 5 2 2
element nonlinearBeamColumn 118 182 183 5 2 2
element nonlinearBeamColumn 119 183 184 5 2 2
element nonlinearBeamColumn 120 184 185 5 2 2
element nonlinearBeamColumn 121 185 186 5 2 2
element nonlinearBeamColumn 122 186 187 5 2 2
element nonlinearBeamColumn 123 187 188 5 2 2
element nonlinearBeamColumn 124 189 190 5 2 2
element nonlinearBeamColumn 125 190 191 5 2 2
element nonlinearBeamColumn 126 191 192 5 2 2
element nonlinearBeamColumn 127 193 194 5 2 2
element nonlinearBeamColumn 128 194 195 5 2 2
element nonlinearBeamColumn 129 195 196 5 2 2
element nonlinearBeamColumn 130 197 198 5 2 2
element nonlinearBeamColumn 131 198 199 5 2 2
element nonlinearBeamColumn 132 199 200 5 2 2
element nonlinearBeamColumn 133 48 41 5 2 3
element nonlinearBeamColumn 134 55 48 5 2 3
element nonlinearBeamColumn 135 62 55 5 2 3
element nonlinearBeamColumn 136 69 62 5 2 3
element nonlinearBeamColumn 137 73 69 5 2 3
element nonlinearBeamColumn 138 77 73 5 2 3
element nonlinearBeamColumn 139 49 42 5 2 3
element nonlinearBeamColumn 140 56 49 5 2 3
element nonlinearBeamColumn 141 63 56 5 2 3
element nonlinearBeamColumn 142 70 63 5 2 3
element nonlinearBeamColumn 143 74 70 5 2 3
element nonlinearBeamColumn 144 78 74 5 2 3
element nonlinearBeamColumn 145 50 43 5 2 3
element nonlinearBeamColumn 146 57 50 5 2 3
element nonlinearBeamColumn 147 64 57 5 2 3
element nonlinearBeamColumn 148 71 64 5 2 3
element nonlinearBeamColumn 149 75 71 5 2 3
element nonlinearBeamColumn 150 79 75 5 2 3
element nonlinearBeamColumn 151 51 44 5 2 3
element nonlinearBeamColumn 152 58 51 5 2 3
element nonlinearBeamColumn 153 65 58 5 2 3
element nonlinearBeamColumn 154 72 65 5 2 3
element nonlinearBeamColumn 155 76 72 5 2 3
element nonlinearBeamColumn 156 80 76 5 2 3
element nonlinearBeamColumn 157 52 45 5 2 3
element nonlinearBeamColumn 158 59 52 5 2 3
element nonlinearBeamColumn 159 66 59 5 2 3
element nonlinearBeamColumn 160 53 46 5 2 3
element nonlinearBeamColumn 161 60 53 5 2 3
element nonlinearBeamColumn 162 67 60 5 2 3
element nonlinearBeamColumn 163 54 47 5 2 3
element nonlinearBeamColumn 164 61 54 5 2 3
element nonlinearBeamColumn 165 68 61 5 2 3
element nonlinearBeamColumn 166 88 81 5 2 3
element nonlinearBeamColumn 167 95 88 5 2 3
element nonlinearBeamColumn 168 102 95 5 2 3
element nonlinearBeamColumn 169 109 102 5 2 3
element nonlinearBeamColumn 170 113 109 5 2 3
element nonlinearBeamColumn 171 117 113 5 2 3
element nonlinearBeamColumn 172 89 82 5 2 3
element nonlinearBeamColumn 173 96 89 5 2 3
element nonlinearBeamColumn 174 103 96 5 2 3
element nonlinearBeamColumn 175 110 103 5 2 3
element nonlinearBeamColumn 176 114 110 5 2 3
element nonlinearBeamColumn 177 118 114 5 2 3
element nonlinearBeamColumn 178 90 83 5 2 3
element nonlinearBeamColumn 179 97 90 5 2 3
element nonlinearBeamColumn 180 104 97 5 2 3
element nonlinearBeamColumn 181 111 104 5 2 3
element nonlinearBeamColumn 182 115 111 5 2 3
element nonlinearBeamColumn 183 119 115 5 2 3
element nonlinearBeamColumn 184 91 84 5 2 3
element nonlinearBeamColumn 185 98 91 5 2 3
element nonlinearBeamColumn 186 105 98 5 2 3
element nonlinearBeamColumn 187 112 105 5 2 3
element nonlinearBeamColumn 188 116 112 5 2 3
element nonlinearBeamColumn 189 120 116 5 2 3
element nonlinearBeamColumn 190 92 85 5 2 3
element nonlinearBeamColumn 191 99 92 5 2 3
element nonlinearBeamColumn 192 106 99 5 2 3
element nonlinearBeamColumn 193 93 86 5 2 3
element nonlinearBeamColumn 194 100 93 5 2 3
element nonlinearBeamColumn 195 107 100 5 2 3
element nonlinearBeamColumn 196 94 87 5 2 3
element nonlinearBeamColumn 197 101 94 5 2 3
element nonlinearBeamColumn 198 108 101 5 2 3
element nonlinearBeamColumn 199 128 121 5 2 3
element nonlinearBeamColumn 200 135 128 5 2 3
element nonlinearBeamColumn 201 142 135 5 2 3
element nonlinearBeamColumn 202 149 142 5 2 3
element nonlinearBeamColumn 203 153 149 5 2 3
element nonlinearBeamColumn 204 157 153 5 2 3
element nonlinearBeamColumn 205 129 122 5 2 3
element nonlinearBeamColumn 206 136 129 5 2 3
element nonlinearBeamColumn 207 143 136 5 2 3
element nonlinearBeamColumn 208 150 143 5 2 3
element nonlinearBeamColumn 209 154 150 5 2 3
element nonlinearBeamColumn 210 158 154 5 2 3
element nonlinearBeamColumn 211 130 123 5 2 3
element nonlinearBeamColumn 212 137 130 5 2 3
element nonlinearBeamColumn 213 144 137 5 2 3
element nonlinearBeamColumn 214 151 144 5 2 3
element nonlinearBeamColumn 215 155 151 5 2 3
element nonlinearBeamColumn 216 159 155 5 2 3
element nonlinearBeamColumn 217 131 124 5 2 3
element nonlinearBeamColumn 218 138 131 5 2 3
element nonlinearBeamColumn 219 145 138 5 2 3
element nonlinearBeamColumn 220 152 145 5 2 3
element nonlinearBeamColumn 221 156 152 5 2 3
element nonlinearBeamColumn 222 160 156 5 2 3
element nonlinearBeamColumn 223 132 125 5 2 3
element nonlinearBeamColumn 224 139 132 5 2 3
element nonlinearBeamColumn 225 146 139 5 2 3
element nonlinearBeamColumn 226 133 126 5 2 3
element nonlinearBeamColumn 227 140 133 5 2 3
element nonlinearBeamColumn 228 147 140 5 2 3
element nonlinearBeamColumn 229 134 127 5 2 3
element nonlinearBeamColumn 230 141 134 5 2 3
element nonlinearBeamColumn 231 148 141 5 2 3
element nonlinearBeamColumn 232 168 161 5 2 3
element nonlinearBeamColumn 233 165 168 5 2 3
element nonlinearBeamColumn 234 182 175 5 2 3
element nonlinearBeamColumn 235 189 182 5 2 3
element nonlinearBeamColumn 236 193 189 5 2 3
element nonlinearBeamColumn 237 197 193 5 2 3
element nonlinearBeamColumn 238 169 162 5 2 3
element nonlinearBeamColumn 239 176 129 5 2 3
element nonlinearBeamColumn 240 183 176 5 2 3
element nonlinearBeamColumn 241 190 183 5 2 3
element nonlinearBeamColumn 242 194 190 5 2 3
element nonlinearBeamColumn 243 198 194 5 2 3
element nonlinearBeamColumn 244 170 163 5 2 3
element nonlinearBeamColumn 245 177 170 5 2 3
element nonlinearBeamColumn 246 184 177 5 2 3
element nonlinearBeamColumn 247 191 184 5 2 3
element nonlinearBeamColumn 248 195 191 5 2 3
element nonlinearBeamColumn 249 199 195 5 2 3
element nonlinearBeamColumn 250 171 164 5 2 3
element nonlinearBeamColumn 251 178 171 5 2 3
element nonlinearBeamColumn 252 185 178 5 2 3
element nonlinearBeamColumn 253 192 185 5 2 3
element nonlinearBeamColumn 254 196 192 5 2 3
element nonlinearBeamColumn 255 200 196 5 2 3
element nonlinearBeamColumn 256 172 165 5 2 3
element nonlinearBeamColumn 257 179 172 5 2 3
element nonlinearBeamColumn 258 186 179 5 2 3
element nonlinearBeamColumn 259 173 166 5 2 3
element nonlinearBeamColumn 260 180 173 5 2 3
element nonlinearBeamColumn 261 187 180 5 2 3
element nonlinearBeamColumn 262 174 167 5 2 3
element nonlinearBeamColumn 263 181 174 5 2 3
element nonlinearBeamColumn 264 188 181 5 2 3
element nonlinearBeamColumn 265 1 41 5 2 1
element nonlinearBeamColumn 266 2 42 5 2 1
element nonlinearBeamColumn 267 3 43 5 2 1
element nonlinearBeamColumn 268 4 44 5 2 1
element nonlinearBeamColumn 269 5 45 5 2 1
element nonlinearBeamColumn 270 6 46 5 2 1
element nonlinearBeamColumn 271 7 47 5 2 1
element nonlinearBeamColumn 272 8 48 5 2 1
element nonlinearBeamColumn 273 9 49 5 2 1
element nonlinearBeamColumn 274 10 50 5 2 1
element nonlinearBeamColumn 275 11 51 5 2 1
element nonlinearBeamColumn 276 12 52 5 2 1
element nonlinearBeamColumn 277 13 53 5 2 1
element nonlinearBeamColumn 278 14 54 5 2 1
element nonlinearBeamColumn 279 15 55 5 2 1
element nonlinearBeamColumn 280 16 56 5 2 1
element nonlinearBeamColumn 281 17 57 5 2 1
element nonlinearBeamColumn 282 18 58 5 2 1
element nonlinearBeamColumn 283 19 59 5 2 1
element nonlinearBeamColumn 284 20 60 5 2 1
element nonlinearBeamColumn 285 21 61 5 2 1
element nonlinearBeamColumn 286 22 62 5 2 1
element nonlinearBeamColumn 287 23 63 5 2 1
element nonlinearBeamColumn 288 24 64 5 2 1
element nonlinearBeamColumn 289 25 65 5 2 1
element nonlinearBeamColumn 290 26 66 5 2 1
element nonlinearBeamColumn 291 27 67 5 2 1
element nonlinearBeamColumn 292 28 68 5 2 1
element nonlinearBeamColumn 293 29 69 5 2 1
element nonlinearBeamColumn 294 30 70 5 2 1
element nonlinearBeamColumn 295 31 71 5 2 1
element nonlinearBeamColumn 296 32 72 5 2 1
element nonlinearBeamColumn 297 33 73 5 2 1
element nonlinearBeamColumn 298 34 74 5 2 1
element nonlinearBeamColumn 299 35 75 5 2 1
element nonlinearBeamColumn 300 36 76 5 2 1
element nonlinearBeamColumn 301 37 77 5 2 1
element nonlinearBeamColumn 302 38 78 5 2 1
element nonlinearBeamColumn 303 39 79 5 2 1
element nonlinearBeamColumn 304 40 80 5 2 1
element nonlinearBeamColumn 305 41 81 5 2 1
element nonlinearBeamColumn 306 42 82 5 2 1
element nonlinearBeamColumn 307 43 83 5 2 1
element nonlinearBeamColumn 308 44 84 5 2 1
element nonlinearBeamColumn 309 45 85 5 2 1
element nonlinearBeamColumn 310 46 86 5 2 1
element nonlinearBeamColumn 311 47 87 5 2 1
element nonlinearBeamColumn 312 48 88 5 2 1
element nonlinearBeamColumn 313 49 89 5 2 1
element nonlinearBeamColumn 314 50 90 5 2 1
element nonlinearBeamColumn 315 51 91 5 2 1
element nonlinearBeamColumn 316 52 92 5 2 1
element nonlinearBeamColumn 317 53 93 5 2 1
element nonlinearBeamColumn 318 54 94 5 2 1
element nonlinearBeamColumn 319 55 95 5 2 1
element nonlinearBeamColumn 320 56 96 5 2 1
element nonlinearBeamColumn 321 57 97 5 2 1
element nonlinearBeamColumn 322 58 98 5 2 1
element nonlinearBeamColumn 323 59 99 5 2 1
element nonlinearBeamColumn 324 60 100 5 2 1
element nonlinearBeamColumn 325 61 101 5 2 1
element nonlinearBeamColumn 326 62 102 5 2 1
element nonlinearBeamColumn 327 63 103 5 2 1
element nonlinearBeamColumn 328 64 104 5 2 1
element nonlinearBeamColumn 329 65 105 5 2 1
element nonlinearBeamColumn 330 66 106 5 2 1
element nonlinearBeamColumn 331 67 107 5 2 1
element nonlinearBeamColumn 332 68 108 5 2 1
element nonlinearBeamColumn 333 69 109 5 2 1
element nonlinearBeamColumn 334 70 110 5 2 1
element nonlinearBeamColumn 335 71 111 5 2 1
element nonlinearBeamColumn 336 72 112 5 2 1
element nonlinearBeamColumn 337 73 113 5 2 1
element nonlinearBeamColumn 338 74 114 5 2 1
element nonlinearBeamColumn 339 75 115 5 2 1
element nonlinearBeamColumn 340 76 116 5 2 1
element nonlinearBeamColumn 341 77 117 5 2 1
element nonlinearBeamColumn 342 78 118 5 2 1
element nonlinearBeamColumn 343 79 119 5 2 1
element nonlinearBeamColumn 344 80 120 5 2 1
element nonlinearBeamColumn 345 81 121 5 2 1
element nonlinearBeamColumn 346 82 122 5 2 1
element nonlinearBeamColumn 347 83 123 5 2 1
element nonlinearBeamColumn 348 84 124 5 2 1
element nonlinearBeamColumn 349 85 125 5 2 1
element nonlinearBeamColumn 350 86 126 5 2 1
element nonlinearBeamColumn 351 87 127 5 2 1
element nonlinearBeamColumn 352 88 128 5 2 1
element nonlinearBeamColumn 353 89 129 5 2 1
element nonlinearBeamColumn 354 90 130 5 2 1
element nonlinearBeamColumn 355 91 131 5 2 1
element nonlinearBeamColumn 356 92 132 5 2 1
element nonlinearBeamColumn 357 93 133 5 2 1
element nonlinearBeamColumn 358 94 134 5 2 1
element nonlinearBeamColumn 359 95 135 5 2 1
element nonlinearBeamColumn 360 96 136 5 2 1
element nonlinearBeamColumn 361 97 137 5 2 1
element nonlinearBeamColumn 362 98 138 5 2 1
element nonlinearBeamColumn 363 99 139 5 2 1
element nonlinearBeamColumn 364 100 140 5 2 1
element nonlinearBeamColumn 365 101 141 5 2 1
element nonlinearBeamColumn 366 102 142 5 2 1
element nonlinearBeamColumn 367 103 143 5 2 1
element nonlinearBeamColumn 368 104 144 5 2 1
element nonlinearBeamColumn 369 105 145 5 2 1
element nonlinearBeamColumn 370 106 146 5 2 1
element nonlinearBeamColumn 371 107 147 5 2 1
element nonlinearBeamColumn 372 108 148 5 2 1
element nonlinearBeamColumn 373 109 139 5 2 1
element nonlinearBeamColumn 374 110 150 5 2 1
element nonlinearBeamColumn 375 111 151 5 2 1
element nonlinearBeamColumn 376 112 152 5 2 1
element nonlinearBeamColumn 377 113 153 5 2 1
element nonlinearBeamColumn 378 114 154 5 2 1
element nonlinearBeamColumn 379 115 155 5 2 1
element nonlinearBeamColumn 380 116 156 5 2 1
element nonlinearBeamColumn 381 117 157 5 2 1
element nonlinearBeamColumn 382 118 158 5 2 1
element nonlinearBeamColumn 383 119 159 5 2 1
element nonlinearBeamColumn 384 120 160 5 2 1
element nonlinearBeamColumn 385 121 161 5 2 1
element nonlinearBeamColumn 386 122 162 5 2 1
element nonlinearBeamColumn 387 123 163 5 2 1
element nonlinearBeamColumn 388 124 164 5 2 1
element nonlinearBeamColumn 389 125 165 5 2 1
element nonlinearBeamColumn 390 126 166 5 2 1
element nonlinearBeamColumn 391 127 167 5 2 1
element nonlinearBeamColumn 392 128 168 5 2 1
element nonlinearBeamColumn 393 129 169 5 2 1
element nonlinearBeamColumn 394 130 170 5 2 1
element nonlinearBeamColumn 395 131 171 5 2 1
element nonlinearBeamColumn 396 132 132 5 2 1
element nonlinearBeamColumn 397 133 173 5 2 1
element nonlinearBeamColumn 398 134 174 5 2 1
element nonlinearBeamColumn 399 135 175 5 2 1
element nonlinearBeamColumn 400 136 176 5 2 1
element nonlinearBeamColumn 401 137 177 5 2 1
element nonlinearBeamColumn 402 138 178 5 2 1
element nonlinearBeamColumn 403 139 179 5 2 1
element nonlinearBeamColumn 404 140 180 5 2 1
element nonlinearBeamColumn 405 141 181 5 2 1
element nonlinearBeamColumn 406 142 182 5 2 1
element nonlinearBeamColumn 407 143 183 5 2 1
element nonlinearBeamColumn 408 144 184 5 2 1
element nonlinearBeamColumn 309 145 185 5 2 1
element nonlinearBeamColumn 410 146 186 5 2 1
element nonlinearBeamColumn 411 147 187 5 2 1
element nonlinearBeamColumn 412 148 188 5 2 1
element nonlinearBeamColumn 413 149 189 5 2 1
element nonlinearBeamColumn 414 150 190 5 2 1
element nonlinearBeamColumn 415 151 191 5 2 1
element nonlinearBeamColumn 416 152 192 5 2 1
element nonlinearBeamColumn 417 153 193 5 2 1
element nonlinearBeamColumn 418 154 194 5 2 1
element nonlinearBeamColumn 419 155 195 5 2 1
element nonlinearBeamColumn 420 156 196 5 2 1
element nonlinearBeamColumn 421 157 197 5 2 1
element nonlinearBeamColumn 422 158 198 5 2 1
element nonlinearBeamColumn 423 159 199 5 2 1
element nonlinearBeamColumn 324 160 200 5 2 1
element pattern Plain 1 Linear {
load 41 0 -1000 0 0 0 0
load 42 0 -1000 0 0 0 0
load 43 0 -1000 0 0 0 0
load 44 0 -1000 0 0 0 0
load 45 0 -1000 0 0 0 0
load 46 0 -1000 0 0 0 0
load 47 0 -1000 0 0 0 0
load 48 0 -1000 0 0 0 0
load 49 0 -1000 0 0 0 0
load 50 0 -1000 0 0 0 0
load 51 0 -1000 0 0 0 0
load 52 0 -1000 0 0 0 0
load 53 0 -1000 0 0 0 0
load 54 0 -1000 0 0 0 0
load 55 0 -1000 0 0 0 0
load 56 0 -1000 0 0 0 0
load 57 0 -1000 0 0 0 0
load 58 0 -1000 0 0 0 0
load 59 0 -1000 0 0 0 0
load 60 0 -1000 0 0 0 0
load 61 0 -1000 0 0 0 0
load 62 0 -1000 0 0 0 0
load 63 0 -1000 0 0 0 0
load 64 0 -1000 0 0 0 0
load 65 0 -1000 0 0 0 0
load 66 0 -1000 0 0 0 0
load 67 0 -1000 0 0 0 0
load 68 0 -1000 0 0 0 0
load 69 0 -1000 0 0 0 0
load 70 0 -1000 0 0 0 0
load 71 0 -1000 0 0 0 0
load 72 0 -1000 0 0 0 0
load 73 0 -1000 0 0 0 0
load 74 0 -1000 0 0 0 0
load 75 0 -1000 0 0 0 0
load 76 0 -1000 0 0 0 0
load 77 0 -1000 0 0 0 0
load 78 0 -1000 0 0 0 0
load 79 0 -1000 0 0 0 0
load 80 0 -1000 0 0 0 0
load 81 0 -1000 0 0 0 0
load 82 0 -1000 0 0 0 0
load 83 0 -1000 0 0 0 0
load 84 0 -1000 0 0 0 0
load 85 0 -1000 0 0 0 0
load 86 0 -1000 0 0 0 0
load 87 0 -1000 0 0 0 0
load 88 0 -1000 0 0 0 0
load 89 0 -1000 0 0 0 0
load 90 0 -1000 0 0 0 0
load 91 0 -1000 0 0 0 0
load 92 0 -1000 0 0 0 0
load 93 0 -1000 0 0 0 0
load 94 0 -1000 0 0 0 0
load 95 0 -1000 0 0 0 0
load 96 0 -1000 0 0 0 0
load 97 0 -1000 0 0 0 0
load 98 0 -1000 0 0 0 0
load 99 0 -1000 0 0 0 0
load 100 0 -1000 0 0 0 0
load 101 0 -1000 0 0 0 0
load 102 0 -1000 0 0 0 0
load 103 0 -1000 0 0 0 0
load 104 0 -1000 0 0 0 0
load 105 0 -1000 0 0 0 0
load 106 0 -1000 0 0 0 0
load 107 0 -1000 0 0 0 0
load 108 0 -1000 0 0 0 0
load 109 0 -1000 0 0 0 0
load 110 0 -1000 0 0 0 0
load 111 0 -1000 0 0 0 0
load 112 0 -1000 0 0 0 0
load 113 0 -1000 0 0 0 0
load 114 0 -1000 0 0 0 0
load 115 0 -1000 0 0 0 0
load 116 0 -1000 0 0 0 0
load 117 0 -1000 0 0 0 0
load 118 0 -1000 0 0 0 0
load 119 0 -1000 0 0 0 0
load 120 0 -1000 0 0 0 0
load 121 0 -1000 0 0 0 0
load 122 0 -1000 0 0 0 0
load 123 0 -1000 0 0 0 0
load 124 0 -1000 0 0 0 0
load 125 0 -1000 0 0 0 0
load 126 0 -1000 0 0 0 0
load 127 0 -1000 0 0 0 0
load 128 0 -1000 0 0 0 0
load 129 0 -1000 0 0 0 0
load 130 0 -1000 0 0 0 0
load 131 0 -1000 0 0 0 0
load 132 0 -1000 0 0 0 0
load 133 0 -1000 0 0 0 0
load 134 0 -1000 0 0 0 0
load 135 0 -1000 0 0 0 0
load 136 0 -1000 0 0 0 0
load 137 0 -1000 0 0 0 0
load 138 0 -1000 0 0 0 0
load 139 0 -1000 0 0 0 0
load 140 0 -1000 0 0 0 0
load 141 0 -1000 0 0 0 0
load 142 0 -1000 0 0 0 0
load 143 0 -1000 0 0 0 0
load 144 0 -1000 0 0 0 0
load 145 0 -1000 0 0 0 0
load 146 0 -1000 0 0 0 0
load 147 0 -1000 0 0 0 0
load 148 0 -1000 0 0 0 0
load 149 0 -1000 0 0 0 0
load 150 0 -1000 0 0 0 0
load 151 0 -1000 0 0 0 0
load 152 0 -1000 0 0 0 0
load 153 0 -1000 0 0 0 0
load 154 0 -1000 0 0 0 0
load 155 0 -1000 0 0 0 0
load 156 0 -1000 0 0 0 0
load 157 0 -1000 0 0 0 0
load 158 0 -1000 0 0 0 0
load 159 0 -1000 0 0 0 0
load 160 0 -1000 0 0 0 0
load 161 0 -1000 0 0 0 0
load 162 0 -1000 0 0 0 0
load 163 0 -1000 0 0 0 0
load 164 0 -1000 0 0 0 0
load 165 0 -1000 0 0 0 0
load 166 0 -1000 0 0 0 0
load 167 0 -1000 0 0 0 0
load 168 0 -1000 0 0 0 0
load 169 0 -1000 0 0 0 0
load 170 0 -1000 0 0 0 0
load 171 0 -1000 0 0 0 0
load 172 0 -1000 0 0 0 0
load 173 0 -1000 0 0 0 0
load 174 0 -1000 0 0 0 0
load 175 0 -1000 0 0 0 0
load 176 0 -1000 0 0 0 0
load 177 0 -1000 0 0 0 0
load 178 0 -1000 0 0 0 0
load 179 0 -1000 0 0 0 0
load 180 0 -1000 0 0 0 0
load 181 0 -1000 0 0 0 0
load 182 0 -1000 0 0 0 0
load 183 0 -1000 0 0 0 0
load 184 0 -1000 0 0 0 0
load 185 0 -1000 0 0 0 0
load 186 0 -1000 0 0 0 0
load 187 0 -1000 0 0 0 0
load 188 0 -1000 0 0 0 0
load 189 0 -1000 0 0 0 0
load 190 0 -1000 0 0 0 0
load 191 0 -1000 0 0 0 0
load 192 0 -1000 0 0 0 0
load 193 0 -1000 0 0 0 0
load 194 0 -1000 0 0 0 0
load 195 0 -1000 0 0 0 0
load 196 0 -1000 0 0 0 0
load 197 0 -1000 0 0 0 0
load 198 0 -1000 0 0 0 0
load 199 0 -1000 0 0 0 0
load 200 0 -1000 0 0 0 0
}
recorder Node -file payam/disp1x.out -time -node 1 -dof 1 disp
recorder Node -file payam/disp1z.out -time -node 1 -dof 3 disp
recorder Node -file payam/node1.out -time -node 1 2 -dof 1 2 3 reaction
recorder Element -file payam/ele11.out -ele 1 section 1 fiber 0 0 stressStrain
recorder Element -file payam/ele12.out -ele 1 section 1 force
recorder display 1 5 5 500 500 -wipe
prp 10 10 10
vup 0 1 0
vpn 0 0 1
display 1 10 10
recorder plot payam/disp1x.out 2 600 200 500 500 -columns 2 1
system BandSPD
numberer RCM
constraints Plain
integrator LoadControl 0.1
algorithm Linear
analysis Static
analyze 10
loadConst -time 0
set a [eigen 3]
set om1 [lindex $a 0]
set om2 [lindex $a 1]
set om3 [lindex $a 2]
set om4 [lindex $a 3]
set Period1 [expr 2*3.14/sqrt($om1)]
set Period2 [expr 2*3.14/sqrt($om2)]
set Period3 [expr 2*3.14/sqrt($om3)]
set Period4 [expr 2*3.14/sqrt($om4)]
puts "$Period1,$Period2,$Period3,$Period4"
#set numEigen 3
#for {set i 1} {$i <= $numEigen} {incr i 1} {
#set eig "eigen $i"
#set recorder "recorder Node -file payam/Eigenvector$i.out -node 3 5 7 -dof 1 {eigen $i}"
#eval $recorder
#}
#analyze 1
set dt 0.01
set AccelSeries "Series -dt $dt -filePath e1.txt -factor 1"
pattern UniformExcitation 100 1 -accel $AccelSeries
recorder Node -file payam/e1 1.out -time -node 1 2 -dof 123 reaction
system SparseGeneral -piv
numberer Plain
constraints Transformation
set Tol 1.e-8 ;
set maxNumIter 10 ;
set printFlag 0 ;
set TestType EnergyIncr ;
test $TestType $Tol $maxNumIter $printFlag ;
set algorithmType Newton ;
algorithm $algorithmType;
set NewmarkGamma 0.5;
set NewmarkBeta 0.25;
integrator Newmark $NewmarkGamma $NewmarkBeta
analysis Transient
set xDamp 0.05;
set lambda [eigen 1];
set omega [expr pow($lambda,0.5)];
set alphaM 0.;
set betaKcurr 0.;
set betaKcomm [expr 2.*$xDamp/($omega)];
set betaKinit 0.;
rayleigh $alphaM $betaKcurr $betaKinit $betaKcomm;
set TmaxAnalysis 25
set DtAnalysis 0.005
set DtMin 1e-6
set DtMax 1e-4
set Nsteps [expr int($TmaxAnalysis/$DtAnalysis)];
set ok [analyze $Nsteps $DtAnalysis];
if {$ok != 0} { ;
set ok 0;
set controlTime [getTime];
while {$controlTime < $TmaxAnalysis && $ok == 0} {
set ok [analyze 1 $DtAnalysis]
set controlTime [getTime]
set ok [analyze 1 $DtAnalysis]
if {$ok != 0} {
puts "Trying Newton with Initial Tangent .."
test NormDispIncr $Tol 1000 0
algorithm Newton -initial
set ok [analyze 1 $DtAnalysis]
test $TestType $Tol $maxNumIter 0
algorithm $algorithmType
}
if {$ok != 0} {
puts "Trying Broyden .."
algorithm Broyden 8
set ok [analyze 1 $DtAnalysis]
algorithm $algorithmType
}
if {$ok != 0} {
puts "Trying NewtonWithLineSearch .."
algorithm NewtonLineSearch .8
set ok [analyze 1 $DtAnalysis]
algorithm $algorithmType
}
}
};
puts "Done.End Time:[getTime]"

BESTS REGARDS

dear vesna ; HELP ME PLEASE

dear vesna ; HELP ME PLEASE