yade-users team mailing list archive

[fengjingyu <question677986@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi,

I changed the properties of the material while adjusting the isotropic compression to the one-dimensional  compression.I track the stresses and strains on the z-axis.But i found that the strain was changed in the second compression because in a one-dimensional compression, the stress should be increasing.

 (-15104.1021154	-0.322694610286
-15040.6585685	-0.322695012996
-14956.6075828	-0.322696878321
-13491.3917129	-6.25046306966e-09
-10620.4527817	-6.32724609453e-09
-7614.76213004	-6.4040291194e-09
-4974.45004156	-6.48081214428e-09
-2838.2818472	-6.55759516916e-09)

What should I do to inherit the stress?

Thanks，

Feng


###########################
# encoding: utf-8

sigmaIso=-15000

from yade import pack, qt, plot

O.periodic=True

O.materials.append(FrictMat(young=64e9,poisson=0.12,density=2650,frictionAngle=0.0005,label='spheres'))
s=O.materials.append(FrictMat(young=64e9,poisson=0.12,density=2650,frictionAngle=0.24))

sp=pack.SpherePack()
sp.makeCloud((0,0,0),(3,3,3),rMean=.1,rRelFuzz=.3,periodic=True)
n=sp.toSimulation(material='spheres')

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),
      maxUnbalanced=.1,relStressTol=1e-3,
      doneHook='compactionFinished()'
   ),
   NewtonIntegrator(damping=.1),
   PyRunner(command='addPlotData()',iterPeriod=1000),
]
O.dt=.5*PWaveTimeStep()

def addPlotData():
   	plot.addData(z=triax.stress[2],c=triax.strain[2],)
	print(triax.stress[2],triax.strain[2])

def compactionFinished():
   for i in n:
	O.bodies[i].material = O.materials[s]
   O.cell.trsf=Matrix3.Identity
   triax.goal=(sigmaIso,sigmaIso,-.5)
   triax.stressMask=3
   triax.maxStrainRate=(0,0,0)
   triax.doneHook='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.