;幫助的 第ftd138.pdf
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 流體的體積模量
;-----以下是命令流
config fluid ;設置流體算法
gen zone brick p1 12 0 0 p2 0 12 0 p3 0 0 12 size 12 12 12 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 ; 設置檔土墻 1 組
group wal2 range x 0 4 y 4 5 z 0 7 ;設置檔土墻 2 組
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 ？？率
oini fdensity 1e3 ; 初始化 水密度
ini fmod 2.0e9 ftens -1e-3 ; 初始化水的體積模量 及 ？？
model fl_null range group excavate ; 開挖部分沒有水流
model fl_null range group wall ;檔土墻 不滲水
jini 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 ;開挖部分也給約束住了，而其他的地方不管
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
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.１
; --- setting --- 設置
set gravity 0 0 10 ;設置重力加速度
solve force 1 ; check initial equilibrium 求解平衡
; --- excavation --- 開挖
set fluid off ;關閉水
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.０
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 relaxsetup ;定義一個 釋放 函數的參數
step0 = step
end

relaxsetup ;調用
def relax ;定義一個釋放函數
rstep = step - step0
if rstep < ncyc then
relax=1.0-(float(rstep)/float(ncyc))
else
relax=0.0
endif
end
set ncyc=1000；先賦值，然后調用
apply nstress 0 grad 0 0 -1e4 his relax range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0
apply nstress 0 grad 0 0 -1e4 his relax range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0
apply nstress -5e4 his relax range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1
cyc ncyc
solve


; --- percolation --- 滲透
set fluid on ;設置水算法
Fluentfix pp 0 range x -0.1 4.1 y -0.1 4.1 z 4.9 5.1 ;施加孔隙水壓力
cyc 9000 ;循環
save exc3.sav ;保存


以上代碼并不是全部，我只是挑了一些覺得重要的，?。。?/span>