3D RC frame Pushover

Forum for OpenSees users to post questions, comments, etc. on the use of the OpenSees interpreter, OpenSees.exe

Moderators: silvia, selimgunay, Moderators

Post Reply
farazmizani
Posts: 7
Joined: Tue Jun 27, 2017 4:06 am
Location: azad tabriz university

3D RC frame Pushover

Post by farazmizani » Sun Feb 07, 2021 3:52 am

Hello everyone,
I modeled a 3d RC structure based on example no.7 codes, the model is ok with eigen analysis and gravity analysis but when I wanna get a displacement-based pushover curve, this error shows up:
*BTW I tried every possible analysis command type but the result was the same*
if you have a solution for this please let me know what my problem is.
WARNING - ForceBeamColumn3d::update - failed to get compatible element forces & deformations for element: 143(dW: << -553.883, dW0: -12.9822)
WARNING - ForceBeamColumn3d::update - failed to get compatible element forces & deformations for element: 140(dW: << 287.511, dW0: 1.43825)
WARNING - ForceBeamColumn3d::update - failed to get compatible element forces & deformations for element: 143(dW: << 24.0492, dW0: -9.70741)
Domain::update - domain failed in update
DisplacementControl::update - model failed to update for new dU
WARNING NewtonRaphson::solveCurrentStep() -the Integrator failed in update()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at load factor 2.01189e+007
OpenSees > analyze failed, returned: -3 error flag
here is my script:
wipe ;
model basic -ndm 3 -ndf 6
###############################################
set g 9.81; #m/s2
set Hstory1 3.2; #m
set Hstory2 6.4; #m
set Hstory3 9.6; #m
set Ax 0.0; #m
set Bx 5.0; #m
set Cx 10.0; #m
set Dx 15.0; #m
set Ex 20.0; #m
set Fx 25.0; #m
set Gx 30.0; #m
set 1y 0.0; #m
set 2y 5.0; #m
set 3y 10.0; #m
set 4y 15.0; #m
set Ec 2500000.0; #tonf/m2
set ConcW 2.5; #tonf/m3
set DL 0.6; #tonf/m2
set LL 0.5; #tonf/m2
set Pu [expr (1.2*$DL)+(1.6*$LL)]; #tonf/m2
set Slab20 [expr 5.0*5.0*0.2*$ConcW]; #ton
puts "==============================================="
puts "****************Start Modeling*****************"
puts "==============================================="
##############################################################################
# START MODELING #
##############################################################################
#
#FIRST PART: NODES DEFINITION AND COORDINATES
#
#
#>>>>>node name digits formula:
#first digit: A=1,B=2,C=3,D=4,E=5,F=6,G=7
#second digit: 1y=1,2y=2,3y=3,4y=4
#third digit: Story1=1,Story2=2,Story3=3<<<<<<<<
#Base
# NAME X Y Z;COORDINATE NAME
node 110 $Ax $1y 0.0;#A1
node 120 $Ax $2y 0.0;#A2
node 130 $Ax $3y 0.0;#A3
node 140 $Ax $4y 0.0;#A4
node 210 $Bx $1y 0.0;#B1
node 220 $Bx $2y 0.0;#B2
node 230 $Bx $3y 0.0;#B3
node 240 $Bx $4y 0.0;#B4
node 310 $Cx $1y 0.0;#C1
node 320 $Cx $2y 0.0;#C2
node 330 $Cx $3y 0.0;#C3
node 340 $Cx $4y 0.0;#C4
node 410 $Dx $1y 0.0;#D1
node 420 $Dx $2y 0.0;#D2
node 430 $Dx $3y 0.0;#D3
node 440 $Dx $4y 0.0;#D4
node 510 $Ex $1y 0.0;#E1
node 520 $Ex $2y 0.0;#E2
node 530 $Ex $3y 0.0;#E3
node 540 $Ex $4y 0.0;#E4
node 610 $Fx $1y 0.0;#F1
node 620 $Fx $2y 0.0;#F2
node 630 $Fx $3y 0.0;#F3
node 640 $Fx $4y 0.0;#F4
node 710 $Gx $1y 0.0;#G1
node 720 $Gx $2y 0.0;#G2
node 730 $Gx $3y 0.0;#G3
node 740 $Gx $4y 0.0;#G4
#STORY1
# NAME X Y Z ;COORDINATE NAME
node 111 $Ax $1y $Hstory1;#A1
node 121 $Ax $2y $Hstory1;#A2
node 131 $Ax $3y $Hstory1;#A3
node 141 $Ax $4y $Hstory1;#A4
node 211 $Bx $1y $Hstory1;#B1
node 221 $Bx $2y $Hstory1;#B2
node 231 $Bx $3y $Hstory1;#B3
node 241 $Bx $4y $Hstory1;#B4
node 311 $Cx $1y $Hstory1;#C1
node 321 $Cx $2y $Hstory1;#C2
node 331 $Cx $3y $Hstory1;#C3
node 341 $Cx $4y $Hstory1;#C4
node 411 $Dx $1y $Hstory1;#D1
node 421 $Dx $2y $Hstory1;#D2
node 431 $Dx $3y $Hstory1;#D3
node 441 $Dx $4y $Hstory1;#D4
node 511 $Ex $1y $Hstory1;#E1
node 521 $Ex $2y $Hstory1;#E2
node 531 $Ex $3y $Hstory1;#E3
node 541 $Ex $4y $Hstory1;#E4
node 611 $Fx $1y $Hstory1;#F1
node 621 $Fx $2y $Hstory1;#F2
node 631 $Fx $3y $Hstory1;#F3
node 641 $Fx $4y $Hstory1;#F4
node 711 $Gx $1y $Hstory1;#G1
node 721 $Gx $2y $Hstory1;#G2
node 731 $Gx $3y $Hstory1;#G3
node 741 $Gx $4y $Hstory1;#G4
#STORY2
# NAME X Y Z ;COORDINATE NAME
node 112 $Ax $1y $Hstory2;#A1
node 122 $Ax $2y $Hstory2;#A2
node 132 $Ax $3y $Hstory2;#A3
node 142 $Ax $4y $Hstory2;#A4
node 212 $Bx $1y $Hstory2;#B1
node 222 $Bx $2y $Hstory2;#B2
node 232 $Bx $3y $Hstory2;#B3
node 242 $Bx $4y $Hstory2;#B4
node 312 $Cx $1y $Hstory2;#C1
node 322 $Cx $2y $Hstory2;#C2
node 332 $Cx $3y $Hstory2;#C3
node 342 $Cx $4y $Hstory2;#C4
node 412 $Dx $1y $Hstory2;#D1
node 422 $Dx $2y $Hstory2;#D2
node 432 $Dx $3y $Hstory2;#D3
node 442 $Dx $4y $Hstory2;#D4
node 512 $Ex $1y $Hstory2;#E1
node 522 $Ex $2y $Hstory2;#E2
node 532 $Ex $3y $Hstory2;#E3
node 542 $Ex $4y $Hstory2;#E4
node 612 $Fx $1y $Hstory2;#F1
node 622 $Fx $2y $Hstory2;#F2
node 632 $Fx $3y $Hstory2;#F3
node 642 $Fx $4y $Hstory2;#F4
node 712 $Gx $1y $Hstory2;#G1
node 722 $Gx $2y $Hstory2;#G2
node 732 $Gx $3y $Hstory2;#G3
node 742 $Gx $4y $Hstory2;#G4
#STORY3
# NAME X Y Z ;COORDINATE NAME
node 113 $Ax $1y $Hstory3;#A1
node 123 $Ax $2y $Hstory3;#A2
node 133 $Ax $3y $Hstory3;#A3
node 143 $Ax $4y $Hstory3;#A4
node 213 $Bx $1y $Hstory3;#B1
node 223 $Bx $2y $Hstory3;#B2
node 233 $Bx $3y $Hstory3;#B3
node 243 $Bx $4y $Hstory3;#B4
node 313 $Cx $1y $Hstory3;#C1
node 323 $Cx $2y $Hstory3;#C2
node 333 $Cx $3y $Hstory3;#C3
node 343 $Cx $4y $Hstory3;#C4
node 413 $Dx $1y $Hstory3;#D1
node 423 $Dx $2y $Hstory3;#D2
node 433 $Dx $3y $Hstory3;#D3
node 443 $Dx $4y $Hstory3;#D4
node 513 $Ex $1y $Hstory3;#E1
node 523 $Ex $2y $Hstory3;#E2
node 533 $Ex $3y $Hstory3;#E3
node 543 $Ex $4y $Hstory3;#E4
node 613 $Fx $1y $Hstory3;#F1
node 623 $Fx $2y $Hstory3;#F2
node 633 $Fx $3y $Hstory3;#F3
node 643 $Fx $4y $Hstory3;#F4
node 713 $Gx $1y $Hstory3;#G1
node 723 $Gx $2y $Hstory3;#G2
node 733 $Gx $3y $Hstory3;#G3
node 743 $Gx $4y $Hstory3;#G4
# define Rigid Floor Diaphragm
node 1111 2.5 2.5 3.2
node 1112 2.5 7.5 3.2
node 1113 2.5 12.5 3.2
node 1114 7.5 2.5 3.2
node 1115 7.5 7.5 3.2
node 1116 7.5 12.5 3.2
node 1117 12.5 2.5 3.2
node 1118 12.5 7.5 3.2
node 1119 12.5 12.5 3.2
node 1120 17.5 2.5 3.2
node 1121 17.5 7.5 3.2
node 1122 17.5 12.5 3.2
node 1123 22.5 2.5 3.2
node 1124 22.5 7.5 3.2
node 1125 22.5 12.5 3.2
node 1126 27.5 2.5 3.2
node 1127 27.5 7.5 3.2
node 1128 27.5 12.5 3.2
node 1129 2.5 2.5 6.4
node 1130 2.5 7.5 6.4
node 1131 2.5 12.5 6.4
node 1132 7.5 2.5 6.4
node 1133 7.5 7.5 6.4
node 1134 7.5 12.5 6.4
node 1135 12.5 2.5 6.4
node 1136 12.5 7.5 6.4
node 1137 12.5 12.5 6.4
node 1138 17.5 2.5 6.4
node 1139 17.5 7.5 6.4
node 1140 17.5 12.5 6.4
node 1141 22.5 2.5 6.4
node 1142 22.5 7.5 6.4
node 1143 22.5 12.5 6.4
node 1144 27.5 2.5 6.4
node 1145 27.5 7.5 6.4
node 1146 27.5 12.5 6.4
node 1147 2.5 2.5 9.6
node 1148 2.5 7.5 9.6
node 1149 2.5 12.5 9.6
node 1150 7.5 2.5 9.6
node 1151 7.5 7.5 9.6
node 1152 7.5 12.5 9.6
node 1153 12.5 2.5 9.6
node 1154 12.5 7.5 9.6
node 1155 12.5 12.5 9.6
node 1156 17.5 2.5 9.6
node 1157 17.5 7.5 9.6
node 1158 17.5 12.5 9.6
node 1159 22.5 2.5 9.6
node 1160 22.5 7.5 9.6
node 1161 22.5 12.5 9.6
node 1162 27.5 2.5 9.6
node 1163 27.5 7.5 9.6
node 1164 27.5 12.5 9.6
puts ">>>>nodes coordinate definition is complete"
#
#
#SECOND PART: MASS CALCULATION AND ASSIGNMENT
#
set HCol 0.5
set BCol 0.5
set LCol 3.2
set LBeam 5.0
set HBeam 0.65
set BBeam 0.45
set Tslab 0.2;
set Lslab 5.0;
set DLfactor 1.2;# scale dead load up a little
set LLfactor 1.6;# scale live load up a little
set GammaConcrete 2.5;# Reinforced-Concrete weight density (weight per volume)
set Qslab [expr $GammaConcrete*$Tslab*$Lslab];
set QdlCol [expr $GammaConcrete*$HCol*$BCol];# self weight of Column, weight per length
set QBeam [expr $GammaConcrete*$HBeam*$BBeam];# self weight of Beam, weight per length
set QdlBeam [expr $Qslab + $QBeam]
set WeightCol [expr $QdlCol*$LCol]
set WeightBeam [expr $QdlBeam*$LBeam]
set FloorWeight3 [expr 28.0*$WeightCol +45.0*$WeightBeam]
set FloorWeight2 [expr 28.0*$WeightCol +45.0*$WeightBeam]
set FloorWeight1 [expr 28.0*$WeightCol +45.0*$WeightBeam]
set Mmid [expr ($WeightCol/2.0+$WeightCol/2.0+$WeightBeam/2.0+$WeightBeam/2.0+$WeightBeam/2.0+$WeightBeam/2.0)/($g)]
set MmidC [expr ($WeightCol/2.0+$WeightCol/2.0+$WeightBeam/2.0+$WeightBeam/2.0)/($g)]
set MmidS [expr ($WeightCol/2.0+$WeightCol/2.0+$WeightBeam/2.0+$WeightBeam/2.0+$WeightBeam/2.0)/($g)]
set Mtop [expr ($WeightCol/2.0+$WeightBeam/2.0+$WeightBeam/2.0+$WeightBeam/2.0+$WeightBeam/2.0)/($g)]
set MtopC [expr ($WeightCol/2.0+$WeightBeam/2.0+$WeightBeam/2.0)/($g)]
set MtopS [expr ($WeightCol/2.0+$WeightBeam/2.0+$WeightBeam/2.0+$WeightBeam/2.0)/($g)]
set WeightTotal [expr $FloorWeight1+$FloorWeight2+$FloorWeight3];# total building weight
set MassTotal [expr $WeightTotal/$g];# total building mass
set sumWiHi [expr $FloorWeight1*3.20+$FloorWeight2*3.20+$FloorWeight3*3.20];# sum of storey weight times height, for lateral-load distribution
set WiHi1 [expr $FloorWeight1*3.20];# storey weight times height, for lateral-load distribution
set WiHi2 [expr $FloorWeight2*3.20];# storey weight times height, for lateral-load distribution
set WiHi3 [expr $FloorWeight3*3.20];# storey weight times height, for lateral-load distribution
set F1 [expr $WiHi1/($sumWiHi*$WeightTotal*18.0)];# lateral load at level
set F2 [expr $WiHi2/($sumWiHi*$WeightTotal*18.0)];# lateral load at level
set F3 [expr $WiHi3/($sumWiHi*$WeightTotal*18.0)];# lateral load at level
#
mass 111 $MmidC $MmidC 0 0 0 0; #A1
mass 141 $MmidC $MmidC 0 0 0 0; #A4
mass 711 $MmidC $MmidC 0 0 0 0; #G1
mass 741 $MmidC $MmidC 0 0 0 0; #G4
mass 121 $MmidS $MmidS 0 0 0 0; #A2
mass 131 $MmidS $MmidS 0 0 0 0; #A3
mass 721 $MmidS $MmidS 0 0 0 0; #G2
mass 731 $MmidS $MmidS 0 0 0 0; #G3
mass 211 $MmidS $MmidS 0 0 0 0; #B1
mass 241 $MmidS $MmidS 0 0 0 0; #B4
mass 311 $MmidS $MmidS 0 0 0 0; #C1
mass 341 $MmidS $MmidS 0 0 0 0; #C4
mass 411 $MmidS $MmidS 0 0 0 0; #D1
mass 441 $MmidS $MmidS 0 0 0 0; #D4
mass 511 $MmidS $MmidS 0 0 0 0; #E1
mass 541 $MmidS $MmidS 0 0 0 0; #E4
mass 611 $MmidS $MmidS 0 0 0 0; #F1
mass 641 $MmidS $MmidS 0 0 0 0; #F4
mass 221 $Mmid $Mmid 0 0 0 0; #B2
mass 231 $Mmid $Mmid 0 0 0 0; #B3
mass 321 $Mmid $Mmid 0 0 0 0; #C2
mass 331 $Mmid $Mmid 0 0 0 0; #C3
mass 421 $Mmid $Mmid 0 0 0 0; #D2
mass 431 $Mmid $Mmid 0 0 0 0; #D3
mass 521 $Mmid $Mmid 0 0 0 0; #E2
mass 531 $Mmid $Mmid 0 0 0 0; #E3
mass 621 $Mmid $Mmid 0 0 0 0; #F2
mass 631 $Mmid $Mmid 0 0 0 0; #F3
mass 112 $MmidC $MmidC 0 0 0 0; #A1
mass 142 $MmidC $MmidC 0 0 0 0; #A4
mass 712 $MmidC $MmidC 0 0 0 0; #G1
mass 742 $MmidC $MmidC 0 0 0 0; #G4
mass 122 $MmidS $MmidS 0 0 0 0; #A2
mass 132 $MmidS $MmidS 0 0 0 0; #A3
mass 722 $MmidS $MmidS 0 0 0 0; #G2
mass 732 $MmidS $MmidS 0 0 0 0; #G3
mass 212 $MmidS $MmidS 0 0 0 0; #B1
mass 242 $MmidS $MmidS 0 0 0 0; #B4
mass 312 $MmidS $MmidS 0 0 0 0; #C1
mass 342 $MmidS $MmidS 0 0 0 0; #C4
mass 412 $MmidS $MmidS 0 0 0 0; #D1
mass 442 $MmidS $MmidS 0 0 0 0; #D4
mass 512 $MmidS $MmidS 0 0 0 0; #E1
mass 542 $MmidS $MmidS 0 0 0 0; #E4
mass 612 $MmidS $MmidS 0 0 0 0; #F1
mass 642 $MmidS $MmidS 0 0 0 0; #F4
mass 222 $Mmid $Mmid 0 0 0 0; #B2
mass 232 $Mmid $Mmid 0 0 0 0; #B3
mass 322 $Mmid $Mmid 0 0 0 0; #C2
mass 332 $Mmid $Mmid 0 0 0 0; #C3
mass 422 $Mmid $Mmid 0 0 0 0; #D2
mass 432 $Mmid $Mmid 0 0 0 0; #D3
mass 522 $Mmid $Mmid 0 0 0 0; #E2
mass 532 $Mmid $Mmid 0 0 0 0; #E3
mass 622 $Mmid $Mmid 0 0 0 0; #F2
mass 632 $Mmid $Mmid 0 0 0 0; #F3
mass 113 $MtopC $MtopC 0 0 0 0; #A1
mass 143 $MtopC $MtopC 0 0 0 0; #A4
mass 713 $MtopC $MtopC 0 0 0 0; #G1
mass 743 $MtopC $MtopC 0 0 0 0; #G4
mass 123 $MtopS $MtopS 0 0 0 0; #A2
mass 133 $MtopS $MtopS 0 0 0 0; #A3
mass 723 $MtopS $MtopS 0 0 0 0; #G2
mass 733 $MtopS $MtopS 0 0 0 0; #G3
mass 213 $MtopS $MtopS 0 0 0 0; #B1
mass 243 $MtopS $MtopS 0 0 0 0; #B4
mass 313 $MtopS $MtopS 0 0 0 0; #C1
mass 343 $MtopS $MtopS 0 0 0 0; #C4
mass 413 $MtopS $MtopS 0 0 0 0; #D1
mass 443 $MtopS $MtopS 0 0 0 0; #D4
mass 513 $MtopS $MtopS 0 0 0 0; #E1
mass 543 $MtopS $MtopS 0 0 0 0; #E4
mass 613 $MtopS $MtopS 0 0 0 0; #F1
mass 643 $MtopS $MtopS 0 0 0 0; #F4
mass 223 $Mtop $Mtop 0 0 0 0; #B2
mass 233 $Mtop $Mtop 0 0 0 0; #B3
mass 323 $Mtop $Mtop 0 0 0 0; #C2
mass 333 $Mtop $Mtop 0 0 0 0; #C3
mass 423 $Mtop $Mtop 0 0 0 0; #D2
mass 433 $Mtop $Mtop 0 0 0 0; #D3
mass 523 $Mtop $Mtop 0 0 0 0; #E2
mass 533 $Mtop $Mtop 0 0 0 0; #E3
mass 623 $Mtop $Mtop 0 0 0 0; #F2
mass 633 $Mtop $Mtop 0 0 0 0; #F3
mass 1111 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1112 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1113 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1114 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1115 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1116 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1117 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1118 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1119 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1120 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1121 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1122 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1123 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1124 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1125 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1126 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1127 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1128 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1129 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1130 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1131 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1132 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1133 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1134 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1135 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1136 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1137 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1138 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1139 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1140 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1141 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1142 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1143 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1144 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1145 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1146 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1147 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1148 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1149 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1150 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1151 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1152 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1153 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1154 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1155 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1156 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1157 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1158 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1159 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1160 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1161 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1162 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1163 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
mass 1164 [expr $Pu*5.0*5.0] [expr $Pu*5.0*5.0] 0 0 0 0
###################################################################
#
puts ">>>>nodes mass definition is complete"
##########################################################################
# BOUNDAREY CONDITIOAN
##########################################################################
#
fix 110 1 1 1 1 1 1;#A1
fix 120 1 1 1 1 1 1;#A2
fix 130 1 1 1 1 1 1;#A3
fix 140 1 1 1 1 1 1;#A4
fix 210 1 1 1 1 1 1;#B1
fix 220 1 1 1 1 1 1;#B2
fix 230 1 1 1 1 1 1;#B3
fix 240 1 1 1 1 1 1;#B4
fix 310 1 1 1 1 1 1;#C1
fix 320 1 1 1 1 1 1;#C2
fix 330 1 1 1 1 1 1;#C3
fix 340 1 1 1 1 1 1;#C4
fix 410 1 1 1 1 1 1;#D1
fix 420 1 1 1 1 1 1;#D2
fix 430 1 1 1 1 1 1;#D3
fix 440 1 1 1 1 1 1;#D4
fix 510 1 1 1 1 1 1;#E1
fix 520 1 1 1 1 1 1;#E2
fix 530 1 1 1 1 1 1;#E3
fix 540 1 1 1 1 1 1;#E4
fix 610 1 1 1 1 1 1;#F1
fix 620 1 1 1 1 1 1;#F2
fix 630 1 1 1 1 1 1;#F3
fix 640 1 1 1 1 1 1;#F4
fix 710 1 1 1 1 1 1;#G1
fix 720 1 1 1 1 1 1;#G2
fix 730 1 1 1 1 1 1;#G3
fix 740 1 1 1 1 1 1;#G4
fix 1111 0 0 1 1 1 0
fix 1112 0 0 1 1 1 0
fix 1113 0 0 1 1 1 0
fix 1114 0 0 1 1 1 0
fix 1115 0 0 1 1 1 0
fix 1116 0 0 1 1 1 0
fix 1117 0 0 1 1 1 0
fix 1118 0 0 1 1 1 0
fix 1119 0 0 1 1 1 0
fix 1120 0 0 1 1 1 0
fix 1121 0 0 1 1 1 0
fix 1122 0 0 1 1 1 0
fix 1123 0 0 1 1 1 0
fix 1124 0 0 1 1 1 0
fix 1125 0 0 1 1 1 0
fix 1126 0 0 1 1 1 0
fix 1127 0 0 1 1 1 0
fix 1128 0 0 1 1 1 0
fix 1129 0 0 1 1 1 0
fix 1130 0 0 1 1 1 0
fix 1131 0 0 1 1 1 0
fix 1132 0 0 1 1 1 0
fix 1133 0 0 1 1 1 0
fix 1134 0 0 1 1 1 0
fix 1135 0 0 1 1 1 0
fix 1136 0 0 1 1 1 0
fix 1137 0 0 1 1 1 0
fix 1138 0 0 1 1 1 0
fix 1139 0 0 1 1 1 0
fix 1140 0 0 1 1 1 0
fix 1141 0 0 1 1 1 0
fix 1142 0 0 1 1 1 0
fix 1143 0 0 1 1 1 0
fix 1144 0 0 1 1 1 0
fix 1145 0 0 1 1 1 0
fix 1146 0 0 1 1 1 0
fix 1147 0 0 1 1 1 0
fix 1148 0 0 1 1 1 0
fix 1149 0 0 1 1 1 0
fix 1150 0 0 1 1 1 0
fix 1151 0 0 1 1 1 0
fix 1152 0 0 1 1 1 0
fix 1153 0 0 1 1 1 0
fix 1154 0 0 1 1 1 0
fix 1155 0 0 1 1 1 0
fix 1156 0 0 1 1 1 0
fix 1157 0 0 1 1 1 0
fix 1158 0 0 1 1 1 0
fix 1159 0 0 1 1 1 0
fix 1160 0 0 1 1 1 0
fix 1161 0 0 1 1 1 0
fix 1162 0 0 1 1 1 0
fix 1163 0 0 1 1 1 0
fix 1164 0 0 1 1 1 0
fixZ 0.0 1 1 1 1 1 1
##########################################################################
# MATERIAL PROPERTIES
##########################################################################
#
set IDconcCore 1; # material ID tag -- confined core concrete
set IDconcCover 2; # material ID tag -- unconfined cover concrete
set IDreinf 3; # material ID tag -- reinforcement
#concrete properties
#unconfined
set wc $ConcW
set fc 2500.0; #tonf/m2 #unconfined compression strength
set eco 0.0023
set ecp 0.005
#confined
set fcc [expr 2.0*$fc] ; #f'cc is the maximum compressive stress of concrete
set ecc 0.01546 ; #ecc is the strain at the maximum compressive stress
set ecu [expr 5.0*$ecc];
uniaxialMaterial Concrete04 $IDconcCover -$fc -$eco -$ecp $Ec
uniaxialMaterial Concrete04 $IDconcCore -$fcc -$ecc -$ecu $Ec
#steel bars
set Fy 24000.0; #tonf/m2
set Es 20000000.0; #tonf/m2
set b 0.013
set R0 18.0
set cR1 0.925
set cR2 0.15
uniaxialMaterial Steel02 $IDreinf $Fy $Es $b $R0 $cR1 $cR2
#-------------------------------------------------------------------------------------------------------------------------------------------
# ------------------------define Rigid Diaphram, dof 2 is normal to floor
set perpDirn 3;
rigidDiaphragm $perpDirn 1111 111 121 221 211
rigidDiaphragm $perpDirn 1112 121 131 231 221
rigidDiaphragm $perpDirn 1113 131 141 241 231
rigidDiaphragm $perpDirn 1114 211 221 321 311
rigidDiaphragm $perpDirn 1115 221 231 331 321
rigidDiaphragm $perpDirn 1116 231 241 341 331
rigidDiaphragm $perpDirn 1117 311 321 421 411
rigidDiaphragm $perpDirn 1118 321 331 431 421
rigidDiaphragm $perpDirn 1119 331 341 441 431
rigidDiaphragm $perpDirn 1120 411 421 521 511
rigidDiaphragm $perpDirn 1121 421 431 531 521
rigidDiaphragm $perpDirn 1122 431 441 541 531
rigidDiaphragm $perpDirn 1123 511 521 621 611
rigidDiaphragm $perpDirn 1124 521 531 631 621
rigidDiaphragm $perpDirn 1125 531 541 641 631
rigidDiaphragm $perpDirn 1126 611 621 721 711
rigidDiaphragm $perpDirn 1127 621 631 731 721
rigidDiaphragm $perpDirn 1128 631 641 741 731
rigidDiaphragm $perpDirn 1129 112 122 222 212
rigidDiaphragm $perpDirn 1130 122 132 232 222
rigidDiaphragm $perpDirn 1131 132 142 242 232
rigidDiaphragm $perpDirn 1132 212 222 322 312
rigidDiaphragm $perpDirn 1133 222 232 332 322
rigidDiaphragm $perpDirn 1134 232 242 342 332
rigidDiaphragm $perpDirn 1135 312 322 422 412
rigidDiaphragm $perpDirn 1136 322 332 432 422
rigidDiaphragm $perpDirn 1137 332 342 442 432
rigidDiaphragm $perpDirn 1138 412 422 522 512
rigidDiaphragm $perpDirn 1139 422 432 532 522
rigidDiaphragm $perpDirn 1140 432 442 542 532
rigidDiaphragm $perpDirn 1141 512 522 622 612
rigidDiaphragm $perpDirn 1142 522 532 632 622
rigidDiaphragm $perpDirn 1143 532 542 642 632
rigidDiaphragm $perpDirn 1144 612 622 722 712
rigidDiaphragm $perpDirn 1145 622 632 732 722
rigidDiaphragm $perpDirn 1146 632 642 742 732
rigidDiaphragm $perpDirn 1147 113 123 223 213
rigidDiaphragm $perpDirn 1148 123 133 233 223
rigidDiaphragm $perpDirn 1149 133 143 243 233
rigidDiaphragm $perpDirn 1150 213 223 323 313
rigidDiaphragm $perpDirn 1151 223 233 333 323
rigidDiaphragm $perpDirn 1152 233 243 343 333
rigidDiaphragm $perpDirn 1153 313 323 423 413
rigidDiaphragm $perpDirn 1154 323 333 433 423
rigidDiaphragm $perpDirn 1155 333 343 443 433
rigidDiaphragm $perpDirn 1156 413 423 523 513
rigidDiaphragm $perpDirn 1157 423 433 533 523
rigidDiaphragm $perpDirn 1158 433 443 543 533
rigidDiaphragm $perpDirn 1159 513 523 623 613
rigidDiaphragm $perpDirn 1160 523 533 633 623
rigidDiaphragm $perpDirn 1161 533 543 643 633
rigidDiaphragm $perpDirn 1162 613 623 723 713
rigidDiaphragm $perpDirn 1163 623 633 733 723
rigidDiaphragm $perpDirn 1164 633 643 743 733
##########################################################################
# COLUMN SECTION DEFINITION
##########################################################################
#sec for 50/16/22
set numBarsTop703228 5; #number of uniformly-distributed longitudinal-reinforcement bars
set numBarsBot703228 5; #number of uniformly-distributed longitudinal-reinforcement bars
set numBarsInt703228 3; #number of uniformly-distributed longitudinal-reinforcement bars
set Barsize 0.022; #m
set barArea [expr $Barsize*$Barsize*3.140/4.0] ; #area of longitudinal-reinforcement bars
set Hsec 0.5; #depth m
set Bsec 0.5; #width m
set cover703228 0.04; #m
set coreY703228 [expr ($Hsec/2.0)-($cover703228+($Barsize/2.0))];
set coreZ703228 [expr ($Bsec/2.0)-($cover703228+($Barsize/2.0))];
set coverY703228 [expr $Hsec/2.0];
set coverZ703228 [expr $Bsec/2.0];
set numSubdivY703228 10;
set numSubdivZ703228 10;
set SecTag703228 10;
section Fiber $SecTag703228 {
patch rect $IDconcCore $numSubdivY703228 $numSubdivZ703228 -$coreY703228 $coreZ703228 $coreY703228 -$coreZ703228
patch quadr $IDconcCover 1 $numSubdivY703228 -$coverY703228 $coverZ703228 -$coreY703228 $coreZ703228 $coreY703228 $coreZ703228 $coverY703228 $coverZ703228
patch quadr $IDconcCover 1 $numSubdivY703228 -$coreY703228 -$coreZ703228 -$coverY703228 -$coverZ703228 $coverY703228 -$coverZ703228 $coreY703228 -$coreZ703228
patch quadr $IDconcCover $numSubdivZ703228 1 -$coverY703228 $coverZ703228 -$coverY703228 -$coverZ703228 -$coreY703228 -$coreZ703228 -$coreY703228 $coreZ703228
patch quadr $IDconcCover $numSubdivZ703228 1 $coreY703228 $coreZ703228 $coreY703228 -$coreZ703228 $coverY703228 -$coverZ703228 $coverY703228 $coverZ703228
layer straight $IDreinf $numBarsInt703228 $barArea -$coreY703228 $coreZ703228 $coreY703228 $coreZ703228; # intermediate skin reinf. +z
layer straight $IDreinf $numBarsInt703228 $barArea -$coreY703228 -$coreZ703228 $coreY703228 -$coreZ703228; # intermediate skin reinf. -z
layer straight $IDreinf $numBarsTop703228 $barArea $coreY703228 $coreZ703228 $coreY703228 -$coreZ703228; # top layer reinfocement
layer straight $IDreinf $numBarsBot703228 $barArea -$coreY703228 $coreZ703228 -$coreY703228 -$coreZ703228; # bottom layer reinforcement
}
set SecTagA703228 101010
set G 1086956.0;
set J [expr (($Hsec*$Bsec)/12.0)*(($Hsec*$Hsec)+($Bsec*$Bsec))] ;
set torsiontag 111111;
set Gc [expr $G*$J/3.20];
uniaxialMaterial Elastic $torsiontag $Gc
section Aggregator $SecTagA703228 $torsiontag T -section $SecTag703228
puts "Nonlinear Materials & Sections are Defined"
#
#
##########################################################################
# ELEMENT PROPERTIES
##########################################################################
#
#Beam elements story 1 vecxzX
set transfTagBx 1
geomTransf Linear $transfTagBx 1 0 0
set transfTagBy 2
geomTransf Linear $transfTagBy 0 -1 0
set Bbeam 0.65
set Hbeam 0.45
set A [expr $Hbeam*$Bbeam]
set J [expr 0.011297*1.0]
set Iy [expr 0.004936*1.0]
set Iz [expr 0.010298*1.0]
set massDens [expr $A*$ConcW]
set numIntgrPts 5
set eleTag1 1
set eleTag2 2
set eleTag3 3
set eleTag4 4
set eleTag5 5
set eleTag6 6
set eleTag7 7
set eleTag8 8
set eleTag9 9
set eleTag10 10
set eleTag11 11
set eleTag12 12
set eleTag13 13
set eleTag14 14
set eleTag15 15
set eleTag16 16
set eleTag17 17
set eleTag18 18
set eleTag19 19
set eleTag20 20
set eleTag21 21
set eleTag22 22
set eleTag23 23
set eleTag24 24
set eleTag25 25
set eleTag26 26
set eleTag27 27
set eleTag28 28
set eleTag29 29
set eleTag30 30
set eleTag31 31
set eleTag32 32
set eleTag33 33
set eleTag34 34
set eleTag35 35
set eleTag36 36
set eleTag37 37
set eleTag38 38
set eleTag39 39
set eleTag40 40
set eleTag41 41
set eleTag42 42
set eleTag43 43
set eleTag44 44
set eleTag45 45
set eleTag46 46
set eleTag47 47
set eleTag48 48
set eleTag49 49
set eleTag50 50
set eleTag51 51
set eleTag52 52
set eleTag53 53
set eleTag54 54
set eleTag55 55
set eleTag56 56
set eleTag57 57
set eleTag58 58
set eleTag59 59
set eleTag60 60
set eleTag61 61
set eleTag62 62
set eleTag63 63
set eleTag64 64
set eleTag65 65
set eleTag66 66
set eleTag67 67
set eleTag68 68
set eleTag69 69
set eleTag70 70
set eleTag71 71
set eleTag72 72
set eleTag73 73
set eleTag74 74
set eleTag75 75
set eleTag76 76
set eleTag77 77
set eleTag78 78
set eleTag79 79
set eleTag80 80
set eleTag81 81
set eleTag82 82
set eleTag83 83
set eleTag84 84
set eleTag85 85
set eleTag86 86
set eleTag87 87
set eleTag88 88
set eleTag89 89
set eleTag90 90
set eleTag91 91
set eleTag92 92
set eleTag93 93
set eleTag94 94
set eleTag95 95
set eleTag96 96
set eleTag97 97
set eleTag98 98
set eleTag99 99
set eleTag100 100
set eleTag101 101
set eleTag102 102
set eleTag103 103
set eleTag104 104
set eleTag105 105
set eleTag106 106
set eleTag107 107
set eleTag108 108
set eleTag109 109
set eleTag110 110
set eleTag111 111
set eleTag112 112
set eleTag113 113
set eleTag114 114
set eleTag115 115
set eleTag116 116
set eleTag117 117
set eleTag118 118
set eleTag119 119
set eleTag120 120
set eleTag121 121
set eleTag122 122
set eleTag123 123
set eleTag124 124
set eleTag125 125
set eleTag126 126
set eleTag127 127
set eleTag128 128
set eleTag129 129
set eleTag130 130
set eleTag131 131
set eleTag132 132
set eleTag133 133
set eleTag134 134
set eleTag135 135
element elasticBeamColumn $eleTag1 111 121 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag2 121 131 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag3 131 141 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag4 211 221 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag5 221 231 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag6 231 241 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag7 311 321 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag8 321 331 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag9 331 341 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag10 411 421 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag11 421 431 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag12 431 441 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag13 511 521 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag14 521 531 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag15 531 541 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag16 611 621 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag17 621 631 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag18 631 641 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag19 711 721 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag20 721 731 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag21 731 741 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag22 112 122 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag23 122 132 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag24 132 142 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag25 212 222 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag26 222 232 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag27 232 242 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag28 312 322 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag29 322 332 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag30 332 342 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag31 412 422 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag32 422 432 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag33 432 442 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag34 512 522 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag35 522 532 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag36 532 542 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag37 612 622 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag38 622 632 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag39 632 642 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag40 712 722 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag41 722 732 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag42 732 742 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag43 113 123 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag44 123 133 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag45 133 143 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag46 213 223 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag47 223 233 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag48 233 243 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag49 313 323 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag50 323 333 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag51 333 343 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag52 413 423 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag53 423 433 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag54 433 443 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag55 513 523 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag56 523 533 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag57 533 543 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag58 613 623 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag59 623 633 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag60 633 643 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag61 713 723 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag62 723 733 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag63 733 743 $A $Ec $Gc $J $Iy $Iz $transfTagBx -mass $massDens
element elasticBeamColumn $eleTag64 111 211 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag65 121 221 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag66 131 231 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag67 141 241 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag68 211 311 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag69 221 321 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag70 231 331 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag71 241 341 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag72 311 411 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag73 321 421 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag74 331 431 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag75 341 441 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag76 411 511 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag77 421 521 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag78 431 531 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag79 441 541 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag80 511 611 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag81 521 621 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag82 531 631 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag83 541 641 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag84 611 711 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag85 621 721 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag86 631 731 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag87 641 741 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag88 112 212 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag89 122 222 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag90 132 232 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag91 142 242 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag92 212 312 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag93 222 322 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag94 232 332 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag95 242 342 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag96 312 412 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag97 322 422 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag98 332 432 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag99 342 442 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag100 412 512 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag101 422 522 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag102 432 532 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag103 442 542 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag104 512 612 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag105 522 622 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag106 532 632 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag107 542 642 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag108 612 712 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag109 622 722 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag110 632 732 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag111 642 742 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag112 113 213 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag113 123 223 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag114 133 233 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag115 143 243 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag116 213 313 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag117 223 323 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag118 233 333 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag119 243 343 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag120 313 413 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag121 323 423 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag122 333 433 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag123 343 443 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag124 413 513 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag125 423 523 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag126 433 533 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag127 443 543 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag128 513 613 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag129 523 623 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag130 533 633 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag131 543 643 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag132 613 713 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag133 623 723 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag134 633 733 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
element elasticBeamColumn $eleTag135 643 743 $A $Ec $Gc $J $Iy $Iz $transfTagBy -mass $massDens
#Column elements vecxzX
set ColTransfTag 3
geomTransf PDelta $ColTransfTag 0 1 0
#####################################
# A-A #
#####################################
set eleTag136 136
set eleTag137 137
set eleTag138 138
set eleTag139 139
set eleTag140 140
set eleTag141 141
set eleTag142 142
set eleTag143 143
set eleTag144 144
set eleTag145 145
set eleTag146 146
set eleTag147 147
set eleTag148 148
set eleTag149 149
set eleTag150 150
set eleTag151 151
set eleTag152 152
set eleTag153 153
set eleTag154 154
set eleTag155 155
set eleTag156 156
set eleTag157 157
set eleTag158 158
set eleTag159 159
set eleTag160 160
set eleTag161 161
set eleTag162 162
set eleTag163 163
set eleTag164 164
set eleTag165 165
set eleTag166 166
set eleTag167 167
set eleTag168 168
set eleTag169 169
set eleTag170 170
set eleTag171 171
set eleTag172 172
set eleTag173 173
set eleTag174 174
set eleTag175 175
set eleTag176 176
set eleTag177 177
set eleTag178 178
set eleTag179 179
set eleTag180 180
set eleTag181 181
set eleTag182 182
set eleTag183 183
set eleTag184 184
set eleTag185 185
set eleTag186 186
set eleTag187 187
set eleTag188 188
set eleTag189 189
set eleTag190 190
set eleTag191 191
set eleTag192 192
set eleTag193 193
set eleTag194 194
set eleTag195 195
set eleTag196 196
set eleTag197 197
set eleTag198 198
set eleTag199 199
set eleTag200 200
set eleTag201 201
set eleTag202 202
set eleTag203 203
set eleTag204 204
set eleTag205 205
set eleTag206 206
set eleTag207 207
set eleTag208 208
set eleTag209 209
set eleTag210 210
set eleTag211 211
set eleTag212 212
set eleTag213 213
set eleTag214 214
set eleTag215 215
set eleTag216 216
set eleTag217 217
set eleTag218 218
set eleTag219 219
element nonlinearBeamColumn $eleTag136 110 111 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag137 111 112 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag138 112 113 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag139 120 121 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag140 121 122 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag141 122 123 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag142 130 131 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag143 131 132 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag144 132 133 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag145 140 141 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag146 141 142 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag147 142 143 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag148 210 211 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag149 211 212 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag150 212 213 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag151 220 221 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag152 221 222 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag153 222 223 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag154 230 231 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag155 231 232 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag156 232 233 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag157 240 241 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag158 241 242 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag159 242 243 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag160 310 311 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag161 311 312 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag162 312 313 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag163 320 321 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag164 321 322 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag165 322 323 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag166 330 331 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag167 331 332 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag168 332 333 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag169 340 341 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag170 341 342 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag171 342 343 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag172 410 411 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag173 411 412 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag174 412 413 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag175 420 421 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag176 421 422 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag177 422 423 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag178 430 431 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag179 431 432 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag180 432 433 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag181 440 441 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag182 441 442 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag183 442 443 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag184 510 511 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag185 511 512 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag186 512 513 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag187 520 521 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag188 521 522 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag189 522 523 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag190 530 531 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag191 531 532 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag192 532 533 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag193 540 541 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag194 541 542 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag195 542 543 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag196 610 611 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag197 611 612 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag198 612 613 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag199 620 621 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag200 621 622 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag201 622 623 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag202 630 631 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag203 631 632 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag204 632 633 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag205 640 641 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag206 641 642 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag207 642 643 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag208 710 711 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag209 711 712 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag210 712 713 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag211 720 721 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag212 721 722 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag213 722 723 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag214 730 731 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag215 731 732 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag216 732 733 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag217 740 741 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag218 741 742 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
element nonlinearBeamColumn $eleTag219 742 743 $numIntgrPts $SecTagA703228 $ColTransfTag ;#
puts "3D Model Section and Elements are Defined"
###################################################################################
###################################################################################
file mkdir Delshad
#
recorder Drift -xml Delshad/ColumneledriftX.txt -time -iNode 110 120 130 140 -jNode 111 121 131 141 -dof 1 -perpDirn 3
recorder Node -xml Delshad/node110reaction.txt -time -node 740 -dof 1 2 3 reaction
recorder Node -xml Delshad/node113disp.txt -time -node 743 -dof 1 2 3 disp
#
# Gravity-analysis parameters -- load-controlled static analysis
set Tol 1.0e-2 ; #convergence tolerance for test
constraints Transformation ; #how it handles boundary conditions
numberer RCM ; #renumber dof's to minimize band-width (optimization), if you want to
system BandGeneral ;
test EnergyIncr $Tol 6 ; #how to store and solve the system of equations in the analysis
algorithm Newton ; #use Newton's solution algorithm: updates tangent stiffness at every iteration
set NstepGravity 1000; #apply gravity in 10 steps
set DGravity [expr 1.0/$NstepGravity]; #first load increment;
integrator LoadControl $DGravity; #determine the next time step for an analysis
analysis Static; #define type of analysis static or transient
analyze $NstepGravity;
puts "Gravity Analysis Done!"
############################################################################
set eigen [eigen 3]
puts "
First Eigen Value is: [eigen 1]
First Mode Period is: [expr 2.0*3.1415/pow([eigen 1],0.5)]s
"
puts "Model Built"
loadConst -time 0.0
pattern Plain 200 Linear {
load 1111 $F1 0. 0 0 0 0
load 1112 $F1 0. 0 0 0 0
load 1113 $F1 0. 0 0 0 0
load 1114 $F1 0. 0 0 0 0
load 1115 $F1 0. 0 0 0 0
load 1116 $F1 0. 0 0 0 0
load 1117 $F1 0. 0 0 0 0
load 1118 $F1 0. 0 0 0 0
load 1119 $F1 0. 0 0 0 0
load 1120 $F1 0. 0 0 0 0
load 1121 $F1 0. 0 0 0 0
load 1122 $F1 0. 0 0 0 0
load 1123 $F1 0. 0 0 0 0
load 1124 $F1 0. 0 0 0 0
load 1125 $F1 0. 0 0 0 0
load 1126 $F1 0. 0 0 0 0
load 1127 $F1 0. 0 0 0 0
load 1128 $F1 0. 0 0 0 0
load 1129 $F2 0. 0 0 0 0
load 1130 $F2 0. 0 0 0 0
load 1131 $F2 0. 0 0 0 0
load 1132 $F2 0. 0 0 0 0
load 1133 $F2 0. 0 0 0 0
load 1134 $F2 0. 0 0 0 0
load 1135 $F2 0. 0 0 0 0
load 1136 $F2 0. 0 0 0 0
load 1137 $F2 0. 0 0 0 0
load 1138 $F2 0. 0 0 0 0
load 1139 $F2 0. 0 0 0 0
load 1140 $F2 0. 0 0 0 0
load 1141 $F2 0. 0 0 0 0
load 1142 $F2 0. 0 0 0 0
load 1143 $F2 0. 0 0 0 0
load 1144 $F2 0. 0 0 0 0
load 1145 $F2 0. 0 0 0 0
load 1146 $F2 0. 0 0 0 0
load 1147 $F3 0. 0 0 0 0
load 1148 $F3 0. 0 0 0 0
load 1149 $F3 0. 0 0 0 0
load 1150 $F3 0. 0 0 0 0
load 1151 $F3 0. 0 0 0 0
load 1152 $F3 0. 0 0 0 0
load 1153 $F3 0. 0 0 0 0
load 1154 $F3 0. 0 0 0 0
load 1155 $F3 0. 0 0 0 0
load 1156 $F3 0. 0 0 0 0
load 1157 $F3 0. 0 0 0 0
load 1158 $F3 0. 0 0 0 0
load 1159 $F3 0. 0 0 0 0
load 1160 $F3 0. 0 0 0 0
load 1161 $F3 0. 0 0 0 0
load 1162 $F3 0. 0 0 0 0
load 1163 $F3 0. 0 0 0 0
load 1164 $F3 0. 0 0 0 0
}
constraints Transformation ; #how it handles boundary conditions
numberer RCM ; #renumber dof's to minimize band-width (optimization), if you want to
system BandGeneral ;
test NormDispIncr $Tol 6 ; #how to store and solve the system of equations in the analysis
algorithm Newton ; #use Newton's solution algorithm: updates tangent stiffness at every iteration
integrator DisplacementControl 743 1 $DGravity; #determine the next time step for an analysis
analysis Static; #define type of analysis static or transient
analyze $NstepGravity;






Thank you,
Faraz Mizani

selimgunay
Posts: 889
Joined: Mon Sep 09, 2013 8:50 pm
Location: University of California, Berkeley

Re: 3D RC frame Pushover

Post by selimgunay » Tue Feb 09, 2021 8:46 am

There is a potential convergence problem at that step. How does the pushover curve look like before this convergence issue happens? You can try other nonlinear solution algorithms, convergence tests and tolerances. Below is an example of this type of adaptive solution algorithm for nonlinear dynamic analysis. You can adapt it for the nonlinear static (pushover) case.

viewtopic.php?f=2&t=67951&p=118077#p118077

farazmizani
Posts: 7
Joined: Tue Jun 27, 2017 4:06 am
Location: azad tabriz university

Re: 3D RC frame Pushover

Post by farazmizani » Sun Feb 21, 2021 12:52 pm

Dear @selimgunay thanks for your guidance my displacement and node reaction is zero in static pushover before converge error and when I use nonlinear dynamic analysis I get converge error flag 3 and I have max displacement 0.0005 meter which is almost zero.
I tried all types of solutions in static pushover and nonlinear dynamic but converge error shows up again.

selimgunay
Posts: 889
Joined: Mon Sep 09, 2013 8:50 pm
Location: University of California, Berkeley

Re: 3D RC frame Pushover

Post by selimgunay » Mon Feb 22, 2021 9:50 am

In this case, there should be something that makes your model not stable. Did you check your eigenvalues, are they reasonable? What are your periods?

Post Reply