sth about stressmask
Hello,
The following code is an example from Yade. I know the first part is about hydrostatic compression. But how to understand "stressMask=
# peri3dControlle
# script, that explains funcionality and input parameters of Peri3dController on the example of
# triaxial compression.
# Firstly, a hydrostatic preassure is applied, then a strain along z axis is increasing
# while x- and y- stress is constant
# The simulation is run on rotated cell to enable localization and strain softening
# (you can also try simulation with command sp.toSimulation() with no rotation,
# in this case there is almost no difference, but in script peri3dControlle
# the cell rotation has significant effect)
from yade import pack,plot,qt
# define material
O.materials.
# create periodic assembly of particles
initSize=1.2
sp=pack.
angle=0
rot=Matrix3(
sp.toSimulation
# plotting
plot.live=False
plot.plots=
def plotAddData():
plot.addData(
iter=
sx=p3d.
syz=p3d.
ex=p3d.
eyz=p3d.
)
O.dt=PWaveTimeS
# define the first part of simulation, hydrostatic compression
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
NewtonIntegrat
Peri3dController( goal=(-
stressMask
nSteps=500,
doneHook=
youngEstim
maxStrain=.5,
label='p3d'
),
PyRunner(
]
O.run(); O.wait()
# second part, z-axis straining and constant transversal stress
O.engines=[
ForceResetter(),
InsertionSortC
InteractionLoop(
[Ig2_
[Ip2_
NewtonIntegrat
Peri3dController( goal=(-
nSteps=1000,
xxPath=
yyPath=
doneHook='print "Simulation with Peri3dController finished."; O.pause()',
maxStrain=.5,
label='p3d',
strain=
stressIdeal=
),
PyRunner(
]
O.run();O.wait()
plot.plot(
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 LuBaoBao for more information if necessary.