← Back to team overview

yade-users team mailing list archive

[Question #434127]: problem in creating the contact law between two polyhedron blocks

 

New question #434127 on Yade:
https://answers.launchpad.net/yade/+question/434127

I want to test slide state between two rigid blocks. For this I tried to run the following simple code with Yade 2016.06a-24-0557faf~trusty and my code is:

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=radians(30),
 isCohesive=True,
 normalCohesion=1e6,
 shearCohesion=1e6,	
 alphaKr=1.8,
 alphaKtw=0,
 etaRoll=1.8,
 momentRotationLaw=True,
 density=2600,
 label='pol')

O.materials.append(matPol)
#############################################################################
O.dt=0.0001
#################### geometry #################################
t1=polyhedron([(3.0,3.0,0.881935),(3.0,3.0,0.0),(0.0,3.0,3.0),(0.0,3.0,0.0),(3.0,0.0,0.0),(0.0,0.0,0.0),(0.820037,3.0,3.0),(0.0,0.0,3.0),(1.01018,0.0,3.0),(3.0,0.0,1.06668)],wire=False,material='pol')
t2=polyhedron([(3.0,0.0,1.06668),(3.0,0.0,3.0),(1.01018,0.0,3.0),(3.0,3.0,3.0),(0.820037,3.0,3.0),(3.0,3.0,0.881935)],wire=False,material='pol')
O.bodies.append(t1)			
O.bodies.append(t2)
##################### bottom surface  ##################
O.bodies.append(utils.wall(0,axis=2,sense=1,material='pol'))  

   
##################### this part should be verified #################

O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Polyhedra_Aabb(),Bo1_Wall_Aabb()], verletDist=0.001),
   
##################  this part seems to have a problem  ##################
InteractionLoop([Ig2_Polyhedra_Polyhedra_PolyhedraGeom(),Ig2_Wall_Polyhedra_PolyhedraGeom()],	
		[Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(label="cohesiveIp")],
[Law2_ScGeom6D_CohFrictPhys_CohesionMoment(label="cohesiveLaw")]
###################################################################

   ),
   NewtonIntegrator(gravity=(0,0,-9.81),damping=0.5)
]

########################################################
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?


-- 
You received this question notification because your team yade-users is
an answer contact for Yade.