第二讲-建模.pptVIP

第二章、结构静力分析的有限单元法 第一节、单元特性的导出方法 进行有限元分析的基本步骤之一是：找出所剖分的单元的刚度矩阵，建立刚度矩阵的基本方法有： 直接方法 虚功原理法 能量变分原理方法 迦辽金方法 虚功原理方法介绍 以平面三角形单元为例： 已知量：三个节点的坐标；未知量：三个节点的位移 1、设定位移函数（形函数、插值函数） 写成矩阵： 单元内“任意点”的位移通过“节点位移”插值得到。对于三角形单元，可假定单元内的位移为 的线性函数。 (1) (2) 能否为二次 或三次函数？ 既然是单元内任意一点的位移表达式，当然也是三个节点上的位移。 可以得到： 这里: 是广义坐标。 (3) (4) C 是与节点坐标相关的矩阵 节点位移 式中： 为不使 A 为负值，图1中的i，j，k必须按逆时针方向标注 任意点位移 任意点坐标 节点几何量 节点位移 将 代入 有： 或写成： 可简写为： 上式建立了单元内“任意点的位移”与“节点位移”的插值关系，称N为形状函数或插值函数，其中： 选取单元位移函数的一般原则： 广义坐标 是由节点场变量确定，因此，其个数必须与节点自由度个数相等； 选取多项式时，常数项（反映刚体位移）和坐标的一次项（反映常应变）必须完备；当单元数目趋于无穷时，单元缩小趋于一点，此时单元应变为常应变； 多项式的选取应由低阶到高阶，尽量选取完全多项式以提高单元精度。 2、由位移函数求应变 由弹性力学的知识： 可得： 或写成： 节点位移 几何矩阵 节点位移 任意点应变 3、由应变求应力 对于平面问题： 对于平面应力： 4、由虚功原理求单元刚度矩阵 给单元节点以任意虚位移： 虚位移、虚应变分别可表示为： 虚功原理：当结构受载荷作用处于平衡状态时， 在任意给定的虚位移下，外力 及内力 所 做的虚功之和为零。 单元节点的外力虚功： 单元内力的虚功： 即 由虚功方程： 是任意的，可消去： 或： 其中 把 B 和 D 代入上式，得平面应力问题的三角形单元刚度矩阵： 是单元厚度 单元刚度矩阵任意元素 的物理意义 当单元的第j (j=1,2,…6)个节点位移为单位位移而其它节点位移为零时，须在单元第 i (i=1,2,3,…6)个节点位移方向上施加的节点力大小。 显然，单元刚性越大，则使节点产生单位位移所需施加的节点力就越大。 单元刚度矩阵的特点 对称性 奇异性 物理解释：单元处于平衡时，节点力相互不是独立的，它们必须满足平衡方程，因此，它们是线性相关的。另一方面，即使给定满足平衡的节点力 f，也不能确定单元节点位移 ，因为单元还有任意的刚体位移。 主元恒正 物理解释：要使节点位移 ，施加在 方向的节点力 必须与位移 同向。 第二节、有限单元法的解题步骤 1、单元剖分和插值函数的确定 根据构件的几何特性、载荷情况及所要求的变形，建立由各种单元组成的计算模型。然后再按单元的性质和精度要求，写出表示单元内任意点的位移函数。 利用节点，可以写出以 表示的节点位移 求 及 ，并代入 它是用节点位移表示单元内任意点位移的插值函数式。 2、单元特性分析 根据位移插值函数，由弹性力学给出的应力应变关系，计算应变为： B 是应变矩阵。相应的变分为： 由物理关系（本构关系） 由虚功原理： 得节点力与节点位移间的关系为： 其中： 3、单元组集 把各单元按“节点”组集成与原结构相似的整体结构，得到整体结构的节点力与节点位移的关系，即整体平衡方程： 对于结构静力分析，载荷列阵 f 可包括： 其中：体积力 表面力： 集中力： 4、解有限元方程 解题前，需对结构平衡方程组进行边界条件处理；（直接代入法；置大数法；对角元素置1 法） 采用不同的方法解有限元方程，得出节点位移。 直接代入法 为待定位移，　为已知位移。 该方法要重新组合方程，组成的方程组阶数降低了，但结点位移的顺序被打破，给编程带来困难。 对角元素置 1 法 对于给定位移值是零位移（铰支座）时，由于　　　　，可以采用对角元素置1 法。　　　　 置大数法 当有结点位移为