| 查看: 833 | 回复: 1 | |||
[交流]
激光熔覆温度场模拟
|
| 各位，最近做ansys激光熔覆模拟时，发现做出来的温度忽高忽低，比如说第一载荷步后节点温度1600，第二个就降到600多，第三个又升到1600多，第四个又降回600多，一直反复，也不知道什么情况？有哪位遇到过吗？ |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
-
国自然面上复盘~欢迎讨论
已经有15人回复
-
中青基了要发朋友圈吗？
已经有3人回复
|
!***********建立热源数组***************** 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 这是求解的部分程序，希望大佬帮我看看哪错了？ |
2楼2021-06-21 20:56:16
10