Interface crack growth for anisotropic plasticity with non-normality effects

A plasticity model with a non-normality plastic flow rule is used to analyze crack growth along an interface between a solid with plastic anisotropy and an elastic substrate. The fracture process is represented in terms of a traction-separation law specified on the crack plane. A phenomenological elastic–viscoplastic material model is applied, using an anisotropic yield criterion, and in each case analyzed the effect of non-normality is compared with results for the standard normality flow rule. Due to the mismatch of elastic properties across the interface the corresponding elastic solution has an oscillating stress singularity, and with conditions of small scale yielding this solution is applied as boundary conditions on the outer edge of the region analyzed. Crack growth resistance curves are calculated numerically, and the effect of the near-tip mode mixity on the steady-state fracture toughness is determined. It is found that the steady-state fracture toughness is quite sensitive to differences in the initial orientation of the principal axes of the anisotropy relative to the interface.     

On non-associative anisotropic finite plasticity fully coupled with isotropic ductile damage for metal forming

In this work, non-associative finite strain anisotropic elastoplasticity fully coupled with ductile damage is considered using a thermodynamically consistent framework. First, the kinematics of large strain based on multiplicative decomposition of the total transformation gradient using the rotating frame formulation, is recalled and different objective derivatives defined. By using different anisotropic equivalent stresses (quadratic and non-quadratic) in yield function and in plastic potential, the evolution equations for all the dissipative phenomena are deduced from the generalized normality rule applied to the plastic potential while the consistency condition is still applied to the yield function. The effect of the objective derivatives and the equivalent stresses (quadratic or non-quadratic) on the plastic flow anisotropy and the hardening evolution with damage is considered. Numerical aspects mainly related to the time integration of the fully coupled constitutive equations are discussed. Applications are made to the AISI 304 sheet metal by considering different loading paths as tensile, shear, plane tensile and bulge tests. For each loading path the effect of the rotating frame, the equivalent stress (quadratic or non-quadratic) and the normality rule (with respect to yield function or to the plastic potential) are discussed on the light of some available experimental results.     

A plane stress anisotropic plastic flow theory for orthotropic sheet metals

A phenomenological theory is presented for describing the anisotropic plastic flow of orthotropic polycrystalline aluminum sheet metals under plane stress. The theory uses a stress exponent, a rate-dependent effective flow strength function, and five anisotropic material functions to specify a flow potential, an associated flow rule of plastic strain rates, a flow rule of plastic spin, and an evolution law of isotropic hardening of a sheet metal. Each of the five anisotropic material functions may be represented by a truncated Fourier series based on the orthotropic symmetry of the sheet metal and their Fourier coefficients can be determined using experimental data obtained from uniaxial tension and equal biaxial tension tests. Depending on the number of uniaxial tension tests conducted, three models with various degrees of planar anisotropy are constructed based on the proposed plasticity theory for power-law strain hardening sheet metals. These models are applied successfully to describe the anisotropic plastic flow behavior of 10 commercial aluminum alloy sheet metals reported in the literature.     

Plastic spin associated with a non-normality theory of plasticity

A plasticity model using a vertex-type plastic flow rule on a smooth yield surface for an anisotropic solid has been proposed recently. This model is here completed by incorporating the effect of plastic spin. Simple shear with a large shear strain is one of the hardest tests on finite strain anisotropic plasticity models, and it is here shown which plastic spin expression is needed to avoid unrealistic oscillatory behavior of the shear stress under large shear strains. The idea of using non-normality with a smooth yield surface originates from a recent proposal of using an abrupt strain path change to determine the subsequent yield surface shape. For this method both polycrystal plasticity calculations and experiments have shown a vertex-type response on the apparently smooth yield surface.     

Nucleation of micro-cracks for polycrystalline metal with anisotropy: micro-stress evaluation

Presented is the local stresses on the crystallographic plane as they are influenced by the metal fracture with anisotropy. It is based on the nucleation of micro-cracks and its unstable equilibrium in a polycrystal with texture. Crystallographic texture causes non-uniform distribution of the crack nucleus orientations owing to their preference to expand or open on certain crystallographic planes. This is the main cause of anisotropy of cleavage fracture stress of textured metal. “Oriented” micro-stresses in textured metal contribute to the anisotropic effect. In view of what has been said, accumulated plastic strains at fracture is analysed.     

弹塑性有限变形的拟流动理论

本文提出一种弹塑性有限变形的拟流动理论。该理论从正交性法则出发，通过引入“拟弹性模量”和模量衰减函数并改进应变率的弹塑性分解，实现了由有限变形Ｐｒａｎｄｔｌ－Ｒｅｕｓｓ流动理论（Ｊ２Ｆ）向基于非正交法则的率形式形变理论（Ｊ２Ｄ）的合理的光滑过渡；并适用于初始及后继各向异性变形分析。在特殊条件下，可退化为Ｊ２Ｆ、Ｊ２Ｄ理论以及由任意各向异性屈服函数描述的流动理论。将该理论用于韧性金属平面应力／应变拉伸失稳与变形局部化的有限元模拟，并与理论分析及实验结果相比较，表明了本文理论的正确性。     

Influence of void nucleation on ductile shear fracture at a free surface

An approximate continuum model of a ductile, porous material is used to study the influence of the nucleation and growth of micro-voids on the formation of shear bands and the occurrence of surface shear fracture in a solid subject to plane strain tension. Bifurcation into diffuse modes is analysed for a plane strain tensile specimen described by these constitutive relations, which account for a considerable plastic dilatancy due to void growth and for the possibility of non-normality of the plastic flow law. In particular, bifurcation into surface wave modes and the possible influence of such modes triggering shear bands is investigated. For solids with initial imperfactions such as a surface undulation, a local material inhomogeneity on an inclusion colony, the inception and growth of plastic flow localization is analysed numerically. Both the formation of void-sheets and the final growth of cracks in the shear bands is described numerically. Some special features of shear band development in the solid obeying non-normality are studied by a simple model problem.     

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.     

A new yield function and a hydrostatic stress-controlled void nucleation model for porous solids with pressure-sensitive matrices

A macroscopic yield function for porous solids with pressure-sensitive matrices modeled by Coulomb's yield function was obtained by generalizing Gurson's yield function with consideration of the hydrostatic yield stress of a spherical thick-walled shell and by fitting the finite element results of the yield stresses of a voided cube. The macroscopic yield function is valid for the negative hydrostatic stress as well as for the positive hydrostatic stress. From the yield function, a plastic potential function for the porous solids was derived either for plastic normality flow or for plastic non-normality flow of the pressure-sensitive matrices. In addition, void nucleation was modeled by a normal distribution function with the macroscopic hydrostatic stress regarded as a controlling stress. This set of constitutive relations was implemented into a finite element code abaqus as a user material subroutine to analyze the cavitation and the deformation behavior of a rubber-modified epoxy around a crack tip under the Mode I plane strain conditions. By comparing the cavitation zone and the plastic zone obtained in the analysis with those observed in an experiment, the mean stress and the standard deviation for the void nucleation model could be determined. The cavitation and the deformation behavior of the rubber-modified epoxy were also analyzed around notches under four-point bending. The size and shape of the cavitation zone and the plastic zone were shown to be in good agreement with those observed in an experiment.     

The ductile fracture of anisotropic materials

A theoretical solution is obtained which adapts the Dugdale model to anisotropic materials. The effect of anisotropy is shown to modify each isotropic field equation by amultiplicative constant which is a function of the anisotropic-material constants. A limit on the validity of a Dugdale-type solution is found using a von Mises' yield criterion which implies that the crack and plastic zone will extend along a direction other than the crack line. As long as the material is within pirescribed limits, the Dugdale finiteness condition is shown to be affected by material anisotropy only in the calculation of a yield stress. Experimental results on crack-opening displacements are presented which indicate that a large degree of anisotropy is necessary to cause significant deviations from isotropic theory. Results are presented which indicate that anisotropy does affect the direction of crack propagation.     

Forming of aluminum alloys at elevated temperatures – Part 1: Material characterization

A temperature-dependent anisotropic material model for use in a coupled thermo-mechanical finite element analysis of the forming of aluminum sheets was developed. The anisotropic properties of the aluminum alloy sheet AA3003-H111 were characterized for a range of temperatures 25–260 °C (77–500 °F) and for different strain rates. Material hardening parameters (flow rule) and plastic anisotropy parameters (R₀, R₄₅ and R₉₀) were calculated using standard ASTM uniaxial tensile tests. From this experimental data, the anisotropy coefficients for the Barlat YLD96 yield function [Barlat, F., Maeda, Y., Chung, K., Yanagawa, M., Brem, J.C., Hayashida, Y., Lege, D.J., Matsui, K., Murtha, S.J., Hattori, S., Becker, R.C., Makosey, S., 1997a. Yield function development for aluminum alloy sheets. J. Mech. Phys. Solids 45 (11/12), 1727–1763] in the plane stress condition were calculated for several elevated temperatures. Curve fitting was used to calculate the anisotropy coefficients of Barlat’s YLD96 model and the hardening parameters as a function of temperature. An analytical study of the accuracy and usability of this curve fitting technique is presented through the calculation of plastic anisotropy R-parameters and yield function plots at different temperatures.     

Modeling of deformation induced anisotropy in free-end torsion

The main purpose of this work is to develop a phenomenological model, which accounts for the evolution of the elastic and plastic properties of fcc polycrystals due to a crystallographic texture development and predicts the axial effects in torsion experiments. The anisotropic portion of the effective elasticity tensor is modeled by a growth law. The flow rule depends on the anisotropic part of the elasticity tensor. The normalized anisotropic part of the effective elasticity tensor is equal to the 4th-order coefficient of a tensorial Fourier expansion of the crystal orientation distribution function. Hence, the evolution of elastic and viscoplastic properties is modeled by an evolution equation for the 4th-order moment tensor of the orientation distribution function of an aggregate of cubic crystals. It is shown that the model is able to predict the plastic anisotropy that leads to the monotonic and cyclic Swift effect. The predictions are compared to those of the Taylor–Lin polycrystal model and to experimental data. In contrast to other phenomenological models proposed in the literature, the present model predicts the axial effects even if the initial state of the material is isotropic.     

Plastic potentials for anisotropic porous solids

The aim of this paper is to incorporate plastic anisotropy into constitutive equations of porous ductile metals. It is shown that plastic anisotropy of the matrix surrounding the voids in a ductile material could have an influence on both effective stress–strain relation and damage evolution. Two theoretical frameworks are envisageable to study the influence of plastic flow anisotropy: continuum thermodynamics and micromechanics. By going through the Rousselier thermodynamical formulation, one can account for the overall plastic anisotropy, in a very simple manner. However, since this model is based on a weak coupling between plasticity and damage dissipative processes, it does not predict any influence of plastic anisotropy on cavity growth, unless a more suitable choice of the thermodynamical force associated with the damage parameter is made. Micromechanically-based models are then proposed. They consist of extending the famous Gurson model for spherical and cylindrical voids to the case of an orthotropic material. We derive an upper bound of the yield surface of a hollow sphere, or a hollow cylinder, made of a perfectly plastic matrix obeying the Hill criterion. The main findings are related to the so-called ‘scalar effect’ and ‘directional effect’. First, the effect of plastic flow anisotropy on the spherical term of the plastic potential is quantified. This allows a classification of sheet materials with regard to the anisotropy factor h; this is the scalar effect. A second feature of the model is the plasticity-induced damage anisotropy. This results in directionality of fracture properties (‘directional effect’). The latter is mainly due to the principal Hill coefficients whilst the scalar effect is enhanced by ‘shear’ Hill coefficients. Results are compared to some micromechanical calculations using the finite element method.     

A simple plasticity model for prediction of non-coaxial flow of sand

A bounding surface plasticity model for non-coaxiality, another aspect of anisotropic behavior of sands under rotation of principal stress axes; is developed in the critical state framework. Numerous experimental evidences exist that corroborate dependence of plastic shear strain rate direction on inherent fabric anisotropy. At first, general expressions for plastic strain rate with respect to possible emerge of non-coaxial flow are obtained. Consequently, using an anisotropy state parameter that is specially developed for this model and accounts for the interaction between imposed loading and soil fabric; effect of anisotropy on plastic flow direction is taken into account. Besides, novel circumstances are proposed for plastic modulus and dilatancy under rotation of principal stress axes. Finally, it is shown that the model is able to simulate successfully the non-coaxial behavior of sands subjected to principal stress axes rotation.     

Micromechanical analysis of friction anisotropy in rough elastic contacts

Computational contact homogenization approach is applied to study friction anisotropy resulting from asperity interaction in elastic contacts. Contact of rough surfaces with anisotropic roughness is considered with asperity contact at the micro scale being governed by the isotropic Coulomb friction model. Application of a micro-to-macro scale transition scheme yields a macroscopic friction model with orientation- and pressure-dependent macroscopic friction coefficient. The macroscopic slip rule is found to exhibit a weak non-associativity in the tangential plane, although the slip rule at the microscale is associated in the tangential plane. Counterintuitive effects are observed for compressible materials, in particular, for auxetic materials.     

Thermoelastic stress intensity factors for a cracked layer between two different anisotropic solids

Steady-state anisotropic thermoelasticity equations are used to obtain the stress intensity factors for a cracked layer sandwiched between two different anisotropic elastic solids. The anisotropy is assumed to arise from discrete fibers whose orientation could alter with reference to the crack edges. A generalized plane deformation prevails in the dissimilar media domain with a line of discontinuity disturbing a uniform heat flow. The flexibility/stiffness matrix approach is used such that the crack problem reduces to solving two sets of singular integral equations. Numerical values of the crack tip stress-intensity factors are obtained for various crack size, crack location, crack surface insulation, fiber volume fraction and orientation angles. The results are displayed graphically.     

A formulation of anisotropic continuum elastoplasticity at finite strains. Part I: Modelling

A constitutive model for anisotropic elastoplasticity at finite strains is developed together with its numerical implementation. An anisotropic elastic constitutive law is described in an invariant setting by use of structural tensors and the elastic strain measure C_e. The elastic strain tensor as well as the structural tensors are assumed to be invariant in relation to superimposed rigid body rotations. An anisotropic Hill-type yield criterion, described by a non-symmetric Eshelby-like stress tensor and further structural tensors, is developed, where use is made of representation theorems for functions with non-symmetric arguments. The model also considers non-linear isotropic hardening. Explicit results for the specific case of orthotropic anisotropy are given. The associative flow rule is employed and the features of the inelastic flow rule are discussed in full. It is shown that the classical definition of the plastic material spin is meaningless in conjunction with the present formulation. Instead, the study motivates an alternative definition, which is based on the demand that such a quantity must be dissipation-free, as the plastic material spin is in the case of isotropy. Equivalent spatial formulations are presented too. The full numerical treatment is considered in Part II.     

基于非关联流动规律的Gotoh屈服准则的参数确定方法

屈服准则对板料成形过程的理论解析、工艺优化和有限元模拟有着重要的影响. 通过提高屈服准则的各向异性表征能力, 可以确保成形过程的可靠性及实际预测的准确性. 本文基于非关联流动法则, 给出了Gotoh屈服准则一套全新的参数求解方法. 在结合常用屈服准则并考虑流动规律的基础上, 分别以5754O铝合金、DP980先进高强钢和SAPH440结构钢作为研究对象, 进行了不同加载路径下各向异性变形行为的预测. 根据Gotoh屈服准则推导的屈服函数、塑性势函数以及基于关联流动的理论函数计算出屈服应力和各向异性指数$r$值随加载角度的分布趋势, 进而针对平面应力状态的屈服轨迹展开分析, 验证了不同屈服准则和流动规律对各向异性屈服行为的预测精度. 理论与实验数据对比结果表明: 不同屈服准则针对同种板料在流动规律一致的情形下其表征各向异性的能力有显著差异; 相同屈服准则基于不同流动规律其表征能力也具有明显差别. 基于非关联流动的屈服准则能极大地提高精度, 各向异性表征能力显著加强. 相关结果能够为各向异性屈服准则在塑性成形领域的实际应用方案提供重要参考.      

Asymmetrical dynamic fracture model of bridging fiber pull-out of unidirectional composite materials

An elastic analysis of an internal crack with bridging fibers parallel to the free surface in an infinite orthotropic anisotropic elastic plane is studied, and asymmetrical dynamic fracture model of bridging fiber pull-out of unidirectional composite materials is presented for analyzing the distributions of stress and displacement with the internal asymmetrical crack under the loading conditions of an applied non-homogenous stress and the traction forces on crack faces yielded by the bridging fiber pull-out model. Thus the fiber failure is ascertained by maximum tensile stress, the fiber ruptures and hence the crack propagation should also appear in the modality of self-similarity. The formulation involves the development of a Riemann-Hilbert problem. Analytical solution of an asymmetrical propagation crack of unidirectional composite materials under the conditions of two increasing loads given is obtained, respectively. In terms of correlative material properties, the variable rule of dynamic stress intensity factor was depicted very well. After those analytical solutions were utilized by superposition theorem, the solutions of arbitrary complex problems could be gained.