机场道面的温度场反问题.pdfVIP

机场道面的温度场反问题 谢永亮方维凤郭跃明 空军后勤学院五系徐州221000 摘要 本文介绍了用有限元反问题的方法在机场道面温度场中的应角，根据道面上两点的测量温度 记录，经有限元离散建立一个最小误差泛函，再用动态规划识别表面热流，从而确定整个温度场． 一、引言 我国北方机场的气温变化大，跑道温度应力影响很大，空军在张家口机场进行这方面的实验 研究工作，但尚未得到完整的理论体系．我们尝试将有限元反问题的方法应用到道面板的温度场 中，实测道面板某些点的温度，识别出边界热流q(t)，由此计算得到板的温度场T(x，t)，为应力场的 计算打下基础． 传统的热传导问题是已知边界热流和其它边初值条件，根据热传导方程来确定温度场．另外 一类问题是：知道某些点的温度和部分边界条件，根据热传导方程来确定另一些边界条件．这属 于热传导反问题·本文提出用反问题的方法来识别边界热流，而确定道面板的温度场．这方法在 工程中有着广泛的应用前景． 二、反问题模型 道面板的热传导应满足于下面的导热微分方／tl： V7(toyT)～Q 觚 pc鲁 (1) 初始条件： T=瓦 whent=0 (2a) 边界条件 On I1 (2b) + — JI 卜wI卜n+q 呱T k 0 on r2 (2c) 上式中V是梯度算子，这里物理量k、p、c仅与位置有关，而与时阈无关，Q(x，t)是热源项。 方程(1)和(2)可用有限元来求解． 根据变分原理，方程(1)的解T(x，t)等价于下面泛函求极值函数。 635 。 将区域Q进行有限元离散，泛函为： I=El(。’ (3) 这里，e’为单元上的泛函，即有 设在单元e上，t时刻的温度T(x，t)可近似地展开为： T。(Kt)=∑N，(x)玉(t)=fN)1{T)。 (5’ 式中N为形函数，{N)为形函数矩阵． l●■■，●■■●●■■， 将式(5)代入式(4)，泛函I就可以用单元上的节点温度表示，即泛函是节点温度的多元函数． I；l(Tl，T2，…，Tn) (6) 泛函，取极值，我们得到： 旦：o 仉 由式(3)、(4)和(5)，上式可写成 ∑厶c科}阿(k甲{N}7{Tj))+Q—pc云({N}-{T／)]dD