안녕하세요?
FLAC 3.4 이하 버젼에서 사면안전율을 구하기 위해서는 전단강도 감소기법을 사용하는 것이 좋을 것 같습니다.
즉, 원지반의 강도정수를 일정한 비율로 감소시켜가면서 사면의 파괴여부를 관찰하다가 사면이 파괴되었을때의 감소된 지반정수와 원지반 강도정수의 비율을 사면의 최소안전율로 채택하는 겁니다.
제 개인적인 생각에는 이 방법은 C-Φ 지반에서 적용하기에는 어려운 점이 있고, 사면의 파괴여부를 판정하는 데에도 어려움이 있는 것이 사실이나 임의의 파괴형상을 고려할 수 있어 기존의 절편법을 보완해줄 수 있는 방법이라 생각됩니다.
위에서 설명한 기본개념을 가지고 아주 간단한 데이터파일을 만들어 첨부하였으니 참고하시기 바랍니다.
(여기서, 사면의 최소안전율은 임의의 안전율 값 범위내에서 탐색하도록 하였으니 사용목적에 맞게 고쳐쓰시면 될 것입니다.)
그럼 좋은 하루 되십시오...
***************************************************************
********** 전단강도 감소기법을 적용한 사면안정해석 예 ***********
NEW
TITLE
Slope Stability Analysis
CONFIG EX=5
*SET LOG ON
GRID 73 73
MO MOHR
GEN -74.0,-20.0 -74.0,54.0 0.0,54.0 0.0,-20.0 I=1,74 J=1,74 RAT 1.0 1.0
def conver
unmax = unbal*0.0001
if unmax < 0.0002
unmax=0.0002
end_if
command
solve force unmax
end_command
end
* Ground Surface modeling
TABLE 1 -74.0,52.3 -56.3,51.3 -54.8,50.6 -50.9,47.4 -49.9,47.4
TABLE 1 -45.9,40.1 -42.9,39.9 -33.0,20.1 -29.9,19.9 -20.0,0.0 0.0,0.0
GEN TABLE 1
* 1st Layer (Soil & Rock)
TABLE 2 -74.0,47.8 -49.3,46.4
TABLE 2 -45.9,40.1 -42.9,39.9 -33.0,20.1 -29.9,19.9 -20.0,0.0 0.0,0.0
GEN TABLE 2
mark i=1 j=72
mark i=1 j=68
ini x -74.0 y 52.5 i 1 j 72
ini x 0.0 y 0.0 i 74 j 21
mar i 74 j 21
m n i 21 j 69,70
m n i 22,23 j 69
m n i 23 j 68,69
m n i 24 j 67,68
m n i 24,25 j 67
m n i 25 j 66,67
m n i 25,26 j 66
m n reg 72,72
;Soil
prop dens=1.9 sh=1.880e4 bu=4.902e4 co=980 fr=33.0 t=980.0 reg 2,70
;Rock
prop dens=2.5 sh=1.008e6 bu=1.603e6 co=980 fr=35.0 t=980.0 reg 2,66
* BOUNDARY CONDITION
fix x i=1
fix x i=74
fix y j=1
* history of unbalance force
his unbal
set grav=9.8
TITLE
Slope Stability Analysis : initial
;
step 1
conver
;
save ini0
*Used k0=1.0
def k0_set
loop i (1,izones)
loop j (1,jzones)
sxx(i,j)=1.0*syy(i,j)
szz(i,j)=1.0*syy(i,j)
end_loop
end_loop
end
k0_set
;
step 1
conver
;
save ini1
; 여기부터는 건기시 해석...
re ini1
ini xdis=0
ini ydis=0
ini xvel=0
ini yvel=0
def prop_ver
loop F (1,5)
FS = 2.1 + (F*0.1)/10.0
coh = coh_org / FS
phi = phi_org / FS
name1='FS = '+ string(FS)
command
title @name1
end_command
command
ini xdis=0 ydis=0 xvel=0 yvel=0
end_command
command
prop dens=1.9 sh=1.880e4 bu=4.902e4 co=coh fr=phi t=coh reg 2,70
end_command
command
step 2000
end_command
command
win -74 -34 24 64
end_command
command
pl hold xv z disp bou
end_command
command
pl hold bou ssr
end_command
command
pl hold hist 1
end_command
end_loop
end
set coh_org=9.8
set phi_org=33.0
prop_ver
save s-dry.sav
;여기부터는 우기시 해석!!!
re ini1
ini xdis=0
ini ydis=0
ini xvel=0
ini yvel=0
def wet_den
loop i (1,izones)
loop j (1,jzones)
if model(i,j) > 1 then
xa=(x(i,j)+x(i+1,j)+x(i+1,j+1)+x(i,j+1))
xc=0.25*xa
ya=(y(i,j)+y(i+1,j)+y(i+1,j+1)+y(i,j+1))
yc=0.25*ya
if yc < table(3,xc) then
density(i,j) = density(i,j) + 1.0
end_if
end_if
end_loop
end_loop
end
def prop_ver
loop F (1,5)
FS = 1.26+ (F*0.1)/10.0
coh = coh_org / FS
phi = phi_org / FS
name1='FS = '+ string(FS)
command
title @name1
end_command
command
ini xdis=0 ydis=0 xvel=0 yvel=0
end_command
command
prop dens=1.9 sh=1.880e4 bu=4.902e4 co=coh fr=phi t=coh reg 2,70
end_command
command
water table 3 den 1.0
end_command
command
TABLE 3 (-74.0,52.3) (-56.3,51.3) (-54.8,50.6) (-50.9,47.4) (-49.9,47.4)
end_command
command
TABLE 3 (-45.9,40.1) (-42.9,39.9) (-33.0,20.1)
end_command
command
TABLE 3 (-29.9,19.9) (-20.0,0.0) (0.0,0.0)
end_command
command
wet_den
end_command
command
plot hold den blo water apply
end_command
command
step 2000
end_command
command
win -74 -34 24 64
end_command
command
pl hold xv z disp bou
end_command
command
pl hold bou ssr
end_command
command
pl hold hist 1
end_command
end_loop
end
set coh_org=9.8
set phi_org=33.0
prop_ver
save s-wet.sav
****************************************************************** 이상 ...끝
┼ 플랙 피시문에 대해...: 조영진(areajin@hanmail.net) ┼
│
│ 안녕하세요^^
│ FLAC/SLOPE(Ver.4.0) 에서 사면의 안전율을 구하는 옵션이 있던데...
│ 3.3 버전을 보유한 관계로 부득이 피시문으로 짜볼 생각임다.
│ ( 논문을 준비중인 학생인데요^^; DAWSON의 논문에 제시되어있는걸로...)
│ 혹 안전율이 어떤식으로 계산이 되는지 아시면 가르쳐주세요.
│
┼ ┼
FLAC 3.4 이하 버젼에서 사면안전율을 구하기 위해서는 전단강도 감소기법을 사용하는 것이 좋을 것 같습니다.
즉, 원지반의 강도정수를 일정한 비율로 감소시켜가면서 사면의 파괴여부를 관찰하다가 사면이 파괴되었을때의 감소된 지반정수와 원지반 강도정수의 비율을 사면의 최소안전율로 채택하는 겁니다.
제 개인적인 생각에는 이 방법은 C-Φ 지반에서 적용하기에는 어려운 점이 있고, 사면의 파괴여부를 판정하는 데에도 어려움이 있는 것이 사실이나 임의의 파괴형상을 고려할 수 있어 기존의 절편법을 보완해줄 수 있는 방법이라 생각됩니다.
위에서 설명한 기본개념을 가지고 아주 간단한 데이터파일을 만들어 첨부하였으니 참고하시기 바랍니다.
(여기서, 사면의 최소안전율은 임의의 안전율 값 범위내에서 탐색하도록 하였으니 사용목적에 맞게 고쳐쓰시면 될 것입니다.)
그럼 좋은 하루 되십시오...
***************************************************************
********** 전단강도 감소기법을 적용한 사면안정해석 예 ***********
NEW
TITLE
Slope Stability Analysis
CONFIG EX=5
*SET LOG ON
GRID 73 73
MO MOHR
GEN -74.0,-20.0 -74.0,54.0 0.0,54.0 0.0,-20.0 I=1,74 J=1,74 RAT 1.0 1.0
def conver
unmax = unbal*0.0001
if unmax < 0.0002
unmax=0.0002
end_if
command
solve force unmax
end_command
end
* Ground Surface modeling
TABLE 1 -74.0,52.3 -56.3,51.3 -54.8,50.6 -50.9,47.4 -49.9,47.4
TABLE 1 -45.9,40.1 -42.9,39.9 -33.0,20.1 -29.9,19.9 -20.0,0.0 0.0,0.0
GEN TABLE 1
* 1st Layer (Soil & Rock)
TABLE 2 -74.0,47.8 -49.3,46.4
TABLE 2 -45.9,40.1 -42.9,39.9 -33.0,20.1 -29.9,19.9 -20.0,0.0 0.0,0.0
GEN TABLE 2
mark i=1 j=72
mark i=1 j=68
ini x -74.0 y 52.5 i 1 j 72
ini x 0.0 y 0.0 i 74 j 21
mar i 74 j 21
m n i 21 j 69,70
m n i 22,23 j 69
m n i 23 j 68,69
m n i 24 j 67,68
m n i 24,25 j 67
m n i 25 j 66,67
m n i 25,26 j 66
m n reg 72,72
;Soil
prop dens=1.9 sh=1.880e4 bu=4.902e4 co=980 fr=33.0 t=980.0 reg 2,70
;Rock
prop dens=2.5 sh=1.008e6 bu=1.603e6 co=980 fr=35.0 t=980.0 reg 2,66
* BOUNDARY CONDITION
fix x i=1
fix x i=74
fix y j=1
* history of unbalance force
his unbal
set grav=9.8
TITLE
Slope Stability Analysis : initial
;
step 1
conver
;
save ini0
*Used k0=1.0
def k0_set
loop i (1,izones)
loop j (1,jzones)
sxx(i,j)=1.0*syy(i,j)
szz(i,j)=1.0*syy(i,j)
end_loop
end_loop
end
k0_set
;
step 1
conver
;
save ini1
; 여기부터는 건기시 해석...
re ini1
ini xdis=0
ini ydis=0
ini xvel=0
ini yvel=0
def prop_ver
loop F (1,5)
FS = 2.1 + (F*0.1)/10.0
coh = coh_org / FS
phi = phi_org / FS
name1='FS = '+ string(FS)
command
title @name1
end_command
command
ini xdis=0 ydis=0 xvel=0 yvel=0
end_command
command
prop dens=1.9 sh=1.880e4 bu=4.902e4 co=coh fr=phi t=coh reg 2,70
end_command
command
step 2000
end_command
command
win -74 -34 24 64
end_command
command
pl hold xv z disp bou
end_command
command
pl hold bou ssr
end_command
command
pl hold hist 1
end_command
end_loop
end
set coh_org=9.8
set phi_org=33.0
prop_ver
save s-dry.sav
;여기부터는 우기시 해석!!!
re ini1
ini xdis=0
ini ydis=0
ini xvel=0
ini yvel=0
def wet_den
loop i (1,izones)
loop j (1,jzones)
if model(i,j) > 1 then
xa=(x(i,j)+x(i+1,j)+x(i+1,j+1)+x(i,j+1))
xc=0.25*xa
ya=(y(i,j)+y(i+1,j)+y(i+1,j+1)+y(i,j+1))
yc=0.25*ya
if yc < table(3,xc) then
density(i,j) = density(i,j) + 1.0
end_if
end_if
end_loop
end_loop
end
def prop_ver
loop F (1,5)
FS = 1.26+ (F*0.1)/10.0
coh = coh_org / FS
phi = phi_org / FS
name1='FS = '+ string(FS)
command
title @name1
end_command
command
ini xdis=0 ydis=0 xvel=0 yvel=0
end_command
command
prop dens=1.9 sh=1.880e4 bu=4.902e4 co=coh fr=phi t=coh reg 2,70
end_command
command
water table 3 den 1.0
end_command
command
TABLE 3 (-74.0,52.3) (-56.3,51.3) (-54.8,50.6) (-50.9,47.4) (-49.9,47.4)
end_command
command
TABLE 3 (-45.9,40.1) (-42.9,39.9) (-33.0,20.1)
end_command
command
TABLE 3 (-29.9,19.9) (-20.0,0.0) (0.0,0.0)
end_command
command
wet_den
end_command
command
plot hold den blo water apply
end_command
command
step 2000
end_command
command
win -74 -34 24 64
end_command
command
pl hold xv z disp bou
end_command
command
pl hold bou ssr
end_command
command
pl hold hist 1
end_command
end_loop
end
set coh_org=9.8
set phi_org=33.0
prop_ver
save s-wet.sav
****************************************************************** 이상 ...끝
┼ 플랙 피시문에 대해...: 조영진(areajin@hanmail.net) ┼
│
│ 안녕하세요^^
│ FLAC/SLOPE(Ver.4.0) 에서 사면의 안전율을 구하는 옵션이 있던데...
│ 3.3 버전을 보유한 관계로 부득이 피시문으로 짜볼 생각임다.
│ ( 논문을 준비중인 학생인데요^^; DAWSON의 논문에 제시되어있는걸로...)
│ 혹 안전율이 어떤식으로 계산이 되는지 아시면 가르쳐주세요.
│
┼ ┼
