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

книги из ГПНТБ / Каландия А.И. Математические методы двумерной упругости

.pdf
Скачиваний:
22
Добавлен:
24.10.2023
Размер:
10.69 Mб
Скачать

Для нахождения первого приближения по формуле (39.10) потребуются значения рІ^ФоОО1)]- Чтобы вычислить эту функ­ цию, достаточно определить из (39.12) Фо(т}) = arg юо(Ç) и под­ ставить найденный аргумент в правую часть соотношения (39.11) вместо Ф. В результате будем иметь

G (#) =

Л'*0/?1 = ,

6 0 1

(39.14)

Задавая

теперь в интервале

(0, я)

п узлов

 

 

q,

=

2v—1

,

0

 

 

Щ

~2я~ я,

v = l , 2, .

 

построим

четный

 

тригонометрический полином

G„_i (f>) =

л—1

 

 

 

 

 

 

= ^ Лѵ cos vù, интерполирующий

G(f>) в данных

узлах. По-

Ѵ=0

лином этот согласно формуле (39.14) не будет содержать чле­

нов с нечетными индексами. Обозначая, далее, через

G„_i (т>)

тригонометрический полином, сопряженный с Gn-ifö),

найдем

первое приближение искомой отображающей функции в окон­

чательном

виде:

 

 

 

 

 

 

<ві (а) =

(ùg (o)gn-i

(о),

g„-i

(а) = G„_, (а) + iG*+l

(Ü). (39.15)

 

Численные

расчеты были проведены для значений

полуосей

 

2

 

9

 

 

 

 

 

а — —т= , о = -^=г ;

было

взято

д = 1 2 . Для первого

прнближе-

ния было найдено

выражение

 

 

 

 

 

 

 

о (£) = Çe**+e V сг 2 ѵ £ 2 ѵ ,

 

(39.16)

где

 

 

 

 

ЙО = — 0,003604,

 

 

а =

0,127706,

а0 = 0,976450,

а., = 0,022810,

ß =

0,016136,

а6 =0,003460, а8 =0,000703,

а1 0 =0,000134.

 

Кривая, являющаяся изображением окружности по первому

приближению

(39.16), отклоняется от заданного эллипса не бо­

лее, чем на 0,006345, а максимальное отклонение для началь­

ного приближения имеет

порядок 0,048524. Первая

итерация

сближает, таким образом,

изображение

нулевого приближения

к истинному контуру примерно в восемь

раз. Отметим,

что ото­

бражение, найденное по методу Л. В. Канторовича в рассмат­

риваемом примере, отходит от

точного примерно на 0,02

(Л. В. Канторович и В. И. Крылов

[1], стр. 441).

Глава пятая

О ПЛОСКИХ ЗАДАЧА Х НЕСИММЕТРИЧНОЙ УПРУГОСТИ

Согласно концепции

братьев Коссера (Е. Cosserat, F. Cosse-

r a t [1])» учитывающей

вращательное

взаимодействие частиц ма­

териала, при изучении

напряженного

состояния твердого дефор­

мируемого континуума необходимо наряду с обычными напря­

жениями (сила на единицу площади) вводить

в рассмотрение

и моментные напряжения (моменты на единицу

площади).

Наличие моментных напряжений предполагает несимметрич­ ность тензора (обычного) напряжения и приводит к рассмотре­ нию градиентов деформации высоких порядков. В случае плос­ кой деформации изотропного тела 1 ) появляется одна дополни­ тельная упругая постоянная, имеющая размерность длины и свя­ зывающая моменты с вызванными ими искривлениями.

Определяющие уравнения и граничные условия, соответст­ вующие тому или иному состоянию равновесия, имеют более сложный вид по сравнению с классическим случаем. Число гра­ ничных условий (геометрических, статических и смешанных), необходимых для определения упругого равновесия тела, равно трем. Число возможных комбинаций задания этих условий зна­ чительно больше числа тех же комбинаций в классической уп­ ругости.

Интерес к теории моментных напряжений значительно возрос за последнее время и различные аспекты ее стали предметом изучения многих авторов2 ). Основная часть имеющихся к нас­ тоящему времени солидного числа научных публикаций посвя­ щена общим теоретическим вопросам обоснования теории.

') Что мы и будем иметь в виду-во всем дальнейшем.

2 ) Укажем здесь

на работы Мпндлниа

[1—3], Миндлина и Тирстеиа [1],

Трусделла и Туппна

(Truesdell and Toupin

[I]) и Тупина (Totipin [I]). Из ра­

бот советских авторов назовем

работы

Э. Л. Аэро и Е. В. Кувшинского [I]

и Г. Н. Савина [3]. Довольно

полную

библиографию по указанным вопро­

сам можно найти в обзорной статье Г. Н. Савина и Ю. Ы. Немиша [1].

261

В настоящей главе мы приведем основные соотношения не­ симметричной теории упругости и укажем применение аппарата теории функций комплексного переменного к решению гранич­ ных задач в простейших случаях.

§ 40. Основные соотношения плоской моментной теории. Постановка граничных задач

1°. К обычным уравнениям равновесия тела, связывающим четыре компоненты напряжения, присоединяется одно уравнение равновесия моментов Мх и Ми. Система уравнений равновесия в случае отсутствия объемных сил здесь имеет вид')

да

дт,,,.

0,

от,.,,

да,.

п

дМ

дМ„ .

XIJX =

п

0.

, л п ..

-л* +

- Г * =

+

— ^ =

0,

-1-

« - I - гХу -

 

40.1)

дх

ду

'

дх '

ду

 

дх

ду

 

 

 

 

Ввиду несимметричности тензора напряжений касательные напряжения имеют симметричную и антисимметричную состав­ ляющие T s и га:

Симметричная составляющая т.„ так же как и в классическом случае, вызывает обычную деформацию сдвига ъ х у ,

х = 2 ^ х у ,

г х у = ± 2 [дх ' ду

(40.3)

где ц. модуль сдвига,

а антисимметричная

составляющая ха,

уравновешиваемая моментнымн напряжениями, поворачивает элемент на угол ы,

Моментные напряжения Мх, Ми искривляют элемент, вследст­ вие чего появляются пропорциональные, по предположению, этим напряжениям кривизны хх, у.„, определяемые формулами

да

<Эш

 

% х ~~ ~дх~ '

КУ ~ ~ду '

(40.5)

где В — постоянная, называемая

модулем изгиба.

 

За новую постоянную материала, характеризующую

степень

влияния моментных напряжений, здесь принята имеющая раз­

мерность длины величина /, определяемая

равенством

1°- = — .

(40.6)

') Содержание настоящего пункта заимствовано

из работМипдлииа [1, 3].

262

Поле напряжений (обычных и моментных) выражается через две функции напряжений U и V по формулам

d2U

д-Ѵ

 

 

_

d4J_ _,

д^Ѵ_

 

 

а * = ду*

дх ду

'

аУ

~~

дх'- г

дх ду '

(40.7)

d*-U - d-W

' ^ _

д'-Ц ^ dW_

 

дх

ду

ду"-'

V

=

дх

ду 1

дх2 ')

 

 

 

дѴ_

 

 

 

 

 

(40.8)

 

'' дх '

 

У

ду

 

 

 

 

 

 

 

Функции напряжении U и V удовлетворяют

системе

урав­

нений

 

 

 

 

 

 

 

 

д-{Ѵ-РЬѴ)

=

2 ( 1 - Ѵ ) Р А Д І / ,

 

дх

 

 

 

 

 

 

 

 

j-(V-l*AV)

=

2 ( \ - v )

l ^ A U .

 

 

Из этих уравнений легко заметить, что функция U — ангармо­ ническая, а V удовлетворяет уравнению

АѴ—12А2Ѵ—0.

Для упругих смещений и, ѵ и напряжений (40.7), (40.8) имеют место следующие формулы общего комплексного пред­ ставления, установленные Р. Д. Миидлиным [3] ' ) :

2р. (и + іѵ) = хер (г) — zw'

(z)

— яр (z) —

 

 

 

 

 

- 8 ( l - v )

I4p"(z)

+

2 i ^ ,

 

 

 

 

 

 

 

dz

 

ox +

(jy

=

2 [w' (г) +

VW)},

 

• ax + i (rxv

+

x,JX)

=

2 [zw"

(г) - f я|/ (г)

+

 

8(1 - v ) / у " ( г ) + 2/ â 2 -

 

(40.9)

(40.10)

+

= 8 t ( l -v)l2w"(z)

+ 2^L,

x = 3 - 4 v . [ ( 4 0 Л 1 )

 

 

dz

j

Здесь ер и яр обозначают потенциалы Колосова — Мусхелишвили,

2U = z~w(z)+zVJ)

+

X(z)

+

W),

%'{*) = №),

(40.12)

v — коэффициент

Пуассона,

Я — решение уравнения

Гельм-

гольца

 

 

 

 

 

 

 

 

АН

— К*Н = 0,

 

% =

~ ;

(40.13)

а — обозначает

операцию

 

 

 

 

 

dz

 

 

 

 

 

 

 

 

д

 

 

 

д

 

(40.14)

 

д'г

= ~\1х~

1'

1 ду

 

 

 

 

') В названной работе имеются опечатки.

263

Наконец, функция V допускает представление (см. там же)

у = - 4t ( 1 - V) Г- [ср' (г) -

+ H.

(40.15)

Выпишем еще известные зависимости между упругими эле­ ментами в полярной H декартовой системах координат:

ог

+ сто = ах

-!- Oy,

 

 

Т г в — Т д г

 

 

(40.16)

о в

сгг - f і (тг* + тег) = К

— ст* -I-1 хѵ - f т„Л )] e2,'»,

 

 

 

 

 

(40.17)

 

ог

+ iü» = (и +

іѵ) е-''*, 2 = x + iij re1

(40.18)

Ha основании этих формул и формул (40.9) — (40.11) по­ лучим следующие формулы комплексного представления реше­ ний уравнений теории моментных напряжений в полярных ко­ ординатах:

 

ar

- f о* =

2 [ср/ (г) +

ср' (г)],

 

с в -

аг + і (хгь +

т*г ) =

2 [г<р" (г) +

і|/ (г)

+

(40.19)

 

- 8 ( 1 - ѵ ) / Ѵ " ( 2 ) - Ь 2 ^ г ,

 

 

 

 

 

8(1 —v) / У ( г ) + 2i

dz

Ol*

(40.20)

2ц (о, + ш») =

[xcp (2) - 2 Ф ' (г) - xp (2) - 8 (1 -

v) /чр" (г)

+

 

 

 

- f 2t. дНdz

 

(40.21)

2°. Подсчитаем компоненты вектора напряжения Jncfs, Fn ds, действующие на элемент дуги какого-либо профиля AB, прове­ денного в плоскости Оху. Имеем

Х„ =

о* cos пх -•- % ѵ COS tt!/, j

(40.22)

 

 

К,, =

Txy COS ПХ + OyCOS /If/,

 

где n обозначает нормаль к кривой AB, проведенную вправо по отношению к положительному на AB направлению. Замечая, что на AB

^

dij

^

dx

COS ПХ =

—T'­

cos n y —

 

 

as

 

 

из формул (40.22) и (40.7) выводим

A " ~ ds 1 ду

 

V

_

dU

dV

дх

dx

ду

1

" ~~

ds

 

 

 

 

264

Складывая предыдущие равенства, после умножения второго из них на і, и интегрируя полученную сумму по s, получим

2 - І - {U — іѴ) =

if (X„ + iYn) ds + const.

dz

•„

Внесем сюда вместо U и V их выражения (40.12) и (40.15). Тогда получим следующее равенство, справедливое вдоль лю­ бого профиля в физической области: •

<p(z) + z с,/ (г) + -ф (z) + 8 ( 1 -

V) /я ср" (г) -

,.

дН

=

 

=г f (X„ 4- iYn) ds 4- const. (40.23} 'о

Тем же путем, исходя из формул (40.8) и используя, пред­ ставление (40.15) для моментиого напряжения М„, действующе­ го на тот же профиль А В,

находим

 

Мп = Мх

cos юс 4- My cos inj,

 

 

(40.24)

 

 

 

 

 

 

 

 

 

 

 

 

Мп

=

- Re f - ^ -

f 8 (1 -

V) /2

<j?4i) - 2i ^

-

] | .

(40.25)-

 

Приведем,

наконец,

формулу

комплексного

представления

для

угла поворота

со, даваемого

равенством (40.4). На

основа­

нии формулы

(40.9)

имеем

 

 

 

 

 

 

со = Im | —

(и + Ц

= - ^ - Im {(и 4- 1)Ф ' (г) + - ^ г я].

(40.26)

 

Рассмотрим

теперь

интеграл

 

 

 

 

 

 

 

 

'

=

f (Х,« и 4- У„ V + Ма со) ds,

 

 

(40.27)

где L обозначает контур физической области 5, предполагаемой

конечной.

 

 

 

 

 

 

 

 

 

 

 

 

Заменив под интегралом Хп, У„, Мп

их значениями из формул

(40.22), (40.24), преобразуем интеграл

по формуле

Грина и вос­

пользуемся уравнениями

равновесия (40.1). Будем

иметь

,

[' С Г да .

дѵ ,

 

дѵ .

du .

 

 

 

 

 

 

4

(%< -

те д ) со +

М, - ^ - 4- Мѵ Щах dij.

(40.28)

 

Учитывая

(40.4)

и (40.3), находим

 

 

 

 

тх„

+ т,, ^

 

+ (%< -

т ж у ) со =

( - * • 4- - І ^ ) =

Ч^%- (40.29>

Сумму первых двух членов под интегралом (40.28) выразим согласно закону Гука через деформации sx, гѵ, примем во.

265

внимание равенство (40.29), а для двух последних членов в ин­ теграле воспользуемся формулами (40.5). Тогда упомянутый интеграл представится в виде

 

/= 2\\Wdxdij,

(40.30)

 

s

 

 

где

 

 

 

2W == а.(е* + еу )2 +

2р{г1 + гІ + 2г%) + 4ß [4

+ xf,). (40.31)

Выражение W — неособенная

положительная

квадратичная

форма совокупности

компонент

деформаций

(как обычных,

так и вызванных моментными напряжениями),— представляет собой потенциальную энергию упругой деформации, заключен­ ную в единице объема.

Отсюда немедленно следует равенство нулю всех шести ком­ понент напряжения ах, ау, ххѵ, хух, Мх, М„ во всей области, если на

ее границе выполняется одна из двух следующих

групп условий:

X,~Yn=Mn

= 0,

(40.32)

и = о = (й =

0.

(40.33)

По условиям (40.32) решение будет определено

с точностью до

жесткого смещения тела как целого,

 

и = — гу-\-а,

w = e.v4-ß,

 

где е, ос, ß — постоянные.

Задачи определения упругого равновесия тела при неодно­ родных условиях (40.32) и (40.33) будем по аналогии с класси­ ческим случаем называть соответственно первой и второй основ­ ными задачами моментной теории упругости.

Выше мы доказали единственность решений этих задач для конечной области. Приведенное рассуждение с известным огра­ ничением поля напряжений на бесконечности справедливо и для случая бесконечных областей.

Так же, как и в классическом случае, для существования ре­ шения первой основной задачи в конечных областях необходимо выполнение условий статики твердого тела. В нашем случае эти условия имеют вид

f (Хя + iYn) ds = 0,

\ (xYn - УХ„ + М„) ds = 0.

(40.34)

L

L

 

Они непосредственно следуют из уравнений равновесия

(40.1).

Указанные выше комплексные представления полей напряже­ ний, смещений и деформации позволяют теперь сформулировать основные граничные задачи теории моментных напряжений в

терминах теории

функций.

 

 

 

 

 

П е р в а я о с н о в н а я

з а д а ч а . Найти

аналитические

в

рассматриваемой

области

функции

cp(z), ^(z)

от

комплексного

266

аргумента

 

z и

регулярное

в

той же

области

решение

H

уравне­

ния

(40.13)

 

по

граничным

 

 

условиям

 

 

 

 

 

 

 

 

ф (t)

+

t Ф' (0 +

Ф (О -+8(1 — V ) Z* ф " (О -

,. дН

=

 

 

 

21

 

 

 

 

 

 

 

 

 

 

 

 

 

= М О + іМ О + const,

(40.35)

 

 

 

 

 

 

! ( 1 - ѵ ) Р Ф " ( 0 - 2 і | '

= М 0 и а L >

 

 

где

fi,

/г,

fo —

заданные

на

 

контуре

L

вещественные

функции

от

t, определяемые

 

через

внешние

усилия

и

моменты,

действующие

на

границу

 

тела,

 

равенствами

 

 

 

 

 

 

 

 

 

 

 

 

 

fi

+ ih

= i Г (Хп

+

îYn)

ds,

f0 =

- M„.

 

(40.36)

 

Второй

 

основной

 

задаче

 

соответствуют

граничные

 

условия

 

 

 

х Ф

(0 - t

ф' (0 — т|з

(0 — 8 (1 ѵ ) / 2

ф " (0 +

 

 

 

 

 

 

 

 

 

 

 

 

 

2*-=- = 2 | i ( g x

- i - igs ),

(40.37)

 

 

 

 

 

 

 

 

 

 

 

 

а/

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Im

 

(х +

1 ) Ф ' ( 0

т

Ш

=

2j.tg-0

на L ,

 

 

 

 

 

 

 

 

2/2

 

 

 

 

 

 

 

go заданные

 

смещения

 

точек

контура

L ,

g0-

•за-

данные

на

L

значения

 

угла

 

 

поворота:

 

 

 

 

 

 

 

 

 

 

u+iv=gl-\-ig2,

 

 

a — go

 

на

L .

 

(40.38)

Наряду с условиями указанных двух основных задач можно, разумеется, ставить и другие граничные условия, например, задавать на границе области и, ѵ, Мп и т. д. Среди этих других задач мы укажем лишь одну, соответствующую соприкасанию тела с жестким профилем заданной формы при отсутствии сил трения (см. § 1). В качестве третьего граничного условия при­ мем здесь наличие углов поворота точек границы и запишем условия задачи в виде

Г = 0 , vn=?g(t),

a = h(t)

на L ,

(40.39)

где Т — касательная составляющая

вектора напряжения, а о , —

нормальное

смещение;

g(t)

и h(t)

обозначают

заданные

па L

функции. Эту задачу мы будем условно называть третьей

задачей

момеитной

теории упругости.

 

 

 

 

Исходя

из формул

(40.35)

и (40.37), можно

без особого

тру­

да построить эффективные решения основных задач для облас­ тей, ограниченным одной окружностью или двумя концентриче­ скими окружностями. Некоторые простейшие задачи такого рода будут рассмотрены в следующих параграфах.

267

§ 41. Решение первой основной задачи для бесконечной плоскости с круговым отверстием1 )

Центр отверстия радиуса R возь'мем в начале координат. Условия (40.35) без несущественной постоянной в правой части примут в нашем случае вид

Ф(0

+tVTt) + W) + m ц

Щ

=

h •!- i-h,

 

R e ^

imp" {t) -

ІЩ-

= - /„,

t = Re'*,

 

 

 

 

ді

 

 

где принято обозначение

 

 

 

(41.1)

 

 

 

 

 

 

ш = 8(1— ѵ)12.

 

 

(41.2)

Мы

будем считать,

что поле

па

бесконечности

однородно2 )

и, кроме того, главный вектор усилий, приложенных к обводу

отверстия, равен

нулю.

 

 

 

Положим

в области вне контура отверстия

 

 

 

 

со

 

со ,

 

 

ф(г) =

Г г - ; - v a * * " * , а|,(г) =

Г'г - f V ак г~к

(Г = Г),

(41.3)

 

 

i

 

о

 

 

H (x,

у) =

v

Ай/с* (Xr)

г - /хм1 ?2 ,

= hk,

(41.4)

 

 

— со

 

 

 

 

где Кп(Кг) обозначает функцию Бесселя второго рода от мни­ мого аргумента (функция Макдопальда).

Выражая операцию —— в полярной системе координат, при­ меняя ее к правой части выражения (41.4) и используя правило дифференцирования функции Кп (см., например, Ватсон [1]), получим

Тг = ~ 2- S Mft-i Kk frr) е^.

(41.5)

—со

Ряды для производных от функции (p(z) мы получим, оче­ видно, дифференцированием соответствующее число раз по z ряда (41.3) для ф(г). Составим все производные нужных поряд­ ков и положим, кроме того, на окружности

U (0 + Ш

=

Ъ

Ake*\

 

 

— СО

(41.6)

 

 

со

— /о (0

=

2

с-л e'"ft!>, cc_fe = a.k. i

') Задача для спободного от усилий отверстия'в однородном поле была решена в цитированной выше работе Миндлина [1]. Решения некоторых ана­ логичных задач имеются в работах Вейтцмана [I] и Хартренфта и Си [1].

2 Ï JVloMeimibie напряжения на бесконечности везде принимаем равными нулю: М[т) = M œ ) = 0.

268

Все предыдущие ряды будем считать равномерно сходящи­ мися на окружности и в этом предположении внесем их в усло­ вия (41.1). Мы будем иметь после некоторых приведении

со

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y i a

k R - k e - i k u

- y

i ( k - 2 ) R - k + 2

[1 — m (ft -

1)

tf~2] ö A _ 2

е і

к Ь 4-

 

i

CO

 

 

з

СО

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

Jj.^LR~K ß'*9 +

X

Уі

Cfe-i

Kk

(Kr) е«» +2ГRe'» 4- P R e - ' »

=

 

 

 

 

 

 

 

 

 

 

 

 

 

v

A,

 

 

 

 

(41.7)

 

 

 

 

 

 

 

 

 

 

l +

 

±j rt-k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l S

tu IKk+i

(Kr)

+

 

 

 

 

 

 

 

 

 

 

Kk-i

(Kr)]e'*>

= H a t e ' f c » ,

 

 

(41.8)

причем

 

 

 

 

Ск

=

t'/lft,

c_/e =

— cA.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сравнение

в

этих равенствах

коэффициентов при е'п *

даст:

 

 

 

 

 

 

 

a'Q — Яс^о = Л0 ,

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a/CiC0 = a0 f

 

f

 

 

 

 

(41.9)

(2 -

л) Я - » + 2

[1 -

m

(я -

1) Я~2 ] Н/ г _2 4-

и;

Д-»

4-

Яля _, /С„

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= Л„

п > 1 ,

(41.10)

//irt(« - fl)/? - » - 8 «„4 - »,[/C n + 2 +/C«Jc„ + 1

=

- 2 i a „ + 1 >

n >

 

1, (41.11)

 

 

 

 

а„ R-"

- Кс~^Кп

=

Л_„,

я > 1 .

 

 

 

(41.12)

п=\

При такой записи систем уравнений необходимо

при

 

первый

член

в

левой

части

соотношения

(41.10),

содержащий

множителем

(и—2), считать равным

нулю,

а

правые

части

Ai

и Л_і равенств

(41.10)

и

(41.12) —заменить

соответственно

на

Ai—

2TR

и Л_, Г ' R .

 

 

 

 

 

 

 

 

 

ап,

 

 

 

 

Система

(41.11), (41.12) определит неизвестные

с п + ь

/ г ^ І ,

а система

(41.10)—все

неизвестные

 

, начиная

с

/ г = 1 . Не­

достающие неизвестные а'0, с0 , ci определяются из формулы (41.9).

Решая совместно системы (41.11) и (41.12), находим

а,, =

(41.13)

д - (И-і) [2(-а„+ 1 / ? 3 - тп (я 4 1) Л _ „ ]

сп+\ =

269

Соседние файлы в папке книги из ГПНТБ