yade-users team mailing list archive

[Bokkisa Srinivas Vivek <question680398@xxxxxxxxxxxxxxxxxxxxx>]

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

Summary changed to:
Inertia Tensor for clumps with Intersecting Spheres

Description changed to:
Hello community,

I was going through the Source code- Clump.cpp

The inertia tensor at local level for the intersecting spheres clump is
calculated from the below equation (line no. 178):

Ig += m*(
x.dot(x)*Matrix3r::Identity()-x*x.transpose())+Matrix3r(Vector3r::Constant(dv*pow(dx,2)/6.).asDiagonal())

In the second part-
"Matrix3r(Vector3r::Constant(dv*pow(dx,2)/6.).asDiagonal())",

I believe it should be m*pow(dx,2)/6 but instead it computes
dv*pow(dx,2)/6, was it a mistake in code or done purposefully?

Reference:
https://github.com/yade/trunk/blob/master/core/Clump.cpp#L178

Thank you

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