The OpenSees Community

Dear all
I am simulating a brige bent retrofitted with a self-centring brace.But I don't know how to model the SCB in the brige bent.I have tried to use tuss element to model SCB,but this method doesn't work.And here is my model code.
model BasicBuilder -ndm 2 -ndf 3
#L column
node 1 -1300 0
node 2 -1300 250
node 3 -1300 650
node 4 -1300 1050
node 5 -1300 1450
node 6 -1300 1850
#R column
node 11 1300 0
node 12 1300 250
node 13 1300 650
node 14 1300 1050
node 15 1300 1450
node 16 1300 1850
#capbeam
node 21 -1800 2050
node 22 -1300 2050
node 23 -1150 2050
node 24 1150 2050
node 25 1300 2050
node 26 1800 2050
#brace
node 31 -1300 0
node 32 1300 2050
#node 33 -1092 164
#node 34 1092 1886
#fix
fix 1 1 1 1
fix 11 1 1 1
#materail
#concrete r6
uniaxialMaterial Concrete02 1 -35.3 -0.002 -7.1 -0.004 0.10 3.70 3530
#concrete
uniaxialMaterial Concrete02 2 -45.9 -0.0026 -9.18 -0.015 0.10 4.22 3531
#bar
uniaxialMaterial Steel02 3 385 190000 1.00E-02 2.00E+01 9.25E-01 1.50E-01 0.01 1.00E+00 0.01 1.00E+00 0.00E+00
#SCB
uniaxialMaterial SelfCentering 4 104500 12100 397000 0.30 0 0 0
#
uniaxialMaterial Elastic 5 1
#uniaxialMaterial Elastic 6 1
section Fiber 1 {
# PatchCirc "con": matTag nfcirc nfrad y z intRad extRad <startAng endAng>
patch circ 1 10 10 +0.000000E+00 +0.000000E+00 +1.250000E+02 +1.5000000E+02 +0.000000E+00 +3.600000E+02
# PatchCirc "un": matTag nfcirc nfrad y z intRad extRad <startAng endAng>
patch circ 2 10 10 +0.000000E+00 +0.000000E+00 +0.000000E+02 +1.250000E+02 +0.000000E+00 +3.600000E+02
# LayerCircular "steel": matTag numBar areaBar y z radius <startAng endAng>
layer circ 3 10 +1.1309734E+02 +0.000000E+00 +0.000000E+00 +1.0700000E+02 +0.000000E+00 +3.600000E+02
}
geomTransf Linear 1
#Lcolumn
element beamWithHinges 1 2 3 1 229.64 1 229.64 37676.82 70685.83471 397607820 1 =-mass 0 =-iter 10 1E-12 =-integration Lobatto
element beamWithHinges 2 3 4 1 229.64 1 229.64 37676.82 70685.83471 397607820 1 =-mass 0 =-iter 10 1E-12 =-integration Lobatto
element beamWithHinges 3 4 5 1 229.64 1 229.64 37676.82 70685.83471 397607820 1 =-mass 0 =-iter 10 1E-12 =-integration Lobatto
element beamWithHinges 4 5 6 1 229.64 1 229.64 37676.82 70685.83471 397607820 1 =-mass 0 =-iter 10 1E-12 =-integration Lobatto
#Rcolumn
element beamWithHinges 11 12 13 1 229.64 1 229.64 37676.82 70685.83471 397607820 1 =-mass 0 =-iter 10 1E-12 =-integration Lobatto
element beamWithHinges 12 13 14 1 229.64 1 229.64 37676.82 70685.83471 397607820 1 =-mass 0 =-iter 10 1E-12 =-integration Lobatto
element beamWithHinges 13 14 15 1 229.64 1 229.64 37676.82 70685.83471 397607820 1 =-mass 0 =-iter 10 1E-12 =-integration Lobatto
element beamWithHinges 14 15 16 1 229.64 1 229.64 37676.82 70685.83471 397607820 1 =-mass 0 =-iter 10 1E-12 =-integration Lobatto
#rigidlink
element elasticBeamColumn 21 1 2 70685.8347 +5.000000E+10 +1.401985E+11 1 +0.000000E+00 +0.000000E+00 -mass +0.000000E+00
element elasticBeamColumn 22 11 12 70685.8347 +5.000000E+10 +1.401985E+11 1 +0.000000E+00 +0.000000E+00 -mass +0.000000E+00
#rigidlink
element elasticBeamColumn 31 6 22 70685.8347 +5.000000E+10 +1.401985E+11 1 +0.000000E+00 +0.000000E+00 -mass +0.000000E+00
element elasticBeamColumn 32 16 25 70685.8347 +5.000000E+10 +1.401985E+11 1 +0.000000E+00 +0.000000E+00 -mass +0.000000E+00
#capbeam
element elasticBeamColumn 41 21 22 +1.60000E+5 +3.096270E+04 2.13E+9 1 +0.000000E+00 +0.000000E+00 -mass +0.000000E+00
element elasticBeamColumn 42 22 23 +1.60000E+5 +5.000000E+10 2.13E+9 1 +0.000000E+00 +0.000000E+00 -mass +0.000000E+00
element elasticBeamColumn 43 23 24 +1.60000E+5 +3.096270E+04 2.13E+9 1 +0 +0 -mass +0
element elasticBeamColumn 44 24 25 +1.60000E+5 +5.000000E+10 2.13E+9 1 +0 +0 -mass +0
element elasticBeamColumn 45 25 26 +1.60000E+5 +3.096270E+04 2.13E+9 1 +0.000000E+00 +0.000000E+00 -mass +0.000000E+00
#define hingle
element zeroLength 61 1 31 -mat 5 -dir 6
element zeroLength 62 32 25 -mat 5 -dir 6
#define hingle
equalDOF 1 31 1 2
equalDOF 25 32 1 2
element elasticBeamColumn 71 31 33 70685.8347 +5.000000E+10 +1.401985E+11 1 +0.000000E+00 +0.000000E+00 -mass +0.000000E+00
element elasticBeamColumn 72 34 32 70685.8347 +5.000000E+10 +1.401985E+11 1 +0.000000E+00 +0.000000E+00 -mass +0.000000E+00
#define SCB
element truss 51 31 33 1 4 -rho +0.000000E+00 -cMass 0 -doRayleigh 0
recorder Node -file node21.out -time -node 21 -dof 1 2 3 disp
recorder Node -file node6.out -time -node 6 -dof 1 2 3 disp
recorder Node -file node32.out -time -node 32 -dof 1 2 3 disp
recorder Node -file node34.out -time -node 34 -dof 1 2 3 disp
recorder Node -file node25.out -time -node 25 -dof 1 2 3 disp
recorder Node -file node33.out -time -node 33 -dof 1 2 3 disp
recorder Node -file disp1.osp -time -node 1 2 3 4 5 6 11 12 13 14 15 16 21 22 23 24 25 26 31 32 33 34 -dof 1 2 3 disp
recorder Element -file ele0.out -time -ele 51 localForce
recorder Element -file ele1.out -ele 1 force
recorder Element -file ele1b.out -ele 1 section 1 deformation
pattern Plain 1 Linear {
load 22 0.0 -203893.2902 0.0
load 25 0.0 -203893.2902 0.0
}
constraints Transformation
numberer Plain
system BandGeneral
test EnergyIncr 1.0e-6 200
algorithm Newton
integrator LoadControl 1
analysis Static
analyze 1
loadConst 0.0
wipeAnalysis
pattern Plain 2 Linear {
load 21 1.0 0.0 0.0
}
source SmartAnalyze.tcl
constraints Transformation
numberer Plain
set protocol { 100 }
SmartAnalyzeStatic 21 1 0.1 $protocol

model BasicBuilder -ndm 2 -ndf 3
node 1 +0.000000E+00 +0.000000E+00 -ndf 3
node 2 +0.000000E+00 +0.000000E+00 -ndf 3
fix 1 1 1 1
fix 2 0 0 1
uniaxialMaterial ElasticPP 1 +1.580000E+03 +2.000000E-03 -2.000000E-03 +0.000000E+00
uniaxialMaterial Elastic 2 0 0
frictionModel VelPressureDep 1 0.085 0.163 7.0686 0.05 0.08 0.77
element flatSliderBearing 1 1 2 1 +1.580000E+03 -P 1 -Mz 2 -orient 0 1 0 -1 0 0
recorder Node -file AnalysisCase01_Node_Recorder01_Dsp.out -time -node 2 -dof 1 2 3 disp
recorder Node -file AnalysisCase01_Node_Recorder01_RFrc.out -time -node 2 -dof 1 2 3 reaction
recorder Node -file AnalysisCase01_Node_Recorder01_RFrcInertia.out -time -node 2 -dof 1 2 3 reactionIncludingInertia
pattern Plain 1 Linear {
load 2 +0.000000E+00 -5.566400E+05 +0.000000E+00
}
constraints Plain
numberer Plain
system BandGeneral
test EnergyIncr 1.0e-12 10
algorithm Newton
integrator LoadControl 0.02
analysis Static
analyze 10
loadConst 0
pattern Plain 2 Linear {
load 2 +1.000000E+00 +0.000000E+00 +0.000000E+00
}
constraints Plain
numberer Plain
system BandGeneral
test EnergyIncr 1.0e-12 10
algorithm Newton
integrator DisplacementControl 2 1 $protocol
set protocol {20 30 40 50 60 70 80 90}
analyze Static

WARNING BandGenLinLapackSolver::solve() -factorization failed, matrix singular U(i,i) = 0, i= 2
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at load factor 0.02
OpenSees > analyze failed, returned: -3 error flag

When I run the analysis, the error has appeared such as bellow,who can help me or give some suggestion? Thank you !

mhscott wrote: ↑Thu Sep 17, 2020 6:43 am There is an error somewhere in your section definition. Try elastic materials and draw out the fibers and back out the process until you find the error.

Thanks for your suggestion,it's very helpful!

This is my model code:
model BasicBuilder -ndm 2 -ndf 3
node 1 -1.300000E+03 +0.000000E+00 -ndf 3
node 2 +1.300000E+03 +0.000000E+00 -ndf 3
node 3 -1.300000E+03 +1.600000E+03 -ndf 3
node 4 +1.300000E+03 +1.600000E+03 -ndf 3
node 5 -1.300000E+03 +1.700000E+03 -ndf 3
node 6 +1.300000E+03 +1.700000E+03 -ndf 3
node 7 -1.150000E+03 +1.700000E+03 -ndf 3
node 8 +1.150000E+03 +1.700000E+03 -ndf 3
node 9 -1.800000E+03 +1.700000E+03 -ndf 3
node 10 +1.800000E+03 +1.700000E+03 -ndf 3
fix 1 1 1 1
fix 2 1 1 1
rigidLink beam 5 3
rigidLink beam 5 7
rigidLink beam 8 6
rigidLink beam 6 4
# Material "»ìÄýÍÁ": matTag fc' epsc0 fcu' epsu
uniaxialMaterial Concrete01 1 -3.530000E+01 -2.000000E-03 -7.100000E+00 -4.000000E-03
# Material "¸Ö½î": matTag Fy Fu E Eh epsh epsu <-GABuck lsr beta r gamma> <-DMBuck lsr <alpha>> <-CMFatigue Cf alpha Cd> <-IsoHard <a1 <limit>>> <-MPCurveParams R1 R2 R3>
uniaxialMaterial ReinforcingSteel 3 +4.550000E+02 +4.980000E+02 +1.900000E+05 +1.075000E+05 +2.900000E-03 +3.300000E-03
# Section "¶ÕÖù": secTag
section Fiber 1 {
# PatchCirc "»ìÄýÍÁ": matTag nfcirc nfrad y z intRad extRad <startAng endAng>
patch circ 1 5 5 +0.000000E+00 +0.000000E+00 +1.250000E+02 +1.500000E+02 +0.000000E+00 +3.600000E+02
# LayerCircular "¸Ö½î": matTag numBar areaBar y z radius <startAng endAng>
layer circ 3 10 +0.000000E+00 +0.000000E+00 +1.190000E+02 +0.000000E+00 +3.600000E+02
}

# Section "¸ÇÁº": secTag
section Fiber 2 {
# PatchQuad "Ô¼Êø": matTag NSIJ NSJK Iy Iz Jy Jz Ky Kz Ly Lz
patch quad 1 5 5 -1.750000E+02 -1.750000E+02 +1.750000E+02 -1.750000E+02 +1.750000E+02 +1.750000E+02 -1.750000E+02 +1.750000E+02
# PatchBox "·ÇÔ¼Êø»ìÄýÍÁ": matTag NSIJ NSJK Iy Iz Jy Jz Ky Kz Ly Lz
patch quad 1 5 5 +1.750000E+02 +2.000000E+02 +1.750000E+02 -2.000000E+02 +2.000000E+02 -2.000000E+02 +2.000000E+02 +2.000000E+02
patch quad 1 5 5 -1.750000E+02 +2.000000E+02 -1.750000E+02 +1.750000E+02 +1.750000E+02 +1.750000E+02 +1.750000E+02 +2.000000E+02
patch quad 1 5 5 -1.750000E+02 -1.750000E+02 -1.750000E+02 -2.000000E+02 +1.750000E+02 -2.000000E+02 +1.750000E+02 -1.750000E+02
patch quad 1 5 5 -2.000000E+02 +2.000000E+02 -2.000000E+02 -2.000000E+02 -1.750000E+02 -2.000000E+02 -1.750000E+02 +2.000000E+02
# LayerStraight "µÚÒ»ÁÐ": matTag numBar areaBar yStart zStart yEnd zEnd
layer straight 3 2 +2.010620E+02 -1.750000E+02 +1.750000E+02 -1.750000E+02 -1.750000E+02
# LayerStraight "µÚÆßÁÐ": matTag numBar areaBar yStart zStart yEnd zEnd
layer straight 3 2 +2.010620E+02 +1.250000E+02 +1.750000E+02 +1.250000E+02 -1.750000E+02
# LayerStraight "µÚÈýÁÐ": matTag numBar areaBar yStart zStart yEnd zEnd
layer straight 3 2 +2.010620E+02 -8.500000E+01 +1.750000E+02 -8.500000E+01 -1.750000E+02
# LayerStraight "µÚ¶þÁÐ": matTag numBar areaBar yStart zStart yEnd zEnd
layer straight 3 2 +2.010620E+02 -1.250000E+02 +1.750000E+02 -1.250000E+02 -1.750000E+02
# LayerStraight "µÚÎåÁÐ": matTag numBar areaBar yStart zStart yEnd zEnd
layer straight 3 2 +2.010620E+02 +2.500000E+01 +1.750000E+02 +2.500000E+01 -1.750000E+02
# LayerStraight "µÚ°ËÁÐ": matTag numBar areaBar yStart zStart yEnd zEnd
layer straight 3 2 +2.010620E+02 +1.750000E+02 +1.750000E+02 +1.750000E+02 -1.750000E+02
# LayerStraight "µÚÁùÁÐ": matTag numBar areaBar yStart zStart yEnd zEnd
layer straight 3 2 +2.010620E+02 +8.500000E+01 +1.750000E+02 +8.500000E+01 -1.750000E+02
# LayerStraight "µÚËÄÁÐ": matTag numBar areaBar yStart zStart yEnd zEnd
layer straight 3 2 +2.010620E+02 -2.500000E+01 +1.750000E+02 -2.500000E+01 -1.750000E+02
}
geomTransf Linear 1
element beamWithHinges 1 1 3 1 248.12 1 248.12 32720 70685.8347 397607820.2 1
element beamWithHinges 2 2 4 1 248.12 1 248.12 32720 70685.8347 397607820.2 1
#element forceBeamColumn $eleTag $iNode $jNode $numIntgrPts $secTag $transfTag
element forceBeamColumn 3 7 8 3 2 1
element forceBeamColumn 4 5 9 3 2 1
element forceBeamColumn 5 6 10 3 2 1
puts "recorder"
recorder Node -file node9.out -time -node 9 -dof 1 2 3 disp
recorder Node -file node3.out -time -node 3 -dof 1 2 3 disp
recorder Element -file ele1.out -time -ele 1 localForce
pattern Plain 1 Linear {
load 5 +0.000000E+00 -3.742810E+05 +0.000000E+00
load 6 +0.000000E+00 -3.742810E+05 +0.000000E+00
}
puts "analysis"
constraints Transformation
numberer Plain
system BandGeneral
test EnergyIncr 1.0e-6 200
algorithm Newton
integrator LoadControl 0.1
analysis Static
analyze 10
loadConst 0.0
pattern Plain 2 Linear {
load 9 +1.000000E+00 +0.000000E+00 +0.000000E+00
}
puts "analysis"
constraints Transformation
numberer Plain
system BandGeneral
set protocol {10 -10 10 -10 15 -15 15 -15 20 -20 20 -20 30 -30 30 -30 35 -35 35 -35 40 -40 40 -40 55 -55 55 -55 57.5 -57.5 57.5 -57.5 60 -60 60 -60 67.5 -67.5 67.5 -67.5 75 -75 75 -75 80 -80 80 -80 90 -90 90 -90 }
When I run the analysis,the warnning is：
Domain::update ForceBeamColumn2d::update() - section failed in setTrial
- domain failed in update
LoadControl::update - model failed to update for new dU
WARNING NewtonRaphson::solveCurrentStep() -the Integrator failed in update()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 2 with domain at load factor 0.3
OpenSees > analyze failed, returned: -3 error flag
analysis

mhscott wrote: ↑Wed Sep 16, 2020 6:25 am Regardless, use constraints Transformation as Selim suggested.

Thank you !I will try.

selimgunay wrote: ↑Tue Sep 15, 2020 10:16 pm You may need to use constraints Transformation.

Are you using equaldof, rigid diaphragm or similar multi-dof constraints?

thanks for your suggestion，actually I am using rigid links to simulate the connection between supscture and bearing.

thanks for your suggestion，actually I am using rigid links to simulate the connection between supscture and bearing.

Hello,recently I am modeling a four span continuous beam bridge,but I just start learning opensees,when I run the analysis,some errors have appeared,and here is my code of this model,please help me fix the code ,thank you!
model BasicBuilder -ndm 3 -ndf 6
# Node tag xCrd yCrd zCrd ndf
node 1 +0.000000E+00 -3.350000E+03 +0.000000E+00 -ndf 6
node 2 +0.000000E+00 +3.350000E+03 +0.000000E+00 -ndf 6
node 3 +0.000000E+00 -3.350000E+03 +6.250000E+03 -ndf 6
node 4 +0.000000E+00 +3.350000E+03 +6.250000E+03 -ndf 6
node 5 +0.000000E+00 -3.350000E+03 +1.300000E+04 -ndf 6
node 6 +0.000000E+00 +3.350000E+03 +1.300000E+04 -ndf 6
node 7 +0.000000E+00 +0.000000E+00 +1.300000E+04 -ndf 6
node 8 +0.000000E+00 -1.675000E+03 +1.300000E+04 -ndf 6
node 9 +0.000000E+00 +1.675000E+03 +1.300000E+04 -ndf 6
node 10 +0.000000E+00 -5.780000E+03 +1.300000E+04 -ndf 6
node 11 +0.000000E+00 +5.780000E+03 +1.300000E+04 -ndf 6
node 12 +0.000000E+00 +0.000000E+00 +1.304000E+04 -ndf 6
node 13 +0.000000E+00 -1.675000E+03 +1.304000E+04 -ndf 6
node 14 +0.000000E+00 +1.675000E+03 +1.304000E+04 -ndf 6
node 15 +0.000000E+00 -3.350000E+03 +1.304000E+04 -ndf 6
node 16 +0.000000E+00 +3.350000E+03 +1.304000E+04 -ndf 6
node 17 -4.000000E+04 -5.780000E+03 +1.300000E+04 -ndf 6
node 18 -4.000000E+04 +5.780000E+03 +1.300000E+04 -ndf 6
node 19 +0.000000E+00 +0.000000E+00 +1.485000E+04 -ndf 6
node 20 +0.000000E+00 -1.675000E+03 +1.485000E+04 -ndf 6
node 21 +0.000000E+00 +1.675000E+03 +1.485000E+04 -ndf 6
node 22 +0.000000E+00 -3.350000E+03 +1.485000E+04 -ndf 6
node 23 +0.000000E+00 +3.350000E+03 +1.485000E+04 -ndf 6
node 24 +0.000000E+00 -5.780000E+03 +1.485000E+04 -ndf 6
node 25 +0.000000E+00 +5.780000E+03 +1.485000E+04 -ndf 6
node 26 -2.000000E+04 -3.350000E+03 +7.900000E+03 -ndf 6
node 27 -2.000000E+04 +3.350000E+03 +7.900000E+03 -ndf 6
node 28 -2.000000E+04 -3.350000E+03 +1.300000E+04 -ndf 6
node 29 -2.000000E+04 +3.350000E+03 +1.300000E+04 -ndf 6
node 30 -2.000000E+04 +0.000000E+00 +1.300000E+04 -ndf 6
node 31 -2.000000E+04 -1.675000E+03 +1.300000E+04 -ndf 6
node 32 -2.000000E+04 +1.675000E+03 +1.300000E+04 -ndf 6
node 33 -2.000000E+04 -5.780000E+03 +1.300000E+04 -ndf 6
node 34 -2.000000E+04 +5.780000E+03 +1.300000E+04 -ndf 6
node 35 -2.000000E+04 +0.000000E+00 +1.304000E+04 -ndf 6
node 36 -2.000000E+04 -1.675000E+03 +1.304000E+04 -ndf 6
node 37 -2.000000E+04 +1.675000E+03 +1.304000E+04 -ndf 6
node 38 -2.000000E+04 -3.350000E+03 +1.304000E+04 -ndf 6
node 39 -2.000000E+04 +3.350000E+03 +1.304000E+04 -ndf 6
node 40 -2.000000E+04 -5.780000E+03 +1.300000E+04 -ndf 6
node 41 -2.000000E+04 +5.780000E+03 +1.300000E+04 -ndf 6
node 42 -2.000000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 43 -2.000000E+04 -1.675000E+03 +1.485000E+04 -ndf 6
node 44 -2.000000E+04 +1.675000E+03 +1.485000E+04 -ndf 6
node 45 -2.000000E+04 -3.350000E+03 +1.485000E+04 -ndf 6
node 46 -2.000000E+04 +3.350000E+03 +1.485000E+04 -ndf 6
node 47 -2.000000E+04 +5.780000E+03 +1.485000E+04 -ndf 6
node 48 -2.000000E+04 -5.780000E+03 +1.485000E+04 -ndf 6
node 49 +2.000000E+04 -3.350000E+03 +5.000000E+03 -ndf 6
node 50 +2.000000E+04 +3.350000E+03 +5.000000E+03 -ndf 6
node 51 +2.000000E+04 -3.350000E+03 +1.300000E+04 -ndf 6
node 52 +2.000000E+04 +3.350000E+03 +1.300000E+04 -ndf 6
node 53 +2.000000E+04 +0.000000E+00 +1.300000E+04 -ndf 6
node 54 +2.000000E+04 -1.675000E+03 +1.300000E+04 -ndf 6
node 55 +2.000000E+04 +1.675000E+03 +1.300000E+04 -ndf 6
node 56 +2.000000E+04 -5.780000E+03 +1.300000E+04 -ndf 6
node 57 +2.000000E+04 +5.780000E+03 +1.300000E+04 -ndf 6
node 58 +2.000000E+04 +0.000000E+00 +1.304000E+04 -ndf 6
node 59 +2.000000E+04 -1.675000E+03 +1.304000E+04 -ndf 6
node 60 +2.000000E+04 +1.675000E+03 +1.304000E+04 -ndf 6
node 61 +2.000000E+04 -3.350000E+03 +1.304000E+04 -ndf 6
node 62 +2.000000E+04 +3.350000E+03 +1.304000E+04 -ndf 6
node 63 +0.000000E+00 -5.780000E+03 +1.300000E+04 -ndf 6
node 64 +0.000000E+00 +5.780000E+03 +1.300000E+04 -ndf 6
node 65 +2.000000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 66 +2.000000E+04 -1.675000E+03 +1.485000E+04 -ndf 6
node 67 +2.000000E+04 +1.675000E+03 +1.485000E+04 -ndf 6
node 68 +2.000000E+04 -3.350000E+03 +1.485000E+04 -ndf 6
node 69 +2.000000E+04 +3.350000E+03 +1.485000E+04 -ndf 6
node 70 +2.000000E+04 -5.780000E+03 +1.485000E+04 -ndf 6
node 71 +2.000000E+04 +5.780000E+03 +1.485000E+04 -ndf 6
node 72 -4.000000E+04 +0.000000E+00 +1.300000E+04 -ndf 6
node 73 -4.000000E+04 -1.675000E+03 +1.300000E+04 -ndf 6
node 74 -4.000000E+04 +1.675000E+03 +1.300000E+04 -ndf 6
node 75 -4.000000E+04 -3.350000E+03 +1.300000E+04 -ndf 6
node 76 -4.000000E+04 +3.350000E+03 +1.300000E+04 -ndf 6
node 77 -4.000000E+04 -5.780000E+03 +1.300000E+04 -ndf 6
node 78 -4.000000E+04 +5.780000E+03 +1.300000E+04 -ndf 6
node 79 -4.000000E+04 +0.000000E+00 +1.304000E+04 -ndf 6
node 80 -4.000000E+04 -1.675000E+03 +1.304000E+04 -ndf 6
node 81 -4.000000E+04 +1.675000E+03 +1.304000E+04 -ndf 6
node 82 -4.000000E+04 -3.350000E+03 +1.304000E+04 -ndf 6
node 83 -4.000000E+04 +3.350000E+03 +1.304000E+04 -ndf 6
node 84 +2.000000E+04 -5.780000E+03 +1.300000E+04 -ndf 6
node 85 +2.000000E+04 +5.780000E+03 +1.300000E+04 -ndf 6
node 86 -4.000000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 87 -4.000000E+04 -1.675000E+03 +1.485000E+04 -ndf 6
node 88 -4.000000E+04 +1.675000E+03 +1.485000E+04 -ndf 6
node 89 -4.000000E+04 -3.350000E+03 +1.485000E+04 -ndf 6
node 90 -4.000000E+04 +3.350000E+03 +1.485000E+04 -ndf 6
node 91 -4.000000E+04 -5.780000E+03 +1.485000E+04 -ndf 6
node 92 -4.000000E+04 +5.780000E+03 +1.485000E+04 -ndf 6
node 93 +4.000000E+04 +0.000000E+00 +1.300000E+04 -ndf 6
node 94 +4.000000E+04 -1.675000E+03 +1.300000E+04 -ndf 6
node 95 +4.000000E+04 +1.675000E+03 +1.300000E+04 -ndf 6
node 96 +4.000000E+04 -3.350000E+03 +1.300000E+04 -ndf 6
node 97 +4.000000E+04 +3.350000E+03 +1.300000E+04 -ndf 6
node 98 +4.000000E+04 -5.780000E+03 +1.300000E+04 -ndf 6
node 99 +4.000000E+04 +5.780000E+03 +1.300000E+04 -ndf 6
node 100 +4.000000E+04 +0.000000E+00 +1.304000E+04 -ndf 6
node 101 +4.000000E+04 -1.675000E+03 +1.304000E+04 -ndf 6
node 102 +4.000000E+04 +1.675000E+03 +1.304000E+04 -ndf 6
node 103 +4.000000E+04 -3.350000E+03 +1.304000E+04 -ndf 6
node 104 +4.000000E+04 +3.350000E+03 +1.304000E+04 -ndf 6
node 105 +4.000000E+04 -5.780000E+03 +1.300000E+04 -ndf 6
node 106 +4.000000E+04 +5.780000E+03 +1.300000E+04 -ndf 6
node 107 +4.000000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 108 +4.000000E+04 -1.675000E+03 +1.485000E+04 -ndf 6
node 109 +4.000000E+04 +1.675000E+03 +1.485000E+04 -ndf 6
node 110 +4.000000E+04 -3.350000E+03 +1.485000E+04 -ndf 6
node 111 +4.000000E+04 +3.350000E+03 +1.485000E+04 -ndf 6
node 112 +4.000000E+04 -5.780000E+03 +1.485000E+04 -ndf 6
node 113 +4.000000E+04 +5.780000E+03 +1.485000E+04 -ndf 6
node 114 -4.000000E+04 +0.000000E+00 +1.200000E+04 -ndf 6
node 115 +4.000000E+04 +0.000000E+00 +1.200000E+04 -ndf 6
node 116 -3.951000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 117 -3.810000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 118 -3.495000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 119 -3.334000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 120 -3.316000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 121 -2.684000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 122 -2.666000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 123 -2.505000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 124 -2.190000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 125 -2.049000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 126 -1.860000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 127 -1.770000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 128 -1.455000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 129 -1.334000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 130 -1.316000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 131 -6.840000E+03 +0.000000E+00 +1.485000E+04 -ndf 6
node 132 -6.660000E+03 +0.000000E+00 +1.485000E+04 -ndf 6
node 133 -5.450000E+03 +0.000000E+00 +1.485000E+04 -ndf 6
node 134 -2.300000E+03 +0.000000E+00 +1.485000E+04 -ndf 6
node 135 -1.400000E+03 +0.000000E+00 +1.485000E+04 -ndf 6
node 136 +3.951000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 137 +3.810000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 138 +3.495000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 139 +3.334000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 140 +3.316000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 141 +2.684000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 142 +2.666000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 143 +2.505000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 144 +2.190000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 145 +2.049000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 146 +1.860000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 147 +1.770000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 148 +1.455000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 149 +1.334000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 150 +1.316000E+04 +0.000000E+00 +1.485000E+04 -ndf 6
node 151 +6.840000E+03 +0.000000E+00 +1.485000E+04 -ndf 6
node 152 +6.660000E+03 +0.000000E+00 +1.485000E+04 -ndf 6
node 153 +5.450000E+03 +0.000000E+00 +1.485000E+04 -ndf 6
node 154 +2.300000E+03 +0.000000E+00 +1.485000E+04 -ndf 6
node 155 +1.400000E+03 +0.000000E+00 +1.485000E+04 -ndf 6
# SPC tag Dx Dy Dz Rx Ry Rz
fix 1 1 1 1 1 1 1
fix 2 1 1 1 1 1 1
fix 26 1 1 1 1 1 1
fix 27 1 1 1 1 1 1
fix 49 1 1 1 1 1 1
fix 50 1 1 1 1 1 1
fix 114 1 1 1 1 1 1
fix 115 1 1 1 1 1 1
mass 19 +1.235976E+04 +1.235976E+04 +1.235976E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 42 +1.015650E+04 +1.015650E+04 +1.015650E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 65 +1.015650E+04 +1.015650E+04 +1.015650E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 86 +2.633167E+03 +2.633167E+03 +2.633167E+03 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 107 +2.633167E+03 +2.633167E+03 +2.633167E+03 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 116 +1.021024E+04 +1.021024E+04 +1.021024E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 117 +2.450457E+04 +2.450457E+04 +2.450457E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 118 +2.557934E+04 +2.557934E+04 +2.557934E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 119 +9.619120E+03 +9.619120E+03 +9.619120E+03 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 120 +3.492976E+04 +3.492976E+04 +3.492976E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 121 +3.492976E+04 +3.492976E+04 +3.492976E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 122 +9.619120E+03 +9.619120E+03 +9.619120E+03 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 123 +2.557934E+04 +2.557934E+04 +2.557934E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 124 +2.450457E+04 +2.450457E+04 +2.450457E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 125 +1.021024E+04 +1.021024E+04 +1.021024E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 126 +1.235976E+04 +1.235976E+04 +1.235976E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 127 +2.176393E+04 +2.176393E+04 +2.176393E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 128 +2.342981E+04 +2.342981E+04 +2.342981E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 129 +7.469596E+03 +7.469596E+03 +7.469596E+03 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 130 +3.492976E+04 +3.492976E+04 +3.492976E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 131 +3.492976E+04 +3.492976E+04 +3.492976E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 132 +7.469596E+03 +7.469596E+03 +7.469596E+03 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 133 +2.342981E+04 +2.342981E+04 +2.342981E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 134 +2.176393E+04 +2.176393E+04 +2.176393E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 135 +1.235976E+04 +1.235976E+04 +1.235976E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 136 +1.021024E+04 +1.021024E+04 +1.021024E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 137 +2.450457E+04 +2.450457E+04 +2.450457E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 138 +2.557934E+04 +2.557934E+04 +2.557934E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 139 +9.619120E+03 +9.619120E+03 +9.619120E+03 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 140 +3.492976E+04 +3.492976E+04 +3.492976E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 141 +3.492976E+04 +3.492976E+04 +3.492976E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 142 +9.619112E+03 +9.619120E+03 +9.619120E+03 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 143 +2.557934E+04 +2.557934E+04 +2.557934E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 144 +2.450457E+04 +2.450457E+04 +2.450457E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 145 +1.021024E+04 +1.021024E+04 +1.021024E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 146 +1.235976E+04 +1.235976E+04 +1.235976E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 147 +2.176393E+04 +2.176393E+04 +2.176393E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 148 +2.342981E+04 +2.342981E+04 +2.342981E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 149 +7.469596E+03 +7.469596E+03 +7.469596E+03 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 150 +3.492976E+04 +3.492976E+04 +3.492976E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 151 +3.392976E+04 +3.392976E+04 +3.392976E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 152 +7.469596E+03 +7.469596E+03 +7.469596E+03 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 153 +2.342981E+04 +2.342981E+04 +2.342981E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 154 +2.176393E+04 +2.176393E+04 +2.176393E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
mass 155 +1.235976E+04 +1.235976E+04 +1.235976E+04 +0.000000E+00 +0.000000E+00 +0.000000E+00
equalDOF 19 12 1 2 3 4 5 6
equalDOF 19 20 1 2 3 4 5 6
equalDOF 20 13 1 2 3 4 5 6
equalDOF 20 22 1 2 3 4 5 6
equalDOF 21 14 1 2 3 4 5 6
equalDOF 21 19 1 2 3 4 5 6
equalDOF 22 15 1 2 3 4 5 6
equalDOF 22 24 1 2 3 4 5 6
equalDOF 23 16 1 2 3 4 5 6
equalDOF 23 21 1 2 3 4 5 6
equalDOF 24 63 1 2 3 4 5 6
equalDOF 25 23 1 2 3 4 5 6
equalDOF 41 47 1 2 3 4 5 6
equalDOF 42 35 1 2 3 4 5 6
equalDOF 42 43 1 2 3 4 5 6
equalDOF 43 36 1 2 3 4 5 6
equalDOF 43 45 1 2 3 4 5 6
equalDOF 44 37 1 2 3 4 5 6
equalDOF 44 42 1 2 3 4 5 6
equalDOF 45 38 1 2 3 4 5 6
equalDOF 45 48 1 2 3 4 5 6
equalDOF 46 39 1 2 3 4 5 6
equalDOF 46 44 1 2 3 4 5 6
equalDOF 47 46 1 2 3 4 5 6
equalDOF 48 40 1 2 3 4 5 6
equalDOF 64 25 1 2 3 4 5 6
equalDOF 65 58 1 2 3 4 5 6
equalDOF 65 66 1 2 3 4 5 6
equalDOF 66 59 1 2 3 4 5 6
equalDOF 66 68 1 2 3 4 5 6
equalDOF 67 60 1 2 3 4 5 6
equalDOF 67 65 1 2 3 4 5 6
equalDOF 68 61 1 2 3 4 5 6
equalDOF 68 70 1 2 3 4 5 6
equalDOF 69 62 1 2 3 4 5 6
equalDOF 69 67 1 2 3 4 5 6
equalDOF 70 84 1 2 3 4 5 6
equalDOF 71 69 1 2 3 4 5 6
equalDOF 78 92 1 2 3 4 5 6
equalDOF 85 71 1 2 3 4 5 6
equalDOF 86 79 1 2 3 4 5 6
equalDOF 86 87 1 2 3 4 5 6
equalDOF 87 80 1 2 3 4 5 6
equalDOF 87 89 1 2 3 4 5 6
equalDOF 88 81 1 2 3 4 5 6
equalDOF 88 86 1 2 3 4 5 6
equalDOF 89 82 1 2 3 4 5 6
equalDOF 89 91 1 2 3 4 5 6
equalDOF 90 83 1 2 3 4 5 6
equalDOF 90 88 1 2 3 4 5 6
equalDOF 91 77 1 2 3 4 5 6
equalDOF 92 90 1 2 3 4 5 6
equalDOF 106 113 1 2 3 4 5 6
equalDOF 107 100 1 2 3 4 5 6
equalDOF 108 101 1 2 3 4 5 6
equalDOF 108 110 1 2 3 4 5 6
equalDOF 109 102 1 2 3 4 5 6
equalDOF 109 107 1 2 3 4 5 6
rigidLink beam 108 110
rigidLink beam 110 103
rigidLink beam 110 112
rigidLink beam 112 105
rigidLink beam 114 72
rigidLink beam 115 93
uniaxialMaterial Elastic 1 +2.900000E+04 +0.000000E+00

uniaxialMaterial Concrete01 2 -2.200000E+01 -4.000000E-03 -3.000000E+01 -1.400000E-02

uniaxialMaterial Concrete01 3 -2.200000E+01 -4.000000E-03 -3.000000E+01 -1.400000E-02

uniaxialMaterial Steel02 4 +3.850000E+02 +2.000000E+05 +5.000000E-02 +1.850000E+01 +9.250000E-01 +1.500000E-01 +0.000000E+00 +1.000000E+00 +0.000000E+00 +1.000000E+00 +0.000000E+00

uniaxialMaterial BoucWen 5 +6.250000E-02 +2.900000E+04 +1.000000E+00 +5.000000E-01 +5.000000E-01 +1.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +1.000000E-08 20

frictionModel Coulomb 1 +5.000000E-02

frictionModel Coulomb 2 +5.000000E-02

section Elastic 1 +2.900000E+04 +1.800000E+02 +4.860000E+03 +1.500000E+03 +1.115400E+04 +3.916000E+03 +8.333333E-01 +8.333333E-01


section Fiber 2 -GJ +5.000000E+00 {
# PatchCirc "concrete": matTag nfcirc nfrad y z intRad extRad <startAng endAng>
patch circ 3 5 5 +0.000000E+00 +0.000000E+00 +1.240000E+03 +1.300000E+03 +0.000000E+00 +3.600000E+02
# LayerCircular "steel": matTag numBar areaBar y z radius <startAng endAng>
layer circ 4 20 +4.908738E+02 +0.000000E+00 +0.000000E+00 +1.250000E+01
}

# Section "contactbeam": secTag
section Fiber 3 -GJ +5.000000E+00 {
# PatchQuad "contactbeam": matTag NSIJ NSJK Iy Iz Jy Jz Ky Kz Ly Lz
patch quad 3 5 5 -5.000000E+02 -6.000000E+02 +5.000000E+02 -6.000000E+02 +5.000000E+02 +6.000000E+02 -5.000000E+02 +6.000000E+02
# LayerStraight "Layer01": matTag numBar areaBar yStart zStart yEnd zEnd
layer straight 4 36 +2.010620E+02 +0.000000E+00 -5.500000E+02 -4.500000E+02 -5.500000E+02
}

geomTransf Linear 1 +1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 2 +0.000000E+00 +1.000000E+00 +0.000000E+00

geomTransf Linear 3 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 4 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 5 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 6 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 7 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 8 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 9 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 10 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 11 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 12 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 13 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 14 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 15 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 16 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 17 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 18 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 19 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 20 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 21 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 22 +0.000000E+00 -1.000000E+00 +0.000000E+00

geomTransf Linear 23 +0.000000E+00 -1.000000E+00 +0.000000E+00

geomTransf Linear 24 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 25 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 26 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 27 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 28 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 29 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 30 +0.000000E+00 -1.000000E+00 +0.000000E+00

geomTransf Linear 31 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 32 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 33 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 34 -1.000000E+00 +0.000000E+00 +0.000000E+00

geomTransf Linear 35 -1.000000E+00 +0.000000E+00 +0.000000E+00

element elasticBeamColumn 1 107 136 +1.405200E+07 +3.000000E+04 +1.200000E+04 +6.688310E+12 +2.025486E+12 +1.416863E+14 2 -mass +1.074762E+01


element elasticBeamColumn 2 136 137 +5.068700E+06 +3.000000E+04 +1.200000E+04 +2.213240E+11 +9.598532E+11 +6.249725E+13 2 -mass +1.074762E+01


element elasticBeamColumn 3 137 138 +4.723150E+06 +3.000000E+04 +1.200000E+04 +1.353550E+11 +1.107110E+12 +5.826396E+13 2 -mass +1.074762E+01


element elasticBeamColumn 4 138 139 +3.993000E+06 +3.000000E+04 +1.200000E+04 +7.538207E+10 +8.375210E+11 +4.922611E+13 2 -mass +1.074762E+01


element elasticBeamColumn 5 139 140 +1.311300E+07 +3.000000E+04 +1.200000E+04 +5.660703E+12 +1.795933E+12 +1.243186E+14 2 -mass +1.074762E+01


element elasticBeamColumn 6 140 141 +3.993000E+06 +3.000000E+04 +1.200000E+04 +7.538207E+10 +8.375210E+11 +4.922611E+13 2 -mass +1.074762E+01


element elasticBeamColumn 7 141 142 +1.311300E+07 +3.000000E+04 +1.200000E+04 +5.660703E+12 +1.795933E+12 +1.243186E+14 2 -mass +1.074762E+01


element elasticBeamColumn 8 142 143 +3.993000E+06 +3.000000E+04 +1.200000E+04 +7.538207E+10 +8.375210E+11 +4.922611E+13 2 -mass +1.074762E+01


element elasticBeamColumn 9 143 144 +4.723150E+06 +3.000000E+04 +1.200000E+04 +1.353550E+11 +1.107110E+12 +5.826396E+13 2 -mass +1.074762E+01


element elasticBeamColumn 10 144 145 +5.068700E+06 +3.000000E+04 +1.200000E+04 +2.213240E+11 +9.598532E+11 +6.249725E+13 2 -mass +1.074762E+01


element elasticBeamColumn 11 145 65 +1.234520E+07 +3.000000E+04 +1.200000E+04 +7.170212E+12 +2.070309E+12 +8.071427E+13 2 -mass +1.074762E+01


element elasticBeamColumn 12 65 146 +1.234520E+07 +3.000000E+04 +1.200000E+04 +7.170212E+12 +2.070309E+12 +8.071427E+13 2 -mass +1.074762E+01


element elasticBeamColumn 13 146 147 +5.068700E+06 +3.000000E+04 +1.200000E+04 +2.213240E+11 +9.598532E+11 +6.249725E+13 2 -mass +1.074762E+01


element elasticBeamColumn 14 147 148 +4.723150E+06 +3.000000E+04 +1.200000E+04 +1.353550E+11 +1.107110E+12 +5.826396E+13 2 -mass +1.074762E+01


element elasticBeamColumn 15 148 149 +3.993000E+06 +3.000000E+04 +1.200000E+04 +7.538207E+10 +8.375210E+11 +4.922611E+13 2 -mass +1.074762E+01


element elasticBeamColumn 16 149 150 +1.311300E+07 +3.000000E+04 +1.200000E+04 +5.660703E+12 +1.795933E+12 +1.243186E+14 2 -mass +1.074762E+01


element elasticBeamColumn 17 150 151 +3.993000E+06 +3.000000E+04 +1.200000E+04 +7.538207E+10 +8.375210E+11 +4.922611E+13 2 -mass +1.074762E+01


element elasticBeamColumn 18 151 152 +1.311300E+07 +3.000000E+04 +1.200000E+04 +5.660703E+12 +1.795933E+12 +1.243186E+14 2 -mass +1.074762E+01


element elasticBeamColumn 19 152 153 +3.993000E+06 +3.000000E+04 +1.200000E+04 +7.538207E+10 +8.375210E+11 +4.922611E+13 2 -mass +1.074762E+01


element elasticBeamColumn 20 153 154 +4.723150E+06 +3.000000E+04 +1.200000E+04 +1.353550E+11 +1.107110E+12 +5.826396E+13 2 -mass +1.074762E+01


element elasticBeamColumn 21 154 155 +5.068700E+06 +3.000000E+04 +1.200000E+04 +2.213240E+11 +9.598532E+11 +6.249725E+13 2 -mass +1.074762E+01


element elasticBeamColumn 22 155 19 +1.234520E+07 +3.000000E+04 +1.200000E+04 +7.170212E+12 +2.070309E+12 +8.071427E+13 2 -mass +1.074762E+01


element elasticBeamColumn 23 19 135 +1.234520E+07 +3.000000E+04 +1.200000E+04 +7.170212E+12 +2.070309E+12 +8.071427E+13 2 -mass +1.074762E+01


element elasticBeamColumn 24 135 134 +5.068700E+06 +3.000000E+04 +1.200000E+04 +2.213240E+11 +9.598532E+11 +6.249725E+13 2 -mass +1.074762E+01


element elasticBeamColumn 25 134 133 +4.723150E+06 +3.000000E+04 +1.200000E+04 +1.353550E+11 +1.107110E+12 +5.826396E+13 2 -mass +1.074762E+01


element elasticBeamColumn 26 133 132 +3.993000E+06 +3.000000E+04 +1.200000E+04 +7.538207E+10 +8.375210E+11 +4.922611E+13 2 -mass +1.074762E+01


element elasticBeamColumn 27 132 131 +1.234520E+07 +3.000000E+04 +1.200000E+04 +7.170212E+12 +2.070309E+12 +8.071427E+13 2 -mass +1.074762E+01


element elasticBeamColumn 28 131 130 +1.311300E+07 +3.000000E+04 +1.200000E+04 +5.660703E+12 +1.795933E+12 +1.243186E+14 2 -mass +1.074762E+01


element elasticBeamColumn 29 130 129 +1.311300E+07 +3.000000E+04 +1.200000E+04 +5.660703E+12 +1.795933E+12 +1.243186E+14 2 -mass +1.074762E+01


element elasticBeamColumn 30 129 128 +3.993000E+06 +3.000000E+04 +1.200000E+04 +7.538207E+10 +8.375210E+11 +4.922611E+13 2 -mass +1.074762E+01


element elasticBeamColumn 31 128 127 +4.723150E+06 +3.000000E+04 +1.200000E+04 +1.353550E+11 +1.107110E+12 +5.826396E+13 2 -mass +1.074762E+01

element elasticBeamColumn 32 127 126 +5.068700E+06 +3.000000E+04 +1.200000E+04 +2.213240E+11 +9.598532E+11 +6.249725E+13 2 -mass +1.074762E+01


element elasticBeamColumn 33 126 42 +1.234520E+07 +3.000000E+04 +1.200000E+04 +7.170212E+12 +2.070309E+12 +8.071427E+13 2 -mass +1.074762E+01


element elasticBeamColumn 34 42 125 +1.234520E+07 +3.000000E+04 +1.200000E+04 +7.170212E+12 +2.070309E+12 +8.071427E+13 2 -mass +1.074762E+01


element elasticBeamColumn 35 125 124 +5.068700E+06 +3.000000E+04 +1.200000E+04 +2.213240E+11 +9.598532E+11 +6.249725E+13 2 -mass +1.074762E+01


element elasticBeamColumn 36 124 123 +4.723150E+06 +3.000000E+04 +1.200000E+04 +1.353550E+11 +1.107110E+12 +5.826396E+13 2 -mass +1.074762E+01


element elasticBeamColumn 37 123 122 +3.993000E+06 +3.000000E+04 +1.200000E+04 +7.538207E+10 +8.375210E+11 +4.922611E+13 2 -mass +1.074762E+01


element elasticBeamColumn 38 122 121 +1.311300E+07 +3.000000E+04 +1.200000E+04 +5.660703E+12 +1.795933E+12 +1.243186E+14 2 -mass +1.074762E+01


element elasticBeamColumn 39 121 120 +1.311300E+07 +3.000000E+04 +1.200000E+04 +5.660703E+12 +1.795933E+12 +1.243186E+14 2 -mass +1.074762E+01


element elasticBeamColumn 40 120 119 +1.311300E+07 +3.000000E+04 +1.200000E+04 +5.660703E+12 +1.795933E+12 +1.243186E+14 2 -mass +1.074762E+01

elasticBeamColumn 41 119 118 +3.993000E+06 +3.000000E+04 +1.200000E+04 +7.538207E+10 +8.375210E+11 +4.922611E+13 2 -mass +1.074762E+01


element elasticBeamColumn 42 118 117 +4.723150E+06 +3.000000E+04 +1.200000E+04 +1.353550E+11 +1.107110E+12 +5.826396E+13 2 -mass +1.074762E+01


element elasticBeamColumn 43 117 116 +5.068700E+06 +3.000000E+04 +1.200000E+04 +2.213240E+11 +9.598532E+11 +6.249725E+13 2 -mass +1.074762E+01


element elasticBeamColumn 44 116 86 +1.405200E+07 +3.000000E+04 +1.200000E+04 +6.688310E+12 +2.025486E+12 +1.416863E+14 2 -mass +1.074762E+01


element dispBeamColumn 45 3 4 3 3 1 -mass +0.000000E+00 -integration Legendre


element elasticBeamColumn 46 93 94 +2.340000E+06 +3.000000E+04 +1.200000E+04 +7.410189E+11 +3.295500E+11 +6.318000E+11 3 -mass +0.000000E+00


element elasticBeamColumn 47 94 96 +2.340000E+06 +3.000000E+04 +1.200000E+04 +7.410189E+11 +3.295500E+11 +6.318000E+11 4 -mass +0.000000E+00


element elasticBeamColumn 48 96 98 +2.340000E+06 +3.000000E+04 +1.200000E+04 +7.410189E+11 +3.295500E+11 +6.318000E+11 5 -mass +0.000000E+00

# Element "damper": eleTag NodeI NodeJ -mat matTag1 matTag2 ... -dir dir1 dir2 ... <-orient x1 x2 x3 y1 y2 y3> <-doRayleigh flag>
element zeroLength 49 78 18 -mat 1 1 1 5 5 -dir 1 2 3 5 6 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +1.000000E+00 -doRayleigh 1

# Element "damper": eleTag NodeI NodeJ -mat matTag1 matTag2 ... -dir dir1 dir2 ... <-orient x1 x2 x3 y1 y2 y3> <-doRayleigh flag>
element zeroLength 50 77 17 -mat 1 1 1 5 5 -dir 1 2 3 5 6 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +1.000000E+00 -doRayleigh 1

# Element "capbeam2": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 51 99 97 +2.340000E+06 +3.000000E+04 +1.200000E+04 +7.410189E+11 +3.295500E+11 +6.318000E+11 6 -mass +0.000000E+00

# Element "capbeam2": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 52 75 77 +2.340000E+06 +3.000000E+04 +1.200000E+04 +7.410189E+11 +3.295500E+11 +6.318000E+11 7 -mass +0.000000E+00

# Element "damper": eleTag NodeI NodeJ -mat matTag1 matTag2 ... -dir dir1 dir2 ... <-orient x1 x2 x3 y1 y2 y3> <-doRayleigh flag>
element zeroLength 53 34 41 -mat 1 1 1 5 5 -dir 1 2 3 5 6 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +1.000000E+00 -doRayleigh 1

# Element "damper": eleTag NodeI NodeJ -mat matTag1 matTag2 ... -dir dir1 dir2 ... <-orient x1 x2 x3 y1 y2 y3> <-doRayleigh flag>
element zeroLength 54 33 40 -mat 1 1 1 5 5 -dir 1 2 3 5 6 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +1.000000E+00 -doRayleigh 1

# Element "damper": eleTag NodeI NodeJ -mat matTag1 matTag2 ... -dir dir1 dir2 ... <-orient x1 x2 x3 y1 y2 y3> <-doRayleigh flag>
element zeroLength 55 11 64 -mat 1 1 1 5 5 -dir 1 2 3 5 6 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +1.000000E+00 -doRayleigh 1

# Element "damper": eleTag NodeI NodeJ -mat matTag1 matTag2 ... -dir dir1 dir2 ... <-orient x1 x2 x3 y1 y2 y3> <-doRayleigh flag>
element zeroLength 56 10 63 -mat 1 1 1 5 5 -dir 1 2 3 5 6 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +1.000000E+00 -doRayleigh 1

# Element "damper": eleTag NodeI NodeJ -mat matTag1 matTag2 ... -dir dir1 dir2 ... <-orient x1 x2 x3 y1 y2 y3> <-doRayleigh flag>
element zeroLength 57 57 85 -mat 1 1 1 5 5 -dir 1 2 3 5 6 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +1.000000E+00 -doRayleigh 1

# Element "damper": eleTag NodeI NodeJ -mat matTag1 matTag2 ... -dir dir1 dir2 ... <-orient x1 x2 x3 y1 y2 y3> <-doRayleigh flag>
element zeroLength 58 56 84 -mat 1 1 1 5 5 -dir 1 2 3 5 6 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +1.000000E+00 -doRayleigh 1

# Element "capbeam2": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 59 97 95 +2.340000E+06 +3.000000E+04 +1.200000E+04 +7.410189E+11 +3.295500E+11 +6.318000E+11 8 -mass +0.000000E+00

# Element "capbeam2": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 60 95 93 +2.340000E+06 +3.000000E+04 +1.200000E+04 +7.410189E+11 +3.295500E+11 +6.318000E+11 9 -mass +0.000000E+00

# Element "bearing2": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 61 97 104 2 +1.658810E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "damper": eleTag NodeI NodeJ -mat matTag1 matTag2 ... -dir dir1 dir2 ... <-orient x1 x2 x3 y1 y2 y3> <-doRayleigh flag>
element zeroLength 62 99 106 -mat 1 1 1 5 5 -dir 1 2 3 5 6 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +1.000000E+00 -doRayleigh 1

# Element "bearing2": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 63 95 102 2 +1.658810E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "damper": eleTag NodeI NodeJ -mat matTag1 matTag2 ... -dir dir1 dir2 ... <-orient x1 x2 x3 y1 y2 y3> <-doRayleigh flag>
element zeroLength 64 98 105 -mat 1 1 1 5 5 -dir 1 2 3 5 6 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +0.000000E+00 +1.000000E+00 -doRayleigh 1

# Element "bearing2": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 65 93 100 2 +1.658810E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "bearing2": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 67 94 101 2 +1.658810E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "bearing2": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 69 96 103 2 +1.658810E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "capbeam": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 71 53 54 +2.700000E+06 +3.000000E+04 +1.200000E+04 +1.018119E+12 +5.062500E+11 +7.290000E+11 10 -mass +0.000000E+00

# Element "capbeam": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 72 54 51 +2.700000E+06 +3.000000E+04 +1.200000E+04 +1.018119E+12 +5.062500E+11 +7.290000E+11 11 -mass +0.000000E+00

# Element "capbeam": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 73 51 56 +2.700000E+06 +3.000000E+04 +1.200000E+04 +1.018119E+12 +5.062500E+11 +7.290000E+11 12 -mass +0.000000E+00

# Element "capbeam": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 84 57 52 +2.700000E+06 +3.000000E+04 +1.200000E+04 +1.018119E+12 +5.062500E+11 +7.290000E+11 13 -mass +0.000000E+00

# Element "capbeam": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 85 52 55 +2.700000E+06 +3.000000E+04 +1.200000E+04 +1.018119E+12 +5.062500E+11 +7.290000E+11 14 -mass +0.000000E+00

# Element "capbeam": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 86 55 53 +2.700000E+06 +3.000000E+04 +1.200000E+04 +1.018119E+12 +5.062500E+11 +7.290000E+11 15 -mass +0.000000E+00

# Element "bearing1": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 87 52 62 1 +2.409744E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "ElasticDefault": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 89 55 60 +1.800000E+02 +2.900000E+04 +1.115400E+04 +3.916000E+03 +1.500000E+03 +4.860000E+03 2 -mass +0.000000E+00

# Element "bearing1": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 91 53 58 1 +2.409744E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "bearing1": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 93 54 59 1 +2.409744E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "bearing1": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 95 51 61 1 +2.409744E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "column": eleTag NodeI NodeJ NIP secTag geoTranTag <-mass massDens> <-cMass>
element dispBeamColumn 97 49 51 5 2 2 -mass +0.000000E+00 -integration Legendre

# Element "column": eleTag NodeI NodeJ NIP secTag geoTranTag <-mass massDens> <-cMass>
element dispBeamColumn 98 50 52 5 2 2 -mass +0.000000E+00 -integration Legendre

# Element "capbeam": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 99 7 8 +2.700000E+06 +3.000000E+04 +1.200000E+04 +1.018119E+12 +5.062500E+11 +7.290000E+11 16 -mass +0.000000E+00

# Element "capbeam": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 100 8 5 +2.700000E+06 +3.000000E+04 +1.200000E+04 +1.018119E+12 +5.062500E+11 +7.290000E+11 17 -mass +0.000000E+00

# Element "capbeam": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 101 5 10 +2.700000E+06 +3.000000E+04 +1.200000E+04 +1.018119E+12 +5.062500E+11 +7.290000E+11 18 -mass +0.000000E+00

# Element "capbeam": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 112 11 6 +2.700000E+06 +3.000000E+04 +1.200000E+04 +1.018119E+12 +5.062500E+11 +7.290000E+11 19 -mass +0.000000E+00

# Element "capbeam": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 113 6 9 +2.700000E+06 +3.000000E+04 +1.200000E+04 +1.018119E+12 +5.062500E+11 +7.290000E+11 20 -mass +0.000000E+00

# Element "capbeam": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 114 9 7 +2.700000E+06 +3.000000E+04 +1.200000E+04 +1.018119E+12 +5.062500E+11 +7.290000E+11 21 -mass +0.000000E+00

# Element "bearing1": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 116 16 6 1 +2.409744E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "bearing1": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 118 14 9 1 +2.409744E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "bearing1": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 120 12 7 1 +2.409744E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "bearing1": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 122 13 8 1 +2.409744E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "bearing1": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 124 15 5 1 +2.409744E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "column": eleTag NodeI NodeJ NIP secTag geoTranTag <-mass massDens> <-cMass>
element dispBeamColumn 125 1 3 5 2 2 -mass +0.000000E+00 -integration Legendre

# Element "column": eleTag NodeI NodeJ NIP secTag geoTranTag <-mass massDens> <-cMass>
element dispBeamColumn 126 3 5 5 2 2 -mass +0.000000E+00 -integration Legendre

# Element "column": eleTag NodeI NodeJ NIP secTag geoTranTag <-mass massDens> <-cMass>
element dispBeamColumn 127 6 4 5 2 22 -mass +0.000000E+00 -integration Legendre

# Element "column": eleTag NodeI NodeJ NIP secTag geoTranTag <-mass massDens> <-cMass>
element dispBeamColumn 128 4 2 5 2 23 -mass +0.000000E+00 -integration Legendre

# Element "capbeam": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 129 30 31 +2.700000E+06 +3.000000E+04 +1.200000E+04 +1.018119E+12 +5.062500E+11 +7.290000E+11 24 -mass +0.000000E+00

# Element "capbeam": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 130 31 28 +2.700000E+06 +3.000000E+04 +1.200000E+04 +1.018119E+12 +5.062500E+11 +7.290000E+11 25 -mass +0.000000E+00

# Element "capbeam": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 131 28 33 +2.700000E+06 +3.000000E+04 +1.200000E+04 +1.018119E+12 +5.062500E+11 +7.290000E+11 26 -mass +0.000000E+00

# Element "capbeam": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 142 34 29 +2.700000E+06 +3.000000E+04 +1.200000E+04 +1.018119E+12 +5.062500E+11 +7.290000E+11 27 -mass +0.000000E+00

# Element "capbeam": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 143 29 32 +2.700000E+06 +3.000000E+04 +1.200000E+04 +1.018119E+12 +5.062500E+11 +7.290000E+11 28 -mass +0.000000E+00

# Element "capbeam": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 144 32 30 +2.700000E+06 +3.000000E+04 +1.200000E+04 +1.018119E+12 +5.062500E+11 +7.290000E+11 29 -mass +0.000000E+00

# Element "bearing1": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 146 39 29 1 +2.409744E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "bearing1": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 148 37 32 1 +2.409744E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "bearing1": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 150 35 30 1 +2.409744E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "bearing1": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 152 36 31 1 +2.409744E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "bearing1": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 154 38 28 1 +2.409744E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "column": eleTag NodeI NodeJ NIP secTag geoTranTag <-mass massDens> <-cMass>
element dispBeamColumn 155 26 28 5 2 2 -mass +0.000000E+00 -integration Legendre

# Element "column": eleTag NodeI NodeJ NIP secTag geoTranTag <-mass massDens> <-cMass>
element dispBeamColumn 156 29 27 5 2 30 -mass +0.000000E+00 -integration Legendre

# Element "capbeam2": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 157 72 73 +2.340000E+06 +3.000000E+04 +1.200000E+04 +7.410189E+11 +3.295500E+11 +6.318000E+11 31 -mass +0.000000E+00

# Element "capbeam2": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 158 73 75 +2.340000E+06 +3.000000E+04 +1.200000E+04 +7.410189E+11 +3.295500E+11 +6.318000E+11 32 -mass +0.000000E+00

# Element "capbeam2": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 170 78 76 +2.340000E+06 +3.000000E+04 +1.200000E+04 +7.410189E+11 +3.295500E+11 +6.318000E+11 33 -mass +0.000000E+00

# Element "capbeam2": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 171 76 74 +2.340000E+06 +3.000000E+04 +1.200000E+04 +7.410189E+11 +3.295500E+11 +6.318000E+11 34 -mass +0.000000E+00

# Element "capbeam2": eleTag NodeI NodeJ A E G J Iy Iz geoTranTag <-mass massDens> <-cMass>
element elasticBeamColumn 172 74 72 +2.340000E+06 +3.000000E+04 +1.200000E+04 +7.410189E+11 +3.295500E+11 +6.318000E+11 35 -mass +0.000000E+00

# Element "bearing2": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 174 83 76 2 +1.658810E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "bearing2": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 176 81 74 2 +1.658810E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "bearing2": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 178 79 72 2 +1.658810E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "bearing2": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 179 73 80 2 +1.658810E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

# Element "bearing2": eleTag NodeI NodeJ frnMdlTag kInit -P matTag -T matTag -My matTag -Mz matTag <-orient <x1 x2 x3> y1 y2 y3> <-shearDist sDratio> <-doRayleigh> <-mass m> <-iter maxIter tol>
element flatSliderBearing 181 75 82 2 +1.658810E+03 -P 1 -T 1 -My 1 -Mz 1 -orient +1.000000E+00 +0.000000E+00 +0.000000E+00 -shearDist +0.000000E+00 -doRayleigh -iter 20 +1.000000E-12

Hello,I am building a model of four span continuous beam bridge recently.When I run the analysis of eigen,some errors has appeared. I have put my tcl.file here,can you give me some suggestions of this model, tank you !
OpenSees > # modify2.tcl
OpenSees > logFile "modify2.log"
OpenSees > set systemTime [clock seconds]
1600222769
OpenSees >
OpenSees > puts "==============================================================================="
===============================================================================
OpenSees > puts "OpenSees Navigator: Version 2.5.7"
OpenSees Navigator: Version 2.5.7
OpenSees > puts "File sourced: E:/µ²¿é/È«ÇÅÄ£Äâ/1/modify/modify2.tcl"
File sourced: E:/碌虏驴茅/脠芦脟脜脛拢脛芒/1/modify/modify2.tcl
OpenSees > puts "Starting OpenSees analysis: [clock format $systemTime -format "%d-%b-%Y %H:%M:%S"]"
Starting OpenSees analysis: 16-Sep-2020 10:19:29
OpenSees > puts "==============================================================================="
===============================================================================
OpenSees >
OpenSees > puts "\n########################"

########################
OpenSees > puts "# Analysis-Sequence 1 #"
# Analysis-Sequence 1 #
OpenSees > puts "########################"
########################
OpenSees >
OpenSees > # Start timing of this analysis sequence
OpenSees > # --------------------------------------
OpenSees > set tStart [clock clicks -milliseconds]
1600222769996
OpenSees >
OpenSees > # Start of model generation
OpenSees > # =========================
OpenSees >
OpenSees > # Create ModelBuilder
OpenSees > # -------------------
OpenSees > model BasicBuilder -ndm 3 -ndf 6
OpenSees >
OpenSees > # Define geometry
OpenSees > # ---------------
OpenSees > source NodeCoord.tcl
OpenSees >
OpenSees > # Define Single Point Constraints
OpenSees > # -------------------------------
OpenSees > source SPConstraint.tcl
42 43 44 45 46 47
OpenSees >
OpenSees > # Define nodal masses
OpenSees > # -------------------
OpenSees > source NodeMass.tcl
OpenSees >
OpenSees > # Define Multi Point Constraints
OpenSees > # ------------------------------
OpenSees > source MPConstraint.tcl
OpenSees >
OpenSees > # Define material(s)
OpenSees > # ------------------
OpenSees > source Materials.tcl
OpenSees >
OpenSees > # Define friction model(s)
OpenSees > # ------------------------
OpenSees > source FrictionModels.tcl
OpenSees >
OpenSees > # Define section(s)
OpenSees > # -----------------
OpenSees > source Sections.tcl
OpenSees >
OpenSees > # Define geometric transformation(s)
OpenSees > # ----------------------------------
OpenSees > source GeoTran.tcl
OpenSees >
OpenSees > # Define element(s)
OpenSees > # -----------------
OpenSees > source Elements.tcl
expected floating-point number but got "-shearDist"
OpenSees >
OpenSees > # Define damping parameters
OpenSees > # -------------------------
OpenSees >
OpenSees > # Define time series
OpenSees > # ------------------
OpenSees > source TimeSeries.tcl
OpenSees >
OpenSees > # Start of analysis generation
OpenSees > # ============================
OpenSees >
OpenSees > # Get Initial Stiffness
OpenSees > # ---------------------
OpenSees > initialize
OpenSees >
OpenSees > puts "o Analysis: EigenDefaultCase"
o Analysis: EigenDefaultCase
OpenSees > # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
OpenSees >
OpenSees > # Define load pattern
OpenSees > # -------------------
OpenSees > source LoadPattern_1.tcl
OpenSees >
OpenSees > # Define recorder(s)
OpenSees > # ------------------
OpenSees > source Recorder_1.tcl
2
OpenSees >
OpenSees > # Define analysis options
OpenSees > # -----------------------
OpenSees > source AnalysisOptn_1.tcl
OpenSees >
OpenSees > # Analyze model
OpenSees > # -------------
OpenSees > set pi [expr acos(-1.0)]
3.141592653589793
OpenSees > set eigFID [open EigenDefaultCase_Node_EigenVectors_EigenVal.out w]
file91a820
OpenSees > set lambda [eigen -fullGenLapack 3]
WARNING PlainHandler::handle() - constraint matrix not identity, ignoring constraint for node 72
non-varyng assumed
WARNING PlainHandler::handle() - constraint matrix not identity, ignoring constraint for node 93
non-varyng assumed
WARNING PlainHandler::handle() - constraint matrix not identity, ignoring constraint for node 103
non-varyng assumed
WARNING PlainHandler::handle() - constraint matrix not identity, ignoring constraint for node 105
non-varyng assumed
WARNING PlainHandler::handle() - constraint matrix not identity, ignoring constraint for node 110
non-varyng assumed
WARNING PlainHandler::handle() - constraint matrix not identity, ignoring constraint for node 112
non-varyng assumed
-479898.10073334811022505164 -678.40598248962487559766 5.30501512295851362921
OpenSees > puts $eigFID " lambda omega period frequency"
OpenSees > foreach lambda $lambda {
set omega [expr sqrt($lambda)]
set period [expr 2.0*$pi/$omega]
set frequ [expr 1.0/$period]
puts $eigFID [format " %+1.6e %+1.6e %+1.6e %+1.6e" $lambda $omega $period $frequ]
}
domain error: argument not in valid range
OpenSees > close $eigFID
OpenSees >
OpenSees > # Record eigenvectors
OpenSees > # -------------------
OpenSees > record
OpenSees >
OpenSees > # Stop timing of this analysis sequence
OpenSees > # -------------------------------------
OpenSees > set tStop [clock clicks -milliseconds]
1600222770588
OpenSees > puts "o Time taken: [expr ($tStop-$tStart)/1000.0] sec"
o Time taken: 0.592 sec
OpenSees >
OpenSees > # Clean up
OpenSees > # --------
OpenSees > wipe
OpenSees > set systemTime [clock seconds]
1600222770
OpenSees > puts "\n==============================================================================="

===============================================================================
OpenSees > puts "Finished OpenSees analysis: [clock format $systemTime -format "%d-%b-%Y %H:%M:%S"]"
Finished OpenSees analysis: 16-Sep-2020 10:19:30