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