ytdongwei的博客

用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