ТРЕЩИНООБРАЗОВАНИЕ В ЛИСТОВОМ ЭЛЕМЕНТЕ КОНСТРУКЦИИ ПРИ НЕОДНОРОДНОМ НАПРЯЖЕННОМ ПОЛЕ

Полный текст

Введение Процесс разрушения конструкционных материалов зависит от особенностей структуры материала, его химического состава, вида напряжения и других. В настоящее время известны различные механизмы зарождения трещины [1-4]. Изучение вопросов разрушения элементов конструкций имеет важное практическое значение. После фундаментальных работ В.М. Мирсалимова [3, 5] проблема зарождения трещины интенсивно исследуется в Азербайджане [6-38]. Постановка задачи Рассмотрим однородную изотропную среду. На бесконечности действуют напряжения, являющиеся полиноминальными функциями декартовых координат x и y. По мере нагружения среды (листового элемента конструкции) в процессе работы силовой нагрузкой в материале будут возникать зоны предразрушения, которые моделируются как области ослабленных межчастичных связей материала. Взаимодействие берегов этих зон моделируется [3] путем введения между берегами зоны предразрушения связей в состоянии пластического течения. Физическая природа таких связей и размеры зон пластического течения зависят от вида материала среды. Поскольку указанные зоны (прослойки) малы по сравнению с остальной частью листового элемента, их можно мысленно заменить разрезами, поверхности которых взаимодействуют по закону, соответтвующему действию удаленного материала. Возникновение зародышевой трещины рассматривается как процесс пере- хода области предразрушения в область разорванных связей между поверхностями материала листового элемента. Исследования [1, 2, 4, 39] возникновения областей с нарушенной структурой материала, показывают, что изначально зоны предразрушения представляют собой узкий вытянутый слой, затем с ростом нагрузки внезапно появляется вторичная система зон, содержащих материал с частично нарушенными связями. Пусть в среде (листовом элементе конструкции) имеется N прямолинейных зон предразрушения (рис. 1). Рассмотрим локальные систем координат , оси которых совпадают с зоной предразрушения и образуют углы с осью х. Начала этих систем координат будут располагаться в центрах зон предразрушения. Длина (k = 1,2,…,N) зон пластического течения заранее неиз- вестна и подлежит определению в процессе решения задачи. Рис. 1. Расчетная схема задачи о трещинообразовании в листовом элементе При действии внешних нагрузок на среду (листовой элемент конструкции) в связях между берегами зон предразрушения, будут возникать нормальные ?s и касательные ?s усилия. Следовательно, к берегам зон предразрушения будут приложены нормальные ?s и касательные ?s напряжения. Граничные условия на берегах зон пластического течения имеют вид: , на , . (1) Для определения внешней нагрузки, при которой происходит трещинообразование, постановку задачи дополним критерием критического раскрытия берегов зоны пластического течения: , (2) где - характеристика сопротивления трещинообразованию материала листового элемента конструкции. Условие (2) является условием появления трещины, т.е. разрыва межчастичных связей материала среды. Метод решения краевой задачи С помощью принципа суперпозиции напряженное состояние в среде с произвольной системой зон предразрушения в состоянии пластического течения при постоянном напряжении представим в вилле суммы двух напряженных со- стояний , , , (3) где , , - компоненты тензора напряжений в листовом элементе конструкции без зон пластического течения, когда на бесконечности действуют напряжения, являющиеся полиноминальными функциями декартовых коорди- нат x и y; , , - компоненты тензора напряжений для среды с произвольной системой зон пластического течения с исчезающими на бесконечности напряжениями. Для компонент , , тензора напряжений имеем: , , (4) , , (5) . Отметим, что функции (5) в зависимости от значений коэффициентов Aj и Bj (j = 0,1,2,…,m) определяют напряженное состояние в листовом элементе без зон пластического течения. С учетом формул (3) краевое условие (1) запишем в следующем виде , на , . (6) Напряжения , , и перемещения , выразим через две кусочно-аналитические функции и . Тогда краевые условия задачи (1) с помощью формул Колосова-Мусхелишвили [40] можно записать в виде граничной задачи для отыскания комплексных потенциалов и : (k = 1,2,…,N), (7) где tk - аффикс точек берегов k-ой зоны пластического течения. Комплексные потенциалы и , дающие решение краевой задачи (7) ищем в виде , (8) . Здесь ; ; (k = 1,2,…,N) - искомые функции, характеризующие раскрытие берегов зон пластического течения: . (9) Определив по формуле напряжения на оси (k = 1,2,…,N) и подставив их в краевые условия (7), по- лучим систему N сингулярных интегральных уравнений в виде (10) (n = 1,2,…,N), , , , . Ядра полученной системы интегральных уравнений (10) регулярны за исключением случая, когда n = k, в этом случае переходит в сингулярное ядро Коши. Систему (10) можно записать в следующем виде: , (n = 1,2,…,N). (11) К системе сингулярных интегральных уравнений (11) для внутренних зон пластического течения следует добавить дополнительные равенства: (k = 1,2,…,N). (12) Эти равенства выражают условие однозначности смещений при обходе контура зоны пластического течения. Если представить неизвестные функции и функцию нагрузки в виде ; , (13) то после отделения действительных и мнимых частей из системы N комплексных сингулярных интегральных уравнений получаем для нахождения и систему 2N действительных сингулярных интегральных уравнений. Чтобы провести алгебраизацию сингулярных интегральных уравнений [39, 40, 41], сначала приведем в системе (11) и дополнительных условиях (12) все интервалы интегрирования к одному интервалу . Сделав затем замену переменных , ( , ), систему уравнений (11) и условия (12) запишем в виде (n = 1,2,…,N), (14) . Для неизвестных функций и свободных членов сохранены прежние обозначения. Решение системы (14) представим в виде [39, 40, 41]: , (15) где - новая неизвестная функция, регулярная в интервале . Использование квадратурных формул типа Гаусса-Чебышева для сингулярного интеграла позволяет свести систему сингулярных интегральных уравнений (14) при дополнительных условиях (12) к конечной системе алгебраических уравнений для определения неизвестных : (16) ( r = 1,2,…,M-1), (n = 1,2,…,N). Здесь значения и определяются формулами: (m = 1,2,…,M), (r = 1, 2,…,M - 1). Переходя в (16) к сопряженным значениям, получаем еще алгебраических уравнений. Для левой части соотношения (9) имеем (k = 1,2,…,N). (17) Полученные системы оказались связанными и должны решаться совместно. Для их замкнутости не хватает комплексных уравнений, определяющих местоположение и размеры зон пластического течения. Так как напряжения в листовом элементе конструкции ограничены, решение сингулярных интегральных уравнений должно искаться в классе всюду ограниченных функций. Такое решение существует при выполнении условий разрешимости сингулярных интегральных уравнений. Поэтому полученные алгебраические системы не являются пока замкнутыми. Записывая условия конечности напряжений у вер- шин зон пластического течения (условия разрешимости), находим еще комплексных уравнений (n = 1,2,…,N), (18) . Для определения предельного состояния, при котором возникает трещина, используем критическое условие (2). Тогда условием, определяющим предельное значение внешней нагрузки, будет равенство . (19) Модуль вектора смещения на берегах зон пластического течения при удобно представить в виде , (k = 1,2,…,N), (20) ; . где - число узловых точек, содержащихся в интервале . Полученная алгебраическая система (19), (25)-(26) из-за неизвестных размеров зон пластического течения является нелинейной. Совместное решение полученных алгебраических систем позволяет найти напряженно- деформированное состояние листового элемента, критическую внешнюю нагрузку, местоположение и размеры зон пластического течения для предельно- равновесного состояния листового элемента конструкции, когда возникает трещина. Объединенная алгебраическая система уравнений из-за неизвестных величин (k = 1,2,…,N) оказывается нелинейной. Для ее решения использовали метод последовательных приближений [39], состоящий в следующем. Решаем алгебраическую систему (19), (25)-(26) при некоторых определенных значениях (k = 1,2,…,N) относительно остальных неизвестных, которые входят в объединенную систему линейным образом. Поскольку значения и соответствующие им значения остальных неизвестных не будут, вообще говоря, удовлетворять уравнениям (27), подбирая значения параметров , многократно повторяем вычисления. Когда уравнения (27) будут удовлетворены с заданной точностью, вычисления прекращаются. В каждом приближении алгебраическая система решалась методом Гаусса с выбором главного элемента. На рис. 2 представлены графики зависимость длины зон пластического течения от безразмерной нагрузки при чистом изгибе. В расчетах было принято M = 30; ? = 0,3; ; . На рис. 3 приведена зависимость безразмерной предельной нагрузки от относительной длины зоны предразрушения . Здесь R характерный размер листового элемента. Рис. 2. Зависимость длины зон пластического течения от безразмерной нагрузки Выводы Разработан эффективный способ решения задач о зарождении трещин в металлическом листовом элементе конструкции при действии внешнего неоднородного напряженного поля. На основе разработанной расчетной модели исследовано трещинообразование в металлическом листовом элементе конструкции с системой произвольно размещенных зон предразрушения при различных силовых нагрузках. Рис. 3. Зависимость предельной нагрузки от относительной длины оны пластического течения

Об авторах

ШАХИН ОГЛЫ ГАСАНОВ

Азербайджанский технический университет

Автор, ответственный за переписку.
Email: iske@mail.ru

д-р техн. наук, профессор

Азербайджан, AZ1073, Баку, пр. Г. Джавида, 25

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