Calculate the radius of the overlaped area between two adjacent particles.
Hi!
I have a quetion with getting the contact radius of the contact particles.
A Hertz-Minglin contact model is used to describe sphere-sphere interactions when I prepared the sample using "TriaxialStress
I want to calculate the radius of the overlaped area between two adjacent particles. But I do not know how to do.
Are there some methods to achieve it? Thanks!
Here is my code to prepare the sample.
##############
from __future__ import print_function
from yade import pack,plot,os,timing
import matplotlib;
import pylab
#Material constants
density=2650
poisson=0.31
young=29e9
damp=0.25
number=5000
frictionangle=0
#Wall constants
walldensity=0
wallfrictionangle=0
wallpoisson=0.5
wallyoung=30000e9
mn=Vector3(
mx=Vector3(
compress=-5000
rate=0.1
O.materials.
O.materials.
wallIds=
# PSD of the particles
psdSizes,
pylab.plot(
#Add particles
sp = pack.SpherePack()
sp.makeCloud(
pylab.plot(
pylab.legend()
sp.toSimulation
#Color
for b in O.bodies:
if isinstance(
if r>0.00025 and r<0.000375:
if r>0.0 and r<0.00025:
#Energy
O.trackEnergy=True
#Defining triaxil engines
triax=TriaxialS
wall_bottom_id = wallIds[2],
wall_top_id = wallIds[3],
wall_left_id = wallIds[0],
wall_right_id = wallIds[1],
wall_back_id = wallIds[4],
wall_front_id = wallIds[5],
thickness=
internalComp
stressMask=7,
goal1=compress,
goal2=compress,
goal3=compress,
)
newton=
###########
O.engines=[
ForceResetter(),
InsertionSort
InteractionLoop(
),
GlobalStiffness
triax,
newton,
TriaxialStat
PyRunner(
PyRunner(
]
#############
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Jan Stránský
- Solved:
- Last query:
- Last reply: