CN114925584B - 变向锻造成形的金属晶粒尺寸和取向的智能协同调控方法 - Google Patents
变向锻造成形的金属晶粒尺寸和取向的智能协同调控方法 Download PDFInfo
- Publication number
- CN114925584B CN114925584B CN202210858281.9A CN202210858281A CN114925584B CN 114925584 B CN114925584 B CN 114925584B CN 202210858281 A CN202210858281 A CN 202210858281A CN 114925584 B CN114925584 B CN 114925584B
- Authority
- CN
- China
- Prior art keywords
- orientation
- metal
- grain size
- forging
- deformation
- 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.)
- Active
Links
Images
Classifications
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B21—MECHANICAL METAL-WORKING WITHOUT ESSENTIALLY REMOVING MATERIAL; PUNCHING METAL
- B21J—FORGING; HAMMERING; PRESSING METAL; RIVETING; FORGE FURNACES
- B21J9/00—Forging presses
- B21J9/10—Drives for forging presses
- B21J9/20—Control devices specially adapted to forging presses not restricted to one of the preceding subgroups
-
- 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
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- 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
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- 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
- 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
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Abstract
本发明公开了一种变向锻造成形的金属晶粒尺寸和取向的智能协同调控方法,属于金属智能制造领域,包括:以任一理想金属晶粒尺寸和取向作为目标,将当前时刻的变向锻造工艺参数输入神经网络模型,输出金属晶粒尺寸和取向的预测值,计算预测值与目标之间的距离;通过调整下一时刻的变向锻造工艺参数,计算前后时刻距离的倒数之差作为奖励,以从当前时刻到最终时刻累积奖励最大时对应的各时刻的变向锻造工艺参数作为最佳工艺路径,在最佳工艺路径下对金属进行变向锻造,最终得到理想金属晶粒尺寸和取向。本发明方法可以通过调整变向锻造工艺参数,对金属材料任意的晶粒尺寸和取向进行智能调控,不需要反复试错和实验,成本较低。
Description
技术领域
本发明属于金属智能制造领域,更具体地,涉及一种变向锻造成形的金属晶粒尺寸和取向的智能协同调控方法。
背景技术
制造业是国民经济的主体,也是“创新驱动,转型升级”的主战场,发展金属智能塑性成形工艺及方法,对于推动我国新时期的经济发展尤为重要。航空航天等关键零部件的加工对材料的服役性能有较高要求,根据性能导向定制化加工工艺,获得金属理想的微观组织特征具有重要应用前景。受加工路径、变形温度、应变速率等加工参数的影响,金属材料会经历复杂的微观组织演变过程,晶粒尺寸和取向特征均会发生改变。如何通过设计变形工艺及参数,调控材料在变形过程中的理想晶粒尺寸和取向特征,是当前研究的难点。
目前关于金属材料晶粒尺寸和取向调控的研究多是采用挤压、轧制或多向锻造等大变形工艺。申请号为202111155892.9,名称为“一种镁及镁合金无缝管材螺旋挤压装置及挤压工艺”的发明专利,其设计的装置通过在模具和螺旋内芯间形成坯料非均匀切变区,实现镁合金坯料在挤压过程中的螺旋流动,细化了晶粒尺寸并对镁及镁合金沿挤压方向的丝织构进行连续弱化以减小材料的拉压各向异性。该方法存在的问题是,只能单一的实现沿挤压方向的晶粒尺寸和取向改变的功能,而不能针对所需材料的性能,对材料任意的晶粒尺寸和取向特征进行智能调控。申请号为202111560971.8,名称为“一种大型饼类锻件的变向锻造工艺”的发明专利,其通过镦粗、横向压扁、压方倒角、滚圆后旋转镦粗四个步骤,实现改善锻件内部疏松,破碎夹杂物等缺陷,提高产品质量的目的。该方法存在的问题是,加工参数的选择依赖于人工经验和实验试错,成本较高,且无法实现对材料微观织构的有效调控。
由此可见,现有技术存在依赖于人工经验和实验试错、成本较高、难以对材料任意的晶粒尺寸和取向特征进行智能调控的技术问题。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种变向锻造成形的金属晶粒尺寸和取向的智能协同调控方法,由此解决通过设计工艺参数对金属进行塑性变形优化微观组织的现有技术存在依赖于人工经验和实验试错、成本较高、难以对材料任意的晶粒尺寸和取向特征进行智能调控的技术问题。
为实现上述目的,按照本发明的一个方面,提供了一种变向锻造成形的金属晶粒尺寸和取向的智能协同调控方法,包括:
以任一理想金属晶粒尺寸和取向作为目标,将当前时刻的变向锻造工艺参数输入神经网络模型,输出金属晶粒尺寸和取向的预测值,计算预测值与目标之间的距离;
通过调整下一时刻的变向锻造工艺参数,计算前后时刻距离的倒数之差作为奖励,以从当前时刻到最终时刻累积奖励最大时对应的各时刻的变向锻造工艺参数作为最佳工艺路径,在最佳工艺路径下对金属进行变向锻造,最终得到理想金属晶粒尺寸和取向;
所述神经网络模型通过如下方式训练得到:
通过仿真模拟获取金属在不同变向锻造工艺参数下,形成的晶粒尺寸和取向,将不同变向锻造工艺参数输入神经网络模型,输出各变向锻造工艺参数下的晶粒尺寸和取向,计算神经网络模型输出的晶粒尺寸和取向的加权和与仿真模拟形成的晶粒尺寸和取向的加权和之间的误差,反向传播更新神经网络模型的权重参数,训练至收敛,得到训练好的神经网络模型。
进一步地,所述变向锻造工艺参数包括变形方向和变形参数,所述神经网络模型的训练数据通过如下方式构造:
通过金属在不同变形方向和不同变形参数下进行压缩的力学响应和微观表征数据,拟合晶体塑性模型和微观组织状态方程;
对金属进行有限元模拟,提取不同变形方向和不同变形参数下锻造形成的金属心部有限单元上积分点的变形梯度,将变形梯度代入晶体塑性模型,计算金属各个晶粒的位错滑移、变形孪晶,得到应力应变分布和位错密度;
将应力应变分布和位错密度代入金属的微观组织状态方程,计算再结晶分数和晶粒尺寸,将位错密度减去用于再结晶的位错密度,得到新的位错密度,将新的位错密度带入晶体塑性模型,得到取向;
将不同变形方向、不同变形参数作为训练数据,将不同变形方向、不同变形参数下锻造金属,形成的晶粒尺寸和取向作为训练数据的模拟输出,用于与神经网络模型的输出进行对比。
进一步地,所述将新的位错密度带入晶体塑性模型,得到用于描述取向的欧拉角,所述取向通过如下方式描述:
其中,为对称广义谐波函数,是对称广义谐波函数的系数,为阶数,为晶体坐标系中的水平方向,为晶体坐标系中的垂直方向,M(l)和N(l)描述了晶体
对称性和金属对称性的线性独立谐波数,为晶体坐标系中的欧拉角,分别表示
章动角,进动角,自转角。
进一步地,所述变形方向包括:轧向、横向和厚度方向,所述变形参数包括:应变量、应变速率和变形温度,所述训练数据以及训练数据的模拟输出均进行归一化处理。
进一步地,所述晶粒尺寸和取向的加权和为:归一化后的晶粒尺寸×权重系数+归一化后的取向×(1-权重系数),其中,权重系数的取值范围是0.4-0.6。
进一步地,所述距离通过如下方式计算:
进一步地,所述奖励通过如下方式计算:
按照本发明的另一方面,提供了一种变向锻造成形的金属晶粒尺寸和取向的智能协同调控系统,包括:
训练模块,用于通过仿真模拟获取金属在不同变向锻造工艺参数下,形成的晶粒尺寸和取向,将不同变向锻造工艺参数输入神经网络模型,输出各变向锻造工艺参数下的晶粒尺寸和取向,计算神经网络模型输出的晶粒尺寸和取向的加权和与仿真模拟形成的晶粒尺寸和取向的加权和之间的误差,反向传播更新神经网络模型的权重参数,训练至收敛,得到训练好的神经网络模型;
调控模块,用于以任一理想金属晶粒尺寸和取向作为目标,将当前时刻的变向锻造工艺参数输入神经网络模型,输出金属晶粒尺寸和取向的预测值,计算预测值与目标之间的距离;通过调整下一时刻的变向锻造工艺参数,计算前后时刻距离的倒数之差作为奖励,以从当前时刻到最终时刻累积奖励最大时对应的各时刻的变向锻造工艺参数作为最佳工艺路径,在最佳工艺路径下对金属进行变向锻造,最终得到理想金属晶粒尺寸和取向。
进一步地,所述变向锻造工艺参数包括变形方向和变形参数,所述训练模块,用于通过金属在不同变形方向和不同变形参数下进行压缩的力学响应和微观表征数据,拟合晶体塑性模型和微观组织状态方程;对金属进行有限元模拟,提取不同变形方向和不同变形参数下锻造形成的金属心部有限单元上积分点的变形梯度,将变形梯度代入晶体塑性模型,计算金属各个晶粒的位错滑移、变形孪晶,得到应力应变分布和位错密度;将应力应变分布和位错密度代入金属的微观组织状态方程,计算再结晶分数和晶粒尺寸,将位错密度减去用于再结晶的位错密度,得到新的位错密度,将新的位错密度带入晶体塑性模型,得到取向;将不同变形方向、不同变形参数作为训练数据,将不同变形方向、不同变形参数下锻造金属,形成的晶粒尺寸和取向作为训练数据的模拟输出,用于与神经网络模型的输出进行对比。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,能够取得下列有益效果:
(1)本发明提出的一种变向锻造成形的金属晶粒尺寸和取向的智能协同调控方法,该方法可以通过调整变向锻造工艺参数,对金属材料任意的晶粒尺寸和取向进行智能调控,定制和优化材料的力学性能。引入神经网络,可以准确快速的进行预测,由此提高效率,降低成本。引入强化学习的方法,可以针对金属理想的晶粒尺寸和取向,智能优化变向锻造工艺参数,找到通向目标结构的最佳工艺路径,准确性高,可靠性强。
(2)本发明在获取训练数据时,通过宏微观跨尺度模拟方法,将宏观有限元模型与晶体塑性模型、材料微观组织状态方程进行迭代耦合,可以实现金属在变向锻造过程中宏观应力应变,微观再结晶分数、晶粒尺寸和取向特征的协同预测,由此获取的数据集进行训练,得到预测结果准确可靠的神经网络。
(3)晶体塑性模型中滑移系和孪晶系的schmid因子反映了多晶体中滑移系和孪晶系的激活程度,由此改变晶体取向的欧拉角。基于此,采用欧拉角作为对称广义谐波函数的输入,通过广义球谐(GSH)矢量对金属的取向分布函数进行逼近,来表示金属的取向特征,可以降低数据维度以提高训练效率。对训练数据进行归一化处理是为了提高模型精度和泛化能力。
(4)本发明可以在金属的不同方向进行变形,通过调整合适的工艺参数,包括变形方向、应变速率、应变量和变形温度,从而促进金属的晶体取向在变形过程中发生指定方向的旋转和指定滑移系、孪晶系的激活,并有效细化晶粒,实现理想晶粒尺寸和取向特征的调控。
(5)本发明设计了晶粒尺寸和取向的加权和表达,可以实现晶粒尺寸和取向的多目标协同优化。本发明在计算距离的时候计算的是初始时刻至某个时刻的累计距离,这样充分考虑了当前时刻和之前时刻的误差,有助于寻找最佳工艺路径。本发明设计的奖励函数,由于充分考虑了前后时刻的误差,得到一个更密集的学习奖励信息,有助于强化学习的准确寻优。
附图说明
图1是本发明实施例提供的方法流程图。
图2中(a)是本发明实施例提供的厚度方向ND锻造示意图。
图2中(b)是本发明实施例提供的轧向RD锻造示意图。
图2中(c)是本发明实施例提供的横向TD锻造示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
针对现有技术在对金属进行塑性变形优化微观组织时,工艺参数的设计多依赖经验和实验试错的方法,实验和人力成本较高,同时面临无法实现对金属任意的晶粒尺寸和取向进行智能化、协同化、精准化调控的技术瓶颈。本发明提供了一种变向锻造成形的金属晶粒尺寸和取向的智能协同调控方法,如图1所示,包括:
以任一理想金属晶粒尺寸和取向作为目标,将当前时刻的变向锻造工艺参数输入神经网络模型,输出金属晶粒尺寸和取向的预测值,计算预测值与目标之间的距离;
通过调整下一时刻的变向锻造工艺参数,计算前后时刻距离的倒数之差作为奖励,以从当前时刻到最终时刻累积奖励最大时对应的各时刻的变向锻造工艺参数作为最佳工艺路径,在最佳工艺路径下对金属进行变向锻造,最终得到理想金属晶粒尺寸和取向;
所述神经网络模型通过如下方式训练得到:
通过仿真模拟获取金属在不同变向锻造工艺参数下,形成的晶粒尺寸和取向,将不同变向锻造工艺参数输入神经网络模型,输出各变向锻造工艺参数下的晶粒尺寸和取向,计算神经网络模型输出的晶粒尺寸和取向的加权和与仿真模拟形成的晶粒尺寸和取向的加权和之间的误差,反向传播更新神经网络模型的权重参数,训练至收敛,得到训练好的神经网络模型。
训练过程的金属和智能协同调控时的金属是相同的,金属材料、金属初始织构、初始晶粒尺寸一致。
神经网络模型的训练包括:
变向锻造工艺参数包括变形方向和变形参数,变形方向(Di)包括:轧向(RD)、横向(TD)、厚度方向(ND)三个方向,变形参数具体指,应变速率,取值范围为:0.001~1s-1;变形温度,取值范围为:金属的再结晶温度~金属的熔点以下100℃;应变量,取值范围为:0~1;
通过金属在不同变形方向和不同变形参数下进行压缩的力学响应和微观表征,拟合晶体塑性模型和微观组织状态方程;
金属在不同变形方向压缩的力学响应具体指金属分别沿轧向(RD)、横向(TD)、厚度方向(ND)三个方向进行切取试样(使得试样满足基础实验要求),并进行单轴压缩实验,从而可以保证和变向锻造过程中材料在不同方向变形的真实力学响应保持一致。
晶体塑性模型可以为平均场也可以为全场形式的晶体塑性模型。
晶体塑性模型的表达形式有多种,本发明以其中一种为例进行介绍:
其中,F表示变形梯度,F e 为弹性应变梯度,F p 为塑性应变梯度,L p 为速度梯度,
为塑性应变梯度随时间的变化率,和分别为滑移系和孪晶系的塑性应变率,和分别为滑移系和孪晶系的schmid因子,为孪晶分数随时间的变化率,N S 和N tw 分别为
金属中滑移系a和孪晶系β的总数;
滑移系和孪晶系的schmid因子反映了多晶体中滑移系和孪晶的激活程度,会改变晶体取向的欧拉角。
拟合前设置初始模型参数以及模型参数的范围,拟合时,以晶体塑性模型的应力预测值与力学响应的应力之间的残差最小且晶体塑性模型的滑移系和孪晶系的schmid因子对应的欧拉角与微观表征所得欧拉角之间的误差最小为目标,对晶体塑性模型的待拟合的模型参数迭代寻优,得到最优模型参数。
微观组织状态方程:
其中,X为再结晶分数,为再结晶分数随时间的变化率,为位错密度,为
再结晶形核的临界位错密度,为塑性应变随时间的变化率(即晶体塑性模型计算的滑移
系和孪晶系的塑性应变率的加权和,通过damask软件计算),分别为微观组
织状态方程中再结晶模型待拟合的第一、第二、第三材料常数,
d drx 为再结晶晶粒尺寸,Z为热激活参数,与变形温度、应变速率相关,a 0和a 1分别为微观组织状态方程中再结晶晶粒尺寸模型待拟合的第一、第二材料常数。
其中,d为晶粒尺寸,d 0为初始晶粒尺寸,通过微观表征得到。
Z通过力学响应数据绘制的应力应变曲线确定,d通过微观表征得到,在已知输入输出的情况下通过最小二乘法拟合方程参数,
或者,拟合前设置初始方程参数和方程参数范围,拟合时,以微观组织状态方程计算的再结晶分数与微观表征的再结晶分数之间误差最小且微观组织状态方程计算的晶粒尺寸微观表征的晶粒尺寸之间的误差最小为目标,对微观组织状态方程的待拟合的方程参数迭代寻优,得到最优方程参数。
对金属进行有限元模拟,提取变向锻造工艺参数下锻造形成的金属心部有限单元上积分点的变形梯度,将变形梯度代入晶体塑性模型,计算金属各个晶粒的位错滑移、变形孪晶,得到应力应变分布和位错密度;
在商用有限元软件(如Abaqus)上进行变形锻造工艺的有限元数值模拟。金属的心部指中心位置,目的是金属在变向锻造过程中,应变分布不均匀,在中心位置处的应变与模拟设置的应变量是一致的。有限单元涉及到有限元法的基本原理,是其最小组成单元。有限元法的思路和原理是以结构力学中的位移法为基础,把复杂的结构或连续体看成有限个单元的组合,各单元在节点处彼此连续而组成整体,把连续体分成有限个单元和节点,称为离散化。积分点即为有限元单元上的节点。变形梯度是在有限元软件后处理中可以直接提取的状态变量。
将应力应变分布和位错密度代入金属的微观组织状态方程,计算再结晶分数和晶粒尺寸,将位错密度减去用于再结晶的位错密度,得到新的位错密度,将新的位错密度带入晶体塑性模型,得到用于描述取向的欧拉角;
将不同变形方向、不同变形参数作为训练数据,将不同变形方向、不同变形参数下锻造金属,形成的晶粒尺寸和取向作为训练数据的模拟输出,用于与神经网络模型的输出进行对比。
模拟出的金属织构采用广义球谐(GSH)矢量对金属的取向分布函数进行逼近,来表示金属的取向特征,可以降低数据维度以提高训练效率。
输入集和输出集均需进行归一化处理在-1到1之间。不同的锻造方向Di分别取值
为-1,0和1。输入集变量:应变量、应变速率和变形温度T,输出集变量:金属晶粒尺寸d
和取向特征t信息均采用如下函数进行归一化:
式中,x为上述输入和输出集变量:应变量、应变速率、变形温度、晶粒尺寸和取向特征。
搭建神经网络模型为多层前馈全连接神经网络模型(一个输入层,N个隐含层,一个输出层)或递归LSTM神经网络。
采用的DQN强化学习算法设置的奖励函数R为:
实施例1
实施例1的材料为轧制态6082铝合金,该合金晶体结构为FCC。初始晶粒尺寸为40μm,目标优化晶粒在20μm以下,并弱化轧向<111>的纤维织构。本发明提供的金属变向锻造成形的晶粒尺寸和取向智能协同调控方法的具体步骤如下:
步骤1、获取轧制态6082铝合金沿轧制方向,横向、厚度方向压缩的力学响应,应变速率设置为0.001s-1-0.1s-1,压缩变形量为60%,变形温度范围为350-500℃。并对压缩后的试样进行EBSD(Electron Backscattered Diffraction,电子背散射衍射)表征,获取晶粒尺寸和取向信息,同时利用压缩应力应变曲线和表征数据,对晶体塑性模型进行参数标定,从而可以实现沿不同变形方向不同变形参数下晶体取向的预测。
步骤2、构建数据集,基于Abaqus有限元平台,对6082铝合金进行不同变形参数交
互作用下的金属变向锻造工艺数值模拟,变形参数及取值范围如下:锻造方向(Di),可选为
轧向(RD)、横向(TD)、厚度方向(ND);应变量(),单道次的应变范围为:0-1;应变速率(
),取值范围为0.001-0.1 s-1,变形温度范围为350-500℃。不同锻造变形方向如图2中(a)、
图2中(b)和图2中(c)所示,Forging direction表示锻造方向,在模拟时主要进行多道次不
同变形方向及路径耦合的影响,共计64组模拟。并提取模拟变形后合金心部代表单元的应
变梯度。晶体塑性模型根据输入的变形梯度,计算金属各个晶粒的位错滑移、变形孪晶,得
到应力应变分布和位错密度等状态变量。同时,将晶体塑性模型得到的状态变量传递给材
料的微观组织状态方程用于计算动态再结晶分数、晶粒尺寸,并传回晶体塑性模型中对各
个晶粒的状态变量进行更新,获取该变形条件下6082铝合金的晶粒尺寸d和取向特征信息t。建立以不同工艺参数为输入集,6082铝合金的晶粒尺寸和取向特征信息为输出集的对应
关系。
步骤3、利用步骤2中获取的数据集,搭建三层全连接神经网络模型,分别包含一个输入层,一个隐含层和一个输出层。训练前对数据集进行归一化处理至-1到1之间。隐含层间选择Relu作为模型的激活函数,采用随机梯度下降法(Adam)来寻找模型参数的最优解,训练得到隐含层节点数为10。并使用该模型预测任意变形条件下6082铝合金的晶粒尺寸d和取向特征信息t。
达式为:。式中,为归一
化后的晶粒尺寸和取向特征,n为权重系数,取值为0.5。采用基于DQN的强化学习算法,针对
细化晶粒尺寸在20μm以下,并弱化轧向<111>的纤维织构为目标, 目标动作值为:,其中,为优化空间域,。其中距离函数
,具体表示为:
奖励函数R为:
式中,T为当前时刻,折扣因子初始取值为0.96。通过在神经网络中进行迭代循环,自动学习出获取输出得到理想6082铝合金晶粒尺寸和取向特征协同调控目标值的最优工艺参数,并保存模型。
步骤5、输出细化6082铝合金晶粒尺寸并弱化沿轧向织构的变向锻造完整工艺路径和参数设置,实现对6082铝合金晶粒尺寸和取向的智能精准调控。
实施例2
实施例2的材料为轧制TA15钛合金,该合金包含α相和β相。初始晶粒尺寸为30μm,目标优化晶粒在15μm以下,并提高其基面织构。本发明提供的金属变向锻造成形的晶粒尺寸和取向智能协同调控方法的具体步骤如下:
步骤1、获取轧制态TA15钛合金沿轧制方向,横向、厚度方向压缩的力学响应,应变速率设置为0.001s-1-0.1s-1,压缩变形量为60%,变形温度范围为800-950℃。并对压缩后的试样进行EBSD表征,获取两相晶粒尺寸和取向信息,同时利用压缩应力应变曲线和表征数据,对晶体塑性模型进行参数标定,从而可以实现沿不同变形方向不同变形参数下晶体取向的预测。
步骤2、构建数据集,基于Abaqus有限元平台,对TA15钛合金进行不同变形参数交
互作用下的金属变向锻造工艺数值模拟,变形参数及取值范围如下:锻造方向(Di),可选为
轧向(RD)、横向(TD)、厚度方向(ND);应变量(),单道次的应变范围为:0-1;应变速率(
),取值范围为0.001-0.1 s-1,变形温度范围为800-950℃。不同锻造变形方向如图2中(a)、
图2中(b)和图2中(c)所示,在模拟时主要进行多道次不同变形方向及路径耦合的影响,共
计64组模拟。并提取模拟变形后合金心部代表单元的应变梯度。晶体塑性模型根据输入的
变形梯度,计算金属各个晶粒的位错滑移、变形孪晶,得到应力应变分布和位错密度等状态
变量。同时,将晶体塑性模型得到的应力应变分布和状态参量传递给材料的微观组织状态
方程用于计算动态再结晶分数、晶粒尺寸以及再结晶,并传回晶体塑性模型中对各个晶粒
的状态变量进行更新,获取该变形条件下轧制态TA15钛合金的晶粒尺寸d和取向特征信息t。建立以不同工艺参数为输入集,轧制态TA15钛合金的晶粒尺寸和取向特征信息为输出集
的对应关系。
步骤3、利用步骤2中获取的数据集,搭建三层全连接神经网络模型,分别包含一个输入层,一个隐含层和一个输出层。训练前对数据集进行归一化处理至-1到1之间。隐含层间选择Relu作为模型的激活函数,采用随机梯度下降法(Adam)来寻找模型参数的最优解,训练得到隐含层节点数为10。并使用该模型预测任意变形条件下轧制态TA15钛合金的两相晶粒尺寸d和取向特征信息t。
步骤4、设计轧制态TA15钛合金晶粒尺寸和取向特征协同调控的多目标函数,函数
表达式为:。式中,为归一化后的晶粒尺寸和取
向特征,n为权重系数,取值为0.4。采用基于DQN的强化学习算法,针对细化晶粒尺寸在15μm
以下,并提高其基面织构为目标, 目标动作值为:,其中, 为
优化空间域,。其中距离函数,具体表示为:
奖励函数R为:
式中,T为当前时刻,折扣因子初始取值为0.92。通过在神经网络中进行迭代循环,自动学习出获取输出得到理想金属晶粒尺寸和取向特征协同调控目标值的最优工艺参数,并保存模型。
步骤5、输出细化TA15钛合金两相晶粒尺寸并提高基面织构的变向锻造完整工艺路径和参数设置,实现对TA15钛合金晶粒尺寸和取向的智能精准调控。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种变向锻造成形的金属晶粒尺寸和取向的智能协同调控方法,其特征在于,包括:
以任一理想金属晶粒尺寸和取向作为目标,将当前时刻的变向锻造工艺参数输入神经网络模型,输出金属晶粒尺寸和取向的预测值,计算预测值与目标之间的距离;
通过调整下一时刻的变向锻造工艺参数,计算前后时刻距离的倒数之差作为奖励,以从当前时刻到最终时刻累积奖励最大时对应的各时刻的变向锻造工艺参数作为最佳工艺路径,在最佳工艺路径下对金属进行变向锻造,最终得到理想金属晶粒尺寸和取向;
所述神经网络模型通过如下方式训练得到:
通过仿真模拟获取金属在不同变向锻造工艺参数下,形成的晶粒尺寸和取向,将不同变向锻造工艺参数输入神经网络模型,输出各变向锻造工艺参数下的晶粒尺寸和取向,计算神经网络模型输出的晶粒尺寸和取向的加权和与仿真模拟形成的晶粒尺寸和取向的加权和之间的误差,反向传播更新神经网络模型的权重参数,训练至收敛,得到训练好的神经网络模型。
2.如权利要求1所述的一种变向锻造成形的金属晶粒尺寸和取向的智能协同调控方法,其特征在于,所述变向锻造工艺参数包括变形方向和变形参数,所述神经网络模型的训练数据通过如下方式构造:
通过金属在不同变形方向和不同变形参数下进行压缩的力学响应和微观表征数据,拟合晶体塑性模型和微观组织状态方程;
对金属进行有限元模拟,提取不同变形方向和不同变形参数下锻造形成的金属心部有限单元上积分点的变形梯度,将变形梯度代入晶体塑性模型,计算金属各个晶粒的位错滑移、变形孪晶,得到应力应变分布和位错密度;
将应力应变分布和位错密度代入金属的微观组织状态方程,计算再结晶分数和晶粒尺寸,将位错密度减去用于再结晶的位错密度,得到新的位错密度,将新的位错密度带入晶体塑性模型,得到取向;
将不同变形方向、不同变形参数作为训练数据,将不同变形方向、不同变形参数下锻造金属,形成的晶粒尺寸和取向作为训练数据的模拟输出,用于与神经网络模型的输出进行对比。
4.如权利要求2所述的一种变向锻造成形的金属晶粒尺寸和取向的智能协同调控方法,其特征在于,所述变形方向包括:轧向、横向和厚度方向,所述变形参数包括:应变量、应变速率和变形温度,所述训练数据以及训练数据的模拟输出均进行归一化处理。
5.如权利要求4所述的一种变向锻造成形的金属晶粒尺寸和取向的智能协同调控方法,其特征在于,所述应变量的取值范围是0~1,所述应变速率的取值范围是0.001~1s-1,所述变形温度的取值范围是金属的再结晶温度~金属的熔点以下100℃。
6.如权利要求4所述的一种变向锻造成形的金属晶粒尺寸和取向的智能协同调控方法,其特征在于,所述晶粒尺寸和取向的加权和为:归一化后的晶粒尺寸×权重系数+归一化后的取向×(1-权重系数),其中,权重系数的取值范围是0.4-0.6。
9.一种变向锻造成形的金属晶粒尺寸和取向的智能协同调控系统,其特征在于,包括:
训练模块,用于通过仿真模拟获取金属在不同变向锻造工艺参数下,形成的晶粒尺寸和取向,将不同变向锻造工艺参数输入神经网络模型,输出各变向锻造工艺参数下的晶粒尺寸和取向,计算神经网络模型输出的晶粒尺寸和取向的加权和与仿真模拟形成的晶粒尺寸和取向的加权和之间的误差,反向传播更新神经网络模型的权重参数,训练至收敛,得到训练好的神经网络模型;
调控模块,用于以任一理想金属晶粒尺寸和取向作为目标,将当前时刻的变向锻造工艺参数输入神经网络模型,输出金属晶粒尺寸和取向的预测值,计算预测值与目标之间的距离;通过调整下一时刻的变向锻造工艺参数,计算前后时刻距离的倒数之差作为奖励,以从当前时刻到最终时刻累积奖励最大时对应的各时刻的变向锻造工艺参数作为最佳工艺路径,在最佳工艺路径下对金属进行变向锻造,最终得到理想金属晶粒尺寸和取向。
10.如权利要求9所述的一种变向锻造成形的金属晶粒尺寸和取向的智能协同调控系统,其特征在于,所述变向锻造工艺参数包括变形方向和变形参数,所述训练模块,用于通过金属在不同变形方向和不同变形参数下进行压缩的力学响应和微观表征数据,拟合晶体塑性模型和微观组织状态方程;对金属进行有限元模拟,提取不同变形方向和不同变形参数下锻造形成的金属心部有限单元上积分点的变形梯度,将变形梯度代入晶体塑性模型,计算金属各个晶粒的位错滑移、变形孪晶,得到应力应变分布和位错密度;将应力应变分布和位错密度代入金属的微观组织状态方程,计算再结晶分数和晶粒尺寸,将位错密度减去用于再结晶的位错密度,得到新的位错密度,将新的位错密度带入晶体塑性模型,得到取向;将不同变形方向、不同变形参数作为训练数据,将不同变形方向、不同变形参数下锻造金属,形成的晶粒尺寸和取向作为训练数据的模拟输出,用于与神经网络模型的输出进行对比。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210858281.9A CN114925584B (zh) | 2022-07-21 | 2022-07-21 | 变向锻造成形的金属晶粒尺寸和取向的智能协同调控方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210858281.9A CN114925584B (zh) | 2022-07-21 | 2022-07-21 | 变向锻造成形的金属晶粒尺寸和取向的智能协同调控方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114925584A CN114925584A (zh) | 2022-08-19 |
CN114925584B true CN114925584B (zh) | 2022-09-20 |
Family
ID=82815882
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210858281.9A Active CN114925584B (zh) | 2022-07-21 | 2022-07-21 | 变向锻造成形的金属晶粒尺寸和取向的智能协同调控方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114925584B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4207762A (en) * | 1978-08-11 | 1980-06-17 | Kelly Vernon A Jr | Method of forming high quality forgings |
CN106055844A (zh) * | 2016-07-06 | 2016-10-26 | 中南大学 | 基于bp神经网络的镍基高温合金微观组织预测控制方法 |
CN106777594A (zh) * | 2016-12-01 | 2017-05-31 | 武汉科技大学 | 一种基于连续曲面的轧制过程数学模型的自学习方法 |
EP3674984A1 (en) * | 2018-12-29 | 2020-07-01 | Dassault Systèmes | Set of neural networks |
CN114417664A (zh) * | 2022-01-04 | 2022-04-29 | 大连理工大学 | 一种基于元胞自动机的钢材热轧微观组织演化在线模拟和可视化方法 |
CN114764525A (zh) * | 2021-01-13 | 2022-07-19 | 上海交通大学 | 针对车轴多道次多工步锻造成形过程的晶粒尺寸预测方法 |
-
2022
- 2022-07-21 CN CN202210858281.9A patent/CN114925584B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4207762A (en) * | 1978-08-11 | 1980-06-17 | Kelly Vernon A Jr | Method of forming high quality forgings |
CN106055844A (zh) * | 2016-07-06 | 2016-10-26 | 中南大学 | 基于bp神经网络的镍基高温合金微观组织预测控制方法 |
CN106777594A (zh) * | 2016-12-01 | 2017-05-31 | 武汉科技大学 | 一种基于连续曲面的轧制过程数学模型的自学习方法 |
EP3674984A1 (en) * | 2018-12-29 | 2020-07-01 | Dassault Systèmes | Set of neural networks |
CN114764525A (zh) * | 2021-01-13 | 2022-07-19 | 上海交通大学 | 针对车轴多道次多工步锻造成形过程的晶粒尺寸预测方法 |
CN114417664A (zh) * | 2022-01-04 | 2022-04-29 | 大连理工大学 | 一种基于元胞自动机的钢材热轧微观组织演化在线模拟和可视化方法 |
Non-Patent Citations (7)
Title |
---|
2024Al/Al18B4O33w复合材料高温变形行为研究;柏阳;《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》;20220215(第02期);B020-241 * |
Behavior and modeling of microstructure evolution during metadynamic recrystallization of a Ni-based superalloy;XuefengTang,BaoyuWang;《Materials Science and Engineering》;20161015;第675卷;192-203 * |
Optimization of Open Die Ironing Process through Artificial Neural Network for Rapid Process Simulation;Silvia Mancini,Luigi Langellotto,etc;《Metals》;20201021;第10卷(第10期);1397 * |
Successful Pass Schedule Design in Open-Die Forging Using Double Deep Q-Learning;Niklas Reinisch,Fridtjof Rudolph,etc;《Processes》;20210622;第9卷(第7期);1084 * |
基于BP神经网络的2.25Cr1Mo0.25V钢高温流变应力模型研究;刘俊林;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20140915(第09期);B022-27 * |
面向大型航空构件形/性控制的局部控流和模具控温模锻工艺研究;姜静;《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》;20200115(第01期);B022-414 * |
齿轮坯预锻件的多目标优化设计;张明月;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20120715(第07期);B022-266 * |
Also Published As
Publication number | Publication date |
---|---|
CN114925584A (zh) | 2022-08-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106096134B (zh) | 基于损伤力学的金属结构疲劳可靠性分析与优化设计方法 | |
Derazkola et al. | Review on modeling and simulation of dynamic recrystallization of martensitic stainless steels during bulk hot deformation | |
Quan et al. | Artificial neural network modeling to evaluate the dynamic flow stress of 7050 aluminum alloy | |
CN113128124B (zh) | 基于改进神经网络的多牌号C-Mn钢力学性能预测方法 | |
CN113156900B (zh) | 一种基于元强化学习的加工变形控制方法 | |
Feng et al. | Optimization of variable blank holder force in deep drawing based on support vector regression model and trust region | |
Shao et al. | A new approach of preform design for forging of 3D blade based on evolutionary structural optimization | |
Heng et al. | Towards intelligent design optimization: Progress and challenge of design optimization theories and technologies for plastic forming | |
Yan et al. | Multi-objective quality control method for cold-rolled products oriented to customized requirements | |
Chen et al. | Mesoscale modeling of dynamic recrystallization: multilevel cellular automaton simulation framework | |
CN116502455A (zh) | 一种激光选区熔化技术的工艺参数确定方法及系统 | |
Kodippili et al. | Data-driven prediction of forging outcome: Effect of preform shape on plastic strain in a magnesium alloy forging | |
CN114925584B (zh) | 变向锻造成形的金属晶粒尺寸和取向的智能协同调控方法 | |
Liu et al. | Hot deformation behavior of near-β titanium alloy Ti-3Mo-6Cr-3Al-3Sn based on phenomenological constitutive model and machine learning algorithm | |
CN110991605A (zh) | 多变量时间序列深度信念网络的低压铸造模温预测方法 | |
CN114564787A (zh) | 用于目标相关翼型设计的贝叶斯优化方法、装置及存储介质 | |
Quan et al. | Modelling the hot flow behaviors of AZ80 alloy by BP-ANN and the applications in accuracy improvement of computations | |
Egorova et al. | Forecasting VT6 titanium alloy rolled bar mechanical properties | |
Mandal et al. | Modeling microstructural evolution during dynamic recrystallization of alloy D9 using artificial neural network | |
Zhang et al. | Rolling force prediction in heavy plate rolling based on uniform differential neural network | |
CN114384931B (zh) | 一种基于策略梯度的无人机多目标最优控制方法和设备 | |
Dong et al. | Just-in-time learning-based soft sensor for mechanical properties of strip steel via multi-block weighted semisupervised models | |
Wei et al. | High-cycle fatigue SN curve prediction of steels based on a transfer learning-guided convolutional neural network | |
Li et al. | A preliminary discussion about the application of machine learning in the field of constitutive modeling focusing on alloys | |
de Melo Filho et al. | Design space exploration for resonant metamaterials using physics guided neural networks |
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 |