Re: [Yade-users] [Question #669371]: Heat exchange
Question #669371 on Yade changed: https://answers.launchpad.net/yade/+question/669371 Ed Carlos Alves Rocha posted a new comment: Thank you Robert, for me it is a great challeng. I will read more information and do something and share with the yade community if I can do anything. Anyone know about any article to get started? Regards, Ed -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
[Yade-users] [Question #669371]: Heat exchange
New question #669371 on Yade: https://answers.launchpad.net/yade/+question/669371 Hello everyone! I have a problem, I would like to do one simulation with heat exchange. The simulation it's about a hot wall, heated wall with resistance/ induction, where the particles heat exchange considering only the conduction. 1) How to give the input of the coefficients of conductivity and specific heat of the particles? 2) How can I give the input of the wall temperature? 3) How to activate heat exchange? -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #652297]: Feng Chen's Yade-OpenFoam coupling model.
Question #652297 on Yade changed: https://answers.launchpad.net/yade/+question/652297 Ed Carlos Alves Rocha posted a new comment: Thanks Robert Caulk for the help. -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #652297]: Feng Chen's Yade-OpenFoam coupling model.
Question #652297 on Yade changed: https://answers.launchpad.net/yade/+question/652297 Ed Carlos Alves Rocha posted a new comment: Hi Feng Chen, Thanks for your answer. I'll try both options. Can I contact you after, by email, to describe the simulation of my master degree and How I can to execute the CFD-DEM coupling? Regards, ED -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
[Yade-users] [Question #652297]: Feng Chen's Yade-OpenFoam coupling model.
New question #652297 on Yade: https://answers.launchpad.net/yade/+question/652297 Hello everyone! In Feng Chen's Yade-OpenFoam coupling model how can I start the Yade loop with OpenFoam? How do I write the FluidDragForceEngine function / engine in the simulation python file on the CFD-DEM coupling? (https://translate.googleusercontent.com/translate_c?act=url&depth=1&hl=pt-BR&ie=UTF8&prev=_t&rurl=translate.google.com.br&sl=en&sp=nmt4&tl=pt-BR&u=https://sites.google.com/site/fchen3/home/yade-openfoam-src&usg=ALkJrhhn1XzfWcXdprHY9pNuv7gs8cng-g) For example, I wrote a python file on yade that simulates the trajectory of a tennis ball (a spheres). Now I want to implement a wind flow transverse to the trajectory of that sphere. How do I call the FluidDragForceEngine function / engine and the OpenFoam? -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #165417]: Fluid Simulation
Question #165417 on Yade changed: https://answers.launchpad.net/yade/+question/165417 Ed Carlos Alves Rocha posted a new comment: Hi Feng Chen? How can I start the Yade Mechanical loop? How can I put FluidDragForceEngine in the simulation file on the yade in the CFD-DEM coupling? -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #647492]: YADE-OpenFOAM Coupling
Question #647492 on Yade changed: https://answers.launchpad.net/yade/+question/647492 Status: Answered => Solved Ed Carlos Alves Rocha confirmed that the question is solved: Thanks Chareyre, that solved my question. -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
[Yade-users] [Question #647492]: YADE-OpenFOAM Coupling
New question #647492 on Yade: https://answers.launchpad.net/yade/+question/647492 I have 2 doubts: 1) My first doubt: I would like export the results about the CFD simulation with OpenFoam to Yade Dem, I talk about the one way simulation, where the OpenFoam export the results to Yade Dem. Is it possible if I use the ParaView to export the results? Or I can export using the OpenFoam? 2) My Second doubt is about the 2 ways simulation, where we have the Yade - OpenFoam coupling model. How can I creat this simulation? Thanks, ED -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #631788]: Flow of sphere clump.
Question #631788 on Yade changed: https://answers.launchpad.net/yade/+question/631788 Status: Answered => Solved Ed Carlos Alves Rocha confirmed that the question is solved: Hello Jan! You are a genius. It worked. Thank you so much. A big hug for you. -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #631788]: Flow of sphere clump.
Question #631788 on Yade changed: https://answers.launchpad.net/yade/+question/631788 Status: Answered => Open Ed Carlos Alves Rocha is still having a problem: Hello Jan! Thank you for your attention. Sorry, I still can not do it. I tried the options below, however, it did not work. # # def inject(): for k in arange(0,4): ss2=O.bodies.appendClumped([utils.sphere((r*i,0,k*4*r),r, material="cana") for i in range(1,9)]) cid,bids=O.bodies.append(ss2) v = (9.1924,9.1924,0) #for s in ss2: s.state.vel = v O.bodies[cid].state.vel = v # # def inject(): for k in arange(0,4): ss2=O.bodies.appendClumped([utils.sphere((r*i,0,k*4*r),r, material="cana") for i in range(1,9)]) cid,bids=(ss2) v = (9.1924,9.1924,0) #for s in ss2: s.state.vel = v O.bodies[cid].state.vel = v # # def inject(): for k in arange(0,4): ss2=O.bodies.appendClumped([utils.sphere((r*i,0,k*4*r),r, material="cana") for i in range(1,9)]) cid,bids=(ss2) v = (9.1924,9.1924,0) for s in ss2: s.state.vel = v O.bodies[cid].state.vel = v # # Injection does not work. Can you help me again? -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #631788]: Flow of sphere clump.
Question #631788 on Yade changed: https://answers.launchpad.net/yade/+question/631788 Status: Solved => Open Ed Carlos Alves Rocha is still having a problem: Hello Jan! I created a sphere clump for the injection, following its tips. However, when I run the particle the error below appears, can you help me with this error? AttributeErrorTraceback (most recent call last) /usr/bin/yade in () /usr/bin/yade in inject() 27 v = (9.1924,9.1924,0) 28 for s in ss2: ---> 29s.state.vel = v 30 O.bodies[cid].state.vel = v 31 AttributeError: 'int' object has no attribute 'state' # # idCana=O.materials.append(FrictMat(young=24.5e6,poisson=.3,frictionAngle=.6, density=4.5e2,label="cana")) idSteel=O.materials.append(FrictMat(young=210e9,poisson=.3,frictionAngle=.8,density=7.85e3,label="steel")) O.bodies.append(geom.facetBox((12,3,1.5),(13,4,3),wallMask=31,material="steel" )) r = 0.125 def inject(): for k in arange(0,4): ss2=O.bodies.appendClumped([sphere((r*i,0,k*4*r),r, material="cana") for i in range(1,9)]) cid,bids = ss2 v = (9.1924,9.1924,0) for s in ss2: s.state.vel = v O.bodies[cid].state.vel = v O.engines = [ ForceResetter(), # sphere, facet, wall InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb(),Bo1_Wall_Aabb()]), InteractionLoop( # the loading plate is a wall, we need to handle sphere+sphere, sphere+facet, sphere+wall [Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom()], [Ip2_FrictMat_FrictMat_FrictPhys()], [Law2_ScGeom_FrictPhys_CundallStrack()] ), # NewtonIntegrator(gravity=(0,0,9.81),damping=0), NewtonIntegrator(gravity=(0,-9.81,0),damping=0.2), PyRunner(iterPeriod=50,command='inject()',initRun=True), ] O.dt=(1e-6)*PWaveTimeStep() #O.dt =1e-7 O.step() # from yade import qt qt.Controller() qt.View() O.saveTmp() # -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #631788]: Flow of sphere clump.
Question #631788 on Yade changed: https://answers.launchpad.net/yade/+question/631788 Status: Answered => Solved Ed Carlos Alves Rocha confirmed that the question is solved: Thanks Jan Stránský, that solved my question. -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #631788]: Flow of sphere clump.
Question #631788 on Yade changed: https://answers.launchpad.net/yade/+question/631788 Ed Carlos Alves Rocha posted a new comment: Thank you Jan. This solves my problem. I wish you a great week. A big hug. cheers Ed -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
[Yade-users] [Question #631788]: Flow of sphere clump.
New question #631788 on Yade:
https://answers.launchpad.net/yade/+question/631788
Hello everyone! I need to create a sphere factory or injector, for example, a
flow of 4 sphere clump. How can I do this programming in YADE?
In Yade there is an example dspherefactorvv. But I do not know how to introduce
a clump of sphere.
"""Simple script which shows how to create an inlet.
"""
shotsId,steelId=O.materials.append([
FrictMat(young=50e9,density=6000,poisson=.2,label='shots'),
FrictMat(young=210e9,density=7800, poisson=.3,label='steel'),
])
## same as
#
# shotsId,steelId=O.materials.index('shots'),O.materials.index('steel')
#
O.bodies.append(geom.facetBox(center=(0,0,0),extents=(30e-3,30e-3,0),wallMask=32,wire=False,material='steel',color=(0,1,.3)))
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()],verletDist=.05*.29e-3),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_MindlinPhys(
# define restitution coefficients between different
pairs of material ids, see the functor's documentation for details
#en=MatchMaker(fallback='zero',matches=((steelId,shotsId,.4),(shotsId,shotsId,1)))
en=MatchMaker(matches=((steelId,shotsId,.4),(shotsId,shotsId,1)))
)],
[Law2_ScGeom_MindlinPhys_Mindlin(label='contactLaw')]
),
NewtonIntegrator(damping=0),
## CircularFactory: disk if length=0 or cylinder if length>0
#CircularFactory(maxParticles=1,radius=8e-3,length=16e-3,center=(0,-15e-3,15e-3),rMin=0.28e-3,rMax=0.29e-3,vMin=100,vMax=100,vAngle=0,massFlowRate=100./60,normal=(0,1.5,-1),label='factory',materialId=shotsId),
## BoxFactory: a line, plane or cuboid
BoxFactory(maxParticles=1,extents=(8e-3,8e-3,8e-3),center=(0,-15e-3,15e-3),rMin=0.28e-3,rMax=0.29e-3,vMin=100,vMax=100,vAngle=0,massFlowRate=100./60,normal=(0,1.5,-1),label='factory',materialId=shotsId),
DomainLimiter(lo=(-30e-3,-30e-3,0),hi=(30e-3,30e-3,60e-3),iterPeriod=200),
#VTKRecorder(recorders=['spheres','facets','velocity'],fileName='/tmp/nozzle-',iterPeriod=500),
# run this every once in a while, to finalize the simulation at some
point
PyRunner(iterPeriod=1,command='if
factory.numParticles>=factory.maxParticles: O.stopAtIter=O.iter+8000;
timing.stats()')
]
# the timestep must be smaller because of high linear velocities of particles
# we cannot use PWaveTimeStep directly, since there are no spheres generated yet
O.dt=SpherePWaveTimeStep(factory.rMin,O.materials[factory.materialId].density,O.materials[factory.materialId].young)
O.saveTmp()
#O.timingEnabled=True
from yade import timing
try:
from yade import qt
# setup 3d view
v=qt.View()
v.upVector=(0,0,1); v.viewDir=(-1,0,-.3); v.center(median=False)
except ImportError: pass
O.run()
--
You received this question notification because your team yade-users is
an answer contact for Yade.
___
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #629486]: Particle Cylinder
Question #629486 on Yade changed: https://answers.launchpad.net/yade/+question/629486 Status: Open => Solved Ed Carlos Alves Rocha confirmed that the question is solved: Thank you Klauss. I'm going to oppen new question. Yes a Clump. I am going to open a new question, ie, a sphere factory or injector, for example, a flow of 4 sphere clump. -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #629486]: Particle Cylinder
Question #629486 on Yade changed:
https://answers.launchpad.net/yade/+question/629486
Status: Solved => Open
Ed Carlos Alves Rocha is still having a problem:
I have another problem. I would like to construct a flow of coupled
sphere particles.. In Yade has the example of the spheres factory.
However, what is the mechanics of the yade where I can build an
injection of particles of coupled spheres?
shotsId,steelId=O.materials.append([
FrictMat(young=50e9,density=6000,poisson=.2,label='shots'),
FrictMat(young=210e9,density=7800, poisson=.3,label='steel'),
])
## same as
#
# shotsId,steelId=O.materials.index('shots'),O.materials.index('steel')
#
O.bodies.append(geom.facetBox(center=(0,0,0),extents=(30e-3,30e-3,0),wallMask=32,wire=False,material='steel',color=(0,1,.3)))
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()],verletDist=.05*.29e-3),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_MindlinPhys(
# define restitution coefficients between different
pairs of material ids, see the functor's documentation for details
#en=MatchMaker(fallback='zero',matches=((steelId,shotsId,.4),(shotsId,shotsId,1)))
en=MatchMaker(matches=((steelId,shotsId,.4),(shotsId,shotsId,1)))
)],
[Law2_ScGeom_MindlinPhys_Mindlin(label='contactLaw')]
),
NewtonIntegrator(damping=0),
## CircularFactory: disk if length=0 or cylinder if length>0
#CircularFactory(maxParticles=1,radius=8e-3,length=16e-3,center=(0,-15e-3,15e-3),rMin=0.28e-3,rMax=0.29e-3,vMin=100,vMax=100,vAngle=0,massFlowRate=100./60,normal=(0,1.5,-1),label='factory',materialId=shotsId),
## BoxFactory: a line, plane or cuboid
BoxFactory(maxParticles=1,extents=(8e-3,8e-3,8e-3),center=(0,-15e-3,15e-3),rMin=0.28e-3,rMax=0.29e-3,vMin=100,vMax=100,vAngle=0,massFlowRate=100./60,normal=(0,1.5,-1),label='factory',materialId=shotsId),
DomainLimiter(lo=(-30e-3,-30e-3,0),hi=(30e-3,30e-3,60e-3),iterPeriod=200),
#VTKRecorder(recorders=['spheres','facets','velocity'],fileName='/tmp/nozzle-',iterPeriod=500),
# run this every once in a while, to finalize the simulation at some
point
PyRunner(iterPeriod=1,command='if
factory.numParticles>=factory.maxParticles: O.stopAtIter=O.iter+8000;
timing.stats()')
]
# the timestep must be smaller because of high linear velocities of particles
# we cannot use PWaveTimeStep directly, since there are no spheres generated yet
O.dt=SpherePWaveTimeStep(factory.rMin,O.materials[factory.materialId].density,O.materials[factory.materialId].young)
O.saveTmp()
#O.timingEnabled=True
from yade import timing
try:
from yade import qt
# setup 3d view
v=qt.View()
v.upVector=(0,0,1); v.viewDir=(-1,0,-.3); v.center(median=False)
except ImportError: pass
O.run()
--
You received this question notification because your team yade-users is
an answer contact for Yade.
___
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #629486]: Particle Cylinder
Question #629486 on Yade changed: https://answers.launchpad.net/yade/+question/629486 Status: Open => Solved Ed Carlos Alves Rocha confirmed that the question is solved: Thanks Klaus Thoeni, that solved my question. -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #629486]: Particle Cylinder
Question #629486 on Yade changed:
https://answers.launchpad.net/yade/+question/629486
Status: Solved => Open
Ed Carlos Alves Rocha is still having a problem:
I have another problem. I would like to construct a flow of coupled sphere
particles.. In Yade has the example of the spheres factory. However, what is
the mechanics of the yade where I can build an injection of particles of
coupled spheres?
shotsId,steelId=O.materials.append([
FrictMat(young=50e9,density=6000,poisson=.2,label='shots'),
FrictMat(young=210e9,density=7800, poisson=.3,label='steel'),
])
## same as
#
# shotsId,steelId=O.materials.index('shots'),O.materials.index('steel')
#
O.bodies.append(geom.facetBox(center=(0,0,0),extents=(30e-3,30e-3,0),wallMask=32,wire=False,material='steel',color=(0,1,.3)))
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()],verletDist=.05*.29e-3),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_MindlinPhys(
# define restitution coefficients between different
pairs of material ids, see the functor's documentation for details
#en=MatchMaker(fallback='zero',matches=((steelId,shotsId,.4),(shotsId,shotsId,1)))
en=MatchMaker(matches=((steelId,shotsId,.4),(shotsId,shotsId,1)))
)],
[Law2_ScGeom_MindlinPhys_Mindlin(label='contactLaw')]
),
NewtonIntegrator(damping=0),
## CircularFactory: disk if length=0 or cylinder if length>0
#CircularFactory(maxParticles=1,radius=8e-3,length=16e-3,center=(0,-15e-3,15e-3),rMin=0.28e-3,rMax=0.29e-3,vMin=100,vMax=100,vAngle=0,massFlowRate=100./60,normal=(0,1.5,-1),label='factory',materialId=shotsId),
## BoxFactory: a line, plane or cuboid
BoxFactory(maxParticles=1,extents=(8e-3,8e-3,8e-3),center=(0,-15e-3,15e-3),rMin=0.28e-3,rMax=0.29e-3,vMin=100,vMax=100,vAngle=0,massFlowRate=100./60,normal=(0,1.5,-1),label='factory',materialId=shotsId),
DomainLimiter(lo=(-30e-3,-30e-3,0),hi=(30e-3,30e-3,60e-3),iterPeriod=200),
#VTKRecorder(recorders=['spheres','facets','velocity'],fileName='/tmp/nozzle-',iterPeriod=500),
# run this every once in a while, to finalize the simulation at some
point
PyRunner(iterPeriod=1,command='if
factory.numParticles>=factory.maxParticles: O.stopAtIter=O.iter+8000;
timing.stats()')
]
# the timestep must be smaller because of high linear velocities of particles
# we cannot use PWaveTimeStep directly, since there are no spheres generated yet
O.dt=SpherePWaveTimeStep(factory.rMin,O.materials[factory.materialId].density,O.materials[factory.materialId].young)
O.saveTmp()
#O.timingEnabled=True
from yade import timing
try:
from yade import qt
# setup 3d view
v=qt.View()
v.upVector=(0,0,1); v.viewDir=(-1,0,-.3); v.center(median=False)
except ImportError: pass
O.run()
--
You received this question notification because your team yade-users is
an answer contact for Yade.
___
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #629486]: Particle Cylinder
Question #629486 on Yade changed: https://answers.launchpad.net/yade/+question/629486 Status: Answered => Solved Ed Carlos Alves Rocha confirmed that the question is solved: Thanks Jan Stránský, that solved my question. -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
Re: [Yade-users] [Question #629486]: Particle Cylinder
Question #629486 on Yade changed: https://answers.launchpad.net/yade/+question/629486 Status: Answered => Open Ed Carlos Alves Rocha is still having a problem: I did the Cylinder example, however, I do not know give the imput of collision detection for this example. Can you take a look and see what impact collision between the cylinders and the wall? Thank you very much for your attention. from yade.gridpfacet import * Parameter L=1.# length of the cylinder element r=0.1 # radius of the cylinder element phi=30. # friction angle E=1e6 # Young's modulus idCana=O.materials.append(FrictMat(young=3.4e6,poisson=.3,frictionAngle=.6, density=4.5e2,label="cana")) idSteel=O.materials.append(FrictMat(young=210e9,poisson=.25,frictionAngle=.8,density=7.85e3,label="steel")) O.bodies.append(geom.facetBox((0,0,0),(3,3,3),wallMask=31,material="steel" )) Engines O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_GridConnection_Aabb(),Bo1_Sphere_Aabb(),Bo1_Facet_Aabb(), Bo1_Wall_Aabb()]), InteractionLoop([ Ig2_GridNode_GridNode_GridNodeGeom6D(), Ig2_GridConnection_GridConnection_GridCoGridCoGeom(), Ig2_Sphere_Sphere_ScGeom(), Ig2_Facet_Sphere_ScGeom(), Ig2_Wall_Sphere_ScGeom() ], [ Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True,setCohesionOnNewContacts=True), Ip2_FrictMat_FrictMat_FrictPhys()], [ Law2_ScGeom6D_CohFrictPhys_CohesionMoment(), Law2_ScGeom_FrictPhys_CundallStrack(), Law2_ScGridCoGeom_FrictPhys_CundallStrack(), Law2_GridCoGridCoGeom_FrictPhys_CundallStrack() ]), NewtonIntegrator(gravity=(-0.,0,-10),damping=0.5,label='newton'), ] 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 Create cylinders nodesIds=[] cylIds=[] cylinder((0,0,0),(L,0,0),radius=r,nodesIds=nodesIds,cylIds=cylIds, fixed=False,color=[1,0,0],intMaterial='cMat',extMaterial='fMat') cylinder((L/4,2*L/3,L),(L/4,-L/3,L),radius=r,nodesIds=nodesIds,cylIds=cylIds, fixed=False,color=[0,1,0],intMaterial='cMat',extMaterial='fMat') cylinder((0,2*L/3,L),(0,-L/3,L),radius=r,nodesIds=nodesIds,cylIds=cylIds, fixed=False,color=[0,1,0],intMaterial='cMat',extMaterial='fMat') cylinder((L/2,L/2,L),(L/2,-L/2,L),radius=r,nodesIds=nodesIds,cylIds=cylIds, fixed=False,color=[0,1,1],intMaterial='cMat',extMaterial='fMat') cylinder((3*L/4,L/3,L),(3*L/4,-2*L/3,L),radius=r,nodesIds=nodesIds,cylIds=cylIds, fixed=False,color=[0,0,1],intMaterial='cMat',extMaterial='fMat') cylinder((L,L/3,L),(L,-2*L/3,L),radius=r,nodesIds=nodesIds,cylIds=cylIds, fixed=False,color=[0,0,1],intMaterial='cMat',extMaterial='fMat') For viewing from yade import qt qt.View() Gl1_Sphere.stripes=True Set a time step O.dt=1e-06 Allows to reload the simulation O.saveTmp() -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
[Yade-users] [Question #629486]: Particle Cylinder
New question #629486 on Yade: https://answers.launchpad.net/yade/+question/629486 I'm new to Yade. I have two particles to insert into a domain. One of them is cylindrical, how to construct for the cylindrical particles in the yade? Another Questio is How do I create a cloud os cylinders? Thanks -- You received this question notification because your team yade-users is an answer contact for Yade. ___ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
