Budżet masy, środek masy i tensor bezwładności pocisku

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

Właściwości masowe jako kontrolowany rekord

Masa, położenie środka masy i tensor bezwładności bywają przedstawiane jako kilka pól w tabeli. Taki skrót jest wygodny, lecz zaciera najważniejszy fakt: liczby nie opisują typu wyrobu ani nawet numeru części, tylko określony egzemplarz lub model w określonej konfiguracji, chwili i układzie odniesienia. Ta sama masa całkowita może należeć do obiektów o innym środku masy, osiach głównych i odpowiedzi obrotowej. Z kolei ten sam egzemplarz zmienia właściwości po zużyciu materiału, przemieszczeniu mechanizmu, separacji, naprawie albo wymianie modułu.

Teza monografii brzmi zatem: właściwości masowe nie są zestawem stałych parametrów, lecz kontrolowanym rekordem modelu fizycznego, którego wiarygodność powstaje w zamkniętym łańcuchu wymaganie → budżet → model rozkładu → pomiar → niepewność → decyzja zgodności → utrzymanie konfiguracji. Każde przerwanie tego łańcucha może pozostawić liczbę poprawną rachunkowo, ale nieprzydatną inżyniersko.

Tekst dotyczy publicznych metod mechaniki, metrologii i zarządzania konfiguracją. Wszystkie studia są jawnie syntetyczne, abstrakcyjne i nieoperacyjne: nie odtwarzają geometrii, tolerancji ani charakterystyk żadnego istniejącego uzbrojenia. Nieznanego, uszkodzonego lub odnalezionego obiektu wojskowego nie wolno podnosić, obracać, zawieszać, opróżniać ani ustawiać na stanowisku. Miejsce zabezpiecza się i przekazuje uprawnionemu zespołowi EOD.

pominięto rzeczywiste rozmieszczenie podzespołów, stany materiałów energetycznych, tolerancje odbiorcze, parametry sterowania, obciążenia oraz procedury prób współczesnych systemów bojowych. W dokumentacji programu miejsca te należy uzupełnić danymi kontrolowanymi i zatwierdzonymi procedurami; przykładów z tego artykułu nie wolno traktować jako zamiennika.

Trzy poziomy opisu jednej bryły

Najprostszy poziom opisuje ruch postępowy. Masa

$$m=\int_V \rho(\mathbf r)\,dV$$

wiąże wypadkową siłę z przyspieszeniem środka masy. Nie mówi jednak, gdzie masa się znajduje. Tę informację pierwszego rzędu wnosi wektor środka masy

$$\mathbf r_C=\frac{1}{m}\int_V \mathbf r\rho(\mathbf r)\,dV, \qquad \mathbf r_C=\frac{\sum_i m_i\mathbf r_i}{\sum_i m_i}$$

dla odpowiednio opisu ciągłego i dyskretnego. Iloczyny $m_i\mathbf r_i$ są pierwszymi momentami masy. Zależność wygląda elementarnie, lecz działa tylko wtedy, gdy wszystkie położenia mają wspólne datum, osie, jednostki i znaki.

Trzeci poziom opisuje rozkład wokół środka. Przyjmijmy prawoskrętny układ $(x,y,z)$ oraz macierz

$$\mathbf I_C= \begin{bmatrix} I_{xx}&I_{xy}&I_{xz}\\ I_{xy}&I_{yy}&I_{yz}\\ I_{xz}&I_{yz}&I_{zz} \end{bmatrix},$$

gdzie używana tutaj konwencja definiuje elementy pozadiagonalne jako $I_{xy}=-\int xy\,dm$, analogicznie dla pozostałych par. W literaturze spotyka się również przechowywanie dodatnich całek iloczynowych i wpisywanie minusów dopiero do macierzy. Nazwa pola Ixy bez definicji nie rozstrzyga konwencji.

Tensor jest symetryczny, ma więc sześć niezależnych składowych. Wiąże prędkość kątową z momentem pędu: $\mathbf H_C=\mathbf I_C\boldsymbol\omega$. Energia obrotu wynosi

$$T_r=\frac12\boldsymbol\omega^T\mathbf I_C\boldsymbol\omega.$$

Dla dowolnej prędkości energia nie może być ujemna, więc ogólnym warunkiem fizyczności tensora jest dodatnia półokreśloność: $\boldsymbol\omega^T\mathbf I_C\boldsymbol\omega\ge0$. Zerowa wartość dla niezerowego $\boldsymbol\omega$ jest możliwa w idealizacji zdegenerowanej, na przykład dla masy punktowej albo idealnej linii obracanej wokół własnej osi. Dodatnia określoność obowiązuje dla niezdegenerowanej bryły przestrzennej o niezerowym rozkładzie masy poza każdą osią przechodzącą przez środek. Nieujemne wartości własne, symetria i nierówności trójkątne między momentami głównymi są użytecznymi kontrolami fizyczności; w przypadku niezdegenerowanym wartości własne muszą być dodatnie. Kontrole te nie dowodzą kompletności modelu.

Znaczenie produktów bezwładności widać w równaniu Eulera dla bryły sztywnej:

$$\mathbf M_C=\mathbf I_C\dot{\boldsymbol\omega}+ \boldsymbol\omega\times(\mathbf I_C\boldsymbol\omega).$$

Jeżeli układ osi nie jest główny, moment zadany nominalnie wokół jednej osi może wywołać odpowiedź sprzężoną. Model z samymi $I_{xx}$, $I_{yy}$ i $I_{zz}$ milcząco zakłada brak tej informacji. Założenie może być dopuszczalne po analizie czułości, ale nie powinno powstawać przez wyzerowanie niewygodnych pól.

Układ odniesienia jest częścią wartości

Rekord powinien definiować prawoskrętny układ bryły, dodatnie kierunki osi, punkt początkowy, jednostki i regułę transformacji. Datum musi być odtwarzalne na wyrobie, w modelu CAD i na stanowisku. „Czubek” lub środek rysunku nie jest dobrym datum, jeśli element jest wymienny albo globalny początek złożenia może zostać przesunięty przez reorganizację modelu.

Wektor punktu przechodzi między układami przez obrót i translację. Dla przyjętej konwencji aktywnego obrotu tensora:

$$\mathbf I'=\mathbf R\mathbf I\mathbf R^T.$$

Czysty obrót zachowuje ślad, wyznacznik i wartości własne tensora. Jest to znakomity test regresji transformacji. Nie zachowuje natomiast poszczególnych elementów macierzy; pojawienie się produktu bezwładności po obrocie bryły diagonalnej jest poprawnym wynikiem, a nie błędem numerycznym.

Przeniesienie tensora między równoległymi osiami opisuje twierdzenie Steinera. Dla wektora $\mathbf d$ od środka masy komponentu do nowego punktu:

$$\mathbf I_O=\mathbf I_C+m\left[(\mathbf d\cdot\mathbf d)\mathbf E- \mathbf d\mathbf d^T\right].$$

Obrót i translacja są różnymi operacjami. Bezpieczny algorytm złożenia najpierw obraca lokalny tensor każdego komponentu do wspólnego układu, następnie wyznacza systemowy środek masy, raz stosuje twierdzenie Steinera i sumuje wkłady:

$$m=\sum_i m_i, \quad \mathbf r_C=\frac{\sum_i m_i\mathbf r_i}{m}, \quad \mathbf I_C=\sum_i\left[ \mathbf R_i\mathbf I_{i,c}\mathbf R_i^T+ m_i\big((\mathbf d_i\cdot\mathbf d_i)\mathbf E-\mathbf d_i\mathbf d_i^T\big) \right].$$

Zastąpienie rozległego komponentu punktem zachowuje masę i środek, ale usuwa jego lokalny tensor. Model skupiony bywa wystarczający dla pierwszego budżetu; nie jest przez to automatycznie odpowiedni do analizy obrotu, drgań lub separacji.

Osie główne i interpretacja dynamiczna

Diagonalizacja $\mathbf I_C\mathbf q_j=\lambda_j\mathbf q_j$ daje trzy główne momenty $\lambda_j$ i odpowiadające im prostopadłe osie $\mathbf q_j$. W bazie złożonej z wektorów własnych produkty bezwładności znikają. Osie te są właściwością rozkładu masy, nie geometrii zewnętrznej. Wiązka, złącze, korekta wyważenia lub niejednakowe aktuatory mogą je obrócić względem osi montażowej.

Gdy dwa momenty główne są niemal równe, kierunek osi w ich płaszczyźnie może być bardzo czuły na małą zmianę danych. Nie oznacza to równie dużej zmiany fizycznej odpowiedzi: wartości własne mogą pozostać stabilne. Dlatego rekord powinien raportować zarówno macierz w ustalonym układzie bryły, jak i wartości oraz orientację osi głównych z właściwą niepewnością. Sam kąt osi bez informacji o rozdzieleniu wartości własnych jest podatny na błędną interpretację.

Promień bezwładności $k=\sqrt{I/m}$ pomaga porównywać rozkłady, ale zawsze wymaga wskazania osi i punktu. Wartość względem datum nie jest wymienna z wartością względem środka masy. Podobnie „środek ciężkości” można utożsamiać ze środkiem masy jedynie w przybliżeniu jednorodnego pola grawitacyjnego; w rekordzie dynamiki przechowuje się wektor środka masy.

Konfiguracja i czas

Najkrótszy poprawny zapis brzmi $\mathcal P=\mathcal P(C,t,s)$: właściwości zależą od konfiguracji $C$, czasu $t$ i stanu $s$. Konfiguracja magazynowa może zawierać osłony i zabezpieczenia transportowe; startowa — wszystkie elementy obecne przy zwolnieniu z platformy; stan po separacji — już nie. Nie wolno przenosić wartości między nimi przez zmianę etykiety.

Materiały zużywalne tworzą historię masy, środka i tensora, a nie tylko malejącą masę skalarną. Element ruchomy dodaje do stanu położenie mechanizmu. Ciecz może wymagać modelu swobodnej powierzchni lub drgań zamiast jednej bryły sztywnej. Każde uproszczenie ma zadeklarowany zakres częstotliwości, przyspieszeń i decyzji, dla których pozostaje użyteczne.

Separacja wymaga dwóch pełnych rekordów — przed zdarzeniem i po nim — ponieważ centrum całego pozostałego złożenia trzeba obliczyć ponownie. Odjęcie tensora odrzuconej części od macierzy zapisanej względem starego środka nie daje macierzy względem nowego środka. Dynamika samego przejścia obejmuje dodatkowo warunki początkowe i oddziaływania między ciałami; nie są one właściwościami masowymi.

Wymaganie, obwiednia i reguła decyzji

Wymaganie musi wskazywać wielkość, konfigurację, datum, układ, punkt odniesienia tensora, jednostki, fazę, obwiednię oraz metodę weryfikacji. „Środek przed stacją X” nie określa znaku, współrzędnych poprzecznych ani sposobu uwzględnienia niepewności. „Bezwładność nominalna” nie mówi, czy ogranicza się element macierzy w układzie bryły, moment główny czy orientację osi.

Obwiednia wielowymiarowa nie powinna powstawać z osobnych minimów i maksimów każdego pola. Połączenie najbardziej przedniego środka z największymi sześcioma elementami tensora może nie należeć do żadnej fizycznej konfiguracji, a nawet naruszać dodatnią półokreśloność lub nierówności trójkątne momentów głównych. Do analiz przekazuje się katalog realizowalnych przypadków, trajektorie zależne od czasu albo wspólny model probabilistyczny zachowujący korelacje.

Kryterium zgodności obejmuje regułę decyzji zatwierdzoną przed poznaniem wyniku. JCGM 106 wiąże specyfikację, wynik pomiaru, niepewność oraz ryzyka błędnej akceptacji i odrzucenia; opisuje też strefy akceptacji zawężane pasem ochronnym.¹³ Pas ochronny (ang. guard band) zmniejsza ryzyko fałszywej akceptacji kosztem większego ryzyka fałszywego odrzucenia. Niepewności nie dodaje się po teście w kierunku wygodnym dla decyzji. NASA-STD-7002B pozostawia szczegółowe wymagania właściwości masowych zależne od misji i łączy ocenę analityczną z pomiarem; nie istnieje jedna publiczna tolerancja odpowiednia dla każdego obiektu.²

Budżet masy jest budżetem wiedzy

Limit systemowy rozdziela się według struktury podziału pracy (WBS), lecz granice WBS nie zawsze pokrywają się ze strukturą złożenia CAD. Potrzebne jest mapowanie oraz właściciel każdego składnika, aby śruba interfejsowa, wiązka lub klej nie zniknęły między zespołami albo nie zostały policzone dwa razy. Lista wyposażenia połączona z bazą właściwości masowych prowadzi od stanu projektowanego do zbudowanego i utrzymywanego.

Należy rozdzielić pojęcia:

alokacja jest przydzielonym fragmentem limitu, nie prognozą;
bieżące najlepsze oszacowanie obejmuje całą znaną treść projektu i podaje pochodzenie danych;
rezerwa na wzrost pokrywa przewidywany przyrost wynikający z niedojrzałości, a jej wartość zależy od klasy danych;
margines jest zdefiniowaną różnicą względem ograniczenia technicznego;
rezerwa zarządcza pozostaje zasobem programu, przydzielanym kontrolowaną decyzją.

Jedno ryzyko nie może równocześnie zwiększać prognozy, rezerwy na wzrost i rezerwy programu. Znany wspornik bez końcowego rysunku jest treścią projektu, a nie „nieprzewidzianą rezerwą”. Klasa dojrzałości rozróżnia analogię, model parametryczny, wstępny CAD, zwolniony rysunek, część wyprodukowaną i wynik pomiaru konkretnej sztuki. Etykieta „zmierzono” obowiązuje tylko dla konfiguracji objętej pomiarem.

Trend powinien zachowywać przyczynę każdego skoku: zmianę fizyczną, poprawkę błędu, zmianę zakresu albo zmianę algorytmu. Gładki wykres nie jest celem. Wczesne ujawnienie przyrostu umożliwia analizę wariantów; ukrycie go do przeglądu odbiera budżetowi funkcję ostrzegawczą. NASA Systems Engineering Handbook traktuje margines masy jako techniczną miarę postępu i wiąże wymagania, konfigurację oraz dowód realizacji.⁵

Studium I: ewolucja syntetycznego obiektu dwuczłonowego

To studium nie opisuje rakiety ani skali rzeczywistego systemu. Cztery abstrakcyjne bryły służą wyłącznie do sprawdzenia matematyki rekordu.

Wymaganie i dane

Wymaganie dydaktyczne R-1 stanowi: dla trzech stanów — początkowego S0, po umownym zużyciu S1 i po odłączeniu tylnego członu S2 — należy wyznaczyć masę, trójwymiarowy środek, pełny tensor w środku masy, momenty główne i odpowiedź pierwszego rzędu na ten sam abstrakcyjny moment. Macierz ma być symetryczna i dodatnio określona. Decyzja nie dotyczy przydatności bojowej; potwierdza jedynie wewnętrzną zgodność modelu demonstracyjnego.

Układ jest prawoskrętny, odległości podano w metrach, masy w kilogramach, a tensory w kg·m². Dane S0 są następujące:

Bryła	$m$	$\mathbf r_i=(x,y,z)$	lokalny tensor centralny	orientacja
A, przednia	12	$(-0{,}80;0{,}02;0)$	diag$(0{,}08;0{,}55;0{,}56)$	obrót $30°$ wokół z
B, środkowa	28	$(0;0;0)$	diag$(0{,}24;3{,}20;3{,}20)$	zgodna z układem
C, tylna	20	$(1{,}00;-0{,}03;0{,}01)$	diag$(0{,}18;1{,}40;1{,}45)$	obrót $20°$ wokół x
D, zużywalna	10	$(0{,}45;0;0)$	diag$(0{,}05;0{,}30;0{,}30)$	zgodna z układem

Model i pełny rachunek S0

Dla bryły A macierz obrotu ma postać

$$\mathbf R_z(30°)= \begin{bmatrix}0{,}8660&-0{,}5000&0\\0{,}5000&0{,}8660&0\\0&0&1\end{bmatrix}.$$

Po obrocie jej tensor centralny wynosi

$$\mathbf I_{A,c}^{B}=\mathbf R_z\mathbf I_{A,c}^{A}\mathbf R_z^T= \begin{bmatrix} 0{,}1975&-0{,}2035&0\\ -0{,}2035&0{,}4325&0\\ 0&0&0{,}5600 \end{bmatrix}.$$

Znak elementu pozadiagonalnego wynika z zadeklarowanej konwencji. Analogiczny obrót bryły C wokół osi x wprowadza człon $I_{yz}$. Suma mas daje $m_0=70$, a suma pierwszych momentów $(14{,}9;-0{,}36;0{,}20)$, zatem

$$\mathbf r_{C,0}=(0{,}21286;-0{,}00514;0{,}00286).$$

Dla A wektor do wspólnego środka to $\mathbf d_A=(-1{,}01286;0{,}02514;-0{,}00286)$. Jej pełny wkład jest równy

$$\mathbf I_{A@C_0}=\mathbf I_{A,c}^{B}+ 12[(\mathbf d_A\cdot\mathbf d_A)\mathbf E-\mathbf d_A\mathbf d_A^T].$$

Taką samą operację wykonuje się raz dla każdej bryły. Suma czterech wkładów daje

$$\mathbf I_{C,0}= \begin{bmatrix} 0{,}68988&0{,}51186&-0{,}15743\\ 0{,}51186&31{,}87321&-0{,}01110\\ -0{,}15743&-0{,}01110&32{,}05853 \end{bmatrix}.$$

Wartości własne wynoszą $(0{,}68069;31{,}88056;32{,}06036)$. Oś najmniejszego momentu jest skierowana w przybliżeniu jak $(0{,}99985;-0{,}01641;0{,}00501)$, czyli nie pokrywa się dokładnie z osią x. Wszystkie wartości są dodatnie; ślad macierzy $64{,}62162$ zgadza się z sumą momentów głównych w granicy zaokrąglenia.

Zmiana stanu, wynik syntetyczny i decyzja

W S1 bryłę D zastępuje abstrakcyjny stan o masie 4, środku $(0{,}25;0;0)$ i tensorze diag$(0{,}02;0{,}12;0{,}12)$. Ponowne złożenie — nie skalarne odjęcie — daje:

Stan	masa	środek masy	główne momenty
S0	70	$(0{,}21286;-0{,}00514;0{,}00286)$	$(0{,}68069;31{,}88056;32{,}06036)$
S1	64	$(0{,}17813;-0{,}00563;0{,}00313)$	$(0{,}64973;31{,}06704;31{,}24629)$
S2	44	$(-0{,}19545;0{,}00545;0)$	$(0{,}46063;10{,}00195;10{,}13258)$

S2 powstaje przez usunięcie C ze stanu S1 i kolejne wyznaczenie środka. Macierze potrzebne do niezależnego odtworzenia odpowiedzi w dwóch zmienionych stanach mają postać

$$\mathbf I_{C,1}=\begin{bmatrix} 0{,}65965&0{,}52436&-0{,}16438\\ 0{,}52436&31{,}05910&-0{,}01119\\ -0{,}16438&-0{,}01119&31{,}24430 \end{bmatrix}, \qquad \mathbf I_{C,2}=\begin{bmatrix} 0{,}46099&-0{,}05843&0\\ -0{,}05843&10{,}00159&0\\ 0&0&10{,}13258 \end{bmatrix}.$$

Dla ilustracji interfejsu z dynamiką przyłóżmy w chwili $\boldsymbol\omega=\mathbf0$ ten sam umowny moment $\mathbf M=(0;10;0)$ N·m. Pełne równanie Eulera ma postać

$$\dot{\boldsymbol\omega}=\mathbf I_C^{-1} \left[\mathbf M-\boldsymbol\omega\times(\mathbf I_C\boldsymbol\omega)\right],$$

a w rozpatrywanej chwili redukuje się dokładnie, nie w przybliżeniu skalarnym, do $\dot{\boldsymbol\omega}=\mathbf I_C^{-1}\mathbf M$. Niezależne rozwiązanie trzech układów liniowych daje:

Stan	$\dot{\boldsymbol\omega}=(\dot\omega_x;\dot\omega_y;\dot\omega_z)$ `rad/s²`	kontrola $\mathbf I_C\dot{\boldsymbol\omega}$, N·m
S0	$(-0{,}235831;0{,}317530;-0{,}001048)$	$(0{,}000000;10{,}000000;0{,}000000)$
S1	$(-0{,}259730;0{,}326351;-0{,}001250)$	$(0{,}000000;10{,}000000;0{,}000000)$
S2	$(0{,}126812;1{,}000582;0{,}000000)$	$(0{,}000000;10{,}000000;0{,}000000)$

Porównanie dowodzi dwóch rzeczy ograniczonych do tego syntetycznego modelu i chwili początkowej: zmiana konfiguracji zmienia odpowiedź na ten sam moment, a produkty bezwładności wywołują składowe przyspieszenia poza osią y i mogą nawet zmienić ich znak. Nie dowodzi przebiegu późniejszego ruchu, stateczności ani jakości układu sterowania; dla $\boldsymbol\omega\ne\mathbf0$ trzeba zachować człon żyroskopowy i całkować pełne równania ruchu.

Niepewność i regresja

W tym studium nie wykonano pomiaru. Syntetyczny model przypisuje masom standardową niepewność względną $0{,}2\%$, położeniom $0{,}5$ mm, a elementom lokalnych tensorów $1\%$; wejścia uznano demonstracyjnie za nieskorelowane. Propagacja Monte Carlo jest częścią rekordu, lecz jej wynik nie stanowi tolerancji żadnego wyrobu. Decyzja: wszystkie trzy stany przyjęto do dalszego ćwiczenia, ponieważ spełniają R-1 i zachowują fizyczność macierzy.

Regresja obejmuje: sumę mas i pierwszych momentów, symetrię, dodatnią półokreśloność oraz — dla zadeklarowanych tu niezdegenerowanych brył — dodatnią określoność, zgodność śladu z sumą wartości własnych, niezmienniczość wartości własnych po czystym obrocie, odwrócenie transformacji oraz niezależne przeliczenie płaskiej listy. Zmiana jednego wejścia wymaga ponownego obliczenia wszystkich dotkniętych stanów i porównania delty z analitycznym pierwszym momentem.

Źródła danych i dojrzewanie modelu

Na początku programu dostępne są analogie i modele parametryczne, później dochodzą modele CAD, dane dostawców, pomiary części i złożeń, a w końcu wynik konkretnego egzemplarza. Nie tworzy to automatycznej hierarchii „pomiar zawsze lepszy”. Pomiar niewłaściwej konfiguracji lub z nieznanym tare może być słabszym dowodem niż zweryfikowany model. Rekord opisuje metodę, zakres, pochodzenie i stosowalność.

Model CAD jest całką po geometrii i przypisanej gęstości. Wiele cyfr w eksporcie nie zmniejsza niepewności gęstości, grubości powłoki ani brakujących elementów. Kontrola kompletności obejmuje bryły ukryte, zdublowane, wiązki, elementy złączne, kleje, uszczelniacze, materiały procesowe, media, pozycje mechanizmów i sprzęt tylko naziemny. Gęstość zastępcza powinna różnić się statusem od wartości z certyfikatu materiałowego.

Dane dostawcy mogą być katalogowe, wstępne, rysunkowe, certyfikowane dla partii albo zmierzone dla sztuki. Numer części nie rozstrzyga, czy masa obejmuje złącza, osprzęt montażowy i opcje. Rozległy moduł wymaga oprócz masy lokalnego środka, tensora i orientacji. Jeśli tych danych brak, model jawnie opisuje przybliżenie i jego wpływ.

Model elementów skończonych powinien zachować masę, środek i tensor, ale także sensowną dystrybucję względem ścieżek obciążeń i modów. Skupienie całej masy wyposażenia w jednym węźle może zamknąć trzy agregaty, a mimo to zafałszować lokalne obciążenie. NASA wymaga kontroli modeli masowych używanych w analizach obciążeń; szczegółowe kryteria z domeny kosmicznej wymagają dostosowania do programu.³^,⁷

Model pomiaru i niepewność

Pomiar zaczyna się od definicji wielkości mierzonej: identyfikatora sztuki, konfiguracji, warunków, czasu, układu, datum i oczekiwanych poprawek. „Masa obiektu” jest opisem potocznym, nie pełną definicją. Model pomiaru łączy wskazania z wynikiem przez wzorcowanie, tare, geometrię, ustawienie, dryf, temperaturę, wypór i oddziaływanie przewodów.

Ocena typu A pochodzi ze statystycznej analizy obserwacji. Seria odczytów bez zdjęcia i ponownego ustawienia obiektu ocenia głównie szum, nie powtarzalność całej procedury. Ocena typu B wykorzystuje świadectwa wzorcowania, rozdzielczość, wcześniejsze dane, tolerancje i wiedzę o stanowisku. Granicy $\pm a$ nie dzieli się mechanicznie przez $\sqrt3$; ten przelicznik wymaga uzasadnionego rozkładu prostokątnego.¹⁵

Dla modelu $y=f(\mathbf x)$ liniowe przybliżenie niepewności ma postać

$$u_c^2(y)=\mathbf c^T\mathbf U_x\mathbf c, \qquad c_i=\frac{\partial f}{\partial x_i},$$

gdzie $\mathbf U_x$ jest macierzą kowariancji wejść. Pierwiastek z sumy kwadratów jest tylko przypadkiem wejść nieskorelowanych. Wyniki wielowymiarowe, takie jak trzy współrzędne środka i sześć składowych tensora, również powinny zachować kowariancję. JCGM 100 opisuje podstawowe ramy GUM, JCGM 101 propagację rozkładów metodą Monte Carlo, a JCGM 102 wiele wielkości wyjściowych.¹⁰^,¹¹^,¹²

Niepewność rozszerzona $U=ku_c$ wymaga podania $k$ i interpretacji pokrycia. Zapis ± bez tej informacji jest niepełny. Rozpoznany błąd systematyczny należy skorygować albo usunąć jego przyczynę; niepewność poprawki pozostaje. Nie wolno chować znanego obciążenia przewodu lub przesunięcia datum w ogólnym „błędzie aparatury”.

Przy silnej nieliniowości, ograniczeniach fizycznych lub asymetrycznych rozkładach propagacja Monte Carlo może być bardziej wiarygodna niż linearyzacja. Losowanie musi jednak zachowywać wspólne błędy wagi, partii i procesu. Niezależne losowanie składników ze wspólnym przesunięciem sztucznie je uśredni. Liczba prób, kryterium zbieżności i polityka ziarna generatora są częścią odtwarzalności.

Ważenie, środek masy i tensor

Ważenie wielopodporowe wykorzystuje sumę reakcji i równowagę momentów. Trzy niekoliniowe podpory wyznaczają rzut środka na ich płaszczyznę; brakującą współrzędną uzyskuje się po zmianie orientacji albo inną metodą. Cztery podpory tworzą układ zależny od sztywności i poziomowania, dlatego większa liczba czujników nie gwarantuje prostszego modelu.

Metoda odwrócenia jest diagnostyczna: rzeczywisty mały mimośród zmienia znak po obrocie o 180°, podczas gdy część przesunięcia kanałów pozostaje. Metoda pochylenia wymaga zmierzonego kąta, rzeczywistych punktów podparcia i kontroli ugięcia. Zawieszenie lub balansowanie na krawędzi może być intuicyjne dla małych obojętnych brył, ale nie jest domyślnie bezpieczną ani dokładną metodą dla dużego sprzętu.

Bezwładność wyznacza się pośrednio z odpowiedzi dynamicznej. Wahadło skrętne wykorzystuje sztywność skrętną, bifilarne i trifilarne — zmianę energii potencjalnej przy wychyleniu, a wahadło fizyczne — znaną odległość osi od środka. Metoda reakcyjna korzysta ze znanego momentu i przyspieszenia. Każdy wynik dotyczy całego poruszającego się zespołu; właściwość obiektu powstaje dopiero po odjęciu przyrządu i transformacji. Podręcznik NASA dla modeli aeroelastycznych opisuje wahadła bi- i trifilarne, odejmowanie bezwładności przyrządu, kalibrację wzorcami oraz konieczność sztywnego połączenia, umieszczenia środka na osi i wyeliminowania ruchu innego niż obrót.¹⁴

Produkty bezwładności można identyfikować z pomiarów względem kilku osi. Zestaw orientacji musi zapewnić obserwowalność, a nadmiarowe orientacje umożliwiają estymację najmniejszych kwadratów i ocenę reszt. Mały produkt jest często różnicą większych wielkości, więc wynik bliski zeru może mieć dużą niepewność względną. Nie należy go zerować bez wykazania nieistotności. NASA TR R-433 omawia doświadczalne wyznaczanie momentów, produktów i nachylenia osi głównej oraz wskazuje jako warunki metody wahadłowej między innymi ograniczenie modów kołysania przez projekt zawieszenia i oddzielenie modów pasożytniczych w danych.¹⁶

Studium II: syntetyczna kampania obojętnego artykułu testowego

Studium jest rachunkiem laboratoryjnym dla abstrakcyjnej, nieenergetycznej bryły. Nie jest instrukcją budowy stanowiska. Rzeczywista próba wymaga zatwierdzonego planu podnoszenia, osprzętu o potwierdzonej nośności, ograniczenia ruchu i przeglądu gotowości.

Wymaganie, dane i model

Wymaganie M-1: dla konfiguracji IT-01/C2 wyznaczyć masę z $U\le0{,}10$ kg, każdą współrzędną środka z $U\le8$ mm oraz moment względem osi $z$ przy $U\le0{,}06$ kg·m²; $k=2$. Akceptacja jest możliwa, jeśli przedział ochronny mieści masę w $49{,}7$–$50{,}3$ kg, $x$ w $0{,}45$–$0{,}51$ m, $|y|<0{,}06$ m, $|z|<0{,}09$ m, a $I_{zz}$ w $1{,}85$–$2{,}15$ kg·m². Granice są wymyślone wyłącznie dla ćwiczenia decyzji.

Trzy podpory leżą w punktach $P_1=(0;-0{,}4)$, $P_2=(0;0{,}4)$ i $P_3=(1{,}2;0)$ m. Po poprawce tare syntetyczne średnie reakcji wyrażone jako równoważne masy wynoszą $(12{,}50;17{,}50;20{,}00)$ kg. Model statyczny daje

$$m=R_1+R_2+R_3=50{,}00,$$

$$x_C=\frac{1{,}2R_3}{m}=0{,}4800,$$

$$y_C=\frac{-0{,}4R_1+0{,}4R_2}{m}=0{,}0400 \text{m}.$$

Po kontrolowanym pochyleniu o $\theta=20°$ zmierzony rzut wynosi $x_\theta=0{,}47328$ m. Dla przyjętej geometrii $x_\theta=x_C\cos\theta+z_C\sin\theta$, więc

$$z_C=\frac{x_\theta-x_C\cos\theta}{\sin\theta}=0{,}0650 \text{m}.$$

Pomiar dynamiczny i tensor

Wahadło mierzy wyłącznie moment wokół pionowej osi $z_p$ przechodzącej przez środek platformy $O_p$. Przed serią obiekt przetłumaczono na platformie tak, aby jego środek miał względem tej osi położenie $\mathbf r_{C/O_p}=(0;0;0{,}120)$ m: rzut środka leży na $z_p$, więc poprawka Steinera do osi równoległej wynosi $m(d_x^2+d_y^2)=0$. Wysokość $0{,}120$ m nie zmienia momentu wokół tej samej prostej. Niepewność resztkowego mimośrodu i nieidealnego ruchu mieści jawny składnik modelowy poniżej.

Syntetyczna platforma trifilarna ma promień zawieszeń $r=0{,}3000$ m i długość $L=0{,}8000$ m. Dla małych oscylacji, sztywnych linek i zespołu poruszającego się jak jedna bryła model daje

$$I=\frac{mgr^2T^2}{4\pi^2L}.$$

Zespół platforma–obiekt o masie $m_f=58{,}00$ kg i okresie $T_f=1{,}1392745$ s ma $I_{z,pełne}=2{,}10449$ kg·m². Pusta platforma o masie $m_p=8{,}00$ kg i okresie $T_p=0{,}65000$ s ma $I_{z,przyrząd}=0{,}09449$ kg·m². Po odjęciu przyrządu wynik względem osi $z_p$ — pokrywającej się z centralną osią z obiektu — wynosi

$$I_{zz}=\frac{gr^2}{4\pi^2L}(m_fT_f^2-m_pT_p^2)=2{,}01000 \text{kg·m}^2.$$

Seria nadmiarowych orientacji i ważone dopasowanie dostarczają pozostałych pięciu niezależnych składowych; wraz z powyższym $I_{zz}$ dają centralny tensor

$$\widehat{\mathbf I}_C= \begin{bmatrix} 0{,}420&0{,}018&-0{,}012\\ 0{,}018&1{,}870&0{,}026\\ -0{,}012&0{,}026&2{,}010 \end{bmatrix} \text{kg·m}^2.$$

Jego momenty główne, obliczone z pokazanej macierzy, to $(0{,}41968;1{,}86560;2{,}01472)$. Oś największego momentu jest nachylona w płaszczyźnie y–z o około $10{,}2°$ względem z; wynik pokazuje, że małe produkty mogą nieznacznie obracać osie bez utraty dodatniej określoności.

Liczbowy budżet niepewności

Dla masy przyjęto standardowe składniki: powtarzalność pełnego ustawienia 0,020 kg, wzorcowanie wspólnej skali 0,025 kg, tare 0,015 kg, rozdzielczość 0,006 kg i warunki środowiskowe 0,010 kg. Przy założeniach studium

$$u_c(m)=\sqrt{0{,}020^2+0{,}025^2+0{,}015^2+0{,}006^2+0{,}010^2}=0{,}0374 \text{kg},$$

zatem $u_c(m)=0{,}0374$ kg i $U(m)=0{,}075$ kg dla $k=2$.

Dla współrzędnych reakcje opisano jako $R_i=(1+s)\bar R_i+\varepsilon_i$, gdzie wspólny względny błąd skali ma $u(s)=0{,}0005$, a niezależne reszty kanałów po pełnym zdjęciu i ponownym ustawieniu mają $u(\varepsilon_i)=0{,}010$ kg. Zatem

$$\begin{aligned} \operatorname{cov}(R_i,R_j) &=\bar R_i\bar R_j u^2(s)\\ &\quad+\delta_{ij}u^2(\varepsilon_i). \end{aligned}$$

Kowariancje od wspólnej skali nie są zerowe, ale ich wkład do ilorazów $x_C$ i $y_C$ znosi się, ponieważ $\sum_i(\partial x_C/\partial R_i)R_i=\sum_i(\partial y_C/\partial R_i)R_i=0$. Reszty $\varepsilon_i$ uznano za niezależne, ponieważ pochodzą z osobnych kanałów i osobnych estymat po wspólnym resecie; lokalne błędy współrzędnych podpór są niezależne na podstawie osobnych odczytów, natomiast wspólne błędy datum zachowano jako wspólne zmienne. Dane i współczynniki wrażliwości są następujące:

Wynik	wejście i standardowa niepewność	współczynnik wrażliwości
$x_C$	$R_1,R_2,R_3$: po $0{,}010$ kg	$(-0{,}0096;-0{,}0096;0{,}0144)$ m/kg
$x_C$	położenie $x_3=1{,}2000$ m, $u=0{,}5$ mm	$R_3/m=0{,}4$
$x_C$	wspólne datum x, $u=0{,}5$ mm	$1$
$y_C$	$R_1,R_2,R_3$: po $0{,}010$ kg	$(-0{,}0088;0{,}0072;-0{,}0008)$ m/kg
$y_C$	lokalne $y_1,y_2$, każde $u=0{,}5$ mm	$(R_1/m;R_2/m)=(0{,}25;0{,}35)$
$y_C$	wspólne datum y, $u=1{,}175$ mm	$1$

Po podstawieniu $u_c(x_C)=0{,}574$ mm i $u_c(y_C)=1{,}200$ mm, a więc $U(x_C)=1{,}15$ mm i $U(y_C)=2{,}40$ mm.

Dla $z_C$ przyjęto $u(x_\theta)=1{,}0$ mm, wyznaczone wyżej $u(x_C)=0{,}574$ mm i $u(\theta)=0{,}05°=8{,}727\cdot10^{-4}$ rad. Wspólne datum x o $u=0{,}5$ mm daje $\operatorname{cov}(x_\theta,x_C)=2{,}50\cdot10^{-7}$ m²; pozostałe składniki tych dwóch estymat pochodzą z odrębnych ustawień i są niezależne. W punkcie pomiarowym współczynniki wynoszą

$$\frac{\partial z}{\partial x_\theta}=2{,}92380,$$

$$\frac{\partial z}{\partial x_C}=-2{,}74748,$$

$$\frac{\partial z}{\partial\theta}=0{,}30144 \text{m/rad}.$$

Wariancja w mm² wynosi $8{,}5486+2{,}4840-4{,}0165+0{,}0692=7{,}0853$, gdzie ujemny wyraz jest podwójnym członem kowariancyjnym. Stąd $u_c(z_C)=2{,}662$ mm i $U(z_C)=5{,}32$ mm.

Dla wyniku wahadła model wejść brzmi $I_{zz}=A(m_fT_f^2-m_pT_p^2)+b_I$, gdzie $A=gr^2/(4\pi^2L)$, a $b_I$ skupia resztkową poprawkę modelu. Pełny budżet jest odtwarzalny z tabeli:

Wejście	estymata	$u$	$\partial I/\partial q$
$m_f$	$58{,}00$ kg	$0{,}04173$ kg	$0{,}036284$ m²
$m_p$	$8{,}00$ kg	$0{,}01552$ kg	$-0{,}011811$ m²
$T_f$	$1{,}1392745$ s	$0{,}003507$ s	$3{,}69444$ kg·m²/s
$T_p$	$0{,}65000$ s	$0{,}002004$ s	$-0{,}290734$ kg·m²/s
$r$	$0{,}3000$ m	$0{,}0005$ m	$13{,}4000$ kg·m
$L$	$0{,}8000$ m	$0{,}0010$ m	$-2{,}5125$ kg·m
$g$	$9{,}810$ m/s²	$0{,}005$ m/s²	$0{,}204893$ kg·m·s²
$b_I$	$0$ `kg·m²`	$0{,}0170$ `kg·m²`	$1$

Niezerowe kowariancje wynoszą $\operatorname{cov}(m_f,m_p)=1{,}16\cdot10^{-4}$ kg² od wspólnej skali masy, $\operatorname{cov}(T_f,T_p)=2{,}962\cdot10^{-8}$ s² od wspólnej skali czasu oraz $\operatorname{cov}(r,L)=2{,}00\cdot10^{-7}$ m², czyli $\rho=0{,}4$, od wspólnego wzorcowania geometrii. Standardowe niepewności mas obejmują odpowiednio niezależne składniki $0{,}030$ i $0{,}015$ kg oraz wspólną względną skalę $0{,}0005$; okresy — niezależne składniki $0{,}0035$ i $0{,}0020$ s oraz wspólną względną skalę $0{,}0002$. $b_I$ obejmuje łącznie granice linearyzacji małego kąta, resztkowy mimośród, ruch pasożytniczy i odjęcie przyrządu. Pozostałe kowariancje przyjęto jako zerowe, ponieważ wielkości pochodzą z niezależnych torów — masowego, czasowego, geometrycznego i lokalnego $g$ — a w ramach każdego toru wspólne skale zapisano jawnie.

Wkłady do $u_c^2(I_{zz})$ od mas, czasu, geometrii, $g$ i $b_I$ wynoszą kolejno $(2{,}2263;168{,}1830;37{,}7357;1{,}0495;289{,}0000)\cdot10^{-6}$ (kg·m²)². Suma $4{,}98195\cdot10^{-4}$ daje $u_c(I_{zz})=0{,}02232$ i $U(I_{zz})=0{,}04464$ kg·m² dla $k=2$.

Decyzja i regresja

Raportuje się $m=(50{,}000\pm0{,}075)$ kg, $\mathbf r_C=(0{,}4800\pm0{,}00115;0{,}0400\pm0{,}00240;0{,}0650\pm0{,}00532)$ m oraz $I_{zz}=(2{,}01000\pm0{,}04464)$ kg·m², wszędzie dla $k=2$. Zastosowana wcześniej reguła wymaga, aby cały symetryczny przedział $[v-U,v+U]$ leżał w granicach. Przedziały wynoszą odpowiednio: masa $[49{,}925;50{,}075]$ kg, x $[0{,}47885;0{,}48115]$ m, y $[0{,}03760;0{,}04240]$ m, z $[0{,}05968;0{,}07032]$ m i $I_{zz}$ $[1{,}96536;2{,}05464]$ kg·m². Każdy mieści się w syntetycznej granicy M-1, a wszystkie wymagane wartości $U$ są nie większe od limitów, więc konfigurację IT-01/C2 akceptuje się wyłącznie dla zadeklarowanego kontekstu użycia. Niepewność niektórych produktów macierzy jest większa względnie; nie wykorzystano ich do tej decyzji bez osobnej kwalifikacji.

Regresja kampanii obejmuje powtórne ustawienie, odwrócenie orientacji, bilans sumy reakcji, zamknięcie momentów, wzorzec kontrolny przed i po serii, analizę reszt okresu względem amplitudy, dodatnią określoność tensora, kontrolę surowych danych drugim algorytmem i zgodność identyfikatora konfiguracji. Ponowny test jest uzasadniony anomalią albo zmianą, nie wybieraniem korzystniejszego wyniku.

Redukcja danych i zamknięcie korelacji

Surowe dane zachowują oryginalne wskazania, czas, status kanału, współczynniki wzorcowania i mapowanie czujników. Filtr, metoda estymacji okresu, poprawki oraz wersja biblioteki numerycznej należą do konfiguracji redukcji. Filtr może przesunąć fazę lub stłumić amplitudę; zgodny wykres po filtracji nie zastępuje kontroli sygnału wejściowego.

Oprogramowanie redukcji wymaga testów jednostkowych dla brył o rozwiązaniu analitycznym, transformacji, jednostek i przypadków granicznych. Drugie uruchomienie tego samego arkusza nie jest kontrolą niezależną. Warto użyć prostszego, osobno utworzonego bilansu reakcji lub drugiej implementacji. Reszty dopasowania powinny być zgodne z modelem szumu i nie wykazywać struktury zależnej od orientacji, amplitudy ani czasu.

Korelacja z CAD zaczyna się od ujednolicenia konfiguracji, datum, układu, punktu i jednostek. Dopiero potem porównuje się masę, pierwsze momenty, tensor i sumy częściowe. Dopasowanie jednolitej gęstości całego modelu może zamknąć masę, jednocześnie pogarszając rozkład. Korekta ma fizyczną przyczynę — brakującą powłokę, zmierzoną gęstość, pozycję lub błąd transformacji — oraz pozostałą niepewność.

Wiarygodność modelu zależy od kontekstu użycia. Model wystarczający do wczesnego bilansu nie musi nadawać się do przypadku kwalifikacyjnego. NASA-STD-7009B opisuje cykl życia, weryfikację, walidację, niepewność i komunikowanie podstaw wiarygodności modeli; podręcznik NASA-HDBK-7009B rozwija wdrożenie tych zasad.⁸^,⁹ Etykieta „zwalidowany” nie przechodzi automatycznie na inną konfigurację ani decyzję.

Wspólny interfejs analiz systemowych

Jedna zmiana rozkładu masy rozchodzi się do wielu odbiorców. W modelu sześciu stopni swobody masa wpływa na translację, środek na ramiona sił, a tensor na odpowiedź obrotową. Algorytm sterowania korzystający ze starych danych może przewidzieć inną odpowiedź niż aktualna konfiguracja, nawet jeśli jego implementacja nie uległa zmianie.

W analizie obciążeń rozkład określa siły bezwładności, ścinanie i zginanie. W analizie drgań wpływa na częstotliwości, postacie modów i udział modalny. W aeroelastyczności liczy się położenie mas względem osi sprężystej. Dla separacji potrzebne są spójne rekordy obu ciał. Dla platformy, transportu i obsługi mogą obowiązywać odrębne złożenia z kontenerem, adapterem lub uchwytami.

Odbiorca danych zapisuje identyfikator dokładnego wydania, nie kopiuje liczb do lokalnego arkusza z ręczną poprawką. Po zmianie bazowej system powinien wskazać analizy zależne: dynamikę, obciążenia, drgania, separację, wyważenie i interfejs platformy. Automatyczny alert inicjuje ocenę; nie zastępuje decyzji inżynierskiej.

Studium III: zamiana modułu bez zmiany masy całkowitej

Studium dotyczy fikcyjnego złożenia ZM-50, nie odpowiada żadnemu wariantowi uzbrojenia. Pokazuje, dlaczego stwierdzenie „masa bez zmian” nie zamyka wniosku o zmianę konfiguracji.

Wymaganie i dane zmiany

Wymaganie Z-1: zamienny moduł musi zachować masę całego złożenia $50{,}00\pm0{,}10$ kg; środek poprzeczny ma leżeć w obwiedni $|y|,|z|<0{,}015$ m; wszystkie główne momenty mają pozostać w zakresie modelu 6DOF-Z3; a zmiana dowolnego produktu bezwładności powyżej $0{,}10$ kg·m² wymaga jawnej ponownej oceny analiz zależnych. Granice są syntetyczne i nie są zaleceniami projektowymi.

Rdzeń o masie 42 kg ma środek $(0;0;0)$ i centralny tensor diag$(0{,}35;5{,}20;5{,}25)$. Stary moduł M-A ma masę 8 kg, środek $(0{,}60;0{,}08;0)$ i tensor diag$(0{,}06;0{,}20;0{,}22)$. Dostawca proponuje M-B o tej samej masie i obrysie, lecz ze środkiem $(0{,}60;-0{,}02;0{,}04)$, tensorze lokalnym diag$(0{,}10;0{,}16;0{,}24)$ i obrocie lokalnych osi o $20°$ wokół z. Różnica wynika z fikcyjnego przemieszczenia ciężkich elementów wewnętrznych.

Model, obrót, Steiner i osie główne

Dla M-B najpierw wykonuje się obrót

$$\mathbf I_{B,c}^{sys}=\mathbf R_z(20°) \operatorname{diag}(0{,}10;0{,}16;0{,}24)\mathbf R_z^T(20°),$$

a dopiero potem translację do systemowego środka. Złożenie stare daje:

$$m_A=50,\qquad \mathbf r_{C,A}=(0{,}0960;0{,}0128;0),$$

$$\mathbf I_{C,A}= \begin{bmatrix} 0{,}45301&-0{,}32256&0\\ -0{,}32256&7{,}81920&0\\ 0&0&7{,}93221 \end{bmatrix}.$$

Jego momenty główne wynoszą $(0{,}43891;7{,}83330;7{,}93221)$, a najmniejsza oś jest odchylona od x o około $2{,}50°$ w płaszczyźnie x–y. Dla nowej konfiguracji:

$$m_B=50,\qquad \mathbf r_{C,B}=(0{,}0960;-0{,}0032;0{,}0064),$$

$$\mathbf I_{C,B}= \begin{bmatrix} 0{,}47046&0{,}06136&-0{,}16128\\ 0{,}06136&7{,}78293&0{,}00538\\ -0{,}16128&0{,}00538&7{,}91189 \end{bmatrix}.$$

Momenty główne zmieniają się na $(0{,}46645;7{,}78333;7{,}91551)$. Najmniejsza oś ma w przybliżeniu kierunek $(0{,}99973;-0{,}00840;0{,}02166)$. Masa całkowita i współrzędna x środka pozostają identyczne, ale $y$ zmienia się o $-16{,}0$ mm, $z$ o $+6{,}4$ mm, $I_{xy}$ o $+0{,}38392$ i pojawia się $I_{xz}=-0{,}16128$ kg·m².

Pomiar syntetyczny, niepewność i wpływ

Odbiorczy wynik syntetyczny dla złożenia z M-B wynosi $m=(50{,}02\pm0{,}08)$ kg, $y_C=(-3{,}0\pm2{,}0)$ mm, $z_C=(6{,}8\pm2{,}4)$ mm, $I_{xy}=(0{,}064\pm0{,}030)$ i $I_{xz}=(-0{,}158\pm0{,}034)$ kg·m², $k=2$. Pomiar jest zgodny z modelem w obrębie łącznej niepewności; nie „kalibruje się” modelu do środka każdego przedziału.

Wpływ przechodzi przez następujący łańcuch:

rekord dostawcy zmienia lokalny środek, tensor i transformację;
złożenie CAD wyznacza nowy środek systemowy i macierz przez obrót oraz twierdzenie Steinera;
model 6DOF otrzymuje pełną macierz, a nie adnotację „0 kg zmiany”;
analizy obciążeń i drgań sprawdzają zmianę rozkładu, nie tylko całkowitej masy;
proces wyważenia ocenia przesunięcie środka oraz obrót osi głównych;
rada kontroli zmian rozpatruje wyniki i aktualizuje stosowalność;
pomiar odbiorczy zamyka korelację konkretnego egzemplarza.

Decyzja konfiguracji i regresja

Stary stan spełniał granicę $|y|<15$ mm tylko z małym zapasem po uwzględnieniu niepewności; nowy stan spełnia obwiednię poprzeczną. Jednak zmiany $I_{xy}$ i $I_{xz}$ przekraczają próg ponownej oceny. Decyzja rady brzmi: nie akceptować zamiany wyłącznie na podstawie równej masy; warunkowo dopuścić M-B po zakończeniu wskazanych analiz i wydaniu nowej bazy ZM-50/B2. Jest to decyzja konfiguracji, nie ocena zdolności systemu bojowego.

Regresja obejmuje porównanie rekordu dostawcy z pomiarem modułu, test śladu i wartości własnych po obrocie, niezależne złożenie systemu, różnicę pierwszych momentów, dodatnią określoność, uruchomienie niezmienionych przypadków referencyjnych 6DOF i obciążeń oraz kontrolę listy odbiorców danych. Wyniki starej bazy pozostają odtwarzalne; nie nadpisuje się ich nowymi liczbami.

Produkcja, wyważanie i utrzymanie

Egzemplarz rozwojowy, kwalifikacyjny, lotny, obojętny i symulator masowy nie są wymienne. Balast może odtworzyć masę i środek, ale nie tensor, sztywność ani lokalną ścieżkę obciążenia. Specyfikacja artykułu testowego powinna wskazywać, które właściwości odwzorowuje i dla jakiej próby.

Zmienność produkcyjna obejmuje tolerancje geometrii, materiału, powłok i montażu. Dwie sztuki o tej samej masie mogą mieć różne środki przez kompensujące się różnice. Składniki z jednej partii lub procesu bywają skorelowane. Plan odbioru może łączyć ważenie każdej sztuki, kontrolę położeń, okresowy pomiar bezwładności i wyważenie, lecz musi wynikać z krytyczności, zdolności procesu i wykrywalności niezgodności.

Wyważenie statyczne sprowadza środek do wymaganej osi. Wyważenie dynamiczne ustawia oś główną; środek na osi nie gwarantuje braku pary niewyważenia. Masa korekcyjna ma kontrolowane położenie, mocowanie i stosowalność. Zmiana pokrywy, śruby albo przewodu po końcowym wyważeniu wymaga oceny wpływu i czasem ponownego zamknięcia.¹

Naprawa, przeróbka i zamiana części aktualizują rekord utrzymywany dla numeru seryjnego. Zgodność formy, montażu i funkcji nie oznacza identycznego rozkładu masy. Zmiany wynikające z zaniku źródeł dostaw kumulują się, dlatego okresowy pełny bilans kontroluje efekt wielu małych decyzji.

Podczas magazynowania utrata medium, korozja, wilgoć lub degradacja uszczelnień może zmieniać masę i jednocześnie sygnalizować stan techniczny. Nadzór eksploatacyjny trenduje wyniki z odniesieniem do sztuki, partii, stanowiska i temperatury. Kontrola przed użyciem może sprawdzać kompletność i nieoczekiwaną zmianę bez powtarzania kwalifikacyjnego pomiaru całego tensora.

Cyfrowy łańcuch i minimalny rekord

Cyfrowy łańcuch nie musi być jedną bazą. Musi natomiast zapewniać stabilne identyfikatory i kontrolowane przejścia między wymaganiem, WBS, listą materiałową, CAD, danymi dostawcy, modelem analitycznym, wynikiem próby i decyzją. Automatyzacja ogranicza przepisywanie, lecz równie sprawnie rozprowadza systematyczny błąd osi lub jednostek. Potrzebuje walidacji schematu, testów zachowania sum i niezależnego przeglądu zmian algorytmu.

Minimalne wydanie właściwości masowych zawiera:

Pole	Znaczenie kontrolne
identyfikator i wersja	jednoznaczne odtworzenie wydania, bez nadpisywania historii
stosowalność	wariant, numer seryjny lub partia oraz stan konfiguracji
czas lub zdarzenie	chwila, przedział albo stan przed/po zmianie
masa i środek	wartość, jednostki, datum i układ
tensor	sześć składowych, punkt odniesienia i konwencja znaków
osie główne	wartości własne, orientacja i informacja o degeneracji
metoda i pochodzenie	pomiar, obliczenie, dostawca, analogia lub założenie
niepewność	model, standardowa lub rozszerzona wartość, $k$, kowariancja
kontekst użycia	decyzje, dla których rekord ma wystarczającą wiarygodność
źródła wejściowe	identyfikatory plików, przyrządów, wzorcowań i kodu
status	roboczy, przejrzany, wydany, zastąpiony albo nieważny
zatwierdzenia	właściciel danych, przegląd i rada zmian

Słownik danych musi rozróżniać zero zmierzone, zero założone, brak pomiaru, wartość nieznaną i pole nie mające zastosowania. Hash wykrywa zmianę bajtów, ale nie zastępuje wersji semantycznej, zatwierdzenia ani informacji, dlaczego wynik się zmienił. Analiza zależna zapisuje identyfikator użytego rekordu i własną konfigurację; tylko wtedy można ustalić zakres ponownej oceny.

Ślad audytowy zachowuje poprzednią i nową wartość, autora, powód, zmianę zakresu, zatwierdzenie i listę dotkniętych konfiguracji. Archiwum obejmuje surowe dane, świadectwa, procedurę, zapis przebiegu próby, wejścia, kod, środowisko obliczeniowe i wydany raport. Otwarty, opisany format tabelaryczny pomaga w długim cyklu życia, lecz nie znosi kontroli dostępu. Dane masowe rzeczywistego systemu mogą ujawniać jego konfigurację; integralność, klasyfikacja, kopie zapasowe i autoryzacja są częścią zarządzania.

Odpowiedzialność i bramka zmiany

Integrator właściwości masowych utrzymuje metodę i bilans systemowy. Właściciele konstrukcji odpowiadają za kompletność elementów, metrologia za ważność wyniku, zarządzanie konfiguracją za stosowalność, a analitycy za użycie właściwego wydania. Jedna osoba nie musi pełnić wszystkich ról; przekazanie odpowiedzialności nie może jednak pozostawiać tabeli bez właściciela.

Każda zmiana inżynierska odpowiada co najmniej na pięć pytań: co fizycznie się zmienia; których stanów dotyczy; jaka jest delta masy, pierwszych momentów i tensora; które analizy i próby zależą od rekordu; jaki dowód zamyka zmianę. Wartość progowa może ustalać skalę oceny, ale zerowa delta masy nigdy nie jest automatycznym zwolnieniem z pozostałych pytań.

Niezgodność opisuje wymaganie, wynik z niepewnością, konfigurację i działanie zabezpieczające. Decyzja „użyć w stanie istniejącym” wymaga systemowej oceny wpływu. Przeróbka przywracająca definicję może dodać klej lub usunąć powłokę; naprawa wprowadza nową konfigurację; wymiana części może zmienić lokalny tensor. Każda z nich ma rzeczywistą deltę i sposób zamknięcia.

Przegląd okresowy nie powinien ograniczać się do masy całkowitej. Pokazuje trend, największych kontrybutorów, dojrzałość, niepewność, margines, otwarte niezgodności, naprawy i świeżość danych. Najmniejszy margines systemowy może dotyczyć poprzecznego środka lub momentu po zdarzeniu, choć całkowita masa wygląda komfortowo.

Bezpieczeństwo badań

Duża bryła na podporach, zawieszeniu albo urządzeniu wirującym ma energię potencjalną lub kinetyczną. Zagrożenia obejmują upadek, przewrócenie, zakleszczenie, zerwanie zawiesia, niekontrolowany obrót i przeciążenie czujnika. Stanowisko wymaga właściwej procedury, oceny zagrożeń, osprzętu o potwierdzonej nośności, strefy wyłączenia, ograniczników, kontroli zmagazynowanej energii i jednoznacznego warunku zatrzymania.

Świadectwo wzorcowania nie dowodzi przydatności przyrządu do danej próby. Należy ocenić zakres, rodzaj obciążenia, rozdzielczość, niepewność, środowisko i czas. W układzie wielopodporowym każda podpora musi mieć zapas dla nierównomiernego rozkładu. Mechaniczny ogranicznik nie może niewidocznie przejąć części reakcji.

Przegląd gotowości potwierdza konfigurację egzemplarza, przyrządy, dokumenty, personel, bezpieczeństwo, dane oczekiwane i uprawnienie do rozpoczęcia. Zapis przebiegu próby obejmuje odstępstwa oraz anomalie w chwili ich wystąpienia. Kontynuowanie dla „dokończenia serii” nie ma pierwszeństwa przed ochroną ludzi i sprzętu.

Skondensowana kontrola wydania

Przed przekazaniem rekordu należy sprawdzić pięć zamknięć:

Fizyczne: lista elementów, materiały zużywalne, mechanizmy, sprzęt naziemny i stan egzemplarza odpowiadają deklarowanej konfiguracji.
Matematyczne: wspólne datum, układ i jednostki; jawna konwencja produktów; jeden obrót i jedna właściwa translacja; symetria, dodatnia półokreśloność, nierówności momentów głównych oraz niezmienniki; dodatnia określoność tylko dla zadeklarowanej niezdegenerowanej bryły przestrzennej.
Metrologiczne: zdefiniowana wielkość mierzona; spójność pomiarowa; tare, geometria, ustawienie i środowisko; niepewność, kowariancja i ustalona wcześniej reguła decyzji.
Analityczne: model ma określony kontekst użycia, wersję, założenia i zakres walidacji; korelacja obejmuje masę, środek, tensor i sumy częściowe; odbiorcy używają tego samego wydania.
Konfiguracyjne: stosowalność, status, zatwierdzenie, ślad zmiany, archiwum i lista analiz do regresji są kompletne.

Jeżeli którekolwiek zamknięcie nie zachodzi, właściwym wynikiem nie jest liczba z większą liczbą miejsc po przecinku, lecz jawny brak dowodu. Precyzja zapisu nie może zastąpić wiarygodności.

Synteza

Masa odpowiada na pytanie o bezwładność postępową. Środek masy dodaje pierwszy moment i ramiona sił. Tensor wraz z osiami głównymi opisuje pełną bezwładność obrotową. Każdy poziom wykrywa klasy błędów niewidoczne dla poprzedniego: zgodna masa nie dowodzi zgodnego środka, a zgodny środek nie dowodzi zgodnego tensora.

Studium dwuczłonowe pokazało, że zużycie i odłączenie wymagają ponownego złożenia wokół nowego centrum, a odpowiedź kątowa może zmienić się wielokrotnie. Kampania obojętnego artykułu połączyła reakcje podpór, zmianę orientacji, wahadło, odjęcie przyrządu, liczbową niepewność i pas ochronny w jedną decyzję. Zamiana modułu udowodniła, że zerowa zmiana masy całkowitej może jednocześnie zmienić środek, produkty bezwładności i osie główne oraz uruchomić regresję wielu analiz.

Budżet, model i pomiar nie są konkurencyjnymi opisami. Budżet zarządza dojrzewaniem wiedzy, model zachowuje rozkład i przewiduje stany, a pomiar ujawnia różnice rzeczywistego egzemplarza. Korelacja nie polega na wymuszaniu zgodności, tylko na znalezieniu fizycznej przyczyny, aktualizacji modelu i zachowaniu pozostałej niepewności.

Najważniejsza reguła praktyczna brzmi więc: nigdy nie przekazuje się samej liczby. Przekazuje się kontrolowany rekord właściwości masowych — z konfiguracją, czasem, układem, datum, jednostkami, metodą, niepewnością, stosowalnością i identyfikatorem wydania — oraz dowód, że wszyscy odbiorcy użyli właściwej wersji.

Źródła i zakres stosowalności

Przypisy

NASA Goddard Space Flight Center, GSFC-STD-7000B — General Environmental Verification Standard, 28 kwietnia 2021, https://standards.nasa.gov/sites/default/files/standards/GSFC/B/0/gsfc-std-7000b_signature_cycle_04_28_2021_fixed_links.pdf. Zakres źródła: sprzęt kosmiczny GSFC, ogólne zasady weryfikacji, właściwości masowych i oceny zmian po wyważeniu; nie jest specyfikacją żadnego pocisku ani uniwersalnym źródłem tolerancji.
↩
NASA, NASA-STD-7002B with Change 1 — Payload Test Requirements, zmiana 1 i ponowne zatwierdzenie 24 marca 2023, https://standards.nasa.gov/standard/NASA/NASA-STD-7002. Zakres źródła: ładunki NASA, relacja oceny analitycznej i pomiarowej; stosowana tu wyłącznie jako zasada planowania weryfikacji zależnej od misji.
↩
NASA Johnson Space Center, JSC-65829 Rev. B — Loads and Dynamics Requirements for Spaceflight Hardware, 2024, https://standards.nasa.gov/sites/default/files/standards/JSC/B/0/JSC-65829-Loads-Dyn_Req_for_SF_Hardware-RevB_signed_page2_footer_corrected_0.pdf. Zakres źródła: sprzęt lotów kosmicznych JSC i kontrola modeli masowych dla analiz obciążeń; wymaga dostosowania poza tą domeną.
↩
Steven R. Hirshorn, Linda D. Voss, Linda K. Bromley, NASA, NASA Systems Engineering Handbook, NASA/SP-2016-6105 Rev. 2, 2016, wydanie elektroniczne 2017, https://www.nasa.gov/wp-content/uploads/2018/09/nasa_systems_engineering_handbook_0.pdf. Zakres źródła: ogólna inżynieria systemów, wymagania, techniczne miary postępu, margines i konfiguracja; nie ustanawia granic dla uzbrojenia.
↩
NASA, NASA-STD-5002B — Load Analyses of Spacecraft and Payloads, 10 marca 2026, https://standards.nasa.gov/standard/NASA/NASA-STD-5002. Zakres źródła: analizy obciążeń statków i ładunków kosmicznych; użyta dla zasady spójnego modelu rozkładu masy, nie jako gotowa procedura programu wojskowego.
↩
NASA, NASA-STD-7009B — Standard for Models and Simulations, 5 marca 2024, https://standards.nasa.gov/standard/nasa/nasa-std-7009. Zakres źródła: wiarygodność modeli i symulacji NASA, kontekst użycia, weryfikacja, walidacja i niepewność; zasady przeniesiono po dostosowaniu terminologii.
↩
NASA, NASA-HDBK-7009B — NASA Handbook for Models and Simulations: An Implementation Guide for NASA-STD-7009B, 3 lutego 2026, https://standards.nasa.gov/standard/nasa/nasa-hdbk-7009. Zakres źródła: przewodnik wdrożeniowy do standardu modeli NASA; nie jest źródłem danych konkretnego systemu.
↩
Joint Committee for Guides in Metrology, JCGM 100:2008 — Evaluation of measurement data — Guide to the expression of uncertainty in measurement, wydanie 2008, korekta 2010; JCGM 100:2008/Amd.1:2026 — Amendment 1: Nonlinearity in measurement models, 2026, https://www.bipm.org/en/committees/jc/jcgm/publications. Zakres źródła: ogólne wyrażanie niepewności pomiaru i uzupełnienie dla nieliniowych modeli; dobór modelu i reguły decyzji pozostaje zależny od zastosowania.
↩
Joint Committee for Guides in Metrology, JCGM 101:2008 — Supplement 1 to the GUM: Propagation of distributions using a Monte Carlo method, 2008, https://www.bipm.org/en/committees/jc/jcgm/publications. Zakres źródła: propagacja rozkładów przez model pomiaru; nie określa rozkładów wejściowych dla konkretnego wyrobu.
↩
Joint Committee for Guides in Metrology, JCGM 102:2011 — Supplement 2 to the GUM: Extension to any number of output quantities, 2011, https://www.bipm.org/en/committees/jc/jcgm/publications. Zakres źródła: wielowymiarowe wyniki i kowariancja, stosowane tu do wektora środka i tensora.
↩
Joint Committee for Guides in Metrology, JCGM 106:2012 — Evaluation of measurement data — The role of measurement uncertainty in conformity assessment, październik 2012, https://www.bipm.org/en/doi/10.59161/jcgm106-2012. Zakres źródła: reguły decyzji zgodności, strefy akceptacji i odrzucenia, pasy ochronne oraz ryzyka konsumenta i producenta; nie narzuca tolerancji właściwości masowych.
↩
John K. Ramsey, NASA Glenn Research Center, NASA Aeroelasticity Handbook, Volume 2: Design Guides, Part 2, NASA/TP-2006-212490/VOL2/PT 2, 1 listopada 2006, https://ntrs.nasa.gov/citations/20070008370. Zakres źródła: rozdział o właściwościach masowych modeli, wahadła bi- i trifilarne, kalibracja, odjęcie przyrządu i założenia ruchu; nie jest procedurą odbiorczą dla uzbrojenia.
↩
Jacek Dusza, Paweł Gąsior, Grzegorz Tarapata, Podstawy pomiarów, 2019; Zbigniew Wrzesiński, Niepewność pomiarów w laboratorium balistyki, 2023. Zakres źródła: tor pomiarowy, typ A/B, propagacja, estymacja i dokumentowanie wyniku. Źródła wspierają metodykę, ale nie są specyfikacją stanowiska właściwości masowych ani źródłem danych uzbrojenia; rachunki studiów są własne i syntetyczne.
↩
Chester H. Wolowicz i Roxanah B. Yancey, NASA Flight Research Center, Experimental Determination of Airplane Mass and Inertial Characteristics, NASA-TR-R-433, 1 października 1974, https://ntrs.nasa.gov/citations/19750001990. Zakres źródła: doświadczalne wyznaczanie środka, momentów, produktów bezwładności i nachylenia osi głównej metodami wahadłowymi wraz z wpływem modów pasożytniczych; metody są stosowalne także do innych brył, ale wymagają kwalifikacji konkretnego stanowiska.
↩