Kn and Ks value for Hertz-Mindilin model
This was an issue I ran into before with the Cundall-Strack law, and I am now also seeing with Hertz-Mindlin law. Hence, I am trying to input specific value for kn and ks for this contact model and the formulation proposed by yade, but I do not know what I can consider them in my simulation.
Is there a way to use these values? Would a simple Ip2_FrictMat_
You can find herein the necessary part of my script.
Thanks in advance.
Sam
#######
young =3e5 # elastic modulus [Pa]
poisson = 0.25
finalFricDegree=30
density = 1950
mn,mx=Vector3(
radius= 1 # Internal Radius
identifier=
output=
rate = strainRate*10000
stabilityThresh
iterper=1000 # Periodicity criterion using step number (deactivated if <= 0)
kn=5e1 # Normal stiffness
ks=5e1 # Shear Stiffness
kr=0 # Rotational stiffness
Bn=0
Bs=0
intRadius=1
#######
Friction=
#frictionless walls
wallmat = O.materials.
walls=aabbWalls
wallIds=
sp = O.bodies.
triax = TriaxialStressC
thickness=0,
stressMask=7,
internalCompac
)
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
),
GlobalStiffnes
triax,
VTKRecorder(
NewtonIntegrat
]
O.dt = 0.3 * utils.PWaveTime
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Sam
- Solved:
- Last query:
- Last reply: