刚才才看到短消息，不好意思，那几个发短消息的兄弟

diyifan001

      本人是在马后炮创建初期，在网上搜索的时候发现炮网的，然后加入的炮网的qq群，感受到了炮网的分享精神。但由于本人经常外面跑，对炮网的贡献不大，也没太注意过短消息。在马甲草堂发过关于编程计算相平衡的帖子（透露一下，本人的第一个宏伟目标是编个只有经常能用到的模块的小的模拟软件-已经搁浅，第二个宏伟目标是实现局部的无政府化，我是个无政府主义者），还给炮网编过一个hysys导出到国内设计院格式execl的小工具（个人认为比网上下到的好用多了，可惜论坛出问题，消失了）。今天才看到几个兄弟发的短消息，说要那个个MATLAB计算相平衡的源代码，可惜的是我不注意保存东西，找不到了。刚开始没贴源代码的原因是，我总觉得对这个感兴趣的人不多。

    不过我的思路还在，如果大家还有兴趣，就跟个贴，人多的话，我就重新按照思路一边发帖子，一边写代码（代码只有matlab或者C#语言可选，其他语言我真不会），然后大家一起讨论。

diyifan001

果然没人理，我觉得我这种低调的人，马哥应该给我开个版，专门讨论extention，炮友们分享自编小工具的版面。

lidaxue1987

回复 diyifan001 的帖子

楼主你好，很高兴看到你的这个帖子以及之前的那个，我现在也是对相平衡有关的泡露点计算很感兴趣，现在的问题是，我在matlab编程，回归二元交互作用参数的时候，发现程序走的不对。唉

lidaxue1987

回复 lidaxue1987 的帖子

我采用的等温泡点计算流程图是陈新志编写的化工热力学  第三版  p103上的，当然  第二版也有，但不知道具体在哪一页。还请楼主留意这个帖子的信息。谢谢啦

lidaxue1987

楼主曾经在别的帖子里提到过的关于压力迭代或者温度迭代的具体格式在陈的书里面都可以找到具体的表达形式。

diyifan001

书上的我也参考过，但是问题是用纯粹的牛顿迭代，有时候会不收敛，所以实践的时候，我在迭代求解过程中采用了阻尼系数，就是迭代目标差值一旦发生正负值的改变，就调整改变量为1/2n   ，N是正负值改变的次数。
关于回归二元交互参数，我没实践过，您说个求解思路，我也试着编编。还请您关注帖子

diyifan001

看了您的介绍，我翻了陈老师的书，给出的框图和其他书上基本差不多，没有给出调整迭代值（T或者P）的方法，还有未给出防止迭代发散的方法，迭代次数应该在比较少的次就应该收敛。因为我刚开始的时候就老收敛到错误的温度。
我在网上找到的书是2001版的，不知道新版的有没有。没找到书。

我同意好友请求，系统老是说我没有权限，晕。呵呵

lidaxue1987

回复 diyifan001 的帖子

至于那个二元交互系数的回归，我们以范德华单流体混合规则为例，附件中的doc是大致的回归流程图，您可以参考一下

lidaxue1987

回复 diyifan001 的帖子

陈老师的书里面 提到了等温泡点的压力及气相组成的迭代格式您可以参考一下

diyifan001

调整迭代值，我用的是割线法加阻尼系数，陈老师的调整方法我没试过，不知道收敛性如何，如果我有时间的话，我用两种调整方法把程序重新写一遍贴上来，咱们在具体讨论下，看那个收敛快些。
感谢您的资料，我正在下，有点慢。网络不好

lidaxue1987

回复 diyifan001 的帖子

好的  感谢楼主。我也在看这部分的知识，静候佳音！

diyifan001

还有，不知道陈老师书里的调整迭代值有没有通用性，比如用在求解泡点温度上。

lidaxue1987

回复 diyifan001 的帖子

根据书上的流程图 陈老师编写的热力学计算程序，楼主可以参考，可惜  是exe可执行程序，查看不到原始代码

diyifan001

ok
我这里下不下来
晚上回去学习学习，然后我重新开个帖子，编点程序就贴上来，到时候您给看看有没有改进的地方。您看用C#行不行，因为最后想做个简单的界面便于对比计算速度。

lidaxue1987

回复 diyifan001 的帖子

C# 我不懂啊，呵呵，感觉还是matlab  还行  呵呵

diyifan001

回复 lidaxue1987 的帖子

ok
我有空就在这个版面重新开个帖子
我一边贴思路，一边贴代码
有优化和改进的地方还请您指出
这个帖子就结了吧
呵呵

小飞侠

好帖，必须顶啊

xaogny

我原来也算过这个的，支持楼主

hebinlin

就是啊！楼主你说你编写过hysys现在还有没？？？能不能传一个上来，在下现在正想学习他，就是没有。。。多谢楼主。。。。。。

diyifan001

回复 hebinlin 的帖子

是个hysys导出工具，导出到个模板，以前就发在这里，论坛出故障，被老马搞没了（这厮）。如果你需要的话，可以为你量身定做一个，把你需要的导出格式放上来，我有空做一下。