流固耦合模态分析实例1

1、圆环在水中（模态分析）【ANSYS】默认分类 2009-11-18 16:48:28 阅读31 评论0 字号：大中小 finish /clear /PREP7 !定义单元类型 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 !材料

2、属性 MP,EX,1,2.068e11 MP,DENS,1,7929 MP,NUXY,1,0 MP,DENS,2,1030 MP,SONC,2,1460 ! 创建四分之一模型 CYL4,0,0,0.254,0,0.26035,90 CYL4,0,0,0.26035,0,0.31242,90 ! 选择属性,网格划分 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 ASEL,S,AREA,2

3、 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 ! 关于Y轴镜像 nsym,x,1000,all ! offset node number by 1000 esym,1000,all ! 关于y轴镜像 nsym,y,2000,all ! offset node number by 2000 esym,2000,all NUMMRG,ALL ! merge all quantities esel,s,type,1 nsle

4、,s esln,s,0 选择附在节点上的单元（0表示单元有一个节点被选中即选中单元）nsle,s esel,inve 反选单元nsle,s emodif,all,type,4 esel,all nsel,all ! 指定无限吸收边界 csys,1 nsel,s,loc,x,0.31242 type,3 real,3 mat,2 esurf esel,all nsel,all ! 标识流固交接面 nsel,s,loc,x,0.26035 esel,s,type,2 sf,all,fsi,1 nsel,all esel,all FINISH /solu antype,modal modopt,da

5、mp,10 mxpand,10,yes solve finish 水箱盛水【ANSYS】默认分类 2009-11-18 16:48:38 阅读61 评论0 字号：大中小 length=2 width=3height=2/prep7 et,1,63 et,2,30 !选用FLUID30单元，用于流固耦合问题 r,1,0.01 mp,ex,1,2e11 mp,nuxy,1,0.3 mp,dens,1,7800 mp,dens,2,1000 !定义Acoustics材料来描述流体材料-水 mp,sonc,2,1400 mp,mu,0, ! block,length,width,height esiz

6、e,0.5 mshkey,1 ! type,1 mat,1 real,1 asel,u,loc,y,width amesh,all alls ! type,2 mat,2 vmesh,all fini /solu antype,2 modopt,unsym,10 !非对称模态提取方法处理流固耦合问题 eqslv,front mxpand,10,1 nsel,s,loc,x, nsel,a,loc,x,length nsel,r,loc,y d,all,ux,uy,uz, nsel,s,loc,y,width, d,all,pres,0 alls asel,u,loc,y,width, sfa,

7、all,fsi !定义流固耦合界面 alls solv fini /post1 set,first plnsol,u,sum,2,1 finiANSYS 声学计算算例 水下圆柱壳体的建模与声学分析 使用有限元软件ANSYS进行计算和分析时水下环肋圆柱壳体有限元模型的建立及结构声学分析主要分为以下一些步骤：1.建立壳体的实体模型（包括有圆柱壳体的建立，给圆柱壳体加环肋）； 2.圆柱壳体外部流体介质的生成； 3.对圆柱壳体和流体介质进行有限元4.设置流固耦合单元，并设置外部声场边界条件；5.在求解器中进行振动模态求解和受激励的谐响应求6.求解结果进行后处理分析。，1.建立壳体的实体模型这个步骤主要

8、是在预处理模块（PREP7）中完成首先是根据要建立的实体模型，进行单元的选取和定义这些单元的物理属 性，水下圆柱壳体半径与壳体壁厚的比超过了20，根据ANSYS中单元的使用原则可以选用Shell63号薄壳单元，这种单元的有限元计算原理在前面已经介绍；环肋选用梁单元，ANSYS提供了多种梁单元的结构形式，其中Beam188号梁单元符合作为壳体加强筋及肋骨的使用，所以在水下圆柱壳体环肋选用T的Beam188号梁单元进行建模；而流体介质根据分析中用途的不同要定义两种，一种是流体介质中的单元Fluid30号流体介质单元，一种是流体与结构接触的流固耦合面的单元选用Fluid30号流固耦合单元，在实际建模

9、操作中还需要定义一种用于平面声场的29号单元（在计算中未用到，但在建模中需使用）；共需要定义五种单元。Shell63壳体单元与Beam188梁单元为同一种材料，所以物理属性相同。而Fluid30流体介质单元与Fluid30流固耦合单元物理属性也相同，及在分析中只需要定义两个物理属性即可。具体的使用APDL命令定义为：  /prep7 ！进行预处理模块et，1，30， ！定义1号单元为Fluid30 流固耦合单元 et，2，29 ！定义2号单元为Fluid29平面流体单元et，3，30 ， ，1 ！定义3号单元为Fluid30流体介质单元et，4，63 ！定义4号单元为Shell63壳

10、体单元et，5，188 ！定义5号单元为Beam188梁单元r，4，0.002 ！定义4号单元的厚度为2 mp，dens，4，7800 ！定义4号物理属性包括有密度mp，ex，4，2.1e11  !杨氏模量、 mp，nuxy，4，0.3 ！泊松比mp，sonc，1，1460 ！设置水中声速 mp，dens，1，1000 ！设置流体密度sectype， 1， beam， T， ! 选取T型梁secoffset,，orig ！设置梁的方向secdata，0.04，0.05，0.002，0.02，0，0，0，0，0，0 、  所建立的圆柱壳体的参数：圆柱长为 50 ，半

11、径为 25 ，壳体的壁厚为 2 ，   cyl4，0，0，0.25，5 ！形成圆面 k，9，0，0，0 ! 定义原点 k，10，0，0，0.5 lstr，9，10  ！通过原点作直线 adrag，5，6，7，8，9 ！通过放样形成圆柱wpoff，0，0，0.1 asel，s，2，5 asbw，all， ！移动工作平面与选取的侧面相切 ！重复上面操作，形成四个环肋面 wpoff，0，0，-0.4 ！工作平面回到原点位置上 k，31，0.2，0，0.1 ！定义环肋的方向点 lsel，s，20 ！选择要划分为环肋的线段latt，4，5，5，31，40，1 ！定义线段物理

12、属性 lesize，20，6 ！划分数目 secnum，1 lmesh，20  ！划分线段 将上述的操作完成以后，壳体的模型基本完成,具体结构如图示 图 2-6 环肋圆柱壳体模型图 2.圆柱壳体外部流体圆柱壳体外部的流体介质主要通过设置好的平面流体单元沿指定的线段进行放样，形成立体的流体介质单元。这里取水体为直径5m的球体，在球体的外边界处声压设置为0（全吸收情况），有实例表明12 这样的近似处理对结构耦合振动的计算影响不大。主要的操作及命令为：  type，2 ! 选取第二种单元lsel，s，1，4 ！选取线段 lesize，all，10 ！线段划分数目lesize，al

13、l，6 amesh，1 ！将通过工作平面的面 1 进行划分esize，0.1，0 ！选取沿面放样的网格大小 mat，1 ！定义放样的形成单元的物理属性vdrag，1，9 ！进行放样形成流体介质   3.对实体进行有限元网格划分壳体环肋的有限元划分在壳体建模后就完成，因为在后面划分十分困难；流体介质的划分在建模时同步就完成了；这部分工作现在主要是对壳体进行网格的划分，选取所有的壳体的面，定义好物质属性和物理常数，选取网格划分的方法和精度，完成对壳体的划分。 元模型如图2-7 示：  流固耦合单元和外部声场边界的设置 网格划分好以后，根据分析的要求

14、要在结构与流体之间设置流固耦合单元， 它是将结构与流体联系在一起的单元，是结构到流体的过渡。它的设置主要有以下要求：必须是设置在所有结构与流体接触的单元上；只能设置一层这样的单元。流固耦合单元是将Fluid30单元的Keyopt（2）设置为0即可，为上面定义的第一种单元。 外部声场的边界条件的设置主要是选取流体介质球体的外部所有的面，在这些面上加载并设定声压为0即可 模拟远场条件asel，s，6，9 !选择壳体面asel，a，2，5 asel，a，14，17nsla，s，1  ！选取面上所有节点esln，s，0， ！选取与这些节点相邻 的所有元素esel，r，type，3

15、！在所选的元素中选取所有3号单元emodif，all，type，1 ！强制转换为1号流固耦合单元sf，all，fsi asel，a，51，52   ！选取最外层流体介质面nsla，s，1 ！选取面上所有节点 d，all，pres，0.0 ！将声压值设为0 5.在求解器中进行求解计算依据上面的参数建立号分析的有限元模型，进入 ANSYS 求解运算对模型进行振动模态分析有限元模型共有元素为20284个FLUID30的流体介质单元；320个FLUID30流固耦合单元；272个SHELL63壳体单元；96个BEAM188梁单元。利用/solu ！进入求解运算器 antype,modal ！进入模态分析 modopt,unsym,10,-100, ！分析方法为不对称法和频率范围 mxpand,10,1,1000 ！模态扩展solve ！求解 为了便于分析及进行比较，求解出（1）同样尺寸的加肋圆柱壳体在真空中 的振 动特性（2）同样尺寸的不加肋圆柱壳体在水中的振动特性，通过数据的比较可以分析圆柱壳体在水中振动的一些特性。通过建立相应的有限元模型，再利用ANSYS计算求解，求解环肋圆柱壳体 有限元模型的前十阶模态频率结果为表1： 在有限元壳体模型上加激励对圆柱壳体在受激励作用下的振动特性进行分 析。 其APDL命令流程序为