Konstrukcja może spełniać klasyczny warunek naprężeniowy i mimo to utracić funkcję wskutek rozwoju pęknięcia. Naprężenie nominalne mniejsze od granicy plastyczności nie rozstrzyga o bezpieczeństwie spoiny, cienkiej ścianki, otworu ani połączenia klejonego. O zachowaniu rzeczywistej części decydują również istniejące nieciągłości, ich orientacja, lokalne pole naprężenia, środowisko, kolejność zdarzeń oraz możliwość wykrycia wady przed osiągnięciem stanu niedopuszczalnego.
Teza monografii brzmi: wiarygodny nadzór odporności na pękanie nie jest pojedynczym obliczeniem odporności materiału, lecz łańcuchem dowodowym łączącym konsekwencję awarii, populację wad, wykrywalność, historię obciążeń i środowiska, model propagacji, wytrzymałość resztkową, konfigurację oraz historię każdej sztuki. Decyzja jest tak wiarygodna, jak najsłabsze z tych ogniw.
Mechanika pękania opisuje relację między wadą, obciążeniem i oporem materiału. Tolerancja uszkodzeń — damage tolerance — odpowiada na pytanie, czy konstrukcja z postulowaną albo wykrytą wadą zachowa wymaganą funkcję do określonej chwili lub wiarygodnej okazji kontroli. Nadzór odporności na pękanie — fracture control — przekształca te ustalenia w program obejmujący klasyfikację części, projekt, materiały, procesy, badania nieniszczące, próby, identyfikowalność, eksploatację i reakcję na zmianę.1,2
Pominięto lokalizacje elementów krytycznych, rzeczywiste widma obciążeń, własności materiałowe, rozmiary wad wykrywalnych, postulowanych i krytycznych, poziomy prób, marginesy, interwały kontroli, limity magazynowania oraz kryteria wycofania dotyczące współczesnych systemów uzbrojenia. Nie przedstawiono procedur demontażu, naprawy ani prób obciążeniowych realnych rakiet, silników, urządzeń inicjujących lub części bojowych. Wszystkie dane liczbowe w studiach I–III są syntetyczne i nie odwzorowują żadnego istniejącego pocisku.
Mechanika pękania, tolerancja uszkodzeń i fracture control
Klasyczna analiza naprężeniowa zakłada zwykle materiał ciągły i pyta, czy naprężenia, odkształcenia oraz przemieszczenia pozostają dopuszczalne. Mechanika pękania wprowadza wadę o określonej geometrii i pyta, czy lokalna siła napędowa jej rozwoju przekracza opór materiału.
Tolerancja uszkodzeń dodaje oś czasu. Stan początkowy wynika z wady postulowanej, wyniku badania nieniszczącego albo modelu populacji wad procesu. Następnie oblicza się wzrost pod wpływem cykli, czasu, temperatury i środowiska oraz sprawdza wytrzymałość resztkową i chronioną funkcję części.
Fracture control obejmuje jeszcze szerszy problem:
- kto ustanawia dane i kryteria;
- jak klasyfikuje się części;
- jak kwalifikuje się materiał, proces i kontrolę;
- co oznacza wynik ujemny NDE;
- kiedy wykonuje się kolejną kontrolę;
- jak wynik przypisuje się do wariantu, partii i konkretnego egzemplarza;
- która zmiana otwiera zatwierdzony dowód ponownie.
NASA-STD-5019A tworzy taki wspólny szkielet dla załogowego sprzętu kosmicznego. Jest cennym źródłem logiki programu, ale nie normą kwalifikacyjną dla rakiet wojskowych.1
Najważniejsze pojęcia należy rozdzielać konsekwentnie:
- wada postulowana jest stanem narzuconym przez wymaganie programu, niezależnie od obserwacji;
- nieciągłość jest cechą materiału albo złącza, która nie musi być niedopuszczalna;
- wskazanie jest odpowiedzią systemu NDE; interpretacja dopiero wiąże je z nieciągłością;
- wymiarowanie — sizing — szacuje geometrię wykrytej cechy wraz z błędem;
- POD — probability of detection jest prawdopodobieństwem wykrycia określonej cechy przez cały zdefiniowany system kontroli;
- wartość dopuszczalna — allowable — ma jawną podstawę statystyczną i dziedzinę zastosowania;
- próba dowodowa — proof test — ma zdefiniowaną funkcję i nie jest automatycznym potwierdzeniem bezwadliwości;
- stan bazowy — baseline — jest zatwierdzoną konfiguracją odniesienia obejmującą geometrię, materiał, proces, model, dane, procedury i ograniczenia;
- nadzór eksploatacyjny — surveillance — sprawdza założenia o starzeniu, populacji i historii;
- granica dowodowa jest stanem, poza którym dostępne dane i modele nie uzasadniają dalszego użycia, nawet jeśli nie przewidują jeszcze fizycznego zniszczenia.
Łańcuch przyczynowy i dowodowy
Argument odporności na pękanie można uporządkować następująco:
funkcja części i konsekwencja jej awarii
│
▼
klasyfikacja krytyczności
│
▼
materiał + proces + geometria wykonana
│
▼
populacja możliwych wad
│
┌─────────┴──────────┐
▼ ▼
NDE, POD i klasyfikacja obciążenia, czas,
morfologii temperatura, środowisko
│ │
└─────────┬──────────┘
▼
warunkowy lub postulowany stan po kontroli
│
▼
lokalne K, J, CTOD, G, FAD oraz da/dN lub da/dt
│
▼
historia propagacji konkretnej sztuki
│
▼
pierwszy utracony warunek funkcjonalny
│
▼
kontrola, resurs, limit kalendarzowy,
ograniczenie zdarzeń, próba lub zmiana projektu
│
▼
konfiguracja + historia sztuki + odstępstwa
│
└────► regresja i ponowne otwarcie dowoduKażda strzałka oznacza założenie wymagające źródła, zakresu ważności, odpowiedzialnej roli i kryterium zamknięcia. Geometria ma reprezentować stan wykonany, nie wyłącznie nominalny model CAD. Widmo powinno pochodzić z pomiaru, zwalidowanej symulacji albo jawnej obwiedni. Dane propagacji muszą odpowiadać materiałowi, orientacji, temperaturze, środowisku i stosunkowi obciążeń. Stan po kontroli powinien wynikać z kwalifikacji kompletnego systemu NDE oraz modelu populacji wad, a nie z nominalnej rozdzielczości aparatury.
NASA-HDBK-5010A podkreśla wielodyscyplinarny charakter fracture control: projekt, analiza, próby, NDE i śledzenie części krytycznych wspólnie ograniczają zagrożenie związane z wadami niewykrytymi.2 Nie oznacza to przeniesienia progów NASA do innej domeny. Program uzbrojenia musi formalnie ustanowić własne wymagania, współczynniki i odpowiedzialności.
Granice transferu między domenami
Trzeba odróżnić fizykę, praktykę inżynierską i wymagania formalne.
| Warstwa | Co można przenieść | Czego nie wolno przenosić automatycznie |
|---|---|---|
| fizyka | równowaga, lokalne pola, $K$, $J$, CTOD, $G$, propagacja, zmiana topologii wady | rozwiązania bez sprawdzenia materiału, skali, geometrii, kontaktu i sposobu obciążenia |
| praktyka inżynierska | identyfikowalność, kwalifikacja NDE, walidacja, kontrola zmian | nazw klas, progów akceptacji i struktur zatwierdzania |
| wymagania domenowe | wymagania formalnie przyjęte lub dostosowane przez właściwy program | współczynników życia, poziomów prób dowodowych, kategorii NASA i interwałów lotniczych |
FAA AC 25.571-1D łączy ocenę tolerancji uszkodzeń, wytrzymałość resztkową, progi i interwały inspekcji z programem utrzymania struktur samolotów transportowych. Jest aktywnym materiałem doradczym FAA, ale nie wymaganiem dla pocisku.10 Polska literatura lotnicza dostarcza wartościowych przykładów związku między widmem, próbą, detalem konstrukcyjnym, badaniem nieniszczącym i współczynnikiem niezawodności. Liczb z programu samolotu nie wolno jednak przenosić do konstrukcji rakietowej bez nowego dowodu.15
Program musi ustanowić co najmniej:
- fazy życia i konfiguracje objęte oceną;
- klasy konsekwencji;
- reprezentatywne obciążenia magazynowania, transportu, gotowości, uruchomienia i lotu;
- kryteria akceptacji niepewności oraz ryzyka technicznego na poziomie systemu;
- wymagany stopień niezależności obliczeń, prób i przeglądu;
- reguły przypisywania wyniku do partii i konkretnej sztuki;
- kryteria kwarantanny, wycofania i ponownego otwarcia kwalifikacji.
Status danych jest częścią wyniku
Liczba bez informacji o pochodzeniu jest niepełnym wejściem.
| Status | Znaczenie | Dopuszczalne wykorzystanie |
|---|---|---|
| wymaganie normatywne | obowiązuje w domenie wskazanego dokumentu | zgodnie z zakresem dokumentu, umową i zatwierdzonym dostosowaniem |
| metoda transferowana | pochodzi z innej domeny, lecz zachowuje użyteczną logikę | po zapisaniu podobieństw, różnic i brakującego dowodu |
| wynik doświadczalny | dotyczy określonej próbki, procesu, środowiska i historii | w granicach wykazanego podobieństwa |
| granica statystyczna | reprezentuje ogon zdefiniowanej populacji | gdy jawne są populacja, metoda, pokrycie i ufność |
| obwiednia deterministyczna | ogranicza zdefiniowany zbiór przypadków | gdy wykazano kompletność zbioru i znak konserwatyzmu |
| wniosek inżynierski | łączy źródła, model, pomiar i osąd | z założeniami, niepewnością i zatwierdzeniem |
| przykład dydaktyczny | wartości dobrano do objaśnienia metody | bez przenoszenia liczb do rzeczywistego projektu |
Status wejścia przechodzi na wynik. Obliczenie wykonane z wartością typową nie staje się wartością projektową przez użycie większej liczby cyfr. Obwiednia obciążeń daje wynik obwiedniowy, nie rozkład prawdopodobieństwa. Największa różnica model–próba w małym zbiorze walidacyjnym nie jest automatycznie statystyczną granicą błędu.
W studiach I–III wszystkie wartości są dydaktyczne. Cyfry robocze umożliwiają odtworzenie rachunku, natomiast decyzje należy zaokrąglać do rozdzielczości uzasadnionej przez dane.
Krytyczność wynika z funkcji
Klasyfikacja krytyczności odpowiada na pytanie o następstwa rozwoju wady. Nie jest klasyfikacją materiału ani prostym podziałem według naprężenia. Ta sama cienka ścianka może być niekrytyczna w jednym układzie i krytyczna w innym, jeżeli zmienia się droga obciążenia, wymagana szczelność, możliwość bezpiecznego zakończenia działania albo sąsiedztwo zespołu magazynującego energię.
| Poziom konsekwencji | Przykładowy skutek | Charakter dowodu |
|---|---|---|
| lokalny | utrata cechy bez zasadniczej utraty funkcji | kontrola jakości, przesiew mechanizmów, zapis odstępstw |
| funkcjonalny | utrata geometrii, sztywności, szczelności lub podsystemu | jawny stan wady, obciążenia, model funkcjonalny, konfiguracja i reakcja |
| krytyczny | pojedyncza wada może spowodować katastrofalną utratę struktury lub niekontrolowane uwolnienie energii | ilościowa wykrywalność, obwiedniowe dane, propagacja, wytrzymałość resztkowa, próby stopniowane, identyfikowalność i niezależny przegląd |
Wzrost konsekwencji rozszerza zakres faz życia objętych analizą i zwiększa rygor danych. NDE wspierające tolerancję uszkodzeń musi wykazać zdolność wykrycia określonej morfologii w rzeczywistej geometrii. Wynik próbki materiałowej nie zastępuje detalu ze spoiną, a próba detalu nie zastępuje podkomponentu, jeśli zmieniają się warunki brzegowe lub droga obciążenia.
Cykl życia i trzy zegary degradacji
Rakieta może długo pozostawać w magazynie, przejść wiele transportów, a następnie wykonać jeden krótki lot. „Jednorazowe użycie” nie oznacza jednego cyklu obciążenia.
| Faza | Wymuszenia | Możliwe mechanizmy |
|---|---|---|
| wytwarzanie i montaż | formowanie, spawanie, utwardzanie, obróbka, wcisk, dokręcanie | pęknięcia skurczowe, brak przetopu, delaminacja, naprężenia resztkowe |
| próby odbiorcze | ciśnienie, drgania, cykle funkcjonalne | wzrost wady, lokalna plastyczność, redystrybucja, rozszczelnienie |
| magazynowanie | temperatura, wilgoć, stałe naprężenie | korozja naprężeniowa, wzrost czasowy, pełzanie kleju, degradacja powłoki |
| transport i przeładunek | drgania szerokopasmowe, udary, zamocowanie | zmęczenie, fretting, utrata zacisku, uszkodzenie udarowe |
| gotowość na platformie | drgania, wiatr, cykle cieplne | pękanie przy otworach, relaksacja połączenia, degradacja interfejsu |
| uruchomienie i start | szybki wzrost obciążenia, ciśnienia i temperatury, udar, akustyka | otwarcie wady, pękanie dynamiczne, utrata drogi obciążenia |
| lot | zginanie, skręcanie, drgania, ciśnienie, gradienty cieplne | zmęczenie termomechaniczne, wzrost szczeliny, odspojenie, utrata stateczności |
Historia egzemplarza powinna zachowywać trzy odrębne zegary:
- liczbę, amplitudę i kolejność cykli;
- czas pozostawania pod obciążeniem;
- czas kalendarzowy w określonym środowisku.
Zmęczenie opisane przez $da/dN$ należy do pierwszego zegara. Korozja naprężeniowa, wzrost wspomagany wodorem albo pełzanie interfejsu mogą zachodzić podczas stałego obciążenia. Starzenie powłoki, chłonięcie wilgoci i zmiana własności kleju mogą zależeć głównie od czasu oraz temperatury. Nie zachodzi więc implikacja:
$$\text{brak cykli} \Rightarrow \text{brak przyrostu uszkodzenia}.$$
Informacja „dziesięć lat magazynowania” jest niepełna bez historii temperatury, wilgotności, naprężenia montażowego, naruszeń powłoki i transportów.
Wada jako zmienna stanu
Nieciągłość materiałowa i matematyczna szczelina nie są synonimami. Por, wtrącenie, rysa powierzchniowa, brak przetopu, podtopienie spoiny, delaminacja i wgniecenie tworzą odmienne pola naprężenia.
Pełny opis stanu wady obejmuje:
- typ i morfologię;
- głębokość, długość i kształt frontu;
- orientację względem naprężeń głównych;
- położenie względem powierzchni, spoiny, otworu, krawędzi i zmiany grubości;
- otwarcie albo kontakt powierzchni;
- materiał i mikrostrukturę przecinane przez front;
- udział modów pękania;
- historię obciążenia i środowiska;
- niepewność detekcji, klasyfikacji i wymiarowania.
Stan można symbolicznie zapisać jako:
$$\mathbf{x}_d= \left[ a,\;c,\;\rho,\;\theta,\;\mathbf r,\; \chi_{\mathrm{kontakt}},\; \psi_I,\psi_{II},\psi_{III},\; \mathcal H \right].$$
Redukcja tego wektora do jednego wymiaru $a$ jest uzasadniona tylko dla kontrolowanej rodziny wad, jeżeli pozostałe zmienne utrzymano w określonym zakresie i wykazano, że nie zmieniają mechanizmu ani decyzji.
Karb nie jest szczeliną
Karb ma skończony promień. W sprężystym modelu jego działanie opisuje współczynnik koncentracji:
$$K_t=\frac{\sigma_{\max}}{\sigma_{\mathrm{nom}}}.$$
Ostra szczelina prowadzi natomiast do osobliwego pola. Maksymalne naprężenie w pojedynczym węźle modelu numerycznego nie jest stabilną wielkością obliczeniową i zwykle rośnie z zagęszczaniem siatki. Mechanika pękania zastępuje je miarą amplitudy pola przy froncie.
Dla idealizowanej szczeliny otwieranej:
$$K_I=Y\,\sigma\sqrt{\pi a}.$$
Współczynnik $Y$ zależy od geometrii, sposobu obciążenia, położenia i kształtu frontu. Inne rozwiązania obowiązują dla szczeliny powierzchniowej, przelotowej, narożnej, przy otworze lub w krótkim ligamencie.
Siła napędowa i opór materiału
Mechanika pękania nie zastępuje analizy konstrukcji. Najpierw trzeba ustalić drogę obciążenia, siły przekrojowe, naprężenia membranowe i zginające, kontakt, podparcie, temperaturę oraz podatność zespołu. Dopiero to pole jest warunkiem brzegowym modelu wady.
Mod I odpowiada otwieraniu. Mod II opisuje ścinanie w płaszczyźnie szczeliny prostopadłe do frontu, a mod III — ścinanie równoległe do frontu. Udział modów może zmieniać się wzdłuż frontu. Pęknięcie może skręcać, rozgałęziać się albo przechodzić między materiałem rodzimym, strefą wpływu ciepła i metalem spoiny.
Współczynnik uwalniania energii $G$ mierzy zmianę energii potencjalnej przypadającą na przyrost powierzchni szczeliny. W jednorodnym materiale liniowo sprężystym, dla modu I:
$$G=\frac{K_I^2}{E'},$$
gdzie:
$$E'=E$$
dla płaskiego stanu naprężenia, a:
$$E'=\frac{E}{1-\nu^2}$$
dla płaskiego stanu odkształcenia.
Osiągnięcie wartości inicjacyjnej nie musi oznaczać natychmiastowej niestabilności. Jeżeli opór materiału rośnie wraz z przyrostem pęknięcia, możliwe jest stabilne rozdzieranie. Stateczność zależy od nachylenia krzywych przyłożonej siły napędowej i oporu, podatności zespołu oraz tego, czy obciążenie jest sterowane siłą, przemieszczeniem czy ciśnieniem.
| Problem | Właściwa rodzina metod |
|---|---|
| jednorodny metal, ostra szczelina, mała strefa plastyczna | liniowo-sprężysta mechanika pękania: $K$, równoważnie $G$ |
| istotna plastyczność i stabilne rozdzieranie | $J$, CTOD, krzywa $J-R$ albo CTOD–$R$ |
| konkurencja pękania i uplastycznienia przekroju | FAD albo analiza sprężysto-plastyczna |
| delaminacja, klej lub interfejs | $G_I/G_{II}/G_{III}$, kryterium mieszanego modu, model kohezyjny i próby |
| wzrost środowiskowy lub pełzanie | $da/dt$ zależne od siły napędowej, temperatury i środowiska |
| szybkie obciążenie przejściowe | analiza dynamiczna i odporność właściwa szybkości obciążenia |
Przybliżony promień strefy plastycznej w płaskim stanie naprężenia wynosi:
$$r_p\approx \frac{1}{2\pi} \left(\frac{K_I}{\sigma_y}\right)^2.$$
Jest to kontrola skali, nie równanie nośności. Podejście $K$ jest uzasadnione, gdy odpowiedź globalna pozostaje zasadniczo liniowa, strefa nieliniowa jest mała względem szczeliny, grubości i ligamentu, a odporność pochodzi z reprezentatywnej próbki.
ASTM E399-24 określa metodę wyznaczania $K_{Ic}$ metali przy zasadniczo liniowo-sprężystej odpowiedzi, silnym skrępowaniu i małej strefie plastycznej.3 Ważność próby materiałowej nie dowodzi automatycznie ważności modelu $K$ dla cienkiej części.
Jeżeli odpowiedź wykazuje istotną nieliniowość i stabilne rozdzieranie, potrzebna jest metoda sprężysto-plastyczna. ASTM E1820-25de2 obejmuje pomiary odporności opisane przez $K$, $J$, CTOD i krzywe oporu w modzie I, przy czym w aktualnym wydaniu wyznaczanie $K_{Ic}$ pozostawiono metodzie E399.4 Trzeba rozdzielić wynik inicjacyjny, ważność wymiarową, krzywą oporu oraz punkt niestabilności konstrukcji.
Diagram oceny uszkodzenia
Diagram FAD zestawia dwa znormalizowane wskaźniki:
$$K_r=\frac{K_I}{K_{\mathrm{mat}}}, \qquad L_r=\frac{\sigma_{\mathrm{ref}}}{\sigma_y}.$$
$K_r$ opisuje wykorzystanie odporności na pękanie, a $L_r$ — zbliżenie do uplastycznienia przekroju. Punkt $(L_r,K_r)$ porównuje się z linią graniczną.
Historyczna linia R6 Revision 3 Option 1 ma postać:
$$K_{r,\lim}(L_r)= (1-0{,}14L_r^2) \left[ 0{,}3+0{,}7\exp(-0{,}65L_r^6) \right],$$
dla:
$$0\leq L_r\leq L_{r,\max}.$$
Przy jednej z powszechnych definicji:
$$L_{r,\max} = \frac{\sigma_y+\sigma_u}{2\sigma_y}.$$
Równanie jest historycznym przykładem procedury R6, a nie deklaracją zgodności z jej bieżącym wydaniem.12 Pomocniczy iloraz:
$$U_{\mathrm{FAD}}= \frac{K_r}{K_{r,\lim}(L_r)}$$
ułatwia porównanie przypadków, lecz kryterium akceptacji musi wynikać z zatwierdzonej procedury.
Naprężenia pierwotne i wtórne nie muszą wchodzić identycznie do obu osi. Naprężenie resztkowe może zwiększać $K_I$, nie będąc jednocześnie obciążeniem powodującym globalne zapadanie przekroju. Sposób superpozycji, relaksacji i redystrybucji pola własnego musi odpowiadać zastosowanej metodzie.
Bramka quasi-statyczna–dynamiczna
Wyniku quasi-statycznego nie wolno automatycznie stosować do startu, udaru ani krótkiego impulsu ciśnienia. Najpierw trzeba ustalić, czy czas narastania obciążenia $\tau_L$ jest długi względem czasów odpowiedzi konstrukcji.
Pomocnicze skale to:
$$\tau_w=\frac{L}{c}, \qquad T_1=\frac{1}{f_1},$$
gdzie $L$ jest charakterystyczną długością drogi propagacji fali, $c$ prędkością fali, a $T_1$ okresem najniższego istotnego modu drgań. Można zbudować wskaźniki:
$$\chi_w=\frac{\tau_w}{\tau_L}, \qquad \chi_m=\frac{T_1}{\tau_L}.$$
Małe wartości obu wskaźników wspierają przybliżenie quasi-statyczne. Wartości porównywalne z jednością wskazują, że fale naprężeń, bezwładność, drgania i faza odpowiedzi mogą wpływać na pole przy szczelinie. Nie istnieje jeden uniwersalny próg $\chi_w$ lub $\chi_m$; granicę należy zwalidować dla geometrii, pasma wymuszenia i wymaganej dokładności.
Druga bramka dotyczy materiału. Trzeba porównać szybkość narastania:
$$\dot K_I=\frac{dK_I}{dt}, \qquad \dot J=\frac{dJ}{dt},$$
oraz szybkość odkształcenia w rejonie frontu z warunkami badania odporności. Szybkość obciążenia może zmienić granicę plastyczności, ciągliwość, mechanizm przełomu i odporność; kierunek tej zmiany zależy od materiału, temperatury, mikrostruktury i środowiska. Statycznego $K_{Ic}$ nie wolno zatem uznawać z definicji za wartość konserwatywną.
ASTM E399-24 zawiera odrębne wymagania dotyczące prób z szybko narastającą siłą i podkreśla, że odporność odnosi się do temperatury, środowiska i szybkości badania.3 Dla odpowiedzi sprężysto-plastycznej potrzebne są odpowiednio zdefiniowane dane dynamiczne $J$ lub CTOD i metoda zachowująca bilans energii. Przy bardzo szybkim zdarzeniu klasyczna całka $J$ wymaga ostrożności z powodu energii kinetycznej i nieustalonego pola.
Bramka dynamiczna ma następującą postać:
historia wymuszenia p(t), F(t), T(t)
↓
pasmo i czas narastania względem L/c oraz T1
↓
czy bezwładność i fale zmieniają lokalne pole?
┌───────┴────────┐
│ │
nie tak
│ │
model quasi-statyczny model przejściowy
│ z bezwładnością, tłumieniem
│ i właściwymi warunkami kontaktu
└───────┬────────┘
↓
K(t), J(t), CTOD(t) i szybkość ich narastania
↓
odporność materiału przy zgodnej szybkości,
temperaturze, środowisku i skrępowaniu
↓
walidacja pomiar–model i decyzjaAnaliza przejściowa powinna zachowywać rzeczywistą historię wymuszenia, masę, podparcie, kontakt i tłumienie. Maksimum obciążenia zewnętrznego nie musi wystąpić jednocześnie z maksimum $K_I(t)$. Zmiana czasu narastania, sztywności mocowania, masy sąsiedniego zespołu, temperatury albo modelu kontaktu otwiera bramkę dynamiczną ponownie.
Propagacja: od widma do decyzji
Dla cyklu liniowo sprężystego:
$$\Delta K=K_{\max}-K_{\min}, \qquad R=\frac{K_{\min}}{K_{\max}}.$$
W środkowym zakresie stabilnego wzrostu często stosuje się prawo Parisa:
$$\frac{da}{dN}=C(\Delta K)^m.$$
Parametry $C$ i $m$ zależą od materiału, orientacji, temperatury, środowiska, stosunku $R$ i procedury badania. Jednostka $C$ zależy od jednostek $a$ oraz $\Delta K$. Prawo Parisa nie opisuje samodzielnie obszaru progowego ani końcowego przyspieszenia.
ASTM E647-24 obejmuje pomiary prędkości wzrostu od rejonu progowego do niestabilności kontrolowanej przez $K_{\max}$. Dokument zwraca uwagę na wpływ naprężeń resztkowych, zamykania szczeliny, stosunku obciążeń, środowiska, historii oraz małych pęknięć.5
Dla dyskretnych bloków:
$$a_{j+1} = a_j+ \sum_k \int_{N_{k,j}}^{N_{k,j+1}} f\!\left[ \Delta K(a),R(a),T,\mathcal E \right]dN + \int_{t_j}^{t_{j+1}} g\!\left[ K(a),T,\mathcal E \right]dt.$$
Nie należy obliczać całego przyrostu z $\Delta K$ wyznaczonym wyłącznie dla początku bloku, jeżeli zmiana $a$ istotnie zmienia $Y$, $\Delta K$ albo prędkość.
Naprężenie resztkowe działa na kilka ekranów
Jeżeli do obciążenia zewnętrznego dochodzi stałe rozciągające naprężenie resztkowe $\sigma_{\mathrm{res}}$, to przy liniowej superpozycji:
$$K_{\max} = Y\left(\sigma_{\max,\mathrm{ext}}+\sigma_{\mathrm{res}}\right) \sqrt{\pi a},$$
$$K_{\min} = Y\left(\sigma_{\min,\mathrm{ext}}+\sigma_{\mathrm{res}}\right) \sqrt{\pi a}.$$
Dla stałego pola własnego:
$$\Delta K = Y\left( \sigma_{\max,\mathrm{ext}}-\sigma_{\min,\mathrm{ext}} \right)\sqrt{\pi a},$$
więc samo $\Delta K$ się nie zmienia. Zmieniają się jednak $K_{\min}$, $K_{\max}$ i:
$$R_{\mathrm{eff}} = \frac{\sigma_{\min,\mathrm{ext}}+\sigma_{\mathrm{res}}} {\sigma_{\max,\mathrm{ext}}+\sigma_{\mathrm{res}}}.$$
Nie wolno dodać pola resztkowego tylko w końcowym ekranie nośności, a jednocześnie stosować danych propagacji właściwych innemu $R$, jeśli nie wykazano niewrażliwości modelu. Podczas postoju należy osobno rozstrzygnąć, czy $\sigma_{\mathrm{res}}$ współtworzy $K_{\mathrm{hold}}$ i czy nie relaksuje się wskutek plastyczności, temperatury lub rozwoju szczeliny.
Sekwencja i krótkie pęknięcia
Kolejność bloków wpływa na wynik przez zmianę rozmiaru przed kolejnym obciążeniem, strefę plastyczną po przeciążeniu, zmianę zamykania frontu i czas postoju przy wysokim $K$. Proste prawo Parisa nie ma pamięci przeciążenia. Zliczanie rainflow zachowuje zakresy cykli, lecz nie zachowuje czasu postoju, temperatury, fazy obciążeń wieloosiowych ani pełnej kolejności.
Próg $\Delta K_{th}$ z badania długiej szczeliny nie jest absolutną granicą braku wzrostu. Małe pęknięcia mogą propagować przy nominalnych wartościach niższych od progu długiej szczeliny.
Diagram Kitagawy–Takahashiego zestawia progowy zakres naprężenia z rozmiarem wady.6 Dla długiej szczeliny:
$$\Delta\sigma_{th}(a)= \frac{\Delta K_{th}} {Y\sqrt{\pi a}}.$$
Długość przecięcia z granicą zmęczeniową detalu $\Delta\sigma_w$ wynosi:
$$a_0= \frac{1}{\pi} \left( \frac{\Delta K_{th}} {Y\Delta\sigma_w} \right)^2.$$
Korekta El Haddada wygładza przejście:
$$\Delta K_{\mathrm{eff}} = Y\Delta\sigma\sqrt{\pi(a+a_0)},$$
$$\Delta\sigma_{th}(a)= \Delta\sigma_w \sqrt{\frac{a_0}{a+a_0}}.$$
Jest to most między granicą zmęczeniową i opisem długiej szczeliny, nie uniwersalne prawo propagacji.7
Proces, POD i stan po kontroli
Trzeba rozdzielić trzy pytania:
- jakie wady może wytworzyć proces;
- z jakim prawdopodobieństwem system NDE wykryje określoną cechę;
- jak opisać stan sztuki po wyniku dodatnim albo ujemnym.
Niech $X_i=(a_i,m_i)$ oznacza stan $i$-tej relewantnej wady w kontrolowanym regionie. Populację przed NDE można opisać rozkładem:
$$p_0(n,x_1,\ldots,x_n),$$
albo deterministycznym zbiorem:
$$\mathcal X_{\mathrm{proc}}.$$
Jeżeli kontrola obejmuje tylko część morfologii:
$$\mathcal X_{\mathrm{proc}} = \mathcal X_{\mathrm{kwal}} \cup \mathcal X_{\mathrm{poza}},$$
gdzie $\mathcal X_{\mathrm{poza}}$ zawiera na przykład wady o niedostępnej orientacji, położone w martwej strefie albo powstałe po kontroli. Kredyt NDE dotyczy wyłącznie zbioru kwalifikowanego.
Co oznacza $a_{90/95}$
Parametr $a_{90/95}$ oznacza rozmiar cechy, dla którego jednostronna dolna 95-procentowa granica ufności krzywej POD osiąga 0,90.8 Nie jest on:
- największą wadą pozostałą po kontroli;
- 95-procentowym kwantylem rozmiaru wady;
- gwarancją wykrycia każdej większej wady;
- błędem wymiarowania;
- opisem częstości występowania wad w procesie.
Krzywa POD należy do całego systemu: części, geometrii, aparatury, procedury, progu, operatora, morfologii i reguły interpretacji. NASA-STD-5009C ustanawia ilościowe wymagania POD dla określonych metalowych elementów NASA wymagających fracture control; nie obejmuje NDE wykonywanego wyłącznie do kontroli procesu ani oceny kompozytów tolerujących uszkodzenia.9
Wynik ujemny dla jednej wady
Dla dokładnie jednej relewantnej wady o stanie $X=(a,m)$:
$$P(\mathrm{ND}\mid a,m)=1-P_D(a,m),$$
a posterior ma postać:
$$p(a,m\mid\mathrm{ND}) = \frac{ [1-P_D(a,m)]\,p_0(a,m) }{ \displaystyle \int\!\!\int [1-P_D(\alpha,\mu)]p_0(\alpha,\mu) \,d\alpha\,d\mu }.$$
Wzór nie jest ogólnym modelem regionu zawierającego nieznaną liczbę wad. Przy $n$ wadach i warunkowo niezależnych odpowiedziach:
$$P(\mathrm{ND}\mid x_1,\ldots,x_n) = \prod_{i=1}^{n} \left[1-P_D(x_i)\right].$$
Jeżeli wspólny stan powierzchni, dostęp, operator lub zakłócenie wpływa na wszystkie odpowiedzi, trzeba wprowadzić zmienną ukrytą $Q$:
$$P(\mathrm{ND}\mid x_1,\ldots,x_n) = \sum_q P(q) \prod_{i=1}^{n} \left[1-P_D(x_i,q)\right].$$
Posterior największej wady otrzymuje się dopiero przez marginalizację pełnego modelu populacji. Samo podstawienie $1-P_D(A)$ jest nieuzasadnione, ponieważ wynik ujemny oznacza przeoczenie wszystkich relewantnych wad.
Jeżeli program ustala limit masy ogona $\varepsilon$, można zdefiniować warunkową granicę:
$$P(A>a_{\mathrm{ND},U}\mid\mathrm{ND}) \leq\varepsilon.$$
Wartość $a_{\mathrm{ND},U}$ zależy równocześnie od populacji procesu i POD. Nie wolno utożsamiać jej z $a_{90/95}$.
Dla wskazania wykrytego:
$$a_{\mathrm{eval},U} = a_{\mathrm{meas}} + e_{\mathrm{sizing},U},$$
o ile $e_{\mathrm{sizing},U}$ jest jednostronną granicą błędu zwalidowaną dla właściwej geometrii, metody i morfologii.
| Wynik | Stan wiedzy | Dalsza czynność |
|---|---|---|
| brak wskazania, wada w domenie POD | warunkowy stan resztkowy | posterior albo zatwierdzona wada postulowana |
| wskazanie wykryte i zwymiarowane | pomiar z niepewnością | niekorzystna granica wymiaru i właściwa morfologia |
| wskazanie wykryte, lecz niewymiarowalne | obecność wady bez wystarczającej geometrii | kwarantanna, dodatkowe badanie albo zatwierdzona obwiednia |
| możliwa wada poza domeną POD | brak wiarygodnego kredytu NDE | metoda komplementarna, próba dowodowa, kontrola procesu albo brak akceptacji |
Powtarzane kontrole i korelacja przeoczeń
Dwie kontrole tej samej części nie tworzą automatycznie dwóch niezależnych dowodów. Trudna morfologia, stan powierzchni, ograniczony dostęp albo trwała cecha egzemplarza mogą korelować kolejne przeoczenia.
Dla ukrytej „inspekcyjności” $Q$:
$$P(\mathrm{ND}_1,\ldots,\mathrm{ND}_k\mid a_1,\ldots,a_k) = \sum_q P(q) \prod_{r=1}^{k} \left[1-P_{D,r}(a_r,q)\right].$$
Odrębny zespół i niezależna ocena mogą ograniczyć niektóre wspólne błędy, lecz nie dowodzą warunkowej niezależności. Model niezależny wolno zastosować dopiero po eksperymencie albo danych terenowych obejmujących powtarzane kontrole tej samej reprezentatywnej populacji.
Pierwszy utracony warunek
Niestabilne pęknięcie nie zawsze jest pierwszym stanem niedopuszczalnym. Wcześniej mogą zostać utracone:
- szczelność;
- sztywność;
- geometria interfejsu;
- stateczność;
- dopuszczalne przemieszczenie;
- ciągłość bariery środowiskowej;
- wymagana częstotliwość własna;
- zdolność mocowania sąsiedniego wyposażenia;
- zapas wobec lokalnego uplastycznienia;
- ważność zwalidowanego modelu albo obwiedni.
Dla ustalonej, jednoparametrowej rodziny wad, w której pozostałe składowe stanu i ścieżka historii są zamrożone, a każdy rozpatrywany warunek pogarsza się monotonicznie wraz ze wzrostem (a), można zapisać:
$$a_{\mathrm{end}} = \min \left[ a_{\mathrm{leak}}, a_{\mathrm{stiff}}, a_{\mathrm{buckling}}, a_{\mathrm{FAD}}, a_{\mathrm{unstable}}, a_{\mathrm{interface}}, a_{\mathrm{evidence}} \right].$$
Równanie nie porządkuje ogólnego wektora stanu wady. Nie wolno go stosować między różnymi morfologiami, orientacjami, historiami obciążenia ani mechanizmami degradacji tylko dlatego, że każdemu z nich przypisano jakiś „rozmiar”. W takim przypadku pierwszy utracony warunek wyznacza się na trajektorii pełnego wektora stanu.
Granica dowodowa nie musi oznaczać fizycznego zniszczenia. Dodatni margines nośności nie unieważnia wcześniejszej utraty funkcji ani utraty ważności argumentu.
Bezpieczny okres użytkowania, bezpieczeństwo po uszkodzeniu i tolerancja uszkodzeń
| Element | Bezpieczny okres użytkowania (safe-life) | Bezpieczeństwo po uszkodzeniu (fail-safe) | Tolerancja uszkodzeń |
|---|---|---|---|
| teza | część zostanie wycofana przed nieakceptowanym uszkodzeniem | po utracie jednej drogi pozostaje wystarczająca funkcja | postulowana wada jest kontrolowana do wykrycia lub limitu |
| główne dane | trwałość, widmo, próby, rozrzut | redystrybucja, nośność resztkowa, wspólne przyczyny | stan początkowy, $da/dN$, $da/dt$, POD i funkcja resztkowa |
| główne zagrożenie | nieoczekiwana wada albo nieznana historia | pozorna redundancja zależna od wspólnego węzła | niewykrywalna morfologia, zły model lub brak dostępu |
| decyzja | limit cykli, zdarzeń lub czasu | bezpieczne zakończenie działania albo kontrola | interwał, limit kalendarzowy, ograniczenie zdarzeń albo wycofanie |
Strategia bezpiecznego okresu użytkowania wymaga reprezentatywnego rozkładu trwałości, widma i kompletnej historii sztuki. Współczynnik niezawodności nie naprawia nieznanego mechanizmu. Strategia bezpieczeństwa po uszkodzeniu wymaga fizycznie niezależnej drogi; dwa elementy oparte na wspólnym okuciu, spoinie lub mechanizmie korozji nie są niezależne. W praktyce potrzebna bywa strategia hybrydowa: dostępna część metalowa podlega okresowemu NDE, a niedostępny węzeł — limitowi kalendarzowemu, nadzorowi środowiska i ograniczeniu liczby transportów.
Próba dowodowa nie zeruje historii
Jeżeli próba dowodowa ma odsiewać wady, trzeba wykazać relację między poziomem obciążenia, rodziną wad, odpowiedzią podczas próby i stanem po jej zakończeniu.
Niech:
- $\mathbf x^-$ oznacza stan przed próbą;
- $L_p(t)$ — historię obciążenia;
- $\mathcal E_p$ — temperaturę, medium i podparcie;
- $\mathcal T_p$ — operator zmiany wady;
- $\mathbf x^+=\mathcal T_p(\mathbf x^-)$ — stan po próbie;
- $\Gamma_p$ — powierzchnię zniszczenia, przecieku albo innego kryterium przerwania.
Zbiór wad odsiewanych:
$$\mathcal S_p= \left\{ \mathbf x^-: \text{trajektoria pod }L_p(t) \text{ osiąga }\Gamma_p \right\}.$$
Zgodnie z definicją zbioru $\mathcal S_p$, przetrwanie wiarygodnie wykonanej i monitorowanej próby wyklucza $\mathcal S_p$, jeżeli historia próby pobudza właściwe miejsce i mod pękania, operator trajektorii oraz powierzchnia $\Gamma_p$ są adekwatne, a osiągnięcie kryterium przerwania zostałoby niezawodnie wykryte. Monotoniczność nie jest potrzebna do samego wykluczenia tak zdefiniowanego zbioru. Jest potrzebna dopiero wtedy, gdy wynik zbiorowy ma zostać zastąpiony prostą skalarną granicą typu „przetrwanie wyklucza każdą wadę większą niż $a_p$” w ustalonej rodzinie wad.
Stan początkowy dalszej analizy jest obwiednią stanów po próbie:
$$\mathcal X_{0,\mathrm{service}} = \left\{ \mathcal T_p(\mathbf x^-): \mathbf x^-\notin\mathcal S_p \right\}.$$
Jeżeli podczas narastania, podtrzymania lub odciążania występuje stabilne rozdzieranie, wzrost środowiskowy albo redystrybucja plastyczna, próba dowodowa może zwiększyć wadę i pozostawić nowe pole naprężeń własnych. Publiczne badania cienkościennych rur aluminiowych i tytanowych pokazują zależność skuteczności takiej próby od materiału, średnicy, grubości, kształtu wady i poziomu obciążenia.13
Niepewność i walidacja
| Rodzaj niepewności | Właściwe traktowanie | Kryterium zamknięcia |
|---|---|---|
| naturalny rozrzut materiału | granica statystyczna | pokrycie materiału, orientacji, grubości, temperatury i środowiska |
| nieznane obciążenie | pomiar, zwalidowana symulacja albo obwiednia | wszystkie zdarzenia należą do zatwierdzonego zbioru |
| kolejność zdarzeń | zachowanie historii lub sekwencje graniczne | żadna wiarygodna kolejność nie zmienia decyzji |
| naprężenia resztkowe | pomiar, walidacja procesu lub obwiednia | rozrzut procesu mieści się w przyjętym polu |
| propagacja | obwiednia danych lub model probabilistyczny | pokryty zakres $\Delta K$, $R$, temperatury i środowiska |
| wykrywalność | POD; osobno wymiarowanie | reprezentatywny system osiąga zatwierdzoną zdolność |
| zależność kontroli | model łączny powtarzanych wyników | dane wykazują strukturę korelacji |
| błąd numeryczny | zbieżność i przypadki odniesienia | dyskretyzacja nie zmienia decyzji |
| błąd modelu | niezależne przypadki walidacyjne | błąd predykcyjny ma granicę właściwą zamierzonemu użyciu |
| brak wiedzy o mechanizmie | konkurencyjne hipotezy | badanie rozstrzygające między hipotezami |
| nieznana historia sztuki | obwiednia albo kwarantanna | historia zostaje odtworzona lub formalnie zastąpiona |
Niepewności mogą być skorelowane. Większa wada może występować w partii o mniej korzystnej mikrostrukturze, większych naprężeniach resztkowych i gorszej wykrywalności. Losowanie każdej wielkości niezależnie może dać zbyt wąski rozkład.
Model powinien przewidywać wynik niewykorzystany do dopasowania. Zgodność z punktami kalibracyjnymi jest kalibracją, nie walidacją. NASA-STD-7009B wiąże wiarygodność modelu z jego zamierzonym użyciem, kryteriami akceptacji, weryfikacją, walidacją, analizą wrażliwości i raportowaniem niepewności.11
Studium I: cienka ścianka i zmiana topologii
Rozpatrywany jest cienkościenny walec o grubości:
$$t=2{,}0 \mathrm{mm}$$
z podłużną półeliptyczną wadą powierzchniową. Przyjęto stały stosunek:
$$\frac{c}{a}=3$$
do chwili penetracji. Chronioną funkcją jest szczelność. Program ćwiczenia postuluje:
$$a_i=0{,}70 \mathrm{mm}.$$
| Wejście | Wartość | Status |
|---|---|---|
| wada $a_i$ | $0{,}70 \mathrm{mm}$ | syntetyczne wymaganie |
| naprężenie pierwotne | $64 \mathrm{MPa}$ | syntetyczna obwiednia deterministyczna |
| naprężenie resztkowe | $36 \mathrm{MPa}$ | syntetyczna obwiednia deterministyczna |
| $\sigma_y$ | $320 \mathrm{MPa}$ | syntetyczna dolna granica deterministyczna |
| $\sigma_u$ | $480 \mathrm{MPa}$ | syntetyczna dolna granica deterministyczna |
| $K_{\mathrm{mat}}$ | $40 \mathrm{MPa}\sqrt{\mathrm m}$ | syntetyczna dolna granica deterministyczna |
Dydaktyczna funkcja geometrii najgłębszego punktu frontu:
$$Y_{\mathrm{surf}}= 1{,}10+0{,}20\frac{a}{t}.$$
Do $K_I$ wprowadzono naprężenie pierwotne i resztkowe:
$$K_I= Y_{\mathrm{surf}} (64+36) \sqrt{\pi a}.$$
Do osi uplastycznienia FAD przyjęto wyłącznie naprężenie pierwotne:
$$L_r=\frac{64}{320}=0{,}200.$$
| $a$ [mm] | $c$ [mm] | $K_I$ [MPa$\sqrt{\mathrm m}$] | $U_{\mathrm{FAD}}$ | Stan |
|---|---|---|---|---|
| 0,70 | 2,10 | 5,49 | 0,14 | wada powierzchniowa |
| 1,20 | 3,60 | 7,49 | 0,19 | wada powierzchniowa |
| 1,80 | 5,40 | 9,63 | 0,24 | wada powierzchniowa |
| 2,00 | 6,00 | 10,30 | 0,26 | penetracja ścianki |
Przy:
$$a=t=2{,}0 \mathrm{mm}$$
front osiąga drugą powierzchnię. Model wady powierzchniowej kończy ważność. Pojawia się szczelina przelotowa, zmieniają się podatność, rozwarcie i możliwość przepływu medium.
Dla niezależnego ekranu po penetracji przyjęto:
$$c=6{,}0 \mathrm{mm}, \qquad Y_{\mathrm{through}}=1{,}10.$$
Wtedy:
$$K_I\approx15{,}1 \mathrm{MPa}\sqrt{\mathrm m}, \qquad \frac{K_I}{K_{\mathrm{mat}}}\approx0{,}38.$$
Jest to wyłącznie ekran składowej pękaniowej. Po zmianie topologii nie zdefiniowano naprężenia odniesienia ani funkcji zapadania przekroju dla szczeliny przelotowej, dlatego nie oblicza się dla niej nowego (L_r) ani wykorzystania FAD. Pełny ekran FAD wymagałby osobnego, zwalidowanego modelu geometrii po penetracji.
Ekran nośności pozostaje poniżej linii granicznej, lecz chroniona funkcja została utracona wcześniej:
$$\text{wzrost wady} \rightarrow \text{penetracja ścianki} \rightarrow \text{utrata szczelności} \quad\text{przed}\quad \text{utratą nośności}.$$
Nie dowodzi to strategii leak-before-break. Potrzebne byłyby dane o rozwarciu, przepływie, medium, progu detekcji, dynamice ciśnienia i czasie reakcji. Studium pokazuje natomiast, że kontynuowanie rozwiązania dla wady powierzchniowej po penetracji jest błędem modelowym, nawet jeśli solver nadal zwraca liczby.
Studium II: od ujemnego NDE do interwału kontroli
Rozpatrywany jest dostępny podczas obsługi metalowy węzeł ze szczeliną krawędziową. Chronioną funkcją jest zachowanie geometrii interfejsu. Studium ogranicza się do jednej dominującej wady o stałej orientacji i położeniu.
Populacja przed kontrolą
Rozmiar $a$ ma rozkład ciągły i jednostajny wewnątrz każdego przedziału:
| Przedział $a$ [mm] | Masa $\pi_j$ |
|---|---|
| 0–0,20 | 0,700 |
| 0,20–0,25 | 0,100 |
| 0,25–0,30 | 0,100 |
| 0,30–0,40 | 0,070 |
| 0,40–0,50 | 0,025 |
| 0,50–0,65 | 0,005 |
Wartości i masy są syntetycznymi danymi dydaktycznymi. Nie reprezentują procesu produkcyjnego rzeczywistego wyrobu.
Skorelowana wykrywalność
Każda sztuka należy do jednej z dwóch trwałych klas inspekcyjności:
$$Q\in\{G,D\},$$
$$P(Q=G)=0{,}80, \qquad P(Q=D)=0{,}20.$$
Klasa $G$ oznacza korzystny dostęp i morfologię, a $D$ — trudny, lecz nadal kwalifikowany przypadek. Krzywe POD mają syntetyczną postać logistyczną:
$$P_D(a\mid q) = \frac{1} {1+\exp[-(a-\mu_q)/s_q]},$$
przy $a$ wyrażonym w milimetrach:
$$\mu_G=0{,}240 \mathrm{mm}, \qquad s_G=0{,}035 \mathrm{mm},$$
$$\mu_D=0{,}330 \mathrm{mm}, \qquad s_D=0{,}050 \mathrm{mm}.$$
Model POD obowiązuje wraz z punktem granicznym:
$$0\leq a\leq0{,}750 \mathrm{mm}.$$
Model propagacji obowiązuje natomiast tylko dla $0\leq a<0{,}750 \mathrm{mm}$. Stan $a\geq0{,}750 \mathrm{mm}$ jest stanem absorbującym „utrata ocenionej obwiedni”. Wartość POD w punkcie $a=0{,}750 \mathrm{mm}$ służy wyłącznie konserwatywnemu warunkowaniu masy tego stanu i nie przedłuża modelu propagacji poza granicę.
Dwuklasowy model $Q\in\{G,D\}$ i jego parametry są autorskim, dydaktycznym modelem zmiennej ukrytej zbudowanym wyłącznie dla tego studium. Nie są wynikiem kwalifikacji konkretnej metody NDE ani modelem populacji rzeczywistego wyrobu.
Po początkowym wyniku ujemnym:
$$P(\mathrm{ND}_0)=0{,}7827386.$$
Kwantyl 99,9% posterioru wynosi:
$$a_{\mathrm{ND0},U}=0{,}393 \mathrm{mm}.$$
Udział klasy trudnej wzrasta do:
$$P(Q=D\mid\mathrm{ND}_0)=0{,}2254.$$
Wynik ujemny zmniejsza prawdopodobieństwo większych wad, lecz jednocześnie zwiększa względny udział egzemplarzy trwale trudniejszych do kontroli.
Obciążenia i dane materiałowe
Dla szczeliny krawędziowej w paśmie o szerokości:
$$W=30 \mathrm{mm}$$
przyjęto:
$$Y(a)= 1{,}12 -0{,}231\frac{a}{W} +10{,}55\left(\frac{a}{W}\right)^2 -21{,}72\left(\frac{a}{W}\right)^3 +30{,}39\left(\frac{a}{W}\right)^4.$$
Jeden blok ma kolejność:
H: 5000 cykli
↓
M: 200 000 cykli
↓
L: 1 000 000 cykli
↓
postój: 365 dniOznaczenia H, M i L są nazwami syntetycznych segmentów amplitudy: odpowiednio wysokiej, średniej i niższej. Litera L nie oznacza „najmniej istotnego” segmentu; ze względu na milion cykli jego wkład może stać się znaczący po wzroście wady.
| Faza | Zakres naprężenia lub obciążenie postoju | Liczba cykli albo czas | Status |
|---|---|---|---|
| H | $\Delta\sigma=170 \mathrm{MPa}$ | 5000 cykli | syntetyczna historia |
| M | $\Delta\sigma=110 \mathrm{MPa}$ | 200 000 cykli | syntetyczna historia |
| L | $\Delta\sigma=70 \mathrm{MPa}$ | 1 000 000 cykli | syntetyczna historia |
| postój | $\sigma_{\mathrm{hold,ext}}=80 \mathrm{MPa}$ | 365 dni | syntetyczna historia |
| naprężenie resztkowe | $\sigma_{\mathrm{res}}=60 \mathrm{MPa}$ | cały blok | syntetyczna obwiednia |
| granica plastyczności | $\sigma_y=320 \mathrm{MPa}$ | — | syntetyczna dolna granica deterministyczna |
| wytrzymałość na rozciąganie | $\sigma_u=480 \mathrm{MPa}$ | — | syntetyczna dolna granica deterministyczna |
| odporność | $K_{\mathrm{mat}}=30 \mathrm{MPa}\sqrt{\mathrm m}$ | — | syntetyczna dolna granica deterministyczna |
Dzięki jawnym wartościom:
$$L_{r,\max} = \frac{320+480}{2\cdot320} = 1{,}25.$$
Dla cykli zewnętrznych:
$$R_{\mathrm{ext}}=0{,}10,$$
$$\sigma_{\max,\mathrm{ext}} = \frac{\Delta\sigma}{1-R_{\mathrm{ext}}}, \qquad \sigma_{\min,\mathrm{ext}} = R_{\mathrm{ext}}\sigma_{\max,\mathrm{ext}}.$$
Pole resztkowe jest w modelu zamrożone: nie relaksuje się i nie redystrybuuje. Wchodzi do $K_{\min}$, $K_{\max}$ i obciążenia postoju, ale znosi się w $\Delta K$. Efektywne stosunki obciążeń wynoszą:
| Faza | $R_{\mathrm{eff}}$ |
|---|---|
| H | 0,317 |
| M | 0,396 |
| L | 0,492 |
Syntetyczny model propagacji
Referencyjną krzywą dla:
$$R_{\mathrm{ref}}=0{,}10$$
interpoluje się liniowo w układzie $\log\Delta K-\log(da/dN)$:
| $\Delta K$ [MPa$\sqrt{\mathrm m}$] | $(da/dN)_{\mathrm{ref}}$ [m/cykl] |
|---|---|
| poniżej 3,0 | 0 w granicach modelu |
| 3,0 | $1{,}0\cdot10^{-12}$ |
| 4,0 | $5{,}0\cdot10^{-11}$ |
| 5,0 | $2{,}0\cdot10^{-10}$ |
| 7,0 | $1{,}0\cdot10^{-9}$ |
| 10,0 | $7{,}0\cdot10^{-9}$ |
| 14,0 | $5{,}0\cdot10^{-8}$ |
| 18,0 | $4{,}0\cdot10^{-7}$ |
W ocenianej domenie nie jest potrzebna ekstrapolacja powyżej ostatniego węzła. Wpływ stosunku obciążeń reprezentuje dydaktyczna korekta:
$$\frac{da}{dN} = \left(\frac{da}{dN}\right)_{\mathrm{ref}} \sqrt{ \frac{1-R_{\mathrm{ref}}} {1-R_{\mathrm{eff}}} }.$$
Podczas postoju:
$$K_{\mathrm{hold}} = Y(a) \left( \sigma_{\mathrm{hold,ext}}+\sigma_{\mathrm{res}} \right) \sqrt{\pi a},$$
$$\frac{da}{dt} = 2{,}0\cdot10^{-8} \left(\frac{K_{\mathrm{hold}}}{4}\right)^2 \mathrm{m/dzień}$$
dla:
$$K_{\mathrm{hold}}>2{,}5 \mathrm{MPa}\sqrt{\mathrm m}.$$
Poniżej progu model przyjmuje zero. Oba prawa są arbitralnymi elementami ćwiczenia, a nie danymi materiałowymi.
Mapa decyzji i jednoznaczny algorytm rachunku
Rachunek prowadzi przez sześć jawnych przekształceń:
| Etap | Wejście | Wynik używany dalej |
|---|---|---|
| 1. populacja | przedziały (a), masy i trwała klasa (Q) | dyskretna aproksymacja rozkładu wspólnego |
| 2. wynik NDE | POD zależna od (a) i (Q) | znormalizowany posterior po wyniku ujemnym |
| 3. propagacja | posterior, historia H–M–L–postój i prawa wzrostu | nowy rozmiar każdego punktu oraz masa w stanie absorbującym |
| 4. granica modelu | przekroczenie (a_{\mathrm{plan}}) | trwały znacznik $\mathsf{ABS}$, bez powrotu do domeny |
| 5. horyzont | kwantyl i limit masy $\varepsilon$ | największa liczba bloków spełniająca kryterium modelowe |
| 6. harmonogram | horyzont oraz wymaganie dwóch okazji | kontrola podstawowa, zapasowa i zakaz rozpoczęcia kolejnego bloku bez wyniku |
Każdy przedział populacji dzieli się na 500 równych komórek. Komórkę reprezentuje punkt środkowy:
$$a_{j,i} = a_{j,\min} + \left(i+\frac12\right) \frac{a_{j,\max}-a_{j,\min}}{500}, \qquad i=0,\ldots,499.$$
Początkowa masa punktu i klasy $q$:
$$w_{j,i,q}^{(0)} = \frac{\pi_j}{500}P(Q=q).$$
Po wyniku ujemnym:
$$\widetilde w_{j,i,q} = w_{j,i,q} \left[1-P_D(a_{j,i}\mid q)\right],$$
$$w_{j,i,q}^{+} = \frac{\widetilde w_{j,i,q}} {\sum_{j,i,q}\widetilde w_{j,i,q}}.$$
Propagację wykonuje się jawną metodą Eulera. Krok cykliczny wynosi najwyżej 100 cykli, a krok postoju dokładnie jeden dzień:
$$a_{n+1}=a_n+ \left.\frac{da}{dN}\right|_{a_n}\Delta N,$$
$$a_{n+1}=a_n+ \left.\frac{da}{dt}\right|_{a_n}\Delta t.$$
Jeżeli w kroku:
$$a_n<a_{\mathrm{plan}}, \qquad a_{n+1}\geq a_{\mathrm{plan}},$$
chwilę przekroczenia zapisuje się przez interpolację liniową:
$$\theta= \frac{a_{\mathrm{plan}}-a_n} {a_{n+1}-a_n}, \qquad 0\leq\theta\leq1.$$
Następnie ustawia się:
$$a=a_{\mathrm{plan}}=0{,}750 \mathrm{mm}$$
i znacznik stanu absorbującego. Taki punkt nie jest dalej propagowany i nie może wrócić do zbioru dopuszczalnego.
Podczas późniejszego, hipotetycznego warunkowania NDE masa absorbująca pozostaje osobnym atomem. Ponieważ rzeczywisty rozmiar poza granicą nie jest już modelowany, jej prawdopodobieństwo wyniku ujemnego oblicza się przy wartości ograniczonej:
$$P(\mathrm{ND}\mid\mathsf{ABS},q) = 1-P_D(a_{\mathrm{plan}}\mid q).$$
Dla monotonicznej krzywej POD jest to niekorzystne założenie względem większych wad: nie przypisuje im mniejszej szansy przeoczenia niż na granicy. Wynik ujemny może zmniejszyć posteriorową masę atomu, ale nie usuwa znacznika $\mathsf{ABS}$ ani nie przywraca ważności dowodu.
Kwantyl 99,9% definiuje się lewostronnie:
$$a_{0{,}999} = \inf\{x:F_a(x)\geq0{,}999\}.$$
Jeżeli masa atomu absorbującego przekracza $10^{-3}$, kwantyl przypada na stan absorbujący.
Pseudokod całego rachunku:
zbuduj punkty środkowe 6 × 500 × 2 i ich masy
pomnóż masy przez 1 − POD(a,q)
znormalizuj
dla każdego bloku:
dla faz H, M, L:
dziel fazę na kroki nie większe niż 100 cykli
dla każdego niezaabsorbowanego punktu:
oblicz Y(a), ΔK(a), Reff i da/dN
wykonaj jawny krok Eulera
jeśli przekroczono aplan:
zapisz interpolowaną chwilę przekroczenia
ustaw a = aplan oraz ABS = prawda
przez 365 kroków jednodniowych:
oblicz Khold(a) i da/dt
wykonaj jawny krok Eulera
obsłuż przekroczenie jak wyżej
zsumuj masę punktów ABS
wyznacz kwantyl 99,9%
przy kolejnej kontroli:
dla punktów żywych użyj 1 − POD(a,q)
dla punktów ABS użyj 1 − POD(aplan,q)
znormalizuj, zachowując znacznik ABSZmniejszenie kroku do 50 cykli i połowy dnia zmienia deterministyczny wynik po czterech blokach o mniej niż $0{,}00006 \mathrm{mm}$ i nie zmienia decyzji. Test ten sprawdza błąd dyskretyzacji, nie błąd modelu propagacji.
Deterministyczny przypadek postulowany
Dla postulowanego stanu, wybranego nieco powyżej kwantyla 99,9% po początkowym wyniku NDE:
$$a_i=0{,}400 \mathrm{mm}$$
otrzymano:
Punkt (0{,}400\ \mathrm{mm}) nie jest obwiednią całej populacji. Rozkład ma niezerową masę do (0{,}650\ \mathrm{mm}), a jego pełny ogon jest zachowany w opisanym dalej rachunku probabilistycznym. Ten pojedynczy przypadek służy jedynie do przejrzystego pokazania kolejności i sprzężenia mechanizmów wzrostu.
| Blok | $a$ na początku [mm] | H [mm] | M [mm] | L [mm] | Postój [mm] | $a$ na końcu [mm] |
|---|---|---|---|---|---|---|
| 1 | 0,400 | 0,0049 | 0,0237 | 0,0000 | 0,0153 | 0,444 |
| 2 | 0,444 | 0,0063 | 0,0338 | 0,0017 | 0,0174 | 0,503 |
| 3 | 0,503 | 0,0090 | 0,0514 | 0,0049 | 0,0204 | 0,589 |
| 4 | 0,589 | 0,0139 | 0,0788 | 0,0194 | 0,0252 | 0,726 |
Granica:
$$a_{\mathrm{plan}}=0{,}750 \mathrm{mm}$$
zostaje osiągnięta w fazie H piątego bloku, po około $4{,}75\cdot10^3$ cykli tej fazy.
Model pokazuje trzy sprzężenia:
- naprężenie resztkowe zwiększa efektywny $R$ i prędkość cykliczną;
- to samo pole zwiększa $K_{\mathrm{hold}}$;
- w miarę wzrostu wady fazy wcześniej progowe zaczynają wnosić mierzalny przyrost.
Granice funkcjonalne i nośności
Syntetyczny model przemieszczenia interfejsu:
$$u_{\mathrm{pred}}(a) = 0{,}0243+0{,}0667a \quad [\mathrm{mm}],$$
przy wymaganiu:
$$u_{\mathrm{int}}\leq0{,}0790 \mathrm{mm}.$$
Nominalna granica funkcjonalna wynosi:
$$a_{\mathrm{interface,nom}}=0{,}820 \mathrm{mm}.$$
Mały syntetyczny zbiór walidacyjny wykazał największą różnicę $0{,}0016 \mathrm{mm}$, a łączny szacunek błędu numerycznego i pomiarowego $0{,}0004 \mathrm{mm}$. Nie są to granice statystyczne. Ćwiczenie nakłada deterministyczną karę:
$$\delta u_{\mathrm{ćw}}=0{,}0020 \mathrm{mm},$$
$$\delta a_{\mathrm{ćw}} = \frac{0{,}0020}{0{,}0667} \approx0{,}030 \mathrm{mm}.$$
Stąd:
$$a_{\mathrm{interface,model}}=0{,}790 \mathrm{mm}.$$
Skorelowana obwiednia rozpatrzonych zmian historii daje:
$$\delta a_{\mathrm{hist}}=0{,}040 \mathrm{mm},$$
a więc:
$$a_{\mathrm{plan}}=0{,}750 \mathrm{mm}.$$
Dla fazy H:
$$\sigma_{\max,\mathrm{ext}} = \frac{170}{1-0{,}10} = 188{,}9 \mathrm{MPa},$$
$$K_I(a) = Y(a) (188{,}9+60) \sqrt{\pi a},$$
$$\sigma_{\mathrm{ref}} = \frac{188{,}9}{1-a/W}.$$
Rozwiązanie:
$$K_r(a)=K_{r,\lim}[L_r(a)]$$
daje:
$$a_{\mathrm{FAD}}\approx2{,}774 \mathrm{mm}.$$
Osobne rozwiązanie:
$$K_I(a)=K_{\mathrm{mat}}$$
daje:
$$a_K\approx3{,}239 \mathrm{mm}.$$
Warunek:
$$L_r=L_{r,\max}=1{,}25$$
zostaje osiągnięty przy:
$$a_{L_r,\max}\approx15{,}833 \mathrm{mm}.$$
Zatem:
$$a_{\mathrm{plan}} < a_{\mathrm{interface,model}} < a_{\mathrm{interface,nom}} < a_{\mathrm{FAD}} < a_K < a_{L_r,\max}.$$
Pierwsza jest granica dowodowa, następnie skorygowana granica funkcjonalna. Pękanie niestabilne i zapadanie przekroju nie sterują decyzją.
Propagacja posterioru
Po początkowym wyniku ujemnym:
| Pełne bloki od $\mathrm{ND}_0$ | Posteriorowa masa $\mathsf{ABS}$ | Kwantyl 99,9% [mm] |
|---|---|---|
| 0 | 0 | 0,393 |
| 1 | $1{,}18\cdot10^{-6}$ | 0,435 |
| 2 | $1{,}08\cdot10^{-5}$ | 0,492 |
| 3 | $2{,}02\cdot10^{-4}$ | 0,572 |
| 4 | $6{,}27\cdot10^{-4}$ | 0,696 |
| 5 | $2{,}48\cdot10^{-3}$ | stan absorbujący |
Program ćwiczenia ustala:
$$\varepsilon=10^{-3}$$
jako największą dopuszczalną posteriorową masę stanów absorbujących:
$$M_{\mathsf{ABS}}(h) = P(\mathsf{ABS}_h\mid\mathcal I).$$
Wielkość $\varepsilon$ nie jest „dopuszczalnym ryzykiem technicznym”. Nie obejmuje konsekwencji awarii, błędu modelu, niepewności danych, wad poza domeną POD ani innych mechanizmów uszkodzenia. Jest wyłącznie limitem warunkowej masy jednego stanu w jednym modelu.
Twardy horyzont modelowy:
$$H_\varepsilon = \max \left\{ h: M_{\mathsf{ABS}}(h)\leq\varepsilon \right\} = 4 \text{bloki}.$$
Dwie okazje kontroli są wymaganiem programu
Podział horyzontu na kontrolę podstawową i zapasową nie wynika z posterioru. Jest odrębnym, syntetycznym wymaganiem organizacyjnym ćwiczenia.
Przyjęto:
- rzeczywista okazja kontroli występuje nie częściej niż co dwa bloki;
- jedna zaplanowana okazja może zostać utracona wskutek niedostępności stanowiska, personelu albo logistyki;
- niewykonana kontrola nie otrzymuje żadnego kredytu;
- wynik i reakcja techniczna muszą zostać zamknięte przed rozpoczęciem następnego bloku;
- prawdopodobieństwo utraty okazji nie jest modelowane ilościowo.
Z tego wymagania organizacyjnego, a nie z równania propagacji, wynika:
- kontrola podstawowa nie później niż po dwóch blokach;
- jeżeli ta okazja zostanie utracona, kontrola zapasowa nie później niż po czterech blokach;
- bez nowego wyniku nie wolno rozpocząć piątego pełnego bloku.
Jeżeli kontrola po drugim bloku zostanie wykonana, kontrola po czwartym bloku jest już kolejną kontrolą planową wynikającą z przyjętego odstępu; nie należy nazywać jej „zapasową” dla wykonanego zdarzenia.
Gałąź A: kontrola po dwóch blokach
Niech $\mathcal H_{0:2}$ oznacza propagację przez pierwsze dwa pełne bloki po wyniku $\mathrm{ND}_0$. Prawdopodobieństwo kolejnego wyniku ujemnego wynosi:
$$P(\mathrm{ND}_1\mid\mathrm{ND}_0,\mathcal H_{0:2})=0{,}877495.$$
Po warunkowaniu:
$$a_{\mathrm{ND1},U}=0{,}370 \mathrm{mm},$$
$$P(Q=D\mid\mathrm{ND}_0,\mathcal H_{0:2},\mathrm{ND}_1)=0{,}2325.$$
Dalsza posteriorowa masa stanów absorbujących:
| Bloki po $\mathrm{ND}_1$ | $M_{\mathsf{ABS}}$ |
|---|---|
| 0 | $2{,}65\cdot10^{-9}$ |
| 1 | $4{,}05\cdot10^{-7}$ |
| 2 | $6{,}14\cdot10^{-6}$ |
| 3 | $9{,}79\cdot10^{-5}$ |
| 4 | $4{,}08\cdot10^{-4}$ |
| 5 | $1{,}04\cdot10^{-3}$ |
Horyzont pozostaje równy czterem pełnym blokom od nowego wyniku. Przy organizacyjnym odstępie dwóch bloków następna kontrola planowa przypada jednak na koniec czwartego bloku całkowitej historii.
Gałąź A: kolejna kontrola na końcu czwartego bloku
Niech $\mathcal H_{2:4}$ oznacza propagację przez bloki trzeci i czwarty po wykonaniu $\mathrm{ND}_1$. Wówczas:
$$P(\mathrm{ND}_2\mid \mathrm{ND}_0,\mathcal H_{0:2},\mathrm{ND}_1,\mathcal H_{2:4}) =0{,}900538.$$
Po warunkowaniu:
$$a_{\mathrm{ND2},U}=0{,}374 \mathrm{mm},$$
$$P(Q=D\mid \mathrm{ND}_0,\mathcal H_{0:2},\mathrm{ND}_1,\mathcal H_{2:4},\mathrm{ND}_2) = 0{,}2358.$$
Kwantyl jest nieco większy niż po NDE$_1$, ponieważ między kontrolami następowała fizyczna propagacja. Wynik ujemny usuwa część masy, ale nie cofa szczeliny.
| Bloki po $\mathrm{ND}_2$ | $M_{\mathsf{ABS}}$ |
|---|---|
| 0 | $1{,}53\cdot10^{-9}$ |
| 1 | $2{,}20\cdot10^{-7}$ |
| 2 | $5{,}15\cdot10^{-6}$ |
| 3 | $4{,}12\cdot10^{-5}$ |
| 4 | $2{,}54\cdot10^{-4}$ |
| 5 | $1{,}20\cdot10^{-3}$ |
Nowy horyzont nadal wynosi:
$$H_{\varepsilon,2}=4 \text{bloki}.$$
Następna kontrola podstawowa przypada po dwóch blokach od NDE$_2$, czyli po szóstym bloku całkowitej historii, a zapasowa nie później niż po ósmym. Bez nowego wyniku nie wolno rozpocząć dziewiątego bloku.
Gałąź B: pierwsza okazja utracona, kontrola zapasowa po czterech blokach
W tej gałęzi nie występuje $\mathrm{ND}_1$. Posterior po $\mathrm{ND}_0$ jest propagowany bez warunkowania przez pełną historię $\mathcal H_{0:4}$. Przed kontrolą masa stanu absorbującego wynosi $6{,}27\cdot10^{-4}$, a kwantyl 99,9% wynosi $0{,}696 \mathrm{mm}$. Obie wartości odpowiadają wierszowi czwartego bloku pierwotnej tabeli.
Prawdopodobieństwo ujemnego wyniku kontroli zapasowej:
$$P(\mathrm{ND}_{B}\mid\mathrm{ND}_0,\mathcal H_{0:4}) =0{,}840898.$$
Po warunkowaniu:
$$a_{\mathrm{NDB},U}=0{,}396 \mathrm{mm},$$
$$P(Q=D\mid\mathrm{ND}_0,\mathcal H_{0:4},\mathrm{ND}_{B}) =0{,}2324,$$
a posteriorowa masa stanu absorbującego spada do $1{,}48\cdot10^{-7}$. Stan absorbujący nie zostaje usunięty; ma jedynie mniejszą masę po wyniku ujemnym. Wszystkie te prawdopodobieństwa są warunkowe względem syntetycznej populacji, krzywych POD, historii obciążeń i reguły atomu absorbującego tego studium; nie są częstościami terenowymi ani oceną ryzyka rzeczywistego wyrobu.
| Bloki po $\mathrm{ND}_{B}$ | $M_{\mathsf{ABS}}$ |
|---|---|
| 0 | $1{,}48\cdot10^{-7}$ |
| 1 | $3{,}34\cdot10^{-6}$ |
| 2 | $3{,}55\cdot10^{-5}$ |
| 3 | $1{,}68\cdot10^{-4}$ |
| 4 | $6{,}73\cdot10^{-4}$ |
| 5 | $2{,}48\cdot10^{-3}$ |
Nowy horyzont wynosi cztery pełne bloki od $\mathrm{ND}_{B}$. Kontrola zapasowa ratuje harmonogram tylko dlatego, że została zakończona przed rozpoczęciem piątego bloku; samo „zaplanowanie” badania nie daje kredytu. Od jej wyniku harmonogram zaczyna nowy odcinek: następna okazja podstawowa przypada po dwóch blokach, a zapasowa nie później niż po czterech.
Model kontrfaktyczny, w którym klasa inspekcyjności jest niezależnie losowana przed każdą kontrolą, daje dłuższy horyzont. Ponieważ niezależność nie została zwalidowana, dodatkowego kredytu nie zaliczono.
Gałęzie decyzji
Wynik ujemny w kwalifikowanej domenie. Stan przechodzi w skorelowany posterior. Nie wolno zastąpić go nominalną czułością aparatury.
Wskazanie wykryte i zwymiarowane. Dla syntetycznego wyniku:
$$a_{\mathrm{meas}}=0{,}46 \mathrm{mm}, \qquad e_{\mathrm{sizing},U}=0{,}03 \mathrm{mm}$$
otrzymuje się:
$$a_{\mathrm{eval},U}=0{,}49 \mathrm{mm}.$$
Obwiedniowa propagacja:
| Stan | $a$ [mm] |
|---|---|
| przy kontroli | 0,490 |
| po 1 bloku | 0,5694 |
| po 2 blokach | 0,6923 |
| podczas 3 bloku | osiągnięcie 0,75 |
Dwa pełne bloki pozostają w ocenionej obwiedni. Przy wymaganiu dwóch okazji kontrolę podstawową planuje się po jednym bloku, a zapasową przed rozpoczęciem trzeciego.
Wskazanie wykryte, ale niewymiarowalne. Nie istnieje liczbowy stan początkowy do propagacji. Decyzją jest kwarantanna techniczna do czasu uzyskania kwalifikowanego wymiarowania albo formalnej obwiedni.
Wynik ujemny poza domeną POD. Jeżeli zmienił się dostęp, stan powierzchni, orientacja, procedura lub aparatura, wynik nie otrzymuje kredytu statystycznego.
Studium pokazuje, że interwał nie jest ilorazem długości życia i współczynnika bezpieczeństwa. Wynika z łącznego modelu procesu, propagacji, wykrywalności, korelacji przeoczeń, granicy funkcjonalnej i rzeczywistych okazji kontroli.
Studium III: pękanie połączenia klejonego
Rozpatrywane jest pojedyncze wewnętrzne odspojenie w kontrolowanej strefie połączenia. Rodzina modelowa zachowuje położenie, kształt frontu, kontakt, odległość od krawędzi, szerokość połączenia i historię kondycjonowania.
Status wszystkich wejść
| Wielkość | Wartość | Status i podstawa niepewności |
|---|---|---|
| zmierzona powierzchnia wskazania | $760 \mathrm{mm^2}$ | syntetyczny wynik pomiaru |
| dodatek powierzchni | $40 \mathrm{mm^2}$ | jednostronna obwiednia deterministyczna ćwiczenia; nie jest kwantylem ani granicą ufności |
| szerokość frontu $b$ | $25 \mathrm{mm}$ | syntetyczna stała rodziny geometrycznej |
| dziedzina $a$ | 8–32 mm | syntetyczna zwalidowana domena modelu |
| $P_I$ | 90 N | syntetyczne obciążenie deterministyczne |
| $P_{II}$ | 200 N | syntetyczne obciążenie deterministyczne |
| $dC_I/da$ | $1{,}11\cdot10^{-3} \mathrm{N^{-1}}$ | syntetyczny wynik zweryfikowanego modelu lokalnego na granicy domeny |
| $dC_{II}/da$ | $3{,}75\cdot10^{-4} \mathrm{N^{-1}}$ | jak wyżej |
| $G_{Ic}$ | $0{,}35 \mathrm{kJ/m^2}$ | syntetyczna obwiednia deterministyczna dolna po kondycjonowaniu; nie kwantyl |
| $G_{IIc}$ | $1{,}10 \mathrm{kJ/m^2}$ | jak wyżej |
| wykładnik $n$ | 1,7 | syntetyczny punktowy współczynnik kalibracji; bez kredytu statystycznego |
„Dolna” oznacza tu minimum narzuconej, zamkniętej macierzy przypadków ćwiczenia, a nie wartość A/B-basis, kwantyl populacji ani dolną granicę ufności. Wynik studium jest warunkowy względem tego zestawu. W rzeczywistym programie potrzebna byłaby statystyczna charakterystyka odporności oraz niepewność współczynnika $n$.
Niekorzystna powierzchnia:
$$A_{d,U}=760+40=800 \mathrm{mm^2}.$$
Modelowa długość:
$$a_U=\frac{A_{d,U}}{b} = \frac{800}{25} = 32 \mathrm{mm}.$$
Jest to górna granica zwalidowanej dziedziny. Zastąpienie powierzchni długością jest dopuszczone wyłącznie dlatego, że rodzina zachowuje kształt, szerokość, położenie i orientację. Odspojenie dochodzące do swobodnej krawędzi byłoby innym przypadkiem.
Energia uwalniania i sprzężenia modalne
Dla wektora obciążenia:
$$\mathbf P= \begin{bmatrix} P_I\\ P_{II} \end{bmatrix}$$
oraz macierzy podatności $\mathbf C(a)$ całkowita energia uwalniania pod obciążeniem sterowanym siłą ma ogólną postać:
$$G_T= \frac{1}{2b} \mathbf P^{T} \frac{d\mathbf C}{da} \mathbf P.$$
Po rozwinięciu:
$$G_T= \frac{1}{2b} \left[ P_I^2\frac{dC_{11}}{da} + 2P_IP_{II}\frac{dC_{12}}{da} + P_{II}^2\frac{dC_{22}}{da} \right].$$
Nie wolno więc sumować dwóch niezależnych pochodnych podatności bez wyjaśnienia wyrazu krzyżowego. W studium przyjęto zwalidowaną partycję modalną, w której współrzędne modalne są ortogonalne energetycznie, a:
$$\frac{dC_{12}}{da}=0$$
w granicach tolerancji numerycznej modelu. Dodatkowo sprawdza się bilans:
$$G_I+G_{II}=G_{\mathrm{total}}$$
wyznaczony niezależnie z całkowitej zmiany energii. Gdyby wyraz krzyżowy nie zanikał albo nie był jawnie rozdzielony między mody, poniższy rachunek byłby nieważny. Udział modów zależy także od przyjętej metody partycji; zmiana metody jest zmianą modelu.
W przyjętej partycji:
$$G_I= \frac{P_I^2}{2b} \frac{dC_I}{da},$$
$$G_{II}= \frac{P_{II}^2}{2b} \frac{dC_{II}}{da}.$$
Na granicy domeny:
$$G_I = \frac{90^2}{2\cdot25} (1{,}11\cdot10^{-3}) \approx0{,}18 \mathrm{N/mm},$$
$$G_{II} = \frac{200^2}{2\cdot25} (3{,}75\cdot10^{-4}) = 0{,}30 \mathrm{N/mm}.$$
Ponieważ:
$$1 \mathrm{N/mm}=1 \mathrm{kJ/m^2},$$
otrzymuje się:
$$G_T=G_I+G_{II}=0{,}48 \mathrm{kJ/m^2},$$
$$\eta= \frac{G_{II}}{G_T} = 0{,}625.$$
Ta sama wartość $G_T$ może prowadzić do innej decyzji przy innym udziale modów.
Odporność mieszanego modu
Kryterium Benzeggagha–Kenane’a ma postać:
$$G_c(\eta) = G_{Ic} + \left(G_{IIc}-G_{Ic}\right)\eta^n.$$
Jest to półempiryczne kryterium, którego parametry muszą pochodzić z danych właściwego systemu materiałowego i sposobu partycji.17
Dla danych ćwiczenia:
$$G_c(0{,}625) = 0{,}35+ (1{,}10-0{,}35)(0{,}625)^{1{,}7} \approx0{,}687 \mathrm{kJ/m^2}.$$
Wskaźnik wykorzystania:
$$g= \frac{G_T}{G_c} = \frac{0{,}48}{0{,}687} \approx0{,}70.$$
Warunkowy wynik brzmi:
$$g<1,$$
zatem model nie przewiduje inicjacji w tej konfiguracji i w granicach syntetycznie ustalonych danych.
ASTM D6671/D6671M-22 rozróżnia inicjację określaną na podstawie odchylenia od liniowości, kryterium 5%/maksimum oraz opcjonalnie pierwszej obserwacji wzrostu. Rozdziela inicjację od wartości propagacyjnych i zwraca uwagę na wpływ mostkowania włókien.16 Podstawowa domena normy obejmuje określoną klasę laminatów jednokierunkowych; wykorzystanie jej logiki do złącza klejonego wymaga wykazania stosowalności.
Gdyby obie składowe siły wzrosły o 25%, energia wzrosłaby z kwadratem obciążenia:
$$G_{T,+25\%} = 0{,}48(1{,}25)^2 = 0{,}75 \mathrm{kJ/m^2},$$
$$g_{+25\%} = \frac{0{,}75}{0{,}687} \approx1{,}09.$$
Taki przypadek przekraczałby kryterium inicjacji. Nie dowodziłby jeszcze niestabilnego zniszczenia całego złącza. O propagacji rozstrzygałyby krzywa oporu, zmiana podatności z $a$, kontakt, sposób sterowania obciążeniem i ewentualna lepkosprężystość.
Granica dziedziny i decyzje
Dla:
$$A_{d,U}=820 \mathrm{mm^2}$$
otrzymano by:
$$a_U=32{,}8 \mathrm{mm}.$$
Równania nadal zwracałyby liczby, lecz przypadek leżałby poza dziedziną $a\leq32 \mathrm{mm}$. Decyzja brzmi:
$$\boxed{\text{dostępny model nie rozstrzyga przypadku}}.$$
Studium rozdziela cztery stany:
- $g<1$ w zwalidowanej dziedzinie — brak przewidywanej inicjacji w modelu;
- $g\geq1$ w zwalidowanej dziedzinie — przewidywana inicjacja, lecz nierozstrzygnięta stateczność propagacji;
- przypadek poza dziedziną — brak rozstrzygnięcia;
- brak reprezentatywnej odporności, partycji modów albo bilansu energii — brak dowodu mimo dostępnego wyniku globalnego.
Zmiana kleju, przygotowania powierzchni, kondycjonowania, grubości spoiny, kształtu frontu, kontaktu, sposobu partycji modów albo współczynnika $n$ otwiera studium ponownie.
Historyczne sprawdzenie łańcucha: obudowa 260‑SL‑1
Eksperymentalna spawana obudowa silnika 260‑SL‑1 o średnicy 260 cali uległa zniszczeniu podczas hydrotestu 11 kwietnia 1965 roku przy ciśnieniu 542 psi, około 56% planowanego ciśnienia próby. Raport NASA TM X‑1194 opublikowano w styczniu 1966 roku. Obudowę wykonano ze stali maraging klasy 250, łączonej głównie automatycznym spawaniem łukiem krytym.14
Źródło pęknięcia znajdowało się w strefie wpływu ciepła podłużnej spoiny. Obszar był wcześniej naprawiany ręcznym spawaniem TIG. Defekt był zanurzony w ściance, zorientowany podłużnie i miał według raportu około 1,4 cala długości oraz około 0,10 cala wymiaru w kierunku grubości. Po awarii znaleziono cztery inne istotne defekty, również pod naprawami TIG.
| Ogniwo | Założenie przed próbą | Co ujawniło dochodzenie |
|---|---|---|
| konsekwencja | krytyczna granica ciśnieniowa | pęknięcie mogło zniszczyć całą obudowę |
| proces | spoiny automatyczne i naprawy ręczne | duże płaskie defekty skupiały się pod naprawami |
| materiał | wysoka wytrzymałość i określona odporność | odporność nie wystarczała dla rzeczywistej populacji wad |
| NDE | zakładano wykrywanie mniejszych nieciągłości | RT i UT miały niższą czułość wobec ciasnych pęknięć |
| sekwencja | kontrola podczas produkcji | między starzeniem a hydrotestem nie wykonano ponownego NDE |
| wada początkowa | większe wady miały zostać odsiane | pozostały większe niewykryte defekty |
| próba dowodowa | miała potwierdzić integralność | obudowa zniszczyła się przed pełnym poziomem próby |
Łańcuch przyczynowy:
naprawy spawalnicze
↓
populacja płaskich defektów pod naprawami
↓
przeceniona zdolność ich wykrywania
↓
brak ponownego NDE po operacji zmieniającej stan materiału
↓
wada większa od tolerowanej przez spoinę
↓
niestabilne pęknięcie podczas narastania ciśnienia próby dowodowejPrzypadek pokazuje, że badanie wykonane przed operacją zmieniającą materiał lub naprężenia nie opisuje automatycznie stanu po tej operacji. Odporność materiału i zdolność NDE muszą być projektowane wspólnie. Proof nie naprawia utraty kontroli procesu. Znalezienie jednego defektu nie kończy analizy, gdy inne podobne wady wskazują na wspólną przyczynę.
Wartością przypadku nie są liczby możliwe do przeniesienia na współczesną obudowę, lecz udokumentowanie kompletnego łańcucha procesu, wady, wykrywalności, odporności, sekwencji kontroli i próby.
Odpowiedzialność organizacyjna
| Rola | Odpowiedzialność | Typowy zapis |
|---|---|---|
| główny konstruktor | funkcja części, droga obciążenia, kryterium końcowe | wymagania i klasyfikacja krytyczności |
| materiały i procesy | własności, środowisko, proces specjalny, partie | wartości projektowe i kwalifikacja procesu |
| obciążenia i dynamika | pełna historia wymuszeń | raport obciążeń i korelacja pomiar–model |
| mechanika pękania | stan początkowy, propagacja i funkcja resztkowa | raport tolerancji uszkodzeń |
| odpowiedzialny poziomu III NDE | procedura, POD, klasyfikacja i wymiarowanie | procedura kwalifikowana, raport POD i dane surowe |
| próby | weryfikacja hipotez i walidacja | plan, konfiguracja, wyniki i korelacja |
| jakość produkcji | wykonanie i identyfikowalność procesu | karta operacji i odstępstwa |
| logistyka i użytkownik | transport, środowisko i zdarzenia | karta sztuki i rejestr przekroczeń |
| zarządzanie konfiguracją | zgodność danych z wariantem, partią i zmianą | stan bazowy i wykaz zastosowalności |
| niezależny organ techniczny | ocena kompletności i wspólnych błędów | niezależny przegląd |
| uprawniony zatwierdzający | akceptacja ograniczeń i ryzyka systemowego | decyzja, limit, interwał albo zakaz użycia |
Braku dokumentu nie należy maskować współczynnikiem. Brak reprezentatywnej POD wymaga kwalifikacji systemu NDE. Nowa morfologia wymaga nowego modelu lub próby. Zaginiona historia sztuki wymaga zatwierdzonej obwiedni albo kwarantanny technicznej.
Regresja dowodu po zmianie
Analiza wrażliwości bada parametry aktualnej konfiguracji. Regresja techniczna ustala, które części zatwierdzonego argumentu pozostają ważne po zmianie.
| Zmiana | Minimalny zakres regresji |
|---|---|
| wymiar lub tolerancja | pole globalne, $K/J/G$, FAD, propagacja i dostęp NDE |
| materiał lub obróbka cieplna | wytrzymałość, odporność, wzrost, krótkie pęknięcia i środowisko |
| dostawca lub proces | populacja wad, próbki technologiczne, NDE i identyfikowalność |
| procedura spawania | mikrostruktura, geometria, naprężenia resztkowe i dane złącza |
| klej lub przygotowanie powierzchni | odporność interfejsu, udział modów, środowisko i wykrywalność |
| widmo lub dodatkowy transport | propagacja, sekwencja, stan sztuki i limity |
| czas narastania lub udar | bramka dynamiczna, bezwładność, $\dot K/\dot J$ i odporność dynamiczna |
| metoda albo próg NDE | POD, fałszywe wskazania, stan po kontroli i decyzja |
| naprawa | nowa geometria, materiał, pole resztkowe, model i ponowne NDE |
| próba dowodowa | stan przed i po próbie, wzrost, pole resztkowe i życie pozostałe |
| wersja solvera lub skryptu | przypadki odniesienia, zbieżność, jednostki i kontrola różnic |
| kryterium funkcjonalne | ponowne przejście do pierwszego stanu niedopuszczalnego |
| klasa konsekwencji | pełna bramka danych, prób, NDE i niezależnego przeglądu |
Pełne ponowne otwarcie jest wymagane po zmianie mechanizmu uszkodzenia, architektury drogi obciążenia, klasy konsekwencji, zasadniczego rachunku, kryterium funkcjonalnego, zdolności wykrywania, modelu korelacji powtarzanych kontroli, dziedziny walidacji, partycji modów albo tezy o niezależności dróg bezpieczeństwa po uszkodzeniu.
Regresja powinna odpowiedzieć:
- czy zmieniło się ogniwo kontrolujące;
- czy zachował się znak konserwatyzmu;
- czy wynik pozostaje w sprawdzonej dziedzinie;
- czy pojawił się nowy mechanizm;
- czy stan konkretnej sztuki odpowiada modelowanemu;
- czy wszystkie wejścia należą do tej samej konfiguracji;
- czy zmiana narusza reprezentatywność próby albo POD;
- czy zmieniła się zależność między kontrolami;
- czy zmienił się pierwszy utracony warunek;
- czy zmienia się interwał, limit kalendarzowy albo liczba zdarzeń.
Sumy kontrolne wejść, wersjonowanie modeli i przypadki odniesienia dowodzą tożsamości artefaktów cyfrowych. Nie dowodzą prawdziwości modelu fizycznego.
Reakcja na wykryte uszkodzenie
Dla potwierdzonego, odłączonego i nieenergetycznego elementu konstrukcyjnego należy przerwać użytkowanie, zabezpieczyć konfigurację, oznaczyć położenie wskazania, zachować dane surowe i przekazać część właściwej organizacji technicznej. Powierzchni przełomu nie należy dopasowywać, agresywnie czyścić ani dotykać przed dokumentacją fraktograficzną.
Jeżeli element może zawierać ciśnienie, energię sprężystą albo ruchome podzespoły, stan zerowej energii może ustanowić wyłącznie uprawniony personel zgodnie z dokumentacją urządzenia. Nie wykonuje się improwizowanego prostowania, nacinania wady ani próby ciśnieniowej.
Uszkodzona lub nieznana rakieta, silnik, inicjator, część bojowa albo zespół pirotechniczny pozostaje w kompetencji uprawnionych służb EOD i personelu danego systemu. Nie należy zbliżać się do znaleziska, dotykać go, przemieszczać, rozkręcać, wiercić, szlifować ani pobudzać w celu „sprawdzenia” pęknięcia. Należy odsunąć ludzi od zagrożenia, powiadomić numer alarmowy 112 i wykonywać polecenia służb.
Warunki zamknięcia dowodu
Argument jest zamknięty dopiero wtedy, gdy niezależny zespół potrafi odtworzyć:
- konfigurację i historię ocenianej sztuki;
- klasyfikację krytyczności;
- mechanizm prowadzący od wady do pierwszego utraconego warunku;
- populację i morfologię wad przed kontrolą;
- domenę POD oraz niepewność wymiarowania;
- warunkowanie stanu po wyniku ujemnym, w tym traktowanie masy absorbującej;
- zależność między powtarzanymi kontrolami;
- źródło, amplitudę i kolejność obciążeń całego cyklu życia;
- trzy zegary degradacji;
- lokalne $K_{\max}$, $K_{\min}$, $\Delta K$, $R$, $J$ lub $G$;
- bramkę między odpowiedzią quasi-statyczną i dynamiczną;
- chwilę każdej zmiany topologii albo mechanizmu;
- relację między próbą dowodową, stanem przed próbą i stanem po próbie;
- pochodzenie każdej rezerwy i jej status statystyczny;
- regułę wyprowadzenia interwału oraz odrębne wymagania organizacyjne;
- limit posteriorowej masy stanów absorbujących, bez nazywania go ryzykiem systemowym;
- decyzję po wyniku ujemnym, dodatnim i niewymiarowalnym;
- wszystkie zmiany otwierające argument ponownie.
Mechanika pękania dostarcza języka siły napędowej i oporu. Tolerancja uszkodzeń przekształca go w wykrywalność, wzrost, czas i funkcję resztkową. Nadzór odporności na pękanie utrzymuje zgodność argumentu z rzeczywistą konfiguracją i historią sztuki. Decyzja powstaje dopiero wtedy, gdy droga od konsekwencji awarii przez proces, stan wady, obciążenie i propagację prowadzi do jawnej granicy, odtwarzalnego terminu kontroli oraz zdefiniowanej reakcji na każdy wynik.