• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.6
• /
• pp.916-922
• /
• 2001

A mechanical face seal is a tribo-element intended to control leakage of working fluid at the interface between a rotating shaft and its housing. Leakage of working fluid decreases drastically as the clearance of the mating seal faces gets smaller. But the very small seal clearance results in an increased reduction of seal life because of high wear and heat generation. Therefore, in the design of mechanical face seals a compromise between low leakage and acceptable seal life is important, and it presents a difficult and practical design problem. A fluid film or sealing dam geometry of the seal clearance affects seal lubrication performance very much, and thereby it is one of the main design considerations. In this study the Reynolds equation for the sealing dam of mechanical face seals is numerically analyzed, using the Galerkin finite element method, which is readily applied to various seal geometries. Film pressures of the sealing dam are analyzed, including the effects of the seal face coning and tilt. Then, lubrication performances of the seals, such as opening forces, restoring moments, leakage, and dynamic coefficients, are calculated, and they are compared to the results obtained by the narrow seal approximation.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.7
• /
• pp.2067-2077
• /
• 1996

In adaptive finite element analysis, r- and h-methods are generally used on the basis of a discretization error estimator. In this paper, an rh-method is proposed as a new adaptive method which can improve the adaptivity performance by using both of them. This suggested rh-method moves nodal coordinates of initially given model to adjust element discretization errors and thereafter performes the h-method tdo obtain the specified accuracy of finite element solutions. Numerical experiments for various plane problems were performed using 4-noded isoparametric quadrilateral elements. As a result, the rh-method has been shown to be an accurate and efficient adaptive analysis method to obtain as improved solution.

• Journal of Mechanical Science and Technology
• /
• v.19 no.spc1
• /
• pp.283-291
• /
• 2005

Accurate seismic analyses of large deformable moving structures are still unsolved problems in the field of earthquake engineering. In order to analyze these problems, the nonlinear finite element method formulated by the absolute nodal coordinate approach is noticed. Because, this formulation has several advantages over the standard procedures on mass matrix, elastic forces and damping forces in the case of large displacement problems. But, it has not been fully studied to build frame structure models by using beam elements in the absolute nodal coordinate formulation. In this paper, we propose the connecting method of the beam elements formulated by the absolute nodal coordinate. The coordinate transformation matrix of this element is introduced into the frame structure. This beam element has the characteristic that the mass matrix and bending stiffiness matrix are constant even if in the case of large displacement problems, and this characteristic is being kept after the transformation. In order to verify the proposed method, we show the numerical simulation results of frame structures for a vibration problem and a large displacement problem.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.9
• /
• pp.2887-2893
• /
• 1996

Governing equations infrictional contact problems are introduced using variational inequality formulations which are regularized to overcome the diffculties of non-differentiability of the friction functional. Also finite element approximations and a priori error estimations are derived based on those formulations. Numerical simulations are performed illustrating the theoretical results.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2005.11a
• /
• pp.66-69
• /
• 2005

Mechanical strength of functionally graded composite plates that composed of ceramic, functionally graded material and metal layers is investigated using 3-D finite element method. In FGM layer, material properties are assumed to be varied continuously in the thickness direction according to a simple power law distribution in terms of the volume fraction of a ceramic and metal. The 3-D finite element model is adopted by using an IS-node solid element to analyze more accurately the variation of material properties in the thickness direction. Numerical results are compared with those of the previous works. In addition, the displacements, the tensile stresses and the compressive stresses are analyzed for the variation of FGM thickness ratio and volume fraction distribution.Mechanical strength of functionally graded composite plates that composed of ceramic, functionally graded material and metal layers is investigated using 3-D finite element method. In FGM layer, material properties are assumed to be varied continuously in the thickness direction according to a simple power law distribution in terms of the volume fraction of a ceramic and metal. The 3-D finite element model is adopted by using an IS-node solid element to analyze more accurately the variation of material properties in the thickness direction. Numerical results are compared with those of the previous works. In addition, the displacements, the tensile stresses and the compressive stresses are analyzed for the variation of FGM thickness ratio and volume fraction distribution.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.1
• /
• pp.163-171
• /
• 1996

For the effective analysis of two dimensional plane problems, Treffiz finite elements and cavity elements have been proposed. These element matrix equaitons were formulated on the basis of hybrid variational principle and Treffiz function sets derived consitstently from the complex theoy of plane elasticity. In order to suggest the accuracy chatacteristics of the proposed Treffiz elements typical plane problems were analyzed and these results were compared with ones obtained by using the conveintional displacement type elements. The accuracy of the proposed elements is less sensitive to the element size and shape than the conventional displacement type elements. These elements, being able to be formed with multi-nodes, give the convenient modeling of an analytic domain. The cavity elements give the comparatively exact values of stress concentration factors of stress intensity factors and can be effectively used for the analysis of mechanical stuctures containing various cavities.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.4
• /
• pp.1133-1145
• /
• 1996

A method is proposed to predict the deformed shape of the structure subjected to the unknown external loads using the signal from the piezoceramic sensors. Such a shape estimation is based on the linear relationship between the deformation of structure and the signal from sensor, which is calculated using finite element method. The deformed shape is, then calculated using the linear matrix and the signals from the piezoceramic sensors attached to the structures. For the purpose, a structural analysis program is developed using a multi-layerd finite element of 8 nodes with 3 displacement and one voltage degrees of freedom at each node. The multiple layers with the different material properties can be layered within the element. The incompatible mode with the element is found to be crucial to catch the bending behavior accurately. The accuracy of the program is, then, verified by being compared with the experimental results performed by Crawley. The proposed shape estimation method is also verified for the different loads and sensor size. It is shown that the results of shape estimation method using the linear matrix well predicts the deflections compared with those of finite element method.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.11
• /
• pp.2922-2931
• /
• 1994

We propose a modified 6-node element, where the sampling point of Gauss quadrature moved in the thickness direction. The modified 6-node element has been applied to static problems and forced motion analyses. In this study, this method is extended to the finite element analysis of the natural frequencies of two dimensional problems. We also propose a modified 16-node element for three dimensional problems, which behaves much like a 20-node element with smaller degree of freedom. The modified 6-node and 16-node elements have been applied to the modal analyses of beams and plates, respectively. The results agree well with the results of the 8-node or 20-node element models.

• Journal of Mechanical Science and Technology
• /
• v.20 no.1
• /
• pp.94-109
• /
• 2006

An improved meshfree method called the Petrov-Galerkin natural element (PG-NE) method is introduced in order to secure the numerical integration accuracy. As in the Bubnov-Galerkin natural element (BG-NE) method, we use Laplace interpolation function for the trial basis function and Delaunay triangles to define a regular integration background mesh. But, unlike the BG-NE method, the test basis function is differently chosen, based on the Petrov-Galerkin concept, such that its support coincides exactly with a regular integration region in background mesh. Illustrative numerical experiments verify that the present method successfully prevents the numerical accuracy deterioration stemming from the numerical integration error.

• Journal of Mechanical Science and Technology
• /
• v.17 no.6
• /
• pp.805-815
• /
• 2003

In order to execute efficiently the free vibration analysis of 2-dimensional structures like plate structures, the author developed the finite element-transfer stiffness coefficient method. This method is based on the combination of the modeling techniques in the FEM and the transfer technique of the stiffness coefficient in the transfer stiffness coefficient method. Numerical results of the simply supported and the elastic supported rectangular plates showed that the present method can be successfully applied to the free vibration analysis of plate structures on a personal computer. We confirmed that, in the case of analyzing the free vibration of rectangular plate structures, the present method is superior to the FEM from the viewpoint of computation time and storage.