Sphere goes through the facet wall in Harmonic vibration
Hi, in my case there are lots of spheres in one box and a piston in the bottom of the box. Then the piston does harmonic vibration. I find when the vibration amplitude goes to high, the spheres go through the piston, which is constructed by facet. I want to know the reason and fix it. In addition, the density is amplified by a factor of 1e9 to increase the calculation speed. And also I want to know how to define the word of "rigid" in yade. In my case, the box is rigid and the spheres are relatively soft. Does it make sense if i set the young's modulus of the box material as 100 times of that of the sphere material?
Following is my code,
# PhysicalParameters
muS = 0.57735
muF = 0.17632
FricAngleS = math.atan(muS)
FricAngleF = math.atan(muF)
densitys_ac = 8150*1e9
densityw_ac = 2700*1e9
yongmodu = 195e9
poisn = 0.3
# material model
matSph = CohFrictMat(
SMat = O.materials.
matFacet = CohFrictMat(
FMat = O.materials.
## box and piston
O.bodies.
Cylinder1IDs=
# spreading process motions first harmonic vibratiion then stop
def change_motion():
if O.time > 0.4:
elif O.time > 0.5:
#define engines:
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
# handle sphere+sphere and facet+sphere collisions
),
PyRunner(
HarmonicMot
NewtonInteg
# VTKRecorder(
]
## generation particles
sp = pack.SpherePack()
sp.makeCloud(
sp.toSimulation
## time step
O.dt = 0.85*utils.
O.usesTimeStepper = False
O.run()
Thanks,
Xuesong
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- gaoxuesong
- Solved:
- Last query:
- Last reply: