Ghostscript wrapper for D:\BBB-ARCH\MTS67\MTS67_t5z1_4_PDF\mts67_t5z3.pdf M E C H AN I K A TEORETYCZNA I STOSOWANA 3, 5 (1967) ROZ WIĄ Z AN IA OSOBLI WE W O G Ó LN E J T E O R I I PŁYT TRÓJWARSTWOWYCH RYSZARD G AN O W I C Z (POZN AŃ ) 1. Wstę p W pracy niniejszej podan e zostan ą rozwią zania nieograniczonej pł yty trójwarstwowej obcią ż onej kolejno sił ą skupion ą oraz m om en tem skupionym . Rozwią zania takie nazywać bę dziemy w dalszym cią gu pracy rozwią zaniami osobliwymi. Analogiczne rozwią zania zn an e są w teorii pł yt cienkich [1, 2 i 3] oraz w teorii pł yt Reissnera [4]. Zastosowanie ich jest bardzo szerokie. W zagadn ien iach technicznych rozwią zania osobliwe są wy- korzystywane przy budowan iu powierzchni wpł ywowych [5 i 6]. Znajdują one także zastosowanie przy rozpatrywan iu problem ów niecią gł ych warun ków brzegowych teorii pł yt. P oza tym w teorii sprę ż ystoś ci i w teorii pł yt rozwią zania osobliwe są podstawą uzyskania zwią zków cał kowych mię dzy rozwią zaniami wewną trz danego obszaru a wiel- koś ciami brzegowymi. Zwią zki takie znane są w teorii sprę ż ystoś ci jako twierdzenie So- m igliano. Otrzymuje się je w oparciu o rozwią zania osobliwe oraz w oparciu o twierdzenie 0 wzajemnoś ci. P odobn e zresztą zwią zki znane są w teorii funkcji harmonicznych [7]. Omówienie tych problem ów dla pł yt reissnerowskich i dla pł yt trójwarstwowych o warst- wach zewnę trznych bez sztywnoś ci n a zginanie po dał autor w pracy [4]. W pracy niniejszej zajmiemy się ogólną teorią zginanych pł yt trójwarstwowych podaną przez H O F F A [8]. W teorii tej przyjmuje się , że warstwy zewnę trzne tych pł yt są jednakowe 1 speł niają wszystkie zał oż enia teorii pł yt cienkich. N atom iast odnoś nie do warstwy ś rod- kowej zakł ada się , że jest on a nieś ciś liwa i pracuje jedynie n a naprę ż enia styczne r xz i T yz . Z podan ych wyż ej zał oż eń wynika, że przy zginaniu odkształ cenia pł yty trójwarstwowej są antysymetryczne wzglę dem jej powierzchni ś rodkowej. Rozwią zania otrzym an e w niniejszej pracy są uogólnieniem wyników podan ych przez au t o ra dla pł yt trójwarstwowych o warstwach zewnę trznych bez sztywnoś ci n a zginanie n a ogólną teorię pł yt trójwarstwowych [4 i 9]. 2. Zasadnicze równania Poniż ej podam y zwią zki mię dzy sił ami wewnę trznymi i przemieszczeniami pł yty trój- warstwowej oraz równ an ia równowagi wyraż one przez przemieszczenia. Zależ noś ci te podam y w oparciu o pracę J. WACH OWIAKA i P . WI LD E G O [10]. Sił y wewnę trzne dla om awian ych pł yt podam y za wyż ej wymienionymi autoram i rozkł adają c je n a sił y tarczowe i n a sił y pł ytowe w warstwach zewnę trznych oraz n a sił y tną ce w wypeł nieniu. 294 RYSZARD G AN OWICZ Sił y wewnę trzne okreś limy przez przemieszczenia n astę pują co: (rys. 1) sił y tarczowe N x = Ed 8u 1 Rys. 1 , N y = (2.1.1) sił y pł ytowe (2.1.2) sił y tną ce w wypeł nieniu (2.1.3) Ed I dv - v2 \ 8y n 8x1' Ed 8u • + • m y = - D '• w; + ROZWIĄ ZAN IA OSOBLIWE W TEORII PŁYT TRÓJWARSTWOWYCH 295 gdzie wprowadzono oznaczenia: E, v stał e materiał owe warstw zewnę trznych, G w = G xz = Gyz moduł odkształ cenia postaciowego warstwy ś rodkowej D = —— — sztywność na zginanie warstw zewnę trznych, w ugię cie pł yty, jednakowe dla wszystkich warstw, u, v przemieszczenia w pł aszczyź nie ś rodkowej warstwy dolnej (równe co do bezwzglę dnej wartoś ci przemieszczeniom w pł aszczyź nie ś rodkowej warstwy górnej, lecz przeciwnie do nich skierowane). W dalszym cią gu pracy zajmiemy się pł ytami trójwarstwowymi poddanymi dział aniu obcią ż enia normalnego p(x, y) oraz obcią ż onymi sił ami n x {x, y), n v {x, y) w warstwach zewnę trznych (rys. 2). Zał oż ymy, że te ostatnie są równomiernie rozł oż one na gruboś ci nx\ | z Rys. 2 warstw zewnę trznych i są przeciwnie skierowane w obu tych warstwach. D ział anie tych sił odpowiada wię c dział aniu n a cał ą pł ytę momentów zginają cych, rozł oż onych w obszarze pł yty. Postę pując podobn ie ja k w wyż ej wymienionej pracy [101 otrzymamy nastę pują cy ukł ad równań dla pł yty poddanej dział aniu obcią ż enia p{x, y), n x (x, y) i n y (x, y): (2.2) v 2 ) 2 2 G w (2h+6) 2 (l~v 2 ) ] G w (2h+d) (l- v*) Idu 8v \ _ l- v 2 4Edh • J 2ESh \ 8x~ t 8y)~ 2EÓ P'Eb G w (2h+d)(l- v 2 ) dw [ 82 1- y d2 G w (l- v 2 )] l+v 8 2 v _ \ ~v 2 dx + [dx 2+ 2 By2 Edh J " + 2 8x8y ~ Ed "*' [ Gw(l- v 2)] 2Edh dx + [dx 2+ 2 By2 Edh J " + 2 8x8y ~ Ed - v 2 ) dw l+v 8 2 u \ 8 2 l- v 8 2 G w (l- v 2 )] l- y 1 " dy + 2 8xdy+l8y2+ 2 dx2 Edh \ V ~ Ed M)I>2Edh dy^ 2 8xdy^ \ _8y2^ 2 8x2 Edh J Ed Powyż szy ukł ad trzech równ ań róż niczkowych moż na sprowadzić do równań na trzy funkcje przemieszczeń [10 i 11]: (2- 3) [ l - ^ d - i - j y j ^ ^ j m - . ^, i = 1,2,3, 296 RYSZARD G AN OWICZ gdzie gi = p(,x, y), qi =• 2nx(pc, y), q* = 2ny(x, y), Ed(2h+Sf 1- 2D oznacza sztywność cał kowitą pł yty trójwarstwowej,2 ( 1- *J ) ESh współ czynnik podatnoś ci wypeł nienia. N atom iast przemieszczenia wyznaczyć m oż na ze zwią zków w = ( 1 —) 2h+8 8 2 1y 2h+6 8 8x (2. 4) G w 2 2h+d 8 8y Łatwo zauważ ymy, że jeż eli przyjmiemy, iż warstwy zewnę trzne pł yty trójwarstwowej nie mają sztywnoś ci n a zginanie (D = 0), to zależ noś ci powyż sze uproszczą się do zależ- noś ci podan ych przez autora w poprzedn io cytowanej pracy [4]. N admienimy pon adto, że przypadek ukł adu równ ań jedn orodn ych (2.3) został prze- dyskutowany w poprzednio cytowanej pracy [10]. 3. Obcią ż enie silą skupioną Zajmiemy się nieograniczoną pł ytą trójwarstwową obcią ż oną sił ą skupioną prosto- padle do powierzchni ś rodkowej. Obcią ż enie dla tego przypadku przedstawimy nastę - pują co: p(x, y) = Pó(x)d(y), n x (x, y) = n y (x, y) = 0. ROZWIĄ ZAN IA OSOBLIWE W TEORII PŁYT TRÓJWARSTWOWYCH 297 Wobec tego, że mamy do czynienia z pł ytą nieograniczoną , interesować nas bę dzie cał ka szczególna ukł adu równ ań (2.3). Przyjmiemy wię c F 2 ~ F 3 = 0. U kł ad równań (2.3) i zależ noś ci (2.4) uproszczą się w tym przypadku do postaci: (3.1) w = ( l —«V' L * 2h+d 8 ̂ L 2 |2 ^ Łatwo zauważ yć, że w zwią zkach powyż szych moż na wył ą czyć operator [10] Ostatecznie otrzymamy wię c nastę pują ce równanie (3.2) DM- x~VW VU = P&(x)d(y). oraz zwią zki mię dzy przemieszczeniami w, u, v a funkcją przemieszczeń V (3 3) W = (1- KV2)U « = 2 / i + ó 8U v- 2h+Ó 8U 2 8x ' 2 by Przejdziemy teraz do rozwią zania naszego problem u. Wykonajmy na równaniu (3.2) nieskoń czoną podwójną transformację F ouriera [12]: — 00 —CO 2n J — 00 cc N(x )- - = —[ f 2n J J (3 . 4 ) 0 0 CO Biorą c pod uwagę , że "p*(tx,, /3) == P / 2JT, otrzymamy nastę pują ce wyraż enie na trans- formatę funkcji przemieszczeń: (3- 5) t/ * ( 8) P o wykorzystaniu zależ noś ci (3.4) funkcję przemieszczeń U(x, y) przedstawimy w po- staci cał kowej: (3.6) U(x,y) = CO co c o co 298 RYSZARD G AN OWICZ Aby przedstawić rozwią zanie U(x,y) (3.6) w postaci wyraź nej, bę dziemy musieli obliczyć cał kę wystę pują cą p o prawej stronie powyż szego wyraż enia. P rzekształ ć my tę cał kę nastę pują co: CO co cc p [ r r e~ Kax+ m ID r (3.7) U(x, y) = - j- ~ 2— I I —r~i—~52i&dtt,df}—x—=—- — c o — c o — c o — c o c o c o — Hax+fly) _| ^ ^ I 2D Di J . 1_Lv (r/ 2- Z _ 0 0 - CO l ^ K p. \ a I f J D wie pierwsze cał ki, wystę pują ce p o prawej stronie wyraż enia (3.7), nie istnieją ja ko cał ki niewł aś ciwe, nie m oż na też wydzielić z nich wartoś ci gł ównej wedł ug C auchy'ego. N ależy je rozumieć w sensie czę ś ci skoń czonej [13, 14 i 4], Wydzieleniem czę ś ci skoń czonej cał ek rozbież nych tego samego typu zajmował się autor w cytowanej poprzednio pracy [4]. Aż eby nie rozszerzać niniejszego opracowan ia, podam y poniż ej jedynie koń cowe wyniki tych obliczeń. — 00 oo co —Hax- \ - py) o i o _ x 2 +y 2 — CO — CO to co „ —i'(ra t- fiy) = —2JI In 1 — CO —CQ Ostatnią cał kę wystę pują cą p o prawej stronie wyraż enia (3.7) ł atwo obliczymy ja ko cał kę niewł aś ciwą: Jeż eli weź miemy pod uwagę , że [15] 2xD t o otrzymamy u. (3.10) i?1 = - ^i- .yr f 2KD J d a = J gdzie K(,(z) oznacza zmodyfikowaną funkcję Bessela drugiego rodzaju, y\ = \ D Z \ 2KD. ROZWIĄ ZAN IA OSOBLIWE W TOERII PŁYT TRÓJWARSTWOWYCH 299 Wstawiają c zależ noś ci (3.8) i (3.10) do zwią zku (3.7) oraz wprowadzają c jedn ostkę dł ugoś ci ra otrzym am y poszukiwan e rozwią zanie osobliwe funkcji przemieszczeń U(x, y) Reszta zadan ia jest już bardzo prosta. Znajdziemy teraz rozwią zania osobliwe prze- mieszczeń. Otrzym am y je wykorzystują c rozwią zanie osobliwe U(x, y) (3.11) oraz zależ- noś ci (3.3) W = - —: ~ . x2+y2 P(2h+d) [ 8 , , - 5 —r\ x\ n ~—\ - x\ - ±——^ - r— I n i / x 2 + / + L Ą J AnD z y\ \ _8x y (3.12) v = — Z auważ my, że jeż eli w zależ noś ciach powyż szych przyjmiemy D = 0, 2hj2h- \ - 6 = 1, to otrzymamy po d an e przez a u t o ia we wcześ niejszej pracy [9] rozwią zania osobliwe dla pł yty trójwarstwowej o warstwach zewnę trznych pracują cych jedynie na sił y tarczowe. P odam y jeszcze rozwią zan ia dla sił wewnę trznych w omawianym przypadku nieogra- niczonej pł yty trójwarstwowej obcią ż onej sił ą skupioną . P o wykorzystaniu rozwią zania dla funkcji przemieszczeń U(x, y) (3.11) oraz zależ noś ci (3.3) i (2.1) otrzymamy przy- kł adowo Ed 2h+d (8 2 U , r82U\ P A i - ,* 2 + 2n{2h+S) 300 RYSZARD G AN OWICZ 8 2 U 8 2 U Id* 82 \ 2 1 8y21 J8x2 dy2 \ 8x2 (3.13) = - - ~ -̂ 67T D z |_ ''o ) - O - ") cos 2tp [- ̂ - —K, (y, r)J J, s m 2 p ( 1 " ) s fc_ _ gdzie oznaczono r2 = x2Ą - y2, sin
0 Jak wiadomo, rozwią zanie dla takiego obcią ż enia otrzymuje się przez róż niczkowanie rozwią zania otrzymanego dla obcią ż enia sił ą skupioną . Wobec tego nie wnosi on o zasadni- czych nowoś ci do rozwią zania podan ego w pun kcie 3 niniejszej pracy. Ciekawsze natom iast jest obcią ż enie m om en tem skupion ym zdefiniowanym jako p ara sił poziomych przył oż onych w warstwach zewnę trznych pł yty trójwarstwowej (rys. 5). Rozwią zanie pł yty obcią ż onej t ak zdefiniowanym m om en tem jest bardziej zł oż one od poprzedn io omówionego i może być podstawą uzyskan ia zwią zków cał kowych mię dzy przemieszczeniami u, v wewną trz danego obszaru a wielkoś ciami brzegowymi [4]. Poniż ej zajmiemy się nieogianiczoną pł ytą poddan ą dział an iu m om en tu skupion ego zdefiniowanego jako p ara sił przył oż onych w warstwach skrajnych. Z ał oż ymy, że sił y ROZWIĄ ZAN IA OSOBLIWE W TEORII PŁ YT TRÓJWARSTWOWYCH 303 te są równom iernie rozł oż one n a gruboś ci warstw zewnę trznych. Obcią ż enie tego typu om ówion e został o w pun kcie 2 niniejszej pracy i uwidocznione jest w postaci obcią ż eń n x i n y w równ an iach (2.2) i (2.3). Przyjmiemy do rozwią zania obcią ż enie m om entem skupionym M x = n x (2h+5) — N 5(x)d(y) {2h+ó). N Rys. 5 P odobn ie ja k przy rozwią zywaniu pł yty nieograniczonej obcią ż onej sił ą skupioną interesować nas bę dzie cał ka szczególna ukł adu równ ań (2.3) przy qi = ą % — 0, qi = = 2N 6(x)6(y). Przyjmiemy więc F\ = F 3 = 0, a n a funkcję przemieszczeń Fi otrzymamy równ an ie róż niczkowe (4.1) i N atom iast przemieszczenia u, v, w m oż na bę dzie wyznaczyć z zależ noś ci (2.4) w = — u = 2 (2h+6) 2 4 (4.2) v = Wykon am y teraz podwójną nieskoń czoną transformację F ouriera (3.4)~równania (4.1), a p o wykon an iu retransform acji otrzym am y poszukiwaną funkcję F 2 w postaci *<* y) = - i^ Dl 1 1 r 2D n r — 1 T - co - cc ( a2+ / ?2)2 1 + « ^ i ( a 2 + ^ ) 1 + i . x ( l _ v ) ( a 2 + / ? ) (4 . 3 ) L £>* J L 2 J C ał kę wystę pują cą p o prawej stron ie powyż szego wyraż enia ł atwo przekształ cimy nastę- pują co: 304 RYSZARD G AN OWICZ - ((coc f fty) X D z ( , 00 CO 2 / 2Z> 1 ., . J J , , 1D z 2 v ' 2 v y v Trzy pierwsze cał ki wystę pują ce w tym wyraż eniu omówiliś my w rozdziale 3 niniejszej pracy. N atom iast czwartą cał kę wyraż enia (4.4) obliczymy ł atwo postę pując podobn ie jak przy obliczaniu cał ki (3.9). Otrzymamy dla niej wyraż enie (4.5) R2= J J j l f J J j 1 £ ^ (J J i i l / 1 \ / 2 i O2\ K i l ^ I - co - co l - f pc(i_ „ ) (az+ / S2) gdzie y2 = ]/ 2/ x(l—v) . Wykorzystując zależ noś ci (3.8), (3.9), (3.10), (4.4), (4.5) ofaz wprowadzając jedn ostkę dł ugoś ci r 0 otrzymamy szukane rozwią zanie osobliwe dla funkcji przemieszczeń F2(x, y) Rozwią zania osobliwe przemieszczeń u, v, w otrzym am y wykonując n a funkcji powyż- szej dział ania przepisane zależ noś ciami (4.2) N (2h+6)ix x 2 +y 2 , x , 1 |' 2x 2y t x 1) \ z i 0 l y x lx +y y x 2 +y 2 l x r i f ' ']\ w— — |2 . , 2 02 - x \ - y\ 8x8y gdzie K v (z) jest zmodyfikowaną funkcją Bessela drugiego rodzaju. R OZ WI Ą Z AN IA OSOBLIWE W TEORII PŁYT TRÓJWARSTWOWYCH 305 Jeż eli w wyraż eniach powyż szych wykonamy przejś cie graniczne D -» 0, to otrzymamy przemieszczenia osobliwe dla pł yt trójwarstwowych o warstwach zewnę trznych bez sztyw- noś ci n a zginanie (por. [4]). Przejdziemy teraz d o wyznaczenia poszczególnych sił wewnę trznych dla nieograniczonej pł yty trójwarstwowej obcią ż onej m om en tem skupionym. Wykorzystują c zależ noś ci (2.1) i (4.7) otrzymamy nastę pują ce wyraż enia dla sił pł yto- wych w warstwach zewn ę trzn ych: N (2h+6) D \ cos2
- ) N (2h+S) D sin2
ix TpexcjioftH Łix njiacTH H . PaccyMCflemm OC H O- i Ha Teoptw npefljTWKeHHoił H . fl>K. Xocb4>OM3 [8]. ROZWIĄ ZAN IA OSOBLIWE W TEORII PŁYT TRÓJWARSTWOWYCH 307 penieH H H 3a«a« o H arpyM cemm cocpeflOToneHHOH CH JIOH neorpaH iweH H ofi TpexcjioHHoft a TaraKe o H arpyscen uK TSKOH r u ia c T a n u cocpeflOToieH H biM MOMemoiw B03HHKaiomaM OT n a p bi ropn3OH TajibH bix CH JI npmio>KeH H brx BO BHeuiiiHX CJIOH X. Bce penieH H H n o jiy^ eu bi B 3aMKHyTOM BH fle. O H H cpaBHHBaioTCH c OCO6WMH pemennHMH Teopn n I O H - KH X njiacTHHj a TaK>Ke c penieH H am n Teop«H TpexcjioftiiBix njiacTH H , B KOTopbnc BH eniH ue CJIOH He oSjia- flaiOT H3rH6HOH HteCTKOCTŁK). S u m m a r y m SIN G U LAR SOLU TION S I N TH E G EN ERAL TH EORY OF TH REE- LAYER PLATES Isotropic three- layer plates are considered in the paper, the solution being based upon the theory given by N . J. Hoff [8], The solution of an infinite three- layer plate is derived in the case when the load consists of a concentra- ted force and a concentrated couple formed by horizontal forces acting in the outer layers. Ali results are given in a closed form. The results given in the paper have been compared with the singular solutions of the theory of thin plates and the theory of three- layer plates in which the outer layers exhibit n o bending rigidity. Praca został a zł oż ona w Redakcji dnia 26 wrześ nia 1966 r. 3 M echanika teoretyczna