!***********建立热源数组*****************
QMAX=P/(3.1415*R**2)
MAX_TIME=len/LSIZE+1
MAX_X=wid/LSIZE+1
MAX_Y=len/ LSIZE+1
!**********施加对流载荷*************
ESEL,S,LIVE
ASEL,S,AREA,,9,12,1
sfa,all,1,conv,150,20
ESEL,S,LIVE
ASEL,S,AREA,,7,7,1
sfa,all,1,conv,150,20
*do,h,1,road,1
m=mod(h,2)
*if,m,eq,1,then
!h=1
!热源数组
*DIM,FLUX1,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME
*DO,K,1,MAX_TIME,1 !循环时间
*DO,I,1,MAX_X,1
FLUX1(I,0,K)=(I-1)*LSIZE
*ENDDO
*DO,J,1,MAX_Y,1
FLUX1(0,J,K)=(J-1)*LSIZE
*ENDDO
*DO,I,1,MAX_X,1
*DO,J,1,MAX_Y,1
XCENTER=(h-1)*weld+weld/2
YCENTER=V*(K-1)*TINC
DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-YCENTER)**2)
*IF,DISTANCE,LE,R,THEN
FLUX1(I,J,K)=2*A*QMAX*EXP(-2*DISTANCE**2/R**2)
*ELSE
FLUX1(I,J,K)=0
*ENDIF
*ENDDO
*ENDDO
FLUX1(0,0,K)=(K-1)*TINC
!激活单元
esel,s,live
nsel,s,loc,x,XCENTER-weld/2 ,XCENTER+weld/2
nsel,r,loc,z,0,hth
nsel,r,loc,y,YCENTER,YCENTER+LSIZE
esln,s,0
ealive,all
allsel,all
!施加热流载荷
ESEL,S,LIVE
NSEL,S,LOC,Z,hth,hth!选择上边界
SF,all,HFLUX,%FLUX1% !热流加载
ALLSEL,ALL
time,(h-1)*(TI+tinc)+(k)*TINC
NSUBST,2
SOLVE
NSEL,S,LOC,Z,hth,hth!选择上边界
SFDELE,ALL,HFLUX !删除上一步加载的热源
/REPLOT
*ENDDO
*else
*DIM,FLUX2,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME
*DO,K,1,MAX_TIME,1 !循环时间
*DO,I,1,MAX_X,1
FLUX2(I,0,K)=(I-1)*LSIZE
*ENDDO
*DO,J,1,MAX_Y,1
FLUX2(0,J,K)=(J-1)*LSIZE
*ENDDO
*DO,I,1,MAX_X,1
*DO,J,1,MAX_Y,1
XCENTER=(h-1)*weld+weld/2
YCENTER=len-V*(K-1)*TINC
DISTANCE=SQRT(((I-1)*LSIZE-XCENTER)**2+((J-1)*LSIZE-YCENTER)**2)
*IF,DISTANCE,LE,R,THEN
FLUX2(I,J,K)=2*A*QMAX*EXP(-2*DISTANCE**2/R**2)
*ELSE
FLUX2(I,J,K)=0
*ENDIF
*ENDDO
*ENDDO
FLUX2(0,0,K)=(K-1)*TINC
!激活单元
esel,s,live
nsel,s,loc,x,XCENTER-weld/2 ,XCENTER+weld/2
nsel,r,loc,z,0,hth
nsel,r,loc,y,YCENTER,YCENTER+LSIZE
esln,s,0
ealive,all
allsel,all
!施加热流载荷
ESEL,S,LIVE
NSEL,S,LOC,Z,hth,hth!选择上边界
SF,all,HFLUX,%FLUX2% !热流加载
ALLSEL,ALL
time,(K)*TINC+(h-1)*(TI+tinc)
NSUBST,2
SOLVE
NSEL,S,LOC,Z,hth,hth!选择上边界
SFDELE,ALL,HFLUX !删除上一步加载的热源
/REPLOT
*ENDDO
*endif
*enddo
FINISH
这是求解的部分程序，希望大佬帮我看看哪错了？ |
回复此楼
