求一个群桩模拟的命令流

googgoodstudy

桩数量非常大的那种，需要用到fish的循环语句。

2226460233

我也想找一个学习学习。

tywan0123

;生成网络
gen zone b size 18 10 30 p0 0 0 0 p1 14 0 0 p2 0 10 0  p3 0 0 30
group clay
model mohr
pro den 1900 bulk 4.17e7  shear 1.92e7 coh 0  fric 30  

fix x range x -0.01 0.01
fix x range x 13.99 14.01
fix y range y -0.01 0.01
fix y range y 9.99 10.01
fix z range z -0.01 0.01

;初始地应力
init szz -5.7e5 grad 0 0 1.9e4
init syy -2.44e5 grad 0 0 0.814e4
init sxx -2.44e5 grad 0 0 0.814e4

set gravity 0 0 -10

plot con szz
plot add axes
solve
save 8p1.sav

;承台
gen zone b size 8 4 2 p0 3 3 31 p1 11 3 31 p2 3 7 31 p3 3 3 33
group ct range x 3 11 y 3 7 z 31 33
model elastic range group ct
pro den 2500 bulk 13.89e9 shear 10.42e9  range group ct


;桩
sel pile id 1 begin 4 4 31  end 4 4 11 nseg 20
sel pile id 1 begin 6 4 31  end 6 4 11 nseg 20
sel pile id 1 begin 8 4 31  end 8 4 11 nseg 20
sel pile id 1 begin 10 4 31 end 10 4 11 nseg 20
sel pile id 1 begin 4 6 31  end 4 6 11 nseg 20
sel pile id 1 begin 6 6 31  end 6 6 11 nseg 20
sel pile id 1 begin 8 6 31  end 8 6 11 nseg 20
sel pile id 1 begin 10 6 31 end 10 6 11 nseg 20

sel pile pro emod 2.5e11 nu 0.2 xcarea 0.19625 &
xcj 0 xciy 0.0625 xciz 0.0625 per 3.14 &
cs_sk 1.3e11 cs_scoh 0 cs_sfric 30.0 &
cs_nk 1.3e11 cs_ncoh 0 cs_nfric 0 cs_ngap off


sel link attach xdir rigid ydir rigid zdir rigid &
xrdir rigid yrdir rigid zrdir rigid range id 1 any id 22 &
any id 43 any id 64 any id  85 any id 106 any id 127 any id 148

;初始置零
ini state 0
ini xdis 0 ydis 0 zdis 0
ini xvel 0 yvel 0 zvel 0


apply szz -1e5 range z 32.9 33.1 group ct

hist id 1 unbal
hist id 2 gp zdis 7 5 33
hist id 3 zone szz 7 5 32.95

plot clear

plot set magn 0.8
plot set plane ori 4 4 0 nor 0 1 0

plot add con zd plane
plot add hist -3 vs -2 xlabel  '承台顶Z位移'  ylabel '承台顶的z分力'
plot show

sel set damp combined
solve

save 8p2.sav

龙一乜

1. def pile7
   
2.    loop n (1,10)
   
3.     x1=3.7+0.5*float(n-1)
   
4.     command
   
5.      sel pile begin (@x1 0.3 23) end (@x1 0.3 26) nseg=6
   
6.      sel pile pro d 1800 emod 2.0e10 nu 0.25 xcarea 0.5027 xciy 20.106e-3 xciz 20.106e-3 xcj
   
7. 
   
8. 40.212e-3 per 2.513 &
   
9.      cs_scoh 18840 cs_sfric 10 cs_sk 50000 cs_ncoh 15000 cs_nfric 10 cs_nk 100000 cs_ng=on slide=on  
   
10.     end_command
    
11.    end_loop
    
12. end
    
13. @pile7

复制代码

仿真cy123

好，谢谢大家的分享啊

757440795

确实不错

guozi001

很好很好

Mortreds

感谢分享