Том 28, № 3 (2024)

Дифференциальные уравнения и математическая физика

Преобразование Халуты, осуществляемое с использованием различных операторов дробной производной

Khalouta A.

Аннотация

Недавно автором было введено и разработано новое интегральное преобразование, которое обобщает множество известных интегральных преобразований. Цель этой работы — расширение данного интегрального преобразования (преобразование Халуты) различными операторами дробной производной. Рассматриваются дробные производные в смысле Римана–Лиувилля, Лиувилля–Капуто, Капуто–Фабрицио, Атанганы–Балеану–Римана–Лиувилля и Атанганы–Балеану–Капуто. Доказаны теоремы, касающиеся свойств преобразования Халуты для решения дробных дифференциальных уравнений с использованием указанных операторов дробной производной. Приведено несколько примеров для проверки надежности и эффективности предложенной техники. Результаты показывают, что преобразование Халуты является эффективным инструментом при работе с дробными дифференциальными уравнениями.

Вестник Самарского государственного технического университета. Серия «Физико-математические науки». 2024;28(3):407-425
pages 407-425 views

Приближение решения уравнения переноса-диффузии в пространстве Гёльдера

Немдили А., Кориши Ф., Фуджита Яшима Х.

Аннотация

Рассматриваются приближенные решения для уравнения переноса-диффузии и их предельная функция, изучается принадлежность предельной функции к пространству Гёльдера, соответствующему регулярности данных. Цель исследования состоит в том, чтобы построить такое приближение решения уравнения переноса-диффузии, чтобы его основное свойство не зависело от величины коэффициента диффузии.
Точнее, рассматривается уравнение переноса-диффузии с постоянным коэффициентом диффузии в целом пространстве $ \mathbb{R}^d $ со свободным членом, который может зависеть от искомой функции. Приближенные решения на каждом шаге дискретизации по времени строятся с использованием ядра теплопроводности и локально  линеаризованного перемещения, соответствующего переносу. Приближенные решения оцениваются  в предположении, что заданные функции и их производные по $ x \in \mathbb{R}^d $ до порядка $m$ включительно ($ m \geqslant 2 $) равномерно ограничены на $[0, \tau]\times \mathbb{R}^d $ для каждого $\tau > 0$ и их производные порядка $m$ непрерывны по Гёльдеру с показателем $ \alpha $, $ {2}/{3} < \alpha \leqslant 1$. Оценки не зависят от величины коэффициента диффузии. На основании этих оценок доказываются равномерная сходимость приближенных решений и их производных по $ x $ до порядка $ m $ включительно на $[0, \tau]\times \mathbb{R}^d $, сходимость их производных порядка $ m $ в пространстве  Гёльдера $ C^{0+\alpha '} (\mathbb{R}^d)$, $ 0 < \alpha ' < \alpha$, для каждого $ t \geqslant 0 $ и непрерывность по Гёльдеру с показателем $\alpha$ производных по $ x $ порядка $ m $ предельной функции, которая удовлетворяет уравнению. То есть показано, что при использовании пространства Гёльдера получается та же дифференцируемость предельной функции, как и дифференцируемость данных, а в предыдущих работах для получения дифференцируемости порядка $ m $ предельной функции предполагалась дифференцируемость порядка $ m +1 $ данных.

Вестник Самарского государственного технического университета. Серия «Физико-математические науки». 2024;28(3):426-444
pages 426-444 views

Механика деформируемого твердого тела

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

Мурашкин Е.В., Радаев Ю.Н.

Аннотация

Рассматриваются вопросы распространения плоских гармонических связанных волн температурного инкремента, трансляционных и спинорных перемещений в полуизотропном термоупругом теле. Получены и проанализированы характеристические уравнения для волновых чисел плоских гармонических связанных термоупругих продольных (бикубическое уравнение) и поперечных волн перемещений (уравнение 8-й степени, естественным образом распадающееся на два алгебраических уравнения 4-й степени). Для продольной волны комплексные амплитуды температурного инкремента, трансляционных и спинорных перемещений оказываются также связанными, в отличие от поперечной волны. С помощью системы символьных вычислений Wolfram Mathematica 13 для волновых чисел поперечных волн получены алгебраические формы, содержащие многозначные комплексные квадратные и кубические радикалы.

Вестник Самарского государственного технического университета. Серия «Физико-математические науки». 2024;28(3):445-461
pages 445-461 views

Упругопластический анализ вращающегося полого цилиндра, жестко посаженного на вал, при наличии температурного градиента

Прокудин А.Н.

Аннотация

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

Вестник Самарского государственного технического университета. Серия «Физико-математические науки». 2024;28(3):462-488
pages 462-488 views

Анализ остаточных напряжений в поверхностно упрочненных вращающихся призматических элементах с полукруглыми надрезами в условиях высокотемпературной ползучести

Радченко В.П., Саушкин М.Н., Шишкин Д.М.

Аннотация

Разработан численный метод для расчета релаксации остаточных напряжений во вращающемся поверхностно упрочненном призматическом образце с полукруглым надрезом в условиях высокотемпературной ползучести. Задача моделирует напряженно-деформированное состояние образца, у которого одна грань закреплена на вращающемся с постоянной скоростью диске.
Методика включает в себя следующие этапы:
– реконструкцию полей остаточных напряжений и пластических деформаций после предварительного поверхностного пластического деформирования;
– расчет релаксации остаточных напряжений в процессе ползучести вращающегося призматического стержня.
Выполнен детальный анализ задачи для призматического образца размером 150 $\times$ 10 $\times$ 10 мм из сплава ЭП742. Одна из граней этого образца была упрочнена с помощью механического ультразвукового упрочнения. Решение задачи рассматривалось для образца с полукруглыми надрезами радиусами 0.1 и 0.3 мм, расположенными на расстоянии 2 и 75 мм от жестко закрепленной грани.
В области надрезов после опережающего поверхностного пластического деформирования задачи решались как в упругой, так и в упруго-пластической постановках. Полученные начальные поля остаточных напряжений и пластических деформаций использовались в качестве исходных данных для решения задачи ползучести.
Анализ влияния радиусов надрезов, их расположения, угловой скорости и начальных полей остаточных напряжений на релаксацию остаточных напряжений проводился при температуре 650 °C на основе феноменологической теории течения, построенной на известных экспериментальных данных для данного сплава.
Результаты показали, что для определения начального напряженно-деформированного состояния после предварительного пластического деформирования для надреза радиусом 0.1 мм необходимо использовать упруго-пластическое решение, тогда как для радиуса 0.3 мм различия между упругими и упруго-пластическими решениями незначительны.
Исследование релаксации остаточных напряжений проводилось при угловых скоростях 1500 и 2000 об/мин в течение 100 часов. Несмотря на заметную релаксацию остаточных напряжений для образцов с надрезами радиусом 0.1 и 0.3 мм, после полной температурно-силовой разгрузки в области надрезов все же наблюдается значительный уровень остаточных сжимающих напряжений. Это свидетельствует о высокой эффективности поверхностного упрочнения при эксплуатации в условиях высокотемпературной ползучести.

Вестник Самарского государственного технического университета. Серия «Физико-математические науки». 2024;28(3):489-515
pages 489-515 views

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

Роговой А.А.

Аннотация

В статье на примере сплава с памятью формы Cu–Al–Ni обоснован подход к описанию на микроструктурном уровне процессов двойникования и раздвойникования мартенситной фазы. Согласованная двойниковая мартенситная структура описывается уравнением совместности Адамара для деформаций, решение которого позволило определить поверхности, вдоль которых происходит сдвиг, направления и величину скольжения в орторомбической кристаллической ячейке, соответствующей рассматриваемому материалу в мартенситном состоянии, приводящих к появлению и исчезновению двойниковой структуры. Показано, что в сплаве с памятью формы одновременно и неразрывно существуют двойники двух видов: деформационный и структурный. Первый связан с деформацией простого сдвига, возникающей в соответствии с условием совместности Адамара в мартенситной пластине, что приводит к излому этой прямой пластины и возникновению двух элементов, повернутых на определенный угол друг относительно друга, которые и формируют этот двойник. Структурный двойник формируется из двух частей, в каждой из которых орторомбические кристаллические ячейки мартенсита одинаково ориентированы, но одна из осей этих ячеек меняет свое направление на 90° при переходе в другую часть двойника. Формирование деформационного двойника инициирует возникновение в среде деформации простого сдвига, а структурного двойника — структурной деформации. Показано, что определенное положение структурного двойника в деформационном приводит к равенству этих деформаций.

Вестник Самарского государственного технического университета. Серия «Физико-математические науки». 2024;28(3):516-542
pages 516-542 views

Способ определения параметров электрического сигнала для управления вынужденными установившимися колебаниями электровязкоупругих тел. Приложение к активному демпфированию колебаний

Севодина Н.В., Юрлова Н.А., Ошмарин Д.А.

Аннотация

При реализации активной стратегии управления динамическим поведением конструкций, в состав которых входят элементы, выполненные из пьезоэлектрических материалов, как правило, используют два пьезоэлемента, один из которых выступает в роли сенсора, а другой – актуатора. При этом проблема заключается в определении величины управляющего сигнала, подаваемого на актуатор, и аппаратной реализации необходимого закона управления. В связи с необходимостью формирования сложных электрических цепей, представляющих собой блок управления, привлекательным становится предварительное моделирование механического отклика на тот или иной управляющий сигнал.
В настоящей работе подход, позволяющий на основе решения задачи о собственных колебаниях электровязкоупругой конструкции получить аналитические выражения для определения величины электрического потенциала, генерируемого в момент резонанса на электродированной поверхности пьезоэлемента при его деформировании на рассматриваемой моде при вынужденных установившихся колебаниях, распространен на случай использования двух пьезоэлементов, выполняющих функции сенсора и актуатора и располагающихся соответствующим образом на поверхности конструкции.
Выведены аналитические выражения для определения величины управляющего сигнала, который подается на актуатор и обеспечивает демпфирование заданной моды колебаний. Управляющий сигнал формируется путем преобразования сигнала, получаемого от сенсора.
Приемлемость предложенного подхода продемонстрирована на примере консольно защемленной пластинки, выполненной из вязкоупругого материала, механическое поведение которого описывается комплексными динамическими модулями. По обе стороны пластинки размещены пьезоэлементы, играющие роли сенсора и актуатора. Численная реализация предложенного подхода осуществляется методом конечных элементов с использованием пакета прикладных программ ANSYS. Продемонстрировано хорошее совпадение результатов, полученных по выведенным формулам, с результатами расчета в ANSYS. Предложенный подход позволяет существенно сократить временные и ресурсные затраты при математическом моделировании активного управления вынужденными установившимися колебаниями электровязкоупругих тел, определить условия, которым должны удовлетворять элементы блока управления при реализации активной стратегии управления динамическим поведением такого рода smart-систем.

Вестник Самарского государственного технического университета. Серия «Физико-математические науки». 2024;28(3):543-561
pages 543-561 views

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

Янковский А.П.

Аннотация

Сформулирована задача термовязкоупругопластического деформирования армированных пологих оболочек при динамическом нагружении. При этом использована уточненная теория их изгиба, простейшим вариантом которой является теория Амбарцумяна. Геометрическая нелинейность учитывается в приближении Кармана. Неупругое поведение материалов фаз композиции описывается уравнениями теории пластического течения; вязкоупругое их деформирование описывается соотношениями модели Максвелла–Больцмана. Температура в трансверсальном направлении аппроксимирована полиномами высоких порядков. Численное интегрирование связанной нелинейной термомеханической задачи осуществляется с использованием явной схемы шагов по времени. Исследовано вязкоупругопластическое изгибное поведение цилиндрической панели из стеклопластика с ортогональной 2D-структурой армирования. Конструкция кратковременно нагружается в поперечном направлении давлением высокой интенсивности. Проведен сравнительный анализ расчетов, выполненных с учетом и без учета температурного отклика в пологой оболочке. Дополнительно исследованы случаи предварительного нагрева панели со стороны одной из лицевых поверхностей. Показано, что для расчета термовязкоупругопластического поведения стеклопластиковых искривленных панелей следует применять уточненную теорию изгиба, а не традиционно используемую теорию Амбарцумяна, которая существенно искажает форму расчетного остаточного прогиба и поля остаточных деформаций компонентов композиции. Продемонстрировано, что неучет температурного отклика в пологой армированной оболочке может привести не только к количественно, но и качественно неверному представлению о расчетной форме ее остаточного прогиба. Наличие предварительного нагрева стеклопластиковой панели приводит к заметному изменению ее остаточного прогиба. При этом важную роль играет то обстоятельство, со стороны какой именно лицевой поверхности осуществляется дополнительный нагрев или охлаждение тонкостенной конструкции.

Вестник Самарского государственного технического университета. Серия «Физико-математические науки». 2024;28(3):562-585
pages 562-585 views

Математическое моделирование, численные методы и комплексы программ

Разработка и сравнительный анализ математических моделей функционирования региональной энергосистемы Самарской области

Зотеев В.Е., Сагитова Л.А., Гаврилова А.А.

Аннотация

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

Вестник Самарского государственного технического университета. Серия «Физико-математические науки». 2024;28(3):586-608
pages 586-608 views