← Back to team overview

yade-users team mailing list archive

Re: [Question #700314]: Hertz Mindlin Contact Law including Rolling

 

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

Description changed to:
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 ?

Also

2) I want to active this contact law after the sample achieve desire
void ratio, since the computation take a long time to achieve desire
void ratio, So How can I do that?

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.