Elastic and elastic–plastic singular fields around a crack-tip in particulate-reinforced composites with progressive debonding damage

This paper deals with elastic and elastic–plastic singular fields around a crack-tip in particulate-reinforced composites with debonding damage of particle-matrix interface. Numerical analyses are carried out on a crack-tip field in elastic-matrix and elastic–plastic-matrix composites reinforced with elastic particles, using a finite element method developed based on an incremental damage theory of particulate-reinforced composites. A particle volume fraction and interfacial strength between particles and matrix of the composites are parametrically changed. In the elastic-matrix composites, a unique elastic singular field is created on the complete damage zone in the vicinity of a crack-tip in addition to the conventional elastic singular field on the no damage zone. The macroscopic stress level around a crack-tip is reduced by the debonding damage while the microscopic stress level of the matrix remains unchanged. In the elastic–plastic-matrix composites, the damage zone develops in addition to the plastic zone due to matrix plasticity, and both the macroscopic and microscopic stress revels around a crack-tip are reduced by the debonding damage. It is concluded from the numerical results that the toughening due to damage could be expected in the elastic–plastic-matrix composites, while it is questionable in the elastic-matrix composites.     

A kinematic elastic nonshakedown theorem for implicit standard materials

Summary  Criteria for a priori recognition of the type of steady-state response induced by cyclic loads and prediction whether a structure will shakedown elastically or not, without the necessity of performing a step-by-step full analysis, have considerable importance. Melan and Koiter theorems provide criteria that guarantee whether elastic shakedown occurs or not under cyclic loads in case of perfect plasticity. However, there remain some aspects of the shakedown theory which deserve further study. One of these, concerned with more realistic nonassociative elastic–plastic constitutive material models, allowing for nonlinear kinematic and isotropic hardening suitable to describe the cyclic plastic behaviour of metallic materials, has strong motivation. Koiter's elastic nonshakedown theorem is reconsidered here, with the objective of extending it to the de Saxcé's implicit standard material class, which contains a wide class of nonassociative elastic–plastic material behaviours. Shakedown analysis is formulated by a kinematic approach based on the plastic accumulation mechanism concept due to Polizzotto. A sufficient condition for elastic nonshakedown and a distinct necessary condition are established. Then, an upper bound to the shakedown multiplier is evaluated. Received 15 February 2001; accepted for publication 18 October 2001     

Elastic–plastic contact force history and response characteristics of circular plate subjected to impact by a projectile

A new elastic–plastic impact–contact model is proposed in this paper. By adopting the principle of minimum acceleration for elastic–plastic continue at finite deformation, and with the aid of finite difference method, the proposed model is applied in the problem of dynamic response of a clamped thin circular plate subjected to a projectile impact centrally. The impact force history and response characteristics of the target plate is studied in detail. The theoretical predictions of the impact force and plate deflection are in good agreements with those of LDA experimental data. Linear expressions of the maximum impact force/transverse deflection versus impact velocity are given on the basis of the theoretical results. The project supported by the National Natural Science Foundation of China (10532020).     

Elastic–plastic behavior of a semicircular frame being pressed against a rigid plane

As a simplified structural model, a semicircular frame is used to study the crashworthiness behavior of an aircraft fuselage. The quasi-static large elastic-plastic deformation of a semicircular frame in the process of its being pressed against a rigid ground is analyzed. First, based on the linear elastic assumption, the quasi-static large deformation contact process of the frame can be divided into three phases, i.e., point contact, line contact and post-buckling. By means of a shooting method, the relations between the displacement and contact force as well as the distribution of bending moment in the three phases are obtained. Then, by assuming an elastic, perfectly-plastic moment-curvature relationship for the semi-circular frame, the contact process is analyzed in detail to reveal the plastic collapse mechanism, the traveling of plastic hinge and the force-displacement relationship. In order to verify the analysis, a preliminary experiment was conducted, in which two types of half rings with clamped ends were pressed by a rigid plate. In addition, a numerical simulation is also conducted by employing ABAQUS to analyze both rectangular cross-sectional beam and I-beam. Finally, the theoretical predictions are compared with the experimental results and numerical solutions, showing that the elastic-plastic analysis can predict the contact process very well.     

On the yield functions in Lagrangian stress space of elastic–plastic materials with internal constraints

In the present paper it is shown that the elastic range in the second Piola–Kirchhoff stress space can be chosen in a hyperplane which is through the origin of Lagrangian stress space and perpendicular to the normal of the constraint manifold at the plastic configuration, if the determinate stress response of the elastic–plastic material with simple internal constraints with some condition is correctly chosen, otherwise, it is in general in a hypersurface and the normal flow rule by Il yushin’s postulate will have an indeterminate part. The choice of determinate stress response is probable because of its indeterminacy. Therefore the yield function should be a function of the second Piola–Kirchhoff stress lying in the hyperplane so that it is more simple and the back stress as the geometric center of the elastic range in general is inside the elastic range. Finally some examples are concerned. The project supported by the National Natural Science Foundation of China (10272055).     

A phenomenological model of fatigue crack growth in perfectly plastic infinite plates under completely reversed uniaxial loading

The fatigue failure of a thin infinite center-cracked plate under completely reversed uniaxial loading is considered. A two-stage fatigue crack model including the incubation and crack propagation stages is constructed. The stress distribution in the vicinity of the crack tip is described using the concept of a conventional elastic crack. The crack-tip plastic zone is simulated by a Dugdale thin plastic zone, and the condition for the movement of the failure front is given by criteria of damage mechanics. It is shown that the fatigue crack growth rate in perfectly plastic materials with a plastic zone of constant length is a power-law function of the stress intensity factor range. This relationship is quadratic when the length of the plastic zone is not constant Published in Prikladnaya Mekhanika, Vol. 41, No. 12, pp. 116–127, December 2005.     

M积分与夹杂/缺陷弹性模量的显式关系

M积分在材料构型力学中表征着缺陷自相似扩展的能量释放率,而有效弹性模量下降量在传统损伤力学中是一个具有内变量属性的损伤参数. 探讨了两者之间的特定关系,以此为材料构型力学与损伤力学搭建桥梁.借助穆斯海里什维利（Muskhelishvili）复势函数方法获取无限大弹性平面含圆形夹杂的弹性场解,根据M 积分的复势函数解析表达式得到M 积分与夹杂弹性模量的显式表达式. 随后通过有限元分析,对含复杂缺陷群的弹塑性材料进行数值模拟,结果表明内部缺陷区域的有效弹性模量下降与M 积分存在着特定关系. 基于此,提出利用材料构型力学中的外变量参数（M 积分）来替代损伤力学中的内变量参数（弹性模量下降量）描述材料的缺陷演化.      

On the Measurement and Prediction of the Out-of-Plane Displacement Surrounding Cold-Expanded Holes

Experimental measurements of the out-of-plane displacement surrounding cold-expanded holes in a 6082-T6 aluminum alloy were made with a 3D optical scanner using the technique of the encoded light-pattern projection in white light. The measured surface profiles have shown the thickness discontinuity along the hole edge due to the effect of the split in the sleeve. An analytical–numerical solution of the out-of-plane displacement is presented based on existing analytical models. Thus, the results given by the analytical model were then compared with the experimental data and with a finite element (FE) model that simulates the cold-expansion process. The location of the elastic–plastic boundary was estimated as the point at which no change in thickness was observed; a good agreement was found in the comparison of measured, FE and analytical results. The measured surface profiles agreed with those predicted by the FE model and analytical solution. The proposed experimental approach can be used together with FE analysis for predicting the radial and circumferential residual stresses in cold-expanded hole. It is quite versatile and can also be used as quality-control technique in the manufacturing processes of cold-expanded holes.     

Sub-Grain Scale Digital Image Correlation by Electron Microscopy for Polycrystalline Materials during Elastic and Plastic Deformation

Damage during loading of polycrystalline metallic alloys is localized at or below the scale of individual grains. Quantitative assessment of the heterogeneous strain fields at the grain scale is necessary to understand the relationship between microstructure and elastic and plastic deformation. In the present study, digital image correlation (DIC) is used to measure the strains at the sub-grain level in a polycrystalline nickel-base superalloy where plasticity is localized into physical slip bands. Parameters to minimize noise given a set speckle pattern (introduced by chemical etching) when performing DIC in a scanning electron microscope (SEM) were adapted for measurements in both plastic and elastic regimes. A methodology for the optimization of the SEM and DIC parameters necessary for the minimization of the variability in strain measurements at high spatial resolutions is presented. The implications for detecting the early stages of damage development are discussed.     

Explosion Dynamics in Saturated Rocks and Solids

Non-elastic pore deformations and crack propagations are the principal causes of dynamic damage in rocks and soils. In the case of downhole blasting from wellbores, these two mechanisms compete with each other. Therefore, to carry out a mechanical analysis of rock blasting, a sufficiently complete model that takes these various mechanisms into account has to be developed. To address this issue, this paper proposes the use of an elastic–plastic model, which includes a yield condition with a non-associated plastic flow rule, the effects of pore fluid saturation, and a brittle failure criterion under extension. The results presented in this paper describe underground explosions with spherical motion (cavity growth under the internal pressure of detonated gases without leakage into the formation), typical for oil or water reservoirs. The governing equations are written in a Cartesian system of coordinates for the case of spatial dynamic medium deformation. For this case, Cartesian coordinates are more convenient than spherical coordinates because they avoid numerical difficulties connected with the non-divergent terms of the non-linear form of the Biot–Frenkel equations. The numerical method uses the Wilkins approach, which has been generalized for the model described in this paper. The dilatancy of the material during plastic deformation is neglected for simplicity. The numerical results show that, when using typical parameters for relatively “soft” porous skeleton, the plastic flow overcomes the brittle failure. An extension zone only appears near the cavity. The results also show the presence of the two Biot P-waves. The second Biot wave, however, is only seen in the case of an extremely high permeability rock. Furthermore, in the case of the first Biot wave, the saturating liquid and the solid skeleton particles are moving with different velocities in a 100 darcy rock and with the same velocity in a 0.01 darcy rock. Calculated radial particle velocities as a function of the scaled radius are close to measured velocities in rigid dense media but are larger than measured ones in clays. It is suggested that the difference is due to different levels of water saturation, assumed full saturation in the calculation, partial saturation in the experiments.     

M积分与夹杂/缺陷弹性模量的显式关系简

M积分在材料构型力学中表征着缺陷自相似扩展的能量释放率,而有效弹性模量下降量在传统损伤力学中是一个具有内变量属性的损伤参数. 探讨了两者之间的特定关系,以此为材料构型力学与损伤力学搭建桥梁.借助穆斯海里什维利（Muskhelishvili）复势函数方法获取无限大弹性平面含圆形夹杂的弹性场解,根据M 积分的复势函数解析表达式得到M 积分与夹杂弹性模量的显式表达式. 随后通过有限元分析,对含复杂缺陷群的弹塑性材料进行数值模拟,结果表明内部缺陷区域的有效弹性模量下降与M 积分存在着特定关系. 基于此,提出利用材料构型力学中的外变量参数（M 积分）来替代损伤力学中的内变量参数（弹性模量下降量）描述材料的缺陷演化.     

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.     

On the plastic wave speeds in rate-independent elastic–plastic materials with anisotropic elasticity

This paper presents a theoretical study of the speeds of plastic waves in rate-independent elastic–plastic materials with anisotropic elasticity. It is shown that for a given propagation direction the plastic wave speeds are equal to or lower than the corresponding elastic speeds, and a simple expression is provided for the bound on the difference between the elastic and the plastic wave speeds. The bound is given as a function of the plastic modulus and the magnitude of a vector defined by the current stress state and the propagation direction. For elastic–plastic materials with cubic symmetry and with tetragonal symmetry, the upper and lower bounds on the plastic wave speeds are obtained without numerically solving an eigenvalue problem. Numerical examples of materials with cubic symmetry (copper) and with tetragonal symmetry (tin) are presented as a validation of the proposed bounds. The lower bound proposed here on the minimum plastic wave speed may also be used as an efficient alternative to the bifurcation analysis at early stages of plastic deformation for the determination of the loss of ellipticity.     

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.     

Modelling and FE-simulation of plastic ductile damage evolution in plates and shells

Summary  The present paper reports on modelling and numerical simulation of thin-walled structures close to failure taking into consideration the effects of both geometrical and physical nonlinearity. The approach accounts for finite displacements and rotations, and the material model adopted includes elastic–plastic behaviour, isotropic and kinematic hardening, and ductile damage. Particular attention is paid to the problems of localised damage, damage progression and final collapse of the structure. Numerical simulation of the nonlinear response of bars, plates and shells to quasistatic monotonic and variable loading illustrates how material damage affects the load-carrying behaviour of structural components. Received 28 November 1999; accepted for publication 29 March 2000     

基于M积分的脆性材料微缺陷等效损伤面积/体积表征

随着脆性材料在工程中的广泛应用,对脆性材料中微缺陷进行统一的损伤水平标定,具有重要的科学研究和工程应用价值. 本研究提出一种基于M积分的材料等效损伤面积/体积标定方法, 以具有相同M积分值的圆孔面积或球孔体积来标定 复杂微缺陷构型的损伤水平,从而实现不同类型微缺陷真实损伤水平的统一表征. 首先,基于Lagrangian能量密度函数推导 了M积分定义式,并简述了M积分的物理意义,基于域积分方法实现二维/三维M积分的数值计算. 随后,提出了基 于M积分的材料缺陷等效损伤面积/体积标定方法,以圆孔面积/球孔体积来标定复杂微缺陷材料系统的损伤水平. 最后, 针对单轴拉伸载荷作用下的二维/三维脆性体含不同缺陷构型,具体计算了椭圆孔、裂纹以及双裂纹、双孔洞、裂纹和孔洞干 涉等复杂缺陷构型情况下的等效损伤面积/体积,并详细分析了缺陷之间的干涉效应及影响因素. 本研究旨在基于M积分等 效方法量化脆性材料中各类微缺陷造成的损伤程度,实现损伤等级标定,有益于工程材料及结构的损伤容限设计及完整性评估.      

帽型试样动态绝热剪切破坏演化分析

利用分离式霍普金森压杆加载Ta2钛合金扁平闭合帽形受迫剪切试样,结合数字图像相关法和“冻结”试样的微观金相观察,研究剪切区剪切应变的演化、绝热剪切带形成条件等。结果显示：受迫剪切试样在动态加载过程,剪切区剪切应变不断集中,形成绝热剪切带,裂纹沿绝热剪切带发展;随加载率提高,绝热剪切起始临界应变减小;进一步利用数字图像相关法DIC场应变分析及金相微观观测对比,利用卸载回复特性对绝热剪切带起始临界条件进行了讨论,计算的绝热剪切带起始时温升仅为86℃。材料软化可能不是绝热剪切带起始的控制条件,相反是由于绝热剪切带形成造成的应变高度集中发展导致温度急剧升高。     

Study on magneto–elastic–plastic deformation characteristics of ferromagnetic rectangular plate with simple supports

In this paper, the magnetic-elastic-plastic deformation behavior is studied for a ferromagnetic plate with simple supports. The perturbation formula of magnetic force is first derived based on the perturbation technique, and is then applied to the analysis of deformation characteristics with emphasis laid on the analyses of modes, symmetry of deformation and influences of incident angle of applied magnetic field on the plate deformation. The theoretical analyses offer explanations why the configuration offer- romagnetic rectangular plate with simple supports under an oblique magnetic field is in-wavy type along the x-direction, and why the largest deformation of the ferromagnetic plate occurs at the incident angle of 45°for the magnetic field. A numerical code based on the finite element method is developed to simulate quantitatively behaviors of the nonlinearly coupled multi-field problem. Some characteristic curves are plotted to illustrate the magneto--elastic-plastic deflections, and to reveal how the deflections can be influenced by the incident angle of applied magnetic field. The deformation characteristics obtained from the numerical simulations are found in good agreement with the theoretical analyses.     

An energy-based crack growth criterion for modelling elastic–plastic ductile fracture

A crack growth criterion is derived based on the Griffith energy concept and the cohesive zone model for modelling fracture in elastic–plastic ductile materials. The criterion is implemented in the finite element context by a virtual crack extension technique. An automatic modelling of the ductile fracture process is realised by combining a local remeshing procedure and the criterion. The validity of the derived criterion is examined by modelling a compact tension specimen.     

Analysis of the plastic zone at the corner point of interface

A symmetric problem of elasticity is formulated to analyze the plastic zone at the corner point of the interface between two isotropic media. The piecewise-homogeneous isotropic body with an interface in the form of angle sides consists of different elastic parts joined by a thin elastoplastic layer. The plastic zone is modeled by discontinuity lines of tangential displacement, which are located at the interface. The exact solution of the problem is found using the Wiener–Hopf method and is then used to determine the length of the plastic zone. The stress at the corner point is analyzed Translated from Prikladnaya Mekhanika, Vol. 45, No. 2, pp. 59–69, February 2009.