New question #707087 on Yade:
https://answers.launchpad.net/yade/+question/707087
I was conducting a triaxial test on the argillite, which contained a matrix and
two inclusions, to investigate its small strain stiffness. I tried to use
spherical grains and JCF model to simulate this kind of rock. And I seperated
the particle into three groups and assigned them different material properties
(Young's modulous and Poisson's ratio). But the results were not desired.
Because the small strain stiffness of the sample is almost equivalent to the
sample containg only one material, which is confusing. Theoretically, the
change of materials in the sample will affect the properties of the rock. How
can I realize this goal? Could you give me some advice?
from yade import ymport, utils , plot
O.load('isocompact200kpa2.yade.bz2')
O.engines[6].dead=True
# ENGINES ARE DEFINED HERE
try:
os.mkdir('post')
except:
pass
PACKING='121_1k.spheres'
OUT=PACKING+'_1MPa_r0.02'
iterMax=5 # maximum number of iterations
saveVTK=iterMax/5 # Vtk files record interval
YOUNG2=100e9
YOUNG3=50e9
ALPHA2=0.06
ALPHA3=0.27
bodies2=[b for b in O.bodies if 5.0 0.001 and e22 < 0.05:
# O.save('save'+'{:.4f}'.format(e22)+'yade.bz2')
# ini_e22a = e22
#if (e22 - ini_e22b) > 0.5:
print ('unbalanced force: %f, mean stress: %f, s11: %f, s22: %f, s33: %f,
coordination number: %f, porosity: %f'
%(unb,mStress,-triax.stress(triax.wall_right_id)[0],-triax.stress(triax.wall_top_id)[1],-triax.stress(triax.wall_front_id)[2],avgNumInteractions(),Porosity))
#Output()
#ini_e22b = e22
plot.saveDataTxt('loadinglog3.txt.bz2')
plot.plots={'e22':('s22',None,'ev')}
plot.plot()
#O.run(iterMax)
--
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