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

Sapunov_Prikladnaya_teoriya_uprugosti_2008

.pdf
Скачиваний:
83
Добавлен:
05.06.2015
Размер:
1.44 Mб
Скачать

Отметим, что в разделе 1.4.1 граничное условие в форме Колосова представлено с использованием угла α , который нормаль n составляет с осью x , а здесь перешли к углу θ. Однако не надо упускать из виду, что направления осей r и θ противоположны направлениям n и t . Соответственно, поскольку α = θ± π, то

e2iα = e2iθ , и вид граничного условия остается прежним. Граничное условие в форме Колосова предполагает использова-

ние комплексных функций

 

 

 

 

(z) .

 

Φ(z)= ϕ (z) ,

 

 

Ψ(z)= ψ

 

В свою очередь функции ϕ(z) и

ψ(z) для бесконечной области

(плоскости) определены соотношениями

 

 

 

 

ϕ(z) = −

X +iY

 

ln z z

 

(z)

,

2π(1+ κ)

0

 

 

 

 

 

 

 

 

κ(X iY )

 

 

 

 

 

 

ψ(z) =

 

 

 

ln z z 0(z) ,

2π(1+ κ)

 

 

 

 

 

 

 

 

где Γ = B + iC и Γ′ = B′+ iC

комплексные постоянные (в даль-

нейшем будем принимать C = 0 ); ϕ0 (z) и ψ0 (z) функции, го-

ломорфные вне окружности

L , представляемые в области S раз-

ложениями вида1

 

 

a1

 

a 2

 

 

 

 

 

ϕ0 (z)= a0

+

+

+ . . . ,

ψ0 (z)= a0

+

a1

+

a2

+ ...

z

z2

z

z2

 

 

 

 

 

 

 

 

Дифференцируем функции ϕ(z) и ψ(z). Будем иметь:

1Для конечной многосвязной области функции ϕ0 (z) и ψ0 (z) голоморфны вне любой окружности, охватывающей все контуры (при достаточно больших значениях z ). Если имеется только один контур (плоскость с отверстием), то эти

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

41

 

X +iY 1

 

 

a1

 

 

 

2a2

 

 

Φ(z)= ϕ (z) = −

 

 

 

 

 

 

+Γ −+

 

 

 

 

 

 

...

2π(1+ κ)

z

z2

 

z3

 

Ψ(z)= ψ′(z) =

κ(X iY )1

+Γ′−

a1

2a2

...

 

2π(1+ κ)

 

z

z2

 

 

z3

 

Перепишем полученные соотношения в форме рядов

Φ(z)= ak zk ,

Ψ(z)= ak zk ,

k =0

k =0

где введены новые обозначения для коэффициентов:

 

a0 = Γ = B ,

a0 = Γ′ = B′+ iC,

 

a1 = −

X +iY

 

,

a1 =

κ(X iY )

и т.д.

 

 

 

 

2π(1+ κ)

2π(1+ κ)

 

 

 

 

Отметим, что коэффициенты разложения a1 и a1связаны между собой условием однозначности перемещений κa1 + a1′ = 0 .

Поскольку соотношения, полученные для комплексных функций Φ(z) и Ψ(z), предназначены для использования в граничном условии при r = R , перепишем их в следующем виде:

Φ(z)= ak rk eikθ ,

Ψ(z)= ak rk eikθ .

k =0

k =0

Правую часть граничного условия также представим в форме ряда

N iT = +Ak eikθ ,

−∞

где коэффициенты Ak разложения известны.

42

Подставляя преобразованные соотношения в граничное условие, получим:

a

 

(1+ k)

 

 

ak

 

 

 

k

 

ak +2

 

 

a1

 

 

 

 

 

 

eikθ +

 

eikθ a0e 2iθ

 

eiθ =

 

 

R k

R k|+2

Rk

R

k =0

 

 

 

 

k =0

 

 

=+Akeikθ . −∞

Сравнивая коэффициенты при различных степенях eiθ , определим неизвестные коэффициенты разложений ak и ak через из-

вестные Ak . Будем иметь:

сравнение постоянных членов

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

(2Γ − A0 ) ,

 

 

 

a2

 

 

 

 

 

 

 

 

a2

 

 

 

a0

 

 

 

+ a0 = A0

 

2a0

 

= A0

 

 

a2 = R

 

 

 

R2

R2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

поскольку a0 = Γ = B действительное число;

 

 

 

 

 

 

 

сравнение коэффициентов при eiθ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a1

a1

= A1 a1 a1′ = A1R ;

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

сравнение коэффициентов при e 2iθ

 

 

 

 

 

 

 

 

 

 

a2

 

a0′ = A2

a2 = R2 (A2 + Γ′) a2 = R2 (

 

2 +

 

) ,

 

 

A

Γ

 

R2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

поскольку a0 = Γ′ = B′+ iC;

 

 

 

 

 

 

 

 

 

 

 

 

сравнение коэффициентов при eikθ , когда k 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ak

 

= Ak

ak =

Ak R

k

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Rk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

43

сравнение коэффициентов при e ikθ, когда k 1

 

ak (1+ k )

 

ak+2

 

 

 

 

= ak (1+ k)R

2

 

k +2

 

 

 

 

= Ak

 

ak

+2

 

Ak R

 

 

Rk

 

Rk +2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или, с заменой индекса k на (k 2),

 

 

 

 

 

 

 

 

 

 

 

 

2

Ak +2R

k

при k 3 .

 

 

 

ak = ak 2 (k 1)R

 

 

 

Для определения коэффициентов a1 и a1имеем два уравнения:

κa1 + a1′ = 0 , a1 a1′ = A1R .

Решение этих уравнений относительно упомянутых неизвестных позволяет получить:

 

1

 

 

 

 

a1

= −

κ

A1R .

a1 =

A1R

,

1+ κ

1+ κ

Окончательно можем выписать все определяемые коэффициенты рядов:

a0 = Γ , a0 = Γ′ ,

 

 

a1 =

1

 

 

 

 

1R ,

a1′ = −

κ

A1R ,

 

 

 

A

 

 

 

1+ κ

1+ κ

 

 

 

a2 = R2 (

 

2 +

 

) ,

a2′ = R2 (2Γ− A0 ) ,

 

 

 

A

Γ

 

ak = Ak R

k

( k 3 ) ,

 

 

 

 

(k 1)R

2

Ak +2R

k

( k 3 ) .

 

 

 

 

ak = ak 2

 

 

 

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

44

Задачи

1.1. Одностороннее растяжение пластины, ослабленной круговым отверстием, имеющим радиус R .

В соответствии с условием задачи имеем, что контур кругового отверстия свободен от нагрузки, а на бесконечности напряжения принимают значения:

σx = p ,

σy = 0 , τxy = 0 .

 

 

Действительные постоянные B ,

Bи C, входящие в соотношения для функ-

ций Φ(z) и Ψ(z) посредством комплексных постоянных Γ и Γ

, имеют простой

 

физический смысл (см. раздел 1.3.3). Рассматривая поведение напряжений на бесконечности, имеем:

σx = 2B B

2B B′ = p

B = p / 4 ,

σy = 2B + B

2B + B′ = 0

B′ = −p / 2 ,

τxy = C

C′ = 0

C′ = 0 .

Соответственно, постоянные

 

Γ и Γ принимают значения:

Γ = p / 4 ,

Γ′ = −p / 2 .

Поскольку контур кругового отверстия не нагружен, имеем N iT = 0 на L и, как следствие, все коэффициенты разложения Ak равны нулю.

Переходим к определению коэффициентов разложения комплексных функций Φ (z) и Ψ (z). Будем иметь:

a0 = p / 4 ,

a0 = −p / 2 ,

a1 = 0 ,

a1 = 0 ,

a2 = −pR2 / 2 ,

a2 = pR2 / 2 ,

a 3 = 0 ,

a3 = 0 ,

a4 = 0 ,

a4 = −3pR4 / 2 .

Все остальные коэффициенты равны нулю.

 

45

С учетом коэффициентов разложений комплексные функции Φ (z) и Ψ (z) принимают вид:

 

Φ(z)=

 

p

 

 

 

 

2R

2

 

Φ(z)=

p

 

 

 

 

2R

2

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

 

 

e 2iθ

 

,

 

 

 

 

 

 

4

 

 

 

 

 

z

2

 

 

 

 

4

 

 

 

 

 

r

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ψ(z)= −

p

 

 

 

R

2

 

 

 

3R

4

 

 

 

Ψ(z)= −

p

 

 

 

 

 

R

2

 

 

 

 

 

 

3R

4

 

 

 

1

 

 

 

 

+

 

 

 

 

 

1

 

 

 

e

2iθ +

 

e

4iθ .

 

 

 

 

2

 

 

 

4

 

 

 

2

 

 

 

4

 

 

 

 

 

z

 

 

 

 

z

 

 

 

2

 

 

 

 

 

r

 

 

 

 

 

 

r

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Найденные соотношения для функций Φ (z) и Ψ (z) определяют решение по-

ставленной задачи в общем виде.

Определим теперь соответствующие компоненты напряжений, используя формулы, полученные в разделе 1.2.5:

σr + σθ = 2 [Φ(z)(z)] ,

σθ −σr + 2 iτr θ = 2e2iθ [z Φ′(z)+ Ψ(z)] .

Подставив соотношения для функций Φ (z) и Ψ (z) в приведенные формулы и проведя некоторые простые преобразования, будем иметь:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2R

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σr + σθ =

p

1

 

 

cos 2θ

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σθ −σr + 2 i τrθ =

 

 

 

 

 

 

 

 

pR

2

 

2 pR

2

 

3 pR

4

 

 

 

 

 

 

 

2 pR

2

 

3 pR

4

 

=

 

+

 

 

p cos 2θ+ i

 

+

 

p sin 2θ .

 

2

 

2

 

 

4

 

 

2

 

 

4

 

 

r

 

r

 

 

r

 

 

 

 

 

 

 

 

r

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку второе уравнение распадается на два, для определения напряжений получаем следующие три уравнения:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2R

2

 

 

 

 

 

 

 

σr + σθ =

p

1

 

 

cos 2θ

,

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pR

2

 

 

 

2 pR

2

 

 

 

3 pR

4

 

 

 

σθ −σr

=

 

 

+

 

 

 

 

p

cos 2

θ ,

 

2

 

 

 

2

 

 

 

 

 

 

4

 

 

 

r

 

 

 

 

r

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

2R

2

 

 

3R

4

 

 

 

τrθ = −

 

 

1

+

 

 

 

 

 

 

sin 2θ .

 

 

 

 

 

 

2

 

 

 

 

4

 

 

 

 

 

 

2

 

 

 

 

 

 

r

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

46

Касательное напряжение τrθ определяется третьим уравнением, а для оты-

скания нормальных напряжений используем первое и второе. Окончательно будем иметь:

 

 

 

p

 

 

 

R

2

 

 

 

 

 

p

 

 

 

 

 

4R

2

 

 

 

3R

4

 

 

σr

=

 

 

1

 

 

 

 

+

1

 

 

 

+

 

cos 2

θ ,

 

 

 

 

2

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

r

 

 

 

 

 

2

 

 

 

 

 

 

r

2

 

 

 

 

 

r

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

R

2

 

 

 

 

p

 

 

 

 

3R

4

 

 

 

 

σθ =

 

1+

 

 

 

 

 

1

+

 

 

cos 2θ ,

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

r

 

2

 

 

 

2

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

2R

2

 

 

 

 

3R

4

 

 

 

 

 

 

 

 

τrθ = −

 

 

1+

 

 

 

 

 

 

 

sin 2θ .

 

 

 

 

 

 

 

2

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

r

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение рассматриваемой задачи в действительных переменных приведено в книге [1]. Там же проведен подробный анализ напряженного состояния с обсуждением его особенностей.

Для определения перемещений vr и vθ точек пластины воспользуемся их комплексным представлением в полярной системе координат (см. раздел 1.2.5):

2μ(vr +ivθ) = eiθ [κϕ(z)zϕ′(z)− ψ(z)] .

Необходимые комплексные функции ϕ(z) и ψ(z) определим по известным функциям Φ (z) и Ψ (z):

ϕ(z)= Φ(z)dz =

p

 

2R

2

 

 

 

 

 

 

 

 

 

ϕ(z)=

 

p

 

2R

2

 

 

 

z +

 

 

 

 

 

 

 

 

reiθ +

 

e

iθ

,

 

4

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ(z)= Ψ (z)dz = −

p

 

 

 

 

R

2

 

R

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z +

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

z

z

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ(z)= −

p

 

 

 

 

R

2

 

 

 

 

R

4

 

 

 

 

 

 

 

 

 

reiθ +

 

 

e

iθ

 

 

 

e

3iθ .

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

r

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя значения комплексных функций в уравнение для перемещений и отделяя действительную и мнимую части, получим:

47

 

p

 

 

 

 

 

R

4

 

 

 

 

 

(κ−1)+ 2R2 +2

 

(κ+1)+ r2

 

 

 

 

 

vr =

 

r2

R2

 

 

cos2θ

,

8μr

r2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vθ = −

p

R2

(κ−1)+ r 2

+

R4

sin 2θ .

 

 

2

 

4μr

 

 

r

 

 

 

 

 

 

 

 

1.2. Всестороннее растяжение пластины, ослабленной круговым отверстием, имеющим радиус R .

Решение задачи о всестороннем растяжении пластины, ослабленной круговым отверстием, строится по схеме, предложенной при решении задачи 1.1.

В соответствии с условием задачи на бесконечности имеем

σx = σy = p , τxy = 0 .

Определение действительных постоянных B , Bи C, входящих в соотношения для функций Φ (z) и Ψ (z), в рассматриваемом случае дает:

σx = 2B B

2B B′ = p

B = p / 2 ,

σy = 2B + B

2B + B′ = p

B′ = 0 ,

τxy = C

C′ = 0

C′ = 0 .

Соответственно, для постоянных Γ и Γ′ получим:

Γ = p / 2 , Γ′ = 0 .

Как и в предыдущей задаче, имеем N iT = 0 на L и Ak = 0 . Переходя к определению коэффициентов разложения комплексных функций Φ (z) и Ψ (z), будем иметь, что только два коэффициента отличны от нуля:

 

 

a0 = p / 2 , a2 = pR2 .

С

учетом найденных

коэффициентов разложений комплексные функции

Φ (z),

Ψ (z) и ϕ(z), ψ(z)

принимают вид:

48

Φ(z)=

p

,

Ψ (z)=

pR 2

;

2

z 2

 

 

 

 

 

 

ϕ(z)=

p

z ,

ψ(z)= −

pR 2

.

 

2

 

 

 

 

z

 

 

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

 

 

 

 

R

2

 

 

 

 

 

R

2

 

 

 

σr

= p

1

 

 

 

,

σθ = p

 

1 +

 

 

,

τrθ = 0

;

 

 

 

2

 

 

 

 

 

 

r

 

 

 

 

 

r

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vr =

 

 

p

 

[r 2 (κ−1)+ 2R2 ] , vθ = 0 .

 

 

 

4μr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

гости справедливы принцип пропорциональности решения заданной нагрузке и принцип независимости действия сил (принцип суперпозиции).

1.3. Равномерное нормальное давление, приложенное к обводу кругового отверстия радиуса R в пластине.

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

Соответственно, постоянные Γ и Γ′

равны нулю ( Γ = Γ′ = 0 a0 = a0 = 0 ),

а проекции N и T внешней нагрузки,

действующей на окружности L контура,

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

N = −p , T = 0 .

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

N iT = +Akeikθ .

−∞

49

В рассматриваемом случае, когда N iT = −p , будем иметь:

A0 = −p , Ak = 0 при k 0 .

Переходим к определению коэффициентов разложения комплексных функций Φ (z) и Ψ (z). В соответствии с имеющимися формулами, получим

a2 = pR2 ,

а все остальные коэффициенты разложений функций Φ (z) и Ψ (z) равны нулю.

С

учетом найденных

коэффициентов разложений, комплексные функции

Φ (z),

Ψ (z) и ϕ(z), ψ(z)

принимают вид:

 

 

 

Φ(z)= 0 ,

Ψ(z)=

 

pR 2

;

 

 

z 2

 

 

 

 

 

 

 

 

ϕ(z)= 0 ,

ψ(z)= −

pR 2

.

 

 

 

 

 

 

 

z

 

 

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

σr

= −

pR 2

, σθ =

pR

2

,

τrθ = 0

;

r 2

 

 

r 2

 

 

 

 

 

 

 

 

 

 

 

 

v

r

=

pR2

,

v

 

= 0 .

 

 

 

 

 

 

 

 

 

 

 

 

2μr

 

θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.4. Сосредоточенная сила, приложенная в точке неограниченной плоскости.

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

Поскольку в поставленной задаче на бесконечности напряжения равны нулю, постоянные Γ и Γ′ также равны нулю ( Γ = Γ′ = 0 a0 = a0 = 0 ).

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

50

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]