yade-users team mailing list archive

[behzad <question265521@xxxxxxxxxxxxxxxxxxxxx>]

New question #265521 on Yade:
https://answers.launchpad.net/yade/+question/265521

Hi,

I have a simulation in which a sinusoidal strain is applied to a cylindrical sample of spheres by a box. I'd like to get the imposed stress history which is supposed to be a sinusoidal curve with a time lag. 
The history of stress on the box with bodyStressTensors()[box_1][2][2] has some noise.
I decided to get the average stress of the spheres in the sample in z direction as the imposed stress. 

This can be done via calling the following function at certain cycles by Pyrunner:

def get2():
	stress=0.0
	sums=0.0
	global a
	num=len(a)
	for i in a:
		sums+=bodyStressTensors()[i][2][2]
	stress=sums/num
	return stress

However, this is not efficient at all (as you can tell, with such a loop in the function). 
So, what's the most efficient way to get the stress in z direction from spherical bodies of a model with iterPeriod of e.g.100.

Shall I use VTKRecorder? 

I wonder if others are using VTKRecorder, that much or not? But, I'm not used to it yet. I appreciate if you can share your ideas or experience in dealing with such a case.

Regards,
Behzad 

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.