CN118378576A - A three-dimensional time-domain hysteresis loop drawing method and system based on hysteresis model - Google Patents

Abstract

The invention discloses a three-dimensional time-domain hysteresis loop drawing method and system based on a hysteresis model. The method of the invention comprises the following steps: constructing a J-A hysteresis model corresponding to a target iron core; obtaining dynamic stray current values of the target iron core at different times, calculating magnetic field strength according to the current value, and superimposing the dynamic stray magnetic field strength and the normal AC magnetic field strength as the total input of the J-A hysteresis model; using a backtracking search optimization algorithm to perform parameter identification on the J-A hysteresis model, and determining the target magnetic induction strength based on the target magnetic field strength and the actual magnetization strength; performing curve fitting on the final magnetic field strength and the magnetic induction strength, and obtaining the three-dimensional hysteresis loop of the target iron core; the invention superimposes the magnetic field strength generated by the dynamic stray current and the normal AC magnetic field strength as the total input of the model, draws a three-dimensional hysteresis loop considering the continuous interference of the dynamic stray current, and demonstrates the time-varying and non-periodicity of the dynamic stray current in affecting the hysteresis loop.

Description

A three-dimensional time-domain hysteresis loop drawing method and system based on hysteresis model

Technical Field

The invention relates to the technical field of iron core magnetization characteristics, and in particular to a three-dimensional time-domain hysteresis loop depiction method and system based on a hysteresis model.

Background technique

With the rapid development of urban rail transit, the dynamic stray current generated by the incomplete insulation of the running rails to the ground is gradually increasing, and its impact on the transformer is increasing. Dynamic stray current, that is, stray current that changes dynamically with time, will cause a series of problems such as increased vibration noise, increased losses, and increased harmonics of the excitation current after flowing into the transformer through the neutral point. Therefore, it is very important to explore the magnetic characteristics of the transformer core under the interference of dynamic stray current.

Hysteresis loop parameter identification is one of the key means to study the magnetic properties of the core. Existing research mainly focuses on parameter identification of hysteresis loops under normal working conditions and changes in material properties, so as to make the description of the hysteresis loop more accurate. The traditional hysteresis loop is used to represent the hysteresis phenomenon of magnetic materials under periodically changing magnetic field strength, indicating that the relationship between the magnetic induction intensity B and the magnetic field intensity H of magnetic materials during repeated magnetization has obvious periodicity. However, the dynamic stray current has a dynamic change in its current value over time, and the magnetic field intensity generated also changes dynamically over time, so the impact on the magnetic properties of the core changes with time and is no longer periodic. Since the traditional hysteresis loop can only represent the changing relationship between H and B on a two-dimensional level, and ignores the changing relationship between H and B in the time domain, the traditional hysteresis loop can no longer accurately describe the changing law of H and B in the time domain, that is, it cannot accurately describe the hysteresis characteristics of magnetic materials under the influence of dynamic stray currents.

Summary of the invention

In view of the deficiencies in the prior art, the purpose of the present invention is to provide a method and system that can accurately describe the changes in the hysteresis characteristics of magnetic materials under the influence of dynamic stray currents, so as to provide assistance for subsequent research on hysteresis losses and the operating temperature of magnetic materials.

To achieve the above object, the present invention provides the following technical solution: a three-dimensional time-domain hysteresis loop drawing method based on a hysteresis model, comprising:

Construct the J-A hysteresis model corresponding to the target core;

The current values of the dynamic stray current of the target core at multiple different times are obtained, and the magnetic field strength generated by the dynamic stray current is calculated based on the obtained current values. The target magnetic field strength is obtained by superimposing the magnetic field strength generated by the dynamic stray current with the AC magnetic field strength generated by the normal power frequency AC current inside the target core. As the total input of JA hysteresis model;

The backtracking search optimization algorithm is used to input the target magnetic field strength The JA hysteresis model is then identified, and the model parameters of the JA hysteresis model obtained after parameter identification are substituted into the JA hysteresis model. The corresponding actual magnetization intensity , using the target magnetic field strength and the actual magnetization Obtain the target magnetic induction intensity B(t);

The target magnetic field intensity and the target magnetic induction intensity are curve fitted to obtain the three-dimensional hysteresis loop of the target core.

Furthermore, the J-A hysteresis model corresponding to the target core is expressed as:

(1);

In the formula, Indicates the actual magnetization Target magnetic field strength The differential of represents a differential operator; The coefficient introduced to prevent non-physical solutions; is the reversible magnetic susceptibility; is the irreversible loss coefficient; is the directional coefficient. When the rate of change of the target magnetic field strength dH(t)/dt>0, , when dH(t)/dt<0, ; t represents time; is the magnetization intensity without hysteresis effect, is the domain wall interaction coefficient.

Furthermore, the hysteresis-free magnetization Expressed as:

(2);

In the formula, is the saturation magnetization; is the magnetization intensity shape factor without hysteresis effect; represents the effective magnetic field strength, where is the domain wall interaction coefficient.

Furthermore, the specific process of calculating the magnetic field strength generated by the dynamic stray current based on the acquired current value is as follows: Solving the magnetic field strength generated by the dynamic stray current based on Ampere's loop law :

(3);

In the formula, is the number of turns of the winding coil; is the dynamic stray current; is the effective magnetic path length.

Furthermore, the magnetic field strength generated by the dynamic stray current The AC magnetic field strength generated by the normal power frequency AC current inside the target core The superimposed target magnetic field strength The total input of the JA hysteresis model is expressed by equation (4):

(4).

Furthermore, the backtracking search optimization algorithm is used to identify the parameters of the JA hysteresis model, and the model parameters obtained after parameter identification are substituted into the JA hysteresis model, and on this basis, the target magnetic field intensity is used to The corresponding actual magnetization intensity , using the target magnetic field strength and the actual magnetization The specific process of obtaining the target magnetic induction intensity B(t) is:

Initialize the parameters of the backtracking search optimization algorithm;

Set the population size and number of iterations of the backtracking search optimization algorithm. At the same time, give the step size factor γ and the random perturbation factor r, and initialize the individuals that meet the constraints within the constraints. Each individual is , , , , There are five parameters in the JA hysteresis model. , , , , The value range of constitutes the constraint condition;

Iteratively update individuals in the population of the backtracking search optimization algorithm;

By changing the individual inclusion , , , , Parameters, to achieve iterative update of individuals within the population;

(5);

In the formula, The number of individuals in a population , , , , The step size of the parameter change, where is the upper bound of the constraint condition, is the lower bound of the constraint; For the Generation population; For the Generation population; It means to generate an array of random numbers uniformly distributed between 0 and 1; Indicates obtaining the number of rows and columns of the target matrix;

For individuals exceeding the constraints , , , , The parameters are controlled by upper and lower limits according to equations (6) and (7):

Exceeding the limit: (6);

Exceeded lower limit: (7);

In the formula, Represents an individual The mth one contained in , , , , parameter; The mth upper limit , , , , parameter; The mth lower limit , , , , parameter;

Substitute all updated individuals in the population into the JA hysteresis model to obtain the actual magnetization intensity at each moment , determine the target magnetic induction intensity B(t) at each moment based on the target magnetic field intensity at each moment and the actual magnetization intensity at each moment;

The fitness function is used to calculate the fitness value of each individual in the population after each iterative update. The fitness function F is expressed as:

(8);

In the formula, is the target magnetic flux density data calculated by the JA hysteresis model. Values; is the jth value in the reference magnetic induction intensity data; is the number of data points;

Determine whether the termination condition is reached. When the termination condition is reached, select the individual with the smallest fitness value as the optimal output. Otherwise, determine whether the number of iterative updates is the specified value. When the number of iterative updates is the specified value, reinitialize. Otherwise, start the next iterative update of individuals in the population.

Furthermore, the target magnetic field intensity and the target magnetic induction intensity are curve fitted to obtain the three-dimensional hysteresis loop of the target iron core. The specific process is as follows: construct a three-dimensional rectangular coordinate system with the x-axis, y-axis, and z-axis as coordinate axes, and transform the target magnetic field intensity and target magnetic induction intensity The three-dimensional hysteresis loop of the target core is drawn by taking the x-axis and y-axis data of the three-dimensional rectangular coordinate system respectively and the time t as the z-axis data of the three-dimensional rectangular coordinate system.

A three-dimensional time-domain hysteresis loop drawing system based on a hysteresis model, comprising:

A construction module is used to construct a J-A hysteresis model corresponding to the target core;

A calculation module is used to obtain the current values of the dynamic stray current of the target iron core at multiple different times, calculate the magnetic field strength generated by the dynamic stray current according to the obtained current value, and superimpose the magnetic field strength generated by the dynamic stray current with the AC magnetic field strength generated by the normal power frequency AC current inside the target iron core to obtain the target magnetic field strength as the total input of the J-A hysteresis model;

An identification module is used to perform parameter identification on the J-A hysteresis model after inputting the target magnetic field strength by using a backtracking search optimization algorithm, substitute the model parameters of the J-A hysteresis model obtained after the parameter identification into the J-A hysteresis model, and obtain the corresponding actual magnetization intensity based on the target magnetic field strength on this basis, and obtain the target magnetic induction intensity by using the target magnetic field strength and the actual magnetization intensity;

The fitting module is used to perform curve fitting on the target magnetic field intensity and the target magnetic induction intensity to obtain the three-dimensional hysteresis loop of the target iron core.

Furthermore, the J-A hysteresis model corresponding to the target core is expressed as:

;

In the formula, Indicates the actual magnetization Target magnetic field strength The differential of represents a differential operator; The coefficient introduced to prevent non-physical solutions; is the reversible magnetic susceptibility; is the irreversible loss coefficient; is the directional coefficient. When the rate of change of the target magnetic field strength dH(t)/dt>0, , when dH(t)/dt<0, ; t represents time; is the magnetization intensity without hysteresis effect, is the domain wall interaction coefficient;

The hysteresis-free magnetization Expressed as:

;

In the formula, is the saturation magnetization; is the magnetization intensity shape factor without hysteresis effect; represents the effective magnetic field strength, where is the domain wall interaction coefficient.

A computer-readable storage medium stores a computer program, which, when executed by a processor, implements a three-dimensional time-domain hysteresis loop drawing method based on a hysteresis model.

Compared with the existing technology, the present invention has the following beneficial effects: the present invention superimposes the magnetic field intensity generated by the dynamic stray current and the normal AC magnetic field intensity as the total input of the model, draws a three-dimensional hysteresis loop considering the continuous interference of the dynamic stray current, and shows the time-varying and non-periodicity of the influence of the dynamic stray current on the hysteresis loop. At the same time, through the three-dimensional hysteresis loop, the change of the hysteresis loop in each normal AC magnetic field intensity change cycle can be observed, and the hysteresis loss and temperature change of the magnetic material when working can be evaluated.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG1 is a flow chart of the present invention for parameter identification of the J-A hysteresis model under dynamic stray current interference.

FIG. 2 is a flow chart of parameter identification of the intelligent algorithm of the present invention.

FIG3 is a diagram showing the measured neutral point current waveform of a substation according to the present invention.

FIG4 is a waveform diagram of the dynamic stray current magnetic field intensity and the AC magnetic field intensity superimposed on each other according to the present invention.

FIG5 is a three-dimensional time domain comparison diagram of the J-A hysteresis model parameter identification of the present invention.

FIG6 is a schematic diagram of the system structure of the present invention.

Detailed ways

As shown in FIG1 , the present invention provides a technical solution: a three-dimensional hysteresis loop drawing method based on a hysteresis model, comprising the following steps:

Step S1: Construct the J-A hysteresis model corresponding to the target core:

(1);

In the formula, Indicates the actual magnetization Target magnetic field strength The differential of represents a differential operator; The coefficient introduced to prevent non-physical solutions; is the reversible magnetic susceptibility; is the irreversible loss coefficient; is the directional coefficient. When the rate of change of the target magnetic field strength dH(t)/dt>0, , when dH(t)/dt<0, ; t represents time; is the domain wall interaction coefficient; is the magnetization intensity without hysteresis effect, expressed as:

(2);

In the formula, is the saturation magnetization; is the magnetization intensity shape factor without hysteresis effect; represents the effective magnetic field strength, where is the domain wall interaction coefficient.

Step S2: Obtain the current values of the dynamic stray current of the target core at multiple different times, as shown in Figure 3, calculate the magnetic field strength generated by the dynamic stray current based on the obtained current value, and superimpose the magnetic field strength generated by the dynamic stray current with the AC magnetic field strength generated by the normal power frequency AC current inside the target core to obtain the target magnetic field strength As the total input of the J-A hysteresis model.

Step S3: Use the backtracking search optimization algorithm to input the target magnetic field strength The JA hysteresis model is then identified, and the model parameters of the JA hysteresis model obtained after parameter identification are substituted into the JA hysteresis model. The corresponding actual magnetization intensity , using the target magnetic field strength and the actual magnetization The target magnetic induction intensity B(t) is obtained, as shown in FIG2 .

Step S4: performing curve fitting on the target magnetic field intensity and the target magnetic induction intensity to obtain a three-dimensional hysteresis loop of the target core.

The specific process of calculating the magnetic field strength generated by the dynamic stray current based on the acquired current value is as follows: Based on Ampere's loop law, the magnetic field strength generated by the dynamic stray current is solved. :

(3);

In the formula, is the number of turns of the winding coil; is the dynamic stray current; is the effective magnetic path length.

Among them, the magnetic field strength generated by the dynamic stray current is The AC magnetic field strength generated by the normal power frequency AC current The superimposed target magnetic field strength The total input of the JA hysteresis model is shown in formula (4), and the waveform of the calculated target magnetic field strength is shown in Figure 4;

(4);

The specific process of step S3 is:

Step S3.1: Initialize the parameters of the backtracking search optimization algorithm;

Set the population size and number of iterations of the backtracking search optimization algorithm. At the same time, give the step size factor γ and the random perturbation factor r, and initialize the individuals that meet the constraints within the constraints. Each individual is , , k, a, c, a total of five parameters in the JA hysteresis model, , The value ranges of k, a, and c constitute the constraints.

Step S3.2: iteratively update the individuals in the population of the backtracking search optimization algorithm;

By changing the individual inclusion , , , , Parameters realize iterative updates of individuals within the population;

(5);

In the formula, The number of individuals in a population , , , , The step size of the parameter change, where is the upper bound of the constraint condition, is the lower bound of the constraint; For the Generation population; For the Generation population; It means to generate an array of random numbers uniformly distributed between 0 and 1; Indicates getting the number of rows and columns of the target matrix.

Step S3.3: For individuals that exceed the constraints , , , , The parameters are controlled by upper and lower limits according to equations (6) and (7):

Exceeding the limit: (6);

Exceeded lower limit: (7);

In the formula, Represents an individual The mth one contained in , , , , parameter; The mth upper limit , , , , parameter; The mth lower limit , , , , parameter.

Step S3.4: Substitute all the updated individuals in the population into the J-A hysteresis model to obtain the actual magnetization intensity M at each moment, and determine the target magnetic induction intensity B(t) at each moment based on the target magnetic field intensity at each moment and the actual magnetization intensity at each moment;

The fitness function is used to calculate the fitness value of each individual in the population after each iteration. The smaller the fitness value, the better the individual. The fitness function F is expressed as:

(8);

In the formula, is the target magnetic flux density data calculated by the JA hysteresis model. Values; is the jth value in the reference magnetic induction intensity data; is the number of data points; the smaller the value of the fitness function F is, the closer the JA hysteresis model parameters identified by the backtracking search optimization algorithm are to the true values.

Step S3.5: Determine whether the termination condition is met. When the termination condition is met, select the individual with the smallest fitness value as the optimal output. Otherwise, determine whether the number of iterative updates is an integer multiple of 5. If the number of iterative updates is an integer multiple of 5, reinitialize. Otherwise, start the next iterative update of individuals in the population.

The specific process of step S4 is: construct a three-dimensional rectangular coordinate system with x-axis, y-axis and z-axis as coordinate axes, and convert the target magnetic field intensity and target magnetic induction intensity The three-dimensional hysteresis loop of the target core is drawn by taking the x-axis and y-axis data of the three-dimensional rectangular coordinate system respectively and the time t as the z-axis data of the three-dimensional rectangular coordinate system.

The three-dimensional hysteresis loop of the target core obtained by using the backtracking search optimization algorithm to identify the parameters of the J-A hysteresis model is shown in Figure 5. It can be seen from Figure 5 that the three-dimensional hysteresis loop of the target core obtained by using the backtracking search optimization algorithm to identify the parameters of the J-A hysteresis model is close to the benchmark value, and this scheme can achieve better results.

As shown in FIG6 , a three-dimensional time-domain hysteresis loop drawing system based on a hysteresis model includes:

A construction module is used to construct a J-A hysteresis model corresponding to the target core;

A calculation module is used to obtain the current values of the dynamic stray current of the target iron core at multiple different times, calculate the magnetic field strength generated by the dynamic stray current according to the obtained current value, and superimpose the magnetic field strength generated by the dynamic stray current with the AC magnetic field strength generated by the normal power frequency AC current inside the target iron core to obtain the target magnetic field strength as the total input of the J-A hysteresis model;

An identification module is used to perform parameter identification on the J-A hysteresis model after inputting the target magnetic field strength by using a backtracking search optimization algorithm, substitute the model parameters of the J-A hysteresis model obtained after the parameter identification into the J-A hysteresis model, and obtain the corresponding actual magnetization intensity based on the target magnetic field strength on this basis, and obtain the target magnetic induction intensity by using the target magnetic field strength and the actual magnetization intensity;

The fitting module is used to perform curve fitting on the target magnetic field intensity and the target magnetic induction intensity to obtain the three-dimensional hysteresis loop of the target iron core.

Among them, the J-A hysteresis model corresponding to the target core is expressed as:

;

In the formula, Indicates the actual magnetization Target magnetic field strength The differential of represents a differential operator; The coefficient introduced to prevent non-physical solutions; is the reversible magnetic susceptibility; is the irreversible loss coefficient; is the directional coefficient. When the rate of change of the target magnetic field strength dH(t)/dt>0, , when dH(t)/dt<0, ; t represents time; is the magnetization intensity without hysteresis effect, is the domain wall interaction coefficient.

Wherein, the hysteresis-free magnetization intensity Expressed as:

;

In the formula, is the saturation magnetization; is the magnetization intensity shape factor without hysteresis effect; represents the effective magnetic field strength, where is the domain wall interaction coefficient.

A computer-readable storage medium stores a computer program, which, when executed by a processor, implements a three-dimensional time-domain hysteresis loop drawing method based on a hysteresis model.

Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present application may adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes. The schemes in the embodiments of the present application may be implemented in various computer languages, for example, object-oriented programming language Java and literal scripting language JavaScript, etc.

The present application is described with reference to the flowcharts and/or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and/or box in the flowchart and/or block diagram, as well as the combination of the processes and/or boxes in the flowchart and/or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate a device for implementing the functions specified in one process or multiple processes in the flowchart and/or one box or multiple boxes in the block diagram.

These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce a manufactured product including an instruction device that implements the functions specified in one or more processes in the flowchart and/or one or more boxes in the block diagram.

These computer program instructions may also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, whereby the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in one or more processes in the flowchart and/or one or more boxes in the block diagram.

Although the preferred embodiments of the present application have been described, those skilled in the art may make other changes and modifications to these embodiments once they have learned the basic creative concept. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments and all changes and modifications falling within the scope of the present application.

Obviously, those skilled in the art can make various changes and modifications to the present application without departing from the spirit and scope of the present application. Thus, if these modifications and variations of the present application fall within the scope of the claims of the present application and their equivalents, the present application is also intended to include these modifications and variations.

Claims (10)

1. A three-dimensional time-domain hysteresis loop drawing method based on a hysteresis model, characterized by comprising: Construct the J-A hysteresis model corresponding to the target core; The current values of the dynamic stray current of the target core at multiple different times are obtained, and the magnetic field strength generated by the dynamic stray current is calculated based on the obtained current values. The target magnetic field strength is obtained by superimposing the magnetic field strength generated by the dynamic stray current with the AC magnetic field strength generated by the normal power frequency AC current inside the target core. As the total input of JA hysteresis model; The backtracking search optimization algorithm is used to input the target magnetic field strength The JA hysteresis model is then identified, and the model parameters of the JA hysteresis model obtained after parameter identification are substituted into the JA hysteresis model. The corresponding actual magnetization intensity , using the target magnetic field strength and the actual magnetization Obtain the target magnetic induction intensity B(t); The target magnetic field intensity and the target magnetic induction intensity are curve fitted to obtain the three-dimensional hysteresis loop of the target core. 2. A three-dimensional time-domain hysteresis loop drawing method based on a hysteresis model according to claim 1, characterized in that: the J-A hysteresis model corresponding to the target core is expressed as: (1); In the formula, Indicates the actual magnetization Target magnetic field strength The differential of represents a differential operator; The coefficient introduced to prevent non-physical solutions; is the reversible magnetic susceptibility; is the irreversible loss coefficient; is the directional coefficient. When the rate of change of the target magnetic field strength dH(t)/dt>0, , when dH(t)/dt<0, ; t represents time; is the magnetization intensity without hysteresis effect; is the domain wall interaction coefficient. 3. A three-dimensional time-domain hysteresis loop drawing method based on a hysteresis model according to claim 2, characterized in that: the hysteresis-free magnetization intensity Expressed as: (2); In the formula, is the saturation magnetization; is the magnetization intensity shape factor without hysteresis effect; represents the effective magnetic field strength, where is the domain wall interaction coefficient. 4. A three-dimensional time-domain hysteresis loop drawing method based on a hysteresis model according to claim 3, characterized in that: the specific process of calculating the magnetic field strength generated by the dynamic stray current according to the acquired current value is: solving the magnetic field strength generated by the dynamic stray current based on Ampere's loop law : (3); In the formula, is the number of turns of the winding coil; is the dynamic stray current; is the effective magnetic path length. 5. A three-dimensional time-domain hysteresis loop drawing method based on a hysteresis model according to claim 4, characterized in that: the magnetic field strength generated by the dynamic stray current The AC magnetic field strength generated by the normal power frequency AC current inside the target core The superimposed target magnetic field strength The total input of the JA hysteresis model is expressed by equation (4): (4). 6. A three-dimensional time-domain hysteresis loop drawing method based on a hysteresis model according to claim 5, characterized in that: a backtracking search optimization algorithm is used to identify the parameters of the JA hysteresis model, the model parameters obtained after the parameter identification are substituted into the JA hysteresis model, and on this basis, the target magnetic field intensity is used to determine the parameters of the JA hysteresis model. The corresponding actual magnetization intensity , using the target magnetic field strength and the actual magnetization The specific process of obtaining the target magnetic induction intensity B(t) is: Initialize the parameters of the backtracking search optimization algorithm; Set the population size and number of iterations of the backtracking search optimization algorithm. At the same time, give the step size factor γ and the random perturbation factor r, and initialize the individuals that meet the constraints within the constraints. Each individual is , , , , There are five parameters in the JA hysteresis model. , , , , The value range of constitutes the constraint condition; Iteratively update individuals in the population of the backtracking search optimization algorithm; By changing the individual inclusion , , , , Parameters, to achieve iterative update of individuals within the population; (5); In the formula, The number of individuals in a population , , , , The step size of the parameter change, where is the upper bound of the constraint condition, is the lower bound of the constraint; For the Generation population; For the Generation population; It means to generate an array of random numbers uniformly distributed between 0 and 1; Indicates obtaining the number of rows and columns of the target matrix; For individuals exceeding the constraints , , , , The parameters are controlled by upper and lower limits according to equations (6) and (7): Exceeding the limit: (6); Exceeded lower limit: (7); In the formula, Represents an individual The mth one contained in , , , , parameter; The mth upper limit , , , , parameter; The mth lower limit , , , , parameter; Substitute all updated individuals in the population into the JA hysteresis model to obtain the actual magnetization intensity at each moment , determine the target magnetic induction intensity B(t) at each moment based on the target magnetic field intensity at each moment and the actual magnetization intensity at each moment; The fitness function is used to calculate the fitness value of each individual in the population after each iterative update. The fitness function F is expressed as: (8); In the formula, is the target magnetic flux density data calculated by the JA hysteresis model. Values; is the jth value in the reference magnetic induction intensity data; is the number of data points; Determine whether the termination condition is reached. When the termination condition is reached, select the individual with the smallest fitness value as the optimal output. Otherwise, determine whether the number of iterative updates is the specified value. When the number of iterative updates is the specified value, reinitialize. Otherwise, start the next iterative update of individuals in the population. 7. A three-dimensional time-domain hysteresis loop drawing method based on a hysteresis model according to claim 6, characterized in that: the target magnetic field intensity and the target magnetic induction intensity are subjected to curve fitting, and the specific process of obtaining the three-dimensional hysteresis loop of the target iron core is as follows: constructing a three-dimensional rectangular coordinate system with the x-axis, the y-axis, and the z-axis as coordinate axes, and fitting the target magnetic field intensity and target magnetic induction intensity The three-dimensional hysteresis loop of the target core is drawn by taking the x-axis and y-axis data of the three-dimensional rectangular coordinate system respectively and the time t as the z-axis data of the three-dimensional rectangular coordinate system. 8. A three-dimensional time-domain hysteresis loop depiction system based on a hysteresis model, characterized by comprising: A construction module is used to construct a J-A hysteresis model corresponding to the target core; A calculation module is used to obtain the current values of the dynamic stray current of the target iron core at multiple different times, calculate the magnetic field strength generated by the dynamic stray current according to the obtained current value, and superimpose the magnetic field strength generated by the dynamic stray current with the AC magnetic field strength generated by the normal power frequency AC current inside the target iron core to obtain the target magnetic field strength as the total input of the J-A hysteresis model; An identification module is used to perform parameter identification on the J-A hysteresis model after inputting the target magnetic field strength by using a backtracking search optimization algorithm, substitute the model parameters of the J-A hysteresis model obtained after the parameter identification into the J-A hysteresis model, and obtain the corresponding actual magnetization intensity based on the target magnetic field strength on this basis, and obtain the target magnetic induction intensity by using the target magnetic field strength and the actual magnetization intensity; The fitting module is used to perform curve fitting on the target magnetic field intensity and the target magnetic induction intensity to obtain the three-dimensional hysteresis loop of the target iron core. 9. A three-dimensional time-domain hysteresis loop drawing system based on a hysteresis model according to claim 8, characterized in that: the J-A hysteresis model corresponding to the target core is expressed as: ; In the formula, Indicates the actual magnetization Target magnetic field strength The differential of represents a differential operator; The coefficient introduced to prevent non-physical solutions; is the reversible magnetic susceptibility; is the irreversible loss coefficient; is the directional coefficient. When the rate of change of the target magnetic field strength dH(t)/dt>0, , when dH(t)/dt<0, ; t represents time; is the domain wall interaction coefficient; is the magnetization intensity without hysteresis effect; The hysteresis-free magnetization Expressed as: ; In the formula, is the saturation magnetization; is the magnetization intensity shape factor without hysteresis effect; represents the effective magnetic field strength, where is the domain wall interaction coefficient. 10. A computer-readable storage medium storing a computer program, characterized in that when the computer program is executed by a processor, it implements the three-dimensional time-domain hysteresis loop depiction method based on the hysteresis model described in any one of claims 1 to 7.