ANSYS单元特性之命令流算例.doc

EX1.1 （LINK1）
（1）进入后处理模块,显示节点位移和杆件内力
MID_NODE = NODE (A/2,-B,0 ) ! 寻找距离位置（A/2,-B,0）最近的点，存入MID_NODE
*GET,DISP,NODE,MID_NODE,U,Y ！提取节点MID_NODE上的位移UY，若果已知要求的节点，直接提取即可。
LEFT_EL = ENEARN (MID_NODE) ! 需找距离节点MID_NODE最近的单元，存入LEFT_EL
ETABLE,STRS,LS,1 ! 用轴向应力SAXL的编号“LS,1”定义单元表STRS
*GET,STRSS,ELEM,LEFT_EL,ETAB,STRS ! 从单元表STRS中提取LEFT_EL单元的应力结果，存入变量STRSS。注意：提取的轴向应力结果具体到指定的单元。
（2）申明数组,提取计算结果,并比较计算误差
*DIM,LABEL,CHAR,2 ！定义2个元素的字符型数组LABEL
*DIM,VALUE,,2,3 ！定义2*3的数值型数组VALUE

LABEL(1) = 'STRS_MPa','DEF_mm' ! 给字符型数组的第1个元素赋值
*VFILL,VALUE(1,1),DATA,1,-0.05498 ! 给其他数值型数组中的元素赋值
*VFILL,VALUE(1,2),DATA,STRSS,DISP
*VFILL,VALUE(1,3),DATA,ABS(STRSS /1 ) ,ABS( DISP /0.05498 )

/OUT,EX1_1,out ！将输出内容重定向到文件EX1_1.out
/COM ! 以注释形式输出内容
/COM,------------------- EX1.1 RESULTS COMPARISON ---------------------
/COM,
/COM, | TARGET | ANSYS | RATIO
/COM,

*VWRITE,LABEL(1),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,' ',F10.3,' ',F10.3,' ',1F5.3)

/COM,----------------------------------------------------------------
/OUT ! 结束数据重定向，关闭输出文件
FINISH

*LIST,EX1_1,out ! 列表显示文件EX1_1.out的内容

EX1.2 （LINK1）
/PNUM, NODE,1 ！打开节点编号显示
/NUMBER, 2 ！只显示编号，不使用色彩

列表显示节点位移和单元的计算结果
PRDISP ! 列表显示节点位移值计算结果
ETABLE, MFORX,SMISC,1 ！以杆单元的轴力为内容，建立单元表MFORX
ETABLE, SAXL, LS, 1 ！以杆的轴向应力为内容，建立单元表SAXL
ETABLE, EPELAXL, LEPEL, 1 ! 以杆单元的轴向应变为内容，建立单元表EPELAXL
PRETAB ! 显示单元表中的计算结果
/NUMBER, 0 ！显示编号，并使用彩色
PLETAB, MFORX ！用色度图显示单元表MFORX中杆件轴力图

EX1.3 （LINK1）
NSEL,S,LOC,Y,1.0 ！选择所有位于Y=1.0位置上的节点
FSUM ！累计叠加选择集中所有节点上的反力
*GET,REAC_1,FSUM,,ITEM,FY ！将累加结果中的FY（Y方向的力）保存到变量REAC_1中

EX1.4 （LINK1）
R,1,65e-6 ！定义第1类实常数，杆件截面面积为65mm^2，在转化为国际单位制时操作

TREF,70 ! 设定参考温度为70度
BFUNIF,TEMP,80 ! 温度从原来的70度均匀上升到80度（TREF+10）

EX1.5 （PLANE42 AND CONTAC26）
ETABLE,STRSX,S,X ！定义X方向的应力为单元表STRSX
*GET,STRSSX,ELEM,3,ETAB,STRSX ！从单元表STRSX中提取3号单元的X向应力，存入STRSSX。注意：已知某一具体单元，在定义完单元表后，就可以直接用*GET语句提取这一具体单元的表中定义的项目。

EX2.1 （BEAM3）
NPLOT ! 显示节点位置图形，但是不显示节点号码
NPLOT, 1 ！显示节点位置图形，同时显示节点号码
NLIST ！列表显示节点在直角坐标系下的坐标值
DSYS, 1 ！改变显示坐标系，列出资料时，转换到圆柱坐标系
NLIST ！列表显示节点在圆柱坐标系下的坐标值

定义以两端弯矩和剪力为内容的单元表，并列出单元的单元表数据
ETABLE,IMOMENT, SMISC, 6 ！以单元I节点弯矩为内容，定义单元表
ETABLE,JMOMENT, SMISC, 12 ！以单元J节点弯矩为内容，定义单元表
ETABLE, ISHEAR, SMISC, 2 ！以单元I节点剪力为内容，定义单元表
ETABLE, JSHEAR, SMISC, 8 ！以单元J节点剪力为内容，定义单元表
PRETAB ！列表显示单元表中单元的计算结果

用单元表数据绘制剪力图
/TITLE, SHEAR FORCE DIAGRAM ! 定义剪力图窗口标题
PLLS, ISHEAR, JSHEAR ！结构剪力分布图

用单元表数据绘制弯矩图
/TITLE, BENDING MOMENT DIAGRAM ！ 定义弯矩图窗口标题
PLLS, IMOMENT, JMOMENT ！结构弯矩分布图

EX2.2 （BEAM3）
E,1,2 ！过1,2节点定义第1个单元
EGEN,4,1,1 ！按照前面的单元模式，生成4次，每次两端节点编号增加1，生成其余3个单元。

SET,1,1 ！读入第1荷载步第1子步的计算结果
PRNSOL,U,COMP ！列表显示节点线位移和阶巅峰线位移的矢量和
PRNSOL,ROT,COMP ！列表显示节点转角和节点角位移的矢量和

EX2.3 （BEAM4）
ET,1,BEAM4,,,,,,1 ！定义单元类型为BEAM4，并控制在单元坐标系中输出弯矩

NSEL,S,,,5 ！选择编号不小于5的所有节点
DSYM,SYMM,X ！定义刚选定的所有节点关于X轴的对称位移边界条件

SFBEAM,ALL,1,PRES,314 ！在所有梁单元表面施加均布荷载

进入后处理模块，显示位移计算结果，提取左边支座节点的转角和跨中挠度
/POST1
NSEL,S,,,1,5,4 ！选中1号和5号节点
PRNSOL,U,Z ！列表显示选择集中所有节点的Z向位移
PRNSOL,ROT,Y ！列表显示选择集中的所有节点的绕Y轴的转角ROTY
NSEL,ALL ！重新选中所有节点
PRRSOL ！列表显示被约束的节点的支反力（Print Reaction Solution）
RGHT_END = NODE (2.54,0,0) ！选择距离(2.54,0,0)最近的节点，并存入RGHT_END
LFT_END = NODE (0,0,0) ！选择距离(0,0,0)最近的节点，并存入LFT_END
*GET,UZ_MX_C2,NODE,RGHT_END,U,Z ！获取节点RGHT_END的线位移UZ
*GET,SLOPE_C2,NODE,LFT_END,ROT,Y ！获取节点LFT_END的转角ROTY
FINISH

在时间历程后处理器中，叠加处理5号节点的弯矩
/POST26
RFORCE,2,RGHT_END,M,Y ！叠加RGHT_END节点的弯矩MY，存入第2个POST26变量
STORE ！保存变量
*GET,M_MX_C2,VARI,2,EXTREM,VMAX ！获取第2个POST26变量的最大值
FINISH

EX2.4 （SOLID5，SOLID95，SOLID98，SOLID147）
smrt,off ! 关闭智能化功能
LESIZE,ALL,,,10 ! 定义所有线的划分段数为10段
ESIZE,25.4 ! 定义默认的单元尺寸为25.4
V,1,2,4,3,5,6,8,7 ! 用关键点连接成体
VMESH,1 ! 执行体单元划分操作

FK,2,FX,4.4483/4 ! 在2号关键点上施加沿着X方向的，大小为4.4483/4的轴向力
*REPEAT,4,2 ! 将前面的命令执行4次，每次操作对象的编号增加2，即完成对2、4、6关键点的荷载定义

FKDELE,ALL,FY ! 删除所有关键点上的轴向力

定义宏程序，完成自由端位移计算结果的处理
*CREATE,MAC ! 定义宏程序计算这些结果
SET,ARG1,1 ! 定义参数ARG1=1
LCDEF,ARG1,ARG1 ! 基于计算结果定义荷载工况
NSEL,S,LOC,X,152.4 ! 选择节点
PRNSOL,U,COMP ! 打印位移结果及其节点位移矢量和
LCFACT,ARG1,ARG2 ! 将缩放因子“AGR2”应用给荷载工况1
LCASE,ARG1 ! 将工况1读入到数据库
PRNSOL,U,COMP ! 打印位移结果及其节点位移矢量和
*END ! 结束宏程序的定义
使用上面定义的宏程序处理所有的3个工况
/COM *** !
*USE,MAC,1,(1/7.6E-4) ! 使用宏程序MAC，使用参数1，(1/7.6E-4)
SET,1,1 ! 读入第1荷载步的计算结果
LCSEL,S,0,0 ! 选择荷载步的子集
N1=NODE(152.4 , 0, 0) ! 选择距离(152.4 , 0, 0)位置最近的节点，命名为N1
*GET,UX1,NODE,N1,U,X ! 获取节点N1的位移UX
处理第2个荷载步的计算结果
*USE,MAC,2,(1/2.4712) ! 使用宏MAC，使用参数2，(1/2.4712)
SET,2,1 !
LCSEL,S,0,0 !
*GET,UY1,NODE,N1,U,Y !
处理第3个荷载步的计算结果
*USE,MAC,3,(1/10.9646) !
SET,3,1 !
LCSEL,S,0,0 !
*GET,UZ1,NODE,N1,U,Z !
定义字符串变量和相应的理论计算结果
*DIM,LABEL,CHAR,3,2 !
*DIM,VALUE,,3,3 !
LABEL(1,1) = 'DEFL X ','DEFL Y ','DEFL Z ' !
LABEL(1,2) = 'mm','mm','mm' !
*VFILL,VALUE(1,1),DATA,7.6E-4,2.471,10.965 !
*VFILL,VALUE(1,2),DATA,UX1,UY1,UZ1 !
*VFILL,VALUE(1,3),DATA,ABS(UX1/(7.6E-4)),ABS(UY1/2.471),ABS(UZ1/10.965)
SAVE, TABLE_1 ! 保存数据到数据库TABLE_1
FINISH

NSEL,S,LOC,X,152.4
NSEL,R,LOC,Y,2.54
NSEL,R,LOC,Z,1.27
*GET,MIDD,NDMX ! 获取当前最大节点号，并存入MIDD

从第1个表中恢复数据，绘制单元应力云图
RESUME,TABLE_1 ! 从数据库中恢复数据
/COM *** CLIPPED AND CAPPED DISPLAY OF STRESS CONTOURS ***
NSEL,ALL ! 选择所有节点
/VIEW,1,2,1,1 ! 设置观察方向为从（2,1,1）指向原点
EPLOT ! 绘制单元图
/TYPE,1,CAP ! 显示梁的剖面图
/DIST,1,5.08 ! 设置观察距离为5.08mm
/FOCUS,1,.3,.15,.09 ! 设置截开截面的位置
/Auto ! 以最佳比例显示
PLNSOL,S,X ! 显示X方向的应力

删除宏文件，临时表文件TABLE_*
/DELETE,MAC ! 删除宏文件
/DELETE,TABLE_1 ! 删除临时性数据库

EX2.5 （PIPE16）
F,2,MZ,1112*914 ! 定义作用在节点2上的绕Z轴大小为M=FD=1112*914N.MM的力偶
F,2,FX,-1112 ! 定义作用在2号节点的沿X轴大小为1112N的力

EX2.6 （SHELL63，BEAM4）
定义耦合自由度，以保证变形对称
CP,1,UZ,2,12 ! 定义 2号节点和12号节点的UZ自由度为耦合自由度，定义为集合1
CP,2,ROTY,2,12 ! 定义 2号节点和12号节点的ROTY自由度为耦合自由度，定义为集合1
CPSGEN,6,1,1,2 ! 按照前面定义的耦合自由度集合1到2，循环6次生成其余的耦合自由度

ETABLE,STRS,S,1 ！用壳单元SHELL63的最大拉应力S1为内容定义单元表STRS
ESORT,STRS ！按最大拉应力对单元表STRS的排序
*GET,SMAX1,SORT,,MAX ！获取排序后的最大拉应力值，存入变量SMAX1
PRNSOL,DOF ！列表显示节点位移计算结果

*DO,I,1,10 ！循环10次定义不同位置的单元不同的截面几何参数
R,I,1,1,(I-1)*76.2*12.7*12.7*12.7/120,6.35,1 ！REAL,1,AREA1,IZ1，IY1，TKZB1，TKYB1,IX1
RMORE,1,1,I*76.2*12.7*12.7*12.7/120,6.35,1 ！AREA2,IZ2，IY2，TKZB2，TKYB2,IX2
*ENDDO
RMODIF,1,3,76.2*12.7*12.7*12.7/120/2 ! 修改实常数
RLIST ！列出实常数

/NOPR ！关闭输入显示
RESUME,TABLE_2 从数据库中恢复数据
/GOPR ！打开输入显示

EX2.7 （BEAM54）
*GET,UY,NODE,1,U,Y ! 提取1号节点的Y方向的线位移UY，并保存到变量UY
*status,parm ！显示内存变量的内容

EX2.8 （BEAM54）
CNVTOL,F,,1E-4 ！设置荷载力的收敛误差为1E-4N
CNVTOL,M,,1E-1 ! 设置集中力偶的收敛误差为0.1N.M

END_NODE = NODE (0,3048/2,0) ! 选择距离(0,3048/2,0)最近的节点，命名为END_NODE
*GET,DEF,NODE,END_NODE,U,X ！提取节点END_NODE的位移UX，存入变量DEF
ETABLE,TENS,NMISC,1 ！以BEAM54的最大拉应力（NMISC,1）为内容定义单元表TENS
ETABLE,COMP,NMISC,2 ！以BEAM54的最大拉应力（NMISC,2）为内容定义单元表COMP
*GET,STS_TENS,ELEM,1,ETAB,TENS ！提取1号单元的最大拉应力，存入变量STS_TENS
*GET,STS_COMP,ELEM,1,ETAB,COMP ！提取1号单元的最大压应力，存入变量STS_COMP

EX3.1 （BEAM3）
*GET,FCR,MODE,1,FREQ ！提取第1阶频率计算结果，并保存到变量FCR.。注意：在屈曲分析中，此处的FCR代表的是压杆的临界失稳荷载。
*status,parm ！显示内存变量

EX3.3 （BEAM3）
FCR=-171.5 ！定义临界荷载
F,11,FY,FCR*1.015 ！定义荷载
F,11,FX,1 ！定义一个很小的扰动力
SOLVE
增加荷载,继续求解
FDEL,11,FX ！删除荷载
F,11,FY,FCR*1.063 ！定义荷载

/POST1
/USER ！定义用户自定义的焦点和观察距离
/FOCUS,,1270,1270 ！焦点定义命令，设置所有图形对象居中
/DIST,,1397 ！设置观察距离
/DSCALE,,1 ！按照模型实际比例显示
SET,1,0 ！读入第1荷载步的计算结果
PLDISP,1 ！显示原来模型图，显示变形后的模型
/NOERASE ！不清楚屏幕，重叠绘制下一个图形
SET,2,0 ！读入第2荷载步的计算结果
PLDISP

EX3.4 （BEAM188，BEAM189）
SECNUM,1 ！设置当前截面编号为1
SECTYPE,1,BEAM,ASEC ！定义梁的截面形状为任意形状，截面几何数据由用户提供
SECDATA,11613,5.61868e8,0.0,224747,0.0,898989.1 ！截面几何参数
以上定义的截面几何参数实际上是实常数
ET,1,BEAM188
MP,EX,1,491.2105
MP,NUXY,1,0.3
*DO,I,1,20,1 ！通过循环来建立单元，其中节点101为控制主轴方向的点
E,I,I+1,101
*ENDDO
*DO,I,21,40,1
E,I,I+1,201
*ENDDO

/OUT,SCRATCH, ! 将输出定向到啊文件SCRATCH，以避免屏幕显示
SOLVE
/OUT ！结束数据重定向
……
/VIEW,1,1,2,3 ！切换视点到（1,2,3）
/ANG,1 ！旋转模型到默认位置
/PBC,F,,1 ！显示力的边界条件
/PBC,U,,1 ！显示位移的边界条件
/PBC,ROT,,1 ！显示转角的边界条件
/ESHAPE,70 ！按照缩放因子70，以实体特征显示模型
EPLOT
/ESHAPE,0 ！关闭单元实体特征显示
/PBC,ALL,OFF ！关闭所有边界条件的显示
/OUT,SCRATCH,,,APPEND ！将中间计算过程输出到草稿文件SCRATCH
SOLVE
/OUT ！结束数据重定向，关闭输出文件

设置求解选项,施加位移约束后求解面内力作用发生平面弯曲的情况
OUTRES,ALL,ALL
ARCLEN,ON ！激活弧长法，默认状态时不激活
ARCTRM,U,1524,41,UZ ！设置弧长法的停止条件：41号节点的UZ向位移超过1524mm
F,41,FX,6.6057
NSUBST,10 ！指定荷载子步为10

NSOL,2,41,U,Z,DISP ！指定2号变量保存41号节点的UZ位移计算结果
RFORCE,3,1,F,X,FORCE ！指定3号变量保存1号节点的X方向的支反力计算结果
PROD,4,3, , ,FORCE , , ,-1.0,1,1, ！将3号变量乘以-1，赋值给4号变量

SET,2,7,1 ！读入第2荷载步，第7子步的计算结果，缩放因子为1
NSEL,S,LOC,X,0
FSUM ！累计叠加节点力和力矩
*GET,CP1,FSUM,,ITEM,FX ！提取计算结果中X向的反力结果，存入变量CP1

设置绘图选项
/TRIAD,OFF ！关闭坐标系图标显示
/PLOPTS,LOGO,0 ！关闭ANSYS图标显示
/PLOPTS,INFO,2 ！关闭图例显示
/PLOPTS,WP,0 ！关闭工作平面显示
/RGB,INDEX,100,100,100,0 ！指定0号颜色的三原色配比为100:100:100

EX4.4 （PLANE42）
OUTPR,NSOL,1 ！输出第1荷载步的节点计算结果
OUTPR,ESOL,1 ！输出第1荷载步的单元计算结果

EX4.6 （PLANE42）
CSYS,1 ！切换到柱坐标系
K,1,89 ！在坐标（89,0,0）处定义关键点1
K,2,89,90 ！在坐标（89,90,0）处定义关键点2
KGEN,2,1,2,1,25.4 ！循环2次，按照1号到2号节点，增量为1的节点范围，每个节点的半径坐标增量为25，生成3、4号节点。

NSEL,R,LOC,X,89 ！选择节点
NSEL,U,LOC,Y,90
PRNSOL,S,COMP ！列表显示已选节点应力计算结果
NSORT,S,Y ！将节点计算结果按曲梁内侧环向应力SY排序
*GET,SI,SORT,,MAX ！提取排序后的最大环向应力，存入变量SI

EX4.7 （PLANE42）
NSEL,S,LOC,Y,90 ！选择节点
CP,1,UY,ALL ！耦合所有选定的节点Y向位移，定义为集合1

EX4.8 （PLANE42）
LESIZE,3,,,4,.25 ！给3号线定义划分段数为4，间距比为0.25
LARC,4,5,6,12.7 ！过关键点4、5，以6号关键点为圆心，以12.7为半径定义圆弧

LSEL,S,LINE,,3,5,2 ！选择3号和5号线
DL,ALL,,SYMM ！对所有的选择集中的线施加对称位移约束

LSEL,S,LINE,,1 ！选择1号线
NSLL,,1 ！选择1号线上的节点
SF,ALL,PRES,-6.895 ！施加在所选节点上的拉力荷载为6.895MPa

NSORT,S,X,,,3 ！将计算结果按SX排序，只保留前3项最高结果
PRNSOL,S,COMP ！列表下式应力分量计算结果

/WINDOW,1,OFF ！关闭1号窗口的显示
/NOERASE ！不擦除图形窗口
/DSCALE,2,1 ！用缩放位移约束命令，缩放位移约束实部2倍，虚部1倍
/WINDOW,2,RTOP ！在图形窗口右上角定义窗口2
PLNSOL,S,X ！图形显示水平方向的节点应力SX分布图

LSEL,S,LINE,,1,2 ！选择1、2号线
NSLL,,1 ！选择1号线上的节点
NWRITE ！将选择的节点写入数据文件Jobanme.NODE

EX4.9 （PLANE42）
/WINDOW,2,-1,2,-2.5,0 ! 定义2号窗口

EX5.2 （SHELL51）
NSEL,S,LOC,X,0 ！选择节点
DSYM,SYMM,X ！施加关于X轴的对称边界条件

ESEL,S,ELEM,,1,1 ！选择单元
ETABLE,CENT,LS,5 ！定义单元表
ESORT,CENT ！排序单元表
*GET,PRSCNT,SORT,,MAX ！提取单元表中的最大值，存入变量
EX7.2 （SHELL51）
再次进入求解模块,进行模态分析
/SOLU ！
ANTYPE,MODAL
MODOPT,SUBSP,3 ！用子空间循环迭代方法计算前3阶振型
MXPAND,3 ！扩展前3阶模态
PSTRES,ON
DDELE,2,UX,13 ！释放2号节点到13号节点的所有X方向的自由度
DDELE,2,UY,13

EX8.1 （BEAM3）
在时间历程后处理器获取频率计算结果
/POST26
*GET,FREQ1,MODE,1,FREQ ! 获取第1阶频率的计算结果，并存入变量FREQ1
*GET,FREQ2,MODE,2,FREQ
*GET,FREQ3,MODE,3,FREQ

EX8.2 （BEAM3）
在一般后处理模块中,获取频率计算结果
/POST1 ！可以再/SOLU、/POST1、/POST26中提取频率
*GET,FREQ1,MODE,1,FREQ ！获取第1阶频率的计算结果，并存入变量FREQ1
*GET,FREQ2,MODE,2,FREQ
*GET,FREQ3,MODE,3,FREQ

EX8.4 （SHELL63）
OUTPR,ALL,1 ！指定输出选项为第1阶频率的所有计算输出参数
OUTRES,ALL,0 ！指定所有输出参数不用保存到数据库

EX9.3 （SHELL63）
CSYS,1
N,1,.0001
N,10,76.2e-3
FILL
NGEN,4,10,1,10,1,,10 ！循环4次，每个节点号增量为10，按照1到10号节点的格式，所有Y坐标（这里是柱坐标系的Theta）增加10度。
NROTAT,ALL ！旋转所有节点的的局部坐标系到当前柱坐标系
E,1,2,12,11 ！过节点1,2,12,11定义第1个单元
！EGEN, ITIME, NINC, IEL1, IEL2, IEINC, MINC, TINC, RINC, CINC, SINC, DX, DY, DZ
EGEN,3,10,-1 ！循环3次，每次单元的节点号增量为10，按照前面1个单元的模式生成单元
EGEN,9,1,-3 ！循环9次，每次单元的节点号增量为1，按照前3个单元的模式生成单元
NSEL,S,LOC,Y,0 ！选择所有Y=0（度）的直线边上的所有节点
CM,RIGHT,NODE ！将这些选中的节点组合成一个部件，叫做RIGHT

EX9.9 （PLANE42，SOLID45）
ET,1,Plane42
ET,2,Solid45
Esize,6.35 ！指定单元划分尺寸为6.35
AMESH,ALL ！将所有面划分为PLANE42单元
指定拉伸选项,通过体拉伸生成体单元
TYPE, 2 ! 切换单元类型为第2类单元SOLID45
EXTOPT,ESIZE,10,0, ！拉伸生成选项，生成10个体单元
EXTOPT,ACLEAR,0 ！拉伸时不清除平面单元的网格
EXTOPT,ATTR,0,0,0 ！保留体单元的特征类型不变
体拉伸命令格式：VEXT, NA1, NA2, NINC, DX, DY, DZ, RX, RY, RZ
VEXT,1, , ,0,0,524,,,, ！按照1号面，执行体拉伸，Z向增量了DZ=254

EX10.1 （PLANE82）
设置单元剖分的尺寸并进行单元剖分
AESIZE,ALL,5 ！划分单元尺寸设置为5mm
AMESH,ALL ！对所有面进行单元划分

EX11.1 （PLANE82）
/POST1 ！进入一般后处理模块
SET,LAST ！读入最后一个荷载步的计算结果
ETABLE,VOLU,VOLU ！用单元体积定义单元表VOLU
PRNSOL,S,PRIN ！列表显示主应力
NSORT,S,1 ！对单元表按第1主应力排序
NSEL,S,LOC,X,0,228.6 ！选择X=0和X=228.6范围内的节点
*GET,STRS,SORT,,MAX ！提取单元表排序后的最大值，即最大拉应力
NSEL,ALL
SSUM ！计算并显示单元表汇总结果，体积
*GET,TVOL,SSUM,,ITEM,VOLU ！将汇总后的总体积保存到变量TVOL
TVOL=TVOL*2 ！将体积乘以2换算成总体积
NSEL,S,LOC,X,253,255 ！选择X=235到X=255范围内的节点，即自由端端点
PRNSOL,U,Y ！列表显示节点位移UY值
NSORT,U,Y,,1 ！对节点计算结果UY排序
PRNSOL,U,Y ！列表显示节点位移UY值
*GET,DEFL,SORT,,MAX ！提取排序后的节点位移结果，保存到变量DEFL
*STATUS,PARM ！显示内存变量的值
DEFL=ABS(DEFL) ！将端部的最大挠度值取绝对值后，存入DEFL

L,1,6 ！过关键点1和6定义直线
*REPEAT,5,1,1 ！重复执行前面定义直线的命令，每次节点编号增加1

*USE,SCRATCH ！指定宏程序SCRATCH
/OPT ！进入优化设计模块
OPANL,SCRATCH ！指定优化循环文件SCRATCH
OPVAR,TVOL,OBJ,,,0.25 ！定义目标函数为总体积TVOL，允许容差为0.25单位
OPVAR,STRS,SV,,207 ！定义应力状态变量DEFL从0到207
OPVAR,TK16,DV,3.81,6.86,0.0254 ！定义设计变量TK16从3.81到6.86，容差0.0254
OPSAVE,INITIAL,OPT ！保存初始设计方案到文件INITIAL.OPT
OPTYPE,SUBP ! 使用子模型近似优化设计方法
OPSUBP,30 ！最大优化设计30个循环
用子模型近似优化设计方法开始优化设计的计算
OPEXE ！用子模型近似优化设计方法开始优化设计的计算
VR1=TVOL ！将体积TVOL保存到变量VR1
PARSAV,,RSET1 ！将所有参数变量保存到文件Jobname.RSET1
OPLIST,ALL,,1 ！显示设计集中的参数
/TITLE,EX11.3, SHAPE OPTIMIZATION OF A CANTILEVER BEAM (SubP Method)
/AXLAB,Y,VOLUME (TVOL)
PLVAROPT,TVOL ！图形显示目标函数收敛情况
STATUS,PARM ！显示内存变量的值
设置一阶方法进行优化设计的参数
OPRESU,INITIAL,OPT ！重新读入初始设计方案
OPVAR,DIF1,DEL ！为一阶方法删除保证锥度的约束DIF1
OPTYPE,FIRST ！指定使用一阶方法
OPFRST,20 ！最大循环次数为20
用一阶方法执行优化设计
OPEXE ！用一阶方法执行优化设计
VR4=TVOL ！将体积TVOL保存到变量VR4
OPLIST,ALL,,1 ！显示设计集中的参数
/TITLE,EX11.3, SHAPE OPTIMIZATION OF A CANTILEVER BEAM (1st Method)
/AXLAB,Y,VOLUME (TVOL)
PLVAROPT,TVOL ！图形显示目标函数收敛情况
STATUS,PARM
FINISH

网架结构优化
Etable,evolume,volu ! 定义单元体积列表
Ssum ! 单元体积列表求和 （计算并打印单元表中每列之和）。
*Get,volume,ssum,,item,evolume ! 定义变量volume存储网架结构消耗材料体积

Etable,smax_e,ls,1 ! 定义单元轴向应力列表
Esort,etab,smax_e,0,1 ! 单元轴向应力降幂排序
*Get,smaxe,sort,,max ! 定义变量smaxe存储网架单元中最大轴向应力

*get,dzmax,node,49,u,z ! 定义变量dzmax存储网架下平面中心49号节点位移UZ
*set,dmax,-1*dzmax ! 定义变量dmax存储网架下平面中心49号节点挠度

Lgwrite,’optimize’,’txt’,,comment

Opanl,’optimize’,’txt’,’’ ! 指定优化命令文件
!定义设计变量、状态变量和目标变量
Opvar,area_1,dv,20,200,0.01 ! 定义设计变量area
Opvar,smaxe,sv,200,210,0.01 ! 定义状态变量smax
Opvar,dmax,sv,0.1,10,0.01 ! 定义状态变量dmax
Opvar,volume,obj,,,20000 ! 定义目标变量volume
!指定优化方法
Optype,firs ! 指定优化方法
Opfrst,100,100,0.2 ! 优化方法设定
!优化分析
Opexe ! Initiates optimization looping.
!查看优化结果
*status

/input,’Jobname’,’MCOM’ ! 读入文件

查看模态求解信息
在ANSYS Output Window中可以查看模态计算时的求解信息。如果想把求解信息保存下来，则需要在求解Solve之前，将输出的信息写入文本中。操作如下：求解之前，在菜单Utility Menu中选择File-Switch Output to-File命令，打开Switch Output to File对话框，接着定义文件名，选择保存路径后，单击“OK”按钮创建文件，然后求解。
求解结束后，在菜单Utility Menu中选择File-Switch Output to-Output Window命令，使信息继续在输出窗口中显示，不再保存到创建的文件中。完整的求解信息中主要包含总质量、结构在各方向的总转动惯量、各种单元质量、各阶频率、周期、参与因数、参与比例及有效质量累积因数等。
/OUTPUT,'AAAA','txt','I:\',
ANTYPE,2
MODOPT,SUBSP,6,0,100, , OFF
SOLVE

/OUTPUT, TERM
声明：本套资料由本人总结概括，如果您在使用过程中发现本套资料有不当或错误之处请联系本人。本人联系QQ：63966955 另：本人空间中有大量ANSYS学习资料，空间地址：