当前课程知识点：有限元分析及应用 > 第5讲 变形体力学方程求解的试函数方法的原理 > 5.1 变形体力学方程求解的主要方法分类及试函数方法 > Video 5.1
同学们好
首先回顾一下上一讲的主要内容
上一讲我们学习了平面问题的几种特殊情况的讨论
一维拉杆问题及平面纯弯梁的描述和求解
空间弹性问题的完整描述
以及关于张量的描述和理解
这一讲我们将学习
变形体也就是弹性体力学方程求解的
方法分类及试函数方法
平面纯弯梁求解的试函数方法
也就是残值处理法
如何降低对试函数的高阶导数和边界条件的要求
平面纯弯梁求解的虚功原理
平面纯弯梁求解的最小势能原理及变分基础
最后我们将学习一般弹性问题的能量原理
变形体特别是弹性变形体的
力学方程求解方法主要有两类
一类是主要针对原始方程进行求解
我们在第一讲中已经作了介绍
主要的方法有解析法、半解析法
还有差分法，还有其它的一些需要技巧的方法
第二类是间接针对原始方程进行求解
主要是进行误差的处理
这也是我们本章节要进行介绍的内容
主要的内容有加权残值法
虚功原理、最小势能原理
还有一些变分的方法
主要的特点是要设置满足边界条件的一些解
主要是试函数
当然试函数中间含有待定的系数
第二个要点就是将试函数代入原来的控制方程
由于是试函数，所以它不是精确解
所以一定有误差
所以说要研究怎么处理误差
通过处理误差来求出待定系数
那么我们针对这个变形体
特别是弹性的变形体
我们详细地来叙述一下
试函数方法的应用背景和基本原理
我们都知道复杂形状的变形体
首先几何形状Ω是非常复杂的
当然，Ω的边界我们叫偏Ω
它包含两个部分
一个是Su，一个是Sp
它也是非常复杂的
那么我们要获取变形体中间的一些力学参量
我们定义了三大类变量
当然三大类变量要通过三大类方程来进行求解
前面说的直接方法也就是直接
通过三大类方程来求解三大类变量
当Ω，也就是几何形状非常复杂的时候
要得到解析解往往是非常困难的
那么试函数方法往往不是直接求解
它通过首先满足一定边界的试函数
先设定一定的试函数
然后再代入控制方程中间来处理误差
这样的话就使得求解的难度大大降低
总的流程基本上是这样的
第一步是设定试函数
也就是含有待定系数，满足边界条件
第二步就是代入控制方程中使得误差最小
可以获得待定系数的一系列方程
那么常用的方法是加权残值法
加权残值法实际上就是对误差进行处理的方法
加权残值法有Galerkin加权残值法
还有残值最小二乘法
那么关于加权残值法的思路可以参见
知识点1.3的微分方程的求解方法
好，我们来具体介绍Galerkin方法
首先我们对这个问题的一般的题法作一个介绍
我们还是以简支梁受均布载荷的
这么一个平面梁的弯曲问题为例
我们要求得它的挠度
首先它有一个控制方程
还有两类边界条件
那么我们抽象地表达这个控制方程是
L()这是一个算子
这个算子是关于挠度的算子
它实际上是一个微分算子
那么加上一个力b一杠等于0
我们还有两类边界条件
我们叫gu，就是关于位移的边界条件
另外还有关于力的边界条件
那么我们把它具体化
对于我们这个平面梁的弯曲问题
那么控制方程我们前面已经推导过了
它是一个EI乘上一个
关于挠度的4阶导数的这么一个方程
边界条件BC(u)就是两边
也就是两个端点的挠度为0
就是说x=0的时候v=0
x=l的时候v=0
另外还是在两个端点上
力的边界条件就是弯矩为0
我们知道弯矩可以表达成挠度的两阶导数
所以说是两阶导数在x=0的地方要等于0
同样，在右端点两阶导数在x=l的地方也要等于0
那么在这个问题里面
算子L()实际是就是EI关于x的4阶导数
载荷b一横就是-px
那么试函数的几个要点
首先我们要设定满足所有边界条件
在我们这个问题里就是BC(u)和BC(p)
的试函数我们叫v(x)
我们给试函数有一个标记就是上面有一个三角
那么把这个试函数代到边界条件里面去
它是完整地满足
那么我们设定试函数主要是
用基底函数的方式进行设定
我们基底函数可以取一系列的
满足BC(u)和BC(p)的基底函数
我们把它进行线性组合
那么它的组合的系数
我们叫c1,c2一直到cn
那么这些系数我们叫待定系数
这个是需要求解的
我们把这个试函数代到控制方程中
也就是说L()这个算子
关于试函数的这么一个算子加上一个b一横
代进去以后，由于它是试函数，它不是精确解
所以说它一定是不等于0的
我们把代入控制方程以后的这个函数叫残值函数
注意这是残值函数，它不是一个值
对于这个梁的弯曲问题
它是关于0到l整个沿着这个梁的长度
它是一个分布，是一个函数
它是关于定义域上的分布函数
由于它是误差的分布函数
所以它有可能为正，也有可能为负
因此需要再定义出一个误差的指标
这个指标应该说是可以处理整个定义域内的
所有的误差的总和
当然刚刚提到了有正有负，那么怎么处理
如果正负抵消了，那么这个误差指标就体现不出来
所以说关于误差指标的定义，这也是一个关键
那么加权残值法具体实施起来有这么几类
一个就是Galerkin加权残值法
它的要点就是把试函数前面取的那些基底函数
也就是φ1，φ2一直到φn
把它作为一个权函数
这个权函数对于刚才所说的残值函数进行加权
在整个域里面进行积分
让它积分后等于0
由此来确定我们试函数中间的一些待定系数
残值最小二乘法主要是要定义一个误差的指标
我们叫Err
它就是把残值函数平方
刚才说了有正有负，平方以后就完全是正的
然后再加权以后对整个域进行积分
所谓积分就是把所有的误差累加起来
让这个累加起来以后的全域上的误差
我们定义的误差指标让它取最小
当然这个里面也有一个权函数
这个权函数一般可以取为1
另外还有一个是配点法
由于这个残值函数是在定义域里面是一个分布函数
配点法就是让这个残值函数
在我们规定的一些关键离散点上让它为0
由此来确定我们试函数中间的一些待定系数
这就是强迫残值函数在规定点上它的余量为0
通过它来求取待定系数
我们重点介绍一下Galerkin加权残值法
要点就是试函数我们取为一系列基底函数的线性组合
基底函数都是满足所有的边界条件
那么线性组合的那些系数是我们的待定系数
也就是说c1,c2一直到cn
那么残值函数就是把试函数代入到控制方程里面去
获得的这么一个函数
是在定义域里的一个分布函数
Galerkin加权残值法就是前面提到的
把针对残值函数前面乘上的这个权函数
也就是加权，取为本身的基底函数
那么你有多少项，我就取多少个
然后在整个域里面加权以后使得误差为0
那这样的话我们就获得一系列的方程
这个方程是关于待定系数
c1,c2一直到cn的线性方程组
通过求解这个方程组
来求出c1,c2一直到cn
残值最小二乘法也是
首先和Galerkin加权残值法一样
设定的试函数也是φ1，φ2一直到φn
然后通过c1,c2一直到cn，进行线性组合
那么同样也有一个残值函数R
定义一个误差的指标叫Err
这个Err刚才说到了
它是残差函数的平方再乘一个权函数，对整个域作积分
一般情况下可以把这个权函数Wt可以取为1
那么取极小值就是通过变化里面的c1,c2一直到cn
使得新定义的这个，我们叫误差指标Err，让它最小
所以说它的一阶导数
关于c1,c2一直到cn的偏导数都为0
同样也可以获得关于c1,c2一直到cn的线性方程组
这样的话也就求出c1,c2一直到cn
那么就可以得到我们所设定试函数的
基底函数基础上的真实的解
-有限的单元 无限的能力
--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
