多效蒸发课程设计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='<hr>';
//進料
$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='一';
include('fs.php');
$word='二';
include('fs.php');
$word='三';
include('fs.php');
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;
?>
?
?
總結
以上是生活随笔為你收集整理的多效蒸发课程设计php程序的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: editor http://www.ev
- 下一篇: PHPStorm 502 / 404