yade-users team mailing list archive

[SayedHessam <question677516@xxxxxxxxxxxxxxxxxxxxx>]

New question #677516 on Yade:
https://answers.launchpad.net/yade/+question/677516

Hey,


I am trying to implement nbTensCracks and nbShearCracks functions [1] & [2], to count the number of breakage after compression of one single particle. BTW, I faced attribute error, the error is:

AttributeError: 'Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM' object has no attribute 'nbTensCracks'.


[1] https://answers.launchpad.net/yade/+question/635871

[2] https://answers.launchpad.net/yade/+question/676014


the MWE is as follows:

from yade import pack, utils, ymport, export
import numpy as np
import yade.timing
import matplotlib.pyplot as plt
import os.path

particleDiameter = 0.6e-3
displacementRate = 1.6667e-6 # m/s
strainRate = displacementRate/particleDiameter # m/s

young =3e8  #elastic modulus [Pa]
poisson = 0.3
density = 5000
relDuctility=30     #relative ductility of bonds in normal direction
mn,mx=Vector3(-0.52e-3, -1e-3,-1e-3),Vector3(0.52e-3, 1e-3, 1e-3) #size of particle
finalFricDegree = 25  #Contact friction angle
intRadius= 1.329   #Internal Radius 
sigmaT = 7e6   #Normal shear stress
identifier=str(particleDiameter)
output='./out/'+identifier
rate = strainRate*10000
stabilityThreshold=0.01
iterper=1000   #Periodicity criterion using step number (deactivated if <= 0)
cohesion=20e6

jCFmat = O.materials.append(JCFpmMat(young=young, cohesion=cohesion, density=density, frictionAngle=radians(finalFricDegree), tensileStrength=sigmaT, poisson=poisson, label='JCFmat'))

#frictionless walls
wallmat = O.materials.append(FrictMat(young=80e9,poisson=.45,frictionAngle=radians(finalFricDegree),density=7000,label='frictionlessWalls'))

walls=aabbWalls([mn,mx], material='frictionlessWalls')
wallIds=O.bodies.append(walls)

sp = O.bodies.append(ymport.textExt('sandGrain_0-6mmdiameter.spheres', 'x_y_z_r',color=(0,0.2,0.7), material='JCFmat'))

triax = TriaxialStressController(	
	thickness=0,
	stressMask=7,
	internalCompaction=False,
	)


O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=intRadius),Bo1_Box_Aabb()]),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom6D(interactionDetectionFactor=intRadius),Ig2_Box_Sphere_ScGeom()],
		[Ip2_FrictMat_FrictMat_FrictPhys(),Ip2_JCFpmMat_JCFpmMat_JCFpmPhys(cohesiveTresholdIteration=10, label='jcf')],
		[Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM(Key=identifier, smoothJoint=True,label='interactionLaw', recordCracks=True, neverErase=True),Law2_ScGeom_FrictPhys_CundallStrack()]
	),	
	GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.5),
	triax,
	VTKRecorder(dead=0,Key=identifier,iterPeriod=iterper,initRun=True,fileName=(output+'-'),recorders=['spheres','intr','stress', 'bstresses','jcfpm','cracks']),
	NewtonIntegrator(damping=0.9)
]

O.dt = 0.3 * utils.PWaveTimeStep()


from yade import qt
v=qt.Controller()
v=qt.View()

from yade import plot
O.engines=O.engines[0:6]+[PyRunner(iterPeriod=iterper/4,command='stressStrainHist()',label='recorder')]+O.engines[6:7]

def stressStrainHist():
	plot.addData(e11=np.exp(-O.engines[4].strain[0])-1, 
		e22=exp(-(O.engines[4].strain[1]))-1, 
		e33=-O.engines[4].strain[2],
		s11=-O.engines[4].stress(0)[0],
		s22=-O.engines[4].stress(2)[1],
		s33=-O.engines[4].stress(4)[2],
		i=O.iter,
		displacement=(exp(-(O.engines[4].strain[1]))-1)*particleDiameter,
		P = -O.engines[4].stress(0)[0]*O.engines[4].height*O.engines[4].depth)
	plot.saveDataTxt('txt/data'+identifier+'.txt',vars= ('displacement','e11','s11','P'))

plot.plots={'displacement':(('P'))} #,None,('Uk'))}
plot.plot()

def addPlotData():
    yade.plot.addData({'i':O.iter,
         'tc':interactionLaw.nbTensCracks,
         'sc':interactionLaw.nbShearCracks,
         'te':interactionLaw.totalTensCracksE,
         'se':interactionLaw.totalShearCracksE,
         'unbF':utils.unbalancedForce()})
	})
    plot.saveDataTxt('txt/data.txt')


#### coordination number verification
numSSlinks=0
numCohesivelinks=0
for i in O.interactions:
    if isinstance(O.bodies[i.id1].shape,Sphere) and isinstance(O.bodies[i.id2].shape,Sphere):
      numSSlinks+=1
    if i.phys.isCohesive :
      numCohesivelinks+=1

print "nblinks=", numSSlinks, " | nbCohesivelinks=", numCohesivelinks, " || Kcohesive=", 2.0*numCohesivelinks/nbSpheres

O.engines=O.engines[0:7]+[PyRunner(iterPeriod=iterper/4,command='stopifDamaged()',label='recorder')]+O.engines[7:8]

def stopifDamaged():
	sigma,eps, P=plot.data['s11'],plot.data['e11'],plot.data['P']
	if O.iter > 5000:
		if  P[-1]>10 and P[-1] < 0.6*max(P):
			print 'failure reached'
			yade.timing.stats()
			O.pause()

triax.stressMask = 0

triax.goal2=rate
triax.goal1=-rate
triax.goal3=rate

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.