巷道喷层shell+拱架beam支护，梁单元模拟拱架无作用，link

lwteng2009

无支护情况下，变形量最大372mm，单独加喷层（shell）情况下348mm，全支护（shell+beam）情况下348mm。
已经处理了LINK的连接。
命令流附上，请大家交流讨论下。


new

;1 圆柱体外环绕放射状网格
gen zone radcylinder size 5 4 22 15 ratio 1 1 1 1 p0 0 0 1.5  p1 5.5 0 1.5 p2 0 1.6 1.5 p3 0 0 7 dimension 2.5 2.5 2.5 2.5

;2
gen zone B size 11 4 7 ratio 1 1 1.0  p0 0 0 7  p1 5.5 0 7 p2 0 1.6 7 p3 0 0 10.5

;3
gen zone B size 10 4 11 ratio 1.1 1 1  p0 5.5 0 1.5  p1 20 0 1.5 p2 5.5 1.6 1.5 p3 5.5 0 7

;4
gen zone B size 10 4 7 ratio 1.1 1 1.0  p0 5.5 0 7  p1 20 0 7 p2 5.5 1.6 7 p3 5.5 0 10.5

;5
gen zone B size 11 4 23 ratio 1 1 1.0  p0 0 0 10.5  p1 5.5 0 10.5 p2 0 1.6 10.5 p3 0 0 22

;6
gen zone B size 10 4 23 ratio 1.1 1 1.0  p0 5.5 0 10.5  p1 20 0 10.5 p2 5.5 1.6 10.5 p3 5.5 0 22



;镜像网格
;gen zone reflect normal(-1 0 0) origin = (0 0 0)

gro 围岩


;巷道圆柱
gen zone c size 10 4 22 ratio 1 1 1  p0 0 0 1.5  p1 2.5 0 1.5 p2 0 1.6 1.5 p3 0 0 4

;巷道矩形
;gen zone B size 10 4 6 ratio 1 1 1  p0 0 0 0  p1 2.5 0 0 p2 0 1.6 0 p3 0 0 1.5

;gen zone reflect normal(-1 0 0) origin = (0 0 0)


gro 巷道 ran gro 围岩 not


gro 1煤 ran z 7 10.5  gro 围岩
gro 2粉细互层 ran z 0 7 gro 围岩

gro 4粉细互层 ran z -1.6 0 gro 围岩
gro 3细砂 ran z 10.5 22 gro 围岩
gro -1泥岩 ran z -1.6 -9.6 gro 围岩
gro -2粉细互层 ran z -9.6 -13.6 gro 围岩
gro -3砂质泥岩 ran z -13.6 -18 gro 围岩

gro 5顶板 ran z 4 7 gro 2粉细互层

;gen zone reflect normal(-1 0 0) origin = (0 0 0)
;gen zone reflect normal(0 0 -1) origin = (0 0 1.5)

attach face



;设置接触面

;gen separate 1煤

;int 1 wrap 2粉细互层  1煤
;int 2 wrap 3细砂 1煤

;interface 1 prop kn 7e9 ks 2e9 c 1e5 f 10 tens 5e4
;interface 2 prop kn 3e9 ks 1e9 c 1e5 f 10 tens 5e4
;----------------------------------------------------------------定义材料属性
def constant
;泥岩
es1=20e8
u1=0.27
  k1=es1/(3.0*(1.0-2.0*u1))
  g1=es1/(2.0*(1.0+u1))
fr1=28           ;fr1=27
c1=10e5           ;c1=2e6
t1=8e5
;细砂
es2=24e8         ;es2=48e9
u2=0.28
  k2=es2/(3.0*(1.0-2.0*u2))
  g2=es2/(2.0*(1.0+u2))
fr2=30
c2=15e5          ;c2=33e6
t2=9e5           ;t2=5e6
;粉细互层
es3=25e8
u3=0.25
  k3=es3/(3.0*(1.0-2.0*u3))
  g3=es3/(2.0*(1.0+u3))
fr3=32
c3=10e5          ;c3=30e6
t3=6e5          ;t3=4e6
;细砂岩=底板
;砂质泥岩
es4=15e8
u4=0.3
  k4=es4/(3.0*(1.0-2.0*u4))
  g4=es4/(2.0*(1.0+u4))
fr4=26           ;fr4=20
c4=9e5           ;c4=14e6     
t4=13.5e5
;煤=改好
es5=10e8          ;es5=4e9
u5=0.35           ;u5=0.4
  k5=es5/(3.0*(1.0-2.0*u5))
  g5=es5/(2.0*(1.0+u5))
fr5=20
c5=6e5
t5=8e5

es6=18e8          ;es5=4e9
u6=0.30           ;u5=0.4
  k6=es5/(3.0*(1.0-2.0*u5))
  g6=es5/(2.0*(1.0+u5))
fr6=27
c6=7e5
t6=4e5

end

constant

;定义模型材料参数
mo mo
;pro bulk k1 she g1  fric fr1 coh c1 ten t1 ran &
;gro -1泥岩 any gro 4粉细互层 any gro 3细砂  any gro 2粉细互层 any gro -2粉细互层 any gro 巷道 any gro -3砂质泥岩 any gro  5顶板

pro bulk k1 she g1  fric fr1 coh c1 ten t1

;pro bulk k3 she g3  fric fr3 coh c3 ten t3 ran &
;gro 2粉细互层 any gro -2粉细互层 any gro 巷道 any

;pro bulk k4 she g4  fric fr4 coh c4 ten t4 ran &
;gro -3砂质泥岩

;pro bulk k5 she g5  fric fr5 coh c5 ten t5 ran &
;gro 1煤

;pro bulk k6 she g6  fric fr6 coh c6 ten t6 ran &
;gro  5顶板

;---------------------------------------------------材料属性定义完成
;---------------------------------------------------施加边界条件
fix x ran x -0.1 0.1
fix x ran x 19.9 20.1
fix y ran y 1.59 1.61
fix y ran y -0.1 0.1
fix z ran z 1.49 1.51
;------------------------------------------------------边界条件施加完成

;------------------------------------------------------施加初始应力
ini dens 2500
set grav 0 0 -10
ini szz -2.62e7 grad 0 0 2.5e4
ini syy -3.85e7 grad 0 0 3.85e4
ini sxx -3.4e7 grad 0 0 3.4e4
app szz -2.57e7 ran z 21.9 22.1
;2.57e7=2.62e7-20*2.5e4
;---------------------------------------------------------------初始应力施加完成

PLOT ADD CONTOUR DISP
PLOT ADD AXES BLACK
PLOT SHOW

SOLVE
SAV INI1.SAV




new
restore INI1.SAV


;---------------------------------------------------------------状态归零


ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
ini state 0
set large


mo null ran gro 巷道


sel shell id=1 range cyl end1 0 0 1.5 end2 0 1.6 1.5  rad 2.5 z 1 4

sel shell prop isotropic 25.5e6 0.2 thickness 0.1 density 2500

;sel delete link  ;连接
sel set link node_tol 0.1   ;只有小于该值，两结点之间才能建立link

def parmater   
   w_y=0.40      ;起始位置
   p=2           ;榀数（横向）
   rbeam=2.480    ;半径
   rshell=2.50   ;(洞室半径)
   det=2.5
end
parmater

;beam半径稍小于shell（洞室）半径
def operate_link
       link_id=10000
       w_y=0.40
       loop m(1,p)
           loop n(23,45)
             x_shellnode=-(rshell*cos((n-1)*0.0714))     ;计算rshell的节点坐标;0.0714为弧度
             y_shellnode=w_y         
             z_shellnode=rshell*sin((n-1)*0.0714)+1.50            

             x_beamnode=-(rbeam*cos((n-1)*0.0714))       ;计算钢拱架的节点坐标
             y_beamnode=w_y                                  ;既然可以计算钢拱架的节点坐标，为什么不可以再这里定义呢？
             z_beamnode=rbeam*sin((n-1)*0.0714)+1.50            

             xx_beamnode=-(rbeam*cos(n*0.0714))       ;计算钢拱架的节点坐标
             yy_beamnode=w_y                                 
             zz_beamnode=rbeam*sin(n*0.0714)+1.50
             if n<45 then
             command
             sel beam id 1   beg x_beamnode y_beamnode z_beamnode end xx_beamnode,yy_beamnode,zz_beamnode
             end_command
             endif            

             sn_pointer=nd_near(x_shellnode,y_shellnode,z_shellnode)
             bn_pointer=nd_near(x_beamnode,y_beamnode,z_beamnode)
             shellnode_id=nd_id(sn_pointer)
             beamnode_id=nd_id(bn_pointer)
             command
                  sel link id link_id  beamnode_id target  node  tgt_num shellnode_id
                  sel link attach  xdir=rigid ydir=rigid  zdir=rigid range id link_id
                  ;sel link attach  xrdir=rigid yrdir=rigid  zrdir=rigid range id link_id
                  sel link attach  xrdir=free  yrdir=free zrdir=free range id link_id
             end_command
             link_id=link_id+1
           end_loop
        w_y=w_y+0.8
      end_loop
end
operate_link



sel beam pro emod 2.1e11 nu 0.3
sel beam pro xcarea 45.69e-4 xciz 929e-8 xciy 1245e-8 xcj 0
;sel beam pro xcarea 45.69e-3 xciz 929e-4 xciy 1245e-4 xcj 0

plot show


solve
;sav kaiwa4-wuchenqi.sav

sunforjean

好人啊，顶起来

sunforjean

还在学习flac3d吗？望交流，qq，1182918643

竹轩蓝韵

我现在也在做这个，可以说跟你是一模一样，遇到的问题都一样，都想了两天了

liujie518620

钢拱架的模拟也可以采用shell单元，shell单元纵向长度设置小一点，相当于拱架，模拟效果会好点

xumeizhuowu

好 不错