Water pressure
Hi
I'm working on the saturated triaxial test. I should have positive water pressure in the simulation and need to have its calculation. I assume 2PFV does not appropriate for this aim.
How can I use water pressure in the triaxial test without use of 2PFV?
I mention my dry-triaxial test code in below.
#######
import matplotlib; matplotlib.
import matplotlib.pyplot as plt
from yade import pack
import pylab
from numpy import *
import numpy as np
from math import *
utils.readParam
from yade.params import table
seed=table.seed
num_spheres=
compFricDegree = table.compFricD
confiningS=-1e5
rate=1
finalFricDegree1=45
## creat a packing with a specific particle side distribution (PSD)
psdSizes,
psdSizesArray=
psdSizesMeter=
sp=pack.
mn,mx=Vector3(
sp.makeCloud(
key='_finalFric
## create material #0, which will be used as default
O.materials.
O.materials.
## create walls around the packing
walls=aabbWalls
wallIds=
O.bodies.
triax=TriaxialS
internalCompac
goal1=confiningS,
goal2=confiningS,
goal3=confiningS,
max_vel=10,
label="triax"
)
newton=
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
),
GlobalStiffnes
triax,
newton
]
while 1:
O.run(1000,True)
unb=unbalance
if unb<0.01 and abs(triax.
break
print ("Particle Distribution is Stable")
print ('Porosity:
#######
## REACH NEW EQU. STATE ###
#######
finalFricDegree = finalFricDegree1 # contact friction during the deviatoric loading
#We move to deviatoric loading, let us turn internal compaction off to keep particles sizes constant
triax.internalC
# Change contact friction (remember that decreasing it would generate instantaneous instabilities)
setContactFrict
while 1:
O.run(1000,True)
unb=unbalance
if unb<0.001 and abs(triax.
break
triax.depth0=
triax.height0=
triax.width0=
print ("Time for the deviator loads")
print ('Porosity:
#######
## Start test ###
#######
AxialstrainMatr
DeviatoricStres
VolumetricStrai
PrincipalStress
MeanStressMatri
triax.stressMask=5
triax.goal2=-rate
triax.goal1=
triax.goal3=
file=open(
while 1:
O.run(1000,True)
print ('eyy='
if triax.strain[1] <-0.05:
break
#######
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 Hossein for more information if necessary.