• 설비공학논문집
• /
• 제3권2호
• /
• pp.97-102
• /
• 1991

Finite element method has been developed recently for the solution of the convection-diffusion problems. Finite element method has several advantages over finite difference method, but its requirement of the larger memory size of the computer has prevented from wide application. In the present study, line-by-line technique has been implemented to finite element method to overcome this disadvantage. Two dimensional laminar natural convection in square cavity was chosen as an example in this study. The numerical result shows good agreement with bench mark solution and the size of the coefficient marix has been reduced drastically.

• Structural Engineering and Mechanics
• /
• 제76권3호
• /
• pp.379-393
• /
• 2020

In this paper, we propose an automatic procedure to improve the accuracy of finite element solutions using enriched 2D solid finite elements (4-node quadrilateral and 3-node triangular elements). The enriched elements can improve solution accuracy without mesh refinement by adding cover functions to the displacement interpolation of the standard elements. The enrichment scheme is more effective when used adaptively for areas with insufficient accuracy rather than the entire model. For given meshes, an error for each node is estimated, and then proper degrees of cover functions are applied to the selected nodes. A new error estimation method and cover function selection scheme are devised for the proposed adaptive enrichment scheme. Herein, we demonstrate the proposed enrichment scheme through several 2D problems.

• Structural Engineering and Mechanics
• /
• 제83권3호
• /
• pp.353-361
• /
• 2022

Linear Elastic Fracture Mechanics (LEFM) has been developed by applying stress analysis to determine the stress intensity factor (SIF, K). The finite element method (FEM) is widely used as a standard tool for evaluating the SIF for various crack configurations. The prediction accuracy can be achieved by applying an adaptive Delaunay triangulation combined with a FEM. The solution can be solved using either direct or iterative solvers. This work adopts the element-by-element preconditioned conjugate gradient (EBE-PCG) iterative solver into an adaptive FEM to solve the solution to heal problem size constraints that exist when direct solution techniques are applied. It can avoid the formation of a global stiffness matrix of a finite element model. Several numerical experiments reveal that the present method is simple, fast, and efficient compared to conventional sparse direct solvers. The optimum convergence criterion for two-dimensional LEFM analysis is studied. In this paper, four sample problems of a two-edge cracked plate, a center cracked plate, a single-edge cracked plate, and a compact tension specimen is used to evaluate the accuracy of the prediction of the SIF values. Finally, the efficiency of the present iterative solver is summarized by comparing the computational time for all cases.

• Structural Engineering and Mechanics
• /
• 제83권2호
• /
• pp.273-282
• /
• 2022

In this paper, we propose a new concept for error estimation in finite element solutions, which we call mode-based error estimation. The proposed error estimation predicts a posteriori error calculated by the difference between the direct finite element (FE) approximation and the recovered FE approximation. The mode-based FE formulation for the recently developed self-updated finite element is employed to calculate the recovered solution. The formulation is constructed by searching for optimal bending directions for each element, and deep learning is adopted to help find the optimal bending directions. Through various numerical examples using four-node quadrilateral finite elements, we demonstrate the improved predictive capability of the proposed error estimator compared with other competitive methods.

• Advances in Computational Design
• /
• 제5권3호
• /
• pp.305-321
• /
• 2020

In this paper is presented the solution method for three-dimensional problem of transversely isotropic body's elastoplastic deformation by the finite element method (FEM). The process of problem solution consists of: determining the effective parameters of a transversely isotropic medium; construction of the finite element mesh of the body configuration, including the determination of the local minimum value of the tape width of non-zero coefficients of equation systems by using of front method; constructing of the stiffness matrix coefficients and load vector node components of the equation for an individual finite element's state according to the theory of small elastoplastic deformations for a transversely isotropic medium; the formation of a resolving symmetric-tape system of equations by summing of all state equations coefficients summing of all finite elements; solution of the system of symmetric-tape equations systems by means of the square root method; calculation of the body's elastoplastic stress-strain state by performing the iterative process of the initial stress method. For each problem solution stage, effective computational algorithms have been developed that reduce computational operations number by modifying existing solution methods and taking into account the matrix coefficients structure. As an example it is given, the problem solution of fibrous composite straining in the form of a rectangle with a system of circular holes.

• East Asian mathematical journal
• /
• 제38권5호
• /
• pp.603-610
• /
• 2022

In [6, 7] they introduced a new finite element method for accurate numerical solutions of Poisson equations with corner singularities. They consider the Poisson equations with homogeneous boundary conditions, compute the finite element solutions using standard FEM and use the extraction formula to compute the stress intensity factor(s), then they posed new PDE with a regular solution by imposing the nonhomogeneous boundary condition using the computed stress intensity factor(s), which converges with optimal speed. From the solution they could get an accurate solution just by adding the singular part. They considered a partial differential equation with the input function f ∈ L²(Ω). In this paper we consider a PDE with the input function f ∈ H¹(Ω) and find the corresponding singular and dual singular functions. We also induce the corresponding extraction formula which are the basic element for the approach.

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

In this study, the low velocity impact behavior of the composite laminates has been described by using 3 dimensional nonlinear finite elements. To describe the geometric nonlinearity due to large deformation, the dynamic contact problem is formulated using the exterior penalty finite element method on the base of Total Lagrangian formulation. The incremental decomposition is introduced, and the converged solution is attained by Newton-Raphson Method. The Newmark's constant-acceleration time integration algorithm is used. To make verification of the finite element program developed in this study, the solution of the nonlinear static problem with occurrence of large deformation is compared with ABAQUS, and the solution of the static contact problem with indentation is compared with the Hertz solution. And, the solution of low velocity impact problem for isotropic material is verificated by comparison with that of LS-DYNA3D. Finally the contact force of impact response from the nonlinear analysis are compared with those from the linear analysis.

• 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.

• 대한기계학회논문집
• /
• 제19권11호
• /
• pp.2971-2980
• /
• 1995

A numerical study of fluid flow in driven cavity was carried out using singular finite element method. The driven cavity problem is known to have infinite velocity gradients as well as dual velocity conditions at the singular points. To overcome such difficulties, a finite element method with singular shape functions was used and a special technique was employed to allow multiple values of velocities at the singular points. Application of singular elements in the driven cavity problem has a significant influence on the stability of solution. It was found the singular elements gave a stable solution, especially, for the pressure distribution of the entire flow field by keeping up a large pressure at the singular points. In the existing solutions of driven cavity problem, most efforts were focused on the study of streamlines and vorticities, and pressure were seldom mentioned. In this study, however, more attention was given to the pressure distribution. Computations showed that pressure decreased very rapidly as the distance from the singular point increased. Also, the pressure distribution along the vertical walls showed a smoother transition with singular elements compared to those of conventional method. At the singular point toward the flow direction showed more pressure increase compared with the other side as Reynolds number increased.

• 소음진동
• /
• 제7권1호
• /
• pp.61-69
• /
• 1997

Despite of the development in the finite element method, it is difficult to get the finite element model describing the dynamic characteristics of the complex structure exactly. Therefore a number of different methods have been developed in order to update the finite element model of a structure using vibration test data. This paper outlines the basic formulation for the frequency response function based updating method. One important advantage of this method is that the intermediate step of performing an eigensolution extraction is unnecessary. Using simulated experimental data, studies are conducted in the case of 10 DOF discrete system. The solution of noisy and incomplete experimental data is discussed. True measured frequency response function data are used for updating the finite element model of a beam and a plate. Its applicability to the joint identification is also considered.