Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési...

20
-1- MISKOLCI EGYETEM GÉPÉSZMÉRNÖKI ÉS INFORMATIKAI KAR TUDOMÁNYOS DIÁKKÖRI DOLGOZAT Ellipszis keresztmetszetű, izotróp rúdba ágyazott, hengeresen anizotróp, körkeresztmetszetű, kompozit rudak Saint-Venant-féle csavarási feladata Gönczi Dávid MSc. 1. évfolyam Alkalmazott Mechanika Szakirány Konzulens: Dr. Ecsedi István egyetemi tanár Mechanika Tanszék Miskolc, 2010

Transcript of Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési...

Page 1: Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési függvény és csúsztatófeszültségek alábbi képletei [1], [2]: 2 2 ... kör keresztmetszetű

-1-

MISKOLCI EGYETEM

GÉPÉSZMÉRNÖKI ÉS INFORMATIKAI KAR

TUDOMÁNYOS DIÁKKÖRI DOLGOZAT

Ellipszis keresztmetszetű, izotróp rúdba ágyazott,

hengeresen anizotróp, körkeresztmetszetű, kompozit

rudak Saint-Venant-féle csavarási feladata

Gönczi Dávid MSc. 1. évfolyam

Alkalmazott Mechanika Szakirány

Konzulens: Dr. Ecsedi István

egyetemi tanár Mechanika Tanszék

Miskolc, 2010

Page 2: Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési függvény és csúsztatófeszültségek alábbi képletei [1], [2]: 2 2 ... kör keresztmetszetű

-2-

1. Összefoglaló

A dolgozat tárgyát az ellipszis keresztmetszetű, izotróp rúdba ágyazott hengeresen

anizotróp, körkeresztmetszetű, kompozit rudak Saint-Venant-féle csavarási

feladatának analitikus vizsgálata adja. Döntően annak az állításnak az igazolásával

foglalkozik, hogy ha az izotróp, ellipszis keresztmetszetű rúdba hengeresen anizotróp,

körkeresztmetszetű rudakat ágyazunk, akkor az izotróp tartomány vetemedési

függvénye nem változik, ha a két komponens rugalmassági állandói bizonyos

kapcsolatban vannak. Az előzőekben megfogalmazott állítás bizonyítását

többféleképpen is el lehet végezni, például Prandtl-féle feszültségfüggvényekkel vagy

a dolgozatban tárgyalt, öblösödési függvényre épített megoldással. A feladat

megoldása a Saint-Venant-féle általános csavarási elmélet ismertetésével kezdődik,

majd koordináta-transzformációkat követően az öblösödési függvény vizsgálata

következik. A közös peremgörbére vonatkozó egyenletek segítségével megállapítható

a kapcsolat a tartományok anyagjellemzői között (csúsztató rugalmassági moduluszok

között). A következő lépésben az ellipszis keresztmetszetű, izotróp rúdba ágyazott

hengeresen anizotróp, körkeresztmetszetű, kompozit rúd csavarási merevségének

levezetése következik az adott esetben. Az utolsó rész az inhomogén tartományokkal

kitöltött ellipszis keresztmetszet vizsgálatát tartalmazza, ami az előzőekben elemzett

eset kiterjesztése.

Page 3: Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési függvény és csúsztatófeszültségek alábbi képletei [1], [2]: 2 2 ... kör keresztmetszetű

-3-

2. Homogén, izotróp anyagú rúd Saint-Venant-féle csavarási feladata

Ebben a fejezetben általánosságban foglalkoztam a homogén, izotróp rudak Saint-

Venant-féle csavarási elméletével, majd a feladatomban szereplő ellipszis

keresztmetszetet külön is megvizsgáltam.

2.1. Általános, homogén, izotróp, prizmatikus rúd csavarása

Vegyünk egy egyszeresen összefüggő keresztmetszetű L hosszúságú prizmatikus

rudat, melynek Mcs csavarónyomatékkal terhelt vázlatát az 1. ábra is szemlélteti.

1.ábra . Csavart egyszeresen összefüggő keresztmetszetű prizmatikus rúd.

Az xyz derékszögű koordinátarendszer O kezdőpontja a z=0 koordinátával kijelölt

keresztmetszeti síkban tetszőleges elhelyezkedésű. Benne értelmezve:

L

z

O

A P

ez

g

x

y ey

ex

Mcs

Mcs

s

n

Mcs

R

Page 4: Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési függvény és csúsztatófeszültségek alábbi képletei [1], [2]: 2 2 ... kör keresztmetszetű

-4-

- a rúd keresztmetszete: A,

- a keresztmetszetet határoló görbe: g,

- a g görbe ívkoordinátája: s,

- R a P anyagi pont helyvektora, ami felírható a bázisvektorokkal:

x y yxR e e , (1)

- az n a normális egységvektor:

x yn n x yn e e , (2)

- és Mcs a rudat terhelő csavarónyomaték.

Ezek alapján az elmozdulásmező az alábbi képletekkel írható le a Saint-Venant

csavarási elméletben [1], [2]:

- az x irányú elmozdulás: u yz , (3)

- az y irányú elmozdulás: v xz , (4)

- a z irányú elmozdulás: ( , )w x y , (5)

ahol a fajlagos elcsavarodás szöge,

( , )x y pedig a csavart keresztmetszet vetemedési (öblösödési) függvénye.

Továbbá a nyúlásokra és a szögtorzulásokra az alábbi összefüggések állnak fenn:

0xu yzx x

(6)

ahol x az x irányú nyúlás, továbbá hasonló módon belátható, hogy

0y z x (7)

adódik a nyúlásokra csavarás során. A szögtorzulásoknál:

( ) ( ( , )) ( )xzu w yz x y yy x y x x

, (8)

( ) ( ( , )) ( )yzv w xz x y xz y z y y

, (9)

( ) ( ) 0xyv u xz yz z zx y x y

, (10)

Ezek után a Hooke-törvényt felhasználva a következőket kapjuk:

Page 5: Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési függvény és csúsztatófeszültségek alábbi képletei [1], [2]: 2 2 ... kör keresztmetszetű

-5-

( )xz xzG G yx

, (11)

( )yz yzG G xy

, (12)

G: csúsztató rugalmassági modulus,

0x y z xy . (13)

Ezen elemekből előállítható a T feszültségtenzor. Majd a mechanikai egyensúlyi

egyenletet írjuk fel ( ,x y A és a térfogati terhelést elhanyagoljuk) [1], [2]:

0yzxz

x y

T 0 , (14)

( )( )0

G xG y yxx y

,

2 2

2 2 0x y

. (15)

ahol ∆ a Laplace operátor, ami a két Hamilton-féle differenciál operátor (x y

-

x,y síkban) skaláris szorzataként adódik. Majd a rúd palástjának terheletlenségét

kifejező egyenletek alapján azt kapjuk, hogy ha ( , )x y g

0xz x yz yn n T n p , (16)

amit kifejtve adódik:

x y x yn n yn xnx y

. (17)

Az ( , )x y öblösödési függvény tehát az előbbi peremérték feladatok egyik

megoldása. Ezek megoldásai egymástól csak egy additív állandóban térnek el.

Továbbá a csavarónyomaték és a kapcsolatát leírva:

( )cs yz xzA

M x y dA , (18)

csM S . (19)

Itt S jelöli a keresztmetszet csavarási merevségét. A (18) egyenletbe behelyettesítve

megkapható a keresztmetszeti csavarási merevségnek az öblösödési függvénnyel

kifejezett alakja:

Page 6: Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési függvény és csúsztatófeszültségek alábbi képletei [1], [2]: 2 2 ... kör keresztmetszetű

-6-

2 2( )A

S G x y x y dAy x

. (20)

2.2. Ellipszis keresztmetszetű, izotróp, prizmatikus rúd csavarása

Majd tekintsük az ellipszis keresztmetszetű prizmatikus rudat (2.ábra). Az erre

levezethető ( , )x y öblösödési függvény és csúsztatófeszültségek alábbi képletei

[1], [2]:

2 2

2 2( , ) b ax y xya b

, (21)

2

2 2( ) 2xza yG y G

x a b

, (22)

2

2 2( ) 2yzb yG x G

y a b

. (23)

Ezekből kiszámítható a (20) egyenlet alapján S:

3 3

2 2

a bS Ga b

. (24)

O

y

x

a a

b

b

2. ábra. Az ellipszis keresztmetszetű prizmatikus rúd paraméterei.

Page 7: Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési függvény és csúsztatófeszültségek alábbi képletei [1], [2]: 2 2 ... kör keresztmetszetű

-7-

3. A hengeres anizotróp, homogén, körkeresztmetszetű rúd komponens mechanikai analízise A következőkben a homogén, hengeresen anizotróp, kör keresztmetszetű

rúdkomponensek mechanikai leírásával foglalkoztam, melyet majd az előzőekben

tárgyalt izotróp, ellipszis keresztmetszetű prizmatikus rúdban helyezünk el.

3.1 A rúd geomatriája, elmozdulásmezője

Vegyük tehát az izotróp, ellipszis keresztmetszetű prizmatikus rúdba ágyazott

anizotróp körkeresztmetszetű rúd vázlatát, ahogy a 3. ábra mutatja.

3. ábra. Ellipszis keresztmetszetű, izotróp rúdba ágyazot,t hengeres anizotróp,

körkeresztmetszetű rúddal erősített kompozit rúd keresztmetszete.

A c sugarú anizotróp körkeresztmetszetű rúd helyzetének leírására a P0

középpontjának x0,y0 koordinátái szolgálnak. A kör egyenlete így felírható a

következő egyenlettel (és a körlap az ellipszisen belül kell hogy elhelyezkedjen):

O

P0(x0,y0) ex

er

φ

P(r,φ)

ruv

b

a

x

y

c

Page 8: Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési függvény és csúsztatófeszültségek alábbi képletei [1], [2]: 2 2 ... kör keresztmetszetű

-8-

2 2 20 0( ) ( )x x y y c . (25)

Az x,y derékszögű koordináták és az r,φ hengeres koordináták között fennálló

kapcsolatokat az alábbi egyenletek fejezik ki (3. ábra):

0 cosx x r , (26)

0 siny y r . (27)

Az anizotróp rúd helyzetét leíró helyvektor:

0 0 0OP x y 0 x yR e e

. (28)

További helyvektorok:

x y x yR e e , (29)

0 rP P r e

. (30)

Az r,φ henger koordináta-rendszer egységvektorait er és eφ jelöli (3.ábra). Majd

tekintsük a Saint-Venant elmélet alapján az elmozdulásokat a kétféle koordi-nátával

felírva. Először a derékszögű koordinátarendszerből indulunk ki [1],[2]:

0( sin )u yz z y r , (31)

0( cos )v xz z x r , (32)

( , )kw r . (33)

Majd hengerkoordináta-rendszerre (rφz) váltunk és felírjuk a megfelelő

elmozdulásokat (3.ábra):

0 0( sin cos )ru z x y , (34)

0 0( sin cos )v z y x zr . (35)

3.2. Az alakváltozási tenzor elemei, feszültségtenzor

A lineáris rugalmasságtan hengerkoordináta-rendszerben érvényes egyenletei [1], [2],

[3] alapján azt kapjuk, hogy a nyúlások értéke zérus:

0r

rur

, (36)

1 0r vur r

, (37)

Page 9: Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési függvény és csúsztatófeszültségek alábbi képletei [1], [2]: 2 2 ... kör keresztmetszetű

-9-

0z

wz

. (38)

A szögtorzulások közül zérus egyedül a r értéke

0 00 0

1

( sin cos )1 ( sin cos ) 0.

rr

v vur r r

z y x rz y x zr r

(39)

A másik két szögtorzulás:

0 0( , ) ( sin cos )kr

rzr z x yuw

r z r z

0 0( , )( sin cos ).k r x yr

(40)

0 0( , ) ( ( sin cos ) )1 1 kz

v r z y x zrwr z r z

0 0( , )1( sin cos )k r y x r

r

. (41)

A képletek egyszerűsítésére bevezetünk egy segédfüggvényt:

0 0( , ) ( sin cos ) ( , )kr r x y r . (42)

Ezáltal a fenti képletekre -(40) és (41)- a következő alakok adódnak:

( , )rz

rr

, 1 ( , )( )z

r rr

. (43)

Alapvetően kétféle anizotrópiáról beszélhetünk, a Cartesian-féle anizotrópiáról és a

nem Cartesian-féle anizotrópiáról. Ez utóbbi egyik esete a hengeres anizotrópia [3].

Jelen feladatban a hengeresen anizotróp, körkeresztmetszetű rúddal foglalkozunk. A

vonatkozó anyagtörvényeket felhasználva az előbbi nyúlásokból és szögtorzulásokból

levezethető, hogy

0r z r . (44)

A hengeresen anizotróp, rugalmas test csavarási feladatának megfogalmazásához

szükségünk van a Grz és Gφz csúsztató rugalmassági moduluszokra (z normálisú síkon

Page 10: Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési függvény és csúsztatófeszültségek alábbi képletei [1], [2]: 2 2 ... kör keresztmetszetű

-10-

az adott r és φ irányokban). Ezeket felhasználva a csúsztató feszültségeket az alábbi

alakban írhatjuk fel [3]:

( , )rz rz

rGr

, (45)

1 ( , )( )z zrG r

r

. (46)

Ezen tagokból tevődik össze hengerkoordináta-rendszerben a feszültségtenzor ( T ).

Írjuk fel rá a mechanikai egyensúlyi egyenletet (feltéve hogy a térfogati terhelés

elhanyagolható):

T 0 1 0zrz rz

r r r

. (47)

Behelyettesítve a (47) egyenletbe a (46) és (45) kifejezéseket az alábbi

differenciálegyenletre jutunk: 2 2

2 2

( , ) ( , ) 1 1 1 ( , )( ( ) 0rz rz zr r rG G G

r r r r r

2

\rz

rG

2 22 2

2 2

( , ) ( , ) ( , ) 0r r rr r kr r

. (48)

itt bevezettük az alábbi jelölést:

2 z

rz

Gk

G . (49)

Page 11: Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési függvény és csúsztatófeszültségek alábbi képletei [1], [2]: 2 2 ... kör keresztmetszetű

-11-

4. A feladat megoldása az öblösödési függvény használatával

Tehát olyan megoldásokat keresünk, melyek során az izotróp, ellipszis alakú

keresztmetszetekben tetszőlegesen elhelyezett anizotróp, körkeresztmetszetű

tartomány nem változtatja meg az izotróp tartomány vetemedési függvényét. Azaz az

izotróp tartományra vonatkozó képleteknek változatlanul kell maradniuk (a 2.2.

fejezetben tárgyalt képletek). Először kapcsolatot teremtünk az ágyazó és az ágyazott

mátrixok között, majd meghatározzuk a kompozit rúd csavarási merevségét.

4.1. Kapcsolat az izotróp és anizotróp tartományok anyagjellemzői között

A (26), (27) képleteket felhasználva felírjuk az öblösödési függvényt r és φ változók

függvényében. Így a (21) egyenlet alapján:

2 2 2

0 0 0 02 2( , ) ( cos sin sin 2 )2

b a rr x y y r x ra b

, (50)

majd ezt a (40) egyenletbe behelyettesítve a következő alakokat kapjuk:

2 2

0 0 0 02 2( ( cos sin sin 2 ) sin cos )rzb a y x r x ya b

. (51)

Ezekután a Hooke törvény felhasználásával az izotróp tartományban felírjuk a ( rz )

csúszatófeszültséget:

2 2 2 20 02 2 (( ) sin 2 2 sin 2 cos )rz rz

GG b a r b x a ya b

. (52)

Keresnénk a (48) parciális differenciálegyenletnek a megoldását, felhasználjuk hozzá

az (50) és (52) kifejezéseket. A megoldást az alábbi formában próbáljuk előállítani:

2 2

0 0 1 2 32 2

1( , ) ( ( ) cos ( )sin ( )sin 2 )2

b ar x y h r h r h ra b

. (53)

Ezt beírva a (48) differenciálegyenletünkbe azt kapjuk, hogy

Page 12: Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési függvény és csúsztatófeszültségek alábbi képletei [1], [2]: 2 2 ... kör keresztmetszetű

-12-

2 22

0 0 1 2 32 22

2

2 2

0 0 1 2 32 2

2 22

0 0 1 2 32 22

2

22 22 1

2 2

1( ( ( )cos ( )sin ( )sin 2 ))2

1( ( ( ) cos ( )sin ( )sin 2 ))2

1( ( ( ) cos ( )sin ( )sin 2 ))2 0

( )( (

b a x y h r h r h ra br

rb a x y h r h r h ra br

rb a x y h r h r h ra bk

h rb a ra b

2232

2 2 2

31 2

21 2 3

( )( ) 1cos sin sin 2 )2

( )( ) ( ) 1( cos sin sin 2 )2

( ( ) cos ( )sin 2 ( )sin 2 )) 0.

h rh rr r r

h rh r h rrr r r

k h r h r h r

Az előző egyenlet alapján eljuthatunk a hi=hi(r) (i=1,2,3) függvényeket tartalmazó

közönséges differenciál egyenletekhez:

2

2 22 0i i

i id h dhr r k hdr dr

0 ,r c (54)

ahol k12=k2

2=k2; k32=4k2. (55)

Az Euler-típusú differenciálegyenlet r=0 helyen korlátos megoldása:

( ) iki ih r C r (i=1,2,3). (56)

Ezeket visszahelyettesítve az (53) kifejezésbe kapjuk

2 2

20 0 1 2 32 2

1( , ) ( cos sin sin 2 )2

k k kb ar x y C r C r C ra b

. (57)

A (42) és (57) egyenletek kombinálásával az alábbi eredmény nyerhető az anizotróp

tartomány öblösödési függvényére: 2 2

20 0 1 0 2 0 32 2

1( , ) ( cos sin sin 2 )2

k k kk

b ar x y C r y C r x C ra b

0 0cos siny r x r . (58)

A képletekben szereplő Ci (i=1,2,3) állandókat az elmozdulásmező és a

feszültségmező folytonosságának feltevéséből kiindulva számítjuk ki. Ugyanis a közös

határgörbén az izotróp és az anizotróp tartományok elmozdulás és feszültség-értékei

meg kell hogy egyezzenek. Tehát, ha r=c (határgörbén vagyunk) és 0 z L (a L

hosszú rúdon belüli z koordinátára):

( , , ) ( , , )u x y z u x y z , ( , , ) ( , , )v x y z v x y z . (59)

Page 13: Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési függvény és csúsztatófeszültségek alábbi képletei [1], [2]: 2 2 ... kör keresztmetszetű

-13-

Ezen feltételek teljesülnek hiszen a határgörbe pontjainak elmozdulását ugyanazon

függvények írják le.

( , ) ( , )kr r , ( , ) ( , )rz rzr r 0 2 . (60)

A fenti egyenlet teljesüléséhez az alábbi 3 egyenletnek kell teljesülnie (az (50) és (58)

egyenletek összevetésével):

1. 2 2 2 2

1 0 0 02 2 2 2kb a b aC y c y c y c

a b a b

, (61)

2 1

1 2 2

2 ka cCb a

, (62)

2. 2 2 2 2

2 0 0 02 2 2 2kb a b aC x c x c y c

a b a b

, (63)

2 1

2 2 2

2 kb cCb a

, (64)

3. 2 23

kC c c , (65)

2(1 )3

kC c . (66)

Az anizotróp tartományban fellépő rz csúsztató feszültséget ezekkel:

2 22

0 0 1 0 2 0 3 0 02 21( ( cos sin sin 2 ) cos sin )2

k k k

rz rz

b a x y C r y C r x C r y r x ra bG

r

= 2 1 1 2 1 1 2 2 2(1 )0 02 2 ( 2 cos 2 sin ( ) sin 2 )k k k k krzG ka c r y kb c r x k b a c

a b

. (67)

Majd a folytonosság ( , ) ( , )rz rzr r , (0 2 , )r c feltétele miatt összevet-jük

egymással a (67) és (52) egyenleteket. Az előírt feltételek csak úgy teljesülnek, ha a

csúsztató rugalmassági modulusokra fennáll a (68) összefüggés.

rzG kG . (68)

Amiből k-t visszaírva az adódik, hogy

rz zG G G . (69)

Page 14: Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési függvény és csúsztatófeszültségek alábbi képletei [1], [2]: 2 2 ... kör keresztmetszetű

-14-

4.2. A csavarási merevség meghatározása

Ezekután az (58) képletbe visszaírjuk a kiszámított Ci állandókat, s az így kapott

öblösödési függvény:

2 2 2 2 2 10 0 02 2

1( , ) ( ) ( ) 2 cosk kk r b a x y a b r a c r y

a b

2 2 2 1 2 2 2(1 ) 20

1( ) 2 sin ( ) sin 2 )2

k k k ka b r b c r x b a c r .(70)

Majd vesszük a csavarási merevség számításához a (18) képletet. A keresztmetszetet

felbontjuk az izotróp és anizotróp tartományokra vonatkozó részekre:

( ) ( ) ( )k k

cs yz xz yz xz yz xzA A A

M x y dA x y dA x y dA . (71)

Majd a (24) egyenlet alapján felírjuk a kifejezés első tagját a már tárgyalt ellipszis

keresztmetszetre: 3 3

2 2( )yz xzA

a bx y dA Ga b

. (72)

Majd vizsgáljuk a (22), (23) alapján a középső tagot:

2 2 2 22 2

2( ) ( )k k

yz xzA A

Gx y dA a y b x dAa b

=

4 4

2 2 2 20 02 2

2 ( ) ( )4 4

G c cc x b c x aa b

2 2 2 22

2 0 02 2( 2 )

2x b y acG c

a b

. (73)

A 3. tag kifejtéséhez szükséges csúsztatófeszültség:

xz yz rz zd z x y r φτ e e e e . (74)

Az ezek alkotta erőrendszert a P0 pontba redukáljuk. Így keletkezik egy eredő erő ( F )

és egy hozzá tartozó nyomaték(,cs zM ). Az eredő erő meghatározásához

A következő kifejezésből indulunk ki:

( ) ( )k k k

zA A A

dA dA dA z zτ R τ R τ R

Page 15: Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési függvény és csúsztatófeszültségek alábbi képletei [1], [2]: 2 2 ... kör keresztmetszetű

-15-

( )k k k

xyxz

A A A

dA dA dAx y

z zR τ τ F , (75)

mivel a R 1 diádikus szorzat értékét és a (14) egyensúlyi egyenletet használtuk ki.

Továbbá a felülvonás mutatja, hogy a differenciál operátor melyik elemekre hat. Majd

a Gauss-féle integrál átalakítási tételt használjuk ki.

2

0

( ) ( )k k

rzA A

dA ds c cd

z z 0 rτ R τ n R R e

2 2 2 2

2 2 2

0 0 0 0

( cos ) ( sin )rz rz rz rzc d c d c d c d

0 r x yR e e e , (76)

ahol rn e . Továbbá a (67) képlet alapján rz φ szerinti integrálásakor az eredetileg

cosinusos tagok 0, a sinusos tagok az integrálás határai miatt kiejtik egymás, így

fennáll:

2

0

( , ) 0rz c d

, (77)

cos sin r x ye e e . (78)

Majd a (67) egyenletet tovább használva kiszámoljuk F eredő erő komponenseit:

2 2

02 2

0

2cosrz rzka yd G

a b

, (79)

2 2

02 2

0

2( , )sinrz rzkb xc d G

a b

. (80)

Majd a (75),…, (80) alapján az eredő erő kiszámítása:

2

2 20 02 2 ( )rz y

cF G a y b xa b

xe e . (81)

Míg a nyomatékot megkapom a (70) képlet használatával

2 2

2 2

0 0 0 0

1 ( , )( )k

c c

cs z z zA

rM r dA r drd G r r drdr

4

2z

cs

G cM . (82)

Ezáltal a hengeresen anizotróp körkeresztmetszet pontjaiban működő csúsztató

feszültségek nyomatéka a z tengelyre

Page 16: Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési függvény és csúsztatófeszültségek alábbi képletei [1], [2]: 2 2 ... kör keresztmetszetű

-16-

0 0( )

k

yz xz cs y xA

x y dA M x F y F (83)

alakba írható fel, ahol

2 2

02 2x

a c yF Ga b

, 2 2

02 2y

b c xF Ga b

. (84)

Behelyettesítve ezeket a (83) egyenletbe azt kapjuk, hogy

4 2 2 2 2

0 02 2( ) ( 2 )

2k

zyz xz

A

G c a c y b c xx y dA Ga b

. (85)

Majd a csavarónyomaték előzőekben kiértékelt tagjait visszaírjuk az eredeti (71)

képletbe, az alábbi összefüggésre jutunk: 42 2 2 2 2 2 2 23 3 2

2 0 0 0 02 2 2 2 2 2( 2 ) ( 2 )

2 2z

cs

G cx b y a a c y b c xa b cM G G c Ga b a b a b

3 3 4

2 2

( 1)( )2cs

a b k cM Ga b

. (86)

Amiből a fenti egyenletből a (19) összefüggés segítségével megkaphatjuk a keresett

csavarási merevséget:

3 3 4

2 2

( 1)( )2

a b k cS Ga b

. (87)

A képletet megvizsgálva azt mondhatjuk, hogy a csavarási merevség nem függ az

anizotróp tartomány helyzetétől (P0(x0,y0)-tól), csak annak geometriájától (c

sugarától), k jellemzőjétől és az ágyazó izotróp anyag paramétereitől (a,b). Ugyanez az

eredmény a [4]. tanulmányban is megtalálható kör alakú keresztmetszetre (a=b).

Page 17: Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési függvény és csúsztatófeszültségek alábbi képletei [1], [2]: 2 2 ... kör keresztmetszetű

-17-

5. Kiterjesztés több anizotróp tartomány esetére

Az előző fejezet eredményeit felhasználva tekintsük azt az esetet mikor n számú

hengeresen anizotróp, kör keresztmetszetű rudat helyezünk el az ellipszisünkben

(4.ábra).

y

b

x

b

a a

4. ábra. Az izotróp, ellipszis keresztmetszetbe ágyazott n darab anizotróp tartomány.

Az előzőekben tárgyalt feltételek teljesülnek az anizotróp elemekre. Ilyen például,

hogy a csúsztató rugalmassági modulusaik között fennáll, hogy

jrz j zG G G (j=1,2,…,n), (88)

j zj

jrz

Gk

G (j=1,2,…,n). (89)

Ekkor az előző fejezetben tárgyalt csavarási merevség képlete a következő:

3 3

42 2

1

1( )

2

nj

jj

ka bS G ca b

. (90)

Majd töltsük ki az ellipszis keresztmetszetet inhomogenitásokkal. Ha feltesszük, hogy

kj=k (j=1,2,…,n), vagyis valamennyi inhomogenitás azonos anyagból készül, akkor a

csavarási merevség képletében kihozható a k a szummajel elé. Majd az izotróp

c1 P1

c2 P2

cn Pn

cj Pj

Page 18: Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési függvény és csúsztatófeszültségek alábbi képletei [1], [2]: 2 2 ... kör keresztmetszetű

-18-

tartomány anizotróp kör alakú rúdkomponensekkel történő kitöltése miatt feltesszük,

hogy n :

2 2 22 2

1

1( ( ) )2

n

j jj

ab kS G a b p qa b

. (91)

Itt bevezettük az alábbi jelöléseket:

jj

cp

a , j

j

cq

b (j=1,2,…,n). (92)

Az izotróp tartomány nagy számú anizotróp kör alakú rúd komponensekkel töltjük ki

úgy, hogy n és max cj0. Ekkor

2

, 0 1

limj

n

jn c jA ab c

. (93)

A (92) és (93) egyenletek kombinálásával nyerjük a (94) egyenletet.

2

1 1 1 1j j j j j j j

j j j jA c ap bq ab p q A p q

(94)

A fenti egyenletből következik, hogy

1

1j jj

p q

max 0, max 0j jp q . (95)

Tekintsük az alábbi egyenlőtenséget:

2

1 1

( ) max( ) max maxj j k k j j j jj j

p q p q p q p q

, 0 k n . (96)

Ennek felhasználásával a következő eredményre jutunk:

3 3

2 2, , 0lim

j jn p q

a bS Ga b

. (97)

Page 19: Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési függvény és csúsztatófeszültségek alábbi képletei [1], [2]: 2 2 ... kör keresztmetszetű

-19-

6. Felhasznált irodalom

[1] Sokolnikoff, I. S. (1956). Mathematical Theory of Elasticity (2nd edition)

Malabar, FL: Krieger Publishing Company.

[2] Timoshenko, S. P. and Goodier, I. N. (1970). Theory of Elasticity. New York:

McGraw-Hill.

[3] Lekhnitskii, S. G. (1981). Theory of Elasticity of an Anisotropic Body.

Moscow, USSR. Mir Publishers.

[4] Y. Benveniste and T. Chen (2003). The Saint-Venant torsion of a circular bar

consisting of a composite cylinder assemblage with cylindrically orthotropic

constituents. International Journal of Solids and Structures (40), page 7093-

7107.

Page 20: Ellipszis keresztmetszetű, izotróp rúdba ágyazott ... · levezethető ( , )xyöblösödési függvény és csúsztatófeszültségek alábbi képletei [1], [2]: 2 2 ... kör keresztmetszetű

-20-

Tartalomjegyzék:

1. Összefoglaló .......................................................................................................................2 2. Homogén, izotróp anyagú rúd Saint-Venant-féle csavarási feladata ....................................3

2.1. Általános, homogén, izotróp, prizmatikus rúd csavarása...............................................3 2.2. Ellipszis keresztmetszetű, izotróp, prizmatikus rúd csavarása.......................................6

3. A hengeres anizotróp, homogén, körkeresztmetszetű rúd komponens mechanikai analízise 7 3.1 A rúd geomatriája, elmozdulásmezője...........................................................................7 3.2. Az alakváltozási tenzor elemei, feszültségtenzor..........................................................8

4. A feladat megoldása az öblösödési függvény használatával ..............................................11 4.1. Kapcsolat az izotróp és anizotróp tartományok anyagjellemzői között........................11 4.2. A csavarási merevség meghatározása.........................................................................14

5. Kiterjesztés több anizotróp tartomány esetére ...................................................................17 6. Felhasznált irodalom.........................................................................................................19