№ 1 (2017)
- Год: 2017
- Статей: 13
- URL: https://ered.pstu.ru/index.php/mechanics/issue/view/18
- DOI: https://doi.org/10.15593/perm.mech/2017.1
Расширение цилиндрической полости в сжимаемой упругопластической среде
Аннотация
Представлена новая постановка задачи о статическом расширении цилиндрической полости в упругопластической сжимаемой среде. Постановка задачи учитывает нелинейную сжимаемость, конечные деформации, закон текучести Мизеса-Шлейхера в зоне пластичности. Особенность настоящей постановки (нелинейная сжимаемость) позволяет рассматривать различные пористые материалы и изучать влияние пористости на параметры задачи. Приведен вывод аналитических соотношений, определяющих сопротивление расширению полости и размер пластической зоны для двух классов материалов - металлов и пористых сред (грунтов). На основе выявленных аналитических соотношений получены качественные выводы о характере и уровне деформаций в среде. В частности, отмечено, что сыпучая среда характеризуется аномально большой зоной пластической деформации. В рамках численного анализа рассмотрен ряд металлов, сплавов и мягких грунтов. Получены оценки области значений наиболее важных параметров. Для металлов это безразмерный объемный модуль упругости , на основе которого с помощью аппроксимирующих зависимостей вычисляется удельное сопротивление расширению и размеры полости и зоны уплотнения. Поскольку мягкие грунты являются материалами с выраженными нелинейными свойствами, то характер зависимости параметров обретает здесь более сложную форму - дополнительное влияние пористости материала на сопротивление расширению и размеры деформации. Результаты численных расчетов позволили провести сравнительный анализ влияния пористости на размер зоны уплотнения в окрестности полости для песчаных и глинистых грунтов. Приведенные оценки сопротивления среды расширению цилиндрической полости полезны для формулировки и уточнения прикладных соотношений теории динамического проникания тонких (заостренных) ударников.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2017;(1):5-23
Влияние состояния разломов литосферных плит на стартовое землетрясение
Аннотация
В более ранней работе авторами была построена модель стартового землетрясения от подготовки до свершения события. Модель основывалась на принятой в ряде работ постановке граничной задачи, представляющей две полубесконечные литосферные плиты с прямыми границами, моделируемыми пластинами Кирхгофа, медленно движущимися навстречу друг другу. Плиты находятся на трехмерном линейно деформируемом основании. Предполагается, что поверхность основания является границей Конрада, разделяющей верхнюю, гранитную, и нижнюю, базальтовую, структуры коры Земли. Для таким образом поставленной граничной задачи исследуется изменение концентрации контактных напряжений под литосферными плитами для случаев наличия и отсутствия дистанции между плитами в зоне сближения литосферных плит. При непосредственной их близости концентрация напряжений в зоне сближения, являющаяся разломом, становится сингулярной. Из практики решения граничных задач для линейно упругих материалов это означает разрушение зоны такой концентрации напряжений, то есть возникновение землетрясения. Для этой модели были изучены различные варианты последствий от воздействия на литосферные плиты. В результате оказалось, что к таким же последствиям приводят реально происходившие сбросовые и разламывающие землетрясения. Исследование осуществлялось для статической задачи при вертикальных внешних воздействиях на литосферные плиты в предположении отсутствия граничных напряжений на торцах. Исследование, основанное на применении метода блочного элемента, выполненное в настоящей работе, направлено на построение более сложной модели в предположении присутствия на торцах плит, т.е. торцевых границах разломов, напряжений - перерезывающих сил и изгибающих моментов. Построены соотношения, позволяющие учитывать присутствие различных вариантов указанных граничных условий на берегах разломов и выяснено их влияние на возможность возникновения землетрясения
Вестник Пермского национального исследовательского политехнического университета. Механика. 2017;(1):24-38
Численное исследование физически нелинейной задачи о продольном изгибе трехслойной пластины с трансверсально-мягким заполнителем
Аннотация
В работе проведено численное исследование физически нелинейной задачи о продольном изгибе бесконечно длинной трехслойной пластины с трансверсально-мягким заполнителем. Предполагается, что в правом торцевом сечении несущие слои жестко защемлены и отсутствует адгезионное соединение заполнителя с опорным элементом, в левом торцевом сечении несущие слои шарнирно оперты на абсолютно жесткие в поперечном направлении диафрагмы, склеенной с торцевым сечением заполнителя. Задача рассматривается в одномерной геометрически нелинейной постановке. Предполагается, что зависимость между касательным напряжением и деформацией поперечного сдвига соответствует идеальной упругопластической модели, т.е. модули касательных напряжений в заполнителе не превосходят некоторого предельного значения. Это условие означает недопущение разрушения конструкции и соответствует учету физической нелинейности в заполнителе по модели идеальной упругопластической модели. Обобщенная постановка сформулирована в виде задачи поиска седловой точки некоторого обобщенного функционала Лагранжа. Исследованы свойства функционала. Доказана выпуклость, полунепрерывность снизу и коэрцитивность по основным переменным (перемещениям точек срединных поверхностей несущих слоев), вогнутость, полунепрерывность сверху и антикоэрцитивность по множителям Лагранжа (касательным напряжениям в заполнителе). Это дало возможность при доказательстве теоремы существования и единственности использовать общую теорию существования седловых точек. Для решения задачи предложен двухслойный итерационный метод типа Удзавы, каждый шаг которого сводится к решению линейной задачи теории упругости и нахождению проекции на выпуклое замкнутое множество. Установлена сходимость метода. С помощью разработанного в среде MatLab комплекса программ проведены численные эксперименты для модельной задачи. Проведен анализ полученных результатов. Результаты численных экспериментов соответствуют физической картине.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2017;(1):39-51
Влияние сжимаемости газа в цилиндре насоса на продольные колебания штанговой колонны
Аннотация
Разработана модель продольных колебаний штанговой колонны глубинно-насосной установки, учитывающая сжимаемость газожидкостной смеси в цилиндре насоса. Был применен подход, использующий формулировку задачи в виде квазивариационного неравенства. Решение этой задачи может быть сведено к последовательности задач негладкой минимизации. Данный подход является достаточно универсальным и может быть использован для колонн в сильно искривленных скважинах. В данной модели газовый компонент смеси подчиняется закону Бойля. Сжимаемость жидкостного компонента, растворение газа в жидкости и выделение газа из жидкости в этой модели не учитываются. На базе предложенной модели рассмотрено влияние содержания газа в смеси, заполняющей цилиндр насоса, на динамограмму штанговой колонны в верхнем сечении. Кроме того, рассмотрено различие поведения штанговой колонны в двух случаях: когда газ равномерно распределен (но не растворен) в жидкости, заполняющей цилиндр насоса и когда газ в цилиндре насоса занимает локализованный объем. Если газ распределен в откачиваемой смеси, то сегменты динамограмм, соответствующие фазе сжатия, имеют более слабый наклон, чем в фазе расширения. Это объясняется тем, что масса газа, а следовательно, и сжимаемость смеси в фазе сжатия больше, чем в фазе расширения. Этот факт также объясняет ослабление свободных колебаний, которое гораздо более значительно для той фазы, когда плунжер движется вниз, чем для фазы, когда плунжер движется вверх. Если газ занимает локализованный объем в цилиндре насоса, то в фазе сжатия имеет ту же форму, что и в фазе растяжения. Повышенного ослабления вибраций при движении плунжера вниз в этом случае не наблюдается.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2017;(1):52-62
Индикация термоупругой неустойчивости скользящего контакта с помощью заглубленной пьезокерамической прослойки
Аннотация
При эксплуатации триботехнических агрегатов и устройств может возникнуть неустойчивость скользящего фрикционного контакта, сопровождающаяся ростом температуры и напряжений, что часто приводит к возникновению нештатных ситуаций, аварий и т.п. Для предупреждения таких явлений предлагается система мониторинга возникновения неустойчивости, основанная на опосредованной индикации температуры и напряжений с помощью измерения и анализа индукции и напряженности электрического поля на встроенной пьезокерамической прослойке. С целью моделирования такой схемы мониторинга термоупругой неустойчивости рассматривается нестационарная динамическая контактная задача термоэлектроупругости о скольжении жёсткого тела в виде полуплоскости по поверхности упругого покрытия. Контакт осуществляется с учётом сил трения и тепловыделения от трения. Пьезокерамическая прослойка располагается между упругим покрытием и недеформируемым основанием и жёстко сцеплена с ними. При этом пьезокерамическая прослойка термоизолирована от термоупругого покрытия, вектор предварительной поляризации направлен вдоль вертикальной оси, ортогональной граням прослойки, на которых располагаются электроды с подведённой разностью потенциалов. Решение поставленной задачи о скользящем термофрикционном контакте строится с помощью интегрального преобразования Лапласа, позволяющего представить основные физические параметры задачи: температуру, напряжения, смещения, индукцию и напряженность электрического поля, электростатический потенциал в виде контурных квадратур обратного преобразования Лапласа или свёртки Лапласа. После изучения свойств подынтегральных функций и изолированных особых точек в комплексной плоскости переменной интегрирования квадратуры вычисляются методами теории функций комплексного переменного. Вычисление квадратур приводит к бесконечным рядам по полюсам подынтегральных функций, являющихся собственными числами задачи. Решение в такой форме позволяет достаточно просто установить области устойчивых и неустойчивых решений задачи на бесконечном временном интервале. Формулы электрического тока и его напряженности, возникающих на пьезокерамической прослойке, совместно с формулами температуры и напряжений на скользящем контакте позволяют во времени диагностировать процесс возникновения и развития термоупругой неустойчивости скользящего контакта.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2017;(1):63-84
Расчетно-экспериментальные исследования прочности корпуса вентилятора турбореактивного двигателя при ударном воздействии оторвавшейся лопатки
Аннотация
Рассмотрена задача по удержанию оторвавшейся лопатки вентилятора турбореактивного авиационного двигателя для обеспечения безопасности пассажиров и экипажа самолета. Указаны основные направления конструирования вентилятора по обеспечению данного требования, одному из которых посвящена данная работа. Представлены расчетные и экспериментальные исследования прочности нескольких вариантов конструкций корпусов вентиляторов. Показано различие результатов расчета и эксперимента для металлов, отличающихся пластичностью. Дан обзор и описание известных методик расчета металлической бронезащиты, указан их общий недостаток, заключающийся в неучете предельной деформации материала. Предложена методика расчета металлического корпуса, учитывающая предельную деформацию материала. При этом указано граничное значение предельной деформации материала, разделяющее металлы на жесткие и пластичные, для которых приведены математические формулы расчета стенки бронезащиты. Показана хорошая сходимость расчета по предложенным соотношениям с результатами эксперимента для титановых и алюминиевых сплавов. Описана конструктивная схема, проведены расчеты прочности и степени поглощения энергии удара оторвавшейся лопатки специально разработанным устройством в виде гофрированного защитного пакета (или чехла) из высокопрочной ткани, который может быть установлен на любую конструкцию корпуса вентилятора как из металлических, так и из полимерных композиционных материалов без принципиальных изменений схемы компоновки. Приведено расчетное обоснование преимущества конструкции гофрированного защитного пакета по сравнению с конструкцией традиционного пакета. Приведена оценка эффективности применения алюминиевых сплавов в качестве материала корпуса вентилятора как альтернативы титановым сплавам. Описана экспериментальная установка для испытания полногабаритных корпусов. Показаны недостатки установки для испытаний модельных корпусов. Описана проблема применения полимерного композиционного материала для силового корпуса вследствие его резания лопатками вентилятора при отрыве одной из лопаток. Представлена проектная оценка применения полимерных композиционных материалов для корпуса вентилятора, указаны пути повышения баллистической стойкости корпуса из полимерных композиционных материалов.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2017;(1):85-103
Экспериментальные исследования влияния повышенных и высоких температур на прочностные и деформационные свойства комбинированных стеклоорганопластиков
Аннотация
Актуальной задачей является исследование и анализ влияния повышенных и пониженных (эксплуатационных) температур на механические свойства и механизмы разрушения композиционных материалов в условиях, максимально приближенных к эксплуатационным режимам, а также установления температурных зависимостей упругих и прочностных характеристик волокнистых композитов, применяемых в ответственных конструкциях. Исследовались механические характеристики композиционных материалов на основе комбинированной стеклоорганоткани и двух типов связующего - полиимидного (КМ-1) и эпоксидного (КМ-2). Целью настоящей работы было исследование деградации механических характеристик композиционных материалов при повышении температуры. Исследования проводились на универсальной электромеханической системе Instron 5882, в состав которой входит термокамера и бесконтактный видеоэкстензометр. Исследования при высокой температуре 500 ºС проводились на сервогидравлической испытательной системе Instron 8850 с использованием высокотемпературной муфельной печи. Получены результаты экспериментальных исследований влияния повышенных и высокой температур на деформационные и прочностные свойства композитов двух рецептур КМ-1 и КМ-2 при растяжении и смятии вдоль утка и вдоль основы. Проведены механические испытания на растяжение и смятие при нормальной и повышенных (100, 250 ºС) температурах для КМ-1 и КМ-2. Для композита КМ-1 также проведены испытания на растяжение при высокой (500 ºС) температуре. Построены диаграммы деформирования, при растяжении и при смятии в условиях повышенных температур выявлена стадия закритического деформирования, определены значения упругих коэффициентов, пределов прочности и прочности при смятии. Построены и проанализированы зависимости деградации механических свойств композитов при повышении температуры. Выявлено преимущество рецептуры КМ-2 над КМ-1 по комплексу механических характеристик, поддающихся сравнению, в диапазоне температур от 22 до 250 ºС при растяжении и температуре 22 ºС при смятии в направлении основы.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2017;(1):104-117
Расчетно-экспериментальная оценка прочности сегмента композитного шпангоута с применением метода in-situ рентгеновского контроля
Аннотация
В настоящей работе на примере конструктивно-подобного элемента типа «сегмент шпангоута», имеющего П-образное сечение и выполненного из углепластика по технологии пропитки под давлением (RTM), с использованием метода микрофокусной рентгенографии проведена качественная и количественная оценка микроповреждений, возникающих в конструкции при проведении механических испытаний. В ходе эксперимента были получены рентгенографические изображения образца при фиксированных уровнях растягивающего усилия. Полученные снимки анализировались на предмет появления расслоений как основного вида повреждений, возникающих в микроструктуре испытываемой детали. Были определены места локализации и ориентировочные значения размеров межслоевых трещин при различных уровнях нагрузки. Отмечено, что после снятия нагрузки межслоевые трещины «смыкаются» и с трудом идентифицируются средствами неразрушающего контроля. Для численного анализа напряженно-деформированного состояния (НДС) рассматриваемой конструкции была разработана трехмерная компьютерная модель с явным описанием её слоистой структуры. Такой подход позволил провести оценку межслоевых отрывных и сдвиговых напряжений, которые определяют разрушение конструкции, зафиксированное в эксперименте. Решение данной задачи осуществлялось методом конечных элементов (МКЭ) в пакете ANSYS Workbench с использованием высокопроизводительного вычислительного комплекса. Сравнение результатов численного моделирования механических испытаний и данных рентгенографического контроля показало хорошее качественное совпадение - области локализации максимальных межслоевых напряжений, определяющих характер разрушения детали, возникают в зонах скругления профиля. Отличие оценок несущей способности, полученных экспериментальными и численными методами, составило порядка 23 %, что может быть скорректировано после проведения натурных испытаний серии подобных образцов и уточнения прочности материала на межслойный отрыв.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2017;(1):118-133
Прерывистая текучесть: модели, основанные на физических теориях пластичности
Аннотация
Прерывистая текучесть как проявление неустойчивости пластического деформирования обнаруживается в ряде металлов и сплавов в определенных диапазонах скоростей деформации и температур при различных нагружениях. Процессы неупругого деформирования и свойства поликристаллических материалов на макроуровне, как показывают многочисленные экспериментальные и теоретические исследования, существенным образом определяются состоянием эволюционирующей мезо- и микроструктуры материала. Для описания пластической неустойчивости необходимо создание математических моделей, позволяющих учитывать самоорганизацию микроструктурных процессов, которые в конечном счете могут привести к спонтанному появлению локализации деформации. В статье (Трусов П.В., Чечулина Е.А. Прерывистая текучесть: физические механизмы, экспериментальные данные, макрофеноменологические модели // Вестник ПНИПУ. Механика. 2014. № 3. С. 186-232) был представлен краткий обзор работ, посвященных описанию физических механизмов и экспериментальных исследований прерывистой пластичности, рассмотрены также наиболее известные макрофеноменологические модели. Однако для корректного описания пластической деформации и ее неоднородности, для учета наиболее значимых реальных физических механизмов, определяющих её и сопутствующих ей, необходимо изучать поведение материала на более низких, чем макроуровень, масштабных уровнях. Пластическая деформация обусловлена неоднородным (во времени и пространстве) движением элементарных носителей пластической деформации - дислокаций. Дислокации и барьеры дислокационной (Ломера-Коттрелла, иммобильные дислокации) и недислокационной природы (облака примесных атомов, твердых частиц, выделений вторичной фазы и т.д.) позволяют описывать многоуровневые модели на мезо- и микроуровне. В настоящей статье предлагается краткий обзор теоретических работ, основанных на физических теориях пластичности, посвященных описанию особенностей деформирования сплавов в температурно-скоростных диапазонах, в которых существенное влияние на поведение материалов оказывают диффузионные процессы. Особое внимание уделено описанию эффекта Портевена-Ле Шателье (ПЛШ), возникновение которого большинство авторов связывают с взаимодействием дислокаций с атмосферами примесных атомов.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2017;(1):134-163
Динамическая осесимметричная задача прямого пьезоэффекта для круглой биморфной пластины
Аннотация
Рассматривается динамическая осесимметричная задача для круглой биморфной конструкции, состоящей из металлической подложки и пьезокерамической аксиально поляризованной пластины. Ее изгибные колебания осуществляются за счет действия на торцевой поверхности механической нагрузки (нормальных напряжений), являющейся произвольной функцией радиальной координаты и времени. Учитывается жесткое и шарнирное закрепление цилиндрической поверхности пластины. Исходные расчетные соотношения сформулированы для пьезокерамического материала с гексагональной кристаллической решеткой класса 6 mm. Для решения задачи теории электроупругости в трехмерной постановке используются конечные интегральные преобразования Ханкеля по аксиальной координате и обобщенное преобразование (КИП) по радиальной переменной. При этом на каждом этапе решения проводится процедура стандартизации, которая позволяет реализовать соответствующий алгоритм преобразования. В первом случае краевые условия представляются в смешанной форме, а во втором приводятся к однородным путем введения вспомогательных функций. Данный подход позволяет получить точные, в рамках используемых моделей, расчетные соотношения в наиболее общем виде. Построенное замкнутое решение позволяет определить частотный спектр собственных осесимметричных колебаний, напряженно-деформированное состояние и характер изменения индуцируемого электрического поля биморфной пластины. Это дает возможность на основании анализа связанности электрических и механических полей напряжений научно обосновать конструктивные решения проектируемых приборов, определить способ фиксации электрического сигнала, подобрать все геометрические, а также физические характеристики типовых элементов пьезокерамических преобразователей. Разработанный алгоритм решения позволяет также решать задачи теории упругости и электроупругости для круглых толстых и тонких пластин с произвольным количеством слоев при наиболее общих условиях загружения без использования кинематических гипотез.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2017;(1):164-180
Стохастическая модель для расчета остаточных напряжений в поверхностно упрочненном полом цилиндре в условиях ползучести
Аннотация
Предложена стохастическая модель для расчета остаточных напряжений в поверхностно упрочненном полом цилиндре в условиях ползучести при осевом растяжении. Выполнен цикл экспериментальных исследований для определения распределения деформации ползучести по длине упрочненного образца из сплава Д16Т при температуре 125 °С при четырех условиях растягивающего напряжения: 353; 385; 406.2; 420 МПа. Установлен существенный разброс деформации ползучести локальных участков, деформации которых могут отличаться друг от друга в 3-4 раза. Предложены феноменологические стохастические уравнения ползучести и выполнена проверка их адекватности экспериментальным данным. Разработана методика численного решения стохастической краевой задачи ползучести поверхностно упрочненных изделий, реализованная для полого цилиндра из сплава Д16Т в условиях осевого растяжения. Теоретически исследована релаксация остаточных напряжений в полом цилиндре и получены статистические оценки для полей остаточных напряжений в процессе ползучести в различные моменты времени. Приводятся доверительные интервалы для компонент тензора остаточных напряжений в процессе ползучести образца в различные временные сечения. Выполнены экспериментальные исследования по влиянию растягивающей нагрузки на релаксацию остаточных напряжений в полом цилиндрическом образце из сплава Д16Т. Приводятся финишные экспериментальные диаграммы для осевой и окружной компонент тензора остаточных напряжений после процесса ползучести длительностью от 83 до 166 часов. Теоретические и экспериментальные исследования позволили сделать вывод о практически полной релаксации остаточных напряжений в исследуемых цилиндрических образцах вследствие ползучести в условиях растяжения в течение 100-160 часов. Обсуждаются вопросы применимости метода колец и полосок для экспериментального определения остаточных напряжений в условиях существенного разброса данных при ползучести и интенсивной релаксации остаточных напряжений.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2017;(1):181-207
О возможностях и ограничениях усредненного описания неупругого поведения хрупких пористых материалов в стесненных условиях
Аннотация
В работе с использованием компьютерного моделирования методом подвижных клеточных автоматов проведено исследование особенностей неупругого деформирования и разрушения микроскопических областей пористых хрупких материалов в условиях неравноосного сжатия. Акцент в исследовании сделан на анализе применимости классических макроскопических критериев пластичности и прочности для интегрального описания механического отклика представительных объемов микроскопического масштаба. Проанализирована связь параметров интегрального механического отклика микроскопических областей с объемным содержанием пор, характером их пространственного распределения в объеме материала и прочностными свойствами материала стенок каркаса. На примере осевого сжатия образцов в условиях постоянного бокового давления исследована стадийность процессов накопления и роста повреждений в стенках пористого каркаса и их связь с интегральным неупругим откликом. Показано, что с ростом величины бокового давления происходит смена характера разрушения пористого материала от упруго-хрупкого к неупругому, локализованному в форме полосы сдвига, и далее к объемному катакластическому. Значения характерных боковых давлений, при которых происходит смена механизма разрушения, существенным образом зависят от чувствительности прочности стенок каркаса к локальному давлению. Анализ результатов моделирования показал, что традиционные условия (критерии) пластичности, учитывающие вклад локального среднего напряжения в линейном приближении, адекватно описывают отклик микроскопических представительных объемов хрупких пористых материалов в стесненных условиях только от начала неупругого деформирования до стадии формирования системы относительно коротких невзаимодействующих трещин. Важно отметить, что разупрочнение представительных микрообъемов хрупких пористых материалов, деформируемых в стесненных условиях, связано не с потерей целостности, но с более поздними процессами формирования полос локализованного сдвига и коллапсом пор в уже фрагментированном материале. Это дает основания предполагать, что экспериментально измеряемые значения прочности стесненных образцов как максимума интегрального сопротивления могут быть существенно завышены по сравнению с истинными значениями (соответствующими фрагментации образцов). Установлено, что условия потери целостности хрупких пористых материалов в стесненных условиях адекватно описываются с использованием «линейных» критериев разрушения, параметры которых определены не из стандартных испытаний на одноосные сжатие/растяжение, но на основе тестирования образцов в стесненных условиях.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2017;(1):208-232
Об учете нелинейных и связанных эффектов тепловой задачи и фазовых переходов при моделировании технологии контактного термосилового поверхностного упрочнения металлических сплавов
Аннотация
Представлена новая постановка задачи о статическом расширении. Рассмотрена постановка контактной термосиловой задачи с учетом фазовых превращений и численное решение её термоструктурной подзадачи, возникающей при моделировании процесса поверхностного упрочнения электромеханической обработкой (ЭМО). Предложена математическая модель формирования структуры металлических сплавов при ЭМО на основе совместного анализа расчетных данных о динамике температурных полей и термокинетической диаграммы на примере титанового сплава Ti6Al2V. Кинетика характерного для ЭМО мартенситного превращения при этом описана с использованием эмпирического уравнения Койстинена-Марбургера. Приведен алгоритм решения тепловой задачи на основе метода конечных элементов в слабой форме Галеркина. При аппроксимации расчетной области, представленной в виде полупространства, использовались специальные конечные элементы типа Зенкевича. На основе серии вычислительных экспериментов исследовано влияние величины временного шага на точность решения задачи. Проведен анализ значимости связанных и нелинейных эффектов, характерных для высокоскоростных высокотемпературных тепловых процессов, в частности изменения теплофизических свойств металла, наличия скрытой теплоты фазовых переходов, теплового излучения и зависимости теплофизических свойств металла от температуры. Учет рассматриваемых эффектов в процессе решения термоструктурной задачи осуществлялся с применением двух методик: решение полностью связанной нелинейной задачи на основе прямого итерационного метода Пикара, использованного совместно с формулой релаксации для ускорения сходимости; решение квазилинейным способом, когда значения нелинейных членов вычисляются на основе распределения температуры, полученного на предыдущем временном шаге. На основе численных экспериментов приведен анализ возникающих нестационарных температурных полей, а также картины распределения структурных областей, характерных для ЭМО титановых псевдоальфа-сплавов.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2017;(1):233-250