CN118378576A - Three-dimensional time domain hysteresis loop drawing method and system based on hysteresis model - Google Patents
Three-dimensional time domain hysteresis loop drawing method and system based on hysteresis model Download PDFInfo
- Publication number
- CN118378576A CN118378576A CN202410841934.1A CN202410841934A CN118378576A CN 118378576 A CN118378576 A CN 118378576A CN 202410841934 A CN202410841934 A CN 202410841934A CN 118378576 A CN118378576 A CN 118378576A
- Authority
- CN
- China
- Prior art keywords
- magnetic field
- hysteresis
- target
- field intensity
- intensity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 84
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical group [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims abstract description 44
- 230000005415 magnetization Effects 0.000 claims abstract description 41
- 230000006698 induction Effects 0.000 claims abstract description 31
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 22
- 238000005457 optimization Methods 0.000 claims abstract description 21
- 230000008569 process Effects 0.000 claims description 35
- 230000000694 effects Effects 0.000 claims description 12
- 230000003993 interaction Effects 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 11
- 230000002427 irreversible effect Effects 0.000 claims description 6
- 230000002441 reversible effect Effects 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000004804 winding Methods 0.000 claims description 3
- 238000009877 rendering Methods 0.000 claims 2
- 238000010586 diagram Methods 0.000 description 10
- 230000006870 function Effects 0.000 description 9
- 239000000696 magnetic material Substances 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000004075 alteration Effects 0.000 description 2
- 230000007935 neutral effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- Mathematical Physics (AREA)
- Biomedical Technology (AREA)
- Artificial Intelligence (AREA)
- Computer Graphics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
The invention discloses a three-dimensional time domain hysteresis loop drawing method and system based on a hysteresis model, wherein the method comprises the following steps: constructing a J-A hysteresis model corresponding to the target iron core; acquiring dynamic stray current values of a target iron core at different moments, calculating magnetic field intensity according to the current values, and superposing the dynamic stray magnetic field intensity and normal alternating current magnetic field intensity to be used as the total input of a J-A hysteresis model; performing parameter identification on the J-A hysteresis model by adopting a backtracking search optimization algorithm, and determining target magnetic induction intensity based on target magnetic field intensity and actual magnetization intensity; curve fitting is carried out on the final magnetic field intensity and the magnetic induction intensity, and a three-dimensional hysteresis loop of the target iron core is obtained; according to the invention, the magnetic field intensity generated by the dynamic stray current is overlapped with the normal alternating current magnetic field intensity to be used as the total input of a model, and the hysteresis loop under the three-dimensional condition of considering the continuous interference of the dynamic stray current is drawn, so that the time variation and the non-periodicity of the influence of the dynamic stray current on the hysteresis loop are shown.
Description
Technical Field
The invention relates to the technical field of iron core magnetization characteristics, in particular to a three-dimensional time domain hysteresis loop drawing method and system based on a hysteresis model.
Background
With the rapid development of urban rail transit, dynamic stray current generated by incomplete insulation of a running rail to the ground is gradually increased, and the influence on a transformer is increasingly increased. The dynamic stray current, namely the stray current which dynamically changes along with time, can cause a series of problems of aggravation of vibration noise, increase of loss, increase of exciting current harmonic wave and the like of the transformer after flowing into the transformer through a neutral point, so that the magnetic property of the transformer core under the interference of the dynamic stray current is very important to be explored.
The parameter identification of the hysteresis loop is used as one of key means for researching the magnetic characteristics of the iron core, and the existing research is mainly aimed at carrying out parameter identification on the hysteresis loop under the working conditions of normal working conditions, changed material characteristics and the like, so that the drawing of the hysteresis loop is more accurate. The traditional hysteresis loop is used for representing hysteresis of magnetic substances under the periodically-changed magnetic field intensity, and shows that the relationship between the magnetic induction intensity B and the magnetic field intensity H of the magnetic substances in the repeated magnetization process has obvious periodicity. However, the magnetic field intensity generated by the dynamic stray current changes dynamically along with time due to the dynamic change of the current value of the dynamic stray current along with time, so that the influence on the magnetic characteristics of the iron core changes along with time, and the magnetic core has no periodicity any more. Because the traditional hysteresis loop can only represent the change relation between H and B on a two-dimensional layer, and neglect the change relation between H and B on a time domain, the traditional hysteresis loop can not accurately express the change rule of H and B on the time domain, namely the hysteresis characteristic of the magnetic material under the influence of dynamic stray current.
Disclosure of Invention
Aiming at the defects of the prior art, the invention aims to provide a method and a system capable of accurately describing the change condition of hysteresis characteristics of a magnetic material under the influence of dynamic stray current, and provides assistance for the subsequent research of hysteresis loss and the working temperature of the magnetic material.
In order to achieve the above purpose, the present invention provides the following technical solutions: a three-dimensional time domain hysteresis loop drawing method based on a hysteresis model comprises the following steps:
constructing a J-A hysteresis model corresponding to the target iron core;
obtaining current values of dynamic stray currents of a target iron core at different moments, calculating the magnetic field intensity generated by the dynamic stray currents according to the obtained current values, and superposing the magnetic field intensity generated by the dynamic stray currents and the alternating current magnetic field intensity generated by normal power frequency alternating current in the target iron core to obtain the target magnetic field intensity As the total input to the J-A hysteresis model;
the retrospective search optimization algorithm is adopted to input the magnetic field intensity of the target The J-A hysteresis model is subjected to parameter identification, the model parameters of the J-A hysteresis model obtained after the parameter identification are substituted into the J-A hysteresis model, and the model parameters are based on the target magnetic field intensityObtain the corresponding actual magnetization intensityBy means of the strength of the target magnetic fieldAnd the actual magnetizationObtaining target magnetic induction intensity B (t);
And 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 iron core.
Further, the J-A hysteresis model corresponding to the target iron core is expressed as:
（1）；
in the method, in the process of the invention, Representing the actual magnetizationFor the target magnetic field intensityIs a derivative of (2); representing a differential operator; coefficients elicited to prevent the occurrence of non-physical solutions; is reversible magnetic coefficient; Is the irreversible loss coefficient; as the direction coefficient, when the change rate dH (t)/dt >0 of the target magnetic field strength, When dH (t)/dt <0,; T represents time; in order to have a magnetization without hysteresis effect, Is the domain wall interaction coefficient.
Further, the magnetization without hysteresis effectExpressed as:
（2）；
in the method, in the process of the invention, Is saturation magnetization; is a magnetization shape factor without hysteresis effect; Indicating the effective magnetic field strength, wherein, Is the domain wall interaction coefficient.
Further, the specific process of calculating the magnetic field intensity generated by the dynamic stray current according to the acquired current value is as follows: solving magnetic field strength generated by dynamic stray current based on ampere loop law：
（3）；
In the method, in the process of the invention,Turns of the winding coil; Is a dynamic stray current; Is the effective magnetic path length.
Further, the magnetic field strength generated by the dynamic stray current is improvedThe intensity of the alternating current magnetic field generated in the target iron core by the normal power frequency alternating currentTarget magnetic field strength obtained by superpositionThe total input as a J-A hysteresis model is represented by equation (4):
（4）。
further, a backtracking search optimization algorithm is adopted to conduct parameter identification on the J-A hysteresis model, model parameters obtained after the parameter identification are substituted into the J-A hysteresis model, and the model parameters are based on the target magnetic field intensity Obtain the corresponding actual magnetization intensityBy means of the strength of the target magnetic fieldAnd the actual magnetizationThe specific process for obtaining the target magnetic induction B (t) is as follows:
initializing parameters of a backtracking search optimization algorithm;
Setting the population size and population iteration times of a backtracking search optimization algorithm, simultaneously giving a step factor gamma and a random disturbance factor r, and initializing in a constraint condition to generate individuals meeting the constraint condition, wherein each individual is composed of 、、、、The parameters in the total five J-A hysteresis models are composed,、、、、The value range of (2) constitutes a constraint condition;
carrying out iterative updating on individuals in the population of the backtracking search optimization algorithm;
altering the inclusion of an individual by formula (5) 、、、、Parameters, realizing iterative updating of individuals in the population;
（5）；
in the method, in the process of the invention, For inclusion in individuals within a population、、、、A step of changing the parameter, wherein,As an upper bound to the constraint condition,Is the lower bound of the constraint condition; Is the first Generating a population; Is the first Generating a population; Representing the generation of an array consisting of random numbers uniformly distributed between 0 and 1; representing the number of rows and columns of the acquisition target matrix;
For individuals who exceed the constraints 、、、、The parameters are controlled to be upper and lower limits according to the formula (6) and the formula (7):
Exceeding the upper limit: （6）；
exceeding the lower limit: （7）；
in the method, in the process of the invention, Representing individualsComprises the m < th > of、、、、Parameters; Mth of the upper limit 、、、、Parameters; Mth of the lower limit 、、、、Parameters;
Substituting the updated individuals in all the populations into the J-A hysteresis model to obtain the actual magnetization intensity at each moment Determining a target magnetic induction B (t) at each moment based on the target magnetic field strength at each moment and the actual magnetization at each moment;
calculating the fitness value of each individual in the population after each iteration update by using a fitness function, wherein the fitness function F is expressed as:
（8）；
in the method, in the process of the invention, To the first of the target magnetic induction data calculated by J-A hysteresis modelA value; A j-th value in the reference magnetic induction data; The number of data points;
judging whether a termination condition is reached, when the termination condition is reached, selecting an individual with the smallest fitness value as an optimal output, otherwise judging whether the iteration update times are specified values, and reinitializing when the iteration update times are specified values, otherwise starting the iteration update of the next individual in the population.
Further, curve fitting is carried out on the target magnetic field intensity and the target magnetic induction intensity, and the specific process for obtaining the three-dimensional hysteresis loop of the target iron core is as follows: constructing a three-dimensional rectangular coordinate system taking an x-axis, a y-axis and a z-axis as coordinate axes, and enabling the target magnetic field intensity to be equal to that of the target magnetic field intensityAnd target magnetic induction intensityRespectively using the time t as the data of the x axis and the y axis of the three-dimensional rectangular coordinate system, and drawing the three-dimensional hysteresis loop of the target iron core by using the time t as the data of the z axis of the three-dimensional rectangular coordinate system.
A three-dimensional time domain hysteresis loop tracing system based on a hysteresis model, comprising:
The construction module is used for constructing a J-A hysteresis model corresponding to the target iron core;
The calculation module is used for obtaining current values of dynamic stray currents of the target iron core at different moments, calculating the magnetic field intensity generated by the dynamic stray currents according to the obtained current values, and taking the target magnetic field intensity obtained by superposing the magnetic field intensity generated by the dynamic stray currents and the alternating current magnetic field intensity generated by normal power frequency alternating current in the target iron core as the total input of the J-A hysteresis model;
The identification module is used for carrying out parameter identification on the J-A hysteresis model after inputting the target magnetic field intensity by adopting a retrospective search optimization algorithm, substituting model parameters of the J-A hysteresis model obtained after parameter identification into the J-A hysteresis model, obtaining corresponding actual magnetization intensity based on the target magnetic field intensity on the basis, and obtaining target magnetic induction intensity by utilizing the target magnetic field intensity and the actual magnetization intensity;
and the fitting module is used for 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 iron core.
Further, the J-A hysteresis model corresponding to the target iron core is expressed as:
；
in the method, in the process of the invention, Representing the actual magnetizationFor the target magnetic field intensityIs a derivative of (2); representing a differential operator; coefficients elicited to prevent the occurrence of non-physical solutions; is reversible magnetic coefficient; Is the irreversible loss coefficient; as the direction coefficient, when the change rate dH (t)/dt >0 of the target magnetic field strength, When dH (t)/dt <0,; T represents time; in order to have a magnetization without hysteresis effect, Is the domain wall interaction coefficient;
said hysteresis-free magnetization Expressed as:
；
in the method, in the process of the invention, Is saturation magnetization; is a magnetization shape factor without hysteresis effect; Indicating the effective magnetic field strength, wherein, Is the domain wall interaction coefficient.
A computer readable storage medium storing a computer program which when executed by a processor implements a method of three-dimensional time domain hysteresis loop tracing based on a hysteresis model.
Compared with the prior art, the invention has the following beneficial effects: according to the invention, the magnetic field intensity generated by the dynamic stray current is overlapped with the normal alternating current magnetic field intensity to be used as the total input of a model, the hysteresis loop under the three-dimensional condition of considering the continuous interference of the dynamic stray current is drawn, the time variation and the non-periodicity of the influence of the dynamic stray current on the hysteresis loop are displayed, meanwhile, the variation condition of the hysteresis loop in each normal alternating current magnetic field intensity variation period can be observed through the hysteresis loop under the three-dimensional condition, and the hysteresis loss and the temperature variation of the magnetic material during working are further evaluated.
Drawings
FIG. 1 is a flow chart of the present invention for parameter identification of J-A hysteresis models under dynamic spurious current interference.
FIG. 2 is a flow chart of the intelligent algorithm parameter identification according to the present invention.
Fig. 3 is a waveform diagram of the neutral point current of the actual measurement transformer substation according to the present invention.
Fig. 4 is a waveform diagram of the superimposed dynamic stray current field strength and ac field strength of the present invention.
FIG. 5 is a three-dimensional time domain comparison chart of the J-A hysteresis model parameter identification of the present invention.
Fig. 6 is a schematic diagram of a system structure according to the present invention.
Detailed Description
As shown in fig. 1, the present invention provides the following technical solutions: a three-dimensional hysteresis loop drawing method based on a hysteresis model comprises the following steps:
Step S1: constructing a J-A hysteresis model corresponding to the target iron core:
（1）；
in the method, in the process of the invention, Representing the actual magnetizationFor the target magnetic field intensityIs a derivative of (2); representing a differential operator; coefficients elicited to prevent the occurrence of non-physical solutions; is reversible magnetic coefficient; Is the irreversible loss coefficient; as the direction coefficient, when the change rate dH (t)/dt >0 of the target magnetic field strength, When dH (t)/dt <0,; T represents time; Is the domain wall interaction coefficient; for magnetization without hysteresis, expressed as:
（2）；
in the method, in the process of the invention, Is saturation magnetization; is a magnetization shape factor without hysteresis effect; Indicating the effective magnetic field strength, wherein, Is the domain wall interaction coefficient.
Step S2: obtaining current values of dynamic stray currents of a target iron core at a plurality of different moments, as shown in fig. 3, calculating the magnetic field intensity generated by the dynamic stray currents according to the obtained current values, and superposing the magnetic field intensity generated by the dynamic stray currents and the alternating current magnetic field intensity generated by normal power frequency alternating current in the target iron core to obtain the target magnetic field intensityAs the overall input to the J-a hysteresis model.
Step S3: the retrospective search optimization algorithm is adopted to input the magnetic field intensity of the targetThe J-A hysteresis model is subjected to parameter identification, the model parameters of the J-A hysteresis model obtained after the parameter identification are substituted into the J-A hysteresis model, and the model parameters are based on the target magnetic field intensityObtain the corresponding actual magnetization intensityBy means of the strength of the target magnetic fieldAnd the actual magnetizationThe target magnetic induction B (t) is obtained as shown in fig. 2.
Step S4: and 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 iron core.
The specific process for calculating the magnetic field intensity generated by the dynamic stray current according to the acquired current value comprises the following steps: solving magnetic field strength generated by dynamic stray current based on ampere loop law：
（3）；
In the method, in the process of the invention,Turns of the winding coil; Is a dynamic stray current; Is the effective magnetic path length.
Wherein the magnetic field strength generated by the dynamic stray currentIntensity of alternating current magnetic field generated by normal power frequency alternating currentTarget magnetic field strength obtained by superpositionThe total input as the J-A hysteresis model is shown in the formula (4), and the waveform diagram of the calculated target magnetic field intensity is shown in the figure 4;
（4）；
the specific process of step S3 is as follows:
step S3.1: initializing parameters of a backtracking search optimization algorithm;
Setting the population size and population iteration times of a backtracking search optimization algorithm, simultaneously giving a step factor gamma and a random disturbance factor r, and initializing in a constraint condition to generate individuals meeting the constraint condition, wherein each individual is composed of 、The parameters in five J-A hysteresis models are composed of k, a and c,、The value ranges of k, a and c constitute constraint conditions.
Step S3.2: carrying out iterative updating on individuals in the population of the backtracking search optimization algorithm;
altering the inclusion of an individual by formula (5) 、、、、The parameters realize iterative updating of individuals in the population;
（5）；
in the method, in the process of the invention, For inclusion in individuals within a population、、、、A step of changing the parameter, wherein,As an upper bound to the constraint condition,Is the lower bound of the constraint condition; Is the first Generating a population; Is the first Generating a population; Representing the generation of an array consisting of random numbers uniformly distributed between 0 and 1; representing the number of rows and columns of the acquisition target matrix.
Step S3.3: for individuals who exceed the constraints、、、、The parameters are controlled to be upper and lower limits according to the formula (6) and the formula (7):
Exceeding the upper limit: （6）；
exceeding the lower limit: （7）；
in the method, in the process of the invention, Representing individualsComprises the m < th > of、、、、Parameters; Mth of the upper limit 、、、、Parameters; Mth of the lower limit 、、、、Parameters.
Step S3.4: substituting the updated individuals in all the populations into a J-A hysteresis model to obtain the actual magnetization M at each moment, and determining the target magnetic induction B (t) at each moment based on the target magnetic field strength at each moment and the actual magnetization at each moment;
calculating the fitness value of each individual in the population after each iteration update by using a fitness function, wherein the smaller the fitness value is, the more excellent the individual is, and the fitness function F is expressed as:
（8）；
in the method, in the process of the invention, To the first of the target magnetic induction data calculated by J-A hysteresis modelA value; A j-th value in the reference magnetic induction data; the number of data points; the smaller the value of the fitness function F, the closer the J-A hysteresis model parameters identified by the retrospective search optimization algorithm to the true values.
Step S3.5: judging whether a termination condition is reached, when the termination condition is reached, selecting an individual with the smallest fitness value as an optimal output, otherwise judging whether the iteration update times are integral multiples of 5, and reinitializing when the iteration update times are integral multiples of 5, otherwise starting the iteration update of the individual in the next population.
The specific process of step S4 is: constructing a three-dimensional rectangular coordinate system taking an x-axis, a y-axis and a z-axis as coordinate axes, and enabling the target magnetic field intensity to be equal to that of the target magnetic field intensityAnd target magnetic induction intensityRespectively using the time t as the data of the x axis and the y axis of the three-dimensional rectangular coordinate system, and drawing the three-dimensional hysteresis loop of the target iron core by using the time t as the data of the z axis of the three-dimensional rectangular coordinate system.
The three-dimensional hysteresis loop of the target iron core obtained by performing parameter identification on the J-A hysteresis model by using the backtracking search optimization algorithm is shown in fig. 5, and as can be seen from fig. 5, the three-dimensional hysteresis loop of the target iron core obtained by performing parameter identification on the J-A hysteresis model by using the backtracking search optimization algorithm is similar to the reference value, and a better effect can be obtained by adopting the scheme.
As shown in fig. 6, a three-dimensional time domain hysteresis loop plotting system based on hysteresis model includes:
The construction module is used for constructing a J-A hysteresis model corresponding to the target iron core;
The calculation module is used for obtaining current values of dynamic stray currents of the target iron core at different moments, calculating the magnetic field intensity generated by the dynamic stray currents according to the obtained current values, and taking the target magnetic field intensity obtained by superposing the magnetic field intensity generated by the dynamic stray currents and the alternating current magnetic field intensity generated by normal power frequency alternating current in the target iron core as the total input of the J-A hysteresis model;
The identification module is used for carrying out parameter identification on the J-A hysteresis model after inputting the target magnetic field intensity by adopting a retrospective search optimization algorithm, substituting model parameters of the J-A hysteresis model obtained after parameter identification into the J-A hysteresis model, obtaining corresponding actual magnetization intensity based on the target magnetic field intensity on the basis, and obtaining target magnetic induction intensity by utilizing the target magnetic field intensity and the actual magnetization intensity;
and the fitting module is used for 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 iron core.
The J-A hysteresis model corresponding to the target iron core is expressed as follows:
；
in the method, in the process of the invention, Representing the actual magnetizationFor the target magnetic field intensityIs a derivative of (2); representing a differential operator; coefficients elicited to prevent the occurrence of non-physical solutions; is reversible magnetic coefficient; Is the irreversible loss coefficient; as the direction coefficient, when the change rate dH (t)/dt >0 of the target magnetic field strength, When dH (t)/dt <0,; T represents time; in order to have a magnetization without hysteresis effect, Is the domain wall interaction coefficient.
Wherein the hysteresis-free magnetizationExpressed as:
；
in the method, in the process of the invention, Is saturation magnetization; is a magnetization shape factor without hysteresis effect; Indicating the effective magnetic field strength, wherein, Is the domain wall interaction coefficient.
A computer readable storage medium storing a computer program which when executed by a processor implements a method of three-dimensional time domain hysteresis loop tracing based on a hysteresis model.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein. The scheme in the embodiment of the application can be realized by adopting various computer languages, such as object-oriented programming language Java, an transliteration script language JavaScript and the like.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
While preferred embodiments of the present application have been described, additional variations and modifications in those embodiments may occur to those skilled in the art once they learn of the basic inventive concepts. It is therefore intended that the following claims be interpreted as including the preferred embodiments and all such alterations and modifications as fall within the scope of the application.
It will be apparent to those skilled in the art that various modifications and variations can be made to the present application without departing from the spirit or scope of the application. Thus, it is intended that the present application also include such modifications and alterations insofar as they come within the scope of the appended claims or the equivalents thereof.
Claims (10)
1. A three-dimensional time domain hysteresis loop drawing method based on a hysteresis model is characterized by comprising the following steps:
constructing a J-A hysteresis model corresponding to the target iron core;
obtaining current values of dynamic stray currents of a target iron core at different moments, calculating the magnetic field intensity generated by the dynamic stray currents according to the obtained current values, and superposing the magnetic field intensity generated by the dynamic stray currents and the alternating current magnetic field intensity generated by normal power frequency alternating current in the target iron core to obtain the target magnetic field intensity As the total input to the J-A hysteresis model;
the retrospective search optimization algorithm is adopted to input the magnetic field intensity of the target The J-A hysteresis model is subjected to parameter identification, the model parameters of the J-A hysteresis model obtained after the parameter identification are substituted into the J-A hysteresis model, and the model parameters are based on the target magnetic field intensityObtain the corresponding actual magnetization intensityBy means of the strength of the target magnetic fieldAnd the actual magnetizationObtaining target magnetic induction intensity B (t);
And 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 iron core.
2. The three-dimensional time domain hysteresis loop plotting method based on hysteresis model as claimed in claim 1, wherein the method comprises the following steps: the J-A hysteresis model corresponding to the target iron core is expressed as:
（1）；
in the method, in the process of the invention, Representing the actual magnetizationFor the target magnetic field intensityIs a derivative of (2); representing a differential operator; coefficients elicited to prevent the occurrence of non-physical solutions; is reversible magnetic coefficient; Is the irreversible loss coefficient; as the direction coefficient, when the change rate dH (t)/dt >0 of the target magnetic field strength, When dH (t)/dt <0,; T represents time; Is magnetization intensity without hysteresis effect; Is the domain wall interaction coefficient.
3. The three-dimensional time domain hysteresis loop plotting method based on hysteresis model as claimed in claim 2, wherein the method comprises the following steps: said hysteresis-free magnetizationExpressed as:
（2）；
in the method, in the process of the invention, Is saturation magnetization; is a magnetization shape factor without hysteresis effect; Indicating the effective magnetic field strength, wherein, Is the domain wall interaction coefficient.
4. A method for three-dimensional time domain hysteresis loop plotting based on hysteresis model as defined in claim 3, wherein: the specific process for calculating the magnetic field intensity generated by the dynamic stray current according to the acquired current value is as follows: solving magnetic field strength generated by dynamic stray current based on ampere loop law：
（3）；
In the method, in the process of the invention,Turns of the winding coil; Is a dynamic stray current; Is the effective magnetic path length.
5. The method for three-dimensional time domain hysteresis loop plotting based on hysteresis model as defined in claim 4, wherein the method comprises the following steps: magnetic field strength generated by dynamic stray currentThe intensity of the alternating current magnetic field generated in the target iron core by the normal power frequency alternating currentTarget magnetic field strength obtained by superpositionThe total input as a J-A hysteresis model is represented by equation (4):
（4）。
6. The method for three-dimensional time domain hysteresis loop plotting based on hysteresis model according to claim 5, wherein the method comprises the following steps: performing parameter identification on the J-A hysteresis model by adopting a backtracking search optimization algorithm, substituting model parameters obtained after the parameter identification into the J-A hysteresis model, and based on the target magnetic field intensity Obtain the corresponding actual magnetization intensityBy means of the strength of the target magnetic fieldAnd the actual magnetizationThe specific process for obtaining the target magnetic induction B (t) is as follows:
initializing parameters of a backtracking search optimization algorithm;
Setting the population size and population iteration times of a backtracking search optimization algorithm, simultaneously giving a step factor gamma and a random disturbance factor r, and initializing in a constraint condition to generate individuals meeting the constraint condition, wherein each individual is composed of 、、、、The parameters in the total five J-A hysteresis models are composed,、、、、The value range of (2) constitutes a constraint condition;
carrying out iterative updating on individuals in the population of the backtracking search optimization algorithm;
altering the inclusion of an individual by formula (5) 、、、、Parameters, realizing iterative updating of individuals in the population;
（5）；
in the method, in the process of the invention, For inclusion in individuals within a population、、、、A step of changing the parameter, wherein,As an upper bound to the constraint condition,Is the lower bound of the constraint condition; Is the first Generating a population; Is the first Generating a population; Representing the generation of an array consisting of random numbers uniformly distributed between 0 and 1; representing the number of rows and columns of the acquisition target matrix;
For individuals who exceed the constraints 、、、、The parameters are controlled to be upper and lower limits according to the formula (6) and the formula (7):
Exceeding the upper limit: （6）；
exceeding the lower limit: （7）；
in the method, in the process of the invention, Representing individualsComprises the m < th > of、、、、Parameters; Mth of the upper limit 、、、、Parameters; Mth of the lower limit 、、、、Parameters;
Substituting the updated individuals in all the populations into the J-A hysteresis model to obtain the actual magnetization intensity at each moment Determining a target magnetic induction B (t) at each moment based on the target magnetic field strength at each moment and the actual magnetization at each moment;
calculating the fitness value of each individual in the population after each iteration update by using a fitness function, wherein the fitness function F is expressed as:
（8）；
in the method, in the process of the invention, To the first of the target magnetic induction data calculated by J-A hysteresis modelA value; A j-th value in the reference magnetic induction data; The number of data points;
judging whether a termination condition is reached, when the termination condition is reached, selecting an individual with the smallest fitness value as an optimal output, otherwise judging whether the iteration update times are specified values, and reinitializing when the iteration update times are specified values, otherwise starting the iteration update of the next individual in the population.
7. The method for three-dimensional time domain hysteresis loop plotting based on hysteresis model as defined in claim 6, wherein the method comprises the following steps: curve fitting is carried out on the target magnetic field intensity and the target magnetic induction intensity, and the specific process for obtaining the three-dimensional hysteresis loop of the target iron core is as follows: constructing a three-dimensional rectangular coordinate system taking an x-axis, a y-axis and a z-axis as coordinate axes, and enabling the target magnetic field intensity to be equal to that of the target magnetic field intensityAnd target magnetic induction intensityRespectively using the time t as the data of the x axis and the y axis of the three-dimensional rectangular coordinate system, and drawing the three-dimensional hysteresis loop of the target iron core by using the time t as the data of the z axis of the three-dimensional rectangular coordinate system.
8. A three-dimensional time domain hysteresis loop plotting system based on hysteresis model, characterized by comprising:
The construction module is used for constructing a J-A hysteresis model corresponding to the target iron core;
The calculation module is used for obtaining current values of dynamic stray currents of the target iron core at different moments, calculating the magnetic field intensity generated by the dynamic stray currents according to the obtained current values, and taking the target magnetic field intensity obtained by superposing the magnetic field intensity generated by the dynamic stray currents and the alternating current magnetic field intensity generated by normal power frequency alternating current in the target iron core as the total input of the J-A hysteresis model;
The identification module is used for carrying out parameter identification on the J-A hysteresis model after inputting the target magnetic field intensity by adopting a retrospective search optimization algorithm, substituting model parameters of the J-A hysteresis model obtained after parameter identification into the J-A hysteresis model, obtaining corresponding actual magnetization intensity based on the target magnetic field intensity on the basis, and obtaining target magnetic induction intensity by utilizing the target magnetic field intensity and the actual magnetization intensity;
and the fitting module is used for 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 iron core.
9. The hysteresis model-based three-dimensional time domain hysteresis loop rendering system of claim 8, wherein: the J-A hysteresis model corresponding to the target iron core is expressed as:
；
in the method, in the process of the invention, Representing the actual magnetizationFor the target magnetic field intensityIs a derivative of (2); representing a differential operator; coefficients elicited to prevent the occurrence of non-physical solutions; is reversible magnetic coefficient; Is the irreversible loss coefficient; as the direction coefficient, when the change rate dH (t)/dt >0 of the target magnetic field strength, When dH (t)/dt <0,; T represents time; Is the domain wall interaction coefficient; Is magnetization intensity without hysteresis effect;
said hysteresis-free magnetization Expressed as:
；
in the method, in the process of the invention, Is saturation magnetization; is a magnetization shape factor without hysteresis effect; Indicating the effective magnetic field strength, wherein, Is the domain wall interaction coefficient.
10. A computer readable storage medium storing a computer program, wherein the computer program when executed by a processor implements the three-dimensional time domain hysteresis loop rendering method based on a hysteresis model as claimed in any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410841934.1A CN118378576B (en) | 2024-06-27 | 2024-06-27 | Three-dimensional time domain hysteresis loop drawing method and system based on hysteresis model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410841934.1A CN118378576B (en) | 2024-06-27 | 2024-06-27 | Three-dimensional time domain hysteresis loop drawing method and system based on hysteresis model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN118378576A true CN118378576A (en) | 2024-07-23 |
| CN118378576B CN118378576B (en) | 2024-10-18 |
Family
ID=91906163
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410841934.1A Active CN118378576B (en) | 2024-06-27 | 2024-06-27 | Three-dimensional time domain hysteresis loop drawing method and system based on hysteresis model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118378576B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107609298A (en) * | 2017-09-26 | 2018-01-19 | 广东电网有限责任公司电力科学研究院 | A kind of Jiles Atherton model parameter identification methods and device |
| CN109444776A (en) * | 2018-11-05 | 2019-03-08 | 国网江苏省电力有限公司电力科学研究院 | Three-phase full-bridge converter iron core remanent magnetism measuring method, system and storage medium |
| CN112733429A (en) * | 2020-12-20 | 2021-04-30 | 南京理工大学 | Method for describing neodymium iron boron hysteresis loop based on Jiles-Atherton model |
| CN114417555A (en) * | 2021-12-14 | 2022-04-29 | 深圳供电局有限公司 | Transformer modeling method and device, computer equipment and storage medium |
| CN115587448A (en) * | 2022-10-31 | 2023-01-10 | 国网四川省电力公司电力科学研究院 | Method, device, equipment and medium for identifying performance of current transformer iron core material |
| CN115758861A (en) * | 2022-10-10 | 2023-03-07 | 三峡大学 | Method for calculating dynamic hysteresis loop of dynamic analytic inverse Preisach model based on R-L type fractional derivative |
| RU2805248C1 (en) * | 2022-12-05 | 2023-10-12 | Роман Александрович Соколов | Device for measuring the magnetic characteristics of a ferromagnet |
| CN117421861A (en) * | 2023-09-15 | 2024-01-19 | 南方电网数字电网科技(广东)有限公司 | Power transmission line multi-working-condition simulation analysis method based on real-time multi-source data fusion |
-
2024
- 2024-06-27 CN CN202410841934.1A patent/CN118378576B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107609298A (en) * | 2017-09-26 | 2018-01-19 | 广东电网有限责任公司电力科学研究院 | A kind of Jiles Atherton model parameter identification methods and device |
| CN109444776A (en) * | 2018-11-05 | 2019-03-08 | 国网江苏省电力有限公司电力科学研究院 | Three-phase full-bridge converter iron core remanent magnetism measuring method, system and storage medium |
| CN112733429A (en) * | 2020-12-20 | 2021-04-30 | 南京理工大学 | Method for describing neodymium iron boron hysteresis loop based on Jiles-Atherton model |
| CN114417555A (en) * | 2021-12-14 | 2022-04-29 | 深圳供电局有限公司 | Transformer modeling method and device, computer equipment and storage medium |
| CN115758861A (en) * | 2022-10-10 | 2023-03-07 | 三峡大学 | Method for calculating dynamic hysteresis loop of dynamic analytic inverse Preisach model based on R-L type fractional derivative |
| CN115587448A (en) * | 2022-10-31 | 2023-01-10 | 国网四川省电力公司电力科学研究院 | Method, device, equipment and medium for identifying performance of current transformer iron core material |
| RU2805248C1 (en) * | 2022-12-05 | 2023-10-12 | Роман Александрович Соколов | Device for measuring the magnetic characteristics of a ferromagnet |
| CN117421861A (en) * | 2023-09-15 | 2024-01-19 | 南方电网数字电网科技(广东)有限公司 | Power transmission line multi-working-condition simulation analysis method based on real-time multi-source data fusion |
Non-Patent Citations (4)
| Title |
|---|
| HESSA AL JUNAID: "Timeless Discretization of Magnetization Slop in the Modeling of Ferromagnetic Hysteresis", 《IEEE》, 31 December 2007 (2007-12-31) * |
| 周求宽;龙国华;刘衍;周友武;王羽;陈斌;: "电压互感器铁芯J-A磁滞参数特性研究", 武汉大学学报(工学版), no. 04, 28 April 2018 (2018-04-28) * |
| 曹祎: "基于改进粒子群算法的电流互感器 J-A 模型参数辨识", 《电测与仪表》, 15 May 2021 (2021-05-15), pages 70 - 76 * |
| 澹台乐琰: "地铁供电系统引起变压器直流偏磁现象的研究与应用", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑 (月刊)》, 15 July 2021 (2021-07-15) * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN118378576B (en) | 2024-10-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106777442A (en) | A kind of permanent-magnet brushless DC electric machine cogging torque Optimization Design | |
| Liang et al. | Hybrid convolutional neural network and projected entangled pair states wave functions for quantum many-particle states | |
| CN117350152A (en) | Electromagnetic field modeling method, system, equipment and medium based on mechanism data fusion | |
| Yong-Yan et al. | Analytical Treatment of the Generalized Hirota‐Satsuma‐Ito Equation Arising in Shallow Water Wave | |
| CN103020433B (en) | A kind of evaluation method of status of electric power | |
| CN114997027B (en) | Method for intelligently solving random signals of axle system | |
| CN102968519A (en) | Model of electronic current transformer with hollow coil in electrical system and modeling method thereof | |
| CN109472114A (en) | A kind of optimum design method and device of magnetic nano-particle simulation test platform | |
| CN114139377A (en) | Power system modeling method and system, readable storage medium and server | |
| CN109344969B (en) | Neural network system, training method thereof, and computer-readable medium | |
| CN118378576B (en) | Three-dimensional time domain hysteresis loop drawing method and system based on hysteresis model | |
| Seebacher et al. | A pseudo density topology optimization approach in nonlinear electromagnetism applied to a 3D actuator | |
| CN113484919A (en) | Magnetic anomaly inversion method, system, terminal and medium based on convolutional neural network | |
| Mousavi et al. | Bi-objective scheduling for the re-entrant hybrid flow shop with learning effect and setup times | |
| Sosnickaya et al. | Modelling the Electromagnetic Processes and Phenomena in Quantum-Sized Systems in the Course of Physical and Mathematical Support of Master’s Programs for the “Electric Power Engineering, Electrical Engineering and Electromechanics Specialty” | |
| Petzold | Modelling, analysis and simulation of multifrequency induction hardening | |
| CN109408927B (en) | Two-dimensional static magnetic field parallel finite element acceleration method based on black box transmission line model | |
| CN113535828B (en) | Aggregation query method, system, equipment and storage medium of time sequence data | |
| CN116822287A (en) | Method and system for reducing turn-to-turn short circuit multi-coupling field model of converter transformer | |
| CN116151092A (en) | Loss metering method and system for extra-high voltage direct current transmission system | |
| CN107944552A (en) | A kind of industrial Internet of Things parameter prediction method based on Elman neutral nets | |
| JP2941516B2 (en) | DC superposition iron loss calculator | |
| CN117592382B (en) | Dynamic response prediction method, system and medium for railway track bridge system | |
| JP2984108B2 (en) | DC superposition inductance calculator | |
| Marcin et al. | Ferromagnetic core coil hysteresis modeling using fractional derivatives |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |