Error...BandGenLinLapackSolver::solve<>

[SunnyPark]

Hi,

I have questions about my error....
here is my error Massege;
=================================================================
WARNING BandGenLinLapackSolver::solve<> -LAPACK routine returned 2
WARNING NewtonRaphson::solveCurrentStep<> -the KnearSysOfEqm failed in solve<>
DirectIntegrationAnalysis::analyze<> - the Algorithm failed at time 22.56
OpenSees > analyze failed, returned: -3 errpr flag
-3
=================================================================

what i wonder is why this error comes out when i analyze my model with trasient analysis..

my model (at below) is analyzed well with static analysis (hysteretic-pushover)
and analyzed by dynamic analysis with EQ input file (with small peak acceleration (like 0.1g)) ; there is no error at all

But!!!!! when above errors comes out with EQ file (with lager peak accel than 0.2g)........................

I search about my error in this community . but i cannot find the clear answer... i don't have any zerolength-element.....and there is no problem with boundary condition.....

Please somebody help me............!!!!!!!!!!!!!



I attatched my model.

Thank you...
==================================================================
# As Built

wipe;
model BasicBuilder -ndm 2 -ndf 3;

# nodal coordinates:
node 1 0 0;
node 2 0 400;
node 3 0 700;
node 4 0 1100;
node 5 0 1400;

# Single point constraints -- Boundary Conditions
fix 1 1 1 1;

# Define ELEMENTS & SECTIONS -------------------------------------------------------------
set ColSecTag 1001;
set coverCol 40;
set barAreaCol 133;

set IDconcCore 1;
set IDconcCover 2;
set IDreinf 3;
set TorMatID 4;
set IDShear 5;

# Material
set dam "energy";

uniaxialMaterial Pinching4 $IDconcCover 6.31 0.0015 6.8 0.003 7.1 0.009 6.8 0.02 -5.81 -0.0015 -6.3 -0.003 -6.6 -0.009 -6.2 -0.02 0.5 0.5 0.05 0.5 0.5 0.05 2.0 2.0 2.5 3.0 2.0 0.0 0.0 1.0 1.0 1.0 0.0 1.0 1.0 1.3 1.3 10 $dam;
uniaxialMaterial Pinching4 $IDconcCore 6.31 0.0015 6.8 0.003 7.1 0.009 6.8 0.02 -5.81 -0.0015 -6.3 -0.003 -6.6 -0.009 -6.2 -0.02 0.5 0.5 0.05 0.5 0.5 0.05 2.0 2.0 2.5 3.0 2.0 0.0 0.0 1.0 1.0 1.0 0.0 1.0 1.0 1.3 1.3 10 $dam;

uniaxialMaterial Pinching4 6 6.3 0.002 7.6 0.005 7.5 0.009 7.0 0.1 -7.0 -0.002 -7.7 -0.005 -7.5 -0.009 -7.0 -0.1 0.1 0.2 0.0 0.7 0.7 0.0 0.0 0.0 0.0 0.0 0.0 0.5 0.5 1.0 1.0 0.0 0.0 1.0 1.5 0.0 0.0 10 $dam;

uniaxialMaterial Hysteretic $IDreinf 370. 0.002 380. 0.005 350. 0.1 -370. -0.002 -380. -0.005 -350. -0.1 1 1 0.025 0.0 0.8;

uniaxialMaterial Elastic $TorMatID 1.0e10;

# SECTION
section Fiber 101 {
patch circ $IDconcCore 32 8 0.0 0.0 0.0 160.0 0.0 360.0
patch circ $IDconcCover 32 2 0.0 0.0 160.0 200.0 0.0 360.0
layer circ $IDreinf 16 $barAreaCol 0.0 0.0 160.0 0.0 360.0
}

section Fiber 102 {
patch circ $IDconcCore 32 8 0.0 0.0 0.0 160.0 0.0 360.0
patch circ 6 32 2 0.0 0.0 160.0 200.0 0.0 360.0
layer circ $IDreinf 16 $barAreaCol 0.0 0.0 160.0 0.0 360.0
}

set ColTransfTag1 1;
set ColTransfTag2 2;
set ColTransfType Corotational;
geomTransf $ColTransfType $ColTransfTag1;
geomTransf Linear $ColTransfTag2 -jntOffset 0 0 0 300.;

section Aggregator 1 $TorMatID T -section 101;
section Aggregator 2 $TorMatID T -section 102;

# element connectivity:
set numIntgrPts 5;
element dispBeamColumn 1 1 2 $numIntgrPts 102 $ColTransfTag1;
element dispBeamColumn 2 2 3 $numIntgrPts 101 $ColTransfTag1;
element dispBeamColumn 3 3 4 $numIntgrPts 101 $ColTransfTag1;
element dispBeamColumn 4 4 5 $numIntgrPts 101 $ColTransfTag2;

puts "column assigned"

set W [expr -4142.077*$N]; # W(N)=Ag*h*2400kg/m3*9.81m/s2
set m [expr -$W/$g];
#set Wg [expr -376991.12*$N]; # Wg(N)=0.1*Ag*fc
set Wg [expr -263893.784*$N]; # Wg(N)=0.1*Ag*fc*0.7
set mg [expr -$Wg/$g];

pattern Plain 1 Constant {

load 3 0.0 [expr $W/2] 0.0
load 5 0.0 [expr $W/2+$Wg] 0.0

}

mass 2 [expr $m/1400.*350.] 0.0 0.0;
mass 3 [expr $m/1400.*350.] 0.0 0.0;
mass 4 [expr $m/1400.*350.] 0.0 0.0;
mass 5 [expr $m/1400.*150.+$mg] 0.0 0.0;

# Define RECORDERS -------------------------------------------------------------
recorder Node -file MA00Displ4GM.out -time -node 4 -dof 1 disp;
recorder Node -file MA00Displ5GM.out -time -node 5 -dof 1 disp;
recorder Node -file MA00AccelGM.out -time -node 5 -dof 1 accel;
recorder Node -file MA00ForceGM.out -time -node 1 -dof 1 reaction;

integrator LoadControl 0.1;

###### Ground motion analysis
set record_length 15000;
set dt 0.01;
set outFile gm1.g3;
set accelSeries "Path -filePath $outFile -dt $dt -factor [expr 3./5.*$g]"

source ReadSMDFile.tcl;
ReadSMDFile gm1.AT2 $outFile dt;
pattern UniformExcitation 2 1 -accel $accelSeries

rayleigh 0.0 0.0 0.0 [expr 2*0.0413/pow([eigen 1],0.5)];

puts "transient load defined";

constraints Transformation
numberer Plain
system BandGeneral
test EnergyIncr 1.0e-8 10
algorithm Newton
integrator Newmark 0.5 0.25
analysis Transient

#Result script
#Record

set lamda [eigen 1]
set omega [expr pow($lamda, 0.5)]
set Tperiod [expr 2.*$pi/$omega]

puts $Tperiod
print node 5


analyze $record_length $dt

=====================================================

Somebody help me to correct this errors???

[vesna]

comment out this line:

integrator LoadControl 0.1;

Solve your analysis in a loop like in this example:
http://opensees.berkeley.edu/wiki/index ... e_Analysis

[SunnyPark]

thanks vesna
i saw the example and fixed analysis command.
but it still have error..and.. error is become more complex....
and i have a question about the command that used in example,
ReadRecord command is not worked at all..
i'm using OpenSees2.2.2 version... should i upgrade it?

i upgraded to version 2.3.2 but there is error.

new error is here...
==========================

WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 5
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 19.86
OpenSees > analyze failed, returned: -3 error flag
regular newton failed .. lets try an initail stiffness for this step
WARNING: CTestNormDispIncr::test() - failed to converge
after: 100 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 19.86
OpenSees > analyze failed, returned: -3 error flag
Transient analysis completed FAILED
ArpackSolver::Error with _saupd info = -9
Starting vector is zero.
WARNING DirectIntegrationAnalysis::eigen() - EigenSOE failed in solve()
eigen values at end of transient:

Node: 5
Coordinates : 0 1400
Disps: -0.422325 -0.75974 0.00024608
Velocities : 2.73214 -0.017282 -0.0024139
commitAccels: 392.516 1223.05 0.180817
unbalanced Load: 6687.2 -268036 0
reaction: -10587.6 -0.0138018 -0.0630025
Mass :
26.9549 0 0
0 0 0
0 0 0

Rayleigh Factor: alphaM: 0
Eigenvectors:
-0.192396 0.00797152
0.000857914 -0.0149773
0.000194138 -0.00356224

ID : 0 1 2

[vesna]

I just noticed that you assigned torsion to your 2D section. 2D section can not have torsion assigned to it. Remove section aggregators from your model.

[SunnyPark]

thank you!
I didn't know that and removed section aggregators!
But!!!!!!!!!!!!!!!!!!!!!still have same error.........
I am just goning carzy..........T^T
How can I fixed this error...???????

[fmk]

you have a nonlinear problem. the 'error' as you call it lies in the fact that your structure has formed a mechanism and, due to your choice of rayleigh coeeficients, the tangent matrix is singular at the point where the analysis fails. You have 2 possible solutions:
1) modify the material parameters a bit so that you don't come to this situation, i.e. they return probably a 0 tangent at the failure point.
2) use some of the initial tangent as part of the damping matrix or some mass proportional damping.
both options effect your model and it's damping characteristics. if either option works you will get an answer, the correctness of which you must determine.

[SunnyPark]

Thank you, fmk!!
i want to choose 2nd slolution, and i've got question....
what is "the initial tangent" mean??
is that mean one of the rayleigh damping constants??
thank you again.

[SunnyPark]

I....still cannot solve this problem... T^T
help me plz........