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.     

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.     

Fatigue crack growth induced by domain switching under electromechanical load in ferroelectrics

Reliability calls for a better understanding of the failure of ferroelectric ceramics. The fracture and fatigue of ferroelectric ceramics under an electric field or a combined electric and mechanical loading are investigated. The small-scale domain-switching model is modified to analyze failure due to fracture and fatigue. Effects of anisotropy and electromechanical load coupling are taken into account. Analytical expressions are obtained for domain-switching regions near the crack tip such that of 90° domain switching can be distinguished from 180° domain switching in addition to different initial poling directions. The crack tip stress intensity variation of ferroelectric ceramics due to the domain switching is analyzed. A positive electric field tends to enhance the propagation of an insulating crack perpendicular to the poling direction, while a negative field impedes it. Fatigue crack growth under various coupling loads and effects of the stress field and electric field on near field stress intensity variation are analyzed. Predicted crack growth versus cyclic electric field agrees well with experiment.     

Pre-kinking of a moving crack in a magnetoelectroelastic material under in-plane loading

A constant moving crack in a magnetoelectroelastic material under in-plane mechanical, electric and magnetic loading is studied for impermeable crack surface boundary conditions. Fourier transform is employed to reduce the mixed boundary value problem of the crack to dual integral equations, which are solved exactly. Steady-state asymptotic fields near the crack tip are obtained in closed form and the corresponding field intensity factors are expressed explicitly. The crack speed influences the singular field distribution around the crack tip and the effects of electric and magnetic loading on the crack tip fields are discussed. The crack kinking phenomena is investigated using the maximum hoop stress intensity factor criterion. The magnitude of the maximum hoop stress intensity factor tends to increase as the crack speed increases.     

起波钢筋高速动态拉伸力学性能研究

针对提高混凝土梁抗冲击爆炸性能的起波配筋新技术,采用理论分析与动态冲击拉伸试验相结合的方法,揭示了起波钢筋的快速拉伸变形作用机理,分析了拉伸速度、起波矢高等因素对起波钢筋抗拉强度的影响规律,确定了起波钢筋静态弹性极限强度计算方法。提出了起波钢筋等效拉伸应变率新概念,建立了弹性极限强度动力放大系数（dynamic increase factors,DIF）计算模型。研究结果表明,预先弯折起波使得钢筋在受力拉直过程中产生截面弯矩,起波钢筋的力学性能存在明显的应变率效应;起波钢筋弹性抗拉极限强度DIF随起波矢高的增高先增大后减小,存在一个最优起波设计矢高,可以使起波钢筋抗拉强度动态放大系数达到最大。研究成果可为进一步推动起波配筋技术在防护工程中的应用提供依据。     

Anti-plane analysis of semi-infinite crack in piezoelectric strip

Using the complex variable function method and the technique of the conformal mapping, the fracture problem of a semi-infinite crack in a piezoelectric strip is studied under the anti-plane shear stress and the in-plane electric load. The analytic solutions of the field intensity factors and the mechanical strain energy release rate are presented under the assumption that the surface of the crack is electrically impermeable. When the height of the strip tends to infinity, the analytic solutions of an infinitely large piezoelectric solid with a semi-infinite crack are obtained. Moreover, the present results can be reduced to the well-known solutions for a purely elastic material in the absence of the electric loading. In addition, numerical examples are given to show the influences of the loaded crack length, the height of the strip, and the applied mechanical/electric loads on the mechanical strain energy release rate.     

Exact solutions of two semi-infinite collinear cracks in piezoelectric strip

Using the complex variable function method and the conformal mapping technique,the fracture problem of two semi-infinite collinear cracks in a piezoelectric strip is studied under the anti-plane shear stress and the in-plane electric load on the partial crack surface.Analytic solutions of the field intensity factors and the mechanical strain energy release rate are derived under the assumption that the surfaces of the crack are electrically impermeable.The results can be reduced to the well-known solutio...     

深埋椭圆形片状裂纹的偏折扩展

基于无限大弹性基体深埋椭圆形片状裂纹的变形场,推导了椭圆形片状裂纹的能量释放率,采用能量平衡方法建立了椭圆形片状裂纹承受拉应力和剪应力时的复合断裂准则. 考虑裂纹在拉-剪应力作用下的偏折扩展,分析了裂纹的偏折方向,提出了椭圆形片状裂纹发生偏折扩展时的初始偏折位置的确定方法.      

The kinked crack solved by Mellin transform

The problem of the evaluation of the stress intensity factors (SIF) at a kinked crack-tip is treated by the use of the most appropriate tool for this problem, the Mellin transform technique. The advantage of this technique lies on the validity of the solution everywhere in the stress field and for every length of the branched crack. Since in the limiting case of the kinked crack, numerical difficulties appear in the evaluation of the kernel functions of the corresponding Fredholm integral-equation system, an asymptotic development, with respect to the length of the kinked crack, is constructed for these functions and therefore for the SIFs at the kinked tip. This analysis yields the interesting result that the SIFs at the crack tip tend to vanish as the kinked crack reduces in length, a result which is similar to the already established behavior of the SIFs at an interfacial tip.     

Anti-plane shear of kinked interface crack in bonded dissimilar anisotropic solid

Complex potentials are derived to describe the anti-plane singular shear stress fields around a kinked crack, the main portion of which is embedded along the interface of two dissimilar anisotropic elastic media. This is accomplished by formulating the problem as singular integral equations with generalized Cauchy kernels. The shear stress singularity at the kink differs from the familiar inverse square root of the local distance; it is found to influence the magnitude of the Mode III crack tip stress intensity factor, K₃. Numerical results of K₃ are obtained and displayed in graphical forms for different degree of material anisotropy and crack dimensions.     

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.     

Strip-electromechanical model solution for piezoelectric plate cut along two semi-permeable collinear cracks

A strip electric saturation and mechanical yielding model solution is proposed for a piezoelectric plate cut along two equal collinear semi-permeable mode-I cracks with electrical polarization reaching a saturation limit and normal stress reaching a yield stress along a line segment in front of the cracks. By using Stroh formalism and complex variable technique, we derived the analytical solution for the field quantities. Three different situations are investigated when developed electrical saturation zone is bigger/smaller or equal to the developed mechanical yield zone. Numerical results show that the effect of different electric boundary conditions on the crack opening displacement and crack opening potential drop is significant and should not be ignored. The influence of electric load displacement on the energy release rate is also investigated for PZT-4, PZT-5H and BaTiO₃ ceramics, and it may assists for the correct choosing of ceramic for specific job.     

Modeling of pre-fracture zones for limited permeable crack in piezoelectric materials

A plane-strain problem for a limited permeable crack in an adhesive thin interlayer between two semi-infinite piezoelectric spaces is considered. The tensile mechanical stress and the electric displacement are applied at infinity. The interlayer is assumed to be softer than the connected materials; therefore, the zones of mechanical yielding and electric saturations can arise at the crack tips on the continuations of the crack. These zones are considered in this work. It was assumed that the length of electric saturation zones is larger than the length of mechanical yield zones. The zones of mechanical yielding are modeled by the crack continuations with normal compressive stresses applied at its faces. The electric saturation zones are modeled by segments at the crack continuations with prescribed saturated electric displacements. These electric displacements can linearly vary along the mechanical yielding zones. The problem is reduced to the Hilbert–Riemann problem of linear relationship, which is solved exactly. The equation for the determination of the yielding zones length, the expressions for the crack-opening displacement jump, electric potential jump, and J-integral is obtained in an analytical form. In case of finite size body, the finite elements method is used and the variation in the fracture mechanical parameters with respect to this size is demonstrated.     

An interface crack with partially electrically conductive crack faces under antiplane mechanical and in-plane electric loadings

An interface crack in a bimaterial piezoelectric space under the action of antiplane mechanical and in-plane electric loadings is analyzed. One zone of the crack faces is electrically conductive while the other part is electrically permeable. All electro-mechanical values are presented using sectionally-analytic vector-functions and a combined Dirichlet-Riemann boundary value problem is formulated. An exact analytical solution of this problem is obtained. Simple analytical expressions for the shear stress, electric field and also for mechanical displacement jump of the crack faces are derived. These values are also presented graphically along the corresponding parts of the material interface. Singular points of the shear stress, electric field and electric displacement jump are found. Their intensity factors are determined as well. Intensity factors variations with respect to the external electric field and different ratios between the electrically conductive and electrically permeable crack face zones are also demonstrated.     

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.     

Closed form solution for two unequal collinear semi-permeable straight cracks in a piezoelectric media

The problem of two unequal collinear straight cracks weakening a poled transversely isotropic piezoelectric ceramic is addressed under semi-permeable electric boundary conditions on the crack faces. The plate has been subjected to combined in-plane normal(to the faces of the cracks) mechanical and electric loads. Problem is formulated employing Stroh formalism and solved using complex variable technique. The elastic field, electric field and energy release rate are obtained in closed analytic form. A case study is presented for poled PZT-5H cracked plate to study the effect of prescribed mechanical load, electric load, inter-crack distance and crack lengths on crack arrest parameters stress intensity factor (SIF), electric displacement intensity factor (EDIF) and mechanical and total energy release rates (ERR). Moreover a comparative study is done of impermeable and semi-permeable crack face boundary conditions on SIF, EDIF and ERR, and results obtained is presented graphically. It is observed that the effect of dielectric medium in the crack gap cannot be ignored.     

A mode III crack crossing the magnetoelectroelastic bimaterial interface under concentrated magnetoelectromechanical loads

A mode III crack cutting perpendicularly across the interface between two dissimilar semi-infinite magnetoelectroelastic solid is studied under the combined loads of a line force, a line electric charge and a line magnetic charge at an arbitrary location. The impermeable conditions are implied on the crack faces. The technique developed in literature for the elastic bimaterial with a crack cutting interface is exploited to treat the magnetoelectroelastic bimaterial. The Riemann-Hilbert problem can be formulated and solved based on complex variable method. Analytical solutions can be obtained for the entire plane. The intensity factors around crack tips can be defined for the elastic, electric and magnetic fields. It shows that, no matter where the load position is, the electric displacement intensity factors (EDIFs), as well as the magnetic induction intensity factors (MIIFs), are identical in magnitude but opposite in sign for both crack tips, on condition that a line force is solely applied. Alternatively, if only a line electric charge is considered, then the stress intensity factors (SIFs) and the MIIFs exhibit the behavior. Likewise, if only a line magnetic charge is applied, it turns to the SIFs and the EDIFs instead. In addition, the dependence of the intensity factors is graphically shown with respect to the location of a line force. It is found that the SIF for a crack tip tends to be infinite if the applied force is approaching the tip itself, but the EDIF, with the complete opposite trend, tends to be vanishing. Finally, focusing on the more practical case of piezoelectric/piezomagnetic bimaterial, variation of the SIF along with the moduli as well as the piezo constitutive coefficients is explored. These analyses may provide some guidance for material selection by minimizing the SIF. It is also believed that the results obtained in this paper can serve as the Green’s function for the dissimilar magnetoelectroelastic semi-infinite bimaterial with a crack cutting the interface under general magnetoelectromechanical loads.     

Crack kinking in piezoelectric materials

A solution is presented for a class of two-dimensional electroelastic branched crack problems. Explicit Green's function for an interface crack subject to an edge dislocation is developed using the extended Stroh formulation allowing the branched crack problem to be expressed in terms of coupled singular integral equations. The integral equations are obtained by the method that models a kink as a continuous distribution of edge dislocations, and the dislocation density function is defined on the line of the branch crack only. Competition between crack extension along the interface and kinking into the substrate is investigated using the integral equations and the maximum energy release rate criterion. Numerical results are presented to show the effect of electric field on the path of crack extension. The work was supported by the Australian Research Council through a Queen Elizabeth II fellowship and by the Australian Academy of Science through the J.G. Russell Award.     

拉伸荷载下分支裂隙破坏机理研究

为研究拉伸荷载下分支裂隙对破坏模式的影响,保持主裂隙参数不变,改变分支裂隙倾角和长度,利用扩展有限元方法模拟了弯折裂隙的动态扩展,总结了分支裂隙参数变化对破坏模式的影响。利用ABAQUS中的轮廓积分计算了分支裂隙尖端应力强度因子,并根据最大周向应力准则计算起裂角。结果表明:拉伸荷载下分支裂隙出现三种破坏模式;分支裂隙倾角和长度均对破坏模式有一定的影响。I型应力强度因子与分支裂隙倾角关系曲线呈斜"S"型,相应II型应力强度因子曲线呈上凸型;由于分支裂隙存在非尖端破坏,利用裂隙尖端应力强度因子判断开裂应结合相应的破坏模式。     

Elastoplastic solution for an eccentric crack loaded by two pairs of point tensile forces

The strain energy density theory and the near crack line analysis method are applied to investigate an eccentric crack loaded by two pairs of tensile point forces in a finite plate. The minimum values of SED in the vicinity of the crack tip are determined, the initial growth orientation of crack are determined. Obtained is the elastic-plastic solution near the crack line of an eccentric crack loaded by two pairs of point tensile forces under large scale yielding condition. More specifically, the near field solution contains the unit normal vector of the elastic-plastic boundary and the elastic-plastic stress field. The length of the plastic zone along the crack line is found to vary with the external load and the bearing capacity of a finite plate with an eccentric crack loaded by two pairs of tensile point forces. Compared with small scale yielding condition, the normalized load obtained is higher than those under small scale yielding condition when the length of the plastic zone is the same.