Simple Question
Moderators: silvia, selimgunay, Moderators
Simple Question
I am sure the answer is yes, but please confirm. When you apply a axial load to a beam element it applies the load through the centroid of the section correct? So if I have some arbitraily shaped section, and this section is applied to a cantilever column. If I apply a vertical axial load to the column, the top of the column should not move in either the local y or z direction, correct? Thanks
That cannot be right define a simple section and put the local 0,0 outside the section or at a corner and place a axial load on it, you don't get any moments in the section. and hte top node doesn't move. The only thing that makes sense is it is through the centroid of the section. I just want someone to confirm I am thinking about this correct.
# Rotation parameters
# theta (rotation angle) 0 degrees CW
set sintheta 0
set costheta 1
# Create reference points for transformation and section generation
# -----------------------
# Confined conrete reference points
set 1y -35.375
set 1z 44.375
set 2y -16.875
set 2z 44.375
set 3y -16.875
set 3z 39.625
set 4y -35.375
set 4z 39.625
set 5y 16.875
set 5z 44.375
set 6y 35.375
set 6z 44.375
set 7y 35.375
set 7z 39.625
set 8y 16.875
set 8z 39.625
set 9y -16.875
set 9z 44.0625
set 10y 16.875
set 10z 44.0625
set 11y 16.875
set 11z 39.9375
set 12y -16.875
set 12z 39.9375
set 13y -2.0625
set 13z 39.9375
set 14y 2.0625
set 14z 39.9375
set 15y 2.0625
set 15z -22.5625
set 16y -2.0625
set 16z -22.5625
set 17y -2.375
set 17z -22.5625
set 18y 2.375
set 18z -22.5625
set 19y 2.375
set 19z -44.375
set 20y -2.375
set 20z -44.375
#Unconfined concrete reference points
set 21y -36.0
set 21z 45.0
set 22y -16.875
set 22z 45.0
set 25y 16.875
set 25z 45.0
set 26y 36.0
set 26z 45.0
set 27y 36.0
set 27z 44.375
set 28y 36.0
set 28z 39.625
set 29y 36.0
set 29z 39.0
set 30y 16.875
set 30z 39.0
set 32y 2.0625
set 32z 39.0
set 33y 3.0
set 33z 39.0
set 34y 3.0
set 34z -22.5625
set 36y 3.0
set 36z -45.0
set 37y 2.375
set 37z -45.0
set 38y -2.375
set 38z -45.0
set 39y -3.0
set 39z -45.0
set 40y -3.0
set 40z -22.5625
set 42y -2.0625
set 42z 39.0
set 43y -3.0
set 43z 39.0
set 44y -16.875
set 44z 39.0
set 46y -36.0
set 46z 39.0
set 47y -36.0
set 47z 39.625
set 48y -36.0
set 48z 44.375
# Create steeel layer reference points
# ------------------------
set ay -34.875
set az 43.875
set by -27.875
set bz 43.875
set cy -24.375
set cz 43.875
set dy -17.375
set dz 43.875
set ey -31.375
set ez 43.9375
set fy -20.875
set fz 43.9375
set gy -34.875
set gz 40.125
set hy -27.875
set hz 40.125
set iy -24.375
set iz 40.125
set jy -17.375
set jz 40.125
set ky -31.375
set kz 40.0625
set ly -20.875
set lz 40.0625
set my 17.375
set mz 43.875
set ny 24.375
set nz 43.875
set oy 27.875
set oz 43.875
set py 34.875
set pz 43.875
set qy 20.875
set qz 43.9375
set ry 31.375
set rz 43.9375
set sy 17.375
set sz 40.125
set ty 24.375
set tz 40.125
set uy 27.875
set uz 40.125
set vy 34.875
set vz 40.125
set wy 20.875
set wz 40.0625
set xy 31.375
set xz 40.0625
set yy -9.00
set yz 43.6875
set zy 9.00
set zz 43.6875
set Ay -6.25
set Az 40.3125
set By 6.25
set Bz 40.3125
set Cy -1.6875
set Cz 36.125
set Dy -1.6875
set Dz -13.875
set Ey 1.6875
set Ez 36.125
set Fy 1.6875
set Fz -13.875
set Gy -1.875
set Gz -26.375
set Hy -1.875
set Hz -33.375
set Iy -1.875
set Iz -36.875
set Jy -1.875
set Jz -43.875
set Ky -1.9375
set Kz -29.875
set Ly -1.9375
set Lz -40.375
set My 1.875
set Mz -26.375
set Ny 1.875
set Nz -33.375
set Oy 1.875
set Oz -36.875
set Py 1.875
set Pz -43.875
set Qy 1.9375
set Qz -29.875
set Ry 1.9375
set Rz -40.375
set Sy -2.0625
set Sz -22.875
set Ty 2.0625
set Tz -22.875
set Uy -2.5
set Uz 43.6875
set Vy 2.5
set Vz 43.6875
# Create rotated points through geometric transformation
set 1yt [expr $1y*$costheta+$1z*$sintheta]
set 1zt [expr -$1y*$sintheta+$1z*$costheta]
set 2yt [expr $2y*$costheta+$2z*$sintheta]
set 2zt [expr -$2y*$sintheta+$2z*$costheta]
set 3yt [expr $3y*$costheta+$3z*$sintheta]
set 3zt [expr -$3y*$sintheta+$3z*$costheta]
set 4yt [expr $4y*$costheta+$4z*$sintheta]
set 4zt [expr -$4y*$sintheta+$4z*$costheta]
set 5yt [expr $5y*$costheta+$5z*$sintheta]
set 5zt [expr -$5y*$sintheta+$5z*$costheta]
set 6yt [expr $6y*$costheta+$6z*$sintheta]
set 6zt [expr -$6y*$sintheta+$6z*$costheta]
set 7yt [expr $7y*$costheta+$7z*$sintheta]
set 7zt [expr -$7y*$sintheta+$7z*$costheta]
set 8yt [expr $8y*$costheta+$8z*$sintheta]
set 8zt [expr -$8y*$sintheta+$8z*$costheta]
set 9yt [expr $9y*$costheta+$9z*$sintheta]
set 9zt [expr -$9y*$sintheta+$9z*$costheta]
set 10yt [expr $10y*$costheta+$10z*$sintheta]
set 10zt [expr -$10y*$sintheta+$10z*$costheta]
set 11yt [expr $11y*$costheta+$11z*$sintheta]
set 11zt [expr -$11y*$sintheta+$11z*$costheta]
set 12yt [expr $12y*$costheta+$12z*$sintheta]
set 12zt [expr -$12y*$sintheta+$12z*$costheta]
set 13yt [expr $13y*$costheta+$13z*$sintheta]
set 13zt [expr -$13y*$sintheta+$13z*$costheta]
set 14yt [expr $14y*$costheta+$14z*$sintheta]
set 14zt [expr -$14y*$sintheta+$14z*$costheta]
set 15yt [expr $15y*$costheta+$15z*$sintheta]
set 15zt [expr -$15y*$sintheta+$15z*$costheta]
set 16yt [expr $16y*$costheta+$16z*$sintheta]
set 16zt [expr -$16y*$sintheta+$16z*$costheta]
set 17yt [expr $17y*$costheta+$17z*$sintheta]
set 17zt [expr -$17y*$sintheta+$17z*$costheta]
set 18yt [expr $18y*$costheta+$18z*$sintheta]
set 18zt [expr -$18y*$sintheta+$18z*$costheta]
set 19yt [expr $19y*$costheta+$19z*$sintheta]
set 19zt [expr -$19y*$sintheta+$19z*$costheta]
set 20yt [expr $20y*$costheta+$20z*$sintheta]
set 20zt [expr -$20y*$sintheta+$20z*$costheta]
set 21yt [expr $21y*$costheta+$21z*$sintheta]
set 21zt [expr -$21y*$sintheta+$21z*$costheta]
set 22yt [expr $22y*$costheta+$22z*$sintheta]
set 22zt [expr -$22y*$sintheta+$22z*$costheta]
set 25yt [expr $25y*$costheta+$25z*$sintheta]
set 25zt [expr -$25y*$sintheta+$25z*$costheta]
set 26yt [expr $26y*$costheta+$26z*$sintheta]
set 26zt [expr -$26y*$sintheta+$26z*$costheta]
set 27yt [expr $27y*$costheta+$27z*$sintheta]
set 27zt [expr -$27y*$sintheta+$27z*$costheta]
set 28yt [expr $28y*$costheta+$28z*$sintheta]
set 28zt [expr -$28y*$sintheta+$28z*$costheta]
set 29yt [expr $29y*$costheta+$29z*$sintheta]
set 29zt [expr -$29y*$sintheta+$29z*$costheta]
set 30yt [expr $30y*$costheta+$30z*$sintheta]
set 30zt [expr -$30y*$sintheta+$30z*$costheta]
set 32yt [expr $32y*$costheta+$32z*$sintheta]
set 32zt [expr -$32y*$sintheta+$32z*$costheta]
set 33yt [expr $33y*$costheta+$33z*$sintheta]
set 33zt [expr -$33y*$sintheta+$33z*$costheta]
set 34yt [expr $34y*$costheta+$34z*$sintheta]
set 34zt [expr -$34y*$sintheta+$34z*$costheta]
set 36yt [expr $36y*$costheta+$36z*$sintheta]
set 36zt [expr -$36y*$sintheta+$36z*$costheta]
set 37yt [expr $37y*$costheta+$37z*$sintheta]
set 37zt [expr -$37y*$sintheta+$37z*$costheta]
set 38yt [expr $38y*$costheta+$38z*$sintheta]
set 38zt [expr -$38y*$sintheta+$38z*$costheta]
set 39yt [expr $39y*$costheta+$39z*$sintheta]
set 39zt [expr -$39y*$sintheta+$39z*$costheta]
set 40yt [expr $40y*$costheta+$40z*$sintheta]
set 40zt [expr -$40y*$sintheta+$40z*$costheta]
set 42yt [expr $42y*$costheta+$42z*$sintheta]
set 42zt [expr -$42y*$sintheta+$42z*$costheta]
set 43yt [expr $43y*$costheta+$43z*$sintheta]
set 43zt [expr -$43y*$sintheta+$43z*$costheta]
set 44yt [expr $44y*$costheta+$44z*$sintheta]
set 44zt [expr -$44y*$sintheta+$44z*$costheta]
set 46yt [expr $46y*$costheta+$46z*$sintheta]
set 46zt [expr -$46y*$sintheta+$46z*$costheta]
set 47yt [expr $47y*$costheta+$47z*$sintheta]
set 47zt [expr -$47y*$sintheta+$47z*$costheta]
set 48yt [expr $48y*$costheta+$48z*$sintheta]
set 48zt [expr -$48y*$sintheta+$48z*$costheta]
set ayt [expr $ay*$costheta+$az*$sintheta]
set azt [expr -$ay*$sintheta+$az*$costheta]
set byt [expr $by*$costheta+$bz*$sintheta]
set bzt [expr -$by*$sintheta+$bz*$costheta]
set cyt [expr $cy*$costheta+$cz*$sintheta]
set czt [expr -$cy*$sintheta+$cz*$costheta]
set dyt [expr $dy*$costheta+$dz*$sintheta]
set dzt [expr -$dy*$sintheta+$dz*$costheta]
set eyt [expr $ey*$costheta+$ez*$sintheta]
set ezt [expr -$ey*$sintheta+$ez*$costheta]
set fyt [expr $fy*$costheta+$fz*$sintheta]
set fzt [expr -$fy*$sintheta+$fz*$costheta]
set gyt [expr $gy*$costheta+$gz*$sintheta]
set gzt [expr -$gy*$sintheta+$gz*$costheta]
set hyt [expr $hy*$costheta+$hz*$sintheta]
set hzt [expr -$hy*$sintheta+$hz*$costheta]
set iyt [expr $iy*$costheta+$iz*$sintheta]
set izt [expr -$iy*$sintheta+$iz*$costheta]
set jyt [expr $jy*$costheta+$jz*$sintheta]
set jzt [expr -$jy*$sintheta+$jz*$costheta]
set kyt [expr $ky*$costheta+$kz*$sintheta]
set kzt [expr -$ky*$sintheta+$kz*$costheta]
set lyt [expr $ly*$costheta+$lz*$sintheta]
set lzt [expr -$ly*$sintheta+$lz*$costheta]
set myt [expr $my*$costheta+$mz*$sintheta]
set mzt [expr -$my*$sintheta+$mz*$costheta]
set nyt [expr $ny*$costheta+$nz*$sintheta]
set nzt [expr -$ny*$sintheta+$nz*$costheta]
set oyt [expr $oy*$costheta+$oz*$sintheta]
set ozt [expr -$oy*$sintheta+$oz*$costheta]
set pyt [expr $py*$costheta+$pz*$sintheta]
set pzt [expr -$py*$sintheta+$pz*$costheta]
set qyt [expr $qy*$costheta+$qz*$sintheta]
set qzt [expr -$qy*$sintheta+$qz*$costheta]
set ryt [expr $ry*$costheta+$rz*$sintheta]
set rzt [expr -$ry*$sintheta+$rz*$costheta]
set syt [expr $sy*$costheta+$sz*$sintheta]
set szt [expr -$sy*$sintheta+$sz*$costheta]
set tyt [expr $ty*$costheta+$tz*$sintheta]
set tzt [expr -$ty*$sintheta+$tz*$costheta]
set uyt [expr $uy*$costheta+$uz*$sintheta]
set uzt [expr -$uy*$sintheta+$uz*$costheta]
set vyt [expr $vy*$costheta+$vz*$sintheta]
set vzt [expr -$vy*$sintheta+$vz*$costheta]
set wyt [expr $wy*$costheta+$wz*$sintheta]
set wzt [expr -$wy*$sintheta+$wz*$costheta]
set xyt [expr $xy*$costheta+$xz*$sintheta]
set xzt [expr -$xy*$sintheta+$xz*$costheta]
set yyt [expr $yy*$costheta+$yz*$sintheta]
set yzt [expr -$yy*$sintheta+$yz*$costheta]
set zyt [expr $zy*$costheta+$zz*$sintheta]
set zzt [expr -$zy*$sintheta+$zz*$costheta]
set Ayt [expr $Ay*$costheta+$Az*$sintheta]
set Azt [expr -$Ay*$sintheta+$Az*$costheta]
set Byt [expr $By*$costheta+$Bz*$sintheta]
set Bzt [expr -$By*$sintheta+$Bz*$costheta]
set Cyt [expr $Cy*$costheta+$Cz*$sintheta]
set Czt [expr -$Cy*$sintheta+$Cz*$costheta]
set Dyt [expr $Dy*$costheta+$Dz*$sintheta]
set Dzt [expr -$Dy*$sintheta+$Dz*$costheta]
set Eyt [expr $Ey*$costheta+$Ez*$sintheta]
set Ezt [expr -$Ey*$sintheta+$Ez*$costheta]
set Fyt [expr $Fy*$costheta+$Fz*$sintheta]
set Fzt [expr -$Fy*$sintheta+$Fz*$costheta]
set Gyt [expr $Gy*$costheta+$Gz*$sintheta]
set Gzt [expr -$Gy*$sintheta+$Gz*$costheta]
set Hyt [expr $Hy*$costheta+$Hz*$sintheta]
set Hzt [expr -$Hy*$sintheta+$Hz*$costheta]
set Iyt [expr $Iy*$costheta+$Iz*$sintheta]
set Izt [expr -$Iy*$sintheta+$Iz*$costheta]
set Jyt [expr $Jy*$costheta+$Jz*$sintheta]
set Jzt [expr -$Jy*$sintheta+$Jz*$costheta]
set Kyt [expr $Ky*$costheta+$Kz*$sintheta]
set Kzt [expr -$Ky*$sintheta+$Kz*$costheta]
set Lyt [expr $Ly*$costheta+$Lz*$sintheta]
set Lzt [expr -$Ly*$sintheta+$Lz*$costheta]
set Myt [expr $My*$costheta+$Mz*$sintheta]
set Mzt [expr -$My*$sintheta+$Mz*$costheta]
set Nyt [expr $Ny*$costheta+$Nz*$sintheta]
set Nzt [expr -$Ny*$sintheta+$Nz*$costheta]
set Oyt [expr $Oy*$costheta+$Oz*$sintheta]
set Ozt [expr -$Oy*$sintheta+$Oz*$costheta]
set Pyt [expr $Py*$costheta+$Pz*$sintheta]
set Pzt [expr -$Py*$sintheta+$Pz*$costheta]
set Qyt [expr $Qy*$costheta+$Qz*$sintheta]
set Qzt [expr -$Qy*$sintheta+$Qz*$costheta]
set Ryt [expr $Ry*$costheta+$Rz*$sintheta]
set Rzt [expr -$Ry*$sintheta+$Rz*$costheta]
set Syt [expr $Sy*$costheta+$Sz*$sintheta]
set Szt [expr -$Sy*$sintheta+$Sz*$costheta]
set Tyt [expr $Ty*$costheta+$Tz*$sintheta]
set Tzt [expr -$Ty*$sintheta+$Tz*$costheta]
set Uyt [expr $Uy*$costheta+$Uz*$sintheta]
set Uzt [expr -$Uy*$sintheta+$Uz*$costheta]
set Vyt [expr $Vy*$costheta+$Vz*$sintheta]
set Vzt [expr -$Vy*$sintheta+$Vz*$costheta]
# ------------------------
# Start model generation
# ------------------------
model basic -ndm 3 -ndf 6
# Create nodes
node 1 0.0 0.0 0.0
node 2 288.0 0.0 0.0
node 3 312.0 0.0 0.0
#node 3 72.0 0.0 0.0
#node 4 144.0 0.0 0.0
#node 5 216.0 0.0 0.0
#node 6 288.0 0.0 0.0
# Fix all DOF
fix 1 1 1 1 1 1 1
fix 2 0 0 0 1 0 0
fix 3 0 0 0 1 0 0
# Create uniaxial material models for the different fibers
# Concrete Material Properties
# ------------------------------
# Confined Concrete (Numbering corresponds to drawing)
uniaxialMaterial Concrete03 1 -6.551 -0.005101 -4.678 -0.021 0.55 0.5303 0.002 .00004 2.5 0.003
#uniaxialMaterial Concrete01 1 -6.551 -0.005101 -4.678 -0.021
uniaxialMaterial Concrete03 2 -7.318 -0.006637 -5.557 -0.028 0.55 0.5303 0.002 .00004 2.5 0.003
#uniaxialMaterial Concrete01 2 -7.318 -0.006637 -5.557 -0.028
uniaxialMaterial Concrete03 3 -5.790 -0.003579 -3.729 -0.013 0.55 0.5303 0.002 .00004 2.5 0.003
#uniaxialMaterial Concrete01 3 -5.790 -0.003579 -3.729 -0.013
uniaxialMaterial Concrete03 4 -6.085 -0.004170 -4.106 -0.016 0.55 0.5303 0.002 .00004 2.5 0.003
#uniaxialMaterial Concrete01 4 -6.085 -0.004170 -4.106 -0.016
# Unconfined Cover Concrete
uniaxialMaterial Concrete03 5 -5.000 -0.003 0.00 -0.008 0.1 0.5303 0.002 .00004 2.5 0.003
#uniaxialMaterial Concrete01 5 -5.00 -0.003 -0.00 -0.008
# Steel Material Properties
# ------------------------------
# Material Properties for #6 bars
uniaxialMaterial Steel02 6 70.0 29000 0.01 18.5 .925 .15 0 5 0 5
# Material Properties for #5 bars
uniaxialMaterial Steel02 7 70.0 29000 0.01 18.5 .925 .15 0 5 0 5
# Material Properties for #3 bars
uniaxialMaterial Steel02 8 70.0 29000 0.01 18.5 .925 .15 0 5 0 5
# Bond Material Properties
# -------------------------------
# Material Properties for #6 bar
uniaxialMaterial StrPen01 9 0.01897 70 0.664 103 0.3 0.8 0.0 1.0 4.35 25.0
# Material Properties for #5 bar
uniaxialMaterial StrPen01 10 0.0169 70 0.593 103 0.3 0.8 0.0 1.0 4.35 25.0
# Material Properties for #3 bar
uniaxialMaterial StrPen01 11 0.0144 70 0.503 103 0.3 0.8 0.0 1.0 4.35 25.0
# Material Properties for interface concrete fibers
uniaxialMaterial Concrete03 12 -6.551 -0.005101 -6.651 -0.04 0.1 0.5303 0.002 .00004 2.5 0.003
uniaxialMaterial Concrete03 13 -7.318 -0.006637 -7.418 -0.04 0.1 0.5303 0.002 .00004 2.5 0.003
uniaxialMaterial Concrete03 14 -5.790 -0.003579 -5.900 -0.04 0.1 0.5303 0.002 .00004 2.5 0.003
uniaxialMaterial Concrete03 15 -6.085 -0.004170 -6.185 -0.04 0.1 0.5303 0.002 .00004 2.5 0.003
uniaxialMaterial Concrete03 16 -5.000 -0.003 -5.100 -0.04 0.1 0.5303 0.002 .00004 2.5 0.003
# Proof of Theory, softening the fange boundary elements
uniaxialMaterial ElasticPPGap 17 12000 63 0.0002
#uniaxialMaterial ElasticPPGap 17 5000 63 0.0
uniaxialMaterial Series 18 6 17
uniaxialMaterial Series 19 7 17
set GJ 1.0e20
# Define cross-section for analysis
# Fibers are approx. 0.25" by 0.25"
section Fiber 1 {
# Create Core Concrete Fibers
patch quad 1 65 17 $4yt $4zt $3yt $3zt $2yt $2zt $1yt $1zt
patch quad 1 65 17 $8yt $8zt $7yt $7zt $6yt $6zt $5yt $5zt
patch quad 3 150 13 $12yt $12zt $11yt $11zt $10yt $10zt $9yt $9zt
patch quad 4 13 273 $16yt $16zt $15yt $15zt $14yt $14zt $13yt $13zt
patch quad 2 17 65 $20yt $20zt $19yt $19zt $18yt $18zt $17yt $17zt
# Create Cover Concrete Fibers
patch quad 5 69 4 $48yt $48zt $2yt $2zt $22yt $22zt $21yt $21zt
patch quad 5 150 6 $9yt $9zt $10yt $10zt $25yt $25zt $22yt $22zt
patch quad 5 69 4 $5yt $5zt $27yt $27zt $26yt $26zt $25yt $25zt
patch quad 5 4 17 $7yt $7zt $28yt $28zt $27yt $27zt $6yt $6zt
patch quad 5 69 4 $30yt $30zt $29yt $29zt $28yt $28zt $8yt $8zt
patch quad 5 69 6 $32yt $32zt $30yt $30zt $11yt $11zt $14yt $14zt
patch quad 5 6 267 $15yt $15zt $34yt $34zt $33yt $33zt $32yt $32zt
patch quad 5 4 69 $37yt $37zt $36yt $36zt $34yt $34zt $18yt $18zt
patch quad 5 17 4 $38yt $38zt $37yt $37zt $19yt $19zt $20yt $20zt
patch quad 5 4 69 $39yt $39zt $38yt $38zt $17yt $17zt $40yt $40zt
patch quad 5 6 267 $40yt $40zt $16yt $16zt $42yt $42zt $43yt $43zt
patch quad 5 69 6 $44yt $44zt $42yt $42zt $13yt $13zt $12yt $12zt
patch quad 5 69 4 $46yt $46zt $44yt $44zt $3yt $3zt $47yt $47zt
patch quad 5 4 17 $47yt $47zt $4yt $4zt $1yt $1zt $48yt $48zt
# Create Reinforcing Steel Fibers
layer straight 6 2 0.44 $ayt $azt $byt $bzt
layer straight 6 2 0.44 $cyt $czt $dyt $dzt
layer straight 7 2 0.31 $eyt $ezt $fyt $fzt
layer straight 6 2 0.44 $gyt $gzt $hyt $hzt
layer straight 6 2 0.44 $iyt $izt $jyt $jzt
layer straight 7 2 0.31 $kyt $kzt $lyt $lzt
layer straight 6 2 0.44 $myt $mzt $nyt $nzt
layer straight 6 2 0.44 $oyt $ozt $pyt $pzt
layer straight 7 2 0.31 $qyt $qzt $ryt $rzt
layer straight 6 2 0.44 $syt $szt $tyt $tzt
layer straight 6 2 0.44 $uyt $uzt $vyt $vzt
layer straight 7 2 0.31 $wyt $wzt $xyt $xzt
layer straight 8 2 0.11 $yyt $yzt $zyt $zzt
layer straight 8 2 0.11 $Ayt $Azt $Byt $Bzt
layer straight 8 5 0.11 $Cyt $Czt $Dyt $Dzt
layer straight 8 5 0.11 $Eyt $Ezt $Fyt $Fzt
layer straight 6 2 0.44 $Gyt $Gzt $Hyt $Hzt
layer straight 6 2 0.44 $Iyt $Izt $Jyt $Jzt
layer straight 7 2 0.31 $Kyt $Kzt $Lyt $Lzt
layer straight 6 2 0.44 $Myt $Mzt $Nyt $Nzt
layer straight 6 2 0.44 $Oyt $Ozt $Pyt $Pzt
layer straight 7 2 0.31 $Qyt $Qzt $Ryt $Rzt
layer straight 8 2 0.11 $Syt $Szt $Tyt $Tzt
layer straight 8 2 0.11 $Uyt $Uzt $Vyt $Vzt
}
set GJ 1.0e20;
uniaxialMaterial Elastic 50 $GJ
section Aggregator 2 50 T -section 1
#geomTransf PDelta 1 1 0 1
geomTransf Linear 1 1 0 1
# Define element
element nonlinearBeamColumn 1 1 2 5 2 1 -iter 1500 1.0e-12
element elasticBeamColumn 2 2 3 100000 3000000 1000000 2000000 2000000 1000000 1
# Create recorder
recorder Node -file NW1-disp(y).out -time -node 2 -dof 2 disp
recorder Node -file NW1-disp(z).out -time -node 2 -dof 3 disp
recorder Node -file NW1-control(y).out -time -node 3 -dof 2 disp
recorder Node -file NW1-control(z).out -time -node 3 -dof 3 disp
# Create the steel fiber stress strain recorders
recorder Element -file NW1section.out -time -ele 1 globalForce
recorder Element -file bar-1stress.out -time -ele 1 section 1 fiber $ayt $azt 6 stressStrain
recorder Element -file bar-2stress.out -time -ele 1 section 1 fiber $gyt $gzt 6 stressStrain
recorder Element -file bar-3stress.out -time -ele 1 section 1 fiber $pyt $pzt 6 stressStrain
recorder Element -file bar-4stress.out -time -ele 1 section 1 fiber $vyt $vzt 6 stressStrain
recorder Element -file bar-5stress.out -time -ele 1 section 1 fiber $Jyt $Jzt 6 stressStrain
recorder Element -file bar-6stress.out -time -ele 1 section 1 fiber $Pyt $Pzt 6 stressStrain
recorder Element -file bar-7stress.out -time -ele 1 section 1 fiber $yyt $yzt 8 stressStrain
recorder Element -file bar-8stress.out -time -ele 1 section 1 fiber $zyt $zzt 8 stressStrain
recorder Element -file bar-9stress.out -time -ele 1 section 1 fiber $cyt $czt 6 stressStrain
recorder Element -file bar-10stress.out -time -ele 1 section 1 fiber $nyt $nzt 6 stressStrain
# Create the concrete fiber stress strain recorders
recorder Element -file conc-Astress.out -time -ele 1 section 1 fiber -35.375 42 1 stressStrain
recorder Element -file conc-Bstress.out -time -ele 1 section 1 fiber 35.375 42 1 stressStrain
recorder Element -file conc-Cstress.out -time -ele 1 section 1 fiber 0 44.0625 3 stressStrain
recorder Element -file conc-Dstress.out -time -ele 1 section 1 fiber 0 -44.375 2 stressStrain
# Define constant axial load
set P -186.5
pattern Plain 1 "Constant" {
load 2 $P 0.0 0.0 0.0 0.0 0.0
}
# Define analysis parameters
integrator LoadControl 0 1 0 0
system SparseGeneral -piv
#system UmfPack
test NormDispIncr 1.0e-5 5000
numberer Plain
constraints Plain
algorithm KrylovNewton
#algorithm Newton
analysis Static
# Do one analysis for constant axial load
analyze 1
fix 3 0 0 1 0 0 0
# Define reference load in Y direction
pattern Plain 3 "Linear" {
# Fx Fy Fz Mx My Mz
load 3 0.0 1.0 0.0 0.0 0.0 24.0
}
# Maximum displacement
set maxU 10.00
set numIncr 500
set dU [expr $maxU/$numIncr]
# Use displacement control at node 3 for analysis
integrator DisplacementControl 3 2 $dU
# Perform the analysis
analyze $numIncr
# theta (rotation angle) 0 degrees CW
set sintheta 0
set costheta 1
# Create reference points for transformation and section generation
# -----------------------
# Confined conrete reference points
set 1y -35.375
set 1z 44.375
set 2y -16.875
set 2z 44.375
set 3y -16.875
set 3z 39.625
set 4y -35.375
set 4z 39.625
set 5y 16.875
set 5z 44.375
set 6y 35.375
set 6z 44.375
set 7y 35.375
set 7z 39.625
set 8y 16.875
set 8z 39.625
set 9y -16.875
set 9z 44.0625
set 10y 16.875
set 10z 44.0625
set 11y 16.875
set 11z 39.9375
set 12y -16.875
set 12z 39.9375
set 13y -2.0625
set 13z 39.9375
set 14y 2.0625
set 14z 39.9375
set 15y 2.0625
set 15z -22.5625
set 16y -2.0625
set 16z -22.5625
set 17y -2.375
set 17z -22.5625
set 18y 2.375
set 18z -22.5625
set 19y 2.375
set 19z -44.375
set 20y -2.375
set 20z -44.375
#Unconfined concrete reference points
set 21y -36.0
set 21z 45.0
set 22y -16.875
set 22z 45.0
set 25y 16.875
set 25z 45.0
set 26y 36.0
set 26z 45.0
set 27y 36.0
set 27z 44.375
set 28y 36.0
set 28z 39.625
set 29y 36.0
set 29z 39.0
set 30y 16.875
set 30z 39.0
set 32y 2.0625
set 32z 39.0
set 33y 3.0
set 33z 39.0
set 34y 3.0
set 34z -22.5625
set 36y 3.0
set 36z -45.0
set 37y 2.375
set 37z -45.0
set 38y -2.375
set 38z -45.0
set 39y -3.0
set 39z -45.0
set 40y -3.0
set 40z -22.5625
set 42y -2.0625
set 42z 39.0
set 43y -3.0
set 43z 39.0
set 44y -16.875
set 44z 39.0
set 46y -36.0
set 46z 39.0
set 47y -36.0
set 47z 39.625
set 48y -36.0
set 48z 44.375
# Create steeel layer reference points
# ------------------------
set ay -34.875
set az 43.875
set by -27.875
set bz 43.875
set cy -24.375
set cz 43.875
set dy -17.375
set dz 43.875
set ey -31.375
set ez 43.9375
set fy -20.875
set fz 43.9375
set gy -34.875
set gz 40.125
set hy -27.875
set hz 40.125
set iy -24.375
set iz 40.125
set jy -17.375
set jz 40.125
set ky -31.375
set kz 40.0625
set ly -20.875
set lz 40.0625
set my 17.375
set mz 43.875
set ny 24.375
set nz 43.875
set oy 27.875
set oz 43.875
set py 34.875
set pz 43.875
set qy 20.875
set qz 43.9375
set ry 31.375
set rz 43.9375
set sy 17.375
set sz 40.125
set ty 24.375
set tz 40.125
set uy 27.875
set uz 40.125
set vy 34.875
set vz 40.125
set wy 20.875
set wz 40.0625
set xy 31.375
set xz 40.0625
set yy -9.00
set yz 43.6875
set zy 9.00
set zz 43.6875
set Ay -6.25
set Az 40.3125
set By 6.25
set Bz 40.3125
set Cy -1.6875
set Cz 36.125
set Dy -1.6875
set Dz -13.875
set Ey 1.6875
set Ez 36.125
set Fy 1.6875
set Fz -13.875
set Gy -1.875
set Gz -26.375
set Hy -1.875
set Hz -33.375
set Iy -1.875
set Iz -36.875
set Jy -1.875
set Jz -43.875
set Ky -1.9375
set Kz -29.875
set Ly -1.9375
set Lz -40.375
set My 1.875
set Mz -26.375
set Ny 1.875
set Nz -33.375
set Oy 1.875
set Oz -36.875
set Py 1.875
set Pz -43.875
set Qy 1.9375
set Qz -29.875
set Ry 1.9375
set Rz -40.375
set Sy -2.0625
set Sz -22.875
set Ty 2.0625
set Tz -22.875
set Uy -2.5
set Uz 43.6875
set Vy 2.5
set Vz 43.6875
# Create rotated points through geometric transformation
set 1yt [expr $1y*$costheta+$1z*$sintheta]
set 1zt [expr -$1y*$sintheta+$1z*$costheta]
set 2yt [expr $2y*$costheta+$2z*$sintheta]
set 2zt [expr -$2y*$sintheta+$2z*$costheta]
set 3yt [expr $3y*$costheta+$3z*$sintheta]
set 3zt [expr -$3y*$sintheta+$3z*$costheta]
set 4yt [expr $4y*$costheta+$4z*$sintheta]
set 4zt [expr -$4y*$sintheta+$4z*$costheta]
set 5yt [expr $5y*$costheta+$5z*$sintheta]
set 5zt [expr -$5y*$sintheta+$5z*$costheta]
set 6yt [expr $6y*$costheta+$6z*$sintheta]
set 6zt [expr -$6y*$sintheta+$6z*$costheta]
set 7yt [expr $7y*$costheta+$7z*$sintheta]
set 7zt [expr -$7y*$sintheta+$7z*$costheta]
set 8yt [expr $8y*$costheta+$8z*$sintheta]
set 8zt [expr -$8y*$sintheta+$8z*$costheta]
set 9yt [expr $9y*$costheta+$9z*$sintheta]
set 9zt [expr -$9y*$sintheta+$9z*$costheta]
set 10yt [expr $10y*$costheta+$10z*$sintheta]
set 10zt [expr -$10y*$sintheta+$10z*$costheta]
set 11yt [expr $11y*$costheta+$11z*$sintheta]
set 11zt [expr -$11y*$sintheta+$11z*$costheta]
set 12yt [expr $12y*$costheta+$12z*$sintheta]
set 12zt [expr -$12y*$sintheta+$12z*$costheta]
set 13yt [expr $13y*$costheta+$13z*$sintheta]
set 13zt [expr -$13y*$sintheta+$13z*$costheta]
set 14yt [expr $14y*$costheta+$14z*$sintheta]
set 14zt [expr -$14y*$sintheta+$14z*$costheta]
set 15yt [expr $15y*$costheta+$15z*$sintheta]
set 15zt [expr -$15y*$sintheta+$15z*$costheta]
set 16yt [expr $16y*$costheta+$16z*$sintheta]
set 16zt [expr -$16y*$sintheta+$16z*$costheta]
set 17yt [expr $17y*$costheta+$17z*$sintheta]
set 17zt [expr -$17y*$sintheta+$17z*$costheta]
set 18yt [expr $18y*$costheta+$18z*$sintheta]
set 18zt [expr -$18y*$sintheta+$18z*$costheta]
set 19yt [expr $19y*$costheta+$19z*$sintheta]
set 19zt [expr -$19y*$sintheta+$19z*$costheta]
set 20yt [expr $20y*$costheta+$20z*$sintheta]
set 20zt [expr -$20y*$sintheta+$20z*$costheta]
set 21yt [expr $21y*$costheta+$21z*$sintheta]
set 21zt [expr -$21y*$sintheta+$21z*$costheta]
set 22yt [expr $22y*$costheta+$22z*$sintheta]
set 22zt [expr -$22y*$sintheta+$22z*$costheta]
set 25yt [expr $25y*$costheta+$25z*$sintheta]
set 25zt [expr -$25y*$sintheta+$25z*$costheta]
set 26yt [expr $26y*$costheta+$26z*$sintheta]
set 26zt [expr -$26y*$sintheta+$26z*$costheta]
set 27yt [expr $27y*$costheta+$27z*$sintheta]
set 27zt [expr -$27y*$sintheta+$27z*$costheta]
set 28yt [expr $28y*$costheta+$28z*$sintheta]
set 28zt [expr -$28y*$sintheta+$28z*$costheta]
set 29yt [expr $29y*$costheta+$29z*$sintheta]
set 29zt [expr -$29y*$sintheta+$29z*$costheta]
set 30yt [expr $30y*$costheta+$30z*$sintheta]
set 30zt [expr -$30y*$sintheta+$30z*$costheta]
set 32yt [expr $32y*$costheta+$32z*$sintheta]
set 32zt [expr -$32y*$sintheta+$32z*$costheta]
set 33yt [expr $33y*$costheta+$33z*$sintheta]
set 33zt [expr -$33y*$sintheta+$33z*$costheta]
set 34yt [expr $34y*$costheta+$34z*$sintheta]
set 34zt [expr -$34y*$sintheta+$34z*$costheta]
set 36yt [expr $36y*$costheta+$36z*$sintheta]
set 36zt [expr -$36y*$sintheta+$36z*$costheta]
set 37yt [expr $37y*$costheta+$37z*$sintheta]
set 37zt [expr -$37y*$sintheta+$37z*$costheta]
set 38yt [expr $38y*$costheta+$38z*$sintheta]
set 38zt [expr -$38y*$sintheta+$38z*$costheta]
set 39yt [expr $39y*$costheta+$39z*$sintheta]
set 39zt [expr -$39y*$sintheta+$39z*$costheta]
set 40yt [expr $40y*$costheta+$40z*$sintheta]
set 40zt [expr -$40y*$sintheta+$40z*$costheta]
set 42yt [expr $42y*$costheta+$42z*$sintheta]
set 42zt [expr -$42y*$sintheta+$42z*$costheta]
set 43yt [expr $43y*$costheta+$43z*$sintheta]
set 43zt [expr -$43y*$sintheta+$43z*$costheta]
set 44yt [expr $44y*$costheta+$44z*$sintheta]
set 44zt [expr -$44y*$sintheta+$44z*$costheta]
set 46yt [expr $46y*$costheta+$46z*$sintheta]
set 46zt [expr -$46y*$sintheta+$46z*$costheta]
set 47yt [expr $47y*$costheta+$47z*$sintheta]
set 47zt [expr -$47y*$sintheta+$47z*$costheta]
set 48yt [expr $48y*$costheta+$48z*$sintheta]
set 48zt [expr -$48y*$sintheta+$48z*$costheta]
set ayt [expr $ay*$costheta+$az*$sintheta]
set azt [expr -$ay*$sintheta+$az*$costheta]
set byt [expr $by*$costheta+$bz*$sintheta]
set bzt [expr -$by*$sintheta+$bz*$costheta]
set cyt [expr $cy*$costheta+$cz*$sintheta]
set czt [expr -$cy*$sintheta+$cz*$costheta]
set dyt [expr $dy*$costheta+$dz*$sintheta]
set dzt [expr -$dy*$sintheta+$dz*$costheta]
set eyt [expr $ey*$costheta+$ez*$sintheta]
set ezt [expr -$ey*$sintheta+$ez*$costheta]
set fyt [expr $fy*$costheta+$fz*$sintheta]
set fzt [expr -$fy*$sintheta+$fz*$costheta]
set gyt [expr $gy*$costheta+$gz*$sintheta]
set gzt [expr -$gy*$sintheta+$gz*$costheta]
set hyt [expr $hy*$costheta+$hz*$sintheta]
set hzt [expr -$hy*$sintheta+$hz*$costheta]
set iyt [expr $iy*$costheta+$iz*$sintheta]
set izt [expr -$iy*$sintheta+$iz*$costheta]
set jyt [expr $jy*$costheta+$jz*$sintheta]
set jzt [expr -$jy*$sintheta+$jz*$costheta]
set kyt [expr $ky*$costheta+$kz*$sintheta]
set kzt [expr -$ky*$sintheta+$kz*$costheta]
set lyt [expr $ly*$costheta+$lz*$sintheta]
set lzt [expr -$ly*$sintheta+$lz*$costheta]
set myt [expr $my*$costheta+$mz*$sintheta]
set mzt [expr -$my*$sintheta+$mz*$costheta]
set nyt [expr $ny*$costheta+$nz*$sintheta]
set nzt [expr -$ny*$sintheta+$nz*$costheta]
set oyt [expr $oy*$costheta+$oz*$sintheta]
set ozt [expr -$oy*$sintheta+$oz*$costheta]
set pyt [expr $py*$costheta+$pz*$sintheta]
set pzt [expr -$py*$sintheta+$pz*$costheta]
set qyt [expr $qy*$costheta+$qz*$sintheta]
set qzt [expr -$qy*$sintheta+$qz*$costheta]
set ryt [expr $ry*$costheta+$rz*$sintheta]
set rzt [expr -$ry*$sintheta+$rz*$costheta]
set syt [expr $sy*$costheta+$sz*$sintheta]
set szt [expr -$sy*$sintheta+$sz*$costheta]
set tyt [expr $ty*$costheta+$tz*$sintheta]
set tzt [expr -$ty*$sintheta+$tz*$costheta]
set uyt [expr $uy*$costheta+$uz*$sintheta]
set uzt [expr -$uy*$sintheta+$uz*$costheta]
set vyt [expr $vy*$costheta+$vz*$sintheta]
set vzt [expr -$vy*$sintheta+$vz*$costheta]
set wyt [expr $wy*$costheta+$wz*$sintheta]
set wzt [expr -$wy*$sintheta+$wz*$costheta]
set xyt [expr $xy*$costheta+$xz*$sintheta]
set xzt [expr -$xy*$sintheta+$xz*$costheta]
set yyt [expr $yy*$costheta+$yz*$sintheta]
set yzt [expr -$yy*$sintheta+$yz*$costheta]
set zyt [expr $zy*$costheta+$zz*$sintheta]
set zzt [expr -$zy*$sintheta+$zz*$costheta]
set Ayt [expr $Ay*$costheta+$Az*$sintheta]
set Azt [expr -$Ay*$sintheta+$Az*$costheta]
set Byt [expr $By*$costheta+$Bz*$sintheta]
set Bzt [expr -$By*$sintheta+$Bz*$costheta]
set Cyt [expr $Cy*$costheta+$Cz*$sintheta]
set Czt [expr -$Cy*$sintheta+$Cz*$costheta]
set Dyt [expr $Dy*$costheta+$Dz*$sintheta]
set Dzt [expr -$Dy*$sintheta+$Dz*$costheta]
set Eyt [expr $Ey*$costheta+$Ez*$sintheta]
set Ezt [expr -$Ey*$sintheta+$Ez*$costheta]
set Fyt [expr $Fy*$costheta+$Fz*$sintheta]
set Fzt [expr -$Fy*$sintheta+$Fz*$costheta]
set Gyt [expr $Gy*$costheta+$Gz*$sintheta]
set Gzt [expr -$Gy*$sintheta+$Gz*$costheta]
set Hyt [expr $Hy*$costheta+$Hz*$sintheta]
set Hzt [expr -$Hy*$sintheta+$Hz*$costheta]
set Iyt [expr $Iy*$costheta+$Iz*$sintheta]
set Izt [expr -$Iy*$sintheta+$Iz*$costheta]
set Jyt [expr $Jy*$costheta+$Jz*$sintheta]
set Jzt [expr -$Jy*$sintheta+$Jz*$costheta]
set Kyt [expr $Ky*$costheta+$Kz*$sintheta]
set Kzt [expr -$Ky*$sintheta+$Kz*$costheta]
set Lyt [expr $Ly*$costheta+$Lz*$sintheta]
set Lzt [expr -$Ly*$sintheta+$Lz*$costheta]
set Myt [expr $My*$costheta+$Mz*$sintheta]
set Mzt [expr -$My*$sintheta+$Mz*$costheta]
set Nyt [expr $Ny*$costheta+$Nz*$sintheta]
set Nzt [expr -$Ny*$sintheta+$Nz*$costheta]
set Oyt [expr $Oy*$costheta+$Oz*$sintheta]
set Ozt [expr -$Oy*$sintheta+$Oz*$costheta]
set Pyt [expr $Py*$costheta+$Pz*$sintheta]
set Pzt [expr -$Py*$sintheta+$Pz*$costheta]
set Qyt [expr $Qy*$costheta+$Qz*$sintheta]
set Qzt [expr -$Qy*$sintheta+$Qz*$costheta]
set Ryt [expr $Ry*$costheta+$Rz*$sintheta]
set Rzt [expr -$Ry*$sintheta+$Rz*$costheta]
set Syt [expr $Sy*$costheta+$Sz*$sintheta]
set Szt [expr -$Sy*$sintheta+$Sz*$costheta]
set Tyt [expr $Ty*$costheta+$Tz*$sintheta]
set Tzt [expr -$Ty*$sintheta+$Tz*$costheta]
set Uyt [expr $Uy*$costheta+$Uz*$sintheta]
set Uzt [expr -$Uy*$sintheta+$Uz*$costheta]
set Vyt [expr $Vy*$costheta+$Vz*$sintheta]
set Vzt [expr -$Vy*$sintheta+$Vz*$costheta]
# ------------------------
# Start model generation
# ------------------------
model basic -ndm 3 -ndf 6
# Create nodes
node 1 0.0 0.0 0.0
node 2 288.0 0.0 0.0
node 3 312.0 0.0 0.0
#node 3 72.0 0.0 0.0
#node 4 144.0 0.0 0.0
#node 5 216.0 0.0 0.0
#node 6 288.0 0.0 0.0
# Fix all DOF
fix 1 1 1 1 1 1 1
fix 2 0 0 0 1 0 0
fix 3 0 0 0 1 0 0
# Create uniaxial material models for the different fibers
# Concrete Material Properties
# ------------------------------
# Confined Concrete (Numbering corresponds to drawing)
uniaxialMaterial Concrete03 1 -6.551 -0.005101 -4.678 -0.021 0.55 0.5303 0.002 .00004 2.5 0.003
#uniaxialMaterial Concrete01 1 -6.551 -0.005101 -4.678 -0.021
uniaxialMaterial Concrete03 2 -7.318 -0.006637 -5.557 -0.028 0.55 0.5303 0.002 .00004 2.5 0.003
#uniaxialMaterial Concrete01 2 -7.318 -0.006637 -5.557 -0.028
uniaxialMaterial Concrete03 3 -5.790 -0.003579 -3.729 -0.013 0.55 0.5303 0.002 .00004 2.5 0.003
#uniaxialMaterial Concrete01 3 -5.790 -0.003579 -3.729 -0.013
uniaxialMaterial Concrete03 4 -6.085 -0.004170 -4.106 -0.016 0.55 0.5303 0.002 .00004 2.5 0.003
#uniaxialMaterial Concrete01 4 -6.085 -0.004170 -4.106 -0.016
# Unconfined Cover Concrete
uniaxialMaterial Concrete03 5 -5.000 -0.003 0.00 -0.008 0.1 0.5303 0.002 .00004 2.5 0.003
#uniaxialMaterial Concrete01 5 -5.00 -0.003 -0.00 -0.008
# Steel Material Properties
# ------------------------------
# Material Properties for #6 bars
uniaxialMaterial Steel02 6 70.0 29000 0.01 18.5 .925 .15 0 5 0 5
# Material Properties for #5 bars
uniaxialMaterial Steel02 7 70.0 29000 0.01 18.5 .925 .15 0 5 0 5
# Material Properties for #3 bars
uniaxialMaterial Steel02 8 70.0 29000 0.01 18.5 .925 .15 0 5 0 5
# Bond Material Properties
# -------------------------------
# Material Properties for #6 bar
uniaxialMaterial StrPen01 9 0.01897 70 0.664 103 0.3 0.8 0.0 1.0 4.35 25.0
# Material Properties for #5 bar
uniaxialMaterial StrPen01 10 0.0169 70 0.593 103 0.3 0.8 0.0 1.0 4.35 25.0
# Material Properties for #3 bar
uniaxialMaterial StrPen01 11 0.0144 70 0.503 103 0.3 0.8 0.0 1.0 4.35 25.0
# Material Properties for interface concrete fibers
uniaxialMaterial Concrete03 12 -6.551 -0.005101 -6.651 -0.04 0.1 0.5303 0.002 .00004 2.5 0.003
uniaxialMaterial Concrete03 13 -7.318 -0.006637 -7.418 -0.04 0.1 0.5303 0.002 .00004 2.5 0.003
uniaxialMaterial Concrete03 14 -5.790 -0.003579 -5.900 -0.04 0.1 0.5303 0.002 .00004 2.5 0.003
uniaxialMaterial Concrete03 15 -6.085 -0.004170 -6.185 -0.04 0.1 0.5303 0.002 .00004 2.5 0.003
uniaxialMaterial Concrete03 16 -5.000 -0.003 -5.100 -0.04 0.1 0.5303 0.002 .00004 2.5 0.003
# Proof of Theory, softening the fange boundary elements
uniaxialMaterial ElasticPPGap 17 12000 63 0.0002
#uniaxialMaterial ElasticPPGap 17 5000 63 0.0
uniaxialMaterial Series 18 6 17
uniaxialMaterial Series 19 7 17
set GJ 1.0e20
# Define cross-section for analysis
# Fibers are approx. 0.25" by 0.25"
section Fiber 1 {
# Create Core Concrete Fibers
patch quad 1 65 17 $4yt $4zt $3yt $3zt $2yt $2zt $1yt $1zt
patch quad 1 65 17 $8yt $8zt $7yt $7zt $6yt $6zt $5yt $5zt
patch quad 3 150 13 $12yt $12zt $11yt $11zt $10yt $10zt $9yt $9zt
patch quad 4 13 273 $16yt $16zt $15yt $15zt $14yt $14zt $13yt $13zt
patch quad 2 17 65 $20yt $20zt $19yt $19zt $18yt $18zt $17yt $17zt
# Create Cover Concrete Fibers
patch quad 5 69 4 $48yt $48zt $2yt $2zt $22yt $22zt $21yt $21zt
patch quad 5 150 6 $9yt $9zt $10yt $10zt $25yt $25zt $22yt $22zt
patch quad 5 69 4 $5yt $5zt $27yt $27zt $26yt $26zt $25yt $25zt
patch quad 5 4 17 $7yt $7zt $28yt $28zt $27yt $27zt $6yt $6zt
patch quad 5 69 4 $30yt $30zt $29yt $29zt $28yt $28zt $8yt $8zt
patch quad 5 69 6 $32yt $32zt $30yt $30zt $11yt $11zt $14yt $14zt
patch quad 5 6 267 $15yt $15zt $34yt $34zt $33yt $33zt $32yt $32zt
patch quad 5 4 69 $37yt $37zt $36yt $36zt $34yt $34zt $18yt $18zt
patch quad 5 17 4 $38yt $38zt $37yt $37zt $19yt $19zt $20yt $20zt
patch quad 5 4 69 $39yt $39zt $38yt $38zt $17yt $17zt $40yt $40zt
patch quad 5 6 267 $40yt $40zt $16yt $16zt $42yt $42zt $43yt $43zt
patch quad 5 69 6 $44yt $44zt $42yt $42zt $13yt $13zt $12yt $12zt
patch quad 5 69 4 $46yt $46zt $44yt $44zt $3yt $3zt $47yt $47zt
patch quad 5 4 17 $47yt $47zt $4yt $4zt $1yt $1zt $48yt $48zt
# Create Reinforcing Steel Fibers
layer straight 6 2 0.44 $ayt $azt $byt $bzt
layer straight 6 2 0.44 $cyt $czt $dyt $dzt
layer straight 7 2 0.31 $eyt $ezt $fyt $fzt
layer straight 6 2 0.44 $gyt $gzt $hyt $hzt
layer straight 6 2 0.44 $iyt $izt $jyt $jzt
layer straight 7 2 0.31 $kyt $kzt $lyt $lzt
layer straight 6 2 0.44 $myt $mzt $nyt $nzt
layer straight 6 2 0.44 $oyt $ozt $pyt $pzt
layer straight 7 2 0.31 $qyt $qzt $ryt $rzt
layer straight 6 2 0.44 $syt $szt $tyt $tzt
layer straight 6 2 0.44 $uyt $uzt $vyt $vzt
layer straight 7 2 0.31 $wyt $wzt $xyt $xzt
layer straight 8 2 0.11 $yyt $yzt $zyt $zzt
layer straight 8 2 0.11 $Ayt $Azt $Byt $Bzt
layer straight 8 5 0.11 $Cyt $Czt $Dyt $Dzt
layer straight 8 5 0.11 $Eyt $Ezt $Fyt $Fzt
layer straight 6 2 0.44 $Gyt $Gzt $Hyt $Hzt
layer straight 6 2 0.44 $Iyt $Izt $Jyt $Jzt
layer straight 7 2 0.31 $Kyt $Kzt $Lyt $Lzt
layer straight 6 2 0.44 $Myt $Mzt $Nyt $Nzt
layer straight 6 2 0.44 $Oyt $Ozt $Pyt $Pzt
layer straight 7 2 0.31 $Qyt $Qzt $Ryt $Rzt
layer straight 8 2 0.11 $Syt $Szt $Tyt $Tzt
layer straight 8 2 0.11 $Uyt $Uzt $Vyt $Vzt
}
set GJ 1.0e20;
uniaxialMaterial Elastic 50 $GJ
section Aggregator 2 50 T -section 1
#geomTransf PDelta 1 1 0 1
geomTransf Linear 1 1 0 1
# Define element
element nonlinearBeamColumn 1 1 2 5 2 1 -iter 1500 1.0e-12
element elasticBeamColumn 2 2 3 100000 3000000 1000000 2000000 2000000 1000000 1
# Create recorder
recorder Node -file NW1-disp(y).out -time -node 2 -dof 2 disp
recorder Node -file NW1-disp(z).out -time -node 2 -dof 3 disp
recorder Node -file NW1-control(y).out -time -node 3 -dof 2 disp
recorder Node -file NW1-control(z).out -time -node 3 -dof 3 disp
# Create the steel fiber stress strain recorders
recorder Element -file NW1section.out -time -ele 1 globalForce
recorder Element -file bar-1stress.out -time -ele 1 section 1 fiber $ayt $azt 6 stressStrain
recorder Element -file bar-2stress.out -time -ele 1 section 1 fiber $gyt $gzt 6 stressStrain
recorder Element -file bar-3stress.out -time -ele 1 section 1 fiber $pyt $pzt 6 stressStrain
recorder Element -file bar-4stress.out -time -ele 1 section 1 fiber $vyt $vzt 6 stressStrain
recorder Element -file bar-5stress.out -time -ele 1 section 1 fiber $Jyt $Jzt 6 stressStrain
recorder Element -file bar-6stress.out -time -ele 1 section 1 fiber $Pyt $Pzt 6 stressStrain
recorder Element -file bar-7stress.out -time -ele 1 section 1 fiber $yyt $yzt 8 stressStrain
recorder Element -file bar-8stress.out -time -ele 1 section 1 fiber $zyt $zzt 8 stressStrain
recorder Element -file bar-9stress.out -time -ele 1 section 1 fiber $cyt $czt 6 stressStrain
recorder Element -file bar-10stress.out -time -ele 1 section 1 fiber $nyt $nzt 6 stressStrain
# Create the concrete fiber stress strain recorders
recorder Element -file conc-Astress.out -time -ele 1 section 1 fiber -35.375 42 1 stressStrain
recorder Element -file conc-Bstress.out -time -ele 1 section 1 fiber 35.375 42 1 stressStrain
recorder Element -file conc-Cstress.out -time -ele 1 section 1 fiber 0 44.0625 3 stressStrain
recorder Element -file conc-Dstress.out -time -ele 1 section 1 fiber 0 -44.375 2 stressStrain
# Define constant axial load
set P -186.5
pattern Plain 1 "Constant" {
load 2 $P 0.0 0.0 0.0 0.0 0.0
}
# Define analysis parameters
integrator LoadControl 0 1 0 0
system SparseGeneral -piv
#system UmfPack
test NormDispIncr 1.0e-5 5000
numberer Plain
constraints Plain
algorithm KrylovNewton
#algorithm Newton
analysis Static
# Do one analysis for constant axial load
analyze 1
fix 3 0 0 1 0 0 0
# Define reference load in Y direction
pattern Plain 3 "Linear" {
# Fx Fy Fz Mx My Mz
load 3 0.0 1.0 0.0 0.0 0.0 24.0
}
# Maximum displacement
set maxU 10.00
set numIncr 500
set dU [expr $maxU/$numIncr]
# Use displacement control at node 3 for analysis
integrator DisplacementControl 3 2 $dU
# Perform the analysis
analyze $numIncr
i just got the following warning, what is it?
OpenSees -- Open System For Earthquake Engineering Simulation
Pacific Earthquake Engineering Research Center -- Version 1.6.2.e
(c) Copyright 1999 The Regents of the University of California
All Rights Reserved
OpenSees > source test.tcl
WARNING could not create uniaxialMaterial StrPen01
OpenSees >
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
SimpleQuestion
I remember in thin wall theory, where the warpping effects are significant, only when axial load is at the shear center (not gravity centroid), will no lateral movement be induced. Actually any sections other than think wall section will have such attribution, but in a much small scale. So I guess if you section is complicated, this might be the reason. 