Understanding why sometimes does not occurs contact between Facet/Spheres using the same script
Dear all,
I am starting a script to reach a specified porosity using REFD method. But right at the beginning, I have found some code behaviors that I still do not understand. In the script below:
1) If I use num=1 in the makeCloud there are no contacts between Facet/Sphere;
2) In the same script, if I activate the gravity=(0,0,-9,81) in NewtonIntegrator, there are contacts between Facet/Sphere;
3) Even without using gravity, but if I use, for example, 100 spheres (num=100 in the makeCloud), there are contacts between Facet/Sphere.
I'm a little unsure of these behaviors, could you tell me what's going on?
Thank you very much.
Marcus.
#############
## Script used ##
#############
from yade import utils,pack,
radius=.05
damp=0.2
O.materials.
O.materials.
O.bodies.
sp=pack.
sp.makeCloud(
sp.toSimulation
triax=TriaxialS
maxMultiplier=
finalMaxMultip
stressMask = 7,
internalCompac
)
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
),
NewtonIntegrat
PyRunner(
PyRunner(
]
O.dt=.5*
O.trackEnergy=True
def radiusExpansion():
for n in range(0,
if isinstance(
if O.bodies[
O.bodies[
def contactForces():
for i in O.interactions:
if not i.isReal: continue
cp = i.geom.contactPoint
nf = i.phys.normalForce
sf = i.phys.shearForce
print cp
def fabricTensor():
utils.
rr=yade.
plot.plots=
' i ':(O.energy.
}
yade.qt.View()
O.saveTmp()
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: