用ansys计算土石坝渗流的代码finish
/TITLE, EARTHDAM SEEPAGE
/PLOPTS,DATE,0
/PREP7
SMRT,OFF
ANTYPE,STATIC
ET,1,PLANE55
MP,KXX,1,1E-4
K,1,24,12
K,2,24,0
K,3,0,0
K,4,28,12
K,5,28,0
K,6,52,0
L,1,3
L,3,2
L,1,2
L,4,5
L,5,6
L,4,6
LESIZE,ALL,,,24
A,1,3,2
A,1,2,5,4
A,4,5,6
MSHK,2
MSHA,0,2D
AMESH,ALL
NUMMRG,NODE
*GET,N_MAX,NODE,,NUM,MAX
*GET,E_MAX,ELEM,,NUM,MAX
*DIM,N_TEMP,ARRAY,N_MAX
*DIM,N_PRE,ARRAY,N_MAX
LSEL,S,LINE,,1
NSLL,S,1
NSEL,R,LOC,Y,0,8
D,ALL,TEMP,8
LSEL,S,LINE,,5
NSLL,S,1
NSEL,R,LOC,X,42.9,52.1
*GET,Nc_NUM,NODE,,COUNT
*get,Nc_min,node,,num,min
*SET,DNN,Nc_min
*DO,I,1,Nc_NUM
D,DNN,TEMP,NY(DNN)
*if,I,LT,Nc_NUM,then
*SET,dnn,ndnext(dnn)
*endif
*ENDDO
ALLSEL,ALL
FINISH
/SOLU
SOLVE
FINISH
SAVE
*SET,MAXCOMP,20
*SET,DD_HEAT,0.01
*DO,COM_NUM,1,MAXCOMP
*SET,DD_H,0
FINI
/POST1
*DO,I,1,N_MAX
*IF,COM_NUM,NE,1,THEN
*SET,DD1,N_TEMP(I)
*IF,ABS(DD1-TEMP(I)),GT,DD_H,THEN
*SET,DD_H,ABS(DD1-TEMP(I))
*ENDIF
*ENDIF
*SET,N_TEMP(I),TEMP(I)
*SET,N_PRE(I),N_TEMP(I)-NY(I)
*ENDDO
*IF,COM_NUM,NE,1,and,DD_H,LE,DD_HEAT,exit
FINISH
/POST1
/CLABEL,,1
/EDGE,,0
/CONTOUR,,8,0,1,8
PLNSOL,TEMP
PLVECT,TG
PLVECT,TF, , , ,VECT,ELEM,ON,0
LSEL,S,LINE,,1
NSLL,S,1
PRRSOL,HEAT
FSUM,HEAT
*GET,Q_day,FSUM,0,ITEM,HEAT
*SET,Q_day,abs(Q_day)
ALLSEL,ALL
*DO,I,1,N_MAX
*SET,N_TEMP(I),TEMP(I)
*SET,N_PRE(I),N_TEMP(I)-NY(I)
DNSOL,I,TEMP,,N_PRE(I)
*ENDDO
PLNSOL,TEMP
FINI