• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2004.10a
• /
• pp.855-859
• /
• 2004

We develop a new four-node flat shell element which is accurate, efficient, and suitable to be used on general purpose. The new element has a hybrid Trefftz element with drilling degrees of freedom as a membrane part. We define the two independent displacement field: the internal displacement field that satisfies governing equations in the domain a priori and the boundary displacement field that is usually used as a conventional finite element method. The hybrid Trefftz variational formulation connects these two displacement fields on the boundary of the domain. To add drilling degrees of freedom, we introduce the Allman's quadratic displacement field to the boundary displacement field. As a result, our flat shell element has 6 degrees of freedom per a node. We also use the well-known DKMQ plate bending element for the plate part of the proposed element. The DKMQ element satisfies Mindlin-Reissner‘s plate theory along the edge of the element and gives proper behavior regardless of the thickness. A series of numerical experiments shows that the performance of the new element such as accuracy, rate of convergence, robustness to mesh quality, and so on.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.6
• /
• pp.829-836
• /
• 2001

A dependable boundary reconstruction technique is proposed. The finite element method is used for the analysis of the direct heat conduction problem to realize the deformable grid system. An appropriate strategy for grid update is suggested. A complete sensitivity analysis is performed to obtain the derivatives required for restoration of the optimal boundary. With the result of the sensitivity analysis, the unknown boundary is sought using the sequential quadratic programming. The method is applied to reconstruction of boundaries with sinusoidal, step, and cavity form. The overall performance of the proposed method is examined by comparison between the estimated the exact boundaries.

• Journal of the Korean Society of Safety
• /
• v.21 no.1 s.73
• /
• pp.15-20
• /
• 2006

The finite element alternating method is a convenient and efficient method to analyze three-dimensional cracks embedded in an infinite or a finite body because the method has the property that the uncracked body and cracks can be modeled independently. In this paper the method was applied for fatigue crack growth simulation. A surface crack in a cylinder was considered as an initial crack and the crack configurations and stress intensity factors during the crack growth were obtained. In this paper the finite element alternating method proposed by Nikishkov, Park and Atluri was used after modification. In the method, as the required solution for a crack in an infinite body, the symmetric Galerkin boundary element method formulated by Li and Mear was used. And a crack was modeled as distribution of displacement discontinuities, and the governing equation was formulated as singularity-reduced integral equations.

• Honam Mathematical Journal
• /
• v.40 no.4
• /
• pp.785-794
• /
• 2018

In [8] they introduced a new finite element method for accurate numerical solutions of Poisson equations with corner singularities. They consider the Poisson equations with homogeneous Dirichlet boundary condition with one corner singularity at the origin, and compute the finite element solution using standard FEM and use the extraction formula to compute the stress intensity factor, then pose a PDE with a regular solution by imposing the nonhomogeneous boundary condition using the computed stress intensity factor, which converges with optimal speed. From the solution they could get an accurate solution just by adding the singular part. This approach uses the polar coordinate and the cut-off function to control the singularity and the boundary condition. In this paper we consider Poisson equations with multiple singular points, which involves different cut-off functions which might overlaps together and shows the way of cording in FreeFEM++ to control the singular functions and cut-off functions with numerical experiments.

• Coupled systems mechanics
• /
• v.1 no.1
• /
• pp.19-37
• /
• 2012

In this work, the iterative coupling of finite element and boundary element methods for the investigation of coupled fluid-fluid, solid-solid and fluid-solid wave propagation models is reviewed. In order to perform the coupling of the two numerical methods, a successive renewal of the variables on the common interface between the two sub-domains is performed through an iterative procedure until convergence is achieved. In the case of local nonlinearities within the finite element sub-domain, it is straightforward to perform the iterative coupling together with the iterations needed to solve the nonlinear system. In particular, a more efficient and stable performance of the coupling procedure is achieved by a special formulation that allows to use different time steps in each sub-domain. Optimized relaxation parameters are also considered in the analyses, in order to speed up and/or to ensure the convergence of the iterative process.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.10
• /
• pp.1742-1756
• /
• 1997

The acoustic characteristics of a direct radiator type loudspeaker has been studied in this paper. The natural modes of the speaker cone vibration analyzed numerically by the finite element method have been verified by comparing them with experimental results. The so-ap-proved finite-element model has been used to calculate the vibration response of the cone excited by the voice coil. The vibration displacement of the speaker cone paper has been converted into the vibration velocity and used as a boundary condition for the acoustic analysis. The frequency characteristics, directivity, and sound pressure distribution of the loudspeaker have been calculated by the boundary element method. The numerical results have been verified by the experiments carried out in an anechoic chamber. The variations of the acoustic characteristics due to the changes of some design parameter values can be examined using the numerical model.

• The Journal of the Acoustical Society of Korea
• /
• v.29 no.7
• /
• pp.448-457
• /
• 2010

Vibrations and wave propagations in waveguide structures can be analysed efficiently by using waveguide finite element (WFE) method. The WFE method only models the 2-dimensional cross-section of the waveguide with finite elements so that the size of the model and computing time are much less than those of the 3-dimensional FE models. For cylindrical shells or pipes which have simple cross-sections, the external coupling with fluids can be treated theoretically. For waveguides of complex cross-sectional geometries, however, numerical methods are required to deal with external fluids. In this numerical approach, the external fluid is modelled by the boundary elements (BEs) and connected to WFEs. In order to validate this WFE/BE method, a pipe submerged in water is considered in this study. The dispersion diagrams and point mobilities of the pipe simulated are compared to those that theoretically obtained. Also the acoustic powers radiated from the pipe are predicted and compared in both cases of air and water as an external medium.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.19 no.9
• /
• pp.922-927
• /
• 2009

The transducers used in active sonar on surface ships are packed in a specific geometry in the array drum in order to meet the requirements such as the source level, directional beam pattern, etc. This paper describes the acoustic characteristics of the cylindrical array which is based on a 64 vertical staves arrangement, each stave composed 5 independent transducers. Firstly, the single transducer on the rigid baffle in the water is analyzed with the Finite Element Method. From the result of the FE analysis nodal velocities on the radiation surface is calculated and used with the boundary conditions of the transducers mounted on the array drum. Then the acoustic pressure is calculated in the field points using the Boundary Element Method and the other acoustic informations, the source level, beam pattern, near field and far-field distance, were acquired.

• Advances in Computational Design
• /
• v.5 no.3
• /
• pp.277-289
• /
• 2020

This article describes the technology for constructing of a multiply-connected three-dimensional area's finite element representation. Representation of finite-element configuration of an area is described by a discrete set that consist of the number of nodes and elements of the finite-element grid, that are orderly set of nodes' coordinates and numbers of finite elements. Corresponding theorems are given, to prove the correctness of the solution method. The adequacy of multiply-connected area topology's finite element model is shown. The merging of subareas is based on the criterion of boundary nodes' coincidence by establishing a simple hierarchy of volumes, surfaces, lines and points. Renumbering nodes is carried out by the frontal method, where nodes located on the outer edges of the structure are used as the initial front.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.32 no.10
• /
• pp.836-843
• /
• 2008

An efficient domain/boundary decomposition method is applied for heat transfer problems with non-linear thermal radiation boundaries. The whole domain of solids or structures is considered as set of subdomains, an interface, and radiation interfaces. In a variational formulation, simple penalty functions are introduced to connect an interface or radiation interfaces with neighboring subdomains that satisfy continuity conditions. As a result, non-linear finite element computations due to the thermal radiation boundaries can be localized within a few subdomains or radiation interfaces. Therefore, by setting up suitable solution algorithms for the governing finite element equations, the computational efficiency can be improved considerably. Through a set of numerical examples, these distinguishing characteristics of the present method are investigated in detail.