wipe
# Opensees dandun
#
#Units:kN, m, sec
# -----------------
# Start of model generation
# -----------------
# Create ModeBulider (with two-dimensions and 3 DOF/node)
model basic -ndm 2 -ndf 3
# -----------------
# tag X Y
node 1 0.0 0.0
node 2 0.0 0.0
node 3 0.0 2.0
辛沪光
node 4 0.0 4.0
node 5 0.0 6.0
node 6 0.0 8.0
node 7 0.0 10.0
ndm
node 8 0.0 12.0
node 9 0.0 14.0
node 10 0.0 16.0
node 11 0.0 18.0
node 12 0.0 20.0
# -----------------
# Fix supports at base of columns
# tag DX DY RZ
fix 1 1 1 1
# ----------------
# Concrete tag fc ec0 fcu ecu
# Core concrete (confined)
uniaxialMaterial Concrete01 1 -25600.0 -0.00219 -17780.0 -0.01
#Cover concrete (unconfined)
uniaxialMaterial Concrete01 2 -23400.0 -0.002 -0.0 -0.006
# STEEL
# Reinforcing steel
set fy 400000.0; #Yield stress
set E 200000000.0;# Young's modulus
# tag fy E0 b
uniaxialMaterial Steel02 3 $fy $E 0.01 18.5 0.925 0.15
uniaxialMaterial Elastic 11 29043600
uniaxialMaterial Elastic 12 12326600
uniaxialMaterial Elastic 13 587247596
#Define cross-section for nonlinear columns
# ---------------------
# set some parameters
set colWidth 8.18
set colDepth 4.28
set cover 0.05
set As 0.00049
# some variables derived from the parameters
set y1 [expr $colDepth/2.0]
set z1 [expr $colWidth/2.0]
section Fiber 1 {
# Create the concrete core fibers
patch rect 1 20 30 [expr $cover-$y1] [expr $cover-$z1] [expr $y1-$cover] [expr $z1-$cover]
# Create the concrete cover fibers (top, bottom, left, right)
patch rect 2 20 5 [expr -$y1] [expr $z1-$cover] $y1 $z1
patch rect 2 20 5 [expr -$y1] [expr -$z1] $y1 [expr $cover-$z1]
patch rect 2 5 10 [expr -$y1] [expr $cover-$z1] [expr $cover-$y1] [expr $z1-$cover]
patch rect 2 5 10 [expr $y1-$cover] [expr $cover-$z1] $y1 [expr $z1-$cover]
# Create the reinforcingfibers (left, middle, right)
layer straight 3 175 $As [expr $y1-$cover] [expr $z1-$cover] [expr $y1-$cover] [expr $cover-$z1]
layer straight 3 175 $As [expr $cover-$y1] [expr $z1-$cover] [expr $cover-$y1] [expr $cover-$z1]
layer straight 3 115 $As [expr $y1-$cover] [expr $z1-$cover] [expr $cover-$y1] [expr $z1-$cover]
layer straight 3 115 $As [expr $y1-$cover] [expr $cover-$z1] [expr $cover-$y1] [expr $cover-$z1]
}
# Define column elements
# ----------------------
# Geometry of column elements
# tag
geomTransf Linear 1
# Number of integration points along length of elements
set np 5
# Create the columns using Beam-column elements
中药制剂# tag ndI ndJ secID transfTag
element nonlinearBeamColumn 2 2 3 $np 1 1
element nonlinearBeamColumn 3 3 4 $np 1 1
element nonlinearBeamColumn 4 4 5 $np 1 1
element nonlinearBeamColumn 5 5 6 $np 1 1
element nonlinearBeamColumn 6 6 7 $np 1 1
牛田洋事件element nonlinearBeamColumn 7 7 8 $np 1 1 一帮一活动
element nonlinearBeamColumn 8 8 9 $np 1 1
element nonlinearBeamColumn 9 9 10 $np 1 1
element nonlinearBeamColumn 10 10 11 $np 1 1
element nonlinearBeamColumn 11 11 12 $np 1 1
equalDOF 1 2 1 2
element zeroLength 111 1 2 -mat 13 -dir 3
set m [expr 3355.7]
set n [expr 175.1]
# tag MX MY RZ
mass 2 $n $n 1e-8
mass 3 $n $n 1e-8
mass 4 $n $n 1e-8
mass 5 $n $n 1e-8
mass 6 $n $n 1e-8
mass 7 $n $n 1e-8
mass 8 $n $n 1e-8
mass 9 $n $n 1e-8
mass 10 $n $n 1e-8
mass 11 $n $n 1e-8
偶合症
mass 12 $m $m 1e-8
# Set a parameter for the axial load
set P 33557.0; # of axial capacity of columns
set Q 1715.5; # of axial capacity of columns
# Create a Plain load pattern with a Linear TimeSeries
pattern Plain 1 "Constant" {
# Create nodal loads at nodes 9
# nd FX FY MZ
load 12 0.0 [expr -$P] 0.0
load 2 0.0 [expr -$Q] 0.0
load 3 0.0 [expr -$Q] 0.0
load 4 0.0 [expr -$Q] 0.0
load 5 0.0 [expr -$Q] 0.0
load 6 0.0 [expr -$Q] 0.0
load 7 0.0 [expr -$Q] 0.0
load 8 0.0 [expr -$Q] 0.0
load 9 0.0 [expr -$Q] 0.0
load 10 0.0 [expr -$Q] 0.0
load 11 0.0 [expr -$Q] 0.0
}
system SparseGeneral -piv
# Create the constraint handler
constraints Transformation
# Create the time integration scheme
# Create the DOF numberer
numberer RCM
test NormDispIncr 1.0e-8 30 5
# Create the solution algorithm
algorithm Newton
# create the transient analysis
integrator LoadControl 1
# -----------------
# End of analysis generation
# -------------------
analysis Static
initialize