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?