yade-users team mailing list archive

[jsonscript <question691560@xxxxxxxxxxxxxxxxxxxxx>]

Question #691560 on Yade changed:
https://answers.launchpad.net/yade/+question/691560

    Status: Answered => Open

jsonscript is still having a problem:
Hello Jan,

I realise I wasn't clear enough in my previous comments. I apologise for
that.

Here is the MWE script (Runtime 2-3 mins) :

#-----------------------------------------------#
from yade import pack
from yade import plot

O.periodic=True
O.cell.setBox(.3,.3,.3)
sp=pack.SpherePack()
radius=5e-3
num=sp.makeCloud(minCorner=(0,0,0),maxCorner=(.3,.3,.3),rMean=radius,rRelFuzz=0.5,num=10000,periodic=True)

O.bodies.append([sphere(s[0],s[1]) for s in sp])

O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb()]),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom()],
		[Ip2_FrictMat_FrictMat_FrictPhys()],
		[Law2_ScGeom_FrictPhys_CundallStrack()]
	),
	PeriTriaxController(dynCell=True,maxUnbalanced=2e-3,relStressTol=2e-4,goal=[-1e5,-1e5,-1e5],stressMask=7,globUpdate=5,maxStrainRate=[1.,1.,1.],doneHook='triaxDone()',label='triax'),
	NewtonIntegrator(damping=.2)
]

def triaxDone():
  print ('Here we are: stress',triax.stress,'strain',triax.strain,'stiffness',triax.stiff)
  for i in O.interactions:
    n=i.geom.normal
    f=open("normals.txt","a")
    print ("%s %s %0.9f %0.9f %0.9f" % (i.id1, i.id2, n[0], n[1], n[2]), file=f)
  f.close()
  O.pause()

O.dt=PWaveTimeStep()
O.run()
O.wait()
#----------------------------------------------------------#

I tried to consolidate a specimen isotropically with 100kPa and  later checked for the isotropy in contact-normal distribution.
>From the contact-normals recorded in the text file "normals.txt", I plotted the spherical histogram using Matlab. It seems to me, the distribution is not even near to isotropic.

I couldn't attach my histogram plot here, respecting "No-external link"
rule of this forum. But, I see the distribution is not isotropic.

Am I doing something wrong here? if not, is there some other preparation
method to generate isotropic sample w.r.t. both stress and contact
normals?

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