剑桥模型提示平均有效应力无效
本帖最后由 gengchsheng 于 2014-5-7 10:34 编辑命令及参数如下，请前辈赐教
;已知土性参数：质量：M=1,42，泊松比：u=0.3,其他参数：λ/(1+e0)=0.0508，κ/(1+e0)=0.0112
;采用剑桥模型和修正剑桥模型分别绘制下述条件三轴压缩与三轴拉伸应力路径固结排水下土的σ-ε关系
;σ1,σ3分别为最大,最小主应力
;第一种情况：主应力σ1=196KPa；第二种情况:主应力σ3=196KPa；第三种情况：平均主应力p=196KPa；
;初始状态的孔隙比取e0 = 0.68，经计算得κ=0.019，λ=0.085
new
title
σ3=196kpa常数三轴压缩
gen zon brip0 0 0 0 p1 7.5 0 0 p2 0 7.5 0 p3 0 0 15 size 10 10 10;三轴试验尺寸：75*75*150 mm
model cam-clay
prop shear 250bulk_bound 10000
prop mm 1.42 lambda 0.085 kappa 0.019 poisson 0.3 mp1 1. mv_l 3.32;赋予材料物理力学参数
fix x y z range z -0.1 0.1;模型底部边界的竖直向速度约束
def sigma3;定义一个sigma变量
sigma3=-196e3
end
;施加围压σ3=196kpa
ini sxx sigma3 ran x -0.1 0.1
ini sxx sigma3 ran x 7.4 7.6
ini syy sigma3 ran y -0.1 0.1
ini syy sigma3 ran y 7.4 7.6
ini szz sigma3 ran z 14.9 15.1
def install_prop
pnt=zone_head
loop while pnt # null
abs_sxx=abs(z_sxx(pnt)) ;|sxx|
abs_syy=abs(z_syy(pnt)) ;|syy|
abs_szz=abs(z_szz(pnt)) ;|szz|
p0_effective=(abs_sxx+abs_syy+abs_szz)/3.0-z_pp(pnt) ;计算有效应力
z_prop(pnt,'cam_cp')=p0_effective
pnt=z_next(pnt)
endloop
end
install_prop
def axi_stress
p_gp=gp_near(0,0,15)
loop n(1,100)
zss_load= sigma3 - float(n)*10e3;加载1000kPa
command
ini szz zss_load ran z 14.9 15.1
ini sxx sigma3 ran x -0.1 0.1
ini sxx sigma3 ran x 7.4 7.6
ini syy sigma3 ran y -0.1 0.1
ini syy sigma3 ran y 7.4 7.6
solve
endcommand
z_dis = -1.0 * gp_zdisp(p_gp) / 15;得到应变
z_load = (sigma3 - zss_load)/sigma3 ;得到主应力比
command
tab 1z_dis z_load ;保存应变和主应力差
endcommand
end_loop
end
axi_stress
你的提示是不是这样的！ 是不是这样的！
页:
[1]