[DEMApplication] New linear contact force model #13008
Conversation
rlrangel
requested review from
AlejandroCornejo and
ChengshunShang1996
and removed request for
a team
|
one question, don't you add any test for the future? |
| const double my_poisson = element->GetPoisson(); | ||
|
|
||
| mKn = my_young * my_radius; // normal | ||
| mKt = my_poisson * mKn; // tangent |
There was a problem hiding this comment.
Hi Rafael. In this line, you are calculating the tangential stiffness using the "poisson_ratio" from the particle. However, in both of your papers, you were using the item "tangential-to-normal stiffness ratio". To me, the physical significance of those two items is different.
There was a problem hiding this comment.
Zhao actually calls it the Poison ration. I decided to call it tangential-to-normal stiffness ratio in my papers exactly because it has a different physical meaning.
| // References using this model: | ||
| // R.L. Rangel et al. (2024). Multiscale data-driven modeling of the thermomechanical behavior of granular media with thermal expansion effects. Comput Geotech, 176:106789. | ||
| // R.L. Rangel et al. (2024). A continuum--discrete multiscale methodology using machine learning for thermal analysis of granular media. Comput Geotech, 168:106118. | ||
| // S. Zhao et al. (2020). Multiscale modeling of thermo-mechanical responses of granular materials: A hierarchical continuum--discrete coupling approach. CMAME, 367:113100. |
There was a problem hiding this comment.
In the paper by Zhao, although they said that the simple linear contact model was used, they did not explicitly give the formulation about how to calculate the normal and tangential stiffness. Should this one be replaced with another piece of literature?
There was a problem hiding this comment.
That's true. I changed the reference.
| CalculateTangentialForceWithNeighbour(normal_contact_force, OldLocalElasticContactForce, LocalElasticContactForce, LocalDeltDisp, sliding, element, wall); | ||
| } | ||
|
|
||
| template<class NeighbourClassType> |
There was a problem hiding this comment.
A template function is suggested to be put in the head file if I am right. Please look at this link.
There was a problem hiding this comment.
That's true. Done it.
I decided to make another PR later exclusively dedicated to adding tests. |
Adding a new simple linear force model, which I have recently used in a research.