Thermomechanical analysis of elastoplastic medium in sliding contact with fractal surface

The effects of mechanical and thermal surface loadings on deformation of elastic–plastic semi-infinite medium were analyzed simultaneously by using the finite element method. Rigid rough surface of a magnetic head and smooth surface of an elastic–plastic hard disk were chosen to perform a comprehensive thermo-elastic–plastic contact analysis at the head–disk interface (HDI). A two-dimensional finite element model of a rigid rough surface characterized by fractal geometry sliding over an elastic–plastic medium was then developed. The evolution of deformation in the semi-infinite medium due to thermomechanical surface loading is interpreted in terms of temperature, von Mises equivalent stress, and equivalent plastic strain. In addition to this, the effects of friction coefficient, sliding, and interference distance on deformation behavior were also analyzed. It is shown that frictional heating increases not only the contact area but also the contact pressure and stresses.     

Pulse loading shape effects on pressure–impulse diagram of an elastic–plastic, single-degree-of-freedom structural model

Pulse loading shape effects on the dynamical response of an elastic–plastic, single-degree-of-freedom (SDOF) structural model are studied in the present paper. Two dimensionless parameters are introduced to classify the studied problem into elastic, elastic–plastic and rigid–plastic structural responses. When structural damage is controlled by maximum structural deflection, a characteristic curve in loading parameter space can be used to define an isodamage curve, i.e., pressure–impulse diagram, which is generally loading-shape-dependent. Dimensionless loading parameters, termed as effective loading parameters, are introduced in the present paper to give unique loading-shape-independent pressure–impulse diagram for each response category.     

A surface flow line model of a scratching tip: apparent and true local friction coefficients

The apparent friction coefficient is the ratio between the tangential force and the normal load applied to a moving scratching tip. It includes a so-called “true local” friction coefficient, which is the scission at the interface between the tip and the surface being scratched, and a “geometrical” friction coefficient, which is the plough effect due to the wave front created ahead of the moving tip and depends on the shape of the tip. Like in any mechanical test, three basic types of behaviour of the material at the interface are observed: purely elastic, elastic–plastic and fully plastic. As is usual in polymers, the material behaviour is time and temperature dependent and may exhibit strain hardening. A surface flow line model is developed here to deduce the geometrical and true friction coefficients at the interface between a moving scratching tip and a surface from the apparent friction coefficient. Using this model, several situations may be simulated to predict the influence of the geometry of the tip on the scratch resistance of the material.     

Anomalous dynamic elastic-plastic response of a galerkin beam model

A study of the anomalous motion of an elastic—plastic beam under short pulse loading is presented. The geometric nonlinearity due to axial end constraints is taken into account. We apply the Galerkin method to the governing partial differential equation of the transverse motion to obtain a general model of n degrees of freedom (nDoF). The results of elastic—plastic deformation analysis and dynamic response for the 2DoF model of a pin-ended beam are presented. The regular and irregular motions of the 2DoF model for the pin-ended beam are examined by various methods including time history, phase diagram, Lyapunov characteristic exponent and power spectral density.     

An elastic–plastic stress analysis steel-reinforced thermoplastic composite cantilever beam

In the present study an analytical elastic–plastic stress analysis is carried out for a low-density homogeneous polyethylene thermoplastic cantilever beam reinforced by steel fibers. The beam is loaded by a constant single force at its free end. The expansion of the region and the residual stress component of σ_x are determined for 0°, 30°, 45°, 60° and 90° orientation angles. Yielding begins for 0° and 90° orientation angles at the upper and lower surfaces of the beam at the same distances from the free end. However, it starts first at the upper surface for 30° and 45° orientation angles. The elastic–plastic analysis is carried out for both the plastic region which spreads only at the upper surface and the plastic region which spreads at the upper and lower surfaces together. The residual stress components of σ_x and τ_xy are also determined. The intensity of the residual stress component is maximum at the upper and lower surfaces of the beam, but the residual stress component of τ_xy is maximum on or around the x-axis. The beam can be strengthened by using the residual stresses. The distance between the plastically collapsed point and the free end is calculated for the same load in the beam for 0°, 30°, 45°, 60° and 90° orientation angles.     

Coupled buckling and plastic instability for tube hydroforming

In this paper, the hydroforming limit of isotropic tubes subjected to internal hydraulic pressure and independent axial load is discussed.Swift's criterion is often used in this case for the prediction of diffuse plastic instability. Here, we first highlight the existence of two different Swift's criteria (for sheets and for tubes).Then, we recall that these types of approaches do not take into account buckling induced by axial loading. In fact, buckling may obviously occur before plastic instability; consequently, Swift's criteria must not be used alone to predict instability in the case of tube hydroforming.Numerical simulation was used to confirm these points and to analyse both the buckling and striction phenomena together. The two types of instability must be treated together in a reasonable approach to the hydroforming process.In this paper, the material verifies a “J2-flow” constitutive rate constitutive law. Jaumann's derivative was chosen and the Prandtl–Reuss equations with von Mises’ yield criterion and the associated flow rule were used. Isotropic hardening was taken into account.     

Plastic deformation and contact area of an elastic–plastic contact of ellipsoid bodies after unloading

This paper presents theoretical and experimental results of the residual or plastic deformation and the plastic contact area of an elastic–plastic contact of ellipsoid bodies after unloading. There are three regime responses of the deformation and contact area: elastic, elastic–plastic and fully plastic. Experimental investigation is presented in order to validate the proposed model. A new technique is introduced to measure the plastic deformation and plastic contact area. Very good correlation is found between the theoretical prediction and the experimental results.     

Numerical study of the flange earring of deep-drawing sheets with stronger anisotropy

A Barlat–Lian anisotropy yield function is introduced into a quasi-flow corner theory of elastic–plastic finite deformation and the elastic–plastic large deformation finite element formulation based on the principle of virtual velocity and the discrete Kirchhoff triangle plate shell element model. The focus of the present researches is on the numerical simulation of the flange earring of deep-drawing process of circular sheets with stronger anisotropy, based on which, the schemes for controlling the flange earring are proposed.     

On the measurement of friction coefficient utilizing the ring compression test

The main objective of this research was to investigate whether generalized friction calibration curves, as recommended in the literature for use with ring compression tests, are applicable to all types of materials and test conditions. Specifically, the effects of material properties, strain-rate sensitivity, and “barreling” on the behavior of friction calibration curves were investigated. To this end, a series of ring compression tests were conducted in order to determine the magnitude of the friction coefficient, μ, as well as the corresponding calibration curves for two types of modeling materials, white and black Plasticine. The experiments were first conducted using the Physical Modeling Technique (PMT) and then simulated via an elastic–plastic finite element code (ABAQUS). In contrast to the results available in the literature, where the same friction calibration curves are recommended for all types of materials and test conditions, the results of this investigation showed that friction calibration curves are indeed affected by the material properties and test conditions and every material possesses its own distinctive friction calibration curve.     

Flywheel energy storage—II: Magnetically suspended superflywheel

This article, the second of a two part paper, describes the general design requirements for a flywheel energy storage system. A new superflywheel energy storage system, using a spokeless, magnetically suspended, composite material pierced disk rotor is proposed. The new system is configured around a permanent magnet (“flux biased”) magnetic suspension system with active control in the radial direction and passive control in the axial direction. The storage ring is used as a moving rotor and electronic commutation of stationary armature coils is proposed. There is no mechanical contact with the rotating ring and long life and low run down losses are projected. A discussion of major components for a 10 kwh system is presented.     

Simulation of rubber isostatic pressing and shape optimization of rubber mold

An elastic–plastic FEM is developed to study the fundamental features of rubber isostatic pressing of powder with a view to performing a net-shape process. Experiments are also performed. It is thus shown that the properties of the rubber, in particular Poisson's ratio, give a great influence on the shape of compact and density distribution. The thickness of the rubber mold is also influential; the thicker the mold, the closer is the shape of compact to the cavity shape. An optimization scheme is also developed, so that one is able to determine the cavity shape that gives a desired shape of compact.     

Numerical simulation of the localization behavior of hydrostatic-stress-sensitive metals

The present paper deals with the numerical simulation of the elastic–plastic deformation and localization behavior of solids which are plastically dilatant and sensitive to hydrostatic stresses. The model is based on a generalized macroscopic theory taking into account macroscopic as well as microscopic experimental data obtained from tests with iron-based metals. It shows that hydrostatic components may have a significant effect on the onset of localization and the associated deformation modes, and that they generally lead to a notable decrease in ductility. The continuum formulation relies on a generalized I₁–J₂–J₃ yield criterion to describe the effect of the hydrostatic stress on the plastic flow properties of metals. In contrast to classical theories of metal plasticity, the evolution of the plastic part of the strain rate tensor is determined by a non-associated flow rule based on a plastic potential function which is expressed in terms of stress invariants and kinematic parameters. Numerical simulations of the elastic–plastic deformation behavior of hydrostatic-stress-sensitive metals show the physical effects of the model parameters and also demonstrate the efficiency of the formulation. Their results are in excellent agreement with available experimental data. A variety of large-strain elastic–plastic problems involving pronounced localizations is presented, and the influence of various model parameters on the deformation and localization behavior of hydrostatic-stress-sensitive metals is discussed.     

Theoretical investigation on the dimple occurrence in the thermal EHL of simple sliding steel–glass circular contacts

The conical depression (surface dimple) phenomena observed by Kaneta et al. (Kaneta M, Nishikawa H, Kameishi K, Sakai T. Effects of elastic moduli of contact surfaces in elastohydrodynamic lubrication. ASME J. Tribol. 1992;114:75–80; Kaneta M, Nishikawa H, Kanada T, Matsuda K. Abnormal phenomena appearing in EHL contacts. ASME J. Tribol. 1996;118:886–892.) in optical interferometry experiments are simulated numerically by a complete solution to the simple sliding circular contact thermal elastohydrodynamic lubrication (TEHL) problem. Good agreement is displayed between the theoretical and experimental results. This agreement is explained by the “temperature–viscosity wedge” mechanism, which was first proposed by Cameron (Cameron A. Hydrodynamic lubrication of rotating discs in pure sliding, a new type of oil film formation. J. Inst. Petrol 1951;37:471.). Effects of the viscosity–pressure coefficient, the ambient viscosity, and the entrainment velocity on the behavior of the surface dimples are discussed.     

Dynamic asymmetrical instability of elastic–plastic beams

Dynamic instability of elastic–plastic beam is investigated by employing a three-degree-of-freedom (3-DoF) beam model. Especially, asymmetrical instability induced by symmetrical load is discussed. The asymmetrical instability is considered as a second-order buckling mode. Four types of perturbations, i.e., geometrical misalignment, material property mismatch, unsymmetry of applied load and disturbance of boundary conditions, are introduced to activate the asymmetrical responses. The asymmetrical response is characterized by a modal participation factor α₂ which corresponds to an asymmetrical mode shape. Phase plane trajectories and Poincaré map are used to illustrate the chaotic characteristics of the beam response. Results show that if the perturbations are small enough, the perturbation type has negligible influence on the critical load for the occurrence of the asymmetrical instability, which implies that the asymmetrical instability is an intrinsic feature of the beam system. However, with the increase of the magnitude of the perturbations, the influence of the asymmetrical vibration is expanded to a large extension of loading parameter.     

Stress distributions in energy generating two-layer tubes subjected to free and radially constrained boundary conditions

Analytical solutions are obtained for thermally induced axisymmetric elastic and elastic–plastic deformations in heat generating composite tubes having a free inner and a radially constrained outer boundaries. Tresca's yield condition and its associated flow rule are used to determine the elastic–plastic response of the assembly. Depending on the physical properties of the materials used, eight different plastic regions with different mathematical forms of the yield condition may occur in the assembly. The closed form solutions for these plastic regions are obtained by assuming linearly hardening material behavior. Various numerical examples are handled using thermal and mechanical properties of real engineering materials.     

New cylinder liner surfaces for low oil consumption

Anticipated emission legislation and reduced fuel consumption are the main driving forces when developing new engines. Optimization of the active surfaces in the piston system is one possible way to meet the above demands. In this study the effects of surface topography and texture direction of the ring/liner contact on oil film thickness and friction were simulated and experimentally tested. “Low wear” results from the experimental wear tests with “glide honed” smooth liner surfaces supported the “low friction” simulation results. In addition a new wear volume sensitive surface roughness parameter, R_ktot, based on the Abbot–Firestone bearing area curve was introduced.     

Pressure between elastic bodies having a slender area of contact and arbitrary profiles

A numerical method is presented for calculating the pressure distribution and contact area shape between two elastic bodies of arbitrary profile which make contact over a slender contact area, i.e. where the relative curvature of the two profiles is much smaller in the longitudinal direction than in the transverse. The pressure distribution is assumed to be piecewise-linear in the longitudinal direction and semi-elliptical in the transverse. No a priori relationship is assumed between the shape of the contact area and the longitudinal variation in pressure; they are found simultaneously from dual integral equations for the compatibility of (a) the normal displacement and (b) the transverse curvature along the longitudinal axis of the contact zone.In cases where the profiles of the contacting bodies are smooth and continuous up to, and beyond, the ends of the contact area, the method gives a very reliable measure of the contact pressure distribution. Where discontinuities in profile are present, at roller ends for example, stress singularities are to be expected and like any numerical method, only approximate values of the stress concentration can be found. In the cases studied, the concentration of pressure associated with a “sharp” edge of contact is found to be very local.The method has been applied to both cylindrical and variously “crowned” rollers, also to a ball “over-riding” the edge of a closely conforming groove.     

Inelastic buckling of columns : The effect of imperfections

In the design of columns of mild steel (idealized as an elastic-perfectly plastic material) it is usual to take account of the effect of possible initial crookedness by means of a “Perry” formula. In contrast, the design of columns of aluminium alloys (and other materials which cannot reasonably be idealized as perfectly plastic) is usually made by means of the “tangent modulus” formula, which is strictly relevant only to initially perfect columns. The paper proposes a way of supplementing this formula for initially imperfect columns, and a simple graphical procedure is devised to generate a family of “column curves” for different degrees of imperfection.It turns out that although the “column curve” based on the tangent-modulus formula is sensitive to the precise shape of the rising stress-strain curve, the curves for the imperfect columns are insensitive to this shape, except for stocky columns. This suggests, paradoxically, a possible design approach using a Perry formula for columns made of aluminium alloys.     

Damage evolution in creep bulging of thin sheet metal

The creeping motion of thin sheet metal, damaged by artificial cavities is observed in bulging tests and simulated ‘semi’-analytically. The sheet metal satisfies Norton’s Law for secondary creep and is subjected to a bi-directional stretch. The stretch is produced by creep bulging through elliptical dies with the virtue of sustaining nearly uniform background stress ratio for each aspect ratio of the die axes. In order to reach large deformations with significant shape evolution of the cavities, the tests were conducted at superplastic conditions. The sheet is double layered (only one layer is cavitated) made of Tin–Lead (50–50 Pb–Sn). The measured damage growth is compared to an approximate simulation. The simulation of the damage evolution, throughout its time history, makes repeated use of a so-called “Green-function solution” for the motion of a single isolated cavity in an infinite viscoplastic continuum. The solution is modified from Muskhelishvili’s elastic solution by replacing the elastic shear modulus by a “viscous-like” variable (“plastic shear modulus”) which depends (non-linearly) on the evolved average strain-rate. Similarly, the stresses in the ligaments between cavities were averaged to approximate the local stress concentrations. Due emphasis is given to the rotation of each elliptical cavity, beside its expansion (contraction) and elongation.     

Substepping algorithms with stress correction for the simulation of sheet metal forming process

The finite element analysis of the sheet metal forming process involves various nonlinearities. To predict accurately the final geometry of the sheet blank and the distribution of strain and stress and control various forming defects, such as thinning, wrinkling and springback, etc., the accurate integration of the constitutive laws over the strain path is essential. Our objective in this paper is to develop an effective and accurate stress integration scheme for the analysis of three-dimensional sheet metal forming problems. The proposed algorithm is based on the explicit “substepping” schemes incorporating with the stress correction scheme. The proposed algorithms have been implemented into ABAQUS/Explicit via User Material Subroutine (VUMAT) interface platform. The algorithms are then employed to analyze a typical deep-cup drawing process and the accuracy of these algorithms has been compared with the implicit “return” algorithm and explicit forward algorithm. The results indicate that the explicit schemes with local truncation error control, together with a subsequent check of the consistency conditions, can achieve the same or even better level of accuracy as “return” algorithm does for integrating large plastic problems like sheet metal forming process.