Том 27, № 3 (2023)
- Год: 2023
- Выпуск опубликован: 15.10.2023
- Статей: 11
- URL: https://journals.eco-vector.com/1991-8615/issue/view/8051
- DOI: https://doi.org/10.14498/vsgtu/v227/i3
Весь выпуск
Дифференциальные уравнения и математическая физика
Устойчивость и сходимость локально-одномерной схемы А. А. Самарского, аппроксимирующей многомерное интегро-дифференциальное уравнение конвекции-диффузии с неоднородными граничными условиями первого рода
Аннотация
Изучена первая начально-краевая задача для многомерного (по пространственным переменным) интегро-дифференциального уравнения конвекции-диффузии. Для приближенного решения поставленной задачи предложена локально-одномерная схема А. А. Самарского с порядком аппроксимации O(h2+τ). Исследование единственности и устойчивости решения проводится с помощью метода энергетических неравенств. Получены априорные оценки решения локально-одномерной разностной схемы, откуда следуют единственность решения, непрерывная и равномерная зависимость решения от входных данных, а также сходимость решения схемы к решению исходной дифференциальной задачи со скоростью, равной порядку аппроксимации разностной схемы. Для двумерной задачи построен алгоритм численного решения, проведены численные расчеты тестовых примеров, иллюстрирующие полученные в работе теоретические результаты.
![pages](/img/style/pages.png)
![views](/img/style/views.png)
![](/img/style/loadingSmall.gif)
Описание спектра одной операторной матрицы четвертого порядка
Аннотация
Рассматривается операторная матрица четвертого порядка A. Этот оператор соответствует гамильтониану системы с несохраняющимся числом и не более четырех частиц на решетке. Показано, что операторная матрица A унитарно эквивалентна диагональной матрице, диагональными элементами которой являются опять операторные матрицы четвертого порядка. Описано местоположение существенного спектра оператора A, т.е. выделены двухчастичная, трехчастичная и четырехчастичная ветви существенного спектра оператора A. Установлено, что существенный спектр операторной матрицы A состоит из объединения отрезков, число которых не больше 14. Построен определитель Фредгольма, такой, что его множество нулей и дискретный спектр операторной матрицы A совпадают, кроме того, доказано, что число простых собственных значений операторной матрицы A, лежащих вне существенного спектра, не превосходит 16.
![pages](/img/style/pages.png)
![views](/img/style/views.png)
![](/img/style/loadingSmall.gif)
О разрешимости одного класса нелинейных двумерных интегральных уравнений типа Гаммерштейна–Немыцкого на плоскости
Аннотация
Работа посвящена изучению вопросов существования и единственности положительного ограниченного и непрерывного решения для одного класса двумерных нелинейных интегральных уравнений с некомпактным оператором Гаммерштейна–Немыцкого на плоскости. Такие уравнения возникают в теории p-адических открытых и открыто-замкнутых струн, в кинетической теории газов, в математической теории географического распространения эпидемических заболеваний. Доказываются конструктивные теоремы существования и единственности ограниченного положительного решения. Исследуется также асимптотическое поведение построенного решения на бесконечности. Приводятся конкретные прикладные примеры указанного класса уравнений.
![pages](/img/style/pages.png)
![views](/img/style/views.png)
![](/img/style/loadingSmall.gif)
Механика деформируемого твердого тела
Упругая составная плоскость с частично оторванным от матрицы межфазным абсолютно жестким тонким включением с учетом проскальзывания на концах
Аннотация
Рассмотрено плоско-деформированное состояние базовой плоскости упругого составного пространства с трещиной конечной длины на линии соединения составляющих полуплоскостей. В один из берегов межфазной трещины под действием сосредоточенной силы вдавливается абсолютно жесткое тонкое включение такой же длины. Для контактирующей стороны включения полагается, что в средней ее части имеет место сцепление с матрицей, а по краям происходит проскальзывание, описываемое законом сухого трения. Задача сформулирована в виде системы сингулярных интегральных уравнений. Исследовано поведение искомых функций в окрестности концов включения-трещины и в точках раздела зон сцепления и проскальзывания. Определяющая система интегральных уравнений решается методом механических квадратур. Найдены законы распределения контактных напряжений, а также длины зон сцепления и проскальзывания в зависимости от коэффициента трения, коэффициентов Пуассона и отношения модулей Юнга материалов полуплоскостей, а также угла наклона внешней силы.
![pages](/img/style/pages.png)
![views](/img/style/views.png)
![](/img/style/loadingSmall.gif)
Влияние поверхностного пластического упрочнения на геометрические параметры круговых концентраторов напряжений в пластинах
Аннотация
Разработана методика изучения влияния упрочняющей обработки на форму концентраторов напряжений в виде сквозных круговых отверстий в пластинах после процедуры поверхностно-пластического деформирования.
Рассмотрены две модельных задачи:
– определение геометрической конфигурации кругового концентратора напряжения, вырезанного в прямоугольной пластине, подвергшейся опережающему поверхностно-пластическому деформированию;
– определение геометрической конфигурации кругового концентратора напряжения в круговой цилиндрической пластине, поверхность которого подверглась поверхностно-пластическому деформированию.
Приведены феноменологические методы восстановления полей остаточных напряжений и пластических деформаций в пластинах после процедуры упрочнения. Краевые задачи реконструкции напряженно-деформированного состояния сведены к корректным задачам фиктивной термоупругости. На модельных расчетах для прямоугольной пластины из сплава ЭП742 и круговой цилиндрической пластины из сплава ЭИ698 проиллюстрирована адекватность предлагаемых подходов.
Получены профили образующих концентраторов напряжений плит. В случае опережающего поверхностного пластического деформирования верхней грани квадратной шарнирно опертой пластины толщиной 10 мм максимальное смещение образующей относительно первоначальной конфигурации составило около 4 мкм. Показано, что с уменьшением толщины пластины максимальное смещение образующей убывает. В случае упрочнения поверхности кругового концентратора напряжений цилиндрической пластины максимальное смещение образующей концентратора напряжений составило около 1.4 мкм для пластин, опертых шарнирно и с жесткой заделкой боковой грани. Показано, что с уменьшением радиуса отверстия смещение образующей возрастает.
![pages](/img/style/pages.png)
![views](/img/style/views.png)
![](/img/style/loadingSmall.gif)
Численное решение задачи о напряженно-деформированном состоянии поверхностно упрочненного призматического образца с надрезом V-образного профиля в упругой и упругопластической постановках
Аннотация
Разработан метод решения задачи расчета напряженно-деформированного состояния в поверхностно упрочненном образце со сквозным поперечным надрезом V-образного профиля при различных значениях угла раскрытия в упругой и упругопластической постановках. Метод базируется на конечно-элементном моделировании и известном начальном напряженно-деформированном состоянии для гладкого упрочненного образца. Выполнено детальное исследование влияния угла раскрытия надреза и его глубины на уровень и характер распределения остаточных напряжений от дна концентратора напряжений по толщине упрочненного слоя для обеих постановок задач. На основании данных расчета обоснована целесообразность исследования поставленной задачи в упругопластической постановке, когда надрез находится полностью или частично в упрочненном слое, так как величины остаточных напряжений при решении задачи в упругой постановке физически нереализуемы, поскольку их значения превосходят по модулю временной предел сопротивления материала в несколько раз.
В этом случае погрешность между решениями в упругой и упругопластической постановках для остаточных напряжений в среднеквадратической норме достигает 100–200 %, а при равномерной оценке (норма Чебышева) — нескольких сотен процентов. Если глубина концентратора превышает величину упрочненного слоя более чем в 1.5 раза, то упругое и упругопластическое решения дают близкие результаты.
![pages](/img/style/pages.png)
![views](/img/style/views.png)
![](/img/style/loadingSmall.gif)
Параметрическое исследование полей, ассоциированных с вершиной трещины, в условиях ползучести с учетом процессов накопления поврежденности с использованием UMAT
Аннотация
Предметом исследования является анализ полей, ассоциированных с вершиной трещины, находящейся в условиях ползучести при принятии во внимание явления накопления повреждений. Целью работы является проведение компьютерного конечно-элементного моделирования одноосного растяжения пластины с центральной горизонтальной и наклонной трещинами в условиях ползучести в плоской постановке задачи и анализ поля сплошности вблизи вершины трещины. При численном моделировании используется степенной закон ползучести Бейли–Нортона. Моделирование выполнено в многофункциональном программном комплексе SIMULIA Abaqus. Проведен анализ окружных распределений напряжений и деформаций ползучести в окрестности вершины трещины.
Cтепенной закон ползучести с помощью пользовательской процедуры UMAT (User Material) пакета SIMULIA Abaqus был дополнен кинетическим уравнением накопления поврежденности Качанова–Работнова в связанной постановке. Примененная подпрограмма UMAT имеет много преимуществ при прогнозировании поврежденности материала и позволяет работать с материалами и определяющими их соотношениями, отсутствующими в библиотеке материалов Abaqus. Подпрограмма UMAT вызывается во всех точках расчета и обновляет напряжения и переменные состояния, зависящие от решения, до их значений в конце приращения. После чего рассчитываются обновленные элементы матрицы Якоби.
Получены распределения напряжений, деформаций и сплошности в условиях ползучести с учетом накопления поврежденности с течением времени. Построены угловые распределения сплошности, напряжений и деформаций с течением времени на различных расстояниях от вершины трещины с применением библиотеки Matplotlib. Проведено сравнение угловых распределений напряжений и деформаций при моделировании без учета поврежденности и в случае учета накопления повреждений. Показано, что наличие поврежденности приводит к большим значениям деформаций ползучести и меньшим значениям напряжений.
![pages](/img/style/pages.png)
![views](/img/style/views.png)
![](/img/style/loadingSmall.gif)
Математическое моделирование, численные методы и комплексы программ
Неоднородное течение Куэтта двухслойной жидкости
Аннотация
Предложено новое точное решение уравнений Навье–Стокса, описывающее установившееся изобарическое изотермическое течение стратифицированной по плотности и/или вязкости несжимаемой двуслойной жидкости. Указанное точное решение принадлежит классу функций, линейных по части пространственных координат, и является обобщением классического течения Куэтта в протяженном горизонтальном слое на случай неодномерных неоднородных течений. В качестве системы краевых условий рассмотрена связка «условие прилипания + воздействие параболического ветра». На общей границе двух слоев заявлено выполнение требования гладкости и непрерывности решения. Построенное для каждого слоя решение было исследовано на предмет возможности описывать возникновение застойных точек поля скорости и генерации противотечений. Строго показано, что указанное решение при определенном граничном управлении и варьировании геометрико-физических характеристик слоя отвечает множественной стратификации как поля скорости, так и порождаемого им поля касательных напряжений.
![pages](/img/style/pages.png)
![views](/img/style/views.png)
![](/img/style/loadingSmall.gif)
Идентификация передаточной функции посредством минимизации рассогласования оценок состояния адаптивного и оптимального фильтров
Аннотация
Статья посвящена дальнейшему развитию активного принципа параметрической идентификации системы в классе линейных, инвариантных во времени, полностью наблюдаемых моделей. В качестве целевой модели идентификации выбран оптимальный фильтр Калмана (ОФК), который не более чем концептуально присутствует в дискретно наблюдаемом отклике системы на обучающее возбуждение типа белого шума. Путем модификации физически заданной структуры в стандартную наблюдаемую модель как в наблюдаемом отклике, так и в адаптивном фильтре Калмана (АФК), строится так называемый обобщенный остаток (ОО), равный рассогласованию между оценками состояния адаптивного и оптимального фильтров плюс независимая от АФК шумовая составляющая. В силу этой модификации средний квадрат ОО становится новым критерием близости модели для этих фильтров. Минимизация этого критерия с помощью обычных практических методов оптимизации дает точно такой же результат (АФК = ОФК), как и минимизация теоретического критерия, который, к сожалению, недостижим для любых методов численной оптимизации АФК. В статье представлена подробная пошаговая процедура, объясняющая вышеуказанное решение в терминах параметризованной передаточной функции. Для наглядности и стимулирования применения подхода в реальном мире в статье используется модель передаточной функции линии витой пары в типичной системе xDSL. Обсуждаются проблемы реализации теоретических положений метода. Вопрос о распространении предложенного подхода на проблемы идентификации линейных моделей для нелинейных систем обозначен в направлениях дальнейших исследований.
![pages](/img/style/pages.png)
![views](/img/style/views.png)
![](/img/style/loadingSmall.gif)
Стереопанорамный анеморумбометр для системы ориентации горизонтально-осевой ветроэнергетической установки
Аннотация
Традиционный подход к ориентации ротора горизонтально-осевой ветроэнергетической установки по ветру приводит к появлению известной дифференциальной ошибки ориентации из-за вращающихся лопастей и периодического отклонения воздушного потока. Для снижения ее величины в традиционном подходе используется флюгер, расположенный сверху гондолы.
В настоящем исследовании предлагается новый подход — использование комплексного или «стереодатчика» в виде двух устройств, симметрично расположенных по обе стороны гондолы (аналогично стереоскопическим устройствам). Для доказательства эффективности подхода были выбраны несколько характерных точек вблизи гондолы для последующего моделирования воздушных потоков в ее области в программе ANSYS® CFX с использованием k–ε модели турбулентности на основе дифференциальных уравнений Навье–Стокса.
В каждой точке была рассчитана средняя величина ошибки угла ориентации при следующих условиях: различных скоростях ветра, значениях быстроходности и углов направления на ветер. В результате выявлены две наиболее подходящие для размещения приборов точки. В численном виде показано преимущество стереопанорамного анеморумбометра перед традиционным на примере расчетного случая с номинальными параметрами. Анализ в Matlab/Simulink показал прирост производительности ветроэнергетической установки за счет повышения достоверности определения направления ветра при применении правильно расположенных датчиков ветрового потока.
Данная статья не дает представления о конструкции датчика, поскольку для определения правильного направления ветра можно использовать любой принцип. Однако авторами рассматривается новый «стереодатчик», который будет более детально исследоваться в следующих работах.
![pages](/img/style/pages.png)
![views](/img/style/views.png)
![](/img/style/loadingSmall.gif)
Краткие сообщения
Исследование напряженно-деформированного состояния упругоподкрепленной сжатой полосы
Аннотация
Проведен анализ непрерывной зависимости функции, описывающей поведение реальной конструкции, от характеристик начальных несовершенств. Получено условие, накладываемое на параметр внешнего воздействия и коэффициент жесткости основания, при нарушении которого форма поперечного сечения полосы уже не будет близкой к прямоугольнику, т.е. полоса теряет устойчивость формы. При проведении исследования параметры внешних воздействий оставались независимыми.
Исходный вариант статьи был опубликован в сборнике Актуальные проблемы прикладной математики, информатики и механики. Воронеж, 2022. С. 1265–1269
![pages](/img/style/pages.png)
![views](/img/style/views.png)
![](/img/style/loadingSmall.gif)