气、固、热耦合计算技术及其进展.ppt

气、热、固耦合计算技术及其进展 目录 1 耦合场分析方法概述 2 气、固、热耦合的物理过程及分析方法 2.1 物理过程分析 2.2 气动及传热分析 2.3 静强度分析 3 气、固、热耦合边界条件 4 气、固、热耦合求解过程分析 4.1 流体与固体传热的耦合分析 4.2 对流换热与固体热应变的耦合分析 4.3 流体动力与固体变形的耦合分析 4.4 流体动力、热应力和弹性变形的耦合分析 目录 5 气、固、热耦合信息的传递 5.1 温度参数空间差值传递方法 5.2 气压载荷的传递 5.3 几何变形的传递 5.4 三种几何变形传递方法 6 基于CBS有限元的流动-传热-变形耦合计算方法 6.1 CBS算法 6.2 基于CBS算法的流动-传热-变形耦合分析 7 小结 1 耦合场分析方法概述 所谓场是指物理量在空间区域上的分布。根据场的数学表达形式可分为数量场和矢量场。而耦合场分析就是指在分析的过程中，考虑了两个或多个工程物理场之间的交叉作用和相互作用。 根据耦合场的分析过程所耦合的物理场，耦合分析方法可分为两大类：顺序耦合和直接耦合。 1.1 顺序耦合 顺序耦合方法包括两个或多个按一定顺序排列的分析，每一种分析属于单一物理场分析。通过将前一分析的结果作为载荷施加到另一分析中的方式进行耦合。顺序耦合分析又分为顺序耦合物理场分析及顺序弱耦合分析。 1.1.1 顺序耦合物理场分析 顺序耦合物理场分析通过将第一个分析的结果做为第二个分析的载荷来耦合两个物理场。载荷传递发生在分析以外，必须通过物理环境指定载荷传递。这类分析的典型实例是顺序热-应力分析，热分析的节点做为体积载荷施加到随后的应力分析中。其实对于分开做的通风散热分析也属于这种形式，先进行通风计算，得到各处的流体速度作为边界条件输入到散热分析，计算出固体表面的散热系数从而进行散热计算。 1 耦合场分析方法概述 1.1.2 顺序弱耦合分析 顺序弱耦合分析一般用于进行流体-结构耦合分析。流体和结构分析的求解顺序发生，穿过流体和结构区域相同或不同的网格界面载荷在内部进行传递。这种分析需要从流体到结构传递流体力和热通量，并将结构的位移和温度传递给流体的结构，进行结构作用分析。 1.2 直接耦合方法 直接耦合方法，一般只进行一次分析，它使用包含多场自由度的耦合单元。通过计算所包含物理量的单元矩阵或载荷向量的方式进行耦合。 1.3 直接耦合法与顺序耦合法的应用场合 对于顺序耦合法，由于每种分析是相对独立的，对于多场的相互作用非线性程度不是很高的耦合场，顺序耦合方法更有效也更灵活。在分析过程中不同物理场间的迭代不断进行直至达到了所期望的收敛准则。顺序耦合可以是双向的，不同物理场之间进行相互耦合分析，直到收敛到一定的精度。 当耦合场之间的相互作用是高度非线性时，直接耦合较具优势。它使用耦合变量一次求解得到所有物理场的结果。 2 气、固、热耦合的物理过程及分析方法 2.1 物理过程分析 流、热、固耦合是指在流体流动、热量传递和固体弹性变形的系统中三者之间的相互影响和相互作用。热应力与流体压力导致固体变形，固体变形导致流场变化，进而使温度场又发生变化，流场、温度场和结构变形交作用直至流场、温度场和结构变形均达到一个稳定的平衡状态（假设无扰动）。 流、热、固三场间的耦合作用是同时发生的、双向相互的，如果单向顺序分析流固、流热、热固各个效应的耦合过程可能与实际情况存在差异。 速度场的发展、热量的传递以及固体的弹性变形三者的发生速度是不一致的。假设t为当前时刻，在极其细小的时间段里，速度场以其自身的发展速度发生着变化，而在这段时间里实际温度场以不同于速度场发展的速度进行着变化，同时，固体结构也在以自身的变形速度变化着，并且影响着流场。经过时间后，流场的约束条件变了，流场将朝着与上一个时间段不同的发展方向去变化，而温度场也已开始在新的时间段、新的条件下发生改变。但是，在数值计算的过程中，如果先分析经过时间后的t+△t时刻的速度场发展，然后，将t +△t时刻速度场视为稳定状态，并在此基础上进行温度场数值计算，那么，所计算得到的温度场不一定与实际经过时间后的t +△t时刻的温度场一致。 2 气、固、热耦合的物理过程及分析方法 另外，实际经过时间后的t+△t时刻的结构变形量也不一定就是在t + △t时刻温度场基础上求解得到的弹性变形量。 但是，数值计算本身就是一种近似求解。在数值求解的过