[Yade-users] [Question #428045]: introduction
New question #428045 on Yade: https://answers.launchpad.net/yade/+question/428045 I am new for this DEM in YADE. so I need basic tutorial and theory about DEM and how to simulate in YADE. -- 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 #428044]: introduction
New question #428044 on Yade: https://answers.launchpad.net/yade/+question/428044 I am new for this DEM in YADE. so I need basic tutorial and theory about DEM and how to simulate in YADE. -- 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
Re: [Yade-users] [Question #428040]: How to use SRM method in order to promote failure in YADE?
Question #428040 on Yade changed: https://answers.launchpad.net/yade/+question/428040 Status: Open => Answered Jérôme Duriez proposed the following answer: Hello, The method is exactly the same as for any other modification of existing interaction parameters in YADE. The corresponding interaction physics parameters are here FnMax [1] and FsMax [2]. If you want to reduce YADE-computed values of these parameters for existing interactions, you have to perform a "manual" loop: for cont in O.interactions: cont.phys.FnMax *= 0.9 # or any other value you would like cont.phys.FsMax *= 0.9 # see Python basics for "*=" meaning [1] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.JCFpmPhys.FnMax [2] https://yade-dem.org/doc/yade.wrapper.html#yade.wrapper.JCFpmPhys.FsMax -- 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 #428040]: How to use SRM method in order to promote failure in YADE?
New question #428040 on Yade: https://answers.launchpad.net/yade/+question/428040 Dear all, I want to use the strength reduction method in order to lead the slope model to failure in trunk-master/examples/jointedCohesiveFrictionalPM/gravityLoading.py example. as refer to this in paper [1], I know that should activate strength reduction process (i.e. decreasing imultaneously the tensile strength and cohesion of the local interaction links) for reaching to unstable condition. How can I use this for my work? Is there any specified function or code for this? I had very searching about this subject in Yade documentation and question, but I couldn't any answer. [1] Modelling progressive failure in fractured rock masses using a 3D discrete element method -Luc Scholtès, Frédéric-Victor Donze - http://www.sciencedirect.com/science/article/pii/S1365160912000391 Thanks a lot for attentions, Alma -- 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
Re: [Yade-users] [Question #406523]: Cylinders and law2_ScGeom_ViscElPhys_Basic
Question #406523 on Yade changed: https://answers.launchpad.net/yade/+question/406523 Bruno Chareyre posted a new comment: Hi Raphael, Quick answer to your last point: it is not a problem to "evaluate the kn, ks, cn and cs from the young, poisson, en and es" between spheres and cylinders since by design (inheritance) the cylinders _are_ spheres. So the problem must be simply that the cylinder has a null mass or something, not that it is not a sphere. What happens between walls/boxes/facets and spheres by the way? They usually have null mass as well. Cheers Bruno -- 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
Re: [Yade-users] [Question #406523]: Cylinders and law2_ScGeom_ViscElPhys_Basic
Question #406523 on Yade changed: https://answers.launchpad.net/yade/+question/406523 Raphaël Maurin confirmed that the question is solved: Thanks Bruno Chareyre, that solved my question. -- 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
Re: [Yade-users] [Question #406523]: Cylinders and law2_ScGeom_ViscElPhys_Basic
Question #406523 on Yade changed: https://answers.launchpad.net/yade/+question/406523 Status: Answered => Solved Raphaël Maurin confirmed that the question is solved: Hi Bruno, Sorry for the late answer, I didn't see your message. Thanks for your answer, I start to understand what is happening and that is getting better. When I use your script, it indeed try to create a cylinder with material viscElMat as asked. However, it prints an error message saying that it didn't succeed to compute the interactions coefficient. This is rather normal as the ip2 contains a function to evaluate the kn, ks, cn and cs from the young, poisson, en and es which has been made for spheres only, and does not work for the cylinder. Therefore, it is possible to make it work by assigning directly the values of kn, ks, cn and cs to define the material. Then, we can observe the viscous dissipation as expected. Here is the script that works: from yade.gridpfacet import * #O.materials.append(ViscElMat(en=0.5, et=0., young=1e7, poisson=0.5, density=2500, frictionAngle=0.4, label='Mat')) O.materials.append(ViscElMat(kn=1e7, ks = 5e6, cn = 1e5,cs = 0., density=2500, frictionAngle=0.4, label='Mat')) O.engines = [ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_GridConnection_Aabb(),]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom(),Ig2_Sphere_GridConnection_ScGridCoGeom()], [Ip2_ViscElMat_ViscElMat_ViscElPhys()], [Law2_ScGeom_ViscElPhys_Basic()], ), NewtonIntegrator(damping=0., gravity = (0,0,-9.81)) ] O.bodies.append(sphere(center = (0,5,5),radius = 1,material = 'Mat')) cylinder(begin = (0,0,0),end =(0,10,0),radius = 1,fixed = True,color = (0,0,1),intMaterial = 'Mat',extMaterial = 'Mat') O.dt = 1e-6 I will use it like that, as it is not a major problem for me. Thanks for everything ! Raphael -- 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
