CN108519569B - JA hysteresis model parameter identification method containing stress - Google Patents
JA hysteresis model parameter identification method containing stress Download PDFInfo
- Publication number
- CN108519569B CN108519569B CN201810426937.3A CN201810426937A CN108519569B CN 108519569 B CN108519569 B CN 108519569B CN 201810426937 A CN201810426937 A CN 201810426937A CN 108519569 B CN108519569 B CN 108519569B
- Authority
- CN
- China
- Prior art keywords
- stress
- magnetic
- experiment
- value
- hysteresis model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/12—Measuring magnetic properties of articles or specimens of solids or fluids
- G01R33/123—Measuring loss due to hysteresis
Abstract
The present invention provides a kind of JA hysteresis model parameter identification method containing stress, comprising steps of S1: magnetism testing device of the setting containing elastic stress；S2: sample to be tested is installed；S3: alternating current is led to primary coil；S4: test sample loading stress is treated；S5: induced current and induced voltage in acquisition secondary coil；S6: calculating and obtain experiment magnetic field strength and experiment magnetic induction intensity, calculates and obtains experiment magnetic hystersis loss value；S7: revised hysteresis model is solved；S8: the peak magnetic induction prediction of strength value, peak field prediction of strength value and the magnetic hystersis loss predicted value that acquire are matched with experiment magnetic field strength, experiment magnetic induction intensity and experiment magnetic hystersis loss value；S9: parameter when using matching error minimum is as optimal value.A kind of JA hysteresis model parameter identification method containing stress of the invention can obtain the parameter distribution with external magnetic field strength and intensity of residual stresses variation, to determine that magnetic hystersis loss provides help.
Description
Technical field
The present invention relates to JA hysteresis model parameter identification field more particularly to a kind of JA hysteresis model parameter knowledges containing stress
Other method.
Background technique
In recent years, non-oriented electrical steel is widely used on the elements such as motor stator, and magnetic property is by more and more
Concern.And influence of the processing and manufacturing process of motor stator for its iron loss (summation of magnetic hystersis loss and eddy-current loss) is very big.
Studies have shown that the iron loss that manufacturing process will lead to electrical sheet material increases above 100%, wherein magnetic hystersis loss accounts for about total losses
60%.And the residual stress that manufacturing process introduces is to influence one of the most important factor of magnetic hystersis loss.
Magnetic hystersis loss be due to ferromagnetic material during magnetized the irreversible motion of neticdomain wall and magnetic domain magnetic away from turn over
Quote.Residual stress will affect the orientation of material internal easy magnetizing axis.According to magnetostrictive relative theory, residual stress
In the presence of can also cause minor change of the material in length and volume, and then affect the power of demagnetizing field.For a kind of mangneto
For the material that coefficient of dilatation is positive, if there is pressure residual stress (σ < 0) along outer magnetic field direction, promote easy magnetizing axis
Influence far from outer magnetic field direction occupies leading role, so that magnetic history becomes difficult, magnetic hystersis loss increases；If along outer
Magnetic direction, which exists, draws residual stress (σ > 0), when residual stress numerical value is smaller, easy magnetizing axis is promoted to tend to outer magnetic field direction
Influence occupy leading role so that magnetic history becomes easy, magnetic hystersis loss reduces, with drawing residual stress to continue to increase
Greatly, it will lead to demagnetizing field enhancing, can magnetic history be become difficult, magnetic hystersis loss increases.
The stress of influence for to(for) magnetic hystersis loss, many scholars have carried out relevant research, wherein 1986 by
It is proposed in the article " Theory of ferromagnetic hysteresis " that D.C.Jiles and D.L.Atherton are delivered
Hysteresis model based on physics can accurately predict the M-H and BH curve of magnetic history, and parameter is less, thus obtain
It is widely applied.Scholar Sablik is based on Jiles-Atherton hysteresis model (hereinafter referred to as JA hysteresis model), has studied
Influence of the stress for magnetic hystersis loss proposes the addition stress influence item in Equivalent Magnetic Field, and gives ferromagnetic material stress
Influence the expression formula of item:
Wherein, M is the intensity of magnetization；σ is residual stress, and positive value represents tensile stress, and negative value represents compression；μ0It is Vacuum Magnetic
Conductance；MsFor saturation magnetization；λ1And λ3It is coefficient.
Parameter in Sablik-Jiles-Atherton hysteresis model shares Ms, C, a, α, k, λ1And λ3.Wherein, MsFor material
The saturation magnetization of material is the property of material itself；C is reversible magnetization coefficient, is definite value；λ1And λ3It also is constant coefficient；And join
Number a, tri- parameters of k, α change with the size of external magnetic field strength and stress.A is anhysteretic susceptibility, determines magnetization
The magnetic induction intensity that later period can reach also is influenced by the pinning effect for magnetizing early period；K be neticdomain wall coefficient of concentration, with rectify
Stupid power is directly proportional；α is main field component, characterizes the slope of hysteresis loop interlude.This three parameters are all σ and Hm(outer magnetic
Peak value) complicated function.There has been no the researchs of this respect at present.Moreover, ground about stress on what magnetic hystersis loss influenced in the past
In studying carefully, experiment compression is up to 50Mpa or so, presses residual stress peak value about 150- much smaller than manufacturing process bring
200Mpa.Relevant research is also lacked for the parameter distribution under biggish pressure residual stress.A set of accurate reasonable parameter is known
Other method is for Sablik-Jiles-Atherton hysteresis model using most important.
In conclusion existing research lacks the hysteresis model parameter identification method containing stress.It is answered therefore it provides one kind contains
The hysteresis model parameter identification method of power is extremely urgent.
Summary of the invention
In view of the deficiency of the prior art, the present invention provides a kind of JA hysteresis model parameter identification side containing stress
Method, in influence process of the residual stress that quantitative analysis manufacturing process is introduced into magnetic hystersis loss, provide one kind can determine compared with
The method of big residual stress lower die shape parameter can obtain the parameter distribution with external magnetic field strength and intensity of residual stresses variation,
To determine that magnetic hystersis loss provides help.
To achieve the goals above, the present invention provides a kind of JA hysteresis model parameter identification method containing stress, including step
It is rapid:
S1: the one magnetism testing device containing elastic stress of setting；The magnetism testing device packet containing elastic stress
It includes: one first sample fixing end, one second sample fixing end, a Stress Control device, a upper magnet yoke, a lower yoke, an insulation
Wire frame and a stress display device；The first sample fixing end and the second sample fixing end are laid along first direction；And
The first sample fixing end can reciprocally be sequentially connected with the Stress Control device along the first direction；On described
Magnetic yoke and the lower yoke are opposed to be set between the first sample fixing end and the second sample fixing end up and down；Institute
It states stress display device and connects the Stress Control device；Two opposed faces of the insulation wire frame are wound with a primary coil respectively
With a secondary coil；
S2: a sample to be tested is installed on the magnetism testing device containing elastic stress；
S3: alternating current is led to the primary coil；
S4: compression or tensile stress are loaded to the sample to be tested by the Stress Control device；
S5: the induced current and induced voltage in the secondary coil are acquired；
S6: calculating and obtain experiment magnetic field strength and experiment magnetic induction intensity, and strong with the reality according to the experiment magnetic field
It tests magnetic induction intensity and calculates acquisition experiment magnetic hystersis loss value；
S7: one revised hysteresis model is solved using MATLAB software, parameter is continued to optimize, acquires a peak value
Magnetic induction intensity predicted value, a peak field prediction of strength value and a magnetic hystersis loss predicted value；
S8: the peak magnetic induction prediction of strength value, the peak field prediction of strength value and the magnetic hystersis loss is pre-
Measured value is matched with the experiment magnetic field strength, the experiment magnetic induction intensity and the experiment magnetic hystersis loss value；
S9: parameter when using matching error minimum is as optimal value.
Preferably, in the S6 step, according to experiment magnetic field strength formula (1) and an experiment magnetic induction intensity formula
(2) it calculates and obtains the experiment magnetic field strength and the experiment magnetic induction intensity:
Wherein, H is experiment magnetic field strength, and B is experiment magnetic induction intensity；N1It is the number of turns of the primary coil, L is described
The length of equivalent magnetic circuit, N in sample to be tested2It is the number of turns of secondary coil, S is the cross-sectional area of the sample to be tested, I1It is
The electric current of the primary coil, U2It is the induced voltage of the secondary coil.
Preferably, the revised hysteresis model are as follows:
DM/dH=((1-C) (M_an-M_irr))/(δ k- α (M_an-M_irr))+C (dM_an)/dH；
Wherein:
M_an=M_S [coth ((H+ α M)/a)-a/ (H+ α M)]；
M=M_irr+C (M_an-M_irr)；
μ _ 0 B=(M+H)；
Wherein, M is the intensity of magnetization, and C is reversible magnetization coefficient, and M_an is the anhysteretic intensity of magnetization, and M_irr is irreversible magnetic
Change intensity, δ is magnetizing parameters, and as dH/dt >=0, δ=1, as dH/dt < 0, δ=- 1, k is neticdomain wall coefficient of concentration, and a is
Anhysteretic susceptibility, α are main field component, and M_S is saturation magnetization, and μ _ 0 is space permeability.
Preferably, the parameter includes anhysteretic susceptibility, neticdomain wall coefficient of concentration and main field component.
Preferably, further include step S10: the parameter is expressed as by stress and magnetic induction intensity according to the optimal value
Function.
Preferably, the S10 further comprises step:
Determine current anhysteretic susceptibility a when a fixed magnetic induction intensity0, current neticdomain wall coefficient of concentration k0With work as
Preceding main field component α0；
Obtain the current anhysteretic susceptibility a0, the current neticdomain wall coefficient of concentration k0With the current main field
Component α0Normalized parameter:
Anhysteretic susceptibility a, the neticdomain wall coefficient of concentration k and the main field component α are expressed as stress and magnetic
The function of induction:
Each coefficient is tested according to practical magnetism testing and is determined in formula (3).
The present invention due to use above technical scheme, make it have it is following the utility model has the advantages that
1, the present invention provides the quick methods for determining hysteresis model parameter, can determine different residual stress, different peaks
It is worth the relevant information of the magnetic hystersis loss and B-H loop under magnetic induction intensity, provides guidance for further control iron loss.Theoretical mould
Type and engineering is practical combines, it is convenient using MATLAB program solution, it is easy to operate.
2, the present invention considers influence of the magnetic induction intensity variation to magnetic hysteresis parameters, is more in line with different inside motor stator
The different actual condition of the magnetic induction intensity of position, can be determined more accurately magnetic hystersis loss.
3, the present invention considers drawing, compression influences the asymmetry of magnetic hystersis loss.
Detailed description of the invention
Fig. 1 is the flow chart of the JA hysteresis model parameter identification method containing stress of the embodiment of the present invention；
Fig. 2 is the structural schematic diagram of the magnetism testing device containing elastic stress of the embodiment of the present invention；
Fig. 3 is the structural schematic diagram of the sample to be tested of the embodiment of the present invention；
Fig. 4 is certain grade electric steel magnetism testing BH curve of the embodiment of the present invention；
Fig. 5~Fig. 7 is each parameter distribution figure of hysteresis model under certain grade electric steel 1.5T of the embodiment of the present invention；
Fig. 8~Figure 10 is certain grade electric steel hysteresis model parameter of the embodiment of the present invention with Magnetic Induction Density Distribution figure.
Specific embodiment
Below according to attached FIG. 1 to FIG. 8, presently preferred embodiments of the present invention is provided, and is described in detail, is enabled more preferable geographical
Solve function of the invention, feature.
Please refer to FIG. 1 to FIG. 3, a kind of JA hysteresis model parameter identification method containing stress of the embodiment of the present invention, including
Step:
S1: the one magnetism testing device 2 containing elastic stress of setting；Magnetism testing device 2 containing elastic stress includes:
One first sample fixing end 21, one second sample fixing end 22, a Stress Control device 23, a upper magnet yoke 241, a lower yoke
242, an insulation wire frame 25 and a stress display device 26；First sample fixing end 21 and the second sample fixing end 22 are along first party
To laying；And first sample fixing end 21 can along first direction reciprocally with Stress Control device 23 be sequentially connected；Upper magnetic
Yoke 241 and about 242 lower yoke are opposed to be set between the first sample fixing end 21 and the second sample fixing end 22；Stress is aobvious
Showing device 26 connects Stress Control device 23；Two opposed faces of insulation wire frame 25 are wound with a primary coil and a secondary wire respectively
Circle；In the present embodiment, Stress Control device 23 uses air compressor.
S2: a sample 1 to be tested is installed on the magnetism testing device 2 containing elastic stress；
S3: alternating current is led to primary coil；
S4: test sample 1 is treated by Stress Control device 23 and loads compression or tensile stress；
S5: induced current and induced voltage in acquisition secondary coil；
S6: calculate obtain experiment magnetic field strength and experiment magnetic induction intensity, and according to experiment magnetic field by force and experiment magnetic induction
Strength co-mputation obtains experiment magnetic hystersis loss value；
In S6 step, calculates and obtain in fact according to experiment magnetic field strength formula (1) and experiment magnetic induction intensity formula (2)
Test magnetic field strength and experiment magnetic induction intensity:
Wherein, H is experiment magnetic field strength, and B is experiment magnetic induction intensity；N1It is the number of turns of primary coil, L is to be measured has a try
The length of equivalent magnetic circuit, N in sample 12It is the number of turns of secondary coil, S is the cross-sectional area of sample 1 to be tested, I1It is primary coil
Electric current, U2It is the induced voltage of secondary coil.
S7: one revised hysteresis model is solved using MATLAB software, parameter is continued to optimize, acquires a peak value
Magnetic induction intensity predicted value, a peak field prediction of strength value and a magnetic hystersis loss predicted value；
Wherein, revised hysteresis model are as follows:
DM/dH=((1-C) (M_an-M_irr))/(δ k- α (M_an-M_irr))+C (dM_an)/dH；
Wherein:
M_an=M_S [coth ((H+ α M)/a)-a/ (H+ α M)]；
M=M_irr+C (M_an-M_irr)；
μ _ 0 B=(M+H)；
Wherein, M is the intensity of magnetization, and C is reversible magnetization coefficient, and M_an is the anhysteretic intensity of magnetization, and M_irr is irreversible magnetic
Change intensity, δ is parameter (δ=1, when dH/dt >=0；δ=- 1, as dH/dt < 0), k is neticdomain wall coefficient of concentration, and a is anhysteretic
Susceptibility, α are main field component, and M_S is saturation magnetization, and μ _ 0 is space permeability.
S8: by peak magnetic induction prediction of strength value, peak field prediction of strength value and magnetic hystersis loss predicted value and experiment magnetic
Field intensity, experiment magnetic induction intensity and experiment magnetic hystersis loss value are matched；
S9: parameter when using matching error minimum is as optimal value.
Wherein, parameter includes anhysteretic susceptibility, neticdomain wall coefficient of concentration and main field component.
Further include step S10 in the present embodiment: parameter is expressed as to the letter of stress and magnetic induction intensity according to optimal value
Number.
S10 further comprises step:
Determine current anhysteretic susceptibility a when a fixed magnetic induction intensity0, current neticdomain wall coefficient of concentration k0With work as
Preceding main field component α0；
Obtain current anhysteretic susceptibility a0, current neticdomain wall coefficient of concentration k0With current main field component α0Normalizing
Change parameter:
Anhysteretic susceptibility a, neticdomain wall coefficient of concentration k and main field component α are expressed as stress and magnetic induction intensity
Function:
Formula (3) each coefficient is tested according to the magnetism testing of actual specific grade electric steel and is determined.
Such as: a length of a, width b are used, with a thickness of the sample to be tested 1 of t；First by sample 1 to be tested from the second sample
In the insert port insertion insulation wire frame 25 of fixing end 22；In the present embodiment, the length of sample 1 to be tested is insulation 25 length of wire frame
Twice, facilitate 1 end of sample to be tested grips and applies pressure.12 internal voids of wire frame insulate having a size of (a/2+
Δ a) × (b+ Δ b) × (t+ Δ t), wherein Δ is small amount, and Δ a=Δ b=Δ t=1mm respectively represents length direction, width
The smaller length for spending direction, thickness direction, primarily to ensuring that the magnetic line of force passes through as far as possible inside sample 1 to be tested.Insulated wire
Two opposed faces of frame 25 are wound with a primary coil and a secondary coil respectively, exciting current are passed through in primary coil, according to electricity
Law of magnetic induction generates the magnetic line of force inside insulation wire frame 25, due to sample 1 to be tested and upper magnet yoke 241 and 242 groups of lower yoke
At circuit magnetic resistance it is minimum, so magnetic line of force almost all passes through inside sample 1 to be tested.Secondary coil is also wrapped in insulated wire
Above frame 25, for detecting the induced voltage generated.Sample 1 to be tested be put into insulation wire frame 25 after, by upper magnet yoke 241 put down with
1 surface of sample contact to be tested, 1 lower surface of sample to be tested is contacted with lower yoke 242.Upper magnet yoke 241 and lower yoke 242 are pair
Claim structure, purpose is also for ensuring that the external magnetic field formed inside sample 1 to be tested is uniform magnetic field, and upper magnet yoke 241 is by leading
Column 27 can move back and forth along the vertical direction, and guide post 27 can connect on a lifting device (not shown)；Guide post 27 is for limiting
The direction that upper magnet yoke 241 processed glides.Upper magnet yoke 241 and lower yoke 242 are in addition to that can form closed magnetic circuit with sample 1 to be tested
Outside, moreover it is possible to play the role of preventing unstability.First sample fixing end 21 and the second sample fixing end 22 are by the two of sample 1 to be tested
End is clamped by rotational fastener knob 28, is then driven the first sample fixing end 21 by air compressor and is generated displacement,
To generate the stress along 1 length direction of sample to be tested, the size of stress can directly read in seif-citing rate display device 26 and obtain
, and magnetism testing can be carried out when keeping stress constant；In the present embodiment, air is can be used in stress display device 26
The included stress display instrument of compressor.The induced voltage that the exciting current and secondary coil of test generate inputs computer software
End, is obtained by calculation magnetic field strength and magnetic induction intensity, obtains B-H loop and magnetic hystersis loss, and the experiment B-H of the present embodiment is bent
Line is as shown in Figure 4.
It is theoretical to be then based on modified JA magnetic hysteresis, is solved using MATLAB software, continuous Optimal Parameters a, k, α,
Solution obtains magnetic field strength H and magnetic induction density B, and the B-H loop predicted.Magnetic hysteresis damage is calculated with the area of B-H loop
The predicted value of consumption.The predicted value of peak magnetic induction intensity, peak field intensity, magnetic hystersis loss is matched with experiment value,
Parameter a, k, α when with error minimum are optimal value.
In the present embodiment, firstly, determining BmThe distribution of parameter when=1.5T, B-H relationship when according to experiment, obtains 1.5T
When magnetic field strength H with residual stress σ variation, it is theoretical based on modified JA magnetic hysteresis, parameter current a is calculated0,k0,α0With answering
The distribution of power, as shown in Fig. 5~Fig. 7.The distribution of parameter under tension and compression stress is obtained by fitting, parameter is expressed as to point of stress
Section function:
a0=-0.75 σ+70.21 (σ < 0)
a0=-0.40 σ+70.21 (σ >=0)
k0=-1.07 σ+50.03 (σ < 0)
k0=-0.43 σ+50.03 (σ >=0)
α0=0.5 (σ≤- 15MPa)
α0=0.067 σ+1.5 (- 15MPa < σ < 0)
α0=1.5 (σ >=0)
As shown in Fig. 8~Figure 10, be under different stress, magnetic hysteresis parameters with magnetic induction intensity variation.Under same stress, magnetic
Stagnant parameter and magnetic induction intensity are in preferable linear relationship, by handling parameter normalizationCurve is oblique
Rate is expressed as the function of stress.In conjunction with front parameter current a0,k0,α0Distribution, can solve to obtain different stress, different magnetic
Parameter distribution under induction, to solve B-H relationship and magnetic hystersis loss.Parameter expression are as follows:
A=a0(-0.0019σ2+0.0428σ+0.2231)(B-1.5)+1；
K=k0(-0.0002σ2+0.0122σ+0.5846)(B-1.5)+1；
α=α0(σ < 0)；
The α of α=0.55640(B-1.5)+1 (σ≥0)；
Sample result shows that revised JA magnetic hysteresis theory in conjunction with mathematical software, can be passed through determination by the present invention
Parameter is expressed as stress and magnetic induction intensity with the changing rule of magnetic induction intensity by parameter distribution and normalized parameter when 1.5T
Function, preferably calculate BH curve and magnetic hystersis loss, for further control magnetic hystersis loss Assistance And Instruction is provided.
The present invention has been described in detail with reference to the accompanying drawings, those skilled in the art can be according to upper
It states and bright many variations example is made to the present invention.Thus, certain details in embodiment should not constitute limitation of the invention, this
Invention will be using the range that the appended claims define as protection scope of the present invention.
Claims (6)
1. a kind of JA hysteresis model parameter identification method containing stress, comprising steps of
S1: the one magnetism testing device containing elastic stress of setting；The magnetism testing device containing elastic stress includes: one
First sample fixing end, one second sample fixing end, a Stress Control device, a upper magnet yoke, a lower yoke, one insulation wire frame and
One stress display device；The first sample fixing end and the second sample fixing end are laid along first direction；And described
One sample fixing end can reciprocally be sequentially connected with the Stress Control device along the first direction；The upper magnet yoke and
The lower yoke is opposed to be set between the first sample fixing end and the second sample fixing end up and down；The stress
Display device connects the Stress Control device；Two opposed faces of the insulation wire frame are wound with a primary coil and primary respectively
Grade coil；
S2: a sample to be tested is installed on the magnetism testing device containing elastic stress；
S3: alternating current is led to the primary coil；
S4: compression or tensile stress are loaded to the sample to be tested by the Stress Control device；
S5: the induced current and induced voltage in the secondary coil are acquired；
S6: calculating and obtain experiment magnetic field strength and experiment magnetic induction intensity, and strong with the experiment magnetic according to the experiment magnetic field
Induction, which calculates, obtains experiment magnetic hystersis loss value；
S7: one revised hysteresis model is solved using MATLAB software, parameter is continued to optimize, acquires a peak induction
Answer prediction of strength value, a peak field prediction of strength value and a magnetic hystersis loss predicted value；
S8: by the peak magnetic induction prediction of strength value, the peak field prediction of strength value and the magnetic hystersis loss predicted value
It is matched with the experiment magnetic field strength, the experiment magnetic induction intensity and the experiment magnetic hystersis loss value；
S9: parameter when using matching error minimum is as optimal value.
2. the JA hysteresis model parameter identification method according to claim 1 containing stress, which is characterized in that the S6 step
In, it calculates according to experiment magnetic field strength formula (1) and experiment magnetic induction intensity formula (2) and obtains the experiment magnetic field strength
With the experiment magnetic induction intensity:
Wherein, H is experiment magnetic field strength, and B is experiment magnetic induction intensity；N1It is the number of turns of the primary coil, L is described to be measured
It has a try the length of equivalent magnetic circuit in sample, N2It is the number of turns of secondary coil, S is the cross-sectional area of the sample to be tested, I1It is described
The electric current of primary coil, U2It is the induced voltage of the secondary coil.
3. the JA hysteresis model parameter identification method according to claim 2 containing stress, which is characterized in that after the amendment
Hysteresis model are as follows:
DM/dH=((1-C) (M_an-M_irr))/(δ k- α (M_an-M_irr))+C (dM_an)/dH；
Wherein:
M_an=M_S [coth ((H+ α M)/a)-a/ (H+ α M)]；
M=M_irr+C (M_an-M_irr)；
μ _ 0 B=(M+H)；
Wherein, M is the intensity of magnetization, and C is reversible magnetization coefficient, and M_an is the anhysteretic intensity of magnetization, and M_irr is that irreversible magnetization is strong
Degree, δ is magnetizing parameters, and as dH/dt >=0, δ=1, as dH/dt < 0, δ=- 1, k is neticdomain wall coefficient of concentration, and a is non-magnetic
Stagnant susceptibility, α are main field component, and M_S is saturation magnetization, and μ _ 0 is space permeability.
4. the JA hysteresis model parameter identification method according to claim 3 containing stress, which is characterized in that the parameter packet
Include anhysteretic susceptibility, neticdomain wall coefficient of concentration and main field component.
5. the JA hysteresis model parameter identification method according to claim 4 containing stress, which is characterized in that further include step
S10: the parameter is expressed as to the function of stress and magnetic induction intensity according to the optimal value.
6. the JA hysteresis model parameter identification method according to claim 5 containing stress, which is characterized in that the S10 into
One step comprising steps of
Determine current anhysteretic susceptibility a when a fixed magnetic induction intensity0, current neticdomain wall coefficient of concentration k0With current master
Want field component α0；
Obtain the current anhysteretic susceptibility a0, the current neticdomain wall coefficient of concentration k0With the current main field component
α0Normalized parameter:
The parameter is expressed as to the function of stress and magnetic induction intensity:
Each coefficient is tested according to practical magnetism testing and is determined in formula (3).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810426937.3A CN108519569B (en) | 2018-05-07 | 2018-05-07 | JA hysteresis model parameter identification method containing stress |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810426937.3A CN108519569B (en) | 2018-05-07 | 2018-05-07 | JA hysteresis model parameter identification method containing stress |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108519569A CN108519569A (en) | 2018-09-11 |
| CN108519569B true CN108519569B (en) | 2019-07-16 |
Family
ID=63429891
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810426937.3A Active CN108519569B (en) | 2018-05-07 | 2018-05-07 | JA hysteresis model parameter identification method containing stress |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108519569B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109086561B (en) * | 2018-10-08 | 2019-09-24 | 重庆邮电大学 | A kind of meter and anisotropic JA hysteresis model parameter extracting method |
| CN109444776A (en) * | 2018-11-05 | 2019-03-08 | 国网江苏省电力有限公司电力科学研究院 | Three-phase full-bridge converter iron core remanent magnetism measuring method, system and storage medium |
| CN112967862B (en) * | 2021-03-23 | 2023-02-07 | 西安热工研究院有限公司 | J-A model-based three-phase transformer demagnetization method |
| CN113158432B (en) * | 2021-03-25 | 2023-03-21 | 潍柴动力股份有限公司 | Stress and strain-containing JA hysteresis model parameter identification method |
| CN114114103A (en) * | 2021-11-10 | 2022-03-01 | 国家石油天然气管网集团有限公司 | Material force magnetic characteristic parameter determination method for strong and weak magnetic detection |
| CN114236433B (en) * | 2021-11-23 | 2022-08-30 | 浙江大学 | Ferromagnetic material magnetic induction intensity online estimation method |
| CN114384449B (en) * | 2022-03-22 | 2022-07-26 | 河北工业大学 | Three-dimensional space vector magnetism measuring device and method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104991994A (en) * | 2015-06-09 | 2015-10-21 | 国网天津市电力公司 | Transformer field current simulation method based on J-A magnetic hysteresis model |
| CN107239643A (en) * | 2017-07-24 | 2017-10-10 | 安徽理工大学 | The parameter identification apparatus and method of super-magnetostrictive drive magnetic hysteresis nonlinear model |
| CN107273573A (en) * | 2017-05-16 | 2017-10-20 | 电子科技大学 | A kind of leakage field emulation mode based on circular current |
| CN107656221A (en) * | 2017-09-08 | 2018-02-02 | 国网江苏省电力公司电力科学研究院 | A kind of transformer core remanent magnetism evaluation method based on minor loop's slope |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5556882B2 (en) * | 2010-03-18 | 2014-07-23 | 富士通株式会社 | Magnetic body simulation method and program |
-
2018
- 2018-05-07 CN CN201810426937.3A patent/CN108519569B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104991994A (en) * | 2015-06-09 | 2015-10-21 | 国网天津市电力公司 | Transformer field current simulation method based on J-A magnetic hysteresis model |
| CN107273573A (en) * | 2017-05-16 | 2017-10-20 | 电子科技大学 | A kind of leakage field emulation mode based on circular current |
| CN107239643A (en) * | 2017-07-24 | 2017-10-10 | 安徽理工大学 | The parameter identification apparatus and method of super-magnetostrictive drive magnetic hysteresis nonlinear model |
| CN107656221A (en) * | 2017-09-08 | 2018-02-02 | 国网江苏省电力公司电力科学研究院 | A kind of transformer core remanent magnetism evaluation method based on minor loop's slope |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108519569A (en) | 2018-09-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108519569B (en) | JA hysteresis model parameter identification method containing stress | |
| US4528856A (en) | Eddy current stress-strain gauge | |
| CN103238064B (en) | Depth of quenching assay method and depth of quenching determinator | |
| CN107085192A (en) | A kind of method and its device that ferromagnetic material hysteresis curve is measured in open-flux path | |
| Zhang et al. | Identification of the ferromagnetic hysteresis simulation parameters using classic non-destructive testing equipment | |
| de La Barrière et al. | A simple compensation method for the accurate measurement of magnetic losses with a single strip tester | |
| CN109388883A (en) | Magnetic flux towards EMTP-current mode JA magnetic hysteresis inductance acquisition methods | |
| Hamzehbahmani | Static hysteresis modeling for grain-oriented electrical steels based on the phenomenological concepts of energy loss mechanism | |
| CN208399665U (en) | A kind of vertical three-phase flexibility wideband rotary magnetic characteristic measuring system | |
| CN106125021B (en) | The measurement method of permeability magnetic material characteristic under a kind of quadrature bias magnetic field | |
| CN105548668B (en) | A kind of method for measuring superconductor critical current density | |
| Deng et al. | Feasibility study of determining axial stress in ferromagnetic bars using reciprocal amplitude of initial differential susceptibility obtained from static magnetization by permanent magnets | |
| Fiorillo et al. | Measuring the hysteresis loop of permanent magnets with the pulsed field magnetometer | |
| CN104375099A (en) | Material detection probe based on initial permeability method | |
| Ravaioli et al. | Fast method to quantify the collective magnetization in superconducting magnets | |
| Ducharne et al. | Electrical steel dynamic behavior quantitated by inductance spectroscopy: Toward prediction of magnetic losses | |
| JP4192333B2 (en) | Method for measuring transformation layer thickness of steel | |
| Soinski | Demagnetization effect of rectangular and ring-shaped samples made of electrical sheets placed in a stationary magnetic field | |
| Zakharov et al. | Evaluating the structure of a ferromagnetic material based on magnetic-field strength between the poles of an attached two-pole magnetizing device | |
| Chetangny et al. | Eddy-Current Loss in a Conductive Material Inserted into a U-Cored Electromagnet Device | |
| CN104697430B (en) | Method and system for measuring steel plate thickness based on dynamic magnetic conductivity feature | |
| Zhang et al. | Static Hysteresis Process Modeling Based on Modified Tellinen Hysteresis Model | |
| Li et al. | A modified jiles-atherton model for estimating the iron loss of electrical steel considering DC bias | |
| CN112731227B (en) | Residual magnetism detection method for current transformer for protection | |
| Charubin et al. | Automatic measurement station for ferrite materials testing |
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 |