Elasto-Plastic Coupled Temperature-Displacement Finite Element Analysis of Two-Dimensional Rolling-Sliding Contact With a Translating Heat Source

1991 ◽  
Vol 113 (1) ◽  
pp. 93-101 ◽  
Author(s):  
S. M. Kulkarni ◽  
C. A. Rubin ◽  
G. T. Hahn

The present paper, describes a transient translating elasto-plastic thermo-mechanical finite element model to study 2-D frictional rolling contact. Frictional two-dimensional contact is simulated by repeatedly translating a non-uniform thermo-mechanical distribution across the surface of an elasto-plastic half space. The half space is represented by a two dimensional finite element mesh with appropriate boundaries. Calculations are for an elastic-perfectly plastic material and the selected thermo-physical properties are assumed to be temperature independent. The paper presents temperature variations, stress and plastic strain distributions and deformations. Residual tensile stresses are observed. The magnitude and depth of these stresses depends on 1) the temperature gradients and 2) the magnitudes of the normal and tangential tractions.

1985 ◽  
Vol 52 (1) ◽  
pp. 67-74 ◽  
Author(s):  
V. Bhargava ◽  
G. T. Hahn ◽  
C. A. Rubin

This paper describes a two-dimensional (plane strain) elastic-plastic finite element model of rolling contact that embodies the elastic-perfectly plastic, cycle and amplitude-independent material of the Merwin and Johnson theory, but is rigorous with respect to equilibrium and continuity requirements. The rolling contact is simulated by translating a semielliptical pressure distribution. Both Hertzian and modified Hertzian pressure distributions are used to estimate the effect of plasticity on contact width and the continuity of the indentor-indentation interface. The model is tested for its ability to reproduce various features of the elastic-plastic indentation problem and the stress and strain states of single rolling contacts. This paper compares the results derived from the finite element analysis of a single, frictionless rolling contact at p0/k = 5 with those obtained from the Merwin and Johnson analysis. The finite element calculations validate basic assumptions made by Merwin and Johnson and are consistent with the development of “forward” flow. However, the comparison also reveals significant differences in the distribution of residual stress and strain components after a single contact cycle.


1988 ◽  
Vol 110 (1) ◽  
pp. 44-49 ◽  
Author(s):  
G. Ham ◽  
C. A. Rubin ◽  
G. T. Hahn ◽  
V. Bhargava

The stresses, strains, and deformations produced by repeated, two-dimensional rolling-sliding contact are analyzed using a modified finite element model developed by Bhargava et al. [1]. Rolling and sliding are simulated by translating an appropriate set of normal and tangential surface tractions across an elastic-perfectly plastic half space. The study examines a peak-pressure-to-shear strength ratio of po/k = 4.5 and normal to tangential force ratios of T/N = 0.20 and T/N = 0.17. The calculations describe the residual stresses, displacements and the continuing cyclic radial, shear and equivalent strains generated at various depths in the rim. The results are compared with previous calculations by Johnson and Jefferis [2] of rolling-sliding contact and with pure rolling. The present work predicts much higher deformations than previously calculated.


1990 ◽  
Vol 57 (1) ◽  
pp. 57-65 ◽  
Author(s):  
S. M. Kulkarni ◽  
G. T. Hahn ◽  
C. A. Rubin ◽  
V. Bhargava

This paper describes a three-dimensional elastoplastic finite element model of repeated, frictionless rolling contact. The model treats a sphere rolling on an elastic-perfectly plastic and an elastic-linear-kinematic-hardening plastic, semi-infinite half space. The calculations are for a relative peak pressure (po/k) = 4.68 (the theoretical shakedown limit for perfect plasticity). Three-dimensional rolling contact is simulated by repeatedly translating a hemispherical (Hertzian) pressure distribution across an elastoplastic semi-infinite half space. The semi-infinite half space is represented by a finite mesh with elastic boundaries. The calculations describe the distortion of the rim, the residual stress-strain distributions, stress-strain histories, and the cyclic plastic strain ranges in the vicinity of the contact.


1985 ◽  
Vol 52 (1) ◽  
pp. 75-82 ◽  
Author(s):  
V. Bhargava ◽  
G. T. Hahn ◽  
C. A. Rubin

This paper presents finite element analyses of two-dimensional (plane strain), elastic-plastic, repeated, frictionless rolling contact. The analysis employs the elastic-perfectly plastic, cycle and strain-amplitude-independent material used in the Merwin and Johnson analysis but avoids several assumptions made by these workers. Repeated rolling contacts are simulated by multiple translations of a semielliptical Hertzian pressure distribution. Results at p0/k = 3.5, 4.35, and 5.0 are compared to the Merwin and Johnson prediction. Shakedown is observed at p0/k = 3.5, but the comparisons reveal significant differences in the amount and distribution of residual shear strain and forward flow at p0/k = 4.35 and p0/k = 5.0. The peak incremental, shear strain per cycle for steady state is five times the value calculated by Merwin and Johnson, and the plastic strain cycle is highly nonsymmetric.


2021 ◽  
Vol 39 (3A) ◽  
pp. 426-439
Author(s):  
Saad A. Al-Taan ◽  
Ayad A. Abdul-Razzak

This paper presents a study on the behavior of fiber reinforced concrete slabsusing finite element analysis. A previously published finite element program is used for the nonlinear analysis by including the steel fiber concrete properties. Concrete is represented by degenerated quadratic thick shell element, which is the general shear deformable eight node serendipity element, and the thickness is divided into layers. An elastic perfectly plastic and strain hardening plasticity approach are used to model the compression behavior of concrete.The reinforcing bars were smeared within the concrete layers and assumed as either an elastic perfectly plastic material or as an elastic-plastic material with linear strain hardening. Cracks initiation is predicted using a tensile strength criterion. The tension stiffening effect of the steel fibers is simulated using a descending parabolic stress degradation function, which is based on the fracture energy concept. The effect of cracking in reducing the shear modulus and the compressive strength of concrete parallel to the crack direction is considered. The numerical results showedgood agreement with published experimental results for two fibrous reinforced concrete slabs.


2011 ◽  
Vol 243-249 ◽  
pp. 2409-2414
Author(s):  
Min Zhou ◽  
Zhong Fu Wang ◽  
Si Wei Wang

In this paper, in order to analyze the capability of squeezed branch pile under different work condition and the cooperation mechanism between the pile and soil, non-liner numerical simulation was carried out using ANSYS. In the finite element model, the elastic-perfectly plastic Drucker-Prager material was assumed for soil. Contact interface elements were placed between the pile and soil. It showed that the squeezed branches took lots of the load, and the ratio it took was related to the load and the elastic modulus of soil; the plastic section of the soil was run-through from bottom to the top; the horizontal displacement of the top soil was moved to the pile, but the horizontal displacement of the soil of the bottom was moved away from the pile; the squeezed branch will break away from the soil above the squeezed branch when the load was at a certain value.


Author(s):  
J. Rodriguez ◽  
M. Him

Abstract This paper presents a finite element mesh generation algorithm (PREPAT) designed to automatically discretize two-dimensional domains. The mesh generation algorithm is a mapping scheme which creates a uniform isoparametric FE model based on a pre-partitioned domain of the component. The proposed algorithm provides a faster and more accurate tool in the pre-processing phase of a Finite Element Analysis (FEA). A primary goal of the developed mesh generator is to create a finite element model requiring only essential input from the analyst. As a result, the generator code utilizes only a sketch, based on geometric primitives, and information relating to loading/boundary conditions. These conditions represents the constraints that are propagated throughout the model and the available finite elements are uniformly mapped in the resulting sub-domains. Relative advantages and limitations of the mesh generator are discussed. Examples are presented to illustrate the accuracy, efficiency and applicability of PREPAT.


2009 ◽  
Vol 44 (6) ◽  
pp. 407-416 ◽  
Author(s):  
P J Budden ◽  
Y Lei

Limit loads for a thick-walled cylinder with an internal or external fully circumferential surface crack under pure axial load are derived on the basis of the von Mises yield criterion. The solutions reproduce the existing thin-walled solution when the ratio between the cylinder wall thickness and the inside radius tends to zero. The solutions are compared with published finite element limit load results for an elastic–perfectly plastic material. The comparison shows that the theoretical solutions are conservative and very close to the finite element data.


1985 ◽  
Vol 107 (1) ◽  
pp. 13-18 ◽  
Author(s):  
B. V. Kiefer ◽  
P. D. Hilton

Capabilities for the analysis of combined viscous and plastic behavior have been added to an existing finite element computer program for two-dimensional elastic-plastic calculations. This program (PAPSTB) has been formulated for elastic-plastic stress and deformation analyses of two-dimensional and axisymmetric structures. It has the ability to model large strains and large deformations of elastic-perfectly plastic, multi-linear hardening, or power-hardening materials. The program is based on incremental plasticity theory with a von Mises yield criterion. Time dependent behavior has been introduced into the PAPSTB program by adding a viscous strain increment to the elastic and plastic strain increment to form the total strain increment. The viscous calculations presently employ a power-law relationship between the viscous strain rate and the effective stress. The finite element code can be easily modified to handle more complex viscous models. The Newmark method for time integration is used, i.e., an input parameter is included which enables the user to vary the time domain approximation between forward (explicit) and backward (implicit) difference. Automatic time stepping is used to provide for stability in the viscous calculations. It is controlled by an input parameter related to the ratio of the current viscous strain increment to the total strain. The viscoplastic capabilities of the PAPSTB program are verified using the axisymmetric problem of an internally pressurized, thick-walled cylinder. The transient viscoplastic case is analyzed to demonstrate that the elastic-perfectly plastic solution is obtained as a steady-state condition is approached. The influence of varying the time integration parameter for transient viscoplastic calculations is demonstrated. In addition, the effects of time step on solution accuracy are investigated by means of the automatic time stepping algorithm in the program. The approach is then applied to a simple forging problem of cylinder upsetting.


2000 ◽  
Author(s):  
Bhavani V. Sankar ◽  
Manickam Narayanan ◽  
Abhinav Sharma

Abstract Nonlinear finite element analysis was used to simulate compression tests on sandwich composites containing debonded face sheets. The core was modeled as an elastic-perfectly-plastic material, and the face-sheet as elastic isotropic. The effects of core plasticity, face-sheet and core thickness, and debond length on the maximum load the beam can carry were studied. The results indicate that the core plasticity is an important factor that determines the maximum load.


Sign in / Sign up

Export Citation Format

Share Document