CN112464430B - Method and system for determining vibration noise of iron core of extra-high voltage shunt reactor - Google Patents
Method and system for determining vibration noise of iron core of extra-high voltage shunt reactor Download PDFInfo
- Publication number
- CN112464430B CN112464430B CN202011117619.2A CN202011117619A CN112464430B CN 112464430 B CN112464430 B CN 112464430B CN 202011117619 A CN202011117619 A CN 202011117619A CN 112464430 B CN112464430 B CN 112464430B
- Authority
- CN
- China
- Prior art keywords
- iron core
- geometric model
- magnetic
- determining
- core
- 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
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical group [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 title claims abstract description 268
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000001133 acceleration Effects 0.000 claims abstract description 8
- 230000005284 excitation Effects 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims description 81
- 238000004804 winding Methods 0.000 claims description 61
- 230000006698 induction Effects 0.000 claims description 42
- 230000035699 permeability Effects 0.000 claims description 42
- 239000004579 marble Substances 0.000 claims description 36
- 239000011159 matrix material Substances 0.000 claims description 36
- 229910000976 Electrical steel Inorganic materials 0.000 claims description 23
- 229910052742 iron Inorganic materials 0.000 claims description 20
- 238000006073 displacement reaction Methods 0.000 claims description 18
- 230000005855 radiation Effects 0.000 claims description 17
- 239000002184 metal Substances 0.000 claims description 12
- 229910052751 metal Inorganic materials 0.000 claims description 12
- 239000002023 wood Substances 0.000 claims description 11
- 230000005415 magnetization Effects 0.000 claims description 6
- 230000001419 dependent effect Effects 0.000 claims description 2
- 238000004088 simulation Methods 0.000 abstract description 8
- 238000013461 design Methods 0.000 abstract description 6
- 230000009471 action Effects 0.000 abstract description 4
- 230000009467 reduction Effects 0.000 abstract description 4
- 230000000694 effects Effects 0.000 abstract description 3
- 238000003475 lamination Methods 0.000 abstract description 3
- 238000009826 distribution Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 10
- 238000004590 computer program Methods 0.000 description 7
- 238000005259 measurement Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000005457 optimization Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/10—Noise analysis or noise optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Magnetic Ceramics (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
The invention provides a method and a system for determining vibration noise of an iron core of an extra-high voltage shunt reactor. The method and the system are characterized in that a vibration simulation experiment model of the shunt reactor is established, magnetostriction force and electromagnetic acting force of an iron core of the simulation geometric model under the action of a magnetic field are calculated, then an acoustic vibration model is adopted, acceleration on the surface of the iron core is used as excitation, sound field distribution around the iron core is calculated, and average sound pressure level around the geometric model is determined. According to the method and the system, the whole iron core is taken as a whole, the Young modulus of an iron core cake is measured to simulate and realize the lamination structure of the iron core, and the pre-strain of the iron core caused by magnetostriction is determined according to the telescopic relation curve of the measured magnetic field and the magnetostriction quantity, so that the influence of the magnetostriction force and the electromagnetic force of the iron core on the vibration noise of the reactor is more accurately distinguished, and the method and the system have an important effect on the theoretical design of the body noise reduction of the reactor.
Description
Technical Field
The invention relates to the field of electrical equipment, in particular to a method and a system for determining vibration noise of an iron core of an extra-high voltage shunt reactor.
Background
The extra-high voltage shunt reactor is generally connected between the tail end of an extra-high voltage transmission line and the ground, mainly plays a reactive compensation role, is one of important equipment for safe operation of an extra-high voltage transformer substation, and is also one of main noise source equipment in the transformer substation. Because the noise energy of the extra-high voltage shunt reactor is mainly distributed on a 1/3 frequency multiplication band with 100Hz as the central frequency, the noise energy is about 80% of the total sound energy, and the phenomenon of sound wave interference exists. The ultra-high voltage shunt reactor has the noise level higher than that of the transformer with the same capacity by about 10dB due to the attraction of alternating magnetic fields existing between the segmented iron core cakes. The vibration mechanism research of the shunt reactor iron core has important significance for the ultra-high voltage shunt reactor body noise reduction technology.
IN order to reduce the body noise of the reactor, measures such as heightening the enclosing wall and adding a Box-IN device to the reactor are mainly adopted at present to reduce the factory boundary noise [3] of the transformer substation. Although these measures reduce factory boundary noise to some extent, since the reactor acoustic energy is mainly concentrated in the low frequency band, the line spectrum characteristics are obvious, so that the effects of these engineering measures are not expected. Therefore, the noise of the reactor is urgently required to be reduced from the body, although some researches on the noise reduction of the body exist in China, more than more, related measures are adopted from experience to reduce the noise, and in general, manufacturers lack comprehensive and careful researches on the vibration and transmission mechanism of the shunt reactor iron core, so that related optimization design methods are lacked. On the other hand, the accuracy of the simulation calculation of the reactor is still to be further improved, so that the fine design guidance on engineering cannot be satisfied. In addition, because the body volume and the weight of the extra-high voltage shunt reactor are too large and the loaded voltage and current parameters are too high, the actual experimental verification of theoretical simulation and calculation results cannot be carried out by adopting the real-type extra-high voltage shunt reactor. In general, the vibration noise of the shunt reactor is calculated relatively more in the prior art, but the error of the vibration noise is relatively large compared with the actual result, and the optimization design of the low-noise reactor cannot be guided accurately.
Disclosure of Invention
In order to solve the technical problems that in the prior art, the calculation error of the vibration noise of the shunt reactor is large and the optimization design of the low-noise reactor cannot be guided accurately, the invention provides a method for determining the vibration noise of an iron core of an extra-high voltage shunt reactor, which comprises the following steps:
applying voltage or current to windings of a 1/2 geometric model of a pre-established shunt reactor, and according to the current density J and the resistivity rho on the windings, presetting the magnetic permeability mu of an acoustic propagation medium 0 And relative permeability mu of iron core 1 Calculating the magnetic vector A of the iron core and the periphery of the iron core of the geometric model, and determining the magnetic induction intensity B of the iron core and the periphery of the iron core according to the rotation of the magnetic vector A;
according to the magnetic induction intensity B, a magnetic field B of a silicon steel sheet adopted by a predetermined iron core under a certain pressure sat And magnetostriction amount epsilon sat Calculating the pre-strain epsilon of the iron core at any position caused by magnetostriction ms ；
According to the magnetic induction intensity B, the Young modulus E of the iron core is determined in advance, the Poisson ratio alpha and the magnetic resistance rate v of the iron core are set in advance 0 And the pre-strain epsilon of the core ms Calculating electromagnetic acting force F of the iron core;
Calculating displacement u of the iron core during vibration according to the electromagnetic acting force F and a mass matrix [ M ] and a rigidity matrix [ K ] of the geometric model, which are determined in advance;
according to the displacement u of the iron core during vibration, the characteristic impedance ρc of the preset sound propagation medium is set, and the sound radiation coefficient k of the j-th surface of the geometric model is divided in advance j The j-th faceArea S of (2) c And the included angle theta between the jth surface and any point i around the geometric model j Determining the total energy W of the radiated sound source at the point i i ；
According to the total energy W of the radiated sound source at the point i i Distance R from point i to geometric model surface i A preset sound field degree of freedom coefficient Q and reference energy W 0 Determining sound pressure level LP i ；
According to the sound pressure level LP i Indeed the average sound pressure level LP around the geometric model A 。
Further, before loading voltage or current on the windings of the 1/2 geometric model of the pre-established shunt reactor, the method further comprises:
establishing a 1/2 geometric model of the extra-high voltage shunt reactor, wherein the geometric model comprises an iron yoke, an iron core, a winding, a marble Dan Diankuai, a winding supporting insulating plate, a wood cushion block, a shielding layer and a metal pull rod device, wherein the iron core and the iron yoke are modeled according to the whole, and the iron yoke provides a magnetic circuit for the reactor to close magnetic force lines; the gap formed by the iron core and the marble cushion blocks forms the inductance of the reactor; the marble cushion block adopts the principle of equal area, and a plurality of marble cushion blocks of each layer of the practical shunt reactor are replaced by an integral ring-shaped marble cushion block for supporting the iron cores and providing air gaps between the iron cores; the winding supporting insulating plate and the wood cushion block are used for supporting the winding; the shielding layer is used for shielding the magnetic field to generate eddy current induction on the shell; the metal pull rod device is used for providing pretightening force for the iron core column consisting of the iron core cake and the marble cushion block; the specific structure of winding turns in the actual shunt reactor is ignored by the winding, modeling is carried out according to the whole size, the geometric model is divided into T surfaces by grid division, and the acoustic radiation coefficient of the jth surface is k j The area of each surface unit is S c ；
Measuring an iron core of the geometric model by adopting a pulse excitation method, and determining the Young modulus E of the iron core;
iron core for measuring geometric model by magnetization and magnetostriction measuring systemMagnetostriction curve of the adopted silicon steel sheet under a certain pressure, wherein the magnetostriction curve is a magnetic field B sat And magnetostriction amount epsilon sat Is a relationship of (2);
determining a mass matrix [ M ] and a stiffness matrix [ K ] of the geometric model;
setting magnetic permeability mu of sound propagation medium at position of geometric model 0 And the characteristic impedance ρc, the core reluctance rate ν 0 Poisson's ratio alpha and relative permeability mu in an acoustic propagation medium 1 And the sound field degree of freedom coefficient Q and the reference energy W 0 。
Further, a voltage or a current is applied to the windings of the geometric model, and the air permeability mu is determined according to the current density J and the resistivity rho of the windings 0 And relative permeability mu of iron core 1 Calculating the magnetic vector position A of the iron core and the periphery thereof, and determining the magnetic induction intensity B of the iron core and the periphery thereof according to the rotation of the magnetic vector position A, wherein the calculation formula is as follows:
further, according to the magnetic induction intensity B, a magnetic field B of a silicon steel sheet adopted by the iron core is determined in advance under a certain pressure sat And magnetostriction amount epsilon sat Calculating the pre-strain epsilon of the iron core at any position caused by magnetostriction ms The calculation formula is as follows:
ε ms ＝1.5×ε sat ×(B/B sat ) 2 。
further, according to the magnetic induction intensity B, the Young modulus E of the iron core is determined in advance, the Poisson ratio alpha of the iron core is set in advance, and the magnetic resistance rate v of the iron core is set in advance 0 And the pre-strain epsilon of the core ms And calculating the electromagnetic acting force F of the iron core, wherein the calculation formula is as follows:
where n represents the vector operator in the normal direction of the magnetic field.
Further, according to the electromagnetic acting force F, a mass matrix [ M ] and a stiffness matrix [ K ] of the geometric model are determined in advance, and the displacement u of the iron core during vibration is calculated, wherein the calculation formula is as follows:
further, according to the acceleration v of the iron core during vibration, the characteristic impedance ρc of the acoustic propagation medium is preset, and the acoustic radiation coefficient k of the j-th surface of the geometric model is divided in advance j Area S of jth face c And the included angle theta between the jth surface and any point i around the geometric model j Determining the total energy W of the radiated sound source at the point i i The calculation formula is as follows:
wherein T represents the number of meshes determined when the geometric model is subjected to meshing.
Further, according to the total energy W of the radiated sound source at the point i i Distance R from point i to geometric model surface i A preset sound field degree of freedom coefficient Q and reference energy W 0 Determining sound pressure level LP i The calculation formula is as follows:
Lp i ＝10lg(W i /W 0 )-lgR i -Q。
further, rootAccording to the sound pressure level LP i Indeed the average sound pressure level LP around the geometric model A The calculation formula is as follows:
wherein N represents a distance R from the surface of the geometric model i Is a total number of arbitrary points.
According to another aspect of the present invention, there is provided a system for determining vibration noise of an extra-high voltage shunt reactor core, the system comprising:
a first calculation unit for applying a voltage or a current to windings of a 1/2 geometric model of a pre-established shunt reactor, the pre-set permeability μ of an acoustic propagation medium being dependent on the current density J and the resistivity ρ on said windings 0 And relative permeability mu of iron core 1 Calculating the magnetic vector A of the iron core and the periphery of the iron core of the geometric model, and determining the magnetic induction intensity B of the iron core and the periphery of the iron core according to the rotation of the magnetic vector A;
a second calculation unit for determining the magnetic field B of the silicon steel sheet adopted by the iron core under a certain pressure according to the magnetic induction intensity B sat And magnetostriction amount epsilon sat Calculating the pre-strain epsilon of the iron core at any position caused by magnetostriction ms ；
A third calculation unit for determining Young's modulus E of the iron core, poisson's ratio alpha and magnetic permeability v of the iron core according to the magnetic induction intensity B 0 And the pre-strain epsilon of the core ms Calculating electromagnetic acting force F of the iron core;
a fourth calculation unit for calculating a displacement u when the iron core vibrates, based on the electromagnetic acting force F, a mass matrix M and a stiffness matrix K of the geometric model determined in advance;
a fifth calculation unit for a geometric model divided in advance according to a preset characteristic impedance ρc of the acoustic propagation medium based on a displacement u of the core during vibrationAcoustic emissivity k of the j-th plane j Area S of jth face c And the included angle theta between the jth surface and any point i around the geometric model j Determining the total energy W of the radiated sound source at the point i i ；
A sixth calculation unit for calculating total energy W of the radiated sound source according to the point i i Distance R from point i to geometric model surface i A preset sound field degree of freedom coefficient Q and reference energy W 0 Determining sound pressure level LP i ；
A seventh calculation unit for calculating a sound pressure level LP based on the sound pressure level LP i Indeed the average sound pressure level LP around the geometric model A 。
Further, the system further comprises:
the model building unit is used for building a 1/2 geometric model of the extra-high voltage shunt reactor, and the geometric model comprises an iron yoke, an iron core, a winding, a marble Dan Diankuai, a winding supporting insulating plate, a wood cushion block, a shielding layer and a metal pull rod device, wherein the iron core and the iron yoke are partially modeled according to the whole, and the iron yoke provides a magnetic circuit for the reactor to enable magnetic lines of force to be closed; the gap formed by the iron core and the marble cushion blocks forms the inductance of the reactor; the marble cushion block adopts the principle of equal area, and a plurality of marble cushion blocks of each layer of the practical shunt reactor are replaced by an integral ring-shaped marble cushion block for supporting the iron cores and providing air gaps between the iron cores; the winding supporting insulating plate and the wood cushion block are used for supporting the winding; the shielding layer is used for shielding the magnetic field to generate eddy current induction on the shell; the metal pull rod device is used for providing pretightening force for the iron core column consisting of the iron core cake and the marble cushion block; the specific structure of winding turns in the actual shunt reactor is ignored by the winding, modeling is carried out according to the whole size, the geometric model is divided into T surfaces by grid division, and the acoustic radiation coefficient of the jth surface is k j The area of each surface unit is S c ；
The first measuring unit is used for measuring the iron core of the geometric model by adopting a pulse excitation method and determining the Young modulus E of the iron core;
a second measuring unit for measuring the magnetostriction curve of the silicon steel sheet adopted by the iron core of the geometric model under a certain pressure by adopting a magnetization and magnetostriction measuring system, wherein the magnetostriction curve is a magnetic field B sat And magnetostriction amount epsilon sat Is a relationship of (2);
a first parameter unit for determining a mass matrix [ M ] and a stiffness matrix [ K ] of the geometric model;
a second parameter unit for setting magnetic permeability mu of the acoustic propagation medium at the position of the geometric model 0 And the characteristic impedance ρc, the core permeability ν 0 Poisson's ratio alpha and relative permeability mu in an acoustic propagation medium 1 And the sound field degree of freedom coefficient Q and the reference energy W 0 。
Further, the first calculation unit loads voltage or current on the windings of the geometric model, and air permeability mu is based on the current density J and the resistivity rho on the windings 0 And relative permeability mu of iron core 1 Calculating the magnetic vector position A of the iron core and the periphery thereof, and determining the magnetic induction intensity B of the iron core and the periphery thereof according to the rotation of the magnetic vector position A, wherein the calculation formula is as follows:
Further, the second calculating unit determines the magnetic field B of the silicon steel sheet adopted by the iron core under a certain pressure according to the magnetic induction intensity B sat And magnetostriction amount epsilon sat Calculating the pre-strain epsilon of the iron core at any position caused by magnetostriction ms The calculation formula is as follows:
ε ms ＝1.5×ε sat ×(B/B sat ) 2 。
further, the third computation sheetThe Young modulus E of the iron core is determined in advance according to the magnetic induction intensity B, the Poisson ratio alpha of the iron core is set in advance, and the magnetic resistance rate v of the iron core is set in advance 0 And the pre-strain epsilon of the core ms And calculating the electromagnetic acting force F of the iron core, wherein the calculation formula is as follows:
where n represents the vector operator in the normal direction of the magnetic field.
Further, the fourth calculation unit calculates the displacement u of the iron core during vibration according to the electromagnetic acting force F and a predetermined mass matrix [ M ] and stiffness matrix [ K ] of the geometric model, and the calculation formula is as follows:
further, the fifth calculation unit sets the characteristic impedance ρc of the acoustic propagation medium in advance according to the acceleration v of the core during vibration, and divides the acoustic radiation coefficient k of the j-th surface of the geometric model in advance j Area S of jth face c And the included angle theta between the jth surface and any point i around the geometric model j Determining the total energy W of the radiated sound source at the point i i The calculation formula is as follows:
further, the sixth calculation unit calculates the total energy W of the radiated sound source according to the point i i Distance R from point i to geometric model surface i A preset sound field degree of freedom coefficient Q and reference energy W 0 Determining soundPressure stage LP i The calculation formula is as follows:
Lp i ＝10lg(W i /W 0 )-lgR i -Q。
further, the seventh calculation unit calculates the sound pressure level LP based on the sound pressure level LP i Indeed the average sound pressure level LP around the geometric model A The calculation formula is as follows:
wherein N represents a distance R from the surface of the geometric model i Is a total number of arbitrary points.
According to the method and the system for determining the vibration noise of the iron core of the ultra-high voltage shunt reactor, provided by the technical scheme of the invention, a vibration simulation experiment model of the shunt reactor is established, magnetostriction force and electromagnetic acting force of the iron core of the simulation geometric model under the action of a magnetic field are calculated, then an acoustic vibration model is adopted, acceleration on the surface of the iron core is utilized as excitation, sound field distribution around the iron core is calculated, and average sound pressure level around the geometric model is determined. According to the method and the system, the Young modulus of the iron core in the simulation geometric model is not the Young modulus of the silicon steel sheet in the prior art, the whole iron core is taken as a whole, the Young modulus of an iron core cake is measured to simulate the lamination structure of the iron core, and the magnetostriction system and the measurement standard are adopted to measure the magnetostriction quantity of the silicon steel sheet used by the iron core under the action of a certain pressure and a fixed magnetic field according to the characteristic that the magnetostriction characteristic of the iron core is greatly influenced by the magnetic field and the compression stress applied to the iron core, and then the prestrain of the iron core caused by magnetostriction is determined according to the telescoping relation curve of the magnetic field and the magnetostriction quantity, so that the influence of which force of the magnetostriction force and the electromagnetic force of the iron core on the vibration noise of the reactor is maximum is more accurately distinguished, and the theoretical design of the body noise reduction of the reactor is important.
Drawings
Exemplary embodiments of the present invention may be more completely understood in consideration of the following drawings:
fig. 1 is a flowchart of a method of determining vibration noise of an extra-high voltage shunt reactor core according to a preferred embodiment of the present invention;
fig. 2 is a schematic diagram of an equivalent scaling model of an extra-high voltage shunt reactor according to a preferred embodiment of the present invention;
fig. 3 is a schematic structural diagram of a system for determining vibration noise of an extra-high voltage shunt reactor core according to a preferred embodiment of the present invention.
Detailed Description
The exemplary embodiments of the present invention will now be described with reference to the accompanying drawings, however, the present invention may be embodied in many different forms and is not limited to the examples described herein, which are provided to fully and completely disclose the present invention and fully convey the scope of the invention to those skilled in the art. The terminology used in the exemplary embodiments illustrated in the accompanying drawings is not intended to be limiting of the invention. In the drawings, like elements/components are referred to by like reference numerals.
Unless otherwise indicated, terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art. In addition, it will be understood that terms defined in commonly used dictionaries should be interpreted as having a meaning that is consistent with their meaning in the context of the relevant art and will not be interpreted in an idealized or overly formal sense.
Fig. 1 is a flowchart of a method of determining vibration noise of an ultra-high voltage shunt reactor core according to a preferred embodiment of the present invention. As shown in fig. 1, a method 100 for determining vibration noise of an extra-high voltage shunt reactor core according to the preferred embodiment starts in step 101.
In step 101, a 1/2 geometric model of the extra-high voltage shunt reactor is established.
Fig. 2 is a schematic diagram of an equivalent scaling model of an extra-high voltage shunt reactor according to a preferred embodiment of the present invention. As shown in FIG. 2, the 1/2 geometric model of the ultra-high voltage shunt reactor comprises an iron yoke 201, an iron core 202, windings 203, a marble Dan Diankuai, winding support insulating plates 205 and woodA head block 206, a shielding layer 207 and a metal pull rod device 208, wherein the iron core 202 and the iron yoke part 201 are integrally modeled, the iron yoke 201 provides a magnetic circuit for the reactor, and magnetic force lines are closed; the gap formed by the iron core 202 and the marble cushion block 204 forms the inductance of the reactor; the marble cushion block 204 adopts the principle of equal area, replaces a plurality of marble cushion blocks of each layer of the practical shunt reactor with an integral ring-shaped marble cushion block for supporting the iron cores and providing air gaps between the iron cores; the winding support insulating plate 205 and the wood pad 206 are used for supporting the winding; the shielding layer 207 is used to shield the magnetic field from eddy current induction on the housing; the metal pull rod device 208 is used for providing pretightening force for the iron core column consisting of the iron core cake 202 and the marble Dan Diankuai; the specific structure of winding turns in the actual shunt reactor is ignored by the winding, modeling is carried out according to the whole size, the geometric model is divided into T surfaces by grid division, and the acoustic radiation coefficient of the jth surface is k j The area of each surface unit is S c . The distance between the iron cores in the geometric model can be adjusted, the characteristics of the model and the actual ultra-high voltage shunt reactor are closer by adjusting the number and the distance of the iron cores, the error of the model in measuring the noise of the reactor can be further reduced, and the measuring precision is improved.
In step 102, a pulse excitation method is adopted to measure the iron core of the geometric model, and the Young modulus E of the iron core is determined. In the prior art, the Young modulus of a common iron core adopts the Young modulus of a silicon steel sheet, but the iron core is of a laminated structure, and the Young modulus of silicon steel cannot be used when the iron core is considered as a whole, so that the Young modulus of an iron core cake is measured firstly; and then setting the Young's modulus parameter of the iron core as the measured Young's modulus parameter of the iron core, thereby simulating the lamination structure of the iron core.
In step 103, a magnetization and magnetostriction measurement system is used to measure the magnetostriction curve of the silicon steel sheet adopted by the iron core of the geometric model under a certain pressure, wherein the magnetostriction curve is a magnetic field B sat And magnetostriction amount epsilon sat Is a relationship of (2). Due to magnetostriction property of iron core and application thereof The magnitude of the added magnetic field and the influence of the compressive stress are larger, so that the IEC60404-3 measurement standard and the magnetostriction measurement system thereof are adopted to fix the pressure P and the fixed magnetic field B of the silicon steel sheet used for the iron core sat Lower magnetostriction amount epsilon sat Obtaining a magnetic field B sat And magnetostriction amount epsilon sat Is a relationship of (2).
In step 104, a mass matrix [ M ] and a stiffness matrix [ K ] of the geometric model are determined.
In step 105, the magnetic permeability μ of the acoustic propagation medium is set at the location of the geometric model 0 And the characteristic impedance ρc, the core reluctance rate ν 0 Poisson's ratio alpha and relative permeability mu in an acoustic propagation medium 1 And the sound field degree of freedom coefficient Q and the reference energy W 0 。
In step 106, a voltage or current is applied to windings of the 1/2 geometric model of the established shunt reactor, and the magnetic permeability mu of the acoustic propagation medium is preset according to the current density J and the resistivity rho on the windings 0 And relative permeability mu of iron core 1 Calculating the magnetic vector A of the iron core and the periphery of the iron core of the geometric model, and determining the magnetic induction intensity B of the iron core and the periphery of the iron core according to the rotation of the magnetic vector A;
in step 107, according to the magnetic induction intensity B, the magnetic field B of the silicon steel sheet adopted by the iron core under a certain pressure sat And magnetostriction amount epsilon sat Calculating the pre-strain epsilon of the iron core at any position caused by magnetostriction ms ；
In step 108, according to the magnetic induction intensity B, young's modulus E of the iron core, poisson's ratio alpha of the iron core and magnetic resistance rate v of the iron core 0 And the pre-strain epsilon of the core ms Calculating electromagnetic acting force F of the iron core;
in step 109, calculating the displacement u of the iron core during vibration according to the electromagnetic acting force F and the mass matrix [ M ] and the rigidity matrix [ K ] of the geometric model;
in step 110, according to the displacement u of the core during vibration, the characteristic impedance ρc of the acoustic propagation medium is preset, and the j-th surface of the geometric model is divided in advanceAcoustic emissivity k j Area S of jth face c And the included angle theta between the jth surface and any point i around the geometric model j Determining the total energy W of the radiated sound source at the point i i ；
In step 111, the total energy W of the sound source is radiated according to the point i i Distance R from point i to geometric model surface i Coefficient of freedom Q of sound field and reference energy W 0 Determining sound pressure level LP i ；
In step 112, according to the sound pressure level LP i Indeed the average sound pressure level LP around the geometric model A 。
Preferably, a voltage or current is applied to the windings of the geometric model, and the air permeability μ is based on the current density J and the resistivity ρ on the windings 0 And relative permeability mu of iron core 1 Calculating the magnetic vector position A of the iron core and the periphery thereof, and determining the magnetic induction intensity B of the iron core and the periphery thereof according to the rotation of the magnetic vector position A, wherein the calculation formula is as follows:
preferably, the magnetic field B of the silicon steel sheet adopted by the predetermined iron core under a certain pressure is based on the magnetic induction intensity B sat And magnetostriction amount epsilon sat Calculating the pre-strain epsilon of the iron core at any position caused by magnetostriction ms The calculation formula is as follows:
ε ms ＝1.5×ε sat ×(B/B sat ) 2 。
preferably, according to the magnetic induction intensity B, the young's modulus E of the iron core is predetermined, the poisson ratio α of the iron core is set in advance, and the magnetic resistance rate v of the iron core is set in advance 0 And the pre-strain epsilon of the core ms Calculating the iron coreThe electromagnetic acting force F of the device is calculated as follows:
where n represents the vector operator in the normal direction of the magnetic field.
Preferably, the displacement u of the iron core during vibration is calculated according to the electromagnetic acting force F and a predetermined mass matrix [ M ] and stiffness matrix [ K ] of the geometric model, and the calculation formula is as follows:
preferably, the acoustic radiation coefficient k of the j-th surface of the geometric model is divided in advance according to the characteristic impedance ρc of the acoustic propagation medium set in advance and the acceleration v when the iron core vibrates j Area S of jth face c And the included angle theta between the jth surface and any point i around the geometric model j Determining the total energy W of the radiated sound source at the point i i The calculation formula is as follows:
wherein T represents the number of meshes determined when the geometric model is subjected to meshing.
Preferably, according to the total energy W of the radiated sound source at the point i i Distance R from point i to geometric model surface i A preset sound field degree of freedom coefficient Q and reference energy W 0 Determining sound pressure level LP i The calculation formula is as follows:
Lp i ＝10lg(W i /W 0 )-lgR i -Q。
preferably, according to said sound pressure level LP i Indeed the average sound pressure level LP around the geometric model A The calculation formula is as follows:
wherein N represents a distance R from the surface of the geometric model i Is a total number of arbitrary points.
In the preferred embodiment, a 1/2 equivalent scaling model of an extra-high voltage shunt reactor is subjected to sound field calculation according to the method of the application. The minimum sound pressure levels under the action of 6.25MPa at the same positions of the two sides of the reactor from the 1.5m plane of the reactor body are 82.79dB respectively, the actual measurement result is 81.2dB, and the error is 1.59dB. The calculation result is relatively close to the actual measurement value.
Fig. 3 is a schematic structural diagram of a system for determining vibration noise of an extra-high voltage shunt reactor core according to a preferred embodiment of the present invention. As shown in fig. 3, a system 300 for determining vibration noise of an extra-high voltage shunt reactor core according to the present preferred embodiment includes:
The model building unit 301 is configured to build a 1/2 geometric model of the extra-high voltage shunt reactor, where the geometric model includes an iron yoke, an iron core, a winding, a marble Dan Diankuai, a winding supporting insulation board, a wood cushion block, a shielding layer and a metal pull rod device, where the iron core and the iron yoke are modeled integrally, and the iron yoke provides a magnetic circuit for the reactor to close magnetic lines; the gap formed by the iron core and the marble cushion blocks forms the inductance of the reactor; the marble cushion block adopts the principle of equal area, and a plurality of marble cushion blocks of each layer of the practical shunt reactor are replaced by an integral ring-shaped marble cushion block for supporting the iron cores and providing air gaps between the iron cores; the winding supporting insulating plate and the wood cushion block are used for supporting the winding; the shielding layer is used for shielding the magnetic field to generate eddy current induction on the shell; the metal pull rod device is used for providing pretightening force for the iron core column consisting of the iron core cake and the marble cushion block; the specific structure that the winding ignores the winding turns in the actual shunt reactor is as a wholeModeling the size, and meshing the geometric model into T surfaces, wherein the acoustic radiation coefficient of the j-th surface is k j The area of each surface unit is S c ；
A first measuring unit 302, configured to measure an iron core of the geometric model by using a pulse excitation method, and determine a young modulus E of the iron core;
a second measuring unit 303 for measuring the magnetostriction curve of the silicon steel sheet used for the iron core of the geometric model under a certain pressure by using a magnetization and magnetostriction measuring system, wherein the magnetostriction curve is a magnetic field B sat And magnetostriction amount epsilon sat Is a relationship of (2);
a first parameter unit 304 for determining a mass matrix [ M ] and a stiffness matrix [ K ] of the geometric model;
a second parameter unit 305 for setting the magnetic permeability μ of the acoustic propagation medium at the location of the geometric model 0 And the characteristic impedance ρc, the core permeability ν 0 Poisson's ratio alpha and relative permeability mu in an acoustic propagation medium 1 And the sound field degree of freedom coefficient Q and the reference energy W 0 。
A first calculation unit 306 for loading a voltage or current on windings of a pre-established 1/2 geometric model of a shunt reactor, a pre-set magnetic permeability μ of the acoustic propagation medium according to the current density J and the resistivity ρ on said windings 0 And relative permeability mu of iron core 1 Calculating the magnetic vector A of the iron core and the periphery of the iron core of the geometric model, and determining the magnetic induction intensity B of the iron core and the periphery of the iron core according to the rotation of the magnetic vector A;
A second calculation unit 307 for determining the magnetic field B of the silicon steel sheet used by the iron core under a certain pressure according to the magnetic induction B sat And magnetostriction amount epsilon sat Calculating the pre-strain epsilon of the iron core at any position caused by magnetostriction ms ；
A third calculation unit 308 for presetting iron according to the magnetic induction intensity B, the Young's modulus E of the iron core, and the preset ironPoisson ratio alpha of core and magnetic permeability v of iron core 0 And the pre-strain epsilon of the core ms Calculating electromagnetic acting force F of the iron core;
a fourth calculation unit 309 for calculating a displacement u of the core during vibration based on the electromagnetic force F, a mass matrix M and a stiffness matrix K of the geometric model determined in advance;
a fifth calculation unit 310 for pre-dividing the acoustic radiation coefficient k of the j-th surface of the geometric model according to the displacement u of the iron core during vibration and the preset characteristic impedance ρc of the acoustic propagation medium j Area S of jth face c And the included angle theta between the jth surface and any point i around the geometric model j Determining the total energy W of the radiated sound source at the point i i ；
A sixth calculation unit 311 for calculating total energy W of the sound source according to the radiation at the i point i Distance R from point i to geometric model surface i A preset sound field degree of freedom coefficient Q and reference energy W 0 Determining sound pressure level LP i ；
A seventh calculation unit 312 for calculating a sound pressure level LP according to the sound pressure level i Indeed the average sound pressure level LP around the geometric model A 。
Preferably, the first computing unit 306 applies a voltage or current to the windings of the geometric model, air permeability μ according to the current density J and resistivity ρ across the windings 0 And relative permeability mu of iron core 1 Calculating the magnetic vector position A of the iron core and the periphery thereof, and determining the magnetic induction intensity B of the iron core and the periphery thereof according to the rotation of the magnetic vector position A, wherein the calculation formula is as follows:
preferably, the secondThe calculating unit 307 determines the magnetic field B of the silicon steel sheet under a certain pressure according to the magnetic induction intensity B sat And magnetostriction amount epsilon sat Calculating the pre-strain epsilon of the iron core at any position caused by magnetostriction ms The calculation formula is as follows:
ε ms ＝1.5×ε sat ×(B/B sat ) 2 。
preferably, the third calculating unit 308 determines the young's modulus E of the iron core, the poisson ratio α of the iron core, and the magnetic reluctance rate v of the iron core according to the magnetic induction intensity B 0 And the pre-strain epsilon of the core ms And calculating the electromagnetic acting force F of the iron core, wherein the calculation formula is as follows:
where n represents the vector operator in the normal direction of the magnetic field.
Preferably, the fourth calculation unit 309 calculates the displacement u of the core during vibration according to the electromagnetic force F, a mass matrix M and a stiffness matrix K of the geometric model, which are predetermined, and the calculation formula is:
preferably, the fifth calculating unit 310 presets the characteristic impedance ρc of the acoustic propagation medium according to the acceleration v of the core during vibration, and presets the acoustic radiation coefficient k of the j-th surface of the geometric model j Area S of jth face c And the included angle theta between the jth surface and any point i around the geometric model j Determining the total energy W of the radiated sound source at the point i i The calculation formula is as follows:
preferably, the sixth calculation unit 311 calculates the total energy W of the radiated sound source according to the point i i Distance R from point i to geometric model surface i A preset sound field degree of freedom coefficient Q and reference energy W 0 Determining sound pressure level LP i The calculation formula is as follows:
Lp i ＝10lg(W i /W 0 )-lgR i -Q。
preferably, the seventh calculation unit calculates the sound pressure level LP based on the sound pressure level i Indeed the average sound pressure level LP around the geometric model A The calculation formula is as follows:
wherein N represents a distance R from the surface of the geometric model i Is a total number of arbitrary points.
The step of determining the sound pressure level around the geometric model by calculating the electromagnetic acting force and the magnetostriction force of the iron core of the simulation geometric model is the same as the step adopted by the method for determining the vibration noise of the iron core of the ultra-high voltage shunt reactor, and the technical effects achieved by the method are the same, and are not repeated here.
The invention has been described with reference to a few embodiments. However, as is well known to those skilled in the art, other embodiments than the above disclosed invention are equally possible within the scope of the invention, as defined by the appended patent claims.
Generally, all terms used in the claims are to be interpreted according to their ordinary meaning in the technical field, unless explicitly defined otherwise therein. All references to "a/an/the [ means, component, etc. ]" are to be interpreted openly as referring to at least one instance of said means, component, etc., unless explicitly stated otherwise. The steps of any method disclosed herein do not have to be performed in the exact order disclosed, unless explicitly stated.
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 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.
Finally, it should be noted that: the above embodiments are only for illustrating the technical aspects of the present invention and not for limiting the same, and although the present invention has been described in detail with reference to the above embodiments, it should be understood by those of ordinary skill in the art that: modifications and equivalents may be made to the specific embodiments of the invention without departing from the spirit and scope of the invention, which is intended to be covered by the claims.
Claims (16)
1. A method for determining vibration noise of an extra-high voltage shunt reactor core, the method comprising:
establishing a 1/2 geometric model of the extra-high voltage shunt reactor, wherein the geometric model comprises an iron yoke, an iron core, a winding, a marble Dan Diankuai, a winding supporting insulating plate, a wood cushion block, a shielding layer and a metal pull rod device, wherein the iron core and the iron yoke are modeled according to the whole, and the iron yoke provides a magnetic circuit for the reactor to close magnetic force lines; the gap formed by the iron core and the marble cushion blocks forms the inductance of the reactor; the marble cushion block adopts the principle of equal area, and a plurality of marble cushion blocks of each layer of the practical shunt reactor are replaced by an integral ring-shaped marble cushion block for supporting the iron cores and providing air gaps between the iron cores; the winding supporting insulating plate and the wood cushion block are used for supporting the winding; the shielding layer is used for shielding the magnetic field to generate eddy current induction on the shell; the metal pull rod device is used for providing pretightening force for the iron core column consisting of the iron core cake and the marble cushion block; the specific structure of winding turns in the actual shunt reactor is ignored by the winding, modeling is carried out according to the whole size, the geometric model is divided into T surfaces by grid division, and the acoustic radiation coefficient of the jth surface is k j The area of each surface unit is S c ；
Measuring an iron core of the geometric model by adopting a pulse excitation method, and determining the Young modulus E of the iron core;
a magnetization and magnetostriction measuring system is adopted to measure the magnetostriction curve of a silicon steel sheet adopted by the iron core of the geometric model under a certain pressure, wherein the magnetostriction curve is a magnetic field B sat And magnetostriction amount epsilon sat Is a relationship of (2);
determining a mass matrix [ M ] and a stiffness matrix [ K ] of the geometric model;
setting magnetic permeability mu of sound propagation medium at position of geometric model 0 And characteristic impedance ρc, core reluctance ratio v 0 Poisson's ratio alpha and relative permeability mu in an acoustic propagation medium 1 And the sound field degree of freedom coefficient Q and the reference energy W 0 ；
Loading voltage or current on windings of a 1/2 geometric model of the established shunt reactor, and presetting the magnetic permeability mu of an acoustic propagation medium according to the current density J and the resistivity rho on the windings 0 And relative permeability mu of iron core 1 Calculating the magnetic vector A of the iron core and the periphery of the iron core of the geometric model, and determining the magnetic induction intensity B of the iron core and the periphery of the iron core according to the rotation of the magnetic vector A;
according to the magnetic induction intensity B, a magnetic field B of a silicon steel sheet adopted by a predetermined iron core under a certain pressure sat And magnetostriction amount epsilon sat Calculating the pre-strain epsilon of the iron core at any position caused by magnetostriction ms ；
According to the magnetic induction intensity B, the Young modulus E of the iron core is determined in advance, the Poisson ratio alpha of the iron core is set in advance, and the magnetic resistance ratio v of the iron core is set in advance 0 And the pre-strain epsilon of the core ms Calculating electromagnetic acting force F of the iron core;
calculating displacement u of the iron core during vibration according to the electromagnetic acting force F and a mass matrix [ M ] and a rigidity matrix [ K ] of the geometric model, which are determined in advance;
according to the displacement u of the iron core during vibration, the characteristic impedance ρc of the preset sound propagation medium is set, and the sound radiation coefficient k of the j-th surface of the geometric model is divided in advance j Area S of jth face c And the j-th faceAn included angle theta between the geometric model and any point i point around the geometric model j Determining the total energy W of the radiated sound source at the point i i ；
According to the total energy W of the radiated sound source at the point i i Distance R from point i to geometric model surface i A preset sound field degree of freedom coefficient Q and reference energy W 0 Determining sound pressure level LP i ；
According to the sound pressure level LP i Indeed the average sound pressure level LP around the geometric model A 。
2. Method according to claim 1, characterized in that a voltage or current is applied to the windings of the geometric model, the air permeability μ being dependent on the current density J and the resistivity ρ on the windings 0 And relative permeability mu of iron core 1 Calculating the magnetic vector position A of the iron core and the periphery thereof, and determining the magnetic induction intensity B of the iron core and the periphery thereof according to the rotation of the magnetic vector position A, wherein the calculation formula is as follows:
3. the method according to claim 1, wherein the magnetic field B of the silicon steel sheet used for the predetermined iron core under a certain pressure is based on the magnetic induction B sat And magnetostriction amount epsilon sat Calculating the pre-strain epsilon of the iron core at any position caused by magnetostriction ms The calculation formula is as follows:
ε ms ＝1.5×ε sat ×(B/B sat ) 2 。
4. the method according to claim 1, characterized in thatCharacterized in that according to the magnetic induction intensity B, the Young modulus E of the iron core is determined in advance, the Poisson ratio alpha of the iron core is set in advance, and the magnetic resistance ratio v of the iron core is set in advance 0 And the pre-strain epsilon of the core ms And calculating the electromagnetic acting force F of the iron core, wherein the calculation formula is as follows:
where n represents the vector operator in the normal direction of the magnetic field.
5. The method according to claim 1, characterized in that the displacement u of the core during vibration is calculated from the electromagnetic force F, a mass matrix [ M ] and a stiffness matrix [ K ] of the geometric model, which are predetermined, with the formula:
6. The method according to claim 1, wherein the acoustic radiation coefficient k of the j-th face of the geometric model divided in advance is based on the characteristic impedance ρc of the acoustic propagation medium set in advance in accordance with the acceleration v at which the core vibrates j Area S of jth face c And the included angle theta between the jth surface and any point i around the geometric model j Determining the total energy W of the radiated sound source at the point i i The calculation formula is as follows:
wherein T represents the number of meshes determined when the geometric model is subjected to meshing.
7. The method according to claim 1, wherein the total energy W of the radiated sound source at the point i i Distance R from point i to geometric model surface i A preset sound field degree of freedom coefficient Q and reference energy W 0 Determining sound pressure level LP i The calculation formula is as follows:
Lp i ＝10lg(W i /W 0 )-lgR i -Q。
8. the method according to claim 1, characterized in that, according to the sound pressure level LP i Indeed the average sound pressure level LP around the geometric model A The calculation formula is as follows:
wherein N represents a distance R from the surface of the geometric model i Is a total number of arbitrary points.
9. A system for determining vibration noise of an extra-high voltage shunt reactor core, the system comprising:
The model building unit is used for building a 1/2 geometric model of the extra-high voltage shunt reactor, and the geometric model comprises an iron yoke, an iron core, a winding, a marble Dan Diankuai, a winding supporting insulating plate, a wood cushion block, a shielding layer and a metal pull rod device, wherein the iron core and the iron yoke are partially modeled according to the whole, and the iron yoke provides a magnetic circuit for the reactor to enable magnetic lines of force to be closed; the gap formed by the iron core and the marble cushion blocks forms the inductance of the reactor; the marble cushion block adopts the principle of equal area, and a plurality of marble cushion blocks of each layer of the practical shunt reactor are replaced by an integral ring-shaped marble cushion block for supporting the iron cores and providing air gaps between the iron cores; the winding supporting insulating plate and the wood cushion block are used for supporting the winding; the shielding layer is used for shielding magnetismThe field creates eddy current induction on the housing; the metal pull rod device is used for providing pretightening force for the iron core column consisting of the iron core cake and the marble cushion block; the specific structure of winding turns in the actual shunt reactor is ignored by the winding, modeling is carried out according to the whole size, the geometric model is divided into T surfaces by grid division, and the acoustic radiation coefficient of the jth surface is k j The area of each surface unit is S c ；
The first measuring unit is used for measuring the iron core of the geometric model by adopting a pulse excitation method and determining the Young modulus E of the iron core;
a second measuring unit for measuring the magnetostriction curve of the silicon steel sheet adopted by the iron core of the geometric model under a certain pressure by adopting a magnetization and magnetostriction measuring system, wherein the magnetostriction curve is a magnetic field B sat And magnetostriction amount epsilon sat Is a relationship of (2);
a first parameter unit for determining a mass matrix [ M ] and a stiffness matrix [ K ] of the geometric model;
a second parameter unit for setting magnetic permeability mu of the acoustic propagation medium at the position of the geometric model 0 And characteristic impedance ρc, core permeability v 0 Poisson's ratio alpha and relative permeability mu in an acoustic propagation medium 1 And the sound field degree of freedom coefficient Q and the reference energy W 0 ；
A first calculation unit for loading voltage or current on windings of a 1/2 geometric model of the established shunt reactor, and presetting magnetic permeability mu of an acoustic propagation medium according to current density J and resistivity rho on the windings 0 And relative permeability mu of iron core 1 Calculating the magnetic vector A of the iron core and the periphery of the iron core of the geometric model, and determining the magnetic induction intensity B of the iron core and the periphery of the iron core according to the rotation of the magnetic vector A;
A second calculation unit for determining the magnetic field B of the silicon steel sheet adopted by the iron core under a certain pressure according to the magnetic induction intensity B sat And magnetostriction amount epsilon sat Calculating the pre-strain epsilon of the iron core at any position caused by magnetostriction ms ；
A third calculation unit for determining Young's modulus E of the iron core, poisson's ratio alpha of the iron core and magnetic permeability v of the iron core according to the magnetic induction intensity B 0 And the pre-strain epsilon of the core ms Calculating electromagnetic acting force F of the iron core;
a fourth calculation unit for calculating a displacement u when the iron core vibrates, based on the electromagnetic acting force F, a mass matrix M and a stiffness matrix K of the geometric model determined in advance;
a fifth calculation unit for dividing the acoustic radiation coefficient k of the j-th surface of the geometric model in advance according to the displacement u of the iron core during vibration and the characteristic impedance ρc of the acoustic propagation medium set in advance j Area S of jth face c And the included angle theta between the jth surface and any point i around the geometric model j Determining the total energy W of the radiated sound source at the point i i ；
A sixth calculation unit for calculating total energy W of the radiated sound source according to the point i i Distance R from point i to geometric model surface i A preset sound field degree of freedom coefficient Q and reference energy W 0 Determining sound pressure level LP i ；
A seventh calculation unit for calculating a sound pressure level LP based on the sound pressure level LP i Indeed the average sound pressure level LP around the geometric model A 。
10. The system according to claim 9, wherein the first calculation unit loads a voltage or current on windings of the geometric model, air permeability μ according to current density J and resistivity ρ on the windings 0 And relative permeability mu of iron core 1 Calculating the magnetic vector position A of the iron core and the periphery thereof, and determining the magnetic induction intensity B of the iron core and the periphery thereof according to the rotation of the magnetic vector position A, wherein the calculation formula is as follows:
11. the system according to claim 9, wherein the second calculation unit calculates a magnetic field B of a silicon steel sheet used for the predetermined iron core under a certain pressure according to the magnetic induction B sat And magnetostriction amount epsilon sat Calculating the pre-strain epsilon of the iron core at any position caused by magnetostriction ms The calculation formula is as follows:
ε ms ＝1.5×ε sat ×(B/B sat ) 2 。
12. the system according to claim 9, wherein the third calculation unit determines, based on the magnetic induction intensity B, a predetermined young's modulus E of the core, a predetermined poisson ratio α of the core, and a core reluctance ratio v 0 And the pre-strain epsilon of the core ms And calculating the electromagnetic acting force F of the iron core, wherein the calculation formula is as follows:
where n represents the vector operator in the normal direction of the magnetic field.
13. The system according to claim 9, wherein the fourth calculation unit calculates the displacement u of the core during vibration from the mass matrix [ M ] and the stiffness matrix [ K ] of the geometric model determined in advance according to the electromagnetic force F, with a calculation formula:
14. the system according to claim 9, wherein the fifth calculation unit sets the characteristic impedance ρc of the acoustic propagation medium in advance based on the acceleration v of the core during vibration, the acoustic radiation coefficient k of the j-th face of the geometric model divided in advance j Area S of jth face c And the included angle theta between the jth surface and any point i around the geometric model j Determining the total energy W of the radiated sound source at the point i i The calculation formula is as follows:
wherein T represents the number of meshes determined when the geometric model is subjected to meshing.
15. The system according to claim 9, wherein the sixth calculation unit calculates the total energy W of the radiated sound source at the point i i Distance R from point i to geometric model surface i A preset sound field degree of freedom coefficient Q and reference energy W 0 Determining sound pressure level LP i The calculation formula is as follows:
Lp i ＝10lg(W i /W 0 )-lgR i -Q。
16. the system of claim 9, wherein the seventh calculation unit is based on the sound pressure level LP i Indeed the average sound pressure level LP around the geometric model A The calculation formula is as follows:
wherein N represents a distance R from the surface of the geometric model i Is a total number of arbitrary points.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011117619.2A CN112464430B (en) | 2020-10-19 | 2020-10-19 | Method and system for determining vibration noise of iron core of extra-high voltage shunt reactor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011117619.2A CN112464430B (en) | 2020-10-19 | 2020-10-19 | Method and system for determining vibration noise of iron core of extra-high voltage shunt reactor |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112464430A CN112464430A (en) | 2021-03-09 |
| CN112464430B true CN112464430B (en) | 2024-03-29 |
Family
ID=74833162
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011117619.2A Active CN112464430B (en) | 2020-10-19 | 2020-10-19 | Method and system for determining vibration noise of iron core of extra-high voltage shunt reactor |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112464430B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113447743B (en) * | 2021-06-28 | 2024-03-12 | 西安西电变压器有限责任公司 | Transformer no-load noise determination method, device and system |
| CN113673188B (en) * | 2021-08-31 | 2024-02-13 | 中国电力科学研究院有限公司 | Method and system for determining air gap structure of shunt reactor |
| CN114216635B (en) * | 2021-11-24 | 2023-08-22 | 广西电网有限责任公司电力科学研究院 | Test device and method for influencing vibration and noise by tensioning screw rod of parallel reactor |
| CN116151169B (en) * | 2022-11-30 | 2024-03-12 | 西安西电变压器有限责任公司 | Method and device for predicting noise intensity of reactor |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2012068072A (en) * | 2010-09-22 | 2012-04-05 | Jfe Steel Corp | Oscillation model determination method for laminate iron core |
| CN105095609A (en) * | 2015-09-21 | 2015-11-25 | 武汉大学 | Transformer electromagnetic vibration noise calculating method based on finite element method |
| CN106570218A (en) * | 2016-10-14 | 2017-04-19 | 中国电力科学研究院 | Method for computing electromagnetic force of transformer iron core based on finite element method |
| CN107195430A (en) * | 2017-05-18 | 2017-09-22 | 华北电力大学 | A kind of Anode saturable reactor of low-vibration noise |
-
2020
- 2020-10-19 CN CN202011117619.2A patent/CN112464430B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2012068072A (en) * | 2010-09-22 | 2012-04-05 | Jfe Steel Corp | Oscillation model determination method for laminate iron core |
| CN105095609A (en) * | 2015-09-21 | 2015-11-25 | 武汉大学 | Transformer electromagnetic vibration noise calculating method based on finite element method |
| CN106570218A (en) * | 2016-10-14 | 2017-04-19 | 中国电力科学研究院 | Method for computing electromagnetic force of transformer iron core based on finite element method |
| CN107195430A (en) * | 2017-05-18 | 2017-09-22 | 华北电力大学 | A kind of Anode saturable reactor of low-vibration noise |
Non-Patent Citations (3)
| Title |
|---|
| Measurement and analysis of sound radiation from coherently vibrating shunt reactors;Yuxing Wang 等;ournal of Mechanical Science and Technology;20191231;第33卷(第1期);全文 * |
| 并联电抗器振动特性及声功率级计算;周兵 等;高电压技术;20191130;第45卷(第11期);全文 * |
| 电抗器铁芯振动噪声的多场耦合分析方法;刘骥 等;电机与控制学报;20160915(第09期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112464430A (en) | 2021-03-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112464430B (en) | Method and system for determining vibration noise of iron core of extra-high voltage shunt reactor | |
| Jingzhu et al. | Electromagnetic vibration noise analysis of transformer windings and core | |
| Li et al. | Core loss calculation based on finite-element method with Jiles–Atherton dynamic hysteresis model | |
| Fleming et al. | Control of resonant acoustic sound fields by electrical shunting of a loudspeaker | |
| CN105095609A (en) | Transformer electromagnetic vibration noise calculating method based on finite element method | |
| CN111324944A (en) | Iron core vibration characteristic calculation method based on COMSOL | |
| JP2013068517A (en) | Noise prediction method for power transformers | |
| Kou et al. | Modeling and analysis of force characteristics for hybrid excitation linear eddy current brake | |
| CN111024343A (en) | Earthquake simulation vibration test method and device for transformer equipment | |
| Li et al. | Structural dynamic vibration absorber using a tuned inerter eddy current damper | |
| Shahaj et al. | A possible method for magnetostrictive reduction of vibration in large electrical machines | |
| CN116136939A (en) | Transformer core magnetic bias vibration calculation method based on ANSYS | |
| Li et al. | Anomalous loss modeling and validation of magnetic materials in electrical engineering | |
| Cornut et al. | From metallurgy to modelling of electrical steels: A multiple approach to their behaviour and use based on physics and experimental investigations | |
| Kim et al. | Load noise prediction of high-voltage transformers by equation Applying 3-D EMCN | |
| Liu et al. | Investigation on the AMDT core model vibration by FEA and validation of the vibration characteristics by the testing platform | |
| Penin et al. | Impact of mechanical deformations of transformer corners on core losses | |
| CN114613453A (en) | Dynamic magnetostriction determination method and system based on JA model | |
| CN216285681U (en) | Magnetic pole and magnetic field generating device | |
| CN114330039A (en) | Finite element simulation method for distribution of vibration noise external field of transformer core | |
| CN115235604A (en) | Method and device for testing sound vibration level of iron core of transformer and electronic equipment | |
| JP2009236904A (en) | Magnetostriction evaluating method of electromagnetic steel sheet for transformer | |
| JP2008082778A (en) | Method and apparatus for measuring natural frequency of iron core of transformer | |
| Haettel et al. | Prediction of transformer core noise | |
| Xian et al. | Numerical estimation and optimization of vibration noise due to magnetostriction and magnetic forces for laminated core structure |
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 |