본 논문에서는 CAD 시스템에서 사용하는 NURBS 기저함수를 사용하는 아이소-지오메트릭 해석(Isogeometric analysis) 방법과 기 하학적으로 엄밀한 빔 모델링(geometrically exact beam model)을 활용하여 회전과 병진 운동이 결합된 새로운 형태의 메타물질 (metamaterial)에 대한 해석을 진행하였다. 이차원 셀 구조는 자유형상변환(Free-form deformation) 법과 적절한 내삽법(Interpolation) 을 통해 원통 위에 입혀졌다. 원통의 치수와 셀 개수가 비틀림 각도에 미치는 영향이 매개변수 연구(parametric study)를 통해 확인되었 다. 비틀림과 병진 운동이 결합된 구조의 메커니즘에 대해 수치 예제를 통해 알아보았다.

본 논문에서는 정전기 흡착패드를 구성하는 곡면형 전극의 기하학적 엄밀성을 고려하기 위해 정전기 문제에 대하여 CAD에서 사 용하는 NURBS 기저함수를 직접 사용하는 아이소-지오메트릭 해석 기법을 도입하였다. 정전기 흡착력을 곡선 접촉면에서 구하는데 법선 벡터의 영향이 크므로 엄밀한 기하형상을 고려하는 아이소-지오메트릭 해석이 강점을 갖는다. 수치 예제를 통해 곡면과 평면에서 반복 구조의 유무에 따른 파라메터 연구를 수행하여 곡면형 전극의 흡착력이 좋은 성능을 가짐을 보였다. 정전기 흡착력의 성분을 분석하였을 때 정전기 흡착력의 차이는 법선 성분 전기장의 증가로 인한 것으로 파악되었다. 결론적으로 곡면형 전극에서도 전극 사이 거리가 가까워지는 아래로 볼록인 경우가 가장 성능이 좋고, 위로 볼록인 경우에는 성능이 가장 낮음을 보였다.

등기하 해석법을 이용한 고유치 해석은 유한요소를 이용한 결과보다 고차 모드에서 더 정확한 결과를 주는 것으로 알려져 있다. 이는 유한요소법이 차수에 상관없이 요소 간에 C0연속성을 보이는 것과 다르게 등기하 해석법은 p차 요소에 대해서 Cp-1의 연속성을 보장하기 때문이다. 본 논문에서는 이러한 장점을 이용하여 등기하 해석법을 이용하여 모드 기반의 축소 모델을 구성하고 동적 거동 해석을 수행하였다. 축소 모델 구성을 위해 Craig-Bampton(CB) 기법을 적용하였다. 수치 예제를 통해 간단한 봉 요소에 대해 등기하 해석법과 유한요소 해석법을 적용하여 요소의 차수에 따른 고유치 해석 결과를 비교 분석하였다. 등기하 해석법에 중첩 노트를 허용하여 요소 간 연속성을 조절하고, 요소 간 연속성이 줄어듦에 따라 고차 모드에서의 수치 오차가 커짐을 확인하였다. 동적 거동 해석을 위한 축소 모델에 높은 차수의 외력이 주어지는 경우 요소간 연속 성이 높은 등기하해석법을 사용하면, 해의 정확도를 높일 수 있다.

본 연구에서는 NURBS 기반 아이소 지오메트릭 쉘 해석을 위해 다중 패치 해석 모델을 정식화하였다. 기존 연구를 통해 개발된 단일 패치로 구성된 전단 변형을 고려한 쉘 요소에 대해 일반 좌표계에서 기하학적으로 엄밀한 쉘 구조물의 아이소지오메트릭 해석 모델을 도입하고 매개변수 연속성을 고려하여 다중 패치 모델로 확장하였다. 인접 곡면의 노트 요소가 결합 경계를 통해 조화를 이루는 경우에 대해 0차와 1차 매개변수 연속성 조건을 고려하였으며, 두 패치 간 마스터-슬레이브 관계를 정립하여 종속된 한 곡면의 자유도를 상대 곡면의 자유도로 표시하여 모델 크기를 줄이면서 두 곡면을 결합하였다. 다중 패치 쉘 예제에 대해 0차와 1차 연속성 조건을 각각 적용하여 구조해석을 수행하여 1차 연속성 조건의 주요한 특성들을 확인하였다. 또한 각 연속성 조건에 대한 해의 수렴 특성을 검토하였으며 결합 경계에서의 두 패치의 연속성을 확인하였다.

본 논문에서는 2개 이상의 기하형상이 순응 또는 비순응 경계면에서 접합된 멀티패치 문제에 대한 아이소-지오메트릭 해석에 대해서 연구하였다. 패치 경계면에서 응력의 연속성을 표현하는 방법으로 Nitsche 방법론과 마스터-슬레이브 방법에 기반한 방법론에 대해서 지배방정식을 유도하고 아이소-지오메트릭 이산화를 수행하였다. 멀티패치 문제에 대해서 두 방법 론의 차이점을 간단하게 비교하였으며, 후처리 과정에서 사용되는 NURBS 곡면 기반의 응력 복원법에 대해서 기술하였다. 수치예제에서 비순응 경계면을 가지는 멀티패치 빔 문제를 통해 Nitshce 방법론을 검증하였으며, 응력집중을 가지는 문제에서 소개된 두 방법론이 유사한 결과를 보이는 것을 확인하였다. 소개된 NURBS 곡면 기반의 응력 복원법을 후처리에서 도입할 경우 멀티패치 문제의 경계면에서 개선된 연속적인 응력을 보임을 알 수 있다.

A study on the vibration and buckling analyses of laminated composite plates is described in this paper. In order to carry out the analyses of laminated composite plates, a NURBS-based isogeometric general plate element based on Reissner-Mindlin (RM) theory is developed. The non-uniform rational B-spline (NURBS) is used to represent the geometry of plate and the unknown displacement field and therefore, all terms required in this element formulation are consistently derived by using NURBS basis function. Numerical examples are conducted to investigate the accuracy and reliability of the present plate element. From numerical results, the present plate element can produce the isogeometric solutions with sufficient accuracy. Finally, the present isogeometric solutions are provided as future reference solutions.

본 논문에서는 아이소-지오메트릭 기법을 기반으로 민들린 후판에 대한 형상 설계민감도 해석법을 제시하였다. 아이소-지오메트릭 기법은 정확한 기하학적 형상의 표현, 요소 사이의 높은 연속성 등 바람직한 강점들을 가지고 있으며 궁극적으로는 해석해로의 빠른 수렴성과 정확한 설계민감도를 제공한다. 선형 형상함수를 사용하는 유한요소법과는 달리 아이소-지오메트릭 기법에서는 높은 차수의 NURBS 기저함수를 활용하여 CAD 형상의 법선벡터와 곡률을 정확하게 고려한다. 전단 잠김(Shear locking) 현상을 극복하기 위해서 선택적 감소적분(Selective reduced integration) 기법을 사용하였다. 이 간단한 방법은 복잡한 정식화 과정 없이 정확한 아이소-지오메트릭 형상 설계민감도 해석을 수행한다. 굽힘 문제에 대한 수치예제를 통하여 제안된 아이소-지오메트릭 해석과 유한요소 해석을 비교하였으며, 유한차분 설계민감도와 비교하여 아이소-지오메트릭 형상 설계민감도는 매우 정확함을 확인하였다.

Finite element analysis is to approximate a geometry model developed in computer-aided design(CAD) to a finite element model, thus the conventional shape design sensitivity analysis and optimization using the finite element method have some difficulties in the parameterization of geometry. However, isogeometric analysis is to build a geometry model and directly use the functions describing the geometry in analysis. Therefore, the geometric properties can be embedded in the NURBS basis functions and control points so that it has potential capability to overcome the aforementioned difficulties. In this study, the isogeometric structural analysis and shape design sensitivity analysis in the generalized curvilinear coordinate(GCC) systems are discussed for the curved geometry. Representing the higher order geometric information, such as normal, tangent and curvature, yields the isogeometric approach to be the best way for generating exact GCC systems from a given CAD geometry. The developed GCC isogeometric structural analysis and shape design sensitivity analysis are verified to show better accuracy and faster convergency by comparing with the results obtained from the conventional isogeometric method.

본 논문은 비정형의 정확한 형상 설계와 구조해석을 위해 넙스 기저함수를 기반으로 한 변수들로 생성된 평연의 등기하 해석법과 해석 결과에 대한 후처리 방법을 제시한다. 제어점, 매듭값, 차수들로 구성되는 변수들을 인터페이스 기법을 통해 변화시킴으로써 다양한 기하 형상을 구축할 수 있다. 등기하 해석에 사용되는 기저함수는 기하형상 구축에 사용되는 함수와 동일하여 형상의 정확한 설계와 해석이 가능하다. 등기하 해석을 위한 요소망 생성을 위해 h-p-k refinement 과정을 수행함으로써 기존 형상의 변형없이 요소망을 생성하여 구조해석을 수행하였다. 해석에 의한 결과값인 제어점의 변위에 대한 시각화를 위해 IGES 포맷과 넙스기반 3D 설계 프로그램 라이노와의 인터페이스 과정을 수행하여 최종 변형 형상 표현을 위한 후처리 방법을 제시한다.

본 논문에서는 등기하 해석법을 이용하여 선형 탄성문제에 대한 형상 최적설계 기법을 개발하였다. 실용적인 공학문제에 대한 많은 최적설계 문제에서는 초기의 데이터가 CAD 모델로부터 주어지는 경우가 많다. 그러나 대부분의 설계 최적화 도구들은 유한요소법에 기초하고 있기 때문에 설계자는 이에 앞서 CAD 데이터를 유한요소 데이터로 변환해야 한다. 이 변환과정에서 기하 모델의 근사화에 따른 수치적 오류가 발생하게 되고, 이는 응답 해석뿐만 아니라 설계민감도 해석에 있어서도 정확도 문제를 발생시킨다. 이러한 점에서 등기하 해석법은 형상 최적설계에 있어서 유망한 방법론 중 하나가 될 수 있다. 등기하 해석법의 핵심은 해석에 사용되는 기저 함수와 기하 모델을 구성하는 함수가 정확히 일치한다는 것이다. 이러한 기하학적으로 정확한 모델은 설계민감도 해석 및 형상 최적설계에 있어서도 사용된다. 이로 인해 높은 정확도의 설계민감도를 얻을 수 있으며, 이는 설계구배 기반의 최적화에 있어서 매우 중요하게 작용한다. 수치 예제를 통하여 본 논문에서 제시된 등기하 해석 기반의 형상 최적설계 방법론이 타당함을 확인하였다.

본 논문에서는 등기하 해석법을 이용하여 평면 탄성문제의 변분식을 유도하였다. 등기하 해석법은 새로이 부각되고 있는 해석법으로서 기저 함수가 NURBS(Non-Uniform Rational B-Splines) 로부터 직접 생성되므로 해 공간은 CAD 모델을 구성하는 함수로써 표현된다. 또한 CAD 모델의 B-Spline 기저 함수를 직접 사용하므로 기하학적으로 엄밀한 형상을 표현할 수 있고 요소망의 재구성 없이 해석모델을 정밀화(Refinement)할 수 있는 강점이 있다. 본 논문에서는 이를 확장하여 연속체 기반의 애드조인트 설계 민감도 해석법을 사용하는 등기하 설계민감도 해석법을 유도하였다. 기존의 유한요소 기반형상 최적설계는 형상의 매개화에 어려움을 겪었으나 등기하 기반 최적설계에서는 기하학적 정보가 이미 B-spline 기저함수와 조정점에 포함되어 있으므로 이러한 어려움을 피할 수 있는 잠재력을 가지고 있다. 몇몇 수치 예제를 통해서 등기하 해석법을 사용한 설계 민감도 해석을 수행하였으며 유한차분 민감도와 비교하여 정확성을 확인하였다.