본 연구는 재료비선형 문제를 다루기 위한 비선형 MLS 차분법의 정식화 과정을 제시한다. MLS 차분법은 절점모델을 기반으로 고속 미분근사식을 활용하여 지배 미분방정식을 직접 이산화 하는데, 변수를 변위로 일원화한 Navier 방정식을 사용하여 탄성재료 문제를 다룬 기존의 MLS 차분법은 재료의 구성방정식을 별도로 고려할 수 없다. 본 연구에서는 비선형 재료의 구성방정식을 반영할 수 있는 강정식화를 위해 1차 미분근사를 반복 사용하는 겹미분근사를 고안했다. 응력의 발산으로 표현되는 평형방정식을 그대로 이산화하고 Newton 방법을 적용하여 반복계산을 통해 수렴해를 찾는 비선형 알고리즘을 제시했다. 응력 계산과 내부변수의 갱신은 return mapping 알고리즘을 활용하였고, 알고리즘 접선계수(algorithmic tangent modulus)의 적용을 통해 빠르고 안정적인 반복계산이 가능하도록 하였다. 재생성 시험을 통해 겹미분근사의 정당성을 검증했고, 비선형재료에 대한 인장문제의 해석을 통해 개발된 비선형 MLS 차분 알고리즘의 정확성과 안정성을 확인하였다.

재료적으로 강한 비선형성을 나타내는 FRP 복합적층판을 진동해석을 수행하였다. 적층판의 적층각도와 순서에 따른 선형해석의 결과와 비선형해석의 결과를 비교, 검토하였다. 본 연구는 FRP 복합적층판의 비선형해석시 동적거동 특성을 예측하기 위한 기초적 연구로써 변위를 비교, 분석한다.

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

철근콘크리트 보에 대해서 인장강화효과의 소성힌지길이를 고려한 재료비선형 해석을 수행하였다. 비선형 해석에서 자유도가 많은 대형구조물에 적용시키기에는 많은 제약이 따르는 복잡한 층상해석기법을 사용하는 대신 단면해석을 통해 미리 구성된 모멘트-곡률 관계를 이용하였으며, 유한요소해석에서 사용요소의 크기에 따른 수치해석상의 오차를 줄이기 위해 인장강화효과와 소성힌지길이 개념을 도입하였다. 마지막으로 제안된 해석 알고리즘의 타당성을 검증하기 위하여 해석결과와 실험결과간의 상호 관계를 비교, 분석하였다.

연속보를 대상으로 한 재료비선형 전달행렬법을 제시하였으며, Gauss-Lobatto 적분법을 사용하여 보의 강성행렬로부터 전달행렬을 도출한다. 전달행렬법에서는 유한요소해석법과는 달리 각 절점의 자유도 수에 상관없이 일정한 미지값만을 가지게 되므로 선형해석뿐만 아니라 비선형해석에서도 빠른 연산속도를 보인다. 연속보에 대한 비선형 해석의 적용 예를 통한 비교 결과, 재료비선형 전달행렬법이 변위-모멘트, 변위-하중, 곡률-모멘트의 관계에서 유한요소해석법에 비해 효율적인 것으로 평가되었다.