发明名称 一种超高强度钢的热成形瞬态成形极限的模拟预测方法
摘要 本发明公开了一种超高强度钢的热成形瞬态成形极限的模拟预测方法,包括如下步骤:确定稳态平衡条件下材料本构模型;在稳态条件下基于材质不匀沟槽假设的集中性失稳理论确定热成形极限模型,其中稳态条件下的热成形极限模型中将采用确定的稳态平衡条件下的板料本构模型;运用所述热成形极限模型进行瞬态工艺条件下热成形极限的模拟预测。
申请公布号 CN102708303B 申请公布日期 2015.11.11
申请号 CN201210192708.2 申请日期 2012.06.12
申请人 湖南大学 发明人 李光耀;李洪周;孙光永;陈响军
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 北京天达知识产权代理事务所(普通合伙) 11386 代理人 王庆海
主权项 一种热成形工艺瞬态成形极限的模拟预测方法,包括如下步骤:步骤1:确定稳态平衡条件下材料本构模型,包括在稳态条件下进行热成形板料高温单轴拉伸实验,获得热成形材料各向异性和流动应力的参数;定义相关屈服准则和硬化模型,通过上述参数回归计算材料本构模型的参数,参数拟合是以本构模型预测的结果与实验数据误差最小为目标,本质上是一个最小化优化问题,所述屈服准则为BBC2005准则,该准则为高阶各向异性屈服函数,其中对处于高温均匀奥氏体状态板料,厚向异性系数采用温度相关的逻辑斯蒂方程描述,其等效应力定义为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>&sigma;</mi><mo>&OverBar;</mo></mover><mo>=</mo><msup><mrow><mo>&lsqb;</mo><mi>a</mi><msup><mrow><mo>(</mo><mi>&Lambda;</mi><mo>+</mo><mi>&Gamma;</mi><mo>)</mo></mrow><mrow><mn>2</mn><mi>k</mi></mrow></msup><mo>+</mo><mi>a</mi><msup><mrow><mo>(</mo><mi>&Lambda;</mi><mo>-</mo><mi>&Gamma;</mi><mo>)</mo></mrow><mrow><mn>2</mn><mi>k</mi></mrow></msup><mo>+</mo><mi>b</mi><msup><mrow><mo>(</mo><mi>&Lambda;</mi><mo>+</mo><mi>&Psi;</mi><mo>)</mo></mrow><mrow><mn>2</mn><mi>k</mi></mrow></msup><mo>+</mo><mi>b</mi><msup><mrow><mo>(</mo><mi>&Lambda;</mi><mo>-</mo><mi>&Psi;</mi><mo>)</mo></mrow><mrow><mn>2</mn><mi>k</mi></mrow></msup><mo>&rsqb;</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn><mi>k</mi></mrow></msup></mrow>]]></math><img file="FDA0000750066960000011.GIF" wi="1333" he="75" /></maths>式中Г=Lσ<sub>11</sub>+Mσ<sub>22</sub><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>&Lambda;</mi><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>N&sigma;</mi><mn>11</mn></msub><mo>-</mo><msub><mi>P&sigma;</mi><mn>22</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msub><mi>&sigma;</mi><mn>12</mn></msub><msub><mi>&sigma;</mi><mn>21</mn></msub></mrow></msqrt></mrow>]]></math><img file="FDA0000750066960000012.GIF" wi="690" he="97" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>&Psi;</mi><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>Q&sigma;</mi><mn>11</mn></msub><mo>-</mo><msub><mi>R&sigma;</mi><mn>22</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msub><mi>&sigma;</mi><mn>12</mn></msub><msub><mi>&sigma;</mi><mn>21</mn></msub></mrow></msqrt></mrow>]]></math><img file="FDA0000750066960000013.GIF" wi="702" he="90" /></maths>其中a,b,L,M,N,P,Q和R为各向异性参数,k是材料晶格相关的参数,面心立方晶体k=3,体心立方晶体k=4;由于热成形中奥氏体态下的板料忽略其面内各向异性,因而确定L=M=N=P=Q=R=1/2,a=1/(1+r),b=r/(1+r),其中r为厚向异性系数,由单向拉伸实验中板料的横向应变与厚向应变的比值计算;所述硬化模型为基于微观特征长度相关的内变量演变的唯象本构模型;步骤2:在稳态条件下基于材质不匀沟槽假设的集中性失稳理论确定热成形极限模型,其中稳态条件下的热成形极限模型中将采用确定的稳态平衡条件下的板料本构模型,包括通过平面应变条件下的实验结果确定所述热成形极限模型的系数;步骤3:运用所述热成形极限模型进行瞬态工艺条件下热成形极限的模拟预测,模拟预测用于预测瞬态非平衡条件下的热成形极限,虽然冷却速度是恒定的,但对不同应变路径其达到临界应变状态的时间不同,导致最终温度不同,因此在瞬态成形极限曲线上,需要标注对应不同应变路径在达到临界应变状态时的最终温度。
地址 410082 湖南省长沙市岳麓区麓山南路1号