• Nie Znaleziono Wyników

Optymalne kształtowanie belki wspornikowej, obciążonej siłami zewnętrznymi i ciężarem własnym, w warunkach pełzania

N/A
N/A
Protected

Academic year: 2021

Share "Optymalne kształtowanie belki wspornikowej, obciążonej siłami zewnętrznymi i ciężarem własnym, w warunkach pełzania"

Copied!
8
0
0

Pełen tekst

(1)

M E C H AN I KA TEORETYCZNA I STOSOWANA

4, 16 (1978)

OPTYMALN E KSZTAŁTOWAN IE BELKI WSPORN IKOWEJ, OBC IĄ Ż ON EJ SI Ł AM I ZEWN Ę TRZN YM I I CIĘ Ż AREM WŁASN YM, W WARU N KACH  P EŁZ AN IA*

M AŁ G O R Z AT A AL B I Ń S K A,  A N T O N I  G AJ E W S K I  ( K R A K Ó W ) 1, Uwagi wstę pne

Elementy konstrukcyjne (np. belki, pł yty) obcią ż one wyłą cznie sił ami zewnę trznymi lub wyłą cznie cię ż arem własnym rzadko wystę pują  w praktyce; na ogół  sił y masowe sta-nowią  dodatkowe obcią ż enie, które jest pewn m uł amkiem obcią ż enia zewnę trznego. W wię kszoś ci opublikowanych dotychczas prac, dotyczą cych optymalnego kształ towania belek (z reguł y liniowo- sprę ż ystych) siły masowe, zależ ą ce od poszukiwanego kształ tu elementu, są  pomijane. W szeregu nowszych prac ([1], [2]) podję to próbę  znalezienia optymalnego kształ tu belki wspornikowej, obcią ż onej wył ą cznie cię ż arem wł asnym. Niestety, przedstawione rozwią zania są  bł ę dne. Poprawne rozwią zanie zagadnienia przed-stawiono w pracy J. M. CIIERNA [3]; otrzymano je w oparciu o zasadę  minimum wzajemnej energii potencjalnej R. T. SHIELDA i W. PRAGERA [4], waż nej tylko w przypadkach materiał u liniowo- sprę ż ystego. W pracy A. GAJEWSKIEGO [5] rozwią zano analogiczny problem opty-malizacji kształ tu belki wspornikowej, znajdują cej się  w niejednorodnym polu sił  grawi-tacyjnych (a wię c obcią ż onej tylko cię ż arem wł asnym). Uwzglę dniono również nielinio-wość fizyczną  materiał u, opisują cą  materiał y nieliniowo- sprę ż yste, sprę ż ysto- plastyczne lub pozostają ce w stanie ustalonego peł zania. W pracy wyznaczono optymalne kształ ty belki przy warunku wyrównania naprę ż eń, i warunku - ustalają cym ugię cie swobodnego koń ca belki. Wykazano, że we wszystkich badanych przypadkach, optymalny kształ t zależy od postaci prawa fizycznego.

Próbę  znalezienia optymalnego kształ tu belki obcią ż onej sił ami zewnę trznymi i cię ż a -rem własnym podję to w pracy A. GAJEWSKIEGO [6], w której przedstawiono przybliż one rozwią zania, otrzymane metodą  mał ego parametru. Obliczenia przeprowadzono przy wa-runku wyrównania naprę ż eń w skrajnych wł óknach belki oraz przy warunku minimali-zacji ugię cia swobodnego koń ca belki. Zał oż ono przy tym, że cię ż ar wł asny stanowi mał y uł amek obcią ż enia zewnę trznego.

W niniejszej pracy przedstawimy rozwią zanie analogicznego zagadnienia optymalizacji kształ tu belki wspornikowej, jednak otrzymany warunek konieczny istnienia ekstremum rozwią ż emy numerycznie dla dowolnego udział u cię ż aru wł asnego w obcią ż eniu cał kowi-tym.

ł )

 Praca wykonana został a w ramach problemu wę zł owego 05.12 pt. „Wytrzymał ość i optymalizacja konstrukcji maszynowych i budowlanych", koordynowanego przez Instytut Podstawowych Problemów Techniki Polskiej Akademii N auk.

(2)

500 M. ALBIŃ SKA, A. G AJEWSKI

U wzglę dnimy również nieliniowość fizyczną  materiał u, zakł adają c nieliniową  zależ ność mię dzy naprę ż eniem i odkształ ceniem. Przyję te dalej nieliniowe, potę gowe prawo fizyczne pozwoli opisać materiał y pozostają ce w stanie ustalonego peł zania [7], nieliniowe sprę -ż yste lub sprę ż ysto- plastyczne (bez odcią ż enia). W przypadku materiał ów pełzają cych, w przyję tym prawie fizycznym symbol „ e " oznacza prę dkość odkształ cenia (e), a ugię cie koń ca belki w(0) należy zastą pić prę dkoś cią ugię cia (w).

2. Sformuł owanie zagadnienia

W niniejszej pracy ograniczymy się  do rozważ ania belek wspornikowych, o prostoką t-nym przekroju poprzecznym, wysokoś ci 2h i szerokoś ci b, wykonanych z materiał u nie-jednorodnego, opisanego potę gowym prawem fizycznym:

• liczba naturalna,

(2.1) "~

l

(a\  . da

w którym przyję to, że funkcje: so(x) i <ro(x) są  znane i okreś lają  niejednorodność wytrzy-mał oś ciową materiał u. W dalszym cią gu pracy bę dziemy opuszczali kropkę  nad s.

Rys. 1

Punktem wyjś cia do dalszych rozważ ań bę dzie podstawowa zależ ność prę dkoś ci krzy-wizny belki od momentu zginają cego (lub zależ ność odwrotna), która w przypadku prawa (2.1) przyjmuje postać [6]:

J

gdzie M* oznacza bezwymiarowy moment: (2 3~) M* -  • { } M  ' — u v / \

Równocześ nie zał oż ymy, że belka obcią ż ona jest sił ami rozł oż onymi o intensywnoś ci

p(x), oraz, że cię ż ar wł aś ciwy belki jest równy q(x) (rys. 1). Wówczas moment zginają cy skł ada się  z dwóch czę ś ci: 1). skł adnika pochodzą cego od obcią ż enia zewnę trznego oraz 2). skł adnika zależ nego od cię ż aru wł asnego, a tym samym od poszukiwanego kształ tu:

(2.4) M(x) =

(3)

OPTYMALNE KSZTAŁTOWANIE 501

Celem pracy jest znalezienie takiego sposobu zmiany wysokoś ci belki wzdł uż dł ugoś ci

h(x) (przy stał ej szerokoś ci b), szerokoś ci b(x) (przy stał ej wysokoś ci) lub wysokoś

ci i sze-rokoś ci (gdy są  do siebie proporcjon aln e: b(x) =  ah(x)), aby okreś lona funkcja celu przyjmował a wartość minimalną , przy warunku ograniczają cym, ustalają cym cał kowity cię ż ar belki:

(2.5) /  q(x)F(x)dx =  W  =  const. o

W dalszym cią gu ograniczymy się  do minimalizacji prę dkoś ci ugię cia swobodn ego koń ca belki. Pewne rozwią zania przy innych funkcjach celu otrzym an o m etodą  m ał ego parametru w pracy [6].

Korzystają c ze wzoru (2.2) prę dkość ugię cia koń ca belki moż emy wyrazić nastę pują cym funkcjonał em:

(2.6) w(0)  = ( ^ ± 1 ) " f

x

\  n } J

( ) f

x

^M

\  n } J h{pć )

W dalszym cią gu zał oż ymy, że obcią ż enie zewnę trzne m a taki sam zwrot jak cię ż ar wł asny, tzn. że m om en t nie zmienia swojego znaku i jest stale dodat n i. Wobec tego we wzorach

(2.2)—(2.6) opuś cimy zn ak m oduł u i przyjmujemy: sgnM =   + 1 .

Jako param etr optymalizacji przyjmiemy pole powierzchni przekroju F(x), co pozwoli ujednolicić tok obliczeń dla wszystkich trzech wymienionych przypadków (zmiany1

 sze-rokoś ci, zmiany wysokoś ci i przypadku gdy: b = a/ i). Po prostych obliczeniach funkcjonał (2.6) m oż na przedstawić w postaci: w której należy przyją ć w poszczególnych przypadkach : 1. h => const., b =  b(x) (zmienna szerokoś ć) 2. b = const., h =  h(x) (zmienna wysokoś ć) 3. b(x) = ah{x) (przekroje powinowate) H =  3/ 2, v = 111, C =

(4)

502 M . ALBIŃ SKA, A. G AIEWSKI

3. Warunek konieczny istnienia ekstremum funkcjonał u (2.7)

W celu znalezienia kształ tu belki F(£), minimalizują cego funkcjonał  (2.7), przy izo-peryraetrycznym warunku ograniczają cym (2.5), utworzymy pomocniczy funkcjonał : (A t — mnoż nik Lagrange'a), i obliczymy jego pierwszą wariację: (3.2) <5/* =  C J ^p%Z~lr l"FdM-  (pn + v) MÓF]d£ + Xtl f gdFdt, o ° . o Obliczając wariację SM ze wzoru (2.4): t (3.3) ÓM =  P f q(rj)dE(V)  ( | - v)drj, o

podstawiając ją do (3.2) i wykonując cał kowanie przez czę ś ci, otrzymujemy:

i o

gdzie: £0 i ^l2 są dowolnymi stał ymi.

Pierwsza wariacja funkcjonał u (3.1) jest zatem równa zeru, gdy przyjmujemy, że stał a lo — 1  °r az gdy speł nione jest równanie:

(3

-

5)

* f

które stanowi warunek konieczny istnienia ekstremum funkcjonał u (3.1).

W celu uproszczenia zapisu, zmienność cię ż aru wł aś ciwego belki q(i) opiszemy za pomocą pewnej znanej bezwymiarowej funkcji K(S), moment zginają cy, pochodzą cy od obcią ż enia zewnę trznego Mp($), za pomocą znanej funkcji m (|), pole powierzchni przekroju za pomocą poszukiwanej funkcji 0(i), oraz wprowadzimy bezwymiarową stał ą e wedł ug nastę pują cych zależ noś ci:

(3.6) g =  q0K(£), Mp =  MpQm(£), F te  Ą * ( l ), e -  l 2

(5)

OPTYMALNE KSZTAŁTOWANIE 503

Wprowadzają c wyraż enia (3.6) do równań (3.5), (2.4), (2.5) otrzymujemy odpowiednio:

(3.8) (3.9) =  Mp 6

j

o = W , gdzie:

Ax—jest mnoż nikiem Lagrange'a.

Uwikł ane równanie cał kowe (3.7) i warunek izoperymetryczny (3.9) pozwalają  na wyz-naczenie poszukiwanego kształ tu belki 0(1) i stał ej Au przy czym rozwią zanie to może być znalezione tylko na drodze numerycznej.

4. Numeryczne rozwią zanie zagadnienia

W celu rozwią zania zagadnienia posł uż ono się  metodą  kolejnych przybliż eń, przedsta-wiają c równanie (3.7) w postaci:

(4- 1) 0(S) -  F[0(S)1, w której F oznacza pewien znany operator.

Zakł adają c funkcję  0o( f) przybliż enia zerowego i podstawiają c ją  do prawej strony

równania (4.1) otrzymujemy pierwsze przybliż enie; nastę pne przybliż enia obliczamy ze wzoru:

(4.2) *i + i -   # 0 |( O J , i =  0, 1, 2, . . . ,

przy czym jako kryterium zbież noś ci rozwią zania przyjmujemy nierównoś ć:

(4.3)

obowią zują cą  dla każ dej wartoś ci f.

Ostatecznie równanie (3.7) przekształ cono do postaci:

(4.4) 0KS) =

limit,) +

w której przyję to, że materiał  belki jest wytrzymał oś ciowo jednorodny i ma stał y cię ż ar właś ciwy tzn.: e0 =  const., a0 =  const., J5T(f) m 1. Równanie to w szczególnym przypadku:

(6)

504 M . ALBIŃ SKA, A. G AJEWSKI

m(C) =  £, ju — 1, v =  0, n =  1 przyjmuje postać równania zamieszczonego w pracy

J. M. Cherna [3].

Obliczenia numeryczne przeprowadzono na maszynie CYBER 72; jako dane wejś ciowe programu przyjmowano: 1. bezwymiarowy cię ż ar belki W , 2. parametr e charaktery-zują cy udział  cię ż aru wł asnego w obcią ż eniu cał kowitym (e =  0 — belka obcią ż ona tylko sił ami zewnę trznymi, E - » oo —belka obcią ż ona tylko cię ż arem wł asnym), 3. parametr n charakteryzują cy prawo fizyczne (n =  1 — materiał  liniowo- sprę ż ysty, n - > oo — materiał sztywno- plastyczny), 4. krok cał kowania: Zlf (po wykonaniu testów dokł adnoś ci przyję to:

Ai =  0,01). '

O 0,1 0,2 0.3 0,4 0.5 0.6 0.7 0,8 0,9 1.0 Rys. 2

0,1 0,2  0,3  0,4  0,5  0,6  0,7  0.8  0,9  1,0

Rys. 3

Stosowana metoda obliczeń okazał a się wystarczają co szybko zbież na, co pozwolił o na wyznaczenie poszukiwanego kształ tu z dużą dokł adnoś cią. Przyjmując np. A = 10~* funkcję  0 ( |) otrzymano po oś miu iteracjach; wówczas czas typowego przebiegu nie prze-kraczał  oś miu minut pracy procesora centralnego.

Obliczenia wykonano dla belki o stał ej wysokoś ci i poszukiwanym sposobie zmiany szerokoś ci b(x), tzn. dla: JX — 1 i v =  0, obcią ż onej: 1. sił ą skupioną /1

 dział ają c ą na swo-bodnym koń cu (rys. 2, 3 i 4) i 2. sił ami rozł oż onymi o stał ej intensywnoś ci p(x) =  const (rys. 5 i 6). W obu przypadkach wyznaczono zależ noś ci kształ tu belki od parametru n przy ustalonej wartoś ci e oraz zależ noś ci od parametru e przy ustalonym n.

N a rys. 2 zilustrowano wpł yw postaci prawa fizycznego (wykł adnika w) na optymalny kształ t belki zginanej sił ą skupioną i cię ż arem wł asnym. D la porównania na rysunku zamieszczono również kształ t belki zginanej wył ą cznie obcią ż eniem zewnę trznym; kształ t ten nie zależy od postaci prawa fizycznego ponieważ minimalizacja ugię ci a w punkcie dzia-ł ania sia w punkcie dzia-ł y skupionej jest równoważ na minimalizacji pracy sia w punkcie dzia-ł  wewnę trznych [6]. Podobny charakter posiadają krzywe przedstawione n a rys. 3, przy czym moż na zauważ yć, że przy wię kszym udziale cię ż aru wł asnego w obcią ż eniu cał kowitym wpł yw wykł adnika n w prawie

(7)

OP TYM ALN E KSZ TAŁ TOWAN I E 505 O 0.1  0 . 2  0 . 3  0,4  0 . 5  0,6  0 . 7  0.B 0,9  1,0 Rys. 4 O 0.1  0.2  0,3  0,4  0.5  0.5  0.7  C.6  0.9  1,0 Rys. 5 0,1 0,2  0,3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 Rys. 6

fizycznym, na optymalny kształ t jest znacznie mniejszy. N a rys. 4 zamieszczono wykresy kształ tu optymalnych belek liniowo- sprę ż ystych (n = 1) dla róż nych wartoś ci parametru e. Uwzglę dnienie cię ż aru wł asnego istotnie wpływa na kształ t optymalny, przy czym wy-kresy otrzymane dla wartoś ci e > 10 praktycznie nie róż nią  się  od wykresu dla e = 10.

Podobny charakter mają  rozwią zania przedstawione na rys. 5 i 6, dotyczą ce belki ob-cią ż onej sił ami równomiernie rozł oż onymi. Dla porównania zamieszczono na rys. 5 kształ t belki liniowo- sprę ż ystej bez uwzglę dnienia cię ż aru wł asnego (e =  0). Kształ t ten zależy od wykł adnika n [6],

(8)

506 M . ALBIŃ SJCA, A. G AIEWSKI Literatura cytowana w tekś cie 1. R. L. BARN ETT, Minimum Deflection Design of a Uniformly Accelerating Cantilever Beam, J. Appl. M ech., 30, (1963), 466—467. 2. L. C. W. D I XON , Pontryagin's maximum principle applied to the profile of a beam, J. of the Aeronautical Society 71, (1967), 513—515. 3. J. M . CH ERN , Optimal structural design for given deflection in presence of body forces, Int. J. Solids Structures, 7, (1971), 373—382.

4. R . T. SH IELD , W. PRAG ER, Optimal Structural Design for Given Deflection, ZAM P, 4, 21, (1970), 513— 523.

5. A. G AJEWSKI, Optymalne kształ towanie belki wspornikowej z materiał u nieliniowego fizycznie, obcią ż onej

cię ż arem wł asnym, Rozpr, Inż ., 3, 24, (1976), 453—467.

6. A. G AJEWSKI, Optymalne kształ towanie wytrzymał oś ciowe w przypadku materiał ów o nieliniowoś

ci fi-zycznej, Zeszyty N aukowe Politechniki Krakowskiej, N r 5, Kraków, 1975.

4. Ju. N . RABOTN OW, Poł zuczest elementów konstrukcji. N auka, Moskwa, 1966.

P e 3 IO M e ,

OI 1TH M AJI Ł H OE OOP M H P OBAH H E KOH COJIfcH OH  EAJIKH  H ATPyjKEH H Ofl: BH EIU H H M H  CH JIAM H  H  COECTBEH H BIM BECOM B yC JIOBH H X

OciroBH oń reMoił  paSoTbi HBJIJKTCJI ormiMajiKiioe (toopMHpoBaHiłe KOHCOJIŁHOH 6ajiKH  c npjiMoyrojib-nonepe^mbiM ce^iesmeM, Harpy>KeHHOH  BHeuiHHMH  pacnpe# ejieHHbiMH  cmiaMtr, a Taiofce CHJiaiwM OT coScTBeH H oro Beca. O6cy>KflaioTcii 6anKH  H3roTOBJieHHBie H3 iwaTepnajia HeoRHOpoftHoro n o conpoTH

B-on acaH H oro cieneHHiiiM dpH SK11601011 *1

 3aK0H0M. 3T O T 3ai- coH  HBjrsieTCH  cymecrBeHHŁijw fljui B COCTOHHMH ycTaHOBHBiiieHCH  nojisyqecTH , H ejiH H efiH o- ynpynix «JIH ynpyro- njiacTH MH tix. B paG oTe o6napy>KeH bi onTHMajiŁHtie cpopMbi SajioK c nepeineHHOH  mapHHOH, Harpy>i<eHHWx cocpefloio-tieHHoii Ha KOHqe cHJioii a TaK>Ke CHuaMK pasH OMepno pacnpeflejiSHHbiiwH  c nocTOHHHOH KHTeHCHBHOcTbio.

H ccjiefloBaH O TaKJKe BjmHHHe noKa3aiejiH  „H" creneH H oro 3aKOHa Ha onTHmanbiibie (bopmbi. Bo n p o c pein eH  c HcnojiB3oBaHH:eM KJiaccircecKHx MeTOflOB Bapnai^iroH H oro Hc^mcjieiiHiH, a pe3yju>TaTbi HyMepHMeaaCc pacieTOB m uiwcrpH poBaiibi pa3JHWHWMH  pncyHKaMH.

S u m m a r y

OPTIM AL D E SI G N  OF  T H E CAN TILEVER BEAM LOAD ED  BY EXTERN AL FORCES AN D  BY ITS OWN  WEIG H T I N  CREEPIN G  CON D ITION S

There was investigated th e optimal design of the cantilever beam with rectangular cross- section, loaded by the uniformly distributed external forces and by its own weight. There were taken into account beams of the inhomogeneous material, described by the power physical law. This law works in the steady state creep and for the materials non- lineary elastic and elastic- plastic. There were found the optimal shapes of beams with varying width, loaded by the force acted on the end of the beam and by the uniformly distri-buted force with constant intensity. There was also examined the influence of an exponent ,,n" in the power law on the optimal shape of cantilever beam. The problem was solved on the basis of classical variational calculus. The results of numerical calculations are presented on numerous figures. POLITECHNIKA KRAKOWSKA INSTYTUT FIZYKI

Cytaty

Powiązane dokumenty

Figures 6, 7 and 8 show the impact of the height, asymmetry angle of the spraying, sampler speed and its position on the degree of spray coverage for the XR 11003 VK spray

10 min, aż osadzająca się miedź przekroczy pewną średnicę, można zaob- serwować skręcanie się coraz bardziej narastających gałęzi w przeciwną stronę, ale już

Stwierdzono, że uwzględnienie sprężystego zamocowania krawędzi ścianki (płyty) wsporniko- wej w segmencie pręta cienkościennego oraz poprzecznej i wzdłużnej zmienności

Marek Cała – Katedra Geomechaniki, Budownictwa i Geotechniki co potwierdza znaczący wpływ kształtu stropu na parametry obudowy kotwiowej koniecznej dla utrzymania stateczności

W ten sposób stosując analizę klasyczną (tzn. niekwan- tową) otrzymaliśmy taką samą wartość i kierunek orbitalnego momentu magne- tycznego, jak w podejściu kwantowym. Być

W anali- zie zastosowano funkcję wielokwadratową, opisano metodę doboru wartości parametru kształtu oraz podano współczynniki korygujące wyniki analizy w przypadku

Analizując drgania belki, własności warstwy MR można opisać, wykorzystując zespolony moduł odkształceń poprzecznych, opisujący tak zmiany jej sztywności jak i

linia ugięcia belki, kąt obrotu belki, warunek sztywności przy zginaniu, równanie różniczkowe linii ugięcia belki, warunki brzegowe, zastosowanie zasady