yade-users team mailing list archive

[enki <question701861@xxxxxxxxxxxxxxxxxxxxx>]

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.