← Back to team overview

yade-users team mailing list archive

Re: [Question #667746]: Segmentation fault (core dumped) with gmshPFacet

 

Question #667746 on Yade changed:
https://answers.launchpad.net/yade/+question/667746

Kun Zeng gave more information on the question:
Hi Jan,

I found for my question #3:

Here is a bug in /yade/ymport.py in line 379. It is now:

if returnConnectivityTable:
		return unvReader.facets, unvReader.nodes, unvReader.elements
	return facets    ####line 379

I think there it should be changed to 'return unvReader.facets' .

1). And can I use the following code to import my file into yade and
apply gravity on the facet? But when I add:

kwm={'color':[1,0,0],'wire':False,'material':0} and run the simulation,
it seems no gravity on the structure I imported.

######################
Here is my code:

from yade import ymport,export
#Add material
O.materials.append(FrictMat(young=10e9,poisson=.25,frictionAngle=0.5,density=1e3))

#Parameters, which will be passed into spheres and facets creators

kwm={'color':[1,0,0],'wire':False,'material':0}

f,n,e = ymport.unv('shell.unv',returnConnectivityTable=True,**kwm)
O.bodies.append(f)
vtk = export.VTKExporter('test')
#vtk.exportFacetsAsMesh(elements=e)

O.engines=[
	#SubdomainBalancer(colorize=True,initRun=True,iterPeriod=100),
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb(),Bo1_Wall_Aabb()],label='collider',verletDist=.1),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom()],
		[Ip2_FrictMat_FrictMat_FrictPhys()],
		[Law2_ScGeom_FrictPhys_CundallStrack()],
	),
	NewtonIntegrator(damping=.1,gravity=(1e-2,1e-2,-9.81)),
]

O.dt=0.0002

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