Solid damping in micro electro mechanical systems

This paper focuses on the problem of the numerical evaluation of dissipation induced by thermoelastic coupling in microelectromechanical systems. An ad hoc conceived, FE based, numerical procedure for the evaluation of the thermoelastic dissipation is proposed and the numerical results are compared with analytical solutions. In order to introduce in the numerical response a dependence on the size of the resonating devices, which is experimentally observed at very small dimensions, a new enhanced non-local coupled thermoelastic model is proposed and the first results are discussed. An erratum to this article can be found at     

三维热弹性应力波传播问题的数值分析

为求解三维动态热弹性问题提出了一种基于特征线法的数值方法，并对矩形截面的热弹性棒在冲击载荷和热冲击作用下的三维应力波传播过程进行了数值模拟。     

Numerical and analytical study of the propagation of thermoelastic waves in a medium with heat-flux relaxation

The thermoelastic problem of laser exposure of metals and dielectrics is studied taking into account the finite speed of propagation of thermal waves and using a numerical finite-difference algorithm. The resulting numerical solution is compared with the analytical one. The problem is solved in coupled and uncoupled formulations. The solutions of the hyperbolic thermoelastic problem are compared with the solutions of the classical problem. Analytical expressions are obtained for the propagation speeds of the thermoelastic wave components. Times are determined at which the difference between the solutions of the hyperbolic and classical thermoelastic problems can be detected experimentally.     

Properties degradation induced by transverse cracks in general symmetric laminates

This paper presents the details of a methodology for predicting the thermoelastic properties degradation in general symmetric laminates with uniform ply cracks in some or all of the 90° layers. First, a stress transfer method is derived by using the concept of state space equation. The laminate can be subjected to any combination of in-plane biaxial and shear loading, and the uniform thermal loading is also taken into account. The method takes into account all independent material constants and guarantees continuous fields of all interlaminar stresses across interfaces between material layers. By this method, a laminate may be composed of an arbitrary number of monoclinic layers and each layer may have different material property and thickness. Second, the concept of the effective thermoelastic properties of a cracked laminate is introduced. Based on the numerical solutions of specially designed loading cases, the effective thermoelastic constants of a cracked laminate can be obtained. Finally, the applications of the methodology are shown by numerical examples and compared with numerical results from other models and experiment data in the literature. It is found that the theory provides good predictions of the thermoelastic properties degradation in general symmetric laminates.     

Analogue solution of two-dimensional thermoelastic problems with internal heat generation and heat transfer

The analogue procedure which the authors have proposed² is expanded to plane thermoelastic fields, with heat generation within the region or heat transfer from both surfaces of the thermoelastic field. A plate similar to the thermoelastic field is prepared and the plate-bending test is carried out with distributed loads over the plate and with a concentrated load on the inner flange. The strains are measured carefully for each bending test and these are combined to transform into the thermal stress of the original thermoelastic field. The results are compared with theoretical or numerical ones and good agreement is confirmed among them.     

Thermoelastic damping in micro-beam resonators

Thermoelastic damping is recognized as a significant loss mechanism at room temperature in micro-scale beam resonators. In this paper, the governing equations of coupled thermoelastic problems are established based on the generalized thermoelastic theory with one relaxation time. The thermoelastic damping of micro-beam resonators is analyzed by using both the finite sine Fourier transformation method combined with Laplace transformation and the normal mode analysis. The vibration responses of deflection and thermal moment are obtained for the micro-beams with simply supported and isothermal boundary conditions. The vibration frequency is analyzed for three boundary condition cases, i.e., the clamped and isothermal, the simply supported and isothermal, and the simply supported and adiabatic. The analytic results show that the amplitude of deflection and thermal moment are attenuated and the vibration frequency is increased with thermoelastic coupling effect being considered. In addition, it can be found from both the analytic results and the numerical calculations that these properties are size-dependent. When the thickness of the micro-beam is larger than its characteristic size, the effect of thermoelastic damping weakens as the beam thickness increases. The size-effect induced by thermoelastic coupling would disappear when the thickness of the micro-beam is over a critical value that depends on the material properties and the boundary conditions.     

ANLYTICAL SOLUTION OF THERMOELASTIC DAMPING IN RECTANGULAR MINDLIN MICRO PLATES 1)

首次给出了四边简支的 Mindlin 矩形微板热弹性阻尼的解析解. 基于考虑一阶剪切变形的 Mindlin 板理论和单向耦合热传导理论建立了微板热弹性耦合自由振动控制微分方程. 忽略温度梯度在面内的变化,在上下表面绝热边界条件下求得了用变形几何量表示的温度场的解析解. 进一步将包含热弯曲内力的结构振动方程转化为只包含挠度振幅的四阶偏微分方程. 利用特征值问题之间在数学上的相似性,在四边简支条件下给出了用无阻尼 Kirchhoff 微板的固有频率表示的 Mindlin 矩形微板的复频率解析解,从而利用复频率法求得了反映热弹性阻尼水平的逆品质因子. 最后,通过数值结果定量地分析了剪切变形、材料以及几何参数对热弹性阻尼的影响 规律. 结果表明,Mindlin 板理论预测的热弹性阻尼小于 Kirchhoff 板理论预测的热弹性阻尼. 两种理论预测的热弹性阻尼之间的差值在临界厚度附近十分显著. 另外,随着微板的边/厚比增大,Mindlin 微板的热弹性阻尼最大值单调增大,而 Kirchhoff 微板的热弹性阻尼最大值却保持不变.     

Thermoelastic surface waves on an exponentially graded half-space

In this paper, we consider the propagation of Rayleigh surface waves in a functionally graded isotropic thermoelastic half-space, in which all thermoelastic characteristic parameters exponentially change along the depth direction. The propagation condition is established in the form of a bicubic equation whose coefficients are complex numbers while the analytical solutions (eigensolutions) of the thermoelastodynamic system are explicitly obtained in terms of the characteristic solutions. The concerned solution of the Rayleigh surface wave problem is subsequently expressed as a linear combination of the three eigensolutions while the secular equation is established in an implicit form. The explicit secular equation is written when an isotropic and homogeneous thermoelastic half-space is considered and some numerical simulations are given for a specific material.     

Mindlin 矩形微板的热弹性阻尼解析解

首次给出了四边简支的 Mindlin 矩形微板热弹性阻尼的解析解. 基于考虑一阶剪切变形的 Mindlin 板理论和单向耦合热传导理论建立了微板热弹性耦合自由振动控制微分方程. 忽略温度梯度在面内的变化,在上下表面绝热边界条件下求得了用变形几何量表示的温度场的解析解. 进一步将包含热弯曲内力的结构振动方程转化为只包含挠度振幅的四阶偏微分方程. 利用特征值问题之间在数学上的相似性,在四边简支条件下给出了用无阻尼 Kirchhoff 微板的固有频率表示的 Mindlin 矩形微板的复频率解析解,从而利用复频率法求得了反映热弹性阻尼水平的逆品质因子. 最后,通过数值结果定量地分析了剪切变形、材料以及几何参数对热弹性阻尼的影响 规律. 结果表明,Mindlin 板理论预测的热弹性阻尼小于 Kirchhoff 板理论预测的热弹性阻尼. 两种理论预测的热弹性阻尼之间的差值在临界厚度附近十分显著. 另外,随着微板的边/厚比增大,Mindlin 微板的热弹性阻尼最大值单调增大,而 Kirchhoff 微板的热弹性阻尼最大值却保持不变.      

A problem for an infinite elastic body with a spherical cavity in the theory of generalized thermoelastic diffusion

In this work, we study a one-dimensional problem in a generalized thermoelastic diffusion in infinite medium with a spherical cavity subjected to a time dependent thermal shock of its internal boundary which is assumed to be traction free. The chemical potential is also assumed to be a known function of time on the bounding cavity. Laplace transform techniques are used. The solution of the problem in the transformed domain is obtained by using a direct approach without the customary use of potential functions. By means of numerical Laplace inversion, the problem is solved in the physical domain. Numerical results predict finite speeds of propagation for thermoelastic and diffusive waves. To investigate the diffusions effects, a comparison is made with the results obtained in the thermoelastic problem.     

Finite amplitude one-dimensional wave propagation in a thermoelastic half-space

The problem of finite wave propagation in a nonlinearly thermoelastic half-space is considered. The surface of the half-space is subjected to a time-dependent thermal and normal mechanical loading. The solution is obtained by a numerical procedure, which is shown to furnish accurate results, and linear dynamic thermoelastic problems are obtained as special cases. The accuracy of the results is checked by comparison with some known analytical solutions which can be obtained in some special cases of both the linear and the nonlinear problems. In those cases where the solution contains shocks, it is shown that the numerical results satisfy the necessary jumps conditions which need to hold across such discontinuities.     

A modified Green–Lindsay thermoelasticity with strain rate to eliminate the discontinuity

Probing the mechanism of ultrafast thermoelastic processes is becoming increasingly important with the development of laser-assisted micro/nano machining. Although thermoelastic models containing temperature rate have been historically proposed, the strain rate has not been considered yet. In this work, a generalized thermoelastic model is theoretically established by introducing the strain rate in Green–Lindsay (GL) thermoelastic model with the aid of extended thermodynamics. Numerically, a semi-infinite one-dimensional problem is considered with traction free at one end and subjected to a temperature rise. The problem is solved using the Laplace transform method, and the transient responses, i.e. displacement, temperature and stresses are graphically depicted. Interestingly, it is found that the strain rate may eliminate the discontinuity of the displacement at the elastic and thermal wave front. Also, the present model is compared with Green–Naghdi (GN) models. It is found that the thermal wave speed of the present model is faster than GN model without energy dissipation, and slower than GN model with energy dissipation. In addition, the thermoelastic responses from the present model are the largest. The present model based upon GL model is free of the jump of GL model in the displacement distribution, and is safer in engineering practices than GN model. The present work will benefit the theoretical modeling and numerical prediction of thermoelastic process, especially for those under extreme fast heating.     

A numerical method based on boundary integral equations and radial basis functions for plane anisotropic thermoelastostatic equations with general variable coefficients

A boundary integral method with radial basis function approximation is proposed for numerically solving an important class of boundary value problems governed by a system of thermoelastostatic equations with variable coefficients. The equations describe the thermoelastic behaviors of nonhomogeneous anisotropic materials with properties that vary smoothly from point to point in space. No restriction is imposed on the spatial variations of the thermoelastic coefficients as long as all the requirements of the laws of physics are satisfied. To check the validity and accuracy of the proposed numerical method, some specific test problems with known solutions are solved.     

压电结构的热弹性比拟建模方法

通过比较逆压电效应本构方程和热弹性本构方程,在线性范围内建立逆压电效应和热弹性效应的比拟关系:将驱动电压比拟为温度载荷,压电应变系数比拟为热膨胀系数,用热弹性有限元方法分析逆压电效应,并采用热弹性比拟方法求解基于长度伸缩逆压电效应和厚度剪切逆压电效应的驱动器问题,给出了平面结构“热-机-电”耦合问题的热弹性比拟求解方法。热弹性比拟方法为快速建立复杂压电结构的有限元模型提供了一种有效的途径,可缩短复杂压电主动结构的设计周期,降低设计成本。最后给出数值算例说明文中方法的有效性。     

THREE-DIMENSIONAL THERMOELASTIC ANALYSIS OF FUNCTIONALLY GRADED PLATE

Three-dimensional thermoelastic analysis is presented for an orthotropic functionally graded rectangular plate,which is simply supported and isothermal on its four lateral edges.With the assumption that material properties have arbitrary dependence on the thickness-coordinate,a Peano-Baker series solution is obtained for the thermoelastic fields of the functionally graded plate subjected to mechanical and thermal loads on its upper and lower surfaces by means of state space method.The correctness of the obtained series solution is validated through numerical examples.The influence of different material properties distributions on the structural response of the plate is also studied.     

A numerical method based on boundary integral equations and radial basis functions for plane anisotropic thermoelastostatic equations with general variable coe?cients

A boundary integral method with radial basis function approximation is proposed for numerically solving an important class of boundary value problems governed by a system of thermoelastostatic equations with variable coe?cients. The equations describe the thermoelastic behaviors of nonhomogeneous anisotropic materials with properties that vary smoothly from point to point in space. No restriction is imposed on the spatial variations of the thermoelastic coe?cients as long as all the requirements of the laws of physics are satis?ed. To check the validity and accuracy of the proposed numerical method, some speci?c test problems with known solutions are solved.     

Analytical Solution for a Free Vibration of a Thermoelastic Hollow Sphere

For a free vibration problem of a thermoelastic hollow sphere into the context of the generalized thermoelasticity theory with one relaxation time, exact analytic solutions are obtained with the use of eigenvalue approach. Both the inner and outer curved surfaces of the sphere are considered stress-free and isothermal surfaces. The dispersion relations for the existence of various types of possible modes of vibrations in the considered hollow sphere are derived. The numerical results have been presented graphically in respect of natural frequencies, thermoelastic damping, and frequency shift.     

Two-dimensional thermoelastic analysis of a functionally graded cylindrical panel due to nonuniform heat supply

This paper is concerned with the theoretical treatment of the steady-state thermoelastic problem of a functionally graded cylindrical panel due to nonuniform heat supply in the circumferential direction. The thermal and thermoelastic constants of the cylindrical panel are expressed as power functions of the radial coordinate. We obtain the exact solution for the two-dimensional temperature change in a steady state, and thermal stresses of a simple supported cylindrical panel under the state of plane strain. Some numerical results are shown in figures and tables. Furthermore, the influence of the nonhomogeneity of the material, the radius ratio and the span angle upon the temperature change, displacements and stresses is investigated.     

Coupled thermoelasticity of functionally graded cylindrical shells

The coupled thermoelastic response of a functionally graded circular cylindrical shell is studied. The coupled thermoelastic and the energy equations are simultaneously solved for a functionally graded axisymmetric cylindrical shell subjected to thermal shock load. A second-order shear deformation shell theory that accounts for the transverse shear strains and rotations is considered. Including the thermo-mechanical coupling and rotary inertia, a Galerkin finite element formulation in space domain and the Laplace transform in time domain are used to formulate the problem. The inverse Laplace transform is obtained using a numerical algorithm. The shell is graded through the thickness assuming a volume fraction of metal and ceramic, using a power law distribution. The results are validated with the known data in the literature.     

Two-dimensional interactions due to moving load in generalized thermoelastic solid with diffusion

The present paper is concerned with the investigation of disturbances in＇a homogeneous, isotropic elastic medium with generalized thermoelastic diffusion, when a moving source is acting along one of the co-ordinate axis on the boundary of the medium. Eigen value approach is applied to study the disturbance in Laplace-Fourier transform domain for a two dimensional problem. The analytical expressions for displacement components, stresses, temperature field, concentration and chemical potential are obtained in the physical domain by using a numerical technique for the inversion of Laplace transform based on Fourier expansion techniques. These expressions are calculated numerically for a copper like material and depicted graphically. As special cases, the results in generalized thermoelastic and elastic media are obtained. Effect of presence of diffusion is analyzed theoretically and numerically.