CN106599395A - 一种油浸变压器的噪声数值仿真计算方法 - Google Patents
一种油浸变压器的噪声数值仿真计算方法 Download PDFInfo
- Publication number
- CN106599395A CN106599395A CN201611052939.8A CN201611052939A CN106599395A CN 106599395 A CN106599395 A CN 106599395A CN 201611052939 A CN201611052939 A CN 201611052939A CN 106599395 A CN106599395 A CN 106599395A
- Authority
- CN
- China
- Prior art keywords
- oil
- type transformer
- immersed type
- model
- modeling
- 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.)
- Pending
Links
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Abstract
本发明涉及一种油浸变压器的噪声数值仿真计算方法,其特征在于,所述方法包括下列步骤:对油浸变压器进行声学建模;对建模后的油浸变压器的油箱壁的每个面施加荷载,采用谐响应分析进行求解;对计算结果进行可视化处理,显示油浸变压器的噪声声场的SPL分布。与现有技术相比,本发明具有求解简单、计算准确以及结果直观等优点。
Description
技术领域
本发明涉及变压器领域,尤其是涉及一种油浸变压器的噪声数值仿真计算方法。
背景技术
噪声是由空气粒子的振动引起的,而粒子振动则是主要由油箱壁振动激励,因此理论上只要通过测量得到油箱壁所有点的振动,就可以计算出空气中的噪声水平,其计算原理具体如下:
其中,ρ0为空气的密度,αn(r')为箱壁上矢径长度为r'点处的法向振动加速度,G(r,r')为格林函数,表示矢径r'处的单位源在矢径r处的声压。
在实际操作中,不可能得到所有点的振动,必须对以上模型进行相应的简化。同时,以上计算比较复杂,得到解析解比较困难,实用性能较差。
发明内容
本发明的目的是针对上述问题提供一种油浸变压器的噪声数值仿真计算方法。
本发明的目的可以通过以下技术方案来实现:
一种油浸变压器的噪声数值仿真计算方法,所述方法包括下列步骤:
1)对油浸变压器进行声学建模;
2)对建模后的油浸变压器的油箱壁的每个面施加荷载,采用谐响应分析进行求解;
3)对步骤2)的计算结果进行可视化处理,显示油浸变压器的噪声声场的SPL分布。
所述声学建模包括结构模型建模、流体模型建模和结构模型与流体模型之间界面的建模。
所述步骤1)具体为:
11)确定声学建模的单元类型和单元尺寸;
12)根据步骤11)确定的单元类型和单元尺寸建立的油浸变压器油箱的外边界模型;
13)通过对称得到整体的油浸变压器油箱的外边界模型;
14)在油浸变压器油箱的外边界模型外侧建立空气模型。
所述单元尺寸不大于声波波长的
所述空气模型包括远场空气模型和近场空气模型。
所述步骤2)具体为:
21)对油浸变压器油箱壁的外边界模型的每个面施加相应速度;
22)对空气模型的外表面施加压力为0的约束。
所述相应速度具体为:
211)在油浸变压器油箱壁的每个面上选取测量点,测量测量点的振动速度;
222)计算测量点的振动速度的平均值,得到油浸变压器油箱壁的外边界模型的每个面的相应速度。
所述测量点的数量为6~12个。
与现有技术相比,本发明具有以下有益效果:
(1)通过对油浸变压器进行声学建模,实现对油浸变压器的简化,通过对模型施加相应的速度,计算得到油浸变压器的周围噪声值,与传统的通过公式计算噪声的方法相比,计算简便,实用性能强。
(2)建模时首先建立的油浸变压器油箱的外边界模型,再通过对称得到整体的油浸变压器邮箱的外边界模型,比起直接建立整体模型相比,操作简便。
(3)建立的单元尺寸不大于声波波长的便于模型的分析和计算。
(4)通过在油浸变压器油箱壁的每个面上选取测量点,测量测量点的振动速度作为仿真计算所施加的相应速度,大大简化了计算的难度,同时也可以通过增加测量点的数目提高计算结果的准确性,实用性能高。
附图说明
图1为本发明的方法流程图;
图2为声学建模所建立的模型示意图;
图3为施加荷载的示意图;
图4为油浸变压器噪声声场的SPL分布图,其中(4a)为100Hz时的噪声声场的SPL分布图,(4b)为200Hz时的噪声声场的SPL分布图;
图5为与油箱壁相距0.3米处的声压计算值与试验比较的结果图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
如图1所示,本实施例的目的是提供一种油浸变压器的噪声数值仿真计算方法,包括以下步骤:
1)对油浸变压器进行声学建模;
2)对建模后的油浸变压器的油箱壁的每个面施加荷载,采用谐响应分析进行求解;
3)对步骤2)的计算结果进行可视化处理,显示油浸变压器的噪声声场的SPL分布。
本实施例中,通过ANSYS软件实现上述步骤,ANSYS是一个大型商用有限元软件,融结构、流体、电场、磁场、声场分析于一体,在声场分析中,既可以考虑结构-流体耦合的情况,也可以考虑不耦合的情况。其基本过程包括建模、施加载荷并求解和对结果的后处理,具体过程如下:
声学建模包括结构模型和流体模型及两者之间的界面,在本实施例中不考虑结构-流体耦合,所以只用建立流体模型。
先建立变压器油箱的外边界,再建立空气部分,包括远场和近场,最后通过布尔操作进行分开。在建模过程中可以先建立1/4模型,再通过对称得到整体模型,如图2所示。
本实施例中,单元采用FLUID221,为10节点四面体单元,对近场部分keyopt(4)=0,远场keyopt(4)=1。单元尺寸与声波波长即频率有关,一般每波长不少于5个单元,本项目中单元尺寸取波长的1/12。对空气,取密度1.21kg/m3,声速343m/s。
对油箱壁的每个面,按该频率下各点速度的平均值近似认为该面整体振动,在油箱壁前后左右四个面的节点上施加相应的速度。该相应的速度,是通过对变压器油箱壁的每个面,测量6~12个点上的振动速度,并认为其平均值能代表整个面的振动平均水平,近似认为该面以相同的速度对空气进行激励而确定的。空气外表层则施加压力为0的约束,如图3所示,由于只对单个频率进行分析,采用谐响应分析进行求解。
在后处理中可以对计算结果进行可视化处理。在本项目中,可以看到空气场中的压力分布和SPL值,既可以用云图显示也可以用等截面图显示。如图4所示,分别为100Hz和200Hz时声场的SPL分布。
在后处理中还可以读出麦克风位置处的声压,以便与试验结果进行对比。图5为油箱正面中间,与油箱壁相距0.3米处的声压计算值与试验比较的结果。
从以上数值计算的结果可以看出,采用以上简化方法可以得到变压器周围声场的分布,与试验值基本吻合,证明了该方法的合理性。
但在个别点上误差较大(300Hz),这是由于油箱壁上可能产生局部振动,此时采用部分点的测量值来代替整个油箱壁具有较大误差。增加测量点的分布应该能提高数值计算的精度。
Claims (8)
1.一种油浸变压器的噪声数值仿真计算方法,其特征在于,所述方法包括下列步骤:
1)对油浸变压器进行声学建模;
2)对建模后的油浸变压器的油箱壁的每个面施加荷载,采用谐响应分析进行求解;
3)对步骤2)的计算结果进行可视化处理,显示油浸变压器的噪声声场的SPL分布。
2.根据权利要求1所述的油浸变压器的噪声数值仿真计算方法,其特征在于,所述声学建模包括结构模型建模、流体模型建模和结构模型与流体模型之间界面的建模。
3.根据权利要求2所述的油浸变压器的噪声数值仿真计算方法,其特征在于,所述步骤1)具体为:
11)确定声学建模的单元类型和单元尺寸;
12)根据步骤11)确定的单元类型和单元尺寸建立的油浸变压器油箱的外边界模型;
13)通过对称得到整体的油浸变压器油箱的外边界模型;
14)在油浸变压器油箱的外边界模型外侧建立空气模型。
4.根据权利要求3所述的油浸变压器的噪声数值仿真计算方法,其特征在于,所述单元尺寸不大于声波波长的
5.根据权利要求3所述的油浸变压器的噪声数值仿真计算方法,其特征在于,所述空气模型包括远场空气模型和近场空气模型。
6.根据权利要求1所述的油浸变压器的噪声数值仿真计算方法,其特征在于,所述步骤2)具体为:
21)对油浸变压器油箱壁的外边界模型的每个面施加相应速度;
22)对空气模型的外表面施加压力为0的约束。
7.根据权利要求6所述的油浸变压器的噪声数值仿真计算方法,其特征在于,所述相应速度具体为:
211)在油浸变压器油箱壁的每个面上选取测量点,测量测量点的振动速度;
222)计算测量点的振动速度的平均值,得到油浸变压器油箱壁的外边界模型的每个面的相应速度。
8.根据权利要求7所述的油浸变压器的噪声数值仿真计算方法,其特征在于,所述测量点的数量为6~12个。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611052939.8A CN106599395A (zh) | 2016-11-25 | 2016-11-25 | 一种油浸变压器的噪声数值仿真计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611052939.8A CN106599395A (zh) | 2016-11-25 | 2016-11-25 | 一种油浸变压器的噪声数值仿真计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106599395A true CN106599395A (zh) | 2017-04-26 |
Family
ID=58593334
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611052939.8A Pending CN106599395A (zh) | 2016-11-25 | 2016-11-25 | 一种油浸变压器的噪声数值仿真计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106599395A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107315886A (zh) * | 2017-07-06 | 2017-11-03 | 国网重庆市电力公司电力科学研究院 | 一种变压器室外部三维空间噪声预测的方法与装置 |
CN110059397A (zh) * | 2019-04-12 | 2019-07-26 | 国网湖南省电力有限公司 | 一种低噪声油浸式配电变压器油箱设计方法、系统及介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104359549A (zh) * | 2014-11-27 | 2015-02-18 | 国家电网公司 | 一种变压器铁芯振动噪声分析方法 |
CN104951635A (zh) * | 2014-03-24 | 2015-09-30 | 国家电网公司 | 一种确定满负荷噪声的变压器噪声预估方法 |
CN104964738A (zh) * | 2015-06-10 | 2015-10-07 | 国家电网公司 | 一种基于噪声分离的电力变压器噪声评估方法 |
CN105095609A (zh) * | 2015-09-21 | 2015-11-25 | 武汉大学 | 基于有限元法的变压器电磁振动噪声计算方法 |
CN106153176A (zh) * | 2016-06-22 | 2016-11-23 | 中国电力科学研究院 | 一种用于对单相油浸式变压器噪声进行预估的系统及方法 |
-
2016
- 2016-11-25 CN CN201611052939.8A patent/CN106599395A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104951635A (zh) * | 2014-03-24 | 2015-09-30 | 国家电网公司 | 一种确定满负荷噪声的变压器噪声预估方法 |
CN104359549A (zh) * | 2014-11-27 | 2015-02-18 | 国家电网公司 | 一种变压器铁芯振动噪声分析方法 |
CN104964738A (zh) * | 2015-06-10 | 2015-10-07 | 国家电网公司 | 一种基于噪声分离的电力变压器噪声评估方法 |
CN105095609A (zh) * | 2015-09-21 | 2015-11-25 | 武汉大学 | 基于有限元法的变压器电磁振动噪声计算方法 |
CN106153176A (zh) * | 2016-06-22 | 2016-11-23 | 中国电力科学研究院 | 一种用于对单相油浸式变压器噪声进行预估的系统及方法 |
Non-Patent Citations (1)
Title |
---|
曹枚根 等: "大型电力变压器振动与噪声数值分析", 《中国电业(技术版)》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107315886A (zh) * | 2017-07-06 | 2017-11-03 | 国网重庆市电力公司电力科学研究院 | 一种变压器室外部三维空间噪声预测的方法与装置 |
CN110059397A (zh) * | 2019-04-12 | 2019-07-26 | 国网湖南省电力有限公司 | 一种低噪声油浸式配电变压器油箱设计方法、系统及介质 |
CN110059397B (zh) * | 2019-04-12 | 2023-02-07 | 国网湖南省电力有限公司 | 一种低噪声油浸式配电变压器油箱设计方法、系统及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kawai et al. | Wall-modeling in large eddy simulation: Length scales, grid resolution, and accuracy | |
CN112765736B (zh) | 一种高超声速钝前缘绕流湍动能入口边界设置方法 | |
Clark et al. | Requirements and challenges for CFD validation within the High-Lift Common Research Model ecosystem | |
CN114662425B (zh) | 一种水轮机启停工况流场仿真预测方法及系统 | |
CN109753681B (zh) | 一种水下流激空腔线谱噪声的综合预报方法 | |
Diskin et al. | Grid convergence for turbulent flows | |
Vatsa et al. | Application of FUN3D solver for aeroacoustic simulation of a nose landing gear configuration | |
CN105468851A (zh) | 一种确定飞机动态重量特性的方法 | |
CN106599395A (zh) | 一种油浸变压器的噪声数值仿真计算方法 | |
Vatsa et al. | Aeroacoustic simulation of nose landing gear on adaptive unstructured grids with FUN3D | |
CN109190232B (zh) | 一种飞机平尾区动能损失计算评估方法 | |
CN110489832A (zh) | 一种用于湍流控制屏单元体气动性能的仿真试验方法 | |
CN108460204A (zh) | 一种通过扬声器振动部件的受力和位移反推其材料动态力学参数的方法 | |
CN116341421B (zh) | 高超声速流场数值模拟方法、系统、电子设备及存储介质 | |
CN102270249B (zh) | 识别零部件特征频率的方法 | |
CN110543677A (zh) | 一种涡特征驱动的旋转湍流pans模型 | |
CN104778321B (zh) | 一种复杂形状声学空间的分析预报方法 | |
Clapp et al. | Validating J-factor as a predictive method for the repeatability of aircraft store separation from unsteady cavity environments | |
Chwalowski et al. | Shock-Buffet Prediction Report in Support of the High Angle Working Group at the Third Aeroelastic Prediction Workshop | |
Szabó et al. | Advanced flutter simulation of flexible bridge decks | |
Kempf et al. | Zonal Hybrid Computational Aeroacoustics Simulation of Trailing Edge Noise Using a High-Order Discontinuous Galerkin Method | |
Zastawny et al. | Validation of intermittency-based transition prediction models | |
Robbins | Detailed Characterization of Flowfields and Uncertainty in a Speed-Bump Turbulent Separated Flow Validation Experiment | |
Jiang et al. | Experimental and Numerical Investigation on the External Aerodynamic Noise of High-Speed Train. | |
Graves et al. | Comparison of computational predictions with experimental data for supersonic cavity geometries |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170426 |
|
RJ01 | Rejection of invention patent application after publication |