基于ANSYS和SYSNOISE的电机噪声仿真分析方法
第15卷第9期
201
电机与控制学报
CONTROL
VoL15No.9
1年9月
ELECTRICMACHINESAND
Sep.2011
基于ANSYS和SYSNOISE的电机噪声仿真分析方法
崔淑梅,
于天达,
宋立伟
(哈尔滨工业大学电气工程及自动化学院,黑龙江哈尔滨150001)
摘要:针对电机噪声分析时有限元软件ANSYS不能直接得到声响特性与实验数据比对的问题,研
究与振动声学分析软件SYSNOISE进行联合仿真的方法。以一台永磁直流电机为例,建立基于AN—
SYS的结构动力模型,计算振动频响特性并作为输入边界激励条件加入SYSNOISE,用SYSNOISE计算出声压频响特性。仿真与实验结果表明,最大相对误差为8.5%,满足工程设计需要。给出了AN.
SYS与SYSNOISE间结构网格和边界条件预处理与导入的具体办法,解决了两个软件的兼容问题。该方法利用两个软件的优势加快了计算速度,仿真结果亦直观,方便与实验结果相比较。关键词:电机;噪声;有限元;边界元;软件接口
中图分类号:TM
351
文献标志码:A文章编号:1007—449X(2011)09—0063—05
Simulationmethodofmotoracousticradiationbased
ANSYSandSYSNOISE
CUIShu—mei,
YUTian—da,
SONGLi—wei
on
(SchoolofElectricalEngineeringandAutomation,Harbin
Institute
ofTechnology,Harbin150001,China)
Abstract:Due
to
belackoftheacousticradiationcharacteristicsandexperimentaldatacompafisioninthe
noiseanalysisofmotorbyANSYS—afiniteelementmethod(FEM)software,the
paperpresented
a
co—sim—
ulationmethodwithANSYSandSYSNOISE-aboundaryelementdynamicmodelof
a
method(BEM)software.The
was
structural
permanentmagnetdirectcurrent(PMDC)motor
as
built,andthevibrationfrequen・
conditionsforcalculating
cyresponseanalysisresultsweretransferredtotheSYSNOISE
the
boundary
theacousticpressurefrequencyresponsefunction.Withthecomparisionofsimulationandexperimentaldata.themaximumrelativeerroris8.5%andtheengineerdesignrequirementsweremet.Thepre-pro・cessingandimportofstructuralgridsand
typroblemsofthetwosoftwareswere
boundary
conditions
were
especiallypresented,andcompatibili・
was
solved.The
calculationspeed
are
acceleratedwiththe
method
by
theadvantagesofthe
two
softwares.Thesimulationresults
intuitiveandconvenienttobecompared
withtheexperimentalresults.Keywords:motor;noise;finite
element;boundaryelement;softwareinterface
o引言
电机气隙磁场引起的电磁激振力作用于电机结
收稿日期:2010—05一19
基金项目:国家自然科学基-全(51077024)
筹言盖蓄_;;:并雾茎萎鲨然蠹徽
采用解析法计算电机噪声辐射He,通常将其近似等
作者简介:崔淑梅(19“一)。女。博士.教授,博士生导师,研究方向为电机及控制等;
于天达(1985一),男.硕士研究生。研究方向为永磁电机的振动与噪声;宋立伟(1970一),男.博士.教授,博士生导师。研究方向为电机设计与工艺等。
万方数据
电机与控制学报
第15卷
效为平面辐射器、无限长圆柱形辐射器或球形辐射器¨。21。然而,实际电机结构复杂,这种方法误差较大,也不能精细地分析与优化电机。数值法可解决复杂结构的计算问题【3。5j。
ANSYS是基于有限元法的一个分析软件,可以分析电机的电磁学特性和热力学特性,也可进
行结构振动分析,但用于电机的噪声分析时,只能
得到振动特性,不能进一步得到电机的声响特性。
尤其是实验一般测得的是电机声响特性,仿真结
果不能直接与实验结果相对应,需要进一步处理。SYSNOISE是一个专业振动与声学分析软件,可以直接得到声响特性,目前其主要应用于船舶、管道等领域哺。8j。
SYSNOISE善于处理复杂结构的振动声学问
题,但没有有限元和边界元网格生成前处理功能;
ANSYS善于进行网格前处理以及结构振动分析。
因此进行复杂结构噪声数值分析时,就需要采用ANSYS和SYSNOISE进行联合仿真,如文献[9]研
究了处理弹性结构水下振动和声辐射问题时AN—
SYS和SYSNOISE之间的数据接口技术,文献[10]
将两软件联合成功地对鱼雷振动及声学特性进行了分析。
由于电机与上述两系统的结构和激振力产生原因不同,因此,以一台永磁直流电机为例,从电机有限元模型的建立、电磁激振力的计算、模型网格导人、边界条件的加载几个方面,研究将ANSYS和
SYSNOISE联合进行电机噪声数值分析的具体
方法。1
基于SYsNOIsE电机噪声分析方法
以某汽车发动机冷却风扇用4极20槽150
W
永磁直流电机为例,探讨结合ANSYS和SYSNOISE
软件平台进行声辐射仿真分析的方法。1.1分析流程
采用SYSNOISE和ANSYS进行声学计算分析
的流程如图l所示。
1)在ANSYS中进行三维建模,并用实体单元(Solid186)和壳单元(shell91)进行网格划分,保存为木.cbd文件。
2)在ANSYS中进行结构的谐响应分析,将振动位移数据保存为宰.fre文件。
3)将网格模型(}.cbd)和结构振动边界条件
(木.fre)导人SYSNOISE,并设置自由边界条件、流
体属性,求解频率范围和步长。
4)声学响应求解。
万方数据
1
图1
ANSYS与SYSNOISE联合仿真流程
Fig.1
Theflowchart
of∞-simulation
in
ANSYSandSYSNOIsE
在图l所示的仿真流程步骤1)中采用实体单
的声学边界元模型;考虑到步骤2)中计算得到的电机的振动频响特性,在声学计算中要作为SYS-NOISE声学边界元模型的输入边界激励条件加入,因此将有限元模型进行抽壳处理,并对整机结构外声场计算结果影响较小的局部特征(如部分附件、
螺栓孔、机壳内永磁体等)进行进一步简化,但节点
不进行任何处理,以保持有限元模型与边界元模型节点的一致性,从而保证数据导人的正确性。另外,
在边界元的噪声辐射计算中,单元的细化程度由结
构总体尺寸和声学流体域中预计算得到的最高频率数值决定,所选单元特征尺寸必须小于由最高关心频率所决定的声波波长的1/6时,声学量结果的计算精确度方可满足要求【6】。综合考虑计算可以达
到的精确度、建模的复杂程度以及求解成功的可能性,本文中机壳选用四边形单元,端盖选择三角形单元实现边界元网格模型的离散,并确定声边界元网格所取的单元长度最大为5mm,以保证所要分析的模型频率精确度控制在661Hz以下。
图2(a)为整机实体模型,图2(b)为结构有限元分析中简化的有限元模型,图2(e)为经抽壳、单元翻转等操作后的边界元模型。边界元模型共有节点8928个,单元3436个,边界元模型能够满足高频计算的要求。
1.2模型建立方法
元网格划分技术,而边界元不能采用实体单元,因此
必须对原有限元模型进行适当处理才可转化成相应
第9期崔淑梅等:基于ANSYS和SYSNOISE的电机噪声仿真分析方法
65
将ANSYS分析所得的位移边界条件通过SYS・
蛾◆
(a)实体模型
(b)有限元横型
(c)边界兀模型
NOISE预设的与有限元软件间的接口导人边界元计
算模型中,即完成了主要边界条件的施加。但SYS-
NOISE能够识别的数据格式与有限元计算结果存在一定的差别,需要编制格式转换程序将响应计算结果转换为SYSNOISE能够识别的数据文件。图4为
格式转换程序框图。SYSNOIsE会根据有限元和边
界元网格间的关系,利用转换后的数据文件自动探测耦合表面,实施正确的模型边界条件传递,这样便实现了ANSYS与sYSNOISE的接力分析。
图2声辐射分析模型的建立
ng.2
EstabUshmenlofmodelforacousticradj血fion
1.3边界条件施加方法
作用在永磁直流电机上的激振力所产生的位移
响应通过ANSYS计算得到。所分析的小型永磁直流电机的噪声主要由电磁激振力产生,其来源有径向电磁力、转矩脉动和不平衡磁拉力,其中径向电磁力主要是由定子主波磁场与转子齿谐波磁场产生;转矩脉动为永磁体与电枢齿间相互作用力的切向分量的波动以及绕组感应电流纹波所引起;不平衡磁
拉力主要由装配误差所致,通常通过工艺手段可以
控制在允许的范围内,与前两项相比较小,计算时暂
不考虑。
图4数据文件转换程序框围
Fig.4
通过分析可知,4极20槽150W永磁直流电机径向电磁力和脉动转矩的主要谐波次数为20次、40次和60次,其中加次谐波成分最大,而且齿槽转矩谐波幅值与转速无关,而电磁转矩谐波幅值随转速正比下降,电磁转矩各次谐波与齿槽转矩相应谐波发生矢量叠加,导致20次谐波彼此增强,成为主要成分。电机电磁力的20次谐波作用在结构上产生
振动变形为4阶,图3为电磁力谐响应分析所得的
erogramflowchartoffileconversion
2仿真计算及结果分析
建立上述边界元模型并施加边界条件后,就可
以利用SYSNOISE软件计算电机的声学特性。图5
为电机机壳表面振速云纹图,表示机壳表面振动速度的分布情况。图6(a)为球面空间的声压云纹图,
黝敢瓯哙
F咤.3
图3电机振型的有限元仿真结果
FEMresultsofmodeshap8forthemotor
rig.5
■蓦_
图5机壳表面振速(430Hz)
Vibration印eedof
_
shell∞rf妇(430Hz
万方数据
电机与控制学报
通过图6(a)可以得到电机在其周围空间噪声辐射
第15卷
的分布情况。计算所选场点为图6(b)所示的球心
位于电机机壳结构中心,半径为0.1m的球面场点。
3实验验证
实验在半消声室中进行,噪声测试装置为BK2238噪声计和南京安正CRAS信号分析系统。为了减少计算时问,仿真时建立的是半径为0.1m的球面并计算该面上的噪声值,实验时测试的也是
距离中心0.1m的相应点。实验中使电机分别在
400~2000
然后可以求取图中任何一点的声压频响函数。图7为图6(a)中场点122处的声压频响函数,通过图7
可以得到电机在其额定转速范围内变化时的噪声辐
射情况。声学频响分析中采用窄频程,计算频率为
133~6612000
Hz,步长为33Hz(对应电机转速为400—r/rain(步长为100r/min)的各转速工作
r/rain,步长为100r/rain)。点上稳定旋转,测试得到了各转速点的噪声频谱。
“…’_…
帆删_¨m船
利用所建立的仿真方法,计算电机的的噪声频谱,图8为举例点600r/rain(对应430Hz)时的情况。求出20次电磁力波所激发的噪声幅值,将各转速点所对应的20次电磁力波所激发的噪声幅值综合到一起,即可得到电机在额定转速范围内变化时噪声的实际辐射情况,如图9所示。与仿真结果相比,实测值与仿真结果存在一定的误差,最大相对误
…一..+。
差为8.5%。产生误差的主要原因为:电机在实际运行中20次力波除了含有文中所给定的径向电磁力以外还有齿槽转矩等力波成分;另外,实验所测得
刷果吲是
。。,-¨…
。缫黜幕
№.6
逛晰“.№甏
一
O
然僳刘
包计接h
含算受_|■电结。.~_
嗍
州惭-h一一一
一
声野等计
一
2
一一一一一1
600
图6场点网格模型及声压(430Hz)
Meshmodelandacousticpressureof
field
频率/kHz
图8
Fig.8
point(430Hz)
r/rain时噪声频谱图
Noiseflpecltrogramat600r/min
/
/
罩出
/
/
,
,一
{出
一
/厂
—F
P
/
韫
,
忸
,少
+仿直结果
图7场点122处频响函数
_毡.7
Frequency
图9不同转速下噪声实测与仿真结果对比
point
122
Fig.9
Thecomparisonofexperiment
response
f.mctionoffield
andsimulation
万方数据
第9期崔淑梅等:基于ANSYS和SYSNOISE的电机噪声仿真分析方法
67
4
结语
本文研究了利用ANSYS与SYSNOISE联合进
行永磁直流电机全三维整体模型振动声学数值计算
的方法。建立了基于有限元技术的结构动力计算模
型。提取了基于边界元技术的噪声预测模型,给出了
能将动力响应计算结果转换为SYSNOISE可识别的
数据文件的转换方法。计算输出了机壳表面振速云纹图、场点声压云纹图和场点声压的频响函数,联合仿真直接给出了噪声频谱和噪声辐射结果。通过仿
真结果与实验结果对比,得到最大相对误差为
8.5%,计算精确度满足了工程设计需要。
该方法既利用了ANSYS较容易建立电机结构模型,得到振动特性的优点,又发挥了SYSNOISE可
以快速得到声响特性的优势,加快了计算速度,仿真
结果亦直观,方便与实验结果相比较。
参考文献:
[1]陈世坤.电机设计[M].北京:机械工业出版社,2000:191
—202.
[2]陈永校,诸自强,应善成.电机噪声的分析和控制[M].杭州:浙江大学出版社。1987:58—67.
[3]LIUzS。LUC,WANGYY,etai.Predictionofnoiseinside
tracked
vehicles[J].AppliedAcoustics,2006,67(1):74—91.
[4]CITARELLA
R,FEDERICOL,CICATIE[.LOA.Modalacoustic
transfer
vecior
approachin
a
FEM・・BEM
vibro-acousticanalysis
【J].E|lgineering
Analysis
with
Boundary
Elements。2007,31
(3):248—258.
AF,HAMILTONDA,HAYESPA.Predictionof
ra-
万方数据
diatednoisefrommachinecomponentsusingtheBEMandthemy—
leish
integral[J].NoiseControlEngineering
Journal。1998,46
(3):77—82.
[6]杨德庆,郑靖明,王德禹,等.基于SYSNOISE软件的船舶振动
声学数值计算[J].中国造船,2002,43(4):32—37.
YANGDeqing,ZHENG
Jingming,WANGDeyu,eta1.Numeri—
cat
analysisofvibro—acousticcharactersofshipwith
SYSNOISE
software[J].Shipbuildingof
China,2002,43(4):32—37.
[7]潘立志,粱亚坤,占敏,等.基于SYSNOISE的轴流式风机管道
噪声辐射仿真分析[J].机械与电子,2009,(9):32—35.
PANLizhi.LIANGYal【un。ZHANMin,ela1.Acousticradiationsimulationofair—intakeandair-exhaustpipeof
an
axis・flowfan
based
on
SYSNOISE[J].Machinery&Electronics,2009,(9):
32—35.
研究[J].机电产品开发与创新,2007,20(4):1—2.
LIWei,YANGJing,CHENBing,eta1.Simulationandanalysis
of
acou鲥cradiationofpipingsystembased
on
SYSNOISE[J].
Development&Innovationof
Machinery&ElectricalProducts.
2007,20(4):1—2.
究[J].应用科技,2004,31(8):32—34.
WANG
Jing。SHANGDejiang.Studyoftheinterfacebet’Ireen
ANYSYSand
SYSNOISE[J].Applied
Science
andTechnology,
2004,31(8):32—34.
肖汉林,于俊卫。。张瑞斌,等.鱼雷电机一艉轴系统振动与声辐射特性分析[J].鱼雷技术.2005。13(4):33—36.
XIAOHanlin,YUJunwei,ZHANGReibin,eta1.Research
on
vibrationandacou鲥eradiationcharacteristicoftorpedoelectric
motor
andstemshaftsystem[J].TorpeIloTechnology,2005,13
(4):33—36.
(编辑:于双)
[8]李伟.杨竞,陈兵,等.基于SYSNOISE管道声辐射仿真分析及
[9]王晶,商德江.ANSYS和SYSNOISE之问的数据接口技术研
[10]
[5]SEYBERT
基于ANSYS和SYSNOISE的电机噪声仿真分析方法
作者:作者单位:刊名:英文刊名:年,卷(期):
崔淑梅, 于天达, 宋立伟, CUI Shu-mei, YU Tian-da, SONG Li-wei哈尔滨工业大学 电气工程及自动化学院,黑龙江哈尔滨,150001电机与控制学报
Electric Machines and Control2011,15(9)
参考文献(10条)
1. 王晶;商德江 ANSYS和SYSNOISE之间的数据接口技术研究[期刊论文]-应用科技 2004(08)
2. 李伟;杨竞;陈兵 基于SYSNOISE管道声辐射仿真分析及研究[期刊论文]-机电产品开发与创新 2007(04)3. 潘立志;梁亚坤;占敏 基于SYSNOISE的轴流式风机管道噪声辐射仿真分析[期刊论文]-机械与电子 2009(09)4. LIU Z S;LU C;WANG Y Y Prediction of noise inside tracked vehicles[外文期刊] 2006(01)5. 陈永校;诸自强;应善成 电机噪声的分析和控制 1987
6. 肖汉林;于俊卫;张瑞斌 鱼雷电机-艉轴系统振动与声辐射特性分析[期刊论文]-鱼雷技术 2005(04)7. 陈世坤 电机设计 2000
8. 杨德庆;郑靖明;王德禹 基于SYSNOISE软件的船舶振动声学数值计算[期刊论文]-中国造船 2002(04)
9. SEYBERT A F;HAMILTON D A;HAYES P A Prediction of radiated noise from machine components using the BEM and therayleigh integral[外文期刊] 1998(03)
10. CITARELLA R;FEDERICO L;CICATIELLO A Modal acoustic transfer vector approach in a FEM-BEM vibro-aconstic analysis 2007(03)
本文链接:http://d.g.wanfangdata.com.cn/Periodical_djykzxb201109011.aspx