工程科学学报，第37卷，第8期：1076-1083,2015年8月 Chinese Journal of Engineering,Vol.37,No.8:1076-1083,August 2015 D0l:10.13374/j.issn2095-9389.2015.08.016:http://journals..ustb.edu.cn 耦合孪生晶体塑性模型硬化参数的灵敏度分析 杨 竞”，黄杰12》，孙朝阳)四，郭宁”，王斌” 1)北京科技大学机械工程学院，北京1000832)宝山钢铁股份有限公司，上海201900 ☒通信作者，Email:suney@usth.cdu.cn 摘要基于经典晶体塑性理论，建立了耦合孪生的晶体塑性本构模型并进行了全隐式积分的数值实现.该本构模型采用 饱和硬化法则，并采用孪生阻力与滑移硬化之间的正比关系来描述孪生对滑移硬化影响及孪生硬化行为，针对该本构模型 的13个参数，结合各参数物理意义提出了参数的分类确定方法.以孪生诱导塑性(TWIP)钢Fe-22M0.6C为例，着重对硬 化参数的局部灵敏度进行了分析，研究了各硬化参数对宏观力学响应、孪生激活和演化的影响，根据变形机制的不同宏观变 形过程可区分为孪生硬化阶段和孪生硬化失效阶段，进而给出了硬化参数确定的步骤及其建议取值范围.结果表明：初始滑 移阻力与屈服极限线性相关，取值范围在80~160MP之间：孪生硬化指数增大使得孪生硬化阶段减弱，其取值范围应在0~ 3之间：李生阻力与滑移阻力比值增大，则孪生增长率降低，硬化率拐点后移，直至拐点消失，其取值范围在1~1.3之间 关键词孪生：晶体塑性：本构模型：灵敏度分析：硬化参数 分类号TG302 Sensitivity analysis of hardening parameters in the crystal plasticity model exhibiting deformation twinning YANG Jing,HUANG Jie,SUN Chao-yang,GUO Ning,WANG Bin 1)School of Mechanical Engineering,University of Science and Technology Beijing,Beijing 100083,China 2)Baoshan Iron Steel Co.Ltd.Shanghai 201900,China Corresponding author,E-mail:suncy@ustb.edu.cn ABSTRACT A crystal plasticity model exhibiting deformation twinning is introduced based on the classical crystal plasticity theory, and its numerical implementation is conducted in which a fully implicit integration procedure is employed.In this constitutive model, the saturation-type hardening law is adopted and the direct proportion relationship between twin resistance and slip hardening is used to describe the effect of twinning on slip hardening and twin hardening.By considering the physical meaning of all the 13 parameters in this model,the classification methodology is presented for these parameters.Taking Fe-22Mn-0.6C twinning induced plasticity (TWIP)steel as an example,a local sensitivity analysis of hardening parameters is investigated emphatically.The effects of hardening parameters on the macro-mechanical response,the activation and evolution of twinning and the strategy for determining the hardening parameters are discussed.According to the difference of deformation mechanisms,the macroscopic deformation are divided into a twin hardening stage and a twin hardening failure stage.And then the determination method and the suggestive ranges for the hardening parameters are given.It is found that the initial slip resistance is linearly associated with the yield limit,and the value of initial resist- ance stress ranges from 80 MPa to 160 MPa.The twin hardening stage weakens when the twin hardening index increases,the range of values allowed for the twin hardening index is 0 to 3.When the ratio between twinning resistance and slip resistance increases,the 收稿日期：2015-01-27 基金项目：国家自然科学基金资助项目(51105029,51575039)：国家自然科学基金委员会-中国工程物理研究院联合基金资助项目 (U1330121):非线性力学国家重点实验室开放基金资助项目(LNM201512):北京市自然科学基金资助项目(3112019)：北京市教委重点基金 资助项目(KZ201010005002)工程科学学报，第 37 卷，第 8 期: 1076--1083，2015 年 8 月 Chinese Journal of Engineering，Vol． 37，No． 8: 1076--1083，August 2015 DOI: 10． 13374 /j． issn2095--9389． 2015． 08． 016; http: / /journals． ustb． edu． cn 耦合孪生晶体塑性模型硬化参数的灵敏度分析 杨 竞1) ，黄 杰1，2) ，孙朝阳1) ，郭 宁1) ，王 斌1) 1) 北京科技大学机械工程学院，北京 100083 2) 宝山钢铁股份有限公司，上海 201900  通信作者，E-mail: suncy@ ustb． edu． cn 摘 要 基于经典晶体塑性理论，建立了耦合孪生的晶体塑性本构模型并进行了全隐式积分的数值实现． 该本构模型采用 饱和硬化法则，并采用孪生阻力与滑移硬化之间的正比关系来描述孪生对滑移硬化影响及孪生硬化行为． 针对该本构模型 的 13 个参数，结合各参数物理意义提出了参数的分类确定方法． 以孪生诱导塑性( TWIP) 钢 Fe--22Mn--0. 6C 为例，着重对硬 化参数的局部灵敏度进行了分析，研究了各硬化参数对宏观力学响应、孪生激活和演化的影响，根据变形机制的不同宏观变 形过程可区分为孪生硬化阶段和孪生硬化失效阶段，进而给出了硬化参数确定的步骤及其建议取值范围． 结果表明: 初始滑 移阻力与屈服极限线性相关，取值范围在 80 ～ 160 MPa 之间; 孪生硬化指数增大使得孪生硬化阶段减弱，其取值范围应在 0 ～ 3 之间; 孪生阻力与滑移阻力比值增大，则孪生增长率降低，硬化率拐点后移，直至拐点消失，其取值范围在 1 ～ 1. 3 之间． 关键词 孪生; 晶体塑性; 本构模型; 灵敏度分析; 硬化参数 分类号 TG302 收稿日期: 2015--01--27 基金项目: 国家自然科学基金资助项目( 51105029，51575039 ) ; 国家自然科学基金委员会--中国工程物理研究院联合基金资助项目 ( U1330121) ; 非线性力学国家重点实验室开放基金资助项目( LNM201512) ; 北京市自然科学基金资助项目( 3112019) ; 北京市教委重点基金 资助项目( KZ201010005002) Sensitivity analysis of hardening parameters in the crystal plasticity model exhibiting deformation twinning YANG Jing1) ，HUANG Jie1，2) ，SUN Chao-yang1)  ，GUO Ning1) ，WANG Bin1) 1) School of Mechanical Engineering，University of Science and Technology Beijing，Beijing 100083，China 2) Baoshan Iron ＆ Steel Co． Ltd． ，Shanghai 201900，China  Corresponding author，E-mail: suncy@ ustb． edu． cn ABSTＲACT A crystal plasticity model exhibiting deformation twinning is introduced based on the classical crystal plasticity theory， and its numerical implementation is conducted in which a fully implicit integration procedure is employed． In this constitutive model， the saturation-type hardening law is adopted and the direct proportion relationship between twin resistance and slip hardening is used to describe the effect of twinning on slip hardening and twin hardening． By considering the physical meaning of all the 13 parameters in this model，the classification methodology is presented for these parameters． Taking Fe--22Mn--0. 6C twinning induced plasticity ( TWIP) steel as an example，a local sensitivity analysis of hardening parameters is investigated emphatically． The effects of hardening parameters on the macro-mechanical response，the activation and evolution of twinning and the strategy for determining the hardening parameters are discussed． According to the difference of deformation mechanisms，the macroscopic deformation are divided into a twin hardening stage and a twin hardening failure stage． And then the determination method and the suggestive ranges for the hardening parameters are given． It is found that the initial slip resistance is linearly associated with the yield limit，and the value of initial resist￾ance stress ranges from 80 MPa to 160 MPa． The twin hardening stage weakens when the twin hardening index increases，the range of values allowed for the twin hardening index is 0 to 3． When the ratio between twinning resistance and slip resistance increases，the