yade-users team mailing list archive

[Pengfei Tan <question686692@xxxxxxxxxxxxxxxxxxxxx>]

Question #686692 on Yade changed:
https://answers.launchpad.net/yade/+question/686692

    Status: Answered => Open

Pengfei Tan is still having a problem:
Thank you very much for your prompt reply, Jan.

1. For your questions:

(1) I already created the cylindrical particles. My code is as follows:


##############################################

# encoding: utf-8
"An example showing how to create cylinders with random length."

from yade.gridpfacet import *
from yade import pack, plot,qt
import gts, os.path, locale

locale.setlocale(locale.LC_ALL, 'en_US.UTF-8')   # Note: gts is locale-
dependent. If, for example, german locale is used, gts.read()-function
does not import floats normally

#### Parameters ####
Lmin=3. # minimum length
Lmax=6. # maximum length 
n=10		# number of cylinders in one direction
r=0.5		# radius of the cylinder element
phi=30. # friction angle
E=1e6   # Young's modulus
radius=1e-02
wire=False
fixed=False
z=-1.2
color=[0,0,1]
Ns=200
R=0.1 #radius of the cylinder
L=5 #length of the cylinder
dx=0.5
dy=0
nx=20
ny=1
z=2
frictionAngle1=radians(15)
frictionAngle2=radians(15)
frictionAngle3=radians(15)

#### Engines ####
O.engines=[
	ForceResetter(),
	InsertionSortCollider([
		Bo1_GridConnection_Aabb(),
		Bo1_PFacet_Aabb(), 
		Bo1_Sphere_Aabb(),
	]),
	InteractionLoop([
		Ig2_GridNode_GridNode_GridNodeGeom6D(),
		Ig2_Sphere_Sphere_ScGeom(),
		Ig2_Sphere_GridConnection_ScGridCoGeom(),
		Ig2_GridConnection_GridConnection_GridCoGridCoGeom(),	# cylinder-cylinder interaction
		Ig2_Sphere_PFacet_ScGridCoGeom(),	# needed for GridNode-pFacet interaction (why is this not included in Ig2_GridConnection_PFacet_ScGeom???)
		Ig2_GridConnection_PFacet_ScGeom(),	# Cylinder-pFcet interaction
	],
	[
		Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True,setCohesionOnNewContacts=False),
		Ip2_FrictMat_FrictMat_FrictPhys()
	],
	[
		Law2_ScGeom_FrictPhys_CundallStrack(),          # contact law for sphere-sphere
		Law2_ScGeom6D_CohFrictPhys_CohesionMoment(),	# contact law for "internal" cylider forces
		Law2_ScGridCoGeom_FrictPhys_CundallStrack(),	# contact law for Cylinder-pFacet interaction
		Law2_GridCoGridCoGeom_FrictPhys_CundallStrack()	# contact law for cylinder-cylinder interaction
	]
	),
	GlobalStiffnessTimeStepper(timestepSafetyCoefficient=0.1,label='ts'), 
	NewtonIntegrator(gravity=(0.,0,-10),damping=0.5,label='newton'),
	RotationEngine(angularVelocity=1,rotationAxis=(0,1,0),ids=[0,1,2])+TranslationEngine(translationAxis=[1,0,0],velocity=0.5,ids=[0,1,2]),
	#VTKRecorder(fileName='3d-vtk-',recorders=['all'],iterPeriod=5000),
	#PyRunner(command='addPlotData()',iterPeriod=100),

]
O.dt=0.1*PWaveTimeStep()

#### Creat materials ####
O.materials.append( CohFrictMat( young=E,poisson=0.3,density=1000,frictionAngle=radians(phi),normalCohesion=1e10,shearCohesion=1e10,momentRotationLaw=True,label='cMat' ) )  #material to create the gridConnections
O.materials.append( FrictMat( young=E,poisson=0.3,density=1000,frictionAngle=radians(phi),label='fMat' ) )  # material for general interactions

#cylindrical roller: the IDs of cylinder must be [0, 1, 2] which are used for translation and rotation engines.
nodesIds=[]
cylIds=[]
cylinder((-8,-8,1.5),(-8,8,1.5),radius=0.5,nodesIds=nodesIds,cylIds=cylIds,color=[1,0,0],fixed=False,intMaterial='cMat',extMaterial='fMat')

#### Spheres ####
O.materials.append(FrictMat(young=4.0e6,poisson=0.5,frictionAngle=frictionAngle1,density=1600,label='spheremat'))
sp=pack.SpherePack()
Ns=sp.makeCloud(Vector3(-5,-5,3.0),Vector3(+5,+5,+6.0),0.3,0.2,Ns,False,0.8)
O.bodies.append([sphere(center,rad,material='spheremat') for center,rad in sp])

nodesIds1=[]
cylIds1=[]
#### Create cylinders ####
for i in range(0,nx):
	x=-5+i*dx
	for j in range(0,ny):
		y=-2.5+j*dy
		cylinder((x,y+L,z),(x,y,z),radius=R,nodesIds=nodesIds1,cylIds=cylIds1,
					 color=[1,0,0],fixed=False,intMaterial='cMat',extMaterial='fMat')

# collect history of data which will be plotted
def addPlotData():
	# each item is given a names, by which it can be the unsed in plot.plots
	# the **O.energy converts dictionary-like O.energy to plot.addData arguments
	plot.addData(i=O.iter,t=O.time,)
	plot.saveDataTxt('bbb.txt.bz2')

####create Pfacet from gts file####
#O.materials.append(CohFrictMat(young=1e8,poisson=0.3,density=2650,frictionAngle=radians(20),normalCohesion=1e100,shearCohesion=1e100,momentRotationLaw=True,label='gridNodeMat'))
#O.materials.append(FrictMat(young=1e8,poisson=0.3,density=2650,frictionAngle=radians(20),label='pFacetMat'))
#nodesIds1,cylIds1,pfIds1 =  gtsPFacet('box.gts',shift=(3.,0.,1.2),scale=1,radius=radius,wire=wire,fixed=False,materialNodes='gridNodeMat',material='pFacetMat',color=color)

#### Creat ground with pFacets ####
color=[255./255.,102./255.,0./255.]
n0=O.bodies.append( gridNode([-10,-10,0],r,wire=False,fixed=True,material='cMat',color=color) )
n1=O.bodies.append( gridNode([10,-10,0],r,wire=False,fixed=True,material='cMat',color=color) )
n2=O.bodies.append( gridNode([10,10,0],r,wire=False,fixed=True,material='cMat',color=color) )
n3=O.bodies.append( gridNode([-10,10,0],r,wire=False,fixed=True,material='cMat',color=color) )
O.bodies.append( gridConnection(n0,n1,r,color=color,material='fMat') )
O.bodies.append( gridConnection(n1,n2,r,color=color,material='fMat') )
O.bodies.append( gridConnection(n2,n0,r,color=color,material='fMat') )
O.bodies.append( gridConnection(n2,n3,r,color=color,material='fMat') )
O.bodies.append( gridConnection(n3,n0,r,color=color,material='fMat') )
O.bodies.append( pfacet(n0,n1,n2,wire=False,material='fMat',color=color) )
O.bodies.append( pfacet(n0,n2,n3,wire=False,material='fMat',color=color) )

#### For viewing ####
qt.View()
Gl1_Sphere.stripes=True

#### Allows to reload the simulation ####
O.saveTmp()

#####################################################

(2). "random distribution" means the positions and orientations.


2. For your solutions:

Because my spherical particles and cylindrical particles have total
different shapes, how to replace the spheres with cylinders?


Thank you a lot!

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.