三维大变形弹塑性有限元分析与边界条件的处理.pdfVIP

第5卷 第4期 塑性工程学报 Vo1．5 No．4 1 998年I 2月 JOURNAI OF PI ASTICITY ENGINEERING Dec． 1 998 三维大变形弹塑性有限元分析与边界条件的处理 TG；7『 ‘ 工业大学 1 30025) 寇淑清 傅沛福 杨慎华 胡 平 李运兴 摘 要 提出适于大应变、大位移与界面摩擦三重非线性的三维弹塑性有限元列式，给出一种处 理三维动态边界摩擦与约束的新方法，并对块体镦挤过程进行了三维数值实验； 关键词 丝塑 塑型!三 ．袭可；，毫段造，踊{ l 前 言 有限元数值模拟技术作为最有前途、应用广泛的工艺仿真工具，在产品优化设计、降低 成本、缩短产品开发周期、预测成型缺陷等方面发挥重要作用 随着精锻工艺的发展，复杂 形状锻件成形过程的三维模拟分析将成为重要的研究方向 ，并且刚塑性有限元在三维体成 形模拟中已取得较大进展。弹塑性有限元考虑了弹性变形和塑性变形的相互关系，残余应力、 残余应变和回弹等问题，对具有载荷大、弹性变形不可忽视等特点的冷锻工艺具有很大的实 用价值 。 。目前，弹塑性有限元法主要用于模具分析及板料成形过程模拟“ 一，在三维体积成 形问题上的应用还很少，原因在于其三重高度非线性性质、传统的界面滑动与摩擦处理造成 的刚度阵非对称性 ，以及求解规模过大等。采用基于T．L法的弹塑性有限元方程处理大变 形问题已有文献报道 ，本文采用逐级更新Lagrange法口 ，建立适宜模拟大变形三维体积成 形的有限元模型，并给出一种处理三维动态边界摩擦与约束的元级代数消元法，为求解大规 模的三维问题奠定基础。 2 大变形三维弹塑性有限元列式 2．1 几何与物理描述 在笛卡尔坐标系下，物体在t。时刻处于状态c。，￡时刻为状态c， + 时刻为状态；。状 态c。时的点Po经位移“．移到状态c的P点，再经时间 移到状态 中的 点，坐标分别为 x ，。， ．，P点的速度为 。于是有如下几何方程： d墨= x 一( +轰 J dx一 一 ． d =(屯+ ． ) l 66 塑性工程学报 第5卷 一 吉(轰+ ,guj+轰 ) 一 吉(轰+轰+瑟 + au． (2) 1( 一 + ) ． ．， 式中 ．——Green应变 ，．— — GIeen应变率 — — 应变率，表示质点处的变形快慢程度 基于J 流动理论的Mises材料本构方程： 一 。f (3) 其中 D嚣 =D (4) 一2G{丢(以 + r ·)一1一r兰f_ ， “) (5) g一一 一。 【 + h】， }一号(E1一击) 式中 、G——分别是泊松比和弹性剪切模量 E一2G (1q-口) E ——单向拉伸真应力与对数应变曲线的切线 一 cauchy应力偏量 ，．— — Kroneker delta符号 一1或o)——分别表示塑性加载和弹性卸载 2．2 单元刚度方程 取t时刻构形为参考构形，虚功率增率原理为： 』 如，