一种基于人工神经元网络的疲劳裂纹扩展速率预测方法
【技术领域】
[0001]本发明涉及人工神经元网络的应用领域，更具体地涉及一种基于人工神经元网络的疲劳裂纹扩展速率预测方法。
【背景技术】
[0002]在当今的航空工业领域，损伤容限理论得到广泛认可并应用于飞行器结构设计中。因此，基于线弹性断裂力学(LEFM)的疲劳裂纹扩展理论和方法被大量运用于材料结构的疲劳寿命的预测。在二十世纪60年代，Paris第一个将裂纹增长速率和应力强度因子幅值联系起来并成功地将LEFM理论应用于金属材料的疲劳裂纹扩展问题。
[0003]自从20世纪80年代以来，由于人工神经元网络对于人脑结构和机能的良好模拟，它开始被广泛的研究应用于各个行业。发展到现在，人工神经元网络已经具有部分与人脑相似的特征，尤其在信息的分布式存储，数据的并行处理和利用外来信息进行自学等方面表现优异。
[0004]由于在线弹性断裂力学下的疲劳裂纹扩展速率的预测可以看成是一个高度非线性的多输入单输出的系统，目前的简单公式对曲线的拟合和主要是对II区的拟合，其中大部分还是建立在将曲线视为对数坐标下线性的情况。然而，实际情况中的疲劳断裂曲线即使是在II区的对数坐标下，也不是完全线性的。而RBF神经网络不仅具有一般神经网络所具有的多维非线性拟合能力，更具有全局逼近能力和分离学习能力。通过RBF神经网络的这些特点，它能够在保证完成对疲劳裂纹扩展速率曲线模拟的前提下，拥有高于其他神经网络的精度和学习速度。
【发明内容】
[0005]本发明的目的是为了解决上述问题，提供一种基于人工神经元网络的疲劳裂纹扩展速率预测方法；
[0006]—种基于人工神经元网络的疲劳裂纹扩展速率预测方法，预测方法的实现平台为mat lab，包括以下几个步骤:
[0007]步骤一:获取载荷Kmax，应力比R，相对应载荷下的疲劳裂纹扩展速率da;
[0008]步骤二:对试件所受载荷与疲劳裂纹扩展速率进行预处理；
[0009]对载荷的预处理包括:若给定载荷形式为应力强度因子，则先将应力强度因子进行对数化，然后以其最大值最小值为限，使用指令mapminmax进行归一化，获得归一化参数psl;若给定载荷形式为应力，则通过裂纹长度与时间几何参数，计算该应力下裂纹尖端所受的应力强度因子，然后再进行归一化；
[0010]对疲劳裂纹扩展速率的预处理包括:将疲劳裂纹扩展速率进行对数化，然后以其最大值最小值为限，使用指令mapminmax进行归一化，获得归一化参数ps2;
[0011]步骤三:将预处理过后的载荷Kmax与应力比作为人工神经元网络的输入项，预处理之后的疲劳裂纹扩展速率作为人工神经元网络的输出项，使用matlab内函数newrbO对人工神经元网络进行训练；
[0012]步骤四:生成与载荷Kmax同等数量级的测试数据txl，与应力比相同数量级的测试数据tx2，输入人工神经元网络，获得人工神经元网络的拟合结果ty，将ty与疲劳裂纹扩展速率da进行对比，对其拟合精度进行检查；
[0013]步骤五:在拟合精度未达到要求时调整人工神经元网络的参数，重复步骤三、四，直至获得满足要求的人工神经元网络netl;
[0014]步骤六:获取需要预测的试件的各项几何参数以及试件所受的载荷F，应力比R2，载荷循环次数cycle;
[0015]步骤七:对获取的载荷进行预处理；
[0016]对获取载荷的预处理包括:若给定的载荷形式为试件所受应力，则先将应力强度因子进行对数化，然后以其最大值最小值为限，使用语句[Kmaxl，psl] =mapminmax(Kmaxl，0，1)进行归一化，获得归一化后的应力强度因子Kmaxl;若给定的载荷形式为试件受力，则通过矩形板的宽和厚，使用公式K = f(Geo，F)将其转化为应力形式，其中Geo为几何修正系数，f()为应力强度因子与试件受力的函数关系；
[0017]步骤八:将预处理过后的载荷Kmaxl与应力比R2输入，调用人工神经元网络netl，具体语句为tyl = sim(netl，[Kmaxl，R1 ])，获得神经元网络的输出结果tyl，并将所得结果进行后处理，获得此次加载下的疲劳裂纹扩展长度dal;
[0018]对神经元网络输出结果的后处理包括:先以步骤二中预处理所使用的上下限，使用语句mapminmax( ‘reverse’，tyl，ps2)对输出结果进行反归一化，再对其进行以e为底的指数化，获得本次疲劳裂纹的扩展长度；
[0019]步骤九:使用步骤八获得的疲劳裂纹扩展长度，对试件的几何参数进行更新；
[0020]将获得的裂纹增长长度累加到裂纹长度上，获得在新一轮载荷下的裂纹长度;并以此裂纹长度为基础，使用步骤七中的公式计算裂纹尖端处的应力强度因子；
[0021]步骤十:循环步骤七至九，直至疲劳裂纹长度达到预设值或者载荷循环次数cycle达到预设值为止；
[0022]步骤十一:得到本次计算所得的至每次加载为止的疲劳裂纹长度，即获得对试件在本加载下的疲劳裂纹扩展速率预测。
[0023]本发明的优点在于:
[0024](1)传统的线弹性锻炼力学将疲劳裂纹扩展速率与应力强度因子的之间的关系视作双对数坐标系下的线性关系，使用指数型模型对疲劳裂纹扩展速率进行拟合，并通过大量的实验数据对公式中的拟合参数进行校正。然而此方法忽视了其非线性特性，同时拟合参数的选择也需要建立在大量实验数据的基础上。
[0025]与传统的线弹性断裂力学下的疲劳裂纹扩展速率预测相比，基于人工神经元网络的疲劳裂纹扩展速率预测在人工神经元网络的基础上，拓展建立了一种全新的算法，对每一次加载下的疲劳裂纹扩展进行计算，之后以此为基础获得每一次加载时的疲劳裂纹长度，最后通过将预测数据与实验数据进行对比来验证预测的精度。除此之外，人工神经元网络的预测精度在训练阶段是可控的，可以通过调整均方误差目标，径向基函数的扩展速度，神经元的最大数目等使人工神经元网络的预测精度最大化。因此使用人工神经元网络建立实验数据之间的关系，在同样规模的实验数据的基础上，可以相对更精确的拟合出疲劳裂纹扩展速率与应力强度因子以及应力比之间的关系；
[0026](2)本发明对系统配置要求不高，方法全，具有可扩展性强，自动化，自学习程度高等特征。在具体操作上简单便捷，只要输入疲劳裂纹扩展速率与应力及应力比的实验数据，就能获得进行预测所需的人工神经元网络。再输入试件的受力，算法即可输出整个裂纹扩展的过程。若同时还输入了实验数据，算法还能一同输出实验数据与预测数据之间的对比；
[0027](3)本发明的拓展性很强，由于人工神经元网络在多元问题方面的优势，可以在预测的过程中加入更多与疲劳裂纹扩展有关的变量，如:裂纹尖端塑性区，裂纹闭合现象，裂纹尖端张开位移，残余应力等等。在加入这些变量之后，训练的出的神经元网络能够更好的拟合加载与裂纹扩展之间潜在的，未被发现的物理规律，更精确的预测疲劳裂纹扩展速率。
【附图说明】
[0028]图1:7075-T6铝合金在不同应力比下的裂纹扩展速率；
[0029]图2:人工神经元网络对实验数据的拟合结果图；
[0030]图3:人工神经元网络的预测结果与实验结果以及Forman公式的预测结果对比。
【具体实施方式】
[0031]下面将结合附图和实施例对本发明作进一步的详细说明。
[0032]本发明是一种基于人工神经元网络的疲劳裂纹扩展速率预测方法，预测方法的实现平台为matlab，包括以下几个步骤:
[0033]步骤一:载入通过具体实验数据得到的载荷Kmax，应力比R，相对应载荷下的疲劳裂纹扩展速率da。
[0034]标准的裂纹扩展速率数据为给定应力比R下裂纹扩展速率da/dn(每次加载时裂纹增长长度)与应力强度因子Kmax之间的关系；
[0035]步骤二:对试件所受载荷与疲劳裂纹扩展速率进行预处理。
[0036]具体为:
[0037]对载荷的预处理包括:若给定载荷形式为应力强度因子，则先将应力强度因子进行对数化，然后以其最大值最小值为限，使用指令mapminmax进行归一化，获得归一化参数pS1。若给定载荷形式为应力，则通过裂纹长度与时间几何参数，计算该应力下裂纹尖端所受的应力强度因子，然后再进行归一化。
[0038]对疲劳裂纹扩展速率的预处理包括:将疲劳裂纹扩展速率进行对数化，然后以其最大值最小值为限，使用指令mapminmax进行归一化，获得归一化参数ps2。
[0039]步骤三:将预处理过后的载荷Kmax与应力比作为人工神经元网络的输入项，预处理之后的疲劳裂纹扩展速率da/dn作为人工神经元网络的输出项，使用matlab内函数newrb()对人工神经元网络进行训练，其具体调用形式为:[net，tr] = newrb(P,T,GOAL，SPEAD，丽，DF)。(P为Q组输入量的R*Q维矩阵，T为Q组目标分类向量组成的S*Q维矩阵，GOAL为均方误差目标，SPREAD为径向基函数扩展速度，MN为神经元的最大数目，默认为Q，DF为两次显示之间所添加的神经元数目)
[0040]步骤四:生成与载荷Kmax同等数量级的测试数据txl，与应力比相同数量级的测试数据tx2，输入人工神经元网络，获得人工神经元网络的拟合结果ty，将ty与实验数据da进行对比，对其拟合精度进行检查。
[0041]用于测试的数据应在实验数据涵括的范围内，可以稍微超出。其拟合精度是否达到要求视具体情况而定。
[0042]步骤五:在拟合精度未达到要求时调整人工神经元网络的参数，重复步骤三、四，直至获得满意的人工神经元网络netl。
[0043]人工神经元网络的可调节参数包括均方误差目标GOAL，径向基函数的扩展速度SPEAD，神经元的最大数目MN等。通过调节神经与网络的这些参数，可以获得相对最好的拟合效果。
[0044]步骤六:获取需要预测的试件的各