Pushover analysis

Forum for OpenSees users to post questions, comments, etc. on the use of the OpenSees interpreter, OpenSees.exe

Moderators: silvia, selimgunay, Moderators

Libra23
Posts: 96
Joined: Tue Nov 07, 2006 9:21 am
Location: University

Pushover analysis

Post by Libra23 »

after running eigen analyis I started with gravity but it gives me this problem :

WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 178
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() -the algorithm failed at iteration :0 with domain at load factor 0.2
analyze failed,returned -3 error flag

so I have also problems with pushover analysis,
if I send you my file colud you look at it?
Thanks
silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia »

it's a good sign, you have convergence problems.
you need to change something about the analysis at this point.
please look at the examples manual, maybe example 5 or so, it has the script you want.
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
Libra23
Posts: 96
Joined: Tue Nov 07, 2006 9:21 am
Location: University

Post by Libra23 »

I took script from example manual and get that problem,this is the script I use:
# Gravity-analysis parameters -- load-controlled static analysis
set Tol 1.0e-8; # convergence tolerance for test
variable constraintsTypeGravity Plain; # default;
if { [info exists RigidDiaphragm] == 1} {
if {$RigidDiaphragm=="ON"} {
variable constraintsTypeGravity Lagrange; # large model: try Transformation
}; # if rigid diaphragm is on
}; # if rigid diaphragm exists
constraints $constraintsTypeGravity ; # how it handles boundary conditions
numberer RCM; # renumber dof's to minimize band-width (optimization), if you want to
system BandGeneral ; # how to store and solve the system of equations in the analysis (large model: try UmfPack)
test EnergyIncr $Tol 6 ; # determine if convergence has been achieved at the end of an iteration step
algorithm Newton; # use Newton's solution algorithm: updates tangent stiffness at every iteration
set NstepGravity 10; # apply gravity in 10 steps
set DGravity [expr 1./$NstepGravity]; # first load increment;
integrator LoadControl $DGravity; # determine the next time step for an analysis
analysis Static; # define type of analysis static or transient
analyze $NstepGravity; # apply gravity

What do you suggest to change?I tried something but get the same!
Thanks
silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia »

no. that one is for gravity, right?
you want to do the pushover.
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
Libra23
Posts: 96
Joined: Tue Nov 07, 2006 9:21 am
Location: University

Post by Libra23 »

Yes I have to do pushover but the problem I get:

Code: Select all

WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 178 
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve() 
StaticAnalysis::analyze() -the algorithm failed at iteration :0 with domain at load factor 0.2 
analyze failed,returned -3 error flag 
is about gravity analysis and then the same for each step of pushover!

I don't know what's wrong!
Libra23
Posts: 96
Joined: Tue Nov 07, 2006 9:21 am
Location: University

Post by Libra23 »

For pushover I use script of example 8.
Libra23
Posts: 96
Joined: Tue Nov 07, 2006 9:21 am
Location: University

Post by Libra23 »

I need help!
Thanks
silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia »

the pushover script in example 8 should help.
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
Libra23
Posts: 96
Joined: Tue Nov 07, 2006 9:21 am
Location: University

Post by Libra23 »

I looked at it,maybe it could be a model problem...if you remember I have modeled walls as non linear beam column elements with large section ,at the edge of section I had two nodes ( 0 )that are connected with nothing...someone suggests me to use beam with elastic section to connect them ( 0 ) with node ( X ) of non linear beamColumn that represent the wall.
I mean :

0-------X--------0


it doens't work the same!
What do you think about it?
silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia »

i am a bit confused, you might want to discuss this with someone who is physically near you, as it is a modelling issue.
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
Libra23
Posts: 96
Joined: Tue Nov 07, 2006 9:21 am
Location: University

Post by Libra23 »

You are right....it's a pity that nobody near me can't use this program and is able to help!
Thanks the same,really, I believe I must help myself!
Libra23
Posts: 96
Joined: Tue Nov 07, 2006 9:21 am
Location: University

Post by Libra23 »

I have read in the manual that ,using eleLoad command ,nonlinearBeamColumn might have problems converging .... could it be this that gives me problems in my model?
silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia »

nope.
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
Libra23
Posts: 96
Joined: Tue Nov 07, 2006 9:21 am
Location: University

Post by Libra23 »

in examples 8 of manual if I want to run pushover do I open example about pushover and call within it example that defines model?
silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia »

first you source the build file,
then you source the analysis file
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
Post Reply