DoubleEllipsoid_ansysworkbench_workbench椭球_

FINISH /PREP7 !双椭球热源模型 Q=ARG1 V=ARG2 a=2e-3 b=5e-3 c1=3e-3 c2=7e-3 f1=2/3 f2=4/3 N=100 *DO,I,1,N BFDELE,ALL,ALL T=I*1E-3/V nsel,s,loc,y,I*1e-3,0.1 CM,FRONT,NODE nsel,s,loc,y,0,I*1e-3 CM,BEHIND,NODE !热源公式1：10.392*f1*Q/a/b/c1/5.5655*exp(-3*({X}^2/a^2+({Y}-{TIME}*V)^2/c1^2+{Z}^2/b^2)) !热源公式2：10.392*f2*Q/a/b/c2/5.5655*exp(-3*({X}^2/a^2+({Y}-{TIME}*V)^2/c2^2+{Z}^2/b^2)) *DEL,_FNCNAME *DEL,_FNCMTID *DEL,_FNC_C1 *DEL,_FNC_C2 *DEL,_FNC_C3 *DEL,_FNC_C4 *DEL,_FNC_C5 *DEL,_FNC_C6 *DEL,_FNCCSYS *SET,_FNCNAME,'FRONT' *DIM,_FNC_C1,,1 *DIM,_FNC_C2,,1 *DIM,_FNC_C3,,1 *DIM,_FNC_C4,,1 *DIM,_FNC_C5,,1 *DIM,_FNC_C6,,1 *SET,_FNC_C1(1),f1 *SET,_FNC_C2(1),Q *SET,_FNC_C3(1),a *SET,_FNC_C4(1),b *SET,_FNC_C5(1),c1 *SET,_FNC_C6(1),V *SET,_FNCCSYS,0 ! /INPUT,.\Desktop\FRONT.func,,,1 *DIM,%_FNCNAME%,TABLE,8,36,1,,,,%_FNCCSYS% ! ! Begin of equation: 10.392*f1*Q/a/b/c1/5.5655*exp(-3*({X}^2/a^2+({Y}-{TIME}*V) ! ^2/c1^2+{Z}^2/b^2)) *SET,%_FNCNAME%(0,0,1), 0.0, -999 *SET,%_FNCNAME%(2,0,1), 0.0 *SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)% *SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)% *SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)% *SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)% *SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)% *SET,%_FNCNAME%(8,0,1), %_FNC_C6(1)% *SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 10.392, 0, 0, 17 *SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, -1, 3, 17 *SET,%_FNCNAME%(0,3,1), 0, -1, 0, 1, -2, 3, 18 *SET,%_FNCNAME%(0,4,1), 0.0, -2, 0, 1, -1, 4, 19 *SET,%_FNCNAME%(0,5,1), 0.0, -1, 0, 1, -2, 4, 20 *SET,%_FNCNAME%(0,6,1), 0.0, -2, 0, 1, -1, 4, 21 *SET,%_FNCNAME%(0,7,1), 0.0, -1, 0, 5.5655, 0, 0, -2 *SET,%_FNCNAME%(0,8,1), 0.0, -3, 0, 1, -2, 4, -1 *SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 0, 0, 0, 0 *SET,%_FNCNAME%(0,10,1), 0.0, -2, 0, 1, 0, 0, -1 *SET,%_FNCNAME%(0,11,1), 0.0, -4, 0, 1, -1, 2, -2 *SET,%_FNCNAME%(0,12,1), 0.0, -1, 0, 3, 0, 0, -4 *SET,%_FNCNAME%(0,13,1), 0.0, -2, 0, 1, -4, 3, -1 *SET,%_FNCNAME%(0,14,1), 0.0, -1, 0, 2, 0, 0, 2 *SET,%_FNCNAME%(0,15,1), 0.0, -4, 0, 1, 2, 17, -1 *SET,%_FNCNAME%(0,16,1), 0.0, -1, 0, 2, 0, 0, 19 *SET,%_FNCNAME%(0,17,1), 0.0, -5, 0, 1, 19, 17, -1 *SET,%_FNCNAME%(0,18,1), 0.0, -1, 0, 1, -4, 4, -5 *SET,%_FNCNAME%(0,19,1), 0.0, -4, 0, 1, 1, 3, 22 *SET,%_FNCNAME%(0,20,1), 0.0, -5, 0, 1, 3, 2, -4 *SET,%_FNCNAME%(0,21,1), 0.0, -4, 0, 2, 0, 0, -5 *SET,%_FNCNAME%(0,22,1), 0.0, -6, 0, 1, -5, 17, -4 *SET,%_FNCNAME%(0,23,1), 0.0, -4, 0, 2, 0, 0, 21 *SET,%_FNCNAME%(0,24,1), 0.0, -5, 0, 1, 21, 17, -4 *SET,%_FNCNAME%(0,25,1), 0.0, -4, 0, 1, -6, 4, -5 *SET,%_FNCNAME%(0,26,1), 0.0, -5, 0, 1, -1, 1, -4 *SET,%_FNCNAME%(0,27,1), 0.0, -1, 0, 2, 0, 0, 4 *SET,%_FNCNAME%(0,28,1), 0.0, -4, 0, 1, 4, 17, -1 *SET,%_FNCNAME%(0,29,1), 0.0, -1, 0, 2, 0, 0, 20 *SET,%_FNCNAME%(0,30,1), 0.0, -6, 0, 1, 20, 17, -1 *SET,%_FNCNAME%(0,31,1), 0.0, -1, 0, 1, -4, 4, -6 *SET,%_FNCNAME%(0,32,1), 0.0, -4, 0, 1, -5, 1, -1 *SET,%_FNCNAME%(0,33,1), 0.0, -1, 0, 1, -2, 3, -4 *SET,%_FNCNAME%(0,34,1), 0.0, -1, 7, 1, -1, 0, 0 *SET,%_FNCNAME%(0,35,1), 0.0, -2, 0, 1, -3, 3, -1 *SET,%_FNCNAME%(0,36,1), 0.0, 99, 0, 1, -2, 0, 0 ! End of equation: 10.392*f1*Q/a/b/c1/5.5655*exp(-3*({X}^2/a^2+({Y}-{TIME}*V)^2 ! /c1^2+{Z}^2/b^2)) !--> *DEL,_FNCNAME *DEL,_FNCMTID *DEL,_FNC_C1 *DEL,_FNC_C2 *DEL,_FNC_C3 *DEL,_FNC_C4 *DEL,_FNC_C5 *DEL,_FNC_C6 *DEL,_FNCCSYS *SET,_FNCNAME,'BEHIND' *DIM,_FNC_C1,,1 *DIM,_FNC_C2,,1 *DIM,_FNC_C3,,1 *DIM,_FNC_C4,,1 *DIM,_FNC_C5,,1 *DIM,_FNC_C6,,1 *SET,_FNC_C1(1),f2 *SET,_FNC_C2(1),Q *SET,_FNC_C3(1),a *SET,_FNC_C4(1),b *SET,_FNC_C5(1),c2 *SET,_FNC_C6(1),V *SET,_FNCCSYS,0 ! /INPUT,.\Desktop\BEHIND.func,,,1 *DIM,%_FNCNAME%,TABLE,8,36,1,,,,%_FNCCSYS% ! ! Begin of equation: 10.392*f2*Q/a/b/c2/5.5655*exp(-3*({X}^2/a^2+({Y}-{TIME}*V) ! ^2/c2^2+{Z}^2/b^2)) *SET,%_FNCNAME%(0,0,1), 0.0, -999 *SET,%_FNCNAME%(2,0,1), 0.0 *SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)% *SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)% *SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)% *SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)% *SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)% *SET,%_FNCNAME%(8,0,1), %_FNC_C6(1)% *SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 10.392, 0, 0, 17 *SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, -1, 3, 17 *SET,%_FNCNAME%(0,3,1), 0, -1, 0, 1, -2, 3, 18 *SET,%_FNCNAME%(0,4,1), 0.0, -2, 0, 1, -1, 4, 19 *SET,%_FNCNAME%(0,5,1), 0.0, -1, 0, 1, -2, 4, 20 *SET,%_FNCNAME%(0,6,1), 0.0, -2, 0, 1, -1, 4, 21 *SET,%_FNCNAME%(0,7,1), 0.0, -1, 0, 5.5655, 0, 0, -2 *SET,%_FNCNAME%(0,8,1), 0.0, -3, 0, 1, -2, 4, -1 *SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 0, 0, 0, 0 *SET,%_FNCNAME%(0,10,1), 0.0, -2, 0, 1, 0, 0, -1 *SET,%_FNCNAME%(0,11,1), 0.0, -4, 0, 1, -1, 2, -2 *SET,%_FNCNAME%(0,12,1), 0.0, -1, 0, 3, 0, 0, -4 *SET,%_FNCNAME%(0,13,1), 0.0, -2, 0, 1, -4, 3, -1 *SET,%_FNCNAME%(0,14,1), 0.0, -1, 0, 2, 0, 0, 2 *SET,%_FNCNAME%(0,15,1), 0.0, -4, 0, 1, 2, 17, -1 *SET,%_FNCNAME%(0,16,1), 0.0, -1, 0, 2, 0, 0, 19 *SET,%_FNCNAME%(0,17,1), 0.0, -5, 0, 1, 19, 17, -1 *SET,%_FNCNAME%(0,18,1), 0.0, -1, 0, 1, -4, 4, -5 *SET,%_FNCNAME%(0,19,1), 0.0, -4, 0, 1, 1, 3, 22 *SET,%_FNCNAME%(0,20,1), 0.0, -5, 0, 1, 3, 2, -4 *SET,%_FNCNAME%(0,21,1), 0.0, -4, 0, 2, 0, 0, -5 *SET,%_FNCNAME%(0,22,1), 0.0, -6, 0, 1, -5, 17, -4 *SET,%_FNCNAME%(0,23,1), 0.0, -4, 0, 2, 0, 0, 21 *SET,%_FNCNAME%(0,24,1), 0.0, -5, 0, 1, 21, 17, -4 *SET,%_FNCNAME%(0,25,1), 0.0, -4, 0, 1, -6, 4, -5 *SET,%_FNCNAME%(0,26,1), 0.0, -5, 0, 1, -1, 1, -4 *SET,%_FNCNAME%(0,27,1), 0.0, -1, 0, 2, 0, 0, 4 *SET,%_FNCNAME%(0,28,1), 0.0, -4, 0, 1, 4, 17, -1 *SET,%_FNCNAME%(0,29,1), 0.0, -1, 0, 2, 0, 0, 20 *SET,%_FNCNAME%(0,30,1), 0.0, -6, 0, 1, 20, 17, -1 *SET,%_FNCNAME%(0,31,1), 0.0, -1, 0, 1, -4, 4, -6 *SET,%_FNCNAME%(0,32,1), 0.0, -4, 0, 1, -5, 1, -1 *SET,%_FNCNAME%(0,33,1), 0.0, -1, 0, 1, -2, 3, -4 *SET,%_FNCNAME%(0,34,1), 0.0, -1, 7, 1, -1, 0, 0 *SET,%_FNCNAME%(0,35,1), 0.0, -2, 0, 1, -3, 3, -1 *SET,%_FNCNAME%(0,36,1), 0.0, 99, 0, 1, -2, 0, 0 ! End of equation: 10.392*f2*Q/a/b/c2/5.5655*exp(-3*({X}^2/a^2+({Y}-{TIME}*V)^2 ! /c2^2+{Z}^2/b^2)) !--> allsel,all /sol ANTYPE,4 TRNOPT,FULL TUNIF,20 BF,FRONT,HGEN,%FRONT% BF,BEHIND,HGEN,%BEHIND% CMDELE,FRONT CMDELE,BEHIND AUTOTS,1 KBC,1 TIME,T allsel,all OUTRES,ALL,ALL SOLVE *ENDDO save