用shell131做的复合材料热分析算例（转贴）

AaronSpark

本算例为4层石墨环氧反对称角铺层方板[45/-45/45/-45]，边长0.08m，厚度0.008m，加载热流于上表面中心区，面积 0.0004m*m，大小为100000W/m*m,材料主轴热导率K11=30W/mK,k22=k33=1.1，对流换热系数a=30W/m*mK.

1. finish
   
2. /clear
   
3.   
   
4. /prep7
   
5. et,1,shell131,,,1,8!热分析用4节点分层壳单元
   
6. sectype,1,shell
   
7. mp,kxx,1,30        !45度方向材料属性   
   
8. mp,kyy,1,1.1
   
9. mp,kzz,1,1.1
   
10. mp,c,1,1000
    
11. mp,dens,1,1389.2
    
12. mp,kxx,2,30        !－45度方向材料属性
    
13. mp,kyy,2,1.1
    
14. mp,kzz,2,1.1
    
15. mp,c,2,1000
    
16. mp,dens,2,1389.2
    
17. secdata,0.001,1,45 !进行各个铺层厚度、材料、方向设定
    
18. secdata,0.001,2,-45
    
19. secdata,0.001,1,45
    
20. secdata,0.001,2,-45
    
21. secdata,0.001,1,45
    
22. secdata,0.001,2,-45
    
23. secdata,0.001,1,45
    
24. secdata,0.001,2,-45
    
25. n,1,               !建模
    
26. n,17,0.08
    
27. fill,1,17
    
28. ngen,17,17,1,17,1,,0.005
    
29. e,1,2,19,18
    
30. egen,16,1,1
    
31. egen,16,17,1,16
    
32. esel,s,,,136,137
    
33. esel,a,,,120,121
    
34. sfe,all,2,hflux,,1e5  !在顶面加载热流载荷
    
35. esel,all
    
36. finish
    
37. !!!!!!!!!!!!!!!!!static
    
38. /solu!稳态
    
39. antype,static
    
40. tref,20
    
41. sfe,all,1,conv,1,30
    
42. sfe,all,1,conv,2,20   
    
43. esel,u,,,136,137
    
44. esel,u,,,120,121      !对未加载单元顶面、所有单元底面  
    
45. sfe,all,2,conv,1,30   !定义边界的空气对流换热
    
46. sfe,all,2,conv,2,20   
    
47. esel,all
    
48. outres,all,all
    
49. solve
    
50. /post1
    
51. plnsol,tbot         !分层显示各个层底温度分布
    
52. plnsol,te2          !由底面到顶面  
    
53. plnsol,te3
    
54. plnsol,te4
    
55. plnsol,te5
    
56. plnsol,te6
    
57. plnsol,te7
    
58. plnsol,te8
    
59. plnsol,ttop
    
60.   
    
61. path,ppp,9         !定义沿X方向半对称轴路径
    
62. ppath,1,145
    
63. ppath,2,146
    
64. ppath,3,147
    
65. ppath,4,148
    
66. ppath,5,149
    
67. ppath,6,150
    
68. ppath,7,151
    
69. ppath,8,152
    
70. ppath,9,153
    
71. pdef,pop,ttop     !显示本路径顶面温度变化
    
72. plpath,pop

复制代码

转自simwe

seashoregang

好东西