已知E, poisson,如何确定kn,ks,fric, pbond参数？

mosquato

看到一些文献上，给出已知E,v，估算kn,ks的公式，但我试了一下，PFC2D仿真出来的E,v和已知数值相差很大。
请教各位高手，究竟如何确定PFC2D中的参数啊？

贴上几种计算kn,ks,pb(kn,ks)的方法，大家看看，为什么我就算不对呢？
method1:
kn=4*r*Eb  ;kn for ball
Kn=2Ec/r    ;kn for pbond
wherein, Eb---ball's Young's Modul, Ec--Pbond's Young Modul...

method2:
kn=3*E*V/(5*Nc*r^2*(2v-1))
kt=((1-2v)/(1+3v))*kn
wherein: E--Young's Modul..; Nc--total number of contacts, v--poisson ratio; V--tot_volume

an example to evaluate method1:  (given parameters as follows
dens=2130
Eb=Ec=64Gpa
Kn/ks=2.0
pb_kn/pb_Ks=2.0
fric=0.65
pb_n=pb_s=124Mpa
;
PFC simulation result:
E=54.2Gpa,  uniaxial comressive strength=100.6Mpa
BTW:这个结果我就是算不出来，大家帮看看，到底是怎么回事？
多谢多谢喽~~~

mengjing120

method1: 是理想的规则组成的理论解，仅作理解参数关系用
method2:不知道你是怎么来的，
没有直接的办法，因为各个颗粒随机生成，组合方式不同，无法建立关系，可以不断调试凑出来

wangpeitao

确定弹模方法的前提是颗粒粒级组成吧，不同颗粒簇对应的弹性模量应该有差别，类似弹簧或者平行键的串并联的结果，这个应该不会有一个统一的确定方法，视颗粒材料而定。

mosquato

我看到FISHtank里面，有一部分是说如何根据宏观参数，确定离散元细观参数的。
但是，我没有fishtank--Biaxial test的*.drv，*.fis文件，
自己试着遍了一个unconfined compression test 程序，但是结果似乎不合理啊？
请教各位，谁有不限制颗粒数的fishtank--biaxial/brizillian模拟程序啊！
多谢喽~~~