高斯分布热源下的两层材料热分析

tammy

试着做了高斯分布热源下的两层材料瞬态热分析，上层为金刚石薄膜，下层为SI底。以下是命令流，如果有不合适的地方，欢迎大家提出、讨论！其中高斯分布体热源是通过GUI方式加的 :

（Parameters -->Functions-->Define/Edit 编辑函数为：
   1e13*exp(-({X}^2+{Y}^2))*exp({Z}-2e-4)*exp(-{TIME})
  再Parameters -->Functions-->Read From File ，并命名为 myload）  

1. /units,si
   
2. /prep7
   
3. ET,1,SOLID90
   
4. ET,2,MESH200,7
   
5. mp,kxx,1,148 !SI参数
   
6. mp,c,1,712
   
7. mp,dens,1,2330
   
8. !mp,kxx,2,2300 !Diamond参数
   
9. !mp,c,2,509
   
10. !mp,dens,2,3500
    
11. mp,kxx,2,170
    
12. mp,c,2,870
    
13. mp,dens,2,2700
    
14. CYLIND,1.25e-3,,2.0e-4,,45
    
15. CYLIND,1.25e-3,,2.0e-4,2.1e-4,,45
    
16. vglue,all
    
17. lesize,6,,,25,2
    
18. lesize,1,,,25,2
    
19. lesize,9,,,20,0.5
    
20. MSHAPE,0,2D
    
21. mshkey,1
    
22. MAT,1
    
23. TYPE,2
    
24. AMESH,4
    
25. allsel
    
26. type,1
    
27. mat,1
    
28. lesize,4,,,10
    
29. VSWEEP,1,4,5
    
30. lesize,15,,,25,2
    
31. lesize,6,,,25,2
    
32. lesize,21,,,2
    
33. MSHAPE,0,2D
    
34. mshkey,1
    
35. MAT,2
    
36. TYPE,2
    
37. AMESH,13
    
38. TYPE,1
    
39. MAT,2
    
40. LESIZE,13,,,10
    
41. VSWEEP,3,13,12
    
42. ALLSEL
    
43. SAVE  
    
44. /SOLU
    
45. AnTYPE,TRANS
    
46. TUNIF,293
    
47. kbc,0
    
48. asel,s,area,,7
    
49. asel,a,area,,11
    
50. asel,a,area,,3
    
51. nsla,s,1
    
52. sf,all,conv,10,293
    
53. allsel
    
54. nsel,s,loc,z,0,2.0e-4
    
55. csys,1
    
56. nsel,r,loc,x,0,1.25e-3
    
57. csys,0
    
58. bf,all,hgen,%myload%
    
59. allsel
    
60. time,1
    
61. deltime,2e-3,2e-3,8e-3
    
62. autots,on
    
63. outres,all,2
    
64. solve

复制代码

[ 本帖最后由 yejet 于 2006-8-12 15:40 编辑 ]

tammy

:)

lixiaoyu0323

程序好像有问题啊，体划分网格的时候自动退出了~

守护蓝天

你好，如果要在一个封闭的环形路径上加载的话，应该如何实现啊？