采用radfrac模拟反应精馏过程，反应选择动力学，但是出现不收敛和质量不守恒的问题..
本帖最后由 soaty 于 2013-5-26 00:41 编辑RT,我用RADFRAC模拟棕榈酸和异丙醇酯化反应指棕榈酸异丙酯的反应精馏过程，先尝试反应部分选择化学平衡模拟，没有错误，再在相同条件下，采用动力学，由于动力学方程的形式和我查到的表达形式不一致，所以用fortran语言编写了动力学方程导入了aspen中，但是运行的时候出现了不守恒和不收敛的问题，我改了迭代次数和迭代方法，以及蒸馏速率，仍然存在这个问题。求问，是什么原因导致的该问题。会不会是fortran那里的代入有问题？如何解决！！！！！！！ 附件里有fortran的文件，要用simulation engine生成obj文件。fortran采用的应该是6.5的版本，如截图里所示。 采用化学平衡模拟的时候也出现了不收敛，我改变Distillate rate或者是D/F的比值后就可以收敛了，弄成动力学减小D/F并不能收敛，但之前减小改值可以收敛。
水和异丙醇会共沸，求大神帮助谢谢！
已经解决了！感谢楚天湘水的帮助！
@楚天湘水
兄弟啊，来看看这个帖子可以解决不。 {:1106_365:}反应精馏以前也模拟过一次，貌似也是出了一些问题 @tdl522 ，呵呵，惭愧啊，其实我也懂的不多。
楼主能不能简单介绍下这个反应，有没有共沸生成？还有能不能把动力学方程贴出来，这样方便大家讨论 楚天湘水 发表于 2013-5-23 09:11
@tdl522 ，呵呵，惭愧啊，其实我也懂的不多。
楼主能不能简单介绍下这个反应，有没有共沸生成？还有能不能 ...
我把动力学传到附件里了，还有异丙醇和水是共沸体系。麻烦你看一下是什么原因。另外fortran语言不是我编写的，是我的导师编的，所以这一块的问题并不是很清楚。
看了楼主的动力学子程序，发现有这样一句话RATNET=RATE*REALB(1)*1.D-03，但我打开楼主的BKP文件里面并没有看到realb这个值，楼主的动力学子程序是从RAMTBE的那个子程序改过来，可能没有读懂那个程序，这个地方没有修改。
如果楼主真的需要用到动力学子程序的话，建议楼主把手册认真读个几遍吧，还有楼主写fortran程序的时候最好用个fortran语言的编辑器，这样可以避免很多语法错误。感觉楼主的PARAMETER(K_IPA=1)这句话跑到前面去了，aspen中用的fortran还是老版本，前面6个字符都不会编译。
今天我的编译器出了问题，没办法调试楼主的程序，楼主修改下那句话，试试吧。 楚天湘水 发表于 2013-5-23 09:58
看了楼主的动力学子程序，发现有这样一句话RATNET=RATE*REALB(1)*1.D-03，但我打开楼主的BKP文件里面并没有 ...
谢谢，我也感觉可能是FORTRAN的问题，不过我完全不懂这一块，是导师编的。这个fortran的版本也是导师电脑上的。那我先安装新版本的fortran试试可否。 soaty 发表于 2013-5-23 10:03
谢谢，我也感觉可能是FORTRAN的问题，不过我完全不懂这一块，是导师编的。这个fortran的版本也是导师电脑 ...
不是安装新版本fortran的问题，aspen中的fortran就是这个样子的 楚天湘水 发表于 2013-5-23 10:16
不是安装新版本fortran的问题，aspen中的fortran就是这个样子的
哦~明白您的意思了。那这个问题怎么解决啊。不懂。我以为是我安装的fortran的问题呢。
和楼主讨论了下，子程序已经修改好了，不收敛是因为共沸的问题，大家都来帮帮忙吧。 楚天湘水 发表于 2013-5-23 10:47
和楼主讨论了下，子程序已经修改好了，不收敛是因为共沸的问题，大家都来帮帮忙吧。
感谢兄弟修改子程序，现已调收敛。
楚天湘水 发表于 2013-5-23 10:47
和楼主讨论了下，子程序已经修改好了，不收敛是因为共沸的问题，大家都来帮帮忙吧。
感谢帮忙，下午的时候已经可以运行错误了。只是很奇怪，生成物水和酯均是从塔底出来，而塔顶出来的水很少，理论上水应该从塔顶出来啊。 我猜想是否与进料温度有关。略微试了下，结果还是不行。恳请回复~~ 请问你用那个版本的fortran编译器可以吗？是不是版本不能太高啊？我用的inter visual fortran2011，总是出现如下图的错误： weihuixiaoniao 发表于 2014-2-27 17:37
请问你用那个版本的fortran编译器可以吗？是不是版本不能太高啊？我用的inter visual fortran2011，总是出 ...
我用的aspen软件是2010年出的aspen plus v7.2 weihuixiaoniao 发表于 2014-2-27 17:43
我用的aspen软件是2010年出的aspen plus v7.2
恩。我用的fortran是以前版本的，好像是6.5的。aspen我用的v7.1 weihuixiaoniao 发表于 2014-2-27 17:43
我用的aspen软件是2010年出的aspen plus v7.2
刚才那个问题解决了，谢谢了，可以配套 借楼主的程序来学习学习 weihuixiaoniao 发表于 2014-3-6 23:29
刚才那个问题解决了，谢谢了，可以配套
我在做编译的时候也出现了您这个问题，请问您是如何解决的？我用的是V8.4 {:1106_365:}
页:
[1]