Axial stress does decrease at the beginning of the deviatoric stress loading in triaxial test
Hi,
I am conducting a simulation of particle breakage in a triaxial test, but at present I'm facing a strange problem: Normally, after the consolidation is completed, the deviatoric loading is entered, but deviatoric pressure is negative at the initial stage, which means that the axial pressure at this stage is less than the confining pressure. In fact, for example, the set confining pressure is 100 kPa, and after the consolidation is completed, the pressures in all three directions are both 100 kpa, then entering deviatoric loading, the axial stress should be increased from 100 kpa. But in my test, after consolidation, confining pressure is 100 kpa in three directions, but in deviatoric loading , axial stress is changed to 80 kpa , as the axial strain increases, the axial stress will gradually increase. When the axial strain is about 5%, the deviatoric stress will be greater than 0.
Do you have any idea about why deciatoric stress is negative in the initial phase of the deviatoric loading?
I urgently need your help . Thanks advance!
Regards!
Cloud
this is my MWE(MWE1.py and MWE2.py are two python file, run MWE1.py first ):
# MWE1.py
from yade import pack,export,
import random
mn=(0,0,0)
mx = (40e-3,40e-3,80e-3)
sp=pack.
sp.makeCloud(
sp.toSimulation()
export.
import random
random.seed(1)
sp = ymport.
colors = [randomColor() for s in sp]
def Children(b):
p = b.state.refPos
r = b.shape.radius
s1 = sphere(
s2 = sphere(
s3 = sphere(
s4 = sphere(
s5 = sphere(
return s1,s2,s3,s4,
for si,i in enumerate(sp):
oriBody = Quaternion(
oriBody *= Quaternion(
oriBody *= Quaternion(
b=O.bodies[si]
children = Children(b)
O.bodies.erase(si)
a=bodiesHandli
ids=O.
for j in ids:
O.bodies[
export.
#MWE2.py
from yade import ymport
key='haxi_'
mn = (0,0,0)
mx = (40e-3,40e-3,80e-3)
young = 300e6
compFricDegree = 30
damp = 0.6
stabilityThreshold = 0.01
rate = -0.02
finalFricDegree = 20
targetPorosity = 0.45
poisson = 0.22
confiningS = 100e3
O.materials.
young=young,
poisson=
frictionAng
density=0,
label='walls')
)
walls = aabbWalls(
wallIds = O.bodies.
O.materials.
young=young,
poisson=
isCohesive=
frictionAng
density=2650,
momentRotat
normalCohes
shearCohesi
alphaKr=0.2, # rolling stiffness
label=
)
attrs = []
sp = ymport.
n = max(int(a[0]) for a in attrs)+1
colors = [randomColor() for _ in xrange(n)]
for s,a in zip(sp,attrs):
aa = int(a[0])
s.agglomerate = aa
s.shape.color = colors[aa]
O.bodies.append(sp) # the last defined material is usedself.
triax=TriaxialS
maxMultiplier=
finalMaxMultip
thickness = 0,
stressMask = 7,
internalCompac
)
factor=1.05
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
),
GlobalStiff
triax,
TriaxialSta
NewtonInteg
]
O.run(1,True)
ig2sss.
for i in O.interactions:
if not isinstance(
continue
b1,b2 = [O.bodies[ii] for ii in (i.id1,i.id2)]
if b1.agglomerate != b2.agglomerate:
else:
i.phys.unp = i.geom.
O.run(1,True)
yade.qt.
#######
### APPLYING CONFINING PRESSURE ###
#######
triax.goal1=
while 1:
O.run(1000, True)
unb=unbalanced
print 'unbalanced force:',unb,' mean stress: ',triax.meanStress
if unb<stabilityTh
break
O.save(
print "### Isotropic state saved ###"
triax.width0=
triax.depth0=
triax.height0=
#######
## DEVIATORIC LOADING ###
#######
setContactFrict
triax.stressMask = 3
triax.goal3= rate
triax.goal1=
triax.goal2=
O.saveTmp()
from yade import plot
## a function saving variables
def history():
plot.addData(
s11=
s22=
s33=
devi = -triax.
i=O.iter)
def stoploading():
axial_strain = abs(triax.
if axial_strain > 0.2:
print 'axial strain reach Limitation, stop loading!'
O.pause()
if 1:
O.engines=
O.engines=
##O.run(100,True)
plot.plots=
# display on the screen (doesn't work on VMware image it seems)
plot.plot()
# In that case we can still save the data to a text file at the the end of the simulation, with:
plot.saveDataTx
#or even generate a script for gnuplot. Open another terminal and type "gnuplot plotScriptKEY.
plot.saveGnuplo
Question information
- Language:
- English Edit question
- Status:
- Expired
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply: