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