毕业设计通流部分高压压气机程序

1、#i nclude stdio.h#in elude math.h #defi ne pai 3.1415 void mai n()/*压气机参数设定*/double d1=0.7;double Mam22=0.8;double cz11=190;double c122=90;double ut仁280;double z=11;double chseff=0.91;/*压气机第一级进口轮毂比*/*压气机第一级动叶进口气流轴向速度*/*压气机出口轴向速度*/*第一级动叶进口处外径圆周速度 */*压气机级数*/*压气机等熵效率*/double Gch; double Tch111;double t

2、ch1;double ach1;double wch1;double um11;double cu11;/*空气流量*/*压气机第一级进口平均直径处圆周速度/*压气机第一级动叶进口气流绝对周向速度*/*/double Tch1;double pch1;double dech1;double arch1;double Dtch1;double Dhch1;double Dmch1;double Tsch1;double Lch1;double omch 仁0.53;double dtchW1;double w11;double c11;double cz12;/*压气机进口静温*/*压气机进口静

3、压*/*压气机进口气体密度/*压气机进口面积*/*压气机进口外径直径/*压气机进口内径直径/*压气机进口中径直径/*压气机第一级级温升/*压气机第一级叶片长度/*压气机第一级反动度*/*/*/*/*/*/*/*压气机第一级平均直径处扭速*/*压气机第一级动叶进口气流相对速度/*压气机第一级动叶进口气流绝对速度/*压气机第一级动叶出口气流轴向速度*/*/*/*压气机第一级动叶出口平均直径处圆周速度/*压气机第一级动叶出口气流绝对周向速度/*压气机第一级动叶出口绝对速度/*压气机第一级动叶出口相对速度/*压气机第一级气流绝对进气角/*压气机第一级气流绝对出气角/*压气机第一级气流相对进气角*/*/

4、*/*/*/*/*/double um12;double cu12;double c12;double w12;double garf11;double garf21;double gbt11;double gbt21;/*压气机第一级气流相对出气角 */double cz21=180;double Tch2;double pch2;double dech2;double arch2;double Dtch2;double Dhch2;double Dmch2;double Tsch2;double d2;double Lch2;double omch2=0.53;double um21;do

5、uble dtchW2;double cu21;double wj21;double cj21;double cz22;double um22;double cu22;double cj22;double wj22;double garf12;double garf22;double gbt12;double gbt22;double cz31=170;double Tch3;double pch3;double dech3;double arch3;double Dtch3;double Dhch3;double Dmch3;/*压气机第二级进口轴向速度*/*压气机第二级进口静温*/*压气机

6、第二级进口静压*/*压气机第二级进口气体密度*/*压气机第二级进口面积*/*压气机第二级进口外径直径*/*压气机第二级进口内径直径*/*压气机第二级进口中径直径*/*压气机第二级级温升*/*压气机第二级进口轮毂比*/*压气机第二级叶片长度*/*压气机第二级发动度*/*压气机第二级进口平均直径处圆周速度 */*压气机第二级平均直径处扭速 */*压气机第二级动叶进口气流绝对周向速度 */*压气机第二级动叶进口气流相对速度*/*压气机第二级动叶进口气流绝对速度*/*压气机第二级动叶出口气流轴向速度*/*压气机第二级动叶出口平均直径处圆周速度*/*压气机第二级动叶出口气流绝对周向速度 */*压气机第二

7、级动叶出口绝对速度*/*压气机第二级动叶出口相对速度*/*压气机第二级气流绝对进气角*/*压气机第二级气流绝对出气角*/*压气机第二级气流相对进气角*/*压气机第二级气流相对出气角*/*压气机第三级进口轴向速度*/*压气机第三级进口静温*/*压气机第三级进口静压*/*压气机第三级进口气体密度*/*压气机第三级进口面积*/*压气机第三级进口外径直径*/*压气机第三级进口内径直径*/*压气机第三级进口中径直径*/*压气机第三级级温升*/*压气机第三级进口轮毂比*/*压气机第三级叶片长度*/*压气机第三级发动度*/*压气机第三级进口平均直径处圆周速度 */*压气机第三级平均直径处扭速 */*压气机第

8、三级动叶进口气流绝对周向速度 */double Tsch3;double d3;double Lch3;double omch3=0.53;double um31;double dtchW3;double cu31;/*压气机第三级动叶进口气流相对速度*/*压气机第三级动叶进口气流绝对速度*/*压气机第三级动叶出口气流轴向速度*/*压气机第三级动叶出口平均直径处圆周速度*/*压气机第三级动叶出口气流绝对周向速度 */double w31;double c31;double cz32;double um32;double cu32;double c32;double w32;double gar

9、f13;double garf23;double gbt13;double gbt23;double cz4仁 160;double Tch4;double pch4;double dech4;double arch4;double Dtch4;double Dhch4;double Dmch4;double Tsch4;double d4;double Lch4;double omch4=0.53;double um41;double dtchW4;double cu41;double w41;double c41;double cz42;double um42;double cu42;do

10、uble c42;double w42;double garf14;double garf24;double gbt14;double gbt24;double cz51=150;double Tch5;double pch5;double dech5;double arch5;/*压气机第三级动叶出口绝对速度*/*压气机第三级动叶出口相对速度*/*压气机第三级气流绝对进气角*/*压气机第三级气流绝对出气角*/*压气机第三级气流相对进气角*/*压气机第三级气流相对出气角*/*压气机第四级进口轴向速度 */*压气机第四级进口静温*/*压气机第四级进口静压*/*压气机第四级进口气体密度*/*压气机

11、第四级进口面积*/*压气机第四级进口外径直径*/*压气机第四级进口内径直径 */*压气机第四级进口中径直径*/*压气机第四级级温升*/*压气机第四级进口轮毂比*/*压气机第四级叶片长度*/*压气机第四级发动度*/*压气机第四级进口平均直径处圆周速度*/*压气机第四级平均直径处扭速*/*压气机第四级动叶进口气流绝对周向速度*/*压气机第四级动叶进口气流相对速度*/*压气机第四级动叶进口气流绝对速度*/*压气机第四级动叶出口气流轴向速度*/*压气机第四级动叶出口平均直径处圆周速度*/*压气机第四级动叶出口气流绝对周向速度 */*压气机第四级动叶出口绝对速度*/*压气机第四级动叶出口相对速度*/*压

12、气机第四级气流绝对进气角*/*压气机第四级气流绝对出气角*/*压气机第四级气流相对进气角*/*压气机第四级气流相对出气角*/*压气机第五级进口轴向速度 */*压气机第五级进口静温*/*压气机第五级进口静压*/*压气机第五级进口气体密度*/*压气机第五级进口面积*/*压气机第五级进口外径直径*/*压气机第五级进口内径直径 */*压气机第五级进口中径直径*/*压气机第五级级温升*/*压气机第五级进口轮毂比 */double Dtch5;double Dhch5;double Dmch5;double Tsch5;double d5;double Lch5;double omch5=0.54;dou

13、ble um51;double dtchW5;double cu51;double w51;double c51;double cz52;double um52;double cu52;double c52;double w52;double garf15;double garf25;double gbt15;double gbt25;double cz61=140;double Tch6;double pch6;double dech6;double arch6;double Dtch6;double Dhch6;double Dmch6;double Tsch6;double d6;dou

14、ble Lch6;double omch6=0.54;double um61;double dtchW6;double cu61;/*压气机第五级叶片长度*/*压气机第五级发动度*/*压气机第五级进口平均直径处圆周速度*/*压气机第五级平均直径处扭速*/*压气机第五级动叶进口气流绝对周向速度*/*压气机第五级动叶进口气流相对速度*/*压气机第五级动叶进口气流绝对速度*/*压气机第五级动叶出口气流轴向速度*/*压气机第五级动叶出口平均直径处圆周速度*/*压气机第五级动叶出口气流绝对周向速度 */*压气机第五级动叶出口绝对速度*/*压气机第五级动叶出口相对速度*/*压气机第五级气流绝对进气角*/*

15、压气机第五级气流绝对出气角*/*压气机第五级气流相对进气角*/*压气机第五级气流相对出气角*/*压气机第六级进口轴向速度 */*压气机第六级进口静温*/*压气机第六级进口静压*/*压气机第六级进口气体密度*/*压气机第六级进口面积 */*压气机第六级进口外径直径*/*压气机第六级进口内径直径 */*压气机第六级进口中径直径*/*压气机第六级级温升*/*压气机第六级进口轮毂比*/*压气机第六级叶片长度*/*压气机第六级发动度*/*压气机第六级进口平均直径处圆周速度*/*压气机第六级平均直径处扭速*/*压气机第六级动叶进口气流绝对周向速度*/*压气机第六级动叶进口气流相对速度*/*压气机第六级动叶

16、进口气流绝对速度*/*压气机第六级动叶出口气流轴向速度*/*压气机第六级动叶出口平均直径处圆周速度*/*压气机第六级动叶出口气流绝对周向速度 */*压气机第六级动叶出口绝对速度*/*压气机第六级动叶出口相对速度*/double w61;double c61;double cz62;double um62;double cu62;double c62;double w62;/*压气机第八级进口平均直径处圆周速度*/*压气机第八级平均直径处扭速*/*压气机第八级动叶进口气流绝对周向速度*/*压气机第八级动叶进口气流相对速度*/*压气机第八级动叶进口气流绝对速度*/double um81;doubl

17、e dtchW8;double cu81;double w81;double c81;/*压气机第八级进口平均直径处圆周速度*/*压气机第八级平均直径处扭速*/*压气机第八级动叶进口气流绝对周向速度*/*压气机第八级动叶进口气流相对速度*/*压气机第八级动叶进口气流绝对速度*/double um81;double dtchW8;double cu81;double w81;double c81;double garf16;double garf26;double gbt16;double gbt26;double cz71=130;double Tch7;double pch7;double

18、dech7;double arch7;double Dtch7;double Dhch7;double Dmch7;double Tsch7;double d7;double Lch7;double omch7=0.54;double um71;double dtchW7;double cu71;double w71;double c71;double cz72;double um72;double cu72;double c72;double w72;double garf17;double garf27;double gbt17;double gbt27;double cz81=120;d

19、ouble Tch8;double pch8;double dech8;double arch8;/*压气机第六级气流绝对进气角*/*压气机第六级气流绝对出气角*/*压气机第六级气流相对进气角*/*压气机第六级气流相对出气角*/*压气机第七级进口轴向速度 */*压气机第七级进口静温*/*压气机第七级进口静压*/*压气机第七级进口气体密度*/*压气机第七级进口面积*/*压气机第七级进口外径直径*/*压气机第七级进口内径直径 */*压气机第七级进口中径直径*/*压气机第七级温升*/*压气机第七级进口轮毂比*/*压气机第七级叶片长度*/*压气机第七级发动度*/*压气机第七级进口平均直径处圆周速度*/

20、*压气机第七级平均直径处扭速*/*压气机第七级动叶进口气流绝对周向速度*/*压气机第七级动叶进口气流相对速度*/*压气机第七级动叶进口气流绝对速度*/*压气机第七级动叶出口气流轴向速度*/*压气机第七级动叶出口平均直径处圆周速度*/*压气机第七级动叶出口气流绝对周向速度 */*压气机第七级动叶出口绝对速度*/*压气机第七级动叶出口相对速度*/*压气机第七级气流绝对进气角*/*压气机第七级气流绝对出气角*/*压气机第七级气流相对进气角*/*压气机第七级气流相对出气角*/*压气机第八级进口轴向速度 */*压气机第八级进口静温*/*压气机第八级进口静压*/*压气机第八级进口气体密度*/*压气机第八级

21、进口面积*/double cz82;double um82;double cu82;double c82;double w82;double garf18;double garf28;double gbt18;double gbt28;/*压气机第八级动叶出口气流轴向速度*/*压气机第八级动叶出口平均直径处圆周速度*/*压气机第八级动叶出口气流绝对周向速度 */*压气机第八级动叶出口绝对速度*/*压气机第八级动叶出口相对速度*/*压气机第八级气流绝对进气角*/*压气机第八级气流绝对出气角*/*压气机第八级气流相对进气角*/*压气机第八级气流相对出气角*/double cz91=110;doub

22、le Tch9;double pch9;double dech9;double arch9;double Dtch9;double Dhch9;double Dmch9;double Tsch9;double d9;double Lch9;double omch9=0.55;double um91;double dtchW9;double cu91;double w91;double c91;double cz92;/*压气机第九级进口轴向速度*/*压气机第九级进口静温*/*压气机第九级进口静压 */*压气机第九级进口气体密度*/*压气机第九级进口面积*/*压气机第九级进口外径直径*/*压气机

23、第九级进口内径直径 */*压气机第九级进口中径直径*/*压气机第九级温升*/*压气机第九级进口轮毂比*/*压气机第九级叶片长度*/*压气机第九级发动度*/*压气机第九级进口平均直径处圆周速度*/*压气机第九级平均直径处扭速*/*压气机第九级动叶进口气流绝对周向速度 */*压气机第九级动叶进口气流相对速度*/*压气机第九级动叶进口气流绝对速度*/*压气机第九级动叶出口气流轴向速度*/double um92; double cu92;double c92;double w92;double garf19; double garf29;double gbt19; double gbt29;/*压气机

24、第九级动叶出口平均直径处圆周速度*/*压气机第九级动叶出口气流绝对周向速度 */*压气机第九级动叶出口绝对速度*/*压气机第九级动叶出口相对速度*/*压气机第九级气流绝对进气角*/*压气机第九级气流绝对出气角*/*压气机第九级气流相对进气角*/*压气机第九级气流相对出气角*/double cz10 仁 100;double Tch10;/*压气机第十级进口轴向速度*/*压气机第十级进口静温*/*压气机第十级进口静压*/*压气机第十级进口气体密度*/*压气机第十级进口面积*/*压气机第十级进口外径直径*/*压气机第十级进口内径直径*/double pchlO; double dechlO;dou

25、ble archlO; double DtchlO;double DhchlO; double DmchlO;double TschlO;double d10;double LchlO;double omch10=0.55; double um101;double dtchWIO; double cu101;double w101; double c101;double cz102;double um102;double cu102;double c102;double w102; double garfllO;double garf210;double gbt11O;double gbt21

26、0;double cz11 仁 90;double Tch11;double pch11;double dechll;double archil;double Dtchll;double Dhchll;double Dmchll;double Tschll;double d11;double Lchll;double omch11=0.56;/*压气机第十级进口中径直径*/*压气机第十级温升*/*压气机第十级进口轮毂比*/*压气机第十级叶片长度*/*压气机第十级发动度*/*压气机第十级进口平均直径处圆周速度*/*压气机第十级平均直径处扭速*/*压气机第十级动叶进口气流绝对周向速度*/*压气机第

27、十级动叶进口气流相对速度*/*压气机第十级动叶进口气流绝对速度*/*压气机第十级动叶出口气流轴向速度*/*压气机第十级动叶出口平均直径处圆周速度*/*压气机第十级动叶出口气流绝对周向速度 */*压气机第十级动叶出口绝对速度*/*压气机第十级动叶出口相对速度*/*压气机第十级气流绝对进气角*/*压气机第十级气流绝对出气角*/*压气机第十级气流相对进气角*/*压气机第十级气流相对出气角*/*压气机第一级进口轴向速度 */*压气机第十一级进口静温*/*压气机第十一级进口静压*/*压气机第一级进口气体密度*/*压气机第十一级进口面积 */*压气机第一级进口外径直径*/*压气机第十一级进口内径直径*/*

28、压气机第一级进口中径直径 */*压气机第十一级温升*/*压气机第一级进口轮毂比*/*压气机第十一级叶片长度*/*压气机第一级发动度*/*压气机第十一级进口平均直径处圆周速度*/*压气机第十一级平均直径处扭速*/*压气机第一级动叶进口气流绝对周向速度 */*压气机第一级动叶进口气流相对速度*/*压气机第一级动叶进口气流绝对速度*/*压气机第一级动叶出口气流轴向速度*/*压气机第十一级动叶出口平均直径处圆周速度*/double um111;double dtchWII;double cu111;double will;double c111;double cz112;double um112;do

29、uble cu112; double c112;double w112;double garflll; double garf211;double gbt111;double gbt211;/*压气机第一级动叶出口气流绝对周向速度*/*压气机第一级动叶出口绝对速度*/*压气机第一级动叶出口相对速度*/*压气机第一级气流绝对进气角*/*压气机第一级气流绝对出气角*/*压气机第十一级气流相对进气角*/*压气机第一级气流相对出气角*/double T92;double p92;double de92;double ar92;double Dt92;double Dh92;double L92;dou

30、ble chspr;double chrcoe;double hsch;double hschz;double hschav;double hsch1;double hsch2;double hsch3;double hsch4;double hsch5;double hsch6;double hsch7;double hsch8;double hsch9;double hsch10;double hsch11;/*压气机出口静温*/*压气机出口静压*/*压气机出口气体密度*/*压气机出口面积*/*压气机出口外径直径*/*压气机出口内径直径*/*压气机末级长度*/*压气机静压比*/*重热系数*

31、/*压气机通流部分静参数等熵压缩功*/*压气机总静参数等熵功*/*压气机静参数平均等熵功*/*第一级等熵功*/*第二级等熵功*/*第三级等熵功*/*第四级等熵功*/*第五级等熵功*/*第六级等熵功*/*第七级等熵功*/*第八级等熵功*/*第九级等熵功*/*第十级等熵功*/*第一级等熵功*/double n ch;double Nch;double hgch;/*风扇转速*/*风扇消耗功率*/*风扇的总耗功*/double G=653;double B=4.84;double Tt22=358.6;double Ka=1.4;double R=287.3;double pt22=112993.4

32、;double Tt3=756.9;double pt3=1242927.0;Gch=G/(B+1);Tch11 仁Tt22-cz11*cz11/(2*R*Ka/(Ka-1);ach仁sqrt(Ka*R*Tch111);wch 仁M am22*ach1;um11=ut1*sqrt(1+d1*d1)/2);cu11=um11-sqrt(wch1*wch1-cz11*cz11);c1仁 sqrt(cu11*cu11+cz11*cz11);tch仁Tt22-c11*c11/(2*R*Ka/(Ka-1);if(-0.02*tch1tch1-Tch10.02*tch1) Tch 仁Tch111;prin

33、 tf(Tch 仁.2* nnn,Tch1);elseprintf(需要重新选取);pch1=pt22*pow(Tch1/Tt22,Ka/(Ka-1);dech 仁pch1/(R*Tch1);arch仁Gch/(cz11*dech1);Dtch 仁sqrt(4*arch1/(1-d1*d1)*pai);Dhch1=Dtch1*d1;Dmch1=(Dtch1+Dhch1)/2;c1仁 cz11;T92=Tt3-c122*c122/(2*Ka*R/(Ka-1);p92=pt3*pow(T92/Tt3,Ka/(Ka-1);de92=p92/(R*T92);ar92=Gch/(cz82*de92);D

34、h92=Dhch1;Dt92=sqrt(Dhch1*Dhch1+(4*arch1)/pai);L92=(Dt92-Dh92)/2;chspr=p92/pch1;chrcoe=(z-1)/z)*(chseff*(pow(chspr,(Ka-1)/(Ka*chseff)-1)/(pow(chspr,(Ka-1)/Ka)-1)-1); hsch=Ka*R*Tch1*(pow(chspr,(Ka-1)/Ka)-1)/(Ka-1);hschz=(1+chrcoe)*hsch;hschav=hschz/z;hsch 仁0.6*hschav;hsch2=1.01*hschav;hsch3=1.10*hsch

35、av;hsch4=1.12*hschav;hsch5=1.12*hschav;hsch6=1.12*hschav;hsch7=1.11*hschav;hsch8=1.11*hschav;hsch9=1.05*hschav;hsch10=1.05*hschav;hsch11=0.9*hschav;n ch=60*ut1/(pai*Dtch1);hgch=Ka*R*(Tt3-Tt22)/(Ka-1);Nch=hgch*Gch;Tsch1=hsch1/(Ka*R*chseff/(Ka-1);Tsch2=hsch2/(Ka*R*chseff/(Ka-1);Tsch3=hsch3/(Ka*R*chsef

36、f/(Ka-1);Tsch4=hsch4/(Ka*R*chseff/(Ka-1);Tsch5=hsch5/(Ka*R*chseff/(Ka-1);Tsch6=hsch6/(Ka*R*chseff/(Ka-1);Tsch7=hsch7/(Ka*R*chseff/(Ka-1);Tsch8=hsch8/(Ka*R*chseff/(Ka-1);Tsch9=hsch9/(Ka*R*chseff/(Ka-1);Tsch10=hsch10/(Ka*R*chseff/(Ka-1);Tsch11=hsch11/(Ka*R*chseff/(Ka-1);Tch2=Tch1+Tsch1;Tch3=Tch2+Tsch2

37、;Tch4=Tch3+Tsch3;Tch5=Tch4+Tsch4;Tch6=Tch5+Tsch5;Tch7=Tch6+Tsch6;Tch8=Tch7+Tsch7;Tch9=Tch8+Tsch8;Tch10=Tch9+Tsch9;Tch11=Tch10+Tsch10;pch2=pch1*pow(Tch2/Tch1,Ka*chseff/(Ka-1);pch3=pch2*pow(Tch3/Tch2,Ka*chseff/(Ka-1);pch4=pch3*pow(Tch4/Tch3,Ka*chseff/(Ka-1);pch5=pch4*pow(Tch5/Tch4,Ka*chseff/(Ka-1);pch

38、6=pch5*pow(Tch6/Tch5,Ka*chseff/(Ka-1);pch7=pch6*pow(Tch7/Tch6,Ka*chseff/(Ka-1);pch8=pch7*pow(Tch8/Tch7,Ka*chseff/(Ka-1);pch9=pch8*pow(Tch9/Tch8,Ka*chseff/(Ka-1); pch10=pch9*pow(Tch10/Tch9,Ka*chseff/(Ka-1); pch11=pch7*pow(Tch11/Tch10,Ka*chseff/(Ka-1);dech2=pch2/(R*Tch2);dech3=pch3/(R*Tch3);dech4=pch4

39、/(R*Tch4);dech5=pch5/(R*Tch5);dech6=pch6/(R*Tch6);dech7=pch7/(R*Tch7);dech8=pch8/(R*Tch8);dech9=pch9/(R*Tch9);dech10=pch10/(R*Tch10);dech11=pch11/(R*Tch11);arch2=Gch/(cz21*dech2);arch3=Gch/(cz31*dech3);arch4=Gch/(cz41*dech4);arch5=Gch/(cz51*dech5);arch6=Gch/(cz61*dech6);arch7=Gch/(cz71*dech7);arch8=

40、Gch/(cz81*dech8);arch9=Gch/(cz91*dech9); arch10=Gch/(cz101*dech10); arch1 仁 Gch/(cz111*dech11);Dhch2=Dhch1;Dhch3=Dhch1;Dhch4=Dhch1;Dhch5=Dhch1;Dhch6=Dhch1;Dhch7=Dhch1;Dhch8=Dhch1;Dhch9=Dhch1;Dhch10=Dhch1;Dhch11=Dhch1;Dtch2=sqrt(Dhch2*Dhch2+4*arch2/pai);Dtch3=sqrt(Dhch3*Dhch3+4*arch3/pai);Dtch4=sqrt(

41、Dhch4*Dhch4+4*arch4/pai);Dtch5=sqrt(Dhch5*Dhch5+4*arch5/pai);Dtch6=sqrt(Dhch6*Dhch6+4*arch6/pai);Dtch7=sqrt(Dhch7*Dhch7+4*arch7/pai);Dtch8=sqrt(Dhch8*Dhch8+4*arch8/pai);Dtch9=sqrt(Dhch9*Dhch9+4*arch9/pai);Dtch10=sqrt(Dhch10*Dhch10+4*arch10/pai);Dtch11=sqrt(Dhch11*Dhch11+4*arch11/pai);d2=Dhch2/Dtch2;

42、d3=Dhch3/Dtch3;d4=Dhch4/Dtch4;d5=Dhch5/Dtch5;d6=Dhch6/Dtch6;d7=Dhch7/Dtch7;d8=Dhch8/Dtch8;d9=Dhch9/Dtch9;d10=Dhch10/Dtch10;d1 仁Dhch11/Dtch11;Lch 仁(Dtch1-Dhch1)/2;Lch2=(Dtch2-Dhch2)/2;Lch3=(Dtch3-Dhch3)/2;Lch4=(Dtch4-Dhch4)/2;Lch5=(Dtch5-Dhch5)/2;Lch6=(Dtch6-Dhch6)/2;Lch7=(Dtch7-Dhch7)/2;Lch8=(Dtch8-

43、Dhch8)/2;Lch9=(Dtch9-Dhch9)/2;Lch10=(Dtch10-Dhch10)/2;Lch11=(Dtch11-Dhch11)/2;Dmch2=(Dhch2+Dtch2)/2;Dmch3=(Dhch3+Dtch3)/2;Dmch4=(Dhch4+Dtch4)/2;Dmch5=(Dhch5+Dtch5)/2;Dmch6=(Dhch6+Dtch6)/2;Dmch7=(Dhch7+Dtch7)/2;Dmch8=(Dhch8+Dtch8)/2;Dmch9=(Dhch9+Dtch9)/2;Dmch10=(Dhch10+Dtch10)/2;Dmch11=(Dhch11+Dtch11

44、)/2;um2仁pai*Dmch2* nch/60; um3仁pai*Dmch3* nch/60; um41=pai*Dmch4* nch/60; um51=pai*Dmch5* nch/60; um61=pai*Dmch6* nch/60; um71=pai*Dmch7* nch/60; um81=pai*Dmch8* nch/60; um91=pai*Dmch9* nch/60; um101=pai*Dmch10* nch/60; um111=pai*Dmch11* nch/60;dtchW 仁 hsch1/(um11*chseff); dtchW2=hsch2/(um21*chseff)

45、; dtchW3=hsch3/(um31*chseff); dtchW4=hsch4/(um41*chseff); dtchW5=hsch5/(um51*chseff); dtchW6=hsch6/(um61*chseff); dtchW7=hsch7/(um71*chseff); dtchW8=hsch8/(um81*chseff); dtchW9=hsch9/(um91*chseff);dtchW10=hsch10/(um101*chseff); dtchW11=hsch11/(um111*chseff);cu11=(1-omch1)*um11-dtchW1/2;cu21=(1-omch2

46、)*um21-dtchW2/2;cu31=(1-omch3)*um31-dtchW3/2;cu41=(1-omch4)*um41-dtchW4/2;cu51=(1-omch5)*um51-dtchW5/2;cu61=(1-omch6)*um61-dtchW6/2; cu71=(1-omch7)*um71-dtchW7/2; cu81=(1-omch8)*um81-dtchW8/2; cu91=(1-omch9)*um91-dtchW9/2; cu10 仁(1-omch10)*um101-dtchW10/2; cu11 仁(1-omch11)*um111-dtchW11/2;c1仁 sqrt(c

47、u11*cu11+cz11*cz11); cj21=sqrt(cu21*cu21+cz21*cz21);c3仁 sqrt(cu31*cu31+cz31*cz31); c4仁 sqrt(cu41*cu41+cz41*cz41);c5仁 sqrt(cu51*cu51+cz51*cz51); c6仁 sqrt(cu61*cu61+cz61*cz61);c7仁 sqrt(cu71*cu71+cz71*cz71);c8仁 sqrt(cu81*cu81+cz81*cz81);c9仁 sqrt(cu91*cu91+cz91*cz91); c101=sqrt(cu101*cu101+cz101*cz101);

48、 c11仁 sqrt(cu111*cu111+cz111*cz111);w11=sqrt(um11-cu11)*(um11-cu11)+cz11*cz11);wj2仁sqrt(um21-cu21)*(um21-cu21)+cz21*cz21); w3仁Sqrt(um31-cu31)*(um31-cu31)+cz31*cz31);W41=sqrt(um41-cu41)*(um41-cu41)+cz41*cz41);W51=sqrt(um51-cu51)*(um51-cu51)+cz51*cz51);w6仁Sqrt(um61-cu61)*(um61-cu61)+cz61*cz61);W71=sqr

49、t(um71-cu71)*(um71-cu71)+cz71*cz71);W81=sqrt(um81-cu81)*(um81-cu81)+cz81*cz81);W91=sqrt(um91-cu91)*(um91-cu91)+cz91*cz91); w101=sqrt(um101-cu101)*(um101-cu101)+cz101*cz101); w11 仁sqrt(um111-cu111)*(um111-cu111)+cz111*cz111);cz12=(cz11+cz21)/2;cz22=(cz21+cz31)/2;cz32=(cz31+cz41)/2;cz42=(cz41+cz51)/2;

50、cz52=(cz51+cz61)/2;cz62=(cz61+cz71)/2;cz72=(cz71+cz81)/2;cz82=(cz81+c91)/2;cz92=(cz91+c101)/2;cz102=(cz101+c111)/2;cz112=(cz111+c111)/2;um12=um11;um22=um21;um32=um31;um42=um41;um52=um51;um62=um61;um72=um71;um82=um81;um92=um91; um102=um101;um112=um111;cu12=(1-omch1)*um12+dtchW1/2;cu22=(1-omch2)*um22+

51、dtchW2/2;cu32=(1-omch3)*um32+dtchW3/2;cu42=(1-omch4)*um42+dtchW4/2;cu52=(1-omch5)*um52+dtchW5/2;cu62=(1-omch6)*um62+dtchW6/2;cu72=(1-omch7)*um72+dtchW7/2;cu82=(1-omch8)*um82+dtchW8/2;cu92=(1-omch9)*um92+dtchW9/2; cu102=(1-omch10)*um102+dtchW10/2; cu112=(1-omch11)*um112+dtchW11/2;c12=sqrt(cu12*cu12+c

52、z12*cz12);cj22=sqrt(cu22*cu22+cz22*cz22);c32=sqrt(cu32*cu32+cz32*cz32);c42=sqrt(cu42*cu42+cz42*cz42);c52=sqrt(cu52*cu52+cz52*cz52);c62=sqrt(cu62*cu62+cz62*cz62);c72=sqrt(cu72*cu72+cz72*cz72);c82=sqrt(cu82*cu82+cz82*cz82);c92=sqrt(cu92*cu92+cz92*cz92);c102=sqrt(cu102*cu102+cz102*cz102);c112=sqrt(cu112*cu112+cz112*cz112);W12=sqrt(um12-cu12)*(um12-cu12)+cz12*cz12);Wj22=sqrt(um22-cu22)*(um22-cu22)+cz22*cz22);W32=sqrt(um32-cu32)*(um32-cu32)+cz32*cz32);W42=sqrt(um42-cu4