ANSYS稳定渗流程序
——不需要迭代多次
/COM,Thermal
/title,3M2M
/filname,dam
/PREP7
H1=1.45
H2=0.75
V=0.1
TT=7
et,1,55
MP,KXX,1,9.59
MP,KXX,2,0.00959
k,1,0,0
k,2,3,0
K,3,3,1.45
k,4,3,2
k,5,0,2
K,6,0,1.45
K,7,0,0.75
a,1,2,3,4,5,6,7
asel,,,,1
aatt,1,,1
allsel,all
ESIZE,0.05,0,
amesh,all
*GET,number1,node,,NUM,MAX
*GET,number2,ELEM,,NUM,MAX
*DIM,NHEAD,,number1,5
/solu
antype,STATIC
time,0.001
nropt, full
nsubst,10,0,0
FLST,2,30,1,ORDE,3
FITEM,2,2
FITEM,2,62
FITEM,2,-90
D,P51X, ,1.45, , , ,TEMP, , , , ,
FLST,2,30,1,ORDE,4
FITEM,2,1
FITEM,2,173
FITEM,2,-200
FITEM,2,162
D,P51X, ,1.45, , , ,TEMP, , , , ,
*DO,i,1,number1
NHEAD(i,1)=i
NHEAD(i,2)=NX(i)
NHEAD(i,3)=NY(i)
NHEAD(i,4)=TEMP(i)
NHEAD(i,5)=TEMP(i)-NY(i)
*IF,NX(i),EQ,0.0,OR,NX(i),EQ,3.0,THEN
*IF,NY(i),GT,1.45,THEN
HEAD=MIN(TEMP(i),NY(i))
FLST,2,1,1,ORDE,1
FITEM,2,i
D,P51X, ,HEAD, , , ,TEMP, , , , ,
*ENDIF
*ENDIF
*IF,NHEAD(i,3),GT,NHEAD(i,4),THEN
NHEAD(i,4)=0
*ENDIF
*ENDDO
*DO,I,1,number2
ZSUM=0
HSUM=0
*DO,J,1,4
*GET,NN1,ELEM,I,NODE,J
ZSUM=ZSUM+NY(NN1)
HSUM=HSUM+TEMP(NN1)
*ENDDO
*IF,ZSUM,GT,HSUM,THEN
mpchg,2,I
*ENDIF
*ENDDO
solve
finish
*CREATE,ansuitmp
*CFOPEN,'3M2M','dat','F:\dam\'
*set,J,1
*VWRITE,NHEAD(J,2),NHEAD(J,3) ,NHEAD(J,5) ,NHEAD(J,4) , , , , , ,
(F10.4,F10.4,F10.4,F10.4)
*CFCLOS
*END
/INPUT,ansuitmp
