algorithm failed, analyze failed..

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

Moderators: silvia, selimgunay, Moderators

Post Reply
melis
Posts: 18
Joined: Fri Dec 31, 2010 6:55 am
Location: Kocaeli Universty

algorithm failed, analyze failed..

Post by melis »

hi;
i have some problem about algorithm and analyze. i modeled cantilever column. i took element properties from "column database"(http://www.ce.washington.edu/~peera1/proprect.htm) my code is failed and my results is different from test results. what is my code problem? i will send my code and please say me what is the solution?

#http://www.ce.washington.edu/~peera1/proprect.htm
#105 Saatcioglu and Ozcebe 1989, U3
#Units; N,m
wipe
model BasicBuilder -ndm 2 -ndf 3
file mkdir DataU3
node 1 0 0 0
node 2 0 0.5 0
node 3 0 1 0
fix 1 1 1 1
set cb 350e-3
set cby [expr $cb/2]
set ch 350e-3
set chy [expr $ch/2]
set cpas 25e-3
set y1 [expr ($cby-$cpas)]
set z1 [expr ($chy-$cpas)]
set A25 [expr (3.14159*25e-3*25e-3)/4]
set detr 10e-3
uniaxialMaterial ConfinedConcrete01 1 S1 -34.8e6 33.e9 -epscu -0.003 -nu 0.2 1.18 10e-3 75.e-3 470.e6 200e9 0.0 1000. 25e-3 -stRatio 0.0321
#uniaxialMaterial Concrete01 2 -20e6 -0.002 -5e6 -0.003
uniaxialMaterial Steel01 3 430e6 200e9 0.0
section Fiber 1 {
patch quad 1 5 5 -$y1 -$z1 $y1 -$z1 $y1 $z1 -$y1 $z1 #1
#patch quad 1 5 5 -$cby -$chy -$y1 -$chy -$y1 $chy -$cby $chy #2
#patch quad 1 5 5 $y1 -$chy $cby -$chy $cby $chy $y1 $chy #3
#patch quad 1 5 5 -$y1 $z1 $y1 $z1 $y1 $chy -$y1 $chy #4
#patch quad 1 5 5 -$y1 -$chy $y1 -$chy $y1 -$z1 -$y1 -$z1 #5
layer straight 3 3 $A25 [expr -($y1-22.5e-3)] [expr -($z1-22.5e-3)] [expr -($y1-22.5e-3)] [expr ($z1-22.5e-3)]
layer straight 3 3 $A25 [expr ($y1-22.5e-3)] [expr -($z1-22.5e-3)] [expr ($y1-22.5e-3)] [expr ($z1-22.5e-3)]
layer straight 3 2 $A25 0 [expr -($z1-22.5e-3)] 0 [expr ($z1-22.5e-3)]
}
puts "section is defined"
geomTransf PDelta 1
set np 5
element dispBeamColumn 1 1 2 $np 1 1
element dispBeamColumn 2 2 3 $np 1 1
recorder Node -file DataU3/disp2.txt -node 3 -dof 1 disp
recorder Node -file DataU3/force1.txt -node 1 -dof 1 reaction
pattern Plain 1 Constant {
load 3 0. -600e3 0.0
}
system BandSPD
constraints Plain
numberer RCM
algorithm Linear
integrator LoadControl 1
analysis Static
analyze 1
loadConst -time 0.0
wipeAnalysis
########################################
set DtAnalysis 0.02
set TmaxAnalysis 20
set maxNumIter 10
set Tol 1.e-8
set TestType NormDispIncr;
set printFlag 0;
set force "Path -filePath forceu3.txt -dt 0.02"
pattern Plain 2 $force {
load 3 1 0.0 0.0
}
system BandGeneral
constraints Plain
numberer Plain
test $TestType $Tol $maxNumIter $printFlag
set algorithmType ModifiedNewton
algorithm $algorithmType
integrator Newmark 0.5 0.25
analysis Transient
set Nsteps [expr int($TmaxAnalysis/$DtAnalysis)];
set ok [analyze $Nsteps $DtAnalysis]; # actually perform analysis; returns ok=0 if analysis was successful

if {$ok != 0} { ; # if analysis was not successful.
# change some analysis parameters to achieve convergence
# performance is slower inside this loop
# Time-controlled analysis
set ok 0;
set controlTime [getTime];
while {$controlTime < $TmaxAnalysis && $ok == 0} {
set ok [analyze 1 $DtAnalysis]
set controlTime [getTime]
set ok [analyze 1 $DtAnalysis]
if {$ok != 0} {
puts "Trying Newton with Initial Tangent .."
test NormDispIncr $Tol 1000 0
algorithm Newton -initial
set ok [analyze 1 $DtAnalysis]
test $TestType $Tol $maxNumIter 0
algorithm $algorithmType
}
if {$ok != 0} {
puts "Trying Broyden .."
algorithm Broyden 8
set ok [analyze 1 $DtAnalysis]
algorithm $algorithmType
}
if {$ok != 0} {
puts "Trying NewtonWithLineSearch .."
algorithm NewtonLineSearch .8
set ok [analyze 1 $DtAnalysis]
algorithm $algorithmType
}
}
}; # end if ok !0


puts "force. End Time: [getTime]"
wipe

thank you...
vesna
Posts: 3033
Joined: Tue May 23, 2006 11:23 am
Location: UC Berkeley

Re: algorithm failed, analyze failed..

Post by vesna »

Show me the OpenSees messages after you run the code.
melis
Posts: 18
Joined: Fri Dec 31, 2010 6:55 am
Location: Kocaeli Universty

Re: algorithm failed, analyze failed..

Post by melis »

OpenSees > analyze failed, returned: -3 error flag
Trying Broyden ..
WARNING: CTestNormDispIncr::test() - failed to converge
after: 10 iterations
ModifiedNewton::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 4.32
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestNormDispIncr::test() - failed to converge
after: 1000 iterations
NewtnRaphson::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 4.32
OpenSees > analyze failed, returned: -3 error flag
Trying Broyden ..
WARNING: CTestNormDispIncr::test() - failed to converge
after: 10 iterations
WARNING: CTestNormDispIncr::test() - failed to converge
after: 10 iterations
WARNING BandGenLinLapackSolver::solve() -LAPACK routine returned 6
WARNING Broyden::solveCurrentStep() -the LinearSysOfEqn failed in solve()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 4.32
OpenSees > analyze failed, returned: -3 error flag
Trying NewtonWithLineSearch ..
WARNING: CTestNormDispIncr::test() - failed to converge
after: 10 iterations
NewtonLineSearch::solveCurrentStep() -the ConvergenceTest object failed in test(
)
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 4.32
OpenSees > analyze failed, returned: -3 error flag
force. End Time: 4.300000
OpenSees >
vesna
Posts: 3033
Joined: Tue May 23, 2006 11:23 am
Location: UC Berkeley

Re: algorithm failed, analyze failed..

Post by vesna »

Try using Newton as a basic algorithm (instead of Modified Newton). Try variable transient instead of transient analysis.
melis
Posts: 18
Joined: Fri Dec 31, 2010 6:55 am
Location: Kocaeli Universty

Re: algorithm failed, analyze failed..

Post by melis »

hi Vesna;
Thank you so much for your quick answer. i tried your suggestion, but it did'nt work. can you say me another solution, because it is really important for me.
do you think my code is wrong or not? what is the main problem?
i wait your reply,thank you..
vesna
Posts: 3033
Joined: Tue May 23, 2006 11:23 am
Location: UC Berkeley

Re: algorithm failed, analyze failed..

Post by vesna »

the fact that your analysis is stopping in the course of the analysis indicates stability of your model. since I do not know much about what you are trying to model I can not say anything about accuracy of your model. I noticed that you are modeling your column with only 2 displacement based elements (DBE). I am afraid that this is not sufficient. I would not model a column with less then 8 DBEs.

Can you tell me what is the second type of analysis that you are performing?
melis
Posts: 18
Joined: Fri Dec 31, 2010 6:55 am
Location: Kocaeli Universty

Re: algorithm failed, analyze failed..

Post by melis »

hi Vesna;
The second type of analysis is dynamic analysis. (105.Saatcioglu and Ozcebe 1989, U3) experimental model has force-displacement test result. i don't know loading type,because test properties table don't have anything about loading procedure(it has only axial load 600kn), but i think this way; if i apply force data (this data is test result file) at free node, i can take displacement result, so i downloaded force.txt file.
then i applied this load like time series .
"set force "Path -filePath forceu3.txt -dt 0.02"
pattern Plain 2 $force {
load 3 1 0.0 0.0
}"
is it true?
if it isn't, what is the correct analysis type?
thank you..
vesna
Posts: 3033
Joined: Tue May 23, 2006 11:23 am
Location: UC Berkeley

Re: algorithm failed, analyze failed..

Post by vesna »

If you expect nonlinear behavior of your column this would not be a right way to define the load for your column. Why don't you try to find ground motion that they used in their test and use that one for your analysis.
melis
Posts: 18
Joined: Fri Dec 31, 2010 6:55 am
Location: Kocaeli Universty

Re: algorithm failed, analyze failed..

Post by melis »

hi Vesna;
after your reply, i checked "http://www.ce.washington.edu/~peera1/" this site, and i didn't find any knowledge abaout loading type. then i wrote "Mike Berry and Marc Eberhard at the University of Washington". they replied me, they said; "The main result in the data file are indeed the applied forces and displacements. In general, the displacements are imposed on the column, and the forces needed to impose those displacements are measured (i.e., displacement control). There is no dynamic effect, since the displacements are applied slowly". i decided to apply force.txt file, and static analysis. i changed second type of analysis, i wrote new code usin tcl command.
###########################################################
set y [open "forceu3.txt" r] #forceu3.txt is data file
for {set j 1} {$j<=1010} {incr j 1} {
gets $y v
set k($j) $v
}
puts "1 $k(1)"
pattern Plain 2 Constant {
load 3 $k(1) 0. 0.
}
pattern Plain 3 Constant {
for {set i 2} {$i<250} {incr i 1} {
set n [expr $i-1]
set m($i) [expr $k($i)-$k($n)]
puts "$i $m($i) "
load 3 $m($i) 0. 0.
system BandSPD
constraints Plain
numberer RCM
test NormDispIncr 1.0e-8 10
algorithm Newton
integrator LoadControl 1
analysis Static
analyze 1
}
}
wipe
##########################
then i took error message,after read number of data 217;
##########################
.
.
216 1800.0
217 4300.0
218 6100.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 6
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
219 4700.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 6
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
220 4800.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 6
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
221 5200.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 6
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
222 4800.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 6
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
223 2700.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 6
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
224 6300.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 6
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
225 2400.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 6
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
226 3500.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 6
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
227 6400.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 5
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
228 9900.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 5
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
229 3100.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 5
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
230 8200.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 5
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
231 0.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 5
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
232 1300.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 5
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
233 -9300.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 5
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
234 -9800.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 5
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
235 -1800.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 5
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
236 -9900.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 6
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
237 -9600.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 6
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
238 -9400.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 6
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
239 -9700.0
WARNING BandSPDLinLapackSolver::solve() - the LAPACK routines returned 6
WARNING NewtonRaphson::solveCurrentStep() -the LinearSysOfEqn failed in solve()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 217
OpenSees > analyze failed, returned: -3 error flag
240 -10700.0
241 -10000.0
242 -10200.0
243 -10300.0
244 -10700.0
245 -10800.0
246 -10500.0
247 -10800.0
248 -9810.0
249 -10670
####################################
i checked two data file, first is applied "forceu3",second is support reaction. both data file are same to number of data 217. then support reaction changed. what is the problem? what is the meanin of this error message? what is your suggestion? i'm waiting your reply Vesna.
Thank you so much,i'm grateful..
vesna
Posts: 3033
Joined: Tue May 23, 2006 11:23 am
Location: UC Berkeley

Re: algorithm failed, analyze failed..

Post by vesna »

If you want to simulate cyclic test you need to use displacement vector and assign it to time series. You can not go with forces. On the message board there are several topics that cover this issue.
Post Reply