본 논문에서는 비교적 최근 정립된 co-rotational 이론을 기반으로 한 4절점 평면요소 정식화를 확장하여 9절점 평면 요소에 적합한 CR 정식화를 제시하고자 한다. 그리고 등방성 재료의 소성 해석을 위해, 선형 경과 규칙(bi-linear hardening rule)을 바탕으로 하는 Newton-Raphson return-mapping 알고리즘을 적용하였다. 이때, von Mises 기준을 적용하여 소성 변형 상태를 예측하였다. Lagrange 승수를 도입하여 2차원 접촉에 대한 구속조건을 부여하였다. 개발한 요소는 상용프로그램인 ABAQUS 해석결과와 비교 검증하였다.