当前课程知识点：有限元分析及应用 > 第3讲 针对复杂几何形状变形体的力学描述(1) > 3.4 平面问题的平衡方程构建 > Video 3.4
好，下面我们来具体构建平面问题的平衡方程
按照刚才说的建模的思路
对于平面的任意的复杂几何形状的变形体
我们从内部取出一个微小的面元
我们叫dxdy
也就是它的x方向的宽度是dx
y方向上的长度是dy
注意这一点的坐标，角点的坐标是x,y
那么这个厚度呢是t
平面问题我们叫等厚度
我们把这个微小面元叫dxdy_t
刚才也提到了，在这个微小面元上
每一个侧面，都有任意的力的作用
我们把它分解以后，两个方向进行分解
并且除上相应的面积
我们就定义相应的应力
那么我们现在把坐标标上
对于最左边的这个侧面，坐标是x,y
那么它的应力呢σxx
沿着y方向，向下的这个呢
第一个下标呢，是表示力的方向是y
然后这个面的法线方向呢是x
所以我们叫σyx
或者叫τyx
那么对应着右边这个侧面
垂直于这个面的力，我们叫σxx
第一个下标表示力的方向
第二个下标表明这个力所作用的面的法线方向是x
注意这一个应力它是有坐标的
它的坐标呢是在右边的这个面
也就是说它的x坐标是x+dx
当然y，还是y的变化
同样呢，沿垂直方向的分解
我们的下标，力呢是y
这个面呢法线方向是x
所以是叫τyx
那么它的坐标是x+dx
也就是说最右边的这个面的坐标
同样y方向的坐标是y
对于下表面
它的应力呢，同样
沿着y方向的这个应力
第一个下标是y
第二个是由于这个侧面它的法线方向是y
所以说它是σyy
那么它的坐标呢是x和y
那么同样呢，沿着x方向的力的分解
我们得到的这个应力，也就是剪应力
它就是τxy
它的坐标呢是x和y
对于上表面，一样
我们得到了σyy
那同样它的坐标要注意一下
x不变，y那就是y+dy的这个表面
另外沿着x方向的这个分力
我们得到的相应的应力呢是τxy
那它的坐标呢是x,y+dy
另外呢我们还要考虑
这个物体呢，它有体积力
那么体积力呢在这个微小的面元里面
我们同样用两个方向的分量来表达
有bx，上面一横，还有by，上面一横
刚才的这些呢，我们简单归纳一下
我们有一个约定
正面正向为正，关于力
负面的负向为正
所以说这样的话，大家在画力的方向的时候呢
我们统一按照这个约定来画
这样就比较简洁一些
另外呢后面我们在构建平衡方程的时候
要注意这么几个点
一个呢是四个侧面
在平衡方程中间一定是考虑合力的平衡
注意，合力就是
刚才的应力一定要乘上相应的面积才能得到合力
那么我们按照合力平衡来构建力的平衡方程
这也是非常重要的一点
另外还有一点要注意的
所有的应力，在不同的侧面，它的坐标不同
由于坐标有dx或者dy的变化
那么它有相应的增量变化
这一点一定要注意
增量怎么表达，我们后面给大家再具体的写出来
另外呢由于我们取出来是dxdy_t的微小的面元
所以说在每个侧面上，这个应力
不管是剪应力还是正应力
我们都认为它是均匀分布
因为dxdy，用极限的概念来说，我们总可以取得很小
那么在很小的范围内，它的变化是非常小
我们可以认为它是均匀分布
那么对于我们所得到的dxdy_t的微小的面元
我们需要建立它的合力的平衡
那么对于一个面元，他实际上也是一个体
也就是说它的厚度为t
那对于体来说，在平面上的一个体
那么我们需要建立它的平衡关系
有三个
一个是x方向所有的合力平衡
还有沿y方向所有合力的平衡
另外所有的合力
关于这个体里面任意一点都有一个力矩的平衡
这是我们三个平面问题的平衡方程
那么我们在构建这个三个平衡方程的时候
要注意一下，刚才提到了
在不同的侧面它由于有dxdy它的坐标的变化
所以说要算力学变量的增量
这个增量呢要用泰勒级数进行展开
我们要计算它的展开项
另外呢就是要算它的合力平衡
关于力学变量的增量的计算
我们要用泰勒级数来展开
我们举个例，比如
σxx这个正应力
它的坐标呢是x+dx，然后这边是y
在最右侧面
我们怎么来算它的这个增量变化引起的应力分量的变化
我们用高等数学的泰勒级数展开
大家要知道，这个x只有一个dx的增量变化
所以说呢我们这个泰勒级数
就针对x做一个针对dx的泰勒级数展开
那这样的话，我们第一项，就是σxx，σxy
这第二项呢就是由于dx的增量变化引起的增量
那这个增量变化由于是x方向的增量变化
所以说它是一个偏x的偏导乘上dx的微量的增量
那二阶的增量变化就是针对x二阶偏导
再乘上一个dx的二次方
同样还可以有高阶的
那么为了简化，刚才我们前面也提到了
就是我们小变形，还有一些这种基本的假定
实际上我们就是要忽略高次的这个级数的增量变化
也就是我们把二阶增量以上的量忽略掉
我们就取一阶的增量
这样的话，我们就可以得到
由于坐标的变化所引起的应力，或者是各个力学分量
它的增量的变化值
好了我们再具体看这个微元体
它呢我们把所有的各个分量都写出来以后
我们来建立x方向的合力平衡
y方向的合力平衡
还有针对任意一点的力矩平衡
那么以x方向的合力平衡为例
我们来看看所有x方向的合力，我们把它标出来
在最右边的这个侧面，我们有σxx
当然呢它的坐标呢是x+dx,y
我们把第一项应力乘上相应的面积，把它写在这儿
这是第一项的合力
沿着x方向的正方向，我们把这个合力为正
第二项是最左边的这个侧面
它也是σxx，它的坐标是x,y
那么它就直接乘上所作用的侧面的面积，叫dy乘上t
由于这是最左边，它是向负的y方向指向的
所以说我们这个力呢是一个负的
那我们再看看最上面的这个剪力
也就是τxy，那么它的坐标呢是x,y+dy
就是在最上面的这个面
那么我们把这个剪应力乘上相应的最上面的这个侧面的面积
叫dx乘t，这个是作为最上边的合力
最底边的这个面，沿着x方向
它的应力呢是τxy，它的坐标是x,y
所以说它乘上相应的最底边的侧面的面积，dx乘t
那么它是负的x方向
所以说这个地方我们取得是这么一个负的合力
另外呢我们还有体积力bx杠
bx杠呢这个乘上相应的体积
我们叫dxdy再乘上厚度t
这样的话就是体积力的合力
我们把这五项加起来，它应该平衡
那么再把刚才说到的基于坐标轴的增量变化
就是说x+dx，这边还有y+dy
用前面提到的泰勒级数展开
把它展开到一阶增量项
这样的话我们重写一下这个平衡方程
那么我们就可以看见有两项
前面σxx（x,y）这两项抵消了
后面这个τxy也抵消了
所以最后剩下三项
我们把这三项再同时除一下这个体的体积dxdyt
这样的话我们就得到这个偏微分方程
我们叫偏σxx，偏x
偏τxy，偏y，加上bx一杠等于零
大家可以看看，这是x方向的合力平衡
我们验证一下
我们前面说了
第一个下标表明是这个应力的力的方向
我们看看第一个下标都是x
表明我们这个方程就是关于x方向的力的平衡方程
同样，我们在y方向
用同样的方式进行所有合力的平衡构建
那么我也可以得到ΣFy等于零
也就是说是y方向的合力为零
也按照刚才的处理方式
由于有坐标轴的增加的面的增量
用泰勒级数展开，我们展开到一阶增量项
然后同样来进行一个处理以后
我们可以得到一个沿着y方向的这么一个平衡方程
我们也做一个验证
这个里面的应力分量的第一项，下标是y
那这表明我们这个方程它是y方向的合力平衡
当我们再看一看这两个偏导，我们也可以解释一下
这个比如第一项，偏yy，偏y
那么第一个下标呢是y方向
第二个呢是对y方向求一个增量
那么后面第二项呢是τyx
这是在yx方向，有一个增量的梯度变化
最后是在y方向有一个力的平衡
这表面前面两项是由于面的增量的变化引起的梯度变化
和我们的体积力要达到一个平衡
就是说这个偏导的含义是由于
不同的面坐标轴引起的梯度变化
所得到的这么一个增量关系
同样，关于合力，也有一个力矩的平衡
也就是说第三个平衡方程
我们叫ΣM0等于零
同样呢，你可以取任意一点
作为简单一点，我们可以取到面元的中心位置
来取它的力矩的平衡
同样呢我们也可以得到四项力矩、两项体积力
因为正好是过这个面元的几何的形心
所以说呢，它的力臂为零
力臂为零么，它的力矩就为零
所以我们只有四个侧面剪应力求矩的这么一个平衡
那么四个面的正应力
我们假设正应力都是均匀的
所以说它的合力也是通过面元的质心的
所以说它的力臂也是为零
所以说也没有力矩
所以说只有四个剪应力的这么一个力矩
那么我们四个剪应力得到的这么一个方程
同样也用泰勒级数展开
关于不同的面的增量用泰勒级数展开
再把展开过后的高阶量给它消掉
剩下的最后就得到τxy等于τyx
这个学过材料力学的同学应该都很清楚
它就是剪应力互等定理
实际上我们关于力矩的平衡就是得到剪应力互等定理
那么这个关系呢以后我们在平衡方程中间可以直接把它代入
就可以用了
综合以上三个平衡方程
也就是说x方向上的合力平衡、y方向上的合力平衡
力矩的平衡
我们就可以写出这三个平衡方程
刚才说到了我们有剪应力互等
我们在应力分量里面有四个
σxx、τxy、σyy、τyx
那么我们注意一下
这个剪应力互等，就可以把这两项代入以后就写一项
这样的话我们把这个方程用一下
把这三个方程变成两个方程
变成两个方程以后呢，剪应力我们就写一个
就是τxy或者τyx
因为剪应力互等，它是相等的
这样我们就比较简单的把平衡方程写成两个
那这个时候的独立的应力分量就是三个了
也就是σxx、σyy，还有一个就是τxy
另外一个τxy呢就是剪应力互等
实际上就是张量对称，它也是对称的这么一个分量
那么我们也可以把它写成指标形式
那怎么写呢？我们看一看
首先第一个下标，表明力的方向
那么第一个方程呢是x方向
第二个方程呢是关于y的方向
那么导数这个也是关于x坐标求一个导数
关于y坐标求一个偏导
那么我们把这个方程呢，就可以写成
σij偏导j加上bi一横等于零
那么我们注意一下
大家可能第一次碰见这么一个偏导的记号
我们这个地方注明一下
也就是说在下标里面带有逗号的这个我们叫作偏导
那么对于这个变量σij,j
就表明σij针对j的方向求偏导
那么写成完整来说应该写成σij偏导xj
那么我们用ij偏导j这个下标来表达
那么由于它要求和
由于jj是哑指标，它要求和
求和的话实际上就对应着前面这两项展开
在平面问题里面就是这两项
那么i呢分别变1、2也就是变x方向y方向
所以说是两个方程
那么这个呢就是我们用指标形式写出来的
比较简单的平衡关系
-有限的单元 无限的能力
--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
