Full Text

Введение В работе рассматривается пространственная задача о равновесном состоянии мягкой сетчатой оболочки при наличии внешней точечной нагрузки, сосредоточенной в некоторой точке. Под сетчатой понимается оболочка, силовой основой которой является сетка, образованная двумя семействами взаимно перекрещивающихся, абсолютно гибких, упругих нитей. Предполагается, что функции, описывающие физические соотношения в нитях, являются непрерывными, неубывающими и имеют линейный рост на бесконечности. Решаемые здесь задачи весьма часто возникают в механике, медицине, при проектировании различного рода конструкций, изготовленных из тканевых или пленочных материалов (см., например, [1–13]). Задачи об определении положения равновесия мягкой сетчатой оболочки рассматривались и ранее (см., например, [14–18]). При этом обобщенные задачи формулировались в виде уравнений или вариационных неравенств с оператором, действующим в случае линейного роста функций, описывающих физические соотношения в нитях, из соболевского пространства в сопряженное с ним, и соответственно рассматривается ситуация, когда функция, описывающая плотность внешних источников, определяет линейный непрерывный функционал на V. В настоящей работе проводится исследование задач теории мягких сетчатых оболочек с менее гладкой правой частью: каждый точечный источник, в неодномерном случае моделируется дельта-функцией Дирака, которая не принадлежит пространству, сопряженному с V. Обойти указанную выше трудность удалось благодаря аддитивному выделению особенности, связанной с дельта-функцией. Обобщенная постановка задачи сформулирована в виде интегрального тождества относительно функции Затем введена вспомогательная задача с правой частью, задаваемой дельта-функцией. Для вспомогательной задачи известно решение в явном виде. Благодаря этому обобщенная постановка свелась к нахождению решения операторного уравнения в V. Установлены свойства оператора, входящего в это уравнение – ограниченность, непрерывность, монотонность, коэрцитивность, что дало возможность применить для доказательства теоремы существования известные результаты теории монотонных операторов (см., например, [19, 20]). Доказано, что множество решений обобщенной задачи не пусто, выпукло и замкнуто. Вообще говоря, решение обобщенной задачи не единственно. Однако установлено, что усилия в нитях определяются однозначно. Отметим, что аналогичный подход был использован при рассмотрении стационарных задач фильтрации несжимаемой жидкости, следующей закону фильтрации с предельным градиентом, при наличии точечного источника [21, 22]. Построены конечномерные аппроксимации задачи, исследована их сходимость. Для решения задач использован разработанный в работе [23] двухслойный итерационный метод. Данный метод был реализован численно. Проведенные для модельных задач численные эксперименты подтвердили эффективность итерационного метода. 1. Постановка задачи Рассмотрим пространственную задачу о равновесном состоянии мягкой сетчатой оболочки при наличии точечной нагрузки. Под сетчатой понимается оболочка, силовой основой которой является сетка, образованная двумя семействами взаимно перекрещивающихся, абсолютно гибких, упругих нитей. Предполагается, что узлы сети фиксированы, материал, заполняющий промежутки между нитями, не сопротивляется деформации, и ни в начальном состоянии, ни в процессе деформации соседние нити не соприкасаются. Ячейки сети считаются малыми и не сопротивляющимися сдвиговым деформациям. Деформации и перемещения допускаются конечными. Введем в пространстве декартову систему координат Считаем, что в недеформированном состоянии оболочка может быть описана поверхностью , где – лагранжевы координаты; – ограниченная область из с непрерывной по Липшицу границей . Предполагаем, что функция удовлетворяет условиям: , для любого справедливо неравенство . Лагранжевы координаты выберем так, что координатные линии сонаправлены с нитями, образующими оболочку. Через обозначим функцию, описывающую поверхность оболочки в деформированном состоянии, через – дискриминант метрического тензора поверхности деформированной оболочки. Здесь использованы обозначения: , ; , и – векторное, скалярное произведения и норма в соответственно. Введем также следующие обозначения (для ): ; – вектора, образующие ковариантный локальный базис на деформированной поверхности; , – параметры Ляме соответственно недеформированной и деформированной поверхности оболочки; – векторы, образующие контрвариантный локальный базис на деформированной поверхности: . Обозначим через внутреннюю силу, действующую на единицу длины -й координатной линии () деформированной оболочки с той стороны оболочки, куда направлен вектор , , через – коэффициенты разложения этой плотности сил по единичным векторам локального базиса: . Тогда ковариантные компоненты тензора напряжений связаны с погонными усилиями соотношениями (см. [24, с. 50]). Для сетчатой оболочки в силу того, что в выбранной лагранжевой системе координат направления осей совпадают с направлениями нитей, имеем [25]: (то есть ячейка сети не оказывает сопротивления повороту нитей в узлах скрепления), , где – относительные степени удлинения. Здесь – функции, характеризующие физические свойства нитей; – количество нитей, сонаправленных с -й координатной осью, на единицу длины -й координатной оси в недеформированном состоянии. Эти функции определены конструкцией сетчатой оболочки, и относительно них считаем выполненными условия , при (1) при *, (2) существуют постоянные такие, что , для всех , . (3) Заметим, что, вообще говоря, направления не являются главными для тензора , хотя смешанные компоненты и равны нулю. Уравнение равновесия оболочки, находящейся под воздействием внешних сил, в декартовой системе координат имеет следующий вид (см. [24, с. 88]): (4) где – векторы плотности соответственно поверхностной и массовой нагрузок; g – плотность материала оболочки в деформированном состоянии. В силу закона сохранения массы , где – заданная плотность материала недеформированной оболочки. Будем считать, что плотность материала равна единице () и . Поверхностная нагрузка предполагается равной нулю: . Поэтому в силу того, что , , уравнение (4) примет следующий вид: (5) Далее будем считать края оболочки закрепленными: (6) Вариационная постановка задачи (5), (6) в перемещениях состоит в нахождении вектор-функции такой, что (7) Предполагаем, что массовая нагрузка, действующая на оболочку, сосредоточена во внутренней точке множества и имеет интенсивность , , . Приведем обобщенную формулировку рассматриваемой задачи. Определим пространство как пополнение пространства по норме . Под обобщенным решением задачи (7) будем понимать функцию такую, что (8) 2. Исследование разрешимости задачи Исследование разрешимости задачи (8) основано на сведении ее к операторному уравнению в гильбертовом пространстве путем введения вспомогательной задачи, снимающей особенность в правой части, связанной с наличием точечной нагрузки, для которой построено точное решение. Определим пространство как пополнение пространства по норме , где , , и обозначим . Рассмотрим вспомогательную задачу поиска функции такой, что (9) где – постоянная из условия (1), d – дельта-функция Дирака. Решение задачи (9) существует, и для каждого (см., например [26, стр. 163]. Поэтому . Вариационная постановка задачи (9) выглядит следующим образом: Найти : (10) Решение задачи (9) будем искать в виде , где – неизвестная функция, а – решение задачи (10). Введем далее операторы , , , по формулам , , , . При этом , , – усилия в нитях. С учетом введенных обозначений задача (8) сводится к поиску функции такой, что (11) Справедлива Теорема 1. Задача (11) имеет непустое, выпуклое, замкнутое множество решений. Если – решение (8), то , где – решение задачи (11); – решение задачи (10). Определим на форму , где , . Нетрудно проверить, что форма непрерывна и линейна по второму аргументу. Поэтому по теореме Рисса-Фишера эта форма порождает оператор , . Таким образом, обобщенная задача (12) эквивалентна операторному уравнению . (12) При выполнении условий, наложенных на функции , можно доказать, что оператор является ограниченным, непрерывным, монотонным и коэрцитивным, поэтому существование решения уравнения (12), выпуклость и замкнутость множества его решений следует из общих результатов теории монотонных операторов (см., например [19, 20]). По аналогии с [21, 22] устанавливается справедливость второго утверждения теоремы 1. Отметим, что, вообще говоря, решение задачи (12) не единственно. Однако имеет место Теорема 2. Пусть выполнены условия (1)–(3), кроме того, функции липшиц-непрерывны: (13) Тогда, если – решение задачи (12), то усилия в нитях , , определяются единственным образом. Справедливость теоремы устанавливается по аналогии с [27]. 3. Конечномерные аппроксимации задачи При построении конечномерных аппроксимаций для простоты будем предполагать, что область является многоугольником. Пусть – некоторое семейство (триангуляция) треугольников , обладающее свойствами: 1) ; 2) 3) для любых множество если не пусто, то является либо стороной, либо вершиной одновременно треугольников ; 4) триангуляция является регулярной [28, 29]: , . Здесь – диаметр наименьшего шара, содержащего ; – наибольшего шара, содержащегося в наибольшем шаре, содержащего в . Определим теперь конечномерные пространства и , ассоциированные c триангуляцией : , где – пространство непрерывных функций, являющихся на каждом из полиномами не выше первого порядка; – пространство непрерывных функций из , равных нулю на границе множества . Очевидно, что , , и, значит, . При этом (см., например, [30]) при . Задаче (12) сопоставим следующую задачу: Найти : . (14) Справедливы следующие результаты: Теорема 3. Пусть выполнены условия (1)–(3). Тогда задача (14) имеет непустое, выпуклое и замкнутое множество решений, множество решений задачи (14) равномерно ограничено по : (15) где константа не зависит от . Если подпоследовательность слабо сходится к в при , то – решение задачи (14). Так же, как и при доказательстве теоремы 1, используя монотонность, непрерывность и коэрцитивность оператора , получаем существование непустого, выпуклого и замкнутого множества решений задачи (14). В силу коэрцитивности оператора множество решений задачи (14) равномерно ограничено по , т.е. справедлива априорная оценка (15). Из этой оценки следует существование подпоследовательности , слабо сходящейся к некоторому элементу в при . Используя монотонность и непрерывность, а значит, псевдомонотонность [19] оператора , получаем стандартным образом (см., например [19, 30]), что – решение задачи (14). Теорема 4. Пусть выполнены условия (1)–(3), (13), решение задачи (12) принадлежит пространству . Тогда справедлива оценка с константой , не зависящей от . 4. Итерационный метод В заключение отметим, что при выполнении условий (1)–(3), (13), наложенных на функции , оператор является потенциальным, липшиц-непрерывным и обратно сильно монотонным [31, 32]. Поэтому для решения задачи (12) можно использовать следующий итерационный процесс [23]. Пусть – произвольный элемент из . Для определим вектор как решение уравнения , (16) где – итерационный параметр. Теорема 5. Пусть выполнены условия (1)–(3), (13), . Тогда итерационная последовательность , построенная согласно (16), сходится слабо в к некоторому решению задачи (12). Указанный алгоритм решения задачи реализован в виде комплекса программ в среде MATLAB. Программа выполнена в соответствии с модульным принципом, что позволило осуществить раздельное программирование, отладку и тестирование составных частей пакета программ, а также простую модернизацию и настройку пакета на решение задач различного уровня сложности. Проведенные для модельных задач численные эксперименты подтвердили эффективность итерационного метода. В качестве области выбирался квадрат , функции полагались равными внешняя нагрузка была сосредоточена в точке , интенсивность нагрузки полагалась В качестве начального приближения в итерационном методе (16) выбирались функции, соответствующие недеформированной оболочке. Рис. 1. Зависимость числа итераций от итерационного параметра τ В результате численных экспериментов было найдено оптимальное (по количеству итераций) значение итерационного параметра (рис. 1). Это значение практически не зависит от числа разбиений (рис. 2). Значение максимального прогиба (в центре квадрата) также не зависит от количества точек сетки (рис. 3). На рис. 4 представлена форма деформированной оболочки. Рис. 2. Зависимость оптимального числа итераций от количества точек сетки Рис. 3. Зависимость прогиба в центральной точке от итерационного параметра Рис. 4. Форма деформированной оболочки

About the authors

Ildar Burkhanovich Badriev

Kazan Federal University

Email: ildar.badriev@ksu.ru
18, Kremlevskaja st., 420008, Kazan, Russian Federation Doctor of Physical and Mathematical Sciences, Professor, Department of Calculation Mathematics

Victor Victorovich Banderov

Kazan Federal University

Email: Victor.banderov@ksu.ru
18, Kremlevskaja st., 420008, Kazan, Russian Federation Ph.D. in Physical and Mathematical Sciences, Deputy Director for Research and Innovations of Institute of Computer Mathematics and Information Technologies

Oleg Anatol’evich Zadvornov

Kazan Federal University

Email: Oleg.Zadvornov@ksu.ru
18, Kremlevskaja st., 420008, Kazan, Russian Federation Doctor of Physical and Mathematical Sciences, Professor, Head of Department of Calculation Mathematics

References

Statistics

Views

Abstract - 113

PDF (Russian) - 48