균열모형을 이용하여 압축 하중하에 있는 암석의 비선형 거동을 예측하는 것은 가능하다. 암석내의 균열의 성장은 암석의 이방성을 가져오며 이러한 이방성의 정도는 탄성계수의 변화로 표현될 수 있다. 본 연구에서는 균열의 성장에 따른 탄성계수의 변화를 이론적인 균열모형을 통해 예측하고, 이를 탄성파 속도 시험을 통해 구한 탄성계수와의 비교를 수행하였다. 또한, 균열모형에 사용되는 초기 암석내 존재하는 균열에 대한 정보는 암석표면의 이미지를 분석하여 구하였다.

암반이나 암석의 역학적 거동을 파악하기 위하여 각종시험 및 수치해석을 수행할 때 대표요소체적(REV)의 크기를 객관적으로 결정할 필요가 있다. 본 연구에서는 영주화강암을 대상으로 암석의 역학적 거동에 영향을 미티는 여러 가지 요소 중 조암광물의 함유율에 주목하여 해석모델의 크기변화에 따른 조암광물의 비율과 탄성제수의 변화를 검토하여 REV를 산정하였다. 조암광물의 분포는 해석모델의 면적이 약 702$ extrm{mm}^2$(900 요소)을 초과하면서부터 일정해 지고 있으며, 탄성계수는 해석모델의 면적이 약 900$ extrm{mm}^2$(1156요소)을 초과하면서부터 일정치를 보여준다.

링 시험펀의 크기와 편심에 따른 응력변화를 FLAC으로 분석하였다. 분석 결과 중공의 크기가 일정한 범위에서만 안정된 값을 얻을 수 있었다. 중공의 크기 오차를 줄이기 위한 상대 공반경은 약 0.3이 적절할 것으로 평가되었다. 한편 중공의 관심 오차가 커지면 응력분포가 불균일하게 되어 정확한 인장강도 산정이 어려운 것으로 사료되었다. 적절한 인장강도를 얻기 위해 중공의 편차는 일정한 정도로 한정하여야 할 것으로 생각된다. 편심으로 인한 오차를 10% 이내로 제한하기 위해서는 중공의 중심은 가압방향에 수평한 방향으로 3%, 수직한 방향으로는 10% 이내가 되어야 하는 것으로 분석되었다.

암석의 파괴인성 측정법으로서 최근에 연구되고 있는 Brazilian 파괴인성 시험을 PFC code를 이용하여 수치해석적으로 시뮬레이션하였다. 이로부터 가압면 형태, 시험편과 모델링 입자의 크기 관계, 가압면 각도 범위, 하중속도 등의 영향 변수에 대한 분석 및 적정 기준에 대해 검 토하였다. 균등한 하중 전달을 위해 가압면이 평면인 Brazilian 시험편을 도입하였는데, 시험편 중앙에서 초기 인장균열을 생성하고 파괴인성 산정에 필요한 안정한 인장균열 진전 시점의 임계하중이 나타나는 하중-변위 곡선을 얻을 수 있는 가압면 각도는 $20^{\circ}C$ -$40^{\circ}C$ 범위로 나타났다. 가압면 각도가 $20^{\circ}C$ 인 경우에, 디스크 시험편의 구성입자 반경이 1 mm 이하이고 하중속도가 0.01 mm/s 이하인 조건에서는 거의 일정한 파괴인성값을 얻을 수가 있었다. 가압면 각도가 $20^{\circ}C$이상이고 하중속도가 0.01 mm/s 이하는 파괴인성 측정에서 중요한 기본 조건인 시험편 중 에서의 인장균열 생성 및 이 인장균열의 안정적인 진전 제어를 위한 조건인 것으로 볼 수가 있다.

암반 절리 조사는 크게 절리 방향 측정과 트레이스 정보 등을 얻기 위한 매핑으로 구분할 수 있다. 본 연구에서는 절리 방향 측정 및 트레이스 매핑을 위한 기존의 수작업을 대체하기 위하여 입체사진측량학과 화상처리 기법 등을 이용한 암반 절리 조사시스템을 개발하였다. 절리 방향 측정을 위해 입체사진측량학을 적용하였으며, 트레이스 길이 및 빈도등의 정보를 얻기 위해 화상 처리 기법을 이용한 두 가지 혼합 접근법(hybrid approach)을 고안하였다. 입체사진측량학을 이용한 절리 방향 측정 기법 및 화상 처리 기법은 접근이 어려운 지역에서의 원거리 측정 및 신속하고 객관적인 매핑이 가능하다는 장점을 가진다.

지하저장공동에 설치되는 콘크리트 플러그는 가능한 모든 하중조건에 대해 안정하게 유지됨은 물론 저장된 유류의 외부 유출과 지하수의 과다 유입을 방지하는 기능을 수행하여야 하므로 지하 유류비축기지의 건설 및 운영 단계에서 매우 중요한 역할을 한다. 플러그가 시공되고 나면 외부와의 접촉이 모두 차단되게 되므로 플러그가 그 기능을 제대로 발휘하기 위해서는 설계단계에서부터 역학적 및 수리적 고려가 이루어져야한다. 따라서 본 연구에서는 2차원 및 3차원 유한 차분프로그램인 FLAC과 FLA $C^{3D}$를 사용하여 플러그에 대한 역학적 및 수리적 해석을 실시하였다. 먼저 측압계수의 변화, 심도의 변화, 플러그의 형상변화, 암반과 콘크리트 접촉면의 상태변화, 굴착손상권의 발생 등 다양한 조건하에서 플러그의 거동 변화에 대한 역학적 해석을 실시하였으며, 2차원 해석조건과 동일한 조건을 가진 3차원 모델을 구성하여 그 결과를 2차원 해석과 비교, 검토하였다. 또한 플러그의 쐐기깊이 변화에 따른 수리해석을 실시하여 그 결과를 비교하였다.다.

방사성 폐기물지하처분이나 열수 및 축열 에너지저장, 지열에너지 개발 등과 같은 대규모 지하공간 프로젝트들이 대두됨에 따라 역학, 수리 및 열적 거동을 동시에 고려한 연구가 필요하게 되었다. 열-수리-역학 상호작용 해석은 열로 교란되고 지하수로 포화된 암반내의 거동을 열 수리, 역학 3가지 지배방정식의 결합을 통해 구현하는 상당히 복잡한 수치 해석 기법 중의 하나이다. 본 연구에서는 기존의 Biot의 압밀이론에 기초한 수식화들을 이용하여 연속체 암반의 열-수리-역학적 상호작용을 모사할 수 있는 유한요소 프로그램을 개발하였다. 개발된 프로그램의 검증을 위해 등온과 비등온 조건하의 일차원 압밀모델에 대한 해석을 실시하여 해석해와 비교하였다. 이타원 압밀에 대한 변수해석을 통하여 포아송비나 수리적 이방성과 같은 인자들이 매체 거동에 미치는 영향을 조사하였다. 앞으로 본 프로그램에 개별체 절리 모델을 통합시켜 보다 일반적인 불연속 암반의 상호작용 거동 해석에 이용할 수 있을 것이다.

불연속 암반의 역학적, 수리적 거동을 현실적으로 해석하기 위한 가장 중요한 요소 중의 하나는 실제 현상을 정확히 모사할 수 있으면서도, 사용하기 쉬운 절리 모델의 개발에 있다. 본 연구에서는 거친 단일 암석 절리의 역학적 거동을 모사하기 위해 기존의 거칢각 손상을 고려한 절리 거동의 탄소성 구성모델을 각각 1차 거칠기와 2차 거칠기의 개념을 사용하여 확장하였다. 제안된 모델을 개별체 유한요소법을 이용하여 수치해석 프로그램을 작성하였다. 거칠기 특성과 손상상수를 달리하여 실시한 변수 해석을 실시한 결과, 제안된 모델은 다양한 형상의 거친 절리의 전단 및 팽창을 수치적으로 잘 근사할 수 있었다. 실험실에서 실시된 단일 절리의 일방향 및 주기전단시험 결과와 수치해석 결과를 비교하였으며, 이를 통해 제안된 모델의 타당성을 검증하였다. 제안된 역학 모델에 충전물 생성계수와 역학 간극과 수리 간극의 관계를 고려한 수리 거동 모델을 도입하였으며, 단일 절리의 전단-수리 시험 결과와의 비교 해석을 실시하였다. 전단에 따른 절리면 사이의 투수계수 중가는 충전물 생성 및 간극 분포의 일정화 경향으로 인해 일정 전단 변위 이후에는 일정해지는 특성을 보였다.

본 논문은 통계적 기법에 의한 균열망 해석 프로그램(NAPSAC)을 사용하여 암반의 수리특성을 평가하기 위해 시도되었다. 암반내 존재하는 절리의 방향성을 고려한 투수계수를 산정하기 위해 마북리 시험터널 주위의 시추자료와 막장 자료로부터 얻어진 균열에 대한 정보로부터 대상지역의 등가 투수계수를 추정하였다. 대상지역의 관찰자료(균열망에 대한 자료, 수리지질학적 자료)로부터 통계적 균열망 해석을 위한 입력자료를 결정하여 해석모델의 신뢰성을 확보하였다. 구현된 모델로부터 모델의 크기를증가함에 따라서 이방성 투수계수 및 투수계수의 변화를계산하였다. 해석결과대상지역의 투수성은 균열군의 방향성에 의해 강한 이방성을 보였다.

지하 채굴에 따른 침하와 수리전도도의 변화 양상에 관하여 조사를 하였다. 채굴 후 암반내 발생하는 변형률과 변형률 의존 수리전도도와의 상호 관계를 나타내는 식을 얻기 위하여 절리 암반을 등가 다공질 매질로 고려하였으며, 단일 절리의 수리전도도는 평행판 이론을 적용하여 결정하였다. 해석 결과에 따르면 채굴 후 암반내 수리전도도의 변화는 침하에 의해 발생한 변형률에 직접적으로 연관이 있는 것으로 나타났다. 채굴적의 폭을 증가시킴에 따라 최대 침하량과 수리전도도의 크기가 증가하였다. 절리 간격은 수리전도도의 변화에 영향을 주는 것으로 나타났다.

고속철도터널 시공전에 폐광된 광산의 채굴공동의 변형거동이 철도터널의 안정성에 미치는 영향을 조사하고자 이지역에서 지질조사, 암반의 공학적 평가 및 평가요소에 대한 다중회귀분석, 수치해석 입력자료의 취득을 위한 많은 조사를 수행하였다. 암반의 공학적 분류결과 절리면에 대한 상태가 RMR 값을 결정하는데 가장 중요한 변수로 나타났고, Q값을 결정하는데는 절리군의 수가 가장 중요한 변수인 것으로 분석되었다. FLAC에 의한 해석 결과 고속철도구간 서측 하부 50 m지점에 위치한 채굴적의 변형거동은 암반의 역학적 특성에 따라 변형양상은 다를 수 있으나 고속철도터널에 침하성 변위를 유발시킬 가능성이 있을 것으로 판단되었다. 그러나 채굴적을 광체와 같은 역학적 특성을 갖는 재료로 충전 시켰을 때 침하성 변위는 나타나지 않는 것으로 나타났다.

터널공학에 있어서 지보시스템에 가해지는 최종 하중을 제어하기 위한 지보의 거동에 관한 많은 연구가실시되었다. 기술적으로 타당한 설계와 안전율이 확보된 경제적인 시공을 위해서는 해석의 신뢰성이 확보되어야 한다. 또한 굴착과 보강의 일련의 시공과정에 대한 역학적인 이해가 필요하며 암반-지보 반응거동에 대한 규명이 이루어 져야 한다. 암반과 지보의 거동에 관한 대부분의 연구는 단순화한 가정에 의한 이론적 해석이 주를 이루고 있다. 또한 터널 주위의 암반 조건에 따른 명확한 기준이 없어 터널 설계시 어려움이 많다. 본 연구에서는 유한차분해석 프로그램인 FLAC을 이용하여 암반조건에 따른 해석을 실시하여 암반-지보 반응곡선을 구하였다. 실제 시공과 유사한 조건을 부여하기 위해 암반등급과 측압계수를 달리하여 해석을 실시하였다. 그 결과 암반조건에 다른 암반-지보 반응곡선의 nomogram을 도출하였으며, 이는 설계 초기에 지보압 및 터널의 허용변위에 대한 효율적인 예측을 실시하는데 있어 유용할 것이다.

절리가 발달한 암반의 거동 특성을 효율적으로 모델링 할 수 있는 수치해석 방법으로서 널리 알려진 개별요소법은 서로 다른 알고리즘에 기초하여 개발된 두 가지 형태의 수치해석 방법으로 발전해 왔다. 본 논문에서는 Plesha 등에 의해 개발된 implicit알고리즘에 기초한 NURBM과 저자에 의해 개발된 explicit 알고리즘에 기초한 CBLOCK을 이용하여 불연속체 암반의 기초적인 거동 특성을 수치해석적으로 고찰하고 모델링상의 문제점들을 비교 분석하였다. CBLOCK 계산 결과에 의하면 절리의 수직 및 접선 강성비는 절리 분포에 따라 불연속체 암반의 안정성을 좌우하는 중요한 지수로서 사용될 수 있는 것으로 나타나고 있으나 NURBM 제산은 상이한 결과를 보였다. 또한 구조물의 안정성과 국부적인 파괴 거동의 해석에서 두 방법이 상이한 결과를 가져올 수 있는 것으로 나타났고 이는 수치해석적 모델링의 제약적 요소로 작용할 수 있음을 보여 주었다.

지하공간 개발에 대한 관심이 높아지면서 국내에서도 안전한 터널시공에 대한 연구가 많이 진행되었다. 본 연구에서는 터널 시공 중 발생 가능한 터널의 붕락 사고의 특성을 조사하여 각종 터널의 붕락사고 유형을 파악하고, 이에 대한 구체적인 원인을 분석하였으며, 사례분석을 통해 터널 붕락의 주요 요인과 개략적인 대책도 함게 제시하였다. 또한, 국내·외 터널보강공법의 기준과 적용현황을 분석하여 세부적이고 합리적인 보강공법의 체계확립을 위한 기초를 마련하고자 하였다. 본 연구에서는 터널 막장의 안전성 등급을 결정하고, 적절한 보강공법을 선정을 위해 퍼지 수량화 이론과퍼지추론시스템을 기반으로 터널정보 데이터베이스를 구축하여 전문가시스템의 모형을 개발하였다. 개발한 시스템의 검증을 위하여 다양한 보강공법을 시공하였던 한강하저터널을 비롯한 국내·외 3곳의 터널 현장 자료를 이용하여 적용한 결과 실제 시행한 보강공범과 근접한 추론결과를 보였다. 정보화시공을 통해 터널 막장기록과 계측자료의 이용을 극대화하고 객관적인 기준의 부재로 인해 일부 전문가의 경험에만 의존하고 있는 국내 보강공법 시공기술을 보다 발전시켜 합리적인 보강공법을 제시하는데 도움이 되고자 하였다.

본 논문에서는 인도네이사 Kalirnantan섬에 위치하는 Pasir탄전의 노천채탄장을 대상으로 기존 지질자료의 분석 및 지질조사, 암반평가, 현지 및 실험실 시험을 통하여 관련 암석들에 대한 물성시험의 수행과 평사투영법과 수치해석에 의한 안정성해석 등을 근거로 채광장내의 전체사면 특히 Roto 북부 및 남북구역의 사면을 적절하게 설계할 수 있는 방안을 재시하고자 하였다. Roto 지역의 주요 불연속면들은 주로 층리로 구성되어 있으며, 층리의 경사는 북쪽의 약 $60^{\circ}$ 에서 남쪽 지역은 $80^{\circ}$ 전후로 남으로 가면서 경사가 급해지며, 많은 사면에서 평면파괴와 전도파괴가 일어나고 있따. Roto 전 지역의 암석물성 시험결과는 이암온 평균 일축압축강도가 9MPa 이고, 사암은 12MPa로 이암보다는 다소 강도가 높다. 현장조사에 의한 암질을 국제암반역학회 분류기준에 따라 분류하면 Roto 북부지역의 암석은 R2~R3 등급에 속할 정도로 대부분 연약하였으며, 석탄층의 강도는 R1~R2 등급에 해당된다. Roto 북부의 4개 지역(동부, 서부,남부 및 북부)과 Roto 남부의 2개 지역(동부 및 서부)의 사면에 대한 상세한 안정성 해석이 수행되었다. 본 연구에서는 최소 사면안전율을 1.2로 결정하여 평사투영법과 한계평형법의 해석을 통해 계산된 30~60$60^{\circ}$의 사면 경사각에 대하여 다시 유한차준법(FLAC)에 의한 안정성 분석을 실시하였다. 최종 사면각으로 Roto 북부지역의 동부사면은 $34^{\circ}$, 서부사면은 $33^{\circ}$, 북부사면은 $32^{\circ}$, 그리고 Roto 남부의 최종사면각은 동부 및 남부사면 모두 $36^{\circ}$로 제안하였다. 본 노천채탄장의 최종사면과 같이 장기적인 유지를 필요로 하는 사면은 사면의 변형거동을 감시하기 위해 보다 체계적인 계측을 필요로 하며, 또한 사면이 많은 비에 노출되기 쉽기 때문에 사면의 풍화와 침식에 대한 대책연구도 수행되어야 한다.

폐갱도를 이용하여 관광공원을 조성할 목적으로 지하에 건설될 예정인 대형 지하광장의 안정성을 분석하기 위하여 FLAC3D를 이용한 3차원 수치해석을 수행하였다. 모델링은 두 단계로 나누어 실시되었다. 첫 단계는 지하광장 자체의 안정성 분석이며, 두 번째 단계는 인근의 채굴적이 미치는 영향분석이다. 첫 번째 단계에서는 모델의 불평형력 수렴 여부, 광장 주변의 소성 또는 인장파괴 영역의 발생상황 그리고 변위 양상으로써 광장의 안정성을 판단하였다. 두 번째 단계에서는 채굴적의 존재 유무에 따른 광장 주변의 응력상태 변화를 분석하였다. 해석결과 대형 지하광장은 역학적으로 아무런 문제가 없으며 채굴적으로 인한 영향도 받지 않는 것으로 나타났다.

일반적으로 암반 구조물의 안정성에 심각한 영향을 미칠 수 있는 것으로 알려진, 암반 내의 rock bridges는 일종의 opening-mode fracture로 간주될 수 있다. 초창기에는 재료의 저항 특성을 고려하여 rock bridge 자체의 거동 특성에 대한 연구가 있었으나, 최근에 와서는 수치 해석적 연구를 통해 파괴역학의 개념을 적용하여 rock bridge의 분포특성이 암반 구조물의 안정성에 미치는 영향 등을 새롭게 연구하고 있다. 본 연구에서는 rock bridges가 수자파쇄 균열의 분포특성에 미치는 영향을 규명하고자, rock bridges의 분포형태에 따른 암반 내 응력 분포 특성과 균열의 bridge현상을 먼저 살펴보았다. 즉, 균열의 길이에 대한 균열의 간격의 비(s/L 비)와 균열의 길이에 대한 균열의 오버랩의 비(d/L 비)를 변화시키면서 응력 분포 특성을 수치 해석적으로 살펴본 결과, 2개의 균열만으로도 여러 개의 균열 분포 특성을 대표할 수 있음을 알았고, s/L 비가 감소할수록 또는 d/L 비가 증가할수록 암반 내 응력 집중은 커짐을 알 수 있었다. 이러한 수치 해석결과를 토대로, 수압파쇄 균열과 rock bridges의 상호 관계를 수치 해석적으로 연구한 결과, rock bridges가 없을 경우 탄성이론에 맞게 발전하던 수압파쇄 균열이, rock bridges의 존재에 의해 수압파쇄 균열의 전파 방향이 왜곡됨은 물론 수압파쇄 균열에 수직으로 작용하는 원거리 응력도 정확히 표현하지 못함을 알았다. 즉, 수압파쇄 균열 선단에서의 mode I 응력집중계수는 s/L 비가 감소할수록 또는 d/L 비가 증가할수륵 커짐을 알 수 있었다. 이러한 연구결과로부터, 균열의 분포특성을 규명하기 위해 노두에서 일반적으로 실시되는 scanline조사에서 얻을 수 있는 s/L비와 d/L비를 통하여 암반 내에 존재하는 응력의 평가 방법이 새로이 제안될 수 있을 것이다. 장기적인 유지를 필요로 하는 사면은 사면의 변형거동을 감시하기 위해 보다 체계적인 계측을 필요로 하며, 또한 사면이 많은 비에 노출되기 쉽기 때문에 사면의 풍화와 침식에 대한 대책연구도 수행되어야 한다. 15ng/$\textrm{cm}^2$로서 90분간 조사로 27ng/$\textrm{cm}^2$량이 생성되었다. 7-DHC은 당초의 123ng/$\textrm{cm}^2$으로부터 계속 감소되어 150분간 조사시 53ng/$\textrm{cm}^2$량까지 내려갔다.하였으며, 그 외의 항목간에는 대동소이하였다.ckarti 와 E. serrulatus가 스파르가눔의 중간숙주가 될 수 있음을 확인하였다. 충란의 배양에서부터 종숙주의 충란 배출까지 약 2개월 정도의 기간이 소요되었고, 우리 나라 자연환경에서는 5일에서 7월에 주로 이 충체의 유충이 발육되고 전파되는 것으로 추측되었다.러 가지 방법들을 적극 적용하여 금후 검토해볼 필요가 있을 것이다.잡은 전혀 삭과가 형성되지 않았다. 이 결과는 종간 교잡종을 자방친으로 하고 그 자방친의 화분친을 사용할 때만 교잡이 이루어지고 있음을 나타내고 있다. 따라서 여교잡을 통한 종간잡종 품종육성 활용방안을 금후 적극 확대 검토해야 할 것이다하였다.함을 보이고 있다.X> , ZnCl$_{3}$$^{-}$같은 이온과 MgCl$^{+}$, MgCl$_{2}$같은 이온종을 형성하기 때문인것 같다. 한편 어떠한 용리액에서던지 NH$_{4}$$^{+}$의 경우 Dv값이 제일 작았다. 바. 본 연구의 목적중의 하나인 인체유해 중금속이온인 Hg(II), Cd(II)등이 NaCl같은 염화물이 함유된 시료용액에 공해이온으로 존재할 경우 흡착에 의한 제거가 가능하다. 한편 이같은 중금속이

대형 구조물 기초로 암반에 근입된 현장타설말뚝의 사용이 현저히 증가하고 있음에도 암반의 공학적 물성과 설계정수와의 관련이 명확히 정립되지 못한실정이다. 이에 본 연구에서는 암반근입말뚝의 거동 특성을 이론적으로 연구하고, 이를 수치적으로 모델링하기 위하여 지반-구조물 상호작용에 관련한 경계면 물성을 기존 연구 결과에 기초하여 합리적으로 산정하였다. 암반의 물성과 경계면 물성간의 관계를 이용하여 암반근입말뚝의 거동을 수치적으로 모사할 수 있으며 , 또한 현장 계측을 통해 얻은 말뚝의 하중전이 양상이 수치적으로 모사된 결과와 잘 일치함을 확인하였다.

본 연구에서는 발파 설계를 쉽게 하기 위하여 RBD(Rock Blasting Design)프로그램을 개발하였다 이 프로그램은 시험 발파 설계, 벤치 발파 설계, 그리고 발파 진동 해석의 기능을 가지며, 지질 특성 값과 발파 상수를 데이터베이스로 제공하고, 여러 상수 값을 재 입력함으로써 보다 빠르고 간편하게 설계할 수 있고 한번 입력한 상수 값을 사용자가 재사용 할 수 있도록 하였다.

발파진동의 전달과 감쇠 특성에는 대상 암반의 특성과 장약량 등이 영향을 미친다. 본 연구에서는 암반의 물성과 장약량의 변화가 발파진동의 전파에 미치는 영향을 파악하기 위해 유한차분 해석프로그램(FLAC)을 이용하여 해석하였다. 대상 지역은 광주시 제2순환도로 공사장이며 설계시 조사된 암반 물성치를 토대로 하였고, 주된 방법으로는 해석의 최대 주요인자를 추정할 수 있는 실험계획법 중 하나인 Taguchi법을 사용하였다. 해당 지반의 각 물성치와 장약량의 변수 레벨을 결정하였고, 이에 의거하여 거리별 발파진동의 전파과정을 고찰하였다. 그 결과 본 연구에서 고려된 암반의 물성치 중에서 발파진동 전달에 가장 영향을 크게 주는 인자는 밀도인 것으로 나타났다.

지반진동을 차단하기 위한 방법으로 방진구가 효과적이지만 대부분의 구조물은 암반 위에 건설되어 있으므로 방진구를 설치하여 진동전파경로를 차단하는 것은 무리가 된다. 그러므로 천공기를 이용하여 방진공을 형성하여 발파진동을 차단하려는 사례가 늘어나고 있는 실정이다. 본 연구에서 2차원 지반모델에 설치된 방진공의 발파진동차단효과를 수치해석적으로 분석하기 위해 이산요소해석프로그램(UDEC)을 이용하여 동적해석을 수행하였으며 이와 비슷한 조건으로 현장실험에 적용하였다.

발파 진동은 보통 지반진동속도로 계측되며, 발파설계나 구조물의 피해 산정에 적절하게 사용된다. 그러나 발파 진동에 의한 인체의 반응이나 가축의 피해 산정 등의 문제에서는 현행 소음진동규제법에서 명시하는 바, 인체의 감각적 감응척도에 충실한 진동레벨을 주로 사용하므로 변환문제가 발생한다. 본 연구에서는 발파 진동 계측기기로 계측한 데이터를 디지털 필터로 처리하여, 진동레벨로 변환하는 프로그램을 개발하고, 이를 단순조화진동과 실제 발파계측결과에 적용하였다. 그 결과, 단순조화진동은 물론이고 국내에서 널리 사용되고 있는 Instantel Inc.의 Blastmate series에서 계측된 진동데이터도 정밀도 높은 진동레벨로 변환이 가능하였다. 비슷한 정밀도로 계측결과를 ASCII file로 송출할 수 있는 기종에서 계측된 자료도 같은 정도로 처리할 수 있을 것으로 판정되었다