[Yade-users] [Question #665130]: Calculating mass or valume of every single fragments after crushing
New question #665130 on Yade: https://answers.launchpad.net/yade/+question/665130 Hello, I am working on modelling particle crushing using bonded particle model in YADE, and I am now little confused on calculating volume of each fragments after crushing. In my model, the 'JCFpmMat' was used to model the agglomerate. Actually, in my simulation, agglomerate can be crushed into several fragments but still exist one or two cohesive bonds connecting neighbouring fragments; thus I cannot determine volumes of them. Can you give me any suggestions on this issue? Thank you in advance, Huihuang -- 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 #648734]: DOFs of facetSphere
Question #648734 on Yade changed: https://answers.launchpad.net/yade/+question/648734 Status: Answered => Solved Huihuang Xia 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
[Yade-users] [Question #648734]: DOFs of facetSphere
New question #648734 on Yade: https://answers.launchpad.net/yade/+question/648734 Hi, I intend to move a facetSphere and get its velocity when it interacts with other elements in my simulation. But to my disappointment, its velocity keeps constant. I have tried to read its source code, and got that facetSphere have a default -- fixed=True. Consequently, I modify this 'fixed=True' to 'fixed=False', but the facetSphere disappears. Besides, I have also tried another method. As facetSphere is fixed, its blockedDOFs = 'xyzXYZ' . Thus, I try to modify this one to "blockedDOFs = 'xzXZ'" (because I intend to get the velocity of facetSphere in y axis). But the facetSphere also disappears in my simulation. Can you give me a hint or solution? Thanks in advance, Huihuang -- 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 #643660]: A simple question on interactionDetectionFactor
Question #643660 on Yade changed: https://answers.launchpad.net/yade/+question/643660 Huihuang Xia confirmed that the question is solved: Thanks Robert Caulk, 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 #643660]: A simple question on interactionDetectionFactor
Question #643660 on Yade changed: https://answers.launchpad.net/yade/+question/643660 Status: Answered => Solved Huihuang Xia confirmed that the question is solved: Thanks Luc and Robert, you solved my problem. -- 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 #643660]: A simple question on interactionDetectionFactor
New question #643660 on Yade: https://answers.launchpad.net/yade/+question/643660 Hi, https://yade-dem.org/doc/yade.wrapper.html?highlight=interaction%20detection#yade.wrapper.Ig2_Sphere_Sphere_ScGeom.interactionDetectionFactor Here gives a parameter that create interaction range for particle packing, but I have no idea on how to determine this parameter ? Because in this paper (A DEM model for soft and hard rocks: Role of grain interlocking on strength), it was stated that if Rmax/Rmin=2, interactionDetectionFactor must not be larger than 1.5. In my simulations, Rmax/Rmin=9, thus, I have no idea to determine the interactionDetectionFactor. I have tried interactionDetectionFactor=1.5 to do some simulations, it got a relatively good result. But this selection (interactionDetectionFactor=1.5) seems go against the theory put in that paper. Can anyone help me solving this confusion ? Thanks in advance, Huihuang -- 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 #632699]: A simple question on selection of alphaKr and alphaKtw
New question #632699 on Yade: https://answers.launchpad.net/yade/+question/632699 Hi, CohFrictMat was used to model breakage of brittle materials, but I have no idea on how to select two parameters: alphaKr and alphaKtw. In my previous simulations, I chose the default value (2.0) of alphaKr and alphaKtw. But the agglomerate created by bonding lots of cohesive sub-particles has a big plastic deformation rather than fragile breakage when this agglomerate impacting on a rigid wall. Can you give me a hint to solve this problem? Thanks in advance, Huihuang -- 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 #631227]: Porosity of a packing with overlapped particles
Question #631227 on Yade changed: https://answers.launchpad.net/yade/+question/631227 Status: Answered => Open Huihuang Xia is still having a problem: Hi Klaus, Thanks for your patience. I have selected two diverse cubic sub-volumes in my irregular packing, but they got a porosity of 1.0 Huihuang Xia -- 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 #631227]: Porosity of a packing with overlapped particles
Question #631227 on Yade changed: https://answers.launchpad.net/yade/+question/631227 Status: Answered => Open Huihuang Xia is still having a problem: Thanks Klaus, I have tried your method, this approach works really good. Here is another question, Can I use this method to calculate the porosity of an irregular packing? Choosing a sub-volume from a irregular packing can give a result as precise as the whole volume? FW your reply, regards, Huihuang -- 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 #631227]: Porosity of a packing with overlapped particles
New question #631227 on Yade: https://answers.launchpad.net/yade/+question/631227 Hi, I defined a function to calculate the porosity of a packing with lots of overlapped particles, because the build-in function "utils.porosity()" always got a negative one. This question (https://answers.launchpad.net/yade/+question/404363) gave me some useful hints. My own function has a serious error, that is my function always got a same porosity whatever the growFactor is. Here is my code segment: from yade import utils from yade import pack from yade import qt mat1 = O.materials.append(FrictMat(young=5.0e7,poisson=0.25,density=2640.0,frictionAngle=radians(18))) pred = pack.inAlignedBox((0,0,0),(0.01,0.01,0.01)) dim = pred.dim() packing = regularHexa(pred, radius=0.0005,gap=0,material=mat1) O.bodies.append(packing) # grow particles in order to get lower porosity growFactor = 1.5 utils.growParticles(growFactor) def overlappedPorosity(totVol,density): totVol = dim[0]*dim[1]*dim[2] mass = sum([b.state.mass for b in O.bodies])*pow(growFactor,-3) return (totVol - mass/density)/totVol print 'The porosity of this packing is:',overlappedPorosity(dim[0]*dim[1]*dim[2],O.materials[0].density) # 3D view and controller try: qt.Controller() qt.View() except: pass In my function, ''mass = sum([b.state.mass for b in O.bodies])*pow(growFactor,-3)'', because after the packing was generated, the radius of each particle in packing was magnified by growFactor using growParticles. Any useful advice will be highly appreciated. Thanks in advance, Huihuang Xia -- 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 #625502]: How to get the minimum porosity of packing using randomDensePack
Question #625502 on Yade changed: https://answers.launchpad.net/yade/+question/625502 Status: Answered => Solved Huihuang Xia confirmed that the question is solved: Thanks Tina, I will try your advice. Regards, Huihuang -- 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 #625502]: How to get the minimum porosity of packing using randomDensePack
Question #625502 on Yade changed: https://answers.launchpad.net/yade/+question/625502 Huihuang Xia confirmed that the question is solved: Thanks Tina Asia, 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 #625502]: How to get the minimum porosity of packing using randomDensePack
Question #625502 on Yade changed: https://answers.launchpad.net/yade/+question/625502 Status: Answered => Open Huihuang Xia is still having a problem: Thanks Robert, I added the following function at the bottom of session1 and using PyRunner call this function, but the exported text file contains nothing. Here is my code: # -*- coding: utf-8 -*- from yade import pack from yade import export ### DEFINING VARIABLES AND MATERIALS ### # The following 5 lines will be used later for batch execution nRead=readParamsFromTable( num_spheres=5000,# number of spheres compFricDegree = 30, # contact friction during the confining phase key='_triax_base_', # put you simulation's name here unknownOk=True ) from yade.params import table num_spheres=table.num_spheres# number of spheres key=table.key targetPorosity = 0.43 #the porosity we want for the packing compFricDegree = table.compFricDegree # initial contact friction during the confining phase (will be decreased during the REFD compaction process) finalFricDegree = 30 # contact friction during the deviatoric loading rate=-0.02 # loading rate (strain rate) damp=0.2 # damping coefficient stabilityThreshold=0.01 # we test unbalancedForce against this value in different loops (see below) young=5e6 # contact stiffness mn,mx=Vector3(0,0,0),Vector3(1,1,1) # corners of the initial packing ## create materials for spheres and plates O.materials.append(FrictMat(young=young,poisson=0.5,frictionAngle=radians(compFricDegree),density=2600,label='spheres')) O.materials.append(FrictMat(young=young,poisson=0.5,frictionAngle=0,density=0,label='walls')) ## create walls around the packing walls=aabbWalls([mn,mx],thickness=0,material='walls') wallIds=O.bodies.append(walls) ## use a SpherePack object to generate a random loose particles packing global sp sp=pack.SpherePack() sp.makeCloud(mn,mx,-1,0.,num_spheres,False, 0.95,seed=1) #"seed" make the "random" generation always the same sps=O.bodies.append([sphere(center,rad,material='spheres') for center,rad in sp]) #or alternatively (higher level function doing exactly the same): #sp.toSimulation(material='spheres') ### DEFINING ENGINES ### triax=TriaxialStressController( ## TriaxialStressController will be used to control stress and strain. It controls particles size and plates positions. ## this control of boundary conditions was used for instance in http://dx.doi.org/10.1016/j.ijengsci.2008.07.002 maxMultiplier=1.+2e4/young, # spheres growing factor (fast growth) finalMaxMultiplier=1.+2e3/young, # spheres growing factor (slow growth) thickness = 0, ## switch stress/strain control using a bitmask. What is a bitmask, huh?! ## Say x=1 if stess is controlled on x, else x=0. Same for for y and z, which are 1 or 0. ## Then an integer uniquely defining the combination of all these tests is: mask = x*1 + y*2 + z*4 ## to put it differently, the mask is the integer whose binary representation is xyz, i.e. ## "100" (1) means "x", "110" (3) means "x and y", "111" (7) means "x and y and z", etc. stressMask = 7, internalCompaction=True, # If true the confining pressure is generated by growing particles ) newton=NewtonIntegrator(damping=damp) O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb()]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()], [Ip2_FrictMat_FrictMat_FrictPhys()], [Law2_ScGeom_FrictPhys_CundallStrack()] ), PyRunner(command='stop()',iterPeriod=10), ## We will use the global stiffness of each body to determine an optimal timestep (see https://yade-dem.org/w/images/1/1b/Chareyre_licensed.pdf) GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.8), triax, TriaxialStateRecorder(iterPeriod=100,file='WallStresses'+table.key), newton ] #Display spheres with 2 colors for seeing rotations better Gl1_Sphere.stripes=0 if nRead==0: yade.qt.Controller(), yade.qt.View() ## UNCOMMENT THE FOLLOWING SECTIONS ONE BY ONE ## DEPENDING ON YOUR EDITOR, IT COULD BE DONE ## BY SELECTING THE CODE BLOCKS BETWEEN THE SUBTITLES ## AND PRESSING CTRL+SHIFT+D ### ### APPLYING CONFINING PRESSURE ### ### #the value of (isotropic) confining stress defines the target stress to be applied in all three directions triax.goal1=triax.goal2=triax.goal3=-1 def stop(): if O.iter>5000: O.pause() # filter this sp into a
Re: [Yade-users] [Question #625502]: How to get the minimum porosity of packing using randomDensePack
Question #625502 on Yade changed:
https://answers.launchpad.net/yade/+question/625502
Status: Needs information => Open
Huihuang Xia gave more information on the question:
Thanks Robert,
https://answers.launchpad.net/yade/+question/246971, gives a good
example. Thus, I want to do a tri-axial compression on a cubic sample
and filter a stone particle shape using pack.filterSpherePack from this
sample. Is this method correct?
According to
https://github.com/yade/trunk/blob/e4e757f2e98a620e3177b7a36a1d10f69f6a6a28/examples
/triax-tutorial/script-session1.py, here is a code was used to filter a
sphere from a cubic sample, but this code failed, because it got the
following error:
Desktop/sphere$ yade sphere.py
Welcome to Yade 2016.06a
TCP python prompt on localhost:9000, auth cookie `ecssdy'
XMLRPC info provider on http://localhost:21000
Running script sphere.py
[[ ^L clears screen, ^U kills line. F12 controller, F11 3d view (use h-key for
showing help), F10 both, F9 generator, F8 plot. ]]
Yade [1]:
---
ValueErrorTraceback (most recent call last)
/home/huihuang/YADE/install/lib/x86_64-linux-gnu/yade-2016.06a/py/yade/__init__.pyc
in refreshEvent(self)
190 def zxySlot(self): self.setViewAxes((0,-1,0),(1,0,0))
191 def refreshEvent(self):
--> 192 self.refreshValues()
193 self.activateControls()
194 def deactivateControls(self):
/home/huihuang/YADE/install/lib/x86_64-linux-gnu/yade-2016.06a/py/yade/__init__.pyc
in refreshValues(self)
263 self.iterLabel.setText('#%ld / %ld, %.1f/s
%s'%(O.iter,stopAtIter,self.iterPerSec,subStepInfo))
264 if t!=float('inf'):
--> 265 s=int(t); ms=int(t*1000)%1000;
us=int(t*100)%1000; ns=int(t*10)%1000
266
self.virtTimeLabel.setText(u'%03ds%03dm%03dμ%03dn'%(s,ms,us,ns))
267 else: self.virtTimeLabel.setText(u'[ ∞ ] ?!')
ValueError: cannot convert float NaN to integer
Here is my code:
from yade import pack
from yade import export
### DEFINING VARIABLES AND MATERIALS ###
# The following 5 lines will be used later for batch execution
nRead=readParamsFromTable(
num_spheres=3000,# number of spheres
compFricDegree = 30, # contact friction during the confining phase
key='_triax_base_', # put you simulation's name here
unknownOk=True
)
from yade.params import table
num_spheres=table.num_spheres# number of spheres
key=table.key
targetPorosity = 0.43 #the porosity we want for the packing
compFricDegree = table.compFricDegree # initial contact friction during the
confining phase (will be decreased during the REFD compaction process)
finalFricDegree = 30 # contact friction during the deviatoric loading
rate=-0.02 # loading rate (strain rate)
damp=0.2 # damping coefficient
stabilityThreshold=0.01 # we test unbalancedForce against this value in
different loops (see below)
young=5e6 # contact stiffness
mn,mx=Vector3(0,0,0),Vector3(1,1,1) # corners of the initial packing
## create materials for spheres and plates
O.materials.append(FrictMat(young=young,poisson=0.5,frictionAngle=radians(compFricDegree),density=2600,label='spheres'))
O.materials.append(FrictMat(young=young,poisson=0.5,frictionAngle=0,density=0,label='walls'))
## create walls around the packing
walls=aabbWalls([mn,mx],thickness=0,material='walls')
wallIds=O.bodies.append(walls)
## use a SpherePack object to generate a random loose particles packing
global sp
sp=pack.SpherePack()
sp.makeCloud(mn,mx,-1,0.,num_spheres,False, 0.95,seed=1) #"seed" make the
"random" generation always the same
O.bodies.append([sphere(center,rad,material='spheres') for center,rad in sp])
### DEFINING ENGINES ###
triax=TriaxialStressController(
## TriaxialStressController will be used to control stress and strain.
It controls particles size and plates positions.
## this control of boundary conditions was used for instance in
http://dx.doi.org/10.1016/j.ijengsci.2008.07.002
maxMultiplier=1.+2e4/young, # spheres growing factor (fast growth)
finalMaxMultiplier=1.+2e3/young, # spheres growing factor (slow growth)
thickness = 0,
## switch stress/strain control using a bitmask. What is a bitmask,
huh?!
## Say x=1 if stess is controlled on x, else x=0. Same for for y and z,
which are 1 or 0.
## Then an integer uniquely defining the combination of all these tests
is: mask = x*1 + y*2 + z*4
## to put it differently, the mask is the integer whose binary
representation is xyz, i.e.
## "100" (1) means "x", "110
Re: [Yade-users] [Question #625502]: How to get the minimum porosity of packing using randomDensePack
Question #625502 on Yade changed:
https://answers.launchpad.net/yade/+question/625502
Huihuang Xia posted a new comment:
Thanks Robert,
https://answers.launchpad.net/yade/+question/246971, gives a good
example. Thus, I want to do a tri-axial compression on a cubic sample
and filter a stone particle shape using pack.filterSpherePack from this
sample. Is this method correct?
According to
https://github.com/yade/trunk/blob/e4e757f2e98a620e3177b7a36a1d10f69f6a6a28/examples
/triax-tutorial/script-session1.py, here is a code was used to filter a
sphere from a cubic sample, but this code failed, because it got the
following error:
Desktop/sphere$ yade sphere.py
Welcome to Yade 2016.06a
TCP python prompt on localhost:9000, auth cookie `ecssdy'
XMLRPC info provider on http://localhost:21000
Running script sphere.py
[[ ^L clears screen, ^U kills line. F12 controller, F11 3d view (use h-key for
showing help), F10 both, F9 generator, F8 plot. ]]
Yade [1]:
---
ValueErrorTraceback (most recent call last)
/home/huihuang/YADE/install/lib/x86_64-linux-gnu/yade-2016.06a/py/yade/__init__.pyc
in refreshEvent(self)
190 def zxySlot(self): self.setViewAxes((0,-1,0),(1,0,0))
191 def refreshEvent(self):
--> 192 self.refreshValues()
193 self.activateControls()
194 def deactivateControls(self):
/home/huihuang/YADE/install/lib/x86_64-linux-gnu/yade-2016.06a/py/yade/__init__.pyc
in refreshValues(self)
263 self.iterLabel.setText('#%ld / %ld, %.1f/s
%s'%(O.iter,stopAtIter,self.iterPerSec,subStepInfo))
264 if t!=float('inf'):
--> 265 s=int(t); ms=int(t*1000)%1000;
us=int(t*100)%1000; ns=int(t*10)%1000
266
self.virtTimeLabel.setText(u'%03ds%03dm%03dμ%03dn'%(s,ms,us,ns))
267 else: self.virtTimeLabel.setText(u'[ ∞ ] ?!')
ValueError: cannot convert float NaN to integer
Here is my code:
from yade import pack
from yade import export
### DEFINING VARIABLES AND MATERIALS ###
# The following 5 lines will be used later for batch execution
nRead=readParamsFromTable(
num_spheres=3000,# number of spheres
compFricDegree = 30, # contact friction during the confining phase
key='_triax_base_', # put you simulation's name here
unknownOk=True
)
from yade.params import table
num_spheres=table.num_spheres# number of spheres
key=table.key
targetPorosity = 0.43 #the porosity we want for the packing
compFricDegree = table.compFricDegree # initial contact friction during the
confining phase (will be decreased during the REFD compaction process)
finalFricDegree = 30 # contact friction during the deviatoric loading
rate=-0.02 # loading rate (strain rate)
damp=0.2 # damping coefficient
stabilityThreshold=0.01 # we test unbalancedForce against this value in
different loops (see below)
young=5e6 # contact stiffness
mn,mx=Vector3(0,0,0),Vector3(1,1,1) # corners of the initial packing
## create materials for spheres and plates
O.materials.append(FrictMat(young=young,poisson=0.5,frictionAngle=radians(compFricDegree),density=2600,label='spheres'))
O.materials.append(FrictMat(young=young,poisson=0.5,frictionAngle=0,density=0,label='walls'))
## create walls around the packing
walls=aabbWalls([mn,mx],thickness=0,material='walls')
wallIds=O.bodies.append(walls)
## use a SpherePack object to generate a random loose particles packing
global sp
sp=pack.SpherePack()
sp.makeCloud(mn,mx,-1,0.,num_spheres,False, 0.95,seed=1) #"seed" make the
"random" generation always the same
O.bodies.append([sphere(center,rad,material='spheres') for center,rad in sp])
### DEFINING ENGINES ###
triax=TriaxialStressController(
## TriaxialStressController will be used to control stress and strain.
It controls particles size and plates positions.
## this control of boundary conditions was used for instance in
http://dx.doi.org/10.1016/j.ijengsci.2008.07.002
maxMultiplier=1.+2e4/young, # spheres growing factor (fast growth)
finalMaxMultiplier=1.+2e3/young, # spheres growing factor (slow growth)
thickness = 0,
## switch stress/strain control using a bitmask. What is a bitmask,
huh?!
## Say x=1 if stess is controlled on x, else x=0. Same for for y and z,
which are 1 or 0.
## Then an integer uniquely defining the combination of all these tests
is: mask = x*1 + y*2 + z*4
## to put it differently, the mask is the integer whose binary
representation is xyz, i.e.
## "100" (1) means "x", "110" (3) means "x and y", "111" (7) mea
Re: [Yade-users] [Question #625502]: How to get the minimum porosity of packing using randomDensePack
Question #625502 on Yade changed: https://answers.launchpad.net/yade/+question/625502 Status: Answered => Open Huihuang Xia is still having a problem: Hi Robert, Thanks for your reply, but in my simulation, the stone particle was scanned by a 3D laser scanner. https://github.com/yade/trunk/blob/e4e757f2e98a620e3177b7a36a1d10f69f6a6a28/examples /triax-tutorial/script-session1.py As for this example, I have tried this code, but all failed, because the stone particle is a irregular one. Huihuang Xia -- 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 #625502]: How to get the minimum porosity of packing using randomDensePack
New question #625502 on Yade: https://answers.launchpad.net/yade/+question/625502 Hi, I am working on a project modeling stone particle using cohesive elements, but in my simulation, randomDensePack always give a low porosity about 0.7. Thus, some numerical simulations using this low porosity packing can not get a satisfactory results. Can you tell me how to insert sub-particles to the profile of a stone particle as many as possible? Huihuang Xia -- 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 #624246]: Call yade-2017.01a after compilation
Question #624246 on Yade changed: https://answers.launchpad.net/yade/+question/624246 Huihuang Xia posted a new comment: Hi Mike, Please open this file folder, myYade/install/bin, in 'bin', can you see a file named 'yade-2017.01a'? If yes, please move to the next step: 2. open a new teminal, and input ' gedit ~/.bashrc' ( there is a blank between t and ~). 3. Copy this line 'alias yade='python SHOME/myYade/install/bin/yade-2017.01a' and place it at the bottom of bashrc, save and exit it. Huihuang Xia -- 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 #611087]: I got a wrong damage ratio
Question #611087 on Yade changed: https://answers.launchpad.net/yade/+question/611087 Status: Answered => Solved Huihuang Xia confirmed that the question is solved: Thanks Luc Scholtès, 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 #611087]: I got a wrong damage ratio
Question #611087 on Yade changed: https://answers.launchpad.net/yade/+question/611087 Status: Answered => Open Huihuang Xia is still having a problem: Hi Robert, Thanks for your patience, but your method also got a wrong result, as damage ratio was greater than 1 and number of broken bonds could not increase gradually, which violated the definition of damage ratio. -- 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 #611087]: I got a wrong damage ratio
Question #611087 on Yade changed: https://answers.launchpad.net/yade/+question/611087 Status: Answered => Open Huihuang Xia is still having a problem: Hi Robert, Thanks for your answer, but put the numBroCohBonds=0 inside the damageRatio() function is reallly a big mistake. You can refer to https://answers.launchpad.net/yade/+question/432617, in this question, it was stated that put the numBroCohBonds=0 inside the damageRatio() function led to a initialization of numBroCohBonds to zero each time. Damage ratio is the fraction of bonds broken relative to the total number of bonds originally created. In each time step, some new bonds broke, and number of bonds broken in damage ratio must be a sum of broken bonds in all time step before current time step. Thus, I still need an answer. Huihuang Xia -- 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 #611087]: I got a wrong damage ratio
New question #611087 on Yade:
https://answers.launchpad.net/yade/+question/611087
Hi,
https://answers.launchpad.net/yade/+question/432617, in this question, damage
ratio was defined and implemented in YADE.
According to this definition, I built a impact crushing simulation using damage
ratio in YADE, but my code got a wrong result, as damage ratio in this
simulation was greater than 1.
Here is my code:
#!/usr/bin/python
# -*- coding: utf-8 -*-
#Created by Huihuang Xia from Huaqiao University,Xiamen,P.R.C
#My E-Mail:huihuang...@sina.com
###
# IMPORT MODULES
###
from yade import pack
from yade import plot
from yade import qt
###
# DEFINE MATERIALS
###
rock=CohFrictMat(young=5.98e7,poisson=0.3,alphaKr=3000,alphaKtw=3000,density=2678,frictionAngle=0.5,isCohesive=True,normalCohesion=7.9e6,shearCohesion=7.9e6,momentRotationLaw=True)
O.materials.append(rock)
steel=CohFrictMat(young=3.06e11,poisson=0.29,density=7861,frictionAngle=0.545,normalCohesion=0,shearCohesion=0)
O.materials.append(steel)
#
# CREATE SAMPLE & RIGID_WALL
#
wall=O.bodies.append(geom.facetBox(center=(0,0,0),extents=(0.015,0.015,0.0005),color=(1,1,0),material=steel,fixed=True))
pred=pack.inSphere(center=(0,0,0.006),radius=0.005)
assembly=pack.randomDensePack(pred,radius=0.0002,rRelFuzz=0.5,spheresInCell=2500,material=rock)
O.bodies.append(assembly)
for b in assembly:
b.state.vel=(0,0,-20)
###
# ENGINES
###
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=1.0),Bo1_Facet_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom6D(interactionDetectionFactor=1.0),Ig2_Facet_Sphere_ScGeom6D()],
[Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True)],
[Law2_ScGeom6D_CohFrictPhys_CohesionMoment()]
),
VTKRecorder(fileName='post/impact-',recorders=['all'],iterPeriod=250),
GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=5,timestepSafetyCoefficient=0.8,defaultDt=PWaveTimeStep()),
NewtonIntegrator(gravity=(0,0,-9.81)),
PyRunner(command='damageRatio()',realPeriod=5),
PyRunner(command='addPlotData()',realPeriod=5),
]
O.trackEnergy=True
O.step()
###
# DEFINE FUNCATIONS
###
#define a function to save damage ratio
global sumCohBonds
sumCohBonds=0
for b in O.interactions:
if b.phys.cohesionBroken==True:
continue
sumCohBonds+=1 # calculate total intact bonds
numBroCohBonds=0
def damageRatio():
global numBroCohBonds
for br in O.interactions:
if br.phys.cohesionBroken==False:
continue
numBroCohBonds+=1 # calculate broken bonds
damageRatio=numBroCohBonds/float(sumCohBonds)
print sumCohBonds,numBroCohBonds
plot.saveDataTxt('data/damageratio.txt')
def addPlotData():
damageRatio=numBroCohBonds/float(sumCohBonds)
plot.addData(t=O.time,dr=damageRatio)
plot.plots={'t':('dr')}
plot.plot()
#set an optimal timestep
O.dt=utils.PWaveTimeStep()
O.usesTimeStepper=True
#3D view and controller
qt.View()
qt.Controller()
and here is the damage ratio vs time:
# drt
0.0 9.49679622677e-09
0.0044209618287 1.60495856089e-06
0.0155272805693 4.9193404254e-06
0.0266335993099 8.39516779415e-06
0.0362303213284 1.10352770982e-05
0.047684063 1.44826140458e-05
0.0585507871469 1.73886336283e-05
0.0679318524908 2.02186788283e-05
0.0788225145568 2.3419099093e-05
0.091654086694 2.71228495603e-05
0.104485658831 3.03992441997e-05
0.117856372655 3.34097285591e-05
0.130687944792 3.48247511657e-05
0.144058658615 3.669561998e-05
0.158291999137 3.75218412351e-05
0.172633167997 3.86994439409e-05
0.190748328661 4.02569184773e-05
0.210696571059 4.18618769955e-05
0.236144058659 4.52902203506e-05
0.278089281863 4.98581792485e-05
0.339120120768 5.47350982335e-05
0.42009920207 5.83200136999e-05
0.520379555747 6.13408059573e-05
0.640931636834 6.41705226598e-05
0.776148371792 6.70709729927e-05
0.946409316368 7.14288492846e-05
1.15074401553 7.49846615977e-05
1.39076989433 7.87721800155e-05
1.673927108048.25338810098e-05
1.99611817986 8.60126620602e-05
2.32769031702 8.71703956744e-05
2.67360362303 8.81834109127e-05
3.04539573 9.07206891265e-05
3.46236791029 9.43729049081e-05
3.92160879879 9.80223033349e-05
4.41718783696 0.000101386989899
4.94899719646 0.000104757207048
5.50323485012 0.000107345725132
6.08561569981 0.000110695198656
6.67565236144 0.000111623153105
7.272697865 0.000112554955106
7.87524261376 0.000113345795924
8.48781539789 0.000116016886332
9.10739702394 0.000118131382236
9.72881173172 0.000120064040351
10.3572352814 0.00012179727
10.9897563080.00012384999
Re: [Yade-users] [Question #517125]: How to set an correct Id of a particle assembly
Question #517125 on Yade changed: https://answers.launchpad.net/yade/+question/517125 Status: Answered => Solved Huihuang Xia 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
[Yade-users] [Question #517125]: How to set an correct Id of a particle assembly
New question #517125 on Yade: https://answers.launchpad.net/yade/+question/517125 Hi, This question bothers me for a long time, i.e. when I create a particle assembly, I do want to use body.state.vel and body.state.blockedDOFs to set a particular motion for the particle assembly. But, I failed for several times. I have tried this two methods, the first one is posted in https://answers.launchpad.net/yade/+question/404127, namely ball_id = O.bodies.append(...) ball = O.bodies[ball_id] but I got an error: Traceback (most recent call last): File "/home/huihuang/YADE/install/bin/yade", line 182, in runScript execfile(script,globals()) File "cube.py", line 8, in cube.state.blockDOFs='xyzXY' AttributeError: 'list' object has no attribute 'state' The second method is: cube=pack.randomDensePack(pred,radius=0.5,rRelFuzz=0.5,spheresInCell=1000) O.bodies.append(cube) cube.state.blockDOFs='xyzXY' cube.state.vel=(0,0,-0.1) and this one also got a same error, AttributeError: 'list' object has no attribute 'state'. Can anyone help me solving this problem? Thanks in advance, Huihuang Xia -- 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 #473506]: How to display force chain with varying width in ParaView
Question #473506 on Yade changed: https://answers.launchpad.net/yade/+question/473506 Status: Open => Solved Huihuang Xia 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 #473506]: How to display force chain with varying width in ParaView
Question #473506 on Yade changed: https://answers.launchpad.net/yade/+question/473506 Huihuang Xia gave more information on the question: Besides, my YADE is a development version (2017.01a). -- 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 #473506]: How to display force chain with varying width in ParaView
Question #473506 on Yade changed: https://answers.launchpad.net/yade/+question/473506 Status: Answered => Open Huihuang Xia is still having a problem: Hello, Jan, thanks for your answer, but [1] only got an equivalent width of all the force chains as I described above. Thus, I still need an answer. Many thanks, Huihuang Xia -- 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 #473506]: How to display force chain with varying width in ParaView
New question #473506 on Yade: https://answers.launchpad.net/yade/+question/473506 Hi, In post-processing, I plan to display force chain in ParaView. Thus, in my intrs files, CellDatatoPointData and Tube were used, but I can not adjust the width of the force chain corresponding to the magnitude of the forces, all the tubes were the same radius. Figure 13b in a publication titled 'YADE-OPEN DEM: an opensource software using a discrete element method to simulate granular material' displayed a varying width of the force chain corresponding to the magnitude of the forces. Can you tell me how to display my results like that? Many thanks, Huihuang Xia -- 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 #464190]: How to choose an appropriate material property
Question #464190 on Yade changed: https://answers.launchpad.net/yade/+question/464190 Status: Answered => Solved Huihuang Xia confirmed that the question is solved: Thanks Jérôme Duriez, 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 #464190]: How to choose an appropriate material property
New question #464190 on Yade: https://answers.launchpad.net/yade/+question/464190 Hi， I have some questions on how to choose a suitable material property, the details are as follows: I plan to model an impact crusher using YADE, rock particles are the raw material crushed in the crusher, and the crusher is made of steel. Thus, I try to use CohFrictMat and bond many sub-particles to model the rock particle, and choose FrictMat to define steel. But in InteractionLoop, I cannot find a suitable Ip2 to set the interactions between rock particles and steel. In a word, can you recommend two appropriate material properties for me to define rock and steel? Many thanks, Huihuang Xia -- 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 #404127]: An impact simulation got an error
Question #404127 on Yade changed: https://answers.launchpad.net/yade/+question/404127 Status: Answered => Solved Huihuang Xia 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
[Yade-users] [Question #404127]: An impact simulation got an error
New question #404127 on Yade: https://answers.launchpad.net/yade/+question/404127 Hi, I simulated a steel ball impacted on a particle assembly at a high speed, but my code got an error, that is the steel ball not move. My codes is as follows, please help me to solve this problem. Many thanks! from yade import pack # define materials for particles and mill # rock=O.materials.append(CohFrictMat(young=5.98e7,poisson=0.3,density=2678,frictionAngle=0.5,normalCohesion=1.2e10,shearCohesion=1.2e9,momentRotationLaw=True)) steel=O.materials.append(FrictMat(young=3.06e11,poisson=0.29,density=7861,frictionAngle=0.545)) # create partciles assembly # pred=pack.inAlignedBox((-0.1,-0.1,-0.1),(0.1,0.1,0.1)) spheres=pack.randomDensePack(pred,radius=0.005,rRelFuzz=0.3,spheresInCell=1500,material=rock,color=(1,0,1)) O.bodies.append(spheres) # create a rigid ball # ball=O.bodies.append(sphere(center=(0,0,0.5),radius=0.1,color=(0,1,1),material=steel)) ball.blockedDOFs='xyXYZ' ball.state.vel=Vector3(0,0,-1) # set velocity along axis-Z O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom6D(interactionDetectionFactor=1.1)], [Ip2_FrictMat_FrictMat_FrictPhys(),Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True)], [Law2_ScGeom_FrictPhys_CundallStrack(),Law2_ScGeom6D_CohFrictPhys_CohesionMoment()] ), NewtonIntegrator(damping=0.3,gravity=(0,0,-9.81)) ] from yade import qt qt.View() qt.Controller() #qt.Renderer().shape=0 #qt.Renderer().intrPhys=1 O.dt=0.001*PWaveTimeStep() -- 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 #394789]: Using radius expansion method to bulid a DEM model
Question #394789 on Yade changed: https://answers.launchpad.net/yade/+question/394789 Huihuang Xia confirmed that the question is solved: Thanks Christian Jakob, 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 #403669]: Uniaxial compression simulation disappeared suddenly in Yade
Question #403669 on Yade changed: https://answers.launchpad.net/yade/+question/403669 Status: Answered => Solved Huihuang Xia confirmed that the question is solved: Thanks Jérôme Duriez, 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 #403669]: Uniaxial compression simulation disappeared suddenly in Yade
New question #403669 on Yade: https://answers.launchpad.net/yade/+question/403669 Hi, I am a beginning learner of Yade, I performed a uniaxial compression in Yade, but this simulation disappeared suddenly. Can you help me to solve this peoblem? Thanks in advance. My codes is as followings: from yade import pack O.periodic=True O.cell.refSize=(0.05,0.05,0.15) #specify materials for wall & particles spheres=O.materials.append(CohFrictMat(young=5.98e7,density=2.7e3,poisson=0.3,frictionAngle=0.6)) wallMat=O.materials.append(CohFrictMat(young=2.06e17,density=7.9e3,poisson=0.29,frictionAngle=0.8)) O.bodies.append( pack.regularHexa(pack.inAlignedBox((0,0,0),(0.05,0.05,0.15)),radius=0.001,gap=0,color=(0,1,0),material=spheres) ) #creat upper & lower platens from yade import utils O.bodies.append(utils.geom.facetBox((0.025,0.025,0.15),(0.055,0.055,0),wallMask=20,material=wallMat,color=(1,0,1),dynamic=1)) O.bodies.append(utils.geom.facetBox((0.025,0.025,0),(0.055,0.055,0),wallMask=21,material=wallMat,color=(0,0,1),dynamic=0,fixed=1)) #Engines O.engines=[ ForceResetter(), InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb()]), InteractionLoop( [Ig2_Sphere_Sphere_ScGeom6D(),Ig2_Box_Sphere_ScGeom6D()], [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True,setCohesionOnNewContacts=False),Ip2_FrictMat_FrictMat_FrictPhys()], [Law2_ScGeom6D_CohFrictPhys_CohesionMoment()] ), TranslationEngine(translationAxis=(0,0,-1),velocity=(0,0,0.005),id=20), NewtonIntegrator(damping=0.2,gravity=[0,0,-9.81]), VTKRecorder(fileName='post/comp-',recorders=['all'],iterPeriod=100) ] #set a suitable timestep for this simulation O.dt=0.01*PWaveTimeStep() #3D view & controller from yade import qt qt.View() qt.Controller() O.run() And here is the error: huihuang@huihuang-Presario-CQ42-Notebook-PC:~$ cd Desktop/test huihuang@huihuang-Presario-CQ42-Notebook-PC:~/Desktop/test$ yade compression_test.py Welcome to Yade 1.20.0 TCP python prompt on localhost:9002, auth cookie `ykasde' XMLRPC info provider on http://localhost:21002 XmbTextListToTextProperty result code -2 XmbTextListToTextProperty result code -2 XmbTextListToTextProperty result code -2 Running script compression_test.py WARN /build/yade-KKgSmd/yade-1.20.0/core/Cell.hpp:115 setRefSize: Setting Cell.refSize is deprecated, use Cell.setBox(...) instead. Traceback (most recent call last): File "/usr/bin/yade", line 182, in runScript execfile(script,globals()) File "compression_test.py", line 31, in TranslationEngine(translationAxis=(0,0,-1),velocity=(0,0,0.005),id=20), TypeError: No registered converter was able to produce a C++ rvalue of type double from this Python object of type tuple [[ ^L clears screen, ^U kills line. F12 controller, F11 3d view (use h-key for showing help), F10 both, F9 generator, F8 plot. ]] [0;34mYade [[1;34m1[0;34m]: [0mXmbTextListToTextProperty result code -2 XmbTextListToTextProperty result code -2 XmbTextListToTextProperty result code -2 python: /build/yade-KKgSmd/yade-1.20.0/pkg/common/InsertionSortCollider.cpp:493: bool InsertionSortCollider::spatialOverlapPeri(Body::id_t, Body::id_t, Scene*, Vector3i&) const: Assertion `maxima[3*id2+axis]-minima[3*id2+axis]<.99*dim' failed. Aborted (core dumped) -- 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 #403219]: randomDensePack got a strange error
Question #403219 on Yade changed: https://answers.launchpad.net/yade/+question/403219 Huihuang Xia posted a new comment: Dear Jérôme Duriez, Thanks for your help, I will read "How to ask questions as a user of Yade" carefully. Regards. -- 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 #403219]: randomDensePack got a strange error
Question #403219 on Yade changed: https://answers.launchpad.net/yade/+question/403219 Status: Answered => Solved Huihuang Xia confirmed that the question is solved: Thanks Jérôme Duriez, 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 #403219]: randomDensePack got a strange error
Question #403219 on Yade changed:
https://answers.launchpad.net/yade/+question/403219
Description changed to:
Hi,
When I conduct simulation using randomDensePack in YADE(vision 1.20.0), a
strange error occur as follows:
Welcome to Yade 1.20.0
TCP python prompt on localhost:9000, auth cookie `caseus'
XMLRPC info provider on http://localhost:21000
[[ ^L clears screen, ^U kills line. F12 controller, F11 3d view (use h-key for
showing help), F10 both, F9 generator, F8 plot. ]]
Yade [1]: from yade import pack
Yade [2]: ## construct the predicate first
Yade [3]:
pred=pack.inHyperboloid(centerBottom=(0,0,-.1),centerTop=(0,0,.1),radius=.05,skirt=.03)
Yade [4]: ## alternatively:
pack.inHyperboloid((0,0,-.1),(0,0,.1),.05,.03)
Yade [5]: ## pack the predicate with spheres (will be explained later)
Yade [6]: spheres=pack.randomDensePack(pred,spheresInCell=2000,radius=3.5e-3)
No suitable packing in database found, running PERIODIC compression
/usr/lib/x86_64-linux-gnu/yade/py/yade/pack.py:296: FutureWarning: The default
behavior will change; specify returnSpherePack=True for the new behavior, and
False to get rid of this warning (your code will break in the future, however).
The returned SpherePack object can be added to the simulation using
SpherePack.toSimulation()
warnings.warn('The default behavior will change; specify
returnSpherePack=True for the new behavior, and False to get rid of this
warning (your code will break in the future, however). The returned SpherePack
object can be added to the simulation using
SpherePack.toSimulation()',category=FutureWarning)
Yade [7]: ## add spheres to simulation
Yade [8]: O.bodies.append(spheres)
-> [8]:
[0,
1,
2,
3,
4,
5,
Finally, Yade got a same results as in the tutorial, but when I using
randomDensePack to conduct other simulation, it got only this error (No
suitable packing in database found, running PERIODIC compression...) without
any particles displayed in the 3D view.
Can you help me to solve this problem?
Many 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
[Yade-users] [Question #403219]: randomDensePack got a strange error
New question #403219 on Yade:
https://answers.launchpad.net/yade/+question/403219
Hi,
When I conduct simulation using randomDensePack in YADE(vision 1.20.0), a
strange error occur as follows:
Welcome to Yade 1.20.0
TCP python prompt on localhost:9000, auth cookie `caseus'
XMLRPC info provider on http://localhost:21000
[[ ^L clears screen, ^U kills line. F12 controller, F11 3d view (use h-key for
showing help), F10 both, F9 generator, F8 plot. ]]
Yade [1]: from yade import pack
Yade [2]: ## construct the predicate first
Yade [3]:
pred=pack.inHyperboloid(centerBottom=(0,0,-.1),centerTop=(0,0,.1),radius=.05,skirt=.03)
Yade [4]: ## alternatively: pack.inHyperboloid((0,0,-.1),(0,0,.1),.05,.03)
Yade [5]: ## pack the predicate with spheres (will be explained later)
Yade [6]: spheres=pack.randomDensePack(pred,spheresInCell=2000,radius=3.5e-3)
No suitable packing in database found, running PERIODIC compression
/usr/lib/x86_64-linux-gnu/yade/py/yade/pack.py:296: FutureWarning: The default
behavior will change; specify returnSpherePack=True for the new behavior, and
False to get rid of this warning (your code will break in the future, however).
The returned SpherePack object can be added to the simulation using
SpherePack.toSimulation()
warnings.warn('The default behavior will change; specify
returnSpherePack=True for the new behavior, and False to get rid of this
warning (your code will break in the future, however). The returned SpherePack
object can be added to the simulation using
SpherePack.toSimulation()',category=FutureWarning)
Yade [7]: ## add spheres to simulation
Yade [8]: O.bodies.append(spheres)
-> [8]:
[0,
1,
2,
3,
4,
5,
Finally, Yade got a same results as in the tutorial, but when I using
randomDensePack to conduct other simulation, it got only this error without any
particles displayed in the 3D view.
Can you help me to solve this problem?
Many 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
Re: [Yade-users] [Question #402952]: Using a gts surface as a shell and filling it using radius expansion method
Question #402952 on Yade changed: https://answers.launchpad.net/yade/+question/402952 Huihuang Xia 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 #402952]: Using a gts surface as a shell and filling it using radius expansion method
Question #402952 on Yade changed: https://answers.launchpad.net/yade/+question/402952 Status: Answered => Solved Huihuang Xia confirmed that the question is solved: Hello Jan, Thanks for your responses, you solved all of my problems and confusions. Best regards, Huihuang Xia -- 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 #402952]: Using a gts surface as a shell and filling it using radius expansion method
Question #402952 on Yade changed: https://answers.launchpad.net/yade/+question/402952 Status: Answered => Open Huihuang Xia is still having a problem: Dear Jan, Thanks for you patience and kindness, I plan to fill a gts surface like a rugby using many sub-particles in YADE. I have saw your examples in github, but the polyhedron was created by some parameters rather than importing a gts or stl surface from a CAD package. I wonder whether rugby like shape packed by sub-particles can be used as a new particle to perform another simulation in YADE?? Many thanks and regards. Huihuang Xia -- 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 #402952]: Using a gts surface as a shell and filling it using radius expansion method
Question #402952 on Yade changed: https://answers.launchpad.net/yade/+question/402952 Status: Needs information => Open Huihuang Xia gave more information on the question: Hi, The polyhedral particle in YADE is crushable?? If yes, this polyhedral particle is a good choice for my simulation, can you give me an example? Thanks, xia -- 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 #402952]: Using a gts surface as a shell and filling it using radius expansion method
New question #402952 on Yade: https://answers.launchpad.net/yade/+question/402952 Hi, I built an irregular gts surface in a CAD software, it looks like a sand or stone particle. I want to use dozens of sub-particle to fill this irregular shell using radius expansion method. After that, I plan to use many these irregular particle to conduct a compression test. I wonder whether YADE can perform this task? If yes, can you give me a similar example for me to copy and study? Many thanks and regards. Huihuang xia -- 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 #394789]: Using radius expansion method to bulid a DEM model
Question #394789 on Yade changed: https://answers.launchpad.net/yade/+question/394789 Status: Answered => Solved Huihuang Xia confirmed that the question is solved: ok, I will try it, many thanks! Huihuang Xia -- 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 #394789]: Using radius expansion method to bulid a DEM model
New question #394789 on Yade: https://answers.launchpad.net/yade/+question/394789 Dear developers of Yade, I am a new user of Yade, I want to build a DEM model using many sub-spheres to fill a given irregular shape with radius expansion method, such as a sand particle. But, I cannot find even an example using this radius expansion method in Yade's tutorial. Can you give me such an example to help me improving my programming skills? Many thanks and regards, Huihuang Xia -- 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
