CN111257807B - Simulation method of magnetic Barkhausen noise signal in stress-containing two-phase ferromagnetic material - Google Patents
Simulation method of magnetic Barkhausen noise signal in stress-containing two-phase ferromagnetic material Download PDFInfo
- Publication number
- CN111257807B CN111257807B CN202010106235.4A CN202010106235A CN111257807B CN 111257807 B CN111257807 B CN 111257807B CN 202010106235 A CN202010106235 A CN 202010106235A CN 111257807 B CN111257807 B CN 111257807B
- Authority
- CN
- China
- Prior art keywords
- stress
- spin
- barkhausen noise
- magnetic barkhausen
- noise signal
- 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
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/12—Measuring magnetic properties of articles or specimens of solids or fluids
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
- Magnetic Heads (AREA)
- Hall/Mr Elements (AREA)
Abstract
The invention discloses a method for simulating a magnetic Barkhausen noise signal in a stress-containing dual-phase ferromagnetic material, which comprises the steps of carrying out mesh division on a spinning plane in an Ising model, and carrying out assignment on different model parameters in a divided sub-spinning plane mesh so as to simulate different magnetic properties of two components. The ratio of the number of meshes occupied by each of the two components can be adjusted to simulate variations in the ratio of the components. Solving the Ising model by using a Monte Carlo algorithm to obtain a magnetic Barkhausen noise signal of the two-phase material with different component ratios. And describing the relation between the spin interaction coefficient and the stress in the Ising model by using a fourth-order polynomial, substituting the relation into a biphase material magnetic Barkhausen noise signal simulation model, and analyzing the change rule of the characteristics of the magnetic Barkhausen noise signal under the combined action of the component ratio and the stress. The influence rule of the component proportion and the stress change on the magnetic Barkhausen noise predicted by the method disclosed by the invention is consistent with the experimental result.
Description
Technical Field
A simulation method of a magnetic Barkhausen noise signal in a stress-containing dual-phase ferromagnetic material belongs to the cross content of the fields of magnetic physics and nondestructive testing, is not only helpful for revealing physical phenomena, but also can be used as a theoretical basis for component proportion and stress nondestructive testing of the dual-phase material.
Background
The magnetic Barkhausen noise detection method has the characteristic of nondestructive detection and is very sensitive to the change of the internal component proportion and the stress of the material. In the experimental process, the magnetic Barkhausen noise signal detected in the two-phase ferromagnetic material is observed to show a bimodal phenomenon, and the bimodal characteristic shows a regular change trend along with the change of the two-phase volume ratio and the stress.
At present, a macroscopic hysteresis model is mostly combined with a linear superposition method to realize a simulation method of a magnetic Barkhausen noise signal of a stress-containing two-phase ferromagnetic material, but the method cannot consider microscopic characteristic information and cannot consider the relation between two single phases. Therefore, a reasonable simulation method of the magnetic barkhausen noise signal of the stress-containing two-phase ferromagnetic material is not available.
The invention creatively provides a method for setting grids by utilizing Ising model spinning plane grid division and equivalent spinning interaction coefficients on the basis of the existing Ising model, realizes the simulation of magnetic Barkhausen noise signals of a two-phase ferromagnetic material, and quantitatively predicts the influence rule of the component proportion and stress of the two-phase material on the bimodal peak of the magnetic Barkhausen noise signals.
Disclosure of Invention
The invention aims to provide a method for simulating a magnetic Barkhausen noise signal in a stress-containing dual-phase ferromagnetic material, so as to realize the prediction of the magnetic Barkhausen noise signal of the dual-phase ferromagnetic material under the action of stress theoretically.
In order to achieve the purpose, the technical scheme provided by the invention is as follows: a simulation method of magnetic Barkhausen noise signals in a stress-containing dual-phase ferromagnetic material is specifically realized by the following steps:
step 1: in the magnetic Barkhausen noise simulation method of the stress-containing two-phase ferromagnetic material, an Ising model is combined with a Monte Carlo algorithm. Firstly, carrying out grid division on spin planes in an Ising model to form N spin sub-planes with equal spin number, wherein each spin sub-plane comprises N spin points, and Q (Q is more than or equal to 1 and less than or equal to N) is selected to set an equivalent spin interaction coefficient to be Je1The spin interaction coefficient in the remaining N-Q sub-spin planes is set to Je2In this case, the component ratio is β ═ Q/N. Equivalent spin interaction coefficient Je1And Je2The expression of (a) is:
in the formula J1And J2Is the spin interaction coefficient in the absence of stress; j. the design is a squareσ1And Jσ2Coefficient of spin interaction due to internal stress σ of the bicomponent lattice, Jσ1And Jσ2The relation with stress sigma is four timesThe polynomial form is expressed as follows:
in the formula p11、p12、p13、p14、p21、p22、p23、p24Is the polynomial coefficient.
In the magnetic barkhausen noise signal simulation of the stress-containing two-phase ferromagnetic material, the energy E of the total spin system of the Ising model considering the component ratio and the stress is expressed as follows:
wherein S represents the two spin states, the value of S is +/-1, the spin-up and the spin-down are represented respectively, and Nxn represents the number of total spin points.
If the influence of the stress (σ ═ 0) is not considered, equation (3) degenerates to a model suitable for the two-phase ferromagnetic material in the unstressed state, and the energy E of the total spin system thus formed is expressed as:
step 2: solving the whole spin plane region of the stress-containing two-phase ferromagnetic material magnetic Barkhausen noise signal simulation Ising model expressed by the formula (3) or the formula (4) by using a Monte Carlo algorithm, and expressing the magnetization M by using the arithmetic mean value of a spin system, wherein the expression is as follows:
in the simulation of the magnetic Barkhausen noise signal of the stress-containing two-phase ferromagnetic material, the magnetic Barkhausen noise signal is calculated, and the expression is as follows:
where t is the simulation time, i.e., the number of solution steps in the Monte Carlo algorithm.
The magnetic Barkhausen noise signal of the two-phase material under the condition of different component volume ratios can be obtained by solving the formula (3). The magnetic Barkhausen noise signal simulation of the two-phase material under the stress action can be obtained by solving the formula (4), and the influence rule of the component ratio and the stress on the magnetic Barkhausen noise signal can be quantitatively predicted.
Drawings
Figure 1Ising model spin-plane plots.
Fig. 2 is a graph of the magnetic barkhausen noise signal of the two-phase material under different component ratios obtained by simulation.
FIG. 3 is a signal envelope diagram of magnetic Barkhausen noise of a two-phase material under different component ratios.
Fig. 4 is a graph of magnetic barkhausen noise signals at a specified stress for a component fraction β of 4%.
In the biphasic material of FIG. 5, represents J1Magnetic Barkhausen noise peaks of the phases are plotted against stress.
In the biphasic material of FIG. 6, representative J2Magnetic Barkhausen noise peaks of the phases are plotted against stress.
Detailed Description
In order to make the objectives, technical solutions and advantages of the present invention more apparent, the present invention will be further described with reference to the accompanying drawings and detailed description of the present invention.
The invention provides a method for simulating a magnetic Barkhausen noise signal in a stress-containing dual-phase ferromagnetic material, which comprises the following steps of:
s1: in the magnetic Barkhausen noise simulation method of the stress-containing two-phase ferromagnetic material, the spinning plane in the Ising model is subjected to equal-size grid division, a two-dimensional 20 multiplied by 20 spinning plane is taken as an example and is divided into 25 4 multiplied by 4 spinning planes, and the spinning interaction coefficients are arranged in Q randomly selected spinning planesIs J1The spin interaction coefficient in the other 25-Q spin sub-planes is J, 0.52Where Q is taken to be 1, 2 and 3, respectively, and the corresponding component proportions β are 4%, 8% and 12%, respectively, the schematic diagram of the spin plan is shown in fig. 1.
S2: and substituting the divided spin planes into the Ising model, and calculating to obtain magnetic Barkhausen noise signals with different component ratios as shown in figure 2. The envelope curve of the magnetic Barkhausen noise signal is extracted by adopting a moving average algorithm, and butterfly curves of the two-phase material under different component ratios are drawn and are shown in figure 3.
S3: the coefficients of the terms in equation (2) are assigned, p11=p21=1.5×10-34MPa-4,p12=p22=-1.1×10-25MPa-3,p13=p23=2.2×10-17MPa-2And p14=p24=5.0×10-15MPa-1Fig. 4 shows a magnetic barkhausen noise signal under a specific stress obtained by simulation by substituting each coefficient into equation (3) and taking β as an example of 4%.
S4: solving is carried out at intervals of 20MPa within the range of simulated stress of 0-200 MPa to obtain magnetic Barkhausen noise signals of the two-phase material with different component ratios under each simulated stress, and the variation of the double peak value of the extracted signals along with the stress is shown in figures 5 and 6.
Claims (3)
1. A simulation method of magnetic Barkhausen noise signals in a stress-containing dual-phase ferromagnetic material is characterized by comprising the following steps: the specific implementation process of the method is as follows,
step 1: in the magnetic Barkhausen noise simulation method of the stress-containing two-phase ferromagnetic material, a mode of combining an Ising model with a Monte Carlo algorithm is adopted; firstly, carrying out grid division on spin planes in an Ising model to form N spin sub-planes with equal spin number, wherein each spin sub-plane comprises N spin points, and Q spin sub-planes are selected to set an equivalent spin interaction coefficient to be Je1The spin interaction coefficient in the remaining N-Q sub-spin planes is set to Je2At this timeThe proportion of the components is beta-Q/N, and Q is more than or equal to 1 and less than or equal to N;
step 2: solving the whole spin plane region of the stress-containing two-phase ferromagnetic material magnetic Barkhausen noise signal simulation Ising model expressed by the formula (3) or the formula (4) by using a Monte Carlo algorithm, and expressing the magnetization M by using the arithmetic mean value of a spin system;
coefficient of equivalent spin interaction Je1And Je2The expression of (a) is:
in the formula J1And J2Is the spin interaction coefficient in the absence of stress; j. the design is a squareσ1And Jσ2Coefficient of spin interaction due to internal stress σ of the bicomponent lattice, Jσ1And Jσ2The relationship to stress σ is expressed in the form of a fourth order polynomial as follows:
in the formula p11、p12、p13、p14、p21、p22、p23、p24Is the coefficient of each polynomial;
in the magnetic barkhausen noise signal simulation of the stress-containing two-phase ferromagnetic material, the energy E of the total spin system of the Ising model considering the component ratio and the stress is expressed as follows:
wherein S represents the two spin states, the value of S is +/-1, the spin states represent spin-up and spin-down respectively, and Nxn represents the number of total spin points;
if the influence of the stress σ ═ 0 is not taken into consideration, equation (3) degenerates to a model suitable for the two-phase ferromagnetic material in the stress-free state, and the energy E of the total spin system thus formed is expressed as:
the expression of M is:
in the simulation of the magnetic Barkhausen noise signal of the stress-containing two-phase ferromagnetic material, the magnetic Barkhausen noise signal is calculated, and the expression is as follows:
where t is the simulation time, i.e., the number of solution steps in the Monte Carlo algorithm.
2. The method for simulating the magnetic Barkhausen noise signal in the stress-containing dual-phase ferromagnetic material according to claim 1, wherein the method comprises the following steps: the method comprises the steps of carrying out equal-size grid division on a spinning plane in an Ising model, setting different spinning interaction coefficients for a spinning array in a grid to simulate components or phases with different magnetic characteristics, wherein for a two-phase material, the spinning plane contains two grid types with different spinning interaction coefficient values, the quantity ratio of the two grids represents the volume ratio of the two components in the two-phase material, and solving the whole spinning plane area of the Ising model through a Monte Carlo algorithm to obtain the magnetic Barkhausen noise signals of the two-phase material under the condition of different component volume ratios.
3. The method for simulating the magnetic Barkhausen noise signal in the stress-containing dual-phase ferromagnetic material according to claim 1, wherein the method comprises the following steps: by means of equivalent spin interaction coefficient JeIntroducing the stress sigma into the model, and analyzing the stress-pair biphase ferromagnetic materialInfluence rule of magnetic Barkhausen noise signals in the material.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010106235.4A CN111257807B (en) | 2020-02-21 | 2020-02-21 | Simulation method of magnetic Barkhausen noise signal in stress-containing two-phase ferromagnetic material |
PCT/CN2021/076119 WO2021164627A1 (en) | 2020-02-21 | 2021-02-09 | Method for simulating magnetic barkhausen noise signal in stress-containing dual-phase ferromagnetic material |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010106235.4A CN111257807B (en) | 2020-02-21 | 2020-02-21 | Simulation method of magnetic Barkhausen noise signal in stress-containing two-phase ferromagnetic material |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111257807A CN111257807A (en) | 2020-06-09 |
CN111257807B true CN111257807B (en) | 2022-05-03 |
Family
ID=70942961
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010106235.4A Active CN111257807B (en) | 2020-02-21 | 2020-02-21 | Simulation method of magnetic Barkhausen noise signal in stress-containing two-phase ferromagnetic material |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN111257807B (en) |
WO (1) | WO2021164627A1 (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111257807B (en) * | 2020-02-21 | 2022-05-03 | 北京工业大学 | Simulation method of magnetic Barkhausen noise signal in stress-containing two-phase ferromagnetic material |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11337527A (en) * | 1998-05-28 | 1999-12-10 | Nippon Steel Corp | Detection method of barkhausen noise |
JP2000340860A (en) * | 1999-03-18 | 2000-12-08 | Read Rite Corp | Transducer and its manufacture |
CN1352790A (en) * | 1999-07-07 | 2002-06-05 | 富士通株式会社 | Magnetoresistance effect head, method for producing magnetoresistance effect head, and informatino reproducing device |
CN1410766A (en) * | 2002-11-02 | 2003-04-16 | 东风汽车公司 | Method of proceeding nondestructive inspection using Barkhausen noise signal |
EP1530707A1 (en) * | 2002-08-20 | 2005-05-18 | EADS Deutschland GmbH | Method for measuring stress/strain using barkhausen noises |
CN102347031A (en) * | 2010-07-23 | 2012-02-08 | 新科实业有限公司 | Method for measuring noise of magnetic head and system thereof |
CN105548924A (en) * | 2016-01-28 | 2016-05-04 | 中国特种设备检测研究院 | Magnetic Barkhausen and magnetic parameter sensor and measurement method of magnetic Barkhausen and magnetic parameter |
CN106248781A (en) * | 2016-07-27 | 2016-12-21 | 南京航空航天大学 | A kind of material magnetic characteristic based on Barkhausen's principle and stress mornitoring method |
CN109933914A (en) * | 2019-03-18 | 2019-06-25 | 北京工业大学 | The modeling method of two-phase ferrimagnet magnetic hysteresis and Barkhausen noise signal |
CN110388996A (en) * | 2019-09-03 | 2019-10-29 | 广东电网有限责任公司 | A kind of barkhausen signal characteristic-acquisition method, device, terminal and storage medium |
CN110749846A (en) * | 2019-09-26 | 2020-02-04 | 南京航空航天大学 | Barkhausen signal detection method based on linear frequency modulation excitation |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4599563A (en) * | 1983-10-04 | 1986-07-08 | American Stress Technologies, Inc. | Barkhausen noise method for analyzing the anisotropic properties of ferromagnetic steel |
CN109190247B (en) * | 2018-09-01 | 2020-11-06 | 刘照森 | Application of optimized quantum Monte-Carlo simulation method in research of complex magnetic system |
CN111257807B (en) * | 2020-02-21 | 2022-05-03 | 北京工业大学 | Simulation method of magnetic Barkhausen noise signal in stress-containing two-phase ferromagnetic material |
-
2020
- 2020-02-21 CN CN202010106235.4A patent/CN111257807B/en active Active
-
2021
- 2021-02-09 WO PCT/CN2021/076119 patent/WO2021164627A1/en active Application Filing
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11337527A (en) * | 1998-05-28 | 1999-12-10 | Nippon Steel Corp | Detection method of barkhausen noise |
JP2000340860A (en) * | 1999-03-18 | 2000-12-08 | Read Rite Corp | Transducer and its manufacture |
CN1352790A (en) * | 1999-07-07 | 2002-06-05 | 富士通株式会社 | Magnetoresistance effect head, method for producing magnetoresistance effect head, and informatino reproducing device |
EP1530707A1 (en) * | 2002-08-20 | 2005-05-18 | EADS Deutschland GmbH | Method for measuring stress/strain using barkhausen noises |
CN1410766A (en) * | 2002-11-02 | 2003-04-16 | 东风汽车公司 | Method of proceeding nondestructive inspection using Barkhausen noise signal |
CN102347031A (en) * | 2010-07-23 | 2012-02-08 | 新科实业有限公司 | Method for measuring noise of magnetic head and system thereof |
CN105548924A (en) * | 2016-01-28 | 2016-05-04 | 中国特种设备检测研究院 | Magnetic Barkhausen and magnetic parameter sensor and measurement method of magnetic Barkhausen and magnetic parameter |
CN106248781A (en) * | 2016-07-27 | 2016-12-21 | 南京航空航天大学 | A kind of material magnetic characteristic based on Barkhausen's principle and stress mornitoring method |
CN109933914A (en) * | 2019-03-18 | 2019-06-25 | 北京工业大学 | The modeling method of two-phase ferrimagnet magnetic hysteresis and Barkhausen noise signal |
CN110388996A (en) * | 2019-09-03 | 2019-10-29 | 广东电网有限责任公司 | A kind of barkhausen signal characteristic-acquisition method, device, terminal and storage medium |
CN110749846A (en) * | 2019-09-26 | 2020-02-04 | 南京航空航天大学 | Barkhausen signal detection method based on linear frequency modulation excitation |
Non-Patent Citations (2)
Title |
---|
Monte Carlo simulation for relationship between magnetic Barkhausen noise and elastic stress of steel;Liang Li,etc;《Advances in Mechanical Engineering》;20161203;43-49 * |
基于巴克豪森效应的无损检测关键技术研究;王文涛;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20171115;C029-28 * |
Also Published As
Publication number | Publication date |
---|---|
WO2021164627A1 (en) | 2021-08-26 |
CN111257807A (en) | 2020-06-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gong et al. | Calculating temperature-dependent properties of Nd 2 Fe 14 B permanent magnets by atomistic spin model simulations | |
CN111257807B (en) | Simulation method of magnetic Barkhausen noise signal in stress-containing two-phase ferromagnetic material | |
Nowik | Mössbauer studies of valence fluctuations | |
CN112259176A (en) | Neural network-based sintered SmCo magnetic performance prediction method and system | |
Garcia-Cervera et al. | Adaptive mesh refinement for micromagnetics simulations | |
Ye et al. | Nanoscale patterning of chemical order induced by displacement cascades in irradiated alloys. I. A kinetic Monte Carlo study | |
CN105720950A (en) | Nuclear pulse signal generating method and nuclear pulse signal generating device | |
CN116885726A (en) | Power equipment operation control method and system based on digital twin technology | |
Andrei et al. | Temperature, stress, and rate dependent numerical implementation of magnetization processes in phenomenological models | |
Levine et al. | A lattice dynamics study of a Langmuir monolayer of monounsaturated fatty acids | |
Todate et al. | Spin correlation in the diluted triangular-lattice antiferromagnet LuFeMgO 4 | |
CN116050135A (en) | Superconducting coil simulation method, system and computer storage medium | |
Kplé et al. | The mixed spin-1/2 and spin-1 Ising system on a two-layer Bethe lattice | |
CN103136401B (en) | The EXPERIMENTAL DESIGN list acquirement method of a kind of man-made microstructure and device | |
Mei-Ling et al. | Critical behavior of a dynamical percolation model | |
Chantrell et al. | Magnetic properties and microstructure of thin-film media | |
Gabay et al. | Dynamics of branched domain structures | |
Jiang et al. | Population spatial migration tendency forecasting with Grey model and Fourier series | |
CN109740582A (en) | A kind of power signal noise filtering method and system for Energy Decomposition | |
Kasai et al. | Hidden mattis phase in annealed system | |
Zhang et al. | Temperature-dependent demagnetization nonlinear Wiener model with neural network for PM synchronous machines in electric vehicle | |
Borštnik et al. | Criticality of D= 2 and D= 3 Ising models: Cluster structure versus populations | |
Schram et al. | Identifying monopoles on a lattice | |
Yang et al. | Magnetization and spin excitations of non-abelian quantum hall states | |
Zhang et al. | Bimerons create bimerons: Proliferation and aggregation induced by currents and magnetic fields: Special Collection: Condensed Matter |
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 |