why in this script the period of structure is very large number? T=8.72Sec in 4 story RC frame
I guess that mass value is large but it is correct!!!!!!
before thanks
#Define units--------------------------------------------
set cm 1;
set kg 1;
set cm2 [expr $cm*$cm];
set ubig 1e10;
set usmall [expr 1/$ubig];
#LibMaterials: define a library of materials---------------------------------
#concrete material
set fpc [expr -250*$kg/$cm2];
set epsco -0.003;
set fpcu [expr -84*$kg/$cm2];
set epsu -0.012;
set concrete 1;
uniaxialMaterial Concrete01 $concrete $fpc $epsco $fpcu $epsu;
#Steel material
set Fy [expr 4000*$kg/$cm2];
set E0 [expr 2.06e6*$kg/$cm2];
set b 0.05;
set steel 2;
uniaxialMaterial Steel01 $steel $Fy $E0 $b;
Define sections---------------------------------------------------------------
proc buildrcrectsection { ID H B cover concrete steel numbartop1 barareatop1 numbartop2 barareatop2 numbartopmid barareatopmid numbarbot1 barareabot1 numbarbot2 barareabot2 numbarbotmid barareabotmid numbarint barareaint nfY nfZ } {
set coverY [expr $H/2];
set coverZ [expr $B/2];
set coreY [expr $coverY-$cover];
set coreZ [expr $coverZ-$cover];
set locaddedbar [expr $coreZ/2];
set D1 [expr 1.128*sqrt($barareatop1)]
set D2 [expr 1.128*sqrt($barareatopmid)]
set h1 [expr $H-($D1/2)-($D2/2)]
set D3 [expr 1.128*sqrt($barareabot1)]
set D4 [expr 1.128*sqrt($barareabotmid)]
set h2 [expr $H-($D3/2)-($D4/2)]
set nfY 20;
set nfZ 20:
section Fiber $ID {
patch quad $concrete 20 20 -$coreY $coreZ -$coreY -$coreZ $coreY -$coreZ $coreY $coreZ
patch quad $concrete 2 20 -$coverY $coverZ -$coreY $coreZ $coreY $coreZ $coverY $coverZ
patch quad $concrete 2 20 -$coreY -$coreZ -$coverY -$coverZ $coverY -$coverZ $coreY $coreZ
patch quad $concrete 20 2 -$coverY $coverZ -$coverY -$coverZ -$coreY -$coreZ -$coreY $coreZ
patch quad $concrete 20 2 $coreY $coreZ $coreY -$coreZ $coverY -$coverZ $coverY $coverZ;
layer straight $steel $numbartop1 $barareatop1 $coreY $coreZ $coreY -$coreZ
layer straight $steel $numbartop2 $barareatop2 $coreY $locaddedbar $coreY -$locaddedbar
layer straight $steel $numbartopmid $barareatopmid $h1 0 $h1 0
layer straight $steel $numbarbot1 $barareabot1 -$coreY $coreZ -$coreY -$coreZ
layer straight $steel $numbarbot2 $barareabot2 -$coreY $locaddedbar -$coreY -$locaddedbar
layer straight $steel $numbarbotmid $barareabotmid -$h2 0 -$h2 0
layer straight $steel $numbarint $barareaint -$coreY $coreZ $coreY $coreZ
layer straight $steel $numbarint $barareaint -$coreY -$coreZ $coreY -$coreZ
};
};
#Build model---------------------------------------------------------
model BasicBuilder -ndm 2 -ndf 3;
source libunits.tcl;
#Define Geometry;
set Lcol [expr 300*$cm];
set Lbeam [expr 500*$cm];
set Nstory 4;
set Nbay 3;
#Define nodal coordinates
for {set level 1} {$level <= [expr $Nstory+1]} {incr level 1} {
set y [expr ($level-1)*$Lcol];
for {set pier 1} {$pier <= [expr $Nbay+1]} {incr pier 1} {
set x [expr ($pier-1)*$Lbeam];
set nodeID [expr $level*10+$pier];
node $nodeID $x $y;
}
};
mass 21 [expr 500*$kg] 0 0
mass 22 [expr 4800*$kg] 0 0
mass 23 [expr 4800*$kg] 0 0
mass 24 [expr 500*$kg] 0 0
mass 31 [expr 500*$kg] 0 0
mass 32 [expr 4800*$kg] 0 0
mass 33 [expr 4800*$kg] 0 0
mass 34 [expr 500*$kg] 0 0
mass 41 [expr 500*$kg] 0 0
mass 42 [expr 4800*$kg] 0 0
mass 43 [expr 4800*$kg] 0 0
mass 44 [expr 500*$kg] 0 0
mass 51 [expr 500*$kg] 0 0
mass 52 [expr 4800*$kg] 0 0
mass 53 [expr 4800*$kg] 0 0
mass 54 [expr 500*$kg] 0 0
equalDOF 21 22 1
equalDOF 21 23 1
equalDOF 21 24 1
equalDOF 31 32 1
equalDOF 31 33 1
equalDOF 31 34 1
equalDOF 41 42 1
equalDOF 41 43 1
equalDOF 41 44 1
equalDOF 51 52 1
equalDOF 51 53 1
equalDOF 51 54 1
fix 11 1 1 1;
fix 12 1 1 1;
fix 13 1 1 1;
fix 14 1 1 1;
#Define sections
source libmaterials.tcl;
source buildrcrectsection.tcl;
buildrcrectsection 44 [expr 35*$cm] [expr 30*$cm] [expr 5.4*$cm] $concrete $steel 3 [expr 2.54*$cm2] 0 0 1 [expr 2.54*$cm2] 3 [expr 2.54*$cm2] 0 0 0 0 2 [expr 2.54*$cm2] nfY nfZ
buildrcrectsection 33 [expr 35*$cm] [expr 30*$cm] [expr 5.6*$cm] $concrete $steel 3 [expr 3.8*$cm2] 2 [expr 2.54*$cm2] 0 0 3 [expr 3.8*$cm2] 0 0 0 0 2 [expr 3.8*$cm2] nfY nfZ
buildrcrectsection 22 [expr 40*$cm] [expr 30*$cm] [expr 5.6*$cm] $concrete $steel 3 [expr 3.8*$cm2] 2 [expr 2.54*$cm2] 0 0 3 [expr 3.8*$cm2] 0 0 0 0 2 [expr 3.8*$cm2] nfY nfZ
buildrcrectsection 11 [expr 40*$cm] [expr 30*$cm] [expr 5.6*$cm] $concrete $steel 3 [expr 3.8*$cm2] 0 0 1 [expr 2*$cm2] 3 [expr 3.8*$cm2] 0 0 0 0 2 [expr 3.8*$cm2] nfY nfZ
buildrcrectsection 12 [expr 55*$cm] [expr 55*$cm] [expr 5.6*$cm] $concrete $steel 5 [expr 3.8*$cm2] 0 0 0 0 5 [expr 3.8*$cm2] 0 0 0 0 5 [expr 3.8*$cm2] nfY nfZ
buildrcrectsection 23 [expr 50*$cm] [expr 50*$cm] [expr 5.6*$cm] $concrete $steel 4 [expr 3.8*$cm2] 0 0 0 0 4 [expr 3.8*$cm2] 0 0 0 0 4 [expr 3.8*$cm2] nfY nfZ
buildrcrectsection 34 [expr 40*$cm] [expr 40*$cm] [expr 5.6*$cm] $concrete $steel 3 [expr 3.8*$cm2] 0 0 0 0 3 [expr 3.8*$cm2] 0 0 0 0 3 [expr 3.8*$cm2] nfY nfZ
#Define transformation
geomTransf Linear 1;
geomTransf PDelta 2;
#Define Beam-Column Elements
#beams
element nonlinearBeamColumn 2122 21 22 4 11 1;
element nonlinearBeamColumn 2223 22 23 4 11 1;
element nonlinearBeamColumn 2324 23 24 4 11 1;
element nonlinearBeamColumn 3132 31 32 4 22 1;
element nonlinearBeamColumn 3233 32 33 4 22 1;
element nonlinearBeamColumn 3334 33 34 4 22 1;
element nonlinearBeamColumn 4142 41 42 4 33 1;
element nonlinearBeamColumn 4243 42 43 4 33 1;
element nonlinearBeamColumn 4344 43 44 4 33 1;
element nonlinearBeamColumn 5152 51 52 4 44 1;
element nonlinearBeamColumn 5253 52 53 4 44 1;
element nonlinearBeamColumn 5354 53 54 4 44 1;
#columns
element nonlinearBeamColumn 1121 11 21 4 12 2;
element nonlinearBeamColumn 1222 12 22 4 12 2;
element nonlinearBeamColumn 1323 13 23 4 12 2;
element nonlinearBeamColumn 1424 14 24 4 12 2;
element nonlinearBeamColumn 2131 21 31 4 23 2;
element nonlinearBeamColumn 2232 22 32 4 23 2;
element nonlinearBeamColumn 2333 23 33 4 23 2;
element nonlinearBeamColumn 2434 24 34 4 23 2;
element nonlinearBeamColumn 3141 31 41 4 34 2;
element nonlinearBeamColumn 3242 32 42 4 34 2;
element nonlinearBeamColumn 3343 33 43 4 34 2;
element nonlinearBeamColumn 3444 34 44 4 34 2;
element nonlinearBeamColumn 4151 41 51 4 34 2;
element nonlinearBeamColumn 4252 42 52 4 34 2;
element nonlinearBeamColumn 4353 43 53 4 34 2;
element nonlinearBeamColumn 4454 44 54 4 34 2;
puts "Model Build"
set pi [expr 2*asin(1)]
set lambda [eigen 1]
set omega [expr pow($lambda,0.5)]
set Tperiod [expr 2*$pi /$omega]; # period (sec)
puts "$Tperiod"