刚塑性有限元分析时计算节点应变值的新方法.pdfVIP

维普资讯 第 18卷 第 3期 太 原 重 型 机 械 学 院 学 报 Vo1．18 NQ3 1997年9月 JOURNALOFTAIYUANHEAVYMACHINERYINSTITUTE Sep．1997 D一己 刚塑性有限元分析时计算节点 ，。^L／ ／ 应变值的新方法 3、 郝南海 刘 岩 (太原重型 军j_力加工系，太原030024) Csg．t岛岛港港务局五公司，秦皇酉066002) 摘 法不仅 美 由 利用刚塑性有限元法模拟金属成形过程时，随着变形的进行，初始规则的有限元网格逐 渐畸变，引起计算上的困难和精度上的降低。当有限元网格畸变到某种程度时，必须停止计 算，重新划分适合于计算的网格。重新划分网格包括两方面的工作，一是在已变形的工件上 确定出规则的新网格，二是把已畸变的旧网格上的信息转换到新网格上。 在有限元计算中，插值函数在插值区域内部精度较好，而在区域边界上，精度较差，计算 时一般取单元形心处的应变值作为单元的平均应变值。新旧网格信息转换时，先根据旧网 格单元的应变值计算出节点的应变值，再由节点应变值通过插值确定新网格节点应变值。 s．I．Oh等人总结了计算节点应变值的三种方法，面积加权平均法，最dx--乘拟合法和 混合法_lJ。本文在分析现有方法优劣的基础上，提出了一种新的计算方法，边界修正法。该 方法不仅容易实现，而且能保证一定的计算精度。 1 面积加权平均法 设某节点 处于 个单元的边界上，则以相关单元的面积作为权数，求其平均值。计算 公式为 ￡= ∑A南／ (1) d 式中 ￡—— 节点 的应变值；e—— 节点 周围单元s的应变值；A—— 单元 s的面积 收稿 日期：1997—03—01 本文第一作者：男。1963年2月生，压力加工景讲师，暗工大在读博士。研究方向为塑性加工 摸报与控嗣。 维普资讯 第 18卷第3期 郝南海等 剐塑性有限元分析时计算节点应变值的新方法 261 ∑—— 对节点 周围所有单元求和 对于内部节点，面积加权平均法能得到较好的精度，对于边界节点，一般周围只有两个 单元，有时甚至只有一个单元．计算误差较大。 2 最小二乘拟合法 构造一个拟合面，使该拟合面与单元应变的给定值之间的均方差为最小，并迫使拟合 面在单元形心处通过给定值，再由拟合面确定出节点应变值，即通过使泛函，最小化来求解 节点的应变值I：’I∑(∑ f—C=)『dAj+KI∑(∑El／A—G)dAj (2) i i ’ i i 式中G——单元J应变给定值；gi——待求的节点i的应变值；K——大的正数，迫使拟合 一 r 面通过给定值；∑——对每一个单元都求和；ld ——在整个区域上求和 』 最dx--乘拟合法所得结果在工件各处都有较好的精度，但计算工作量较大，需求解线性 方程组。 3 混合法 s．I．Oh等人建议，为了提高计算效率并保持一定的精度，可对内部节点采用面积加权 平均法，对边界节点采用最dx--乘拟合法。笔者经实际计算发现，采用混合法计算时，由于 边界单元所涉及到的内部节点其应变值由面积加权平均法首先给出，为在边界单元内使泛 函最小，有时边界节点的应变值计算精度会低于单纯采用