版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、7 流体流动问题的有限元法 本章应用加权余量法,将求解域的微分方程,转化为积分表达式,然后通过求积分的极值,找到原问题的解。 7-1 问题的提出 如果可以找到与微分方程对应的泛函积分表达式,则可以通过变分原理,建立有限元格式。但是以下两种情况没有泛函: (1)解函数在求解域不连续或不可导; (2)无法找到与微分方程相对应的泛函表达式。,描述稳态不可压缩流体流动的微分方程就没有与之对应的泛函表达式。而流体流动问题却是工程中经常遇到的问题 (1)车辆高速运行时的气动稳定性; (2)两列高速运行的列车会车时的压力波动; (3)列车进入隧道时的压力波动; (4)建筑物的风载荷; (5)室内的通风与空调
2、; (6)桥梁的风致振动; (7)船舶的运行阻力; (8)飞机的升力、阻力。,不论是那种原因,如果找不到与微分方程对应的泛函表达式,那么就无法利用变分原理建立有限元的计算格式。这时我们只有寻求另外的途径。 这个途径就是:加权余量法。 二 加权余量法 加权余量法的基本思想:通过使试探函数与真值的加权误差在求解域内的总和为零,以求得满足微分方程的近似解。 设某物理问题的控制微分方程及其边界条件分别为,为待求函数。如果无法或不易直接求解,可选一个试探函数 式中 ci待定常数; i试探函数项。 将试探函数带入控制微分方程及其边界条件,一般来讲不可能正好满足方程,在域内和边界S上会产生误差,即 式中R和
3、Rb称为余量(或残数,残差,残值)。,加权余量法的基本思想:在域内 和/或 边界S上寻找n个线性无关的函数Wi(i=1,2,n),使余量R和Rb在加权求和的意义上等于零,即 这里Wi称为权函数。 加权余量法所假设的试探函数并不能满足微分方程及其边界条件,但是当加权的试探函数与真值的误差(余量)在求解域上积分为零时,那么试探函数就在总体上满足微分方程及其边界条件。当n足够大时,试探函数就趋近于真解。,介绍两种常用的权函数。 1 最小二乘加权余量法 设有满足边界条件的试探函数 带入控制微分方程将产生余量 如果希望余量R在最小二乘的意义下为最小(即令R的平方和为最小),则构造,使 比较 可知,权函数
4、 通过求解 可求出ci,进而得到 。,2 伽辽金加权余量法 如果选用试探函数 中的试探函数项i作为权函数Wi,就成为伽辽金加权余量法。即 在许多物理问题控制微分方程的有限元法求解过程中,都采用伽辽金加权余量法推导有限元计算格式。,7-2 二维流体流动的有限元计算格式 二维稳态可不压缩流体流动方程由连续方程和动量方程描述 方程中的待求变量为流体速度u,v和压力p。,根据有限元法的计算思路,首先选取插值函数来近似描述速度u,v和压力p在单元内的变化情况。 式中N单元形状函数; ui, vi, pi单元节点处的速度和压力值。,这里的插值函数,就作为加权余量法中的试探函数;其中的形状函数,就作为加权余
5、量法中的权函数。,经过推导和简化,可得单元方程为 或 其中 分别为单元刚度矩阵、单元节点列向量和单元节点受到的来自“铰链”的“节点力”。,将 在求解区域内分别按节点号叠加,就可以构成整个流场的有限元计算的总体方程,7-3 流场有限元分析的几个特殊问题 1 速度和压力插值函数的阶次。 速度插值函数要高于压力一阶,否则方程会出现“病态”。 2 主对角线元素为零。 采用罚函数法,将压力用速度表示。求出速度后,再计算压力。 3 刚度矩阵不对称。 原来介绍的压缩存储方法全部没有用。 4 非线性方程组 需要迭代计算。,小结: (1)本章讨论了利用有限元法求解非结构问题的又一个例子流体流动问题的计算。所用方法为加权余量法,通过将试探函数带入控制微分方程,基于使所产生的误差(余量)在加权平均的意义上等于零的思想,来推导该控制方程的有限元计算格式。 (2)本章简要介绍了加权余量法的基本概念,最小二乘加权余量法中权函数的计算,以及伽辽金加权余量法中权函数的确定。对于无法利用变分原理,即找不到等价的泛函极值问题的控制微分方程有限元求解问题,一般来讲,都可以利用加权余量法推导其有限元计算格式。,(3)本章简述了二维稳定流场问题的有限元计算格式的推导思路,同时涉及了流场有限
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 互联网在线审批系统技术设计方案
- 2021年中国痛风与高尿酸血症临床诊治指南
- 2024年红细胞类诊断抗原项目建议书
- 物流营销与客户关系 课件全套 张广敬 项目1-9 初识物流营销和客户关系-评估物流营销绩效
- 电子线路装调 课件 A1-5能正确识别晶闸管
- 2024年各类型谱仪(含多道系统)项目建议书
- 2024年金融私募项目发展计划
- “双减”背景下区域教研如何赋能学科教师
- 2024年PP-R再生料项目合作计划书
- 2024年加气混凝土砌块项目合作计划书
- 核潜艇发展史
- 巨大儿难产的护理查房
- 2024年安徽江东控股集团有限责任公司招聘笔试参考题库含答案解析
- 兰州地产项目分析报告
- 平价饰品市场分析报告
- 交警大队业务培训方案
- 高光谱成像行业发展趋势报告
- 2024年1月安徽省普通高等学校招生考试适应性测试物理试题(解析版)
- 习作单元精读课教学策略探析以《搭船的鸟》为例
- 2024年湖南省建筑工程集团招聘笔试参考题库含答案解析
- 30题纪检监察位岗位常见面试问题含HR问题考察点及参考回答
评论
0/150
提交评论