有限单元法应用中的若干实际考虑要点.ppt

单元的自由度凝聚： （5.5.7） 由上式中的第二式可解出： （5.5.9） 将上式的第一式，有 整理，消去?e , 有 令： （5.5.10） 说明： （1）方程（5.5.10）包括了附加内结点位移项得到的单元刚度阵和载荷列阵。 （2）方程（5.5.10）阶数与原线性单元相同。消去附加自由度 ?1~?4 的过程，称为自由度凝聚。 （5.5.10） 说明： （1）方程（5.5.10）包括了附加内结点位移项得到的单元刚度阵和载荷列阵。 （2）方程（5.5.10）阶数与原线性单元相同。消去附加自由度 ?1~?4 的过程，称为自由度凝聚。 （3）若不存在体积力（ f ≡ 0）, 则有 进一步略去 中的第二项，则 —— 与原线性单元中相同 实践证明， 作以上处理后，计算量大大减少，且对精度影响不太大。 例： 4050 3000 101.3 100.0 非协调元网格2 4050 3000 101.5 100.0 非协调元网格1 2954 2188 72.3 70.6 协调元网格2 2945 2182 70.1 68.1 协调元网格1 4050 3000 103.0 100.0 理论解 载荷B 载荷A 载荷B 载荷A i 点的弯曲应力 j 结点的位移 悬臂梁受载荷A和载荷B作用，如图所示。采用4结点矩形单元计算。 其中： —— 附加的内部自由度 1 2 3 4 存在的问题： （5.5.4） 单元边界上位移分布： —— Wilson元 不满足协调条件，为非协调单元。 相邻单元边界的位移不能连续。 Wilson非协调元单元收敛性如何？ 实践证明： 对于 C0 类型问题，若在单元尺寸趋于零（即单元应变趋于常应变）时，其位移的连续性能得到恢复，则非协调元的解仍能趋于精确解。 检验非协调元是否收敛性的条件为： （1）位移模式能否描述常应变？ （2）在常应变的条件下，能否自动地保证位移的连续性？ 检验方法： 分片试验 —— 检验采用非协调元任意网格单元时能否达到上述连续性的要求。 —— 能通过分片试验的非协调元，其有限元解一定收敛于精确解。 * 有限单元法应用中的若干实际考虑 第 5 章 要点： （1）结点应力的计算与修正； （2）结构特点的考虑； （3）非协调单元简介 主 要 内 容 5.1 引 言 5.2 应力计算结果的性质与处理 5.3 子结构法 5.4 结构对称性与周期性的利用 5.6 小 结 5.5 非协调元与分片试验 5. 1 引 言 1. 有限单元法的求解过程 （1）划分单元，输入结点和单元信息； —— 前处理器 （2）单元分析：N、Ke、Pe； （3）整体分析： 引入位移边界条件，得到 （4）求解方程 —— 得解a ； （5）对位移 a 结果进行有关整理、计算单元或结点的应力、应变。 —— 求解器 —— 后处理器 的可视化表示。 2. 目前存在的问题： （1） 的精度较低。 如何由应力、应变结果的特点改善其精度？ （2） 如何利用结构的几何特点、受力特点简化计算，减少工作量，提高计算效率？ （3） 如：结构与受力的对称性、周期性物质等； 子结构法； 非协调元概念与应用 （Wilson非协调元）。 5.2 应力计算结果的性质与处理 应力、应变的计算： —— 精度较低。 误差的原因： （1） 单元内平衡方程不能精确满足； （2） 单元交界面上应力不连续； （3） 边界上边界条件不能得到精确满足等； 5.2.1 应力近似解的性质 1. 位移解 a?的性质 a?——有限元近似解， a ——真实解； 由最小位能原理，可知，a?具有下限的性质： 原因： 单元离散等相当于加大了原结构的刚度。 2. 应力、应变解 ??、 ??的性质 设u? 、? ?、 ??——近似解， u、? 、 ? ——真实解，有 近似解对应的位能： ——?P(u)实际的总位能 —— ??P(u)=0 —— ?2?P(u) 2. 应力、应变解 ??、 ??的性质 在线弹性下，有 对于一具体问题， ?P(u)应为一定值， 则 ?P(u*)的极值问题归结为： ?2?P(u)的极小值问题。 将?2?P(u)表示成单元位能泛函的形式，有 （5.2.5） （5.2.6） （5.2.5） （5.2.6） 上式表明： ?2?P(u)的极小值问题 求解 的加权 二乘最小值问题。即 ??、 ??为?、 ? 在加权（D、C）最小二乘意义下的近似解。 ??、 ??的特点： （1）??、 ?? 在真正解?、 ? 上下振荡； （2）在某些点上有：?? =? 、 ?? =?，即存在最佳应力点。 利用??、 ?? 的上述特点，作适当处理，可提高应力、应变结果的精度。 5.2.2 等