De volgcode was gemaakt om de tijd te berekenen om de TE te vergroten ... (2024)

functietime_to_increase_temperature = calculate_time_to_increase_temperature (...

rho_concrete, rho_air, k_concrete, k_air, ~, cp_concrete,...

cp_air, t_initial, t_surface, l1, l2, l3, volume_p1_p3, ~,...

NX, NT, Delta_T, Heat_power)

% Bereken de totale dikte

Total_thickness = l1 + l2 + l3;

% Thermische diffusiviteiten

alpha_concrete = k_concrete / (rho_concrete * cp_concrete);% m²/s

alpha_air = k_air / (rho_air * cp_air);% m²/s

% Discretisatieparameters

dx = total_thickness / nx;% Ruimtelijke stapgrootte

dt = 6 / nt;% Tijdstapgrootte

% Preallocate temperatuurarrays

T = t_initial * ones (nx, 1);

T_new = t;

% Definieer materiaaleigenschappen op elk punt

len_p1 = round (nx * l1 / total_thickness);% Lengte van betonlaag P1

len_p2 = round (nx * l2 / total_thickness);% Lengte van luchtlaag P2

alleen_p3 = nx - alleen_p1 - alleen_p2;% Lengte van betonlaag P3

% Zorgt ervoor dat LEN_P3 correct is aangepast als afronding problemen heeft veroorzaakt

alsalleen_p1 + alleen_p2 + alleen_p3 ~ = nx

alleen_p3 = nx - alleen_p1 - alleen_p2;

einde

% Preallocate materiële array

materiaal = nullen (nx, 1);

materiaal (1: Len_p1) = alpha_concrete;

materiaal (Len_P1+1: Len_P1+Len_P2) = alpha_air;

materiaal (LEN_P1+LEN_P2+1: nx) = alpha_concrete;

tijd om de temperatuur te verhogen = in f;% Initialiseer om de tijd op te slaan om de temperatuur te verhogen

% Bereken warmte -invoer per ruimtelijk punt

hitte_input_per_point = heat_power * dt / (rho_concrete * cp_concrete * volume_p1_p3 / nx);

% Tijdintegratie met behulp van de eindige verschilmethode (FDM)

tijd = 0;% Initialiseer de tijd

voorn = 1: nt

voorI = 2: NX-1

T_new (i) = t (i) + materiaal (i) * dt / dx^2 * (t (i + 1) - 2 * t (i) + t (i -1));

einde

% Grensvoorwaarden

T_new (1) = t (1) + alpha_concrete * dt / dx^2 * (t (2) - t (1)) + heat_input_per_point * (t_surface - t (1));

T_new (nx) = t (nx-1);

% Updatetemperatuur

T = t_new;

tijd = tijd + dt;

% Controleer of de temperatuur aan de rechterkant van P3 T_INITIAL + DELTA_T bereikt

alsT (nx)> = t_initial + delta_t

time_to_increase_temperature = time;

pauze;

einde

einde

% Retourneer de tijd om de temperatuur met 1 K aan de rechterkant van P3 te verhogen

einde

functiemain_process ()

% Parameters met vaste waarden

rho_concrete = 2500;% kg/m³

rho_air = 1293;% kg/m³

k_concrete = 1,35;% W/(m · k)

k_air = 0.024;% W/(m · k) uitgaande van typische waarde voor lucht

cp_concrete = 0,92 * 1000;% J/(kg · k)

Cp_air = 1.005 * 1000;% J/(kg · k)

T_initial = 300;% Initiële temperatuur in K

T_surface = 411;% Oppervlaktetemperatuur in k

volume_p1_p3 = 0.00488;% Volume van elke betonnen plaat in m³

Nx = 100;% Aantal ruimtelijke punten

NT = 1000;% Aantal tijdstappen

delta_t = 1;% Temperatuurstijging in K

% Vraag de gebruiker om specifieke invoerparameters

% L1 = input ('Voer de dikte van de eerste betonlaag in (m):');

% L2 = input ('Voer de dikte van de luchtlaag in (m):');

% L3 = input ('Voer de dikte van de tweede betonlaag in (m):');

% heat_power = input ('Voer het warmtevermogen in (w):');

L1 = 1;

Nee = 1;

L = 1;

heat_power = 1;

% Bel de functie

time_to_increase_temperature = calculate_time_to_increase_temperature (...

rho_concrete, rho_air, k_concrete, k_air, [], cp_concrete,...

cp_air, t_initial, t_surface, l1, l2, l3, volume_p1_p3, [],...

NX, NT, Delta_T, Heat_power);

% Geeft het resultaat weer

fprintf ('Tijd om de temperatuur met 1 K aan de rechterkant van P3 te verhogen: %.2f seconden \ n', time_to_increase_temperature);

% Plot het resultaat

Total_thickness = l1 + l2 + l3;

x = linspace (0, total_thickness, nx);% Ruimtelijke punten

y = linspace (0, time_to_increase_temperature, nx);% Tijdelijke punten voor het plotten

Plot (X, Y,'Weergavenaam',,'Tijd om de temperatuur met 1 K te verhogen');

xlabel ('Positie (M)');

YLabel ('Keer)');

legende('show');

titel('Tijd om de temperatuur met 1 K aan de rechterkant van P3 te verhogen');

einde

main_process ()

⋮