ansys声学分析.docx

什么是声学声学研究声压波在流体介质中的产生、传播、吸收和反射。声学有如下的应用：声声纳一声学上雷达的对应物 设计音乐厅，希望声压均匀分布。 减小机器厂房内的噪音 汽车中的噪声消除 水下声学 设计扬声器、音箱、声滤、消音器及其他类似装置。 地球物理探测声场分析的类型只有在 ANSYS/Multiphysics 和 ANSYS/Mechanical 中能进行声场分析，通常包括对流体介质及其周围结构的建模。典型感兴趣的是不同频率的声波在流体中的压力分布、 压力梯度、 粒子速度、 声压级及声波的散射、 衍射、 传输、 辐射、衰减和散射。耦合的声场分析将考虑流体- 结构的相互作用。非耦合的声场分析模型只考虑流体而忽略任何流体-结构的相互作用。ANSY邻序假定流体是可压的,但只允许压力与平均压力相比有较小的变化。而且，流体假定为非流动并且无粘的（即粘性不引起耗散作用）。假定平均密度和平均压力不变，压力求解偏离平均压力而不是绝对压力。求解声学问题通过执行一个谐波响应分析可以解决许多声学问题。分析计算流体-结构界面上的谐波载荷 （正弦变化） 引起流体中的压力分布。 通过指定载荷的频率范围，可以观察到在不同的频率时压力的分布。 可以执行模态和瞬态的声学分析。 （参见 ANSYS Structural Analysis Guide 中关于这种分析更详细的叙述。）谐波声场分析的过程包括以下三个主要步骤：建立模型。施加边界条件和载荷并获得求解。查看结果。建立模型在此步骤中，用户指定工作名称和分析标题，然后用PREP砌处理器定义单元类型，单元实常数，材料属性和模型几何尺寸。这些任务与多数分析相同，在 ANSYS Basic Analysis Guide 中有叙述。谐波声场分析准则对一个谐波声场分析，考虑以下几点：单元类型一ANSY*场分析指定了四种单元类型：对二维和三维模型的流体 部分分别使用 Fluid29 和 Fluid30 单元，Fluid129 和 Fluid130 与 FLUID29和 FLUID30单元一起使用，用来构造包围 Fluid29和Fluid30单元的无限外壳。禾 用这些单元单型可以构造流体部分的模型， 然后利用相应的结构单元（ PLANE4、2SOLID45）构造固体模型。只有Fluid29和Fluid30单元才能与结构单元相接 触（在结构的内部或外部）；Fluid129和Fluid130单元只能与Fluid29和Fluid30 单元相接触，而不能直接与结构单元接触。FLUID29 与 FLUID30 单元对与固体相接触的声单元，要确保使用KEYOPT（2）=0缺省的设置允许流体-结构的相互作用。UX, UY UZ和PRES乍为自由度引起单元矩阵的不对称。对所有 其它的声单元，设置KEYOPT2） =1,致使带有PRESS由度的单元矩阵的对称。（见图 5-1 ） 对称矩阵需要的内存和计算时间更少， 因此只要可能就应该使用它。关于流体-结构的相互作用的详细信息参见ANSYS,Inc. Theory Reference 。图 5-1 二维声模型的例子（流体在结构的内部）FLUID129 和 FLUID130 单元对无限的吸收压力波的声单元，模拟在FLUID29和FLUID30单元之外无限延 伸域的输出效果。FLUID129和FLUID130单元提供了第二级的吸收边界条件，所 以输出的压力波到达模型的边界以最小的反射吸收到流体域内。FLIUD129单元用来建立二维流体区域的边界和诸如线单元。FLIUD130单元用来建立三维流体区域的边界和诸如平面表面单元。材料属性一声单元需要密度 （DEN$和声速（SONC作为材料属性（FLUID129 和FLUID130只需要SONC。如果在流体-结构界面存在声的吸收，利用标记 MU 来指定边界导纳B （吸收系数）。值B通常由实验来测定。对结构单元，指定 杨氏模量（EX）、密度（DENS和泊松比（PRXYE NUXY。常常数一当用FLUID129和FLUID130单元时，里面的有限元网格边界必须是 圆形的（二维及轴对称）或球形的（三维），而且圆形或球形边界的有限区域的半径必须指定为实常的RAD （见图5-2）圆或球的中心也必须用实常数指定：R,3,RAD,X0,Y0!REAL set 3 for FLUID129R,3,RAD,X0,Y0,Z0 !REAL set 3 for FLUID130如果二维轴对称情况圆的中心坐标（ X0， Y0） 或三维情况球的中心坐标（ X0，Y0, Z0）不是通过实常数指定的，ANSY赧定中心为总体坐标系的原点。5-2 应用吸收单元的例子5-3 浸在水中的圆筒对模型进行网格划分下面为使用二维无限的声单元进行网格划分的典型程序。 对三维单元的程序与此相同。如果还有结构部分，结构单元必须与FLUID29单元相邻，而不要与无 限流体单元（FLUID129相邻。这个过程自动在有限区域的边界加入FLUID129单元。这里环形结构用PLANE4站构单元进行网格划分。与PLANE42I元接触的流体单元层用带有 UX 和UY自由度的FLUID29单元划分网格同时打开流体-结构接触面。流体的外层用 没有UX和UY自由度的FLUID29单元建模。在 X0=Y0=Oft放置FLUID129单元的 半径为（参见§节）。可用下列方法定义FLUID129单元：命令： ESURFGUI: Main Menu>Preprocessor>Create>Elements>Inf Acoustic步骤一：内部流体区域的网格划分用FLUID29单元对圆形或球形边界（PLANE42所包围的内部流体区域划分 网格。5-4 对流体区域划分网格步骤二：生成无限的声单元按以下步骤：1选择圆形或球形边界上的节点：命令： NSELGUI: Utility Menu>Select>Entities2.指定FLUID129单元作为与FLUID29单元相联系的无限流体单元。命令： TYPE, REALGUI: Main Menu>Preprocessor>Attributes>Default AttribsMain Menu>Preprocessor>Real Constants无限单元对低频及高频激励都执行得很好。 数值实验已确定吸收单元远离结构或振动源区域以外大约人放置能产生准确的结果。这里人=c/f压力波的主波长。 c 是流体中的声速(SON)，Cf 是压力波的主频。例如，对浸没在水中的圆盘或球壳的直径D,封闭边界的半径RADS少应为D/2+入。而且，对一般的声分 析，网格必须足够的细致以能分辨最小的主频。3在边界生成吸收单元(FLUID129)。命令： ESURFGUI: Main Menu>Preprocessor>Create>Elements>Inf Acoustic图 5-5 在边界加入吸收单元步骤三：指定流体 - 结构接触面指定流体 - 结构接触面：1 选择界面上的所有节点。命令： NSELGUI: Utility Menu>Select Entities2选择附属于这些节点上的流体单元。命令： ESELGUI: Utility Menu>Select>Entities3指定所选择的节点作为流体- 结构接触面上的节点。命令： SFGUI: Main Menu>Preprocessor>Loads>Apply>Fluid-Struct>On Nodes 注意 在对分析进行求解前确保重选所有的节点。图 5-6 指定流体 -结构界面施加载荷并获得求解在这个步骤里，用户定义分析类型和选项，施加载荷，指定载荷步选项，并开始有限元求解。下面几节解释如何做这些工作。步骤一：进入 SOLUTION解器通过选择GUI途彳全Main Menu>Solution 或执行/SOLU命令进入SOLUTION： 解器。步骤二：定义分析类型利用GUI途径或命令集，定义分析类型和分析选项。用下列方法定义分析类型：命令： ANTYPE,HARMICGUI: Main Menu>Solution>New Analysis必须选择新的分析 NEW Analysis 因为重启动对谐波响应分析无效。如果需要施加另外的谐波载荷，每次作一个新的分析（或用 ANSYS Basic AnalysisGuide中叙述的“部分求解”程序）。步骤三：定义分析选项用下列方法指定求解方法：命令： HROPTGUI: Main Menu>Solution>Analysis Options尽管全量、减缩或模态叠加方法都是可选的方法，选择全量方法因为它自己 就可以处理非对称矩阵。用下列方法定义求解列表格式：命令：HROUTGUI: Main Menu>Solution>Analysis Options这个选项决定ANSY9口何在打印输出（）中对谐波自由度结果进行列表。用下列方法指定方程求解器：命令：EQSLVGUI: Main Menu>Solution>Analysis Options可以选择波前求解器（缺省），Jocobi共腕梯度（JCG求解器，或不完全 的Cholesky共腕梯度（ICCG求解器。对大多数模型推荐使用 JCG求解器。步骤四：在模型上施加载荷由谐波分析的定义可知，假定任何施加的载荷随时间是简谐（正弦）变化的. 在声学分析中完整地指定一个谐波载荷， 通常需要两条信息：强迫力的频率和振 幅。振幅是载荷的最大值，如表5-2所示用命令指定振幅。强迫力的频率是谐波 载荷的频率（循环/时间）。可用HARFR命令或GUI途彳全Main Menu> Solution>Time/Frequenc>Freq & Substeps在随后的载荷步中指定它。参见§节步骤5。表5-1列出了所有对谐波声分析可用的载荷及定义、列表和删除它们的命令。 除了惯性载荷，用户既可以在实体模型（关键点、线和面）上也可以在有限元模 型（节点和单元）上定义载荷。对实体模型载荷与有限元载荷综合的讨论参见 ANSYSS本分析程序指南中加载与求解的有关章节。表5-1声学分析中可用的载荷载荷类型种类命令族菜单途径位移 (UX,UY,UZ),压力(PRES)约 束DMainMenu>Preprocessor>-Loads->Apply>Displacement or PotentialMain Menu>Solution>Apply>Displacement or PotentialMain力(FX,FX,FZ),Menu>Preprocessor>-Loads->力矩Apply>Force/Moment(MX,MY,MZ),力FMain流体载荷Menu>Solution>-Loads-Apply>Force/Moment压力(PRES)表MainMenu>Preprocessor>-Loads->阻抗（IMPD）面 载SFApply>load type流体-结构相互作用标记（FSI）荷Main Menu>Solution>Apply>load typeACEL,Main惯Menu>Preprocessor>-Loads->向性OMEGA,Apply>load type旋转等轨荷DOMEGA, CGLOCCGOMEGA, DCGOf/Main Menu>Solution>Apply>loadIRLFtype在一个分析中，用户可以施加、移去、操作或列表载荷用GUI施加载荷通过一系列的层叠菜单可以得到除列表（见下面）以外的所有载荷操作。从 求解菜单中选择操作（应用等），载荷类型（位移、力等），然后是施加载荷的 目标（关键点等）。例如，按以下GUI途径，给一条线施加位移载荷：Main Menu>Solution>Apply>Displacement>On Lines用GUI途径对载荷进行列表：Utility Menu>List>Loads>load type利用命令施加载荷表5-2列出在声分析中所有可用来加载的命令。表5-2在声分析中施加载荷的命令载荷类型实体模型或 有限元图元施加删除列表操作施加设 置位移压力实体模型关键占八、DKDKDELEKLISTDTRAN实体模型线DLDLDELEDLLISTDTRAN实体模型面DADADELEALISTDTRAN有限元节点DDDELEDLISTDSCALE:DCUM力力矩实体模型关键占八、FKFKDELEFKLISTFTRAN有限元节点FFDELEFLISTFSCALEFCUM压力阻抗流体结构相互 作用标记实体模型线SFLSFLDELIESFLLISTSFTRANSFGRAD实体模型面SFASFADEL5FALISTSFTRANSFGRAD有限元节点SFSFDELESFLISTSFSCALSFCUMSFGRAD有限元单元SFESFEDEL5FELISTSFSCALSFBEAMSFCUMSFFUNSFGRAD惯性ACEL CGLOCCGOMGDCGOMOMEGADOMEGAIRLF卜卜载荷类型位移（UX UY UZZ和压力（PRES这些是DOF自由度）约束。例如，在刚性的流体-结构接触面上指定零位移。 也可以指定非零位移，但应记住它们被假定为简谐的。在自由流面边界通常指定 零压力（此处流体不封闭，诸如一个开口）。力（FX, FY, FZ）和力矩（MX MY MZ通常在模型的固体部分指定载荷以激励流体。当得知一个节点处的声载荷，可以通过指定 FLOW：标记来指定流动载荷：命令：F,FLOWGUI: Main Menu>Solution>Apply>-Structural-Force/Moment>On Nodes压力（PRES可以在固体部分指定表面载荷来代替力和力矩。阻抗（IMPD实际上显示不是载荷而是表面吸收声。指定声音的吸收度作为材料的属性 MU （边界导纳或吸收系数）。流体-结构交互作用标记（FSI）这表示模型的流体-结构部分之间的接触表面。步骤五：指定载荷步选项对谐波声分析有如下选项：表5-3对谐波声分析的载荷步选项选项命令菜单途径动力学选项强迫频率范围HARFR工Main Menu>Solution>-Load StepOpts-Time/Frequenc> Freq & Substeps一般选项谐波数NSUBSMain Menu>Solution>-Load Step Opts-Time/ Frequenc>Freq & Substeps阶跃载荷或斜 坡载荷KBCMain Menu>Solution>-Load Step Opts-Time/ Frequenc> Freq & Substeps输出控制打印输出控制OUTPFMain Menu>Solution>Output Ctrls>Solu Printout数据库和结果 文件输出OUTREMain Menu>Solution>Output Ctrls>DB/ Results File结果的外推ERESXMain Menu>Solution>Output Ctrls>Integration Pt动力学选项在此范畴内唯一有效的选项即强迫力的频率范围，在谐波分析中是必须定义的（循环/时间）。在此范围内，然后指定要求解的解的个数。（见“常规选项”）常规选项可以请求任意数目的待计算谐波求解数（通过 NSUBS命令或其相应的GUI 途径）。解（或子步）将在当前指定的频率范围内平均分布。例如，如果在 30 到40Hz内指定10个解，程序就会计算在31, 32, 33,，39和40Hz处的响应。 在低端频率处的响应将不进行计算。载荷可以是阶梯式的或斜坡式的（通过 KB0T令或其GUI途径）。缺省地， 它们是斜坡式的；即载荷值在每一个子步逐渐地增加。 通过阶跃载荷，可以得到 在频率范围内的所有相同的载荷值。输出控制如果在打印输出文件中要包括任何结果数据用 OUTP威其相应的GUI途径（）。OUTRES其GUI途径控制结果文件（）中的数据。ERES双其GUI途径允 许通过拷贝到节点来查看单元积分点的结果而不是通过外推得到（缺省）。注意-缺省地程序只将每个载荷步的最后一个子步写到结果文件中。如果想 要将所有子步（在所有频率处的解）写入结果文件，要保证用OUTRESFREQ选项设置为ALL （或1）。步骤六：备份数据库用Utility Menu>File>Save as或SAV瑜令保存数据库的备份拷贝为一个命名的文件中。这样做使用户在求解失败时能恢复模型。（要恢复模型，重新进入 ANSY和执行 RESUME令或选择 Utility Menu>File>Resume.）步骤七：施加另外的载荷步（可选择）如果要施加另外的载荷步，重复步骤 5和6。步骤8:完成求解进行分析的求解然后结束命令：SOLVEGUI: Main Menu>Solution>-Solve-Current LS命令： FINISHGUI: Main Menu>Finish查看结果ANSY邻序将谐波声分析的结果写到结构的结果文件中。结果包括下列数据, 所有的结果在每一个的强迫频率处都按简谐变化：? 初始数据- 节点压力- 节点位移? 得出数据- 节点和单元压力梯度- 节点和单元应力- 单元力- 节点反力可用POST恢POST2而看此信息。流体 -结构相互作用流体和结构在网格界面处的相互作用引起声压施加给结构一个强迫力， 并且 结构运动产生一个有效的“流体载荷”。有限元的控制矩阵方程变为：（ 1）（ 2）R是一个耦合矩阵代表与流体-结构界面（FSI）上的节点相联系的有效表面面积。耦合矩阵R 也考虑进了组成接触表面的每一对重合的流体和结构单元面的法线矢量方向。ANSYS§序使用的法线矢量的正方向定义为由流体网格以外朝向结构的方向。结构和流体载荷量都是定义在流体- 结构的界面处并为节点自由度的未知函数。将未知的载荷量放在方程的左边并且将两方程合并为一个方程，产生如下结果：3）方程（3）表明着流体- 结构界面处的节点包括位移和压力自由度。应用示例以下两个问题为声单元应用的例子。 第一个例子说明声吸收单元在构造远场问题模型的应用，并通过 GUI和批处理工具进行示范。第二个例子是近场问题， 说明了在一个封闭空间中驻波预测的运用。例 1：流体 - 结构耦合声场分析（命令方法）在此例中，将确定一个浸没在水中延伸到无限远的环孔的第一阶弯曲模态（卵形模态）。使用谐波分析法通过对34Hz和38Hz频率之间进行频率扫描。环 的中心到无限单元的距离至少应等于（D/2） +,这里D是环的外径。l=c/f是压 力波的主波长。 用作为环的外径， 声速为1460， 主频率为 36000/2p ， 给出 （ D/2）+=+（1460）（2p）/36000= 。距离相应的结果频率为。/BATCH,LIST/VERIFY,EV129-1S/PREP7/TITLE,AMA,EV129-1S,FLUID129,HARMONIC ANALYSISET,1,PLANE42! structural elementET,2,FLUID29! acoustic fluid element with ux & uyet,3,129! acoustic infinite line elementr,3,0,0ET,4,FLUID29,1,0! acoustic fluid element without ux & uy! material propertiesMP,EX,1,MP,DENS,1,7929MP,NUXY,1,0MP,DENS,2,1030MP,SONC,2,1460! create inner and outer quarter circlesCYL4,0,0,0,90CYL4,0,0,0,90! select, assign attribute to and mesh area 1ASEL,S,AREA,1AATT,1,1,1,0LESIZE,1,16,1LESIZE,3,16,1LESIZE,2,1,1LESIZE,4,1,1MSHKEY,1MSHAPE,0,2D! mapped quad meshAMESH,1! select, assign attribute to and mesh area 2ASEL,S,AREA,2AATT,2,1,2,0LESIZE,5,16,1LESIZE,7,16,1LESIZE,6,5LESIZE,8,5MSHKEY,0MSHAPE,0,2D! mapped quad meshAMESH,2reflect quarter circle into semicircle about x-axisnsym,x,1000,all! offset node number by 1000esym,1000,all! reflect semicircle into full circle about y-axis nsym,y,2000,all! offset node number by 2000esym,2000,allNUMMRG,ALL! merge all quantities! modify outer 2 layers of el29 into type 4esel,s,type,1nsle,sesln,s,0nsle,sesel,invensle,semodif,all,type,4esel,allnsel,all! define el129 line elementcsys,1nsel,s,loc,x,type,3real,3mat,2esurfesel,allnsel,all! flag interface as fluid-structure interfacensel,s,loc,x,esel,s,type,2sf,all,fsi,1nsel,allesel,allFINISH! enter solution module/SOLUANTYPE,harmic! select harmonic analysishropt,fullf,19,fx,1000.f,1019,fx,-1000.harfrq,34.,38.nsubst,100kbc,1SOLVEFINISH/post26plcplx,0nsol,2,1,u,x,d1uxstoreconjug,3,2prod,4,2,3sqrt,5,4*get,uxmx,vari,5,extrem,tmax/COM/COMExpected Result:/COM/COM The following "uxmx" should equal Hz./COM*status,uxmxfinish例 2 ：室内声场分析(命令方法)此例的问题示例了利用FLUID30单元预测一个典型盒子表示房间的声学驻波 模式。 吸声材料放在盒子的下表面， 圆筒表面的振动结构放在盒子的右上角。 此 问题将确定结构在激励频率为80Hz 时振动的声压波模式。/batch,list/com, Harmonic Analysis - Room Acoustics/PREP7/TITLE,Room Acoustic AnalysisANTYPE,HARMET,1,30! Acoustic elements in contact with walls and vibrating surfaceET,2,30,1Acoustic elements in interior (not in contact with walls)! Set parameters for mesh generationXDIV=29! Number of divisions along x-axisYDIV=19! Number of divisions along y-axisZDIV=1! Number of divisions along z-axisCDIV=2! Number of divisions along radius! Dimensions of the roomLEN=27HGT=20RAD=! Mesh generationK,1K,2,LENK,3,LEN,HGTK,4,HGTK,5,-1K,6,LEN,-1K,7,LEN,HGT,-1K,8,HGT,-1L,1,5,1L,2,6,1L,3,7,1L,4,8,1CIRC,3,RAD,7,2,90,2ADRAG,5,6, 3PIO4=ATAN(1)LENC=COS(PIO4)LENC=LENC*RADHGTC=HGT-LENCLENC=LEN-LENCK,15,HGTCK,16,HGTC,-1K,17,LENCK,18,LENC,-1L,1 ,17,XDIVL,10,15,XDIVL,11,4 ,XDIVL,17,10,YDIVL,15, 1,YDIVL,2 ,9 ,YDIVL,5 ,18,XDIVL,13,16,XDIVL,14,8 ,XDIVL,18,13,YDIVL,16,5 ,YDIVL,6 ,12,YDIVESIZE,CDIVV,1 ,17,10,15, 5,18,13,16V,15,10,11,4 ,16,13,14,8V,17,2 ,9 ,10,18,6 ,12,13VMESH,ALL! Material propertiesMP,DENS,1,! (Rho) density of air (lb/ft*3)MP,SONC,1,! (C) speed of sound in air (ft/sec)MP,MU,1,! (Beta) absorption coefficient of the walls! ("Beta" should be between 0 - 1)MP,DENS,2,MP,SONC,2,MP,MU,2,! (Beta) absorption coefficient of! the absorbing material on the floor! ("Beta" should be between 0 -1)MP,DENS,3,MP,SONC,3,MP,MU,3,! (Beta) zero absorption coefficient for interior elementsNSEL,S,LOC,Y,NSEL,R,LOC,X,12,15ESLNMAT,2EMODIF,ALL!Elements which have the absorbing material on the floor! Boundary conditionsALLSNSEL,S,LOC,X,NSEL,A,LOC,Y,NSEL,A,LOC,X,LENNSEL,A,LOC,Y,HGTD,ALL,UX,UY,UZ! Constrain all displacements to zero at the wallsLOCAL,11,1,27,20NSEL,A,LOC,X,RADESLNESEL,INVETYPE,2MAT,3EMODIF,ALL! Interior elements are specified as Type=2 & material=3ALLS! Fluid-Structure Interface (FSI)NSEL,S,LOC,X,RAD ! Select interface (FSI) surface nodesESLN! Select elements attached to interface surfaceSF,ALL,FSI! Specify vibrating surface as Fluid-structure interfaceNROTAT,ALLD,ALL,UX,.01! Radial vibration amplitue of Vibrating surfaceD,ALL,UY,UZ! Impedance Surface (IMPD)CSYS,0ALLSNSEL,S,LOC,X,NSEL,A,LOC,X,LENNSEL,A,LOC,Y,NSEL,A,LOC,Y,HGTSF,ALL,IMPD,1! Specify the walls as Impedance Surface flag! to activate absorptionALLSFINISH/SOLU! Excitation Frequency for Harmonic AnalysisHARF,80,80! Frequency of excitation = 80 HzSOLVEFINISH! Plot the Standing Wave Pattern (f = 80 Hz)/POST1/SHOW,ENCL1,GRPH,1SET,1,1! plot the real part of pressure responseEPLOT/EDGE,1/TITLE,-ROOM ACOUSTICS-* REAL PART OF PRESSUREPLNS,PRESSET,1,1,1! plot the imag. part of pressure response/TITLE,-ROOM ACOUSTICS-* IMAG. PART OF PRESSUREPLNS,PRESFINISH