Dynamic magnetostriction determination method and system based on JA model

Abstract

The present invention relates to a dynamic magnetostriction determination method and system based on a JA model, belonging to the research field of magnetostriction characteristics of ferromagnetic materials. The total loss of ferromagnetic materials under dynamic conditions is decomposed into hysteresis loss, eddy current loss and residual loss, and a total magnetic field intensity calculation model under dynamic conditions is obtained by using a field separation method. A magnetostriction model is established in combination with a secondary magnetic domain rotation model. In practical application, the magnetostriction is determined only according to the magnetic induction intensity of the ferromagnetic material to be measured under static conditions, thereby avoiding the defect of calculation deviation caused by using measurement data, and effectively improving the accuracy of dynamic magnetostriction analysis.

Description

Technical Field

[0001] The invention relates to the research field of magnetostrictive characteristics of ferromagnetic materials, and in particular to a dynamic magnetostriction determination method and system based on a JA model. Background Art

[0002] Magnetostriction is one of the important inherent characteristics of magnetic materials. The magnetostrictive characteristics of different magnetic materials are mainly expressed by the butterfly curve of the material. Magnetic materials exist in almost all electrical equipment, such as the iron core of the transformer in the power system, the inductor coil in the generator and electronic circuit, etc. Magnetostriction is the main source of vibration noise of the iron core of power equipment. The international community has been discussing how to effectively reduce the vibration noise of electrical equipment such as transformers and reactors. Therefore, it is necessary to explore the magnetostrictive characteristics of iron core silicon steel sheets, so as to further solve the problem of vibration reduction and noise reduction of electrical equipment.

[0003] At present, laser measurement methods have been widely used in the magnetostriction measurement of silicon steel sheets. Many scholars have conducted certain research on the magnetostriction and anisotropy of oriented and non-oriented silicon steel sheets under stress, harmonics, DC bias and other working conditions, but they are basically based on the analysis of measurement data. Therefore, it is necessary to establish a magnetostriction model so that it can obtain accurate and effective magnetostriction using as little experimental data as possible. Summary of the invention

[0004] The purpose of the present invention is to provide a method and system for determining dynamic magnetostriction based on a JA model, which determines magnetostriction only according to the magnetic induction intensity under static conditions, thereby improving the accuracy of dynamic magnetostriction analysis.

[0005] To achieve the above object, the present invention provides the following solutions:

[0006] A method for determining dynamic magnetostriction based on a JA model, the method comprising:

[0007] Based on the loss statistics theory, the total loss of ferromagnetic materials under dynamic conditions is decomposed into magneto-induced loss, eddy current loss and residual loss.

[0008] Determine the eddy current loss calculation model and the residual loss calculation model respectively;

[0009] According to the eddy current loss calculation model and the residual loss calculation model, based on the field separation theory, the dynamic total magnetic field intensity calculation model is determined;

[0010] According to the dynamic total magnetic field intensity calculation model, combined with the secondary magnetic domain rotation model, a dynamic magnetostriction model is determined; the dynamic magnetostriction model is a relationship model between magnetostriction and magnetic induction intensity;

[0011] Obtaining the magnetic induction intensity of the ferromagnetic material to be tested;

[0012] According to the magnetic induction intensity of the ferromagnetic material to be tested, the magnetostriction of the ferromagnetic material to be tested is obtained by using the dynamic magnetostriction model.

[0013] Optionally, the eddy current loss calculation model is:

[0014]

[0015] Among them, W ed is the eddy current loss, d is the thickness of the silicon steel sheet, ρ is the resistivity, β is the shape parameter, B is the magnetic induction intensity, and t is the time.

[0016] Optionally, the method for determining the residual loss calculation model is:

[0017] The linear relationship between the difference between the total loss and the eddy current loss and the root mean square of the frequency is determined by numerical fitting method;

[0018] According to the slope of the linear relationship, using the formula Determine the statistical parameters that characterize the local magnetic field distribution of the magnet; where K is the slope, H0 is the statistical parameter that characterizes the local magnetic field distribution of the magnet, G is the dimensionless coefficient, w is the lamination width, B m is the AC peak magnetic flux density;

[0019] Substituting the statistical parameters into the formula The residual loss calculation model is obtained; where W an The remaining loss.

[0020] Optionally, the dynamic total magnetic field strength calculation model is determined based on the eddy current loss calculation model and the residual loss calculation model and field separation theory, specifically including:

[0021] According to the eddy current loss calculation model and the residual loss calculation model, based on the field separation theory, the total loss model is obtained as follows: Where W is the total loss, W hy is the magneto-induced loss, H hy is the static magnetic field strength;

[0022] According to the total loss model, the dynamic total magnetic field strength calculation model is determined as follows: Among them, H total is the dynamic total magnetic field intensity, and γ is the coefficient.

[0023] Optionally, the dynamic magnetostriction model is determined based on the dynamic total magnetic field intensity calculation model in combination with the secondary magnetic domain rotation model, specifically including:

[0024] Based on the dynamic total magnetic field strength calculation model, combined with the secondary magnetic domain rotation model and B = μ0(H total +M), the dynamic magnetostrictive model is determined as

[0025] Among them, λ is magnetostriction, λ s is the saturation magnetostriction, M is the magnetization intensity, M s is the saturation magnetization, μ0 is the vacuum permeability.

[0026] Optionally, obtaining the magnetostriction of the ferromagnetic material to be tested by using the dynamic magnetostriction model according to the magnetic induction intensity of the ferromagnetic material to be tested specifically includes:

[0027] According to the magnetic induction intensity of the ferromagnetic material to be tested, the static magnetic field intensity is obtained by using the static JA hysteresis model;

[0028] The magnetostriction of the ferromagnetic material to be tested is obtained by using the dynamic magnetostriction model according to the static magnetic field intensity and the magnetic induction intensity of the ferromagnetic material to be tested.

[0029] Optionally, the static magnetic field strength is obtained by using a static JA hysteresis model according to the magnetic induction strength of the ferromagnetic material to be measured, specifically including:

[0030] Based on the particle swarm optimization algorithm, the parameters of the static JA hysteresis model are extracted and the static JA hysteresis model is determined as Among them, δ M To prevent the occurrence of non-physical solutions, M an is the hysteresis-free magnetization intensity, c is the reversible magnetization parameter, k is the pinning parameter between magnetic domains, δ is the direction coefficient, H e is the effective magnetic field intensity, α is the mean field parameter of the coupling inside the magnetic domain;

[0031] According to the magnetic induction intensity of the ferromagnetic material to be tested, the static JA hysteresis model is solved by the fourth-order Runge-Kutta method to obtain the magnetization intensity;

[0032] According to the magnetization intensity, using the formula B = μ0 (H hy +M), determine the static magnetic field strength.

[0033] A dynamic magnetostriction determination system based on a JA model, the system comprising:

[0034] The loss decomposition module is used to decompose the total loss of ferromagnetic materials under dynamic conditions into magneto-induced loss, eddy current loss and residual loss based on the loss statistics theory;

[0035] A calculation model determination module, used to determine an eddy current loss calculation model and a residual loss calculation model respectively;

[0036] A dynamic total magnetic field strength calculation model determination module is used to determine the dynamic total magnetic field strength calculation model based on the eddy current loss calculation model and the residual loss calculation model and the field separation theory;

[0037] A dynamic magnetostriction model determination module is used to determine a dynamic magnetostriction model based on a dynamic total magnetic field intensity calculation model in combination with a secondary magnetic domain rotation model; the dynamic magnetostriction model is a relationship model between magnetostriction and magnetic induction intensity;

[0038] A magnetic induction intensity acquisition module is used to obtain the magnetic induction intensity of the ferromagnetic material to be tested;

[0039] The magnetostriction obtaining module is used to obtain the magnetostriction of the ferromagnetic material to be tested according to the magnetic induction intensity of the ferromagnetic material to be tested and using the dynamic magnetostriction model.

[0040] Optionally, the dynamic total magnetic field strength calculation model determination module specifically includes:

[0041] The total loss model acquisition submodule is used to obtain the total loss model based on the eddy current loss calculation model and the residual loss calculation model based on the field separation theory: Where W is the total loss, W hy is the magneto-induced loss, H hy is the static magnetic field strength;

[0042] The dynamic total magnetic field strength calculation model acquisition submodule is used to determine the dynamic total magnetic field strength calculation model according to the total loss model: Among them, H total is the dynamic total magnetic field intensity, and γ is the coefficient.

[0043] Optionally, the dynamic magnetostriction model determination module specifically includes:

[0044] Dynamic magnetostrictive model determination submodule, used to calculate the model based on the dynamic total magnetic field strength, combined with the secondary magnetic domain rotation model and B = μ0(H total +M), the dynamic magnetostrictive model is determined as

[0045] Among them, λ is magnetostriction, λ s is the saturation magnetostriction, M is the magnetization intensity, M s is the saturation magnetization, μ0 is the vacuum permeability.

[0046] According to the specific embodiments provided by the present invention, the present invention discloses the following technical effects:

[0047] The present invention discloses a method and system for determining dynamic magnetostriction based on a JA model. The method decomposes the total loss of a ferromagnetic material under dynamic conditions into hysteresis loss, eddy current loss and residual loss, and uses a field separation method to obtain a total magnetic field intensity calculation model under dynamic conditions. A magnetostriction model is established in combination with a secondary magnetic domain rotation model. In practical applications, the magnetostriction is determined only according to the magnetic induction intensity of the ferromagnetic material to be measured under static conditions, thereby avoiding the defect of calculation deviation caused by using measurement data, and effectively improving the accuracy of dynamic magnetostriction analysis. BRIEF DESCRIPTION OF THE DRAWINGS

[0048] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.

[0049] Figure 1 A flow chart of a method for determining dynamic magnetostriction based on a JA model provided by the present invention;

[0050] Figure 2 A schematic diagram of a method for determining dynamic magnetostriction based on a JA model provided by the present invention;

[0051] Figure 3 This is a flow chart of the particle swarm optimization algorithm provided by the present invention. DETAILED DESCRIPTION

[0052] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0053] The purpose of the present invention is to provide a method and system for determining dynamic magnetostriction based on a JA model, which determines magnetostriction only according to the magnetic induction intensity under static conditions, thereby improving the accuracy of dynamic magnetostriction analysis.

[0054] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the present invention is further described in detail below with reference to the accompanying drawings and specific embodiments.

[0055] The present invention provides a method for determining dynamic magnetostriction based on the JA model, such as Figure 1-2 As shown, the method includes:

[0056] Step 101, based on the loss statistics theory, the total loss of the ferromagnetic material under dynamic conditions is decomposed into magneto-induced loss, eddy current loss and residual loss.

[0057] Step 102, respectively determine an eddy current loss calculation model and a residual loss calculation model.

[0058] (1) The eddy current loss calculation model is:

[0059]

[0060] Among them, W ed is the eddy current loss, d is the thickness of the silicon steel sheet, ρ is the resistivity, β is the shape parameter, B is the magnetic induction intensity, and t is the time.

[0061] The expression of eddy current loss is further obtained as

[0062] (2) The residual loss expression is:

[0063]

[0064] In the residual loss expression, the value of H0 is related to the AC peak magnetic flux density B m It is related and cannot be solved directly. total and hysteresis loss W hy , eddy current loss W ed , Residual loss W an The difference between the total loss and the eddy current loss is W. total -W ed and the root mean square of the frequency f 0.5 There is a linear relationship between them. The intersection of the function corresponding to this linear relationship and the ordinate is the hysteresis loss W hy Therefore, after finding the slope, the size of H0 can be quickly calculated.

[0065] The method for determining the residual loss calculation model is:

[0066] The linear relationship between the difference between the total loss and the eddy current loss and the root mean square of the frequency is determined by numerical fitting method;

[0067] According to the slope of the linear relationship, using the formula Determine the statistical parameters that characterize the local magnetic field distribution of the magnet; where K is the slope, H0 is the statistical parameter that characterizes the local magnetic field distribution of the magnet, G is the dimensionless coefficient, w is the lamination width, B m is the AC peak magnetic flux density;

[0068] Substituting statistical parameters into the formula The residual loss calculation model is obtained; where W anThe remaining loss.

[0069] Step 103, determining a dynamic total magnetic field intensity calculation model according to the eddy current loss calculation model and the residual loss calculation model based on the field separation theory.

[0070] Specifically include:

[0071] According to the eddy current loss calculation model and the residual loss calculation model, based on the field separation theory, the total loss model is obtained as follows: Where W is the total loss, W hy is the magneto-induced loss, H hy is the static magnetic field strength;

[0072] According to the total loss model, the dynamic total magnetic field strength calculation model is determined as follows: Among them, H total is the dynamic total magnetic field intensity, and γ is the coefficient.

[0073] Step 104, determining a dynamic magnetostriction model based on the dynamic total magnetic field intensity calculation model and the secondary magnetic domain rotation model; the dynamic magnetostriction model is a relationship model between magnetostriction and magnetic induction intensity.

[0074] Specifically include:

[0075] Based on the dynamic total magnetic field strength calculation model, combined with the secondary magnetic domain rotation model and B = μ0(H total +M), the dynamic magnetostrictive model is determined as

[0076] Among them, λ is magnetostriction, λ s is the saturation magnetostriction, M is the magnetization intensity, M s is the saturation magnetization, μ0 is the vacuum permeability.

[0077] Step 105, obtaining the magnetic induction intensity of the ferromagnetic material to be tested;

[0078] Step 106, obtaining the magnetostriction of the ferromagnetic material to be tested by using a dynamic magnetostriction model according to the magnetic induction intensity of the ferromagnetic material to be tested.

[0079] Specifically include:

[0080] Step 106-1, according to the magnetic induction intensity of the ferromagnetic material to be tested, using the static JA hysteresis model, obtain the static magnetic field intensity, the implementation process is as follows:

[0081] Based on the particle swarm optimization algorithm, the parameters of the static JA hysteresis model are extracted and the static JA hysteresis model is determined as Among them, δ M To prevent the occurrence of non-physical solutions, Man is the hysteresis-free magnetization intensity, c is the reversible magnetization parameter, k is the pinning parameter between magnetic domains, δ is the direction coefficient, H e is the effective magnetic field intensity, α is the mean field parameter of the coupling inside the magnetic domain;

[0082] According to the magnetic induction intensity of the ferromagnetic material to be tested, the fourth-order Runge-Kutta method is used to solve the static JA hysteresis model to obtain the magnetization intensity;

[0083] According to the magnetization intensity, using the formula B = μ0 (H hy +M), determine the static magnetic field strength.

[0084] Reference Figure 3 , Based on the particle swarm optimization algorithm, the method of extracting the parameters of the static JA hysteresis model is as follows:

[0085] a: Initialize the population, including the particle position X(a, α, M s , k, c) and speed V (a, α, M s , k, c). a represents the shape parameter of the hysteresis-free magnetization curve.

[0086] b: Evaluate the fitness of each particle according to the fitness function, compare its current fitness with the fitness corresponding to its individual historical best position pb, and if the current fitness is higher, the current position will be used to update the historical best position.

[0087] c: For each particle, compare its current fitness value with the fitness value corresponding to the global best position gb. If the current fitness is higher, the global best position will be updated with the current particle position.

[0088] d: Update the speed and position of each particle according to equations (1) and (2). If the end condition is not met, return to step b. If the end condition is met, the algorithm stops.

[0089] v i =ωv i +c1random(0,1)(pb i -x i )+c2random(0,1)(gb-x i )

[0090] x i =x i +v i

[0091] Among them: ω is the inertia factor, c1 and c2 are learning factors, v i is the velocity of the ith particle, random(0,1) represents a random number between (0,1), pb iis the best historical position of the ith particle, x i is the current position of the ith particle.

[0092] Step 106 - 2 , obtaining the magnetostriction of the ferromagnetic material to be tested by using a dynamic magnetostriction model according to the static magnetic field intensity and the magnetic induction intensity of the ferromagnetic material to be tested.

[0093] The present invention decomposes the total loss of ferromagnetic materials under dynamic conditions into hysteresis loss, eddy current loss and residual loss, and uses a field separation method to obtain the total magnetic field intensity under dynamic conditions, and establishes a magnetostriction model in combination with a secondary magnetic domain rotation model, and the algorithm is simple; secondly, the present invention only requires the magnetic induction intensity under static conditions, and requires less experimental data, thus avoiding complex experimental design and measurement work and the inefficiency and errors caused thereby, and effectively improving the accuracy of generating magnetostriction characteristic simulation results.

[0094] The present invention also provides a dynamic magnetostriction determination system based on the JA model, the system comprising:

[0095] The loss decomposition module is used to decompose the total loss of ferromagnetic materials under dynamic conditions into magneto-induced loss, eddy current loss and residual loss based on the loss statistics theory;

[0096] A calculation model determination module, used to determine an eddy current loss calculation model and a residual loss calculation model respectively;

[0097] A dynamic total magnetic field strength calculation model determination module is used to determine the dynamic total magnetic field strength calculation model based on the eddy current loss calculation model and the residual loss calculation model and the field separation theory;

[0098] A dynamic magnetostriction model determination module is used to determine a dynamic magnetostriction model based on a dynamic total magnetic field intensity calculation model and a secondary magnetic domain rotation model; the dynamic magnetostriction model is a relationship model between magnetostriction and magnetic induction intensity;

[0099] A magnetic induction intensity acquisition module is used to obtain the magnetic induction intensity of the ferromagnetic material to be tested;

[0100] The magnetostriction obtaining module is used to obtain the magnetostriction of the ferromagnetic material to be tested by using a dynamic magnetostriction model according to the magnetic induction intensity of the ferromagnetic material to be tested.

[0101] The dynamic total magnetic field strength calculation model determination module specifically includes:

[0102] The total loss model acquisition submodule is used to obtain the total loss model based on the eddy current loss calculation model and the residual loss calculation model based on the field separation theory: Where W is the total loss, W hy is the magneto-induced loss, H hyis the static magnetic field strength;

[0103] The dynamic total magnetic field strength calculation model acquisition submodule is used to determine the dynamic total magnetic field strength calculation model according to the total loss model: Among them, H total is the dynamic total magnetic field intensity, and γ is the coefficient.

[0104] The dynamic magnetostrictive model determination module specifically includes:

[0105] Dynamic magnetostriction model determination submodule, used to calculate the model based on the dynamic total magnetic field strength, combined with the secondary magnetic domain rotation model and B = μ0(H total +M), the dynamic magnetostrictive model is determined as

[0106] Among them, λ is magnetostriction, λ s is the saturation magnetostriction, M is the magnetization intensity, M s is the saturation magnetization, μ0 is the vacuum permeability.

[0107] In this specification, each embodiment is described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the embodiments can be referred to each other. For the system disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant parts can be referred to the method part.

[0108] The principles and implementation methods of the present invention are described in this article using specific examples. The description of the above embodiments is only used to help understand the method and core idea of ​​the present invention. At the same time, for those skilled in the art, according to the idea of ​​the present invention, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as limiting the present invention.

Claims

1. A method for determining dynamic magnetostriction based on a JA model, characterized in that: The method comprises: Based on the loss statistics theory, the total loss of ferromagnetic materials under dynamic conditions is decomposed into hysteresis loss, eddy current loss and residual loss. Determine the eddy current loss calculation model and the residual loss calculation model respectively; According to the eddy current loss calculation model and the residual loss calculation model, based on the field separation theory, the dynamic total magnetic field intensity calculation model is determined; According to the dynamic total magnetic field intensity calculation model, combined with the secondary magnetic domain rotation model, a dynamic magnetostriction model is determined; the dynamic magnetostriction model is a relationship model between magnetostriction and magnetic induction intensity; Obtaining the magnetic induction intensity of the ferromagnetic material to be tested; According to the magnetic induction intensity of the ferromagnetic material to be tested, the magnetostriction of the ferromagnetic material to be tested is obtained by using the dynamic magnetostriction model; The dynamic magnetostriction model is determined based on the dynamic total magnetic field intensity calculation model and the secondary magnetic domain rotation model, and specifically includes: Based on the dynamic total magnetic field strength calculation model, combined with the secondary magnetic domain rotation model and B = μ0(H total +M), the dynamic magnetostrictive model is determined as Among them, λ is magnetostriction, λ s is the saturation magnetostriction, M is the magnetization intensity, M s is the saturation magnetization, μ0 is the vacuum magnetic permeability; B is the magnetic induction intensity, H total is the dynamic total magnetic field strength, H hy is the static magnetic field intensity, d is the thickness of the silicon steel sheet, ρ is the resistivity, β is the shape parameter, γ is the coefficient, G is the dimensionless coefficient, w is the lamination width, and H0 is the statistical parameter characterizing the local magnetic field distribution of the magnet.

2. The method for determining dynamic magnetostriction based on the JA model according to claim 1, characterized in that: The eddy current loss calculation model is: Among them, W ed is the eddy current loss, d is the thickness of the silicon steel sheet, ρ is the resistivity, β is the shape parameter, B is the magnetic induction intensity, and t is the time.

3. The method for determining dynamic magnetostriction based on the JA model according to claim 2, characterized in that: The method for determining the residual loss calculation model is: The linear relationship between the difference between the total loss and the eddy current loss and the root mean square of the frequency is determined by numerical fitting method; According to the slope of the linear relationship, using the formula Determine the statistical parameters that characterize the local magnetic field distribution of the magnet; where K is the slope, H0 is the statistical parameter that characterizes the local magnetic field distribution of the magnet, G is the dimensionless coefficient, w is the lamination width, B m is the AC peak magnetic flux density; Substituting the statistical parameters into the formula The residual loss calculation model is obtained; where W an The remaining loss.

4. The method for determining dynamic magnetostriction based on the JA model according to claim 3, characterized in that: The dynamic total magnetic field strength calculation model is determined based on the eddy current loss calculation model and the residual loss calculation model and the field separation theory, and specifically includes: According to the eddy current loss calculation model and the residual loss calculation model, based on the field separation theory, the total loss model is obtained as follows: Where W is the total loss, W hy is the hysteresis loss, H hy is the static magnetic field strength; According to the total loss model, the dynamic total magnetic field strength calculation model is determined as follows: Among them, H total is the dynamic total magnetic field intensity, and γ is the coefficient.

5. The method for determining dynamic magnetostriction based on the JA model according to claim 4, characterized in that: The method of obtaining the magnetostriction of the ferromagnetic material to be tested by using the dynamic magnetostriction model according to the magnetic induction intensity of the ferromagnetic material to be tested specifically includes: According to the magnetic induction intensity of the ferromagnetic material to be tested, the static magnetic field intensity is obtained by using the static JA hysteresis model; The magnetostriction of the ferromagnetic material to be tested is obtained by using the dynamic magnetostriction model according to the static magnetic field intensity and the magnetic induction intensity of the ferromagnetic material to be tested.

6. The method for determining dynamic magnetostriction based on the JA model according to claim 5, characterized in that: The method of obtaining the static magnetic field strength according to the magnetic induction strength of the ferromagnetic material to be measured by using the static JA hysteresis model specifically includes: Based on the particle swarm optimization algorithm, the parameters of the static JA hysteresis model are extracted and the static JA hysteresis model is determined as Among them, δ M To prevent the occurrence of non-physical solutions, M an is the hysteresis-free magnetization intensity, c is the reversible magnetization parameter, k is the pinning parameter between magnetic domains, δ is the direction coefficient, H e is the effective magnetic field intensity, α is the mean field parameter of the coupling inside the magnetic domain; According to the magnetic induction intensity of the ferromagnetic material to be tested, the static JA hysteresis model is solved by the fourth-order Runge-Kutta method to obtain the magnetization intensity; According to the magnetization intensity, using the formula B = μ0 (H hy +M), determine the static magnetic field strength.

7. A dynamic magnetostriction determination system based on JA model, characterized in that: The system comprises: Loss decomposition module, which is used to decompose the total loss of ferromagnetic materials under dynamic conditions into hysteresis loss, eddy current loss and residual loss based on loss statistics theory; A calculation model determination module, used to determine an eddy current loss calculation model and a residual loss calculation model respectively; A dynamic total magnetic field strength calculation model determination module is used to determine the dynamic total magnetic field strength calculation model based on the eddy current loss calculation model and the residual loss calculation model and the field separation theory; A dynamic magnetostriction model determination module is used to determine a dynamic magnetostriction model based on a dynamic total magnetic field intensity calculation model in combination with a secondary magnetic domain rotation model; the dynamic magnetostriction model is a relationship model between magnetostriction and magnetic induction intensity; A magnetic induction intensity acquisition module is used to obtain the magnetic induction intensity of the ferromagnetic material to be tested; A magnetostriction obtaining module, used to obtain the magnetostriction of the ferromagnetic material to be tested by using the dynamic magnetostriction model according to the magnetic induction intensity of the ferromagnetic material to be tested; The dynamic magnetostrictive model determination module specifically includes: Dynamic magnetostrictive model determination submodule, used to calculate the model based on the dynamic total magnetic field strength, combined with the secondary magnetic domain rotation model and B = μ0(H total +M), the dynamic magnetostrictive model is determined as Among them, λ is magnetostriction, λ s is the saturation magnetostriction, M is the magnetization intensity, M s is the saturation magnetization, μ0 is the vacuum magnetic permeability; B is the magnetic induction intensity, H total is the dynamic total magnetic field strength, H hy is the static magnetic field intensity, d is the thickness of the silicon steel sheet, ρ is the resistivity, β is the shape parameter, γ is the coefficient, G is the dimensionless coefficient, w is the lamination width, and H0 is the statistical parameter characterizing the local magnetic field distribution of the magnet.

8. The dynamic magnetostriction determination system based on the JA model according to claim 7, characterized in that: The dynamic total magnetic field strength calculation model determination module specifically includes: The total loss model acquisition submodule is used to obtain the total loss model based on the eddy current loss calculation model and the residual loss calculation model based on the field separation theory: Where W is the total loss, W hy is the hysteresis loss, H hy is the static magnetic field strength; The dynamic total magnetic field strength calculation model acquisition submodule is used to determine the dynamic total magnetic field strength calculation model according to the total loss model: Among them, H total is the dynamic total magnetic field intensity, and γ is the coefficient.

Images