一个不很成熟的变参数批量分析的ANSYS程序

!为了批量研究变化参数得到的结果变化，进行批量分析（采用的是不断新建模型的方法）
!里面有不少参数还没有变为DO和ENDDO的方法
!还需继续改进
!华南理工大学 建筑学院 qjchen
!http://chenqj.blogspot.com
D=400
!b=100
!h=100
num1=10
num11=num1+1
num21=num11+num1
num31=num21+num1
num3=num1*3
num2=num1*2
cnum=4


*CFopen,myout,out
*DO,I,1,cnum,1
*DO,J,1,cnum,1
dh=I/2*D
db=J/2*D
/PREP7
ET,1,SHELL63
R,1,0.5
MP,EX,1,2E11
MP,NUXY,1,0.2
*DO,ROW,0,num1,1
*DO,COL,0,num3,1
*IF,COL,LT,num11,THEN
N,ROW*num31+COL+1,db*COL/num1-D*0.5-db,ROW*dh/num1,0
*ELSEIF,COL,LT,num21
N,ROW*num31+COL+1,-D*0.5+D*(COL-num1)/num1,ROW*dh/num1,0
*ELSE
N,ROW*num31+COL+1,db*(COL-num2)/num1+D*0.5,ROW*dh/num1,0
*ENDIF
*ENDDO
*ENDDO
*DO,ROW,0,9,1
*DO,NNUM,31*ROW+1,31*ROW+30,1
E,NNUM,NNUM+1,32+NNUM,31+NNUM
*ENDDO
*ENDDO
NAME=STRCAT(STRCAT('H',CHRVAL(I/10)),STRCAT('B',CHRVAL(J/10)))
!SAVE,NAME,DB
FINISH
/SOLU
D,11,,,,,,UX
D,11,,,,,,UY
D,21,,,,,,UY
D,12,,,,,,UY
D,13,,,,,,UY
D,14,,,,,,UY
D,15,,,,,,UY
D,16,,,,,,UY
D,17,,,,,,UY
D,18,,,,,,UY
D,19,,,,,,UY
D,20,,,,,,UY
F,321,FY,-100
F,322,FY,-100
F,323,FY,-100
F,324,FY,-100
F,325,FY,-100
F,326,FY,-100
F,327,FY,-100
F,328,FY,-100
F,329,FY,-100
F,330,FY,-100
F,331,FY,-100
NAME1=STRCAT(NAME,'_FORCE')
!SAVE,NAME1,DB
!/output,dd,myout,,
SOLVE
/post1
set,first
*do,jjj,156,186,1
nodenow=jjj
*GET,sz,NODE,jjj,S,Z
*GET,sx,NODE,jjj,S,X
*GET,sy,NODE,jjj,S,Y
*vwrite,('b:'),db,(' h:'),dh,(' i:'),jjj,(' sx:'),sx,(' sy:'),sy,(' sz:'),sz
(A3,F6.0,A3,F6.0,A3,F6.0,A4,F10.2,A4,F10.2,A4,F10.2)
*enddo
NAME1=STRCAT(NAME,'_RESULT')
!SAVE,NAME1,DB

FINISH
/PREP7
EDELE,ALL
NDELE,ALL
FINISH
*ENDDO
*ENDDO
*cfclos