이 논문에서는 경계반력법을 이용한 비선형 지반-구조물 상호작용 해석을 위해 LS-DYNA나 MIDAS/Civil 등의 유한요소해석 프로그램과 연계하는 방법을 제시하였다. 경계반력법 적용시 유한요소프로그램에서 구조물과 지반은 선형 또는 비선형 유한요소를 이용하여 모델링하였다. 유한요소의 해석모델 외부의 무한영역으로 전달되는 탄성파를 최대한 흡수하기위해 유한요소 모델의 외측에 LS-DYNA의 경우에는 PML(Perfectly Matched Layer) 요소를, MIDAS/Civil의 경우에는 점성감쇠-스프링 요소를 적용하였다. 비선형 유한요소는 구조물영역에만 적용되는 것으로 가정하였다. 이 연구에서는 입사지진파에 의한 경계반력은 KIESSI-3D 프로그램을 이용하여 계산하였다. 선형 지반-구조물 상호작용 문제에 대해 일반적인 KIESSI-3D의 해석결과와 BRM해석결과를 비교하여 제시된 방법의 효율성을 제시하였다. 또한 수치적 비교를 통해 비선형 구조에 대해 보수적인 응답을 보이는 선형 SSI문제에 대하여 얻은 경계반력이 비선형 지반-구조물 상호작용해석에 효과적으로 적용 가능함을 알 수 있었다.