求大牛帮小弟我分析一下。第二次发帖了，这个有关问题很急多谢大家

求大牛帮我分析一下。。第二次发帖了，这个问题很急，谢谢大家

m_fQj=m_fQj-50;

  for(int i=0;i<100;i++)/*喷淋水*/

  {

    m_ft1=m_ft1-5;

for(int j=0;j<10;j++)

{

      m_ft2=m_ft2-5;

  for(int k=0;k<10;k++)

  {

       m_fuof=m_fuof-2;

for(int h=0;h<40;h++)

{

           m_fts1=m_fts1-2;

 for(int zz=0;zz<4;zz++)

 {	

if(m_jiezhi2==0)


{

     CString strSQL;

 int a=(m_fts1+0.5);

 float b=a;

strSQL.Format(_T("select * from 水的物性参数 where 温度=%f"),b);

pRecordsethb.CreateInstance((_uuidof(Recordset)));

try

{

pRecordsethb->Open((_bstr_t)strSQL,connh2.GetInterfacePtr(),adOpenDynamic,adLockOptimistic,adCmdText);

}

catch(_com_error& error)

{

CString erroMsg;

erroMsg.Format("error:%s",(LPTSTR)error.Description());

AfxMessageBox(erroMsg);	

}

_variant_t vt5; 

_variant_t vt6;

_variant_t vt7;

_variant_t vt8;

vt5=pRecordsethb->GetCollect("密度");

vt6=pRecordsethb->GetCollect("比热容");

vt7=pRecordsethb->GetCollect("动力黏度");

vt8=pRecordsethb->GetCollect("导热系数");

CString str1;

str1=(LPCSTR)_bstr_t(vt5);

pw1=atof(str1);

str1=(LPCSTR)_bstr_t(vt6);

cpw1=atof(str1);

str1=(LPCSTR)_bstr_t(vt7);

uw1=atof(str1);

str1=(LPCSTR)_bstr_t(vt8);

rw1=atof(str1);

    

}

//*************水膜温度ts下的饱和蒸汽压************//

int ff=1;

float m11,m12,m13,Lt,qms,qvs,Rsw,z,dm,qma;

z=log(m_fdo/m_fdi);

dm=(m_fdo-m_fdi)/z;

Rsw=m_fdo*(m_fdo-m_fdi)/2/dm/m_fDaoxi;

float m1,m2,m3,m4,A,AS,W,Bp,us,A1,A2,A3,m41,m42,m43,ma1,ma2,ma3,Ka1,Ka2,Ka3;

 CString stra,strb,strc;

 int nn ,nn1;          /*每排管数*/

 m11=m_fdo/m_fdi/ai1;

 m12=m_fdo/m_fdi/ai2;

 m13=m_fdo/m_fdi/ai3;

 nn=m_iN/m_nb+0.99999;

 Bp=(nn-1)*m_fBt*1000+m_fdo*1000;

 nn1=(m_iN+1)/2;

     us=2*nn1*m_fQj*1000/pw1/Bp/3600;

     m1=pow(us,0.56);

 m2=pow(m_fuof,0.52);

 m3=pow(m_ft1,1.38);

 m4=100.6*m1*m2*m3;/*拟合公式*/


     m41=m11+Rsw+m_fRi+m_fRo+1/m4;

 m42=m12+Rsw+m_fRi+m_fRo+1/m4;

 m43=m13+Rsw+m_fRi+m_fRo+1/m4;

 A1=Q1*m41*1000/tm1;

 A2=Q2*m42*1000/tm2;

 A3=Q3*m43*1000/tm3;

 A=A1+A2+A3;

     if (j==0&&i==0&&k==0&&h==0&&zz==0)

     {

 AA=A;

     }

     if (AA<=A)

 AA=AA;

 else AA=A;

  K1=1/m41;

  K2=1/m42;

  K3=1/m43;

  ma1=m11+Rsw+1/m4;

  ma2=m12+Rsw+1/m4;

  ma3=m13+Rsw+1/m4;

  Ka1=1/ma1;

  Ka2=1/ma2;

  Ka3=1/ma3;

  float RY=(K1+K2+K3)/3;

  if (KK<=RY)

  KK=(K1+K2+K3)/3;

  else

    KK=KK;

 

  Lt=A/m_iN/3.1415926/m_fdo/m_inp;

  if (Lt>=1.2&&Lt<=1.56) Lt=1.56;

 else if (Lt>1.56&&Lt<=2.475) Lt=2.475;

 else if (Lt>2.475&&Lt<=3.39) Lt=3.39;

 else if (Lt>3.39&&Lt<=5.22) Lt=5.22;