各向异性复合材料中界面屈服/脱粘对裂纹扩展的影响

本研究针对层状复合材料中正交于层合面的裂纹,研究裂纹前方层状界面发生屈服或脱粘现象对该裂纹前沿应力场的扰动。通过利用叠加原理,借用滑移型位错密度表征界面的屈服或脱粘。利用Chebyshev数值积分法求解相应的位错密度的奇异积分方程,得到沿界面屈服/脱粘区域的位错密度分布及裂端区应力场。结果表明,若层状复合材料界面为发生屈服或脱粘,将减弱独立层裂尖的应力奇异性,进而抑制独立层中裂纹的扩展。     

Cracking characteristics of mixed mode dislocations near a lip-like mode crack

This work is concerned with the cracking characteristics of mixed mode dislocations near a lip-like mode crack, stress intensity and strain energy density factor are obtained by using conformal mapping, singularity analysis and Cauchy integrals. Shielding effect generated by screw dislocation near a lip-like mode crack decreases with the increment of the distance between screw dislocation and crack tip. Larger distance between two faces of the crack leads to the shielding effect waning. The strain energy density factor of mode III decreases with the increment of the distance between dislocation and crack tip. Larger distance between two faces of lip-like mode crack also leads to the strain energy density factor waning and encourages crack initiation; the shielding effects generated by edge dislocation near the crack decrease with the increment of the distance between edge dislocation and crack tip.     

Crack growth resistance of shape memory alloys by means of a cohesive zone model

Crack growth resistance of shape memory alloys (SMAs) is dominated by the transformation zone in the vicinity of the crack tip. In this study, the transformation toughening behavior of a slowly propagating crack in an SMA under plane strain conditions and mode I deformation is numerically investigated. A small-scale transformation zone is assumed. A cohesive zone model is implemented to simulate crack growth within a finite element scheme. Resistance curves are obtained for a range of parameters that specify the cohesive traction-separation constitutive law. It is found that the choice of the cohesive strength t₀ has a great influence on the toughening behavior of the material. Moreover, the reversibility of the transformation can significantly reduce the toughening of the alloy. The shape of the initial transformation zone, as well as that of a growing crack is determined. The effect of the Young's moduli ratio of the martensite and austenite phases is examined.     

压电双材料界面钝裂纹附近螺型位错的屏蔽效应与发射条件

研究了压电双材料界面钝裂纹附近螺型位错的屏蔽效应与发射条件.应用保角变换技术,得到了复势函数与应力场的封闭形式解,讨论了位错方位、双材料电弹常数及裂纹钝化程度对位错屏蔽效应和发射条件的影响.结果表明,Burgers矢量为正的螺型位错可以降低界面钝裂纹尖端的应力强度因子(屏蔽效应),屏蔽效应随位错方位角及位错与裂纹尖端距离的增大而减弱,压电双材料中螺型位错对裂纹的屏蔽效应强于相应弹性双材料中螺型位错对裂纹的屏蔽效应;位错发射所需的临界无穷远加载或电位移随位错方位角及裂纹钝化程度的增加而增大;最可能的位错发射角度为零度即位错最可能沿裂纹前方发射.论文解答的特殊情况与已有文献一致.     

A multiscale model for the ductile fracture of crystalline materials

In this paper, a multiscale model that combines both macroscopic and microscopic analyses is presented for describing the ductile fracture process of crystalline materials. In the macroscopic fracture analysis, the recently developed strain gradient plasticity theory is used to describe the fracture toughness, the shielding effects of plastic deformation on the crack growth, and the crack tip field through the use of an elastic core model. The crack tip field resulting from the macroscopic analysis using the strain gradient plasticity theory displayes the 1/2 singularity of stress within the strain gradient dominated region. In the microscopic fracture analysis, the discrete dislocation theory is used to describe the shielding effects of discrete dislocations on the crack growth. The result of the macroscopic analysis near the crack tip, i.e. a new K-field, is taken as the boundary condition for the microscopic fracture analysis. The equilibrium locations of the discrete dislocations around the crack and the shielding effects of the discrete dislocations on the crack growth at the microscale are calculated. The macroscopic fracture analysis and the microscopic fracture analysis are connected based on the elastic core model. Through a comparison of the shielding effects from plastic deformation and the discrete dislocations, the elastic core size is determined.     

Dislocation behaviours ahead of crack tip

In this paper, a unified mechanics model for dislocation nucleation, emission and dislocation free zone is proposed based on the Peierls framework. Three regions are identified ahead of the crack tip. The emitted dislocations within the plastic zone in the form of an inverse pile up are treated as discrete elastic edge dislocations. Between that zone and the cohesive zone immediately ahead of the crack tip, there is a dislocation free zone. With the stress field and the dislocation density field in the cohesive zone, respectively, expressed in the first and second Chebyshev polynomial series, and the opening and slip displacements in trigonometric series, a set of nonlinear governing equations are obtained which take into account for the interaction between the emitted dislocations and cohesive zone and the nonlinear interaction between sliding displacement and the opening displacement. After discretization, the governing equations are transformed into a set nonlinear algebraic equations which are solved with Newton-Raphson Method. The results of calculation for pure shearing and combined tension and shear loading after dislocation emission are given in detail. Finally, the process of dislocation nucleation and emission on a pair of symmetric slip planes of angle α with respect to the crack plane under pure mode I load is analysed. The equilibrium positions and the number of emitted dislocation are determined. Several possible competition behaviors of dislocation emission vs cleavage are revealed.     

Arbitrarily oriented crack near interface in piezoelectric bimaterials

Arbitrarily oriented crack near interface in piezoelectric bimaterials is considered. After deriving the fundamental solution for an edge dislocation near the interface, the present problem can be expressed as a system of singular integral equations by modeling the crack as continuously distributed edge dislocations. In the paper, the dislocations are described by a density function defined on the crack line. By solving the singular integral equations numerically, the dislocation density function is determined. Then, the stress intensity factors (SIFs) and the electric displacement intensity factor (EDIF) at the crack tips are evaluated. Subsequently, the influences of the interface on crack tip SIFs, EDIF, and the mechanical strain energy release rate (MSERR) are investigated. The J-integral analysis in piezoelectric bimaterals is also performed. It is found that the path-independent of J₁-integral and the path-dependent of J₂-integral found in no-piezoelectric bimaterials are still valid in piezoelectric bimaterials.     

内聚力模型裂纹问题分析的解析奇异单元

内聚力模型已经被广泛应用于需要考虑断裂过程区的裂纹问题当中,然而常用的数值方法应用于分析内聚力模型裂纹问题时还存在着一些不足,比如不能准确地给出断裂过程区的长度、需要网格加密等.为了克服这些缺点,论文构造了一个新型的解析奇异单元,并将之应用于基于内聚力模型的裂纹分析当中.首先将虚拟裂纹表面处的内聚力用拉格拉日插值的方法近似表示为多项式形式,而多项式表示的内聚力所对应的特解可以被解析地给出.然后利用一个简单地迭代分析,基于内聚力模型的裂纹问题就可以被模拟出来了.最后,给出二个数值算例来证明论文方法的有效性.     

Interaction between an interface crack and subinterface microcracks in metal/piezoelectric bimaterials

The J-integral analysis is presented for the interaction problem between a semi-infinite interface crack and subinterface matrix microcracks in dissimilar anisotropic materials. After deriving the fundamental solutions for an interface crack subjected to different loads and the fundamental solutions for an edge dislocation beneath the interface, the interaction problem is deduced to a system of singular integral equations with the aid of a superimposing technique. The integral equations are then solved numerically and a conservation law among three values of the J-integral is presented, which are induced from the interface crack tip, the microcracks and the remote field, respectively. The conservation law not only provides a necessary condition to confirm the numerical results derived, but also reveals that the microcrack shielding effect in such materials could be considered as a redistribution of the remote J-integral. It is this redistribution that does lead to the phenomenological shielding effect.     

Exact solutions for piezoelectric screw dislocations near two asymmetrical edge cracks emanating from an elliptical hole

The interaction between piezoelectric screw dislocations and two asymmetrical interfacial cracks emanating from an elliptic hole under combined mechanical and electric load at infinity is dealt with. The closed-form solutions are derived for complex potentials and generalized stress fields. In the limiting cases, some well-known results can be obtained from the present solutions. Moreover, some new exact solutions are shown. The stress intensity factor and the energy release rate at the right tip due to a screw dislocation near the right interfacial crack are also calculated. The results show that the shielding effect of dislocation on crack expanding decreases with the increase in dislocation azimuth angle and the distance between the dislocation and the crack tip, and the repulsion acting on the dislocation from the other half plane demotes crack propagation. The increasing of the length of the other crack promotes crack growth, but the increasing of the minor semi-axis demotes it.     

Subsonic and intersonic mode II crack propagation with a rate-dependent cohesive zone

A recent experimental study has demonstrated the attainability of intersonic shear crack growth along weak planes in otherwise homogeneous, isotropic, linear elastic solids subjected to remote loading conditions (Rosakis et al., Science 284 (5418) (1999) 1337). The relevant experimental observations are summarized briefly here and the conditions governing the attainment of intersonic crack speeds are examined. Motivated by experimental observations, subsonic and intersonic mode II crack propagation with a rate-dependent cohesive zone is subsequently analyzed. A cohesive law is assumed, wherein the cohesive shear traction is either a constant or varies linearly with the local sliding rate. Complete decohesion is assumed to occur when the crack tip sliding displacement reaches a material-specific critical value. Closed form expressions are obtained for the near-tip fields. With a cohesive zone of finite size, it is found that the dynamic energy release rate is finite through out the intersonic regime. Crack tip stability issues are addressed and favorable speed regimes are identified. The influence of shear strength of the crack plane and of a rate parameter on crack propagation behavior is also investigated. The isochromatic fringe patterns predicted by the analytical solution are compared with the experimental observations of Rosakis et al. (1999) and comments are made on the validity of the proposed model.     

一维六方准晶中平行螺型位错与有限长裂纹的相互作用

论文采用解析函数理论,利用Muskhelishvili著名的复变函数方法,给出了无限大一维六方准晶中多条平行螺型位错与有限长直裂纹相互作用的解析解,并得到了裂纹尖端的场强度因子和作用在位错上的像力.数值算例讨论了镜像力、场强度因子随位错位置的变化,及位错的位置与分布对广义应力场的影响.     

Simulation of local instabilities during crack propagation in the ductile–brittle transition region

The crack propagation for a cohesive zone model within an elastic-plastic material under small-scale yielding conditions is simulated numerically. The resulting crack growth resistance curves show local instabilities, so-called pop-ins even for homogeneous cohesive properties if the cohesive strength lies sufficiently close to the maximum stress of the corresponding blunting solution. The formation of secondary cracks and unloading zones in front of the actual crack tip is identified as the underlying mechanism. It is found that the shape of the cohesive law has a considerable influence on the crack arrest behavior. Furthermore, requirements to the mesh resolution are derived.     

Comparative analysis of extrinsic and intrinsic cohesive models of dynamic fracture

A comparative analysis of intrinsic and extrinsic cohesive models has been performed for the case of spontaneous and steady-state dynamic crack propagation. Spontaneous crack propagation was simulated using a spectral form of the elastodynamic boundary integral equation, while steady-state solutions were obtained by numerically integrating the governing Cauchy singular integral equation. Spontaneous crack propagation results showed that intrinsic models are less numerically stable than the extrinsic ones. Under steady-state propagation conditions, some intrinsic cohesive models lead to unrealistic results as the crack opening velocity becomes negative at the cohesive zone tip. By imposing a positive crack opening acceleration at the cohesive zone tip, the envelope of the required minimum initial strength has been calculated.     

位错、裂纹与断裂

本文评述近些年来在位错、裂纹与断裂方面的研究进展.内容主要包括:①以屈服强度温度依赖性为例,讨论力学性能的结构敏感性,使其对位错芯结构与塑性之间的关系有一定的认识;②纯弹性断裂总是极少的,因之在加载应力作用下,裂纹顶端总是存在着或大或小的塑性区,存在着裂纹与位错的交互作用.这里以我们自己的工作为主,重点讨论裂端位错过程的研究结果.      

Finite deformation analysis of mechanism-based strain gradient plasticity: torsion and crack tip field

A finite deformation theory of mechanism-based strain gradient (MSG) plasticity is developed in this paper based on the Taylor dislocation model. The theory ensures the proper decomposition of deformation in order to exclude the volumetric deformation from the strain gradient tensor since the latter represents the density of geometrically necessary dislocations. The solution for a thin cylinder under large torsion is obtained. The numerical method is used to investigate the finite deformation crack tip field in MSG plasticity. It is established that the stress level around a crack tip in MSG plasticity is significantly higher than its counterpart (i.e. HRR field) in classical plasticity.     

A solution to the thermo-elastic interface crack branching in dissimilar anisotropic bi-material media

An interface crack or delamination may often branch out of the interface in a laminated composite due to thermal stresses developing around the delamination/crack tip when the media is exposed to heat flow induced by environmental events such as a sudden short-duration fire. In this paper, the thermo-elastic problem of interface crack branching in dissimilar anisotropic bi-media is studied by using the theory of Stroh’s dislocation formalism, extended to thermo-elasticity in matrix notation. Based on the complex variable method and the analytical continuation principle, the thermo-elastic interface crack/delamination problem is examined and a general solution in compact form is derived for dissimilar anisotropic bi-media. A set of Green’s functions is proposed for the dislocations (conventional dislocation and thermal dislocation/heat vortex) in anisotropic bi-media. These functions may be more suitable than those which have appeared in the literature on addressing thermo-elastic interface crack branching in dissimilar anisotropic bi-materials. Using the contour integral method, a closed form solution to the interaction between the dislocations and the interface crack is obtained. Within the scope of linear fracture mechanics, the thermo-elastic problem of interface crack branching is then solved by modelling the branched portion as a continuous distribution of dislocations. The influence of thermal loading and thermal properties on the branching behavior is examined, and criteria for predicting interface crack branching are suggested, based on the extensive numerical results from the study of various cases.     

Singularity-free dislocation dynamics with strain gradient elasticity

The singular nature of the elastic fields produced by dislocations presents conceptual challenges and computational difficulties in the implementation of discrete dislocation-based models of plasticity. In the context of classical elasticity, attempts to regularize the elastic fields of discrete dislocations encounter intrinsic difficulties. On the other hand, in gradient elasticity, the issue of singularity can be removed at the outset and smooth elastic fields of dislocations are available. In this work we consider theoretical and numerical aspects of the non-singular theory of discrete dislocation loops in gradient elasticity of Helmholtz type, with interest in its applications to three dimensional dislocation dynamics (DD) simulations. The gradient solution is developed and compared to its singular and non-singular counterparts in classical elasticity using the unified framework of eigenstrain theory. The fundamental equations of curved dislocation theory are given as non-singular line integrals suitable for numerical implementation using fast one-dimensional quadrature. These include expressions for the interaction energy between two dislocation loops and the line integral form of the generalized solid angle associated with dislocations having a spread core. The single characteristic length scale of Helmholtz elasticity is determined from independent molecular statics (MS) calculations. The gradient solution is implemented numerically within our variational formulation of DD, with several examples illustrating the viability of the non-singular solution. The displacement field around a dislocation loop is shown to be smooth, and the loop self-energy non-divergent, as expected from atomic configurations of crystalline materials. The loop nucleation energy barrier and its dependence on the applied shear stress are computed and shown to be in good agreement with atomistic calculations. DD simulations of Lomer–Cottrell junctions in Al show that the strength of the junction and its configuration are easily obtained, without ad-hoc regularization of the singular fields. Numerical convergence studies related to the implementation of the non-singular theory in DD are presented.     

纳米晶体铝在单轴拉伸下的位错堆积对微裂纹扩展的影响

摘要：针对纳米晶体材料,研究了单轴拉伸载荷作用下纳米晶体铝中的裂纹与裂纹尖端发射的位错所形成的滑移面之间的相互作用。通过分布位错法,将裂纹和滑移面等效为均匀分布的连续位错,获得了裂纹面上应力场。并引入裂纹尖端的无位错区,研究了裂纹尖端无位错区对微裂纹的萌生和主裂扩展之间的影响。结果表明,不考虑裂纹尖端无位错区时,裂纹长度较短,会先在晶界处形成微裂纹,主裂纹较长时,主裂纹会直接穿晶扩展。滑移面与裂纹尖端夹角较大时,会增加裂纹尖端发射的位错个数,从而抑制主裂纹的扩展。考虑裂纹尖端无位错区时,无位错区先于晶界处出现微裂纹,通过主裂纹与微裂纹之间位错的相互发射,导致裂纹与尖端处微裂纹汇合,有效加速了主裂纹的扩展。