yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #17104
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.