(19)中华人民共和国国家知识产权局
(12)发明专利申请
(10)申请公布号 (43)申请公布日 (21)申请号 202011426575.1
(22)申请日 2020.12.09
(71)申请人 中南大学
地址 410083 湖南省长沙市岳麓区麓山南
路932号
(72)发明人 蔺永诚 姜玉强 王冠强
(51)Int.Cl.
G06F 30/23(2020.01)
G06F 30/27(2020.01)
G16C 60/00(2019.01)
G06F 119/14(2020.01)
(54)发明名称一种确定晶体塑性有限元模型材料参数的方法(57)摘要本发明公开了一种确定晶体塑性有限元模型材料参数的方法。该方法如下：(1)设置不同的材料参数组合A；(2)采用参数组合A，进行晶体塑性有限元计算，得到不同参数下的真应力真应变；(3)将参数组合A作为输入，把第二步得到的真应力真应变作为输出，建立深度置信神经网络模型；(4)选取额外的材料参数组合B，进行晶体塑性有限元和深度置信神经网络模型计算，得到不同模型的真应力真应变，并进行比较，验证深度置信神经网络模型的准确性；(5)设置多组材料参数组合C，进行深度置信网络模型计算，将计算得到的真应力真应变和试验值比较，若误差大于设定值，继续修改材料参数组合C，直至误差小于设定值，该参数组合C即为晶体塑性有限元模 型的材料参数。权利要求书1页 说明书4页 附图1页CN 112651153 A 2021.04.13
C N 112651153
A
1.一种确定晶体塑性有限元模型材料参数的方法，其特征在于：该方法通过深度置信神经网络模型，采用较少的晶体塑性有限元计算量，获得准确的晶体塑性有限元模型的材料参数，包括如下步骤：
步骤1：设置不同的材料参数组合A，进行晶体塑性有限元计算，得到不同参数下的真应力真应变；
本步骤所述的晶体塑性有限元模型为唯象学幂率模型：式中：为第α个滑移系剪切应变率，是参考剪切应变率(数值一般为0.001)，τα为第α个滑移系的分解剪切应力，为临界剪切应力，m是滑移速率敏感系数，为滑移系初始强度，h αβ
为滑移硬化模量，表征滑移系β的单位剪切变形对滑移系α滑移阻力的影响，为第β个滑移系剪切应变率，h 0为自硬化模量，τs
是滑移系饱和强度，τβ是滑移系β的分解剪切应力，a为应力指数，q αβ
是描述自硬化系数和潜在硬化系数的矩阵，
可以表示为：其中，q是潜在硬化系数与自硬化系数的比值，对共面滑移系，比值为1，其他情况为
1.4，A为三阶单位矩阵，参数m，h αβ
和a为待确定的材料参数；步骤2：将参数组合A作为输入，把第二步得到的真应力真应变作为输出，建立深度置信神经网络模型；
步骤3：选取额外的材料参数组合B，进行晶体塑性有限元和深度置信神经网络模型计算，得到不同模型的真应力真应变，并进行比较，以验证深度置信神经网络模型的准确性；
步骤4：设置多组材料参数组合C，进行深度置信网络模型计算，将计算得到的真应力真应变和试验值比较，若误差大于设定值，继续修改材料参数组合C，直至误差小于设定值，该参数组合C即为晶体塑性有限元模型的材料参数。
权 利 要 求 书1/1页CN 112651153 A
一种确定晶体塑性有限元模型材料参数的方法
技术领域：
[0001]本发明属于金属材料加工工程技术领域，涉及一种确定晶体塑性有限元模型材料参数的方法。
背景技术：
[0002]高强韧钛合金具有比强度高、冲击韧性好、抗疲劳和抗腐蚀能力强等优点，被广泛用来制造飞机起落架，航空发动机叶片等关键结构件。通过热变形，β相因动态再结晶而细化，进而可以得到细小均匀的组织，这有利于提高材料的力学性能。然而，由于钛合金中β相的取向不同，导致材料会发生不均匀变形和不均匀的再结晶现象。
[0003]晶体塑性有限元是模拟材料不均匀变形及织构演变的重要方法。而晶体塑性有限元模拟的计算周期长，材料参数多，这限制了该方法的应用。传统的确定晶体塑性有限元模型材料参数的方法为试错法。该方法具有不确定性、迭代次数多、计算量极大、收敛性差，且不容易得到精确的材料参数值的缺陷。因此，需要一种准确确定晶体塑性有限元模型材料参数的方法，以建立材料准确的晶体塑性有限元模型。
发明内容：
[0004]本发明的目的在于提供一种确定晶体塑性有限元模型材料参数的方法，该方法通过深度置信神经网络模型，采用较少的晶体塑性有限元计算量，获得准确的晶体塑性有限元的材料参数，解决了晶体塑性有限元材料参数试错法带来的不确定性、迭代次数多、计算量大、收敛性差及精度低的难题。
[0005]本发明解决上述难题的方案是：
[0006]步骤1：设置不同的材料参数组合A，进行晶体塑性有限元计算，得到不同参数下的真应力真应变；本步骤所述的晶体塑性有限元模型为唯象学幂率模型：
[0007]
[0008]式中：为第α个滑移系剪切应变率，是参考剪切应变率，τα为第α个滑移系的分解剪切应力，为临界剪切应力，m是滑移率速率敏感系数，为滑移系初始强度，h
αβ
为滑移硬化模量，表征滑移系β的单位剪切变形对滑移系α滑移阻力的影响，为第β个滑移系剪切
应变率，h
0为自硬化模量，τ
s
是滑移系饱和强度，τβ是滑移系β的分解剪切应力，a为应力指
数，q
αβ
是描述自硬化系数和潜在硬化系数的矩阵，可以表示为：
[0009]
[0010]其中，q是潜在硬化系数与自硬化系数的比值，对共面滑移系，比值为1，其他情况
和a为待确定的材料参数；
为1.4，A为三阶单位矩阵，参数m，h
αβ
[0011]步骤2：将参数组合A作为输入，把第二步得到的真应力真应变作为输出，建立深度置信神经网络模型；
[0012]步骤3：选取额外的材料参数组合B，进行晶体塑性有限元和深度置信神经网络模型计算，得到不同模型的真应力真应变，并进行比较，以验证深度置信神经网络模型的准确性；
[0013]步骤4：设置多组材料参数组合C，进行深度置信网络模型计算，将计算得到的真应力真应变和试验值比较，若误差大于设定值，继续修改材料参数组合C，直至误差小于设定值，该参数组合C即为晶体塑性有限元模型的材料参数。
附图说明：
[0014]图1深度置信神经网络模型预测精度验证；
[0015]图2实施例1获得的晶体塑性有限元模型预测精度；
[0016]表1材料参数组合A；
[0017]表2材料参数组合B；
[0018]表3晶体塑性有限元模型材料参数；
具体实施方式：
[0019]下面结合附图和具体实施案例对本发明进行详细说明。
[0020]本发明为一种可以快速获得晶体塑性有限元材料参数的方法，下面结合TC18钛合金的晶体塑性有
限元模型的建立过程为例，详细介绍本发明的具体实施细节，其方法包括：[0021]步骤1：设置不同的材料参数组合A，材料参数组合A如表1所示，进行晶体塑性有限元计算，得到不同参数下的真应力真应变；本步骤所述的晶体塑性有限元模型为唯象学幂率模型：
[0022]
[0023]式中：为第α个滑移系剪切应变率，是参考剪切应变率，τα为第α个滑移系的分
为滑移解剪切应力，为临界剪切应力，m是滑移率速率敏感系数，为滑移系初始强度，h
αβ
硬化模量，表征描述滑移系β的单位剪切变形对滑移系α滑移阻力的影响，为第β个滑移系
剪切应变率，h
0为自硬化模量，τ
s
是滑移系饱和强度，τβ是滑移系β的分解剪切应力，a为应力
指数，q
αβ
是描述自硬化系数和潜在硬化系数的矩阵，可以表示为：
[0024]
[0025]其中，q是潜在硬化系数与自硬化系数的比值，对共面滑移系，比值为1，其他情况
为1.4，A为三阶单位矩阵，参数m，h
αβ
和a为待确定的材料参数；
[0026]表1材料参数组合A
[0027]
[0028]步骤2：将参数组合A作为输入，把第二步得到的真应力真应变作为输出，建立深度置信神经网络模型；本实施例建立的深度置信神经网络模型具有三个隐含层，每一层的单元数量分别为10‑22‑19；
[0029]步骤3：选取额外的材料参数组合B，进行晶体塑性有限元和深度置信神经网络模型计算，得到不同模型的真应力真应变，并进行比较，以验证深度置信神经网络模型的准确性。表2为材料参数组合B；
[0030]表2材料参数组合B
[0031]
[0032]步骤4：设置多组材料参数组合C，进行深度置信网络模型计算，将计算得到的真应力真应变和试验值比较，若误差大于设定值，继续修改材料参数组合C，直至误差小于设定值，该参数组合C即为晶体塑性有限元模型的材料参数。
[0033]实施例建立的深度置信神经网络模型预测精度如图1所示，其中可以看到，深度置
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议