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

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

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.209-212
• /
• 2002

Since Berkhoff proposed the mild-slope equation in 1972, it has widely been used for calculation of shallow water wave transformation. Recently, it was extended to give an extended mild-slope equation, which includes the bottom slope squared term and bottom curvature term so as to be capable of modeling wave transformation on rapidly varying topography. These equations were derived by integrating the Laplace equation vertically. In the present study, we develop a finite element model to solve the Laplace equation directly while keeping the same computational efficiency as the mild-slope equation. This model assumes the vertical variation of wave potential as a cosine hyperbolic function as done in the derivation of the mild-slope equation, and the Galerkin method is used to discretize . The computational domain was discretized with proper finite elements, while the radiation condition at infinity was treated by introducing the concept of an infinite element. The upper boundary condition can be either free surface or a solid structure. The applicability of the developed model was verified through example analyses of two-dimensional wave reflection and transmission. .

• Journal of applied mathematics & informatics
• /
• 제26권5_6호
• /
• pp.1057-1069
• /
• 2008

We consider singularly perturbed Boundary Value Problems (BVPs) for third and fourth order Ordinary Differential Equations(ODEs) of convection-diffusion type with discontinuous source term and a small positive parameter multiplying the highest derivative. Because of the type of Boundary Conditions(BCs) imposed on these equations these problems can be transformed into weakly coupled systems. In this system, the first equation does not have the small parameter but the second contains it. In this paper a computational method named as 'An asymptotic finite element method' for solving these systems is presented. In this method we first find an zero order asymptotic approximation to the solution and then the system is decoupled by replacing the first component of the solution by this approximation in the second equation. Then the second equation is independently solved by a fitted mesh Finite Element Method (FEM). Numerical experiments support our theoritical results.

• 대한조선학회지
• /
• 제19권4호
• /
• pp.19-29
• /
• 1982

Galerkin's finite element method is applied to a two-dimensional heat convection-diffusion problem arising in the hydrodynamic lubrication of thrust bearings used in naval vessels. A parabolized thermal energy equation for the lubricant, and thermal diffusion equations for both bearing pad and the collar are treated together, with proper juncture conditions on the interface boundaries. it has been known that a numerical instability arises when the classical Galerkin's method, which is equivalent to a centered difference approximation, is applied to a parabolic-type partial differential equation. Probably the simplest remedy for this instability is to use a one-sided finite difference formula for the first derivative term in the finite difference method. However, in the present coupled heat convection-diffusion problem in which the governing equation is parabolized in a subdomain(Lubricant), uniformly stable numerical solutions for a wide range of the Peclet number are obtained in the numerical test based on Galerkin's classical finite element method. In the present numerical convergence errors in several error norms are presented in the first model problem. Additional numerical results for a more realistic bearing lubrication problem are presented for a second numerical model.

• Structural Engineering and Mechanics
• /
• 제39권1호
• /
• pp.21-46
• /
• 2011

The natural frequencies of continuous systems depend on the governing partial differential equation and can be numerically estimated using the finite element method. The accuracy and convergence of the finite element method depends on the choice of basis functions. A basis function will generally perform better if it is closely linked to the problem physics. The stiffness matrix is the same for either static or dynamic loading, hence the basis function can be chosen such that it satisfies the static part of the governing differential equation. However, in the case of a rotating beam, an exact closed form solution for the static part of the governing differential equation is not known. In this paper, we try to find an approximate solution for the static part of the governing differential equation for an uniform rotating beam. The error resulting from the approximation is minimized to generate relations between the constants assumed in the solution. This new function is used as a basis function which gives rise to shape functions which depend on position of the element in the beam, material, geometric properties and rotational speed of the beam. The results of finite element analysis with the new basis functions are verified with published literature for uniform and tapered rotating beams under different boundary conditions. Numerical results clearly show the advantage of the current approach at high rotation speeds with a reduction of 10 to 33% in the degrees of freedom required for convergence of the first five modes to four decimal places for an uniform rotating cantilever beam.

• East Asian mathematical journal
• /
• 제33권5호
• /
• pp.511-525
• /
• 2017

We introduce an extrapolated higher order characteristic finite element method to construct approximate solutions of a Sobolev equation with a convection term. The higher order of convergence in both the temporal direction and the spatial direction in $L^2$ normed space is established and some computational results to support our theoretical results are presented.

• Journal of applied mathematics & informatics
• /
• 제22권1_2호
• /
• pp.223-235
• /
• 2006

By applying the Landau-type transformation, we transform a Stefan problem with nonlinear free boundary condition into a system consisting of a parabolic equation and the ordinary differential equations. Fully discrete finite element method is developed to approximate the solution of a system of a parabolic equation and the ordinary differential equations. We derive optimal orders of convergence of fully discrete approximations in $L_2,\;H^1$ and $H^2$ normed spaces.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1992년도 가을 학술발표회 논문집
• /
• pp.23-28
• /
• 1992

The ultimate capacity of end-plate connection is investigated through nonlinear finite element analysis. The example models are divided into stiffened case and unstiffened one. The refined finite element models are analyzed by utilizing a general purpose structural analysis computer program ADINA and the moment-rotation relationships of the connection are determined. The results are compared with the regression equation deduced by Krishnamurthy. It is planned to deduce a bilinear regression equation through a parametric study on various dimensions of the connection.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 춘계학술대회
• /
• pp.507-512
• /
• 2003

A finite element analysis for near-net-shape forming of Al6061 powder was performed under warm pressing. The advantages of warm compaction by rubber isostatic pressing were discussed to obtain parts with better density distributions. To simulate densification and deformed shape of a powder compact during warm pressing, the elastoplastic constitutive equation based on yield function of Shima-Oyane was implemented into a finite element program(ABAQUS). The hyperelastic constitutive equation based on the Ogden strain energy potential was employed to analyze nonlinear elastic response of rubber. Finite element results were compared with experimental data for Al6061 powder compacts under warm pressing.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제20권3호
• /
• pp.243-259
• /
• 2016

The Richards equation for water movement in unsaturated soil is highly nonlinear partial differential equations which are not solvable analytically unless unrealistic and oversimplifying assumptions are made regarding the attributes, dynamics, and properties of the physical systems. Therefore, conventionally, numerical solutions are the only feasible procedures to model flow in partially saturated porous media. The standard Finite element numerical technique is usually coupled with an Euler time discretizations scheme. Except for the fully explicit forward method, any other Euler time-marching algorithm generates nonlinear algebraic equations which should be solved using iterative procedures such as Newton and Picard iterations. In this study, lumped mass and distributed mass in the frame of Picard and Newton iterative techniques were evaluated to determine the most efficient method to solve the Richards equation with finite element model. The accuracy and computational efficiency of the scheme and of the Picard and Newton models are assessed for three test problems simulating one-dimensional flow processes in unsaturated porous media. Results demonstrated that, the conventional mass distributed finite element method suffers from numerical oscillations at the wetting front, especially for very dry initial conditions. Even though small mesh sizes are applied for all the test problems, it is shown that the traditional mass-distributed scheme can still generate an incorrect response due to the highly nonlinear properties of water flow in unsaturated soil and cause numerical oscillation. On the other hand, non oscillatory solutions are obtained and non-physics solutions for these problems are evaded by using the mass-lumped finite element method.