• Tribology and Lubricants
• /
• v.10 no.3
• /
• pp.54-61
• /
• 1994

The regimes of elastohydrodynamic lubrication at the points where line contacts occur between the vane tip and camring in an oil hydraulic vane pump is studied. A study of the contact conditions in vane pump provided most of the early interest in the possibility of fluid film lubrication in highly loaded contacts. The variation of viscosity with pressure and the elastic deformation associated with the high pressures generated in the contact region are the major causes of the complexity attributed to lubrication behavior. Therefore a numerical solutions to the problem of elastohydrodynamic lubrication of line contact are obtained by using a finite-difference formulation.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.3
• /
• pp.389-396
• /
• 2001

Direct numerical simulation of bubble growth and merger process on a single nucleation site during partial nucleate boiling is performed. The equations governing conservation of mass, momentum and energy are solved using a finite difference method combined with a level set method for capturing the vapor-liquid interface. The level set method is modified to include the effects of phase change at the interface and contact angle at the wall. Also, a simplified formulation for predicting the evaporative heat flux in a thin liquid micro-layer is developed and incorporated into the level set formulation. Based on the numerical results, the bubble growth and merger pattern and its effect on the heat transfer are discussed.

• Journal of the Korea Institute of Information and Communication Engineering
• /
• v.12 no.12
• /
• pp.2253-2258
• /
• 2008

To design micro flow sensor varying depending on temperature of driving heater in the detector of Oxide semiconductor, Markov chain MCM(MCMCM), which is a kind of stochastic and microscopic method, was introduced. The formulation for the thermal transfer equation based on the FDM to obtain the MCMCM solution was performed and investigated, in steady state case. MCMCM simulation was successfully applied, so that its application can be expanded to a three-dimensional model with inhomogeneous material and complicated boundary.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.6
• /
• pp.1465-1477
• /
• 1993

Powder Injection Molding(PM) is an advanced and complicated technology for manufacturing ceramic or metal products making use of a conventional injection molding process, which is generally used for plastic products. Among many technologies involved in the successful PIM, injection molding process is one of the key steps to form a desired shape out of powder/binder mixtures. Thus, it is of great importance to have a numerical tool to predict the powder injection molding filling process. In this regard, a finite element analysis system has been developed for numerical simulations of filling process of powder injection molding. Powder/polymer mixtures during the filling pro cess of injection molding can be rheologically characterized as Non-Newtonian fluids with a so called yield phenomena and have a peculiar feature of apparent slip phenomena on the wall boundaries surrounding mold cavity. Therefore, in the present study, a physical modeling of the filling process of powder/polymer mixtures was developed to take into account both the yield stress and slip phenomena and a finite element formulation was developed accordingly. The numerical analysis scheme for filling simulation is accomplished by combining a finite element method with control volume technique to simulate the movement of flow front and a finite difference method to calculate the temperature distribution. The present study presents the modeling, numerical scheme and some numerical analysis results showing the effect of the yield stress and slip phenomena.

• Steel and Composite Structures
• /
• v.30 no.5
• /
• pp.483-492
• /
• 2019

In the current study, the influence of the initial lateral (sweep) shape and the cross-sectional twist imperfection on the lateral torsional buckling (LTB) response of doubly-symmetric steel I-beams was investigated. The material imperfection (residual stress) was not considered. For this objective, standard European IPN 300 beam with different unbraced span was numerically analyzed for three imperfection cases: (i) no sweep and no twist (perfect); (ii) three different shapes of global sweep (half-sine, full-sine and full-parabola between the end supports); and (iii) the combination of three different sweeps with initial sinusoidal twist along the beam. The first comparison was done between the results of numerical analyses (FEM) and both a theoretical solution and the code lateral torsional buckling formulations (EC3 and AISC-LRFD). These results with no imperfection effects were then separately compared with three different shapes of global sweep and the presence of initial twist in these sweep shapes. Besides, the effects of the shapes of initial global sweep and the inclusion of sinusoidal twist on the critical buckling load of the beams were investigated to unveil which parameter was considerably effective on LTB response. The most compatible outcomes for the perfect beams was obtained from the AISC-LRFD formulation; however, the EC-3 formulation estimated the $P_{cr}$ load conservatively. The high difference from the EC-3 formulation was predicted to directly originate from the initial imperfection reduction factor and high safety factor in its formulation. Due to no consideration of geometric imperfection in the AISC-LFRD code solution and the theoretical formulation, the need to develop a practical imperfection reduction factor for AISC-LRFD and theoretical formulation was underlined. Initial imperfections were obtained to be more influential on the buckling load, as the unbraced length of a beam approached to the elastic limit unbraced length ($L_r$). Mode-compatible initial imperfection shapes should be taken into account in the design and analysis stages of the I-beam to properly estimate the geometric imperfection influence on the $P_{cr}$ load. Sweep and sweep-twist imperfections led to 10% and 15% decrease in the $P_{cr}$ load, respectively, thus; well-estimated sweep and twist imperfections should considered in the LTB of doubly-symmetric steel I-beams.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.7
• /
• pp.1738-1750
• /
• 1994

Both cylindrical cup drawing and square cup drawing are analyzed using membrane analysis as well as shell analysis by the elastic-plastic finite element method. An incremental formulation incorporating the effect of large deformation and normal anisotropy is used for the analysis of elastic-plastic non-steady deformation. The computed results are compared with the existing experimental results to show the validity of the analysis. Comparisons are made in the punch load and distribution of thickness strain between the membrane analysis and the shell analysis for both cylindrical and square cup drawing processes. In punch load, both analyses show very little difference and also show generally good agreement with the experiment. For the cylindrical cup deep drawing, the computed thickness strain of a membrane analysis, however, shows a wide difference with the experiment. In the shell analysis, the thickness strain shows good agrement with the experiment. For the square cup deep drawing, both membrane and shell analyses show a wide difference with experiment, this may be attributable to the ignorance of the shear deformation. Concludingly, it has been shown that the membrane approach shows a limitation for the deep drawing process in which the effect of bending is not negligible and more exact information on the thickness strain distribution is required.

• Journal of Korea Foundry Society
• /
• v.13 no.4
• /
• pp.323-332
• /
• 1993

An implicit finite difference formulation with three methods of latent heat treatment, such as equivalent specific heat method, temperature recovery method and enthalpy method, was applied to solidification analysis. The Neumann problem was solved to compare the numerical results with the exact solution. The implicit solutions with the equivalent specific heat method and the temperature recovery method were comparatively consistent with the Neumann exact solution for smaller time steps, but its error increased with increasing time step, especially in predicting the solidification beginning time. Although the computing time to solve energy equation using temperature recovery method was shorter than using enthalpy method, the method of releasing latent heat is not realistic and causes error. The implicit formulation of phase change problem requires enthalpy method to treat the release of latent heat reasonably. We have modified the enthalpy formulation in such a way that the enthalpy gradient term is not needed, and as a result of this modification, the computation stability and the computing time were improved.

• Ocean Systems Engineering
• /
• v.11 no.1
• /
• pp.59-81
• /
• 2021

The nonlinear formulation using the principle of virtual work-energy for free vibration of a large-sag extensible catenary riser in two dimensions is presented in this paper. A support at one end is hinged and the other is a free-sliding roller in the horizontal direction. The catenary riser has a large-sag configuration in the static equilibrium state and is assumed to displace with large amplitude to the motion state. The total virtual work of the catenary riser system involves the virtual strain energy due to bending, the virtual strain energy due to axial deformation, the virtual work done by the effective weight, and the inertia forces. The nonlinear equations of motion for two-dimensional free vibration in the Cartesian coordinate system is developed based on the difference between the Euler's equations in the static state and the displaced state. The linear and nonlinear stiffness matrices of the catenary riser are obtained and the eigenvalue problem is solved using the Galerkin finite element procedure. The natural frequencies and mode shapes are obtained. The results are validated with regard to the reference research addressing the accuracy and efficiency of the proposed nonlinear formulation. The numerical results for free vibration and the effect of the nonlinear behavior for catenary riser are presented.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.15 no.10 s.89
• /
• pp.977-982
• /
• 2004

In this paper, a dispersive anisotropic perfectly matched layer(APML) is proposed using piecewise linear recursive convolution(PLRC) for finite difference time domain(FDTD) methods. This proposed APML can be utilized for the analysis of a nonlinear dispersive medium as absorbing boundary condition(ABC). The formulation is simple modification to the original AMPL and can be easily implemented. Also it has some advantages of the PLRC approach-fast speed, low memory cost, and easy formulation of multiple pole susceptibility. We applied this APML to 2-D propagation problems in dispersive media such as Debye and Lorentz media The results showed good absorption at boundaries.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.123-130
• /
• 2002

A general formulation for shape design sensitivity analysis over three dimensional beam structure is developed based on a variational formulation of the beam in linear elasticity. Sensitivity formula is derived based on variational equations in cartesian coordinates using the material derivative concept and adjoint variable method for the displacement and Von-Mises stress functionals. Shape variation is considered for the beam shape in general 3-dimensional direction as well as for the orientation angle of the beam cross section. In the sensitivity expression, the end points evaluation at each beam segment is added to the integral formula, which are summed over the entire structure. The sensitivity formula can be evaluated with generality and ease even by employing piecewise linear design velocity field despite the bending model is fourth order differential equation. For the numerical implementation, commercial software ANSYS is used as analysis tool for the primal and adjoint analysis. Once the design variable set is defined using ANSYS language, shape and orientation variation vector at each node is generated by making finite difference to the shape with respect to each design parameter, and is used for the computation of sensitivity formula. Several numerical examples are taken to show the advantage of the method, in which the accuracy of the sensitivity is evaluated. The results are found excellent even by employing a simple linear function for the design velocity evaluation. Shape optimization is carried out for the geometric design of an archgrid and tilted bridge, which is to minimize maximum stress over the structure while maintaining constant weight. In conclusion, the proposed formulation is a useful and easy tool in finding optimum shape in a variety of the spatial frame structures.