Question #701861 on Yade changed:
https://answers.launchpad.net/yade/+question/701861
Description changed to:
In my case, the consolidation is controlled by the def func():
However, when I checked the mean stress p = (sigmax+sigmay+sigmaz)/3 it
did not steadily increase with iterPeriod . So my function doesn't work
as my assumption in my simulation. But I can run the whole script with
no error.
>From terminal, I can get:
fabs(getStress()[0,0]+getStress()[1,1]+getStress()[2,2])/3.0
Out[75]: 6.432893541478686
In [76]: fabs(getStress()[0,0]+getStress()[1,1]+getStress()[2,2])/3.0
Out[76]: 11.011277803935808
In [77]: fabs(getStress()[0,0]+getStress()[1,1]+getStress()[2,2])/3.0
Out[77]: 19.46480982572519
In [78]: fabs(getStress()[0,0]+getStress()[1,1]+getStress()[2,2])/3.0
Out[78]: 18.54078400999076
In [79]: fabs(getStress()[0,0]+getStress()[1,1]+getStress()[2,2])/3.0
Out[79]: 15.016490138434548
In [80]: fabs(getStress()[0,0]+getStress()[1,1]+getStress()[2,2])/3.0
Out[80]: 6.468796221436976
In [81]: fabs(getStress()[0,0]+getStress()[1,1]+getStress()[2,2])/3.0
Out[81]: 17.7652627151387
The mean stress is jumping...
#limited consolidation pressure
finalconso_p=35.0
flag=0
def checkStress():
global flag
if
fabs(getStress()[0,0]+getStress()[1,1]+getStress()[2,2])/3.0>=finalconso_p:
fabric()
addPlotData()
print("mean stress
p=",fabs(getStress()[0,0]+getStress()[1,1]+getStress()[2,2])/3.0)
Finished()
--
You received this question notification because your team yade-users is
an answer contact for Yade.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
