yade-users team mailing list archive

[xjin <question622060@xxxxxxxxxxxxxxxxxxxxx>]

New question #622060 on Yade:
https://answers.launchpad.net/yade/+question/622060

I am trying to run a packing model,and try to use TriaxialTest() to obtain the strain-changing-field. I write the script following the example: tesselationWrapper.py .

My specific script is like tihs:
-----------------------------------------------------------------------------------------------------------------
O.reset()
O.bodies.append(geom.facetBox((.5,.5,.5),(.5,.5,.5),wallMask=63))
from yade import pack
sp=pack.SpherePack()
sp.makeCloud((0,0,0),(1,1,1),rMean=0.02, rRelFuzz=0.1)
sp.toSimulation()
O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
   InteractionLoop(
      [Ig2_Sphere_Sphere_L3Geom(),Ig2_Facet_Sphere_L3Geom()],
      [Ip2_FrictMat_FrictMat_FrictPhys()],
      [Law2_L3Geom_FrictPhys_ElPerfPl()]
   ),
   NewtonIntegrator(gravity=(0,0,-9.81),damping=0.4),
   
]
O.dt=0.05*PWaveTimeStep()
O.run(5)

TW=TesselationWrapper()
TW.triangulate()
TW.computeVolumes()
TW.volume(20)
TW.setState(0)
O.run(10,True)
TW.setState(1)
TW.defToVtk("strain.vtk")
-----------------------------------------------------------------------------------------------------------------------------------------
when I run the script, the result is:
-----------------------------------------------------------------------------------------------------------------------------------------
Triangulated Grains : 7432
sym_grad_u_total_g (wrong averaged strain):
0 0 -1.0106e-17 
0 0 -1.0839e-17 
-1.0106e-17 -1.0839e-17 -7.99784e-17 

Total volume = 0.87657, grad_u = 
0 0 0 
0 0 0 
-2.0212e-17 -2.16781e-17 -7.99784e-17 

sym_grad_u (true average strain): 
0 0 -1.0106e-17 
0 0 -1.0839e-17 
-1.0106e-17 -1.0839e-17 -7.99784e-17 

Macro strain : 
0 0 0 
0 0 0 
0 0 0 

Cannot open file [strain.vtk]
-----------------------------------------------------------------------------------------------------------------------------------------
Can someone help me? For I am new to YADE,I cannot understand the meaning of the result.
Thanks very much!!!!

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.