Ferromagnetic cavity physical force magnetic coupling modeling method

Abstract

The invention discloses a ferromagnetic cavity physical and magnetic coupling modeling method, and relates to the technical field of electromagnetic property analysis, the method adopts a COMSOL Multiphysics finite element simulation technology, and combines hysteresis loop measurement of a ferromagnetic material, a Jiles-Atherton parameter recognition algorithm, a tensile test, residual magnetism measurement and a force and magnetic coupling test, so that the ferromagnetic cavity physical and magnetic coupling modeling method is realized. According to the method, the magnetization behaviors of the ferromagnetic cavity under the action of different stresses and magnetic fields are accurately simulated, the high-precision force-magnetic coupling simulation model is established by automatically identifying five Jiles-Atherton parameters of the ferromagnetic material and combining mechanics and residual magnetism data, space magnetic field distribution under the specified working condition is rapidly calculated, and the method has the advantages of being high in modeling efficiency and high in calculation precision.

Description

Technical Field

[0001] The present invention relates to the technical field of electromagnetic characteristic analysis, and in particular to a method for modeling force-magnetic coupling of a ferromagnetic cavity. Background Art

[0002] In the deep sea, there are ferromagnetic cavities such as functional pipelines and various underwater vehicles. Due to their concealment underwater, magnetism can generally be used to evaluate their status. The magnetization characteristics of ferromagnetic materials are affected by stress in addition to their own magnetism and the earth's magnetic field. Therefore, studying the magnetomechanical coupling characteristics of ferromagnetic cavities is important for analyzing their underwater magnetic field characteristics and then evaluating their status.

[0003] The magnetomechanical coupling effect was first studied in 1945, which found that magnetization is a lagging process that contains irreversible components. Subsequent experiments showed that the effect of stress on magnetic domains is discontinuous. Based on the experimental results, the proximity principle was proposed and verified by subsequent experiments. Based on the proximity principle, researchers proposed the Smith uniform energy model, the Stoner-Wolffatrth model, and the currently widely used Jiles-Atherton model.

[0004] At present, the magnetic field modeling methods for ferromagnetic cavities mainly include: 1) establishing a mapping model of the magnetic field inside and outside the ferromagnetic cavity, and then directly calculating the unknown external magnetic field from the known internal magnetic field; 2) solving the Laplace boundary value problem in the half space through the vector magnetic potential distribution, and constructing a magnetic field prediction model for the ferromagnetic cavity; 3) using the magnet simulation method, the induced magnetic field of the ferromagnetic cavity is regarded as the superposition of three orthogonal magnetic dipoles, and according to the change law of the ellipsoid model of the ferromagnetic cavity in the external field, the analytical expression is obtained; 4) for the magnetic field simulation problem of thin-shell magnetic targets such as ferromagnetic cavities, a multi-model magnetic field simulation method based on COMSOL Multiphysics is proposed, the ferromagnetic cavity is equivalent to a hollow cylindrical model, and the size and distribution of its magnetic field are obtained through simulation calculation; 5) using the boundary integral equation to establish a mathematical model of the ferromagnetic cavity extension, when determining the model parameters from the magnetic field measurement values, the boundary conditions of the model parameters are considered, and the Tikhonov regularization method is used to solve the model parameters to effectively suppress the interference of measurement noise on the calculated values.

[0005] However, the above methods do not consider the influence of the magnetomechanical coupling effect of the ferromagnetic cavity on the spatial magnetic field, so the calculation accuracy is low under the force condition of the ferromagnetic cavity. Summary of the invention

[0006] The purpose of the present invention is to solve the problem that the magnetomechanical coupling effect is not considered in the existing ferromagnetic cavity magnetic field modeling method.

[0007] The technical solution of the present invention is to provide a ferromagnetic cavity force-magnetic coupling modeling method, the method comprising:

[0008] S1. Construct a finite element magnetoelastic coupling model framework, including a model input module, a magnetoelastic coupling parameter automatic identification module, a residual magnetism inversion module, a magnetoelastic coupling calculation module and a model output module; wherein the model input module, the magnetoelastic coupling parameter automatic identification module and the residual magnetism inversion module all receive external data and input the results into the magnetoelastic coupling calculation module, the magnetoelastic coupling calculation module inputs the results into the model output module, and the model output module outputs the magnetoelastic coupling model of the ferromagnetic cavity;

[0009] S2, measuring the geometric parameters of the ferromagnetic cavity, inputting the measured geometric parameters into the model input module, and establishing a geometric model of equal scale in COMSOL Multiphysics finite element software;

[0010] S3, measuring the saturation hysteresis loop of the material sample of the ferromagnetic cavity, bringing the measured saturation hysteresis loop result into the Jiles-Atherton parameter identification algorithm, and inputting the five finally obtained Jiles-Atherton parameters into the model input module;

[0011] S4, measuring the Young's modulus of the material sample of the ferromagnetic cavity, and inputting the measured Young's modulus into the model input module;

[0012] S5, performing a residual magnetism measurement test on the ferromagnetic cavity in an environment where the background magnetic field is 0, and inputting the result of the residual magnetism measurement test into a residual magnetism inversion module;

[0013] S6. Performing a force-magnetic coupling test on the ferromagnetic cavity in an environment where the background magnetic field is 0, and inputting the force-magnetic coupling test result into the force-magnetic coupling parameter automatic identification module;

[0014] S7, the magnetomechanical coupling calculation module outputs the magnetic field of a specified point under a specified working condition, the magnetic field cloud map of a specified surface, and the magnetization intensity distribution cloud map of the ferromagnetic cavity according to the data from the model input module, the magnetomechanical coupling parameter automatic identification module, and the residual magnetism inversion module;

[0015] The specified working condition is determined by the pressure function y=press(t) and the time variable t, the specified background magnetic field is input through the model, and the specified point and the specified surface are determined by the spatial coordinates;

[0016] S8. The model output module generates and outputs the ferromagnetic cavity magnetomechanical coupling model according to the input results of the magnetomechanical coupling calculation module.

[0017] In any of the above technical solutions, further, the method for measuring the saturation hysteresis loop in step S3 includes: cutting the ferromagnetic material into cubes with an edge length of 1 mm and sticking them to the center of the test rod, adjusting the scanning step size and other parameters of the superconducting quantum interference instrument so that the maximum magnetic field exceeds the magnetic field size of the material saturation, and using the superconducting quantum interference instrument to measure the saturation hysteresis loop of the sample.

[0018] In any of the above technical solutions, further, the Jiles-Atherton parameter identification algorithm specifically includes:

[0019] The magnetization formula of the Jiles-Atherton model is:

[0020]

[0021] Among them, M s is the saturation magnetization, k is the internal coupling coefficient, c is the reversibility coefficient, and M is the magnetization of the material under an external magnetic field H;

[0022] M irr is the irreversible magnetization intensity, which can be expressed as:

[0023]

[0024] Among them, H eff =H-αM is the effective magnetic field, α is the shape factor, which represents the microscopic properties of the material, and b is the elastic coefficient of the magnetic domain wall;

[0025] The total magnetization consists of a reversible and an irreversible part:

[0026] M=M irr +c·(M s -M irr );

[0027] The particle swarm optimization algorithm is used to adjust M s , α, k, c, b are used to make the model calculation results match the experimentally measured hysteresis loop. The objective function E is defined as:

[0028]

[0029] Among them, M exp (H i ) represents the experimentally measured magnetization intensity, M calc (H i ) represents the magnetization intensity calculated by the model;

[0030] The objective function E is to minimize the experimental data M exp (H i ) and the calculated result M calc (Hi ) to adjust the model parameters so that the calculated results of the model match the experimental data.

[0031] In any of the above technical solutions, further, the method for measuring the Young's modulus in step S4 includes: cutting the ferromagnetic material into long strip samples of 1 cm×1 cm×10 cm, using a tensile testing machine to stretch both ends of the 10 cm edge of the sample to obtain the Young's modulus of the sample.

[0032] In any of the above technical solutions, further, the residual magnetism measurement test step in step S5 includes: adjusting the compensation current to make the background magnetic field 0, arranging a three-phase fluxgate sensor at a specified point, measuring the fixed magnetic field generated by the residual magnetism of the ferromagnetic body, and obtaining the fixed magnetic field distribution of the ferromagnetic cavity.

[0033] In any of the above technical solutions, further, the remanent magnetism inversion module adopts the "conversion matrix method" to measure the remanent magnetism of the ferromagnetic body, and the relationship between the magnetization intensity matrix M0 and the magnetic induction intensity matrix B0 measured by the sensor can be expressed as:

[0034] M0=AB0;

[0035] Where A is the transformation matrix, expressed as:

[0036] A=B -1 M;

[0037] Among them, B -1 is the magnetic field calculated in COMSOL Multiphysics finite element software, and M is the magnetization intensity matrix obtained using random function in COMSOL Multiphysics finite element software.

[0038] In any of the above technical solutions, further, the steps of the magnetomechanical coupling test in step S5 include: after adjusting the background magnetic field to 0, pressurizing the ferromagnetic cavity and recording the spatial magnetic field once every 2MPa until it reaches 10Mpa; then decompressing the ferromagnetic cavity and recording the spatial magnetic field once every 2MPa until the pressure reaches 0.

[0039] In any of the above technical solutions, further, in the magneto-mechanical coupling parameter automatic identification module, the magneto-mechanical coupling test results are input, and the equivalent field effect parameter expression to be optimized is as follows:

[0040]

[0041] Among them, α σx , α σy , α σz is the equivalent field effect parameter in the x, y, and z directions, α is the field effect parameter of the ferromagnetic material, σ x , σy , σ z is the stress magnitude in the x, y, and z directions, calculated by COMSOL Multiphysics finite element software, μ0 is the vacuum permeability, γ 1x ~γ 1z ’ are the six parameters to be optimized. Similarly, the particle swarm algorithm is used to optimize the six parameters to achieve matching with the magnetomechanical coupling test data.

[0042] The beneficial effects of the present invention are:

[0043] The technical solution of the present invention adopts the Jiles-Atherton parameter automatic identification algorithm, and the Jiles-Atherton parameters of the ferromagnetic material can be automatically identified by inputting the saturation hysteresis loop measurement results; the residual magnetism inversion module and the magnetomechanical coupling parameter automatic identification module only need to input the test results to automatically obtain the residual magnetism distribution and equivalent field effect parameters of the ferromagnetic cavity;

[0044] Using innovative theories, high calculation accuracy: The traditional Jiles-Atherton theory is only used to describe the hysteresis loop of ferromagnetic materials, and does not consider the magnetomechanical coupling effect of ferromagnets. This modeling method uses the improved Jiles-Atherton theory based on equivalent field effect parameters to realize the finite element simulation of magnetomechanical coupling of ferromagnetic cavities, effectively solving the problem that magnetomechanical coupling effect of ferromagnetic materials cannot be realized, and has high calculation accuracy.

[0045] Using module packaging, it is highly convenient: This modeling method divides the complex magnetoelastic coupling modeling process into four modules, namely, model input module, magnetoelastic coupling parameter automatic identification module, residual magnetism inversion module and magnetoelastic coupling calculation module. Parameters are directly transferred between different modules. Users only need to complete the necessary measurement tests and input the results and target working conditions into the finite element model to obtain the target magnetic field under the specified working conditions. Therefore, the operation is simple and easy to use. BRIEF DESCRIPTION OF THE DRAWINGS

[0046] The advantages of the above and additional aspects of the present invention will become apparent and easily understood in the description of the embodiments in conjunction with the following drawings, in which:

[0047] Figure 1 is a finite element force-magnetic coupling model framework diagram of a ferromagnetic cavity force-magnetic coupling modeling method according to an embodiment of the present invention;

[0048] Figure 2 is a flowchart of the steps of a method for modeling ferromagnetic cavity force-magnetic coupling according to an embodiment of the present invention;

[0049] Figure 3is a schematic diagram of a force-magnetic coupling model of an oxygen cylinder according to a ferromagnetic cavity force-magnetic coupling modeling method according to an embodiment of the present invention;

[0050] Figure 4 is a schematic diagram of the fixed magnetic field distribution of an oxygen cylinder according to a ferromagnetic cavity force magnetic coupling modeling method according to an embodiment of the present invention;

[0051] Figure 5 It is a schematic diagram comparing simulation and experimental magnetic field data of a ferromagnetic cavity force-magnetic coupling modeling method according to an embodiment of the present invention. DETAILED DESCRIPTION

[0052] In order to more clearly understand the above-mentioned purpose, features and advantages of the present invention, the present invention is further described in detail below in conjunction with the accompanying drawings and specific embodiments. It should be noted that the embodiments of the present invention and the features in the embodiments can be combined with each other without conflict.

[0053] In the following description, many specific details are set forth to facilitate a full understanding of the present invention. However, the present invention may also be implemented in other ways different from those described herein. Therefore, the protection scope of the present invention is not limited to the specific embodiments disclosed below.

[0054] like Figure 1 and Figure 2 As shown, this embodiment provides a ferromagnetic cavity force-magnetic coupling modeling method, the method comprising:

[0055] S1. Construct a finite element magnetoelastic coupling model framework, including a model input module, a magnetoelastic coupling parameter automatic identification module, a residual magnetism inversion module, a magnetoelastic coupling calculation module and a model output module; wherein the model input module, the magnetoelastic coupling parameter automatic identification module and the residual magnetism inversion module all receive external data and input the results into the magnetoelastic coupling calculation module, the magnetoelastic coupling calculation module inputs the results into the model output module, and the model output module outputs the magnetoelastic coupling model of the ferromagnetic cavity.

[0056] S2. Measure the geometric parameters of the ferromagnetic cavity, input the measured geometric parameters into the model input module, and establish a geometric model of equal scale in COMSOL Multiphysics finite element software.

[0057] S3. Measure the saturation hysteresis loop of the material sample in the ferromagnetic cavity. The measurement method includes: cutting the ferromagnetic material into cubes with an edge length of 1 mm and sticking them to the center of the test rod, adjusting the scanning step size and other parameters of the superconducting quantum interference device so that the maximum magnetic field exceeds the magnetic field size of the material saturation, and using the superconducting quantum interference device to measure the saturation hysteresis loop of the sample.

[0058] The measured saturation hysteresis loop results are introduced into the Jiles-Atherton parameter identification algorithm. The Jiles-Atherton parameter identification algorithm modifies the traditional Jiles-Atherton equation according to the format of the fourth-order Runge-Kutta method. The characteristic points of the test results are selected, and then the particle swarm algorithm is used to adjust the M of the ferromagnetic material. s , α, k, c, b five Jiles-Atherton parameters, to achieve the matching of the calculation results with the characteristic point test results, the specific steps include:

[0059] The magnetization formula of the Jiles-Atherton model is:

[0060]

[0061] Among them, M s is the saturation magnetization, k is the internal coupling coefficient, c is the reversibility coefficient, and M is the magnetization of the material under an external magnetic field H.

[0062] M irr is the irreversible magnetization intensity, which can be expressed as:

[0063]

[0064] Among them, H eff =H-αM is the effective magnetic field, α is the shape factor, which represents the microscopic properties of the material, and b is the elastic coefficient of the magnetic domain wall.

[0065] The total magnetization consists of a reversible and an irreversible part:

[0066] M=M irr +c·(M s -M irr );

[0067] The particle swarm optimization algorithm is used to adjust M s , α, k, c, b are used to make the model calculation results match the experimentally measured hysteresis loop. The objective function E is defined as:

[0068]

[0069] Among them, M exp (H i ) represents the experimentally measured magnetization intensity, M calc (H i ) represents the magnetization intensity calculated by the model.

[0070] The objective function E is to minimize the experimental data M exp (H i ) and the calculated result M calc (H i) to adjust the model parameters so that the calculated results of the model match the experimental data.

[0071] The final five Jiles-Atherton parameters are input into the model input module.

[0072] S4. Measure the Young's modulus of the material sample of the ferromagnetic cavity. The measuring method includes: cutting the ferromagnetic material into long strips of 1 cm×1 cm×10 cm, using a tensile testing machine to stretch both ends of the 10 cm edge of the sample to obtain the Young's modulus of the sample.

[0073] Enter the measured Young's modulus into the model input module.

[0074] S5. Perform a residual magnetism measurement test on the ferromagnetic cavity in an environment where the background magnetic field is 0. Arrange a three-phase fluxgate sensor at a specified point to measure the fixed magnetic field generated by the residual magnetism of the ferromagnetic body, and obtain the fixed magnetic field distribution of the ferromagnetic cavity. The residual magnetism inversion module uses the "conversion matrix method" to measure the residual magnetism of the ferromagnetic body. Specifically, the relationship between the magnetization intensity matrix M0 and the magnetic induction intensity matrix B0 measured by the sensor can be expressed as:

[0075] M0=AB0;

[0076] Where A is the transformation matrix, expressed as:

[0077] A=B -1 M;

[0078] Among them, B -1 is the magnetic field calculated in COMSOL Multiphysics finite element software, and M is the magnetization intensity matrix obtained using random function in COMSOL Multiphysics finite element software.

[0079] The results of the remanent magnetization measurement test are input into the remanent magnetization inversion module.

[0080] The background magnetic field is made zero by adjusting the compensation current.

[0081] S6. Conduct a magnetomechanical coupling test on the ferromagnetic cavity in an environment where the background magnetic field is 0. The test steps include: after adjusting the background magnetic field to 0, pressurizing the ferromagnetic cavity and recording the spatial magnetic field every 2MPa until it reaches 10Mpa; then depressurizing the ferromagnetic cavity and recording the spatial magnetic field every 2MPa until the pressure reaches 0.

[0082] The results of the magnetomechanical coupling test are input into the magnetomechanical coupling parameter automatic identification module. The equivalent field effect parameter expression to be optimized is as follows:

[0083]

[0084] Among them, α σx , α σy , α σz is the equivalent field effect parameter in the x, y, and z directions, α is the field effect parameter of the ferromagnetic material, σ x , σ y , σ z is the stress magnitude in the x, y, and z directions, calculated by COMSOL Multiphysics finite element software, μ0 is the vacuum permeability, γ 1x ~γ 1z ’ are the six parameters to be optimized. Similarly, the particle swarm algorithm is used to optimize the six parameters to achieve matching with the magnetomechanical coupling test data.

[0085] S7, the magnetomechanical coupling calculation module outputs the magnetic field of a specified point under specified working conditions, the magnetic field cloud map of a specified surface, and the magnetization intensity distribution cloud map of the ferromagnetic cavity based on the data from the model input module, the magnetomechanical coupling parameter automatic identification module, and the residual magnetism inversion module.

[0086] The specified working condition is determined by the pressure function y=press(t) and the time variable t, the specified background magnetic field is input through the model, and the specified point and the specified surface are determined by the space coordinates.

[0087] S8. The model output module generates and outputs the ferromagnetic cavity magnetomechanical coupling model according to the input results of the magnetomechanical coupling calculation module.

[0088] In another embodiment of the present invention, Figure 3 As shown in the figure, a magnetomechanical coupling model of an oxygen cylinder is established in the Comosl Multiphysics finite element simulation software. The oxygen cylinder is made of 37Mn. A tensile test is performed on the 37Mn material and its saturation hysteresis loop is measured. Parameter identification is performed based on the test data, and the Jiles-Atherton parameters of 37Mn are obtained as follows: M s =1.8×10 6 A / m, a=100801.22A / m, k=30.41A / m, c=0.99, α=1.12×10 -3 At room temperature, the relationship between its Young's modulus and temperature is: T = 2.11 × 10 11 +3.57×10 7 *T – 1.06×10 5 *T 2 .

[0089] The residual magnetism measurement test and force-magnetic coupling test were carried out on the oxygen cylinder, and the fixed magnetic field distribution on the keel line was obtained as follows: Figure 4 As shown, the magnetomechanical coupling parameter of 37Mn is: γ 1x(0) = 5.1 × 10 -17 m 2 / A 2 , γ 1x ’ (0) = -1.5 × 10 -24 m 2 Pa / A 2 , γ 1y (0) = 5 × 10 -17 m 2 / A 2 , γ 1y ’ (0) = -1.7 × 10 -24 m 2 Pa / A 2 , γ 1z (0) = 5.1 × 10 -17 m 2 / A 2 , γ 1z ’ (0) = -2.3 × 10 -24 m 2 Pa / A 2 .

[0090] The above parameters are brought into the finite element model, and the simulation results show that the magnetic field data on the keel line are compared with the experimental data when the pressure is 8.1 MPa. Figure 5 As shown, after calculation, the model calculation accuracy exceeds 90%.

[0091] In summary, the present invention proposes a ferromagnetic cavity force-magnetic coupling modeling method, comprising:

[0092] S1. Construct a finite element magnetoelastic coupling model framework, including a model input module, a magnetoelastic coupling parameter automatic identification module, a residual magnetism inversion module, a magnetoelastic coupling calculation module and a model output module; wherein the model input module, the magnetoelastic coupling parameter automatic identification module and the residual magnetism inversion module all receive external data and input the results into the magnetoelastic coupling calculation module, the magnetoelastic coupling calculation module inputs the results into the model output module, and the model output module outputs the magnetoelastic coupling model of the ferromagnetic cavity.

[0093] S2. Measure the geometric parameters of the ferromagnetic cavity, input the measured geometric parameters into the model input module, and establish a geometric model of equal scale in COMSOL Multiphysics finite element software.

[0094] S3. Measure the saturation hysteresis loop of the material sample of the ferromagnetic cavity, bring the measured saturation hysteresis loop result into the Jiles-Atherton parameter identification algorithm, and input the five finally obtained Jiles-Atherton parameters into the model input module.

[0095] S4. Measure the Young's modulus of the material sample of the ferromagnetic cavity, and input the measured Young's modulus into the model input module.

[0096] S5. Perform a remanent magnetism measurement test on the ferromagnetic cavity in an environment where the background magnetic field is 0, and input the result of the remanent magnetism measurement test into the remanent magnetism inversion module.

[0097] S6. Perform a magnetomechanical coupling test on the ferromagnetic cavity in an environment where the background magnetic field is 0, and input the magnetomechanical coupling test results into the magnetomechanical coupling parameter automatic identification module.

[0098] S7, the magnetomechanical coupling calculation module outputs the magnetic field of a specified point under specified working conditions, the magnetic field cloud map of a specified surface, and the magnetization intensity distribution cloud map of the ferromagnetic cavity based on the data from the model input module, the magnetomechanical coupling parameter automatic identification module, and the residual magnetism inversion module.

[0099] The specified working condition is determined by the pressure function y=press(t) and the time variable t, the specified background magnetic field is input through the model, and the specified point and the specified surface are determined by the space coordinates.

[0100] S8. The model output module generates and outputs the ferromagnetic cavity magnetomechanical coupling model according to the input results of the magnetomechanical coupling calculation module.

[0101] The steps in the present invention can be adjusted in order, combined or deleted according to actual needs.

[0102] The units in the device of the present invention can be combined, divided and deleted according to actual needs.

[0103] Although the present invention is disclosed in detail with reference to the accompanying drawings, it should be understood that these descriptions are merely exemplary and are not intended to limit the application of the present invention. The scope of protection of the present invention is defined by the appended claims and may include various modifications, alterations and equivalents made to the invention without departing from the scope and spirit of the present invention.

Claims

1. A ferromagnetic cavity force-magnetic coupling modeling method, characterized in that: The method comprises: S1. Construct a finite element magnetoelastic coupling model framework, including a model input module, a magnetoelastic coupling parameter automatic identification module, a residual magnetism inversion module, a magnetoelastic coupling calculation module and a model output module; wherein the model input module, the magnetoelastic coupling parameter automatic identification module and the residual magnetism inversion module all receive external data and input the results into the magnetoelastic coupling calculation module, the magnetoelastic coupling calculation module inputs the results into the model output module, and the model output module outputs the magnetoelastic coupling model of the ferromagnetic cavity; S2, measuring the geometric parameters of the ferromagnetic cavity, inputting the measured geometric parameters into the model input module, and establishing a geometric model of equal scale in COMSOL Multiphysics finite element software; S3, measuring the saturation hysteresis loop of the material sample of the ferromagnetic cavity, bringing the measured saturation hysteresis loop result into the Jiles-Atherton parameter identification algorithm, and inputting the five finally obtained Jiles-Atherton parameters into the model input module; S4, measuring the Young's modulus of the material sample of the ferromagnetic cavity, and inputting the measured Young's modulus into the model input module; S5, performing a residual magnetism measurement test on the ferromagnetic cavity in an environment where the background magnetic field is 0, and inputting the result of the residual magnetism measurement test into a residual magnetism inversion module; S6. Performing a force-magnetic coupling test on the ferromagnetic cavity in an environment where the background magnetic field is 0, and inputting the force-magnetic coupling test result into the force-magnetic coupling parameter automatic identification module; S7, the magnetomechanical coupling calculation module outputs the magnetic field of a specified point under a specified working condition, the magnetic field cloud map of a specified surface, and the magnetization intensity distribution cloud map of the ferromagnetic cavity according to the data from the model input module, the magnetomechanical coupling parameter automatic identification module, and the residual magnetism inversion module; The specified working condition is determined by the pressure function y=press(t) and the time variable t, the specified background magnetic field is input through the model, and the specified point and the specified surface are determined by the spatial coordinates; S8. The model output module generates and outputs the ferromagnetic cavity magnetomechanical coupling model according to the input results of the magnetomechanical coupling calculation module.

2. The ferromagnetic cavity force-magnetic coupling modeling method according to claim 1, characterized in that: The method for measuring the saturation hysteresis loop in step S3 includes: cutting the ferromagnetic material into cubes with an edge length of 1 mm and sticking them to the center of the test rod, adjusting parameters such as the scanning step length of the superconducting quantum interference instrument so that the maximum magnetic field exceeds the magnetic field size of the material saturation, and using the superconducting quantum interference instrument to measure the saturation hysteresis loop of the sample.

3. The ferromagnetic cavity force-magnetic coupling modeling method according to claim 2, characterized in that: The Jiles-Atherton parameter identification algorithm specifically includes: The magnetization formula of the Jiles-Atherton model is: Among them, M s is the saturation magnetization, k is the internal coupling coefficient, c is the reversibility coefficient, and M is the magnetization of the material under an external magnetic field H; M irr is the irreversible magnetization intensity, which can be expressed as: Among them, H eff =H-αM is the effective magnetic field, α is the shape factor, which represents the microscopic properties of the material, and b is the elastic coefficient of the magnetic domain wall; The total magnetization consists of a reversible and an irreversible part: M=M irr +c·(M s -M irr ); The particle swarm optimization algorithm is used to adjust M s , α, k, c, b are used to make the model calculation results match the experimentally measured hysteresis loop. The objective function E is defined as: Among them, M exp (H i ) represents the experimentally measured magnetization intensity, M calc (H i ) represents the magnetization intensity calculated by the model; The objective function E is to minimize the experimental data M exp (H i ) and the calculated result M calc (H i ) to adjust the model parameters so that the calculated results of the model match the experimental data.

4. The ferromagnetic cavity force-magnetic coupling modeling method according to claim 1, characterized in that: The method for measuring the Young's modulus in step S4 includes: cutting the ferromagnetic material into long strip samples of 1 cm×1 cm×10 cm, using a tensile testing machine to stretch both ends of the 10 cm edge of the sample to obtain the Young's modulus of the sample.

5. The ferromagnetic cavity force-magnetic coupling modeling method according to claim 1, characterized in that: The residual magnetism measurement test in step S5 includes: adjusting the compensation current to make the background magnetic field zero, arranging a three-phase fluxgate sensor at a specified point, measuring the fixed magnetic field generated by the residual magnetism of the ferromagnetic body, and obtaining the fixed magnetic field distribution of the ferromagnetic cavity.

6. The ferromagnetic cavity force-magnetic coupling modeling method according to claim 5, characterized in that: The remanent magnetism inversion module uses the "conversion matrix method" to measure the remanent magnetism of the ferromagnetic body. The relationship between the magnetization intensity matrix M0 and the magnetic induction intensity matrix B0 measured by the sensor can be expressed as: M0=AB0; Where A is the transformation matrix, expressed as: A=B -1 M; Among them, B -1 is the magnetic field calculated in COMSOL Multiphysics finite element software, and M is the magnetization intensity matrix obtained using random function in COMSOL Multiphysics finite element software.

7. The ferromagnetic cavity force-magnetic coupling modeling method according to claim 1, characterized in that: The steps of the magnetomechanical coupling test in step S5 include: after adjusting the background magnetic field to 0, pressurizing the ferromagnetic cavity and recording the spatial magnetic field every 2MPa until it reaches 10MPa; then depressurizing the ferromagnetic cavity and recording the spatial magnetic field every 2MPa until the pressure reaches 0.

8. The ferromagnetic cavity force-magnetic coupling modeling method according to claim 7, characterized in that: In the automatic identification module of magnetomechanical coupling parameters, the magnetomechanical coupling test results are input, and the equivalent field effect parameter expression to be optimized is as follows: Among them, α σx , α σy , α σz is the equivalent field effect parameter in the x, y, and z directions, α is the field effect parameter of the ferromagnetic material, σ x , σ y , σ z is the stress magnitude in the x, y, and z directions, calculated by COMSOL Multiphysics finite element software, μ0 is the vacuum permeability, γ 1x ~γ 1z ’ are the six parameters to be optimized. Similarly, the particle swarm algorithm is used to optimize the six parameters to achieve matching with the magnetomechanical coupling test data.

Images