第三讲平面问题有限元分析.pptVIP

同理，第31页，共72页，星期日，2025年，2月5日例题：在均质、等厚的三角形单元ijm的ij边上作用有沿x方向按三角形分布的载荷，求移置后的结点载荷。取局部坐标s，在i点s=0，在j点s=L，L为ij边的长度。在ij边上，以局部坐标表示的插值函数为， ，，第32页，共72页，星期日，2025年，2月5日载荷为 第33页，共72页，星期日，2025年，2月5日例：总载荷的2/3移置到结点i，1/3移置到结点j，与原载荷同向第34页，共72页，星期日，2025年，2月5日第35页，共72页，星期日，2025年，2月5日2-5单元应力矩阵本节利用几何方程、物理方程，实现用结点位移表示单元的应变和单元的应力。用结点位移表示单元的应变的表达式为 ，[B]矩阵称为几何矩阵。 第36页，共72页，星期日，2025年，2月5日2-5单元应力矩阵由物理方程，可以得到单元的应力表达式为应力矩阵 第37页，共72页，星期日，2025年，2月5日2-6单元刚度矩阵讨论单元内部的应力与单元的结点力的关系，导出用结点位移表示结点力的表达式。由应力推算结点力，需要利用平衡方程。第一章中已经用虚功方程表示出平衡方程。 第38页，共72页，星期日，2025年，2月5日2-6单元刚度矩阵考虑上图三角形单元的实际受力，结点力和内部应力为： 任意虚设位移，结点位移与内部应变为第39页，共72页，星期日，2025年，2月5日2-6单元刚度矩阵令实际受力状态在虚设位移上作虚功，外力虚功为第40页，共72页，星期日，2025年，2月5日2-6单元刚度矩阵计算内力虚功时，从弹性体中截取微小矩形，边长为dx和dy，厚度为t，图示微小矩形的实际应力和虚设变形。第41页，共72页，星期日，2025年，2月5日2-6单元刚度矩阵微小矩形的内力虚功为整个弹性体的内力虚功为第42页，共72页，星期日，2025年，2月5日2-6单元刚度矩阵根据虚功原理，得这就是弹性平面问题的虚功方程，实质是外力与应力之间的平衡方程。虚应变可以由结点虚位移求出：代入虚功方程第43页，共72页，星期日，2025年，2月5日2-6单元刚度矩阵接上式，将应力用结点位移表示出有令则建立了单元的结点力与结点位移之间的关系，称为单元刚度矩阵。它是6*6矩阵，其元素表示该单元的各结点沿坐标方向发生单位位移时引起的结点力，它决定于该单元的形状、大小、方位和弹性常数，而与单元的位置无关，即不随单元或坐标轴的平行移动而改变。第44页，共72页，星期日，2025年，2月5日2-6单元刚度矩阵由于[D]中元素是常量，而在线性位移模式下，[B]中的元素也是常量，且因此可以进一步得出平面应力问题和平面应变问题中的单元刚度矩阵。第45页，共72页，星期日，2025年，2月5日2-7单元刚度矩阵的物理意义及其性质已经求出了下列关系第46页，共72页，星期日，2025年，2月5日2-7单元刚度矩阵的物理意义及其性质结点力和结点位移的关系：(以简单平面桁架为例)平面问题中，离散化的单元组合体极为相似，单元组合体在结点载荷的作用下，结点对单元、单元对结点都有作用力与反作用力存在，大小相等方向相反，统称为结点力。结点力和结点位移的关系前面已经求出：第47页，共72页，星期日，2025年，2月5日2-7单元刚度矩阵的物理意义及其性质单元刚度矩阵的物理意义：将写成分块矩阵写成普通方程其中表示结点S(S=i,j,m)产生单位位移时，在结点r(r=i,j,m)上所需要施加的结点力的大小。第48页，共72页，星期日，2025年，2月5日2-7单元刚度矩阵的物理意义及其性质单元刚度矩阵的物理意义：将结点力列矩阵与结点位移列矩阵均展开成(6*1)阶列矩阵，单元刚度矩阵相应地展开成(6*6)阶方阵：元素K的脚码，标有“-”的表示水平方向，没有标“-”的表