关于ansys

常见问题1、如何导出ansys 命令流文本？你的模型是自己建立的还是在其他地方建立的？如果是在ansys里面建立的，就会有你的命令流！大概有三种方式可以提取：1，在工作目录里面后缀为.logd的文件里面自动保存了你的所有操作，当然包括命令流了，可以用记事本打开查看，估计你得有基本的命令流常识，不然看不懂，下同！2，1楼说的也是一种；应该是在Main Menu>Preprocessor>Session Editor这样的GUI路径，打开以后可以使用ctrl+A全选，然后再ctrl+c复制，这样就复制了这里面的命令流，记住要使用键盘哦，复制好了以后就可以建立一个记事本，然后这时候粘贴就行了，保存记事本即可！当然使用另存为也可以！3，还有一种是通过一个命令：命令是：LGWRITE, Fname, Ext, --, KeditGUI是：Utility Menu>File>Write DB Log File （用户名就是Fname，生成文件后缀就是ext,自己定义，比如LGWRITE,me,txt,就生成了文件me.txt，里面包括了你在这个模型已经操作的命令流）这个命令在默认情况下与第一种方法生成的基本上一样，而且是单独的一个案例的命令，在第一种方法里面，如果你没有及时清除你的工作空间里面的所有的文件，可能是命令流的一种累加，包括了你操作了的好几个的案例、也就是例子的命令流，可能需要你去分辨了；而且第三种方法还可以对命令流里面进行简单的筛选，当然是你要熟悉这个命令以后，不然默认就行了！2、什么意思？这里是选择你要使用的是那种ansys的产品，比如是结构的，lsdyna的，电磁的，等等，还有可以设定内存的需要啊，改变工作空间等！比如：ansys mechanical/ls-dyna是说的在一般的机械分析中使用，是有一定的机械分析优势，但是在很多的单元使用中是局限性的ansys mechanical/emag/ls-dyna是说的包括了电磁的一些东西，当然也包括机械方面的ansys mechanical/flotran/ls-dyna这个就偏重流体，也就是CFD方面的综合使用了3、ANSYS文件类型ANSYS在分析过程中需要读写文件，文件格式为jobname.ext，其中jobname是设定的工作文件名文件类型内容plane.db 二进制数据库文件plane.dbb 二进制数据库备份文件（当非线性分析不正常终止时产生）plane.emat 二进制单元矩阵plane.err 文本错误或警告信息plane.esav 二进制单元存储数据（当非线性分析不能向上兼容时产生）plane.full 二进制装配的整体刚度和质量矩阵plane.ldhi 文本载荷步中载荷和边界条件plane.log 文本命令行输入历史记录plane.mntr 二进制监视文件plane.opt 文本优化数据plane.osav 二进制单元存储文件的备份plane.rdb 二进制第一载荷步第一子步起始时的数据状态plane.rst 二进制结构或耦合场分析得到的结果文件文件名.ext是由ANSYS定义的扩展名，用于区分文件的用途和类型，默认的工作文件名是file。

ANSYS分析中有一些特殊的文件，其中主要的几个是数据库文件jobname.db、记录文件jobname.log、输出文件jobname.out、错误文件jobname.err、结果文件jobname.rxx及图形文件jobname.grph2009-06-09 12:104、残差曲线震荡问题一. 残差波动的主要原因：1、高精度格式；2、网格太粗；3、网格质量差；4、流场本身边界复杂，流动复杂；5、模型的不恰当使用。

二. 问：在进行稳态计算时候，开始残差线是一直下降的，可是到后来各种残差线都显示为波形波动，是不是不收敛阿？答：有些复杂或流动环境恶劣情形下确实很难收敛。

计算的精度（2 阶），网格太疏，网格质量太差，等都会使残差波动。

经常遇到，一开始下降，然后出现波动，可以降低松弛系数，我的问题就能收敛，但如果网格质量不好，是很难的。

通常，计算非结构网格，如果问题比较复杂，会出现这种情况，建议作网格时多下些功夫。

理论上说，残差的震荡是数值迭代在计算域内传递遭遇障碍物反射形成周期震荡导致的结果，与网格亚尺度雷诺数有关。

例如，通常压力边界是主要的反射源，换成OUTFLOW 边界会好些。

这主要根据经验判断。

所以我说网格和边界条件是主要因素。

三. 1、网格问题：比如流场内部存在尖点等突变，导致网格在局部质量存在问题，影响收敛。

2、可以调整一下courant number，courant number实际上是指时间步长和空间步长的相对关系，系统自动减小courant数，这种情况一般出现在存在尖锐外形的计算域，当局部的流速过大或者压差过大时出错，把局部的网格加密再试一下。

在fluent中，用courant number来调节计算的稳定性与收敛性。

一般来说，随着courant number的从小到大的变化，收敛速度逐渐加快，但是稳定性逐渐降低。

所以具体的问题，在计算的过程中，最好是把courant number从小开始设置，看看迭代残差的收敛情况，如果收敛速度较慢而且比较稳定的话，可以适当的增加courant number的大小，根据自己具体的问题，找出一个比较合适的courant number，让收敛速度能够足够的快，而且能够保持它的稳定性。

本人觉得可以重点参考第四个回答。

另外，如果出现连续方程残差很高收敛慢的情况，首先应该检查的是网格质量；由于现在大量使用分块网格，这时要看看两相邻块处的网格大小是不是相差较大，也就是看看有没有出现cell jump的情况，相邻网格的大小最好不要超过2倍的关系，这时出现高连续方程残差的一个主要原因，这需要在划分网格时做好规划6、网格独立性数值计算的与实验值之间的误差来源只要有这几个：物理模型近似误差（无粘或有粘，定常与非定常，二维或三维等等），差分方程的截断误差及求解区域的离散误差（这两种误差通常统称为离散误差），迭代误差（离散后的代数方程组的求解方法以及迭代次数所产生的误差），舍入误差（计算机只能用有限位存储计算物理量所产生的误差）等等。

在通常的计算中，离散误差随网格变细而减小，但由于网格变细时，离散点数增多，舍入误差也随之加大。

由此可见，网格数量并不是越多越好的。

再说说网格无关性的问题，由上面的介绍，我们知道网格数太密或者太疏都可能产生误差过大的计算结果，网格数在一定的范围内的结果才与实验值比较接近，这样在划分网格时就要求我们首先依据已有的经验大致划分一个网格进行计算，将计算结果（当然这个计算结果必须是收敛的）与实验值进行比较（如果没有实验值，则不需要比较，后面的比较与此类型相同），再酌情加密或减少网格，再进行计算，再与实验值进行比较，并与前一次计算结果比较，如果两次的计算结果相差较小（例如在2%），说明这一范围的网格的计算结果是可信的，说明计算结果是网格无关的。

再加密网格已经没有什么意义（除非你要求的计算精度较高）。

但是，如果你用粗网格也能得到相差很小的计算结果，从计算效率上讲，你就可以完全使用粗网格去完成你的计算。

加密或者减少网格数量，你可以以一倍的量级进行。

7、TDMA算法（c#描述）// <summary>/// 求解三对角矩阵/// </summary>/// <param name="a">系数矩阵对角下方的系数</param>/// <param name="b">系数矩阵对角线系数</param>/// <param name="c">系数矩阵对角上方系数</param>/// <param name="res">结果矩阵</param>/// <returns>返回求解后的结果，存放于数组中</returns>public static double [] SolveTDMA(double []a,double []b,double[]c,double []res){int dimension = b.Length; //结果阶数和b相同double []retMatrix = new double[dimension];double []beta = new double[dimension];if(b[0] != 0){beta[0] = c[0]/b[0];for(int i =1;i<dimension;i++){beta = c/(b - a*beta);}double []y = new double[dimension];y[0] = res[0]/b[0];for(int i = 1;i<dimension;i++){y = (res-a*y)/(b-a*beta);}retMatrix[dimension-1] = y[dimension-1];for(int i=dimension-2;i>=0;i--){retMatrix = y - beta* retMatrix;}}else{//出错，无法计算下去，因为b[0]是要做分母的}return retMatrix;}8、gambit网格质量判断具体内容可以参考Gambit Documentation中的Quality Type Definitions章节。

判断网格质量的方面有：Area单元面积，适用于2D单元，较为基本的单元质量特征。

Aspect Ratio长宽比，不同的网格单元有不同的计算方法，等于1是最好的单元，如正三角形，正四边形，正四面体，正六面体等；一般情况下不要超过5：1.Diagonal Ratio对角线之比，仅适用于四边形和六面体单元，默认是大于或等于1的，该值越高，说明单元越不规则，最好等于1，也就是正四边形或正六面体。

Edge Ratio长边与最短边长度之比，大于或等于1，最好等于1，解释同上。

EquiAngle Skew通过单元夹角计算的歪斜度，在0到1之间，0为质量最好，1为质量最差。

最好是要控制在0到0.4之间。

EquiSize Skew通过单元大小计算的歪斜度，在0到1之间，0为质量最好，1为质量最差。

2D质量好的单元该值最好在0.1以内，3D单元在0.4以内。

MidAngle Skew通过单元边中点连线夹角计算的歪斜度，仅适用于四边形和六面体单元，在0到1之间，0为质量最好，1为质量最差。