有限元 | 弹性支座

考虑一个有限元模型的势能泛函其中，是整体节点位移矩阵，是整体刚度矩阵。已知沿支座的自由度1的方向产生位移为，即。那么求势能泛函(1)的极值变成了在约束条件下的极值问题。用罚函数将有约束问题转化为无约束问题。引入一个很大的正参数，构造新的泛函令得到新的平衡方程这里可以看到，为处理，需要对以上方程进行修正，即将一个大数加到的第1个对角元上以及将加到上，式(4)的解即为节点位移列阵。节点1处的支反力为需要说明的是:这里所说的罚函数法只是一种近似的方法，最终求解的精度，特别是支反力的求解精度，取决于的选取。考察方程组(4)中的第一个方程(6)两边除以C从(7)式中我们可以发现，如果比刚度系数大得多时，那么有。是施加在支座处的外力(如果存在的话)，而且通常是一个很小的值。例1图1中，有一个外力在作用在杆中点，求杆的位移、应力和支座支反力(取)。▲图1在这个问题中，我们首先应该确定杆与墙之间是否发生接触。假设墙不存在，那么点B的位移为.从这个结果可以看出接触是存在的，因为边界条件发生了变化，即点B的位移是给定的，所以需要重新求解。建立两个单元的有限元模型如图1b所示，边界条件为和，结构刚度矩阵为整体节点力列阵为按照上述方法，取,则修正后的方程如下解得支反力为精确解为,。由罚函数法得到的结果有一个小的近似误差。例2▲图2如图2所示的四杆桁架结构。对于每个单元，给定,。节点2竖直下移。用先处理法建立如图3所示的有限元模型，得到的平衡方程组为▲图3考虑支座移动的影响，需要修正平衡方程组。在与指定位移值的自由度相对应的结构刚度矩阵对角线元素上加上一个大常数。一般来讲，可以取为修正前刚度矩阵的最大对角线元素的倍。此外，载荷列阵的对应位置还要加上力，其中为指定的位移值。在本例当中，对于自由度2来说，，因而载荷列阵的第2行上要加上一个大小为的力。修正后的有限元方程组为来源：数值分析与有限元编程