yade-users team mailing list archive

[Alexander <question271394@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Answered => Open

Alexander is still having a problem:
Hello, Luc.  Thank you for your response.

So now i have 3 questions:

1) Do i understand correct that in the incremental form we have
following computation formulas for moments:

dO = (w_1 - w_2)dt - is the rotational vector during the time increment
dt between two simulation steps, w_i - are rotational velocuties at time
t_0 + dt

dO_t =  (d_n*n)n - is twisting component, which is the projection of
dOmega to the normal of the contact at time t_0 + dt

dO_b = dO - dO_n - is the bending component, which lies in the contact
plane:

twisting increment dM_t = k_t * dO_t - total twisting at time t+dt: M_t
= M_t + dM_t

bending increment dM_b= k_b * dO_b - total bending at time t+dt: M_b =
M_b+ dM_b

2) How these formulas look like for non-uncremental form. Is it just simple replacing  M with dM, i mean at step t+dt
we just compute M_t and M_b directly?

M_t = k_t * dO_t  - twisting

M_t = k_t * dO_t - bending

3) If  Law2_ScGeom6D_CohFrictPhys_CohesionMoment::momentRotationLaw
=False then there is now angular velocities (w_i=0) of particles and
shear displacement is computed without part of angular velocitiues (only
displacement of sphere centers)?

with regards, Alexander

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