최 등1)은 total lagrangian formulation에 근거한 증분 평형방정식을 적용하고, 강도행렬 산정시 회전각의 2차항을 포함시켜 기하학적 비선형 해석시 해의 수렴성을 향상시켰다. 또한 등매개 쉘 유한요소의 단점인 전단구속 현상과 제로 에너지 모드가 발생하는 문제를 극복하기 위하여 가정 변형률장을 적용하여 보강된 판 및 쉘 구조의 비선형 해석법을 개발하였다. 본 연구에서는 잔류응력을 고려한 쉘구조의 극한강도 해석을 수행하기 위하여, 대변형거동과 함께 소성붕괴거동을 추적할 수 있는 알고리즘을 제시한다. 잔류응력을 고려한 증분평형방정식에 return mapping algorithm을 이용한 탄소성 해석법을 결합시켜서 보강된 판 및 쉘구조의 극한거동을 파악한다. 수치해석 예제를 통하여 본 연구에서 제시된 유한요소 및 비선형 해석 알고리즘에 대한 효율성 및 적용성을 확인하였다.

보강된 판 및 쉘구조의 안정성 및 후좌굴을 포함하는 기하학적 비선형 해석을 수행하기 위하여, total Lagrangian formulation에 근거한 연속체의 증분평형방정식으로부터 변형된 쉘요소인 유한요소이론을 제시하였다. 쉘구조의 곡률이 불연속적으로 변하거나 쉘부재들이 유한한 각도로 만나는 보강된 판 및 쉘구조의 비선형 해석이 가능하도록 주부재와 보강재 간의 연결점에 대한 일반적인 변환관계를 제시하였으며 좌굴해석 및 기하학적 비선형해석의 경우에 해의 정확성 및 수렴성을 개선시키기 위하여 접선강도행렬 산정시 회전각의 2차항을 포함시켰다. 또한, shear locking 현상을 극복하기 위하여 감차적분을 적용하였고 쉘구조의 좌굴해석에서는 power method를 적용하여 해석의 효율을 높였으며, 후좌굴해석에서는 변위 및 하중증분법을 적절히 결합시켜 보강된 쉘구조의 후좌굴 거동추적이 용이하였다. 또한, 입력자료를 손쉽게 준비하고 좌굴모드 및 후좌굴거동을 효율적으로 분석하기 위하여 전, 후 처리 프로그램을 개발하였고 다양한 해석예제를 통하여 다른 문헌의 해석결과를 비교함으로써 본 연구에서 개발된 유한요소 해석프로그램의 타당성 및 정확성을 입증하였다.