yade-users team mailing list archive

[behzad <question262922@xxxxxxxxxxxxxxxxxxxxx>]

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

behzad posted a new comment:
Yes, this value came from PFC!

I changed it to zero. Now the ratio of the applied torque to the obtained angular acceleration is correct.
However, the problem is the calculation of inertia is not correct. 

we have:

upper_plate =
O.bodies.append(box((0,0,0),(1,1,0.1),fixed=True,material=CohMat))

]: O.bodies[upper_plate].state.mass
 ->  [2]: 1600.0


So from

body->state->inertia=Vector3r(mass*(4*extents[1]*extents[1]+4*extents[2]*extents[2])/12.,mass*(4*extents[0]*extents[0]+4*extents[2]*extents[2])/12.,mass*(4*extents[0]*extents[0]+4*extents[1]*extents[1])/12.);

we must have inertia as (538.66,538.66,1066.66)

But we get:

 O.bodies[upper_plate].state.inertia
 ->  [3]: Vector3(8080,8080,16000)

Where does it come from?

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