Question #656214 on Yade changed:
https://answers.launchpad.net/yade/+question/656214
Status: Open => Answered
Robert Caulk proposed the following answer:
Hello Jesse,
>Since the clumps are still not being compressed, I believe the MWE I
posted before should still be valid.
Just because you believe it, does not make it true ;-) [1].
Here is a script that uses makeclumpcloud together with
periTriaxController.
from yade import pack,qt
O.periodic=True
O.cell.hSize=Matrix3(1, 0, 0,
0 ,1, 0,
0, 0, 1)
sp=pack.SpherePack()
sp=pack.SpherePack()
c1=pack.SpherePack([((0,0,0),.1),((.15,0,0),.05),((0,.1,0),.05)])
sp.makeClumpCloud(Vector3().Zero,O.cell.refSize,[c1],periodic=True)
sp.toSimulation()
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()]
),
#PeriTriaxController(maxUnbalanced=0.01,relStressTol=0.02,goal=[-1e4,-1e4,0],stressMask=3,globUpdate=5,maxStrainRate=[1.,1.,1.],doneHook='triaxDone()',label='triax'),
#using cell inertia
PeriTriaxController(dynCell=True,mass=0.2,maxUnbalanced=0.01,relStressTol=0.02,goal=(-1e4,-1e4,0),stressMask=3,globUpdate=5,maxStrainRate=(1.,1.,1.),doneHook='triaxDone()',label='triax'),
NewtonIntegrator(damping=.2),
]
O.dt=PWaveTimeStep()
O.run();
qt.View()
phase=0
def triaxDone():
global phase
if phase==0:
print 'Here we are:
stress',triax.stress,'strain',triax.strain,'stiffness',triax.stiff
print 'Now εz will go from 0 to .2 while σx and σy will be kept
the same.'
triax.goal=(-1e4,-1e4,-0.2)
phase+=1
elif phase==1:
print 'Here we are:
stress',triax.stress,'strain',triax.strain,'stiffness',triax.stiff
print 'Done, pausing now.'
O.pause()
[1]https://yade-dem.org/wiki/Howtoask
--
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 : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
