Interaction between heat dipole and circular interracial crack

The heat dipole consists of a heat source and a heat sink. The problem of an interfacial crack of a composite containing a circular inclusion under a heat dipole is investigated by using the analytical extension technique, the generalized Liouville theo-rem, and the Muskhelishvili boundary value theory. Temperature and stress fields are formulated. The effects of the temperature field and the inhomogeneity on the interracial fracture are analyzed. As a numerical illustration, the thermal stress intensity factors of the interfacial crack are presented for various material combinations and different po-sitions of the heat dipole. The characteristics of the interfacial crack depend on the elasticity, the thermal property of the composite, and the condition of the dipole.     

螺型位错偶极子与界面钝裂纹的干涉效应

研究了螺型位错偶极子与界面钝裂纹的干涉效应.应用保角变换技术,得到了复势函数与应力场的封闭解析解,讨论了位错偶极子方位、臂长及裂纹钝化程度对位错偶板子屏蔽效应和发射条件的影响.结果表明,与单个螺型位错不同,螺型位错偶极子与x轴夹角在一定范围内时才可以降低界面钝裂纹尖端的应力强度因子(屏蔽效应),屏蔽效应随偶板子臂长的增大而增强,随裂纹钝化程度的增大而增强,屏蔽区域也随裂纹钝化程度的增大而增大;位错偶极子发射所需的临界无穷远加载随偶极子臂长的增加而减小,随位错方位角及裂纹钝化程度的增加而增大;最可能的位错偶极子发射角度为0.螺型位错偶极子的发射比单个螺型位错的发射要困难.本文解答的特殊情况与相关文献给出的解答一致.     

Perturbation analysis of Mode III interfacial cracks advancing in a dilute heterogeneous material

The paper addresses the problem of a Mode III interfacial crack advancing quasi-statically in a heterogeneous composite material, that is a two-phase material containing elastic inclusions, both soft and stiff, and defects, such as microcracks, rigid line inclusions and voids. It is assumed that the bonding between dissimilar elastic materials is weak so that the interface is a preferential path for the crack. The perturbation analysis is made possible by means of the fundamental solutions (symmetric and skew-symmetric weight functions) derived in Piccolroaz et al. (2009). We derive the dipole matrices of the defects in question and use the corresponding dipole fields to evaluate “effective” tractions along the crack faces and interface to describe the interaction between the main interfacial crack and the defects. For a stable propagation of the crack, the perturbation of the stress intensity factor induced by the defects is then balanced by the elongation of the crack along the interface, thus giving an explicit asymptotic formula for the calculation of the crack advance. The method is general and applicable to interfacial cracks with general distributed loading on the crack faces, taking into account possible asymmetry in the boundary conditions.The analytical results are used to analyse the shielding and amplification effects of various types of defects in different configurations. Numerical computations based on the explicit analytical formulae allows for the analysis of crack propagation and arrest.     

Thermoelastic problem of steady-state heat flows disturbed by a crack at an arbitrary angle to the graded interfacial zone in bonded materials

Plane thermoelasticity solutions are presented for the problem of a crack in bonded materials with a graded interfacial zone. The interfacial zone is treated as a nonhomogeneous interlayer having spatially varying thermoelastic moduli between dissimilar, homogeneous half-planes. The crack is assumed to exist in one of the half-planes at an arbitrary angle to the graded interfacial zone, disturbing uniform steady-state heat flows. The Fourier integral transform method is employed in conjunction with the coordinate transformations of field variables in the basic thermoelasticity equations. Formulation of the current nonisothermal crack problem lends itself to the derivation of two sets of Cauchy-type singular integral equations for heat conduction and thermal stress analyses. The heat-flux intensity factors and the thermal-stress intensity factors are defined and evaluated in order to quantify the singular characters of temperature gradients and thermal stresses, respectively, in the near-tip region. Numerical results include the variations of such crack-tip field intensity factors versus the crack orientation angle for various combinations of material and geometric parameters of the dissimilar media bonded through the thermoelastically graded interfacial zone. The dependence of the near-tip thermoelastic singular field on the degree of crack-surface partial insulation is also addressed.     

Stress intensity factor for an edge crack in two bonded dissimilar materials under convective cooling

The transient thermal stress crack problem for two bonded dissimilar materials subjected to a convective cooling on the surface containing an edge crack perpendicular to the interface is considered. The problem is solved using the principle of superposition and the uncoupled quasi-static thermoelasticity. The crack problem is formulated by applying the transient thermal stresses obtained from the uncracked medium with opposite sign on the crack surfaces to be the only external loads. Fourier integral transform is used to solve the perturbation problem resulting in a singular integral equation of Cauchy type in which the derivative of the crack surface displacement is the unknown function. The numerical results of the stress intensity factors are calculated for both the edge crack and the crack terminating at the interface using two different composite materials and illustrated as a function of time, crack length, coefficient of heat transfer, and the thickness ratio.     

Stress intensity factor for a kinked interfacial crack in dissimilar anisotropic materials under antiplane deformation

The asymptotic problem of a kinked interfacial crack in dissimilar anisotropic materials under antiplane deformation is investigated. The linear transformation method for the problem of the anisotropic bimaterial with a straight interface is proposed. The stress intensity factor for the kinked interfacial crack in the anisotropic composite is obtained from the solution of the transformed problem of the kinked interfacial crack in the isotropic bimaterial based on the linear transformation method. The effects of the material parameters as well as the kink angle on the stress intensity factor are discussed from numerical results of the stress intensity factor. The finite element analysis is carried out to verify the stress intensity factor obtained by using the linear transformation. The influence of the material orientations on the stress intensity factor is investigated for the kinked crack in the bimaterial consisting of dissimilar inclined orthotropic materials.     

基于纯声子辐射模型的超薄绝缘夹层结构

在复合材料结构中起绝热、增韧作用的绝缘夹层,其加工厚度现在已达到纳米量级,原有的傅立叶热传导定律已无法描述其热能的传递行为,需从分子动力学、量子力学从发,针对不同研究对象建立相应的热传导模型．针对超薄绝缘夹层结构,将纯声子辐射模型和傅立叶热传导模型相结合数值求解热冲击条件下的温度场,并作为热载荷,用于求解结构上表面应力和夹层裂纹驱动力,其结果与只采用傅立叶热传导模型计算的结果相比较, 分析了物理参数对温度、应力和裂纹驱动力的影响．结果表明：与只采用傅立叶热传导模型计算的结果相比,按EPRT计算的热传导明显变慢,其表面剥离应力偏大,而夹层裂纹驱动力偏小．同时随着松弛时间增大和声子速度的降低,热传导减缓,表面横向剥离应力增大,超薄绝缘夹层内裂纹尖端驱动力减小．     

Theoretical model of crack branching in magnetoelectric thermoelastic materials

Thermomagnetoelectroelastic crack branching of magnetoelectro thermoelastic materials is theoretically investigated based on Stroh formalism and continuous distribution of dislocation approach. The crack face boundary condition is assumed to be fully thermally, electrically and magnetically impermeable. Explicit Green’s functions for the interaction of a crack and a thermomagnetoelectroelastic dislocation (i.e., a thermal dislocation, a mechanical dislocation, an electric dipole and a magnetic dipole located at a same point) are presented. The problem is reduced to two sets of coupled singular integral equations with the thermal dislocation and magnetoelectroelastic dislocation densities along the branched crack line as the unknown variables. As a result, the formulations for the stress, electric displacement and magnetic induction intensity factors and energy release rate at the branched crack tip are expressed in terms of the dislocation density functions and the branch angle. Numerical results are presented to study the effect of applied thermal flux, electric field and magnetic field on the crack propagation path by using the maximum energy release rate criterion.     

INTERACTION OF A DISLOCATION DIPOLE WITH A MODE Ⅲ INTERFACE CRACK

Interaction between a screw dislocation dipole and a mode Ⅲ interface crack is investigated. By using the complex variable method, the closed form solutions for complex potentials are obtained when a screw dislocation dipole lies inside a medium. The stress fields and the stress intensity factors at the tip of the interface crack produced by the screw dislocation dipole are given. The influence of the orientation, the dipole arm and the location of the screw dislocation dipole as well as the material mismatch on the stress intensity factors is discussed. zThe image force and the image torque acting on the screw dislocation dipole center are also calculated. The mechanical equilibrium position of the screw dislocation dipole is examined for various material property combinations and crack geometries. The results indicate that the shielding or anti-shielding effect on the stress intensity factor increases abruptly when the dislocation dipole approaches the tip of the crack. Additionally, the disturbation of the interface crack on the motion of the dislocation dipole is also significant.     

INTERACTION OF A DISLOCATION DIPOLE WITH A MODE III INTERFACE CRACK

Interaction between a screw dislocation dipole and a mode III interface crack is investigated. By using the complex variable method, the closed form solutions for complex potentials are obtained when a screw dislocation dipole lies inside a medium. The stress fields and the stress intensity factors at the tip of the interface crack produced by the screw dislocation dipole are given. The influence of the orientation, the dipole arm and the location of the screw dislocation dipole as well as the material mismatch on the stress intensity factors is discussed. The image force and the image torque acting on the screw dislocation dipole center are also calculated. The mechanical equilibrium position of the screw dislocation dipole is examined for various material property combinations and crack geometries. The results indicate that the shielding or anti-shielding effect on the stress intensity factor increases abruptly when the dislocation dipole approaches the tip of the crack. Additionally, the disturbation of the interface crack on the motion of the dislocation dipole is also significant.     

Axially symmetric thermal stress of an external circular crack under general thermal conditions

Summary The paper presents a solution for the linear thermoelastic problem of determining axisymmetric stress and displacement fields in an isotropic elastic solid of infinite extent weakened by an external circular crack under general mechanical loadings and general thermal conditions. The mechanical loadings and thermal conditions applied on the crack faces are axisymmetric, being non-symmetric about the crack plane. In similar lines of [7], equations of equilibrium of an elastic solid conducting heat have been solved using Hankel transforms and Abel operators of the first kind. Expressions for stress, displacement, temperature and heat flux functions are obtained in terms of Abel transforms of the first kind of the jumps of stress, displacement, temperature and heat flux at the crack plane. Two types of thermal conditions, that is, general surface temperatures and general heat flux on faces of the crack are considered. In both the cases, closed form solutions have been obtained for the unknown functions solving Abel type of integral equations. Explicit expressions for stresses, displacements, temperature fields, stress intensity factors have been obtained. Two special cases of thermal conditions in which: (i) crack faces are subjected to constant non-symmetric temperatures over a circular ring area, (ii) crack faces are subjected to constant non-symmetric heat flux over a circular ring area, have been considered. In some special cases, results have been compared with those from the literature.     

Impact of thermal loads on interfacial debonding in FRP strengthened beams

The effect of thermal loads on the debonding mechanisms in beams strengthened with externally bonded composite materials is analytically investigated. The analytical approach adopts a high-order stress analysis model and a fracture mechanics model that uses the concept of the energy release rate through the thermo-mechanical form of the J-integral. The two models are combined to synthesize the relation between the energy release rate, the mechanical loads, the thermal loads, and the interfacial crack length simulating the thermo-mechanical debonding process. The model is supported through comparison with experimental results taken from the literature. The comparison quantifies and explains various phenomena observed in the experiments and mainly the non-monotonic dependency of the debonding failure load on the temperature. The impact of the temperature on the interfacial stresses and on the stability of the debonding process is also studied. Finally, the effect of an uniform thermal load on the debonding behavior of a strengthened beam is studied revealing the impact of the thermal load on the debonding stability and strength characteristics.     

A rigid triangular inclusion partially bonded in an elastic matrix

The two-dimensional problem of a rigid rounded-off angle triangular inclusion partially bonded in an infinite elastic plate is studied. The unbonded part of the inclusion boundary forms an interfacial crack. Based on the complex variable method for curvilinear boundaries, the problem is reduced to a non-homogeneous Hilbert problem and the stress and displacement fields in the plate are obtained in closed form. Special attention is paid in the investigation of the stress field in the vicinity of the crack tip. It is found that the stresses present an oscillatory singularity and the general equations for the local stresses are derived. The singular stress field is coupled with the maximum circumferential stress and the minimum strain energy density criteria to study the fracture characteristics of the composite plate. Results are given for the complex stress intensity factors, the local stresses, the crack extension angles and the critical applied loads for unstable crack growth from its more vulnerable tip or two types of interfacial cracks along the inclusion boundary.     

位错偶极子和Ⅲ型界面裂纹的干涉效应

研究了多晶体材料中螺型位错偶极子和界面裂纹的弹性干涉作用.利用复变函数方法,得到了该问题复势函数的封闭形式解答.求出了由位错偶极子诱导的应力场和裂纹尖端应力强度应子,分析了偶极子的方向,偶臂和位置以及材料失配对应力强度因子的影响.推导了作用在螺型位错偶极子中心的像力和力偶矩,并讨论了界面裂纹几何条件和不同材料特征组合对位错偶极子平衡位置的影响规律.结果表明,裂纹尖端的螺型位错偶极子对应力强度因子会产生强烈的屏蔽或反屏蔽效应.同时,界面裂纹对螺型位错偶极子在材料中运动有很强的扰动作用.     

含埋藏椭圆形裂纹金属构件电磁热止裂时热应力场分析

对带有椭圆埋藏裂纹金属构件在脉冲放电瞬间的热应力场进行了理论分析。在热应力场的求解中,利用了热传导和非定常热应力理论,由求解的热应力场公式发现：在脉冲放电瞬间,电磁热在裂纹尖端形成热压应力场,热压力场可有效地抑制裂纹的扩展。以Cr12MoV模具钢中埋藏椭圆裂纹止裂为例,具体计算了脉冲放电瞬间的热应力场分布情况,为空间裂纹电磁热止裂技术的实际应用提供了理论基础。     

A technique for studying interaction between a screw dislocation dipole or a concentrated load and a mode III crack crossing an interface

A method of potentially wide application is developed for deriving analytical expressions of the elastic interaction between a screw dislocation dipole or a concentrated force and a crack cutting perpendicularly across the interface of a bimaterial. The cross line composed of the interface and the crack is mapped into a line, and then the complex potentials are educed. The Muskhelishvili method is extended by creating a Plemelj function that matches the singularity of the real crack tips, and eliminates the pseudo tips’ singularity induced by the conformal mapping. The stress field is obtained after solving the Riemann–Hilbert boundary value problem. Based on the stress field expressions, crack tip stress intensity factors, dislocation dipole image forces and image torque are formulated. Numerical curves show that both the translation and rotation must be considered in the static equilibrium of the dipole system. The crack tip stress intensity factor induced by the dipole may rise or drop and the crack may attract or reject the dipole. These trends depend not only on the crack length, but also on the dipole location, the length and the angle of the dipole span. Generally, the horizontal image force exerted at the center of the dislocation dipole is much smaller than the vertical one. Whether the dipole subjected to clockwise torque or anticlockwise torque is determined by whether the Burgers vector of the crack-nearby dislocation of the dipole is positive or negative. A concentrated load induces no singularity to crack tip stress fields as the load is located at the crack line. However, as the concentrated force is not located on the crack line but approaches the crack tip, the nearby crack tip stress intensity factor K_IIIu increases steeply to infinity.     

纤维复合材料损伤过程的数值模拟

利用界面断裂力学和有限元法数值模拟纤维增强复合材料的细观损伤过程，研究各种主要破坏模式之间的相互转变和影响，指出以断裂能和混合度表示的界面性能是控制复合材料损伤过程的主要细观参数。分析了界面韧度对破坏性能的影响，探讨了基于破坏模式控制的复合材料韧度设计的新途径。     

Thermal stress intensity factors for a crack in an anisotropic half plane

On the basis of the two-dimensional theory of anisotropic thermoelasticity, a solution is given for the thermal stress intensity factors due to the obstruction of a uniform heat flux by an insulated line crack in a generally anisotropic half plane. The crack is replaced by continuous distributions of sources of temperature discontinuity and dislocations. First, the particular thermoelastic dislocation solutions for the half plane are obtained; then the corresponding isothermal solutions are superposed to satisfy the traction-free conditions on the crack surfaces. The dislocation solutions are applied to calculate the thermal stress intensity factors, which are validated by the exact solutions. The effects of the uniform heat flux, the ply angle and the crack length are investigated.     

Stress intensity factor analysis of a three-dimensional interfacial corner between anisotropic bimaterials under thermal stress

A numerical method using a path-independent H-integral based on the conservation integral was developed to analyze the singular stress field of a three-dimensional interfacial corner between anisotropic bimaterials under thermal stress. In the present method, the shape of the corner front is smooth. According to the theory of linear elasticity, asymptotic stress near the tip of a sharp interfacial corner is generally singular as a result of a mismatch of the materials’ elastic constants. The eigenvalues and the eigenfunctions are obtained using the Williams eigenfunction method, which depends on the anisotropic materials’ properties and the geometry of an interfacial corner. The order of the singularity related to the eigenvalue is real, complex or power-logarithmic. The amplitudes of the singular stress terms can be calculated using the H-integral. The stress and displacement around an interfacial corner for the H-integral are obtained using finite element analysis. In this study, a proposed definition of the stress intensity factors of an interfacial corner, which includes those of an interfacial crack and a homogeneous crack, is used to evaluate the singular stress fields. Asymptotic solutions of stress and displacement around an interfacial corner front are uniquely obtained using these stress intensity factors. To prove the accuracy of the present method, several different kinds of examples are shown such as interfacial corners or cracks in three-dimensional structures.     

Stress intensity factors around a cylindrical crack in an interfacial zone in composite materials

Axisymmetric stresses around a cylindrical crack in an interfacial cylindrical layer between an infinite elastic medium with a cylindrical cavity and a circular elastic cylinder made of another material have been determined. The material constants of the layer vary continuously from those of the infinite medium to those of the cylinder. Tension surrounding the cylinder and perpendicular to the axis of the cylinder is applied to the composite materials. To solve this problem, the interfacial layer is divided into several layers with different material properties. The boundary conditions are reduced to dual integral equations. The differences in the crack faces are expanded in a series so as to satisfy the conditions outside the crack. The unknown coefficients in the series are solved using the conditions inside the crack. Numerical calculations are performed for several thicknesses of the interfacial layer. Using these numerical results, the stress intensity factors are evaluated for infinitesimal thickness of the layer.