A Finite Element to Perform 3D Shakedown Analysis for Limited Kinematic Hardening Materials

The use of the Design by Analysis (DBA) route is a modern trend in pressure vessel and piping international codes in mechanical engineering. However, to apply the DBA to structures under variable mechanical and thermal loads, it is necessary to assure that the plastic collapse modes, alternate plasticity and incremental collapse (with instantaneous plastic collapse as a particular case), be precluded. The tool available to achieve this target is the shakedown theory. Unfortunately, the practical numerical applications of the shakedown theory result in very large nonlinear optimization problems with nonlinear constraints. Precise, robust and efficient algorithms and finite elements to solve this problem in finite dimension has been a more recent achievements. However, to solve real problems in an industrial level, it is necessary also to consider more realistic material properties as well as to accomplish 3D analysis. Limited kinematic hardening, is a typical property of the usual steels and it should be considered in realistic applications. In this paper, a new finite element with internal thermodynamical variables to model kinematic hardening materials is developed and tested. This element is a mixed ten nodes tetrahedron and through an appropriate change of variables is possible to embed it in a shakedown analysis software developed by Zouain and co-workers for elastic ideally-plastic materials, and then use it to perform 3D shakedown analysis in cases with limited kinematic hardening materials

The Influence of Limited Kinematic Hardening in Shakedown Analysis

The use of the Design by Analysis concept is a trend in modern pressure vessel and piping calculations. DBA flexibility allow us to deal with unexpected configurations detected at in-service inspections. It is also important, in life extension calculations, when deviations of the original standard hypotesis adopted initially in Design by Formula, can happen. To apply the DBA to structures under variable mechanic and thermal loads, it is necessary that, alternate plasticity and incremental collapse (with instantaneous plastic collapse as a particular case), be precluded. These are two basic failure modes considered by ASME or European Standards in DBA. The shakedown theory is the tool available to achieve this goal. In order to apply it, is necessary only the range of the variable loads and the material properties. Precise, robust and efficient algorithms to solve the very large nonlinear optimization problems generated in numerical applications of the shakedown theory is a recent achievement. Zouain and co-workers developed one of these algorithms for elastic ideally-plastic materials. But, it is necessary to consider more realistic material properties in real practical applications. This paper shows an enhancement of this algorithm to dealing with limited kinematic hardening, a typical property of the usual steels. This is done using internal thermodynamic variables. A discrete algorithm is obtained using a plane stress, mixed finite element, with internal variable. An example, a beam encased in an end, under constant axial force and variable moment is presented to show the importance of considering the limited kinematic hardening in a shakedown analysis.

Shakedown Analysis Considering Limited Kinematic Hardening Materials

In design or safety assessment of mechanical structures, the use of the Design by Analysis (DBA) route is a modern trend. However, for making possible to apply DBA to structures under variable loads, two basic failure modes considered by ASME or European Standards must be precluded. Those modes are the alternate plasticity and incremental collapse (with instantaneous plastic collapse as a particular case). Shakedown theory is a tool that permit us to assure that those kinds of failures will be avoided. However, in practical applications, very large nonlinear optimization problems are generated. Due to this facts, only in recent years have been possible to obtain algorithms sufficiently accurate, robust and efficient, for dealing with this class of problems. In this paper, one of these shakedown algorithms, developed for dealing with elastic ideally-plastic structures, is enhanced to include limited kinematic hardening, a more realistic material behavior. This is done in the continuous model by using internal thermodynamic variables. A corresponding discrete model is obtained using an axisymmetric mixed finite element with an internal variable. A thick wall sphere, under variable thermal and pressure loads, is used in an example to show the importance of considering the limited kinematic hardening in the shakedown calculations

Constitutive modeling of kinematic-hardening and damage in solder joints

Solder constitutive models are important as they are widely used in FEA simulations to predict the lifetime of soldered assemblies. This paper briefly reviews some common constitutive laws to capture creep in solder and presents work on laws capturing both kinematic hardening and damage. Inverse analysis is used to determine constants for the kinematic hardening law which match experimental creep curves. The mesh dependence of the damage law is overcome by using volume averaging and is applied to predict the crack path in a thermal cycled resistor component

On the definition of an kinematic hardening effect graph for sheet metal forming process simulations

The objective of this work is to develop a kinematic hardening effect graph (KHEG) which can be used to evaluate the effect of kinematic hardening on the model accuracy of numerical sheet metal forming simulations and this without the need of complex material characterisation. The virtual manufacturing process design and optimisation depends on the accuracy of the constitutive models used to represent material behaviour. Under reverse strain paths the Bauschinger effect phenomenon is modelled using kinematic hardening models. However, due to the complexity of the experimental testing required to characterise this phenomenon in this work the KHEG is presented as an indicator to evaluate the potential benefit of carrying out these tests. The tool is validated with the classic three point bending process and the U-channel width drawbead process. In the same way, the capability of the KHEG to identify effects in forming processes that do not include forming strain reversals is identified.

Anisotropic hardening model based on non-associated flow rule and combined nonlinear kinematic hardening for sheet materials

 A material model for more effective analysis of plastic deformation of sheet materials is presented in this paper. The model is capable of considering the following aspects of plastic deformation behavior of sheet materials: the anisotropy in yielding stresses in different directions by using a quadratic yield function (based on Hill’s 1948 model and stress ratios), the anisotropy in work hardening by introducing non-constant flow stress hardening in different directions, the anisotropy in plastic strains in different directions by using a quadratic plastic potential function and non-associated flow rule (based on Hill’s 1948 model and plastic strain ratios, r-values), and finally some of the cyclic hardening phenomena such as Bauschinger’s effect and transient behavior for reverse loading by using a coupled nonlinear kinematic hardening (so-called Armstrong-Frederick-Chaboche model). Basic fundamentals of the plasticity of the model are presented in a general framework. Then, the model adjustment procedure is derived for the plasticity formulations. Also, a generic numerical stress integration procedure is developed based on backward-Euler method (so-called multistage return mapping algorithm). Different aspects of the model are verified for DP600 steel sheet. Results show that the new model is able to predict the sheet material behavior in both anisotropic hardening and cyclic hardening regimes more accurately. By featuring the above-mentioned facts in the presented constitutive model, it is expected that more accurate results can be obtained by implementing this model in computational simulations of sheet material forming processes. For instance, more precise results of springback prediction of the parts formed from highly anisotropic hardened materials or that of determining the forming limit diagrams is highly expected by using the developed material model.

A non-associated plasticity model with anisotropic and nonlinear kinematic hardening for simulation of sheet metal forming

A material model for more thorough analysis of plastic deformation of sheet materials is presented in this paper. This model considers the following aspects of plastic deformation behavior of sheet materials: (1) the anisotropy in yield stresses and in work hardening by using Hill's 1948 quadratic yield function and non-constant stress ratios which leads to different flow stress hardening in different directions, (2) the anisotropy in plastic strains by using a quadratic plastic potential function and non-associated flow rule, also based on Hill's 1948 model and r-values, and (3) the cyclic hardening phenomena such as the Bauschinger effect, permanent softening and transient behavior for reverse loading by using a coupled nonlinear kinematic hardening model. Plasticity fundamentals of the model were derived in a general framework and the model calibration procedure was presented for the plasticity formulations. Also, a generic numerical stress integration procedure was developed based on backward-Euler method, so-called multi-stage return mapping algorithm. The model was implemented in the framework of the finite element method to evaluate the simulation results of sheet metal forming processes. Different aspects of the model were verified for two sheet metals, namely DP600 steel and AA6022 aluminum alloy. Results show that the new model is able to accurately predict the sheet material behavior for both anisotropic hardening and cyclic hardening conditions. The drawing of channel sections and the subsequent springback were also simulated with this model for different drawbead configurations. Simulation results show that the current non-associated anisotropic hardening model is able to accurately predict the sidewall curl in the drawn channel sections.

A large strain anisotropic elastoplastic continuum theory for nonlinear kinematic hardening and texture evolution

In this paper we present a continuum theory for large strain anisotropic elastoplasticity based on a decomposition of the modified plastic velocity gradient into energetic and dissipative parts. The theory includes the Armstrong and Frederick hardening rule as well as multilayer models as special cases even for large strain anisotropic elastoplasticity. Texture evolution may also be modelled by the formulation, which allows for a meaningful interpretation of the terms of the dissipation equation

Ratchetting of railhead in the vicinity of the gap of the insulated rail joints

Insulated Rail Joints (IRJs) are designed to electrically isolate two rails in rail tracks to control the signalling system for safer train operations. Unfortunately the gapped section of the IRJs is structurally weak and often fails prematurely especially in heavy haul tracks, which adversely affects service reliability and efficiency. The IRJs suffer from a number of failure modes; the railhead ratchetting at the gap is, however, regarded as the root cause and attended to in this thesis. Ratchetting increases with the increase in wheel loads; in the absence of a life prediction model, effective management of the IRJs for increased wagon wheel loads has become very challenging. Therefore, the main aim of this thesis is to determine method to predict IRJs' service life. The distinct discontinuity of the railhead at the gap makes the Hertzian theory and the rolling contact shakedown map, commonly used in the continuously welded rails, not applicable to examine the metal ratchetting of the IRJs. Finite Element (FE) technique is, therefore, used to explore the railhead metal ratchetting characteristics in this thesis, the boundary conditions of which has been determined from a full scale study of the IRJ specimens under rolling contact of the loaded wheels. A special purpose test set up containing full-scale wagon wheel was used to apply rolling wheel loads on the railhead edges of the test specimens. The state of the rail end face strains was determined using a non-contact digital imaging technique and used for calibrating the FE model. The basic material parameters for this FE model were obtained through independent uniaxial, monotonic tensile tests on specimens cut from the head hardened virgin rails. The monotonic tensile test data have been used to establish a cyclic load simulation model of the railhead steel specimen; the simulated cyclic load test has provided the necessary data for the three decomposed kinematic hardening plastic strain accumulation model of Chaboche. A performance based service life prediction algorithm for the IRJs was established using the plastic strain accumulation obtained from the Chaboche model. The predicted service lives of IRJs using this algorithm have agreed well with the published data. The finite element model has been used to carry out a sensitivity study on the effects of wheel diameter to the railhead metal plasticity. This study revealed that the depth of the plastic zone at the railhead edges is independent of the wheel diameter; however, large wheel diameter is shown to increase the IRJs' service life.

Failure of discontinuous railhead edges due to plastic strain accumulation

Railhead is perhaps the highest stressed civil infrastructure due to the passage of heavily loaded wheels through a very small contact patch. The stresses at the contact patch cause yielding of the railhead material and wear. Many theories exist for the prediction of these mechanisms of continuous rails; this process in the discontinuous rails is relatively sparingly researched. Discontinuous railhead edges fail due to accumulating excessive plastic strains. Significant safety concern is widely reported as these edges form part of Insulated Rail Joints (IRJs) in the signalling track circuitry. Since Hertzian contact is not valid at a discontinuous edge, 3D finite element (3DFE) models of wheel contact at a railhead edge have been used in this research. Elastic–plastic material properties of the head hardened rail steel have been experimentally determined through uniaxial monotonic tension tests and incorporated into a FE model of a cylindrical specimen subject to cyclic tension load- ing. The parameters required for the Chaboche kinematic hardening model have been determined from the stabilised hysteresis loops of the cyclic load simulation and imple- mented into the 3DFE model. The 3DFE predictions of the plastic strain accumulation in the vicinity of the wheel contact at discontinuous railhead edges are shown to be affected by the contact due to passage of wheels rather than the magnitude of the loads the wheels carry. Therefore to eliminate this failure mechanism, modification to the contact patch is essential; reduction in wheel load cannot solve this problem.

First Measurement of the b-Jet Cross Section in Events with a W Boson in pp̅ Collisions at √s=1.96 TeV

The cross section for jets from b quarks produced with a W boson has been measured in ppbar collision data from 1.9/fb of integrated luminosity recorded by the CDF II detector at the Tevatron. The W+b-jets process poses a significant background in measurements of top quark production and prominent searches for the Higgs boson. We measure a b-jet cross section of 2.74 +- 0.27(stat.) +- 0.42(syst.) pb in association with a single flavor of leptonic W boson decay over a limited kinematic phase space. This measured result cannot be accommodated in several available theoretical predictions.

Constitutive Response of Passivated Copper Films: Experiments and Analyses

The constitutive behavior of passivated copper films is studied. Stresses in copper films of thickness ranging from 1000 nm to 40 nm, passivated with silicon oxide on a quartz or silicon substrate, were measured using the curvature method. The thermal cycling spans a temperature range from - 196 to 600°C. It is seen that the strong relaxation at high temperatures normally found in unpassivated films is nonexistent for passivated films. The copper film did not show any rate-dependent effect over a range of heating/cooling rate from 5 to 25°C/min. Further analyses showed that significant strain hardening exists during the course of thermal loading. In particular, the measured stress- temperature response can only be fitted with a kinematic hardening model, if a simple constitutive law within the continuum plasticity framework is to be used. The analytic procedures for extracting the film properties are presented. Implications to stress modeling of copper interconnects in actual devices are discussed.

Experimental and theoretical comparison between static and dynamic torsional soil tests

The dynamic properties of dry Leighton Buzzard sand have been investigated using a resonant column test apparatus. These data are compared with very low frequency cyclic tests on identical specimens of sand. The comparison indicates that the properties of dry sand are independent of frequency. A simple one-dimensional model of kinematic hardening plasticity is used to predict the dynamic behaviour of the sand. The input parameters for this model are based on the results of static tests. These may be conducted on standard laboratory equipment with only minor modifications. The predictions are in good agreement with the measured data.