ADINA流固耦合软件使用手册.pdf-资料库

cwang1117-3180686-4744302543377986577.pdf-第1页.png

第1页 / 共300页

cwang1117-3180686-4744302543377986577.pdf-第2页.png

第2页 / 共300页

cwang1117-3180686-4744302543377986577.pdf-第3页.png

第3页 / 共300页

cwang1117-3180686-4744302543377986577.pdf-第4页.png

第4页 / 共300页

cwang1117-3180686-4744302543377986577.pdf-第5页.png

第5页 / 共300页

cwang1117-3180686-4744302543377986577.pdf-第6页.png

第6页 / 共300页

cwang1117-3180686-4744302543377986577.pdf-第7页.png

第7页 / 共300页

cwang1117-3180686-4744302543377986577.pdf-第8页.png

第8页 / 共300页

第一章计算流体力学概述第一章计算流体力学概述 1.1 计算流体力学概述计算流体力学（Computational Fluid Dynamics，简称 CFD）是以计算机作为模拟手段，运用一定的计算技术寻求流体力学各种复杂问题的离散化数值解的计算方法。计算流体力学可以看作是对基本守恒方程（质量守恒方程、动量守恒方程、能量守恒方程）控制下的流动过程进行数值模拟。通过这种数值模拟，可以得到极其复杂问题的流场内的各个位置上的基本物理量（如速度、压力、温度、浓度等）分布，以及这些物理量随时间的变化情况。计算流体力学、理论流体力学、实验流体力学是流体力学研究工作的三种主要手段。理论分析具有普遍性，各种影响因素清晰可见、为实验和计算研究提供依据。对于非线性情况，只有少数问题能给出解析解。实验研究仍是研究工作的基石，数值研究的许多方面都密切依赖于实验研究提供数据；计算结果需由实验验证；观察实验现象分析实验数据以建立计算模型等等。数值模拟是特殊意义下的实验，也称数值实验，它比起实验研究，经济效益极为显著。三种手段既互相独立又相辅相成。近年来，由于实际工程设计对于流体计算提出越来越高的要求，计算流体力学在明显地突破传统的单纯流体的观念，各种涉及到复杂物理现象的流体问题求解方法是计算流体力学发展的主要趋势，这些复杂现象是涉及热传递、多物质流动、相变、流固耦合体系求解、变边界（变流动区域）、湍流模拟等等。从工程角度看，流体力学研究的起因通常是基于对各种工程结构的设计需要，例如分析飞机机翼在气流作用下随机摆动问题的目的，是要求流场计算结果要对机翼的非稳态振动、强度特性提出明确的力学设计指标。因此可以说，对于相当多的流体计算问题，实际上我们需要知道的是一个耦合力学系统的响应特性－流固耦合体系特性，尽管这是一个更为复杂的计算体系。流固耦合（Fluid-Structure Interaction，简称 FSI）计算方法的开发和应用是目前工程计算流体力学发展的重点领域，也是计算流体力学指导工程设计的直接途径。计算流体力学的发展和计算机硬件求解能力、工程设计需求高速增长密切相关，可以肯定地说计算流体力学在未来的研究领域和工程领域，都会越来越走向实用化，越来越发挥不可或缺的作用。 1

第一章计算流体力学概述 1.2 基本概念流体是气体和液体的总称。在人们的生活和生产活动中随时随地都可遇到流体，大气和水是最常见的两种流体。 1.2.1 理想流体和粘性流体流体在静止时虽不能承受切向力，但在运动时，任意相邻两层流体之间却是有相互抵抗作用力的，这种相互抵抗的作用力称为剪切力，流体所具有的这种抵抗两层流体相对滑动速度的性质称为流体的粘性。粘性是流体的固有属性之一，不论流体处于静止还是流动，都具有粘性。它是流体状态（压力、温度、组成）的函数。气体的粘性随温度的升高而增大，液体的粘性随温度的升高而减小。自然界中存在的流体都具有粘性，具有粘性的流体统称为粘性流体。完全没有粘性的流体称为理想流体。自然界中并不存在真正的理想流体，它只是为便于处理某些流动问题所做的假设而已。 1.2.2 牛顿流体和非牛顿流体牛顿内摩擦定律： τ μ = lim n 0 Δ → u Δ n Δ = u ∂ μ n ∂ 其中 τ表示流体内摩擦应力， nΔ 为法线方向的距离增量； uΔ 为对应于的流体速度 nΔ 增量。牛顿内摩擦定律表示流体内摩擦应力和单位距离上的两层流体间的相对速度成正比。比例系数μ称为流体的动力粘度，简称粘度。牛顿流体：是指μ为常数的流体，即遵循牛顿内摩擦定律的流体称为牛顿流体。非牛顿流体：不符合上述条件的均称为非牛顿流体。所有气体和大多数低分子量液体均属牛顿型流体，如水、空气等；而某些高分子溶液、油漆、血液等则属于非牛顿流体。下面用一个简单的例子来说明牛顿流体和非牛顿流体的差别。ADINA 的操作步骤参见第二部分。 2

例 1 分别考虑两个平行板间的牛顿流体和非牛顿流体，如下图所示，给出流体速度大小为第一章计算流体力学概述 10m/s。 v 先考虑牛顿流体流体参数为常参数模型，密度为 1 /kg m 3 ，粘度为 0.04 /N s m⋅ 2 。查看流场中纵坐标为 0.3 处的节点的速度大小，如下图所示： m 非牛顿流体非牛顿流体模型和牛顿流体模型的差别在于材料的定义。这里非牛顿流体使用 Power Law 模型。流体参数：密度为 1 /kg m 3 ，常数 μ = 0 0.04, A = 1, n = − 0.75 。查看非牛顿流体模型的计算结果，仍然画相同节点的处的速度图像，如下图所示，可以看出，非牛顿流体的速度增长较牛顿流体稍缓一些。 3

第一章计算流体力学概述 1.2.3 可压缩流体和不可压缩流体根据密度 ρ是否为常数，流体分为可压缩流体（compressible ）与不可压缩流体 (incompressible)两大类。在温度不变的情况下，当密度 ρ为常数时，流体为不可压缩流体，否则为可压缩流体。空气为可压缩流体，水为不可压缩流体。水的可压缩性是很小的，压强每增加一个大气压，其体积变化不到万分之一。工程中常用的其它工作液体，如液压油，机械油等，其体积模量数值也都很大，在一般工程计算中，可以忽略其可压缩性，将其看作是不可压缩流体。气体的可压缩性与液体相比则大很多，因此在研究气流场的时候，当流速较低 (<0.3 马赫)时可以认为是不可压缩流体；当流速较高时需要考虑压强对气体密度的影响，甚至温度、压强对体积和密度的影响都必须考虑。 1.2.4 定常与非定常流动根据流体的物理量（如速度、压力、温度等）是否随时间变化，将流动分为定常（steady ）与非定常（unsteady）两大类。当流动的物理量不随时间变化，即又称稳态流动；当流动的物理量随时间变化，即 () ∂ t ∂ ≠ 0 动或瞬态流动。 1.2.5 层流和湍流 () ∂ t ∂ = 0 时，为定常流动，时，为非定常流动，又称非稳态流自然界中的流体流动状态主要分为两种形式－层流和湍流。在平行管的流动试验中，从流动质点的运动状态看，当仅发生层流时，流体质点互不干扰，流体流动呈线性或层状，且平行于管道轴线；当湍流（又称紊流）发生后，质点运动杂乱无章，除了平行于管道轴线的 4

第一章计算流体力学概述运动以外，还存在着剧烈的横向运动。层流和湍流是两种不同性质的流态。层流时，流体流速较低，质点受粘性制约，不能随意运动，粘性力起主导作用；湍流时，液体流速较高，粘性的制约作用减弱，惯性力起主导作用。液体流动时，究竟是层流还是湍流，通常要用雷诺数来判定。雷诺数（Reynolds）是表征流体流动特性的一个重要参数。它是由管内的平均流速、 u 液体的运动粘度 ν、管径三个参数决定的，即 d Re = 。又有运动粘度ν与动力粘度μ的 ud ν 关系， ν = ，ρ为密度，所以雷诺数又可以表示为 Re μ ρ = udρ μ 。从层流向湍流过渡发生在雷诺数约等于 2300，即临界雷诺数的时候。当流体的实际流动时的雷诺数小于临界雷诺数时为层流，反之则为湍流。对于非圆截面管道来说，做简单等效，可用 R 代替上式中的， R 为通流截面的水力 d 半径。它等于液流的有效截面积 A 和它的湿周χ（通流截面上与液体接触的固体壁面的周长）之比，即 AR = 。 x 关于层流和湍流的区别，可用典型的后台阶流动模型进行说明。这部分内容请参考第四章。 1.3 流体力学方程组流体流动满足三个基本方程：（1）质量守恒方程。它反映的是物质不生不灭这一最自然的物理定律。它用以描述流体密度 ρ的变化规律。它不需要补充任何其他关系式，也就是说，质量守恒方程是物质不生不灭的最直观的体现，方程的形式十分简单。质量守恒方程也称连续性方程。（2）动量守恒方程。它反映的是牛顿定律，即物体在力的作用下做加速运动。具体说，物质体所受的合力等于其质量与加速度的积，也可以理解为流体微团所受的力等于其动量变化率。因此，只要能求出合力，便可以得到动量守恒方程。合力包括体积力 fρ 和面积力 nP 。体积力可以是多相物质的相互作用力、远程力（如重力与电磁力）和惯性力。面积力主要是压力和粘性应力。 5

第一章计算流体力学概述（3）能量守恒方程。能量守恒属于经典的热力学定律。流体微团单位质量的能量即总能 E 包括内能 e 与动能 1 2 V V⋅ 。合力所做的功、热传导、（由化学反应等引起的）生成热都引起总能的变化。 1.3.1 预备知识流体速度 =u u u u ( , 1 3 , 2 ) ， u u u , 1 , 2 分别表示 x,y,z 方向上的速度分量。 3 ρ是密度，p 是压强，T 是热力学温度。若有过点 =x ( x x x , 1, 2 3 ) 的面积微元，单位法向量为 n 。 dS 在[ , t t dt+ ] 内沿 n 方向流过 dS 的流体体积为 ⋅u n dSdt 在[ , t t dt+ ] 内沿 n 方向流过 dS 的流体质量为 ρ ⋅u n dSdt 在[ , t t dt+ ] 内沿 n 方向流过 dS 的流体动量为其中 u ⊗ u = 2 u 1 u u 2 1 u u 3 1 ⎛ ⎜ ⎜ ⎜ ⎝ u u 1 2 u 2 2 u u 3 2 u u 1 3 u u 2 3 u 2 3 ⎞ ⎟ ⎟ ⎟ ⎠ ρ ( u u n) ⋅ dSdt = ( ρ u ⊗ ) u n dSdt 面积微元受到 n 正向一侧的流体压力为 dS − n （压力方向与法线方向相反） p dS 热力学中，有状态方程 p Tρ= , ( f ) （1.1）若 p RTρ= ，则为理想气体。用 e 表示单位质量流体的内能，单位体积中流体的能量为 2 （ u = 2 u 1 2 + u 2 2 + u 3 ）在[ , t t dt+ ] 内沿 n 方向流过的流体能量为 dS ( 1.3.2 理想流体力学方程组（1）质量守恒方程 u ρ ρ+ e 1 2 2 ) ⋅u n dSdt eρ ρ+ 1 2 2 u ，质量守恒定律表述为：单位时间内流体微元中质量的增量，等于同一时间间隔内流入该微元的净质量，即 6

资料库

ADINA流固耦合软件使用手册.pdf

相关推荐

课程资源

热门标签

最新资料