Hi,
I'm having a problem where the spherical particles simply pass through my
meshed facet geometry. It's as if the mesh was never converted to facets
Using GMSH, I drew up a vessel with impeller geometry & meshed it in 3D. I then
saved it as a mesh file under ASCII format & physical entity element tag. I
then imported it into my yade script using ymport & assigned it the material
'steel'. But when I ran the script, the particles generated simply fell through
the impeller & vessel portion of the geometry instead of bouncing off. The same
thing happens when I save the mesh using a partition element tag. I'm having
trouble pinpointing the problem, but having looked at some examples of
ymport.gsmh, perhaps I'm doing something wrong in gmsh
My script are as follows:
from yade import utils,pack,plot,timing,ymport
from numpy import linspace
#import granulator geometry & convert to facets
O.bodies.append(yade.ymport.gmsh(meshfile='try_1_mesh.mesh', shift=Vector3(0,
0, 0), scale=1.0, orientation=Quaternion((1, 0, 0), 0)))
#define sphere material & facet material
stuff=O.materials.append(FrictMat(young=3e10,poisson=.2,density=1000))
steel=O.materials.append(FrictMat(young=210e9,poisson=.25,frictionAngle=.8,density=8000))
#pack volume with spheres with mean radius 0.5 & rRelFuzz=.3
sp=pack.SpherePack()
sp.makeCloud(Vector3(-10,-10,31),Vector3(10,10,35),.5,rRelFuzz=.3)
O.bodies.append([utils.sphere(c,r,material=stuff) for c,r in sp])
#open file for writing plot data
u=open('hsg1_vel.txt','w')
#define functions used in commands
def addPlotData(): #position & velocity of
particle 5
yade.plot.addData({'t':O.time,'i':O.iter,'position':O.bodies[5].state.pos,'velocity':O.bodies[5].state.vel[1]})
#write plot data into open file
def writedata():
u.write(str(O.time))#plot time int text file
u.write(" %s"%str(O.bodies[5].state.vel[0])) #plot x-velocity of sphere
4 into text file
u.write(" %s"%str(O.bodies[5].state.vel[1])) #plot y-velocity of sphere
4 into text file
u.write(" %s"%str(O.bodies[5].state.vel[2])) #plot z-velocity of sphere
4 into text file
u.write("\n")
#define constants
intradius=1.2 #bounding interaction radius
dampingmag=0.4 #damping for NewtonIntegrator
O.dt=5e-6 #timestep size=5e-6s
#mesh to facet keyword parameter
kwMeshes={'color':[1,1,0],'wire':True,'dynamic':False,'material':steel}
#define engine components
O.engines=[
ForceResetter(),
BoundDispatcher([Bo1_Sphere_Aabb(aabbEnlargeFactor=intradius)]),
InsertionSortCollider(),
InteractionDispatchers(
[Ig2_Sphere_Sphere_Dem3DofGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_Dem3DofGeom_FrictPhys_Basic()], #include Law2?
),
GravityEngine(gravity=(0,0,-1e4)),
NewtonIntegrator(damping=dampingmag),
PeriodicPythonRunner(iterPeriod=3,command='addPlotData()'),
PeriodicPythonRunner(iterPeriod=3,command='writedata()'),
VTKRecorder(iterPeriod=100,recorders=['spheres','facets','colors'],fileName='/home/kenmendoza/yade_results/vtu/p1')
]
#launch simulation viewing window & player
from yade import qt
qt.Controller()
qt.View()
O.run(50000)
I'd appreciate any help/comments on this matter.
Thanks,
Ken
This email is confidential and may be privileged. If you are not the intended
recipient, please delete it and notify us immediately.
Please do not copy it or use it for any purpose, or disclose its contents to
any other person. Thank you.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
