problem in creating the contact law between two polyhedron blocks
I want to test slide state between two rigid blocks. For this I tried to run the following simple code with Yade 2016.06a-
from yade import plot, polyhedra_utils, utils
from yade import qt
from yade import pack, export
#######
#######
# Materials
matPol=CohFrictMat(
young=1e9,
poisson=0.25,
frictionAngle=
isCohesive=True,
normalCohesion
shearCohesion=1e6,
alphaKr=1.8,
alphaKtw=0,
etaRoll=1.8,
momentRotation
density=2600,
label='pol')
O.materials.
#######
O.dt=0.0001
#######
t1=polyhedron(
t2=polyhedron(
O.bodies.append(t1)
O.bodies.append(t2)
#######
O.bodies.
#######
O.engines=[
ForceResetter(),
InsertionSor
################## this part seems to have a problem ##################
InteractionLoop
[Ip2_
[Law2_ScGeom6D_
#######
),
NewtonIntegr
]
#######
from yade import qt
qt.Controller()
V = qt.View()
According to the code, I expect that the upper block slides on the bottom block by decreasing cohesion in contact face, but I got the following error:
Segmentation fault (core dumped)
I think the cohesion parameter isn't work correctly.Is this parameter defined in polyhedron materials? and as second question about this, how the contact law works between polyhedrons?
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 alma for more information if necessary.