본 논문에서는 합성단면을 갖는 구조물의 극한 거동 해석에서 요구하는 재료 및 기하학적 비선형 해석을 수행하기 위한 보 요소를 제시하였다. 제안된 요소는 기하학적 비선형성을 효과적으로 모사할 수 있는 co-rotational 정식화를 통해 도출되 었으며, 다양한 합성단면의 저항성능을 재현할 수 있도록 화이버 단면법이 요소의 내력 및 강성을 산정하는데 활용되었다. 제안된 방법을 구현할 수 있도록 해석 프로그램이 개발되었으며, 호장법을 적용하여 최대내력 발생 이후의 연성거동뿐만 아 니라 심한 비선형 응답(snap-through 또는 snapback)까지 추적해낼 수 있도록 하였다. 본 연구에서 제안된 요소 정식화와 해석 프로그램의 정확성을 검증을 위해 몇 가지 수치예제가 수행되었고, 해석결과는 제안된 요소의 정확성과 효율성을 보이 기 위해 3차원 연속체 모델 및 기존 연구의 결과와 비교되었다. 추가로 합성단면을 갖는 골조 구조물에 대한 수치예제를 통 해, 합성단면을 구성하는 재료의 탄성계수 비 및 강도 비에 따른 영향을 분석하였다. 해석결과는 외층 재료의 탄성계수가 증 가됨에 따라 준취성 거동이 나타났으며, 외층 재료의 항복강도가 높을수록 선형 거동하는 기하적 비선형 응답과 유사한 응 답을 보였다.

플랜트의 유황 저장 탱크는 강재로 구성되며, 탱크 저면은 앵커볼트에 의해 Ring Wall 형상의 콘크리트 기초와 연결된다. 탱크 내 유황이 내부 열원에 의해 고온상태를 유지하기 때문에, 유황 저장 탱크는 상온의 유체를 저장하는 다른 탱크에 비해 큰 체적팽창을 겪게 된다. 일반적으로 탱크 기초의 구조설계는 기초의 내외부의 온도차를 하중으로 적용한 구조해석이 수행 되는데, 이 방법은 탱크의 열팽창 특성이 앵커볼트에 의해 집중하중 형태로 콘크리트 기초에 전달되는 현상을 고려할 수 없 다. 이는 온도하중의 영향을 과소평가하게 되며, 앵커볼트에 인접한 콘크리트의 균열을 야기한다. 본 연구는 앵커볼트에 의 한 온도 하중전달 메커니즘을 고려한 하중 평가식을 제안함으로써, 콘크리트 기초에 작용하는 하중을 보다 합리적으로 결정 하고자 한다. 이를 위해 탱크 바닥판과 앵커볼트가 포함된 유한요소모델을 이용해 앵커볼트 개수 증감에 따른 온도하중의 변화를 분석하였으며, 분석결과를 이론해와 결합해 앵커볼트에 의해 전달되는 하중을 평가할 수 있는 명시적인 형태로 해를 제시하였다. 제안된 식의 유효성을 확인하기 위해 실제 플랜트 현장의 유황 저장 탱크 설계에 적용하였으며, 실무적으로 사 용 가능함을 보였다.

국내외의 실험을 통해 알려진 바와 같이 RC교각의 연성능력은 소성힌지구간에 위치한 주철근 겹침이음 비율 및 횡방향 구속철근이 보유한 횡구속력에 따라 차이를 보이고 있다. 내진설계가 반영되지 못한 기존 교각의 경우 소성힌지구간의 겹침이음 비율에 따라 강도 및 연성능력의 저하에 미치는 영향이 크다. 우리나라에서는1992년 내진설계가 도입된 이후 철근콘크리트 교각의 주철근 겹침 이음에 대한 규정은 없었으나, 2005년 도로교 설계기준에서 주철근 겹침이음을 50% 이내에서 허용하고 있다. 기존 교량의 내진성능 평가 요령에 있어서도 이를 반영하여 교각의 연성능력을 평가하고 있지만 주철근 겹침이음 비율을 단순 유 무에 따라서 구분하여 평가하고 있고 평가기법 또한 명확하게 정립되어 있지 못하다. 따라서 본 연구에서는 비내진 교각의 겹침이음 비율별 연성능력평가를 위하여 현재까지 국내에서 수행된 연구실험결과를 분석하고, 교각의 비선형성 및 주철근의 겹침이음을 고려하기 위하여 섬유요소를 이용한 단면해석으로 RC교각의 연성능력을 산정시 적절한 콘크리트의 극한변형률을 제시하였다.

합리적인 도로 설계를 통해 운전자들의 기대를 만족시켜, 운전자들이 도로의 조건을 충분히 인지하여 주행하도록 하는 것은 매우 중요하다. 이러한 도로의 경우 우리는 흔히 선형의 일관성이 확보되어 운전자들이 안전하고 편안한 차량 주행을 할 수 있는 것이라고 말한다. 위에서 언급한 설계 일관성의 확보 유무를 평가하기 위해 국내 외에서 많은 연구가 이루어지고 있다. 본 연구에서는 기존의 국내 외 설계 일관성 평가 방법을 유형별로 구분하여 정리해 보고, 국내현실에 가장 적합하다고 판단되는 방법론을 응용하여, 일반국도의 과거 조사 자료를 바탕으로 현실적인 설계 일관성기법을 개발하고자 한다. 본 연구에서 개발한 모형은 곡선구간 상류부 100m지점의 속도(Vt)를 변수로 포함하는 예측식으로 양방향 2차로는 R≤200m인 경우와 R〉200m인 경우로 나누어 예측식을 개발하였다. R≤200m인 경우에는 곡선반경(R)과 Vt의 영향이 있고, R〉200m인 구간에서는 기하구조의 영향보다는 곡선 상류부 100m지점의 속도(Vt)가 영향이 있는 것으로 판단된다. 양방향 4차로의 경우도 R값을 400m라는 값에서 두 구간으로 구분하여 예측하였지만 양방향 2차로와는 달리 두 구간에서 모두 기하구조의 영향보다는 Vt의 영향이 더 큰 것으로 판단되었다. 향후 보다신뢰도가 높은 예측모형식 개발을 위해서는 더욱 정확하고 다양한 기하구조의 자료의 수집이 필요하다. 그리고 곡선 상류부 직선구간의 여러 지점의 속도를 수집하면 보다 신뢰도가 높은 예측 모형식 개발이 가능할 것이다.

유체자유수면의 동적거동을 합리적으로 예측하기 위해서는 비선형 특성을 보이는 자유수면의 동역학적 경계조건을 고려해야할 뿐만 아니라 시간에 따라 변화하는 자유수면의 위치변화에 따른 운동학적 경계조건을 고려하여야 한다. 이러한 문제는 대상구조물이 3차원이 될 경우 더욱 복잡해지므로 3차원 비선형 유체자유수면의 해석은 이론해의 도출이 어려우며 수치해석 방법을 이용하는 것이 효과적이다. 본 연구에서는 수치해석 안정성이 높고 3차원 문제에서도 하나의 변수로 유체거동을 모사할 수 있는 arbitrary Lagrangian-Eulerian approach 를 경계요소에 적용하여 효율적이며 안정적인 유체 대변형 해석기법을 개발하였다. 개발된 기법은 향후 자유수면의 비선형 효과를 고려한 유체-구조물 상호작용 해석에 효과적으로 적용할 수 있을 것으로 판단된다.

본 연구에서는 기 개발된 유한요소-경계요소 조합을 통한 지반구조물 상호작용해석기법을 이용하여 주파수 특성이 다른 여러 가지 지진파를 이용한 수치해석을 통하여 지진파에 따른 거동특성을 분석하였다. 사용한 기본해의 검증을 위하여 적용된 다층 반무한 해를 Estorff 등의 연구결과와 비교하였으며, 기본해를 이용하여 개발된 지반-구조물 상호작용해석기법의 검증을 위하여 자유장해석을 수행하였다. 자유장 해석결과는 1차원 파전달 이론에 의하여 개발된 자유장응답해석 프로그램인 SHAKE의 결과와 비교하여 그 타당성을 검증하였다. 검증된 해석기법을 이용하여 특성이 다른 3종류의 지진파가 적용된 2차원 평면상의 지반-구조물 상호작용해석을 수행함으로써 지진파와 말뚝유무에 따른 지반-구조물 상호작용거동특성을 분석하였다. 해석결과 지진이 작용할 때는 말뚝기초를 사용하는 것이 반드시 유리한 결과를 주지는 않는다는 것을 알았으며, 지반의 비선형성을 고려 할때는 말뚝의 유무가 지진응답결과에 큰 영향을 미치지 않는다는 것을 알 수 있었다.

본 연구에서는 지하구조물의 내진 설계를 위한 수정된 응답변위법을 제안하였다. 먼저 기존 응답변위법의 적용성을 검증하기 위해서 매설심도와 지반조건을 변수로 하여 국내의 특성이 반영된 경우에 대한 해석을 수행하였으며 정밀 동적 해석의 결과와 비교하였다. 최대휨모멘트 및 상대 변위를 비교한 결과, 응답변위법의 정확도에 영향을 크게 미지는 속도응답스펙트럼과 지반반력계수 산정의 필요성이 확인되었다. 따라서 국내의 실정에 맞도록 수정된 속도응답스펙트럼과 새로운 지반반력계수 산정식 제안을 위한 방법을 고찰하였다. 제안된 응답변위법의 타당성 및 적용성을 검증하기 위하여, 지반의 조건 및 구조물의 크기 매설심도, 기반암의 위치를 변화시키면서 해석을 수행하였다. 응답변위법의 해석결과와 정밀 동적 해석법의 결과를 비교ㆍ분석하여 국내 지하구조물의 내진 설계를 위한 응답변위법의 적용 가능성을 검증하였다.

지진 발생시 매립지반과 같은 연약지반에 위치하게 되는 지중매설관로는 액상화-영구지반변형에 의해 심각한 규모의 구조적인 손상이 발생할 수 있는 특징을 가지고 있다. 이 경우 매설관로의 거동특성 해석은 수치모형 및 해석적 관계식의 적용을 통해 주로 수행되는데, 특히 횡방향 지반변형을 받는 지중매설관로에 대한 해석적 관계식의 경우. 횡방향 지반변형의 폭에 따라서 해석적 관계식 자체가 이원화되는 단점을 가지고 있는 동시에 액상지반의 특징인 지반강성의 감소로 야기되는 다양한 지반변형의 형상을 반영하지 못하는 단점을 가지고 있다. 그러므로 본 연구에서는 먼저 기존의 해석적 관계식을 개략적으로 살펴본 후, 유한요소 해석과의 비교를 통해 기존의 해석적 관계식이 가지고 있는 적용적 한계성을 검토하였으며, 전체적인 매설관로의 거동을 케이블과 빔의 조합된 형태로 고려하고 지반변형의 다양한 형상을 대변하는 상호작용 형상계수를 도입함으로써 지반변형의 폭과 무관하게 적용될 수 있는 동시에 다양한 지반변형의 형상을 반영할 수 있는 개선된 형태의 해석적 관계식을 제안하였다. 제안된 해석식의 합리적 적용성과 객관적 타당성을 검증하기 위해 지반변형의 크기와 형상 변화에 따른 수치해석을 수행하고 이를 유한요소 해석결과와 비교하였으며, 지반변형의 크기 및 형상변화에 따른 매설관의 거동특성 변화 및 그 의의를 분석하였다.

Response analysis of buried pipeline subjected to permanent ground deformation(PGD) due to liquefaction is mainly executed by use of numerical analysis or semi-analytical relationship, When applying these methods, so called interfacial pipelineㆍsoil interaction force plays an dominant part. Currently used interaction force is mode up of indispensable mechanical and physical components for the response analysis of buried pipeline. However, it has somewhat limited applicability to the liquefied region since it is based on the experimental results for the non-liquefied region. Therefore, in this study, improved type of pipelineㆍsoil interaction force is proposed based on the existing interaction force and experimental research accomplishments. Above all, proposed interaction force includes various patterns of PGD or spatial distributions of interaction force caused by the decrease of soil stiffness. Through the comparison of numerical results using the proposed and the existing interaction force, relative influences of interaction force on the response of pipeline are evaluated and noticeable considerations in the application of semi-analytical relationship are discussed. Moreover, analyses due to the change of pipe thickness and burial depth are performed.

본 연구에서는 액상화-종방향 영구지반변형에 대한 지중매설관로의 거동특성을 해석하기 위하여 수치해석 알고리즘을 개발하였다. 기존의 연구결과가 간략한 해석식의 제안을 중심으로 진행되어 왔으며 영구지반변형의 형상과 폭에 따라 해석방법이 달라지는 단점을 가지고 있었던 것을 고려한다면, 개발된 수치해석 기법은 다양한 영구지반변형의 형상과 폭을 단일한 알고리즘 내에서 처리할 수 있다는 특징을 가지고 있다. 이를 위해 본 연구에서는 연속관 형태의 지중매설관로와 주변지반을 보요소와 등가지반강성으로 표현되는 탄-소성 지반 스프링을 이용하여 모형화하였으며, 지진발생시 실측된 지반변형에 기초하여 영구지반변형의 형상을 5가지의 대표적인 형태로 이상화하여 고려하였다. 국내 계기지진피해사례의 부족으로 인하여 영구지반변형의 크기와 지반변형의 폭은 기존의 연구결과를 참조하여 설정하였으며, 국내에서 사용되는 일반적인 강관을 대상으로 지반변형의 형상과 크기 및 폭, 매설관로의 관경, 관두께 등을 변화시켜 가면서 다양한 수치해석을 수행하였다. 수치해석 결과, 종방향 영구지반변형에 대한 매설관로의 거동에 미치는 주요 인자들의 영향정도를 평가할 수 있었다.

본 연구에서는 2차원 평면상에서 자유장응답 해석을 위하여 유한요소-경계요소 조합에 의한 수치해석기법을 개발하였다. 전체 계를 외부영역과 내부영역으로 구분하였다. 외부영역은 동적 다층반무한 기본해를 이용한 경계요소로 모형화되고 내부영역은 유한요소로 모형화하여 조합하였다. 다층지반의 외부에서 입사하는 지진에 의한 지진응답해석을 수행하기 위하여 동적기본해를 이용한 자유장응답해석을 수행하였다. 지진응답해석에서는 지반의 전단병형률에 따라 변화하는 비선형특성을 모형화하기 위해 등가선형화기법을 적용하였다. 지진응답해석의 검증에 의하여 해석결과를 상용프로그램의 결과와 비교하였다. 결과적으로 지진응답해석을 효과적으로 수행할 수 있는 수치해석기법을 개발하였고 구조물이 있는 경우로의 확장돠 가능하게 되었다.

본 연구에서는 외부의 동하중에 의한 다층 지반-말뚝 상호작용계의 해석을 위한 동적 유한요소-경계요소 조합 주파수 응답해석 알고리즘을 개발하였다. 전체 상호작용계를 내부영역과 외부영역으로 나누고, 내부영역에 보요소를 도입하여 말뚝을 모형화 하고 평면변형률 요소로 모형화된 지반과 조합하였다. 말뚝머리 절점에 집중질량을 이용하여 상부구조물을 고려하므로써 전체 지반-말뚝 상호작용계의 내부영역을 형상화하였다. 외부영역에 동적 기본해를 이용한 경계요소 해석을 도입하고 유한요소로 구성된 내부영역과 조합하므로써 반무한체에 대한 방사조건을 만족시키고 내부의 복잡한 기하학적 성질과 다양한 물성의 고려가 가능한 수치해석 기법을 개발하였다. 개발된 지반-말뚝-구조물계의 상호작용 해석법에 대한 타당성을 알아보기 위해 다층반무한 지반에 근입되어진 말뚝에 조화하중을 가하여 동적 응답해석을 실시하고 기존의 연구결과 및 실험값과 비교 검증하였다. 또한 상호작용계의 주요 인자들의 변화를 통한 다양한 해석을 수행하므로써 다층 반무한 지반에 근입되어진 말뚝의 동적 특성을 고찰하였다.

지진파 전파로 인한 매설관에 작용하는 지진하중은 지진특성 및 지반조건에 따른 지반변형률로부터 산정되어야 한다. 그러나. 기존에 사용되고 있는 경험적인 방법에 의해 계산된 지반변형률 모형은 지진 및 지반의 지역적 특수성을 고려할 수 없는 문제점을 내포하고 있다. 따라서, 본 연구에서는 이러한 문제점을 개선하기 위하여 지진특성 및 지반조선을 반영할 수 있는 수정된 지반변형률 모형을 제안하고 개발된 모형을 매설관로의 지진해석에 지진하중으로 적용하였다. 여기서, 지반변형률을 예측하기 위한 지진판 전파속도는 지반조건을 고려할 수 있도록 파 에너지분포에 근거한 분산곡선을 제안하여 산정하였다. 이러한 과정을 통해 얻어진 지반변형률 산정방법에 타당성을 파악하기 위해 예측한 지반변형률과 과거 지진으로 실측된 지반변형률을 비교하였다. 타당성이 입증된 지반변형률 모형을 매설관의 하중으로 적용하여 지진해석을 실시하였으며, 계산결과는 범용 유한요소해석을 통한 동해석 및 응답변위법에 의한 결과와 비교하였다. 이를 통해 지반 변형률 모형을 적용한 매설관 지진해석의 타당성을 검증하였다. 또한, 지진 및 지반환경이 다른 다양한 관의 특성을 반영하기 위해, 지진 지반 및 관의 영향 인자에 대해 매개변수 해석에 실시되었으며, 이로써 본 연구의 활용성을 검토하였다.

지하구조물의 건전성을 평가하기 위한 비파괴시험으로써 탄성응력파를 이용한 충격반향탐사법을 수치해석적인 방법을 통하여 수행하였다. 즉, 일면만으로 접근 가능한 터널 면에서의 충격가진과 동적응답의 측정으로 이질면을 포함한 내부의 상태를 예측할 수 있다. 연구의 수행은 탄성거동을 하는 매질 내부에서 전파되는 탄성응력파의 특성을 이해하고, 이를 동적 유한요소해석으로 모형화하여 충격반향탐사법을 수치해석적으로 수행한다. 이질재료가 2개의 층을 이루고 있는 경우 표면층의 두께를 쉽게 측정할 수 있었으며, 구조물의 병진운동, 휨운동과 구조물 내에서 다중반사되는 탄성응력파에 의한 복합적인 영향을 받는 터널과 같은 원통형 구조물에서 동적응답의 주파수 특성으로부터 터널라이닝 내부에 형성된 공동의 위치와 크기의 예측이 가능하였다. 수치해석적인 방법과 병행하여 다양한 형태의 경계조건을 가지는 구조물에 대한 충격반향탐사법의 실험을 수행할 경우 실제적인 문제에 적용, 건전성 평가의 지표를 마련할 수 있을 것으로 사료된다.

터널 등과 같은 지하구조계를 유한요소법 등의 수치적 방법으로 해석할 경우 인위적인 경제에서 파의 반사가 발생하게 되어 실제 결과의 큰 차이를 발생시킨다. 따라서 동역학적 하중을 받는 지하구조계는 실질적인 반무한 구조계로 고려되어야 한다. 특히 지하구조계는 실제 다층구조로 구성되어 있으므로 이러한 다층문제를 고려할 수 있어야 한다. 이를 위해 본 연구에서는 외부영역 경제적문제로 해석하기 위한 동적 수치기본해를 개발하였다. 주파수영역의 정적인 경우에 대한 엄밀 적분해와 Bessel 함수의 점근식을 이용한 적분을 통해 축대치문제를 2차원 문제로 보다 쉽게 적용할 수 있도록 하였다. 이와 같이 개발된 동적 수치기본해를 경제 적분 방정식에 적용하여 해석한 결과와 기존 해석결과와의 비교를 통해 그 효율성을 입증하였다. 또한 다층지반내 지하구조물에 대해 지반매체의 각 물성 및 공동의 깊이에 따른 민감도분석을 수해하여 지하구조계의 동적 거동특성 파악의 적용성을 다루었다.

본 논문은 철근 콘크리트(RC) 격자보위에 설치된 access floor 에 발생하는 진동 중 수평방향 진동을 제어하기 위해 정밀 스프링 댐퍼를 이용한 실험을 이용한 실험을 통하여 수평진동 방진 시스템을 개발하였고, 대상 구조물의 진동해석을 위한 모델링과 그에 따른 해석을 수행하였다. 설계된 스프링 댐버 시스템의 방진효과를 알아보기 위하여 모형구조물에 댐퍼를 설치하지 않은 경우, 설피시 댐버를 pedestal과 pedestal에 연결하여 설치한 경우 그리고 pedestal과 격자보에 연결하여 설치한 경우로 나누어 실험을 실시하였다. 각각의 경우에 대해 충격 가진 및 외부 진동을 가해 슬래브와 access floor에서의 가속도 응답을 측정하였다. 실험결과 댐퍼를 설치한 경우에 공진 응답은 댐퍼가 없는 경우에 비해 응답크기가 감소하고, 공진 최대치도 부분적으로 저진동수 대역으로 이동하는 경향으로 나타났다. 또한 저진동수 대역에서 보다는 고진동수 대역에서 가속도 성분의 감소가 크게 나타나고 있고, 특히 외부 진동에 대해서는 상당한 효과가 있다는 것을 알 수 있었다. 대상 구조물의 진동해석을 유한요소법에 의해 실시 하였으며 해석을 통해서도 스프링 댐퍼 시스템의 방진효과를 확인할 수 있었고, 실험과 해석의 결과가 잘 일치함을 알 수 있었다. 따라서 본 연구는 정밀 방.제진 시스템 및 미진동 제어 콘크리트 구조물 설계를 위해 유용하게 이용될 수 있으리라 판단된다.

지하터널은 그 경계가 반무한영역에서 설정되고 재료나 형상의 복잡성을 갖고 있기 때문에, 동적하중에 대하여 정확한 거동을 해석하기 위해서는 3차원 동적해석이 필요하다. 이때 일반적인 수치해석기법인 유한요소만을 이용한 방법은 인위적 경계에서의 파의 반사, 입력자료의 방대함 등으로 인하여 효율적이지 못하게 된다. 본 연구는 이러한 점을 고려하여 지하터널에 직접 가해지는 동적하중에 대한 효율적인 해석기법을 개발하는데 그 목적이 있다. 개발된 프로그램에서 지반의 반무한성은 3차원 경계요소로 고려되었으며, 구조물에는 3차원 동적해석을 수행한 결과 기존의 2차원 터널해석에서 고려가 곤란했던 차량의 진행하중으로 인한 반복효과가 합리적으로 반영되는 것으로 분석되었다.

지하구조물의 주위지반은 일반적으로 퇴적층의 형성 또는 지각의 변동에 의해 다층구조를 가지게 되므로, 구조물 및 주위지반의 거동을 정확히 예측하기 위해서는 해석에 다층구조의 영향을 반영해야 한다. 본 연구에서는 다층으로 구성된 지하구조계를 대상으로 하여 구조물과 그 주변에는 비선형 유한요소를 사용하고, 비선형성이 상대적으로 미약한 주변 다층지반에는 선형 경계요소를 사용하여 재료의 비선형성과 비균질성을 고려한 효율적인 조합해석방법을 개발하고자 한다. 반무한영역에 설정되는 다층구조계를 경계요소로 해석할 경우 그 기본해가 제한되어 있으므로, 본 연구에서는 기존의 무한기본해를 이용하는 방법을 사용하였다. 무한기본해를 이용하는 내부영역문제의 경우 각각의 균질한 층을 부영역(subdomain)으로 분할하고 계방정식을 구성한 뒤에 접합면에 대하여 평형조건과 적합조건을 만족시켜 주는 방법을 사용하여 비균질성을 고려한다. 부영역으로 층을 분할한 내부영역문제의 경계요소해석 결과는 선형 유한요소해석 결과와 비교하여 검증하였고, 검증된 경계요소 프로그램을 비선형 유한요소 프로그램과 조합하였다. 조합해석 결과, 굴착부 주변의 응력집 중부에는 비선형 유한요소를 사용하고, 비선형의 영향이 미소한 주변의 다층지반에 대해서는 부영역에 의한 선형 경계요소를 사용하는 조합해석방법이 합리적이고 효율적임을 알 수 있었다.

최근 원자력의 사용이 증가함에 따라 핵폐기물을 효과적으로 처리하는 문제에 관심이 집중되고 있다. 이러한 핵폐기물을 지층내에 저장할 경우 고온의 열에 의해 핵폐기물 구조체에 지대한 영향을 미치므로 지반의 열력학적 거동을 분석할 필요성이 요구된다. 본 연구는 지반내에 처분된 고온의 사용후 핵연료에 의한 열역학적인 응력이 집중되어 비선형 거동이 예상되는 저장구조체 주변에는 비선형 유한요소를 적용하고 선형거동이 예상되는 무한영역에는 선형경계요소를 사용하여, 일반적인 역학적 계와 동일한 방법으로 비선형 유한요소와 경계요소를 조합한 프로그램을 개발하였다. 사용후 핵연료 폐기구조체와 같이 국부적인 비선형거동이 예상되는 구조물에서는 조합방법이 전 영역을 비선형 유한요소로 모형화하여 해석하는 것보다 효율적임을 알 수 있었다. 또한, 지층내 지반에 영향 미치는 주요 지반계수를 변화시킨 경우, 터널경계의 변위에 이러한 계수들이 어떠한 영향을 미치는가를 개발된 방법을 사용하여 검토하였다. 검토결과, 다른 계수들의 변화보다 열팽창계수의 변화가 터널주위의 변위에 상당한 영향을 미침을 알 수 있었다.