PFC手册例子耦合FLAC3D

1．P146，通过ball命令生成边界，球的半径一定。

2．A finite wall iscreated when the face keyword appears on the WALL command 有限平面是WALL command follow face keyword

流固耦合手册现场翻译

1 流体－固体耦合与单相渗流 1.1介绍 FLAC3D模拟了流体流过可渗透的介质，例如土体。渗流模型可以独立于通常FLAC3D的固体力学计算，而只考虑渗透；或者为了描述流体和固体的耦合特性，与固体模型并行计算。例如，固结是一类流固耦合的现象，在固结过程中孔隙压力逐渐消散，从而导致了土体的位移。这种行为包含了两种力学效应。其一，孔隙水压力的改变导致了有效应力的改变，有效应力的改变影响了固体的力学性能，例如有效应力的降低可能引发塑性屈服；其二，土体中的流体对孔隙体积的变化产生反作用，表现为孔隙水压力的变化。 本程序可以不仅可以解决完全饱和土体中的渗流，也可以分析有浸润线定义的饱和与非饱和区的渗流计算。该条件下，浸软面以上的土体的孔隙水压力为零，气体是被动的（气体的压力考虑成负的）。这种方法用于颗粒比较粗的毛细现象可以忽略土体。 渗流分析中有如下的特征： 1.对应于渗流各向同性和各向异性渗透性采用两种不同的流体传输定律。渗流区域中的不可渗透 的区域用流体的null材料定义。 2.不同的zone可以赋予不同的渗流模型(isotropic, anisotropic or null)和属性。 3.流体压力，涌入量，渗漏量和不可渗透边界都可以定义。 4.土体中可以加入抽水井，考虑成点源或者体积源。 5.计算完全饱和土体中的渗流问题，可以采用显式差分法或者隐式差分法；而非饱和渗流问题只 能采用显式差分法。 6.渗流模型可以和固体力学模型和传热模型一起使用。在耦合问题中，可以考虑饱和材料的压缩 和热膨胀。 7.流体和固体的耦合程度依赖于土体颗粒（骨架）的压缩程度，用Biot系数表示颗粒的可压缩 程度。（即用biot系数确定颗粒变形，来模拟力学与流体的耦合） 8.用线性热膨胀系数和非排干导热系数来考虑热耦合计算。 9.对于热-流体-流动计算是基于线性理论，假设材料属性为常数，不考虑对流。流体和固体温度 局部平衡。如果要模拟非线性行为，需要制定孔隙压力和用fish函数确定材料属性。 由于循环荷载引起的动水压力和液化问题也可以用FLAC3D模拟。FLAC3D不考虑毛细现象，土体颗粒间的电、化学作用力。然而，可以根据土体的局部饱和度，孔隙率，或者其他的变量，通过编写一段FISH语言来考虑这种力。类似的，由于液体中溶解了空气而引起的液体刚度变化，也不能显式的模拟，而通过FISH将液体刚度表示为压力，时间和其他变量的函数。 这以章节可以分为七个主要部分： 1.数学模型描述（章节1.2）和相应的数值方法（章节1.3）（单相渗流和流固耦合计算）。 2.计算模式和渗流分析相关的命令。（章节1.4） 3.渗流分析所需要的材料属性的单位。（章节1.5） 4.不同边界条件，初始条件，和流体的源和汇的的描述。（章节1.6） 5.求解单相渗流问题和流固耦合问题的推荐方法。（章节 1.7）建议先练习这里的例子再做自己 的流体分析。 6.（章节1.8）提供了几个验证算例，演示了用falc3d的流体计算方法的准确性。 7.（章节1.9）总结了与流体计算的相关命令。 用户在尝试解决流固耦合问题之间，强烈建议先熟悉一下FLAC3D求解固体力学问题的步骤。流固耦合的力学行为通常非常复杂，需要用户对计算结果是否真确需要远见和判断。在开始计算一个大的项目支前，有必要在一个网格比较少的例子上作试验，尝试不同的边界条件和建模策略。“浪费”在试验上的时间，必定可以通过计算时间来弥补回来。

Ansys CFX 流固耦合分析

流固耦合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)到工作区

三个流固耦合分析实例

length=2 !定义体各种变量参数，长宽高 width=3 height=2 /prep7 et,1,63 !选用壳模型 et,2,30 !选用FLUID30单元，用于流固耦合问题r,1,0.01 增加实常数，壳厚为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 !建立长方体 esize,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,all,,fsi !定义流固耦合界面

alls !选择所有项 solv !求解 fini /post1 !后处理 set,first plnsol,u,sum,2,1 !显示图形 fini /PREP7 !定义壳材料与性质 !壳元素与材料 ET,1,shell63 $MP,EX,1,201E9 $MP,prxy,1,0.26 $MP,dens,1,7.85E3 $r,1,0.006 !流体元素与材料 ET,2,FLUID80 $MP,EX,2,1.5e9 $MP,DENS,2,0.84e3 $mp,visc,2,1.0e-10 !以下这个keyoption怎么用? 如过用1，就会显示[Element 877 may not have a positive Z coordinate IF KEYOPT(2) = 1.]，显示这个错误代表要做什么修正吗?所以我暂时用KEYOPT(2) = 0就可以跑。 KEYOPT,2,2,0 !建立壳关键点 K,1,10,0,0 $K,2,10,0,12 !建立中心线关键点 k,3,0,0,0 $k,4,0,0,20 !定义壳壁线 L,1,2 $L,1,3 !以关键点3,4为中心线旋转360度生成壳体 AROTAT,all,,,,,,3,4,360 !划分壳体网格 AATT,1,1,1 $esize,2 $mshape,0,3D $mshkey,2 $amesh,all $alls !延伸出水位体积 VEXT,2,8,2,0,0,10,0,0,0 $vglue,all

ansys14workbench血管流固耦合分析实例

Ansys14 workbench血管流固耦合实例 根据收集的一些资料，进行学习后，试着做了这个ansys14workbench的血管流固耦合模拟，感觉能够耦合上，仅是熟悉流固耦合分析过程，不一定正确，仅供参考，希望大家多讨论。谢谢！ 1、先在proe5中建立血管与血液流体区的模型（两者装配起来），或者直接在workbench中建模。 图1 模型图 2、新建工程。在workbench中toolbox中选custom system，双击FSI: FluidFlow(fluent)->static structure. 图2 计算工程 3、修改engineering data，因为系统缺省材料是钢，需要构建血管材料，如图3所示。先复制steel，而后修改密度1150kg/m3，杨氏模量4.5e8Pa，泊松比0.3，重新命名，最后在主菜单中点击“update project”保存.

图3 修改工程材料 4、模型导入，进入gemetry模块，import外部模型文件。 图4 模型导入图 5、进入FLUENT网格划分。 在workbench工程视图中的Mesh上点击右键，选择Edit…，如图5所示，进入网格划分meshing界面，如图6所示。我们这里需要去掉血管部分，只保留血液几何。

图5 进入网格划分

图6 禁用血管模型 6、设置网格方法。 默认是采用ICEM CFD进行网格划分，设置方式如图7所示，截面圆弧边分为12份，纵截面的边均分为10份，网格结果如图8所示。另外在这个界面中要设置边界的几何面，如inlet、outlet、symmetry 图7 设置网格划分方式 图8 最终出网格

流固耦合分析成功的条件

流固耦合分析成功的条件 作者：mbrslydddd 本人近来查看好多CFD分析的论坛，发现好多网友都在为流固耦合而头大，且网上也没有完整的、系统的案例可供参考，目前小弟也在做流固耦合的问题，即为流体输送固体颗粒的过程的模拟，好不容易找到了完整的解决办法，不敢独享，特奉献于此。 首先在建模和条件设臵方面要按照这样的设臵顺序： 1)选取流体单元，（打开keyopt（4）选项）,建立流体模型，注意此处挖去固体所占 的空间，然后分区划分流体场网格（好像在ls_dyna里面不要挖去固体所占空间），注意靠近挖去空间的部分网格应该细小些，还有若要采用remesh在计算中重划网 格，一定要使用三角形单元（所有流体场）； 2)流体场模型建立完成后，首先要在流固耦合的边界上施加流体耦合标签FSI，然后 在在流体场区域施加必要的边界条件，诸如位移约束，速度、压力等等。然后设臵 求解流体场的时间步长、求解时间、流体属性,打开ALE选项（瞬态分析）网格重 画属性等等； 3)这样的工作完成后，进入/prep7，加入固体单元，设臵固体材料属性，在挖去的部 分建立固体模型，划分固体网格，在固体网格与流体场接触的固体边界上施加流体 耦合标签FSI，注意要和前面的number相同； 4)为固体实际必要的约束条件（看是固体推动流体还是流体推动固体）； 5)设臵固体求解的时间步长和求解结束时间，设臵流固耦合属性，（看是固体推动流 体还是流体推动固体），求解时间步长和求解时间，收敛准则，迭代次数等等； 6)保存求解。 总之，在流固耦合分析中，你最好要按着先流体后固体再耦合的属性设臵顺序，流固耦合标签FSI要分别加在流固耦合边界的流体边界上和固体边界上，加在的顺序要按照上面所述。在实际的建模中，流固耦合的边界上由于建模的原因会出现节点的重合现象，注意一定不要使用捏合节点的命令来将重合的节点变成一个，这个很重要。固体单元一定要设臵求解的时间步长和求解中止时间，时间步长一般和求解流体场和流固耦合的时间步长相等。 ALE＋remesh选项是解决瞬态流固耦合问题的一个很重要的方法，流固耦合一般要伴随着流体的形状改变和位臵的移动，因此首先启动ALE选项使流体与固体的耦合边界保持一致并规则化流场内部由于挤压而畸变的网格，其次若网格畸变的过于严重，就要启动remesh 选项重新划分网格单元。 看一个流固耦合时间步长设臵的是否合理的方法是在求解中若发现软件频频调动remesh进行网格重画，那这个时间步长基本上处于最大的时间步长上，若很久出现一次或不出现，说明你的时间步长太小，可以适当提高，但是正常的情况下，软件只允许在一个求解过程中最多有100次网格重画，若超过，将退出求解。 先写到这些，匆匆草就，可能有很多缺点和错误，大家指正（注意若固体推动流体且又不考虑固体的形变和应力等等，可以直接按照帮助文档上ALE一章） 还有，若流体推动固体的分析中，固体发生了大位移，那么一定要在/solu中打开大变形分析选项，否则求解不收敛而退出。

流固耦合 - 基础知识

流固耦合概念 流固耦合，是研究可变形固体在流场作用下的各种行为以及固体变形对流场影响这二者相互作用的一门科学。它是流体力学(CFD) 与固体力学(CSM) 交叉而生成的一门力学分支，同时也是多学科或多物理场研究的一个重要分支。流固耦合力学的重要特征是两相介质之间的相互作用，变形固体在流体载荷作用下会产生变形或运动。变形或运动又反过来影响流体运动，从而改变流体载荷的分布和大小，正是这种相互作用将在不同条件下产生形形色色的流固耦合现象。当你研究的问题，不仅涉及到了流场的分析，还涉及到了结构场的分析，而且二者之间存在着明显的相互作用的时候，你就考虑进行流固耦合分析。 流固耦合求解方法 流固耦合问题分析根据流体域和固体域之间物理场耦合程度的不同，可分为强耦合和弱耦合，对应的求解方法分别为直接解法和分离解法。直接解法通过将流场和结构场的控制方程耦合到同一方程矩阵中求解，即在同一求解器中同时求解流固控制方程，理论上非常先进，适用于大固体变形、生物隔膜运动等。但在实际应用中，直接法很难将现有的计算流体动力学和计算固体力学技术真正结合到一起。 另外，考虑到同步求解的收敛难度以及耗时问题，直接解法目前主要应用于模拟分析热－结构耦合和电磁－结构耦合等简单问题中，对于流体－结构耦合只进行了一些非常简单的研究，还难以应用在实际工程问题中。而弱流固耦合的分离解法是分别求解流体和固体的控制方程，通过流固耦合交界面进行数据传递。该方法对计算机性能的需求大幅降低，可用来求解实际的大规模问题。 目前的商业软件中，流固耦合分析基本都采用分离解法。ANSYS很早便开始进行流固耦合的研究和应用，目前ANSYS中的流固耦合分析算法和功能