yade-dev team mailing list archive

[noreply@xxxxxxxxxxxxx]

------------------------------------------------------------
revno: 1933
committer: Bruno Chareyre <bchareyre@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 <eudoxos@xxxxxxxx>
+"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);