小弟最近在做焊接温度场的模拟,使用的软件是ANSYS,在加载双椭球热源的时候遇到了点问题,就是热源的前半球和后半球的位置总是相反的,应该是前半球的半径短点,可是我不知道哪里出了问题,总是出现这种情况。
( )
附上代码
/SOL
*DO,I,0,3
*SET,DISY,I*0.002
*SET,T,(DISY+0.002)/0.01
ANTYPE,4
TRNOPT,FULL
TUNIF,293,
OUTRES,BASIC,LAST
TIME,T
NROPT,FULL,,ON
AUTOTS,1
NSUBST,5,,,1
KBC,1
lnsrch,1
TSRES,ERASE
TIMINT,1,THERM
NSEL,S,EXT
NSEL,U,LOC,X,0
SF,ALL,CONV,-1,293
ALLSEL,ALL
LOCAL,11,2,0,DISY,0.009,,,,1,15 //前半球
NSEL,S,LOC,X,0,0.004
NSEL,R,LOC,Y,0,90
NSEL,R,LOC,Z,-90,0
CM,NODEF,NODE
CMSEL,,NODEF
NPLOT
ALLSEL,ALL
LOCAL,12,0,0,DISY,0.009,,,,1,1 !!!考虑焦点位置变化坐标
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNCCSYS
*SET,_FNCNAME,\\\\\\\'tb\\\\\\\'
*DEL,\\\\\\\'tb\\\\\\\'
*SET,_FNCCSYS,12
*DIM,%_FNCNAME%,TABLE,6,26,1,,,,%_FNCCSYS%
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 10, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 12, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 17, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 0.4176614, 0, 0, -3
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, -3
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -1, 2, -3
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 3, 0, 0, -4
*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 1, -4, 3, -1
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 2, 0, 0, 2
*SET,%_FNCNAME%(0,12,1), 0.0, -4, 0, 1, 2, 17, -1
*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,14,1), 0.0, -5, 0, 1, 3, 17, -1
*SET,%_FNCNAME%(0,15,1), 0.0, -1, 0, 1, -4, 1, -5
*SET,%_FNCNAME%(0,16,1), 0.0, -4, 0, 2, 0, 0, 4
*SET,%_FNCNAME%(0,17,1), 0.0, -5, 0, 1, 4, 17, -4
*SET,%_FNCNAME%(0,18,1), 0.0, -4, 0, 1, -1, 1, -5
*SET,%_FNCNAME%(0,19,1), 0.0, -1, 0, 1, -3, 3, -4
*SET,%_FNCNAME%(0,20,1), 0.0, -3, 0, 0.002, 0, 0, 0
*SET,%_FNCNAME%(0,21,1), 0.0, -4, 0, 2, 0, 0, -3
*SET,%_FNCNAME%(0,22,1), 0.0, -5, 0, 1, -3, 17, -4
*SET,%_FNCNAME%(0,23,1), 0.0, -3, 0, 1, -1, 4, -5
*SET,%_FNCNAME%(0,24,1), 0.0, -1, 7, 1, -3, 0, 0
*SET,%_FNCNAME%(0,25,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,26,1), 0.0, 99, 0, 1, -3, 0, 0
BF,NODEF,HGEN,%tb%
CMDELE,NODEF
ALLSEL,ALL
LOCAL,13,2,0,DISY,0.009,,,,6,15 //后半球
NSEL,S,LOC,X,0,0.004
NSEL,R,LOC,Y,270,360
NSEL,R,LOC,Z,0,90
CM,NODEB,NODE
CMSEL,,NODEB
NPLOT
ALLSEL,ALL
LOCAL,14,0,0,DISY,0.009,,,,1,1 !!!考虑焦点位置变化坐标
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNCCSYS
*SET,_FNCNAME,\\\\\\\'tbb\\\\\\\'
*DEL,\\\\\\\'tbb\\\\\\\'
*SET,_FNCCSYS,14
*DIM,%_FNCNAME%,TABLE,6,26,1,,,,%_FNCCSYS%
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 10, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 12, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 17, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 0.2361406, 0, 0, -3
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, -3
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -1, 2, -3
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 3, 0, 0, -4
*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 1, -4, 3, -1
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 2, 0, 0, 2
*SET,%_FNCNAME%(0,12,1), 0.0, -4, 0, 1, 2, 17, -1
*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,14,1), 0.0, -5, 0, 1, 3, 17, -1
*SET,%_FNCNAME%(0,15,1), 0.0, -1, 0, 1, -4, 1, -5
*SET,%_FNCNAME%(0,16,1), 0.0, -4, 0, 2, 0, 0, 4
*SET,%_FNCNAME%(0,17,1), 0.0, -5, 0, 1, 4, 17, -4
*SET,%_FNCNAME%(0,18,1), 0.0, -4, 0, 1, -1, 1, -5
*SET,%_FNCNAME%(0,19,1), 0.0, -1, 0, 1, -3, 3, -4
*SET,%_FNCNAME%(0,20,1), 0.0, -3, 0, 0.002, 0, 0, 0
*SET,%_FNCNAME%(0,21,1), 0.0, -4, 0, 2, 0, 0, -3
*SET,%_FNCNAME%(0,22,1), 0.0, -5, 0, 1, -3, 17, -4
*SET,%_FNCNAME%(0,23,1), 0.0, -3, 0, 1, -1, 4, -5
*SET,%_FNCNAME%(0,24,1), 0.0, -1, 7, 1, -3, 0, 0
*SET,%_FNCNAME%(0,25,1), 0.0, -3, 0, 1, -2, 3, -1
*SE T,%_FNCNAME%(0,26,1), 0.0, 99, 0, 1, -3, 0, 0
BF,NODEB,HGEN,%tbb%
CMDELE,NODEB
ALLSEL,ALL
SOLVE
BFDELE,ALL,HGEN
*ENDDO
SAVE |