OpenSEES方钢管混凝土滞回性能模拟的疑问
kqai_16546
kqai_16546 Lv.7
2015年09月10日 14:30:00
只看楼主

内田混凝土方钢管悬臂柱高800mm,横截面长、宽均200mm,钢管壁厚3mm;内填C50混凝土(fc=38.5 MPa,Ec0=35310 MPa)。钢屈服强度238.4 MPa,弹模210000 MPa。轴力650 kN恒定。混凝土本构采用修正的Mander约束混凝土模型。以下为OpenSEES分析的命令流:# set up -----------------------------------------------------------------------

内田混凝土方钢管悬臂柱高800mm,横截面长、宽均200mm,钢管壁厚3mm;内填C50混凝土(fc=38.5 MPa,Ec0=35310 MPa)。钢屈服强度238.4 MPa,弹模210000 MPa。
轴力650 kN恒定。
混凝土本构采用修正的Mander约束混凝土模型。

以下为OpenSEES分析的命令流:

# set up -----------------------------------------------------------------------
wipe;
model BasicBuilder -ndm 2 -ndf 3;
#set dataDir Data;
#file mkdir $dataDir/;

# define units -----------------------------------------------------------------
set N 1;
set mm 1;
set sec 1;

# nodal coordinates ----------------------------------------------------------
node 1 0 0;
node 2 0 800;

# boundary conditions -------------------------------------------------------
fix 1 1 1 1;

# define materials ------------------------------------------------------------
set TAGsteel 1;
set TAGcon 2;
set fy 238.4; #yield strength of steel
set Es 210000;
set b 0.01; #strain hardening ratio of steel
set R0 18.5;
set cR1 0.925;
set cR2 0.15;
uniaxialMaterial Steel02 $TAGsteel $fy $Es $b $R0 $cR1 $cR2;
set fcc -53.529; #compressive strength of confined concrete
set ep -0.0053; #compressive strain of confined concrete
set fcu -42.656; #ultimate compressive strength of confined concrete
set epu -0.016; #ultimate compressive strain of confined concrete
set lambda 0.15;
set ftu 5.353; #ultimate tensile strength of confined concrete
set Ets 15512.09;
uniaxialMaterial Concrete02 $TAGcon $fcc $ep $fcu $epu $lambda $ftu $Ets;

# section discretization ----------------------------------------------------------
set TAGsec 3;
section Fiber $TAGsec {
patch quad $TAGsteel 3 1 -97 -97 -97 -100 100 -100 100 -97
patch quad $TAGsteel 1 3 97 -97 100 -97 100 100 97 100
patch quad $TAGsteel 3 1 -100 97 97 97 97 100 -100 100
patch quad $TAGsteel 1 3 -100 97 -100 -100 -97 -100 -97 97

patch quad $TAGcon 3 3 -97 -97 97 -97 97 97 -97 97
};

# geometric-transformation considering P-Delta effect ------------------
set TAGcrdtran 5;
geomTransf PDelta $TAGcrdtran;

# define nonlinear beam-column element ---------------------------------
set TAGele 6;
element nonlinearBeamColumn $TAGele 1 2 5 $TAGsec $TAGcrdtran;

# record results ---------------------------------------------------------------
recorder Node -file out.txt -node 1 -dof 1 reaction -node 2 -dof 1 disp;

# define loads ----------------------------------------------------------------
set TAGF 8;
pattern Plain $TAGF Linear {
load 2 1 -650000 0;
};

# problem solving ------------------------------------------------------------
constraints Plain;
numberer Plain;
system UmfPack;
test EnergyIncr 1.0e-4 5000;
algorithm Newton;
analysis Static;
integrator DisplacementControl 2 1 0.5;
analyze 100;

运行结果是一堆警告:
WARNING analysis Static - no Integrator specified,
StaticIntegrator default will be used
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 1
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 14311.9
OpenSees > analyze failed, returned: -3 error flag

应该说不是很复杂的一个题,怎么就得不出结果呢?
请各位高手不灵赐教!谢谢!
免费打赏
温故而知新
2015年10月16日 14:59:09
2楼
:handshake:handshake
回复
lrywjy
2015年11月04日 10:40:12
3楼
楼主解决了么?我也碰到了这个问题。。

abby菠萝

请问你现在解决了吗?我的模型也是这样不知道怎么办。

回复

相关推荐

APP内打开