热传导 流体流动 薄板弯曲.ppt

第七章　热传导问题的有限元法 一、问题的提出 热传导方程 热量传递有三种方式：传导、对流和辐射。热量通过固体传递的方式称为传导；由于流体流动而引起的热量传递成为对流，通过电磁波进行热量传递成为辐射。 在热传导过程中，结构中每一点都有一个温度值，他们构成具有物体形状的温度场，场变量就是温度T。温度可以使随着时间变化的，称为瞬态温度场，这时候T=T（x,y,z,t），也可能与时间无关，成为稳态温度场， T=T（x,y,z） 对均质物体，物体内温度不随时间而变化，温度分布函数，应满足拉普拉斯方程 热边界条件： 总结: 1.将描述物理现象的微分方程及其边界条件通过变分原理转化为某一泛函求极值问题.再通过一定的推导得到有限元计算格式.这一方法可以应用到其它场问题的有限元求解中.重要我们找到与微分方程及其边界条件等价的泛函,就可以用变分原理. 2.泛函求极值问题与微分方程边值问题比较,有边界条件自动满足,导数阶次降低,推导有限元格式方便等优点. 3.推导方法也可用于不同单元节点数的二维或三维等参数元计算格式,也可用于其它场问题的有限元计算格式推导. 习题: 1.在稳定温度场有限元计算格式推导过程中,泛函求极值问题是如何转换成多元函数求极值问题的? 2.简述利用变分原理推导稳定温度场有限元计算的主要思路和基本过程. 第八章 流体流动问题的有限元方法 一 问题的提出 1. 采用变分原理可以推导某些偏微分方程所描述的非结构问题的有限元格式.但是泛函求极值问题和微分方程求解问题只是在一定的条件下等价,既要求求解函数连续且可导; 2. 给出一个泛函形式,通常可以找到相对应的微分方程;反之,若给出一个微分方程,并不一定能够找到与之对应泛函极值问题的具体形式. 要用有限元方法求解某一场问题的微分方程,如果找不到与之相对应的泛函,则无法用变分原理来推导相应的有限元计算格式. 第八章 流体流动问题的有限元方法 很多问题都不能用变分法 流体流动问题 高速车辆运行时的空气阻力问题 列车进入隧道时的压力波动问题 两高速运动的列车会车时的压力波动问题 桥梁的风致震动问题 汽车高速行驶时的气动稳定性问题 方程为非线性偏微分方程 第八章 流体流动问题的有限元方法 求解方法:加权余量法 思想: 某物理问题的控制微分方程及边界条件分别为 在域 内 在域 边界S上 为待求函数. 选择一试探函数 待求常数 试探函数项 第八章 流体流动问题的有限元方法 将试探函数代入控制方程和边界条件，在域内和边界上会产生误差 在域 内 在域 边界S上 为余量（残数、残差、残值）。 基本思想：在域内或边界上寻找n个线性无关和函数 使余量R和Rb在加权平均的意义上等于零 第八章 流体流动问题的有限元方法 尽管试探函数本身不能满足控制方程和边界条件。但是当其雨量与许多线性无关的权函数相乘并积分十，这个余量在总体上接近于零，也就是说试探函数在积分意义上满足微分方程式和边界条件。当n足够大时，试探函数就趋近于真实解。 基本方法：1.最小二乘加权余量法 首先找到一个已经满足边界条件式的试探函数 代入控制方程会产生余量 第八章 流体流动问题的有限元方法 基本方法：1.最小二乘加权余量法 代入控制方程会产生余量 简单地说，最小二乘的思想就是要使得观测点和估计点的距离的平方和达到最小.这里的“二乘”指的是用平方来度量观测点与估计点的远近（在古汉语中“平方”称为“二乘”），“最小”指的是参数的估计值要保证各个观测点与估计点的距离的平方和达到最小. 希望余量在最小二乘的意义下为最小. 第八章 流体流动问题的有限元方法 方差为: 由此可知道: 方差最小的条件: 通过求解线性方程组,可得到待定系数. 第八章 流体流动问题的有限元方法 控制方程 最小二乘加权余量法: 有一个待定常数 代入控制方程 对于复杂工程问题,不需要求精确解 第八章 流体流动问题的有限元方法 控制方程 两个未知数 代入控制方程 第八章 流体流动问题的有限元方法 解得: 基本方法：2.伽辽金加权余量法 如果在 即称为伽辽金加权余量法. 权函数选用试探函数,即 许多物理问题控制微分方程的有限元求解,都采用伽辽金加权余量法. 第八章 流体流动问题的有限元方法 控制方程 伽辽金加权余量法: 有一个待定常数 代入控制方程 第八章 流体流动问题的有限元方法 两个未知数 代入控制方程 伽辽金加权余量法: 控制方程 根据边界条件,设 边界条件 取前2项 为形函数. 代入 代入可得到 第9章 薄板的弯曲有限元分析 第1节 薄板弯曲理论基础 一、薄板弯曲问题应满足的条件 1.几何条件 要求结构属于薄板.工程中常将厚度尺寸小于其他两个方向尺寸的结