传热学多媒体教学辅助系统（赵镇南）例题中的数值计算源程序.docVIP

下载本文档

4
0
约1.23万字
约 11页
2017-10-13 发布于浙江
举报
版权申诉

传热学多媒体教学辅助系统（赵镇南）例题中的数值计算源程序.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

例题中的数值计算源程序目录例题3-10：热力发动机壁面的冷却例题3-11：椭圆管矩形翅片的效率例题3-12：梯形散热肋的温度分布与肋效率例题4-11：一维对流边界条件下平壁的非稳态导热源程序例题3-10：热力发动机壁面的冷却 c program for ex3-10 dimension t (65),tt (65),a (65) data dx,dy,cond,tf1,h1,tf2,h2,tf3,h3/.002,.001,4.,1600.,500., * 450.,40.,350.,1500./ data eps,N/0.01,200/ write(*,100) 100 format(1x,input cond:F6.2) c read(*,200) cond 200 format(f6.2) write(*,101)dx,dy,cond,tf1,h1,tf2,h2,tf3,h3 101 format(1x,dx/dy=,2f6.3,2x,cond=,f6.2/ * 1x,tfi/hi:,6(f6.1,1x)) Bi1=h1*dx/cond Bi2=h2*dx/cond Bi3=h3*dx/cond it=0 MN=65 c area of boundary for node points a(1)=0.25 do 51 J=2,65 51 a(J)=1.0 do 52 J=2,6 52 a(J)=0.5 do 53 J=7,35,7 53 a(J)=0.5 do 54 J=8,29,7 54 a(J)=0.5 do 55 J=30,32 55 a(J)=0.5 a(33)=0.75 a(36)=0.5 a(39)=0.5 a(42)=0.5 do 56 J=38,44,3 56 a(J)=0.5 a(45)=0.25 do 57 J=46,48 57 a(J)=0.5 a(49)=0.75 a(51)=0.5 a(58)=0.5 a(52)=0.5 a(59)=0.25 do 58 J=60,64 58 a(J)=0.5 a(65)=0.25 c assign initial temp. value do 10 I=1,MN 10 t(I)=750.0 12 do 20 J=1,MN tt(J)=t(J) 20 continue c nodal eq. 14 continue t(1)=(0.5*t(2)+2.*t(8)+Bi1*tf1)/(2.5+Bi1) do 2 i=2,6 2 t(i)=(t(i-1)+t(i+1)+8.*t(i+7)+4.*Bi1*tf1)/(10.+4.*Bi1) t(7)=(0.5*t(6)+2.*t(14)+Bi1*tf1)/(2.5+Bi1) t(8)=(t(1)+t(15)+0.5*t(9))/2.5 do 3 i=9,13 3 t(i)=(0.5*t(i-1)+0.5*t(i+1)+2.*t(i-7)+2.*t(i+7))/5. t(14)=(t(7)+t(21)+0.5*t(13))/2.5 t(15)=(t(8)+t(22)+0.5*t(16))/2.5 do 4 i=16,20 4 t(i)=(0.5*t(i-1)+0.5*t(i+1)+2.*t(i-7)+2.*t(i+7))/5. t(21)=(t(14)+t(28)+0.5*t(20))/2.5 t(22)=(t(15)+t(29)+0.5*t(23))/2.5 do 5 i=23,27 5 t(i)=(0.5*t(i-1)