| 飞鱼游鹰 |
2010-03-16 09:08 |
高手帮我看看这个蠕变程序 （FLAC2D）
大家好 以下是我编的三轴试验蠕变的程序
config axis creep
g 20 20
gen 0 0 0 0.2 0.05 0.2 0.05 0
model cvisc
prop den 1060 dilation 27 bu 1.2e10 she=1.1e10
prop coh 2.72e5 friction 44 ten 2e5
prop k_visc=1.1e10 k_shear_mod=1.1e10
def confine ;定义围压
confine=-con
app_x=con
end
set con=1.5e5 ;给围压赋初值(400 kPa)
def load ;定义轴向压力
sum2=0.0
loop i (1,igp)
sum2=sum2+yforce(i,jgp)
end_loop
lo=sum2/(x(igp,jgp)-x(1,jgp))
sdiff=lo-con ;定义偏差应力
load=lo
end
def vstrain ;定义轴向应变
vst = -ydisp(1,jgp) / y(1,jgp)
if vst < 0.0 then
vst = 0.0
end_if
vstrain = vst
end
def superstep
loop n (1,n_steps)
app_y=n*(5e4)
command
app syy app_y j=jgp
step 10000
end_command
end_loop
end
set n_steps=8
fix y j 1 ;固定下边界y方向的位移
fix y j 21 ;固定上边界y方向的位移
ini yvel -1e-7 j=21;施加上边界y方向速度
ini sxx confine syy confine szz confine ;初始化单元应力
apply pressure app_x i 11 ;施加右边界应力(400 kPa)
set crdt 2.592
step 1000000
hist ydis i=1,j=21
hist vstrain
hist crtime
plot hist 2 vs 3
plot hist 1 vs 3
我做的是三轴的蠕变试验 采用分级加载的（每级荷载下维持72个小时）方法 这个程序怎么执行一下就没有了呢？我编的有什么问题？我的QQ 914108739 邮箱是qin20054184@126.com 麻烦高手指点 ^_^
|
|
