当前课程知识点：有限元分析及应用 > 第10讲 有限元方法中的基本性质 > 10.4 位移函数构造与收敛性要求 > Video 10.4
对于单元来说
位移模式将直接影响它的收敛性
什么叫收敛
也就是说当节点数或者单元插值位移的项数
趋近于无穷大时
也就是说当单元尺寸趋近于零时
最后得到的结果是无限地逼近准确解
也就是计算量越大、划的单元越多
我们的结果越好
我们看看这个收敛曲线
横坐标是节点数，纵坐标是势能
这是真实的解
那么1号曲线，随着我们节点数增加
它逐渐收敛到真实的解
这个是收敛的
2号曲线，我们看看
2号曲线它是收敛于真实解
但是收敛速度比较慢
3号曲线，它虽然收敛，但是没有收敛到真实解
4号曲线，它是一个非单调的收敛
它不构成真实解的上界和下界
5号曲线，它完全是一个发散的
那怎么才能保证我们所构建的有限元单元是收敛的呢
主要看单元的内部，还有单元和单元之间的连接
只有保证单元内部和单元的连接是收敛的
那么我们整个组装出来的这个模型
才能保证它是收敛的
由于我们在有限元里面的基本原理是要算势能
我们看看这个势能是要计算每个单元的应力、应变
具体表达出来就是
应变的转置乘上弹性系数矩阵，再乘上一个应变
这边是体积力乘上一个位移
面积力乘上一个位移
我们把它作一个展开
也就是说考虑当单元尺寸趋向于0的时候的展开
这个函数我们把它展开
第一项就是展开过后的常数项
第二项就是一次项，当然还有高次项
同样，对于位移这边也是一样
把它展开，第一项是常数项
第二项是一次项，还有高次项
同样，这个后面也是这样的
当单元尺寸趋向于0的时候
这个常数项应该要存在
这个一次项以及后面的高次项，它要趋向于0
这种情况下我们就可以说
单元有这种性能，它就是收敛的
为了得到这种性能
也就是说位移函数本身在单元上应该是连续
而且至少要使得单元位移函数和对应的应变
都为常数的项，也就是说要包含常位移项和常应变项
上面我们可以看到，也就是说这个常数项一定要存在
那么同样还有单元之间位移函数应该保证足够的连续
至少要使得应变的计算能够存在
这样我们就针对单元的内部有准则1
前面提到了，要使得我们的应变能或者势能存在
至少要存在常位移项和常应变项
我们准则1就是针对单元内部的，我们也叫完备性要求
具体描述是这样的
如果在势能泛函中所出现位移函数的
最高阶导数是m阶
那么有限元解答收敛的条件之一就是
单元内部的位移场函数的试函数
也就是我们的位移模式至少是m阶完全多项式
这个叙述起来比较繁琐
实际上它的意思就是要具有常位移项和常应变项
针对单元之间，前面我们说了就是要保证足够的连续
要使得应变的计算能够存在
那么怎么保证呢
我们就有一个协调性的要求，就是准则2
同样如果势能泛函里的位移函数出现的最高阶导数是m阶
那么要求在单元交界面上的必须有直至m-1阶的连续导数
也就是说具有Cm-1的连续性
那关于准则1、准则2的具体应用我们给出一个后面的讨论
对于平面问题，势能泛函是这样的
应力乘应变
我们看看这个应变
由几何方程我们可以得到这么一个偏导的关系
这个里面的最高阶导数明显能看出来就是1
所以说m=1
那么由准则1，形状函数至少要包含完整的直至m次
在这个里面m是等于1的
所以是一次多项式
也就是说对于二维问题
一定要包含这么一个完全的一次项
由它我们求偏导以后，可以得到
它是包含常位移项，还有常应变项
由准则2，位移函数要求
刚才是m=1，所以是要求Cm-1阶连续，也就是C0连续
什么叫C0呢
就是要求函数本身连续
也就是说仅采用节点位移进行插值
就可以使得应变的计算能够存在
那么我们再考查一下平面梁的弯曲问题
同样它的势能泛函写出来
平面梁呢，它是关于挠度的二阶导数
这是它的几何方程，是这样表达的
由准则1，也就是说这个时候它的最高阶导数是2
也就是说m=2
准则1就是要求m阶多项式要存在
这样我们就得到，它的挠度函数是
这是它可以保证常位移项和常应变项
你可以把它求两阶导数以后，你就可以发现
它的常应变是存在的
由准则2，它的位移函数的连续性是Cm-1连续
也就是C1连续
也就是说要求在单元之间的位移函数
至少要求一阶导数连续
这是什么意思呢
就是在节点位移上，我们既要有节点位移，也就是挠度
同时还有一阶导数，也作为它的节点位移
也就是θ，也就是转角
只有通过这两种
就是说节点位移和节点转角来进行插值以后
才可以使得应变的计算成立
也就是说前面提到的v(x)的常应变项
一定要取到x的平方
那必须要有4个待定系数
左节点要有2个位移条件，右节点有2个位移条件
由这4个位移条件才可以确定4个待定系数
这样才可以保证前面所说的常应力和常应变的存在
这样就可以得到插值函数为
对它进行两阶导数的求导以后
我们就可以说它包含有常应变项
同时还有应变的一次项
在常用的单元中间，能够保证常位移项
和常应变项的多项式，我们把它列出来
对于一维问题的杆单元
它是一个常数，还有一个x
这是杆单元
平面单元，要保证常位移和常应变的项是1,x,y
空间单元，也一样，是1,x,y,z
平面梁单元
要保证常应力和常位移的项就是1,x,x的平方
对于平面弯曲板单元，它要求
1,x,y,x的平方,xy,y的平方
对于单元的位移模式的选择
我们前面在构造单元的位移模式的时候已经提到了
第一个是唯一确定性原则
也就是说我们待定系数的个数应该与
节点位移的自由度数要相等
而且要选择多项式是从低阶到高阶
尽量选择完全多项式，这样可以提高单元精度
同样，我们在选择位移模式的时候
还要考虑单元内部的完备性，也就是收敛准则1
它必须要有位移的常数项，还有完备的一次项
对于单元之间的协调性，我们有收敛性准则2
对于C0型单元，它完全是可以保证的
对于C1型单元，往往很难保证我们所取的位移模式
比较难以保证单元之间的协调性
对于2D问题和3D问题我们来构造这个多项式
从低阶到高阶
我们可以参考Pascal三角形
再参考上面的一些准则来选择相应的函数的项次
-有限的单元 无限的能力
--Video
-课程大纲
--课程大纲
-第一章第一节测试题
-1.1 力学的分类：质点、刚体、变形体的力学
-第一章第二节测试题
-1.2 变形体力学的要点
-第一章第三节测试题
-1.3 微分方程求解的方法
--勘误
-第一章第四节测试题
-1.4 关于函数逼近的方式
-第一章第五节测试题
-1.5 针对复杂几何域上的函数表征及逼近
-第一章第六节测试题
-1.6 有限元的核心：针对复杂几何域的分片函数逼近
-第一章第七节测试题
-1.7 有限元发展的历史和软件
-课后讨论
--讨论题
-课后练习
-第二章第一节测试题
-2.1 弹簧的力学分析原理
-第二章第二节测试题
-2.2 弹簧单元与杆单元的比较
-第二章第三节测试题
-2.3 杆单元的坐标变换
-第二章第四节测试题
-2.4 一个四杆结构的实例分析
-2.5 四杆结构的ANSYS实例分析
--ANSYS
-课后讨论
--讨论题
-第三章第一节测试题
-3.1 力学描述的基本思路及关于变形体材料的基本假设
-第三章第二节测试题
-3.2 指标记法
-第三章第三节测试题
-3.3 关于三大变量及三大方程的思路
-第三章第四节测试题
-3.4 平面问题的平衡方程构建
-第三章第五节测试题
-3.5 平面问题的几何方程构建
--勘误
-第三章第六节测试题
-3.6 平面问题的物理方程构建
-第三章第七节测试题
-3.7 两类边界条件
-课后讨论
--讨论题
-第四章第一节测试题
-4.1 几种特殊情况的讨论
--勘误
-第四章第二节测试题
-4.2 简单拉杆问题的完整弹性力学求解
-第四章第三节测试题
-4.3 平面纯弯梁的描述及求解
-第四章第四节测试题
-4.4 空间弹性问题的完整描述
--勘误
-第四章第五节测试题
-4.5 关于张量的描述及理解
-课后讨论
--讨论题
-第五章第一节测试题
-5.1 变形体力学方程求解的主要方法分类及试函数方法
-第五章第二节测试题
-5.2 平面弯曲梁求解的试函数方法-残值处理法
-第五章第三节测试题
-5.3 如何降低对试函数的高阶导数的要求
-第五章第四节测试题
-5.4 平面弯曲梁求解的虚功原理
-第五章第五节测试题
-5.5 平面弯曲梁求解的最小势能原理的变分基础
-第五章第六节测试题
-5.6 一般弹性问题的能量原理
-课后讨论
--讨论题
-第六章第一节测试题
-6.1 基于试函数的经典方法与有限元方法
-第六章第二节测试题
-6.2 有限元方法中的自然离散与逼近离散
-第六章第三节测试题
-6.3 有限元方法中的基本步骤
-第六章第四节测试题
-6.4 经典方法及有限元方法的比较
-课后讨论
--讨论题
-第七章第一节测试题
-7.1 局部坐标系中的杆单元构建及MATLAB编程
-第七章第二节测试题
-7.2 局部坐标系中的平面纯弯梁单元构建及MATLAB编程
-第七章第三节测试题
-7.3 局部坐标系中的一般梁单元构建(组装)
-第七章第四节测试题
-7.4 梁单元的坐标变换
-第七章第五节测试题
-7.5 分布力的处理
-7.6 门型框架结构的实例分析及MATLAB编程
--【知识点7.6】 门型框架结构的实例分析及MATLAB编程(PDF)
-7.7 门型框架结构的ANSYS实例分析
--【知识点7.7】 ANSYS算例-门型框架结构有限元GUI操作与命令流(PDF)
-课后讨论
--讨论题
-第八章第一节测试题
-8.1 平面3节点三角形单元及MATLAB编程
-第八章第二节测试题
-8.2 平面4节点矩形单元及MATLAB编程
-第八章第三节测试题
-8.3 轴对称单元
-第八章第四节测试题
-8.4 分布力的处理
-8.5 平面矩形薄板分析的MATLAB编程
--【知识点8.5】 平面矩形薄板分析的MATLAB编程(PDF)
-8.6 平面矩形薄板的ANSYS实例分析
--【知识点8.6】 ANSYS算例-平面矩形薄板有限元GUI操作与命令流(PDF)
-课后讨论
--讨论题
-第九章第一节测试题
-9.1 空间4节点四面体单元及MATLAB编程
-第九章第二节测试题
-9.2 空间8节点正六面体单元及MATLAB编程
-第九章第三节测试题
-9.3 参数单元的原理
-第九章第四节测试题
-9.4 数值积分
-9.5 典型空间问题的MATLAB编程
--【知识点9.5】 典型空间问题的MATLAB编程(PDF)
-9.6 典型空间问题的ANSYS分析实例
--【知识点9.6】 ANSYS算例-典型空间问题有限元GUI操作与命令流(PDF)
-课后讨论
--讨论题
-第十章第一节测试题
-10.1 节点编号与存储带宽
-第十章第二节测试题
-10.2 形状函数矩阵与刚度矩阵的性质
-第十章第三节测试题
-10.3 边界条件的处理与支反力的计算
-第十章第四节测试题
-10.4 位移函数构造与收敛性要求
-第十章第五节测试题
-10.5 C0单元与C1单元
-第十章第六节测试题
-10.6 单元的拼片试验
-第十章第七节测试题
-10.7 有限元分析数值解的精度与性质
-第十章第八节测试题
-10.8 单元应力计算结果的误差与平均处理
-第十章第九节测试题
-10.9 控制误差和提高精度的h方法和p方法
-课后讨论
--讨论题
-第十一章第一节测试题
-11.1 1D高阶单元
-第十一章第二节测试题
-11.2 2D高阶单元
-第十一章第三节测试题
-11.3 3D高阶单元
-第十一章第四节测试题
-11.4 基于薄板理论的弯曲板单元
-第十一章第五节测试题
-11.5 子结构与超级单元
-课后讨论
--讨论题
-第十二章第一节测试题
-12.1 结构振动的有限元分析：基本原理
-第十二章第二节测试题
-12.2 结构振动的有限元分析实例
-第十二章第三节测试题
-12.3 弹塑性问题的有限元分析：基本原理
-第十二章第四节测试题
-12.4 弹塑性问题的有限元分析：非线性方程求解
-课后讨论
--讨论题
-第十三章第一节测试题
-13.1 传热问题的有限元分析：基本原理
-第十三章第二节测试题
-13.2 传热问题的有限元分析实例
-第十三章第三节测试题
-13.3 热应力问题的有限元分析：基本原理
-第十三章第四节测试题
-13.4 热应力问题的有限元分析实例
-课后讨论
--讨论题
-【基本建模Project1】2D问题：带孔平板的有限元分析
--Doc I-1
-【基本建模Project2】3D问题：花型卡盘网格划分的控制
--Doc I-2
-【应用建模Project3】振动模态分析：斜拉桥的模态分析
--Doc I-3
-【应用建模Project4】弹塑性分析：厚壁圆筒受内压的弹塑性分析
--Doc I-4
-【应用建模Project5】传热分析：钢制圆柱冷却过程温度场的瞬态问题
--Doc I-5
-【应用建模Project6】热应力分析：桁架结构的温度及装配应力分析
--Doc I-6
-【高级建模Project7】结构的概率：大型液压机机架的概率设计分析
--Doc I-7
-【高级建模Project8】p方法的建模与应用：平面问题的p型单元建模与分析
--Doc I-8
