Author: richefeu Date: 2009-07-02 15:09:34 +0200 (Thu, 02 Jul 2009) New Revision: 1832
Added: trunk/scripts/test/SpherePadder.py Log: Add an example of python script for SpherePadder Added: trunk/scripts/test/SpherePadder.py =================================================================== --- trunk/scripts/test/SpherePadder.py 2009-07-02 11:49:10 UTC (rev 1831) +++ trunk/scripts/test/SpherePadder.py 2009-07-02 13:09:34 UTC (rev 1832) @@ -0,0 +1,37 @@ +# encoding: utf-8 +# © 2009 Vincent Richefeu <[email protected]> +"""Simple example to show the SpherePadder usage in the built-in python interpreter of yade. +The sample is show with qt.View() but no computation can be done.""" + +# For the moment, the module packing (ie. file packing.so) +# must copyed manually in the folder gui (of the dynamic libs tree) +import packing + +mesh = packing.TetraMesh() +padder = packing.SpherePadder() + +mesh.read_gmsh('bench_1000.msh') +padder.plugTetraMesh(mesh) + +padder.setRadiusRatio(4.0,0.125) +padder.setMaxOverlapRate(1.0e-4) +padder.setVirtualRadiusFactor(100.0) + +padder.pad_5() +padder.place_virtual_spheres() +padder.setMaxNumberOfSpheres(8500) +padder.densify() +padder.save_mgpost("mgp.out.001") + +from yade import utils +o=Omega() + +lst = packing.getSphereList(padder) +for Sph in lst: + o.bodies.append(utils.sphere(center=[Sph[0],Sph[1],Sph[2]],radius=Sph[3])) + +o.bodies.append(utils.box(center=[0.0,-0.1,0.0],extents=[3.0,0.05,3.0])) + +from yade import qt +qt.View() + _______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
