剑桥模型原存固结应力的计算方法讨论

guihonglong

cam-clay模型计算时需要先生成初始应力场与初始的固结应力pc0。ycding、 math等人都提供过一个计算pc0的方法。思路是按单元循环分别计算其p、q值，再计算得到pc0。该方法避免了手动计算pc0值，只需提供一些基本的参数，即可建模计算，而且考虑到了pc0沿深度的变化，是个很好的思路。我曾把这个方法试用于手册中的路基填筑的例子。手册中的pc0赋值的方法是给一个定值160kpa。与两们大牛的赋值方法，与手册中的pc0明显要大一些。采用ycding与math的方法赋值后，则无法完成计算，提示有效应力为负。对此非常不理解，想向两位大牛及各们高手讨教是何原因。

guihonglong

;-----------------------------------------------------------
; emc3.dat
; Loading by an embankment on a Cam-clay foundation
;-----------------------------------------------------------
new
title
loading by an embankment on a Cam-clay foundation
config fluid
gen zone brick p0 (0,0,0) P1 (20,0,0) P2 (0,1,0) P3 (0,0,10) &
size 20 1 10
; --- model and properties ---
model cam-clay
prop poisson .3 bulk_bound 5e6
prop mm 0.888 lambda 0.161 kappa 0.062
prop mp1 1e3 mv_l 2.858
;prop mpc 160e3
;
ini density 2000
;
model fl_iso
set fluid biot off
prop perm 1e-12 por .3
ini fmod 2e7 ftens -1e2 fdens 1e3
; --- fish function ---
def camclay_ini_p
   nccz = 0
   p_z = zone_head
   loop while p_z # null
   if z_model(p_z) = 'cam-clay' then
      
     s1=z_sxx(p_z)
     s2=z_syy(p_z)
     s3=z_szz(p_z)

     p0 = -(s1+s2+s3)/3.0 - z_pp(p_z)
     z_prop(p_z,'cam_cp') = p0
     
                 q0=sqrt(((s1-s2)^2+(s2-s3)^2+(s3-s1)^2)*0.5)
                 temp1=q0/(z_prop(p_z,'mm')*p0)
                 pc0=p0*(1.0+temp1^2)*1        ;先期有效固结压力，用于确定屈服面
                 z_prop(p_z,'mpc')=pc0
                 ;z_prop(p_z,'mpc')=160e3

   end_if
   p_z = z_next(p_z)
   endloop
end
def ramp
   ramp = min(1.0,float(step)/_nsteps)
end
; --- initial conditions ---
ini sat 1
ini sxx -1.6e5 grad 0 0 1.6e4
ini syy -1.6e5 grad 0 0 1.6e4
ini szz -2.3e5 grad 0 0 2.3e4
ini pp 1e5 grad 0 0 -1e4
;(the CamClay property of mean effective stress must be initialized)
camclay_ini_p
;pause
; --- boundary conditions ---
fix x range x=0
fix x range x=20
fix y
fix x y z range z=0
fix pp 0 range z=10
set _nsteps=5000
apply nstress=-5e4 his ramp range x 0 4 y 0 1 z 9.9 10.1
; --- settings ---
set gravity 0 0 -10
set mech on fluid off
; --- histories ---
hist nstep 100
hist unbal
hist ramp
hist gp zdisp 0 0 10
hist gp zdisp 2 0 10
hist gp zdisp 4 0 10
hist gp zdisp 6 0 10
hist zone szz 0.5 0.5 9.5
hist zone sxx 0.5 0.5 9.5
hist zone pp 0.5 0.5 9.5
hist zone pp 1.5 0.5 7.5
; --- undrained response ---
plot disp
step _nsteps
solve
save emc0.sav
; --- drained response ---
set fluid on
hist fltime
set mech force 0 ratio 5e-5
set fluid substep 100
set mech substep 10 auto
solve age 1e5
set mech substep 50 auto
solve age 1.e8
save emc.sav
plot his 3 4 5 6 vs 11
ret

guihonglong

自己顶一下。
FLAC中用剑桥的人似乎还是有点少啊。参数取值确实困难一些。希望大家一起讨论。

tanglj

你好 请问你解决了这个问题了吗？

lxh684536

我也碰到这个问题了

jiaochiyu

我也碰到同样的问题，不过提示的是，应力为0，无法开方，我想对多层土赋修正剑桥模型参数可是没法实现啊。