第三章平面问题的有限单元法.pptVIP

即：其中：为形函数在集中力作用点处的值。平面问题的有限单元法3-27第61页，共121页，星期日，2025年，2月5日y0xijm设在单元ijk上受有分布体力。取微元体，则此微元体可看成一个集中力：体力等效结点力：由（3-27）式，并对三角形单元体积积分，即可得到体力位置到结点上的等效结点力列阵：平面问题的有限单元法2.体积力的等效载荷列阵{R}第62页，共121页，星期日，2025年，2月5日3、面力的移置：设在单元的某一个边界上作用有分布的面力,单位面积上的面力为{P}＝[XY]，在此边界上取微面积tds，将微面积上的面力视为集中载荷，利用(3—27)式，对整个边界面积分，得到:第63页，共121页，星期日，2025年，2月5日4、特殊载荷的移置：y0xijm已知在ij边受有面力q，则移置到i、j结点上的等效结点力为：y0xijm平面问题的有限单元法当某一边上有三角形分布的面力时，可由刚体静力等效直接写出第64页，共121页，星期日，2025年，2月5日如果任意三角形单元ijk的重心c上受有自重则按刚体静力等效原理可把W直接移置到i，j，m三个结点上而组成：第65页，共121页，星期日，2025年，2月5日等效结点力、载荷列阵小结第66页，共121页，星期日，2025年，2月5日例:图示三角形单元的ijm的jm边上作用有如图所示线性均匀分布面载荷，三角形单元厚度为t，求结点载荷向量。解：运用叠加原理三角形单元jm边长度为：Xm(5,6)j(6,3)i(2,2)q1q2OY第67页，共121页，星期日，2025年，2月5日习题1、设有三角形单元ijm受自重作用，容重为Pg。求出其等效结点载荷。2．设有：三角形单元ijm的ij边受线性分布的法向压力作用，ij两点的压力分别为Pi和Pj试求其等效结点载荷。第68页，共121页，星期日，2025年，2月5日3.4整体刚度矩阵平面问题的有限单元法讨论了单元的力学特性之后，就可转入结构的整体分析。假设弹性体被划分为N个单元和n个结点，对每个单元按前述方法进行分析计算，便可得到N组形如（3-21）式的方程。将这些方程集合起来，就可得到表征整个弹性体的平衡关系式。为此，我们先引入整个弹性体的结点位移列阵{?}2n×1，它是由各结点位移按结点号码以从小到大的顺序排列组成，即其中子矩阵(j)(i=1,2,…,n)(k)是结点i的位移分量。第69页，共121页，星期日，2025年，2月5日平面问题的有限单元法继而再引入整个弹性体的载荷列阵{R}2n×1，它是移置到结点上的等效结点载荷依结点号码从小到大的顺序排列组成，即(l)其中子矩阵(i=1,2,…,n)(m)是结点i上的等效结点载荷。第70页，共121页，星期日，2025年，2月5日平面问题的有限单元法现将各单元的结点力列阵{R}e6×1加以扩充，使之成为2n×1阶列阵其中，子矩阵(n)(i,j,m轮换)(o)是单元结点i上的等效结点力。(n)式中的省略号处的元素均为零，矩阵号上面的i,j,m表示在分块矩阵意义下Ri所占的列的位置。此处假定了i,j,m的次序也是从小到大排列的、并且与结点号第71页，共121页，星期日，2025年，2月5日平面问题的有限单元法码的排序一致。各单元的结点力列阵经过这样的扩充之后就可以进行相加，把全部单元的结点力列阵叠加在一起，便可得到(l)式所表示的弹性体的载荷列阵，即这是由于相邻单元公共边内力引起的等效结点力，在叠加过程中必然会全部相互抵消，所以只剩下载荷所引起的等效结点力。同样，将单元刚度矩阵的六阶方阵[k]加以扩充，使之成为2n阶的方阵(p)返回第72页，共121页，星期日，2025年，2月5日