第5章通用后处理器(POST1)
静力分析
5.1概述
使用POST1通用后处理器可观察整个模型或模型的一部分在某一时间点（或频率）上针对指定载荷组合时的结果。POST1有许多功能，包括从简单的图象显示到针对更为复杂数据操作的列表，如载荷工况的组合。
要进入ANSYS通用后处理器，输入/POST1命令（Main Menu>General Postproc）.
5.2将数据结果读入数据库
POST1中第一步是将数据从结果文件读入数据库。要这样做，数据库中首先要有模型数据（节点，单元等）。若数据库中没有模型数据，输入RESUME命令(Utility Menu>File>Resume Jobname.db)读入数据文件Jobname.db。数据库包含的模型数据应该与计算模型相同，包括单元类型、节点、单元、单元实常数、材料特性和节点座标系。
注：数据库中被选来进行计算的节点和单元组应和模型中的节点和单元组属于相同组，否则会出现数据不匹配。有关数据不匹配的详细资料见5.2.2.3章。
一旦模型数据存在数据库中，输入SET,SUBSET或APPEND命令均可从结果文件中读入结果数据。
5.2.1 读入结果数据
输入SET命令（Main Menu>General PostProc>datatype）,可在一特定的载荷条件下将整个模型的结果数据从结果文件中读入数据库，覆盖掉数据库中以前存在的数据。边界条件信息（约束和集中力）也被读入，但这仅在存在单元节点载荷或反作用力的情况下，详情请见OUTRES命令。若它们不存在，则不列出或显示边界条件，但约束和集中载荷可被处理器读入，而且表面载荷和体积载荷并不更新，并保持它们最后指定的值。如果表面载荷和体积载荷是使用表格指定的，则它们将依据当前的处理结果集，表格中相应的数据被读入。加载条件靠载荷步和子步或靠时间（或频率）来识别。命令或路径方式指定的变元可以识别读入数据库的数据。例如：SET，2，5读入结果，表示载荷步为2，子步为5。同理，SET，，，，，3.89表示时间为3.89时的结果（或频率为3.89，取决于所进行分析的类型）。若指定了尚无结果的时刻，程序将使用线性插值计算出该时刻的结果。
结果文件（Jobname.RST）中缺省的最大子步数为1000，超出该界限时，需要输入SET，Lstep，LAST引入第1000个载荷步，使用/CONFIG增加界限。
注：对于非线性分析，在时间点间进行插值常常会降低精度。因此，要使解答可用，务必在可求时间值处进行后处理。
115
对于SET命令有一些便捷标号：
·SET，FIRST 读入第一子步，等价的GUI方式为First Set。
·SET，NEXT 读入第二子步，等价的GUI方式为Next Set。
·SET，LAST 读入最后一子步，等价的GUI方式为Last Set。
·SET命令中的NSET字段（等价的GUI方式为Set Number）可恢复对应于特定数据组号的数据，而不是载荷步号和子步号。当有载荷步和子步号相同的多组结果数据时，这对FLOTRAN的结果非常有用。因此，可用其特定的数据组号来恢复FLOTRAN的计算结果。SET 命令（或GUI中的List Results）LIST选项列出了其对应的载荷步和子步数，可在接下来的SET命令的NSET字段输入该数据组号，以申请处理正确的一组结果。
·SET命令中的ANGLE字段规定了谐调元的周边位置（结构分析－PLANE25，PLANE83和SHELL61；温度场分析－PLANE75和PLANE78）。
5.2.2 其他用于恢复数据的选项
其他GUI路径和命令也可恢复结果数据。
5.2.2.1定义待恢复的数据
POST1中的命令INRES(Main Menu>General PostProc>Data&File Opts)与PREP7和SOLUTION处理器中的OUTRES命令是姐妹命令，OUTRES命令控制写入数据库和结果文件的数据，而INRES命令定义要从结果文件中恢复的数据类型，通过命令SET，SUBSET和APPEND 等命令写入数据库。尽管不须对数据进行后处理，但INRES命令限制了恢复和写入数据库的数据量。因此，对数据进行后处理也许占用的时间更少。
5.2.2.2读入所选择的结果信息
为了只将所选模型部分的一组数据从结果文件读入数据库，可用SUBSET命令（Main Menu>General Postproc>By characteristic）。结果文件中未用INRES命令指定恢复的数据，将以零值列出。
SUBSET命令与SET命令大致相同，除了差别在于SUBSET只恢复所选模型部分的数据。用SUBSET命令可方便地看到模型的一部分的结果数据。例如，若只对表层的结果感兴趣，可以轻易地选择外部节点和单元，然后用SUBSET命令恢复所选部分的结果数据。
5.2.2.3 向数据库追加数据
每次使用SET，SUBSET命令或等价的GUI方式时，ANSYS就会在数据库中写入一组新数据并覆盖当前的数据。APPEND命令（Main Menu>General Postproc>By characteristic）从结果文件中读入数据组并将与数据库中已有的数据合并（这只针对所选的模型而言）。已有的数据库并不清零（或重写全部），而允许将被查询的结果数据并入数据库。
116
可用SET，SUBSET，APPEND命令中的任一命令从结果文件将数据读入数据库。命令方式之间或路径方式之间的唯一区别是所要恢复的数据的数量及类型。追加数据时，务必不要造成数据不匹配。例：请看下一组命令：
/ POST1
INRES,NSOL !节点DOF求解的标志数据
NSEL,S,NODE,,1,5 !选节点1至5
SUBSET,1 !从载荷步1开始将数据写入数据库
此时载荷步1内节点1到5的数据就存在于数据库中了。
NSEL,S,NODE,,6,10 ! 选节点6至10
APPEND,2 ! 将载荷步2的数据并入数据库中
NSEL,S,NODE,,1,10 ! 选节点1至10
PRNSOL,DOF ! 打印节点DOF求解结果
数据库当前就包含有载荷步1和载荷步2的数据。这样数据就不匹配。使用PRNSOL命令（Main Menu>General PostProc> List Results>Nodal Solution）时，程序将通知从第二个载荷步中取出数据，而实际上数据是从现存于数据库中的两个不同的载荷步中取得的。程序列出的载荷步仅为与此同时最近一次存入的载荷步相对应的载荷步。当然，若希望将不同载荷步的结果进行对比，将数据加入数据库中是很有用的。但若有目的地混合数据，要极其注意跟踪追加数据的来源。
在求解曾用不同单元组计算过的模型子集时，为避免出现数据不匹配，按下列任一方法进行。
·不要重选在后处理的当前解答中未被选中的任何单元
·从ANSYS数据库中删除以前的解答。可在多步的求解过程中从每步的求解后退出ANSYS或在求解中间存储数据库。
详见ANSYS命令参考中对NSEL，APPEND，PRNSOL，SUBSET命令的描述。
若想清空数据库中所有以前的数据，使用下列任一方式：
命令：LCZERO
GUI：Main Menu>General PostProc>Load Case>Zero Load Case
上述两种方法均会将数据库中所有以前的数据置零。因而可重新进行数据存储。若在向数据库追加数据之前将数据库置零，假如SUBSET和APPEND命令中的变元等价，其结果与使用SUBSET命令（或等价的GUI路径）一样。
注：SET命令可用的全部选项对SUBST命令和APPEND命令完全可用。
缺省情况下，SET，SUBSET和APPEND命令将寻找这些文件中的一个：Jobname.RST, Jobname.RTH,Jobname.RMG,Jobname.RFL。在使用SET，SUBSET和APPEND命令之前用FILE 命令可指定其它文件名（Main Menu>General Postproc>Data&File Opts）。
117
5.2.3 创建单元表
ANSYS程序中单元表有两个功能：第一，它是在结果数据中进行数学运算的工具。第二，它能够访问其他方法无法直接访问的单元结果。例如：从结构一维元派生的数据（尽管SET，SUBSET和APPEND命令将所有申请的结果项读入数据库中，但并非所有的数据均可直接用PRNSOL命令和PLESON等命令访问）。
将单元表作为扩展表，每行代表一单元，每列则代表单元的特定数据项。例如：一列可能包含单元的平均应力SX，而另一列则代表单元的体积，第三列则包含各单元质心的Y座标。
使用下列任一命令创建或删除单元表：
命令：ETABLE
GUI：Main Menu>General Postproc>Element Table>Define Table or Erase Table
5.2.3.1填上按名字来识别变量的单元表
为识别单元表的每列，在GUI方式下使用Lab字段或在ETABLE命令中使用Lab变元给每列分配一个标识，该标识将作为所有的以后的包括该变量的POST1命令的识别器。进入列中的数据靠Item名和Comp名以及ETABLE命令中的其它两个变元来识别。例如：对上面提及的SX应力，SX是标识，S将是Item变元，X将是Comp变元。
有些项，如单元的体积，不需Comp。这种情况下，Item为VOLU, 而Comp为空白。按Item和Comp（必要时）识别数据项的方法称为填写单元表的“元件名”法。使用“元件名”法访问的数据通常是那些针对大多数单元类型和单元类型组而经过计算的数据。
ETABLE命令的文档通常列出了所有的Item和Comp的组合情况。要清楚何种组合有效，见ANSYS单元参考手册中每种单元描述中的“单元输出定义”。表5—1是一个对于BEAM4的列表示例，可在表中的NAME列中的冒号后面使用任意名称作为名字，通过“Compmnent Name”法将单元表填满。冒号前面的名字部分应输入作为ETABLE命令的Item变元，冒号后的部分（如果有的话）应输入作为ETABLE命令的Comp变元，O列与R列表示在Jobname.OUT 文件（O）中或结果文件（R）中该项是否可用：Y表示该项总可用，数字则表示表的脚注，描述何时该项在一定条件下可用，而a则显示该项不可用。
表5—1三维BEAM4单元输出定义
118
1．项目经I端、中间位置(见KEYOPT(9))及J端重复进行。
2．若KEYOPT(6)=1
5.2.3.2填充按序号识别变量的单元表
可对每个单元加上不平均的或非单值载荷，将其填入单元表中。该数据类型包括积分点的数据、从结构一维单元（如杆，梁，管单元等）和接触单元派生的数据、从一维温度单元派生的数据、从层状单元中派生的数据等。这些数据将列在“Item and Sequence Numbers for the ETABLE and ESOL Commands”表中，而每一单元类型在ANSYS单元参考指南中都有描述。表5—2是BEAM4的示例。
表中的数据分成项目组（如LS,LEPEL,SMISC等），项目组中每一项有一用于识别的序号。将项目组（如LS,LEPEL,SMISC等）作为ETABLE命令的Item变元，将序列号作为Comp 变元，将数据装入单元列表中，称之为填写单元表的“序列号”法。例如：BEAM4单元的J 点处的最大应力为Item=NMISC及Comp=3。而单元（ E）的初始轴向应变（EPINAXL）为Item=LEPYH，Comp=11 。
表5-2梁单元（KEYOPT（9）=0）的ETABLE和ESOL命令中的项目及序号
119
120
对于某些一维单元，如BEAM4，KEYOPT设置控制了计算数据的量。这能改变特定数据项的序号。因此，在这些情况下提供了有关每个KEYOPT设置的表格。表5-3如表5-2一样显示了关于BEAM4的相同信息，但列出的为KEYOPT（9）=3时的序号（3个中间计算点）。例如：当KEYOPT（9）=0时，单元J端Y向的力矩（MMOM）在表5-2中是序号11（SMISC项），而当KEYOPT（9）=3时，其序号（表5-3）为29。
表5-3ETABLE命令和ESOL命令的BEAM4(KEYOPT(9)=3) 的项目名和序号
121
5.2.3.3 定义单元表的注释
·ETABLE命令仅对选中的单元起作用，即只将所选单元的数据送入单元表中。在ETABLE 命令中改变所选单元，可以有选择地填写单元表的行。
·相同序号的组合表示对不同单元类型有不同数据。例如：组合SMISC，1对梁单元表示MFOR(X)（单元X向的力），对SOLID45单元表示P1（面1上的压力），对CONTACT48单元表示FNTOT（总的法向力）。因此，若模型中有几种单元类型的组合，务必要在使用ETABLE 命令前选择一种类型的单元（用ESEL命令或Utility Menu >Select>Entities）
·ANSYS程序在读入不同组的结果（例如对不同的载荷步）或在修改数据库中的结果（例如在组合载荷工况），不能自动刷新单元表。例如：假定模型由提供的样本单元组成，在POST1中发出下列命令：
SET,1 ! 读入载荷步1结果
ETABLE,ABC,1S,6！在以ABC开头的列下将J端（KEYOPT（9）=0的SDIR移入单元表中
SET,2 ！读入载荷步2中结果
此时，单元表“ABC”列下仍含有载荷步1的数据。用载荷步2中的数据更新该列数据时，应用命令ETABLE，REFL或通过GUI方式指定更新项。
·可将单元表当作一“工作表”，对结果数据进行计算。详见5.4章对该特征的描述。
·使用POST1中的SAVE,FNAME,EXT或/EXIT,ALL在退出ANSYS程序时，可以对单元表进行存盘（若使用GUI方式，选择UtilityMenu>File>Save as或Utility Menu>File>Exit 后按照对话框内的提示进行）。这样可将单元表及其余数据存到数据库文件中。
·为从内存中删除整个单元表，用ETABLE，ERASE(Main Menu>General Postproc >Element Table>Erase Table),或用ETABLE，LAB, ERASE删去单元表中的Lab 列。用RESET命令(Main Menu>General Postproc>Reset)可自动删除内存中的单元表。
5.2.4 对主应力的专门研究
在POST1中，SHELL61单元的主应力不能直接得到，缺省情况下，可得到其它单元的主应力，除以下两种情况之一：
·在SET命令中要求进行时间插值或定义了某一角度
122
·执行了载荷工况操作时
在上述情况（对SHELL61单元包括所有情况）下，必须用Main Menu>General Postproc>Load Case>Line Elem Stress或用LCOPER,LPRIN命令以计算主应力。然后通过ETABLE命令或用其他适当的打印或绘图命令访问该数据。
5.2.5 读入FLOTRAN的计算结果
使用命令FLREAD(Main Menu>General Postproc>FLOTRAN 2.1A)可以将结果从FLOTRAN 的“剩余”文件读入数据库。FLOTRAN的计算结果(Jobname.RFL)可以用普通的后处理函数或命令（例如SET命令）（Utility Menu>List>Results>Load Step Summary）读入。
5.2.6 数据库复位
RESET命令（Main Menu>General Postproc> Reset）可在不脱离POST1情况下初始化POST1命令的数据库缺省部分，该命令在离开或重新进入ANSYS程序时的效果相同。
5.3 在POST1中观察结果
一旦所需结果存入数据库，可通过图象显示和表格方式观察。另外，可映射沿某一路径的结果数据（详见5.3.4节）。
5.3.1图象显示结果
图象显示可能是观察结果的最有效方法。POST1可显示下列类型图象：
·云图显示
·变形后的形状显示
·矢量图显示
·路径绘图
·反作用力显示
·粒子流轨迹
5.3.1.1云图显示
云图显示表现了结果项（如应力、温度、磁场磁通密度等）在模型上的变化。云图显示中有四个可用命令：
命令：PLNSOL
GUI：Main Menu>General Postproc> Plot Results>Nodal Solu
命令：PLESOL
GUI：Main Menu>General Postproc>Plot Results>Element Solu
123
命令：PLETAB
GUI：Main Menu>General Postproc> Plot Results>Elem Table
命令：PLLS
GUI：Main Menu>General Postproc> Plot Results>Line Elem Res
PLNSOL命令生成连续的过整个模型的云图。该命令或GUI方式可用于原始解或派生解。对典型的单元间不连续的派生解，在节点处进行平均，以便可显示连续的云图。下面将举出原始解（TEMP）和派生解（TGX）梯度显示的示例。
PLNSOL, TEMP!原始解：自由度TEMP
图5-1：使用PLNSOL得到的原始解的云图
若有PowerGraphics(性能优化的增强型RISC体系图形)，可用下列任一命令中控制派生数据的平均。
命令：AVRES
GUI：Main Menu>General Postproc> Options for Outp
Utility Menu>List>Results>Options
上述任一命令均可确定在材料及（或）实常数不连续的单元边界上是否对结果进行平均。详见第10章“性能优化的增强型RISC体系图形”。
当心：若PowerGraphics无效（对大多数单元类型而言，这是默认值），不能用AVRES 命令去控制平均计算。平均算法则不管连接单元的节点属性如何，均会在所选单元上的所有节点处进行平均操作。这对材料和几何形状不连续处是不合适的。当对派生数据进行云图显示时（这些数据在节点处已做过平均），务必选择相同材料、相同厚度、（对板单元），相同坐标系等下的单元。
PLNSOL,TG,X！派生数据：温度梯度函数TGX
124
详见PLNSOL命令的描述。
图5-2：用PLNSOL命令对派生数据进行云图显示
PLESOL命令在单元边界上生成不连续的云图，该显示主要用于派生的解数据。如：X
PLESOL，TG，
图5-4：平均的PLETAB云图
PLETAB命令云图显示单元表中的数据。在PLETAB命令中的AVGLAB字段，提供了是否对节点处数据进行平均的选择项（对连续云图，平均；缺省状态，对不连续云图，不平均）。下例假设采用SHELL99单元（层状壳）模型，对结果进行平均和不平均之间的差异如图所示。
125
ETABLE,SHEARXZ,SMISC,9 !在第二层底部存在层内剪切(ILSXZ)
PLETAB,SHEARXZ,AVG ! SHEARXZ的平均云图。
PLETAB，SHEARXZ,NOAVG ! SHEARXZ的未平均(默认值)的云图。
图5-5：未平均的PLETAB云图。
PLLS命令用云图的形式显示一维单元的结果。该命令也要求数据存储在单元表中。该形式的显示常用于梁分析中对剪力和力矩图进行表示。下例假定有KEYOPT(9)=1梁模型BEAM3(二维的梁):
ETABLE,IMOMENT,SMISC,6 !I端的弯矩，命名为IMOMENT
ETABLE,JMOMENT,SMISC,18 !J端的弯矩命名为JMOMENT
PLLS,IMOMENT,JMOMENT !显示IMOMENT,JMOMENT结果。
图5-6用PLLS命令显示的弯矩图
PLLS命令只会在单元的I,J节点处的数值间画一直线，不管结果沿单元长度如何变化。可用负的比例因子将图形倒过来。
126
注释：
·可用/CTYPE命令（Utility Menu>PlotCtrls>Style>Contours>Contour Style）首先设置KEY为1来生成等轴侧的云图显示。见第6章有关等轴侧的信息。
·平均主应力：缺省情况下，各节点处的主应力根据平均分应力计算。也可反过来做，首先计算每个单元的主应力，然后在各节点处平均。步骤如下：
命令：AVPRIN
GUI：Main Menu>General Postproc> Options for Outp
Utility Menu>List >Results>Options
该法不常用，但在特定情况下很有用。在不同材料的结合面处不应采用平均算法。
·矢量求和数据：与主应力的做法相同。缺省情况下，在每个节点处的矢量和的模（平方和的开方）是按平均过的分量来求的。用AVPRIN命令，可反过来计算，先计算每单元矢量和的模，然后在节点处进行平均。
·壳单元或分层壳单元：缺省情况下，壳单元或分层壳单元得到的计算结果假定为其上表面处的结果。要显示上表面、中部或下表面的结果，用SHELL命令（Main Menu>General Postproc> Options for Outp）。对于分层单元，使用LAYER命令(Main Menu>General Postproc> Options for Outp)指明需显示的层号。
·Von Mises当量应力（EQV）:使用命令AVPRIN可以改变用来计算这些量的有效泊松比。
命令：AVPRIN
GUI：Main Menu>General Postproc>Plot Results>-Contour Plot-Nodal Solu
Main Menu>General Postproc>Plot Results>-Contour Plot-Element Solu
Utility Menu>Plot>Results >Contour Plot>Elem Solution
典型情况下，对弹性当量应变（EPEL,EQV），可将有效泊松比设为输入泊松比，对非弹性应变（EPPL,EQV或EPCR,EQV），设为0.5。对于整个当量应变（EPTOT,EQV），应在输入泊松比和0.5之间选用一有效泊松比。另一种方法是，用命令ETABLE存储当量弹性应变，使有效泊松比等于输入泊松比；在另一张表中用0.5作为有效泊松比存储当当量塑性应变，然后用SADD命令将两张表合并，得到整个当量应变。
·/EFACET命令的效果：当你看连续的云图时(PLNSOL),通过不同的/EFACET命令设置,你将看到不同的线图。当设置/EFACET,1,时,区域中间位置的梯度值由平均过后的相邻角点的梯度值插值得到;当/EFACET,2,时,先计算将每个单元中点节点的梯度值,由未平均过的角点处的梯度值插值得到,然后所有中点节点的梯度值由(PLNSOL)平均再描绘出线图;当
127
/EFACET,4,时,ANSYS利用形状函数(高阶P单元除外)在每个单元的边界处计算三个网格交叉节点的梯度值,交叉节点值首先在每个单元内计算然后再由(PLNSOL)一起平均。因此，区域中间位置的梯度值将在不同的/EFACET设置下不同。当你选用一个不是全局坐标系时(KCN=1,2等)，梯度值首先在全局坐标系中平均，然后将该平均值转化为你指定的结果坐标系值。
大多数情况下，不管/EFACET设置如何，PLNSOL命令将得到相同的云图。但是，若在连接处存在着除KCN=0的其他RSYS设置时，改变/EFACET的设置时使用PLNSOL你将看到不同的云图。
5.3.1.2 变形后的形状显示
在结构分析中可用这些显示命令观察结构在施加的载荷下的变形情况。用下列方法之一可产生变形形状显示：
命令：PLDISP
GUI：Utility Menu>Plot>Results>Deformed Shape
Main Menu>General Postproc>Plot Results> Deformed Shape
例如：可用以下PLDISP命令：
图5-7 PLDISP图样
PLDISP，1 ！变形后的形状与原始形状叠加在一起
可用命令/DSCALE来改变位移比例因子。
在你进入POST1时，当心所有载荷符号被自动关闭。以后再次进入PREP7或SOLUTION 处理器时仍不会见到这些载荷符号。若在POST1中打开所有载荷符号，结果会在变形形状图上显示载荷。
128
5.3.1.3 矢量显示
矢量显示用箭头显示模型中某个矢量大小和方向的变化。平移（U）、转动（ROT）、磁力矢量势（A）、磁通密度（B）、热通量（TF）、温度梯度（TG）、液流速度（V）、主应力（S）等都是矢量的例子。
用下列方法可产生矢量显示：
命令：PLVECT
GUI： Main Menu>General Postproc>Plot Results>Predefined or User-Defined
可用下列方法改变矢量箭头长度比例：
命令：/VSCALE
GUI：Utility Menu>PlotCtrls>Style>Vector Arrow Scaling
PLVECT,B ! 磁通密度（B）的矢量显示
图5-8磁场强度的PLVECT矢量图
在PLVECT命令中定义两个或两个以上分量，你可生成自己所需的矢量值。
5.3.1.4 路径图
路径图是显示某个量沿过模型的某一预定路径的变化图。要产生路径图，执行下述步骤：
用命令PATH定义路径属性(Main Menu>General Postproc>Path Operations> Define Path>Defined Paths)。
用命令PPATH定义路径点(Main Menu>General Postproc>Path Operations> Define Path>Modify Path)。
用命令PDEF将所需的量映射到路径上(Main Menu>General Postproc>Map Onto Path)。
129
用命令PLPATH和PLPAGM显示结果(Main Menu>General Postproc>Path Operations>Plot Path Items)。
有关路径图的进一步内容见5.3.4节。
5.3.1.5 反作用力显示
与边界条件显示相同，用命令/PBC下的RFOR或RMOM来激活。以后的任何显示（由NPLOT,EPLOT或PLDISP命令生成）将在定义了DOF约束的点处显示反作用力符号。约束方程中某一自由度节点力之和不应包含过该方程的力。详见15.8.2节中的ANSYS理论参考。
如反作用力一样，也可用命令/PBC(Utility Menu>Plotctrls>Symbols)中的NFOR或NMOM项显示节点力，这是单元在其节点上施加的力。每一节点处这些力之和通常为0，约束点处或加载点除外。
缺省情况下，打印出的或显示出的力（或力矩的）的数值代表合力（静力、阻尼力和惯性力的总和）。FORCE命令(Main Menu>General Postproc>Options For Outp)可将合力分解成各分力。
5.3.1.6 粒子流和带电粒子轨迹
粒子流轨迹是一种特殊的图象显示形式，用于描述流动流体中粒子的运动情况。带电粒子轨迹是显示带电子粒子在电、磁场中如何运动的图象显示。详见12章图形显示及15章粒子轨迹动画的信息。见ANSYS理论参考的第15章对电磁粒子轨迹的简化假设。
粒子流或带电粒子轨迹需要两个函数：
TRPOIN命令(Main Menu>General Postproc>Plot Results>Defi Trace Pt)。任一方式都在路径轨迹上定义一点（起点、终点或两点中的任意一点）。
PLTRAC命令(Main Menu>General Postproc>Plot Results>Plot Flow Tra)。任一方式都在单元上显示上产生流动轨迹，能同时定义和显示多达50点。
PLTRAC图样如下所示：
图5-9 粒子流轨迹示例
PLTRAC命令中的Item字段和comp字段能使用户看到某一特定项的变化情况（如：对粒子流动轨迹为速度、压力和温度或对带电粒子轨迹为电荷）。项目的变化情况沿路径用按梯度变化的色带显示出来，其他命令有：
130
·TRPLIS命令（Main Menu>General Postproc>Plot Results>List Trace Pt）：列出轨迹点。
·TRPDEL命令（Main Menu>General Postproc>Plot Results>Dele Trace Pt）：删除轨迹点。
·TRTIME命令（Main Menu>General Postproc>Plot Results>Time Interval）：定义流动轨迹时间间隔。
·ANFLOW命令（Main Menu>General Postproc>Plot Results>Paticle Flow）：生成粒子流的动画序列。
注意：
·在描绘粒子流轨迹时会生成三个数组参数：TRACPOIN, TRACDATA TRACLABL.这些参数可以用来得到以粒子的速度和已经过去的时间为变量的粒子轨迹方程。利用这些参数来得到一个命名为PATHNAME的轨迹的示例如下：
*get,npts,PARM,TRACPOIN,DIM,x
PATH,PATHNAME,npts,9,1
PAPUT,TRACPOIN,POINTS
PAPUT,TRACDATA,TABLES
PAPUT,TRACLABL,LABELS
PRPATH,S,T_TRACE,VX_TRACE,VY_TRACE,VZ_TRACE,VS_TRACE
·粒子流轨迹偶尔会无明显原因地停止。在靠近管壁处的静止流体区域，或当粒子沿单元边界运动时，会出现这种情况。为解决这个问题，在流线交叉方向轻微调整子粒初始点。
·对带电粒子轨迹，用TRPON命令（Main Menu>General Postproc>Plot Results>Defi Trace Pt）输入的变量Chrg和Mass在MKS单位制中具有相应的单位“库仑”和“千克”。
·粒子轨迹跟踪算法会导致死循环，例如某一带电粒子轨迹会导致无限循环。要避免出现死循环，PLTRAC命令的MXLOOP变元设置了极限值。
5.3.1.7 破碎图
131
图5-10具有裂缝的混凝土梁
若在模型中有SOLID65单元，你可用PLCRACK命令(Main Menu>General Postproc>Plot Results>Crack/Crush)确定那些单元已断裂或碎开。以小圆圈标出已断裂，以小八边形表示混凝土已压碎（见图5-10）。在使用不隐藏矢量显示的模式下，可见断裂和压碎的符号。为指定这一设备，用命令/DEVICE,VECTOR,ON(Utility Menu>Plotctrls>Device Options)。
5.3.2合成表面结果
INTSRF命令(Main Menu>General Postproc>Nodal Calcs>Pressure Integrl)允许合成外表面处的节点结果，首先必须要在外表面选择所需要合成的节点。
利用INTSRF计算升力和阻力时，当表面为流-固表面时，仅需选择流体表面节点来进行合成。利用NSEL命令(Utility Menu>Select>Entities)和EXT选项选择节点。
利用INTSRF计算升力和阻力时，必须指定一个结果坐标系，其X方向为流体区域的引入方向，其Y方向为重力的方向。以便使的升力方向为X方向，阻力方向为Y方向。然后利用INTSRF，PRES和INTSRF，TAUW分别得到升力和阻力，也可以利用INTSRF，FLOW得到分开的升力和阻力。结果将写入输出文件中(jobname.OUT)。
合成结果为结果坐标系下的数值(见RSYS命令)，结果坐标系的类型必须和分析中使用的坐标系类型相匹配，但可以对力和力矩进行所需的转换和旋转。利用*GET命令得到结果。
********PAUSEPAUSEPAUSEPAUSEPAUSEPAUSEPAUSE*******
5.3.3 用表格形式列出结果
将结果存档的有效方法（例如：报告、呈文等）是在POST1中制表。列表选项对节点、单元、反作用力等求解数据可用。
132
PRESOL,ELEM的样表
PRINT ELEM ELEMENT SOLUTION PER ELEMENT
***** POST1 ELEMENT SOLUTION LISTING *****
LOAD STEP1SUBSTEP=1
TIME=1.0000LOAD CASE=0
EL=1NODES=13MAT=1
BEAM3
TEMP =0.000.000.000.00
LOCATIONSDIRSBYTSBYB
1 (I)0.00000E+00130.00-130.00
2 (J)0.00000E+00104.00-104.00
LOCATIONSMAXSMIN
1 (I)130.00-130.00
2 (J)104.00-104.00
LOCATIONEPELDIREPELBYTEPELBYB
1 (I)0.0000000.000004-0.000004
2 (J)0.0000000.000003-0.000003
LOCATIONEPTHDIREPTHBYTEPTHBYB
1 (I)0.0000000.0000000.000000
2 (J)0.0000000.0000000.000000
EPINAXL =0.000000
EL=2NODES=34MAT=1
BEAM3
TEMP =0.000.000.000.00
LOCATIONSDIRSBYTSBYB
1 (I)0.00000E+00104.00-104.00
2 (J)0.00000E+0078.000-78.000
LOCATIONSMAXSMIN
1 (I)104.00-104.00
2 (J)78.000-78.000
LOCATIONEPELDIREPELBYTEPELBYB
1 (I)0.0000000.000003-0.000003
2 (J)0.0000000.000003-0.000003
LOCATIONEPTHDIREPTHBYTEPTHBYB
1 (I)0.0000000.0000000.000000
2 (J)0.0000000.0000000.000000
EPINAXL =0.000000
5.3.3.1 列出节点、单元求解数据
用下列方式可以列出指定的节点求解数据（原始解及派生解）：
命令：PRNSOL
GUI：Main Menu>General Postproc>List Results>Nodal Solution
133
用下列方式可以列出所选单元的指定结果：
命令：PRNSEL
GUI：Main Menu>General Postproc>List Results>Element Solution
要获得一维单元的求解输出，在PRNSOL命令中指定ELEM选项，程序将列出所选单元的所有可行的单元结果。
PRNSOL，S的样表
PRINT SNODAL SOLUTION PER NODE
***** POST1 NODAL STRESS LISTING *****
LOAD STEP=5SUBSTEP=2
TIME=1.0000LOAD CASE=0
THE FOLLOWING X,Y,Z VALUES ARE IN GLOBAL COORDINATES
NODESXSYSZSXYSYZSXZ
1148.01-294.54.00000E+00 -56.256.00000E+00.00000E+00
2144.89-294.83.00000E+0056.841.00000E+00.00000E+00
3241.8473.743.00000E+00 -46.365.00000E+00.00000E+00
4401.98-18.212.00000E+00 -34.299.00000E+00.00000E+00
5468.15-27.171.00000E+00.48669E-01.00000E+00.00000E+00
6401.46-18.183.00000E+0034.393.00000E+00.00000E+00
7239.9073.614.00000E+0046.704.00000E+00.00000E+00
8-84.741-39.533.00000E+0039.089.00000E+00.00000E+00
93.2868-227.26.00000E+0068.563.00000E+00.00000E+00
10-33.232-99.614.00000E+0059.686.00000E+00.00000E+00
11-520.81-251.12.00000E+00.65232E-01 .00000E+00.00000E+00
12-160.58-11.236.00000E+0040.463.00000E+00.00000E+00
13-378.5555.443.00000E+0057.741.00000E+00.00000E+00
14-85.022-39.635.00000E+00 -39.143.00000E+00.00000E+00
15-378.8755.460.00000E+00 -57.637.00000E+00.00000E+00
16-160.91-11.141.00000E+00 -40.452.00000E+00.00000E+00
17-33.188-99.790.00000E+00 -59.722.00000E+00.00000E+00
183.1090-227.24.00000E+00 -68.279.00000E+00.00000E+00
1941.81151.777.00000E+00 -66.760.00000E+00.00000E+00
20-81.0049.3348.00000E+00 -63.803.00000E+00.00000E+00
21117.64-5.8500.00000E+00 -56.351.00000E+00.00000E+00
22-128.2130.986.00000E+00 -68.019.00000E+00.00000E+00
23154.69-73.136.00000E+00.71142E-01.00000E+00.00000E+00
24-127.64-185.11.00000E+00.79422E-01.00000E+00.00000E+00
25117.22-5.7904.00000E+00 56.517.00000E+00.00000E+00
26-128.2031.023.00000E+0068.191.00000E+00.00000E+00
134
压电变换器的自振频率分析及详细过程 1.模态分析的定义及其应用 模态分析用于确定设计结构或机器部件的振动特性(固有频率和振型)，即结构的固有频率和振型，它们是承受动态载荷结构设计中的重要参数。同时，也可以作为其它动力学分析问题的起点，例如瞬态动力学分析、谐响应分析和谱分析，其中模态分析也是进行谱分析或模态叠加法谐响应分析或瞬态动力学分析所必需的前期分析过程。 ANSYS的模态分析可以对有预应力的结构进行模态分析和循环对称结构模态分析。前者有旋转的涡轮叶片等的模态分析，后者则允许在建立一部分循环对称结构的模型来完成对整个结构的模态分析。 ANSYS提供的模态提取方法有:子空间法(subspace)、分块法(block lancets),缩减法(reduced/householder)、动态提取法(power dynamics)、非对称法(unsymmetric),阻尼法(damped), QR阻尼法(QR damped)等，大多数分析都可使用子空间法、分块法、缩减法。 ANSYS的模态分析是线形分析，任何非线性特性，例如塑性、接触单元等，即使被定义了也将被忽略。 2.模态分析操作过程 一个典型的模态分析过程主要包括建模、模态求解、扩展模态以及观察结果四个步骤。 (1).建模 模态分析的建模过程与其他分析类型的建模过程是类似的，主要包括定义单元类型、单元实常数、材料性质、建立几何模型以及划分有限元网格等基本步骤。 (2).施加载荷和求解 包括指定分析类型、指定分析选项、施加约束、设置载荷选项，并进行固有频率的求解等。 指定分析类型，Main Menu- Solution-Analysis Type-New Analysis,选择Modal。 指定分析选项，Main Menu-Solution-Analysis Type-Analysis Options，选择MODOPT(模态提取方法〕，设置模态提取数量MXPAND. 定义主自由度，仅缩减法使用。 施加约束，Main Menu-Solution-Define Loads-Apply-Structural-Displacement。 求解，Main Menu-Solution-Solve-Current LS。 (3).扩展模态 如果要在POSTI中观察结果，必须先扩展模态，即将振型写入结果文件。过程包括重新进入求解器、激话扩展处理及其选项、指定载荷步选项、扩展处理等。 激活扩展处理及其选项，Main Menu-Solution-Load Step Opts-Expansionpass-Single Expand-Expand modes。 指定载荷步选项。 扩展处理，Main Menu-solution-Solve-Current LS。 注意:扩展模态可以如前述办法单独进行，也可以在施加载荷和求解阶段同时进行。本例即采用了后面的方法 (4).查看结果 模态分析的结果包括结构的频率、振型、相对应力和力等
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 /TITILE,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,,ALL SFL,3,PRES,1000 SFL,2,PRES,1000 ...... SOLVE !求解标识 FINISH !求解模块结束标识 /POST1 !进入通用后处理器标识 ...... /POST26 !进入时间历程后处理器 …… /EXIT,SAVE !退出并存盘 以下是日志文件中常出现的一些命令的标识说明,希望能给大家在整理LOG文件时有所帮助 /ANGLE !指定绕轴旋转视图 /DIST !说明对视图进行缩放 /DEVICE !设置图例的显示,如:风格,字体等 /REPLOT !重新显示当前图例 /RESET !恢复缺省的图形设置 /VIEW !设置观察方向 /ZOOM !对图形显示窗口的某一区域进行缩放
514. TALLOW，TEMP1，TEMP2，TEMP3，TEMP4，TEMP5，TEMP6（定义允许应力值相应的温度） 515. TB，Lab，MAT，NTEMP，NPTS，TBOPT，EOSOPT（在非线性材料属性或特殊单元输入中激活一单元表格） 516. TBDATA，STLOC，C1，C2，C3，C4，C5，C6（定义单元表格中的数据） 517. TBLIST，Lab，MAT（列表显示材料非线性特性） 518. TBPLOT，Lab，MAT，TBOPT，TEMP，SEGN（图形显示非线性材料的应力-应变曲线）519. TBPT, oper, x,y（在应力-应变曲线上定义一个点） 【注】oper: defi 定义一个点 dele 删除一个点 x,y：坐标 520. TCHG，ELEM1，ELEM2，ETYPE2（将四面体退化单元转化为非退化单元） 521. TIME，TIME（通过时间定义载荷步） 522. TIMP，ELEM，CHGBND，IMPLEVEL（对不附属于体的四面体单元进行改进） 523. /TLABEL，XLOC，YLOC，Text（使用文字注释） 524. TOFFST，VALUE（选择温度的单位） 525. TORQ2D（计算磁场中物体上的扭矩） 526. TORQC2D，RAD，NUMN，LCSYS（计算磁场中物体上环行路径的扭矩） 527. TORQSUM，Cnam1，Cnam2，…，Cnam8，Cnam9（对2-D平面问题中单元上的电磁麦克斯韦和虚功扭矩求和） 528. TORUS，RAD1，RAD2，RAD3，THETA1，THETA2（生成环体） 【注】RAD1，RAD2，RAD3中最大直径为主半径，最小为内半径，中间值为外半径。529. TRANSFER，KCNTO，INC，NODE1，NODE2，NINC（将节点模式转换到另一坐标系中）530. TREF，TREF（定义参考温度） 531. /TRIAD，Lab（控制是否显示整体坐标系标志，并对其位置进行定义） 【注】Lab=ORIG（在原点显示坐标系）、OFF（关闭显示）、LBOT（在左下角显示坐标系）、RBOT（在右下角显示坐标系）、LTOP（在左上角显示坐标系）、RTOP（在右上角显示坐标系）。532. /TRLCY，Lab，TLEVEL，N1，N2，NINC（透明显示） 533. TRPDEL，NTRP1，NTRP2，TRPINC（删除轨迹点） 534. TRPLIS，NTRP1，NTRP2，TRPINC（列表显示轨迹点信息） 535. TRPOIN，X，Y，Z，VX，VY，VZ，CHRG，MASS（定义粒子流轨迹上的点） 536. TRTIME，TIME，SPACING，OFFSET，SIZE，LENGTH（定义流动轨迹时间间隔） 537. /TSPEC，TCOLOR，TSIZE，TXTHIC，PANGLE，IANGLE（定义文字标注属性） 538. TUNIF，TEMP（定义结构中所有节点的温度）。 【注】适用于均匀温度负载时使用） 539. /TXTRE，Lab，NUM，N1，N2，NINC（为所选项选择纹理） /TXTRE，VOLU，NUM，N1，N2，NINC（为体选择纹理） /TXTRE，ON（激活纹理显示） 540. /TYPE，WN，Type（定义显示类型） 541. TYPE，ITYPE（指定单元类型） 542. /UDOC，Wind，Class，Key（指定图例栏中图例和文本在窗口中的位置） 543. UIMP，MAT，Lab1，Lab2，Lab3，VAL1，VAL2，VAL3（求解过程中修改材料特性）544. /UNITS，Label，LENFACT，MASSFACT，TIMEFACT，TEMPFACT，TOFFSET，CHARGEFACT，FORCEFACT，HEATFACT（选择单位制）
（转）ANSYS学习也有一个来月的时间了，可是还是什么都不会！郁闷！整理了一些ANSYS 常用的命令；但深知自己的水平，还不敢保证完全正确；给大家一些参考，望指正： 1. A，P1，P2，…，P17，P18（以点定义面） 2. AADD，NA1，NA2，…NA8，NA9（面相加） 3. AATT，MAT，REAL，TYPE，ESYS，SECN（指定面的单元属性） 【注】ESYS为坐标系统号、SECN为截面类型号。 4. *ABBR，Abbr，String（定义一个缩略词） 5. ABBRES，Lab，Fname，Ext（从文件中读取缩略词） 6. ABBSAVE，Lab，Fname，Ext（将当前定义的缩略词写入文件） 7. ABS，IR，IA，--，--，Name，--，--，FACTA（取绝对值） 【注】************* 8. ACCAT，NA1，NA2（连接面） 9. ACEL，ACEX，ACEY，ACEZ（定义结构的线性加速度） 10. ACLEAR，NA1，NA2，NINC（清除面单元网格） 11. ADAMS，NMODES，KSTRESS，KSHELL 【注】************* 12. ADAPT, NSOLN, STARGT, TTARGT, FACMN, FACMX, KYKPS, KYMAC 【注】************* 13. ADD，IR, IA, IB, IC, Name, --,-- , FACTA, FACTB, FACTC（变量加运算） 14. ADELE，NA1，NA2，NINC，KSWP（删除面） 【注】KSWP =0删除面但保留面上关键点、1删除面及面上关键点。 15. ADRAG，NL1，NL2，…，NL6，NLP1，NLP2，…，NLP6（将既有线沿一定路径拖拉成面） 16. AESIZE，ANUM，SIZE（指定面上划分单元大小） 17. AFILLT，NA1，NA1，RAD（两面之间生成倒角面） 18. AFSURF，SAREA，TLINE（在既有面单元上生成重叠的表面单元） 19. *AFUN, Lab（指定参数表达式中角度单位） 20. AGEN, ITIME, NA1, NA2, NINC, DX, DY, DZ, KINC, NOELEM, IMOVE（复制面） 21. AGLUE，NA1，NA2，…，NA8，NA9（面间相互粘接） 22. AINA，NA1，NA2，…，NA8，NA9（被选面的交集） 23. AINP，NA1，NA2，…，NA8，NA9（面集两两相交） 24. AINV，NA，NV（面体相交） 25. AL，L1，L2，…，L9，L10（以线定义面） 26. ALIST，NA1，NA2，NINC，Lab（列表显示面的信息） 【注】Lab=HPT时，显示面上硬点信息，默认为空。 27. ALLSEL，LabT，Entity（选择所有实体） 【注】LabT=ALL（指定实体及其所有下层实体）、BELOW（指定实体及其下一层实体）；Entity=ALL、VOLU、AREA、LINE、KP、ELEM、NODE。 28. AMESH，NA1，NA2，NINC（划分面生成面单元） AMESH，AREA，KP1，KP2，KP3，KP4（通过点划分面单元） 29. /AN3D，Kywrd，KEY（三维注释） 30. ANCNTR，NFRAM，DELAY，NCYCL（在POST1中生成结构变形梯度线的动画） 31. ANCUT，NFRAM，DELAY，NCYCL，QOFF，KTOP，TOPOFF，NODE1，NODE2，NODE3（在POST1中生成等势切面云图动画） 32. ANDATA，DELAY，NCYCL，RSLTDAT，MIN，MAX，INCR，FRCLST，AUTOCNTRKY（生成某一
第5章通用后处理器(POST1) 静力分析 5.1概述 使用POST1通用后处理器可观察整个模型或模型的一部分在某一时间点（或频率）上针对指定载荷组合时的结果。POST1有许多功能，包括从简单的图象显示到针对更为复杂数据操作的列表，如载荷工况的组合。 要进入ANSYS通用后处理器，输入/POST1命令（Main Menu>General Postproc）. 5.2将数据结果读入数据库 POST1中第一步是将数据从结果文件读入数据库。要这样做，数据库中首先要有模型数据（节点，单元等）。若数据库中没有模型数据，输入RESUME命令(Utility Menu>File>Resume Jobname.db)读入数据文件Jobname.db。数据库包含的模型数据应该与计算模型相同，包括单元类型、节点、单元、单元实常数、材料特性和节点座标系。 注：数据库中被选来进行计算的节点和单元组应和模型中的节点和单元组属于相同组，否则会出现数据不匹配。有关数据不匹配的详细资料见5.2.2.3章。 一旦模型数据存在数据库中，输入SET,SUBSET或APPEND命令均可从结果文件中读入结果数据。 5.2.1 读入结果数据 输入SET命令（Main Menu>General PostProc>datatype）,可在一特定的载荷条件下将整个模型的结果数据从结果文件中读入数据库，覆盖掉数据库中以前存在的数据。边界条件信息（约束和集中力）也被读入，但这仅在存在单元节点载荷或反作用力的情况下，详情请见OUTRES命令。若它们不存在，则不列出或显示边界条件，但约束和集中载荷可被处理器读入，而且表面载荷和体积载荷并不更新，并保持它们最后指定的值。如果表面载荷和体积载荷是使用表格指定的，则它们将依据当前的处理结果集，表格中相应的数据被读入。加载条件靠载荷步和子步或靠时间（或频率）来识别。命令或路径方式指定的变元可以识别读入数据库的数据。例如：SET，2，5读入结果，表示载荷步为2，子步为5。同理，SET，，，，，3.89表示时间为3.89时的结果（或频率为3.89，取决于所进行分析的类型）。若指定了尚无结果的时刻，程序将使用线性插值计算出该时刻的结果。 结果文件（Jobname.RST）中缺省的最大子步数为1000，超出该界限时，需要输入SET，Lstep，LAST引入第1000个载荷步，使用/CONFIG增加界限。 注：对于非线性分析，在时间点间进行插值常常会降低精度。因此，要使解答可用，务必在可求时间值处进行后处理。
FX MX UX ROTX VX AX ACLX OMGX TEMP RBFX RBMX RBUX RBRX RBVX RBOX PRESS DCURVE DCURVE Option,LCID,Par1,Par2 Option----ADD,DELE, LIST, PLOT LCID---- Par1,Par2---[ *DIM *DIM Par,Type,IMAX,JMAX,KMAX,Var1,Var2,Var3 ! Par--- Type--- ARRAY IMAX,JMAX,KMAX--- *SET *SET Par,V ALUE! Par--- V ALUE--- EDLOAD ~ EDLOAD Option, Lab, KEY, Cname, Par1,Par2,PHASE,LCID,SCALE,BTIME,DTIME Option---ADD,DELE,LIST Lab--- FX UY PRSSURE KEY--- PRESSURE KEY ID EDLCS CID Cname--- Par1,Par2--- PHASE--- 0= =1 =2 LCID--- SCALE--- BTIME DTIME--- GUI Preprofessor>LS-DYNA Options>Loading Options>Specify Loads Solution>Loading Options>specify Loads EDFPLOT EDFPLOT KEY KEY--- ON 1 OFF 0 GUI Main Menu>Preprofessor>LS-DYNA Options>Loading Options>Show Forces EDVEL
ANSYS最常用命令流+中文注释 VSBV, NV1, NV2, SEPO, KEEP1, KEEP2 —Subtracts volumes from volumes，用于2个solid相减操作，最终目的是要nv1-nv2=?通过后面的参数设置，可以得到很多种情况：sepo项是2个体的边界情况，当缺省的时候，是表示2个体相减后，其边界是公用的，当为sepo的时候，表示相减后，2个体有各自的独立边界。keep1与keep2是询问相减后，保留哪个体？当第一个为keep时，保留nv1,都缺省的时候，操作结果最终只有一个体，比如：vsbv,1,2,sepo,,keep,表示执行1-2的操作，结果是保留体2，体1被删除，还有一个1-2的结果体，现在一共是2个体（即1-2与2），且都各自有自己的边界。如vsbv,1,2,,keep,,则为1-2后，剩下体1和体1-2，且2个体在边界处公用。同理，将v换成a 及l是对面和线进行减操作！ mp,lab, mat, co, c1,…….c4 定义材料号及特性 lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens） ex: 弹性模量 nuxy: 小泊松比 alpx: 热膨胀系数 reft: 参考温度 reft: 参考温度 prxy: 主泊松比 gxy: 剪切模量 mu: 摩擦系数 dens: 质量密度 mat: 材料编号（缺省为当前材料号） co: 材料特性值，或材料之特性，温度曲线中的常数项 c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数 定义DP材料： 首先要定义EX和泊松比：MP，EX，MA T，…… MP，NUXY，MAT，…… 定义DP材料单元表（这里不考虑温度）：TB，DP，MA T 进入单元表并编辑添加单元表：TBDATA，1，C TBDATA，2，ψ TBDATA，3，…… 如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下：MP，EX，1，1E8 MP，NUXY，1，0.3 TB，DP，1 TBDATA，1，27 TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,deg VSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Type，是选择的方式，有选择（s）,补选（a），不选(u)，全选(all)、反选（inv）等,其余方式不常用 Item, Comp 是选取的原则以及下面的子项 如volu 就是根据实体编号选择， loc 就是根据坐标选取，它的comp就可以是实体的某方向坐标！ 其余还有材料类型、实常数等 MIN, VMAX, VINC，这个就不必说了吧！ ,例：vsel,s,volu,,14 vsel,a,volu,,17,23,2 上面的命令选中了实体编号为14，17，19，21，23的五个实体 VDELE, NV1, NV2, NINC, KSWP: 删除未分网格的体 nv1:初始体号 nv2:最终的体号 ninc:体号之间的间隔 kswp=0:只删除体 kswp=1:删除体及组成关键点,线面 如果nv1=all,则nv2,ninc不起作用 其后面常常跟着一条显示命令VPLO,或aplo,nplo,这个湿没有参数的命令，输入后直接回车，就可以显示刚刚选择了的体、面或节点，很实用的哦！ Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备 Type: S: 选择一组新节点（缺省） R: 在当前组中再选择 A: 再选一组附加于当前组 U: 在当前组中不选一部分 All: 恢复为选中所有 None: 全不选 Inve: 反向选择 Stat: 显示当前选择状态 Item: loc: 坐标 node: 节点号
一般的接触分类 (2) ANSYS接触能力 (2) 点─点接触单元 (2) 点─面接触单元 (2) 面─面的接触单元 (3) 执行接触分析 (4) 面─面的接触分析 (4) 接触分析的步骤： (4) 步骤1：建立模型，并划分网格 (4) 步骤二：识别接触对 (4) 步骤三：定义刚性目标面 (5) 步骤4：定义柔性体的接触面 (8) 步骤5：设置实常数和单元关键字 (10) 步骤六： (21) 步骤7：给变形体单元加必要的边界条件 (21) 步骤8：定义求解和载步选项 (22) 第十步：检查结果 (23) 点─面接触分析 (25) 点─面接触分析的步骤 (26) 点－点的接触 (35) 接触分析实例（GUI方法） (38) 非线性静态实例分析（命令流方式） (42) 接触分析 接触问题是一种高度非线性行为，需要较大的计算资源，为了进行实为有效的计算，理解问题的特性和建立合理的模型是很重要的。 接触问题存在两个较大的难点：其一，在你求解问题之前，你不知道接触区域，表面之间是接触或分开是未知的，突然变化的，这随载荷、材料、边界条件和其它因素而定；其二，大多的接触问题需要计算摩擦，有几种摩擦和模型供你挑选，它们都是非线性的，摩擦使问题的收敛性变得困难。
一般的接触分类 接触问题分为两种基本类型：刚体─柔体的接触，半柔体─柔体的接触，在刚体─柔体的接触问题中，接触面的一个或多个被当作刚体，（与它接触的变形体相比，有大得多的刚度），一般情况下，一种软材料和一种硬材料接触时，问题可以被假定为刚体─柔体的接触，许多金属成形问题归为此类接触，另一类，柔体─柔体的接触，是一种更普遍的类型，在这种情况下，两个接触体都是变形体（有近似的刚度）。 ANSYS接触能力 ANSYS支持三种接触方式：点─点，点─面，平面─面，每种接触方式使用的接触单元适用于某类问题。 为了给接触问题建模，首先必须认识到模型中的哪些部分可能会相互接触，如果相互作用的其中之一是一点，模型的对立应组元是一个结点。如果相互作用的其中之一是一个面，模型的对应组元是单元，例如梁单元，壳单元或实体单元，有限元模型通过指定的接触单元来识别可能的接触匹对，接触单元是覆盖在分析模型接触面之上的一层单元，至于ANSTS使用的接触单元和使用它们的过程，下面分类详述。 点─点接触单元 点─点接触单元主要用于模拟点─点的接触行为，为了使用点─点的接触单元，你需要预先知道接触位置，这类接触问题只能适用于接触面之间有较小相对滑动的情况（即使在几何非线性情况下） 如果两个面上的结点一一对应，相对滑动又以忽略不计，两个面挠度（转动）保持小量，那么可以用点─点的接触单元来求解面─面的接触问题，过盈装配问题是一个用点─点的接触单元来模拟面─与的接触问题的典型例子。 点─面接触单元 点─面接触单元主要用于给点─面的接触行为建模，例如两根梁的相互接触。 如果通过一组结点来定义接触面，生成多个单元，那么可以通过点─面的接触单元来模拟面─面的接触问题，面即可以是刚性体也可以是柔性体，这类接触问题的一个典型例子是插头到插座里。
ansys学习－耦合与约束方程 1 耦合 当需要迫使两个或多个自由度取得相同（但未知）值，可以将这些自由度耦合在一起。耦合自由度集包含一个主自由度和一个或多个其它自由度。 典型的耦合自由度应用包括： "模型部分包含对称； "在两重复节点间形成销钉、铰链、万向节和滑动连接； "迫使模型的一部分表现为刚体。 如何生成耦合自由度集 1.在给定节点处生成并修改耦合自由度集 命令：CP GUI: Main Menu>Preprocessor>Coupling / Ceqn>Couple DOFs 在生成一个耦合节点集之后，通过执行一个另外的耦合操作（保证用相同的参考编号集）将更多节点加到耦合集中来。也可用选择逻辑来耦合所选节点的相应自由度。用CP命令输入负的节点号来删除耦合集中的节点。要修改一耦合自由度集（即增、删节点或改变自由度标记）可用CPNGEN命令。（不能由GUI直接得到CPNBGEN命令）。 2.耦合重合节点。 CPINTF命令通过在每对重合节点上定义自由度标记生成一耦合集而实现对模型中重合节点的耦合。此操作对“扣紧”几对节点（诸如一条缝处）尤为有用。 命令：CPINTF GUI: Main Menu>Preprocessor>Coupling / Ceqn>Coincident Nodes 3.除耦合重复节点外，还可用下列替换方法迫使节点有相同的表现方式： o如果对重复节点所有自由度都要进行耦合，常用NUMMRG命令（GUI：Main Menu>Preprocessor>Numbering Ctrls>Merge Items）合并节点。 o可用EINTF命令（GUI：Main Menu> Preprocessor>Create> Elements >At Coincid Nd）通在重复节点对之间生成2节点单元来连接它们。 o用CEINTF命令（GUI：Main Menu>Preprocessor> Coupling/Ceqn >Adjacent Regions）将两个有不相似网格模式的区域连接起来。这项操作使一个区域的选定节点与另一个区域的选定单元连接起来生成约束方程。 生成更多的耦合集 一旦有了一个或多个耦合集，可用这些方法生成另外的耦合集： 1.用下列方法以相同的节点号但与已有模式集不同的自由度标记生成新的耦合集。 命令：CPLGEN GUI: Main Menu>Preprocessor>Coupling / Ceqn>Gen w/Same Nodes 2.用下列方法生成与已有耦合集不同（均匀增加的）节点编号但有相同的自由度标记的新的耦合集： 命令：CPSGEN GUI: Main Menu>Preprocessor>Coupling / Ceqn>Gen w/Same DOF 使用耦合注意事项 1.每个耦合的节点都在节点坐标系下进行耦合操作。通常应当保持节点坐标系的一致性。 2.自由度是在一个集内耦合而不是集之间的耦合。不允许一个自由度出现在多于一个耦合集中。 3.由D或共它约束命令指定的自由度值不能包括在耦合集中。
161. EMF（电磁场分析中计算沿路径的电动势和电压降） 162. EMID，Key，Edges（增加或删除中间节点） 163. EMODIF，IEL，STLOC，I1，I2，I3，I4，I5，I6，I7，I8（调整单元坐标系方向）164. EMORE，Q，R，S，T，U，V，W，X（单元节点超过个时，在E命令后使用）165. EMUNIT, Lab, V ALUE(定义磁场单位) 166. EN，IEL，IJ，K，L，M，N，O，P（通过节点生成指定单元） 167. ENGEN，IINC，ITIME，NINC，IEL1，IEL2，IEINC，MINC，TINC，RINC，CINC，SINC，DX，DY，DZ（元素复制：用户自己进行编号） 168. ENORM，ENUM（重新定义壳单元的法线方向） 169. ENSYM，IINC，--，NINC，IEL1，IEL2，IEINC（镜像生成新单元：用户自己进行编号） 170. EPLOT（元素显示） 171. ERASE（擦除当前图形窗口显示的内容） 172. EREFINE，NE1，NE2，NINC，LEVEL，DEPTH，POST，RETAIN（将单元附近的单元网格细化） 173. ERESX，Key（控制单元积分点解的外推方式） Key=DEFA（线形材料单元节点解由积分点解外推得到） YES（节点解由积分点解外推得到） NO（节点解由积分点解拷贝得到） 174. ERNORM，Key（定义是否进行误差估计） 175. ERRANG，EMIN，EMAX，EINC（从文件读入单元数据） 176. ESEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS（选择单元子集） 177. /ESHAPE，SCALE（显示单元形状） 178. ESIZE，SIZE，NDIV（指定线划分单元的默认数目） 179. ESLA, Type（选择已选面上的单元） 180. ESLL, Type（选择已选线上的单元） 181. ESLN, Type, EKEY, NodeType（选择已选节点上的单元） 182. ESORT，Item，Lab，ORDER，KABS，NUMB（对单元数据指定新的排序方式）183. ESURF，XNODE，Tlab，Shape（在既有单元表面生成表面单元） 184. ESYM，--，NINC，IEL1，IEL2，IEINC（镜像生成新单元：自动编号） 185. ESYS，KCN（定义单元坐标系。【注】只能通过局部坐标系定义） 186. ET，ITYPE，Ename，KOPT1，KOPT2，KOPT3，KOPT4，KOPT5，KOPT6，INOP R（定义单元） 【注】KOPT1～KOPT6为元素特性编码，BEAM3的KOPT6=1时，表示分析后的结果可输出节点的力或力矩。 187. ETABLE，Lab，Item，Comp（将单元某项结果作成表格） 【注】Lab为字段名，最多8个字符；Item，Comp分别为单元输出表中的名称和分量。
面与面接触实例：插销拨拉问题分析 定义单元类型 Element/add/edit/delete 定义材料属性 Material Props/Material Models Structural/Linear/Elastic/Isotropic 定义材料的摩擦系数 … 建立几何模型 Modeling/Create/Volumes/Block/By Dimensions X1=Y1=0，X2=Y2=2，Z1=，Z2=
Modeling/Create/Volumes/Cylinder/By Dimensions Modeling/Operate/Booleans/Subtract/Volumes 先拾取长方体，再拾取圆柱体。 Modeling/Create/Volumes/Cylinder/By Dimensions 、 划分掠扫网格 Meshing/Size Cntrls/ManualSize/Lines/Picked Lines 拾取插销前端的水平和垂直直线，输入NDIV=3再拾取插座前端的曲线，输入NDIV=4
PlotCtrls/Style/Size and Shape，在Facets/element edge列表中选择2 facets/edge 建立接触单元 : Modeling/Create/Contact pair，弹出Contact Manager对话框，如图所示。 单击最左边的按钮，启动Contact Wizard（接触向导），如图所示。
单击Pick Target，选择目标面。 选择接触面 定义位移约束 施加对称约束，Define Loads/Apply/Structural/Displacement/Symmetric On Areas，选择对称面。 再固定插座的左侧面。 ） 设置求解选项 Analysis Type/Sol’s Control
260. L，P1，P2，NDIV，SPACE，XV1，YV1，ZV1，XV2，YV2，ZV2（定义线） 261. L2ANG，NL1，NL2，ANG1，ANG2，PHIT1，PHIT2（生成直线与两直线均成一定角度）262. L2TAN，NL1，NL2（生成直线与两直线均相切） 263. LANG，NL1，P3，ANG，PHIT，LOCAT（生成直线与已知直线成一定角度） 264. LARC，P1，P2，PC，RAD（生成弧线） 265. /LARC，XCENTER，YCENTER，XLRAD，ANGLE1，ANGLE2（使用弧线注释） 266. LAREA，P1，P2，NAREA（在面上两关键点之间生成一条最短的线） 267. LATT，MAT，REAL，TYPE，--，KB，KE，SECNUM（指定线的单元属性） 268. LCABS，LCNO，KABS（指定是否对载荷工况取绝对值） 269. LCASE，LCNO（将载荷工况读入） 270. LCDEF，LCNO，LSTEP，SBSTEP，KIMG（从结果文件中定义载荷工况） LCDEF，LCNO，ERASE（删除一载荷工况） 271. LCFACT，LCNO，FACT（指定载荷工况的比例因子） 272. LCFILE，LCNO，Fname，Ext，--（从载荷工况文件中定义载荷工况） 273. LCLEAR，NL1，NL2，NINC（清除线单元网格） 274. LCOMB，NL1，NL2，KEEP（线线合并） 275. LCOPER，Oper1，LCASE，Oper2，LCASE2（载荷工况的组合运算） 【注】Oper1=ADD（加）、SUB（减）、SQUA（平方）、SQRT（平方根）、SRSS（平方和求平方根）、MIN（比较存储最小值）、MAX（比较存储最大值）。 LCOPER，LPRIN（重新计算线单元的主应力） 276. LCSEL，Type，LCMIN，LCMAX，LCINC（选择载荷工况） 277. LCWRITE，LCNO，Fname，Ext，--（将当前载荷工况写入载荷工况文件中） 278. LCZERO（清空数据库中以前的数据） 279. LDELE，NL1，NL2，NINC，KSWP（删除线） 【注】KSWP=0删除线但保留线上关键点、1删除线及线上关键点。 280. LDIV，NL1，RATIO，PDIV，NDIV，KEEP（将线分割） 281. LDRAG，NK1，NK2，…，NK6，NL1，NL2，…，NL6（将一组既有关键点按一定路径拖拉成线） 282. LDREAD，Lab，LSTEP，SBSTEP，TIME，KIMG，Fname，Ext，--（施加耦合场载荷）283. LESIZE，NL1，SIZE，ANGSIZ，NDIV，SPACE，KFORC，LAYER1，LAYER2，KYNDIV（指定所选线上单元数） 284. LEXTND，NL1，NK1，DIST，KEEP（将线延伸） 285. LFILLT，NL1，NL2，RAD，PCENT（两条相交线生成倒角） 286. LFSURF，SLINE，TLINE（在既有面单元上生成重叠的表面单元） LGEN,ITIME,NL1,NL2,NINC,DX,DY,DZ,KING,NOELEM,IMOVE(从一条线或多条线生成新的线) 287. LGLUE，NL1，NL2，…，NL8，NL9（线间相互粘接） 288. /LIGHT，WN，NUM，INT，XV，YV，ZV，REFL（为模型添加光源） 289. LINA，NL，NA（线面相交） 290. /LINE，X1，Y1，X2，Y2（使用线注释） 291. LINL，NL1，NL2，…NL8，NL9（被选线的交集） 292. LINP，NL1，NL2，…NL8，NL9（线集两两相交）
CERIG命令详解 2011-12-21 21:27:42| 分类：ANSYS | 标签：|字号大中小订阅 ansys学习－耦合与约束方程 1 耦合 当需要迫使两个或多个自由度取得相同（但未知）值，可以将这些自由度耦合在一起。耦合自由度集包含一个主自由度和一个或多个其它自由度。 典型的耦合自由度应用包括： "模型部分包含对称； "在两重复节点间形成销钉、铰链、万向节和滑动连接； "迫使模型的一部分表现为刚体。 如何生成耦合自由度集 1.在给定节点处生成并修改耦合自由度集 命令：CP GUI: Main Menu>Preprocessor>Coupling / Ceqn>Couple DOFs 在生成一个耦合节点集之后，通过执行一个另外的耦合操作（保证用相同的参考编号集）将更多节点加到耦合集中来。也可用选择逻辑来耦合所选节点的相应自由度。用CP命令输入负的节点号来删除耦合集中的节点。要修改一耦合自由度集（即增、删节点或改变自由度标记）可用CPNGEN命令。（不能由GUI直接得到CPNBGEN命令）。 2.耦合重合节点。 CPINTF命令通过在每对重合节点上定义自由度标记生成一耦合集而实现对模型中重合节点的耦合。此操作对“扣紧”几对节点（诸如一条缝处）尤为有用。 命令：CPINTF GUI: Main Menu>Preprocessor>Coupling / Ceqn>Coincident Nodes 3.除耦合重复节点外，还可用下列替换方法迫使节点有相同的表现方式： o如果对重复节点所有自由度都要进行耦合，常用NUMMRG命令（GUI：Main Menu>Preprocessor>Numbering Ctrls>Merge Items）合并节点。 o可用EINTF命令（GUI：Main Menu> Preprocessor>Create> Elements >At CoincidNd）通在重复节点对之间生成2节点单元来连接它们。 o用CEINTF命令（GUI：Main Menu>Preprocessor> Coupling/Ceqn>Adjacent Regions）将两个有不相似网格模式的区域连接起来。这项操作使一个区域的选定节点与另一个区域的选定单元连接起来生成约束方程。 生成更多的耦合集 一旦有了一个或多个耦合集，可用这些方法生成另外的耦合集： 1.用下列方法以相同的节点号但与已有模式集不同的自由度标记生成新的耦合集。 命令：CPLGEN GUI: Main Menu>Preprocessor>Coupling / Ceqn>Gen w/Same Nodes 2.用下列方法生成与已有耦合集不同（均匀增加的）节点编号但有相同的自由度标记的新的耦合集：
第一天目标： 熟悉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 中定义接触通常有两种方法： 1. 用户自己手工创建接触单元和目标单元。这种方法，在定义接触和目标单元时还比较简单，但是在设置或修改单元属性和定义实常数时却比较复杂。需要用户对接触有较深刻的理解和通过实践积累丰富的经验。 2. 使用接触管理器中的接触向导定义接触对：使用接触管理器 (接触向导) 定义接触对(即接触单元和目标单元) 时，可以定义除了点-点接触以外的各种接触类型；它可以自动生成接触单元和目标单元，并提供了一组默认的单元属性和实常数值。使用这些默认的设置，加上适当的求解设置，对于多数接触问题都能够获得收敛的结果。而且，如果使用默认设置时，计算不收敛或对结果不太满意，也可以通过接触管理器(接触向导) 对单元属性和实常数方便的进行修改和调整。 因此，我们推荐，在可能的情况下，尽量使用接触管理器(接触向导) 来定义接触。本文将通过一个实例介绍接触管理器的基本使用方法。 所使用的例子如下： 两块平板，中间夹一个圆球。上面平板的上表面承受压力，分析模型的变形和应力随压力的变化。 两块平板，尺寸都是(100*100*20)，相距100。中间夹一个半径50 的圆球。两个平板分别与圆球的上下边缘接触。尺寸单位为mm。几何模型如图1。
图 1 中，为了能够划分映射网格，分别对体积进行了切割材料属性为：两块平板： E = 201000 Mpa；μ= 0.3 圆球： E = 70100 Mpa；μ= 0.33 接下来对各个Volumes 划分网格，单元类型采用solid186 (20 节点六面体)，单元边长统一取 6 mm。网格划分结果如图 2 所示：