r * 0.02;
lambtp2=NumVar * 5;=0;=3;=zeros (Nx, Ny); _ g=zeros (Nx, Ny);=zeros (Nx, Ny);=zeros (Nx, Ny, Nz); _ etta_gamma_ksi=zeros (Nx, Ny, Nz);
U=zeros (Nx, Ny);
. Розрахунок даних:
Розрахунок управління:=1;=1509.977;=1.39;
E2=934.1678;=0.00099;=80.765158;=57.231289;=34.8664826;=9.278370;=0.05944;=zeros (5,5); _ prev=zeros (5,5); lmb=1:1 * 10 ^ 5=t + delt; etta=2: Nx - 1gamma=2: Ny - 1_g (etta, gamma)=10; etta=1: Nxgamma=1: Ny (etta, gamma)=T_g (etta, gamma) - T (etta, gamma, 4) etta=2: Nx - 1gamma=2: Ny - 1 (etta, gamma)=( F (etta - 1, gamma) - 2 * F (etta, gamma) + F (etta + 1, gamma))/delX ^ 2 + ...
(F (etta, gamma - 1) - 2 * F (etta, gamma) + F (etta, gamma + 1))/delY ^ 2; _prev (etta, gamma)=E4 * (((N4-1)/N4) * F (etta, gamma) - (1/N4) * lapl (etta, gamma)) * delt; (etta, gamma)=SUMM (etta, gamma) + SUMM_prev (etta , gamma);
% U (etta, gamma)=1000 * cos (pi * (i - 0.5) * etta/(Nx - 1)) * sin (pi * gamma * (i - 0.5)/( Ny - 1)); (etta, gamma)=E1 * (((N1-1)/N1) * F (etta, gamma) - (1/N1) * lapl (etta, gamma)) + SUMM (etta, gamma) + E2 * ((F (etta, gamma) -Fold (etta, gamma))/delt);
end
Теплові процеси, що протікають в печі
for etta=2: Nx - 1gamma=2: Ny - 1ksi=2: 6_etta_gamma_ksi (etta, gamma, ksi)=a1 * delt * ((T (etta - 1, gamma, ksi)- 2 * T (etta, gamma, ksi) + T (etta + 1, gamma, ksi))/delX ^ 2 + ...
(T (etta, gamma - 1, ksi) - 2 * T (etta, gamma, ksi) + T (etta, gamma + 1, ksi))/delY ^ 2 +.
(T (etta, gamma, ksi - 1) - 2 * T (etta, gamma, ksi) + T (etta, gamma, ksi + 1))/delZ1 ^ 2); etta= 2: Nx - 1gamma=2: Ny - 1ksi=8_etta_gamma_ksi (etta, gamma, ksi)=a2 * delt * ((T (etta - 1, gamma, ksi) - 2 * T (etta, gamma, ksi) + T (etta + 1, gamma, ksi))/delX ^ 2 + ...
(T (etta, gamma - 1, ksi) - 2 * T (etta, gamma, ksi) + T (etta, gamma + 1, ksi))/delY ^ 2 +
(T (etta, gamma, ksi - 1) - 2 * T (etta, gamma, ksi) + T (etta, gamma, ksi + 1))/delZ2 ^ 2); etta= 2: Nx - 1gamma=2: Ny - 1ksi=10: 14_etta_gamma_ksi (etta, gamma, ksi)=a3 * delt * ((T (etta - 1, gamma, ksi) - 2 * T (etta, gamma, ksi) + T (etta + 1, gamma, ksi))/delX ^ 2 + ...
(T (etta, gamma - 1, ksi) - 2 * T (etta, gamma, ksi) + T (etta, gamma + 1, ksi))/delY ^ 2 +
(T (etta, gamma, ksi - 1) - 2 * T (etta, gamma, ksi) + T (etta, gamma, ksi + 1))/delZ3 ^ 2);
end
Граничні умови при переході з середовища в середу
for etta=2: Nx - 1gamma=2: Ny - 1 (etta, gamma, 7)=(lambtp1 * T (etta, gamma, 6) + lambtp2 * T (etta, gamma, 8))/(lambtp1 + lambtp2); (etta, gamma, 9)=(lambtp1 * T (etta, gamma, 10) + lambtp2 * T (etta, gamma, 8))/(lambtp1 + lambtp2); ksi= 2: Nz - 1gamma=2: Ny - 1etta=2: Nx - 1 (etta, gamma, ksi)=T (etta, gamma, ksi) + delT_etta_gamma_ksi (etta, gamma, ksi);
Вхідні воздействіяetta=2: Nx - 1gamma=2: Ny - 1
% T (etta, gamma, 1)=delZ * utp/lambtp1 + T (etta, gamma, 2); (etta, gamma, 1)=delZ * U (etta, gamma)/ lambtp1 + T (etta, gamma, 2);
Граничні умови: etta=1: Nxgamma=1: Nyksi=1: Nz (Nx, gamma, ksi)=T (Nx - 1, gamma, ksi); (etta, Ny, ksi) =T (etta, Ny - 1, ksi); (etta, gamma, Nz)=T (etta, gamma, Nz - 1); (Nx, gamma)=T (Nx - 1, gamma); (etta, Ny )=T (etta, Ny - 1); etta=1: Nxgamma=1: Ny (etta, gamma)=F (etta, gamma); (lmb, :)=T (4,2,6); (1 ), plot (lmb, T (4,4,12)), hold on (2), plot (lmb, U), hold on
Розрахунок характеристик:
n=length (ForPlotT);
M=ForPlotT (n); _ vv=U (4,4);=M/M_vv_x=pi * (i - 0.5)/Lx; _y=pi * (i - 0.5)/Ly ;=Psi_x ^ 2 + Psi_y ^ 2
