6第六章用有限单元法解平面问题2011s讲述.pptVIP

以平面问题为例，利用坐标变换可以将任意形状的四边形映射到自然坐标系下的边长为2的正方形，如图所示，如果坐标变换和位移函数插值采用相同的结点，并采用相同的插值函数，则这种变换称为等参变换。采用等参变换的单元称为等参单元。如果坐标变换结点多于位移函数插值的结点数，则称为超参变换，反之称为次参变换。 第五节 单元的结点力列阵与劲度矩阵 为了推导单元的节点力和节点位移之间的关系，可应用虚位移原理对单元e进行分析。单元e是在等效节点力的作用下处于平衡的，而这种结点力可采用列阵表示为 假设在单元的三个节点i、j、m 的虚位移为 单元内的虚应变{? *}为 于是，作用在单元体上的外力在虚位移上所做的功可写为 而单元内的应力在虚应变上所做的功为 根据虚位移原理，得到单元的虚功方程为 注意到虚位移是任意的，所以等式两边与相乘的项应该相等，即得 推导可得 即 记 则有 上式就是表征单元的节点力和节点位移之间关系的刚度方程，[k]e就是单元劲度矩阵，也称单元刚度矩阵。如果单元的材料是均质的，那么矩阵[D] 中的元素就是常量，并且对于三角形常应变单元，[B]矩阵中的元素也是常量。当单元的厚度也是常量时， [k]e =[B]T [D][B] A 与前面讨论过的情况类似，单元刚度矩阵[k]中任一列的元素分别等于该单元的某个节点沿坐标方向发生单位位移时，在各节点上所引起的节点力。单元的刚度取决于单元的大小、方向和弹性常数，而与单元的位置无关，即不随单元或坐标轴的平行移动而改变。 平面应力问题中三角形单元的刚度矩阵的分块形式为： 其中 ( r = i、j、m；s = i、j、m ) 对于平面应变问题，只要将上式中的E、?分别换成E / 1-? 2 和? / 1-? 即可。于是 ( r = i、j、m；s = i、j、m ) 第六节 载荷向结点移置 单元的结点载荷列阵 为了简化每个单元的受力情况，便于分析，我们把单元所受的外力载荷都向结点移置而成为结点载荷。这种移置必须按照静力等效的原则进行。对于变形体，包括弹性体内，所谓静力等效，是指原载荷与结点载荷在任意虚位移上的虚功都相等。 设单元 在坐标为的任意一点M，在单位厚度上受有集中载荷 ，分力为 ， 。将此集中力移置到单元的结点处，转换为结点载荷，并用单元结点载荷列阵表示为， 设单元的各结点发生了虚位移 由前面的位移模式，得到相应于集中力 的作用点（x,y）的虚位移为 按静力等效原则，结点载荷在结点虚位移上的虚功，应当等于原载荷集中力在其作用点的虚位移上的虚功，即 将 代入得 由于虚位移可以是任意的，从而矩阵 也是任意的，于是得 即， 其中， 应当是它们在M点的函数值。 设上述单元受有分布体力 ， 可将微分体积 上的体力当作集中力， 则 即 第七节结构的整体分析 结点平衡方程组 讨论了单元的力学特性之后，就可转入结构的整体分析。假设弹性体被划分为N个单元和n个节点，对每个单元按前述方法进行分析计算，便可得到N组形如 的方程。将这些方程集合起来，就可得到表征整个弹性体的平衡关系式。为此，我们先引入整个弹性体的节点位移列阵 {?}2n×1 ，它是由各节点位移按节点号码以从小到大的顺序排列组成，即 其中子矩阵 (i =1,2, …, n ) 是节点i的位移分量。 继而再引入整个弹性体的载荷列阵{F}2n×1 ，它是移置到节点上的等效节点载荷依节点号码从小到大的顺序排列组成，即 其中子矩阵 (i =1,2, …, n ) 是节点i上的等效节点载荷。 同样，将六阶方阵[k]加以扩充，使之成为2n阶的方阵 不难看出2×2阶子矩阵[ki j ]将处于上式中的第i双行、第j双列中。 考虑到[k]扩充以后，除了对应的i, j, m 双行和双列 上的九个子矩阵之外，其余元素均为零，故式 中的单元位移列阵{?}e2n×1 便可用整体的位移列阵 {?}2n×1 来替代。这样 式可改写为 把上式对N个单元进行求