# How to connect the material parameters with contact stiffness for polyhedra?

Hey

#########################################
Real Kna = mat1->young;
Real Knb = mat2->young;
Real Ksa = mat1->young*mat1->poisson;
Real Ksb = mat2->young*mat2->poisson;
Real frictionAngle = math::min(mat1->frictionAngle,mat2->frictionAngle);
contactPhysics->tangensOfFrictionAngle = math::tan(frictionAngle);
contactPhysics->kn = Kna*Knb/(Kna+Knb);
contactPhysics->ks = Ksa*Ksb/(Ksa+Ksb);
#########################################

I am wondering why the value of young could be assigned to the Kn directly? I am speculating that the name "young" does not have any specific meaning here,just represents the stiffness. Is that right?
If so, is there any method implemented in Yade actually connects the stiffness with the material parameters(such as E,shear modulus G) for a polyhedron? (distinguished from calculating the stiffness using Young's moduluus, diameter and poison's ratio for sphere contacts)

Thanks

Jasper

## Question information

Language:
English Edit question
Status:
For:
Assignee:
No assignee Edit question
Last query:
2020-03-11
2020-03-11
 Jan Stránský (honzik) said on 2020-03-11: #1

Hello,

> I am wondering why the value of young could be assigned to the Kn directly?

why not? :-)

> just represents the stiffness.

yes

> the name "young" does not have any specific meaning here ...

The specific meaning is that it represents the normal stiffness :-)

if you mean that "young" should be Young's modulus, than yes, it may not be the case here.
The name is inherited from FrictMat and has the meaning of stiffness in normal direction, not necessarily Young's modulus.

have a look at a similar discussion [1]

> If so, is there any method implemented in Yade actually connects the stiffness with the material parameters(such as E,shear modulus G) for a polyhedron? (distinguished from calculating the stiffness using Young's moduluus, diameter and poison's ratio for sphere contacts)

as usually in DEM, you have to calibrate the parameters to obtain desired response (if I got the paragraph correctly)

cheers
Jan

 Jérôme Duriez (jduriez) said on 2020-03-11: #2

If you read source code, [*] may be relevant as well.

You will see that PolyhedraMat.young and PolyhedraPhys.kn may even have different "units" than usual (and possibly variable ones, depending on the Law2*.volumePower)

 Chareyre (bruno-chareyre-9) said on 2020-03-11: #3

> is there any method implemented in Yade actually connects the stiffness
with the material parameters(such as E,shear modulus G) for a polyhedron?

Hi,
Hertzian theory of contact is the answer to your question in the case of
spheres. And (surprise surprise) it leads to non linear elasticity and
vanishing stiffness in unconfined materials. So even for spheres there is
no answer to your question (no single stiffness value due to non-linearity).
There is no established theory for polyhedral particles AFAIK.
Cheers
Bruno

On Wed, 11 Mar 2020 at 14:09, Jan Stránský <

> Question #689263 on Yade changed:
>
>
> Jan Stránský proposed the following answer:
> Hello,
>
> > I am wondering why the value of young could be assigned to the Kn
> directly?
>
> why not? :-)
>
> > just represents the stiffness.
>
> yes
>
> > the name "young" does not have any specific meaning here ...
>
> The specific meaning is that it represents the normal stiffness :-)
>
> if you mean that "young" should be Young's modulus, than yes, it may not
> be the case here.
> The name is inherited from FrictMat and has the meaning of stiffness in
> normal direction, not necessarily Young's modulus.
>
> have a look at a similar discussion [1]
>
> > If so, is there any method implemented in Yade actually connects the
> stiffness with the material parameters(such as E,shear modulus G) for a
> polyhedron? (distinguished from calculating the stiffness using Young's
> moduluus, diameter and poison's ratio for sphere contacts)
>
> as usually in DEM, you have to calibrate the parameters to obtain
> desired response (if I got the paragraph correctly)
>
> cheers
> Jan
>
>
> --
>
> _______________________________________________
> Post to : <email address hidden>
>

--
--
_______________
Bruno Chareyre
Associate Professor
ENSE³ - Grenoble INP
Lab. 3SR
BP 53
38041 Grenoble cedex 9
Tél : +33 4 56 52 86 21
________________

Email too brief?
Here's why: email charter