ansys第六讲.pptVIP

假定该系统中节点位移u1 和u2分别被指定为u1 ?1，u2 ?2，当引入这些节点的已知位移之后，方程 a 就变成 然后，就用这组维数不变的方程来求解所有的节点位移。显然，其解答仍为原方程 a 的解答。 保持方程组为2n×2n系统，仅对[K]和 R 进行修正。例如，若指定节点i在方向y的位移为vi ，则令[K]中的元素k2i, 2i 为1，而第2i行和第2i列的其余元素都为零。 R 中的第2i个元素则用位移vi 的已知值代入， R 中的其它各行元素均减去已知节点位移的指定值和原来[K]中该行的相应列元素的乘积。 * 将[K]中与指定的节点位移有关的主对角元素乘上一个大数，如1015，同时将 R 中的对应元素换成指定的节点位移值与该大数的乘积。实际上，这种方法就是使[K]中相应行的修正项远大于非修正项。 若把此方法用于上面的例子，则方程 a 就变成 事实上，该方程组的第一个方程为 返回 * 图1所示为一厚度t 1cm的均质正方形薄板，上 下受均匀拉力q 106N/m，材料弹性模量为E，泊松比 ，不记自重，试用有限元法求其应力分量。 1 2 3 4 2 1 x 图 2 y 2m y x q 106N/m 图 1 例1 9、计算实例 * 解： １.力学模型的确定 ２.结构离散 　　由于此结构长、宽远大于厚度，而载荷作用于板平面内，且沿板厚均匀分布，故可按平面应力问题处理，考虑到结构和载荷的对称性，可取结构的1/4来研究。 　　该1/4结构被离散为两个三角形单元，节点编号，单元划分及取坐标如图2所示，其各节点的坐标值见表1。 ３.求单元的刚度矩阵 计算单元的节点坐标差及单元面积 　　单元１（i、j、m　　1，2，3） 节点 坐标 １ ２ ３ ４ x y ０ ０ １ ０ １ １ ０ １ 表1 * 计算各单元的刚度矩阵 先计算用到的常数 代入可得: * 所以单元1的刚度矩阵为: 1 2 3 1 2 3 * 由于单元2若按341对应单元1的123排码时,则这两个单元刚度矩阵内容完全一样，故有: 3 4 1 3 4 1 * 组集整体刚度矩阵 由于[Krs] [Ksr]T，又单元1和单元2的节点号按123对应341，则可得： 按刚度集成法可得整体刚度矩阵为： * * 所以组集的整体刚度矩阵为： * 5. 引入约束条件，修改刚度方程并求解 根据约束条件：u1 v1 0；v2 0；u4 0和等效节点力列阵： ，并代入刚度方程： ，划去[K]中与0位移相对应的1，2，4，7的行 和列，则刚度方程变为： 求解上面方程组可得出节点位移为： 所以 * 先求出各单元的应力矩阵[S]1、[S]2，然后再求得各单元的应力分量： 6. 计算各单元应力矩阵，求出各单元应力 单元应力可看作是单元形心处的应力值。 * 例2 图中所示为一平面应力问题离散化以后的结构图，其中图（a）为离散化后的总体结构，图（b）为单元1，2，3，4的结构，图（c）为单元3的结构。用有限元法计算节点位移、单元应变及单元应力（为简便起见，取泊松比 ，单元厚度t 1）。 x y 1 2 3 4 6 5 1 2 3 4 a 3 i j m a a a 1,2,4 i j m 图 计算实例2的结构图 * 　　首先求确定各单元刚度所需的系数 及面积A，对于单元1，2，4有： 解： 对于单元3有: * 其次，求出各单元的单元刚度矩阵。对于1，2，4单元，其单元刚度矩阵为： i j m i j m * 各单元的节点编号与总体结构的总编号之间的对应关系见表2。 对于单元3，其单元刚度矩阵为： i j m i j m * 各单元节点号与总体节点号对应表 单元号 1 2 3 4 节点号 节 点 总 编 号 I 1 2 2 3 j 2 4 5 5 m 3 5 3 6 表2 * 将各单元刚度矩阵按节点总数及相应的节点号关系扩充成12*12矩阵,分别如下: * * * * 将扩充后的各单元刚度矩阵相加,得总体刚度矩阵K，即： * 所以结构总方程为： 其中 考虑到边界条件： * 用对角元乘大数法消除奇异性后的结构总体方程为： * 由以上方程解得的各节点的位移为： * 然后将相应的节点位移代入公式，可分别求得各单元的应 变和应力。 对于单元1： * 对于单元2： * 对于单元3： * 对于单元4： * 思考题 简述有限元方法进行问题求解的基本步骤？ 三角形常应变单元位移模式和各矩阵的公式、维数及意义？ 形函数的性质？ 半带宽的求解？ 如何使用刚体静力等效原理进行单元载荷的移置？ 如何进行边界条件的处理及整体刚度矩阵的修正？ 例题的求解步骤和注意事项？ * * * 图 6 a * 图 6 b 半带宽B （相邻节点号的最大差值D+1）*2 * 若第r双行的第一