------------------------------------------------------------ revno: 1933 committer: Bruno Chareyre <bchare...@r1arduina> branch nick: trunk timestamp: Thu 2009-12-31 17:34:47 +0100 message: A funny script, and no crash! added: scripts/test/periodic-simple-shear.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 'scripts/test/periodic-simple-shear.py' --- scripts/test/periodic-simple-shear.py 1970-01-01 00:00:00 +0000 +++ scripts/test/periodic-simple-shear.py 2009-12-31 16:34:47 +0000 @@ -0,0 +1,55 @@ +# coding: utf-8 +# 2009 © Václav Šmilauer <[email protected]> +"Test and demonstrate use of PeriTriaxController." +from yade import * +from yade import pack,log,qt +log.setLevel('PeriTriaxController',log.TRACE) +O.periodic=True +O.cell.refSize=Vector3(.1,.1,.1) +#O.cell.Hsize=Matrix3(0.1,0,0, 0,0.1,0, 0,0,0.1) +sp=pack.SpherePack() +radius=5e-3 +num=sp.makeCloud(Vector3().ZERO,O.cell.refSize,radius,.2,500,periodic=True) # min,max,radius,rRelFuzz,spheresInCell,periodic +O.bodies.append([utils.sphere(s[0],s[1]) for s in sp]) + + +O.engines=[ + ForceResetter(), + BoundDispatcher([Bo1_Sphere_Aabb()]), + InsertionSortCollider(nBins=5,sweepLength=.05*radius), + InteractionDispatchers( + [Ig2_Sphere_Sphere_Dem3DofGeom()], + [SimpleElasticRelationships()], + [Law2_Dem3Dof_Elastic_Elastic()] + ), + PeriTriaxController(dynCell=True,mass=0.2,maxUnbalanced=0.01,relStressTol=0.02,goal=[-1e4,-1e4,0],stressMask=3,globUpdate=5,maxStrainRate=[1.,1.,1.],doneHook='triaxDone()',label='triax'), + NewtonIntegrator(damping=.2), +] + +phase=0 +def triaxDone(): + global phase + if phase==0: + print 'Here we are: stress',triax['stress'],'strain',triax['strain'],'stiffness',triax['stiff'] + print 'Now shearing.' + O.cell.velGrad[1,2]=3.0 + triax.stressMask=7 + triax['goal']=[-1e4,-1e4,-1e4] + phase+=1 + elif phase==1: + print 'Here we are: stress',triax['stress'],'strain',triax['strain'],'stiffness',triax['stiff'] + #print 'Done, pausing now.' + #O.pause() + +O.dt=utils.PWaveTimeStep() +O.run(10000); +qt.View() +O.wait() + +O.cell.velGrad[1,2]=0 +O.cell.velGrad[2,1]=-3 +O.run(10000); + +O.cell.velGrad[1,2]=-3 +O.cell.velGrad[2,1]=3 +O.run(10000);
_______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
