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