yade-users team mailing list archive

[Gavin He <question680307@xxxxxxxxxxxxxxxxxxxxx>]

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

Recenetly I have tried to run periodic triaxial test on polyhedras. Now I can use polyhedra_utils.fillbox function to create a box of polyhedras, but I cannot run PeriTriaxController on them. Could anyone help me to solve this problem?

____Code______________________________________

sigmaIso=-1e5

from yade import pack, qt, plot, polyhedra_utils

gravel = PolyhedraMat()
gravel.density = 2600 #kg/m^3 
gravel.young = 1E7 #Pa
gravel.poisson = 20000/1E7
gravel.frictionAngle = 0.5 #rad

O.periodic=True

polyhedra_utils.fillBox((0,0,0), (2,2,2),gravel,sizemin=[0.2,0.2,0.2],sizemax=[0.2,0.2,0.2],seed=4)


O.engines=[
      ForceResetter(),
      InsertionSortCollider([Bo1_Polyhedra_Aabb(),Bo1_Wall_Aabb(),Bo1_Facet_Aabb()]),
      InteractionLoop(
      [Ig2_Wall_Polyhedra_PolyhedraGeom(), Ig2_Polyhedra_Polyhedra_PolyhedraGeom(), Ig2_Facet_Polyhedra_PolyhedraGeom()], 
      [Ip2_PolyhedraMat_PolyhedraMat_PolyhedraPhys()], # collision "physics"
      [Law2_PolyhedraGeom_PolyhedraPhys_Volumetric()]   # contact law -- apply forces
   ),
	PeriTriaxController(label='triax',
		# specify target values and whether they are strains or stresses
		goal=(sigmaIso,sigmaIso,sigmaIso),stressMask=7,
		# type of servo-control
		dynCell=True,maxStrainRate=(10,10,10),
		# wait until the unbalanced force goes below this value
		maxUnbalanced=.1,relStressTol=1e-3,
		# call this function when goal is reached and the packing is stable
		doneHook='compactionFinished()'
	),
	NewtonIntegrator(damping=.2),
	PyRunner(command='addPlotData()',iterPeriod=100),
]
O.dt=.5*PWaveTimeStep()

_______________________________________________________________

Before simulation starts:

generated 253 polyhedrons
WARN  /build/yade-KKgSmd/yade-1.20.0/pkg/dem/Shop_01.cpp:457 PWaveTimeStep: PWaveTimeStep has not found any suitable spherical body to calculate dt. dt is set to 1.0
yade.plot: creating fake plot, since there are no y-data yet
/usr/lib/python2.7/dist-packages/matplotlib/axes/_axes.py:519: UserWarning: No labelled objects found. Use label='...' kwarg on individual plots.
  warnings.warn("No labelled objects found. "

After simulation starts:

ERROR /build/yade-KKgSmd/yade-1.20.0/pkg/common/InsertionSortCollider.cpp:244 action: verletDist is set to 0 because no spheres were found. It will result in suboptimal performances, consider setting a positive verletDist in your script.
yade.plot: creating new line for elastPotential
yade.plot: creating new line for gravWork
yade.plot: creating new line for kinetic
yade.plot: creating new line for plastDissip
yade.plot: creating new line for velGradWork
yade.plot: creating new line for nonviscDamp



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