[Yade-users] New version of Yade, 2022. Release plan
Dear Yade users and developers, as always at the beginning of January we want to release a new Yade version. Release process takes some time, so please commit all your planned features till the end of the December 2021, so we can prepare tarball, test it on all supported architectures and upload it into the package archives. The version 2022.01 should go into the next Long-term-support Ubuntu Release, which is planned to be released in April 2022 and will be supported till 2027, and even with Extended Security Maintenance till 2032. Please plan your work accordingly. Thanks and best regards Anton ___ 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 #699362]: The uniaxial compression experiment fails to reach the critical state
Question #699362 on Yade changed: https://answers.launchpad.net/yade/+question/699362 Status: Needs information => Answered Robert Caulk proposed the following answer: Hello, Please, no external links [1]. >>Usually, the critical state can be reached at 25% of the axial strain. Why is this? What is "usually?" Does that mean you are comparing to a common "triaxial" lab test with flexible membrane? Keep in mind that the walls are stiff and frictional in your simulation here so you should expect different behavior. Cheers, Robert [1] https://www.yade-dem.org/wiki/Howtoask -- 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 #699362]: The uniaxial compression experiment fails to reach the critical state
Question #699362 on Yade changed: https://answers.launchpad.net/yade/+question/699362 Status: Open => Needs information Jan Stránský requested more information: Hello, > The uniaxial compression ... > After isotropic consolidation, the confining pressure remains unchanged in > the X and Y directions ... usually "uniaxial compression", uniaxial stress etc. refers to a scenario, where the lateral stress is zero. This I would call triaxial compression. > Usually, the critical state can be reached at 25% of the axial strain. Please be more specific. What is "usually"? In experiments? simulations? > Why is this? Because of the combination of packing, material model and loading :-) Cheers Jan -- 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 #699346]: How to simulate the material with big elastoplastic property such as PLA
Question #699346 on Yade changed: https://answers.launchpad.net/yade/+question/699346 Wang Benxin posted a new comment: Hi Robert, I think this would be of great help for the problem. And I'd like to try. Thanks and best regards, Benxin -- 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 #699345]: How to calibrate the magnitude of local stress field in 2D specimen
Question #699345 on Yade changed: https://answers.launchpad.net/yade/+question/699345 Wang Benxin posted a new comment: Hi Robert, In experiment the stress field was calculated base on the damage evolution constitutive equation. And it's magnitude is 1e6 Pa. The numerical test uses the Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM. And the local stress in numerical model is exported by VTK and showed in ParaView. -- 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 #699362]: The uniaxial compression experiment fails to reach the critical state
New question #699362 on Yade:
https://answers.launchpad.net/yade/+question/699362
After isotropic consolidation, the confining pressure remains unchanged in the
X and Y directions, and the axial compression in the Z-axis direction is 40%,
which has not yet reached the critical state. Usually, the critical state can
be reached at 25% of the axial strain. Why is this? The procedure is as
follows：
from __future__ import print_function
from yade import pack, qt, plot
from math import *
nRead=readParamsFromTable(
## model parameters
num_spheres=3,
targetPorosity= .32,
confiningPressure=-5,
## material parameters
compFricDegree=25,#contact friction during the confining phase
finalFricDegree=30,#contact friction during the deviatoric loading
young=2e8,
poisson=.2,
density=2650,
alphaKr=7.5,
alphaKtw=0,
competaRoll=.1,
finaletaRoll=.1,
etaTwist=0,
normalCohesion=0,
shearCohesion=0,
## fluid parameters
fluidDensity=1000,
dynamicViscosity=.001,
## control parameters
damp=0,
stabilityThreshold=.001,
## output specifications
filename='suffusion',
unknowOk=True
)
from yade.params.table import *
O.load('compactedStatesuffusion25.yade.gz')
##__third section, deviatoric
loading__
triax.width0=aabbDim()[0]
triax.height0=aabbDim()[1]
triax.depth0=aabbDim()[2]
# change control type: keep constant confinement in x,y, 40% compression in z
triax.stressMask=3
triax.goal3=-.1
triax.goal1=confiningPressure
triax.goal2=confiningPressure
# enable energy tracking in the code
O.trackEnergy=True
# define function to record history
def history():
sig_xx=-triax.stress(triax.wall_right_id)[0]
sig_yy=-triax.stress(triax.wall_top_id)[1]
sig_zz=-triax.stress(triax.wall_front_id)[2]
plot.addData(unbalanced=unbalancedForce(),i=O.iter,exx=-triax.strain[0],
eyy=-triax.strain[1], ezz=-triax.strain[2],
sxx=sig_xx,syy=sig_yy,szz=sig_zz,
ev=-triax.strain[0]-triax.strain[1]-triax.strain[2],
dsratio=(-sig_zz+sig_xx)/((-sig_xx-sig_yy-sig_zz)/3.0),
porosity=porosity(),Etot=O.energy.total(),**O.energy# save all
available energy data
)
O.engines=O.engines+[PyRunner(command='history()',iterPeriod=500)]
t0=O.time
while True:
O.run(1000,True)
print('running',-triax.strain[2])
unb=unbalancedForce()
if unbhttps://ibb.co/Yc7jbbQ
Obviously the volumetric strain is still changing and has not reached the
critical state
--
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 #699347]: How to use periodic boundaries for triaxial test simulation after erosion
Question #699347 on Yade changed: https://answers.launchpad.net/yade/+question/699347 Status: Needs information => Solved Zhicheng Gao confirmed that the question is solved: > This feature is also what I am researching. In other words, is it best to > delete the wall directly and impose periodic boundary conditions? I am studying the effect of erosion on granular materials. -- 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 #699345]: How to calibrate the magnitude of local stress field in 2D specimen
Question #699345 on Yade changed: https://answers.launchpad.net/yade/+question/699345 Robert Caulk requested more information: I guess this question depends entirely on how you are estimating the local stress field. Please indicate this information if you wish to receive accurate assistance. -- 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 #699346]: How to simulate the material with big elastoplastic property such as PLA
Question #699346 on Yade changed: https://answers.launchpad.net/yade/+question/699346 Status: Needs information => Answered Robert Caulk proposed the following answer: Hello, >>'PLA' is the polylactic acid. https://sv.wikipedia.org/wiki/Polylaktid I am unaware of any contact models in yade which were designed particularly for plastic materials. However, CohFrictMat does allow for perfectly plastic behavior with fragile=True [1]. Cheers, Robert [1]https://yade- dem.org/doc/yade.wrapper.html#yade.wrapper.CohFrictMat.fragile -- 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 #699345]: How to calibrate the magnitude of local stress field in 2D specimen
Question #699345 on Yade changed: https://answers.launchpad.net/yade/+question/699345 Wang Benxin posted a new comment: Hi Robert, >>Why do you presume that the micro-parameters should be of the same order of magnitude as the macro-parameters? I got the local stress field in experiment. Hence, I want to make the comparison and verification of local stress field between simulation and experiment. >>Have you considered increasing the interaction range to simulate grain interlocking better? [1] Yes, I had increased the interaction range as the paper saying. Thanks, Benxin -- 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 #699342]: generate two different radius particle in a box
Question #699342 on Yade changed: https://answers.launchpad.net/yade/+question/699342 Jan Stránský proposed the following answer: Hello, > I've just started using yade recently welcome :-) you can try something like (not tested): ### sp = pack.randomDensePack(...,returnSpherePack=True) # big particles sp.makeCloud(...) # small particles sp.toSimulation() ### Cheers Jan -- 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 #699342]: generate two different radius particle in a box
Question #699342 on Yade changed:
https://answers.launchpad.net/yade/+question/699342
Robert Caulk proposed the following answer:
Hello,
No, you cannot use randomDensePack to do this. But you can very easily
replicate randomDensePack action by increasing particle radii yourself
according to the stresses on the walls.
Add a pyrunner to your engines list:
PyRunner(command='addStress()',dead=1,iterPeriod=10,label='loadSample')
Then add the following function to your script:
## Set mn,mx. They should be the values you used to create the aabbWalls.
'walls' should be a list of the ids of your walls.
multiplier = 1.001
def addStress():
global multiplier
preStress = 5e3
maxMultiplier=1.001
growParticles(multiplier)
pressure = 0
for i in walls:
f = O.bodies[i]
a = mn[0]*mn[1]
force = O.forces.f(f.id,sync=True)
pressure = pressure + np.linalg.norm(force)/a
stabCoeff = abs((abs(pressure/len(walls))-preStress)/preStress)
facetStress = abs(pressure/len(walls))
print('unbalanced force:',unbalancedForce(),'avg top pressure
',facetStress ,'stability coeff ',stabCoeff)
multiplier = 1+(preStress-facetStress)/preStress * (maxMultiplier-1.)
if multiplier <= 0.95: multiplier=0.95
print('multiplier adjusted up to ', multiplier)
if stabCoeff<0.001:
print('packing finished, export it to file ;)')
loadSample.dead = 1
O.pause()
I hope it is obvious that this should not be active during any scientific
simulation, this is just to achieve a dense packing. Best practice would be to
export this packing to file using export.textExt() and then import that packing
into whatever you need to do. I leave the export/import up to you as homework
;).
Cheers,
Robert
--
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 #699342]: generate two different radius particle in a box
Question #699342 on Yade changed: https://answers.launchpad.net/yade/+question/699342 Status: Open => Answered Karol Brzezinski proposed the following answer: Hi, I don't think that you can do this with randomDensePack (I am not sure though). I think that the simplest approach would be to: - generate packing with makeCloud() inside a bigger cube, - generate six walls around the packing and compress the packing until some condition is reached (stress or porosity, etc.) - export your packing to an external file. If you need a more specific answer, you should share more details [1]. Especially, after you try to implement some solution by yourself. Cheers, Karol [1] https://www.yade-dem.org/wiki/Howtoask -- 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 #699346]: How to simulate the material with big elastoplastic property such as PLA
Question #699346 on Yade changed: https://answers.launchpad.net/yade/+question/699346 Wang Benxin posted a new comment: Hi Robert, Thank you for your interest in my question. >>I don't quite understand. What is "stable after the yield stage". Do you mean the material should be perfectly plastic? I guess yes, because you say "no strain softening." Yes, you are right. The material is perfectly plastic. Its property is like the Fig.9 in the literature https://link.springer.com/content/pdf/10.1007/s10409-015-0524-4.pdf >>What is "big plastic characteristic?" Can you describe that a bit more technically? "big", "small" what does it mean? are those relative values to some other behavior? "big", "small" are use to describe the behavior of the material like 'perfectly' you said. >>What is "PLA?" 'PLA' is the polylactic acid. https://sv.wikipedia.org/wiki/Polylaktid Cheers, Benxin -- 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 #699345]: How to calibrate the magnitude of local stress field in 2D specimen
Question #699345 on Yade changed: https://answers.launchpad.net/yade/+question/699345 Status: Open => Answered Robert Caulk proposed the following answer: Hello, >>However the bonding parameters such as cohesion and tensile strength are in the magnitude about 1e10 Pa which is too high if I want to get a macro stress-strain curve and UCS with the unit or magnitude MPa like the experimental tests. Why do you presume that the micro-parameters should be of the same order of magnitude as the macro-parameters? Have you considered increasing the interaction range to simulate grain interlocking better? [1] Cheers, Robert [1]https://www.sciencedirect.com/science/article/abs/pii/S0022509612002268 -- 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 #699346]: How to simulate the material with big elastoplastic property such as PLA
Question #699346 on Yade changed: https://answers.launchpad.net/yade/+question/699346 Status: Open => Needs information Robert Caulk requested more information: Hello, >>stable after the yield stage showing that it has a big plastic characteristic and there is not a strain softening stage. The material is just like PLA I don't quite understand. What is "stable after the yield stage". Do you mean the material should be perfectly plastic? I guess yes, because you say "no strain softening." What is "big plastic characteristic?" Can you describe that a bit more technically? "big", "small" what does it mean? are those relative values to some other behavior? What is "PLA?" Cheers, Robert -- 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 #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Status: Expired => Needs information Robert Caulk changed the question status: Avoiding launchpad janitor -- 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 #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Robert Caulk posted a new comment: Hello, Please confirm that updating your yade installation has solved your problem. If not, please indicate the exact changes you have made compared to the example script, and which change in particular is the one that results in your issue. Thank you, Robert -- 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 #699072]: About the examples/ThermalEngine
Question #699072 on Yade changed: https://answers.launchpad.net/yade/+question/699072 Status: Open => Expired Launchpad Janitor expired the question: This question was expired because it remained in the 'Open' state without activity for the last 15 days. -- 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
