Implementing Bounding Cam Clay

Forum for asking and answering questions related to use of the OpenSeesPy module

Moderators: silvia, selimgunay, Moderators

Post Reply
Trish
Posts: 7
Joined: Thu Feb 24, 2022 1:44 am

Implementing Bounding Cam Clay

Post by Trish » Mon May 16, 2022 2:02 am

Hello,
I am trying to implement nDMaterial('BoundingCamClay'). I am getting an error as follows:
CTestNormDispIncr::test() - iteration: 1 current Norm: 64.3004 (max: 1e-06, Norm deltaR: nan)
CTestNormDispIncr::test() - iteration: 2 current Norm: nan (max: 1e-06, Norm deltaR: nan)
.......
WARNING: CTestNormDispIncr::test() - failed to converge
after: 20 iterations current Norm: nan (max: 1e-06, Norm deltaR: nan)
NewtonRaphson::solveCurrentStep() -the ConvergenceTest object failed in test()
DirectIntegrationAnalysis::analyze() - the Algorithm failed at time 100
OpenSees > analyze failed, returned: -3 error flag
My code is

Code: Select all

from openseespy.opensees import *
wipe() 
model('basic','-ndm',3,'-ndf',4)
node(1,1,0,0)
node(2,1,1,0)
node(3,0,1,0)
node(4,0,0,0)
node(5,1,0,1)
node(6,1,1,1)
node(7,0,1,1)
node(8,0,0,1)

fix(1,1,1,1,1)
fix(2,1,1,1,1)
fix(3,1,1,1,1)
fix(4,1,1,1,1)
fix(5,1,1,0,1)
fix(6,1,1,0,1)
fix(7,1,1,0,1)
fix(8,1,1,0,1)

nDMaterial('BoundingCamClay',1,1.8,1,4e3,1,5.4e4,0,0.13,0,1.5)

element('SSPbrickUP',1,1,2,3,4,5,6,7,8,1,2.2e6,1,1e-9,1e-9,1e-9,0.7,6e-5)

test('NormDispIncr',1e-6,20,1)
algorithm('Newton')
analysis('Transient')

t = 60*60
timeSeries('Path',1,'-time',0,24*t,48*t,72*t,96*t,120*t,144*t,'-values',0,10,20,40,80,20,5,'-factor',1)

p = 1/4
pattern('Plain',1,1)
load(5,0,0,-p,0)
load(6,0,0,-p,0)
load(7,0,0,-p,0)
load(8,0,0,-p,0)

analyze(5184,100)
Can any one make me realize, what is the mistake.The analysis is not converging from second step.
The model is working fine with nDMaterial-ElastoIsotropic, Manzari Dafalias and DruckerPrager. With Bounding Cam Clay I am facing this issue.

Thank you.

Post Reply