무한영역에서 진행하는 탄성파를 유한영역에서 수치적으로 해석하기 위해 많은 흡수경계조건들이 제안되어져 왔다. Paraxial 경계조건은 흡수경계조건의 하나로서 스칼라 및 탄성파 방정식의 paraxial 근사화를 통해 얻어지며, 그 성능이 우수하고 수치해석시 계산적 부담을 주지 않는다. 그러나 경계조건이 복잡한 편미분 방정식으로 표현되어 있어 유한요소해석으로의 적용이 어렵다. 본 논문에서는 penalty function method를 이용하여 전체 에너지 범함수와 paraxial 경계조건을 함께 변분정식화 함으로써 유한요소해석을 수행하였다. 유한요소해석에 가장 적용이 용이하며, 많이 사용되어지는 Lysmer-Kuhlemeyer의 흡수경계조건과 성능을 비교함으로써 연구결과의 타당성을 입증하였다.