No 4 (2025)
- Year: 2025
- Articles: 12
- URL: https://ered.pstu.ru/index.php/geo/issue/view/480
ARTICLES
A STATISTICAL CRYSTAL PLASTICITY MODEL BASED ON CONTINUUM DISLOCATION DYNAMICS AND ITS APPLICATION TO THE COMPREHENSIVE ANALYSIS OF MICROSTRUCTURE EVOLUTION IN COPPER DURING EQUAL-CHANNEL ANGULAR PRESSING
Abstract
В ходе технологических процессов (ТП) обработки давлением металлов и сплавов, особенно с применением методов интенсивной пластической деформации, происходит изменение их зеренной структуры. При низких гомологических температурах доминирующим процессом перестроения внутренней структуры материала является фрагментация, которая реализуется за счет образования блоков ячеек и постепенной их эволюции в субзерна и новые мелкие зерна за счет разворотов друг относительно друга в ходе деформирования. Развитие вычислительных средств в последние десятилетия позволяет все более эффективно применять для анализа и совершенствования реальных ТП аппарат математического моделирования. При этом для исследования происходящих в ходе ТП изменений внутри материала целесообразно использовать многоуровневые конститутивные модели (КМ), основанные на принципах физической теории пластичности, за счет явного учета в них эволюции внутренней структуры материала и реализующихся механизмов деформирования. В работе предлагается новая многоуровневая статистическая КМ на базе континуальной дислокационной динамики для описания фрагментации зерен металлов с ГЦК-решеткой. Разработанная КМ комплексно описывает процесс фрагментации: учитывает эволюцию плотности статистически накопленных и геометрически необходимых дислокаций, а также её влияние на изменение средних значений углов разориентации, размеров фрагментов и зерен в поликристалле. Верификация результатов моделирования осуществлялась на опытных данных для четырех проходов равноканального углового прессования меди при комнатной температуре по маршруту Bc (с поворотом заготовки). Полученные в КМ результаты предсказания изменения напряженно-деформированного состояния и ключевых параметров структуры (размеров зерен, текстуры, углов разориентации и размеров фрагментов, плотностей дислокаций) находятся в удовлетворительном соответствии с данными экспериментов.
Perm Journal of Petroleum and Mining Engineering. 2025;(4):5-23
![]() | ![]() |
Identification and verification of visco-hyperelastic constants of defining relations for shock absorbing rubber
Abstract
В настоящее время в конструкциях виброзащитных устройств широкое применение нашли высокоэластические материалы, проявляющие ярко выраженные реологические свойства, способствующие снижению пиковых амплитуд колебаний в резонансных режимах. Проектирование подобных конструкций в современных условиях требует построения математической модели, наиболее ответственной частью которой являются определяющие соотношения для описания вязкоупругого и гиперупругого поведения материала и адекватная методика идентификации материальных констант. Цель данного исследования заключается в построении методики определения вязко-гиперупругих свойств материала на основании данных квазистатических испытаний. При описании гиперупругого поведения материала использовались потенциалы энергии деформации Муни-Ривлина с двумя константами и нео-Гука с учетом сжимаемости, а при учете вязкоупругости функция релаксации раскладывалась в ряды Прони. Получены полуаналитические выражения зависимости напряжений от главных растяжений и времени с учетом вязко-гиперупругого поведения материала в случаях одноосного напряженного и одноосного деформированного состояний. Для выбора оптимального распределения параметров аппроксимации относительных модулей сдвига и диапазона времен релаксации было оценено влияние изменения вязкоупругих параметров на величину гистерезиса и максимальных напряжений в цикле при различной скорости деформирования. Для вычисления параметров аппроксимации материальных свойств рассматривались результаты натурных циклических испытаний на свободное растяжение и сжатие, а также стесненное сжатие. Разработана процедура определения вязко-гиперупругих свойств материала, которая заключается в минимизации разницы между результатами натурных испытаний и вычисленного с помощью построенных полуаналитических выражений значений отклика с использованием метода Нелдера-Мида. В результате были определены гиперупругие и вязкоупругие свойства материала, которые с достаточной точностью согласуются с экспериментальными кривыми. Показано, что потенциал нео-Гука при своей более простой форме позволяет описать данные натурных испытаний наравне с потенциалом Муни-Ривлина с двумя константами. Вычисленные свойства материала были верифицированы на основании сравнения данных испытаний амортизатора и результатов численного моделирования деформирования демпфера.
Perm Journal of Petroleum and Mining Engineering. 2025;(4):24-37
![]() | ![]() |
Optimisation of the first natural frequency for heterogeneous circular plates
Abstract
Проведено исследование спектральной задачи максимизации первой собственной частоты неоднородной упругой и вязкоупругой круглой пластинки переменной жесткости при различных условиях закрепления края. В качестве управляющих функций выступают механические упругие или вязкоупругие свойства материала. Вязкоупругие свойства учитываются в рамках принципа соответствия. Рассмотрены несколько подходов к решению спектральной оптимизационной задачи: аналитические, полуаналитические и численные схемы. В модельной задаче для круглой пластинки, шарнирно опертой по краю, получено аналитическое решение – оптимальное распределение вязкоупругих модулей – длительного и мгновенного модулей, с учетом условия на среднее распределение жесткости пластинки по объему. Построено соотношение Релея, а также обобщенный функционал на базе метода множителей Лагранжа, учитывающий дополнительное условие. Найдено условие оптимальности, которое состоит в постоянстве энергетической характеристики. В случае круглой пластинки, жестко закрепленной по краю, для нахождения оптимального распределения комплексных модулей предложена численная схема оптимального управления характеристиками материала в классе полиномов, с учетом их свойств положительности и четности. В основе предлагаемой схемы лежит метод Галеркина и метод наименьших квадратов для определения коэффициентов полиномиального разложения управляющих функций. Предложена также численная схема решения спектральной задачи поиска коэффициентов полиномиального разложения управляющих функций, которая запускает эвристический поиск соответствующих параметров. Проведены расчеты по каждой из предлагаемых схем решения задач оптимизации. Выявлены недостатки и преимущества каждого подхода. Проведен сравнительный анализ результатов на примере модельных задач.
Perm Journal of Petroleum and Mining Engineering. 2025;(4):38-49
![]() | ![]() |
STRUCTURAL-PHENOMENOLOGICAL MODEL OF FATIGUE DAMAGE ACCUMULATION AND FAILURE OF INHOMOGENEOUS SOLIDS
Abstract
На основе предложенного Ю.В. Соколкиным подхода к математическому моделированию процессов усталостного разрушения структурно неоднородных сред, в рамках концепции поврежденности Ю.Н. Работнова и Л.М. Качанова в данной работе разработана двухуровневая структурно-феноменологическая модель усталостного накопления повреждений и разрушения неоднородного тела. Подробно рассмотрен переход от наиболее простых критериев и моделей усталости и моделей накопления повреждений к предложенной модели, учитывающей изменение деформационных и прочностных характеристик структурных элементов по мере роста числа циклов воздействия. Приведена постановка краевой задачи деформирования и разрушения структурно неоднородного тела с учетом усталостного накопления повреждений. Разработан алгоритм решения поставленной краевой задачи. Отмечена перспективность применения разработанной модели для прогнозирования долговечности конструкций, их механического поведения на макроуровне после предварительных циклических воздействий, а также рассмотрения кинетики процессов накопления повреждений и разрушения.
Perm Journal of Petroleum and Mining Engineering. 2025;(4):50-60
![]() | ![]() |
STRESS-STRAIN STATE OF CYLINDRICAL PARTS WHEN STRENGTHENED BY FLAT WEDGE PLATES
Abstract
Повышение качества поверхностного слоя деталей является актуальной задачей современного машиностроения, от эффективного решения которой зависит надежность работы изделиий. Формирование благоприятного напряжённо-деформированного состояния в поверхностных слоях позволяет повысить износостойкость, усталостную прочность и эксплуатационные характеристики деталей машин. На сегодняшний день существуют различные методы упрочнения деталей поверхностным пластическим деформированием, среди которых можно выделить накатывание, обкатку роликами, ультразвуковую обработку и тд. Выбор конкретного способа зависит от назначения детали, её геометрии и требований к эксплуатационным характеристикам. Особые технологические проблемы возникают при обработке цилиндрических деталей малых диаметров типа осей, силовых шпилек, валиков, которые легко изгибаются. Для нежёстких деталей, характеризующихся отношением длины к диаметру более 15, разработан способ упрочнения с использованием плоских клиновых плит. Данный метод представлен и исследован в настоящей работе. Цель работы заключается в оценке напряженно-деформированного состояния в очаге деформации и остаточных напряжений в упрочненных цилиндрических деталях в зависимости от геометрической формы рабочего инструмента и основных параметров процесса упрочнения плоскими клиновыми плитами. Методом исследования служит моделирование в программной среде ANSYS Workbench с учетом билинейной диаграммы деформирования материала образца. Проведены численные эксперименты при различных углах клинового участка β (от 30° до 90°) и степенях относительного обжатия Q (от 0,5% до 4%), что позволило определить закономерности распределения временных и остаточных напряжений, а также интенсивности пластической деформации как по длине, так и по радиусу заготовки. Установлено, что угол β клинового участка плоских плит равный 45° является наиболее рациональным, так как обеспечивает равномерное распределение остаточных напряжений, предотвращает локальные перегрузки и способствует упрочнению без повреждений поверхностного слоя. Рекомендуемая величина относительного обжатия при упрочнении плоскими клиновыми плитами должна находится в интервале Q = (1,3-1,5)%.
Perm Journal of Petroleum and Mining Engineering. 2025;(4):61–71
![]() | ![]() |
Prediction of the Polymer Reinforced Pipeline Stress-Strain State taking into Account Technological Residual Stresses
Abstract
В статье приводятся постановки и результаты численной реализации трех этапов прогнозирования напряженно-деформированного состояния полимерных армированных трубопроводов (ПАТ) под действием внутреннего давления с учетом остаточных напряжений, возникающих в процессе изготовления. Для этого на первом этапе построены процедуры прогнозирования эффективных характеристик армированного слоя; на втором этапе решена задача о прогнозировании остаточных технологических напряжений; на третьем этапе решена задача об отыскании напряженно-деформированного состояния под действием внутреннего давления с учетом остаточных технологических напряжений. Решение поставленных задач осуществлено численно, методом конечных элементов. Данное исследование носит поисковый характер, направленный на выявление проблемных зон конструкции ПАТ, формирование и отработку методологии проведения моделирования поведения под действием различных нагрузок. В рамках представленной работы рассмотрены упругие характеристики материалов и стационарные формулировки задач. Рассмотренные задачи позволили провести качественную приближенную оценку напряженно-деформированного состояния ПАТ в процессе изготовления и эксплуатации, выявить закономерности деформирования и зоны возможного разрушения элементов конструкции. Для поиска и формулирования отбракованных характеристик, рекомендаций по монтажу и эксплуатации необходимо реализовать переход к более сложным моделям описания поведения материалов, учитывающих температурные и временные зависимости механических характеристик полиэтилен высокой плотности ПЭ 100 и композиционного материала на его основе.
Perm Journal of Petroleum and Mining Engineering. 2025;(4):72–85
![]() | ![]() |
ACOUSTIC PROPERTIES OF PIEZOACTIVE HETEROSTRUCTURES MADE OF THIN LITHIUM NIOBATE PLATES
Abstract
Разработана и исследована математическая модель возбуждения поверхностных акустических волн (ПАВ) на поверхности пьезоактивной гетероструктуры «пластинка W-среза ниобата лития – подложка кремния», где W – срез, нормаль к поверхности которого имеет направления наиболее распространенных промышленных срезов Z, X, Y, 41Y, 64Y, 128Y. Наряду с ориентацией среза модель учитывает направление распространения волны и толщину пластинки по отношению к длине волны, что позволяет осуществлять выбор кристаллографических и геометрических параметров для эффективного возбуждения ПАВ. Проведено исследование фазовой скорости распространения волн Рэлея и Лява на поверхности гетероструктуры для двух случаев: когда ее поверхность свободна и когда она металлизирована и заземлена. Это позволило вычислить коэффициент электромеханической связи (КЭМС), характеризующий эффективность возбуждения в гетероструктуре как релеевских, так и горизонтально поляризованных сдвиговых ПАВ. Исследование показало, что на пластинках ниобата лития всех типах срезов КЭМС волны Лява достигает более высоких значений по сравнению с КЭМС рэлеевских волнами. Из всех промышленных типов срезов пластинок ниобата лития на данной гетероструктуре наиболее эффективными для возбуждения волн Лява оказываются X-, Y- и 128Y-срез, на которых КЭМС достигает значений 0,4, 0,38 и 0,39 соответственно. Дополнительное исследование позволило обнаружить срез 160Y, на котором КЭМС достигает абсолютного максимума. Полученные результаты предоставляют возможность существенно упростить процесс проектирования приборов и устройств на ПАВ за счет выбора оптимальных геометрических параметров гетероструктуры c пластиной ниобата лития заданного среза.
Perm Journal of Petroleum and Mining Engineering. 2025;(4):86-97
![]() | ![]() |
VISCOELASTIC-PLASTIC MODEL OF POLYLACTIDE FOR ANALYZING DEFORMATION OF BONE LATTICE SCAFFOLDS AT DIFFERENT STRUCTURAL LEVELS
Abstract
Решетчатые структуры из биоразлагаемых полимеров (скаффолды) используются в биомедицинских приложениях для замещения поврежденной костной ткани. Для понимания связанных деформационных и биологических процессов, происходящих в микронеоднородном материале, необходимо использование адекватной модели упруговязкопластичности материала. С этой целью предпринято экспериментальное исследование полилактида, включающего следующие программы нагружения: а) несколько циклов одноосного стесненного сжатия с различными скоростями деформирования, б) одноосное свободное сжатие с различными скоростями деформирования, в) одноосное свободное сжатие по программе нагружения с циклами нагружения и разгрузки с выдержками со ступенчатыми повышением максимального напряжения, а также г) одноосное свободное сжатие по треугольным отнулевым циклам нагружения при различных значениях максимального напряжения. Для численного исследования особенностей неоднородного распределения деформаций и напряжений при деформировании решетчатого скаффолда из полилактида сжатием в зависимости от скорости деформирования использована модель упругопластичности с законом упрочнения Джонсона – Кука, константы которой были определены экспериментально по данным испытаний а)-б). Обнаружено, что с увеличением деформаций сжатия образца скаффолда от 2,7 до 10,7 % начинает заметно проявляться зависимость пространственного распределения эквивалентных напряжений и скоростей деформаций от скорости деформации сжатия, приложенной к образцу, в диапазоне изменения от 1,11‧10-4 до 1,11‧10-1 с-1. Это показывает важность учета зависимости упруговязкопластических свойств полилактида от скорости деформации для корректного описания локализации пластической деформации и концентрации напряжений, влияющих на микромеханику различных физико-химических и механических процессов в скаффолдах в условиях эксплуатации. Также установлено, что использованная для описания одноосного нагружения модель упругопластичности качественно описывает поведение сплошных образцов полилактида при их нагружении: релаксацию напряжений на этапах выдержки при превышении напряжения текучести, а также накопление вязкопластических деформаций в циклах (рэтчетинг). Данные эффекты представляются важными для описания деформирования скаффолдов из полилактида в процессах отнулевого циклического нагружения с выдержками при максимальных значениях, повторяющих ходьбу человека.
Perm Journal of Petroleum and Mining Engineering. 2025;(4):98-112
![]() | ![]() |
ELASTIC-PLASTIC BENDING OF RECTANGULAR PLATES IN A 3D SETTING
Abstract
В работе построены математические модели физически нелинейных пластин в трехмерной постановке. Математическая модель строится на основании трехмерных уравнений равновесия Коши. Уравнения состояния получены в перемещениях. Используется теория малых упругопластических деформаций. При определении пластических деформаций используется критерий пластичности Мизеса. Для решения полученной нелинейной задачи используется метод переменных параметров упругости И.А. Биргера, сводящий решение нелинейных уравнений к последовательному решению линейных уравнений на каждом шаге итерационной процедуры метода. Численное исследование проводится на основе метода конечных элементов с использованием пакета Comsol Multiphysics. Исследовалась сходимость решений в зависимости от количества конечных элементов. В работе анализируется напряженно-деформированное состояние квадратной пластинки, жестко заделанной по краям, в зависимости от отношения линейного размера к толщине. Для каждого из рассмотренных случаев получены зависимости нагрузка-прогиб, распределение пластических зон по объему пластинки и распределение значений перемещения w по толщине пластинки.
Perm Journal of Petroleum and Mining Engineering. 2025;(4):113-123
![]() | ![]() |
EVALUATING THE STRESS-STRAIN BEHAVIOR OF A LARGE COMPOSITE SHELL DURING CURING IN FREE SPACE
Abstract
Изучается напряженно деформированное состояние (НДС) и формоизменение крупногабаритной цилиндрической оболочки из композитного материала, отверждаемой под действием прямого излучения от Солнца и отраженного излучения от Земли в условиях открытого космоса на околоземных орбитах. Особенностью рассматриваемой задачи является изменение механических свойств связующего материала конструкции, обусловленное двумя факторами: изменением температуры и степенью отверждения связующего. Зависимость свойств материала от изменения температуры при фиксированной степени отверждения обратима. Зависимость свойств от степени отверждения необратима. Оба фактора имеют место одновременно и могут быть описаны независимыми кинетическими уравнениями. Целью настоящей работы является разработка численной методики исследования НДС конструкций, учитывающая оба эти фактора. Предлагаемая методика является обобщением модели мгновенного стеклования И.И. Бугакова на случай малых приращений по времени и необратимого изменения материальных характеристик полимера. Методика демонстрируется на изучении НДС крупногабаритной цилиндрической композитной оболочки. Рассматриваются состояния оболочки на низкой околоземной полярной орбите, когда учитывается тепловой поток от Солнца и от Земли, и случай высокой орбиты, когда учитывается нагрев только от Солнца. В обоих случаях из-за необратимого изменения свойств связующего вследствие его отверждения появляются остаточные напряжения и деформации, обусловливающие формоизменение оболочки в виде удлинения или укорочения оси оболочки, ее изгиба и искажения формы поперечного сечения. Предложенная в работе методика исследования НДС конструкций из материалов, свойства которых в процессе нагружения необратимо изменяются, может использоваться при расчете конструкций из композитов не только в условиях космоса, но и при исследовании остаточных напряжений и формоизменения изделий в технологических задачах.
Perm Journal of Petroleum and Mining Engineering. 2025;(4):124-138
![]() | ![]() |
EXPERIMENTAL DETERMINATION OF ARTIFICIAL NOTCH OPENING DISPLACEMENTS FOR TENSION OF COMPOSITE PLATE WITH A HOLE
Abstract
Представлены результаты определения величин раскрытия искусственного надреза, моделирующего трещину в реальном композиционном материале, с помощью метода электронной спекл-интерферометрии (ESPI). Главная особенность представленного подхода заключается в том, что прецизионные оптико-интерференционные измерения величин раскрытия проводятся непосредственно на берегах надреза. Визуализированы картины интерференционных полос высокого качества, вызванные локальным удалением материала между двумя экспозициями. Полученные данные обеспечивают количественное описание полей чисто деформационных перемещений, как в «ближнем поле», так и в дальней окрестности вершины надреза. На этой основе построены распределения раскрытия надреза по его длине при растяжении прямоугольной композитной пластины с отверстием. Надрез состоит из двух ветвей, которые направлены в противоположных направлениях от точек пересечения контура отверстия и короткой оси симметрии прямоугольной пластины. Растягивающая нагрузка действует перпендикулярно линии надреза. Показано, что последовательность нанесения двустороннего надреза практически не влияет на результаты измерений компонент перемещений. Экспериментальные данные получены для надрезов различной длины при различных уровнях внешней нагрузки. Построены зависимости средних значений раскрытия надреза в его начальных точках от усредненной длины надреза и величин внешних растягивающих напряжений. Полученные зависимости, демонстрируют явную нелинейность процесса деформирования в нерегулярной зоне образцов, как при изменении длины надреза, так и внешней нагрузки. Приведенные данные представляют значительный интерес с точки зрения верификации расчетных методов, используемых при определении параметров механики разрушения для трещин в композитных конструкциях.
Perm Journal of Petroleum and Mining Engineering. 2025;(4):139-153
![]() | ![]() |
STRESS-STRAIN ANALYSIS OF A COMPOSITE CERAMIC COATING B4C-BN DEPOSITED ON A STEEL MATRIX BY LASER CLADDING
Abstract
Покрытия на основе формируют сверхтвердый поверхностный слой со сверхнизким коэффициентом трения и перспективны для использовании в парах трения с высокими относительными скоростями трения скольжения. Использование метода лазерной наплавки при нанесении покрытия обеспечивает его прочный адгезионный контакт с металлической матрицей за счет высоких температур обработки. В работе проведен анализ напряженно-деформированного состояния композиционного поверхностного слоя, формирующегося при охлаждении с высокой (до ) скоростью охлаждения. Проанализирована структура покрытия, имеющего аморфно-кристаллическое строение, и ее влияние на напряженное состояние. Установлено, что в аморфной матрице имеются включения кристаллического нитрида бора с характерными размерами 30-50 нм. Далее методами математического моделирования проведен расчет уровня остаточных напряжений для всего цикла термического нагружения. По результатам моделирования выявлен следующий механизм, что последовательная наплавка дорожек (элементов) покрытия приводит к сложному термомеханическому воздействию. В покрытии и в приповерхностном слое металлической матрицы формируются знакопеременные по времени и пространству напряжения. Это приводит к возникновению концентраторов напряжений на границе между отдельными дорожками и развитию поперечных трещин. При этом не прогнозируется значительных напряжений в направлении, перпендикулярном поверхности. Этим объясняется отсутствие трещин вдоль границы «покрытие – матрица» и наблюдаемая в триботехнических испытаниях высокая адгезия покрытия к матрице. Дополнительно исследован вопрос релаксации напряжений и упругой энергии в результате формирования трещин. Прогнозируется снижение упругой энергии на 60–70%, что снижает остаточные напряжения до уровня 400–500 МПа. Полученные результаты применимы для оптимизации режимов лазерной обработки с целью получения покрытий с наименьшим уровнем остаточных напряжений.
Perm Journal of Petroleum and Mining Engineering. 2025;(4):154–166
![]() | ![]() |


