G03：Itasca（FLAC、UDEC/3DEC、PFC）-ggjj,命令,请教-ggjj,请教我的命令流哪里有问题？，{Simwe仿真论坛(forum.simwe.com)，CAE/CAD/CAM/，FEA/FEM/有限元分析论坛-(手机验证注册)}

handsomejun 发表于 2009-1-14 21:48:01

ggjj,请教我的命令流哪里有问题？

new

gen zone radcylp0 0 0 0 p1 5 0 0 p2 0 20 0 p3 0 0 5 dim 3 3 3 3 size 2 10 4 2ratio 1.0 1.0 1.0 1.2 fill group tunnel
gen zone brick p0 0 0 5 p1 5 0 5 p2 0 20 5 p3 0 0 15 ratio 1.0 1.0 1.0size 2 10 4
gen zone brick p0 5 0 0 p1 25 0 0 p2 5 20 0 p3 5 0 15 ratio 1.0 1.0 1.0size 8 10 6
group soil range group tunnel not
gen zone reflect normal 0 0 -1
; Create named range synonyms
range name=soil group soil
; Assign models to groups
model mohr range soil

; Assign properties
prop shear 1e8 bulk 2e8 fric 30 range soil
ini den 2000

; Settings
set large
set grav 0,0,-10
plot create displ
plot set back white
plot set rot 0 0 40
plot set mag 1
plot add axes red
plot add contour disp out on
plot show

;加位移边界条件
fix x range x 0
fix xrange x 25
fix y range y 0
fixyrange y 20
fixzrange z-15
ini xdis=0 ydis=0 zdis=0                            ;位移归零

;hist fltime                   ;1
;14m横向沉降槽
hist id=2gp zd 0,10,15          ;2
hist id=3gp zd 1,10,15          ;3
hist id=4gp zd 2,10,15          ;4
hist id=5gp zd 3,10,15          ;5
hist id=6gp zd 4,10,15          ;6
hist id=7gp zd 5,10,15          ;7
hist id=8gp zd 6,10,15          ;8
hist id=9gp zd 7,10,15          ;9
hist id=10 gp zd 8,10,15          ;10

;纵向沉降曲线
hist id=11 gp zd 0,1,15             ;11
hist id=12 gp zd 0,2,15             ;12
hist id=13 gp zd 0,3,15             ;13
hist id=14 gp zd 0,4,15             ;14
hist id=15 gp zd 0,5,15             ;15
hist id=16 gp zd 0,6,15             ;16
hist id=17 gp zd 0,7,15             ;17
hist id=18 gp zd 0,8,15             ;18
hist id=19 gp zd 0,9,15             ;19
hist id=20 gp zd 0,10,15          ;20
hist id=21 gp zd 0,11,15          ;21
hist id=22 gp zd 0,12,15          ;22
hist id=23 gp zd 0,13,15          ;23
hist id=24 gp zd 0,14,15          ;24
hist id=25 gp zd 0,15,15          ;25
hist id=26 gp zd 0,16,15          ;26
hist id=27 gp zd 0,17,15          ;27
hist id=28 gp zd 0,18,15          ;28
hist id=29 gp zd 0,19,15          ;29
hist id=30 gp zd 0,20,15          ;30
;孔压
;横向孔压分布
hist id=31 gp pp 3.0,10,0          ;31
hist id=32 gp pp 3.2,10,0          ;32
hist id=33 gp pp 3.4,10,0          ;33
hist id=34 gp pp 3.6,10,0          ;34
hist id=35 gp pp 3.8,10,0          ;35
hist id=36 gp pp 4.0,10,0          ;36
hist id=37 gp pp 4.5,10,0          ;37
hist id=38 gp pp 5.0,10,0          ;38
hist id=39 gp pp 6.0,10,0          ;39
hist id=40 gp pp 7.0,10,0          ;40
hist id=41 gp pp 9.0,10,0          ;41
hist id=42 gp pp 12.0,10,0          ;42
;纵向下部孔压分布
hist id=43 gp pp 0,10,-3.0          ;43
hist id=44 gp pp 0,10,-3.2          ;44
hist id=45 gp pp 0,10,-3.4          ;45
hist id=46 gp pp 0,10,-3.6          ;46
hist id=47 gp pp 0,10,-3.8          ;47
hist id=48 gp pp 0,10,-4.0          ;48
hist id=49 gp pp 0,10,-4.5          ;49
hist id=50 gp pp 0,10,-5.0          ;50
hist id=51 gp pp 0,10,-6.0          ;51
;纵向上部孔压分布
hist id=52 gp pp 0,10,3.0          ;52
hist id=53 gp pp 0,10,3.2          ;53
hist id=54 gp pp 0,10,3.4          ;54
hist id=55 gp pp 0,10,3.6          ;55
hist id=56 gp pp 0,10,3.8          ;56
hist id=57 gp pp 0,10,4.0          ;57
hist id=58 gp pp 0,10,4.5          ;58
hist id=59 gp pp 0,10,5.0          ;59
hist id=60 gp pp 0,10,6.0          ;60
hist id=61 gp zd 0,10,3.0          ;61
hist id=62 gp zd 0,10,3.2          ;62
hist id=63 gp zd 0,10,3.4          ;63
hist id=64 gp zd 0,10,3.6          ;64
hist id=65 gp zd 0,10,3.8          ;65
hist id=66 gp zd 0,10,4.0          ;66
hist id=67 gp zd 0,10,4.5          ;67
hist id=68 gp zd 0,10,5.0          ;68
hist id=69 gp zd 0,10,6.0          ;69


def sup_stress
lumda=0.9
o_press=-208009.55
o_grad=1.39e4
s_press=o_press*lumda
s_grad=o_grad*lumda
end
sup_stress
;开挖第一环
def excate_step1
n=1
a1=2*n-2
a2=2*n
b1=2*n-0.01
b2=2*n+0.01
n1=16*n-1
n2=16*n
t=4*3600*n
end
excate_step1

model null range cyl end1 0 a1 0 end2 0 a2 0 rad 3.0

;添加管片支护
sel shell id=1 range cyl end1 0 a1 0 end2 0 a2 0 rad 3.0 group soil
sel node local xdir 1 0 0 ydir 0 -1 0          range x -0.1 0.1
sel node fix lsys                               range x -0.1 0.1
sel node fix x yr zr                            range x -0.1 0.1
sel node local xdir 0 0 -1 ydir 0 -1 0          range y -0.1 0.1
sel node fix lsys                               range y -0.1 0.1
sel node fix y xr zr                            range y -0.1 0.1
sel shell id=1 prop iso=(3.45e10,0.3) thick=0.35 density 2450
;ret
sel shell apply press -160e3 range cid n1 n2

apply nstress s_press grad 0 0 s_grad range cyl end1 0 b1 0 end2 0 b2 0 rad 3.0
solve

handsomejun 发表于 2009-1-14 21:53:20

附上txt程序代码，谢谢大家，请多指教！

附上程序代码，谢谢大家，请多指教！

jfchai 发表于 2009-1-15 09:13:48

中间要开挖的圆形土芯缺少本构模型、力学参数等等。
所以不能计算。
不过值得称赞，你的模型建的不错。

jfchai 发表于 2009-1-15 09:14:19

中间要开挖的圆形土芯tunnel缺少本构模型、力学参数等等。
所以不能计算。
不过值得称赞，你的模型建的不错。

handsomejun 发表于 2009-1-15 19:14:32

回复 4# jfchai 的帖子

tunnel的土体都是要开挖掉的，所以不设置本构模型，力学参数应该也可以算阿！？

lookcity 发表于 2009-1-16 10:42:48

ini xdis=0 ydis=0 zdis=0 ;位移归零
位移清零太早了，应该是模型开挖前计算平衡了之后在清零

lookcity 发表于 2009-1-16 10:44:32

还有位移边界条件的设置：
fix x range x 0
应该是：
fix x range x -0.1 0.1
其他边界条件类似

shenning139883 发表于 2010-4-17 22:42:41

有限元学得太离散，离散元学得太有限。

shenning139883 发表于 2010-4-17 22:43:28

有限元学得太离散，离散元学得太有限。

shenning139883 发表于 2010-4-17 22:44:01

如何给圆柱形边界加边界应力呢？

aizhongluo 发表于 2010-4-23 16:31:06

你的问题太多了，边界条件不对，模型赋值不对

页: [1]

Simwe仿真论坛(forum.simwe.com)，CAE/CAD/CAM/，FEA/FEM/有限元分析论坛---(手机验证注册)'s Archiver

ggjj,请教我的命令流哪里有问题？

附上txt程序代码，谢谢大家，请多指教！

回复 4# jfchai 的帖子