yade-dev team mailing list archive

[noreply@xxxxxxxxxxxxx]

------------------------------------------------------------
revno: 1934
committer: Bruno Chareyre <bchareyre@r1arduina>
branch nick: trunk
timestamp: Thu 2009-12-31 18:05:15 +0100
message:
  -Fix the loading path, this example shows how shear create rotations.
modified:
  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.

=== modified file 'scripts/test/periodic-simple-shear.py'
--- scripts/test/periodic-simple-shear.py	2009-12-31 16:34:47 +0000
+++ scripts/test/periodic-simple-shear.py	2009-12-31 17:05:15 +0000
@@ -32,7 +32,7 @@
 	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
+		O.cell.velGrad[1,2]=6.0
 		triax.stressMask=7
 		triax['goal']=[-1e4,-1e4,-1e4]
 		phase+=1
@@ -42,14 +42,17 @@
 		#O.pause()
 		
 O.dt=utils.PWaveTimeStep()
-O.run(10000);
+O.run(7000);
 qt.View()
+r=qt.Renderer()
+r['Background_color']=1,1,1
 O.wait()
 
 O.cell.velGrad[1,2]=0
-O.cell.velGrad[2,1]=-3
-O.run(10000);
+O.cell.velGrad[2,1]=-6
+O.run(5000);
+O.wait()
 
-O.cell.velGrad[1,2]=-3
-O.cell.velGrad[2,1]=3
-O.run(10000);
+O.cell.velGrad[1,2]=6
+O.cell.velGrad[2,1]=0
+O.run(5000);