Oedometric consolidation test
Hi friends :)
I am trying to conduct a simulation of a laboratory oedometric test used for testing the consolidation of a soil sample.
I went through the example in the Yade documentation already (https:/
Here is what I am intending to do:
(1) Create a cylinder and fill it with spherical particles
from yade import pack,plot
pred=
sp=
O.bodies.
(2) For the sampler (container): create a facet cylinder as follows (the cylinder is to be closed from all sides -( top, bottom and curved)
facets = geom.facetCylin
O.
(3) Engines
O.
(4) Time step for simulation:
O.
(5) Check Unbalanced forces:
if O.iter<5000: return
if unbalancedForce
(6) Now as in the experiment for consolidation, I wish to load the top plate with a constant stress/ force for a certain period of time
Can I define the top wall/facet of the facet cylinder as the loading plate? If yes, how ?
Say. the object is named as "plate". In order to subject this plate to a certain stress/force for a certain period of time, do I use the following:
if O.time < 3600 ( say the time for the experiment is 1 hour)
preStress= 2e4
def addForces():
for f in facets:
n = f.shape.normal
a = f.shape.area
O.forces.
the above code applies for all the facets uniformly. How do I apply the above only for the top facet (which acts as the loading plate) and keep the rest stable?
I think I have to do something like:
for f in facets:
f.state.
f.state.mass = 0.5
(Uncertain about the application here)
(7) Call Pyrunner to save the data for the plotting:
O.
(8) Defining the data to be saved and plotted:
I intend to plot the displacement of the top plate against square root of time (in mins)
I think it should be as follows:
import math
def addPlotData():
t= O.time
plot. addData(
plot.plots= { 'Tx':('Dz')}
plot.plot()
O.run
Please correct me if I am wrong anywhere. I am just a beginner :)
Thanks..
Question information
- Language:
- English Edit question
- Status:
- Answered
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Swapnil for more information if necessary.