基于阶级阶级理论的钢板渐进折弯三维数值模拟.docxVIP

基于阶级阶级理论的钢板渐进折弯三维数值模拟 大、金属大板件的零件只能通过单通道弯曲法采用多功复杂的加工和制造工艺。多通道弯是一种柔性的板料形成工艺。计算机数字控制每个插入的板的进给量，worm模式下的位移参数被抑制，以渐进形成高精度零件。渐弯形成是一种非常复杂的物理过程。采用传统方法难以分析和解决板材料的渐进式弯曲变形和回弹问题。作者采用了比例素描分析法建立了板材料回弹数学模型，但该模型仅适用于单一通道的折叠，并且有限面积内的形状与回弹非线性相交的多通道弯曲有其自身的局限性。 有限元方法已成为预测板料成形回弹的重要手段.文献分别基于显式和隐式的有限元算法模拟了板材成形回弹过程,但预测精度不是很理想.其主要原因是目前在板料折弯成形领域,还是沿袭传统的描述各向同性的Mises屈服准则,视各向同性材料模型为材料本构模型.这使得对各向异性板材成形的有限元预测成形和回弹过程中的应力场精度不会很高.由于本课题研究的大幅面钢板具有较明显的平面各向异性和厚向各向异性,研究金属的变形只能用各向异性塑性理论来描述弹塑性变形体中应力(应力率)与应变(应变率)之间的本构关系.Hill和Barlat在这一领域的研究成果最具代表性.目前,有限元求解中一般采用描述厚向异性的Hill48屈服准则和用于平面应力的Barlat89各向异性屈服准则.Hill48各向异性屈服准则在描述ˉγγˉ(各向异性系数)值较大(ˉγ1)(γˉ1)的金属钢板的各向异性行为时比较准确,而Barlat各向异性屈服准则能精确描述ˉγγˉ较小(ˉγ1)的铝合金板料的屈服行为.遵循上述规则,本文设想在Hill各向异性屈服准则、非线性随动强化模型、平面应变假设条件下,构建金属钢板折弯时反映其真实应力应变关系的本构方程,作为大幅面钢板渐进折弯成形制造的理论基础.最后通过编程并数值模拟大尺度U形板材工件的多道次渐进折弯成形过程. 1 材料各向异性初始屈服准则 Hill提出用一个二次屈服函数来描述各向异性材料的塑性行为.为简单起见,只考虑每一点上具有3个互相垂直的对称平面的各向异性,其交线叫各向异性体的主轴.其屈服函数为 f(σij,εij)=F(σy-σz)2+G(σz-σx)2+Η(σx-σy)2+2Lτ2yz+2Μτ2zx+2Ντ2xy=2ˉσ2e.(1) 式中:F、G、H、L、M、N为各向异性瞬态特征参数;σx、σy、σz分别为正交对称应力轴(x、y、z)上的单向拉伸屈服应力;τyz、τzx、τxy分别为yz-、zx-、xy-坐标平面上的剪切屈服应力;ˉσe是等效应力. 各正交对称应力轴方向上应变增量关系表示为 dεx=dλ[Η(σx-σy)+G(σx-σz)],dγyz=dλLτyz,dεy=dλ[F(σy-σz)+Η(σy-σx)],dγzx=dλΜτzx,dεz=dλ[G(σz-σx)+F(σz-σy)],dγxy=dλΝτxy. 对于各向异性瞬态特征参数,可由下式确定: 从式(2)可以看出,要完整描述一个单元体中的各向异性状态,就需要知道各主轴的方位以及6个互相独立的屈服应力值.迄今还不能把屈服应力与微观结构(如择优方位的程度)联系起来,因此只能通过实验来确定材料各向异性系数.经拉伸试验可得到材料各向异性系数(ˉγ)?r0、r45、r90分别为相对于冷轧方向 0°、45°、90°的各向异性参数,其关系式为 ˉγ=r0+2r45+r904. 在板材折弯成形时,平面应变条件下的Hill各向异性屈服准则中特征参数可通过下式求解: F=r90τyz2(r0+r90)(2r45+1);G=r0r90τzx2(r0+r90)(2r45+1);Η=r0τxy2(r0+r90)(2r45+1).} 根据式(1),省去高价项和忽略剪切效应的影响,Hill各向异性理论的屈服条件可表述为 f(σij,εij)=σ2x+σ2y-2ˉγ1+ˉγσxσy+21+ˉγ(σ2z-σxσz-σyσz)-ˉσe2=0.(3) 对于板料折弯成形过程,如采用平面应变假设,即εx=0,使得dεx=dλ?f(σij,εij)?σy=0(dλ是一系数).与式(3)组合求解得到ˉσe,其表示为 ˉσe=√1+2ˉγ|σx-σy|, (4) ˉσe=√1+2ˉγ1+ˉγ|σy-σz|. (5) 根据塑性功公式,等效应变可表示为 ˉεe=1+ˉγ√1+2ˉγ|εy|. (6) 根据非线性随动硬化模型,等效应力与等效应变的关系表示为 ˉσe=Κ(ˉεo+ˉεe)n. (7) 式中,ˉε0为初始应变,ˉεe为等效应变. 假定在板料长度方向(轧制方向)的应变近似为零、折弯过程中板材体积保持不变及折弯工件的应力中性层与几何中性层相重合,即折弯工件的中性层曲率半径rn可表示为rn=r+t/2(见图1)