ansys 声学分析Word文件下载.docx

ansys 声学分析Word文件下载.docx

《ansys 声学分析Word文件下载.docx》由会员分享，可在线阅读，更多相关《ansys 声学分析Word文件下载.docx（18页珍藏版）》请在冰豆网上搜索。

关于流体-结构的相互作用的详细信息参见《ANSYS,Inc.TheoryReference》。

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

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!

REALset3forFLUID129

R,3,RAD,X0,Y0,Z0!

REALset3forFLUID130

如果二维轴对称情况圆的中心坐标（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（参见§

命令：

ESURF

GUI:

MainMenu>

Preprocessor>

Create>

Elements>

InfAcoustic

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

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

按以下步骤：

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

NSEL

UtilityMenu>

Select>

Entities

2．指定FLUID129单元作为与FLUID29单元相联系的无限流体单元。

TYPE,REAL

Attributes>

DefaultAttribs

MainMenu>

RealConstants

无限单元对低频及高频激励都执行得很好。

数值实验已确定吸收单元远离结构或振动源区域以外大约0.2λ放置能产生准确的结果。

这里λ=c/f压力波的主波长。

c是流体中的声速（SONC），f是压力波的主频。

例如，对浸没在水中的圆盘或球壳的直径D，封闭边界的半径RAD至少应为D/2+0.2λ。

而且，对一般的声分析，网格必须足够的细致以能分辨最小的主频。

3．在边界生成吸收单元（FLUID129）。

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

指定流体-结构接触面：

1．选择界面上的所有节点。

SelectEntities

2．选择附属于这些节点上的流体单元。

ESEL

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

SF

Loads>

Apply>

Fluid-Struct>

OnNodes

注意—在对分析进行求解前确保重选所有的节点。

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

5.5施加载荷并获得求解

在这个步骤里，用户定义分析类型和选项，施加载荷，指定载荷步选项，并开始有限元求解。

下面几节解释如何做这些工作。

通过选择GUI途径MainMenu>

Solution或执行/SOLU命令进入SOLUTION求解器。

利用GUI途径或命令集，定义分析类型和分析选项。

用下列方法定义分析类型：

ANTYPE,HARMIC

Solution>

NewAnalysis

必须选择新的分析NEWAnalysis因为重启动对谐波响应分析无效。

如果需要施加另外的谐波载荷，每次作一个新的分析（或用《ANSYSBasicAnalysisGuide》中叙述的“部分求解”程序）。

用下列方法指定求解方法：

HROPT

AnalysisOptions

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

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

HROUT

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

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

EQSLV

可以选择波前求解器（缺省），Jocobi共轭梯度（JCG）求解器，或不完全的Cholesky共轭梯度（ICCG）求解器。

对大多数模型推荐使用JCG求解器。

由谐波分析的定义可知，假定任何施加的载荷随时间是简谐（正弦）变化的。

在声学分析中完整地指定一个谐波载荷，通常需要两条信息：

强迫力的频率和振幅。

振幅是载荷的最大值，如表5-2所示用命令指定振幅。

强迫力的频率是谐波载荷的频率（循环/时间）。

可用HARFRQ命令或GUI途径MainMenu>

Solution>

Time/Frequenc>

Freq&

Substeps在随后的载荷步中指定它。

参见§

表5-1列出了所有对谐波声分析可用的载荷及定义、列表和删除它们的命令。

除了惯性载荷，用户既可以在实体模型（关键点、线和面）上也可以在有限元模型（节点和单元）上定义载荷。

对实体模型载荷与有限元载荷综合的讨论参见ANSYS基本分析程序指南中加载与求解的有关章节。

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

载荷类型

种类

命令族

菜单途径

位移（UX,UY,UZ）,

压力（PRES）

约束

D

-Loads->

DisplacementorPotential

力（FX,FX,FZ）,

力矩（MX,MY,MZ）,

流体载荷

力

F

Force/Moment

-Loads-Apply>

阻抗（IMPD）

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

表面载荷

loadtype

重力

旋转等

惯性载荷

ACEL,

OMEGA,

DOMEGA,CGLOC,CGOMEGA,DCGOM,IRLF

在一个分析中，用户可以施加、移去、操作或列表载荷

通过一系列的层叠菜单可以得到除列表（见下面）以外的所有载荷操作。

从求解菜单中选择操作（应用等），载荷类型（位移、力等），然后是施加载荷的目标（关键点等）。

例如，按以下GUI途径，给一条线施加位移载荷：

Displacement>

OnLines

用GUI途径对载荷进行列表：

UtilityMenu>

List>

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

表5-2在声分析中施加载荷的命令

实体模型或有限元

图元

施加

删除

列表

操作

施加设置

位移

压力

实体模型

关键点

DK

DKDELE

DKLIST

DTRAN

线

DL

DLDELE

DLLIST

面

DA

DADELE

DALIST

有限元

节点

DDELE

DLIST

DSCALE

DCUM

力矩

FK

FKDELE

FKLIST

FTRAN

FDELE

FLIST

FSCALE

FCUM

阻抗

流体结构相互作用标记

SFL

SFLDELE

SFLLIST

SFTRAN

SFGRAD

SFA

SFADELE

SFALIST

SFDELE

SFLIST

SFSCALE

SFCUM

单元

SFE

SFEDELE

SFELIST

SFBEAM

SFFUN

惯性

ACELCGLOC

CGOMGA

DCGOM

OMEGA

DOMEGA

IRLF

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

这些是DOF（自由度）约束。

例如，在刚性的流体-结构接触面上指定零位移。

也可以指定非零位移，但应记住它们被假定为简谐的。

在自由流面边界通常指定零压力（此处流体不封闭，诸如一个开口）。

力（FX，FY，FZ）和力矩（MX，MY，MZ）

通常在模型的固体部分指定载荷以激励流体。

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

F,,FLOW

-Structural-Force/Moment>

压力（PRES）

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

阻抗（IMPD）

实际上显示不是载荷而是表面吸收声。

指定声音的吸收度作为材料的属性MU（边界导纳或吸收系数）。

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

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

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

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

选项

命令

动力学选项

强迫频率范围

HARFRQ

-LoadStepOpts-Time/Frequenc>

Freq&

Substeps

一般选项

谐波数

NSUBST

-LoadStepOpts-Time/Frequenc>

阶跃载荷或斜坡载荷

KBC

输出控制

打印输出控制

OUTPR

OutputCtrls>

SoluPrintout

数据库和结果文件输出

OUTRES

DB/ResultsFile

结果的外推

ERESX

IntegrationPt

在此范畴内唯一有效的选项即强迫力的频率范围，在谐波分析中是必须定义的（循环/时间）。

在此范围内，然后指定要求解的解的个数。

（见“常规选项”）

可以请求任意数目的待计算谐波求解数（通过NSUBST命令或其相应的GUI途径）。

解（或子步）将在当前指定的频率范围内平均分布。

例如，如果在30到40Hz内指定10个解，程序就会计算在31，32，33，…，39和40Hz处的响应。

在低端频率处的响应将不进行计算。

载荷可以是阶梯式的或斜坡式的（通过KBC命令或其GUI途径）。

缺省地，它们是斜坡式的；

即载荷值在每一个子步逐渐地增加。

通过阶跃载荷，可以得到在频率范围内的所有相同的载荷值。

如果在打印输出文件中要包括任何结果数据用OUTPR或其相应的GUI途径（Jobname.OUT）。

OUTRES和其GUI途径控制结果文件（Jobname.RST）中的数据。

ERESX及其GUI途径允许通过拷贝到节点来查看单元积分点的结果而不是通过外推得到（缺省）。

注意-缺省地程序只将每个载荷步的最后一个子步写到结果文件中。

如果想要将所有子步（在所有频率处的解）写入结果文件，要保证用OUTRES将FREQ选项设置为ALL（或1）。

用UtilityMenu>

File>

Saveas或SAVE命令保存数据库的备份拷贝为一个命名的文件中。

这样做使用户在求解失败时能恢复模型。

（要恢复模型，重新进入ANSYS并执行RESUME命令或选择UtilityMenu>

Resume.）

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

进行分析的求解然后结束。

SOLVE

-Solve-CurrentLS

FINISH

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,HARMONICANALYSIS

ET,1,PLANE42!

structuralelement

ET,2,FLUID29!

　acousticfluidelementwithux&

uy

et,3,129!

acousticinfinitelineelement

r,3,0.31242,0,0

ET,4,FLUID29,,1,0!

　acousticfluidelementwithoutux&

!

materialproperties

MP,EX,1,2.068e11

MP,DENS,1,7929

MP,NUXY,1,0

MP,DENS,2,1030

MP,SONC,2,1460

createinnerandouterquartercircles

CYL4,0,0,0.254,0,0.26035,90

CYL4,0,0,0.26035,0,0.31242,90

select,　assignattributetoandmesharea1

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!

mappedquadmesh

AMESH,1

select,　assignattributetoandmesharea2

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

AMESH,2

　reflectquartercircleintosemicircleaboutx-axis

nsym,x,1000,all!

offsetnodenumberby1000

esym,,1000,all

　reflectsemicircleintofullcircleabouty-axis

nsym,y,2000,all!

offsetnodenumberby2000

esym,,2000,all

NUMMRG,ALL!

mergeallquantities

modifyouter　2layersofel29intotype4

esel,s,type,,1

nsle,s

esln,s,0

esel,inve

emodif,all,type,4

esel,all

nsel,all

defineel129lineelement

csys,1

nsel,s,loc,x,0.31242

type,3

real,3

mat,2

esurf

flaginterfaceasfluid-structureinterface

nsel,s,loc,x,0.26035

esel,s,type,,2

sf,all,fsi,1

entersolutionmodule

/SOLU

ANTYPE,harmic!

selectharmonicanalysis

hropt,full

f,19,fx,1000.

f,1019,fx,-1000.

harfrq,34.,38.

nsubst,100

kbc,1

/post26

plcplx,0

nsol,2,1,u,x,d1ux

store

conjug,3,2

p