Thermoelastic analysis for an opening crack in an orthotropic material

The thermoelastic analysis of an opening crack embedded in an orthotropic material is made under applied uniform heat flux and mechanical loadings. To simulate the case of an opening crack filled with a medium, a thermal-medium crack model is proposed. The thermally permeable and impermeable cracks are the limiting ones of the proposed thermal-medium one. The crack-tip thermoelastic fields induced by a crack in an orthotropic material are determined in closed forms. The elastic T-stress can be also obtained explicitly. The effects of applied mechanical loadings and the thermal conductivity of crack interior on the heat flux at the crack surfaces and the mode-II stress intensity factor are investigated through numerical computations. The obtained results reveal that an increase of the thermal conductivity of crack interior decreases the mode-II stress intensity factor. And when an applied mechanical loading is increasing, the mode-II stress intensity factor is rising.     

Scattering of the harmonic anti-plane shear waves by a crack in functionally graded piezoelectric materials

In the present paper, the dynamic behavior of a Griffth crack in the functionally graded piezoelectric material (FGPM) is investigated. It is assumed that the elastic stiffness, piezoelectric constant, dielectric permittivity and mass density of the FGPM vary continuously as an exponential function, and that FGPM is under the anti-plane mechanical loading and in-plane electrical loading. By using the Fourier transform and defining the jumps of displacement and electric potential components across the crack surface as the unknown functions, two pairs of dual integral equations are derived. To solve the dual integral equations, the jumps of the displacement and electric potential components across the crack surface are expanded in a series of Jacobi polynomial. Numerical examples are provided to show the effects of material properties on the stress and the electric displacement intensity factors.     

在反平面变形下功能梯度压电板条的应力强度因子分析

利用Fourier变换技术将混合边值问题化为对偶积分方程,求解对偶积分方程得到应力强度因子以及电位移强度因子的表达式。最后通过数值计算讨论了材料参数、载荷条件以及裂纹的几何参数等对功能梯度压电材料中裂纹尖端应力强度因子的影响。     

Non-local theory solution of a mode-I crack in a piezoelectric/piezomagnetic composite material plane

The non-local theory solution of a mode-I permeable crack in a piezoelectric/piezomagnetic composite material plane was given by using the generalized Almansi’s theorem and the Schmidt method in this paper. The problem was formulated through Fourier transform into two pairs of dual integral equations, in which the unknown variables are the displacement jumps across the crack surfaces. To solve the dual integral equations, the displacement jumps across the crack surfaces were directly expanded as a series of Jacobi polynomials. Numerical examples were provided to show the effects of the crack length and the lattice parameter on the stress field, the electric displacement field and the magnetic flux field near the crack tips. Unlike the classical elasticity solutions, it is found that no stress, electric displacement and magnetic flux singularities are present at the crack tips in piezoelectric/piezomagnetic composite materials. The non-local elastic solution yields a finite hoop stress at the crack tip, thus allowing us to use the maximum stress as a fracture criterion.     

The extended finite element method in thermoelastic fracture mechanics

The extended finite element method (XFEM) is applied to the simulation of thermally stressed, cracked solids. Both thermal and mechanical fields are enriched in the XFEM way in order to represent discontinuous temperature, heat flux, displacement, and traction across the crack surface, as well as singular heat flux and stress at the crack front. Consequently, the cracked thermomechanical problem may be solved on a mesh that is independent of the crack. Either adiabatic or isothermal condition is considered on the crack surface. In the second case, the temperature field is enriched such that it is continuous across the crack but with a discontinuous derivative and the temperature is enforced to the prescribed value by a penalty method. The stress intensity factors are extracted from the XFEM solution by an interaction integral in domain form with no crack face integration. The method is illustrated on several numerical examples (including a curvilinear crack, a propagating crack, and a three‐dimensional crack) and is compared with existing solutions. Copyright © 2007 John Wiley & Sons, Ltd.     

Micro-crack initiation at the tip of a semi-infinite rigid line inhomogeneity in piezoelectric solids

A dislocation emission mechanism for micro-crack initiation at the tip of a semi-infinite rigid line inhomogeneity in a piezoelectric solid is proposed in the present paper. For a rigid line inhomogeneity embedded in a piezoelectric matrix, dislocations of one sign are driven away from the tip due to high stress level, while the stationary dislocations of the opposite sign are left behind near the tip of the inhomogeneity. As a result, a micro-Zener–Stroh crack is initiated ahead of the line inhomogeneity. In the current study, a dislocation pileup mechanism for micro-crack initiation at the inhomogeneity tip is proposed. An interesting result is that the critical stress intensity factors for a line inhomogeneity perpendicular to the poling direction can be related to the fracture toughness of a conventional crack in the same material. Analytical solutions show that the critical plane shear stress intensity factor depends on the plane shear mechanical and displacement loadings, and the critical opening stress and electric displacement intensity factors depend on not only the mechanical and displacement loadings, but the electric field and displacement loadings as well.     

Hybrid finite element analysis of transient thermoelastic fracture problems subjected to general heat transfer conditions

This paper presents the singular characteristics of heat flux in the vicinity of the crack-tip for two dimensional transient thermoelastic fracture problems subjected to general heat transfer conditions at crack surfaces. Based on a restricted variational principle, a rigorous hybrid finite element procedure is then developed to perfectly describe the singularities of heat flux and thermal stress induced at the crack-tip. For verification purposes, the examples of transient thermoelastic problems with insulated crack surfaces are first analyzed. Excellent agreements between the computed results and referenced solutions can be drawn. To evaluate the influence of heat convection and radiation on the computation of temperature distributions and thermal stress intensity factors, several numerical examples are also presented.     

The penny-shaped crack problem in micropolar thermoelasticity

The axisymmetry problem of a penny-shaped crack opened out by thermal loads is studied. The linear theory of micropolar elasticity is employed. Two types of thermal loads are considered—prescribed temperature on the crack faces and prescribed heat flux across the faces. It is shown that, in both the cases, the problem is equivalent to the isothermal problem of the crack opened out by suitable normal tractions on the crack faces. The stress intensity factors are found to depend on, in addition to the usual parameters, two parameters N and M; N is a number characterising the coupling of the displacement field with the microtation field and M is the ratio N/τ where τ is a non-dimensional material characteristic length. The classical values of the stress intensity factors are recovered as a limiting case. Numerical results are presented for the case of constant heat flux across the crack faces. These results show that the presence of couple stresses elevates the values of the stress intensity factors.     

A magneto-electro-elastic material with a penny-shaped crack subjected to temperature loading

Summary The analysis of intensity factors for a penny-shaped crack under thermal, mechanical, electrical and magnetic boundary conditions becomes a very important topic in fracture mechanics. An exact solution is derived for the problem of a penny-shaped crack in a magneto-electro-thermo-elastic material in a temperature field. The problem is analyzed within the framework of the theory of linear magneto-electro-thermo-elasticity. The coupling features of transversely isotropic magneto-electro-thermo-elastic solids are governed by a system of partial differential equations with respect to the elastic displacements, the electric potential, the magnetic potential and the temperature field. The heat conduction equation and equilibrium equations for an infinite magneto-electro-thermo-elastic media are solved by means of the Hankel integral transform. The mathematical formulations for the crack conditions are derived as a set of dual integral equations, which, in turn, are reduced to Abel's integral equation. Solution of Abel's integral equation is applied to derive the elastic, electric and magnetic fields as well as field intensity factors. The intensity factors of thermal stress, electric displacement and magnetic induction are derived explicitly for approximate (impermeable or permeable) and exact (a notch of finite thickness crack) conditions. Due to its explicitness, the solution is remarkable and should be of great interest in the magneto-electro-thermo-elastic material analysis and design.     

Thermal stresses near a penny-shaped crack in an elastic sphere embedded in an infinite elastic space

This paper gives an analysis of the distribution of thermal stresses in a sphere which is bonded to an infinite elastic medium. The thermal and the elastic properties of the sphere and the elastic infinite medium are assumed to be different. The penny-shaped crack lies on the diametral plane of the sphere and the centre of the crack is the centre of the sphere. By making a suitable representation of the temperature function, the heat conduction problem is reduced to the solution of a Fredholm integral equation of the second kind. Using suitable solution of the thermoelastic displacement differential equation, the problem is then reduced to the solution of a Fredholm integral equation, in which the solution of the earlier integral equation arising from heat conduction problem occurs as a known function. Numerical solutions of these two Fredholm integral equations are obtained. These solutions are used to evaluate numerical values for the stress intensity factors. These values are displayed graphically.     

Stress intensity factors for curvilinear cracks loaded under anti-plane strain (mode III) conditions

A general solution to the elastic and thermoelastic problems with a rigid circular-arc inclusion is presented. The proposed analysis is based upon the complex variable theory dealing with sectionally holomorphic functions which is reduced to the solution of the Hilbert problem. It is indicated that both the stress and thermal stress fields near the inclusion tip possess a square-root singularity similar to that for the corresponding crack problem. In analogy to the stress intensity factors defined for crack problem, stress singularity coefficients are introduced in this paper to characterize the near tip fields. Complete stress fields and the corresponding stress singularity coefficients as the circular-arc inclusion are under uniform remote load, concentrated force and uniform heat flux are given. Failure initiation of an infinite plate embedded with a rigid arc inclusion under different loading conditions is also discussed.     

Transient response of a center crack in a functionally graded piezoelectric strip under electromechanical impact

The dynamic fracture problem for a functionally graded piezoelectric strip containing a center crack parallel to the free boundaries is considered in this study. It is assumed that the electroelastic properties of the medium vary continuously in the thickness direction, and that the strip is under in-plane mechanical and electric impact loadings. Integral transform techniques and dislocation density functions are employed to reduce the problem to the solutions of a system of singular integral equations. The dynamic stress and electric displacement intensity factors versus time are presented for various values of dimensionless parameters representing the crack size, the material nonhomogeneity and the loading combination.     

An Efficiency Study of Transient Wave Generation in a Thermoelastic Layer with a Pulsed Laser

The efficiency of transient wave generation in a thermoelastic silicon layer excited by a pulsed laser is considered. First a principle-based transfer matrix formulation with relaxation effect, also referred to as the generalized dynamic theory of linear thermoelasticity, is used in obtaining transfer functions between the input heat field and the elements of the thermoelastic state vector. The second sound effect, through this relaxation time term, is included to eliminate the thermal wave travelling with infinite velocity as predicted by the diffusion heat transfer model. By employing the fast Fourier transform (FFT) algorithm, the transient response of a silicon thermoelastic layer under a thermal excitation (by a pulsed laser) is investigated to quantify the conversion efficiency from thermal to mechanical energy. The transient acceleration, stress, heat, temperature, and mechanical power flux responses are presented. The pulse duration of the laser excitation is submicrosecond level and, consequently, a large number of modes of motion are excited. Rigid body singularities are eliminated by considering the higher order time derivatives of the state variables. A layer made of bulk silicon under this laser excitation is considered and it is found that the amplitude ratio of the applied heat field to the propagating heat flux at the data points is in the order of 10°. The ratio of the applied power (heat flux) to the generated mechanical power flux is in the order of 10°. The resulting rigid body motion of the layer due to the laser excitation is excluded in calculating the mechanical power.     

Axisymmetric thermoelastic shape sensitivity analysis and its application to turbine disc design

A general method for shape design sensiti vityi analysis (SDSA) as applied to an axisymmetric thermoelasticity problem is presented using the material derivative concept and the adjoint variable method. The sensitivity of a general functional composed of thermal and mechanical quantities is considered. The method for deriving the sensitivity formula is based on standard direct thermal and elastic boundary integral equation formulation. It is then applied to obtain explicit formulas for a representative displacement and stress constraint imposed on a sector of the boundary. Results of numerical implementation are presented for weight minimization of a turbine disc under thermomechanical loading. The sensitivities of the displacement and stress constraint calculated by the formulas are compared with those by finite differences. Optimum shape obtained under the thermomechanical loading is discussed with that under the mechanical loading only, clearly showing the practical importance of the SDSA of thermoelastic systems.     

非均匀复合材料中反平面裂纹的动态断裂力学研究

对于非均匀复合材料中多个裂纹的动态断裂力学问题, 提出了一种分析方法, 假设复合材料为正交各向异性并含有多个垂直于厚度方向的裂纹, 材料参数沿厚度方向为变化的, 沿该方向将复合材料划分为许多单层, 假设单层材料参数为常数, 应用柔度矩阵/刚度矩阵方法及Fourier变换法, 在L aplace 域内推导出了控制问题的奇异积分方程组, 并用虚位移原理求解, 给出了应力强度因子及能量释放率的表达式, 然后利用Laplace 数值反演, 得出了裂纹尖端的动态应力强度因子和能量释放率。作为算例, 研究了带有两个裂纹的功能梯度结构, 分析了材料参数的优化对降低应力强度因子的意义。      

Thermal stress singularities at tips of a crack in a semi-infinite medium under uniform heat flow

In the framework of plane thermoelastic problems is discussed the thermal stress field near the tips of an arbitrarily inclined crack in an isotropic semi-infinite medium with the thermally insulated edge surface under uniform heat flow. The crack is replaced by continuous distributions of quasi-Volterra dislocations corresponding to line heat sources and edge dislocations, and we obtain a set of simultaneous singular integral equations for dislocation density functions, whose solution is given in the forms of series in terms of Tchebycheff polynomials of the first kind. By means of this method, the thermal stress singularities at the crack tips are estimated exactly and the stress intensity factors can be readily evaluated. Numerical results are given for the particular case where the surface of the inclined crack is maintained at constant temperature and the heat supplied across the surface of the crack vanishes as a whole. The effects of the distance from the crack tip to the edge surface of the semi-infinite medium and the angle of inclination of the crack on the stress intensity factors and the initial direction of crack extension are shown graphically.     

Dynamic stress intensity factors of two 3D rectangular permeable cracks in a transversely isotropic piezoelectric material under a time‐harmonic elastic P‐wave

The dynamic behavior of two 3D rectangular permeable cracks in a transversely isotropic piezoelectric material is investigated under an incident harmonic stress wave by using the generalized Almansi's theorem and the Schmidt method. The problem is formulated through double Fourier transform into three pairs of dual integral equations with the displacement jumps across the crack surfaces as the unknown variables. To solve the dual integral equations, the displacement jumps across the crack surfaces are directly expanded as a series of Jacobi polynomials. Finally, the relations among the dynamic stress field and the dynamic electric displacement filed near the crack edges are obtained, and the effects of the shape of the rectangular crack, the characteristics of the harmonic wave, and the distance between two rectangular cracks on the stress and the electric intensity factors in a piezoelectric composite material are analyzed. Copyright © 2015 John Wiley & Sons, Ltd.     

Dual boundary element method for three-dimensional thermoelastic crack problems

This paper describes the formulation and numerical implementation of the three-dimensional dual boundary element method (DBEM) for the thermoelastic analysis of mixed-mode crack problems in linear elastic fracture mechanics. The DBEM incorporates two pairs of independent boundary integral equations; namely the temperature and displacement, and the flux and traction equations. In this technique, one pair is applied on one of the crack faces and the other pair on the opposite one. On non-crack boundaries, the temperature and displacement equations are applied. This revised version was published online in July 2006 with corrections to the Cover Date.     

Dynamic crack extension along the interface of materials that differ in their thermal properties: with and without thermal relaxation

Summary Moving surface stresses cause crack extension along the interface of perfectly bonded thermoelastic materials at a constant sub-critical speed. The materials differ only in their thermal properties, and are governed by coupled thermoelastic equations that admit as special cases Fourier heat conduction as well as thermal relaxation with one or two relaxation times. A dynamic steady state of plane strain is assumed. The exact transform solution for a propagating displacement and temperature discontinuity is used to find solutions to the interface crack valid away from the crack edge for low extension speeds and solutions valid at the crack edge for high speeds. Results show that Fourier heat conduction dominates the former case, but solution behavior in the latter is dependent upon the particular thermal model. Thermal mismatch is seen to by itself cause a solution behavior similar to that for bonded dissimilar isothermal elastic solids. In particular, the two-relaxation time solution exhibits both oscillatory and non-oscillatory terms, and the interface temperature at the crack edge is finite.