Precrack
Hi,
I want to put a precrack in my model. Currently, my model consists of 5 by 6 particles in the xy plane. Each particle is linked with its neighbouring particle. I want to zero out the bonds (a precrack) between particle 11 and 16, 12 and 17, and 13 and 18. How do I do this?
Thanks!
MY CODE:
#MATERIAL PROPERTIES
intR = 1.00 #connectivity parameter
DENS = 2500.00 #Density (kg/mm^3)
YOUNG = 1800.00 #Young's modulus (N/mm^2)
FRICT = 7.00 #Friction angle (degrees)
ALPHA = 0.10 #Poisson's ratio
TENS = 100.00 #Tensile strength (N/mm^2)
COH = 100.00 #Shear strength (N/mm^2)
iterMax = 1.00 #Max Itteration
rate = 0.00 #pull rate
O.materials.
# PARTICLES
O.bodies.append([
sphere(center=( 1.00, 1.00,0),radius= 0.50,material=
sphere(center=( 2.00, 1.00,0),radius= 0.50,material=
sphere(center=( 3.00, 1.00,0),radius= 0.50,material=
sphere(center=( 4.00, 1.00,0),radius= 0.50,material=
sphere(center=( 5.00, 1.00,0),radius= 0.50,material=
sphere(center=( 1.00, 2.00,0),radius= 0.50,material=
sphere(center=( 2.00, 2.00,0),radius= 0.50,material=
sphere(center=( 3.00, 2.00,0),radius= 0.50,material=
sphere(center=( 4.00, 2.00,0),radius= 0.50,material=
sphere(center=( 5.00, 2.00,0),radius= 0.50,material=
sphere(center=( 1.00, 3.00,0),radius= 0.50,material=
sphere(center=( 2.00, 3.00,0),radius= 0.50,material=
sphere(center=( 3.00, 3.00,0),radius= 0.50,material=
sphere(center=( 4.00, 3.00,0),radius= 0.50,material=
sphere(center=( 5.00, 3.00,0),radius= 0.50,material=
sphere(center=( 1.00, 4.00,0),radius= 0.50,material=
sphere(center=( 2.00, 4.00,0),radius= 0.50,material=
sphere(center=( 3.00, 4.00,0),radius= 0.50,material=
sphere(center=( 4.00, 4.00,0),radius= 0.50,material=
sphere(center=( 5.00, 4.00,0),radius= 0.50,material=
sphere(center=( 1.00, 5.00,0),radius= 0.50,material=
sphere(center=( 2.00, 5.00,0),radius= 0.50,material=
sphere(center=( 3.00, 5.00,0),radius= 0.50,material=
sphere(center=( 4.00, 5.00,0),radius= 0.50,material=
sphere(center=( 5.00, 5.00,0),radius= 0.50,material=
sphere(center=( 1.00, 6.00,0),radius= 0.50,material=
sphere(center=( 2.00, 6.00,0),radius= 0.50,material=
sphere(center=( 3.00, 6.00,0),radius= 0.50,material=
sphere(center=( 4.00, 6.00,0),radius= 0.50,material=
sphere(center=( 5.00, 6.00,0),radius= 0.50,material=
])
nbSpheres = 30.00
# BOUNDARY CONDITIONS
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
O.bodies[
# FUNCTIONAL COMPONENTS
def addPlotData1():
plot.addData(
# PRINTING
from yade import plot
from pprint import pprint
# SIMULATION LOOP
O.engines=[
PyRunner(
ForceResetter(),
InsertionSortCo
InteractionLoop(
[Ig2_Sphere_
[Ip2_JCFpmMat_
[Law2_ScGeom_
),
GlobalStiffness
NewtonIntegrato
]
#TIME STEP
O.dt=0.
#PLOTTING
plot.plots=
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 Najm for more information if necessary.