ADINA流固耦合实例

合集下载

TG导入到adina进行流固耦合分析

TG导入到adina进行流固耦合分析

TG导入到ADINA进行流固耦合分析的基本步骤
Wry618
第一步:在tg中定义nastran的输出格式
第二步:建立相应的模型
在建立模型过程中,在建立体后,对自己想要建立约束、施加载荷或定义流固耦合面的地方需要定义相应的面。

由于经常定义的的约束、施加载荷的面或定义流固耦合的面不止一个,需要在tg中灵活应用mate、mti命令流。

建立完模型之后,输入命令流
merge
write
默认的文件名为truegrdo,将其拷贝到adina的工作文件目录下,这样便于查找(可以重命名)。

第三步:模型导入到ADINA
在ADINA file中import NASTRAN,导入truegrdo,
在导入时应注意,在creat boundary cells(BCELL)from shell element下选择replace。

第四步:在adina中定义约束和流固耦合界面
在adina模型树中就可以看见定义的shell和实体单元,此时的壳就是在tg中定义的面。

然后在mesh的nodeset或elmentset中定义相应的节点集或单元面集。

在method中选择from element group,在element group中选择相应的单元组,这样就能够完成节点集、单元集和单元面集的定义。

最后在施加相应的载荷,约束或流固耦合面,就可以得到
Tg入门相对较难,界面操作不是很方便,但是六面体网格划分能力确实令人佩服,在cnshijun兄的指导和单单的合作下,完成了tg到ADINA的转化的真正入门,谢谢他们。

希望大家以后多多交流,共同进步。

ADINA流固耦合建模方法

ADINA流固耦合建模方法

ADINA流固耦合建模方法ADINA流固耦合建模方法是一种综合考虑流体和固体相互作用的建模方法。

它结合了计算流体力学(CFD)和有限元力学(FEM),能够模拟和分析各种流体与固体相互作用的现象,如流体对结构的冲击、振动和与固体结构的热传导等。

1.定义流体区域:首先,需要在模型中定义流体的几何形状和流体域。

可以使用ADINA提供的几何建模工具或者导入已有的CAD文件来创建流体区域。

2.定义流体边界条件:在流体区域中定义流体的边界条件,如流体的入口速度、出口压力、壁面摩擦等。

可以通过给定边界条件来模拟各种流体流动情况。

3.网格划分:将流体区域划分为离散的网格单元,以便进行数值计算。

ADINA提供了自动划分网格的工具,也可以手动调整网格单元的大小和形状。

4.定义固体区域:在流体区域中定义固体的几何形状和固体域。

可以使用ADINA提供的几何建模工具或者导入已有的CAD文件来创建固体区域。

5.定义固体的边界条件:在固体区域中定义固体的边界条件,如固体的材料属性、固体的初始应力等。

根据具体问题,可以指定不同的边界条件。

6.载荷施加:在固体区域中施加外部载荷,如重力载荷、声压载荷等。

这些载荷将影响固体结构和流体流动的耦合过程。

7.运行求解器:通过ADINA的求解器对流固耦合建模进行求解。

求解器将同时考虑流体流动和固体结构的相互作用,求解固体受力平衡、流体流动动量方程等。

8.分析结果:根据求解结果,可以分析固体结构的变形、应力分布,以及流体流动的速度、压力等。

ADINA提供了丰富的结果分析工具,如绘制流线、应力云图等。

ADINA流固耦合建模方法能够模拟和分析多种流固耦合问题,如流体力学冲击载荷下的结构响应、流体流动对结构振动的影响、流体流动中的温度变化等。

它在航空航天、汽车工程、水利工程等领域具有广泛的应用。

通过ADINA流固耦合建模方法,可以提前发现和解决流固耦合问题,优化设计方案,提高产品的可靠性和性能。

流体流固耦合分析手册

流体流固耦合分析手册

5
第一章 计算流体力学概述
(3)
能量守恒方程。能量守恒属于经典的热力学定律。流体微团单位质量的能量 (由化学反应 即总能 E 包括内能 e 与动能 V ⋅ V 。合力所做的功、热传导、 等引起的)生成热都引起总能的变化。
1 2
1.3.1 预备知识 流体速度 u = (u1 , u2 , u3 ) ,u1 , u2 , u3 分别表示 x,y,z 方向上的速度分量。 ρ 是密度, p 是压强, T 是热力学温度。 若有过点 x = ( x1, x2 , x3 ) 的面积微元 dS ,单位法向量为 n 。 在 [t , t + dt ] 内沿 n 方向流过 dS 的流体体积为 u ⋅ ndSdt 在 [t , t + dt ] 内沿 n 方向流过 dS 的流体质量为 ρ u ⋅ ndSdt 在 [t , t + dt ] 内沿 n 方向流过 dS 的流体动量为 ρ u(u ⋅ n)dSdt =
下面用一个简单的例子来说明牛顿流体和非牛顿流体的差别。ADINA 的操作步骤参见 第二部分。
2
第一章 计算流体力学概述
例 1 分别考虑两个平行板间的牛顿流体和非牛顿流体,如下图所示,给出流体速度大小为 10m/s。
v
先考虑牛顿流体 流体参数为常参数模型,密度为 1 kg / m ,粘度为 0.04 N ⋅ s / m 2 。
流体流固耦合分析 手册
ADINA 中国代表处 Copyright © 2007
第一章 计算流体力学概述
第一章 计算流体力学概述
1.1 计算流体力学概述
计算流体力学(Computational Fluid Dynamics,简称 CFD)是以计算机作为模拟手段, 运用一定的计算技术寻求流体力学各种复杂问题的离散化数值解的计算方法。 计算流体力学可以看作是对基本守恒方程(质量守恒方程、动量守恒方程、能量守恒方 程)控制下的流动过程进行数值模拟。通过这种数值模拟,可以得到极其复杂问题的流场内 的各个位置上的基本物理量(如速度、压力、温度、浓度等)分布,以及这些物理量随时间 的变化情况。 计算流体力学、理论流体力学、实验流体力学是流体力学研究工作的三种主要手段。理 论分析具有普遍性, 各种影响因素清晰可见、 为实验和计算研究提供依据。 对于非线性情况, 只有少数问题能给出解析解。 实验研究仍是研究工作的基石, 数值研究的许多方面都密切依 赖于实验研究提供数据; 计算结果需由实验验证; 观察实验现象分析实验数据以建立计算模 型等等。数值模拟是特殊意义下的实验,也称数值实验,它比起实验研究,经济效益极为显 著。三种手段既互相独立又相辅相成。 近年来,由于实际工程设计对于流体计算提出越来越高的要求,计算流体力学在明显 地突破传统的单纯流体的观念, 各种涉及到复杂物理现象的流体问题求解方法是计算流体力 学发展的主要趋势,这些复杂现象是涉及热传递、多物质流动、相变、流固耦合体系求解、 变边界(变流动区域) 、湍流模拟等等。 从工程角度看,流体力学研究的起因通常是基于对各种工程结构的设计需要,例如分 析飞机机翼在气流作用下随机摆动问题的目的,是要求流场计算结果要对机翼的非稳态振 动、强度特性提出明确的力学设计指标。因此可以说,对于相当多的流体计算问题,实际上 我们需要知道的是一个耦合力学系统的响应特性-流固耦合体系特性, 尽管这是一个更为复 杂的计算体系。流固耦合(Fluid-Structure Interaction,简称 FSI)计算方法的开发和应用是 目前工程计算流体力学发展的重点领域,也是计算流体力学指导工程设计的直接途径。 计算流体力学的发展和计算机硬件求解能力、工程设计需求高速增长密切相关,可以 肯定地说计算流体力学在未来的研究领域和工程领域, 都会越来越走向实用化, 越来越发挥 不可或缺的作用。

基于ADINA的混流式水轮机流固耦合分析

基于ADINA的混流式水轮机流固耦合分析

合 界 面 ,建 立 了 混 流 式 水 轮 机数 值 模 型 。通 过 A I A 软件 平 台 ,采 用 流 固耦 合 模 块 ,实 现 了 转 轮 叶 片 与 内部 流场 耦 DN 合 的模 拟 ,得 到 了 转 轮 结 构 的应 力 应 变 、流场 的压 力 、流 速 和 k s分 布 。 - 关 键 词 :R G 湍 流 :流 固耦 合 :混 流式 水 轮 机 AN
(col f n r n ni n n, iu nvr t C e guS h a 10 9 S ho o eg adE vr metXh aU iesy hn d i un60 3) E y o i, c
Absr t t ac :The mate ai a d lo h F a i ur i e wa v lpe b u i g c n iuiy qu to h m tc lmo e ft e r ncs t b n s de eo d y sn o tn t e a in,m o n u e a in, me t m qu to
Re n l s sr s q a in n NG u b ln e 一 d la h o e nn q ai n .a d t e fud s u t r n e a t n y od te s e u t s a d RA o t r u e c s mo e s t e g v r i g e u to s n h l i —t cu e i tr c i r o me h d t t c h i tr c in s r c .B s d n t e t o o r e t e n e a t u f e a e o h ADI s f a e p af r a o a NA ot r l t m,t e f w f l o r n i u b n s w o h o ed f F a c s t r i e wa l i

基于ADINA的梭式止回阀流固耦合性能分析及优化

基于ADINA的梭式止回阀流固耦合性能分析及优化

( e a oa r f et gT c n l yf n fc r gPoes ns d ct n S h o o a ua tr g K y brt o s n eh oo r L o y T i g o Ma uat i rcs Mii r o E u ai ,co l f n fcu n un t f y o M i S i c n n ier g Suh et nvri f c n eadT c n l y Mi yn 2 0 0 C ia ce ea dE gn e n ,o t s U i s yo i c n eh o g , a a g6 1 1 , hn ) n i w e t Se o n
, 1 e ^ ● 、 . t § . § 、 t一 ● 、" § 1 喜 . 、. , — t 1 e斤 喜 \| 量 — " 、毫 e 斤 ’ 、 、 — 一、 | 重十 0 ” " 、 . ~
l 【 要】 摘 梭式止回阀作为管道非能动技术的关键组成部分 , 可以在无人力或外加动力驱动下 自 动实 i现管道系统的安全保障作用, 可以有效地降低输流管道中意外事故的发生。 围绕大口 径梭式止回阀的流固 ;耦合性能数值分析展开, DN 以A IA软件为平台, 建立梭式止回阀的流固耦合有限元模型, 对管道流场进行 l动态分析, 出了 给 数值模拟计算数据。 根据分析结果 , 对梭式止回阀的结构设计提 出了 优化方案 , 使梭式止 》回阀在管道流体中的运动流场更趋于平稳, 从而有利于提高梭式止回阀的使用寿命及安全性。 ; 关键词 : 梭式止回阀; DN ; A IA 流固耦合 ; 输流管道
O f ADI NA-b s d s ut p h c av a e h t e t e c e kv le l y

基于ADINA的隔膜流固耦合分析研究

基于ADINA的隔膜流固耦合分析研究

生产产量 。
【 参考文献】
[] 1 曹世海. 加热炉燃烧控制 系统 的几项改进措施叨. 自动化, 0 , ) 8 5 . 冶金 2 2( : — 0 0 44 [] 2张晶涛, 钱晓龙, 王伟, 张莉, 天佑. 柴 步进式加热 炉燃烧控制 的新 方法嗍. 制 控 与决策. 0 . 1 ( . 2 1 . 3 0 56) [] 3张元福加 热炉空, 比自寻优模糊控制器研究[. 燃 J山东冶金, 0 ,) 】 2 05. 0 (
压力低时。 保护步进梁, 为 煤气也能紧急 切断。 煤气切断 后自 动进行氮
在实际操作 中 , 当出现下列情况之一时 , 为了保证人员和设备 在 炉温 控制系统里 , 考虑 到控制功能 的连续性 , 由于烘炉 时温 度 气吹扫。 较低 , 升温时间短 , 常规 的温度控 制不稳定 , 反应 慢 , 在该项 目里增 加 的安全 , 燃烧控制 自动切 断或者 自 动锁定或继续。 另外 , 一个在操作 台 了程序 升温功能 . 通过设定升温的方式 、 升温 的温度 、 温的时间来 进 上会启动语音报警及声光报警 升 行程序升温控制 , 通过程序来 自动调节各 区域 的煤气流量 , 通过 温度 切断条件 : 检测来进行调整 。 控制功能如图 3所示 , 中,I : 其 TC炉温控制器 , FC G I: () 1系统跳电。 煤气 流量控 制器 。 FC空气流量控制器。 A I: () 2 动力气源压 力低 。 () 3 燃烧空气 压力低 。 温眉 () 4 煤气压力低。 T 2 () 5 冷却水压 力低。 () 6 操作急停。 其他报警信号: / ’ / 如果在操作中出现以下报警信号 . 以下的控制如燃烧控制 、 2 N 压 挂傩戕杰 , 动 1 库彳; 动 ' 年韵 自 陧 自 力控 制 、 压控 制 、 炉 换热 器人 E废 气温度控 制将继续 , l 报警 只在 C T R r 啬i f责 M M A 上显示 。 FT 冉 膏 r - M A f ) 2压力低 1N 青 M C () 2 燃烧空气温度高。 () 3 空气换热器上游温度高。 () 4 冷却水温度高。 缝 o O 干 确 v FV打 选择程序升温 ( ) 内温度高。 5炉 选择 FC启动 I () 6 煤气换热器下游的燃烧 空气温度高 。 图 3程 序 升 温 的控 制 po rmmi e t gc nrl rga n h ai o t g n o 3结 束语 . 2 其它控 制及安全连锁 . 4 通过对炉温的多种控制方式 的应用 .提高 了系统 的响应速度 . 使 241区域空气及煤气流量控制 .. 钢坯保温时间缩短 . 使氧化烧损 由原来的 1 %降低 燃烧控制 是建 立在煤气和空气 的流量控制上 的 . 两个 流量控制 炉 内温度控制均匀 , . 在一定程度上提高 了轧前水 除鳞效果 . 8 避免 了轧件 的打滑 现 是靠从 区域温度控制的输 出来 执行的。煤气和空气 的流速是 通过双 到 O %. 象, 减少 了废钢的几率 , 了生产 的顺行 , 保证 提高 了产 品成材率 和型钢 交叉 限幅燃 烧模式来维持空气, 煤气 比例值 在一个特 定的范 围 . 即使

adina热-流-固耦合建模过程

adina热-流-固耦合建模过程

基于adina热-流-固耦合建模过程热-流-固耦合作用是存在高度非线性的复杂耦合作用。

有关这三场的耦合作用研究在地石油工程、热资源开发、地下核废料存储安全、采矿工程等很多领域有着非常重要的应用价值。

由于研究对象的不同,热流固耦合模型的形式存在差异,建立符合实际问题的三场耦合模型十分困难,文中在国内外学者对三场耦合模型理论研究的进展状况的基础上,通过一个例子,介绍了用adina建立模型的过程。

1三场耦合理论模式介绍在三场耦合尤其是三场耦合机制的研究过程中,人们根据各自对三场耦合的认识提出了不同的三场耦合作用模式。

1995年前有关三场耦合作用模式的研究在场与场之间的联系关系上主要是以速度等变量为桥梁,如HART、Jing提出的作用模式,其中Jing主要描述的核储存库三场耦合模式,后来作用模式发展为主体为物理现象,它们之间的相互联系是以场作用或物理作用为桥梁的,如Guvanasen、柴军瑞的作用模式,前者同样以核废料储库库围岩三场耦合作用研究为主,后者为一般模式。

Jing等描述了核废料贮库围岩裂隙岩体中的热-液-力耦合过程,如图1所示。

H art等提出了如图2所示的三场耦合作用模式。

柴军瑞从岩体渗流-应力-温度三者两两之间的相互关系出发,建立了如图3的作用模式。

图中:口渗透水流对岩体固相的力学作用,一般应用有效应力原理来反映;a’为应力引起裂隙岩体空隙率和渗透特性变化,目前有经验关系式(如Lours负指数关系式)和理论关系式(包括各种概化情况下和各种概化模型下的理论关系式)两大类表示方法;b为温度引起热应变(力)及与温度有关的岩体固相力学特性变化;b’为岩体固相力学变形引起热力学特性变化及岩体固相内部热耗散;c为水流的热对流及与岩体固相的热交换;c’为温度势梯度引起水份运动及与温度有关的水特性变化。

图1裂隙岩体中的热液力耦合过程(据Jing等。

1995年)图2三场耦合模式(Hart)图3渗流-应力-温度之间的相互关系图2热流固耦合理论的提出三场耦合理论是由流-固两场耦合理论发展而来的,在流-固耦合理论中,有的假设温度场是恒定的,或者是不考虑温度场的变化与流体流动、岩石变形间的耦合作用。

ADINA在重力坝地震响应分析中的应用

ADINA在重力坝地震响应分析中的应用

连续性方程 运动方程
∂ui = 0
(2)
∂ xi
dui dt
+ cj
∂ui ∂x j
=
∂σ ij ∂x j
+
fi
(3)本构方程σFra bibliotekij=

p ρ
δ
ij
+ v( ∂ui ∂x j
+
∂u j ∂ xi
)
(4)
边界条件
⎧ ⎪

ui = ui
on SW
⎪⎨σ ijm j

=σi
on
Sf
(5)


⎪ p= p ⎩
2. ADINA 中 ALE 方法简介
ADINA 程序中采用 ALE 方法解决流固耦合问题。ALE 方法通过追踪自由液面的真实位
置,固体和液体之间的相互作用通过耦合面相互传递,实现流固耦合问题的分析.基本方程如
下:
ALE 描述下的水体导数可写为
f
= ∂f ∂t
x
+
∂f ∂ xi
ci
=
df dt
+ c∇f
图 3 水工建筑物抗震设计反应谱
在反应谱法计算中,设计加速度反应谱最大值 β max 选取:对重力坝可取 βmax = 2.0 。
一般来说,结构在地震作用下的响应中我们最关心的是位移和应力,而应力中又以最大主 应力和最小主应力最为重要。本文取坝体和坝基上若干关键点(见图)来考察最大主应力和最 小主应力,且只考察它们的最值,结果见表 2-3。取 A-I 关键点来考察位移, 坝体与坝基关键点 的示意图如下:
从表 1 可以看出,两种模型重力坝结构的自振频率的基频基本上一致的,差别相当小,说 明用附加质量模型代替流固耦合模型进行自振特性计算也比较精确。

adina15例子

adina15例子

橡胶支座的大变形分析
E7. 涡轮叶片的动力学分析 采用 ADINA 进行涡轮的频率和振形分 析,考虑涡轮高速旋转的离心力效应。
(挪威 Kvaemer Energy a.s.)
涡轮叶片的固有模态和振形
E8. 光导纤维束加工过程仿真 采用 ADINA 对多根纤维扭转 720 度而搓结 成束过程进行仿真,分析涉及到材料大变形和多 体接触摩擦等复杂因素。 (英国 Akzo Nobel 公司)
Ford 汽车 FMVSS 216 项标准规定的准静态试验
承载力分析结果和试验结果的比较
ADINA 的仿真结果
E2. 汽车 ABS 系统流动分析
ADINA 软 件 的 流 固 耦 合 ( Fluid Structure Interaction,简称 FSI)求解功能是真正意义上的多 场耦合计算。流固耦合计算中的一个突出的难点是移 动壁面的问题。通常,较小的结构变形可以采取 ALE 单元描述实现结构节点的移动并修正流场空间。 但是, 如果结构壁面移动范围很大, 仅靠 ALE 技术不能实现。 ADINA 采用 ALE 和网格自动重划分技术 (Automatic Remeshing)顺利解决这一难题。以下是汽车刹车防 抱死机构的示意图。刹车导致钢活塞突然发生上下移 动,液压油在机构内腔流动,并推动钢珠运动,与钢 珠相连的弹簧或松或紧控制刹车阻力。钢珠表面是流 场中的移动壁面,即随着钢珠的移动,流场空间在改 变。
美国芝加哥重点医疗咨询中心(MC3)采 用 ADINA 开发人工肺产品。人工肺的材料是 毛细纤维束,MC3 选用 ADINA 的多孔介质材 料模拟人工肺,外部覆盖一层结构材料。通过 外部结构运动使气体流入、流出人工肺。分析 得到人工肺中血压分布和血液流动特征。 (美国 MC3 研究中心)

ADINA流固耦合实例

ADINA流固耦合实例

实例3 隧道内具有柔性结构的流固耦合分析问题:隧道内具有柔性结构的流固耦合如图3-1所示。

图3-1 流体-固体结构示意图一、目的1. 掌握流固耦合作用FSI在Adina-AUI中的操作过程。

2. 掌握用伸缩比例因子画流固耦合模型。

3. 定义引导点(leader-follower points)。

二、定义模型主控数据1. 定义标题:选Control→Heading→敲入标题“exe03: Fluid flow over a flexible structure in a channel, ADINA input”→and click OK。

2. FSI分析:在右边Analysis Type区选FSI按钮。

3. 主控自由度选Control→Degrees of Freedom→不选X-Translation, X-Rotation, Y-Rotation andZ-Rotation按钮→and click OK。

4. 分析假设:大位移,小应变。

选Control→Analysis Assumptions→Kinematics→设置“Displacements/Rotations”为Large→ click OK。

(注:非常薄的结构,因此为小应变)。

三、力学模型1. 柔性结构建立模型1). 柔性结构几何模型坐标点如表3-1,几何结构如图3-2所示。

其几何面见表3-2所示。

①选Define Points 图标→按表3-1输入几何点坐标→ click OK .②选Define Surfaces 图标→设置TYPE 为Vertex → click OK(如图3-2所示)。

2). 施加固定边界条件和流-固边界条件①.图3-2中,在L2线上施加固定约束,其过程可用Adina-AUI 完成。

②. 流-固边界,选Model →Boundary Conditions →FSIBoundary →add FSI boundary number 1→在表中头两行敲入流固边界线编号1和 3 and click OK 。

ADINA压力容器多场耦合技术

ADINA压力容器多场耦合技术

第一部分 ADINA多场耦合求解功能介绍在很多实际工程问题中,我们需要知道的是流固耦合力学系统的整体响应特性,而不是单纯流体或者固体的响应特性,尽管系统响应特性与后者都具有相关性。

相对单一物理场求解,流固耦合面对一个更为复杂的计算体系,不仅需要求解固体问题和流体问题,还要求解固体和流体的相互作用,甚至热量在流体、固体之间进行传递的三场耦合。

ADINA-FSI 是全球领先的流固耦合求解器。

由于ADINA的结构求解器(ADINA-Sturctures)和流体求解器(ADINA-CFD)都是在同一内核基础上开发的,ADINA-FSI很容易将ADINA-Sturctures 和ADINA-CFD的功能完全地融合在一起,实现流体-固体耦合的高级分析。

流固耦合力学的重要特征是两相介质之间的相互作用,变形固体在流体载荷作用下会产生变形或运动,固体的变形或运动又反过来影响流体,从而改变流体载荷的分布和大小,正是这种相互作用将在不同条件下产生形形色色的流固耦合现象。

流固耦合问题可由其耦合方程定义,这组方程的定义域同时有流体域和固体域。

而未知变量含有描述流体现象的变量和含有描述固体现象的变量,一般而言具有以下两点特征:1)流体域与固体域均不可单独求解;2)无法显式地削去描述流体运动的独立变量及描述固体变形的独立变量一般的,我们概括ADINA流固耦合的特点如下:∙FCBI及FCBI-C单元提供了最大的稳定性,且适用于从低雷诺数到高雷诺数的各种问题。

∙ FSI分析可以用于各种流体,包括不可压流、轻微可压流、低速可压流和高速可压流。

另外,所有的流体材料模型包括非牛顿流体,湍流模型、相变模型和VOF法都可以用在FSI分析中。

∙ ADINA允许流体模型和结构模型使用任意的网格。

而且,流体和结构的网格在流固耦合界面上不必完全匹配。

∙热和多孔介质的耦合可以用在流体和结构模型之间。

∙所有的结构单元类型、接触算法和结构材料模型(如弹性、粘弹性、超弹、塑性、徐变或蠕变、复合材料、记忆合金、垫片材料等)都可以用于FSI求解。

液压悬置的流固耦合模拟

液压悬置的流固耦合模拟

液压悬置的流固耦合模拟
作者:ADINA
关键词:流固耦合,液压悬置,液力发动机悬置,强耦合,ALE 坐标,直接耦合求解器
这里我们演示汽车工业中用到的ADINA 流固耦合分析。

模拟的部件是图1所示的液压悬置(hydromount )。

分析的目的是计算液压悬置在不同频率下的刚度。

液压悬置的橡胶部分经历大变形以及它与金属悬置(metal mount )和底座(plate )之间的接触。

流体是粘性不可压缩的,用ALE 坐标描述。

这个问题用到了直接耦合的求解器,因为流体和结构之间是强耦合的。

对于液压悬置,首先做施加了预应力的静力分析,然后进行不同激励频率下的瞬态动力分析。

上边的动画显示的是两个不同激励频率10Hz 和400Hz (注意:动画显示的是相同的速度)下,橡胶上的最大主应力和流体速度。

图2显示了频率在10~1000Hz 之间,液压悬置的刚度变化。

图1 液压悬置
图2 液压悬置刚度曲线这个结果说明ADINA流固耦合技术有着广泛的应用。

基于ADINA的多孔介质流固耦合分析

基于ADINA的多孔介质流固耦合分析

基于ADINA的多孔介质流固耦合分析朱庆杰;任瑞;贾海波;徐伟彬;张涛【摘要】基于流固耦合的研究理论,根据多孔介质(过滤板)中流体与结构之间的耦合关系,以有限元分析软件ADINA为平台,建立多孔介质(过滤板)的三维有限元模型,计算得到多孔介质(过滤板)在流固耦合作用下应力分布及变形情况.计算结果表明流固耦合作用对多孔介质(过滤板)变形破坏有重要的影响.针对计算结果,为过滤板的设计提出几点建议.【期刊名称】《常州大学学报(自然科学版)》【年(卷),期】2015(027)001【总页数】5页(P81-85)【关键词】多孔介质;流固耦合;有限元;ADINA【作者】朱庆杰;任瑞;贾海波;徐伟彬;张涛【作者单位】常州大学石油工程学院,江苏常州213016;常州大学石油工程学院,江苏常州213016;常州大学石油工程学院,江苏常州213016;常州大学石油工程学院,江苏常州213016;常州大学石油工程学院,江苏常州213016【正文语种】中文【中图分类】TE992.2多孔介质流固耦合作用存在于诸多工程问题中,如水利水电工程的渗流与控制、水库诱发地震、核废料深埋处理、煤与瓦斯突出、岩坡和坝基的稳定性、热能开发、地下水资源的利用等领域[1]。

在油田污水处理中的精细过滤装置中,由于过滤介质一般选用多孔介质材料,所以在含油污水处理过程中也会出现多孔介质流固耦合作用。

如精细过滤过程中的微孔过滤[2],其过滤原理是在压力的推动下,料浆穿过过滤介质,并依靠这种过滤介质的拦截作用实现过滤。

其中,过滤介质多采用金属丝网烧结的微孔过滤材料,也可称之为多孔介质(过滤板)。

含油污水通过过滤介质(过滤板)必然会存在流固耦合作用,而流固耦合作用会影响过滤介质(过滤板)的使用寿命。

因此有必要采用有限元分析方法对过滤介质(过滤板)的工作过程进行数值模拟,研究过滤介质(过滤板)工作时的受力及变形情况。

国外学者已将有限元数值模拟方法应用到工程分析和设计中[3]。

ADINA流固耦合建模方法

ADINA流固耦合建模方法

ADINA流固耦合建模方法尽量采用几何模型做为最初的模型信息输入,而不是单元和网格;•所有的载荷、边界条件、初始条件施加在几何模型上,而不是节点或单元上面;将流固耦合等边界条件定义在单元上的工作量要远大于几何元素的情况;•结构模型和流体模型分布建立;结构模型可以包括ADINA Structure模块提供的所有特性;流体模型包括ADINA-Fluid所有特性;ADINA提供了各种流体专用网格生成方法;将结构模型和流体模型分别以dat文件的形式写出;同时将两个dat问题提交给ADINA-FSI求解器进行求解;•后处理中可以同时或者分开查看结果;原则:绝不能包含无关紧要的模型细节;•用正确的边界条件(或其它手段)补偿截断的流场空间;例如没有包括模拟障碍物后面的回流区,需要延长下游流场空间;如果不能包括足够大的流场空间,则需要采用瞬态方式进行求解;•当结构小到不能影响主要流场特性时,尽量抹除。

一个模型存在错误的原因可能非常多,此时将无法判断从什么方面进行模型的修正。

从节省时间的角度,应该按照下面的进行模型测试过程:1. 用的你经验或者试验现象分析模型,确定2D/3D?有无可简化部分?可压缩性?技术难点或者无法把握的问题?数值稳定性、存储空间、CPU占用时间的估计;2. 如可能,用简化模型先测试。

例如采用3D模型前先采用2D模型计算;3. 采用粗网格并使用能使模型快速求解的材料数据和载荷(例如,高粘度值,低速度,低压力,放松收敛精度等)。

此测试是确认模型具有合理性。

一旦模型出现问题,可以很快查出原因;4. 如果出现错误,可以查看*.out and *.log文件中的信息;5. 采用细网格通常更容易收敛;也有可能出现截断误差带来的影响,这样可以通过增大迭代次数、减小载荷增量、使用CFL数获得收敛。

5. 在进行流固耦合计算前,首先分别测试结构模型和流场计算模型;结构模型测试-在FSI边界施加相当于流体作用的压力;流体模型测试-将流固耦合边界定义为Wall或者移动的Wall;6. 当结构和流场模型能够正常求解,采用FSI进行求解;a)瞬态分析采用合理的初始条件;如果第一步不收敛,关注初始条件可能的影响;b)有预应力结构,控制预应力施加在FSI耦合之前;(第一步不施加流体载荷)c)先做大步长稳态计算,之后重启动瞬态计算;d)定义合理的时间步长。

高级结构非线性及流固耦合计算系统ADINA

高级结构非线性及流固耦合计算系统ADINA

有限元方法 ( 1 i e 1 m n M t o F n t E e e t e h d),最早是 由 中国工业近十年的高速发展和对有限元软件的功能需求引起 美国加州大学伯克利分校 ( 1e s t 0 a 1 o n , U v 1 f C ]f r 1 n m y a 用。在现在 ,有 限元 技术 不仅 在商业上有广阔的发展空 间, 在 理论和技术上也有更大 发展 。有限元方法概念的提 出,引 出 了美国加州大学伯 克利分校 有限元技术研究小组最为辉煌 的十年历程 。在这十年里 ,第一个 得到开发并公开发布 的自 动计算系统 ,称之为S P( tu t r n lss Po rm)。 A Sr cue A ay i rg a
AD A 为 当 今 最 为 可 靠 的 结 构 非 线 性 、 流 固 耦 合 计 算 系 种复合材料失效准则 ) I 成 N 。
D n m c n r m n a o ln a n lss 首 字 母 缩 写 。 y a is I ce e t l nie r A ay i的 N
AI U——前后处理模块包括 :基于P r s ld a ao i建模 内核、
这基本表达 了软 件开发者的根本 目标 ,PAD A 了求解线 P r s 1d . I 除 P N a a o 1几何接 口、有限元N s r n at a 文件 接 口、I E 通 用几 GS 性 问题外 ,还要具备分析 非线 性问题的强大功能——求解结 何传 输 、S 通 用几何 传输 、 点集数据 读入 、 自动 网格 划 T L 构 以及涉及结构 场外 的多场耦 合问题。结构非线性和 多场耦 分 、加载和边界条件 、模型列表 、结果列表 、等值线显示 、 合 的求解 ,是A IA D N 开发一贯坚持 的 目标 ,也是A I A D N 软件广 向量 显示 、流 场粒 子流显示 、动 画生成 、输 出格式 ( m 、 B P 大用户推崇 的特点之一 。 Je 等 ) pg 、用户 白定义图标 、在 线帮助文档 、宏语言 以及二 S r c ue t u t r—— 结构分 析模块 ,包括结构线性 ( 力、 静 1 8 年 以前 ,A I A 96 D 软件 的源代 码是 公开 的 ,这主要 次开 发资源库 。 N 包括两个著 名版本 ,RA I A 8 版 和A I A 4 。当时 , P D N 1 D N 8 版

油箱流固耦合模态分析

油箱流固耦合模态分析

做了油箱一个流固耦合模态分析的例子,共享给大家。

汽车油箱流固耦合模态分析实际的油箱几何结构很复杂,这是他提供的一个简化的模型。

几何参数:油箱容积42L,油液装载体积:21L;油箱材料参数:密度0.934g\cm3 ,弹性模量1100MP ,泊松比0.4 ,厚度5mm,边界条件为底部四边全约束。

油液参数:密度680kg/m3,体积模量1.3E9N/m2。

1. 启动ADINA,选择,2.点击,选择红色部分,设置箱体材料参数点击OK。

然后点击红色部分设置势流体油液,设置如下:点击OK。

关闭材料设置选项卡。

3.点击,如下设置4.点击,设置如下5.点击,进行如下设置面:6.点击设置拉伸体:7.显示如下8.通过面6继续拉伸体9.显示如下10.划分网格,进行如下操作点击OK。

11.点击,如下设置点击OK。

12.点击,如下设置连续两次点击OK。

13.设置然后进行如下设置:14.设置自由面15.加重力g。

点击红色define设置:最后设置16.保存ms.idb。

然后另存一个名为mm.idb。

17.静力计算,打开ms.idb,点击,求解ms。

18.17步求解结束后,关闭,然后打开mm.idb。

进行如下设置。

19.选择,点击,进行如下设置:分析前100阶模态,选用Determinant-search法求解流固耦合模态:20.重启动设置。

21.点击,输入mm,点击保存,提示选择重启动文件,选择ms.res,点击copy,程序求解。

22后处理第一阶第二阶第四阶第五阶第六阶MODENUMBER FREQUENCY1 8.16987E-012 1.49813E+003 1.49813E+004 1.78437E+005 2.12175E+006 2.12289E+007 2.24613E+008 2.24613E+009 2.52886E+0010 2.60776E+0011 2.60776E+0012 2.67783E+00第一步静力计算的放大云图。

ADINA 流固耦合FSI

ADINA 流固耦合FSI

ADINA 和流固耦合(FSI )ADINA 独具的流固耦合求解功能可以在单一系统ADINA 中模拟流体和因大变形、非弹性、接触及温度而经历明显的非线性响应的结构之间完全耦合的物理现象。

一个完全耦合的流固耦合模型意味着结构的变形影响流体区域,反过来流体的作用力也会施加到结构上。

从流体的角度看,Navier-Stokes 流体可以是不可压的,轻微可压的,低速和高速可压的。

从结构的角度看,各种结构单元类型都可以参与FSI 过程(即壳单元,2D 和3D 结构单元,梁单元,等参梁单元,接触面等),支持各种材料模型、支持各种非线性物理过程如材料失效、单元生死、结构失稳、相变等等。

此外,ADINA 还提供了针对流体是势流理论的完全耦合的流固耦合模型。

但由于势流体计算理论相对简单,不是本文主要讨论的内容。

ADINA FSI 是如何工作的?ADINA 在一个单一系统中组合了结构和流体动力学方程,获得这个系统的统一方程组,并对其进行求解。

对流体模型可以选择基于节点的FCBI (Flow-Condition-Based Interpolation)算法和基于单元的FCBI-C 算法进行单元的定义。

• FCBI 单元算法:基于速度自由度的FCBI 算法是用来提供稳定性的。

有限元方程可以通过Newton-Raphson 迭代计算一致的Jacobian 矩阵来进行求解。

因此流固耦合系统中能建立一致的刚度阵可以解决极为复杂的非线性问题。

• FCBI-C 单元算法:所有的解变量定义在单元的中心,速度和压力间的耦合是迭代地处理的。

因此在使用FCBI-C 单元算法的FSI 分析中,结构模型和流体模型之间的耦合也是迭代地处理的。

这种算法可以用来解决很大计算规模的实际问题。

这些算法可以适用于从低雷诺数到高雷诺数的各种流体。

一旦计算区域的任何一部分发生变形,对流体的Eulerian 描述就不再可用了。

因此,ADINA 求解流体的控制方程使用Arbitrary-Lagrangian-Eulerian(ALE)表示。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

实例3 隧道内具有柔性结构的流固耦合分析问题:隧道内具有柔性结构的流固耦合如图3-1所示。

图3-1 流体-固体结构示意图一、目的1. 掌握流固耦合作用FSI在Adina-AUI中的操作过程。

2. 掌握用伸缩比例因子画流固耦合模型。

3. 定义引导点(leader-follower points)。

二、定义模型主控数据1. 定义标题:选Control→Heading→敲入标题“exe03: Fluid flow over a flexible structure in a channel, ADINA input”→and click OK。

2. FSI分析:在右边Analysis Type区选FSI按钮。

3. 主控自由度选Control→Degrees of Freedom→不选X-Translation, X-Rotation, Y-Rotation andZ-Rotation按钮→and click OK。

4. 分析假设:大位移,小应变。

选Control→Analysis Assumptions→Kinematics→设置“Displacements/Rotations”为Large→ click OK。

(注:非常薄的结构,因此为小应变)。

三、力学模型1. 柔性结构建立模型1). 柔性结构几何模型坐标点如表3-1,几何结构如图3-2所示。

其几何面见表3-2所示。

①选Define Points 图标→按表3-1输入几何点坐标→ click OK .②选Define Surfaces 图标→设置TYPE 为Vertex → click OK(如图3-2所示)。

2). 施加固定边界条件和流-固边界条件①.图3-2中,在L2线上施加固定约束,其过程可用Adina-AUI 完成。

②. 流-固边界,选Model →Boundary Conditions →FSIBoundary →add FSI boundary number 1→在表中头两行敲入流固边界线编号1和 3 and click OK 。

3). 定义材料特性弹性模量1.0×106(dyne/cm 2),泊松比0.3。

(线弹性问题)选Model →Materials →Elastic_Isotropic →add material 1, 设置弹性模量1.0E6→泊松比 0.3 and click OK .4).定义单元和单元划分(1). 2-D 实体单元,此问题属平面应变问题。

Element group : 选 Meshing →Element Groups → 增加单元组号 1→ 设置 the Type to 2-D Solid →设置 the Element柔性结构图3-2 几何模型表3-1 模型几何点坐标 几何点 X1 X2 X3 坐标系 130.025 15.0 0 230.0 0.0 0 330.050.0图3-3 结构网格Sub-Type to Plane Stain and click OK .(2). 网格密度划分所有点设置均匀密度:①. Choose Meshing →Mesh Density →Line →选线编号1→设置划分为5段→在表中第一行敲入线编号3→ and click OK . ②. 选 Mesh Surface 图标→设置 “Number of Nodes per Element ” 为 9, 选择“Triangular Surfaces Treated as Degenerate ” 按钮→在表中第一行敲入面编号1→click OK . (如图3-3所示)。

③. 用伸缩比例因子显示网格模型,选Modify Meshplot 图标→click the View… 按钮→设置X 伸缩比例为100.0 →在对话框中click OK 两次(结果如图3-3所示)。

5). 保存结果选 Data File/Solution 图标, 敲入文件名exe03_a →不选Run ADINA 按钮→ click Save 。

2. 流体建立模型(Adina-F)1). 在菜单中选 New → 建立一个新的 ADINA-IN 数据库。

2). 从模块下拉式菜单中选ADINA-F 。

3). 设置FSI 分析:在右边的分析类型中选定FSI 选项。

4). 流体假设:选Model →Flow Assumptions →设置 Flow Dimensions 为2D (在YZ 平面) →设置the Default Upwinding Type 为 Finite Element →去掉Includes Heat Transfer 选项→and click OK 。

5). 时间步长和时间函数: 垂直入口分20步计算。

①. 选 Control →Time Step →在表中第一行敲入20→ click OK 。

②. 选 Control →Time Function →编辑表3-3→click OK 。

6). 流体几何模型建立坐标点如表3-4,几何结构如图3-4所示。

选Define Points 图标→按表3-4敲入坐标数据→click OK 。

(没有按比例化)图3-4 流体几何模型表3-3 时间函数 time value 0.0 0.0 20.00.04表3-4 模型几何点坐标 几何点 X1 X2 X3 1 500 40 2 30.025 40 3 0 40 4 500 15 5 30.025 15 6 0 15 7 500 0 8 30.05 0 9 30 0 10选Define Surfaces 图标→设置TYPE 为Vertex→ 按表3-5建立几何面→click OK (如图3-4所示)。

3). 定义材料特性选Model →Materials →Constant →add material 1→设置速度 Viscosity 为 1.7E-4→设置 Density 为0.001 →and click OK 。

4). 定义边界条件和载荷①. 壁面边界条件:在上下壁上施加无滑移(no-slip)壁面边界条件。

选Model →Special Boundary Conditions →add special boundary condition 1 →检查类型为Wall 。

在Line#表中头四行分别敲入壁面线编号1, 5, 9 and 12 →Click OK 。

②.流固耦合边界条件(FSI):选 Model →Special Boundary Conditions →add special boundary condition 2→ Type to Fluid Structure Interface ,且确保Fluid-Structure Boundary 为 1。

在Line#表中头两行分别敲入壁面线编号8 and 13 →Click OK 。

③. 载荷:在隧道入口施加垂直载荷选 Apply Load 图标→设置Load Type → Normal Traction , 并在右边Load Number field 按 Define... 按钮→ add Normal Traction 1, 设置 Magnitude 为1.0 →and click OK 。

在Apply Usual Boundary Conditions/Loads 对话框表的头两行Site#分别敲入线编号 6 和11→click OK 。

图3-5 模型载荷及其边界条件5). 定义流体导向关系为了保持好的网格质量,需要几何点2跟随几何点5沿隧道壁面移动。

选Meshing→ALE Mesh Constraints→Leader-Follower。

在表的第一行设置Label # to 1,设置Leader Point # to 5→设置Follower Point # to 2→ and click OK。

6). 定义网格单元①选Meshing→Element Groups→ add element group 1→确保单元类型为2-D Fluid→设置Element Sub-Type为Planar and click OK。

②选Meshing→Mesh Density→Surface→每个面划分的段数如表3-6所示→and click OK(如图3-6所示)。

注: 在图3-6中沿柔性结构有6个流体单元,在图3-3中只有5个固体单元,而且流体单元为3-节点线性单元,固体单元为9-节点四边形单元。

 ③网格划分:选Mesh Surface图标→在表中头四行分别敲入面积编号1, 2, 3, 4→ click OK。

7). 保存结果选Data File/Solution图标, 敲入文件名exe03_f→不选Run ADINA-F按钮→ click Save。

四、运行ADINA-FSI 和其后处理1.运行ADINA-FSI选 Solution →Run ADINA-FSI →选exeo3_f 文件,然后按住Ctrl 建选择exe03_a 文件→ click Start 。

2. 后处理ADINA-PLOT①. 从模块菜单中选ADINA-PLOT ;②. 从主菜单中打开数据结果文件,调入后处理文件exe03_f.por 和exe03_a.por 。

其处理结果如图3-7所示。

注: 图3-7中,在结构上面的垂直单元边界仍然保持垂直,因为在几何点5和2直接定义了导向关系。

 ③. 速度矢量:选 Model Outline →Load Plot 图标显示载荷→Quick Vector Plot 图标(AUI同时显示流体速度和结构应力)→Modify Vector Plot 图标→确保Element Vector Quantity 为Stress →选 Delete 按钮→click Yes → click Cancel orClose (如图3-8所示)。

④. 压力场:选Clear Vector Plot 图标→ Quick Bandplot 图标(AUI 同时显示流体压力和结构等效应力) → Modify Bandplot 图标→确保Band Plot Variable 为 (Stress:EFFECTIVE_STRESS )→Delete 按钮→click Yes →Cancel or Close(如图3-9所示)。

可用按Previous Timestep, Next Timestep, First Timestep and Last Timestep图标显示每步的结果。

图3-7 移动网格图3-8 流体速度矢量分部图3-9 压力场分布⑤.画结构移动: 只显示ADINA 固体模型→click Clear 图标→Display Zone 图标→选ADINA →and click OK 。

相关文档
最新文档