Идея построения специального тензора для описания параметров структурно-неоднородных материалов возникла из целого ряда попыток количественно охарактеризовать микроструктуру упругого пористого материала. Использование специальных тензорных величин для описания стереометрических характеристик структурно-анизотропных материалов позволяет в компактном виде выразить значимые структурные параметры исследуемых объектов. Преимущественная ориентация пор внутри образца хорошо описывается тензором структуры и, алгебраически связанным с ним тензором анизотропии. Приведенные в работе математические выкладки, позволили формализовать процесс вычисления всех необходимых для построения тензора структуры параметров. Алгоритмизация метода определения среднего расстояния между порами, легла в основу разработанного специализированного программного обеспечения для расчета компонент тензора структуры. Верификация программного модуля была осуществлена путем проведения стереологического исследования ряда идеализированных тестовых структур и образца пористого материала, для которого тензор структуры был известен заранее. Полученные результаты не противоречили природной действительности, совпадали с ранее полученными данными и описывали степень анизотропии исследованных структур с высокой степенью точности. В качестве демонстрации практического использования разработанного программного комплекса в работе представлены результаты исследования образца трабекулярной костной ткани шейки бедренной кости человека и образца автоклавного газобетона. Проведены вычисления всех необходимых параметров и приведены изображения эллипса структуры исследованных пористых материалов. Из полученных результатов видно, что тензор структуры способен описывать стереометрические характеристики натуральных и искусственных пористых структур, а пакет проблемно-ориентированных программ позволяет автоматизировать процесс определения всех необходимых параметров.
Идентификаторы и классификаторы
- SCI
- Информатика
Существует большое количество структурно-неоднородных материалов как естественного (горные породы, трабекулярная костная ткань), так и искусственного происхождения (пеностекло, газобетон, высокопористый пироуглерод), физические свойства которых во многом определяются их внутренним строением [1–3]. Методы количественной стереологии широко применяются для описания структурных особенностей зернистых и пористых материалов [1; 4–11]. Идея построения специального тензора для описания параметров структурнонеоднородных материалов возникла из целого ряда попыток количественно охарактеризовать микроструктуру упругого пористого материала. В 1978 г. Cowin обратил внимание на то, что эллипсоид структуры эквивалентен тензору второго ранга, отмечая при этом, что для пористых материалов эллипсоид структуры является наилучшей мерой описания структурного распределения материала [1; 12]. В 1980 г. Oda, Konishi и Nemat-Nasser использовали концепцию эллипсоида структуры, аргументируя это тем, что для зернистых материалов эллипсоид структуры является второй наилучшей мерой описания структуры материала после пористости
Список литературы
1. Cowin, S.C. Bone mechanics handbook / S.C. Cowin. - Second edition. - New York: CRC Press. - 2001. - 1136 p.
2. Киченко, А.А. Математическое моделирование структуры губчатой костной ткани с использованием тензора структуры / А.А. Киченко, В.М. Тверье, А.В. Сотин // Математические методы в технологиях и технике. - 2023. - № 9. - С. 41-44. EDN: HXABXL
3. Няшин, Ю.И. Экспериментальные методы в биомеханике / Ю.И. Няшин, Р.М. Подгаец. - Пермь: Изд-во Перм. гос. техн. ун-та, 2008. - 400 с.
4. Harrigan, T.P. Characterization of microstructural anisotropy in orthotropic materials using a second rank tensor / T.P. Harrigan, R.W. Mann //j. Mater. Sci. - 1984. - Vol. 19. - P. 761-767. EDN: RGUBXZ
5. Telega, J.J. Fabric tensor in bone mechanics /j.J. Telega, S. Jemiolo //j. Engineering Transactions. - 1998. - Vol. 46. - P. 3-26.
6. Underwood, E. Quantitative stereology / E. Underwood. - Mass.: Addision Wesley, 1970. - 274 p.
7. Whitehouse, W.J. The quantitative morphology of anisotropic trabecular bone / W.J. Whitehouse //j. Microscopy. - 1974. - Vol. 101. - P. 153-168.
8. Whitehouse, W.J. A stereological method for calculating the internal surface areas in structures which have become anisotropic as the result of linear expansions or contractions / W.J. Whitehouse //j. Microscopy. - 1974. - Vol. 101. - P. 169-176.
9. Экспериментальное определение тензора структуры трабекулярной костной ткани / А.А. Киченко, В.М. Тверье, Ю.И. Няшин, А.А. Заборских // Российский журнал биомеханики. - 2011. - Т. 15, № 4. - С. 78-93. EDN: OOQWFD
10. О приложении теории перестройки трабекулярной костной ткани / А.А. Киченко, В.М. Тверье, Ю.И. Няшин, М.А. Осипенко, В.А. Лохов // Российский журнал биомеханики. - 2012. - Т. 16, № 4. - С. 53-72. EDN: PPMLBP
11. Experimental construction of the fabric tensor for trabecular bone tissue / V. Tverier, A. Kichenko, Y. Nyashin, V. Lokhov // Series on Biomechanics. - 2015. - Vol. 29, no. 4. - P. 33-38. EDN: XNBNIO
12. Martin, R.B. Skeletal tissue mechanics. Second edition / R.B. Martin, D.B. Burr, N.A. Sharkey. - New York: Springer-Verlag. - 1998. - 392 p.
13. Oda, M. Fabrics and their effects on the deformation behaviors of sand / M. Oda. - Saitama University: Dept. of Foundation Eng., 1976. - 59 p.
14. Cowin, S.C. Fabric dependence of an anisotropic strength criterion / S.C. Cowin //j. Mech. Materials. - 1986. - Vol. 5. - P. 251-260.
15. Cowin, S.C. Wolff’s law of trabecular architecture at remodeling equilibrium / S.C. Cowin //j. Biomech. Engng. - 1986. - Vol. 108. - P. 83-88.
16. The fabric dependence of the orthotropic elastic constants of cancellous bone / C.H. Turner, S.C. Cowen, J.Y. Rho, R.B. Ashman, J.C. Rice //j. Biomechanics. - 1990. - Vol. 23. - P. 549-561.
17. Cowin, S.C. An evolutionary Wolff’s law for trabecular architecture / S.C. Cowin //j. Biomech. Engng. - 1992. - Vol. 114. - P. 129-136.
18. Постановка начально-краевой задачи о перестройке трабекулярной костной ткани / А.А. Киченко, В.М. Тверье, Ю.И. Няшин, М.А. Осипенко, В.А. Лохов // Российский журнал биомеханики. - 2012. - Т. 16, № 4. - С. 36-52. EDN: PPMKWP
19. Mathematical modelling of trabecular bone tissue remodelling under load / V. Tverier, A. Kichenko, Y. Nyashin, V. Lokhov // Series on Biomechanics. - 2015. - Vol. 29, No. 4. - P. 39-44. EDN: XNHYCD
20. Биомеханическое моделирование трабекулярной костной ткани в состоянии равновесия / Т.Н. Чикова, А.А. Киченко, В.М. Тверье, Ю.И. Няшин // Российский журнал биомеханики. - 2018. - Т. 22, № 3. - С. 282-291. EDN: OKJHJU
21. Моделирование перестройки трабекулярной костной ткани в ветви нижней челюсти человека / Т.Н. Чикова, А.А. Киченко, В.М. Тверье, Ю.И. Няшин // Российский журнал биомеханики. - 2018. Т. 22, № 3. - С. 292-300. EDN: USBIDU
22. Киченко, А.А. Перестройка структуры губчатой костной ткани: математическое моделирование / А.А. Киченко // Российский журнал биомеханики. - 2019. - Т. 23, № 3. - С. 336-358. EDN: JNRCKN
23. Киченко, А.А. Математическое моделирование адаптации губчатой костной ткани применительно к зубочелюстной системе человека / А.А. Киченко // Прикладная математика и вопросы управления. - 2020. - № 3. - С. 35-48. EDN: RVIGZD
24. Киченко, А.А. Описание адаптации трабекулярной костной ткани посредством тензора структуры на примере зубочелюстной системы человека / А.А. Киченко // XIII Всероссийский съезд по теоретической и прикладной механике: сборник тезисов докладов: в 4 т., Санкт-Петербург, 21-25 августа 2023 г. - Т. 4. Материалы симпозиумов и исторической сессии. - С. 58-60. EDN: DPTGMC
Выпуск
Другие статьи выпуска
Рассматриваются нестационарные эффекты при работе двухканальной RQ система массового обслуживания (СМО) с обменом заявками в процессе обслуживания. Исследование направлено на выявление особенностей функционирования гибридного кол-центра с параллельным обслуживанием человеком-диспетчером и интеллектуальным голосовым ботом. Заявки, не получившие обслуживания, переходят на «орбиты», откуда осуществляются повторные звонки. По истечении случайного времени, определяемого терпеливостью клиентов, заявки покидают «орбиты». Интеллектуальный бот обладает более высокой скоростью обслуживания в сравнении с диспетчером-человеком. Поток заявок на входе распределяется между каналами обслуживания случайным образом. Учитывается, что часть клиентов в ходе обслуживания может предпочесть переход к альтернативному диспетчеру. В процессе работы у диспетчера-человека имеются короткие перерывы для отдыха и длинный перерыв на обед. В эти временные промежутки обслуживание полностью осуществляется ботом. Прерывистость режимов обслуживания приводит к нестационарным эффектам в виде текущего изменения пропускной способности СМО. Поступающие заявки описываются пуассоновским процессом, когда обслуживание открыто. Численные расчеты основаны на дискретно-событийном имитационном моделировании системы. На каждом шаге таймер модельного времени увеличивается на фиксированную величину шага. Состояние системы за временной шаг изменяется случайным образом, как марковский процесс, т. е. вероятность перехода зависит только от текущего состояния СМО и не учитывает эффектов памяти. Предложен алгоритм имитационного моделирования, основанный на разбиении всего периода работы СМО на малые интервалы, в течение каждого из которых вероятности изменений малы. Наличие альтернатив в работе СМО приводит к ветвлению процесса. Статистический характер функционирования СМО учитывается в ансамбле реализаций компьютерной модели и характеризуется дисперсией результатов. В качестве примера проведено моделирование кол-центра жилищной управляющей компании, показавшее перспективность повышения доли интеллектуальных ботов в обслуживании звонков.
Применение математических методов в различных прикладных областях играет большую роль при принятии управленческих решений. Оптимизационные модели являются неотъемлемой частью математического аппарата, используемого как различными государственными институтами, так и бизнесом для помощи лицам, принимающим решение в сложных условиях с целью проведения полного и объективного анализа предметной деятельности. Рассматривается оптимизационный подход к решению проблемы определения перечня площадок распространения информации в средствах массовой коммуникации. Сформулированы новые постановки задач целочисленного линейного программирования и многокритериальной оптимизации для моделирования распространения информации. Имплементирован алгоритм обработки статистических данных для формирования матрицы объектов-признаков. Реализованы и апробированы методы решения сформулированных задач оптимизации в задаче определения перечня площадок распространения информации. Проведен анализ чувствительности в задаче многокритериальной оптимизации, рассмотрены результаты численного моделирования при различных входных параметрах, сделаны соответствующие выводы и замечания. Актуальность продиктована нарастающей ролью информационных площадок и потребностью оптимизировать процесс принятия решений в области управления информацией.
Рассматривается актуальная проблема поиска закономерностей в больших объемах статистических данных. Инструментом анализа данных выступает регрессионный анализ. При построении регрессионных моделей исследователи зачастую стремятся только к их высокому качеству аппроксимации. Но, как отмечено в современных научных работах, одной такой метрики недостаточно. Поэтому сегодня активно развивается интерпретируемое машинное обучение. Ранее автором было предложено определение вполне интерпретируемой линейной регрессии, а задача ее построения была формализована в виде задачи частично-булевого линейного программирования. Исследования выявили высокую эффективность разработанного математического аппарата при решении задач обработки больших данных. Поэтому было принято решение расширить предложенную технологию для построения квазилинейных регрессий. В статье дано определение вполне интерпретируемой квазилинейной регрессии, включающее 6 условий. Разработан алгоритм интерпретации влияния в оцененной квазилинейной регрессии монотонно преобразованных объясняющих переменных на зависимую переменную. Задача построения вполне интерпретируемой квазилинейной регрессии формализована в виде задачи частично-булевого линейного программирования. Показано, как в этой задаче выбирать допустимые границы параметра M. Для демонстрации работоспособности предложенного математического аппарата решена задача моделирования прочности бетона на сжатие по данным, содержащим более 1000 наблюдений. Для этого использовалась программа «ВИнтер-2». В построенную модель вошли следующие преобразованные переменные: цементно-водное отношение, шлак доменной печи, пластификатор и возраст бетона. Построенная регрессия оказалась лучше по качеству аппроксимации и проще по структуре существующей модели. Дана интерпретация построенной квазилинейной регрессии. Влияние объясняющих переменных на прочность бетона в ней согласуется как с содержательным смыслом задачи, так и с другими существующими математическими моделями. Предложенная в статье технология построения вполне интерпретируемых квазилинейных регрессий обладает высоким потенциалом для решения задач обработки больших данных в различных предметных областях.
Рассматриваются проблемы зависимости сверточных нейронных сетей от цветовых параметров изображений. Выдвигаются гипотезы о том, что обучение нейронных сетей в задачах компьютерного зрения является зависимым от цветов объектов на изображениях обучающей выборки. Разработаны и проведены специальные эксперименты по обучению нейронных сетей на синтетических изображениях для подтверждения данных гипотез. Проведен анализ результатов экспериментов, который показывает наличие зависимости обучения сверточной сети от цветовых признаков объектов. Сформулированы методы преодоления некоторых выявленных проблем и задачи для дальнейших исследований.
Развивается схема построения систем прямого адаптивного управления с неявной эталонной моделью для линейных непрерывных объектов. Ключевым условием сходимости настраиваемых параметров к их идеальным значениям является постоянство «богатого» по гармоникам возбуждения системы. В реальных системах это требование трудно выполнимо, поскольку задающие воздействия должны быть отработаны системой на ограниченном интервале времени. Стремление получить в системах адаптивного управления быструю параметрическую сходимость является двигателем исследований в области гибридных адаптивных систем. В данной работе гибридность понимается в том смысле, что объект управления является непрерывным, а регулятор или какая-либо его часть - дискретным. Предлагается схема синтеза гибридного контура адаптации, в которой непрерывная часть алгоритма использует в качестве параметра дискретную оценку матрицы идеальной настройки регулятора. Ставится задача получить модель параметрической ошибки в виде уравнения линейной регрессии относительно параметров идеальной настройки регулятора. Вводится квадратичный критерий параметрического рассогласования. Соответствующая задача минимизации по искомым оценкам приводит к алгоритму наименьших квадратов, который может быть использован как в одношаговой, так и в рекуррентной форме. Проводится анализ устойчивости методом функций Ляпунова. Определены условия, при которых замкнутая система является диссипативной. Цель данной работы - показать преимущество гибридного контура адаптации перед традиционным. Для этого теоретические результаты дополняются примером численного моделирования. Сравнивается динамика системы с традиционным непрерывным алгоритмом адаптации и с гибридным алгоритмом. В обоих случаях исследуемая функция Ляпунова является невозрастающей. Это отражает устойчивость системы. Однако в традиционной системе скорость сходимости существенно меньше, что приводит к большей координатной ошибке. В системе же с гибридным контуром адаптации имеет место относительно быстрая как координатная, так и параметрическая сходимость в допустимо малую зону диссипации.
Объектом исследования является камера инертная с внутрикамерным технологическим оборудованием участка пирохимического передела - этапа технологии замкнутого ядерного топливного цикла, предназначенная для отработки технологического процесса переработки отработавшего ядерного топлива пирохимического передела модуля переработки опытно-демонстрационного энергетического комплекса. Приводится описание данных, полученных в ходе стендовых испытаний экспериментального образца инертной камеры, конструкции экспериментального образца и численных математических моделей его теплового состояния. Выполнен анализ исходной геометрии камеры и обоснованы принятые упрощения. Обоснован принцип моделирования конструктивных элементов ограждающих конструкций, основанный на применении оболочек нулевой толщины. При этом толщина стенок и теплотехнические свойства задаются математически с учетом свойств как поперек, так и вдоль оболочки. Описана методика разработки численной математической модели инертной камеры на базе метода конечных объемов с использованием программного пакета ANSYS Fluent. Приведены результаты расчетов теплового состояния исследуемой камеры как в стационарной, так и в нестационарной постановках, а также валидация разработанной численной модели на базе результатов натурных экспериментальных исследований. Показано, что разработанная численная модель обладает высокой адекватностью (отклонение от экспериментальных данных не более 5 %).
Развивается теория робастного оценивания параметров статистических моделей с привлечением аппарата теории информации. Рассмотрен подход А. М. Шурыгина, основанный на модели серии выборок со случайным точечным засорением (модели байесовского точечного засорения). Из предложенных А. М. Шурыгиным оценок, пожалуй, наиболее интересными свойствами обладают стойкие оценки. Хотя данный подход можно связать с подходом Ф. Хампеля к робастному оцениванию, необходимость при нахождении стойких оценок постулировать параметрический вид распределения засоряющей точки не позволяет считать это робастной процедурой. В первой части нашей работы предложен непараметрический способ выбора указанного распределения - посредством максимизации энтропии Шеннона в окрестности модельного распределения, ограниченной величиной дивергенции Кульбака - Лейблера. Такой способ нахождения плотности распределения засоряющей точки позволяет рассматривать получаемые оценки как робастные, причем обладающие свойством оптимальности. Полученные оценки мы называем обобщенными радикальными, поскольку их частным случаем являются радикальные оценки А. М. Шурыгина. Обобщенные радикальные оценки широко известны в зарубежных публикациях как оценки минимума логарифмической дивергенции степени плотности (гамма-дивергенции), при этом вопрос их оптимальности там не исследуется. К обобщенным радикальным относятся некоторые популярные оценки параметра сдвига: оценки Мешалкина (Уэлша), Эндрюса, Смита, Бернулли, бивес-оценка Тьюки, оценка Хьюбера типа урезанного среднего, обобщенные оценки Шарбонье. Также в первой части работы предложено использовать функционал перекрестной энтропии. Перекрестная энтропия, применяемая в качестве оптимизируемого функционала вместо энтропии Шеннона, позволяет получить семейство оценок с наиболее широким диапазоном значений параметра, задающего это семейство. К задаче максимизации перекрестной энтропии сводится задача максимизации математического ожидания функции потерь оценок максимального правдоподобия в модели байесовского точечного засорения. По этой причине обобщенные радикальные оценки могут интерпретироваться как защищенные от намеренного искажения оценок максимального правдоподобия. Во второй части работы получено другое оптимальное решение на основе формализма А. Реньи (или эквивалентного с точки зрения нашей задачи формализма К. Цаллиса), дающее новое семейство оценок, частными случаями которого также являются некоторые известные оценки. Для выбора одной оценки из семейства, определяемого разными ограничениями на дивергенцию, предложен оптимизационный подход, аналогичный таковому, приводящему к стойким оценкам, но, в отличие от последнего, остающийся непараметрическим. Основные теоретические результаты, полученные в работе, иллюстрируются во второй ее части на примере оценивания параметра сдвига косинусного распределения.
Рассматриваются задачи генерации уровней для однопользовательских 3D игр и предлагается метод, направленный на ускорение и удешевление процесса разработки игровых уровней, раскрываются подходы и технологии, примененные для решения задач генерации уровней, реализации алгоритмов генерации комнат и путей между ними. После проведения анализа существующих методов генерации уровней выбран метод генерации BSP-tree, который может создавать уникальные уровни на основе входных переменных, позволяющий сократить сроки разработки игровых уровней. Создание бесконечного уровня - сложная задача, однако с использованием некоторых полезных советов и техник становится гораздо проще. Первым шагом для организации бесконечного уровня является создание пустого объекта, который служит основой для уровня, затем можно добавлять различные элементы окружения. Для достижения эффекта бесконечного уровня предлагается использовать технику «прокрутки». Это означает, что когда игрок движется в одном направлении, то объекты в уровне перемещаются в противоположном направлении. Это создает иллюзию бесконечности и позволяет игроку продолжать исследовать новые области уровня.
Исследована устойчивость по начальной функции линейного автономного функционально-дифференциального уравнения нейтрального типа. Анализируется устойчивость по Ляпунову, асимптотическая и сильная асимптотическая, а также экспоненциальная устойчивость уравнения и их взаимосвязь. Определения всех типов устойчивости формулируются в терминах функции Коши - функции, позволяющей в явном виде записать общее решение уравнения. Основное внимание уделено исследованию устойчивости по начальной функции из пространств суммируемых функций. Используется известное представление решения функционально-дифференциального уравнения с помощью интегрального оператора, ядром которого является функция Коши. Вопросы устойчивости исследуются для уравнения с кратным запаздыванием при производной и распределенным запаздыванием при неизвестной функции. Показано, что для такого уравнения сохраняются все свойства, ранее доказанные для уравнения с кратным запаздыванием при неизвестной функции. А именно показано, что сильная асимптотическая устойчивость рассматриваемого уравнения с начальной функцией из пространства L 1 эквивалента экспоненциальной оценки функции Коши, кроме того, из любого из этих свойств следует экспоненциальная устойчивость по начальной функции в любом из пространств L p при 1≤p≤∞. При этом, как и для уравнения с кратными запаздываниями, сильная асимптотическая устойчивость в пространстве L p для некоторого p>1 может не быть равносильной экспоненциальной устойчивости.
Издательство
- Издательство
- ПНИПУ
- Регион
- Россия, Пермь
- Почтовый адрес
- 614990, Пермский край, г. Пермь, Комсомольский проспект, д. 29
- Юр. адрес
- 614990, Пермский край, г. Пермь, Комсомольский проспект, д. 29
- ФИО
- ТАШКИНОВ АНАТОЛИЙ АЛЕКСАНДРОВИЧ (ИСПОЛНЯЮЩИЙ ОБЯЗАННОСТИ РЕКТОРА)
- E-mail адрес
- rector@pstu.ru
- Контактный телефон
- +7 (342) 2198067
- Сайт
- https://pstu.ru