yade-users team mailing list archive

[Hossein <question700314@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi every body

I need to rolling resistance in order to increase the stiffness of the sample during a triaxial test. I have already probe with below parameters such as young modules and Poisson ratio and the result is the same with a simulation including rolling stiffness (krot) and coefficent of bending moment (eta) was performed. So my question is that :

1) How to should I assign rolling and twisting in a simulation with considering Hertz Mindlin contact law ?   

Here in below I mention my code where the contact law must be activated:

O.materials.append(FrictMat(young=9e10,poisson=.25,frictionAngle=radians(compFricDegree),density=2600,label='spheres'))
O.materials.append(FrictMat(young=9e10,poisson=.25,frictionAngle=0,density=0,label='frictionless'))


## create walls around the packing
walls=aabbWalls((mn,mx),thickness=0,material='frictionless')
wallIds=O.bodies.append(walls)


O.bodies.append([utils.sphere(center,rad,material='spheres') for center,rad in sp])


triax=TriaxialStressController(
	internalCompaction=False,
	goal1=confiningS,
	goal2=confiningS,
	goal3=confiningS,
	max_vel=10,
	label="triax"
)

newton=NewtonIntegrator(damping=0.7)

O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb()]),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Box_Sphere_ScGeom()],
		[Ip2_FrictMat_FrictMat_MindlinPhys(krot=.4,eta=.2)],
		[Law2_ScGeom_MindlinPhys_Mindlin(includeAdhesion=False,includeMoment=True)]
	),
	GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.8),
	VTKRecorder(iterPeriod=1000,fileName="./export",recorders=['spheres','intr','color'],dead=1,label='VTK'),
	triax,
	newton
]

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