yade-users team mailing list archive

[Oguz Cebeci <question241338@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi,
i´m building up a simulation. I want to run about a day. But i have one problem. I get a segmentation fault.
I think it happens, when some colision appear. especially when spheres collider with Facets.

Here is my code. Yade version 1.05

Did i miss something?
Best Regards


from yade import pack
from yade import ymport


facetMat=O.materials.append(FrictMat(density=7850,young=5e5,poisson=.25))
sphereMat=O.materials.append(FrictMat(density=7850,young=5e5,poisson=.25))
fctIds=O.bodies.append(ymport.stl('bottich_triangles_unit_meter-binary.stl',color=(1,1,1),wire=True,fixed=True,material=facetMat))
fctIds2=O.bodies.append(ymport.stl('blisk_triangles_unit_meter-binary.stl',color=(1,1,1),wire=True,fixed=True,material=facetMat))

amplitude=((2*pi)/360)*.75



sp=pack.SpherePack()
sp.makeCloud((-0.2,-.2,.1),(.2,.2,.25),rMean=.01,rRelFuzz=1,num=100)
sp.toSimulation()


O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
   InteractionLoop(
      # handle sphere+sphere and facet+sphere collisions
      [Ig2_Sphere_Sphere_L3Geom(),Ig2_Facet_Sphere_L3Geom()],
      [Ip2_FrictMat_FrictMat_FrictPhys()],
      [Law2_L3Geom_FrictPhys_ElPerfPl()]
   ),
   NewtonIntegrator(gravity=(0,0,-9.81),damping=0.75),
   HarmonicRotationEngine(ids=fctIds+fctIds2,rotationAxis=[0,1,0],zeroPoint=(0,0,0),rotateAroundZero=True,A=amplitude,f=30,dead=False,label='auslenkung',fi=pi/2),

   CombinedKinematicEngine(ids=fctIds+fctIds2,label='comb',dead=True) + HarmonicRotationEngine(ids=fctIds,rotationAxis=[1,0,0],zeroPoint=(0,0,0),rotateAroundZero=True,A=amplitude,f=30,dead=False,label='xRot',fi=0) 
+ HarmonicRotationEngine(ids=fctIds,rotationAxis=[0,1,0],zeroPoint=(0,0,0),rotateAroundZero=True,A=amplitude,f=30,dead=False,label='yRot',fi=pi/2),

   
   VTKRecorder(fileName='paraview/3d-vtk-',recorders=['all'],virtPeriod=0.01),
   PyRunner(command='finish()',realPeriod=9000),
   PyRunner(command='status()',iterPeriod=1000),
   PyRunner(command='finish_auslenkung()',virtPeriod=0.0083, nDo=2),

]
O.dt=5*utils.PWaveTimeStep()




def status():
	print 'iteration:' + str(O.iter)

def finish():
	print 'finished'
	O.pause()
	

def finish_auslenkung():
	print 'auslenkung'	
	auslenkung.dead=True
	comb.dead=False
	




print 'run'
O.run()




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