Ansys CFX 流固耦合分析
- 格式:doc
- 大小:812.00 KB
- 文档页数:22
ANSYS 14.0中Workbench提供了进行流固耦合(FSI)分析的模块,可以十分方便的对轴流叶轮机械进行气动载荷分析,包括最大变形量和等效应力分布。
1.进入ANSYS14.0 Workbench界面。
2.在左下角中的custom system模块中选择第一个流固耦合模块FSI:Fluid Flow(CFX)-staticstructural,双击。
3.屏幕中出现了FSI模块。
4.右击A5(solution)选择import solution,导入已经计算完毕的CFX结果.res文件。
5.导入结果后的界面如下图所示。
CFX部分已经完成了计算,所以不需要额外的设置。
6.双击B3(Geometry)进入结构分析的几何单元,初始单位选择meter。
7.导入一个叶片的几何实体,可以选择的几何文件类型很多,x_t、iges等等都可以。
在CFX中,我们通常计算的都是多个转子,多个叶片,但是在分析流固耦合时,只需导入自己关心的那个叶片就可以了。
8.然后点击Generate,就可以看到生成的叶片实体了。
8.关闭Geometry窗口回到Workbench截面,可以看到此时B3(Geometry)后已经变成了绿色的√,说明生成正确。
9.双击B4(model)进入。
可以看到Geometry、coordinate system、connections等项目前面已经是绿色的对号,不需要再进行设置。
10.单击mesh,在左下角的Details of mesh,如图进行设置。
10.右击mesh,选择generate mesh生成网格。
11.生成的叶片网格如图所示。
12.点击static structural ,选择工具栏中的support 下的fixed support,为叶片根部添加约束。
13.选中叶根面,点击左下角中的Apply,完成约束添加。
14.点击上工具栏中units,选择转速单位为RPM.15.如图所示添加转速16.按自己的算例输入转速。
ANSYS流固耦合分析成功的条件1.首先在建模和条件设置方面要按照这样的设置顺序:(1)选取流体单元,(打开keyopt(4)选项),建立流体模型,注意此处挖去固体所占的空间,然后分区划分流体场网格(好像在ls_dyna里面不要挖去固体所占空间),注意靠近挖去空间的部分网格应该细小些,还有若要采用remesh在计算中重划网格,一定要使用三角形单元(所有流体场)(2)流体场模型建立完成后,首先要在流固耦合的边界上施加流体耦合标签FSI,然后在在流体场区域施加必要的边界条件,诸如位移约束,速度、压力等等。
然后设置求解流体场的时间步长、求解时间、流体属性,打开ALE选项(瞬态分析)网格重画属性等等(3)这样的工作完成后,进入/prep7,加入固体单元,设置固体材料属性,在挖去的部分建立固体模型,划分固体网格,在固体网格与流体场接触的固体边界上施加流体耦合标签FSI,注意要和前面的number相同。
(4)为固体实际必要的约束条件(看是固体推动流体还是流体推动固体)(5)设置固体求解的时间步长和求解结束时间(6)设置流固耦合属性,(看是固体推动流体还是流体推动固体),求解时间步长和求解时间,收敛准则,迭代次数等等。
(7)保存求解。
-----------------------------总之,在流固耦合分析中,你最好要按着先流体后固体再耦合的属性设置顺序,流固耦合标签FSI要分别加在流固耦合边界的流体边界上和固体边界上,加在的顺序要按照上面所述。
在实际的建模中,流固耦合的边界上由于建模的原因会出现节点的重合现象,注意一定不要使用捏合节点的命令来将重合的节点变成一个,这个很重要。
固体单元一定要设置求解的时间步长和求解中止时间,时间步长一般和求解流体场和流固耦合的时间步长相等。
ALE+remesh选项是解决瞬态流固耦合问题的一个很重要的方法,流固耦合一般要伴随着流体的形状改变和位置的移动,因此首先启动ALE选项使流体与固体的耦合边界保持一致并规则化流场内部由于挤压而畸变的网格,其次若网格畸变的过于严重,就要启动remesh选项重新划分网格单元。
ANSYS流体(CFX)/结构(Structure)耦合计算流程本人最近在学习这方面的知识,对流固耦合问题有了初步的认识,现发在这里,和大家分享,并请求指正!在ANSYS的早期版本,ANSYS与CFX之间的流固耦合计算是单向耦合的,而从ANSYS10.0开始,ANSYS可以和CFX进行双向的流固耦合计算,即对一个包含固体和流体计算域的模型可以分别在ANSYS和CFX中同时进行计算,数据进行时时交换耦合;对于从ANSYS 传来的网格位移,CFX中可以自动进行网格变形。
一般单向耦合适合于结构形状对流体影响不大的情况,而当结构形状对流体影响很显著时就得用双向耦合。
在ANSYS和CFX之间进行流固耦合计算的过程如下:分别在ANSYS中建立结构域模型和在CFX中建立流体域模型,并对结构域模型和流体域模型分别划分有限元网格以及物理定义,之后会在CFX中针对流体域会生成*.def文件,在ANSYS中针对结构域生成*.in文件。
有了这两个文件后,启动ANSYS/CFX,分别指定*.def 文件和*.in文件开始ANSYS和CFX之间的双向耦合计算,在流固耦合计算中,定义流固界面,程序自动进行在流固界面进行平衡迭代,完成稳态和瞬态流固耦合分析。
图片附件: 流程.JPG (2006-4-4 17:53, 62.23 K)上图为流固耦合以及与sysnoise声学软件的耦合解决方案,其中红色框为ansys10.0提供的流固耦合的流程。
下图就是10。
0中进行流固耦合时的启动界面。
分别指定*.in文件和*.def文件,就可以进行双向的流固耦合计算。
程序同时启动ansys和CFX进行计算。
这与以前版本的流固耦合具有非常大的区别。
图片附件: 启动界面.JPG (2006-4-4 18:00, 92.11 K)。
CFX_流固双向耦合的实现实现流固双向耦合需要以下几个步骤:1. 网格生成:首先需要生成流体和固体模型的网格。
对于流体,可以使用常规的CFD网格生成软件(如Ansys ICEM-CFD)生成适当的流体网格。
对于固体,可以使用CAD软件生成固体模型,并通过网格生成软件(如Ansys Meshing)将其转换为固体网格。
2. 物理模型设定:根据实际情况,选择合适的流体和固体模型进行设定。
对于流体,可以选择使用Navier-Stokes方程来描述流体的运动。
对于固体,可以选择使用弹性力学方程进行模拟。
3.边界条件设定:对于流体和固体的边界条件进行设定。
对于流体,包括入口流速、出口压力、壁面摩擦等边界条件。
对于固体,包括固体的位移、力或者应力等边界条件。
4. 数值求解:根据设定的物理模型和边界条件,使用CFX软件进行数值求解。
CFX使用有限体积法对Navier-Stokes方程进行离散化,同时使用显式或隐式方法求解弹性力学方程。
5.耦合求解:在流固双向耦合中,流体和固体之间的相互作用需要通过迭代的方式求解。
首先,在给定流体的边界条件下,使用CFX求解流体部分的问题。
然后,在给定固体的边界条件下,使用CFX求解固体部分的问题。
接着,将固体的变形信息传递给流体,影响流体的边界条件。
再次使用CFX求解流体的问题,得到新的流场分布。
重复这个过程,直到流体和固体的解收敛。
6.结果分析:对求解得到的结果进行分析和后处理。
可以通过CFX提供的后处理工具,如应力和变形分布、速度和压力分布等来评估流固耦合模拟的效果。
值得注意的是,流固双向耦合模拟的实现通常需要较高的计算资源和时间。
同时,由于流固耦合问题的复杂性,对物理模型的设定以及边界条件的设定也需要经验和专业知识。
综上所述,CFX流固双向耦合的实现可以分为网格生成、物理模型设定、边界条件设定、数值求解、耦合求解和结果分析等几个步骤。
通过迭代的方式求解流固双向耦合问题,可以模拟流体和固体之间的相互作用,为工程实践提供有价值的参考。
流固耦合方法要实现流固耦合首先要确定固体模型和流体模型。
本设计的流体和固体模型都是在ANSYS Workbench中建立的(也可以用其他的三维建模工具)。
建模后将流体模型输出为IGES格式,然后导入到ICEM中化分网格备用。
固体模型则直接在Workbench中的Simulation中定义固支面和流固耦合面等,并划分网格。
最后输出inp文件。
最后利用CFX进行双向流固耦合计算。
具体步骤在下面分别叙述。
流体部分1.流体模型的建立启动ANSYS Workbench,点击New geometry开启DesignModeler。
在这里可以采用多种建模方法,我用的是直接在Create下拉菜单中选择最底部的Primitives 中的bend选项,直接建立一系列的扇形环柱体。
即可组成流体域的模型,在这里不再赘述。
2.流体网格的划分建模后将流体模型输出为IGES格式,然后导入到ICEM中。
在ICEM中通过点和线处理工具删掉一些多于的线并补上必要的点(分块的时候用来固定节点)。
然后创建体,在创建的提示栏中选择所有的面。
注意体形成后要多还几个角度观察,以确定solid点在整个模型的内部。
接下来就要分块,点击blocking按钮,选择整个实体为对象创建块。
之后将块沿边界切割成若干子块,删掉不需要的子块。
然后把快上的节点固定在相应的实体的点上,把块的边界线固定在相应的实体线上。
然后在块的特征边上设定节点数,化分网格。
3.网格输出网格划分完之后,需要将网格从ICEM CFD导出到计算软件ANSYS CFX中做计算,其具体的操作步骤如下:(1)选择File > Mesh > Load from Blocking;(2)选择File > Blocking > Save Multiblock Mesh,在出来的对话框上选择V olume。
(3)选择Output > Select solver,在出来的对话框上选择计算所用的软件ANSYS CFX,点击Okay。
CFX—Ansys单向流固耦合分析——散热器热分析1.计算模型模型一般选用Iges格式,单位mm,去掉直径小于φ30的孔及半径小于R10等细小特征,将散热元器件接触面分割出来。
2.元器件发热量IGBT损耗IGBT器件损耗包括IGBT、FWD稳态损耗及开关损耗。
根据I C=f(V CE)查出V CE 及I C的两个值,分别输入表一中,即可求出IGBT的稳态损耗P IGBT-DC,根据表二求出FWD的稳态损耗P FWD-DC 。
根据表三、表四求出IGBT的开关损耗,表五求出FWD的开关损耗。
根据已知的电阻、电流查出损耗值,然后根据实际电阻及已知电阻的比例关系求出实际的损耗值。
整流桥损耗根据公式求出I dc,然后查表得出整流桥损耗,其中I o为电流值,其他参数不变。
具体参见“IGBT、整流桥功耗计算(ECM※).xlsx” 及元器件手册。
热量输入采用热流密度,单位W/m2。
3.创建文件夹a)建立一个总文件夹,如“radiator(ECM43)”b)在radiator(ECM43)里建立IGS、hm、cfx及wb文件夹,IGS文件夹里放入模型数据,hm里放入网格数据,cfx里放入流体分析文件,wb文件夹里放入热分析文件 4.网格划分采用Hypermesh划分网格,必须使散热器网格与流体网格坐标一致。
4.1.启动“Hypermesh”在桌面上双击图标开始Æ所有程序ÆAltair Hyperworks 8.0 sr1Æhypermesh4.2.导入模型FileÆImportÆGeometryÆIGESGeometry ColorÆBy TopoShaded Geometry and Surface Edges4.3.程序参数设置数据格式选择PreferencesÆUser Profiles…选择PreferencesÆOptions,4.4.体网格划分3DÆtetrameshÆvolume tetraEnclosed volume:,Hint:Shift+左键拖动选择散热器全部表面点击开始网格划分,看到提示信息表示网格划分结束。
流固耦合FSI分析分析原理:流场采用CFX12,固体采用ANSYS12分别计算,通过界面耦合。
流体网格:流体部分采用HyperMesh9.0分网,按照流体分网步骤即可,没有特殊要求。
网格导出:CFX可以很好的支持Fluent的.cas格式。
直接导出这个格式即可。
流体的其余设置都在CFX-PRE中设置。
固体网格即设置:HyperMesh9.0划分固体网格。
设置边界条件,载荷选项,求解控制,导出.cdb文件。
实例练习:以CFX12实例CFX tutorial 23作为练习。
为节省时间,将计算时间缩短为2s。
网格划分:提取CFX tutorial 23中的实体模型到hm中,分别划分流体,固体网格。
分别导出为fluent的.cas格式和ansys的cdb格式。
流体网格如下:网格文件见:fluid.cas固体网格为:特别注意:做FSI分析时,ANSYS固体部分必须在BATCH下运行(即将.cdb文件导入ansys不需要任何操作就能直接计算出结果),所以导出的.CDB文件需要添加一个命令,在hm建立FSIN_1的set,以方便在.cdb中手动添加命令SF,FSIN_1,FSIN,1,具体位置在定义了节点集合FSIN_1之后。
另一个set:pressure用于施加压强。
这里还设置了一些控制卡片用于分析,当然也可以直接修改.cdb文件详细.cdb文件请参看plate.cdb将固体部分在ansys中计算一下,以确定没有问题。
通过ansys计算检查最大位移:最上面的点x向变形曲线至此,固体部分的计算文件已经准备好,流体网格需要导入CFX以进一步设置求解选项和耦合选项。
以下在CFX-PRE中进行设置由于固体模型已经生成,故不需要利用workbench,所以不必按照指南的做法。
启动workbench,拖动fluid flow(CFX)到工作区直接双击setup进入CFX-PRE 导入流体网格然后设置分析选项:注意:mechanical input file即是固体部分网格。
ANSY S -CFX 单向流固耦合分析的方法刘志远 郑 源 张文佳 司佳钧摘 要 在用Ansys 软件对风轮进行结构静力分析的过程中,无法从流体计算软件F LUE NT 中直接获取叶片在流场中所受的压力,即风施加在叶片上的瞬态压强值。
此类问题的研究属于流固耦合的范畴,也是目前流、固体力学研究领域比较前沿的课题。
通过研究ANSY S -CFX 组合软件,发现了分析单向流固耦合问题的方法,从而在3D 软件Ansys W orkbench 中实现了对风轮受力变形更合理、更精确的数值模拟。
关键词 ANSY S -CFX 流固耦合 垂直轴风力机 静力分析中图分类号 T V73411 文献标识码 A 文章编号 100726980(2009)022******* 耦合场分析是考虑两个或两个以上工程学科(物理场)间相互作用的分析。
例如流体与结构的耦合分析,即流固耦合(Fluid S tructure Interaction ),流体流动的压力作用于结构,结构将产生变形,而结构的变形又影响了流体的流道,因此是相互作用的问题。
目前,在工程学科中,特别是流体动力学领域中,越来越多的实际问题需要进行耦合场的模拟分析,例如水轮机的叶片在水流中的变形情况,风机的叶片在风场中的变形情况等。
因为受软件开发水平的限制,很多软件只能完成单一物理场的模拟,而不能够实现多物理场的耦合。
为了能够实现对某一物理模型的多场耦合分析,国内的一些高校和研究机构通常对相关软件进行二次开发,例如对Ansys [1]和F LUE NT 进行程序接口的二次开发,来解决不同软件之间的数据交换问题,但这种方法不仅要求开发者具有相当高的编程水平,同时也需要耗费大量的时间,而且这些机构开发出来的程序也往往只适用于他们自己所研究的领域,所以在推广上具有很大的局限性。
虽然实现流固耦合分析的软件很多,方法也不少,但由于受方方面面因素的制约,国内在这方面的资料却很稀缺,本文以一螺旋S 型风力机的叶片在风场中旋转时的某一瞬时状态为例,来介绍如何通过Ansys W orkbench [2](ANSY S -CFX )来实现单向流固耦合(FSI )分析的方法。
ANSYS流固耦合分析9.2 问题的描述(2)进入前处理器,定义所需的单元类型。
对于该问题而言,模型中只含有Plane 42一种单元,由于本计算是针对二维平面应变问题,因而需将Plane 42的KEYOPT(3)设置为2。
对应的命令流如下:/prep7et,1,plane42KEYOPT,1,3,2 !平面应变选项(3)定义模型中的材料参数。
模型中共有7种材料,包括基岩和不同饱和度下的滑体材参。
定义材参的命令如下:mp,ex,1,20e9 !岩体弹性模型mp,prxy,1,0.22 !岩体泊松比mp,dens,1,2600 !岩体密度mp,ex,2,3e9 !滑体初始参数mp,prxy,2,0.35mp,dens,2,2140mp,ex,3,2.7e9 !饱和度为35%mp,prxy,3,0.35mp,dens,3,2180mp,ex,4,2.4e9 !饱和度为50%mp,prxy,4,0.35mp,dens,4,2210mp,ex,5,2.1e9 !饱和度为65%mp,prxy,5,0.35mp,dens,5,2240mp,ex,6,1.8e9 !饱和度为80%mp,prxy,6,0.35mp,dens,6,2270mp,ex,7,1.5e9 !饱和度为100%的弹性模量mp,prxy,7,0.35 !饱和度为100%的泊松比mp,dens,7,2300 !饱和度为100%的密度save !存储数据库(4)建立几何模型并划分单元。
由于模型的几何形状较为复杂,故采用自底向上的建模方式,即先生成点,然后生成面,最后划分单元。
k,1,0,0,0 !根据坐标生成关键点k,2,0,505,0k,3,562,505,0k,,621,459,0k,,658,433,0k,,693,393,0k,,802,313,0k,,850,303,0k,,892,273,0k,,913,253,0k,,930,243,0k,,1034,233,0k,,1186,228,0k,,1216,223,0k,,1232,208,0*do,i,1,14 !以循环的方式创建直线,共画了14条l,i,i+1*enddok,16,660,387,0 !再创建关键点k,,770,248,0k,,827,230,0k,,888,217,0k,,930,216,0k,,1034,213,0k,,1186,209,0l,4,16 !再根据关键点创建直线l,16,17l,17,18l,18,19l,19,20l,20,21l,21,22l,22,15nummrg,all !合并所有元素numcmp,all !压缩所有元素编号k,23,1243,198,0 !再创建关键点k,,612,380,0k,,745,214,0k,,884,181,0k,,1242,174,0l,3,24 !再根据关键点创建直线l,24,25l,25,26l,26,27l,27,23l,23,15nummrg,all !合并所有元素numcmp,all !压缩所有元素编号k,28,1610,198,0 !再创建关键点k,,1610,0,0k,,1243,0,0k,,884,0,0k,,745,0,0k,,612,0,0l,6,16 !再根据关键点创建直线l,7,17l,8,18l,9,19l,11,20l,13,22l,16,24l,17,25l,19,26l,24,33l,25,32l,26,31l,27,30l,1,33l,30,31l,31,32l,32,33l,23,28l,28,29l,29,30saveal,4,5,15,29 !根据线创建面al,29,6,16,30al,7,17,30,31al,8,18,31,32al,9,10,19,32,33al,11,12,20,21,33,34al,13,14,22,34al,3,15,35,23al,35,16,24,36al,36,17,18,37,25al,19,20,21,22,26,27,28,37al,1,2,23,38,42al,38,39,24,45al,39,40,25,44al,40,41,26,43al,46,47,48,27,41allsel,all !选择所有元素type,1 !选择要划分单元的类型mat,1 !选择要划分单元的材料模型mshape,1,2d !设置成平面三角形单元mshkey,0 !设置成自由方式ESIZE,10,0, !设置单元大小amesh,1,11,1 !对面1到面11划分单元ESIZE,,, !设置单元大小为默认值allsel,all !选择所有元素amesh,12,16,1 !划分面12到16 savefinish!划分好后的有限元网格如图9-2所示。
CFX user Fortran: 在每个时间步输出结果问题是这样的,用CFX user CEL完成变量的同时,还希望CFX输出每个时间步的变量值(比如,每个时间步的压力值)到一个文本文件。
虽然在程序中作了一些控制(以避免在每个coefficient loop都输出结果),但是CFX总是在一个时间步输出同样的变量多次。
User CEL的功能其实是和CFX的内置CEL一样的,只是用户化的CEL。
User CEL会在“需要时”,被实时引用。
如果user CEL是用来定义边界条件,在一个coefficient loop,user CEL会被运行多次。
这样结果也会被多次重复输出。
其实CFX user Fortran有两种,第一种是user CEL,它的主要功能是表达式的计算;第二种是junction box routine,它的主要功能是控制计算流程。
user CEL会被实时引用,而junction box routine只在用户设定的运行点运行。
下图是从CFX帮助里复制的流程图。
要完成上面所说的功能,只要定义个junction box routine,然后把运行点设成End of Time Step (transient only)就可以了。
/?p=116}Fri, 14 Dec 2007 18:34:57 +0000Fluent: 计算叶轮功率问题是这样的,Fluent可以不可以计算驱动搅拌器叶轮所需要的功率;如果可以的话,是不是要做面积分。
按照物理学上的定义,对于旋转的系统,功率=力矩*转速。
通过表面积分(把每个壁面单元,wall face,的功率相加)可以得到整个系统的功率。
但是,这个方法比较麻烦,需要用UDF计算,可能比较费时间。
既然功率=力矩*转速。
转速我们是已经知道的。
力矩可以在Fluent->Report->Forces…上计算得到。
因此,如果叶轮转动轴是和x,y,z轴其中之一方向一致的话,我们可以直接用Fluent图形界面进行计算。
在航空中ANSYS CFX流固耦合模拟的应用如今计算流体力学(CFD)已经发展成为分析工业设备外部和内部流动的可靠工具,其所面临的新挑战是对于涉及不同物理现象的多物理场的模拟。
一个重要的例子是流动与周围固体结构的干扰。
干扰可以是流体作用力与固体变形的力学耦合,也可以是流固界面之间温度和热通量的热耦合。
一个典型的例子是机翼或叶片颤振力学耦合系统的数值模拟。
在本文的研究中,ANSYS公司的两个软件包ANSYS和CFX被用于结构和流动力学耦合的模拟。
ANSYS是多用途非线性的有限元求解器,用于计算固体结构和非固体结构(例如,静电场、静磁场、声学)。
CFX是通用的CFD代码,以高鲁棒性和高精度的流动数值算法、高级湍流模式和多种复杂物理模型而著称。
ANSYS CFX软件介绍在CFX中,NS方程组采用守恒形式的有限体积法来离散,时间采用隐格式,可以计算混合网格和非结构网格,网格单元可以是六面体型、棱柱型、楔型和四面体型。
在每个网格节点周围构造控制体,通量通过位于两个控制体界面上的结合点来计算。
离散方程采用有界的高精度对流格式来求解。
通过Rhie和Chow 的算法来计算质量流量,以保证压力速度耦合。
离散方程组通过由Raw发展的耦合代数多重网格法求解。
该方法的数值能力随参与计算的网格节点的数量增加而线性增加。
定常计算采用时间迭代法,直到达到用户指定的收敛标准。
对于非定常计算,迭代程序在每个时间步内更新非线性系数,而时间步由外层循环来推进。
由于力学耦合,将导致流体和固体之间的界面发生移动。
因此,离散方程必须被拓展以允许网格移动和网格变形。
这种拓展通过空间守恒律来实现。
壁面网格节点界面的移动需要重新计算求解域内部网格节点的位置,可以通过求解描述网格变形的拉普拉斯方程来实现,这类似于网格光顺所做的操作。
它是描述动网格运动的经典粘弹动力学方程的简化形式。
如果网格严重变形,光顺网格的方法不足以提供高质量的网格。
在这种情况下,必须建立拥有不同网格拓扑结构的新网格,在下一时间步,通过二阶插值将求解变量插值到新网格节点上。
CFX 流固双向耦合的实现2010.6.4 精仪 设计所 王昊森○流固耦合基本问题流固耦合问题一般分为两类,一类是流‐固单向耦合,一类是流‐固双向耦合。
单向耦合应用于流畅对固体最用后,固体变形不大,即流畅的边界形貌改变很小,不影响流畅分布的,可以使用流固单向耦合。
先计算出流畅分布,然后将其中的关键参数作为载荷加载到固体结构上。
典型应用比如小型飞机按刚性体设计的机翼,机翼有明显的应力受载,但是形变很小,对绕流不产生影响。
当固体结构变形比较大,导致流畅的边界形貌发生改变后,流畅分布会有明显变化时,单向耦合显然是不合适的,因此需要考虑固体变形对流畅的影响。
两者相互作用,最终达到一个平衡状态(稳态问题中)。
比如大型客机的机翼,上下跳动量可以达到5米,以及一切机翼的气动弹性问题,都是因为两者相互影响产生的。
因此在解决这类问题时,需要进行流固双向耦合计算。
○双向流固耦合的基本方式双向流固耦合的基本思路是:分别计算流场和固体结构,然后通过中间平台交换耦合量。
每次大迭代中,分别进行一次流体计算和固体计算,并交换2次数据(每个方向各一次),直到最终收敛。
这种通过交换耦合边界的耦合量实现的双向耦合称为松耦合方式(即非耦合方程方式求解)。
(1) MPCCIMPCCI(Mesh-based parallel Code Coupling Interface)是德国法兰克福SCAI研究中心开发的针对耦合问题的工程软件。
本身没有计算功能,是作为数据交换和控制台软件存在。
具体是通过计算机内部的口通讯的监听和操作来实现数据的交换和求解器的控制。
其基本框架与CFX-ANSYS耦合方式相同。
下图为MPCCI的框架结构,示一般流固双向耦合的实现方式。
MPCCI框架结构MPCCI3.06版本可以使用ANSYS11.0以及FLUENT6.3版本。
MPCCI4.0可以使用ANSYS12.0和FLUENT12版本。
(2) CFX-ANSYS(12.0版本)CFX是集成于ANSYS中的流体计算软件,虽然功能较fluent差,但是目前只有CFX可以不借助第三方软件与ANSYS实现双向流固耦合的计算。
流固耦合FSI分析
分析原理:流场采用CFX12,固体采用ANSYS12分别计算,通过界面耦合。
流体网格:流体部分采用HyperMesh9.0分网,按照流体分网步骤即可,没有特殊要求。
网格导出:CFX可以很好的支持Fluent的.cas格式。
直接导出这个格式即可。
流体的其余设置都在CFX-PRE中设置。
固体网格即设置:HyperMesh9.0划分固体网格。
设置边界条件,载荷选项,求解控制,导出.cdb文件。
实例练习:
以CFX12实例CFX tutorial 23作为练习。
为节省时间,将计算时间缩短为2s。
网格划分:提取CFX tutorial 23中的实体模型到hm中,分别划分流体,固体网格。
分别导出为fluent的.cas格式和ansys的cdb格式。
流体网格如下:
网格文件见:fluid.cas
固体网格为:
特别注意:
做FSI分析时,ANSYS固体部分必须在BATCH下运行(即将.cdb文件导入ansys不需要任何操作就能直接计算出结果),所以导出的.CDB文件需要添加一个命令,在hm建立FSIN_1
的set,以方便在.cdb中手动添加命令SF,FSIN_1,FSIN,1,具体位置在定义了节点集合FSIN_1之后。
另一个set:pressure用于施加压强。
这里还设置了一些控制卡片用于分析,当然也可以直接修改.cdb文件
详细.cdb文件请参看plate.cdb
将固体部分在ansys中计算一下,以确定没有问题。
通过ansys计算检查最大位移:最上面的点x向变形曲线
至此,固体部分的计算文件已经准备好,流体网格需要导入CFX以进一步设置求解选项和耦合选项。
以下在CFX-PRE中进行设置
由于固体模型已经生成,故不需要利用workbench,所以不必按照指南的做法。
启动workbench,拖动fluid flow(CFX)到工作区
直接双击setup进入CFX-PRE 导入流体网格
然后设置分析选项:
注意:mechanical input file即是固体部分网格。
再新建一个流体,取名fluid。
设置domain
添加边界条件
取名为interface设置流固耦合界面,对应为abc。
这就是流固耦合界面的设置过程。
同理,建立sym1
Sym2
这个选项默认为no slip 的wall,最普通的那种,不必特殊设置初始化:
求解控制
输出控制:
Output variables list
看清楚字母,别搞错了!
Monitor板最上面的一个角点,也是ansys中计算关心的点
如这个十字号
保存以上设置即可。
返回到workbench,双击solution启动计算器。
不用更改设置,直接计算即可。
以上过程小心仔细,一般不会出现错误的。
至此,几分钟后即可计算完成。
先看收敛曲线
再看monitor点
曲线跟单独的板的振动十分相似,说明计算是正确的。
后处理:
再次返回workbench,双击results,启动后处理。
先查看板的应力
是对的
以前我用CFX11做过同样的例子,所以以下提供CFX11中的结果,与cfx12类似。
可用查看速度矢量。
从以上结果来看,分析是完全正确的!
动画制作略。
关于后处理,大家可以自己参考官方帮助做一下。