材料成形回弹分析.pdfVIP

中国机械工程第 13 卷第 17 期 2002 年 9 月上半月 文章编号: 1004- 132 (2002) 17- 1458- 04 板材成形回弹数值分析的静力隐式方法 蔡中义　李明哲　李湘吉 　　　　摘要: 给出了板材成形加载过程、卸载过程数值模拟的静力隐式有限元 方法, 通过平衡迭代求解有限元方程, 并采用隐式算法进行应力积分及接触 压力计算。这种方法应力计算准确, 因而回弹计算结果具有较高的精度。数 值计算表明, 纯弹性板大变形卸载后能够完全恢复到原始形状, 说明该方法 对几何非线性问题的分析是非常准确的。最后对板材柱面、球面成形的回弹 进行了分析。 蔡中义　副教授 关键词: 板材成形; 回弹; 静力隐式; 弹塑性有限元 中图分类号: T G386; TB 115　　　文献标识码: A 　　回弹是板材成形中不可避免的现象, 是决定 1　静力隐式有限元理论 零件最终形状的重要因素。如果回弹量不能有效 地控制, 将严重影响零件的几何精度, 成为成形缺 1. 1　有限元方程 陷。回弹计算是板材成形数值模拟的一个十分重 板材成形是一种拟静态的变形过程, 这一变 要课题[ 1, 2 ]。传统的回弹数值模拟通常采用动力显 形过程由各物质点在空间的运动构成, 各物质点 [ 3 ] 式有限元方法 , 由于显式时间积分格式不形成 的运动可由位移来描述。在任意时刻, 物质点都处 总刚度矩阵, 不进行力的平衡迭代, 应力计算不准 于平衡状态。采用有限元格式空间离散, 并考虑到 确, 因而无法给出准确的回弹计算结果。 工具与板材间的接触作用后, 平衡方程可简单表 回弹计算是在加载过程应力分析的基础上进 示为 行的。为精确计算加载应力必须采用静力隐式有 R in t + R ci = R ex t + R ce ( 1) [ 4 ] 式中, R in t 为节点的内力向量; R ex t 为作用在节点的外力向 限元格式 。这种分析格式采用迭代方法求解非 线性有限元方程, 迭代过程中不断更新切线刚度, 量; R ci 及R ce 分别为由接触作用产生的节点内力向量及外 直至节点不平衡力及位移增量小于给定值为止。 力向量。这些向量都是由各单元对节点的作用累加后得 到的。 由于每一增量步都进行平衡迭代, 因而加载过程 在板材成形分析中, 最令人感兴趣的是变形。 的应力计算精度高, 以此为基础进行的卸载回弹 计算也比较准确。 因此, 需要从上述方程中求解出位移U 。然而, 由 于R in t 及R ci 是位移的复杂非线性函数, 方程不能 收稿日期: 2001—01—16 直接求解, 必须将外载荷划分成一系列的载荷增 ( ) 量, 将非线性方程组线性化。 基金项目: “十五”国家科技攻关课题2001 203 11 ; 教育部跨 BA B 世纪优秀人才专项基金项目 从 t 时刻给定的平衡状态出发, 施加载荷增 量后产生变形, t + t 时刻新构形下的平衡方