有限元方法预案.pptVIP

北京机械工业学院 * 步骤： 代入：已知外力F1x， F1y， F2x； 支反力F4x， F4y， F5y， F6y以外诸力为零； 已知位移约束条件u4，v4，v5，v6＝0 消元：消去对应于零位移（ u4，v4，v5，v5 ）的行和列，同时消去位移列 阵和外载列阵{F}中相应的元素。 （26） ＝0 ＝0 ＝0 ＝0 0= 0= 0= 0= 0= 0= u1 v1 v2 u3 v3 u4 v4 u5 v5 u6 v6 u2 北京机械工业学院 * （27） 则方程成为： 则方程（27）成为可解。 先解位移ui，vi，再计算应变、应力。 （3）整体刚度矩阵为稀疏矩阵 北京机械工业学院 * 因每结点的刚度项只与相邻结点有关，故两不相邻的结点间kij为零。 结论： 划分结点时，要尽可能使同一单元的几个结点号码尽可能接近。 北京机械工业学院 * 五、有限元建模 有限元模型包括：单元、结点、边界条件、外载荷、材料特性 模型质量的影响：分析精度、效率、计算机CPU时间和存储量 （1）常用单元类型及网格划分 一维 二维 二维 三维 网格划分：应力变化梯度或应力集中处单元划分密。 北京机械工业学院 * 结点编号： 同一单元的几个结点号码尽可能接近。 北京机械工业学院 * （2）结点外载荷：据力的等效原理，将外载荷简化到结点上。 （3）边界约束条件：平面问题至少三个约束 1、将支承转化为约束 2、构件和载荷为对称时，可取对称轴线的一侧计算，并给出相应约束 北京机械工业学院 * 六、有限元法解题步骤 北京机械工业学院 * 七、有限元软件 第一代：人工完成前、后处理，如结点号、单元号及结点 坐标的计算及逐点输入，例SAP 半自动：半自动生成网格，给出关键点及区域，例Super SAP 全自动：自动生成网格 各种图形显示，例Ansys 北京机械工业学院 * 有限元解的性质和收敛性 1、收敛准则的物理意义（以平面问题为例） 其中 ，常数项代表单元 的刚体位移模式 是u、v的一次导数，用m＝1表示导数次数 . . . . . . 位移 插值函数 应变 常数项代表单元的常应变状态 位移插值函数必须至少是x、y的一次完全多项式 北京机械工业学院 * 安全性要求（完备性）： 由插值函数构成的有限元解必须能反映单元的刚体位移和常应变状态，否则单元内部将产生非零或非常值的应变，则有限元解不可能收敛于真正解。 如场函数的最高阶导数是m阶,则单元插值函数至少是m次完全多项式(即高至m阶导数都有常数项)。 连续性要求（协调性）： 位移函数自身在单元交界面上是连续的，否则，交界面上引起无限大应变，有限元不能收敛于真正解。 如场函数的最高阶导数是m阶，则单元插值函数在交界面上必须具有（m-1）阶的连续导数。 3结点三角形的插值函数满足以上二要求。 2、有限元解的收敛准则 如采用完全多项式作为单元的插值函数，则有限元解在一个有限尺寸的单元内可以精确地和真正解一致，但实际上只取有限项多项式，故结果是近似的。 北京机械工业学院 * 设取p次多项式为插值函数，m阶导数是非零的，则须 ，有 常数项 可知，p次完全多项式 的高至m阶导数的表达式中都有常数项，当单元尺寸趋于零时，每单元内 及其至m阶导数将趋于它的精确值，即趋于常数。如插值函数还满足连续性要求，则有限元阶趋于精确解，即收敛。 北京机械工业学院 * 准则1：完备性要求 如果场函数的最高阶导数是m阶，则单元插值函数至少是m次完全多项式（即函数本身和高至m阶导数都有常数项）——收敛条件之一。 如应变、应力（磁、温度场） 准则2：协调性要求 如果场函数的最高阶导数是m阶，则单元插值函数在单元交界面上必须具有函数直至（m－1）阶的连续导数——收敛条件之二。 系数齐全 有限元解收敛条件： 选取的单元既完备又协调。即当单元尺寸趋于零时，有限元解趋于真正解。 （如：梁、板、壳的变形能含有位移二阶导数，故位移至少是二次函数，梁设为三次函数） 北京机械工业学院 * 3、收敛速度和精度估计 单元尺寸为h，则x、y与h同量级，则位移u的误差是 如一维 误差估计： 如平面问题中位移u是p阶完全多项式： 泰勒级数： 如平面问题采用线性插值函数： p＝1，则u的误差是 ，应变误差是 应变由位移的m阶导数得到，则其误差是 误