recorder elements for shear wall
Moderators: silvia, selimgunay, Moderators
-
- Posts: 17
- Joined: Thu Dec 29, 2011 1:10 am
recorder elements for shear wall
hi
i have wrote a script for apply IDA method on Simple Steel Framing with Reinforced Concrete Shear Wall (5,8 and 12 story). by using concrete01 steel02 material, disp beamcolumnInt for shear wall and nonlinear beamcolumn for beam and column elements, i want get stress or strain response, for beam and column elements it gives output that first column is time second is stress and third is strain by this command:
recorder Element -file $RecordDir/Accel/ss/$acc/1.out -time -ele 1 section 1 fiber 0 0 1 stressStrain
i use this command for shearwall elements but it gives only time column! i have used stresses instead of stressStrain but still give only time column!
is the reason of this problem the using of DBE elements?
i wait your suggests. . .
i have wrote a script for apply IDA method on Simple Steel Framing with Reinforced Concrete Shear Wall (5,8 and 12 story). by using concrete01 steel02 material, disp beamcolumnInt for shear wall and nonlinear beamcolumn for beam and column elements, i want get stress or strain response, for beam and column elements it gives output that first column is time second is stress and third is strain by this command:
recorder Element -file $RecordDir/Accel/ss/$acc/1.out -time -ele 1 section 1 fiber 0 0 1 stressStrain
i use this command for shearwall elements but it gives only time column! i have used stresses instead of stressStrain but still give only time column!
is the reason of this problem the using of DBE elements?
i wait your suggests. . .
-
- Posts: 35
- Joined: Thu Apr 12, 2012 10:09 am
- Location: National Technical University of Athens
Re: recorder elements for shear wall
How did you use disp beamcolumnInt for shear wall ( as a big column??) and how did you connect it with the frame?
-
- Posts: 17
- Joined: Thu Dec 29, 2011 1:10 am
Re: recorder elements problem
gerber wrote:
> How did you use disp beamcolumnInt for shear wall ( as a big column??) and
> how did you connect it with the frame?
hi gerber
i have used disp beamcolumnInt similar what in opensees manual said for using this command. i have not defined connection shear wall whit frame. how i can defin this? can you help me that solve my problem
it is apart of my script that consist of modeling, defining recorders and static analysis parameters:
# --------------------------------------------------------------------------------------------------
#
# SET UP ----------------------------------------------------------------------------
wipe; # clear memory of all past model definitions
model BasicBuilder -ndm 2 -ndf 3; # Define the model builder, ndm=#dimension, ndf=#dofs
set RecordDir $Record; # set up name of data directory (can remove this)
file mkdir $RecordDir;
#file mkdir $RecordDir/drifts
file mkdir $RecordDir/Accel/drift/$acc
file mkdir $RecordDir/Accel/disp/$acc
file mkdir $RecordDir/Accel/envelopdrif/$acc
file mkdir $RecordDir/Accel/stressStrain/$acc
file mkdir $RecordDir/Accel/stressStrain1003/$acc
# create data directory
set GMdir "C:/tcleditor/bin/GMfiles"; # ground-motion file directory
#source LibUnits.tcl; # define units
source DisplayPlane.tcl; # procedure for displaying a plane in model
source DisplayModel2D.tcl; # procedure for displaying 3D perspectives of model
# define GEOMETRY -------------------------------------------------------------
# define NODAL COORDINATES
node 1 0.0 0.0
node 2 5.0 0.0
node 3 7.5 0.0
node 4 10.0 0.0
node 5 15.0 0.0
node 6 0.0 3.0
node 7 5.0 3.0
node 8 7.5 3.0
node 9 10.0 3.0
node 10 15.0 3.0
node 11 0.0 6.5
node 12 5.0 6.5
node 13 7.5 6.5
node 14 10.0 6.5
node 15 15.0 6.5
node 16 0.0 10.0
node 17 5.0 10.0
node 18 7.5 10.0
node 19 10.0 10.0
node 20 15.0 10.0
node 21 0.0 13.5
node 22 5.0 13.5
node 23 7.5 13.5
node 24 10.0 13.5
node 25 15.0 13.5
node 26 0.0 17.0
node 27 5.0 17.0
node 28 7.5 17.0
node 29 10.0 17.0
node 30 15.0 17.0
node 31 0.0 20.5
node 32 5.0 20.5
node 33 7.5 20.5
node 34 10.0 20.5
node 35 15.0 20.5
node 36 0.0 24.0
node 37 5.0 24.0
node 38 7.5 24.0
node 39 10.0 24.0
node 40 15.0 24.0
node 41 0.0 27.5
node 42 5.0 27.5
node 43 7.5 27.5
node 44 10.0 27.5
node 45 15.0 27.5
node 46 0.0 31.0
node 47 5.0 31.0
node 48 7.5 31.0
node 49 10.0 31.0
node 50 15.0 31.0
node 51 0.0 34.5
node 52 5.0 34.5
node 53 7.5 34.5
node 54 10.0 34.5
node 55 15.0 34.5
node 56 0.0 38.0
node 57 5.0 38.0
node 58 7.5 38.0
node 59 10.0 38.0
node 60 15.0 38.0
node 61 0.0 41.5
node 62 5.0 41.5
node 63 7.5 41.5
node 64 10.0 41.5
node 65 15.0 41.5
# Nodes for shear wall in story 1:
node 101 7.5 0.3
node 102 7.5 0.6
node 103 7.5 0.9
node 104 7.5 1.2
node 105 7.5 1.5
node 106 7.5 1.8
node 107 7.5 2.1
node 108 7.5 2.4
node 109 7.5 2.7
# Nodes for shear wall in story 2:
node 201 7.5 3.3
node 202 7.5 3.6
node 203 7.5 3.9
node 204 7.5 4.2
node 205 7.5 4.5
node 206 7.5 4.8
node 207 7.5 5.1
node 208 7.5 5.4
node 209 7.5 5.7
node 210 7.5 6.0
node 211 7.5 6.3
# Nodes for shear wall in story 3:
node 301 7.5 6.8
node 302 7.5 7.1
node 303 7.5 7.4
node 304 7.5 7.7
node 305 7.5 8.0
node 306 7.5 8.3
node 307 7.5 8.6
node 308 7.5 8.9
node 309 7.5 9.2
node 310 7.5 9.5
node 311 7.5 9.8
# Nodes for shear wall in story 4:
node 401 7.5 10.3
node 402 7.5 10.6
node 403 7.5 10.9
node 404 7.5 11.2
node 405 7.5 11.5
node 406 7.5 11.8
node 407 7.5 12.1
node 408 7.5 12.4
node 409 7.5 12.7
node 410 7.5 13.0
node 411 7.5 13.3
# Nodes for shear wall in story 5:
node 501 7.5 13.8
node 502 7.5 14.1
node 503 7.5 14.4
node 504 7.5 14.7
node 505 7.5 15.0
node 506 7.5 15.3
node 507 7.5 15.6
node 508 7.5 15.9
node 509 7.5 16.2
node 510 7.5 16.5
node 511 7.5 16.8
# Nodes for shear wall in story 6:
node 601 7.5 17.3
node 602 7.5 17.6
node 603 7.5 17.9
node 604 7.5 18.2
node 605 7.5 18.5
node 606 7.5 18.8
node 607 7.5 19.1
node 608 7.5 19.4
node 609 7.5 19.7
node 610 7.5 20.0
node 611 7.5 20.3
# Nodes for shear wall in story 7:
node 701 7.5 20.8
node 702 7.5 21.1
node 703 7.5 21.4
node 704 7.5 21.7
node 705 7.5 22.0
node 706 7.5 22.3
node 707 7.5 22.6
node 708 7.5 22.9
node 709 7.5 23.2
node 710 7.5 23.5
node 711 7.5 23.8
# Nodes for shear wall in story 8:
node 801 7.5 24.3
node 802 7.5 24.6
node 803 7.5 24.9
node 804 7.5 25.2
node 805 7.5 25.5
node 806 7.5 25.8
node 807 7.5 26.1
node 808 7.5 26.4
node 809 7.5 26.7
node 810 7.5 27.0
node 811 7.5 27.3
# Nodes for shear wall in story 9:
node 901 7.5 27.8
node 902 7.5 28.1
node 903 7.5 28.4
node 904 7.5 28.7
node 905 7.5 29.0
node 906 7.5 29.3
node 907 7.5 29.6
node 908 7.5 29.9
node 909 7.5 30.2
node 910 7.5 30.5
node 911 7.5 30.8
# Nodes for shear wall in story 10:
node 1001 7.5 31.3
node 1002 7.5 31.6
node 1003 7.5 31.9
node 1004 7.5 32.2
node 1005 7.5 32.5
node 1006 7.5 32.8
node 1007 7.5 33.1
node 1008 7.5 33.4
node 1009 7.5 33.7
node 1010 7.5 34.0
node 1011 7.5 34.3
# Nodes for shear wall in story 11:
node 1101 7.5 34.8
node 1102 7.5 35.1
node 1103 7.5 35.4
node 1104 7.5 35.7
node 1105 7.5 36.0
node 1106 7.5 36.3
node 1107 7.5 36.6
node 1108 7.5 36.9
node 1109 7.5 37.2
node 1110 7.5 37.5
node 1111 7.5 37.8
# Nodes for shear wall in story 12:
node 1201 7.5 38.3
node 1202 7.5 38.6
node 1203 7.5 38.9
node 1204 7.5 39.2
node 1205 7.5 39.5
node 1206 7.5 39.8
node 1207 7.5 40.1
node 1208 7.5 40.4
node 1209 7.5 40.7
node 1210 7.5 41.0
node 1211 7.5 41.3
fix 1 1 1 1
fix 2 1 1 1
fix 3 1 1 1
fix 4 1 1 1
fix 5 1 1 1
set NStory 12
# determine support nodes where ground motions are input, for multiple-support excitation
set iSupportNode "1 2 3 4 5"
# calculated MODEL PARAMETERS, particular to this model
# Set up parameters that are particular to the model for displacement control
set IDctrlNode 63.; # node where displacement is read for displacement control
set IDctrlDOF 1; # degree of freedom of displacement read for displacement control
#set LBuilding [expr 36*$in] # total building height
#material for frame Elements :
# Define Steel property variables:
set Fy 240e5
set E 210e8
set b 0.03; #strain-hardening ratio (ratio between post-yield tangent and initial elastic tangent)
set cR0 20; #$R0=between 10 and 20
set cR1 0.925
set cR2 0.15
set a1 0.0005
set a2 0.01
set a3 0.0005
set a4 0.01
#
# ----- Columns & Beams Material -------- #
#
uniaxialMaterial Steel02 1 $Fy $E $b $cR0 $cR1 $cR2 $a1 $a2 $a3 $a4
#
#material for shear wall reinforcement:
#
uniaxialMaterial Steel02 1003 400e5 $E $b $cR0 $cR1 $cR2 $a1 $a2 $a3 $a4
#
# --------------------------------------- #
# Define Concrete property variables:
#
set fpc -210e4
set epsc0 -0.003
set fpcu -180e4
set epsU -0.004
# --------- ShearWall Material --------- #
#
uniaxialMaterial Concrete01 2 $fpc $epsc0 $fpcu $epsU
# Define Column & Beams Sections :
# --------------------------------------- #
#
#
# I Section
proc Wsection { secID matID d tw bf tf nfdw nftw nfbf nftf} {
set dw [expr ($d-2*$tf)*1e-2]
set y1 [expr (-$d*0.5)*1e-2]
set y2 [expr (-$dw*0.5)*1e-2]
set y3 [expr ($dw*0.5)*1e-2]
set y4 [expr ($d*0.5)*1e-2]
set z1 [expr (-$bf*0.5)*1e-2]
set z2 [expr (-$tw*0.5)*1e-2]
set z3 [expr ($tw*0.5)*1e-2]
set z4 [expr ($bf*0.5)*1e-2]
#
section Fiber $secID {
patch quad $matID $nftf $nfbf $y1 $z1 $y2 $z1 $y2 $z4 $y1 $z4
patch quad $matID $nfdw $nftw $y2 $z2 $y3 $z2 $y3 $z3 $y2 $z3
patch quad $matID $nftf $nfbf $y3 $z1 $y4 $z1 $y4 $z4 $y3 $z4
}
}
#
# Beam section:
# tag mat H tw bf tf nfdw nftw nfbf nftf
Wsection 1 1 25 0.10 20.0 1.00 6 2 8 4; #PG1
#
# --------------------------------------- #
proc Isection { secID1 matID d tw bf tf nfdw nftw nfbf nftf} {
set dw [expr ($d-2*$tf)]
set z1 [expr (-$d/2)*1e-2]
set z2 [expr (-$dw/2)*1e-2]
set z3 [expr ($dw/2)*1e-2]
set z4 [expr ($d/2)*1e-2]
set y1 [expr (-$bf/2)*1e-2]
set y2 [expr (-$tw/2)*1e-2]
set y3 [expr ($tw/2)*1e-2]
set y4 [expr ($bf/2)*1e-2]
#
section fiberSec $secID1 {
patch quadr $matID $nfbf $nftf $y1 $z3 $y4 $z3 $y4 $z4 $y1 $z4
patch quadr $matID $nftw $nfdw $y2 $z2 $y3 $z2 $y3 $z3 $y2 $z3
patch quadr $matID $nfbf $nftf $y1 $z1 $y4 $z1 $y4 $z2 $y1 $z2
}
}
#
# Column section:
# tag mat H tw bf tf nfdw nftw nfbf nftf
Isection 2 1 12 0.65 12.0 1.10 6 2 8 4; #HE120B
#
Isection 3 1 14 0.7 14.0 1.20 6 2 8 4; #HE140B
#
Isection 4 1 16 0.8 16.0 1.30 6 2 8 4; #HE160B
#
Isection 5 1 18 0.85 18.0 1.40 6 2 8 4; #HE180B
#
Isection 6 1 20 0.9 20.0 1.50 6 2 8 4; #HE200B
#
Isection 7 1 22 0.95 22.0 1.60 6 2 8 4; #HE220B
#
Isection 8 1 24 1.0 24.0 1.70 6 2 8 4; #HE240B
#
Isection 9 1 28 1.05 28.0 1.80 6 2 8 4; #HE160B
#
Isection 10 1 32 1.15 30.0 2.05 6 2 8 4; #HE180B
#
Isection 11 1 34 1.2 30.0 2.15 6 2 8 4; #HE200B
#
Isection 12 1 36 1.25 30.0 2.25 6 2 8 4; #HE200B
#
#shear wall section
section FiberInt 15 -NStrip 2 0.5 10 0.5 2 0.5 {
#vertical fibers
fiber -2.5 0 0.015 2; fiber -2.5 0 0.00165 1003;
fiber -2.2 0 0.02 2; fiber -2.2 0 0.00165 1003;
fiber -1.8 0 0.02 2; fiber -1.8 0 0.00165 1003;
fiber -1.4 0 0.02 2; fiber -1.4 0 0.0022 1003;
fiber -1.0 0 0.02 2; fiber -1.0 0 0.0022 1003;
fiber -0.6 0 0.02 2; fiber -0.6 0 0.0022 1003;
fiber -0.2 0 0.01 2; fiber -0.2 0 0.0022 1003;
fiber 0.2 0 0.01 2; fiber 0.2 0 0.0022 1003;
fiber 0.6 0 0.02 2; fiber 0.6 0 0.0022 1003;
fiber 1.0 0 0.02 2; fiber 1.0 0 0.0022 1003;
fiber 1.4 0 0.02 2; fiber 1.4 0 0.0022 1003;
fiber 1.8 0 0.015 2; fiber 1.8 0 0.00165 1003;
fiber 2.2 0 0.02 2; fiber 2.2 0 0.00165 1003;
fiber 2.5 0 0.015 2; fiber 2.5 0 0.00165 1003;
#horiz. reinf.
Hfiber 0 0 0.000915 1003;
}
#element dispBeamColumnInt 111 2 111 1 5 3 0.1
geomTransf Linear 1; # beam
geomTransf Linear 2; # column
geomTransf LinearInt 3; # shear wall
# column elements
element nonlinearBeamColumn 1 1 6 5 8 1
element nonlinearBeamColumn 2 2 7 5 12 1
element nonlinearBeamColumn 3 4 9 5 12 1
element nonlinearBeamColumn 4 5 10 5 8 1
element nonlinearBeamColumn 5 6 11 5 8 1
element nonlinearBeamColumn 6 7 12 5 12 1
element nonlinearBeamColumn 7 9 14 5 12 1
element nonlinearBeamColumn 8 10 15 5 8 1
element nonlinearBeamColumn 9 11 16 5 7 1
element nonlinearBeamColumn 10 12 17 5 11 1
element nonlinearBeamColumn 11 14 19 5 11 1
element nonlinearBeamColumn 12 15 20 5 7 1
element nonlinearBeamColumn 13 16 21 5 7 1
element nonlinearBeamColumn 14 17 22 5 11 1
element nonlinearBeamColumn 15 19 24 5 11 1
element nonlinearBeamColumn 16 20 25 5 7 1
element nonlinearBeamColumn 17 21 26 5 6 1
element nonlinearBeamColumn 18 22 27 5 10 1
element nonlinearBeamColumn 19 24 29 5 10 1
element nonlinearBeamColumn 20 25 30 5 6 1
element nonlinearBeamColumn 21 26 31 5 6 1
element nonlinearBeamColumn 22 27 32 5 10 1
element nonlinearBeamColumn 23 29 34 5 10 1
element nonlinearBeamColumn 24 30 35 5 6 1
element nonlinearBeamColumn 25 31 36 5 5 1
element nonlinearBeamColumn 26 32 37 5 9 1
element nonlinearBeamColumn 27 34 39 5 9 1
element nonlinearBeamColumn 28 35 40 5 5 1
element nonlinearBeamColumn 29 36 41 5 5 1
element nonlinearBeamColumn 30 37 42 5 9 1
element nonlinearBeamColumn 31 39 44 5 9 1
element nonlinearBeamColumn 32 40 45 5 5 1
element nonlinearBeamColumn 33 41 46 5 3 1
element nonlinearBeamColumn 34 42 47 5 6 1
element nonlinearBeamColumn 35 44 49 5 6 1
element nonlinearBeamColumn 36 45 50 5 3 1
element nonlinearBeamColumn 37 46 51 5 3 1
element nonlinearBeamColumn 38 47 52 5 6 1
element nonlinearBeamColumn 39 49 54 5 6 1
element nonlinearBeamColumn 40 50 55 5 3 1
element nonlinearBeamColumn 41 51 56 5 2 1
element nonlinearBeamColumn 42 52 57 5 4 1
element nonlinearBeamColumn 43 54 59 5 4 1
element nonlinearBeamColumn 44 55 60 5 2 1
element nonlinearBeamColumn 45 56 61 5 2 1
element nonlinearBeamColumn 46 57 62 5 4 1
element nonlinearBeamColumn 47 59 64 5 4 1
element nonlinearBeamColumn 48 60 65 5 2 1
# beam elements
element nonlinearBeamColumn 100 6 7 5 1 1
element nonlinearBeamColumn 101 7 8 5 1 1
element nonlinearBeamColumn 102 8 9 5 1 1
element nonlinearBeamColumn 103 9 10 5 1 1
element nonlinearBeamColumn 104 11 12 5 1 1
element nonlinearBeamColumn 105 12 13 5 1 1
element nonlinearBeamColumn 106 13 14 5 1 1
element nonlinearBeamColumn 107 14 15 5 1 1
element nonlinearBeamColumn 108 16 17 5 1 1
element nonlinearBeamColumn 109 17 18 5 1 1
element nonlinearBeamColumn 110 18 19 5 1 1
element nonlinearBeamColumn 111 19 20 5 1 1
element nonlinearBeamColumn 112 21 22 5 1 1
element nonlinearBeamColumn 113 22 23 5 1 1
element nonlinearBeamColumn 114 23 24 5 1 1
element nonlinearBeamColumn 115 24 25 5 1 1
element nonlinearBeamColumn 116 26 27 5 1 1
element nonlinearBeamColumn 117 27 28 5 1 1
element nonlinearBeamColumn 118 28 29 5 1 1
element nonlinearBeamColumn 119 29 30 5 1 1
element nonlinearBeamColumn 120 31 32 5 1 1
element nonlinearBeamColumn 121 32 33 5 1 1
element nonlinearBeamColumn 122 33 34 5 1 1
element nonlinearBeamColumn 123 34 35 5 1 1
element nonlinearBeamColumn 124 36 37 5 1 1
element nonlinearBeamColumn 125 37 38 5 1 1
element nonlinearBeamColumn 126 38 39 5 1 1
element nonlinearBeamColumn 127 39 40 5 1 1
element nonlinearBeamColumn 128 41 42 5 1 1
element nonlinearBeamColumn 129 42 43 5 1 1
element nonlinearBeamColumn 130 43 44 5 1 1
element nonlinearBeamColumn 131 44 45 5 1 1
element nonlinearBeamColumn 132 46 47 5 1 1
element nonlinearBeamColumn 133 47 48 5 1 1
element nonlinearBeamColumn 134 48 49 5 1 1
element nonlinearBeamColumn 135 49 50 5 1 1
element nonlinearBeamColumn 136 51 52 5 1 1
element nonlinearBeamColumn 137 52 53 5 1 1
element nonlinearBeamColumn 138 53 54 5 1 1
element nonlinearBeamColumn 139 54 55 5 1 1
element nonlinearBeamColumn 140 56 57 5 1 1
element nonlinearBeamColumn 141 57 58 5 1 1
element nonlinearBeamColumn 142 58 59 5 1 1
element nonlinearBeamColumn 143 59 60 5 1 1
element nonlinearBeamColumn 144 61 62 5 1 1
element nonlinearBeamColumn 145 62 63 5 1 1
element nonlinearBeamColumn 146 63 64 5 1 1
element nonlinearBeamColumn 147 64 65 5 1 1
# shear wall elements in story 1
element dispBeamColumnInt 1001 3 101 1 15 3 0.4
element dispBeamColumnInt 1002 101 102 1 15 3 0.4
element dispBeamColumnInt 1003 102 103 1 15 3 0.4
element dispBeamColumnInt 1004 103 104 1 15 3 0.4
element dispBeamColumnInt 1005 104 105 1 15 3 0.4
element dispBeamColumnInt 1006 105 106 1 15 3 0.4
element dispBeamColumnInt 1007 106 107 1 15 3 0.4
element dispBeamColumnInt 1008 107 108 1 15 3 0.4
element dispBeamColumnInt 1009 108 109 1 15 3 0.4
element dispBeamColumnInt 1010 109 8 1 15 3 0.4
# shear wall elements in story 2
element dispBeamColumnInt 2001 8 201 1 15 3 0.4
element dispBeamColumnInt 2002 201 202 1 15 3 0.4
element dispBeamColumnInt 2003 202 203 1 15 3 0.4
element dispBeamColumnInt 2004 203 204 1 15 3 0.4
element dispBeamColumnInt 2005 204 205 1 15 3 0.4
element dispBeamColumnInt 2006 205 206 1 15 3 0.4
element dispBeamColumnInt 2007 206 207 1 15 3 0.4
element dispBeamColumnInt 2008 207 208 1 15 3 0.4
element dispBeamColumnInt 2009 208 209 1 15 3 0.4
element dispBeamColumnInt 2010 209 210 1 15 3 0.4
element dispBeamColumnInt 2011 210 211 1 15 3 0.4
element dispBeamColumnInt 2012 211 13 1 15 3 0.4
# shear wall elements in story 3
element dispBeamColumnInt 3001 13 301 1 15 3 0.4
element dispBeamColumnInt 3002 301 302 1 15 3 0.4
element dispBeamColumnInt 3003 302 303 1 15 3 0.4
element dispBeamColumnInt 3004 303 304 1 15 3 0.4
element dispBeamColumnInt 3005 304 305 1 15 3 0.4
element dispBeamColumnInt 3006 305 306 1 15 3 0.4
element dispBeamColumnInt 3007 306 307 1 15 3 0.4
element dispBeamColumnInt 3008 307 308 1 15 3 0.4
element dispBeamColumnInt 3009 308 309 1 15 3 0.4
element dispBeamColumnInt 3010 309 310 1 15 3 0.4
element dispBeamColumnInt 3011 310 311 1 15 3 0.4
element dispBeamColumnInt 3012 311 18 1 15 3 0.4
# shear wall elements in story 4
element dispBeamColumnInt 4001 18 401 1 15 3 0.4
element dispBeamColumnInt 4002 401 402 1 15 3 0.4
element dispBeamColumnInt 4003 402 403 1 15 3 0.4
element dispBeamColumnInt 4004 403 404 1 15 3 0.4
element dispBeamColumnInt 4005 404 405 1 15 3 0.4
element dispBeamColumnInt 4006 405 406 1 15 3 0.4
element dispBeamColumnInt 4007 406 407 1 15 3 0.4
element dispBeamColumnInt 4008 407 408 1 15 3 0.4
element dispBeamColumnInt 4009 408 409 1 15 3 0.4
element dispBeamColumnInt 4010 409 410 1 15 3 0.4
element dispBeamColumnInt 4011 410 411 1 15 3 0.4
element dispBeamColumnInt 4012 411 23 1 15 3 0.4
# shear wall elements in story 5
element dispBeamColumnInt 5001 23 501 1 15 3 0.4
element dispBeamColumnInt 5002 501 502 1 15 3 0.4
element dispBeamColumnInt 5003 502 503 1 15 3 0.4
element dispBeamColumnInt 5004 503 504 1 15 3 0.4
element dispBeamColumnInt 5005 504 505 1 15 3 0.4
element dispBeamColumnInt 5006 505 506 1 15 3 0.4
element dispBeamColumnInt 5007 506 507 1 15 3 0.4
element dispBeamColumnInt 5008 507 508 1 15 3 0.4
element dispBeamColumnInt 5009 508 509 1 15 3 0.4
element dispBeamColumnInt 5010 509 510 1 15 3 0.4
element dispBeamColumnInt 5011 510 511 1 15 3 0.4
element dispBeamColumnInt 5012 511 28 1 15 3 0.4
# shear wall elements in story 6
element dispBeamColumnInt 6001 28 601 1 15 3 0.4
element dispBeamColumnInt 6002 601 602 1 15 3 0.4
element dispBeamColumnInt 6003 602 603 1 15 3 0.4
element dispBeamColumnInt 6004 603 604 1 15 3 0.4
element dispBeamColumnInt 6005 604 605 1 15 3 0.4
element dispBeamColumnInt 6006 605 606 1 15 3 0.4
element dispBeamColumnInt 6007 606 607 1 15 3 0.4
element dispBeamColumnInt 6008 607 608 1 15 3 0.4
element dispBeamColumnInt 6009 608 609 1 15 3 0.4
element dispBeamColumnInt 6010 609 610 1 15 3 0.4
element dispBeamColumnInt 6011 610 611 1 15 3 0.4
element dispBeamColumnInt 6012 611 33 1 15 3 0.4
# shear wall elements in story 7
element dispBeamColumnInt 7001 33 701 1 15 3 0.4
element dispBeamColumnInt 7002 701 702 1 15 3 0.4
element dispBeamColumnInt 7003 702 703 1 15 3 0.4
element dispBeamColumnInt 7004 703 704 1 15 3 0.4
element dispBeamColumnInt 7005 704 705 1 15 3 0.4
element dispBeamColumnInt 7006 705 706 1 15 3 0.4
element dispBeamColumnInt 7007 706 707 1 15 3 0.4
element dispBeamColumnInt 7008 707 708 1 15 3 0.4
element dispBeamColumnInt 7009 708 709 1 15 3 0.4
element dispBeamColumnInt 7010 709 710 1 15 3 0.4
element dispBeamColumnInt 7011 710 711 1 15 3 0.4
element dispBeamColumnInt 7012 711 38 1 15 3 0.4
# shear wall elements in story 8
element dispBeamColumnInt 8001 38 801 1 15 3 0.4
element dispBeamColumnInt 8002 801 802 1 15 3 0.4
element dispBeamColumnInt 8003 802 803 1 15 3 0.4
element dispBeamColumnInt 8004 803 804 1 15 3 0.4
element dispBeamColumnInt 8005 804 805 1 15 3 0.4
element dispBeamColumnInt 8006 805 806 1 15 3 0.4
element dispBeamColumnInt 8007 806 807 1 15 3 0.4
element dispBeamColumnInt 8008 807 808 1 15 3 0.4
element dispBeamColumnInt 8009 808 809 1 15 3 0.4
element dispBeamColumnInt 8010 809 810 1 15 3 0.4
element dispBeamColumnInt 8011 810 811 1 15 3 0.4
element dispBeamColumnInt 8012 811 43 1 15 3 0.4
# shear wall elements in story 9
element dispBeamColumnInt 9001 43 901 1 15 3 0.4
element dispBeamColumnInt 9002 901 902 1 15 3 0.4
element dispBeamColumnInt 9003 902 903 1 15 3 0.4
element dispBeamColumnInt 9004 903 904 1 15 3 0.4
element dispBeamColumnInt 9005 904 905 1 15 3 0.4
element dispBeamColumnInt 9006 905 906 1 15 3 0.4
element dispBeamColumnInt 9007 906 907 1 15 3 0.4
element dispBeamColumnInt 9008 907 908 1 15 3 0.4
element dispBeamColumnInt 9009 908 909 1 15 3 0.4
element dispBeamColumnInt 9010 909 910 1 15 3 0.4
element dispBeamColumnInt 9011 910 911 1 15 3 0.4
element dispBeamColumnInt 9012 911 48 1 15 3 0.4
# shear wall elements in story 10
element dispBeamColumnInt 10001 48 1001 1 15 3 0.4
element dispBeamColumnInt 10002 1001 1002 1 15 3 0.4
element dispBeamColumnInt 10003 1002 1003 1 15 3 0.4
element dispBeamColumnInt 10004 1003 1004 1 15 3 0.4
element dispBeamColumnInt 10005 1004 1005 1 15 3 0.4
element dispBeamColumnInt 10006 1005 1006 1 15 3 0.4
element dispBeamColumnInt 10007 1006 1007 1 15 3 0.4
element dispBeamColumnInt 10008 1007 1008 1 15 3 0.4
element dispBeamColumnInt 10009 1008 1009 1 15 3 0.4
element dispBeamColumnInt 10010 1009 1010 1 15 3 0.4
element dispBeamColumnInt 10011 1010 1011 1 15 3 0.4
element dispBeamColumnInt 10012 1011 53 1 15 3 0.4
# shear wall elements in story 11
element dispBeamColumnInt 11001 53 1101 1 15 3 0.4
element dispBeamColumnInt 11002 1101 1102 1 15 3 0.4
element dispBeamColumnInt 11003 1102 1103 1 15 3 0.4
element dispBeamColumnInt 11004 1103 1104 1 15 3 0.4
element dispBeamColumnInt 11005 1104 1105 1 15 3 0.4
element dispBeamColumnInt 11006 1105 1106 1 15 3 0.4
element dispBeamColumnInt 11007 1106 1107 1 15 3 0.4
element dispBeamColumnInt 11008 1107 1108 1 15 3 0.4
element dispBeamColumnInt 11009 1108 1109 1 15 3 0.4
element dispBeamColumnInt 11010 1109 1110 1 15 3 0.4
element dispBeamColumnInt 11011 1110 1111 1 15 3 0.4
element dispBeamColumnInt 11012 1111 58 1 15 3 0.4
# shear wall elements in story 12
element dispBeamColumnInt 12001 58 1201 1 15 3 0.4
element dispBeamColumnInt 12002 1201 1202 1 15 3 0.4
element dispBeamColumnInt 12003 1202 1203 1 15 3 0.4
element dispBeamColumnInt 12004 1203 1204 1 15 3 0.4
element dispBeamColumnInt 12005 1204 1205 1 15 3 0.4
element dispBeamColumnInt 12006 1205 1206 1 15 3 0.4
element dispBeamColumnInt 12007 1206 1207 1 15 3 0.
element dispBeamColumnInt 12008 1207 1208 1 15 3 0.4
element dispBeamColumnInt 12009 1208 1209 1 15 3 0.4
element dispBeamColumnInt 12010 1209 1210 1 15 3 0.4
element dispBeamColumnInt 12011 1210 1211 1 15 3 0.4
element dispBeamColumnInt 12012 1211 63 1 15 3 0.4
# --------------------------------------------------------------------------------------------------------------------------------
# Define GRAVITY LOADS, weight and masses
set iFloorWeight "93653 93653 93653 93653 93653 93653 93653 93653 93653 93653 93653 89169"
set WeightTotal 1489480
mass 6 1592.7 0 0
mass 7 3185.5 0 0
mass 9 3185.5 0 0
mass 10 1592.7 0 0
mass 11 1592.7 0 0
mass 12 3185.5 0 0
mass 14 3185.5 0 0
mass 15 1592.7 0 0
mass 16 1592.7 0 0
mass 17 3185.5 0 0
mass 19 3185.5 0 0
mass 20 1592.7 0 0
mass 21 1592.7 0 0
mass 22 3185.5 0 0
mass 24 3185.5 0 0
mass 25 1592.7 0 0
mass 26 1592.7 0 0
mass 27 3185.5 0 0
mass 29 3185.5 0 0
mass 30 1592.7 0 0
mass 31 1592.7 0 0
mass 32 3185.5 0 0
mass 34 3185.5 0 0
mass 35 1592.7 0 0
mass 36 1592.7 0 0
mass 37 3185.5 0 0
mass 39 3185.5 0 0
mass 40 1592.7 0 0
mass 41 1592.7 0 0
mass 42 3185.5 0 0
mass 44 3185.5 0 0
mass 45 1592.7 0 0
mass 46 1592.7 0 0
mass 47 3185.5 0 0
mass 49 3185.5 0 0
mass 50 1592.7 0 0
mass 51 1592.7 0 0
mass 52 3185.5 0 0
mass 54 3185.5 0 0
mass 55 1592.7 0 0
mass 56 1592.7 0 0
mass 57 3185.5 0 0
mass 59 3185.5 0 0
mass 60 1592.7 0 0
mass 61 1592.7 0 0
mass 62 3185.5 0 0
mass 64 3185.5 0 0
mass 65 1592.7 0 0
#mass 4 $MassNode $MassNode $MassNode 0. 0. 0.; # define mass
set MassTotal [expr $WeightTotal/9.806];
# --------------------------------------------------------------------------------------------------------------------------------
# Define RECORDERS -------------------------------------------------------------
recorder Node -file $RecordDir/DFree.out -time -node 63 -dof 1 2 3 accel;
recorder Drift -file $RecordDir/Accel/drift/$acc/1.out -iNode 3 -jNode 8 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/1.out -node 8 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/1.out -iNode 3 -jNode 8 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/1.out -time -ele 1006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/1.out -time -ele 1006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/2.out -iNode 8 -jNode 13 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/2.out -node 13 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/2.out -iNode 8 -jNode 13 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/2.out -time -ele 2006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/2.out -time -ele 2006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/3.out -iNode 13 -jNode 18 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/3.out -node 18 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/3.out -iNode 13 -jNode 18 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/3.out -time -ele 3006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/3.out -time -ele 3006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/4.out -iNode 18 -jNode 23 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/4.out -node 23 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/4.out -iNode 18 -jNode 23 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/4.out -time -ele 4006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/4.out -time -ele 4006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/5.out -iNode 23 -jNode 28 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/5.out -node 28 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/5.out -iNode 23 -jNode 28 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/5.out -time -ele 5006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/5.out -time -ele 5006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/6.out -iNode 28 -jNode 33 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/6.out -node 33 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/6.out -iNode 28 -jNode 33 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/6.out -time -ele 6006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/6.out -time -ele 6006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/7.out -iNode 33 -jNode 38 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/7.out -node 38 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/7.out -iNode 33 -jNode 38 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/7.out -time -ele 7006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/7.out -time -ele 7006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/8.out -iNode 38 -jNode 43 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/8.out -node 43 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/8.out -iNode 38 -jNode 43 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/8.out -time -ele 8006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/8.out -time -ele 8006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/9.out -iNode 43 -jNode 48 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/9.out -node 48 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/9.out -iNode 43 -jNode 48 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/9.out -time -ele 9006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/9.out -time -ele 9006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/10.out -iNode 48 -jNode 53 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/10.out -node 53 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/10.out -iNode 48 -jNode 53 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/10.out -time -ele 10006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/10.out -time -ele 10006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/11.out -iNode 53 -jNode 58 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/11.out -node 58 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/11.out -iNode 53 -jNode 58 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/11.out -time -ele 11006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/11.out -time -ele 11006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/12.out -iNode 58 -jNode 63 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/12.out -node 63 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/12.out -iNode 58 -jNode 63 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/12.out -time -ele 12006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/12.out -time -ele 12006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/globaldrift.out -iNode 63 -jNode 3 -dof 1 -perpDirn 2
# Define DISPLAY -------------------------------------------------------------
#DisplayModel2D DeformedShape ; # options: DeformedShape NodeNumbers ModeShape
# GRAVITY -------------------------------------------------------------
# define GRAVITY load applied to beams and columns -- eleLoad applies loads in local coordinate axis
#pattern Plain 101 Linear {
#
pattern Plain 1 "Linear" {
eleLoad -ele 100 -type -beamUniform -6243.5
eleLoad -ele 101 -type -beamUniform -6243.5
eleLoad -ele 102 -type -beamUniform -6243.5
eleLoad -ele 103 -type -beamUniform -6243.5
eleLoad -ele 104 -type -beamUniform -6243.5
eleLoad -ele 105 -type -beamUniform -6243.5
eleLoad -ele 106 -type -beamUniform -6243.5
eleLoad -ele 107 -type -beamUniform -6243.5
eleLoad -ele 108 -type -beamUniform -6243.5
eleLoad -ele 109 -type -beamUniform -6243.5
eleLoad -ele 110 -type -beamUniform -6243.5
eleLoad -ele 111 -type -beamUniform -6243.5
eleLoad -ele 112 -type -beamUniform -6243.5
eleLoad -ele 113 -type -beamUniform -6243.5
eleLoad -ele 114 -type -beamUniform -6243.5
eleLoad -ele 115 -type -beamUniform -6243.5
eleLoad -ele 116 -type -beamUniform -6243.5
eleLoad -ele 117 -type -beamUniform -6243.5
eleLoad -ele 118 -type -beamUniform -6243.5
eleLoad -ele 119 -type -beamUniform -6243.5
eleLoad -ele 120 -type -beamUniform -6243.5
eleLoad -ele 121 -type -beamUniform -6243.5
eleLoad -ele 122 -type -beamUniform -6243.5
eleLoad -ele 123 -type -beamUniform -6243.5
eleLoad -ele 124 -type -beamUniform -6243.5
eleLoad -ele 125 -type -beamUniform -6243.5
eleLoad -ele 126 -type -beamUniform -6243.5
eleLoad -ele 127 -type -beamUniform -6243.5
eleLoad -ele 128 -type -beamUniform -6243.5
eleLoad -ele 129 -type -beamUniform -6243.5
eleLoad -ele 130 -type -beamUniform -6243.5
eleLoad -ele 131 -type -beamUniform -6243.5
eleLoad -ele 132 -type -beamUniform -6243.5
eleLoad -ele 133 -type -beamUniform -6243.5
eleLoad -ele 134 -type -beamUniform -6243.5
eleLoad -ele 135 -type -beamUniform -6243.5
eleLoad -ele 136 -type -beamUniform -6243.5
eleLoad -ele 137 -type -beamUniform -6243.5
eleLoad -ele 138 -type -beamUniform -6243.5
eleLoad -ele 139 -type -beamUniform -6243.5
eleLoad -ele 140 -type -beamUniform -6243.5
eleLoad -ele 141 -type -beamUniform -6243.5
eleLoad -ele 142 -type -beamUniform -6243.5
eleLoad -ele 143 -type -beamUniform -6243.5
eleLoad -ele 144 -type -beamUniform -6243.5
eleLoad -ele 145 -type -beamUniform -6243.5
eleLoad -ele 146 -type -beamUniform -6243.5
eleLoad -ele 147 -type -beamUniform -6243.5
}
#
# apply GRAVITY-- # apply gravity load, set it constant and reset time to zero, load pattern has already been defined
puts goGravity
# Gravity-analysis parameters -- load-controlled static analysis
set Tol 1.e-4; # convergence tolerance for test
variable constraintsTypeGravity Plain; # default;
constraints $constraintsTypeGravity ; # how it handles boundary conditions
numberer RCM; # renumber dof's to minimize band-width (optimization), if you want to
system BandGeneral ; # how to store and solve the system of equations in the analysis (large model: try UmfPack)
test EnergyIncr $Tol 6 ; # determine if convergence has been achieved at the end of an iteration step
algorithm Newton; # use Newton's solution algorithm: updates tangent stiffness at every iteration
set NstepGravity 10; # apply gravity in 10 steps
set DGravity [expr 1./$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; # apply gravity
# ------------------------------------------------- maintain constant gravity loads and reset time to zero
loadConst -time 0.0
# -------------------------------------------------------------
#puts "Model Built"
set negen 1
set lambda [eigen [expr $negen]]
puts "lambda= $lambda"
set lambda [expr abs($lambda)]
set omega [expr pow($lambda,0.5)]
set tperiod [expr 2*3.1416/$omega]
puts " T= $tperiod "
i have two problem, first my script failed to converge in all iteration in low level of scaled acceleration and second dont show steress strain response for displacemet based element of shear wall.
what is your opinion about my problem?
> How did you use disp beamcolumnInt for shear wall ( as a big column??) and
> how did you connect it with the frame?
hi gerber
i have used disp beamcolumnInt similar what in opensees manual said for using this command. i have not defined connection shear wall whit frame. how i can defin this? can you help me that solve my problem
it is apart of my script that consist of modeling, defining recorders and static analysis parameters:
# --------------------------------------------------------------------------------------------------
#
# SET UP ----------------------------------------------------------------------------
wipe; # clear memory of all past model definitions
model BasicBuilder -ndm 2 -ndf 3; # Define the model builder, ndm=#dimension, ndf=#dofs
set RecordDir $Record; # set up name of data directory (can remove this)
file mkdir $RecordDir;
#file mkdir $RecordDir/drifts
file mkdir $RecordDir/Accel/drift/$acc
file mkdir $RecordDir/Accel/disp/$acc
file mkdir $RecordDir/Accel/envelopdrif/$acc
file mkdir $RecordDir/Accel/stressStrain/$acc
file mkdir $RecordDir/Accel/stressStrain1003/$acc
# create data directory
set GMdir "C:/tcleditor/bin/GMfiles"; # ground-motion file directory
#source LibUnits.tcl; # define units
source DisplayPlane.tcl; # procedure for displaying a plane in model
source DisplayModel2D.tcl; # procedure for displaying 3D perspectives of model
# define GEOMETRY -------------------------------------------------------------
# define NODAL COORDINATES
node 1 0.0 0.0
node 2 5.0 0.0
node 3 7.5 0.0
node 4 10.0 0.0
node 5 15.0 0.0
node 6 0.0 3.0
node 7 5.0 3.0
node 8 7.5 3.0
node 9 10.0 3.0
node 10 15.0 3.0
node 11 0.0 6.5
node 12 5.0 6.5
node 13 7.5 6.5
node 14 10.0 6.5
node 15 15.0 6.5
node 16 0.0 10.0
node 17 5.0 10.0
node 18 7.5 10.0
node 19 10.0 10.0
node 20 15.0 10.0
node 21 0.0 13.5
node 22 5.0 13.5
node 23 7.5 13.5
node 24 10.0 13.5
node 25 15.0 13.5
node 26 0.0 17.0
node 27 5.0 17.0
node 28 7.5 17.0
node 29 10.0 17.0
node 30 15.0 17.0
node 31 0.0 20.5
node 32 5.0 20.5
node 33 7.5 20.5
node 34 10.0 20.5
node 35 15.0 20.5
node 36 0.0 24.0
node 37 5.0 24.0
node 38 7.5 24.0
node 39 10.0 24.0
node 40 15.0 24.0
node 41 0.0 27.5
node 42 5.0 27.5
node 43 7.5 27.5
node 44 10.0 27.5
node 45 15.0 27.5
node 46 0.0 31.0
node 47 5.0 31.0
node 48 7.5 31.0
node 49 10.0 31.0
node 50 15.0 31.0
node 51 0.0 34.5
node 52 5.0 34.5
node 53 7.5 34.5
node 54 10.0 34.5
node 55 15.0 34.5
node 56 0.0 38.0
node 57 5.0 38.0
node 58 7.5 38.0
node 59 10.0 38.0
node 60 15.0 38.0
node 61 0.0 41.5
node 62 5.0 41.5
node 63 7.5 41.5
node 64 10.0 41.5
node 65 15.0 41.5
# Nodes for shear wall in story 1:
node 101 7.5 0.3
node 102 7.5 0.6
node 103 7.5 0.9
node 104 7.5 1.2
node 105 7.5 1.5
node 106 7.5 1.8
node 107 7.5 2.1
node 108 7.5 2.4
node 109 7.5 2.7
# Nodes for shear wall in story 2:
node 201 7.5 3.3
node 202 7.5 3.6
node 203 7.5 3.9
node 204 7.5 4.2
node 205 7.5 4.5
node 206 7.5 4.8
node 207 7.5 5.1
node 208 7.5 5.4
node 209 7.5 5.7
node 210 7.5 6.0
node 211 7.5 6.3
# Nodes for shear wall in story 3:
node 301 7.5 6.8
node 302 7.5 7.1
node 303 7.5 7.4
node 304 7.5 7.7
node 305 7.5 8.0
node 306 7.5 8.3
node 307 7.5 8.6
node 308 7.5 8.9
node 309 7.5 9.2
node 310 7.5 9.5
node 311 7.5 9.8
# Nodes for shear wall in story 4:
node 401 7.5 10.3
node 402 7.5 10.6
node 403 7.5 10.9
node 404 7.5 11.2
node 405 7.5 11.5
node 406 7.5 11.8
node 407 7.5 12.1
node 408 7.5 12.4
node 409 7.5 12.7
node 410 7.5 13.0
node 411 7.5 13.3
# Nodes for shear wall in story 5:
node 501 7.5 13.8
node 502 7.5 14.1
node 503 7.5 14.4
node 504 7.5 14.7
node 505 7.5 15.0
node 506 7.5 15.3
node 507 7.5 15.6
node 508 7.5 15.9
node 509 7.5 16.2
node 510 7.5 16.5
node 511 7.5 16.8
# Nodes for shear wall in story 6:
node 601 7.5 17.3
node 602 7.5 17.6
node 603 7.5 17.9
node 604 7.5 18.2
node 605 7.5 18.5
node 606 7.5 18.8
node 607 7.5 19.1
node 608 7.5 19.4
node 609 7.5 19.7
node 610 7.5 20.0
node 611 7.5 20.3
# Nodes for shear wall in story 7:
node 701 7.5 20.8
node 702 7.5 21.1
node 703 7.5 21.4
node 704 7.5 21.7
node 705 7.5 22.0
node 706 7.5 22.3
node 707 7.5 22.6
node 708 7.5 22.9
node 709 7.5 23.2
node 710 7.5 23.5
node 711 7.5 23.8
# Nodes for shear wall in story 8:
node 801 7.5 24.3
node 802 7.5 24.6
node 803 7.5 24.9
node 804 7.5 25.2
node 805 7.5 25.5
node 806 7.5 25.8
node 807 7.5 26.1
node 808 7.5 26.4
node 809 7.5 26.7
node 810 7.5 27.0
node 811 7.5 27.3
# Nodes for shear wall in story 9:
node 901 7.5 27.8
node 902 7.5 28.1
node 903 7.5 28.4
node 904 7.5 28.7
node 905 7.5 29.0
node 906 7.5 29.3
node 907 7.5 29.6
node 908 7.5 29.9
node 909 7.5 30.2
node 910 7.5 30.5
node 911 7.5 30.8
# Nodes for shear wall in story 10:
node 1001 7.5 31.3
node 1002 7.5 31.6
node 1003 7.5 31.9
node 1004 7.5 32.2
node 1005 7.5 32.5
node 1006 7.5 32.8
node 1007 7.5 33.1
node 1008 7.5 33.4
node 1009 7.5 33.7
node 1010 7.5 34.0
node 1011 7.5 34.3
# Nodes for shear wall in story 11:
node 1101 7.5 34.8
node 1102 7.5 35.1
node 1103 7.5 35.4
node 1104 7.5 35.7
node 1105 7.5 36.0
node 1106 7.5 36.3
node 1107 7.5 36.6
node 1108 7.5 36.9
node 1109 7.5 37.2
node 1110 7.5 37.5
node 1111 7.5 37.8
# Nodes for shear wall in story 12:
node 1201 7.5 38.3
node 1202 7.5 38.6
node 1203 7.5 38.9
node 1204 7.5 39.2
node 1205 7.5 39.5
node 1206 7.5 39.8
node 1207 7.5 40.1
node 1208 7.5 40.4
node 1209 7.5 40.7
node 1210 7.5 41.0
node 1211 7.5 41.3
fix 1 1 1 1
fix 2 1 1 1
fix 3 1 1 1
fix 4 1 1 1
fix 5 1 1 1
set NStory 12
# determine support nodes where ground motions are input, for multiple-support excitation
set iSupportNode "1 2 3 4 5"
# calculated MODEL PARAMETERS, particular to this model
# Set up parameters that are particular to the model for displacement control
set IDctrlNode 63.; # node where displacement is read for displacement control
set IDctrlDOF 1; # degree of freedom of displacement read for displacement control
#set LBuilding [expr 36*$in] # total building height
#material for frame Elements :
# Define Steel property variables:
set Fy 240e5
set E 210e8
set b 0.03; #strain-hardening ratio (ratio between post-yield tangent and initial elastic tangent)
set cR0 20; #$R0=between 10 and 20
set cR1 0.925
set cR2 0.15
set a1 0.0005
set a2 0.01
set a3 0.0005
set a4 0.01
#
# ----- Columns & Beams Material -------- #
#
uniaxialMaterial Steel02 1 $Fy $E $b $cR0 $cR1 $cR2 $a1 $a2 $a3 $a4
#
#material for shear wall reinforcement:
#
uniaxialMaterial Steel02 1003 400e5 $E $b $cR0 $cR1 $cR2 $a1 $a2 $a3 $a4
#
# --------------------------------------- #
# Define Concrete property variables:
#
set fpc -210e4
set epsc0 -0.003
set fpcu -180e4
set epsU -0.004
# --------- ShearWall Material --------- #
#
uniaxialMaterial Concrete01 2 $fpc $epsc0 $fpcu $epsU
# Define Column & Beams Sections :
# --------------------------------------- #
#
#
# I Section
proc Wsection { secID matID d tw bf tf nfdw nftw nfbf nftf} {
set dw [expr ($d-2*$tf)*1e-2]
set y1 [expr (-$d*0.5)*1e-2]
set y2 [expr (-$dw*0.5)*1e-2]
set y3 [expr ($dw*0.5)*1e-2]
set y4 [expr ($d*0.5)*1e-2]
set z1 [expr (-$bf*0.5)*1e-2]
set z2 [expr (-$tw*0.5)*1e-2]
set z3 [expr ($tw*0.5)*1e-2]
set z4 [expr ($bf*0.5)*1e-2]
#
section Fiber $secID {
patch quad $matID $nftf $nfbf $y1 $z1 $y2 $z1 $y2 $z4 $y1 $z4
patch quad $matID $nfdw $nftw $y2 $z2 $y3 $z2 $y3 $z3 $y2 $z3
patch quad $matID $nftf $nfbf $y3 $z1 $y4 $z1 $y4 $z4 $y3 $z4
}
}
#
# Beam section:
# tag mat H tw bf tf nfdw nftw nfbf nftf
Wsection 1 1 25 0.10 20.0 1.00 6 2 8 4; #PG1
#
# --------------------------------------- #
proc Isection { secID1 matID d tw bf tf nfdw nftw nfbf nftf} {
set dw [expr ($d-2*$tf)]
set z1 [expr (-$d/2)*1e-2]
set z2 [expr (-$dw/2)*1e-2]
set z3 [expr ($dw/2)*1e-2]
set z4 [expr ($d/2)*1e-2]
set y1 [expr (-$bf/2)*1e-2]
set y2 [expr (-$tw/2)*1e-2]
set y3 [expr ($tw/2)*1e-2]
set y4 [expr ($bf/2)*1e-2]
#
section fiberSec $secID1 {
patch quadr $matID $nfbf $nftf $y1 $z3 $y4 $z3 $y4 $z4 $y1 $z4
patch quadr $matID $nftw $nfdw $y2 $z2 $y3 $z2 $y3 $z3 $y2 $z3
patch quadr $matID $nfbf $nftf $y1 $z1 $y4 $z1 $y4 $z2 $y1 $z2
}
}
#
# Column section:
# tag mat H tw bf tf nfdw nftw nfbf nftf
Isection 2 1 12 0.65 12.0 1.10 6 2 8 4; #HE120B
#
Isection 3 1 14 0.7 14.0 1.20 6 2 8 4; #HE140B
#
Isection 4 1 16 0.8 16.0 1.30 6 2 8 4; #HE160B
#
Isection 5 1 18 0.85 18.0 1.40 6 2 8 4; #HE180B
#
Isection 6 1 20 0.9 20.0 1.50 6 2 8 4; #HE200B
#
Isection 7 1 22 0.95 22.0 1.60 6 2 8 4; #HE220B
#
Isection 8 1 24 1.0 24.0 1.70 6 2 8 4; #HE240B
#
Isection 9 1 28 1.05 28.0 1.80 6 2 8 4; #HE160B
#
Isection 10 1 32 1.15 30.0 2.05 6 2 8 4; #HE180B
#
Isection 11 1 34 1.2 30.0 2.15 6 2 8 4; #HE200B
#
Isection 12 1 36 1.25 30.0 2.25 6 2 8 4; #HE200B
#
#shear wall section
section FiberInt 15 -NStrip 2 0.5 10 0.5 2 0.5 {
#vertical fibers
fiber -2.5 0 0.015 2; fiber -2.5 0 0.00165 1003;
fiber -2.2 0 0.02 2; fiber -2.2 0 0.00165 1003;
fiber -1.8 0 0.02 2; fiber -1.8 0 0.00165 1003;
fiber -1.4 0 0.02 2; fiber -1.4 0 0.0022 1003;
fiber -1.0 0 0.02 2; fiber -1.0 0 0.0022 1003;
fiber -0.6 0 0.02 2; fiber -0.6 0 0.0022 1003;
fiber -0.2 0 0.01 2; fiber -0.2 0 0.0022 1003;
fiber 0.2 0 0.01 2; fiber 0.2 0 0.0022 1003;
fiber 0.6 0 0.02 2; fiber 0.6 0 0.0022 1003;
fiber 1.0 0 0.02 2; fiber 1.0 0 0.0022 1003;
fiber 1.4 0 0.02 2; fiber 1.4 0 0.0022 1003;
fiber 1.8 0 0.015 2; fiber 1.8 0 0.00165 1003;
fiber 2.2 0 0.02 2; fiber 2.2 0 0.00165 1003;
fiber 2.5 0 0.015 2; fiber 2.5 0 0.00165 1003;
#horiz. reinf.
Hfiber 0 0 0.000915 1003;
}
#element dispBeamColumnInt 111 2 111 1 5 3 0.1
geomTransf Linear 1; # beam
geomTransf Linear 2; # column
geomTransf LinearInt 3; # shear wall
# column elements
element nonlinearBeamColumn 1 1 6 5 8 1
element nonlinearBeamColumn 2 2 7 5 12 1
element nonlinearBeamColumn 3 4 9 5 12 1
element nonlinearBeamColumn 4 5 10 5 8 1
element nonlinearBeamColumn 5 6 11 5 8 1
element nonlinearBeamColumn 6 7 12 5 12 1
element nonlinearBeamColumn 7 9 14 5 12 1
element nonlinearBeamColumn 8 10 15 5 8 1
element nonlinearBeamColumn 9 11 16 5 7 1
element nonlinearBeamColumn 10 12 17 5 11 1
element nonlinearBeamColumn 11 14 19 5 11 1
element nonlinearBeamColumn 12 15 20 5 7 1
element nonlinearBeamColumn 13 16 21 5 7 1
element nonlinearBeamColumn 14 17 22 5 11 1
element nonlinearBeamColumn 15 19 24 5 11 1
element nonlinearBeamColumn 16 20 25 5 7 1
element nonlinearBeamColumn 17 21 26 5 6 1
element nonlinearBeamColumn 18 22 27 5 10 1
element nonlinearBeamColumn 19 24 29 5 10 1
element nonlinearBeamColumn 20 25 30 5 6 1
element nonlinearBeamColumn 21 26 31 5 6 1
element nonlinearBeamColumn 22 27 32 5 10 1
element nonlinearBeamColumn 23 29 34 5 10 1
element nonlinearBeamColumn 24 30 35 5 6 1
element nonlinearBeamColumn 25 31 36 5 5 1
element nonlinearBeamColumn 26 32 37 5 9 1
element nonlinearBeamColumn 27 34 39 5 9 1
element nonlinearBeamColumn 28 35 40 5 5 1
element nonlinearBeamColumn 29 36 41 5 5 1
element nonlinearBeamColumn 30 37 42 5 9 1
element nonlinearBeamColumn 31 39 44 5 9 1
element nonlinearBeamColumn 32 40 45 5 5 1
element nonlinearBeamColumn 33 41 46 5 3 1
element nonlinearBeamColumn 34 42 47 5 6 1
element nonlinearBeamColumn 35 44 49 5 6 1
element nonlinearBeamColumn 36 45 50 5 3 1
element nonlinearBeamColumn 37 46 51 5 3 1
element nonlinearBeamColumn 38 47 52 5 6 1
element nonlinearBeamColumn 39 49 54 5 6 1
element nonlinearBeamColumn 40 50 55 5 3 1
element nonlinearBeamColumn 41 51 56 5 2 1
element nonlinearBeamColumn 42 52 57 5 4 1
element nonlinearBeamColumn 43 54 59 5 4 1
element nonlinearBeamColumn 44 55 60 5 2 1
element nonlinearBeamColumn 45 56 61 5 2 1
element nonlinearBeamColumn 46 57 62 5 4 1
element nonlinearBeamColumn 47 59 64 5 4 1
element nonlinearBeamColumn 48 60 65 5 2 1
# beam elements
element nonlinearBeamColumn 100 6 7 5 1 1
element nonlinearBeamColumn 101 7 8 5 1 1
element nonlinearBeamColumn 102 8 9 5 1 1
element nonlinearBeamColumn 103 9 10 5 1 1
element nonlinearBeamColumn 104 11 12 5 1 1
element nonlinearBeamColumn 105 12 13 5 1 1
element nonlinearBeamColumn 106 13 14 5 1 1
element nonlinearBeamColumn 107 14 15 5 1 1
element nonlinearBeamColumn 108 16 17 5 1 1
element nonlinearBeamColumn 109 17 18 5 1 1
element nonlinearBeamColumn 110 18 19 5 1 1
element nonlinearBeamColumn 111 19 20 5 1 1
element nonlinearBeamColumn 112 21 22 5 1 1
element nonlinearBeamColumn 113 22 23 5 1 1
element nonlinearBeamColumn 114 23 24 5 1 1
element nonlinearBeamColumn 115 24 25 5 1 1
element nonlinearBeamColumn 116 26 27 5 1 1
element nonlinearBeamColumn 117 27 28 5 1 1
element nonlinearBeamColumn 118 28 29 5 1 1
element nonlinearBeamColumn 119 29 30 5 1 1
element nonlinearBeamColumn 120 31 32 5 1 1
element nonlinearBeamColumn 121 32 33 5 1 1
element nonlinearBeamColumn 122 33 34 5 1 1
element nonlinearBeamColumn 123 34 35 5 1 1
element nonlinearBeamColumn 124 36 37 5 1 1
element nonlinearBeamColumn 125 37 38 5 1 1
element nonlinearBeamColumn 126 38 39 5 1 1
element nonlinearBeamColumn 127 39 40 5 1 1
element nonlinearBeamColumn 128 41 42 5 1 1
element nonlinearBeamColumn 129 42 43 5 1 1
element nonlinearBeamColumn 130 43 44 5 1 1
element nonlinearBeamColumn 131 44 45 5 1 1
element nonlinearBeamColumn 132 46 47 5 1 1
element nonlinearBeamColumn 133 47 48 5 1 1
element nonlinearBeamColumn 134 48 49 5 1 1
element nonlinearBeamColumn 135 49 50 5 1 1
element nonlinearBeamColumn 136 51 52 5 1 1
element nonlinearBeamColumn 137 52 53 5 1 1
element nonlinearBeamColumn 138 53 54 5 1 1
element nonlinearBeamColumn 139 54 55 5 1 1
element nonlinearBeamColumn 140 56 57 5 1 1
element nonlinearBeamColumn 141 57 58 5 1 1
element nonlinearBeamColumn 142 58 59 5 1 1
element nonlinearBeamColumn 143 59 60 5 1 1
element nonlinearBeamColumn 144 61 62 5 1 1
element nonlinearBeamColumn 145 62 63 5 1 1
element nonlinearBeamColumn 146 63 64 5 1 1
element nonlinearBeamColumn 147 64 65 5 1 1
# shear wall elements in story 1
element dispBeamColumnInt 1001 3 101 1 15 3 0.4
element dispBeamColumnInt 1002 101 102 1 15 3 0.4
element dispBeamColumnInt 1003 102 103 1 15 3 0.4
element dispBeamColumnInt 1004 103 104 1 15 3 0.4
element dispBeamColumnInt 1005 104 105 1 15 3 0.4
element dispBeamColumnInt 1006 105 106 1 15 3 0.4
element dispBeamColumnInt 1007 106 107 1 15 3 0.4
element dispBeamColumnInt 1008 107 108 1 15 3 0.4
element dispBeamColumnInt 1009 108 109 1 15 3 0.4
element dispBeamColumnInt 1010 109 8 1 15 3 0.4
# shear wall elements in story 2
element dispBeamColumnInt 2001 8 201 1 15 3 0.4
element dispBeamColumnInt 2002 201 202 1 15 3 0.4
element dispBeamColumnInt 2003 202 203 1 15 3 0.4
element dispBeamColumnInt 2004 203 204 1 15 3 0.4
element dispBeamColumnInt 2005 204 205 1 15 3 0.4
element dispBeamColumnInt 2006 205 206 1 15 3 0.4
element dispBeamColumnInt 2007 206 207 1 15 3 0.4
element dispBeamColumnInt 2008 207 208 1 15 3 0.4
element dispBeamColumnInt 2009 208 209 1 15 3 0.4
element dispBeamColumnInt 2010 209 210 1 15 3 0.4
element dispBeamColumnInt 2011 210 211 1 15 3 0.4
element dispBeamColumnInt 2012 211 13 1 15 3 0.4
# shear wall elements in story 3
element dispBeamColumnInt 3001 13 301 1 15 3 0.4
element dispBeamColumnInt 3002 301 302 1 15 3 0.4
element dispBeamColumnInt 3003 302 303 1 15 3 0.4
element dispBeamColumnInt 3004 303 304 1 15 3 0.4
element dispBeamColumnInt 3005 304 305 1 15 3 0.4
element dispBeamColumnInt 3006 305 306 1 15 3 0.4
element dispBeamColumnInt 3007 306 307 1 15 3 0.4
element dispBeamColumnInt 3008 307 308 1 15 3 0.4
element dispBeamColumnInt 3009 308 309 1 15 3 0.4
element dispBeamColumnInt 3010 309 310 1 15 3 0.4
element dispBeamColumnInt 3011 310 311 1 15 3 0.4
element dispBeamColumnInt 3012 311 18 1 15 3 0.4
# shear wall elements in story 4
element dispBeamColumnInt 4001 18 401 1 15 3 0.4
element dispBeamColumnInt 4002 401 402 1 15 3 0.4
element dispBeamColumnInt 4003 402 403 1 15 3 0.4
element dispBeamColumnInt 4004 403 404 1 15 3 0.4
element dispBeamColumnInt 4005 404 405 1 15 3 0.4
element dispBeamColumnInt 4006 405 406 1 15 3 0.4
element dispBeamColumnInt 4007 406 407 1 15 3 0.4
element dispBeamColumnInt 4008 407 408 1 15 3 0.4
element dispBeamColumnInt 4009 408 409 1 15 3 0.4
element dispBeamColumnInt 4010 409 410 1 15 3 0.4
element dispBeamColumnInt 4011 410 411 1 15 3 0.4
element dispBeamColumnInt 4012 411 23 1 15 3 0.4
# shear wall elements in story 5
element dispBeamColumnInt 5001 23 501 1 15 3 0.4
element dispBeamColumnInt 5002 501 502 1 15 3 0.4
element dispBeamColumnInt 5003 502 503 1 15 3 0.4
element dispBeamColumnInt 5004 503 504 1 15 3 0.4
element dispBeamColumnInt 5005 504 505 1 15 3 0.4
element dispBeamColumnInt 5006 505 506 1 15 3 0.4
element dispBeamColumnInt 5007 506 507 1 15 3 0.4
element dispBeamColumnInt 5008 507 508 1 15 3 0.4
element dispBeamColumnInt 5009 508 509 1 15 3 0.4
element dispBeamColumnInt 5010 509 510 1 15 3 0.4
element dispBeamColumnInt 5011 510 511 1 15 3 0.4
element dispBeamColumnInt 5012 511 28 1 15 3 0.4
# shear wall elements in story 6
element dispBeamColumnInt 6001 28 601 1 15 3 0.4
element dispBeamColumnInt 6002 601 602 1 15 3 0.4
element dispBeamColumnInt 6003 602 603 1 15 3 0.4
element dispBeamColumnInt 6004 603 604 1 15 3 0.4
element dispBeamColumnInt 6005 604 605 1 15 3 0.4
element dispBeamColumnInt 6006 605 606 1 15 3 0.4
element dispBeamColumnInt 6007 606 607 1 15 3 0.4
element dispBeamColumnInt 6008 607 608 1 15 3 0.4
element dispBeamColumnInt 6009 608 609 1 15 3 0.4
element dispBeamColumnInt 6010 609 610 1 15 3 0.4
element dispBeamColumnInt 6011 610 611 1 15 3 0.4
element dispBeamColumnInt 6012 611 33 1 15 3 0.4
# shear wall elements in story 7
element dispBeamColumnInt 7001 33 701 1 15 3 0.4
element dispBeamColumnInt 7002 701 702 1 15 3 0.4
element dispBeamColumnInt 7003 702 703 1 15 3 0.4
element dispBeamColumnInt 7004 703 704 1 15 3 0.4
element dispBeamColumnInt 7005 704 705 1 15 3 0.4
element dispBeamColumnInt 7006 705 706 1 15 3 0.4
element dispBeamColumnInt 7007 706 707 1 15 3 0.4
element dispBeamColumnInt 7008 707 708 1 15 3 0.4
element dispBeamColumnInt 7009 708 709 1 15 3 0.4
element dispBeamColumnInt 7010 709 710 1 15 3 0.4
element dispBeamColumnInt 7011 710 711 1 15 3 0.4
element dispBeamColumnInt 7012 711 38 1 15 3 0.4
# shear wall elements in story 8
element dispBeamColumnInt 8001 38 801 1 15 3 0.4
element dispBeamColumnInt 8002 801 802 1 15 3 0.4
element dispBeamColumnInt 8003 802 803 1 15 3 0.4
element dispBeamColumnInt 8004 803 804 1 15 3 0.4
element dispBeamColumnInt 8005 804 805 1 15 3 0.4
element dispBeamColumnInt 8006 805 806 1 15 3 0.4
element dispBeamColumnInt 8007 806 807 1 15 3 0.4
element dispBeamColumnInt 8008 807 808 1 15 3 0.4
element dispBeamColumnInt 8009 808 809 1 15 3 0.4
element dispBeamColumnInt 8010 809 810 1 15 3 0.4
element dispBeamColumnInt 8011 810 811 1 15 3 0.4
element dispBeamColumnInt 8012 811 43 1 15 3 0.4
# shear wall elements in story 9
element dispBeamColumnInt 9001 43 901 1 15 3 0.4
element dispBeamColumnInt 9002 901 902 1 15 3 0.4
element dispBeamColumnInt 9003 902 903 1 15 3 0.4
element dispBeamColumnInt 9004 903 904 1 15 3 0.4
element dispBeamColumnInt 9005 904 905 1 15 3 0.4
element dispBeamColumnInt 9006 905 906 1 15 3 0.4
element dispBeamColumnInt 9007 906 907 1 15 3 0.4
element dispBeamColumnInt 9008 907 908 1 15 3 0.4
element dispBeamColumnInt 9009 908 909 1 15 3 0.4
element dispBeamColumnInt 9010 909 910 1 15 3 0.4
element dispBeamColumnInt 9011 910 911 1 15 3 0.4
element dispBeamColumnInt 9012 911 48 1 15 3 0.4
# shear wall elements in story 10
element dispBeamColumnInt 10001 48 1001 1 15 3 0.4
element dispBeamColumnInt 10002 1001 1002 1 15 3 0.4
element dispBeamColumnInt 10003 1002 1003 1 15 3 0.4
element dispBeamColumnInt 10004 1003 1004 1 15 3 0.4
element dispBeamColumnInt 10005 1004 1005 1 15 3 0.4
element dispBeamColumnInt 10006 1005 1006 1 15 3 0.4
element dispBeamColumnInt 10007 1006 1007 1 15 3 0.4
element dispBeamColumnInt 10008 1007 1008 1 15 3 0.4
element dispBeamColumnInt 10009 1008 1009 1 15 3 0.4
element dispBeamColumnInt 10010 1009 1010 1 15 3 0.4
element dispBeamColumnInt 10011 1010 1011 1 15 3 0.4
element dispBeamColumnInt 10012 1011 53 1 15 3 0.4
# shear wall elements in story 11
element dispBeamColumnInt 11001 53 1101 1 15 3 0.4
element dispBeamColumnInt 11002 1101 1102 1 15 3 0.4
element dispBeamColumnInt 11003 1102 1103 1 15 3 0.4
element dispBeamColumnInt 11004 1103 1104 1 15 3 0.4
element dispBeamColumnInt 11005 1104 1105 1 15 3 0.4
element dispBeamColumnInt 11006 1105 1106 1 15 3 0.4
element dispBeamColumnInt 11007 1106 1107 1 15 3 0.4
element dispBeamColumnInt 11008 1107 1108 1 15 3 0.4
element dispBeamColumnInt 11009 1108 1109 1 15 3 0.4
element dispBeamColumnInt 11010 1109 1110 1 15 3 0.4
element dispBeamColumnInt 11011 1110 1111 1 15 3 0.4
element dispBeamColumnInt 11012 1111 58 1 15 3 0.4
# shear wall elements in story 12
element dispBeamColumnInt 12001 58 1201 1 15 3 0.4
element dispBeamColumnInt 12002 1201 1202 1 15 3 0.4
element dispBeamColumnInt 12003 1202 1203 1 15 3 0.4
element dispBeamColumnInt 12004 1203 1204 1 15 3 0.4
element dispBeamColumnInt 12005 1204 1205 1 15 3 0.4
element dispBeamColumnInt 12006 1205 1206 1 15 3 0.4
element dispBeamColumnInt 12007 1206 1207 1 15 3 0.
element dispBeamColumnInt 12008 1207 1208 1 15 3 0.4
element dispBeamColumnInt 12009 1208 1209 1 15 3 0.4
element dispBeamColumnInt 12010 1209 1210 1 15 3 0.4
element dispBeamColumnInt 12011 1210 1211 1 15 3 0.4
element dispBeamColumnInt 12012 1211 63 1 15 3 0.4
# --------------------------------------------------------------------------------------------------------------------------------
# Define GRAVITY LOADS, weight and masses
set iFloorWeight "93653 93653 93653 93653 93653 93653 93653 93653 93653 93653 93653 89169"
set WeightTotal 1489480
mass 6 1592.7 0 0
mass 7 3185.5 0 0
mass 9 3185.5 0 0
mass 10 1592.7 0 0
mass 11 1592.7 0 0
mass 12 3185.5 0 0
mass 14 3185.5 0 0
mass 15 1592.7 0 0
mass 16 1592.7 0 0
mass 17 3185.5 0 0
mass 19 3185.5 0 0
mass 20 1592.7 0 0
mass 21 1592.7 0 0
mass 22 3185.5 0 0
mass 24 3185.5 0 0
mass 25 1592.7 0 0
mass 26 1592.7 0 0
mass 27 3185.5 0 0
mass 29 3185.5 0 0
mass 30 1592.7 0 0
mass 31 1592.7 0 0
mass 32 3185.5 0 0
mass 34 3185.5 0 0
mass 35 1592.7 0 0
mass 36 1592.7 0 0
mass 37 3185.5 0 0
mass 39 3185.5 0 0
mass 40 1592.7 0 0
mass 41 1592.7 0 0
mass 42 3185.5 0 0
mass 44 3185.5 0 0
mass 45 1592.7 0 0
mass 46 1592.7 0 0
mass 47 3185.5 0 0
mass 49 3185.5 0 0
mass 50 1592.7 0 0
mass 51 1592.7 0 0
mass 52 3185.5 0 0
mass 54 3185.5 0 0
mass 55 1592.7 0 0
mass 56 1592.7 0 0
mass 57 3185.5 0 0
mass 59 3185.5 0 0
mass 60 1592.7 0 0
mass 61 1592.7 0 0
mass 62 3185.5 0 0
mass 64 3185.5 0 0
mass 65 1592.7 0 0
#mass 4 $MassNode $MassNode $MassNode 0. 0. 0.; # define mass
set MassTotal [expr $WeightTotal/9.806];
# --------------------------------------------------------------------------------------------------------------------------------
# Define RECORDERS -------------------------------------------------------------
recorder Node -file $RecordDir/DFree.out -time -node 63 -dof 1 2 3 accel;
recorder Drift -file $RecordDir/Accel/drift/$acc/1.out -iNode 3 -jNode 8 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/1.out -node 8 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/1.out -iNode 3 -jNode 8 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/1.out -time -ele 1006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/1.out -time -ele 1006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/2.out -iNode 8 -jNode 13 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/2.out -node 13 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/2.out -iNode 8 -jNode 13 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/2.out -time -ele 2006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/2.out -time -ele 2006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/3.out -iNode 13 -jNode 18 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/3.out -node 18 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/3.out -iNode 13 -jNode 18 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/3.out -time -ele 3006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/3.out -time -ele 3006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/4.out -iNode 18 -jNode 23 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/4.out -node 23 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/4.out -iNode 18 -jNode 23 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/4.out -time -ele 4006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/4.out -time -ele 4006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/5.out -iNode 23 -jNode 28 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/5.out -node 28 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/5.out -iNode 23 -jNode 28 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/5.out -time -ele 5006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/5.out -time -ele 5006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/6.out -iNode 28 -jNode 33 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/6.out -node 33 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/6.out -iNode 28 -jNode 33 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/6.out -time -ele 6006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/6.out -time -ele 6006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/7.out -iNode 33 -jNode 38 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/7.out -node 38 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/7.out -iNode 33 -jNode 38 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/7.out -time -ele 7006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/7.out -time -ele 7006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/8.out -iNode 38 -jNode 43 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/8.out -node 43 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/8.out -iNode 38 -jNode 43 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/8.out -time -ele 8006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/8.out -time -ele 8006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/9.out -iNode 43 -jNode 48 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/9.out -node 48 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/9.out -iNode 43 -jNode 48 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/9.out -time -ele 9006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/9.out -time -ele 9006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/10.out -iNode 48 -jNode 53 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/10.out -node 53 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/10.out -iNode 48 -jNode 53 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/10.out -time -ele 10006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/10.out -time -ele 10006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/11.out -iNode 53 -jNode 58 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/11.out -node 58 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/11.out -iNode 53 -jNode 58 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/11.out -time -ele 11006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/11.out -time -ele 11006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/Accel/drift/$acc/12.out -iNode 58 -jNode 63 -dof 1 -perpDirn 2
recorder Node -file $RecordDir/Accel/disp/$acc/12.out -node 63 -dof 1 disp
recorder EnvelopeDrift -file $RecordDir/Accel/envelopdrif/$acc/12.out -iNode 58 -jNode 63 -dof 1 -perpDirn 2
recorder Element -file $RecordDir/Accel/stressStrain/$acc/12.out -time -ele 12006 section 15 fiber 2.5 0 2 stressStrain
recorder Element -file $RecordDir/Accel/stressStrain1003/$acc/12.out -time -ele 12006 section 15 fiber 2.5 0 1003 stressStrain
recorder Drift -file $RecordDir/globaldrift.out -iNode 63 -jNode 3 -dof 1 -perpDirn 2
# Define DISPLAY -------------------------------------------------------------
#DisplayModel2D DeformedShape ; # options: DeformedShape NodeNumbers ModeShape
# GRAVITY -------------------------------------------------------------
# define GRAVITY load applied to beams and columns -- eleLoad applies loads in local coordinate axis
#pattern Plain 101 Linear {
#
pattern Plain 1 "Linear" {
eleLoad -ele 100 -type -beamUniform -6243.5
eleLoad -ele 101 -type -beamUniform -6243.5
eleLoad -ele 102 -type -beamUniform -6243.5
eleLoad -ele 103 -type -beamUniform -6243.5
eleLoad -ele 104 -type -beamUniform -6243.5
eleLoad -ele 105 -type -beamUniform -6243.5
eleLoad -ele 106 -type -beamUniform -6243.5
eleLoad -ele 107 -type -beamUniform -6243.5
eleLoad -ele 108 -type -beamUniform -6243.5
eleLoad -ele 109 -type -beamUniform -6243.5
eleLoad -ele 110 -type -beamUniform -6243.5
eleLoad -ele 111 -type -beamUniform -6243.5
eleLoad -ele 112 -type -beamUniform -6243.5
eleLoad -ele 113 -type -beamUniform -6243.5
eleLoad -ele 114 -type -beamUniform -6243.5
eleLoad -ele 115 -type -beamUniform -6243.5
eleLoad -ele 116 -type -beamUniform -6243.5
eleLoad -ele 117 -type -beamUniform -6243.5
eleLoad -ele 118 -type -beamUniform -6243.5
eleLoad -ele 119 -type -beamUniform -6243.5
eleLoad -ele 120 -type -beamUniform -6243.5
eleLoad -ele 121 -type -beamUniform -6243.5
eleLoad -ele 122 -type -beamUniform -6243.5
eleLoad -ele 123 -type -beamUniform -6243.5
eleLoad -ele 124 -type -beamUniform -6243.5
eleLoad -ele 125 -type -beamUniform -6243.5
eleLoad -ele 126 -type -beamUniform -6243.5
eleLoad -ele 127 -type -beamUniform -6243.5
eleLoad -ele 128 -type -beamUniform -6243.5
eleLoad -ele 129 -type -beamUniform -6243.5
eleLoad -ele 130 -type -beamUniform -6243.5
eleLoad -ele 131 -type -beamUniform -6243.5
eleLoad -ele 132 -type -beamUniform -6243.5
eleLoad -ele 133 -type -beamUniform -6243.5
eleLoad -ele 134 -type -beamUniform -6243.5
eleLoad -ele 135 -type -beamUniform -6243.5
eleLoad -ele 136 -type -beamUniform -6243.5
eleLoad -ele 137 -type -beamUniform -6243.5
eleLoad -ele 138 -type -beamUniform -6243.5
eleLoad -ele 139 -type -beamUniform -6243.5
eleLoad -ele 140 -type -beamUniform -6243.5
eleLoad -ele 141 -type -beamUniform -6243.5
eleLoad -ele 142 -type -beamUniform -6243.5
eleLoad -ele 143 -type -beamUniform -6243.5
eleLoad -ele 144 -type -beamUniform -6243.5
eleLoad -ele 145 -type -beamUniform -6243.5
eleLoad -ele 146 -type -beamUniform -6243.5
eleLoad -ele 147 -type -beamUniform -6243.5
}
#
# apply GRAVITY-- # apply gravity load, set it constant and reset time to zero, load pattern has already been defined
puts goGravity
# Gravity-analysis parameters -- load-controlled static analysis
set Tol 1.e-4; # convergence tolerance for test
variable constraintsTypeGravity Plain; # default;
constraints $constraintsTypeGravity ; # how it handles boundary conditions
numberer RCM; # renumber dof's to minimize band-width (optimization), if you want to
system BandGeneral ; # how to store and solve the system of equations in the analysis (large model: try UmfPack)
test EnergyIncr $Tol 6 ; # determine if convergence has been achieved at the end of an iteration step
algorithm Newton; # use Newton's solution algorithm: updates tangent stiffness at every iteration
set NstepGravity 10; # apply gravity in 10 steps
set DGravity [expr 1./$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; # apply gravity
# ------------------------------------------------- maintain constant gravity loads and reset time to zero
loadConst -time 0.0
# -------------------------------------------------------------
#puts "Model Built"
set negen 1
set lambda [eigen [expr $negen]]
puts "lambda= $lambda"
set lambda [expr abs($lambda)]
set omega [expr pow($lambda,0.5)]
set tperiod [expr 2*3.1416/$omega]
puts " T= $tperiod "
i have two problem, first my script failed to converge in all iteration in low level of scaled acceleration and second dont show steress strain response for displacemet based element of shear wall.
what is your opinion about my problem?
Re: recorder elements for shear wall
your recorders fail because you don't have 15 sections in your element!
-
- Posts: 1
- Joined: Sun Jul 03, 2011 6:33 am
- Location: Iran University of Science and Technology
Re: recorder elements problem
alihaji1365 wrote:
>
> #shear wall section
> section FiberInt 15 -NStrip 2 0.5 10 0.5 2 0.5 {
> #vertical fibers
> fiber -2.5 0 0.015 2; fiber -2.5 0 0.00165 1003;
> fiber -2.2 0 0.02 2; fiber -2.2 0 0.00165 1003;
>
> fiber -1.8 0 0.02 2; fiber -1.8 0 0.00165 1003;
> fiber -1.4 0 0.02 2; fiber -1.4 0 0.0022 1003;
> fiber -1.0 0 0.02 2; fiber -1.0 0 0.0022 1003;
> fiber -0.6 0 0.02 2; fiber -0.6 0 0.0022 1003;
> fiber -0.2 0 0.01 2; fiber -0.2 0 0.0022 1003;
> fiber 0.2 0 0.01 2; fiber 0.2 0 0.0022 1003;
> fiber 0.6 0 0.02 2; fiber 0.6 0 0.0022 1003;
> fiber 1.0 0 0.02 2; fiber 1.0 0 0.0022 1003;
> fiber 1.4 0 0.02 2; fiber 1.4 0 0.0022 1003;
> fiber 1.8 0 0.015 2; fiber 1.8 0 0.00165 1003;
>
> fiber 2.2 0 0.02 2; fiber 2.2 0 0.00165 1003;
> fiber 2.5 0 0.015 2; fiber 2.5 0 0.00165 1003;
>
> #horiz. reinf.
> Hfiber 0 0 0.000915 1003;
> }
In your code, I think area of fiber for concrete should be 0.15 instead of 0.015, or 0.2 instead of 0.02.
>
> #shear wall section
> section FiberInt 15 -NStrip 2 0.5 10 0.5 2 0.5 {
> #vertical fibers
> fiber -2.5 0 0.015 2; fiber -2.5 0 0.00165 1003;
> fiber -2.2 0 0.02 2; fiber -2.2 0 0.00165 1003;
>
> fiber -1.8 0 0.02 2; fiber -1.8 0 0.00165 1003;
> fiber -1.4 0 0.02 2; fiber -1.4 0 0.0022 1003;
> fiber -1.0 0 0.02 2; fiber -1.0 0 0.0022 1003;
> fiber -0.6 0 0.02 2; fiber -0.6 0 0.0022 1003;
> fiber -0.2 0 0.01 2; fiber -0.2 0 0.0022 1003;
> fiber 0.2 0 0.01 2; fiber 0.2 0 0.0022 1003;
> fiber 0.6 0 0.02 2; fiber 0.6 0 0.0022 1003;
> fiber 1.0 0 0.02 2; fiber 1.0 0 0.0022 1003;
> fiber 1.4 0 0.02 2; fiber 1.4 0 0.0022 1003;
> fiber 1.8 0 0.015 2; fiber 1.8 0 0.00165 1003;
>
> fiber 2.2 0 0.02 2; fiber 2.2 0 0.00165 1003;
> fiber 2.5 0 0.015 2; fiber 2.5 0 0.00165 1003;
>
> #horiz. reinf.
> Hfiber 0 0 0.000915 1003;
> }
In your code, I think area of fiber for concrete should be 0.15 instead of 0.015, or 0.2 instead of 0.02.
-
- Posts: 1
- Joined: Fri Nov 02, 2012 1:13 am
- Location: Ain Shams University
Re: recorder elements for shear wall
i think there is something wrong with fiber area ie firs strip area is 0.15 not 0.05 according to your wall thickness t=0.5 , similarly second strip area 0.02