№ 1 (2018)
- Год: 2018
- Статей: 10
- URL: https://ered.pstu.ru/index.php/mechanics/issue/view/14
- DOI: https://doi.org/10.15593/perm.mech/2018.1
Математическая модель развития трещины гидроразрыва пласта в трехмерной пороупругой среде
Аннотация
В настоящее время гидроразрыв пласта (ГРП) является одним из самых распространенных методов увеличения нефтеотдачи, применяемым при промышленной разработке нефтегазовых месторождений. В то же время используемые модели развития трещины ГРП являются во многом упрощенными: геометрия трещины предполагается плоской и заранее заданной, практически всегда отсутствует полноценный учет геомеханических эффектов, эволюция трещины зачастую описывается эмпирическими критериями. Несмотря на успешность применения таких моделей, их возможностей недостаточно для решения ряда важных задач разработки. В работе представлена полностью трехмерная самосогласованная физико-математическая модель развития крупномасштабной трещины гидроразрыва пласта. Модель включает в себя несколько групп уравнений, включая неизотермическую пороупругую модель Био для описания поведения вмещающей трещину среды, двумерные уравнения смазочного слоя, описывающие течение в трещине, соответствующие условия согласования на границе «трещина-среда». Геометрическая модель трещины предполагает, что она является произвольной гладкой поверхностью с краем. Эволюция поверхности трещины в ходе ее роста описывается физически обоснованными критериями на основе векторной формы J-интеграла Райса-Черепанова. Модель пригодна как для анализа непосредственно процесса развития трещины ГРП, так и для анализа вызванных наличием трещины фильтрационных и геомеханических эффектов в ходе эксплуатации скважины с ГРП. Основное назначение предложенной модели - согласованное описание процесса развития трещины ГРП в достаточно общей постановке с минимумом априорных допущений о характере протекания процесса и вместе с тем пригодное для использования на практике с применением современных вычислительных подходов. В связи с этим в работе дан краткий обзор вычислительных алгоритмов, пригодных для практической реализации модели.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2018;(1):5-17
Численное исследование термомеханического поведения кристаллизующейся полимерной среды с учетом больших деформаций
Аннотация
Предложены феноменологические определяющие соотношения для кристаллизующейся полимерной среды, полученные в рамках нелинейной механики деформируемого твердого тела. Соотношения основаны на представлении среды в виде композиции расплавленного и полностью закристаллизованного материала с учетом истории непрерывного зарождения и деформирования новой фазы в интервале температур фазовых превращений. Выполнена общая постановка эволюционной краевой задачи нелинейной механики полимерных материалов в условиях протекания фазовых переходов с использованием предложенных определяющих соотношений. Рассмотрены алгоритмы численной реализации теплокинетической задачи и задачи определения напряженно-деформированного состояния затвердевающей системы для случая плоского деформированного состояния. Разработана процедура линеаризации, удобная с точки зрения построения численных алгоритмов решения поставленных эволюционных краевых задач, использующая предположение о близости каждой промежуточной конфигурации к текущей, что соответствует процедуре наложения на конечные деформации кристаллизующейся среды малых деформаций закристаллизовавшихся частиц. Процедура линеаризации исходной постановки задачи механики реализована с учетом температурных и структурных деформаций. Разработан и реализован численный алгоритм решения поставленной плоской краевой задачи с целью исследования особенностей эволюции напряженно-деформированного состояния в полимерной конструкции. Предложена и реализована методика построения дискретного аналога поставленных краевых задач, основанная на использовании процедуры метода Галеркина с выбором базисных функций с компактным носителем по методу конечных элементов. При этом в качестве узловых неизвестных приняты приращения функций перемещений на текущем временном шаге. Установлены закономерности формирования дефекта типа раковины в кристаллизующемся полимерном цилиндре.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2018;(1):18-28
Пластическое деформирование материалов, чувствительных к виду напряженного состояния
Аннотация
Рассматриваются основные положения и уравнения теории пластичности материалов, чувствительных к виду напряженного состояния, т.е. материалов, которые имеют различные кривые пластического деформирования при одноосном растяжении, сжатии, кручении (сдвиге). Таким образом, для таких материалов не существует единой кривой пластического деформирования при лучевых (простых) процессах нагружения. Рассматриваемая теория пластичности относится к теориям пластического течения при комбинированном упрочнении, в которой радиус поверхности нагружения принимается зависящим от первого инварианта тензора напряжений и параметра вида активного напряженного состояния, а определяющие функции эволюционного уравнения для смещения поверхности нагружения - от параметра вида добавочного напряженного состояния (состояния микронапряжений). Параметр вида определяется отношением третьего инварианта ко второму инварианту в степени 3/2 соответствующих девиаторов и равен при сжатии -1, при растяжении +1, а при сдвиге 0. В рамках этой теории рассматривается пластическое изменение объема (разрыхление) в случае зависимости радиуса поверхности нагружения от первого инварианта тензора напряжений. Для описания процессов накопления повреждений приводится кинетическое уравнение, базирующееся на работе микронапряжений на поле пластических деформаций. В этом уравнении энергия разрушения принимается зависящей от первого инварианта тензора напряжений и параметра вида состояния микронапряжений. Приводятся материальные функции, замыкающие теорию, и метод их определения. Анализируются результаты теоретических и экспериментальных исследований упругопластического деформирования образцов из алюминиевого сплава Д16Т по двухзвенным траекториям деформаций в виде «веера», а также образцов из стали 30XГСА при нагружении по двухзвенным ортогональным траекториям напряжений. Получено удовлетворительное соответствие расчетных и экспериментальных результатов. Рассмотрен эффект «расщепления деформаций», приводящий к тому, что лучевым траекториям деформаций (напряжений) могут отвечать нелучевые траектории напряжений (деформаций), а плоским траекториям - неплоские.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2018;(1):29-39
Экспериментальное и теоретическое исследование взаимосвязи фазовой и структурной деформаций в сплавах с памятью формы
Аннотация
Структура ориентированного мартенсита, определяющая макроскопическую деформацию в сплавах с памятью формы (СПФ), может быть сформирована двумя способами: непосредственно из аустенитной фазы в результате прямого фазового превращения под нагрузкой, а также из хаотического мартенсита при его структурном превращении в процессе изотермического нагружения. Деформация, приобретенная первым способом, называется фазовой, вторым - структурной, однако различие в названиях отражает лишь разность механизмов их формирования, тогда как конечный продукт - ориентированный мартенсит - одинаков для обоих типов деформации. В связи с этим некоторые феноменологические модели для СПФ учитывают единообразие этих двух составляющих деформации, определяя их взаимосвязь через диаграммы прямого превращения F1 (для фазовой деформации) и мартенситной неупругости F2 (для структурной деформации). Лежащее в основе этих моделей теоретическое представление содержит гипотезу о независимости процесса дальнейшего деформирования ориентированного мартенсита от способа его формирования и включает три материальные функции: F1, F2, а также связывающую их функцию f . В настоящей работе проведено экспериментальное исследование на проволочных образцах из никелида титана, целью которого являлось подтверждение гипотезы и установление материальных функций, используемых данным теоретическим представлением. Предложен новый способ определения функции f , который может использоваться в качестве проверочного эксперимента. Определены пределы выполнимости гипотезы как для дальнейшего изотермического деформирования образцов с начальными фазовой и структурной деформациями, так и для процессов, связанных с их последующим нагревом. Экспериментально показано совпадение диаграмм доориентации для таких образцов, что служит иллюстрацией эффекта перекрестного упрочнения.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2018;(1):40-57
НЕОДНОРОДНОСТИ ПОЛЕЙ ДЕФОРМАЦИЙ В ЗЕРНАХ ПОЛИКРИСТАЛЛИЧЕСКИХ МАТЕРИАЛОВ И ЗАДАЧА ЭШЕЛБИ
Аннотация
Предложен и реализован способ вычисления неоднородных полей деформаций в зернах поликристаллических материалов. Вычисления основаны на разработанном ранее методе решения краевой задачи механики неоднородных поликристаллических тел с помощью оригинального варианта теории возмущений, основанного на аналогиях c квантовой теорией поля. Краевая задача для неоднородных полей деформирования в дифференциальной форме преобразуется в интегральное уравнение для тензора деформаций. Решение интегрального уравнения строится в виде ряда по интенсивности взаимодействия деформаций. Это позволяет интерпретировать неоднородную деформацию в какой-либо точке зерна как суперпозицию макродеформации, обусловленной граничными условиями, и двух составляющих, обусловленных внутризеренным и межзеренным взаимодействием. Показано, что в нетекстурированных поликристаллах, несмотря на дальнодействующий характер упругого взаимодействия, при оценке влияния межзеренного взаимодействия на неоднородность деформаций в выделенном зерне можно ограничиться учетом взаимодействия только с ближайшими и вторыми по удалению зернами-соседями. Вклады от взаимодействия с более далекими зернами взаимно компенсируют друг друга. Неоднородное в пределах одного зерна поле деформаций аппроксимируется кусочно-постоянной функцией. Для этого каждое зерно разбивается на большое число малых субзерен, в пределах которых поля деформаций принимаются однородными. Такая аппроксимация сводит интегральные уравнения для локальных деформаций к линейным алгебраическим, которые решаются численно. Применение метода к классической задаче вычисления деформаций в сферическом включении, погруженном в неограниченную матрицу, дает решение Эшелби. На модельных поликристаллах цинка выполнена численная оценка неоднородных деформаций. Вблизи границ сферического зерна экстремальные значения деформаций, обусловленные межзеренным взаимодействием, на 30 % превосходят средние. В материалах с более низкой упругой симметрией зерен концентрация деформаций существенно выше.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2018;(1):58-72
Математическая модель импульсного сканирования давления по длине пьезоэлектролюминесцентного оптоволоконного датчика
Аннотация
Разработана математическая модель локации неоднородностей давления по длине оптоволоконного пьезоэлектролюминесцентного датчика с использованием локационного сканирующего электрического видеоимпульса с пошаговым изменением его величины. Разработан алгоритм нахождения функции распределения давления по локальному участку и по всей длине датчика по результатам замеряемой на торцевом сечении датчика интенсивности исходящего из оптоволокна света для случая нелинейной «функции свечения» - зависимости интенсивности света от действующего на электролюминесцентный элемент электрического напряжения; задача сведена к решению интегрального уравнения Фредгольма 1-го рода с разностным ядром, зависящим от управляющего и информативного передаточных коэффициентов датчика и заданной функции свечения электролюминесцентного элемента. Получены аналитические решения для функций распределения давления по длине датчика для частных случаев, когда ядро или сама функция плотности распределения выражаются через дельта-функцию и интегральное уравнение Фредгольма сводится к алгебраическим. Определены области допустимых значений управляющего напряжения датчика для различных режимов диагностики распределения давления. Представлены результаты численных решений прямой и обратной задач для неоднородного распределения давления посредством «точечного» сканирования этого давления предельно узким импульсом управляющего электрического напряжения. В прямой задаче найдены функции свечения на выходе из оптоволокна для различных моментов времени и значений величины импульса управляющего электрического напряжения с учетом заданной функции свечения электролюминесцентного элемента; в обратной задаче найдено распределение давления по значениям функции интенсивности свечения в различные моменты времени.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2018;(1):73-82
Модель процесса кристаллизации с учетом смены фаз при формировании металлического материала методом лазерного сплавления
Аннотация
Разрабатываемые теоретические основы моделей кристаллизующихся тел применяются в настоящий момент для моделирования технологического процесса селективного лазерного сплавления в части описания взаимодействия расплав-твердое тело и определения коробления тела вследствие появления остаточных напряжений в остывающей заготовке. При этом одной из актуальных задач области исследования является построение определяющих соотношений, позволяющих описывать взаимодействие расплав-твердое тело. В данной работе рассмотрено применение известных определяющих соотношений вязкоупругого растущего тела для затвердевающего металлического сплава. При этом процесс кристаллизации протекает в широком температурном диапазоне и сопровождается структурными изменениями. Под процессом кристаллизации понимается переход металлического материала из жидкого состояния в твердое. В рамках исследования рассмотрена общая постановка краевой задачи механики кристаллизующегося тела. Реализация определяющих соотношений, описывающих процесс кристаллизации, выполнена на двух задачах: задаче о неравномерном управляемом остывании стержня и задаче об остывании пластины, начиная с температуры выше точки плавления исследуемого материала. Таким образом, для одномерной и плоской постановок получены численные модели кристаллизации изотропных тел с учетом фазовых переходов на основе метода конечных элементов. Результаты применения метода конечных элементов на основе полученных новых определяющих соотношений не противоречат физике процесса кристаллизации и могут быть применены при моделировании процесса селективного лазерного сплавления металлических материалов с учетом смены фаз. Для каждой модельной задачи получены картины напряженно-деформированного состояния конструкции в целом: поля распределения перемещений, деформаций и напряжений. Также в рамках исследования выполнен анализ сходимости численного решения, проанализировано выполнение естественных граничных условий, получены поля температур и степени кристаллизации.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2018;(1):83-92
ОЦЕНКА ЖЕСТКОСТИ РАЗВЕРТЫВАЕМОЙ ВНУТРЕННИМ ДАВЛЕНИЕМ ЦИЛИНДРИЧЕСКОЙ КОМПОЗИТНОЙ ОБОЛОЧКИ НА НАЧАЛЬНОМ ЭТАПЕ ПОЛИМЕРИЗАЦИИ СВЯЗУЮЩЕГО
Аннотация
Рассматривается экспериментально-расчетный метод оценки жесткостных свойств цилиндрической оболочечной конструкции из композиционного материала с полимерной матрицей на начальном этапе ее отверждения. Оценка жесткостных параметров композиции на этом этапе полимеризации вследствие различного физического состояния армирующих элементов и связующего методами механики композитных материалов приводит к плохо обусловленным матрицам жесткости, непригодным для достоверного описания механического поведения конструкции. Актуальность проводимого исследования связана с изучением технологии изготовления крупногабаритных пневматических конструкций на основе композиций, подвергаемых отверждению в условиях космоса. В предлагаемом методе экспериментально определяется давление развертывания цилиндра (давление, при котором диаметр цилиндра принимает номинальное значение), соответствующее текущей степени полимеризации связующего. Степень полимеризации характеризуется измеряемыми визкозиметром вязкостью и динамическим модулем полимера. Приводится описание установки, используемой в экспериментах, приборов, осуществляющих фиксацию измеряемых параметров состояния, и порядок проведения опытов. Расчетным путем на основе геометрически нелинейной упругой модели методом последовательных приближений определяется эффективный модуль упругости материала цилиндра, соответствующий его жесткостным характеристикам в течение нагружения внутренним давлением до давления развертывания. Время развертывания мало по сравнению с временем полной полимеризации связующего. Посредством сопоставления экспериментальных и расчетных данных устанавливается зависимость эффективного модуля упругости от параметров отверждения связующего. Выявлена практически линейная зависимость давления развертывания цилиндра от эффективного модуля упругости, что позволяет экстраполировать результаты исследования на значения параметров связующего, не подкрепленные опытом. Приведенные результаты позволяют проводить оценку внутреннего давления, необходимого для развертывания цилиндрических оболочек из композитов с частично отвержденным связующим, путем решения задач механики деформируемого твердого тела.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2018;(1):93-99
Моделирование образования новых материальных поверхностей в процессах адгезионного расслоения композита
Аннотация
Построена модель адгезионного расслоения композиционного материала. Расслоение рассматривается как термомеханический процесс, в котором напряженное состояние одной из границ адгезионного слоя изменяется до нуля при нарушении связей с соединяемым телом. В результате расслоения прекращается взаимодействие между частью композита, включающей адгезивный слой, и остальным телом. Получена система двух вариационных условий равновесности в скоростях, описывающая докритическое деформирование и процесс расслоения. Осреднение напряженно-деформированного состояния (НДС) адгезионного слоя позволяет избежать сингулярности в тупиковой точке образуемого математического разреза, и продвижение по границам слоя поверхности разрыва не приводит к ее образованию. При решении задачи докритического деформирования выделяется малая δ-поверхность на границе адгезива, где достигнут критерий отслоения. Распределение нагрузки (узловых сил) на δ-поверхности определяется путем повторного решения задачи докритического деформирования с известным из первоначального решения законом движения границы адгезивного слоя. Решается задача при простой разгрузке δ-поверхности тела и сохранении внешней нагрузки, соответствующей началу процесса расслоения. В результате НДС тела в момент начала локальной разгрузки отличается от его состояния по окончании δ-разгрузки. Для линейно-упругого материала проведено сравнение решений задачи в рамках данной модели и предложенной ранее модели когезионного расслоения, в которой предполагается полное разрушение связующего слоя. Установлено существенное различие граничных перемещений основных слоев композита в процессе его разрушения при увеличении поверхности разрыва контакта адгезионного слоя и основного материала.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2018;(1):100-109
Прикладная модель измельчения твердой частицы простой формы ударом о жесткую поверхность
Аннотация
Сложность явлений, происходящих при измельчении и разрушении твердых тел, затрудняет теоретическое описание этого процесса. В связи с этим представляется актуальным установление взаимосвязей между параметрами, определяющими показатели процесса измельчения, выявление степени их влияния друг на друга, создание и анализ модели процесса измельчения с учетом определяющих его параметров измельчителя и физико-механических свойств материала и на этой основе совершенствование формализованных методов расчета и обоснования рациональных параметров измельчителей, обеспечивающих эффективность их применения при эксплуатации. На основе анализа состояния реального материала большой группой ученых был создан ряд теорий, объясняющих условия и механизмы разрушения в твердых материалах. Однако практическое применение существующих теорий для расчета процессов измельчения довольно затруднительно. Поэтому до сих пор существует необходимость разработки новой простой и удобной для практического применения теории. Авторами предлагается новый способ теоретического описания процесса разрушения материалов. На основе упрощенной энергетической гипотезы и прикладной технической теории распространения волн в упругой сплошной среде было получено новое уточненное решение фундаментальной динамической задачи механики упругодеформируемого твердого тела о продольном соударении бруса постоянного произвольного поперечного сечения, моделирующего частицу материала, с абсолютно жесткой поверхностью, моделирующей рабочий орган измельчителя, учитывающее параметр времени и линейный размер движущегося стержневого элемента (частицы). Разработанная уточненная механико-математическая модель, доведенная до практически применимых расчетных аналитических зависимостей и проиллюстрированная характерными численными примерами, позволяет количественно оценивать прочность частицы твердого вещества в процессе ее разрушения и измельчения, дает возможность осуществить комплексный подход к описанию динамического процесса разрушения частиц материала путем регулирования и подбора оптимальных физико-геометрических характеристик, обеспечивающих требуемое качество измельчения, а также прогнозировать протекание разрушения частицы в зависимости от параметров этого процесса.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2018;(1):110-120