角接焊过程热分析，加载请教！！高手指点
我在做焊接温度场的分析，对于角接焊运用高斯移动热源加载时，我分别定义了两个数组，分别加载在两个面上，两面垂直相交（腹板和面板），可是显示的计算结果，只在一个面上加载成功了，另一个面没有加上！那位高手作过类似的热分析，帮忙指点一下！谢谢啦！下面是我加载部分的命令流，是我数组定义错误还是加载的部分有问题？帮忙指点一下！谢谢啦
/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 !! 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 这个问题很想知道，那位高手帮忙啊！
谢谢啦 解决了吗，呵呵
页:
[1]