各向异性屈服准则的UMAT子程序二次开发研究
|
摘要
通过ABAQUS提供的UMAT子程序二次开发接口,结合完全隐式向后Euler图形返回算法,用Fortran语言编程,将Yld2004-18p各向异性屈服准则本构模型嵌入ABAQUS软件,并提出了一个通用的、柔性的本构模型二次开发结构模式。将Yld2004-18p屈服准则退化到Mises各向同性屈服准则,用于实例有限元分析,计算结果与ABAQUS自带的Mises屈服准则计算结果做比较,验证了二次开发过程的正确性,为后续嵌入更高级、更强健、更特殊的屈服准则提供研究思路和方法,使之用于精密塑性有限元分析。
The secondary development interface of UMAT subroutine provided by ABAQUS is combined with the complete implicit backward Euler graph return algorithm. The Fortran language programming has been used to embed the Yld2004-18 p anisotropic yield criterion constitutive model into ABAQUS software. A general and flexible secondary development structure model for constitutive model has been put forward. The Yld2004-18 p yield criterion has been degraded to the Mises isotropic yield criterion for the example finite element analysis. The calculation results have been compared with the calculation results of the Mises yield criterion of ABAQUS. The correctness of the secondary development process has been verified for subsequent embedding. It provides research ideas and methods for subsequent embedding the more advanced, robust and specific yield criteria in order for the application of precision plastic finite element analysis.
The secondary development interface of UMAT subroutine provided by ABAQUS is combined with the complete implicit backward Euler graph return algorithm. The Fortran language programming has been used to embed the Yld2004-18 p anisotropic yield criterion constitutive model into ABAQUS software. A general and flexible secondary development structure model for constitutive model has been put forward. The Yld2004-18 p yield criterion has been degraded to the Mises isotropic yield criterion for the example finite element analysis. The calculation results have been compared with the calculation results of the Mises yield criterion of ABAQUS. The correctness of the secondary development process has been verified for subsequent embedding. It provides research ideas and methods for subsequent embedding the more advanced, robust and specific yield criteria in order for the application of precision plastic finite element analysis.
引文
[1]庄茁,张帆,岑松,等.ABAQUS非线性有限元分析与实例[M].北京:科学出版社,2005.
[2]蒋向华,郑金桥,王义林,等.基于有限元分析的覆盖件拉延筋设计与优化[J].锻压装备与制造技术,2005,40(2):77-80.
[3]黄志超,包忠诩,黄克坚,等.基于ANSYS二次开发的三维金属成形数值模拟研究[J].锻压装备与制造技术,2004,39(1):73-75.
[4]王晓东,余心宏.TC4钛合金U形带筋板件超塑性成形的有限元数值模拟[J].锻压装备与制造技术,2010,45(6):79-82.
[5]李平,庄茁,李殿中,等.板带钢热连轧及微观组织演变非线性三维仿真[J].力学与实践,2002,24:35-38.
[6]杨曼娟.ABAQUS用户材料子程序开发及应用[D].武汉:华中科技大学,2005.
[7]唐榕蔚.统一强度理论的ABAQUS二次开发及其应用[D].北京:北京交通大学,2012.
[8]乔顺成,吴建军,郭绍唱,等.基于PCL语言的螺栓连接结构有限元分析系统研究[J].锻压装备与制造技术,2016,51(1):118-122.
[9]王平,崔礼春,马国礼,等.铝合金顶盖充液成形工艺研究[J].锻压装备与制造技术,2017,52(5):58-61.
[10]吴建军,周维贤.板料成形性基本理论[M].西安:西北工业大学出版社,2010.
[11]倪向贵,吴恒安,王宇,等.各向异性本构关系在板料成形数值模拟中的应用[J].计算力学学报,2003,20(2):231-235.
[12]Barlat F,Aretz H,Yoon J W,et al.Linear transfomation-based anisotropic yield functions[J].International Journal of Plasticity,2005,21(5):1009-1039.
[13]Belytschko T,Liu W K,Moran B,et al.Nonlinear finite elements for continua and structures[M].John wiley&sons,2013.
[2]蒋向华,郑金桥,王义林,等.基于有限元分析的覆盖件拉延筋设计与优化[J].锻压装备与制造技术,2005,40(2):77-80.
[3]黄志超,包忠诩,黄克坚,等.基于ANSYS二次开发的三维金属成形数值模拟研究[J].锻压装备与制造技术,2004,39(1):73-75.
[4]王晓东,余心宏.TC4钛合金U形带筋板件超塑性成形的有限元数值模拟[J].锻压装备与制造技术,2010,45(6):79-82.
[5]李平,庄茁,李殿中,等.板带钢热连轧及微观组织演变非线性三维仿真[J].力学与实践,2002,24:35-38.
[6]杨曼娟.ABAQUS用户材料子程序开发及应用[D].武汉:华中科技大学,2005.
[7]唐榕蔚.统一强度理论的ABAQUS二次开发及其应用[D].北京:北京交通大学,2012.
[8]乔顺成,吴建军,郭绍唱,等.基于PCL语言的螺栓连接结构有限元分析系统研究[J].锻压装备与制造技术,2016,51(1):118-122.
[9]王平,崔礼春,马国礼,等.铝合金顶盖充液成形工艺研究[J].锻压装备与制造技术,2017,52(5):58-61.
[10]吴建军,周维贤.板料成形性基本理论[M].西安:西北工业大学出版社,2010.
[11]倪向贵,吴恒安,王宇,等.各向异性本构关系在板料成形数值模拟中的应用[J].计算力学学报,2003,20(2):231-235.
[12]Barlat F,Aretz H,Yoon J W,et al.Linear transfomation-based anisotropic yield functions[J].International Journal of Plasticity,2005,21(5):1009-1039.
[13]Belytschko T,Liu W K,Moran B,et al.Nonlinear finite elements for continua and structures[M].John wiley&sons,2013.