02有限元法的基本概念-引例.pptVIP

第二章 一维问题有限元法——引例 设有一根受自重作用的等截面直杆，上端固定，下端自由，如图所示。单位杆长的重力为q，杆长为L，横截面积为A，材料弹性模量为E。试求该直杆各横截面上的应力。 材料力学的精确解答： 用有限元法求解。 第一步：用若干个分点将该直杆分割成为许多个有限长度的小微段（这里各小段的长度不一定相等），并将每个小段的重力等效地移置到分点上去。分点称为节点，小段称为单元，移置到节点后的荷载称为节点荷载。本例中分成3个小段。 由材料力学解答可知，无论怎样分割该直杆，每个单元的位移都是x的二次函数。但若节点较多，单元较短，对每个单元来讲，我们就可以用线性函数近似地描述实际的单元位移。 我们用eij表示任一单元，其中e为单元编号①、②、③，ij为其两端点的节点编号12、23、34。 设线性单元位移函数为 其中，?1、?2为待定常数，可由该单元节点的位移值ui、uj来确定。单元位移函数也适于单元节点，则 解得 代回到前式中，并记{f} = u，得 其中，(xj-x)/(xj-xi)和(x-xi)/(xj-xi)与x成线性关系，反映了单元的位移形态，称为形函数，并令 Ni = (xj-x)/(xj-xi)，Nj = (x-xi)/(xj-xi) 记形函数矩阵为 [N]=[Ni Nj] 节点位