当前课程知识点：有限元分析及应用 > 第3讲 针对复杂几何形状变形体的力学描述(1) > 3.5 平面问题的几何方程构建 > Video 3.5
下面我们构建平面问题的几何方程
那么这是变形前的取出来的一个微元
叫dxdy，那么这一点的坐标叫x,y
那么它受力以后它会变形
变形以后，变成像这样一个菱形一样
我们分别定义它的几个位置
A点、P点、B点
变形以后，A’、P’、B’
也就是说对应点它相对移动了一个位置
那么移动的前提条件就是有一个位移场
u（x,y)、v（x,y）
那么这个位移场是表明这个变形体任何一个位置
它移动的位移，x方向上的分量和y方向上的分量
那就是u和v
那么如果变化x、y这个坐标就可以得到这个变形体
任意一个位置它沿两个方向的位移分量
好，我们要构建几何方程，就是要考虑两个问题
一个呢就是考虑变形以后的长度变化
也就是相对伸长量
另外还有一个，它形状改变以后有一个夹角的变化
那么我们定义沿着x方向的相对伸长量
定义为εxx，也叫x方向的正应变
那么我们来看
这个定义叫相对伸长量
也就是变形以后的绝对长度P'A'，这个长度
减去原来这个长度PA，就是绝对伸长量
再除上一个原长PA
这样的话就是我们的相对伸长量
那么我们具体已知它的位移场
怎么来做它的相对伸长量
我们来推导一下
我看一看，P一点是xy这个坐标
由于变形的描述
也就是说它的位移描述是u和v
那我们把这一点的坐标代进去
我们就得到了P'
也就是移动位移在x方向上的分量就是u
y方向上的分量就是v
这一点的位移，我们就得到了
A这一点的变形前变形后的移动量
同样把A这一点的坐标应该是x坐标是x+dx
y坐标就是y
我们把这个坐标代到u和v的这个函数里面去
就可以得到A一点的x方向的位移移动量和y方向的移动量
因为我们这个位移场么，不管你把什么坐标代进去
它都要满足这个位移场，用函数来表达
那同样，在x方向由于它有增量
所以说呢我们用泰勒级数展开
按照前面的展开方式，它就有一个
沿x方向的位移加上x方向的一阶梯度变化乘上增量
这就是A点的x方向的由于dx这一点的
坐标变化引起的x方向的位移的增量变化
那同样，A这一点变到A'这一点在y方向的位移量
把这一点的坐标，代到这个位移场函数里面去
那么在y方向的这个分量呢，用泰勒级数展开
写出它的一阶增量
就得到了A这一点它变形以后的y方向的分量
B这一点相对于P这一点呢
它是在y坐标有一个dy的变化
那同样呢我们把这一点的坐标
也就是说x,y+dy代入到位移场函数里面去u和v
那么也同样用泰勒级数展开
我们保留一阶增量项
同样也可以得到
B这一点变到B'这一点，x方向的位移
x方向的位移量就是u加上偏u偏ydy
y方向上的分量呢就是v加偏v dy，偏y dy
好，有了这三个点具体的位移的描述
那我们下面我们来计算P’A'
当然我们这个角度比较小
虽然是斜的，但我们就认为是一个沿x方向的水平段
那么P'A'我们怎么算呢
就用这个大的长段，减去前面这个小段
就可以得到P'A'
这个大的长段怎么做呢
就是用前面这一段PAdx，再加上这一段，就是大的长段
这一段大的长段我们做出来就是dx加上u
加上个偏u偏x dx
然后减去PP'
PP'就是u
那么我们做一个化简就是u和u消掉了
我们就得到了dx加上偏u偏x dx
这就是P'A'
然后还要减去一个PA
PA呢就是dx，就是原来的原长
所以说把dx去掉
所以说我们这一段上面的绝对伸长量就是偏u偏x dx
原长呢就是dx
那么把dx消掉以后就是偏u偏x
那同样，在y方向，我们定义相对伸长量
也是按照刚才的这种方式
就是P'B'减去PB，再除上一个PB
同样我们也可以最后得到一个偏v偏y
好了我们两个方向的相对伸长量，我们就得到了
下面我们来定义一下夹角的变化
我们看一下，原来是一个直角
现在变形以后，变成一个相当于菱形一样的东西
那么我们和x轴的这个夹角呢叫α角
和y轴的这个夹角呢我们叫β角
那同样我们来算一算α角和β角
那么我们看看α角
α角怎么算呢
我们来算这一段长度减去原来的这段长度
也就是说这一条边，再比上一个原长
这个dx变化以后的这个长度
和变形以前的这个长度应该是接近的
也就是说，可能它加上一个高阶小量
那么它呢可以去掉
那么上面v和v消掉以后
剩下的dx/dx去掉
我们就得到一个偏v偏x
由于这个角度比较小，所以tanα就是α
同样，β角最后得到的是偏u偏y
我们把这两个角度加起来，得到一个总的角度
也就是α加β
我们把它写成剪应变，这个剪应变就是γxy
基于位移场来表达它呢就是偏u偏y和偏v偏x
我们归纳一下上面所建立的几何方程
也就是说x方向的相对伸长量我们定义为εxx
y方向的相对伸长量定义为εyy
原来是个直角，变化以后的角度变化就是α加β
我们定义为剪应变γxy，它是偏u偏y，偏v偏x
那同样我们把它写成指标形式
可以写成εij等于二分之一ui偏导j加上uj偏导i
那这个地方注意两点
一个是ui偏导j这个偏导，这个下标
还是前面应力的偏导记号一样
我们这个逗号表示偏导
比如ui,j就表明是偏ui对xj进行求偏导
那么另外呢，由于我们写成指标形式的这个量
要和原来的分量形式要完全对应
那么我们可以验证一下
当i和j相等的时候
也就是i等于j分别等于1的时候
那我们就是εxx
当i等于j分别等于2的时候，就是εyy
当i和j不相等，分别等于1、2
也就是对应着x和y的时候
那这个时候我们得到的是两个的偏导呢
和上面γxy这个偏导呢差了二分之一
那么这个地方我们就要约定一下
我们由于是用的εij
所以说εij要等于二分之一γij
这样的话就可以把我们的指标形式和分量形式完全的对应起来
那么我们在构建几何方程的时候是基于位移分量
对于我们平面问题就是两个位移分量
就是uxy，vxy
我们分别构造了三个形变分量，也就是应变分量
也就是x方向的相对伸长量、y方向的相对伸长量
还有夹角的变化γxy
也就是说有两个独立的位移分量
我们给出了三个形变分量，也就是应变分量
这是由前面的公式可以进行计算的
那反过来，如果我们已知三个形变分量或者应变分量
能不能唯一的确定两个位移分量
从理论上来讲不可以唯一确定
除非这三个分量有一个相关条件
也就是说这三个分量里面
只有两个是独立的，另外一个是相关的
所以这相关条件我们就叫做变形协调条件
由这三个分量来推出变形协调条件
那么我们来看，εxx，εyy，还有γxy
分别都是基于u和v表达出来的
那么我们分别把εxx对y方向求两阶偏导数
把εyy对x方向求两阶偏导数
把它加到一块儿
同样呢把它写出关于u和v的偏导的表达（此处2应为偏导符号）
就可以把这个加起来的偏导呢
可以写成对γxy的两阶偏导，分别对偏x和偏y
这样的话我们构建的这个方程
就是三个应变分量必须满足的这么一个关系
也就是三个应变分量只有两个是独立的
其中要满足一个协调方程，也就是一个耦合方程
是一个非独立方程
那么这个方程呢，我们叫变形协调条件
那么变形协调条件从物理上来解释的话
就是说它只有满足变形协调方程所得到的位移
u和v才是真实的变形场的位移场，也就是u和v
如果是不满足变形协调方程得到的这个应变
所确定出来的u和v
它是不协调的，也就是不是真实的
那么在工程上或者在物理上我们来看
真实的协调的变形
那么这个是变形协调，这是变形前的情况
变形以后的协调方程应该是连续的
变化是无缝的
如果出现中间有间隙，有缝
甚至出现材料之间的搭接
那就说明三个应变分量出现了不协调
也就是说三个应变分量不满足变形协调方程
这样的话就不是真实的变形场
那么由它不能够确定出真实的位移场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
