土木在线论坛 \ 建筑结构 \ 结构设计软件 \ OpenSEES方钢管混凝土滞回性能模拟的疑问

OpenSEES方钢管混凝土滞回性能模拟的疑问

发布于:2015-09-10 14:30:10 来自:建筑结构/结构设计软件 [复制转发]
内田混凝土方钢管悬臂柱高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

应该说不是很复杂的一个题,怎么就得不出结果呢?
请各位高手不灵赐教!谢谢!

全部回复(3 )

只看楼主 我来说两句抢地板
  • lrywjy
    lrywjy 沙发
    楼主解决了么?我也碰到了这个问题。。
    2015-11-04 10:40:04

    回复 举报
    赞同0
  • 温故而知新
    :handshake:handshake
    2015-10-16 14:59:16

    回复 举报
    赞同0
这个家伙什么也没有留下。。。

结构设计软件

返回版块

41.5 万条内容 · 265 人订阅

猜你喜欢

阅读下一篇

OpenSees框架拟静力模拟的问题

最近做框架拟静力模拟,遇到几个问题,苦于身边没有用这个软件的朋友,只好拿到这来请教高手。1.有没有办法实现先荷载控制到屈服-后位移控制到破坏?因为官网里面看到有关加载例子只有三种:重力、pushover、地震动。唯一的1个cyclic load位移控制的例子也没有说清屈服位移怎么确定。2.我用下列代码做拟静力位移控制加载,结果报错expected floating-point number but got “$a”

回帖成功

经验值 +10