№ 2 (2022)
- Год: 2022
- Статей: 12
- URL: https://ered.pstu.ru/index.php/mechanics/issue/view/299
- DOI: https://doi.org/10.15593/perm.mech/2022.2
Plane-Strain Extrusion of a Green Type Porous Plastic Material through a Wedge-Shaped Die
Аннотация
This paper presents the solutions for the plane-strain extrusion of porous material. We consider the problem of a stationary plastic flow through a wedge shaped die. We neglect friction between the die and the deformed material since it is rather a negative effect and should be avoided in manufacturing. The elliptic Green type yield condition and its piecewise-linear approximation are adopted for this problem. In the last case, we obtain analytical solution that links extrusion pressure and area reduction to the initial and final density of the porous material. For elliptic Green yield condition the problem reduced to nonlinear ODE that integrated numerically. The results are compared with known solution for Gurson model. The extrusion pressure predicted by the piecewise-linear model is lower than what obtained by the elliptic Green model. In turn, the pressure predicted by elliptic Green model is lower than the pressure obtained in the frame of Gurson model. At low values of area reduction, all three models predict approximately the same extrusion pressure. With a small initial porosity of the material, the Gurson model gives results that are close to the elliptic Green model, and with a large initial porosity, to the piecewise-linear Green model.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2022;(2):5-9
Контурные интегралы в нелинейной механике разрушения при смешанных формах деформирования
Аннотация
Современные знания в области механики разрушения являются ключевыми при решении проблем безопасности и прочности объектов с трещиноподобными повреждениями различной природы происхождения. Нелинейная механика разрушения при анализе напряженно-деформированного состояния в области вершины трещины базируется на одно- и двухпараметрическом подходах. Классические однопараметрические исследования предполагают изучение сингулярных величин, включающих контурный, независящий от пути интегрирования J -интеграл, коэффициент интенсивности напряжений (КИН) и т.д. Значения КИН и J -интеграла являются взаимозависимыми. Большой популярностью пользуются комбинированные методы на основе союза численных, экспериментальных и аналитических расчетов, которые позволяют получить максимально четкое описание параметров механики разрушения. Вычисление J -интеграла в конечно-элементных моделях методом реакций или напряжений является весьма эффективным, но для этого необходимы достаточно точные аналитические представления контурного J -интеграла. Существуют определенные ограничивающие условия при получении подобных формул. В многочисленных научных работах доказано, что J -интеграл в большинстве случаев не зависит от пути интегрирования, но весьма зависим от способов описания параметров напряженно-деформированного состояния, а также их производных, от размерности задачи и от степени удаления контура интегрирования от вершины трещины. В данной работе проведен обзор и приведены авторские исследования контурных интегралов в нелинейной механике разрушения для трех случаев: классическое решение Хатчинсона - Розенгрена - Райса (ХРР), контурные интегралы в градиентной теории пластичности и вычисление J -интеграла для общего случая, когда компоненты напряжений и перемещений являются функциями трех декартовых координат. Выведен обобщенный контурный J -интеграл, используемый для характеристики нелинейного амплитудного множителя.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2022;(2):10-24
Энергетические условия формирования гидрида магния
Аннотация
Предложен новый подход к анализу кинетики гидридного превращения в магнии, основанный на учете вклада механических факторов: работы внешних сил, энергии упругих деформаций и энергии, необходимой на образование единицы объема новой фазы. Энергетическое требование устойчивого состояния механической системы или инициированного изменением термодинамических условий фазового перехода состоит в том, что полная энергия системы стремится принять минимальное значение. Анализ дополнительных условий, когда зародыш гидридной фазы оказывается устойчивым (при неизменной температуре и концентрации водорода), показал, что контролирующими параметрами являются следующие величины: объемный эффект фазового превращения, отношение упругих модулей металла и гидридной фазы, работа, затраченная на образование единичного объема гидрида. В рамках данного подхода показано, что форма зародыша гидрида в виде вытянутого эллипсоида является энергетически более выгодной. Более того, наиболее устойчивым будет зародыш с существенной разницей в величине полуосей эллипсоида. Полученные в работе энергетические критерии появления новой фазы могут быть использованы в качестве дополнительного условия для оценки стабильности критического зародыша и кинетики гидридного превращения. Проведен расчет напряженно-деформированного состояния по обе стороны от границы раздела «металл - гидрид». Показано, что приграничная область фазового перехода является областью аккумулирования неоднородных интенсивных напряжений, которые могут способствовать двум параллельно протекающим процессам. Во-первых, увеличению концентрации водорода в искаженной области. Во-вторых, накоплению вблизи области некогерентной границы раздела двух фаз напряжений, превышающих более чем в два раза предел текучести при сдвиге. Наличие этих напряжений способно как инициировать появление «свежих» дефектов (дислокаций, микротрещин), так и привести к уменьшению фракционного состава порошкового материала.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2022;(2):25-38
Анализ предельных состояний цилиндрических упругопластических оболочек при растяжении и комбинированном нагружении внутренним давлением и растяжением
Аннотация
Теоретически и экспериментально исследуется упругопластическое деформирование, предельные состояния и закритическое поведение цилиндрических оболочек при растяжении и комбинированном нагружении внутренним давлением и растяжением до разрушения. Данная задача характеризуется возникновением больших деформаций, формоизменений и, как следствие, неоднородного напряженно-деформированного состояния. При численном решении подобных задач возникает проблема построения истинных диаграмм деформирования материалов. В связи с этим для исследования деформационных и прочностных свойств материалов актуально использование экспериментально-расчетного подхода, позволяющего без принятия упрощающих гипотез учесть неодноосность и неоднородность напряженно-деформированного состояния. В работе представлен новый эффективный алгоритм построения истинной диаграммы деформирования, в основе которого лежит процедура нелинейной экстраполяции диаграммы. Такой алгоритм в процессе прямого численного решения задачи последовательно выполняет построение диаграммы деформирования без использования повторных прямых расчетов, что существенно (в разы) повышает его эффективность. На основе экспериментально-расчетного подхода определены истинные диаграммы деформирования для сплошных стержней и оболочек из сталей 10ХСНД и 10Г2ФБЮ при растяжении и комбинированном нагружении внутренним давлением и растяжением до разрушения. Разрушение оболочек при растяжении происходит при меньших (в разы) значениях истинных деформаций, чем сплошных стержней. Значительные отличия в истинных напряжениях и деформациях в момент разрушения обусловлены различной локализацией деформирования сплошных стержней и оболочек после потери устойчивости пластического деформирования при растяжении. Численно и экспериментально показано, что после потери устойчивости пластического деформирования по Консидеру при растяжении цилиндрическая оболочка содержит две формы потери устойчивости до момента разрушения. Первая форма потери устойчивости, как и в сплошных стержнях, характеризуется локализацией деформаций по диаметру образца в виде «шейки», а вторая - локализацией деформаций по толщине образца, которая и определяет финальную стадию разрушения. При действии на оболочку внутреннего давления происходит вырождение первой формы потери устойчивости пластического деформирования с образованием шейки вовнутрь оболочки и наблюдается только форма потери устойчивости, вызванная локализацией деформаций по толщине оболочки.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2022;(2):39-48
Решение плоской задачи теории упругости об изгибе шарнирно закреплённой многослойной панели с круговой осью
Аннотация
Предложен метод решения плоской задачи теории упругости об изгибе шарнирно закрепленной многослойной панели с круговой осью, основанный на полиноминальной аппроксимации перемещений по толщине панели. В отличие от известных решений этой задачи, в данном случае коэффициенты аппроксимирующих полиномов находятся из условий равновесия и равенства перемещений и поперечных напряжений при переходе через границу раздела слоёв и выполнения дифференциальных уравнений равновесия в нескольких точках по толщине слоёв. Окончательно задача сводится к решению системы линейных уравнений относительно коэффициентов аппроксимирующих полиномов. Достоверность данного метода продемонстрирована путем сравнения результатов расчётов, полученных на его основе, и результатов, полученных с помощью эталонной конечно-элементной модели. Задача решена в два этапа. На первом этапе для однослойной панели выполнено исследование зависимости степени полинома от отношения среднего радиуса панели к её толщине и отношения модуля поперечного сдвига к модулю продольной упругости, характеризующих нелинейность изменения перемещений. На втором этапе на примере трёхслойной панели рассмотрено применение предлагаемого метода для расчёта многослойных панелей. При этом результаты, полученные на первом этапе, использованы при выборе начальной степени полиномов, аппроксимирующих перемещения по толщине слоёв. Предложенный в работе метод позволяет получить аналитическое решение без введения упрощающих гипотез о характере перемещения слоёв и их упругих характеристик в широком диапазоне изменения геометрических размеров и упругих характеристик слоёв панелей. Данный метод может найти применение как для верификации численных моделей, так и при выполнении расчётов на прочность многослойных панелей.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2022;(2):49-57
Моделирование процессов упругопластического деформирования монокристаллических сплавов
Аннотация
Разработана и верифицирована микромеханически мотивированная модель упругопластического деформирования двухфазных монокристаллических сплавов на никелевой основе, предсказывающая их поведение при высокотемпературных механических воздействиях с учетом наличия γ- и γ'-фаз. Использование модели актуально при выполнении уточненных расчетов напряженно-деформированного состояния охлаждаемых монокристаллических лопаток газотурбинных установок. Формулировка определяющих уравнений для каждой из фаз учитывала анизотропию упругих и пластических свойств, наличие октаэдрических систем скольжения, особенности кубической сингонии, различные механизмы упрочнения, включающие кинематическое, изотропное и латентное. Идентификация упругих и пластических констант материала для γ- и γ'-фаз производилась на основе известных диаграмм деформирования для каждой фазы. Определение эффективных свойств и диаграмм деформирования монокристаллического сплава с учетом наличия γ- и γ'-фаз производилось как на основе конечно-элементной гомогенизации для элементарного представительного объема, так и с использованием простейших реологических (структурных) моделей материала, рассматривающих последовательное и параллельное соединение фаз. Путем многовариантных вычислительных экспериментов и аналитических оценок определены зависимости упругопластических свойств двухфазных монокристаллических сплавов на никелевой основе от объемной доли γ'-фазы. В целях определения оптимальной стратегии решения рассматриваемого класса задач проведены многовариантные вычислительные эксперименты для различных видов граничных условий задачи гомогенизации, числа ячеек периодичности, форм включения γ'-фазы, объемных долей γ'-фазы, видов упрочнения, вариантов реологических моделей и даны соответствующие рекомендации. Результаты моделирования с использованием предложенной двухуровневой микроструктурной модели материала демонстрируют хорошее соответствие с экспериментальными данными для монокристаллического жаропрочного сплава CMSX-4.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2022;(2):58-72
Моделирование влияния вибраций на поверхностное натяжение капли жидкости с применением бессеточных методов
Аннотация
Применение вибрационных воздействий для целенаправленного влияния на такие процессы, как формирование капли, формирование ванны расплава и кристаллизация наплавляемого валика, позволяет управлять тепломассопереносом в жидкости, процессом кристаллизации и формой валика в технологических процессах наплавки. Рассматривается влияние вибрационных воздействий на характер движения жидкости в капле, которое отражается в изменении значения коэффициента поверхностного натяжения. Предложена математическая модель течения жидкости с учетом силы поверхностного натяжения в формализме метода гидродинамики сглаженных частиц. Этот метод позволяет напрямую учитывать влияние вибрационных воздействий через введение дополнительных граничных условий. Верификация разработанной математической модели проведена в сравнении с натурными экспериментами, в которых определялась зависимость величины коэффициента поверхностного натяжения от амплитуды скорости вибрационных воздействий. Для определения коэффициента поверхностного натяжения реализовано два метода: метод висячей капли и сталагмометрический метод. Реализованная модель удовлетворительно описывает эффект снижения коэффициента поверхностного натяжения для воды. Проведена серия численных экспериментов по определению влияния вибрационных воздействий на значение коэффициента поверхностного натяжения для стали марки 12Х18Н10Т. Выявлено, что при вибрациях с амплитудой скорости, равной 2,0 м/с, наблюдается снижение значения коэффициента поверхностного натяжения на 30 %. Снижение коэффициента поверхностного натяжения должно облегчить реализацию непрерывного стекания металла с проволоки, что может положительно сказаться на формировании металла при проволочной наплавке. Таким образом предложенная математическая модель позволяет явно моделировать эффект влияния вибрационных воздействий на значение коэффициента поверхностного натяжения и позволит в дальнейшем исследовать влияние вибрационных воздействий при аддитивном производстве.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2022;(2):73-84
Решение смешанной неосесимметричной задачи теории упругости для анизотропных тел вращения
Аннотация
Разработана методика решения смешанных неосесимметричных задач теории упругости для ограниченных тел вращения из трансверсально-изотропного материала, находящихся под действием поверхностных сил, заданных по циклическому закону. Методика предполагает развитие энергетического метода граничных состояний, основу которого составляют понятия пространств внутренних и граничных состояний, сопряженных изоморфизмом, что позволяет установить взаимно однозначное соответствие между элементами этих пространств. Во внутреннее состояние входят компоненты тензора напряжений, деформаций и вектора перемещений. В граничное состояние входят усилия и перемещения на границе тела. Доказан изоморфизм пространств состояний, позволяющий отыскание внутреннего состояния, что сводится к исследованию изоморфного ему граничного состояния. Базис формируется на основе общего решения краевой задачи эластостатики для трансверсально-изотропного тела вращения. Проводится ортогонализация пространств состояний, где в качестве скалярных произведений в пространстве внутренних состояний используется внутренняя энергия упругого деформирования; в пространстве граничных состояний используется работа внешних сил. Окончательно отыскание внутреннего состояния сводится к решению бесконечной системы линейных алгебраических уравнений относительно коэффициентов Фурье. Представлено решение задачи со смешанными граничными условиями для кругового в плане цилиндра из трансверсально-изотропного алевролита крупного темно-серого с осью анизотропии, совпадающей с геометрической осью симметрии. Заданные поверхностные силы изменяются по закону синуса и косинуса. Решение является аналитическим и характеристики напряженно-деформированного состояния имеют полиномиальный вид. Представлены явные и косвенные признаки сходимости решения задачи и графическая визуализация результатов.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2022;(2):85-97
О предварительной статистической обработке информации исследования усталостной прочности деталей машин
Аннотация
Предлагается обобщенный подход статистической обработки результатов усталостных параметров статической прочности с использованием методов математической статистики и обычной статистики при помощи среднего квадратичного отклонения, начального момента, центрального момента дисперсии и коэффициента вариации ряда распределения физико-механических характеристик материала детали, анализ исследования усталостной прочности на начальном этапе согласно статистической обработке информации, графического оформления ряда распределения напряжений, статистического анализа напряжений, применяется графический подход в виде гистограммы ряда распределения, полигона частостей, полигона накопленных частостей, сопоставляется подбор теоретических законов распределения напряжений с эмпирическим их подтверждением с более строгой оценкой соответствия законов распределения, которая производится при помощи специальных критериев согласия, например критерия Пирсона, предлагается параллельная проверка правильности выбранного подхода с использованием классических зависимостей сопротивления материалов. Произведена оценка суммарной ошибки в виде методической и непосредственного расхождения теоретических и экспериментальных значений напряжений и температуры физико-механического процесса по экспериментальным и теоретическим нормальным и касательным напряжениям, возникающим в ходе эксплуатации детали, используется температурно-временная суперпозиция в виде функции долговечности флуктуационной кинетической теории прочности, в которой температура как линейная функция может заменяться на удельную энергию, относительную деформацию, нормальные или касательные напряжения и любой другой энергетический и силовой фактор
Вестник Пермского национального исследовательского политехнического университета. Механика. 2022;(2):98-104
Контактные задачи для трансверсально-изотропного слоя
Аннотация
Рассматриваются две пространственные, одна осесимметричная и две плоские контактные задачи для трансверсально-изотропного упругого слоя, одна грань которого находится в условиях скользящей заделки. В пространственных и плоских контактных задачах плоскости изотропии могут быть параллельны или перпендикулярны граням слоя. В случае осевой симметрии плоскости изотропии параллельны граням слоя. При помощи интегрального преобразования Фурье контактные задачи сводятся к интегральным уравнениям относительно контактного давления, из которых предельными переходами можно получить известные уравнения соответствующих задач для изотропного слоя. Для решения пространственных задач с неизвестной областью контакта применяется метод нелинейных граничных интегральных уравнений, который позволяет одновременно определить контактные давления и область контакта. При выделении главной части ядра интегрального уравнения пространственной задачи, когда плоскости изотропии перпендикулярны граням слоя, используется полученное ранее в форме, свободной от квадратур, ядро интегрального уравнения соответствующей контактной задачи для трансверсально-изотропного полупространства. Интегральное уравнение осесимметричной задачи при помощи метода парных уравнений сводится к интегральному уравнению Фредгольма второго рода, для численного решения которого применяется метод механических квадратур. Получены замкнутые решения плоских контактных задач, основанные на специальных аппроксимациях символов ядер. Точность аппроксимаций возрастает с увеличением анизотропии. При этом уровень анизотропии можно охарактеризовать отличием отношения корней характеристического уравнения от единицы, соответствующей изотропному случаю. Расчеты механических характеристик и погрешностей аппроксимаций сделаны для известных трансверсально-изотропных материалов.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2022;(2):105-113
Эволюция зеренной структуры металлов и сплавов при интенсивном пластическом деформировании: многоуровневые модели
Аннотация
Хорошо известно, что эксплуатационные свойства изделий из металлов и сплавов определяются главным образом мезо- и микроструктурой последних. Структура материалов формируется и претерпевает существенные изменения в процессах изготовления деталей и конструкций методами термомеханической обработки. Важным параметром, определяющим физико-механические характеристики материалов, является зеренная структура (размер, форма, взаиморасположения зерен и включений различных фаз). Особое внимание в последние десятилетия в связи с этим уделяется процессам интенсивной пластической деформации (ИПД), позволяющим получать зеренную структуру субмикро- и нанокристаллического строения, обеспечивающую существенное повышение эксплуатационных свойств изделиям из металлов и сплавов. Разработка технологий ИПД в современных условиях немыслима без математического моделирования рассматриваемых процессов; важнейшей составляющей при разработке такого «инструментария» являются определяющие соотношения (или - более широко - конститутивные модели). В связи с вышесказанным последние должны обладать возможностью описывать эволюционирующую структуру на различных масштабных уровнях. До настоящего времени в практике деятельности разработчиков технологий переработки материалов преобладает применение макрофеноменологических моделей, основанных на классических континуальных теориях пластичности, вязкопластичности, ползучести. Со второй половины ХХ в. и до настоящего времени предлагаются различные усовершенствования конститутивных моделей указанного выше класса, в которые вводятся дополнительные параметры и кинетические уравнения для них, описывающие те или иные характеристики структуры материалов. Как правило, такие модели позволяют получить адекватную картину изменяющейся структуры, однако - для конкретных материалов и способов термомеханической обработки. При этом подобные модели, к сожалению, не обладают необходимой универсальностью, при изменении материала или способа обработки их приходится существенно «подгонять» под конкретные условия, вплоть до полного изменения входящих в модель соотношений. Краткий обзор работ, посвященных созданию и применению моделей данного класса, приведен в предыдущей статье авторов. Наиболее перспективными и обладающими значительной степенью универсальности, по мнению авторов, в настоящее время являются многоуровневые конститутивные модели, основанные на введении внутренних переменных и физических теориях пластичности (упруговязкопластичности). Обзор работ, в которых рассматриваются различные аспекты построения, модификации, численной реализации и применения таких моделей, предлагается в настоящей статье. Основное внимание уделено моделям, ориентированным на описание изменения структуры материалов, обусловленных дислокационно-дисклинационными механизмами; приведена краткая справка о моделях, в которых учитываются термически активируемые диффузионные механизмы, за счет которых реализуются процессы возврата и рекристаллизации.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2022;(2):114-146
Реализация модели поврежденности Lemaitre с кинематическим упрочнением в конечно-элементном комплексе ANSYS
Аннотация
В настоящее время одним из популярных направлений в механике разрушения является учет поврежденности и ее влияние на напряженно-деформированное состояние элементов конструкций. В данной работе проведена интеграция модели поврежденности Lemaitre с учетом комбинированного закона упрочнения, сочетающего в себе кинематический закон Armstrong - Frederick и изотропный закон упрочнения Voce, в программный комплекс расчетов по методу конечных элементов ANSYS. Модель внедрена в виде динамически подключаемой библиотеки пользовательского материала для трехмерных задач, которая апробирована на цилиндрическом образце с внешним кольцевым надрезом как в осесимметричной постановке, так и в трехмерной. Приведены модельные представления перечисленных законов, которые были объединены в общую структуру разрешающих уравнений и представлены в качестве обобщенного алгоритма реализации пользовательского материала. Данная работа демонстрирует лишь один из двух этапов верификации созданной программы - сравнение полей поврежденности при монотонном нагружении с известными в литературе данными - и не учитывает проверку поцикловой кинетики накопления пластических деформаций с экспериментальными данными на малоцикловую усталость. Результат верификации показал, что стадии и поля накопленных повреждений в численном эксперименте при осесимметричной и трехмерной постановках совпадают с аналогичными численными экспериментами известными в литературе. Помимо этого, была найдена аналогия с результатами, полученные с использованием TSL закона когезионного подхода механики разрушения. Несмотря на то что в модели когезионной прочности и Lemaitre используются два разных типа конститутивных уравнений, физический смысл данных уравнений состоит в одном - визуализация и идентификация механизмов и координат повреждений.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2022;(2):147-157