ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ СОБСТВЕННЫХ КОЛЕБАНИЙ СЛОИСТЫХ УСЕЧЕННЫХ КОНИЧЕСКИХ ОБОЛОЧЕК, ЗАПОЛНЕННЫХ ЖИДКОСТЬЮ
Бочкарев С.А., Лекомцев С.В.
Аннотация
Представлены результаты численных исследований собственных колебаний усечённых прямых слоистых конических оболочек вращения, полностью заполненных идеальной сжимаемой жидкостью. Поведение упругой конструкции и жидкой среды описывается в рамках классической теории оболочек, основанной на гипотезах Кирхгофа – Лява, и уравнений Эйлера. Уравнения движения оболочки совместно с соответствующими геометрическими и физическими соотношениями сводятся к системе обыкновенных дифференциальных уравнений относительно новых неизвестных. Акустическое волновое уравнение, записанное относительно гидродинамического давления, преобразуется к системе дифференциальных уравнений с помощью метода обобщённых дифференциальных квадратур. Решение сформулированной краевой задачи осуществляется методом ортогональной прогонки Годунова с численным интегрированием дифференциальных уравнений методом Рунге – Кутты четвёртого порядка точности. Для вычисления собственных частот колебаний используется сочетание пошаговой процедуры с последующим уточнением найденных значений в полученном диапазоне методом Мюллера. Достоверность получаемых результатов подтверждена сравнением с известными численными решениями. В работе детально проанализированы зависимости низших частот колебаний от угла конусности и угла армирования свободно опёртых, жёстко закреплённых и консольных двухслойных и трёхслойных композитных конических оболочек. Оценена возможность изменения низших частот и соответствующих им форм колебаний в зависимости от заданной комбинации угла конусности, условий закрепления краёв конструкции, схемы укладки и угла армирования композиционного материала. Проведённая обширная серия вычислений позволила выявить наличие конфигураций, для которых низшие частоты превышают значения, соответствующие эквивалентной слоистой круговой цилиндрической оболочке.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2025;(1):5-19
РЕШЕНИЕ ЗАДАЧИ О ПОПЕРЕЧНОМ ИЗГИБЕ ЭЛЕКТРОМАГНИТОУПРУГОЙ ПОЛУПЛОСКОСТИ С ОТВЕРСТИЯМИ И ТРЕЩИНАМИ
Калоеров С.А., Сероштанов А.В.
Аннотация
С использованием комплексных потенциалов теории электромагнитоупругого изгиба тонких плит решена задача об изгибе пьезоплиты в виде полуплоскости с отверстиями и трещинами. При этом функции, голоморфные вне контуров отверстий и трещин, разлагаются в ряды Лорана, а функции, голоморфные в нижних полуплоскостях, методом интегралов типа Коши выражаются через функции, сопряженные к указанным функциям. При таком подходе полученные суммарные функции точно удовлетворяют граничным условиям на прямолинейной границе полуплоскости, а для определения неизвестных коэффициентов рядов Лорана используются граничные условия на контурах отверстий и трещин, которые в работе удовлетворяются обобщенным методом наименьших квадратов, приводящим задачу к переопределенной системе линейных алгебраических уравнений, решаемой методом сингулярного разложения. Описаны результаты численных исследований электромагнитоупругого состояния полуплоскости с круговым отверстием или трещиной, с круговым отверстием и внутренней трещиной в перемычке, с круговым отверстием, имеющим краевую трещину в перемычке. Установлены закономерности изменения электромагнитоупругого состояния плиты в зависимости от ее материала и геометрических характеристик отверстий и трещин, их взаиморасположения. Установлено, что с приближением отверстия или трещины к прямолинейной границе значения моментов в точках перемычки резко возрастают, незначительно изменяясь в других зонах. Большая концентрация моментов наблюдается и в точках прямолинейной границы вблизи перемычки. Значения этих моментов особенно велики в задаче для полуплоскости с круговым отверстием, имеющим краевую трещину в перемычке. На значения изгибающих моментов значительно влияет учет пьезосвойств материала, особенно в зонах высокой концентрации изгибающих моментов, поэтому в этих случаях нельзя ограничиваться решением задачи теории упругости об изгибе плиты, а нужно решать задачу электромагнитоупругости.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2025;(1):20-33
Остаточные напряжения в зоне повреждения пластины из композиционного материала при статическом и динамическом воздействии стального индентора
Елеонский С.И., Писарев В.С., Чернов А.В.
Аннотация
Разработан и реализован новый подход, обеспечивающий определение главных компонент остаточных напряжений, которые возникают в результате, как статического, так и ударного контактного взаимодействия сферического индентора и плоской поверхности пластины, изготовленной из полимерного композиционного материала с продольно-поперечной укладкой. Экспериментальная методика включает сверление зондирующего отверстия и последующее измерение приращений диаметров этого отверстия в направлении главных остаточных деформаций методом спекл-интерферометрии. Картины интерференционных полос высокого качества, необходимые для реализации этой процедуры, визуализированы, как внутри контактной вмятины, так и в ближайшей внешней окрестности контура углубления, которое расположено на длинной оси симметрии прямоугольных образцов. Данные интерференционных измерений тангенциальных компонент перемещений используются для определения главных компонент остаточных напряжений на основе единственного решения корректно сформулированной обратной задачи. Это обстоятельство обеспечивает минимально возможную погрешность конечного результата. Ключевым моментом разработанного подхода является то, что остаточные напряжения определяются количественным образом по результатам прямых физических измерений приращения диаметров зондирующих отверстий, выполненных в различных зонах контактного взаимодействии. На этой основе установлено наличие значительных по величине главных компонент остаточных напряжений, которые возникают в зоне контактного взаимодействия стального индентора и поверхности композиционной пластины. Выявлено, что распределения главных компонент остаточных напряжений вдоль контура контактной вмятины по отношению к ее центру имеет радиально-симметричный характер. При этом тангенциальные компоненты являются растягивающими напряжениями, а радиальные напряжения являются сжимающими в обоих случаях. Наибольшие сжимающие компоненты в обоих случаях возникают в центре контактной вмятины. Представлено сравнение величин главных компонент остаточных напряжений, возникающих при статическом вдавливании сферического индентора и ударном воздействии, которые характеризуются практически одинаковым диаметром контактной вмятины. Полученные данные свидетельствуют, что остаточные напряжения можно использовать в качестве репрезентативного параметра, который способен установить корреляцию между снижением прочности композитной пластины вследствие ударного воздействия и количественными характеристиками поля остаточных напряжений в окрестности контактного углубления.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2025;(1):34-44
Применение экспериментально-расчетного подхода для идентификации деформационных и прочностных характеристик упругопластических материалов при квазистатическом и динамическом нагружениях
Баженов В.Г., Константинов А.Ю., Куканов С.С., Осетров Д.Л., Рябов А.А.
Аннотация
В работе представлены численные и экспериментальные исследования процессов деформирования и разрушения сталей 10Г2ФБЮ и Ст35 при различных видах квазистатического и динамического нагружения: растяжение сплошного стержня и болтов М8 с гладкой и резьбовой рабочей частью, а также сжатии компактных цилиндрических образцов. Подробно описаны и схематично показаны разработанные экспериментально-расчетные методики построения квазистатических и динамических истинных диаграмм деформирования в экспериментах на растяжение и сжатие. В результате применения этих методик определены квазистатические и динамические истинные диаграммы деформирования для сталей 10Г2ФБЮ и Ст35. Полученные квазистатические диаграммы деформирования стали 10Г2ФБЮ при растяжении и сжатии практически совпадают. Установлено, что при построении истинных диаграмм деформирования для стали 10Г2ФБЮ существенное влияние оказывает скорость деформаций при значениях выше 500 1/c. Для исследования и описания процессов деформирования и разрушения упругопластических материалов предложена и реализована связанная модель разрушения на основе кинетического уравнения накопления повреждений в сочетании с критерием прочности типа Писаренко-Лебедева. На основе экспериментально-расчетного подхода разработана новая методика по оснащению материальными параметрами для предложенной модели разрушения. Соответственно, определены материальные параметры связанной модели разрушения для стали Ст35. Результатами численного моделирования процесса разрушения болтов М8 с гладкой и резьбовой рабочей частью при квазистатическом растяжении показано количественное и качественное соответствие вида и характера разрушения с экспериментом. Совместное применение экспериментально-расчетного подхода и предлагаемой связанной модели разрушения позволяет более точно описывать процессы деформирования и разрушения упругопластических материалов.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2025;(1):45-58
Эффекты слагаемых высшего порядка малости в обобщенных рядах для представления полей, ассоциированных с вершиной трещины в анизотропных средах. Часть 1. Сравнение точного и асимптотического решений определения напряжений у вершины трещины
Степанова Л.В., Мушанкова К.А.
Аннотация
Статья посвящена анализу обобщенных на случай анизотропных линейно упругих сред асимптотических рядов, представляющих поля перемещений, деформаций и напряжений вблизи кончика острой трещины в анизотропных средах. Ключевой целью исследования является оценка механических полей в ближайшей окрестности вершины трещины в анизотропных материалах с простейшей кубической сингонией упругих свойств среды. Асимптотические ряды построены на основании классической теории упругости анизотропного тела. С помощью анализа рядов, ассоциированных с полями у вершины трещины, показано, что высшие приближения с коэффициентами, получившими название обобщенных коэффициентов интенсивности напряжений (ОКИН), оказывают значительное влияние на целостное представление поля напряжений при увеличении расстояния от кончика трещины. На примере плоской задачи для бесконечной анизотропной плоскости с кубической симметрией свойств (с различными тензорами упругих модулей, имеющими три независимых элемента) при различной ориентации трещины по отношению к осям симметрии упругих свойств показано, что в обобщенных рядах следует сохранять помимо первых двух слагаемых (содержащих коэффициенты интенсивности напряжений (КИН) и Т-напряжения) члены более высоких порядков малости. Построены угловые ( -распределения) зависимости компонент тензора напряжений в плоской задаче для двух различных расстояний от острия трещины, с принятием во внимание различное количество слагаемых ряда. Простое визуальное сопоставление и анализ зависимостей компонент тензора от полярного угла, выстроенных с учетом различного числа слагаемых ряда, ясно указывают на необходимость удержания высших приближений ряда. В целях расширения области, в которой справедливо решение в рядах, необходимо сохранять большее количество слагаемых. Все вычисления проведены для реальных материалов, упругие постоянные которых определены с помощью метода молекулярной динамики для монокристаллических веществ с гранецентрированной кубической (ГЦК) решеткой.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2025;(1):59-81
Метод Ритца при дискретной аппроксимации перемещений для расчета плит мостовых сооружений, подкрепленных ребрами различной формы
Афанасьева Е.О.
Аннотация
В статье рассматриваются плиты, подкрепленные ребрами жесткости, два противоположных конца которых закреплены жестко, а два другие свободны. Ребра расположены в одном направлении, параллельно свободным граням плиты. Предложена методика расчета таких конструкций, объединяющая понижение мерности задачи при помощи метода Л. В. Канторовича и метода Ритца при дискретной аппроксимации перемещений. Рассмотрены способы задания ребер жесткости при помощи единичных столбчатых функций. Выявлено, что единичные столбчатые функции являются удобным способом моделирования ребер, в том числе при необходимости учета ребер в краевых условиях. Рассмотрены ребра жесткости в виде короба, тавра, двутавра и сплошного прямоугольного сечения. Приведены их геометрические характеристики, используемые при построении функционала полной потенциальной энергии с учетом их дискретного расположения. Выполнено сравнение результатов расчетов, полученных предложенным методом с решением, полученным с помощью метода конструктивной анизотропии при решении краевой задачи для обыкновенных дифференциальных уравнений. В численном эксперименте установлено число конечных элементов по ширине ребра, необходимое и достаточное для решения задачи. Построены графики перемещений центральной части плит, на основе которых сделан вывод, что предложенный метод имеет достаточно хорошую сходимость с методом конструктивной анизотропии. Однако, плиты, посчитанные двумя этими способами, деформируются различно и метод Ритца при дискретной аппроксимации дает более точную картину деформирования: плита меньше деформируется в зоне ребер и больше между ними. Из результатов расчетов можно сделать вывод, что представленный метод является наиболее оптимальным с точки зрения трудоемкости и точности. Однако, при отношении суммарной ширины ребер к ширине плиты, равном 1:3 и более, можно использовать метод конструктивной анизотропии, как более простой.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2025;(1):82-91
НОВЫЙ ПОДХОД ДЛЯ ЧИСЛЕННОЙ ОЦЕНКИ РАЗРУШЕНИЯ КОМПОЗИТНОГО ЗАМКОВОГО СОЕДИНЕНИЯ «ЛАСТОЧКИН ХВОСТ» ПРИ РАСТЯГИВАЮЩЕЙ НАГРУЗКЕ
Гусейнов К.А., Кудрявцев О.А., Сапожников С.Б.
Аннотация
This paper proposes and tests a new approach to assessing stress-strain states and determining failure regions of thick-walled tapered composite elements. The mesostructure features of composites, such as ply drops and resin pockets, are not explicitly considered in this approach. The constitutive model based on the multiphase finite element approach was extended to describe the three-dimensional stress-strain state. The model consisted of orthotropic linear-elastic and isotropic elastic-plastic sub-elements which simulate the properties of the fibres and matrix, respectively. The assumption of independence of the shear curve from the type of stress state was adopted to describe the nonlinear deformation response in the model. The calibration of the interlaminar nonlinear response of the constitutive model was performed using the test results of V-notched specimens under combined loading. Then, the verified model was used to determine the delamination load of the dovetail specimens. The delamination load of the dovetail specimens was estimated with the Nouthwestern (NU-Daniel) and the Hashin failure criteria. Finite element analysis of the influence of interlaminar strength and taper angle on the failure load of the dovetail specimens was performed. Based on the results obtained, we proposed the method for determining the rational parameters of the dovetail specimen. It was shown that the new approach could be effective for strength assessment and nonlinear behaviour analysis of tapered thick-walled composite structures at interlaminar shear strains up to 6%.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2025;(1):92-103
МОДИФИЦИРОВАННАЯ ТЕОРИЯ НЕУПРУГОСТИ
Бондарь В.С., Абашев Д.Р., Шарова В.И.
Аннотация
Рассматриваются основные положения и уравнения модифицированной теории неупругости, относящейся к классу теорий течения при комбинированном упрочнении. Модифицированная теория неупругости является простейшим вариантом теории неупругости, который интегрирован в конечноэлементный комплекс для практических расчетов выработанного и остаточного ресурса материала конструкций в условиях повторности и длительности воздействия термомеханических нагрузок. Тензор деформации представляется в виде суммы тензоров упругой и неупругой деформаций, т.е. здесь нет условного разделения необратимой (неупругой) деформации на деформации пластичности и ползучести. Упругая деформация следует закону Гука, обобщенному на неизотермическое нагружение. В пространстве составляющих тензора напряжений вводится поверхность нагружения, которая изотропно расширяется или сужается и смещается в процессе нагружения. Для радиуса поверхности нагружения (изотропное упрочнение) формулируется эволюционное уравнение, обобщенное на неизотермическое нагружение и процессы возврата механических свойств при отжиге. Смещение поверхности нагружения (анизотропное упрочнение) описывается на основе эволюционного уравнения с трехчленной структурой, обобщенного на неизотермическое нагружение и процессы снятия микронапряжений (смещения) при отжиге. Для разделения процессов монотонного и циклического деформирования в пространстве тензора неупругой деформации вводится поверхность памяти, ограничивающая область циклического деформирования. Для описания процессов посадки и вышагивания (ratcheting) петли неупругого деформирования при несимметричных циклических нагружениях вводится модификация теории неупругости. Модификация теории неупругости сводится к тому, что при формулировке эволюционного уравнения для микронапряжений определяющий (материальный) параметр уравнения для микронапряжений первого типа принимается зависящим от накопленной неупругой деформации на основе разных соотношений как при циклическом деформировании, так и монотонном. Для определения тензора скоростей неупругой деформации используется ассоциированный (градиентальный) закон течения. Сформулированы условия упругого и неупругого состояний. Для описания нелинейных процессов накопления повреждений вводится кинетическое уравнение накопления повреждений, базирующееся на работе микронапряжений второго типа на поле неупругих деформаций. Кинетическое уравнение обобщено на неизотермическое нагружение и процессы охрупчивания и залечивания. Выделяются материальные параметры и функции, замыкающие теорию, формулируется базовый эксперимент и метод их определения. Приводятся материальные параметры и функции бронзового сплава БрХ08-Ш при температурах 20, 400, 500, 600 °С. Проводится верификация теории при циклическом изотермическом деформировании и разрушении (малоцикловая прочность) в условиях высоких температур. Рассматривается также ползучесть и длительная прочность при изотермических и неизотермических нагружениях. Результаты расчетов сопоставляются с результатами экспериментов.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2025;(1):104-116
О ВЛИЯНИИ МЕХАНИЗМА УПРОЧНЕНИЯ НА ПРИСПОСОБЛЯЕМОСТЬ СОСУДА ПОД ДАВЛЕНИЕМ В УСЛОВИЯХ ТЕРМОЦИКЛИЧЕСКОГО НАГРУЖЕНИЯ
Федоренко Р.В., Лукин А.В., Муртазин И.Р.
Аннотация
В классических трудах по исследованию приспособляемости сосуда под давлением в условиях термоциклического нагружения материал конструкции принимается упруго-пластическим без упрочнения (идеальная пластичность). Решение задачи в данной постановке может быть получено аналитически и, как правило, приводится в виде так называемой «диаграммы Бри», построенной в осях механических и тепловых напряжений. В то же время поведение реального материала в условиях циклического термо-силового нагружения требует применения математических моделей изотропного, кинематического или смешанного изотропно-кинематического упрочнения, что оказывает значительное влияние на приспособляемость конструкции. Современные отечественные и зарубежные нормы проектирования (ПНАЭ, ГОСТ 59115, ASME, RCC-MR) высокотемпературных реакторных установок атомных станций (реакторы с жидкометаллическим теплоносителем, высокотемпературные газоохлаждаемые реакторы и др.) допускают возникновение пластических деформаций в материале, но ограничивают уровень их накопления в течение всего жизненного цикла конструкции. В частности, действующие стандарты регламентируют использование как упрощенной методики, основанной на классическом решении Бри для идеально-пластического материала, так и прямого конечно-элементного расчета жизненного цикла конструкции в условиях пластического деформирования и высокотемпературной ползучести. В настоящей работе рассматривается задача о термоциклическом нагружении сосуда под давлением с применением различных моделей упрочнения. Разработан численный алгоритм построения и анализа эволюции диаграмм приспособляемости в зависимости от числа циклов нагружения. Проводится конкретный числовой расчет процесса циклического термосилового нагружения сосуда под давлением, изготовленного из стали с известными параметрами моделей упруго-пластического деформирования. в контексте прогнозирования ресурса конструкции на основе действующих норм и правил проектирования корпусов и оборудования реакторных установок.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2025;(1):117-128
Сравнение энергетических спектров сигналов акустической эмиссии по фрактальным размерностям и сравнительным диаграммам
Волков А.Е., Черняева Е.В., Казаринов Н.А., Волкова Н.А.
Аннотация
Рассмотрены проблемы анализа сигналов акустической эмиссии (АЭ), направленного на исследование изменения состояния материала образцов в результате их испытаний. Ввиду того, что получить информацию непосредственно о характеристиках источников АЭ очень сложно ввиду искажений сигналов из-за наличия дисперсии, неодинакового затухания на разных частотах, отражений от свободных поверхностей образца, искажений, создаваемых датчиком, волноводом и усилителем электрических колебаний, в данной работе предложено сравнивать характеристики сигналов, полученных на "свежих" (контрольных) образцах и образцах, испытавших механические воздействия. Рассмотрено выявление различий сигналов путем построения сравнительной диаграммы спектров, представляющей модуль отношения усредненных Фурье-образов сигналов. Другой способ анализа сигналов основан на поиске таких их характеристик, которые не изменяются под действием многих из перечисленных искажений. Поскольку график спектральной плотности мощности имеет сложную изломанную форму, он может быть рассмотрен, как фрактальная кривая. Ее важнейшей характеристикой является фрактальная размерность, которая определяется условиями формирования сигналов и, следовательно, может служить признаком для их классификации. Для примера исследовали образец из стали 20, который подвергли циклическому нагружению от напряжения σ = 0 до напряжения σmax =1,2σт, (σт - предел текучести) и разгрузке до нуля, с частотой f = 20 Гц. Испытание прекратили на 8851-м цикле, когда сформировалась «шейка», а относительное удлинение составило 15%. Сравнивали сигналы АЭ, возникающие при индентировании образца в областях, расположенных на различных расстояниях от шейки. Фрактальная размерность спектров мощности при приближении к зоне шейки снижалась с 0.72 до 0.62. На сравнительной диаграмме спектров выделялись пики вблизи частот 270 кГц и 680 кГц. Таким образом, рассмотренные способы сравнения сигналов позволяют оценить степень изменения состояния образцов в результате механических испытаний.
Вестник Пермского национального исследовательского политехнического университета. Механика. 2025;(1):129-138