基于流固耦合算法多孔介质的有限元建模分析
有限元语言在CFD与流固耦合领域的应用成果
有限元语言在CFD与流固耦合领域的应用成果
1.NACA0012跨音速机翼颤振流固耦合计算
利用有限元自动生成系统(有限元语言)计算机
翼振动的位移和速度，FLUENT软件计算机翼外
流场的速度和压力，两者之间的接口程序实现数
据交互，实现了流体和固体的耦合计算，求解了
标准翼型NACA0012。
从结果分析来看，两个软
件结合起来求解机翼振颤问题是有效、方便的，
并能够扩展至其它的流固耦合问题。
机翼颤振计算
2.热管内液气两相流计算
热管是当今最理想的传热元件之一。
由于其设计、制造简单，广泛应用于航空航天、能源、电子、核能和太阳能等工程领域。
热管的计算涉及不可压蒸汽流动计算，多孔介质内液体流动计算，两部分流场之间的衔接，汽/液界面的相变，速度场与温度场的耦合等复杂物理现象。
用有限元语言对先进的传热元件——热管内流动与传热作数值模拟显示了软件在处理复杂多场耦合问题上的便易及潜在优势。
热管结构示意图
3.二维圆柱绕流计算
丹麦气象局访问学者提出一种新的湍流模型，并采用算子分裂算法，把N-S扩展方程的物理过程分解为扩散和对流两个过程，分别计算这两个过程，以模拟圆柱绕流的普遍规律。
二维圆柱绕流计算。
与应变状态相关的岩体双重孔隙介质 流-固耦合的有限元计算
dU dt
+ S2
dP1 dt
=
f2
(9)
式中：各系数的表达式为
图 1 双重介质模型简图 Fig.1 Schematic representation of double porosity model
∫ K = BT DBdV
(10)
∫ L1 = −α1 BTmNdV1
(11)
根据有效应力定律、达西定律、流体-固体的质 量守恒定律，利用格林公式经过数学运算最终可得
(21)
dt
式中：A 为刚度矩阵，B 为质量矩阵，X 为解矢量。
采用隐式变量θ 方法，可得
X t+1 = θX t+1 + (1 − θ ) X t
(22)
dX = X t+1 − X t
(23)
将式(22)，(23)代入式(21)并与时间相乘得
[B + θΔtA]X t+1 = [B − (1 −θ )ΔtA]X t + F t+1 (24)
∂2 p2 ∂ xi∂ xi
= a22
∂ p2 ∂t
−α2
∂ 2ui ∂ t∂ xi
−
Γ ρf
(3)
Γ ρf
= γ ( p1 − p2 )
(4)
式中：λ，G 为拉梅常数，可由弹性模量和泊松比确 定；α1 ，α 2 为耦合系数，定义为单元代表体中孔隙、 裂隙各自的体积变形与总体积变形的比率；a11 ，a22 分别为孔隙、裂隙岩体的压缩性，可由孔隙基质孔
第 22 卷 第 10 期
吉小明等. 与应变状态相关的岩体双重孔隙介质流-固耦合的有限元计算
• 1637·
论的孔隙-裂隙岩体流-固直接耦合的数学模型及有 限元公式，给出了由岩体分类指标确定渗透系数与 应变状态相关的渗透系数的计算表达式，编制出相 应的二维有限元程序，最后给出了算例以验证其可 行性。
流固耦合作用下油气储层变形的有限元分析
Finite Element Analysis of Oil-Gas Reservoir Deformation under Fluid Structure InteractionA Dissertation Submitted toChangzhou UniversityByJia Hai-boPetroleum and Natrual Gas EngineeringDissertation Supervisor: Prof. Qing-jie ZhuApril，2015常州大学学位论文原创性声明本人郑重声明：所呈交的学位论文是本人在导师指导下独立进行的研究工作及取得的研究成果。
除文中已经注明引用的内容外，本论文不含任何其他个人或集体已经发表或撰写过的作品成果。
对本文的研究做出重要贡献的个人和集体，均已在论文中以明确方式标明。
本人已完全意识到本声明的法律结果由本人承担。
作者签名：签字日期：年月日学位论文版权使用授权的说明本学位论文作者完全了解常州大学有关保留、使用学位论文的规定，即：研究生在校攻读学位期间论文工作的知识产权单位属常州大学。
学校有权保留并向国家有关部门或机构送交论文的复印件和磁盘，允许论文被查阅和借阅。
学校可以公布学位论文的全部或部分内容，可以采用影印、缩印或扫描等复制手段保存、汇编本学位论文。
保密论文注释：本学位论文属于保密范围，在年解密后适用本授权书。
非保密论文注释：本学位论文不属于保密范围，适用本授权书。
学位论文作者签名：签字日期：年月日导师签名：签字日期：年月日中文摘要油气储层变形是引起油气井套管变形和破坏的重要影响因素，尤其是在稠油油藏热采过程中更加突出，所以不可避免是套管的损坏率会非常的高。
通过近年辽河油田锦州采油厂对大量热采油井进行的调查，发现新投产的锦607区块套损严重，在统计的套损井中发生在油气储层部位的套损情况最为严重，这对油田开发生产构成极大的威胁。
油气储层段的地层变形是引起套损情况在储层段多发的重要原因，通过研究发现多孔介质流固耦合作用与储层发生变形是甚为相关联的，所以通过运用多孔介质流固耦合作用机理来研究注汽热采中储层的变形是相当必要和有研究价值的。
多孔介质流_固耦合渗流数学模型研究_李培超
第23卷第16期岩石力学与工程学报23(16)：2842～2842 2004年8月Chinese Journal of Rock Mechanics and Engineering Aug.，2004 多孔介质流－固耦合渗流数学模型研究＊ 李培超1，2(1长庆石油勘探局博士后科研工作站西安 710021) (2石油大学博士后流动站北京 102200)博士学位论文摘要有效应力是多孔介质力学中一个十分重要的基础概念之一，很多多孔介质力学问题的分析都建立在有效应力的概念之上，因此，如何正确地确定有效应力则是多孔介质有关力学理论成败的关键。
Terzaghi因其在土力学和地基处理工程领域中的卓越工作被称作土力学之父。
他的重要贡献之一就是第一个提出了有效应力的计算公式，或者说建立了多孔介质有效应力原理。
然而该有效应力原理在实际应用过程中出现了一些偏差，其原因可能在于该原理是在研究土力学时提出的，它可能仅适用于大孔隙、固体颗粒以点接触方式存在的土壤；而对于低孔隙度的土壤或致密岩石，颗粒之间的接触面积通常不能忽略时，仍需要在细观上考虑其结构性因素。
论文第1章在建立多孔介质的有效应力原理时提出，必须考虑其细观结构，并解析推导出了基于多孔介质之上的饱和多孔介质有效应力原理。
该原理包含了孔隙度φ这一表征多孔介质结构特性的重要参数；并体现了孔隙流体与固体骨架对总应力的分担作用，分担比例为φ∶(1－φ)。
第2章讲述的是该有效应力原理的一些工程应用研究，通过该研究取得如下成果：(1) 建立了上覆岩层压力、孔隙流体压力、基岩应力之间的解析关系式；(2) 推导出饱和土体孔隙度与孔隙水压之间的非线性关系式，该关系式与多孔介质的初始状态及初始孔隙水压力有关；(3) 建立了基于多孔介质的岩石剪切强度理论和岩石破裂压力计算公式。
第3章则将该有效应力原理引入流-固耦合渗流中，并根据平衡条件得出了应力场方程；充分分析流-固耦合渗流的物理特性，建立起孔隙度和渗透率动态模型；依据流体力学连续性方程，考虑流-固耦合情形下多孔介质骨架变形特性和流体的可压缩性，得到了孔隙流体的连续性方程。
基于有限元方法的流固耦合分析研究
基于有限元方法的流固耦合分析研究随着人们对材料和结构性能的需求不断提高，流固耦合问题研究越来越受到关注。
流固耦合分析主要研究流体流动与固体变形之间的相互作用和耦合效应。
在流固耦合分析中，有限元方法成为最常用的分析方法之一。
本文将从三个方面对基于有限元方法的流固耦合分析进行讨论。
一. 基本原理在流固耦合分析中，有限元方法起到了至关重要的作用。
在有限元方法应用中，流场模拟和固体变形分析的研究已经非常成熟。
但是，在将流场和固体变形结合起来进行研究时，需要考虑流场对固体进行的压力和剪切载荷以及固体对流体流动的影响。
这是一种非线性问题，需要进行高效的数值分析.有限元方法的基本原理是将运动的物体划分成有限个小部分，然后用有限元模型离散处理每个小部分，通过简单的微分方程组成了均衡方程，然后使用数值法求解。
这种方法是一种数值分析方法，主要用来解决固体强度、稳定性、弹性和塑性分析问题。
同时，它还能够被用于研究流体流动、热传递和电磁场问题的解决方案。
二. 动力学中的应用在动力学中，有限元方法是一种广泛应用的方法，可以帮助人们准确预测材料和结构的力学性能。
基于有限元方法的流固耦合分析使得我们能够更为准确地预测材料结构的变形。
通过对其物理和机械属性的模拟，我们能够更好地了解物体的反应和行为。
流固耦合分析可用于模拟一般材料和技术性成分的结构，如飞机飞行时飞行表面的变形，以及汽车发动机在运行时的振动和变形。
同时，它还可以用于研究核反应堆的材料和组件，以及火箭发动机的设计。
三. 工业中的应用有限元方法的工业应用很广泛。
对于自动化工业而言，这种方法可以减少费用和时间，同时提高产品的质量。
例如，有限元模拟可以使用软件来模拟产品的变形，如汽车的碰撞试验，振动，加速度等。
有限元方法还可以用于模拟钢铁、铝和塑料的加工，同时还可用于气动设计和水动力学分析。
总结：流固耦合分析是一项复杂的工程技术，而有限元方法则是一种解决流固耦合问题的重要方法。
多孔介质流固耦合模型
多孔介质流固耦合模型全文共四篇示例，供读者参考第一篇示例：多孔介质是一种具有孔隙结构的介质，其内部空间可以充满气体或液体。
在工程应用中，多孔介质广泛存在于地下水层、土壤、岩石等领域，其流体运动受到物理和化学环境的复杂影响，因此需要建立多孔介质流固耦合模型来研究其流动特性。
多孔介质流固耦合模型是将多孔介质内部的固相颗粒和流体相耦合起来，通过数学方程描述多孔介质内部的流体流动和固体颗粒的运动规律。
多孔介质流固耦合模型的构建涉及多孔介质内部流动的描述、固相颗粒的运动规律、力学和流体动力学方程的建立等多个方面。
在多孔介质内部的流动描述方面，常用的方法包括达西定律和布里渊方程。
达西定律是描述多孔介质内部流体流动速度与渗透率之间的关系的经典定律，其表示为：\[ v = -K \frac{∇p}{μ} \]v为流速，K为渗透率，p为压力，μ为粘度。
布里渊方程则是描述多孔介质内部流体流动速度随位置变化的变化规律，其表示为：ρ为流体密度，g为重力加速度，z为高度。
达西定律和布里渊方程可以结合使用，建立多孔介质内部流动的数学模型。
固相颗粒的运动规律是多孔介质流固耦合模型的另一个关键方面。
通常情况下，多孔介质内部的固相颗粒会受到流体作用力和固相颗粒之间的相互作用力的影响，其运动规律可以用牛顿第二定律和达西定律描述。
牛顿第二定律表示为：\[ F = ma \]F为作用力，m为质量，a为加速度。
达西定律描述了固相颗粒受到的流体作用力与流速之间的关系。
通过以上分析，我们可以看出，多孔介质流固耦合模型是描述多孔介质内部流动和固相颗粒运动的重要数学模型，其建立需要考虑多个方面的因素。
在实际工程应用中，多孔介质流固耦合模型可以用于研究地下水流动、土壤固结、岩石力学等问题，为工程设计和科学研究提供重要的参考依据。
希望通过本文的介绍，读者能够对多孔介质流固耦合模型有更深入的了解，并在实际应用中取得更好的研究成果。
第二篇示例：多孔介质是一种具有孔隙结构的材料，它可以允许流体通过或在其中存储流体。
流体固体动力耦合分析的有限元法
第22卷第6期 工 程 力 学 V ol.22 No.6 2005年 12 月ENGINEERING MECHANICSDec. 2005———————————————收稿日期：2003-11-13；修改日期：2004-04-19 基金项目：国家自然科学基金资助(50209014)作者简介：*刘云贺(1968)，男，辽宁人，副教授，博士，中国水利水电科学院博士后，从事流固耦合及水工结构抗震研究(E-mail: liuyhe@)； 俞茂宏(1934)，男，浙江人，教授，博士生导师，从事工程力学和结构强度理论研究；陈厚群(1932)，男，江苏人，教授，中国工程院院士，中国水利学会副理事长，从事水工结构抗震研究.文章编号：1000-4750(2005)06-0001-06流体固体动力耦合分析的有限元法*刘云贺1,2，俞茂宏3，陈厚群2(1. 西安理工大学，西安 710048；2. 中国水利水电科学研究院，北京 100044；3. 西安交通大学，西安 710049)摘 要：应用有限元法探讨了流体、固体接触界面由无限接触点对组成，并以接触点对的瞬态接触内力作为待定变量的流体固体动力耦合模型的数值求解方法。
分析了流体、固体域插值函数的特点，用二维八节点等参元及流固接触面上的接触点对单元，对流固耦合系统进行了离散化处理；并采用变分原理推导了反映流体固体动力相互作用机理的接触约束矩阵(或称动力耦合矩阵)，建立了有限元控制方程，给出了完整的数值计算方法, 研编了动力耦合系统的分析程序。
数值计算结果与经典理论解误差很小，验证了动力耦合模型和有限元求解方法的正确性及其较高的计算精度。
关键词：流固耦合；数值分析；接触约束矩阵；瞬态动力；变分原理 中途分类号：TV312 文献标识码：AFINITE ELEMENT METHOD FOR TRANSIENT ANALYSIS OFFLUID-STRUCTURE COUPLING PROBLEMLIU Yun-he 1,2, YU Mao-hong 3 , CHEN Hou-qun 2(1. Xi’an University of Technology, Xi’an 710048, China; 2. China Institute of Water Resources and Hydropower Research, Beijing 100044, China;3. Xi’an Jiaotong University, Xi’an 710049, China)Abstract: The transient analysis of fluid-structure coupling system is conducted using finite element method, in which the transient contact force of the unlimited contact points coupling in the contact interface between fluid and structure is taken as unknown variable. The 2-dimensional 8-node element and the contact point couple elements are used, based on the interpolating function feature in the fluid and structure domain. The fluid-structure coupling system is divided into elements, and the finite element equation with the contact constrain matrix that reflects the mechanism of interaction is derived from a variational principle. A computer program is developed for the transient analysis of fluid-solid coupling system. Numerical results are in excellent agreement with those of the available analytical solution.Key words: fluid-solid coupling; numerical analysis; contact constraint matrix; transient dynamic; variationalprinciple流体－固体动力耦合是广泛存在于水利、船舶、海洋以及航空航天等许多个工程领域中十分重要和多学科交叉的研究课题。
流固耦合动力学仿真方法及工程应用
流固耦合动力学仿真方法及工程应用流固耦合动力学仿真方法及工程应用包括:1. 流固耦合仿真方法:该方法通常使用流体力学和材料力学的基础原理来建模流固耦合现象。
将流体与固体材料紧密耦合,考虑流体的运动和材料表面的应力、应变和变形协调。
该方法的应用范围非常广泛,包括机械工程、航空航天工程、土木工程和水利工程等领域。
2. 有限元分析(FEA):有限元分析是一种计算机模拟方法,通过将整个系统划分为许多小部分,并在每个小部分中求解复杂的问题。
流固耦合动力学仿真通常涉及将流体和固体材料划分为许多单元,并对每个单元进行求解。
这种方法常用于机械工程、航空航天工程和土木工程等领域。
3. 数值模拟方法:数值模拟方法是一种通过计算机程序计算数值模型的方法。
流固耦合动力学仿真通常涉及对系统进行数值模拟,并使用计算机程序进行求解。
这种方法常用于机械工程、航空航天工程和土木工程等领域。
4. 基于物理模型的方法:基于物理模型的方法是一种将实际物理过程建模为数学模型的方法。
流固耦合动力学仿真通常涉及将实际物理过程建模为数学模型,并将其与材料力学和流体力学的基础原理进行耦合。
这种方法常用于机械工程、航空航天工程和水利工程等领域。
流固耦合动力学仿真方法及工程应用的主要方法如下:1. 有限元分析(FEA):该方法可以用于模拟流固耦合现象,将流体和固体材料划分为许多小部分,并对每个小部分进行求解。
2. 数值模拟方法:该方法可以用于模拟流固耦合现象,使用计算机程序进行求解。
3. 基于物理模型的方法:该方法可以用于将实际物理过程建模为数学模型,并将其与材料力学和流体力学的基础原理进行耦合。
4. 混合方法:混合方法是一种结合多种方法的方法,将不同方法结合起来,以获得更准确的结果。
多孔介质热流固耦合有限元分析
２ 年 ２月 ０ １ １
Ｆｅ ． ｂ２０１ １
文 章 编 号 ：１ ７－２ ２２ ｌ）ｌ０ １－６ ６４ｏ ６ （Ｏ Ｏ 一０ ２０ １
多孑 介质 热流 固耦 合 有 限元分 析 Ｌ
杨
（ 河北 理工 大学
通 ，刘
峰
唐 山 ０ ３０ ） ６ ０９
建 筑工程 学 院 ，河 北
关 键词 ：多孔介质 ；Ａ Ｎ ＤＩ Ａ；热 力采 油 ；耦 合
摘
要 ：探 讨 了应 用有 限元分析 软件 ＡＤ ＮＡ 进 行 多孔介 质 热流 固耦合 的分 析方 法。 以石 Ｉ
油生产 中的热力 采油为例 ，对 多孔介质 热 流 固耦合 问题 进行 数值模 拟。分析 了在 热力采 油
பைடு நூலகம்
过程 中 由于注蒸 汽以及流 体流 动等 因素对储 油层 的 影响 。对开采 选址 与提 高石 油产 量有 一
定 的参 考意 义。 中图分 类号 ：Ｔ ３ ２ Ｅ １ 文献标 志码 ：Ａ
０ 引言
在 热流 固耦合 问题研 究领域 ，最早 是对 地热 资源 的开发 研究 ，探讨地 热 区域 内地应 力 、温度变 化和 岩 石 渗透率 变化 的影 响 ；近年来 ， 内学 者研究 了非等 温条件 下煤层 中瓦斯渗流 的规 律 ， 国 并建 立 了数学 模型 ，
于研究 热流 固耦 合有很 大启发 ，流 固耦合 理论 中 的岩石 变形方 程可 以通过 改变本 构关 系推广 到热 流耦合 理 论 中的岩石 变形方 程等 。 Ｊ
１ 热 流 固耦合 的两种模 型
对 于热 流 固耦 合理 论 目前 还处 于探索 阶段 ，所 以对 该理 论 的深入研究 主要 以流 固耦合 理论 为平 台 ，借 鉴其 中的一些 思想 和方程 。现 在对该 问题 的研究 主要 有两 个方 向 ，并 提 出了两类模 型 ： Ｊ 第一类 称 为非 完全耦 合模 型 ，主要 是 以流 固耦 合理 论 为基础 进行研究 。在此 模 型中 ，认 为骨架 温度 的
基于多孔介质模型的污水过滤器有限
Value Engineering0引言在工业污水处理中，过滤是最为重要的环节，过滤环节的本质在于借助过滤介质的筛选和截留功能，对二级出水做进一步的去质，使其达到所要求的水质再进行排放[1]。
现代过滤理论研究认为，在过滤过程中，悬浮颗粒的去除，主要是由颗粒与滤料之间以及颗粒与颗粒之间的吸附（粘附）作用。
其中引起颗粒迁移的作用力有重力、惯性力或水流紊动等，对应可产生拦截、沉淀、惯性、扩散和水动力等多种过程[2]。
Ohio 国立大学的Young Woon Kang [3]采用砂滤池处理高浓度有机废水，结果表明能去90%以上的固体悬浮颗粒（SS ）。
Odegaard 等人[4]用粗填料处理城市污水，在滤速为20~30m/h 的情况下SS 的去除率能达到80%~90%。
通常的工程污水过滤的参数依靠人工经验。
随着有限元等仿真软件的发展，可事先通过模拟过滤相关参数，从而获得最优的设计结果。
中交第三航务工程局有限公司厦门分公司的薛宏伟[5]、赖银波[6]等工程师借助有限元软件构建了大直径污水排海顶管基坑的三维数值模型，并动态模拟了顶管施工的力学效应。
李明辉等人[7]基于ABAQUS 软件，研究了一种水平井多簇压裂的有限元模拟方法，结果表明该模拟数值结果与全耦合裂缝扩展模型结果基本相符。
董刚[8]使用间断有限元法计算油藏渗流问题，通过数值直接计算得到压力梯度，给出相应的压力分布图，饱和度分布图。
综上所述，目前，国内外对于常用的过滤方法、机理、滤料的种类进行了分析，提出其中的一些关键性参数，建立了一些简单的模型。
但是，污水整个过滤过程是非常复杂、非稳态、非线性的，基于简单的模型或者经验的设计结果是不精准的。
基于此，本文基于本研究基于有限元软件Workbench 的CFX 流体动力学模块，以固体悬浮颗粒（SS ）为指标，模拟污水在过滤器中的动态过程，研究和获得固定过滤器参数下的最大过滤流量。
1污水过滤器有限元分析过程本文首先建立了过滤器的三维模型，并结合软件Workbench 中的CFX 流体动力学模块，运用有限元法对过滤器模型以及过滤过程中的流体动力学进行分析，对过滤器进行大量运算以弥补人工计算速度和精度的不足。
流-固耦合问题边界元-有限元耦合方法分析
流-固耦合问题边界元-有限元耦合方法分析
张效松;叶天麒
【期刊名称】《石家庄铁道大学学报（自然科学版）》
【年(卷),期】2000(013)002
【摘要】利用边界元-有限元耦合方法对流-固耦合振动问题进行了分析,假设流体控制方程为Laplace方程,利用非连续边界元对流体域进行离散,从而有效地解决了边界元分析中的"角点效应"问题.固体以平面梁为模型,采用有限元进行离散,对非连续边界元和有限元的耦合问题进行了分析,通过对悬臂梁在一侧受液体作用时的瞬态响应分析的数值解同解析解的比较,表明了本文所给方法的有效性,同时为利用边界元-有限元耦合方法对流-固耦合分析奠定了基础.
【总页数】4页(P6-9)
【作者】张效松;叶天麒
【作者单位】石家庄铁道学院基础部石家庄 050043;西北工业大学飞机系西安710072
【正文语种】中文
【中图分类】TB115
【相关文献】
1.Oseen子边界元法在粘性流－固耦合运动中的应用 [J], 司马灿;张效慈;吴有生
2.边界元——有限元耦合方法分析 [J], 张效松
3.分析流-固耦合问题的ALE有限元平衡迭代算法 [J], 刘力菱;易丽清;魏泳涛
4.流-固耦合问题的ALE有限元分析 [J], 魏泳涛;刘力菱;Philippe H.Geubelle
5.非连续边界元-有限元耦合方法分析 [J], 张效松;叶天麒
因版权原因，仅展示原文概要，查看原文内容请购买。
U形渡槽流固耦合有限元建模与响应谱分析
U形渡槽流固耦合有限元建模与响应谱分析U形渡槽流固耦合有限元建模与响应谱分析韩晓育1，2，王俊1，2（1.黄河水利职业技术学院，河南开封475004；2.小流域水利河南省高校工程技术研究中心，河南开封475004）摘要：渡槽结构的振动特性分析是进行结构抗震设计与安全评价的重要内容之一。
结合景泰川电力提灌二期工程的特点，考虑实际工程输水过程中水体与槽体间的耦合作用，采用ANSYS软件中融合的流固耦合模块建立了渡槽三维有限元模型，并对3种不同水位工况进行模态计算，得出相应的固有频率和主振型特征，再将正常水位工况下的计算结果与用HHT方法辨识的结果进行对比，误差范围为0.2%～4.4%，能够满足实际工程要求，说明有限元模型的适用性较强，同时对其主振型特征进行分析总结并提出合理的建议，最后利用反应谱理论对渡槽结构进行动力响应计算，得出不同工况下渡槽结构的位移、应力响应最值及变化规律，为渡槽结构的抗震设计、后续研究和该工程的加固改造提供理论依据。
关键词：渡槽；流固耦合；有限元建模；模态参数辨识；模态分析；响应谱分析近年来，随着国家对水利建设投资的不断加大，渡槽在各地水利工程建设中得到了广泛应用［1］，为满足工程结构的实际需求，加强渡槽结构的分析研究是十分必要的。
常用于渡槽结构理论计算的方法有结构力学法和有限元模拟法，进行结构动力计算时往往选用有限元模拟法，其能通过图像来方便、直观地反映出结构的振动特征。
目前，渡槽抗震设计计算研究中考虑更多的是槽体与水体之间的流固耦合作用［2－4］，且不同模型的计算结果差别很大，动力计算模型的选取还有待统一。
实际工程抗震设计中只有合理地选取流固耦合仿真模型，其动力响应结果才能作为抗震设计的依据。
为此，笔者以景泰川电力提灌二期工程（以下简称景电工程）中大型渡槽为研究对象，结合流固耦合理论，使用ANSYS有限元软件建立了槽体－水体－排架－基础－地基为一体的三维有限元动力模型，对渡槽结构进行了3种不同水位工况的模态计算，通过将正常水位工况下渡槽结构的自振频率与用HHT［5－7］方法得到的辨识结果进行对比、分析，证实了FSI系统模型的适用性，结果精度较高，可在实际工程抗震设计计算中推广应用；最后利用基于FSI系统的有限元模型对景电工程渡槽结构进行了动力响应谱分析，其计算结果可为该工程的后期改造、加固提供理论依据。
隧道开挖的双重介质水力耦合模型有限元分析
隧道开挖的双重介质水力耦合模型有限元分析
吉小明
【期刊名称】《岩土力学》
【年(卷),期】2006(0)S1
【摘要】基于混合物理论孔隙-裂隙岩体的双重孔隙介质水力耦合计算的微分方程,利用伽辽金有限元法提出了相应的有限元公式,并基于岩体分类指标(RQD,RMR)提出了与岩体应变状态相关的渗透系数计算公式。
编制了裂隙岩体双重介质流固耦合的2-D 有限元程序,给出的验证算例表明,该程序是合理和实用的。
同时将该程序用于隧道开挖的模拟计算,探讨渗流效应对开挖隧道围岩变形与渗流场的影响。
计算结果表明,在隧道设计中不考虑渗流的影响是偏于不安全的。
【总页数】6页(P305-310)
【关键词】双重孔隙介质;水力耦合;有限元法;隧道开挖
【作者】吉小明
【作者单位】广东工业大学建设学院
【正文语种】中文
【中图分类】U451.2
【相关文献】
1.考虑裂隙的几何-力学特性的双重孔隙介质水-应力耦合模型及其有限元分析 [J], 张玉军;张维庆
2.双重孔隙介质THMM耦合模型及其有限元分析 [J], 张玉军;张维庆
3.岩体水力学基础（六）——岩体渗流场与应力场耦合的双重介质模型 [J], 仵彦卿
4.三维双重孔隙-裂隙介质热-水-应力-迁移耦合模型及其有限元分析 [J], 张玉军;张维庆
因版权原因，仅展示原文概要，查看原文内容请购买。
动力固结流_固耦合模型的求解方法
.
李尧臣， 男， 教授， 现从事固体力学方面教学和科研工作。 !&,+ 年生，
!"#$%&’() *&#")(+&,- .,% ( .,"-/(+&,- %$&-.,%’$#$-+ 0%,’$/"%$ "*&-1 /2-(#&’ ’,-*,)&/(+&,- +$’3-&4"$
!" #$%&’()*+ , -./0 12)+ , 3.45 6(7*&(7$8
成的向量； ! 为总体结点自由度向量。在计算中采用 对时域的增量法， 时间步长取 %;%%’ <。 设圆柱形土体的底部和侧面是固定的，且不透水， 顶部除夯击部分外是自由表面， 夯锤底部为不透水界 面 （图 ’） 。多孔介质的总体积为 # ， 总表面面积为 .， 于是可得土体的边界条件为： 在位移已知的边界 3* 上 在应力已知的边界 3% 上 在压力已知的边界 3, 上 在速度已知的边界 37 上 式中 量在边界上的给定值。 （"）
5 &
*% 为位移分量； #%& 为应变张量分量； # %% 为体积应变； %%&
为应力张量分量； ) % 为体积力分量； , 为孔隙水压力； 1% 为水的体力分量； 7 % 为渗流速度； - %& 为渗透系 数 张 量 分量； / 为储存系数； $为 *+,- 有效应力系数； $ %&’( 为应 力应变关系张量； 4% 为边界上的已知力； # 6 为边界上 已知的法向流速。 $和 / 可用 ./012-,3 孔隙压力系数 排水和不排水情况下的泊松比 ’ 、 8、 ’* 和多孔介质的 弹性常数来表示。初始位移场和孔隙压力场为零， 对 于实际地基土中存在初始地应力和初始孔隙水 压 力 的问题， 可采用迭加法进行处理。
饱和多孔介质流固耦合渗流的数学模型
Ξ
收稿日期 : 2002212202 基金项目 : 国家自然科学基金资助项目 ( 10102020) ; 宁波 — 中国科学院高新技术种子资金项目 作者简介 : 李培超 ( 1976～) ,男 ,山东寿光人 ,博士研究生 。 © 1994-2010 China Academic Journal Electronic Publishing House. All rights reserved.
Mathematical modeling of flow in saturated porous media on account of fluid2structure coupling effect
L I Pei2chao , KON G Xiang2yan , L U De2tang
A 辑第 18 卷第 4 期 水 动 力 学 研 究 与 进 展 Se r . A , V ol . 18 , N o. 4 2003 年 7 月 J O U R N AL O F H YD RO D YN A M I CS J uly , 2003 文章编号 :100024874 ( 2003 ) 0420419208
420
水 动 力 学 研 究 与 进 展 2003 年第 4 期
1 引言
经典渗流力学一般假定流体流动的多孔介质
( 比如岩石 、 土壤等 ) 是完全刚性的 , 即在孔隙流体
三维固结理论开展的 , 只是假设的多孔介质的应力 应变本构关系有所不同 ( 比如有的假设多孔介质为 弹塑性 , 有的为粘弹性等等) 或者孔隙流体假设为多 相流体或单相流体的差别 ( 有的称为非饱和和多孔 介质或者饱和多孔介质) 。 综上所述 , 目前所应用的流固耦合渗流模型几 乎都是建立在 Te rza ghi 有效应力原理之上的 。而 Te rza ghi 有效应力原理有其不足之处 , 本文作者从 渗流力学的观点出发 , 在前人研究的基础上 , 推导出 了基于多孔介质的有效应力原理[ 9 ,10 ] 。该有效应 力原理包含了多孔介质的结构参数 — 孔隙度 φ, 从 而代替了以前诸多有效应力公式中广为使用的经验 参数 ( 各个模型的经验参数范围各异 , 取决于其适 用对 象 和 使 用 条 件 , 比 如 比 较 常 用 的 Bi ot 常 数 [ 5 ,6 ,7 ] 等) ; 而且该有效应力公式是解析的 , 理论上 是严格的 , 它可以适用于各种不同的多孔介质 , 而无 须区分不同的适用对象和特定条件 。 基于多孔介质的有效应力原理为 :
基于ADINA的多孔介质流固耦合分析
基于ADINA的多孔介质流固耦合分析朱庆杰;任瑞;贾海波;徐伟彬;张涛【摘要】基于流固耦合的研究理论,根据多孔介质(过滤板)中流体与结构之间的耦合关系,以有限元分析软件ADINA为平台,建立多孔介质(过滤板)的三维有限元模型,计算得到多孔介质(过滤板)在流固耦合作用下应力分布及变形情况.计算结果表明流固耦合作用对多孔介质(过滤板)变形破坏有重要的影响.针对计算结果,为过滤板的设计提出几点建议.【期刊名称】《常州大学学报（自然科学版）》【年(卷),期】2015(027)001【总页数】5页(P81-85)【关键词】多孔介质;流固耦合;有限元;ADINA【作者】朱庆杰;任瑞;贾海波;徐伟彬;张涛【作者单位】常州大学石油工程学院,江苏常州213016;常州大学石油工程学院,江苏常州213016;常州大学石油工程学院,江苏常州213016;常州大学石油工程学院,江苏常州213016;常州大学石油工程学院,江苏常州213016【正文语种】中文【中图分类】TE992.2多孔介质流固耦合作用存在于诸多工程问题中，如水利水电工程的渗流与控制、水库诱发地震、核废料深埋处理、煤与瓦斯突出、岩坡和坝基的稳定性、热能开发、地下水资源的利用等领域[1]。
在油田污水处理中的精细过滤装置中，由于过滤介质一般选用多孔介质材料，所以在含油污水处理过程中也会出现多孔介质流固耦合作用。
如精细过滤过程中的微孔过滤[2]，其过滤原理是在压力的推动下，料浆穿过过滤介质，并依靠这种过滤介质的拦截作用实现过滤。
其中，过滤介质多采用金属丝网烧结的微孔过滤材料,也可称之为多孔介质(过滤板)。
含油污水通过过滤介质(过滤板)必然会存在流固耦合作用，而流固耦合作用会影响过滤介质(过滤板)的使用寿命。
因此有必要采用有限元分析方法对过滤介质(过滤板)的工作过程进行数值模拟，研究过滤介质(过滤板)工作时的受力及变形情况。
国外学者已将有限元数值模拟方法应用到工程分析和设计中[3]。
多孔介质流_固_热三场全耦合数学模型及数值模拟_盛金昌
变形、能量传输、流体流动 3 个相互耦合的过程， 多物理场之间的交叉耦合还包括材料性质与独立变 量之间的耦合关系)，然后利用 FEMLAB 软件(基于 偏微分方程组的多物理过程模拟工具)作为平台，成 功地求解了这一全耦合偏微分方程组，避免了松散 耦合法求解多场耦合问题带来的误差，实现了同时 求解多物理场耦合过程。本文的求解方法是一个全 耦合算法，在理论上它能给出最真实的结果。通过 对一个具有已知解析解和数值解的算例的计算分析 来证明本文耦合模型及求解方法的正确性：一维砂 柱的等温固结和非等温固结问题。最后模拟分析了 通过井孔向岩体中注入冷水时流、固、热全耦合过 程，详细地分析了全耦合作用对井壁围岩应力的影 响。
·3030·
岩石力学与工程学报
2006 年
由于应变、流体压力和温度所引起的流体体积的变
化，等号左边的最后 1 项代表由压力梯度和重力作
用而引起的流体流量。
2.2 能量守恒方程 固体骨架和流体共同存在于同一个体积空间，
但它们具有不同的热动力学特性：如比热容和热传
导系数等。因此固体骨架和流体的能量守恒方程需
式中：φ 为岩体孔隙率， ρl 为流体的密度，t 为时 间，Vl 为流体速度矢量， Q 为流体的源汇项。
根据流体流动的动量方程可得 Darcy 定律：
Vl
=
−
k µl
(∇P
−
ρl g)
(2)
式中： µl 为流体的动力黏滞系数，k 为孔隙介质的 渗透率， P 为孔隙压力，g 为重力加速度矢量。
将式(2)代入式(1)，并加上固体骨架的变形项，
基于多孔介质流固耦合模型的地质灾害预测系统和方法[发明专利]
![基于多孔介质流固耦合模型的地质灾害预测系统和方法[发明专利]](https://img.taocdn.com/s3/m/20c05f6c814d2b160b4e767f5acfa1c7ab00824d.png)
(10)申请公布号 (43)申请公布日 2013.12.04C N 103425877 A (21)申请号 201310329385.1(22)申请日 2013.07.31G06F 19/00(2011.01)G06T 17/30(2006.01)(71)申请人湖北思高科技发展有限公司地址430070 湖北省武汉市东湖开发区创业街73号文化路武汉留学生创业园E 栋3176(72)发明人陈晓辉 付磊(74)专利代理机构北京轻创知识产权代理有限公司 11212代理人杨立(54)发明名称基于多孔介质流固耦合模型的地质灾害预测系统和方法(57)摘要本发明涉及基于多孔介质流固耦合模型的地质灾害预测系统和方法，包括以下步骤：采集待预测地理区域的地质参数；根据待预测地理区域的地理特征和地质参数，建立水文地质概念模型；根据要预测的地质灾害类型和待预测区域的地理特征，采用有限元法将所述水文地质概念模型转化为对应的三维地下水数学模型，并得到待预测地理区域的位移和应力随时空分布情况的模拟结果；采用图形和/或数据文字形式输出所述模拟结果，所述模拟结果用于描述地下水在空隙介质间的流动和地下土壤的沉降情况；预测地质灾害发生的时间、地理范围及形成机理，并给出解决方案。
本发明可以预测地陷、地基沉降、山体滑坡、泥石流等地质灾害，运算速度快且准确性高。
(51)Int.Cl.权利要求书2页 说明书7页 附图7页(19)中华人民共和国国家知识产权局(12)发明专利申请权利要求书2页 说明书7页 附图7页(10)申请公布号CN 103425877 A*CN103425877A*1.基于多孔介质流固耦合模型的地质灾害预测方法，包括以下步骤：a、采集待预测地理区域的地质参数；b、根据待预测地理区域的地理特征和所述地质参数，建立水文地质概念模型；c、根据要预测的地质灾害类型和待预测区域的地理特征，采用有限元法将所述水文地质概念模型转化为对应的三维地下水数学模型，并向所述数学模型输入所述地质参数得到待预测地理区域的位移和应力随时空分布情况的模拟结果；d、采用图形和/或数据文字形式输出所述模拟结果，所述模拟结果用于描述地下水在空隙介质间的流动和地下土壤的沉降情况；e、根据输出的模拟结果，预测地质灾害发生的时间、地理范围及形成机理，并给出解决方案。
基于流固耦合算法的多孔铜爆炸压实数值模拟
基于流固耦合算法的多孔铜爆炸压实数值模拟
张杰;赵铮;杜长星
【期刊名称】《科学技术与工程》
【年(卷),期】2010(010)032
【摘要】利用LS-DYNA有限元程序及流固耦合算法,对多孔铜的爆炸压实过程进行数值模拟研究.建立了多孔铜爆炸压实的三维四分之一模型,模型由空气、炸药和多孔铜三部分组成.炸药、空气与多孔铜的相互作用通过流固耦合算法来实现.根据得到的压力云图,对球面爆轰波和滑移爆轰的传播,多孔铜侧壁压力变化等问题进行了分析.结果表明,爆炸压实过程中多孔铜侧壁压力随滑移爆轰向下传播而逐渐增加,当夹角增加到90度后,角度及压力不再变化,形成稳定的滑移爆轰,多孔铜侧壁压力达到最大值,约为炸药爆压的85%.
【总页数】4页(P8001-8004)
【作者】张杰;赵铮;杜长星
【作者单位】中船重工713研究所,郑州,450052;南京理工大学能源与动力工程学院,南京,210094;南京理工大学能源与动力工程学院,南京,210094
【正文语种】中文
【中图分类】O347
【相关文献】
1.多孔覆盖层水下爆炸流固耦合分析 [J], 殷彩玉;金泽宇;谌勇;华宏星
2.爆炸冲击载荷作用下流固耦合数值模拟 [J], 郭攀;刘君;武文华;
3.基于流固耦合算法多孔介质的有限元建模分析 [J], 张杰;卜祥英;杨通
4.基于流固耦合的立式拱顶储罐油气爆炸数值模拟 [J], 欧益宏;李润;袁广强;李国庆;王世茂
5.变形多孔介质流固耦合模型及数值模拟研究 [J], 褚卫江;徐卫亚;苏静波
因版权原因，仅展示原文概要，查看原文内容请购买。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
震 、 坝稳 定 性 、 堤 煤层 瓦斯 渗 流等 ； 技 术上 与不 同工程 领域 如 土木 、 空 航天 、 从 航 船舶 、 动力 、 洋 、 海 石化 、 械 、 机
介 质在 流 固耦 合作 用 下的 变形 情况 。 中 图分 类 号 ： ＴＵ３ ２ ３ 文献标 志码 ： ９． Ａ
０ 引 言
传 统 的输 液 管道 系统 振动 分析 不考 虑流 体流 动 的多孔 介 质对 管 道 的影 响 ， 及认 为 多 孑 介质 骨 架 是 完 全 Ｌ 刚性 的 ， 在孔 隙 流体 压力 变化 过程 中 ， 固体骨 架不 产生 任 何弹 性或 塑 性 的变 形 ， 对输 液 管 道 振 动 分析 不 产 生 影 响 。这 是将 此 问题 作为 非耦 合 问题 来 研 究 ， 然 得 到 简化 ， 不切 合 实 际 ， 到 的 结 果 也是 不 精 确 的［ 。 虽 却 得 １ ］
１ ＡＤＩ 有 限元 建 模 ＮＡ
１ １ 建 模 方 法 介 绍 ．
ＡＤ ＮＡ 提供 了两 种建 模 的基本 方 法—— ｎ ｔ ｅ Ｓｍｐｅ 和 ｐ ｒｓ ｌ （ － ｐ 建模 Ｌ 。Ｎａｉｅ 模 是 一 Ｉ ａｉ （ ｉ ｌ） ａａｏｉ Ｂ Ｒｅ ） ｖ ｄ ４ ］ ｔ 建 ｖ 种 传统 的建模 方 法 ， 几何 元 素包 括 几 何 点 （ ｏｎ ） 几何 线 （ｉｅ 、 何 面 （ ｕｆｃ） 几 何 体 （ ｌｍｅ 等 。 其 ｐ ｉｔ 、 １ ）几 ｎ Ｓ ｒａｅ 、 Ｖｏｕ ） 每 个线 、 面和 体形 状 比较简 单 ， 能够 独立 进行 网格 剖分 。Ｎ ｔ ｅ 何 建模 方 式 是一 种从 底 向上 的建 模 方式 ， ａｉ 几 ｖ 由点 生成 线 ， 由线 生成 面 ， 由面 生成 体 。Ｐ ｒｓｌ （ － ｐ 建模 是 一 种 最 为广 泛 的 ３ 复 杂实 体 模 型 建 模 再 ａａｏｉ Ｂ Ｒｅ ） ｄ Ｄ 技术 ， 使用 边 界 曲面表 达法 描 述每个 实体 ， 个实 体都 典 型 的描述 了一个 要 进行 网格 划 分 的完 整 几 何 结构 。 每 其几何 元 素 为 Ｂ ｄ 、 ａ ｅＥｄ ｅ Ｐ ｉｔ其 中 Ｆ ｃ ｏ ｙ Ｆ ｃ 、 ｇ 、 ｏｎ ， ａｅ和 Ｅ ｇ 能作 为 Ｂ ｄ ｄ ｅ只 ｏ ｙ的组 成 部 分 存 在 ， 不 能 单 独 存 而
收稿 日期 ：０ ００ —８ ２ １ －９２
１６ ２
河北 理工 大学学报 （ 自然科拟 的 问题是水 道 中间有 块板 ， 为多孔 介质 属性 ， 可 以从 板渗透 过去 ； 水 由于 板两侧 有压 差存 在 ， 将 导致 板发 生变 形 。 １ ２ 建模 过 程 ． １ ２ １ 建立 结构模 型 ． ．
Ａｕ ． ０ ｇ ２ １ １
２１ ０ １年 ８月
文章 编 号 ：６ ４０ ６ （０ １０ —１ ５０ １ ７—２ ２ ２ １ ） ３０ ２ —４
基 于流 固耦 合 算 法 多孔 介 质 的有 限元 建 模 分 析
张杰 ， 卜祥 英 ， 通 杨
（ 北联合大学 河北省地震工程研究中心 ， 北 唐 山 ０３０） 河 河 ６ ０ ９
在 。Ｐ ｒｓｌ ａａｏｉ ｄ提供 了一系列 基本 几何 体 ， 过 布尔 操 作将 几 何 体 相加 减 ， 立 分 析需 要 的几 何 模 型 。两 种 通 建 方法 建 立 的几 何 模型 在几 何形 状和 尺寸 上完 全相 同 ， 但在 网格划 分 上有 所 不 同 。Ｎａｉｅ建 立 的模 型 在 划 分 ｔ ｖ 网格时 生成 的 网格质 量 比 ｐ ｒｓｌ ａａｏｉ 成 的 网格 质 量 要较 好 些 ， ｄ生 网格 更 加 规 则［ 。Ｐ ｒｓｌ ５ ］ ａａｏｉ 立 的 模 型 在 ｄ建 运算 时 容 易 出现 网格 扭 曲 ， 而导致 模 型不能 运行 。 从
翻
・
定 义几 何模 型 ： 采用 ｎ ｔ ｅ 模方 法 ， ｅｉ 建 ｖ 通过 点 、 、 的方式 建立 出多孔 介质 属性 的水 道 中间 板 ， 寸选 线 面 尺
第 ３卷 ３
第 ３期
河北 理工大 学 学报 （ 自然科学 版 ）
Ｊ ｕｎ ｌ ｆ ｂ ｉ ｏｙｅｈ ｉ Ｕ ｉ ｒｉ Ｎａｕ ａ Ｓ ｉｎ ｅＥ ｉ ｏ ） ｏ ｒ ａ ｏ ｅ Ｐ ｌｔ ｎｃ ｎｖ ｓ ｙ（ ｔ ｒｌ ｃ ｃ ｄｔ ｎ Ｈｅ ｃ ｅ ｔ ｅ ｉ
Ｖｏ． ３ Ｎ ． １３ ｏ ３
关 键词 ： 多孔介 质 ； 固耦 合 ； Ｉ 流一 ＡＤ ＮＡ
．
摘 要 ： 多孔介 质 流 固耦 合 的研 究 已被 广 泛 学者重视 ， 回顾 国内外 学者对 流 固耦 舍模 型 理论 在 进 展状 况 的基础 上 ， 主要 介 绍 了 多孔介 质 流 固耦 合 的基础 理论 ， 以有 限元 分析 软件 ＡＤ ＮＡ 为 Ｉ 平 台， 考虑 流 固耦合 分析 算 法 ， 建立 了 多孔介 质 的 三 维有 限 元模 型 ， 过 分析 计 算 得 到 了多孔 通
核 动力 、 地质 、 物工 程等 均有 关 系［ 。流 固耦 合 问题也 广泛 存在 于石 油工 业 中 ， 生 ２ ］ 如注 采过 程 中 的油 藏 渗 流 、 井 壁稳 定 、 井 和套 管 的损 坏 、 油 产层 出砂 、 水力 压 裂 、 层失 稳坍 塌等臼 。 地 ］ 流 固耦合 作 用对埋 地 管道 破坏 的影 响分 析 的研 究 日趋成 熟 ， 广泛 应 用工 程 实 际 ， 取得 了实 际 效 益 。 已 并 多 孔介 质 流 固耦 合 的研 究 也 吸引 了大批 工作 者逐 步深 入 ， 文 的工作 以有 限元 分析 软 件 ＡＤＩ 为平 台 ， 本 ＮＡ 介 绍 建模 方 法 ， 建立 多孔 介质 流 固耦合 有 限元模 型 。