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_ELETABLE,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的数值型数组VALUELABEL(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 ！以杆单元的轴力为内容，建立单元表MFORXETABLE, SAXL, LS, 1 ！以杆的轴向应力为内容，建立单元表SAXLETABLE, EPELAXL, LEPEL, 1 ! 以杆单元的轴向应变为内容，建立单元表EPELAXLPRETAB ! 显示单元表中的计算结果/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类实常数，杆件截面面积为65mm2，在转化为国际单位制时操作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 ！在所有梁单元表面施加均布荷载进入后处理模块，显示位移计算结果，提取左边支座节点的转角和跨中挠度/POST1NSEL,S,1,5,4 ！选中1号和5号节点PRNSOL,U,Z ！列表显示选择集中所有节点的Z向位移PRNSOL,ROT,Y ！列表显示选择集中的所有节点的绕Y轴的转角ROTYNSEL,ALL ！重新选中所有节点 PRRSOL ！列表显示被约束的节点的支反力（Print Reaction Solution）RGHT_END = NODE (2.54,0,0) ！选择距离(2.54,0,0)最近的节点，并存入RGHT_ENDLFT_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的转角ROTYFINISH在时间历程后处理器中，叠加处理5号节点的弯矩/POST26 RFORCE,2,RGHT_END,M,Y ！叠加RGHT_END节点的弯矩MY，存入第2个POST26变量STORE ！保存变量*GET,M_MX_C2,VARI,2,EXTREM,VMAX ！获取第2个POST26变量的最大值FINISHEX2.4 （SOLID5，SOLID95，SOLID98，SOLID147）smrt,off ! 关闭智能化功能LESIZE,ALL,10 ! 定义所有线的划分段数为10段ESIZE,25.4 ! 定义默认的单元尺寸为25.4V,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=1LCDEF,ARG1,ARG1 ! 基于计算结果定义荷载工况NSEL,S,LOC,X,152.4 ! 选择节点PRNSOL,U,COMP ! 打印位移结果及其节点位移矢量和LCFACT,ARG1,ARG2 ! 将缩放因子“AGR2”应用给荷载工况1LCASE,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_1FINISHNSEL,S,LOC,X,152.4NSEL,R,LOC,Y,2.54NSEL,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自由度为耦合自由度，定义为集合1CP,2,ROTY,2,12 ! 定义 2号节点和12号节点的ROTY自由度为耦合自由度，定义为集合1CPSGEN,6,1,1,2 ! 按照前面定义的耦合自由度集合1到2，循环6次生成其余的耦合自由度ETABLE,STRS,S,1 ！用壳单元SHELL63的最大拉应力S1为内容定义单元表STRSESORT,STRS ！按最大拉应力对单元表STRS的排序*GET,SMAX1,SORT,MAX ！获取排序后的最大拉应力值，存入变量SMAX1PRNSOL,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*ENDDORMODIF,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-4NCNVTOL,M,1E-1 ! 设置集中力偶的收敛误差为0.1N.MEND_NODE = NODE (0,3048/2,0) ! 选择距离(0,3048/2,0)最近的节点，命名为END_NODE*GET,DEF,NODE,END_NODE,U,X ！提取节点END_NODE的位移UX，存入变量DEFETABLE,TENS,NMISC,1 ！以BEAM54的最大拉应力（NMISC,1）为内容定义单元表TENSETABLE,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_COMPEX3.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荷载步的计算结果PLDISPEX3.4 （BEAM188，BEAM189）SECNUM,1 ！设置当前截面编号为1SECTYPE,1,BEAM,ASEC ！定义梁的截面形状为任意形状，截面几何数据由用户提供SECDATA,11613,5.61868e8,0.0,224747,0.0,898989.1 ！截面几何参数以上定义的截面几何参数实际上是实常数ET,1,BEAM188MP,EX,1,491.2105MP,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 ！将中间计算过程输出到草稿文件SCRATCHSOLVE/OUT ！结束数据重定向，关闭输出文件设置求解选项,施加位移约束后求解面内力作用发生平面弯曲的情况OUTRES,ALL,ALL ARCLEN,ON ！激活弧长法，默认状态时不激活ARCTRM,U,1524,41,UZ ！设置弧长法的停止条件：41号节点的UZ向位移超过1524mmF,41,FX,6.6057 NSUBST,10 ！指定荷载子步为10NSOL,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子步的计算结果，缩放因子为1NSEL,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:100EX4.4 （PLANE42）OUTPR,NSOL,1 ！输出第1荷载步的节点计算结果OUTPR,ESOL,1 ！输出第1荷载步的单元计算结果EX4.6 （PLANE42）CSYS,1 ！切换到柱坐标系K,1,89 ！在坐标（89,0,0）处定义关键点1K,2,89,90 ！在坐标（89,90,0）处定义关键点2KGEN,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 ！提取排序后的最大环向应力，存入变量SIEX4.7 （PLANE42）NSEL,S,LOC,Y,90 ！选择节点CP,1,UY,ALL ！耦合所有选定的节点Y向位移，定义为集合1EX4.8 （PLANE42）LESIZE,3,4,.25 ！给3号线定义划分段数为4，间距比为0.25LARC,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.895MPaNSORT,S,X,3 ！将计算结果按SX排序，只保留前3项最高结果PRNSOL,S,COMP ！列表下式应力分量计算结果/WINDOW,1,OFF ！关闭1号窗口的显示/NOERASE ！不擦除图形窗口/DSCALE,2,1 ！用缩放位移约束命令，缩放位移约束实部2倍，虚部1倍/WINDOW,2,RTOP ！在图形窗口右上角定义窗口2PLNSOL,S,X ！图形显示水平方向的节点应力SX分布图LSEL,S,LINE,1,2 ！选择1、2号线NSLL,1 ！选择1号线上的节点NWRITE ！将选择的节点写入数据文件Jobanme.NODEEX4.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,13EX8.1 （BEAM3）在时间历程后处理器获取频率计算结果/POST26*GET,FREQ1,MODE,1,FREQ ! 获取第1阶频率的计算结果，并存入变量FREQ1*GET,FREQ2,MODE,2,FREQ*GET,FREQ3,MODE,3,FREQEX8.2 （BEAM3）在一般后处理模块中,获取频率计算结果/POST1 ！可以再/SOLU、/POST1、/POST26中提取频率*GET,FREQ1,MODE,1,FREQ ！获取第1阶频率的计算结果，并存入变量FREQ1*GET,FREQ2,MODE,2,FREQ*GET,FREQ3,MODE,3,FREQEX8.4 （SHELL63）OUTPR,ALL,1 ！指定输出选项为第1阶频率的所有计算输出参数OUTRES,ALL,0 ！指定所有输出参数不用保存到数据库EX9.3 （SHELL63）CSYS,1N,1,.0001N,10,76.2e-3FILLNGEN,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, DZEGEN,3,10,-1 ！循环3次，每次单元的节点号增量为10，按照前面1个单元的模式生成单元EGEN,9,1,-3 ！循环9次，每次单元的节点号增量为1，按照前3个单元的模式生成单元NSEL,S,LOC,Y,0 ！选择所有Y=0（度）的直线边上的所有节点CM,RIGHT,NODE ！将这些选中的节点组合成一个部件，叫做RIGHTEX9.9 （PLANE42，SOLID45）ET,1,Plane42ET,2,Solid45Esize,6.35 ！指定单元划分尺寸为6.35AMESH,ALL ！将所有面划分为PLANE42单元指定拉伸选项,通过体拉伸生成体单元TYPE, 2 ! 切换单元类型为第2类单元SOLID45EXTOPT,ESIZE,10,0, ！拉伸生成选项，生成10个体单元EXTOPT,ACLEAR,0 ！拉伸时不清除平面单元的网格EXTOPT,ATTR,0,0,0 ！保留体单元的特征类型不变体拉伸命令格式：VEXT, NA1, NA2, NINC, DX, DY, DZ, RX, RY, RZVEXT,1, , ,0,0,524, ！按照1号面，执行体拉伸，Z向增量了DZ=254EX10.1 （PLANE82）设置单元剖分的尺寸并进行单元剖分AESIZE,ALL,5 ！划分单元尺寸设置为5mmAMESH,ALL ！对所有面进行单元划分EX11.1 （PLANE82）/POST1 ！进入一般后处理模块SET,LAST ！读入最后一个荷载步的计算结果ETABLE,VOLU,VOLU ！用单元体积定义单元表VOLUPRNSOL,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 ！将汇总后的总体积保存到变量TVOLTVOL=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) ！将端部的最大挠度值取绝对值后，存入DEFLL,1,6 ！过关键点1和6定义直线*REPEAT,5,1,1 ！重复执行前面定义直线的命令，每次节点编号增加1*USE,SCRATCH ！指定宏程序SCRATCH/OPT ！进入优化设计模块OPANL,SCRATCH ！指定优化循环文件SCRATCHOPVAR,TVOL,OBJ,0.25 ！定义目标函数为总体积TVOL，允许容差为0.25单位OPVAR,STRS,SV,207 ！定义应力状态变量DEFL从0到207OPVAR,TK16,DV,3.81,6.86,0.0254 ！定义设计变量TK16从3.81到6.86，容差0.0254OPSAVE,INITIAL,OPT ！保存初始设计方案到文件INITIAL.OPTOPTYPE,SUBP ! 使用子模型近似优化设计方法OPSUBP,30 ！最大优化设计30个循环用子模型近似优化设计方法开始优化设计的计算OPEXE ！用子模型近似优化设计方法开始优化设计的计算VR1=TVOL ！将体积TVOL保存到变量VR1PARSAV,RSET1 ！将所有参数变量保存到文件Jobname.RSET1OPLIST,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 ！为一阶方法删除保证锥度的约束DIF1OPTYPE,FIRST ！指定使用一阶方法OPFRST,20 ！最大循环次数为20用一阶方法执行优化设计OPEXE ！用一阶方法执行优化设计VR4=TVOL ！将体积TVOL保存到变量VR4OPLIST,ALL,1 ！显示设计集中的参数/TITLE,EX11.3, SHAPE OPTIMIZATION OF A CANTILEVER BEAM (1st Method)/AXLAB,Y,VOLUME (TVOL)PL