激光脉冲加热的二维广义热弹性问题研究

基于L-S广义热弹性理论,研究了半无限大板局部受到激光脉冲加热时的广义热弹性问题.为避免常规积分变换方法求解带来的精度丢失,采用有限元法直接在时间域进行求解,得到了激光脉冲加热时板中的温度、位移及应力的变化规律.结果表明,直接求解方法可以准确描述热在介质中以有限的速度传播,同时发现,激光脉冲加热过后,结构的最高温度随着时间的推移逐渐降低,且最高温度的位置总在热波波前附近,此处的应力也明显高于其他区域.     

直接有限元法求解广义磁热弹二维旋转问题

为了验证直接有限元法求解广义磁热弹耦合旋转问题的有效性及准确性,该文基于Lord和Shulman(L-S)广义热弹性理论,采用直接有限元方法,求解了置于磁场中的旋转半无限大体受热冲击作用的动态响应问题.文中给出了L-S型广义磁热弹耦合旋转问题的控制方程,建立了L-S型广义磁热弹旋转问题的虚位移原理,推导得到了相应的有限...     

属性与温度相关材料的广义热弹性问题有限元分析

计及材料特性与温度的相关性,基于Lord和Shulman(L-S)广义热弹性理论,建立了此类问题的有限元控制方程. 由于材料属性的温度相关性,温度控制方程具有非线性,积分变换求解方法难以采用,因而将有限元方程直接在时间域求解. 利用所建立方法研究了材料特性与温度相关、带有孔洞的无限大体在热冲击和机械冲击作用下的广义热弹性问题. 分析表明,在时间域直接求解材料属性与温度相关的广义热弹性问题是可行的,所得结果具有很高的精度,热的波动性得到充分的展现. 同时发现,热冲击载荷作用时,材料属性与温度的相关性对结构的机械响应影响显著,对温度响应影响很小;机械载荷作用时,材料参数与温度的相关性对所有响应影响都很小. 因此,研究热冲击载荷作用的机械响应时,必须考虑材料属性的温度相关性,而研究温度响应时,无论何种冲击载荷,都可以不考虑材料属性的温度相关性.      

无限长圆柱体中广义电磁热弹耦合问题研究

基于L-S广义热弹性理论,研究了处于磁场中无限长理想圆柱导体在边界受热冲击作用时的电磁热弹耦合问题.建立了广义电磁热弹耦合的有限元方程,为避免积分变换方法求解带来的精度丟失.采用将有限元方程直接在时间域求解的方法,得到了圆柱体中的温度、位移、应力、感应磁场和感应电场的分布规律,反映了热的波动性及电磁热弹的耦合效应.结果表明,将有限元方程直接在时间域求解,可以获得各物理量的准确分布.得到温度在热波波前处的阶跃,准确地反应热波的波动效应.     

半无限长压电杆的瞬态热冲击问题

采用具有一个热松弛时间的L－S广义热弹性理论，研究了一维半无限长压电杆在一端受到热冲击时的边值问题。借助拉普拉斯正，反变换技术，在所考虑时间非常短的情况下，对问题进行了求解，得到了压电杆上的位移，压力及温度分布的近似解析解，发现应力及温度分布中分别存在两个阶跃点，并给出了算例。     

多层压电材料层合板的精确解

抛弃有关位移和应力的所有假设，直接从三维弹性力学理论的静电学理论，先导出正交各向异性压电材料板的状态方程，由此得到四边简支压电材料板的状态主程，再根据矩阵分析理论，建立了单层压电材料板的上下表面状态量之间的关系，进一步建立了多层压电板上，下表面状态量之间关系式，利用上下表面已知状态量，得到上表面未知状态的求解方程解。通过求解方程组，便得上表面未知状态量，最终可以得到任意位置处状态量，最后，同时给出了四边简支，两层不同压电材料组成，不同纵横比的层合板受正弦分布载荷作用下的精确解，其结果与现有解比较，吻合较好。     

A GENERALIZED PIEZOELECTRIC-THERMOELASTIC THEORY WITH STRAIN RATE AND ITS APPLICATION 1)

工程中大量材料的形变介于弹性与黏性之间, 既具有弹性固体特性, 又具有黏性流体特点, 即为黏弹性. 黏弹性使得材料出现很多力学松弛现象, 如应变松弛、滞后损耗等行为. 在研究受热载荷作用的多场耦合问题的瞬态响应时, 考虑此类问题中的热松弛和应变松弛现象, 对准确描述其瞬态响应尤为重要. 针对广义压电热弹问题的瞬态响应, 尽管已有学者建立了考虑热松弛的广义压电热弹模型, 但迄今, 尚未计入应变松弛. 本文中, 考虑到材料变形时的应变松弛, 通过引入应变率, 在Chandrasekharaiah广义压电热弹理论的基础之上, 经拓展, 建立了考虑应变率的广义压电热弹理论. 借助热力学定律, 给出了理论的建立过程并得到了相应的状态方程及控制方程. 在本构方程中, 引入了应变松弛时间与应变率的乘积项, 同时, 分别在本构方程和能量方程中引入了热松弛时间因子. 其后, 该理论被用于研究受移动热源作用的压电热弹一维问题的动态响应问题. 采用拉普拉斯变换及其数值反变换, 对问题进行了求解, 得到了不同应变松弛时间和热源移动速度下的瞬态响应, 即无量纲温度、位移、应力和电势的分布规律, 并重点考察了应变率对各物理量的影响效应, 将结果以图形形式进行了表示. 结果表明: 应变率对温度、位移、应力和电势的分布规律有显著影响.     

考虑应变率的广义压电热弹理论及其应用

工程中大量材料的形变介于弹性与黏性之间, 既具有弹性固体特性, 又具有黏性流体特点, 即为黏弹性. 黏弹性使得材料出现很多力学松弛现象, 如应变松弛、滞后损耗等行为. 在研究受热载荷作用的多场耦合问题的瞬态响应时, 考虑此类问题中的热松弛和应变松弛现象, 对准确描述其瞬态响应尤为重要. 针对广义压电热弹问题的瞬态响应, 尽管已有学者建立了考虑热松弛的广义压电热弹模型, 但迄今, 尚未计入应变松弛. 本文中, 考虑到材料变形时的应变松弛, 通过引入应变率, 在Chandrasekharaiah广义压电热弹理论的基础之上, 经拓展, 建立了考虑应变率的广义压电热弹理论. 借助热力学定律, 给出了理论的建立过程并得到了相应的状态方程及控制方程. 在本构方程中, 引入了应变松弛时间与应变率的乘积项, 同时, 分别在本构方程和能量方程中引入了热松弛时间因子. 其后, 该理论被用于研究受移动热源作用的压电热弹一维问题的动态响应问题. 采用拉普拉斯变换及其数值反变换, 对问题进行了求解, 得到了不同应变松弛时间和热源移动速度下的瞬态响应, 即无量纲温度、位移、应力和电势的分布规律, 并重点考察了应变率对各物理量的影响效应, 将结果以图形形式进行了表示. 结果表明: 应变率对温度、位移、应力和电势的分布规律有显著影响.      

功能梯度压电圆板自由振动问题的三维精确分析

本文对周边为广义刚性滑动和广义简支两种边界条件下的功能梯度压电材料圆板自由振动问题进行分析。根据轴对称横观各向同性压电材料基本方程，并利用有限Hankel变换得到了功能梯度压电材料圆板的状态空间方程。假设材料的机械和电学性质均沿板厚方向按统一的指数函数形式梯度分布，从而获得了周边为广义刚性滑动和广义弹性简支两种边界条件下功能梯度压电圆板自由振动问题的三维精确频率方程，该方程是一个关于自由振动频率的超越方程，通过求解该超越方程可得到在不同板厚以及不同的材料性质梯度变化情况下的圆板自由振动频率值，结果表明在相同的材料性质梯度变化情况下频率均随着板厚增加而增大，而在相同的板厚情况下频率则随材料性质梯度变化指数的增大而减小的结论。     

分数阶广义热弹性理论下中空柱热弹性分析

基于Laplace变换及特征值法,推导并给出了分数阶广义热弹性理论下中空柱内表面作用有热冲击情况的解析解,通过Laplace数值逆变换法求解得到了位移场、温度场、应力场的分布规律。结果表明:特征值法能准确给出Laplace域内方程组的解;分数阶参数对温度场和应力场有较大影响,对位移场影响较小。作为广义热弹性理论的一种推广,在处理热传导问题时,通过分数阶广义热弹性理论进行研究更科学、全面。     

Variation principle of piezothermoelastic bodies, canonical equation and homogeneous equation

Combining the symplectic variations theory,the homogeneous control equa- tion and isoparametric element homogeneous formulations for piezothermoelastic hybrid laminates problems were deduced.Firstly,based on the generalized Hamilton variation principle,the non-homogeneous Hamilton canonical equation for piezothermoelastic bod- ies was derived.Then the symplectic relationship of variations in the thermal equilibrium formulations and gradient equations was considered,and the non-homogeneous canoni- cal equation was transformed to homogeneous control equation for solving independently the coupling problem of piezothermoelastic bodies by the inceusement of dimensions of the canonical equation.For the convenience of deriving Hamilton isoparametric element formulations with four nodes,one can consider the temperature gradient equation as constitutive relation and reconstruct new variation principle.The homogeneous equa- tion simplifies greatly the solution programs which are often performed to solve non- homogeneous equation and second order differential equation on the thermal equilibrium and gradient relationship.     

The effect of rotation on generalized micropolar thermoelasticity for a half-space under five theories

A model of the equations of a two-dimensional problem in a micropolar thermoelastic medium for a half-space whose surface is free and subjected to an instantaneous thermal point source is studied. The entire elastic medium is rotating with a uniform angular velocity. The formulation is applied under five theories of the generalized thermoelasticity: Lord–Shulman with one relaxation time, Green–Lindsay with two relaxation times, Green–Naghdi theory (of type II) without energy dissipation and Chandrasekharaiah–Tzou theory with dual-phase-lag, as well as the coupled theory. The normal mode analysis is used to obtain the exact expressions for the considered variables. The distributions of the considered variables are illustrated graphically. Comparisons are made with the results predicted by the five theories in the presence and absence of rotation.     

A direct finite element method study of generalized thermoelastic problems

By using the principle of virtual work, the finite element equations corresponding to the generalized thermoelasticity with two relaxation times (i.e., the G–L theory) are derived. In order to improve the accuracy of integral transform method, especially in two/three-dimension, the equations are solved directly in time-domain. As numerical examples, the developed method is used to investigate the generalized thermoelastic behavior of a slim strip and a semi-infinite plate subject to thermal shock loading. The results demonstrate that the developed method can faithfully predict the deformation of the structure and most importantly the delicate second sound effect of heat conduction in both one and two-dimensional solids whilst it is usually difficult to model, especially in two-dimensional, by using the transform method. This method can also be used to study generalized piezothermoelastic and magnetothermoelastic problems.     

Propagation characteristics of thermoelastic waves in piezoelectric (6 mm class) rotating plate

The propagation of Lamb waves in a homogeneous, transversely isotropic (6 mm class), piezothermoelastic plate rotating with uniform angular velocity about normal to its boundary has been investigated. The generalized (non-classical) theories of thermoelasticity in contrast to Sharma and Pal [Sharma, J.N., Pal, M., 2004. Lamb wave propagation in transversely isotropic piezothermoelastic plate. J. Sound Vib. 270, 587–610] have been used to investigate the problem. The surfaces of the plate are subjected to stress free, thermally insulated/isothermal and electrically shorted boundary conditions. Secular equations for wave propagation modes in the plate are derived from a coupled system of governing partial differential equations of linear piezothermoelasticity. After obtaining the complex characteristic roots with the help of Descartes' algorithm, the transcendental secular equations have been solved by functional iteration numerical technique to compute phase velocity and attenuation coefficient. Finally, in order to illustrate the analytical development, numerical solution of secular equations is carried out for PZT-5A piezo-thermoelastic material. The corresponding simulated results of various physical quantities such as phase velocity, attenuation coefficients, specific loss factor of energy dissipation, thermo-mechanical coupling factor and relative frequency shifts have been presented graphically for both rotating and non-rotating plates for comparison purpose. There is a scope for extension of the present work to other classes of piezo/pyroelectric crystals. The study will be useful in design and construction of gyroscope, rotation sensors, temperature sensors and other pyro/piezoelectric surface acoustic wave (SAW) devices.     

厚板的高阶剪切变形理论研究

已有多种厚板理论和高阶剪切变形模型,但仍需要进一步研究以更加完善.首先根据平均转角及上下表面剪应力自由这两个条件,提出了具有统一高阶剪切变形模型的中面位移模式,并将之表示为正交分解形式.根据正交特性,定义了板的广义应力;运用板问题应变能密度表示的等价性,提出了与广义应力功共轭的广义应变表示形式,建立了板的本构关系.证明了不同转角定义时虚功原理板理论表示的客观性,以及与三维弹性理论表示的等价性.运用虚功原理,建立了变分自洽的高阶厚板理论和变分渐近的低阶厚板理论,推导了相应的平衡方程及边界条件,分析了与已有板理论的异同.以广义应力形式建立了厚板理论的平衡方程,厘清了不同转角表示时板理论间的关系、低阶厚板理论与高阶厚板理论间的关系以及剪切系数计算等若干基本问题.对圣维南扭转问题的求解证明了该理论的正确性.      

A hierarchy of dynamic equations for solid isotropic circular cylinders

This work considers homogeneous isotropic circular cylinders adopting a power series expansion method in the radial coordinate. Equations of motion together with consistent sets of end boundary conditions are derived in a systematic fashion up to arbitrary order using a generalized Hamilton’s principle. Time domain partial differential equations are obtained for longitudinal, torsional, and flexural modes, where these equations are asymptotically correct to all studied orders. Numerical examples are presented for different sorts of problems, using exact theory, the present series expansion theories of different order, and various classical theories. These results cover dispersion curves, eigenfrequencies and the corresponding displacement and stress distributions, as well as fix frequency motion due to prescribed end displacement or lateral distributed forces. The results illustrate that the present approach may render benchmark solutions provided higher order truncations are used, and act as engineering cylinder equations using low order truncation.     

Generalized thermoelastic analysis of the stress-focusing effect in a solid cylinder under rotationally asymmetrical instantaneous thermal loading

Summary  This paper investigates the stress-focusing effect in an infinitely long cylinder under rotationally asymmetrical instantaneous thermal loading on the basis of the generalized thermoelastic Lord–Shulman (L-S) and Green–Lindsay (G-L) theories. Combined forms of the governing equations of both theories are given in a cylindrical coordinate system. The two-dimensional generalized thermoelastic problems are solved by numerical inversion of Laplace transform. Calculations have been performed to find distributions of thermal stresses on the basis of the L-S theory. Stress-focusing phenomena under different heating conditions are presented. The effects of thermomechanical coupling and relaxation time on the stress-focusing phenomena as well as the singularity of stresses are discussed. Received 15 November 2000; accepted for publication 15 November 2001     

Transient piezothermoelasticity for a cylindrical composite panel composed of angle-ply and piezoelectric laminae

This paper is concerned with the theoretical treatment of transient piezothermoelastic problem is developed for a cylindrical composite panel composed of angle-ply laminae and piezoelectric material of crystal class mm2, subject to non-uniform heat supply in the circumferential direction. We obtain the exact solution for the two-dimensional temperature change in a transient state, and transient piezothermoelastic response of a simple supported cylindrical composite panel under the state of generalized plane deformation. As an example, numerical calculations are carried out for an angle-ply laminated composite panel made of alumina fiber reinforced aluminum composite, associated with a piezoelectric layer of a cadmium selenide solid. Some numerical results for temperature change, displacement, stress and electric potential distributions in a transient state are shown in figures. Furthermore, the influence of thickness of the angle-ply laminate is investigated.     

Effect of rotation on plane waves in generalized thermo-microstretch elastic solid with a relaxation time

The present paper is aimed at studying the effect of rotation on the general model of the equations of generalized thermo-microstretch for a homogeneous isotropic elastic half-space solid whose surface is subjected to a Mode-I Crack problem considered. The problem is in the context of the generalized thermoelasticity Lord-?hulman??s (L-S) theory with one relaxation time, as well as the classical dynamical coupled theory (CD) The normal mode analysis is used to obtain the exact expressions for the displacement components, force stresses, temperature, couple stresses and microstress distribution. The variations of the considered variables through the horizontal distance are illustrated graphically. Comparisons are made with the results in the presence and absence of rotation and in the presence and absence of microstretch constants between the two theories.