Failed to triangulate body
Hi,
I use the FlowEngine and get some problems. I use the pack module to generate particles, when I set the seed=2 or 3, I get such warning:
Failed to triangulate body with id=47 Point=-1976.72 1535.93 -26.6935 rad=0.03
Failed to triangulate body with id=57 Point=1343.75 -3820.83 -831.97 rad=0.03
.........
NEGATIVE DETERMINANT
NEGATIVE DETERMINANT
Thank you for your attention to this matter !
May code is shown as follow:
#======
from yade import pack,qt,plot
num_spheres=1000# number of spheres
young=1e6
compFricDegree = 3 # initial contact friction during the confining phase
finalFricDegree = 30 # contact friction during the deviatoric loading
mn,mx=Vector3(
flag=1
flag2=0
O.materials.
O.materials.
walls=aabbWalls
wallIds=
sp=pack.
sp.makeCloud(
sp.makeCloud(
sp.toSimulation
triax=TriaxialS
thickness = 0,
stressMask = 7,
max_vel = 0.005,
internalCompac
# wall_bottom_
)
newton=
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
),
FlowEngine(
GlobalStiffnes
triax,
newton,
]
triax.goal1=
while 1:
O.run(1000,1)
unb=unbalance
if unb<0.001 and abs(-10000-
break
setContactFrict
## ______________ Oedometer section _________________
triax.stressMask=7
triax.goal1=
#B. Activate flow engine and set boundary conditions in order to get permeability
flow.dead=0
flow.defToleran
flow.meshUpdate
flow.useSolver=3
flow.permeabili
flow.viscosity=10
flow.bndCondIsP
flow.bndCondVal
flow.boundaryUs
flow.updateTria
#O.dt=0.1e-4
#O.dynDt=False
newton.damping=0
#triax.
O.run(2,1)
csdList=
print(len(
flow.dead=1
for i in range(10,900):
O.bodies.
O.run(1000,1)
O.engines=
O.step()
O.step()
O.engines=
flow2.defTolerance = -1
flow2.meshUpdat
flow2.useSolver = 3
flow2.permeabil
flow2.viscosity = 10
flow2.bndCondIs
flow2.bndCondValue = [0, 0, 0, 0, 100, 0]
flow2.boundaryU
flow2.updateTri
O.run(2,1)
csdList=
print(len(
Question information
- Language:
- English Edit question
- Status:
- Expired
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply: