铜管游动芯头拉拔过程应力应变场的数值模拟.pdfVIP

铜管游动芯头拉拔过程应力、应变场的数值模拟 林鹏 张世英 辽宁工程技术大学材料科学与工程系，辽宁阜新（123000） E-mail ：songbo1984923@163.com 摘 要：采用国际上通用的非线性有限元软件 Marc 对铜管游动芯头拉拔过程进行了数值模 拟和分析研究。通过模拟分析，清楚的再现了游动芯头拉拔过程中铜管的塑性变形过程，得 出了应力、应变分布并进行了分析，同时合理的揭示了拉拔过程中铜管缺陷产生的原因，对 提高铜管质量有指导意义。 关键词： 铜管 游动芯头拉拔 数值模拟 marc 0．引言 随着空调、冰箱市场的逐年扩大，对铜管的需求增多的同时对铜管的质量要求也越来越 高。由于现场条件与观测条件的限制，人们很难实时、定量地把握金属成形过程，也就很难 找出在拉拔过程中产生缺陷的真正原因，为能有效的解决拉拔过程中出现的问题，对游动芯 头拉拔过程进行数值模拟分析研究具有非常重要的意义[1-3] 。 本文采用国际上通用的非线性有限元软件 Marc 对铜管游动芯头拉拔过程进行了数值模 拟和分析研究。通过模拟分析，清楚的再现了游动芯头拉拔过程中铜管的塑性变形过程，得 出了应力、应变分布并进行了分析，同时合理的揭示了拉拔过程中铜管缺陷产生的原因，对 提高铜管质量有指导意义。 1．有限元模型的建立 根据其轴对称性，可将铜管模型简化为二维轴对称模型，模型都用四边形表示。铜管单 元体可取长 15mm ，宽等于实际加工的铜管初始壁厚 1.65mm 的平面单元，将其划分 720 个 四边形四节点单元，总节点数 819。铜管材料为磷脱氧铜，牌号为TP2，弹性模量为 1.172 ×1011N/m ,泊松比为 0.33 。模具的模型建立通常是把它建成刚性接触体 （不变形接触体）。 本文为了也能够体现模具的应力、应变情况，所以将模具也设置为单元体 （可变形接触体）， 模具的外形并不十分规整，在 Marc 中用AUTOMESH （通用的自动划分面网格和体网格的 工具）对其划分为外模，162 个四节点单元，总节点数 207;芯头，132 个四节点单元，总节 点数 172 。模具的材料为硬质合金 （YG8 ）,弹性模量为 7 ×1011N/mm2 ，泊松比为 0.21 。整 体模型如图 1 图 1 整体模型 材料的物理性能 （热传导系数、热容及热膨胀系数））是随温度变化的，变形抗力是热 -1- 力学参数 （变形程度，变形速度和变形温度）的函数，这些数据均可在MARC 的材料库中读 取[4-5] 。为了能够得出更加准确的模拟结果，本文将采用Marc软件的热—机耦合分析方法， 即温度场的求解和应力场的求解合在一起的大变形热弹塑性有限元来分析。对于接触摩擦问 题的处理，在接触计算中采用了直接约束法,为能保证精度采用了修正的剪切摩擦模型： σ 2 V ⎛ ⎞ σ ≤−m arctan r •t f ⎜ ⎟ 3 π rv ⎝ ⎠ rv 表示发生相对滑动时接触体之间的临界相对速度 σ 等效应力 σf 切向摩擦应力 V 管与模具的相对滑动速度向量 r m 摩擦因子 t 相对滑动速度方向上的切向单位向量 2．模拟结果及分析 通过数值模拟可以使铜管在拔制过程中应力场分布可视