恳请各位版主看看，这个热电分析为什么不收敛？
以前发过，当时没有得到解决方案。现在在很多地方完善了程序，观察不收敛的结果，更加符合理论值。问题是死活没有办法收敛。还请各位能够看看。不胜感激。弧焊电弧模拟，如果有同样经验的人看看条件更好。
材料只用了1和3，2是以后备用的，一起写上算了。
Finish
/clear
/PLOPTS,INFO,3
/PLOPTS,LOGO,0
/PLOPTS,WP,0
/PLOPTS,DATE,2
/TRIAD,OFF
/REPLOT
/RGB,INDEX,100,100,100, 0
/RGB,INDEX, 80, 80, 80,13
/RGB,INDEX, 60, 60, 60,14
/RGB,INDEX, 0, 0, 0,15
/TITLE,Arc_Welding_Simulation
/REPLOT
/PREP7
ET,1,SOLID69
MPTEMP,1,273,2500,5000,7500,10000,12500
MPTEMP,7,15000,17500,20000
MPDATA,DENS,1,1,1.67,0.2,0.1,0.07,0.04,0.03
MPDATA,DENS,1,7,0.02,0.015,0.014
MPDATA,KXX,1,1,0,0.1,0.2,0.3,0.6,1.6
MPDATA,KXX,1,7,2.4,2.0,2.2
MPDATA,RSVX,1,1,1000,100,1,1.3e-3,3.3e-4,1.6e-4
MPDATA,RSVX,1,7,1.2e-4,9.8e-5,8.5e-5
MPDATA,C,1,1,521,521,521,600,1250,6000
MPDATA,C,1,7,10000,4000,2500
MPDATA,VISC,1,1,2e-5,7e-5,1e-4,2.2e-4,2.7e-4,2e-4
MPDATA,VISC,1,7,6e-5,3e-5,2.6e-5
MPDATA,DENS,2,,7200
MPDATA,KXX,2,,20
MPDATA,RSVX,2,,1/770000
MPDATA,DENS,3,,19350
MPDATA,RSVX,3,,0.0000000528
MPDATA,KXX,3,,173
MPDATA,MURX,1,,1
MPDATA,MURX,3,,1
/REPLOT,RESIZE
TOFFST,0
CYL4,0,0,0.01,0,0.01,90,0.005+0.0014*sqrt(3)
CYL4,0,0,0.0002,0,0.0002,90,0.005
CYL4,0,0,0.0016,0,0.0016,90,0.005+0.0014*sqrt(3)
CONE,0.0002,0.0016,0.005,0.005+0.0014*sqrt(3),0,90,
VOVLAP,ALL
VGLUE,ALL
AGLUE,ALL
LGLUE,ALL
NUMCMP,ALL
/PNUM,KP,1
/PNUM,AREA,1
/VIEW,1,-1,-1,1
/VUP,1,Z
APLOT
VSEL,S,VOLU,,2
VATT,3,,1,0,
VSEL,S,VOLU,,1
VSEL,A,VOLU,,3,4
VATT,1,,1,0,
ALLSEL,ALL
ESIZE,0.0005,0,
MSHKEY,1
VMESH,ALL
FINISH
/SOL
TUNIF,10000,
DA,13,TEMP,1000,0
DA,1,TEMP,1000,0
ASEL,S,AREA,,9
ASEL,A,AREA,,2
DA,ALL,TEMP,3000,0
ASEL,S,AREA,,5
ASEL,A,AREA,,16
ASEL,A,AREA,,12
DA,ALL,TEMP,5000,0
DA,ALL,VOLT,0,0
ASEL,S,AREA,,8
NSLA,S,0
CP,next,VOLT,ALL
*GET,NMIN,NODE,,NUM,MIN
F,NMIN,AMPS,-160
ALLSEL,ALL
!cnvtol,amps,,,,-1
SOLVE
SAVE
FINISH 电流不输入吗？
页:
[1]