NPLOT ! 只显示节点位置,不显示节点号码
NPLOT,1 !显示节点位置,显示节点号码
NLIST !列出节点在直角坐标系下的坐标值
PRDISP !列表显示节点位移值计算结果
PRDISP !列表显示节点位移值计算结果
PRETAB ! 显示单元表中的计算结果
PLETAB, MFORX ! 用色度图显示单元表MFORX中杆件轴力图
5 复制和填充
FILL,1,11 !在节点1和节点11之间均匀填充节点(2到10)
EGEN,10,1,1,1,1 !元素复制命令,复制10次,相对应节点增量为1,后三位数表示要复制的元素
EGEN,10,1,1 ! 按照前面单元的模式循环10次,生成其余的9个单元 NGEN,2,11,1,11,,12.7 ! 循环2次, 节点号增量11, 按1号节点到11号节点的范围, X坐标增加12.7mm。生成12号到22号节点
NGEN命令的格式是“NGEN, ITIME,INC, NODE1, NODE2, NINC, DX, DY, DZ, SPACE”表示循环ITIME次(包括原来的那一批节点),每次节点号增加INC
个,按照从NODE1到NODE2,增量为NINC的节点范围,每次循环时3个坐标增量为DX,DY和DZ,“SPACE”是间隔比例因子
EGEN, ITIME, NINC, IEL1, IEL2, IEINC, MINC,TINC, RINC, CINC, SINC, DX, DY, DZ,各项参数的意义是按照增量NINC循环ITIME次,按照从单元IEL1到IEL2,增量为IEINC范围内的单元
EGEN,3,10, -1 !循环3次,每次单元的节点号增量为10,按照前1个单元的模式生成
单元,这里“-1”表示前面刚定义的前1个单元
AGEN,ITIME,NA1,NA2,NINC,DX,DY,DZ,KINC,NOELEM,IMOVE面积复制命令 AGEN,2,1,2,1,,62 !循环两次,从1号面到2号面,每次循环面号增加1,向Y方向平移62mm
6 提取指定位置的节点和单元,定义单元表并获取弯曲应力
MID_NODE = NODE (2,,, ) ! 选择距离(2,0,0)位置最近的节点,命名为MID_NODE
*GET,DISP,NODE,MID_NODE,U,Y ! 提取名称为MID_NODE的节点的竖向位移UY,存入变量DISP
MID_ELM = ENEARN (MID_NODE) ! 选择距离MID_NODE节点最近的单元, 命名为MID_ELM
ETABLE,STRS,LS,3 ! 以LS,3(-Y一侧的弯曲应力SBYB)为内容定义单元表STRS
*GET,STRSS,ELEM,MID_ELM,ETAB,STRS ! 提取单元表STRS中单元MID_ELM的应力,存入 STRSS
7 弯矩和剪力单元表的提取
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 ! 结构弯矩分布图
8 移去前面施加的(横向)荷载和删除约束:
FDEL,11,FX
DDELE,ALL !删除所有节点的所有约束
9 显示节点应力:
/DEVICE,VECTOR,1 ! 切换显示风格为矢量线方式,便于单元网格的显示。 /Device,VECTOR,0 ! 关闭图形矢量显示开关 PLNSOL,S,X ! 在图形窗口显示节点上的X方向正应力 PLNSOL,S,XY ! 在图形窗口显示节点上的XY方向的剪应力 PLNSOL,S,Y ! 在图形窗口显示节点上的Y方向正应力
PLNSOL,S,1 ! 在图形窗口显示节点上的第一主应力(拉应力) PLNSOL,S,3 ! 在图形窗口显示节点上的第三主应力(压应力) PLNSOL,S,EQV ! 在图形窗口显示所有节点Mises应力
PRNSOL,S,PRIN !列表显示节点主应力(包括第一,二,三主应力和剪应力,切应力)
ETABLE,SIGY,S,Y ! 以Y方向的正应力为内容,定义单元表SIGY PRETAB,SIGY ! 显示单元表中的单元应力计算结果
PRNSOL,S,COMP ! 显示节点应力计算结果
PRNSOL,U,COMP ! 列表显示结点位移值
PRNSOL,U,Y !列表显示节点位移UY值
PRNSOL,DOF !列表显示节点位移分量结果
FSUM ! 统计汇总节点上的合力(Nsel,S,LOC,Z,0 ! 通过位置选择Z=0位置的节点)
NSORT,S,1 !对单元表按第一主应力排序 NSORT,U,Y,,1 对节点计算结果UY 排序
10 显示内存变量: *status,parm ! 显示内存变量的值
11 循环生成节点:
*DO,i,1,20,1 ! 水平方向节点的循环 *DO,j,1,10,1 ! 高度方向节点循环
n,i+(j-1)*20,(i-1)*5/19,(j-1)*1/9 ! 依次计算节点号,节点水平位置和铅直位置
*ENDDO ! 结束循环 *ENDDO ! 结束循环
*DO,i,1,19,1 ! 水平方向单元循环 *DO,j,1,9,1 ! 高度方向单元循环
E,i+(j-1)*20,i+(j-1)*20+1,i+(j-0)*20+1,i+(j-0)*20 ! 定义单元的四节点 *ENDDO ! 结束循环
*ENDDO ! 结束循环
12 圆弧命令
LARC,4,5,6,12.7 ! 以6号点为圆心,过4号点和5号点,以12.7mm为半径定义圆弧(LARC,P1,P2,PC,RAD)
13 节点的选择
NSEL,S,LOC,X,0 !选择X=0位置处的节点
NSEL,S,LOC,X,152.4 !选择位于X=152.4面上的所有节点
NSEL,R,LOC,Y,2.54 ! 重新(继续)选择Y=2.54mm的所有节点
LSEL,S,LINE,,1 ! 选择1号线
NSLL,,1 !选择1号线上的节点
NSEL,S,LOC,X,0,228.6 !选择X=0到X=228.6mm范围内的节点
14 映射网格划分
MSHK,2 !如果可能,采用映射网格剖分
MSHA,0,2D !使用四边形网格
MSHKEY,1 !设置单元网格划分类型为映射网格