increasing the applied stress every 5000 iteration
Hello,
I realise a triaxial compression test.
I want to increase the major stress every 5000 iterations
For example
0<= i < 5000 : sigma1=2MPa
5000 <= i <10000 : sigma1 = 3MPa
10000 <= i <15000 : sigma1 = 4 MPa etc...
The problem is that my code increase the major stress once and the major stress becomes constant.
Best regards
Jabrane
This is my code:
from yade import ymport, utils , plot
import math
PACKING=
OUT=PACKING+
l
DAMP=0.4
saveData=100
iterMax=2000000
saveVTK=10000
confinement=-2e6
uniaxial_
delta_stress=-1e6
stress_max=-20e6
intR=1.47618
DENS=4000
YOUNG=65e9
FRICT=10
ALPHA=0.4
TENS=8e6
COH=160e6
def sphereMat(): return JCFpmMat(
def wallMat(): return JCFpmMat(
O.bodies.
dim=utils.
xinf=dim[0][0]
xsup=dim[1][0]
X=xsup-xinf
yinf=dim[0][1]
ysup=dim[1][1]
Y=ysup-yinf
zinf=dim[0][2]
zsup=dim[1][2]
Z=zsup-zinf
R=0
Rmax=0
numSpheres=0.
for o in O.bodies:
if isinstance(
numSpheres+=1
R+=o.
if o.shape.
Rmax=
Rmean=R/numSpheres
O.reset()
mn,mx=Vector3(
walls=utils.
wallIds=
O.bodies.
for o in O.bodies:
if isinstance(
o.shape.
O.engines=[
InteractionLoop(
[Ig2_
[Ip2_
[Law2_
),
]
tensCks=
e10=e20=e30=0
def recorder():
global tensCks, shearCks, e10,e20,e30
tensCks=0
shearCks=0
for o in O.bodies:
tensCks+
shearCks+
yade.
,i=O.iter
,e1=
,e2=
,e3=
,s1=
,s2=
,s3=
,tc=
)
plot.
plot.plots=
plot.plot()
O.step();
SSgeom.
Saabb.aabbEnlar
numSSlinks=0
numCohesivelinks=0
numFrictionalLi
for i in O.interactions:
if not i.isReal : continue
if isinstance(
numSSlinks+=1
if i.phys.isCohesive :
numCohesi
else :
numFricti
print "nbSpheres=", numSpheres," | coordination number =", 2.0*numCohesive
triax.stressMask=7
triax.goal1=
triax.goal2=
triax.goal3=
triax.max_vel=0.001
while 1:
if confinement==0:
O.run(
break
O.run(100,True)
unb=unbalance
meanS=
print 'unbalanced force:',unb,' mean stress: ',meanS
if unb<0.005 and abs(meanS-
O.run(
e10=
e20=
e30=
break
O.bodies[
O.bodies[
triax.stressMask=7
triax.goal1=
triax.goal3=
triax.goal2=
while 1:
O.run(5000,True)
triax.
print 'stress', triax.goal2
break
O.run(iterMax)
Question information
- Language:
- English Edit question
- Status:
- Expired
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
This question was reopened
- by Yor1