[工学]有限元简介3.pptVIP

3.1 引 言 有限元法是由一位工程师提出来的数值方法，其基本思想就是将一个连续体划分为有限个在节点处相连接的小单元，然后利用各单元内所假设的近似函数通过分片逼近整个求解域上的未知场函数。这些单元可以是形状不同，但应以不同的连接方式组合，因此有限元法适合分析几何形状复杂的问题，包括结构工程、流体工程、电磁场等。 如下图(a)所示中间开有一圆孔的板，可以用三角形单元和四边形单元离散为图(b)所示。 商业有限元软件有丰富的单元，足以解决工程中所需单元的选取和应用。 3.2 轴力杆单元 1.单元形函数 如图所示长为 的一维轴力杆元，单元有二个节点i和j，节点的位移分量分别为 和 ，所以该单元具有二个自由度。这两个自由度属待定参数，它可用一个线性函数确定。单元内的位移场可假定为： (1) 式中 为待定常数。由于位移场u(x)在单元节点i和j应等于单元的节点位移 和 ，即 由上式可解得待定常数 为： 将上式代入式(1)，并注意到 ，则有： (2) 其中 (3) 称 为形函数或插值函数。 2.形函数的性质 1) 在节点i 处的值为1，在节点 j 的值为零， 亦类似。 即 2.形函数的性质 2) 形函数必须满足条件 。这一条件可由将式(3)相加证明，它表示当单元发生刚体位移 时单元内任意点的位移均应等于 ,将该条件代入式(2) 可得： 由于相邻单元公共节点的位移相等，因此轴力杆单元的位移函数在单元间连续，称这类单元为协调元。轴力杆单元的位移函数 是线性的，其一阶导数 在单元内为常数，因此在单元内连续。 式(2)可写为矩阵形式，即 (4) 称为插值函数矩阵，而 称为单元节点位移列阵。 3.应变矩阵与应力矩阵 轴力杆单元任意点的应变为: (5) 称为单元的应变矩阵或几何矩阵。 根据虎克定律，轴力杆单元内任意点的应力为 ， 于是由式(5)可得： (6) 式中 为单元材料的弹性模量， 称为单元应力矩阵。 4.单元刚度矩阵 轴力杆单元的应变能为： 式中 为轴力杆元的横截面积， 称为单元刚度矩阵。 轴力杆单元的刚度矩阵具体为： (7) 它还可写为分块形式： 5.单元等效节点载荷 因为有限元是对节点建立平衡方程进行求解，因此必须将作用于单元上的载荷等效到节点上。对于轴力杆单元，如果在其上作用有集中力F，单位长度上的外力q(x)其外力势能为： 将式(4)代入上式可得： 式中