箱梁 温度梯度的模拟

现在在做一个箱梁温度分析的作业，但是有一些不明白的地方，
对于温度荷载的施加通常有两种方式：
其一是直接法，就是利用既可用于结构分析又可用于热分析的单元直接求解得出相应的应力；
其二是间接法：先采用热力学单元将温度转化为应力，然后将其施加到结构上(相当于外荷)，并将相应的热力学单元用结构力学单元替代，并求解得出相应的结构内力(应力)；
对于梁单元，如188，189我已经知道如何加温度效应，但是BFE只能施加线性的温度梯度，对于实体或者板壳单元的箱梁温度梯度的施加不是很明白，特别是如何施加非线性温度梯度，希望高手能发个例子上来，指导一二。
现贴一个简单的188单元，施加线性温度梯度荷载的例子（参看《基于ANSYS的桥梁结构分析》），以抛砖引玉。
fini
/cle
/prep7
height=1.5
uptemp=30
botemp=730
midtemp=(30+730)/2
temp_gradient=(uptemp-midtemp)/(height/2)
temp_z1=midtemp+temp_gradient
et,1,beam188
sectype,1,beam,i,,0
secoffset,cent
secdata,1,1,height,0.1,0.1,0.1,0,0,0,0
k,1
k,2,10
k,3,,10
l,1,2
latt,,,1,,3,,1
esize,1
lmesh,all
/eshape,1
/view,1,1,2,3
bfe,all,temp,1,midtemp,midtemp,temp_z1,midtemp
bfe,all,temp,5,midtemp,temp_z1
/pbf,temp,,1
eplo

[ 本帖最后由 fengpei001 于 2008-12-13 21:38 编辑 ]

怎么没有反应呢，自己先顶一个。

可以采用ANSYS耦合场求解能力。首先进行热分析求解温度分布，再进行结构分析就可以读入热分析求解任意时刻的温度分布作为结构的温度荷载。也就是上述的间接法。可是，如何才能得到自己想要的温度分布模式。现在学习中。。。。。。希望。。。。先发一个《ANSYS分析指南》，以便共同学习！

请问楼主，日照温度 是不是可以采用体荷载，本人在施加温度荷载时，却只能加到一个侧面上（如x方向），为什么不能同时加到y方向上呢？

可以用表格输入的方发来实现任何温度模式的输入哦～～
附件是我对一个小箱梁进行的折线温度加载，温度加载实现了折线分布
不过结构应力～～我还没有分析明白。     

下面是命令流  请多指教～～

1. fini
   
2. /CLEAR
   
3. /prep7
   
4. /title,Anlyse of kongxinban-side
   
5. !单位系统采用mm，N，MPa
   
6. et,1,42
   
7. ET,2,65
   
8. et,3,45 !模拟支座
   
9. MP,EX,1,3.45e4 !C50
   
10. MP,PRXY,1,0.2
    
11. MP,DENS,1,2600e-12
    
12. MP,ALPX,1,1e-5
    
13. MP,EX,2,2e5 !模拟支座
    
14. MP,PRXY,2,0.3
    
15. MP,DENS,2,7800e-12
    
16. MP,ALPX,2,1e-5
    
17. r,1
    
18. K,1, 542.5 , 0
    
19. K,2, 782.5 , 0
    
20. K,3, 782.5 , 230
    
21. K,4, 782.5 , 575
    
22. K,5, 702.5 , 655
    
23. K,6, 702.5 , 695
    
24. K,7, 757.5 , 750
    
25. K,8, 757.5 , 800
    
26. K,9, -1182.5 , 800
    
27. K,10, -1182.5 , 690
    
28. K,11, -782.5 , 600
    
29. K,12, -782.5 , 230
    
30. K,13, -782.5 , 0
    
31. K,14, -542.5 , 0
    
32. K,15, 542.5 , 110
    
33. K,16, 662.5 , 230
    
34. K,17, 662.5 , 440
    
35. K,18, 512.5 , 640
    
36. K,19, 0 , 680.009
    
37. K,20, -512.5 , 640
    
38. K,21, -662.5 , 440
    
39. K,22, -662.5 , 230
    
40. K,23, -542.5 , 110
    
41. 
    
42. L,1,2
    
43. *repeat,13,1,1
    
44. L,14,1
    
45. L,15,16
    
46. *repeat,3,1,1
    
47. larc,18,20,19
    
48. L,20,21
    
49. *repeat,3,1,1
    
50. L,23,15
    
51. L,1,15
    
52. L,16,3
    
53. L,18,5
    
54. L,11,20
    
55. L,12,22
    
56. L,14,23
    
57. 
    
58. AL,15,16,17,18,19,20,21,22
    
59. AL,14,23,22,28
    
60. AL,11,27,20,19,26
    
61. AL,25,4,3,24,16,17
    
62. AL,27,21,28,13,12
    
63. AL,24,2,1,23,15
    
64. AL,5,6,7,8,9,10,26,18,25
    
65. ALLSEL,ALL
    
66. 
    
67. esize,50
    
68. asel,s,,,5,6
    
69. asel,a,,,2
    
70. LCCAT,12,13
    
71. LCCAT,1,2
    
72. mshape,0,2d
    
73. mshkey,1
    
74. amesh,all
    
75. 
    
76. ALLSEL,ALL
    
77. LSEL,S,,,19,21,1
    
78. LCCAT,ALL
    
79. LSEL,S,,,15,17,1
    
80. LCCAT,ALL
    
81. amesh,1
    
82. 
    
83. ALLSEL,ALL
    
84. LCCAT,20,19
    
85. LCCAT,16,17
    
86. LCCAT,3,4
    
87. amesh,3,4
    
88. 
    
89. ALLSEL,ALL
    
90. LSEL,S,,,5,7
    
91. LCCAT,ALL
    
92. LSEL,S,,,10,18,8
    
93. LSEL,A,,,25,26
    
94. LCCAT,ALL
    
95. amesh,7
    
96. ALLSEL,ALL
    
97. LDELE,29,37,1
    
98. 
    
99. KGEN,2,14,,,,-60
    
100. L,14,24
     
101. ADRAG,14,,,,,,29
     
102. TYPE,1
     
103. amesh,8
     
104. allsel,all
     
105. 
     
106. k,100
     
107. k,101,,,-500
     
108. k,102,,,7500
     
109. lstr,100,101
     
110. lstr,100,102
     
111. asel,s,,,2,6,1
     
112. cm,a1,area
     
113. asel,s,,,1,7,1
     
114. cm,a2,area
     
115. /VIEW,1,1,1,1
     
116. 
     
117. type,2
     
118. mat,1
     
119. ESIZE,200
     
120. allsel,all
     
121. vdrag,7,,,,,,34
     
122. cmsel,s,a1
     
123. vdrag,all,,,,,,34
     
124. ESIZE,100
     
125. cmsel,s,a2
     
126. vdrag,all,,,,,,33
     
127. type,3
     
128. mat,2
     
129. ESIZE,100
     
130. allsel,all
     
131. vdrag,8,,,,,,33
     
132. 
     
133. allsel,all
     
134. esel,s,type,,1
     
135. aclear,all
     
136. 
     
137. allsel,all
     
138. NUMCMP,ALL
     
139. NUMMRG,ALL,,,,LOW
     
140. eplot
     
141. /dist,1,0.729,1
     
142. /triad,ltop
     
143. /image,save,16m_kongxinban_3d.bmp
     
144. SAVE
     
145. 
     
146. /solu
     
147. time,1
     
148. tref,0 !定义原始温度
     
149. *dim,htemp,table,3,1,1,Y,TEMP, !定义表格数据模拟梯度温度
     
150. htemp(1)=0,0,15
     
151. htemp(1,0)=0,400,800
     
152. 
     
153. NSEL,S,LOC,Y,0,800
     
154. BF,ALL,TEMP, %HTEMP% !调用表格数据加载梯度温度
     
155. /PBf,temp,,1 !显示温度荷载
     
156. /VIEW,1,1,1,1
     
157. /ANG,1
     
158. /REP,FAST
     
159. 
     
160. allsel,all
     
161. nsel,s,loc,y,-60
     
162. nsel,r,loc,z,-200,-500
     
163. d,all,all !定义支座约束
     
164. nsel,s,loc,Z,7500
     
165. d,all,uz !定义跨中对称约束
     
166. allsel,all
     
167. outpr,basic,last
     
168. solve
     
169. /dscale,1,1 !显示图形中不显示变形

复制代码

终于找到了一个日照温度的，哈哈