


№ 1 (2024)
Статьи
Моделирование изгибно-гравитационных волн в ледяном покрове на упругих пленках
Аннотация
Отмечено, что на основании наблюдений за колебаниями ледяного покрова в натурных условиях под действием движущихся нагрузок, т.е. при возбуждении изгибно-гравитационных волн (ИГВ), последний ведет себя аналогично упругой изотропной пластине. На этом основании предложено новое направление в моделировании некоторых задач деформирования ледяного покрова ИГВ на упругих пленках в обычных опытовых бассейнах. Возможность такой технологии подтверждена результатами сопоставления записей деформирования движущимися нагрузками упругого модельного слоя и натурного ледяного покрова. На основании теории подобия и размерностей получены зависимости для пересчета данных модельных испытаний на натуру. Отмечается, что затраты на проведение подобных модельных экспериментов несоизмеримо меньше с затратами при проведении опытов в ледовых бассейнах. Перечисляются ледотехнические задачи, при решении которых может быть использована разработанная методика моделирования ИГВ.



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



Совершенствование метода фазовой плоскости для исследования влияния эффекта “бифуркационной памяти” на динамику корабля
Аннотация
Изучение эффекта “бифуркационной памяти” играет важную роль в исследовании динамических особенностей реальных систем. Практический интерес заключается в изучении возможности прогнозирования временного снижения реакции на управление, что может существенно повысить безопасность навигации. Эффект “бифуркационной памяти” заключается во временном снижении фазовой скорости изображающей точки при прохождении определенной области (“фазового пятна”) на фазовой плоскости. “Фазовое пятно” возникает вблизи исчезнувшего при бифуркации состояния равновесия. За почти полувековую историю изучения этой динамической особенности, предложено очень мало методов, позволяющих однозначно и с достаточной точностью идентифицировать эффект “бифуркационной памяти”. В данной работе предложен усовершенствованный метод фазовой плоскости, заключающийся в построении годографа фазовой скорости. Отличительной чертой предлагаемого метода является не только то, что он охватывает обе фазовые координаты, а также дает адекватный результат при любых начальных условиях. Метод вполне универсален и может использоваться для изучения эффекта “бифуркационной памяти” в различных динамических системах. Информация ограничных значениях параметра – угла перекладки руля, при которых начинает (заканчивает) проявляться эффект “бифуркационной памяти” может быть использована, например, в задаче оптимизации конструкции корпуса и рулей или при создании алгоритмов управления.



Кватернионная регуляризация дифференциальных уравнений возмущенного центрального движения и регулярные модели орбитального (траекторного) движения: обзор и анализ моделей, их приложение
Аннотация
В обзорной статье кратко излагается предложенная нами общая кватернионная теория регуляризующих и стабилизирующих преобразований ньютоновских дифференциальных уравнений возмущенного движения материальной точки в центральном силовом поле, потенциал которого полагается произвольной дифференцируемой функцией расстояния от точки до центра поля. Точка находится также под действием возмущающего потенциала, полагаемого произвольной функцией времени и декартовых координат местоположения точки, и под действием возмущающего ускорения, полагаемого произвольной функцией времени, радиус-вектора и вектора скорости точки. Рассмотрены условия приводимости излагаемых кватернионных уравнений возмущенного центрального движения к осцилляторному виду с помощью использования трех регуляризующих функций, содержащих расстояние до центра поля. Приведены различные дифференциальные кватернионные уравнения возмущенного центрального движения в осцилляторной и нормальной формах, построенные с помощью этой теории, в том числе регулярные уравнения, в которых используются четырехмерные параметры Эйлера (Родрига–Гамильтона) или четырехмерные переменные Кустаанхеймо–Штифеля или их модификации, предложенные нами. Рассмотрены регулярные кватернионные уравнения пространственного невозмущенного центрального движения материальной точки, связи используемых четырехмерных переменных с элементами орбиты, униформизированное решение пространственной задачи невозмущенного центрального движения. В качестве приложения изложены регуляризованные дифференциальные кватернионные уравнения движения искусственного спутника в гравитационном поле Земли в четырехмерных переменных Кустаанхеймо–Штифеля, а также в наших модифицированных четырехмерных переменных и в параметрах Эйлера. Приведен анализ изложенных регулярных кватернионных уравнений возмущенного центрального движения, показывающий, что кватернионный метод регуляризации, основанный на использовании параметров Эйлера или переменных Кустаанхеймо–Штифеля или их модификаций уникален в совместной регуляризации, линеаризации и увеличении размерности для трехмерных кеплеровских систем и центрального движения.
Изложенные регуляризованные (в отношении ньютоновской силы притяжения) дифференциальные кватернионные уравнения движения искусственного спутника в гравитационном поле Земли в наших модифицированных четырехмерных переменных имеют перед кватернионными уравнениями в переменных Кустаанхеймо–Штифеля преимущества, указанные в статье. В изложенных дифференциальных кватернионных уравнениях движения спутника, построенных с использованием четырехмерных параметров Эйлера, регуляризуются слагаемые уравнений, содержащих отрицательные степени расстояния до центра Земли четвертого порядка включительно. Во всех этих регуляризованных уравнениях в описании гравитационного поля Земли учитываются не только центральная (ньютоновская), но и зональные, тессеральные и секториальные гармоники потенциала поля тяготения Земли (учитывается несферичность Земли).



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



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



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



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



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



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



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



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



Два способа управления cферическим роботом маятникового типа на подвижной платформе в задаче преследования
Аннотация
Рассматривается проблема управления сферическим роботом с маятниковым приводом, катающимся по платформе, которая способна двигаться поступательно в горизонтальной плоскости абсолютного пространства. На сферический робот наложены голономные и неголономные ограничения. Некоторая точечная цель движется на уровне геометрического центра сферического робота и не касается самой подвижной платформы. Программа движения, позволяющая сферическому роботу преследовать цель, задается посредством двух сервосвязей. Движение за целью робот может осуществлять из любого положения и с любыми начальными условиями. Предлагаются два способа управления данной системой в абсолютном пространстве: посредством управления вынужденным движением платформы (колебания маятника свободные) и посредством управления крутящим моментом маятника (платформа неподвижна или колеблется несогласованно со сферическим роботом). Построены уравнения движения системы. В случае свободных колебаний маятника система уравнений движения обладает первыми интегралами и при необходимости может быть редуцирована на фиксированный уровень этих интегралов. При движении сферического робота по прямой для редуцированной на уровень интегралов системы построены фазовые кривые, графики расстояния от геометрического центра сферического робота до цели, траектории движения выбранной точки платформы при управлении платформой, квадрат управляющего крутящего момента при управлении маятниковым приводом. При движении робота по криволинейной траектории интегрирование ведется в исходных переменных. Строятся графики квадратов угловой скорости маятника и самого сферического робота, а также траектории движения робота в абсолютном пространстве и на подвижной платформе. Численные эксперименты выполнялись в программном пакете Maple.



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



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


