Conservation Law and Application of J-Integral in Multi-Materials

The conservation law of J-integral in two-media with a crack paralleling to the interface of the two media was firstly proved by analytical and numerical finite element method. Then a schedule model was established that an interface crack is inserted in four media. According to the J-integral conservation law on multi-media, the energy release ratio of Ⅰ-type crack was considered to be conservation when the middle medium layers are very thin. And the conservation law was also convinced by numerical method. By means of the dimension analysis on the model, the asymptotic results and formula calculating the energy release ratio and complex stress intensity factor are presented.     

涂层式裂纹监测系统中基体裂纹穿越行为研究

建立了智能涂层的两相模型与三相模型,基于能量准则分别用这两种模型研究了基体裂纹达到涂层界面后的穿越/偏转行为. 用有限元法分析了相对裂纹扩展长度、弹性错配参数及界面层厚度对偏转裂纹与穿越裂纹能量释放率之比的影响,结果表明当基体裂纹到达驱动层与基体界面时,能量释放率之比不仅与基体和驱动层之间的弹性错配相关,而且当驱动层较薄时对驱动层与传感层之间的弹性错配亦有较强的依赖性. 此外,随着驱动层厚度的增加,能量释放率之比对驱动层与传感层之间的弹性错配的依赖性逐渐降低. 通过与实验结果相比,建立的模型能够较好的解释基体裂纹在界面的扩展行为,可用于智能涂层裂纹传感器的优化设计.      

Numerical solution of polarization saturation/dielectric breakdown model in 2D finite piezoelectric media

The polarization saturation (PS) model [Gao, H., Barnett, D.M., 1996. An invariance property of local energy release rates in a strip saturation model of piezoelectric fracture. Int. J. Fract. 79, R25–R29; Gao, H., Zhang, T.Y., Tong, P., 1997. Local and global energy release rates for an electrically yielded crack in a piezoelectric ceramic. J. Mech. Phys. Solids 45, 491–510], and the dielectric breakdown (DB) model [Zhang, T.Y., Zhao, M.H., Cao, C.F., 2005. The strip dielectric breakdown model. Int. J. Fract. 132, 311–327] explain very well some experimental observations of fracture of piezoelectric ceramics. In this paper, the nonlinear hybrid extended displacement discontinuity-fundamental solution method (NLHEDD-FSM) is presented for numerical analysis of both the PS and DB models of two-dimensional (2D) finite piezoelectric media under impermeable and semi-permeable electric boundary conditions. In this NLHEDD-FSM, the solution is expressed approximately by a linear combination of fundamental solutions of the governing equations, which includes the extended point force fundamental solutions with sources placed at chosen points outside the domain of the problem under consideration, and the extended Crouch fundamental solutions with extended displacement discontinuities placed on the crack and the electric yielding zone. The coefficients of the fundamental solutions are determined by letting the approximated solution satisfy certain conditions on the boundary of the domain, on the crack face and the electric yielding zone. The zero electric displacement intensity factor in the PS model or the zero electric field strength intensity factor in the DB model at the outer tips of the electric yielding zone is used as a supplementary condition to determine the size of the electric yielding zone. Iteration approaches are adopted in the NLHEDD-FSM. The electric yielding zone is determined, and the extended intensity factors and the local J-integral are calculated for center cracks in piezoelectric strips. The effects of finite domain size, saturation property and different electric boundary conditions, as well as different models on the electric yielding zone and the local J-integral, are studied.     

A NEW CYCLIC J-INTEGRAL FOR LOW-CYCLE FATIGUE CRACK GROWTH

The constitutive equation under the low-cycle fatigue (LCF) was discussed, and a two-dimensional (2-D) model for simulating fatigue crack extension was put forward in order to propose a new cyclic J-integral. The definition, primary characteristics, physical interpretations and numerical evaluation of the new parameter were investigated in detail. Moreover, the new cyclic J-integral for LCF behaviors was validated by the compact tension (CT) specimens. Results show that the calculated values of the new parameter can correlate well with LCF crack growth rate, during constant-amplitude loading. In addition, the phenomenon of fatigue retardation was explained through the viewpoint of energy based on the concept of the new parameter.     

Energy release rate and contact zone in a cohesive and an interface crack by hypersingular integral equations

Solution of Cauchy-type singular integral equations permits the evaluation of the fracture parameters at the crack tips very accurately. However, it does not permit the determination of the crack opening and sliding displacements while ensuring no crack surface interpenetration unless the location of the contact zone is known a priori. In order to circumvent this shortcoming, this study presents a solution method based on the Hadamard-type singular integral equations to obtain the crack opening and sliding displacements directly while enforcing the appropriate conditions to prevent interpenetration. Furthermore, the crack opening displacements are physically more meaningful and readily validated against the finite element analysis predictions. The numerical solutions of the hypersingular integral equations provide not only crack opening and sliding displacements directly but also the stress intensity factors and energy release rates. Also, the behavior of the energy release rate is examined as the cohesive crack located parallel to the interface approaches the interface from either the soft or the stiff side of the interface. The limiting value of the energy release rate is established by considering an interface crack. As the cohesive crack approaches the interface from either side of the interface, the energy release rate approaches to that of the interface crack. However, the length of contact zone between the cohesive crack surfaces under uniform shear loading does not approach to that of the interface crack.     

Interaction between an interface crack and subinterface microcracks in metal/piezoelectric bimaterials

The J-integral analysis is presented for the interaction problem between a semi-infinite interface crack and subinterface matrix microcracks in dissimilar anisotropic materials. After deriving the fundamental solutions for an interface crack subjected to different loads and the fundamental solutions for an edge dislocation beneath the interface, the interaction problem is deduced to a system of singular integral equations with the aid of a superimposing technique. The integral equations are then solved numerically and a conservation law among three values of the J-integral is presented, which are induced from the interface crack tip, the microcracks and the remote field, respectively. The conservation law not only provides a necessary condition to confirm the numerical results derived, but also reveals that the microcrack shielding effect in such materials could be considered as a redistribution of the remote J-integral. It is this redistribution that does lead to the phenomenological shielding effect.     

Driving forces and kinking of an oblique crack in bonded nonhomogeneous materials

Summary  Plane elasticity solutions are presented for the problem of an oblique crack in two bonded media. The material model under consideration consists of a homogeneous half-plane with an arbitrarily oriented crack and a nonhomogeneous half-plane. The Fourier integral transform method is employed in conjunction with the coordinate transformations of field variables in the basic elasticity equations. Formulation of the crack problem results in having to solve a system of singular integral equations for arbitrary crack surface tractions. A crack perpendicular to or along the bonded interface between the homogeneous and nonhomogeneous constituents arises as a limiting case. In the numerical results, the values of mixed-mode stress intensity factors are provided for various combinations of relevant geometric and material parameters of the bonded media. Subsequently, the infinitesimal kinks from the tips of a main crack are presumed, with the corresponding local driving forces being evaluated in terms of the stress intensities of the main crack. The criterion of maximum energy release rate is applied with the aim of making some conjectures concerning the likelihood of kinking and the probable kink direction based on the approximation of local homogeneity and brittleness of the crack-tip behavior. Received 25 September 2001; accepted for publication 13 February 2002     

基于构型力内变量的界面损伤问题研究

本文基于材料构型力的基本理论和损伤力学中含有内变量的热力学框架,提出了新的损伤变量定义方式,为研究界面损伤问题提供了一种新思路。首先,基于双相弹性体的能量分析,给出界面材料构型力表达式,通过构型力的离散化方法,实现了其在有限元中的数值计算。其次,定义构型力为界面损伤内变量,进而提出一种新的损伤演化模型,并采用刚度劣化的方法,对该界面损伤模型进行数值实现。最后,通过对复合材料界面损伤问题（有裂纹或无裂纹）进行数值模拟,分析了其界面损伤发展趋势,探讨了此模型的合理性和优越性。基于构型力内变量的界面损伤模型,可为研究复合材料的界面损伤失效问题提供一种普适性的方法。     

Boundary element analysis of cracked film–substrate media

This study evaluates the stress behavior of a cracked film–substrate medium by applying the multi-region boundary element method. Four problems addressed herein are the crack tip within a film, the crack tip terminating at the interface, interface debonding, and the crack penetrating into the substrate. The multi-region boundary element method is initially developed and, then, the stress intensity factors or the energy release rates are evaluated according to the different stress singularities of the four considered problems. These results indicate that the stress intensity factors or the energy release rates of the four problems rely not only on the different elastic mismatches and crack lengths, but also on the thickness ratio of the film and the substrate.     

Investigation on behavior of crack penetration/deflection at interfaces in intelligent coating system

Based on the three-phase model, the propagation behavior of a matrix crack in an intelligent coating system is investigated by an energy criterion. The effect of the elastic mismatch parameters and the thickness of the interface layer on the ratio of the energy release rate for infinitesimal deflected and penetrated crack is evaluated with the finite element method. The results show that the ratio of the energy release rates strongly depends on the elastic mismatch α1between the substrate and the driving layer.It also strongly depends on the elastic mismatch α2between the driving layer and the sensing layer for a thinner driving layer when a primary crack reaches an interface between the substrate and the driving layer. Moreover, with the increase in the thickness of the driving layer, the dependence on α2gradually decreases. The experimental observation on aluminum alloys monitored with intelligent coating shows that the established model can better explain the behavior of matrix crack penetration and can be used in optimization design of intelligent coating.     

A numerical method for crack growth analysis in linearly viscoelastic materials

This paper presents a numerical method for the analysis of crack initiation and extension in linearly viscoelastic materials undergoing Mode I plane stress deformation. Plastic deformation near the crack tip is considered by a strip-yielding model. The crack initiation and growth are taken to follow a critical energy release rate fracture criterion, and the plastic work per unit crack extension is included in the calculation of the total energy release rate. The resulting numerical method is presented in algorithmic form and two example problems are solved to demonstrate its application. The two example problems represent unstable and stable crack propagations respectively. The results obtained lend insight into the effect of creep deformation on the initiation and growth of a crack.     

Interface crack between isotropic Kirchhoff plates

In this work Kirchhoff plate theory is used to calculate the energy release rate function in delaminated isotropic plates. The approximation is based on the consideration of the equilibrium equations and the displacement continuity between the interface plane of a double-plate model. It is shown that the interface shear stresses are governed by a fourth order partial differential equation system. As an example, a simply supported delaminated plate subjected to a point force is analyzed adopting Lévy plate formulation and the mode-II and mode-III energy release rate distributions along the crack front were calculated by the J-integral. To confirm the analytical results the 3D finite element model of the delaminated plate was created, the energy release rates were calculated by the virtual crack-closure technique and the J-integral. The results indicate a good agreement between analysis and numerical computation.     

Crack deflection at an interface between dissimilar elastic materials: Role of residual stresses

A crack intersecting an interface between two dissimilar materials may advance by either penetrating through the interface or deflecting into the interface. The competition between deflection and penetration can be assessed by comparison of two ratios: (i) the ratio of the energy release rates for interface cracking and crack penetration; and (ii) the ratio of interface to material fracture energies. Residual stresses caused by thermal expansion misfit can influence the energy release rates of both the deflected and penetrating crack. This paper analyses the role of residual stresses. The results reveal that expansion misfit can be profoundly important in systems with planar interfaces (such as layered materials, thin film structures, etc.), but generally can be expected to be of little significance in fiber composites. This paper corrects an earlier result for the ratio of the energy release rate for the doubly deflected crack to that for the penetrating crack in the absence of residual stress.     

弹性T项对裂尖参数和裂纹扩展路径稳定性的作用

研究了弹性T项在主裂纹与近裂尖微空洞干涉问题中对裂尖参数和裂纹扩展路径稳定性的影响作用.利用“伪力”法,并通过解决主裂纹与近裂尖微空洞干涉问题,对远场纯1型载荷和远场弹性T项载荷下,该干涉问题中弹性T项的影响作用进行了分析从数值结果可以看出：由于空洞的存在;释放了弹性T项所引起的应力,弹性T项对裂尖参数;应力强度因子和J积分都有直接显著的影响,因而,它对该载荷下的裂纹扩展路径的稳定性有控制作用。     

MESHLESS METHOD FOR 2D MIXED-MODE CRACK PROPAGATION BASED ON VORONOI CELL

A meshless method integrated with linear elastic fracture mechanics (LEFM) is presented for 2D mixed-mode crack propagation analysis. The domain is divided automatically into sub-domains based on Voronoi cells, which are used for quadrature for the potential energy. The continuous crack propagation is simulated with an incremental crack-extension method which assumes a piecewise linear discretization of the unknown crack path. For each increment of the crack extension, the meshless method is applied to carry out a stress analysis of the cracked structure. The J-integral, which can be decomposed into mode I and mode II for mixed-mode crack, is used for the evaluation of the stress intensity factors (SIFs). The crack-propagation direction, predicted on an incremental basis, is computed by a criterion defined in terms of the SIFs. The flowchart of the proposed procedure is presented and two numerical problems are analyzed with this method. The meshless results agree well with the experimental ones, which validates the accuracy and efficiency of the method.     

The derivation of J-integral conservation law from the principle of virtual work

The general form of J-integral conservation law in three di-mensions has been derived from the principle of virtual work.The necessary conditions for J-integral conservation are alsogiven.     

固体炸药爆轰与惰性介质相互作用的一种扩散界面模型

提出一种保持热力学一致性的扩散界面模型，用来数值模拟固体炸药爆轰与惰性介质的相互作用问题。基于混合网格内各组分物质间可以达到力学平衡状态而不能达到热学平衡状态的假设，由混合网格能量守恒以及压力相等条件，推导出每种组分物质的体积分数演化方程。由此获得的扩散界面模型包括组分物质的质量守恒方程、混合物质的动量及总能量守恒方程，同时包括组分物质的体积分数演化方程和混合物质的压力演化方程。该扩散界面模型的主要特点是考虑了化学反应以及热学非平衡的影响。提出的扩散界面模型在物质界面附近不会出现物理量的非物理振荡现象、适用于任意表达形式的物质状态方程以及任意数目的惰性介质。