목차
1.1 UDEC(Universal Distinct Element Code)
1.2 UDEC의 입력 Data
1.3 UDEC 해석 따라하기
1.2 UDEC의 입력 Data
1.3 UDEC 해석 따라하기
본문내용
=xforce2
xtable(yforce2,n)=yforce2
n=n+1
endif
endif
igp=gp_next(igp) ; next gridpoint
endloop
ib= b_next(ib) ; next block in list
endloop
end
relax
STEP 100;STEP 5000
save 39015-2.sav
;-----------------------------------------------------------------;
; ** hard shotcrete **
;-----------------------------------------------------------------;
prop mat 6 st_d 2.5E-3 st_p 0.2 st_ym 15e3
prop mat 6 if_kn 4.46e3 if_ks 2.23e3 cfric 45 if_c 1e0 if_t 1e0
;prop mat 6 st_d=2.5E-3 st_pr=0.2 st_ym=1.5E4 st_yc=21e0 st_y=3.62e-1
;;prop mat 6 if_kn 4.46E7 if_ks 2.23E7 cfric 45 if_c 1E6 if_t 1E6
def remove
ib=block_head ; start of block list
loop while ib #0 ; loop through all blocks
igp=b_gp(ib) ; start of gridpoint list for block ib
loop while igp # 0 ; loop through all gridpoints
ibou=gp_bou(igp) ; index of boundary corner associated with gridpoint
if(ibou) < 0 then ; if address is negative then it is interior
ibou2=abs(ibou)
if (imem(ibou2+2)) = 1 then
imem(ibou2+2)= 1 ; force boundary
imem(ibou2+3)= 1 ; force boundary
ini_force(1,n)=0.0
ini_force(2,n)=0.0
fmem(ibou2+4)= 0.0 ; apply 0% reaction force
fmem(ibou2+5)= 0.0 ; apply 0% reaction force
xforce3=0.0
yforce3=0.0
xforce3=fmem(ibou2+4)
yforce3=fmem(ibou2+5)
ii=out('ibou2 '+string(ibou2))
ii=out('xforce3 '+string(xforce3))
ii=out('yforce3 '+string(yforce3))
xtable(xforce3,ibou2)=fmem(ibou2+4)
xtable(yforce3,ibou2)=fmem(ibou2+5)
endif
endif
igp=gp_next(igp) ; next gridpoint
endloop
ib= b_next(ib) ; next block in list
endloop
end
remove
STEP 100;STEP 5000
save 39015-3.sav
; ** down Excavation **
;del range ann 96.040,60.1172 0,7.6 xr 88.775,103.305 yr 58.999,67.6972
del block 10618,8278,10968,2004,5159
;**************************************************
;하중 분담율 50-30-20
;**************************************************
cycle 1
;하중을 꺼꾸로 주기위한 fish 함수
;터널 벽면의 절점 수에 맞게 ini_force(2,300)의 300를 변경시켜야 한다.
;ini_force(1,300)는 xforce, ini_force(2,300)는 yforce가 저장됨
find
STEP 100;STEP 2000
save 39015-4.sav
; ** soft shotcrete & rockbolt **
struct gen xc 60 yc 28.83 np 100 thi 0.20 mat 6 fang -205.8 the 231
;--------------------------------------------------------
;rockbolt 좌표, 굴진장 1.0m 5.067e-4/굴진장
;--------------------------------------------------------
cable 54.885,29.610 51.250,27.950 10 5 5.067e-4 5 ;connect
cable 54.520,30.751 50.597,29.990 10 5 5.067e-4 5 ;connect
cable 54.440,31.940 50.440,32.140 10 5 5.067e-4 5 ;connect
cable 65.570,31.940 69.560,32.140 10 5 5.067e-4 5 ;connect
cable 65.480,30.751 69.400,29.990 10 5 5.067e-4 5 ;connect
cable 65.110,29.610 68.750,27.950 10 5 5.067e-4 5 ;connect
prop mat 6 st_d 2.5E-3 st_p 0.2 st_ym 5e3
prop mat 6 if_kn 4.46e3 if_ks 2.23e3 cfric 45 if_c 1e0 if_t 1e0
prop mat 5 cb_de 7.85E-3 cb_yi 250e0 cb_yc 25e0 cb_ym 200e3
prop mat 5 cb_kb 1.5e3 cb_sb 0.8e1
relax
STEP 100
;STEP 5000
save 39015-5.sav
; ** hard shotcrete **
prop mat 6 st_d 2.5E-3 st_p 0.2 st_ym 15e3
prop mat 6 if_kn 4.46e3 if_ks 2.23e3 cfric 45 if_c 1e0 if_t 1e0
remove
STEP 100
;STEP 5000
save 39015-6.sav
xtable(yforce2,n)=yforce2
n=n+1
endif
endif
igp=gp_next(igp) ; next gridpoint
endloop
ib= b_next(ib) ; next block in list
endloop
end
relax
STEP 100;STEP 5000
save 39015-2.sav
;-----------------------------------------------------------------;
; ** hard shotcrete **
;-----------------------------------------------------------------;
prop mat 6 st_d 2.5E-3 st_p 0.2 st_ym 15e3
prop mat 6 if_kn 4.46e3 if_ks 2.23e3 cfric 45 if_c 1e0 if_t 1e0
;prop mat 6 st_d=2.5E-3 st_pr=0.2 st_ym=1.5E4 st_yc=21e0 st_y=3.62e-1
;;prop mat 6 if_kn 4.46E7 if_ks 2.23E7 cfric 45 if_c 1E6 if_t 1E6
def remove
ib=block_head ; start of block list
loop while ib #0 ; loop through all blocks
igp=b_gp(ib) ; start of gridpoint list for block ib
loop while igp # 0 ; loop through all gridpoints
ibou=gp_bou(igp) ; index of boundary corner associated with gridpoint
if(ibou) < 0 then ; if address is negative then it is interior
ibou2=abs(ibou)
if (imem(ibou2+2)) = 1 then
imem(ibou2+2)= 1 ; force boundary
imem(ibou2+3)= 1 ; force boundary
ini_force(1,n)=0.0
ini_force(2,n)=0.0
fmem(ibou2+4)= 0.0 ; apply 0% reaction force
fmem(ibou2+5)= 0.0 ; apply 0% reaction force
xforce3=0.0
yforce3=0.0
xforce3=fmem(ibou2+4)
yforce3=fmem(ibou2+5)
ii=out('ibou2 '+string(ibou2))
ii=out('xforce3 '+string(xforce3))
ii=out('yforce3 '+string(yforce3))
xtable(xforce3,ibou2)=fmem(ibou2+4)
xtable(yforce3,ibou2)=fmem(ibou2+5)
endif
endif
igp=gp_next(igp) ; next gridpoint
endloop
ib= b_next(ib) ; next block in list
endloop
end
remove
STEP 100;STEP 5000
save 39015-3.sav
; ** down Excavation **
;del range ann 96.040,60.1172 0,7.6 xr 88.775,103.305 yr 58.999,67.6972
del block 10618,8278,10968,2004,5159
;**************************************************
;하중 분담율 50-30-20
;**************************************************
cycle 1
;하중을 꺼꾸로 주기위한 fish 함수
;터널 벽면의 절점 수에 맞게 ini_force(2,300)의 300를 변경시켜야 한다.
;ini_force(1,300)는 xforce, ini_force(2,300)는 yforce가 저장됨
find
STEP 100;STEP 2000
save 39015-4.sav
; ** soft shotcrete & rockbolt **
struct gen xc 60 yc 28.83 np 100 thi 0.20 mat 6 fang -205.8 the 231
;--------------------------------------------------------
;rockbolt 좌표, 굴진장 1.0m 5.067e-4/굴진장
;--------------------------------------------------------
cable 54.885,29.610 51.250,27.950 10 5 5.067e-4 5 ;connect
cable 54.520,30.751 50.597,29.990 10 5 5.067e-4 5 ;connect
cable 54.440,31.940 50.440,32.140 10 5 5.067e-4 5 ;connect
cable 65.570,31.940 69.560,32.140 10 5 5.067e-4 5 ;connect
cable 65.480,30.751 69.400,29.990 10 5 5.067e-4 5 ;connect
cable 65.110,29.610 68.750,27.950 10 5 5.067e-4 5 ;connect
prop mat 6 st_d 2.5E-3 st_p 0.2 st_ym 5e3
prop mat 6 if_kn 4.46e3 if_ks 2.23e3 cfric 45 if_c 1e0 if_t 1e0
prop mat 5 cb_de 7.85E-3 cb_yi 250e0 cb_yc 25e0 cb_ym 200e3
prop mat 5 cb_kb 1.5e3 cb_sb 0.8e1
relax
STEP 100
;STEP 5000
save 39015-5.sav
; ** hard shotcrete **
prop mat 6 st_d 2.5E-3 st_p 0.2 st_ym 15e3
prop mat 6 if_kn 4.46e3 if_ks 2.23e3 cfric 45 if_c 1e0 if_t 1e0
remove
STEP 100
;STEP 5000
save 39015-6.sav
소개글