Сосредоточенные силы и моменты в упругой полуплоскости с отверстием Текст научной статьи по специальности «Физика»

ВЕСТНИК САНКТ-ПЕТЕРБУРГСКОГО УНИВЕРСИТЕТА

Сер. 10. 2009. Вып. 2

УДК 539.3 Ю. Г. Пронина

СОСРЕДОТОЧЕННЫЕ СИЛЫ И МОМЕНТЫ В УПРУГОЙ ПОЛУПЛОСКОСТИ С ОТВЕРСТИЕМ *)

1. Введение. Аналитические решения задач линейной теории упругости о сосредоточенных силах и моментах получены лишь для немногих областей с несложной границей, таких как бесконечная плоскость, полуплоскость, круг, клин, две соединенных полуплоскости [1, 2]. В случае неодносвязной области задача существенно усложняется. Теоретические и экспериментальные исследования напряженного состояния тел с отверстиями под действием нагрузки, приложенной на их границе, проводились многими учеными. Результаты и методы этих исследований представлены, например, в работах [3-7]. В общем случае задача приводится к решению граничных интегральных уравнений (ГИУ) [1, 2, 6-8]. В отдельных ситуациях решения построены с помощью конформного отображения или криволинейных координат [3, 9, 10]. В частности, в монографии [3] указаны вытекающие из формул G. Jeffery [9] решения некоторых авторов о сосредоточенных силах на границе полубесконечной пластины с круговым вырезом. В ряде задач используется аналитическое продолжение потенциалов Г. В. Колосова в верхнюю полуплоскость или их представление в виде специальных рядов [10, 11].

В отличие от упомянутых работ в данной статье изучено тело с отверстием, загруженное не только усилиями, приложенными на его сложной границе и на бесконечности, но также силами и моментами, действующими во внутренних точках. Исключение составляет публикация [10] о сосредоточенной силе в плоскости с двумя равными круговыми вырезами. Выведены интегральные уравнения Фредгольма первого рода, аналогичные уравнениям, полученным М. П. Савруком без учета сосредоточенных воздействий [6]. Помимо этого, задача непосредственно приведена к системе линейных алгебраических уравнений с помощью введения дискретной фиктивной нагрузки.

2. Постановка задачи. Рассмотрим упругую полуплоскость у ^ 0 с отверстием, описываемым простым гладким замкнутым контуром Г (рис. 1). Будем предполагать, что контур Г имеет кривизну, удовлетворяющую условию Гельдера [1].

При этом

1) на бесконечности к полуплоскости приложены постоянные напряжения аХХХ = qX° = const, а™ = °, а™ = 0;

2) прямолинейная граница у = 0 полуплоскости загружена нормальными qyy (x) и касательными qxy(x) усилиями, причем qyy(x) и qxy(x) — функции абсциссы х, удовлетворяющие при больших \x\ условиям [1]

qyy(x) = o(x-1) , qxy (x) = o (x-1) ;

Пронина Юлия Григорьевна — доцент кафедры вычислительных методов механики деформируемого тела факультета прикладной математики-процессов управления Санкт-Петербургского государственного университета. Количество опубликованных работ: 40. Научные направления: сингулярная плоская задача теории упругости, разрушение тел в коррозионных средах. E-mail: DearJuly@gmail.com.

+ ) Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (гранты №06-01-0171, 08-01-00394).

© Ю. Г. Пронина, 2009

Рис. 1. Представление исходной задачи в виде суперпозиции двух вспомогательных задач

3) на контуре Г отверстия действуют известные нагрузки с нормальной qnn (£) и касательной qnt(£) составляющими, непрерывными на Г, £ Е Г;

4) в точках гр, гш полуплоскости с отверстием находятся сосредоточенная сила Р = X + гУ и момент М соответственно.

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

Цель работы: построение численно-аналитического решения поставленной задачи.

Вопросы существования и единственности решения и при указанных условиях освещены в монографиях [1, 7].

3. Вывод основных разрешающих уравнений. Искомое поле напряжений представим в виде [12]

ахх(г) = (z), ауу(х) = (z), аху(г) = а£) (г)+ аХ) (z), (1)

где слагаемые с индексом (1) суть напряжения в полуплоскости с прямолинейной границей, загруженной усилиями qyy(x), qxy(x), qX° и сосредоточенными воздействиями Р и М (задача 1); индексом (2) отмечены аналогичные величины в сплошной полуплоскости (без отверстия), внутри которой на линии, соответствующей контуру Г, действуют неизвестные заранее и подлежащие определению усилия р(£) = рх(С) + гру(С), С Е Г (задача 2) — рис. 1; z = х + гу. Очевидно, что если внешняя нагрузка задана только на контуре отверстия, то все компоненты напряжений в задаче 1 равны нулю.

Компоненты напряжений в каждой из этих задач (р = 1, 2) удовлетворяют формулам Г. В. Колосова [1]

^) + аУУ ^ = 4ИеФ(р)^), аУУ ^) - аХ]^) + 21*%)^) = 2^)'^) + Ф(р)^)]. (2)

Черта сверху означает комплексное сопряжение, штрих — производную по аргументу. Функции Ф(1)^) и Ф(1)^), согласно [1, 13], определяются выражениями

ф<1,м =

Чуу(т) - 'Щху(т) ^ _

г — г 4

Х + гУ

27г(1 + х)

+

X -гУ

гМ

2п(1 + к) ^ — гр)2 2п(z — хт)2 '

qy у(т ) + кху(т)

в,т +

(3)

+ _±_ [ Чуу(т) - кху(т) тЛт_Ч^__ Х + 1У

2пг

(т — ^2

2 2п(1 + к)

+

_(z — zp)2 (z — гр)2 _

+

+

Х-гУ 27г(1 + х)

1,М ~2тт

+

2 Zp(^Zp гр)

(г - гр)2 (г - гр)3

+

+

2г

где интегрирование производится по вещественной оси (т - локальная переменная интегрирования); к = 3 — 4v в случае плоской деформации и к = (3 — V)/(1 + V) при обобщенном плоском напряженном состоянии; V - коэффициент Пуассона материала полуплоскости.

Функции Г. В. Колосова для второй задачи запишем в форме [14]

Ф(2)(^

р(С)

2п(1 + к)

1к

+ ■

—

р(С)

С ~ С 1

—— ав,

27г(1 + ж) (г - С)2

1

к

р

р

z — z

z — z

р

р.

т — z

р

р

1

к

z — z

z — z

р

р

1

1

3

Ъ(2)(г) = - !

Р(0

2п(1 + к)

+

><

[(г-С)2 (г - С)2]

ів +

(4)

+

р(С )

2п(1 + к)

+

с - с 2С(С - С)

г - С г-С {г - С)2 (г - ()3

ів.

Отметим, что р(С) означает величину, комплексно сопряженную к р(С), <1з - дифференциал дуги контура Г, £ Є Г, £ ^ Г.

При таком выборе функций Ф(р) (г) и ^(р)(г) краевые условия удовлетворяются на прямолинейной границе полуплоскости и на бесконечности. На контуре отверстия они имеют вид

яПП^+^ПП (і) = Чии(і), яПА^+^Пл (і) = Чиг(і)

І2)

(!)/•

12),

аПП(і) - іаПі!(і) = дпп (і) - ідпі(і) - аПП)(І) + іаПі)(і), і Є г

(5)

Здесь апП (£) и а(£) - нормальная и касательная компоненты напряжения на кривой Г в задаче 1 (р =1) и задаче 2 (р = 2). Следуя [1], правомерно написать

4п}(-г) - і-аїй (*) = 211еФ^(,г) - \гФ^'{г) + ^р\г)

Лгд

(р =1, 2),

(6)

где в - угол между нормалью к контуру Г и осью 0х.

Используя формулы (3), (4), (6), краевое условие (5) можно привести к интегральному уравнению относительно неизвестной нагрузки р(£):

J р(С)к(г, () (1в + ! р(С)Д(і, С) Лв = /(і), і є г,

(7)

в котором

_ с-с

2тг (1 + я) 1 і - С ' і - С (і - С)2

1к

+

+

і-С — С)

|_(*-С)2 (і - С)2

„2ів

R(t, С) = —

і к с ~ с

=—= + =—г + , +

2тг (1 + к) І І - С і - С (І-СУ

(8)

+

I* - С і-с (і - с)2

(і - С)3

Лів

1 [ Чуу(т) — Щху(т)

т — і

іт —

1 [ Чууіг) + ідху(т) _ с _

2тті І т — і 2

с

1

к

или

1

и

1 Г й»(т)-«/,»(т) (г_т),;т+ 1

2тп У (г — £)2 27Г*

9та(т) + г^(т) <7^

г 2

+

Мг

~2тт

1

1

(г - гт)‘2 (г - гт)2

М'1

~2тг

_1_________2(2 - гт)

(£ — гт)2 (4 — ит)2 (4 — ит)3

+

в21в -

- (X + гУ)К(г, гр) - (X - гУ)Я(г, гр).

При вычислении интегралов в выражениях (4), (7) можно пользоваться следующими соотношениями [14]. Пусть, например, «верхняя» и «нижняя» части контура Г описываются уравнениями типа

С = £ + «7(0, С = £-«7(0, се Г, с^г.

Поскольку

^ = У1 + Ь?£(£)]2 ^

< = [1 + Ц(0] ^ ^С = [1 - «7^(0] ^

то для «верхней» и «нижней» частей контура справедливы выражения

% =

ас

ЙС

лЛ + ^ОР

<1в = -------5-— <1С,

1 + *^(0

лЛ + [^(ОР _

^ = ----------гг^г- с^С-

1 - «?£(0

4. Периодические воздействия. Рассмотрим теперь ситуацию, когда в полуплоскости с отверстием действуют периодические системы сил и моментов. При этом, конечно, следует иметь в виду, что тело должно находиться в равновесии. Бесконечный ряд сосредоточенных сил необходимо уравновесить периодической или равномерной нагрузкой на прямолинейной границе у = 0, усилиями на бесконечности, либо другой системой или системами сил, приложенных во внутренних точках тела. При само-уравновешенной нагрузке условие, накладываемое на функции дуу(х) и дху(х) в п. 2, становится необязательным. Будем считать, что условия равновесия выполняются.

Итак, пусть в полуплоскости с отверстием в точках = гр + к1р, = гт + к1т

(к = 0, ±1, ±2,...) действуют периодические системы одинаковых сосредоточенных сил Р = X+ гУ и моментов М соответственно (рис. 2). Для полноты исследований положим, что нагрузка на прямолинейной границе также является периодической:

Яуу(х) Яуу(Х + к1д ^, Чху(х) Яху(х + к1д \ х , b\, \Ь а\ ^ 1д.

В таком случае, согласно [14, 15], функции Г. В. Колосова для задачи 1 имеют вид

ф(1)(» = -7^- Jкуу(Т) - Щху(т)] — <1т + Ц- -

и

Х + гУ 21р{\ -\- х)

п(г - гр) п(г - гр)

—— + х ——

п(X - гУ) гр - гр

Мпг

21р(1 + к) вт2 [п(г - 2р)/1р\ Ыт вт2 [п(г - Гт)/1т\’

п(Т - г) гп Г дуу (т) - гдху (т)

212г У вт2 [п (т - г)/1я\

с1т-¥г +

+

Х + гУ 2/р(1 -\- х)

п (г — гр) п (г — гр)

^-Ц—^ + х с^-Ц——

гр гр + г

+

к г

1р\вт [п(г - гр)/1р\ вт [п(г - Гр)/1р\

+

+

X -гУ

21 р (1 + к)

к ctg

п(г - гр) п(г - гр)

1р

+ <*£-

+

+

т(гр - гр)

2пг ctg[п(г - Гр)/1р\

вт [п(г - гр)/1р\ 1р вт [п(г - гр)/1р\

Миг

2ттг ^[тт(г - гт)/1т\ 2т)/1т\ [тг(г ^т)/1т\ 1т si.Il [тг(^ ^т)/^т]

+

Комплексные потенциалы Ф(2) (г) и Ф(2) (г) остаются без изменений. В силу этого, ядра К(г,() и Я(г,С) интегрального уравнения (7) определяются прежними выражениями (8), а его правая часть / (г) принимает вид

о

/0) = Чпп{Ь) - гды(г) + —!-т J[qyy(т) - 1дху{т)}

г(т - г)

ат —

о

1

у- ] [%у(т) +кху(т)} (*£■

Г(т - Г)

р

р

1

о

я

я

р

р

п

1

I

р

я

I

я

b _ b i (<,„{ r) *■+f Ч:)г'“ь,(:), dr+%

ln

2/2*

sin [n(r — t)/lq ] 2

e2i0 +

+

+

+

X + iY 2/p(l + я)

X-iY 2/p(l + x)

X + *У

п (t — г„) п (t — z„) п

ctg —-—-—^ + x ctg v

2lp(1 + к) п к (t — t)

ctg

ctg

lp

7Г (t — ZP)

lp

7Г (t — ZP)

+ к ctg + к ctg

lp

7Г (t — ZP)

lp

7Г (t — Zp)

+ -r

lp sin [n(t — Zp)/lp]_ Zp Zp

lp sin [n(t — Zp)/lp]_

+

+

Ip sin [n(t — Zp)/lp] _ ZZp - Zp

e2ie +

X-iY 2/p(l + я)

lp sin2 [n(t — Zp)/lp]

n(t — zp) п (t — Zp)

>c ctg-Ц—21 + Ctg-——^ +

п + —

п2 2 (Zp — Zp)(t — t) ctg[n(t — Zp)/lp]

lp sin [n(t — Zp)/lp]

l2

lp

sin [n(t — Zp)/lp] _

e2if> +

p

p

p

p

+

Miri

Ж

1

1

sin [n(t — Zm )/lm] sin [n(t — Zm)/lm]_

Miri

1

sin [n(t — Zm)/lm]

1

+

2n(t — t) ctg[n(t — Zm )/lm]

„2i6

(10)

sin2 [n(t — Zm)/lm] lm sin2 [n(t — Zm)/lm]_

Следует сказать, что аналогичные формулы нетрудно вывести для сил и моментов, распределенных по некоторым кривым, лежащим внутри тела, или для других точечных воздействий (краевых дислокаций, силовых диполей, центров расширения-сжатия). Комплексные потенциалы Г. В. Колосова для них представлены, например, в работах [2, 14, 16, 17].

Различные методы решения ГИУ описаны в [6, 7].

5. Приведение задачи к системе алгебраических уравнений. Минуя приближенное решение полученных выше ГИУ, задача может быть непосредственно сведена к системе алгебраических уравнений с помощью введения дискретной неизвестной нагрузки [8, 15]. С этой целью в качестве второй вспомогательной задачи возьмем задачу о сплошной полуплоскости, загруженной в точках Zj, лежащих справа от контура Г (при обходе контура по часовой стрелке), фиктивными сосредоточенными силами Pj = P(Zj) = p'j + ipVy (j = 1,...,N) — рис. 3.

Расположение точек Zj в зоне отверстия на удалении от контура Г обусловлено тем, что в окрестности мест приложения сосредоточенных сил компоненты напряжения и смещения перестают быть ограниченными. Решения задач о точечных воздействиях правомерно применять только к частям тела, не слишком близким к местам их приложения [1]. Размещение некоторых точек Zj на самом контуре Г возможно только в том случае, если на контуре отверстия действуют заданные сосредоточенные силы или же контур имеет угловые точки. Чтобы познакомиться с методом решения задач теории упругости для многосвязных областей с угловыми точками, отсылаем читателя к исследованиям [6]. Отметим лишь, что при наличии угловых точек или участков контуров отверстий с достаточно малым радиусом кривизны целесообразно использовать

Рис. 3. Дискретизация неизвестной нагрузки

не только фиктивные сосредоточенные силы, но также и фиктивные точечные воздействия более высокого порядка, такие как двойные силы и центры расширения-сжатия [16, 17].

Итак, в таком случае вместо формул (4) будем иметь соотношения [13]

' [ Рз ( 1 \ к_ \ | Рз 0 - С;

~ н) \г ~ Сз г - Сі ) 27г(1 -ь

'( Рз ( 0 , <3

ф|Я<г) - + т^)+ 2,(1 + *) <2 - с,р

N

Р7 ( * і Сі ~Сі 2Сі(^і~Сі)

+

2тг(1 + лг) ^ - 0 * - С,' (* - О)2 (* - Сі)3 ) ] '

Внося в граничное условие (5) выражения (6), (3), (11) и разделяя затем его на вещественную и мнимую части, после соответствующих преобразований получим систему двух вещественных уравнений

N N

Ч-

£рХ Кг (і,-) + £ рУ Яг (і,-) = /і (і),

З = г -=г

N N

Т,Рз К2 (ЪСз ) + Ц Ру Я (1'Сз ) = /2 (г), (12)

3 = 1 3 = 1

£ е Г.

Здесь

К1 (г, з) = и.е(к(г, з) + я(г, 3)), К2(г, 3) = 1ш(к(г, Сз) + я(г, Сз)), я1 (г,Сз) = 1ш(-к(г,Сз) + я(г,Сз)), Я2(г,Сз) = ^е(к(г,Сз) - я(г,Сз

/1(г) = Я,е/(г) = дПи(г) - стПП(г), /2(г) = 1ш/(г) = -диг(г) + (г),

где ядра K(г,(),Я(г,() определяются формулами (8), а правая часть /(г) - соотношениями (9) или (10).

Заставим уравнения (12) удовлетворяться в N точках е Г (к = 1,..., Ы). Тогда краевые условия сведутся к системе 2Ы линейных алгебраических уравнений относи-

х У

тельно неизвестных рз и рз:

N N

^2р'3 К1(ги, Сз) + ^2 Ру я^к, Сз) = ЛЫ,

з=1 з=1

N N

^2рз К2(гк, Сз) + ^2 Ру я2(гк, Сз) = /2 (г к), (13)

з=1 з=1

гк е Г\ к = ^^.^ Ы.

6. Вычисление напряжений и результаты расчетов. Решив систему уравнений (13) относительно неизвестных рх и рУ, по формулам (11) несложно определить функции Ф(2) и Ф(2). Окончательно комплексные потенциалы Г. В. Колосова для исходной задачи выражаются суммой

Ф(г) = Ф(1) (г) + Ф(2)(г), Ф(г) = Ф(1)(г) + Ф(2)(г).

Зная функции Ф и Ф, по формулам Г. В. Колосова (2), (6) нетрудно вычислить компоненты напряжений в любой точке тела. Или, иначе, можно сначала с помощью формул Г. В. Колосова найти компоненты напряжений в задачах 1 и 2, а затем уже, просуммировав их, согласно условию (1), получить распределение напряжений для исходной задачи.

Отметим, что метод решения системы уравнений (12) (метод коллокаций) не является сходящимся. При значительном увеличении числа фиктивных сил, а следовательно, и размерности системы уравнений (13) ее матрица становится плохо обусловленной. Тип фиктивных особенностей (силы, диполи, ...) и места их расположения следует выбирать из условий минимальной разности между заданными на контуре отверстия усилиями цпп,Чт и вычисляемыми в ходе решения задачи напряжениями опп, оП1 на контуре Г в точках, отличных от узлов коллокации г к, где решение удовлетворяется автоматически. Указанная разность характеризует погрешность метода. Фиктивную нагрузку необходимо подбирать для каждой ситуации индивидуально: она зависит от кривизны контура и усилий на нем. Для тестирования данного метода было исследовано несколько известных задач, при этом получено хорошее соответствие результатов вычислений и имеющихся теоретических и экспериментальных данных [3, 5]. Расчеты показали, что при непрерывной нагрузке на контуре Г оптимальнее выбирать расстояние от точек Сз приложения фиктивных сил до контура Г не менее чем удвоенное расстояние между соседними точками Сз. В таком случае суперпозиция полей напряжений от действия всех фиктивных сил приводит к достаточно гладкому их распределению на контуре и далее внутри тела. При соблюдении этого условия изменение расстояния от точек Сз до контура Г или размерности системы (13) даже в несколько раз (пока обусловленность матрицы системы не становится критической) в отдельных случаях практически не влияет на погрешность расчетов. В протестированных модельных задачах при оптимальном выборе расчетных параметров погрешность не превышала 3%.

На рис. 4 приведен пример результатов расчетов, проведенных по предложенной схеме. На нем для сравнения показано распределение безразмерных напряжений охх/Р,

вызванных действием сосредоточенной силы Р = гУ в полуплоскости с круговым вырезом (а) и без него (б). Здесь Р - модуль силы Р. Графики построены для следующих исходных данных: точка приложения силы Р гр = — 2г; отверстие — круг радиусом 0.3 с центром в точке гс = —0.9г. Распределение напряжений представлено в области {х Е [—2, 2], у Е [—1.9, 0]} с исключенным кругом \г + 0.9г| < 0.3.

■2 -2 -2 -2

Рис. 4. Распределение напряжений в полуплоскости с отверстием (а) и без него (б)

Из рис. 4 видно, что наличие отверстия повышает уровень напряжений. Дальнейшие расчеты показали, что чем ближе вырез к месту приложения силы Р, тем выше значения напряжений на его границе и вблизи нее. Однако стоит сказать, что любая граница, будь то внешняя поверхность тела или контур отверстия, приводит к росту одних компонент напряжений и уменьшению других. Определенное влияние на рост напряжений оказывает кривизна границы. При увеличении радиуса кривизны максимальные напряжения на контуре отверстия стремятся к аналогичным величинам на прямолинейной кромке полуплоскости при равных расстояниях а от точки приложения силы Р до исследуемой границы. Для относительно больших отверстий такое влияние незначительно: при увеличении радиуса кривизны от г = а до г = го напряжения возрастают лишь на 20%. Это, конечно, относится к случаям, когда силы направлены так, как в рассмотренном выше примере. В некоторых других ситуациях (как, например, при растяжении) уменьшение радиуса кривизны отдельных участков контура отверстия приводит к повышению уровня напряжений. Близость выреза к прямолинейной кромке полуплоскости также способствует росту напряжений. При уменьшении расстояния между отверстием и прямолинейной границей от Н = го до Н = г (где г - радиус кругового отверстия) напряжения возрастают примерно на 30%. Дальнейшее сближение указанных границ приводит к более существенному росту напряжений, особенно в перемычке между ними (где в пределе при Н 0 они стремятся к бесконечности). Отметим, что предложенный метод дает хорошие результаты и для близко расположенных границ области. В заключение укажем еще на тот очевидный факт, что чем меньше максимальный линейный размер отверстия, тем меньшее искажение оно вносит в поле напряжений, вызванное сосредоточенными воздействиями.

1. Мусхелишвили Н. И. Некоторые основные задачи математической теории упругости. М.: Наука, 1966. 707 с.

2. Греков М. А. Сингулярная плоская задача теории упругости. СПб.: Изд-во С.-Петерб. ун-та, 2001. 192 с.

3. Савин Г. Н. Распределение напряжений около отверстий. Киев: Наукова думка, 1968. 887 с.

4. Уфлянд Я. С. Интегральные преобразования в задачах теории упругости. Изд. 2-е, доп. Л.: Наука, 1967. 402 с.

5. Петерсон Р. Коэффициенты концентрации напряжений: Графики и формулы для расчета конструктивных элементов на прочность / пер. с англ. И. А. Нечая и др. М.: Мир, 1977. 304 с.

6. Саврук М. П., Осив П. Н., Прокопчук И. В. Численный анализ в плоских задачах теории

трещин. Киев: Наукова думка, 1989. 248 с.

7. Линьков А. М. Комплексный метод граничных интегральных уравнений теории упругости. СПб.: Наука, 1999. 382 с.

8. Пронина Ю. Г., Чистякова Е. О. Первая основная краевая задача для упругой полуплоскости с отверстием // Нелинейные проблемы механики и физики деформируемого твердого тела / под ред. К. Ф. Черныха. СПб.: Соло, 2004. Вып. 8. С. 232—240.

9. Jeffery G. B. Plane stress and plane strain in bipolar coordinates // Philos. Trans. Roy. Soc. London. Ser. A. 1921. Vol. 221. P. 265-293.

10. Sekhri J. N. Some problems relating to stresses in an infinite plate with holes. IV. The effect of

a concentrated force // Proc. Nat. Inst. Sci. India. Ser. A. 1961. Vol. 27, N 4. P. 295-304.

11. Араманович И. Г. О распределении напряжений в упругой полуплоскости, ослабленной подкрепленным круговым отверстием // Докл. АН СССР. 1955. Т. 104, № 3. С. 372-375.

12. Даль Ю. М. Напряженно-деформированное состояние упругой полуплоскости с краевой выемкой // Вопросы механики и процессов управления / Труды ЛГУ им. А. А. Жданова. Л.: Изд-во Ленингр. ун-та, 1977. Вып. 13: Академик В. В. Новожилов — ученый, педагог, гражданин. 1990. С. 77-78.

13. Даль Ю. М., Пронина Ю. Г. Сосредоточенные силы и моменты у границы упругой полуплоскости // Изв. РАН. Сер. МТТ. 1998. № 5. С. 78-87.

14. Веселков С. Ю., Даль Ю. М., Пронина Ю. Г. Сосредоточенные силы и моменты в некоторых двумерных задачах теории упругости. СПб.: Науч.-исслед. ин-т химии С.-Петерб. ун-та, 1998. 90 с.

15. Пронина Ю. Г. Концентрация напряжений в упругой полуплоскости с краевыми выемками // Изв. РАН. Сер. МТТ. 1998. № 1. С. 103-109.

16. Пронина Ю. Г. Центры расширения-сжатия в упругой полуплоскости // Вестн. С.-Петерб. ун-та. Сер. 1: Математика, механика, астрономия. 2007. Вып. 2. С. 140-149.

17. Пронина Ю. Г. Некоторые сингулярные решения теории упругости для полуплоскости // Актуальные проблемы прочности: Материалы XLVI Междунар. конференции. 15-17 окт. 2007 г. Витебск, Беларусь. Ч. 2. Витебск: УО «ВГТУ», 2007. С. 243-248.

Статья рекомендована к печати проф. Ю. М. Далем.

Статья принята к печати 25 декабря 2008 г.