Link: OLE-Object-Data Duktilní deformace, část 4 Deformace eliptické nebo elipsoidální částice je popsána vztahem: kde A je matice elipsy (či elipsoidu) před deformací A‘ je matice elipsy (či elipsoidu) a D je matice deformace. Pouze z tvaru a orientace deformovaných eliptických částic nelze určit velikost deformace! Pro takové určení je nutné vyslovit doplňující předpoklad, který blíže specifikuje celkovou stavbu celého souboru částic. Deformační analýza se pak soustředí nikoli na změnu parametrů popisujících geometrii jedné částice, ale na změnu parametrů popisujících celkovou stavbu. Deformace eliptických částic – všesměrná stavba Jedním z nejjednodušších předpokladů je předpoklad původně všesměrné stavby. Předpokládáme tedy, že částice měly původně eliptický tvar, jejich dlouhé osy byly ale orientovány chaoticky do všech směrů (neměly přednostní orientaci). Budeme-li uvažovat jednotnou původní elipticitu R[i], pak se nám tato stavba znázorní v R[f]/f grafu jako vodorovná přímka (tvořena body popisujícími stav jednotlivých objektů). Body popisující stav objektů po deformaci pak bude v R[f]/f grafu vymezovat křivku - při dostatečně velké deformaci jde o uzavřenou křivku cibulovitého tvaru a nazývá se proto „cibulová“ křivka. Tvar cibulové křivky (její „protažení“ ve směru paralelním s hlavním směrem deformace) závisí na velikosti (elipticitě) deformace. Šířka cibulové křivky je funkcí původní elipticity částic R[i] a elipticity deformace R[s] a nazývá se fluktuace F: Opět si všimněme vztahu pro R[f] a podívejme se, jak bude vztah zjednodušen pro elipsy, jejichž dlouhé osy byly původně paralelní se směrem maximálního protažení, nebo byly naopak na tento směr kolmé (q=0°, q=90°). Předpokládáme koaxiální deformaci, pak také dlouhé osy těchto částic po deformaci jsou buď paralelní se směrem maximálního protažení, nebo jsou na tento směr kolmé (f=0°, f=90°): q=0°: q=90°: Lze ukázat, že: q=0° pak f=0° q=90° a R[s]R[i]: ... f=0° Pak tedy pro q=0° nabývá vztah pro R[f] tvaru: pro q=90° je nutno vztah upravit (vzhledem k aplikaci funkce tangens), lze ale ukázat, že platí: Je-li R[s]R[i]: Předpokládejme jednotnou původní elipticitu. Pak částice s původní dlouhou osou paralelní se směrem maximálního protažení bude mít po deformaci největší elipticitu (R[f-maximum]), částice s původní dlouhou osou kolmou na směr maximálního protažení bude mít po deformaci nejmenší elipticitu (R[f-minimum]). Přitom jsme před chvílí určili jednoduché vztahy pro velikosti těchto konečných elipticit částic: tedy: Těchto vztahů využívá tzv. metoda extrémních tvarů. Podmínky platnosti této metody extrémních tvarů jsou: - Soubor analyzovaných částic musí zahrnovat jak částici, jejíž dlouhá osa byla přibližně paralelní se směrem maximálního protažení, tak i částici, jejíž dlouhá osa byla na tento směr kolmá (tj. nejlépe původně všesměrná orientace částic) - Původní elipticita částic s extrémními tvary, které jsou použity k analýze, byla vzájemně podobná; soubor neobsahuje částici s původně abnormálně eliptickým tvarem, případně je taková částice známa a proto vyloučena z dalších analýz (tj. nejlépe původně tvarově uniformní částice) Průměrná matice Pouze z tvaru a orientace deformovaných eliptických částic nelze určit velikost deformace! Pro takové určení je nutné vyslovit doplňující předpoklad, který blíže specifikuje celkovou stavbu celého souboru částic. Můžeme-li pak tento soubor částic jednoduše matematicky popsat (nějakou rovnicí), lze pak deformaci tohoto souboru částic vyjádřit také jako transformaci daného matematického popisu. Toshihiko Shimamoto a Yukio Ikeda v roce 1976 navrhli popsat soubor eliptických či elipsoidálních částic jejich průměrnou maticí, která representuje tzv. průměrnou elipsu či průměrný elipsoid (averaged ellipse, averaged ellipsoid). Průměrná matice je definována jako matice, jejíž členy jsou aritmetickým průměrem příslušných členů všech matic popisujících jednotlivé objekty. A ... průměrná matice; a[k] ... matice popisující k-tý objekt Dlouhá osa průměrné elipsy odpovídají směru přednostní orientace jednotlivých eliptických částic. Elipticita průměrné elipsy je úměrná jednak elipticitě průměrných částic a jednak míře jejich uspořádání. Při paralelním uspořádání všech dlouhých os jednotlivých částic je elipticita průměrné matice maximální - nemůže být ale větší, než elipticita nejvíce eliptické částice. Při všesměrném uspořádání jednotlivých částic je elipticita průměrné matice minimální - tvar průměrné elipsy odpovídá nebo se alespoň blíží kružnici. Všimněme si opět transformační rovnici popisující deformaci eliptických či elipsoidálních částic: Průměrnou matici A‘[t] popisující celý soubor deformovaných objektů sestavíme jako aritmetický průměr jednotlivých matic deformovaných objektů A‘[k]. Aritmetický průměr matic jednotlivých objektů před deformací je ovšem průměrnou maticí A[t] popisující celý soubor objektů před deformací. Pak tedy platí jednoduchá transformace: kde A‘[t] je tzv. průměrná konečná matice, A[t] je tzv. průměrná počáteční matice a D je matice deformace. Pokud je původní stavba všesměrná, pak její průměrná matice representuje kruh (2D) nebo povrch koule (3D) a lze jde tedy o C-násobek jednotkové matice: Pokud byla původní stavba všesměrná (nebo také „náhodná stavba“ - random fabric), pak průměrná konečná matice této stavby má tvar a směr deformační elipsy (či elipsoidu)! Ovšem pozor: Pokud původní stavba nebyla všesměrná, pak průměrná konečná matice této stavby nemá ani tvar ani směr shodný s tvarem a směrem deformační elipsy (či elipsoidu)! Testy symetrie Soubory částic, které měly před deformací všesměrnou stavbu, mají tendenci vykazovat v R[f]/f grafu symetrické uspořádání kolem osy průměrné elpticity a osy průměrné orientace dlouhé osy. Naopak soubory částic, které měly před deformací výrazně přednostně uspořádanou stavbu, mají někdy (ne vždy!!!) tendenci vykazovat v R[f]/f grafu asymetrické uspořádání kolem osy průměrné elpticity a osy průměrné orientace dlouhé osy. Asymetrické uspořádání bodů reprezentujících tvar a orientaci deformovaných objektů v R[f]/f grafu tedy vylučuje možnost aplikaci předpokladu původně všesměrné stavby. Symetrické uspořádání bodů reprezentujících tvar a orientaci deformovaných objektů v R[f]/f grafu tedy nevylučuje možnost aplikaci předpokladu původně všesměrné stavby. Nedokazuje však, že stavba původně všesměrná byla - pouze nevylučuje, že stavba původně všesměrná mohla být! Přesto je velmi užitečné provést test symetrie uspořádání bodů v R[f]/f grafu, abychom se vyhnuli nesprávné aplikaci předpokladu původní všesměrné stavby v případech, kdy tento předpoklad je zcela jistě nesprávný a vede k chybným výsledkům. Jednoduchým testem symetrie je tzv. I[sym] test. Test je založený na odvození parametru symetrie I[sym]. R[f]/f graf je rozdělen osou harmonického průměru elipticit částic a osou vektorového průměru směrů dlouhých os částic na čtyři pole (A, B, C a D). Parametr I[sym] je pak definován vztahem: Parametr symetrie I[sym]. nabývá hodnot od 0 do 1, čím vyšší je jeho hodnota, tím více je distribuce bodů v R[f]/f grafu symetrická. Pokud je hodnota parametru I[sym] větší než kritická hodnota, pak je distribuce bodů v R[f]/f grafu symetrická. Pokud je hodnota parametru I[sym] menší než kritická hodnota, pak je distribuce bodů v R[f]/f grafu symetrická. Kritické hodnoty parametru Isym závisí na velikosti deformace (na elipticitě deformace Rs) a na počtu objektů (n) tvořících soubor. počet objektů: 35 60 100 200 20 Rs: 1.5 0.3 0.51 0.60 0.74 0.82 (0.4) (0.63) (0.67) (0.78) (0.85) 2.0 0.5 0.63 0.73 0.80 0.86 (0.5) (0.63) (0.77) (0.82) (0.88) 3.0 0.5 0.63 0.73 0.80 0.87 (0.6) (0.63) (0.77) (0.82) (0.88) 5.0 0.5 0.63 0.73 0.82 0.87 (0.6) (0.63) (0.77) (0.82) (0.88) 10.0 0.6 0.63 0.73 0.82 0.87 (0.6) (0.63) (0.77) (0.84) (0.89) kritické hodnoty Isym pro 95% (90%) pravděpodobnost asymetrie. Důležitým testem symetrie je také tzv. test q-distribuce. Je založený na sledování hodnoty c^2 a jejím porovnání s kritickou hodnotou, která plyne z předpokladu náhodné distribuce v původní orientaci dlouhých os částic. O je počet částic pozorovaný v dané části grafu E je počet částic předpokládaný ve stejné části grafu na základě modelu Rf/f graf je rozdělen na části, které mají v případě splnění podmínky o původně náhodné stavbě obsahovat při správně určené velikosti deformace (ohraničení částí závisí na Rs) shodné počty prvků. Hodnota c^2 by tak měla být nízká. Je-li hodnota pro některé elipticity deformace Rs podkritická, lze tuto elipticitu deformace pokládat za řešení deformační analýzy.