• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 춘계학술대회논문집
• /
• pp.597-600
• /
• 2008

Preconditioners for a two-dimensional combined finite element formulation have been devised and tested for fluid-structure interaction (FSI) problems. The FSI code simulating the interaction of a elastic body with an unsteady flow is based on P2P1 finite element combined formulation. It has been shown that two preconditioners among them perform well with respect to computational memory and convergence for a bench-mark problem. Based on the verification of the preconditioners for the two-dimensional combined formulation, four preconditioners are proposed for the problem of an elastic body interacting with a flow.

• Interaction and multiscale mechanics
• /
• 제6권2호
• /
• pp.237-255
• /
• 2013

In this paper, a meshfree-enriched finite element method (ME-FEM) is introduced for the large deformation analysis of nonlinear path-dependent problems involving contact. In linear ME-FEM, the element formulation is established by introducing a meshfree convex approximation into the linear triangular element in 2D and linear tetrahedron element in 3D along with an enriched meshfree node. In nonlinear formulation, the area-weighted smoothing scheme for deformation gradient is then developed in conjunction with the meshfree-enriched element interpolation functions to yield a discrete divergence-free property at the integration points, which is essential to enhance the stress calculation in the stage of plastic deformation. A modified variational formulation using the smoothed deformation gradient is developed for path-dependent material analysis. In the industrial metal forming problems, the mortar contact algorithm is implemented in the explicit formulation. Since the meshfree-enriched element shape functions are constructed using the meshfree convex approximation, they pose the desired Kronecker-delta property at the element edge thus requires no special treatments in the enforcement of essential boundary condition as well as the contact conditions. As a result, this approach can be easily incorporated into a conventional displacement-based finite element code. Two elasto-plastic problems are studied and the numerical results indicated that ME-FEM is capable of delivering a volumetric locking-free and pressure oscillation-free solutions for the large deformation problems in metal forming analysis.

• Advances in aircraft and spacecraft science
• /
• 제1권2호
• /
• pp.145-175
• /
• 2014

The results of a series of numerical experiments are presented to verify some of the important developments made in the first part of this paper. Firstly, the static solution of an algebraic system obtained through Strong Formulation Finite Element Method (SFEM) is presented. Secondly, the stress and strain recovery procedure is descripted for the present technique. It will be clear that the present approach is suitable for any strong formulation finite element methodology, due to the presented general approach based on the unknown displacements and on the elasticity equations. Thirdly, the numerical solutions for some classical and other numerical results found in literature are exposed. Finally, an arbitrarily shaped composite plate is solved and good agreement is observed for all the presented cases.

• 한국자동차공학회논문집
• /
• 제5권6호
• /
• pp.141-147
• /
• 1997

For the effective finite element analysis of the structures including material interfaces or contact surfaces, interface elements are proposed. Most of early works in this problem require not only iterative computation but also complex formulation because of the kinematic nonlinearities caused from the discontinuous behavior and the stress concentration phenomena. The proposed elements, however, are consistently formulated using relative displacements and tractions between top and bottom regular finite elements. The effectiveness of these elements are shown by solving various numerical sample problems including an leaf spring and comparing with results of general finite element analysis. As a result, more stable solutions are conveniently obtaines using interface elements than regular finite elements.

• Interaction and multiscale mechanics
• /
• 제6권2호
• /
• pp.83-105
• /
• 2013

In this paper, a multi-scale meshfree-enriched finite element formulation is presented for the analysis of acoustic wave propagation problem. The scale splitting in this formulation is based on the Variational Multi-scale (VMS) method. While the standard finite element polynomials are used to represent the coarse scales, the approximation of fine-scale solution is defined globally using the meshfree enrichments generated from the Generalized Meshfree (GMF) approximation. The resultant fine-scale approximations satisfy the homogenous Dirichlet boundary conditions and behave as the "global residual-free" bubbles for the enrichments in the oscillatory type of Helmholtz solutions. Numerical examples in one dimension and two dimensional cases are analyzed to demonstrate the accuracy of the present formulation and comparison is made to the analytical and two finite element solutions.

• 한국지반공학회논문집
• /
• 제16권4호
• /
• pp.191-199
• /
• 2000

본 연구에서는 피에조콘 시험의 유한요소해석을 수행하였다. 이를 위하여 점탄소성 bounding surface 모델, 가상일의 방정식(virtual work equation) 및 theory of mixtures를 Updated Lagrangian reference frame에서 formulation하였다. 결과적으로 구성된 유한요소 formulation을 컴퓨터 프로그래밍 하였으며 유한요소해석에서 얻은 콘 저항치, 과잉간극수압 및 간극수압소산 등의 결과를 실험치와 비교 분석하였으며 피에조콘 주변의 응력, 변형율 및 과잉간극수압의 contour를 유한요소해석에서 구하여 이를 고찰하였다. 비등방성 및 점성이 추가된 구성모델을 사용함으로서 응력의 비등방성 및 관입속도를 효과적으로 simulation할 수 있었다. 유한요소 Formulation 과정은 'I' 결과는 'II'에서 설명된다.

• 대한기계학회논문집
• /
• 제13권5호
• /
• pp.820-830
• /
• 1989

본 논문에서는 타이어의 각 부분의 물성치 계산을 위한 식을 유한요소법에 적용할 수 있도록 제안하였다. 이 식은 강철 코드의 굽힙효과를 고려 하였으며, 특히 각 요소에서 전단변형이 일어나는 동안의 굽힘효과를 고려하였다. 유한요소 공식화는 가상일의 원리에 의하여 평형 방정식으로부터 유도하였고, Updated refer- ence coordinate에 대해 증분해석을 적용하여 Updated Lagrangian공식화를 하였다. 그리고 차량하중에 의하여 타이어가 노면에 접지될때의 응력상태를 게산할 수 있도록 접촉문제 공식화를 유한요소 공식화에 첨가 하였다.

• 한국CDE학회논문집
• /
• 제15권6호
• /
• pp.411-417
• /
• 2010

In this paper, the boom of a floating crane is modeled as a 3-dimensional elastic beam in order to analyze the dynamic response of the crane and its cargo. The boom is divided into more than two elements based on finite element formulation, and deformation of each element is expressed in terms of shape matrix and nodal coordinates. The equations of motion for the elastic boom consist of a mass matrix, a stiffness matrix, and a quadratic velocity vector that contains the gyroscopic and Coriolis forces. The size and complicity of the matrices increase in proportion with the number of elements. Therefore, it is not possible to derive the equations of motion explicitly for different number of elements. To overcome this difficulty, matrices for one 3-dimensional element are expressed with elementary sub-matrices. In particular, the quadratic velocity vector is derived as a product of a shape matrix and a 3-dimensional rotation matrix. By using the derived matrices, the equations of motion for the multi-element boom are automatically constructed. To verify the implementation of the elastic boom based on finite element formulation, we simulated a simple vibration of the elastic boom and compared the average deformation with the analytic solution. Finally, heave motion of the floating crane and surge motion of the cargo are presented as application examples of the elastic boom.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 1992년도 추계학술대회논문집; 반도아카데미, 20 Nov. 1992
• /
• pp.117-123
• /
• 1992

In recent years, the finite element method has become one of the most popular numerical technique for obtaining solutions of engineering science problems. However, there exist various uncertainties in modeling the problems, such as the dimensions(geometry shape), the material properties, boundary conditions, etc. The consideration for the uncertainties inherent in the problems can be made by understanding the influences of uncertain parameters[1]. Determining the influences of uncertainties as statistical quantities using the standard finite element method requires enormous computing time, while the probabilistic finite element method is realized as an efficient scheme[2,3] yielding statistical solution with just a few direct computations. In this paper, a formulation of the probabilistic fluid-structure interaction problem accounting for the first order perturbation of geometric shape is derived, and especially probabilistical acoustic pressure scattering from the structure with surrounding fluid is focused on. In Section 2, governing equations for the fluid-structure problems are given. In Section 3, a finite element formulation, based on the functional, is presented. First order perturbation of geometric shape with randomness is incorporated into the finite element formulation in conjunction with discretization of the random fields in Section 4 and 5. Finally, the proposed formulation is applied to a acoustic pressure scattering problem from an infinitely long cylindrical shell structure with randomness of radial perturbation.

• Structural Engineering and Mechanics
• /
• 제69권6호
• /
• pp.615-626
• /
• 2019

The 3-noded metric Timoshenko beam element with an offset of the internal node from the element centre is used here to demonstrate the best-fit paradigm using function space formulation under locking and mesh distortion. The best-fit paradigm follows from the projection theorem describing finite element analysis which shows that the stresses computed by the displacement finite element procedure are the best approximation of the true stresses at an element level as well as global level. In this paper, closed form best-fit solutions are arrived for the 3-noded Timoshenko beam element through function space formulation by combining field consistency requirements and distortion effects for the element modelled in metric Cartesian coordinates. It is demonstrated through projection theorems how lock-free best-fit solutions are arrived even under mesh distortion by using a consistent definition for the shear strain field. It is shown how the field consistency enforced finite element solution differ from the best-fit solution by an extraneous response resulting from an additional spurious force vector. However, it can be observed that when the extraneous forces vanish fortuitously, the field consistent solution coincides with the best-fit strain solution.