a simple example

For developers writing C++, Fortran, Java, code who have questions or comments to make.

Moderators: silvia, selimgunay, Moderators

Post Reply
tangzhen
Posts: 30
Joined: Thu May 30, 2013 4:58 am
Location: HIT

a simple example

Post by tangzhen » Wed Oct 09, 2013 7:11 pm

i used the changed opensees in which i added a new element and a new integrator,"trussforenergy" and "enegryconserving" respectively, then i got the failed information below. i am wondering how to deal with this problem.

WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 2
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 0.01
OpenSees > analyze failed, returned: -3 error flag
-3

the example code :
wipe
model BasicBuilder -ndm 2 -ndf 2
node 1 0.0 0.0
node 2 144.0 0.0
node 3 168.0 0.0
node 4 72.0 96.0
fix 1 1 1
fix 2 1 1
fix 3 1 1
uniaxialMaterial Elastic 1 3000
element trussforenergy 1 1 4 10.0 1
element trussforenergy 2 2 4 5.0 1
element trussforenergy 3 3 4 5.0 1
pattern Plain 1 "Linear" {
load 4 100 -50
}
system BandGeneral
numberer RCM
constraints Plain
integrator EnergyConserving
test EnergyIncr 1.0e-12 6 0
algorithm Newton
analysis Transient
analyze 100 0.01

fmk
Site Admin
Posts: 5883
Joined: Fri Jun 11, 2004 2:33 pm
Location: UC Berkeley
Contact:

Re: a simple example

Post by fmk » Thu Oct 10, 2013 8:12 am

the solver failed .. the tangent is not correct in the element, the material, or both. in your code use opserr << theTangent or *theTangent or whatever to spit out the tangent you are returing from the element or material to find the problem. put this just BEFORE the return statement.

Post Reply