ansys 声学分析

第五章 声学

5.1什么是声学？

声学研究声压波在流体介质中的产生、传播、吸收和反射。声学有如下的应用：

·声纳—声学上雷达的对应物 ·设计音乐厅，希望声压均匀分布。 ·减小机器厂房内的噪音 ·汽车中的噪声消除 ·水下声学

·设计扬声器、音箱、声滤、消音器及其他类似装置。 ·地球物理探测

5.1.1声场分析的类型

只有在ANSYS/Multiphysics 和 ANSYS/Mechanical中能进行声场分析，通常包括对流体介质及其周围结构的建模。典型感兴趣的是不同频率的声波在流体中的压力分布、压力梯度、粒子速度、声压级及声波的散射、衍射、传输、辐射、衰减和散射。耦合的声场分析将考虑流体-结构的相互作用。非耦合的声场分析模型只考虑流体而忽略任何流体-结构的相互作用。

ANSYS程序假定流体是可压的，但只允许压力与平均压力相比有较小的变化。而且，流体假定为非流动并且无粘的（即粘性不引起耗散作用）。假定平均密度和平均压力不变，压力求解偏离平均压力而不是绝对压力。

5.2求解声学问题

通过执行一个谐波响应分析可以解决许多声学问题。分析计算流体-结构界面上的谐波载荷（正弦变化）引起流体中的压力分布。通过指定载荷的频率范围，可以观察到在不同的频率时压力的分布。可以执行模态和瞬态的声学分析。（参见《ANSYS Structural Analysis Guide》中关于这种分析更详细的叙述。）

谐波声场分析的过程包括以下三个主要步骤： ·建立模型。

·施加边界条件和载荷并获得求解。

·查看结果。

5.3建立模型

在此步骤中，用户指定工作名称和分析标题，然后用PREP7前处理器定义单元类型，单元实常数，材料属性和模型几何尺寸。这些任务与多数分析相同，在《ANSYS Basic Analysis Guide》中有叙述。

5.3.1谐波声场分析准则

对一个谐波声场分析，考虑以下几点：

单元类型—ANSYS声场分析指定了四种单元类型：对二维和三维模型的流体部分分别使用Fluid29和Fluid30单元，Fluid129和Fluid130与FLUID29和FLUID30单元一起使用，用来构造包围Fluid29和Fluid30单元的无限外壳。利用这些单元类型可以构造流体部分的模型，然后利用相应的结构单元（PLANE42、SOLID45等）构造固体模型。只有Fluid29和Fluid30单元才能与结构单元相接触（在结构的内部或外部）；Fluid129和Fluid130单元只能与Fluid29和Fluid30单元相接触，而不能直接与结构单元接触。

5.3.1.1 FLUID29与FLUID30单元

对与固体相接触的声单元，要确保使用KEYOPT(2)=0，缺省的设置允许流体-结构的相互作用。UX，UY，UZ和PRES作为自由度引起单元矩阵的不对称。对所有其它的声单元，设置KEYOPT（2）=1，致使带有PRES自由度的单元矩阵的对称。（见图5-1）对称矩阵需要的内存和计算时间更少，因此只要可能就应该使用它。关于流体-结构的相互作用的详细信息参见《ANSYS, Inc. Theory Reference》。

图5-1二维声模型的例子（流体在结构的内部）

5.3.1.2 FLUID129和FLUID130单元

对无限的吸收压力波的声单元，模拟在FLUID29和FLUID30单元之外无限延伸域的输出效果。FLUID129和FLUID130单元提供了第二级的吸收边界条件，所以输出的压力波到达模型的边界以最小的反射吸收到流体域内。

FLIUD129单元用来建立二维流体区域的边界和诸如线单元。FLIUD130单元用来建立三维流体区域的边界和诸如平面表面单元。

材料属性—声单元需要密度（DENS）和声速（SONC）作为材料属性（FLUID129和FLUID130只需要SONC）。如果在流体-结构界面存在声的吸收，利用标记MU来指定边界导纳β（吸收系数）。值β通常由实验来测定。对结构单元，指定杨氏模量（EX）、密度（DENS）和泊松比（PRXY或NUXY）。

实常数—当用FLUID129和FLUID130单元时，里面的有限元网格边界必须是圆形的（二维及轴对称）或球形的（三维），而且圆形或球形边界的有限区域的半径必须指定为实常的RAD。（见图5-2）圆或球的中心也必须用实常数指定：

R,3,RAD,X0,Y0!REAL set 3 for FLUID129 R,3,RAD,X0,Y0,Z0!REAL set 3 for FLUID130

如果二维轴对称情况圆的中心坐标（X0，Y0）或三维情况球的中心坐标（X0，Y0，Z0）不是通过实常数指定的，ANSYS假定中心为总体坐标系的原点。

图5-2应用吸收单元的例子

图5-3浸在水中的圆筒

5.4对模型进行网格划分

下面为使用二维无限的声单元进行网格划分的典型程序。对三维单元的程序与此相同。如果还有结构部分，结构单元必须与FLUID29单元相邻，而不要与无限流体单元（FLUID129）相邻。

这个过程自动在有限区域的边界加入FLUID129单元。这里环形结构用PLANE42结构单元进行网格划分。与PLANE42单元接触的流体单元层用带有UX和UY自由度的FLUID29单元划分网格同时打开流体-结构接触面。流体的外层用没有UX和UY自由度的FLUID29单元建模。在X0=Y0=0处放置FLUID129单元的半径为0.31242（参见§5.4.2节）。可用下列方法定义FLUID129单元：

命令：ESURF

GUI: Main Menu>Preprocessor>Create>Elements>Inf Acoustic

5.4.1步骤一：内部流体区域的网格划分

用FLUID29单元对圆形或球形边界（PLANE42）所包围的内部流体区域划分网格。

图5-4 对流体区域划分网格

5.4.2步骤二：生成无限的声单元

按以下步骤：

1．选择圆形或球形边界上的节点： 命令：NSEL

GUI: Utility Menu>Select>Entities

2．指定FLUID129单元作为与FLUID29单元相联系的无限流体单元。 命令：TYPE, REAL

GUI: Main Menu>Preprocessor>Attributes>Default Attribs Main Menu>Preprocessor>Real Constants

无限单元对低频及高频激励都执行得很好。数值实验已确定吸收单元远离结构或振动源区域以外大约0.2λ放置能产生准确的结果。这里λ=c/f压力波的主波长。c是流体中的声速（SONC），f是压力波的主频。例如，对浸没在水中的圆盘或球壳的直径D，封闭边界的半径RAD至少应为D/2+0.2λ。而且，对一般的声分析，网格必须足够的细致以能分辨最小的主频。 3．在边界生成吸收单元（FLUID129）。 命令：ESURF

GUI: Main Menu>Preprocessor>Create>Elements>Inf Acoustic

图5-5 在边界加入吸收单元

5.4.3步骤三：指定流体-结构接触面 指定流体-结构接触面： 1．选择界面上的所有节点。 命令：NSEL

GUI: Utility Menu>Select Entities 2．选择附属于这些节点上的流体单元。 命令：ESEL

GUI: Utility Menu>Select>Entities

3．指定所选择的节点作为流体-结构接触面上的节点。 命令：SF

GUI: Main Menu>Preprocessor>Loads>Apply>Fluid-Struct>On Nodes 注意—在对分析进行求解前确保重选所有的节点。

图5-6 指定流体-结构界面

5.5施加载荷并获得求解

在这个步骤里，用户定义分析类型和选项，施加载荷，指定载荷步选项，并开始有限元求解。下面几节解释如何做这些工作。

5.5.1步骤一：进入SOLUTION求解器

通过选择GUI途径Main Menu>Solution或执行/SOLU命令进入SOLUTION求解器。

5.5.2步骤二：定义分析类型

利用GUI途径或命令集，定义分析类型和分析选项。 用下列方法定义分析类型： 命令：ANTYPE,HARMIC

GUI: Main Menu>Solution>New Analysis

必须选择新的分析NEW Analysis因为重启动对谐波响应分析无效。如果需要施加另外的谐波载荷，每次作一个新的分析（或用《ANSYS Basic Analysis Guide》中叙述的“部分求解”程序）。

5.5.3步骤三：定义分析选项

用下列方法指定求解方法： 命令：HROPT

GUI: Main Menu>Solution>Analysis Options

尽管全量、减缩或模态叠加方法都是可选的方法，选择全量方法因为它自己就可以处理非对称矩阵。

用下列方法定义求解列表格式： 命令：HROUT

GUI: Main Menu>Solution>Analysis Options

这个选项决定ANSYS如何在打印输出（Jobname.OUT）中对谐波自由度结果进行列表。

用下列方法指定方程求解器： 命令：EQSLV

GUI: Main Menu>Solution>Analysis Options

可以选择波前求解器（缺省），Jocobi共轭梯度（JCG）求解器，或不完全的Cholesky共轭梯度（ICCG）求解器。对大多数模型推荐使用JCG求解器。

5.5.4步骤四：在模型上施加载荷

由谐波分析的定义可知，假定任何施加的载荷随时间是简谐（正弦）变化的。在声学分析中完整地指定一个谐波载荷，通常需要两条信息：强迫力的频率和振幅。振幅是载荷的最大值，如表5-2所示用命令指定振幅。强迫力的频率是谐波载荷的频率（循环/时间）。可用HARFRQ命令或GUI途径Main Menu>

Solution>Time/Frequenc>Freq & Substeps在随后的载荷步中指定它。参见§5.5.5节步骤5。

表5-1列出了所有对谐波声分析可用的载荷及定义、列表和删除它们的命令。除了惯性载荷，用户既可以在实体模型（关键点、线和面）上也可以在有限元模型（节点和单元）上定义载荷。对实体模型载荷与有限元载荷综合的讨论参见ANSYS基本分析程序指南中加载与求解的有关章节。

表5-1声学分析中可用的载荷

载荷类型 种命令族 菜单途径 类 Main Menu>Preprocessor>-Loads-> 位移 (UX,UY,UZ), 约束 压力 (PRES) D Apply>Displacement or Potential Main Menu>Solution>Apply> Displacement or Potential Main 力(FX,FX,FZ), 力矩 (MX,MY,MZ), 流体载荷 Menu>Preprocessor>-Loads-> Apply>Force/Moment 力 F Main Menu>Solution>-Loads-Apply> Force/Moment 压力 (PRES) 表阻抗 (IMPD) 面载流体-结构相互作荷 用标记(FSI) ACEL, 重力 旋转等 惯性载荷 OMEGA, Apply>load type DOMEGA, CGLOC, CGOMEGA, DCGOM, Main Menu>Solution>Apply>load IRLF type SF Main Menu>Preprocessor>-Loads-> Apply>load type Main Menu>Solution>Apply>load type Main Menu>Preprocessor>-Loads-> 在一个分析中，用户可以施加、移去、操作或列表载荷

5.5.4.1用GUI施加载荷

通过一系列的层叠菜单可以得到除列表（见下面）以外的所有载荷操作。从求解菜单中选择操作（应用等），载荷类型（位移、力等），然后是施加载荷的目标（关键点等）。

例如，按以下GUI途径，给一条线施加位移载荷： Main Menu>Solution>Apply>Displacement>On Lines

用GUI途径对载荷进行列表： Utility Menu>List>Loads>load type

5.5.4.2利用命令施加载荷

表5-2列出在声分析中所有可用来加载的命令。 表5-2在声分析中施加载荷的命令

载荷类型 实体模型或有限元 实体模型 位移 压力 图元 施加 关键点 DK DL DA D FK F SFL SFA SF 删除 列表 操作 施加设置 DKDELE DKLIST DTRAN DLDELE DLLIST DTRAN DADELE DALIST DTRAN 实体模型 线 实体模型 面 有限元 节点 力 实体模型 力矩 关键点 DDELE DLIST DSCALE DCUM FKDELE FKLIST FTRAN 有限元 节点 实体模型 线 实体模型 面 FDELE FLIST FSCALE FCUM SFLDELE SFLLIST SFTRAN SFGRAD SFADELE SFALIST SFTRAN SFGRAD SFCUM SFDELE SFLIST SFSCALE SFGRAD SFBEAM SFCUM 压力 阻抗 流体结构相互作用标记 有限元 节点 有限元 单元 SFE SFEDELE SFELIST SFSCALE SFFUN SFGRAD ACEL CGLOC CGOMGA 惯性 DCGOM OMEGA DOMEGA IRLF 5.5.4.3载荷类型

位移（UX，UY，UZ）和压力（PRES）

这些是DOF（自由度）约束。例如，在刚性的流体-结构接触面上指定零位移。也可以指定非零位移，但应记住它们被假定为简谐的。在自由流面边界通常指定零压力（此处流体不封闭，诸如一个开口）。

力（FX，FY，FZ）和力矩（MX，MY，MZ） 通常在模型的固体部分指定载荷以激励流体。

当得知一个节点处的声载荷，可以通过指定FLOW力标记来指定流动载荷： 命令：F,,FLOW

GUI: Main Menu>Solution>Apply>-Structural-Force/Moment>On Nodes

压力（PRES）

可以在固体部分指定表面载荷来代替力和力矩。

阻抗（IMPD）

实际上显示不是载荷而是表面吸收声。指定声音的吸收度作为材料的属性MU（边界导纳或吸收系数）。

流体-结构交互作用标记（FSI）

这表示模型的流体-结构部分之间的接触表面。

5.5.5步骤五：指定载荷步选项

对谐波声分析有如下选项：

表5-3对谐波声分析的载荷步选项

选项 命令 菜单途径 动力学选项 强迫频率范围 HARFRQ Main Menu>Solution>-Load Step Opts-Time/Frequenc> Freq & Substeps 一般选项 谐波数 NSUBST Main Menu>Solution>-Load Step Opts-Time/

Frequenc>Freq & Substeps 阶跃载荷或斜坡载荷 KBC Main Menu>Solution>-Load Step Opts-Time/ Frequenc> Freq & Substeps 输出控制 打印输出控制 OUTPR Main Menu>Solution>Output Ctrls>Solu Printout 数据库和结果文件输出 OUTRES Main Menu>Solution>Output Ctrls>DB/ Results File 结果的外推 ERESX Main Menu>Solution>Output Ctrls>Integration Pt 5.5.5.1动力学选项

在此范畴内唯一有效的选项即强迫力的频率范围，在谐波分析中是必须定义

的（循环/时间）。在此范围内，然后指定要求解的解的个数。（见“常规选项”）

5.5.5.2常规选项

可以请求任意数目的待计算谐波求解数（通过NSUBST命令或其相应的GUI途径）。解（或子步）将在当前指定的频率范围内平均分布。例如，如果在30到40Hz内指定10个解，程序就会计算在31，32，33，?，39和40Hz处的响应。在低端频率处的响应将不进行计算。

载荷可以是阶梯式的或斜坡式的（通过KBC命令或其GUI途径）。缺省地，它们是斜坡式的；即载荷值在每一个子步逐渐地增加。通过阶跃载荷，可以得到在频率范围内的所有相同的载荷值。

5.5.5.3输出控制

如果在打印输出文件中要包括任何结果数据用OUTPR或其相应的GUI途径（Jobname.OUT）。OUTRES和其GUI途径控制结果文件（Jobname.RST）中的数据。ERESX及其GUI途径允许通过拷贝到节点来查看单元积分点的结果而不是通过外推得到（缺省）。

注意-缺省地程序只将每个载荷步的最后一个子步写到结果文件中。如果想要将所有子步（在所有频率处的解）写入结果文件，要保证用OUTRES将FREQ选项设置为ALL（或1）。

5.5.6步骤六：备份数据库

用Utility Menu>File>Save as 或SAVE命令保存数据库的备份拷贝为一个命名的文件中。这样做使用户在求解失败时能恢复模型。（要恢复模型，重新进入ANSYS并执行RESUME命令或选择Utility Menu>File>Resume.）

5.5.7步骤七：施加另外的载荷步（可选择）

如果要施加另外的载荷步，重复步骤5和6。

5.5.8步骤8：完成求解

进行分析的求解然后结束。 命令：SOLVE

GUI: Main Menu>Solution>-Solve-Current LS 命令：FINISH

GUI: Main Menu>Finish

5.6查看结果

ANSYS程序将谐波声分析的结果写到结构的结果文件Jobname.RST中。结果包括下列数据，所有的结果在每一个的强迫频率处都按简谐变化：

?初始数据 - 节点压力 - 节点位移 ?得出数据

- 节点和单元压力梯度 - 节点和单元应力 - 单元力 - 节点反力

可用POST1或POST26查看此信息。

5.7流体-结构相互作用

流体和结构在网格界面处的相互作用引起声压施加给结构一个强迫力，并且结构运动产生一个有效的“流体载荷”。有限元的控制矩阵方程变为：

（1）

（2）

[R]是一个耦合矩阵代表与流体-结构界面（FSI）上的节点相联系的有效表面面积。耦合矩阵[R]也考虑进了组成接触表面的每一对重合的流体和结构单元面的法线矢量方向。ANSYS程序使用的法线矢量的正方向定义为由流体网格以外朝向结构的方向。结构和流体载荷量都是定义在流体-结构的界面处并为节点自由度的未知函数。将未知的载荷量放在方程的左边并且将两方程合并为一个方程，产生如下结果：

（3）

方程（3）表明着流体-结构界面处的节点包括位移和压力自由度。

5.8应用示例

以下两个问题为声单元应用的例子。第一个例子说明声吸收单元在构造远场问题模型的应用，并通过GUI和批处理工具进行示范。第二个例子是近场问题，说明了在一个封闭空间中驻波预测的运用。

5.9例1：流体-结构耦合声场分析（命令方法）

在此例中，将确定一个浸没在水中延伸到无限远的环孔的第一阶弯曲模态（卵形模态）。使用谐波分析法通过对34Hz和38Hz频率之间进行频率扫描。环的中心到无限单元的距离至少应等于（D/2）+0.2l，这里D是环的外径。l=c/f是压力波的主波长。用0.26035作为环的外径，声速为1460，主频率为36000/2p，给出（D/2）+0.2l=0.26035+(0.2)(1460)(2p)/36000=0.311。0.31242距离相应的结果频率为35.240Hz。

/BATCH,LIST /VERIFY,EV129-1S /PREP7

/TITLE,AMA,EV129-1S,FLUID129,HARMONIC ANALYSIS ET,1,PLANE42! structural element

ET,2,FLUID29! acoustic fluid element with ux & uy et,3,129! acoustic infinite line element r,3,0.31242,0,0

ET,4,FLUID29,,1,0! acoustic fluid element without ux & uy

! material properties MP,EX,1,2.068e11 MP,DENS,1,7929 MP,NUXY,1,0 MP,DENS,2,1030 MP,SONC,2,1460

! create inner and outer quarter circles CYL4,0,0,0.254,0,0.26035,90 CYL4,0,0,0.26035,0,0.31242,90

! select, assign attribute to and mesh area 1 ASEL,S,AREA,,1 AATT,1,1,1,0 LESIZE,1,,,16,1 LESIZE,3,,,16,1 LESIZE,2,,,1,1 LESIZE,4,,,1,1 MSHKEY,1

MSHAPE,0,2D! mapped quad mesh AMESH,1

! select, assign attribute to and mesh area 2 ASEL,S,AREA,,2 AATT,2,1,2,0 LESIZE,5,,,16,1 LESIZE,7,,,16,1

LESIZE,6,,,5 LESIZE,8,,,5 MSHKEY,0

MSHAPE,0,2D! mapped quad mesh AMESH,2

! reflect quarter circle into semicircle about x-axis nsym,x,1000,all! offset node number by 1000 esym,,1000,all

! reflect semicircle into full circle about y-axis nsym,y,2000,all! offset node number by 2000 esym,,2000,all

NUMMRG,ALL! merge all quantities

! modify outer 2 layers of el29 into type 4 esel,s,type,,1 nsle,s esln,s,0 nsle,s esel,inve nsle,s

emodif,all,type,4 esel,all nsel,all

! define el129 line element csys,1

nsel,s,loc,x,0.31242 type,3 real,3 mat,2 esurf esel,all nsel,all

! flag interface as fluid-structure interface nsel,s,loc,x,0.26035 esel,s,type,,2 sf,all,fsi,1 nsel,all esel,all FINISH

! enter solution module /SOLU

ANTYPE,harmic! select harmonic analysis hropt,full f,19,fx,1000. f,1019,fx,-1000. harfrq,34.,38. nsubst,100 kbc,1 SOLVE

FINISH /post26 plcplx,0

nsol,2,1,u,x,d1ux store conjug,3,2 prod,4,2,3 sqrt,5,4

*get,uxmx,vari,5,extrem,tmax

/COM ------------------------------------------------------------- /COMExpected Result: /COM

/COM The following \

/COM ------------------------------------------------------------- *status,uxmx finish

5.10 例2：室内声场分析（命令方法）

此例的问题示例了利用FLUID30单元预测一个典型盒子表示房间的声学驻波模式。吸声材料放在盒子的下表面，圆筒表面的振动结构放在盒子的右上角。此问题将确定结构在激励频率为80Hz时振动的声压波模式。

/batch,list

/com, Harmonic Analysis - Room Acoustics /PREP7

/TITLE,Room Acoustic Analysis

ANTYPE,HARM

ET,1,30! Acoustic elements in contact with walls and vibrating surface

ET,2,30,,1Acoustic elements in interior (not in contact with walls) ! Set parameters for mesh generation XDIV=29! Number of divisions along x-axis YDIV=19! Number of divisions along y-axis ZDIV=1! Number of divisions along z-axis CDIV=2! Number of divisions along radius ! Dimensions of the room LEN=27 HGT=20 RAD=0.9

! Mesh generation K,1 K,2,LEN K,3,LEN,HGT K,4,,HGT K,5,,,-1 K,6,LEN,,-1 K,7,LEN,HGT,-1 K,8,,HGT,-1 L,1,5,1 L,2,6,1

L,3,7,1 L,4,8,1

CIRC,3,RAD,7,2,90,2 ADRAG,5,6,,,,, 3 PIO4=ATAN(1) LENC=COS(PIO4) LENC=LENC*RAD HGTC=HGT-LENC LENC=LEN-LENC K,15,,HGTC K,16,,HGTC,-1 K,17,LENC K,18,LENC,,-1 L,1 ,17,XDIV L,10,15,XDIV L,11,4 ,XDIV L,17,10,YDIV L,15, 1,YDIV L,2 ,9 ,YDIV L,5 ,18,XDIV L,13,16,XDIV L,14,8 ,XDIV L,18,13,YDIV L,16,5 ,YDIV

L,6 ,12,YDIV ESIZE,,CDIV

V,1 ,17,10,15, 5,18,13,16 V,15,10,11,4 ,16,13,14,8 V,17,2 ,9 ,10,18,6 ,12,13 VMESH,ALL

! Material properties

MP,DENS,1,2.35E-3! (Rho) density of air (lb/ft**3) MP,SONC,1,1100.0! (C) speed of sound in air (ft/sec) MP,MU,1,0.04! (Beta) absorption coefficient of the walls ! (\MP,DENS,2,2.35E-3 MP,SONC,2,1100.0

MP,MU,2,0.70! (Beta) absorption coefficient of ! the absorbing material on the floor ! (\MP,DENS,3,2.35E-3 MP,SONC,3,1100.0

MP,MU,3,0.0! (Beta) zero absorption coefficient for interior elements

NSEL,S,LOC,Y,0.0 NSEL,R,LOC,X,12,15 ESLN MAT,2

EMODIF,ALL! Elements which have the absorbing material on the floor ! Boundary conditions ALLS

NSEL,S,LOC,X,0.0 NSEL,A,LOC,Y,0.0 NSEL,A,LOC,X,LEN NSEL,A,LOC,Y,HGT

D,ALL,UX,,,,,UY,UZ! Constrain all displacements to zero at the walls LOCAL,11,1,27,20 NSEL,A,LOC,X,RAD ESLN ESEL,INVE TYPE,2 MAT,3

EMODIF,ALL! Interior elements are specified as Type=2 & material=3 ALLS

! Fluid-Structure Interface (FSI)

NSEL,S,LOC,X,RAD ! Select interface (FSI) surface nodes ESLN! Select elements attached to interface surface

SF,ALL,FSI! Specify vibrating surface as Fluid-structure interface NROTAT,ALL

D,ALL,UX,.01! Radial vibration amplitue of Vibrating surface D,ALL,UY,,,,,UZ

! Impedance Surface (IMPD)

CSYS,0 ALLS

NSEL,S,LOC,X,0.0 NSEL,A,LOC,X,LEN NSEL,A,LOC,Y,0.0 NSEL,A,LOC,Y,HGT

SF,ALL,IMPD,1! Specify the walls as Impedance Surface flag ! to activate absorption ALLS FINISH /SOLU

! Excitation Frequency for Harmonic Analysis HARF,80,80! Frequency of excitation = 80 Hz SOLVE FINISH

! Plot the Standing Wave Pattern (f = 80 Hz) /POST1

/SHOW,ENCL1,GRPH,1

SET,1,1! plot the real part of pressure response EPLOT /EDGE,1

/TITLE,-ROOM ACOUSTICS-* REAL PART OF PRESSURE PLNS,PRES

SET,1,1,,1! plot the imag. part of pressure response

/TITLE,-ROOM ACOUSTICS-* IMAG. PART OF PRESSURE PLNS,PRES FINISH

本文来源：https://www.bwwdw.com/article/7ohp.html