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.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
