• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2214-2219
• /
• 2008

The thermal contact resistance(TCR) of nanosized contact spots has been investigated through a multiscale analysis which considers the resolution of surface topography. A numerical simulation is performed on the finite element model of rough surfaces. Especially, as the contact size decreases below the phonon mean free path, the size dependent thermal conductivity is considered to calculate the TCR. In our earlier model which follows an elastic material, the TCR increases without limits as the number of nanosized contact spots increases in the process of scale variation. However, the elastoplastic contact induces a finite limit of TCR as the scale varies. The results are explained through the plastic behavior of the two contacting models. Furthermore, the effect of air conduction in nanoscale is also investigated.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.1
• /
• pp.1-6
• /
• 2017

This paper presents development of the nine-node co-rotational(CR) planar element applicable for elastoplastic and contact analysis. The CR formulation is one of the efficient geometrically nonlinear formulations. It is based on the assumptions of small strain and large displacement. Further, it is extended to both elastoplastic analysis and contact analysis in this paper. For accurate plastic analysis, nine-node quadrilateral element, which can provide accurate stress prediction, is chosen. Bi-linear hardening rule based on Newton- Raphson return-mapping is employed. Also, Lagrange multiplier is used in order for constraints regarding the contact analysis. The present development is validated via the time transient problems. The present results are compared with those obtained by the other existing software.

• Tribology and Lubricants
• /
• v.40 no.1
• /
• pp.8-16
• /
• 2024

This study conducts numerical modeling and eigen-analysis of a rod-fastened rotor, which is mainly used in aircraft gas turbine engines in which multiple disks are in contact through curvic coupling. Nayak's theory is adopted to calculate surface parameters measured from the tooth profile of the curvic coupling gear. Surface parameters are important design parameters for predicting the stiffness between contact surfaces. Based on the calculated surface parameters, elastoplastic contact analysis is performed according to the interference between two surfaces based on the Greenwood-Williamson model. The equivalent bending stiffness is predicted based on the shape and elastoplastic contact stiffness of the curvic coupling. An equation of motion of the rod-fastened rotor, including the bending stiffness of the curvic coupling, is developed. Methods for applying the bending stiffness of a curvic coupling to the equation of motion and for modeling the equation of motion of a rotor that includes both inner and outer rotors are introduced. Rotordynamic analysis is performed through one-dimensional finite element analysis, and each element is modeled based on Timoshenko beam theory. Changes in bending stiffness and the resultant critical speed change in accordance with the rod fastening force are predicted, and the corresponding mode shapes are analyzed.

• Transactions of Materials Processing
• /
• v.8 no.3
• /
• pp.306-312
• /
• 1999

In the recent sheet metal forming simulations, it increases to adopt the dynamic explicit method for an effective computation and the elastoplastic formulation for stress recovery. It is inevitable in the dynamic explicit method that some noises occur, which sometimes partly spoil results of simulations. This phenomenon becomes severer when complicate contact conditions are included in simulations. In commercial dynamic codes, the concept of contact damping is introduced. However, the formulation process of it is not revealed well. In this paper, a contact damping method is formulated in order for effectively suppressing noises occurring due to complicated contact conditions. This is checked by analyzing a simple sheet metal stamping process (U-draw bending). From the computational results, it is shown that the contact damping can effectively control the noises due to contacts, especially when considering the sheet thickness, and help to develop more reliable internal stress states, which result in more realistic shapes after springbank.

• Proceedings of the KSME Conference
• /
• 2008.11a
• /
• pp.144-149
• /
• 2008

A transient finite element simulation is developed for the two-dimensional stationary elastoplastic layer between sliding layers, to investigate thermoelastoplastic instability(TEPI) due to frictional heating in the material. The analysis will show some differences between the case of thermoelastic instability and TEPI, especially according to the contact pressure above yield stress. A transient behavior of contact pressure is captured to explain the behavior of thermoplasticity of contact with different sliding velocity. The instability of contact pressure in the long range of braking time will be explored to understand the generation mechanism of hot spots.

• The Journal of Engineering Geology
• /
• v.8 no.3
• /
• pp.215-224
• /
• 1998

In this study, total strain was measured by LVDTs and local strains on the surface of specimens were measured by strain gauges. And axi-symmetrically elastoplastic FEM analyses was carried out for cylindrical specimens. Considering the influence of the restraint induced by the loading platen, in the case of H/D=1, the strain distribution on the side of a specimen is obviously affected by the condition of platen contact. Furthermore, it is clear that the larger H/D ratio becomes, the smaller the influence to the strain distribution is. For the smooth contact condition, the strain on the side is not influenced by the stiffness of the specimen, the shape and the scale effect, the strain distribution coincides with the nominal total strain. Whereas, in the case of rough contact condition, the strain distribution is remarkably affected. It is made clear that strain responses of hard rock specimens may more sensitive than these of soft rock specimens as a results of interaction between loading platens and specimen and the uniaxial strength of specimens may strongly depends on this interaction and stress-strain relation is affected by the contact condition.

• Journal of Advanced Marine Engineering and Technology
• /
• v.34 no.7
• /
• pp.943-950
• /
• 2010

An efficient meshfree method based on a stabilized conforming nodal integration method is developed for elastoplastic contact analysis of metal forming processes. In this approach, strain smoothing stabilization is introduced to eliminate spatial instability in Galerkin meshfree methods when the weak form is integrated by a nodal integration. The gradient matrix associated with strain smoothing satisfies the integration constraint for linear exactness in the Galerkin approximation. Strain smoothing formulation and numerical procedures for path-dependent problems are introduced. Applications of metal forming analysis are presented, from which the computational efficiency has been improved significantly without loss of accuracy.

• Journal of the Korean Society for Precision Engineering
• /
• v.20 no.6
• /
• pp.37-46
• /
• 2003

This paper is to predict quality deterioration resulting from a caulking process of yoke which is a part of automotive steering system. The caluking is a plastic deformation process involving such as impact of high speed tool, contacts between part and fixtures and strain rate sensitivity of the part material. Elaborate application of finite element method is neccesary to calculate changes of part dimensions because they fall into a level of tolerances. Simple work hardening and strain rate sensitive model is proposed fur the material and applied for the simulation by using Abaqus which is able to cater for elastoplastic rate sensitive material and contacts. Numerical results of test models that represent tensile bar and tensile plate are compared with material data inputs. Dimensional changes for the yoke are calculated from simulations and compared to the mesurements and they show good agreement. The method presented here with the material model proved to be valuable to assess quality deterioration for similar metal forming processes.

• Advances in concrete construction
• /
• v.5 no.3
• /
• pp.221-240
• /
• 2017

The modeling of the mechanical behavior of quasi-brittle materials is still a challenge task, mainly in failure processes when fracture and plasticity phenomena become important actors in dissipative processes which occur in materials like concrete, as instance. Many homogenization-based approaches have been proposed to deal with heterogeneous materials in the last years. In this context, a computational homogenization modeling for concrete is presented in this work using the concept of Representative Volume Element (RVE). The material is considered as a three-phase material consisting of interface zone (ITZ), matrix and inclusions-each constituent modeled by an independent constitutive model. The Representative Volume Element (RVE) consists of inclusions idealized as circular shapes symmetrically and nonsymmetrically placed into the specimen. The interface zone is modeled by means of cohesive contact finite elements. The inclusion is modeled as linear elastic and matrix region is considered as elastoplastic material. A set of examples is presented in order to show the potentialities and limitations of the proposed modeling. The consideration of the fracture processes in the ITZ is fundamental to capture complex macroscopic characteristics of the material using simple constitutive models at mesoscopic level.

• Steel and Composite Structures
• /
• v.13 no.4
• /
• pp.397-408
• /
• 2012

This paper describes the refined 3-D finite element (FE) modeling of composite frames composed of concrete-filled steel tubular (CFST) columns and steel-concrete composite beams based on the test to get a better understanding of the seismic behavior of the steel-concrete composite frames. A number of material nonlinearities and contact nonlinearities, as well as geometry nonlinearities, were taken into account. The elastoplastic behavior, as well as fracture and post-fracture behavior, of the FE models were in good agreement with those of the specimens. Besides, the beam and panel zone deformation of the analysis models fitted well with the corresponding deformation of the specimens. Parametric studies were conducted based on the refined finite elememt (FE) model. The analyzed parameters include slab width, slab thickness, shear connection degree and axial force ratio. The influences of these parameters, together with the presence of transverse beam, on the seismic behavior of the composite frame were studied. And some advices for the corresponding seismic design provisions of composite structures were proposed.