P3-有限单元法-3.pptVIP

* * 回顾有限单元法的分析步骤：三大步。第二大步单元特性分析包括四小步。 * 四步骤的依据：收敛准则、几何方程、物理本构关系、虚功原理。 基础：单元位移模式 目标：节点位移与节点力之间的关系 核心：单元刚度矩阵 * * 位移模式如何选取或确定？ 位移模式的选取必须具备两个条件： (1)从物理上应满足收敛条件； (2)从数学上既能较好地反映单元实际位移，又便于数学处理。 * 收敛：随着单元的不断细化，其极值收敛于真实解。有限单元法作为一种数值方法，其解答应随着网格的逐步细分，收敛于问题的精确解。 1. 收敛：0.99999……＝1 ？！ 2.实践证明，在给定载荷之下，有限单元法计算模型的变形要比实际变形小，即当单元网格逐步细化时，位移的近似解将由下方收敛于精确解，得到真实解的下界。这是因为人为地给定位移模式，相当于对结构多加了约束，限制了变形、提高结构的刚度。 3.收敛与位移模式的选取关系重大。单元形状确定之后，位移模式的选取拭关键，因为应力、应变及刚度矩阵的建立都依赖于位移模式。 条件1：当节点位移是由某个刚体位移引起时，单元内不会有应力和应变。引入刚体位移，不但具有描述单元自身变形的能力，也具备了描述由于其它单元变形而连带引起的刚体位移的能力，这种能力是必要的，因为多数单元的刚体位移成分很大，甚至是主要的。例如：杂靠近悬臂梁的自由端处，单元应变很小，其位移成分主要是由于其它单元变形而引起的刚体位移； 条件2：每个单元的实际应变一般包含两个部分，一部分是与该单元各点的位置坐标有关的变应变，另一部分是与位置坐标无关的常应变。当单元尺寸较小时，单元应变趋于均匀，常应变就成为应变的主要部分。如位移模式中不包含常应变，就不可能正确反映单元的位移形态。 条件3：单元内位移是连续的，单元间变形是协调的，不会出现两部分相互叠合或相互脱离。单元内位移的连续性是容易保证的，只要将位移函数取为坐标的单值连续函数即可。而在单元间保证连续性较困难。可以证明：当单元界面上的位移仅由该界面上的节点位移确定时，单元间的变形协调性要求即可满足。 4.满足条件1和条件2的单元称为完备单元；满足条件3的单元称为协调单元或保续单元。同时满足三个条件的单元为完备协调单元。 * 普遍选择多项式作为位移模式，原因有二： 任一闭区间上的光滑函数，都可以用多项式逼近，项数愈多、愈精确 ； 多项式便于数学处理(微分和积分)，而且很容易满足收敛条件 。 位移模式应该与局部坐标系的方位无关：几何各向同性，根据巴斯卡三角形来选择二维多项式各项即可满足要求。 * 位移模式应该与局部坐标系的方位无关：几何各向同性，根据巴斯卡三角形来选择二维多项式各项即可满足要求。 * 高次单元 * 三节点三角形单元有三个节点，每个节点有两个位移分量，单元自由度数为6。 * 1.两个方向的位移，分别写成该点坐标的函数。 2.既然假设的是单元内任意点的位移，所以该位移模式对单元节点也应该适用，代入即可获得关于六个待定系数的方程组。 3.从中可以看出为什么必须待定系数的个数应等于单元的自由度数了。 * 由解析几何知，上式中的det是三角形ijm的面积，为使求得的面积值不为负值，节点ijm的次序必须逆时针排列。 * 所有这些转换都是为了将来计算机处理的方便性？ Ni反映单元的位移形态，称为形态函数。 * * * 当选取了单元位移函数，并用节点位移确定待定系数之后，就可以通过几何方程和弹性方程求得用节点位移表示的单元应变和应力。 * 当选取了单元位移函数，并用节点位移确定待定系数之后，就可以通过几何方程和弹性方程求得用节点位移表示的单元应变和应力。 由于其元素均是仅与单元几何性质有关的常量，所以在每一个单元中，各个应变分量为一常数。 * 单元应力只与单元材料性能和几何参数有关，而与单元上点的位置坐标无关，也就是说单元应力再一个单元的内部也是常量。当然，相邻单元一般将具有不同的应力，所以在它们的共用边界上，应力将不连续而有突变。不过随着单元尺寸逐步取小，这种应力突变将会减小，最终收敛于正确解。 * * * 节点力是指连续体离散化后，外载、约束和其它单元通过节点作用在某一单元上的力。对于已从整体结构中取出来的单元来说，作用在其上的节点力就是外力。 这些节点力在单元内部会引起相应的内力。当整体处于平衡状态时，单元在节点力作用下也处于平衡状态。 平面问题：分别用U和V加节点号下标和单元号上标表示该节点水平和垂直分量。 * * 单元节点位移和节点虚位移列阵时常量，即与单元内点的位置坐标无关，可以提到积分号的前面和后面。虚位移是任意的，因此等式两边与它相乘的矩阵应相等，即得。 具有普遍意义的平面问题单元的节点力和节点位移关系式。 * 单元节点位移和节点虚位移列阵时常量，