CN106153176A - A kind of system and method for single-phase oil immersion type transformer noise is estimated - Google Patents
A kind of system and method for single-phase oil immersion type transformer noise is estimated Download PDFInfo
- Publication number
- CN106153176A CN106153176A CN201610460417.5A CN201610460417A CN106153176A CN 106153176 A CN106153176 A CN 106153176A CN 201610460417 A CN201610460417 A CN 201610460417A CN 106153176 A CN106153176 A CN 106153176A
- Authority
- CN
- China
- Prior art keywords
- solid
- matrix
- coupling
- acoustic
- unit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
- G01H9/002—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means for representing acoustic field distribution
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H11/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by detecting changes in electric or magnetic properties
- G01H11/06—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by detecting changes in electric or magnetic properties by electric means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
- G01H9/004—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means using fibre optic sensors
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a kind of noise prediction system, for estimating the noise of phase oil-filled transformer, this system includes: transformator component finite element unit, determines the solid unit stiffness matrix K relevant to solidS, damping matrix CSAnd mass matrix MS.Hydroacoustics finite element unit, according to casing extraneous air parameter and the oil body parameter of inside transformer fluid field, sets up hydroacoustics finite element equation.Solid and fluid coupling member, set up solid and fluid coupling kinetic equation.Test cell, test cell includes fibre optic accelerometer, Acquisition Instrument and sound meter.Electromagnetic excitation simulation input block, applies the excitation vector of the multiple harmonic component of solid structure in electromagnetic excitation simulation input block.Data correction unit, couples solid and fluid kinetic equation as prediction model, utilizes acoustic mode and mathematic interpolation coupling kinetic equation to obtain the acoustic pressure of sensitive spot.
Description
Technical field
The present invention relates to noise prediction field, more particularly, to one for single-phase oil immersion type transformer noise is entered
The system and method that row is estimated.
Background technology
Along with fast development especially people's growth in the living standard in city, per capita household electricity consumption is substantially increased, and uses
Electricity demand is in the most zooming situation, and supply of electric power is the most in short supply, and transformer substation construction quantity demand increases year by year.And city
Densely inhabited district is exactly the area that load is high in fact, and therefore the addressing of most transformer stations must go deep into load center, thus
Cause the environmental sensitive area noise figures such as resident living area, school, hospital to exceed standard, have a strong impact on work and the life of nearby residents
Live.Substation operation noise just becomes the most serious environmental problem to the impact of surrounding enviroment.
Power transformer is the Main Noise Sources of transformer station, from product design, produce dispatch from the factory, bid and purchase, planning and designing
And install, until the links that puts into operation is all considered as the radiated noise limit value of transformator.In addition to existing transformer station noise
Control and administer to need also exist for the noise of transformator being estimated and testing.Transformer noise is not carried out pre-by prior art
The system estimated or method.
Summary of the invention
In order to solve the problems referred to above, the invention provides a kind of system, for the noise of single-phase oil immersion type transformator is entered
Row is estimated, and this system includes:
Transformator component finite element unit, according to iron core parameter, winding parameter and the casing of single-phase oil immersion type transformator
Parameter, determines the solid unit stiffness matrix K relevant to solidS, damping matrix CSAnd mass matrix MS;
Hydroacoustics finite element unit, according to casing extraneous air parameter and the oil body parameter of inside transformer fluid field,
Set up hydroacoustics finite element equation, and obtain the parameters,acoustic relevant to fluid, wherein inside transformer is closed in space
Sound field be separated into more than one little sound field unit;
Solid and fluid coupling member, set up the coupling contact face of solid and fluid, to carry out between solid and fluid
Power is conducted, according to solid unit stiffness matrix KS, damping matrix CSAnd mass matrix MSSet up the solid transformator structure after coupling
The kinetics equation of part finite element and set up the hydroacoustics equation after coupling, by consolidating after described coupling according to parameters,acoustic
Hydroacoustics equation after the kinetics equation of body transformator component finite element and coupling merges computing, thus sets up solid
Kinetic equation is coupled with fluid;
Test cell, described test cell includes fibre optic accelerometer, Acquisition Instrument and sound meter, and described test cell is installed
In described iron core and the outside of winding, described fibre optic accelerometer is for testing the ferrum of described transformator component finite element unit
The vibration parameters of core, winding and tank surface, described sound meter is for measuring the sound pressure level outside described casing;
Electromagnetic excitation simulation input block, applies the multiple harmonic wave of solid structure in described electromagnetic excitation simulation input block
The excitation vector of component;
Data correction unit, couples kinetic equation as prediction model using described solid and fluid, utilize acoustic mode and
Mathematic interpolation coupling kinetic equation is to obtain the acoustic pressure of sensitive spot;Described prediction model is by described in described test optical fiber unit
The data of fibre optic accelerometer and the test of described sound meter carry out parameters revision, to improve the degree of accuracy of described prediction model.
Preferably, described hydroacoustics finite element equation is:
(Kf+jωCf-ω2Mf)·{pi}={ Ffi,
K in formulafFor acoustics stiffness matrix, CfFor acoustic damping matrix, MfFor acoustic mass matrix;piIt it is the i-th little sound field list
The acoustic pressure of unit, FfiFor acoustic excitation vector, j is complex representation, and ω is circular frequency.
Preferably, the sound field closed in space is separated into more than one little sound field by described hydroacoustics FEM (finite element) model
Unit, according to the formula described little sound field element length L of calculating:
In formula, c is that sound teaches speed, f in liquid mediummaxFor computation model pass frequency of heart
Big value.
Preferably, described casing extraneous air parameter longest edge 1/2nd based on described fuel tank outer surface length range
Interior air fluid.
Preferably, the kinetics equation of the described solid transformator component finite element set up after coupling is:
(Ks+jωCs-ω2μMs)·{ui}+Kc{pi}={ Fsi,
In formula, KsFor described winding and the stiffness matrix of described iron core, CsFor interface damping matrix, MsFor described winding
With the combination quality matrix of described iron core, KcFor Coupling stiffness matrix;FsiFor load excitation vector, μ is winding described in transformator
With the mass ratio of described iron core, uiFor construction unit displacement, piBeing the acoustic pressure of the i-th little sound field unit, j is complex representation, and ω is circle
Frequency.
Preferably, the described hydroacoustics equation set up after coupling includes: described winding and described core structure vibration speed
Degree can input as the addition speed of acoustic pressure, and hydroacoustics equation becomes:
(Kf+jωCf-ω2Mf)·{pi}-ω2Mc{ui}={ Ffi,
KfFor acoustics stiffness matrix, CfFor acoustic damping matrix, MfFor acoustic mass matrix;piIt it is the i-th little sound field unit
Acoustic pressure, FfiFor acoustic excitation vector, j is complex representation, and ω is circular frequency, w2For circular frequency square, McFor coupling mass matrix,
uiFor construction unit displacement.
Preferably, described set up solid and fluid coupling kinetic equation be:
P in formula0For the initial acoustic pressure of fluid, Ks is described winding and the stiffness matrix of described iron core, KfFor acoustics rigidity square
Battle array, Kc is Coupling stiffness matrix, and j is complex representation, and w is circular frequency, and Cs is interface damping matrix, CfFor acoustic damping square
Battle array, w2For circular frequency square, u is winding described in transformator and the mass ratio of described iron core, and Ms is described winding and described iron core
Combination quality matrix, McFor coupling mass matrix, MfFor acoustic mass matrix, uiFor construction unit displacement, piIt it is the i-th fluid
The acoustic pressure of micro unit, FsiFor load excitation vector, FfiFor acoustic excitation vector.
Preferably, described test cell, use the sound meter outside A weighted sound pressure level of test transformator simultaneously.
Preferably, in described electromagnetic excitation simulation input model, apply the excitation vector of the multiple harmonic component of solid structure
For:
F in formulasiFor load excitation vector, n is the frequency multiplication number of times chosen, and m is unit
Quality, f50nFor harmonic frequency, t is the time, ∑ A50nFor harmonic component amplitude.
Preferably, for estimating single-phase oil immersion type transformer noise, the method includes:
Set up transformator component finite element unit, according to the iron core parameter of single-phase oil immersion type transformator, winding parameter and
Box body parameter, determines the solid unit stiffness matrix K relevant to solidS, damping matrix CsAnd mass matrix MS;
Set up hydroacoustics finite element unit, join according to the oil body of casing extraneous air parameter and inside transformer fluid field
Number, sets up hydroacoustics finite element equation, and obtains the parameters,acoustic relevant to fluid, wherein inside transformer is closed space
In sound field be separated into more than one little sound field unit;
Set up solid and fluid coupling member, set up the coupling contact face of solid and fluid, with carry out solid and fluid it
Between power conduction, according to solid unit stiffness matrix KS, damping matrix CSAnd mass matrix MSSet up the solid transformation after coupling
The kinetics equation of device component finite element and set up the hydroacoustics equation after coupling, after described coupling according to parameters,acoustic
The kinetics equation of solid transformator component finite element and coupling after hydroacoustics equation merge computing, thus set up
Solid and fluid coupling kinetic equation;
Setting up test cell, described test cell includes fibre optic accelerometer, Acquisition Instrument and sound meter, described test cell
Being arranged on the outside of described iron core and winding, described fibre optic accelerometer is for testing described transformator component finite element unit
The vibration parameters of iron core, winding and tank surface, described sound meter is for measuring the sound pressure level outside described casing;
Set up electromagnetic excitation simulation input block, in described electromagnetic excitation simulation input block, apply solid unit multiple
The excitation vector of harmonic component;
Set up data correction unit, described solid and fluid are coupled kinetic equation as prediction model, utilize acoustic mode
State and mathematic interpolation coupling kinetic equation are to obtain the acoustic pressure of sensitive spot;Described prediction model is by described in described test cell
The data of fibre optic accelerometer and the test of described sound meter carry out parameters revision, to improve the degree of accuracy of described prediction model.
The present invention is a kind of single-phase oil immersion type transformer noise Prediction System and method.Use the present invention can set up single-phase
Oil-filled transformer noise prediction system, analyzes the distribution obtaining transformator inside and outside radiated sound field, accurately estimates transformator
Radiated noise, thus apply for newly-built transformer station Noise Design and noise reduction technology and estimation results accurately is provided.
Accompanying drawing explanation
By reference to the following drawings, the illustrative embodiments of the present invention can be more fully understood by:
Fig. 1 is the structural representation of the noise prediction system according to embodiment of the present invention;And
Fig. 2 is the schematic flow sheet of the noise prediction method according to embodiment of the present invention.
Detailed description of the invention
With reference now to accompanying drawing, introduce the illustrative embodiments of the present invention, but, the present invention can be by many different shapes
Formula is implemented, and is not limited to embodiment described herein, it is provided that these embodiments are at large and fully disclose
The present invention, and the scope of the present invention is fully passed on to person of ordinary skill in the field.For representing showing in the accompanying drawings
Term in example embodiment is not limitation of the invention.In the accompanying drawings, identical cells/elements uses identical attached
Figure labelling.
Except as otherwise noted, person of ordinary skill in the field is had by term used herein (including scientific and technical terminology)
Common understand implication.Further it will be understood that the term limited with normally used dictionary, it should it is understood to and it
The linguistic context of association area has consistent implication, and is not construed as Utopian or the most formal meaning.
Fig. 1 is the structural representation of the noise prediction system according to embodiment of the present invention.The present invention is a kind of single phase
Immersion transformer noise Prediction System, is used for estimating single-phase oil immersion type transformer noise, analyzes and obtains the radiation of transformator inside and outside
The distribution of sound field, accurately estimates the radiated noise of transformator.Oil-filled transformer is mainly by iron core, folder, winding, oil body, case
The multiple appurtenances of body and inside and outside form.Transformer noise mainly magnetostriction and winding electric power by iron core causes, logical
The coupling crossing transformator oil body propagates to casing, and then produces the radiated noise of external box of microvave.The present invention is based on solid and fluid
Coupling communication theory, establishes sound field coupled wave equation between solid and fluid, uses the method for acoustics finite element to obtain transformation
The sound pressure level of device inside and outside each unit node, the interpolated precise irradiation noise being calculated all sensitive spots, calculate required ginseng
Number can be beneficial to optical fiber sensing technology and install optical fiber acceleration transducer on Transformer Winding, iron core and oil tank wall and measure, and becomes
Depressor radiated noise uses sound meter to record, and obtains single-phase oil immersion type transformer noise accurately eventually through Modifying model and estimates
System, thus apply for newly-built transformer station Noise Design and noise reduction technology and estimation results accurately is provided.
As it is shown in figure 1, the system for estimating the noise of single-phase oil immersion type transformator includes: transformator component has
Limit unit's unit 101,2 hydroacoustics finite element unit 102,3 solid and fluid coupling member 103,4 test cell 104,5 electromagnetism
Excitation simulation input block 105 and 6 data correction unit 106.The system of the present invention can obtain oil immersed type transformation accurately
Device noise prediction value.
Preferably, transformator component finite element unit 101, according to iron core parameter, the winding ginseng of single-phase oil immersion type transformator
Number and box body parameter, determine the solid unit stiffness matrix K relevant to solidS, damping matrix CSAnd mass matrix MS.Single
Phase oil-filled transformer is mainly made up of the multiple Accessory Members in iron core, folder, oil body, casing and inside and outside.Setting up transformator
During component finite element unit, the appurtenances such as negligible internal cabling, switch, small casing, simplify transformator member unit, Gu
Body unit is defined as iron core, winding and three components of casing, is divided into finite element unit, determines solid unit stiffness matrix KS、
Damping matrix CSAnd mass matrix MS。
Preferably, hydroacoustics finite element unit 102, according to casing extraneous air parameter and inside transformer fluid field
Oil body parameter, sets up hydroacoustics finite element equation, and obtains the parameters,acoustic relevant to fluid, wherein inside transformer sealed
Close the fluid sound field in space and be separated into a number of little sound field unit.The present invention fluctuates according to according to acoustics Helmholtz
Equation, is separated into a number of little sound field unit by the sound field closed in space.By shape function and boundary condition difference meter
Calculate the acoustic pressure distribution obtaining all site.Hydroacoustics finite element equation can be converted into:
(Kf+jωCf-ω2Mf)·{pi}={ Ffi,
K in formulafFor acoustics stiffness matrix, CfFor acoustic damping matrix, MfFor acoustic mass matrix;piIt it is the i-th little sound field list
The acoustic pressure of unit, FfiFor acoustic excitation vector;J is complex representation, and ω is circular frequency.The present invention uses acoustic mode method to solve stream
Body acoustics finite element equation.Preferably, casing extraneous air parameter longest edge 1/2nd based on described fuel tank outer surface is long
Air fluid in the range of degree.Preferably, the sound field closed in space is separated into more than one by hydroacoustics FEM (finite element) model
Little sound field unit, according to the formula described little sound field element length L of calculating:
In formula, c is that sound teaches speed, f in liquid mediummaxPass frequency of heart for computation model is maximum
Value, the present invention is during calculating, and the audio parameter choosing input has various frequency, and wide range, and the present invention is at meter
During calculation, choosing frequency range of interest is 50Hz to 2000Hz, it is preferable that be concerned about in frequency selection purposes range set close
1000Hz。
Preferably, solid and fluid coupling member 103, set up the coupling contact face of solid and fluid, with carry out solid and
Power conduction between fluid.At coupling boundary, the vibration velocity of surface of solids normal direction is equal with the vibration velocity of fluid,
Sound acts on the acoustic pressure load on the solid structures such as inner core winding can be as additional normal plane load, after setting up coupling
The kinetics equation of solid transformator component finite element be:
(Ks+jωCs-ω2μMs)·{ui}+Kc{pi}={ Fsi,
In formula, KsFor described winding and the stiffness matrix of described iron core, CsFor interface damping matrix, MsFor winding and ferrum
The combination quality matrix of core, KcFor Coupling stiffness matrix;FsiFor load excitation vector, μ is the quality of Transformer Winding and iron core
Ratio, uiFor construction unit displacement, piBeing the acoustic pressure of the i-th little sound field unit, j is complex representation, and ω is circular frequency.
Preferably, the vibration velocity of winding and iron core inputs as the addition speed of acoustic pressure, and hydroacoustics equation becomes:
(Kf+jωCf-ω2Mf)·{pi}-ω2Mc{ui}={ Ffi,
KfFor acoustics stiffness matrix, CfFor acoustic damping matrix, MfFor acoustic mass matrix;piIt it is the i-th little sound field unit
Acoustic pressure, FfiFor acoustic excitation vector, j is complex representation, and ω is circular frequency, w2For circular frequency square, McFor coupling mass matrix,
uiFor construction unit displacement.
Hydroacoustics equation after the kinetics equation of the solid transformator component finite element after coupling and coupling is carried out
Union operation, thus set up solid and fluid coupling kinetic equation be:
P in formula0For the initial acoustic pressure of fluid, Ks is described winding and the stiffness matrix of described iron core, KfFor acoustics rigidity square
Battle array, Kc is Coupling stiffness matrix, and j is complex representation, and ω is circular frequency, and Cs is interface damping matrix, CfFor acoustic damping square
Battle array, w2For circular frequency square, u is winding described in transformator and the mass ratio of described iron core, and Ms is described winding and described iron core
Combination quality matrix, McFor coupling mass matrix, MfFor acoustic mass matrix, uiFor construction unit displacement, piIt it is the i-th fluid
The acoustic pressure of micro unit, FsiFor load excitation vector, FfiFor acoustic excitation vector.
Preferably, test cell 104, test cell includes fibre optic accelerometer, collecting fiber instrument and sound meter, and test is single
Unit is arranged on the outside of iron core and winding, fibre optic accelerometer and collecting fiber instrument and is used in combination, and is used for testing transformator structure
The vibration parameters of the iron core of part finite element unit, winding and tank surface, sound meter is for measuring the sound outside described casing
Arbitrarily downgrading, in the present invention, sound pressure level is mainly used in the environment noise of measuring transformer near zone.The survey that the embodiment of the present invention is chosen
Test ring border is quiet environment, and sound pressure meter input parameter is 0.
Preferably, electromagnetic excitation simulation input block 105, in electromagnetic excitation simulation input block, apply solid structure many
Plant the excitation vector of harmonic component.The winding that transformer noise is mainly caused by magnetostriction and the load current of iron core silicon-steel sheet
Vibration is produced, and the present invention, in order to simplify electromagnetic exciting, applies solid in the noise prediction system of single-phase oil immersion type transformator
The excitation vector of the multiple harmonic component of structure is:
F in formulasiFor load excitation vector, n is the frequency multiplication number of times chosen, and m is unit
Quality, f50nFor harmonic frequency, t is the time, Σ A50nFor harmonic component amplitude.
The test environment that the embodiment of the present invention is chosen is quiet environment, and sound pressure meter input parameter is 0.Acoustic excitation vector
FfiIt is 0.
Preferably, data correction unit 106, solid and fluid are coupled kinetic equation as prediction model, utilize acoustics
Mode and mathematic interpolation coupling kinetic equation are to obtain the acoustic pressure of sensitive spot.Prediction model is by light described in test optical fiber unit
The data of fine accelerometer and the test of described sound meter carry out parameters revision, to improve the degree of accuracy of described prediction model.
Fig. 2 is according to another embodiment of the present invention, and the present invention provides a kind of predictor method, for becoming single-phase oil immersion type
Depressor noise is estimated, and the method includes:
In step 201, set up transformator component finite element unit, according to the iron core parameter of single-phase oil immersion type transformator, around
Group parameter and box body parameter, determine the solid unit stiffness matrix K relevant to solidS, damping matrix CsAnd mass matrix
MS;
In step 202, set up hydroacoustics finite element unit, according to casing extraneous air parameter and inside transformer fluid
The oil body parameter of field, sets up hydroacoustics finite element equation, and obtains the parameters,acoustic relevant to fluid, wherein by transformator
Sound field in closing space, portion is separated into more than one little sound field unit;
In step 203, set up solid and fluid coupling member, set up the coupling contact face of solid and fluid, solid to carry out
Power conduction between body and fluid, according to solid unit stiffness matrix KS, damping matrix CSAnd mass matrix MSAfter setting up coupling
Solid transformator component finite element kinetics equation and according to parameters,acoustic set up coupling after hydroacoustics equation, will
Hydroacoustics equation after the kinetics equation of the solid transformator component finite element after described coupling and coupling merges fortune
Calculate, thus set up solid and fluid coupling kinetic equation;
In step 204, setting up test cell, described test cell includes fibre optic accelerometer, Acquisition Instrument and sound meter, institute
Stating test cell and be arranged on the outside of described iron core and winding, described fibre optic accelerometer is used for testing described transformator component to be had
The limit unit iron core of unit, winding and the vibration parameters of tank surface, described sound meter is for measuring the sound outside described casing
Arbitrarily downgrade;
In step 205, set up electromagnetic excitation simulation input block, apply solid in described electromagnetic excitation simulation input block
The excitation vector of the multiple harmonic component of body unit;
In step 206, set up data correction unit, described solid and fluid coupled kinetic equation as prediction model,
Utilize acoustic mode and mathematic interpolation coupling kinetic equation to obtain the acoustic pressure of sensitive spot;Described prediction model is by described test
The data of fibre optic accelerometer described in unit and the test of described sound meter carry out parameters revision, to improve described prediction model
Degree of accuracy.
Preferably, described hydroacoustics finite element equation is:
(Kf+jωCf-ω2Mf)·{pi}={ Ffi,
K in formulafFor acoustics stiffness matrix, CfFor acoustic damping matrix, MfFor acoustic mass matrix;piIt it is the i-th little sound field list
The acoustic pressure of unit, FfiFor acoustic excitation vector, j is complex representation, and ω is circular frequency.
Preferably, the sound field closed in space is separated into more than one little sound field by described hydroacoustics FEM (finite element) model
Unit, according to the formula described little sound field element length L of calculating:
In formula, c is that sound teaches speed, f in liquid mediummaxFor computation model pass frequency of heart
Big value.
Preferably, described casing extraneous air parameter longest edge 1/2nd based on described fuel tank outer surface length range
Interior air fluid.
Preferably, the kinetics equation of the described solid transformator component finite element set up after coupling is:
(Ks+jωCs-ω2μMs)·{ui}+Kc{pi}={ Fsi,
In formula, KsFor described winding and the stiffness matrix of described iron core, CsFor interface damping matrix, MsFor described winding
With the combination quality matrix of described iron core, KcFor Coupling stiffness matrix;FsiFor load excitation vector, μ is winding described in transformator
With the mass ratio of described iron core, uiFor construction unit displacement, piBeing the acoustic pressure of the i-th little sound field unit, j is complex representation, and ω is circle
Frequency.
Preferably, the described hydroacoustics equation set up after coupling includes: described winding and described core structure vibration speed
Degree can input as the addition speed of acoustic pressure, and hydroacoustics equation becomes:
(Kf+jωCf-ω2Mf)·{pi}-ω2Mc{ui}={ Ffi,
KfFor acoustics stiffness matrix, CfFor acoustic damping matrix, MfFor acoustic mass matrix;piIt it is the i-th little sound field unit
Acoustic pressure, FfiFor acoustic excitation vector, j is complex representation, and ω is circular frequency, w2For circular frequency square, McFor coupling mass matrix,
uiFor construction unit displacement.
Preferably, described set up solid and fluid coupling kinetic equation be:
P in formula0For the initial acoustic pressure of fluid, Ks is described winding and the stiffness matrix of described iron core, KfFor acoustics rigidity square
Battle array, Kc is Coupling stiffness matrix, and j is complex representation, and w is circular frequency, and Cs is interface damping matrix, CfFor acoustic damping square
Battle array, w2For circular frequency square, u is winding described in transformator and the mass ratio of described iron core, and Ms is described winding and described iron core
Combination quality matrix, McFor coupling mass matrix, MfFor acoustic mass matrix, uiFor construction unit displacement, piIt it is the i-th fluid
The acoustic pressure of micro unit, FsiFor load excitation vector, FfiFor acoustic excitation vector.
Preferably, described test cell, use the sound meter outside A weighted sound pressure level of test transformator simultaneously.
Preferably, in described electromagnetic excitation simulation input model, apply the excitation vector of the multiple harmonic component of solid structure
For:
F in formulasiFor load excitation vector, n is the frequency multiplication number of times chosen, and m is unit
Quality, f50nFor harmonic frequency, t is the time, ∑ A50nFor harmonic component amplitude.
By describing the present invention with reference to a small amount of embodiment.But, known in those skilled in the art, as
Subsidiary Patent right requirement is limited, except the present invention other embodiment disclosed above falls equally in the present invention
In the range of.
Normally, all terms used in the claims are all solved in the usual implication of technical field according to them
Release, unless additionally defined clearly wherein.All of reference " one/described/be somebody's turn to do [device, assembly etc.] " is all opened ground
It is construed at least one example in described device, assembly etc., unless otherwise expressly specified.Any method disclosed herein
Step all need not be run, unless explicitly stated otherwise with disclosed order accurately.
Claims (10)
1. a noise prediction system, for estimating the noise of single-phase oil immersion type transformator, this system includes:
Transformator component finite element unit, according to iron core parameter, winding parameter and the box body parameter of single-phase oil immersion type transformator,
Determine the solid unit stiffness matrix K relevant to solidS, damping matrix CSAnd mass matrix MS;
Hydroacoustics finite element unit, according to casing extraneous air parameter and the oil body parameter of inside transformer fluid field, sets up
Hydroacoustics finite element equation, and obtain the parameters,acoustic relevant to fluid, wherein inside transformer is closed the sound in space
Field is separated into more than one little sound field unit;
Solid and fluid coupling member, set up the coupling contact face of solid and fluid, to carry out the power biography between solid and fluid
Lead, according to solid unit stiffness matrix KS, damping matrix CSAnd mass matrix MSSet up the solid transformator component after coupling to have
Limit the kinetics equation of unit and set up the hydroacoustics equation after coupling according to parameters,acoustic, the solid after described coupling is become
Hydroacoustics equation after the kinetics equation of depressor component finite element and coupling merges computing, thus sets up solid and stream
Body coupling kinetic equation;
Test cell, described test cell includes that fibre optic accelerometer, Acquisition Instrument and sound meter, described test cell are arranged on institute
State the outside of iron core and winding, described fibre optic accelerometer for test described transformator component finite element unit iron core, around
Group and the vibration parameters of tank surface, described sound meter is for measuring the sound pressure level outside described casing;
Electromagnetic excitation simulation input block, applies the multiple harmonic component of solid structure in described electromagnetic excitation simulation input block
Excitation vector;
Data correction unit, couples described solid and fluid kinetic equation as prediction model, utilizes acoustic mode and difference
Calculate coupling kinetic equation to obtain the acoustic pressure of sensitive spot;Described prediction model is accelerated by optical fiber described in described test cell
The data of degree meter and the test of described sound meter carry out parameters revision, to improve the degree of accuracy of described prediction model.
System the most according to claim 1, described hydroacoustics finite element equation is:
(Kf+jωCf-ω2Mf)·{pi}={ Ffi,
K in formulafFor acoustics stiffness matrix, CfFor acoustic damping matrix, MfFor acoustic mass matrix;piIt it is the i-th little sound field unit
Acoustic pressure, FfiFor acoustic excitation vector, j is complex representation, and ω is circular frequency.
System the most according to claim 1, the sound field closed in space is separated into by described hydroacoustics FEM (finite element) model
More than one little sound field unit, according to the formula described little sound field element length L of calculating:
In formula, c is that sound teaches speed, f in liquid mediummaxCare frequency maxima for computation model.
System the most according to claim 1, described casing extraneous air parameter longest edge based on described fuel tank outer surface
Air fluid in 1/2nd length ranges.
System the most according to claim 1, the kinetics side of the solid transformator component finite element after described foundation coupling
Cheng Wei:
(Ks+jωCs-ω2μMs)·{ui}+Kc{pi}={ Fsi,
In formula, KsFor described winding and the stiffness matrix of described iron core, CsFor interface damping matrix, MsFor described winding and institute
State the combination quality matrix of iron core, KcFor Coupling stiffness matrix;FsiFor load excitation vector, μ is winding described in transformator and institute
State the mass ratio of iron core, uiFor construction unit displacement, piBeing the acoustic pressure of the i-th little sound field unit, j is complex representation, and ω is circle frequency
Rate.
System the most according to claim 1, the hydroacoustics equation after described foundation coupling includes: described winding and described
Core structure vibration velocity can input as the addition speed of acoustic pressure, and hydroacoustics equation becomes:
(Kf+jωCf-ω2Mf)·{pi}-ω2Mc{ui}={ Ffi,
KfFor acoustics stiffness matrix, CfFor acoustic damping matrix, MfFor acoustic mass matrix;piIt is the acoustic pressure of the i-th little sound field unit,
FfiFor acoustic excitation vector, j is complex representation, and ω is circular frequency, w2For circular frequency square, McFor coupling mass matrix, uiFor
Construction unit displacement.
System the most according to claim 1, described set up solid and fluid coupling kinetic equation be:
P in formula0For the initial acoustic pressure of fluid, Ks is described winding and the stiffness matrix of described iron core, KfFor acoustics stiffness matrix, Kc
For Coupling stiffness matrix, j is complex representation, and w is circular frequency, and Cs is interface damping matrix, CfFor acoustic damping matrix, w2For
Circular frequency square, u is winding described in transformator and the mass ratio of described iron core, and Ms is the combination matter of described winding and described iron core
Moment matrix, McFor coupling mass matrix, MfFor acoustic mass matrix, uiFor construction unit displacement, piIt it is the i-th fluid micro unit
Acoustic pressure, FsiFor load excitation vector, FfiFor acoustic excitation vector.
System the most according to claim 1, described test cell, use the sound meter outside A weighted of test transformator simultaneously
Sound pressure level.
System the most according to claim 1, applies solid structure multiple humorous in described electromagnetic excitation simulation input model
The excitation vector of wave component is:
F in formulasiFor load excitation vector, n is the frequency multiplication number of times chosen, and m is element quality,
f50nFor harmonic frequency, t is the time, ∑ A50nFor harmonic component amplitude.
10. a predictor method, for estimating single-phase oil immersion type transformer noise, the method includes:
Set up transformator component finite element unit, according to iron core parameter, winding parameter and the casing of single-phase oil immersion type transformator
Parameter, determines solid unit stiffness matrix Ks, the damping matrix C relevant to solidsAnd mass matrix MS;
Set up hydroacoustics finite element unit, according to casing extraneous air parameter and the oil body parameter of inside transformer fluid field,
Set up hydroacoustics finite element equation, and obtain the parameters,acoustic relevant to fluid, wherein inside transformer is closed in space
Sound field be separated into more than one little sound field unit;
Set up solid and fluid coupling member, set up the coupling contact face of solid and fluid, to carry out between solid and fluid
Power is conducted, according to solid unit stiffness matrix KS, damping matrix CSAnd mass matrix MSSet up the solid transformator structure after coupling
The kinetics equation of part finite element and set up the hydroacoustics equation after coupling, by consolidating after described coupling according to parameters,acoustic
Hydroacoustics equation after the kinetics equation of body transformator component finite element and coupling merges computing, thus sets up solid
Kinetic equation is coupled with fluid;
Setting up test cell, described test cell includes fibre optic accelerometer, Acquisition Instrument and sound meter, and described test cell is installed
In described iron core and the outside of winding, described fibre optic accelerometer is for testing the ferrum of described transformator component finite element unit
The vibration parameters of core, winding and tank surface, described sound meter is for measuring the sound pressure level outside described casing;
Set up electromagnetic excitation simulation input block, in described electromagnetic excitation simulation input block, apply the multiple harmonic wave of solid unit
The excitation vector of component;
Set up data correction unit, described solid and fluid coupled kinetic equation as prediction model, utilize acoustic mode and
Mathematic interpolation coupling kinetic equation is to obtain the acoustic pressure of sensitive spot;Described prediction model is by optical fiber described in described test cell
The data of accelerometer and the test of described sound meter carry out parameters revision, to improve the degree of accuracy of described prediction model.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610460417.5A CN106153176B (en) | 2016-06-22 | 2016-06-22 | System and method for predicting noise of single-phase oil-immersed transformer |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610460417.5A CN106153176B (en) | 2016-06-22 | 2016-06-22 | System and method for predicting noise of single-phase oil-immersed transformer |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106153176A true CN106153176A (en) | 2016-11-23 |
CN106153176B CN106153176B (en) | 2020-04-24 |
Family
ID=57353714
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610460417.5A Active CN106153176B (en) | 2016-06-22 | 2016-06-22 | System and method for predicting noise of single-phase oil-immersed transformer |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106153176B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106599395A (en) * | 2016-11-25 | 2017-04-26 | 国网上海市电力公司 | Numerical simulation calculation method for noise of oil immersed transformer |
CN108362966A (en) * | 2018-02-12 | 2018-08-03 | 广东电网有限责任公司电力科学研究院 | A kind of oil-immersed type transformer high-precision noise on-line monitoring method and system |
CN108760030A (en) * | 2018-03-22 | 2018-11-06 | 国网湖南省电力有限公司 | A kind of self-cooled transformer iron core and winding noise level detection method and system |
CN109443523A (en) * | 2018-10-23 | 2019-03-08 | 国网天津市电力公司电力科学研究院 | A kind of 35kV is indoor or underground substation structural noise indirect measurement method |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013068517A (en) * | 2011-09-22 | 2013-04-18 | Jfe Steel Corp | Noise prediction method for power transformers |
CN104359549A (en) * | 2014-11-27 | 2015-02-18 | 国家电网公司 | Transformer core vibration noise analysis method |
CN104951635A (en) * | 2014-03-24 | 2015-09-30 | 国家电网公司 | Transformer noise estimating method for determining full load noise |
CN104964738A (en) * | 2015-06-10 | 2015-10-07 | 国家电网公司 | Noise separation-based power transformer noise assessment method |
CN105095609A (en) * | 2015-09-21 | 2015-11-25 | 武汉大学 | Transformer electromagnetic vibration noise calculating method based on finite element method |
-
2016
- 2016-06-22 CN CN201610460417.5A patent/CN106153176B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013068517A (en) * | 2011-09-22 | 2013-04-18 | Jfe Steel Corp | Noise prediction method for power transformers |
CN104951635A (en) * | 2014-03-24 | 2015-09-30 | 国家电网公司 | Transformer noise estimating method for determining full load noise |
CN104359549A (en) * | 2014-11-27 | 2015-02-18 | 国家电网公司 | Transformer core vibration noise analysis method |
CN104964738A (en) * | 2015-06-10 | 2015-10-07 | 国家电网公司 | Noise separation-based power transformer noise assessment method |
CN105095609A (en) * | 2015-09-21 | 2015-11-25 | 武汉大学 | Transformer electromagnetic vibration noise calculating method based on finite element method |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106599395A (en) * | 2016-11-25 | 2017-04-26 | 国网上海市电力公司 | Numerical simulation calculation method for noise of oil immersed transformer |
CN108362966A (en) * | 2018-02-12 | 2018-08-03 | 广东电网有限责任公司电力科学研究院 | A kind of oil-immersed type transformer high-precision noise on-line monitoring method and system |
CN108760030A (en) * | 2018-03-22 | 2018-11-06 | 国网湖南省电力有限公司 | A kind of self-cooled transformer iron core and winding noise level detection method and system |
CN109443523A (en) * | 2018-10-23 | 2019-03-08 | 国网天津市电力公司电力科学研究院 | A kind of 35kV is indoor or underground substation structural noise indirect measurement method |
Also Published As
Publication number | Publication date |
---|---|
CN106153176B (en) | 2020-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Pelat et al. | The acoustic black hole: A review of theory and applications | |
Zhang et al. | A two-way interactive nesting procedure with variable terrain resolution | |
CN106153176A (en) | A kind of system and method for single-phase oil immersion type transformer noise is estimated | |
Xue et al. | Design, fabrication, and preliminary characterization of a novel MEMS bionic vector hydrophone | |
Barney et al. | Fluid flow in a dynamic mechanical model of the vocal folds and tract. I. Measurements and theory | |
Barre et al. | Experimental study of a normal shock/homogeneous turbulence interaction | |
CN108416127A (en) | Submarine high pressure blows down pipe-line system multiple physical field coupling simulation method and system | |
Hou et al. | Design and modeling of a magnetic-coupling monostable piezoelectric energy harvester under vortex-induced vibration | |
CN104408295B (en) | A kind of Loads of Long-span Bridges substructure stormy waves coupling load method for numerical simulation | |
Suryadi et al. | Wall pressure spectra on a DU96-W-180 profile from low to pre-stall angles of attack | |
Lee et al. | Prediction of far-field wind turbine noise propagation with parabolic equation | |
Hafizh et al. | Vortex induced vibration energy harvesting using magnetically coupled broadband circular-array piezoelectric patch: Modelling, parametric study, and experiments | |
Robertson et al. | Planar analysis of a quasi-zero stiffness mechanism using inclined linear springs | |
Lai et al. | Nonhydrostatic three-dimensional method for hydraulic flow simulation. II: Validation and application | |
Zhang et al. | Experimental investigation of tip vortex formation noise produced by wall-mounted finite airfoils | |
CN110543677A (en) | vortex characteristic driven rotational turbulence PANS model | |
Staubs | Real airfoil effects on leading edge noise | |
Wang et al. | A miniature hydro-energy generator based on pressure fluctuation in Kármán vortex street | |
Zhou et al. | Wind loads and wind-induced responses of Guangzhou New TV Tower | |
Liu et al. | Effect of a variable-bend slat on tones due to the cove’s self-excited oscillation | |
Gatti | Spatially-varying multi-degree-of-freedom electromagnetic energy harvesting | |
Yuksel et al. | Active noise control in a duct with flow | |
Di Lorenzo | Operational Modal Analysis for rotating machines | |
Chakrabarti | Studying shocks in model astrophysical flows | |
Cotte et al. | Towards a semi-empirical trailing edge noise model valid for attached and separated turbulent boundary layers |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |