Hello everyone!
I'm trying to model a earth dam with opensees. To do TCL file, I used GiD ProblemTypes for 2D Slope Analysis Input File Generation.
When I do dynamic analysis with small earthquakes the model works, but with longer earthquakes and higher accelerations I can't reach convergence and I get the following warning:
WARNING: CTestNormDispIncr::test() - failed to converge
after: 200 iterations current Norm: 0.00296852 (max: 0.001, Norm deltaR: 9.76904e+07)
NewtnRaphson::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 0.08125
OpenSees > analyze failed, returned: -3 error flag
Finished with dynamic analysis...
Analysis execution time: 1927 seconds
To solve this problem I tried to change the constitutive law of the material to the base of the dam passing from PressureIndependMultiYield Material to Elastic isotropic, the system gives the following warning:
damping coefficients: a_0 = 0.04976780441980199; a_1 = 0.00031515830311111683
number of steps in analysis: 7990
analysis time step: 0.005
WARNING: MaterialStageParameter::setDomain() - no effect with material tag 1
WARNING: MaterialStageParameter::setDomain() - no effect with material tag 4
CTestNormDispIncr::test() - iteration: 1 current Norm: 1.19313e+06 (max: 0.001, Norm deltaR: 9.38208e+07)
CTestNormDispIncr::test() - iteration: 2 current Norm: 0.42086 (max: 0.001, Norm deltaR: 9.38208e+07)
Is there someone that can help me?
Hi David,
Changing your algorithm to Krylov-Newton is one of the ways of solving the convergence problem.
- Badges
- Science topic
- Similar topics
- Filing