双轴压缩载荷下裂纹动态传播的实验研究

本文研究了双轴压缩载荷下裂纹的起裂和非共面传播问题，用动态光弹性方法成功地记录下裂纹的动态扩展过程。获得了双轴压缩载荷下扩展裂纹的传播速度，临界开裂载荷以及开裂角第一系重要断裂参数。研究表明，在双轴压缩载荷下不闭合裂纹前缘应力场是一个混合型的奇异场，裂纹的扩展过程为稳态过程并且表现出显著的阶段性特点。     

EXPERIMENTAL STUDIES ON DYNAMIC PROPAGATION OF CRACKS UNDER BIAXIAL COMPRESSION

This paper deals with the initiation and non-coplaner propagation ofcracks under biaxial compression.The running process of the propagating cracks issuccessfully recorded with the method of dynamic photoclasticity.A series of importantfracture parameters such as cracks' propagating velocity,critical opening loading andinitial angle of cracks' growth under hiaxial compression is calculated.Results showthat the stress field in front of the non-closed crack tip under biaxial compression is amixed-mode singular field,the dynamic propagation of the new tensile cracks is juststable and shows distinct stages,but the unstable crack appears under high loads andresults in the rupture of the matrix.     

Dynamic crack deflection and penetration at interfaces in homogeneous materials: experimental studies and model predictions

We examine the deflection/penetration behavior of dynamic mode-I cracks propagating at various speeds towards inclined weak planes/interfaces of various strengths in otherwise homogeneous isotropic plates. A dynamic wedge-loading mechanism is used to control the incoming crack speeds, and high-speed photography and dynamic photoelasticity are used to observe, in real-time, the failure mode transition mechanism at the interfaces. Simple dynamic fracture mechanics concepts used in conjunction with a postulated energy criterion are applied to examine the crack deflection/penetration behavior and, for the case of interfacial deflection, to predict the crack tip speed of the deflected crack. It is found that if the interfacial angle and strength are such as to trap an incident dynamic mode-I crack within the interface, a failure mode transition occurs. This transition is characterized by a distinct, observable and predicted speed jump as well as a dramatic crack speed increase as the crack transitions from a purely mode-I crack to an unstable mixed-mode interfacial crack.     

On the higher order asymptotic analysis of a non-uniformly propagating dynamic crack along an arbitrary path

Transient mixed-mode elastodynamic crack growth along arbitrary smoothly varying paths is considered. Asymptotically, the crack tip stress field is square root singular with the angular variation of the singular term depending weakly on the instantaneous values of the crack tip speed and on the mode-I and mode-II stress intensity factors. However, for a material particle at a small distance away from the moving crack tip, the local stress field will depend not only on the instantaneous values of the crack tip speed and stress intensity factors, but also on the past history of these time dependent quantities. In addition, for cracks propagating along curved paths the stress field is also expected to depend on the nature of the curved crack path. Here, a representation of the crack tip fields in the form of an expansion about the crack tip is obtained in powers of radial distance from the tip. The higher order coefficients of this expansion are found to depend on the time derivative of crack tip speed, the time derivatives of the two stress intensity factors as well as on the instantaneous value of the local curvature of the crack path. It is also demonstrated that even if cracks follow a curved path dictated by the criterion K ₁₁ ^d =0, the stress field may still retain higher order asymmetric components related to non-zero local curvature of the crack path.     

Life prediction of thermally cracked railway wheels: Growth estimation of cracks with arbitrary shape

A computer program for railway wheel life prediction has been developed incorporating an elastoplastic calculation of both the residual stress field induced by shoe-braking and the superposed alternating rolling contact stresses. The possibilities of unstable propagation of surface cracks, crack arrest, fatigue propagation through a complex stress field and crack detention on reaching the threshold value have been accounted for. An approximate method for crack growth prediction is used providing bounds for instantaneous crack front position and estimations of crack shape within short computing times. Some examples of application for an UIC R7 wheel are presented.     

A domain independent integral expression for the crack extension force of a curved crack in three dimensions

An integral expression that is domain independent in curvilinear coordinates and compatible with zero divergence of Eshelby's (Phil. Trans. Roy. Soc. (London) 244 (1951) 87.) energy momentum tensor was obtained from the principle of virtual work. By applying Eshelby's definition of the force on a material defect a general expression of the crack extension force for a curved crack in three dimensions, here called the F-integral, was derived from the domain independent integral expression. The F-integral is given explicitly for a number of curved cracks and found to be in agreement with previously known solutions, when available. The influence of crack surface and crack front curvature upon the various forms of the F-integral is discussed. The F-integral presented in this work is a generalisation of the J-integral (Rice, J. Appl. Mech. 35 (1968) 379.) to curved cracks in orthogonal curvilinear coordinates.     

Stability of an advancing crack to small perturbation of its path

In this work we investigate the stability of a nominally straight two-dimensional quasistatically growing crack to a small perturbation of its path. Formulae for perturbations of stress intensity factors induced by slight deviation of the crack trajectory were developed by Movchan et al. (Int. J. Solids Struct. 35, 3419) Their solution is exploited to derive an equation for the perturbation of the crack path on the assumption that the crack advances in pure “opening” mode (i.e. local K_II=0). Various types of loading conditions are considered, including a cracked body loaded by a pair of point body forces and a crack whose faces are subjected to given tractions acting in the direction normal to the crack boundary. The body is also subjected to a remotely maintained uniaxial stress, aligned with the direction of the unperturbed crack. The loading is assumed to advance as the crack advances, to maintain the critical value of Mode I stress intensity factor. Numerical computations of possible crack paths have been performed, extending results on crack stability obtained by Cotterell and Rice (Int. J. Fract. 16, 155). The results show that in the case of loading by point body forces the stability of the crack path depends on the positions of the points of application of the applied forces and the magnitude of the applied stress acting parallel to the crack. There exists a critical value of this stress such that the crack path is stable for values less than critical and unstable otherwise. It is shown that the crack is always unstable in the case of point force tractions applied normal to the crack faces.     

压电薄板板边导电裂纹尖端的力电场分析

采用有限元方法，分析了压电薄板板边不同长度导电裂纹尖端的力电场分布规律，发现导电裂纹尖端的应力场和电场强度存在明显的集中和奇异现象，集中和奇异的程度与裂纹长度有关。而且，在裂纹延长线上分别存在两点，这里的应力和电场对裂纹长度不太敏感，总等于无裂纹时薄板的均匀应力和均匀电场强度；同时，还研究了导电裂纹尖端的应力强度因子和电场强度因子对裂纹长度的依赖关系，发现在线性本构的前提下，导电裂纹尖端的应力强度因子与电场强度因子之间具有近似的线性关系。     

Effect of electric fields on fracture behavior of ferroelectric ceramics

The asymptotic problem of a semi-infinite crack perpendicular to the poling direction in a ferroelectric ceramic subjected to combined electric and mechanical loading is analyzed to investigate effect of electric fields on fracture behavior. Electromechanical coupling induced by the piezoelectric effect is neglected in this paper. The shape and size of the switching zone is shown to depend strongly on the relative magnitude between the applied electric field and stress field as well as on the ratio of the coercive electric field to the yield electric field. A universal relation between the crack tip stress intensity factor and the applied intensity factors of stress and electric field under small-scale conditions is obtained from the solution of the switching zone. It is found that the ratio of the coercive electric field to the yield electric field plays a significant role in determining the enhancement or reduction of the crack tip stress intensity factor. The fracture toughness variation of ferroelectrics under combined electric and mechanical loading is also discussed.     

Partition of unity-based thermomechanical meshfree method for two-dimensional crack problems

We present a partition of unity-enriched element-free Galerkin method for thermoelastic two-dimensional crack problems. Therefore, the displacement field is enriched by the step enrichment. In the vicinity of the crack tip, the asymptotic branch enrichment functions commonly used in linear elastic fracture mechanics are employed. The same enrichment strategy is employed for the temperature field. Level set functions are used in order to model the crack surface. The accuracy of the method is demonstrated for three examples, one involves the crack propagation due to temperature and mixed traction-temperature loading conditions with complex curved crack paths.     

Brittle fracture dynamics with arbitrary paths. II. Dynamic crack branching under general antiplane loading

The dynamic propagation of a bifurcated crack under antiplane loading is considered. The dependence of the stress intensity factor just after branching is given as a function of the stress intensity factor just before branching, the branching angle and the instantaneous velocity of the crack tip. The jump in the dynamic energy release rate due to the branching process is also computed. Similar to the single crack case, a growth criterion for a branched crack is applied. It is based on the equality between the energy flux into each propagating tip and the surface energy which is added as a result of this propagation. It is shown that the minimum speed of the initial single crack which allows branching is equal to 0.39c, where c is the shear wave speed. At the branching threshold, the corresponding bifurcated cracks start their propagation at a vanishing speed with a branching angle of approximately 40°.     

Effect of electric displacement saturation on the hysteretic behavior of ferroelectric ceramics and the initiation and propagation of cracks in piezoelectric ceramics

This paper presents a computational investigation of a proposed simplified account for electric displacement saturation on the hysteretic behavior of initially unpoled ferroelectric ceramics as well as on the initiation and propagation of cracks in poled ferroelectric ceramics within the linear regime of piezoelectricity. For the latter case, experimental observations suggest an odd dependency of the onset of crack initiation in these brittle materials on the orientation of the applied electric field with respect to their poling direction which contradicts theoretical results which propose an even dependency of the energy release rate on the applied electric field within the framework of anisotropic linear piezoelectricity. Electric non-linearities arising at regions of inhomogeneities such as inclusions or at the crack tip are proposed in the literature to avoid this discrepancy. Electric displacement saturation is one such non-linear effect which is investigated in this work. A simplified account of this effect is proposed based on an exponential saturation model of the identified material parameters which can be related to this non-linearity. Its advantage over the superposition of a complex function onto the singular solution of a crack within the framework of linear piezoelectricity lies in the straightforward extension of the proposed approach to problems where no analytical solutions exist. This is outlined based on its incorporation into a rate-dependent ferroelectric model accounting for polarization switching as well as based on its incorporation into a finite element framework capable of simulating the initiation and propagation of cracks in piezoelectric ceramics through strong discontinuities in the displacement field and the electric potential. It is shown that besides the determination of the crack initiation onset also the crack propagation direction is influenced by the appearance of saturation zones arising at the crack tip normal to the polarization direction. The numerically obtained crack paths are found to be close to the experimentally reported results.     

Coupled solution for anisotropic piezoelectric materials with cracks

The coupled elastic and electric fields for anisotropic piezoelectric materials with electrically permeable cracks are analyzed by using Stroh formula in anisotropic elasticity. It is shown from the solution that the tangent component of the electric field strength and the normal component of the electric displacement along the faces of cracks are all constants, and the electric field intensity and electric displacement have the singularity of type (1/2) at the crack tip. The energy release rate for crack propagation depends on both the stress intensity factor and material constants. The electric field intensity and electric displacement inside electrically permeable cracks are all constants.     

Coplanar propagation paths of 3D cracks in infinite bodies loaded in shear

Bower and Ortiz, recently followed by Lazarus, developed a powerful method, based on a theoretical work of Rice, for numerical simulation of planar propagation paths of mode 1 cracks in infinite isotropic elastic bodies. The efficiency of this method arose from the need for the sole 1D meshing of the crack front. This paper presents an extension of Rice’s theoretical work and the associated numerical scheme to mixed-mode (2 + 3) shear loadings. Propagation is supposed to be channeled along some weak planar layer and to remain therefore coplanar, as in the case of a geological fault for instance. The capabilities of the method are illustrated by computing the propagation paths of cracks with various initial contours (circular, elliptic, rectangular, heart-shaped) in both fatigue and brittle fracture. The crack quickly reaches a stable, almost elliptic shape in all cases. An approximate but accurate analytic formula for the ratio of the axes of this stable shape is derived.     

Parametric instability of a Jeffcott rotor with rotationally asymmetric inertia and transverse crack

Both the rotationally asymmetric inertia and transverse crack frequently appear in the rotor system. The parametric excitations induced by this two features cause instability and severe vibration under certain operating conditions. Thus, the parametric instability of a Jeffcott rotor with asymmetric disk and open transverse crack is studied analytically. The vibration equations of four degrees-of-freedom of the system are established, and the stiffness coefficients of cracked rotor shaft are derived based upon the compliance method and strain energy release rate method. Then, utilizing the harmonic balance method and Taylor expansion technique, the unstable widths of simple and combination instability regions (SIR and CIR) are solved approximately. For a practical rotor system, the approximate unstable widths are verified by the Floquet numerical analysis. The effects of crack depth and position upon the unstable widths are discussed, and the conditions for zero unstable points (ZUPs) are given: Besides the asymmetric angle should be π/2 (for SIR) or 0 (for CIR), the relationships between the inertia asymmetry and crack parameters (depth and position) are also presented analytically. These results would be useful for crack detection and instability control of the asymmetric rotor-bearing system.     

各向异性压电材料平面裂纹的耦合场分析

用Stroh方法分析了各向异性压电材料电导通型裂纹问题的耦合场。结果表明,裂纹面上的切向电场强度和法向电位移均为常数,在裂纹尖端有由弹性场的耦事作用产生的奇异电导通裂纹模型中的静电场对裂纹尖端扩展的能量释放率不作贡献。     

混凝土双K断裂参数计算的半解析有限元法

混凝土裂缝扩展的双$K$断裂准则，用于描述混凝土结构裂缝的起裂、稳定扩 展和失稳断裂. 其相应的双$K$断裂参数(起裂断裂韧度$K_{\rm IC}^{\rm ini} $和失 稳断裂韧度$K_{\rm IC}^{\rm un}$)一般通过简便的试验和基于虚拟裂缝扩展粘 聚力的解析方法确定. 利用平面扇形域哈 密顿体系的方程，通过分离变量法及共轭辛本征函数向量展开法，以解析的方法推导出基于混 凝土虚拟裂缝扩展线性粘聚力模型的平面裂缝解析元列式. 将该解析元与有限元相结合，构成 半解析的有限元法，可求解任意结构几何形状的混凝土平面裂缝双$K$断裂参数的计算问题. 数值计算结果表明半解析有限元法对该类问题的求解是十分有效的.     

Evolving crack patterns in thin films with the extended finite element method

This paper develops the extended finite element method (XFEM) to evolve patterns of multiple cracks, in a brittle thin film bonded to an elastic substrate, with a relatively coarse mesh, and without remeshing during evolution. A shear lag model describes the deformation in three dimensions with approximate field equations in two-dimensions. The film is susceptible to subcritical cracking, obeying a kinetic law that relates the velocity of each crack to its energy release rate. At a given time, the XFEM solves the field equations and calculates the energy release rate of every crack. For a small time step, each crack is extended in the direction of maximal hoop stress, and by a length set by the kinetic law. To confirm the accuracy of the XFEM, we compare our simulation to the exiting solutions for several simple crack patterns, such as a single crack and a set of parallel cracks. We then simulate the evolution of multiple cracks, initially in a small region of the film but of different lengths, showing curved crack propagation and crack tip shielding. Starting with multiple small cracks throughout the film, the XFEM can generate the well-known mud crack pattern.     

压电压磁复合材料中正六边形孔边裂纹的反平面断裂问题

基于线性磁电弹性理论,利用Schwarz-Christoffel(CS)变换技术和Stroth公式,首次系统研究了压电压磁复合材料中含带两个不对称裂纹的正六边形孔口问题在部分渗透磁电边界条件下的解析解.当忽略磁场时,磁电非渗透裂纹和磁电渗透裂纹两种极端情况下的解析解答可退化为文献已有研究结果.数值结果揭示了正六边形孔口尺寸、裂纹长度以及力电载荷和磁载荷对能量释放率的影响规律.研究结果表明：减小孔口边长和裂纹长度可以提高材料的可靠性;机械载荷总是促进裂纹扩展;在磁电非渗透和磁电部分渗透边界条件下,负电场和负磁场会延缓裂纹的扩展,而正电场可以增强或阻碍裂纹的扩展,这取决于所施加的电场和磁场的强度以及机械载荷的水平;在磁电渗透边界条件下,电场和磁场对裂纹的扩展没有影响.     

冲击载荷下两裂纹间的连通规律

脆性材料内部含有大量裂纹，当某一裂纹扩展时，其他裂纹会对扩展裂纹产生影响。为了研究冲击载荷下，脆性材料内两裂纹的相互影响、连通规律及裂纹尖端应力强度因子的变化规律，利用有机玻璃板制作了含非平行双裂纹的实验试件，利用落板冲击设备进行了中低速冲击实验，结合有限元分析软件ABAQUS计算出裂纹尖端应力强度因子，利用有限差分软件AUTODYN进行了动态数值模拟研究，并将其模拟结果与实验结果进行对比分析。实验及模拟结果表明:裂纹破坏形态与AUTODYN数值模拟破坏形态基本一致；试件的断裂形态随着两裂纹间距不同而不同；裂纹间的相互影响程度随着裂纹间间距增大而减小；裂纹尖端应力强度因子K_I随着裂纹间距的增大而减小，而K_II随着裂纹间距增大而增大。