• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.7 s.94
• /
• pp.1658-1667
• /
• 1993

In this research, a numerical algorithm has been developed, which can be applied to the large deformation and large displacement contact problems between two deformable bodies. The contact conditions expressed in terms of the rate of angle change have been proposed considering the change in geometric shape and rate of contact force. A set of linear simultaneous equations is constructed by adding the geometric shape change and contact conditions to the original stiffness matrix. A new method to determine time increment has been proposed based on Euler method, in which the condition to prevent the contact bodies from penetrating and overrunning each other has been taken into consideration. Practical application to contact problem is extrusion in which bodies are sliding along the contact boundary.

• Tribology and Lubricants
• /
• v.20 no.5
• /
• pp.284-291
• /
• 2004

We present a three-dimensional average flow model considering elastic deformation of pad asperities for chemical mechanical planarization. To consider the contact deformation of pad asperities in the calculation of the flow factor, three-dimensional contact analysis of a semi-infinite solid based on the use of influence functions is conducted from computer generated three dimensional roughness data. The average Reynolds equation and the boundary condition of both force and momentum balance are used to investigate the effect of pad roughness and external pressure conditions on film thickness and wafer position angle.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1994.03a
• /
• pp.47-54
• /
• 1994

The study is concerned with the coupled analysis of heat transfer and thermoelastoplastic deformation. The thermoelastoplastic model is very useful for the analysis of residual stress and the analysis of thermal stress as well as the analysis of metal forming. Heat of deformation, phase transformation and contact heat transfer boundary are considered. The contact heat transfer boundary is treated by the interpolation of shape function. The analysis of deformation and the analysis of heat transfer are carried out for the could upsetting and the hot rolling. The computed results are found to be in good agreement with the experimental results.

• Structural Engineering and Mechanics
• /
• v.6 no.1
• /
• pp.63-76
• /
• 1998

This paper presents dynamic responses of structures with sliding base which limits the translation of external loads from ground excitation. A discrete element model based on the discontinuous deformation analysis method is proposed to study this sliding boundary problem. The sliding base is simulated using sets of fictitious contact springs along the sliding interface. The set of contact spring is to translate friction force from ground to superstructure. Validity of the proposed model is examined by the closed-form solutions of an idealized mass-spring structural model subjected to harmonic ground excitation. This model is also applied to a problem of a three-story structural model subjected to the ground excitation of 1940 El Centro earthquake. Analyses of both sliding-base and fixed-base conditions are performed as comparisons. This study shows that using this model can simulate the dynamic response of a sliding structure with frictional cut-off quite accurately. Results reveal that lowering the frictional coefficient of the sliding joint will reduce the peak responses. The structure responses in little deformation, but it displaces at the end of excitation.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.237-255
• /
• 2013

In this paper, a meshfree-enriched finite element method (ME-FEM) is introduced for the large deformation analysis of nonlinear path-dependent problems involving contact. In linear ME-FEM, the element formulation is established by introducing a meshfree convex approximation into the linear triangular element in 2D and linear tetrahedron element in 3D along with an enriched meshfree node. In nonlinear formulation, the area-weighted smoothing scheme for deformation gradient is then developed in conjunction with the meshfree-enriched element interpolation functions to yield a discrete divergence-free property at the integration points, which is essential to enhance the stress calculation in the stage of plastic deformation. A modified variational formulation using the smoothed deformation gradient is developed for path-dependent material analysis. In the industrial metal forming problems, the mortar contact algorithm is implemented in the explicit formulation. Since the meshfree-enriched element shape functions are constructed using the meshfree convex approximation, they pose the desired Kronecker-delta property at the element edge thus requires no special treatments in the enforcement of essential boundary condition as well as the contact conditions. As a result, this approach can be easily incorporated into a conventional displacement-based finite element code. Two elasto-plastic problems are studied and the numerical results indicated that ME-FEM is capable of delivering a volumetric locking-free and pressure oscillation-free solutions for the large deformation problems in metal forming analysis.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1992.03a
• /
• pp.11-26
• /
• 1992

In this research, a numerical algorithm has been developed, which can be applied to the large deformation and large displacement contact problems between angle change have been proposed considering the change in geometric shape and rate of contact force. A set of linear simultaneous equations is constructed by adding the geometric shape change and contact conditions to the original stiffness matrix. A new method to determine time increment has been proposed based on Euler method, in which the condition to prevent the contact bodies from penetrating and overrunning each other has been taken into consideration. Practical application to contact problem is extrusion in which bodies are sliding along the contact boundary.

• Structural Engineering and Mechanics
• /
• v.25 no.5
• /
• pp.501-518
• /
• 2007

In this work, a new approach is developed for dynamic analysis of a composite beam with an interply crack, based on finite element solution of partial differential equations with the use of the COMSOL Multiphysics package, allowing for fast and simple change of geometric characteristics of the delaminated area. The use of COMSOL Multiphysics package facilitates automatic mesh generation, which is needed if the problem has to be solved many times with different crack lengths. In the model, a physically impossible interpenetration of the crack faces is prevented by imposing a special constraint, leading to taking account of a force of contact interaction of the crack faces and to nonlinearity of the formulated boundary value problem. The model is based on the first order shear deformation theory, i.e., the longitudinal displacement is assumed to vary linearly through the beam's thickness. The shear deformation and rotary inertia terms are included into the formulation, to achieve better accuracy. Nonlinear partial differential equations of motion with boundary conditions are developed and written in the format acceptable by the COMSOL Multiphysics package. An example problem of a clamped-free beam with a piezoelectric actuator is considered, and its finite element solution is obtained. A noticeable difference of forced vibrations of the delaminated and undelaminated beams due to the contact interaction of the crack's faces is predicted by the developed model.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.2
• /
• pp.153-161
• /
• 2009

A domain/boundary decomposition method is applied for efficient analyses of thermo-elasto-viscoplastic damage and contact problems under the assumption of infinitesimal deformation. For the decomposition of a whole domain and contact boundaries, all the equality constraints on the interface and contact interfaces are restated with simple penalty functional. Therefore, the non-linearity of the problem is localized within finite element matrices in a few subdomains and on contact interfaces. By setting up suitable solution algorithms, the computational efficiency can be improved considerably. The general tendency of the computational efficiency is illustrated with some numerical experiments.

• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 2004.11a
• /
• pp.331-338
• /
• 2004

We present a three-dimensional average flow model considering elastic deformation of pad asperities for chemical mechanical planarization. To consider the contact deformation of pad asperities in the calculation of the flow factor, three-dimensional contact analysis of a semi-infinite solid based on the use of influence functions is conducted from computer generated three dimensional roughness data. The average Reynolds equation and the boundary condition of both force and momentum balance are used to investigate the effect of pad roughness and external pressure conditions on film thickness and wafer position angle.

• Transactions of Materials Processing
• /
• v.12 no.1
• /
• pp.34-42
• /
• 2003

This paper presents a rigid-plastic finite element method to handle the frictional contact problem between two deformable bodies experiencing large deformation. The variational formulation combined with incremental quasi-static model is employed for treating the contact boundary condition. The frictional behavior of the model obeys Coulomb's law of friction. The proposed contact algorithms are classified into two categories, one for searching contacting nodes and the other for calculating contact forces at the contact surface. A slave node and master contact segment are defined using the geometric condition of finite elements on the contact interface. The penalty parameter is used to limit the penetration between contacting bodies, and the finite elements are coupled with contact boundary elements.us gates and cavity thicknesses. Through this study we have observed that the jetting is related to the die swell of material. This means that the jotting is strongly affected by the elastic flow property rather than the viscous flow property in viscoelastic characteristics of molten polymer. Different resins have different elastic properties, and elastic flow behavior depends on the shear rate of flow, i.e. injection speed. Large die swell would eliminate jetting however, the retardation of die swell would stimulate jetting. In the point of mole design, reducing the thickness ratio of cavity to gate can reduce or eliminate jetting regardless of amount of elasticity of polymer melt.