- 1、本文档共12页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
6. 稳态热传导问题的有限元法
本章的内容如下:
6.1热传导方程与换热边界
6.2稳态温度场分析的一般有限元列式
6.3三角形单元的有限元列式
6.4温度场分析举例
6.1热传导方程与换热边界
在分析工程问题时,经常要了解工件内部的温度分布情况,例如发动机的工作温度、金属工件在热处理过程中的温度变化、流体温度分布等。物体内部的温度分布取决于物体内部的热量交换,以及物体与外部介质之间的热量交换,一般认为是与时间相关的。物体内部的热交换采用以下的热传导方程(Fourier方程)来描述,
(6-1)
式中为密度,kg/m3; 为比热容,;为导热系数,;T为温度,℃;t为时间,s;为内热源密度,w/m3。
对于各向同性材料,不同方向上的导热系数相同,热传导方程可写为以下形式,
(6-2)
除了热传导方程,计算物体内部的温度分布,还需要指定初始条件和边界条件。初始条件是指物体最初的温度分布情况,
(6-3)
边界条件是指物体外表面与周围环境的热交换情况。在传热学中一般把边界条件分为三类。
给定物体边界上的温度,称为第一类边界条件。
物体表面上的温度或温度函数为已知,
或 (6-4)
给定物体边界上的热量输入或输出,称为第二类边界条件。
已知物体表面上热流密度,
或 (6-5)
给定对流换热条件,称为第三类边界条件。
物体与其相接触的流体介质之间的对流换热系数和介质的温度为已知。
(6-6)
其中h为换热系数,W/(m2 K);是物体表面的温度;是介质温度。
如果边界上的换热条件不随时间变化,物体内部的热源也不随时间变化,在经过一定时间的热交换后,物体内各点温度也将不随时间变化,即
这类问题称为稳态(Steady state)热传导问题。稳态热传导问题并不是温度场不随时间的变化,而是指温度分布稳定后的状态,我们不关心物体内部的温度场如何从初始状态过渡到最后的稳定温度场。随时间变化的瞬态(Transient)热传导方程就退化为稳态热传导方程,三维问题的稳态热传导方程为,
(6-7)
对于各向同性的材料,可以得到以下的方程,称为Poisson方程,
(6-8)
考虑物体不包含内热源的情况,各向同性材料中的温度场满足Laplace方程,
(6-9)
在分析稳态热传导问题时,不需要考虑物体的初始温度分布对最后的稳定温度场的影响,因此不必考虑温度场的初始条件,而只需考虑换热边界条件。计算稳态温度场实际上是求解偏微分方程的边值问题。温度场是标量场,将物体离散成有限单元后,每个单元结点上只有一个温度未知数,比弹性力学问题要简单。进行温度场计算时有限单元的形函数与弹性力学问题计算时的完全一致,单元内部的温度分布用单元的形函数,由单元结点上的温度来确定。由于实际工程问题中的换热边界条件比较复杂,在许多场合下也很难进行测量,如何定义正确的换热边界条件是温度场计算的一个难点。
6.2稳态温度场分析的一般有限元列式
在前面我们已经介绍了有限元方法可以用来分析场问题,稳态温度场计算是一个典型的场问题。我们可以采用虚功方程建立弹性力学问题分析的有限元格式,推导出的单元刚度矩阵有明确的力学含义。在这里,介绍如何用加权余量法(Weighted Residual Method)建立稳态温度场分析的有限元列式。
微分方程的边值问题,可以一般地表示为未知函数u满足微分方程组,
(在域内) (6-10)
未知函数u还满足边界条件,
(在边界上) (6-11)
如果未知函数u是上述边值问题的精确解,则在域中的任一点上u都满足微分方程(6-10),在边界的任一点上都满足边界条件(6-11)。对于复杂的工程问题,这样的精确解往往很难找到,需要设法寻找近似解。所选取的近似解是一族带有待定参数的已知函数,一般表示为
(6-12)
其中为待定系数,为已知函数,被称为试探函数。试探函数要取自完全的函数序列,是线性独立的。由于试探函数是完全的函数序列,任一函数都可以用这个序列来表示。
采用这种形式的近似解不能精确地满足微分方程和边界条件,所产生的误差就称为余量。
微分方程(6-10)的余量为,
(6-13)
边界条件(6-11)的余量为,
(6-14)
选择一族已知的函数,使余量的加权积分为零,强迫近似解所产生的余量在某种平均意义上等于零,
(6-15)
称为权函数,通过公式(6-15)可以选择待定的参数。
这种采用使余量的加权积分为零来求得微分方程近似解的方法称为加权余量法。对权函数的不同选择就得到了不同的加权余量法,常用的方法包括配点法、子域法、最小二乘法、力矩法和伽辽金法(Galerkin method)。在很多情况下,采用Galer
您可能关注的文档
- 水中本底成分对催化臭氧氧化降解富里酸影响.pdf
- 水中隐孢子虫检测与活性鉴定.pdf
- 水资源—多目标粒子群算法在梯级水库联合防洪调度中应用研究.pdf
- 水资源价值定量探讨.pdf
- 水资源转移利用利益补偿测算:模型构建与应用.pdf
- 顺酐装置溶剂吸收塔方案模拟及分析.pdf
- 瞬时无功功率补偿和全频谱滤波.pdf
- 瞬态热线法导热系数测量数值模拟.pdf
- 说说诗歌中意象和意境.doc
- 丝瓜瓤固定简青霉吸附废水中Pb_2_和Cu_2_机理.pdf
- 中国国家标准 GB 42125.11-2025测量、控制和实验室用电气设备的安全要求 第11部分:用于处理医用材料的灭菌器和清洗消毒器的特殊要求.pdf
- 《GB 42125.11-2025测量、控制和实验室用电气设备的安全要求 第11部分:用于处理医用材料的灭菌器和清洗消毒器的特殊要求》.pdf
- 《GB/T 15237-2025术语工作及术语科学 词汇》.pdf
- GB/T 15237-2025术语工作及术语科学 词汇.pdf
- GB/T 22454-2025企业建模与体系结构 企业建模构件.pdf
- 中国国家标准 GB/T 15237-2025术语工作及术语科学 词汇.pdf
- 中国国家标准 GB/T 22454-2025企业建模与体系结构 企业建模构件.pdf
- 《GB/T 22454-2025企业建模与体系结构 企业建模构件》.pdf
- 中国国家标准 GB/T 18910.64-2025液晶显示器件 第6-4部分:测试方法 带动态背光的液晶显示模块.pdf
- 《GB/T 18910.64-2025液晶显示器件 第6-4部分:测试方法 带动态背光的液晶显示模块》.pdf
文档评论(0)