Re: [Yade-users] [Question #698444]: REFD method with clumps
Question #698444 on Yade changed: https://answers.launchpad.net/yade/+question/698444 Status: Answered => Solved panjinhong confirmed that the question is solved: Thanks Jérôme Duriez, 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 #679515]: How to have a constant-p simulations with periodic boundary condition?
Question #679515 on Yade changed: https://answers.launchpad.net/yade/+question/679515 Status: Answered => Solved panjinhong confirmed that the question is solved: Thanks 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
[Yade-users] [Question #698444]: REFD method with clumps
New question #698444 on Yade: https://answers.launchpad.net/yade/+question/698444 Hello everybody! I am trying to use the TriaxialStressController to prepare clumps samples using radius-expansion friction-decrease method. However, I found that the ratio between the distance of the spheres in clumps and the radius of spheres is changing during the process of radius-expansion. I notice that Marcus has experienced the same problem as me [1], and this question was solved with the help of Bruno. But the solution in [2] is not accessible for me now. Is there any solution to this problem or example using REFD method with clumps? Thank you very much. Jinhong [1] https://answers.launchpad.net/yade/+question/661926 [2] https://github.com/yade/trunk/commit/243ee130ba97e628763d36c1b2442700e38d9185 -- 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 #682985]: A strange phenomenon in constant-p(mean stress) simulations
Question #682985 on Yade changed: https://answers.launchpad.net/yade/+question/682985 Description changed to: Hello, everyone! I'm trying to take a constant-p simulation under different b(intermediate principal stress ratio) with Yade but the result seems to be a bit strange. I used non-spherical clumped particle consisting of two identical spheres and periodic boundary condition. In the Stress ratio-Deviation strain curve with different b(sorry I don't know how to upload the picture), When b=0.0 and 1.0, the curve is quite normal whose stress ratio has a peak and then drops, becomes steady. But when b=0.25 and 0.25, after the peak, the stress ratio drops and then rises. The residual value of it is even higher than the stress ratio when deviation strain is 0.3. In my mind, this phenomenon would not happen in reality. Actually, some researchers also do constant-p simulation with PFC-3D or Liggghts, but their results didn't show this phenomenon. Looking forward to your help, Thanks! This is parts of my code: O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb()]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom()], [Ip2_FrictMat_FrictMat_FrictPhys()], [Law2_ScGeom_FrictPhys_CundallStrack()] ), PeriTriaxController(label='triax', # specify target values and whether they are strains or stresses goal=(sigmaIso1,sigmaIso1,sigmaIso1),stressMask=7, # type of servo-control dynCell=True,maxStrainRate=(10,10,10), # wait until the unbalanced force goes below this value maxUnbalanced=0.01,relStressTol=1e-3, # call this function when goal is reached and the packing is stable doneHook='compactionFinished1()' ), NewtonIntegrator(damping=.2), PyRunner(command='check()',iterPeriod=100), ] O.dt=.5*PWaveTimeStep() def compactionFinished(): O.cell.trsf=Matrix3.Identity triax.goal=(sigmaIso,sigmaIso,-.8) triax.stressMask=3 triax.maxStrainRate=(1.,1.,.01) triax.doneHook='triaxFinished()' O.engines = O.engines + [PyRunner(command='history()',iterPeriod=200,label='recorder')] O.engines = O.engines + [PyRunner(command='changeGoalStress()',iterPeriod=50,label='checker')] triax.maxUnbalanced=10 def changeGoalStress(): s11=-triax.stress[0] s22=-triax.stress[1] s33=-triax.stress[2] r=-((b + 1.0) * s33 -1 * 3 * 10) / (b - 2.0) s=-((1.0 - 2 * b ) * s33 + (b - 1) * 3 * 10) / (b - 2.0) triax.goal=(r,s,-0.8) -- 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 #682985]: A strange phenomenon in constant-p(mean stress) simulations
New question #682985 on Yade: https://answers.launchpad.net/yade/+question/682985 Hello, everyone! I'm trying to take a constant-p simulation under different b(intermediate principal stress ratio) with Yade but the result seems to be a bit strange. I used non-spherical clumped particle consisting of two identical spheres and periodic boundary condition. And this is the Stress ratio-Deviation strain curve with different b. As you can see, when b=0.0 and 1.0, the curve is quite normal whose stress ratio has a peak and then drops, becomes steady. But when b=0.25 and 0.25, after the peak, the stress ratio drops and then rises. The residual value of it is even higher than the stress ratio when deviation strain is 0.3. In my mind, this phenomenon would not happen in reality. Actually, some researchers also do constant-p simulation with PFC-3D or Liggghts, but their results didn't show this phenomenon. Looking forward to your help, Thanks! This is parts of my code: O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb()]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom()], [Ip2_FrictMat_FrictMat_FrictPhys()], [Law2_ScGeom_FrictPhys_CundallStrack()] ), PeriTriaxController(label='triax', # specify target values and whether they are strains or stresses goal=(sigmaIso1,sigmaIso1,sigmaIso1),stressMask=7, # type of servo-control dynCell=True,maxStrainRate=(10,10,10), # wait until the unbalanced force goes below this value maxUnbalanced=0.01,relStressTol=1e-3, # call this function when goal is reached and the packing is stable doneHook='compactionFinished1()' ), NewtonIntegrator(damping=.2), PyRunner(command='check()',iterPeriod=100), ] O.dt=.5*PWaveTimeStep() def compactionFinished(): O.cell.trsf=Matrix3.Identity triax.goal=(sigmaIso,sigmaIso,-.8) triax.stressMask=3 triax.maxStrainRate=(1.,1.,.01) triax.doneHook='triaxFinished()' O.engines = O.engines + [PyRunner(command='history()',iterPeriod=200,label='recorder')] O.engines = O.engines + [PyRunner(command='changeGoalStress()',iterPeriod=50,label='checker')] triax.maxUnbalanced=10 def changeGoalStress(): s11=-triax.stress[0] s22=-triax.stress[1] s33=-triax.stress[2] r=-((b + 1.0) * s33 -1 * 3 * 10) / (b - 2.0) s=-((1.0 - 2 * b ) * s33 + (b - 1) * 3 * 10) / (b - 2.0) triax.goal=(r,s,-0.8) -- 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 #679516]: How to get the Volume in the simulation with periodic boundary condition?
Question #679516 on Yade changed: https://answers.launchpad.net/yade/+question/679516 Status: Answered => Solved panjinhong confirmed that the question is solved: Thank you! -- 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 #679516]: How to get the Volume in the simulation with periodic boundary condition?
Question #679516 on Yade changed: https://answers.launchpad.net/yade/+question/679516 Status: Answered => Open panjinhong is still having a problem: Thank you, Jan. By the second question, I mean that when I output the position of all the particle (by state.pos) and observed it in paraview, it become strange. But I want it to like a cuboid, because I want to calculate the stress inside and the unevenness of the boundary cause lots of troubles. Jinhong -- 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 #679516]: How to get the Volume in the simulation with periodic boundary condition?
New question #679516 on Yade: https://answers.launchpad.net/yade/+question/679516 Hello, everyone. I use PeriTriaxController to have a simulation with periodic boundary condition. But I don't know how to measure the volume change. It seems a silly problem, but I need your help! Actually, when I use a Cube with 30,000 particles to do the simulation, I can not keep the boundary's shape (I want to keep it in as a flat at each surface). Is there any way to achieve my demand? Looking forward to your help. 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 #679515]: How to have a constant-p simulations with periodic boundary condition?
New question #679515 on Yade: https://answers.launchpad.net/yade/+question/679515 Hello, everyone. I am looking forward to a way to do a constant-p simulation with periodic boundary condition. I have done some simulations with TriaxialStressController. At that case, I can change the stress after a few steps. Because the stress it provide by the wall would change immediately, I can get a constant-p simulation. But, when I use periodic boundary condition (PeriTriaxController), I could only change the Goal Stress. It means that the stress would not change to the target stress I want immediately, which result in a shaking of the stress (I change the stress quiet often to keep a constant p ). Is there some good ways to achieve this kind of simulations? Looking forward to your help, 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