各向异性屈服准则对铝合金成形过程预测的影响.pdfVIP

摘要 采用完全隐式向后 Euler 回映算法，通过 ABAQUS 软件中的用户材料子程序 UMAT ，将Hill 1948 ，Barlat YLD91 ，Barlat YLD2004-18p ，CMTP 各向异性屈服 准则嵌入到商业有限元软件中。以圆杯拉深为例，采用 Hill 1948 ，Barlat YLD91 和 CMTP，对深冲过程中产生的起皱进行了预测；基于Hill 1948 ，Barlat YLD91 ， Barlat YLD2004-18p 和 CMTP 各向异性屈服准则，预测了深冲过程中的制耳行为， 并和实测值、解析方法进行了比较。分析了深冲过程中等效应力、等效应变的分 布情况，此外考虑了工模具之间的摩擦和压边力对深冲制耳的影响，主要工作和 结论如下： （1）本文推导了 Mises 、Hill 1948 、Barlat YLD91 、Barlat YLD2004-18p 和 CMTP 屈服函数的“三参数”，其皆可实现于ABAQUS 软件中。 2 f (σ) σ UMAT （）对任意形如 的屈服函数，编写了通用的 用户子程序接 口，该接口具有普适、精度高、收敛速度快的特点。采用此接口，将 Hill 1948 ， Barlat YLD91 Barlat YLD2004-18p CMTP ABAQUS ， ， 屈服函数在 软件中成功地实 现。 （3 ）利用Hill 1948 ，Barlat YLD91 ，CMTP 屈服函数对 3104 铝合金板深冲过 程进行模拟。模拟结果表明，Hill 1948 ，Barlat YLD91 和 CMTP 都可以预测到深 冲过程中的起皱现象，但采用不同的屈服准则模拟的趋势有所不同，采用屈服准 则 Hill 1948 和 Barlat YLD91 预测的起皱趋势稍差，而 CMTP 屈服准则预测结果具 有较高的起皱趋势。 （4 ）采用 CMTP 屈服函数可将 2090 铝合金板中的织构信息和由其产生的制 耳行为直接联系起来；而 Barlat YLD2004-18p 则输入更多的材料实验数据，如沿 轧制方向每隔 15°的屈服应力、r 值和平衡双轴屈服应力、双轴 r 值，以此来捕获 2090 铝合金的各向异性从而预测制耳。模拟结果表明：深冲过程中，法兰区周向 受压而增厚，凹模圆角处变薄，杯壁减薄，杯底厚度基本不变，等效应力、等效 应变最大值首先出现在凹模圆角处。深冲完成后，杯壁径向由上而下厚度逐渐的 变薄，等效应力、等效应变最大值出现在法兰区的外缘。 （5 ）Barlat YLD2004-18p 屈服函数模拟的制耳轮廓和实测值比较吻合，并且 Barlat YLD2004-18p 屈服函数模拟的精度要比 CMTP，Barlat YLD91 和 Hill 1948 屈服函数要高，解析方法通过 r 值（塑性应变比）可以预测到制耳的个数，但是预 测制耳轮廓的精度不高，只可以作为一次近似。 （6 ）随圆片和凹模、凸模和压边圈之间摩擦系数及压边圈所施加的压边力的 - I - 增加，圆片深冲变形过程中径向受到的拉应力相应增加，从而导致制耳轮廓高度 也增加。 关键词：各向异性，深冲，有限元模拟，制耳， 起皱 - II - Abstract The Hill48, Barlat YLD91, Barlat YLD2004-18P, CMTP theories and the finite element method are combined using a full implicit backward