Реферат диссертация 152 стр, 75 рис, 19 табл, 116 библ. Ключевые слова. Тепловой анализ, параметрическая генерация цепных моделей, электрические схемы замещения, неоднородность магнитного поля, распределительные трансформаторы, системы инженерного анализа cae системы.
 Скачать 6.09 Mb.
|
|
Solidworks. Компьютерное моделирование в инженерной практике / А.А. Алямовский и др. – СПб.: БХВ-Петербург, 2005 – 800 с. 139 Иностранные источники 105. Arslan, S. Investigation of current density, magnetic flux density and ohmic losses for single-veined, litz and foil structured conductors at different fre- quences / S. Arslan, İ. Tarimer, M. Emin Güven // Pamukkale University Jour- nal of Engineering Sciences, Vol.19, 2013. pp.195-200. 106. Different Approach to Thermal Modeling of Transformers - a comparison of methods / V. Madžarević [and others] // International Journal of Energy and Environment, Issue 5, Volume 5, 2011 107. Ghareh, M. Thermal Modeling of Dry-Transformers and Estimating Tem- perature Rise / M. Ghareh, L. Sepahi // World Academy of Science, Engineer- ing and Technology, 21, 2008 108. Holtshausen, C.B Transformer thermal modeling, load curve development and life estimation / R & D Journal of the South Africa institution of mechani- cal engineering, 2015, Volume 31, pp. 12-16. 109. Hot spot studies for sheet wound transformers windings / S. Kennedy [and others] // Doble Engineering Company – 80th International Conference of Doble Clients, 2013 110. Mamizadech, A. Analyzing and Comparing the Hot-spot Thermal Models of HV/LV Prefabricated and Outdoor Oil-Immersed Power Transformers / A. Mamizadech, I. Iskander // International Journal of Electronics and Electric- al Engineering 6, 2012 111. Mullineux, N. Current distribution in sheet- and foil-wound transformers / N. Mullineux, J.R. Reed, I.J. Whyte // The Proceedings of the Institution of Electrical Engineers,1969. Vol. 116, Issue 1, p. 127 – 129. 112. Nigam, M. Multilayer Barrel-wound Foil Winding Design / M. Nigam, C.R. Sullivan // IEEE Industry Applications Society Annual Meeting, Oct. 2008. 140 113. Reddy, A. Hottest spot and life evaluation of power transformer design using finite element method / A. Reddy, M. Vijaykumar // Journal of Theoretical and Applied Information Technology, 2008. 114. Reducing losses in distribution transformers / Juan C. Olivares [and others] // IEEE Transactions on Power Delivery, Vol.18, No.7, 2003. 115. Susa, D. Dynamic thermal modelling of power transformers / Doctoral Dis- sertation. – Helsinki University of Technology, Espoo. – 2005. 116. Swift, G. A Fundamental Approach to Transformer Thermal Modeling – Part I: Theory and Equivalent Circuit / G. Swift, T. Molinski, W. Lehn // IEEE Transactions on Power Delivery, Vol. 16, No.2, April 2001. 141 ПРИЛОЖЕНИЕ 1 ПРОГРАММНЫЙ КОД ПОДСИСТЕМЫ ТЕПЛОВОГО АНАЛИЗА SHAPT В СРЕДЕ MATLAB %************************************************************* % Подсистема теплового расчета распределительного трансформатора на основе электрических схем замещения % %************************************************************* % Этап 1. Подготовка схемы (расчет числа ветвей и узлов соответсвтующих типов = evalin ( 'base' , 'NconL' ); NsecL = evalin ( 'base' , 'NsecL' ); [p,q,p_lmb,p_alfa,p_air,p_is,p_power,q_0,q_lmb,q_alfa,q_can, q_is,p_gor,p_lmb_down,p_z,p_sector] = PrepareNet (NconL,NsecL); % p - число ветвей схемы замещения q - число узлов схемы замещения p_lmb - число ветвей схемы замещения, характеризующих теплопередачу путем теплопроводности в обмотке p_alfa - число ветвей схемы замещения, характеризующих теплопередачу путем теплоотдачи с поверхности p_air - число ветвей схемы замещения, характеризующих теплопередачу в канале охлаждения p_is - число ветвей схемы замещения, характеризующих теплопередачу в бортовой изоляции путем совместного действия теплопроводности и теплоотдачи p_power - число ветвей схемы замещения, характеризующих источники тепловыделения q_0 - нулевой узел q_lmb - число узлов, характеризующих начало ветвей теплопроводности q_alfa - число узлов, характеризующих начало ветвей теплоотдачи q_can - число узлов, характеризующих начало ветвей теплопередачи в канале охлаждения q_is - число узлов, характеризующих начало ветвей бортовой изоляции p_gor - число горизонтальных ветвей водной секции p_lmb_down - число вертикальных ветвей теплопроводности в обмотке снизу первой секции p_z - число вертикальных ветвей теплопроводности в обмотке p_sector - число ветвей в секторе (сектор - совокупность ветвей, охватывающих все концентры в пределах одной секции NconL - число концентров обмотки НН % NsecL - число секций обмотки НН function [p,q,p_lmb,p_alfa,p_air,p_is,p_power,q_0,q_lmb,q_alfa,q_can, q_is,p_gor,p_lmb_down,p_z,p_sector] = PrepareNet (NconL,NsecL) p_lmb = NconL*(3*NsecL+1); p_alfa = 2*NconL*NsecL; p_air = NsecL*(NconL+3); p_is = 2*(2*NconL+1); p = p_lmb+p_alfa+p_air+p_is; p_power = NconL*NsecL; q_0 = 1; q_lmb = NconL*NsecL; q_alfa = 2*(p_alfa-NsecL*(NconL-1)); q_can = NsecL*(NconL-1); q_is = 3*NconL+1; q =q_0+ q_lmb+ q_alfa+ q_can+q_is; p_gor = 4*NconL+2; p_lmb_down = NconL; 142 p_z = 2*NconL + 1; p_sector = p_gor + p_z; %------------------------------------------------------------------------------------------------------------------ - assignin ( 'base' , 'q' ,q); assignin ( 'base' , 'p' ,p); assignin ( 'base' , 'p_lmb' ,p_lmb); assignin ( 'base' , 'p_alfa' ,p_alfa); assignin ( 'base' , 'p_air' ,p_air); assignin ( 'base' , 'p_is' ,p_is); assignin ( 'base' , 'p_power' ,p_power); assignin ( 'base' , 'q_0' ,q_0); assignin ( 'base' , 'q_lmb' ,q_lmb); assignin ( 'base' , 'q_alfa' ,q_alfa); assignin ( 'base' , 'q_can' ,q_can); assignin ( 'base' , 'q_is' ,q_is); assignin ( 'base' , 'p_gor' ,p_gor); assignin ( 'base' , 'p_lmb_down' ,p_lmb_down); assignin ( 'base' , 'p_z' ,p_z); assignin ( 'base' , 'p_sector' ,p_sector); % Этап 2. Расчет тепловых сопротивлений схемы замещения = evalin ( 'base' , 'NsLt' ); bpL = evalin ( 'base' , 'bpL' ); apL = evalin ( 'base' , 'apL' ); DvLt = evalin ( 'base' , 'DvLt' ); DnLt = evalin ( 'base' , 'DnLt' ); BcLt = evalin ( 'base' , 'BcLt' ); alfa_rad = 0.00000001; %ek_oil = 19; %19; [R_lmb_r_in, R_lmb_r_out, R_lmb_z,R_alfa_conv_in,R_alfa_conv_out, R_alfa_rad_in, R_alfa_rad_out, R_oil_can, Rbort] = ResistanceCalcTMG (NconL,NsLt,bpL,apL,DvLt,DnLt,alfa_rad, NsecL,dTms,BcLt,Ts,P1,P2,P1dob,P2dob); % Входные параметры - число концентров обмотки НН %NsLt - число слоев концентров обмотки НН %bpL - осевой размер неизолированного проводника НН %apL - радиальный размер неизолированного проводника НН %DvLt - внутренний диаметр концентра обмотки НН %DnLt - наружный диаметр концентра обмотки НН %alfa_rad - коэффициент радиационной теплоотдачи - число секций обмотки НН %dTms - превышение температуры средних слоев масла над температурой окружающей среды - ширина канала охлаждения - температура поверхности концентра %P1 - потери в обмотке ВН %P2 - потери в обмотке НН %P1dob - добавочные потери ВН %P2dob - добавочные потери НН % Выходные параметры % R_lmb_r_in - тепловое сопротивление теплопроводности обмотки в радиальном направлении в сторону внутренней поверхности концентра % R_lmb_r_out - тепловое сопротивление теплопроводности обмотки в радиальном направлении в сторону наружной поверхности концентра % R_lmb_z - тепловое сопротивление теплопроводности обмотки в осевом направлении 143 % R_alfa_conv_in - тепловое сопротивление конвективной теплоотдачи с поверхности обмотки в сторону внутренней поверхности концентра % R_alfa_conv_out - тепловое сопротивление конвективной теплоотдачи с поверхности обмотки в сторону наружной поверхности концентра % R_alfa_rad_in - тепловое сопротивление радиационной теплоотдачи с поверхности обмотки в сторону внутренней поверхности концентра % R_alfa_rad_out - тепловое сопротивление радиационной теплоотдачи с поверхности обмотки в сторону наружной поверхности концентра % R_oil_can - тепловое сопротивление теплопередачи в канале охлаждения Rbort - тепловое сопротивление теплопередачи в бортовой изоляции function [R_lmb_r_in, R_lmb_r_out, R_lmb_z,R_alfa_conv_in,R_alfa_conv_out, R_alfa_rad_in, R_alfa_rad_out, R_oil_can, Rbort] = ResistanceCalcTMG (NconL,NsLt,bpL,apL,DvLt,DnLt,alfa_rad, NsecL,dTms,BcLt,Ts,P1,P2,P1dob,P2dob) R_lmb_r_in = zeros (1,NconL); % сопротивления теплопроводности в радиальном направлении внутрь обмотки = zeros (1,NconL); % сопротивления теплопроводности в радиальном направлении наружу из обмотки = zeros (1,NconL); % сопротивления теплопроводности в осевом направлении = zeros (1,NconL); % сопротивления теплоотдачи конвекцией с поверхности в канал внутри обмотки = zeros (1,NconL); % сопротивления теплоотдачи конвекцией с поверхности в канал снаружи обмотки = zeros (1,NconL); % сопротивления теплоотдачи излучением с поверхности в канал внутри обмотки = zeros (1,NconL); % сопротивления теплоотдачи излучением с поверхности в канал снаружи обмотки = zeros (1,NconL-1); % тепловое сопротивление масляного канала с учетом конвекции масла = zeros(1,NconL); S_in = zeros (1,NconL); S_out = zeros (1,NconL); Scan = zeros (1,NconL-1); lmb_al = 230; lmb_oil = 0.15; lmb_eq = lmb_al*lmb_oil/(lmb_al+lmb_oil); %Pp=P1+P2+P1dob+P2dob; % Rbort = 10^-9*Pp^2-3*10^-5*Pp+0.2595; % Rkretnn = 2*10^-7*Pp^2+0.0036*Pp-20.71; Rbort = 0.085; Rkretnn=10; g=9.81; %T = 65; %Ts = 60; ro_20 = 860; alfa1 = 0.000686; lmb_0 = 0.15; h_con = 0.47; beta=0.0006; Co=0.75; %Tcoil=473; nu = 41.456*exp(-0.029*dTms)*10^-6; % кинематическая вязкость ro = ro_20*(1-alfa1*(dTms-20))/1000; % плотность cp = 4.1868*(1/sqrt(ro))*(0.403+0.00081*dTms)*1000; % удельная теплоемкость lmb = lmb_0*(1-(0.018/200)*dTms); % теплопроводность 144 Pr = nu*ro*cp*1000/lmb; Gr =((g*h_con^3)/nu^2)*beta*(Ts-dTms); Ra = Gr*Pr; Nu = Co*Ra^0.25; alfa_TMG = Nu*lmb/h_con; if Ra > 10^3 if Ra < 10^6 ek_oil = 0.105*Ra^0.3; end end if Ra > 10^6 if Ra < 10^12 ek_oil = 0.4*Ra^0.2; end end for i =1:1:NconL DavLt(i) = 0.5*(DvLt(i)+DnLt(i)); S_in(i) = pi*DvLt(i)*bpL*10^-3; S_out(i) = pi*DnLt(i)*bpL*10^-3; % 1. В радиальном направлении внутрь теплопроводность R_lmb_r_in (i) =NsLt(i)*(1/(2*pi*lmb_eq))*log(DavLt(i)/DvLt(i)); % 2. В радиальном направлении наружу теплопроводность R_lmb_r_out (i) =NsLt(i)*(1/(2*pi*lmb_eq))*log(DnLt(i)/DavLt(i)); % 3. В осевом направлении R_lmb_z (i) = (bpL*10^-3/NsecL)/(lmb_al*(pi/4)*((DvLt(i)+apL*10^-3*NsLt(i))^2-DvLt(i)^2)); % 4. Конвекция R_alfa_conv_in (i)= 1/(2*alfa_TMG*S_in(i)); R_alfa_conv_out (i) = 1/(alfa_TMG*S_out(i)); % 5. Излучение R_alfa_rad_in = 1/(alfa_rad*S_in(i)); R_alfa_rad_out = 1/(alfa_rad*S_out(i)); end % 6. Теплопроводность с конвекцией в масле for i =1:1:NconL - 1 Scan(i) = (pi/4)*(DvLt(i+1)^2-DnLt(i)^2); R_oil_can(i) = ((bpL*10^-3)/NsecL)/(lmb_oil*ek_oil*Scan(i))/(BcLt(i)*1000*Rkretnn); end assignin ( 'base' , 'R_lmb_r_in' ,R_lmb_r_in); assignin ( 'base' , 'R_lmb_r_out' ,R_lmb_r_out); assignin ( 'base' , 'R_lmb_z' ,R_lmb_z); assignin ( 'base' , 'R_alfa_conv_in' ,R_alfa_conv_in); assignin ( 'base' , 'R_alfa_conv_out' ,R_alfa_conv_out); assignin ( 'base' , 'R_alfa_rad_in' ,R_alfa_rad_in); assignin ( 'base' , 'R_alfa_rad_out' ,R_alfa_rad_out); assignin ( 'base' , 'R_oil_can' ,R_oil_can); % Этап 3. Построение матрицы соединений А, заполнение векторов Т и N00 [A00,N00,T00] = Matrix- ANet(p,q,p_is,p_gor,p_sector,NconL,NsecL,R_lmb_r_in,R_lmb_r_out,R_lmb_z,R_alfa_conv_in, 145 R_alfa_conv_out,R_alfa_rad_in,R_alfa_rad_out,R_oil_can,Rbort, P2, P2dob); assignin ( 'base' , 'A00' ,A00); assignin ( 'base' , 'N00' ,N00); assignin ( 'base' , 'T00' ,T00); % А - матрица соединений тепловой цепи электрической схемы замещения N00 - матрица номиналов элементов тепловой цепи электрической схемы замещения Т - матрица типов элементов тепловой цепи электрической схемы замещения function [A00M,N00,T00] = Matrix- ANet(p,q,p_is,p_gor,p_sector,NconL,NsecL,R_lmb_r_in,R_lmb_r_out,R_lmb_z,R_alfa_conv_in, R_alfa_conv_out,R_alfa_rad_in,R_alfa_rad_out,R_oil_can,Rbort, P2, P2dob) A00M = zeros(q,p); T00 = zeros(1,p+NsecL*NconL); N00 = zeros(1,p+NsecL*NconL); %A00 = zeros(q,p); % заполняем матрицу типов for i=1:1:p T00(i) = 3; end for i = p + 1:1: p + NconL * NsecL T00(i) = 5; end P_vet = zeros (1,p+NconL+NsecL*NconL); Q_uzl = zeros(1,q+1); for i =1:1:p+NsecL*NconL P_vet(i) = i; end for j=1:1:q+1 Q_uzl(j) = начала источников тока for i = p + 1:1:p + NconL* NsecL A00M(Q_uzl(1), P_vet(i) ) = 1; end % определяем начала и концы ветвей снизу - готово j =1; for i = 1 :1: 2 * NconL + 1 if mod (i,2) == 0 A00M(Q_uzl(0.5 * i + 1), P_vet(i) ) = 1; N00(P_vet(i))=Rbort; else A00M( Q_uzl(NconL + 2 * i), P_vet(i) ) = 1; N00(P_vet(i))=R_oil_can(j); % проверить индексы (учесть канал между сердечником, канал ОНН и главный канал A00M(Q_uzl(1), P_vet(i) ) = -1; A00M(Q_uzl(1), P_vet(i) ) = -определяем начала и концы ветвей сверху for i = 1:1: p_is / 2 A00M(Q_uzl(q - p_is / 2 + i), P_vet(p - p_is / 2 + i)) = 1; A00M(Q_uzl(1), P_vet(p - p_is / 2 + i)) = -1; if mod(NconL,2) == 0 if mod(i,2) == 0 N00(P_vet(p - p_is / 2 + i)) = Rbort; 146 else N00(P_vet(p - p_is / 2 + i)) = R_oil_can(j); % проверить индексы end else if mod(i,2) == 0 N00(P_vet(p - p_is / 2 + i)) = Rbort; else N00(P_vet(p - p_is / 2 + i)) = R_oil_can(j); % проверить индексы end end определяем начала и концы нижних ветвей (их число равно числу концентров) j = 1; for i = 2 * NconL + 2:1:3 * NconL + 1 A00M(Q_uzl(NconL + 4 * j), P_vet(i) ) = 1; A00M(Q_uzl(j + 1), P_vet(i) ) = -1; N00( P_vet(i)) = R_lmb_z(j) / 2; j = j + 1; end % определяем начала и концы ветвей слева и справа v =1; % временно для запуска программы j = 3 * NconL + 2; k = NconL + 2; k1 = 5 * NconL + 2; for i = 1:1: NsecL A00M( Q_uzl(1), P_vet(j) ) = -1; конец air ветви слева A00M(Q_uzl(k), P_vet(j) ) = 1; начало горизонтальной air ветви слева A00M(Q_uzl(k), P_vet(j + 4 * NconL + 2)) = 1; начало вертикальной air ветви слева A00M(Q_uzl(k + 4 * NconL + 1), P_vet(j + 4 * NconL + 2)) = -1; конец вертикальной air ветви слева A00M(Q_uzl(1), P_vet(j + p_gor - 1)) = -1; %' конец air ветви справа A00M(Q_uzl(k1), P_vet(j + p_gor - 1)) = 1; %' начало air ветви справа A00M(Q_uzl(k1), P_vet(j + p_gor - 1 + 2 * NconL + 1)) = 1; %' начало вертикальной air ветви справа N00(P_vet(j)) = R_oil_can(v); N00(P_vet(j + 4 * NconL + 2)) = R_oil_can(v); N00(P_vet(j + p_gor - 1)) = R_oil_can(v); N00(P_vet(j + p_gor - 1 + 2 * NconL + 1)) = R_oil_can(v); if k1 + 4 * NconL + 1 < q A00M(Q_uzl(k1 + 4 * NconL + 1), P_vet(j + p_gor - 1 + 2 * NconL + 1)) = -1; % ' конец вертикальной ветви справа else A00M(Q_uzl(k1 + 2 * NconL + 1), P_vet(j + p_gor - 1 + 2 * NconL + 1)) = -1; %' конец вертикальной ветви справа end j = j + p_sector; k = k + 4 * NconL + 1; k1 = k1 + 4 * NconL + Определяем характеристические узлы k = 0; k1 = 0; num = 1; for j = 1:1:NsecL for i = 1:1:NconL h = 5 * NconL + 4 * i + 1 + k; A00M(Q_uzl(NconL + 4 * i + k - 1), P_vet(3 * NconL + 4 * i + k1 - 1)) = 1; %alfa ветвь слева начало A00M(Q_uzl(NconL + 4 * i + k - 1 - 1), P_vet(3 * NconL + 4 * i + k1 - 1)) = -1; %alfa ветвь слева конец A00M(Q_uzl(NconL + 4 * i + k), P_vet(3 * NconL + 4 * i + k1)) = 1; %lambda ветвь слева начало 147 A00M(Q_uzl(NconL + 4 * i + k), P_vet(p + num)) = -1; источника тока ветвь конец' lambda ветвь слева конец A00M(Q_uzl(NconL + 4 * i + k), P_vet(7 * NconL + 3 + 2 * i + k1)) = 1; % ' lambda ветвь вверх начало N00(P_vet(3 * NconL + 4 * i + k1 - 1) ) = R_alfa_conv_in(i)*R_alfa_rad_in(v)/(R_alfa_rad_in(v)+R_alfa_conv_in(i)); % !!!!!!!!!!!! N00(P_vet(3 * NconL + 4 * i + k1)) = R_lmb_r_in(i) / 2 ; N00(P_vet(7 * NconL + 3 + 2 * i + k1)) = R_lmb_z(i) / 2; N00(P_vet(3 * NconL + 4 * i + 1 + k1)) = R_lmb_r_out(i) / 2; N00(P_vet(3 * NconL + 4 * i + 1 + k1 + 1)) = R_alfa_conv_out(i)*R_alfa_rad_out(v)/(R_alfa_rad_out(v)+R_alfa_conv_out(i)); % !!!!!!!!!!!! if h < q - (p_is / 2 - 1) A00M(Q_uzl(NconL + 4 * i + 4 * NconL + 1 + k), P_vet(7 * NconL + 3 + 2 * i + k1)) = -1; %' lamb- da ветвь вверх конец else A00M(Q_uzl(NconL + 4 * NconL + 2 * (i - 1) + 4 + k), P_vet(7 * NconL + 3 + 2 * i + k1)) = -1; % ' lambda ветвь вверх конец end A00M(Q_uzl(NconL + 4 * i + k), P_vet(3 * NconL + 4 * i + 1 + k1)) = 1; %' lambda ветвь справа начало %' lambda ветвь справа конец A00M(Q_uzl(NconL + 4 * i + k + 1), P_vet(3 * NconL + 4 * i + 1 + k1 + 1)) = 1; %' alfa ветвь справа начало A00M(Q_uzl(NconL + 4 * i + k + 1 + 1), P_vet(3 * NconL + 4 * i + 1 + k1 + 1)) = -1; %' alfa ветвь справа конец num = num + 1; end k = k + 4 * NconL + 1; k1 = k1 + 6 * NconL + 3; end % 'Строим ветви в канале nk = NconL - 1; if nk > 0 k = 0; k1 = 0; for j = 1:1:NsecL for i = 1:1:nk h = 5 * nk + 4 * i + k + 8; A00M(Q_uzl(nk + 4 * i + 3 + k), P_vet(7 * nk + 13 + 2 * (i - 1) + k1)) = 1; T00(P_vet(7 * nk + 13 + 2 * (i - 1) + k1)) = 6; N00(P_vet(7 * nk + 13 + 2 * (i - 1) + k1) ) = R_oil_can(v); if h < q - (p_is / 2 - 2) A00M(Q_uzl(nk + 4 * i + 3 + k + 4 * nk + 5), P_vet(7 * nk + 13 + 2 * (i - 1) + k1)) = -1; else A00M(Q_uzl(nk + 4 * i + 3 + k + 4 * nk + 3 - 2 * (i - 1)), P_vet(7 * nk + 13 + 2 * (i - 1) + k1)) = -1; end end k = k + 4 * nk + 5; k1 = k1 + 6 * nk + 9; end end for i = p + 1:1: p + NconL*NsecL N00(i) = (P2+P2dob)/NsecL; end [m,n] = size(N00); T0 = T00'; TN0 = zeros(n,1); for i=1:1:n 148 switch T0(i) case 6 T0(i) = 3; TN0(i) = 10000000; end end N0 = zeros(n,n); for i=1:1:n N0(i,i) = N00(i); end A0 = A00; [A,T,N,Tur,p,q,nx,y1,y2,z1,z2,NCepi,NPrBase,NPrOld] = AnalisNet(A0,T0,N0); assignin ( 'base' , 'A' ,A); assignin ( 'base' , 'T' ,T); assignin ( 'base' , 'N' , N); assignin ( 'base' , 'Tur' , Tur); assignin ( 'base' , 'p' , p); assignin ( 'base' , 'q' , q); assignin ( 'base' , 'nx' , nx); assignin ( 'base' , 'y1' ,y1); assignin ( 'base' , 'y2' ,y2); assignin ( 'base' , 'z1' ,z1); assignin ( 'base' , 'z2' ,z2); assignin ( 'base' , 'NCepi' , NCepi); assignin ( 'base' , 'NPrBase' , NPrBase); assignin ( 'base' , 'NPrOld' , NPrOld); [G12,G21,Dy,Dyt,Cz,Czt] = GetPrewMatrix(A,p,q,y1,y2,z1,z2,NCepi); [ni,nr,np,Gr,Grp,Gp,Xp,Xr] = InitCalcNet(T,Tur,p,nx); I = zeros(p,1); [G,Y,Z] = GetPostMatrix(Dy,Dyt,Cz,Czt,G12,G21,NPrOld,T0,TN0,N0,N,I,p,y1,y2); tn = 0; tk = 0.03; dt = 0.01; nt = (tk - tn) / dt + 1; X = zeros(nx,1); %Xt = zeros(nx,nt); содержит результаты расчета вовсе моменты времени = zeros(p,nt); содержит результаты расчета вовсе моменты времени = zeros(p,nt); содержит результаты расчета вовсе моменты времени = zeros(p,nt); содержит результаты расчета вовсе моменты времени t = tn; iter = 0; while t < tk X0 = X; [Fi,U,I] = GetUI(A0,X,p,q,nx,y1,NPrBase,N0,T0); X = BuildX(Tur,Y,Z,Xp,Xr,nx,ni); iter = iter + 1; %Xt(:,iter) = X(:); Fit(:,iter) = Fi(:); Ut(:,iter) = U(:); It(:,iter) = I(:); %for i=2:1:nx %line([tn,t],[X0(i),X(i)],'EraseMode','none', 'color','blue','LineWidth',1); %end [G,Y,Z] = GetPostMatrix(Dy,Dyt,Cz,Czt,G12,G21,NPrOld,T0,TN0,N0,N,I,p,y1,y2); [Xp,Xr] = StepCalcNet(G,Gr,Grp,Gp,Xp,X,Tur,nx,nr,np,dt); tn = t; t = t + dt; end X = BuildX(Tur,Y,Z,Xp,Xr,nx,ni); [Fi,U,I] = GetUI(A0,X,p,q,nx,y1,NPrBase,N0,T0); iter = iter + 1; Fit(:,iter) = Fi(:); 149 Ut(:,iter) = U(:); It(:,iter) = I(:); T_uzl_c = zeros(NsecL,NconL); T_uzl_int = zeros(NsecL,NconL); T_uzl_ext = zeros(NsecL,NconL); T_uzl_air = zeros(NsecL,NconL); for j = 1:1:NconL for i = 1:1:NsecL T_uzl_c(i, j) = Fi(NconL + 4 * j + (4 * NconL + 1) * (i - 1)); T_uzl_int(i, j) = Fi(NconL + 4 * j + (4 * NconL + 1) * (i - 1) - 1); T_uzl_ext(i, j) = Fi(NconL + 4 * j + (4 * NconL + 1) * (i - 1) + 1); T_uzl_air(i, j) = Fi(NconL + 4 * j + (4 *NconL + 1) * (i - 1) + 2); end end subplot(2,3,1); plot(T_uzl_c, 'Color' , 'r' ); hold Центр тепловыделения' ); xlabel( ' Номер секции' ); ylabel( 'Превышение температуры on ; subplot(2,3,2); plot(T_uzl_int, 'Color' , 'g' ); hold Внутренние поверхности' ); xlabel( 'Номер секции' ); ylabel( 'Превышение температуры on ; subplot(2,3,3); plot(T_uzl_ext, 'Color' , 'b' ); hold Наружные поверхности' ); xlabel( 'Номер секции' ); ylabel( 'Превышение температуры on ; subplot(2,3,4); plot(T_uzl_air, 'Color' , 'k' ); title( 'Канал' ); xlabel( 'Номер секции' ); ylabel( 'Превышение температуры on ; tc1 = mean (T_uzl_c); tc = mean(tc1); ti1 = mean (T_uzl_int); ti = mean (ti1); te1= mean (T_uzl_ext); te= mean(te1); tair1= mean (T_uzl_air); tair= mean(tair1); err = ti - Ts; if err < 0 Ts = Ts - dterr; else Ts = Ts + dterr; end %end assignin ( 'base' , 'tc' ,tc); assignin ( 'base' , 'ti' ,ti); assignin ( 'base' , 'te' ,te); assignin ( 'base' , 'tair' ,tair); 150 ПРИЛОЖЕНИЕ 2 151 ПРИЛОЖЕНИЕ 3 152 ПРИЛОЖЕНИЕ 4 |