yade-users team mailing list archive

[张文卿 <question562171@xxxxxxxxxxxxxxxxxxxxx>]

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

I’m using PeriTriaxController to get packing with desired porosity.But when I used ‘triax.stress' to get average stresses ,it showed (0,0,0).And it's weird this script can work.So where should I change?Here is my script:
r=1.75e-05
O.periodic=True
O.cell.hSize=Matrix3(41*r,0,0,
		    0,41*r,0,
		    0,0,41*r)

O.materials.append(FrictMat(young=1.1e11,poisson=0.3,density=8900e12,frictionAngle=0,label='sphere'))
sp=pack.SpherePack()
sp.makeCloud((0,0,0),(40*r,40*r,40*r),r,0,4000,False)
O.bodies.append([sphere(center,rad,material='sphere') for center,rad in sp])


############################triax########################################################
triax=PeriTriaxController(goal=(-2e4,-2e4,-2e4),stressMask=7)
############################Engine#######################################################
O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb()],verletDist=0.05*r),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom()],
		[Ip2_FrictMat_FrictMat_MindlinPhys()],
		[Law2_ScGeom_MindlinPhys_HertzWithLinearShear()]
	),
	PyRunner(command='controlporosity()',iterPeriod=1),
	triax,
	NewtonIntegrator(damping=.2)
]

O.dt=PWaveTimeStep()
def controlporosity():
	targateporosity=0.36
	a=utils.porosity()
	if a<targateporosity:
		O.pause()



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