uncertainty in the shear behavior of material
Hello All,
I was encountered to a problem in my dem modeling. I want to simulate the bulk behavior of material. Due to reach this goal, I have to choose a number of particles which could represent the bulk behavior as similar as experiment. I did two type of simulation.
in the first one, I put the particles in the box and the the box faces compress the paricles, after in the x and y direction, stress remain constant and in the z direction, we increase stress and plot the strain-stress in z direction. Due to initial conditions (randomness of initial position of particles) we will reach to different results in each realization. The code is as below:
from yade import utils, plot
from yade import pack, qt
from datetime import datetime
from yade import export
O.periodic=True
O.cell.
O.materials.
nums=['t']
mats=['Coke']
coke=(5e-3,3000)
nums=pack.
nums.makeCloud(
O.bodies.
sigmaIso=-1e5
O.engines=[
ForceReset
InsertionS
Interactio
),
PeriTriaxC
goal=
dynCell=
maxUnbalanced
doneHook=
),
NewtonInte
PyRunner(
]
O.dt=1e-6
O.saveTmp()
def compactionFinis
O.cell.
triax.
triax.stressMask=3
triax.
triax.
triax.
dataOld = plot.data
plot.data = {}
plot.plot()
def history():
plot.addData(
sxx=-
exx=-
print 'porosity:', porosity()
print 'iteration', O.iter
O.run(8000000,True)
plot.saveDataTx
plot.plots=
plot.plot()
def triaxFinished():
print 'Compaction finished'
O.pause()
In the second one, I made the sample first with the code below:
from yade import utils, plot
from yade import pack, qt
from datetime import datetime
from yade import export
O.periodic=True
O.cell.
O.materials.
nums=['t']
mats=['Coke']
coke=(5e-3,3000)
nums=pack.
nums.makeCloud(
O.bodies.
sigmaIso=-1e5
O.engines=[
ForceReset
InsertionS
Interactio
),
PeriTriaxC
goal=
dynCell=
maxUnbalanced
doneHook=
),
NewtonInte
PyRunner(
]
O.dt=1e-6
O.saveTmp()
def compactionFinis
print 'Finished'
O.pause()
plot.data = {}
O.run(1500000,True)
and after that, I compress the sample of "Pre-compact300
from yade import utils, plot
from yade import pack, qt
from datetime import datetime
from yade import export
O.load(
sigmaIso=-1e5
def compactionFinis
O.cell.
triax.
triax.stressMask=3
triax.
triax.
triax.
dataOld = plot.data
# plot.saveDataTx
plot.data = {}
plot.plot()
def history():
plot.addData(
sxx=-
exx=-
print 'porosity:', porosity()
print 'iteration', O.iter
O.run(5000000,True)
plot.saveDataTx
plot.plots=
plot.plot()
def triaxFinished():
print 'Compaction finished'
O.pause()
However, when I run the code again, I reach to a different results. But we know that the initial sample is the same.
Where does this uncertainty come from?Why the stress-strain curves are different for the same initial sample?
Thank you very much for your patient and help in advance.
Best Reagrds
Alireza
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Robert Caulk
- Solved:
- Last query:
- Last reply: