ANSYS单元特性之命令流算例

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