← Back to team overview

yade-users team mailing list archive

[Question #693640]: utils.wall acts like air...

 

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

hi,

I wanna generate a 'wall' object to act as infinite axes-aligned plane, there is a example 'wall.py' in example file which runs as what i thought, somehow it wont work as it should be when i add it to my code...wall acted as air and sphere just went through it. please help me with this problem. below is two scripts 'wall.py' the example and my script 'air-wall.py'.

#--------------------------------------------wall.py---------------------------------------
O.materials.append(FrictMat(young=30e9,density=1000,poisson=.2,frictionAngle=.5))
O.bodies.append([
	wall(1,axis=2,sense=-1),
	wall(-5,axis=0,sense=1),
	wall(1,axis=1),
	wall((1,0,0),0),
	sphere([0,0,0],.5),
	sphere([-4,-4,-3],.5)
])
Gl1_Wall(div=10)

from yade import qt
qt.Controller()
qt.View()


O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb(),Bo1_Wall_Aabb()]),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom()],
		[Ip2_FrictMat_FrictMat_FrictPhys()],
		[Law2_ScGeom_FrictPhys_CundallStrack()],
	),
	NewtonIntegrator(damping=0.01,gravity=[1e2,1e2,1e2]),
	]

O.dt=PWaveTimeStep()
O.save('/tmp/a.xml')
O.saveTmp()
#---------------------------------------------------------------------------------------------------------



#------------------------------------------------------air-wall.py-------------------------------------
#sphere
O.materials.append(FrictMat(young=4.0e6,poisson=.3,frictionAngle=30,density=2630e+6,label='sphereMat'))#for sphere
s=utils.sphere(center=[0,100e-3,100e-3],radius=20e-3)
O.bodies.append(s)
	        
#### Wall ####
O.materials.append(FrictMat(young=1.0e6,poisson=0.2,density=2.60e3,frictionAngle=20,label='wallmat'))
O.bodies.append(utils.wall(0,axis=2))

from yade import qt
qt.Controller()
qt.View()


O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb(),Bo1_Wall_Aabb()]),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom()],
		[Ip2_FrictMat_FrictMat_FrictPhys()],
		[Law2_ScGeom_FrictPhys_CundallStrack()],
	),
	NewtonIntegrator(damping=0.01,gravity=[0,0,-1e2]),
	]
#---------------------------------------------------------------------------------------------

thanks!!!!!!!

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