• 터널과지하공간
• /
• 제20권1호
• /
• pp.65-72
• /
• 2010

암반발파에서 평활한 파단면을 확보하고 굴착손상영역을 제어하기 위하여 노치장약공을 이용한 균열제어공법이 제안되었다. 본 연구에서는 노치형상과 장약조건이 균열발생 및 성장에 미치는 영향을 살펴보기 위하여 날개형 노치장약공을 갖는 발파모델을 구축하고 동적 파괴과정 해석법을 이용한 암반 파괴과정 해석을 수행하였다. 그 결과, 노치 길이가 증가함에 균열의 성장 길이가 증가하며 파단면의 거칠기가 감소하고 장약공 상하부에 손상균열의 발생이 억제되는 경향을 보였다. 해석결과로부터 노치 길이 및 개구 폭에 따른 응력집중계수의 변화 및 균열발생 양상을 비교 분석하여 균열제어에 미치는 영향인자에 대하여 고찰하였다.

• 화약ㆍ발파
• /
• 제27권2호
• /
• pp.33-41
• /
• 2009

비장약량 또는 장약계수는 어떤 발파에서 파괴대상이 되는 암석의 총 부피 당 또는 총 무게 당 폭약 소비량으로 정의된다. 암석 톤당 또는 입방미터 당 폭약 소비량의 변화는 언제나 암질변화에 대한 좋은 지표가 된다. 광산현장에서는 통상 광석(ore) 톤당 폭약 소비량을 암석에 대한 발파 용이성의 척도로 사용하는 반면, 건설현장에서는 암석 입방미터 당 폭약 소비량을 사용한다. 본 논문에서는 터널발파를 대상으로 하므로 건설현장에서 사용하는 비장약량의 정의를 채택하였다. 지금까지 다양한 터널발파 설계법들이 제안되어 있지만 이런 방법들을 현장에 적용하였을 때 잘 맞지 않는 경우가 많다. 그 이유는 무엇보다 각 나라나 지역별로 암질조건이 서로 상이하기 때문인 것으로 보이며, 이러한 문제는 발파의 설계자나 시공자에게 기술적으로 상당한 부담이 될 수도 있다. 그런데 만일 우리가 주어진 암석에 대한 적정 비장약량을 알고 있다면 발파설계가 매우 쉬워질 수 있을 것이다. 이런 측면에서 본 논문에서는 사전에 결정된 비장약량이 있는 경우 그 비장약량에 맞추어서 터널발파를 설계할 수 있는 알고리즘을 제안하였다. 이 알고리즘은 터널단면 상의 다양한 영역에 서로 다른 구속도를 부여하는 개념을 토대로 하고 있으며, 기존에 알려진 터널발파 설계방법들과 조합하여 사용할 수 있는 경험적인 설계방법이다.

• 화약ㆍ발파
• /
• 제9권1호
• /
• pp.3-13
• /
• 1991

발파에 의한 지반진동의 크기는 화약류의 종류에 따른 화약의 특성, 장약량, 기폭방법, 전새의 상태와 화약의 장전밀도, 자유면의 수, 폭원과 측간의 거리 및 지질조건 등에 따라 다르지만 지질 및 발파조건이 동일한 경우 특히 측점으로부터 발파지점 까지의 거리와 지발당 최대장약량 (W)간에 깊은 함수관계가 있음이 밝혀졌다. 즉 발파진동식은 $V=K{\cdot}(\frac{D}{W^b})^n{\;}{\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots}$ (1) 여기서 V ; 진동속도, cm /sec D ; 폭원으로부터의 거리, m W ; 지발 장약량, kg K ; 발파진동 상수 b ; 장약지수 R ; 감쇠지수 이 발파진동식에서 b=1/2인 경우 즉 $D{\;}/{\;}\sqrt{W}$를 자승근 환산거리(Root scaled distance), $b=\frac{1}{3}$인 경우 즉 $D{\;}/{\;}\sqrt[3]{W}$를 입방근환산거리(Cube root scaled distance)라 한다. 이 장약 및 감쇠지수와 발파진동 상수를 구하기 위하여 임의거리와 장약량에 대한 진동치를 측정, 중회귀분석(Multiple regressional analysis)에 의해 일반식을 유도하고 Root scaling과 Cube root scaling에 대한 회귀선(regression line)을 구하여 회귀선에 대한 적합도가 높은 쪽을 택하여 비교, 검토하였다. 위 (1)식의 양변에 log를 취하여 linear form(직선형)으로 바꾸어 쓰면 (2)式과 같다. log V=A+BlogD+ClogW ----- (2) 여기서, A=log K B=-n C=bn (2)식은 다시 (3)식으로 표시할 수 있다. $Yi=A+BXi_{1}+CXi_{2}+{\varepsilon}i{\;}{\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots}$(3) 여기서, $Xi_{1},{\;}Xi_{2} ;(두 독립변수 logD, logW의 i번째 측정치. Yi ; ($Xi_1,{\;}Xi_2$)에 대한 logV의 측정치 ${\varepsilon}i$ ; error term 이다. (3)식에서 n개의 자료를 (2)식의 회귀평면으로 대표시키기 위해서는 $S={\sum}^n_{i=1}\{Yi-(A+BXi_{1}+CXi_{2})\}\^2$을 최소로하는 A, B, C 값을 구하면 된다. 이 방법을 최소자승법이 라 하며 S를 최소로 하는 A, B, C의 값은 (4)식으로 표시한다. $\frac{{\partial}S}{{\partial}A}=0,{\;}\frac{{\partial}S}{{\partial}B}=0,{\;}\frac{{\partial}S}{{\partial}C}=0{\;}{\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots}$ (4) 위식을 Matrix form으로 간단히 나타내면 식(5)와 같다. [equation omitted] (5) 자료가 많아 계산과정이 복잡해져서 본실험의 정자료들은 전산기를 사용하여 처리하였다. root scaling과 Cube root scaling의 경우 각각 $logV=A+B(logD-\frac{1}{2}W){\;}logV=A+B(logD-\frac{1}{3}W){\;}\}{\;}{\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots}$ (6) 으로 (2)식의 특별한 형태이며 log-log 좌표에서 직선으로 표시되고 이때 A는 절편, B는 기울기를 나타낸다. $\bullet$ 측정치의 검토 본 자료의 특성을 비교, 검토하기 위하여 지금까지 발표된 국내의 몇몇 자료를 보면 다음과 같다. 물론, 장약량, 폭원으로 부터의 거리등이 상이하지만 대체적인 경향성을 추정하는데 참고할수 있을 것이다. 금반 총실측자료는 총 88개이지만 환산거리(5.D)와 진동속도의 크기와의 관계에서 차이를 보이고 있어 편선상 폭원과 측점지점간의 거리에 따라 l00m말만인 A지역과 l00m이상인B지역으로 구분하였다. 한편 A지역의 자료 56개중, 상하로 편차가 큰 19개를 제외한 37개자료와 B지역의 29개중 2개를 낙외한 27개(88개 자료중 거리표시가 안된 12월 1일의 자료3개는 원래부터 제외)의 자료를 computer로 처리하여 얻은 발파진동식은 다음과 같다. $V=41(D{\;}/{\;}\sqrt[3]{W})^{-1.41}{\;}{\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots}$ (7) (-100m)(R=0.69) $V=124(D{\;}/{\;}\sqrt[3]{W})^{-1.66){\;}{\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots}$ (8) (+100m)(R=0.782) 식(7) 및 (8)에서 R은 구한 직선식의 적합도를 나타내는 상관계수로 R=1인때는 모든 측정자료가 하나의 직선상에 표시됨을 의미하며 그 값이 낮을수록 자료가 분산됨을 뜻한다. 본 보고에서는 상관계수가 자승근거리때 보다는 입방근일때가 더 높기 때문에 발파진동식을 입방근($D{\;}/{\;}\sqrt[3]{W}$)으로 표시하였다. 특히 A지역에서는 R=0.69인데 비하여 폭원과 측점지점간의 거리가 l00m 이상으로 A지역보다 멀리 떨어진 B지역에서는 R=0.782로 비교적 높은 값을 보이는 것은 진동성분중 고주파성분의 상당량이 감쇠를 당하기 때문으로 생각된다.