CN107037383B - A kind of magnetic domain Deflection Model method for building up and device - Google Patents

A kind of magnetic domain Deflection Model method for building up and device Download PDF

Info

Publication number
CN107037383B
CN107037383B CN201710313713.7A CN201710313713A CN107037383B CN 107037383 B CN107037383 B CN 107037383B CN 201710313713 A CN201710313713 A CN 201710313713A CN 107037383 B CN107037383 B CN 107037383B
Authority
CN
China
Prior art keywords
formula
scheduled
magnetization
stress
magnetic
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.)
Expired - Fee Related
Application number
CN201710313713.7A
Other languages
Chinese (zh)
Other versions
CN107037383A (en
Inventor
王亚文
谷爱昱
严柏平
江伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong University of Technology
Original Assignee
Guangdong University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN201710313713.7A priority Critical patent/CN107037383B/en
Publication of CN107037383A publication Critical patent/CN107037383A/en
Application granted granted Critical
Publication of CN107037383B publication Critical patent/CN107037383B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/12Measuring magnetic properties of articles or specimens of solids or fluids
    • G01R33/18Measuring magnetostrictive properties

Abstract

The embodiment of the invention discloses a kind of magnetic domain Deflection Model method for building up and devices, the magnetization characteristic that load acts on lower material can preferably be described for solving original magnetic domain angular deflection model, but model is simple to the description of external macroscopic properties, solution procedure is complicated, and model needs the technical issues of carrying out parameters revision.A kind of magnetic domain Deflection Model method for building up of the embodiment of the present invention includes: S1: under scheduled stress and magnetic fields, obtaining the scheduled free energy formula of single domain particle in giant magnetostrictive material;S2: the intensity of magnetization formula on scheduled direction is obtained by scheduled free energy formula;S3: the derivative operation to stress and magnetic field is carried out to the intensity of magnetization formula respectively, obtains corresponding partial derivative;S4: the partial derivative is substituted into scheduled magnetisation density increment formula, obtains scheduled magnetic field incremental value, obtains the calculated result of magnetisation density increment.

Description

A kind of magnetic domain Deflection Model method for building up and device
Technical field
The present invention relates to magnetic domain angular deflection model field more particularly to a kind of magnetic domain Deflection Model method for building up and dress It sets.
Background technique
Giant magnetostrictive material has complicated machine magnetic parameter coupled characteristic, and super mangneto is stretched under stress and magnetic field load effect The machine magnetic coupling properties of compression material can be studied by the deflection of magnetic domain angle.Take a single domain in giant magnetostrictive material Research object of the grain as magnetic domain Deflection Model (S-W model).
The effect of stress and magnetic field load will make magnetization M in single domain particle that the deflection of angle occur, and then influence The magnetization characteristic of giant magnetostrictive material.By the solution of single domain particle free energy Nonlinear System of Equations extreme-value problem in material, The rule of magnetic domain angular deflection in giant magnetostrictive material can be studied, and then describes magnetization and the machine magnetic coupling properties of material.
Summary of the invention
The embodiment of the invention provides a kind of magnetic domain Deflection Model method for building up and devices, solve original magnetic domain angle Deflection Model can preferably describe the magnetization characteristic that load acts on lower material, but model is simple to the description of external macroscopic properties, Solution procedure is complicated, and model needs the technical issues of carrying out parameters revision.
A kind of magnetic domain Deflection Model method for building up provided in an embodiment of the present invention, comprising:
S1: under scheduled stress and magnetic fields, the scheduled freedom of single domain particle in giant magnetostrictive material is obtained It can formula;
S2: the intensity of magnetization formula on scheduled direction is obtained by scheduled free energy formula;
S3: the derivative operation to stress and magnetic field is carried out to the intensity of magnetization formula respectively, obtains corresponding local derviation letter Number;
S4: the partial derivative is substituted into scheduled magnetisation density increment formula, scheduled magnetic field incremental value is obtained, obtains Take the calculated result of magnetisation density increment.
Preferably, scheduled free energy formula are as follows:
Wherein, H is magnetic field;Ms is saturation magnetization;K1And K2For cubic crystal magnetic anisotropy constant;σ is stress Amount;λ100And λ111For magnetostriction parameter;For field direction;For stress direction,For the direction of magnetization;μ0 For magnetic conductivity.
Preferably, the step S2 is specifically included:
Energy of magnetization formula, magnetocrystalline anisotropy energy formula, stress anisotropy energy are obtained by scheduled free energy formula Formula is obtained pre- by the energy of magnetization formula, the magnetocrystalline anisotropy energy formula, the stress anisotropy energy formula Fixed magnetization intensity vectorIntensity of magnetization formula on direction:
Wherein, σ is scheduled stress;H is scheduled magnetic field;For scheduled magnetization intensity vector;For The direction of magnetization;EfieldFor the energy of magnetization;EanFor magnetocrystalline anisotropy energy;EσFor stress anisotropy energy;K is normaliztion constant;ω For Potential distribution parameter.
Preferably, scheduled magnetisation density increment formula are as follows:
Wherein, Δ B, Δ H, Δ σ are respectively magnetisation density increment, magnetic field increment and stress increment.
Preferably, magnetisation density is equal to the product of the intensity of magnetization and relative permeability.
A kind of magnetic domain Deflection Model provided in the embodiment of the present invention establishes device, comprising:
First acquisition unit, for obtaining single domain in giant magnetostrictive material under scheduled stress and magnetic fields The scheduled free energy formula of grain;
Second acquisition unit, for obtaining the intensity of magnetization formula on scheduled direction by scheduled free energy formula;
Local derviation unit, for carrying out the derivative operation to stress and magnetic field, acquisition pair respectively to the intensity of magnetization formula The partial derivative answered;
Third acquiring unit obtains predetermined for substituting into the partial derivative in scheduled magnetisation density increment formula Magnetic field incremental value, obtain magnetisation density increment calculated result.
Preferably, scheduled free energy formula are as follows:
Wherein, H is magnetic field;Ms is saturation magnetization;K1And K2For cubic crystal magnetic anisotropy constant;σ is stress Amount;λ100And λ111For magnetostriction parameter;For field direction;For stress direction,For the direction of magnetization;μ0 For magnetic conductivity.
Preferably, the second acquisition unit specifically includes:
Subelement is obtained, is specifically used for obtaining energy of magnetization formula, magnetocrystalline anisotropy energy by scheduled free energy formula Formula, stress anisotropy energy formula, it is each by the energy of magnetization formula, the magnetocrystalline anisotropy energy formula, the stress Anisotropy energy formula, obtains scheduled magnetization intensity vectorIntensity of magnetization formula on direction:
Wherein, σ is scheduled stress;H is scheduled magnetic field;For scheduled magnetization intensity vector;For The direction of magnetization;EfieldFor the energy of magnetization;EanFor magnetocrystalline anisotropy energy;EσFor stress anisotropy energy;K is normaliztion constant;ω For Potential distribution parameter.
Preferably, scheduled magnetisation density increment formula are as follows:
Wherein, Δ B, Δ H, Δ σ are respectively magnetisation density increment, magnetic field increment and stress increment.
Preferably, magnetisation density is equal to the product of the intensity of magnetization and relative permeability.
As can be seen from the above technical solutions, the embodiment of the present invention has the advantage that
A kind of a kind of magnetic domain Deflection Model method for building up and device provided in the embodiment of the present invention, wherein magnetic domain deflection Method for establishing model includes: S1: under scheduled stress and magnetic fields, single domain particle is pre- in acquisition giant magnetostrictive material Fixed free energy formula;S2: the intensity of magnetization formula on scheduled direction is obtained by scheduled free energy formula;S3: to institute It states intensity of magnetization formula and carries out derivative operation to stress and magnetic field respectively, obtain corresponding partial derivative;S4: by the local derviation Function substitutes into scheduled magnetisation density increment formula, obtains scheduled magnetic field incremental value, obtains the calculating of magnetisation density increment As a result.In the present embodiment, by the solution of single domain particle free energy Nonlinear System of Equations extreme-value problem in material, simplification was solved External macroscopic properties is described in journey, and solving original magnetic domain angular deflection model can preferably describe under load effect The magnetization characteristic of material, but model is simple to the description of external macroscopic properties, and solution procedure is complicated, and model needs to carry out parameter The technical issues of amendment.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this Some embodiments of invention without any creative labor, may be used also for those of ordinary skill in the art To obtain other attached drawings according to these attached drawings.
A kind of process of one embodiment of the magnetic domain Deflection Model method for building up provided in Fig. 1 embodiment of the present invention is illustrated Figure;
A kind of magnetic domain Deflection Model provided in Fig. 2 embodiment of the present invention establishes the structural representation of one embodiment of device Figure;
The structure for another embodiment that a kind of magnetic domain Deflection Model provided in Fig. 3 embodiment of the present invention establishes device is shown It is intended to;
Magnetization curve schematic diagram under a kind of uniaxial mechanical stress provided in Fig. 4 embodiment of the present invention;
The schematic diagram of calculation result of the magnetic field partial differential d1 provided in Fig. 5 embodiment of the present invention a kind of;
The two-dimentional cloud atlas result schematic diagram of the magnetic field partial differential d1 provided in Fig. 6 embodiment of the present invention a kind of;
The schematic diagram of calculation result of the compression partial differential d2 provided in Fig. 7 embodiment of the present invention a kind of;
The two-dimentional cloud atlas result schematic diagram of the compression partial differential d2 provided in Fig. 8 embodiment of the present invention a kind of;
A kind of schematic diagram of calculation result of the flux density variable quantity provided in Fig. 9 embodiment of the present invention;
A kind of calculating two dimension cloud atlas result schematic diagram of the flux density variable quantity provided in Figure 10 embodiment of the present invention;
A kind of abscissa provided in Figure 11 embodiment of the present invention is the ratio of magnetic field increment and stress increment, and ordinate is The change of magnetic field strength amount schematic diagram of stress increment.
Specific embodiment
The embodiment of the invention provides a kind of magnetic domain Deflection Model method for building up and devices, for solving original magnetic domain angle Degree Deflection Model can preferably describe the magnetization characteristic that load acts on lower material, but description letter of the model to external macroscopic properties Single, solution procedure is complicated, and model needs the technical issues of carrying out parameters revision.
In order to make the invention's purpose, features and advantages of the invention more obvious and easy to understand, below in conjunction with the present invention Attached drawing in embodiment, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that disclosed below Embodiment be only a part of the embodiment of the present invention, and not all embodiment.Based on the embodiments of the present invention, this field Those of ordinary skill's all other embodiment obtained without making creative work, belongs to protection of the present invention Range.
Referring to Fig. 1, a kind of one embodiment packet of the magnetic domain Deflection Model method for building up provided in the embodiment of the present invention It includes:
101, under scheduled stress and magnetic fields, the scheduled freedom of single domain particle in giant magnetostrictive material is obtained It can formula;
102, the intensity of magnetization formula on scheduled direction is obtained by scheduled free energy formula;
103, the derivative operation to stress and magnetic field is carried out to the intensity of magnetization formula respectively, obtains corresponding local derviation letter Number;
104, the partial derivative is substituted into scheduled magnetisation density increment formula, obtains scheduled magnetic field incremental value, Obtain the calculated result of magnetisation density increment.
The above is that a kind of magnetic domain Deflection Model method for building up is described in detail, and will deflect mould to a kind of magnetic domain below The process of type method for building up is described in detail, a kind of magnetic domain Deflection Model method for building up provided in the embodiment of the present invention Another embodiment includes:
201, under scheduled stress and magnetic fields, the scheduled freedom of single domain particle in giant magnetostrictive material is obtained Energy formula, scheduled free energy formula are as follows:
Wherein, H is magnetic field;Ms is saturation magnetization;K1And K2For cubic crystal magnetic anisotropy constant;σ is stress Amount;λ100And λ111For magnetostriction parameter;For field direction;For stress direction,For the direction of magnetization;μ0 For magnetic conductivity;
Giant magnetostrictive material has complicated machine magnetic parameter coupled characteristic, and super mangneto is stretched under stress and magnetic field load effect The machine magnetic coupling properties of compression material can be studied by the deflection of magnetic domain angle.Take a single domain in giant magnetostrictive material Research object of the grain as magnetic domain Deflection Model (S-W model).The effect of stress and magnetic field load will be so that magnetic in single domain particle Change the deflection that angle occurs for intensity M, and then influences the magnetization characteristic of giant magnetostrictive material.It is free by single domain particle in material The solution of energy Nonlinear System of Equations extreme-value problem, can study the rule of magnetic domain angular deflection in giant magnetostrictive material, and then retouch State magnetization and the machine magnetic coupling properties of material.Under external load (stress σ and magnetic field H) effect, anisotropic influence is considered, The free energy of single domain particle may be expressed as:
In formula (1), wherein H is magnetic field;Ms is saturation magnetization;K1And K2For cubic crystal magnetic anisotropy constant;σ For stress tensor;λ100And λ111For magnetostriction parameter;For field direction;For stress direction,For magnetic Change direction;μ0For magnetic conductivity, μ0=4 Ω × 10-7H/m;
E is the energy of magnetization, magnetocrystalline anisotropy energy, the summation of stress anisotropy energy.
The energy of magnetization is equal to the intensity H in magnetic field, the angle between saturation magnetization Ms and magnetization and the magnetic vector of application The product of degree, the formula of the energy of magnetization are as follows:
Magnetocrystalline anisotropy energy is the energy of the per unit volume discharged when magnetization is rotated away from the hard direction of crystals It measures, wherein K1And K2It is cubic crystal magnetic anisotropy constant, the formula of magnetocrystalline anisotropy energy are as follows:
Stress anisotropy can be generated by the interaction of the stress tensor σ and magnetostrictive strain tensor that are applied, Middle λ100And λ111For magnetostriction parameter, the formula of stress anisotropy energy are as follows:
202, energy of magnetization formula, magnetocrystalline anisotropy energy formula, stress are obtained respectively to different by scheduled free energy formula Performance formula is obtained by the energy of magnetization formula, the magnetocrystalline anisotropy energy formula, the stress anisotropy energy formula Take scheduled magnetization intensity vectorIntensity of magnetization formula on direction:
Wherein, σ is scheduled stress;H is scheduled magnetic field;For scheduled magnetization intensity vector;For The direction of magnetization;EfieldFor the energy of magnetization;EanFor magnetocrystalline anisotropy energy;EσFor stress anisotropy energy;K is normaliztion constant;ω For Potential distribution parameter;
DirectionOn superposed magnetization M can by from all possible magnetization orientation in this direction Magnetization contribution is overlapped and obtains:
K normaliztion constant is to make the function for a spy for the constant contained by any section of any nonnegative function The integral for determining section is exactly equal to 1, obtains the magnetization curve under Fig. 4 single shaft mechanical stress, and wherein abscissa is magnetic field H, indulges and sits It is designated as magnetization M, every curve represents the magnetization curve under different stress.
203, the derivative operation to stress and magnetic field is carried out to the intensity of magnetization formula respectively, obtains corresponding local derviation letter Number;
Local derviation is asked to obtain Fig. 5 the intensity of magnetization, 6,7,8.
Wherein B is magnetisation density, the relationship with magnetization M are as follows: B=μ0M, Δ B, Δ H, Δ σ are respectively to magnetize Density increment, magnetic field increment and stress increment, magnetic field partial differential d1 result (result under unit 1kA/m): if Fig. 5, Fig. 6 are magnetic field For the two-dimentional cloud atlas of partial differential d1 as a result, wherein abscissa is magnetic field H, ordinate is stress, compression partial differential d2 result (unit Result under 1MPa): if the two-dimentional cloud atlas that Fig. 7, Fig. 8 are compression partial differential d2 is as a result, abscissa is magnetic field H, ordinate is to answer Power σ.
204, the partial derivative is substituted into scheduled magnetisation density increment formula, obtains scheduled magnetic field incremental value, Obtain the calculated result of magnetisation density increment, scheduled magnetisation density increment formula are as follows:
Wherein, Δ B, Δ H, Δ σ are respectively magnetisation density increment, magnetic field increment and stress increment, and magnetisation density is equal to magnetic Change the product of intensity and relative permeability.
B is magnetisation density, the relationship with magnetization M are as follows: B=μ0M, Δ B, Δ H, Δ σ are respectively magnetisation density Increment, magnetic field increment and stress increment, magnetic field partial differential d1 result (result under unit 1kA/m): if Fig. 5, Fig. 6 are that magnetic field is partially micro- Divide the two-dimentional cloud atlas of d1 as a result, wherein abscissa is magnetic field H, ordinate is stress, compression partial differential d2 result (unit 1MPa Lower result): if the two-dimentional cloud atlas that Fig. 7, Fig. 8 are compression partial differential d2 is as a result, wherein abscissa is magnetic field H, ordinate is to answer Power σ substitutes into the partial derivative acquired in formula 8:
Since the energy generated under stress and magnetic field independent role need to reach 1:1, i.e. Δ b1: Δ b2=1:1.According to the mankind It is -10MPa that walking, which generates stress Δ σ value,.Δ H value in (2.8~12.5km/A) range can be released.It is quasi- first to take one Initial value Δ H=8.8km/A is substituted into formula 8.It is calculated Fig. 9,10, the calculated result of flux density variable quantity: such as Fig. 9, Tu10Wei For the calculating two dimension cloud atlas of flux density variable quantity as a result, wherein abscissa is magnetic field H, ordinate is stress σ.It, can by observing Figure 10 Preliminary value σ0=-5MPa, H0=30000km/A is independent variable calculation formula 8, formula with Δ σ, Δ H to examine result above 8 can be changed into formula 9:
Figure 11 is obtained by calculation formula (9), the abscissa of Figure 11 is the ratio of magnetic field increment and stress incrementOrdinate is stress increment Δ σ, and Figure 11 is shown as Δ σ=- 8MPa,When, the intensity of magnetization reaches Maximum variable quantity.It thus can more than reverse-direction derivation preliminary value be suitable.Therefore σ0=-5MPa, H0=30000km/A is Best value.
In the present embodiment, directionOn integral magnetization M can be by from all possible magnetization orientation Magnetization contribution in this direction is integrated and is obtained:
The present embodiment can take a kind of easier method to realize the integral process in above formula, and θ is enabled to take in 0- π 360 values,180 values are taken in 0-2 π, by eachThe opposite sex under value can be carried out superposition, obtain
WhereinK is assessed under each applied field and stress condition Normaliztion constant, ω is Potential distribution parameter, and solution procedure is easier, formulaFirst item and Section 2, which is respectively described, individually applies stress and magnetic field intensity of magnetization increment achieved to magnetostriction materials, to external macro The description for seeing characteristic is more specific, and the magnetic that load acts on lower material can preferably be described by solving original magnetic domain angular deflection model Change characteristic, but model is simple to the description of external macroscopic properties, solution procedure is complicated, and model needs to carry out the skill of parameters revision Art problem.
Referring to Fig. 2, a kind of magnetic domain Deflection Model provided in the embodiment of the present invention establishes one embodiment packet of device It includes:
First acquisition unit 301, for obtaining single domain in giant magnetostrictive material under scheduled stress and magnetic fields The scheduled free energy formula of particle;
Second acquisition unit 302, for obtaining the public affairs of the intensity of magnetization on scheduled direction by scheduled free energy formula Formula;
Local derviation unit 303 is obtained for carrying out the derivative operation to stress and magnetic field respectively to the intensity of magnetization formula Corresponding partial derivative;
Third acquiring unit 304 obtains pre- for substituting into the partial derivative in scheduled magnetisation density increment formula Fixed magnetic field incremental value obtains the calculated result of magnetisation density increment.
The above is that a kind of each unit that magnetic domain Deflection Model establishes device is described in detail, below will be to a kind of magnetic Each extra cell that farmland Deflection Model establishes device is described in more detail, referring to Fig. 3, provided in the embodiment of the present invention A kind of another embodiment that magnetic domain Deflection Model establishes device includes:
First acquisition unit 401, for obtaining single domain in giant magnetostrictive material under scheduled stress and magnetic fields The scheduled free energy formula of particle, scheduled free energy formula are as follows:
Wherein, H is magnetic field;Ms is saturation magnetization;K1And K2For cubic crystal magnetic anisotropy constant;σ is stress Amount;λ100And λ111For magnetostriction parameter;For field direction;For stress direction,For the direction of magnetization;μ0 For magnetic conductivity, μ0=4 Ω × 10-7H/m;
Second acquisition unit 402, for obtaining the public affairs of the intensity of magnetization on scheduled direction by scheduled free energy formula Formula;
The second acquisition unit 402 specifically includes:
Subelement 4021 is obtained, is specifically used for different by scheduled free energy formula acquisition energy of magnetization formula, magnetocrystalline Performance formula, stress anisotropy energy formula by the energy of magnetization formula, the magnetocrystalline anisotropy energy formula, described are answered Power anisotropy energy formula, obtains scheduled magnetization intensity vectorIntensity of magnetization formula on direction:
Wherein, σ is scheduled stress;H is scheduled magnetic field;For scheduled magnetization intensity vector;For The direction of magnetization;EfieldFor the energy of magnetization;EanFor magnetocrystalline anisotropy energy;EσFor stress anisotropy energy;K is normaliztion constant;ω For Potential distribution parameter.
Local derviation unit 403 is obtained for carrying out the derivative operation to stress and magnetic field respectively to the intensity of magnetization formula Corresponding partial derivative;
Third acquiring unit 404 obtains pre- for substituting into the partial derivative in scheduled magnetisation density increment formula Fixed magnetic field incremental value obtains the calculated result of magnetisation density increment, scheduled magnetisation density increment formula are as follows:
Wherein, Δ B, Δ H, Δ σ are respectively magnetisation density increment, magnetic field increment and stress increment, and magnetisation density is equal to magnetic Change the product of intensity and relative permeability.
It is apparent to those skilled in the art that for convenience and simplicity of description, the system of foregoing description, The specific work process of device and unit, can refer to corresponding processes in the foregoing method embodiment, and details are not described herein.
In several embodiments provided herein, it should be understood that disclosed system, device and method can be with It realizes by another way.For example, the apparatus embodiments described above are merely exemplary, for example, the unit It divides, only a kind of logical function partition, there may be another division manner in actual implementation, such as multiple units or components It can be combined or can be integrated into another system, or some features can be ignored or not executed.Another point, it is shown or The mutual coupling, direct-coupling or communication connection discussed can be through some interfaces, the indirect coupling of device or unit It closes or communicates to connect, can be electrical property, mechanical or other forms.
The unit as illustrated by the separation member may or may not be physically separated, aobvious as unit The component shown may or may not be physical unit, it can and it is in one place, or may be distributed over multiple In network unit.It can select some or all of unit therein according to the actual needs to realize the mesh of this embodiment scheme 's.
It, can also be in addition, the functional units in various embodiments of the present invention may be integrated into one processing unit It is that each unit physically exists alone, can also be integrated in one unit with two or more units.Above-mentioned integrated list Member both can take the form of hardware realization, can also realize in the form of software functional units.
If the integrated unit is realized in the form of SFU software functional unit and sells or use as independent product When, it can store in a computer readable storage medium.Based on this understanding, technical solution of the present invention is substantially The all or part of the part that contributes to existing technology or the technical solution can be in the form of software products in other words It embodies, which is stored in a storage medium, including some instructions are used so that a computer Equipment (can be personal computer, server or the network equipment etc.) executes the complete of each embodiment the method for the present invention Portion or part steps.And storage medium above-mentioned includes: USB flash disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic or disk etc. are various can store journey The medium of sequence code.
The above, the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations;Although referring to before Stating embodiment, invention is explained in detail, those skilled in the art should understand that: it still can be to preceding Technical solution documented by each embodiment is stated to modify or equivalent replacement of some of the technical features;And these It modifies or replaces, the spirit and scope for technical solution of various embodiments of the present invention that it does not separate the essence of the corresponding technical solution.

Claims (6)

1. a kind of magnetic domain Deflection Model method for building up characterized by comprising
S1: under scheduled stress and magnetic fields, the scheduled free energy for obtaining single domain particle in giant magnetostrictive material is public Formula, wherein scheduled free energy formula are as follows:
Wherein, H is magnetic field;Ms is saturation magnetization;K1And K2For cubic crystal magnetic anisotropy constant;σ is stress tensor;λ100 And λ111For magnetostriction parameter;For field direction;For stress direction,For the direction of magnetization;μ0For magnetic conductance Rate;
S2: obtaining energy of magnetization formula, magnetocrystalline anisotropy energy formula, stress anisotropy by scheduled free energy formula can be public Formula is obtained predetermined by the energy of magnetization formula, the magnetocrystalline anisotropy energy formula, the stress anisotropy energy formula Magnetization intensity vectorIntensity of magnetization formula on direction:
Wherein, σ is scheduled stress;H is scheduled magnetic field;For scheduled magnetization intensity vector;For magnetization Direction;EfieldFor the energy of magnetization;EanFor magnetocrystalline anisotropy energy;EσFor stress anisotropy energy;K is normaliztion constant;ω is electricity Bit distribution parameter;
S3: the derivative operation to stress and magnetic field is carried out to the intensity of magnetization formula respectively, obtains corresponding partial derivative;
S4: the partial derivative is substituted into scheduled magnetisation density increment formula, obtains scheduled magnetic field incremental value, obtains magnetic Change the calculated result of density increment.
2. magnetic domain Deflection Model method for building up according to claim 1, which is characterized in that scheduled magnetisation density increment is public Formula are as follows:
Wherein, Δ B, Δ H, Δ σ are respectively magnetisation density increment, magnetic field increment and stress increment.
3. magnetic domain Deflection Model method for building up according to claim 2, which is characterized in that magnetisation density is equal to the intensity of magnetization With the product of relative permeability.
4. a kind of magnetic domain Deflection Model establishes device characterized by comprising
First acquisition unit, under scheduled stress and magnetic fields, obtaining single domain particle in giant magnetostrictive material Scheduled free energy formula;
Wherein, scheduled free energy formula are as follows:
Wherein, H is magnetic field;Ms is saturation magnetization;K1And K2For cubic crystal magnetic anisotropy constant;σ is stress tensor;λ100 And λ111For magnetostriction parameter;For field direction;For stress direction,For the direction of magnetization;μ0For magnetic conductance Rate;
Second acquisition unit, for obtaining the intensity of magnetization formula on scheduled direction by scheduled free energy formula;
The second acquisition unit specifically includes:
Obtain subelement, be specifically used for by scheduled free energy formula acquisition energy of magnetization formula, magnetocrystalline anisotropy energy formula, Stress anisotropy energy formula passes through the energy of magnetization formula, the magnetocrystalline anisotropy energy formula, the stress anisotropy Energy formula, obtains scheduled magnetization intensity vectorIntensity of magnetization formula on direction:
Wherein, σ is scheduled stress;H is scheduled magnetic field;For scheduled magnetization intensity vector;For magnetization Direction;EfieldFor the energy of magnetization;EanFor magnetocrystalline anisotropy energy;EσFor stress anisotropy energy;K is normaliztion constant;ω is electricity Bit distribution parameter;
Local derviation unit obtains corresponding for carrying out the derivative operation to stress and magnetic field respectively to the intensity of magnetization formula Partial derivative;
Third acquiring unit obtains scheduled magnetic for substituting into the partial derivative in scheduled magnetisation density increment formula Field incremental value, obtains the calculated result of magnetisation density increment.
5. magnetic domain Deflection Model according to claim 4 establishes device, which is characterized in that scheduled magnetisation density increment is public Formula are as follows:
Wherein, Δ B, Δ H, Δ σ are respectively magnetisation density increment, magnetic field increment and stress increment.
6. magnetic domain Deflection Model according to claim 5 establishes device, which is characterized in that magnetisation density is equal to the intensity of magnetization With the product of relative permeability.
CN201710313713.7A 2017-05-05 2017-05-05 A kind of magnetic domain Deflection Model method for building up and device Expired - Fee Related CN107037383B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710313713.7A CN107037383B (en) 2017-05-05 2017-05-05 A kind of magnetic domain Deflection Model method for building up and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710313713.7A CN107037383B (en) 2017-05-05 2017-05-05 A kind of magnetic domain Deflection Model method for building up and device

Publications (2)

Publication Number Publication Date
CN107037383A CN107037383A (en) 2017-08-11
CN107037383B true CN107037383B (en) 2019-08-09

Family

ID=59538206

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710313713.7A Expired - Fee Related CN107037383B (en) 2017-05-05 2017-05-05 A kind of magnetic domain Deflection Model method for building up and device

Country Status (1)

Country Link
CN (1) CN107037383B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4094003A (en) * 1976-03-29 1978-06-06 Canadian Patents And Development Limited Sonic magnetic domain sensor
CN1303090A (en) * 1999-11-26 2001-07-11 株式会社日立制作所 Perpendicular magnetic recording medium, magnetic recording device
JP4029785B2 (en) * 2003-06-30 2008-01-09 株式会社ジェイテクト Rolling bearing unit with sensor
CN101852775A (en) * 2010-06-08 2010-10-06 浙江大学 Device for regulating magnetostrictive detection of longitudinal static magnetic field
CN102288673A (en) * 2011-05-10 2011-12-21 中国科学院电工研究所 Quantitative detection method adopting metal magnetic memory

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4094003A (en) * 1976-03-29 1978-06-06 Canadian Patents And Development Limited Sonic magnetic domain sensor
CN1303090A (en) * 1999-11-26 2001-07-11 株式会社日立制作所 Perpendicular magnetic recording medium, magnetic recording device
JP4029785B2 (en) * 2003-06-30 2008-01-09 株式会社ジェイテクト Rolling bearing unit with sensor
CN101852775A (en) * 2010-06-08 2010-10-06 浙江大学 Device for regulating magnetostrictive detection of longitudinal static magnetic field
CN102288673A (en) * 2011-05-10 2011-12-21 中国科学院电工研究所 Quantitative detection method adopting metal magnetic memory

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"Electromagnetic energy harvester by using buried NdFeB";S. Miki, 等;《2012 IEEE 25th International Conference on Micro Electro Mechanical Systems (MEMS)》;20121231;1221-1224 *
"超磁致伸缩执行器磁机耦合模型及自感知应用研究";严柏平;《中国博士学位论文全文数据库 信息科技辑》;20141215;I140-37 *

Also Published As

Publication number Publication date
CN107037383A (en) 2017-08-11

Similar Documents

Publication Publication Date Title
Garcia-Gonzalez Magneto-visco-hyperelasticity for hard-magnetic soft materials: theory and numerical applications
Muhanna Bohr's phenomenon in subordination and bounded harmonic classes
Fukuyama Theory of orbital magnetism of Bloch electrons: Coulomb interactions
Mehar et al. Free vibration and bending behaviour of CNT reinforced composite plate using different shear deformation theory
Delduc et al. On q-deformed symmetries as Poisson–Lie symmetries and application to Yang–Baxter type models
Allahyarov et al. Magnetomechanical response of bilayered magnetic elastomers
Nguyen-Van et al. Geometrically nonlinear analysis of composite plates and shells via a quadrilateral element with good coarse-mesh accuracy
Klevtsov Lowest Landau level on a cone and zeta determinants
Ishiki et al. Kähler structure in the commutative limit of matrix geometry
Ichimura et al. A numerical exploration of signal detector arrangement in a spin-wave reservoir computing device
Wang et al. Enhancing the nonlocal advantage of quantum coherence by local parity-time (PT)-symmetric operation
CN107037383B (en) A kind of magnetic domain Deflection Model method for building up and device
Gran et al. All superalgebras for warped AdS2 and black hole near horizon geometries
Augusiak et al. Checking the optimality of entanglement witnesses: An application to structural physical approximations
Dong et al. Pitchfork bifurcation and circuit implementation of a novel Chen hyper-chaotic system
Daghigh et al. Gravitational and electromagnetic radiation from an electrically charged black hole in general nonlinear electrodynamics
O’Dwyer et al. Investigation of strain hardening effects under in-plane shear of unidirectional composite materials
Masuku et al. De Alfaro, Fubini and Furlan from multi matrix systems
Chirka Orthogonal complex structures in
Mizel Entanglement versus gap, quantum teleportation, and the AKLT model
Yin et al. Quantum correlations dynamics of three-qubit states coupled to an XY spin chain: Role of coupling strengths☆
Zhang et al. A meshfree cell‐based smoothed point interpolation method for solid mechanics problems
Mahapatra et al. Effects of hygrothermal conditions on free vibration behaviour of laminated composite structures
Trallori et al. Excitation spectrum of an antiferromagnetic film with a non-uniform ground state
Fu Canonical bases for modified quantum gln and q-Schur algebras

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190809