非傅里叶热弹性的时域间断迦辽金有限元方法

基于Lord-Shulman非傅里叶热弹性模型,提出了采用修正的时域间断迦辽金有限元方法(time discontinuousGalerkin finite element method, DGFEM)求解方法. DGFEM对温度场、位移场基本未知向量及其时间导数向量在时域中分别插值;在最终的求解公式中,引入了人工阻尼. 数值结果显示所发展的DGFEM 较好地捕捉了波的间断并消除了热冲击作用下虚假的数值振荡,能够良好地模拟热弹性问题并具有较高的精度.     

分数阶理论下实心球体受热冲击作用的渐进分析

基于分数阶广义热弹性理论,针对实心球体在外表面受均匀热冲击作用下的一维广义热弹性问题进行研究分析. 利用热冲击的瞬时特征,借助于Laplace 正、反变换技术及柱函数的渐进性质,推导了热冲击作用周期内位移场、温度场和应力场的渐进表达式. 通过计算,得到了不同传热能力下受热冲击作用时热波、热弹性的传播规律以及位移场、温度场及应力场的分布规律. 结果表明:分数阶参数取值的不同,热波、热弹性波的传播以及各物理场的分布均有所不同,分数阶参数可视为延迟时间的影响因子,通过改变延迟效应对热弹性行为的影响来改变热冲击的作用效果.      

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

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

材料强度对层板抗弹性能和损伤特性的影响

根据纤维增强复合材料宏细观结构特征,基于ABAQUS软件平台,建立了层合板高速冲击损伤三维有限元分析模型。该模型在复合材料层间引入界面单元模拟层间分层,采用三维粘弹性本构,结合Hashin失效准则模拟单层板面内损伤.利用该模型,深入研究了复合材料层板的抗弹性能和损伤特性,数值分析结果与实验结果吻合良好,证明了该方法的合理有效性。通过数值分析,详细探讨了材料强度参数对层板抗弹性能和损伤特性的影响规律,获得了一些有价值的结论。     

热冲击下物性与温度相关性对热弹性响应的渐进分析

计及材料物性与温度的相关性,基于Green-Naghdi能量无耗散广义热弹性理论(G-N II理论),对热冲击下具有变物性特征材料的热弹性响应进行了求解分析。借助Laplace正、反变换技术以及Krichhoff变换,在热物性参数随真实温度呈线性规律的前提下,推导了半无限大体受热冲击作用时热弹性响应的解析表达式,通过求解分析,得到了热冲击下热波、热弹性波的传播规律,位移场、温度场以及应力场的分布情况,以及物性随温度相关性对热弹性响应的影响效果。结果表明：当考虑材料物性随温度的变化时,热波、热弹性波的传播以及各物理场的分布均受到不同程度的影响,且物性随温度相关性对热弹性响应的作用效果将受到材料热-力耦合特性的影响。     

热冲击下温度相关物性对热弹性响应的渐进分析

计及材料物性与温度的相关性,基于Green-Naghdi能量无耗散广义热弹性理论(G-N II理论),对热冲击下具有变物性特征材料的热弹性响应进行了求解分析.借助Laplace正、反变换技术以及Krichhoff变换,在热物性参数随真实温度呈线性规律的前提下,推导了半无限大体受热冲击作用时热弹性响应的解析表达式,通过求解分析,得到了热冲击下热波、热弹性波的传播规律,位移场、温度场以及应力场的分布情况,以及物性随温度相关性对热弹性响应的影响效果.结果表明:当考虑材料物性随温度的变化时,热波、热弹性波的传播以及各物理场的分布均受到不同程度的影响,且物性随温度相关性对热弹性响应的作用效果将受到材料热-力耦合特性的影响.     

实心圆柱外表面受热冲击作用的轴对称平面应变热弹性动力学解

利用分离变量法，成功地给出了各向同性实心圆柱外表面受热冲击作用的轴对称平面应变热弹性动力学问题的解析解．运用此方法，可以避免积分变换且易于实现数值计算．在数值结果中，与拟静态问题和均匀热冲击情形进行了比较，表明惯性项对中心处应力响应影响不大，而热传导的过程能大大减低中心处的热应力．     

圆柱外表面受热冲击问题的广义热弹性分析

基于 L-S 广义热弹性理论, 针对实心圆柱体在外表面受均匀热冲击作用下的一维广义热弹性问题进行研究分析. 利用热冲击的瞬时特征, 借助于 Laplace 正、反变换技术及柱函数的渐近性质, 推导了热冲击作用周期内温度场、位移场和应力场的渐近表达式. 通过计算, 得到了热冲击条件下各物理场的分布规律以及延迟效应和耦合效应对热弹性响应的影响规律. 结果表明: 当考虑延迟效应和耦合效应时, 热扰动将以两组速度不同的波的形式向前传播, 延迟效应和耦合效应对各物理场的建立时间, 阶跃间隔和阶跃峰值均产生影响, 且延迟效应和耦合效应均在一定程度上削弱了热冲击的作用效果.      

求解三维弹性力学问题的快速耦合方法

将多种数值方法耦合,充分利用各种方法的优点建立新的数值方法,是求解三维复杂问题的有效途径之一．本文将无单元Galerkin (Element-Free Galerkin, EFG)方法、有限元法和维数分裂法耦合,提出了求解三维弹性力学问题的快速耦合方法(Fast Hybrid Method, FHM)．将三维弹性力学问题分裂为若干个二维平面问题,对于每个二维问题采用罚函数法施加边界条件,并推导其相应的积分弱形式,引入Shepard基函数的移动最小二乘法建立形函数,进而推导二维平面问题的离散方程．第三个方向上采用有限元法将这些二维离散方程进行耦合,可以得到原三维弹性力学问题的快速耦合方法数值解的求解公式．通过数值算例验证了本文快速耦合方法求解三维弹性力学问题的收敛性,将数值解与解析解对比,说明了本文方法求解三维弹性力学问题的有效性．     

二维弹性结构入水冲击过程中的流固耦合效应

描述了一个研究弹性结构入水冲击过程中水弹性效应的数值方法,在弹性结构入水冲击过程中,流体域作用在结构上的水动力载荷由边界元法获得,而结构的弹性动力响应则由有限元方法求解,通过线性给离散Ｂｅｒｎｏｕｌｌｉ方程将有限元方程和边界元方程耦合到一起,从而获得了求解流场和结构动力响应的相互耦合的运动方程。在数值考虑了自由表面的非线性边界条件,通过引入射流单元以及最大射流厚度,较好地处理了冲击引起的射流问题。     

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     

A NUMERICAL ANALYSIS OF THREE-DIMENSIONAL THERMOELASTIC STRESS WAVE PROPAGATION PROBLEM

This paper is concerned with the numerical technique based on the method of character-istics for three-dimensional dynamic thermoelastic problems. A numerical example for the three-di-mensional stress wave propagation in a thermoelastic bar of square cross section subjected to both animpact loading and a thermal shock is presented.     

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.