yade-dev team mailing list archive

[noreply@xxxxxxxxxxxxx]

------------------------------------------------------------
revno: 4140
committer: Anton Gladky <gladky.anton@xxxxxxxxx>
timestamp: Tue 2014-08-12 21:46:58 +0200
message:
  Add checkClumpHopper autotest.
added:
  examples/clumps/clump-hopper-test.py
  scripts/checks-and-tests/checks/checkClumpHopper.py


--
lp:yade
https://code.launchpad.net/~yade-pkg/yade/git-trunk

Your team Yade developers is subscribed to branch lp:yade.
To unsubscribe from this branch go to https://code.launchpad.net/~yade-pkg/yade/git-trunk/+edit-subscription

=== added file 'examples/clumps/clump-hopper-test.py'
--- examples/clumps/clump-hopper-test.py	1970-01-01 00:00:00 +0000
+++ examples/clumps/clump-hopper-test.py	2014-08-12 19:46:58 +0000
@@ -0,0 +1,79 @@
+# -*- coding: utf-8
+
+from yade import pack,export,qt,geom
+import itertools
+from numpy import *
+
+
+kinEnergyMax = 100000
+# Parameters
+tc=0.001# collision time 
+en=.3  # normal restitution coefficient
+es=.3  # tangential restitution coefficient
+frictionAngle=radians(35)# 
+density=2700
+# facets material
+facetMat=O.materials.append(ViscElMat(frictionAngle=frictionAngle,tc=tc,en=en,et=es))
+# default spheres material
+dfltSpheresMat=O.materials.append(ViscElMat(density=density,frictionAngle=frictionAngle,tc=tc,en=en,et=es)) 
+
+O.dt=.05*tc # time step
+
+Rs=0.05 # particle radius
+
+# Create geometry
+
+x0=0.; y0=0.; z0=0.; ab=.7; at=2.; h=1.; hl=h; al=at*3
+
+zb=z0;   x0b=x0-ab/2.; y0b=y0-ab/2.; x1b=x0+ab/2.; y1b=y0+ab/2.
+zt=z0+h; x0t=x0-at/2.; y0t=y0-at/2.; x1t=x0+at/2.; y1t=y0+at/2.
+zl=z0-hl;x0l=x0-al/2.; y0l=y0-al/2.; x1l=x0+al/2.; y1l=y0+al/2.
+
+vibrationPlate = O.bodies.append(geom.facetBunker((x0,y0,z0),
+  dBunker=ab*4.0,dOutput=ab*1.7,hBunker=ab*3,hOutput=ab,hPipe=ab/3.0,
+  wallMask=4,segmentsNumber=10,material=facetMat))
+
+# Create engines
+O.engines=[
+	ForceResetter(),
+	InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
+	InteractionLoop(
+		[Ig2_Sphere_Sphere_ScGeom(), Ig2_Facet_Sphere_ScGeom()],
+		[Ip2_ViscElMat_ViscElMat_ViscElPhys()],
+		[Law2_ScGeom_ViscElPhys_Basic()],
+	),
+	DomainLimiter(lo=(-4,-4,-1),hi=(4,4,4),iterPeriod=5000,label='domLim'),
+	NewtonIntegrator(damping=0, gravity=[0,0,-9.81]),
+	PyRunner(iterPeriod=6500,command='addBodies()',nDo=7,label='addb'),
+	PyRunner(iterPeriod=1000,command='state()',label='state'),
+]
+
+numSphereGen = 0
+def addBodies():
+	global numSphereGen
+	# Create clumps...
+	clumpColor=(0.0, 0.5, 0.5)
+	for k,l in itertools.product(arange(0,10),arange(0,10)):
+		clpId,sphId=O.bodies.appendClumped([sphere(Vector3(x0t+Rs*(k*4+2),y0t+Rs*(l*4+2),i*Rs*2+zt+ab*3),Rs,color=clumpColor,material=dfltSpheresMat) for i in xrange(4)])
+		numSphereGen += len(sphId)
+
+	# ... and spheres
+	spheresColor=(0.4, 0.4, 0.4)
+	for k,l in itertools.product(arange(0,9),arange(0,9)):
+		sphAloneId=O.bodies.append( [sphere( Vector3(x0t+Rs*(k*4+4),y0t+Rs*(l*4+4),i*Rs*2.3+zt+ab*3),Rs,color=spheresColor,material=dfltSpheresMat) for i in xrange(4) ] )
+		numSphereGen += len(sphAloneId)
+
+def state():
+	global numSphereGen
+	print "Iter %d: Total number of generated spheres %d, removed particles %d, current particles %d, kinEnergy %g"%(O.iter, numSphereGen, domLim.nDeleted, numSphereGen-domLim.nDeleted, utils.kineticEnergy())
+	if (utils.kineticEnergy() > kinEnergyMax):
+		print "Kinetic energy is over a threshold value! Error!"
+
+addBodies()
+
+
+from yade import qt
+qt.View()
+O.saveTmp()
+
+#O.run(50001)

=== added file 'scripts/checks-and-tests/checks/checkClumpHopper.py'
--- scripts/checks-and-tests/checks/checkClumpHopper.py	1970-01-01 00:00:00 +0000
+++ scripts/checks-and-tests/checks/checkClumpHopper.py	2014-08-12 19:46:58 +0000
@@ -0,0 +1,77 @@
+# -*- coding: utf-8
+# This script is used in --check to keep the functionality of adding/removal
+# of particles and clumps functional
+
+import itertools
+from numpy import *
+
+
+kinEnergyMax = 100000
+# Parameters
+tc=0.001# collision time 
+en=.3  # normal restitution coefficient
+es=.3  # tangential restitution coefficient
+frictionAngle=radians(35)# 
+density=2700
+# facets material
+facetMat=O.materials.append(ViscElMat(frictionAngle=frictionAngle,tc=tc,en=en,et=es))
+# default spheres material
+dfltSpheresMat=O.materials.append(ViscElMat(density=density,frictionAngle=frictionAngle,tc=tc,en=en,et=es)) 
+
+O.dt=.05*tc # time step
+
+Rs=0.05 # particle radius
+
+# Create geometry
+
+x0=0.; y0=0.; z0=0.; ab=.7; at=2.; h=1.; hl=h; al=at*3
+
+zb=z0;   x0b=x0-ab/2.; y0b=y0-ab/2.; x1b=x0+ab/2.; y1b=y0+ab/2.
+zt=z0+h; x0t=x0-at/2.; y0t=y0-at/2.; x1t=x0+at/2.; y1t=y0+at/2.
+zl=z0-hl;x0l=x0-al/2.; y0l=y0-al/2.; x1l=x0+al/2.; y1l=y0+al/2.
+
+vibrationPlate = O.bodies.append(geom.facetBunker((x0,y0,z0),
+  dBunker=ab*4.0,dOutput=ab*1.7,hBunker=ab*3,hOutput=ab,hPipe=ab/3.0,
+  wallMask=4,segmentsNumber=10,material=facetMat))
+
+# Create engines
+O.engines=[
+	ForceResetter(),
+	InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
+	InteractionLoop(
+		[Ig2_Sphere_Sphere_ScGeom(), Ig2_Facet_Sphere_ScGeom()],
+		[Ip2_ViscElMat_ViscElMat_ViscElPhys()],
+		[Law2_ScGeom_ViscElPhys_Basic()],
+	),
+	DomainLimiter(lo=(-4,-4,-1),hi=(4,4,4),iterPeriod=5000,label='domLim'),
+	NewtonIntegrator(damping=0, gravity=[0,0,-9.81]),
+	PyRunner(iterPeriod=6500,command='addBodies()',nDo=7,label='addb'),
+	PyRunner(iterPeriod=10000,command='state()',label='state'),
+]
+
+numSphereGen = 0
+def addBodies():
+	global numSphereGen
+	# Create clumps...
+	clumpColor=(0.0, 0.5, 0.5)
+	for k,l in itertools.product(arange(0,10),arange(0,10)):
+		clpId,sphId=O.bodies.appendClumped([sphere(Vector3(x0t+Rs*(k*4+2),y0t+Rs*(l*4+2),i*Rs*2+zt+ab*3),Rs,color=clumpColor,material=dfltSpheresMat) for i in xrange(4)])
+		numSphereGen += len(sphId)
+
+	# ... and spheres
+	spheresColor=(0.4, 0.4, 0.4)
+	for k,l in itertools.product(arange(0,9),arange(0,9)):
+		sphAloneId=O.bodies.append( [sphere( Vector3(x0t+Rs*(k*4+4),y0t+Rs*(l*4+4),i*Rs*2.3+zt+ab*3),Rs,color=spheresColor,material=dfltSpheresMat) for i in xrange(4) ] )
+		numSphereGen += len(sphAloneId)
+
+def state():
+	global numSphereGen
+	global resultStatus
+	print "Iter %d: Total number of generated spheres %d, removed particles %d, current particles %d, kinEnergy %g"%(O.iter, numSphereGen, domLim.nDeleted, numSphereGen-domLim.nDeleted, utils.kineticEnergy())
+	if (utils.kineticEnergy() > kinEnergyMax):
+		print "Kinetic energy is over a threshold value! Error!"
+		resultStatus += 1
+		O.pause()
+
+addBodies()
+O.run(50002, True)