Все выпуски

В статье рассматривается модель антропоморфного механизма типа экзоскелета со звеньями переменной длины. Комплексно рассмотрены четыре модели звеньев переменной длины: модель звена экзоскелета переменной длины с упругим элементом и абсолютно твердым весомым стержнем, модель телескопического звена; модель звена с массами в шарнирах-суставах и между ними, модель звена с произвольным количеством масс. Составлены дифференциальные уравнения движения в форме уравнений Лагранжа второго рода. На основе проведенного анализа дифференциальных уравнений движения для многозвенных стержневых механических систем типа экзосклета выявлена их структура, позволившая представить их в векторно-матричном виде. Впервые установлены общие закономерности построения матриц и получены обобщения выражений для элементов матриц в двухмерном случае. Приводятся новые рекуррентный и матричный методы составления дифференциальных уравнений движения. Предлагается единый подход к построению дифференциальных уравнений движения экзоскелета на основе разработанных рекуррентного и матричного методов записи дифференциальных уравнений движения экзоскелета. Проведено сопоставление времени составления дифференциальных уравнений движения предложенными методами, в сравнении с уравнениями Лагранжа второго рода, в системе компьютерной математики Mathematica. Осуществлено аналитическое исследование модели экзоскелета. Установлено, что для механизмов с $n$ подвижными звеньями решение задачи Коши для систем дифференциальных уравнений движения при любых начальных условиях существует, единственно и неограниченно продолжаемо. Управление экзоскелетом осуществляется с помощью крутящих моментов, расположенных в шарнирах-суставах в местах соединения звеньев и моделирующих управляющие воздействия. Выполнено численное исследование модели экзоскелета, проведено сопоставление результатов расчетов для экзоскелетов с различными моделями звеньев. Для численного исследования использованы эмпирические данные о человеке и его движениях. Установлено, что при выборе конструкции экзоскелета модель с сосредоточенными массами является предпочтительной, нежели модель с абсолютно твердым весомым стержнем, так как экзоскелет, обеспечивающий комфортабельные передвижения человека в нем, должен повторять свойства опорно-двигательного аппарата.

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

В настоящей работе предложен вариант импорта в систему ANSYS Mechanical APDL модели поведения гибких тканых композиционных материалов с армирующей тканью полотняного переплетения при статическом растяжении вдоль нитей армирования. Импорт осуществлен при помощи использования, разработанного авторами и представленного в текущей работе модуля интеграции, основанного на использовании аналитической модели деформирования исследуемого материала, представленной в опубликованных ранее статьях и учитывающей изменения геометрической структуры, происходящие в армирующем слое материала в процессе деформирования, образование необратимых деформаций и взаимодействие накрест лежащих нитей армирующей ткани. Во введении кратко описаны основные вводимые характеристики полотняного переплетения армирующей ткани и аналитической модели, импортируемой в ANSYS. Аналитическая модель основана на описании процессов деформирования элементарной периодической ячейки исследуемого класса материалов. Входными параметрами модуля являются механические характеристики материалов, входящих в состав композита (связующее и материал нитей армирования), геометрические характеристики переплетения армирующей ткани. Алгоритм импорта модели основан на вычислении и передачи в ANSYS расчетных точек диаграммы деформирования материала при одноосном растяжения вдоль направления армирования и использовании вложенных в систему ANSYS пользовательской моделей материала Multilinear Kinematich Hardening. Аналитическая модель, импортируемая при помощи представленного модуля, позволяет моделировать композиционный материал с армирующей тканью без детального описания геометрии переплетения нитей при моделировании материала в целом. Выполнена верификация импортированной модели. Для верификации были проведены натурные экспериментальные исследования и численное моделирование растяжения образцов из гибких тканых композитов. В экспериментах принимали участие образцы трех марок: VP4126, VP6131 и VP6545. Погрешность импортированной в ANSYS модели составила менее 10 % относительно экспериментальных исследований для всех марок материала. Анализ полученных результатов показал хорошее качественное и количественное согласование расчетов в системе ANSYS с применением импортированной модели и натурных испытаний до величин предельных деформаций, соответствующих разрушению образцов материала для всех исследуемых марок гибких тканых композитов, что позволяет сделать вывод о возможности применения предложенного модуля при моделировании процессов деформирования гибких тканых композитов и конструкций, созданных из таких материалов при статическом одноосном растяжении вдоль нитей армирования.

Гибкие тканые композиты относят к классу высокотехнологичных инновационных материалов. Благодаря совмещению различных компонентов наполнителя и элементов армирования такие материалы применяют в строительстве, оборонной промышленности, судо- и авиастроении и др. В отечественной литературе уделено недостаточное внимание к тканым композитам, изменяющим свою геометрическую структуру армирующего слоя в процессе деформирования. В настоящей работе приводится анализ предложенного ранее комплексного подхода к моделированию поведения гибких тканых композитов при статическом одноосном растяжении для дальнейшего обобщения подхода на двухосное растяжение. Работа нацелена на качественное и количественное описание механических деформационных процессов, протекающих в структуре исследуемых материалов при растяжении, к которым относится распрямления нитей армирующего слоя и увеличение величины взаимного надавливания накрест лежащих нитей армирования. В начале процесса деформирования распрямление нитей и увеличение взаимного надавливания нитей наиболее интенсивны. С увеличением уровня нагрузки изменение указанных параметров замедляется. Например, изгиб нитей армирования переходит в центральное растяжение, а величина нагрузки от взаимного надавливания более не увеличивается (стремится к константе). Для моделирования описанных процессов вводятся основные геометрические и механические параметры материала, влияющие на процесс формоизменения, приводятся необходимая терминология и описание характеристик. В связи с высокой геометрической нелинейностью все процессы описаны в приращениях, так как на начальных значениях нагрузки происходит значительное формоизменение армирующего слоя. Для количественного и качественного описания механических деформационных процессов, протекающих в армирующем слое, выведены аналитические зависимости, позволяющие определить приращение угла распрямления нитей армирования и нагрузки, вызванной взаимным надавливанием накрест лежащих нитей на каждом шаге приращения нагрузки. Для апробации выведенных зависимостей приведен пример их применения для гибких тканых композиционных материалов марок VP4126, VP6131 и VP6545. Результаты моделирования подтвердили предположения о процессах выпрямления нитей и замедления увеличения взаимного надавливания нитей. Приведенные в данной работе результаты и зависимости имеют непосредственное отношение к дальнейшему обобщению предложенных ранее аналитических моделей для двухосного растяжения, так как растяжение в двух направлениях существенно уменьшит выпрямление нитей и увеличит величину взаимного надавливания при аналогичных нагрузках.

В данной статье рассматривается метод построения прогнозной нейросетевой модели временного ряда, основанный на определении состава входных переменных, построения обучающей выборки и самого обучения с использованием метода обратного распространения ошибки. Традиционные методы построения прогнозных моделей временного ряда (авторегрессионной модели, модели скользящего среднего или модели авторегрессии – скользящего среднего) позволяют аппроксимировать временной ряд линейной зависимостью текущего значения выходной переменной от некоторого количества ее предыдущих значений. Такое ограничение, как линейность зависимости, приводит к значительным ошибкам при прогнозировании.

Технологии интеллектуального анализа с применением нейросетевого моделирования позволяют аппроксимировать временной ряд нелинейной зависимостью. Причем процесс построения нейросетевой модели (определение состава входных переменных, числа слоев и количества нейронов в слоях, выбор функций активации нейронов, определение оптимальных значений весов связей нейронов) позволяет получить прогнозную модель в виде аналитической нелинейной зависимости.

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

В работе предлагается метод определения состава входных переменных нейросетевых моделей для стационарных и нестационарных временных рядов, базирующийся на построении и анализе автокорреляционных функций. На основе предложенного метода разработаны алгоритм и программа в среде программирования Python, определяющая состав входных переменных прогнозной нейросетевой модели — персептрона, а также строящая саму модель. Осуществлена экспериментальная апробация предложенного метода на примере построения прогнозной нейросетевой модели временного ряда, отражающего потребление электроэнергии в разных регионах США, открыто опубликованной компанией PJM Interconnection LLC (PJM) — региональной сетевой организацией в Соединенных Штатах. Данный временной ряд является нестационарным и характеризуется наличием как тренда, так и сезонности. Прогнозирование очередных значений временного ряда на ос- нове предыдущих значений и построенной нейросетевой модели показало высокую точность аппроксимации, что доказывает эффективность предлагаемого метода.

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

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

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

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

Одной из перспективных технологий повышения нефтеотдачи при разработке нетрадиционных нефтяных пластов является метод термогазового воздействия. Метод основан на закачке в пласт кислородосодержащей смеси и ее трансформации в высокоэффективный смешивающийся с пластовой нефтью вытесняющий агент за счет самопроизвольных внутрипластовых окислительных процессов. В ряде случаев этот метод обладает большим потенциалом по сравнению с другими способами повышения нефтеотдачи. В данной работе рассматриваются некоторые вопросы распространения волн внутрипластового горения. В зависимости от параметров коллектора и закачиваемой смеси такие волны могут распространяться в различных режимах. В данной работе рассматривается только прямоточный инверсный режим распространения. В этом режиме волна горения распространяется в направлении течения окислителя и фронт реакции отстает от тепловой волны, в которой вещество (углеводородные фракции, пористый скелет и др.) прогреваются до температур, достаточных для протекания реакции окисления. В работе представлены результаты аналитического исследования и численного моделирования структуры инверсной волны внутрипластового горения при двухфазном течении в пористом слое. Сделаны упрощающие предположения о теплофизических свойствах флюидных фаз, которые позволяют, с одной стороны, сделать модель внутрипластового горения обозримой для анализа, а с другой — передать основные особенности этого процесса. Рассмотрено решение типа «бегущая волна» и указаны условия его реализации. Выделено два режима распространения инверсных волн внутрипластового горения: гидродинамический и кинетический. Численное моделирование распространения волны внутрипластового горения проводилось с помощью термогидродинамического симулятора, разработанного для численного интегрирования неизотермических многокомпонентных фильтрационных течений, сопровождающихся фазовыми переходами и химическими реакциями.

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

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