刚性圆形压头作用下热电材料半平面的无摩擦接触问题

热电材料可以将热能转化为电能,反之亦然,这一优良的性质将有助于研发更具成本效益的设备和器件。本文研究了刚性圆形压头作用在热电材料半平面的无摩擦接触问题。假定压头为电导体、热导体,且压头压入深度及与材料的接触区域宽度未知。首先求解电场和温度场,利用傅里叶变换得到了电势函数、温度、电流密度和能量通量的解析表达式。然后求解弹性场,利用积分变换和边界条件,将该热弹性接触问题转化为第一类奇异积分方程并数值求解。数值结果讨论了压头半径和热电载荷对法向接触应力、电流强度因子和能量通量强度因子的影响。结果表明,对于圆压头,热电材料的法向电流密度、法向能量通量在接触边缘表现出奇异性,而表面法向接触应力在接触边缘为零。本文建立的研究模型有助于更深层次的了解热电材料的接触行为。     

Thermo-electro-mechanical contact behavior of a finite piezoelectric layer under a sliding punch with frictional heat generation

The thermal contact problem of a piezoelectric strip with heat supply generated by the frictional tangential traction under the action of a rigid sliding punch is investigated. The inertial effects are considered. It is convenient to introduce the Galilean transform. Whole cases of the root distribution of the corresponding characteristic equation are detailed. Appropriate fundamental solutions that can lead to real solutions of the thermo-electro-mechanical quantities are derived for the piezoelectric governing equation. The stated problem is reduced to Cauchy singular integral equation of the second kind finally. Numerical results are also presented. The solutions have a reduced dependence on the material properties. The singular behaviors at the edges of the punch are revealed. The stress distribution and temperature distribution above the punch with the variations of the relative sliding speed, the frictional coefficient and the thickness are plotted. The effects of the material constants on the stress distribution and temperature distribution above the punch are presented.     

A solution for an isotropic sector under anti-plane shear loadings

In this study, the problem of an isotropic sector subjected to anti-plane shear loadings is investigated. The loadings were applied to the arc of the sector, and the radial edges of the sector were under traction-free or fixed conditions. Depending on these conditions, three problems, namely, free–free, fixed–free, and fixed–fixed edges were studied. A procedure using the finite Mellin transform combined with the Laplace transform was proposed for solving these problems. Explicit closed form solutions for the displacement and stress fields throughout the sector were obtained. The stress intensity factor (SIF) for each problem was analyzed using the obtained stress fields. It was determined that the SIF disappeared under the special condition of a fixed–fixed edge. Other special cases having anti-symmetric conditions were deduced from the derived solutions, and the results of these verified those cited in the literature as well as those obtained using finite element analysis (FEA).     

A generalised stress intensity approach to characterising the process zone in complete fretting contacts

The known analytical contact solution for the stress field induced by a rigid, square-ended punch, sliding on an elastic half-plane defines the stress state everywhere in the half-plane. An asymptotic approach is then used to determine the characteristic stress field at the edge of the contact, which is matched with the contact solution. Hence, the regions over which the asymptotic solution is valid are found. Using a method analogous to the crack-tip stress field, a generalised stress intensity factor is defined, with the aim of providing a single variable characterisation of the stress state at the punch corner. The crack initiation process zone for a fretting fatigue crack is therefore captured, and the conditions for small scale yielding explicitly found.     

Thermoelastic transient response of an infinitely long annular multilayered cylinder

Thermoelastic transient response of multilayered annular cylinders of infinite lengths subjected to known temperature at traction-free inner and outer surfaces are considered. A method based on the Laplace transformation and finite difference method has been developed to analyze the thermoelasticity problem. Using the Laplace transform with respect to time, the general solutions of the governing equation are obtained in transform domain. The solution is obtained by using the matrix similarity transformation and inverse Laplace transform. Solutions for the temperature and thermal stress distributions in a transient state were obtained. It was found that the temperature distribution, the displacement and the thermal stresses change slightly as time increases. There is no limit of number of annular layers of the cylinder in the presented computational procedures.     

Electro-mechanical frictionless contact behavior of a functionally graded piezoelectric layered half-plane under a rigid punch

The frictionless contact problem of a functionally graded piezoelectric layered half-plane in-plane strain state under the action of a rigid flat or cylindrical punch is investigated in this paper. It is assumed that the punch is a perfect electrical conductor with a constant potential. The electro-elastic properties of the functionally graded piezoelectric materials (FGPMs) vary exponentially along the thickness direction. The problem is reduced to a pair of coupled Cauchy singular integral equations by using the Fourier integral transform technique and then is numerically solved to determine the contact pressure, surface electric charge distribution, normal stress and electric displacement fields. For a flat punch, the normal stress intensity factor and electric displacement intensity factor are also given to quantitatively characterize the singularity behavior at the punch ends. Numerical results show that both material property gradient of the FGPM layer and punch geometry have a significant influence on the contact performance of the FGPM layered half-plane.     

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

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

Transient response of a piezoelectric ceramic strip with an eccentric crack under electromechanical impacts

The transient response of a piezoelectric strip with an eccentric crack normal to the strip boundaries under applied electromechanical impacts is considered. By using the Laplace transform, the mixed initial-boundary-value problem is reduced to triple series equations, then to a singular integral equation of the first kind by introducing an auxiliary function. The Lobatto–Chebyshev collocation technique is adopted to solve numerically the resulting singular integral equation. Dynamic field intensity factors and energy release rate are obtained for both a permeable crack and an impermeable crack. The effects of the crack position and the material properties on the dynamic stress intensity factor are examined and numerical results are presented graphically.     

Impact of a plane die on an elastic orthotropic half-plane

To study the process of impact of a rigid body on the surface of an elastic body made of a composite material, we consider a nonstationary dynamic contact problem about the impact of a plane rigid die on an elastic orthotropic half-plane. The problem is reduced to solving an integral equation of the first kind for the Laplace transform of the contact stresses under the die base. An approximate solution of the integral equation is constructed with the use of a special approximation to the symbol of the kernel of the integral equation in the complex plane. The inverse Laplace transform of the solution results in determining the scalar contact stress field on the die base, the force exerted by the die on the elastic medium, and the vertical displacement field of the free surface of the orthotropic medium out side the die. The solutions thus obtained permit studying specific features of the process of die penetration into an orthotropic medium and the strain properties of the medium.     

Transient dynamic response of a coated piezoelectric strip with a vertical crack

In this study, the dynamic response of a coated piezoelectric strip containing a crack vertical to the interfaces under normal impact load is considered. Based on the superposition principle and the integral transform techniques, the solution in the Laplace transformed plane is obtained in terms of a singular integral equation. The order of stress singularity around the tip of the terminated crack is also obtained. The singular integral equation is solved by using the Gauss–Jacobi integration formula, and the numerical Laplace inversion is then carried out to obtain the resulting dynamic stress and electric displacement intensities. The effects of the material properties and the geometric parameters on the dynamic stress intensity factor and the dynamic energy density factors are shown graphically.     

Torsional impact of transversely isotropic solid with functionally graded shear moduli and a penny-shaped crack

The torsional impact response of a penny-shaped crack in an unbounded transversely isotropic solid is considered. The shear moduli are assumed to be functionally graded such that the mathematics is tractable. Laplace transform and Hankel transform are used to reduce the problem to solving a Fredholm integral equation. The crack tip stress fields are obtained. Investigated are the influence of material nonhomogeneity and orthotropy on the dynamic stress intensity factor. The peak value of the dynamic stress intensity factor can be suppressed by increasing the shear moduli's gradient and/or increasing the shear modulus in a direction perpendicular to the crack surface.     

Viscoelastic composite with a cracked layer under anti-plane shear

A viscoelastic cracked sandwich composite under anti-plane load is presented via the application of the Laplace transform and the complex variable formulation. Using a special technique of analytical continuation associated with the method of image, the triple-layered problem is reduced to a singular function by simple algebraic manipulations. The superposition method and the singular integral equation is applied to deal with the concentrated load and crack problem. Finally, some typical viscoelastic models and its corresponding stresses and stress intensity factors are also discussed. The results show that the time dependent stress intensity factor is affected by the relative strength of creep compliances of each layer.     

功能梯度材料裂纹尖端动态应力场

研究受反平面剪切作用的功能梯度材料动态裂纹问题,通过积分变换-对偶积分方程方法推出了裂纹尖端动态应力场,时间域内的动态应力强度因子由Laplace数值反演获得,研究结果表明功能梯度材料的梯度越大,相应的裂纹问题的动态应力强度因子值越低。     

The evolution of the process zone when a complete contact is subject to cyclically varying loads

A rigid punch pressed against an incompressible half-plane by a constant normal load is subject to both a cyclically varying shearing force, and a synchronously varying underlying bulk tension. The evolution of the slip/stick distribution along the contact interface is found during a loading cycle. An asymptotic approach is subsequently used to determine the order of singularity, and hence spatial distribution of the local stress state. A generalised stress intensity factor is then defined to scale the corresponding asymptotic solution, and hence map out the evolution of the process zone. The analysis shows that there is a smooth transition in the form of the asymptotic solution in moving from an ‘adhered’ contact edge to a ‘slipping’ contact edge, and vice versa. The implications of the results to fretting fatigue are discussed.     

Thermo-mechanical contact behavior of a finite graded layer under a sliding punch with heat generation

The problem of a rigid punch contacting with a finite graded layer on a rigid substrate is investigated within the framework of steady-state plane strain thermoelasticity, in which heat generated by contact friction is considered with a constant friction coefficient and inertia effects are neglected. The material properties of the graded layer vary according to an exponential function in the thickness direction. Fourier integral transform method and transform matrix approach are employed to reduce the current thermocontact problem to the second kind of Cauchy-type singular integral equation. Distributions of the contact pressure and the in-plane stress under the prescribed thermoelastic environment with different parameter combinations, including ratio of shear moduli, relative sliding speed, friction coefficient and thermal parameters are obtained and analyzed, as well as the stress singularity and the stress intensity factors near the contact edges. The results should be helpful for the design of surfaces with strong wear resistance and novel graded materials for real applications.     

Thermoelastic contact mechanics for a flat punch sliding over a graded coating/substrate system with frictional heat generation

The problem of thermoelastic contact mechanics for the coating/substrate system with functionally graded properties is investigated, where the rigid flat punch is assumed to slide over the surface of the coating involving frictional heat generation. With the coefficient of friction being constant, the inertia effects are neglected and the solution is obtained within the framework of steady-state plane thermoelasticity. The graded material exists as a nonhomogeneous interlayer between dissimilar, homogeneous phases of the coating/substrate system or as a nonhomogeneous coating deposited on the substrate. The material nonhomogeneity is represented by spatially varying thermoelastic moduli expressed in terms of exponential functions. The Fourier integral transform method is employed and the formulation of the current thermoelastic contact problem is reduced to a Cauchy-type singular integral equation of the second kind for the unknown contact pressure. Numerical results include the distributions of the contact pressure and the in-plane component of the surface stress under the prescribed thermoelastic environment for various combinations of geometric, loading, and material parameters of the coated medium. Moreover, in order to quantify and characterize the singular behavior of contact pressure distributions at the edges of the flat punch, the stress intensity factors are defined and evaluated in terms of the solution to the governing integral equation.     

三维动态裂纹问题的超奇异积分方程法

基于弹性材料的动态基本方程,结合广义Betti-Rayleigh互易等式与时域下的边界积分方程,推导得到时域下的超奇异积分方程组。引入Laplace域下的动态基本解,将经过主部分析的积分核函数分解为静态和动态部分,其中动态积分核不具有奇异性。在裂纹前沿附近单元,采用与理论分析一致的平方根位移模型。结合Lubich时间卷积实现拉氏变换,采用配置点法计算超奇异积分,获得问题的数值解。并针对椭圆裂纹算例编写Fortran程序,得到冲击荷载作用下张开型裂纹的动态应力强度因子变化规律,数值结果稳定且收敛速度快。     

Galin’s problem for a spatial elastic wedge

We study a three-dimensional contact problem on the indentation of an elliptic punch into a face of a linearly elastic wedge. The wedge is characterized by two parameters of elasticity and its edge is subjected to the action of an additional concentrated force. The other face wedge is free from stresses. The problem is reduced to an integral equation for the contact pressure. An asymptotic solution of this equation is obtained which is effective for a given contact region fairly remote from the edge. Calculations are performed that allow one to evaluate the effect of a force applied outside the contact region on the contact pressure distribution. The problem under study is a generalization of L. A. Galin’s problem on a force applied outside a circular punch on an elastic half-space [1, 2]. In a special case of a wedge with an opening angle of 180° and zero contact ellipse eccentricity, the obtained asymptotic relation coincides with the expansion of Galin’s exact solution in a series. Problems of indentation of an elliptic punch into a spatial wedge with the face not loaded outside the contact region have been studied previously. For example, the paper [3] dealt with the case of a known contact region (asymptotic method) and the paper [4] considered the case of an unknown contact region (numerical method). The solution of Galin’s problem allowed the authors of [2] to reduce the contact problem on the interaction of several punches applied to a half-space to a system of Fredholm integral equations of the second kind (Andreikin-Panasyuk method). A topical direction in contact mechanics is the model of discrete contact as well as related problems on the interaction of several punches [2, 5–8]. The interaction of several punches applied to a face of a wedge can be treated in a similar manner and an asymptotic solution can be obtained for the case where a concentrated force is applied at an arbitrary point of this face beyond the contact region rather than on the edge.     

Impact response of a piezoelectric layered composite plate with a crack

The dynamic response of a central crack in a piezoelectric layered composite plate under normal impact is analyzed. The crack is oriented normally to the interfaces. The Laplace and Fourier transform techniques are used to formulate the problem in terms of a singular integral equation. The order of stress singularity around the tip of the terminated crack is also obtained. Numerical calculations are carried out, and the main results presented are the variations of the dynamic stress intensity factor and the dynamic energy density factor versus time as functions of the geometric parameters and the piezoelectric material properties of the layered composite plate.     

沿厚度非均匀复合材料的动态断裂力学研究

对于非均匀复合材料中多个裂纹的动态断裂力学问题，提出了一种分析方法，假设复合材料为正交各向异性并含有多个垂直于厚度方向的裂纹，材料参数沿厚度方向为变化的，沿该方向将复合划分为许多单层，假设单层材料参数为常数，Ｆｏｕｒｉｅｒ变换法，在Ｌａｐｌａｃｅ域内推导出了控制问题的奇异积分方程组并用虚位移原理求解，然后利用Ｌａｐｌａｃｅ数值反得刺裂纹尖端的动态应力强度因子和能量释放率，作为算例，研究了带有两个裂