• Structural Engineering and Mechanics
• /
• 제51권4호
• /
• pp.601-625
• /
• 2014

Based on continuum mechanics and the principle of virtual displacements, incremental total Lagrangian formulation (T.L.) and incremental updated Lagrangian formulation (U.L.) were presented. Both T.L. and U.L. considered the large displacement stiffness matrix, which was modified to be symmetrical matrix. According to the incremental updated Lagrangian formulation, small strain, large displacement, finite rotation of three dimensional Timoshenko fiber beam element tangent stiffness matrix was developed. Considering large displacement and finite rotation, a new type of tangent stiffness matrix of the beam element was developed. According to the basic assumption of plane section, the displacement field of an arbitrary fiber was presented in terms of nodal displacement of centroid of cross-area. In addition, shear deformation effect was taken account. Furthermore, a nonlinear finite element method program has been developed and several examples were tested to demonstrate the accuracy and generality of the three dimensional beam element.

• 한국전산구조공학회논문집
• /
• 제25권3호
• /
• pp.195-202
• /
• 2012

기하학적 비선형성을 고려한 두 개의 비선형 프레임요소의 co-rotational 정식화 과정을 제시한다. 운동학적으로 엄밀한 첫 번째 프레임요소는 변형된 상태의 총 변형성분으로부터 부재력을 산정하며, 정확한 접선강성행렬을 적용한다. 아울러 total Lagrangian 및 updated Lagrangian 정식화에 따른 첫 번째 요소의 엄밀한 접선강성행렬이 동일하다는 것을 보인다. 이에 반하여 두 번째 프레임요소는 절점과 절점사이의 변형을 무시하고 직선으로 가정하여 근사적인 접선강성행렬을 산정하고, 반복계산 시 증분변위로 부터 증분부재력을 구하여 총부재력을 산정한다. 두 개의 수치예제를 통해 첫 번째 프레임 요소가 기하비선형 거동을 추적하는데 있어서 더 정확하고 성능이 우수하다는 것을 입증한다. 특히 케이블부재의 비선형해석 예제를 통하여 첫 번째 프레임 요소가 휨강성을 고려한 케이블요소로 사용할 수 있음을 보인다.

• 대한토목학회논문집
• /
• 제10권3호
• /
• pp.29-38
• /
• 1990

본(本) 연구(硏究)에서는 종래(從來)의 일반적(一般的)인 비선형(非線型) 해석(解析) 방법(方法)인 Lagrangian방법(方法)을 개선(改善)한 새로운 Pseudo Lagrangian방법(方法)에 의한 기하학적(幾何學的) 선형(非線型) 보요소(要素)의 정식화(Formulation) 과정을 제시(提示)하고 예제해석(例題解析)을 통한 대변형(大變形) 해석결과(解析結果)를 고찰(考察)하여 Pseudo Lagrangian방법(方法)에 의한 보요소(要素)의 대변형(大變形) 해석(解析)의 타당성(妥當性)과 정확도(正確度)를 검증(檢證)하였다. Pseudo Lagrangian방법(方法)에서는 변형전(變形前) 또는 변형후(變形後)의 상태(狀態)를 기준(基準)로 변형(變形)을 나타내는 Total Lagrangian방법(方法)과는 달리 어떠한 물리적(物理的)인 의미(意味)가 없는 임의(任意)의 Pseudo(가상(假想))변형상태(變形狀態)를 기준(基準)으로 변형(變形)을 나타낸다. 유한요소법(有限要素法)에 의한 비선형(非線形) 해석시(解析時) 일반적(一般的)인 Lagrangian방법(方法)은 각(各) 구조요소(構造要素)에 대한 유한요소(有限要素) Mapping과 변형(變形) Mapping을 따로 수행(遂行)하여야 하는데 비하여 Pseudo Lagrangian방법(方法)에서는 한번의 직접적(直接的) Mapping으로 구조요소(構造要素)의 변형상태(變形狀態)를 나타낸다. 예제분석(例題分析) 결과(結果)는 이 PLF방법(方法)에 의한 비선형(非線型) 정식화가 종래(從來)의 Lagrangian방법(方法)보다 적용성(適用性)이 양호(良好)하며 유도과정(誘導過程)이 간편(簡便)함을 보여주고 있으며 PLF방법(方法)에 의한 보요소(要素)의 비선형(非線型) 해석(解析)이 매우 정확(正確)한 결과(結果)를 나타내고 있음을 입증(立證)하고

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 추계학술대회
• /
• pp.474-479
• /
• 2004

According to our previous study, it is confirmed that the Petrov-Galerkin natural element method (PGNEM) completely resolves the numerical integration inaccuracy in the conventional Bubnov-Galerkin natural element method (BG-NEM). This paper is an extension of PG-NEM to two-dimensional nonlinear dynamic problem. For the analysis, a constant average acceleration method and a linearized total Lagrangian formulation is introduced with the PG-NEM. At every time step, the grid points are updated and the shape functions are reproduced from the relocated nodal distribution. This process enables the PG-NEM to provide more accurate and robust approximations. The representative numerical experiments performed by the test Fortran program, and the numerical results confirmed that the PG-NEM effectively and accurately approximates the nonlinear dynamic problem.

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

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

• 한국전산구조공학회논문집
• /
• 제18권2호
• /
• pp.123-131
• /
• 2005

기존의 부브노프-갤러킨 자연요소법(BG-NEM)에서 발생하는 수치적분의 부정확성을 페트로프-갤러킨 자연요소법(PG-NEM)에서 완벽히 해결할 수 있음을 저자들의 이전 논문에서 확인하였다. 본 논문에서는 PG-NEM을 확장하여 2차원 기하학적 비선형 문제를 다룬다. 해석을 위해 선형화된 토탈 라그랑지 정식화를 도입하고 PG-NEM을 적용하여 근사화한다. 각 하중 단계마다 절점은 새로운 위치로 갱신되며, 재분포된 절점을 바탕으로 형상함수를 새롭게 구성한다. 이러한 과정은 PG-NEM이 더 정확하고 안정적인 근사함수를 제공하는 것을 가능하게 한다. 개발된 포트란 시험 프로그램을 이용하여 대표적인 수치 예제를 수행하였으며, 수치결과로부터 PG-NEM이 효율적이고 정확하게 대변형 문제를 근사화하는 것을 확인하였다.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2005년도 추계학술대회논문집
• /
• pp.788-791
• /
• 2005

Nodal displacements are referred to the Initial configuration in the total Lagrangian formulation and to the last converged configuration in the updated Lagrangian formulation. This research proposes a relative nodal displacement method to represent the position and orientation for a node in truss structures. Since the proposed method measures the relative nodal displacements relative to its adjacent nodal reference frame, they are still small for a truss structure undergoing large deformations for the small size elements. As a consequence, element formulations developed under the small deformation assumption are still valid fer structures undergoing large deformations, which significantly simplifies the equations of equilibrium. A structural system is represented by a graph to systematically develop the governing equations of equilibrium for general systems. A node and an element are represented by a node and an edge in graph representation, respectively. Closed loops are opened to form a spanning tree by cutting edges. Two computational sequences are defined in the graph representation. One is the forward path sequence that is used to recover the Cartesian nodal displacements from relative nodal displacements and traverses a graph from the base node towards the terminal nodes. The other is the backward path sequence that is used to recover the nodal forces in the relative coordinate system from the known nodal forces in the absolute coordinate system and traverses from the terminal nodes towards the base node. One closed loop structure undergoing large deformations is analyzed to demonstrate the efficiency and validity of the proposed method.

• 한국복합재료학회:학술대회논문집
• /
• 한국복합재료학회 2001년도 춘계학술발표대회 논문집
• /
• pp.163-167
• /
• 2001

A two-dimensional progressive failure analysis method is presented for the strength characterization of the composite joints under pin loading. The eight-nodes laminated she]1 element is utilized based on the updated Lagrangian formulation. The criteria by Yamada-Sun, Tsai-Wu, and the maximum stress are used for the failure estimation. The stiffness of failed layer is degraded by the complete unloading method. No factor depending on test is included in the finite element analysis except for the material strength and stiffness. Total 20 plate specimens with and without hole are tested to validate the finite element prediction. The Tsai-Wu failure criterion most conservatively estimates the strength of laminate, and the maximum stress criterion yields the highest strength because it does not consider the coupling of the failure modes. The strength by Yamada-Sun method neglecting the matrix failure effect are located between other two methods and shows best agreement with test result for laminate with hole.

• Selected Papers of The Society of Naval Architects of Korea
• /
• 제1권1호
• /
• pp.91-100
• /
• 1993

A displacement-based finite element method Is presented for the geometrically nonlinear analysis of eccentrically stiffened plates. A nonlinear degenerated shell element and a nonlinear degenerated eccentric isoparametric beam (isobeam) element are formulated on the basis of Total Agrangian and Updated Lagrangian descriptions. In the formulation of the isobeam element, some additional local decrees of freedom are implementd to describe the stiffener's local plate buckling modes. Therefore this element can be effectively employed to model the eccentric stiffener with fewer D.O.F's than the case of a degenerated shell element. Some detailed buckling and nonlinear analyses of an eccentrically stiffened plate are performed to estimate the critical buckling loads and the post buckling behaviors including the local plate buckling of the stiffeners discretized with the degenerated shell elements and the isobeam elements. The critical buckling loads are found to be higher than the analytical plate buckling load but lower than Euler buckling load of the corresponding column, i.e, buckling strength requirements of the Classification Societies for the stiffened plates.

• 대한기계학회논문집A
• /
• 제29권4호
• /
• pp.534-539
• /
• 2005

Nodal displacements are referred to the initial configuration in the total Lagrangian formulation and to the last converged configuration in the updated Lagrangian furmulation. This research proposes a relative nodal displacement method to represent the position and orientation for a node in truss structures. Since the proposed method measures the relative nodal displacements relative to its adjacent nodal reference frame, they are still small for a truss structure undergoing large deformations for the small size elements. As a consequence, element formulations developed under the small deformation assumption are still valid for structures undergoing large deformations, which significantly simplifies the equations of equilibrium. A structural system is represented by a graph to systematically develop the governing equations of equilibrium for general systems. A node and an element are represented by a node and an edge in graph representation, respectively. Closed loops are opened to form a spanning tree by cutting edges. Two computational sequences are defined in the graph representation. One is the forward path sequence that is used to recover the Cartesian nodal displacements from relative nodal displacement sand traverses a graph from the base node towards the terminal nodes. The other is the backward path sequence that is used to recover the nodal forces in the relative coordinate system from the known nodal forces in the absolute coordinate system and traverses from the terminal nodes towards the base node. One open loop and one closed loop structure undergoing large deformations are analyzed to demonstrate the efficiency and validity of the proposed method.