


Том 89, № 1 (2025)
Статьи
Софья Васильевна Ковалевская (К 175-летию со дня рождения)



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



Преобразование нестационарных уравнений Навье–Стокса вязкой сжимаемой жидкости при произвольном конформном отображении
Аннотация
Показано, что при произвольном конформном отображении области двумерного течения вязкой сжимаемой жидкости циркуляция скорости и расход жидкости на любом замкнутом или незамкнутом контуре сохраняются. Выведены преобразованные нестационарные уравнения Навье–Стокса, неразрывности и баланса тепла, которым подчиняются аэродинамические параметры в отображенной области.



Моделирование эксперимента ONERA с дозвуковым предварительно перемешанным турбулентным горением в канале с обратным уступом
Аннотация
Рассматривается предварительно перемешанное дозвуковое турбулентное горение метано-воздушной смеси в канале с обратным уступом (P. Magre и др., ONERA, 1975–1989). В экспериментах воспроизведены базовые физические механизмы, характерные для процессов горения в газотурбинных установках. Дан краткий обзор предыдущих работ по численному моделированию этих экспериментов. Представлены новые результаты численного исследования режима стабилизированного горения в данной экспериментальной установке. Описан выбор модели турбулентности и ее влияние на структуру течения. Несколько подходов к описанию турбулентного горения из класса PaSR (Partially Stirred Reactor – модели реактора частичного перемешивания) сравниваются с квазиламинарным подходом. Даны рекомендации по выбору между глобальным и многоступенчатым кинетическими механизмами в сочетании с различными подходами к описанию турбулентного горения. Описано влияние модели переменных турбулентных чисел Прандтля и Шмидта на моделирование данного течения. Сформулированы идеи для дальнейших исследований.



О свойствах решения контактных задач с трением для штампа в виде четверти плоскости, контактирующего со слоистым основанием
Аннотация
В работе построено точное решение статической контактной задачи о действии с трением жесткого клиновидного штампа, занимающего первый квадрант, на слой изотропного композитного материала.
В отличие от многочисленных, в основном безуспешных, попыток решения этой и подобных задач аналитическими или численными методами, позволявшими выявлять лишь частичные свойства решения этой задачи, метод блочного элемента дал возможность вскрыть полную структуру ее решения. В работе доказано, что полученное решение точно удовлетворяет двумерному уравнению Винера–Хопфа для произвольной правой части. Расчет показателя особенности концентрации контактных напряжений в этой точке близок к значениям, рассчитанным численными методами в ряде работ. Доказано, что зона вблизи вершины штампа имеет превосходящую податливость при внедрении штампа в среду, по сравнению с удаленными зонами. Развитый метод применим для линейно упругих материалов и кристаллов, допускающих построение функции Грина и сводящихся к двумерным интегральным уравнениям Винера–Хопфа.
Установление общего вида решения рассмотренной контактной задачи открывает возможность изучения предвестников повышения сейсмичности в горных территориях, а также совершенствования численных методов для получения более точных решений усложненных контактных задач в инженерной практике.



Свободные колебания тонкой упругой ортотропной цилиндрической панели с шарнирно закрепленной граничной образующей
Аннотация
Используя систему уравнений соответствующей классической теории ортотропных цилиндрических оболочек, исследуются свободные колебания ортотропной тонкой упругой цилиндрической панели с шарнирно закрепленной граничной образующей. Для расчета собственных частот и идентификации соответствующих собственных мод используется обобщенный метод сведения к обыкновенным дифференциальным уравнениям Канторовича–Власова. Получены дисперсионные уравнения для нахождения собственных частот возможных типов колебаний. Установлена асимптотическая связь между дисперсионными уравнениями рассматриваемой задачи и аналогичной задачи для ортотропной прямоугольной пластины. Приводится механизм, с помощью которого расчленяются возможные типы краевых колебаний. На примере ортотропной цилиндрической панели получены приближенные значения безразмерной характеристики собственных частот колебаний.



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



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



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



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



Консервативная численная схема для решения уравнения Кана–Хилларда
Аннотация
В данной статье представлен консервативный численный алгоритм для решения уравнения Кана–Хилларда. Предложен способ линеаризации уравнения Кана–Хилларда, построена численная схема на основе метода контрольного объема. Подробно описана реализация предложенного численного алгоритма. Консервативность предложенной дискретной схемы проверена путем численного моделирования. Проведены численные эксперименты.


