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. _______________________________________________ Mailing list: https://launchpad.net/~yade-users Post to : yade-users@lists.launchpad.net Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp