Termika i zarządzanie cieplne pocisku

Publikacja: 2026-09-04 · Aktualizacja: 2026-09-11

Stan cieplny pocisku jest pamięcią całego jego życia. Temperatura baterii przy uruchomieniu zależy nie tylko od pogody w chwili startu, lecz także od magazynowania, transportu, czasu na wyrzutni i wcześniejszego poboru mocy. Podobnie temperatura detektora nie jest temperaturą jego obudowy, a temperatura złącza półprzewodnika nie jest temperaturą powietrza w sekcji. Projektant nie rozwiązuje więc zadania „odporności na zakres temperatur”, ale odtwarza zmienne w czasie pole temperatur i sprawdza, czy każda potrzebna funkcja ma margines w odpowiedniej fazie misji.

Teza tej monografii brzmi: poprawność cieplną pocisku można wykazać tylko wtedy, gdy zachowuje się ciągły ślad od profilu życia i przepływu energii, przez temperatury konkretnych węzłów oraz ich wpływ na funkcje, po model, próbę i konfigurację seryjnego wyrobu. Ten ślad jest ważniejszy niż imponująca mapa barwna z symulacji. Mapa bez stanu początkowego, mocy, warunków brzegowych, niepewności i limitów funkcjonalnych nie odpowiada jeszcze na żadne pytanie eksploatacyjne.

Opis pozostaje na poziomie fizyki i inżynierii systemowej. Nie podaje progów bezpieczeństwa, parametrów wykonawczych ani procedur obsługi współczesnych środków bojowych. W przypadku przegrzanego, uszkodzonego lub nieznanego obiektu obowiązują wyłącznie autoryzowane instrukcje systemu, zabezpieczenie miejsca i wezwanie właściwych służb lub personelu usuwania materiałów wybuchowych (EOD).

Jeden łańcuch dowodowy zamiast tabeli temperatur

Praca termiczna zaczyna się od osi czasu. Dla każdej fazy określa się środowisko zewnętrzne, stan zasilania, źródła wewnętrzne i konfigurację mechaniczną. Wynik fazy staje się warunkiem początkowym następnej. Dopiero z tego opisu powstają temperatury węzłów, a z nich — odpowiedź baterii, czujników, elektroniki, urządzeń wykonawczych, materiałów i połączeń. Na końcu porównuje się przewidywaną odpowiedź, wraz z niepewnością, z wymaganiem na konkretną funkcję.

Pełny ślad ma postać:

$$\text{profil życia}\rightarrow\text{obciążenia cieplne}\rightarrow \text{drogi przepływu}\rightarrow T_i(t)\rightarrow \text{parametr funkcjonalny}\rightarrow\text{margines}\rightarrow\text{dowód}.$$

Brak któregokolwiek ogniwa prowadzi do typowego błędu. Znany zakres temperatury otoczenia bez dróg przepływu nie wyznacza temperatury podzespołu. Znana temperatura obudowy bez modelu złącza nie wyznacza temperatury krzemu. Udana próba jednego egzemplarza bez identyfikacji jego materiałów, oprogramowania i historii nie kwalifikuje całej populacji.

Standard ECSS-E-ST-31 traktuje definicję, analizę, projekt, wytwarzanie, weryfikację i eksploatację kontroli cieplnej jako jedną dyscyplinę.¹ Nie jest to norma uzbrojenia, ale logika ciągłości dowodu pozostaje użyteczna po dopasowaniu do krótkiej misji atmosferycznej, długiego magazynowania, powtarzanego przenoszenia oraz właściwych wymagań wojskowych.

Różne klasy pocisków, różne pamięci cieplne

Nie istnieje jeden reprezentatywny „profil rakiety”. O dominującym zjawisku decydują czas, sposób przechowywania, platforma, źródło napędu i stan wyposażenia przed startem.

Klasa i sytuacja	Dominujący stan początkowy	Typowe źródła podczas misji	Co najłatwiej przeoczyć
pocisk kontenerowy	wielomiesięczne składowanie i mikroklimat kontenera	krótki test, elektronika i krótkotrwały lot	różnicę między temperaturą kontenera a wnętrzem pocisku
pocisk przenoszony przez samolot	wielogodzinne przenoszenie, wysokość i sąsiedztwo płatowca	głowica naprowadzająca, tor radiowy, komputer, nagrzewanie aerodynamiczne	powrót bez odpalenia jako kolejny cykl resursowy
pocisk bardzo szybki	temperatura skóry przed startem	intensywny, krótki strumień zewnętrzny i napęd	opóźnione maksimum temperatury wnętrza po maksimum strumienia
pocisk manewrujący	stan wyrzutni i przygotowania	długotrwały napęd, elektronika, RF i otoczenie	narastanie temperatury do stanu bliskiego ustalonemu
efektor wysokościowy	stan naziemny	szybko zmieniająca się konwekcja, ciśnienie i nagrzewanie lotne	zmianę mechanizmu wymiany wraz z wysokością

MIL-STD-810H porządkuje badania przez wpływy występujące w całym cyklu życia, a nie przez bezmyślne zastosowanie jednego katalogu poziomów.⁶ Dla termiki oznacza to rozdzielenie gorącego magazynowania, zimnego uruchomienia, udaru cieplnego, promieniowania słonecznego, wilgoci i pracy. Próba obejmująca jedno z tych zjawisk nie staje się automatycznie dowodem dla pozostałych.

Cztery granice: przetrwanie, praca, jakość i trwałość

Jedna wartość graniczna nie opisuje zachowania podsystemu. Wymagania powinny rozróżniać:

Rodzaj granicy	Pytanie dowodowe	Przykład skutku
przetrwanie	czy nie powstanie nieodwracalne uszkodzenie?	pęknięcie, rozszczelnienie, trwała degradacja materiału
dopuszczenie pracy	czy podsystem można bezpiecznie uruchomić?	rozruch komputera, zdolność baterii, ruch mechanizmu
jakość funkcji	czy działanie ma wymaganą dokładność?	błąd inercyjnej jednostki pomiarowej (IMU), szum detektora, transmisja radomu, stabilność zegara
trwałość	ile czasu i cykli można zaakceptować?	zmęczenie lutów, starzenie kleju, utrata własności baterii

Każda granica wymaga miejsca pomiaru i czasu. Sformułowanie „elektronika pracuje do temperatury X” jest niepełne, dopóki nie wiadomo, czy chodzi o powietrze, obudowę, płytkę czy złącze. Zakres dopuszczenia pracy może być szerszy niż zakres pełnej dokładności, a zakres przetrwania szerszy od obu. Układ może się uruchomić, lecz potrzebować czasu na osiągnięcie dokładności; może też przetrwać gorące składowanie, ale zużyć przy tym część trwałości.

Wniosek o trwałości wymaga również przedziału, nie tylko jednej krzywej. Rozważmy czysto hipotetyczny mechanizm starzenia, dla którego badania materiałowe potwierdziły w danym zakresie model przyspieszenia

$$a(T;b)=\exp\left[b\left(\frac{1}{T_{ref}}-\frac{1}{T}\right)\right].$$

Nie jest to uniwersalny model baterii, kleju ani elektroniki: wolno go zastosować dopiero po wykazaniu, że w badanym przedziale nie zmienia się dominujący mechanizm. Niech $T_{ref}=313 \mathrm K$, czas ekspozycji wynosi 100 godzin, a budżet odniesienia 1000 godzin. Dla niepewności $T\in[331;335] \mathrm K$ i $b\in[3200;4200] \mathrm K$ współczynnik przyspieszenia mieści się w przedziale $a\in[1{,}744;2{,}414]$, więc przyrost znormalizowanego uszkodzenia wynosi $\Delta D^*\in[0{,}174;0{,}241]$.

Jeżeli wcześniejsza historia daje $D_0^*\in[0{,}60;0{,}68]$, łączny wynik to

$$D^*=D_0^*+\Delta D^*\in[0{,}774;0{,}921], \qquad R^*=1-D^*\in[0{,}079;0{,}226].$$

Wartość nominalna, obliczona dla $D_0^*=0{,}64$, $T=333 \mathrm K$ i $b=3700 \mathrm K$, daje $R^*=0{,}157$. Planowany kolejny cykl o zapotrzebowaniu $0{,}10$ mieści się w wartości nominalnej, lecz nie w dolnej granicy $0{,}079$. Poprawna decyzja nie polega więc na zaokrągleniu do wyniku pozytywnego: trzeba zawęzić niepewność temperatury, parametru materiałowego lub wcześniejszej historii, a jeśli nie jest to możliwe — zastosować konserwatywną granicę z dokumentacji właściwego wyrobu. Liczby nie opisują żadnego rzeczywistego resursu.

Bilans energii i skale czasu

Punktem wyjścia jest zachowanie energii:

$$\frac{dU}{dt}=\dot Q_{in}-\dot Q_{out}+\dot Q_{gen}.$$

Dla skupionego węzła o masie $m$, cieple właściwym $c_p$, temperaturze $T$ i mocy strat $P(t)$ można zapisać

$$mc_p\frac{dT}{dt}=\sum_j G_j(T_j-T)+P(t),$$

gdzie $G_j$ jest przewodnością do sąsiedniego węzła. Iloraz pojemności cieplnej i przewodności daje skalę czasu

$$\tau=\frac{mc_p}{G}.$$

Jeżeli faza trwa znacznie krócej niż $\tau$, węzeł zachowuje pamięć poprzedniego stanu. Jeżeli trwa wiele stałych czasowych, zbliża się do równowagi. To rozróżnienie wyjaśnia pozorny paradoks: krótki, wysoki strumień może najbardziej obciążyć cienką powłokę, a dłuższe oczekiwanie z włączoną elektroniką — głęboko położony procesor. Maksimum strumienia, całkowita dostarczona energia i maksimum temperatury nie muszą wystąpić jednocześnie.

NASA opisuje analogiczny bilans dla małych statków kosmicznych, gdzie energia pochłaniana, generowana, wypromieniowana i magazynowana musi zamknąć się przy ograniczonej masie, objętości i mocy.³ W atmosferycznym pocisku inne są warunki brzegowe, ale bilans i pamięć cieplna pozostają tą samą fizyką.

Drogi przepływu: materiał, styk, gaz i promieniowanie

Przewodzenie w jednorodnej warstwie opisuje prawo Fouriera:

$$\dot Q=-kA\frac{dT}{dx},\qquad R_{th}=\frac{L}{kA}.$$

W prawdziwej konstrukcji często nie gruba część, lecz cienki styk dominuje opór. Nacisk, płaskość, chropowatość, powłoka, klej, materiał interfejsowy i cykle zmieniają rzeczywistą powierzchnię kontaktu. Śruba może być jednocześnie drogą strukturalną i mostkiem cieplnym; wiązka może odprowadzać mało energii z całej sekcji, lecz wystarczająco dużo, by zaburzyć precyzyjny czujnik. Parametr styku skorelowany na prototypie nie reprezentuje serii, jeżeli proces montażu nie kontroluje docisku i materiałów.

Konwekcję przybliża zależność

$$\dot Q=hA(T_s-T_\infty),$$

lecz współczynnik $h$ zależy od geometrii, ruchu i właściwości płynu. Komora z wymuszonym obiegiem może ogrzewać szybciej niż magazyn, a próżnia usuwa konwekcję całkowicie. Żadne z tych środowisk nie jest zawsze „bardziej surowe”; trzeba porównać odpowiedź krytycznych węzłów.

Promieniowanie ma zależność nieliniową

$$\dot Q=\varepsilon\sigma A(T_s^4-T_{sur}^4).$$

Znaczenie mają współczynniki widoku, emisyjność, pochłanianie promieniowania słonecznego i stan powierzchni. Farba, zabrudzenie, wilgoć lub korozja mogą zmienić bilans. Izolacja nie usuwa energii: tylko spowalnia jej przepływ. Może chronić przed krótkim nagrzaniem zewnętrznym, a jednocześnie zatrzymywać moc wydzielaną w środku.

Od środowiska do sieci węzłów

Sieć cieplna zamienia geometrię w układ węzłów o pojemności i połączeń o przewodności. Rozdzielczość wynika z decyzji. Jeden węzeł może wystarczyć do oszacowania średniej temperatury masywnej baterii, ale nie do wyznaczenia temperatury złącza procesora, gradientu radomu albo lokalnej spoiny.

Abstrakcyjny przekrój pokazuje, dlaczego nazwa „temperatura elektroniki” ukrywa kilka stanów:

      promieniowanie                 moc układu Pj
            ↓                              ↓
otoczenie ⇄ [SKÓRA] ── Gm ── [OBUDOWA] ── Gi ── [PŁYTKA] ── Rj ── [ZŁĄCZE]
   konwekcja  Cs                    Ch                 Cp                  Tj
                                      │
                                      ├── [CZUJNIK] ──► kanał pomiarowy
                                      │       Tc
                                      └── Gw ──► wiązka / sąsiednia sekcja

$G_m$ obejmuje mocowanie skóry, $G_i$ — rzeczywisty interfejs obudowa–płytka, $R_j$ — drogę od elementu do złącza, a $G_w$ — odpływ przez wiązkę lub inną konstrukcję. Czujnik o temperaturze $T_c$ obserwuje własny węzeł i własną stałą czasową. Nawet bez błędu elektrycznego może więc reagować później niż $T_j$ albo pokazywać przeciwny kierunek krótkotrwałej zmiany, jeżeli źródło znajduje się po drugiej stronie sieci. Schemat jest autorski i funkcjonalny; nie przedstawia układu żadnego konkretnego pocisku.

Każdy węzeł powinien mieć fizyczne znaczenie, masę, własności, źródło mocy i ślad do konfiguracji. Każde połączenie potrzebuje materiału, geometrii lub uzasadnionego parametru styku. Warunki brzegowe obejmują nie tylko temperaturę powietrza, lecz także promieniowanie słoneczne, mocowanie do platformy, przepływ, ciśnienie i stan sąsiednich sekcji. ECSS Thermal Analysis Handbook porządkuje takie modele, przypadki, niepewności i raportowanie.²

Źródła zewnętrzne to klimat, słońce, sąsiednia platforma, nagrzewanie aerodynamiczne i gorąca struga napędu. Źródła wewnętrzne to bateria, przetwornice, komputer, tor radiowy (RF), głowica naprowadzająca, grzałki i urządzenia wykonawcze. Podręcznik programu rakiet sondażowych NASA dostarcza jawnego przykładu powiązania środowiska cieplnego, nagrzewania aerodynamicznego, montażu i sprawdzeń przed lotem, choć nie stanowi źródła wymagań dla uzbrojenia.⁷ Oprogramowanie ustala kolejność i czas aktywacji odbiorników. Zmiana częstotliwości procesora, okresu transmisji albo czasu pracy czujnika jest więc zmianą cieplną, mimo że rysunek mechaniczny pozostaje ten sam.

Temperatura staje się błędem funkcjonalnym

Model jest użyteczny dopiero po przełożeniu temperatury na funkcję. Wspólna macierz zapobiega tworzeniu oddzielnych katalogów dla podzespołów:

Węzeł	Wielkość funkcjonalna	Mechanizm wpływu temperatury	Kompensacja lub ograniczenie	Dowód
bateria	napięcie i dostępna moc	rezystancja wewnętrzna, samonagrzewanie, starzenie	przygotowanie cieplne, budżet obciążenia, nadzór stanu	próba obciążeniowa po reprezentatywnej historii
IMU	błąd stały, skala i ustawienie osi	dryf czujnika, gradient mocowania, opóźnienie czujnika	kalibracja dynamiczna i estymacja	próba funkcjonalna podczas zmian temperatury
detektor podczerwieni, w tym obrazujący (IR/IIR)	szum, przesunięcie poziomu, niejednorodność i czas gotowości	temperatura detektora, optyki i okna	kontrola stanu, korekcja i kwalifikowana sekwencja	scena optyczna połączona z profilem cieplnym
tor RF i radom	wzmocnienie, faza, strata i błąd osi	własności dielektryczne, odkształcenie, moc strat	kalibracja, rozpraszanie ciepła, sterowanie cyklem pracy	pomiar RF w stanie zimnym i gorącym
urządzenie wykonawcze	moment, prędkość i tarcie	rezystancja uzwojeń, lepkość, luzy, samonagrzewanie	dobór mechaniczny, ograniczenie mocy, diagnostyka	próba pod obciążeniem w temperaturze
elektronika	czas, moc obliczeniowa i trwałość	temperatura złącza, upływność, zmiana parametrów	droga złącze–obudowa–struktura i zarządzanie mocą	pomiar mocy i temperatury dla znanej wersji oprogramowania
polimery i spoiny	szczelność, sztywność i geometria	mięknięcie, relaksacja, różnica rozszerzalności, wilgoć	dobór procesu, kontrola powierzchni i geometrii	kupony, badanie elementu i próba zespołu

Wiersze optyczny i radiowy wymagają dwóch różnych rodzajów dowodu. NASA TM-82394 wyprowadza dla detektora piroelektrycznego odpowiedź cieplną i elektryczną oraz składowe szumu ograniczające czułość.¹⁶ Potwierdza to metodę wiązania temperatury, przewodności cieplnej i pasma z mierzalnym parametrem detektora, ale nie dostarcza charakterystyki innej technologii ani wojskowego układu obrazującego. Z kolei NASA TM-110321 opisuje bezpośredni pomiar jednokierunkowej straty transmisji próbki radomu w funkcji częstotliwości i kąta padania oraz osobny pomiar jej własności dielektrycznych.¹⁷ Raport nie bada temperatury, dlatego podpiera wybór wielkości wyjściowych, lecz nie dowodzi ich cieplnej zależności. Tę zależność trzeba wykazać przez pomiar tych samych wielkości w kontrolowanych stanach cieplnych albo przez zwalidowany model materiałowy.

Kompensacja programowa nie usuwa fizyki. Tabela kalibracyjna IMU nie naprawi nieobserwowanego gradientu ani histerezy. Korekcja detektora nie usunie zmiany transmisji rozgrzanego okna. Próg nadtemperaturowy chroni sprzęt tylko wtedy, gdy czujnik reprezentuje chroniony węzeł, ma znane opóźnienie, a reakcja nie tworzy groźniejszego stanu systemowego.

Materiały muszą być oceniane w czasie, wilgoci i cyklach. NASA-STD-6016 wskazuje w swojej domenie temperaturę pracy, zmęczenie cieplne i mechaniczne, temperaturę zeszklenia, niedopasowanie współczynników rozszerzalności oraz wilgoć jako powiązane przesłanki doboru.⁸ Podstawy przewodnictwa, dielektryków i starzenia materiałów elektrycznych opisuje Celiński,¹³ a zależność właściwości polimerów od temperatury, wilgotności i metody badania — Garbacz, Tor-Świątek i Samujło.¹⁴

Architektura cieplna jest kompromisem systemowym

Środki pasywne — rozpraszacze, kontrolowane styki, taśmy cieplne, powłoki, bariery, masa cieplna lub materiały przemiany fazowej — nie wymagają aktywnego sterowania, ale zależą od warunku początkowego i wykonania. Masa chłodna na początku pochłonie krótki impuls; ta sama masa po długim wygrzaniu nie daje tej rezerwy. Dobra ścieżka od procesora może ogrzać sąsiedni czujnik. Zmiana powłoki może pomóc w słońcu i zaszkodzić sygnaturze lub pracy anteny.

Środki aktywne — grzałki, termostaty, chłodnice termoelektryczne, wentylatory lub pompy — dodają zapotrzebowanie na energię, sterowanie, czujniki i tryby uszkodzeń. Chłodnica termoelektryczna przenosi ciepło i sama wydziela moc; bez drogi po stronie gorącej może podnieść temperaturę całego modułu. Grzałka wymaga analizy zaniku oraz niekontrolowanego załączenia, ale szczegółowe progi i reakcje pozostają elementem dokumentacji konkretnego systemu.

Konfliktów nie rozwiązuje maksymalizacja jednego marginesu. Priorytetem jest bezpieczeństwo, następnie dostępność wymaganych funkcji, ich dokładność i dopiero potem optymalizacja masy, energii oraz kosztu. Rozwiązanie poprawiające temperaturę jednego węzła jest poprawą systemu tylko wtedy, gdy nie odbiera marginesu funkcji bardziej krytycznej.

Studium I: gotowość po zimnym magazynowaniu

Rozważmy kontenerowy efektor, którego bateria, IMU i mechanizm wykonawczy mają różne pojemności cieplne i minimalne stany potrzebne do pełnej funkcji. Aby pokazać metodę bez ujawniania danych wyrobu, użyjmy temperatury znormalizowanej

$$\theta_i=\frac{T_i-T_c}{T_h-T_c}$$

oraz czasu bezwymiarowego $s=t/t_{ref}$. Dla węzła ogrzewanego przez otoczenie o $\theta_e=1$ i bez źródła wewnętrznego rozwiązanie równania pierwszego rzędu ma postać

$$\theta_i(s)=1-\exp\left(-\frac{G_i^*}{C_i^*}s\right).$$

Poniższe współczynniki są wyłącznie syntetycznym przykładem dydaktycznym:

Węzeł	$C_i^*$	$G_i^*$	minimalne $\theta_i$ dla funkcji	wynikowy $s_{got}$
bateria	1,60	0,35	0,40	2,34
IMU	0,80	0,55	0,60	1,33
mechanizm	1,20	0,25	0,50	3,33

Czas wynika z $s_{got}=-(C_i^*/G_i^*)\ln(1-\theta_{min})$. W tym modelu granicą gotowości nie jest bateria, lecz wolno nagrzewający się mechanizm. Gdyby sprawdzono wyłącznie napięcie źródła, system uznano by za gotowy zbyt wcześnie. Gdyby zaś czujnik umieszczono na szybciej reagującej obudowie, wynik pomiaru nie dowodziłby temperatury wnętrza mechanizmu.

Następny krok nie polega na przyjęciu tych liczb jako projektu, lecz na zastąpieniu ich danymi konkretnej konfiguracji, dodaniu sprzężeń między węzłami i mocy podczas sprawdzenia przedstartowego. Próba powinna zmierzyć zarówno warunki brzegowe i pobór mocy, jak i odpowiedź węzłów. Korelacja modelu musi następnie przewidzieć inny przebieg niż ten użyty do strojenia. Dopiero taki pakiet pozwala uzasadnić czas gotowości i wskazać podsystem ograniczający.

Studium II: przenoszenie lotnicze i powrót bez użycia

Pocisk podwieszony pod samolot przechodzi przez gorący postój, wznoszenie, chłodne otoczenie wysokościowe, nagrzewanie zależne od prędkości oraz pracę własnej elektroniki. Temperatura skóry może szybko podążać za lotem, masywny moduł zachowuje stan z postoju, a zasilana głowica naprowadzająca tworzy lokalne źródło. Po lądowaniu zimne powierzchnie spotykają się z cieplejszym i wilgotnym powietrzem. Ta sama misja zawiera więc ogrzewanie, chłodzenie, gradient, cykl zmęczeniowy i ryzyko kondensacji.

Nie wystarczy zapisać minimum i maksimum. Dla degradacji zależnej od czasu można prowadzić znormalizowaną dawkę cieplną

$$D_T=\int \max[0,\theta(t)-\theta_{ref}]\,ds,$$

oddzielnie od liczby oraz amplitudy cykli. Wielokrotne loty bez odpalenia zwiększają oba składniki, choć licznik użyć bojowych pozostaje zerowy. Historia powinna być związana z numerem egzemplarza i konfiguracją platformy.

Dawka nie zastępuje jednak bilansu. Rozdzielmy skórę $s$, masywny moduł $m$ i zasilany detektor $d$. W jednostkach bezwymiarowych ich sprzężenie można zapisać jako

$$\begin{aligned} C_s^*\dot\theta_s&=G_{se}^*(\theta_e-\theta_s)+G_{sm}^*(\theta_m-\theta_s)+q_a^*,\\ C_m^*\dot\theta_m&=G_{sm}^*(\theta_s-\theta_m)+G_{md}^*(\theta_d-\theta_m)+P_m^*,\\ C_d^*\dot\theta_d&=G_{md}^*(\theta_m-\theta_d)+P_d^*. \end{aligned}$$

Kropka oznacza tu pochodną względem czasu bezwymiarowego: $\dot\theta=d\theta/ds$. W dydaktycznym rachunku przyjęto $C^*=(0{,}4;1{,}6;0{,}3)$, $G_{se}^*=0{,}8$, $G_{sm}^*=0{,}25$, $G_{md}^*=0{,}2$ i wspólny stan początkowy $\theta=0{,}4$. Kolejne fazy mają syntetyczne warunki: postój $(\Delta s,\theta_e,q_a^*,P_m^*,P_d^*)=(1{,}0;0{,}8;0;0{,}05;0{,}10)$, lot $(1{,}2;-0{,}6;0{,}2;0{,}10;0{,}25)$ i powrót $(0{,}8;0{,}5;0;0{,}02;0)$. Całkowanie układu daje:

Koniec fazy	$\theta_s$ — skóra	$\theta_m$ — moduł	$\theta_d$ — detektor
postój	0,693	0,471	0,661
lot	-0,105	0,544	1,270
powrót	0,439	0,572	0,975

Minimum skóry po locie nie jest minimum modułu, a wyłączenie detektora przy powrocie nie usuwa zgromadzonej energii. Zewnętrzny kanał może już pokazywać ocieplenie, gdy źródło wewnętrzne wciąż pozostaje najgorętszym węzłem. Wartość $\theta_d>1$ nie oznacza przekroczenia limitu konkretnego urządzenia — skala $T_c$–$T_h$ służy tu tylko normalizacji. Studium nie zawiera danych istniejącego pocisku.

Dowód funkcjonalny łączy co najmniej cztery obserwacje: temperaturę okna lub radomu, temperaturę wewnętrznego czujnika, stan toru RF i czas. Jeżeli po lądowaniu pojawi się nietypowy odczyt, przyczyną może być rzeczywisty gradient, kondensacja, opóźnienie czujnika albo zmiana kalibracji. Decyzji o dalszym użyciu nie wolno wyprowadzić z ogólnego artykułu; podejmuje się ją według autoryzowanej dokumentacji, danych konkretnego egzemplarza i analizy resursowej. Studium pokazuje natomiast, dlaczego „nieodpalony” nie znaczy „nieobciążony”.

Studium III: długotrwałe źródło napędowe i zmiana oprogramowania

W pocisku manewrującym z miniaturowym silnikiem turboodrzutowym napęd działa na tyle długo, że sekcja gorąca, wlot, paliwo, struktura i elektronika tworzą sprzężoną sieć. Zmniejszanie skali zwiększa udział powierzchni w stosunku do objętości, dlatego wymiana ciepła i straty nie skalują się tak jak ciąg czy masa. Gieras omawia to zagadnienie wraz z materiałami wysokotemperaturowymi, chłodzeniem, CFD i metrologią małego napędu.¹²

Pierwszy model może mieć trzy węzły: sekcję napędową, strukturę pośrednią i moduł elektroniki, połączone z otoczeniem jako warunkiem brzegowym. Powinien zamknąć bilans oraz wskazać dominującą ścieżkę. Model przestrzenny dodaje lokalne maksimum temperatury i gradient mocowania. Obliczeniowa mechanika płynów (CFD) jest uzasadniona dopiero wtedy, gdy niepewność przepływu zmienia decyzję. Historyczne badanie swobodnego lotu konstrukcji szybowcowo-rakietowej pokazuje, że lokalny transfer ciepła silnie zależy od geometrii; jego wartości nie są jednak danymi wejściowymi współczesnego systemu.¹¹

Załóżmy teraz, że nowy obraz oprogramowania zwiększa średnią moc modułu o $\Delta P$. W liniowym przybliżeniu wzrost ustalonej temperatury wynosi

$$\Delta T_{ss}\approx \Delta P\,R_{eq},$$

a po czasie $t$ osiąga część $1-\exp(-t/\tau)$. Ta prosta relacja rozstrzyga, czy krótka misja zobaczy pełny skutek i czy potrzebny jest model przestrzenny. Zmiana „tylko programowa” wpływa na węzeł elektroniki, potem przez strukturę na sąsiednie wyposażenie, a ostatecznie na margines i kwalifikację. Jeżeli pobór mocy, wersja modelu i wynik próby nie są związane z tą samą wersją kompilacji, dowód przestaje być ciągły.

Od rachunku kontrolnego do modelu predykcyjnego

Hierarchia modeli powinna rosnąć wraz z pytaniem, nie z ambicją zespołu. Rachunek bilansowy sprawdza znak, skalę czasu i możliwy zakres odpowiedzi. Sieć skupiona identyfikuje stany przejściowe i dominujące ścieżki. Model elementów skończonych opisuje gradient oraz lokalną geometrię. CFD służy konwekcji i sprzężeniu przepływu, a analiza termostrukturalna przenosi pole temperatury na odkształcenie, naprężenie i ustawienie osi. Podstawę związku odkształceń termicznych z kompatybilnością, siłami i sztywnością przedstawia Misiak.¹⁵

Wzorcowy proces ma siedem kroków:

wymaganie i decyzja określają potrzebną wielkość wyjściową;
rachunek niskiego rzędu wyznacza skalę i dominanty;
parametry otrzymują źródła, zakresy i wersje;
model przestrzenny rozwija tylko obszary wpływające na decyzję;
plan instrumentacji rozróżnia konkurencyjne hipotezy;
próba służy korelacji w fizycznych granicach;
inny przypadek sprawdza zdolność predykcji wraz z niepewnością.

NASA-STD-7009 rozdziela weryfikację implementacji, walidację wobec rzeczywistości, zamierzony sposób użycia i wiarygodność modelu.¹⁰ Zbieżność procedury numerycznej nie gwarantuje poprawnej fizyki. Idealny styk, uśredniona moc lub zbyt bliska granica o stałej temperaturze mogą dać gładki, lecz błędny wynik. Niezależny rachunek energii pozostaje wartościowy nawet przy rozbudowanej symulacji.

Instrumentacja, korelacja i rozpoznawanie przyczyn

Czujnik mierzy własny stan po zamontowaniu, a nie abstrakcyjną temperaturę punktu. Warstwa kleju, przewody, emisyjność i bezwładność zmieniają odpowiedź. Termopara potrzebuje kompensacji złącza odniesienia, czujnik rezystancyjny — kontroli toru pomiarowego, a kamera podczerwieni — znanej emisyjności i odbić. Kalibracja powinna obejmować cały kanał.

Rozmieszczenie wynika z hipotez: lokalne maksimum temperatury, najzimniejszy węzeł, warunek brzegowy, styk o niepewnej przewodności i punkt ograniczający funkcję. Kanały temperatury trzeba zsynchronizować z mocą, stanem oprogramowania, ciśnieniem, prędkością i zdarzeniami. Bez wspólnego czasu nie wiadomo, czy wzrost spowodowała aktywacja toru RF, zmiana otoczenia czy błąd montażu czujnika.

Korelacja nie jest dowolnym dopasowaniem wielu oporów do wielu punktów. Modyfikuje tylko fizycznie uzasadnione parametry w znanych granicach, zachowuje wyniki przed zmianą i ma przewidzieć niezależny przypadek. GSFC-STD-7000 opisuje próbę bilansu cieplnego jako źródło danych o temperaturach i mocy do korelacji modelu używanego później do predykcji.⁵ Dla pocisku atmosferycznego stanowisko będzie inne, ale musi zachować tę samą kontrolę wejść.

Kwalifikacja: co można przenieść ze źródeł kosmicznych

Źródła NASA i ECSS dostarczają dojrzałego języka bilansu, modeli, korelacji, konfiguracji i rozdzielenia celów prób. Nie dostarczają gotowych poziomów dla uzbrojenia. Różnicę porządkuje tabela:

Element podejścia	Można przenieść	Trzeba zdefiniować dla pocisku
bilans i ciągłość energii	równania, kontrolę zamknięcia i niepewność	profil atmosferyczny, platformę, napęd i czas
wiarygodność modelu	weryfikację, walidację, zamierzony sposób użycia	krytyczne funkcje i kryteria akceptacji
próba bilansowa i cykle	rozdzielenie korelacji od wykazania trwałości	medium, gradienty, stany pracy i poziomy
konfiguracja dowodu	ślad model–próba–wyrób	wariant bojowy, oprogramowanie, partie i resurs
próżnia cieplna	metodę kontroli promieniowania i przewodzenia	uzasadnienie, czy brak konwekcji reprezentuje fazę misji

NASA-STD-7002 wiąże przypadki gorące i zimne z wymaganiami prób oraz marginesami w swojej domenie,⁴ ale próżnia nie jest automatycznie surowszym odpowiednikiem atmosfery. Pocisk wysokościowy zmienia ciśnienie w czasie, zaś wyrób magazynowany i przenoszony zwykle zależy od konwekcji. Równoważność trzeba wykazać przez temperatury i gradienty krytycznych węzłów, nie nazwę komory.

Jawne źródło Army Test and Evaluation Command domyka zakres domenowy. Aktywna procedura TOP-5-2-582 dotyczy wprost prób temperatura–wysokość systemu rakietowego i jego wyposażenia podsystemowego. Obejmuje zdolność pracy, przetrwanie ekspozycji, gwałtowne zmiany ciśnienia oraz wykrywanie degradacji, a środowisko wiąże z całym cyklem życia, w tym transportem i magazynowaniem bez opakowania ochronnego.¹⁸ Nie ustanawia przez to jednego uniwersalnego profilu dla każdego pocisku: testowany wyrób, fazy życia, funkcje i kryteria nadal muszą należeć do tej samej podstawy kwalifikacyjnej.

Raport DOT&E za rok fiskalny 2020 pokazuje praktyczną granicę takiego dowodu na przykładzie JAGM. Program zakończył badania środowiskowe w komorze, ale nie wykonał jeszcze prób lotnych w skrajnym zimnie; DOT&E wskazało, że reprezentatywne środowisko arktyczne może ujawnić ograniczenia niewidoczne w komorze statycznej.¹⁹ Nie przeczy to wartości próby laboratoryjnej. Oznacza, że kontrolowany test odpowiada na pytanie o zadane wymuszenie, a lot — o sprzężony przebieg środowiska, platformy i działającego systemu.

Próba kwalifikacyjna wykazuje projekt z odpowiednią rezerwą. Próba odbiorcza ma wychwycić wady egzemplarza bez niepotrzebnego zużywania trwałości. Badania łączone są potrzebne, gdy fizycznie uzasadniona interakcja — na przykład gorący klej pod drganiami lub zimne uszczelnienie w ruchu — może ujawnić inny mechanizm niż osobne próby. NASA-STD-5001 przypomina w swojej domenie o łącznym uwzględnianiu obciążeń mechanicznych i termicznych.⁹

Produkcja, oprogramowanie i dokładna konfiguracja

Konfiguracja cieplna obejmuje geometrię, masę, materiały, powłoki, styki, kleje, docisk, komponenty, tryby mocy, oprogramowanie układowe, programowalny układ logiczny (FPGA), tablice kalibracyjne, procedurę numeryczną, siatkę, przypadki oraz stanowisko badawcze. Zmiana koloru, śruby, układu płytki albo harmonogramu zadań może być termicznie istotna. „Ten sam kształt i funkcja” nie dowodzi podobieństwa, jeżeli zmieniły się pojemność, lokalne maksimum temperatury, przewodność styku lub moc.

Dowód przez podobieństwo powinien porównać wszystkie drogi oraz źródła i wskazać, który przypadek pozostaje ograniczający. Zmiana pozornie korzystna również wymaga sprawdzenia: lepsze odprowadzanie z procesora może podnieść temperaturę IMU, a lżejsza obudowa szybciej reaguje na krótki impuls.

Dane odbiorcze mają wartość w analizie trendu. Czas odpowiedzi, różnica temperatur, pobór mocy i wyznaczony opór można porównywać między egzemplarzami. Dryf mieszczący się jeszcze w prostym limicie może ujawnić zmianę procesu. Powtórzenie próby aż do wyniku pozytywnego bez zapisu niezgodności niszczy informację o rzeczywistym stanie wykonania.

Nadzór eksploatacyjny i rekonstrukcja anomalii

Długie magazynowanie zmienia baterie, czujniki, powłoki, styki, uszczelnienia i materiały interfejsowe. Program nadzoru powinien reprezentować partie, zakłady, wiek, mikroklimat i historię transportu. Świeży egzemplarz zapasowy nie opisuje automatycznie najstarszej populacji. Wyniki testów funkcjonalnych można uzupełniać inspekcją i badaniami niszczącymi próbek wycofanych.

Po anomalii trzeba najpierw zabezpieczyć surowe kanały temperatury i mocy, wspólny czas, dziennik zdarzeń, warunki otoczenia, historię wcześniejszych faz, numery seryjne, kalibracje, model i fotografie. Szybki wzrost temperatury może oznaczać większą moc, utratę styku, zmianę warunku brzegowego lub błąd czujnika. Poprawienie jednego współczynnika tak, aby wykres się zgadzał, nie rozstrzyga przyczyny.

Wspólnym narzędziem jest reszta bilansu dla węzła $i$:

$$r_i(t)=C_i\dot T_i-\sum_jG_{ij}(T_j-T_i)-P_i-\dot Q_{b,i}.$$

Składnik $\dot Q_{b,i}$ zbiera wymianę z warunkami brzegowymi, przy zachowaniu tej samej konwencji znaków co w modelu. W poprawnym zakresie modelu $r_i$ powinno być zgodne z zerem w granicach niepewności. Niezerowa reszta nie wskazuje automatycznie uszkodzonej części, ponieważ błąd pojemności, styku, mocy, warunku brzegowego lub czujnika może dać podobny wykres. Hipotezę wybiera się przez dodatkową obserwację, która powinna zmienić się tylko dla części mechanizmów:

Hipoteza	Wzorzec w danych	Próba rozcinająca	Co ją osłabia
wzrost mocy	zgodny wzrost prądu poprzedza temperaturę; odpowiedź przechodzi znanymi drogami	odtworzenie obrazu oprogramowania na bezpiecznym emulatorze obciążenia	brak zmiany energii wejściowej przy poprawnej synchronizacji
utrata przewodności styku	rośnie różnica złącze–obudowa i zmienia się stała czasowa przy podobnej mocy	mały, kontrolowany impuls cieplny oraz niezależny pomiar po obu stronach styku	jednakowa odpowiedź przed i po ponownym kontrolowanym montażu
zły warunek brzegowy	kilka węzłów odchyla się wspólnie i koreluje z otoczeniem lub uchwytem	dodatkowy pomiar otoczenia, mocowania i strumienia na granicy	anomalia ograniczona do jednego kanału wewnętrznego
błąd toru pomiarowego	skok jednego kanału nie zamyka się energetycznie i nie ma odpowiedzi sąsiadów	wzorzec kalibracyjny albo drugi czujnik o innej zasadzie działania	zgodna odpowiedź elektryczna, cieplna i funkcjonalna kilku kanałów

Próba rozcinająca ma być bezpieczna i możliwie niskoenergetyczna. Nieznanego uzbrojenia nie uruchamia się w celu „sprawdzenia wykresu”; rekonstrukcja wykorzystuje zatwierdzone stanowisko, emulator albo próbkę materiałową właściwą dla hipotezy.

Wartością modelu po anomalii jest przewidywanie niezależnej obserwacji. Hipoteza utraty przewodności styku powinna wyjaśnić nie tylko temperaturę procesora, lecz także zmianę opóźnienia i odpowiedź sąsiedniej obudowy. Hipoteza błędnej mocy powinna znaleźć potwierdzenie w prądzie lub stanie oprogramowania. Decyzję o przedłużeniu resursu podejmuje się na podstawie spójności całego pakietu, a nie pojedynczego kanału.

Jak oceniać dowód cieplny

Rzetelny raport pozwala przejść w obie strony między wymaganiem a dowodem. Powinien wskazywać wariant i egzemplarz, profil życia, miejsca i rodzaje limitów, warunki brzegowe, historię mocy, wersję modelu, źródła parametrów, niepewność, margines, konfigurację próby i zakres korelacji. Wykres temperatur potrzebuje osi czasu z fazami i stanami systemu; mapa potrzebuje skali, geometrii i warunków brzegowych.

Najbardziej użyteczne pytanie kontrolne brzmi nie „czy pocisk wytrzyma daną temperaturę?”, lecz: który węzeł i która funkcja, po jakiej historii, przy których źródłach i drogach przepływu, z jaką niepewnością oraz na podstawie jakiej reprezentatywnej próby?

Trzy studia pokazują konsekwencję tej zasady. W efekcie kontenerowym czas gotowości ograniczył inny węzeł niż pierwszy sprawdzany. W przenoszeniu lotniczym brak odpalenia nie oznaczał braku zużycia cieplnego. W pocisku z długotrwałym napędem zmiana oprogramowania przeniosła się przez moc i strukturę aż do ważności kwalifikacji. We wszystkich przypadkach o wartości analizy decydowała ciągłość śladu, a nie liczba obliczeń.

Źródła

Przypisy

European Cooperation for Space Standardization, ECSS-E-ST-31C — Thermal control, 2008, https://ecss.nl/standard/ecss-e-st-31c-thermal-control/, dostęp: 2026-09-04; domena kosmiczna.
↩
European Cooperation for Space Standardization, ECSS-E-HB-31-03A — Thermal analysis handbook, 2016, https://ecss.nl/home/ecss-e-hb-31-03a-15november2016/, dostęp: 2026-09-04.
↩
NASA Small Spacecraft Systems Virtual Institute, State-of-the-Art of Small Spacecraft Technology — Thermal Control, edycja 2026, https://www.nasa.gov/smallsat-institute/sst-soa/thermal-control/, dostęp: 2026-09-04, oficjalny przegląd bilansu oraz pasywnych i aktywnych technologii; część danych bazuje na informacjach producentów.
↩
NASA, Payload Test Requirements, NASA-STD-7002B w/Change 1, https://standards.nasa.gov/standard/NASA/NASA-STD-7002, dostęp: 2026-09-04.
↩
NASA Goddard Space Flight Center, General Environmental Verification Standard, GSFC-STD-7000B, 2021, https://standards.nasa.gov/sites/default/files/standards/GSFC/B/0/gsfc-std-7000b_signature_cycle_04_28_2021_fixed_links.pdf, dostęp: 2026-09-04; domena programów GSFC.
↩
U.S. Department of Defense / DLA ASSIST, Environmental Engineering Considerations and Laboratory Tests, MIL-STD-810H Change 1, https://quicksearch.dla.mil/qsDocDetails.aspx?ident_number=35978, dostęp: 2026-09-04.
↩
NASA Sounding Rocket Program, NASA Sounding Rockets User Handbook, NASA/TP-20230006855, 2023, https://ntrs.nasa.gov/citations/20230006855, dostęp: 2026-09-04; rakiety sondażowe, nie uzbrojenie.
↩
NASA, Standard Materials and Processes Requirements for Spacecraft, NASA-STD-6016C w/Change 1, https://standards.nasa.gov/standard/nasa/nasa-std-6016, dostęp: 2026-09-04.
↩
NASA, Structural Design and Test Factors of Safety for Spaceflight Hardware, NASA-STD-5001B w/Change 3, https://standards.nasa.gov/standard/NASA/NASA-STD-5001, dostęp: 2026-09-04; wymaga dopasowania.
↩
NASA, Standard for Models and Simulations, NASA-STD-7009B, 2024, https://standards.nasa.gov/standard/NASA/NASA-STD-7009, dostęp: 2026-09-10.
↩
NASA/NACA, Free-Flight Investigation of Aerodynamic Heat Transfer to a Simulated Glide-Rocket Shape at Mach Numbers up to 10, rekord NTRS 19630010650, https://ntrs.nasa.gov/citations/19630010650, dostęp: 2026-09-04; bez przenoszenia wartości na aktywne systemy.
↩
Marian Gieras, Miniaturowe silniki turboodrzutowe, 2025, rozdz. 3–6.
↩
Zdzisław Celiński, Materiałoznawstwo elektrotechniczne, 2018, rozdz. 1–5.
↩
Tomasz Garbacz, Aneta Tor-Świątek i Bronisław Samujło, Właściwości mechaniczne i cieplne tworzyw polimerowych. Ćwiczenia laboratoryjne, 2017; nie jest normą projektową.
↩
Jan Misiak, Wytrzymałość materiałów z przykładami obliczeń, 2025, rozdziały o odkształceniach, układach niewyznaczalnych, dynamice i powłokach; nie jest normą pocisku.
↩
R. L. Kroes, D. Reiss, The Pyroelectric Properties of TGS for Application in Infrared Detection, NASA-TM-82394, 1981, https://ntrs.nasa.gov/citations/19810010326, dostęp: 2026-09-10; nie jest charakterystyką głowicy naprowadzającej.
↩
Robin L. Cravey, Pacita I. Tiemsin, W-Band Transmission Measurements and X-Band Dielectric Properties Measurements for a Radome Material Sample, NASA-TM-110321, 1997, https://ntrs.nasa.gov/citations/19970017316, dostęp: 2026-09-10; bez zależności temperaturowej i bez przenoszenia wartości na radomy pocisków.
↩
U.S. Army Test and Evaluation Command, TOP-5-2-582 — Temperature-Altitude Tests for Missile System and Missile Subsystem Equipment, 2 czerwca 2011, https://quicksearch.dla.mil/qsDocDetails.aspx?ident_number=278816, dostęp: 2026-09-10; poziomy i kryteria wymagają dopasowania do badanego wyrobu.
↩
Director, Operational Test & Evaluation, FY 2020 Annual Report, część „Joint Air-to-Ground Missile (JAGM)”, s. 96, 2021, https://www.dote.osd.mil/Portals/97/pub/reports/FY2020/other/2020DOTEAnnualReport.pdf, dostęp: 2026-09-10.
↩