yade-users team mailing list archive

[Sina Jafari <question247136@xxxxxxxxxxxxxxxxxxxxx>]

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

hi all,
I tried to create polar interaction histogram using interaction-histogram.py for my simulation but the problem is that nothing is plotted in the output file. I'm using utils.plotDirections() at the same time and it works fine but I need to plot it myself. Can you please help me with this? thanks in advance.

I've defined a function to plot histogram at certain iterperiods. it is as follows:

def Intrhisto():
	axis=2 # x, i.e. plot the xy plane
	ax1,ax2=(axis+0)%3,(axis+1)%3 ## get the other two indices, i.e. 0 and 1 in this case
	angles,forces=[],[]
	for z in O.interactions:
		if not z.isReal: continue
		norm=z.geom.normal
		if norm[ax1]==0:
			angle=90
			force=z.phys.normalForce.norm()
		else:
			angle=np.arctan(norm[ax2]/norm[ax1])
			force=z.phys.normalForce.norm()
		angles.append(angle)
		forces.append(force)
	# polar histogram
	pylab.figure()
	## prepare data
	values,bins=numpy.histogram(angles,weights=forces,bins=20)
	## prepare polar plot
	pylab.subplot(111,polar=True);
	pylab.plot()
	pylab.savefig('Norm-Force-Histogram.JPEG')

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