Software error in update material function (updateMaterialStage)

[francojj]

Cheers,
I'm doing a Tcl to Py model of a site effect column but during static analysis updating the elastic material to plastic produces an error in the program:

WARNING: CTestNormDispIncr::test() - failed to converge
after: 30 iterations current Norm: 3.54162e-05 (max: 1e-05, Norm deltaR: 1744.74)
AcceleratedNewton::solveCurrentStep() -The ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 10000.1
OpenSees > analyze failed, returned: -3 error flag
Error de convergencia en análisis estático de modeloPerfil

Code: Select all

              
#----------------------------------------------------------------------------------------
#  7. GRAVITY ANALYSIS
#-----------------------------------------------------------------------------------------
        ops.model('basic', '-ndm', 3, '-ndf', 4)
        for k in range(1, capas+1):
            ops.updateMaterialStage('-material', k, '-stage', 0)
               
        print('Inicio de análisis estático elástico\n\n')       
        ops.start()
        ops.analyze(20, 5.0e2)    
        print('Fin de análisis estático elástico\n\n')
        

        
        
        #WHEN TRANSFORMING FROM ELASTIC TO PLASTIC HERE THE ERROR OCCURS:
        # update materials to consider plastic behavior
        for k in range(1, capas+1):
            ops.updateMaterialStage('-material', k, '-stage', 1)
           
           
           
           
           
           
        # plastic gravity loading
        print('Inicio de análisis estático plástico\n\n')
        ok=ops.analyze(40, 5e-2)         
        if ok != 0:
            error = 'Error de convergencia en análisis estático de modelo' + str(perfil) + '\n\n'
            print(error)
            break
        print('Fin de análisis estático plástico\n\n')    

[mhscott]

Can you provide a more complete minimum working example? The code you posted does not define any materials.

[francojj]

Thank you very much for your answer, below I send you the source code (OpenseesPy) the error exists in line # 464.
https://colab.research.google.com/drive ... niqifier=1



I also upload the code that I am transforming into TCL format:

https://colab.research.google.com/drive ... sp=sharing


REGARDS

[francojj]

The results in the design spectrum without executing the material update command exactly match me from the TCL source code with PY.

Code: Select all

ops.updateMaterialStage('-material', k, '-stage', 0)

While when I activate the material update the results vary in the design spectrum

Code: Select all

ops.updateMaterialStage('-material', k, '-stage', 1)