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

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

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

приведении

4 (to +

?о) N0 (t0) + ЬѢ=Ы

j

HfMi

+

 

 

 

 

 

 

 

V

0

 

 

 

 

 

+

(2p0-?o)jWo jM*+ R e ^[qoFo( g _

 

где

 

 

 

-

/"of i

( Ш = h (t0)

на у,

(20.11)

 

 

 

 

 

 

 

 

 

A„ =

M-i («о I) llo (Xi 1),

Po =

ï*i («о 4- I) +

h> («a4-1), 1

<7o =

l l i («o +

1),

 

mo =

p-o К

+ 1 ) ,

 

1

 

 

 

 

 

 

 

 

 

(20.12)

a F 0 и Fy — известные функции, даваемые

формулами (20.6).

Это и есть искомое интегральное уравнение для неизвестного

давления N0(t).

К этому уравнению

следует, как и прежде, при­

соединить

для определения

дуги контакта

соотношение

(19.24):

 

 

 

$N0(t)dt

= - P ,

 

(20.13)

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

 

 

 

 

 

N0

(**) = N0

( - "О =

0,

 

t4 = re'N

 

(20.14)

где

t^w

 

 

— аффиксы

концевых

точек линии ч, лежащие от

оси у справа

и слева соответственно.

 

 

 

 

 

 

 

 

Подстановка

(19.36)

при p — г

приводит эти уравнения к виду

 

m

м

,е,

,

J _

ç

Nfàdn

 

,

m ß 2 |

ç

N [n)di\

_ Ea +

ß 2 p

( p .

g « _ ß » i V

W

't

2n

J

 

 

"t" я

( | * - ß * ) J_ uf +

ß»

~~ S 2

- ß a

™ '

 

 

 

 

 

 

_

1

 

 

 

_

1

 

 

 

(20.15)

 

 

 

 

 

 

 

 

§ f w » N {

l ) d

b

= 1'

 

 

 

( 2 0 Л 6 )

где

 

 

 

 

 

—î

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

(g) =

±{h(t0)-Ke[t0(q0Fo(to)-m0F1(W),

 

 

 

 

 

 

\ ( 2

Q 1 ? )

*

=

m

=

-2»-,

N{l)

=

N0(t0),

 

N{1)

= N

(-1)^-0.

J

 

 

 

Po

 

 

Po

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, для определения неизвестных N и ß мы име­

ем сингулярное интегральное уравнение (20.15) при дополни­ тельном условии (20.16).

После определения N и ß функции напряжения Фо{г) и <Di(z) будут даны формулами (20.4) и (20.5). Что касается функций 4f0(z) и x Fi(z), то они, как и в предыдущем парагра-

130

фе, непосредственно находятся по функциям

Ф0 и

Ч / 0 из

гра­

ничных условий задачи.

 

 

 

 

 

 

 

 

Сингулярное уравнение (20.15) относится к типу уравнений,

рассмотренных в § 131 ) , и согласно условию

N(l)

=N(—1)

= ;

= 0, к нему применима

схема

численного решения, указанная

там же, в п. 3°.

 

 

 

 

 

 

 

 

 

Уравнение

(20.15) запишем в виде (13.22). Имеем

 

а Ш N(I) +

j Ш * ! +

J *( Е , „ ) Д/(•)•))Д і

=

/

(20.18)

где

- 1

 

 

 

- 1

 

 

 

 

 

 

 

 

 

 

2mß2

 

 

 

 

V

Ь{1)=\,

 

Ä(g,i]) =

 

 

 

 

 

( ё 2 - ß2 ) Of +

ß 2 ) '

 

Po

/1

г

,

P

+

 

 

 

(20.19)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ я («о +

 

mßj>£ (R-r)

(ga + ß2 )

 

 

 

 

 

1) ( I 2 -

ß 2 ) [І 2

(« - 2r)» +

ß*Ä«J-

 

Согласно изложенному

в п. 3 § 13^З узламии

интерполяции в

нашем случае будут

точки

 

/ля

 

 

 

 

 

 

costfm ,

flm

(m = l ,

...,n),

 

 

 

 

 

 

 

л + 1

 

 

 

 

 

а система линейных уравнений для дискретных значений иско­

мой функции

(система

(13.10))

будет иметь вид

 

 

 

 

 

fm

+

2J ат,ѵЛ^ѵ

(m = 1, . . . , /г),

 

(20.20)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n + l

cos •&

— cos

 

 

-j-k(C0S®m,

C O S # v

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

*

 

 

 

 

 

 

 

 

 

 

6m.v —

f 1 при \m — v\ — 1,3, . . .

 

 

 

 

[0

при \m—v\

0, 2, . . .

 

 

 

Заметим,

наконец,

что по условиям задачи функция Л^(|)

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

в системе

(20.20)

сократит­

ся вдвое. Она будет содержать всего

[ — ]

неизвестных.

 

 

 

§

21. Числовой

пример

 

 

 

 

 

Будем для примера

считать, что круглая

шайба

изготовлена

из

стали, а

окружающее тело — из

чугуна.

В

соответствии

с

этим для упругих

постоянных

примем

следующие

значения

 

ѵо—0,33, .

ѵі=0,27,•

хо=3—4ѵо=1,67,

 

 

хі=3 — 4 ѵ і = 1,92,'

Но

f—^ = 8,1 - 105

кГ/см*,

н

= 2 ( 1

^ V i

)

= 4,5 • Ю5

кГІсмК

2(1Е п+

 

 

 

 

 

 

 

 

 

 

 

 

 

') См. замечание в конце параграфа.

131

Возьмем,

далее, Р=\ кГ

и г—3 см. Вычисления мы про­

ведем в

трех

случаях,

когда

дробь

'7R равна,

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

0,9; 0,5 и 0.

 

 

 

 

п=8,

 

 

 

В системе

(20.20)

возьмем

т. е. будем решать систе­

му из

четырех

уравнений

относительно Nu

N2, N3 и

Nit

Ov = - ^ - ( v =

l , 2 , 3,4).

 

 

 

 

 

 

Результаты

вычислений

приведем в виде

следующей

таб­

лицы:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а 1

 

 

 

 

 

0,9

 

 

0,5

0

 

 

Ni

 

—1,00402

 

—2,24803

—3,18046

 

 

N2

 

—1,89353

 

—4,23317

—5,98690

 

 

N3

 

—2,55914

 

—5,71031

—8,07370

 

 

Nt

 

—2,91538

 

—6,49854

-9,18564

 

 

 

 

—85°52'20"

—88°09'20"

—88°41'40"

 

Для найденных Nv и ß') левая часть выражения (20.16), обозначаемая через c(ß), будучи вычислена с помощью квадра­ турной формулы типа Гаусса:

f M _ 2rnß ^ sin ö y ( c o s a fly— ß°) ^Ѵу

v=l

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

c(ß) = 1,00143 при /7Я = 0,9;

с(ß) =0,99996 при г/Я = 0,5;

c(ß) =0,9999 7 при r/R = 0.

Отсюда видно, что равенство (20.16) выполняется во всех трех случаях с высокой точностью.

Как было отмечено в начале параграфа, рассматриваемая задача при известных допущениях может быть решена в замк­ нутой форме. Давление Л7 (я) и участок контакта f, представля­ ющий собой на сей раз отрезок [—I, I] оси х, даются при этом формулами (Н. И. Мусхелишвили [1J, стр. 438)2 )

- М ( г \ - У ' 2 - * 2 Г ^(0 d t

Г «) dt _ к р

Напомним, что ß определяется через по формуле (19.36). Формулы мы приводим в симметричном случае.

132

причем

'

Но в нашем случае, согласно формуле (20.3а)

*w=4(-f-4-)*.

ввиду чего предыдущие интегралы вычисляются элементарно. Будем иметь

Для сравнения результатов приведем численные замечания напряжений и полярных углов, полученные из предыдущих формул при тех же параметрах задачи.

 

 

 

Т а б л и ц а 2

 

0.9

0.5

0

Ni

—1,00729

—2,25223

—3,18532

 

—1,89324

—4,23315

—5,98690

ІѴз

—2,55063

—5,70303

—8,06575

*. .

—2,90053

—6,48538

—9,17223

—85°52'19"

—88°09'14"

—88°4П4"

Разница между данными таблиц 1 и 2 ничтожно мала.

§ 22. Об обобщении метода

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

<»(£)

=

+

(22.1)

где R , m и п •— некоторые

постоянные. Подходящим

подбором

этих постоянных можно при помощи формулы (22.1)

рассматри­

вать отверстия в форме круга, эллипса, криволинейного четырех­ угольника и др.

При решении задачи о жестком включении .в этих случаях следует, предварительно воспользовавшись отображением (22.1),

133

применить способ, изложенный в § 14, что приведет к сингуляр­ ному уравнению вида (14.12), правда, несколько более сложно­ м у 1 ) . К решению этого уравнения может быть с успехом приме­ нен метод Мультоппа, или же можно его решить каким-либо другим приближенным способом. Этим путем задача о жестком эллиптическом включении решена в работах Уилсона-младшего (Wilson, Jr. [1] и Гори (Goree [1]) .

') Сведение задачи к нитегро-днффереициальному уравнению можно осу­ ществить и несколько иначе,— применением метода И. Н. Карцпвадзе (см. примечание в конце § 14).

Глава четвертая

СПЕЦИАЛЬНЫЕ МЕТОДЫ РЕШЕНИЯ ПЛОСКИХ ЗАДАЧ

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

I. РЕШЕНИЕ ОСНОВНЫХ ЗАДАЧ

ДЛЯ ОДНОГО КЛАССА ОДИОСВЯЗНЫХ ОБЛАСТЕЙ

Класс областей, который' мы будем рассматривать в настоя­ щем отделе, определяется условием, чтобы функция и ( ^ ) , даю­ щая конформное отображение области 5 на полукруг (на полуплоскость с вырезом в форме полукруга), допускала анали­ тическое продолжение на всю область единичного круга (на бесконечную область вне кругового отверстия). После продол­ жения соответствующим образом одного из искомых комплекс­ ных потенциалов, плоская задача для области 5 сводится к не­ которой смешанной задаче для определения одной-единственной функции, голоморфной внутри единичного круга (в плоскости с круговым отверстием). Для таких областей мы укажем спо­ соб решения основных плоских задач, позволяющий их сведение к бесконечной системе линейных алгебраических уравнении.

§ 23. Способ решения

Пусть Б + и Б~ обозначают верхний и нижний полукруги еди­ ничного радиуса на плоскости £ (с центром в начале координат), Ч+ и.7~ — дуговые части их границ соответственно,.а / — отрезок (—1,1) на той же плоскости.

135

После отображения физической области 5 на нижний полу­ круг £~ граничное условие плоской задачи представится в виде (гл. I , формула (1.14)).

ср (а) +

ю(ст)

( о г ) + яр (а) = / (а) на

у - + I,

(23.1)

 

«' (о)

 

 

 

где /(а) —заданная функция. От этой функции

будем

требовать,

чтобы ее значения на отрезке / были вещественны и представля­ ли собой граничные значения некоторой функции fo(z), голо­

морфной в 2~. Предполагается,

что условие (23.1) выполняется

всюду на

за

исключением, быть

может,

угловых

точек.

Распространим

определение

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

потенциала ф(£)

на область 2 +

равенством

 

 

 

 

Ф(С) =

- Ш ( £ ) | І І § _ ^ ( Е ) + / 0 ( Е )

при

£ в 2 + ,

(23.2)

где использовано обозначение (Н. И. Мусхелишвили [1], § 76)

F(Ç)=Tffî. (23.3)

Легко убедиться, что функция ф(£), определяемая равенством (23.2), голоморфна в области S+ и согласно (23.1) удовлетво­ ряет условию

Ф + ( о ) = ф - ( 0 )

н а / .

(23.4)

Заменим в формуле (23.2) переменную £ на £ (считая, что £ принадлежит Б + ) н перейдем в обеих частях равенства к сопря­ женным значениям. Будем иметь

^ ( £ ) = - Ф ( £ ) - т ( £ ) , ф ' 1 £ ) + f ° ® при S e s - (23.5)

Предыдущее равенство позволяет определить потенциал яр(£;) че­ рез функцию ф(£), голоморфную согласно условию (23.4) во всем единичном круге.

Граничное условие (23.1) на отрезке I нами уже использовано. Если в оставшейся части этого условия заменить функцию •ф(о) ее значением из формулы (23.5), получим

Ф (а) - ф (а) +

[со (а) - œ (а)] Ф (а) =

/* (а)

на -у- (23.6>

где

 

 

 

 

/Лет) = /

(а) - / о ( а ) ,

Ф(ог)

Ф' (<0

=

^

Таким образом, для определения одной-единственной функ­ ции ф(£), голоморфной в-круге | £ | < 1 , мы имеем условие (23.6) на нижней полуокружности. Условие на верхней полуокружности

136

получится из равенства (23.6), если заменить в нем точку о на а. Будем иметь

Ф (а) - Ф (а) + [со (а) - а (а)] Ф (а) = - f * (о) на 7+ . (23.6а)' Объединив формулы (23.6) и (23.6а) в одно граничное условие,

приходим к следующей смешанной задаче теории

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

функций.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Найти

голоморфную

в

единичном

круге

 

функцию

ф(£) по

граничному

 

 

условию

 

 

 

 

 

 

 

 

 

 

 

 

Ф ( а ) - ф ( а ) + (а) -Cû(ô)]G(a) =g (а) на

 

у

(у=у~+

у+), (23.7)

 

r

. .

) =

f Ф(ОТ на

у - ,

, .

//* (о)

на

у-,

 

(23.8)v

 

G

(ѵ °

'

7 К / (

J-

s ^w ( а

) =

l

f

 

г \

 

4-

7+,

 

 

 

 

[ Ф(а)

на

7+,

 

 

 

— / *

(ст) на

 

где

(0) —заданная функция на 4-,

а ш(о) —граничные

значе­

ния заданной функции со(£), голоморфной внутри f-

 

 

К

граничной

задаче

(23.7)

при известных

условиях

можно

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

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

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

§ 24. Основная смешанная задача для полукруга1 )

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

Радиус полукруга примем равным единице и расположим упругое тело со штампом в плоскости z=x-\-iy так, чтобы уп­ ругая среда занимала нижнюю половину круга с центром в на­ чале координат. В рассматриваемом случае, очевидно,

z = c a ( £ ) = £ .

Верхний и нижний полукруги в плоскости z будем обозначать соответственно через S+ и S~, а для их дуговых границ, а также для отрезка — 1 < г < 1 , оставим прежние обозначения.

') Настоящий параграф воспроизводит статью автора [10] почти без из­ менения.

137

Предположим, что на линии соприкасания двух тел (уп­ ругого полукруга и жесткого основания) отсутствует скольже­ ние и отставание1 ). Тогда для комплексных потенциалов qp(z) и ip(z), голоморфных в области S~, будем иметь граничные условия

 

Ф (0 + t^Jt) + Ш = f (0

на у-,

(24.1 )

 

х Ф

(t) — Up' (t) — ар (t) = О

на I.

(24.2)

Здесь f(t)

—заданная функция

точки

на нижней полуокружно­

сти, непрерывная на закрытой

дуге,

а х — постоянная, равная

3—4ѵ в

случае

плоской деформации

и ^ ~ ^

в случае

плоского напряженного состояния. Граничные условия (24.1), (24.2) выполняются всюду на соответствующих разомкнутых кривых, за исключением, быть может, концевых точек z = ± l .

Следуя приему, указанному в предыдущем параграфе, опре­ делим в области 5 + функцию q>(z) соотношением

хф (z) = z ф' (z) 4- яр" (z) при z в S+.

(24.3)

Если в формуле (24.3) заменить z на z (считая, что z при­ надлежит S~) и перейти к сопряженным значениям в обеих ча­ стях равенства, то из этого равенства определится ip(z) через функцию cp(z), голоморфную в обоих полукругах:

ір (z) = хф (z) — гф' (г) при г в 5 _ .

(24.4)

Легко убедиться, что граничное условие (24.2) не что иное, как условие аналитической продолжимости функции cp(z) через отрезок /. Для этой функции ф, голоморфной в силу сказанного выше, во всем единичном круге, мы будем иметь на основании формул (24.1) и (24.4) следующее условие:

Ф (0 4- щ

(f) + (t-t)

ф' (t) =

f (t)

на y-.

(24.5)

Продолжив, как и раньше, равенство

(24.5) на верхнюю полу­

окружность ч+, мы

получим

соотношения

для определения

функции cp(z) на всей окружности у:

 

 

 

ф (0 + иф (f) +

(t - 1 ) c(77Ö = g(t)+k

(t) Si (t)

(24.6)

g(t)=f(t)

на y-,

g (*)=/(/)

на y+,

(24.7)

Q (t) = (« -

1)[Ф Çt) - Ф (/)] 4- (* - 7)[фЧО + ф' (01, (24.8)

 

1

на

 

 

 

') Предлагаемый способ принципиально пригоден

и в других

случаях, ес­

ли только коэффициент трения сохраняет постоянное значение вдоль контакт­ ной линии.

138

К решению смешанной задачи (24.6) будем применять метод степенных рядов. Положим в единичном круге

СО с о

 

 

 

Ф (Z ) =

2

akz\

 

Ф' (z) =

У> kctkzk-K

 

(24.9)

Предполагая

 

k=0

 

 

 

 

ft=l

 

 

разложений

в ря­

законность

соответствующих

ды Фурье, положим еще

 

 

 

 

 

 

 

 

 

 

 

 

 

g(t)=

 

2

Akt\

Q ( 0 =

2

Akt"

 

(t = e*). (24.10)

 

 

 

fe=—со

 

 

 

 

k=—CO

 

 

 

 

 

 

 

Величины A h известны. Коэффициенты

же Фурье A h неизве­

стной функции Q(t) представляются согласно (24.8)

формулами

Л л = ^ . [ { ( Х - І ) [ Ф ( о - Ф ( 0 ] +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

(*-*)[ф'(0

+ V'(t)}}t-n-ldt.

 

(24.11)

Характер

искомого

решения

в замкнутой

области

должен

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

энергии.

Для

выполнения

этого

условия

достаточно

 

предположить, что

поле

упругих

смещений

 

непрерывно

в замкнутой

 

области,

а

ком­

поненты

напряоісенил

 

также

непрерывны

 

вплоть

до

 

границы

всюду,

за

исключением

 

точек z = ± l ;

в окрестности

этих

точек

напряжения

 

могут

допустить

 

сингулярность

 

 

порядка

ниже

еди­

ницы

(см., например,

Н. И. Мусхелишвнли

 

[2], §§ 113—115).

При

этих

условиях

можно

утверждать, что производные

ф'(г),

ф"(г)

непрерывны в любой части замкнутой

области, не содер­

жащей точек z = ±

l ,

а вблизи указанных

точек

справедлива

оценка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

І Ф ' ( 2 ) | < С | г ±

 

 

 

( 0 < а <

1, C = const).

(24.12)

Вернемся

к

решению задачи

(24.6).

Составим

выражение

Q{t)

на основании разложений

(24.9). Имеем

 

 

 

 

 

 

 

 

Й (0 =

2 и* [a;

 

- *-*) +

аг

(t -

t~%

(24.13)

 

 

 

 

 

k=i

 

 

 

 

 

 

 

 

 

 

 

 

 

fife [a;

x] = — [(x — 1) ak + (k - f 2) a k + 2

На основании формулы

(24.11) имеем, далее,

 

 

Л 0

= 2

afeufife[û; «],

1

 

fe=i

 

 

œ

 

Л я = ^(sgn/iK^ + Q|„, [a; x]} J - Sßfti fife [a- x],

Л„ = j(sgn n) Q,„, [A; x] - f aUl ax

+

 

 

с о

 

— kâk].

(24.14)

n ~ ± 1, (24.15)

+ 2'a/m fife [я; x], n = ± 2, ± 3, . . . (24.16)

k=i

139

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