------------------------------------------------------------ revno: 1743 committer: Václav Šmilauer <vac...@flux> branch nick: trunk timestamp: Fri 2009-08-28 00:37:02 +0200 message: Added script for http://www.youtube.com/watch?v=KUv26xlh89I added: examples/mill.py
-- lp:yade https://code.launchpad.net/~yade-dev/yade/trunk Your team Yade developers is subscribed to branch lp:yade. To unsubscribe from this branch go to https://code.launchpad.net/~yade-dev/yade/trunk/+edit-subscription.
=== added file 'examples/mill.py' --- examples/mill.py 1970-01-01 00:00:00 +0000 +++ examples/mill.py 2009-08-27 22:37:02 +0000 @@ -0,0 +1,74 @@ +# encoding: utf-8 +from yade import pack +from numpy import linspace +# geometry parameters +bumpNum=20 +bumpHt,bumpTipAngle=0.07,60*pi/180 +millRad,millDp=1,1 # radius and depth (cylinder length) of the mill +sphRad,sphRadFuzz=0.03,.8 # mean radius and relative fuzz of the radius (random, uniformly distributed between sphRad*(1-.5*sphRadFuzz)â¦sphRad*(1+.5*sphRadFuzz)) +dTheta=pi/24 # circle division angle + + + +### +### mill geometry (parameteric) +### +bumpPeri=2*bumpHt*tan(.5*bumpTipAngle) # length of a bump on the perimeter of the mill +bumpAngle=bumpPeri/millRad # angle of one bump from the axis of the mill +interBumpAngle=2*pi/bumpNum +bumpRad=millRad-bumpHt +pts=[]; thMin=0 +for i in range(0,bumpNum): + thMin+=interBumpAngle + thMax=thMin+interBumpAngle-bumpAngle + thTip=thMax+.5*bumpAngle + # the circular parts spanning from thMin to thMax + for th0 in linspace(thMin,thMax,interBumpAngle/dTheta,endpoint=True): + pts.append(Vector3(-.5*millDp,millRad*cos(th0),millRad*sin(th0))) + # tip of the bump + pts.append(Vector3(-.5*millDp,bumpRad*cos(thTip),bumpRad*sin(thTip))) +# close the curve +pts+=[pts[0]] +# make the second contour, just shifted by millDp; ppts contains both +ppts=[pts,[p+Vector3(millDp,0,0) for p in pts]] +mill=pack.sweptPolylines2gtsSurface(ppts,threshold=.01*min(dTheta*millRad,bumpHt))#,capStart=True,capEnd=True) +millIds=O.bodies.append(pack.gtsSurface2Facets(mill,color=(1,0,1))) # add triangles, save their ids +# make the caps less comfortably, but looking better as two triangle couples over the mill +mrs2=millRad*sqrt(2) +cap1,cap2=[Vector3(0,0,mrs2),Vector3(0,-mrs2,0),Vector3(0,0,-mrs2)],[Vector3(0,0,mrs2),Vector3(0,0,-mrs2),Vector3(0,mrs2,0)] # 2 triangles at every side +for xx in -.5*millDp,.5*millDp: millIds+=O.bodies.append([utils.facet([p+Vector3(xx,0,0) for p in cap1],color=(0,0,0)),utils.facet([p+Vector3(xx,0,0) for p in cap2],color=(0,0,0))]) + +# define domains for initial cloud of red and blue spheres +packHt=.8*millRad # size of the area +bboxes=[(Vector3(-.5*millDp,-.5*packHt,-.5*packHt),Vector3(.5*millDp,0,.5*packHt)),(Vector3(-.5*millDp,0,-.5*packHt),Vector3(.5*millDp,.5*packHt,.5*packHt))] +colors=(1,0,0),(0,0,1) +for i in (0,1): # red and blue spheres + sp=pack.SpherePack(); bb=bboxes[i]; vol=(bb[1][0]-bb[0][0])*(bb[1][1]-bb[0][1])*(bb[1][2]-bb[0][2]) + sp.makeCloud(bb[0],bb[1],sphRad,sphRadFuzz,int(.25*vol/((4./3)*pi*sphRad**3)),False) + O.bodies.append([utils.sphere(s[0],s[1],density=3000,color=colors[i]) for s in sp]) + +print "Numer of grains",len(O.bodies)-len(millIds) + +O.dt=utils.PWaveTimeStep() + +O.engines=[ + BexResetter(), + BoundingVolumeMetaEngine([InteractingSphere2AABB(),InteractingFacet2AABB(),MetaInteractingGeometry2AABB()]), + InsertionSortCollider(), + InteractionDispatchers( + [ef2_Sphere_Sphere_Dem3DofGeom(),ef2_Facet_Sphere_Dem3DofGeom()], + [SimpleElasticRelationships()], + [Law2_Dem3Dof_Elastic_Elastic()], + ), + GravityEngine(gravity=(0,0,-3e4)), # gravity artificially high, to make it faster going ;-) + RotationEngine(rotateAroundZero=True,zeroPoint=(0,0,0),rotationAxis=(1,0,0),angularVelocity=-120,subscribedBodies=millIds), + SnapshotEngine(iterPeriod=30,fileBase='/tmp/mill-',viewNo=0,label='snapshooter'), + NewtonsDampedLaw(damping=.3), +] + +O.saveTmp() +from yade import qt +v=qt.View() +v.eyePosition=(3,.8,.96); v.upVector=(-.4,-.4,.8); v.viewDir=(-.9,-.25,-.3); v.axes=True; v.sceneRadius=1.9 +O.run(10000); O.wait() +utils.encodeVideoFromFrames(snapshooter['savedSnapshots'],out='/tmp/mill.ogg',fps=30)
_______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
