Elastic–plastic stresses in linearly hardening rotating solid disks of variable thickness

The distribution of stress, displacement and plastic strain in a rotating elastic–plastic solid disk of variable thickness in a power function form is investigated. The analysis is based on Tresca's yield condition, its associated flow rule and linear strain hardening material behavior. An analytical solution is obtained and numerical results are presented for different values of the geometric parameters. The validity of the solution is demonstrated by comparing the results with those for a uniform thickness disk available in the literature.     

Elastoplastic analysis of nonlinearly hardening variable thickness annular disks under external pressure

Based on von Mises’ yield criterion, deformation theory of plasticity and Swift’s hardening law, elasto-plastic deformation of variable thickness annular disks subjected to external pressure is studied. A nonlinear shooting method using Newton’s iterations with numerically approximated tangent is designed for the solution of the problem. Considering a thickness profile in the form of a general parabolic function, the condition of occurrence of plastic deformation at the inner and outer edges of the annular disk is investigated. A critical disk profile is determined and the corresponding elastic–plastic stresses as well as the residual stress distribution upon removal of the applied pressure are computed and discussed.     

Analysis of functionally graded rotating disks with variable thickness

Elastic solutions for axisymmetric rotating disks made of functionally graded material with variable thickness are presented. The material properties and disk thickness profile are assumed to be represented by two power-law distributions. In the case of hollow disk, based on the form of the power-law distribution for the mechanical properties of the constituent components and the thickness profile function, both analytical and semi-analytical solutions are given under free–free and fixed-free boundary conditions. For the solid disk, only semi-analytical solution is presented. The effects of the material grading index and the geometry of the disk on the stresses and displacements are investigated. It is found that a functionally graded rotating disk with parabolic or hyperbolic convergent thickness profile has smaller stresses and displacements compared with that of uniform thickness. It is seen that the maximum radial stress for the solid functionally graded disk with parabolic thickness profile is not at the centre like uniform thickness disk. Results of this paper suggest that a rotating functionally graded disk with parabolic concave or hyperbolic convergent thickness profile can be more efficient than the one with uniform thickness.     

Numerical analysis of elastic–plastic rotating disks with arbitrary variable thickness and density

A unified numerical method is developed in this article for the analysis of deformations and stresses in elastic–plastic rotating disks with arbitrary cross-sections of continuously variable thickness and arbitrarily variable density made of nonlinear strain-hardening materials. The method is based on a polynomial stress–plastic strain relation, deformation theory in plasticity and Von Mises’ yield condition. The governing equation is derived from the basic equations of the rotating disks and solved using the Runge–Kutta algorithm. The proposed method is applied to calculate the deformations and stresses in various rotating disks. These disks include solid disks with constant thickness and constant density, annular disks with constant thickness and constant density, nonlinearly variable thickness and nonlinearly variable density, linearly tapered thickness and linearly variable density, and a combined section of continuously variable thickness and constant density. The computed results are compared to those obtained from the finite element method and the existing approaches. A very good agreement is found between this research and the finite element analysis. Due to the simplicity, effectiveness and efficiency of the proposed method, it is especially suitable for the analysis of various rotating disks.     

Limit angular velocities of variable thickness rotating disks

Elastic and plastic limit angular velocities are calculated for rotating disks of variable thickness in power function form. Analytical solution is obtained and used to calculate elastic limit angular velocities of variable thickness rotating annular disks and annular disks with rigid inclusion. An efficient numerical solution procedure is designed and used to obtain the elastic limit angular velocities of variable thickness rotating solid disks. Von Mises yield criterion and its flow rule is used to estimate plastic limit angular velocities. Both linear and nonlinear hardening material behaviors are treated numerically. The results are verified by comparing with those of uniform thickness rotating solid disks available in the literature. Elastic and plastic limit angular velocities are found to increase with the reduction of the disk thickness at the edge as well as the reduction in the disk mass due to the shape of the profile.     

A centrally cracked thin circular disk, Part I: 3D Elastic–plastic finite element analysis

A full field solution, based on small deformation, three-dimensional elastic–plastic finite element analysis of the centrally cracked thin disk under mode I loading has been performed. The solution for the stresses under small-scale yielding and lo!cally fully plastic state has been compared with the HRR plane stress solution. At the outside of the 3D zone, within a distance of rσ_o/J=18, HRR dominance is maintained in the presence of a significant amount of compressive stress along the crack flanks. Ahead of this region, the HRR field overestimate the stresses. These results demonstrate a completely reversed state of stress in the near crack front compared to that in the plane strain case. The combined effect of geometry and finite thickness of the specimen on elastic–plastic crack tip stress field has been explored. To the best of our knowledge, such an attempt in the published literature has not been made yet. For the qualitative assessment of the results some of the field parameters have been compared to the available experimental results of K, gives a fair estimate of the crack opening stress near the crack front at a distance of order 10⁻² in. On the basis of this analysis, the Linear Elastic Fracture Mechanics approach has been adopted in analyzing the fatigue crack extension experiments performed in the disk (Part II).     

Von mises yield criterion and nonlinearly hardening variable thickness rotating annular disks with rigid inclusion

A computational model is developed to investigate inelastic deformations of variable thickness rotating annular disks mounted on rigid shafts. The von Mises yield condition and its flow rule are combined with Swift’s hardening law to simulate nonlinear hardening material behavior. An efficient numerical solution procedure is designed and used throughout to handle the nonlinearities associated with the von Mises yield condition and the boundary condition at the shaft–annular disk interface. The results of the computations are verified by comparison with an analytical solution employing Tresca’s criterion available in the literature. Inelastic stresses and deformations are calculated for rotating variable thickness disks described by two different commonly used disk profile functions i.e. power and exponential forms. Plastic limit angular velocities for these disks are calculated for different values of the geometric and hardening parameters. These critical angular velocities are found to increase as the edge thickness of the disk reduces. Lower plastic limit angular velocities are obtained for disks made of nonlinearly hardening materials.     

A model for the through-thickness elastic–plastic behaviour of paper

An elastic–plastic material model for the out-of-plane mechanical behaviour of paper is presented. This model enables simulation the elastic–plastic behaviour under high compressive loads in the through-thickness direction (ZD). Paper does not exhibit a sharp transition from elastic to elastic–plastic behaviour. This makes it advantageous to define critical stress states based on failure stresses rather than yield stresses. Moreover, the failure stress in out-of-plane shear is strongly affected by previous plastic through-thickness compression. To cover these two features, a model based on the idea of a bounding surface that grows in size with plastic compression is proposed. Here, both the bounding and the yield surfaces are suggested as parabolas in stress space. While the bounding surface is open for compressive loads, the yield surface is bordered by the maximum applied through-thickness compression.     

Residual stress field and reduction of stress intensity factors in cold-worked holes

Closed-form and semi-analytical solutions are obtained for the residual stress distributions in a plate caused by pressure acting on a central circular hole, representing the cold-work process. The material is elastic–perfectly plastic. Both Tresca and von Mises yield criteria are used and the corresponding residual stress distributions are compared. The relation between the dimension of the plastic zone and the value of internal pressure is presented. The relation between the magnitude of the residual stresses and the remote uniform tensile stress required to open symmetrical radial cracks is also presented. The reduction of the stress intensity factors of cracked open and riveted holes as a function of the internal pressure applied (or mandrel radial displacement) is investigated using numerical models for both an elastic–perfectly plastic material and for an Al 2024-T3 Alclad aluminum alloy.     

Unit cell debonding analyses for arbitrary orientation of plastic anisotropy

Debonding of rigid inclusions embedded in the elastic–plastic aluminum alloy Al 2090-T3 is analyzed numerically using a unit cell model taking full account of finite strains. The cell is subjected to overall biaxial plane strain tension and periodical boundary conditions are applied to represent arbitrary orientations of plastic anisotropy. Plastic anisotropy is considered using two phenomenological anisotropic yield criteria, namely Hill [Proceedings of the Royal Society of London A 193 (1948) 281] and Barlat et al. [International Journal of Plasticity 7 (1991) 693]. For this material plastic anisotropy delays debonding compared to plastic isotropy except for the case of Hill’s yield function when the tensile directions coincided with the principal axes of anisotropy. For some inclinations of the principal axes of anisotropy relative to the tensile directions, the stress strain responses are identical but the deformation modes are mirror images of each other.     

Heterogeneity of Plastic Flow of Zirconium Alloys with a Parabolic Law of Strain Hardening

The specific features of plastic–strain macrolocalization at the stage of the parabolic law of strain hardening in samples from industrial zirconium–based alloys are considered. It is shown that in predeformed blanks, zones with a different character of plastic–strain localization are formed. It is also shown that the strain–localization macropattern can be used as a characteristic of the susceptibility of a material to further plastic form–changing, for example, upon tube rolling. The sign of fracture of alloys upon plastic deformation is revealed. The scale effect in the formation of localizedplastic–flow zones is shown and studied.     

Elastic–plastic and elastic anti-plane shear of two parallel cracks

Two infinite interacting parallel cracks in an elastic–plastic and in an elastic body under anti-plane strain (mode III) loading conditions are considered. The body is subjected to vanishing remote loading and the cracks are traction free. Closed-form solution is found for the elastic–plastic problem in terms of elementary functions, where the shape of the plastic boundary is obtained. The complete stress distribution is obtained in an inverse form i.e. physical coordinates are functions of stresses.     

Elastic–plastic FE analysis of a notched shaft under multiaxial nonproportional synchronous cyclic loading

An elastic–plastic finite element analysis is presented for a notched shaft subjected to multiaxial nonproportional synchronous cyclic tension/torsion loading. The elastic–plastic material property is described by the von Mises yield criterion and the kinematic hardening rule of Prager/Ziegler. The finite element program system ABAQUS is used to solve the boundary value problem. Special emphasis is given to explore the effects of the stress amplitude, the mean-stress, and the mutual interactions on the local stress–strain responses at the notch root.     

Displacement Field and Strain Distribution in a Rotating Annular Disk

Abstract

The displacement field and strain distribution in a thin rotating disk with constant thickness and density are found based on Mises’ yield criterion and its associated flow rule. The material of the disk is elastic-perfectly plastic and the assumption of plane stress is adopted. The solution is illustrated by an example.     

A finite thickness band method for ductile fracture analysis

We present a finite element method with a finite thickness embedded weak discontinuity to analyze ductile fracture problems. The formulation is restricted to small geometry changes. The material response is characterized by a constitutive relation for a progressively cavitating elastic–plastic solid. As voids nucleate, grow and coalesce, the stiffness of the material degrades. An embedded weak discontinuity is introduced when the condition for loss of ellipticity is met. The resulting localized deformation band is given a specified thickness which introduces a length scale thus providing a regularization of the post-localization response. Also since the constitutive relation for a progressively cavitation solid is used inside the band in the post-localization regime, the traction-opening relation across the band depends on the stress triaxiality. The methodology is illustrated through several example problems including mode I crack growth and localization and failure in notched bars. Various finite element meshes and values of the thickness of the localization band are used in the calculations to illustrate the convergence with mesh refinement and the dependence on the value chosen for the localization band thickness.     

Interaction of microcracks with a macrocrack yielded in a narrow strip

Plastic zone growth of collinear cracks has had a longstanding interest in ductile fracture. This work further considers yield zone growth in an isotropic, homogeneous elastic–perfectly plastic infinite plate containing a macrocrack with several neighboring microcracks. Normal loading is considered at distances far away from the cracks. The strip yield is adopted where the plastic zone is assumed to be confined to two narrow strips extending from the ends of a finite length crack while the microcracks are assumed to be elastic. The plastic zone length and crack opening displacement are found from asymptotic solution and compared with finite element solution.     

Near crack line elastic–plastic analysis for a infinite plate loaded by two pairs of point tensile forces

The near crack line analysis method is used to investigate a center crack loaded by two pairs of point tensile forces in an infinite plate in an elastic–perfectly plastic solid, and the analytical solutions are obtained in this paper. These solutions include: the elastic–plastic stress field near the crack line, the law that the length of the plastic zone along the crack line is varied with an external loads and the bearing capacity of an infinite plate with a center crack. The results of this paper are sufficiently precise near the crack line because the assumptions of the small scale yielding theory have not been used and no other assumptions have been taken.     

Plasticity aspects of interfacial fracture for polymer/glass specimens

Experimental results suggest that the interfacial fracture resistance is minimal for approximate near tip Mode I accompanied by positive and negative near tip Mode II. Finite-strain FE analysis is made for an elastic–plastic medium bonded to an ideally elastic medium with an interface crack. Small-scale plasticity conditions are invoked and examined in relation to the elastic–plastic stress distribution along the bond line. Plasticity engenders a tendency to turn near tip biaxiality towards pure Mode I regardless of the mixed-mode loading. High levels of hydrostatic stress are attained. For different mode mixities of the applied load, the dependence of the elastic–plastic normal bond stress on load level is examined. It is found that under positive Mode II loading, the normal bond stress σ_yy tends to saturate as the load level rises. This does not occur for Mode I and negative Mode II loading. In addition, deformation patterns inside the plastic zone are examined for mixed-mode situations. A displacement criterion based on the normal bond crack opening suggests a dependence of the critical load level on the extent of mixed mode. Under positive mode II fracture, traces of the ductile material are found at the top of the elastic substrate. Some of these conclusions appear to be consistent with the fracture patterns observed for LD-polyethylene/glass interfacial mixed-mode fracture.     

J-integral and crack driving force in elastic–plastic materials

This paper discusses the crack driving force in elastic–plastic materials, with particular emphasis on incremental plasticity. Using the configurational forces approach we identify a “plasticity influence term” that describes crack tip shielding or anti-shielding due to plastic deformation in the body. Standard constitutive models for finite strain as well as small strain incremental plasticity are used to obtain explicit expressions for the plasticity influence term in a two-dimensional setting. The total dissipation in the body is related to the near-tip and far-field J-integrals and the plasticity influence term. In the special case of deformation plasticity the plasticity influence term vanishes identically whereas for rigid plasticity and elastic-ideal plasticity the crack driving force vanishes. For steady state crack growth in incremental elastic–plastic materials, the plasticity influence term is equal to the negative of the plastic work per unit crack extension and the total dissipation in the body due to crack propagation and plastic deformation is determined by the far-field J-integral. For non-steady state crack growth, the plasticity influence term can be evaluated by post-processing after a conventional finite element stress analysis. Theory and computations are applied to a stationary crack in a C(T)-specimen to examine the effects of contained, uncontained and general yielding. A novel method is proposed for evaluating J-integrals under incremental plasticity conditions through the configurational body force. The incremental plasticity near-tip and far-field J-integrals are compared to conventional deformational plasticity and experimental J-integrals.