:--------------------------------------------------------------------------------------------------------------------
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 1
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis:: analyze() - the Algorithm failed at iteration: 0 with domain at load factor 0.1
OpenSees > analyze failed, returned: -3 error flag
--------------------------------------------------------------------------------------------------------------------
i've tried to change the type of the solve system but it didn't work.here is my scripts,hope that you could give my some suggestions,Thank you!
---------------------------------------------------------------------------------------------------------------------
#Units:KN,M
wipe;
model Basic -ndm 2 -ndf 3;
#set some parameters for the model geomtry
set Length [expr 1.5*$m];
set Length1 [expr 0.5*$m];
set Length2 [expr 0.75*$m];
set Length3 [expr 1.0*$m];
node 1 0 0;
node 2 $Length 0;
node 3 $Length1 0;
node 4 $Length2 0;
node 5 $Length3 0;
fix 1 1 1 1;
fix 2 1 1 1;
#Define the Materials for the Beam
uniaxialMaterial Concrete01 1 -23.4 -0.002 8.6 -0.003;
uniaxialMaterial Steel01 2 400 2.1e5 0.01;
#parameters of the section
set BeaWidth 120*$mm;
set BeaDepth 200*$mm;
set Cover 25*$mm;
set As [expr 50.3*$mm2];
set y1 [expr $BeaWidth/2.0];
set z1 [expr $BeaDepth/2.0];
section Fiber 1 {
patch rect 1 1 10 [expr $Cover-$y1] [expr $Cover-$z1] [expr $y1-$Cover] [expr $z1-$Cover];
patch rect 2 7 1 [expr -$y1] [expr $z1-$Cover] $y1 $z1;
patch rect 2 7 1 [expr -$y1] [expr -$z1] $y1 [expr $Cover-$z1];
patch rect 2 1 10 [expr -$y1] [expr $Cover-$z1] [expr $Cover-$y1] [expr $z1-$Cover];
patch rect 2 1 10 [expr $y1-$Cover] [expr $Cover-$y1] $y1 [expr $z1-$Cover];
layer straight 2 2 $As [expr $Cover-$y1] [expr $z1-$Cover] [expr $Cover-$y1] [expr $Cover-$z1];
layer straight 2 2 $As [expr $y1-$Cover] [expr $z1-$Cover] [expr $y1-$Cover] [expr $Cover-$z1];
};
geomTransf Linear 1;
element forceBeamColumn 1 1 2 1 Lobatto 1 5;
set F 20.43;
timeSeries Linear 1
pattern Plain 1 1 {
load 3 0.0 [expr -$F] 0.0
load 5 0.0 [expr -$F] 0.0
}
recorder Node -file node4.out -time -node 4 -dof 2 disp;
recorder Element -file Element1_Steel.out -time -ele 1 section 1 fiber 0.035 -0.075 stressstrain;
recorder Element -file Element1_Concrete.out -time -ele 1 section 1 fiber -0.035 -0.003 stressstrain;
system BandGeneral;
numberer Plain;
constraints Plain;
test NormDispIncr 1.0e-12 10 3;
algorithm Newton;
integrator LoadControl 0.1;
analysis Static;
analyze 100;