《一种动态再结晶模型参数反分析识别方法.pdf》由会员分享，可在线阅读，更多相关《一种动态再结晶模型参数反分析识别方法.pdf（8页珍藏版）》请在专利查询网上搜索。
1、10申请公布号CN102339344A43申请公布日20120201CN102339344ACN102339344A21申请号201110143339322申请日20110525G06F17/50200601G06N3/0220060171申请人深圳大学地址518060广东省深圳市南山区南海大道3688号72发明人娄燕李落星彭怀伟吴文华74专利代理机构深圳市中知专利商标代理有限公司44101代理人张皋翔54发明名称一种动态再结晶模型参数反分析识别方法57摘要本发明公开了一种动态再结晶模型参数反分析识别方法，基于遗传神经网络元胞自动机，其特征在于依次有以下步骤1确定参数的取值范围，设计构造实验方。
2、案；2进行模拟计算，构建学习样本库；3构建BP神经网络；4训练BP神经网络，建立非线性映射关系；5随机选取至少一组参数，预测流变应力；6比较预测值与实测值，直至误差达到设计要求；7保存误差达到设计要求的参数最佳值。有效克服了对动态再结晶参数取值的随意性，用正交试验、遗传算法、神经网络和数值模拟优化计算出动态再结晶参数的取值，解决了反分析优化计算时不稳定性、强非线性、计算量大及学习样本不均匀等问题，明显提高再结晶参数取值的准确性和快速性。51INTCL19中华人民共和国国家知识产权局12发明专利申请权利要求书1页说明书5页附图1页CN102339357A1/1页21一种动态再结晶模型参数反分析识。
3、别方法，基于遗传神经网络元胞自动机，其特征在于依次有以下步骤1确定参数的取值范围，设计构造实验方案依据实际问题，确定材料动态再结晶元胞自动机模型参数的取值范围，并依据正交试验设计构造实验方案；2进行模拟计算，构建学习样本库采用有限元模拟对构造的每一组实验方案进行模拟计算，获得每一组实验方案关键点的流变应力值，并将每一组实验方案与对应的流变应力值构成一组学习样本，构建学习样本库；3构建BP神经网络构建BP神经网络结构；4训练BP神经网络，建立非线性映射关系用所述学习样本对BP神经网络进行训练，建立待反分析的动态再结晶参数与流变应力之间的非线性映射关系；5随机选取至少一组参数，预测流变应力随机选取。
4、至少一组动态再结晶参数，用所述训练后的BP神经网络预测流变应力；6比较预测值与实测值，直至误差达到设计要求将预测的所述流变应力值与对应的实测流变应力值进行比较，计算误差，如果误差达不到设计要求，则用遗传算法进行优化计算，推出新的一代群体，所述新的一代群体是一组新的、优化的动态再结晶参数，将所述新的动态再结晶参数用所述训练后的神经网络预测流变应力，再将预测的所述流变应力值与实测的流变应力值进行比较，计算误差，直至误差达到设计要求；7保存误差达到设计要求的参数最佳值保存误差达到设计要求的动态再结晶参数最佳值。2如权利要求1所述的动态再结晶模型参数反分析识别方法，其特征在于所述步骤3的BP神经网络结。
5、构是对神经网络的隐含层层数和每个隐含层节点数从网络精度和训练时间上综合考虑建立的。3如权利要求1或2所述的动态再结晶模型参数反分析识别方法，其特征在于所述步骤6计算误差的公式为式中BP为关键点的神经网络预测应力值；0为对应的实测应力值。4如权利要求3所述的动态再结晶模型参数反分析识别方法，其特征在于所述步骤6的误差如果达不到设计要求，则将误差E设置为遗传算法的个体适应度值，将所述随机选取至少一组动态再结晶参数设置为初始群体，然后用遗传算法的选择、交叉、变异推出新的一代群体。权利要求书CN102339344ACN102339357A1/5页3一种动态再结晶模型参数反分析识别方法技术领域0001本。
6、发明涉及材料动态再结晶微观模拟的参数反分析识别技术，特别是涉及一种基于遗传神经网络元胞自动机的动态再结晶模型参数反分析识别方法。背景技术0002动态再结晶行为是改善材料热变形微观组织及性能的重要机制。定量研究动态再结晶过程的物理冶金原理，形成预报动态再结晶过程的数值方法，对控制材料组织演变和获得良好力学性能具有重要意义。然而，现有用元胞自动机模拟动态再结晶过程，是随意选取元胞自动机模型参数，只能通过反复调用有限元模拟，寻找计算流变应力值和最接近的实测流变应力值，以获得较理想的动态再结晶数值模拟过程。虽然用有限元模拟代替实际试验过程，节约了成本，但模拟时的网格重划降低了计算效率；有些用数值模拟和。
7、高斯牛顿迭代算法对材料参数进行反分析识别，但此回归拟合方法只能解决简单的线性和非线性映射关系，不能解决动态再结晶的形核参数与宏观流变应力之间复杂的非线性关系；有些用反向传播BACKPROPAGATION，缩略词为BP神经网络对铸件与铸型间界面热换系数进行反分析识别，能解决复杂的非线性关系，在调整学习样本数据时采用试凑法，却无法快速优化计算出所需参数值，计算量大，网络收敛速度慢；还有些采用自适应响应面方法反求低碳钢动态再结晶形核参数的全局最优值，虽然计算量比较小，速度也比较快，但预测精度差，样本数据容易失衡。发明内容0003本发明所要解决的技术问题是弥补上述现有技术的缺陷，提供一种基于遗传神经网。
8、络元胞自动机的动态再结晶模型参数反分析识别方法。0004本发明基于正交试验、遗传算法、BP神经网络和有限元数值模拟技术优化计算材料动态再结晶元胞自动机模型参数，可以有效解决动态再结晶模型参数反分析求解过程的不稳定性、强非线性、计算量大及学习样本不平衡等问题，实现动态再结晶模型参数的快速、准确求解。0005本发明的技术问题通过以下技术方案予以解决。0006这种动态再结晶模型参数反分析识别方法，基于遗传神经网络元胞自动机。0007这种动态再结晶模型参数反分析识别方法的特点是0008依次有以下步骤00091确定参数的取值范围，设计构造实验方案0010依据实际问题，确定材料动态再结晶元胞自动机模型参数。
9、的取值范围，并依据正交试验设计构造实验方案；00112进行模拟计算，构建学习样本库0012采用有限元模拟对构造的每一组实验方案进行模拟计算，获得每一组实验方案关键点的流变应力值，并将每一组实验方案与对应的流变应力值构成一组学习样本，构建学说明书CN102339344ACN102339357A2/5页4习样本库；00133构建BP神经网络0014构建BP神经网络结构；00154训练BP神经网络，建立非线性映射关系0016用所述学习样本对BP神经网络进行训练，建立待反分析的动态再结晶参数与流变应力之间的非线性映射关系；00175随机选取至少一组参数，预测流变应力0018随机选取至少一组动态再结晶参。
10、数，用所述训练后的BP神经网络预测流变应力；00196比较预测值与实测值，直至误差达到设计要求0020将预测的所述流变应力值与对应的实测流变应力值进行比较，计算误差，如果误差达不到设计要求，则用遗传算法进行优化计算，推出新的一代群体，所述新的一代群体是一组新的、优化的动态再结晶参数，将所述新的动态再结晶参数用所述训练后的神经网络预测流变应力，再将预测的所述流变应力值与实测的流变应力值进行比较，计算误差，直至误差达到设计要求；00217保存误差达到设计要求的参数最佳值0022保存误差达到设计要求的动态再结晶参数最佳值。0023本发明的技术问题通过以下进一步的技术方案予以解决。0024所述步骤1构。
11、造的实验方案中的数据样本分布平衡、均匀。0025所述步骤3的BP神经网络结构是对神经网络的隐含层层数和每个隐含层节点数从网络精度和训练时间上综合考虑建立的。0026所述步骤6计算误差的公式为0027式中0028BP为关键点的神经网络预测应力值；00290为对应的实测应力值。0030所述步骤6的误差达不到设计要求，则将误差E设置为遗传算法的个体适应度值，将所述随机选取几组动态再结晶参数设置为初始群体，然后用遗传算法的选择、交叉、变异推出新的一代群体。0031本发明与现有技术对比的有益效果是0032本发明方法有效克服了对动态再结晶参数取值的随意性，用正交试验、遗传算法、神经网络和数值模拟优化计算出。
12、动态再结晶参数的取值，解决了反分析优化计算时不稳定性、强非线性、计算量大及学习样本不均匀等问题，明显提高再结晶参数取值的准确性和快速性。附图说明0033图1是本发明具体实施方式的微观组织模拟金相图；0034图2是本发明具体实施方式的微观组织实测金相图。具体实施方式说明书CN102339344ACN102339357A3/5页50035下面结合具体实施方式并对照附图对本发明进行说明。0036一种如图12所示的基于遗传神经网络元胞自动机的动态再结晶模型参数反分析识别方法，对AZ31镁合金试样以下简称样品圆柱热压缩变形时具体的动态再结晶模型参数进行反分析识别，步骤如下00371确定参数的取值范围，设。
13、计构造实验方案0038依据实际问题，确定试样材料动态再结晶元胞自动机模型参数的取值范围，其模型参数有4个，取值范围为0039初始位错密度CR128091014，149761014M2；0040形核率常数NC5000，250000M1S1；0041初始软化参数R06000，8900；0042初始硬化参数H0100001015，129001015；0043依据正交试验表设计构造三十二组实验方案，其中三组实验方案见下表10044表1004500462进行模拟计算，构建学习样本库0047采用有限元模拟对构造的每一组实验方案进行模拟计算，模拟时试样材料参数取值见下表20048表200490050模拟时压。
14、缩温度为300，应变速率为03S1，材料与模具的摩擦系数取1。0051模拟结束后，获得每一组实验方案流变应力峰值BP，并将每组实验方案与对应的流变应力峰值BP构成一组学习样本，即输入值为4个，输出值为1个，从而构建学习样本库。三组实验方案对应的学习样本见下表30052表30053说明书CN102339344ACN102339357A4/5页600543构建BP神经网络0055通过反复计算，从网络精度和训练时间上综合考虑，BP神经网络采用双隐层结构，具体各层节点数为435371。此时网络均方误差00016，训练时间64367S。00564训练BP神经网络，建立非线性映射关系0057用所述学习样本。
15、对BP神经网络对进行训练，即采用435371结构的BP神经网络建立待反分析的动态再结晶参数与流变应力之间的非线性映射关系。00585随机选取十组参数，预测流变应力0059随机选取十组动态再结晶参数，用所述训练后的BP神经网络进行预测，得到相应十个预测应力峰值，其中的三组动态再结晶参数和预测应力峰值见下表40060表4006100626比较预测值与实测值，直至误差达到设计要求0063将预测的流变应力峰值与实测的流变应力峰值进行比较，计算误差，误差公式如下0064式中0065BP为神经网络预测应力峰值；00660为实测应力峰值。0067实测流变应力试验是在GLEEBLE3500热模拟试验机上进行等。
16、应变速率下热压缩实验，压缩前，将试样原挤压棒材加工成10MM12MM的圆柱体试样，石墨薄片作为润滑剂放在样品和压头间，以减小试样与压头间的摩擦影响，样品在变形过程中的温度变化用直径为008MM的热电偶记录，测量反应时间为01S。样品加热到预设温度的速率为10/S，保温时间为60S，以保证温度的均匀化，然后由高度12MM压缩到高度44MM，即达到真应变1，整个试验都在氮气氛下进行，压缩温度为300，应变速率为03S1，实测应力峰值为说明书CN102339344ACN102339357A5/5页7906786MPA，则对应的三组参数误差见下表50068表500690070计算出的十组误差都达不到设。
17、计要求E0001，则用遗传算法的选择、交叉、变异推出新的一代群体，遗传算法参数如下变量个数4，种群大小10，最大迭代次数100，交叉概率04，变异概率02，采用实数编码；将误差E设置为遗传算法的个体适应度值，将随机选取十组动态再结晶参数编码为初始群体，然后用遗传算法的选择、交叉、变异推出新的一代群体，将其解码形成新的动态再结晶参数，即CR143201014M2；NC170000M1S1；R072532；H01231015。将这些新的动态再结晶参数用训练后的神经网络预测其流变应力峰值为858760MPA，再与实测的流变应力峰值进行比较，计算出的误差E005300001仍达不到设计要求，则将此组新。
18、的动态再结晶参数与随机选取的十组动态再结晶参数中误差小的九组参数，再组成十组动态再结晶参数，仍然编码为初始群体，然后用遗传算法的选择、交叉、变异优化，推出新的又一代群体，直至误差达到设计要求E0001。00717保存误差达到设计要求E0001的动态再结晶参数值CR4524921014M2，NC129300M1S1，0511297，H01681015。0072上述动态再结晶参数值是压缩温度为300、应变速率为03S1时的动态再结晶参数值。0073图1是压缩温度为300、应变速率03S1时的试样热变形的微观组织模拟的金相图，图2是与之对应的微观组织实测的金相图。对比表明试样热变形时，在晶界产生动态再结晶“项链”结构，晶内大量出现了细小再结晶晶粒。用DM1500金相仪测量图1和图2的晶粒大小，分别为5776M和5782M，非常接近，达到材料动态再结晶精确数值模拟的要求。0074以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下做出如果干等同替代或明显变型，而且性能或用途相同，都应当视为属于本发明由所提交的权利要求书确定的专利保护范围。说明书CN102339344ACN102339357A1/1页8图1图2说明书附图CN102339344A。