Можете подсказать что не так.Я недавно начал работать в матлабе,и еще не все знаю

@Andreii · Регистрация: 04.04.2015

Author24 — интернет-сервис помощи студентам

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
Gapa=4900/3600; %Debitul de apa,[kg/s]
ti=17; %Temperatura de intrare a apei,[grade Celsius]
te=74; %Temperatura de iesire a apei,[grade Celsius]
tv=100;%Temperatura vaporilor,[grade Celsius]
ta=45.5;%temperatura apei,[grade Celsius]
u=0.04; %Umiditatea
r=2260*10^3; %Caldura latenta a aburului,[J/kg]
cp_apa=4180.2; %Caldura specifica a apei,[J/(kg*K)]
Qu=Gapa*cp_apa*(te-ti);%Caldura utila pentru incalzirea apei,[W]
Qp=0.07*Qu; %Pierderi de caldura,[W]
Gabur=(Qp+Gapa*cp_apa*(te-ti))/(r*(1-u)); %Debitul de abur,[kg/s]
Qc=Gabur*r*(1-u); %Debitul de caldura consumat, [W]
k=1100; %coeficient total de transfer termic propus,[W/(m2*K)]
Tm=(tv-ti-(tv-te))/log((tv-ti)/(tv-te));%Temperatura medie logaritmica 
                                             % a apei,[grade Celsius]
Apres=Qc/(k*Tm);%Aria presupusa,[m2]  
disp('Apres');
disp(Apres);
ASTAS=6.33; %Aria STAS,[m2]
Di=0.2;     %Diametrul interior al mantalei,[m]
s=0.0045;    %Grosimea mantalei,[m]
beta=0.002; %Grosimea tevilor,[m]
de=0.02;    %Diametrul exterior al tevilor,[m]
di=0.016;    %Diametrul interior al tevilor,[m]
n=28;        %Numarul de tevi
Nt=4;        %Numarul de treceri
l=4;         %Lungimea unei tevi,[m]
dechiv=0.065;%Diametrul echivalent,[m]
dm=(de-di)/log(de/di);%Diametrul mediu logaritmic al tevilor,[m]
 
lambdaOL=46.5; %Conductivitatea termica a otelului,[W/(m*K)]
% Calculam coeficientul de transfer termic in tevi alfa apa
ro_apa=990;%Densitatea apei,[kg/m3]
Gvapa=Gapa/ro_apa;%[kg/s]
w=4*Gvapa*Nt/(pi*di^2*n);%Viteza de curgere a apei prin tevi,[m/s]
disp('Vieza');
disp(w);
ita_apa=59.082*10^-5;%Viscozitatea cinematica a apei,[kg/(m*s)]
lambda_apa=0.64;%Conductivitatea termica a apei,[W/(m*K)]
 Re=w*di*ro_apa/ita_apa;
 disp('Reynolds');
 disp(Re);
 Pr=cp_apa*ita_apa/lambda_apa;
 
   g=9.81; %Aceleratia gravitationala, [m/s]
    
% Determinam coeficientul de transfer termic dintre tevi
 
    
tol=0.01;
   
for tp1=0:0.0001:tv
    
    lambda_cr=1.5;    %Conductivitatea termica a crustei
    r2=0.00053;   %Rezistenta termica a crustelor(apa)
    r1=0.00009;   %Rezistenta termica a crustelor(vapori)
    beta_1= lambda_cr * r1;   %Grosimea de pe exteriorul tevii [m]
    beta_2= lambda_cr *r2;    %Grosimea din interiorul tevii [m]
    de_cr1 = de+2*beta_1;        
    di_cr2 = di-2*beta_2;
    dm_cr1 = (de_cr1 - de)/ log( de_cr1/de);   %Diametrul mediu logaritmic al crustei1 [m]
    dm_cr2 = (di - di_cr2)/ log( di/di_cr2);    %Diametrul mediu logaritmic al crustei2 [m]
    
     
    alfa_v=0.943*((r*ro_apa^2*lambda_apa^3*g)/(ita_apa*de*(tv-tp1))^(1/4));
    k2=(alfa_v*(tv-tp1))/(tv - ta);
 
    tp2=tp1 - (alfa_v*de*(tv - tp1)*beta_1)/(lambda_cr*dm_cr1);
    tp3=tp2 - (alfa_v*de*(tv - tp1)*beta_1)/(lambdaOL*dm);
    tp4=tp3 - (alfa_v*de*(tv - tp1)*beta_2)/(lambda_cr*dm_cr2);
    
    ita_p=0.0015*exp(-0.018*tp2);
    
if Re<2300
    B=Re*Pr*di/l;
         Nu = (3.65+0.0668*B/(1+0.045*B^0.66))*(ita_apa/ita_p)^0.14;
    if 2300<Re<7000
         Nu = 0.008*Re^0.9*Pr^0.43;
    else Re>10^4
               Nu = 0.023*Re^0.8*Pr^0.4;
     
           disp('Nuselt');
           disp(Nu);  
    end
end
 
    alfa_a = Nu*lambda_apa/di;
 
    k1=(alfa_a*di*(tp2-ta))/(de*(tv - ta));
 
 
    if abs(k1 - k2)< tol
        disp( ' Valorile obtinute sunt:' );
        disp( ' tp1     tp4    k1     k2  ');
        disp([tp1,tp4,k1,k2]);
        break
   
    end
end
 
 
kcalc=1/((1/alfa_v)*(de/de_cr1)+(r1*(de/dm_cr1))+(beta/lambdaOL)*(de/dm)+(r2*(de/dm_cr2))+(1/alfa_a)*(de/di_cr2));
disp('Kcalc');
disp(kcalc);
Acalc=Qc/(kcalc*Tm);
E=(ASTAS-Acalc)/ASTAS*100;
disp('Acalc         E       ');
disp([Acalc,E]);
disp([Acalc,E]);

@Norwall · 10.05.2015, 11:35

Andreii, немного поправил, считает:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
Gapa=4900/3600; %Debitul de apa,[kg/s]
ti=17; %Temperatura de intrare a apei,[grade Celsius]
te=74; %Temperatura de iesire a apei,[grade Celsius]
tv=100;%Temperatura vaporilor,[grade Celsius]
ta=45.5;%temperatura apei,[grade Celsius]
u=0.04; %Umiditatea
r=2260*10^3; %Caldura latenta a aburului,[J/kg]
cp_apa=4180.2; %Caldura specifica a apei,[J/(kg*K)]
Qu=Gapa*cp_apa*(te-ti);%Caldura utila pentru incalzirea apei,[W]
Qp=0.07*Qu; %Pierderi de caldura,[W]
Gabur=(Qp+Gapa*cp_apa*(te-ti))/(r*(1-u)); %Debitul de abur,[kg/s]
Qc=Gabur*r*(1-u); %Debitul de caldura consumat, [W]
k=1100; %coeficient total de transfer termic propus,[W/(m2*K)]
Tm=(tv-ti-(tv-te))/log((tv-ti)/(tv-te));%Temperatura medie logaritmica
% a apei,[grade Celsius]
Apres=Qc/(k*Tm);%Aria presupusa,[m2]
disp('Apres');
disp(Apres);
ASTAS=6.33; %Aria STAS,[m2]
Di=0.2;     %Diametrul interior al mantalei,[m]
s=0.0045;    %Grosimea mantalei,[m]
beta=0.002; %Grosimea tevilor,[m]
de=0.02;    %Diametrul exterior al tevilor,[m]
di=0.016;    %Diametrul interior al tevilor,[m]
n=28;        %Numarul de tevi
Nt=4;        %Numarul de treceri
l=4;         %Lungimea unei tevi,[m]
dechiv=0.065;%Diametrul echivalent,[m]
dm=(de-di)/log(de/di);%Diametrul mediu logaritmic al tevilor,[m]
 
lambdaOL=46.5; %Conductivitatea termica a otelului,[W/(m*K)]
% Calculam coeficientul de transfer termic in tevi alfa apa
ro_apa=990;%Densitatea apei,[kg/m3]
Gvapa=Gapa/ro_apa;%[kg/s]
w=4*Gvapa*Nt/(pi*di^2*n);%Viteza de curgere a apei prin tevi,[m/s]
disp('Vieza');
disp(w);
ita_apa=59.082*10^-5;%Viscozitatea cinematica a apei,[kg/(m*s)]
lambda_apa=0.64;%Conductivitatea termica a apei,[W/(m*K)]
Re=w*di*ro_apa/ita_apa;
disp('Reynolds');
disp(Re);
Pr=cp_apa*ita_apa/lambda_apa;
 
g=9.81; %Aceleratia gravitationala, [m/s]
 
% Determinam coeficientul de transfer termic dintre tevi
 
 
tol=0.01;
 
for tp1=0:0.0001:tv
    
    lambda_cr=1.5;    %Conductivitatea termica a crustei
    r2=0.00053;   %Rezistenta termica a crustelor(apa)
    r1=0.00009;   %Rezistenta termica a crustelor(vapori)
    beta_1= lambda_cr * r1;   %Grosimea de pe exteriorul tevii [m]
    beta_2= lambda_cr *r2;    %Grosimea din interiorul tevii [m]
    de_cr1 = de+2*beta_1;
    di_cr2 = di-2*beta_2;
    dm_cr1 = (de_cr1 - de)/ log( de_cr1/de);   %Diametrul mediu logaritmic al crustei1 [m]
    dm_cr2 = (di - di_cr2)/ log( di/di_cr2);    %Diametrul mediu logaritmic al crustei2 [m]
    
    
    alfa_v=0.943*((r*ro_apa^2*lambda_apa^3*g)/(ita_apa*de*(tv-tp1))^(1/4));
    k2=(alfa_v*(tv-tp1))/(tv - ta);
    
    tp2=tp1 - (alfa_v*de*(tv - tp1)*beta_1)/(lambda_cr*dm_cr1);
    tp3=tp2 - (alfa_v*de*(tv - tp1)*beta_1)/(lambdaOL*dm);
    tp4=tp3 - (alfa_v*de*(tv - tp1)*beta_2)/(lambda_cr*dm_cr2);
    
    ita_p=0.0015*exp(-0.018*tp2);
    
    if Re<2300
        B=Re*Pr*di/l;
        Nu = (3.65+0.0668*B/(1+0.045*B^0.66))*(ita_apa/ita_p)^0.14;
    elseif Re > 2300 && Re<7000
        Nu = 0.008*Re^0.9*Pr^0.43;
    else
        Nu = 0.023*Re^0.8*Pr^0.4;
        
        disp('Nuselt');
        disp(Nu);
    end
end
 
alfa_a = Nu*lambda_apa/di;
 
k1=(alfa_a*di*(tp2-ta))/(de*(tv - ta));
 
 
if abs(k1 - k2)< tol
    disp( ' Valorile obtinute sunt:' );
    disp( ' tp1     tp4    k1     k2  ');
    disp([tp1,tp4,k1,k2]);
    break
    
end
 
kcalc=1/((1/alfa_v)*(de/de_cr1)+(r1*(de/dm_cr1))+(beta/lambdaOL)*...
(de/dm)+(r2*(de/dm_cr2))+(1/alfa_a)*(de/di_cr2));
disp('Kcalc');
disp(kcalc);
Acalc=Qc/(kcalc*Tm);
E=(ASTAS-Acalc)/ASTAS*100;
disp('Acalc         E       ');
disp([Acalc,E]);
disp([Acalc,E]);

Ошибка была здесь:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
if Re<2300
    B=Re*Pr*di/l;
         Nu = (3.65+0.0668*B/(1+0.045*B^0.66))*(ita_apa/ita_p)^0.14;
    if 2300<Re<7000
         Nu = 0.008*Re^0.9*Pr^0.43;
    else Re>10^4
               Nu = 0.023*Re^0.8*Pr^0.4;
     
           disp('Nuselt');
           disp(Nu);  
    end
end

Только вывод на экран числа Нуссельта замедляет вычисления:

Matlab M
1
2
disp('Nuselt');
disp(Nu);

@Andreii · 10.05.2015, 12:03 **[ТС]**

спасибо большое

,сейчас все исправлю

Новые блоги и статьи Все статьи Все блоги /
sshd restrictions, ssh access limitations jigi33 26.03.2025 sshd restrictions \| ssh access limitations рестрикции доступа на сервер sshd статья: https:/ / www. golinuxcloud. com/ restrict-allow-ssh-certain-users-groups-rhel	Компиляция C++ с Clang API NullReferenced 24.03.2025 Компиляторы обычно воспринимаются как черные ящики, которые превращают исходный код в исполняемые файлы. Мы запускаем компилятор командой в терминале, и вуаля — получаем бинарник. Но что если нужно. . .	Многопоточное программирование в C#: Класс Thread UnmanagedCoder 24.03.2025 Когда запускается приложение на компьютере, операционная система создаёт для него процесс - виртуальное адресное пространство. В C# этот процесс изначально получает один поток выполнения — главный. . .	SwiftUI Data Flow: Передача данных между представлениями mobDevWorks 23.03.2025 При первом знакомстве со SwiftUI кажется, что фреймворк предлагает избыточное количество механизмов для передачи данных: @State, @Binding, @StateObject, @ObservedObject, @EnvironmentObject и другие. . . .	Моки в Java: Сравниваем Mockito, EasyMock, JMockit Javaican 23.03.2025 Как протестировать класс, который зависит от других сложных компонентов, таких как базы данных, веб-сервисы или другие классы, с которыми и так непросто работать в тестовом окружении? Для этого и. . .
Архитектурные паттерны микросервисов: ТОП-10 шаблонов ArchitectMsa 22.03.2025 Популярность микросервисной архитектуры объясняется множеством важных преимуществ. К примеру, она позволяет командам разработчиков работать независимо друг от друга, используя различные технологии и. . .	Оптимизация рендеринга в Unity: Сортировка миллиона спрайтов GameUnited 22.03.2025 Помните, когда наличие сотни спрайтов в игре приводило к существенному падению производительности? Время таких ограничений уходит в прошлое. Сегодня геймдев сталкивается с задачами совершенно иного. . .	Образование и практика Igor3D 21.03.2025 Добрый день А вот каково качество/ эффективность ВУЗовского образования? Аналитическая геометрия изучается в первом семестре и считается довольно легким курсом, что вполне справедливо. Ну хорошо,. . .	Lazarus. Таблица с объединением ячеек. Massaraksh7 21.03.2025 Понадобилась представление на экране таблицы с объединёнными ячейками. И не одной, а штук триста, и все разные. На Delphi я использовал для этих целей TStringGrid, и то, кривовато получалось. А в. . .	Async/await в Swift: Асинхронное программирование в iOS mobDevWorks 20.03.2025 Асинхронное программирование долго было одной из самых сложных задач для разработчиков iOS. В течение многих лет мы сражались с замыканиями, диспетчеризацией очередей и обратными вызовами, чтобы. . .

@Andreii 0 / 0 / 0 Регистрация: 04.04.2015 Сообщений: 6
	10.05.2015, 12:03 [ТС]
	спасибо большое ,сейчас все исправлю 0

Опции темы