英语翻译/com obtain nodes' information*GET,NTOT,NODE,0,COUNT*GET

英语翻译
/com obtain nodes' information
*GET,NTOT,NODE,0,COUNT
*GET,Nmin,NODE,0,num,min
*GET,Nmax,NODE,0,num,max
*DIM,XY,ARRAY,NTOT,2
*dim,nv,array,ntot
*dim,nmold,array,ntot
*dim,nheat,array,ntot
*dim,ncenter,array,ntot
*dim,ncon1,array,ntot
*dim,ncon2,array,ntot
*dim,nnum,array,ntot
nheat0=0
nmold0=0
ncenter0=0
ncon10=0
ncon20=0
*DO,i,1,NTOT,1
NI=I+nmin-1
nnum(I)=NI
XY(I,1)=NX(nnum(I))
XY(I,2)=NY(nnum(I))
*GET,nbc,NODE,nnum(I),D,TEMP
*if,nbc,eq,100,then
nmold0=nmold0+1
nmold(nmold0)=nnum(I)
*elseif,nbc,eq,10,then
nheat0=nheat0+1
nheat(nheat0)=nnum(I)
*elseif,nbc,eq,2,then
ncon20=ncon20+1
ncon2(ncon20)=nnum(I)
*elseif,nbc,eq,1,then
ncon10=ncon10+1
ncon1(ncon10)=nnum(I)
*else
ncenter0=ncenter0+1
ncenter(ncenter0)=nnum(I)
*endif
*ENDDO
*MWRITE,XY(1,1),NCOOR,DAT,JIK,2,NTOT,
(2F10.6)
*MWRITE,nmold(1),nmold,DAT,JIK,1,nmold0,
(F10.0)
*MWRITE,ncenter(1),ncenter,DAT,JIK,1,ncenter0,
(F10.0)
*MWRITE,nheat(1),nheat,DAT,JIK,1,nheat0,
(F10.0)
*MWRITE,ncon1(1),ncon1,DAT,JIK,1,ncon10,
(F10.0)
*MWRITE,ncon2(1),ncon2,DAT,JIK,1,ncon20,
(F10.0)
*MWRITE,nnum(1),nnum,DAT,JIK,1,ntot,
(F10.0)
*GET,NELM,ELEM,0,COUNT
*GET,EL,ELEM,0,NUM,MIN
*DIM,ELEMS,ARRAY,NELM,4
*do,i,1,nelm,1
II=NELEM(I,1)-NMIN+1
JJ=NELEM(I,2)-NMIN+1
KK=NELEM(I,3)-NMIN+1
LL=NELEM(I,4)-NMIN+1
*CFOPEN,ELEMS,DAT
*VWRITE,II,JJ,KK,LL
(2X,F8.0,2X,F8.0,2X,F8.0,2X,F8.0)
*enddo
这是第一部分 有兴趣的再帮忙看看第二部分

好费劲啊!
你参考下!
关于屈曲分析的.
/solu !
antype,0 !
eqslv,spar !选择spar求解器
pstres,on !打开预应力选项
solve !
finish
/solu !再一次求解,前一次的主要是考虑预应力效应,这个屈曲分析中比较常见
antype,1 !屈曲分析
bucopt,lanb,1,0,0 !指定屈曲分析的分析选项,后面的三个数字是默认值
mxpand,1,0,0,1,0.001 !模态扩展分析扩展分析,并制定计算单元和反作用力求解信息,后面的0.001是默认值
solve
finish
/post1 !进入后处理器
set,first !选择第一载荷步
set,list !列出第一载荷步的求解信息结果
pldisp,1 !显示模态结果,也是位移,包括单元和轮廓线
finish
/config,nres,100 !允许存在结果文件中的最大号码（一般指的是载荷步）
/prep7 !进入预处理器
tb,biso,1,1,2 !指定属性数据表,用于同向双线型分析,塑性的分布情况等,指定1好材料属性,再1指定温度号码
tbtemp,0 !为数据表定义温度
tbdata,2.0e8,0,!为刚才定义的属性数据表定义数据,从默认开始,这个要结合所选择单元才能看明白,第一个省略说明是默认选项一次定义6个数据
upgeom,0.01,1,1,'qiangdu','rst'," !第一个是比例因子,向后面定义模型,小于1,说明缩小,其次两个1是载荷步,和子步；后面的是新生成的文件的名字和后缀将分析所得的位移加到有限元模型的 节点上并更新有限元模型的几何形状
finish
/solu
antype,0 !稳态分析
nlgeom,1 !指出包含大变形
outres,all,all !输出所有的结点和单元的求解信息
arclen,1,0 !激活非线性分析的弧长法,不过最大半径设为0
arctrm,u,0.5,78,uy !弧长求解终止条件,当一次位移等于或者大于设定值时停止分析,位移设定为0.5,以78结点的位移为标准,Uy方向有效,也就是78结点爱UY方向有0.5位移是就停止分析!
nsubst,200,1 !指定本载荷步为200,并把当前载荷步的最后时间步作为起始时间步
solve
finish
/post26 !进入时间后处理器
nsol,2,78,u,y,deflection !把78结点的在Y方向的位移结果储存在deflection标识名下,2为代号（特征号,这个求解的）
rforce,3,1861,f,x,reactionf !把1861结点在X的结构力储存在reactionf标识下,3为代号
/axlab,x,deflection !表明图形显示中X轴标签的显示内容,下面为Y轴的
/axlab,y,reactionf
xvar,2 !指定X轴的2代号的信息,图形,
plvar,3 !显示3代号的信息,这里也就是在Y轴的显示图形,最多可以显示10个变量,一次
finish !推出处理器