多效蒸发课程设计php程序

 
从邮件里面整理到一些大学里的资料，做个blog备份

多效蒸发是用于工业中的浓缩方案。

1.多效蒸发本身就存在很多不确定因素，结果仅供参考！
2.感谢老师的指导！
3.同时特别感谢易晓霆同学在本程序编写过程中的指导和纠错！

后话：

网络确实很有趣，需要大量的时间去接触，然而但你深入的同时，或许已经身不由己了。无法接受这个残酷的现实。当你在忙忙碌碌中突然醒悟的同时，一切都太晚了。失去的东西却是那么多，令人感到心痛，然而悲剧还在一次次重演。曾看到一个故事：一个经历了很多坎坷的老人对一个刚刚步入社会的年轻人说：别走这条路了，这条路不值得去走！年轻人说：这条路我要走！我可以走的更好！当这个年轻人执拗的走完时，却发现自己只是重复了老人的历史。然而网络也是这样，选择了网络就回不了头了！

部分程序

<?php
function get_t($p_now,$p_1,$p_2,$get_1,$get_2){
        $ap1=$p_2-$p_1;
        $ap2=$p_now-$p_1;
        $ap3=$ap2/$ap1;
        $ap4=$get_2-$get_1;
        $ap5=$ap4*$ap3;
        return $get_1+$ap5;
}
function do_ai($x){
        return 1.78*$x+6.22*$x*$x;
        }

$br=&apos;<hr>&apos;;
//进料
$f =1000*1000*1.05/24; //请在左边设置原始数据 （单位：kg/h）
echo "f=$f".$br;
//总蒸发水量
$x0=0.13;$x4=0.64;   //请在左边设置原始数据
$w=$f*(1-$x0/$x4);
echo  "w=$w".$br;
//各效抽用额外蒸汽量
$e1=0.06*$f;
$e2=0.07*$f;
$e3=0.08*$f;
echo "e1=$e1   e2=$e2  e3=$e3".$br;
$w4=($w-$e1-2*$e2-3*$e3)/4;
$w3=$w4+$e3;
$w2=$w3+$e2;
$w1=$w2+$e1;
echo "w1=$w1   w2=$w2  w3=$w3 w4=$w4".$br;
$x1=$f*$x0/($f-$w1);
$x2=$f*$x0/($f-$w1-$w2);
$x3=$f*$x0/($f-$w1-$w2-$w3);
$x4=$f*$x0/($f-$w1-$w2-$w3-$w4);
echo "x1=$x1   x2=$x2  x3=$x3 x4=$x4".$br;
//假设等压强分配
$p0=300;$p4=20;
$ap=$p0-$p4;
$p1=$p0-$ap/4;
$p2=$p1-$ap/4;
$p3=$p2-$ap/4;
$p4=$p3-$ap/4;
echo "p0=$p0 p1=$p1  p2=$p2  p3=$p3 p4=$p4".$br;
//查表
$tg0=get_t(300,270.25,313.11,130,135);
$tg1=get_t(230,198.64,232.19,120,125);
$tg2=get_t(160,143.31,169.11,110,115);
$tg3=get_t(90,84.556,101.33,95,100);
$tg4=get_t(20,19.923,25.014,60,65);
echo "tg0=$tg0 tg1=$tg1  tg2=$tg2  tg3=$tg3 tg4=$tg4".$br;//0c
$h0=get_t(300,270.25,313.11,2723.9,2731.0);
$h1=get_t(230,198.64,232.19,2708.9,2716.4);
$h2=get_t(160,143.31,169.11,2693.4,2701.3);
$h3=get_t(90,84.556,101.33,2668.7,2677.0);
$h4=get_t(20,19.923,25.014,2606.3,2615.5);
echo "h0=$h0 h1=$h1  h2=$h2  h3=$h3 h4=$h4".$br;   //kj/kg
$r0=get_t(300,270.25,313.11,2177.6,2163.3);
$r1=get_t(230,198.64,232.19,2205.2,2191.8);
$r2=get_t(160,143.31,169.11,2232.0,2219.0);
$r3=get_t(90,84.556,101.33,2270.9,2258.4);
$r4=get_t(20,19.923,25.014,2355.1,2343.4);
echo "r0=$r0 r1=$r1  r2=$r2  r3=$r3 r4=$r4".$br; //kj/kg
//计算公式 沸点升高引起的温度差损失
$da1=do_ai($x1);
$da2=do_ai($x2);
$da3=do_ai($x3);
$da4=do_ai($x4);
echo "da1=$da1  da2=$da2  da3=$da3 da4=$da4".$br; //0c
//暂取l=3m ,正常液位为1/3l=1m,查表
/*-------------------人工查询--------------------------------*/
$daa1=1.6;
$daa2=get_t($tg2,100,110,2.9,2.2);
$daa3=get_t($tg3,90,100,4.1,3.0);
$daa4=get_t($tg4,60,70,11.0,8.1);
echo "daa1=$daa1  daa2=$daa2  daa3=$daa3 daa4=$daa4".$br; //0c
/*-------------------人工查询--------------------------------*/
$t1=$tg1+$da1+$daa1;
$t2=$tg2+$da2+$daa2;
$t3=$tg3+$da3+$daa3;
$t4=$tg4+$da4+$daa4;
echo "t1=$t1  t2=$t2  t3=$t3 t4=$t4".$br; //0c
//由x0=0.13 t1=126.810c查表
$cpw=4.187;
$cp0=get_t($x0,0.1,0.2,0.963,0.927)*$cpw;
echo "cp0=$cp0".$br;  //kj/kg0c
$af1=1;
$af2=1;
$af3=1;
$af4=1;
$beita1=0;
$beita2=($t1-$t2)/($h2-$cpw*$t2);
$beita3=($t2-$t3)/($h3-$cpw*$t3);
$beita4=($t3-$t4)/($h4-$cpw*$t4);
echo "beita1=$beita1  beita2=$beita2  beita3=$beita3 beita4=$beita4".$br; //0c
$y1=$y2=$y3=$y4=0.97;
$a1=$af1*$y1;
$a2=($af2*$a1-$cpw*$beita2*$a1)*$y2;
$temp=$a1+$a2;
$a3=($af3*$a2-$cpw*$beita3*$temp)*$y3;
$temp=$a1+$a2+$a3;
$a4=($af4*$a3-$cpw*$beita4*$temp)*$y4;
echo "a1=$a1  a2=$a2  a3=$a3 a4=$a4".$br; //0c
$b1=$f*$cp0*$beita1*$y1;
$temp=$b1;
$b2=$y2*$af2*($b1-$e1)+$y2*$beita2*($f*$cp0-$cpw*$temp);
$temp=$b1+$b2;
$b3=$y3*$af3*($b2-$e2)+$y3*$beita3*($f*$cp0-$cpw*$temp);
$temp=$b1+$b2+$b3;
$b4=$y4*$af4*($b3-$e3)+$y4*$beita4*($f*$cp0-$cpw*$temp);
echo "b1=$b1  b2=$b2  b3=$b3 b4=$b4".$br; //0c
$biga=$a1+$a2+$a3+$a4;
$bigb=$b1+$b2+$b3+$b4;
$d1=($w-$bigb)/$biga;
$_w1=$a1*$d1+$b1;
$_w2=$a2*$d1+$b2;
$_w3=$a3*$d1+$b3;
$_w4=$a4*$d1+$b4;
echo "bigb=$bigb biga=$biga _w1=$_w1  _w2=$_w2  _w3=$_w3 _w4=$_w4".$br; //0c
$_x1=$f*$x0/($f-$_w1);
$_x2=$f*$x0/($f-$_w1-$_w2);
$_x3=$f*$x0/($f-$_w1-$_w2-$_w3);
$_x4=$f*$x0/($f-$_w1-$_w2-$_w3-$_w4);
echo "_x1=$_x1   _x2=$_x2  _x3=$_x3 _x4=$_x4".$br;
$daa_all=$tg0-$tg4-($da1+$da2+$da3+$da4)-($daa1+$daa2+$daa3+$daa4)-3;
echo "daa_all=$daa_all".$br; //0c
//加热蒸汽消耗量
$d2=$_w1-$e1;
$d3=$_w2-$e2;
$d4=$_w3-$e3;
echo "d1=$d1   d2=$d2  d3=$d3 d4=$d4".$br;
//带入数据
$k1=3000;$k2=1900;$k3=1000;$k4=520;
//有效温度差分配
$drk1=$d1*$r0*1000/($k1*3600);
$drk2=$d2*$r1*1000/($k2*3600);
$drk3=$d3*$r2*1000/($k3*3600);
$drk4=$d4*$r3*1000/($k4*3600);
echo "drk1=$drk1   drk2=$drk2  drk3=$drk3 drk4=$drk4".$br;
$drk_all=$drk1+$drk2+$drk3+$drk4;//m2*k
echo "drk_all=$drk_all".$br; //0c
$at1=$drk1*$daa_all/$drk_all;
$at2=$drk2*$daa_all/$drk_all;
$at3=$drk3*$daa_all/$drk_all;
$at4=$drk4*$daa_all/$drk_all;
echo "at1=$at1   at2=$at2  at3=$at3 at4=$at4".$br;
$adt1=$drk1/$at1;
$adt2=$drk2/$at2;
$adt3=$drk3/$at3;
$adt4=$drk4/$at4;
echo "adt1=$adt1   adt2=$adt2  adt3=$adt3 adt4=$adt4".$br;
echo "根据_x查表 _x1=$_x1   _x2=$_x2  _x3=$_x3 _x4=$_x4
这里根据公式".$br;
$dap1=do_ai($_x1);
$dap2=do_ai($_x2);
$dap3=do_ai($_x3);
$dap4=do_ai($_x4);
echo "dap1=$dap1 dap2=$dap2 dap3=$dap3 dap4=$dap4".$br;
/*-------------------人工查询--------------------------------*/
$tgp4=$tg4;
$daap4=get_t($tgp4,70,60,8.1,11.0);
$tgp3=$tgp4+1+$at4+$daap4+$dap4;
echo "查表tgp3=$tgp3 _x3=$_x3".$br;
$daap3=get_t($tgp3,100,90,3.0,4.1);
$tgp2=$tgp3+1+$at3+$daap3+$dap3;
echo "查表tgp2=$tgp2 _x2=$_x2".$br;
$daap2=get_t($tgp2,120,110,1.7,2.2);
$tgp1=$tgp2+1+$at2+$daap2+$dap2;
echo "查表tgp1=$tgp1 _x1=$_x1".$br;
$daap1=get_t($tgp1,130,120,1.3,1.6);
$tgp0=$tgp1+$at1+$daap1+$dap1;
echo "tgp0=$tgp0".$br;
//echo "$tgp3=$tgp4+1+$at4+$daap4+$dap4;";
//$tgp3=$tgp4+1+$at4+$daap4+$dap4;
echo "tgp1=$tgp1 tgp2=$tgp2 tgp3=$tgp3 tgp4=$tgp4".$br;
echo "daap1=$daap1 daap2=$daap2 daap3=$daap3 daap4=$daap4".$br;
$asd1=$daap1+$dap1;
$asd2=$daap2+$dap2;
$asd3=$daap3+$dap3;
$asd4=$daap4+$dap4;
echo "asd1=$asd1   asd2=$asd2  asd3=$asd3 asd4=$asd4".$br;
/*-------------------人工查询--------------------------------*/
$last_t4=$t4;
$last_t3=$tgp3+$daap3+$dap3;
$last_t2=$tgp2+$daap2+$dap2;
$last_t1=$tgp1+$daap1+$dap1;
echo "last_t1=$last_t1 last_t2=$last_t2
last_t3=$last_t3 last_t4=$last_t4".$br;
$word=&apos;一&apos;;
include(&apos;fs.php&apos;);
$word=&apos;二&apos;;
include(&apos;fs.php&apos;);
$word=&apos;三&apos;;
include(&apos;fs.php&apos;);
function hd($w,$u,$p){
        $temp=3.14*$u*$p*3600;
        return sqrt($temp);
        }
        $hd1=hd($w1,48,1.2425);
        $hd2=hd($w2,46,0.8213);
        $hd3=hd($w3,44,0.4163);
        $hd4=hd($w4,42,0.1307);
        echo "hd1=$hd1   hd2=$hd2  hd3=$hd3 hd4=$hd4".$br;
$p1=get_t(126.5884,125,130,1.296,1.494);
$p2=get_t(113.8908,110,115,0.8254,09635);
$p3=get_t(96.5069,95,100,0.5039,0.5970);
$p4=get_t(74.7406,70,75,0.1979,0.2416);
 echo "t1=$t1   t2=$t2  t3=$t3 t4=$t4".$br;
  echo "p1=$p1   p2=$p2  p3=$p3 p4=$p4".$br;
  function d($d,$u,$p){
        $temp=3.14*$p*3600*30;
        return sqrt($d*4/$temp);
        }
        $d1=d(1.3890,30,48,1.2425);
        $d2=d($d1,46,0.8213);
        $d3=d($w3,44,0.4163);
        $d4=d($w4,42,0.1307);

?>

 

 

<?php
//查表入录
/*
$tg0=get_t($tg0,120,110,130,135);
$tg1=get_t($tgp1,198.64,232.19,120,125);
$tg2=get_t(160,143.31,169.11,110,115);
$tg3=get_t(90,84.556,101.33,95,100);
$tg4=get_t(20,19.923,25.014,60,65);
echo "tg0=$tg0 tg1=$tg1  tg2=$tg2  tg3=$tg3 tg4=$tg4".$br;//0c
*/

echo "------------------进入第 $word 次复算！！！！！----------------".$br;
echo "tgp1=$tgp1 tgp2=$tgp2 tgp3=$tgp3 tgp4=$tgp4".$br;
$h0=get_t($tg0,130,135,2723.9,2731.0);
$h1=get_t($tgp1,120,125,2708.9,2716.4);
$h2=get_t($tgp2,110,115,2693.4,2701.3);
$h3=get_t($tgp3,90,95,2659.9,2668.7);
$h4=get_t($tgp4,60,65,2606.3,2615.5);
echo "h0=$h0 h1=$h1  h2=$h2  h3=$h3 h4=$h4".$br;   //kj/kg
$r0=get_t($tg0,130,135,2177.6,2163.3);
$r1=get_t($tgp1,120,125,2205.2,2191.8);
$r2=get_t($tgp2,110,115,2232.0,2219.0);
$r3=get_t($tgp3,90,95,2283.1,2270.9);
$r4=get_t($tgp4,60,65,2355.1,2343.4);
echo "r0=$r0 r1=$r1  r2=$r2  r3=$r3 r4=$r4".$br; //kj/kg
/*
$t1=$tg1+$da1+$daa1;
$t2=$tg2+$da2+$daa2;
$t3=$tg3+$da3+$daa3;
$t4=$tg4+$da4+$daa4;
echo "t1=$t1  t2=$t2  t3=$t3 t4=$t4".$br; //0c
//由x0=0.13 t1=$last_t1 0c查表
$cpw=4.187;
$cp0=get_t($x0,0.1,0.2,0.963,0.927)*$cpw;
echo "cp0=$cp0".$br;  //kj/kg0c
 */
$af1=1;
$af2=1;
$af3=1;
$af4=1;
$beita1=0;
$beita2=($last_t1-$last_t2)/($h2-$cpw*$last_t2);
$beita3=($last_t2-$last_t3)/($h3-$cpw*$last_t3);
$beita4=($last_t3-$last_t4)/($h4-$cpw*$last_t4);
echo "beita1=$beita1  beita2=$beita2  beita3=$beita3 beita4=$beita4".$br; //0c
$y1=$y2=$y3=$y4=0.97;
$a1=$af1*$y1;
$a2=($af2*$a1-$cpw*$beita2*$a1)*$y2;
$temp=$a1+$a2;
$a3=($af3*$a2-$cpw*$beita3*$temp)*$y3;
$temp=$a1+$a2+$a3;
$a4=($af4*$a3-$cpw*$beita4*$temp)*$y4;
echo "a1=$a1  a2=$a2  a3=$a3 a4=$a4".$br; //0c
$b1=$f*$cp0*$beita1*$y1;
$temp=$b1;
$b2=$y2*$af2*($b1-$e1)+$y2*$beita2*($f*$cp0-$cpw*$temp);
$temp=$b1+$b2;
$b3=$y3*$af3*($b2-$e2)+$y3*$beita3*($f*$cp0-$cpw*$temp);
$temp=$b1+$b2+$b3;
$b4=$y4*$af4*($b3-$e3)+$y4*$beita4*($f*$cp0-$cpw*$temp);
echo "b1=$b1  b2=$b2  b3=$b3 b4=$b4".$br; //0c
$biga=$a1+$a2+$a3+$a4;
$bigb=$b1+$b2+$b3+$b4;
$d1=($w-$bigb)/$biga;
$_w1=$a1*$d1+$b1;
$_w2=$a2*$d1+$b2;
$_w3=$a3*$d1+$b3;
$_w4=$a4*$d1+$b4;
echo "bigb=$bigb biga=$biga _w1=$_w1  _w2=$_w2  _w3=$_w3 _w4=$_w4".$br; //0c
$_x1=$f*$x0/($f-$_w1);
$_x2=$f*$x0/($f-$_w1-$_w2);
$_x3=$f*$x0/($f-$_w1-$_w2-$_w3);
$_x4=$f*$x0/($f-$_w1-$_w2-$_w3-$_w4);
echo "_x1=$_x1   _x2=$_x2  _x3=$_x3 _x4=$_x4".$br;
$daa_all=$tg0-$tg4-($da1+$da2+$da3+$da4)-($daa1+$daa2+$daa3+$daa4)-3;
echo "daa_all=$daa_all".$br; //0c
//加热蒸汽消耗量
$d2=$_w1-$e1;
$d3=$_w2-$e2;
$d4=$_w3-$e3;
echo "d1=$d1   d2=$d2  d3=$d3 d4=$d4".$br;
//带入数据
$k1=3000;$k2=1900;$k3=1000;$k4=520;
//有效温度差分配
$drk1=$d1*$r0*1000/($k1*3600);
$drk2=$d2*$r1*1000/($k2*3600);
$drk3=$d3*$r2*1000/($k3*3600);
$drk4=$d4*$r3*1000/($k4*3600);
echo "drk1=$drk1   drk2=$drk2  drk3=$drk3 drk4=$drk4".$br;
$drk_all=$drk1+$drk2+$drk3+$drk4;//m2*k
echo "drk_all=$drk_all".$br; //0c
$at1=$drk1*$daa_all/$drk_all;
$at2=$drk2*$daa_all/$drk_all;
$at3=$drk3*$daa_all/$drk_all;
$at4=$drk4*$daa_all/$drk_all;
echo "at1=$at1   at2=$at2  at3=$at3 at4=$at4".$br;
$adt1=$drk1/$at1;
$adt2=$drk2/$at2;
$adt3=$drk3/$at3;
$adt4=$drk4/$at4;
echo "adt1=$adt1   adt2=$adt2  adt3=$adt3 adt4=$adt4".$br;
echo "根据_x查表 _x1=$_x1   _x2=$_x2  _x3=$_x3 _x4=$_x4
这里根据公式".$br;
$dap1=do_ai($_x1);
$dap2=do_ai($_x2);
$dap3=do_ai($_x3);
$dap4=do_ai($_x4);
echo "dap1=$dap1 dap2=$dap2 dap3=$dap3 dap4=$dap4".$br;
/*-------------------人工查询--------------------------------*/
$tgp4=$tg4;
$daap4=get_t($tgp4,70,60,8.1,11.0);
$tgp3=$tgp4+1+$at4+$daap4+$dap4;
echo "查表tgp3=$tgp3 _x3=$_x3".$br;
$daap3=get_t($tgp3,100,90,3.0,4.1);
$tgp2=$tgp3+1+$at3+$daap3+$dap3;
echo "查表tgp2=$tgp2 _x2=$_x2".$br;
$daap2=get_t($tgp2,120,110,1.7,2.2);
$tgp1=$tgp2+1+$at2+$daap2+$dap2;
echo "查表tgp1=$tgp1 _x1=$_x1".$br;
$daap1=get_t($tgp1,130,110,1.6,2.1);
$tgp0=$tgp1+$at1+$daap1+$dap1;
echo "tgp0=$tgp0".$br;
//$tgp3=$tgp4+1+$at4+$daap4+$dap4;
echo "tgp1=$tgp1  tgp2=$tgp2  tgp3=$tgp3 tgp4=$tgp4".$br; //0c
/*-------------------人工查询--------------------------------*/
$last_t4=$t4;
$last_t3=$tgp3+$daap3+$dap3;
$last_t2=$tgp2+$daap2+$dap2;
$last_t1=$tgp1+$daap1+$dap1;
echo "dap1=$dap1 dap2=$dap2 dap3=$dap3 dap4=$dap4".$br;
echo "daap1=$daap1 daap2=$daap2 daap3=$daap3 daap4=$daap4".$br;
$asd1=$daap1+$dap1;
$asd2=$daap2+$dap2;
$asd3=$daap3+$dap3;
$asd4=$daap4+$dap4;
echo "asd1=$asd1   asd2=$asd2  asd3=$asd3 asd4=$asd4".$br;
echo "last_t1=$last_t1 last_t2=$last_t2
last_t3=$last_t3 last_t4=$last_t4".$br;
?>