№ 2 (2016)
- Год: 2016
- Статей: 11
- URL: https://ered.pstu.ru/index.php/mechanics/issue/view/21
- DOI: https://doi.org/10.15593/perm.mech/2016.2
Моделирование технологических дефектов и оценка их влияния на статическую прочность композитных фланцев
Аннотация
Предложена приближенная методика оценки влияния дефектов на статическую прочность композитных фланцев для высоконагруженных корпусных деталей и узлов авиационной техники из полимерных композиционных материалов (ПКМ), основанная на современных методах численного моделирования. В процессе формования фланцев в зонах перегиба возможно появление различных дефектов: искривление слоев, смоляные карманы, поры, расслоения и т.п. Таким образом, задача оценки влияния технологических дефектов на прочность данного изделия является актуальной. Определены основные типы дефектов в слоистых композитах, а также проведен обзор основных работ, посвященных исследованию вопросов прочности слоистых конструкций из ПКМ с дефектами. С использованием двумерной осесимметричной модели в среде конечно-элементного анализа ANSYS проведено численное моделирование механического поведения стеклопластикового фланца с основными типами дефектов при силовом нагружении. Задача решалась в упругой постановке с учетом анизотропии механических свойств материала. В разработанной модели предполагается параметризация размеров и конфигурации дефектов, а также возможная вариация места расположения дефектов по выбранной подобласти композиционного материала. Для сохранения исходного контура конструкции общая толщина пакета слоев в локальной зоне дефекта уменьшалась на величину дефекта. Таким образом, оценка прочности конструкции с дефектом осуществлялась с запасом. Оценка запаса статической прочности рассматриваемого фланца проводилась по компонентам напряжений с использованием критерия максимальных напряжений. Модель позволяет анализировать как нормальные напряжения в плоскости слоя, так и межслоевые отрывные и сдвиговые напряжения. Сравнение полученных значений с результатами для неповрежденной конструкции фланца позволяет сделать вывод о степени влияния рассматриваемого дефекта или их сочетаний на статическую прочность изделия.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2016;(2):5-21
Особенности численного моделирования упругопластического выпучивания полусферических оболочек при нагружении жестким индентором
Аннотация
Рассматривается трехмерная задача упругопластического деформирования и выпучивания полусферических оболочек при контактном взаимодействии с жесткими телами. Для описания деформирования оболочки применяется текущая лагранжева формулировка. Уравнение движения выводится из баланса виртуальных мощностей. В качестве уравнений состояния используются соотношения теории течения с изотропным упрочнением. Контактное взаимодействие жесткого тела и оболочки моделируется исходя из условия непроникания. Для дискретизации определяющей системы уравнений по пространственным переменным применяется 8-узловой изопараметрический конечный элемент с полилинейными функциями формы. Решение задачи при заданных граничных и начальных условиях основано на моментной схеме метода конечных элементов и явной конечно-разностной схеме интегрирования по времени типа «крест». Выполнены расчеты упругопластического деформирования и выпучивания полусферической стальной оболочки, расположенной на неподвижной плите при квазистатическом вдавливании недеформируемого индентора. Индентор имеет вид цилиндрической оболочки с продольными вырезами в области контакта. Количество вырезов в расчетах варьировалось. Проведено численное исследование влияния формы индентора на уровень пластических деформаций и величину критической нагрузки. Как показали результаты расчетов, полусферическая оболочка в процессе нагружения претерпевает значительные локальные формоизменения в области контакта с индентором, характеризуемые большими смещениями и углами поворота конечных элементов как жесткого целого. Достоверность результатов численного решения задачи подтверждается хорошим соответствием с экспериментальными данными других авторов. Показано, что при малых деформациях поперечного сдвига и больших углах поворота применение моментной схемы метода конечных элементов с определением скоростей деформаций и напряжений в локальном базисе вполне обоснованно при использовании явной конечно-разностной схемы интегрирования типа «крест» с малыми шагами по времени.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2016;(2):22-33
Моделирование процесса разделения композита с адгезионным слоем
Аннотация
Рассматривается модель деформирования композиционного материала с тонким адгезионным слоем. Напряженное состояние слоя предлагается рассматривать на основе связи средних по толщине слоя напряжений и напряжений по границам слоя. Средние деформации слоя выражаются через его граничные перемещения. Использование средних напряжений и деформаций позволяет избежать зависимости напряженно-деформированного состояния слоя от формы его торцевых поверхностей. В рамках малых деформаций получено вариационное условие равновесия тел, соединенных посредством адгезионного слоя. Задача рассматривается в рамках линейной теории упругости. Для сопрягаемых тел законом Гука непосредственно связываются поля деформаций и напряжений. В результате связанная система вариационных уравнений сводится к уравнениям относительно полей перемещений в сопрягаемых телах, в том числе и на границах контакта со слоем. В качестве параметра система вариационных уравнений в перемещениях содержит толщину адгезионного слоя. Существенно, что данная система уравнений не является дискретной, так как поля перемещений полагаются непрерывными. Для получения приближенного решения можно использовать различные аппроксимации перемещений. В частности, применялся метод конечного элемента с квадратичной аппроксимацией полей перемещений для случая плоской деформации. Исследовалось влияние характерного размера конечного элемента на сходимость решения. Установлено, что если отношение грани конечного элемента к толщине слоя равно четырем и более, то имеет место численная сходимость. В силу отсутствия сингулярности напряжений в точках сопряжения адгезионного слоя с телами предлагаемый подход позволяет использовать известные локальные критерии разрушения. Проведен анализ возможных видов разрушения композита как за счет разрушения материала слоя, так и посредством разрыва связей между слоем и смежными материалами.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2016;(2):34-44
НЕОСЕСИММЕТРИЧНАЯ ПОТЕРЯ УСТОЙЧИВОСТИ ПРИ ОСЕСИММЕТРИЧНОМ НАГРЕВЕ КРУГЛОЙ ПЛАСТИНЫ
Аннотация
Представлены результаты простых экспериментов, демонстрирующие неосесимметричную краевую потерю устойчивости с разным числом волн по окружности при равномерном нагревании тонких пластиковых круглых пластин, края которых зафиксированы на жёстких кольцах. При этом формы потери устойчивости с большим числом волн по окружности наблюдались для пластин с меньшей толщиной. Для описания такого типа многоволновой потери устойчивости использована классическая модель устойчивости круглой пластины под действием радиальных сжимающих усилий, распределённых по контуру пластины, которая включает в том числе температурное выпучивание пластины. Получена зависимость минимальной критической нагрузки от числа волн по окружности в форме потери устойчивости пластины. Выполнено сопоставление многоволновой формы потери устойчивости пластины с похожими на неё формами собственных колебаний, локализованными у края пластины. Сопоставление проводилось по положению переходной линии, математически отделяющей область с активной осцилляцией от плато с практически недеформированной центральной областью пластины, смещающейся как жёсткое целое, а также - по расположению экстремумов формы потери устойчивости и формы колебаний относительно центра пластины. По этим же параметрам выполнено сравнение наблюдаемых в эксперименте неосесимметричных форм потери устойчивости пластины с результатами расчётов по теоретической модели при условиях жесткого, шарнирного и упругого закреплений края пластины. В последнем случае показана возможность определения жёсткости опоры исходя из условий совпадения теоретических и экспериментальных значений радиуса переходной линии и радиуса окружности, на которой располагаются экстремумы формы потери устойчивости пластины. Приведенные данные иллюстрируют тенденцию сдвига экстремумов формы потери устойчивости к границе пластины при ослаблении краевых условий.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2016;(2):45-53
Обоснование использования метода конструктивной анизотропии при расчете пологих оболочек двоякой кривизны, ослабленных вырезами
Аннотация
Приводится геометрически нелинейная математическая модель деформирования изотропных пологих оболочек двоякой кривизны, ослабленных вырезами. Модель основывается на гипотезах теории оболочек Кирхгофа-Лява и представлена в виде геометрических соотношений, физических соотношений и функционала полной потенциальной энергии. Также приводятся выражения для усилий и моментов. Рассмотрено два способа введения вырезов: дискретно и методом конструктивной анизотропии, который позволяет наиболее точно «размазать» нулевую жесткость вырезов по полю оболочки. Для минимизации функционала применяется метод Ритца, что сводит задачу к решению системы нелинейных алгебраических уравнений, которая решается методом Ньютона. Алгоритм реализован в среде аналитических вычислений Maple 2015. Проводится анализ устойчивости пологих оболочечных конструкций двоякой кривизны, выполненных из стали, при действии на них внешней равномерно распределенной поперечной нагрузки и шарнирно-неподвижном способе закрепления контура оболочки. Расчеты производились при наличии разного числа вырезов, при этом фиксировался коэффициент отношения общей площади вырезов к площади оболочки. Таким образом, при увеличении числа вырезов уменьшался их размер. Распределение вырезов по оболочке делалось двумя разными способами. Для всех исследованных конструкций приводятся значения критических нагрузок потери устойчивости. Проводится сравнение значений, полученных при дискретном введении вырезов и методом конструктивной анизотропии. Для нескольких вариантов конструкций показаны графики «нагрузка - прогиб». Для одного варианта оболочки, ослабленной большим числом вырезов, приводятся поля прогибов до и после потери устойчивости, цветом также показана интенсивность напряжений. На основании полученных данных показано, что при увеличении числа вырезов дискретность их ввода теряется и становится возможным использование специально разработанного метода конструктивной анизотропии. Таким образом, сделано обоснование использования данного метода при расчете устойчивости пологих оболочек, ослабленных большим числом вырезов.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2016;(2):54-68
Методика расчета механики систем связанных тонких упругих стержней по дифференциальной модели
Аннотация
Предложена методика расчета геометрически нелинейного динамического деформирования системы связанных между собой стержней по дифференциальной модели тонкого упругого стержня. Деформирование каждого стержня описывается 12 дифференциальными уравнениями с краевыми условиями. Краевые условия для каждого стержня формируются из условий закрепления, нагружения и условий стыковки стержней между собой. Полученная нелинейная краевая задача решается численно, для этого требуется многократно решать большие СЛАУ, матрицей которых является матрица Якоби совокупности дифференциальных уравнений, краевых условий и условий стыковки для всех стержней, входящих в систему. На главной диагонали такой матрицы стоят блоки от дифференциальных уравнений и краевых условий для каждого отдельного стержня, и эти блоки связаны между собой в остальной части матрицы элементами, полученными из уравнений связи между стержнями в системе. В статье предлагается метод решения таких блочных разреженных СЛАУ, на главной диагонали которых стоят невырожденные блоки, связанные с другими блоками относительно небольшим количеством уравнений. Решение исходной СЛАУ распадается на решение нескольких систем меньшего размера (для каждого блока) и решения СЛАУ размером, равным числу ненулевых строк, связывающих матричные блоки. Поэтому предлагаемый метод эффективен при относительно небольшом по сравнению с размерами матрицы количестве ненулевых строк, связывающих выделенные матричные блоки, стоящие на главной диагонали. В статье приведены тестовые примеры расчета больших перемещений нагруженных рам, результаты расчета по предложенной методике сравниваются с результатами расчета в конечно-элементном пакете ANSYS. Результаты расчетов по методам практически совпали, точность совпадения зависит от дискретизации моделей.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2016;(2):69-88
Процессы трансформации и взаимодействия микротрещин в металле под воздействием высокоэнергетического импульсного электромагнитного поля
Аннотация
Моделирование процессов, происходящих в металлах при обработке их кратковременными импульсами электрического тока высокой плотности, представляет интерес прежде всего в связи с изучением явления электропластичности, физический механизм которого до сих пор остается непонятным для исследователей. Одно из объяснений этого явления заключается в эффекте залечивания микродефектов в металле. В настоящей работе исследуются процессы трансформации и взаимодействия дефектов типа плоских микротрещин с линейными размерами порядка 10 мкм, протекающие в металле при обработке его импульсами тока. Исследование осуществляется численно на основе связанной модели воздействия интенсивным электромагнитным полем на предварительно поврежденный термоупругопластический материал с дефектами, которая учитывает плавление и испарение металла, а также зависимость всех его физико-механических свойств от температуры. Решение получающейся системы уравнений ищется методом конечных элементов на подвижных сетках с использованием смешанного эйлеро-лагранжева метода. Расчеты показали, что под действием импульсов тока происходит сварка трещины и залечивание микродефекта. Залечивание происходит путем одновременного уменьшения длины, выброса расплавленного металла внутрь трещины и смыкания берегов, что приводит к тому, что берега трещины начинают контактировать со струей расплавленного материала и в конце этих процессов струя оказывается полностью зажатой берегами трещины. В статье исследуется влияние расстояния между трещинами и их взаимного расположения относительно друг друга на процессы деформирования и залечивания микродефектов, а также обсуждается выбор области интегрирования и условий на ее границах. Численное моделирование показало, что при изучении процесса залечивания микротрещин можно без потери точности ограничиться рассмотрением в качестве области интегрирования одного представительного элемента (или одной четверти осесимметричного представительного элемента), задавая на ее границах, не являющихся осями симметрии, разность потенциалов, определенную для элемента без дефекта (в состоянии, «не возмущенном» наличием микротрещины). При расстояниях между трещинами, превышающих 5-6 их длин, процессы залечивания будут протекать одинаково, независимо от того, моделируем мы их в области интегрирования, состоящей из одного или нескольких представительных элементов. С увеличением расстояния между трещинами наблюдается уменьшение влияния взаимного расположения микротрещин на процессы их залечивания, так что при расстояниях между микротрещинами, превышающих шесть их длин, время залечивания микротрещин оказывается фактически одинаковым при любом расположении трещин относительно друг друга. Взаимодействие между микротрещинами начинает заметно сказываться на процессе их залечивания, когда расстояния между ними сокращается до 5-6 длин микротрещин. Если же расстояние между трещинами превышает шесть и более их длин, то процессы залечивания микротрещин, становятся практически независящими как от расстояния между дефектами, так и от расположения дефектов относительно друг друга. Уменьшение расстояния между трещинами вплоть до 1-2 их линейных размеров (с учетом изменения их взаимного расположения) качественно не меняет описанный процесс залечивания, однако приводит к его существенному замедлению: выброс расплавленного материала в трещину сохраняется, но уменьшение трещины особенно в поперечном направлении значительно сокращается.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2016;(2):89-110
КОМПЛЕКСНЫЙ ПОДХОД К ИССЛЕДОВАНИЮ ХАРАКТЕРИСТИК ДИНАМИЧЕСКОГО ДЕФОРМИРОВАНИЯ НА ПОВЕРХНОСТИ НЕЖЕСТКИХ ДОРОЖНЫХ ОДЕЖД С ИСПОЛЬЗОВАНИЕМ МЕТОДОВ НЕРАЗРУШАЮЩЕГО КОНТРОЛЯ
Аннотация
Предложен комплексный подход к исследованию характеристик динамического деформирования на поверхности нежестких дорожных одежд с использованием методов неразрушающего контроля, включающий в себя как анализ расчетных и экспериментальных чаш максимальных динамических прогибов, так и анализ характеристик динамического отклика нежесткой дорожной конструкции в частотной области при кратковременном ударном нагружении. Представлены результаты численного эксперимента по моделированию влияния потери сцепления между отдельными элементами нежестких дорожных одежд на форму амплитудно-частотной характеристики (АЧХ) ускорения на поверхности дорожной одежды. Установлено что при полном межслойном сцеплении между всеми элементами нежесткой дорожной конструкции в амплитудно-частотной характеристике ускорения присутствует единственный частотный экстремум. При потере межслойного сцепления между отдельными элементами нежесткой дорожной конструкции на амплитудно-частотной характеристике проявляется несколько локальных частотных экстремумов. Адекватность результатов численного моделирования обосновывается хорошим соответствием между формами амплитудно-частотных характеристик ускорения, рассчитанных с использованием аналитической модели динамического напряженно-деформированного состояния, и амплитудно-частотными характеристиками, зарегистрированными на поверхности покрытия нежестких дорожных одежд в натурных условиях. Рассмотрены вопросы построения чаш максимальных динамических прогибов на поверхности нежесткой дорожной одежды при воздействии кратковременного ударного нагружения. Приведены примеры применения разработанного комплексного подхода для решения задач оценки модулей упругости конструктивных элементов нежестких дорожных одежд, устроенных на участках автомобильных дорог. Определены перспективы развития предлагаемого подхода, неразрывно связанные с совершенствованием методов и средств анализа поведения дорожной конструкции в частотной области при динамическом нагружении.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2016;(2):111-130
О расслоении полосы по границе раздела упругих свойств Часть 2. Случай сдвиговой трещины
Аннотация
Рассматривается задача о полосе, образованной двумя изотропными упругими полосами меньшей и в общем случае различной толщины, обладающими различными упругими свойствами, с полубесконечной трещиной, проходящей по границе раздела упругих свойств и нагружаемой на бесконечности системой усилий. В первой части работы [1] представлена математическая формулировка, и с помощью преобразования Лапласа задача сведена к однородной задаче Римана-Гильберта с матричным коэффициентом. В предположении возможности пренебрежения влиянием нормальных напряжений на сдвиговые смещения и сдвиговых напряжений на нормальные смещения задача сведена к двум скалярным задачам Римана-Гильберта. Данная постановка может рассматриваться как приближенная для общего случая (данное приближение при этом заведомо не хуже общепринятого приближения, заключающегося в рассмотрении узкого слоя в рамках теории балок либо стержней) и как точная - для случая, когда прилегающие слои могут скользить друг относительно друга, но удерживаются силами адгезии от нормального отрыва. Путем факторизации получено точное аналитическое решение одной из сформулированных в [1] задач, а именно задачи о сдвиге. Получены асимптотические выражения для смещений берегов трещины вдали от ее вершины. Показано, что ведущие члены асимптотики смещений берегов трещины вдали от вершины соответствуют смещению стержня при граничных условиях типа упругой заделки, т.е. условиях пропорциональности смещения в точке заделки приложенному усилию. Для данного коэффициента пропорциональности получено аналитическое выражение. Также получены асимптотические выражения для поля напряжений вблизи вершины трещины (коэффициент интенсивности напряжений и скорость высвобождения энергии).
Вестник Пермского национального исследовательского политехнического университета. Механика. 2016;(2):131-142
Экспериментально-теоретическое исследование механического поведения 3D-композитов при квазистационарном разрушении
Аннотация
В конструкции всех современных самолетов и авиационных двигателей применяются композиционные материалы. Разработка таких материалов является очень трудоёмким, длительным и дорогостоящим процессом ввиду необходимости проведения большого количества экспериментов. Также существуют определенные ограничения, делающие невозможным определение некоторых свойств опытным путем. Для сокращения времени и стоимости разработки таких материалов применяются различные численные методики, такие как построение виртуального прототипа, проведение ряда численных экспериментов и применение различных разновидностей метода рассеянных повреждений. Метод механики рассеянных повреждений (МРП) - это популярный инструмент для моделирования развития повреждений в армированных композитах. Он основан на предположении, что материал с дефектами может быть заменён воображаемым материалом без дефектов, но с пониженным модулем упругости. Таким образом, наличие дефектов отражается только на упругих свойствах материала, и развитие повреждений регистрируется через потерю этих свойств. Для реализации одной из разновидностей данного метода на языке APDL был разработан макрос для пакета ANSYS. Рассматривается методика получения идеализированного виртуального прототипа трехмерно-армированного композиционного материала по результатам томографии. Представлена методика определения эффективных упругих свойств трехмерно-армированного композиционного материала путем проведения ряда численных экспериментов. Метод механики рассеянных повреждений применяется для нахождения параметров прочности. Показаны графики зависимости акустической эмиссии и кривые деформирования образцов при испытании на растяжение и сдвиг. Результаты моделирования сравниваются с полученными экспериментальными данными.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2016;(2):143-158
О геометрически нелинейных определяющих соотношениях упругого материала
Аннотация
Вестник Пермского национального исследовательского политехнического университета. Механика. 2016;(2):159-160