erase rigid boundaries and add periodic boundaries
Hi,
I met a question when building the model.
I have two scripts and they should run in sequence.
The first script is used to generate particles. The particles expand until the wallstress reach 90kPa. The walls are rigid.
In the second script, I want to erase the rigid walls and replace them by the periodic boundaries, but I don’t know how to add the periodic boundaries after executing the following command
###erase the six rigid boundaries###
for b in range(6):
O.bodies.erase(b)
The first script:
#######
### Defining parameters and variables ###
#######
#Material constants
Density = 3000
FrictionAngle = 1.5
PoissonRatio = 0.5
Young = 300e6
Damp = 0.5
AvgRadius = 0.1
N_particles = 4000
#Wall constants
WDensity = 0
WFrictionAngle = 0.0
WPoissonRatio = 0.5
WYoung = 50e9
#Packing variables
mn = Vector3(1,0,0)
mx = Vector3(101,20,1)
#Confining variables
ConfPress1 = -90000 #pre-compression
ConfPress = -1.0e5
#time calculation
startT = O.time
endT = O.time
timeSpent = endT - startT
#######
#import necessary packages
from yade import pack,plot,os,timing
import matplotlib; matplotlib.
import pylab
#######
### Sample Preparing ###################
#######
#Create materials for spheres and plates
SphereMat = O.materials.
WallMat = O.materials.
#Create walls for packing
wallIds = O.bodies.
sp = pack.SpherePack()
sp.makeCloud(
sp.toSimulation
O.usesTimeStepp
O.trackEnerty=True
#######
#####Defining triaxil engines####
#######
###first step: compression#######
triax1=
wall_
wall_
thickness = 0.001,
maxMultipli
finalMaxMul
internalCom
stressMask = 7,
computeStre
goal1 = ConfPress1,
goal2 = ConfPress1,
)
newton=
###engine
O.engines=[
ForceResett
InsertionSo
Interaction
[Ig2_
[Ip2_
[Law2_
),
GlobalStiff
triax1,
newton,
PyRunner(
]
# Simulation stop conditions defination
def checkUnbalanced():
unb=
if unb<0.001:
O.pause()
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Jan Stránský
- Solved:
- Last query:
- Last reply: