提取计算结果.docx
《提取计算结果.docx》由会员分享，可在线阅读，更多相关《提取计算结果.docx（18页珍藏版）》请在冰豆网上搜索。
提取计算结果
Ansys中NSORT命令和ETABLE命令详解
用NSORT命令获取最大应力点数据
在ANSYS的结果后处理时，实体单元的应力结果显示时默认会用"MX"和"MIN"在图形中标出最大最小应力点位置，要知道最大（或最小）应力点的节点编号的话可通过PRNSOL命令列表查得。
一般来讲，结构在不同工况下的最大应力点一般不在同一位置，如果工况较多或者要反复查询的话，这种方法就需要人为地反复记录最大（或最小应力节点号），比较繁杂。
如果能够把最大（或最小）应力点的节点号动态地符给一个变量，通过APDL语言就可以随时求得最大（或最小）应力点的数据了。
方法1：
/POST1
NSEL,S,ALL
NNUM=NDINQR（0,13）!
当前选择集中的节点数目
*GET,N,NODE,,NUM,MIN!
获取选择集中最小节点号
N_MAX=N;!
当前应力最大节点号
*GET,SEQV_MAX,NODE,N,S,EQV!
当前最大应力
*DO,I,2,NNUM!
开始循环
*GET,TEMP,NODE,NDNEXT（N）,S,EQV
*IF,TEMP,GT,SEQV_MAX,THEN!
比较大小
N_MAX=NDNEXT（N）
SEQV_MAX=TEMP
*ENDIF
N=NDNEXT（N）!
下一个节点
*ENDDO
方法2：
上面的命令流，由于需要反复对每一个节点进行遍历比较，如果节点数量较大时运算时间会比较长。
其实，ANSYS中还提供了一个对节点数据进行排序的命令NSORT，借助它可以快速地实现上面命令流的功能，而且如果只查为了得到一个最大应力点的话，运算速度会很快。
命令流如下：
/POST1
NSEL,S,ALL
NSORT,S,EQV,0,,1,SELECT!
结果按应力大小进行排序
*GET,N_MAX,NODE,0,NUM,MAX!
取得应力最大点的节点编号
NUSORT!
恢复缺省的排序方式
*GET,SEQV_MAX,NODE,N_MAX,S,EQV!
取得最大应力值
关于NSORT命令的格式：
NSORT,Item,Comp,ORDER,KABS,NUMB,SEL
其中：
Item,Comp——排序的依据项标识，如S、LOC等；
ORDER——排序方式。
默认为0,按降序，若为1按升序；
KABS——是否按绝对值进行排序。
默认为0,不按绝对值排序，若为1按绝对值进行排序；
NUMB——排序后结果记录的个数。
默认为所有节点，此数越少排序速度越快；
SEL——控制结果结点是否可以进行再次选择。
如为空则不可选择（默认），如果排序后使用了NSEL等节点选择命令，则排序结果自动恢复为默认的排序方式；如为SELECT则当前节点选择集自动变为排序的结果记录结点，此时可用NSEL命令从当前选择集中进行进一步选择，如果节点不在排序后的选择集中，则也恢复默认排序方
用ANSYS的APDL方式直接获得最大应力和最大应变的节点编号及其数值的有好多...
用ANSYS的APDL方式直接获得最大应力和最大应变的节点编号及其数值的有好多种方法：
如下面的这种，先对节点的值进行SORT，在提取最大的值即可。
最大应力节点编号及其数值的提取：
ALLSEL
NSORT,S,EQV,0,0,ALL
*GET,MAX_SEQV,SORT,0,IMAX
*GET,MAXSEQV,NODE,MAX_EQV,S,EQV
依次类推，可以最大应变节点编号及其值数值的提取：
ALLSEL
NSORT,EPTO,EQV,0,0,ALL
*GET,MAX_EPTOEQV,SORT,0,IMAX
*GET,MAXEPTOEQV,NODE,MAX_EPTOEQV,S,EQV
最大位移节点编号及其数值的提取：
ALLSEL
NSORT,U,SUM,0,0,ALL
*GET,MAX_U,SORT,0,IMAX
*GET,MAXU,NODE,MAX_U,U,SUM
还有其他方法：
如把所有的节点应力应变数值提取然后找最大值
*GET,NUMALL,NODE,0,COUNT
*GET,NUMSTART,NODE,0,NUM,MIN
SMAX=0
INUSE=NUMSTART
*DO,I,1,NUMALL,1
*GET,SSUM,NODE,INUSE,S,EQV
*IF,SMAX,GE,SSUM,THEN
SMAX=SMAX
*ELSEIF,SMAX,LT,SSUM,THEN
SMAX=SSUM
*ENDIF
INUSE=NDNEXT（INUSE）
*ENDDO
还有一种是首先在窗口上显示应力应变云图，然后利用GET命令得到最大值。
如提取最大应力为：
PLNSOL,S,EQV,0
*GET,PAR,PLNSOL,0,MAX
或将所有的节点定义成一个路径，*GET,PAR,PATH,0,MAX，
另：
获取结果各个时间点的最大应力值
*DIM,MAXSEQVNODE,TABLE,63,1,1,,,!
定义数组
*DO,ISET1,1,63,1!
根据结果存储的数目而定,如63
/POST1
SET,1,,1,,,,ISET1,
PLNSOL,S,EQV,0
*GET,MAXSEQV_NODEI,PLNSOL,0,MAX!
获取结果,这里是最大MISES等效应力
*SET,MAXSEQV_NODE（ISET-1,0,1）,（ISET-1）*4E-4!
确定对应的时间
*SET,MAXSEQV_NODE（ISET-1,1,1）,MAXSEQV_NODEI!
写入最大等效应力值
*ENDDO
。
。
。
。
。
。
。
。
。
。
。
。
。
用ETABLE命令提取计算结果
1、ETABLE简介
原来ansys求解完问题之后，会把计算结果保存在一个表里面，表的行是单元的编号，表的列是单元的计算结果，如节点位移、节点力、应力、应变等等。
ETABLE这个命令就是把表中的某个列取出来，赋值给某个自定义的向量，再通过*GET命令可以指定某个具体的单元，就可以把该单元的对应计算结果提取出来了。
ETABLE,Lab,Item,Comp将单元值形成一个表以便进一步的处理。
2、命令选项解释：
1）Lab
用户自定义的表名，用于后续命令或输出的标题，最多可使用8个字母，不可与预定义的表名称重复。
默认的表名是Item和Comp项的前四个字母组合而成的8个字母。
如果与用户之前定义的表名相同，本次结果将被包括在同一表中。
最多可定义200个不同的表名。
以下表名是ANSYS预定义的，不可用作用户自定表名：
REFL,STAT,和ERAS.
Lab=REFL以ETABLE的最新选项重写所有ETABLE命令预定义的表，但保留字段将被忽略，这个命令在载荷步改变后重写表时很方便。
Lab=STAT将显示储存的表的值。
Lab=ERAS将删除整个表。
2）Item
选项名称。
常用的选项名称见后表。
某些选项需要栏目名。
Item=ERAS将删除表中的某一栏。
3）Comp
选项的栏目名（如果需要的话）。
常用的栏目名见后表。
4）说明
定义单元值的表以便后续处理。
单元值表可以被认为是工作表，其行代表所有被选择的单元。
其列代表通过ETABLE命令输入表中的单元值。
每一列数据有一个用户定义的名称，用于列表和显示。
将数据输入单元表后，你不仅可以列出和显示你的数据，还可以对数据进行许多操作，例如列相加或列相乘[SADD,SMULT]，为安全计算定义允许的应力[SALLOW]，或者将一列数据和另一列相乘[SMULT]。
更多的细节请看ANSYSBasicAnalysisGuide。
有很多不同类型的数据结果可以被存在单元表中。
例如，许多单元的选项只有一个值（也就是说，每一个单元对应此选项只有一个值）单一值的选项包括：
SERR,SDSG,TERR,TDSG,SENE,TENE,KENE,JHEAT,JS,VOLU和CENT.其余的选项是多个值的（也就是说，这些值在单元中是变化的，每个节点有不同的值）。
因为每个单元只能有一个值存在单元表中，多值的选项存入的是平均值（视节点数而定）。
例外的是FMAG和所有的单元力选项，它们存入的是相关节点值的和。
（这段话的意思是说，单值的单元选项，如单元体积，存入表中的就是这个值；而在单元不同位置有不同值的选项，如应力，写入表中的是单元的平均值。
根本原因在于一个单元只能对应表中的一个数据。
）
5）访问方法
ETABLE命令中可以使用两种数据访问方法，视你想储存的数据不同而不同。
一些结果只用通用名就可以访问（要素名法），而另一些结果需要一个标志名和标志数（序列数法）。
要素名法用于访问常用的单元数据（也就是说，绝大部分单元类型都有的数据）。
所有的单值选项以及一部分多值选项可以用要素名法访问。
不同的单元值视计算方法不同和选择集不同而不同。
（AVPRIN,RSYS,LAYER,SHELL,和ESEL）
尽管节点值不用单元表也可以很容易地列出和显示，你仍然可能需要利用单元表储存这些节点数据以便后需的操作。
要素名法的选项名和栏目名见后表。
序列数法可以使你访问那些非平均值的结果（例如节点的压力，连接点的温度等等），或者是那些不宜用普通格式描述的数据（例如结构线单元和接触单元的导出数据，热线单元、层单元的导出数据等）。
描述不同单元的这些选项（如LS,LEPEL,LEPTH,SMISC,NMISC,SURF等）和对应的序列数的表请见ANSYSElementsReference.
一些单元表数据是基于结果坐标系的，这些数据包括所有的要素结果（例如：
位移UX,UY，应变SX,SY等）。
求解器把要素结果依照求解坐标系写入数据库。
当你使用ETABLE的时候，这些结果在写入表中之前被转换为结果坐标系。
默认的结果坐标系是global坐标[RSYS,0]。
所有其他的数据在从数据库中提取出来写入表中时没有经过坐标转换。
6）后处理方法
使用PRETAB,PLETAB,orETABLE,STAT命令可以显示存储的表值。
使用ETABLE,ERAS命令删除整个表，使用ETABLE,Lab,ERAS命令删除表中名为Lab的栏。
在GUI界面下，如果对话框将某一单元表的DELETE命令写入日志文件（Jobname.LOGorJobname.LGW），你会发现明令行中的Lab为空缺，Item=ERASE，而Comp是一个整数。
在这种情况下，GUI给Comp指定了一个值，这个值对应于列表框中被选择的变量名，这并不表明你也要在ANSYS中给Comp输入这样一个值。
然而，包括如此由GUI产生的ETABLE命令的文件可以用作批处理输入或用于/INPUT命令。
7）ETABLE－常用选项和选项栏目表
Item选项
Comp选项栏目
Description描述
U
X，Y，Z
X,Y,Z方向的位移
ROT
X，Y，Z
X,Y,Z方向的旋转
TEMP
温度
PRES
压力
VOLT
电压
MAG
磁梯度位
V
X，Y，Z
X,Y,Z方向流体速度
A
X，Y，Z
X,Y,Z方向磁矢量差
EMF
电动势降
ENKE
紊流动能
ENDS
紊流能量损失
SP0n
核质量分率
FLOTRAN节点结果可用选项
TTOT
总温度
HFLU
热流量
HFLM
热传导系数
COND
层流传导率
PCOE
压力系数
PTOT
总压力
MACH
马赫数
STRM
流量函数
DENS
流体密度
VISC
层流粘性
EVIS
流体有效粘性
ECON
流体有效传导率
YPLU
Y+参数
TAUW
壁剪切应力
LMDn
核素片质量扩散率
EMDn
核素有效质量扩散率
单元结果可用选项
S
X,Y,Z,XY,YZ,XZ
Componentstress.
各方向应力
1,2,3
Principalstress.
主应力
INT
Stressintensity.
应力强度
EQV
Equivalentstress.
相当应力
EPEL
X,Y,Z,XY,YZ,XZ
Componentelasticstrain.
各方向弹性应变
1,2,3
Principalelasticstrain.
主弹性应变
INT
Elasticstrainintensity.
弹性应变强度
EQV
Elasticequivalentstrain.
弹性等效应变
EPTH
X,Y,Z,XY,YZ,XZ
Componentthermalstrain.
各方向热应变
1,2,3
Principalthermalstrain.
主热应变
INT
Thermalstrainintensity.
热应变强度
EQV
Thermalequivalentstrain.
热相当应变
EPPL
X,Y,Z,XY,YZ,XZ
Componentplasticstrain.
各方向塑性应变
1,2,3
Principalplasticstrain.
主塑性应变
INT
Plasticstrainintensity.
塑性应变强度
EQV
Plasticequivalentstrain.
塑性相当应变
EPCR
X,Y,Z,XY,YZ,XZ
Componentcreepstrain.
各方向蠕变应变
1,2,3
Principalcreepstrain.
主蠕变应变
INT
Creepstrainintensity.
蠕变应变强度
EQV
Creepequivalentstrain.
蠕变相当应变
EPSW
Swellingstrain.
膨胀应变
EPTO
X,Y,Z,XY,YZ,XZ
Componenttotalmechanicalstrain（excludingthermal）（EPEL+EPPL+EPCR）.
各方向总机械应变（不包括热应变），即EPEL+EPPL+EPCR
1,2,3
Principaltotalmechanicalstrain.
总主机械应变
INT
Totalmechanicalstrainintensity.
总机械应变强度
EQV
Totalequivalentmechanicalstrain.
总机械相当应变
EPTT
X,Y,Z,XY,YZ,XZ
Componenttotalstrainincludingthermal（EPEL+EPTH+EPPL+EPCR）.
各方向总应变（包括热应变），即EPEL+EPTH+EPPL+EPCR
1,2,3
Principaltotalstrain.
总主应变
INT
Totalstrainintensity.
总应变强度
EQV
Totalequivalentstrain.
总相当应变
NL
SEPL
Equivalentstress（fromstress-straincurve）.
相当应变（基于应力－应变图）
SRAT
Stressstateratio.
应力状态比率
HPRES
Hydrostaticpressure.
流体静力学压力
EPEQ
Accumulatedequivalentplasticstrain.
累积塑性相当应变
PSV
Plasticstatevariable.VISCO106,VISCO107,andVISCO108only.
塑性状态变量（仅适用于VISCO106,VISCO107,andVISCO108）
PLWK
Plasticwork/volume.VISCO106,VISCO107,andVISCO108only.
塑性体积（仅适用于VISCO106,VISCO107,andVISCO108）
SEND
ELASTIC
Elasticstrainenergydensity.
弹性应变能量密度
PLASTIC
Plasticstrainenergydensity.
塑性应变能量密度
CREEP
Creepstrainenergydensity.
蠕变应变能量密度
TG
X,Y,Z,SUM
Componentthermalgradientorvectorsum.
各方向热梯度或矢量和
TF
X,Y,Z,SUM
Componentthermalfluxorvectorsum.
各方向热通量或矢量和
PG
X,Y,Z,SUM
Componentpressuregradientorvectorsum.
各方向气体压力梯度或矢量和
EF
X,Y,Z,SUM
Componentelectricfieldorvectorsum.
各方向电场或矢量和
D
X,Y,Z,SUM
Componentelectricfluxdensityorvectorsum.
各方向电通量密度或矢量和
H
X,Y,Z,SUM
Componentmagneticfieldintensityorvectorsum.
各方向磁场强度或矢量和
B
X,Y,Z,SUM
Componentmagneticfluxdensityorvectorsum.
各方向磁感应密度或矢量和
FMAG4
X,Y,Z,SUM
Componentmagneticforcesorvectorsum.
各方向磁力或矢量和
SERR
Structuralerrorenergy.Anisotropicelasticmaterials（TB,ANEL）arenotsupported.
结构误差能量
SDSG
Absolutevalueofmaximumvariationofanynodalstresscomponent.
节点应力最大改变量的绝对值
TERR
Thermalerrorenergy.
热误差能量
TDSG
Absolutevalueofthemaximumvariationofanynodalthermalgradientcomponent.
节点热梯度最大改变量的绝对值
F
X,Y,Z
Componentstructuralforce.Sumofelementnodalvalues.
各方向力，为单元各节点力的和
M
X,Y,Z
Componentstructuralmoment.Sumofelementnodalvalues.
各方向力矩，为单元各节点力矩的和
HEAT
Heatflow.Sumofelementnodalvalues.
热流速，为单元各节点热流速的和
FLOW
Fluidflow.Sumofelementnodalvalues.
流体流速，为单元各节点流速的和
AMPS
Currentflow.Sumofelementnodalvalues.
电流，为单元各节点电流的和
FLUX
Magneticflux.Sumofelementnodalvalues.
磁通量，为单元各节点磁通量的和
VF
X,Y,Z
Componentfluidforce.
各方向流体力
CSG
X,Y,Z
Componentmagneticcurrentsegment.
磁流片段
SENE
"Stiffness"energyorthermalheatdissipation（appliestoallelementswheremeaningful）.SameasTENE.
硬度能量或热消散（取决于单元），同TENE
AENE
Artificialenergyoftheelement.Thisincludesthesumofhourglasscontrolenergyandenergygeneratedbyin-planedrillingstiffnessfromshellelements（appliestoallelementswheremeaningful）.TheenergyisusedforcomparisonstoSENEenergytopredictthesolutionerrorduetoartificialstiffness.SeetheTheoryReferenceforANSYSandANSYSWorkbench.
单元的人工能量
TENE
Thermalheatdissipationor"stiffness"energy（appliestoallelementswheremeaningful）.SameasSENE.
热量消散或刚度能量。
同SENE
KENE
Kineticenergy（appliestoallelementswheremeaningful）.
动能（使用所有有意义的单元）
STEN
Elementalenergydissipationduetostabilization.
JHEAT
ElementJouleheatgeneration.
单元焦耳热
JS
X,Y,Z,SUM
Sourcecurrentdensityforlow-frequencymagneticanalyses.Totalcurrentdensity（sumofconductionanddisplacementcurrentdensities）inlowfrequencyelectricanalyses.Components（X,Y,Z）andvectorsum（SUM）.
低频磁分析的源电流密度。
低频磁分析的总源电流密度（SUM）
JT
X,Y,Z,SUM
Totalmeasureablecurrentdensityinlow-frequencyelectromagneticanalyses.（Conductioncurrentdensityinalow-frequencyelectricanalysis.）Components（X,Y,Z）an
升级会员 