当前课程知识点：有限元分析及应用 > 第10讲 有限元方法中的基本性质 > 10.3 边界条件的处理与支反力的计算 > Video 10.3
我们得到了整体刚度矩阵为
当然我们前面已经知道了
所装配的刚度方程，这个K矩阵是奇异的
这就是因为它有刚体位移
所以我们要处理边界条件
我们常见的边界条件有两种
一种是零边界条件，就是比如固定的约束
另外还有一种给定具体数值的位移边界条件
我们把它表达成
那么前面我们也提到了
我们对刚度方程进行分块
也就是说把位移分成已知的qa一横
还有一个未知的qb
同样把力也分成已知的Pb一横，还有未知的Pa
前面已经提到了
已知的位移和未知的力，它形成这么一个互补的关系
我们要对这个分块的刚度方程进行边界条件的处理
首先我们讨论直接法
直接法是针对零边界条件，也就是说
那么在这个分块的刚度方程里面
我们就把qa一横对应的行和列划掉
剩下的这个我们就可以直接得到
关于未知节点位移的这么一个方程
这个方程我们就可以直接把它求出来
同样，针对给定数值的位移边界条件
也就是给定u一横的这种情况
我们对刚才分块的这个刚度方程
我们对于第二个方程，我们看看
我们把给定的qa一横代进去
我们由方程b，就是第二个方程
同样也可以直接把未知的节点位移qb算出来
当得到了qb以后，我们再由方程a就可以求出Pa
Pa也就是支反力
那么直接法，我们看看
它的处理过程是比较直观的
求解的矩阵的规模变小了，因为划行划列划掉了
适合于手工处理
但是由于划行划列以后
刚度矩阵的节点编号和排序改变了
不利于计算机的规范化处理
特别是大规模的计算
因为秩序改变了过后要重新排序
这个计算量和工作量是相当大的
另一个方法是置1法
置1法我们主要是针对零边界条件的情况
假定在r这个自由度的位置有一个零边界
那么置1法就是把r这一行这一列作一个处理
也就是把krr置为1，它的行和列其它元素置为0
然后再把所对应的载荷位置pr置为0
这就是置1法
置1法对不对呢
我们要考查一下对原始的边界条件是不是能反映出来
同时也要看置1法对其它节点位移的影响是不是合理的
我们看看，我们考查第r行的方程
我们把它拿出来
由于其它的元素都为0
只有krr还有qr，我们把它写到这
然后还有一个载荷pr
由于我们置1法,krr=1
这个对应的pr=0
所以只好得到
那这正好就是我们的零边界条件所需要求的
当然，我们得到节点位移qb以后
我们同样也可以由前面这个分块方程计算支反力
置1法的特点是
只能处理零边界条件
它可以保持原矩阵的规模，就说不需要重新排序
而且使得总刚度矩阵的对称性得到了保持
有利于计算机的规范化处理
乘大数法
乘大数法它是针对非零的边界条件和零边界条件都可以
我们看看，对于一个非零的边界条件u一横
它同样也是在第r个自由度的位置
我们对第r行，我们把它的系数krr乘上一个大数α
然后再把对应的载荷pr换一换
把它换成α乘krr再乘上一个u一横
这就是第r行的情况
那么我们考查一下第r行的一个关系
我们把这个乘了大数，又把对应的载荷作了一个置换
的这么一个情况，我们把这个方程拿出来
我们就得到这么一个方程
由于α这个系数很大
使得和它周边的这些元素相比
因为它比较大，所以其它的元素就可以省略
所以前面这一项就得到了
我们力的这一项，我们刚才已经置了
这样我们就把前面两项α乘krr消掉
剩下就是
这正好就是我们所针对的边界条件
当然我们这是一个近似的
为什么呢
因为前面在左边这个式子得到的时候
已经把其它的元素的影响去掉了
为什么去掉呢，就是因为α这个系数很大
所以α的取值将影响该方法的精度
一般情况下，我们可以把α取为刚度系数
就是kii里面最大的，再乘上一个10的4次方
取这么一个值往往就能得到一个比较好的精度
当然我们得到了节点位移以后
同样也可以进行支反力的计算
乘大数法的特点就是
它既可以处理零边界条件又可以处理指定边界条件
它保持原矩阵的规模，不需要重新排序
同时它也保持总刚度矩阵的对称性
便于计算机的规范化处理
也就是说只要改变这两个数
完全按照整个刚度方程的常规的求解来做就行了
下面再介绍一下罚函数法
罚函数法它针对非零的这么一个边界
它的思想就是用一个大刚度的弹簧
来支撑所对应自由度的这么一个边界
我们假定这个弹簧的刚度为C
这样我们就把这个支撑边界的位置算一下它的应变能
弹簧的应变能就是
也就是它的相对伸长量
整个系统的总势能
就是原来的势能再加上这个弹簧的势能
这个弹簧的势能刚才已经给了
由势能最小原理，我们对ui进行求导数
那么这个ui是分别对所有的自由度，从1到n
我们求了导数后就得到这么一个方程
这个方程就是刚度方程
当然这个刚度方程是带有弹簧的这么一个刚度方程
当然这个弹簧是我们人为加上去的
我们对这个刚度方程的第一行我们拿来分析一下
第一行这个地方就多了一个弹簧C
后面这个地方，同样由于求了偏导以后
后面的这个载荷这里也多了一个C
那么我们把C分别除一下就得到这么一个表达
由于C比较大，所以其它的那些项都为0
这一项也为0
这一项也是很小，它和1比很小
就得到了u1等于u1一横
当C足够大的时候那么这个可以说
是我们指定边界下的一种近似
当然我们这个C一般是取
所有刚度系数里最大的再乘上10的4次方
当然在近似的得到u1的位移以后
我们还可以直接求支反力
我们还由第一个方程来做，那么我们可以得到R1
我们把其它的消掉
大家可以看看
正因为u1和u1一横这中间有一个差
因为我们得到的是一个近似的解
正因为有这个近似的解，它有一个差
这个差应该是比较小的
再乘上这个刚度系数，这个C是比较大的
乘上以后恰好就可以得到这个支反力
罚函数法的最大好处就是可以直接计算出
位移边界上的支反力
-有限的单元 无限的能力
--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
