Re: [Yade-users] [Question #676012]: Questions about facet and CohesiveFrictionalContactLaw
Question #676012 on Yade changed: https://answers.launchpad.net/yade/+question/676012 Gastbye Liu gave more information on the question: Hi everyone, For example, I think I can add shear cohesion between facets and spheres and I use the material CohFricMat. However, I don't know how to add normal cohesion. When facets and spheres are not physically contact (split), the normal direction and contact point is hard to define? How to calculate the tensile force between facets and spheres? Thanks! Best regards, Guangyu -- 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 #676012]: Questions about facet and CohesiveFrictionalContactLaw
Question #676012 on Yade changed: https://answers.launchpad.net/yade/+question/676012 Status: Answered => Open Gastbye Liu is still having a problem: Dear Bruno, Thanks for your reply! I am not sure about the name of the model. You know, the cohesion exists in spherical particles, in yade, the material is CohFrictMat. However, between facet and particles, I only can set Coulomb friction in yade. But I want to add cohesion between facet and spheres. How to do that? Use which kind of material? Best regards, Guangyu -- 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 #676012]: Questions about facet and CohesiveFrictionalContactLaw
New question #676012 on Yade: https://answers.launchpad.net/yade/+question/676012 Hi everyone, I am intertested about the cohesion between particles and I manage to get the algorithm of that. But now, I want to add cohesion between the sphere particle and facet. I am concerned about the contact normal, penetration depth and contact point when facet and sphere are split but there are still tensile force between sphere and facet. In Yade, how does that solve? The code is in this file CohesiveFrictionalContactLaw.cpp? Best wishes, Guangyu -- 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 #667267]: Questions about CohesiveFrictionalContactLaw
Question #667267 on Yade changed: https://answers.launchpad.net/yade/+question/667267 Status: Answered => Solved Gastbye Liu 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 #667267]: Questions about CohesiveFrictionalContactLaw
Question #667267 on Yade changed: https://answers.launchpad.net/yade/+question/667267 Status: Answered => Open Gastbye Liu is still having a problem: Ohhh I got it. Thanks! So it doesn't consider the torque in CohesiveFrictionalContactLaw? -- 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 #667267]: Questions about CohesiveFrictionalContactLaw
Question #667267 on Yade changed: https://answers.launchpad.net/yade/+question/667267 Status: Answered => Open Gastbye Liu is still having a problem: In https://github.com/yade/trunk/blob/master/pkg/dem/CohesiveFrictionalContactLaw.cpp I can see: Vector3r moment = phys->moment_twist + phys->moment_bending; But when calculating moment, the tangential force would affect moment. M = Fs * R. So... I thought in CohesiveFrictionalContactLaw, it doesn't consider the moment Fs * R ? -- 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 #667267]: Questions about CohesiveFrictionalContactLaw
Question #667267 on Yade changed: https://answers.launchpad.net/yade/+question/667267 Status: Answered => Open Gastbye Liu is still having a problem: Thanks for your reply! By the way, in CohesiveFrictionalContactLaw, does it consider the the moment caused by tangential force? -- 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 #667267]: Questions about CohesiveFrictionalContactLaw
Question #667267 on Yade changed: https://answers.launchpad.net/yade/+question/667267 Gastbye Liu posted a new comment: Thanks for your reply! By the way, in CohesiveFrictionalContactLaw, does it consider the the moment caused by tangential force? -- 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 #667267]: Questions about CohesiveFrictionalContactLaw
Question #667267 on Yade changed: https://answers.launchpad.net/yade/+question/667267 Gastbye Liu posted a new comment: And in cohesion model how does the cohesion break? Through reading the code, I think when cohesion break because of two situlations. One: the tensile force between particles exceeds the normal adhesion. Then there would be no force between particles. Second: The shear force exceed the shear adhesion and there would be no force between particles. Am I right? Thanks! -- 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 #667267]: Questions about CohesiveFrictionalContactLaw
New question #667267 on Yade: https://answers.launchpad.net/yade/+question/667267 Hi everyone, I am using the CohesiveFrictionalContactLaw and trying to understand how it works. I noticed: contactPhysics->normalAdhesion = normalAdhPreCalculated*pow(std::min(geom->radius2, geom->radius1),2); contactPhysics->shearAdhesion = shearAdhPreCalculated*pow(std::min(geom->radius2, geom->radius1),2); 1) Does that mean normalAdhesion = tensileStrength * pow(min(R1,R2),2) shearAdhesion= shearStrength * pow(min(R1,R2),2)? But I saw in papers the normalAdhesion should be equal to tensileStrength * pow(min(R1,R2),2) * Pi. I was wondering that whether I neglect Pi in some place... 2)If I set fragile = true, so the value unp is not used. Right? Best Wishes, Liu -- 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 #631723]: Add Force every iterPeriod
Question #631723 on Yade changed: https://answers.launchpad.net/yade/+question/631723 Status: Answered => Solved Gastbye Liu confirmed that the question is solved: Dear Jan: Many Thanks! That's solved my question. Liu -- 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 #631723]: Add Force every iterPeriod
New question #631723 on Yade: https://answers.launchpad.net/yade/+question/631723 Hi everyone: I want to apply force on the polyhedron and the force grows as the time grows. Question1: This is my code. Is that correct? Question2: Also I don't know how to print my added force in every iterPeriod. Could you give me some advice? from yade import * from yade import polyhedra_utils matF = PolyhedraMat() matF.density = 10e15 matF.young = 2e11 matF.poisson = 0.1 matF.Ks = 1e14 matF.Kn = 1e15 matF.frictionAngle = 0.2 yade.O.materials.append(matF) # basic YADE engines fo tetrahedra contact O.bodies.append(polyhedra_utils.polyhedra(matF,v=[[-0.025,-0.025,-0.005],[-0.025,0.025,-0.005],[0.025,-0.025,-0.005],[0.025,0.025,-0.005],[-0.025,-0.025,0],[-0.025,0.025,0],[0.025,-0.025,0],[0.025,0.025,0]])) O.bodies.append(polyhedra_utils.polyhedra(matF,v=[[-0.025,-0.025,0.1],[-0.025,0.025,0.1],[0.025,-0.025,0.1],[0.025,0.025,0.1],[-0.025,-0.025,0.105],[-0.025,0.025,0.105],[0.025,-0.025,0.105],[0.025,0.025,0.105]])) O.engines = [ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Polyhedra_Aabb(),Bo1_Wall_Aabb()],verletDist=0.1), InteractionLoop( [Ig2_Wall_Polyhedra_PolyhedraGeom(), Ig2_Polyhedra_Polyhedra_PolyhedraGeom(), Ig2_Facet_Polyhedra_PolyhedraGeom()], [Ip2_PolyhedraMat_PolyhedraMat_PolyhedraPhys()], # collision "physics" [Law2_PolyhedraGeom_PolyhedraPhys_Volumetric()] # contact law -- apply forces ), #GravityEngi NewtonIntegrator(damping=0.2), PyRunner(command='addForce()',iterPeriod=1) ] print "utils.PWaveTimeStep= "+str(utils.PWaveTimeStep()) def addForce(): O.forces.addF(0,Vector3(0,0,5*O.iterPeriod)) O.forces.addF(1,Vector3(0,0,-5*O.iterPeriod)) ### Many Thanks! Liu -- 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 #574421]: meaning of Ig2_facet_sphere.cpp
Question #574421 on Yade changed: https://answers.launchpad.net/yade/+question/574421 Status: Answered => Solved Gastbye Liu confirmed that the question is solved: Hi Jan: Thank you so much! Yeah, I mean "Inscribing circle radius" instead of "circular radius in facet". This solved my problem absolutely! How kind you are! Guangyu -- 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 #574421]: meaning of Ig2_facet_sphere.cpp
Question #574421 on Yade changed: https://answers.launchpad.net/yade/+question/574421 Gastbye Liu gave more information on the question: And what's the meaning of "circular radius in Facet.hpp"? Why facet has circular radius? -- 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 #574421]: meaning of Ig2_facet_sphere.cpp
New question #574421 on Yade: https://answers.launchpad.net/yade/+question/574421 hi everyone: "Vector3r cl " represents the center of sphere and the the center of facet? The coordinate of facet's normal belongs to global coordinate system? And the " Vector3r cl" belongs to local coordinate system? What's the meaning of "shrinkFactor" ? Do you have some introduction to help me understand that? Thanks! Guangyu -- 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 #562190]: Change about Ig2_Facet_Sphere_ScGeom.cpp
Question #562190 on Yade changed: https://answers.launchpad.net/yade/+question/562190 Gastbye Liu posted a new comment: hi Jan Thanks! I would try! -- 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 #562190]: Change about Ig2_Facet_Sphere_ScGeom.cpp
Question #562190 on Yade changed: https://answers.launchpad.net/yade/+question/562190 Status: Solved => Open Gastbye Liu is still having a problem: Hi: Many Thanks to your reply! You are very lovely! I have compared many code. For example, comparing "Ig2_Facet_Sphere_ScGeom.cpp" (Yade's version is 1.07) and "Ig2_Facet_Sphere_ScGeom.cpp"(Yade's version is 2017.07a) , it's all the same. And then I changed the function content in "Ig2_Facet_Sphere_ScGeom_new.cpp" (Yade's version is 2017.07a) and made the funtion's code are just the same like the code in " Ig2_Facet_Sphere_ScGeom.cpp" . The values don't have "nan" at all. It comes out to me that maybe my new.cpp have bugs. Next I plan to make some output(such as 'Vector3r cl' and so on) to terminal in the new.cpp to check every step. Jan,what do you mean? What's the meaning of define some scenarios? Change different algothrim in new.cpp? -- 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 #562190]: Change about Ig2_Facet_Sphere_ScGeom.cpp
Question #562190 on Yade changed: https://answers.launchpad.net/yade/+question/562190 Gastbye Liu posted a new comment: you are right, I need to change that. -- 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 #562190]: Change about Ig2_Facet_Sphere_ScGeom.cpp
Question #562190 on Yade changed: https://answers.launchpad.net/yade/+question/562190 Status: Open => Solved Gastbye Liu 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 #562190]: Change about Ig2_Facet_Sphere_ScGeom.cpp
Question #562190 on Yade changed: https://answers.launchpad.net/yade/+question/562190 Status: Answered => Open Gastbye Liu is still having a problem: Thanks for your reply! After your advise, I made some output on datafile, which contains O.iterations.id1 and O.iterations.id2. I compared the version output on the same condition. On the first time step, the difference is clear. The two version's output data is same but the order changed. ( The first step don't use my Ig2_Facet_Sphere_ScGeom.cpp. ) When I use my new.cpp do calculation, then the second step has wrong data and the wrong data is more and more as the timestep grows. What should I do? Make change on my new.cpp ? -- 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 #562190]: Change about Ig2_Facet_Sphere_ScGeom.cpp
Question #562190 on Yade changed: https://answers.launchpad.net/yade/+question/562190 Gastbye Liu posted a new comment: Thanks for your reply! After your advise, I made some output on datafile, which contains O.iterations.id1 and O.iterations.id2. I compared the version output on the same condition. On the first time step, the difference is clear. The two version's output data is same but the order changed. ( The first step don't use my Ig2_Facet_Sphere_ScGeom.cpp. ) When I use my new.cpp do calculation, then the second step has wrong data and the wrong data is more and more as the timestep grows. What should I do? Make change on my new.cpp ? -- 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 #562190]: Change about Ig2_Facet_Sphere_ScGeom.cpp
Question #562190 on Yade changed: https://answers.launchpad.net/yade/+question/562190 Status: Answered => Open Gastbye Liu is still having a problem: How cound I use the patch to be applied on trunk? -- 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 #562190]: Change about Ig2_Facet_Sphere_ScGeom.cpp
New question #562190 on Yade: https://answers.launchpad.net/yade/+question/562190 Because I want to do FEM*DEM coupling project, I have to add functions about the geo-information for the contact between facet and sphere. My senior had made a new cpp named Ig2_Facet_Sphere_ScGeom_new.cpp, it worked on the version Yade-1.07. Just put the new.cpp on the folder of pkg/dem/, then compile it and make. But now I use the new version Yade-2017.01a and in the same way compile it and make and make install. It seems that there is no strange thing happened. But when I use that to calculate, I found problem. After initialization and O.engine started working, the vtk-output said some data are "nan" . In other words, the sphere's position are not right. But I don't use the new.cpp, it works fine. I tried but it also don't work. I am wondering that it need another change (in other file) beacause of the version changing. Please help me. Hope your reply. Many Thanks. -- 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