yade-users team mailing list archive

[Arthur Mercier <question154219@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi,

I use Yade Bzr2811 and I'd like to model a concrete beam. I tried to put in the CpmMat and CpmPhys. There is no problem when I compil, I mean no error occurs. But when I try to run it, Yade closes automatically. Here is my script, it works with Frictmat and FrictPhys in comment:

"import numpy, plot, pack, geom
rad,num=.1,15; dist=1.9999*rad


# create rectangular box from facets
O.bodies.append(geom.facetBox((.5,-.5,.25),(.5,.5,.5),wallMask=63))
O.bodies.append(geom.facetBox((5,-.5,.25),(.5,.5,.5),wallMask=63))
# one arm
O.bodies.append([utils.sphere((x,.1,0),rad,wire=True) for x in numpy.arange(0,55*rad,dist)])

O.dt=.03*utils.PWaveTimeStep() 

O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
	InteractionLoop(
		#[Ig2_Sphere_Sphere_L6Geom(distFactor=-1),Ig2_Facet_Sphere_L3Geom()],
		#[Ip2_FrictMat_FrictMat_FrictPhys()],
		#[Law2_L6Geom_FrictPhys_Linear(charLen=1),Law2_L3Geom_FrictPhys_ElPerfPl()]),
		[Ig2_Sphere_Sphere_Dem3DofGeom(),Ig2_Facet_Sphere_Dem3DofGeom()],
		[Ip2_CpmMat_CpmMat_CpmPhys()],
		[Law2_Dem3DofGeom_CpmPhys_Cpm()]),	
	GravityEngine(gravity=(0,-20,0)),
	ForceEngine(force=(0,-5000,0),ids=[38]),
	NewtonIntegrator(damping=0.1),
]

O.saveTmp()"

Where is my mistake?

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.