New question #680414 on Yade:
https://answers.launchpad.net/yade/+question/680414
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
"PeriTriaxController run on aperiodic simulation". I set up the FCC structure
in the rock sample by specifying the centroid position and radius of every
sphere particle, however I'm 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.materials.append(FrictMat(young=15e6,poisson=.4,frictionAngle=0,density=0,label='frictionless'))
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:
sphereArr.append((L*i,m*L/2,j*L/2))
else:
sphereArr.append((L*(i+(.5)),m*L/2,j*L/2))
O.bodies.append([sphere(center, r,material='spheres',color = (0.42,0.06,0.51)
)for center in sphereArr])
## create walls around the packing
walls=aabbWalls([(0,0,0),(20,30,20)],thickness=1e-10,material='frictionless')
wallIds=O.bodies.append(walls)
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()]
),
PeriTriaxController(label='triax',
goal=(sigmaIso,sigmaIso,sigmaIso),stressMask=7,
dynCell=True, maxStrainRate=(10,10,10),
# wait until the unbalanced force goes below this value
maxUnbalanced=.1,relStressTol=1e-3,
# call this function when goal is reached and the packing is stable
doneHook='compactionFinished()'
),
NewtonIntegrator(damping=.2),
PyRunner(command='GlobalData()',iterPeriod=500),
]
O.dt=.5*PWaveTimeStep()
def GlobalData():
print 'porosity '+utils.porosity()
print 'coord #' +utils.avgNumInteractions()
def compactionFinished():
# set the current cell configuration to be the reference one
O.cell.trsf=Matrix3.Identity
# change control type: keep constant confinement in x,y, 20%
compression in z
triax.goal=(sigmaIso,sigmaIso,-.2)
triax.stressMask=3
# allow faster deformation along x,y to better maintain stresses
triax.maxStrainRate=(1.,1.,.1)
# next time, call triaxFinished instead of compactionFinished
triax.doneHook='triaxFinished()'
# do not wait for stabilization before calling triaxFinished
triax.maxUnbalanced=10
def triaxFinished():
print 'Finished'
O.pause()
--
You received this question notification because your team yade-users is
an answer contact for Yade.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
