yade-users team mailing list archive

Thread
Date

[Question #295169]: Problem with MatchMaker

To: yade-users@xxxxxxxxxxxxxxxxxxx
From: Jonathan Pergoli <question295169@xxxxxxxxxxxxxxxxxxxxx>
Date: Fri, 10 Jun 2016 21:27:40 -0000
Reply-to: question295169@xxxxxxxxxxxxxxxxxxxxx
Sender: bounces@xxxxxxxxxxxxx

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

Hi everybody!

I'm using YADE to simulate an impact between two body with differents material's properties. Untill few days ago my code worked perfectly. Now it gives me the following problem:

terminate called after throwing an instance of 'std::invalid_argument'
  what():  MatchMaker: no match for (1,0), and values required for algo computation 'avg' not specified.
Aborted (core dumped)

The code is:
# Material
mli=FrictMat(density=643,frictionAngle=0.1489,label="MLI")
mat1=O.materials.append(mli)
gravel=FrictMat(density=1700,frictionAngle=0.9188,label="gravel")
mat2=O.materials.append(gravel)

# Ground
s=ymport.textExt('rest',format='x_y_z_r')
sphere=O.bodies.append(s)
print len(O.bodies)

# Create a vector o spheres to eliminate those that have COG above the container
for i in O.bodies:
	if isinstance(i.shape,Sphere):
		print i.state.mass
		if i.state.pos[2]>.4:
			O.bodies.erase(i.id)
for i in O.bodies:
	if isinstance(i.shape,Sphere):
		if i.state.pos[2]<0:
			O.bodies.erase(i.id)
for i in O.bodies:
	if isinstance(i.shape,Sphere):
		x=i.state.pos[0]
		y=i.state.pos[1]
		r=m.sqrt(x**2+y**2)
		if r>.75:
			O.bodies.erase(i.id)
print len(O.bodies)

# Cylinder
h=.4
c=geom.facetCylinder((0,0,.2),radius=.75,height=h,segmentsNumber=100,wallMask=6,material="gravel")
O.bodies.append(c)

p=utils.polyhedron((v1,v2,v3,v4,v5,v6,v7,v8),fixed=False,color=(.6,.45,0),material="MLI",wire=False)

# Engines
O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb(),Bo1_Polyhedra_Aabb()]),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom(),Ig2_Sphere_Polyhedra_ScGeom()],
		[Ip2_FrictMat_FrictMat_MindlinPhys(en=MatchMaker(matches=((mat1,mat2,.6),(mat2,mat2,.55))),es=MatchMaker(matches=((mat1,mat2,.6),(mat2,mat2,.55))),krot=3)],
		[Law2_ScGeom_MindlinPhys_Mindlin(includeMoment=True)]
	),
	NewtonIntegrator(gravity=(0,0,-2.5e-4),damping=0.0)
]

O.trackEnergy=True
O.dt=.5*PWaveTimeStep()

I use yadedaily on Ubuntu 14.04 LTS server.

Cheers,

Jonathan


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