请ycding博士进来看看指导一下

guihonglong

http://forum.simwe.com/viewthread.php?tid=935935&highlight=%BD%A3%C7%C5

以前的一个贴子，尚未找到解决的办法。请丁博等各位高手指导一下。

guihonglong

或者哪位能上传一个用剑桥模型计算的例子吗？

guihonglong

或者哪位能上传一个用剑桥模型计算的例子吗？

lmm

找丁博算是找对人了

guihonglong

是在论坛上下载了丁博的有关剑桥模型初应力计算的一个命令流试用了一下，结果有些问题。
百思不得其解，所以再向丁博及各位请教一下。请知道的大侠赐教，

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

命令流再发一下。

guihonglong

丁博做点售后服务，人民版主为人民。

guihonglong

看来还是没动静啊。我现在并没有相关的项目，只是突然对此感兴趣，想弄明白。
突发奇想，没事的时间我就隔几天来顶一下这个贴子，并邀请有兴趣的同仁们帮忙顶一下，直到ycding博士版主或哪位高手答复为止。
这应该算是FLAC中的一个比较不容易搞明白的事吧，反正我是没明白。一样与我不明白的请帮着顶一下，搞明白的就答复一下。

guihonglong

完整的描述一下问题
针对路基填筑做一个算例，采取如下方案
第一步，采用mc，流固耦合，完成计算，没问题。
第二步，采用camclay，流固耦合，提示yield envelop can not be reached 或者 effective porepressure   为负（除换本构外不做任何修改）
第三步，采用camclay，不考虑流体，仅做弹塑性分析，完成计算，而且量值与mc固结的情况的结果比较符合。

请各位分析一下camclay情况下出错的可能性？

坚持顶，顶到几位版主或其它大侠们出手为止。
怎么没人响应呢？没有人遇到这样的问题吗？看来是我太笨了？