热点耦合分析出错！急求大家帮忙！

lilichunyan

原来是一个热分析，现在的是热电耦合的，不收敛，请大家帮忙看看原因吧！感激不尽！

1. 
   
2. Finish
   
3. /clear
   
4. /filname,redian,1
   
5. /UNITS,MKS
   
6. /PREP7
   
7. ET,1,solid69
   
8. ET,2,solid69
   
9. MPTEMP,1,273,2500,5000,7500,10000,12500
   
10. MPTEMP,7,15000,17500,20000  
    
11. MPDATA,RSVX,1,1,1000,100,1,1.3e-3,3.3e-4,1.6e-4
    
12. MPDATA,RSVX,1,7,1.2e-4,9.8e-5,8.5e-5
    
13. MPDATA,KXX,1,1,0,0.1,0.2,0.3,0.6,1.6
    
14. MPDATA,KXX,1,7,2.4,2.0,2.2
    
15. MPDATA,C,1,1,500,500,500,600,1250,6000
    
16. MPDATA,C,1,7,10000,4000,2500
    
17. MPDATA,DENS,1,1,1.1,0.2,0.1,0.07,0.04,0.03
    
18. MPDATA,DENS,1,7,0.02,0.015,0.014
    
19. !mp,KXX,2,,173   
    
20. !MP,RSVX,2,,0.0000000528
    
21. MPTEMP,,,,,,,,  
    
22. MPTEMP,1,0
    
23. MPDATA,KXX,2,,173   
    
24. MPTEMP,,,,,,,,  
    
25. MPTEMP,1,0  
    
26. MPDATA,RSVX,2,,0.0000000528
    
27. !建模
    
28. wpro,,-90   
    
29. CYL4,0,0,0.006,90, , ,0.006  
    
30. WPAVE,0,0.004,0  
    
31. CONE,0.00024,0.0014,0,0.002,0,90,
    
32. VSBW,all,,dele!分割平面
    
33. WPAVE,0,0,0  
    
34. CYL4,0,0,0.00024,90, , ,0.004  
    
35. CYL4,0,0,0.0014,90, , ,0.006   
    
36. CYL4,0,0,0.003,90, , ,0.006
    
37. VOVLAP,all
    
38. VGLUE,ALL
    
39. AGLUE,ALL
    
40. LGLUE,ALL
    
41. !NUMMRG,NODE
    
42. NUMCMP,all
    
43. /VIEW,1,1,2,3
    
44. WPCSYS,-1,0        !工作平面坐标归位
    
45. WPSTYLE,,,,,,,,0
    
46. 
    
47. !划分网格
    
48. VSEL,,,,2    !分配钨极材料属性
    
49. VATT,2,1,2
    
50. vsel,inve    !选择其他所有的，分配气体材料属性
    
51. VATT,1,1,1
    
52. allsel,all
    
53. AESIZE,7,0.0001,
    
54. AESIZE,12,0.0001,
    
55. AESIZE,18,0.0001,
    
56. AESIZE,22,0.0002,
    
57. AESIZE,26,0.0002,
    
58. AESIZE,32,0.0002,
    
59. AESIZE,34,0.0002,
    
60. LESIZE,2, , ,16, , , , ,1
    
61. MSHAPE,0,3d
    
62. MSHKEY,1
    
63. VMESH,all
    
64. finish
    
65. save
    
66. /sol
    
67. VSEL,S, , ,1
    
68. VSEL,A, , ,4
    
69. VSEL,A, , ,7
    
70. TUNIF,10000,
    
71. ALLSEL,ALL
    
72. asel,s,,,2
    
73. asel,a,,,9
    
74. DA,all,TEMP,3000
    
75. asel,s,,,5
    
76. asel,a,,,6
    
77. asel,a,,,25
    
78. asel,a,,,33  
    
79. DA,all,TEMP,1000
    
80. asel,s,loc,y,0
    
81. DA,all,TEMP,5000
    
82. DA,all,VOLT,0,0
    
83. ASEL,S,AREA,,4
    
84. NSLA,S,1
    
85. CP,next,VOLT,ALL
    
86. *GET,NMIN,NODE,,NUM,MIN                                 
    
87. F,NMIN,AMPS,-150               
    
88. ALLSEL,ALL  
    
89. solve
    
90. 
    

复制代码

[ 本帖最后由 lilichunyan 于 2009-6-6 20:47 编辑 ]

16443

ALLSEL,ALL
asel,s,,,2
asel,a,,,9
DA,all,TEMP,3000
asel,s,,,5
asel,a,,,6
asel,a,,,25
asel,a,,,33  
DA,all,TEMP,1000

asel,s,loc,y,0
DA,all,TEMP,5000

是否应该在每次asel,s,..之前增加一个allsel,all

jxxansys

"VSEL,S, , ,1
VSEL,A, , ,4
VSEL,A, , ,7
TUNIF,10000,  
ALLSEL,ALL "

1。tunif命令为整个模型施加初始温度，所以选择命令无效
2。你的命令流我执行并收敛了，附件为计算结果。
3。根据ANSYS HELP，SOLID70单元的材料特性为
Material Properties KXX, KYY, KZZ, DENS, C, ENTH, VISC, MU (VISC and MU used only if KEYOPT(7) = 1. Do not use ENTH with KEYOPT(8) = 1).

[ 本帖最后由 jxxansys 于 2009-6-6 15:28 编辑 ]

lilichunyan

这个倒是不用，我试过！

lilichunyan

呵呵，还真收敛了，发帖的时候稍微改了一下，但只是去掉的没用的啊！呵呵，很奇怪！
但是这个结果是不对的，那我怎么能为部分区域加上初始温度10000啊？

lilichunyan

原来的程序是我把热电部分分开了，现在是整个程序，不收敛的，您看看会不会？