©2010-2026 宁波自信网络信息技术有限公司 版权所有
客服电话：0574-28810668 投诉电话：18658249818
浙ICP备2021020529号-1 | 浙B2-20240490
关注我们 ：
1、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,
2、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
3、),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),VALU
4、E(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 ！只显示编号，不使用色彩 列表显示节点位
5、移和单元的计算结果 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
6、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
7、定义X方向的应力为单元表STRSX *GET,STRSSX,ELEM,3,ETAB,STRSX ！从单元表STRSX中提取3号单元的X向应力，存入STRSSX。注意：已知某一具体单元，在定义完单元表后，就可以直接用*GET语句提取这一具体单元的表中定义的项目。 EX2.1 （BEAM3） NPLOT ! 显示节点位置图形，但是不显示节点号码 NPLOT, 1 ！显示节点位置图形，同时显示节点号码 NLIST ！列表显示节点在直角坐标系下的坐标值 DSYS, 1 ！改变显示坐标系，列出资料时，转换到圆柱坐标系 NLIST ！列表显示节点在圆柱坐标系下的坐
8、标值 定义以两端弯矩和剪力为内容的单元表，并列出单元的单元表数据 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 ! 定义剪力图
9、窗口标题 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,RO
10、T,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 ！列表
11、显示选择集中所有节点的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
12、 *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 ! 关闭智能化功能
13、 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
14、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
15、 *** ! *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
16、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
17、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个表
18、中恢复数据，绘制单元应力云图 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 ! 以最佳比例显示
19、 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
20、自由度为耦合自由度，定义为集合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 ！
21、循环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 从数据库
22、中恢复数据 /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
23、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,F
24、CR,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,
25、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,116
26、13,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 ！结束数据重定向 …… /VI
27、EW,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 ！结束数据重定向，关闭输出文件 设置
28、求解选项,施加位移约束后求解面内力作用发生平面弯曲的情况 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
29、 ！将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号颜色的
30、三原色配比为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,8
31、9 ！选择节点 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,1
32、2.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号窗
33、口的显示 /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）
34、 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
35、 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,MOD
36、E,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坐标（这里是柱坐
37、标系的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（度）的直线边上的
38、所有节点 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,
39、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,V
40、OLU,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
41、选择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,SCRATC
42、H ！指定宏程序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 ! 使用子模型近似优化设计方法
43、 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 ！图形显示目标函数收敛
44、情况 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 CANTILEVE
45、R 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 !
46、单元轴向应力降幂排序 *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 ! 定
47、义设计变量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’ !
48、 读入文件 查看模态求解信息 在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学习资料，空间地址：
©2010-2026 宁波自信网络信息技术有限公司 版权所有
客服电话：0574-28810668 投诉电话：18658249818
浙ICP备2021020529号-1 | 浙B2-20240490
关注我们 ：
