用ANSYS软件分析压电换能器入门

用ANSYS软件分析压电换能器入门
A：分析过程基本步骤
一：问题描述（草稿纸上完成）
1：画出换能器几何模型，包括尺寸
2：选定材料
3：查材料手册确定材料参数
二：建立模型
1：根据对称性确定待建模型的维数
2：根据画出的几何模型确定关键点坐标，给关键点编好号码
3：建立一个文件夹用于当前分析
4：启动ANSYS软件，指定路径到建立的文件夹，
5：定义单元类型
压电换能器分析使用的单元类型：
solid5：8个节点3D六面体耦合场单元（也可缩减为三角柱形单元或四面体单元）。

无实常数。

plane13：4个节点2D四边形耦合场单元（也可缩减为三角形单元）。

无实常数。

solid98：10个节点3D四面体耦合场单元。

无实常数。

Fluid30：8个节点3D六面体声学流体单元（也可缩减为三角柱形单元或四面体单元）。

应用于近场水和远场水。

实常数为参考声压，可缺省。

Fluid130：4个节点面无穷吸收水声学流体单元（也可缩减为三角形面单元）。

实常数：半径，球心X，Y，Z坐标值。

6：定义材料参数
对一般均匀各向同性材料要给出材料密度，杨氏模量，泊松系数。

（静态分析不用密度）
对压电材料：
一般使用的压电方程：e型压电方程，因此输入的常数为
注意！一般顺序为：XX，YY，ZZ，YZ，XZ，XY。

在ANSYS中为XX，YY，ZZ，XY，YZ，XZ。

因此，前两矩后三行和后三列要做相应变化。

7：建立关键点
8：把关键点连成线
9：把线段围成面
10：通过适当的方法生成体
11：指定单元类型和材料参数
12：划分线段
13：划分体单元
14：坐标转换，（转换到柱坐标系下）
15：节点转换
三：加载约束条件
1：加载边界约束条件
2：电极上加电压
四：求解
1：模态分析
2：谐响应分析
五：查看结果
1：查看模态分析结果，计算导纳。

2：各模态的动态演示
3：查看谐响应分析结果，计算导纳、发射与接收响应。

六：生成命令流文件
1：给程序分块，添加适当的注释
2：把相应参数具体值改成变量，同时给变量赋值
B：空气中建模过程
一：问题描述
弯曲式换能器实体模型为轴对称结构。

结构尺寸如图所示。

弯曲式换能器由压电陶瓷，‘金属壳’粘结而成，粘结材料为环氧树脂。

由于粘结层为一薄层，建模的时候可以忽略。

本例中‘金属壳’材料采用铝，压电材料选择压电陶瓷PZT5。

二：建立模型
1．指定工作空间，设定分析作业名和标题
（1）首先建立一个文件夹用于当前分析，如D:\ANSYSfiles ，路径当中不允许有中文出现。

（2）启动ANSYS9.0软件，界面如下图，界面包括实用菜单，命令输入窗口，ANSYS工具栏，主菜单栏，图形窗口，图形调整工具栏，标准工具栏。

（3）指定工作空间：实用菜单－》File－》Change Directory… ,出现对话框，选择刚才建立的目录D:\ANSYSfiles，以后操作生成的文件都在此文件夹下。

（4）设定标题：实用菜单－》File－》Change Title…，出现对话框，输入‘daocheqi Analysis ’作为标题。

点击ok。

（5）保存db（数据库文件）文件：File－》Save as ,出现对话框，保存文件名为title.db。

（6）定义单位：在ANSYS主界面输入窗口输入“/UNITS，Label”其中，Label 指定单位制标签有：
USER：用户自定义单位系统。

SI：国际单位，即MKS单位。

CGS：CGS单位系统，即cm,g,s,c 等。

MPA：MPA单位系统，即mm,Kg,s,c 等。

BFT ：用英尺的单位系统，即ft,slug,s, F
等。

BIN ：用英寸的单位系统，即in,lbm,s, F
等。

一般采用国际单位：即 然后按Enter 即可。

如果不改变单位制，以后默认。

2．定义单元类型
金属壳采用Solid45 单元，压电陶瓷采用Solid5 单元。

Solid45为结构单元里面的三维块模型单元。

Solid5为三维耦合场单元，因为压电陶瓷涉及到电－力耦合问题，因此采用此单元类型。

(1）
主菜单栏－》Preprocessor －》Element Type －》Add/Edit/Delete , 出现对话框，点击Add …键，出现对话框，点击框图左边Solid , 选择右边框图中Brick 8node 45，点击Apply 键，选定了一种Solid45单元。

继续点击框图左边Coupled Field , 选择右边框图中Scalar Brick 5，选定了Solid5单元，点击ok 键。

从对话框中可以看出我们已经定义了solid45和solid5两种单元类型。

点击close 键关闭窗口。

（2）保存文件：File- Save as ,出现对话框，保存文件名为 element.db
3．定义材料参数
查阅材料手册我们可以得到如下材料参数。

铝材料参数：
密度3/2790m kg =ρ，杨氏模量2
10/1015.7m N E ⨯=，泊松比34.0=σ。

压电陶瓷PZT5的材料参数：（ANSYS 输入参数） 密度3
7750/kg m ρ=
刚度矩阵（也可输入柔顺矩阵） 压电应力常数矩阵 介电常数矩阵 也可输入绝对值
（1） 输入黄铜材料参数：主菜单栏－》Preprocessor->>Material Props －》Material Models ,
出现对话框，点击左边Material Model Number 1 ,点击右边Favorites －》Linear Static —》Density ，出现密度对话框，输入密度2790，点击ok 。

同样继续点击Linear Isotropic （线性各项同性）输入泊松比和杨氏模量。

输入7.15e10和0.34，点击ok 。

（2）输入PZT5材料参数：定义完一种材料参数之后要定义第二中材料，选择对话框‘Define Material Model Behavior’ 的菜单Material －》New Model …，出现对话框，材料添2（默认值也是2），点击ok 。

定义第二中材料，首先定义PZT5 的密度。

如图点击左边Material Model Number 2 ,点击右边Structural —》Density ，输入7750，点击ok 。

定义刚度矩阵，点击右边Structural —》Linear －》Elastic －》Anisotropic,输入参数 D11=12.1e10，D12=7.54e10，D13=7.52e10，D22=12.1e10，D23=7.52e10，D33=11.1e10，D44=2.26e10，D55=2.11e10，D66=2.11e10
定义介电常数矩阵，点击右边Electromagnetics —》Relative Permittivity －》Orthotropic,
输入参数PERX=8.11e-9 ,PERY=8.11e-9 ,PZRZ=7.35e-9
定义压电矩阵：点击右边Piezoelectrics—》Piezoelectric matrix,输入如下图
（3）保存文件matiral.db。

4．建立关键点
（1）首先在草稿纸上画出几何模型，指定好关键点的编号以及结构尺寸。

主菜单栏－》Preprocessor－》Modeling －》Create－》Keypoints－》In Actice CS，出现对话框如下图，根据指定好的编号即结构尺寸输入坐标。

（注意，由于最后要转换到柱坐标系下分析，所以y即为旋转方向）
输入坐标为：
1 (0, 0, 0)
2 (0, 0, 0.5e-3)
3 (0, 0, 1.0e-3)
4 (10e-3, 0,0)
5 (10e-3, 0, 0.5e-3)
6 (10e-3, 0, 1.0e-3)
7 (17e-3, 0, 0)
8 (17e-3, 0, 0.5e-3,)
9 (20e-3, 0, 0)
10 (20e-3, 0, 0.5e-3)
11 (20e-3, 0, 10.5e-3)
12 (20e-3, 0, 13.5e-3)
13 (17e-3, 0, 10.5e-3)
14 (17e-3, 0, 13.5e-3)
15 (0, 0, 10.5e-3)
16 (0, 0, 13.5e-3)
输入完之后我们可以看到图形窗口中我们所定义的关键点
（2）保存文件：keypoints.db
5. 把关键点连成线
主菜单栏－》Preprocessor－》Modeling －》Create－》lines－》Straight Line
通过鼠标选取两个关键点就可以生成一条线，如此操作，把所有的线都连好，点击ok。

保存文件：line.db
6.把线段围成面
主菜单栏－》Preprocessor－》Modeling －》Create－》Areas－》By Lines
通过鼠标顺时针（逆时针）方向选择线段围成一个面积，如此操作，把所有的线都围成面。

保存文件：area.db
7. 旋转生成体积，并粘接
（1）主菜单栏－》Preprocessor－》Modeling －》Operate－》Extrude－》Areas－》About Axis, 出现一个对话框，点击Pick All , 表示全部面都要旋转。

之后又出现对话框，要求选择旋转轴，通过鼠标选择Z轴上的两个点，比如第1，3两个点，点击ok，出现对话框，要求说明旋转的度数，这里选择30度模型，输入30。

点击ok，。

生成如下图的体。

（2）粘接各部分体积：主菜单栏－》Preprocessor－》Modeling －》Operate－》Boolean－》Glue－》Volumes，
出现对话框，要求选择要粘接的体积，这里选择Pick All 表示所有体积都粘在一起。

（3）保存文件：volume.db
8．指定各部分单元类型以及材料属性
（1）指定金属壳的单元及材料属性：主菜单栏－》Preprocessor－》Meshing－》Mesh Attributes －》Picked Volumes
用鼠标点中金属壳的体积单元
出现对话框，材料参数选择1，单元类型选择1 SOLID45
（2）指定压电材料的单元及材料属性：主菜单栏－》Preprocessor－》Meshing－》Mesh Attributes－》Picked Volumes
用鼠标点中压电材料的体积单元
出现对话框，材料参数选择2，单元类型选择2 SOLID5
（3）点击ok ，保存文件material.db
9 .划分线段
为了网格单元数可以控制，我们采用手动划分网格，而不采用自动网格划分。

（1）首先对‘线段’划分，由于体积图看不清一些线的位置，先显示线段图，
实用菜单－》Plot－》Lines
(2) 主菜单－》Preprocessor－》Meshing－》Size Cntrls－》ManualSize－》Lines－》Picked Lines
选择这些线段，（比如压电片连同金属壳在压电片半径内划分为4单元），出现对话框，划
分单元数填4，表示这些线段都划分为4个单元，点击ok，
采用相同的方法划分其他线段的单元数目。

（为了采用6面体划分体积单元，注意几个划分10格的线段划分时要求单元数目一致，当然这些线段也可以划分成10格，但要一致）
10.划分体积单元
主菜单－》Preprocessor－》Meshing－》Mesh－》Volumes－》Mapped－》4 to 6 sided 我们采用6面体单元划分体，这样我们划分的网格大小就比较‘整齐’
出现对话框，点击Pick All，表示所有体积都按6面体单元划分。

（4）保存文件：meshvolume.db
11．转换坐标系
把坐标系转换到柱坐标系下，有利于加载约束条件。

实用菜单－》WorkPlane－》Change Active CS to －》Global Cylindrical
保存文件：changecys.db
12．节点转换
把划分网格后生成的节点转换到柱坐标系下，使得后续的求解过程计算速度更快。

主菜单栏－》Preprocessor－》Modeling －》Create－》Nodes－》Rotate Node CS－》To Active CS ，出现对话框，点击Pick All
保存文件：changenodes.db
三：加载约束条件
1.加载边界约束条件
（1）旋转方向固定
主菜单栏－》Preprocessor－》Loads－》Define Loads－》Apply－》Structural－》Displacement －》On Areas
用鼠标选定如下的几个面，点击ok
出现对话框，点击UY，表示在Y方向，在柱坐标系下也即旋转方向固定，点击ok
(2) 旋转轴固定
主菜单栏－》Preprocessor－》Loads－》Define Loads－》Apply－》Structural－》Displacement －》On Lines
选择轴线，点击ok。

出现如下对话框，选择UX，表示固定其X方面，即在径向上不移动。

保存文件：loaddisp.db
2 .加载电压
给压电陶瓷加上电压，我们给定陶瓷下底面为0V ，陶瓷上底面为1V 。

主菜单栏－》Preprocessor－》Loads－》Define Loads－》Apply－》Electric－》Boundary－》Voltage－》On Areas
选择陶瓷下底面，如下
出现对话框，要求输入电压值，电压值填0V
用同样的方法选择陶瓷上底面，加1V的电压。

保存文件：loadele.db
四求解
1 .模态分析
（1）选择分析类型：主菜单－》Solution－》Analysis Type－》New Analysis
点击选择模态分析选项‘Modal’，点击ok
(2)分析参数设定：
主菜单－》Solution－》Analysis Type－》Analysis Options, 设定求解模态数目为8
点击ok，出现对话框，频率范围可以不设定，直接点击ok
（3）进入求解：主菜单－》Solution－》Solve－》Current LS ，表示采用当前约束条件进行求解
出现下面两个提示，分别点击ok和yes。

程序进行求解，求解完成之后出现求解完毕提示，表示求解完成。

保存文件:solve.db
五查看结果
1 查看模态结果
主菜单－》General Postproc－》Results Summary，我们可以看出求出的前8阶模态分别为：0Hz，4268.9Hz，19476Hz，21172Hz，51261Hz，58121Hz，58445Hz ，73133Hz。

2.动态演示模态结果
(1）读取模态结果读取第1阶模态结果：主菜单－》General Postproc－》Read Results－》First Set
实际上，由于模型未在Z方向加位移限制条件，故第阶1模态频率为零，只作平动。

（2）动态演示第2阶模态的振动形态
读取模态结果读取第2阶模态结果：主菜单－》General Postproc－》Read Results－》Next Set ,实用菜单－》PlotCtrls－》Animate－》Mode Shape…
点击之后出现对话框对动画进行参数设置，我们采用默认参数，
点击ok将会出现第2阶模态的振动形态动画。

也可选择Main Menu>General Postproc>Plot Results>Deformed Shape 命令，弹出对话框
选择[Def+undef edge] ,单击[OK]得到变形图
六生成命令流文件
以上步骤都是在图形操作下完成的，但由于采用‘图形操作’来建模效率太低，而且不易于改变参数（结构参数，材料参数和网格参数等），因此我们有必要采用文本编程的方法来建模和求解。

ANSYS软件提供了这方法，但是要完全掌握ANSYS编程语言相对比较困难，而采用‘图形操作’建模过程中，软件会按照图形操作过程生成相应的文本命令，称之为命令流文件（.lgw）。

因此我们在‘图形操作’建模之后生成命令流文件，然后对命令流文件进行编辑整
理，分块，‘具体参数变量化’(如轴对称模型只分析圆柱的一定角度，电极给定电压，划分线网格数)等等操作之后，再由软件读入该文件就可以得到我们所要求解的东西。

1 生成.lgw文件
实用菜单－》File－》Write DB log file…，输入docheqi.lgw文件名。

（用蓝色字表示注释，红色标注提请注意）
2 编辑命令流文件
注意：所有以“！”开始的语句为注释语句，并不参加执行，可以去掉。

/BATCH
/TITLE,normal cymbal
/PREP7
ET,1,SOLID45！定义第一种单元
ET,2,SOLID5！定义第二种单元
!*
!*
! SAVE, element,db,D:\ANSYSfiles\docheqi Analisis\
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@material1金属壳材料
mat1dens=2790 ！铝密度
mat1yang=7.15e10 ！铝杨氏模量
mat1person=0.34 ！铝泊松比
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@material3压电体材料PZT5
mat2dens=7750 ！密度
c11=12.1e10 ！弹性常数矩阵
c12=7.54e10
c13=7.52e10
c22=12.1e10
c23=7.52e10
c33=11.1e10
c44=2.26e10
c55=2.11e10
c66=2.11e10
e13=-5.4 ！压电常数矩阵
e23=-5.4
e33=15.8
e52=12.3
e61=12.3
fx=0.811e-8 ！介电常数矩阵
fy=0.811e-8
fz=0.735e-8
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@structural几何结构参数
r1=10e-3 ！陶瓷片半径
r2=17e-3 ！铝壳内半径
r3=20e-3 ！铝壳外半径
tj=0.5e-3 ！双迭片中铝片厚度
h=10e-3 ！铝壳内腔深度
ty=0.5e-3 ！压电片厚度
td=3e-3 ！铝壳底厚度
theta=30 ！分析换能器的旋转角
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ meshline网格分线控制tcbj=4
lvbj=4
bhow=2
hudu=10
dibj=10
dihow=4
shen=8
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@volut control电压控制voltvalue=1
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,mat1yang ！铝杨氏模量MPDATA,PRXY,1,,mat1person ！铝泊松比MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,mat1dens ！铝密度
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,2,,mat2dens！压电材料密度
TB,ANEL,2,1,21,0
TBTEMP,0
TBDATA,,c11,c12,c13,,, ！压电材料弹性系数TBDATA,,c22,c23,,,,c33
TBDATA,,,,,c44,,
TBDATA,,c55,,c66,,,
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,PERX,2,,fx！压电材料介电常数MPDATA,PERY,2,,fy
MPDATA,PERZ,2,,fz
TB,PIEZ,2,,,0
TBMODIF,1,1,
TBMODIF,1,2,
TBMODIF,1,3,e13 ！压电材料压电常数TBMODIF,2,1,
TBMODIF,2,2,
TBMODIF,2,3,e23
TBMODIF,3,1,
TBMODIF,3,2,
TBMODIF,3,3,e33
TBMODIF,4,1,
TBMODIF,4,2,
TBMODIF,4,3,
TBMODIF,5,1,
TBMODIF,5,2,e52
TBMODIF,5,3,
TBMODIF,6,1,e61
TBMODIF,6,2,
TBMODIF,6,3,
! SAVE, matiral constant,db,D:\ANSYSFILES\docheqi Analisis\
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
K,1,0,0,0, ！定义关键点
K,2,0,0,tj,
K,3,0,0,tj+ty,
K,4,r1,0,0,
K,5,r1,0,tj,
K,6,r1,0,tj+ty,
K,7,r2,0,0,
K,8,r2,0,tj,
K,9,r3,0,0,
K,10,r3,0,tj,
K,11,r3,0,tj+h,
K,12,r3,0,tj+td+h,
K,13,r2,0,tj+h,
K,14,r2,0,tj+td+h,
K,15,0,0,tj+h,
K,16,0,0,tj+td+h,
!SAVE,keypoint,db,D:\ANSYSFILES\docheqi Analisis\
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ LSTR, 1, 2 ！关键点连线
LSTR, 2, 3
LSTR, 4, 5
LSTR, 5, 6
LSTR, 7, 8
LSTR, 9, 10
LSTR, 10, 11
LSTR, 11, 12
LSTR, 8, 13
LSTR, 13, 14
LSTR, 15, 16
!!@@@@@@@@@@@@@@@@@@@@@@@2
LSTR, 1, 4
LSTR, 4, 7
LSTR, 2, 5
LSTR, 5, 8
LSTR, 3, 6
LSTR, 7, 9
LSTR, 8, 10
LSTR, 15, 13
LSTR, 13, 11
LSTR, 16, 14
LSTR, 14, 12
! SAVE, line,db,D:\ANSYSFILES\docheqi Analisis\
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ FLST,2,4,4 ！线围成面
FITEM,2,1
FITEM,2,12
FITEM,2,3
FITEM,2,14
AL,P51X
FLST,2,4,4
FITEM,2,2
FITEM,2,14
FITEM,2,4
FITEM,2,16
AL,P51X
FLST,2,4,4
FITEM,2,13
FITEM,2,5
FITEM,2,15
FITEM,2,3
AL,P51X
FLST,2,4,4
FITEM,2,5
FITEM,2,17
FITEM,2,6
FITEM,2,18
AL,P51X
FLST,2,4,4
FITEM,2,9
FITEM,2,18
FITEM,2,7
FITEM,2,20
AL,P51X
FLST,2,4,4
FITEM,2,10
FITEM,2,20
FITEM,2,8
FITEM,2,22
AL,P51X
FLST,2,4,4
FITEM,2,11
FITEM,2,19
FITEM,2,10
FITEM,2,21
AL,P51X
! SAVE, area,db,D:\ANSYSFILES\docheqi Analisis\
FLST,2,7,5,ORDE,2 ！面旋转成体
FITEM,2,1
FITEM,2,-7
FLST,8,2,3
FITEM,8,1
FITEM,8,16
VROTAT,P51X, , , , , ,P51X, , theta, , ！旋转角度
! SAVE, volume,db,D:\ANSYSFILES\docheqi Analisis\
FLST,2,7,6,ORDE,2 !粘结体各部分
FITEM,2,1
FITEM,2,-7
VGLUE,P51X
! SAVE, glue,db,D:\ANSYSFILES\docheqi Analisis\
FLST,5,6,6,ORDE,3
FITEM,5,1
FITEM,5,3
FITEM,5,-7
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT, 1, , 1, 0 ！指定金属壳材料为第一种单元，第一种材料属性CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
CM,_Y,VOLU
VSEL, , , , 2
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT, 2, , 2, 0！指定压电材料为第二种单元，第二种材料属性CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
! SAVE, matiral,db,D:\ANSYSFILES\docheqi Analisis\ FLST,5,6,4,ORDE,6
FITEM,5,12
FITEM,5,14
FITEM,5,16
FITEM,5,23
FITEM,5,25
FITEM,5,28
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,tcbj, , , , ,1 ！划分线段网格FLST,5,4,4,ORDE,4
FITEM,5,13
FITEM,5,15
FITEM,5,31
FITEM,5,33
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,lvbj, , , , ,1 ！划分线段网格FLST,5,8,4,ORDE,8
FITEM,5,17
FITEM,5,-18
FITEM,5,20
FITEM,5,22
FITEM,5,36
FITEM,5,38
FITEM,5,42
FITEM,5,47
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,bhow, , , , ,1 ！划分线段网格FLST,5,11,4,ORDE,11
FITEM,5,26
FITEM,5,-27
FITEM,5,30
FITEM,5,34
FITEM,5,-35
FITEM,5,39
FITEM,5,-40
FITEM,5,44
FITEM,5,-45
FITEM,5,49
FITEM,5,-50
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,hudu, , , , ,1 ！划分线段网格
FLST,5,4,4,ORDE,4
FITEM,5,19
FITEM,5,21
FITEM,5,51
FITEM,5,-52
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,dibj, , , , ,1 ！划分线段网格
FLST,5,5,4,ORDE,5
FITEM,5,8
FITEM,5,10
FITEM,5,-11
FITEM,5,46
FITEM,5,48
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,dihow, , , , ,1 ！划分线段网格FLST,5,4,4,ORDE,4
FITEM,5,7
FITEM,5,9
FITEM,5,41
FITEM,5,43
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,shen, , , , ,1 ！划分线段网格
! SAVE, meshline,db,D:\ANSYSFILES\docheqi Analisis\
FLST,5,7,6,ORDE,2 ！划分体网格
FITEM,5,1
FITEM,5,-7
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
MSHAPE,0,3d
MSHKEY,1
VMESH,_Y1
MSHKEY,0
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
! SAVE, meshvolume,db,D:\ANSYSFILES\docheqi Analisis\ CSYS,1！转换坐标系到柱坐标系下
! SAVE, changecys,db,D:\ANSYSFILES\docheqi Analisis\ FLST,2,1697,1,ORDE,2 ！转换节点到柱坐标系下FITEM,2,1
FITEM,2,-1697
NROTAT,P51X
! SAVE, changenodes,db,D:\ANSYSFILES\docheqi Analisis\ !!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
NSEL,S,LOC,Y,0 ！边界约束条件
DSYM,SYMM,Y,1
NSEL,S,LOC,Y,30
DSYM,SYMM,Y,1
NSEL,S,LOC,X,0
DSYM,SYMM,X,1
! SAVE, loaddisp,db,D:\ANSYSFILES\docheqi Analisis\
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
NSEL,S,LOC,Z,tj+ty ！给压电材料表面加电压NSEL,R,LOC,X,0,r1
CP,1,VOLT,ALL
*GET,N1,NODE,,NUM,MIN
D,N1,VOLT,voltvalue ！给定一电极电压值
NSEL,S,LOC,Z,tj
NSEL,R,LOC,X,0,r1
D,ALL,VOLT,0 ！给定另一电极电压值
NSEL,ALL
! SAVE, loadele,db,D:\ANSYSFILES\docheqi Analisis\ FINISH
/SOLU ！求解
ANTYPE,2 ！进行模态分析
MSAVE,0
MODOPT,LANB,8！计算前8个模态
EQSLV,SPAR
MXPAND,0, , ,0
LUMPM,0
PSTRES,0
MODOPT,LANB,8,0,0, ,OFF
SOLVE
! SAVE, solve,db,D:\ANSYSFILES\docheqi Analisis\
C：水中建模过程（兰色字部分操作方法见后面附录）水中建模过程与空气中建模过程基本相同，由于换能器置于无限介质水中，ANSYS设置了近场水、远场水和边界水，其中，近场水只需设置在换能器近处周围，将换能器包围，远场水必须满足远场条件，边界水处于远场外边界，为全吸收层。

例如：双叠片换能器在水中的谐响应分析命令流文件如下：
/BATCH
/TITLE,head
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@material1金属材料
mat1dens=2790
mat1yang=7.15e10
mat1person=0.34
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@material3压电体材料PZT4
mat2dens=7550
c11=13.9e10
c12=7.78e10
c13=7.43e10
c22=13.9e10
c23=7.43e10
c33=11.5e10
c44=3.06e10
c55=2.56e10
c66=2.56e10
e13=-5.2
e23=-5.2
e33=15.1
e52=12.7
e61=12.7
fx=0.328e-8
fy=0.328e-8
fz=0.562e-8
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@material3近场水
mat3dens=1000 ！水密度
sonc3=1480 ！水中声速
mu3=0 ！吸收系数
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@material4远场水
mat4dens=1000
sonc4=1480
mu4=0
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@structura5边界水
mat5dens=1000
sonc5=1480
mu5=1
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ 结构参数
h1=0.002
h2=0.010
h3=0.001
h4=0.0005
h5=0.0005
h6=0.007
r1=0.012
r2=0.015
r3=0.0075
r4=0.011
r5=0.025
r6=0.050
h7=r6*(h1+h2+h3+h4-h6)/r2
h8=r6*(h1+h2+h3+h4+r5-h6)/r2
h9=r6*(h1+h2+h3+h4+r5-h6)/r1
h10=r6*(h1+h2+h3+h4+r5-h6)/r4
h11=r6*(h1+h2+h3+h4+r5-h6)/r3
h12=r6*(h1-h6)/r2
h13=-r6*h6/r2
h14=-r6*h6/r1
theta=45 ！绕轴旋转45度
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@电压控制参数voltvalue=1
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@网格划分控制参数tcbj=20
ouco=1
sdpd=3
dbj=20
bhw=4
wqs=18
wqx=12
hudu=20
dihw=6
sdps=12
sdpx=2
tchw=4
ys=20
js=12
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@单元类型
/PREP7
ET,1,SOLID45
ET,2,SOLID5
ET,3,FLUID30
KEYOPT,3,2,0 ！单元类型号，任意值，近场水0远场水1 ET,4,FLUID30
KEYOPT,4,2,1
ET,5,FLUID130
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@实参量
R,1,20e-6, ！近场水实常数，任意号，参考压力
R,2,20e-6, ！远场水实常数，任意号，参考压力
R,3,r6,0,0,0, ！边界水实常数，半径，中心坐标X，Y，Z值
! SAVE, definerealconstants,db,D:\ANSYSFILES\1030WATERNEW\ !!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@金属材料MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,mat1dens
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,mat1yang
MPDATA,PRXY,1,,mat1person
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@压电片MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,2,,mat2dens
TB,ANEL,2,1,21,0
TBTEMP,0
TBDATA,,c11,c12,c13,,,
TBDATA,,c22,c23,,,,c33
TBDATA,,,,,c44,,
TBDATA,,c55,,c66,,,
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,PERX,2,,fx
MPDATA,PERY,2,,fy
MPDATA,PERZ,2,,fz
TB,PIEZ,2,,,0
TBMODIF,1,1,
TBMODIF,1,2,
TBMODIF,1,3,e13
TBMODIF,2,1,
TBMODIF,2,2,
TBMODIF,2,3,e23
TBMODIF,3,1,
TBMODIF,3,2,
TBMODIF,3,3,e33
TBMODIF,4,1,
TBMODIF,4,2,
TBMODIF,4,3,
TBMODIF,5,1,
TBMODIF,5,2,e52
TBMODIF,5,3,
TBMODIF,6,1,e61
TBMODIF,6,2,
TBMODIF,6,3,
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@近场水MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,3,,mat3dens MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,SONC,3,,sonc3 MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,MU,3,,mu3
!!@@@@@@@@@@@@@@@@@@@@@@@@@@远场水MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,4,,mat4dens MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,SONC,4,,sonc4 MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,MU,4,,mu4
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@边界水MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,5,,mat5dens MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,SONC,5,,sonc5 MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,MU,5,,mu5
！！@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
K,1,0,0,-h6,
K,2,0,0,h1-h6,
K,3,0,0,h1+h2-h6,
K,4,0,0,h1+h2+h3-h6,
K,5,0,0,h1+h2+h3+h4-h6,
K,6,0,0,h1+h2+h3+h4+h5-h6,
K,7,r3,0,h1+h2-h6,
K,8,r3,0,h1+h2+h3-h6,
K,9,r3,0,h1+h2+h3+h4-h6,
K,10,r3,0,h1+h2+h3+h4+h5-h6,
K,11,r4,0,h1+h2+h3-h6,
K,12,r4,0,h1+h2+h3+h4-h6,
K,13,r4,0,h1+h2+h3+h4+h5-h6,
K,14,r1,0,-h6,
K,15,r1,0,h1-h6
K,16,r1,0,h1+h2+h3+h4-h6,
K,17,r1,0,h1+h2+h3+h4+h5-h6,
K,18,r2,0,-h6,
K,19,r2,0,h1-h6,
K,20,r2,0,h1+h2+h3+h4-h6,
K,21,r2,0,h1+h2+h3+h4+h5-h6,
K,22,0,0,r5,
K,23,0,0,r6,
K,24,r5,0,0,
K,25,r6,0,0,
K,26,r6,0,h7,
K,27,r6,0,h8,
K,28,r6,0,h9,
K,29,r6,0,h10,
K,30,r6,0,h11,
K,31,0,0,-r5,
K,32,0,0,-r6,
K,33,0,0,0,
K,34,r1,0,0,
K,35,r2,0,0,
K,36,r6,0,h12,
K,37,r6,0,h13,
K,38,r6,0,h14,
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ LSTR, 1, 2
LSTR, 2, 33
LSTR, 3, 4
LSTR, 4, 5
LSTR, 5, 6
LSTR, 7, 8
LSTR, 8, 9
LSTR, 9, 10
LSTR, 11, 12
LSTR, 12, 13
LSTR, 14, 15
LSTR, 15, 34
LSTR, 34, 16
LSTR, 16, 17
LSTR, 18, 19
LSTR, 19, 35
LSTR, 35, 20
LSTR, 20, 21
LSTR, 1, 14
LSTR, 14, 18
LSTR, 2, 15
LSTR, 15, 19
LSTR, 3, 7
LSTR, 4, 8
LSTR, 5, 9
LSTR, 6, 10
LSTR, 8, 11
LSTR, 9, 12
LSTR, 10, 13
LSTR, 12, 16
LSTR, 13, 17
LSTR, 16, 20
LSTR, 17, 21
LSTR, 6, 22
LSTR, 22, 23
LSTR, 33, 34
LSTR, 34, 35
LSTR, 35, 24
LSTR, 24, 25
LARC,22,24,34,r5 ！圆弧线，起始关键点号，终点关键点号，圆心号，半径。

LARC,23,25,34,r6
LSTR, 20, 26
LSTR, 21, 27
LSTR, 17, 28
LSTR, 13, 29
LSTR, 1, 31
LSTR, 31, 32
LSTR, 19, 36
LSTR, 18, 37
LSTR, 14, 38
LARC,24,31,34,r5
LARC,25,32,34,r6
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!分割线，操作如下：
！Main Menu>Preprocessor>Modeling>Operate>Booleans>Partition>Lines FLST,2,12,4,ORDE,4
FITEM,2,41
FITEM,2,-47
FITEM,2,50
FITEM,2,-54
LPTN,P51X
！删除线,操作如下：
!!Main Menu>Preprocessor>Modeling>Delete>Lines Only
FLST,2,8,4,ORDE,4
FITEM,2,69
FITEM,2,-73
FITEM,2,76
FITEM,2,-78
LDELE,P51X
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
FLST,2,4,4
FITEM,2,1
FITEM,2,22
FITEM,2,12
FITEM,2,20
AL,P51X
FLST,2,4,4
FITEM,2,12
FITEM,2,23
FITEM,2,16
FITEM,2,21
AL,P51X
FLST,2,4,4
FITEM,2,13
FITEM,2,38
FITEM,2,17
FITEM,2,23
AL,P51X
FITEM,2,14 FITEM,2,33 FITEM,2,18 FITEM,2,38 AL,P51X FLST,2,4,4 FITEM,2,4 FITEM,2,25 FITEM,2,7 FITEM,2,24 AL,P51X FLST,2,4,4 FITEM,2,5 FITEM,2,26 FITEM,2,8 FITEM,2,25 AL,P51X FLST,2,4,4 FITEM,2,6 FITEM,2,27 FITEM,2,9 FITEM,2,26 AL,P51X FLST,2,4,4 FITEM,2,8 FITEM,2,29 FITEM,2,10 FITEM,2,28 AL,P51X FLST,2,4,4 FITEM,2,9 FITEM,2,30 FITEM,2,11 FITEM,2,29 AL,P51X FLST,2,4,4 FITEM,2,11 FITEM,2,32 FITEM,2,15 FITEM,2,31 AL,P51X FLST,2,4,4 FITEM,2,15
FITEM,2,19 FITEM,2,33 AL,P51X FLST,2,4,4 FITEM,2,35 FITEM,2,63 FITEM,2,55 FITEM,2,27 AL,P51X FLST,2,4,4 FITEM,2,55 FITEM,2,79 FITEM,2,56 FITEM,2,30 AL,P51X FLST,2,4,4 FITEM,2,58 FITEM,2,80 FITEM,2,56 FITEM,2,32 AL,P51X FLST,2,4,4 FITEM,2,58 FITEM,2,81 FITEM,2,62 FITEM,2,34 AL,P51X FLST,2,4,4 FITEM,2,62 FITEM,2,82 FITEM,2,61 FITEM,2,19 AL,P51X FLST,2,4,4 FITEM,2,18 FITEM,2,61 FITEM,2,65 FITEM,2,39 AL,P51X FLST,2,4,4 FITEM,2,17 FITEM,2,39 FITEM,2,66
AL,P51X FLST,2,4,4 FITEM,2,16 FITEM,2,60 FITEM,2,95 FITEM,2,59 AL,P51X FLST,2,4,4 FITEM,2,21 FITEM,2,59 FITEM,2,96 FITEM,2,57 AL,P51X FLST,2,4,4 FITEM,2,20 FITEM,2,57 FITEM,2,74 FITEM,2,48 AL,P51X FLST,2,4,4 FITEM,2,36 FITEM,2,64 FITEM,2,91 FITEM,2,63 AL,P51X FLST,2,4,4 FITEM,2,91 FITEM,2,83 FITEM,2,90 FITEM,2,79 AL,P51X FLST,2,4,4 FITEM,2,90 FITEM,2,84 FITEM,2,89 FITEM,2,80 AL,P51X FLST,2,4,4 FITEM,2,89 FITEM,2,85 FITEM,2,88 FITEM,2,81 AL,P51X
FITEM,2,88
FITEM,2,86
FITEM,2,87
FITEM,2,82
AL,P51X
FLST,2,4,4
FITEM,2,65
FITEM,2,87
FITEM,2,68
FITEM,2,40
AL,P51X
FLST,2,4,4
FITEM,2,66
FITEM,2,40
FITEM,2,67
FITEM,2,92
AL,P51X
FLST,2,4,4
FITEM,2,95
FITEM,2,92
FITEM,2,97
FITEM,2,93
AL,P51X
FLST,2,4,4
FITEM,2,96
FITEM,2,93
FITEM,2,98
FITEM,2,94
AL,P51X
FLST,2,4,4
FITEM,2,74
FITEM,2,94
FITEM,2,75
FITEM,2,49
AL,P51X
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ FLST,2,31,5,ORDE,2
FITEM,2,1
FITEM,2,-31
FLST,8,2,3
FITEM,8,32
FITEM,8,23
VROTAT,P51X, , , , , ,P51X, ,theta, ,
FLST,2,31,6,ORDE,2
FITEM,2,1
FITEM,2,-31
VGLUE,P51X
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ FLST,5,10,6,ORDE,4 ！指定材料类型FITEM,5,1
FITEM,5,-4
FITEM,5,6
FITEM,5,-11
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT, 1, 1, 1, 0 CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
CM,_Y,VOLU
VSEL, , , , 5
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT, 2, 1, 2, 0 CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
FLST,5,10,6,ORDE,2
FITEM,5,12
FITEM,5,-21
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT, 3, 1, 3, 0 CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
FLST,5,10,6,ORDE,2
FITEM,5,22
FITEM,5,-31
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
CMSEL,S,_Y1
VATT, 4, 2, 4, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
FLST,5,10,5,ORDE,10
FITEM,5,101
FITEM,5,104
FITEM,5,107
FITEM,5,110
FITEM,5,113
FITEM,5,116
FITEM,5,120
FITEM,5,122
FITEM,5,125
FITEM,5,128
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y
CMSEL,S,_Y1
AATT, 5, 3, 5, 0,
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
FLST,5,8,4,ORDE,8 ！线网格划分，注意事项见附录FITEM,5,20
FITEM,5,22
FITEM,5,41
FITEM,5,43
FITEM,5,74
FITEM,5,-75
FITEM,5,152
FITEM,5,180
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,dbj, , , , ,1
FLST,5,18,4,ORDE,18
FITEM,5,23
FITEM,5,33
FITEM,5,-34
FITEM,5,38
FITEM,5,46
FITEM,5,50
FITEM,5,54
FITEM,5,73
FITEM,5,81
FITEM,5,85
FITEM,5,96
FITEM,5,98
FITEM,5,122
FITEM,5,134
FITEM,5,149
FITEM,5,162
FITEM,5,177
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,bhw, , , , ,1 FLST,5,8,4,ORDE,8
FITEM,5,31
FITEM,5,-32
FITEM,5,80
FITEM,5,84
FITEM,5,118
FITEM,5,120
FITEM,5,132
FITEM,5,159
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,ouco, , , , ,1 FLST,5,10,4,ORDE,9
FITEM,5,28
FITEM,5,-30
FITEM,5,79
FITEM,5,83
FITEM,5,110
FITEM,5,112
FITEM,5,128
FITEM,5,156
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,sdpd, , , , ,1 FLST,5,12,4,ORDE,10
FITEM,5,24
FITEM,5,-27
FITEM,5,63
FITEM,5,-64
FITEM,5,99
FITEM,5,101
FITEM,5,104
FITEM,5,107
FITEM,5,125
FITEM,5,153
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,tcbj, , , , ,1 FLST,5,8,4,ORDE,8
FITEM,5,14
FITEM,5,18
FITEM,5,65
FITEM,5,68
FITEM,5,72
FITEM,5,76
FITEM,5,140
FITEM,5,168
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,wqs, , , , ,1 FLST,5,13,4,ORDE,13
FITEM,5,6
FITEM,5,9
FITEM,5,11
FITEM,5,15
FITEM,5,19
FITEM,5,86
FITEM,5,108
FITEM,5,116
FITEM,5,119
FITEM,5,123
FITEM,5,137
FITEM,5,165
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,sdps, , , , ,1 FLST,5,8,4,ORDE,8
FITEM,5,13
FITEM,5,17
FITEM,5,53
FITEM,5,66
FITEM,5,-67
FITEM,5,69
FITEM,5,143
FITEM,5,172
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,wqx, , , , ,1 FLST,5,9,4,ORDE,9
FITEM,5,1
FITEM,5,12
FITEM,5,16
FITEM,5,42
FITEM,5,47
FITEM,5,95
FITEM,5,97
FITEM,5,146
FITEM,5,174
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,dihw, , , , ,1 FLST,5,5,4,ORDE,5
FITEM,5,5
FITEM,5,10
FITEM,5,105
FITEM,5,111
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,sdpx, , , , ,1 FLST,5,3,4,ORDE,3
FITEM,5,4
FITEM,5,7
FITEM,5,100
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,tchw, , , , ,1 FLST,5,35,4,ORDE,35
FITEM,5,44
FITEM,5,-45
FITEM,5,51
FITEM,5,-52
FITEM,5,70
FITEM,5,-71
FITEM,5,77
FITEM,5,-78
FITEM,5,102
FITEM,5,-103
FITEM,5,106
FITEM,5,109
FITEM,5,113
FITEM,5,-114
FITEM,5,117
FITEM,5,121
FITEM,5,124
FITEM,5,127
FITEM,5,130
FITEM,5,133
FITEM,5,136
FITEM,5,139
FITEM,5,142
FITEM,5,145
FITEM,5,148
FITEM,5,155
FITEM,5,158
FITEM,5,161
FITEM,5,164
FITEM,5,167
FITEM,5,170
FITEM,5,173
FITEM,5,176
FITEM,5,179
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,hudu, , , , ,1 FLST,5,20,4,ORDE,14
FITEM,5,35
FITEM,5,39
FITEM,5,48
FITEM,5,55
FITEM,5,-62
FITEM,5,126
FITEM,5,129
FITEM,5,131
FITEM,5,135
FITEM,5,138
FITEM,5,141
FITEM,5,144
FITEM,5,147
FITEM,5,150
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,js, , , , ,1 FLST,5,20,4,ORDE,14
FITEM,5,36
FITEM,5,40
FITEM,5,49
FITEM,5,87
FITEM,5,-94
FITEM,5,154
FITEM,5,157
FITEM,5,160
FITEM,5,166
FITEM,5,169
FITEM,5,171
FITEM,5,175
FITEM,5,178
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,ys, , , , ,1
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
FLST,5,31,6,ORDE,2 ！体网格划分FITEM,5,1
FITEM,5,-31
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
MSHAPE,0,3d
MSHKEY,1
VMESH,_Y1
MSHKEY,0
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ FLST,5,10,5,ORDE,10 ！边界水网格划分FITEM,5,101
FITEM,5,104
FITEM,5,107
FITEM,5,110
FITEM,5,113
FITEM,5,116
FITEM,5,120
FITEM,5,122
FITEM,5,125
FITEM,5,128
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
MSHKEY,1
AMESH,_Y1
MSHKEY,0
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
CSYS,1 ！转换柱坐标系下
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
FLST,2,17436,1,ORDE,2 ！转换所有节点到柱坐标下FITEM,2,1
FITEM,2,-17436
NROTAT,P51X
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
NSEL,S,LOC,Y,0 ！加载边界约束条件DSYM,SYMM,Y,1
NSEL,S,LOC,Y,45
DSYM,SYMM,Y,1
NSEL,S,LOC,X,0
DSYM,SYMM,X,1
! SAVE, loaddisp,db,D:\ANSYSFILES\docheqi Analisis\ !!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
NSEL,S,LOC,Z,h1+h2-h6 ！加载电边界条件
NSEL,R,LOC,X,0,r3
CP,1,VOLT,ALL
*GET,N1,NODE,,NUM,MIN
D,N1,VOLT,voltvalue
NSEL,S,LOC,Z,h1+h2+h3-h6
NSEL,R,LOC,X,0,r3
D,ALL,VOLT,0
NSEL,ALL
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
FLST,2,10,5,ORDE,10 ！换能器与流体分界面的耦合FITEM,2,34
FITEM,2,37
FITEM,2,-38
FITEM,2,42
FITEM,2,46
FITEM,2,55
FITEM,2,62
FITEM,2,65
FITEM,2,69
FITEM,2,-70
SFA,P51X, ,FSI。