THERMODYNAMIC BASIS AND TRANSIENT RESPONSE OF GENERALIZED THERMOELASTICITY
-
摘要: 微纳科技的快速发展与超短脉冲激光技术的广泛运用, 对描述微纳尺度超快热冲击的广义热传导及其热弹耦合理论提出迫切需求. 基于拓展热力学原理, 本文建立了考虑热传导双相滞后效应和高阶热流率的广义热弹耦合理论. 类比于力学领域黏弹性本构关系的串联、并联模型, 并受Green-Naghdi (GN)广义热传导模型启发, 本文提出了热学“弹性”单元和“黏性”单元模型, 并采用串联、并联方法实现了Cattaneo-Vernotte (CV)、GN、双相滞后(DPL)和Moore-Gibson-Thompson (MGT) 热传导模型的重构. 理论推导进一步表明, 本文新建模型对应于热学Burgers模型, 并得到了新模型中各相位滞后中松弛时间之间的比例关系. 运用拉普拉斯变换方法, 研究了一维结构受边界热冲击和移动热源作用下的瞬态响应, 计算结果表明: 新模型克服了热波速度无限大的悖论; 仅有边界热冲击载荷时, 新模型得到的响应结果均较大, 响应范围最小; 相比于无热源作用情形, 受移动热源作用时, 新模型会产生更大的峰值响应. 新模型与经典弹性理论耦合构建了广义热弹性理论, 运用该理论, 可以清晰观察到在热波和弹性波波前的应力突变. 理论方面, 本文推动了拓展热力学与连续介质力学的结合, 对于远离平衡态极端力学基础理论问题的研究具有启发意义; 应用方面, 本文研究结果可为激光等移动热源作用下材料的瞬态响应分析提供理论基础和数值方法.Abstract: The rapid development of micro/nano technology and the wide application of ultrashort pulsed laser technology have put forward an urgent need for generalized heat conduction and thermoelastic coupling theory to describe ultrafast thermal shock at micro/nano scale. Based on extended thermodynamic principle, a generalized thermoelastic coupling theory considering the dual-phase-lagging effect of heat conduction and the rate of higher order heat flux is established. Inspired by Green-Naghdi (GN) generalized heat conduction model, thermal "elastic" and "viscous" element models are proposed, which are similar to the series and parallel models of viscoelastic constitutive relations in the field of mechanics. The Cattanoe-Vernotte (CV), GN, dual-phase-lag (DPL) and Moore-Gibson-Thompson (MGT) heat conduction models were obtained by series and parallel methods. Theoretical derivation further shows that the newly formulated model corresponds to the Burgers model of heat conduction. In these models, the proportional relationship between the relaxation time of each phase lag is also obtained. Laplace transform method is used to study the transient response of one-dimensional structure under thermal shock and moving heat source. The results show that the present model overcomes the paradox of infinite thermal wave velocity. When the boundary thermal shock load is applied, the results obtained by the new model have higher peak and the smallest affected region. And under the effect of moving heat source, the new model can generate a larger peak response. The new model could coupled with the classical elastic theory and built a generalized thermoelasticity. With this theory, the jump of stress at wavefront of thermal wave and elastic wave can be clearly observed. Theoretically, this paper promotes the combination of extended thermodynamics and continuum mechanics, which is of enlightening significance to the study of fundamental theoretical problems far from equilibrium of extreme mechanics. For applications, this work can provide theoretical basis and numerical method for the transient response analysis under the moving heat sources.
-
引 言
热传导理论是热弹耦合理论的重要基础, 经典傅里叶热传导定律为
$$ {\boldsymbol{q}}({\boldsymbol{r}},t) = - k\nabla T({\boldsymbol{r}},t) $$ (1) 其中,
$k$ 为热传导系数,${\boldsymbol{q}}$ 为热流,$T$ 为温度. 傅里叶定律及其热弹耦合理论被广泛用于描述热力共同作用下材料/结构的热传导和变形行为. 同时, 傅里叶定律表明[1-3]: 热以无限大速度传导, 即一点的温度变化能瞬间传导到无穷远处. 解决这一理论“悖论”成为不断推动热传导方程深入发展的动力. 同时, 飞速发展的激光、集成电路等技术领域, 具有作用载体和作用环境极小、作用时间极短的特点, 使得器件在极小空间、极短时间内产生大量热量. 建立在宏观、近平衡态的经典热传导和热力耦合关系不再适用于此类以微纳尺度和超快热冲击为代表的极端非平衡态问题[4-6].为此, 国内外学者提出了多种广义热传导模型, 并建立了相应的广义热弹性理论. Cattaneo[7]和Vernotte[8]提出了CV (Cattaneo-Vernotte)热传导模型
$$ {\boldsymbol{q}}({\boldsymbol{r}},t) + \tau {\dot{\boldsymbol q}}({\boldsymbol{r}},t) = - k\nabla T({\boldsymbol{r}},t) $$ (2) 其中,
$\tau $ 为松弛时间,${\dot{\boldsymbol q}}$ 为热流率, 松弛时间与热流率的引入表征了热流的滞后效应. Tzou[9]引入了热流与温度梯度的松弛时间, 得到了双相滞后(DPL)模型$$ {\boldsymbol{q}}({\boldsymbol{r}},t + {\tau _q}) = - k\nabla T({\boldsymbol{r}},t + {\tau _T}) $$ (3) 其中,
${\tau _q}$ ,${\tau _T}$ 分别为热流相、温度梯度相的松弛时间. 若左右两端分别按一阶泰勒级数展开, 可得到$$ {\boldsymbol{q}}({\boldsymbol{r}},t) + {\tau _q}{\dot{\boldsymbol q}}({\boldsymbol{r}},t) = - k\left[ {\nabla T({\boldsymbol{r}},t) + {\tau _T}\nabla \dot T({\boldsymbol{r}},t)} \right] $$ (4) Green等[10-11]提出了GN (Green-Naghdi)热传导模型
$$ {\boldsymbol{q}}\left( {{\boldsymbol{r}},t} \right) = - \left[ {{k^*}\nabla \chi \left( {{\boldsymbol{r}},t} \right) + k\nabla T\left( {{\boldsymbol{r}},t} \right)} \right] $$ (5) 其中,
${k^*}$ 为导热比,$\chi $ 为热位移:$\dot \chi = T$ . 当${k^*} = 0$ 时, 该式退化为傅里叶定律. 对式(5)求导, 可得$$ {\dot{\boldsymbol q}}({\boldsymbol{r}},t) = - \left[ {{k^*}\nabla T({\boldsymbol{r}},t) + k\nabla \dot T({\boldsymbol{r}},t)} \right] $$ (6) 通过量纲分析, 引入关系式
$ {k^*} = {k \mathord{\left/ {\vphantom {k \tau }} \right. } \tau } $ , 则GN热传导方程可进一步改写为$$ \tau {\dot{\boldsymbol q}}({\boldsymbol{r}},t) = - \left[ {k\nabla T({\boldsymbol{r}},t) + k\tau \nabla \dot T({\boldsymbol{r}},t)} \right] $$ (7) 在GN模型(5)中也考虑热流率和松弛时间, 可以得到MGT (Moore-Gibson-Thompson)热传导方程[12-13]
$$ {\boldsymbol{q}}({\boldsymbol{r}},t) + \tau {\dot{\boldsymbol q}}({\boldsymbol{r}},t) = - \left[ {{k^*}\nabla \chi ({\boldsymbol{r}},t) + k\nabla T({\boldsymbol{r}},t)} \right] $$ (8) 同样地, 考虑到关系式
$ {k^*} = {k \mathord{\left/ {\vphantom {k \tau }} \right. } \tau } $ , 则有$$ \tau {\dot{\boldsymbol q}}({\boldsymbol{r}},t) + {\tau ^2}{\boldsymbol{\ddot q}}({\boldsymbol{r}},t) = - \left[ {k\nabla T({\boldsymbol{r}},t) + k\tau \nabla \dot T({\boldsymbol{r}},t)} \right] $$ (9) 其他的广义热传导模型包括: Cao等[14]基于相对论理论, 并计及声子质量, 建立的热质热传导方程; Kuang[15]通过引入惯性熵的概念, 建立的惯性熵理论等.
广义热传导的研究直接推动了广义热弹耦合理论的发展. Lord等[16]在CV热传导模型的基础上, 建立了LS广义热弹性模型; Bazarra等[17]在该模型的基础上探讨了材料孔隙率对瞬态响应的影响; Green等[11]建立了无耗散的GN热弹性理论; El-Karamany等[18]考虑了滞后效应, 建立了GN模型Ⅱ型和Ⅲ型的单相滞后、双相滞后模型; Alizadeh Hamidi等[19]探讨了GN模型在Euler-Bernoulli梁的运用; Tzou[20-21]由理论推导建立了DPL热弹耦合理论, 针对瞬态响应和滞后效应开展了大量实验研究; Singh[22]模拟了平面波在横观各向同性双相滞后广义热弹性固体中的传播, 以图形方式显示了双相滞后热弹性和LS广义热弹性情况下的传播角度; Quintanilla[23]发展了MGT热弹性理论, 并考虑了双温下的扩展, 证明了解的多项式衰减; Youssef[24]基于热质热传导方程, 构建了双温度广义热弹性理论; 文献[25-27]建立并发展了GL热弹耦合理论, 讨论了应变率的影响, 证明了解的唯一性定理. 此外, 考虑空间尺度效应和时间记忆效应的热弹耦合理论也已取得较大进展, 如分数阶微积分[28-31]和非局部效应[32-33]的引入等.
综上, 广义热传导模型的形式多样, 对应的热弹耦合理论也很丰富. 因此, 存在的问题有: (1) 如何从连续介质力学视角, 建立起相关理论模型的热力学基础; (2) 各模型之间存在怎样的关联, 如何揭示之.
基于拓展热力学原理, 本文建立了广义热弹耦合理论的新模型, 是现有相关模型的统一形式, 可通过参数赋值实现模型的退化. 此外, 提出了热学“弹性”单元和“黏性”单元, 通过串并联模式实现了上述模型的重构, 从而揭示了已有模型间的关系. 结合数值分析结果, 研究了各模型中松弛时间等对瞬态响应结果的影响规律.
1. 拓展热力学和广义热弹耦合理论
根据连续介质力学, 热力学第一定律为
$$ \rho \dot U = - \nabla \cdot{\boldsymbol{q}} + r + {\boldsymbol{\sigma}} :{\dot{\boldsymbol \varepsilon }} $$ (10) 其中,
$\;\rho $ 为密度,$U$ 为比内能,$r$ 为内热源,${\boldsymbol{\sigma }}$ 为应力张量,${\boldsymbol{\varepsilon }}$ 为应变张量. 热力学第二定律为$$ \rho \dot \eta + \nabla \cdot {{\boldsymbol{J}}^s} - \dfrac{r}{T} \geqslant 0 $$ (11) 其中,
$\eta $ 为比熵,${{\boldsymbol{J}}^s}$ 为熵流矢量. 引入比亥姆霍兹自由能$\psi $ $$ \rho \psi = \rho {{U}} - \rho T\eta $$ (12) 由式(10)~式(12), 有
$$ - \rho \dot \psi + {\boldsymbol{\sigma }}:{\dot{\boldsymbol \varepsilon }} - \rho \dot T\eta - \nabla \cdot {\boldsymbol{q}} + T\nabla \cdot {{\boldsymbol{J}}^s} \geqslant 0 $$ (13) 根据拓展热力学, 现引入高阶热流[34]
$$ -\nabla \cdot {{\boldsymbol{Q}}}^{\left(n\right)} = {\dot{{\boldsymbol{Q}}}}^{(n-1)}\; (n = 2,3) $$ (14) 其中一阶热流
${{\boldsymbol{Q}}^{\left( 1 \right)}}$ 是经典热流${\boldsymbol{q }}={{\boldsymbol{Q}}^{(1)}}$ . 经典热力学中, 熵$\;\rho \eta $ 是比内能和应变的函数. 对于拓展热力学, 假设其也是各阶热流、温度梯度的函数, 则拓展的吉布斯方程为$$ \begin{split} & \rho \dot \eta \left( {U,{\boldsymbol{\varepsilon }},{{\boldsymbol{Q}}^{\left( 1 \right)}},{{\boldsymbol{Q}}^{\left( 2 \right)}},{{\boldsymbol{Q}}^{\left( 3 \right)}},\nabla T} \right)= \\ &\qquad \rho \left( {\dfrac{{\partial \eta }}{{\partial U}}} \right)\dot U + \rho \left( {\dfrac{{\partial \eta }}{{\partial {\boldsymbol{\varepsilon }}}}} \right):{\dot{\boldsymbol \varepsilon }}+ \\ &\qquad \mathop \sum \limits_{n = 1}^3 - {\alpha _n}{{\boldsymbol{Q}}^{\left( n \right)}}\mathop \otimes \limits^n {{\dot {\boldsymbol{Q}}}^{\left( n \right)}} - \varUpsilon \nabla T \cdot \nabla \dot T \end{split} $$ (15) 其中,
$\mathop \otimes \limits^n $ 表示n阶内积, 可得$$ \rho \left( {\dfrac{{\partial \eta }}{{\partial {{\boldsymbol{Q}}^{\left( n \right)}}}}} \right) = - {\alpha _n}{{\boldsymbol{Q}}^{\left( n \right)}}\text{, }\rho \left( {\dfrac{{\partial \eta }}{{\partial \nabla T}}} \right) = - \varUpsilon \nabla T $$ 由式(12)进一步可知, 比亥姆霍兹自由能
$\psi $ 是温度、应变、各阶热流和温度梯度的函数, 故有$$ \begin{split} & \rho \dot \psi \left( {T,{\boldsymbol{\varepsilon }},{{\boldsymbol{Q}}^{\left( 1 \right)}},{{\boldsymbol{Q}}^{\left( 2 \right)}},{{\boldsymbol{Q}}^{\left( 3 \right)}},\nabla T} \right) = \\ &\qquad \rho \left( {\dfrac{{\partial \psi }}{{\partial T}}} \right)\dot T + \rho \left( {\dfrac{{\partial \psi }}{{\partial {\boldsymbol{\varepsilon }}}}} \right):{\dot{\boldsymbol \varepsilon }} + \\ &\qquad \mathop \sum \limits_{n = 1}^3 \rho \left( {\dfrac{{\partial \psi }}{{\partial {{\boldsymbol{Q}}^{\left( n \right)}}}}} \right)\mathop \otimes \limits^n {{\dot {\boldsymbol{Q}}}^{\left( n \right)}} + \rho \left( {\dfrac{{\partial \psi }}{{\partial \nabla T}}} \right) \cdot \nabla \dot T \end{split} $$ (16) 其中
$$ \rho \left( {\dfrac{{\partial \psi }}{{\partial {{\boldsymbol{Q}}^{\left( n \right)}}}}} \right) = \left( {\dfrac{{\partial \psi }}{{\partial \eta }}} \right)\rho \left( {\dfrac{{\partial \eta }}{{\partial {{\boldsymbol{Q}}^{\left( n \right)}}}}} \right) = T{\alpha _n}{{\boldsymbol{Q}}^{\left( n \right)}} $$ $$ \rho \left( {\dfrac{{\partial \psi }}{{\partial \nabla T}}} \right) = \left( {\dfrac{{\partial \psi }}{{\partial \eta }}} \right)\rho \left( {\dfrac{{\partial \eta }}{{\partial \nabla T}}} \right) = T\varUpsilon \nabla T \quad$$ 代入式(16), 有
$$ \begin{split} & \rho \dot \psi \left( {T,{\boldsymbol{\varepsilon }},{{\boldsymbol{Q}}^{\left( 1 \right)}},{{\boldsymbol{Q}}^{\left( 2 \right)}},{{\boldsymbol{Q}}^{\left( 3 \right)}},\nabla T} \right) = \\ &\qquad \rho \left( {\dfrac{{\partial \psi }}{{\partial T}}} \right)\dot T + \rho \left( {\dfrac{{\partial \psi }}{{\partial {\boldsymbol{\varepsilon }}}}} \right):{\dot{\boldsymbol \varepsilon }} + \\ &\qquad \mathop \sum \limits_{n = 1}^3 T{\alpha _n}{{\boldsymbol{Q}}^{\left( n \right)}}\mathop \otimes \limits^n {{\dot {\boldsymbol{Q}}}^{\left( n \right)}} + T\varUpsilon \nabla T \cdot \nabla \dot T \end{split} $$ (17) 由式(13)和式(17)可得
$$ \begin{split} & - \rho \left( {\eta + \dfrac{{\partial \psi }}{{\partial T}}} \right)\dot T + \left( {{\boldsymbol{\sigma }} - \rho \dfrac{{\partial \psi }}{{\partial {\boldsymbol{\varepsilon }}}}} \right):{\dot{\boldsymbol \varepsilon }}- \\ & \left( \mathop \sum \limits_{n = 1}^3 T{\alpha _n}{{\boldsymbol{Q}}^{\left( n \right)}}\mathop \otimes \limits^n {{\dot {\boldsymbol{Q}}}^{\left( n \right)}} + T\varUpsilon \nabla T \cdot \nabla \dot T {\text{ + }}\nabla \cdot {\boldsymbol{q}} - T\nabla \cdot {{\boldsymbol{J}}^s} \right) \geqslant 0 \end{split} $$ (18) 可得
$$ {\boldsymbol{\sigma }} = \rho \dfrac{{\partial \psi }}{{\partial {\boldsymbol{\varepsilon }}}} $$ (19) $$ \rho \eta = - \rho \dfrac{{\partial \psi }}{{\partial T}} $$ (20) $$ \begin{split} &- \Biggr( \mathop \sum \limits_{n = 1}^3 T{\alpha _n}{{\boldsymbol{Q}}^{\left( n \right)}}\mathop \otimes \limits^n {{\dot {\boldsymbol{Q}}}^{\left( n \right)}} + T\varUpsilon \nabla T \cdot \nabla \dot T+ \\ &\qquad \nabla \cdot {\boldsymbol{q}} - T\nabla \cdot {{\boldsymbol{J}}^s} \Biggr) \geqslant 0 \end{split}$$ (21) 经典热力学的熵流为
${{\boldsymbol{J}}^s} = {T^{ - 1}}{\boldsymbol{q}}$ , 通过考虑高阶热流项, 将其拓展为$$ {{\boldsymbol{J}}^s} = {T^{ - 1}}{\boldsymbol{q}} + {\beta _1}{{\boldsymbol{Q}}^{\left( 2 \right)}} \cdot {\boldsymbol{q}} + {\beta _2}{{\boldsymbol{Q}}^{\left( 3 \right)}}:{{\boldsymbol{Q}}^{\left( 2 \right)}} $$ (22) 上式代入式(21), 并由正定性得
$$ - T\varUpsilon \nabla \dot T - {A_1}{T^{ - 1}}{\boldsymbol{q}} = {\chi _0}\nabla T $$ (23) $$ - {A_2}{T^{ - 1}}\nabla T - T{\alpha _1}{\dot{\boldsymbol q}} + T{\beta _1}\nabla \cdot {{\boldsymbol{Q}}^{\left( 2 \right)}} = {\chi _1}{\boldsymbol{q}} $$ (24) $$ T{\beta _1}\nabla {\boldsymbol{q}} - T{\alpha _2}{\dot {\boldsymbol{Q}}^{\left( 2 \right)}} + T{\beta _2}\nabla \cdot {{\boldsymbol{Q}}^{\left( 3 \right)}} = {\chi _2}{{\boldsymbol{Q}}^{\left( 2 \right)}} $$ (25) $$ T{\beta _2}\nabla {{\boldsymbol{Q}}^{\left( 2 \right)}} - T{\alpha _3}{\dot {\boldsymbol{Q}}^{\left( 3 \right)}} = {\chi _3}{{\boldsymbol{Q}}^{\left( 3 \right)}} $$ (26) 其中,
$ {A_1} $ 和$ {A_2} $ 满足$ {A_1} + {A_2} = 1 $ ,${\chi _i}\;\left( {i = 0,1,2,3} \right) \geqslant 0$ . 考虑式(14), 并由式(26)可得$$ \begin{split} & {{\boldsymbol{Q}}^{\left( 2 \right)}} = \dfrac{{T{\beta _1}}}{{{\chi _2}}}\nabla {\boldsymbol{q}} - \dfrac{{T{\alpha _2}}}{{{\chi _2}}}{{\dot {\boldsymbol{Q}}}^{\left( 2 \right)}} +\\ &\qquad \dfrac{{{T^2}{\beta _2}}}{{{\chi _2}}}\left( {\dfrac{{{\beta _2}}}{{{\chi _3}}}\nabla \cdot \nabla {{\boldsymbol{Q}}^{\left( 2 \right)}} - \dfrac{{{\alpha _3}}}{{{\chi _3}}}\nabla \cdot {{\dot {\boldsymbol{Q}}}^{\left( 3 \right)}}} \right) \end{split} $$ (27) 将上式代入式(24), 有
$$ \begin{split} & {\boldsymbol{q}} + \dfrac{{T{\alpha _1}}}{{{\chi _1}}}{\dot{\boldsymbol q}} = - \dfrac{{{A_2}}}{{{\chi _1}}}{T^{ - 1}}\nabla T + \\ &\qquad \dfrac{{T{\beta _1}}}{{{\chi _1}}}\dfrac{{T{\beta _1}}}{{{\chi _2}}}{\nabla ^2}{\boldsymbol{q}} + \dfrac{{T{\beta _1}}}{{{\chi _1}}}\dfrac{{T{\alpha _2}}}{{{\chi _2}}}{\boldsymbol{\ddot q}}- \\ &\qquad \dfrac{{T{\beta _1}}}{{{\chi _1}}}\dfrac{{T{\beta _2}}}{{{\chi _2}}}\dfrac{{T{\beta _2}}}{{{\chi _3}}}{\nabla ^2}{\dot{\boldsymbol q}} - \dfrac{{T{\beta _1}}}{{{\chi _1}}}\dfrac{{T{\beta _2}}}{{{\chi _2}}}\dfrac{{T{\alpha _3}}}{{{\chi _3}}}{\dddot{\boldsymbol q}} \end{split} $$ (28) 考虑到
$ {A_1} + {A_2} = 1 $ , 则式(23)可改写为$$ {A_2}\nabla T = - \dfrac{{T\varUpsilon }}{{{\chi _0}}}\nabla \dot T - \dfrac{{{A_1}{T^{ - 1}}}}{{{\chi _0}}}{\boldsymbol{q}} - {A_1}\nabla T $$ (29) 代入式(28)
$$ \begin{split} & \left( {1 - \dfrac{{{A_1}{T^{ - 2}}}}{{{\chi _0}{\chi _1}}}} \right){\boldsymbol{q}} + \dfrac{{T{\alpha _1}}}{{{\chi _1}}}{\dot{\boldsymbol q}} = \dfrac{{{A_1}{T^{ - 1}}}}{{{\chi _1}}}\nabla T{\text{ + }}\dfrac{\varUpsilon }{{{\chi _0}{\chi _1}}}\nabla \dot T + \\ &\qquad \dfrac{{T{\beta _1}}}{{{\chi _1}}}\dfrac{{T{\beta _1}}}{{{\chi _2}}}{\nabla ^2}{\boldsymbol{q}} + \dfrac{{T{\beta _1}}}{{{\chi _1}}}\dfrac{{T{\alpha _2}}}{{{\chi _2}}}{\boldsymbol{\ddot q}} - \\ &\qquad \dfrac{{T{\beta _1}}}{{{\chi _1}}}\dfrac{{T{\beta _2}}}{{{\chi _2}}}\dfrac{{T{\beta _2}}}{{{\chi _3}}}{\nabla ^2}{\dot{\boldsymbol q}} - \dfrac{{T{\beta _1}}}{{{\chi _1}}}\dfrac{{T{\beta _2}}}{{{\chi _2}}}\dfrac{{T{\alpha _3}}}{{{\chi _3}}}{{\dddot{\boldsymbol{ q}}}} \end{split} $$ (30) 若取
$$ T{\alpha _i} = {\tau _q}{\left( {\dfrac{{{\tau _q}}}{l}} \right)^{2n - 2}}{\chi _1}\left( {1 - \dfrac{{{A_1}{T^{ - 2}}}}{{{\chi _0}{\chi _1}}}} \right),{\text{ }}i = 1,2,3 $$ (31) $$ T{\beta _i} = - {\tau _q}{\left( {\dfrac{{{\tau _q}}}{l}} \right)^{2n - 2}}{\chi _1}\left( {1 - \dfrac{{{A_1}{T^{ - 2}}}}{{{\chi _0}{\chi _1}}}} \right),{\text{ }}i = 1,2 $$ (32) $$ {\chi _i} = {\left( {\dfrac{{{\tau _q}}}{l}} \right)^{2n - 2}}{\chi _1}\left( {1 - \dfrac{{{A_1}{T^{ - 2}}}}{{{\chi _0}{\chi _1}}}} \right),{\text{ }}i = 2,3 $$ (33) 以及
$$ k = - \dfrac{{{A_1}{T^{ - 1}}}}{{{\chi _1}}}{\left( {1 - \dfrac{{{A_1}{T^{ - 2}}}}{{{\chi _0}{\chi _1}}}} \right)^{ - 1}} \text{, } {\tau _T} = \dfrac{\varUpsilon }{{{A_1}{T^{ - 1}}{\chi _0}}} $$ 式(30)可更新为
$$ \begin{split} & {\boldsymbol{q}} + {\tau _q}{\dot{\boldsymbol q}} = - k\nabla T - k{\tau _T}\nabla \dot T + \\ &\qquad {l^2}{\nabla ^2}{\boldsymbol{q}} - \tau _q^2{\boldsymbol{\ddot q}}{\text{ + }}{\tau _q}{l^2}{\nabla ^2}{\dot{\boldsymbol q}} - \tau _q^3{\dddot{\boldsymbol q}} \end{split} $$ (34) 至此, 得到了广义热传导方程. 虽然, 本文中自由能
$ \rho \psi $ 与应变、温度以及各阶热流相关, 但由式(19)和式(20)可看出, 应力和熵本构关系具有经典热弹耦合理论的形式. 广义热传导方程(34)同时考虑了时空尺度效应, 若忽略空间尺度效应, 可得$$ {\boldsymbol{q}} + {\tau _q}{\dot{\boldsymbol q}} + \tau _q^2{\boldsymbol{\ddot q}} + \tau _q^3{{\dddot {\boldsymbol{q}}}} = - k\nabla T - k{\tau _T}\nabla \dot T $$ (35) 显然, 上式可退化为CV模型(2)、DPL模型(4)、GN模型(7)和MGT模型(9).
2. 热学“弹性”单元和“黏性”单元及广义热传导方程的重构
上节中, 基于拓展热力学, 通过考虑热弹耦合得到了应力、熵本构关系和广义热传导方程. 若将其与运动学方程、能量守恒方程和应变−位移关系结合, 则可建立广义热弹耦合模型. 同时发现, 本文得到的广义热传导模型可退化为CV, DPL, GN, MGT等模型.
为了进一步阐释各模型之间的内在关联, 类比于黏弹性力学, 本节提出热学“弹性”单元和热学“黏性”单元, 并通过串并联方式, 实现各广义热传导模型的重构. 为此, 在一维问题中, 考虑GN模型(5)中热位移
$\chi $ 的概念($\dot \chi = T$ ), 热位移梯度记为$\omega = \nabla \chi $ , 则傅里叶定律(1)可写作$$ q = - k\dot \omega $$ (36) GN模型(5)中若忽略
$ k\nabla T $ 项, 可得$$ q = - {k^*}\omega $$ (37) 类比于力学, 称式(37)为热学“弹性”单元, 式(36)为热学“黏性”单元. 其中
$q$ 为热流,${k^*}$ 和$k$ 分别为热学“弹性”单元和热学“黏性”单元的传热系数, 满足关系式$ {k^*} = {k \mathord{\left/ {\vphantom {k \tau }} \right. } \tau } $ . 基于以上单元, 现通过串并联模型, 揭示广义热传导模型间的关联. 基本模型如图1所示.2.1 CV模型
如图1(a)示, 热弹性单元和热黏性单元串联组合, 其结构与黏弹性理论中的Maxwell模型一致. 模型两端的热流相等, 即
$$ q = {q_1} = {q_2} $$ (38) 而模型的总热位移
$\chi $ 为各单元热位移之和, 则$$ \omega = {\omega _1} + {\omega _2} $$ (39) 对式(39)求导
$$ \dot \omega = {\dot \omega _1} + {\dot \omega _2} $$ (40) 其中
$$ {\dot \omega _1} = - \dfrac{{\dot q}}{{{k^{*}}}}, \qquad{\dot \omega _2} = - \dfrac{q}{k} $$ (41) 因此
$$ \dot \omega = - \dfrac{q}{k} - \dfrac{{\dot q}}{{{k^{\text{*}}}}} $$ (42) 整理后可得CV热传导模型
$$ q + \tau \dot q = - k\nabla T $$ (43) 2.2 GN模型
如图1(b)示, 热弹性单元和热黏性单元并联组合, 其结构与黏弹性理论中的Kelvin模型一致. 模型的热位移与各单元的热位移相等, 而总热流为各单元热流之和
$$ \omega = {\omega _1} = {\omega _2}, \qquad q = {q_1} + {q_2} $$ (44) 可得到
$$ \dot q = {\dot q_1} + {\dot q_2} = - k{\ddot \omega _1} - {k^*}{\dot \omega _2} $$ (45) 整理后可得GN热传导模型
$$ \tau \dot q = - k\nabla T - k\tau \nabla \dot T $$ (46) 2.3 DPL模型
如图1(c)所示, 两端总热流和各部分的热流相等, 总热位移为各单元热位移之和
$$ q = {q_1} = {q_2},\qquad \omega = {\omega _1} + {\omega _2} $$ (47) 对第1部分有
$$ {q_1} = - {k_1}{\dot \omega _1} $$ (48) 对第2部分有
$$ {\dot q_2} = - {k_2}{\ddot \omega _2} - k_2^{\text{*}}{\dot \omega _2} $$ (49) 联立式(47)~式(49), 可得
$$ \ddot \omega = - \dfrac{{\dot q}}{{{k_1}}} - \dfrac{{\dot q}}{{{k_2}}} - \dfrac{1}{{{\tau _2}}}\left(\dot \omega + \dfrac{q}{{{k_1}}}\right) $$ (50) 整理可得DPL型热传导方程
$$ q + (1 + \dfrac{{{k_1}}}{{{k_2}}}){\tau _2}\dot q = - {k_1}\dot \omega - {k_1}{\tau _2}\ddot \omega $$ (51) 进一步, 若设定第1和2部分的单元参数一致, 即
$k = {k_1} = {k_2}$ ,${\tau _2}{\text{ = }}\tau $ , 方程改写为$$ q + 2\tau \dot q = - k\nabla T - k\tau \nabla \dot T $$ (52) 值得指出的是, 该方法推导出的DPL模型是特殊形式, 和一般形式(具有两个松弛时间, 热流相松弛时间
${\tau _q}$ 和温度梯度相松弛时间${\tau _T}$ )存在差异, 这是因为推导过程中假定了两个黏性单元参数保持一致.2.4 MGT模型
如图1(d)示, 两端总热流和各部分的热流相等, 总热位移为两部分之和
$$ q = {q_1} = {q_2}, \qquad \omega = {\omega _1} + {\omega _2} $$ (53) 对第1部分有
$$ {q_1} = - k_1^*{\omega _1} $$ (54) 对第2部分有
$$ {q_2} = - k_2^*{\omega _2} - {k_2}{\dot \omega _2} $$ (55) 联立式(53)~式(55)可得
$$ \omega = {\omega _1} + {\omega _2} = - \dfrac{q}{{k_1^{\text{*}}}} - \dfrac{1}{{k_2^*}}(q + {k_2}{\dot \omega _2}) $$ (56) 其中
$$ {\dot \omega _2} = \dot \omega - {\dot \omega _1} = \dot \omega + \dfrac{{\dot q}}{{k_1^*}} $$ (57) 联立式(56)和式(57)有
$$ k_1^*\omega = - q - \dfrac{{k_1^*}}{{k_2^*}}q - k_1^*{\tau _2}\dot \omega - {\tau _2}\dot q $$ (58) 同样地, 若
$k = {k_1} = {k_2}$ ,${\tau _2}{\text{ = }}\tau $ ,$k_1^* = k_2^*$ , 则式(58)整理为$$ 2\tau \dot q + {\tau ^2}\ddot q = - k\nabla T - k\tau \nabla \dot T $$ (59) 虽然式(59)和MGT方程的一般形式仍有差别, 但两式都是由
$ \dot{q},\;\ddot{q},\;\nabla T,\;\nabla \dot{T} $ 等分别作为函数构成的等式, 仅仅是系数的差异. 而这一问题可以通过对常系数的赋值加以解决. 事实上, 区分热传导模型也可以通过对比所含的函数这一方法来完成.2.5 新模型——热场Burgers模型
基于以上的组合思想, 自然可以形成一种新的模型. 如图1 (f)示, 该模型在结构上是由Maxwell模型和Kelvin模型串联组合形成, 与黏弹性理论中的Burgers模型的结构一致. 在前文推导的基础上, 可以直接得出第1和第2部分的方程分别为
$$ {q_1} + {\tau _1}{\dot q_1} = - {k_1}{\dot \omega _1} $$ (60) $$ {q_2} = - k_2^*{\omega _2} - {k_2}{\dot \omega _2} $$ (61) 而对于该模型, 两端总热流和各部分两端热流仍然相等, 总热位移为各部分热位移之和
$$ {q_1} = {q_2} = q, \qquad \omega = {\omega _1} + {\omega _2} $$ (62) 联立式(60)~式(62)并求导可得
$$ \begin{split} & q + \left({\tau _1} + {\tau _2} + \dfrac{{{k_1}}}{{{k_2}}}{\tau _2}\right)\dot q + {\tau _1}{\tau _2}\ddot q = - {k_1}\dot \omega - {k_1}{\tau _2}\ddot \omega \end{split} $$ (63) 假设模型中的弹性单元、黏性单元各自的参数均保持一致, 即
$k = {k_1} = {k_2}$ ,$\tau = {\tau _1} = {\tau _2}$ ,${k^{\text{*}}} = k_1^{\text{*}} = k_2^{\text{*}}$ , 则可将式(63)整理为$$ q + 3\tau \dot q + {\tau ^2}\ddot q = - k\nabla T - k\tau \nabla \dot T $$ (64) 至此, 通过热学“弹性”和“黏性”单元的串并联组合, 重构了CV, GN, DPL, MGT热传导模型, 并通过Burgers模型得到了第1节基于拓展热力学的广义热传导模型. 特别指出的是, 串并联模型给出了广义热传导方程中各松弛时间的关系. 综合2.1~2.5节的推导, 广义热传导各模型的控制方程如表1第2列所示.
表 1 广义热传导模型的变换和统一化Table 1. Transformation and unification of generalized heat conduction modelsModel Heat conduction equation Transformed equation $\bar M$ $\bar N$ CV $ q + \tau \dot q = - k\nabla T $ $ q + \tau sq = - k\nabla T $ $ 1 + \tau s $ 1 GN $ \tau \dot q = - k\nabla T - k\tau \nabla \dot T $ $ \tau sq = - k\nabla T - k\tau s\nabla T $ $ \tau s $ $ 1 + \tau s $ DPL $ q + 2\tau \dot q = - k\nabla T - k\tau \nabla \dot T $ $ q + 2\tau sq = - k\nabla T - k\tau s\nabla T $ $ 1 + 2\tau s $ $ 1 + \tau s $ MGT $ 2\tau \dot q + {\tau ^2}\ddot q = - k\nabla T - k\tau \nabla \dot T $ $ 2\tau sq + {\tau ^2}{s^2}q = - k\nabla T - k\tau s\nabla T $ $ 2\tau s + {\tau ^2}{s^2} $ $ 1 + \tau s $ new model $ q + 3\tau \dot q + {\tau ^2}\ddot q = - k\nabla T - k\tau \nabla \dot T $ $ q + 3\tau sq + {\tau ^2}{s^2}q = - k\nabla T - k\tau s\nabla T $ $ 1 + 3\tau s + {\tau ^2}{s^2} $ $ 1 + \tau s $ 3. 算例计算
考虑一维半无限的材料受边界热冲击和移动热源的作用, 假设无穷远处不受扰动, 且左侧边界无应力, 因此, 边界条件可表达为
$$ \begin{aligned} & \sigma {(x,t)_{x = 0}} = 0, \qquad \theta {(x,t)_{x = 0}} = {T_0}H(t) \\ & u{(x,t)_{x \to \infty }} \to 0, \qquad \theta {(x,t)_{x \to \infty }} \to 0 \end{aligned} $$ 其中,
$\theta = T - {T_0}$ 为温度变化值,${T_0}$ 为起始温度,$H(t)$ 为Heaviside单位阶跃函数,$u$ 为位移. 一维问题下, 运动平衡方程为$$ \dfrac{{\partial \sigma }}{{\partial x}} = \rho \ddot u $$ (65) 考虑了热力耦合的应力本构方程为
$$ \sigma {\text{ = }}\left( {\lambda + 2\mu } \right)\varepsilon - \left( {3\lambda + 2\mu } \right){\alpha _\theta }\theta $$ (66) 其中,
$\lambda $ 和$\mu $ 是拉梅常数,${\alpha _\theta }$ 是热膨胀系数. 几何方程为$$ \varepsilon = \dfrac{{\partial u}}{{\partial x}} $$ (67) 能量守恒方程为
$$ \dfrac{{\partial q}}{{\partial x}} = Q - \rho {T_0}\dot \eta $$ (68) 其中
$Q$ 为热源. 熵本构方程为$$ \rho {T_0}\eta = \left( {3\lambda + 2\mu } \right){\alpha _\theta }{T_0}\varepsilon + \rho {c_E}\theta $$ (69) 为简单起见, 表1中的广义热传导方程统一写为
$$ Mq = - Nk\dfrac{{\partial \theta }}{{\partial x}} $$ (70) 其中, M和N是微分算子, 如表1第2列所示: 对于CV模型, 有
$ M = 1 + \tau {\partial _t} $ ,$ N = 1 $ ; 对本文新建模型, 有$ M = 1 + 3\tau {\partial _t} + {\tau ^2}{\partial _{tt}} $ ,$ N = 1 + \tau {\partial _t} $ . 热传导方程经拉普拉斯变换后的变换方程如表1第3列所示.联立式(65)~式(70), 可得位移和温度的控制方程
$$ (\lambda + 2\mu ){\nabla ^2}u - \gamma \nabla \theta = \rho \ddot u $$ (71) $$ kN{\nabla ^2}\theta = M({T_0}\gamma \nabla \dot u + \rho {c_E}\dot \theta - Q) $$ (72) 其中
$\nabla = \dfrac{\partial }{{\partial x}}$ 为微分算子,$\gamma = (3\lambda + 2\mu ){\alpha _\theta }$ .定义无量纲变量如下
$$ \left\{\begin{split} & (\tilde x,\tilde u) = {n_1}{n_2}(x,u)\\ &(\tilde t,\tilde \tau ) = n_1^2{n_2}(t,\tau )\\ &\tilde \sigma = \dfrac{\sigma }{\mu }\\ & \tilde \theta = \dfrac{\theta }{{{T_0}}}\\ & \tilde Q = \dfrac{Q}{{k{T_0}n_1^2{n_2}}} \end{split}\right. $$ 其中
${n_1} = \sqrt {\dfrac{{\lambda + 2\mu }}{\rho }} $ ,${n_2} = \dfrac{{\rho {c_E}}}{k}$ . 式(66)、式(71)和式(72)无量纲化后为(为了简洁, 省略了波浪线)$$ \qquad\qquad \sigma = {\beta ^2}\nabla u - b\theta $$ (73) $$ \qquad\qquad{\beta ^2}{\nabla ^2}u - b\nabla \theta = {\beta ^2}\ddot u $$ (74) $$ N{\nabla ^2}\theta = M(g\nabla \dot u + \dot \theta - Q) $$ (75) 其中
$$ \left\{ \begin{aligned} & {\beta ^2}{\text{ = }}\dfrac{{\lambda {\text{ + }}2\mu }}{\mu }\\ & b = \dfrac{{(3\lambda + 2\mu ){\alpha _\theta }{T_0}}}{\mu }\\ & g = \dfrac{{(3\lambda + 2\mu ){\alpha _\theta }}}{{k{n_2}}}\end{aligned}\right.$$ 本文运用拉普拉斯变换方法求解这一问题, 为此对式(73)~式(75)做拉普拉斯变换, 可得
$$ \qquad\qquad \bar \sigma = {\beta ^2}\nabla \bar u - b\bar \theta $$ (76) $$ \qquad\qquad{\beta ^2}{\nabla ^2}\bar u - b\nabla \bar \theta = {\beta ^2}{s^2}\bar u $$ (77) $$ \qquad\qquad\bar N{\nabla ^2}\bar \theta = \bar M(gs\nabla \bar u + s\bar \theta - \bar Q) $$ (78) 其中, 变换后的M和N(即
$\bar M$ 和$\bar N$ )如表1第4、5列所示. 联立式(77)和式(78)可得$$ A{\nabla ^4}\bar u - B{\nabla ^2}\bar u + C\bar u = D $$ (79) 其中
$$ \left\{\begin{aligned} &A = \bar N{\beta ^2}\\ &B = \bar N{\beta ^2}{s^2} + \bar Mgbs + \bar M{\beta ^2}s\\ &C = \bar M{\beta ^2}{s^3}\\ &D = - \bar Mb\nabla \bar Q\end{aligned}\right.$$ 至此, 可以通过设定不同的热源函数求解不同条件下的式(79).
3.1 热源函数
$\bar Q = 0$ 式(79)退化为
$$ A{\nabla ^4}\bar u - B{\nabla ^2}\bar u + C\bar u = 0 $$ (80) 该方程的特征方程有4个根, 考虑到一维半无限问题下, 无穷远处不受扰动, 故舍去两个正根. 位移、温度、应力解的形式为
$$ \left.\begin{array}{l} \bar u = \displaystyle\sum\limits_{i = 1}^2 {{{\bar u}_i}} {{\rm{e}}^{ - {k_i}\bar x}} \\ \bar \theta = \displaystyle\sum\limits_{i = 1}^2 {{{\bar \theta }_i}} {{\rm{e}}^{ - {k_i}\bar x}} \\ \bar \sigma = \displaystyle\sum\limits_{i = 1}^2 {{{\bar \sigma }_i}} {{\rm{e}}^{ - {k_i}\bar x}} \end{array}\right\} $$ (81) 其中
${k_i} = \sqrt {\dfrac{{B \pm \sqrt {{B^2} - 4 AC} }}{{2 A}}}\; (i = 1,2)$ .将式(81)代入式(76)~式(78)中, 可解得
$$ \;\;\;\;\;{\bar u_i} = \bar u_i^\theta {\bar \theta _i}, \qquad\bar u_i^\theta = b{k_i}/({\beta ^2}{s^2} - {\beta ^2}k_i^2) $$ (82) $$ {\bar \sigma _i} = \bar \sigma _i^\theta {\bar \theta _i}, \qquad \bar \sigma _i^\theta = - ({\beta ^2}{k_i}\bar u_i^\theta + b) $$ (83) 边界条件经拉普拉斯变换后为
$$ \sum\limits_{i = 1}^2 {{{\bar \sigma }_i}} = 0,\qquad \sum\limits_{i = 1}^2 {{{\bar \theta }_i}} = \dfrac{1}{s} $$ (84) 联立式(82)~式(84), 可得方程组
$$ \left. \begin{gathered} {{\bar \sigma }_1}{\text{ + }}{{\bar \sigma }_2}{\text{ = 0}} \\ {{\bar \theta }_1}{\text{ + }}{{\bar \theta }_2} = \dfrac{1}{s} \\ \end{gathered} \right\}$$ (85) 联立式(83)和式(85), 可解得温度表达式为
$$\left. \begin{array}{l} {\bar \theta _1} = \dfrac{{\bar \sigma _2^\theta }}{{\left( {\bar \sigma _2^\theta - \bar \sigma _1^\theta } \right)s}} \\ {\bar \theta _2} = \dfrac{{\bar \sigma _1^\theta }}{{\left( {\bar \sigma _1^\theta - \bar \sigma _2^\theta } \right)s}} \end{array} \right\}$$ (86) 由式(82)可知,
$ \bar u_i^\theta $ 已知, 可由式(83)进一步解得$ \bar \sigma _i^\theta $ , 又由式(86)可解得${\bar \theta _i}$ , 再回代至式(82)和式(83), 分别解得${\bar u_i}$ 和${\bar \sigma _i}$ . 至此, 求得响应的拉普拉斯域解答.3.2 移动热源
$\bar Q = {{L}}(Q) = \dfrac{{{Q_0}}}{\nu }{{\rm{e}}^{ - \tfrac{s}{\nu }x}}$ 式(79)具体表示为
$$ A{\nabla ^4}\bar u - B{\nabla ^2}\bar u + C\bar u = D{{\rm{e}}^{ - \tfrac{{sx}}{v}}} $$ (87) 其中,
$D = \dfrac{{bs}}{{{v^2}}}\bar M{Q_0}$ . 同样地, 舍去特征方程的两个正根, 方程特解形式为${u_Q}{{\rm{e}}^{ - \tfrac{{sx}}{v}}}$ , 则位移方程为$$ \bar u = {\bar u_1}{{\rm{e}}^{ - {r_1}\bar x}} + {\bar u_2}{{\rm{e}}^{ - {r_2}\bar x}} + {\bar u_Q}{{\rm{e}}^{ - \tfrac{{s\bar x}}{v}}} $$ (88) 其中
$$\left\{ \begin{gathered} \;\;{r_1} = \sqrt {\left(B - \sqrt {{B^2} - 4AC} \right)/(2A)} \\ \;\;{r_2} = \sqrt {\left(B + \sqrt {{B^2} - 4AC} \right)/(2A)} \\ {{\bar u}_Q} = \dfrac{D}{{A{{(s/v)}^4} - B{{(s/v)}^2} + C}} \end{gathered} \right.$$ 拉普拉斯域中的温度、应力解同样可表示为
$$ \bar \theta = {\bar \theta _1}{{\rm{e}}^{ - {r_1}\bar x}} + {\bar \theta _2}{{\rm{e}}^{ - {r_2}\bar x}} + {\bar \theta _Q}{{\rm{e}}^{ - \tfrac{{s\bar x}}{v}}} $$ (89) $$ \bar \sigma = {\bar \sigma _1}{{\rm{e}}^{ - {r_1}\bar x}} + {\bar \sigma _2}{{\rm{e}}^{ - {r_2}\bar x}} + {\bar \sigma _Q}{{\rm{e}}^{ - \tfrac{{s\bar x}}{v}}} $$ (90) 将式(88)和式(89)代入式(76), 可得
$$ {\bar u_i} = \bar u_i^\theta {\bar \theta _i}, \qquad \bar u_i^\theta = b{r_i}/({\beta ^2}{s^2} - {\beta ^2}r_i^2) $$ (91) $$ {\bar u_Q} = \bar u_Q^\theta {\bar \theta _Q},\qquad \bar u_Q^\theta = \dfrac{{bv}}{{{\beta ^2}s({v^2} - 1)}} $$ (92) 联立式(76)、式(88)、式(91)和式(92), 可解得
$$ {\bar \sigma _i} = \bar \sigma _i^\theta {\bar \theta _i}, \qquad \bar \sigma _i^\theta = - ({\beta ^2}{r_i}\bar u_i^\theta + b) $$ (93) $$ {\bar \sigma _Q} = \bar \sigma _Q^\theta {\bar \theta _Q},\qquad \bar \sigma _Q^\theta = - [(s/v){\beta ^2}\bar u_Q^\theta + b] $$ (94) 边界条件经拉普拉斯变换为
$$ \left. \begin{gathered} {{\bar \sigma }_1} + {{\bar \sigma }_2} + {{\bar \sigma }_Q} = 0 \\ {{\bar \theta }_1} + {{\bar \theta }_2} + {{\bar \theta }_Q} = F(s) = \dfrac{1}{s} \end{gathered} \right\} $$ (95) 联立式(94)和式(95), 可得
$$ \left.\begin{gathered} {{\bar \theta }_1} = \dfrac{{\bar \sigma _2^\theta - \bar \sigma _Q^\theta }}{{\bar \sigma _1^\theta - \bar \sigma _2^\theta }}{{\bar \theta }_Q} - \dfrac{{\bar \sigma _2^\theta }}{{\bar \sigma _1^\theta - \bar \sigma _2^\theta }}F(s) \\ {{\bar \theta }_2} = \dfrac{{\bar \sigma _1^\theta - \bar \sigma _Q^\theta }}{{\bar \sigma _2^\theta - \bar \sigma _1^\theta }}{{\bar \theta }_Q} - \dfrac{{\bar \sigma _1^\theta }}{{\bar \sigma _2^\theta - \bar \sigma _1^\theta }}F(s) \\ \end{gathered}\right\} $$ (96) 其中
$$ \left\{ \begin{aligned} &{\bar \theta _Q} = \dfrac{D}{{\bar u_Q^\theta [A{{(s/v)}^4} - B{{(s/v)}^2} + C]}}\\ &\bar \sigma _Q^\theta = - \left(\dfrac{b}{{{v^2} - 1}} + b\right)\end{aligned}\right.$$ 进一步,
${\bar \sigma _1}$ ,${\bar \sigma _2}$ ,${\bar u_1}$ ,${\bar u_2}$ 可由式(89)、式(90)和式(92)联立得到. 将相关参数代入拉普拉斯空间中的应力、温度、位移方程即得解拉普拉斯域下的应力、温度、位移响应.4. 结果与分析
在上一节中对问题进行了描述并给出了数值方法. 本节将对所得数据进行讨论. (1)在无热源的情况下, 将对四个基本模型和新模型的温度、位移、应变响应进行对比; (2)在移动热源下的各模型温度、位移、应变响应对比; (3)对比不同热源移速下
$x = 0.12$ 处DPL模型和新模型温度、位移、应变的时间演化特征. 材料常数如表2示.表 2 材料常数表(铜)Table 2. Material constants (copper)$\lambda/ {{\text{GPa}}} $ $\mu /{{\text{GPa}}} $ ${\alpha _\theta }/\left( {{\rm{m\cdot K^{-1}}}} \right)$ $\rho/\left( {{\text{kg}}\cdot{{\text{m}}^{-3}}} \right)$ ${c_E}/\left( { {\text{J} }\cdot { {\rm{kg^{-1}\cdot K^{-1} } } } } \right)$ ${T_0}/{\text{K}} $ 77.6 38.6 1.78 × 10−5 8945 381 293 通过对
$M$ 和$N$ 的不同赋值可以实现不同模型的转化, 由此可以计算得到无热源状态下各模型$t = 0.06$ 时刻的空间响应如图2所示. 图2(a)为无热源状态下各模型位移响应的空间分布特征. 新模型表现出相比于其他模型更为激进的变化, 而GN模型则更趋保守. 具体为: 新模型会发生更大的峰值位移, 同时位移在空间上递减更快, 至热波边界$x = 0.3$ 处为0, 变形区域更小. 相反地, GN模型产生最小的峰值位移响应, 位移递减的曲线也最为平缓, 变形区域最大. 图2(b)为无热源状态下各模型应力响应的空间分布特征. 各模型均在约产生突变, CV, MGT和新模型在约$x = 0.3$ , 即热波边界处产生突变,$x = 0.3$ 处以外的区域热波尚未到达, 不产生应力响应. 而GN和DPL模型并没有清晰的波前突变, 反映出光滑连续的分布特点. 新模型相比于CV和MGT模型在波前产生的应力更小. 图2(c)为无热源状态下各模型温度的空间分布特征. 与应力特征类似, CV, MGT和新模型在约$x = 0.3$ , 即波前处的温度并不连续, 热边界条件对$x = 0.3$ 以外的区域不产生温度影响. 新模型在波前的温度响应更小, 而GN和DPL模型的温度分布光滑连续, 逐步归于无响应. 综合以上信息, 新模型在热波到达时, 位移与应力响应峰值更大, 但响应区域更小. 新模型和CV, MGT模型均有明显的波前突变, 在热波未达区不产生响应, 但新模型的波前突变更小. GN和DPL模型则呈现光滑连续的变化, 响应区域更大.![]()
图 2 无热源下各模型位移、应力、温度响应Figure 2. Displacement, stress and temperature responses of each model without heat source在移动热源移速
$v = 2$ 的作用下, 各模型位移、应力、温度的空间分布如图3所示. 图3(a)为移动热源下各模型位移响应的空间分布, 相比于无热源状态, 各模型的位移均有较大增加, 尤其在$x = 0.12$ 处即移动热源作用处, 位移达到最大值. 新模型仍然具有各模型中最大的位移, 同时位移在空间上的递减也更快, 变形区域最小. 图3(b)为移动热源下各模型应力响应的空间分布, 热源作用下的应力峰值相比于无热源状态明显增大. 无热源状态下的应力峰值约为0.05, 移动热源下的峰值应力已跃升至约0.125, 且应力峰值出现在移动热源作用处. 新模型表现出最大的峰值应力, 而GN模型应对热源作用产生的应力变化则最小. 图3(c)为移动热源下各模型温度的空间分布. 热源作用下必然会引起各模型温度峰值的上升, 其中尤为明显的是CV, MGT和新模型. 其温度曲线也在, 即波前处发生突变, 回复为常温状态. GN和DPL模型的温度曲线是光滑连续的, 温度沿空间逐步降至常温, 故温度变化的区域更广. 由于移动热源的作用, 导致各模型无论在波前、波后均产生更大的位移、应力和温度响应, 同时在移动热源作用点$x = 0.12$ 处相比无热源情况明显升高. 同无热源状态的响应特点相似的是: CV, MGT和新模型在波前的响应也存在突变, 新模型的突变更小. 与无热源状态的响应特点不同的是: CV, MGT和新模型在热源作用点$x = 0.12$ 处也产生了突变, GN和DPL模型无论在热源作用点和波前处, 均呈现光滑连续的变化特点.![]()
图 3 移动热源下各模型的位移、应力、温度响应Figure 3. Displacement, stress and temperature responses of each model under moving heat sourceDPL模型和新模型在
$x = 0.12$ 处受$v = 2$ 和$v = 4$ 的移动热源作用, 移动热源到达该点的时间分别为$t = 0.06$ 和$t = 0.03$ , 其时间演化特征如图4所示. 图4 (a)为该作用环境下位移随时间演化的规律. 两模型均在经历过移动热源的加热后产生较大位移, 其中$v = 2$ 热源作用下的位移响应更大.![]()
图 4 移动热源下DPL模型和新模型位移、应力、温度的时间演化特征Figure 4. Evolution laws of displacement, stress and temperature of DPL model and new model under moving heat source新模型的位移响应较DPL模型峰值更高. 图4(b)为DPL和新模型的应力时间演化曲线. 由于移动热源的加热作用, DPL和新模型均在热源到达时取得极值. 新模型具有比DPL模型更大的峰值应力. 在移动热源离开该点后, 新模型的应力响应迅速回复至与DPL模型相近的水平. 图4(c)为该作用环境下温度的时间演化曲线. 与图4(b)类似, DPL和新模型经历热源作用后, 温度均有明显上升. 在温度上升阶段, 新模型表现出比DPL模型更大更快的变化, 拥有更高的峰值温度. 新模型的温度回复也较快, 与应力时间演化曲线呈现的特点一致, 新模型在受热源作用前后产生了突变. 离开移动热源作用后, 新模型的温度响应迅速回复至与DPL模型相近的数值, 之后两模型以较为近似的曲线缓慢下降. 热源移速的变化, 使得各项数据的达峰时间随移速的增大而提前. 在时间演化特征曲线上, 可以明显看到新模型在热源作用时产生的响应突变, 而DPL模型则呈现为光滑连续的变化, 这一点和空间响应呈现的特征是一致的.
5. 总 结
超快热冲击等极端热作用环境下, 以傅里叶定律为代表的经典热传导理论不再适用, 为发展非平衡态热弹耦合新理论提出迫切需求. 学者们通过拓展经典热传导方程, 已提出多种广义热传导模型, 并通过热力耦合建立了广义热弹性理论. 为了阐清已有各模型间的内在关联, 并进一步发展极端热弹耦合理论, 本文基于拓展热力学原理, 建立了考虑高阶项的广义热传导理论, 其可退化为已有相关理论. 同时, 类比于黏弹性力学理论, 本文抽象出热学“弹性”单元和“黏性”单元模型, 基于串并联方式实现了已有各广义热传导方程的重构, 并通过Burgers模型得到了本文新建立的热传导方程, 又利用拓展热力学原理, 在原理层面上推导支撑了新模型的理论基础. 广义热传导方程的重构过程也揭示了部分模型中, 热流率相与温度梯度率相的松弛时间的比例关系. 如对DPL模型而言, 热流率
$\dot q$ 相与温度梯度率$\nabla \dot T$ 相的松弛时间的比例为2:1. 数值结果表明: 无移动热源作用时, 新模型得到的位移、应力、温度响应峰值更大, 影响范围较小; 在移动热源作用下, 新模型预测的响应峰值明显升高, 且高于其他模型. 对比DPL模型, 新模型在拥有DPL模型双相滞后特点的同时, 解决了热波速度无限大的悖论, 在响应曲线上反映出清晰的波前突变. 本研究将推动远离平衡态极端热弹耦合理论的发展, 并为相关问题的求解提供了分析方法. -
![]()
图 2 无热源下各模型位移、应力、温度响应
Figure 2. Displacement, stress and temperature responses of each model without heat source
![]()
图 3 移动热源下各模型的位移、应力、温度响应
Figure 3. Displacement, stress and temperature responses of each model under moving heat source
![]()
图 4 移动热源下DPL模型和新模型位移、应力、温度的时间演化特征
Figure 4. Evolution laws of displacement, stress and temperature of DPL model and new model under moving heat source
表 1 广义热传导模型的变换和统一化
Table 1 Transformation and unification of generalized heat conduction models
Model Heat conduction equation Transformed equation $\bar M$ $\bar N$ CV $ q + \tau \dot q = - k\nabla T $ $ q + \tau sq = - k\nabla T $ $ 1 + \tau s $ 1 GN $ \tau \dot q = - k\nabla T - k\tau \nabla \dot T $ $ \tau sq = - k\nabla T - k\tau s\nabla T $ $ \tau s $ $ 1 + \tau s $ DPL $ q + 2\tau \dot q = - k\nabla T - k\tau \nabla \dot T $ $ q + 2\tau sq = - k\nabla T - k\tau s\nabla T $ $ 1 + 2\tau s $ $ 1 + \tau s $ MGT $ 2\tau \dot q + {\tau ^2}\ddot q = - k\nabla T - k\tau \nabla \dot T $ $ 2\tau sq + {\tau ^2}{s^2}q = - k\nabla T - k\tau s\nabla T $ $ 2\tau s + {\tau ^2}{s^2} $ $ 1 + \tau s $ new model $ q + 3\tau \dot q + {\tau ^2}\ddot q = - k\nabla T - k\tau \nabla \dot T $ $ q + 3\tau sq + {\tau ^2}{s^2}q = - k\nabla T - k\tau s\nabla T $ $ 1 + 3\tau s + {\tau ^2}{s^2} $ $ 1 + \tau s $ 
下载: 导出CSV
表 2 材料常数表(铜)
Table 2 Material constants (copper)
$\lambda/ {{\text{GPa}}} $ $\mu /{{\text{GPa}}} $ ${\alpha _\theta }/\left( {{\rm{m\cdot K^{-1}}}} \right)$ $\rho/\left( {{\text{kg}}\cdot{{\text{m}}^{-3}}} \right)$ ${c_E}/\left( { {\text{J} }\cdot { {\rm{kg^{-1}\cdot K^{-1} } } } } \right)$ ${T_0}/{\text{K}} $ 77.6 38.6 1.78 × 10−5 8945 381 293
下载: 导出CSV
-
[1] Zhmakin AI. Heat conduction beyond the Fourier law. Technical Physics, 2021, 66(1): 1-22 doi: 10.1134/S1063784221010242
[2] 李吉伟, 何天虎. 考虑应变率的广义压电热弹理论及其应用. 力学学报, 2020, 52(5): 1267-1276 (Li Jiwei, He Tianhu. A generalized piezoelectric-thermoelastic theory with strain rate and its application. Chinese Journal of Theoreticaland and Applied Mechanics, 2020, 52(5): 1267-1276 (in Chinese) doi: 10.6052/0459-1879-20-120 [3] 李妍, 何天虎, 田晓耕. 超短激光脉冲加热薄板的广义热弹扩散问题. 力学学报, 2020, 52(5): 1255-1266 (Li Yan, He Tianhu, Tian Xiaogeng. A generalized thermoelastic diffusion problem of thin plate heated by the ultrashort laser pulses. Chinese Journal of Theoretical and Applied Mechanics, 2020, 52(5): 1255-1266 (in Chinese) doi: 10.6052/0459-1879-20-118 [4] 张培, 何天虎. 考虑非局部效应和记忆依赖微分的广义热弹问题. 力学学报, 2018, 50(3): 508-516 (Zhang Pei, He Tianhu. A generalized thermoelastic problem with nonlocal effect and memory-dependent derivative. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(3): 508-516 (in Chinese) doi: 10.6052/0459-1879-18-079 [5] Li SN, Cao BY. Anomalous heat diffusion from fractional Fokker–Planck equation. Applied Mathematics Letters, 2020, 99: 105992 doi: 10.1016/j.aml.2019.07.023
[6] Yu YJ, Deng ZC. New insights on microscale transient thermoelastic responses for metals with electron-lattice coupling mechanism. European Journal of Mechanics-A/Solids, 2020, 80: 103887 doi: 10.1016/j.euromechsol.2019.103887
[7] Cattaneo C. A form of heat equation which eliminates the paradox of instantaneous propagation. Compete Rendus, 1958, 247: 431-433
[8] Vernotte P. Paradoxes in the continuous theory of the heat conduction. Compte Rendus, 1958, 246: 3154-3155
[9] Tzou DY. The generalized lagging response in small-scale and high-rate heating. International Journal of Heat and Mass Transfer, 1995, 38(17): 3231-3240
[10] Green AE, Naghdi PM. On undamped heat waves in an elastic solid. Journal of Thermal Stresses, 1992, 15(2): 253-264
[11] Green AE, Naghdi PM. Thermoelasticity without energy dissipation. Journal of Elasticity, 1993, 31(3): 189-208
[12] Quintanilla RXFN. Instability and non-existence in the nonlinear theory of thermoelasticity without energy dissipation. Continuum Mechanics and Thermodynamics, 2001, 13(2): 121-129
[13] Quintanilla R. Moore−Gibson−Thompson thermoelasticity. Mathematics and Mechanics of Solids, 2019, 24(12): 4020-4031 doi: 10.1177/1081286519862007
[14] Cao BY, Guo ZY. Equation of motion of a phonon gas and non-Fourier heat conduction. Journal of Applied Physics, 2007, 102(5): 053503 doi: 10.1063/1.2775215
[15] Kuang Z. Variational principles for generalized dynamical theory of thermopiezoelectricity. Acta Mechanica, 2009, 203(1-2): 1-11 doi: 10.1007/s00707-008-0039-1
[16] Lord HW, Shulman Y. A generalized dynamical theory of thermoelasticity. Journal of the Mechanics and Physics of Solids, 1967, 15(5): 299-309 doi: 10.1016/0022-5096(67)90024-5
[17] Bazarra N, Fernández JR, Quintanilla R. Lord−Shulman thermoelasticity with microtemperatures. Applied Mathematics & Optimization, 2021, 84(2): 1667-1685
[18] El-Karamany AS, Ezzat MA. On the phase−lag Green–Naghdi thermoelasticity theories. Applied Mathematical Modelling, 2016, 40(9-10): 5643-5659
[19] Alizadeh Hamidi B, Hosseini SA, Hassannejad R, et al. An exact solution on gold microbeam with thermoelastic damping via generalized Green-Naghdi and modified couple stress theories. Journal of Thermal Stresses, 2020, 43(2): 157-174 doi: 10.1080/01495739.2019.1666694
[20] Tzou DY. Experimental support for the lagging behavior in heat propagation. Journal of Thermophysics and Heat Transfer, 1995, 9(4): 686-693
[21] Tzou DY. Macro-to Microscale Heat Transfer: The Lagging Behavior. John Wiley & Sons, 2014: 388-391
[22] Singh B. Wave propagation in dual-phase-lag anisotropic thermoelasticity. Continuum Mechanics and Thermodynamics, 2013, 25(5): 675-683 doi: 10.1007/s00161-012-0261-x
[23] Quintanilla R. Moore-Gibson-Thompson thermoelasticity with two temperatures. Applications in Engineering Science, 2020, 1: 100006 doi: 10.1016/j.apples.2020.100006
[24] Youssef HM. A novel theory of generalized thermoelasticity based on thermomass motion and two-temperature heat conduction. Journal of Thermal Stresses, 2021, 44(2): 133-148 doi: 10.1080/01495739.2020.1838247
[25] Green AE, Lindsay KA. Thermoelasticity. Journal of Elasticity, 1972, 2(1): 1-7 doi: 10.1007/BF00045689
[26] Yu YJ, Xue ZN, Tian XG. A modified Green–Lindsay thermoelasticity with strain rate to eliminate the discontinuity. Meccanica, 2018, 53(10): 2543-2554 doi: 10.1007/s11012-018-0843-1
[27] Marin M, Craciun EM, Pop N. Some results in Green–Lindsay thermoelasticity of bodies with dipolar structure. Mathematics, 2020, 8(4): 497 doi: 10.3390/math8040497
[28] Yu YJ, Zhao LJ. Fractional thermoelasticity revisited with new definitions of fractional derivative. European Journal of Mechanics-A/Solids, 2020, 84: 104043 doi: 10.1016/j.euromechsol.2020.104043
[29] Yu YJ, Deng ZC. Fractional order theory of Cattaneo-type thermoelasticity using new fractional derivatives. Applied Mathematical Modelling, 2020, 87: 731-751 doi: 10.1016/j.apm.2020.06.023
[30] Yu YJ, Deng ZC. Fractional order thermoelasticity for piezoelectric materials. Fractals, 2021, 29(4): 2150082 doi: 10.1142/S0218348X21500821
[31] Yu YJ, Li SS, Deng ZC. Unified theory of 2n + 1 order size-dependent beams: Mathematical difficulty for functionally graded size-effect parameters solved. Applied Mathematical Modelling, 2020, 79: 314-340 doi: 10.1016/j.apm.2019.10.038
[32] Abouelregal AE. A novel model of nonlocal thermoelasticity with time derivatives of higher order. Mathematical Methods in the Applied Sciences, 2020, 43(11): 6746-6760 doi: 10.1002/mma.6416
[33] Yu YJ, Tian XG, Xiong QL. Nonlocal thermoelasticity based on nonlocal heat conduction and nonlocal elasticity. European Journal of Mechanics-A/Solids, 2016, 60: 238-253 doi: 10.1016/j.euromechsol.2016.08.004
[34] Machrafi H, Lebon G. General constitutive equations of heat transport at small length scales and high frequencies with extension to mass and electrical charge transport. Applied Mathematics Letters, 2016, 52: 30-37
-
期刊类型引用(3)
1. 郭颖，杜翠，马建军. 基于分数阶理论的两端固定杆热弹动力响应. 河南科技大学学报(自然科学版). 2024(01): 68-78+8-9 . 
百度学术
2. 尉亚军，吴华，徐成辉，谷旭东，邓子辰. 连续介质力学的类比思维教学法——黏弹性力学vs.广义热传导1). 力学与实践. 2024(03): 625-630 . 百度学术
3. Hua WU，Xinyi LI，Yajun YU，Zichen DENG. A novel electron-phonon coupling thermoelasticity with Burgers electronic heat transfer. Applied Mathematics and Mechanics(English Edition). 2023(11): 1927-1940 . 必应学术
其他类型引用(2)
计量
- 文章访问数: 1513
- HTML全文浏览量: 415
- PDF下载量: 198
- 被引次数: 5
