[Yade-users] [Question #688491]: Error when manually saving a state in yade.gz
New question #688491 on Yade:
https://answers.launchpad.net/yade/+question/688491
Hi there,
I am using Yade 1.20.0 on Ubuntu 16.04 LTS (Xenial).
I used O.pause() followed by O.save('filename') in my script to save a
simulation state in yade.gz format, which works well and saved with no errors.
However, when I use the GUI pause and save tabs to manually pause and save the
simulation, I obtain the following error:
RuntimeError Traceback (most recent call last)
/usr/lib/x86_64-linux-gnu/yade/py/yade/__init__.py in saveSlot(self)
152 f=str(f)
153 if not f: return # cancelled
--> 154 O.save(f)
155 def reloadSlot(self):
156 self.deactivateControls()
RuntimeError: output stream error
Appreciate your feedback on how to resolve this. 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 #688400]: Energy of the body is very high, defying gravity and normal stress - cpm material
Question #688400 on Yade changed: https://answers.launchpad.net/yade/+question/688400 Status: Needs information => Open Akm gave more information on the question: I think that might have been the problem. I had given the tensile strength as too high and that has created a lot of problems for me. How do I know which version of Yade I am using? I had downloaded the package from https://github.com/yade/trunk. So it shows the version as Yade unknown. Should I be downloading it from somewhere else? Many thanks, Arun -- 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 #688489]: applying force boundary condition
New question #688489 on Yade: https://answers.launchpad.net/yade/+question/688489 Hi, To apply a velocity, I did the following: O.bodies[63].state.blockedDOFs = "xyzXYZ" O.bodies[63].state.vel = (0,-0.000100,0) O.bodies[63].state.angVel = (0,0,0) I want to apply a force to element 63 as a function of time. For example, 1 N/s in the x-direction. Thank you for your 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
[Yade-users] [Question #688488]: more than 255 arguments
New question #688488 on Yade: https://answers.launchpad.net/yade/+question/688488 Hi, I have a model where I define a function in which there are 300 variables that are save (mainly the poistion of some elements). I get an error saying that more than 255 arguments. i tried to make two different functions but that did not work as well. Any idea of how I can record positions of 300 elements? Thanks! # FUNCTIONAL COMPONENTS def addPlotData1(): plot.addData(t=O.time,Displacement = O.bodies[2450].state.pos[1] - 50, Force = -1*(O.forces.f(2450)[1]+ O.forces.f(2451)[1] + O.forces.f(2452)[1] + O.forces.f(2453)[1] + O.forces.f(2454)[1] + O.forces.f(2455)[1] + O.forces.f(2456)[1] + O.forces.f(2457)[1] + O.forces.f(2458)[1] + O.forces.f(2459)[1] + O.forces.f(2460)[1] + O.forces.f(2461)[1] + O.forces.f(2462)[1] + O.forces.f(2463)[1] + O.forces.f(2464)[1] + O.forces.f(2465)[1] + O.forces.f(2466)[1] + O.forces.f(2467)[1] + O.forces.f(2468)[1] + O.forces.f(2469)[1] + O.forces.f(2470)[1] + O.forces.f(2471)[1] + O.forces.f(2472)[1] + O.forces.f(2473)[1] + O.forces.f(2474)[1] + O.forces.f(2475)[1] + O.forces.f(2476)[1] + O.forces.f(2477)[1] + O.forces.f(2478)[1] + O.forces.f(2479)[1] + O.forces.f(2480)[1] + O.forces.f(2481)[1] + O.forces.f(2482)[1] + O.forces.f(2483)[1] + O.forces.f(2484)[1] + O.forces.f(2485)[1] + O.forces.f(2486)[1] + O.forces.f(2487)[1] + O.forces.f(2488)[1] + O.forces.f(2489)[1] + O.forces.f(2490)[1] + O.forces.f(2491)[1] + O.forces.f(2492)[1] + O.forces.f(2493)[1] + O.forces.f(2494)[1] + O.forces.f(2495)[1] + O.forces.f(2496)[1] + O.forces.f(2497)[1] + O.forces.f(2498)[1] + O.forces.f(2499)[1] ) ,DX1 = O.bodies[1250].state.pos[0], DY1 = O.bodies[1250].state.pos[1], DX2 = O.bodies[1251].state.pos[0], DY2 = O.bodies[1251].state.pos[1], DX3 = O.bodies[1252].state.pos[0], DY3 = O.bodies[1252].state.pos[1], DX4 = O.bodies[1253].state.pos[0], DY4 = O.bodies[1253].state.pos[1], DX5 = O.bodies[1254].state.pos[0], DY5 = O.bodies[1254].state.pos[1], DX6 = O.bodies[1255].state.pos[0], DY6 = O.bodies[1255].state.pos[1], DX7 = O.bodies[1256].state.pos[0], DY7 = O.bodies[1256].state.pos[1], DX8 = O.bodies[1257].state.pos[0], DY8 = O.bodies[1257].state.pos[1], DX9 = O.bodies[1258].state.pos[0], DY9 = O.bodies[1258].state.pos[1], DX10 = O.bodies[1259].state.pos[0], DY10 = O.bodies[1259].state.pos[1], DX11 = O.bodies[1260].state.pos[0], DY11 = O.bodies[1260].state.pos[1], DX12 = O.bodies[1261].state.pos[0], DY12 = O.bodies[1261].state.pos[1], DX13 = O.bodies[1262].state.pos[0], DY13 = O.bodies[1262].state.pos[1], DX14 = O.bodies[1263].state.pos[0], DY14 = O.bodies[1263].state.pos[1], DX15 = O.bodies[1264].state.pos[0], DY15 = O.bodies[1264].state.pos[1], DX16 = O.bodies[1265].state.pos[0], DY16 = O.bodies[1265].state.pos[1], DX17 = O.bodies[1266].state.pos[0], DY17 = O.bodies[1266].state.pos[1], DX18 = O.bodies[1267].state.pos[0], DY18 = O.bodies[1267].state.pos[1], DX19 = O.bodies[1268].state.pos[0], DY19 = O.bodies[1268].state.pos[1], DX20 = O.bodies[1269].state.pos[0], DY20 = O.bodies[1269].state.pos[1], DX21 = O.bodies[1270].state.pos[0], DY21 = O.bodies[1270].state.pos[1], DX22 = O.bodies[1271].state.pos[0], DY22 = O.bodies[1271].state.pos[1], DX23 = O.bodies[1272].state.pos[0], DY23 = O.bodies[1272].state.pos[1], DX24 = O.bodies[1273].state.pos[0], DY24 = O.bodies[1273].state.pos[1], DX25 = O.bodies[1274].state.pos[0], DY25 = O.bodies[1274].state.pos[1], DX26 = O.bodies[1275].state.pos[0], DY26 = O.bodies[1275].state.pos[1], DX27 = O.bodies[1276].state.pos[0], DY27 = O.bodies[1276].state.pos[1], DX28 = O.bodies[1277].state.pos[0], DY28 = O.bodies[1277].state.pos[1], DX29 = O.bodies[1278].state.pos[0], DY29 = O.bodies[1278].state.pos[1], DX30 = O.bodies[1279].state.pos[0], DY30 = O.bodies[1279].state.pos[1], DX31 = O.bodies[1280].state.pos[0], DY31 = O.bodies[1280].state.pos[1], DX32 = O.bodies[1281].state.pos[0], DY32 = O.bodies[1281].state.pos[1], DX33 = O.bodies[1282].state.pos[0], DY33 = O.bodies[1282].state.pos[1], DX34 = O.bodies[1283].state.pos[0], DY34 = O.bodies[1283].state.pos[1], DX35 = O.bodies[1284].state.pos[0], DY35 = O.bodies[1284].state.pos[1], DX36 = O.bodies[1285].state.pos[0], DY36 = O.bodies[1285].state.pos[1], DX37 = O.bodies[1286].state.pos[0], DY37 = O.bodies[1286].state.pos[1], DX38 = O.bodies[1287].state.pos[0], DY38 = O.bodies[1287].state.pos[1], DX39 = O.bodies[1288].state.pos[0], DY39 = O.bodies[1288].state.pos[1], DX40 = O.bodies[1289].state.pos[0], DY40 = O.bodies[1289].state.pos[1], DX41 = O.bodies[1290].state.pos[0], DY41 = O.bodies[1290].state.pos[1], DX42 = O.bodies[1291].state.pos[0], DY42 = O.bodies[1291].state.pos[1], DX43 = O.bodies[1292].state.pos[0], DY43 = O.bodies[1292].state.pos[1], DX44 = O.bodies[1293].state.pos[0], DY44 = O.bodies[1293].state.pos[1], DX45 = O.bodies[1294].state.pos[0], DY45 = O.bodies[1294].state.pos[1], DX46 = O.bodies[1295].state.pos[0], DY46 = O.bodies[1295].state.pos[
Re: [Yade-users] [Question #688485]: Cycling AddToClump example fails
Question #688485 on Yade changed:
https://answers.launchpad.net/yade/+question/688485
Status: Needs information => Open
Zoheir Khademian gave more information on the question:
Hi Jan,
Thanks for the quick reply.
My Yade version is 20191220-2907
and Ubuntu version is 18.04.3
and Python version is 2.7.17
Here is the MWE, which is only a "O.run()" different from the trunk example.
from __future__ import print_function
from builtins import range
from yade import pack,export,qt
#define material for all bodies:
#define engines:
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()]
),
NewtonIntegrator(damping=0.7,gravity=[0,0,-10])
]
#create a box:
id_box = O.bodies.append(box((0,0,0),(2,2,.1),fixed=True))
show how to use appendClumped():
#create 2 clumps:
clump1=O.bodies.appendClumped([\
sphere([0,0,1], radius=0.5),\
sphere([0.2,0,1], radius=0.5)\
])
clump2=O.bodies.appendClumped([\
sphere([3,1,2], radius=0.5),\
sphere([3.2,1,2], radius=0.5)\
])
#get clump ids:
id_clump1 = clump1[0]
id_clump2 = clump2[0]
#definition for getting informations from all clumps:
def getClumpInfo():
for b in O.bodies:
if b.isClump:
print('Clump ',b.id,' has following members:')
keys = list(b.shape.members.keys())
for ii in range(0,len(keys)):
print('- Body ',keys[ii])
print('inertia:',b.state.inertia)
print('mass:',b.state.mass,'\n')
show how to use addToClump():
#create a new sphere:
id_new=O.bodies.append(sphere([0,0.2,1], radius=0.5))
print('\nSTATE before adding sphere to clump ')
getClumpInfo()
#add a sphere to the clump:
O.bodies.addToClump([id_new],id_clump1)
print('\nSTATE after adding sphere to clump ')
getClumpInfo()
#add a clump to a clump:
O.bodies.addToClump([id_clump2],id_clump1)
print('\nSTATE after adding the second clump to clump ')
getClumpInfo()
#try to add clump member to a clump (should give error message):
#O.bodies.addToClump(1,id_clump1)
#try to add clump to the same clump (should give error message):
#O.bodies.addToClump(id_clump1,id_clump1)
O.dt=1e-6
print('\nPress Play button ... ')
renderer = qt.Renderer()
qt.View()
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 #685707]: Periodic flow and permeability
Question #685707 on Yade changed:
https://answers.launchpad.net/yade/+question/685707
Status: Answered => Open
Jibril Coulibaly is still having a problem:
Hi,
Thank you very much rcaulk and bruno-chareyre for your answers. Using
flow.imposePressure() and flow.averageVelocity() provide sensible
results.
I have a couple remaining questions:
- (1) : If I try to run the script in a loop to obtain permeability in
all 3 directions, I get sensible results on the first run, then I get
zero or nan on subsequent runs. I did not get these issues when using
walls (see script below)
- (2) : Running a similar script, but reading the particles coordinates from an
external file instead, I obtain the error message: "GS did not converge in 20k
iterations (maybe because the reference pressure is zero?)".
I tried to figure it out from the code but I am not sure what this message is
warning me about.
I appreciate your help,
Jibril B. Coulibaly
###
from yade import pack,Vector3
O.periodic=True
O.cell.hSize=Matrix3(0.1,0,0, 0,0.1,0, 0,0,0.1)
pos_center = Vector3(0.05,0.05,0.05)
sp=pack.SpherePack()
radius=5e-3
num=sp.makeCloud((0,0,0),(.1,.1,.1),radius,.2,1000,periodic=True)
O.bodies.append([sphere(s[0],s[1]) for s in sp])
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb()],verletDist=.05*radius),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()]
),
GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.8),
NewtonIntegrator(damping=0.),
PeriodicFlowEngine(dead=0,label="flow")
]
O.dt=0.1e-8
O.dynDt=False
flow.defTolerance=0.3
flow.meshUpdateInterval=300
flow.useSolver=0 # Solver >0 not supported for PBC
flow.permeabilityFactor=1
flow.viscosity=1
flow.imposePressure(pos_center,1.0)
for i in [0,1,2]:
if i==0:
gradPress = Vector3(1,0,0)
elif i==1:
gradPress = Vector3(0,1,0)
elif i==2:
gradPress = Vector3(0,0,1)
flow.updateTriangulation=True
flow.gradP=gradPress
flow.updateTriangulation=True
O.run(1,1)
aveVel = flow.averageVelocity()
print("average velocity=(",aveVel[0],aveVel[1],aveVel[2],")")
###
--
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 #688487]: constant head permeability code
New question #688487 on Yade: https://answers.launchpad.net/yade/+question/688487 Hello, Is it possible to get the code that was used to produce this paper? Tong, A. T., Catalano, E., & Chareyre, B. (2012). Pore-scale flow simulations: model predictions compared with experiments on bi-dispersed granular assemblies. Oil & Gas Science and Technology–Revue d’IFP Energies nouvelles, 67(5), 743-752. Thank you, Othman -- 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 #688485]: Cycling AddToClump example fails
Question #688485 on Yade changed: https://answers.launchpad.net/yade/+question/688485 Status: Open => Needs information Jan Stránský requested more information: Hello Zoheir, Your MWE and info about Yade and OS version would be appreciated [1] cheers Jan [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
[Yade-users] [Question #688485]: Cycling AddToClump example fails
New question #688485 on Yade: https://answers.launchpad.net/yade/+question/688485 Hello all, I was able to successfully run addToClump-example.py example at "https://github.com/yade/trunk/blob/master/examples/clumps/addToClump-example.py " Hopwver, when I try to settle the clumps under gravity, I get this error: python3.6: /builds/yade-dev/trunk/deb/yadedaily/core/Clump.cpp:61: void yade::Clump::addForceTorqueFromMembers(const yade::State*, yade::Scene*, yade::Vector3r&, yade::Vector3r&): Assertion `member->isClumpMember()' failed. Aborted (core dumped) Your help in this would be appreciated Zoheir -- 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 #688340]: Problem about running yade in parallel
Question #688340 on Yade changed: https://answers.launchpad.net/yade/+question/688340 Huang peilun posted a new comment: That's great. 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 #688340]: Problem about running yade in parallel
Question #688340 on Yade changed: https://answers.launchpad.net/yade/+question/688340 Janek Kozicki posted a new comment: And we will not have this bug ever again. Because a test for it was added to the pipeline ;) -- 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 #688340]: Problem about running yade in parallel
Question #688340 on Yade changed: https://answers.launchpad.net/yade/+question/688340 Status: Answered => Solved Huang peilun 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 #688400]: Energy of the body is very high, defying gravity and normal stress - cpm material
Question #688400 on Yade changed: https://answers.launchpad.net/yade/+question/688400 Status: Open => Needs information Jan Stránský requested more information: > It's 1.20.0. > I have only changed three lines in the source code. Comparing the provided file ConcretePM.cpp and the one from 1.20.0, there is much more than three lines changes.. > Were you able to figure out any discrepancy in my code that led to the > overall behaviour? > epsCrackOnset = 5e-1, > epsCrackOnset = 5e-2 Where these values come from? this epsCrackOnset would lead to non-realistically high strength, very roughly f=E*e=32e9*5e-1=16 GPa / 1,6 GPa. Then the material behaves as elastic and you will not get any damage.. The elastic sample cannot do anything else (damage is not possible and stiffness is too high to somehow deform dramatically) than to ''FLY'' up.. 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 #688400]: Energy of the body is very high, defying gravity and normal stress - cpm material
Question #688400 on Yade changed: https://answers.launchpad.net/yade/+question/688400 Jérôme Duriez posted a new comment: Yade 1.20.0 is extremely old (before June 2016), see https://launchpad.net/yade/+download -- 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 #688376]: Which kind of default distribution in makeCloud
Question #688376 on Yade changed: https://answers.launchpad.net/yade/+question/688376 Jérôme Duriez posted a new comment: Also note that referring to github as in your initial post is no longer relevant to understand current Yade. Yade development is now reflected on gitlab.com -- 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 #688340]: Problem about running yade in parallel
Question #688340 on Yade changed: https://answers.launchpad.net/yade/+question/688340 Status: Needs information => Answered Jérôme Duriez proposed the following answer: I think this issue is known and fixed / workarounded as of yadedaily 20200202-3270~7acb24d~bionic1: https://gitlab.com/yade-dev/trunk/issues/138 -- 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 #688400]: Energy of the body is very high, defying gravity and normal stress - cpm material
Question #688400 on Yade changed: https://answers.launchpad.net/yade/+question/688400 Status: Needs information => Open Akm gave more information on the question: Hi Jan, It's 1.20.0. Were you able to figure out any discrepancy in my code that led to the overall behaviour? The conditions are valid and experiments have been performed under those conditions. I was looking for an approximate behaviour like that. Please let me know if you found any issues out here. Many thanks, Arun -- 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 #688458]: The mismatach of the sample after using export.text()
Question #688458 on Yade changed: https://answers.launchpad.net/yade/+question/688458 Status: Answered => Solved jasperb 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 #688388]: On the interaction range of the JCFPm
Question #688388 on Yade changed: https://answers.launchpad.net/yade/+question/688388 Status: Answered => Solved jasperb 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 #688439]: effect of rolling stiffness on peak friction angle of sample
Question #688439 on Yade changed: https://answers.launchpad.net/yade/+question/688439 Luc Sibille proposed the following answer: > But my problem is about the amount of peak friction angle before considering > rolling stiffness (I mean without considering CohesionMomment) and after > that. when I use just FrictPhys_CundallStrack as the contact law, with the > Inter-Particle friction angle of 30 degree, I can't get peak friction angle > of experimental data that is about 40 degree. This is a classical result: with spherical particles and without rolling resistance, internal peak friction angle relatively high (typically 40 deg) cannot be reproduced. > Because of this problem, I tried to model a sample with more stiffness until > getting more peak friction angle. After using CohesionMoment and setting > alpaKr (for example alphaKr=0.5 or 1), I couldn't see any changes. Again, stiffnesses, as normal, tangential or rolling stiffnesses, if their are chosen sufficient high should not affect the failure properties as the peak friction angle. So it may be normal that peak friction angle is not affected by such stiffnesses. Then peak friction angle is more strongly affected by the contact friction angle and plastic rolling moment coefficient (or rolling friction coefficient) etaRoll (and not alphaKr). Besides, be sure that the rolling resistance is activated. It is not because you use Law2_ScGeom6D_CohFrictPhys_CohesionMoment that rolling resistance is necessarily activated: rolling resistance can be enable or disable depending on some arguments in Law2_ScGeom6D_CohFrictPhys_CohesionMoment and/or CohFrictMat Luc -- Luc Sibille Université Grenoble Alpes / IUT1 de Grenoble Laboratoire 3SR: Sols, Solides, Structures, Risques Tel lab.: +33 (0)4 76 82 63 48 Tel IUT: +33 (0)4 76 82 53 36 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 #688400]: Energy of the body is very high, defying gravity and normal stress - cpm material
Question #688400 on Yade changed: https://answers.launchpad.net/yade/+question/688400 Status: Open => Needs information Jan Stránský requested more information: >>what version of Yade you use? > It is 1.2.0 1.20.0? > I have only changed three lines in the source code. in this case, you could have post the three lines instead of complete file.. > We know that the physical nature of the experiment cannot be simulated here in the simulations. then you cannot expect physical nature of the results.. you can try some simple hand calculations, maybe you just imposed such conditions that real gravity and plate pressure simply does not have enough effect.. 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 #688458]: The mismatach of the sample after using export.text()
Question #688458 on Yade changed: https://answers.launchpad.net/yade/+question/688458 Status: Open => Answered Jan Stránský proposed the following answer: I guess you are using O.periodic=True. In this case, the particles' positions are virtually wrapped to the cell for both computational and visualization purposes, even they are "exploded" in reality. Before exporting, do ## for b in O.bodies: b.state.pos = O.cell.wrapPt(b.state.pos) export.text(...) ## 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 #688439]: effect of rolling stiffness on peak friction angle of sample
Question #688439 on Yade changed: https://answers.launchpad.net/yade/+question/688439 Elham Hosseinkhani posted a new comment: Hello Luc Sibille, and thanks for your comment. I've studied paper [1] before and I've observed that point you said (peak friction angle does not depend on alphaKr when alphaKr is sufficiently high). But my problem is about the amount of peak friction angle before considering rolling stiffness (I mean without considering CohesionMomment) and after that. when I use just FrictPhys_CundallStrack as the contact law, with the Inter-Particle friction angle of 30 degree, I can't get peak friction angle of experimental data that is about 40 degree. Because of this problem, I tried to model a sample with more stiffness until getting more peak friction angle. After using CohesionMoment and setting alpaKr (for example alphaKr=0.5 or 1), I couldn't see any changes. I should say that I know about using more Inter-Particle friction angle( higher than 30 Degree like 50, 60,...) and getting more peak friction angle in the macroscopic response. I've tried this dependency an observed exact results. But I want to see effects of rolling parameters like stiffness and strength in macroscopic parameters. Thanks Luc Sibille for your suggestion for more studying. I also will study paper [2]. >>[1] Discrete numerical modeling of loose soil with spherical particles and interparticle rolling friction RA Hosn, L Sibille, N Benahmed, B Chareyre >>[2] Quantitative prediction of discrete element models on complex loading paths Luc Sibille, Pascal Villard, Félix Darve, Rodaina Aboul Hosn -- 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
