Facet in the model disappears
Hello everyone!
I am trying to compact a mixed specimen of ball and polyhedra. In this process, I generate a facet as a compaction plate above the specimen. But when the platen was in contact with the specimen, it disappeared immediately.
Here is the script:
from yade import pack,qt,
from yade import *
import numpy
from pprint import pprint
import random
import math
from random import uniform
from random import randint
from math import *
matP = PolyhedraMat()
matP.density = 2650 #kg/m^3
matP.young = 5.5E9#Pa 5.5E9
#matP.young = 2.75 #Pa 1E7
matP.poisson = 0.25 # 20000/1E7
matP.frictionAngle = 0.7 #rad0.7
O.materials.
steel= PolyhedraMat()
steel.density = 2600 #kg/m^3
steel.young = 15E9 #inital steel was 15E9
steel.poisson = 0.21
steel.frictionAngle = 0.2 #rad
O.materials.
sphereballast = CohFrictMat()#
sphereballast.
#sphereballast.
sphereballast.young = 3E8
sphereballast.
sphereballast.
sphereballast.
sphereballast.
sphereballast.
O.materials.
h=0
b=0.15 #if the shape or the size of the boundary need to change, dont forget the contact number in wall's velocity.
mpi=math.pi
O.bodies.
O.bodies.
O.bodies.
O.bodies.
O.bodies.
O.bodies.
O.bodies.
O.bodies.
O.bodies.
O.bodies.
O.bodies.
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
),
NewtonIntegrat
PyRunner(
PyRunner(
]
global a
a=1
def Compact():
global a
if a == 1:
ldpltheight=
fa1=utils.
O.bodies.
fa2=utils.
O.bodies.
a=a+1
elif a == 2:
Lastnum1=
O.bodies[
O.bodies[
O.bodies[
Lastnum=
O.bodies[
O.bodies[
O.bodies[
O.dt=0.000001
The vanishing facets I described above are fa1 and fa2. But as you can see, the boundaries in the model are also facets, and these facets work fine.
How did the above problem arise?
xue
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Karol Brzezinski
- Solved:
- Last query:
- Last reply: