ANSYS Workbench LS-DYNA流固耦合方法应用
ANSYS Workbench LS-DYNA流固耦合方法应用贮液容器（含塑料瓶）广泛应用于化工、食品包装、储运等领域。由于容器（含塑料瓶）在运输和使用过程中常常会因为跌落或碰撞冲击导致破损而造成损失和污染，因此，研究贮液容器（含塑料瓶）在跌落碰撞过程中的力学行为，对认识容器（含塑料瓶）跌落碰撞损伤机理，优化容器（含塑料瓶）结构，提高其安全性和使用价值意义重大。．
贮液容器的跌落是一个典型的流固耦合问题，可采用LS-DYNA的ALE算法（任意拉格朗日欧拉算法）进行模拟。下面以一个封闭的装水水箱为例，介绍ANSYS Workbench LS-DYNA 分析此类型跌落问题的方法和步骤：
1.建立几何模型
调用ANSYS Workbench中的LS-DYNA模块，如图1所示。然后使用ANSYS的CAD工具DesignModeler建立几何模型，如图2所示。
图1 调用Workbench LS-DYNA图2 DesignModeler中建立几何模型
2.生成K文件
双击进入“Model”后，对模型进行网格划分、边界条件设置、速度设置和分析设置，如图3所示。设置完成后点击“solve”求解，生成K文件，如图4所示。
图3调用Workbench LS-DYNA图4DesignModeler中建立几何模型
3.编辑K文件
通过Workbench LS-DYNA生成的K文件中关键字是不够完善的，并不能直接递交LS-DYNA求解器进行求解。K文件中所欠缺的一些关键字，在流固耦合分析中是必不可少的，如空材料的定义、跟随坐标系的定义、空白域的定义以及状态方程的定义等。
3.1重要关键字释义
（1）LS-DYNA程序提供了运动的多物质ALE网格，可以方便地为多物质ALE算法定义跟随坐标系
*ALE_REFERENCE_SYSTEM_NODE
*ALE_REFERENCE_SYSTEM_GROUP
（2）定义空材料和状态方程的关键字
*MAT_NULL*EOS
（3）初始化空白域的关键字
*INITIAL_VOID_PART
（4）结构和流体之间耦合的关键字
*CONSTRAINED_LAGRANGE_IN_SOLID
（5）单元算法定义（单点积分的单物质加空白材料）的关键字
*SECTION_SOLID_ALE ELF0RM=12
（6）在重力作用下产生下落的关键字
*LOAD_BODY
……
3.2关键字编辑方法
关键字的编辑或修改一般有两种方法，一种是直接在ls-prepost中对关键字进行编辑设置，如图5所示；另一种是在文本编辑器UltraEdit中对关键字进行编辑或修改，如图6所示。当然也可以使用TEXT文本编辑器来编辑K文件，不过TEXT在此方面的专业性不及UltraEdit，尤其对于较大的K文件，使用TEXT打开将会很慢。
图5 ls-prepost中编辑关键字图6UltraEdit中编辑关键字
3.3关键字检查
关键字编辑完成后，很难确保所有关键字编写无误，若直接递交LS-DYNA求解器求解，计算一段时间后可能会报错。这样往往会浪费我们的时间，可以在求解之前先对关键字进行检查，有针对性的将错误提前排除掉，具体检查方法如图7所示。图中提示关键字LOAD存在错误，对LOAD进行重点查看，错误排查后如图8所示，然后再提交至LS-DYNA求解器计
算。
图7关键字错误检查图8 关键错误排查
4.后处理查看结果
在ls-prepost中查看计算结果，可将计算得到的d3plot文件直接拖入至ls-prepost，勾选“Assembly and Select Part”中的“Fluid(Ale)”，便可查看装水水箱整个跌落过程中液体水的流动情况，如图9所示。
图9 后处理结果查看
流固耦合应用研究进展
文章编号:1671-3559(2004)02-0123-04 收稿日期:2003-12-03 基金项目:山东省科学技术发展计划资助项目(012050107);山 东省自然科学基金资助项目(Y 2002F19) 作者简介:郭术义(1971-),男,山东济南人,山东大学机械工 程学院博士研究生。 流固耦合应用研究进展 郭术义,陈举华 (山东大学机械工程学院,山东济南250061) 摘 要:流固耦合力学是一门新兴学科。本文简要介绍了该学科的典型应用进展情况,总结了各种研究中的典型方程、数值解法,展望了进一步发展的趋势。关键词:流固耦合;数值模拟;展望中图分类号:O35112;O34717 文献标识码:A 流固耦合力学是一门比较新的力学边缘分支, 是流体力学与固体力学二者相互交叉而生成的。它的研究对象是固体在流场作用下的各种行为以及固体变形或运动对流场的影响。流固耦合力学的重要特征是两相介质之间的相互作用:固体在流体动载荷作用下产生变形或运动,而固体的变形或运动又反过来影响到流场,从而改变流体载荷的分布和大小。总体上,流固耦合问题按耦合机理可分为两大类:一类的特征是流固耦合作用仅仅发生在流、固两相交界面上,在方程上耦合是由两相耦合面的平衡及协调关系引入的;另一类的特征是流、固两相部分或全部重叠在一起,耦合效用通过描述问题的微分方程来实现。本文就流固耦合问题的两大分类中三种基本情况进行了讨论。 1 流固耦合典型应用 流固耦合作用的研究在航空、航天、水利、建筑、石油、化工、海洋以及生物领域都有着十分重要的意义。如液体晃动对火箭飞行稳定性的影响,大型贮液管在地震激励作用下产生的流固耦合作用,液体湍振对输液管道的影响。本文就如下三个大方面进行了总结。1.1 输流管道流固耦合 流体引起输流管道振动的研究最初来源于横跨 阿拉伯输油管道振动的分析[1]。管道在众多的工业领域中应用十分广泛,作用极其重要。但是,在管道 内流体流动状态的微弱变化往往引起在工作过程中的湍振现象,诱发流体、管道之间的耦合振动,动力学行为相当复杂。这使得人们很早就开始了这方面的研究,Paidoussis M P [2]是其中最具有代表性的。输流管道的振动问题之所以能引起学者的兴趣,除因为该问题的广泛工业背景和现实意义之外,还因为输流管道虽然是最简单的流固耦合系统,但它却涉及了流固耦合的大多数问题,并且它的物理模型简单,系统比较容易实现,因而便于理论与试验的相互协同。 考虑因素侧重面的不同,输液管道非线性运动方程有几种类型[3-5],它们之间有一定的差别。它们的基本假设都是:流体无粘且不可压;管道作为梁模型来处理;管道只是在平面内振动。尽管输流管道的非线性动力问题受到50多年极为广泛的研究,但至今尚没有一个公认的模型。文[6]建立的4个独立变量(轴向位移、横向位移、流速和压力)的全耦合模型(耦合形式包含摩擦耦合、P oiss on 耦合、结合部耦合以及管道轴向和横向运动的耦合)在众多的非线性分析模型中是一个较为完整的模型。 m ¨u +m f [ υf (1+u ′)+2υf u ′+υ2 f u ″+ ωυ′f ]+ P (υf + u )/c 2F -[(1-2υ)P (1+u ′)]′+4f ρf ρ′?υ2f /DK -gm f (1-2υ)(1+u ′)ω′-EI (7ω″ω +ω′ ω )-E A p (2u ″+6u ′u ″+2ω′ω″ )/2=0(1)m ¨ω+m f [ υf (1+ω′)+2υf ω′+υ2f u ″+ω″υ2 f ]+ P (υf + ω)/c 2F -[(1-2υ)P ω′]′-gm +EI ω″″-EI (u ′ω′+6u ″ω +4u ′ω ′)-E A p (u ″ω′+u ′ω″ )=0(2) P /c 2F +m f [(1-2υ)( u +υf )u ″- u ′+υ′f ]-m f (1-2υ)( u ′+u ′ u ′+ω′ ω′ )=0(3)P ′+m f (¨u + υf )+m f ¨ωω′+gm f ω′+Df ρf υ2 f /2=0 (4)随着对输流管道问题研究的深入,各种不同的 分析计算方法也相继被提出。其中有限元法(FE M ) 第18卷第2期2004年6月 济南大学学报(自然科学版) JOURNA L OF J I NAN UNI VERSITY (Sci.&T ech 1) V ol.18 N o.2 Jun.2004
中国科学院流固耦合系统力学重点实验室
中国科学院流固耦合系统力学 重点实验室 Key Laboratory for Mechanics in Fluid Solid Coupling Systems Institute of Mechanics, Chinese Academy of Sciences 季报 2019年第1期（总第17期） 目录 中科院流固耦合系统力学重点实验室现场评估工作顺利完成 (2) 中科院流固耦合系统力学重点实验室召开2019年室务会 (3) 中国航空学会空气动力学分会飞行载荷专业工作会在扬州召开 (6) 圆柱阵列波浪力幅值的波动现象和预报公式 (8) 轻质金属点阵圆柱壳结构制备与力学性能研究进展 (9) 力学所提出一种大幅提升3D打印点阵结构力学性能的新方法 (11) 雾化稠油掺稀降粘技术研究进展 (12) 南海天然气水合物试采安全评价研究进展 (14) 油气水多相流量计研究进展 (15) 空化致板间液滴界面稳定性研究获得多个奖项 (16) 空泡与柔性膜的流固耦合研究获得2019度中国力学大会优秀墙报奖. 18
中科院流固耦合系统力学重点实验室现场评估工作顺利完成 7月15日，中科院前沿科学与教育局、中科院重点实验室现场评估专家组一行14人莅临中科院力学所，对依托力学所建设的流固耦合系统力学重点实验室进行现场评估。专家组组长顾逸东院士主持了评估会议并宣布了现场评估的议程安排。力学所所长秦伟，党委书记、副所长刘桂菊，副所长魏宇杰，副所长尹明及流固耦合系统力学重点实验室学术委员会主任、实验室主任参加会议。 实验室主任黄晨光做实验室主任工作报告，围绕发展定位与研究方向、科研任务与代表性成果、队伍建设与人才培养、开放交流与运行管理等方面，向专家组汇报了评估期内的发展成果和工作成效。杨国伟研究员、王展研究员分别做“高速列车气动设计与流固耦合动力学特性研究”和“极端海洋环境及其与工程结构的流固耦合理论”代表性成果报告。专家组肯定了实验室取得的成绩以及工作亮点，并就汇报和自评估报告中的存疑事项进行了交流。 现场评估专家组还查看了高速列车动模型试验平台、海洋流固土耦合实验室、多相流体力学实验室、冲击与耦合效应实验室的科研仪器建设、大型科研仪器设备使用共享等情况，同时，参观了实验室的展板窗口。在此基础上，专家组召开会议，根据现场考核情况对实验室进行打分，并初步形成了评估意见。 经过努力，实验室顺利完成了此次中科院重点实验室现场评估工作，并在评估中充分展现了自身的优势和特色，最终取得良好的评估成绩。 在国家科技创新基地优化整合的背景下，实验室将积极适应新形势和新要求，进一步加强实验室建设和运行管理工作，全面提升科研平台建设水平和运行效率，为加快科技创新提供良好的条件支撑。 （流固耦合系统力学重点实验室供稿）
流固耦合手册现场翻译
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 workbench 流固耦合计算实例
Oscillating Plate with Two-Way Fluid-Structure Interaction Introduction This tutorial includes: ?Features ?Overview of the Problem to Solve ?Setting up the Solid Physics in Simulation (ANSYS Workbench) ?Setting up the Fluid Physics and ANSYS Multi-field Settings in ANSYS CFX-Pre ?Obtaining a Solution using ANSYS CFX-Solver Manager ?Viewing Results in ANSYS CFX-Post If this is the first tutorial you are working with, it is important to review the following topics before beginning: ?Setting the Working Directory ?Changing the Display Colors Unless you plan on running a session file, you should copy the sample files used in this tutorial from the installation folder for your software (
流固耦合
流固耦合定义：它是研究变形固体在流场作用下的各种行为以及固体位形对流场影响这二者相互作用的一门科学。流固耦合力学的重要特征是两相介质之间的相互作用,变形固体在流体载荷作用下会产生变形或运动。变形或运动又反过来影响流，从而改变流体载荷的分布和大小，正是这种相互作用将在不同条件下产生形形色色的流固耦合现象。 (一)流固耦合动力学：求解方法与基本理论---张阿漫，戴绍仕 ●有限元法 ●边界元法 ●SPH法与谱单元法 ●瞬态载荷作用下流固耦合分析方法 ●小尺度物体的流固耦合振动 ●水下气泡与边界的耦合效应 按耦合机理分两大类： 1 耦合作用只发生在两相交界面---界面耦合（场间不相互重叠与渗透），耦合作用通过界面力（包括多相流的相间作用力等）起作用。它的计算只要满足耦合界面力平衡，界面相容就可以了（其耦合效应是通过在方程中引入两相耦合面边界条件的平衡及协调关系来实现的）。如气动弹性，水动弹性等。 按照两相间相对运动的大小及相互作用分为三类： (1)流体和固体结构之间有大的相对运动问题"最典型的例子是飞机机翼颤振和 悬索桥振荡中存在的气固相互作用问题,一般习惯称为气动弹性力学问题" (2)具有流体有限位移的短期问题"这类问题由引起位形变化的流体中的爆炸或 冲击引起"其特点是:我们极其关心的相互作用是在瞬间完成的,总位移是有限的,但 流体的压缩性是十分重要的" (3)具有流体有限位移的长期问题"如近海结构对波或地震的响应!噪声振动的 响应!充液容器的液固耦合振动!船水响应等都是这类问题的典型例子"对这类问题, 主要关心的是耦合系统对外加动力荷载的动态响应" 2 两域部分或全部重叠在一起，难以明显的分开，使描述物理现象的方程，特别是本构方程需要针对具体的物理现象来建立，其耦合效应应通过建立与不同单相介质的本构方程等微分方程来体现。 按耦合求解方法分两大类： 1 直接耦合求解：直接耦合是在一个求解器中同时求解不同物理场的所有变量，需要针对具体的物理现象来建立本构方程，其耦合效应通过描述问题的微分方程来体现。 2 间接耦合求解：而间接耦合不需要重写本构方程，仅只利用当前比较成熟的单物理场求解器求解各自相域，并实现不同的物理场之间的信息交换。 范例（一个经典的间接耦合求解范例步骤）：利用CFX 进行全三维非定常粘性数值模拟，利用ANSYS 进行结构瞬态动力分析，其耦合面数据交换以MFX-ANSYS/CFX为平台，在每个物理时间步上进行耦合迭代，各自收敛后再瞬态向前推进，结构变形引起的流场网格位移由CFX内部的动网格技术来处理，整个耦合过程充分考虑了流场的三维非定常性和结构响应的瞬态变化。https://www.360docs.net/doc/5e19056934.html,/s/blog_6817db3a0100ju4s.html） 迭代求解,也就是在流场,结构上分别求解,在各个时间步之间耦合迭代,收敛后再向前推进.好
流固耦合概述及应用研究进展
流固耦合概述及应用研究进展 摘要 流固耦合力学是流体力学与固体力学交叉而生成的一门力学分支。顾名思义，它是研究变形固体在流场作用下的各种行为以及固体位形对流场影响这二者交互作用的一门科学。流固耦合力学的重要特征是两相介质之间的交互作用(fluid．solid interaction)：变形固体在流体载荷作用下会产生变形或运动，而变形或运动又反过来影响流场，从而改变流体载荷的分布和大小。总体上 , 流固耦合问题按耦合机理可分为两大类:一类的特征是流固耦合作用仅仅发生在流、固两相交界面上 ,在方程上耦合是由两相耦合面的平衡及协调关系引入的;另一类的特征是流、固两相部分或全部重叠在一起 ,耦合效用通过描述问题的微分方程来实现。 1 流固耦合概述 1．1引言 历史上，人们对流固耦合现象的早期认识源于飞机工程中的气动弹性问题。Wright兄弟和其它航空先驱者都曾遇到过气动弹性问题。直到1939年二战前夕，由于飞机工业的迅猛发展，大量出现的飞机气动弹性问题的需要，有一大批科学家和工程师投入这一问题的研究。从而，气动弹性力学开始发展成为一门独立的力学分支。如果将与飞机颤振密切相关的气动弹性研究作为流固耦合的第一次高潮的话，则与风激振动及化工容器密切相关的研究可作为流固耦合研究的第二次高潮。 事实上，从美国ASME应用力学部召开的历次流固耦合研讨会上可以看出，流固耦合问题涉及到很多方面。比如：空中爆炸及响应，噪声相互作用问题，气动弹性，水弹性问题，充液结构内的爆炸分析，管道中的水锤效应，充液容器的晃动及毛细流中血细胞的变形，沉浸结构的瞬态运动，流固相互冲击，板的颤振及流体引起的振动，圆柱由于热交换引起支持附件松动的非线性流固耦合系统，声音与结构的相互作用，涡流与结构的相互作用，机械工程中的机械气动弹性问题等等。 1.2流固耦合力学定义和特点 流固耦合力学是流体力学与固体力学交叉而生成的--I'l力学分支。顾名思义，它是研究变形固体在流场作用下的各种行为以及固体位形对流场影响这二者交互作用的一门科学。流固耦合力学的重要特征是两相介质之间的交互作用（fluid-solid interaction)．变形固体在流体载荷作用下会产生变形或运动，而变形或运动又反过来影响流场，从而改变流体载荷的分布和大小。正是这种相互作用将在不同条件下产生形形色色的流固耦合现象。流固耦合问题可由其耦合方程来定义，这组方程的定义域同时有流体域与固体域，而未知变量含有描述流体现象的变量及描述固体现象的变量，一般而言，具有以下两点特征： a)流体域或固体域均不可能单独地求解； b)无法显式地消去描述流体运动的独立变量或描述固体运动的独立变量。 1．3流固耦合力学涉及领域及分类 流固耦合问题涉及到很多方面。比如：工程实际中所涉及到的流固耦合问题，
(完整版)5流固耦合
第五章 轴流泵的流固耦合 5-1 流固耦合概论 流固耦合问题一般分为两类，一类是流‐固单向耦合，一类是流‐固双向耦合。单向耦合 应用于流场对固体作用后，固体变形不大，即流场的边界形貌改变很小，不影响流场分布的， 可以使用流固单向耦合。先计算出流场分布，然后将其中的关键参数作为载荷加载到固体结 构上。典型应用比如小型飞机按刚性体设计的机翼，机翼有明显的应力受载，但是形变很小， 对绕流不产生影响。当固体结构变形比较大，导致流场的边界形貌发生改变后，流场分布会 有明显变化时，单向耦合显然是不合适的，因此需要考虑固体变形对流场的影响，即双向耦 合。比如大型客机的机翼，上下跳动量可以达到5 米，以及一切机翼的气动弹性问题，都是 因为两者相互影响产生的。因此在解决这类问题时，需要进行流固双向耦合计算。下面简单 介绍其理论基础。 连续流体介质运动是由经典力学和动力学控制的，在固定产考坐标系下，它们可以被表 达为质量、动量守恒形式: ()0v t ρρ?+??=? (1) ()B v vv f t ρρτ?+??-=? (2) 式中，ρ为流体密度；v 为速度向量；B f 流体介质的体力向量；τ为应力张量；在旋 转的参考坐标系下，控制方程变为： ()0r v v t ρρ?+??=? （3） (-)+B r r c v v v f f t ρρτ?+??=? （4） 形式和固定坐标系下基本相同，只是速度变成了相对速度，另外就是增加了附加力项 c f 。 固体有限元动力控制方程为： []{}[]{}{}...[]{}M u C u K u F ++= （5） 式中，[]M ，[]C ，[]K 分别是质量矩阵，阻尼矩阵以及刚度矩阵，{}F 为载荷矩阵。 流固耦合遵循最基本的守恒原则，所以在流固耦合交界面处，应满足流体与固体应力、 位移、热流量、温度等变量的相等或守恒，即满足如下四方程： f f s s n n ττ?=? (6) f s d d = (7) f s q q = (8) f s T T = (9) 5-2 单向流固耦合
流固耦合ansys命令流
new config fluid title 基于流固耦合作用下的双龙富水隧道稳定性研宄 set fluid off set log on set logfile yang 1 .log genzonradcylpOOOOpl 9.00 0p2 0 50 0 p3 0 0 8 size4 20 64 dim6 5 6 5 rat 1 1 I 1 group 围岩 gen zon cshell pOOOOpl 6.0 0 0 p2 0 50 0 p3 0 0 5.0 size 4 2064 dim 5.6 4.6 5.6 4.6 rat 1 1 1 1 group初期支护 gen zon cshell pO 0 0 0 pi 5.6 0 0 p2 0 50 0 p3 0 0 4.6 size 4 2064 dim 5.0 4.0 5.0 4.0 rat 1 1 1 1 group 二次衬砲fill group原岩 gen zon radcyl pOOOOpl 0 0 ?8.0 p2 0 50 0 p3 9.0 0 0 size 4 20 64 dim 3 6 3 6 rat 1 1 1 1 group围岩2 gen zon cshell pOOO Opl 0 0 -3.0 p2 0 50 0 p3 6.0 0 0 size 4 2064 dim 2.6 5.6 2.6 5.6 rat 1 1 1 1 group仰拱初期支护 gen zon cshell pO 0 0 Opl 0 0 -2.6 p2 0 50 0 p3 5.6 0 0 size 4 2064 dim2 5 2 5 rat 1 1 1 1 group仰拱二次衬砲fill group仰拱原岩 gen zone reflect normal -10 0 gen zone radtun pO 0 0 0 pi 45 0 0 p2 0 50 0 p3 0 0 20 size 3 20 3 12 dim 9 8 9 8 rat 1 1 1 1.1 group围岩3 gen zon reflect dip 0 ori 0 0 0 range x09y0 50z8 20 gen zon reflect dip 0 ori 0 0 0 range x 9 45 y 0 50 z 0 20 gen zon reflect dip 90 dd 270 ori 0 0 0 range x09y0 50z8 20 gen zon reflect dip 90 dd 270 ori 0 0 0 range x 0 9 y 0 50 z -8 -20 gen zon reflect dip 90 dd 270 ori 0 0 0 range x 9 45 y 0 50 z -20 20 save shuitun一model.sav model fl_iso prop perm 1.23e-9 poro 0.45 range z 4.5 20 prop perm 4.70e-10 poro 0.4 range z -20 4.5 set fl biot off ini fdensity le3 ini sat 1.0 ini food 2.0e9 ftens -le-3 ;假设围岩岩体符合mohr-coulomb本构模型,给围岩陚参数命令流如下, ;mohr-coulomb model model mohr def derive s_modl=E一modl/(2.0*(l .0+p—ratio 1)) b_modl=E_modl/(3.0*(1.0-2.0*p_ratiol)) s_mod2=E—mod2/(2.0*(l .0+p_ratio2)) bjmod2~E_mod2/(3.0*( 1.0-2.0*p_ratio2)) end
流固耦合文献总结
小论文拟采用DP模型，在应力较高的土体中，比Mohr-coulomb理想弹塑性模型的数值计算结果更精确。设定DP模型需要输入3个特殊参数，粘聚力，内摩擦角，膨胀角，其中的膨胀角是用来控制体积膨胀的大小的。在岩土工程中，一般密实的砂土和超强固结土在发生剪切的时候会出现体积膨胀，因为颗粒重新排列了；而一般的砂土或者正常固结的土体，只会发生剪缩。在使用DP模型的时候，对于一般的土，膨胀角设置为0度比较符合实际。渗流耦合分析拟采用的边界条件是全地基边界，即把要分析的模型所有的区域看成是一个封闭的整体。在计算渗流应力耦合分析时，考虑基坑空间效应，建立三维实体模型，不仅考虑施工降水耦合，也考虑施工间歇变形耦合。最终通过支护结构桩和锚杆的变形以及基坑的变形，得出以下两条结论：（1）采用渗流应力耦合理论计算的基坑工程变形形态符合实际情况，随着基坑开挖深度增加，基坑变形规律也符合实际情况。（2）渗流应力耦合情况下基坑变形与不考虑渗流耦合影响下基坑变形曲线相比，数值较大，可见，分析基坑变形时不考虑渗流耦合影响是偏不安全的，耦合分析对基坑变形的影响不能忽视。 1、基于渗流场-应力场耦合作用下的深基坑降水支护结构的位移研究工程勘察2012 本文采用大型通用岩土工程有限元软件PLAXIS对复合土钉支护进行分析,模型采用平面应变模型,土体采用Mohr-coulomb理想弹塑性模型且具有对称性，故取一半对其分析，模型底部为固定约束，侧面只限制水平位移，上表面为自由边界。 本工程的数值模拟主要为比较在有降水作用下和未考虑地下水两种情况下的支护结构体系的位移，为此，首先进行了在未考虑地下水条件下的模拟，即不考虑孔隙水压，地下水位线默认为基坑底部。其次依据实际工程的地下水位线-7.24m，进行了数值模拟，以便找到降水作用对支护结构体系位移的影响。 2、考虑流-固耦合效应的基坑水土压力计算工程勘察2011 针对地下水绕过围护墙渗流情况，分析了传统的水土压力分算、合算及考虑土体渗流-固结变形方法计算土压力的区别，并利用实测数据进行对比。 流过耦合分析，PLAXIS程序采用水土分算的方法，通过输入地下水水头执行地下水渗流程序进行计算，利用单元应力点上的压力水头求得孔隙水压力，将围护墙与土体接触界面上的有效压力与孔隙水压力值相加，得到基坑围护墙上总的水土压力分布。 3、考虑流固耦合作用的深基坑有限元分析地下空间与工程学报2012 利用FLAC流固耦合模型对复杂地质条件下深基坑降水开挖过程中深基坑的时间效应进行研究。建立考虑参数变化的弹塑性流固耦合数值模型，分析基坑开挖及降水作用下地表沉降、水压力、基底隆起随时间变化的规律。平面应变模型，土体采用修正的剑桥模型模拟，只是在理论上提出考虑基坑开挖过程中渗透系数随孔隙比变化的现象，未应用在模型模拟中。 4、考虑渗流-应力耦合基坑开挖降水数值分析广东工业大学学报2013 本文运用通用软件MIDAS/GTS考虑渗流应力耦合作用下模拟基坑开挖降水的详细过程，分析了不同阶段渗流情况,同时探讨了止水帷幕、渗透系数与不同降水深度对基坑支护特性的影响，以期为基坑降水和支护结构优化提供理论参考。采用的摩尔库伦土体模型，基坑较小，应力水平较低，平面应变模型，未考虑基坑的空间效应。 5、深基坑工程降水与地面沉降耦合数值模拟研究中国市政工程2012 采用基坑降水与地面沉降耦合模型分析，四周边界取为定水头边界，其中，求解地下水问题简化为求解地下水在多孔介质中流动的问题，建立相适应的地下水三维非稳定渗流数学模型为 地面沉降模型为 方程的求解条件为： 利用建立的三维渗流沉降模型预测抽水减压期间对水位降深和区域沉降影响。计算结果
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)到工作区
流固耦合问题研究进展及展望
流固耦合问题研究进展及展望 摘要：天然岩体大多数为多相不连续介质，岩体内充满着诸如节理、裂隙、断层、接触带、剪切带等各种各样的不连续面，为地下水提供了储存和运动的场所。地下水的渗流以渗透应力作用于岩体，影响岩体中应力场的分布；同时岩体应力场的改变使裂隙产生变形，从而影响了裂隙的渗透性能，因此，流固耦合问题研究主要考虑流体在固体中的变化规律，尤其是流体渗流与和岩体应力之间的耦合作用,通过对国内外相关文献的分析与整理，从流固耦合的研究现状、特点、研究方法及展望这四个方面进行了论述。 关键词：流固耦合；岩体；地下水；研究方法；渗流 中图分类号：X523文献标识码：A 文章编号： 天然岩石不只是单一固相介质，尚有固相、液相和气相并存的多孔介质组合，岩石经历了漫长的成岩和改造历史，其内部富含各种缺陷，包括微裂纹、孔隙以及节理裂隙等宏观非连续面，它们的存在为地下水提供了储存和运动的场所。地下水的渗流还以渗透应力作用于岩体，影响岩体中应力场的分布，同时岩体应力场的改变往往使裂隙产生变形，影响裂隙的渗透性能，所以渗流场随着裂隙渗透性的变化重新分布，因此，在许多情况下必须考虑流体,包括液体(油或水)、气体(天然气、煤矿瓦斯等)在多孔介质中的流动规律及其对岩体本身的变形或强度造成的影响，即应考虑岩体内应力场与渗流场之间的相互耦合作用。 近年来，流固耦合问题越来越受到人们的重视，这方面的研究涉及许多领域，在采矿领域，涉及地热开发，石油开采中的流固耦合渗流，采矿围岩突水问题等。在建筑工程领域，包括地下水抽取引起的地面沉降问题，基坑渗流引起变形问题，坝基渗流及稳定性问题，隧道建设等。在环境工程领域涉及地下核废料存储，城市垃圾废弃物处理等以及生物医学工程等领域，这一问题的研究对促进科技进步和解决实际工程技术问题有着重要意义。 1 国内外研究现状 关于岩体和流体相互作用研究最早见诸K.Terzaghi对有关地面沉降研究，其内容主要限于考虑一维弹性孔隙介质中饱和流体流动时的固结，提出了著名的有效应力公式，迄今该公式仍是研究岩体和流体相互作用的基础公式之一。二十世纪中期Biot(1941,1956)进一步研究了三向变形材料与孔隙压力的相互作用，并在一些假设，如材料为各向同性、线弹性小变形，孔隙流体是不可压缩的且充满固体骨架的孔隙空间，而流体通过孔隙骨架的流动满足达西定律的基础上，建立了比较完善的三维固结理论。在此基础上，进一步发展了多相饱和渗流与孔隙介质耦合作用的理论模型，并在连续介质力学的系统框架内建立了多相流体运移和变形空隙介质耦合问题的理论模型。 Lous等（1974）运用单裂隙试件进行单向水流的室内模型，综合研究了天
流固耦合分析相关关键字
流固耦合分析相关关键字单元算法定义 *SECTION_SOLID *SECTION_SOLID_ALE *INITIAL_VOID_OPTIONS 多物质单元定义 *ALE_MUL TI_MA TERIAL_GROUP 多物质材料ALE网格控制 *ALE_REFERENCE_SYSTEM_CURVE *ALE_REFERENCE_SYSTEM_GROUP *ALE_REFERENCE_SYSTEM_NODE *ALE_REFERENCE_SYSTEM_SWITCH 流固耦合定义 *CONSTRAINED_LAGRANGE_IN_SOLID ALE算法控制 *CONTROL_ALE *ALE_SMOOTHING 材料本构及状态方程定义 *MA T_NULL(空气、水等材料) *MA T_V ACUUM *MA T_OPTION(结构材料) *EOS_OPTION（流体、结构材料的状态方程）
爆炸分析相关关键字 材料模型 *MA T_HIGH_EXPLOSIVE_BURN（炸药材料） *MA T_ELASTIC_PLASTIC_HYDRO（推进剂） *MA T_NULL（空气、水等材料） *MA T_OPTION（结构材料） 状态方程 *EOS_JWL（各种炸药） *EOS_IGNITION_AND_GROWTH_OF_REACTION_IN_HE（推进剂燃烧） *EOS_JWLB（各种炸药） *EOS_SACK_TUESDA Y（炸药材料） *EOS_OPTION（结构材料的状态方程） *EOS_LINEAR_POL YNOMIAL（空气） *EOS_GRUNEISEN（水、油等） 接触类型（Lagrange方法） *CONTACT_2D_AUTOMA TIC_SURFACE_TO_SURFACE *CONTACT_2D_SLIDING_ONL Y *CONTACT_SLIDING_ONL Y_OPTIONS *CONTACT_SURFACE_TO_SURFACE *CONTACT_ERODING_SURFACE_TO_SURFACE 起爆方式（单点、多点、线起爆）
adina热-流-固耦合建模过程
基于adina热-流-固耦合建模过程 热-流-固耦合作用是存在高度非线性的复杂耦合作用。有关这三场的耦合作用研究在地石油工程、热资源开发、地下核废料存储安全、采矿工程等很多领域有着非常重要的应用价值。由于研究对象的不同，热流固耦合模型的形式存在差异，建立符合实际问题的三场耦合模型十分困难，文中在国内外学者对三场耦合模型理论研究的进展状况的基础上，通过一个例子，介绍了用adina建立模型的过程。 1三场耦合理论模式介绍 在三场耦合尤其是三场耦合机制的研究过程中，人们根据各自对三场耦合的认识提出了不同的三场耦合作用模式。1995年前有关三场耦合作用模式的研究在场与场之间的联系关系上主要是以速度等变量为桥梁，如HART、Jing提出的作用模式，其中Jing主要描述的核储存库三场耦合模式，后来作用模式发展为主体为物理现象，它们之间的相互联系是以场作用或物理作用为桥梁的，如Guvanasen、柴军瑞的作用模式，前者同样以核废料储库库围岩三场耦合作用研究为主，后者为一般模式。 Jing等描述了核废料贮库围岩裂隙岩体中的热-液-力耦合过程，如图1所示。H art等提出了如图2所示的三场耦合作用模式。柴军瑞从岩体渗流-应力-温度三者两两之间的相互关系出发，建立了如图3的作用模式。图中：口渗透水流对岩体固相的力学作用，一般应用有效应力原理来反映;a’为应力引起裂隙岩体空隙率和渗透特性变化，目前有经验关系式(如Lours负指数关系式)和理论关系式(包括各种概化情况下和各种概化模型下的理论关系式)两大类表示方法;b为温度引起热应变(力)及与温度有关的岩体固相力学特性变化;b’为岩体固相力学变形引起热力学特性变化及 岩体固相内部热耗散;c为水流的热对流及与岩体固相的热交换;c’为温度势梯度引起水份运动及与温度有关的水特性变化。 图1裂隙岩体中的热液力耦合过程(据Jing等。1995年)
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 最终出网格
几个ansys流固耦合的例子
一般说来，ANSYS的流固耦合主要有4种方式： 1，sequential 这需要用户进行APDL编程进行流固耦合 sequentia指的是顺序耦合 以采用MpCCI为例,你可以利用ANSYS和一个第三方CFD产品执行流固耦合分析。在这个方法中，基于网格的平行代码耦合界面(MpCCI) 将ANSYS和CFD程序耦合起来。即使网格上存在差别，MpCCI也能够实现流固界面的数据转换。ANSYS CD中包含有MpCCI库和一个相关实例。关于该方法的详细信息，参见ANSYS Coupled-Field Analysis Guide中的Sequential Couplin 2，FSI solver 流固耦合的设置过程非常简单，推荐你使用这种方式 3，multi-field solver 这是FSI solver的扩展，你可以使用它实现流体，结构，热，电磁等的耦合 4，直接采用特殊的单元进行直接耦合，耦合计算直接发生在单元刚度矩阵 一个流固耦合的例子 length=2 width=3 height=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 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
三个流固耦合分析实例
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