箱梁 温度梯度的模拟

fengpei001 发表于 2008-12-13 21:37:38 | 显示全部楼层 | 阅读模式
现在在做一个箱梁温度分析的作业,但是有一些不明白的地方,
对于温度荷载的施加通常有两种方式:
其一是直接法,就是利用既可用于结构分析又可用于热分析的单元直接求解得出相应的应力;
其二是间接法:先采用热力学单元将温度转化为应力,然后将其施加到结构上(相当于外荷),并将相应的热力学单元用结构力学单元替代,并求解得出相应的结构内力(应力);
对于梁单元,如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 编辑 ]
精彩评论倒序浏览

7615查看5评论

fengpei001 发表于 2008-12-16 13:29:00
怎么没有反应呢,自己先顶一个。
举报 回复
fengpei001 发表于 2008-12-16 21:47:28
可以采用ANSYS耦合场求解能力。首先进行热分析求解温度分布,再进行结构分析就可以读入热分析求解任意时刻的温度分布作为结构的温度荷载。也就是上述的间接法。可是,如何才能得到自己想要的温度分布模式。现在学习中。。。。。。希望。。。。先发一个《ANSYS分析指南》,以便共同学习!

ANSYS热分析指南.rar

93.5 KB, 下载次数: 61

举报 回复
我是飞沙... 发表于 2009-5-10 21:24:11
请问楼主,日照温度 是不是可以采用体荷载,本人在施加温度荷载时,却只能加到一个侧面上(如x方向),为什么不能同时加到y方向上呢?
举报 回复
murongcao 发表于 2010-4-27 09:12:26
本帖最后由 murongcao 于 2010-4-27 09:18 编辑

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

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

  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. L,1,2
  42. *repeat,13,1,1
  43. L,14,1
  44. L,15,16
  45. *repeat,3,1,1
  46. larc,18,20,19
  47. L,20,21
  48. *repeat,3,1,1
  49. L,23,15
  50. L,1,15
  51. L,16,3
  52. L,18,5
  53. L,11,20
  54. L,12,22
  55. L,14,23

  56. AL,15,16,17,18,19,20,21,22
  57. AL,14,23,22,28
  58. AL,11,27,20,19,26
  59. AL,25,4,3,24,16,17
  60. AL,27,21,28,13,12
  61. AL,24,2,1,23,15
  62. AL,5,6,7,8,9,10,26,18,25
  63. ALLSEL,ALL

  64. esize,50
  65. asel,s,,,5,6
  66. asel,a,,,2
  67. LCCAT,12,13
  68. LCCAT,1,2
  69. mshape,0,2d
  70. mshkey,1
  71. amesh,all

  72. ALLSEL,ALL
  73. LSEL,S,,,19,21,1
  74. LCCAT,ALL
  75. LSEL,S,,,15,17,1
  76. LCCAT,ALL
  77. amesh,1

  78. ALLSEL,ALL
  79. LCCAT,20,19
  80. LCCAT,16,17
  81. LCCAT,3,4
  82. amesh,3,4

  83. ALLSEL,ALL
  84. LSEL,S,,,5,7
  85. LCCAT,ALL
  86. LSEL,S,,,10,18,8
  87. LSEL,A,,,25,26
  88. LCCAT,ALL
  89. amesh,7
  90. ALLSEL,ALL
  91. LDELE,29,37,1

  92. KGEN,2,14,,,,-60
  93. L,14,24
  94. ADRAG,14,,,,,,29
  95. TYPE,1
  96. amesh,8
  97. allsel,all

  98. k,100
  99. k,101,,,-500
  100. k,102,,,7500
  101. lstr,100,101
  102. lstr,100,102
  103. asel,s,,,2,6,1
  104. cm,a1,area
  105. asel,s,,,1,7,1
  106. cm,a2,area
  107. /VIEW,1,1,1,1

  108. type,2
  109. mat,1
  110. ESIZE,200
  111. allsel,all
  112. vdrag,7,,,,,,34
  113. cmsel,s,a1
  114. vdrag,all,,,,,,34
  115. ESIZE,100
  116. cmsel,s,a2
  117. vdrag,all,,,,,,33
  118. type,3
  119. mat,2
  120. ESIZE,100
  121. allsel,all
  122. vdrag,8,,,,,,33

  123. allsel,all
  124. esel,s,type,,1
  125. aclear,all

  126. allsel,all
  127. NUMCMP,ALL
  128. NUMMRG,ALL,,,,LOW
  129. eplot
  130. /dist,1,0.729,1
  131. /triad,ltop
  132. /image,save,16m_kongxinban_3d.bmp
  133. SAVE

  134. /solu
  135. time,1
  136. tref,0 !定义原始温度
  137. *dim,htemp,table,3,1,1,Y,TEMP, !定义表格数据模拟梯度温度
  138. htemp(1)=0,0,15
  139. htemp(1,0)=0,400,800

  140. NSEL,S,LOC,Y,0,800
  141. BF,ALL,TEMP, %HTEMP% !调用表格数据加载梯度温度
  142. /PBf,temp,,1 !显示温度荷载
  143. /VIEW,1,1,1,1
  144. /ANG,1
  145. /REP,FAST

  146. allsel,all
  147. nsel,s,loc,y,-60
  148. nsel,r,loc,z,-200,-500
  149. d,all,all !定义支座约束
  150. nsel,s,loc,Z,7500
  151. d,all,uz !定义跨中对称约束
  152. allsel,all
  153. outpr,basic,last
  154. solve
  155. /dscale,1,1 !显示图形中不显示变形
复制代码

评分

参与人数 1堡币 +2 威望 +1 收起 理由
水幽寒 + 2 + 1

查看全部评分

举报 回复
shunmei_sm... 发表于 2015-3-9 21:42:22
终于找到了一个日照温度的,哈哈
举报 回复
您需要登录后才可以回帖 登录 | 注册

本版积分规则

 
 
  • QQ:56984982
  • 点击这里给我发消息
    电话:13527553862
    站务咨询群桥头堡站务咨询桥梁专业交流群:
    中国桥梁专业领袖群
    工作时间
    8:00-18:00