opensees模拟滞回

wipeputs "System"model basic -ndm 2 -ndf 3puts "restraint"node 1 0 0 node 2 0 700node 3 0 1400node 4 0 700node 5 475 700node 6 1250 700node 7 1500 700node 8 1650 700puts "rigidDiaphragm"puts "mass"#mass 1 3.361E-002 3.361E-002 0.000E+000 0.000E+000 0.000E+000 0.000E+000#mass 2 8.130E-002 8.130E-002 0.000E+000 0.000E+000 0.000E+000 0.000E+000#mass 3 3.361E-002 3.361E-002 0.000E+000 0.000E+000 0.000E+000 0.000E+000#mass 4 4.314E-002 4.314E-002 0.000E+000 0.000E+000 0.000E+000 0.000E+000#mass 5 3.845E-002 3.845E-002 0.000E+000 0.000E+000 0.000E+000 0.000E+000#mass 6 9.379E-003 9.379E-003 0.000E+000 0.000E+000 0.000E+000 0.000E+000puts "node"fix 1 1 1 0;fix 3 1 0 0;fix 8 1 0 0;puts "Equal DOF" equalDOF 2 4 1 2puts "material"#钢筋本构uniaxialMaterial Steel02 1 582 289850 0.01 15 0.925 0.15 0.04 1 0.04 1 0 #uniaxialMaterial Steel01 1 582 289850 0.01#柱核心混凝土本构 #01本构 #uniaxialMaterial Concrete01 2 -57.22 -0.0052 -28.61 -0.0445 #（mander本构模型） #uniaxialMaterial Concrete02 2 -57.22 -0.0052 -28.61 -0.0445 0.1 4.08 1000 uniaxialMaterial Concrete07 2 -57.22 -0.0052 32939 4.08 0.00025 2 8.58 1.5 #（modifed kent-park本构模型） #uniaxialMaterial Concrete02 2 -53.38 -0.0025 -10.68 -0.0502 0.1 4.08 1000#柱保护层混凝土本构 #01本构 uniaxialMaterial Concrete01 3 -43.4 -0.002 -8.47 -0.008 # mander本构模型） #uniaxialMaterial Concrete02 3 -43.4 -0.002 -8.47 -0.008 0.1 4.08 1000 #uniaxialMaterial Concrete07 3 -43.4 -0.002 32939 4.08 0.00025 2 2.3 2.93 #（过镇海本构模型） #uniaxialMaterial Concrete02 3 -43.4 -0.0027 -8.47 -0.008 0.1 4.08 1000#梁加密区（1）核心混凝土本构 #01本构 #uniaxialMaterial Concrete01 4 -47.13 -0.0029 -23.56 -0.011 #（mander本构模型） #uniaxialMaterial Concrete02 4 -47.13 -0.0029 -23.56 -0.011 0.1 4.08 1000 uniaxialMaterial Concrete07 4 -47.13 -0.0029 32939 4.08 0.00025 2 3.85 2 #（modifed kent-park本构模型） #uniaxialMaterial Concrete02 4 -45.36 -0.0021 -9.07 -0.0155 0.1 4.08 1000#梁加密区（1）保护层混凝土本构 #01本构 uniaxialMaterial Concrete01 5 -43.4 -0.002 -8.47 -0.008 #（mander本构模型） #uniaxialMaterial Concrete02 5 -43.4 -0.002 -8.47 -0.008 0.1 4.08 1000 #uniaxialMaterial Concrete07 5 -43.4 -0.002 32939 4.08 0.00025 2 2.3 2.93 #（过镇海本构模型） #uniaxialMaterial Concrete02 5 -43.4 -0.0027 -8.47 -0.008 0.1 4.08 1000#梁非加密区核心混凝土本构 #01本构 #uniaxialMaterial Concrete01 6 -43.95 -0.0021 -21.98 -0.0055 #（mander本构模型） #uniaxialMaterial Concrete02 6 -43.95 -0.0021 -21.98 -0.0055 0.1 4.08 1000 uniaxialMaterial Concrete07 6 -43.95 -0.0021 32939 4.08 0.00025 2 2.59 2.68 #（modifed kent-park本构模型） #uniaxialMaterial Concrete02 6 -45.24 -0.0021 -9.05 -0.0082 0.1 4.08 1000#梁非加密区保护层混凝土本构 #01本构 uniaxialMaterial Concrete01 7 -43.4 -0.002 -8.47 -0.008 #（mander本构模型） #uniaxialMaterial Concrete02 7 -43.4 -0.002 -8.47 -0.008 0.1 4.08 1000 #uniaxialMaterial Concrete07 7 -43.4 -0.002 32939 4.08 0.00025 2 2.3 2.93 #（过镇海本构模型） #uniaxialMaterial Concrete02 7 -43.4 -0.0027 -8.47 -0.008 0.1 4.08 1000 #梁加密区（2）核心混凝土本构 #01本构 #uniaxialMaterial Concrete01 8 -47.13 -0.0029 -23.56 -0.011 #（mander本构模型） #uniaxialMaterial Concrete02 8 -47.13 -0.0029 -23.56 -0.011 0.1 4.08 1000 uniaxialMaterial Concrete07 8 -47.13 -0.0029 32939 4.08 0.00025 2 3.85 2 #（modifed kent-park本构模型） #uniaxialMaterial Concrete02 8 -45.36 -0.0021 -9.07 -0.0155 0.1 4.08 1000#梁加密区（2）保护层混凝土本构 #01本构 uniaxialMaterial Concrete01 9 -43.4 -0.002 -8.47 -0.008 #（mander本构模型） #uniaxialMaterial Concrete02 9 -43.4 -0.002 -8.47 -0.008 0.1 4.08 1000 #uniaxialMaterial Concrete07 9 -43.4 -0.002 32939 4.08 0.00025 2 2.3 2.93 #（过镇海本构模型） #uniaxialMaterial Concrete02 9 -43.4 -0.0027 -8.47 -0.008 0.1 4.08 1000 #板核心混凝土本构 #01本构 #uniaxialMaterial Concrete01 10 -43.95 -0.0021 -21.98 -0.0055 #（mander本构模型） #uniaxialMaterial Concrete02 10 -43.95 -0.0021 -21.98 -0.0055 0.1 4.08 1000 uniaxialMaterial Concrete07 10 -43.95 -0.0021 32939 4.08 0.00025 2 2.59 2.68 #（modifed kent-park本构模型） #uniaxialMaterial Concrete02 10 -45.24 -0.0021 -9.05 -0.0082 0.1 4.08 1000#板保护层混凝土本构 #01本构 uniaxialMaterial Concrete01 11 -43.4 -0.002 -8.47 -0.008 #（mander本构模型） #uniaxialMaterial Concrete02 11 -43.4 -0.002 -8.47 -0.008 0.1 4.08 1000 #uniaxialMaterial Concrete07 11 -43.4 -0.002 32939 4.08 0.00025 2 2.3 2.93 #（过镇海本构模型） #uniaxialMaterial Concrete02 11 -43.4 -0.0027 -8.47 -0.008 0.1 4.08 1000 #零长度单元 #（钢筋滑移本构） uniaxialMaterial Bond_SP01 100 582 0.34027 878.2 11.9094 0.25 0.25 #uniaxialMaterial BarSlip 100 30 582 289850 855 2898.5 200 8 4 125 250 weak beambot #核心混凝土本构 #01本构 #uniaxialMaterial Concrete01 200 -47.13 -0.0029 -23.56 -0.011 #（mander本构模型） #uniaxialMaterial Concrete02 200 -47.13 -0.0029 -23.56 -0.011 0.1 4.08 1000 uniaxialMaterial Concrete07 200 -47.13 -0.0029 32939 4.08 0.00025 2 3.85 2 #（modifed kent-park本构模型） #uniaxialMaterial Concrete02 200 -45.36 -0.0021 -9.07 -0.0155 0.1 4.08 1000 #保护层混凝土本构 #01本构 uniaxialMaterial Concrete01 300 -43.4 -0.002 -8.47 -0.008 #（mander本构模型） #uniaxialMaterial Concrete02 300 -43.4 -0.002 -8.47 -0.008 0.1 4.08 1000 #uniaxialMaterial Concrete07 300 -43.4 -0.002 32939 4.08 0.00025 2 2.3 2.93 #（过镇海本构模型） #uniaxialMaterial Concrete02 300 -43.4 -0.0027 -8.47 -0.008 0.1 4.08 1000#柱截面抗剪抗扭uniaxialMaterial Elastic 201 3.447E+008uniaxialMaterial Elastic 301 3.447E+008uniaxialMaterial Elastic 401 2.330E+012#梁截面抗剪抗扭uniaxialMaterial Elastic 202 3.228E+008uniaxialMaterial Elastic 302 1.67856E+008uniaxialMaterial Elastic 402 1.03635E+012#柱子纤维模型NC200X200section Fiber 1 {#核心混凝土patch rect 2 5 5 -75 -75 75 75#cover concrete fiberspatch rect 3 5 1 -75 75 75 100patch rect 3 1 7 75 -100 100 100patch rect 3 5 1 -75 -100 75 -75patch rect 3 1 7 -100 -100 -75 100#reinforcing fiberslayer straight 1 3 50.27 -75 75 75 75 layer straight 1 2 50.27 -75 0 75 0layer straight 1 3 50.27 -75 -75 75 -75}#梁加密区（1）纤维模型section Fiber 2 {#核心混凝土patch rect 4 3 8 -37.5 -100 37.5 100#cover concrete fiberspatch rect 5 3 1 -37.5 100 37.5 125patch rect 5 1 10 37.5 -125 62.5 125patch rect 5 3 1 -37.5 -125 37.5 -100patch rect 5 1 10 -62.5 -125 -37.5 125#板核心区混凝土patch rect 10 10 1 62.5 85 422.5 105patch rect 10 10 1 -422.5 85 -62.5 105 #板保护层区混凝土patch rect 11 10 1 62.5 65 422.5 85patch rect 11 10 1 62.5 105 422.5 125patch rect 11 10 1 -422.5 65 -62.5 85patch rect 11 10 1 -422.5 105 -62.5 125#reinforcing fiberslayer straight 1 4 50.27 -37.5 -100 37.5 -100layer straight 1 4 78.54 -37.5 100 37.5 100layer straight 1 4 12.57 -422.5 85 -62.5 85layer straight 1 4 12.57 -422.5 105 -62.5 105layer straight 1 4 12.57 62.5 85 422.5 85layer straight 1 4 12.57 62.5 105 422.5 105}#梁非加密区纤维模型section Fiber 3 {#核心混凝土patch rect 6 3 8 -37.5 -100 37.5 100#cover concrete fiberspatch rect 7 3 1 -37.5 100 37.5 125patch rect 7 1 10 37.5 -125 62.5 125patch rect 7 3 1 -37.5 -125 37.5 -100patch rect 7 1 10 -62.5 -125 -37.5 125#板核心区混凝土patch rect 10 10 1 62.5 85 422.5 105patch rect 10 10 1 -422.5 85 -62.5 105 #板保护层区混凝土patch rect 11 10 1 62.5 65 422.5 85patch rect 11 10 1 62.5 105 422.5 125patch rect 11 10 1 -422.5 65 -62.5 85patch rect 11 10 1 -422.5 105 -62.5 125#reinforcing fiberslayer straight 1 4 50.27 -37.5 -100 37.5 -100layer straight 1 2 78.54 -37.5 100 37.5 100layer straight 1 4 12.57 -422.5 85 -62.5 85layer straight 1 4 12.57 -422.5 105 -62.5 105layer straight 1 4 12.57 62.5 85 422.5 85layer straight 1 4 12.57 62.5 105 422.5 105}#梁加密区（2）纤维模型section Fiber 4 {#核心混凝土patch rect 8 3 8 -37.5 -100 37.5 100#cover concrete fiberspatch rect 9 3 1 -37.5 100 37.5 125patch rect 9 1 10 37.5 -125 62.5 125patch rect 9 3 1 -37.5 -125 37.5 -100patch rect 9 1 10 -62.5 -125 -37.5 125#板核心区混凝土patch rect 10 10 1 62.5 85 422.5 105patch rect 10 10 1 -422.5 85 -62.5 105 #板保护层区混凝土patch rect 11 10 1 62.5 65 422.5 85patch rect 11 10 1 62.5 105 422.5 125patch rect 11 10 1 -422.5 65 -62.5 85patch rect 11 10 1 -422.5 105 -62.5 125#reinforcing fiberslayer straight 1 4 50.27 -37.5 -100 37.5 -100layer straight 1 2 78.54 -37.5 100 37.5 100layer straight 1 4 12.57 -422.5 85 -62.5 85layer straight 1 4 12.57 -422.5 105 -62.5 105layer straight 1 4 12.57 62.5 85 422.5 85layer straight 1 4 12.57 62.5 105 422.5 105}#零长度单元纤维模型section Fiber 5 {#核心混凝土patch rect 200 3 8 -37.5 -100 37.5 100#cover concrete fiberspatch rect 300 3 1 -37.5 100 37.5 125patch rect 300 1 10 37.5 -125 62.5 125patch rect 300 3 1 -37.5 -125 37.5 -100patch rect 300 1 10 -62.5 -125 -37.5 125#reinforcing fiberslayer straight 100 4 50.27 -37.5 -100 37.5 -100layer straight 100 4 78.54 -37.5 100 37.5 100}#柱子section Aggregator 1001 201 Vy 301 Vz 401 T -section 1#加密区（1）section Aggregator 1002 202 Vy 302 Vz 402 T -section 2#非加密区section Aggregator 1003 202 Vy 302 Vz 402 T -section 3#加密区（2）section Aggregator 1004 202 Vy 302 Vz 402 T -section 4#零长度section Aggregator 1005 202 Vy 302 Vz 402 T -section 5puts "section"puts "transformation"geomTransf PDelta 1 #geomTransf Linear 1 geomTransf Linear 2 puts "element"set np 4#element dispBeamColumn 1 1 2 $np 1001 1#element dispBeamColumn 2 2 3 $np 1001 1#element dispBeamColumn 3 4 5 $np 1002 2#element dispBeamColumn 4 5 6 $np 1003 2#element dispBeamColumn 5 6 7 $np 1004 2#element dispBeamColumn 6 7 8 $np 1004 2element nonlinearBeamColumn 1 1 2 $np 1001 1element nonlinearBeamColumn 2 2 3 $np 1001 1element nonlinearBeamColumn 3 4 5 $np 1002 2element nonlinearBeamColumn 4 5 6 $np 1003 2element nonlinearBeamColumn 5 6 7 $np 1004 2element nonlinearBeamColumn 6 7 8 $np 1004 2element zeroLengthSection 7 2 4 1005puts "recorder"recorder Node -file node7.out -time -node 7 -dof 2 dispputs "gravity"## Load Case = DEADset p -1.66E+005#set p 0pattern Plain 1 Linear {load 3 0.000E+000 $p 0.000E+000}puts "analysis"constraints Plainnumberer Plainsystem BandGeneraltest NormDispIncr 1.0e-1 200algorithm Newtonintegrator LoadControl 1.000E-001analysis Staticanalyze 10loadConst -time 0.0puts " pushover"## Load Case = PUSHpattern Plain 2 Linear {load 7 0.000E+000 1.000E+005 0.000E+000 }puts "analysis"constraints Plain numberer RCMsystem UmfPacktest NormDispIncr 1.0e-1 200algorithm Newtonanalysis Staticintegrator DisplacementControl 7 2 -0.002analyze 1000integrator DisplacementControl 7 2 0.006analyze 1000integrator DisplacementControl 7 2 -0.012analyze 1000integrator DisplacementControl 7 2 0.018analyze 1000integrator DisplacementControl 7 2 -0.022analyze 1000integrator DisplacementControl 7 2 0.028analyze 1000integrator DisplacementControl 7 2 -0.036analyze 1000integrator DisplacementControl 7 2 0.040analyze 1000integrator DisplacementControl 7 2 -0.040analyze 1000integrator DisplacementControl 7 2 0.040analyze 1000integrator DisplacementControl 7 2 -0.050analyze 1000integrator DisplacementControl 7 2 0.060analyze 1000integrator DisplacementControl 7 2 -0.060analyze 1000integrator DisplacementControl 7 2 0.060analyze 1000integrator DisplacementControl 7 2 -0.080analyze 1000integrator DisplacementControl 7 2 0.100analyze 1000integrator DisplacementControl 7 2 -0.100analyze 1000integrator DisplacementControl 7 2 0.100analyze 1000integrator DisplacementControl 7 2 -0.105analyze 1000integrator DisplacementControl 7 2 0.110analyze 1000integrator DisplacementControl 7 2 -0.110analyze 1000integrator DisplacementControl 7 2 0.110analyze 1000integrator DisplacementControl 7 2 -0.121analyze 1000integrator DisplacementControl 7 2 0.132analyze 1000integrator DisplacementControl 7 2 -0.132analyze 1000integrator DisplacementControl 7 2 0.132analyze 1000integrator DisplacementControl 7 2 -0.143analyze 1000integrator DisplacementControl 7 2 0.143analyze 1000integrator DisplacementControl 7 2 -0.143analyze 1000integrator DisplacementControl 7 2 0.143analyze 1000integrator DisplacementControl 7 2 -0.167analyze 1000integrator DisplacementControl 7 2 0.180analyze 1000integrator DisplacementControl 7 2 -0.180analyze 1000integrator DisplacementControl 7 2 0.180analyze 1000integrator DisplacementControl 7 2 -0.190analyze 1000integrator DisplacementControl 7 2 0.210analyze 1000integrator DisplacementControl 7 2 -0.220analyze 1000integrator DisplacementControl 7 2 0.220analyze 1000integrator DisplacementControl 7 2 -0.220analyze 1000integrator DisplacementControl 7 2 0.220analyze 1000integrator DisplacementControl 7 2 -0.220analyze 1000integrator DisplacementControl 7 2 0.230analyze 1000