修正剑桥模型程序解读

xldpig

多谢分享！！

xldpig

谢谢~~~分享~~

zxb911

楼主，关于剑桥您还有研究吗

chenpf09

学习了。。。。

烟火里的尘埃

我在运行着段命令的时候，最后提示某一个zone requires density if gravity is none-zero，请问这是什么情况

prosper9715

thank you,hao dongxi  yiding  yao  fenxiang

科费卡

请问在用修正剑桥时候，按你的算例怎么总是出现
*** Bad conversion of parameter number 2 (p_ratio).
源代码如下：


new
gen z b p0 0 0 188.35 p1 2000 0 188.35 p2 0 2000 188.35 p3 0 0 200.5 size 100 100 2 ratio 1 1 1 group diceng1
gen z b p0 0 0 175 p1 2000 0 175 p2 0 2000 175 p3 0 0 188.35 size 100 100 2 ratio 1 1 1 group diceng2
gen z b p0 0 0 165.1 p1 2000 0 165.1 p2 0 2000 165.1 p3 0 0 175 size 50 50 1 ratio 1 1 1 group diceng3
gen z b p0 0 0 154.8 p1 2000 0 154.8 p2 0 2000 154.8 p3 0 0 165.1 size 50 50 1 ratio 1 1 1 group diceng4
gen z b p0 0 0 146 p1 2000 0 146 p2 0 2000 146 p3 0 0 154.8 size 50 50 1 ratio 1 1 1 group diceng5
gen z b p0 0 0 91 p1 2000 0 91 p2 0 2000 91 p3 0 0 146 size 10 10 3 ratio 1 1 1 group diceng6
gen z b p0 0 0 84.7 p1 2000 0 84.7 p2 0 2000 84.7 p3 0 0 91 size 10 10 1 ratio 1 1 1 group diceng7
gen z b p0 0 0 60.5 p1 2000 0 60.5 p2 0 2000 60.5 p3 0 0 84.7 size 5 5 1 ratio 1 1 1 group diceng8
gen z b p0 0 0 49.3 p1 2000 0 49.3 p2 0 2000 49.3 p3 0 0 60.5 size 5 5 1 ratio 1 1 1 group diceng9
gen z b p0 0 0 44.5 p1 2000 0 44.5 p2 0 2000 44.5 p3 0 0 49.3 size 5 5 1 ratio 1 1 1 group diceng10
gen z b p0 0 0 16.1 p1 2000 0 16.1 p2 0 2000 16.1 p3 0 0 44.5 size 1 1 1 ratio 1 1 1 group diceng11
gen z b p0 0 0 10.5 p1 2000 0 10.5 p2 0 2000 10.5 p3 0 0 16.1 size 1 1 1 ratio 1 1 1 group diceng12
gen z b p0 0 0 0 p1 2000 0 0 p2 0 2000 0 p3 0 0 10.5 size 1 1 1 ratio 1 1 1 group diceng13
attach face range z 0 201
Initial density 2000
Fix x range x=-0.1 0.1
Fix x range x=2000
fix y range y=-0.1 0.1
fix y range y=2000
fix z range z=0
set gravity 0 0 -9.8

model cam-clay
;lProperty shear=10.06e6 bulk_bound=30e6
roperty mm=0.75358 lambda=0.070981 kappa=0.007617
;Property mpc=0.453045e6
;Property mp1=1.0 mv_l=2.3

;model fl_iso
;prop perm 2.66e-13
;ini fdens 1000
;ini fmod 2e9
;ini sat 1.0

def install_prop
pnt=zone_head
loop while pnt # null
  OCR=1.2
  abs_sxx=abs(z_sxx(pnt))
  abs_syy=abs(z_syy(pnt))
  abs_szz=abs(z_szz(pnt))
  p0=(abs_sxx+abs_syy+abs_szz)/3.0
  p0_effective=p0-z_pp(pnt)
  q0=sqrt(((abs_sxx-abs_syy)^2+(abs_syy-abs_szz)^2+(abs_szz-abs_sxx)^2)*0.5)
  z_prop(pnt,'mm')=6.0*sin(fai*degrad)/(3.0-sin(fai*degrad))
  temp1=q0/(z_prop(pnt,'mm')*p0_effective)
  ;temp1=q0/(mm*p0_effective)
  pc0=p0_effective*(1.0+temp1^2)*OCR
  v0=1.0+_e0
  z_prop(pnt,'cam_cp')=p0_effective
  z_prop(pnt,'mpc')=pc0
  z_prop(pnt,'poisson')=p_ratio
  z_prop(pnt,'lambda')=_lambda
  z_prop(pnt,'kappa')=_kappa
  z_prop(pnt,'mp1')=_mp1
  z_prop(pnt,'mv_l')=v0+_lambda*ln(2.0*_cu/(z_prop(pnt,'mm')*_mp1))+(_lambda-_kappa)*ln(2.0)
  z_prop(pnt,'bulk_bound')=100*40e6
  ;z_prop(pnt,'bulk_bound')=100*(s_mod+4.0/3.0*s_mod)
  pnt=z_next(pnt)
endloop
end
set p_ratio=0.35 fai=34.5 _lambda=0.14 _kappa=0.01 _mp1=1e3 _e0=1.2 _cu=10e3 OCR=1.2
@install_prop

solve
save diceng4.sav

科费卡

茳上风 发表于 2015-3-30 23:38
补充：trying to divide by zero 可能是你的zone里的应力太小或者没有计算。

请问这个问题怎么解决？

aaa479111747

谢谢分享！

科费卡

茳上风 发表于 2015-3-30 22:04
补充：yield envelope cannot be reach  与OCR 、mp1和cu有关，当OCR较大时，mp1和cu可以取得比较小，当OCR ...

按照你的建议调整之后还是没能解决。

wht3911

受教了  谢谢楼主分享

我想说给你听

学习了，谢谢分享

仿真cy123

谢谢分享，大家进步

xspxspxsp

学习学习，高手这么多

xspxspxsp

茳上风 发表于 2015-4-2 18:16
土层输入参数对变形的影响：
参考徐中华博士论文第二章有，
1、孔隙比e，越大，K也越大，土体卸载时的回弹 ...

我做修正剑桥遇到一个问题就是算单层的时候可以算，但是一旦算多层的剑桥模型就出问题，现有FLAC书中都没有提多层剑桥模型，一直显示property bad，一直没有搞出来，你对此有过研究么，我刚到论坛没有这里的钱币，所以想先请教一下

xspxspxsp

科费卡 发表于 2016-8-1 16:26
请问在用修正剑桥时候，按你的算例怎么总是出现
*** Bad conversion of parameter number 2 (p_ratio).
源 ...

赋值前都要加@的

xspxspxsp

科费卡 发表于 2016-8-1 16:46
请问这个问题怎么解决？

我是施加一个数值为1的力，就行了，这句话说分母不能为零

暗度陈仓

急需，谢楼主

暗度陈仓

好好研究下

zaiyunduan

科费卡 发表于 2016-8-1 16:26
请问在用修正剑桥时候，按你的算例怎么总是出现
*** Bad conversion of parameter number 2 (p_ratio).
源 ...

楼主 你这个问题解决了吗