Re: [Yade-users] [Question #683001]: Physic interaction betweem sphere and polyhedra

[Jan Stránský]

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

Status: Open => Answered

Jan Stránský proposed the following answer:
> What puzzles me is that polyhedra do not have radius r, then how to
calculate normal stiffness between polyhedra and sphere

the same way as it does for walls, facets... just considers that the
other body has some artificial radius related to the sphere's radius
(e.g. equal) [1,2,3]

cheers
Jan

[1] https://gitlab.com/yade-dev/trunk/blob/master/pkg/dem/FrictPhys.cpp#L19
[2] https://gitlab.com/yade-dev/trunk/blob/master/pkg/dem/Polyhedra_Ig2.cpp#L419
[3] https://gitlab.com/yade-dev/trunk/blob/master/pkg/dem/Polyhedra_Ig2.cpp#L335

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

___
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help   : https://help.launchpad.net/ListHelp

[Yade-users] [Question #683001]: Physic interaction betweem sphere and polyhedra

[lip]

New question #683001 on Yade:
https://answers.launchpad.net/yade/+question/683001

hello
 i know that Ig2_Sphere_Polyhedra_ScGeom can calculate the ig between 
polyhedra and sphere, if the material of polyhedra is polyhedraMat and sphere 
is frictMat, the ip between them can be calculate by 
Ip2_FrictMat_PolyhedraMat_FrictPhys, and then use cundall law, the interaction 
between polyhedra and sphere can be solved.
 Ip2_FrictMat_PolyhedraMat_FrictPhys  describe in Polyhedra.cpp is as below:
void Ip2_FrictMat_PolyhedraMat_FrictPhys::go(const shared_ptr& pp1, 
const shared_ptr& pp2, const shared_ptr& interaction){
const shared_ptr& mat1 = YADE_PTR_CAST(pp1);
const shared_ptr& mat2 = YADE_PTR_CAST(pp2);
Ip2_FrictMat_FrictMat_FrictPhys().go(mat1,mat2,interaction);
}

and stiffness culculate in Ip2_FrictMat_FrictMat_FrictPhys() is like that:
Real Kn = 2*Ea*Ra*Eb*Rb/(Ea*Ra+Eb*Rb);
Real Ks = 2*Ea*Ra*Va*Eb*Rb*Vb/(Ea*Ra*Va+Eb*Rb*Vb);

What puzzles me is that polyhedra do not have radius r, then how to calculate 
normal stiffness between polyhedra and sphere use  
Ip2_FrictMat_PolyhedraMat_FrictPhys?
lookforward for your answer.

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

___
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help   : https://help.launchpad.net/ListHelp