• Structural Engineering and Mechanics
• /
• v.8 no.1
• /
• pp.19-26
• /
• 1999

Superior performance of field consistent eight-node hexahedron element in static bending problems has already been demonstrated in literature. In this paper, its performance in free vibration is investigated. Free vibration frequencies of typical test problems have been computed using this element. The results establish its superior performance in free vibration, particularly in thin plate application and near incompressibility regimes, demonstrating that shear locking, Poisson's stiffening and volumetric locking have been eliminated.

• Economic and Environmental Geology
• /
• v.41 no.6
• /
• pp.731-746
• /
• 2008

Most multiphase flow simulators following fractional flow approach assume incompressibility of fluid and matrix or consider only two phase flow (water and air, water and NAPL). However, in this study, mathematical governing equations were developed for fully compressible three-phase flow using fractional flow based approach. Also, fully compressible multiphase flow simulator (CMPS) considering compressibilities of matrix and fluid was developed using the mathematical governing equations. In order to verify CMPS, the CMPS were compared with analytical solution and the existing multiphase flow simulator, MPS, which had been developed for simulating incompressible multiphase flow (Suk and Yeh 2007; Suk and Yeh 2008). According to the results, solutions of CMPS and MPS and analytical solutions are well matched each other. Thus, it is found that CMPS has the capability of simulating compressible three phase flow phenomena assuming compressibilities of fluids and matrix.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.18 no.2
• /
• pp.83-90
• /
• 2010

The windshield wiper consists of 4 parts: a blade, an arm, a linkage and a motor. The wiper blade makes contact with the windshield and is designed to be operated normally at an angle of 30~50 degrees to the front glass. If the contact pressure between the wiper blade and windshield surface is too high, noise and wear of the rubber will result. On the other hand, if the contact pressure is too low, the performance will do badly, since foreign substances such as dust and stains will not be removed well. The pressure and friction of the wiper blade has a great influence on its effectiveness in cleaning the front window. This is due to the contact of the rubber with the window. This paper presents the dynamic analysis method to estimate the performance of the flat type blade of the wiper system. The blade has a nonlinear characteristic since the rubber is an incompressible hyper-elastic and visco-elastic material. Thus, Structural dynamic analysis using a complex contact model for the blade is performed to find the characteristics of the blade. The flexible multi-body dynamic model is verified by the comparison between test and analysis result. Also, the optimization using the central composite design table is performed.

• Transactions of Materials Processing
• /
• v.8 no.4
• /
• pp.399-406
• /
• 1999

It is more appropriate to treat that the semi-solid mixture as a single phase material that obeys incompressibility in the global sense and to analyze the liquid flow only locally than the approach based on compressible yield criteria. In the present study, a numerical algorithm of updating the solid volume fraction based on mixture theory has been developed. Finite element analysis of simple upsetting was carried out using the proposed algorithm to investigate the degree of macro-segregation according to friction conditions and compressive strain rates under the isothermal condition. The simulation results were compared to experimental results available in reference to test the validity of the currently proposed algorithm. Since the comparison results show a good agreement it is construed that the proposed algorithm can contribute to the development of numerical analysis of determining the solid volume fraction semi-solid processing.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.22 no.1
• /
• pp.121-131
• /
• 1998

This paper applies the FE analysis procedure, developed in the Part I of the companion article, to the three-dimensional rubber pad deformation during rubber-pad forming process. Effects of different algorithms corresponding to incompressibility constraint and time integration methods on numerical solution responses are investigated. Laboratory scale experiments support the validity of the developed FE procedure an demonstrate the accuracy of the numerical models. Full scale model responses are also predicted using the reasonable method and parameters obtained in laboratory modeling.

• Journal of Advanced Marine Engineering and Technology
• /
• v.27 no.4
• /
• pp.494-502
• /
• 2003

The near-tip field of mode-I dynamic cracks steadily propagating in a strain softening material is investigated under plane strain conditions. The material is assumed to be incompressible and its deformation obeys the $J_2$ flow theory of plasticity. A power-law stress-strain relation with strain softening is adopted to account for the damage behavior of materials near the dynamic crack tip. By assuming that the stresses and strain have the same singularity at the crack tip. this paper obtains a fully continuous dynamic crack-tip field in the damage region. Results show that the stress and strain components the same logarithmic singularity of (In(R/r))$\delta$, and the angular variations of filed quantities are identical to those corresponding to the dynamic cracks in the elastic-perfectly plastic material.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.22 no.2
• /
• pp.365-371
• /
• 1998

In the present study, we attempted to add the incompatible functions as additional variable terms to the conventional u-p formulation. It is derived from the four-field generalized variational principle that encompasses velocity, pressure, velocity strains and stress fields as independent interpolated variables. As a severe test of the present formulation, we have investigated the driven cavity with the corner velocity singularity like leaky lid. Through the test, the present element performs very well without unstable oscillation of pressure distribution.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.11
• /
• pp.2312-2321
• /
• 2002

The first-order least-squares meshfree method for linear elasticity is presented. The conventional and the compatibility-imposed least-squares formulations are studied on the convergence behavior of the solution and the robustness to integration error. Since the least-squares formulation is a type of mixed formulation and induces positive-definite system matrix, by using shape functions of same order for both primal and dual variables, higher rate of convergence is obtained for dual variables than Galerkin formulation. Numerical examples also show that the presented formulations do not exhibit any volumetric locking for the incompressible materials.

• Journal of the Korean Mathematical Society
• /
• v.38 no.1
• /
• pp.1-23
• /
• 2001

This paper is concerned with an optimal shape control problem for the stationary Navier-Stokes system. A two-dimensional channel flow of an incompressible, viscous fluid is examined to determine the shape of a bump on a part of the boundary that minimizes the viscous drag. by introducing an artificial compressibility term to relax the incompressibility constraints, we take the penalty method. The existence of optima solutions for the penalized problem will be shown. Next, by employing Lagrange multipliers method and the material derivatives, we derive the shape gradient for the minimization problem of the shape functional which represents the viscous drag.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.112-122
• /
• 1996

A finite element scheme using the concept of PISO method has been developed to solve the Navier-Stokes viscous flows in all speed range. This scheme includes development of new pressure equation that retains both the hyperbolic term related with the density variation and the elliptic term reflecting the incompressibility constraint. The present method is applied to the incompressible two-dimensional driven cavity flow problems(Re=100, 400 and 1,000). For compressible flows, the Carter plate problem(M=3 and Re=1,000) is computed. Finally, we have simulated the shock-boundary layer interaction(M=2 and Re=2.96×10/sup 5/), a more difficult problem, and compared its results with the experiment to demonstrate the shock capturing capability of the present solution algorithm.