Driving force acting on a running crack

The asymptotic analysis of the dynamic crack problem for the anti‐plane shear mode is provided. The field near the crack tip is studied in detail for a nonlinear elastic incompressible material whose stored energy behaves asymptotically as a power of the first invariant of the strain tensor at large strains. It is shown that the hardening parameter characterizes fully the singularity degree of the near‐crack‐tip field. Based on the latter knowledge the driving force acting on the crack tip is calculated. Possible scenarios of the crack propagation are discussed.     

A generalized cohesive element technique for arbitrary crack motion

A computational method for arbitrary crack motion through a finite element mesh, termed as the generalized cohesive element technique, is presented. In this method, an element with an internal discontinuity is replaced by two superimposed elements with a combination of original and imaginary nodes. Conventional cohesive zone modeling, limited to crack propagation along the edges of the elements, is extended to incorporate the intra-element mixed-mode crack propagation. Proposed numerical technique has been shown to be quite accurate, robust and mesh insensitive provided the cohesive zone ahead of the crack tip is resolved adequately. A series of numerical examples is presented to demonstrate the validity and applicability of the proposed method.     

Crack propagation calculation using trial cracks

We present a method to perform finite element calculations for crack propagation problems with arbitrary crack directions in two dimensions. The crack direction (angle of propagation) is determined by inserting small “trial cracks” at the crack tip. For each trial crack, the domain is remeshed to allow crack propagation between elements. The trial cracks are then opened and the energy release rate is measured. The optimum crack direction (i.e., the crack direction with maximum energy release) is determined by an optimisation procedure. Although the method is computationally expensive due to the need to perform several calculations for each crack increment, it has the advantage that the energy release rate can be calculated even in cases where other methods fail. After explaining the method, it is applied to some test examples. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

弹性-幂硬化蠕变性材料Ⅱ型界面裂纹准静态扩展的渐近分析

建立了弹性-幂硬化蠕变性材料Ⅱ型界面裂纹准静态扩展的力学模型,求得了在裂纹表面自由和裂纹面有摩擦接触两种情况下,裂纹尖端应力场分离变量形式的渐近解．求解结果表明:Ⅱ型界面裂纹问题的应力、应变具有相同的奇异性;Ⅱ型界面裂纹尖端场不存在振荡奇异性;材料的幂硬化指数n和弹性模量比对裂纹尖端应力场幂硬化蠕变性材料区有着显著的影响,而弹性区仅受幂硬化指数n的影响,当n很大时,蠕变变形占主导地位,应力场趋于稳定,不随n的变化而变化;泊松比对裂纹尖端应力场的影响不明显．     

平面应变和反平面应变复合型裂纹尖端的各向异性塑性应力场

在裂纹尖端的理想塑性应力分量都只是θ的函数的条件下,利用平衡方程,各向异性塑性应力应变率关系、相容方程和Hill各向异性屈服条件,本文导出了平面应变和反平面应变复合型裂纹尖端的各向异性塑性应力场的一般解析表达式.将这些一般解析表达式用于复合型裂纹,我们就可以得到Ⅰ-Ⅲ、Ⅱ-Ⅲ及Ⅰ-Ⅱ-Ⅲ复合型裂纹尖端的各向异性塑性应力场的解析表达式.     

A Yoffe-type moving crack in one-dimensional hexagonal piezoelectric quasicrystals

A Yoffe-type moving crack in one-dimensional hexagonal piezoelectric quasicrystals is considered. The Fourier transform technique is used to solve a moving crack problem under the action of antiplane shear and inplane electric field. Full elastic stresses of phonon and phason fields and electric fields are derived for a crack running with constant speed in the periodic plane. Obtained results show that the coupled elastic fields inside piezoelectric quasicrystals depend on the speed of crack propagation, and exhibit the usual square-root singularity at the moving crack tip. Electric field and phason stresses do not have singularity and electric displacement and phonon stresses have the inverse square-root singularity at the crack tip for a permeable crack. The field intensity factors and energy release rates are obtained in closed form. The crack velocity does not affect the field intensity factors, but alters the dynamic energy release rate. Bifurcation angle of a moving crack in a 1D hexagonal piezoelectric quasicrystal is evaluated from the viewpoint of energy balance. Obtained results are helpful to better understanding crack advance in piezoelectric quasicrystals.     

线性硬化材料中稳恒扩展裂纹尖端场的粘塑性解

采用弹粘塑性力学模型,对线性硬化材料中平面应变扩展裂纹尖端场进行了渐近分析．假设人工粘性系数与等效塑性应变率的幂次成反比,通过量级匹配表明应力和应变均具有幂奇异性,奇异性指数由粘性系数中等效塑性应变率的幂指数唯一确定．通过数值计算讨论了Ⅱ型动态扩展裂纹尖端场的分区构造随各材料参数的变化规律．结果表明裂尖场构造由硬化系数所控制而与粘性系数基本无关．弱硬化材料的二次塑性区可以忽略,而较强硬化材料的二次塑性区和二次弹性区对裂尖场均有重要影响．当裂纹扩展速度趋于零时,动态解趋于相应的准静态解;当硬化系数为零时便退化为HR(Hui-Riedel)解．     

Evaluation of the stress singularities of plane V-notches in bonded dissimilar materials

According to the linear theory of elasticity, there exists a combination of different orders of stress singularity at a V-notch tip of bonded dissimilar materials. The singularity reflects a strong stress concentration near the sharp V-notches. In this paper, a new way is proposed in order to determine the orders of singularity for two-dimensional V-notch problems. Firstly, on the basis of an asymptotic stress field in terms of radial coordinates at the V-notch tip, the governing equations of the elastic theory are transformed into an eigenvalue problem of ordinary differential equations (ODEs) with respect to the circumferential coordinate θ around the notch tip. Then the interpolating matrix method established by the first author is further developed to solve the general eigenvalue problem. Hence, the singularity orders of the V-notch problem are determined through solving the corresponding ODEs by means of the interpolating matrix method. Meanwhile, the associated eigenvectors of the displacement and stress fields near the V-notches are also obtained. These functions are essential in calculating the amplitude of the stress field described as generalized stress intensity factors of the V-notches. The present method is also available to deal with the plane V-notch problems in bonded orthotropic multi-material. Finally, numerical examples are presented to illustrate the accuracy and the effectiveness of the method.     

关于平面应变和反平面应变复合型裂纹尖端的理想塑性应力场

文献[1]在裂纹尖端的理想塑性应力分量都只是θ的函数的条件下,利用平衡方程、应力应变率关系、相容方程和屈服条件导出了平面应变和反平面应变复合型裂纹尖端的理想塑性应力场的一般解析表达式。但文献[1]对应力应变率关系式中的比例因子λ(r,θ)作了很多限制,即假定λ与θ无关,并假定λ=c或cr-1。本文取消了对λ的这些限制。而文献[1]所研究的λ=crn(n=0或-1)的情形,只是本文的一个特殊情况。     

圆形界面刚性线夹杂的反平面问题

研究了在反平面集中力和无穷远纵向剪切作用下,不同弹性材料圆形界面上有多条刚性线夹杂的问题．运用Riemann-Schwarz解析延拓技术与复势函数奇性主部分析方法,首次获得了该问题的一般解答,求出了几种典型情况的封闭解,并给出了刚性线夹杂尖端的应力场分布A·D2结果表明,在反平面加载的情况下圆形界面刚性线夹杂尖端应力具有平方根奇异性,无奇异性应力振荡;应力场与刚性线夹杂的形状,加载方式和材料性质有关．退化结果与已有的解答完全吻合．     

弹性-粘塑性材料Ⅰ型动态扩展裂纹尖端场的渐近解

提出了一种新的弹性-粘塑性模型用于分析Ⅰ型动态扩展裂纹尖端的应力应变场.给出了适当的位移模式,推导了渐近方程并且给出了数值解.分析和计算表明:对于低粘性情况,裂纹尖端场具有对数奇异性;对于高粘性情况,渐近方程无解.分析比较表明该结果具有高玉臣提出的单参数解的所有优点,并且消除了粘性区随裂纹扩展而移动的不足.     

正交异性双材料的Ⅱ型界面裂纹尖端场

通过引入含16个待定实系数和两个实应力奇异指数的应力函数,再借助边界条件,得到了两个八元非齐次线性方程组.求解该方程组,在双材料工程参数满足适当条件下,确定了两个实应力奇异指数.根据极限唯一性定理,求出了全部系数,得到了应力函数的表示式.代入相应的力学公式,推出了当特征方程组两个判别式都小于0时,每种材料的裂纹尖端应力强度因子、应力场和位移场的理论解.裂纹尖端附近的应力和位移有混合型断裂特征,但没有振荡奇异性和裂纹面相互嵌入现象作为特例,当两种正交异性材料相同时,可以推出正交异性单材料Ⅱ型断裂的应力奇异指数、应力强度因子公式、应力场、位移场表示式.     

Ⅰ型动态扩展裂纹尖端的弹粘-理想塑性场

由于材料在扩展裂纹尖端的粘性效应的存在,考虑粘性效应并假设粘性系数与塑性等效应变率的幂次成反比,对理想塑性材料中平面应变扩展裂纹尖端场进行了弹粘塑性渐近分析,得到了不含间断的连续解,并讨论了Ⅰ型裂纹数值解的性质随各参数的变化规律．分析表明,应力和应变均具有幂奇异性,通过分析使尖端场的弹、粘、塑性可以合理匹配．对于Ⅰ型裂纹,裂尖场不含弹性卸载区．趋于极限情况时,裂纹尖端处于一种超粘性状态,并积聚了大量的能量,在各个受压应力状态下裂纹扩展．     

裂纹尖端场弹塑性分析的加权残数法及塑性应力强度因子的计算

本文以幂强化材料,平面应变情形为例,系统地提出了裂纹尖端场弹塑性分析的加权残数法,并根据此法,得出了裂纹尖端场的解析式弹塑性近似解.在此基础上.对整个裂纹区域,构造了弹塑性解叠加非线性有限元计算塑性应力强度因子的方法,从而为裂纹尖端场和整个裂纹体的分析和计算,提供了一个方法.     

Modeling crack micro-branching using finite elements with embedded strong discontinuities

The emphasis of this work lies in the development of a numerical method which is capable of representing the complex physical phenomena arising in the case of crack branching in brittle materials. In particular, the formation of crack micro-branches needs to be accounted for when it comes to the prediction of the propagation pattern of crack macro-branches which will ultimately lead to the failure of the material. This is achieved by numerically modeling the failure zones within the individual finite elements based on the concept of the embedded finite element method, where all the information with regard to the geometry of the failure zone is stored locally on the element level leading to a very efficient methodology capable of discretely resolving the failure zone. The main feature of the current work is the redundancy of the branching criterion based on crack tip velocity and that both, micro- as well as macro-branches can be modeled. Whether a micro-crack develops into a macro-crack solely depends on the local state of the material as it is outlined based on the application of the proposed numerical scheme on a rectangular block with a pre-existing notch set under tension. A comparison of the oscillatory behavior of the obtained crack tip velocity every time a micro-crack develops with experimental results from the literature is provided. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Asymptotics of eigenvalues of the nonlinear eigenvalue problem arising from the near mixed-mode crack-tip stress-strain field problems

In the present paper, approximate analytical and numerical solutions to nonlinear eigenvalue problems arising in nonlinear fracture mechanics in studying stress-strain fields near a crack tip under mixed-mode loading are presented. Asymptotic solutions are obtained by the perturbation method (the artificial small parameter method). The artificial small parameter is the difference between the eigenvalue corresponding to the nonlinear eigenvalue problem and the eigenvalue related to the linear “undisturbed” problem. It is shown that the perturbation technique is an effective method of solving nonlinear eigenvalue problems in nonlinear fracture mechanics. A comparison of numerical and asymptotic results for different values of the mixity parameter and hardening exponent shows good agreement. Thus, the perturbation theory technique for studying nonlinear eigenvalue problems is offered and applied to eigenvalue problems arising in fracture mechanics analysis in the case of mixed-mode loading.     

两压电介质之间的界面夹杂问题

应用Stroh理论,研究了两压电介质之间的刚性介电线夹杂问题。首先该问题被化为Hilbert问题,然后分别给出了压电介质内的复势函数解、夹杂内的电场解和夹杂尖端场的解析表达式。结果表明,在夹杂尖端附近,所有的场变量均呈现奇异性和振荡性,且其强度取决于介质的材料常数和无限场远处的应变。此外,结果还表明,当从夹杂内部趋近夹杂尖端时,夹杂内的电场也呈现奇异性和振荡性。     

An experimental-analytical hybrid model for the elasto-plastic stresses at a crack

The primary obstacle preventing the analytical determination of physically sensible stresses at a crack tip is the presence of a mathematical singularity there. This singularity is best known in its elastic form; however it persists even in elasto-plastic crack-tip stresses. To overcome the difficulty we adopt the following strategy: we attempt to capture initial elastic stresses experimentally, than track subsequent elasto-plastic stress distributions analytically.We infer a finite stress at a crack tip from the experimental behaviour of cracked specimens at fracture when the specimens are made of a truly brittle material. Given a size-independent result, we argue that the crack-tip stress at fracture must equal the ultimate stress for such a material; thus dividing by the applied stress at the same point gives a measure of the stress concentration factor, K_T. The approach is checked for size independence and against hole configurations with known theoretical, yet physically reasonable, K_T. Then the effective experimental KT are taken as inputs for the second phase of the study in which we model the crack as being a smooth notch having the same stress concentration factor as found experimentally. In this way our configuration initially shares the same stresses at the crack tip as we inferred physically. Next we track effects of incremental plastic flow on a set of finite element grids. Satisfactory resolution in return for modest computational effort is obtained by employing a substructuring method. The accuracy in both the elastic and the elasto-plastic regime is checked against trial problems with exact solutions. Thereafter, physically interpretable stress distributions ahead of the crack are determined for a range of materials and for varying load levels.     

压电陶瓷板中非电渗透型反平面裂纹的电弹性场

对受4种机电载荷的内含裂纹的压电陶瓷板的电弹性行为进行了分析。利用积分变换方法将非电渗透型反平面裂纹问题化为对偶积分方程组,求解这些方程组可以获得裂纹线上电弹性场的明显解析表达式,及裂尖处一些量的强度因子和机械应变能释放率。当板的厚度趋近于无穷大时,所得结果还原为熟知结果。     

On corner singularities in Reissner-Mindlin's theory of elastic plates

In the framework of linear elasticity, singularities occur in domains with non-smooth boundaries. Particularly in Fracture Mechanics, the local stress field near stress concentrations is of interest. In this work, singularities at re-entrant corners or sharp notches in Reissner-Mindlin plates are studied. Therefore, an asymptotic solution of the governing system of partial differential equations is obtained by using a complex potential approach which allows for an efficient calculation of the singularity exponent λ. The effect of the notch opening angle and the boundary conditions on the singularity exponent is discussed. The results show, that it can be distinguished between singularities for symmetric and antisymmetric loading and between singularities of the bending moments and the transverse shear forces. Also, stronger singularities than the classical crack tip singularity with free crack faces are observed. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)