• Interaction and multiscale mechanics
• /
• 제6권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.

• 한국공작기계학회논문집
• /
• 제12권5호
• /
• pp.24-33
• /
• 2003

This study proposed a technique for inverse problem, linear approximation of contact position and loading in single and double meshing of transmission contact element using 2-dimension model considered near the tooth by root stress. Determination of root stress is carried out far the gear tooth by finite element method and boundary element method. Boundary element discretization near contact point is carefully performed to keep high computational accuracy. The predicted results of boundary element method are good accordance with that of finite element method.

• Tribology and Lubricants
• /
• 제36권5호
• /
• pp.262-266
• /
• 2020

Herein, we present a numerical investigation of wear analysis of sliding systems with a constant speed subjected to Archard's wear law. For this investigation, we compared two methods: eigen-wear analysis and adaptive meshing technique. The eigen-wear analysis is advantageous to predict the evolution of contact pressure due to wear using the initial contact pressure and contact stiffness. The adaptive meshing technique in finite element analysis is employed to obtain transient wear behavior, which needs significant computational resources. From the eigen-wear analysis, we can determine the appropriate element size required for finite element analysis and the time increment required for wear evolution by a dimensionless variable above a certain value. Since the prediction of wear depends on the maximum contact pressure, the finite element model should have a reasonable representation of the maximum contact pressure. The maximum contact pressure and wear amount according to this dimensionless variable shows that the number of fine meshes in the contact area contributes more to the accuracy of the wear analysis, and the time increment is less sensitive when the number of contact nodes is significantly larger. The results derived from a two-dimensional wear model can be applied to a three-dimensional wear model.

• Structural Engineering and Mechanics
• /
• 제54권4호
• /
• pp.607-622
• /
• 2015

This paper presents a comparative study of analytical method and finite element method (FEM) for analysis of a continuous contact problem. The problem consists of two elastic layers loaded by means of a rigid circular punch and resting on semi-infinite plane. It is assumed that all surfaces are frictionless and only compressive normal tractions can be transmitted through the contact areas. Firstly, analytical solution of the problem is obtained by using theory of elasticity and integral transform techniques. Then, finite element model of the problem is constituted using ANSYS software and the two dimensional analysis of the problem is carried out. The contact stresses under rigid circular punch, the contact areas, normal stresses along the axis of symmetry are obtained for both solutions. The results show that contact stresses and the normal stresses obtained from finite element method (FEM) provide boundary conditions of the problem as well as analytical results. Also, the contact areas obtained from finite element method are very close to results obtained from analytical method; disagree by 0.03-1.61%. Finally, it can be said that there is a good agreement between two methods.

• 한국주조공학회지
• /
• 제40권6호
• /
• pp.146-150
• /
• 2020

In general, hot metal castings contract and molds expand during the cooling step of a casting process. Therefore, it is important to consider both the casting and mold at the same time in a casting process analysis. For a more accurate analysis that includes the contact characteristics, matching each node of the casting and mold in the contact area is recommended. However, it is very difficult to match the nodes of the casting and the mold when generating elements due to the geometric problem of CAD model data. The present study proposes a mesh generation technique that considers mechanical contact between the casting and the mold in a casting analysis (finite element analysis). The technique focuses on the fact that the mold surrounds the casting. After generating the 3D elements for the casting, the surface elements of the casting in contact with the mold are projected inside the mold to create contact elements that coincide with the contact surface of the casting. It was confirmed that high-quality contact element information and a 3D element net can be automatically generated by the method proposed in this study.

• Structural Engineering and Mechanics
• /
• 제64권4호
• /
• pp.505-513
• /
• 2017

In this paper, a receding contact problem for an elastic layer resting on a half plane is considered. The layer is pressed by two rectangular stamps placed symmetrically. It is assumed that the contact surfaces are frictionless and only compressive traction can be transmitted through the contact surfaces. In addition the effect of body forces is neglected. Firstly, the problem is solved analytically based on theory of elasticity. In this solution, the problem is reduced into a system of singular integral equations in which half contact length and contact pressures are unknowns using boundary conditions and integral transform techniques. This system is solved numerically using Gauss-Jacobi integral formulation. Secondly, two dimensional finite element analysis of the problem is carried out using ANSYS. The dimensionless quantities for the contact length and the contact pressures are calculated under various stamp size, stamp position and material properties using both solutions. The analytic results are verified by comparison with finite element results.

• Structural Engineering and Mechanics
• /
• 제76권3호
• /
• pp.325-336
• /
• 2020

In this study, the continuous and discontinuous contact problems of functionally graded (FG) layer resting on a rigid foundation were considered. The top of the FG layer was loaded by a distributed load. It was assumed that the shear modulus and the density of the layer varied according to exponential functions along the depth whereas the the Poisson ratio remained constant. The problem first was solved analytically and the results were verified with the ones obtained from finite element (FE) solution. In analytical solution, the stress and displacement components for FG layer were obtained by the help of Fourier integral transform. Critical load expression and integral equation for continuous and discontinuous contact, respectively, using corresponding boundary conditions in each case. The finite element solution of the problem was carried out using ANSYS software program. In continuous contact case, initial separation distance and contact stresses along the contact surface between the FG layer and the rigid foundation were examined. Separation distances and contact stresses were obtained in case of discontinuous contact. The effect of material properties and loading were investigated using both analytical and FE solutions. It was shown that obtained results were compatible with each other.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 추계학술대회
• /
• pp.1286-1290
• /
• 2003

This paper presents a newly implemented parallel finite element procedure for contact-impact problems. Three sub-algorithms are includes in the proposed parallel contact-impact procedure, such as a parallel Belytschko-Lin-Tsay (BLT) shell element generation, a parallel explicit time integration scheme, and a parallel contact search algorithm based on the master slave slide-line algorithm. The underlying focus of the algorithms is on its effectiveness and efficiency for inclusion in future finite element systems on parallel computers. Throughout this research, a prototype code, named GT-PARADYN, is developed on the IBM SP2, a distributed-memory computer. Some numerical examples are provided to demonstrate the timing results of the procedure, discussing the accuracy and efficiency of the code.

• 산업기술연구
• /
• 제28권A호
• /
• pp.81-88
• /
• 2008

Generally, Hertz theory is used to analyze the contact problem of two bodies. It is simple derivation of solution in the contact part. And calculation time is short Moreover, it can mean well that many wear occurs relatively. However, material property becomes plastic deformation when large perpendicular pressure acts on a small contact surface product. In this case, Hertz theory is inapplicable. Therefore this thesis carried the finite element analysis in consideration of material elasticitystrain and the shape of the geometric from contact point. And it compared with Hertz theory that change of the contact surface and contact pressure.

• 대한기계학회논문집A
• /
• 제24권5호
• /
• pp.1238-1245
• /
• 2000

One of the main causes of severe wear or crack initiation in wheel and rail is the contact stress due to wheel-rail contact. First, we obtain contact stress due to the rail mounting slope using the finite element method. Second, the shape design based on more reasonable contact stress analysis rather th~n a general Hertzian contact theory is investigated in order to reduce the contact stress. The optimum -design is performed using the simple 2-D finite element model and its results are verified by 311) finite element analysis.