Влияние кругового отверстия на напряженное состояние оболочки произвольной гауссовой кривизны

ВЕСТНИК ПНИПУ

2014 Механика № 1

УДК 539.3

Е.Н. Довбня, Н.А. Крупко

Донецкий национальный университет, Донецк, Украина

ВЛИЯНИЕ КРУГОВОГО ОТВЕРСТИЯ НА НАПРЯЖЕННОЕ СОСТОЯНИЕ ОБОЛОЧКИ ПРОИЗВОЛЬНОЙ ГАУССОВОЙ КРИВИЗНЫ

Работа посвящена определению напряженно-деформированного состояния изотропной оболочки произвольной гауссовой кривизны с круговым отверстием, расположенным в центре конструкции. Оболочка находится под действием осевого растяжения или внутреннего давления. Использовались уравнения теории пологих изотропных оболочек, которые совпадают с уравнениями теории изотропных оболочек с большим показателем изменяемости. Были задействованы интегральное преобразование Фурье, теория обобщенных функций. В результате задача сведена к решению системы граничных интегральных уравнений. Одно из преимуществ использования метода граничных интегральных уравнений для исследования напряженно-деформированного состояния оболочек, ослабленных отверстием, состоит в возможности определять искомые величины непосредственно на контуре отверстия, не вычисляя их на всей поверхности оболочки. Для получения ядер сингулярных интегральных уравнений были использованы интегральные представления перемещений и фундаментальные решения уравнений статики пологих изотропных оболочек. В качестве неизвестных функций использовались комбинации перемещений, углов поворота и их производных. Аналитические выкладки существенно упрощаются, если считать неизвестными на контуре не четыре функции, как это принято, а пять. В данной работе в качестве пятого уравнения используется дифференциальное уравнение, связывающее неизвестные функции. При численном решении задачи для сведения системы интегральных уравнений к системе линейно-алгебраических уравнений использовались специальные квадратурные формулы для интегралов типа Коши, если неизвестные функции имели корневую особенность на концах промежутка интегрирования. Для сведения дифференциального уравнения к линейно-алгебраическому уравнению использовался метод конечных разностей. Приведены результаты значений коэффициентов концентрации напряжений в зависимости от кривизны изотропной оболочки. Также было произведено сравнение результатов с другими исследователями.

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

E.N. Dovbnya, N.A. Krupko

Donetsk National University, Donetsk, Ukraine

INFLUENCE OF CIRCULAR HOLE ON THE SHELL STRESS STATE FOR ARBITRARY GAUSSIAN CURVATURE

The work is devoted to determining isotropic shell of stress-strain state for arbitrary Gaussian curvature with a circular hole, located in the center of the structure. An axial tension or an internal pressure is applied to the surface of the shell. The isotropic shallow shell theory equations were used, which coincide with the isotropic shell theory equations with a large measure of variability. The integral Fourier transformation and the theory of generalized functions were applied. As a result the problem was reduced to solving the system of boundary integral equations. One benefit of using the method of boundary integral equations for the study of shell stress-strain state weakened by a hole is the ability to define the unknown quantities directly on the contour of the hole, not evaluating them on the whole surface of the shell. To obtain the kernels of the singular integral equations, integral representations of displacements and shallow isotropic shells static equations fundamental solutions were used. As the unknown functions, a combination of displacements, rotation angles and their derivatives were used. Analytical calculations are considerably simplified if it is assumed to take into account not four unknown functions on the contour, as it is customary, but five. In this paper it is chosen to use the differential equation which relates the unknown functions as the fifth equation of the boundary integral equations system. In order to obtain numerical solution of the problem, the method of mechanical quadratures for systems of integral equations and finite difference method for the fifth differential equation were used to reduce a problem to a system of linear algebraic equations. The stress concentration factors values depending on the isotropic shell curvature are given. Also the results were compared with other researchers.

Keywords: circular hole, isotropic shell, stress concentration factor, Fourier transformation, method of boundary integral equations.

Введение

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

Первым, кто направил свои усилия для решения задачи о напряженном состоянии вокруг кругового отверстия в круговой цилиндрической оболочке, был А.И. Лурье [1]. Оболочка была нагружена осевым растяжением и внутренним давлением. Его расчеты основывались на методе возмущений. Подобное исследование было проведено

Ю.А. Шевляковым [2]. Он исследовал напряжения в сферическом днище резервуара, ослабленном малым круговым отверстием в центре при осесимметрической нагрузке. Считалось, что контур отверстия нагружен перерезающими усилиями, которые компенсируют действие внутреннего давления, или подкреплен абсолютно жестким патрубком. Решения обоих авторов действительны только для малых значений параметра r0/Rh, где r0 - радиус отверстия; h - толщина оболочки; R -радиус кривизны срединной поверхности оболочки. V. Dyke [3], J.G. Lekkerkerker [4] исследовали круговую цилиндрическую оболочку с отверстием при различных условиях нагрузки. Использовав численные методы, каждый из них расширяет результаты, которые были получены предыдущими авторами, для больших значений параметра r0/Rh. Цилиндрическим оболочкам с разными видами подкрепления

отверстий были посвящены работы И.М. Пирогова [5].

Л.Б. Именитов рассмотрел сферическую оболочку в виде купола с одним жестко заделанным краем и произвольно расположенным непод-крепленным круговым отверстием. В его работе был предложен метод расчленения напряженного состояния [6]. А. А. Сяський, Е.И. Лунь занимались решением задачи о напряженном состоянии изотропной сферической оболочки с круговым отверстием, край которого подкреплен широким упругим кольцом [7].

В дальнейшем наиболее существенные результаты были получены в работах Г.Н. Савина, который сформулировал общую постановку задачи о напряженном состоянии оболочки с отверстием, которое задается в виде гладкого криволинейного контура [8].

Достаточно полный анализ проведенных исследований по концентрации напряжений вблизи отверстий в оболочках приведен в монографии [9]. Влияние поперечного сдвига оболочки с отверстием и ряда других нагрузок, приложенных к конструкциям, описывается в работах [10-16]. Конические оболочки с круговым вырезом рассматриваются в [17].

Оболочки положительной гауссовой кривизны с отверстием были рассмотрены в [18-20]. Задачи сводились к решению системы граничных интегральных уравнений. Одно из преимуществ использования метода граничных интегральных уравнений для исследования напряженно-деформированного состояния оболочек, ослабленных отверсти-

ем, состоит в возможности определять искомые величины непосредственно на контуре отверстия, не вычисляя их на всей поверхности оболочки. В качестве неизвестных функций использовались комбинации перемещений, углов поворота и их производных. Аналитические выкладки существенно упрощаются, если считать неизвестными на контуре не четыре функции, как это принято, а пять. В результате количество неизвестных превышает количество уравнений. Поэтому для обеспечения единственности решения необходимо добавить еще одно уравнение на контуре. В работе [20] предлагались разные варианты выбора пятого уравнения, но численные расчеты проводились только для случая, если задавалось пятое уравнение как компонент главного вектора. В данном исследовании в качестве пятого уравнения использовалось дифференциальное уравнение, связывающее неизвестные функции. Таким образом, первые четыре уравнения - интегральные, пятое - дифференциальное.

1. Постановка задачи

Рассматривается тонкая изотропная оболочка произвольной гауссовой кривизны и постоянной толщины h, которая находится под действием симметричной внешней нагрузки. Оболочка отнесена к системе ортогональных координат x, ус осями вдоль направлений главных кривизн оболочки, а ось z направлена по нормали к ней (рис. 1). Конструкция ослаблена неподкрепленным круговым отверстием радиусом r0 (контур L), расположенным в центре системы координат, уравнение отверстия имеет вид

x=r0 cos0, у=r0 sin 0,0e (0;2п].

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

Предполагается, что контур отверстия свободен от нагрузки на L.

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

2 I

Рис. 1. Тонкая оболочка с круговым отверстием

Граничные условия для возмущенного напряженного состояния имеют вид

Тп =-С Б„ =-С, Мп =-М*, Ит =-И'„, Ып =-N1,

где Тп,- мембранные усилия на контуре отверстия; Мп,Ипт - изгибающий и крутящий момент соответственно; Nn - перерезывающая сила.

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

Как показывают результаты экспериментальных исследований, возмущенное напряженное состояние возле отверстий имеет локальный характер и быстро затухает при удалении от отверстия. Поэтому в свое время Г.Н. Савин [8] предложил для описания возмущенного напряженно-деформированного состояния использовать уравнения теории оболочек с большим показателем изменяемости [22], которые совпадают с уравнениями пологих оболочек [23].

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

постоянными [24].

Основные соотношения теории пологих оболочек в случае изотропного материала принимают следующий вид [20]:

1) уравнения состояния исходя из гипотезы Кирхгофа-Лява записываются как

ЭТ дБ ^ дБ дТ2 -1+-= - X , -+-2:

Эх ду Эх ду

ЭМ1 ЭН лг дН ЭМ -^-~-Ы1 --1--:

-+

дх ду

+-

дх ду

---У.

-( Т + к 2Т2 )=-7, 1 1 1 2 2> дх ду

где X,У,7 - проекции внешней нагрузки на оси координат; Т1,Т2,Б -мембранные усилия; М1,М2,Н - изгибающие и крутящий моменты, Ы1,Ы2 - перерезывающие силы;

2) геометрические соотношения выражаются следующим обра-

зом:

ди

дд

е1 =--I к, —, е2 =--I к2 —, е12 =—

1 дх 12 ду 2 12 2

1 ( ди дд ^

ду дх

01 =—э—, 02 =-

= де1 =

Х1 = -Ч , X 2 =

дх де2

э— ¥

, Х12 ~

1 (де1 де 2 ^ =--—2

ду дх

дх ду 2

где е1, е2, е12 - компоненты тангенциальной деформации срединной поверхности оболочки; х1, X2, Х12 - компоненты изгибающей деформации; и,у, — - перемещения в направлениях х, у, г соответственно; е1,е 2 -

углы поворота нормали;

3) уравнение неразрывности деформации принимает вид

к0Х + к0Х I Э 2е1 — 2 Э 2е12 , Э 2е2 = 0. к1 Х1 + к2Х 2 + Л 2 2 Л -л + Л 2 =0;

Эх Эхду ду

4) соотношения упругости запишутся как

Ек Ек ЕН 71=-+ 2 ), Т2 = "-+е2 ), Б =-е1

1 1-у2

Н =В(1 V)Х12, М1 = В(Х1 +vx2), М2 +Х2),

где Б=-

ЕИ3

; Е - модуль Юнга; V - коэффициент Пуассона.

12

Дальнейшее исследование основано на использовании теории обобщенных функций, двумерного интегрального преобразования Фурье [20].

2. Система сингулярных интегральных уравнений

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

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

Обозначим через часть контура Ь с

концом в текущей точке 5 (рис. 2), где п, т

нормаль и касательная к контуру Ь5 в точке 5.

Граничные условия для возмущенного состояния запишем с учетом задания на контуре компонентов главного

— —

вектора Р и главного момента М усилий, которые действуют вдоль части контура Ь [18],

Рис. 2. Часть контура Ь

Рх =-Р*; Ру =-Р*; М« =М)!>; му) =-М<-;^ на Ь.

Неизвестные функции задаются через скачок функции при переходе через контур со стороны внешней нормали

¥ = ЕИ ГМ-ъЦЛ ¥2 = ЕИ(п М1; ¥з = ЕИГШ

V д ^2 У

V д Я у

я2 V д(

ЕИ ГЭ[02]^ ЕИ , т ^

¥4 =—

где п1, п2 - направляющие косинусы вектора внешней нормали к Ь

Я

вдоль осей х, у; Я1,Я2 - радиусы главных кривизн оболочки; А,=—-.

Я,

Система граничных интегральных уравнений для данной задачи, учитывая работу [20], где подробно изложено сведение подобных задач к системам граничных интегральных уравнений, запишется следующим образом:

\jiKj

Lj=I

где / - полудлина контура.

— — и

Компоненты главного вектора Р и главного момента М усилий приведены ниже:

Рх (s)=| Тхйу-Sdx; Ру (s)=18йу-Тгдх; Рг (s)=|Ы^+к2Р/у+кхРхсх,

Ьз Ь Ъ3

М^^)=-{Pzdy+кхМ[s-1dX; М«

Ь* Ls

М[s) ^ ) = | Р^У - Pydx.

С учетом (3) правая часть в (2) будет задана через производные выражений в (3):

F1 Рх'; F2 ^ )=Ру'; Fъ (s)=MX; F4 ^)=М'у.

Но при вычислении компонент главного вектора и момента в правой части системы (2) возникают 4 дополнительных константы с{ ,i=1,4, для определения которых используются дополнительные соотношения, следующие из условий однозначности перемещения и углов поворота при обходе контура Ь.

| (V: - С2 ( (У2 -Лх2) V з - 2Ьуу 4+2 У у 5

где с-

- с(ху 4 -Лууз))=0; |(хуз + ууА))=0,

|2 + "2"^

с

ds=0,

Ядра системы (2) имеют особенность типа Коши. Подробная методика получения ядер описывается в работе [21].

Таким образом, задача сведена к решению системы граничных интегральных уравнений относительно пяти неизвестных функций. Для решения поставленной задачи необходимо выбрать соответствующее количество уравнений, которое бы обеспечило единственность решения. Данный вопрос был затронут во многих работах [18-21]. В частности, в работе [20] предлагалось выбрать интегральные уравнения, соответствующие компонентам главного вектора Рх,Ру и главного

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

Численная реализация в [20] была проведена для случая, когда в качестве пятого уравнения задавались компоненты главного вектора Рг.

Существенным отличием данного исследования от предыдущих работ [18-21] является выбор в качестве пятого уравнения дифференциального уравнения и проведение численных расчетов с его использованием. Дифференциальное уравнение приведено ниже:

Г\

у С ёу3 ёу ^

ё у

5=от,.

ё е2 0

т Л I

+ 2(-П2У3 4

1 ёе 2 ёе \ 2тз 14

V У

у 5.

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

Численный расчет системы (2) основывался на использовании специальных квадратурных формул для интегралов типа Коши в случае, если неизвестные функции имеют корневую особенность на концах промежутка интегрирования [24]. Таким образом, каждое уравнение системы (2) сводилось к N линейно-алгебраическим уравнениям относительно значений неизвестных подынтегральных функций в выбранных точках (узлах интерполяционного полинома). Для сведения

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

Для обеспечения единственности решения были использованы еще два алгебраических уравнения, которые были получены из следующих интегральных уравнений:

=0; |=0.

Ь Ь

В результате получена система линейных алгебраических уравнений размерности (5N+4)х(5N+4).

3. Результаты численного исследования

Для определения напряженного состояния около кругового отверстия в оболочке были вычислены коэффициенты концентрации напряжения

ЕЪ 1

КТв = Те/РоЪ; Тв=-—2-ех; ех=—(¥1-П¥2); (РоЪ=рК/2).

1-У ЕЪ

Как показали результаты проведенных исследований, при выборе узлов разбиения достаточно ограничиться п=20, так как дальнейшее увеличение количества узлов разбиения практически не сказывается на результате.

Для сравнения с известными результатами рассматривалась оболочка нулевой гауссовой кривизны под действием осевого растяжения. Результаты сравнивались с работами [4, 9] и имеют хорошую сходимость. В работе [4] Ю. Ьеккегкегкег осуществил эксперимент для образца в виде цилиндрической оболочки с круговым отверстием, сделанной из стали. Нагрузка - осевое растяжение.

В табл. 1 приводятся результаты работы [4] в сравнении с коэффициентами концентрации напряжений, полученными в данном исследовании. Необходимо учитывать, что для моделирования задачи использовались различные теории, но, несмотря на это, погрешность вычислений не превосходит и 5 %.

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

нагружен перерезывающими усилиями, уравновешивающими действие внутреннего давления, имеют следующий вид:

Тп\ь =-рЯ/2; Мп\ь = 0; Кп\ь =-ръ/2.

Таблица 1

Сравнение коэффициентов концентрации напряжений с работой [4]

0 Результаты эксперимента [4] Теоретические результаты [4] Результаты, полученные в данной работе

0° 3,61 3,66 3,59

22° 2,62 2,78 2,73

90° -1,27 -1,2 -1,24

135° 0,87 0,91 0,88

180° 3,7 3,66 3,59

225° 0,67 0,91 0,88

270° -1,27 -1,2 -1,24

292° -0,59 -0,65 -0,72

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

в=сг0.

Таблица 2

Сравнение коэффициентов концентрации KTQ для оболочки,

находящейся под действием внутреннего давления, с результатами

других исследователей

в K'T , полученные в работе [2] K'T , полученные в работе [10] Результаты данного исследования

0,10 2,008 2,008 2,008

0,20 2,031 2,034 2,034

0,45 2,184 2,177 2,178

1,00 2,716 2,744 2,792

1,10 2,844 2,867 2,931

1,20 3,044 2,994 3,087

1,35 3,243 3,188 3,345

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

Рис. 3. Зависимость коэффициентов концентрации напряжений от кривизны оболочки

К 2,6

2,5

2,4

2,3

2,2

2,1

2

1,9

х = -

¿0,75

,-0,5

К 2,8

2,6

2,4

2,2

2

1,8

г1

0,75

~0;5

А,25

Х = 0

0,15 0,47 0,785 1,25 0 о,15 0,47 0,785 1,25

а

б

Рис. 4. Зависимость коэффициентов концентрации напряжений вдоль контура отверстия при заданных значениях параметра X: а - X < 0, б - X > 0

На рис. 3 показано влияние гауссовой кривизны на напряженное состояние оболочки. Здесь параметр X = 1 соответствует сферической оболочке; X = 0- цилиндрической оболочке; X = -1 - псевдосфере. Данные и последующие результаты приводятся для случая, когда контур отверстия нагружен перерезывающими усилиями, в результате действия которых Р2 = 0. Приведенные результаты соответствуют значению 0=0°.

На рис. 4 представлены результаты для случая в = 1.

Выводы

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

1) Максимальные значения коэффициентов концентрации напряжения KT при 0 = 0 достигаются для случая X = 1, что соответствует

сферической оболочке. Минимальные значения соответствуют случаю X = -1, что соответствует псевдосферической оболочке.

2) При увеличении параметра в значения коэффициентов концентрации напряжения возрастают. После значения X = 0 наблюдается монотонное возрастание .

3) Наибольшие значения KJ при в = 1 в точках 0 = 0 и 0 = п/2 соответствуют случаю X = 1 и принимают значение KJ = 2,79. Наименьшее значение KTQ в точке 0 = 0 соответствует случаю X = -1, где

= 2,35, а в точке 0 = п/2 соответствует случаю X = -0,25, где KTe = 2,00.

Библиографический список

1. Лурье А.И. Статика тонкостенных упругих оболочек. - Л.: Гос-техиздат, 1947. - 252 с.

2. Шевляков Ю.А. Напряжение в сферическом днище, ослабленном круговым вырезом // Инж. сб. - М: Наука, 1956. - Вып. 24. -С.226-230.

3. Dyke P. Stresses about a circular hole in a cylindrical shell // ATAA Journ. - 1967. - Vol. 5. - No. 9. - P. 87-91.

4. Lekkerkerker J.G. On the stress distribution in cylindrical shells weakened by a circular hole. - Delft: Vitgeverij Waltman, 1965.

5. Пирогов И.М. Влияние кривизны на распределение напряжений около отверстия в цилиндрической оболочке // Прикл. механика. -1965. - Т.1, № 12. - С. 116-119.

6. Именитов Л.Б. Задача о сферической оболочке с неподкреп-ленным отверстием // Инженерный журнал. - 1963. - Т. 3, № 1. -С. 93-100.

7. Сяський А.А., Лунь Е.И. Напряженное состояние изотропной сферической оболочки с круговым отверстием // Прикладная механика. -1974. - Т. 10, № 8. - С. 98-102.

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

9. Методы расчета оболочек: в 5 т. Т.1. Теория тонких оболочек, ослабленных отверстиями / А.Н. Гузь [и др.]. - Киев: Наук. думка, 1980. - 636 с.

10. Закора С.В., Чехов В.Н., Шнеренко К.И. О концентрации напряжений около кругового отверстия в трансверсально-изотропной сферической оболочке // Прикл. механика. - 2004. - Т. 40, № 12. -С. 99-106.

11. Loo Wen-da, Cheng Yao-shun. The effect of transverse shear deformation on stress concentration factors for shallow shells with a small circular hole // Applied Mathematics and Mechanics. - 1991. - Vol. 12(2). -Р.195-202.

12. Шевченко В.П., Закора С.В. О влиянии сдвиговой жесткости на напряженное состояние в транстропной сферической оболочке с двумя круговыми отверстиями при их сближении // Доповвд Нащонально! академп наук Украши. - 2010. - № 12. - С. 56-62

13. Хома И.Ю., Стрыгина О. А. Напряженное состояние пологой сферической оболочки с круговым отверстием, на поверхности которого заданы касательные напряжения // Теоретическая и прикладная механика. - 2010. - № 1 (47). - С. 62-68.

14. Nitin K.J. Analysis of Stress Concentration and Deflection in Isotropic and Orthotropic Rectangular Plates with Central Circular Hole under Transverse Static Loading // World Academy of Science, Engineering and Technology. - 2009. - Vol. 36. - P. 407-413.

15. Reissner E., Wan F.Y.M. Further considerations of stress concentration problems for twisted or sheared shallow spherical shells // International Journal Solids Structures. - 1994. - Vol. 31. - No. 16. - P. 2153-2165.

16. Hsu Chin-yun, Yeh Kai-yuan. The General Solution of Bending of a Spherical Thin Shallow Shell with a Circular Hole at the Center under Arbitrary Transverse Loads // Applied Mathematics and Mechanics. - 1980. -Vol. 1. - No. 3. - P. 341-356.

17. Чернышенко И.С., Сторожук Е.А., Харенко С.Б. Упругоплас-тическое состояние гибких конических оболочек с круговым вырезом

при осевом растяжении // Прикладная механика. - 2011. - Т. 47, № 6. -С. 86-92.

18. Довбня Е.Н. Деформация контура кругового отверстия в ор-тотропных оболочках при растягивающей нагрузке // Вюник Донец. ун-ту. Сер. А. - 2000. - № 1. - C. 51-55.

19. Довбня Е.Н. Напряженно-деформированное состояние ортот-ропной оболочки с круговым отверстием // Теорет. и прикл. механика. -1997. - № 27. - С. 154-158.

20. Довбня К.М. Розвиток методу граничних штегральних р1внянь в теорп ортотропних оболонок з розр1зами та отворами: дис. ... д-ра ф1з.-мат. наук. - Донецьк, 2001. - 362 с.

21. Довбня Е.Н. Система граничных интегральных уравнений для ортотропных оболочек нулевой и отрицательной кривизн, ослабленных разрезами и отверстиями // Вюник Донец. ун-ту. Сер. А: природничi науки. - 1998. - Вип. 2. - С. 45-52.

22. Гольденвейзер А. Л. Теория упругих тонких оболочек. - М.: Наука. - 1976. - 512 с.

23. Власов В.З. Избранные труды: в 3 т. Т. 1. - М.: Изд-во АН СССР, 1962. - 528 с.

24. Панасюк В.В., Саврук М.П., Дацышин А.П. Распределение напряжений около трещин в пластинах и оболочках. - Киев: Наук. думка. - 1976. - 444 с.

References

1. Lur'e A.I. Statics of Thin-walled Elastic Shells. State Publishing House of Technical and Theoretical Literature. Moscow, 1947; translation, AEC-tr-3798, Atomic Energy Commission, 1959.

2. Shevliakov Iu.A. Napriazhenie v sfericheskom dnishche, oslablen-nom krugovym vyrezom [Stresses in the spherical bottom weakened by a circular cut]. Inzhenernyj sbornik, 1956, vol. 24, pp. 226-230.

3. Dyke P. Stresses about a circular hole in a cylindrical shell. ATAA Journ., 1967, vol. 5, no. 9, pp. 87-91.

4. Lekkerkerker J.G. On the stress distribution in cylindrical shells weakened by a circular hole. Delft: Vitgeverij Waltman, 1965.

5. Pirogov I.M. Vliyanie krivizny na raspredelenie napriazhenij okolo otverstiya v tsilindricheskoj obolochke [Effect of curvature on the stress

distribution near the hole in the cylindrical shell]. Applied Mechanics, 1965, vol.1, no. 12, pp. 116-119.

6. Imenitov L.B. Zadacha o sfericheskoj obolochke s nepodkreplen-nym otverstiem [The problem of a spherical shell with an unsupported hole]. Inzhenernyi zhurnal, 1963, vol. 3, no. 1, pp. 93-100.

7. Syaskii A.A., Lun E.I. State of stress in an isotropic spherical shell with a circular hole. Prikladnaya mekhanika, 1974, vol. 10, no. 8, pp. 884887.

8. Savin G.N. Raspredelenie napryazhenij okolo otverstij [Stress distribution around holes]. Kiev: Naukova dumka, 1968. 888 p.

9. Guz A.N., Chernyshenko I.S., Chekhov Val.N., Chekhov Vik.N., Shnerenko K.I. Metody rascheta obolochek: v 5 t. T.1. Teoriya tonkikh obolochek, oslablennykh otverstiyami [Methods of calculation of shells. In 5 vols. Vol. 1. The theory of thin shells weakened by holes]. Kiev: Naukova dumka, 1980, 636 p.

10. Zakora S.V., Chekhov V.N., Shnerenko K.I. Stress concentration around a circular hole in a transversely isotropic spherical shell. Prikladnaya mekhanika, 2004, vol. 40, no. 12, pp. 1391-1397.

11. Loo Wen-da, Cheng Yao-shun. The effect of transverse shear deformation on stress concentration factors for shallow shells with a small circular hole. Prikladnaya matematika i mekhanika, 1991, no. 12(2), pp. 195-202.

12. Shevchenko V.P., Zakora S.V. O vliyanii sdvigovoj zhestkosti na napryazhennoe sostoyanie v transtropnoj sfericheskoj obolochke s dvumya krugovymi otverstiyami pri ikh sblizhenii [On the effect of shear stiffness on the stress state in transtropic spherical shell with two circular holes on their approach]. Dopovidi Natsional'noi akademii nauk Ukraini, 2010, no. 12, pp. 56-62.

13. Khoma I.Yu., Starygina O.A. Napryazhennoe sostoyanie pologoj sfericheskoj obolochki s krugovym otverstiem, na poverkhnosti kotorogo zadany kasatelnye napryazheniya [Stress state of shallow spherical shell with a circular hole where shear stresses are given on the surface]. Teo-reticheskaya iprikladnaya mekhanika, 2010, no. 1 (47), pp. 62-68.

14. Nitin K.J. Analysis of Stress Concentration and Deflection in Isotropic and Orthotropic Rectangular Plates with Central Circular Hole under Transverse Static Loading. World Academy of Science, Engineering and Technology, 2009, vol. 36, pp. 407-413.

E.H. floeÓHX, H.A. KpynKO

15. Reissner E., Wan F.Y.M. Further considerations of stress concentration problems for twisted or sheared shallow spherical shells. International Journal Solids Structures, 1994, vol. 31, no. 16, pp. 2153-2165.

16. Hsu Chin-yun, Yeh Kai-yuan. The General Solution of Bending of a Spherical Thin Shallow Shell with a Circular Hole at the Center under Arbitrary Transverse Loads. Applied Mathematics and Mechanics, 1980, vol. 1, no. 3, pp. 341-356.

17. Chernyshenko I.S., Storozhuk E.A., Kharenko S.B. Elastoplastic state of flexible conical shells with a circular hole under axial tension. International Applied Mechanics, 2011, vol. 47, no. 6, pp. 679-684.

18. Dovbnya E.N. Deformatsiya kontura krugovogo otverstiya v ortotropnykh obolochkakh pri rastyagivaiushchej nagruzke [Deformation of the contour of a circular hole in orthotropic shells under tensile load]. Visnik Donetskogo universitetu. Ser.A, 2000, no. 1, pp. 51-55.

19. Dovbnya E.N. Napryazhenno-deformirovannoe sostoyanie orto-tropnoy obolochki s krugovym otverstiem [Stress-strain state of orthotropic shell with a circular hole]. Teoreticheskaya iprikladnaya mekhanika, 1997, no. 27, pp. 154-158.

20. Dovbnya K.M. Rozvitok metodu granichnikh integral'nikh rivnian' v teorii ortotropnikh obolonok z rozrizami ta otvorami [The development of the method of boundary integral equations in the theory of orthotropic shells with holes and cuts]. Thesis of doctors degree dissertation physics and mathematics sciences. Donetsk, 2001, 362 p.

21. Dovbnya E.N. Sistema granichnykh integral'nykh uravnenii dlia ortotropnykh obolochek nulevoi i otritsatel'noi krivizn, oslablennykh raz-rezami i otverstiiami [System of boundary integral equations for orthotropic shells with zero and negative curvature, weakened by cuts and holes]. Visnik Donets'kogo universitetu. Ser. A: prirodnichi nauki, 1998, vol. 2, pp. 45-52.

22. Goldenveizer A.L. Teoriya uprugikh tonkikh obolochek [Theory of elastic thin shells]. Moscow: Nauka, 1976. 512 p.

23. Vlasov V.Z. Izbrannye trudy: v 3 t., T. 1. [Selected Works, in 3 vols, vol. 1]. Moscow: Akademiya nauk SSSR, 1962. 528 p.

24. Panasiuk V.V., Savruk M.P., Datsyshin A.P. Raspredelenie na-pryazhenii okolo treshchin v plastinakh i obolochkakh [The stress distribution around cracks in plates and shells]. Kiev: Naukova dumka, 1976. 444 p.

Об авторах

Довбня Екатерина Николаевна (Донецк, Украина) - доктор физико-математических наук, профессор, профессор кафедры прикладной механики и компьютерных технологий Донецкого национального университета (83055, Украина, г. Донецк, ул. Университетская, 24).

Крупко Наталья Андреевна (Донецк, Украина) - аспирант кафедры прикладной механики и компьютерных технологий Донецкого национального университета (83055, Украина, г. Донецк, ул. Университетская, 24, e-mail: nataliekrupko@gmail.com).

About the authors

Dovbnya Ekaterina Nikolaevna (Donetsk, Ukraine) - Doctor of Physics and Mathematics Sciences, Professor, Professor of the Department of Applied Mechanics and Computer Technology Donetsk National University (24, Universitetskaya str., 83055, Donetsk, Ukraine).

Krupko Nataliia Andreevna (Donetsk, Ukraine) - Doctoral student Department of Applied Mechanics and Computer Technology Donetsk National University (24 Universitetskaya str., 83055, Donetsk, Ukraine, e-mail: nataliekrupko@gmail.com).

Получено 10.02.2014

Просьба ссылаться на эту статью в русскоязычных источниках следующим образом:

Довбня Е.Н., Крупко Н.А. Влияние кругового отверстия на напряженное состояние оболочки произвольной гауссовой кривизны // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2013. -№ 1. - С. 108-125.

Please cite this article in English as:

Dovbnya E.N., Krupko N.A. Influence of circular hole on the shell stress state for arbitrary Gaussian curvature. PNRPUMechanics Bulletin. 2013. No. 1. P. 108-125.