------------------------------------------------------------ revno: 3040 committer: Bruno Chareyre <[email protected]> branch nick: trunk timestamp: Tue 2012-03-06 22:13:31 +0100 message: - example script for large bodies in periodic BCs added: examples/periodicSandPile.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/periodicSandPile.py' --- examples/periodicSandPile.py 1970-01-01 00:00:00 +0000 +++ examples/periodicSandPile.py 2012-03-06 21:13:31 +0000 @@ -0,0 +1,51 @@ +#!/usr/bin/python +# -*- coding: utf-8 -*- +# © 2012 Bruno Chareyre <[email protected]> +"""Script showing how large bodies can be combined with periodic boundary conditions using InsertioSortCollider::allowBiggerThanPeriod=True (1).""" + +from yade import pack +from pylab import * +from yade import qt + +O.periodic=True +length=0.4 +height=0.6 +width=0.2 +thickness=0.01 + +O.cell.hSize=Matrix3(length, 0, 0, + 0 ,3.*height, 0, + 0, 0, width) + +O.materials.append(FrictMat(density=1,young=1e5,poisson=0.3,frictionAngle=radians(30),label='boxMat')) +lowBox = utils.box( center=(length/2.0,height-thickness/2.0,width/2.0), extents=(length*1000.0,thickness/2.0,width*1000.0) ,fixed=True,wire=False) +O.bodies.append(lowBox) + +radius=0.01 +O.materials.append(FrictMat(density=1000,young=1e4,poisson=0.3,frictionAngle=radians(30),label='sphereMat')) +sp=pack.SpherePack() +sp.makeCloud((0.*length,height+1.2*radius,0.25*width),(0.5*length,2*height-1.2*radius,0.75*width),-1,.2,2000,periodic=True) +O.bodies.append([utils.sphere(s[0],s[1],color=(0.6+0.15*rand(),0.5+0.15*rand(),0.15+0.15*rand())) for s in sp]) + +O.dt=0.2*utils.PWaveTimeStep() +O.usesTimeStepper=True +newton=NewtonIntegrator(damping=0.6,gravity=(0,-10,0)) + +O.engines=[ + ForceResetter(), + #(1) This is where we all big bodies, else it would crash due to the very large bottom box: + InsertionSortCollider([Bo1_Box_Aabb(),Bo1_Sphere_Aabb()],allowBiggerThanPeriod=True), + InteractionLoop( + [Ig2_Sphere_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()], + [Ip2_FrictMat_FrictMat_FrictPhys()], + [Law2_ScGeom_FrictPhys_CundallStrack()] + ), + GlobalStiffnessTimeStepper(timeStepUpdateInterval=100,timestepSafetyCoefficient=0.8,defaultDt=O.dt), + newton +] + +from yade import qt +v = qt.Controller() +Gl1_Sphere.stripes=1 + +
_______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
