clc;

clear;

D=4;

D0 = 0.012;

qm = 0.06;

I = 380;

G = 810;

Tfi = 323;

Ta = 285;

V = 4.2;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

if D>70

D=70;

else

D=D;

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

M = 0.85;%透过率τ

N = 0.9; %吸收率α

L = 2.3*10^(-3); %厚度L

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

K = M*N/(1-(1-M)*0.16)+0.27*(1-(exp(1))^(-4*L))  ; %2-8A

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

A = 5.67*10^(-8); %σ

B = 0.1; %ε_p吸热板发射率

C = 0.9; %ε_c盖板发射率

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%假设Tp

Tpp=70+273;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

hw=5.7+3.8*V; %

f=(1+0.0892*hw-0.1166*hw*B)*(1+0.07866); %

c=520*(1-0.000051*D^2); %用D表示倾斜角度β

Z = 0.43*(1-Tpp);    %e=0.43(1-100?T_p )%

Ut=(1/(c/Tpp*((Tpp-Ta)/(1+f))^Z)+1/hw)^(-1)+A*(Tpp+Ta)*(Tpp^2+Ta^2)/((B+0.00591*hw)^(-1)+(1+f+0.133*B)/C-1);

%公式2-39       A=5.67*10^(-8)     B=0.1   %吸热版发射功率          C=    %盖板发射率

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

E= 0.04       ;% λ_(p-b),λ_(p-e)隔热层导热系数

H= 25*10^(-3) ;%L_(p-b)底部厚度

J= 25*10^(-3) ;% L_(p-e) 侧壁厚度

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Ub=1/(H/E+1/(5.7+3.8*V))  ;

%2-43a

Ue=1/(J/E+1/(5.7+3.8*V))  ;

%2-43b  E= 0.04  ;%隔热层导热系数        H= 25*10^(-3) ;%地步厚度      J= 25*10^(-3) ;%侧壁厚度

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Ae = 0.5;      %侧壁面积  A_e

At = 1.82;     %采光面积  A_t

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Ul=Ut+Ub+Ae/At*Ue  ; %2-14A_t采光面积;A_e侧面面积;V环境风速

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

tp=70;

ta=12;

Cp=4.185*10^3 ;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Ql = At*Ul*(tp-ta);   %2-3 Ql

Qu = At*G*K-At*Ul*(tp-ta);   %2-4  Qu

Tf = Qu/(2*Cp*qm)+Tfi ; %2-61c

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

W = 135*10^(-3);  %翅片宽度  W

Di = 10*10^(-3);   %内径D_i

O=1000;%水的密度

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Ww =qm/((pi*Di^2/4)*O*floor(1/W))  ; %2-51

tf=Tf-273;

hi =(1430+23.3*tf^2)*Ww^0.8*Di^(-0.2);   %2-50

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

W = 135*10^(-3);  %翅片宽度  W

I = 380         ;   %导热系数λ_p

J = 0.35*10^(-3) ;  %翅片厚度  δ

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

m=sqrt(Ul/(I*J))     ;             %2-49  I  J

F=tanh(m*(W-D0)/2)/(m*(W-D0)/2) ;  %2-48

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

a=2.5*10^(-3);%焊缝宽度

b=0.3*10^(-3);%焊缝高度

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Cb=I*a/b;

Ff = 1/Ul/(W*(1/(Ul*((W-D0)*F+D0)))+1/Cb+1/(hi*Di*pi)); %2-46

Fr = qm*Cp/(At*Ul)*(1-exp(-At*Ul*Ff/(qm*Cp)));  %2-54

Tp = Tfi+(Qu/(At*Fr*Ul))*(1-Fr) ;   %

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%if abs(Tp-Tp')<10^(-5)I=Fr*(K-Ul*(Tfi-Ta)/G)  %η=F_R [(τα)_e-U_L  (T_fi-T_a)/G]elseTpp = Tp;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

if  abs(Tp-Tpp)<10^(-5)

S = Fr*(K-Ul*(Tfi-Ta)/G) ;

else

while  abs(Tp-Tpp)>10^(-5)

hw = 5.7+3.8*V; %

f=(1+0.0892*hw-0.1166*hw*B)*(1+0.07866); %

c=520*(1-0.000051*D^2); %用D表示倾斜角度β

Z = 0.43*(1-Tp);    %e=0.43(1-100?T_p )%

Ut=(1/(c/Tp*((Tp-Ta)/(1+f))^Z)+1/hw)^(-1)+A*(Tp+Ta)*(Tp^2+Ta^2)/((B+0.00591*hw)^(-1)+(1+f+0.133*B)/C-1);

%公式2-39       A=5.67*10^(-8)     B=0.1   %吸热版发射功率          C=    %盖板发射率

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

E= 0.04       ;% λ_(p-b),λ_(p-e)隔热层导热系数

H= 25*10^(-3) ;%L_(p-b)底部厚度

J= 25*10^(-3) ;% L_(p-e) 侧壁厚度

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Ub=1/(H/E+1/(5.7+3.8*V))  ;

%2-43a

Ue=1/(J/E+1/(5.7+3.8*V))  ;

%2-43b  E= 0.04  ;%隔热层导热系数        H= 25*10^(-3) ;%地步厚度      J= 25*10^(-3) ;%侧壁厚度

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Ae = 0.5;      %侧壁面积  A_e

At = 1.82;     %采光面积  A_t

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Ul=Ut+Ub+Ae/At*Ue  ; %2-14A_t采光面积;A_e侧面面积;V环境风速

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

tp=70;

ta=12;

Cp=4.185*10^3 ;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Ql = At*Ul*(tp-ta);   %2-3 Ql

Qu = At*G*K-At*Ul*(tp-ta);   %2-4  Qu

Tf = Qu/(2*Cp*qm)+Tfi ; %2-61c

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

W = 135*10^(-3);  %翅片宽度  W

Di = 10*10^(-3);   %内径D_i

O=1000;%水的密度

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Ww =qm/((pi*Di^2/4)*O*floor(1/W))  ; %2-51

tf=Tf-273;

hi =(1430+23.3*tf^2)*Ww^0.8*Di^(-0.2);   %2-50

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

W = 135*10^(-3);  %翅片宽度  W

I = 380         ;   %导热系数λ_p

J = 0.35*10^(-3) ;  %翅片厚度  δ

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

m=sqrt(Ul/(I*J))     ;             %2-49  I  J

F=tanh(m*(W-D0)/2)/(m*(W-D0)/2) ;  %2-48

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

a=2.5*10^(-3);%焊缝宽度

b=0.3*10^(-3);%焊缝高度

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Cb=I*a/b;

Ff = 1/Ul/(W*(1/(Ul*((W-D0)*F+D0)))+1/Cb+1/(hi*Di*pi)); %2-46

Fr = qm*Cp/(At*Ul)*(1-exp(-At*Ul*Ff/(qm*Cp)));  %2-54

Tp = Tfi+(Qu/(At*Fr*Ul))*(1-Fr) ;   %

end

end

S = Fr*(K-Ul*(Tfi-Ta)/G)

matlab怎么实现循环,matlab怎么实现直到型循环相关推荐

  1. python循环语句-python语句中Python循环语句(附代码)

    python语句多如牛毛,对于很多初学者来说,不知道该如何下手.今天本文将着重讲述python语句中for语句和while语句.都知道这2种语句都属于循环语句,for语句属于遍历循环,while语句属 ...

  2. c语言第六章条件型循环结构,C语言课件(第六章 循环结构)

    C语言课件 第六章 循环结构 在许多复杂的问题中,常常需要做大量的类似的计算处理,这些类似的相同的处理过程存在重复性过程,即在某种条件成立的情况下需要处理相同的操作时,就要使用循环结构来完成. 类似于 ...

  3. 分类预测 | MATLAB实现WOA-CNN-GRU鲸鱼算法优化卷积门控循环单元数据分类预测

    分类预测 | MATLAB实现WOA-CNN-GRU鲸鱼算法优化卷积门控循环单元数据分类预测 分类效果 基本描述 1.Matlab实现WOA-CNN-GRU多特征分类预测,多特征输入模型,运行环境Ma ...

  4. 数学建模与MATLAB计算之007.MATLAB中的循环语句

    007.MATLAB中的循环语句 文章目录 007.MATLAB中的循环语句 一.for循环 (一)格式1: (二)格式2: 二.while循环 三.break语句和continue语句 四.循环的嵌 ...

  5. matlab程序循环,matlab循环程序只得到一个结果

    如题.matlab编程,for循环语句,本来要它循环100次,但是得到的结果却只有一个,本想得到100个的. 求牛人指点.整个程序较简单,很多类似语句,目的是求7种方法得到的列向量的7个离差平方和,这 ...

  6. do matlab的 while循环_PHP while和do while循环

    循环语句可以在满足条件的情况下反复执行某一个操作.在 PHP 中提供了 4 个循环控制语句,分别是 while.do while.for 和 foreach 循环语句.本节我们主要介绍一下 while ...

  7. mysql matlab for循环嵌套循环_matlab里面的for循环嵌套

    展开全部 函数思路:先对1到e5a48de588b632313133353236313431303231363533313334313732331000的数K进行循环,从2到K-1进行j循环,如果K能 ...

  8. matlab中step_使用MATLAB中bertool进行通信仿真

    MATLAB中的Communication System Toolbox中的bertool图像化工具bertool提供了极为方便的通信仿真环境以测试仿真模型的误比特率性能,并且可以很好地和理论值进行对 ...

  9. Matlab学习一本通,matlab基础教程

    链接:https://pan.baidu.com/s/1uTCbiRfIxcrt6lmiy6_QlQ  提取码:f2dn  Matlab学习一本通,matlab基础教程 <MATLAB R201 ...

  10. matlab dsp实验报告,matlab实验报告14.pdf

    matlab实验报告14 Matlab实验报告 院系名称 :信息科学与工程学院 专业班级 :通信工程 1303 指导老师 :陈科文 ,支国明 ,张金焕 ,周扬 学生姓名 : 学号 : 目录 实验一 熟 ...

最新文章

  1. linux怎么添加ubuntu源,ubuntu/linuxmint如何添加和删除PPA源
  2. (二)nimlang web开发 hello world jester
  3. Flask-模板渲染
  4. 手把手实现YOLOv3(三)
  5. 搜索框中“请输入搜索keyword”
  6. SQL SERVER数据库备份时出现“操作系统错误5(拒绝访问)。BACKUP DATABASE 正在异常终止。”错误的解决办法...
  7. 大数据与Hadoop有什么关系?大数据Hadoop入门简介
  8. 3.企业应用架构模式 --- 映射到关系数据库
  9. 关于人生的10个真相, 希望你越早知道越好
  10. 访问服务器根目录文件,可道云 访问服务器根目录
  11. C#数组和集合类详解
  12. 寻星时卫星数字电视接收机的信号检测功能
  13. 【毕业设计】STM32智能药箱系统 - 单片机 嵌入式 物联网
  14. IDEA使用Git大全
  15. CentOS7图形界面启动报错unable to connect to X server
  16. hp服务器装系统键盘没反应,u盘装系统f12没反应怎么回事|u盘装系统按f12没反应怎么解决...
  17. [wine5.0] 解决wine所有软件乱码(中文不显示)问题
  18. 一点资讯signature分析
  19. 京东再次传来裁员消息,IT行业又有什么波动呢!
  20. omv检查硬盘坏道_技术员们有福了,关于怎么用MHDD修复硬盘坏道,现在免费给大家了!!!...

热门文章

  1. TFS2008 基本安装
  2. [Vuex系列] - Vuex中的getter的用法
  3. django模板中使用JQ代码实现瀑布流显示效果
  4. WebStorm错误--无法显示文件夹目录
  5. 判断点是否在一个任意多边形中
  6. TableCellRenderer TableCellEditor(三)
  7. 优化tableView性能(针对滑动时出现卡的现象)
  8. 中国女足在2022女足亚洲杯决赛中击败韩国队,时隔16年再夺亚洲杯冠军。
  9. 不停歇的 Java 即将发布 JDK 16,新特性速览!
  10. Visual C++语言编程开发详解(孙鑫老师)