Ghostscript wrapper for D:\BBB-ARCH\MTS68\MTS68_t6z1_4_PDF\mts68_t6z1.pdf M E C H AN I K A TEORETYCZNA I STOSOWANA 1,6 (1968) O E F E KC I E SKALI CIAŁA KRU CH EG O WYTRZYM U JĄ CEGO U STALON Ą KON C EN TRAC JĘ M IKRO- U SZ KOD Z EŃ JAN U SZ M U RZEWSKI, JÓZEF SOJKA (KRAKÓW) 1. Model oś rodka i mechanizm zniszczenia Jedn ostkowa gran ica wytrzymał oś ci oś rodka kruchego zależy od absolutnych rozmia- rów ciał a, w tym sensie, że jej wartość oczekiwana jest maleją cą funkcją obję toś ci ciał a. P rawidł owość ta, zwan a efektem skali, został a zbadan a doś wiadczalnie i wytł umaczona n a gruncie rach u n ku prawdopodobień stwa przez W. WE I BU LI A [12], J. I. F REN KIELA [5], B. B. C Z E C Z U LI N A [2] i in n ych autorów [1, 3]. P odstawą heurystyczną teorii efektu skali jest t ak zwan a h ipoteza «najsł abszego ogniwa w ł ań cuchu», wedł ug której noś ność graniczną cał ego ciał a okreś la się wytrzymał oś cią miejscową jego najsł abszego elementu oraz hipoteza, że w m ateriale znajdują się przypadkowe, bezł adnie rozmieszczone uszko- dzen ia lokaln e, czyli miejsca o zaniż onej wytrzymał oś ci. M etody rozważ ań, jakie są przedstawicn e n a podstawie tych hipotez, moż na by z grub- sza rzecz biorą c podzielić n a dwa kierun ki. D o kierun ku pierwszego zaliczymy te prace, w których wprowadza się cią głe funkcje rozkł adu minimalnych wartoś ci lokalnej wytrzyma- ł oś ci m ateriał u. W. WEIBTJLL ja ko pierwszy t a k wł aś nie postą pił [12], a zapropon owan e przez niego prawo rozkł adu wartoś ci m inim alnych okazał o się trafne i zyskał o duże roz- powszechnienie [13]. D alszy rozwój tego kierun ku, dla którego charakterystyczne jest traktowan ie m ateriał u ja ko oś rodka cią gł ego, polegać chyba bę dzie na zastosowaniach teorii stacjon arn ych funkcji stochastycznych [8]. D rugi kierun ek — to uż ycie klasycznych schematów losowania i «dyskretnych» funkcji rozkł adu, n p . funkcji ro zkł ad u Bernoulliego. Stosują c tę m etodę należy sobie wyobrazić ciał o ja ko zbiór czą stek, w których mogą wystę pować defekty. Z reguł y uwzglę dnia się n iejedn orodn ość tych defektów i charakteryzuje sieje funkcją prawdopodobień stw G aussa, P earson a lub inną . Ten p u n kt widzenia został przedstawiony n p . w pracy T. A. K O N T O - ROWEJ i J. I. F R EN KI ELA [5]. W tym uję ciu poszukiwan ia teoretyczne zwrócone są w kierun ku okreś len ia prawdopodobień stwa rzadko wystę pują cych, najbardziej osł abionych czą stek m ateriał u, bo te determinują wytrzymał ość cał ego zbioru. Jednakże w znanych n am pracach nie spotkaliś my konsekwentnie zastosowanego schematu rzadkich zdarzeń i prawa prawdopodobień stw P oisson a. N awet S. D . WO Ł K Ó W [14], który przedstawia w gruncie rzeczy P oissonowski schem at losowania, dokonuje ostatecznie bł ę dnego przejś cia gra- nicznego i aproksym uje koń cowe wyniki G aussowskim rozkł adem asymptotycznym (str. 93 m on ografii [14]). 44 JAN U SZ M U RZEWSKI, JÓZEF SOJKA W niniejszej pracy zjawisko kruchego pę knię cia analizować bę dziemy w oparciu o uogól- nioną (tak jak u S. D . WOŁ KOWA) hipotezę najsł abszego ogniwa w ł ań cuchu i hipotezę przypadkowego rozł oż enia defektów przy uwzglę dnieniu dyskretnej struktury m ateriał u. Z arodek m akropę knię cia w oś rodku traktować bę dziemy ja ko zdarzen ie rzadkie, wywoł ane nadzwyczajną lokalną kumulacją mikrorys, i prowadzą ce ostatecznie do zastosowan ia prawa prawdopodobień stw P oissona. Wyniki, które przedstawimy, bę dą ogólniejsze od klasycznych wzorów opisują cych efekt skali i dlatego ich zakres stosowalnoś ci w praktyce może być szerszy. Przy analizie efektu skali ograniczymy się do oś rodków quasi- jednorodnych, sprę - ż ysto- kruchych, począ tkowo doskonale izotropowych. P rzez quasi- jednorodność rozu- miemy jedn orodn ość w skali makroskopowej przy jednoczesnej niejednorodnoś ci w skali mikroskopowej. Sprę ż ysto- krucha wł asność m ikroelem entu polega n a tym , że p o od- kształ ceniu sprę ż ystym i osią gnię ciu przez m ikronaprę ż enie s odpowiedniej granicy P nastę puje momentalnie opadnię cie mikro- naprę ż enia do zera. Sprę ż ysto- krucha wł asność makroelementu polega n a tym, że aż do granicy sprę ż ystoś ci Q m ateriał podlega prawu H ooke'a, nastę pnie mogą zachodzić nieliniowe odkształ cenia n a skutek spę kan ia, czyli tzw. plastycznoś ci destrukcyjnej i dopiero po osią gnię ciu przez naprę ż enie er granicy wy- trzymał oś ci R nastę puje rozerwanie ciał a. Rozpatrywany jest stan naprę ż enia doskon ale jednoosiowy i quasi- jednorodny, a wię c s jest zmienną losową a a — wielkoś cią w peł n i okreś loną. Przy zał oż eniu doskonał ej izotropii począ tkowej m ateriał u m ikrorysy wystę pować bę dą w pł aszczyznach prostopadł ych do kierun ku dział an ia n aprę ż eń rozcią gają cych; N a skutek takiego spę kania m ateriał nabywa wł asnoś ci an izotropowych. Z agadn ien ie nabytej anizotropii «destrukcyjnej» w przestrzennym stanie n aprę ż en ia był o już przedm io- tem badań doś wiadczalnych i teoretycznych. Przy uproszczon ym podejś ciu do tego zagad- nienia, przedstawionym w pracy [9] i ograniczeniu się do jednoosiowego stan u n aprę ż en ia jedynym efektem tej anizotropii bę dą wyż ej wspom niane odstę pstwa od prawa H ooke'a.. Z akł adamy, że m ateriał ma m ikrostrukturę podobn ą do tej, jaką opisywali w swoich pracach J. M U R Z EWSKI [9] i S. D . WO Ł K Ó W [14], tzn . polegają cą n a podziale oś rodka n a elementy obję toś ciowe m akroskopowe ii (I rzę du) i m ikroskopowe O (I I rzę du), czyli n a tzw. punkty fizyczne pierwszego i drugiego rodzaju, w stosun ku do których stosuje się róż ne prawa fizyczne. Wymiary elementów m akroskopowych szacowane są [9] n a 0,1 m m , czyli 10~2 cm, tak ż eby podlegać badan iom laboratoryjnym , elementów zaś m ikroskopo- wyc h — n a 103 A, czyli 10~5 cm, t ak ż eby zachować jeszcze wł asnoś ci m aterii cią gł ej. Przy tak pomyś lanym geometrycznym podziale oś rodka ilość elementów O w elemencie Q wynosi ilość zaś pun któw pierwszego rodzaju jest tego samego rzę du ju ż przy stosun kowo m ał ych rozm iarach elementów konstrukcyjnych, bo ~ 10 cm, (1.2) N- EF EKT SKALI PRZY KONCENTRACJI MIKROUSZKODZEŃ 45 M oż na się zatem spodziewać, że przy tak duż ej liczbie elementów zachodzą prawa wielkich liczb, które zastosujemy w dalszych rozważ aniach tej pracy. Wielkoś ci mikroskopowe O i makroskopowe Q, a zatem i M bę dziemy w dalszym cią gu traktować jako stał e dla danego oś rodka kruchego, ale nie bę dziemy konkretyzować ich wartoś ci. Bę dą to wię c parametry zależ ne od rodzaju m ateriał u. W zwią zku z t ak przyję tym modelem oś rodka rozróż niamy zniszczenie oś rodka mikro- skopowe, czyli drugiego rodzaju, i makroskopowe, czyli pierwszego rodzaju, przy czym t o ostatnie przy jedn orodn ym stanie naprę ż eń jest równoznaczne z utratą noś noś ci cał ego ciał a. O wytrzymał oś ci materiał u decydują losowo rozł oż one mikrouszkodzenia (mikro- rysy) wewnę trzne, powstają ce w nim w procesie obcią ż enia. Element obję toś ciowy drugiego rzę du traci wytrzymał ość z chwilą pojawienia się w nim jednej mikrorysy. N atomiast pę knię cie elementu pierwszego rzę du nastę puje przy pewnej granicznej liczbie mikrorys r uważ anej za jedną m akrorysę . Liczba n aturaln a r traktowana jest jako specyficzna stał a materiał owa. Z akł adam y, że ciał o o obję toś ci V znajduje się w stanie mikronaprę ż eń 5 i że s k , k = — 1, 2, ..., L , oznaczają naprę ż enia panują ce w elementach mikroskopowych O k . Jeż eli granica mikro wytrzymał oś ci elementu O k wynosi P k , to mikroskopowy warunek pę knię cia tego elementu zapisujemy w postaci (1.3) s k >P k - Jest to warunek wystarczają cy dla powstania stacjonarnej mikrorysy. Jeś li w &- tym ele- mencie m ikroskopowym nie był o mikrorysy, to jest to również warunek konieczny. Wa- runek mikrospójnoś ci m a wówczas postać nierównoś ci przeciwnej (1.4) s k
1, która charakteryzuje kruche zachowanie się
materiał u, jest podstawową hipotezą tej pracy. N ie jest to hipoteza nowa, gdyż m. in.
S. D . WOŁ K ÓW [14], ja k już wspomniano, w ten sposób wł aś nie uogólnił koncepcję naj-
sł abszego ogniwa w ł ań cuchu. Istnieje wiele argumentów przemawiają cych za celowoś cią
wprowadzenia takiego uogólnienia. D o nich zaliczyć należy powszechnie znane fakty
doś wiadczalne, że pojedyncze mikropę knię cia rozcią ganego oś rodka kruchego, wykrywane
n p . przez osł uchiwanie stetoskopem, nie powodują cał kowitego rozerwania ciał a. Prze-
mawia za tą hipotezą także fakt, że materiał y kruche n a granicy wytrzymał oś ci z reguł y
wykazują odstę pstwa od prawa H ooke'a, mimo że odkształ cenia plastyczne typu poś lizgo-
wego lub lepkiego nie wchodzą w rachubę z uwagi n a niskie naprę ż enie i krótki czas próby.
Jeś li doś wiadczalnie stwierdza się w tym przypadku, że odkształ calność graniczna charak-
teryzuje wytrzymał ość m ateriał u, t o ze statystycznej teorii plastycznoś ci destrukcyjnej [9]
wynika, że m aksym alna koncentracja spę kania jest specyficzną stał ą oś rodka, a stą d z kolei
wniosek, który bę dzie jeszcze w tej pracy szczegół owo wyprowadzony, że istnieje liczba r
46 JAN U SZ M U RZ EWSKI, JÓZEF SOJKA
ograniczają ca ilość m ikrorys w jedn ym makroelemencie. N adm ien ić jedn ak należ y, że
hipoteza o stał ej liczbie r nie jest jedyn a, która pozwala wytł umaczyć fakt, że granica wy-
trzymał oś ci może przewyż szać granicę proporcjonalnoś ci ciał a kruchego. D la m ateriał ów
kruchych o silnie nieliniowym prawie odkształ cenia m oż na in terpretować utratę wytrzy-
mał oś ci jako niestateczność stanu naprę ż enia [9], ską d wcale nie wynika r = con st.
W ogólnym przypadku zarówno m ikro wytrzymał ość P ja k też m ikron aprę ż en ie s
są zmiennymi losowymi. D la uproszczenia obliczeń i przejrzystoś ci rozważ ań tylko jedn ą
z tych wielkoś ci bę dziemy uważ ać za zmienną losową o skoń czonej wariancji, a drugą —
ustalimy. Kwestia, którą wielkość przyją ć za losową przy wyprowadzan iu statystycznego
kryterium pę knię cia, czy mikrowytrzymał ość P, jak czę sto zakł ada w swych pracach pierw-
szy z autorów, n p. [7, 9], czy też — mikronaprę ż enie s, ja k to czyni S. D . WO Ł K Ó W [14],
pozostaje otwartą . W pracy rozpatrzymy obydwa warianty, jakkolwiek bardziej uzasadn io-
nym wydaje się przyjmowanie za zmienną losową mikrowytrzymał oś ci, gdyż rozrzut m ikro-
naprę ż eń zależy gł ównie od zmiennoś ci mikroskopowych m oduł ów sprę ż ystoś ci, które n ie
podlegają tak duż ym fluktuacjom jak cechy wytrzymał oś ciowe.
Z akł adam y, że rozkł ady prawdopodobień stw m ikron aprę ż en ia s i logarytm u m ikro-
wytrzymał oś ci In P są typu gaussowskiego. P rawo G aussa bywa z reguł y przyjm owane
dla scharakteryzowania rozkł adu m ikron aprę ż eń lub mikrodefektów w pracach teoretycz-
nych na tem at efektu skali [5, 14]. Potwierdzają je ostatn io przeprowadzon e m ikroskopowe
badan ia eksperymentalne [4]. Logarytmiczno n orm aln e prawo rozkł adu dla m ikrowytrzy-
mał oś ci jest czę ś ciej stosowane niż n orm aln e prawo rozkł adu, gdyż usuwa prawdopodobień -
stwa ujemnych wartoś ci granicy mikrowytrzymał oś ci, nie mają ce sensu realnego. P oza tym
róż nice mię dzy normalną i logarytmiczno norm alną funkcją rozkł adu przy nieduż ych
współ czynnikach zmiennoś ci są bardzo m ał e. W ogóle m oż na mówić, że rozkł ad n orm aln y
i log- normalny jest rozkł adem asymptotycznym dla wielu symetrycznych i niesymetrycznych
rozkł adów prawdopodobień stw [11] i stosowanie jego m a w pierwszym przybliż eniu zn a-
czenie do pewnego stopnia uniwersalne.
2. Wpł yw skali na granicę sprę ż ystoś ci
Korzystamy z nastę pują ch oznaczeń:
V obję tość ciał a,
O obję tość elementu makroskopowego,
O obję tość elementu mikroskopowego,
L ilość elementów O w ciele o obję toś ci V,
N ilość elementów O w ciele o obję toś ci V,
M ilość elementów O w elemencie o obję toś ci Q,
I ogólna ilość mikrorys w ciele o obję toś ci V,
n losowa ilość makropę knięć w ciele o obję toś ci V,
m losowa ilość mikrorys w elemencie Q,
m ś rednia ilość mikrorys w elemencie Q,
r graniczna ilość mikrorys w elemencie O,
mię dzy którym i zachodzą zwią zki:
(2.1) L = M- N ,
(2.2) Q = M- O,
(2.3) V = N - Q = M- N - O = L - O,
(2.4) l=m- N .
E F E K T SKALI P R Z Y KON C E N TR AC JI M I KR OU Ś Z KOD Z EŃ 47
G eometryczne prawdopodobień stwo m ikrorys X, czyli ś rednia koncentracja zniszczo-
nych m ikroelem entów O w obję toś ci ciał a V wyraża się wzorem
W pracach [7, 9] propon uje się nazwać prawdopodobień stwo X spę kaniem. Spę kanie
m ateriał ów czysto- kruchych, nie wykazują cych efektów cią gliwoś ci (plastycznoś ci poś li-
zgowej itp.), równ a się wytę ż eniu w probabilistycznym sensie tego sł owa.
Przyjmujemy najpierw za zmienną losową m ikronaprę ż enie s i charakteryzujemy je
n orm aln ym prawem rozkł adu N (s, / A), O wartoś ci ś redniej ~s i odchyleniu standardowym fi.
Wartość ś rednia s równ a się w przybliż eniu m akroskopowem u naprę ż eniu gł ównemu a,
co wynika z prawa akcji i reakcji [9, str. 268],
(2.6) . * = r = T K a -
Powyż szą przybliż oną równ ość wprowadzam y dlatego, że X w naszych rozważ aniach jest
liczbą bardzo m ał ą, pon ieważ m jest liczbą skoń czoną, a M jest bardzo duż e. Wariancja
m ikron aprę ż eń w oś rodku quasi- jednorodnym wg twierdzenia S. D . WOŁ KOWA [14, str. 29]
jest proporcjon aln a do wł aś ciwej energii sprę ż ystej, a więc
gdzie E
o
> 0 jest stał ą m ateriał ową, okreś lają cą stopień niejednorodnoś ci sprę ż ystej
oś rodka.
Energia poten cjaln a odkształ cen ia sprę ż ystego dla oś rodka izotropowego w rozpatry-
wanym jedn oosiowym przypadku n aprę ż eń wynosi
(2.8) 0 = j - ,
gdzie E jest ustalon ym m oduł em Youn ga. A więc odchylenie stan dardowe mikronaprę ż eń
wynosi
(2.9)
Wprowadzając sym bol v dla oznaczenia współ czynnika zmiennoś ci dostajemy wzór
(2.10) v = - ~- s= *|/ - £• = co n st .
W pierwszym warian cie obliczeń traktujemy mikrowytrzymał ość P
k
mikroelementu a
k
ja ko jedn akową wielkość dla każ dego k,
(2.11) P
k
cn P = c o n st .
Biorąc p o d uwagę (1.4), (2.6) i (2.11) oraz rozkł ad n orm aln y, wyraż amy prawdopodobień-
stwo m ikrospójnoś ci, czyli koncentrację m ikroelem entów niezniszczonych nastę pują cym
wzorem :
(2.12) 9(s P .
W ogólnym przypadku natomiast
(2.16) / = o dla a ! dla er 5
EFEKT SKALI PRZY KONCENTRACJI MIKROUSZKODZEŃ 49
Symbolem Q oznaczyliś my granicę sprę ż ystoś ci. Rozgranicza ona stan doskonale
spójny, liniowo sprę ż ysty, od mikrozarysowania, czyli plastycznoś ci destrukcyjnej materia-
ł u. D la skoń czonej liczby L , ilość mikrorys /, a także granica sprę ż ystoś ci Q są zmiennymi
losowymi. M edianę Q oblicza się z równania
(2.17) ( l - Ai )L = y ,
gdzie \ wyraża się wzorem (2.14) dla a = Q, a modę Q — metodami podanymi p o raz
pierwszy przez T. A. KON TOROWĄ i J. I. FREN KIELA [5]. D la wielkiej liczby L — M- N ,
wg (2.1) i (2.2), i ustalonej ś redniej iloś ci mikrorys I — 1- L wyprowadzamy wzór asympto-
tyczny, z którego okreś limy kres górny Q', realizowany z prawdopodobień stwem równym
jednoś ci,
(2.18) l - ( l - J t ; )L - > l, stą d X[- >\ ,
N admieniamy przy tym, że dla wielkiej liczby L wartość ś rednia At dą ży do tej samej gra-
nicy
(2.19) Ź L, =L j h(l- X)L- 1dX =
1 J_
L +l ~ L
A wię c teoretyczna granica liniowej sprę ż ystoś ci, czyli granica proporcjonalnoś ci, jest
znikomo mał a dla ciał duż ej obję toś ci, jak wynika z równań (2.14) i (2.18),
(2.20) g ' - >0 dla L ~ > o o .
W praktyce okreś la się Q" > Q', jako naprę ż enie konwencjonalne, wywoł ują ce umowną
odchył kę As od prawa H ooke'a (n p. As. = 0,0002). Ponieważ przy plastycznoś ci destruk-
cyjnej odkształ cenie nieliniowe wywoł ane jest redukcją powierzchni przekroju spójnego
(przenoszą cego naprę ż enia) [9], wię c praktycznie
(2.21) l"=%, gdzie « " - X.
S .Er
G raniczny warunek proporcjonalnoś ci formuł ujemy jak nastę puje:
(2.22) %\ = l- Fit') lub t' = W {\ - X[),
gdzie t' — —^ , a W (x) jest funkcją odwrotną do dystrybuanty G aussa, tzn . y = W (x)
v
gdy . = F(y), (rys. 2).
Wartoś ci x bliskie jednoś ci, odpowiadają ce duż ym, niestabelaryzowanym wartoś ciom
argumentu y, wyznacza się ze wzoru asymptotycznego [6, str. 595], rys. 2,
(2.23)
Oznaczają c symbolem Q'
o
granicę proporcjonalnoś ci próbki o obję toś ci V
o
oraz
4 Mechanika teoretyczna
50 JAN U SZ M U RZ EWSKI, JÓZEF SOJKA
przy czym t'
o
= W (l~P^ ), dostajemy zależ ność granicy proporcjonalnoś ci Q' od obję toś ci
ciał a V,
(2.24) Q' = Q
exp(- t
2
m
- « —« - 20 - Zł0 _ z - 4 - ff ^ f l HÓ - « ^
Rys. 2
Zależ ność (2.24) dla róż nych wartoś ci param etru v i przykł adowej wartoś ci X'
o
= l/ L =
= 10~18 przedstawiono na rys. 3.
N atomiast konwencjonalna granica sprę ż ystoś ci Q", okreś lona równaniami (2.21)
i (2.14) dla a == Q", nie zależy od obję toś ci ciał a.
1,08
1,04
1,00
0,98
no?
Hi
r—
ti
Qo 1+VV(1- A'O
uia
)
|
*
- 3 - 2 0
Rys. 3
3
ig v/ v
0
Przyjmijmy teraz w drugim wariancie obliczeń, że m ateriał quasi- jednorodny charak-
teryzuje się zmienną losowo mikrowytrzymał oś cią P, a mikronaprę ż enie s
k
mikroelementu
O
k
moż na w przybliż eniu zastę pować w peł ni okreś loną wartoś cią,
(2.25) s
k
K ? « a,
EF EKT SKALI PRZY KONCENTRACJI MIKROUSZKODZEN 51
i zał óż my, t ak jak t o p o d an o w p . 1, że rozkł ad zmiennej losowej P opisuje prawo loga-
rytm iczn on orm aln e
(2.26)
gdzie P oznacza m edian ę m ikrowytrzym ał oś ci, v logarytmiczny wskaź nik zmiennoś ci,
P o = vP « vP jest odchyleniem stan dardowym .
Jeż eli uwzglę dnimy m ikroskopowy warunek pę knię cia (1.3), to dystrybuanta okreś lona
wzorem,
(2.27)
a
- 0,30-f
- 2
- - 12
- 18
f
//
f - <
H
i
/
/
/
/
l\
ł
i
I/
1
\
1
Ą4
/
—-—
//
iS
C i
1/
j
/
1/
\\
1
i
1
*̂ -
G
\
i
1
P
!
—»•
V
Rys. 4
okreś la prawdopodobień stwo m ikrorys dla n aprę ż en ia a. P orównują c ze sobą formuł y
(2.5) i (2.27) otrzymujemy zależ ność mię dzy n aprę ż en iami a i koncentracją m ikrode*
fektów I,
(2.28) A = l - F ( 0,
przy czym stan daryzowan a zm ien n a losowa m a w tym przypadku postać,
_ In P/er
v
52 JAN U SZ M U RZEWSKI, JÓZEF SOJKA
Zależ ność I od