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

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

广义热弹模型的热力学基础与瞬态响应

微纳科技的快速发展与超短脉冲激光技术的广泛运用, 对描述微纳尺度超快热冲击的广义热传导及其热弹耦合理论提出迫切需求. 基于拓展热力学原理, 本文建立了考虑热传导双相滞后效应和高阶热流率的广义热弹耦合理论. 类比于力学领域黏弹性本构关系的串联、并联模型, 并受Green-Naghdi (GN)广义热传导模型启发, 本文提出了热学“弹性”单元和“黏性”单元模型, 并采用串联、并联方法实现了Cattaneo-Vernotte (CV)、GN、双相滞后(DPL)和Moore-Gibson-Thompson (MGT) 热传导模型的重构. 理论推导进一步表明, 本文新建模型对应于热学Burgers模型, 并得到了新模型中各相位滞后中松弛时间之间的比例关系. 运用拉普拉斯变换方法, 研究了一维结构受边界热冲击和移动热源作用下的瞬态响应, 计算结果表明: 新模型克服了热波速度无限大的悖论; 仅有边界热冲击载荷时, 新模型得到的响应结果均较大, 响应范围最小; 相比于无热源作用情形, 受移动热源作用时, 新模型会产生更大的峰值响应. 新模型与经典弹性理论耦合构建了广义热弹性理论, 运用该理论, 可以清晰观察到在热波和弹性波波前的应力突变. 理论方面, 本文推动了拓展热力学与连续介质力学的结合, 对于远离平衡态极端力学基础理论问题的研究具有启发意义; 应用方面, 本文研究结果可为激光等移动热源作用下材料的瞬态响应分析提供理论基础和数值方法.      

由运动内热源引起的磁热黏弹性问题的研究

在具有两个热松弛时间的广义热弹性理论下, 研究了处于定常磁场中的均布各向同性黏弹性半空间中, 由以均匀速度运动的线热源引起的瞬态波问题. 通过引入黏弹性向量势和热黏弹性标量势,问题退化为求解3个偏微分方程. 运用Laplace变换(对时间变量)和Fourier变换(对一个空间变量), 得到了变换域内应力和位移的解析表达式. 采用级数展开法, 得到了边界位移在小时间范围内的近似解, 给出了解的近似范围, 同时还研究了两种特例：(1)热源静止不动, (2)不考虑热松弛时间的影响. 最后对于丙烯酸塑料介质给出了数值结果.     

短纤维增强三元乙丙橡胶横观各向同性黏-超弹性本构模型

短纤维增强三元乙丙橡胶包覆薄膜,是一种应用于固体火箭发动机缠绕包覆装药的新型复合材料.为了描述其在工作过程中受振动、冲击等载荷作用时的力学行为,基于黏弹性理论和纤维增强连续介质力学理论,提出了一种考虑应变率强化效应的横观各向同性黏-超弹本构模型.模型中应变能函数被分解为超弹性应变能和黏性应变能,其中超弹性应变能包括表征各向同性的橡胶基体应变能和表征各向异性的纤维拉伸应变能,黏性应变能采用表征橡胶和纤维黏性响应的宏观唯象模型.选取表征各应变能的函数形式,经过数学变换、替代、叠加,求解确定最终的应力应变形式,明确模型参数获取的具体步骤,将预测结果与实验结果对比分析,准确性较高.研究表明:该模型能有效预测材料在低应变率下纤维方向为0?～45?的非线性率相关力学特性;模型形式易于实现有限元开发,对固体火箭发动机装药结构完整性分析具有参考价值.     

NEPE推进剂低高应变率下改进的黏-超弹本构模型

为研究低高应变率条件下NEPE推进剂的力学特性,通过电子万能试验机和分离式霍普金森杆装置,对NEPE推进剂进行了准静态和冲击实验,得到了不同应变率下（1.667×10?4～4 500 s?1）的应力-应变曲线。实验结果表明NEPE推进剂具有明显的非线性弹性和应变率敏感性,随着应变率的增加,材料的强度、屈服应力和弹性模量显著增加,与低应变率相比,高应变率条件下材料的应变率敏感性更高。在高速冲击下材料内部瞬间产生大量热量无法及时散发出去,使得材料内部温度升高,导致材料出现软化效应,力学性能降低。本文建立了一个非线性黏超弹本构模型,其中采用Rivlin应变能函数来描述稳态超弹响应部分,采用积分型本构模型来描述材料的动态黏弹性响应部分,考虑到松弛时间具有应变率相关性,本文采用了一个率相关松弛函数来替代传统的Prony级数形式。使用极慢速压缩实验数据对本构模型中的超弹部分进行拟合获得超弹参数,然后用准静态和动态实验数据对本构模型进行拟合得出其他参数。不同应变率下的预测曲线与实验曲线具有较好的重合度,证明了该模型可以很好地描述低高应变率下NEPE推进剂的力学特性。     

坡形加热下的二维广义磁热黏弹性问题研究

运用具有一个热松弛时间的广义热黏弹性理论，研究了处于均布磁场中的二维磁热黏弹 性问题. 运用Laplace变换(对时间变量)和Fourier变换(对于一个空间变量)，得到了变 换域内场量的精确表达式，并把结果应用到表面受到坡形加热的半空间问题. 应用 数值逆变换得到了时间-空间域内场量的解，对丙烯酸塑料 给出场量的响应图. 并把运用广义热黏弹性理论所得的结果与传统热黏弹性理 论及热弹性理论下的结果进行了比较.     

高应变率下航空透明聚氨酯的动态本构模型

采用低阻抗分离式霍普金森压杆对航空透明聚氨酯进行了高应变率下的动态力学性能测试,得到的应力应变曲线表现出了显著的非线性黏弹性特征。基于本构理论和实验数据,构建了航空透明聚氨酯的松弛时间应变相关的超黏弹性本构形式。该本构模型由2部分组成:一部分表征准静态下的超弹性行为,另一部分描述非线性应变率的相关特性。利用超黏弹性本构模型对不同应变率下航空透明聚氨酯的动态应力应变曲线进行拟合,拟合曲线与实验曲线一致性良好。     

半无限粘弹杆的广义磁热粘弹问题

应用Lord-Shulman(L-S)和Green-Lindsay(G-L)广义热弹性理论,研究了在磁场中受移动热源作用的半无限长均质各向同性粘弹杆的磁热粘弹动态响应,并与经典耦合理论进行了对比.给出了杆的广义磁热粘弹耦合的控制方程,借助拉普拉斯积分变换及其数值反变换对控制方程进行了求解,计算得到了杆内温度、应力和位移的分布规律.研究结果表明:时间、热源移动速度和磁场大小对以上分布规律都有一定的影响.     

广义热弹性问题研究进展

本文总结了广义热弹性问题最近10年的研究进展, 包括不同类型广义热弹耦合问题的研究、考虑磁\!--\!电多场耦合的广义电磁热弹耦合问题研究以及计及扩散效应和黏弹性效应的广义热弹性理论的发展、广义热弹性问题基本求解方法等, 通过总结, 使读者对广义热弹性问题的研究现状及发展趋势有较全面的认识, 帮助研究人员进一步开展广义热弹性问题更高层次的研究.      

超短激光脉冲加热薄板的广义热弹扩散问题

由于超短激光脉冲具有功率密度高、持续时间短、加工精度高等优势, 近年来被广泛应用于超精细加工、光学储存和微电子器件制造等领域. 本文基于L-S型广义热弹扩散理论, 建立了考虑材料记忆依赖效应和空间非局部效应的记忆依赖型非局部广义热弹扩散耦合理论, 它能够准确预测几何尺寸与内部特征尺寸相近结构的热弹扩散瞬态响应. 推导了所建理论的控制方程, 并基于拉普拉斯积分变换获得了控制方程的解. 作为算例, 利用所建理论和求解方法研究了半无限大薄板受非高斯激光脉冲加热和化学冲击联合作用下的热弹扩散瞬态响应问题, 得到了薄板的温度、化学势、位移、应力和浓度等随非局部参数、热时间迟滞因子和扩散时间迟滞因子等参数变化的分布规律. 结果表明: 传热对传质影响显著, 传质对传热影响甚微; 非局部参数对位移、应力影响显著, 对温度、化学势和浓度几乎没有影响. 该理论及求解方法的建立, 旨在实现材料在机械、热、化学势等冲击作用下传热传质瞬态响应的准确预测.      

Thermo-electro-elastic transient responses in piezoelectric hollow structures

The paper presents an analytical method to solve thermo-electro-elastic transient response in piezoelectric hollow structures subjected to arbitrary thermal shock, sudden mechanical load and electric excitation. Volterra integral equation of the second kind caused by interaction between elastic deformation and electric field is solved by using an interpolation method. Thus, the exact expressions for the transient responses of displacement, stresses, electric displacement and electric potential in the piezoelectric hollow structures are obtained by means of Hankel transform, Laplace transform, and their inverse transforms. In Section 2, based on spherical coordinates, the governing equation of thermo-electro-elastic transient responses in a piezoelectric hollow sphere is found and the associated numerical results are carried out. In Section 3, based on cylindrical coordinates, the governing equation of thermo-electro-elastic transient responses in a non-homogeneous piezoelectric hollow cylinder is found and the corresponding numerical results are carried out. The results carried out may be used as a reference to solve other transient coupled problems of thermo-electro-elasticity in piezoelectric structures.     

Transient thermopiezoelectric response of a one-dimensional functionally graded piezoelectric medium to a moving heat source

The thermopiezoelectricity problem of a one-dimensional (1-D), finite length, functionally graded medium excited by a moving heat source is investigated in this paper. The Lord and Shulman theory of generalized coupled thermoelasticity is employed to account for both the finite speed of thermal waves and coupling of temperature field with displacement and electric fields. Except thermal relaxation time and specific heat, which are taken to be constant for simplicity, all other properties are assumed to vary exponentially along the length through an arbitrary non-homogeneity index. Laplace transform has been used to eliminate the time effect, and three coupled fields, namely, displacement, temperature, and electric fields are obtained analytically in the Laplace domain. The solutions are then inverted to time domain using a numerical Laplace inversion method. Numerical examples are displayed to illustrate the effects of non-homogeneity index, length and thermal relaxation time on the results. When the medium is homogeneous, the results of the current paper are reduced to exactly the same results available in the literature.     

Exact Electromagnetothermoelastic Solution for a Transversely Isotropic Piezoelectric Hollow Sphere Subjected to Arbitrary Thermal Shock

This paper presents analytical study for electromagnetothermoelastic transient behavior of a transversely isotropic hollow sphere, placed in a uniform magnetic field, subjected to arbitrary thermal shock. Exact solutions for the transient responses of stresses, perturbation of magnetic field vector, electric displacement and electric potential in the transversely isotropic piezoelectric hollow sphere are obtained by means of the Hankel transform, the Laplace transform and their inverse transforms. An interpolation method is used to solve the Volterra integral equation of the second kind caused by interactions among electric, magnetic, thermal and elastic fields. From the sample numerical calculations, it is seen that the present method is suitable for the transversely isotropic hollow sphere, placed in a uniform magnetic field, subjected to arbitrary thermal shock. Finally, the result can be used as a reference to solve other transient coupling problems of electromagnetothermoelasticity.     

热冲击载荷作用下的含球形空腔的广义热弹性功能梯度球形各向同性体

This paper is concerned with the determination of thermoelastic displacement, stress and temperature in a functionally graded spherically isotropic infinite elastic medium having a spherical cavity, in the context of the linear theory of generalized thermoelasticity with two relaxation time parameters （Green and Lindsay theory）. The surface of cavity is stress-free and is subjected to a time-dependent thermal shock. The basic equations have been written in the form of a vector-matrix differential equation in the Laplace transform domain, which is then solved by an eigenvalue approach. Numerical inversion of the transforms is carried out using the Bellman method. Displacement, stress and temperature are computed and presented graphically. It is found that variation in the thermo-physical properties of a material strongly influences the response to loading. A comparative study with a corresponding homogeneous material is also made.     

物性与温度相关材料广义热弹性模型及渐近分析

计及材料物性与温度的相关性,基于Clausius不等式和L-S广义热弹性理论,通过对自由能公式的高阶展开,构建了具有变物性特征的广义耦合热弹性动力学模型。推导了各向同性材料表面受热冲击问题的线性化控制方程组,利用热冲击的瞬时特征,借助于Laplace正、逆变换技术及其极限性质,给出了变物性条件下一维热冲击问题的温度场、位移场和应力场的渐近表达式。通过算例,得到了热冲击作用下各物理场的分布规律以及材料物性与温度相关性对于热弹性响应的影响规律。结果表明：材料物性与温度相关性对于各物理场的阶跃位置、阶跃间隔以及阶跃峰值均产生影响,但值得注意的是,相比于位移场和应力场的显著影响,其对温度场的影响效果并不明显。     

An adaptation of finite linear viscoelasticity theory for rubber-like viscoelasticity by use of a generalized strain measure

In this paper a simplified three-dimensional constitutive equation for viscoelastic rubber-like solids is derived by employing a generalized strain measure and an asymptotic expansion similar to that used by Coleman and Noll (1961) in their derivation of finite linear viscoelasticity (FLV) theory. The first term of the expansion represents exactly the time and strain separability relaxation behavior exhibited by certain soft polymers in the rubbery state and in the transition zone between the glassy and rubbery states. The relaxation spectra of such polymers are said to be deformation independent. Retention of higher order terms of the asymptotic expansion is recommended for treating deformation dependent spectra.Certain assumptions for the solid theory are relaxed in order to obtain a constitutive equation for uncross-linked liquid materials which exhibit large elastic recovery properties.Apart from the strain energyW(I₁,I ₂), which alternatively characterizes the long-time elastic response of solids or the instantaneous elastic response of elastic liquids, only the linear viscoelastic relaxation modulus is required for the first-order theory. Both types of material functions can be obtained, in theory, from simple laboratory testing procedures. The constitutive equations for solids proposed by Chang, Bloch and Tschoegl (1976) and a special form of K-BKZ theory for elastic liquids are shown to be particular cases of the first-order theory.Previously published experimental data on a cross-linked styrene-butadiene rubber (SBR) and an uncross-linked polyisobutylene (PIB) rubber is used to corroborate the theory.     

Nonlinear dependence of viscosity in modeling the rate-dependent response of natural and high damping rubbers in compression and shear: Experimental identification and numerical verification

The rate-dependent behavior of filled natural rubber (NR) and high damping rubber (HDR) is investigated in compression and shear regimes. In order to describe the viscosity-induced rate-dependent effects, a constitutive model of finite strain viscoelasticity founded on the basis of the multiplicative decomposition of the deformation gradient tensor into elastic and inelastic parts is proposed. The total stress is decomposed into an equilibrium stress and a viscosity-induced overstress by following the concept of the Zener model. To identify the constitutive equation for the viscosity from direct experimental observations, an analytical scheme that ascertains the fundamental relation between the inelastic strain rate and the overstress tensor of the Mandel type by evaluating simple relaxation test results is proposed. Evaluation of the experimental results using the proposed analytical scheme confirms the necessity of considering both the current overstress and the current deformation as variables to describe the evolution of the rate-dependent phenomena. Based on this experimentally based motivation, an evolution equation using power laws is proposed to represent the effects of internal variables on viscosity phenomena. The proposed evolution equation has been incorporated in the finite strain viscoelasticity model in a thermodynamically consistent way. Simulation results for simple relaxation tests, multi-step relaxation tests and monotonic tests at different strain rates using the developed model show an encouraging correlation with the experiments conducted on HDR and NR in both compression and shear regimes. Finally, an approach to extend the proposed evolution equation for rate-dependent cyclic processes is proposed. The simulation results are critically compared with the experiments.     

Magneto–thermo–electro–elastic transient response in a piezoelectric hollow cylinder subjected to complex loadings

The article presents an analytical solution for magneto–thermo–electro–elastic problems of a piezoelectric hollow cylinder placed in an axial magnetic field subjected to arbitrary thermal shock, mechanical load and transient electric excitation. Using an interpolation method solves the Volterra integral equation of the second kind caused by interaction among magnetic, thermal, electric and mechanical fields, the electric displacement is determined. Thus, the exact expressions for the transient responses of displacement, stresses, electric displacement, electric potential and perturbation of the magnetic field vector in the piezoelectric hollow cylinder are obtained by means of Hankel transforms, Laplace transforms, and inverse Laplace transforms. From sample numerical calculations, it is seen that the present method is suitable for a piezoelectric hollow cylinder subjected to arbitrary thermal shock, mechanical load and transient electric excitation, and the result carried out may be used as a reference to solve other transient coupled problems of magneto–thermo–electro–elasticity.