직교이방성 적층평판해석을 위해 퇴화 쉘요소에 기초를 둔 p-version 유한요소법이 제안되었다. 이 모델의 비선형 정식화과정에서 기하비선형의 경우 von Karman의 대변형-소변형률 가정을 설명하기 위해 Total Lagrangian 방법이 채택되었으며, 재료비선형의 경우 Huber-Mises의 항복기준과 변형률경화 항복함수에 근거를 둔 Prandtl-Reuss 유동법칙이 사용되었다. 재료모델은 이방성을 표현하는 매개변수에 의해 이방겅재료를 고려할 수 있도록 하였다. 적층평판이론으로는 전단변형 효과를 고려할 수 있는 등가단출이론(ESL Theory)에 기초를 두었기 때문에 두 적층간 계면에서의 전단변형률은 연속이라는 조건을 갖게된다 적분형 르장드르 다항식이 형상함수로 사용되었으며 형상함수의 차수는 1차에서 10차까지 변화시킬 수 있다. 또한, Causs-Lobatto 수치적될법을 사용하기 때문에 기존의 가우스 적분점에서 계산되던 응력값은 이 적분법의 적분점이 절점에 위치하므로 절점에서 바로 응력값이 산출되도록 하였다 극한하중 수렴성, 비선형 효과, 소성역의 형상 등의 비교관점을 통해 p-version 유한요소 모델의 적정성을 보이고자 하였다.

선형탄성 파괴해석은 균열을 갖는 변형률 경화재료의 파괴거동을 예측하는데 불충분하기 때문에 최근에는 균열 선단 부에서 대규모 소성 역을 갖는 균열 체에 적용할 수 있는 많은 파괴역학개념이 제안되고 있다. 따라서, 본 연구에서는 대규모항복 조건하의 연성파괴를 보이는 평판을 정확하게 해석할 수 있는 새로운 유한요소모델을 제시하고자 한다. 균열 선단 부의 응력 장을 정의하는데 가장 지배적인 파괴매개변수인 J-적분 값과 소성 역의 크기 및 형상을 J-적분법과 등가영역적분법을 통해 파괴거동을 설명할 수 있도록 증분소성이론에 기초를 둔 p-version 유한요소해석이 채택되었다. 제안된 유한요소모델에 의한 수치해석결과는 이론 해와 h-version 유한요소해석과 비교되었다.

p-version 유한요소법에 의한 고정밀해석은 응력특이가 발생하는 선형탄성 문제에 매우 적합한 방법으로 인식되고 있다. 해석 결과의 정확도, 모델링의 단순성, 입력자료에 대한 통용성 및 사용자와 CPU 시간의 절감 등 여러장점이 선형탄성 문제에 적용되어 우수성이 입증되었지만, 탄소성 해석분야는 아직 적용이 이루어지지 않고 있다. 그러므로 본 논문에서는 일-경화재료에 대한 구성방정식을 이용하여 정식화된 증분소성이론과 소성유동법칙에 근거한 재료비선형 p-version 유한요소모델이 제안되었다. 비선형방정식을 풀기 위해 Newton-Raphson법과 초기강성도법 등의 반복법이 모색되었다. 제안된 모델을 이용하여 개구부를 가진 사각형 평판과 내압을 받는 두꺼운 실린더, 그리고 등분포하중을 받는 원판해석 등의 수치실험이 수행되었다. 한편, p-version 모델에 의한 해석결과는 문헌의 이론값과 상용유한요소프로그램인 ADINA의 해석결과와 비교 검증되었다.

p-version 유한요소법을 사용한 바닥 슬래브의 탄성해석은 어떤 종류의 요형모서리, 개구 그리고 손상단면을 갖는 점에서 응력특이성을 수반하게 된다. Reissner-Mindlin의 평판이론에 근거한 C.deg.- 평판 계층요소를 사용한 결과가 이론치 및 참고문헌에 발표된 수치해석값과 비교되었다. h-, p-와 hp-version의 수렴속도는 전체적 차원에서의 자유도 증가에 따른 에너지 노름값을 사용하여 예측할 수 있다. 만약에 자유도의 항으로 나타내지는 정확도를 여러 해석방법을 비교하는 기준으로 삼으면 본 연구에서 새로 제안되는 p-version 유한요소해석법의 근사해는 종래의 h-version에 근거하여 현재 까지 발표된 어떤 것 보다 훨씬 효율적 접근방법이라 하겠다. 해석예로는 150.deg. 둔각을 갖는 마름모꼴평판과 손상단면을 갖는 정방형평판이 사용되었다.