ПРОЧНОСТЬ И ОПТИМИЗАЦИЯ ЭЛЕМЕНТОВ КОНСТРУКЦИЙ: СОВРЕМЕННЫЕ ПОДХОДЫ И ВЫЗОВЫ

DURABILITY AND OPTIMIZATION OF STRUCTURAL ELEMENTS: MODERN APPROACHES AND CHALLENGES

Popov Oleg Alekseevich,

Student, Kuban State Agrarian University named after I.T. Trubilin,

Russia, Krasnodar

Nikolenko Alexander Yuryevich,

scientific supervisor. Assistant, Kuban State Agrarian University named after I.T. Trubilin,

Russia, Krasnodar

АННОТАЦИЯ

В статье рассматриваются фундаментальные аспекты обеспечения прочности элементов конструкций в контексте задач топологической и параметрической оптимизации, причём особое внимание уделяется противоречию между снижением материалоёмкости и сохранением несущей способности. На основе анализа классических критериев прочности (Губера–Мизеса, Мора–Кулона) и численных методов (МКЭ, генетические алгоритмы) показано, что внедрение оптимизационных процедур требует пересмотра традиционных коэффициентов запаса. Актуальность темы обоснована необходимостью проектирования лёгких и одновременно надёжных конструкций в авиастроении, ракетно-космической технике и энергетическом машиностроении.

ABSTRACT

The article examines the fundamental aspects of ensuring the strength of structural elements in the context of topological and parametric optimization problems, with special attention being paid to the contradiction between reducing material consumption and maintaining load-bearing capacity. Based on the analysis of classical strength criteria (Huber–Mises, Mohr–Coulomb) and numerical methods (FEM, genetic algorithms), it is shown that the introduction of optimization procedures requires a revision of traditional stock coefficients. The relevance of the topic is justified by the need to design lightweight and at the same time reliable structures in the aircraft industry, rocket and space technology and energy engineering.

Ключевые слова: прочность, оптимизация топологии, метод конечных элементов, концентрация напряжений, критерии разрушения, многоцелевая оптимизация, аддитивные технологии.

Keywords: strength, topology optimization, finite element method, stress concentration, fracture criteria, multipurpose optimization, additive technologies.

Современный этап развития технических наук характеризуется всё более жёсткими требованиями к массогабаритным показателям несущих систем, причём снижение массы конструкции, как правило, вступает в прямое противоречие с требованием её прочности и долговечности. Поскольку в таких отраслях, как аэрокосмическая техника, автомобилестроение и мостостроение, каждый килограмм материала непосредственно влияет на эксплуатационные расходы и энергоэффективность, проблема поиска компромисса между минимальной материалоёмкостью и максимальной несущей способностью приобретает ключевое значение. Актуальность данной темы дополнительно обостряется переходом к аддитивным технологиям, которые позволяют реализовывать сложные топологические формы, недостижимые при традиционной обработке резанием; однако такие формы требуют новых критериев оценки прочности, поскольку классические методы расчёта стержневых и тонкостенных систем не всегда адекватно описывают поведение конструкций с развитой внутренней архитектурой.

Несмотря на то, что основы механики деформируемого твёрдого тела были заложены ещё в XIX веке трудами Коши, Сен-Венана и Кирхгофа, систематическое применение оптимизационных алгоритмов к задачам распределения материала в упругих телах началось лишь во второй половине XX века благодаря работам Бендсё и Сигмунда. Тем не менее, и в настоящее время сохраняется разрыв между теоретически оптимальными решениями (например, формами, полученными методом асимптотической гомогенизации) и реальными конструкциями, где необходимо учитывать технологические ограничения, усталостную прочность и чувствительность к начальным несовершенствам. Следовательно, разработка методов оптимизации, которые не только снижают массу, но и гарантируют отсутствие хрупкого или вязкого разрушения при всех регламентных нагрузках, представляет собой одну из центральных задач вычислительной механики.

Для корректной постановки задачи оптимизации необходимо выбрать функцию отклика, которая количественно характеризует достижение предельного состояния. В теории упругости наибольшее распространение получил критерий энергии формоизменения (фон Мизеса), согласно которому пластическое течение начинается в тот момент, когда интенсивность касательных напряжений достигает критического значения, определяемого пределом текучести при одноосном растяжении. Данный критерий хорошо зарекомендовал себя для пластичных материалов, таких как малоуглеродистые стали и алюминиевые сплавы, но его использование в задачах топологической оптимизации сталкивается с проблемой недифференцируемости функции напряжений в областях с нулевой плотностью материала.

Альтернативные подходы, включая критерий Мора–Кулона для сыпучих сред и критерий максимального главного напряжения для хрупких материалов (например, высокопрочных чугунов или керамики), требуют модификации при интегрировании в градиентные методы оптимизации, поскольку поверхности текучести в пространстве напряжений имеют угловые точки, что порождает неединственность направления спуска. В этой связи, как показывают исследования последних лет, более продуктивным оказывается использование релаксированных критериев, в которых предельная поверхность аппроксимируется гладкой функцией, дифференцируемой во всей области допустимых состояний.

Наиболее эффективным инструментом решения сформулированной выше задачи признаётся метод конечных элементов (МКЭ) в сочетании с методами оптимальности (Optimality Criteria, OC) или методами подвижных асимптот (Method of Moving Asymptotes, MMA). При этом каждый шаг оптимизации требует вычисления поля напряжений во всей сетке конечных элементов, что сопряжено с высокой вычислительной стоимостью — особенно для трёхмерных задач с числом степеней свободы порядка 10⁶–10⁷.

Ключевой проблемой является зависимость напряжений от плотности материала в промежуточных элементах (серых зонах), где 0 < ρ < 1. Согласно традиционной модели SIMP (Solid Isotropic Material with Penalization), жёсткость интерполируется как E(ρ) = E_min + ρ^p (E₀ – E_min), а напряжения обычно вычисляются через деформации по закону Гука: σ(ρ) = E(ρ)·ε. Однако такая интерполяция приводит к занижению напряжений в элементах с промежуточной плотностью, что побуждает оптимизатор использовать эти зоны для обхода ограничений по прочности. Для устранения этого артефакта исследователями предложены модифицированные модели, в частности модель RAMP (Rational Approximation of Material Properties), а также релаксация ограничений с использованием p-нормы напряжений.

Важно подчеркнуть, что даже при использовании корректных интерполяционных схем сходимость итерационного процесса остаётся медленной из-за невыпуклости целевой функции и наличия локальных минимумов. В этой связи многие авторы обращаются к эвристическим методам (генетические алгоритмы, роевой интеллект), которые, хотя и требуют на порядки больше вычислений, позволяют выходить из зон локального оптимума, характерных для градиентных методов.

Отдельным, исключительно сложным направлением является оптимизация элементов конструкций, работающих в условиях циклического нагружения, поскольку критерии статической прочности здесь оказываются недостаточными. Усталостное разрушение, инициируемое в зонах локальной концентрации напряжений (галтелях, отверстиях, сварных швах), требует введения ограничений на амплитуду циклов или на размах напряжений. В рамках оптимизационных постановок это обычно реализуется через метод эквивалентных напряжений по Содербергу или Гудману, где учитывается влияние среднего напряжения.

Однако фундаментальное противоречие заключается в том, что оптимальная по жёсткости топология часто создаёт области с высокой градиентностью поля напряжений, что негативно сказывается на усталостной долговечности. Например, классическая MBB-балка, оптимизированная для минимизации податливости, демонстрирует ярко выраженные концентраторы напряжений в узлах решётки, которые становятся очагами зарождения усталостных трещин. Следовательно, для обеспечения ресурса требуется либо введение дополнительных ограничений на максимальный градиент напряжений, либо переход к многокритериальной оптимизации по Парето, где наряду с жёсткостью минимизируется пиковое эквивалентное напряжение.

Даже при получении теоретически оптимального распределения материала, удовлетворяющего всем прочностным ограничениям, инженер сталкивается с проблемой его технологической реализуемости. Так, топологии, содержащие замкнутые полости, внутренние консоли и тонкие перемычки, не могут быть изготовлены методами литья или механической обработки без использования специальной оснастки, что резко удорожает производство. Ситуация кардинально меняется с внедрением аддитивных технологий (селективное лазерное сплавление, электронно-лучевая плавка), которые позволяют воспроизводить структуры с внутренней ячеистостью — так называемые топологически оптимизированные «скелеты».

При этом возникает новая научная проблема: прочностные свойства материала, полученного методом послойного синтеза, анизотропны и зависят от ориентации детали при печати, что должно быть учтено в оптимизационной модели. Экспериментальные данные свидетельствуют, что модуль упругости и предел текучести образцов, напечатанных вертикально, могут на 15–25% отличаться от горизонтально ориентированных образцов. Следовательно, задача оптимизации прочности должна быть дополнена ограничениями, связанными с технологическими углами наклона элементов, что превращает её в многослойную иерархическую задачу, решаемую методами параллельного проектирования.

Проведённый анализ позволяет утверждать, что современные методы оптимизации элементов конструкций по критериям прочности достигли уровня, при котором возможно проектирование деталей с массой, близкой к теоретическому минимуму, однако при условии учёта ряда дополнительных факторов. К числу наиболее значимых направлений, требующих дальнейшей разработки, следует отнести: создание эффективных алгоритмов сглаживания поля напряжений в задачах топологической оптимизации; разработку гибридных методов, сочетающих градиентный спуск и эволюционные стратегии для преодоления локальных оптимумов; а также построение верифицированных моделей усталостной прочности для аддитивно производимых структур с развитой пористостью.

Кроме того, остаётся открытым вопрос о нормировании коэффициентов запаса для оптимальных конструкций: поскольку традиционные нормативные коэффициенты (например, 1,5 для пластичных материалов) были эмпирически установлены для конструкций с «избыточным» материалом, их прямое перенесение на тонкие оптимизированные топологии может оказаться либо излишне консервативным, либо, напротив, недостаточным. В связи с этим представляется перспективным вероятностный подход, в котором запас прочности связывается с распределением дефектов материала, выявляемых методами неразрушающего контроля на этапе изготовления.

Таким образом, синтез идей механики разрушения, вычислительной топологической оптимизации и цифрового производства образует тот междисциплинарный фронт, где в ближайшие десятилетия будет формироваться новая парадигма проектирования несущих систем — систем, которые одновременно легки, прочны и технологичны.

Список литературы:

1. Бендсё М.П., Сигмунд О. Топологическая оптимизация. Теория, численные методы и приложения. — М.: Мир, 2008. — 384 с.
2. Леонтьев Н.В., Соколова А.И. Учёт усталостной прочности при топологической оптимизации авиационных кронштейнов // Вестник МГТУ им. Н.Э. Баумана. Сер. «Машиностроение». — 2021. — №4. — С. 67–81.
3. Rozvany G.I.N. A critical review of established methods of structural topology optimization // Structural and Multidisciplinary Optimization. — 2009. — Vol. 37. — P. 217–237.
4. Зеленцов В.Б., Миронов А.А. О сходимости методов подвижных асимптот в задачах с ограничениями по напряжениям // Вычислительная механика сплошных сред. — 2022. — Т.15. — №2. — С. 144–159.