角接焊过程热分析，加载请教！！高手指点

jy424

我在做焊接温度场的分析，对于角接焊运用高斯移动热源加载时，我分别定义了两个数组，分别加载在两个面上，两面垂直相交（腹板和面板），可是显示的计算结果，只在一个面上加载成功了，另一个面没有加上！
那位高手作过类似的热分析，帮忙指点一下！谢谢啦！下面是我加载部分的命令流，是我数组定义错误还是加载的部分有问题？帮忙指点一下！谢谢啦
/SOLU

V=0.0073
T=L/V
TINC=LSIZE1/V
U=28
I=270
Efficiency=0.85
KK=3
R1=0.015
R2=0.02
Qmax1=3*U*I*Efficiency/(3.1415926*R1**2)
Qmax2=3*U*I*Efficiency/(3.1415926*R2**2)


MAX_TIME=L/LSIZE1+1

MAX_Z=L/LSIZE1+1

MAX_X=L1/LSIZE1+1

*DIM,FLUX2,TABLE,MAX_X,MAX_Z,MAX_TIME,X,Z,TIME

*DO,K,1,MAX_TIME,1

*DO,J,1,MAX_Z,1

FLUX2(0,J,K)=(J-1)*LSIZE1

*ENDDO

*DO,I,1,MAX_X,1

FLUX2(I,0,K)=H-(I-1)*LSIZE1

*ENDDO

*DO,J,1,MAX_Z,1

*DO,I,1,MAX_X,1

ZCENTER=V*(K-1)*TINC

XCENTER=H

DISTANCE=SQRT(((J-1)*LSIZE1-ZCENTER)**2+((H-(I-1)*LSIZE1)-XCENTER)**2)

*IF,DISTANCE,LE,R1,THEN

FLUX2(I,J,K)=Qmax1/EXP(KK*DISTANCE**2/(R1**2))

*ELSE

FLUX2(I,J,K)=0

*ENDIF

*ENDDO

*ENDDO

FLUX2(0,0,K)=(K-1)*TINC

*ENDDO



MAX_TIME=L/LSIZE1+1
MAX_Z=L/LSIZE1+1
MAX_Y=L1/LSIZE1+1

*DIM,FLUX3,TABLE,MAX_Y,MAX_Z,MAX_TIME,Y,Z,TIME

*DO,K,1,MAX_TIME,1

*DO,J,1,MAX_Z,1

FLUX3(0,J,K)=(J-1)*LSIZE1

*ENDDO

*DO,I,1,MAX_Y,1

FLUX3(I,0,K)=B+(I-1)*LSIZE1

*ENDDO

*DO,J,1,MAX_Z,1

*DO,I,1,MAX_Y,1

ZCENTER=V*(K-1)*TINC

YCENTER=B

DISTANCE=SQRT(((J-1)*LSIZE1-ZCENTER)**2+(B+(I-1)*LSIZE1-YCENTER)**2)

*IF,DISTANCE,LE,R2,THEN

FLUX3(I,J,K)=Qmax2/EXP(KK*DISTANCE**2/(R2**2))

*ELSE

FLUX3(I,J,K)=0

*ENDIF

*ENDDO

*ENDDO

FLUX3(0,0,K)=(K-1)*TINC

*ENDDO

NSEL,,LOC,X,H
NSEL,R,LOC,Y,B,B
sf,all,hflux,%flux3%
NSEL,ALL

NSEL,,LOC,Y,B
NSEL,R,LOC,X,H-L1,H-LSIZE1
sf,all,hflux,%flux2%
NSEL,ALL

jinqinghai

!! Interesting,  follow you!

/SOLU

V=0.0073
T=L/V
TINC=LSIZE1/V
U=28
I=270
Efficiency=0.85
KK=3
R1=0.015
R2=0.02
Qmax1=3*U*I*Efficiency/(3.1415926*R1**2)
Qmax2=3*U*I*Efficiency/(3.1415926*R2**2)


MAX_TIME=L/LSIZE1+1

MAX_Z=L/LSIZE1+1

MAX_X=L1/LSIZE1+1

*DIM,FLUX2,TABLE,MAX_X,MAX_Z,MAX_TIME,X,Z,TIME

*DO,K,1,MAX_TIME,1

*DO,J,1,MAX_Z,1

FLUX2(0,J,K)=(J-1)*LSIZE1

*ENDDO

*DO,I,1,MAX_X,1

FLUX2(I,0,K)=H-(I-1)*LSIZE1

*ENDDO

*DO,J,1,MAX_Z,1

*DO,I,1,MAX_X,1

ZCENTER=V*(K-1)*TINC

XCENTER=H

DISTANCE=SQRT(((J-1)*LSIZE1-ZCENTER)**2+((H-(I-1)*LSIZE1)-XCENTER)**2)

*IF,DISTANCE,LE,R1,THEN

FLUX2(I,J,K)=Qmax1/EXP(KK*DISTANCE**2/(R1**2))

*ELSE

FLUX2(I,J,K)=0

*ENDIF

*ENDDO

*ENDDO

FLUX2(0,0,K)=(K-1)*TINC

*ENDDO



MAX_TIME=L/LSIZE1+1
MAX_Z=L/LSIZE1+1
MAX_Y=L1/LSIZE1+1

*DIM,FLUX3,TABLE,MAX_Y,MAX_Z,MAX_TIME,Y,Z,TIME

*DO,K,1,MAX_TIME,1

*DO,J,1,MAX_Z,1

FLUX3(0,J,K)=(J-1)*LSIZE1

*ENDDO

*DO,I,1,MAX_Y,1

FLUX3(I,0,K)=B+(I-1)*LSIZE1

*ENDDO

*DO,J,1,MAX_Z,1

*DO,I,1,MAX_Y,1

ZCENTER=V*(K-1)*TINC

YCENTER=B

DISTANCE=SQRT(((J-1)*LSIZE1-ZCENTER)**2+(B+(I-1)*LSIZE1-YCENTER)**2)

*IF,DISTANCE,LE,R2,THEN

FLUX3(I,J,K)=Qmax2/EXP(KK*DISTANCE**2/(R2**2))

*ELSE

FLUX3(I,J,K)=0

*ENDIF

*ENDDO

*ENDDO

FLUX3(0,0,K)=(K-1)*TINC

*ENDDO

NSEL,,LOC,X,H
NSEL,R,LOC,Y,B,B
sf,all,hflux,%flux3%
NSEL,ALL

NSEL,,LOC,Y,B
NSEL,R,LOC,X,H-L1,H-LSIZE1
sf,all,hflux,%flux2%
NSEL,ALL

jy424

这个问题很想知道，那位高手帮忙啊！
谢谢啦

sara_eva

解决了吗，呵呵