Title : Error message : can only concatenate list (not "RotationEngine") to list
Hello everyone,
I am using a program in which particles rotate in a drum (by the rotation of the drum).
I would like to stop the rotation of the drum at a certain point and rotate with respect to the rotation axis of the drum , only the set of all the particles by a certain angle dtetab for a single time step, in order to get the image of the flow after this rotation( angle dtetab).
The code :
from __future__ import print_function
from yade import qt
from yade import plot
from yade import utils, plot, qt
#******
#Unité de distance
#Diamètre d'une particule en m
d= 0.008
#Angular Velocity (rd/s)
angVel = 40* (2 * pi / 60)
# Diamètre du tambour
rD = 0.5*(50*d)
# Diamètre du cylindre intérieur
rcylint = 0.5*(10*d)
# Profondeur de la cellule
rP = 10*d
# hauteur couche supérieure
alpha=0
rh=alpha*
#Coefficient de restitution
e= 15/16
#Period (s)
TangVel = 2*pi/angVel
#******
O.load(
#******
id_SphereMat= O.materials.
SphereMat = O.materials[
id_BMat = O.materials.
BMat = O.materials[
def unitSquare():
"""Return square composed of 2 facets in the xy plane, centered at zero with unit extents."""
import gts
vv = [gts.Vertex(1,1,0), gts.Vertex(1,-1, 0), gts.Vertex(
ee = [gts.Edge(vv[0], vv[1]), gts.Edge(vv[1], vv[2]), gts.Edge(vv[2], vv[3]), gts.Edge(vv[3], vv[0]), gts.Edge(vv[0], vv[2])]
surf = gts.Surface()
surf.
surf.
return surf
def unitCylinder(
import numpy
from yade import pack
thetas = numpy.linspace(0, 2 * pi, nDiv, endpoint=True)
ptsBase = [Vector3(
ptsTop = [p + Vector3(0,0,2) for p in ptsBase]
return pack.sweptPolyl
from yade import pack, timing
cyl = unitCylinder()
sq1 = unitSquare()
sq1.translate(
cyl.copy(sq1)
sq2 = unitSquare()
sq2.translate(
cyl.copy(sq2)
cyl.scale(rD,rD,rP)
cylint = unitCylinder()
cylint.
cylIds = O.bodies.
cylintIds = O.bodies.
O.dt = 0.95*PWaveTimeS
itround=
O.engines = [
]),
]
pi=3.141592654
dtetab=pi/4
wtetab=dtetab/O.dt
def checkUnbalanced():
if O.iter<1*itround :
return
else:
for b in O.bodies:
if isinstance(
yade.qt.View()
O.run()
I get this error : can only concatenate list (not "RotationEngine") to list
Can you please help me to solve this error or suggest me another method ? Thanks for your help.
pegbehe
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Karol Brzezinski
- Solved:
- Last query:
- Last reply: