内田混凝土方钢管悬臂柱高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

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