有限元法基础-6应用考虑答案.pptVIP

* * * * * * * * * * * * * * * * * * * * * 6. 5 子结构 有限元法基础 * 为解决计算机求解能力的问题，尤其是对大型结构分析时，提出了总体-局部方法、子结构、超级单元法，其目的是提高分析效率，并能够在硬件条件有限的情况下，能够分析大型结构问题。 6. 5 子结构 有限元法基础 * 总体-局部方法主要用于应力集中问题 先由较粗的网格进行整体分析，然后在关心的局部加密网格将总体分析的位移解施加在局部模型的边界。 6. 5 子结构 有限元法基础 * 构建超级单元的目的是减少“计算时间”上的重复重复性，超级单元实际上也是一种子结构。 6. 5 子结构 有限元法基础 * 对于几何空间上有重复的大型结构，构建子结构减少重复计算量。 6. 5 子结构 有限元法基础 * 实施步骤 1）取具有重复性的结构作为子结构（可以有多级子结构） 2）对底层子结构形成刚度矩阵，将其内部DOF凝聚 边界DOF 内部DOF 6. 5 子结构 有限元法基础 * 凝聚内部DOF 由第二组方程 6. 5 子结构 有限元法基础 * 6. 5 子结构 有限元法基础 * 3）将子结构装配，形成上一级子结构（采用坐标变换并再凝聚） 4）组装子结构，并形成整体刚度矩阵，求解 5）将结果进行回代，再求解子结构内部的节点位移和其他量 6. 5 子结构 有限元法基础 * 采用Gauss消去法进行内部DOF的凝聚 将矩阵消元 第一方程组为子结构的凝聚后的矩阵 第二方程组为子结构回代求解方程 6. 6 自适应分析 有限元法基础 * 如何对有限元分析结果的误差进行判断，其划分的网格是否合适？ 自适应有限元分析是基于前一次计算结果，进行事后误差估计，并以提高精度为目的，由分析软件自动对计算方案重新调整，重新计算，直至满足精度的计算分析方法。 6. 6 自适应分析 有限元法基础 * 调整方式 1） h方法（h-refinement method）：逐步细化网格，减小单元尺度 2） p方法（p-refinement method）：在不改变有限元网格剖分的基础上，提高单元插值函数阶次 3） r方法 （r-refinement method）：不改变单元类型和节点数目，调整网格节点布局 4） h-p混合方法 6. 6 自适应分析 有限元法基础 * 几种细化方法 6. 6 自适应分析 有限元法基础 * 自适应分析涉及的关键技术 1）误差估计 误差估计是其理论基础 2）单元升阶 提供足够的有效的高阶单元 3）网格再生成 自动网格剖分的计算方法 6. 6 自适应分析 有限元法基础 * 误差估计（Zienkiewicz和Zhu于1987提出） 定义误差 定义应变能范数 定义应变能误差范数 经过平滑后的 单元中直接计算 6. 6 自适应分析 有限元法基础 * 相对误差 在每个单元里计算，可接受的值 实际计算中可简化为 6. 6 自适应分析 有限元法基础 * 自适应网格划分 目标为所由单元的 不超过允许值 。在误差估计区有m个单元，其允许的误差 在一个典型单元i中 6. 6 自适应分析 有限元法基础 * 当 时，表明单元比需要的小； 当 时，表明则个区域需要细化网格。 细化网格的希望的尺度为 p为单元场插值的完全多项式的最高阶次， 为奇异点的强度。 6. 6 自适应分析 有限元法基础 * 6. 6 自适应分析 有限元法基础 * 自适应分析的过程 * * * * * * * * * * * * * 6. 2 应力计算结果的处理 有限元法基础 * 节点应力平均 节点i 有m个单元与之相关 6. 2 应力计算结果的处理 有限元法基础 * 整体应力修匀 应力常在全域是不连续的，可采用整体应力修匀的方法，以得到应力场全域连续。 构造一个改进的应力解 Ni 是插值函数， 改进后的节点应力值，ne是单元节点数。 6. 2 应力计算结果的处理 有限元法基础 * 改进的应力与有限元解 应满足加权最小二乘的原则 变分取驻值 代入 表达式，由 的任意性 N为全部单元的节点数。 6. 2 应力计算结果的处理 有限元法基础 * 可得NxS个线性方程组