当前课程知识点：有限元分析及应用 > 第3讲 针对复杂几何形状变形体的力学描述(1) > 3.6 平面问题的物理方程构建 > Video 3.6
下面我们构建平面问题的物理方程
这个物理方程实际上就是广义虎克定理
一维的问题我们已经做过实验了
我们做过一个单拉实验，在线弹性
我们可以得到应力和应变的关系
也就是说ε等于σ除上E
注意一下我们应力和应变都是两个下标εxx，σxx
它两个下标相同的时候有时候就简写成εx
σxx呢往往我们也可以简写是σx
不要看成一个下标就是一阶张量
实际上它就是两阶张量，也就是两个下标
那么我们要构建平面问题的物理方程呢
还是基于实验
那么这个实验呢应该是个广义的
除了沿拉伸方向构建的虎克定理以外
我们还要考虑另外一个方向的变化
那么同样我们也来做一个拉伸实验
那么主方向上的伸长，就是刚才提到的单拉实验
它就是这么一个线性关系
那么我们在主方向拉长的情况下
另外一个方向，一定要缩短
我们把这个缩短的比例，就是应变
相对伸长量的这个变化比例呢定义为泊松效应
它的两个之比，一个伸长一个缩短的之比呢
加个负号，我们把它定义为泊松比
那么有了主伸长和另外方向的缩短用泊松效应来描述以外
我们就可以给出平面问题在任意情况下的虎克定理
我们叫广义虎克定理
那我们看一下平面的受力情况
每一个侧面分别有两个力
一个是正应力，另外一个是剪应力
我们把这么几组力分解成x方向上的一组拉伸也叫σxx
y方向上的一组拉伸，σyy
还有四个剪切，我们分成这三组情况
那我们看看x方向的拉伸就相当于是个单拉实验
主方向相对伸长
那么伸长力和应变的关系，比例系数就是我们的弹性模量
这就是我们单拉主方向的虎克定理
另外它拉长以后，另外一个方向它就缩短
在y方向上的缩短量呢，我们用泊松系数来表达
它就是主方向拉伸量的泊松比相乘，再加个负号
因为它是缩短
那么我们把σxx拉伸引起的两个应变
就是一个伸长一个缩短，分别写到这儿
另外y方向应力σyy，在y方向的主伸长拉长
这是εyy，和σyy，这个弹性模量比例系数
另外它同样也有一个在x方向的泊松效应
因为y方向是主方向
所以x方向就是它的泊松效应的方向
那我们同样也可以用泊松比来写
得到εxx等于负的μσyy除上一个E
另外剪切，也就是τxy
这四个剪应力分量引起的变形的情况
我们考察一下
这个剪切呢，它只会引起夹角的变化
不会引起两个方向伸长的变化
所以说呢它只有剪切角的变化
我们都知道剪切角α加β我们用剪应变γxy来表达
所以它呢和剪应力之间也是满足一个线性关系
这个线性关系呢它的比例系数就是我们的剪切模量
那把这三种情况叠加
分别沿着x方向的伸长量、y方向的伸长量还有夹角的变化
分别找到对应的位置进行叠加
我们就可以得到合成以后的平面问题的本构方程
也就是广义虎克定理
它分别等于εxx，等于E分之一（σxx-μσyy）
那么我们把所得到的平面问题的物理方程进行一个汇总
就得到这样的方程
注意这里面的材料系数
E为弹性模量，G为剪切模量，μ为泊松比
这三个材料当中呢，只有两个是独立的
所以还要满足三个材料系数之间的一个关系
那么我们把这个物理方程写成一个逆形式
也就是说把应力写到前面，把应变写到这个关系式里面
那么同样我们可以得到逆形式的物理方程
那么前面我们得到的是应力和应变之间的关系
那一维情况最简单
二维情况我们得到的是一个方程组
那么按照指标形式来写
我们可以把应变写成εij
把应力呢也写成一个指标形式σkl，就是两个下标
两个下标不同
那么按照指标的运算规则
就两个二阶下标的指标形式的物理量之间
一定通过一个四阶下标的物理量建立它的关系
这就是我们指标规则里所得到的这么一个关系
那么下面我们来看看如果用指标形式来写
那么我们怎么来写出这个比例系数Cijkl
也就说是四阶张量，四个下标
那么我们怎么来确定这个Cijkl的这个系数呢
我们还是对应着前面所构建的平面问题的物理方程
我们看看第一项
第一项的物理方程呢是εxx，xx我们取1、1
那这样就是i=1、j=1
那么我们第一项这个系数呢是对应着σxx
σxx呢我们也是x对应着1
那这样就是k等于1，l等于1
那这样我们就可以得到ijkl分别等于1的情况
就是c1111它的对应的值就是E分之一
同样我们来写σyy的这一项
那同样xx还是1、1，也就是i等于1，j等于1
那么这个yy呢，y是表明是2，对应着x2轴
这样就是分别对应着k等于2，l等于2
那么这一项写出来就是i等于1，j等于1
k等于2，l等于2
对过来这个系数呢就是E分之负μ
所以写出来这个系数呢是C1122
就是负的E分之μ
那么其他的对应关系也是一样
我们分别来变下标
然后找出这两个下标分别对应的关系
来找出C这个对应下标的关系
那同样，这个物理方程的另外一个的逆形式
也是这样写的
那么σij要和εkl建立关系
同样有一个四阶张量的系数，我们叫Dijkl
同样我们用σxx，也就是i等于1、j等于1对过来
那么kl分别对应εxx也是等于1
1、1对应的是k等于1、l等于1
这样的话我们就对应着D1111
它对应的系数就是1减μ平方分之E
那同样，我们要做第二项
第二项对过来是εyy，也就是y呢对应的是2
那这样就是k等于2，l等于2
这样我们得到的系数就是i等于1，j等于1
k等于2，l等于2
D系数呢就是1122
它对过来这个系数就是1减μ平方分之μE
那么其它的参量也可以按照这种方式找出对应的关系
可以看出如果用指标形式来写
它是二阶张量之间通过一个四阶张量来进行转换
非常的复杂
要写出这个系数，要变这个下标都是非常的复杂
那这个时候我们就可以用Voigt的这个规则
把它这个阶次降低，来写出它的矩阵形式
那这样的话就是，对于σij等于Dijkl εkl
这么一个指标形式的物理方程
我们把ij，压缩成一个下标
也就是说通过Voigt这个规则
把一个两阶的矩阵变成一个列阵
那么同样对于εkl，把kl也变成一个下标，变成q
那这样D里面的ij就变成p，kl就变成q
那这样的话我们就得到了一个方程就是σp和εq之间的关系
那这样的话前面的系数原来是四阶的张量
就变成这么一个两个下标的降阶以后的张量
那么我们来写出四阶张量四个下标的
和两个下标的矩阵之间的对应关系
那么我们可以看，我们把pq分别变1、2、3
我们可以得到D一横11、12、13
D一横21、22、23
注意一下，这里面下标的第一个，表明是p的变化
第二个，表明是q的变化
那么p的变化又对着D矩阵前面的两个ij
这个q的变化又对着原来这个D下标的后面两个kl
也就是前面两个对着p，后面对着是q
那这样我们来做相应的变化
那么p取1的时候
也就是说我们σp，第一个叫σ1
也就是σxx、σyy、τxy
也就是说，σ1就对着σxx
也就是p取1的时候，就意味着ij分别取1、1
那么q呢，也要取1，对应着ε1
我们是这个εq
那么q要取1，也就对应着ε1，也就对应着εxx
也就是说，这个时候的kl也是取1、1
那么这个数就是D1111，它对着我们前面的D一杠11
这两个是一个对应关系
那我们看看第二项
我们的p还是取1
也就对应着σx，也就是说ij分别为1、1
然后q呢这个时候要取2，也就在这儿
这个时候q要取2
q取2就意味着εq等于2
εq等于2就对应着εyy，这个yy就对应着kl分别取2、2
那这个时候，我们所对应过来，D呢就是1122
这就是我们所对应的关系
那么我们第三项也是一样
我们的p还是取1，对应着σ1，也就说i等于1，j等于1
第二项呢就是q等于3
我们来看看εq等于3，就相当于是第三项
ε就是γ12，就对应着k等于1，l等于2
那这样的话就是第三项D，前面是11，后面就是12
那么我们就得到了这么一个对应的
四阶张量和二阶张量Voigt变换以后的这么一个对应关系
那么其它的，也都可以写出来
那么这个降阶的这九个量，分别对过来，就是这么具体的数字
这个是D1111，也就是这个降阶以后的D一杠11
这个呢对过来，是D1122，也对着我们这个D一杠12
那么其他也是这样的
那这就是用Voigt ij把它降为p一个下标和kl降成一个下标以后
就可以把原来很复杂的指标形式的本构关系
可以写成一个矩阵表达的这么一个形式
但注意一点，我们只是上工程上表达成矩阵
但是用指标形式可不能写成这个矩阵形式
因为指标形式一定是二阶张量的指标
这本来就是一个二阶矩阵
那么四阶的这个的话应该是一个四个下标变化的
那么这个矩阵是在纸面上是写不出来的
只能用下标来进行表达
-有限的单元 无限的能力
--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
