ANSYS流体分析CFD

第一章 FLOTRAN 计算流体动力学(CFD)分析概述

FLOTRAN CFD 分析的概念

ANSYS程序中的FLOTRAN CFD分析功能是一个用于分析二维及三维流体流动场的先进的工具，使用ANSYS中用于FLOTRAN CFD分析的FLUID 141和FLUID 142 单元，可解决如下问题：

?作用于气动翼(叶)型上的升力和阻力

?超音速喷管中的流场

?弯管中流体的复杂的三维流动

同时，FLOTRAN还具有如下功能：

?计算发动机排气系统中气体的压力及温度分布

?研究管路系统中热的层化及分离

?使用混合流研究来估计热冲击的可能性

?用自然对流分析来估计电子封装芯片的热性能

?对含有多种流体的(由固体隔开)热交换器进行研究

FLOTRAN 分析的种类

FLOTRAN可执行如下分析：

?层流或紊流

?传热或绝热

?可压缩或不可压缩

?牛顿流或非牛顿流

?多组份传输

这些分析类型并不相互排斥，例如，一个层流分析可以是传热的或者是绝热的，一个紊流分析可以是可压缩的或者是不可压缩的。

层流分析

层流中的速度场都是平滑而有序的，高粘性流体（如石油等）的低速流动就通常是层流。

紊流分析

紊流分析用于处理那些由于流速足够高和粘性足够低从而引起紊流波动的流体流动情况，ANSYS中的二方程紊流模型可计及在平均流动下的紊流速度波动的影响。如果流体的密度在流动过程中保持不变或者当流体压缩时只消耗很少的能量，该流体就可认为是不可压缩的，不可压缩流的温度方程将忽略流体动能的变化和粘性耗散。

热分析

流体分析中通常还会求解流场中的温度分布情况。如果流体性质不随温度而变，就可不解温度方程。在共轭传热问题中，要在同时包含流体区域和非流体区域（即固

体区域）的整个区域上求解温度方程。在自然对流传热问题中，流体由于温度分布的不均匀性而导致流体密度分布的不均匀性，从而引起流体的流动，与强迫对流问题不同的是，自然对流通常都没有外部的流动源。

可压缩流分析

对于高速气流，由很强的压力梯度引起的流体密度的变化将显著地影响流场的性质，ANSYS对于这种流动情况会使用不同的解算方法。

非牛顿流分析

应力与应变率之间成线性关系的这种理论并不能足以解释很多流体的流动，对于这种非牛顿流体，ANSYS程序提供了三中粘性模式和一个用户自定义子程序。

多组份传输分析

这种分析通常是用于研究有毒流体物质的稀释或大气中污染气体的传播情况，同时，它也可用于研究有多种流体同时存在（但被固体相互隔开）的热交换分析。

第二章 FLOTRAN分析基础

FLOTRAN单元的特点

ANSYS中的FLOTRAN单元，即FLUID141和FLUID142，用于解算单相粘性流体的二维和三维流动、压力和温度分布。对于这些单元，ANSYS通过质量、动量和能量三个守恒性质来计算流体的速度分量、压力、以及温度。

FLUID141单元

FLUID141单元具有下列特征：

维数：二维

形状：四节点四边形或三节点三角形

自由度：速度、压力、温度、紊流动能、紊流能量耗散、多达六种流体的各自质量所占的份额

FLUID142单元

FLUID142单元具有下列特征：

维数：三维

形状：四节点四面体或八节点六面体

自由度：速度、压力、温度、紊流动能、紊流能量耗散、多达六种流体的各自质量所占的份额

FLUID141单元 FLUID142单元

FLUID单元的其他特征

FLUID单元的其他特征包括：

?用于模拟紊流的二方程紊流模式

?有很多推导结果，诸如：流场分析中的马赫数、压力系数、总压、剪应力、壁面处的y-plus、以及流线函数；热分析中的热流、热交换（膜）系数等。

?流体边界条件，包括：速度、压力、紊流动能以及紊流能量耗散率。用户无需提供流场进口处紊流项的边界条件，因FLOTRAN对此提供的缺省值适用

于绝大多数分析。

?热边界条件，包括：温度、热流、体积热源、热交换（膜）系数。

用户可使用的坐标系有：的卡尔坐标系、柱坐标系、极坐标系和轴对称坐标系。如果所计算的问题是轴对称的，激活旋转（swirl）选项即可算出垂直于对称平面的速度分量。

使用FLOTRAN单元的一些限制及注意事项

FLOTRAN单元的一些局限性：

?在同一次分析中不能改变求解的区域

?单元不支持自由流面边界条件

?ANSYS程序的某些特征不能同FLOTRAN单元一起使用

?使用FLOTRAN单元时不能使用某些命令或菜单

?当使用ANSYS的图形用户界面时，程序将只能显示那些在菜单和对话框中的 FLOTRAN SetUp部分要求了的特征和选项。

FLOTRAN单元使用中的一些限制

当使用FLOTRAN单元时，要避免使用ANSYS的某些特征和命令，至少，要注意到在使用FLOTRAN单元时与别的分析稍微有些不同，当使用了无效的命令时，程序会给出相应的警告或错误信息。使用FLOTRAN单元要注意如下几点（下面所列命令相应的菜单路径请参见ANSYS命令手册或联机帮助中的“Commands and Their Location in the GUI”

?FLOTRAN单元不能和其他单元联合使用。

?节点坐标系必须与总体坐标系一致。

?/CLEAR命令并不破坏业已存在的FLOTRAN结果文件（Jobname.RFL），这有助于防止用户不小心破坏那些花了很多时间和精力才求得的结果，用户必

须在操作系统里才能删除那些无用的结果文件。

?CP命令通过对自由度进行耦合来形成周期边界条件，ANSYS命令手册对CP 命令的描述是可以只对某些自由度进行耦合，但作FLOTRAN分析时，周期

边界的所有自由度都将被耦合。

用户不能对同一个单元中的节点进行耦合，相邻单元间节点耦合也很困难。

?ADAPT命令不适用于FLOTRAN分析。

?不能用ANTYPE命令来引入FLOTRAN的瞬态分析。

?FLOTRAN分析不支持自动时间步长功能，详见“FLOTRAN瞬态分析”。

?如果用户通过BFCUM、BFDELE或BFUNIF来定义节点热源，则ANSYS会在内部用BFE命令来代替。

?不能使用LDREAD,FORC命令来电磁载荷转换到FLOTRAN分析中，而必须使用相应的宏来进行转换。

?FLOTRAN分析不能使用CE、CECMOD、CEDELE和DEINTF命令。

?FLOTRAN分析不能使用CNVTOL命令来设置收敛容差。

?不能用DSYM命令来定义FLOTRAN的对称和反对称边界条件。

?FLOTRAN不支持旋转坐标系中的角加速度向量。

?在FLOTRAN分析中，使用FLDATA4,TIME命令而不是DELTIM命令来定义一个载荷步的时间。

?对于FLOTRAN分析，不能使用DESOL命令或PRESOL命令来修改节点的热（HEAT）、流（FLOW）、或流密（FLUX）结果。

?FLOTRAN分析不允许将积分点结果外推到节点上（ERESX命令）。

?FLOTRAN分析不允许通过镜象操作来形成单元。

?FLOTRAN分析不能用KBC命令来施加渐变载荷，而必须用多个载荷步来逐渐改变载荷。

?FLOTRAN分析不允许用LCCALC、LCDEF、LCFA、LCFI等命令来作不同载荷状况之间的运算操作。

?NCNV命令中的收敛工具不能用于FLOTRAN中相互独立的求解器。

?FLOTRAN分析不允许用NEQUIT命令来定义非线性分析的平衡迭代数。

?FLOTRAN分析要求节点坐标系必须是总体的卡尔坐标系，故不能用N命令、NMODIF命令、和NROTAT命令的旋转（rotational）域。

?FLOTRAN分析中，FLDATA2,ITER命令和FLDATA4,TIME命令用来控制一个载荷步中的总体迭代数。

?FLOTRAN分析不允许用户自定义单元。

?FLOTRAN分析中，FLDATA2,ITER命令和FLDATA4,TIME命令用来对写入数据库中的结果进行控制。

?PRNLD命令不能用于FLOTRAN分析中，因其不能将边界条件作为可打印的单元节点载荷来保存。

?FLOTRAN分析中没有节点反力解。

?部分和预定义求解选项（PSOLVE命令定义）不适用于FLOTRAN各自独立的求解器。

?TIME命令不能用在FLOTRAN分析中。

?FLOTRAN用FLDATA1,SOLU命令而不是TIMINT命令来定义瞬态载荷步。

?FLOTRAN用FLDATA4,TIME命令而非TRNOPT命令来定义瞬态分析选项。

FLOTRAN分析的主要步骤

一个典型的FLOTRAN分析有如下七个主要步骤：

1. 确定问题的区域。

2. 确定流体的状态。

3. 生成有限元网格。

4. 施加边界条件。

5. 设置FLOTRAN分析参数。

6. 求解。

7. 检查结果。

第一步：确定问题的区域

用户必须确定所分析问题的明确的范围，将问题的边界设置在条件已知的地方，如果并不知道精确的边界条件而必须作假定时，就不要将分析的边界设在靠近感兴趣区域的地方，也不要将边界设在求解变量变化梯度大的地方。有时，也许用户并不知道自己的问题中哪个地方梯度变化最大，这就要先作一个试探性的分析，然后再根据结果来修改分析区域。这些在后面章节中都有详述。

第二步：确定流体的状态

用户在此需要估计流体的特征，流体的特征是流体性质、几何边界以及流场的速度幅值的函数。FLOTRAN能求解的流体包括气流和液流，其性质可随温度而发生显著变化，FLOTRAN中的气流只能是理想气体。用户须自己确定温度对流体的密度、粘性、和热传导系数的影响是否是很重要，在大多数情况下，近似认为流体性质是常数，即不随温度而变化，都可以得到足够精确的解。

通常用雷诺数来判别流体是层流或紊流，雷诺数反映了惯性力和粘性力的相对强度，详见第四章。

通常用马赫数来判别流体是否可压缩，详见第七章。流场中任意一点的马赫数是该点流体速度与该点音速之比值，当马赫数大于0.3时，就应考虑用可压缩算法来进行求解；当马赫数大于0.7时，可压缩算法与不可压缩算法之间就会有极其明显的差异。

第三步：生成有限元网格

用户必须事先确定流场中哪个地方流体的梯度变化较大，在这些地方，网格必须作适当的调整。例如：如果用了紊流模型，靠近壁面的区域的网格密度必须比层流模型密得多，如果太粗，该网格就不能在求解中捕捉到由于巨大的变化梯度对流动造成的显著影响，相反，那些长边与低梯度方向一致的单元可以有很大的长宽比。

为了得到精确的结果，应使用映射网格划分，因其能在边界上更好地保持恒定的网格特性，映射网格划分可由命令MSHKEY,1或其相应的菜单Main Menu>Preproce ssor > -Meshing-Mesh>-entity-Mapped来实现。

第四步：施加边界条件

可在划分网格之前或之后对模型施加边界条件，此时要将模型所有的边界条件都考虑进去，如果与某个相关变量的条件没有加上去，则该变量沿边界的法向值的梯度将被假定为零。求解中，可在重启动之间改变边界条件的值，如果需改变边界条件的值或不小心忽略了某边界条件，可无须作重启动，除非该改变引起了分析的不稳定。

第五步：设置FLOTRAN分析参数

为了使用诸如紊流模型或求解温度方程等选项，用户必须激活它们。诸如流体性质等特定项目的设置，是与所求解的流体问题的类型相关的，该手册的其他部分详细描述了各种流体类型的所建议的参数设置。

第六步：求解

通过在观察求解过程中相关变量的改变率，可以监视求解的收敛性及稳定性。这些变量包括速度、压力、温度、动能 (ENKE自由度) 和动能耗散率 (ENDS自由度) 等紊流量以及有效粘性（EVIS）。一个分析通常需要多次重启动。

第七步：检查结果

可对输出结果进行后处理，也可在打印输出文件里对结果进行检查，此时用户应使用自己的工程经验来估计所用的求解手段、所定义的流体性质、以及所加的边界条件的可信程度。

FLOTRAN分析中产生的一些文件

在ANSYS中进行的大多数流体分析都是通过多次中断和重启动来完成的，通常, 分析人员需要在各个重启动之间改变诸如松弛系数等参数或开关某些项（如求解温度方程的开关）。每当用户继续一个分析时，ANSYS程序会自动将数据附加在所有的由FLOTRAN单元产生的文件中。下面将对FLOTRAN单元产生的所有文件进行说明: ?结果文件， Jobname. RFL，包含节点结果。

?打印文件， Jobname.PFL，包含各量的收敛记录及进/出口状态(如流量等)。

?壁面文件， Jobname.RSW，包含壁面剪切应力以及Y-Plus信息。

?残差文件， Jobname.RDF，包含节点残差。

?调试文件， Jobname.DBG，包含数学求解器的有关信息。

?结果备份文件， Jobname.RFO，包含结果文件数据的一个拷贝。

?重启动文件， Jobname.CFD，包含FLOTRAN的数据结构。

结果文件

FLOTRAN分析的结果并不自动保存在ANSYS的数据库中，在每次求解之后，程序会将一个结果集附加在结果文件Jobname.RLF中。用户可对结果文件的内容及程序对结果文件的更新频率进行控制，ANSYS命令手册中对FLDATA5,OUTP命令的介绍就详细说明了结果文件会基于用户的选择而保存些什么内容。

在一个稳态FLOTRAN分析中，结果文件能保存多少个结果集是没有限制的，在求解的初期多保存几个结果有很多好处：可以比较各结果集之间的变化、可以使用不同的选项或松弛系数来从一个分析的较早状态重新开始分析。

当开始一个新分析时（在其第一次迭代之前），ANSYS程序会保存一个结果, 然后在当中断发生时保存再保存结果，在这些事件之间，用户还可通过设置将一些中间结果附在结果文件里，这样就可以从较早的分析状态开时，通过激活一些不同的选项和特征来重新分析，例如，可以通过这种方式来提高分析的稳定性。

使用ANSYS 的覆盖频率选项是一个明智的方法，它就可以周期性地保存和更新一个临时的结果集，这样，当由于断电或其他系统原因而发生求解中断时，总可以有一个可用的结果集用于重新开始分析。设置覆盖频率的方式如下：

命令：FLDATA2,ITER,OVER,value

菜单：Main Menu>Solution>FLOTRAN Set Up>Execution Ctrl

设置附加频率的方式如下：

命令：FLDATA2,ITER,APPE,value

菜单：Main Menu>Solution>FLOTRAN Set Up>Execution Ctrl

打印文件（ Jobname.PFL）

Jobname.PFL文件包含了所有FLOTRAN输入参数的完整记录，该信息每在发出一个求解命令时就保存一次以完整地记录整个分析历程。同时，所有激活了的变量的收敛过程也记录了下来，还有一个对结果的总结，即每个性质和自由度的最大最小值，这些记录的频率都由用户自己设定。所记录的其他量还有：各记录量的平均值、质量流的边界、质量平衡的计算、所有热传导和热源的相关信息。

节点残差文件

节点残差文件，即Jobname.RDF，显示了当前解的收敛好坏程度。在求解过程的每一个阶段，流场、性质场、温度场都用于对每个自由度计算系数矩阵和强迫函数，如果解完全收敛，这些矩阵和强迫函数将会生成一个与产生它们的速度场一样的速度场，同时，矩阵方程的残差也会变得很小。要得到一个残差文件，必须至少执行一次迭代。

当求解过程发生振荡时，残差的幅值将显示分析的错误所在。（矩阵的主对角元素对残差作归一化处理）这种归一化使用户可对自由度的值及其残差作比较。

对每一个激活了的自由度计算残差并将其存入残差文件的方式如下：

命令：FLDATA5,OUTP,RESI,TRUE

菜单：Main Menu>Solution>FLOTRAN Set Up>Additional Out>Residual File

要读取残差文件，可通过菜单Main Menu>General Postproc>FLOTRAN 2.1A或命令FLREAD来实现。

重启动文件

通常，FLOTRAN在一个重启动的起始处计算数据结构，对于一个大模型，这种计算将消耗大量的时间，为了避免这种重新计算，可要求FLOTRAN将数据结构保存在重启动文件 Jobname.CFD中，FLOTRAN从ANSYS的数据库中产生该文件。

对 Jobname.CFD文件的读和写的方式如下：

命令： FLDATA32,REST,RFIL,T

菜单： Main Menu>Preprocessor>FLOTRAN Set Up>Restart Options>CFD Restart File

可将RFIL状态设置为开(ON)或关(OFF)，若设为开，则FLOTRAN开始执行分析时将读入重启动文件，若此时重启动文件不存在，则将产生一个重启动文件。

如果在改变了边界条件之后再进行重启动分析，则必须覆盖掉业已存在的.CFD 文件以使得ANSYS能用新的边界条件进行重新分析，覆盖.CFD文件的方式如下：命令： FLDATA32,REST, WFIL,T

菜单： Main Menu>Preprocessor>FLOTRAN Set Up>Restart Options>CFD Restart File

这就使FLOTRAN在下一载荷步产生一个新的重启动文件，并自动将RFIL状态设置为关闭。当新的重启动文件产生之后，用FLDATA32,REST,RFIL,T命令使随后的重启动能使用新的重启动文件。

FLOTRAN重启动分析（续算）

用户可在结果文件Jobname.RFL中任意一个解集的基础上开始一个重启动分析，重启动位置的设置可基于解集号(NSET)、迭代数(ITER)、载荷步/子步号(LSTP)或瞬态分析的时间(TIME)，方式如下：

命令： FLDATA32,REST,lable,value

其中，lable为上面的NSET、ITER、LSTP、TIME等

菜单：Main Menu>Preprocessor>FLOTRAN Set Up>Restart Options> Restart/Iteratio（或Restart/Load step, Restart/Set, 等）

当重启动一个分析时，ANSYS将原始的结果文件拷贝到Jobname.RSO中并将重启动点、所有在重启动点之前的结果集、所有的后续结果集放在新的结果文件中。如果在FLDATA32,REST命令中的value值是一个负值，则将不产生Jobname.RSO文件，而重启动的点将由value的绝对值来指定。

提高收敛性和稳定性的常用的工具

ANSYS程序提供几个有助于收敛和求解稳定的工具，理论手册对其机理有详述。

松弛系数

松弛系数是一个其值介于0和1之间的小数，它表示旧结果与附加在旧结果上以形成新结果的最近一次计算量之间的变化量。设置松弛系数的方式如下：命令：FLDATA25,RELX,lable,value

菜单： Main Menu>Preprocessor>FLOTRAN SetUp>Relax/Stab/Cap>DOF Relaxation

Main Menu>Preprocessor>FLOTRAN SetUp>Relax/Stab/Cap>Prop Relaxation Main Menu>Solution>FLOTRAN SetUp>Relax/Stab/Cap>DOF Relaxation

Main Menu>Solution>FLOTRAN SetUp>Relax/Stab/Cap>Prop Relaxation

注：命令手册中对该命令的自由度和性质量有详述。

惯性松弛

对某个自由度的方程组的惯性松弛就是使其矩阵的主对角占优以保持求解的稳定性。如果当一个解在收敛过程中没有发生舍入误差，则惯性释放的值不会影响到求解的最终结果。但是通常的求解过程都会发生舍入误差，故惯性松弛可能对结果产生影响。用户可对动量方程（MOME）、紊流方程（TURB）。压力方程（PRES) 和温度方程（TEMP）施加惯性松弛，其方式如下：

命令：FLDATA26,STAB,lable,value

菜单：Main Menu>Preprocessor>FLOTRAN SetUp>Relax/Stab/Cap>Stability Parms Main Menu>Solution>FLOTRAN SetUp>Relax/Stab/Cap>Stability Parms

惯性松弛系数是以所加项的分母的形式出现的，故其值越小，所起作用越大，其典型值介于1.0（作用中等）到1.0×10－7（作用很大）之间。

人工粘性

人工粘性用于在梯度较大的区域平抑速度解。它有助于可压缩问题的收敛，也有助于对有分布阻力的不可压缩问题的速度解进行平抑。对于不可压缩问题，应使人工粘性的幅值与有效粘性的幅值处于相同的数量级。施加人工粘性的方式如下：命令：FLDATA26,STAB,VISC,value

菜单：Main Menu>Preprocessor>FLOTRAN SetUp>Relax/Stab/Cap>Stability Parms Main Menu>Solution>FLOTRAN SetUp>Relax/Stab/Cap>Stability Parms

速度限值

速度限值使所求解量不能超出用户所定义的值，可对速度、压力和温度自由度进行限制（VX、VY、VZ、PRES、TEMP），方式如下：

命令：FLDATA31,CAPP

菜单：Main Menu>Preprocessor>FLOTRAN SetUp> Relax/Stab/Cap>Results Capping Main Menu>Solution>FLOTRAN SetUp>Relax/Stab/Cap>Results Capping

速度限值可消除速度尖峰的不利影响，这种速度尖峰通常发生在收敛过程中的较早阶段。它还特别适合用于可压缩流分析，因这类分析中速度尖峰通常使动能项大到产生负的静温。

当对压力进行限值时，所限的值是由压力方程解算出来的压力而不是松弛后的压力，故当限值后作重启动时，压力值仍有可能超出限值。

注意：当有速度限值时，质量有可能不守恒。

面积积分阶次（Quadrature Order）

缺省的用于计算单元面积积分的阶次是单点积分，用户可对其进行控制。对于轴对称问题，求解时，该值自动设为2，因为当面积积分阶次为2时，可使含有异常形状单元的问题收敛到更精确的解。用下面的方式改变动量、压力、热或紊流项的面积积分阶次：

命令：FLDATA30,QUAD,lable,value

其中，lable为要改变的单元积分，value为积分点的数目。

菜单：Main Menu>Preprocessor>FLOTRAN SetUp> Mod Res/Quad Ord> CFD Quad Orders

Main Menu>Solution>FLOTRAN SetUp>Mod Res/Quad Ord>CFD Quad Orders

FLOTRAN分析过程中应处理的问题

确定总体迭代的数目

FLOTRAN分析是一个非线性的序列求解过程，故每次分析首先得确定要让程序执行多少次迭代。一次总体迭代就是对所有相关的控制方程按序列进行求解，并且在求解过程中流体性质会随时更新。在瞬态分析中，时间步循包含了总体迭代循环。在一个总体迭代中，程序首先获得动量方程的近似解，再在质量守恒的基础上将动量方程的解作为强迫函数来求解压力方程，然后用压力解来更新速度，以使速度场保持质量守恒。如果要求了程序求解温度，则程序会同时求解温度方程并更新与温度相关的流体性质。最后，如果激活了紊流模型，则程序将求解紊流方程并用紊流动能及其耗散率来计算有效粘性和热传导系数，有效粘性和热传导系数将分别代替层流粘性和热传导系数以在平均流上模拟紊流的影响。用下面的方式定义总体迭代的数目：命令：FLDATA2,ITER,EXEC,value（value即为迭代数）

菜单：Main Menu>Preprocessor>FLOTRAN SetUp> Execution Ctrl

Main Menu>Solution>FLOTRAN SetUp>Execution Ctrl

收敛监测

在FLOTRAN求解过程中，程序在每一个总体迭代里对每一个自由度计算出一个收敛监测量，这些自由度包括：速度(VX、VY、VZ)、压力(PRES)、温度（TEMP)、紊流动能（ENKE）、动能耗散率（ENDS）、以及激活了的多组份传输方程（SP01 ～SP06）。收敛监测量就是两次迭代之间结果改变量的归一化值，若以 表示任一自由度，则该自由度的收敛监测量可由下式表示：

收敛监测量表示变量在当前迭代（k th）的结果和前一次迭代（(k-1)th）的结果之间差值的总和除以当前值的总和，这种求和是在所有节点上进行的，并且使用的是差值的绝对值。

在批处理或交互式运行过程中，当求解进行时, 程序的“图形求解跟踪（GST）” 功能将实时显示出所计算的收敛监测量，GST的缺省值在交互运行时是开（ON）, 而在批处理运行时是关（OFF）。用户可用下面的方式定义其开关：

命令：/GST

菜单：Main Menu>Solution>Output Ctrls>Grph Solu Track

图2－1是两个典型的GST图形。图2－1(b)是一个FLOTRAN的瞬态分析过程，图中的每一个尖峰表示了一个新时间步的开始。

在初始阶段可能出现的一些振荡之后，收敛监测量的大小将随着分析过程的收敛而逐渐减小，但其减小的程度将依赖于几个因素，诸如：

?几何边界的复杂程度

?高梯度区域有限元网格的精度

?紊流的严重程度（由雷诺数确定）

?出口边界处流场的发展是否充分

当使用图形求解跟踪（GST）功能时，还应注意：

?不单是FLOTRAN分析有GST功能，非线性的结构分析、非线性的热分析和非线性的电磁场分析都有GST功能。详见各自的分析指导手册。

?GST可同时显示多达10条的跟踪曲线，如果用户的模型有多于10个的自由度, 则GST将只显示前10个自由度的收敛跟踪曲线。

?当GST开始显示时，程序会弹出一个带STOP按钮的对话框，用户可在任意时刻通过点取该STOP按钮来中断求解过程，而后要进行重启动分析时，可通过

执行命令SOLVE或其相应的菜单Main Menu>Solution>Run Flotran来实现。

图2－1 由GST显示的收敛监测量

（a）稳态求解

（b）瞬态求解

中断一个FLOTRAN分析

用户可以定义一个基于压力和温度收敛监测量的目标值来中断一个FLOTRAN分析，定义方式如下：

命令：FLDATA3,TERM,PRES,value

FLDATA3,TERM,TEMP,value

菜单：Main Menu>Preprocessor>FLOTRAN SetUp> Execution Ctrl

Main Menu>Solution>FLOTRAN SetUp>Execution Ctrl

压力和温度的收敛缺省值都是1.0×10-8，如果没有激活温度方程的求解，则程序只检测压力的收敛值是否满足要求，而若同时激活了流体方程和温度方程的求解，则二者的收敛标准都必须同时满足。在满足了压力和温度的收敛条件或总体平衡迭代数达到了所要求的值后，FLOTRAN求解过程就自动中断。

要中断一个正在以批处理方式或后台方式执行的FLOTRAN分析，则需在当前工作目录下生成一个Jobname.ABT文件，该文件的第一行应含有terminate字样，且该字样的起始位置应是第一行的第一列。在每一次总体迭代之前，FLOTRAN都会在当前目录下搜寻Jobname.ABT文件，如果程序找到该文件并发现其含有terminate字样,则立即完成该次总体迭代并正常中断程序的执行，而且将结果写入结果文件中。

对一个FLOTRAN分析进行评价

分析员必须回答的两个问题是：

1. 所作的分析是何时结束的？

2. 所作的分析是否是正确的？

这两个问题是相互关联的，因为，如果没有正确地设置和正确地分析一个流体问

题，它一般都是不会收敛的。

如果所输入的初始参数和所有的边界条件都是正确的，则当所有变量的收敛监测

量都停止增长，以及所有求解量的平均、最大、最小值都不再升降时，求解过程就算

是完成了。然而，这并不能保证所求解的结果是唯一正确的，因为自然界本身并不保

证存在唯一解。振荡问题（例如：柱体绕流的旋涡脱落问题）用稳态或瞬态求解技术

都不能得到一个稳定的解。要验证一个分析是稳定的或是振荡的，可以通过对它执行

大量的迭代求解来实现。

ANSYS将求解变量的平均、最小、最大值保存在文件Jobname.RFL中，该文件同

时还保存了FLOTRAN的输入数据和计算出的收敛监测量、所有自由度的结果总结、

层流特性和有效特性。可用下面的方式来规定ANSYS进行结果总结的频率：命令：FLDATA5,OUTP,SUMF,value

菜单：Main Menu>Preprocessor>FLOTRAN SetUp>Additional Out>RFL Out Derived Main Menu>Solution>FLOTRAN SetUp>Additional Out>RFL Out Derived

验证结果

验证求解结果的可靠性是所有分析人员的责任，如果一个FLOTRAN分析得到了

非预期的结果，则应进行下列所示的一些操作，这些操作的大部分都可以在开始一个

分析前完成。即使只进行了零次迭代，ANSYS也会生成一个Jobname.RFL文件并检查

所有的输入数据。

1. 检查作为结果总结的一部分而打印出来的质量平衡情况。内部检查将确定是

否有任何的可能会通过模型的质量流，允许质量流的边界条件是：

?确定的速度边界条件

?确定的压力边界条件

?未确定的边界(这有可能是由于用户忘了施加边界条件而致) ANSYS会将进口和出口边界编号列表，而这些应与所希望的条件相对应。

2. 在ANSYS里检查边界条件，以保证其正确性。

3. 检查所定义的流体性质及其随温度的可变动性正确与否，这可在.RFL文件中

方便地检查。

4. 检查用以建立模型的单位制与用以定义流体性质的单位制是否一致。

5. 有时，还需确认与所选选项相联系的方程的求解是否正确（例如：可压缩流

中的压力方程）。

6. 如果求解发散，可能的原因还有：有限元网格不够精细、或者邻近出口处流

场梯度太大，要解决这些问题，可以使用一些诸如惯性松弛等有助于收敛的

手段，本手册的后面将详述各种松弛技术。

7. 如果仅仅只有某个特定的量产生发散，则可将该量重新初始化到一个单值，

并作重启动分析，方式如下：

命令：FLDATA29,MODV

菜单：Main Menu>Preprocessor>FLOTRAN SetUp> Mod Res/Quad

Ord>Modify Results

Main Menu>Solution>FLOTRAN SetUp> Mod Res/Quad Ord>

Modify Results

第三章FLOTRAN设置命令

本章有如下设置

1 FLOTRAN求解控制命令15

2 FLOTRAN执行及输出控制命令?稳态控制参数设置16

3 FLOTRAN执行及输出控制命令?瞬态控制参数设置17

4 FLOTRAN输出及保存文件控制21

5 FLOTRAN输出收敛监测量的控制24

6 FLOTRAN流体类型及其特性的可变性控制25

7 FLOTRAN流体性质参数设定27

8 FLOTRAN分析的坐标系统的指定29

9 FLOTRAN分析参考条件的设置30

10 指定FLOTRAN分析的旋转坐标系统31

11 指定FLOTRAN分析的重力加速度32

12 设置有助于FLOTRAN求解稳定的参数33

13 设定FLOTRAN自由度松弛系数34

14 设定FLOTRAN流体性质松弛因子35

15 设置FLOTRAN分析的自由度限值36

16 选择FLOTRAN各自由度相应的求解器38

17 对FLOTRAN各求解器的控制39

18 设置FLOTRAN湍流模型的一些常数42

19 重新设定FLOTRAN各分析参数的值45

20 控制FLOTRAN面积积分的阶次48

21 FLOTRAN多组份疏运分析的设置及控制49

22 定义FLOTRAN的重启动（续算）控制55

23 设置并执行一个零迭代FLOTRAN分析58

1 FLOTRAN求解控制命令

命令：FLDATA1, SOLU, Lable, Value

功能：设定求解控制

菜单：Main Menu>Preprocessor>FLOTRAN Set Up>Solution Options Main Menu>Solution>FLOTRAN Set Up>Solution Option

其中，Lable的选项及其各自含义如下：

TRAN：求解稳态或非稳态选项，缺省为稳态

FLOW：是否求解流动方程选项，缺省为求解流动方程

TEMP：是否求解温度方程选项，缺省为绝热(不求解温度方程) TURB：层流或湍流选项，缺省为层流

COMP：不可压缩或可压缩选项，缺省为不可压缩流

SWRL：求解轴对称旋流选项，缺省为否

SPEC：求解多组份疏运选项，缺省为否

Value的值为“是”或“否”二者之一。

该命令的菜单形式如下：

2 FLOTRAN执行及输出控制命令 稳态控制参数设置

命令：FLDATA2,ITER, Lable, Value

功能：设置稳态流的迭代及输出控制

其中， Lable的选项及其各自含义如下：

EXEC：总体迭代次数，缺省为10次

OVER：结果文件覆盖频率(每隔Value次迭代)，缺省为0

APPE：结果文件附加频率(每隔Value次迭代)，缺省为0

Value的值以上Lable各自的迭代次数

命令：FLDATA3,TERM, Lable, Value

功能：设置压力和温度的收敛准则

其中， Lable的选项及其各自含义如下：

PRES：设置压力收敛准则，缺省为10-8

TERM：设置温度收敛准则，缺省为10-8

Value为以上两项各自的收敛值

命令：FLDATA5, OUTP, Lable, Value

功能：设置输出控制

其中， Lable的选项及其各自含义如下：

SUMF：输出总结的频率(每隔Value次迭代)，缺省为10

Value为上面SUMF的迭代数

以上三条命令的菜单路径如下：

菜单：Main Menu>Preprocessor>FLOTRAN Set Up>Execution Ctrl

Main Menu>Solution>FLOTRAN Set Up>Execution Ctrl

以上三条命令的菜单形式如下：

3 FLOTRAN执行及输出控制命令 瞬态控制参数设置

命令：FLDATA4, TIME, Lable, Value

功能：基于瞬态时间来设置瞬态分析的求解及输出控制

命令：FLDATA4A, STEP, Lable, Value

功能：基于时间步来设置瞬态分析的输出控制

菜单：Main Menu>Preprocessor>FLOTRAN Set Up> Execution Ctrl

Main Menu>Solution>FLOTRAN Set Up>Execution Ctrl

在点取上面的菜单路径时，首先会弹出一个下图所示的对话框，该对话框的第二项即是激活上面的第二条命令，即用户需要选择到底是基于瞬态分析的时间值、或是时间步、或是二者兼顾来对求解及输出进行控制。同时，该对话框的第一项是用以定义到底是由程序还是由用户自己来对瞬态分析的时间步长进行控制，也就是下面将要讲到的STEP项。这两个选项的值直接影响了随后的弹出菜单的内容，本手册在此是合起来讲的，谨请注意。

第一条命令的 Lable选项及其各自含义如下：

STEP：定义时间步长，当在上面的对话框中设置为用户自定义时，该选项的Value值应为一正值，表示真正的时间步长，命令输入就在Value处输一

真正时间即可；当在上面的对话框中设置为Advection时，即为由程序

定义时间步长，命令输入时可在Value处输-1、-2、-3或-4四值之一，其

含义分别如下：

-1：时间步长会小到在单一的任何一个时间步长内，流场中任意一点的运动距离都不会大于一个单元的长度

-2：只用于可压缩流分析，它使时间步长会小到在单一的任何一个时间步长内，流场中压力信号的传输距离都不会大于一个单元的长度-3：只用于可压缩流分析，它使时间步长取上面-1和-2项中的最小值

-4：只用于纯传热分析（不计算流场方程），它使时间步长会小到在单一的任何一个时间步长内，任何一个“热点”的传导或对流距离

都不会大于一个单元的长度

ISTEP：指定初始时间步长，仅在时间步长由程序定义时有效，缺省为0 NUMB：指定时间步的数目，缺省为10

TEND：指定瞬态分析的终止时间，缺省为1×10

GLOB：指定每一时间步的总体迭代数，缺省为20

PRES：指定每一时间步的压力收敛准则，缺省为1×10-6

TEMP：指定每一时间步的温度收敛准则，缺省为1×10-6

OVER：指定结果文件的覆盖频率 (每隔Value时间覆盖一次)，缺省为0，即不覆盖

APPE：指定结果的输出频率(每隔Value时间输出一次结果到结果文件中)，缺省为1×106

SUMF：指定结果总结的输出频率(每隔 Value时间输出一次结果总结)，缺省为

1×106

BC：指定瞬态边界条件的变化形式，缺省为阶跃变化

其Value的含义如上面各项所述

第二条命令的 Lable选项及其各自含义如下：

OVER：指定结果文件的覆盖频率(每隔Value个时间步覆盖一次)，缺省为0，即不覆盖

APPE：指定结果的输出频率(每隔 Value个时间步输出一次结果到结果文件

中)，缺省为10

SUMF：指定结果总结的输出频率(每隔Value个时间步输出一次结果总结)，缺省为10

下页所示为该二命令的菜单形式，第一幅图为用户自定义时间步长，第二幅图为程序定义时间步长（上述对话框的第一项设定），同时，该二图所显示的都是设定由时间值和时间步同时控制求解及输出选项时（上述对话框的第二项设为“both”）的菜单情况：

ANSYS非稳态热分析及实例详解解析

本章向读者介绍非稳态热分析的基本知识， 主要包括非稳态热分析的应用、 非稳态热分析的基本步骤。 非稳态导热的基本概念 非稳态热分析的应用 非稳态热分析单元 分析的基本步骤 丄本章案例 钢球非稳态传热过程分析 不同材料金属块水中冷却的非稳态传热过程分析 高温铜导线冷却过程分析 7.1 非稳态热分析概述 物体的温度随时间而变化的导热过程称为非稳态导热。 根据物体温度随着时间的推移而变化的 特性可本章要点 非稳态热分析单兀、

以区分为两类非稳态导热：物体的温度随时间的推移逐渐趋于恒定的值以及物体的温度随时间而作周期性的变化。无论在自然界还是工程实际问题中，绝大多数传热过程都是非稳态的。许多工程实际问题需要确定物体内部的温度场随时间的变化，或确定其内部温度达到某一限定值所需要的时间。例如：在机器启动、停机及变动工况时，急剧的温度变化会使部件因热应力而破坏，因此需要确定物体内部的瞬时温度场；钢制工件的热处理是一个典型的非稳态导热过程，掌握工件中温度变化的速率是控制工件热处理质量的重要因素。再例如，金属在加热炉内加热时，需要确定它在加热炉内停留的时间，以保证达到规定的中心温度。可见，非稳态热分析是有相当大的应用价值的。 ANSYS 11.0 及其相关的下属产品均支持非稳态的热分析。非稳态热分析确定了温度以及其它随时间变化的热参数。 7.1.1 非稳态热分析特性 瞬态热分析用于计算一个系统的随时间变化的温度场及其它热参数。在工程上一般用瞬态热分析计算温度场，并将之作为热载荷进行应力分析。 瞬态热分析的基本步骤与稳态热分析类似。主要的区别是瞬态热分析中的载荷是随时间变化的。为了表达随时间变化的载荷，首先必须将载荷 - 时间曲线分为载荷步。对于每一个载荷步，必须定义载荷值及时间值，同时必须选择载荷步为渐变或阶越。

ansys旋转经典命令流

1 旋转摩擦 (1) 2. 电磁三d命令流实例（论坛看到） (11) 3. 帮助感应加热例子induction heating of a solid cylinder billet (15) 4. 感应加热温度场的数值模拟（论文）inducheat30命令流 (19) 5. 如何施加恒定的角速度？Simwe仿真论坛 (24) 6. 旋转一个已经生成好的物体 (27) 7. 产生这样的磁力线 (28) 8. 旋转摩擦生热简单例子（二维旋转） (32) 8.1. 原版 (32) 8.2. 部分gui操作 (35) 9. VM229 Input Listing (39) 10 轴承---耦合+接触分析 (47) 11. 板的冲压仿真 (52) 1 旋转摩擦 FINISH /FILNAME,Exercise24 !定义隐式热分析文件名 /PREP7 !进入前处理器 ET,1,SOLID5 !选择单元类型 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,7800 !定义材料1的密度 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,C,1,,460 !定义材料1的比热 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,KXX,1,,66.6 !定义材料1的热传导系数 MPTEMP,,,,,,,, MPTEMP,1,0 UIMP,1,REFT,,,30 !定义材料1的热膨胀系数的参考温度 MPDATA,ALPX,1,,1.06e-5 !定义材料1的热膨胀系数MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,206e9 !定义材料1的弹性模量 MPDATA,PRXY,1,,0.3 !定义材料1的泊松比 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,2,,8900 !定义材料2的密度 MPTEMP,,,,,,,, MPTEMP,1,0

ANSYS热分析指南与经典案例

第一章简介 一、热分析的目的 热分析用于计算一个系统或部件的温度分布及其它热物理参数，如热量的获取或损失、热梯度、热流密度（热通量〕等。 热分析在许多工程应用中扮演重要角色，如内燃机、涡轮机、换热器、管路系统、电子元件等。 二、ANSYS的热分析 ?在ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/FLOTRAN、ANSYS/ED五种产品中包含热分析功能,其中 ANSYS/FLOTRAN不含相变热分析。 ?ANSYS热分析基于能量守恒原理的热平衡方程，用有限元法计算各节点的温度，并导出其它热物理参数。 ?ANSYS热分析包括热传导、热对流及热辐射三种热传递方式。此外，还可以分析相变、有内热源、接触热阻等问题。 三、ANSYS 热分析分类 ?稳态传热：系统的温度场不随时间变化 ?瞬态传热：系统的温度场随时间明显变化 四、耦合分析 ?热－结构耦合 ?热－流体耦合 ?热－电耦合 ?热－磁耦合 ?热－电－磁－结构耦合等

第二章 基础知识 一、符号与单位 W/m 2-℃ 3 二、传热学经典理论回顾 热分析遵循热力学第一定律，即能量守恒定律: ● 对于一个封闭的系统（没有质量的流入或流出〕 PE KE U W Q ?+?+?=- 式中: Q —— 热量; W —— 作功; ?U ——系统内能; ?KE ——系统动能； ?PE ——系统势能； ● 对于大多数工程传热问题：0＝＝PE KE ??； ● 通常考虑没有做功：0=W , 则：U Q ?=； ● 对于稳态热分析：0=?=U Q ，即流入系统的热量等于流出的热量； ● 对于瞬态热分析：dt dU q = ，即流入或流出的热传递速率q 等于系统内能的变化。 三、热传递的方式 1、热传导 热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温度梯度而引起的内能的交换。热传导遵循付里叶定律：dx dT k q -=''，式中''q 为热流

ansys命令流

第一天目标： 熟悉ANSYS基本关键字的含义k --> Keypoints关键点l --> Lines线a --> Area 面v --> Volumes体e --> Elements单元n --> Nodes节点cm --> component组元et --> element type单元类型mp --> material property材料属性r --> real constant实常数d --> DOF constraint约束f --> Force Load集中力sf --> Surface load on nodes 表面载荷bf --> Body Force on Nodes体载荷ic --> Initial Conditions初始条件第二天目标： 了解命令流的整体结构,掌握每个模块的标识!文件说明段/BATCH/TILE,test analysis!定义工作标题/FILENAME,test!定义工作文件名/PREP7!进入前处理模块标识!定义单元,材料属性,实常数段ET,1,SHELL63!指定单元类型ET,2,SOLID45!指定体单元MP,EX,1,2E8!指定弹性模量MP,PRXY,1, 0.3!输入泊松比MP,DENS,1, 7.8E3!输入材料密度R,1, 0.001!指定壳单元实常数-厚度......!建立模型K,1,0,0,,!定义关键点 K,2,50,0,,K,3,50,10,,K,4,10,10,,K,5,10,50,,K,6,0,50,,A,1,2,3,4,5,6,!由关键点生成面......!划分网格ESIZE,1,0,AMESH, 1......FINISH!前处理结束标识/SOLU!进入求解模块标识!施加约束和载荷DL,5,,ALLSFL,3,PRES,1000SFL,2,PRES, 1000......SOLVE!求解标识FINISH!求解模块结束标识/POST1!进入通用后处理器标识....../POST26!进入时间历程后处理器……/EXIT,SAVE!退出并存盘以下是日志文件中常出现的一些命令的标识说明,希望能给大家在整理LOG文件时有所帮助/ANGLE!指定绕轴旋转视图/DIST!说明对视图进行缩放/DEVICE!设置图例的显示,如： 风格,字体等/REPLOT!重新显示当前图例/RESET!恢复缺省的图形设置/VIEW!设置观察方向/ZOOM!对图形显示窗口的某一区域进行缩放第三天生成关键点和线部分 1.生成关键点K,关键点编号,X坐标,Y坐标,Z坐标例：

ANSYS稳态热分析的基本过程和实例

ANSYS稳态热分析的基本过程 ANSYS热分析可分为三个步骤： ?前处理：建模、材料和网格 ?分析求解：施加载荷计算 ?后处理：查看结果 1、建模 ①、确定jobname、title、unit； ②、进入PREP7前处理，定义单元类型，设定单元选项； ③、定义单元实常数； ④、定义材料热性能参数，对于稳态传热，一般只需定义导热系数，它可 以是恒定的，也可以随温度变化； ⑤、创建几何模型并划分网格，请参阅《ANSYS Modeling and Meshing Guide》。 2、施加载荷计算 ①、定义分析类型 ●如果进行新的热分析： Command: ANTYPE, STATIC, NEW GUI: Main menu>Solution>-Analysis Type->New Analysis>Steady-state ●如果继续上一次分析，比如增加边界条件等： Command: ANTYPE, STATIC, REST GUI: Main menu>Solution>Analysis Type->Restart ②、施加载荷 可以直接在实体模型或单元模型上施加五种载荷(边界条件) ： a、恒定的温度 通常作为自由度约束施加于温度已知的边界上。 Command Family: D GUI：Main Menu>Solution>-Loads-Apply>-Thermal-Temperature b、热流率 热流率作为节点集中载荷，主要用于线单元模型中(通常线单元模型不能施加对流或热流密度载荷)，如果输入的值为正，代表热流流入节点，即单元获取热量。如果温度与热流率同时施加在一节点上则ANSYS读取温度值进行计算。 注意：如果在实体单元的某一节点上施加热流率，则此节点周围的单元要

Ansys 第 例瞬态热分析实例一水箱

第33例瞬态热分析实例——水箱 本例介绍了利用ANSYS进行瞬态热分析的方法和步骤、瞬态热分析时材料模型所包含的内容，以及模型边界条件和初始温度的施加方法。 33.1概述 热分析是计算热应力的基础，热分析分为稳态热分析和瞬态热分析，稳态热分析将在后面两个例子中介绍，本例介绍瞬态热分析。 33.1.1 瞬态热分析的定义 瞬态热分析用于计算系统随时间变化的温度场和其他热参数。一般用瞬态热分析计算温度场，并找到温度梯度最大的时间点，将此时间点的温度场作为热载荷来进行应力计算。 33.1.2 嚼态热分析的步骤 瞬态热分析包括建模、施加载荷和求解、查看结果等几个步骤。 1．建模 瞬态热分析的建模过程与其他分析相似，包括定义单元类型、定义单元实常数、定义材料特性、建立几何模型和划分网格等。 注意：瞬态热分析必须定义材料的导热系数、密度和比热。 2.施加载荷和求解 (1)指定分析类型， Main Menu→Solution→Analysis Type→New Analysis,选择 Transient。 (2)获得瞬态热分析的初始条件。 定义均匀的初始温度场：Main Menu→Solution→Define Loads→Settings→Uniform Temp，初始温度仅对第一个子步有效，而用Main Menu →Solution→Define Loads→Apply→Thermal→Temperature命令施加的温

度在整个瞬态热分析过程中均不变，应注意二者的区别。 定义非均匀的初始温度场：如果非均匀的初始温度场是已知的，可以用Main Menu→Solution→Define Loads→Apply→Initial Condit'n→Define 即IC命令施加。非均匀的初始温度场一般是未知的，此时必须先进行行稳态分析确定该温度场。该稳态分析与一般的稳态分析相同。 注意：要设定载荷（如已知的温度、热对流等），将时间积分关闭，选择Main Menu→Solution→Load Step Opts→Time/Frequenc→Time Integration→Amplitude Decay;设定只有一个子步，时间很短（如(0.01s)的载荷步， Main Menu→Solution→Load Step Opts→Time/Frequenc→Time →Time Step。 (3)设置载荷步选项。 普通选项包括每一载荷步结束的时间、每一载荷步的子步数、阶跃选项等，选择Main Menu→Solution→Load Step Opts→Time/Frequenc→Time-Time Step. 非线性选项包括：迭代次数（默认25），选择Main Menu→Solution→Load Step Opts→Nonlinear→Equilibrium Iter;打开自动时间步长，选择Main Menu→Solution→Load Step Opts→Time/Frequenc→Time→Time Step：将时间积分打开，选择Main Menu→Solution→Load Step Opts→Time/Frequenc→Time Integration→Amplitude Decay. 输出选项包括：控制打印的输出，选择Main Menu→Solution→Load Step Opts→Output Ctrls→Solu Printout; 结果文件的输出，选择Main Menu →Solution→Load Step Opts→Output Ctrls→DB/Results File.

ANSYS-结构稳态(静力)分析之经典实例-命令流格式

ANSYS 结构稳态(静力)分析之经典实例-命令流格式.txt两人之间的感情就像织毛衣，建立 的时候一针一线，小心而漫长，拆除的时候只要轻轻一拉。。。。/FILNAME,Allen-wrench,1 ! Jobname to use for all subsequent files /TITLE,Static analysis of an Allen wrench /UNITS,SI ! Reminder that the SI system of units is used /SHOW ! Specify graphics driver for interactive run; for batch ! run plots are written to pm02.grph ! Define parameters for future use EXX=2.07E11 ! Young's modulus (2.07E11 Pa = 30E6 psi) W_HEX=.01 ! Width of hex across flats (.01m=.39in) *AFUN,DEG ! Units for angular parametric functions定义弧度单位 W_FLAT=W_HEX*TAN(30) ! Width of flat L_SHANK=.075 ! Length of shank (short end) (.075m=3.0in) L_HANDLE=.2 ! Length of handle (long end) (.2m=7.9 in) BENDRAD=.01 ! Bend radius of Allen wrench (.01m=.39 in) L_ELEM=.0075 ! Element length (.0075 m = .30 in) NO_D_HEX=2 ! Number of divisions on hex flat TOL=25E-6 ! Tolerance for selecting nodes (25e-6 m = .001 in) /PREP7 ET,1,SOLID45 ! 3维实体结构单元；Eight-node brick element ET,2,PLANE42 ! 2维平面结构；Four-node quadrilateral (for area mesh) MP,EX,1,EXX ! Young's modulus for material 1；杨氏模量 MP,PRXY,1,0.3 ! Poisson's ratio for material 1；泊松比 RPOLY,6,W_FLAT ! Hexagonal area创建规则的多边形 K,7 ! Keypoint at (0,0,0) K,8,,,-L_SHANK ! Keypoint at shank-handle intersection K,9,,L_HANDLE,-L_SHANK ! Keypoint at end of handle L,4,1 ! Line through middle of hex shape L,7,8 ! Line along middle of shank L,8,9 ! Line along handle LFILLT,8,9,BENDRAD ! Line along bend radius between shank and handle! 产生 一个倒角圆，并生成三个点 /VIEW,,1,1,1 ! Isometric view in window 1 /ANGLE,,90,XM ! Rotates model 90 degrees about X! 不用累积的旋转 /TRIAD,ltop /PNUM,LINE,1 ! Line numbers turned on LPLOT

几个ansys经典实例(长见识)

平面问题斜支座的处理 如图5-7所示，为一个带斜支座的平面应力结构，其中位置2及3处为固定约束，位置4处为一个45o的斜支座，试用一个4节点矩形单元分析该结构的位移场。 (a)平面结构(b)有限元分析模型 图5-7 带斜支座的平面结构 基于ANSYS平台，分别采用约束方程以及局部坐标系的斜支座约束这两种方式来进行处理。 (7) 模型加约束 左边施加X,Y方向的位移约束 ANSYS Main Menu: Solution →Define Loads →Apply →-Structural→Displacement On Nodes →选取2，3号节点→OK →Lab2: All DOF(施加X,Y方向的位移约束) →OK 以下提供两种方法处理斜支座问题，使用时选择一种方法。 ?采用约束方程来处理斜支座 ANSYS Main Menu：Preprocessor →Coupling/ Ceqn →Constraint Eqn ：Const :0, NODE1:4, Lab1: UX,C1:1,NODE2:4,Lab2:UY,C2:1→OK 或者?采用斜支座的局部坐标来施加位移约束 ANSYS Utility Menu：WorkPlane →Local Coordinate System →Create local system →At specified LOC + →单击图形中的任意一点→OK →XC、YC、ZC分别设定为2，0，0,THXY：45 →OK ANSYS Main Menu：Preprocessor →modeling →Move / Modify →Rotate Node CS →To active CS → 选择4号节点 ANSYS Main Menu：Solution →Define Loads →Apply →Structural →Displacement On Nodes →选取4号节点→OK →选择Lab2:UY(施加Y方向的位移约束) →OK 命令流； !---方法1 begin----以下的一条命令为采用约束方程的方式对斜支座进行处理 CE,1,0,4,UX,1,4,UY,-1 !建立约束方程(No.1): 0=node4_UX*1+node_UY*(-1) !---方法1 end --- !--- 方法2 begin --以下三条命令为定义局部坐标系，进行旋转，施加位移约束 !local,11,0,2,0,0,45 !在4号节点建立局部坐标系 !nrotat, 4 !将4号节点坐标系旋转为与局部坐标系相同 !D,4,UY !在局部坐标下添加位移约束 !--- 方法2 end

一个经典的ansys热分析实例(流程序)

/PREP7 /TITLE,Steady-state thermal analysis of pipe junction /UNITS,BIN ! 英制单位；Use U. S. Customary system of units (inches) ! /SHOW, ! Specify graphics driver for interactive run ET,1,90 ! Define 20-node, 3-D thermal solid element MP,DENS,1,.285 ! Density = .285 lbf/in^3 MPTEMP,,70,200,300,400,500 ! Create temperature table MPDATA,KXX,1,,8.35/12,8.90/12,9.35/12,9.80/12,10.23/12 ! 指定与温度相对应的数据材料属性；导热系数；Define conductivity values MPDATA,C,1,,.113,.117,.119,.122,.125 ! Define specific heat values（比热） MPDATA,HF,2,,426/144,405/144,352/144,275/144,221/144 ! Define film coefficient;除144是单位问题，上面的除12也是单元问题 ! Define parameters for model generation RI1=1.3 ! Inside radius of cylindrical tank RO1=1.5 ! Outside radius Z1=2 ! Length RI2=.4 ! Inside radius of pipe RO2=.5 ! Outside pipe radius Z2=2 ! Pipe length CYLIND,RI1,RO1,,Z1,,90 ! 90 degree cylindrical volume for tank WPROTA,0,-90 ! 旋转当前工作的平面；从Y到Z旋转-90度；；Rotate working plane to pipe axis CYLIND,RI2,RO2,,Z2,-90 ! 角度选择在了第四象限；90 degree cylindrical volume for pipe WPSTYL,DEFA ! 重新安排工作平面的设置；另外WPSTYL,STAT to list the status of the working plane；；Return working plane to default setting BOPT,NUMB,OFF ! 关掉布尔操作的数字警告信息；Turn off Boolean numbering warning VOVLAP,1,2 ! 交迭体；Overlap the two cylinders /PNUM,VOLU,1 ! 体编号打开；Turn volume numbers on /VIEW,,-3,-1,1

ANSYS_结构稳态(静力)分析之经典实例-命令流格式

/FILNAME,Allen-wrench,1 ! Jobname to use for all subsequent files /TITLE,Static analysis of an Allen wrench /UNITS,SI ! Reminder that the SI system of units is used /SHOW ! Specify graphics driver for interactive run; for batch ! run plots are written to pm02.grph ! Define parameters for future use EXX=2.07E11 ! Young's modulus (2.07E11 Pa = 30E6 psi) W_HEX=.01 ! Width of hex across flats (.01m=.39in) *AFUN,DEG ! Units for angular parametric functions定义弧度单位 W_FLAT=W_HEX*TAN(30) ! Width of flat L_SHANK=.075 ! Length of shank (short end) (.075m=3.0in) L_HANDLE=.2 ! Length of handle (long end) (.2m=7.9 in) BENDRAD=.01 ! Bend radius of Allen wrench (.01m=.39 in) L_ELEM=.0075 ! Element length (.0075 m = .30 in) NO_D_HEX=2 ! Number of divisions on hex flat TOL=25E-6 ! Tolerance for selecting nodes (25e-6 m = .001 in) /PREP7 ET,1,SOLID45 ! 3维实体结构单元；Eight-node brick element ET,2,PLANE42 ! 2维平面结构；Four-node quadrilateral (for area mesh) MP,EX,1,EXX ! Young's modulus for material 1；杨氏模量 MP,PRXY,1,0.3 ! Poisson's ratio for material 1；泊松比 RPOLY,6,W_FLAT ! Hexagonal area创建规则的多边形 K,7 ! Keypoint at (0,0,0) K,8,,,-L_SHANK ! Keypoint at shank-handle intersection K,9,,L_HANDLE,-L_SHANK ! Keypoint at end of handle L,4,1 ! Line through middle of hex shape L,7,8 ! Line along middle of shank L,8,9 ! Line along handle LFILLT,8,9,BENDRAD ! Line along bend radius between shank and handle! 产生一个倒角圆，并生成三个点 /VIEW,,1,1,1 ! Isometric view in window 1 /ANGLE,,90,XM ! Rotates model 90 degrees about X! 不用累积的旋转 /TRIAD,ltop /PNUM,LINE,1 ! Line numbers turned on LPLOT ! Line numbers off !

ansys热分析例题

问题描述：一个30公斤重、温度为70℃的铜块,以及一个20公斤重、温度为80℃的铁块，突然放入温度为20℃、盛满了300升水的、完全绝热的水箱中,如图所示。过了一个小时，求铜块与铁块的最高温度(假设忽略水的流动)。 材料热物理性能如下：热性能单位制 铜铁水 导热系数W/m℃ 383 37 密度Kg/m 8889 7833 996 比热J/kg℃ 390 448 4185 菜单操作过程： 一、设置分析标题 1、选择“Utility Menu>File>Change Jobname”，输入文件名Transient1。 2、选择“Utility Menu>File>Change Title”输入Thermal Transient Exercise 1。 二、定义单元类型 1、选择“Main Menu>Preprocessor”，进入前处理。 2、选择“Main Menu>Preprocesor>Element Type>Add/Edit/Delete”。选择热平面单元plane77。 三、定义材料属性 1、选择“Main Menu>Preprocessor>Material Props>Material Models”，在弹出的材料定义窗口中顺序双击Thermal选项。 2、点击Conductivity，Isotropic，在KXX框中输入383；点击Density，在DENS框中输入8898；点击Specific Heat，在C框中输入390。 3、在材料定义窗口中选择Material>New Model，定义第二种材料。 4、点击Conductivity，Isotropic，在KXX框中输入70；点击Density，在DENS框中输入7833；点击Specific Heat，在C框中输入448。 5、在材料定义窗口中选择Material>New Model，定义第三种材料。 6、点击Conductivity，Isotropic，在KXX框中输入.61；点击Density，在DENS框中输入996；点击Specific Heat，在C框中输入4185。 四、创建几何模型 1、选择“Main Menu>Preprocessor>-Modeling->Create>-Areas->Retangle>By Dimensions”，输入X1=0, Y1=0, X2=, Y2=, 点击Apply；输入X1=, Y1=, X2= ,Y2=, 点击Apply；输入X1= Y1=, X2= Y2=+, 选择OK。 2、选择“Main Menu>Preprocessor>-Modeling->Operate>Booleans>Overlap”，选择Pick All。 3、选择“Utility Menu>Plotctrls>Numbering>Areas, on”。 4、选择“Utility Menu>Plot>Areas”。 五、划分网格 1、选择“Main Menu>Preprocessor>-Attributes->Define->All Areas”，选择材料1。 2、选择“Main Menu>Preprocessor>Meshing->Size Cntrls->-Manualsize->-Global->Size”，输入单元大小。 3、选择“Main Menu>Preprocessor>Meshing->Mesh->-Areas->Mapped>3 or 4 sided”，选择铜块。 4、选择“Main Menu>Preprocessor>-Attributes->Define->All Areas”，选择材料2。 5、选择“Main Menu>Preprocessor>Meshing->Mesh->-Areas->Mapped>3 or 4 sided”，选

命令流ansys经典实例

ansys钢筋混凝土建模实例 finish /clear /prep7 et,1,solid65,,,,,,,1 et,2,link8 et,3,185 et,4,solid45 !************************定义材料属性***************************** !混凝土材料属性 mp,ex,1,1.596e10 mp,prxy,1,0.2 mp,dens,1,2400 fc=1.68e7 !c30混凝土轴心抗压强度设计值 ft=1.86e6 tb,concr,1 tbdata,,0.5,0.95,ft,-1 tb,miso,1,,11 tbpt,,0.0002,fc*0.19 tbpt,,0.0004,fc*0.36 tbpt,,0.0006,fc*0.51 tbpt,,0.0008,fc*0.64 tbpt,,0.0010,fc*0.75 tbpt,,0.0012,fc*0.84 tbpt,,0.0014,fc*0.91 tbpt,,0.0016,fc*0.96 tbpt,,0.0018,fc*0.99 tbpt,,0.002,fc tbpt,,0.0033,fc !57到68 v,6,7,56,54,10,11,60,58 vgen,4,57,57,1,0,0.38,0 vgen,2,57,60,1,0.25,0,0 vgen,2,57,60,1,-0.25,0,0 !69到80 v,189,190,234,233,191,192,238,237 vgen,4,69,69,1,0,0.38,0

vgen,3,69,72,1,0,0,0.25 !121到136 vgen,2,105,120,1,0,0,0.107 !193到208 k,1178,0.81,0,0.107 k,1179,0.823,0,0.107 k,1180,0.81,0,0.417 k,1181,0.81,0.16,0.107 k,1182,0.823,0.16,0.107 k,1183,0.81,0.16,0.417 k,1184,0.81,0.38,0.107 k,1185,0.823,0.38,0.107 k,1186,0.81,0.38,0.417 v,1178,1179,187,1180,1181,1182,189,1183 v,1181,1182,189,1183,1184,1185,191,1186 vgen,4,193,194,1,0,0.38,0 vgen,2,193,200,1,0.107,0,0 !209到232 k,1271,-0.02,0,-0.013 k,1272,-0.02,0,0 k,1273,-0.02,0.16,-0.013 k,1274,-0.02,0.16,0 k,1275,-0.02,0.38,-0.013 k,1276,-0.02,0.38,0 v,677,1,1272,1271,679,5,1274,1273 v,679,5,1274,1273,681,9,1276,1275 vgen,4,209,210,1,0,0.38,0 vgen,3,209,216,1,0.25,0,0 !233到256 vgen,2,209,232,1,0,0,0.107 vgen,2,257,280,1,0.107,0,0 k,5001,-0.1,1.3,0 k,5002,-0.02,1.3,0 k,5003,-0.02,1.3,0.094

ANSYS热分析指南

ANSYS热分析指南 第一章简介 1.1热分析的目的 热分析用于计算一个系统或部件的温度分布及其它热物理参数，我们一般关心的参数有： 温度的分布 热量的增加或损失 热梯度 热流密度 热分析在许多工程应用中扮演着重要角色，如内燃机、涡轮机、换热器、管路系统、电子元件等等。通常在完成热分析后将进行结构应力分析，计算由于热膨胀或收缩而引起的热应力。 1.2ANSYS中的热分析 ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Professional、 ANSYS/FLOTRAN四种产品中支持热分析功能。ANSYS热分析基于由能量守恒原理导出的热平衡方程，有关细节，请参阅《ANSYS Theory Reference》。ANSYS使用有限元法计算各节点的温度，并由其导出其它热物理参数。 ANSYS可以处理所有的三种主要热传递方式：热传导、热对流及热辐射。1.2.1对流 热对流在ANSYS中作为一种面载荷，施加于实体或壳单元的表面。首先需要输入对流换热系数和环境流体温度，ANSYS将计算出通过表面的热流量。如果对流换热系数依赖于温度，可以定义温度表，以及在每一个温度点处的对流换热系数。 1.2.2辐射 ANSYS提供了四种方法来解决非线性的辐射问题： 辐射杆单元（LINK31） 使用含热辐射选项的表面效应单元（SURF151-2D,或SURF152-3D）

在AUX12中，生成辐射矩阵，作为超单元参与热分析 使用Radiosity求解器方法 有关辐射的详细描述请阅读本指南第四章。 1.2.3特殊的问题 除了前面提到的三种热传递方式外，ANSYS热分析还可以解决一些诸如：相变（熔融与凝固）、内部热生成（如焦耳热）等的特殊问题。例如，可使用热质点单元MASS71模拟随温度变化的内部热生成。 1.3热分析的类型 ANSYS支持两种类型的热分析： 1．稳态热分析确定在稳态的条件下的温度分布及其他热特性，稳态条件指热量随时间的变化可以忽略。 2．瞬态热分析则计算在随时间变化的条件下，温度的分布和热特性。 1.4耦合场分析 ANSYS中可与热分析进行耦合的方式有热—结构、热－电磁等。耦合场分析可以使用ANSYS中的矩阵耦合单元，或者在独立的物理环境中使用序惯荷载耦合。有关耦合场分析的详细描述，请参阅《ANSYS Coupled-Field Analysis Guide》。 1.5关于菜单路径和命令语法 在本指南中，您将会看到相关的ANSYS命令及其等效的菜单路径。这些参考的命令仅仅包括命令名，因为并不总是需要指定所有的参数，而且不同的参数组合会有不同的作用。有关ANSYS命令的更多的叙述，请参考《ANSYS Commands Reference》。 菜单路径将近可能完整得列出。对于多数情况，选择菜单就能够完成所需要的功能；但还有一些情况，选择文中所示菜单后会弹出一个菜单或是对话框，由此定义其他的选项来执行一些特定的任务。 第二章基础知识 2.1符号与单位

实例ANSYS

ANSYS应用实例—圆轴扭转分析 摘要：本文建立了圆轴的模型，并对其施加了扭转应力，应用ANSYS软件，对圆轴的扭转进行了分析，分析表明，ANSYS的分析结果与实际理论计算的结果是十分接近的。 关键词：圆轴扭转ANSYS 问题描述:设等直圆轴截面直径D=50mm，长度120mm，作用在圆轴两端上的转矩M=1.5X103N.m，试分析该圆轴的扭转情况。 分析步骤： (1)改变工作名：拾取菜单Utility Menu→File→Change Jobname，在弹出的文本框中输入适当的文件名，单击OK即可。 (2) 创建单元类型：拾取菜单Utility Menu→Preprocessor→Element Type→ Add/Edit/Delete。在弹出的对话框中单击Add按钮，在左侧列表中选Structural Solid，在右侧列表中选Quad 8node 183，单击Apply按钮，再在右侧列表中选择Brick 20node 186，单击OK按钮，然后单击对话框中的Close按钮。 图1 单元类型对话框 （3 ）定义材料类型：拾取菜单Utility Menu→Preprocessor→Material Props→Material Model。弹出如图2所示的对话框，在右侧列表中一次双击Structural、Linear、Elastic、Isotropic，弹出如图3所示的对话框，在EX文本框中输入弹性模量2.08e11，在PRXY中输入泊松比0.3，单击OK，关闭对话框即可。 图3 材料特性对话框

（4 ）创建矩形面：拾取菜单Utility Menu→Preprocessor→Modeling→Create→Areas→Rectangle→By Dimensions。弹出如图4所示的对话框，在“X1，X2“文本框中输入0,0.025，在”Y1，Y2”文本框中输入0，0.12，单击OK按钮即可。 图4 创建矩形面对话框 （5 ）划分单元：拾取菜单Utility Menu→Preprocessor→Meshing→Mesh Tool。弹出如图6所示的对话框，单击Size Controls区域中的Lines后的Set按钮，弹出拾取窗口，拾取矩形面的任意短边，单击OK，在如图7中的弹出对话框在NDIV 文本框中输入5，单击Apply按钮，在此弹出拾取窗口，拾取矩形面的任一长边，单击OK按钮，再次弹出对话框，在NDIV文本框中输入8，单击OK按钮。在Mesh区域，选择单元形状为Quad，选择划分单元的方法为Mapped，单击Mesh 按钮，弹出拾取窗口，拾取面，单击OK，然后关闭窗口即可。 图5 划分单元工具对话框

Ansys经典例子：如何使用Ansys帮助文件

Ansys经典例子：如何使用Ansys帮助文件 很多网友都曾觉得ANSYS使用起来有一定的难度，经常会遇到这样或那样的问题，但市面上的参考书又不尽如人意，那究竟有没有比较好的参数书? 有的，个人认为ANSYS的帮助文件就是一本不错的参数书。接下来就ANSYS在线帮助的使用做一些基本的介绍，希望能对初学者有所帮助。 ANSYS的帮助文件包括所有ANSYS命令解释及所有的GUI解释，还包括ANSYS各模块的分析指南，实例练习等。 一.进入帮助系统 可以通过下列三种方式进入: 1.进入ANSYS的操作界面后，在应用菜单中选取Help进入; 2.在ANSYS程序组中选取Help System进入：Start Menu > Programs >ANSYS XX>Help System; 3.在任何对话框中选取Help。 二.帮助系统的内容安排： 点击帮助系统的目录，就看到如下的ANSYS帮助系统的整体内容安排： 1.前面4个部分是与软件版本，安装，注册相关的信息，只需作相应的了解即可，如下： ※Release Notes ※ANSYS Installation and Configuration Guide for UNIX ※ANSYS Installation and Configuration Guide for Windows ※ANSYS, Inc. Licensing Guide 2.接下来两个部分是比较重要的部分，ANSYS的命令和单元手册，对用到的命令和单元应作详细的了解和掌握。 ※ANSYS Commands Reference ※ANSYS Element Reference 3.下面四个部分是ANSYS相关的操作手册，说明如下： ※Operations Guide 基本界面，操作指南 ※Basic Analysis Procedures Guide 基础分析指南 ※Advanced Analysis Techniques Guide 高级分析指南 ※Modeling and Meshing Guide 建模与分网指南 4.以下几个部分则是ANSYS分模块的分析指南，如下： ※Structural Analysis Guide 结构分析指南 ※Thermal Analysis Guide 热分析指南 ※CFD FLOTRAN Analysis Guide 流体分析指南 ※Electromagnetic Field Analysis Guide 电磁场分析指南 ※Coupled-Field Analysis Guide 耦合场分析指南 5.为更好的使用ANSYS方便，快捷的解决更多的工程实际问题，建议仔细学习以下几个部分： ※APDL Programmer's Guide：APDL操作手册 ※ANSYS Troubleshooting Guide：ANSYS错误信息指南 ※Mechanical Toolbar：机械工具栏 ※ANSYS/LS-DYNA User's Guide：ANSYS/LS-DYNA操作指南 ※ANSYS Connection Users Guide：接口技术指南