Full Text

Введение В работе предлагается новый аналитический подход к выводу напряжений цилиндрической оболочки с круговым отверстием при растяжении вдоль образующей. Постановка задачи и решение методом разложения по малому параметру β, характеризующим отношение радиуса отверстия, толщины оболочки и радиуса кривизны, принадлежат А.И. Лурье [1, 25] (40-е гг. XX в.). Далее в середине 60-70-х гг. всплеск интереса к этой проблеме происходит не только среди советских, но и среди зарубежных ученых, нашедших ошибку в постановке краевых условий на границе кругового отверстия и пересмотревших эту задачу либо тем же способом (Пирогов, Murthy, Naghdi, Eringen [4, 6, 17-23, 26, 30]), либо численным методом коллокаций (Lekkerkerker, Eringen, Naghdi, VanDyke [3, 5, 7]). Однако предложенное аналитическое решение было крайне трудоемким и работало для диапазона значений параметра β, мало отличающегося от плоской задачи Кирша, а результаты, полученные методом коллокаций, разнились [9]. Также предпринимались попытки подойти к этой задаче энергетическим методом [23, 29] и методом комплексных переменных для произвольных отверстий [16, 27, 28]. Ресурсы рассмотренных методов себя исчерпали, не дав удобного решения, а альтернативных методов пока предложено не было. Последующие работы опираются на компьютерное моделирование, в частности основанное на методе конечных элементов [10-15]. Однако по сей день актуальность и применяемость этой задачи остаются высокими [31-37], особенно в области авиапромышленности, и именно аналитически найденное поле напряжений в области отверстия может дать толчок для фундаментального изучения вопросов разрушения и устойчивости. В настоящей работе представлена новая идея решения, которая позволяет легко находить численные значения напряжений и открывает перспективы их аналитического исследования, а также снимает математические ограничения на значения параметра β. В работе особое внимание уделено технике решения задачи и строгой математический постановке. 1. Постановка задачи Рассматривается цилиндрическая оболочка с круговым отверстием, растягиваемая вдоль образующей х некоторой силой р на бесконечности. В задаче будут использоваться следующие обозначения: параметр - основной параметр, отвечающий за отношение геометрических и физических параметров, в том числе за кривизну оболочки; r0 - радиус отверстия (не умаляя общности, можно принять r0 за единицу измерения, т.е. r0 = 1); R, h - радиус кривизны и толщина оболочки соответственно; v - коэффициент Пуассона; E - модуль Юнга. Заметим, что предельный случай при приводит нас к задаче Кирша. Также вводится функция , зависящая от прогиба w и функции напряжений U [1, 25].Связь между усилиями T и функцией U задается следующим образом: а напряжение срединной поверхности тонкой оболочки как В [1, 2] показано, что система уравнений равновесия оболочки сводится к следующему уравнению: (1) Сформулируем условие задачи полностью: найти функцию, удовлетворяющую уравнению (1) и всем граничным условиям - на бесконечности (2) - на границе кругового отверстия в полярной системе координат (3) Здесь - момент, - обобщенное граничное условие на свободном крае [1]. 2. Решение Несмотря на то что метод решения уравнения (1) хорошо известен [1, 2, 4, 6, 9, 24, 25], некоторым техническим деталям не было уделено должного внимания. Рассмотрим два коммутирующих линейных оператора и , где Тогда уравнение (1) можно записать в виде Другими словами, задача сводится к нахождению ядра произведения Из того, что операторы коммутируют, следует Нахождение по отдельности решений уравнений и с последующей возможностью найти их сумму существенно упрощает решение исходного уравнения, так как понижает его порядок. Однако важно отметить, что сумма ядер не совпадает с множеством всех решений уравнения (1), что может приводить к потере решений. Поэтому такой метод может быть использован для доказательства существования решения и его конструктивного нахождения, а исследование единственности должно быть проведено отдельно. С помощью замены и разделения переменных легко установить [1, 2, 25], что решениями являются функции Выбор функций Ханкеля при построении решения связан с тем, что это единственные из бесселевых функций, стремящиеся к нулю в бесконечно далекой комплексной точке [8]: Поскольку имеет аргумент функции заведомо не удовлетворяют граничным условиям, так как прогиб на бесконечности. В то же время, руководствуясь [8], можно вывести, что Заметим, что первые три граничных условия системы (2) ставятся не относительно функции U, а относительно ее вторых производных. Поэтому необходимо убедиться в том, что не только потенциал, но и напряжения стремятся к нулю при больших r. Это верно, так как из рекуррентных соотношений для функций Бесселя [8] следует, что является линейной комбинацией функций такого же вида Кроме того, непосредственные вычисления показывают, что является линейной комбинацией функций вида Тем самым решение удовлетворяет всем граничным условиям на бесконечности. Такую проверку не удалось найти ни в одной из известных работ. Отсутствие полной четкой формулировки задачи и недостаточное внимание к поведению функций на бесконечности привело к тому, что полученные в работах [1, 6] решения граничным условиям не удовлетворяют. Если перейти к тригонометрической форме записи и учесть симметрию кругового отверстия, получается, что решение задачи (1)-(3) можно искать в следующем виде для четных и нечетных n [1, 2]: (4) Функция Ф является решением уравнения математической физики (1), при этом Ф удовлетворяет граничным условиям (2). Остается лишь найти коэффициенты и из граничных условий (3). На этом этапе авторы предшествующих работ столкнулись с наибольшими затруднениями [подробнее см. 9]. Именно метод поиска коэффициентов является основной содержательной частью данной работы. Как уже отмечалось выше, ранее предлагались три различных метода поиска коэффициентов: 1. В работе [1] и разложены в ряд по степеням Это позволило разделить переменные r и Однако вычисления получились слишком громоздкими, и появилось сильное ограничение на параметр Этот подход не оставляет возможности аналитического исследования напряжений и не позволяет получить численных значений напряжений для всех 2. Численный поиск коэффициентов (методом коллокаций). Имеет мало перспектив для аналитического исследования напряжений. Полученные численные результаты отличались друг от друга (причем даже в статьях одного автора). Напряжения для всех были получены только в работе [7]. 3. Подставить решение в виде (4) в граничные условия. Получившиеся крайне громоздкие соотношения не позволяют провести аналитическое исследование, а возможность получения численных результатов представляется сомнительной. 3. Новый подход Основная идея заключается в разделении переменных r и в каждой базисной функции. Только в отличие от [1] для достижения этой цели предлагается разложение в тригонометрический ряд Фурье. Из известного разложения в ряд Лорана [8] для получается Замена n на -n приводит к Если сложить оба полученных равенства, получится Теперь можно заменить α на -α в последней формуле и сложить оба равенства при четных n, а при нечетных n вычесть из одного другое: В результате четные и нечетные базисные функции можно записать одной общей формулой (5) где - целая часть числа и в знаменателе введен нормирующий коэффициент Последнее сделано для того, чтобы численные значения неизвестных коэффициентов имели умеренные значения, с которыми удобно работать. Далее для удобства вводятся обозначения для коэффициентов Фурье в тригонометрическом разложении базисной функции Тогда (5) примет вид (6) Теперь решение (4) может быть записано в виде (7) удобном для подстановки в граничные условия (3). 4. Граничные условия Первое граничное условие в полярных координатах приводит к уравнению (8) Коэффициенты при косинусах дают нам следующую систему: (8а) Второе граничное условие дает уравнение (9) Коэффициенты при синусах: (9а) Из третьего граничного условия получаем (10) Коэффициенты при косинусах: (10а) Из четвертого граничного условия получаем (11) Коэффициенты при косинусах: (11а) 5. Исследование системы Во второй паре систем (10a-11a) можно домножить выражения, стоящие под знаком суммы, на Тогда во все уравнения четырех систем будет входить только мнимая часть от суммы. Из всех полученных систем составим общую линейную систему с бесконечным числом неизвестных и уравнений. Предварительно сделаем некоторые элементарные преобразования: 1. Для уравнения систем (8a) и (9а) запишем следующим образом: Получилась однородная система относительно неизвестных и , имеющая определитель, отличный от нуля а следовательно, имеющая только тривиальное решение, т.е. (12) 2. Для те же уравнения дают неоднородную систему Разрешая ее относительно неизвестных, получаем (13) 3. Введем обозначения Для того чтобы избавиться от в выражении можно вычесть из последнего уравнения системы (11а) последнее уравнение (10а): (14) Важно отметить, что первое уравнение системы (8а), как доказано авторами, является следствием четырех уравнений: двух уравнений системы (13) и двух уравнения системы (14) при Доказательство строится на свойствах функций Бесселя [8] и его идея изложена в аппендиксе. Факт линейной зависимости уравнений, но непонимание - каких именно ставили в тупик авторов предшествующих работ [6]. Все уравнения систем (12)-(14) можно записать в матричной форме (табл. 1). Таким образом, получена линейная система с бесконечным числом уравнений и неизвестных. При этом элементы бесконечной матрицы системы, существенно отличающиеся от нуля, расположены вблизи главной диагонали (аналогия с диагональным преобладанием), так как значения от которых зависят все элементы матрицы системы, с ростом k становятся очень малы, а именно рядом с главной диагональю расположены элементы, у которых индекс близок к нулю. Пример для представлен в табл. 2. Следовательно, подматрица, составленная из первых 4N строк и столбцов, будет иметь ненулевой определитель, и при решении системы однозначно будут найдены коэффициенты для первых 2N базисных функций. При этом, как показывают расчеты, при увеличении N первые найденные коэффициенты практически не изменяются, а коэффициенты при базисных функциях с большими индексами стремятся к нулю. Этот метод не имеет математических ограничений на значения основного параметра β. С точки зрения механики эта модель применима для диапазона [2]. Таблица 1 Системы (12)-(14) в матричной форме Table 1 Systems (12)-(14) in the matrix form n l 0 Im Re 1 Im Re 2 Im Re 3 Im Re Неизв. Свободн. 0 … 0 0 0 1 1 1 0 1 0 2 0 2 0 … Таблица 2 Численный пример для Table 2 The numerical solution for n: 0 1 2 3 4 5 l: Im Re Im Re Im Re Im Re Im Re Im Re 0 0,003 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0,107 -0,011 0 0 0 1 -0,214 0,022 0 0,001 -0,001 1 -0,44 0 0,047 0,002 0,002 0 1 -0,44 0 0,047 0,002 0,002 0 2 0 0 0 0,107 -0,022 1 0,105 0,107 0 2 0 0 -0,214 0,06 -4 -0,42 -0,428 0 2 0 0 -0,063 14 0,014 1,493 -1,475 0 2 0 -8,49 41,99 1,528 4,525 -4,378 0 3 0 0 0 0 0 0 0,0003 -0,011 0 0,105 0,107 3 0 0 0 0 0 0 0 0 0,022 0 -0,42 -0,428 0 6. Обсуждение результатов Найденные коэффициенты можно подставить в (4). При этом любая конечная частичная сумма является точным решением уравнения матфизики в области (в отличие, например, от решения данной задачи методом Ритца) и удовлетворяет граничным условиям на бесконечности. Причем, как показывают расчеты, краевые условия на границе отверстия выполняются достаточно точно для любого уже при 18 базисных функциях, коэффициенты которых найдены из редуцированной системы. Так, максимальное отклонение граничных условий от нуля для составляет не более , а при не более Естественно, что с ростом β максимальное отклонение увеличивается: для большей точности можно взять 24 базисные функции для больших значений β, и тогда отклонение составит не более При возрастании β увеличивается число базисных функций, существенно влияющих на ответ, т.е. растут базисные коэффициенты при больших n. Полученная система позволяет построить графики зависимости коэффициентов от β. Рис. 1. Графики коэффициентов an и bn в зависимости от Fig. 1. Graphs of coefficients an and bn depending on По графикам рис. 1 видно, как при возрастании β начинают себя проявлять коэффициенты например, коэффициенты при равны нулю, а далее начинают влиять на ответ (однако их порядок гораздо меньше, чем у первых коэффициентов). В табл. 3 указано количество базисных функций для получения напряжения с точностью до 10-3 при различных значениях β. Таблица 3 Количество базисных функций при различных значениях β Table 3 Number of basic functions with different values of β Результаты, изображенные на графике рис. 2, полностью совпадают с результатами, полученными в работе [7] методом коллокаций. Ранее в работах разных авторов результаты разнились, и оставалось неясным, на какие результаты опираться. Теперь удалось найти аналитический метод, простой в реализации и дающий достоверные результаты и возможности дальнейшего исследования напряжений. Рис. 2. Сравнение результатов, полученными авторами (слева) и в работе [7] (справа) Fig. 2. Comparison of results obtained by the authors (on the left) and in work [7] (on the right) Приводим сравнение результатов с известным асимптотическим анализом задачи при малых значениях параметра β и (табл. 4). Ниже представлены графики зависимости напряжений от r, по которым видно, что напряжения меняются в пределах значений даже для большого значения (рис. 3, 4). Таблица 4 Сравнение результатов по с интервалом 10° Table 4 Comparison of results according to with an interval 10° Угол Наш метод Метод Лурье* Наш метод Метод Лурье* 0 10° 20° 30° 40 50° 60° 70° 80° 90° -1,01434 -0,892895 -0,543185 -0,00732409 0,650131 1,34993 2,00767 0,54396 2,89405 3,01564 -1,01434 -0,892897 -0,543185 -0,00732314 0,65013 1,34993 0,00766 2,54396 2,89405 3,01564 -1,09003 -0,965589 -0,606287 -0,0529143 0,630768 1,36436 2,05954 2,63051 3,00529 3,13585 -1,08888 -0,964564 -0,6056 -0,0526668 0,630601 1,3639 2,05891 2,62984 3,00462 3,13519 *В данной работе разложение посчитано до β6. а б Рис. 3. Зависимость напряжений от r для β = 4 при (а) и (б) Fig. 3. Dependence of stresses on r for β = 4 when (a) and (b) Рис. 4. Поле напряжений от координат при и Fig. 4. Stress field from coordinates when and

About the authors

S. V Kashtanova

Institute for Problems in Mechanical Engineering of the Russian Academy of Sciences

A. V Rzhonsnitskiy

Saint-Petersburg State Institute of Technology

References

Statistics

Views

Abstract - 337

PDF (Russian) - 234