有限元法及应用总结.pptVIP

位移函数在相邻单元的公共边界上必须协调。对一般单元而言，协调性是指相邻单元在公共节点处有相同的位移，而且沿单元边界也有相同的位移，也就是说，要保证不发生单元的相互脱离开裂和相互侵入重叠。要做到这一点，就要求函数在公共边界上能由公共节点的函数值唯一确定。对一般单元，协调性保证了相邻单元边界位移的连续性。但是，在板壳的相邻单元之间，还要求位移的一阶导数连续，只有这样，才能保证结构的应变能是有界量。有限元法的收敛性概念与收敛条件（续）有限元法的收敛性概念与收敛条件（续）总的说来，协调性是指在相邻单元的公共边界上满足连续性条件。前三条又叫完备性条件，满足完备条件的单元叫完备单元；第四条是协调性要求，满足协调性的单元叫协调单元；否则称为非协调单元。完备性要求是收敛的必要条件，四条全部满足，构成收敛的充分必要条件。在实际应用中，要使选择的位移函数全部满足完备性和协调性要求是比较困难的，在某些情况下可以放松对协调性的要求。有限元法的收敛性概念与收敛条件（续）需要指出的是，有时非协调单元比与它对应的协调单元还要好，其原因在于近似解的性质。假定位移函数就相当于给单元施加了约束条件，使单元变形服从所加约束，这样的替代结构比真实结构更刚一些。但是，这种近似结构由于允许单元分离、重叠，使单元的刚度变软了，或者形成了（例如板单元在单元之间的绕度连续，而转角不连续时，刚节点变为铰接点）对于非协调单元，上述两种影响有误差相消的可能，因此利用非协调单元有时也会得到很好的结果。在工程实践中，非协调元必须通过“小片试验后”才能使用。应力的单元平均或节点平均处理方法?最简单的处理应力结果的方法是取相邻单元或围绕节点各单元应力的平均值。取相邻单元应力的平均值这种方法最常用于3节点三角形单元中。这种最简单而又相当实用的单元得到的应力解在单元内是常数。可以将其看作是单元内应力的平均值，或是单元形心处的应力。由于应力近似解总是在精确解上下振荡，可以取相邻单元应力的平均值作为此两个单元合成的较大四边形单元形心处的应力。如2单元的情况下，取平均应力可以采用算术平均，即平均应力=（单元1的应力+单元2的应力）/2。也可以采用精确一些的面积加权平均，即平均应力=[单元1应力×单元1的面积+单元2应力×单元2面积]/（单元1面积+单元2面积）当相邻两单元面积相差不大时，两者的结果基本相同。在单元划分时应避免相邻两单元的面积相差太多，从而使求解的误差相近。应力的单元平均或节点平均处理方法?（续）0102应力的单元平均或节点平均处理方法?（续）一般而言，3节点三角形单元的最佳应力点是单元的中心点，此点的应力具有1阶的精度。取围绕节点各单元应力的平均值首先计算围绕该节点（i）周围的相关单元在该节点出的应力值，然后以他们的平均值作为该节点的最后应力值，即其中，1~m是围绕在i节点周围的全部单元。取平均值时也可进行面积加权。11.有限元法求解问题的基本步骤对整个结构进行离散化，将其分割成若干个单元，单元间彼此通过节点相连；1.结构离散化是由单元节点位移量求单元节点力向量的转移矩阵，其关系式为2.求出各单元的刚度矩阵0201集成总体刚度矩阵[K]并写出总体平衡方程：总体刚度矩阵[K]是由整体节点位移向量求整体节点力向量的转移矩阵，其关系式为，此即为总体平衡方程。引入支撑条件，求出各节点的位移节点的支撑条件有两种：一种是节点n沿某个方向的位移为零，另一种是节点n沿某个方向的位移为一给定值。求出各单元内的应力和应变。12.单元刚度矩阵的特性单元刚度矩阵无论在局部坐标系中还是在整体坐标系中都具有相同的三个特性：称性由材料力学中的位移互等定理可知，对一个构件，作用在点j的力引起点i的绕度等于有同样大小而作用于点i的力引起的点j的绕度，即，表明单元刚度矩阵是一个对称矩阵。logo单元刚度矩阵的特性（续）2）异性无逆阵的矩阵就叫做奇异矩阵，其行列式的值为0，即，这一点可以从例题直接得到验证。其物理意义是引入支撑条件之前，单元可平移。12.单元刚度矩阵的特性（续）块性有前面所讲的内容可以看出，矩阵可以用虚线分成四块，因此可写成如下的分块形式，式中——局部坐标系中单元(e)按局部码标记的节点m、n之间的刚度子矩阵。*13.求图中所示刚架中各单元在整体坐标系中的单元刚度矩阵(1)。设两杆的长度与截面尺寸彼此相等。（空心杆）其中，L=200cm,D=5cm,d=4cm,E=2×107N/cm2。*13：求图所示结构的节点位移向量。已知节点1处承受外载