Considering the non-linear behavior of structure and soil when evaluating a nuclear power plant's seismic safety under a beyond-design basis earthquake is essential. In order to obtain the nonlinear response of a nuclear power plant structure, a time-domain SSI analysis method that considers the nonlinearity of soil and structure and the nonlinear Soil-Structure Interaction (SSI) effect is necessary. The Boundary Reaction Method (BRM) is a time-domain SSI analysis method. The BRM can be applied effectively with a Perfectly Matched Layer (PML), which is an effective energy absorbing boundary condition. The BRM has a characteristic that the magnitude of the response in far-field soil increases as the boundary interface of the effective seismic load moves outward. In addition, the PML has poor absorption performance of low-frequency waves. For this reason, the accuracy of the low-frequency response may be degraded when analyzing the combination of the BRM and the PML. In this study, the accuracy of the analysis response was improved by adjusting the PML input parameters to improve this problem. The accuracy of the response was evaluated by using the analysis response using KIESSI-3D, a frequency domain SSI analysis program, as a reference solution. As a result of the analysis applying the optimal PML parameter, the average error rate of the acceleration response spectrum for 9 degrees of freedom of the structure was 3.40%, which was highly similar to the reference result. In addition, time-domain nonlinear SSI analysis was performed with the soil's nonlinearity to show this study's applicability. As a result of nonlinear SSI analysis, plastic deformation was concentrated in the soil around the foundation. The analysis results found that the analysis method combining BRM and PML can be effectively applied to the seismic response analysis of nuclear power plant structures.

In this study, a numerical approach based on mid-point integrated finite elements and a viscous boundary is proposed for time-domain wave-propagation analyses in infinite poroelastic media. The proposed approach is accurate, efficient, and easy to implement in time-domain analyses. In the approach, an infinite domain is truncated at some distance. The truncated domain is represented by mid-point integrated finite elements with real element-lengths and a viscous boundary is attached to the end of the domain. Given that the dynamic behaviors of the proposed model can be expressed in terms of mass, damping, and stiffness matrices only, it can be implemented easily in the displacement-based finite-element formulation. No convolutional operations are required for time-domain calculations because the coefficient matrices are constant. The proposed numerical approach is applied to typical wave-propagation and soil-structure interaction problems. The model is verified to produce accurate and stable results. It is demonstrated that the numerical approach can be applied successfully to nonlinear soil-structure interaction problems.

이 연구에서는 3축 방향 지반운동이 작용하는 지반-구조물 상호작용계의 비선형 지진응답 해석을 수행한다. 비선형 거동 이 예상되는 구조물과 지반의 근역은 비선형 유한요소에 의해 모형을 구성한다. 기하학적 형상과 재료 성질이 균일하고 선 형 거동을 가정하는 원역지반은 무한 영역으로의 에너지 방사를 정확히 고려할 수 있는 3차원 perfectly matched discrete layer에 의해 수치 모형을 구성한다. 이와 같은 지반-구조물 상호작용계의 수치모형을 사용하여 3축 방향 지반운동이 작용 하는 비선형 지진-구조물 상호작용계의 지진응답해석을 수행한다. 3축 방향 지반운동이 작용하는 경우에는 입력 지반운동의 특성에 따라 시스템의 응답이 우세하게 발현되는 방향이 존재하고 그 수준 또한 정밀한 지진응답해석을 통해 산정하여야 한 다. 이 연구의 해석기법은 구조물과 지반의 재료 비선형 거동, 기초와 지반 경계면에서의 경계 비선형 거동 등 다양한 비선 형 지반-구조물 상호작용 해석에 확장 적용할 수 있을 것이다.

이 연구에서는 3축 방향 지반운동이 작용하는 지반-구조물 상호작용계의 비선형 지진응답 해석을 수행한다. 비선형 거동이 예상되는 구조물과 지반의 근역은 비선형 유한요소에 의해 모형을 구성한다. 기하학적 형상과 재료 성질이 균일하고 선형 거동을 가정하는 원역지반은 무한 영역으로의 에너지 방사를 정확히 고려할 수 있는 3차원 perfectly matched discrete layer에 의해 수치 모형을 구성한다. 이와 같은 지반-구조물 상호작용계의 수치모형을 사용하여 3축 방향 지반운동이 작용하는 비선형 지진-구조물 상호작용계의 지진응답해석을 수행한다. 3축 방향 지반운동이 작용하는 경우에는 입력 지반운동의 특성에 따라 시스템의 응답이 우세하게 발현되는 방향이 존재하고 그 수준 또한 정밀한 지진응답해석을 통해 산정하여야 한다. 이 연구의 해석기법은 구조물과 지반의 재료 비선형 거동, 기초와 지반 경계면에서의 경계 비선형 거동 등 다양한 비선형 지반-구조물 상호작용 해석에 확장 적용할 수 있을 것이다.

This paper presents a detailed procedure for a nonlinear soil-structure interaction of a seismically isolated NPP(Nuclear Power Plant) structure using the boundary reaction method (BRM). The BRM offers a two-step method as follows: (1) the calculation of boundary reaction forces in the frequency domain on an interface of linear and nonlinear regions, (2) solving the wave radiation problem subjected to the boundary reaction forces in the time domain. For the purpose of calculating the boundary reaction forces at the base of the isolator, the KIESSI-3D program is employed in this study to solve soil-foundation interaction problem subjected to vertically incident seismic waves. Wave radiation analysis is also employed, in which the nonlinear structure and the linear soil region are modeled by finite elements and energy absorbing elements on the outer model boundary using a general purpose nonlinear FE program. In this study, the MIDAS/Civil program is employed for modeling the wave radiation problem. In order to absorb the outgoing elastic waves to the unbounded soil region, spring and viscous-damper elements are used at the outer FE boundary. The BRM technique utilizing KIESSI-3D and MIDAS/Civil programs is verified using a linear soil-structure analysis problem. Finally the method is applied to nonlinear seismic analysis of a base-isolated NPP structure. The results show that BRM can effectively be applied to nonlinear soil-structure interaction problems.

이 논문에서는 비선형 지반-구조물 상호작용해석을 위한 새로운 시간-주파수영역 복합법을 제시하였다. 제안한 방법은 등가선형 지반-구조물 상호작용해석 프로그램과 범용 비선형 유한요소해석 프로그램을 동시에 사용하는 실용적인 방법이다. 이 방법에서는 먼저 주파수영역에서 등가선형 지반-구조물 상호작용해석을 수행하여 유한요소 영역의 경계면에서 응답을 구한 다음, 이를 범용 비선형 유한요소해석 프로그램에 의한 비선형 동적해석의 시간의존 경계조건으로 입력한다. 제안된 방법의 검증을 위하여 2차원 지하철 정거장 구조물에 대한 지진해석을 수행하였다. 이를 위하여 등가선형 지반-구조물 상호작용해석 프로그램 KIESSI-2D와 비선형 유한요소해석 프로그램 ANSYS를 사용하였다 수치적인 해석결과로부터 이 연구에서 제안한 방법의 타당성을 확인할 수 있었다.

교대-인접토체사이의 상호작용으로 인한 비선형 교대거동이 교량구조물의 전체지진응답에 미치는 영향을 다양한 인자들을 고려할 수 있도록 개발된 이상화된 교량 해석모형을 이용하여 분석하였다. 교대의 비선형 운동은 교대의 강성저하를 반영하는 비선형 스프링으로 모형화하였으며, 비선형효과를 분석하기 위하여 현행 도로교 설계기준에서 제시하고 있는 일정강성을 적용한 선형스프링을 이용한 상대적인 선형모형과 결과를 비교하였다. 분석결과로부터 전체적인 교량구조물의 지진응답은 교대진동계의 모형화 방법 밑 인접한 토체의 조건에 따라 다양하게 나타나며, 교대진동계는 교량구조물의 지진응답에 중요한 영향을 미치는 것으로 분석되었다. 인접진동계간 최대상대거리는 비선형 모델을 적용한 경우가 상당히 증가하는 것으로 나타났으며, 특히 전체 교량구조물에서 낙교의 발생가능성이 가장 큰 위치에서 최대 30%, 정도까지도 증가하는 것으로 분석되었다. 또한 촘촘한 모래를 갖는 토체조건 하에서는 경간수가 증가할수록 교대의 비선형 거동에 따른 영향은 증가하는 것으로 평가되었다. 따라서 교량구조물의 지진거동 분석시 교대의 거동특성을 보다 실제적으로 반영하기 위해서는 교대의 비선형거동이 합리적으로 고려되어야 할 것으로 판단된다.

본 논문에서는 국제공동연구원 대형지진시험구조물의 강세진동시험결과 대한 상관해석와 지진응답해석에 관해 연구하였다. 지반-구조물 상호작용을 위해서 구조물과 근영지반은 유한요소로 모형화하고 원역지반은 무한요소로 모형화하는 직적법을 사용하였으며, 지진응답은 부분구조법에 근거한 파 입력기법을 사용하여 해석하였다. 시험후 상관해석을 통해 각 지반영역의 물성이 강제진동 시험에서 계측된 구조물 응답과 일치하도록 보정하였다. 보정된 지반물성을 초기 선형값으로 사용하고 등가선형화기법을 적용하여 지진에 관한 구조물의 응답을 예측하였다. 지반의 비선형거동을 고려하여 얻어진 구조물 응답은 계측된 결과와 매우 잘 일치한 반면, 초기 선형물성치를 사용한 응답결과는 상당한 차이를 보이고 있어서, 지반 비선형 거동의 영향이 중요함을 알 수 있었다.

파랑제어 및 표사제어라고 하는 해안·항만구조물의 원래의 기능 이외에 해양·해안환경의 유지·개선의 기능까지도 갖춘 투과성잠제에 의한 비선형파랑변형을 해석하였다. 해석법으로서 Stokes 3차파에 기초한 섭동전개법과 경계요소법을 병용하는 주파수영역해석법을 처음으로 투과성잠제의 파랑변형에 적용하였다. 투과층내의 유체운동의 저항에는 Dupuit-Forchheimer의 저항식을 적용하며, Lorentz의 등가일원리에 기초하여 얻어지는 등가선형마찰계수의 산정식을

2차의 섭동법과 경계요소법에 기초한 시간영역해석법은 불규칙파의 파동장에 있어서 파-구조물의 비선형간섭을 해석할 수 있는 해석법이지만. 파와 구조물의 운동이 정상상태에 도달하기까지 시간스텝으로 계산을 수행하여야 하므로 계산시간이 매우 길어지고, 각 성분파와 그에 의한 운동요소를 평가하는 것이 어렵다. 반면에 주파수영역해석법은 계산시간이 짧고, 각 성분요소들의 변화특성을 쉽게 판단할 수 있지만, 불규칙파동장으로의 적용이 현실적으로 어렵다는 단점을 가진다. 본 연구에서는 잠제 등에 대해서 전개되어 있는 주파수영역해석법을 임의형상의 부체 구조물에 대해 새롭게 수식의 전개를 수행하고, 압축공기주입 부체구조물에 적용하여 실험 및 이론해석결과로부터 그의 타당성을 확인한다. 이 때 압축공기의 거동은 Boyle법칙을 사용하여 평가한다.