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

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

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

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

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

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

不同理论下广义压电热弹性问题的有限元求解

基于G-L和L-S广义压电热弹性理论研究了无限大厚压电板在上下表面受到条带状热冲击时的广义压电热弹性问题。在时间非常短的情况下，为避免积分变换求解带来的精度丢失，采用有限元方法对问题在时间域进行直接求解，获得压电板在热冲击作用下的温度、位移、应力及电势等，并将结果与经典压电热弹性理论进行比较。结果表明，直接求解方法可以准确描述热在介质中以有限的速度传播。     

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

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

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

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

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

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

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

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

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

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

应变控制的热机械疲劳行为的数值模拟

根据高温合金材料的力学性能，以弹粘塑性本构模型为基础，用数值模拟方法研究材料的热机械疲劳循环特性，模型将应变分为弹性应变、温度应变和粘塑性应变三部分，认为材料在高温循环载荷下呈现明显的弹粘塑性特征，根据虚位移原理建立轴对称体的弹粘塑性计算有限元格式，对于循环机械载荷和循环温度载荷，程序中采用了增量法迭代求解，在非线性项中不仅考虑了机械载荷增量的影响，同时也考虑了温度增量的影响，根据应变控制热机械疲劳的特点，发展了应变增量法的有限元计算方法、通过数值模拟，得到材料在各种循环载荷下的应力—应变响应，数值模拟较好地反映了粘塑性变形过程以及温度变化的效应，所描述的不可逆系统在某一时刻的状态完全由当时的状态参数、内变量、承载时间及塑性应变累积量决定，对带缺口试件的模拟结果显示了程序对复杂轴对称结构进行热机械疲劳计算的有效性。     

EFFECT OF TEMPERATURE-DEPENDENT PROPERTIES ON THERMOELASTIC PROBLEMS WITH THERMAL RELAXATIONS

Based on the generalized thermoelasticity proposed by Green and Lindsay, the dynamic response of generalized thermoelastic problems with temperature-dependent material properties is investigated. The governing equations are formulated and found to be nonlinear because of the temperature-dependence of properties. Owing to the nonlinearity of the governing equations, the finite element method is resorted to for solution. The results obtained show that the temperature-dependent properties influence the variables considered by reducing their magnitudes. This indicates that taking the temperature-dependence of properties into consideration in the investigation of generalized thermoelastic problems is necessary and practical for accurately predicting the thermoelastic behavior.     

Nonlinear transient thermal stress analysis of thick-walled FGM cylinder with temperature-dependent material properties

In this paper, the problem of generalized thermoelasticity in a thick-walled FGM cylinder with one relaxation time is presented. The material properties are taking as function of temperature and graded in the radial. Due to the nonlinearity of the governing equations, finite element method is adopted to solve such problem. Both the inner and outer curved surfaces of the cylinder are stress free while the inner surface is subjected to thermal shock, while the outer surface is thermally isolated. The effects of temperature-dependent properties, volume fraction parameter and the thermal relaxation time on the physical quantities behavior are evaluated. Results confirm the efficiency of the present algorithm and reveal the significant effects of the temperature-dependent of the material properties.     

Generalized thermo-piezoelectric problems with temperature-dependent properties

The model of the equation of generalized thermo-piezoelectricity in an isotropic elastic medium with temperature-dependent mechanical properties is established. The modulus of elasticity is taken as a linear function of reference temperature. The state-space approach is adopted for the solution of one-dimensional problems in the absence or presence of heat sources. A numerical technique is employed to obtain the solution in the physical domain. The results are given and illustrated graphically. A comparison is made with results obtained in case of temperature-independent modulus of elasticity.     

Hybrid numerical scheme for nonlinear two-dimensional phase-change problems with the irregular geometry

A hybrid numerical scheme combining the Laplace transform and control-volume methods is presented to solve nonlinear two-dimensional phase-change problems with the irregular geometry. The Laplace transform method is applied to deal with the time domain, and then the control-volume method is used to discretize the transformed system in the space domain. Nonlinear terms induced by the temperature-dependent thermal properties are linearized by using the Taylor series approximation. Control-volume meshes in the solid and liquid regions during simulations are generated by using the discrete transfinite mapping method. The location of the phase-change interface and the isothermal distributions are determined. Comparison of these results with previous results shows that the present numerical scheme has good accuracy for two-dimensional phase-change problems. Received on 17 October 1996     

Methods of fluid–structure coupling in frequency and time domains using linearized aerodynamics for turbomachinery

Two methods of fluid–structure coupling for turbomachinery are presented, the first one in the frequency domain and the second in both frequency and time domains. In both methods, the structure and the fluid are assumed to have circumferential cyclic symmetric properties and the unsteady aerodynamic forces are assumed to be linear in terms of the structural displacements. The motion equation of the reference sector in the travelling wave coordinates is projected on the complex eigenmodes for each phase number. The generalized unsteady aerodynamic forces are computed by solving the Euler equations and by assuming the structural motion to be harmonic with a constant phase angle between two adjacent sectors. In the frequency domain, the complex, nonlinear eigenvalue problem for the aeroelastic stability analysis is solved iteratively either by the double scanning method or by using Karpel's minimum state smoothing of the aerodynamic coefficient matrix. In the time domain, Karpel's smoothing method is used to obtain an approximation of the generalized unsteady aerodynamic forces by means of auxiliary state variables. These coupling methods are tested on a compressor blade row and the good agreement obtained between their results and those of the direct coupling method shows that the proposed numerical methods, already used in aircraft applications, are adapted to turbomachinery.     

Thermally induced vibration and stability of laminated composite plates with temperature-dependent properties

This paper presents a study on the buckling and vibration of initially stressed composite plates with temperature-dependent material properties in thermal environments. The initial stress is taken to be a combination of a pure bending stress and an axial stress. The temperature distribution in the plate is assumed to be uniform and linear in the transverse direction. The governing equations including the transverse shear deformation effects are established using the variational method. The effects of various parameters on the buckling and vibration behaviors of laminated plates with respective temperature-dependent and temperature-independent material properties are investigated. The buckling load and natural frequency are sensitive to the thermal stresses and initial stresses. Numerical results reveal that temperature-dependent material properties should be considered in the buckling and vibration analysis for laminated plates under thermal conditions.     

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.     

Application of a fourth-order relaxation scheme to hyperbolic systems of conservation laws

A fourth-order relaxation scheme is derived and applied to hyperbolic systems of conservation laws in one and two space dimensions. The scheme is based on a fourth-order central weighted essentially nonoscillatory (CWENO) reconstruction for one-dimensional cases, which is generalized to two-dimensional cases by the dimension-by-dimension approach. The large stability domain Runge-Kutta-type solver ROCK4 is used for time integration. The resulting method requires neither the use of Riemann solvers nor the computation of Jacobians and therefore it enjoys the main advantage of the relaxation schemes. The high accuracy and high-resolution properties of the present method are demonstrated in one- and two-dimensional numerical experiments. The project supported by the National Natural Science Foundation of China (60134010) The English text was polished by Yunming Chen.