Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Механика сплошной среды. Т.4 Основы механики твёрдых сред

.pdf
Скачиваний:
19
Добавлен:
19.11.2023
Размер:
41.7 Mб
Скачать

§ 2.8. Двумерные задачи для тел с прямолинейной анизотропией

151

 

'

2

 

о

 

 

 

 

Е

 

+ pfa = О,

 

 

 

/3=1

2

 

 

 

 

 

&аа =

 

а 12 = 2 CQQ£\2 ,

 

 

Е

 

 

<

 

/3=1

 

(2.8.47)

 

 

^а(3

(1 /2 ) (Да,/3

^(3,а) > Ol, [3

1 ?2,

 

 

 

=

Е

 

а

 

 

 

Д^ск/З^^ —^пе’

где

и Са — части контура £, соответствующие поверхностям Е„ и Е^-.

 

Задачу (2.8.47)

называют задачей о плоской

деформации в перемеще­

ниях. После решения этой задачи напряжение <73 вычисляют по формуле (2.8.44), где константу AQ находят из соотношения (2.8.46а).

2.8.9. Метод функций комплексного переменного для плоской деформации

Если тело является изотропным, т. е.

П 11 = П22 = П33 =

1/£7,

П12

= П |з = П23 = —и/Е,

П13П

и{

1

+и)

 

 

 

1

13J-J-23

 

.Bfifi = П«; = ——=

В \2 Пщ — П;

 

Е

Е

 

 

 

2G

■33

 

nf

1 -

Л

 

 

 

 

+ -Вбб,

 

-В11 = В 22 = Пп ———=

Е

= В \2

 

 

 

-LJ-33

 

 

 

Ь\2 = —и(\ - и)/Е,

О

О

и,

(2.8.48)

В\

= В 2 =

тогда уравнение (2.8.45) при отсутствии массовых сил (% = 0) становится бигармоническим. Добавив к нему граничные условия (2.8.35) на контуре, получим следующую задачу:

Д Д Ф 0

в Е 1,

 

 

Ф2 2 = р1 +

Ф0 2 ,

4 > 2 1 =

р1 + Фщ на Д«7 >

(2.8.49)

U\ U\e5

U2

U2е

На

 

где

Д Д Ф = Фд 111 + 2Ф д 122+ Ф,2222. Д Ф = Ф .11+ Ф,22> Р[ = tieds, I = 1,2;

о

Са — часть контура £, на которой заданы усилия; Си — часть контура, на которой заданы перемещения.

Для решения задачи (2.8.49) могут быть эффективно использованы мето­ ды функций комплексного переменного, разработанные еще Г. В. Колосовым [33], Н. И. Мусхелишвили [44], А. Лявом [38], согласно которым вводятся комплексная переменная z и ее сопряженная z:

152

Глава 2. Упругие среды с малыми деформациями

2

- .

1

x2=^-(z — z), (2.8.50)

z = x L+ ix*,

z = x l —ix2,

i = л/-—I , x l = - ( z + z),

 

 

2 V~ ' ~7’

~ 2

где i — мнимая единица, и рассматриваются функции комплексной переменной f(z).

В частности, функцию напряжений Ф можно рассматривать как Ф(^) = = Ф(х1,х2), тогда, ввиду определения (2.8.50), ее производная по z связана с производными по х1 следующим образом:

Ф' =

дФ(г)

=

ж дх1

+ ^

дх1

1f

Л

 

 

a7i

*' l e r

e ?

=

2^

| - 1* '2>

 

дФ

 

дх1

дх2

=

1

,

 

\

(2.8.51)

_

 

= ф ,

_ + ф , _

 

- ( ф , , + 1ф 2).

Отсюда находим

 

 

 

 

 

 

 

,/дФ

дФ\

 

,

 

<9Ф

<9Ф

,

 

 

(2.8.52a)

 

л _

аТ + Hi'

Ф<> = 1Vdz

dz )'

 

 

 

, дх1

 

 

 

 

дх2

 

 

- f — ^ = - , ( * п - 1 Ф . 1 2 & + Ц Ф .1 2 - 1 Ф , 2 2 ) ^ =

 

d z \d z J

 

 

dz

 

 

 

 

dz

 

 

 

 

 

=

 

1Ф 12 + 1(Ф,12 -

1Ф,2г) = ^ АФ.

(2.8.526)

Следовательно, бигармоническое уравнение для Ф можно представить в виде

ААФ = Г1

дАФ = о.

(2.8.53)

16

dz2dz2

 

Решением этого уравнения является функция

 

2Ф = z<p(z) + zip\ (z) + £(z) + £i (z)

(2.8.54)

(при дифференцировании (2.8.54) переменные z и z рассматриваются как независимые).

Поскольку сама функция Ф должна быть вещественнозначной, то она должна совпадать со своей комплексно-сопряженной функцией: Ф = Ф, где Ф^) получаем из Ф^) заменой z на z и всех постоянных коэффициентов на сопряженные к ним величины. Например, для линейной функции

Lp{z) = a + bz, Lp{z) = a + bz, Lp{z) = a + bz, (p(z) = a + bz = +{z).

Тогда в (2.8.54) должны выполняться следующие условия: cp\(z) = p(z) и £1(2) = £(z), и сама формула (2.8.54) принимает вид

Ф = zcp(z) + zcp(z) + £(z) + £(z)

(2.8.55)

и называется формулой Гурса. Согласно этой формуле, бигармоническая функция Ф выражается через две функции комплексного переменного p(z)

И Щ ).

Вычисляя первые производные по х 1 от (2.8.55) с помощью (2.8.52), получаем

§ 2.8. Двумерные задачи для тел с прямолинейной анизотропией

153

 

2Фд = zp'(z) + p(z) + (3{z) + <p(z) + zp'(z) + j3(z),

 

2Ф> = i(z<p'(z) + ip(z) + P(z) - ip(z) + zip'(z) -

(3(z)),

(2.8.56)

так как — p(z) = д ф )

Здесь введено обозначение /3(z) = £'(z).

 

a z

a z

 

 

 

Вычислим вторые производные:

 

 

02

1

 

<p'(z)) + 0{z)

+ (3’{z)),

= Ф,11 = -^(z(f"(z) + Zip"(z) + 2(ip'(z) +

ai = Ф>22 =

^ ( - z p '(z )

- Ztp"(z) + 2(y/(z) + tp'(z)) -

0{z) - P'(z)),

(2.8.57)

0-12 = -Ф .12 = ~^(z<p"(z) - zip"(z) + (3\z) - !3'(z)) = Im {zp"{z) + (3'{z)).

Из (2.8.57) вытекают формулы Колосова:

 

 

 

р = СТ1+ ст2 = 2(y/(z) + <^(z)) = 4 Re (ip'(z)),

(2.8.58)

 

cr2 -

ai + 2i <7i2 = 2(ztp”(z) + j3'(z)).

(2.8.59)

Если ввести функцию

 

 

 

q = 41m (<p'(z)) = 2i(p'(z) p(z)),

(2.8.60)

являющуюся мнимой частью p'{z), то несложно проверить, что она удовле­ творяет соотношениям Коши — Римана

ЯЛ =

- Р . 2 ,

Я.2

= Р л -

(2.8.61)

Поскольку

 

 

 

 

Д<? = Ял 1 +

Я,22 =

~ Р Л 2

+ Р , 21 =

О,

Ар = Д(<тц + сг22) = Д(Ф,ц + Ф,22) = ДДФ = О,

то обе функции р и q являются гармоническими. Из (2.8.59) и (2.8.60) следует формула

4<p'(z) = p + iq.

(2.8.62)

Если обозначить действительную и мнимую части самой функции р:

p(z) = V + iQ,

(2.8.63)

то, дифференцируя (2.8.63) по г с учетом (2.8.51) и (2.8.61), находим соот­ ношения между р, q и Р, Q:

p = 4 R i = 4 Q 2,

q = —4'Р2 = 4 Q j.

(2.8.64)

Используя формулы (2.8.64), можно выразить перемещения щ и г<2 через функции комплексного переменного. Для этого воспользуемся соотношения­ ми Коши (2.8.23а), (2.8.236) и (2.8.23г), в которые подставим определяющие соотношения (2.8.10) для изотропных сред:

154 Глава 2. Упругие среды с малыми деформациями

«1,1 =£\ = 2 ^ (0 - U)(J\ - va2 ) - vA0,

uh 2 + «2,1

= V2e6 = ^<712,

«2,2 = £2 = 2 ^ (_гУСГ1+ 0 - u)a2 ) -

vAo, 2 G

=

J (2.8.65)

Здесь учтены соотношения (2.8.48).

Подставляя выражения

(2.8.44), преобразуем соотношения (2.8.65) следу­

ющим образом:

 

 

 

 

 

 

 

 

 

2Gun =

(1 -

v)(cr

1+ сг2) -

02

- нА0 =

(1 -

v)p -

Ф,п -

1уА0,

2G U 2j2 =

(1 -

v)((T\ + 0*2) -

°\

~ v A 0 =

(1 -

v)p -

Ф 22 -

vAo,

 

 

 

G { li\ 2 + u 2p) = Фд2*

 

 

(2.8.66)

Подставляя в первые два соотношения (2.8.66) выражения (2.8.64), полу­ чаем

2 Gu\p = 4(1 - v)V,\ - Фд 1- vA0,

2GU2,2 4(1 —Z/)Q,2 —Ф,22 —VAQ.

(2.8.67)

Проинтегрируем эти соотношения:

 

 

2Gu\ = 4(1 —v)V Фд —vA§xx + /i(x 2),

 

2GU2 = 4(1 —u)Q Ф 2 —UAQX2 + /^(x1).

(2.8.68)

Подставляя выражения (2.8.68) в третье соотношение (2.8.66), имеем

2(1 -

v)(V,2 + Q,i) +

f[(x2) + /Дж1)) = 0,

(2.8.69)

но, в силу (2.8.64), Vf2

+ Qд = 0, и из (2.8.69) находим, что функции f\ и / 2

могут быть только линейными:

 

 

/ i = 2 G ( w i x 2 + г*ю),

/ 2 = 2 G ( - w i ® 1 + и 2о),

 

где сщ, г^ю и г^о — константы, определяющие перемещение тела как жесткого целого. Тогда выражения (2.8.68) принимают окончательный вид:

2Gu\ = 4(1 v)V —Фд VAQX1 + 2G(CJ\X2 + U\Q),

 

2GU2 = 4(1 —u)Q Ф 2 —UAQX2 + 2G(—сщж1+ IA2Q)

(2.8.70)

и называются формулами Лява.

Составляя из этих формул комплексную комбинацию и учитывая то, что

на основании (2.8.51) и (2.8.55)

 

Ф,1 + 1Ф,2 = 2 ^ = ip(z) + ZLp'(z) + j3(z),

(2.8.71)

получаем формулу Колосова для перемещений:

 

2G{u\ + ш2) = (3 —4z/)y?(z) z(p'(z) — fd{z) —

 

UAQZ T 26r(icjiz T (г/до T i^ 2 o))*

(2 .8 .72)

§ 2.8. Двумерные задачи для тел с прямолинейной анизотропией

155

Обратимся к граничным условиям в задаче (2.8.49). Используя представ­ ление (2.8.56) для производных Фд и и образуя комплексную величину Ф 1+ 1Ф,2, находим граничное условие на контуре Са:

( p ( z ) + z i p ' ( z ) + (3{z) = —р2 + ip \ + (Ф01+ 1Ф02)

(2.8.73)

На части контура Cu, где заданы перемещения, с учетом (2.8.72) получаем граничное условие

2^((3 - 4u)ip(z) - zip'(z) - f3(z)) - uA0z - iwiz+

+ (M10 + i«2o) = (u le + iM2e)- (2.8.74)

Таким образом, задача (2.8.49) свелась к нахождению двух функций ком­ плексного переменного p(z) и £f(z), удовлетворяющих граничным условиям (2.8.73), (2.8.74).

Константы и\ и щ о, що, Що находим из условий (2.8.2) закрепления тела в одной точке (XQ,0), число которых с учетом (2.8.41) равно четырем.

Подставляя (2.8.72) в условие (2.8.2а), получаем

2^((3 - A v ) i p ( z o ) - z 0 L p ' ( z 0 ) - (3( Z Q ) - U A Q ZQ ) + \ ш \ z 0 + и10 + ш20 = u l0 + I U Q ,

Що = uo>

где ZQ = XQ + ixg — комплексная координата точки закрепления; uJQ — коор­ динаты заданного вектора перемещений в точке закрепления.

Вычислим ненулевую компоненту ротора вектора перемещений, используя

(2.8.70):

2

{

 

Щ,2 ~ Щ,\ =

^ 0 - ^)(^,2 - Я,0 -

^ ф, 12 + 2сщ.

(2.8.76)

Подставляя в это выражение формулы (2.8.57) и (2.8.64), (2.8.60) для Фд2 и <2д = q/4, а затем получившееся выражение в (2.8.26), находим окончатель­ ное уравнение для сщ:

— Im (—P '(ZQ) + ZQP"(ZQ) + (3'(ZQ)) + 2cu\ = 2щзо.

(2.8.77)

Пример 2.8.1. Рассмотрим применение метода комплексного переменного в задаче Кирше, в которой сечение тела имеет форму прямоугольника раз­ мерами h°{ и и содержит круглое отверстие радиусом а с центром в начале координат (рис. 2.8.3). Поверхность отверстия и две противоположные боковые грани х2 = ±Л°/2 тела свободны от нагрузок, а на другой паре боковых граней х 1 = ±Лд/2 задано растягивающее усилие:

х 1 = =Ь/гд/2: од =

ст12 = 0;

х2 = ± /^ /2: од = <Ji2 = 0;

г2 = (х1)2 + (х2)2 = a2: fii<7i + П2& \2 = 0, ni<Ji2 + Щ(У2 0*

(2.8.78)

§ 2.8. Двумерные задачи для тел с прямолинейной анизотропией

157

t

Оп2

tl

п2

п2

(2.8.85)

=

+ — )>

£(*) = —f (* + - + - - + ) •

Подставляя выражения (2.8.85) в формулы (2.8.57) и используя представ­ ление комплексных чисел в цилиндрической системе координат:

z = re1^ = r(cos0 + isin</>), z z = г 2 , z = z = 2Re(z),

приходим к следующим выражениям для напряжении:

ai =

 

3+

 

2 ^ )c0s4

^

f t ( 1 -

COS 4<^ -

^2 COS 2^) •

Л )

 

v

2г27

г

2rz

Из формул (2.8.86) следует, что напряжения убывают отверстия как (а/г)2.

На поверхности отверстия при г = а имеем

(2.8.86)

с удалением от

г = а: р — (J1+ 02 —0>г +

= 4(1 —2 cos 20).

(2.8.87)

Здесь учтено, что Ji = сг • • Е = од + од +

03 = а г г + Сфф + a z z —

инвариант

при любых ортогональных преобразованиях, а 0-3 = azz.

 

При ф = ± 7г/2 среднее напряжение р имеет максимальное значение, кото­

рое в три раза превышает напряжение од на бесконечности:

 

Р |

г=а

®ФФтах

max

3te,

(2.8.88)

 

0=±7г/2

 

 

 

 

(поскольку 0-2! _ = 0).

 

 

 

 

ф=±%/2

 

 

 

 

 

Эпюры напряжений

од и од вдоль оси

Ох2

(ф = к /2, г = ж2),

а также

эпюра среднего напряжения р вдоль поверхности отверстия приведены на рис. 2.8.3. □

2.8.10. Изгиб моментами

Вернемся к общему случаю цилиндрических тел, рассмотренному в пи. 2.8.2-2.8.7, и снова будем считать выполненным условие 4 из и. 2.8.8. По­ ложим, что на торцах заданы изгибающие моменты М\, М 2 и поверхностная

сила

а крутящий момент МI и поверхностные усилия t\, 4

на боковой

поверхности отсутствуют:

 

 

 

 

M l = 0,

tl = fe = 0,

со30 = 0.

(2.8.89)

Выберем, как и в и. 2.8.8, оси декартовой системы координат совпадаю­ щими с осями ортотропии Oci, тогда в общем решении, представленном в пи. 2.8.5-2.8.7, можно принять

д = 0, Ф = 0, 7 = 0, сщ = 0.

(2.8.90)

158 Глава 2. Упругие среды с малыми деформациями

Напряжения и деформации в теле в этом случае определяются по формулам

<73 = Е 3{ А \ Х Х+ А 2Х2 + А 0), 0\ = <72 =

<712 =

а 13 = &23 = О,

£3 = А \ Х { + А 2Х2 + AQ,

Е\2 =

£13 =

е 23 = О,

£1 = -г/31£3,

£2 = —гУ32£з-

(2.8.91)

Подставляя выражение для <т3 в условия на торце (2.8.37а), (2.8.376), в

силу допущения 4 о симметрии

сечения Е 1 тела относительно осей Ох1

и

Ох2, находим три константы AQ, А\

и А2:

 

 

Ао =

,

 

Ml

,

 

Ml

lae Ха d£, а = 1, 2,

 

А] = —

 

Ао —

Eshe

 

 

ЕЛIS

ЕЛ\е

 

 

 

 

где 1 ае — моменты инерции сечения И1 относительно оси Оха.

 

Если

же подставить

деформации

е\, е2,

£12 = 0 в уравнения (2.8.30),

а

£13 = £23 —0 — в уравнения (2.8.31), а затем проинтегрировать их, то получим следующие выражения для W\, W2 и W3 (определяем их с точностью до перемещений тела как жесткого целого, так как такие перемещения уже учтены в общем представлении перемещений (2.8.32)):

W\ =

\ д ^ - ( х х)2 + А 2х 1х 2 + Л ож 1) + и32^ { х 2)2,

 

w 2 = - щ2 (А\х хх2 + А 2 2 Г - + AQX2) + Vz\^-(x1)2, W3 = 0.

(2.8.93)

Если рассматривать жесткое закрепление в точке О — в начале коорди­

нат, то граничные условия

(2.8.26) на торце ж3 = 0 примут вид

 

u lQ=

?1Q = ?1Q = 0, UJ\о = CJ20 = 0

(2.8.94)

(в точке О исключаются не только перемещения, но и малые повороты как жесткого целого в плоскостях Оххх2, Ох2х 3). Тогда, подставляя (2.8.32) в

(2.8.94), находим

 

 

 

 

 

що = и2о = щзо = 0,

и\ = и 2 = 0.

(2.8.95)

Формулы (2.8.32) совместно с (2.8.93) принимают окончательный вид:

щ = - 4р-(Ч3)2 “

^31Ц 1)2 + VЗ2(ж2)2) -

^(Л гж 'ж 2 + Л0ж'),

 

и2 = - ^ ( ( ж 3)2 -

v3

2{x2 ) 2 + г'зЩ 1)2) -

и3 2{А\ххх2 + Л0ж2),

 

щ =

(A Q + A i x 1 +

А.2ж2)ж3.

(2 .8 .96)

Полученное решение, очевидно, представляет собой суперпозицию реше­

ний трех отдельных задач: изгиба моментом М \ (когда JF3 = 0, М 2

= 0), из­

гиба моментом М 2 (JF3 = 0, М \ = 0) и растяжения силой JF3 (М х =

М 2 = 0).

В частности, при чистом изгибе моментом М\ из (2.8.91) и (2.8.96) полу­ чаем следующие выражения для напряжений, деформаций и перемещений:

160

Глава 2. Упругие среды с малыми деформациями

 

 

 

7 =

# 7 ,

 

 

(2.8.103)

тогда для 7 получим следующую задачу (В44

1/ 2G23, -В55

1/2G13):

 

1

1

 

 

в Е 1,

(2.8.104)

 

2С?гз7 , п

+ 2Gi3 7,22

=

 

 

 

 

 

7 = 0

на

С,

 

 

а для константы $ из (2.8.103) имеем выражение

 

 

д = M l/С е,

Се = 2

7

(2.8.105)

X1 константу Се называют жесткостью тела при кручении.

Напряжения и деформации в задаче кручения, согласно (2.8.14) и (2.8.10), запишем следующим образом:

СГ1 = СГ2 = СГ12 =

0,

СГ13 =

^7)2,

023 =

- $ 7,ь

 

(2.8.106)

А

 

 

А

 

#7,2

 

£23 =

$7,1

£3 = 0,

е 1 = £ 2 =

£12 = 0,

£13

= ^Ст13

,

^Ст23>

а перемещения, согласно (2.8.30) и (2.8.32), —

 

 

 

 

щ = —$ж2ж3,

U2

= Ле1#3,

г^з = rdW{xI).

 

(2.8.107)

Здесь введена новая функция W: W3 = $147, для которой, согласно (2.8.31),

имеем следующую задачу:

 

 

 

 

 

 

 

 

W i = ^ T

+ x 2 ,

W 2

= X L

- X 1

B E 1;

W \

= 0 .

 

(2.8.108)

G13

 

G23

 

 

 

a

 

 

Для изотропных сред задача (2.8.104) представляет собой внутреннюю задачу Дирихле для уравнения Пуассона, при решении которой могут быть эффективно применены методы теории функций комплексного переменного. Для простейших форм сечений И1 тела решение можно получить без приме­ нения специальных методов.

Пример 2.8.2. Для цилиндрического изотропного тела с круговой формой поперечного сечения И1 решение задачи (2.8.104) имеет вид

7 = (G/2)(a2 —г2), г2 = (ж1)2 + (ж2)2 ^ а2, (2.8.109)

где а — радиус кругового сечения. Подставляя это выражение в (2.8.109), убеждаемся, что Wt\ = П72 = 0, и решением задачи (2.8.108) является нулевая функция W = 0.

Вычислим жесткость при кручении кругового цилиндра Се из (2.8.105):

Се = 2 7 dZ = 2TiG (а2 — r2)r dr = TrGa4 $ = 2Ml

( 2.8. 110)

7rGa4

x1