yade-users team mailing list archive

[Klaus Thoeni <question177203@xxxxxxxxxxxxxxxxxxxxx>]

Question #177203 on Yade changed:
https://answers.launchpad.net/yade/+question/177203

    Status: Open => Answered

Klaus Thoeni proposed the following answer:
Hi Ning,

you have to check out and compile the trunk version in order to get new 
commits. Yade version 0.70.0 will not be updated!

However, I tried it and there is still a problem with the 2d version.

@Bruno, did you test your implementation? Yade is calculating and calculating 
but I can't get any packing.

I guess 2d remains still experimental ;-)

Klaus

On Wed, 2 Nov 2011 02:35:44 PM ceguo wrote:
> New question #177203 on Yade:
> https://answers.launchpad.net/yade/+question/177203
> 
> Hi,
> 
> In the discussion of #175606, Bruno committed a change of sp.makeCloud to
> accommodate true 2d tests. Where is the new version code? I recompiled
> code using version 0.70.0, but still cannot generate packing with 0
> thickness.  The error info:
> 
> num =
> sp.makeCloud(minCorner=(0,0,0),maxCorner=(0.01,0.01,0),rMean=0.0002,rRelFu
> zz=.5,periodic=True) ValueError: The box defined as null volume. Define at
> least maxCorner of the box, or hSize if periodic.
> 
> Also with this zero-thickness packing, how to conduct biaxial test? Still
> using Peri3dController is OK? Do we need to block the corresponding DOFs?
> Will the following code work smoothly?
> 
> from yade import pack
> 
> O.materials.append()
> sp = pack.SpherePack()
> num =
> sp.makeCloud(minCorner=(0,0,0),maxCorner=(1,1,0),rMean=0.02,rRelFuzz=.5,pe
> riodic=True) O.bodies.append([utils.sphere(s[0],s[1]) for s in sp])
> print num
> for p in O.bodies:
>    p.state.blockedDOFs = ['z','rx','ry']
> 
> O.dt = .5*utils.PWaveTimeStep()
> O.engines = [
>    ForceResetter(),
>    InsertionSortCollider([Bo1_Sphere_Aabb()]), # detect new collisions
>    InteractionLoop(
>       [Ig2_Sphere_Sphere_Dem3DofGeom()],
>       [Ip2_FrictMat_FrictMat_FrictPhys()],
>       [Law2_Dem3DofGeom_FrictPhys_CundallStrack()]
>    ),
>    Peri3dController(
>       goal=(-2e5,-2e5,0,0,0,0),
>       stressMask=3,
>       nSteps=20000,
>       doneHook='print "Consolidation finished."; O.pause()',
>       maxStrain=.5,
>       label='triax'
>    ),
>    NewtonIntegrator()
> ]
> O.run(); O.wait()
> 
> O.engines = [
>    ForceResetter(),
>    InsertionSortCollider([Bo1_Sphere_Aabb()]), # detect new collisions
>    InteractionLoop(
>       [Ig2_Sphere_Sphere_Dem3DofGeom()],
>       [Ip2_FrictMat_FrictMat_FrictPhys()],
>       [Law2_Dem3DofGeom_FrictPhys_CundallStrack()]
>    ),
>    Peri3dController(
>       goal=(2.e-2,-2.e-2,0,0,0,0),
>       stressMask=0,
>       nSteps=20000,
>       doneHook='print "Shear finished."; O.pause()',
>       maxStrain=.5,
>       label='triax'
>    ),
>    NewtonIntegrator(),
> ]
> O.run(); O.wait()
> 
> Ning

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.