CN116312902A - 一种确定材料Johnson-Cook动态本构模型的方法 - Google Patents

一种确定材料Johnson-Cook动态本构模型的方法 Download PDF

Info

Publication number
CN116312902A
CN116312902A CN202310594115.7A CN202310594115A CN116312902A CN 116312902 A CN116312902 A CN 116312902A CN 202310594115 A CN202310594115 A CN 202310594115A CN 116312902 A CN116312902 A CN 116312902A
Authority
CN
China
Prior art keywords
johnson
cook
constitutive model
test
dynamic constitutive
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202310594115.7A
Other languages
English (en)
Other versions
CN116312902B (zh
Inventor
戴兰宏
张洧菡
李统
乔禹
陈艳
汪海英
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Mechanics of CAS
Original Assignee
Institute of Mechanics of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Mechanics of CAS filed Critical Institute of Mechanics of CAS
Priority to CN202310594115.7A priority Critical patent/CN116312902B/zh
Publication of CN116312902A publication Critical patent/CN116312902A/zh
Application granted granted Critical
Publication of CN116312902B publication Critical patent/CN116312902B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computing Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明提供了一种确定材料Johnson‑Cook动态本构模型的方法，首先获得待测材料在每组工况下的拉伸试验数据，然后读取待测材料在对应工况下的真实塑形段曲线，确定Johnson‑Cook动态本构模型中的参数，设置收敛阈值并迭代，直至Johnson‑Cook动态本构模型中的参数与上一代参数相比小于等于设定的收敛阈值时判定迭代结束，输出结果，完成拟合，最后对得出的Johnson‑Cook动态本构模型参数进行评估，得出可准确预测该材料的力学行为。本发明构思合理，基于多种高低应变率和温度的五组拉伸试验数据，能够低成本、高效、准确的预测出材料在高应变率的力学行为。

Description

一种确定材料Johnson-Cook动态本构模型的方法
技术领域
本发明涉及材料力学行为的模拟仿真技术领域，具体涉及一种确定材料Johnson-Cook动态本构模型的方法。
背景技术
在汽车、航空航天、军事国防等许多重要领域都涉及材料动态力学行为，其中的关键部件在冲击变形下常常伴随高温、高压等瞬态过程。然而，选取的材料在这一极端条件下是否满足要求是一个亟待解决的问题。相对于传统实验方法的成本高、周期长，通过模拟仿真技术可以快速高效的得到答案，其中的基础和前提是找到一个能够准确描述材料在高应变率和宽温度范围内力学行为的本构关系。
Johnson-Cook（J-C）动态本构模型是Johnson和Cook于1983年建立的一种经验模型。J-C动态本构关系由于形式简单、参数容易标定而被广泛应用，并被嵌入ABAQUS、ANSYS等商用有限元软件中。J-C动态本构关系认为金属的等效流动应力可表示为等效应变函数、等效应变率函数和温度函数的乘积，具体形式为：
Figure SMS_1
（1）；
式（1）中：
Figure SMS_12
为等效流动应力；/>
Figure SMS_4
为等效塑性应变；/>
Figure SMS_8
为无量纲塑性应变率，其中/>
Figure SMS_16
为塑性应变率，/>
Figure SMS_19
为参考应变率；/>
Figure SMS_18
为无量纲温度，其中/>
Figure SMS_20
为当前温度，
Figure SMS_13
为熔化温度，/>
Figure SMS_17
为参考温度；/>
Figure SMS_2
、/>
Figure SMS_7
、/>
Figure SMS_10
、/>
Figure SMS_14
和/>
Figure SMS_11
为材料常数，其中，/>
Figure SMS_15
为参考温度和参考应变率下的真实屈服应力，/>
Figure SMS_3
为应变硬化系数，/>
Figure SMS_6
为应变硬化指数，/>
Figure SMS_5
和/>
Figure SMS_9
分别代表应变率硬化系数和热软化系数的材料常数。该本构关系具有明确的物理意义，等号右边第一个括号内表达式为应变对材料的强化作用，第二个括号内表达式为应变率对材料的强化作用，第三个括号内的表达式为温度对材料的软化作用。
目前，在已发表的文献中可以找到许多改进方法，其中广泛用于标定Johnson-Cook动态本构模型参数的方法可以概括为如下的“三步走”：
1）第一步在参考条件下，式（1）简化为：
Figure SMS_21
（2）；
根据材料在参考试验工况下的真实屈服应力，确定参数
Figure SMS_22
；然后对塑性段曲线采用最小二乘法拟合确定/>
Figure SMS_23
、/>
Figure SMS_24
。
2）第二步在参考温度下，式（1）简化为：
Figure SMS_25
（3）；
其中参数
Figure SMS_26
、/>
Figure SMS_27
、/>
Figure SMS_28
已知，对参考温度下的高应变率真实塑性段曲线采用最小二乘法拟合确定参数/>
Figure SMS_29
。
3）第三步在参考应变率下，式（1）简化为：
Figure SMS_30
（4）；
与第二步相类似，对参考应变率下的高温真实塑性段曲线采用最小二乘法拟合确定参数
Figure SMS_31
。
然而，实践证明这一方法确定的参数能够很好的描述材料准静态力学行为，但是并不能准确预测材料的动态力学行为。分析后可以发现有两方面原因：
①一方面，在“第二步”确定应变率硬化系数
Figure SMS_32
中有时只采用一组动态工况，有时采用的多组动态工况中所涉及应变率范围较小，据此所得到的参数往往不能准确预测材料在大应变率范围下的力学行为；
②另一方面，材料在高应变率或者高温环境下的应力-应变曲线和在室温准静态下的应力-应变曲线往往不是一个常数关系，但是在以上标定过程中的“单向”默认了这一关系。例如在“第二步”中，应变硬化系数
Figure SMS_33
和应变硬化指数/>
Figure SMS_34
已由第一步的参考试验工况确定，再通过动态工况标定应变率硬化系数/>
Figure SMS_35
中自然更倾向于能够拟合上低应变率，而不能很准确的拟合高应变率。
Johnson-Cook动态本构模型是材料爆炸与冲击动力学及相关工程应用领域使用最广泛的本构模型之一。然而，经典的Johnson-Cook动态本构模型参数确定方法是依据一组高应变率和多组准静态应变率实验数据拟合确定的，导致本构理论预测与实验偏差较大。
综上所述，有必要对现有技术做进一步创新。
发明内容
针对上述背景技术中存在的问题，本发明提出了一种确定材料Johnson-Cook动态本构模型的方法，其构思合理，基于多种高低应变率和温度的六组拉伸试验数据，通过在传统“三步走”基础上搭建了反馈通道的方法标定出材料的Johnson-Cook动态本构模型参数，能够低成本、高效、准确的预测出材料在高应变率的力学行为。
为解决上述技术问题，本发明提供的一种确定材料Johnson-Cook动态本构模型的方法，首先获得待测材料在每组工况下的拉伸试验数据，然后读取待测材料在对应工况下的真实塑形段曲线，确定Johnson-Cook动态本构模型中的参数，设置收敛阈值并迭代，直至Johnson-Cook动态本构模型中的参数与上一代参数相比小于等于设定的收敛阈值时判定迭代结束，输出结果，完成拟合，最后对得出的Johnson-Cook动态本构模型参数进行评估，得出可准确预测该材料的力学行为。
所述确定材料Johnson-Cook动态本构模型的方法，其中，所述方法具体包括以下步骤：
1）获得待测材料在六组试验工况下的应力-应变曲线，且六组试验工况分别为：第一组试验工况为应变率1.0×10-3s-1、温度25℃,第二组试验工况为应变率1.0×10-3s-1、温度0.15Tm,第三组试验工况为应变率1.0×10-3s-1、温度0.30Tm,第四组试验工况为应变率5.0×102s-1、温度25℃,第五组试验工况为应变率4.0×103s-1、温度25℃,来标定Johnson-Cook动态本构参数，第六组试验工况为应变率1.5×102s-1、温度25℃；
2)设定第一组试验工况为参考试验工况,在参考试验工况下将下式（1）中的Johnson-Cook动态本构模型简化为式（2），读取待测材料在参考试验工况下的真实屈服应力，确定式（2）中Johnson-Cook动态本构模型中的参数
Figure SMS_36
；读取参考试验工况中的塑形段，通过最小二乘法拟合出Johnson-Cook动态本构模型中的应变硬化系数/>
Figure SMS_37
和应变硬化指数/>
Figure SMS_38
，记为/>
Figure SMS_39
和/>
Figure SMS_40
；
Figure SMS_41
（1）；
Figure SMS_42
（2）；
其中上式（1）中
Figure SMS_44
为等效流动应力,/>
Figure SMS_46
为等效塑性应变；/>
Figure SMS_49
为无量纲塑性应变率，其中/>
Figure SMS_43
为塑性应变率，/>
Figure SMS_47
为参考应变率；/>
Figure SMS_50
为无量纲温度，其中/>
Figure SMS_52
为当前工况中的温度，/>
Figure SMS_45
为材料熔化温度，/>
Figure SMS_48
为参考试验工况中的温度,/>
Figure SMS_51
为参考试验工况下的屈服应力;
3)在参考试验工况的参考温度25℃下，上式（1）中的Johnson-Cook动态本构模型简化为下式（3）;读取待测材料在第四组试验工况和第五组试验工况下的真实应力-应变曲线，对第四组试验工况和第五组试验工况中的两条塑形段采用最小二乘法拟合出Johnson-Cook动态本构模型中的应变率硬化系数
Figure SMS_53
，记为/>
Figure SMS_54
；
Figure SMS_55
（3）;
4）在参考试验工况的参考应变率1.0×10-3s-1下，将上式（1）中的Johnson-Cook动态本构模型简化为下式（4）;读取待测材料在第二组试验工况和第三组试验工况下的真实应力-应变曲线，对第二组试验工况和第三组试验工况中的两条塑形段采用最小二乘法拟合出Johnson-Cook动态本构模型中的热软化系数
Figure SMS_56
，记为/>
Figure SMS_57
；
Figure SMS_58
（4）；
5）读取待测材料在第一组至第五组试验工况下的整条塑形段以及应变率硬化系数
Figure SMS_60
和热软化系数/>
Figure SMS_63
，通过最小二乘法重新对上式(2)中Johnson-Cook动态本构模型中的应变硬化系数/>
Figure SMS_65
和应变硬化指数/>
Figure SMS_61
进行拟合，记为/>
Figure SMS_62
、/>
Figure SMS_64
，并替换上述步骤2)所得到的
Figure SMS_66
、/>
Figure SMS_59
；
6）重复以上步骤3)至步骤5)，设置收敛阈值为1.0×10-6；当上式（1）中Johnson-Cook动态本构模型参数
Figure SMS_83
、/>
Figure SMS_86
、/>
Figure SMS_89
、/>
Figure SMS_68
、/>
Figure SMS_72
迭代至第/>
Figure SMS_75
代得到此时的五个参数值/>
Figure SMS_79
、/>
Figure SMS_69
、/>
Figure SMS_73
、/>
Figure SMS_77
、/>
Figure SMS_81
与第/>
Figure SMS_70
代的参数值/>
Figure SMS_71
、/>
Figure SMS_76
、/>
Figure SMS_80
、/>
Figure SMS_85
、/>
Figure SMS_88
相比均相差小于等于1.0×10-6，判定迭代结束，参数/>
Figure SMS_91
、/>
Figure SMS_93
、/>
Figure SMS_67
、/>
Figure SMS_74
、/>
Figure SMS_78
作为Johnson-Cook动态本构模型参数/>
Figure SMS_82
、/>
Figure SMS_84
、/>
Figure SMS_87
、/>
Figure SMS_90
、/>
Figure SMS_92
的输出结果,完成拟合,得到待测材料的Johnson-Cook动态本构参数;
7）通过第六组试验工况的试验数据对上述步骤1)-6)得出的Johnson-Cook动态本构模型在参考试验工况下的真实屈服应力
Figure SMS_94
、应变硬化系数/>
Figure SMS_95
、应变硬化指数/>
Figure SMS_96
、应变率硬化系数/>
Figure SMS_97
和热软化系数/>
Figure SMS_98
进行评估，得到可以准确预测待测材料的力学行为。
所述确定材料Johnson-Cook动态本构模型的方法，其中，所述步骤1）具体步骤为：通过所述第一组至第五组试验工况来标定Johnson-Cook动态本构参数，通过所述第六组试验工况来评判所标定的Johnson-Cook动态本构参数，以获得待测材料在六组试验工况下的应力-应变曲线。
所述确定材料Johnson-Cook动态本构模型的方法，其中：所述步骤7）在进行评估过程中当平均误差处于15%以内，认为可准确预测该待测材料的力学行为。
采用上述技术方案，本发明具有如下有益效果：
本发明确定材料Johnson-Cook动态本构模型的方法构思合理,其作为一套标准的标定材料Johnson-Cook动态本构模型参数方案，基于多种高低应变率和温度的六组拉伸试验数据，通过在传统“三步走”基础上搭建了反馈通道的方法标定出材料的Johnson-Cook动态本构参数,可以实现准确预测材料在高应变率下的力学行为，且具有低成本、高效、准确特点。
本发明综合利用第二、三组试验工况和第一、四、五组试验工况的实验数据，同时构建一种依据高低应变率实验数据融合且具有反馈机制的整体优化策略,依据新策确定的Johnson-Cook动态本构模型预测与实与实验数据符合令人满意，并获多种材料检验证实。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明确定材料Johnson-Cook动态本构模型的方法中在标定Johnson-Cook动态本构模型参数的传统“三步走”基础上增加迭代的流程图；
图2为本发明确定材料Johnson-Cook动态本构模型的方法中锻造的TC4材料在六组试验工况下的工程应力-应变曲线图；
图3为本发明确定材料Johnson-Cook动态本构模型的方法中3D打印的TC4材料在六组试验工况下的工程应力-应变曲线图；
图4为本发明确定材料Johnson-Cook动态本构模型的方法中锻造的TC4材料的J-C动态本构方程曲线与试验数据曲线对比图；
图5为本发明确定材料Johnson-Cook动态本构模型的方法中3D打印的TC4材料的J-C动态本构方程曲线与试验数据曲线对比图；
图6为本发明确定材料Johnson-Cook动态本构模型的方法中锻造的TC4材料的预测曲线与试验数据曲线对比图；
图7为本发明确定材料Johnson-Cook动态本构模型的方法中3D打印的TC4材料的预测曲线与试验数据曲线对比图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
下面结合具体的实施方式对本发明做进一步的解释说明。
具有双相结构的TC4材料具有高强度、低密度和耐腐蚀等特点，并广泛应用于航空航天、汽车、生物医学等领域。此外，与传统的加工相比，增材制造技术因其制造和修复具有复杂几何形状的产品的能力而受到越来越多的关注。然而，目前尚缺乏对增材制造得到的TC4材料在高应变率和宽温度范围内力学行为的研究。本发明可以低成本、高效、准确的解决这一问题，其中分别对锻造和3D打印得到的TC4材料进行了Johnson-Cook动态本构模型参数标定。
如图1所示，本实施例提供的确定材料Johnson-Cook动态本构模型的方法，首先获得待测材料在每组工况下的拉伸试验数据，然后读取待测材料在对应工况下的真实塑形段曲线，确定Johnson-Cook动态本构模型中的参数，设置收敛阈值并迭代，直至Johnson-Cook动态本构模型中的参数与上一代参数相比小于等于设定的收敛阈值时判定迭代结束，输出结果，完成拟合，最后对得出的Johnson-Cook动态本构模型参数进行评估，得出可准确预测该材料的力学行为。
本发明确定材料Johnson-Cook动态本构模型的方法，具体包括以下步骤：
S100、分别获得锻造和3D打印得到的TC4材料在表1所示的六组试验工况下的应力-应变曲线。其中通过表1中的前五组试验工况即第一组试验工况为应变率1.0×10-3s-1、温度25℃,第二组试验工况为应变率1.0×10-3s-1、温度0.15Tm,第三组试验工况为应变率1.0×10-3s-1、温度0.30Tm,第四组试验工况为应变率5.0×102s-1、温度25℃,第五组试验工况为应变率4.0×103s-1、温度25℃,来标定Johnson-Cook动态本构参数;第六组试验工况为应变率1.5×102s-1、温度25℃，用于评判所标定的Johnson-Cook动态本构参数；两种状态材料的熔点相近均为1665℃；得到的两种状态材料在六组试验工况下的工程应力-应变曲线如图2和图3所示。
表1用于标定并评判Johnson-Cook动态本构参数的六组试验
Figure SMS_99
S200、设定第一组试验工况为参考试验工况,在参考试验工况下将下式（1）中的Johnson-Cook动态本构模型简化为下式（2），读取待测材料在参考试验工况下的真实屈服应力，确定式（2）中Johnson-Cook动态本构模型中的参数
Figure SMS_100
；读取参考试验工况中的塑形段，通过最小二乘法拟合出Johnson-Cook动态本构模型中的应变硬化系数/>
Figure SMS_101
和应变硬化指数/>
Figure SMS_102
，记为/>
Figure SMS_103
和/>
Figure SMS_104
；
Figure SMS_105
（1）；
Figure SMS_106
（2）；
其中上式（1）中
Figure SMS_107
为等效流动应力；/>
Figure SMS_110
为等效塑性应变；/>
Figure SMS_113
为无量纲塑性应变率，其中/>
Figure SMS_109
为塑性应变率，/>
Figure SMS_111
为参考应变率；/>
Figure SMS_114
为无量纲温度，其中/>
Figure SMS_116
为当前工况中的温度，/>
Figure SMS_108
为材料熔化温度，/>
Figure SMS_112
为参考试验工况中的温度,/>
Figure SMS_115
为参考试验工况下的屈服应力。
S300、在参考试验工况的参考温度25℃下,将上式（1）中的Johnson-Cook动态本构模型简化为下式（3）;读取待测材料在第四组试验工况和第五组试验工况下的真实应力-应变曲线，对第四组试验工况和第五组试验工况中的两条塑形段采用最小二乘法拟合出Johnson-Cook动态本构模型中的应变率硬化系数
Figure SMS_117
，记为/>
Figure SMS_118
；
Figure SMS_119
（3）。
S400、在参考试验工况的参考应变率1.0×10-3s-1下，将上式（1）中的Johnson-Cook动态本构模型简化为下式（4）;读取待测材料在第二组试验工况和第三组试验工况下的真实应力-应变曲线，对第二组试验工况和第三组试验工况中的两条塑形段采用最小二乘法拟合出Johnson-Cook动态本构模型中的热软化系数
Figure SMS_120
，记为/>
Figure SMS_121
；
Figure SMS_122
（4）。
S500、读取待测材料在表1所示的前五组即第一组至第五组试验工况下的整条塑形段以及应变率硬化系数
Figure SMS_124
和热软化系数/>
Figure SMS_126
，通过最小二乘法重新对式(2)中Johnson-Cook动态本构模型中的应变硬化系数/>
Figure SMS_128
和应变硬化指数/>
Figure SMS_125
进行拟合，记为/>
Figure SMS_127
、/>
Figure SMS_129
，并替换上述步骤S200所得到的/>
Figure SMS_130
、/>
Figure SMS_123
。
S600、重复以上步骤S300至步骤S500，设置收敛阈值为1.0×10-6；当式（1）中Johnson-Cook动态本构模型参数
Figure SMS_148
、/>
Figure SMS_151
、/>
Figure SMS_154
、/>
Figure SMS_132
、/>
Figure SMS_136
迭代至第/>
Figure SMS_140
代，得到此时的五个参数值
Figure SMS_144
、/>
Figure SMS_134
、/>
Figure SMS_138
、/>
Figure SMS_141
、/>
Figure SMS_145
与第/>
Figure SMS_147
代的参数值/>
Figure SMS_150
、/>
Figure SMS_153
、/>
Figure SMS_156
、/>
Figure SMS_149
、/>
Figure SMS_152
相比均相差小于等于1.0×10-6，判定迭代结束，参数/>
Figure SMS_155
、/>
Figure SMS_157
、/>
Figure SMS_131
、/>
Figure SMS_137
、/>
Figure SMS_142
为作为Johnson-Cook动态本构模型参数/>
Figure SMS_146
、/>
Figure SMS_133
、/>
Figure SMS_135
、/>
Figure SMS_139
、/>
Figure SMS_143
的输出结果，J-C动态本构方程曲线与试验数据曲线如图4、5所示，表2为完成拟合后锻造和3D打印得到TC4材料的Johnson-Cook动态本构参数。
表2标定出两种材料的JOHNSON-COOK动态本构参数取值
Figure SMS_158
S700、通过第六组试验工况的试验数据对步骤S100-S600得出的Johnson-Cook动态本构模型在参考试验工况下的真实屈服应力
Figure SMS_159
、应变硬化系数/>
Figure SMS_160
、应变硬化指数/>
Figure SMS_161
、应变率硬化系数/>
Figure SMS_162
和热软化系数/>
Figure SMS_163
进行评估，预测曲线与试验数据曲线如图6、7所示，表3为计算得到的平均误差，均处于15%以内，认为所得到的Johnson-Cook动态本构参数可以准确预测该待测材料的力学行为。
表3两种材料的J-C动态本构方程预测曲线与试验数据曲线的平均误差。
Figure SMS_164
本发明构思合理，基于多种高低应变率和温度的六组拉伸试验数据，通过在传统“三步走”基础上搭建了反馈通道的方法标定出材料的Johnson-Cook动态本构参数，能够低成本、高效、准确的预测出材料在高应变率的力学行为。
本发明包括但不局限于以上的工程问题，凡是在本发明的精神和原则之下，进行的任何局部改进、修改，都将视为在本发明的保护范围之内。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (4)

1.一种确定材料Johnson-Cook动态本构模型的方法，其特征在于，首先获得待测材料在每组工况下的拉伸试验数据，然后读取待测材料在对应工况下的真实塑形段曲线，确定Johnson-Cook动态本构模型中的参数，设置收敛阈值并迭代，直至Johnson-Cook动态本构模型中的参数与上一代参数相比小于等于设定的收敛阈值时判定迭代结束，输出结果，完成拟合，最后对得出的Johnson-Cook动态本构模型参数进行评估，得出可准确预测该材料的力学行为。
2.如权利要求1所述的确定材料Johnson-Cook动态本构模型的方法，其特征在于，所述方法具体包括以下步骤：
1）获得待测材料在六组试验工况下的应力-应变曲线，且六组试验工况分别为：第一组试验工况为应变率1.0×10-3s-1、温度25℃,第二组试验工况为应变率1.0×10-3s-1、温度0.15Tm,第三组试验工况为应变率1.0×10-3s-1、温度0.30Tm,第四组试验工况为应变率5.0×102s-1、温度25℃,第五组试验工况为应变率4.0×103s-1、温度25℃,来标定Johnson-Cook动态本构参数，第六组试验工况为应变率1.5×102s-1、温度25℃；
2)设定第一组试验工况为参考试验工况,在参考试验工况下将下式（1）中的Johnson-Cook动态本构模型简化为式（2），读取待测材料在参考试验工况下的真实屈服应力，确定式（2）中Johnson-Cook动态本构模型中的参数
Figure QLYQS_1
；读取参考试验工况中的塑形段，通过最小二乘法拟合出Johnson-Cook动态本构模型中的应变硬化系数/>
Figure QLYQS_2
和应变硬化指数/>
Figure QLYQS_3
，记为/>
Figure QLYQS_4
和/>
Figure QLYQS_5
；
Figure QLYQS_6
（1）；
Figure QLYQS_7
（2）；
其中上式（1）中
Figure QLYQS_9
为等效流动应力,/>
Figure QLYQS_12
为等效塑性应变；/>
Figure QLYQS_14
为无量纲塑性应变率，其中/>
Figure QLYQS_10
为塑性应变率，/>
Figure QLYQS_13
为参考应变率；/>
Figure QLYQS_16
为无量纲温度，其中/>
Figure QLYQS_17
为当前工况中的温度，/>
Figure QLYQS_8
为材料熔化温度，/>
Figure QLYQS_11
为参考试验工况中的温度,/>
Figure QLYQS_15
为参考试验工况下的屈服应力;
3)在参考试验工况的参考温度25℃下，上式（1）中的Johnson-Cook动态本构模型简化为下式（3）;读取待测材料在第四组试验工况和第五组试验工况下的真实应力-应变曲线，对第四组试验工况和第五组试验工况中的两条塑形段采用最小二乘法拟合出Johnson-Cook动态本构模型中的应变率硬化系数
Figure QLYQS_18
，记为/>
Figure QLYQS_19
；
Figure QLYQS_20
（3）;
4）在参考试验工况的参考应变率1.0×10-3s-1下，将上式（1）中的Johnson-Cook动态本构模型简化为下式（4）;读取待测材料在第二组试验工况和第三组试验工况下的真实应力-应变曲线，对第二组试验工况和第三组试验工况中的两条塑形段采用最小二乘法拟合出Johnson-Cook动态本构模型中的热软化系数
Figure QLYQS_21
，记为/>
Figure QLYQS_22
；
Figure QLYQS_23
（4）；
5）读取待测材料在第一组至第五组试验工况下的整条塑形段以及应变率硬化系数
Figure QLYQS_26
和热软化系数/>
Figure QLYQS_28
，通过最小二乘法重新对上式(2)中Johnson-Cook动态本构模型中的应变硬化系数/>
Figure QLYQS_30
和应变硬化指数/>
Figure QLYQS_24
进行拟合，记为/>
Figure QLYQS_27
、/>
Figure QLYQS_29
，并替换上述步骤2)所得到的/>
Figure QLYQS_31
、/>
Figure QLYQS_25
；
6）重复以上步骤3)至步骤5)，设置收敛阈值为1.0×10-6；当上式（1）中Johnson-Cook动态本构模型参数
Figure QLYQS_49
、/>
Figure QLYQS_52
、/>
Figure QLYQS_54
、/>
Figure QLYQS_33
、/>
Figure QLYQS_37
迭代至第/>
Figure QLYQS_41
代得到此时的五个参数值/>
Figure QLYQS_45
、/>
Figure QLYQS_36
、/>
Figure QLYQS_39
、/>
Figure QLYQS_42
、
Figure QLYQS_46
与第/>
Figure QLYQS_34
代的参数值/>
Figure QLYQS_38
、/>
Figure QLYQS_43
、/>
Figure QLYQS_47
、/>
Figure QLYQS_48
、/>
Figure QLYQS_51
相比均相差小于等于1.0×10-6，判定迭代结束，参数/>
Figure QLYQS_55
、/>
Figure QLYQS_57
、/>
Figure QLYQS_32
、/>
Figure QLYQS_35
、/>
Figure QLYQS_40
作为Johnson-Cook动态本构模型参数/>
Figure QLYQS_44
、/>
Figure QLYQS_50
、
Figure QLYQS_53
、/>
Figure QLYQS_56
、/>
Figure QLYQS_58
的输出结果,完成拟合,得到待测材料的Johnson-Cook动态本构参数;
7）通过第六组试验工况的试验数据对上述步骤1)-6)得出的Johnson-Cook动态本构模型在参考试验工况下的真实屈服应力
Figure QLYQS_59
、应变硬化系数/>
Figure QLYQS_60
、应变硬化指数/>
Figure QLYQS_61
、应变率硬化系数/>
Figure QLYQS_62
和热软化系数/>
Figure QLYQS_63
进行评估，得到可以准确预测待测材料的力学行为。
3.如权利要求2所述的确定材料Johnson-Cook动态本构模型的方法，其特征在于，所述步骤1）具体步骤为：通过所述第一组至第五组试验工况来标定Johnson-Cook动态本构参数，通过所述第六组试验工况来评判所标定的Johnson-Cook动态本构参数，以获得待测材料在六组试验工况下的应力-应变曲线。
4.如权利要求2所述的确定材料Johnson-Cook动态本构模型的方法，其特征在于：所述步骤7）在进行评估过程中当平均误差处于15%以内，认为可准确预测该待测材料的力学行为。
CN202310594115.7A 2023-05-25 2023-05-25 一种确定材料Johnson-Cook动态本构模型的方法 Active CN116312902B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310594115.7A CN116312902B (zh) 2023-05-25 2023-05-25 一种确定材料Johnson-Cook动态本构模型的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310594115.7A CN116312902B (zh) 2023-05-25 2023-05-25 一种确定材料Johnson-Cook动态本构模型的方法

Publications (2)

Publication Number Publication Date
CN116312902A true CN116312902A (zh) 2023-06-23
CN116312902B CN116312902B (zh) 2023-08-04

Family

ID=86785508

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310594115.7A Active CN116312902B (zh) 2023-05-25 2023-05-25 一种确定材料Johnson-Cook动态本构模型的方法

Country Status (1)

Country Link
CN (1) CN116312902B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117454673A (zh) * 2023-12-22 2024-01-26 中南大学 考虑降温速率的金属材料本构模型构建方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111554362A (zh) * 2020-03-30 2020-08-18 西安建筑科技大学 一种建立hbprc动态损伤本构模型的方法
WO2021037036A1 (zh) * 2019-08-28 2021-03-04 肖锋 一种测定材料真实应力应变曲线的试验与计算方法
WO2022077204A1 (zh) * 2020-10-13 2022-04-21 华东理工大学 高温结构强度和寿命分析的时间相关局部应力应变方法及工具软件
CN115831295A (zh) * 2023-01-06 2023-03-21 中铝材料应用研究院有限公司 材料本构方程参数标定方法、装置及计算机设备
CN116030923A (zh) * 2023-03-28 2023-04-28 深圳大学 获取材料动态本构关系的方法、装置、设备及存储介质

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021037036A1 (zh) * 2019-08-28 2021-03-04 肖锋 一种测定材料真实应力应变曲线的试验与计算方法
CN111554362A (zh) * 2020-03-30 2020-08-18 西安建筑科技大学 一种建立hbprc动态损伤本构模型的方法
WO2022077204A1 (zh) * 2020-10-13 2022-04-21 华东理工大学 高温结构强度和寿命分析的时间相关局部应力应变方法及工具软件
CN115831295A (zh) * 2023-01-06 2023-03-21 中铝材料应用研究院有限公司 材料本构方程参数标定方法、装置及计算机设备
CN116030923A (zh) * 2023-03-28 2023-04-28 深圳大学 获取材料动态本构关系的方法、装置、设备及存储介质

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
DING-NI ZHANG 等: "A modified Johnson–Cook model of dynamic tensile behaviors for 7075-T6 aluminum alloy", 《JOURNAL OF ALLOYS AND COMPOUNDS》；, vol. 619, pages 186 - 194 *
M BURLEY 等: "Johnson-Cook parameter evaluation from ballistic impact data via iterative FEM modelling", 《INTERNATIONAL JOURNAL OF IMPACT ENGINEERING》, vol. 112, pages 180 - 192, XP085269507, DOI: 10.1016/j.ijimpeng.2017.10.012 *
YU QIAO 等: "Enhancing the dynamic temperature stability of epoxy with graphene oxide", 《MECHANICS OF MATERIALS》, vol. 150, pages 103593 *
YUNTIAN WANG 等: "Modified Johnson-Cook constitutive model of metallic materials under a wide range of temperatures and strain rates", 《RESULTS IN PHYSICS》, vol. 27, pages 104498 *
杨曾宇 等: "Ti-6Al-4V合金动态本构模型参数敏感性及其优化", 《钛工业进展》, vol. 34, no. 5, pages 12 - 17 *
柳爱群 等: "高应变率变形的Johnson—Cook动态本构模型参数识别方法", 《应力数学和力学》, vol. 35, no. 2, pages 219 - 225 *
鄢阿敏 等: "高熵合金药型罩射流成型与稳定性", 《力学学报》, vol. 54, no. 8, pages 2119 - 2130 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117454673A (zh) * 2023-12-22 2024-01-26 中南大学 考虑降温速率的金属材料本构模型构建方法及装置
CN117454673B (zh) * 2023-12-22 2024-03-12 中南大学 考虑降温速率的金属材料本构模型构建方法及装置

Also Published As

Publication number Publication date
CN116312902B (zh) 2023-08-04

Similar Documents

Publication Publication Date Title
CN116312902B (zh) 一种确定材料Johnson-Cook动态本构模型的方法
US8346693B2 (en) Method for hammerstein modeling of steam generator plant
CN104809273B (zh) 一种描述蠕变变形的方法
Cho et al. A study on the hot-deformation behavior and dynamic recrystallization of Al–5 wt.% Mg alloy
Zhang et al. Identification of friction coefficients and strain-compensated Arrhenius-type constitutive model by a two-stage inverse analysis technique
CN113673030B (zh) 一种金属材料韧性断裂耦合失效仿真分析方法
CN111310348A (zh) 一种基于pso-lssvm的材料本构模型预测方法
CN113609719B (zh) 一种高强钛管数控弯曲滞后回弹补偿控制的研究方法
CN115034114B (zh) 一种基于数值模拟的挤压有限元仿真模型优化方法
CN112347574B (zh) 一种增材制造金属材料的弹塑性损伤失效预测方法
CN107633106A (zh) 一种基于全局差分法的非均匀温度场下热模态灵敏度分析方法
CN106446390A (zh) 金属材料稳态蠕变速率拟合方程的计算方法
CN116611348A (zh) 基于机器学习的金属材料统一无参数疲劳寿命预测方法
CN110488888B (zh) 一种基于自适应迭代学习的电阻加热炉温度控制方法
JP2003194686A (ja) 応力−ひずみ関係シミュレート方法および除荷過程における降伏点を求める方法
CN111723513A (zh) 一种通过机器学习神经网络进行反演模拟参数的方法
KR102457539B1 (ko) 구조물 성능기반 재료의 역설계 방법
JP2009125768A (ja) 温度分布履歴推定方法
CN112651153B (zh) 一种确定晶体塑性有限元模型材料参数的方法
CN107363142B (zh) 复合金属板的热冲压成型控制方法
JP2007265266A (ja) 粘弾性応答性能予測方法、ゴム製品設計方法、及び粘弾性応答性能予測装置
JP2001349883A (ja) 金属材料の特性予測方法
CN116720283A (zh) 融合径向基函数与克里金模型的高维代理模型构建方法
CN113935164A (zh) 一种考虑晶格长度的金属材料本构方程拟合方法
CN109408926B (zh) 求解复杂结构多维瞬态非线性热传导反问题的方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant