_14-2014中文第14章划分网格的一般原则及计算精度分析概要.pptVIP

网格加密 §14.3 网格加密的h,p,r方法 在应力梯度比较高的地方使用细网格(通常位于不连续处)。 网格密度要平稳过渡，防止相邻网格尺寸相差悬殊。 网格密度要在计算精度和计算量之间折衷。 三种加密网格的方法 1. h-refinement (改变单元尺寸大小) 2. p-refinement (使用插值次数较高的高阶单元) 3. r-refinement (移动节点位置) h-adaptivity p-adaptivity ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? Eight-node elt. ? ? ? ? ? 一般来说，具有不同DOF的不同类型单元不能连接而共享同一自由度。 (例如不能把三维梁单元与平面应力单元连接) 有时候允许某些不同类型的单元连接在同一节点。 §14.4 不同类型单元之间的连接 例：平面Q4与杆单元之间的连接 基本思想：使单元的变形在连接处协调。 x y ? ? ? ? ? 1 4 3 2 5 Bilinear plane element (dofs: ui, vi) Truss element (dofs: ui, vi) Compatible dofs: no problem y ? ? ? ? ? 1 4 3 2 5 Bilinear plane element (dofs: ui, vi) Frame element (dofs: ui, vi, ?i) ? 6 a b L 问题: ? 节点 5 不与节点 2 和 3同坐标 ? 自由度不同: - 1和4节点仅有平动自由度 - 5和6节点有平动和转动自由度 方法: 把刚架节点5上的三个自由度 ( u5, v5 and ?5) 转换为节点 2和 3上的自由度 (u2, v2, u3, v3). 板单元节点上有转角自由度，而实体单元没有。把板单元的转角与实体单元的节点位移协调起来。 u1 u2 q a (u1 - u2 )/a - q = 0 例：三维实体块单元与板单元之间的连接 平衡和协调 节点上的力和力矩平衡条件满足，因为 KU=F 在节点上变形是协调的，这是由刚度阵的组装方式决定的。 在单元之间的边界上，力不平衡，但是随着网格加密，不平衡程度逐渐减小。 单元中心处的应力精度比单元边上的应力精度要高。 在某些高阶单元以及不同类型单元的连接处可能存在单元边界上不协调情况，随着网格加密这种情况的严重程度逐渐减小。 在单元内部，不满足平衡方程。这是由于有限元方程 KU = F 不等价于微元体的平衡方程。 在单元内部满足协调性条件 ，这是由于有限元假定的单元位移函数是连续、单值的。 应力分析结果验证 进行变形、支反力、应力的初步计算。 查看变形形状和主应力向量。 查看等值线的变化。 应力通常依赖于形状和载荷，而与材料特性无关。 位移与材料特性相关。 支反力(力矩)的合力通常与载荷平衡。 对于分布力载荷，计算其合力，检查是否与支反力平衡。 如果软件提供了误差功能，检查计算误差。 §14.5 应力计算结果的性质和处理 §14.5.1 应力近似解的性质 位移元的有限单元法 得到各节点的位移值 得到e=Bu, s=De=DBu 应变、应力的精度比位移低一阶。 B=LN 应力解的近似性体现在： 在单元内部一般不满足平衡方程； 单元与单元之间一般不连续； 在力的边界条件上一般不满足力的边界条件。 §14.5.2 应力精度的改进 1.应力计算公式 应力在单元内分布，对于p 完备的位移模式，应力为p-1次分布 Gauss积分点应力精度最好。 2. 应力磨平： 1 单元应力磨平 (最简单的一种方法，3节点三角形) 可视单元内应力平均值，或形心处应力。 所以，平均应力= 得到的平均应力是矩形大单元形心的值。 2节点应力磨平 3 加权磨平 面积加权 单元面积占单元片面积的权重 4整体磨平 由于位移元得到的应力场不连续，用整体应力磨平改进计算结果，以使应力场全域连续。 (a)构造一个改进后的应力解， ne是单元节点数， 是第i点的插值函数 是改进后的节点应力， 是改进后的应力场函数。 由于 变分得到 为由节点位移计算得到的应力场。 注意到 于是有： 这是一个未知数为 线性代数方程组，求解后可以得到改进后的节点应力 ，从而由插值得到全域内连续的应力场。 得到各个节点改进的应力值，单元交界面上是连续。 显然整体磨平工作量很大，相当于又进行一次有限元计算。 5 利用应力泛函的单元内磨平： 应力在单元内磨平，减少工作量。 原全域应力场 现单元应力场 单元足够小时，令C=I，有 对于等参元