Small density sphere in micron size makes segment fault (core dumped)
In my case, the sphere used is in micron size. When i use the real material property, especially the density, the calculation crashes after some time(Segment fault (core dumped)). Then i try to magnify the density and find the run goes well when the density amplified by one million times.
Following is the testing code where the parameter masscof is the density amplification factor,
from yade import geom,pack
import math
from yade import qt
import os
os.system('rm -rf plt;mkdir ./plt/')
os.system('rm -rf out;mkdir ./out/')
# PhysicalParameters
muS = 0.57735 # Friction coefficient (30 degree)
FricAngleS = math.atan(muS)
### mass amplification factor ###
masscof = 1e6
# PhysicalParameters
matSph = CohFrictMat(
density = 7990*masscof,
young = 193e9,
poisson = 0.3,
frictionAngle = FricAngleS,
momentRotationLaw = True)
SMat = O.materials.
# create empty sphere packing
sp = pack.SpherePack()
# generate spheres1 with a certain diameter distribution
sp.makeCloud(
# add the sphere pack to the simulation
sp.toSimulation
# create rectangular box1(Left) from facets
O.bodies.
O.bodies.
O.dt = 0.85*utils.
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
# handle sphere+sphere and facet+sphere collisions
),
NewtonIntegrat
qt.
VTKRecorder
]
#show geometry:
qtr = qt.Renderer()
qtr.bgColor = [1,1,1]
qt.Controller()
v = qt.View()
#v.center()
v.axes=False
v.viewDir=
v.eyePosition=
#O.run()
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: