yade-users team mailing list archive

[Ehsan Izadi <question240920@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Answered => Open

Ehsan Izadi is still having a problem:
Thanks, Here is my code, I tried to make it shorter: 
----------------------------------------------------------------------------------------------
box=ymport.stl('/home/ehsan/Desktop/try3.stl')
O.bodies.append(box)

sp=pack.SpherePack()
sp.makeCloud((-0.2,-0.15,0),(0.2,0.15,0.1),rRelFuzz=0.5, num=20000)
sp.toSimulation(material='particles')

triax=TriaxialStressController(
	maxMultiplier=1.+2e4/1e6,
	finalMaxMultiplier=1.+2e3/1e6, 
	thickness = 0,
	stressMask = 7,
	max_vel = 0.005,
	internalCompaction=True 
)


O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
		[Ip2_FrictMat_FrictMat_FrictPhys()],
		[Law2_ScGeom_FrictPhys_CundallStrack()]
	),
	NewtonIntegrator(gravity=(0,0,-9.81),damping=0.2),
	triax,
	GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.8)
	]


triax.goal1=triax.goal3=triax.goal2=1000


O.saveTmp()
----------------------------------------------------------------------------------

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