← Back to team overview

yade-users team mailing list archive

Re: About the moment calculation in cohesion-friction model



Thanks for your warm help. I have tried this model. However, the new problem is coming. You know, the CFphPhys just can be matched to Dem3DofGeom. There are no the Ig2_Box_Sphere_Dem3DofGeom. If I must choose the Ig2_Facet_Sphere_Dem3DofGeom instead of Ig2_Box_Sphere_Dem3DofGeom. then the ThreeDTriaxialEngine can not approve the Facet as the boundary. If I want to change Facet instead of the box in ThreeDTriaxialEngine, I cannot give the only one Facet id to impose the strain on it. You know, The rectangle should be composed of two facet.
So, how can I deal with it? Thank you!


At 2010-09-17 10:49:31,"luc scholtes" <lscholtes63@xxxxxxxxx> wrote:
Try to have a look at CohesiveFrictionalPM with CFpmPhys and CFpmMat which allow cohesion, tensile strength and moment at contact.



2010/9/17 焦利青<jiaoliqing@xxxxxxx>

Yes, you are all right. I mentioned the cohesion-friction model is Law2_ScGeom_CohFrictPys_ElasticPlastic. I am sorry, I should point out clearly. As your means, If I would like to consider the effect of moment, this one is not suitable, isn't it? I know the code of moment algorithm has been realized in the Cohesionless one. But, the effect of cohesion must be considered in my model. In the conclusion, if both the moment and cohesion should be considered, how can I do it? May I have your help? thank you.


At 2010-09-16 19:32:15,"Jerome Duriez" <duriez@xxxxxxxxxxxxxxx> wrote:

>You can find this plastic threshold for the moment in 
>Law2_ScGeom_NormalInelasticityPhys_NormalInelasticity (lines 193 to 202 
>of the .cpp). For what concerns the "cohesion-friction model", I assume 
>that you speak of Law2_ScGeom_CohFrictPhys_ElasticPlastic, which I know 
>less. But a quick look to the .cpp gives me the same feeling than yours, 
>that the moment is always elastic.
>Le 16/09/2010 12:40, 焦利青 a écrit :
>> Hello,
>> I found the moment calculation in cohesion-friction model is different 
>> from the traditional algorithm which is listed in the 
>> paper(Jean-Patrick Plassiard etl, A spherical discrete element 
>> model:calibration procedure and incremental response, 209; N. Belheine 
>> etl, Numerial simulation of drained triaxial test using 3d discrete 
>> element modeling; Rolling resistance at contacts in simulation of 
>> shear band development by DEM). In the traditional algorithms, the 
>> slider for rolling starts working if the moment exceeds a value of 
>> "eta* Normal force". However, I cannot find the related code described 
>> it. If no any condition to relax the moment. then the moment will act 
>> all the time. I think it should be wrong. Maybe I haven't found the 
>> right file. May I have your help? Thank you.
>> Best
>> Liqing

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

Follow ups