yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #25208
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.