№ 1 (2013)
- Год: 2013
- Статей: 18
- URL: https://ered.pstu.ru/index.php/mechanics/issue/view/33
- DOI: https://doi.org/10.15593/perm.mech/2013.1
О механизмах возбуждения автоколебаний в потоке газа подводом теплоты
Аннотация
Установлены механизмы возбуждения термоакустических автоколебаний, связанные со своеобразной перестройкой структуры тензоров вязких напряжений и диссипации тепловой энергии, в случае когда к движущемуся газу локально подводится теплота. В одномерном потоке тензор диссипации тепловой энергии определяет тепловое сопротивление, для которого при произвольном политропном подводе теплоты найдено аналитическое представление. Это позволило определить области образования отрицательного теплового сопротивления в зависимости от показателя политропы, что составляет новый механизм возбуждения термоакустических автоколебаний. Б.В. Раушенбахом было предложено аппроксимировать зону локального подвода теплоты некоторой поверхностью, называемой поверхностью теплоподвода. Данная поверхность является поверхностью разрывов первого рода гидродинамических и термодинамических параметров сплошной среды. В данной работе найдены граничные условия, которые должны выполняться на поверхности теплоподвода. Эти граничные условия являются следствиями выполнения на поверхности теплоподвода основных законов сохранения: массы, импульса массы и энергии и могут рассматриваться как уравнения для определения компонент тензора диссипации тепловой энергии, который является ассоциированным с поверхностью теплоподвода. Для описания рассматриваемых термоакустических автоколебаний получена математическая модель в виде диссипативной распределенной динамической системы. В качестве примера применения данной математической модели была рассмотрена задача построения продольных термоакустических автоколебаний в трубе Рийке при расположении источника теплоты непосредственно на ее входе. Показано, что причиной возбуждения автоколебаний феномена Рийке в этом случае является N-образная зависимость напряжения силы вязкостного трения по длине трубы как функции средней по сечению трубы скорости движения нагретого воздуха.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2013;(1):7-24
Оценка напряженно-деформированного состояния вращающегося длинного цилиндра
Аннотация
Объект исследования – длинный полый цилиндр, находящийся в сложных условиях деформирования под действием собственного веса и вращения крутящим моментом. Исследуется напряженное состояние по сечению цилиндра и кривизна его продольной оси. В качестве основных допущений, принимаемых для решения рассматриваемой задачи, рассматриваются обычные для механики материалов и инженерных приближений гипотеза о линейности физических соотношений между напряжениями и деформациями (линейная теория упругости) и предположение о малых деформациях. Решение задачи строится на основе дифференциального уравнения упругого изгиба центральной оси цилиндра и основных соотношений между кривизной этой оси, приложенными нагрузками, деформациями и напряжениями по сечению цилиндра. Граничная задача для обыкновенного дифференциального уравнения четвертого порядка с соответствующими граничными условиями решается с помощью метода вариаций произвольных постоянных. Это позволило получить точное решение задачи изгиба длинного цилиндра, вращающегося вокруг продольной оси. Это позволило определить зависимость кривизны центральной оси цилиндра от продольной координаты и найти характер распределения напряжения по сечению цилиндра. Выполнена оценка вклада в напряженное состояние от каждого из факторов, действующих на рассмотренное изделие. Решение поставленной задачи позволило определить эквивалентное напряжение в периферийных слоях цилиндра как результат воздействия всех рассмотренных силовых факторов.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2013;(1):25-38
Применение методов оптимизации к численному анализу температурных полей заготовок кварцевого оптического волокна
Аннотация
Предложена методика определения величины теплового нагружения при высоких температурах, позволяющая с высокой точностью найти температуру нагрева поверхности пламенем газовой горелки, что необходимо для обеспечения точности дальнейших исследований температурных полей и напряженно-деформированного состояния изделий. Рассматривается задача определения температуры поверхности заготовки оптического волокна из кварцевого стекла в зоне разогрева газовой горелкой и дальнейшее исследование температурных полей в заготовке в процессе производства. Методика определения температуры поверхности заготовки в зоне нагревания газовой горелкой включает проведение эксперимента, позволяющего измерить температуру в нескольких точках поверхности заготовки, удаленных от пламени горелки, и последующее решение задачи оптимизации для нахождения температуры поверхности непосредственно в зоне нагрева. В процессе эксперимента температура поверхности трубки измерялась двумя термопарами, расположенными на расстоянии от зоны нагрева, а затем решалась задача оптимизации. В качестве целевой функции выбрана сумма квадратов отклонений расчетных значений температуры в заданных точках от значений, полученных в ходе эксперимента. Разработанная методика проверена на ряде тестовых задач. Единственность решения задачи оптимизации была подтверждена путем «спуска» из нескольких различных начальных значений. Исследованы температурные поля в заготовке кварцевого оптического волокна, нагреваемой с помощью равномерно движущейся газовой горелки. Для решения нестационарной задачи теплопроводности использован метод конечных элементов, задача оптимизации решалась методом золотого сечения.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2013;(1):39-51
Сверхмногоцикловое усталостное разрушение титановых дисков компрессора
Аннотация
Описана процедура расчета напряженного состояния диска компрессора газотурбинного двигателя в полетном цикле нагружения (режим малоцикловой усталости – МЦУ) и при низкоамплитудных вибрациях лопаток (режим сверхмногоцикловой усталости – СВМУ). Исследуются критерии и модели многоосного разрушения в условиях малоцикловой усталости и в условиях сверхмногоцикловой усталости. Определены параметры моделей на основе экспериментальных данных одноосных усталостных испытаний при различных коэффициентах асимметрии цикла. С использованием рассчитанного напряженного состояния с помощью моделей многоосного усталостного разрушения получены оценки долговечности диска компрессора для альтернативных механизмов усталости МЦУ и СВМУ. С целью подтверждения полученных результатов решены две задачи теории упругости для кольцевого диска, моделирующие центробежные нагрузки от лопаток (аналог режима МЦУ) и нагрузки от кручения лопаток (аналог режима СВМУ). Оценки долговечности также указывают, что усталостное разрушение в режимах МЦУ и СВМУ может происходить за одинаковый период реального времени.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2013;(1):52-67
Теоретическое моделирование работы атомно-силового микроскопа при исследовании поверхностей со сложной наноструктурой
Аннотация
Стандартное математическое обеспечение, поставляемое для расшифровки результатов атомно-силового сканирования (АСМ), базируется в основном на моделях, использующих классическое решение задачи Герца о контакте двух линейно-упругих сфер (или сферы и плоского полупространства, если одна из них имеет бесконечно большой радиус). В большинстве случаев этого вполне достаточно. Однако существуют такие ситуации, когда решение Герца следует применять с большой осторожностью. Теоретическому исследованию этих вариантов и посвящена данная работа. Представлены результаты численного моделирования контактного взаимодействия зонда атомно-силового микроскопа и поверхности со сложной наноструктурой. Исследования вели для двух типов материалов: 1) упругая анизотропная среда (зубная эмаль); 2) нелинейно-упругий конечно-деформируемый полимер. Это два класса материалов, которые принципиально различаются по своему механическому поведению. Соответственно, для правильной расшифровки экспериментальных данных требуются разные теоретические модели. Для материалов первого типа решена задача внедрения зонда АСМ в трансверсально-изотропную упругую поверхность. Построены расчетные зависимости силы реакции на инденторе от глубины вдавливания и степени анизотропии материала. Для материалов второго типа (эластомеров) проведено компьютерное моделирование контактного взаимодействия зонда АСМ с эластомерными нанотяжами, которые могут образовываться в полимере в вершине трещины. Тяж представлялся в виде длинной продольной выпуклости, лежащей на плоской упругой поверхности. Задача решалась в трехмерной постановке.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2013;(1):68-80
Численный алгоритм решения задачи о больших упругопластических деформациях МКЭ
Аннотация
Работа посвящена методике исследования конечных упругопластических деформаций. В качестве тензоров, описывающих деформацию и скорость деформации, используются левый тензор Коши–Грина, тензор пространственного градиента скорости и тензор деформации скорости. Вводится удельная потенциальная энергия деформации, которая зависит от левого тензора Коши–Грина. Рассмотрен изотропный материал. Напряженное состояние описывается тензором истинных напряжений Коши–Эйлера, который определяется в актуальном состоянии. Получены линеаризованные определяющие соотношения упругого деформирования в виде зависимости производной Трузделла тензора напряжении Коши–Эйлера от деформации скорости. В рамках теории течения используются аддитивное представление для полной деформации скорости. Предполагается справедливость ассоциированного закона течения. Критерием упругого деформирования является условие Мизеса–Губера. Алгоритм исследования основан на методе последовательных нагружений. В качестве базового уравнения принимается уравнение мощностей в актуальном состоянии. После линеаризации получена разрешающая система линейных уравнений, где неизвестным является приращение перемещений в текущем временном слое. При моделировании пластических деформаций применяется метод проецирования напряжений на поверхность текучести с итерационным уточнением текущего напряженно-деформированного состояния, основанным на введение в разрешающие уравнение мощности дополнительных напряжений. В качестве примера рассмотрено построение алгоритма решения для материала второго порядка. Выбрано соответствующее выражение потенциала упругих деформаций, критерием пластического течения служит условие Губера–Мизеса с изотропным упрочнением. Получены линеаризированные определяющие соотношения. Численная реализация основана на методе конечных элементов. Используется восьмиузловой конечный элемент. Созданный алгоритм исследования больших упругопластических деформаций опробован на решении тестовой задачи о растяжении круглого стержня с образованием шейки. Приводятся результаты решения и сравнение с результатами, полученными другими авторами. Также исследовалось деформирование квадратной плиты под действием внутреннего давления.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2013;(1):81-93
Исследование эффекта Баушингера и границы текучести при упругопластическом деформировании металлов
Аннотация
Представлены результаты опыта на многократное знакопеременное нагружение при растяжении-сжатии тонкостенного трубчатого образца из стали 45 на автоматизированном испытательном комплексе на сложное нагружение СН-ЭВМ. Дана экспериментальная оценка эффекта Баушингера и влияния на него различных допусков на остаточную деформацию. Исследуется влияние допуска на остаточную деформацию при определении вторичных пределов текучести, а также на радиус-функцию и положение центра гипотетической поверхности текучести в пространстве напряжений, используемой в теории пластического течения. По результатам проведенного исследования установлено, что с ростом допуска на остаточную деформацию параметр, характеризующий эффект Баушингера, и радиус-функция поверхности текучести увеличиваются, а смещение ее центра уменьшается. Показано, что параметр, характеризующий эффект Баушингера, уменьшается с ростом длины дуги пластического деформирования и стремится к некоторому стационарному значению. Экспериментально установлено, что для стали 45 смещение центра поверхности текучести увеличивается с ростом длины дуги пластического деформирования. Также установлено, что радиус сферической поверхности текучести совершает временное понижение, а затем с ростом длины дуги пластического деформирования увеличивается. При этом в некоторых математических моделях теории пластического течения временное убывание радиуса поверхности текучести (сужение поверхности текучести) в начале процесса пластического деформирования не связывают с определением пределов текучести по допуску на остаточную деформацию.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2013;(1):94-105
Оптимизация прессового инструмента при производстве низкотемпературных сверхпроводников
Аннотация
Предложена методика определения оптимальных углов технологического инструмента при прессовании сборных композиционных заготовок. В основу оптимизации положено напряжение прессования, обеспечивающее минимальные энергозатраты в ходе процесса. Актуальной является задача совершенствования и оптимизации технологии производства низкотемпературных композиционных сверхпроводников для достижения необходимых производственных объемов выпуска при соответствии требованиям качества. Важным является разработка теоретических основ и методик проектирования технологических процессов применительно к технологии производства низкотемпературных композиционных сверхпроводников, научно обосновывающих выбор технологических режимов и технологической оснастки для повышения качества низкотемпературных композиционных сверхпроводников и обеспечения высоких технико-экономических показателей их производства. Процесс прессования находит широкое применение при обработке металлов давлением. Сущность процесса прессования заключается в выдавливании материала, помещенного в замкнутый объем, через канал, образованный прессовым инструментом. Достоинством процесса прессования является благоприятная схема напряженного состояния с преобладающим влиянием сжимающих напряжений, обеспечивающих повышенную пластичность прессуемого материала. Поэтому процесс прессования широко используется при обработке давлением малопластичных труднодеформируемых металлов и сплавов.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2013;(1):106-120
Об устойчивости конвективных течений при изотермическом диффузионном смешении трехкомпонентных газовых смесей
Аннотация
В рамках линейной теории устойчивости проведен анализ экспериментальных данных по изучению неустойчивости механического равновесия при диффузионном смешении трехкомпонентных газовых смесей в изотермических условиях. Рассматривается канал в виде вертикального цилиндра конечной высоты при граничных условиях, предполагающих отсутствие переноса компонентов через стенки диффузионного канала. Результаты теории сравниваются с полученным в опытах положением границ устойчивости для системы 0,4722 He + 0,5278 Ar – N2 при варьировании диаметра канала. Показано, что в условиях развитой конвекции существование максимумов интенсивности переноса компонентов в зависимости от термодинамических параметров связано с модой возмущений, определяющей различные виды конвективных течений. Структурные образования, движущиеся навстречу друг другу, число которых в поперечном сечении диффузионного канала определяется модой возмущения, взаимодействуют и тем самым определяют волнообразное изменение интенсивности неустойчивого процесса. Для рассматриваемой системы первый максимум интенсивности неустойчивого процесса в зависимости от диаметра, возникающий при d = 6 мм, характеризуется модой возмущений n = 3 и критическим числом Рэлея R3 = 972,7. При этих условиях происходит стабилизация конвективного переноса, характеризующаяся определенным типом движения конвективных формирований, т.е. движением по шести каналам (три тока и три противотока). Также результаты показали, что увеличение моды возмущений и диаметра диффузионного канала приводит к нелинейному увеличению критического числа Рэлея для всех компонентов, участвующих в переносе. Приведенные результаты хорошо согласуются с экспериментальными данными для системы 0,4722 He + 0,5278 Ar – N2.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2013;(1):121-135
Параметры процесса механической обработки с учетом температурных режимов
Аннотация
В работе рассмотрены основные процессы механической обработки и их влияние на прочность цилиндрических изделий, рассчитаны температурные режимы при резании труб, параметры процесса резания выбраны с учетом величины и знака температурных остаточных напряжений. Под механической обработкой понимают различные виды резания: точение, сверление, фрезерование, шлифование, хонингование, полирование. В связи с этим показано, как процессы механической обработки влияют на качество трубных изделий. Целью обработки резанием является придание деталям желаемой формы. Также рассмотрены параметры процесса резания, выбранные с учетом величины и знака температурных остаточных напряжений. Остаточные напряжения при резании металлов образуются в результате неравномерности пластической деформации и значительного нагрева поверхностных слоев. Решая задачу термоупругости, можно определить напряжения, возникающие при нагреве. Исходя из условия не наступления пластических деформаций на поверхности детали после прохождения инструмента можно оптимизировать параметры резания при обработке различных осесимметричных металлоизделий. Таким образом, зная соотношения между этими параметрами, а также механические свойства материала заготовки и геометрию, можно оценить величину интенсивности напряжений в зависимости от начальной температуры и параметров резания при обработке трубы. Далее, на основе полученных результатов построены соотношения и зависимости важнейших параметров резания для трубных изделий из различных материалов и представлены предельные соотношения параметров механической обработки трубы. В работе предложена методика определения предельных соотношений параметров процесса резания для осесимметричных изделий с учетом температурного разогрева поверхности за счет трения при обточке. Методика позволяет определять напряжения на поверхности заготовки и рационально подбирать параметры процесса в зависимости от геометрии детали и обрабатываемого материала.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2013;(1):136-147
Эффективность использования идеологии мастер-процесса при параллельной реализации итерационных процедур решения линейных систем
Аннотация
В работе рассматривается эффективность использования идеологии мастер-процесса при построении параллельной реализации процедуры решения систем линейных алгебраических уравнений. Исследования проведены применительно к итерационным алгоритмам методов сопряженных градиентов и Якоби. Использовался разреженный формат RR(C)U хранения матрицы коэффициентов алгебраической системы. Исследования проводились применительно к системе линейных алгебраических уравнений, сформированной в ходе решения двумерной краевой задачи линейной теории упругости методом конечных элементов. В качестве критерия количественного анализа использовалось «ускорение», равное отношению времени выполнения последовательного алгоритма ко времени выполнения параллельного алгоритма. Анализ проведен для алгебраических систем от 100 до 50000 уравнений с использованием ЭВМ на базе шестиядерного процессора AMD® Phenom II X6 1075T. Программная реализация итерационных алгоритмов решения системы линейных алгебраических уравнений была произведена на языке программирования C#. Для обеспечения взаимодействия параллельных процессов использовался стандарт MPI 2.0. Полученные данные позволяют сделать вывод, что использование мастер-процесса при параллельной реализации метода сопряженных градиентов приводит к незначительному, от 10 до 15 процентов, снижению величины коэффициента ускорения по отношению к варианту реализации, не использующему мастер-процесс. В случае использования метода Якоби эффект выражен намного слабее. Принимая во внимание структурное удобство использования идеологии мастер-процесса, можно говорить о допустимости применения данного архитектурного решения при реализации рассмотренных методов.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2013;(1):148-160
Моделирование напряженно-деформированного состояния тонких композитных покрытий на основе решения плоской задачи градиентной теории упругости для слоя
Аннотация
В работе получена и исследована постановка задачи плоской деформации в рамках градиентной модели упругости для композитного слоя. На основе синтеза аналитических и численных методов решена плоская задача градиентной теории упругости для бесконечного слоя с целью учета распределения напряжений в плоскости покрытий для сверхтонких структур. Исследована задача о воздействии распределенной поверхностной сжимающей нормальной нагрузки на бесконечный слой, находящийся на жестком основании. Для решения задачи используется интегральное преобразование Фурье, при этом обратное преобразование вычисляется с использованием численной процедуры. В работе показано, что предложенные модели позволяют прогнозировать эффекты локализации напряжений в окрестности межслойных зон в покрытии и учитывать влияние неклассических масштабных факторов – толщины слоев покрытия и градиентных параметров моделей. Для моделирования привлекается наиболее простой вариант градиентной теории упругости – прикладная модель межфазного слоя, содержащая единственный дополнительный физический параметр, определяющий «градиентность» среды и протяженность межфазных зон в области границ материала. Этот параметр является дополнительной физической константой, характеризующей контакт разнородных материалов. Для численных вычислений в работе используются гипотетические значения градиентного параметра. Прикладное значение решенной задачи связано с возможностью достоверного моделирования и оптимизации микроструктурного строения ультратонких защитных композитных покрытий, применяемых в авиакосмической отрасли. Также построенное решение может быть использовано для идентификации дополнительных физических параметров градиентной теории упругости на основе сопоставления результатов моделирования и экспериментальных данных по индентированию тонкослойных структур.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2013;(1):161-181
Смешанные дискретные модели в анализе упругих трехмерных неоднородных тел сложной формы
Аннотация
Предложена процедура построения смешанных дискретных моделей для анализа напряженного состояния трехмерных упругих тел, конструкций сложной формы, имеющих неоднородную (композитную) структуру. Смешанные модели состоят из однородных односеточных трехмерных конечных элементов (КЭ) первого порядка формы куба и двухсеточных конечных элементов (ДвКЭ) неоднородной структуры формы прямоугольного параллелепипеда, т. е. состоят из КЭ различной сеточной структуры. В окрестности крепления тела или сложной формы границы используем мелкое разбиение, состоящее из односеточных КЭ и учитывающее неоднородную структуру и сложную форму тела, в остальной части тела – крупное, представленное ДвКЭ. Мелкое и крупное разбиения склеиваем с помощью связующих КЭ, построенных на основе ДвКЭ. Для построения ДвКЭ используем две вложенные сетки: мелкую и крупную. Область ДвКЭ представляем базовым (мелким) разбиением на КЭ первого порядка, которое учитывает его неоднородную структуру и порождает мелкую сетку. На базовом разбиении определяем в матричной форме функционал полной потенциальной энергии ДвКЭ, который (с помощью аппроксимаций, построенных на крупной сетке) проецируем на крупную сетку. Из условия минимизации полученного функционала по узловым перемещениям крупной сетки находим формулы для вычисления матрицы жесткости и вектора узловых сил ДвКЭ. Достоинства ДвКЭ состоят в том, что с помощью базового (мелкого) разбиения учитывается неоднородная структура ДвКЭ, они образуют дискретные модели малой размерности и порождают решения с заданной погрешностью. Погрешность решения варьируется с помощью соотношения шагов мелкой и крупной сеток ДвКЭ. Достоинства смешанных дискретных моделей состоят в том, что они имеют малую размерность, учитывают сложную форму тел, неоднородную структуру и порождают решения с заданной погрешностью. Приведен пример расчета.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2013;(1):182-195
Моделирование волн на поверхности цилиндрической конфигурации магнитной жидкости, окружающей длинное пористое ядро
Аннотация
Построена и исследована математическая модель распространения и неустойчивости волн на поверхности цилиндрического столба магнитной жидкости бесконечной длины, окружающей коаксиально расположенное, бесконечно длинное ядро (из пористого материала) круглого сечения. Найдены условия, при которых возмущения поверхности жидкого столба становятся неустойчивыми и приводят к его распаду на цепочку из соединенных капель. Учитывается наличие поверхностного натяжения. Сила тяжести предполагается отсутствующей. Ось пористого цилиндра совпадает с осью коаксиально расположенного соленоида, создающего однородное магнитное поле. Задача решается в цилиндрической системе координат ( r , θ, z ), в которой жидкий столб покоится. Ось z направлена по оси соленоида. Записаны уравнения движения магнитной жидкости внутри и вне пористой среды, а также уравнения для магнитного поля в пористой среде, жидкости и воздушном зазоре. Сформулированы граничные условия для гидродинамических и магнитных величин на поверхностях раздела сред. Возмущенное (в связи с распространением волны) магнитное поле ищется внутри и вне пористой среды, а также в воздушном зазоре соленоида. Найдено полное решение краевой задачи для гидродинамических и магнитных величин. Проведен численный анализ полученного дисперсионного уравнения, описывающего распространение поверхностных волн. Рассмотрены различные частные случаи. Найдены условия, при которых возмущения поверхности жидкого столба устойчивы (затухающие волны) либо неустойчивы (что приводит к нарастанию возмущений и распаду цилиндра на цепочку капель). Показано, что размер образующихся при распаде капель увеличивается с ростом магнитного поля, т.е. магнитное поле оказывает стабилизирующее влияние на распад жидкого столба.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2013;(1):196-209
Эволюция полей технологических напряжений в цилиндрическом силовом стержне для заготовки оптоволокна типа Panda в процессе отжига
Аннотация
Исследованы закономерности эволюции полей технологических напряжений в цилиндрическом силовом стержне для заготовки оптоволокна типа Panda в процессе отжига. Реализована математическая модель формирования технологических напряжений, построенная на основе соотношений линейной термовязкоупругости, в температурных режимах, соответствующих технологическому процессу отжига силового стержня. Установлены количественные характеристики релаксации напряжений в различных условиях протекания процесса отжига. Выявлены режимы, приводящие к максимально возможному снижению опасных растягивающих нормальных напряжений и интенсивности тензора напряжений. Установлено, что для стержней, закон легирования в которых близок к равномерному, процесс отжига является более эффективным, чем для стержней с законами легирования, близкими к параболическому.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2013;(1):210-220
Конечно-элементная реализация метода геометрического погружения на основе вариационного принципа Кастильяно для плоской задачи теории упругости
Аннотация
Рассмотрен вариант метода геометрического погружения для плоских задач теории упругости, основанного на методе конечных элементов в напряжениях в рамках принципа минимума дополнительной работы упругой системы. Суть метода геометрического погружения заключается в сведении исходной задачи для линейно-упругого тела произвольной формы к итерационной последовательности задач теории упругости на некоторой канонической области. Сформулирована итерационная процедура для решения вариационного уравнения метода геометрического погружения, а также процедура построения его дискретного аналога с помощью метода конечных элементов в напряжениях для плоской задачи теории упругости в декартовой системе координат. Использован вариант конечного элемента в терминах функции напряжений для удовлетворения аппроксимирующих выражений уравнениям равновесия. Продемонстрировано практическое применение метода на примере решения плоской задачи для упругой пластины с прямоугольным вырезом. Получено достаточно хорошее соответствие результатов определения полей напряжений в сравнении с традиционным методом конечных элементов в перемещениях. Установлена практическая сходимость итерационной процедуры метода геометрического погружения. Уделено внимание способам задания статических граничных условий, являющихся главными для данной вариационной формулировки. Использован способ модификации матрицы податливости системы конечных элементов и метод множителей Лагранжа.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2013;(1):221-234
Метод вычисления поверхностной энергии и энергии адгезии упругих тел
Аннотация
Существует ряд задач, для успешного решения которых необходимо знать поверхностную энергию, энергию и силы адгезии и когезии. Примером является задача о расслоении композитов. В выражение энергетического критерия прочности входит суммарная поверхностная энергия контактирующих элементов или энергия адгезии. Эти величины определить методами классической теории упругости невозможно. Силовые критерии прочности используют значения предельных напряжений – сил адгезии и когезии. Обычно их определяют экспериментально, что не всегда возможно или экономически невыгодно. В работе предложен метод расчета поверхностной энергии и энергии адгезии упругих тел, находящихся в состоянии адгезии. Учтено следующее. Распределенная по границе тела суммарная поверхностная энергия равна распределенному по его объему изменению свободной энергии, произошедшему при образовании границы. При адгезии двух тел вдоль поверхности контакта формируется переходный слой. В нем физические и термодинамические свойства одного тела непрерывно переходят в свойства другого тела. Метод базируется на варианте градиентной модели сплошной упругой среды. В ее основе лежит предположение о многочастичном потенциальном нелокальном взаимодействии бесконечно малых частиц, составляющих среду. Дополнительные к классическим характеристики упругого состояния вычисляются с помощью дифференцирования известного выражения объемной плотности распределения свободной энергии. Оно строится на основании дополнительных гипотез о составе среды: упругая часть, фононный и электронный (для металлов) газы.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2013;(1):235-259
К теории процесса гидратообразования газовых пузырей в условиях Мирового океана
Аннотация
В работе предложена и построена теоретическая модель миграции метановых пузырей в условиях образования гидрата в вертикальном канале, находящегося над газовым источником. Изучено влияние различных параметров на процесс образования гидратных пузырей. Получены минимальные значения массовых расходов газа и воды, необходимых для процесса образования гидрата. Установлено, что при миграции газовых пузырей в вертикальном канале возможны два режима протекания процесса гидратообразования в зависимости от значения массового расхода воды. Если начальное значение массового расхода воды больше критического, то при данном массовом расходе газа, достаточного для гидратообразования, газовые пузырьки полностью превращаются в гидратные, при этом температура воды в канале не достигает равновесной температуры образования гидрата. Если начальное значение массового расхода воды меньше критического, то газовые пузырьки покрываются гидратной оболочкой, а температура воды в канале достигает равновесной и процесс гидратообразования завершается. Получено, что с увеличением начального массового расхода воды высота всплытия гидратных пузырьков растет до некоторой максимальной величины, которая соответствует его критическому значению, а при дальнейшем увеличении – падает. В работе проведен анализ влияния на процесс гидратообразования различных глубин, на которых расположены газовые источники: 800, 1200 и 1500 м. Получено, что если происходит процесс частичного гидратообразования, то чем глубже расположен источник, тем на большую высоту всплывают пузырьки. Если происходит процесс полного образования гидрата, то картина совершенно иная: чем выше расположен источник газа – тем больше высота гидратообразования.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2013;(1):260-273