Creating a FCC structued packing and applying confining stress

Asked by Anqi H on 2019-04-24

Hey everyone!

I've been trying to create a cubic rock sample of FCC structure, and apply a 0.5 MPa confining stress to prepare the sample for other simulations. I tried following the periodic triaxial testing script but something is not correct in my setup as I'm getting an error message "body larger than half of the cell size encountered". I set up the FCC structure in the rock sample by specifying the centroid position and radius of every sphere particle, however I'm not sure if this is the correct way to create the packing. I have attached the my code below, if anyone could show me some pointers I would really appreciate it.

from yade import pack, qt, plot

r = 0.45
L = 2*sqrt(2)*r
sigmaIso = -5e5

# rock material
rockID = O.materials.append(CpmMat(young=15e6,poisson=.4,frictionAngle=radians(atan(0.5)),density=2650,label='spheres'))
#wall material
O.periodic = True

sphereArr = []

for m in range(0,46): #y axis
 for j in range(0,30): #z axis
  for i in range(0,15): #x axis
   if j%2==0:

O.bodies.append([sphere(center, r,material='spheres',color = (0.42,0.06,0.51) )for center in sphereArr])

## create walls around the packing

 dynCell=True, maxStrainRate=(10,10,10),
 # wait until the unbalanced force goes below this value
 # call this function when goal is reached and the packing is stable


def GlobalData():
 print 'porosity '+utils.porosity()
 print 'coord #' +utils.avgNumInteractions()

def compactionFinished():
 # set the current cell configuration to be the reference one
 # change control type: keep constant confinement in x,y, 20% compression in z
 # allow faster deformation along x,y to better maintain stresses
 # next time, call triaxFinished instead of compactionFinished
 # do not wait for stabilization before calling triaxFinished

def triaxFinished():
 print 'Finished'

Question information

English Edit question
Yade Edit question
No assignee Edit question
Solved by:
Anqi H
Last query:
Last reply:
Jérôme Duriez (jduriez) said : #1


A first step would to be to correct O.periodic = Trus ?

Anqi H (analoq) said : #2

You're right Jerome, casualties of typing code intead of copying. The script I run doesn't have this typo.

Anqi H (analoq) said : #4

Hi Robert, pack.regularOrtho is what I'm looking for, thank you for your help!