含应力双相铁磁材料中磁巴克豪森噪声信号的模拟方法，属于磁性物理学与无损检测领域的交叉内容，既有助于揭示物理现象，又可以作为双相材料组分占比和应力无损检测的理论依据。
背景技术：
磁巴克豪森噪声检测方法具有无损检测的特点，对材料内部组分占比和应力的变化非常敏感。实验过程中观察到双相铁磁性材料中检测得到的磁巴克豪森噪声信号往往呈现双峰现象，且双峰特征随双相体积占比及应力变化而呈现规律性的变化趋势。
目前，对于含应力双相铁磁性材料磁巴克豪森噪声信号的模拟方法大都采用宏观磁滞模型结合线性叠加方法实现，但该方法无法考虑微观特征信息，也无法考虑两单相间的关系。因此，尚缺乏含应力双相铁磁性材料磁巴克豪森噪声信号的合理模拟方法。
本发明在现有ising模型基础上，创新性地提出利用ising模型自旋平面网格划分和等效自旋相互作用系数分网格设置的方法，实现了双相铁磁性材料的磁巴克豪森噪声信号模拟，定量预测了双相材料组分占比和应力对磁巴克豪森噪声信号双峰峰值的影响规律。
技术实现要素：
本发明的目的在于提供了一种含应力双相铁磁材料中磁巴克豪森噪声信号的模拟方法，以实现在理论上对应力作用下的双相铁磁性材料磁巴克豪森噪声信号进行预测。
为实现上述发明目的，本发明提供的技术方案如下：一种含应力双相铁磁材料中磁巴克豪森噪声信号的模拟方法，该方法的具体实现过程如下：
步骤1：在含应力双相铁磁性材料的磁巴克豪森噪声模拟方法中，采用ising模型结合蒙特卡洛算法的方式。首先对ising模型中的自旋平面进行网格划分，形成n个自旋数目相等的子自旋平面，每个子自旋平面中包含n个自旋点，选择其中q(1≤q≤n)个子自旋平面设置等效自旋相互作用系数为je1，其余n-q个子自旋平面内的自旋相互作用系数设为je2，此时组分占比计为β＝q/n。等效自旋相互作用系数je1和je2的表达式为：
式中j1和j2为无应力时的自旋相互作用系数；jσ1和jσ2为双组分网格内应力σ所致的自旋相互作用系数，jσ1和jσ2与应力σ的关系用四次多项式形式表述如下：
式中p11、p12、p13、p14、p21、p22、p23、p24为多项式各项系数。
在含应力双相铁磁性材料的磁巴克豪森噪声信号模拟中，考虑组分占比和应力的ising模型的总自旋系统的能量e表示为：
式中s表示这两种自旋状态，取值为±1，分别代表自旋向上与自旋向下，n×n表示总自旋点数。
如果不考虑应力(σ＝0)的影响，则式(3)退化为适用于无应力状态下双相铁磁材料的模型，由此形成的总自旋系统的能量e表示为：
步骤2：利用蒙特卡洛算法对式(3)或式(4)所表达的含应力双相铁磁性材料磁巴克豪森噪声信号模拟ising模型整体自旋平面区域进行求解，用自旋系统的算数平均值表述磁化强度m，表达式为：
在含应力双相铁磁性材料的磁巴克豪森噪声信号模拟中，计算磁巴克豪森噪声信号，其表达式为：
式中t为模拟时间，即在蒙特卡罗算法中的求解步数。
不同组分体积占比条件下双相材料的磁巴克豪森噪声信号可通过对式(3)的求解得到。双相材料在应力作用下的磁巴克豪森噪声信号模拟可通过对式(4)的求解得到，并可定量预测组分占比与应力对磁巴克豪森噪声信号的影响规律。
附图说明
图1ising模型自旋平面划分图。
图2模拟得到的不同组分占比下双相材料的磁巴克豪森噪声信号图。
图3不同组分占比下的双相材料磁巴克豪森噪声信号包络图。
图4组分占比β＝4％时特定应力下的磁巴克豪森噪声信号图。
图5双相材料中代表j1相的磁巴克豪森噪声峰值随应力的变化图。
图6双相材料中代表j2相的磁巴克豪森噪声峰值随应力的变化图。
具体实施方式
为使本发明实施的目的、技术方案和优点更加清楚，下面将结合本发明附图和具体实施方式对本发明作进一步说明。
本发明提供一种含应力双相铁磁材料中磁巴克豪森噪声信号的模拟方法，所述建模方法包括以下步骤：
s1：在含应力双相铁磁性材料的磁巴克豪森噪声模拟方法中，对ising模型中的自旋平面进行等尺寸网格划分，以二维20×20自旋平面为例，划分为25个4×4子自旋平面，对其中随机选择的q个子自旋平面内设置自旋相互作用系数为j1＝0.5，其它的25-q个子自旋平面内的自旋相互作用系数设为j2＝0.9，这里取q分别为1、2和3，对应的组分占比β分别为4％、8％和12％，自旋平面示意图如图1所示。
s2：将划分后的自旋平面代入ising模型，计算得到不同组分占比下的磁巴克豪森噪声信号如图2所示。采用滑动平均算法提取磁巴克豪森噪声信号包络曲线，绘制不同组分占比下双相材料的蝶形曲线如图3所示。
s3：对公式(2)中的各项系数进行赋值，p11＝p21＝1.5×10-34mpa-4，p12＝p22＝-1.1×10-25mpa-3，p13＝p23＝2.2×10-17mpa-2和p14＝p24＝5.0×10-15mpa-1，将各项系数代入公式(3)，以β＝4％为例，模拟得到特定应力下的磁巴克豪森噪声信号如图4所示。
s4：在模拟应力为0到200mpa范围内以间隔20mpa进行求解，得到各模拟应力下不同组分占比双相材料的磁巴克豪森噪声信号，提取信号的双峰峰值随应力的变化如图5与图6所示。
1.一种含应力双相铁磁材料中磁巴克豪森噪声信号的模拟方法，其特征在于：该方法的具体实现过程如下，
步骤1：在含应力双相铁磁性材料的磁巴克豪森噪声模拟方法中，采用ising模型结合蒙特卡洛算法的方式；首先对ising模型中的自旋平面进行网格划分，形成n个自旋数目相等的子自旋平面，每个子自旋平面中包含n个自旋点，选择其中q个子自旋平面设置等效自旋相互作用系数为je1，其余n-q个子自旋平面内的自旋相互作用系数设为je2，此时组分占比计为β＝q/n，1≤q≤n；
步骤2：利用蒙特卡洛算法对式(3)或式(4)所表达的含应力双相铁磁性材料磁巴克豪森噪声信号模拟ising模型整体自旋平面区域进行求解，用自旋系统的算数平均值表述磁化强度m。
2.根据权利要求1所述的一种含应力双相铁磁材料中磁巴克豪森噪声信号的模拟方法，其特征在于：等效自旋相互作用系数je1和je2的表达式为：
式中j1和j2为无应力时的自旋相互作用系数；jσ1和jσ2为双组分网格内应力σ所致的自旋相互作用系数，jσ1和jσ2与应力σ的关系用四次多项式形式表述如下：
式中p11、p12、p13、p14、p21、p22、p23、p24为多项式各项系数；
在含应力双相铁磁性材料的磁巴克豪森噪声信号模拟中，考虑组分占比和应力的ising模型的总自旋系统的能量e表示为：
式中s表示这两种自旋状态，取值为±1，分别代表自旋向上与自旋向下，n×n表示总自旋点数；
如果不考虑应力σ＝0的影响，则式(3)退化为适用于无应力状态下双相铁磁材料的模型，由此形成的总自旋系统的能量e表示为：
3.根据权利要求1所述的一种含应力双相铁磁材料中磁巴克豪森噪声信号的模拟方法，其特征在于：m表达式为：
在含应力双相铁磁性材料的磁巴克豪森噪声信号模拟中，计算磁巴克豪森噪声信号，其表达式为：
式中t为模拟时间，即在蒙特卡罗算法中的求解步数。
4.根据权利要求1所述的一种含应力双相铁磁材料中磁巴克豪森噪声信号的模拟方法，其特征在于：对ising模型中的自旋平面进行等尺寸网格划分，网格内自旋阵列设置不同的自旋相互作用系数以模拟具有不同磁特性的组分或者相，对于双相材料而言，自旋平面内共包含两种自旋相互作用系数取值不同的网格类型，两种网格的数量比值代表双相材料内两种组分的体积占比，通过蒙特卡洛算法对ising模型整体自旋平面区域进行求解，即可得到不同组分体积占比条件下双相材料的磁巴克豪森噪声信号。
5.根据权利要求1所述的一种含应力双相铁磁材料中磁巴克豪森噪声信号的模拟方法，其特征在于：借助于等效自旋相互作用系数je将应力σ引入模型，分析应力对双相铁磁材料中磁巴克豪森噪声信号的影响规律。
技术总结