Question #622060 on Yade changed:
https://answers.launchpad.net/yade/+question/622060
Description changed to:
I am trying to run a packing model,and try to use TesselationWrapper()
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.
_______________________________________________
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
