1. 渗流示例
;------------------------------参数部分
bulk modulus, K 390 MPa 体积 模量（土质比较硬）
shear modulus, G 280 MPa 剪切模量
soil dry density, ρd 1200 kg/m3 土的密度
water density, ρw 1000 kg/m3 水的密度
wall density, ρwal 1500 kg/m3 档土墙的密度
permeability, k 10?12 m2/Pa-s 渗透系数
porosity, n 0.3 孔隙率
fluid bulk modulus, Kf 2.0 GPa
;-------------------------------------
;
-----以下是命令流
;------------------------------------------------------
; Excavation in a saturated soil
;------------------------------------------------------
config fluid ;
; --- geometrical model --- 建立模型
gen zone brick p1 18 0 0 p2 0 18 0 p3 0 0 18 size 18 18 18 rat 1 1 1
group soil ；设置土组
作弊检讨group excavate range x 0 4 y 0 4 z 0 5 ；设置开挖范 围组
group wal1 range x 4 5 y 0 5 z 0 7 ;
group wal2 range x 0 4 y 4 5 z 0 7 ;
group wall range group wal1 any group wal2 any ；档土墙1 和2 合并为一个档土墙组
本科英语
; --- fluid flow model ---
model fl_iso ；各项同性流体
prop perm 1e-12 poro 0.3 ；设置渗透系数和？？比 or ？？率
ini fdensity 1e3 ; 初始化水密 度
ini fmod 2.0e9 ftens -1e-3 ; 初始化水的体积模量 及 ？？
model fl_null range group excavate ; 开挖部分没有水流
model fl_null range group wall ；档土墙 不渗水
ini pp 0 grad 0 0 1e4 ；初始化 孔隙水压力
fix pp range z -0.1 0.1 ；约束 应该是整个模型的底部吧？
fix pp range x -0.1 4.1 y -0.1 4.1 z 4.9 5.1 ；开挖部分也给约束住了，而其他的地方不管
性生活多久一次才正常;
--- mechanical model --- 力学模型
model elas ; 弹塑性
prop bul 3.9e6 shea 2.8e6 ；设置体积模量和 剪切模量，这两个必须，有一次我只设置了体积模量，就错的一塌糊涂
model null range group excavate ；力学模型的开挖
ini density 1.2e3 ；初始化密度 ，应该是土的密度
ini density 1.5e3 range group wall ；初始化墙（档土墙）的密度，
fix x range x -.1 .1 ；施加约 束，据说差分法约束0 的时候，必须取一个-0.1<0<0.1，范 围
网络游戏情侣名fix x range x 11.9 12.1 ；施加约束，或者叫做边界条件
fix y range y -.1 .1 ；施加约 束
儿童画画入门fix y range y 11.9 12.1
fix z range z 11.9 12.1
; initial total stress 出世总应力
ini szz 0 grad 0 0 -1.5e4 ;z方向的应力大一些
ini sxx 0 grad 0 0 -1.2e4 ；水平方向的，小一些天龙屯堡古镇
ini syy 0 grad 0 0 -1.2e4 ;
apply nstress 0 grad 0 0 -1.2e4 range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0
apply nstress 0 grad 0 0 -1.2e4 range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0
apply nstress -7.5e4 range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1
; --- tting --- 设置
t gravity 0 0 10 ；设置重力加速度
; --- initial state --- ；初 始化状态
solve force 1 ; check initial equilibrium 求解平衡
;
--- histories --- 记录
t hist_rep 40 ；记录步为40
hist fltime ；记录时间
hist gp pp 0 0 6 ；记录某点的孔隙 水压力
hist gp xdis 4 0 0 ；水平位移
hist gp xdis 4 0 2 ；记
录水平位移
hist gp xdis 4 2 0
hist gp zdis 0 0 5 ；记录竖向位移
hist gp zdis 2 0 5
hist gp zdis 4 0 5
hist gp zdis 2 2 5
hist gp zdis 4 2 5
hist gp zdis 4 4 5
hist gp zdis 10 0 1
hist gp zdis 10 0 2
伊索寓言
;
; --- excavation --- 开挖
t fluid off ；关闭水
; apply pore pressure at walls 在墙上施加孔隙水压力
apply nstress 0 grad 0 0 -1.e4 range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0
apply nstress 0 grad 0 0 -1.e4 range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0
apply nstress -5.e4 range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1
solve ;force 1 ；求解
save exc1.sav ; 保存
;
; --- drainage ---
apply remove nstress ；撤掉刚才的压力
def relaxtup ；定义一个释放 函数的参数
step0 = step
end
relaxtup ；调用
def relax ；定义一个释放函数
rstep = step - step0
if rstep < ncyc then
relax=1.0-(float(rstep)/float(ncyc))
el
relax = 0.0
endif
end ；结束
t ncyc = 1000 ；先赋值，随 后调用
apply nstress 0 grad 0 0 -1.e4 hist relax &
range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0
apply nstress 0 grad 0 0 -1.e4 hist relax &
range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0
apply nstress -5.e4 hist relax &
range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1
cyc ncyc
solve ；求解
save exc2.sav ；保存
;
; --- percolation --- 渗透
t fluid on ；设置水算法
fix pp 0 range x -0.1 4.1 y -0.1 4.1 z 4.9 5.1 ；施加孔隙水压力
cyc 9000 ；循环
save exc3.sav ；保存
plot create excav ；建立一个 显示的视图 一下内容都是显示内容
plot t rot 200 0 195
plot t cent 6 6 6
plot t dist 39.18
plot t magn 0.8
plot t plane ori 0 0 0 normal 0 -1 0
爱你是我做过最好的事plot add cont pp out on
plot add block group lgra range group wall
plot add cont pp int 10000 max 110000 range x 0 4 y 0 4 z 5 7
plot add flow plane
plot show
ret
