【达尔整理】ANSYS流固耦合分析实例命令流
- 格式:doc
- 大小:306.00 KB
- 文档页数:53
达尔文档DareDoc分享知识传播快乐ANSYS流固耦合分析实例命令流本资料来源于网络,仅供学习交流2015年10月达尔文档|DareDoc整理目录ANSYS流固耦合例子命令流............................................................................. 错误!未定义书签。
ANSYS流固耦合的方式 (3)一个流固耦合模态分析的例子1 (3)一个流固耦合模态分析的例子2 (4)一个流固耦合建模的例子 (7)一加筋板在水中的模态分析 (8)一圆环在水中的模态分析 (10)接触分析实例---包含初始间隙 (14)耦合小程序 (19)流固耦合练习 (21)一个流固耦合的例子 (22)使用物理环境法进行流固耦合的实例及讲解 (23)针对液面晃动问题,ANSYS/LS-DYNA提供三种方法 (30)1、流固耦合 (30)2、SPH算法 (34)3、ALE(接触算法) (38)脱硫塔于浆液耦合的分析 (42)ANSYS坝-库水流固耦合自振特性的例子 (47)空库时的INP文件 (47)满库时的INP文件 (49)计算结果 (52)ANSYS流固耦合的方式一般说来,ANSYS的流固耦合主要有4种方式:1,sequential这需要用户进行APDL编程进行流固耦合sequentia指的是顺序耦合以采用MpCCI为例,你可以利用ANSYS和一个第三方CFD产品执行流固耦合分析。
在这个方法中,基于网格的平行代码耦合界面(MpCCI) 将ANSYS和CFD程序耦合起来。
即使网格上存在差别,MpCCI也能够实现流固界面的数据转换。
ANSYS CD中包含有MpCCI库和一个相关实例。
关于该方法的详细信息,参见ANSYS Coupled-Field Analysis Guide中的Sequential Couplin2,FSI solver流固耦合的设置过程非常简单,推荐你使用这种方式3,multi-field solver这是FSI solver的扩展,你可以使用它实现流体,结构,热,电磁等的耦合4,直接采用特殊的单元进行直接耦合,耦合计算直接发生在单元刚度矩阵一个流固耦合模态分析的例子1这是一个流固耦合模态分析的典型事例,采用ANSYS/MECHANICAL可以完成。
处理过程中需要注意以下几个方面的问题:1、单元的选择;2、流体材料模式;3、流固耦合关系的定义;4、模态提取方法。
length=2width=3height=2/prep7et,1,63et,2,30 !选用FLUID30单元,用于流固耦合问题r,1,0.01mp,ex,1,2e11mp,nuxy,1,0.3mp,dens,1,7800mp,dens,2,1000 !定义Acoustics材料来描述流体材料-水mp,sonc,2,1400mp,mu,0,!block,,length,,width,,heightesize,0.5mshkey,1!type,1mat,1real,1asel,u,loc,y,widthamesh,allalls!type,2mat,2vmesh,allfini/soluantype,2modopt,unsym,10 !非对称模态提取方法处理流固耦合问题eqslv,frontmxpand,10,,,1nsel,s,loc,x,nsel,a,loc,x,lengthnsel,r,loc,yd,all,,,,,,ux,uy,uz,nsel,s,loc,y,width,d,all,pres,0allsasel,u,loc,y,width,sfa,all,,fsi !定义流固耦合界面allssolvfini/post1set,firstplnsol,u,sum,2,1fini一个流固耦合模态分析的例子2一实例,水箱采用SHELL63单元,水箱中的水采用FLUID30单元,以下即为整个流固耦合模态计算的命令流文件:length=1width=0.6height=0.8/prep7et,1,63et,2,30 !选用FLUID30单元,用于流固耦合问题r,1,0.01mp,ex,1,2e11mp,nuxy,1,0.3mp,dens,1,7800mp,dens,2,1000 !定义Acoustics材料来描述流体材料-水mp,sonc,2,1400mp,mu,2,!block,,length,,width,,heightesize,0.1mshkey,1!type,1mat,1real,1asel,u,loc,y,widthamesh,allalls!type,2mat,2vmesh,allfini/soluantype,2modopt,unsym,10 !非对称模态提取方法处理流固耦合问题eqslv,frontmxpand,10,,,1nsel,s,loc,x,nsel,a,loc,x,lengthnsel,r,loc,yd,all,,,,,,ux,uy,uz,nsel,s,loc,y,width,d,all,pres,0allsasel,u,loc,y,width,sfa,all,,fsi !定义流固耦合界面allssolvfini/post1set,firstplnsol,u,sum,2,1fini我得出的结果是:SET TIME/FREQ LOAD STEP SUBSTEP CUMULATIVE1 29.195 1 1 12 0.0000 1 1 13 33.202 1 2 24 0.0000 1 2 25 37.598 1 3 36 0.0000 1 3 37 44.592 1 4 48 0.0000 1 4 49 47.917 1 5 510 0.0000 1 5 511 88.015 1 6 612 0.0000 1 6 613 88.565 1 7 714 0.0000 1 7 715 97.133 1 8 816 0.0000 1 8 817 109.49 1 9 918 0.0000 1 9 919 109.81 1 10 1020 0.0000 1 10 10为了比较我还计算了水箱不装水的情况,(水箱体模型的建立和约束与上面的相同,只是把流体部分去掉),计算结果是:SET TIME/FREQ LOAD STEP SUBSTEP CUMULATIVE1 45.108 1 1 12 52.591 1 2 23 68.922 1 3 34 76.373 1 4 45 91.398 1 5 56 126.49 1 6 67 140.43 1 7 78 152.15 1 8 89 156.77 1 9 910 160.79 1 10 10我想请问:1。
水箱装水时,怎么计算结果中有0频率啊?其中不等于0的频率是水箱体的频率么?2。
水箱不装水时的结果和装水时的相差怎么这么大啊?还是我理解错了么?企盼指教一个流固耦合建模的例子/prep7 !进行预处理模块et,1,30, !定义1号单元为Fluid30 流固耦合单元et,2,29 !定义2号单元为Fluid29平面流体单元et,3,30 , ,1 !定义3号单元为Fluid30流体介质单元et,4,63 !定义4号单元为Shell63壳体单元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 ㎝,半径为25 ㎝,壳体的壁厚为2 ㎝,cyl4,0,0,0.25,,5 !形成圆面k,9,0,0,0 ! 定义原点k,10,0,0,0.5lstr,9,10 !通过原点作直线adrag,5,6,7,8,,,9 !通过放样形成圆柱wpoff,0,0,0.1asel,s,,,2,5asbw,all,,, !移动工作平面与选取的侧面相切!重复上面操作,形成四个环肋面wpoff,0,0,-0.4 !工作平面回到原点位置上k,31,0.2,0,0.1 !定义环肋的方向点lsel,s,,,20 !选择要划分为环肋的线段latt,4,5,5,,31,40,1 !定义线段物理属性lesize,20,,,6 !划分数目secnum,1lmesh,20 !划分线段type,2 ! 选取第二种单元lsel,s,,,1,4 !选取线段lesize,all,,,10 !线段划分数目lesize,all,,,6amesh,1 !将通过工作平面的面1 进行划分esize,0.1,0 !选取沿面放样的网格大小mat,1 !定义放样的形成单元的物理属性vdrag,1,,,,,,9 !进行放样形成流体介质一加筋板在水中的模态分析再给大家一个实例!考虑结构在水中的自振频率:例子是一加筋板在水中的模态分析。
命令流如下:FINISH/CLEAR/FILENAME,plane/UNITS,SI/TITLE,plane/PREP7!*********ELEMENT DEFINE********ET,63,63ET,4,beam4et,30,fluid30!****MATERIAL DEFINE*********MP,EX,1,2.10E11MP,DENS,1,7850MP,NUXY,1,0.3mp,dens,30,1025mp,sonc,30,1500mp,mu,30,0.5!*******REAL CONSTANT***********r,30,1e-06r,50,0.05r,75,0.375e-02,0.78125e-06,0.000016406k,1k,4,1kfill,1,4,2,,1kgen,4,1,4,1,,1/3,,10a,1,2,12,11*do,i,0,2*do,j,0,2*10,10a,1+i+j,2+i+j,12+i+j,11+i+j*enddo*enddo!***************************fluid element****************k,100,-14.5,-14.5k,101,-14.5,15.5k,102,15.5,15.5k,103,15.5,-14.5k,140,-14.5,-14.5,30k,141,-14.5,15.5,30k,142,15.5,15.5,30k,143,15.5,-14.5,30a,100,101,102,103,4,14,24,34,33,32,31,21,11,1 a,1,2,3,4,103,100a,140,141,142,143a,100,101,141,140a,101,102,142,141a,142,143,103,102a,140,143,103,100a,14,24,34,33,32,31,21,11,1,2,3,4asel,u,,,1,FLST,2,8,5,ORDE,2FITEM,2,10FITEM,2,-17V A,P51Xnummrg,allallsMSHKEY,0MSHAPE,0esize,1lsel,s,loc,y,1/3lsel,r,loc,x,0,1lsel,r,loc,z,0latt,1,75,4lmesh,alllsel,s,loc,y,2/3lsel,r,loc,x,0,1lsel,r,loc,z,0latt,1,75,4lmesh,alllsel,s,loc,x,1/3lsel,r,loc,y,0,1lsel,r,loc,z,0latt,1,75,4lmesh,alllsel,s,loc,x,2/3lsel,r,loc,y,0,1lsel,r,loc,z,0latt,1,75,4lmesh,allasel,s,,,1,9aatt,1,50,63amesh,allallsMSHAPE,1,3desize,3vsel,s,,,1type,30 $mat,30 $real,30vmesh,allallsFINISH/solualls!**** 求解***********!*********************ANTYPE,MODAL!MODOPT,lanb,25,0MODOPT,UNSYMM,25,0SOLVEFINISH总是出现error 说矩阵不对称,不可以用lanb计算。