线性三角形单元.pptVIP

刚度矩阵*(1)刚度矩阵中每个元素有明确的物理意义。例如，Kij表示当节点位移中第j个元素为1(dj=1)其余元素为零时，引起的单元力中的第i个节点力fi。把平衡方程写开主对角线上元素Kii(i=1,Ns)恒为正值：位移和作用力同向刚度矩阵*(2)K的每一行或每一列元素之和为零以上式中第i行为例:f2i-1=0f2i=0f2j-1=0f2j=0f2m-1=0f2m=0ijmxyijm11当所有节点沿x向或y向都产生单位位移时，单元作平动运动，无应变，也无应力，因而单元结点力为零（不含初应力）。所以有即，K的每一行元素之和为零。由于对称性，每一列元素之和也为零。刚度矩阵刚度矩阵是在系统处于平衡状态的前提下得出的。作用在它上面的外力必定是平衡力系。然而，研究系统平衡时没有引入约束。承受平衡力系作用的无约束系统，其变形是确定的，但位移不是确定的。所以出现性质（2）中的“平动问题”，即可以发生任意的刚体运动。从数学上讲，系统平衡方程的解不是唯一的或不能确定的。由此，系统刚度矩阵一定是奇异的。单元刚度矩阵也一定是奇异的。(3)系统刚度矩阵是奇异矩阵（即K的行列式为零）01系统的节点力和节点位移成线性关系是基于弹性理论的结果。(4)系统刚度矩阵是常量矩阵02位移边界条件的处理*系统刚度矩阵是奇异矩阵，其物理原因是结构缺少刚性位移的约束，实际的工程结构都受有足够的支承约束，排除了发生任何刚体位移的可能性，因此，必须引入位移约束。有限元中，位移约束都设置在节点处。这里，只讨论刚性约束情况，即被约束的位移分量为零。设讨论的结构有Nn个节点，每个节点有ndf个自由度。则系统的总自由度为Ns，且节点总位移列向量d中共包含Ns个分量为了引入位移约束，把节点总位移列向量d分成两部分。一部分是不受约束的位移分量，记为df。另一部分是受刚性约束的位移分量，记为dr。不失一般性，设1～N号位移分量是不受约束的；N+1～N+Nr共Nr个分量是受刚性约束的。即：位移边界条件的处理显然不受约束的节点位移的总数N为N=Ns-Nr对方程中的刚度矩阵K和节点荷载向量列阵f也作相应分割，则得到式中，ff是已知力边界，fr是约束反力。位移边界条件的处理12按矩阵乘法规则得3每个受刚性支承约束的位移分量都等于零，即4从而得到位移边界条件的处理01Kff——引入约束后的约化的系统刚度矩阵。这是一个非奇异矩阵，它的逆矩阵Kff-1是存在的。02引入约束后的约化的系统平衡方程03在分析计算时，从无约束的系统刚度矩阵K中删去与受约束位移号对应的行和列，再将矩阵压缩排列成N×N阶方阵，即为约化后的结构刚度矩阵。04位移边界条件的处理位移边界条件的处理*00010置一法显然乘大数法显然位移边界条件的处理节点位移和约束反力*通过求解平衡方程即可解出全部未知的节点位移：约束反力把解出的d代入未经修改的平衡方程，即可得到约束反力：关于上述方程的解算方法，一般不采用求逆的方法求解，而是直接采用高斯消元法等求解线性方程组的方法求解求解。施加边界条件后，得到修改后的平衡方程(未约化的)(约化的)或节点位移或单元应变和应力*根据三角形节点的位移，求出单元应力应变为ijk如何求系统应变能和节点应力？三角形单元*引言杆梁结构：由于有自然的连接关系，可以凭一种直觉将其进行自然的离散。连续体：它的内部没有自然的连接节点，必须完全通过人工的方法进行离散。三维问题平面问题平面应力平面应变平面问题平面应力平面应变离散三节点平面三角形单元*节点1的位移节点2的位移节点3的位移三节点三角形单元的位移函数可假设为：“位移函数”也称“位移模式”，是单元内部位移变化的数学表达式，是坐标的函数。有限元分析必须事先给出（设定）位移函数。一般而论，位移函数选取会影响甚至严重影响计算结果的精度。弹性力学中，恰当选取位移函数不是一件容易的事情。有限单元法中当单元划分得足够小时，把位移函数设定为简单的多项式也可得到相当精确的结果。这正是有限单元法具有的重要优势之一。引入位移函数的概念：平面三角形单元*显然，三角形三个节点的的位移可由下列方程给出，在各节点上的水平位移方程为：u1=?1+?2x1+?3y1u2=?1+?2x2+?3y2