wipe
model basic -ndm 2 -ndf 3
if { [file exists oupt] == 0} {
file mkdir output
}
node 1 0 0
node 2 0 2
node 3 0 4
node 4 0 6
fix 1 1 1 1
uniaxialMaterial Concrete01 1 -26.8e6 -0.002 -10e6 -0.0033
uniaxialMaterial Concrete02 2 -40e6 -0.002 -10e6 -0.005 0.1 2.68 1000
uniaxialMaterial Steel02 3 400 206000 0.01 18 0.925 0.15 0 1 0 1
section Aggregator
section Fiber 1 {
patch rect 1 15 15 -0.675 -0.675 0.675 0.675
patch rect 2 15 1 -0.75 0.675 0.75 0.75
patch rect 2 15 1 -0.75 -0.75 0.75 -0.675
patch rect 2 2 13 -0.75 -0.675 -0.675 0.675
patch rect 2 2 13 0.675 -0.675 0.75 0.675
layer straight 3 11 8.01e-4 0.675 0.675 0.675 -0.675
layer straight 3 11 8.01e-4 0.64 0.64 0.64 -0.64
layer straight 3 11 8.01e-4 -0.675 0.675 -0.675 -0.675
layer straight 3 11 8.01e-4 -0.64 0.64 -0.64 -0.64
layer straight 3 9 2.01e-4 0.54 -0.675 -0.54 -0.675
layer straight 3 9 2.01e-4 0.54 -0.64 -0.54 -0.64
layer straight 3 9 2.01e-4 0.54 0.675 -0.54 0.675
layer straight 3 9 2.01e-4 0.54 0.64 -0.54 0.64
}
geomTransf Linear 1
element nonlinearBeamColumn 1 1 4 5 1 1
recorder Node -file output/disp_4.out -time -node 4 -dof 1 2 3 disp
recorder Node -file output/reaction_1.out -time -node 1 -dof 1 2 3 reaction
recorder Element -file output/force_1.txt -time -ele 1 globalForce
recorder Element -file output/defocolsec_1.txt -time -ele 1 section 1 deformation
recorder Element -file output/forcecolsec_1.txt -time -ele 1 section 1 force
pattern Plain 1 Linear {
load 4 0.0 -653860 0.0
}
constraints Plain
numberer Plain
system BandGeneral
test NormDispIncr 1.0e-6 6 2
algorithm Newton
integrator LoadControl 0.01
analysis Static
analyze 100
puts "重力分析完成..."