CTS试件中复合型疲劳裂纹扩展

针对复合型循环载荷作用下的金属构件中的裂纹扩展问题进行了实验分析和理论建模. 首先 采用紧凑拉剪试件(CTS)和 Richard研制的复合型载荷加载装置，对承受复合型循环载荷的裂纹进行了实验研究. 实验选择了两种金属材料试件，分别承受3种形式的复合型循环载荷的作用，在裂纹尖端具 有相同的初始应力场强度的条件下考察复合型循环载荷对裂纹扩展规律的影响. 实验结果表明，疲劳裂纹的扩展速率与加载角度有关. 对于同样金属材料的试件，当裂尖处 初始应力场强度相等时，载荷越接近于II型，裂纹增长速率越快. 采用等效应力强度 因子(I型和II型应力强度因子的组合)、裂纹扩展速率及复合强度等参数，以实验数据为 基础，建立了一个疲劳裂纹扩展模型，用来预测裂纹在不同模式疲劳载荷作用下的扩展速率. 为验证其有效性，该模型被应用于钢制试件的数值模拟计算中. 实验结果与模拟计算曲线保 持一致，表明该模型可以用来估算带裂纹金属构件的寿命.     

Electroelastic analysis of a penny-shaped crack in a piezoelectric ceramic under mode I loading

Following the theory of linear piezoelectricity, we consider the electroelastic problem for a piezoelectric ceramic with a penny-shaped crack under mode I loading. The problem is formulated by means of Hankel transform and the solution is solved exactly. The stress intensity factor, energy release rate and energy density factor for the exact and impermeable crack models are expressed in closed form and compared for a P-7 piezoelectric ceramic. Based on current findings, we suggest that the energy release rate and energy density factor criteria for the exact crack model are superior to fracture criteria for the impermeable crack model.     

压电材料中椭圆片状裂纹问题精确解的一种方法

在线性压电陶瓷本构关系和裂纹边界绝缘的框架下,用超奇异积分方程的方法对椭圆类片状裂纹问题进行了重新研究．超奇异积分方程中的未知位移间断和电势间断近似地表示为基本密度函数与多项式之积,其中基本密度函数反映了椭圆片状裂纹前沿电弹性场的奇异性,而多项式在均布载荷作用下可用一个常数来表达．引入椭球坐标系后,得到了均布载荷作用下未知位移间断和电势间断的解析解．使用这些解析解和电弹性场强度的定义,得到了裂纹前沿Ⅰ型、Ⅱ型和Ⅲ型应力强度因子以及电位移强度因子的精确表达式．法向均布载荷作用下的结果与现有精确解完全一致,切向均布载荷作用下的结果则尚未见有其它报道．     

Multiple parallel cracks interaction problem inpiezoelectric ceramics

This paper has two goals. First, we propose the pseudo-traction–electric displacementmethod for solving the interaction problem of multiple parallel cracks in transversely isotropicpiezoelectric ceramics. Second, we present a fundamental understanding for the role that theelectric displacement loading plays in the interaction problem. Detailed comparisons between theresults under the compound mechanical–electric loading conditions and those derived underpurely mechanical loading conditions are performed. It is shown that the mechanical fractureparameters such as the stress intensity factors are no longer independent of the electric loading asthey would be in single crack problems. Quite contrary, the electric displacement loading has asignificant influence on the stress intensity factors, the total potential energy release rate and themechanical strain energy release rate. This important conclusion is mainly due to the interactioneffect, i.e., one of the multiple cracks releases the stresses and disturbs the electric fields near theother crack. It is also found that there are some special relative locations for the multiple parallelcracks at which the electric displacement loading has no effect on the Mode I stress intensityfactor. However, the mechanical strain energy release rate has no such a property.     

Crack propagation in an elastic solid subjected to general loading— IV. Obliquely incident stress pulse

The stress intensity factors of a half-plane crack extending nonuniformly in an isotropic elastic solid subjected to stress wave loading are considered. A plane stress pulse is obliquely incident on the crack, and the wavefront strikes the crack at some initial time. At some arbitrary later time, the crack begins to extend at a nonuniform rate. It is found that the mode I and mode II stress intensity factors each have the form of the product of a universal function of instantaneous cracktip speed with the stress intensity factor for an equivalent stationary crack. An energy-rate balance fracture criterion is applied to obtain an equation of motion for the crack tip and to determine the actual delay time between the arrival of the incident wave and the onset of fracture as a function of angle of incidence of the loading wave.     

A computational study of local stress intensity factor solutions for kinked cracks near spot welds in lap-shear specimens

In this paper, the local stress intensity factor solutions for kinked cracks near spot welds in lap-shear specimens are investigated by finite element analyses. Based on the experimental observations of kinked crack growth mechanisms in lap-shear specimens under cyclic loading conditions, three-dimensional and two-dimensional plane-strain finite element models are established to investigate the local stress intensity factor solutions for kinked cracks emanating from the main crack. Semi-elliptical cracks with various kink depths are assumed in the three-dimensional finite element analysis. The local stress intensity factor solutions at the critical locations or at the maximum depths of the kinked cracks are obtained. The computational local stress intensity factor solutions at the critical locations of the kinked cracks of finite depths are expressed in terms of those for vanishing kink depth based on the global stress intensity factor solutions and the analytical kinked crack solutions for vanishing kink depth. The three-dimensional finite element computational results show that the critical local mode I stress intensity factor solution increases and then decreases as the kink depth increases. When the kink depth approaches to 0, the critical local mode I stress intensity factor solution appears to approach to that for vanishing kink depth based on the global stress intensity factor solutions and the analytical kinked crack solutions for vanishing kink depth. The two-dimensional plane-strain computational results indicate that the critical local mode I stress intensity factor solution increases monotonically and increases substantially more than that based on the three-dimensional computational results as the kink depth increases. The local stress intensity factor solutions of the kinked cracks of finite depths are also presented in terms of those for vanishing kink depth based on the global stress intensity factor solutions and the analytical kinked crack solutions for vanishing kink depth. Finally, the implications of the local stress intensity factor solutions for kinked cracks on fatigue life prediction are discussed.     

Analytical representation of the non-square-root singular stress field at a finite angle sharp notch

The stress field near the tip of a finite angle sharp notch is singular. However, unlike a crack, the order of the singularity at the notch tip is less than one-half. Under tensile loading, such a singularity is characterized by a generalized stress intensity factor which is analogous to the mode I stress intensity factor used in fracture mechanics, but which has order less than one-half. By using a cohesive zone model for a notional crack emanating from the notch tip, we relate the critical value of the generalized stress intensity factor to the fracture toughness. The results show that this relation depends not only on the notch angle, but also on the maximum stress of the cohesive zone model. As expected the dependence on that maximum stress vanishes as the notch angle approaches zero. The results of this analysis compare very well with a numerical (finite element) analysis in the literature. For mixed-mode loading the limits of applicability of using a mode I failure criterion are explored.     

In-plane perturbation of a system of two coplanar slit-cracks – I: Case of arbitrarily spaced crack fronts

In order to lay the grounds for a future study of the deformation of the fronts of coplanar cracks during their final coalescence, we consider the model problem of a system of two coplanar, parallel, identical slit-cracks loaded in mode I in some infinite body. The first, necessary task is to determine the distribution of the stress intensity factors along the crack fronts resulting from some small but otherwise arbitrary in-plane perturbation of these fronts. This is done here in the case where the distances between the various crack fronts are arbitrary and fixed.The first order expression of the local variation of the stress intensity factor is provided by a general formula of Rice (1989) in terms of some “fundamental kernel” tied to the mode I crack face weight function. In the specific case considered, this fundamental kernel reduces to six unknown functions; the problem is to determine them. This is done by using another formula of Rice (1989) which provides the variation of the fundamental kernel in a similar way. This second formula is applied to special perturbations of the crack fronts preserving the shape and relative dimensions of the cracks while modifying their absolute size and orientation. The output of this procedure consists of nonlinear integro-differential equations on the functions looked for, which are transformed into nonlinear ordinary differential equations through Fourier transform in the direction of the crack fronts, and then solved numerically.     

Stable crack growth along a brittle\ductile interface—II. Small scale yielding solutions and interfacial toughness predictions

The problem of a crack growing steadily and quasi-statically along a brittle\ductile interface under plane strain, mixed mode, and small scale yielding conditions is considered. The ductile material is assumed to be characterized by the J₂-flow theory of plasticity with linear strain hardening, while the brittle material is assumed to be linear elastic. A displacement-based finite element method, exploiting the convective nature of the problem, is utilized to solve the relevant boundary value problem. In Part I of this work, the corresponding asymptotic problem was solved. This paper addresses the full-field problem in order to validate the asymptotic solutions, and to explore the physical implications of the results. The numerical full-field results are found to be in good agreement with the analytical asymptotic solutions. In particular, the full-field results strongly suggest that the stress fields in the vicinity of the crack tip are variable-separable of the power singular type; and also that the mode mix of the near-tip stress fields is, to a large extent, independent of the applied elastic mode mix. The amplitude (the plastic stress intensity factor) and the regions of validity of the asymptotic fields are estimated from the full-field results, and are observed to be strongly dependent on the applied mode mix. The remote elastic loading fields appear to influence the near-tip fields, primarily, through the plastic stress intensity factor. The present work also explores the suggestion made by Bose and Ponte Castaneda, 1992 that the solutions to the small scale yielding problem may be used in the context of a standard crack growth criterion, requiring that continued growth take place with a fixed near-tip crack opening profile, to obtain theoretical predictions for the dependence of interfacial toughness on the applied mode mix. Based on the numerical results, predictions for mixed mode toughness of the brittle\ductile interface are reported. The results, which are in qualitative agreement with available experimental data and also with some recent theoretical results, predict a strong dependence of interfacial toughness on mode mix. This suggests that ductility provides the main operating mechanism for explaining the dependence of interfacial toughness on the mode mix of the applied loading fields, during steady crack growth.     

NiTi形状记忆合金裂纹尖端热力耦合行为研究

本文对NiTi形状记忆合金I型裂纹尖端热力耦合行为进行了数值仿真分析和实验验证。建立了包含相变和热力耦合的本构模型,通过有限元计算得到了裂纹尖端附近的纵向应变、马氏体体积分数和温度场分布,依据马氏体相变情况对裂纹尖端有效应力强度因子进行了修正,揭示了加载速率对形状记忆合金裂纹尖端有效应力强度影子的影响规律。参数研究表明,随着加载频率的增加,裂纹尖端附近温度逐渐升高,马氏体相变区域逐渐缩小,有效应力强度因子呈下降趋势,形状记忆合金表现出增韧效应,有助于减缓裂纹扩展。本研究结果对于揭示热力耦合作用下超弹性形状记忆合金疲劳裂纹扩展规律具有重要参考意义。     

Mode I fracture analysis of a piezoelectric strip based on real fundamental solutions

In existing papers, mode I crack problems of piezoelectric ceramics are generally solved in complex domain because of the complex fundamental solutions of in-plane piezoelectric governing equations. In fact, these problems can alternatively be analyzed in real number field by recasting the solutions in real form instead. The main purpose of the present work is to develop such real fundamental solutions by detailed eigenvalue and eigenvector analyses. As an example of application, the widely studied fracture problem of a piezoelectric strip with a center-situated crack under mode I loading condition is then revisited based on the real fundamental solutions. Mixed boundary value conditions of the crack are transformed into Cauchy singular integral equations, which are then solved numerically to get fracture parameters including the energy release rate and intensity factors. Convergence behaviors of the kernel functions are surveyed. Theoretical derivation and computation are validated by the exact solution in a special case. The effect of a combined geometrical parameter on the crack is also investigated.     

Development of fracture facets from a crack loaded in mode I+III: Solution and application of a model 2D problem

It is experimentally well-known that a crack loaded in mode I+III propagates through formation of discrete fracture facets inclined at a certain tilt angle on the original crack plane, depending on the ratio of the mode III to mode I initial stress intensity factors. Pollard et al. (1982) have proposed to calculate this angle by considering the tractions on all possible future infinitesimal facets and assuming shear tractions to be zero on that which will actually develop. In this paper we consider the opposite case of well-developed facets; the stress field near the lateral fronts of such facets becomes independent of the initial crack and essentially 2D in a plane perpendicular to the main direction of crack propagation.To determine this stress field, we solve the model 2D problem of an infinite plate containing an infinite periodic array of cracks inclined at some angle on a straight line, and loaded through uniform stresses at infinity. This is done first analytically, for small values of this angle, by combining Muskhelishvili's (1953) formalism and a first-order perturbation procedure. The formulae found for the 2D stress intensity factors are then extended in an approximate way to larger angles by using another reference solution, and finally assessed through comparison with some finite element results.To finally illustrate the possible future application of these formulae to the prediction of the stationary tilt angle, we introduce the tentative assumption that the 2D mode II stress intensity factor is zero on the lateral fronts of the facets. An approximate formula providing the tilt angle as a function of the ratio of the mode III to mode I stress intensity factors of the initial crack is deduced from there. This formula, which slightly depends on the type of loading imposed, predicts somewhat smaller angles than that of Pollard et al. (1982).     

Dynamic crack propagation in a viscoelastic strip

The dynamic propagation of a semi-infinite crack in a finite linear viscoelastic strip subjected to Mode I loading is investigated. Through the use of integral transforms the problem is reduced to solving a Wiener-Hopf equation. The asymptotic properties of the transforms are exploited to establish the stress intensity factor. Plane-stress and plane-strain stress intensity factors as a function of crack speed for both fully-clamped and shear-free lateral boundaries are presented for the standard linear viscoelastic solid. Comparisons are made with previously obtained asymptotic stress intensity factors and with stress intensity factors for the equivalent elastic strips.     

层状磁电复合材料界面共线裂纹平面问题分析

本文研究了面内电磁势载荷作用下双层压电压磁复合材料中共线界面裂纹问题．考虑了压电材料的导磁性质和压磁材料的介电性质,引入了界面电位移和磁感强度的连续性条件．利用Fourier 变换得到一组第二类Cauchy 型奇异积分方程．进一步导出了相应问题的应力强度因子、电位移强度因子和磁感强度强度因子的表达式,给出了应力强度因子的数值结果．结果表明电磁载荷会导致界面裂纹尖端I、II 混合型应力奇异性,同时还伴随着电位移和磁感强度的奇异性．比较了双裂纹左右端的应力强度因子,发现在面内极化方向上施加面内磁势载荷时共线裂纹内侧尖端区域的两个法向应力场发生互相干涉增强．     

裂纹在冲击载荷作用下起裂的临界载荷面

采用有限元方法研究裂纹在I型短脉冲载荷作用下应力强度因子随时间的变化 ,用应力强度因子的初始上升时间Tr 对时间坐标无量纲化 ,对应力强度因子初始上升段进行曲线拟合 ,得到了上升段的曲线表达式。运用简单弹性梁理论和Lagrangian运动方程 ,获得载荷与时间对裂纹作用的关系式 ,结合有限元的结果 ,得到了上升时间Tr 的计算表达式 ,并进一步推出了裂纹在冲击载荷作用下起裂的临界载荷面。     

Contact Zone Approach for a Moving Interface Crack in a Piezoelectric Bimaterial under Thermoelectromechanical Loading

An interface crack of a finite length moving with a constant subsonic speed v along an interface of two semi-infinite piezoelectric spaces is considered. It is assumed that the bimaterial compound is loaded by a remote mixed mode mechanical loading and a thermoelectrical field and that a frictionless contact zone arises at the leading crack tip. Electrically permeable and electrically insulated cases of the open part of the crack are involved into the consideration. By introducing a moving coordinate system at the crack tip the problem is reduced to a combined Dirichlet–Riemann boundary value problem which is solved exactly. For both cases of the electrical conditions the transcendental equations are obtained for the determination of the real contact zone length, and moreover, the associated closed form asymptotic formulas are found for small values of this parameter. Variations of the contact zone length and the stress intensity factor with respect to the crack speed and the loading have been investigated both for electrically permeable and electrically insulated cases.     

Perturbation of a dynamic planar crack moving in a model viscoelastic solid

Weight functions, which give stress intensity factors in terms of applied loading, are constructed, for three-dimensional time-dependent loading of a semi-infinite crack, propagating at uniform speed. Both a model problem, governed by a scalar wave equation, and the full vectorial problem for Mode I loading, are considered. The medium through which the crack propagates is viscoelastic; the approach is general but explicit formulae are given when the medium is a Maxwell fluid. The weight functions are exploited to develop formulae for the first-order perturbations of stress intensity factors when the crack edge is no longer straight but becomes slightly wavy. Implications for stability, and for “crack front waves” in the case of the Mode I problem, are discussed.     

Analytic solutions to crack tip plastic zone under various loading conditions

Based on stress field equations and Hill yield criterion, the crack tip plastic zone is determined for orthotropic materials and isotropic materials under small-scale yielding condition. An analytical solution to calculating the crack tip plastic zone in plane stress states is presented. The shape and size of the plastic zone are analyzed under different loading conditions. The obtained results show that the crack tip plastic zones present “butterfly-like” shapes, and the elastic–plastic boundary is smooth. The size of the plastic zone for orthotropic composites is less at the crack tip for various loading conditions, compared with the case of isotropic materials. Crack inclination angle and loading conditions affect greatly the size and shape of crack tip plastic zone. The mode I crack has a crucial effect on the plastic zone for mixed mode case in plane stress state. The plastic zone for pure mode I crack and pure mode II crack have a symmetrical distribution to the initial crack plane.     

Full field analysis of an anti-plane interface crack subjected to dynamic body forces

In this study, the transient full field response of an interface crack between two different media subjected to dynamic body force at one material is investigated. For time t < 0, the bimaterial medium is stress free and at rest. At t = 0, a concentrated anti-plane dynamic point loading is applied at the medium as shown in Fig. 1. The total wave field is due to the effect of this point loading and the scattering of the incident waves by the interface crack. An alternative methodology that is different from the conventional superposition method is used to construct the reflected, refracted and diffracted wave fields. A useful fundamental solution is proposed in this study and the full field solution is determined by superposition of the fundamental solution in the Laplace transform domain. The proposed fundamental problem is the problem of applying an exponentially distributed traction (in the Laplace transform domain) on the interfacial crack faces. The Cagniard–de Hoop method of Laplace inversion is used to obtain the transient solution in time domain. Exact transient closed form solutions for stresses and stress intensity factors are obtained. Numerical results for the time history of stresses and stress intensity factors during the transient process are discussed in detail.     

Loading effect on ACPD of a crack in ferromagnetic material

The loading effect on alternating current potential drop (ACPD) for a ferromagnetic material containing a two-dimensional surface crack was investigated under opening mode loading without shear (mode I). The change in potential drop due to load was obtained with and without a magnetic field around the specimen. To remove the magnetic field from the circumference of the specimen, a new measuring system using the characteristic of coaxial transmission line was made. When the magnetic field does not exist around the specimen, the change in potential drop with load is governed by the change in electromagnetic properties near the crack tip. The results obtained by using the new measuring system are the basis for an application of the ACPD technique to the experimental determination of the stress intensity factor, since they are independent of the arrangement of the measuring probe lines and the current supply lines. The relationship between the change in potential drop and the change in load is linearized by demagnetization. The change in potential drop per unit change in the stress intensity factor is independent of the crack length.