• 대한조선학회논문집
• /
• 제30권1호
• /
• pp.145-156
• /
• 1993

대변형, 대회전, 대변형도 문제를 고려한 탄소성-손상 유한요소 정식화 과정을 연구함으로써 구조물의 모든 비선형 거동 및 손상을 합리적으로 예측할 수 있는 수치모형을 개발하였다. 재료의 소성 변형과정에서 발생되는 손상을 합리적으로 고려하기 위하여 연속체 손상역학의 접근방법을 이용하여 구성방정식을 정식화하였으며 Updated Lagrangian 정식화방법, 호장증분법 등의 비선형 강성방정식 해법을 적용하여 2차원 평면문제를 대상으로 하는 탄소성-손상 유한요소해석 프로그램을 구성하였다. 여러가지 예제 계산을 통하여 이 수치모형의 적용성 및 타당성을 검토한 결과 대변형 문제, 손상을 포함하는 재료 비선형문제 공히 합리적인 해석결과를 제시하고 있슴을 확인할 수 있었다.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2011년도 정기 학술대회
• /
• pp.760-763
• /
• 2011

본 논문의 목적은 효율적인 3절점 판 유한요소를 개발하는 것이다. Hellinger-Reissner 범함수에 근거한 혼합정식화(mixed formulation)를 사용한다. 잠김현상을 일으키는 전단변형률장을 독립적으로 분리한 후, MITC(Mixed Interpolation of Tensorial Components)방법을 이용하여 대체전단변형률장(assumed transverse shear strain field)을 구성한다. 추가적으로 회전된 반변기저벡터(contravariant base vector)로 정의된 근사전단변형률장(approximated transverse shear strain field)에 미지수(unknowns)를 도입하여 혼합정식화를 완성시키고 정적응축(static condensation)을 통해 최종 강성행렬을 구성한다. 거짓영에너지모드시험(spurious zero energy mode test), 조각시험(patch test), 등방성시험(isotropic test) 등을 실시하였으며, 4변 완전구속 정사각형 판 구조물과 60도 기울어진 단순지지 판 구조물 등 예제들을 해석하여 MITC3판 유한요소와 수렴성능을 비교하였다.

• 한국전산구조공학회논문집
• /
• 제20권3호
• /
• pp.303-309
• /
• 2007

무한영역에서 진행하는 탄성파를 유한영역에서 수치적으로 해석하기 위해 많은 흡수경계조건들이 제안되어져 왔다. Paraxial 경계조건은 흡수경계조건의 하나로서 스칼라 및 탄성파 방정식의 paraxial 근사화를 통해 얻어지며, 그 성능이 우수하고 수치해석시 계산적 부담을 주지 않는다. 그러나 경계조건이 복잡한 편미분 방정식으로 표현되어 있어 유한요소해석으로의 적용이 어렵다. 본 논문에서는 penalty function method를 이용하여 전체 에너지 범함수와 paraxial 경계조건을 함께 변분정식화 함으로써 유한요소해석을 수행하였다. 유한요소해석에 가장 적용이 용이하며, 많이 사용되어지는 Lysmer-Kuhlemeyer의 흡수경계조건과 성능을 비교함으로써 연구결과의 타당성을 입증하였다.

• 한국전산구조공학회논문집
• /
• 제16권2호
• /
• pp.165-172
• /
• 2003

본 연구에서는 축대칭 쉘구조물의 정동적해석을 효과적으로 수행할 수 있는 새로운 유한요소를 제안하였다. 본 유한요소는 축대칭 쉘의 전단변형률을 고려하였으며, 쉘의 경계에서 기술할 수 있는 변수들만으로 표현되는 효율적인 형태의 수정된 혼합 변분이론에 바탕하여 유한요소정식화를 수행하였다. 또한, 변위장에 대해 무절점 자유도를 추가적으로 도입하여 요소의 수치적 성능을 크게 향상시켰다 계산의 효율성을 위해, 요소정식화의 최종단계에서 정치조건으로부터 응력매개변수들을 제거하고, 동적축약을 통하여 무절점 자유도 성분들 또한 최종적인 유한요소방정식에서 제거되어짐으로써, 일반적인 변위기저 요소와 같은 크기의 유한요소방정식을 얻을 수 있다. 몇 가지 수치예제들에 대한 해석을 통하여, 무절점 자유도와 변위장에 일치하는 적절한 응력매개변수를 가지는 제안된 혼합 축대칭 쉘요소가 정동적해석에서 대단히 정확하고 효율적임을 확인할 수 있었다

• 한국전산구조공학회논문집
• /
• 제28권5호
• /
• pp.485-490
• /
• 2015

구조의 기하학적 비선형해석을 위해 대표적으로 Total Lagrangian, Updated Lagrangian 정식화 기법이 있다. 이러한 고전적인 정식화 과정은 요소의 변형률을 가정하는 방법에 따라 그리고 요소의 절점 수에 따라 추가의 수학적 정식화 과정이 요구된다. 하지만 비교적 최근에 정립된 Co-roational(CR) 이론은 기 존재하는 보, 판, 쉘 요소에 독립적으로 요소 절점자유도에 따라 일정하게 적용하여 대변위, 작은 변형률을 갖는 구조의 기하비선형 해석을 가능케 한다. 본 논문에서는 회전자유도를 갖는 삼각평면요소에 대한 CR 기법을 정식화하였고 동적해석으로 확장하여 이를 상용프로그램과 검증하였다. 해석에 사용한 삼각평면요소는 OPtimal Triangular(OPT) 평면요소이다.

• 한국항공우주학회지
• /
• 제31권5호
• /
• pp.55-62
• /
• 2003

복합재료 구조물의 동적 접촉 문제를 효율적으로 해석하기 위하여 부영역과 공유면에 기반을 둔 변분 정식화 과정을 제안하였다. 벌칙 함수법을 이용하여 접촉면에서의 부등식 구속 조건은 물론, 유한요소 부영역과 공유면의 연결을 위한 등식 적합 조건까지 만족하게 하였다. 이에 따라 구조 형상이 복잡한 경우라도 공유면에서의 절점 연속성을 별도로 고려하지 않고 전체 영역긍 분할한 후, 분할된 부영역별로 독립적인 유한요소로 모델링하여 필요한 수치 연산을 수행할 수 있다. 개발된 컴퓨터 코드를 이용한 수치 해석을 통하여 제안된 정식화에 대한 여러 특성을 고찰하였다.

• 한국전산구조공학회논문집
• /
• 제12권2호
• /
• pp.119-128
• /
• 1999

본 연구는 등방탄성체에 고속충격하중이 작용하는 경우에 대한 유한요소해석에 관한 것으로 대상구조는 여러 가지 모양의 2차원 평면구조를 택하였다. Galerkin 방법을 이용하여 유한요소 정식화하였으며 직접시간적분법에 의해 수치해를 구하였다. 본 해석에서는 균열이 없는 평판으로 수치해와 이론해를 비교하여 수치해의 신뢰성을 확인하였으며, 0°, 30°, 45°경사 균열이 없는 평판에 적용한 3가지 예를 분석하였다. 수치해석 결과는 이론해의 결과와 상호 잘 일치하였다.

• 한국공간구조학회지
• /
• 제3권4호
• /
• pp.11-17
• /
• 2003

• 한국전산구조공학회논문집
• /
• 제13권1호
• /
• pp.37-49
• /
• 2000

본 논문은 구조물의 동역학 및 열탄성 연성문제 해석을 위한 통합된 유한요소법을 개발하는데 초점을 두고있다. 첫째로, 열전도 방정식에 열변위라는 물리량을 도입하여 동역학의 운동 방정식과 유사하도록 유도한 후, 변분법과 일반좌표계를 이용하여 시간영역에서 정식화하였다. 둘째로, 두 방정식에 라플라스 변환을 동시에 도입하고, 공간변수만을 갖는 형상함수와 가중잔여법을 적용하여 유한요소식을 변환영역에서 표현하였다. 연성된 방정식을 문제의 특성에 따라서 분류하였고 정식화 과정을 검증하였다. 또한 수치해석 알고리듬이 갖는 수치 역 변환의 정성적인 경향에 대하여 검토하였다.

• 대한토목학회논문집
• /
• 제26권2A호
• /
• pp.383-390
• /
• 2006

본 연구에서는 급수전개를 이용한 추계론적 유한요소해석법의 개선을 위한 등가몬테카를로 추계장함수를 제안하고 1차 Taylor전개를 이용한 추계론적 유한요소해석법인 가중적분법에 적용하였다. 일반적으로 1차 Taylor전개를 이용하는 수치해석법에서의 응답변화도는 고려하고 있는 추계장의 분산계수에 대하여 선형거동을 보인다. 그러나 몬테카를로 해석의 경우 추계장 분산계수에 대하여 비선형 거동을 나타낸다. 이는 급수전개법의 1차 Taylor전개에 따른 선형특성에 기인한다. 따라서, 가중적분법에서 사용되는 Taylor전개된 변위벡터와 몬테카를로 해석에서의 변위벡터를 비교하고 이들 두 변위벡터 사이에 상호 불일치 하는 점을 고찰하여 몬테카를로 해석에서의 변위벡터와 등가의 변위벡터를 구성하고 이를 가중적분법에 적용하였다. 제안한 등가몬테카를로 추계장은 본래의 추계장 함수에 대한 고차함수로 주어진다. 평면구조에 대한 수치해석을 통하여 제안한 등가몬테카를로 추계장을 이용한 정식화의 타당성을 고찰하였다 새로운 정식화는 기존의 l차 가중적분법을 위한 정식화 과정과 유사하게 수행할 수 있었다.