高等工程热力学-编程.docVIP

下载本文档

23
0
约1.3万字
约 22页
2018-02-02 发布于河北
举报
版权申诉

高等工程热力学-编程.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文档。上传文档

高等工程热力学-编程.doc

第一题题目：编写一个P－R方程已知 p、T 求比体积 v 的计算机程序。计算R134a的比体积，并将结果与文献值进行比较。将计算结果与文献值列表，并计算相对偏差；将计算值与文献值画在p –v 图上。 1、比较饱和线上的比体积； 2、比较等温线t=90℃的比体积； R134a: Tc=374.18K，pc=4.056MPa，ω=0.326， M=102.03 一、计算过程 P-R方程为：将P-R方程整理，用牛顿迭代法求v值 f = RgT(v2 + 2bv - b2) - a(v - b) - P(v - b)(v2 + 2bv - b2) f1 = RgT(2v + 2b) - a - P(3v2 - 3b2 + 2bv) 根据v(k+1) = v(k) - f / f1经过数次迭代后既可求得P、T所对应的v值二、源程序计算并比较饱和线上的比体积 Private Sub Command1_Click() Dim b, k, Rg As Double, i As Integer Dim p(10), t(10), f(10), f1(10), v(10), v1(10), v2(10), a(10), at(10), tr(10) Picture1.Cls: Picture2.Cls Tc = 374.18: Pc = 4056000: w = 0.326 k = 0.37464 + 1.5422 * w - 0.26992 * w ^ 2 t(1) = 246.7: t(2) = 262.96: t(3) = 273.73: t(4) = 282.08: t(5) = 288.87 t(6) = 294.7: t(7) = 299.87: t(8) = 304.47: t(9) = 308.65: t(10) = 312.54 For i = 1 To 10 p(i) = i * 10 ^ 5 tr(i) = t(i) / Tc a(i) = (1 + k * (1 - Sqr(tr(i)))) ^ 2 Rg = 8.3145 / 0.10203 at(i) = (0.45724 * Rg ^ 2 * Tc ^ 2 * a(i)) / Pc b = (0.0778 * Rg * Tc) / Pc v(i) = 1 For n = 1 To 5000 f(i) = Rg * t(i) * (v(i) ^ 2 + 2 * b * v(i) - b ^ 2) - at(i) * (v(i) - b) - p(i) * (v(i) - b) * (v(i) ^ 2 + 2 * b * v(i) - b ^ 2) f1(i) = Rg * t(i) * (2 * v(i) + 2 * b) - at(i) - p(i) * (3 * v(i) ^ 2 - 3 * b ^ 2 + 2 * b * v(i)) v1(i) = v(i) - f(i) / f1(i) If Abs(v1(i) - v(i)) 0.00001 Then v(i) = v1(i) End If Next n Picture1.Print Format(v1(i), ##0.000000) Picture1.Print Next i v2(1) = 0.189737: v2(2) = 0.098326: v2(3) = 0.066694: v2(4) = 0.050444: v2(5) = 0.040612 v2(6) = 0.03387: v2(7) = 0.029081: v2(8) = 0.025428: v2(9) = 0.022569: v2(10) = 0.020228 For i = 1 To 10 v2(i) = Abs((v1(i) - v2(i)) / v2(i)) * 100 Picture2.Print Format(v2(i), ##0.000) Picture2.Print Next i End Sub 计算并比较等温线t=80℃的比体积 Private Sub Command1_Click() Dim b, k, Rg As Double, i As Integer Dim p(10) As Double, t(10) As Double, f(10) As Double, f1(10) As Double, v(10) As Double, v1(10)