当前课程知识点：有限元分析及应用 > 第4讲 针对复杂几何形状变形体的力学描述(2) > 4.1 几种特殊情况的讨论 > Video 4.1
同学们好
首先简单回顾一下上一讲的基本内容
上一讲我们介绍了
力学描述的思路及变形体材料的基本假设
介绍了三大类变量及三大类方程
同时也介绍了复杂变形体描述的建模思路
最后介绍了位移和力的边界条件
这一讲我们重点要介绍
平面问题的几种特殊情况的讨论
针对1D也就是一维的拉杆问题
及平面纯弯曲问题的梁进行描述和求解
还要介绍空间弹性问题的完整描述
最后我们对张量的描述和理解进行一个讨论
那么简单回顾一下
我们在二维建模里面
我们给定的几何Ω
在Ω里面我们给定体积力bi一横
同样我们在Ω的表面
也就是偏Ω呢，我们有两类边界条件
一个是Su，还有一个是Sp
在Su上面我们给定了位移的这个值
我们叫位移边界条件Ui一横
同样在Sp上面，在力的边界条件上呢
我们给定了边界的外力值Pi一横
这就是我们给定的情况
那么我们基于内部的建模和表面的描述呢
要定义三大类变量
第一呢就是位移变量
我们在二维问题里面有x方向的位移
和y方向的位移u和v
还有呢我们独立的应力的变量
我们有三个分量σxx，σyy，τxy
应变呢我们也一样
我们有三个独立的应变分量εxx，εyy，γxy
那么所有的三大类变量呢都是坐标点的函数
也就是说xy的函数，在二维问题里面
那么我们构建的相应的平衡方程呢
独立的平衡方程呢有两个
一个是x方向的合力平衡，还有y方向的合力平衡
几何方程呢，我们得到了三个
一个是x方向的相对伸长量的表达
还有y方向的相对伸长量的表达
以及夹角的变化
那么我们也构建了材料的物理方程
我们可以写出来应力和应变之间的相对的关系
这个关系呢是一个线性关系
下面，我们还有呢是两类边界条件
关于位移的边界条件
就是在Su上面给出u一横，还有v一横
那么在Sp上面呢，我们给定了力的分布
Px一横，还有Py一横
同样我们有两个的边界上的力的平衡方程
下面要就几种情况来进行讨论
这几种情况是一些特殊情况
比如平面应力、平面应变，还有刚体位移
那么平面应力问题呢，我们首先看一下对象
这个对象首先是个二维问题
是xy坐标系里面的一个平面问题
它是一个很薄的等厚度的薄板
那么所有的变量呢在z方向都不变化
z方向板的厚度是t
那么我们可以看出来
在薄板的内表面和外表面
也就是说在z等于正负二分之t的两个内外表面上
所有应力都为零，它是自由表面
所以说σzz，还有τxz，τyz
在两个内外表面上都是为零的
那么由于板比较薄，所以t呢很小
所以我们就近似的认为在板的中间
刚才的这三个分量也都为零
也就是σzz，τxz，τyz处处为零
当然这是一个近似的假定
那么由此呢我们看看
所对应的剪应力的这一块我们可以得到剪应变为零
也就是γxz，还有γyz等于零
那么对应着σzz的对应它的应变
比如εzz，是不是等于零
或者它不等于零，它等于多少
那么这个需要由物理方程来进行确定
也许大家会想，σzz等于零
那应力为零了，那这个方向的应变也应该为零呀
但是要知道，我们广义虎克定理里面两个特点
一个是主方向拉长，另外一个方向上会有一个泊松效应
也就是说，它会收缩
所以说，由于泊松效应会引起这个方向
也就是z方向的应变有可能不为零
所以说我们要用物理方程来进行推导
那么简单归纳一下我们平面问题的基本变量
也就是说我们的位移是
u(x,y)，v(x,y)两个分量
应力呢，σxx(x,y)，σyy(x,y)，τxy(x,y)
应变呢是εxx(x,y)，εyy(x,y)，γxy(x,y)
由于近似认为z方向处处为零
那么我们由z方向的物理方程
当然这个物理方程式三维的物理方程
这个三维的物理方程就是[方程]
可以明显看出来，这是z方向施加的力引起的应变
同样还有x方向和y方向作用力引起的泊松效应
在这个地方
那么刚才σzz等于零
我们看看由这个方程就可以得到[方程]
那么我们可以看出来，确实是在平面应力问题里面
z方向的应变是不为零的
它是等于由σxx、σyy泊松效应引起的
另外反方向的一个应变
那么第二个特殊情况讨论呢就是平面应变问题
那么我们同样也看一看平面应变问题的对象
它的对象是一个无限长得等截面的柱形体
比如这个水坝，当然我们是要等截面
并且所有的外载荷都不随着z变化
当然理想化的这种无限长的是没有的
我们一般来说，相对于横截面的尺寸
它比较长，我们就可以近似认为它是一个平面应变问题
那么从中间呢我们可以任意取出一个截面
这个截面也是一个薄片
那么它呢也可以看成是一个平面问题
但是在z方向，它和刚才的平面应力问题完全不一样
也就是说，由于它是一个无限长的等截面的柱形体
并且所有的外载荷不随z变化
那么任何一个在z任意位置上截下来的截面
它都是对称面
那么对称面就意味着它z方向的位移和应变都为零
也就是说w，就是z方向上的位移它是为零的
那同样呢，和z相关的应变它也是为零
γyz，γxz，它也是为零
那么由z方向的相对伸长量的计算
也就是说三维问题的我们计算它的应变
由于它的位移为零，就是说它完全被限制死了
它z方向不能任何的移动
那么同样，z方向的相对伸长量也就是应变它也为零
那么从应力上我们来看看找对应关系
由于和z相关的这些剪切应变为零
那么得到的和z方向相关的相应的剪切的应力它也是为零
那么在z方向的应力是不是为零
那我们要用物理方程来进行讨论
好，我们把平面应变问题的三大类变量稍微归纳一下
首先所有的力学变量都是xy的函数
不随着厚度z方向的变化而变化
那么我们可以得到位移同样也是x方向的位移u(x,y)
还有y方向的位移分量v(x,y)
同样呢也有σxx(x,y)，σyy(x,y)，τxy(x,y)
应变呢也是εxx(x,y)，εyy(x,y)，γxy(x,y)
由于刚才我们提到的平面应变的特征
就是z方向的应变相对伸长量为零
那么我们由z方向的物理方程我们来看εzz等于零
那么我们用物理方程就可以得到
σzz等于μ乘上σxx加上σyy
那么明显能看出来，在平面应变问题
z方向的主应力，我们也叫正应力，它是不为零的
σzz是不为零的
那么我们进行两个平面问题的比较
从前面的基本变量和方程我们可以看
我们把三维问题就平面问题进行简化
就是说z方向的变化我们不考虑
我们都可以得到相应的平面应力和平面应变的三大类方程
应该说平衡方程和几何方程式完全一样的
那么在物理方程有一些不一样
那么我们看平面应力问题的特征
就是z方向的特征是平面应力它的z方向的应力为零
有σzz为零，但是应变不为零
它有泊松效应，εzz等于负的μ除上E
然后乘上σxx加上σyy（此处有误，字幕为准）
那么平面应变呢，z方向的特征是应变为零，应力不为零
正好这两个是互补的，一个应力为零一个应变为零
那么我们分别由三维问题由这两个特征
分别得到相应的物理方程
我们平面应力得到的方程是这样的
那么平面应变得到的是这样的
我们比较一下这两个物理方程
我们就发现，将平面应力问题里面的弹性模量E
换为E除上1减μ平方
把平面应力问题里面的μ泊松比换成μ除上1减μ
这样我们就从平面应力问题的物理方程转化为
平面应变问题的物理方程
反过来，也可以由平面应变问题的物理方程转化成
平面应力问题的物理方程
我们第三个特殊问题的讨论就是平面刚体位移
什么叫刚体位移
也就是说这个物体只有刚体的移动或转动
它不会产生任何内部的应变或应力
那么这样的位移是什么样子，它怎么表达
那首先，我们基于前面的三大类变量
描述变形
我们有应变分量
我们都知道平面问题三个应变分量
它的物理含义，就是x方向的相对伸长量εxx
y方向的相对伸长量εyy
还有夹角的变化γxy
那么我们让这三个变形分量
就是描述变形的分量为零
也就是它不产生变形
这就是刚体运动的情况
那我们看看这个时候确定出来的刚体位移
也就是u和v，是什么样的形式
那我们得到的这三个应变为零的方程
我们分别来求一下
我们对第一个，它是偏u偏x等于零
那这样的话，我们就可以说
u函数，因为它本身是xy的函数
那么它要满足偏x等于零，它一定是一个y的函数
我们把这个函数取为f1（y）
那同样，对于v函数，它对y偏导为零
也就是说
这个v(x,y)的这个函数呢它一定是一个关于x的函数
我们把这个函数设定为f2(x)
把u函数和v函数对应的f1和f2
代到γxy这里面的u和v里面去
由于这个里面的偏导因为我们都变成
这两个函数都成为一个自变量的函数
所以我们就可以得到一个不是偏导的微分方程
也就是[方程]
也就是说剪应变要为零
这两个呢我们把它移一下项
我们把第一项移到右边去
最后我们就可以得到这么一个方程
我们分析一下，这个方程，左边是一个导数，一阶导数
但是我们看看它的自变量就是x
右边呢也是一个导数，但是它的自变量是y
那么是一个x的函数要等于一个y的函数
而且在我们所分析的对象Ω这个物体里呢要处处满足关系
由于它是刚体运动么，处处都要满足这个关系
它是一个恒等式
那一个x的函数要等于一个y的函数，要处处满足
那只能说让它们这个函数分别等于一个常数
这样才能做一个恒等式
那我们令它等于一个常数ω0
这样我们就可以分别就
f2和f1这两个一阶的常微分方程
来进行求解
这个求解结果非常简单[方程]
那么这个f1和f2实际上
和我们前面设定的这个一看一对照
我们就知道，它就是关于u和v的刚体位移的函数表达
也就是说u(x,y)它的函数表达是这样的
v(x,y)，就是y方向的位移分量的函数表达就是ω0x加v0
这就是我们平面刚体位移的函数表达式
注意这里面的三个常数
一个是u0 v0加ω0
它是表征平面刚体位移的常数
好，我们对这三个表征平面刚体位移的常数
分别要分析一下它的物理含义
那么我们分别看一看这三个表示平面刚体位移的常数
那么如果我们令ω0等于0，v0等于0
这个时候的u(x,y)就等于u0
那表明什么意思呢
表明处处x方向的位移等于一个常数
所以说u0的物理含义就是为
整个物体在x方向的刚体平移量
同样，v0也是，为整个物体在y方向的刚体平移量
那我们看看ω0，它的物理含义是什么
那么我们分别来分析一下
那么对于物体里面的任意一点这个坐标
这个坐标是xoy坐标系
对于它任意一点，xy这个坐标在这儿
这个坐标呢有刚体位移
刚体位移的描述呢，我们先不考虑u0和v0
把这两个去掉
那我们看看x方向的位移，就是负的ω0乘上y
也就是说，它是负方向的负ω0y，这是u
那么这一点的y方向的移动量是ω0x
那么这是两个分量
那么真正这一点的移动是这两个分量的合成
那么合成以后就等于这么一个我们叫d
那么d的这两个平方和相加再开根号
得到的呢就是r乘上ω0
那么这个r正好就是从原点到这一点的坐标的距离
这就可以看出来，那么它的位移量
就正好是距离r的这一点乘上一个ω0
小变形、转动的情况下乘上的一个转动的刚体位移
所以呢，ω0就是整个物体刚体转动的角度
已知u和v同样根据这个关系呢
可以求出来它的刚体的转动量
那么我们由前面这两个方程
我们u0和v0都不考虑的情况下、等于零的情况下
我们就可以u(x,y)对于y求偏导
它就等于负的ω0
对于v函数呢，对于x求偏导
它就等于ω0
把这两个一减再除个2刚好就是刚体位移的转动量
也就是说我们已知u和v
可以通过这么一个求偏导相减的关系呢
可以把刚体位移的转动把它求出来
-有限的单元 无限的能力
--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
