有限元基础理论教程lecture05.pptVIP

2.4 单元刚度矩阵 2.5单元刚度矩阵的性质与物理意义 假设单元的结点位移如下： 单元刚度矩阵的性质： 1）对称性 2）奇异性 3）对角线上的主元恒正 2.6整体分析 整体分析包括以下4个步骤： 1）建立整体刚度矩阵， 2）根据约束条件修改整体刚度矩阵， 3）解方程组，求出结点的位移， 4）根据结点位移，求出单元的应变和应力。 建立整体刚度矩阵的基本方法是刚度集成法，即整体刚度矩阵是单元刚度矩阵的集成。 编程实现整体刚度矩阵集成 NEC(i,1)，表示单元(i)中I结点的整体编号； NEC(i,2)，表示单元(i)中J结点的整体编号； NEC(i,3)，表示单元(i)中M结点的整体编号。 分块矩阵[Krs]中的系数按以下地址加入整体刚度矩阵： 在编制有限元程序时，可以用一个二维数组来保存单元结点的局部编号与该结点在整体结构中的总体编号之间的关系。假定最大单元数目为NELEM，用Fortran语言可以定义以下数组， DIMENSION NEC(NELEM,3) 数组的行与单元对应，数组的第1列为单元中I结点的整体编号，第2列为单元中J结点的整体编号，第一列为单元中M结点的整体编号。例如， NEC(2,1)=5，表示单元(2)中I结点的整体编号为5； NEC(2,2)=2，表示单元(2)中J结点的整体编号为2； NEC(2,3)=4，表示单元(2)中M结点的整体编号为4。 假定弹性体总共划分为m个单元，n个结点。 * 2）体力的移置 令单元所受的均匀分布体力为， (2-17) 虚功相等， 例题2.4、设有均质等厚的三角形单元ijm，受到沿y方向的均布载荷qy的作用。求均布体力移置到各结点的载荷。 三角形中的一点P可以用子三角形面积定义的自然坐标来确定。面积坐标定义为， 点P表示为， 3）分布面力的移置 设在单元的边上分布有面力， 虚功相等， 面积坐标在三角形全面积上的积分为 例题2.5、在均质、等厚的三角形单元ijm的ij边上作用有沿x方向按三角形分布的载荷，求移置后的结点载荷。 取局部坐标s，在i点s=0，在j点s=l，L为ij边的长度。在ij边上，以局部坐标表示的插值函数为， 载荷为 设ij边的长度为L，先把分布面力等效为作用在距 i结点2/3L处P点的集中力，再移置到结点上。 由几何方程可以得到单元的应变表达式， (2-19) 记为 [B]矩阵称为几何矩阵。 [B]矩阵可以表示为分块矩阵的形式， (2-20) 由物理方程，可以得到单元的应力表达式， [D]称为弹性矩阵，对于平面应力问题， 应力矩阵 应用虚功原理可以建立单元结点位移与结点力的关系矩阵，即单元刚度矩阵。 虚功原理：在外力作用下处于平衡状态的弹性体，如果发生了虚位移，则所有外力在虚位移上做的虚功等于内应力在虚应变上做的虚功。 单元的结点力 单元的虚应变 单元结点力虚功 单元的内力虚功 用结点位移表示应力 单元刚度矩阵 在3结点等厚三角形单元中[B]和[D]的分量均为常量，则单元刚度矩阵可以表示为， （2-24） （2-25） 单元刚度矩阵表示为分块矩阵： 单元刚度矩阵的分块为 例2.6、属于平面应力问题的弹性体被划分成3个单元、5个结点。结点坐标如下：1(0,2a)，2(0,a)， 3(a,a)，4(0,0)，5(a,0)。 单元结点的局部编号顺序如下， e1(1,2,3); e2(2,4,5); e3(2,5,3) 试求单元1的单元刚度矩阵。 单元1结点的局部编号如下， i - 1, j - 2, m - 3 单元1的单元刚度矩阵： 得到结点力如下： 表示s结点在水平方向产生单位位移时，在结点r的垂直方向上需要施加的结点力。 表示s结点在水平方向产生单位位移时，在结点r的水平方向上需要施加的结点力。 选择不同的单元结点位移，可以得到单元刚度矩阵中每个元素的物理含义： 表示s结点在垂直方向产生单位位移时，在结点r的水平方向上需要施加的结点力。 表示s结点在垂直方向产生单位位移时，在结点r的垂直方向上需要施加的结点力。 单元刚度矩阵中每个元素都可以理解为刚度系数，即在结点发生单位位移时所需要施加的结点力。 首先证明单元刚度矩阵的对称性，即 假定单元产生了x方向的单位位移的刚体移动， 再证明单元刚度矩阵是奇异的，即 对应的单元结点力为零， 可以得到，在单元刚度矩阵中1，3，5列中对应行的系数相加为零，由行列式的性质可知，单元刚度矩阵是奇异的。 如何说明单元刚度矩阵对角线上的元素恒正？ 刚度集成法的物理意义 结构中的结点力是相关单元结点力的叠加，整体刚度矩阵的元素是相关单元的单元刚度矩阵元素的集成。 刚度矩阵集成的规则： 1）将单元刚度矩阵中的每个分块放到在整体刚度矩阵中的对应位置上，得到单元的扩大刚度矩