当前课程知识点：有限元分析及应用 > 第7讲 杆、梁结构的有限元分析 > 7.6 门型框架结构的实例分析及MATLAB编程 > Video 7.6
下面对门型框架结构的有限元分析
给出完整的有限元的计算的流程
我们看一下
这是一个门型的框架结构
三个梁，我们分两个垂直梁，上面一个横梁
横梁上面作用一个均布载荷
我们把这个节点分为
1号节点、2号节点、3号节点、4号节点
在1号节点还作用有水平方向的力
3号节点和4号节点是固定的
结构的参数有
弹性模量、惯性矩、横截面面积
这三个梁都一样
首先对这个结构进行离散化和编号
我们对这个门型结构
首先把最上面这个分为①号单元
左边这个是②号单元，右边这个是③号单元
那么①号单元对应的节点编号是1、2
②号单元是3、1
③号单元是对应着4、2
首先我们对①号单元，由于它是分布力
我们把这个分布力作一个等效
按照我们上一个知识点所讲的
我们可以分别得到1号节点、2号节点
y方向的等效力和它的力矩
这样我们就可以来看看节点的位移列阵
它是一个由4个节点组成的
每一个节点3个自由度
所以它是由12个自由度组成的一个位移列阵
分别是u1,v1,θ1一直到u4,v4,θ4
对于节点外载荷列阵
我们把刚才的分布力已经等效到节点上
所以我们这样列出来就是
1号节点它有x方向的力、y方向的力还有一个力矩
2号节点同样也是这样的，我们列出来
那3号节点、4号节点
因为它是约束，所以没有外力
但我们还有一个支反力的列阵
我们全部把它列出来
当然1号节点、2号节点肯定没有支反力
那3号节点、4号节点分别有x方向和y方向的支反力
还有转角，也就是力矩
那这样的话我们就有6个支反力
我们把节点的外载和支反力合到一块
写一个总的节点载荷列阵
列出来就是这样的
下面对各个单元进行描述
单元①局部坐标和整体坐标是一致的
它的刚度矩阵是一个6X6的
因为它是一般平面问题的梁单元
我们查询一下前面梁单元的刚度矩阵的一些元素
那么我们分别算出来
②号单元和③号单元实际上情况相同
它们都是垂直的，只是节点编号不同
所以它的局部坐标系下的单元刚度矩阵也是一个6X6的
它们的各个元素都是相同的
坐标转换矩阵
由于它是垂直的一般的梁单元
所以我们把②号单元和③号单元的方向余弦分别来做一做
我们把这个6X6的转换阵也写出来
那么对②号单元和③号单元
我们分别前乘一个坐标转换阵的转置
后乘一个坐标转换阵
这样得到的是在整体坐标系下的②号单元和③号单元
各个矩阵里面的元素的系数
但是要注意，对于单元②
它对应的节点位移是u3,v3,θ3,u1,v1,θ1
对于单元③来说
它是对应着u4,v4,θ4,u2,v2,θ2
也就说它的节点是不相同的，元素是相同的
我们把这三个单元进行组装
就可以得到整体的刚度方程
那么我们得到的刚度方程是
其中总体刚度阵分别由三个单元的刚度阵进行组装
对于节点位移实际上也是由三个单元的节点位移
把它写成一个总的
虽然我们表达成一个Σ的关系
实际上写成一个总的
对于总体坐标系下的所有的节点的载荷
包括施加的外力和支反力
同样我们也把它写成一个总的
装好以后的总的刚度方程
我们在这个基础上要处理边界条件
我们看一下这个结构
节点3和节点4它是固定的
也就是说u3,v3,θ3以及u4,v4,θ4
是等于0的
那么对应着总体刚度阵的
这几个自由度的行和列我们把它划掉
剩下的就是我们处理边界条件以后的刚度方程
这个时候得到的刚度方程就是这样的
这个节点位移是
节点1是u1,v1,θ1
节点2是u2,v2,θ2
所有的这些载荷都是所施加的外载荷
所有的支反力都和刚才的位移边界条件约束对应的地方
都已经去掉了
我们先求这个节点位移
我们求出来后就得到了u1,v1,θ1,u2,v2,θ2
相应的这些值
按照我们前面的流程
这个时候我们就可以求支反力
还可以求单元的其它物理量
那么这个地方我们省略
基于MATLAB编程
我们对这个门型框架结构进行分析
同样，节点的离散化和编号也和前面一样
各个单元的描述
我们在MATLAB环境下
需要输入弹性模量、横截面面积、惯性矩、长度
由于单元②和单元③的刚度矩阵相同
所以我们针对单元①，还有单元②和单元③
分别调用两次Stiffness来计算单元的刚度阵
那么每个单元的刚度阵都是6X6的
分别计为K1和K2
这是MATLAB窗口的输入和输出情况
这是弹性模量、这是惯性矩、横截面面积、长度
然后调用两次，k1、k2，这样得到的刚度矩阵
那么我们还需要对单元②和单元③进行一个转换
把它转成整体坐标下的形式
这样我们来组装整体刚度方程
由于这个结构一共有4个节点
总的自由度是12个
所以我们要构建的总的刚度矩阵KK是12X12的
我们首先对KK进行清零
然后3次调用Assembly进行刚度矩阵的组装
那么我们看看，这个T是坐标转换矩阵
我们对k2前乘一个T矩阵的转置
后乘一个T矩阵，就得到k3
那么我们对总的刚度矩阵，12X12的先清零
然后对k1调用一次
因为k3是单元②和单元③相同的刚度阵
所以我们对k3分别调用两次
但是它的组装的位置有所不同
这样就得到组装以后的刚度阵12X12的情况
边界条件是3号节点和4号节点固定
也就是相应的x,y方向还有转角都为0
我们只要针对1号节点和2号节点的自由度进行求解
我们从总体刚度矩阵里的前6列和前6行
我们把它取出来
取出来后置给k
然后再生成对应的载荷列阵
采用高斯消元法进行求解
那么我们看一看
首先这个总体刚度阵的kk
我们把它前6行和前6列取出来
取到k里面
同样我们生成一个载荷的矩阵
然后对取出来的这个k矩阵
进行高斯消元法
就是k\p
这就是对k矩阵进行高斯消元
最后我们得到的u就得到了节点的位移
将前面得到的1号节点和2号节点的节点位移
以及3号节点和4号节点为0的情况
也就是说它们是固定的BC(u)的情况
我们把它写到一块，写成12X1的
再代回原来的整体刚度方程
这样算出来就是所有力的节点力
也就是P（12X1）的
这样我们计算的结果就是这样的
我们就可以按照总节点位移和支反力对应的关系
也就是说在哪个位置就可以找到我们支反力的情况
在下一个知识点
我们将基于ANSYS平台
对于门型框架结构进行一个分析
-有限的单元 无限的能力
--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
