← Back to team overview

yade-users team mailing list archive

Re: [Question #696037]: contact forces in VTKrecorder

 

Question #696037 on Yade changed:
https://answers.launchpad.net/yade/+question/696037

    Status: Answered => Open

Eleni Gerolymatou is still having a problem:
Thanks for the quick reply Jan and sorry for the lack of an MWE, I missed the instructions. 
I am only partly interested in the visualization and was hoping to avoid using the exporter and the recorder, but of course that is fine. What I need is a list of the contacts with contact force and torque vectors and the contacting particles. Am I right that recorder cannot provide this and I should use exporter instead? Where can I find what 'dictionary' includes?

As a workaround I found that
namefor='force')+'.txt'
a=open(namefor,'w')
for b in O.interactions:
	a.write(str(b.id1)+' '+str(b.id2)+' '+str(b.phys.normalForce[0])+' '+str(b.phys.normalForce[1])+' '+str(b.phys.normalForce[2])+' '+str(b.phys.shearForce[0])+' '+str(b.phys.shearForce[1])+' '+str(b.phys.shearForce[2])+' '+str(b.phys.moment_bending[0])+' '+str(b.phys.moment_bending[1])+' '+str(b.phys.moment_bending[2])+' '+str(b.phys.moment_twist[0])+' '+str(b.phys.moment_twist[1])+' '+str(b.phys.moment_twist[1])+'+\n')
a.close()  

works, but maybe this is not optimal?

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.