RC柱低周反复分析问题-混凝土抗压强度影响，求助！

westinglie

刚刚发帖失败了，再发下，求助！
最近学习了OpenSEES模拟RC柱的低周反复滞回曲线，现在才刚开始做对比柱就出现问题了。混凝土本构关系设置中，抗压强度fc对计算结果很大，我把fc设置为35MPa时，计算结果好，抗压强度我只要调整一点点，就算不出来了，如调为34MPa或36MPa，滞回曲线就不像样了，滞回曲线不对称、或直接乱了，不知道什么原因。求解！如果碰到滞回曲线很不规则，一般要怎么调整相关参数设置？使滞回曲线漂亮些。 谢谢！
#=======tcl文件如下==============
# 单位制：N , mm
wipe;
model BasicBuilder -ndm 2 -ndf 3;
file mkdir Data;


node 1 0 0;
node 2 0 0;
node 3 0 1175;
fix 1 1 1 1;
fix 2 1 1 0;

set Fy 350;
set Fu [expr 1.2*$Fy];
set Es 200000;
set Esh [expr 0.0255*$Es];
set esh 0.01152;
set eult 0.1105;
uniaxialMaterial ReinforcingSteel 1 $Fy $Fu $Es $Esh $esh $eult;

set fc -35;      # 混凝土抗压强度设置
set Ec 30744.12;
# unconfined concrete
set fc1u $fc;
set eps1u -0.002;
set eps2u -0.004;
set Ecu $Ec;
uniaxialMaterial Concrete04 2 $fc1u $eps1u $eps2u $Ecu;
# confined concrete
set fc1c [expr 1.09*$fc];
set eps1c -0.0034;
set eps2c -0.062;
set Ecc 32826.12;
uniaxialMaterial Concrete04 3 $fc1c $eps1c $eps2c $Ecc;

# zerolength section element reinforcement
set Bzerolength 0.4;
set Rzerolength 0.6;
set zerolengthec -0.0075;
set zerolengthecu -0.15;
uniaxialMaterial Bond_SP01 4 $Fy 0.56 $Fu 22.4 $Bzerolength $Rzerolength;

section fiberSec 101 {
     patch quadr 3 4 16 -105 105 -105 -105 105 -105 105 105;
     patch quadr 2 1 16 -125 125 -105 105 105 105 125 125;
     patch quadr 2 1 16 -105 -105 -125 -125 125 -125 105 -105;
     patch quadr 2 4 1 -125 125 -125 -125 -105 -105 -105 105;
     patch quadr 2 4 1 105 105 105 -105 125 -125 125 125;
     layer straight 1 3 153.94 -105 105 -105 -105;
     layer straight 1 3 153.94 105 105 105 -105;
};

# zeroLengthSection
section fiberSec 102 {
     patch quadr 3 4 16 -105 105 -105 -105 105 -105 105 105;
     patch quadr 2 1 16 -125 125 -105 105 105 105 125 125;
     patch quadr 2 1 16 -105 -105 -125 -125 125 -125 105 -105;
     patch quadr 2 4 1 -125 125 -125 -125 -105 -105 -105 105;
     patch quadr 2 4 1 105 105 105 -105 125 -125 125 125;
     layer straight 4 3 153.94 -105 105 -105 -105;
     layer straight 4 3 153.94 105 105 105 -105;
};

geomTransf Linear 1;
element nonlinearBeamColumn 1 2 3 5 101 1;
element zeroLengthSection 2 1 2 102;

recorder Node -file Data/DFree.out -time -node 3 -dof 1 2 3 disp;
recorder Node -file Data/Base.out -time -node 2 -dof 1 2 3 reaction;
recorder Drift -file Data/Drift.out -time -iNode 3 -jNode 2 -dof 1   -perpDirn 2 ;
recorder plot Data/Drift.out test 0 0 400 400 -columns 2 1;

puts "Model Built";

pattern Plain 1 Linear {
     load 3 0 -150000 0;
};
constraints Plain;
numberer Plain;
system BandGeneral;
test NormDispIncr 1e-8 6;
algorithm Newton;
integrator LoadControl 0.1;
analysis Static;
analyze 10;

loadConst -time 0;

pattern Plain 2 Linear {
      load 3 150000 0 0;
};

constraints Plain;
numberer Plain;
system SparseGeneral -piv;
test NormDispIncr 1e-4 1000;
algorithm KrylovNewton;
analysis Static;

# Compute displacement increment
set dU1 0.05
set dU2 -0.05

# Perform the analysis
set IDctrlNode 3;
set IDctrlDOF 1;
integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU1
analyze 50
integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU2
analyze 100
integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU1
analyze 50

integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU1
analyze 100
integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU2
analyze 200
integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU1
analyze 100

integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU1
analyze 150
integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU2
analyze 300
integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU1
analyze 150

integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU1
analyze 200
integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU2
analyze 400
integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU1
analyze 200

integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU1
analyze 250
integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU2
analyze 500
integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU1
analyze 250

integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU1
analyze 300
integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU2
analyze 600
integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU1
analyze 300

integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU1
analyze 350
integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU2
analyze 700
integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU1
analyze 350

integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU1
analyze 400
integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU2
analyze 800
integrator DisplacementControl $IDctrlNode   $IDctrlDOF $dU1
analyze 400
# =======================================================

westinglie

刚摸了OpenSEES一段时间，对RC构件的模拟才刚开始，现在还不知道怎么有效的调试参数，碰到收敛问题，调试参数还没成功过，希望路过的大侠，知道一些解决RC构件模拟中，出现的不收敛、计算出错、结果不理想（如滞回曲线不规则）等问题，可以留下经验什么，现在一抹黑，谢谢！

ysy

为什么我ID禁止发言了？

kmani

楼主真悲剧，求助没有解决，反而还被禁言了，嘿嘿。

alvestian

关注一下。

406663989

我也有问题

tangkkktang

不仅仅是材料强度问题 还有后后面荷载等联系起来调整

jianku

问题比较多啊

宁静的sinx

持续关注，

lx3221

请问楼主的滞回曲线是怎么提取的？直接提取加载点的水平力和位移吗？