← Back to team overview

yade-users team mailing list archive

Re: [Question #682859]: how to calculate anisotropy of contact normal vector

 

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

    Status: Open => Answered

William Chevremont proposed the following answer:
Hi,

PDFEngine do that for you. It compute the average of quantities like
stress, distance, and probability of interaction for each orientation.
For now, it works only for spherical particles.

Just add 
PDFEngine(filename="PDF.txt", iterPeriod=1000, numDiscretizeAnglePhi=9,numDiscretizeAngleTheta=13, label="PDFRecorder")

in your engines. This will output a PDF.txt containing the density
functions, with column name like NormalStress_xy(theta, phi)

Formulas: 
x = r sin theta sin phi
y = r cos phi
z = r sin theta cos phi

This is a work in progress (not well documented for now) but it work as
is.

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