Why use ymport.textPolyhedra command make the simulation to slow?
I first prepare some polyhedra particle in my model and used 'export.
#######
# INPUTS
#######
gravity = 100
# box dimensions
widthl = 1.3
widthr = 1.3
widthc = 1.3
height = 1
thick = 3
deep = -.2
# size of grains
sizeMin = 40e-3
sizeMax = 95e-3
frictionAngle = .5
young = 45e9 # stiffness
density = 16000 #kg/m^3
dt = 1e-4 # time step
# how much larger the initial make box should be
fillBoxHFactor = 1
#######
from yade import polyhedra_utils, export, ymport
width = widthl+
# mat, engines, ...
mat = PolyhedraMat(
O.materials.
O.engines=[
ForceResetter(),
InsertionSor
InteractionLoop(
[
[
[
),
NewtonIntegr
]
O.dt = dt
# box
p000 = Vector3(0,0,0)
p100 = Vector3(widthl,0,0)
p200 = Vector3(
p300 = Vector3(
pxs = (p000,p100,
p001,p101,p201,p301 = [p+Vector3(
p010,p110,p210,p310 = [p+Vector3(
p011,p111,p211,p311 = [p+Vector3(
p00b,p10b,p20b,p30b = [p+Vector3(
p01b,p11b,p21b,p31b = [p+Vector3(
def rect(vs,**kw):
v1,v2,v3,v4 = vs
return [
facet(
facet(
]
rects = (
(p000,
(p000,
(p300,
(p000,
(p010,
)
rects = sum((rect(r) for r in rects),[])
O.bodies.
O.bodies.
Question information
- Language:
- English Edit question
- Status:
- Answered
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Kun Zeng for more information if necessary.