热对流作用下筒壁涂层的边裂行为

边裂(边缘开裂)是涂层热致损伤的主要模式之一. 边缘裂纹穿透涂层后,常导致界面脱粘从而驱使涂层与基体剥离,最终丧失对基体的保护作用. 本文以热应力强度因子表征边缘裂纹的扩展驱动力,研究筒壁涂层在热对流作用下的边裂行为. 首先,利用拉普拉斯变换法,得到了瞬态温度场及热应力场的封闭解. 其次,运用Fett等的三参数法确定了筒壁涂层边缘裂纹的权函数. 最后,基于叠加原理和权函数方法计算了边缘裂纹的热应力强度因子. 探讨了无量纲时间、边缘裂纹深度、基体/涂层厚度比、热对流强度等参数对热应力强度因子的影响规律. 结果表明:热应力强度因子的峰值既非发生在热载荷初始时刻,也非发生在热稳态时刻,而出现在时间历程的中间时刻;增大热对流强度不仅可提高热应力强度因子的峰值,而且使峰值提前出现;其他条件相同时,热应力强度因子随着边缘裂纹长度的增大而降低;增大涂层厚度或减小基体厚度可增强涂层抵抗瞬态热载荷的能力.      

Crack propagation in glass coatings under expanding spherical contact

The growth of transverse cracks under expanding spherical contact in a model system consisted of soda-lime glass bonded to a polycarbonate substrate is observed in situ from below or from the polished edge of the bilayer. Abrasion or chemical etching is employed on the coating surfaces to control the initial fracture. In the limit case of monoliths, the crack mouth becomes fully engulfed by the expanding contact, which results in a much steeper crack angle compared to the classical Hertzian cone case. As the coating thickness is reduced, flexure stresses are set in the coating which drive the cone crack to well away from the contact circle and initiate semi-elliptical-like radial cracks at the subsurface, right under the contact. Common to all three fracture modes is an initial unstable propagation phase following by a stable growth, with detrimental failure associated with severe damage to the top surface and/or delamination at the coating/substrate interface taking place at loads several times the fracture initiation loads.LEFM in conjunction with a large-strain FEM contact code is used to study the post-initiation fracture, with the crack path controlled by the principal stress trajectory or zero-mode II S.I.F. The analysis exposes the leading geometric and material parameters in each fracture mode, which may be useful in the design of bilayer structures for optimal mechanical performance. The well-known Auerbach law governing the initial fracture of monoliths is found to apply also to the bilayer crack systems within a certain range of the problem parameters. The numerical prediction for the crack profiles and the fracture envelopes generally collaborate well with the tests.     

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.     

Debonding of graded coatings under in-plane compression

Graded materials are multiphase composites with continuously varying thermophysical properties. The concept provides material scientists and engineers with an important tool to develop new materials tailored for some specific applications. One such application of this new class of materials is as top coats or interfacial regions in thermal barrier systems. A widely observed failure mode in these layered materials is known to be interfacial cracking that leads to spallation. In many cases it is the buckling instability of coating under mechanically or thermally induced compressive stresses that triggers spallation. Under in-plane loading since the linear elastic small deformation theory gives only a trivial solution, in this study the plane strain interface crack problem for a graded coating bonded to a homogeneous substrate is formulated by using a kinematically nonlinear continuum theory. Both the instability and the postbuckling problems are considered. The main objective of the study is the investigation of the influence of material nonhomogeneity, kinematic nonlinearity and plate approximation on the critical instability load and on such fracture mechanics parameters as strain energy release rate, stress intensity factors and crack opening displacements.     

Periodic cracks in a functionally graded piezoelectric layer bonded to a piezoelectric half-plane

Studied is the problem of a periodic array of cracks in a functionally graded piezoelectric strip bonded to a homogeneous piezoelectric material. The properties of the functionally graded piezoelectric strip, such as elastic modulus, piezoelectric constant and dielectric constant, are assumed in exponential forms and vary along the crack direction. The crack surface condition is assumed to be electrically impermeable or permeable. Integral transform and dislocation density functions are employed to reduce the problem to the solution of a system of singular integral equations. The effects of the periodic crack spacing, material constants and the geometry parameters on the stress intensity factor, the energy release ratio and the energy density factor are studied.     

Buckling of an orthotropic graded coating with an embedded crack bonded to a homogeneous substrate

To simulate buckling of nonuniform coatings, we consider the problem of an embedded crack in a graded orthotropic coating bonded to a homogeneous substrate subjected to a compressive loading. The coating is graded in the thickness direction and the material gradient is orthogonal to the crack direction which is parallel with the free surface. The elastic properties of the material are assumed to vary continuously along the thickness direction. The principal directions of orthotropy are parallel and perpendicular to the crack orientation. The loading consists of a uniform compressive strain applied away from the crack region. The graded coating is modeled as a nonhomogeneous medium with an orthotropic stress–strain law. Using a nonlinear continuum theory and a suitable perturbation technique, the plane strain problem is reduced to an eigenvalue problem describing the onset of buckling. Using integral transforms, the resulting plane elasticity equations are converted analytically into singular integral equations which are solved numerically to yield the critical buckling strain. The Finite Element Method was additionally used to model the crack problem. The main objective of the paper is to study the influence of material nonhomogeneity on the buckling resistance of the graded layer for various crack positions, coating thicknesses and different orthotropic FGMs.     

Multiple axial cracks in a coated hollow cylinder due to thermal shock

The transient thermal stress problem of an inner-surface-coated hollow cylinder with multiple pre-existing surface cracks contained in the coating is considered. The transient temperature, induced thermal stress, and the crack tip stress intensity factor (SIF) are calculated for the cylinder via finite element method (FEM), which is exposed to convective cooling from the inner surface. As an example, the material pair of a chromium coating and an underlying steel substrate 30CrNi2MoVA is particularly evaluated. Numerical results are obtained for the stress intensity factors as a function of normalized quantities such as time, crack length, convection severity, material constants and crack spacing.     

An axisymmetric problem of an embedded crack in a graded layer bonded to a homogeneous half-space

In an attempt to simulate non-uniform coating delamination, the elasto-static problem of a penny shaped axisymmetric crack embedded in a functionally graded coating bonded to a homogeneous substrate subjected to crack surface tractions is considered. The coating’s material gradient is parallel to the axisymmetric direction and is orthogonal to the crack plane. The graded coating is modeled as a non-homogeneous medium with an isotropic constitutive law. Using Hankel transform, the governing equations are converted into coupled singular integral equations, which are solved numerically to yield the crack tip stress intensity factors. The Finite Element Method was additionally used to model the crack problem. The main objective of this paper is to study the influence of the material non-homogeneity and the crack position on the stress intensity factors for the purpose of gaining better understanding on the behavior of graded coatings.     

Dynamic crack growth under Rayleigh wave

A nonuniform crack growth problem is considered for a homogeneous isotropic elastic medium subjected to the action of remote oscillatory and static loads. In the case of a plane problem, the former results in Rayleigh waves propagating toward the crack tip. For the antiplane problem the shear waves play a similar role. Under the considered conditions the crack cannot move uniformly, and if the static prestress is not sufficiently high, the crack moves interruptedly. For fracture modes I and II the established, crack speed periodic regimes are examined. For mode III a complete transient solution is derived with the periodic regime as an asymptote. Examples of the crack motion are presented. The crack speed time-period and the time-averaged crack speeds are found. The ratio of the fracture energy to the energy carried by the Rayleigh wave is derived. An issue concerning two equivalent forms of the general solution is discussed.     

Elastic thermal stresses in a hollow circular overlay/substrate system

This paper investigates the thermal elastic fields in the hollow circular overlay fully bonded to a rigid substrate, which is subjected to a temperature change. Following our previous work for a solid circular overlay/substrate system (Yuan and Yin, Mech. Res. Commun. 38, 283–287, 2011), this paper presents a closed form approximate solution to the axisymmetric boundary value problem using the plane assumption, whose accuracy is verified by the finite element models. When the inner radius is reduced to zero, the present solution recovers the previous solution. When the outer radius approaches infinite, the solution provides the elastic fields for a tiny hole in the overlay. The effects of thickness and width of the overlay are investigated and discussed. When a circular crack initiates in a solid circular overlay, the fracture energy release rate is investigated. This solution is useful for thermal stress analysis of hollow circular thin film/substrate systems and for fracture analysis of spiral cracking in the similar structures.     

Thermal shock analyses for a strip containing an infinite row of periodically distributed cracks normal to its edge

The transient thermal stress problem of a semi-infinite plate containing an infinite row of periodically distributed cracks normal to its edge is investigated in this paper. The elastic medium is assumed to be cooled suddenly on the crack-containing edge. By the superposition principle, the formulation leads to a mixed boundary value problem, with the negating tractions arisen from the thermal stresses for a crack-free semi-infinite plate. The resulting singular integral equation is solved numerically. The effects on the stress intensity factors due to the presence of periodically distributed cracks in a semi-infinite plate are illustrated. For both the edge crack and the embedded crack arrays, the stress intensity factors increase, due to the reduction of the shielding effect, as the stacking cracks are more separated. For the case of embedded crack array, one has the further conclusion that the stress intensity factors decline as the crack array shifts from the plate edge.     

Interaction of cracks with bonds between their faces in an isotropic medium reinforced by a regular system of stringers

A problem of an elastic isotropic medium with a system of foreign (transverse with respect to crack alignment) rectilinear inclusions is considered. The medium is assumed to be attenuated by a periodic system of rectilinear cracks with zones where the crack faces interact with each other. These zones are assumed to be adjacent to the crack tips, and their sizes can be commensurable with the crack size. Interaction between the crack faces in the tip zone is modeled by introducing bonds (adhesion forces) between the cracks with a specified strain diagram. The boundary-value problem of the equilibrium of a periodic system of cracks with bonds between their faces under the action of external tensile loads and forces in the bonds is reduced to a nonlinear singular integrodifferential equation with a kernel of the Cauchy kernel type. The condition of critical equilibrium of the cracks with the tip zones is formulated with allowance for the criterion of critical tension of the bonds. A case of a stress state of the medium containing zones where the crack faces interact with each other is considered.     

An embedded crack in a functionally graded orthotropic coating bonded to a homogeneous substrate under a frictional Hertzian contact

In this paper, we consider the elasto-static problem of an embedded crack in a graded orthotropic coating bonded to a homogeneous substrate subject to statically applied normal and tangential surface loading. The crack direction is parallel to the free surface. The coating is graded in the thickness direction and is orthogonal to the crack direction. This coating is modelled as a non-homogeneous medium with an orthotropic stress–strain law. The equivalent crack surface stresses are first obtained and substituted in the plane elasticity equations. Using integral transforms, the governing equations are converted into singular integral equations which are solved numerically to yield the displacement field as well as the crack-tip stress intensity factors. This study presents a complete theoretical formulation for the problem in the static case. A numerical predictive capability for solving the singular integral equations and computing the crack-tip stress intensity factors is proposed. Since the loading is compressive, a previously developed crack-closure algorithm is applied to avoid interpenetration of the crack faces. The main objective of the paper is to investigate the effects of the material orthotropy and non-homogeneity of the graded coating on the crack-tip stress intensity factors, with and without using the crack-closure algorithm, for the purpose of gaining better understanding on the behavior and design of graded coatings.     

热压电介质中的渗透型周期裂纹问题

利用Stroh公式,研究了含共线周期裂纹热的压电介质的广义二维问题。该工作有两个特征：一是裂纹被建模为具有渗透表面的缝隙,并假设为跨越上下表面时,电场的切向分量和电位移的法向分量是连续的;另一个特征是,机－电载荷和热载荷被假设作用在无限远处,而不是在裂纹表面。基于这两个假设,我们获得了有关场强因子,以及裂纹内电场的相当简洁的表达式。结果表明：①在裂纹内电场是线性变化的,②电位移的奇异性总是取决于应力的奇异性.③所有场的奇异性与所加的电载荷无关。     

Crack spacing effect for a piezoelectric cylinder under electro-mechanical loading or transient heating

This paper investigates the fracture problem of a piezoelectric cylinder with a periodic array of embedded circular cracks. An electro-mechanical fracture mechanics model is established first. The model is further used to the thermal fracture analysis of a piezoelectric cylinder subjected to a sudden heating on its outer surface. The temperature field and the associated thermal stresses and electric displacements are obtained and are added to the crack surface to form a mixed-mode boundary value problem for the electro-mechanical coupling fracture. The stress and stress intensities are investigated for the effect of crack spacing. Strength evaluation of piezoelectric materials under the transient thermal environment is made and thermal shock resistance of the medium is given.     

Stress analysis of an elastic cracked layer bonded to a viscoelastic substrate

The effect of a viscoelastic substrate on an elastic cracked layer under an in-plane concentrated load is solved and discussed in this study. Based on a correspondence principle, the viscoelastic solution is directly obtained from the corresponding elastic one. The elastic solution in an anisotropic trimaterial is solved as a rapidly convergent series in terms of complex potentials via the successive iterations of the alternating technique in order to satisfy the continuity condition along the interfaces between dissimilar media. This trimaterial solution is then applied to a problem of a thin layer bonded to a half-plane substrate. Using the standard solid model to formulate the viscoelastic constitutive equation, the real-time stress intensity factors can be directly obtained by performing the numerical calculations. The results obtained in this paper are useful in studying the problem with bone defects where a crack is assumed to exist in an elastic body made of the cortical bone that is bonded to a viscoelastic substrate made of the cancellous bone.     

Fracture mechanics analysis of thin coatings under plane-strain indentation

A combined experimental/analytical work is carried out to elucidate the fracture resistance of a thin, hard coating bonded to a semi-infinite substrate due to indentation by a cylindrical surface. The bending of the coating under the softer substrate induces concentrated tensile stress regions at the lower and upper surfaces of the coating, from which cracks may ensue. The evolution of such damage in a model transparent system (glass/polycarbonate) is viewed in situ from below and from the side of the specimen. The critical load needed to initiate a crack on the lower coating surface generally increase proportionally to the coatings thickness, d. An interesting departure from this trend occurs for thin coatings, where the fracture load, although marred by a large scatter, increases somewhat with decreasing d. The fracture data for the upper coating surface are limited to relatively thick coatings due to the recurrence of premature failure from the coating edges. The behavior in this range is similar to that for the lower surface crack, albeit with an order of magnitude greater fracture resistance.A fracture mechanics analysis in conjunction with FEM is performed to elucidate the stress intensity factors responsible for crack propagation. A crack normal to the coating surface is assumed to emanate either from the lower or upper surface of the coating. A major feature of the solution is the occurrence of a bending-induced compression stress field over a region ahead of the crack tip. This effect, which become more dominant as the ratio between the contact length and the coating thickness is increased, tends to delay the onset of crack propagation, especially for the lower surface crack. Consequently, in applications associated with large indenters, thin and/or tough coatings and stiff substrates, cracking from the upper coating surface may precede that from the lower surface. An interesting feature of this crack shielding mechanism is that when the coating surface contains a distribution of flaws rather than a single crack, small flaws in this population may be more detrimental than large ones. Incorporation of these aspects into the analysis leads to a good correlation with the test results. In the special case of line loading, which constitutes a lower bound for the critical loads, a closed-form, approximate solution for the stress intensity factors or the critical loads are obtained.Plane-strain indentation, although less common than spherical indentation, allows for characterizing the fracture resistance of opaque films through observation from the specimen edge. This approach is not easily implemented to thin films (i.e., less than about a hundred microns), however.     

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.     

功能梯度材料涂层平面裂纹分析

研究粘接于均质基底材料上功能梯度涂层平面裂纹问题. 假设功能梯度材料剪切模量的倒数为坐标的线性函数,而泊松比为常数. 采用Fourier变换和传递矩阵法将该混合边值问题化为奇异积分方程组,通过数值求解获得 应力强度因子. 考察了材料梯度变化形式、结构几何尺寸和材料梯度参数对裂纹应力强度因子的影响,发现 功能梯度材料涂层尺寸、裂纹长度以及材料梯度参数均对应力强度因子有显著影响.     

Out-of-plane interface cracks in dissimilar piezoelectric materials

Summary This paper investigates the problem of an anti-plane interfacial crack between two dissimilar piezoelectric material layers. A single crack is first considered. The effect of interaction of two collinear cracks in the medium on the field intensity factors is investigated. The solutions of several particular cases, including an infinite piezoelectric bi-material and a piezoelectric material bonded to an elastic medium, are given. The bi-material constants governing the behavior of the crack tip fields are identified. By considering the crack as a notch of finite thickness, it is shown that the thickness of the notch has a pronounced influence on the crack tip field. The results for the assumption of a permeable crack represent the limit case where the notch thickness is reduced to zero.BLW would like to thank the National Science Foundation of China (#10102004) and the City University of Hong Kong (DAG #7100219) for the support of this work. YGS also thanks the Multidiscipline Scientific Research Foundation Project (HIT. MD 2001. 39) of the Harbin Institute of Technology and the SRF for ROCS, SEM.accepted for publication 3 April 2003