当前位置:文档之家› ANSYS单元特性之命令流算例

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,V ALUE,,2,3!定义2*3的数值型数组V ALUELABEL(1) = 'STRS_MPa','DEF_mm' ! 给字符型数组的第1个元素赋值*VFILL,V ALUE(1,1),DATA,1,-0.05498 ! 给其他数值型数组中的元素赋值*VFILL,V ALUE(1,2),DATA,STRSS,DISP*VFILL,V ALUE(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),V ALUE(1,1),V ALUE(1,2),V ALUE(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类实常数,杆件截面面积为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!在所有梁单元表面施加均布荷载进入后处理模块,显示位移计算结果,提取左边支座节点的转角和跨中挠度/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号节点的弯矩/POST26RFORCE,2,RGHT_END,M,Y !叠加RGHT_END节点的弯矩MY,存入第2个POST26变量STORE!保存变量*GET,M_MX_C2,V ARI,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,V ALUE,,3,3 !LABEL(1,1) = 'DEFL X ','DEFL Y ','DEFL Z '!LABEL(1,2) = 'mm','mm','mm' !*VFILL,V ALUE(1,1),DATA,7.6E-4,2.471,10.965 !*VFILL,V ALUE(1,2),DATA,UX1,UY1,UZ1!*VFILL,V ALUE(1,3),DATA,ABS(UX1/(7.6E-4)),ABS(UY1/2.471),ABS(UZ1/10.965) SA VE, 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自由度为耦合自由度,定义为集合1 CPSGEN,6,1,1,2! 按照前面定义的耦合自由度集合1到2,循环6次生成其余的耦合自由度ETABLE,STRS,S,1 !用壳单元SHELL63的最大拉应力S1为内容定义单元表STRS ESORT,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)为内容定义单元表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_COMPEX3.1 (BEAM3)*GET,FCR,MODE,1,FREQ!提取第1阶频率计算结果,并保存到变量FCR.。

相关主题