考虑流固耦合的某海底隧道的渗流分析

nietaoyi

太贵了吧

nietaoyi

下下来看看

nietaoyi

太贵了，建议出售价格2仿真比

anningjhj

keng神坑，超级坑，巨坑

sjzdh

楼主先介绍一下啦

taulia

太贵了。。。。还是不看了吧

sjzdh

下载学习学习

taulia

没钱。。。。。。。。

周毅

群里面早就有的

jk在路上

楼主好人…。。。

zqwyx2012

不错 要看看 不过仿真比不够用啊 亲

zqwyx2012

怎么没有模型 哎 楼主能不能发一个啊

745257116

回复一个，赚点币回来！

a570705800

顶一下，看看自己多少钱了。

whx1002102308

太贵了。。。。。。。。。。。

wystyler

看看是啥

zaijianjack

楼主是个大坑货！没法运行，上来就rest model.sav，尼玛又不把model.sav给出来，真不厚道！
new
rest model.sav
title
Calculation for External water pressure on lining of subsea tunnels
def setup
  c_perm=2.0e-11         ermeability of rock (m/s)/(Pa/m)
  g_perm=4.0e-13         ;Permeability of grouting loop (m/s)/(Pa/m)
  l_perm=2.0e-12         ;Permeability of liner (m/s)/(Pa/m)
  
  c_por=0.2              ;porosity of rock
  g_por=0.03             ;porosity of grouting loop
  l_por=0.02             ;porosity of liner
  
  den_water=1000

  c_dens=2000            ;denstity of rock
  g_dens=2400            ;denstity of grouting loop
  l_dens=2500            ;denstity of liner

  c_fric=24              ;friction of rock
  g_fric=30              ;friction of grouting loop
  l_fric=35              ;friction of liner

  c_coh=0.72e6           ;cohesion
  g_coh=1.2e6            ;cohesion
  l_coh=1.2e6            ;cohesion

  c_dil=0                ;dilation
  g_dil=0                ;dilation
  l_dil=0                ;dilation  

  c_tens=0.24e6          ;tension
  g_tens=0.8e6           ;tension
  l_tens=1.4e6           ;tension

  c_E=5.9e9
  g_E=6.5e9
  l_E=30e9

  c_v=0.312
  g_v=0.3
  l_v=0.3

  c_bulk=c_E/(3*(1-2*c_v))
  g_bulk=g_E/(3*(1-2*g_v))
  l_bulk=l_E/(3*(1-2*l_v))

  c_shear=c_E/(2*(1+c_v))
  g_shear=g_E/(2*(1+g_v))
  l_shear=l_E/(2*(1+l_v))

  c_k=c_v/(1-c_v)
  g_k=c_v/(1-g_v)
  l_k=c_v/(1-l_v)

  den_ef=c_dens+(c_por-1)*den_water
  str_ver=(den_ef+den_water)*10
  str_hor=(c_k*den_ef+den_water)*10
  str_hor_c=0.8*str_ver  

end
setup

config fluid

; -------------- mechanical model ----------------

model elas

ini dens c_dens
prop bulk c_bulk shear c_shear  

;;;; --- assign boundary conditions ---
fix x         range x    99.9 100.1
fix x         range x  -100.1 -99.9
fix x y z     range z  -150.1 -149.9
fix y         range y  -0.1 0.1
fix y         range y   0.9 1.1

;; --- initial stress state ---
set grav 0 0 -10
ini szz -45e4  grad 0 0  str_ver
ini sxx -45e4  grad 0 0  str_hor_c
ini syy -45e4  grad 0 0  str_hor_c

apply sxx -45e4 grad 0 0  str_hor_c range x   -100.1  -99.9
apply sxx -45e4 grad 0 0  str_hor_c range x    99.9  100.1
apply syy -45e4 grad 0 0  str_hor_c range y   -0.1   0.1
apply syy -45e4 grad 0 0  str_hor_c range y    0.9   1.1
apply szz -45e4 grad 0 0  0         range z   -0.1   0.1
apply szz -45e4 grad 0 0  str_ver   range z   -150.1 -149.9

; --------------fluid flow model----------------
model fl_iso
prop perm c_perm poro c_por ; biot_c 1

ini fmod 0
ini fdensity den_water

ini pp  45e4   grad 0 0 -10000
fix pp         range z -0.1 0.1
fix pp         range z -150.1 -149.9
fix pp         range x -100.1 -99.9
fix pp         range x 99.9 100.1

;; --- initial -----

set mech on
set fluid off     ; 在config 状态下，是on by defaunlt

set large
hist unbal
plot add his 1
;solve               ; check initial equilibrium

plot
plot creat view_pp
plot add cont pp
plot show

sav ini_1.sav

ini xdis 0 ydis 0 zdis 0
ini xvel 0.0 yvel 0.0 zvel 0.0
ini state 0

ini dens g_dens range group 9
ini dens l_dens range group 1 any group 4 any group 6 any group 8 any

model mo
prop bulk c_bulk shear c_shear coh c_coh fri c_fric tens c_tens  range group 10 any group 11 any
prop bulk c_bulk shear c_shear coh c_coh fri c_fric tens c_tens  range group 2 any group 3 any group 5 any group 7 any
prop bulk g_bulk shear g_shear coh g_coh fri g_fric tens g_tens  range group 9
prop bulk l_bulk shear l_shear coh l_coh fri l_fric tens l_tens  range group 1 any group 4 any group 6 any group 8 any

model fl_iso range group 9
prop perm g_perm poro g_por range group 9 ; biot_c 1
model fl_iso range group 1 any group 4 any group 6 any group 8 any
prop perm l_perm poro l_por range group 1 any group 4 any group 6 any group 8 any  ; biot_c 1

ini fmod 2e7
ini fdensity den_water

;; --- excavation -----
model null range group 2 any group 3 any group 5 any group 7 any
model fl_null range group 2 any group 3 any group 5 any group 7 any

ini pp 0 range group 2 any group 3 any group 5 any group 7 any


set mech on
set fluid on

hist id 10 gp pp id 1845      ;注浆圈顶pp
hist id 11 gp pp id 2095      ;注浆圈腰pp
hist id 12 gp pp id 2313      ;注浆圈底pp

hist id 13 gp zdisp id 751    ;衬砌顶zdisp   
hist id 14 gp xdisp id 523    ;衬砌腰xdisp   
hist id 15 gp zdisp id 1401   ;衬砌底zdisp   
hist id 16 gp xdisp id 99     ;衬砌腰xdisp

hist id 17 gp pp id 1335      ;衬砌背顶pp
hist id 18 gp pp id 713       ;衬砌背腰pp
hist id 19 gp pp id 1563      ;衬砌背底pp

hist id 100 fltime

def day
  day=fltime/24/3600
end
day

hist id 101 day


hist id 102 qflac

set mech force  1.5e3  
set mech substep 1000 auto
set fluid substep 20

solve age 3.456e6

save result.sav

cumtsafe

无图无真相！！1:(:(:(

echrome

10仿真币，我还得多挣点。。。

zzw1990

你好啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊