Full Text

Введение Тема расчетной оценки прочности элементов конструкций (в том числе сварных) из малопластичных материалов с острыми концентраторами напряжений сохраняет свою актуальность уже на протяжении многих десятилетий [1-7], поскольку появляются новые высокопрочные материалы, совершенствуются технологии производства и методы расчета, учитывающие все особенности геометрии концентраторов. Для сварных соединений острые концентраторы (непровары, грат, резкие изменения толщины и др.) являются опасными дефектами, учет которых при оценке несущей способности соединений обязателен [8-11]. Традиционные подходы конструкционной прочности основаны на использовании локального условия прочности в виде неравенства max σ < σB, (1) в котором в левой части стоит максимальное действующее эквивалентное напряжение (или соответствующая деформация для деформационного критерия) в опасной точке, а в правой - предел прочности (или разрушающая деформация, определяемая на стандартных образцах), и применяются в конструкциях, имеющих концентраторы известного и достаточно большого радиуса. Определение «достаточно большой» оценивает многократное превышение величины радиуса над размером элемента структуры материала, определяющего его физическую неоднородность: зерна в металлах, гравий в бетоне, волокна в композитах, глобулы или пачки в полимерах и т.п. Это является следствием допущения об однородности материала в классической механике деформируемого твердого тела [5, 11]. В случае острых концентраторов (с неизвестными малыми или даже нулевыми радиусами в вершине) применение подходов конструкционной прочности невозможно, и здесь включается в работу механика разрушения, оперирующая, в частности, для малопластичных (хрупких) материалов понятием коэффициента интенсивности напряжений в вершине острого концентратора: max K1 < K1С, (2) где в левой части - максимальное значение коэффициента интенсивности напряжений (в вершине) при отрыве, а в правой - его предельное значение, определенное из опытов на образцах с трещинами известной длины. Напряженное состояние вблизи вершины трещины (рис. 1, r ® 0 на продолжении линии трещины) характеризуется при этом сингулярностью со степенью m = 0,5, формула (3): (3) Рис. 1. Напряженное состояние у вершины трещины Fig. 1. The stress state at the top of the crack где σx(r) - напряжение в направлении оси x в точке, удаленной от вершины трещины на расстояние r; K1 - коэффициент интенсивности напряжений, который для короткой трещины длиной L при одноосном растяжении напряжениями σ определяется как K1 = σ(πL)0,5. (4) Для трехточечного изгиба образцов с трещиной нормального отрыва (двугранный угол v = 0) есть численное решение Броека [12] для коэффициента интенсивности напряжений K1 в зависимости от нагрузки Р, размеров поперечного сечения образца (b, h), расстояния между опорами (l) и длины трещины L (рис. 2): (5) Однако для геометрии, в которой имеется ненулевой двугранный угол (пример - образец Шарпи по ГОСТ Р ИСО 148-1-2013, ASTM E23 и др.), в литературе нет выражений, описывающих зависимость коэффициента интенсивности напряжений от величины двугранного угла v и степени сингулярности m. Имеются данные [14, 15], что степень сингулярности m снижается с увеличением двугранного угла v ® π вплоть до нуля по зависимости (6) а коэффициент интенсивности определяется выражением (7) В формуле (7) функция πx(r) заранее неизвестна и может быть найдена, например, из решения соответствующей задачи теории упругости методом конечных элементов [20, 21] для вырезов с заданным двугранным углом v. При этом коэффициент интенсивности напряжений имеет различную размерность для различных двугранных углов v (МПа×мm). Классические подходы линейной механики разрушения к оценке прочности тел с острыми надрезами необходимо дополнить работами, в которых рассматривалась удельная энергия деформирования в заданном объеме у вершины надреза, которая сравнивалась с предельной энергией, полученной из испытаний стандартных образцов до разрушения [13, 17-19]. Ключевым словом здесь является «заданный объем», иногда это «заданная площадь» или «заданное расстояние». Иначе говоря, вводится размерный фактор, который является параметром в критерии разрушения [16, 18, 19]. Работы по анализу полей напряжений в окрестности острых надрезов публиковались рядом авторов [23-28] на протяжении последних 30 лет, в них разрабатывались аналитические и численные подходы, рассматривались упругие и малопластичные материалы с позиции линейной механики разрушения с использованием сингулярных представлений о распределении компонент напряжений в окрестности вершины выреза. К отмеченным выше можно добавить работы, в которых удерживаются дополнительно несингулярные члены, уточняющие поля напряжений у вершины надреза (расширенные ряды М.Уильямса) и применяемые в анализе коэффициентов интенсивности напряжений, оценке устойчивости направления движения трещины, усталостной прочности деталей с надрезом [29-34]. Наличие несингулярных членов иллюстрируется конечно-элементными расчетами, в которых наличие этих членов проявляется на удалении от вершины острого выреза и может быть существенным при оценке прочности деталей из материалов с крупнозернистой (крупноячеистой) структурой (бетон, композиты, пеноматериалы и др.). Влияние на распределение напряжений на удалении от вершины трещины высших приближений в асимптотическом разложении М.Уильямса подробно рассмотрено в обзоре [39], где, в частности, точно решена задача о двух коллинеарных трещинах конечной длины в пластине при двух разных случаях нагружения. Однако результаты такого рода работ пока не доведены до включения в методики оценки хрупкой прочности. Важно отметить, что линейная механика разрушения для прогноза прочности требует проведения экспериментальных исследований прочности образцов с острыми вырезами с конкретным двугранным углом v, вычисления критического коэффициента интенсивности и сравнения текущих значений с критическим значением , что влечет за собой существенные временные и материальные затраты. Учет высших приближений в асимптотическом разложении дополнительно усложняет эту проблему. Иной подход предложен в работе [35], где полагается, что в окрестности вершины v-образного надреза возникает зона рассеянных микроповреждений, снижающих модуль упругости по приближении к вершине выреза из-за почти равноосного растяжения в этой зоне, а по достижении определенной удельной энергии деформации в этой зоне возникает макротрещина, которая и приводит к полному разрушению. На аморфном графите показано хорошее соответствие расчетных и экспериментальных нагрузок разрушения. К недостаткам следует отнести необходимость определения четырех параметров, определяющих нагрузку разрушения. В недавней статье [36] отмечается, что двух- (и более) параметрические критерии прочности имеют шанс успешного соответствия экспериментальным данным. И, как правило, одним из них является размерный параметр, например размер области интенсивного высвобождения энергии, представительный размер конечного элемента, размер зоны микрорастрескивания и т.п. В связи с этим в данной работе для оценки прочности элементов конструкции с острым v-образным вырезом предлагается рассмотреть двухпараметрический нелокальный силовой критерий прочности (8) в котором в зоне вершины выреза производится вычисление максимальных осредненных на заданной базе А напряжений и сравнение их с предельным значением SIG. Идейной базой этого критерия являются подходы Новожилова, Нейбера и Бажанта [5, 37, 38]. Преимущество данного критерия в том, что параметры А и SIG, определенные на образцах с двумя-тремя двугранными углами, могут использоваться для всех других конфигураций образцов из того же материала. В данной статье при трехточечном изгибе образцов с острыми v-образными вырезами (двугранные углы 0-180°) различной глубины с помощью метода конечных элементов (пакет ANSYS Workbench [20]) рассмотрены особенности напряженно-деформированного состояния, определены коэффициенты интенсивности напряжений и показатели степени сингулярности, проведены расчеты нагрузок разрушения с использованием нелокального критерия прочности (8), проведено их сопоставление с экспериментальными данными для квазистатического трехточечного изгиба образцов с надрезами из органического стекла (ПММА). 1. Экспериментальное исследование нагрузок разрушения при изгибе образцов с острыми надрезами Из листа органического стекла типа ТОСП (ГОСТ 17622-72) толщиной b = 3,9 мм были вырезаны образцы высотой h = 10 мм и длиной 75 мм. Специальными фрезами с двугранными углами 25, 90 и 120° и малым радиусом в вершине (менее 10 мкм) с водяным охлаждением без сколов были выполнены вырезы глубиной L = 0,25 … 5 мм. На испытательной машине INSTRON 5882 с силоизмерителем 5кН были проведены испытания на трехточечный изгиб с расстоянием между опорами l = 60 мм и радиусом опор 5 мм (рис. 2). Рис. 2. Схемы нагружения образца на трехточечный изгиб (a) и КЭ-модели (б) Fig. 2. Scheme of loading of a sample on a three-point bending (a) and FEM (b) Скорость нагружения была 2 мм/мин (квазистатическое нагружение). Механические характеристики ПММА: модуль упругости Е » 3500 ГПа, коэффициент Пуассона μ » 0,34, предел прочности при изгибе σBи » 120…130 МПа, а при растяжении σBp » 90…100 МПа, относительное удлинение при разрушении δ » 4…5 %. В работах [24, 25] характеристики ПММА близки к приведенным выше. Диаграммы «усилие - прогиб» (рис. 3), полученные при испытаниях образцов с надрезами, имели линейный характер вплоть до разрушения, что свидетельствовало о хрупком поведении ПММА. Рис. 3. Диаграммы изгиба образцов с надрезами (v = 25°) глубиной 0,25-5 мм Fig. 3. Diagrams of bending samples with notches (v = 25°) with the depths of 0.25-5 mm На рис. 4 приведены фотографии поверхностей разрушения образцов с двугранным углом v = 25° и глубиной L = 0,25…5 мм. Стрелками показаны зоны надрезов. L = 0,5 мм L = 1 мм L = 1,5 мм L = 2,5 мм L = 3 мм L = 3,5 мм L = 4,5 мм L = 5 мм Рис. 4. Фотографии поверхности разрушения Fig. 4. Pictures of fractured surfaces Отметим, что во время испытаний образцы с неглубокими надрезами (L £ 2 мм) разрушались на две части с хорошо слышимой акустической эмиссией, тогда как в образцах с более глубокими вырезами развитие трещины от надреза было достаточно медленным и заметным невооруженным глазом. У надрезов глубиной 0,5-1,5 мм на поверхностях разрушения имеется бороздчатая часть (быстрое прерывистое движение трещины). Длина бороздчатой зоны уменьшается с увеличением глубины надреза. Для глубоких надрезов (L > 2,5 мм) поверхность разрушения практически гладкая. На этих фотографиях можно отметить место начала разрушения на дне выреза (см. рис. 4, стрелки). Шаг бороздок в образцах практически неизменен и составляет (0,253 ± 0,01) мм. В зоне бороздок поверхность разрушения имеет существенно большую площадь по сравнению с гладкими картинами (L > 2,5 мм), что свидетельствует о рассеянии значительно большей удельной энергии деформации в таких образцах. Этот факт подтверждается диаграммами изгиба (см. рис. 3), на которых на разрушение образцов с надрезами меньшей глубины потрачена намного большая энергия, чем на образцы с большой глубиной вырезов. Определенный по картинам на рис. 4 шаг бороздок оказался тесно связан с размером площадки осреднения напряжений для нелокального критерия прочности. 2. Анализ напряженно-деформированного состояния образца с острым v-образным надрезом при изгибе На рис. 2, a приведена схема изгиба образцов с v-образным надрезом. Глубина выреза L изменялась в диапазоне 0,25-5 мм, двугранный угол v = 0…180°, высота h = 10 мм, толщина 3,9 мм. Расчет напряженного состояния проводили в плоской постановке, рассматривали лишь половину образца ввиду симметрии (рис. 2, б). Во всех примерах сетка КЭ имела сгущение в области вершины выреза (рис. 5) радиусом 0,25 мм с размером КЭ 0,005 мм. Рис. 5. Сгущение сетки КЭ для выреза v = 90° глубиной L = 0,5 мм Fig. 5. Densification of the FE grid for the notch v = 90° with depth of L = 0.5 mm Поскольку сингулярные задачи показывают стремление максимальных напряжений к бесконечности с уменьшением размера КЭ, стандартное доказательство сходимости можно осуществить, лишь сравнивая раcчетные значения коэффициента интенсивности и степени сингулярности с теоретическими значениями по формуле (5). Это сделано ниже, в п.3, с использованием решения Броека [12] и формул (10). Стоит отметить, что для обычных трещин (с нулевым двугранным углом) в коммерческих пакетах, реализующих метод конечных элементов (ANSYS, ABAQUS [20, 21]), имеются специальные сингулярные КЭ со степенью сингулярности 0,5. Однако для других случаев, с ненулевыми двугранными углами, таких КЭ нет. 3. Обработка данных расчетов Для вычисления коэффициента интенсивности К1v и степени сингулярности m в образцах с острыми надрезами и произвольным двугранным углом v в работах [24-26] предложена достаточно сложная в реализации методика вычисления G-интеграла у надреза и его связь с коэффициентом интенсивности напряжений. Проблема в том, что степень сингулярности здесь не равна 0,5. В данной статье используется выражение (3) для напряжений σx(r) у вершины острого надреза на расстоянии r. Две переменные и m предлагается искать из сопоставления напряженных состояний, полученных МКЭ в нескольких узлах на продолжении разреза в окрестности его вершины с теоретическим выражением (3). При этом анализ показывает, что в одном-двух первых, наиболее напряженных конечных элементах погрешность вычислений МКЭ еще велика, так как функция формы КЭ лишь квадратичная (плоский элемент второго порядка типа PLANE183, опция plane strain [20]) и сингулярность корректно отобразить не может. А на достаточно удаленном диапазоне расстояний зависимость (3) работает вполне адекватно (рис. 7). Уместно добавить, что в работе [34] использован аналогичный прием, в котором для вычисления коэффициентов интенсивности напряжений с учетом несингулярных слагаемых использованы узлы сетки конечных элементов, далеко отстоящие от вершины выреза. Минимизируя среднеквадратичные отклонения расчетных (3) и численных МКЭ (пакет ANSYS Workbench) значений, можно найти оптимальное сочетание и m с использованием программы, написанной в пакете MathCAD [22]. Для этой цели в окрестности вершины разреза выделена зона радиусом 0,25 мм, в которой сетка конечных элементов была более мелкой (размер 5 мкм, рис. 6, а). На линии - продолжении надреза - были выведены полученные расчетом МКЭ значения напряжений σx в узлах сетки КЭ (20 точек на рис. 6, б, из которых лишь 10-15 учитывались при оптимизации. Там же показана зависимость σx(r) по формуле (3) при различных сочетаниях m и (см. полосу разброса при изменении оптимальных параметров на ±10%). При меньшем количестве учитываемых точек получаемые оптимальные параметры сильно осциллировали, а при большем количестве точек на результат оптимизации (резкое снижение степени сингулярности) начинал влиять способ нагружения - изгиб - со снижением номинальных напряжений растяжения при удалении от вершины надреза. Таким образом, учет 10-15 точек давал наиболее стабильный результат оптимизации параметров напряженного состояния у выреза. Рис. 6. Сетка КЭ в зоне выреза (а) и распределение напряжений (б) вдоль линии 1-2. Глубина надреза 0,5 мм Fig. 6. FE grid in the area of the notch (a) and stress distribution (b) along the line 1-2. The depth of notch is 0.5 mm Рис. 7. Примеры распределения напряжений при найденных оптимальных m и для образцов с глубиной надреза L = 5 мм и разными двугранными углами v = 10…160° Fig. 7. Examples of stress distribution for given optimal m and for samples with a depth of cut L = 5 mm and different two-edge angles v = 10...160° На рис. 6 и 7 линии - расчет по (3) с оптимальными значениями и m. Получено, что степень сингулярности m не зависит от глубины выреза и хорошо согласуется со значениями, полученными из решения уравнения (6) (рис. 8). Рис. 8. Степень сингулярности m в функции двугранного угла и глубины надреза (L = 0,5…5 мм) Fig. 8. Degree of singularity of m in a function of an two-edge angle and a depth of notch (L = 0,5...5 mm) Здесь сплошная линия - численное решение трансцендентного уравнения (6) с помощью пакета MathCAD, прерывистая - удобная его аппроксимация для дальнейших расчетов m(v) = 0,5(1-7,486×10-9v3,603), (9) в которую угол v подставляется в градусах. Погрешность аппроксимации - менее 1 %. Коэффициенты интенсивности напряжений для различных глубин вырезов h и двугранных углов v при нагрузке Р = 200 Н (см. рис. 2) приведены на рис. 9 (значки - МКЭ, сплошные линии - аппроксимации, прерывистые - расчет по формуле (5)). Для расчета коэффициентов интенсивности предлагается использовать следующую зависимость, из которой получается известная формула Броека (5) заменой m(v) на 0,5: (10) В подтверждение правомочности использования (10) отметим, что в работе [24] показано, что для трехточечного изгиба образца с v-образным вырезом коэффициент интенсивности может быть выражен через номинальное напряжение в гладком образце σn = 1,5Pl/bh2, множитель hm и функцию геометрии образца F(x,…) откуда появляется что с точностью до константы согласуется с выражением (10). Отметим, что функция F(x,…) предполагает включение кроме относительной глубины x и угол раскрытия v граней выреза. Однако в данной работе использование функции Броека (5) с одним аргументом x для v-образных вырезов оказывается экспериментально оправданным. Рис. 9. Коэффициенты интенсивности напряжений для различных глубин вырезов L и двугранных углов v Fig. 9. The stress intensity factors for different depths of notches L and two-edge angles v 4. Расчеты на прочность Условие прочности (8) - нелокальный критерий прочности - требует вычисления наибольшего значения осредненного на базе А напряжения σx: откуда, используя (8), получаем нагрузку разрушения как функцию величины двугранного угла v и относительной глубины выреза x = L/h: (11) Результаты расчета нагрузок разрушения P образцов приведены на рис.10 вместе с экспериментальными данными для наилучшего сочетания локального предела прочности ПММА (SIG = 95 МПа) и базы осреднения A = 0,25 мм. Этот размер весьма близок к шагу 0,253 мм бороздок на поверхности разрушения (см. рис. 4), а предел прочности SIG - к пределу прочности на растяжение ПММА, полученных для образцов с двугранными углами 25 и 90°. Прогноз на двугранный угол 120° показан в виде кривой синего цвета на рис. 10. Совпадение с экспериментальными данными можно считать вполне удовлетворительным. Рис. 10. Результаты расчета нагрузок разрушения образцов по двухпараметрическому критерию (а) и по линейной механике разрушения (б) в функции глубины выреза Fig. 10. Results of calculation of the samples’ failure loads on two-parametric criterion (a) and to the linear mechanics of failure (b) in a function of a notch depth Использование линейной механики разрушения приводит к необходимости расчета критических значений коэффициентов интенсивности для каждого двугранного угла v. Так, при малых глубинах надрезов (h = 0,5…1 мм) для угла v = 25° = 60,7 МПа·мм0,5; для угла 90° - 75,4 МПа·мм0,459, а для угла 120° - 145 МПа·мм0,384, расчетные кривые и экспериментальные точки показаны на рис. 10, б. Здесь погрешности прогноза прочности меньше, чем в случае (а), так как для каждой кривой имеется собственный независимый параметр - критический коэффициент интенсивности напряжений. Выводы В работе рассмотрен элемент конструкции из хрупкого материала ПММА с острым надрезом с произвольным углом между сторонами, подвергнутый трехточечному изгибу. Для расчетной оценки прочности предложен двухпараметрический нелокальный критерий прочности, который избавляет от учета сингулярности в вершине надреза и сводит задачу механики разрушения к задаче конструкционной прочности. Нагрузки разрушения образцов с двугранными углами от 25 до 120° получены расчетом по двухпараметрическому критерию и экспериментально хорошо согласуются. При этом параметры предложенного критерия прочности оказались тесно коррелированными с шагом бороздок на поверхностях разрушения и пределом прочности на растяжение использованного ПММА соответственно. Это совпадение можно было бы посчитать случайным, так как в параметры A и SIG изначально не закладывались конкретные физические величины. Однако, если обратиться к работам, в которых для упругопластического материала рассматривалась предельная удельная энергия деформирования в качестве правой части критерия разрушения, то для линейно-упругого материала она пропорциональна квадрату предела прочности, т.е. полностью соответствует выражению нелокального критерия (8) данной работы. Использование МКЭ с оптимальной аппроксимацией зависимости напряжений от координаты в окрестности выреза позволило найти степень сингулярности и коэффициенты интенсивности напряжений для вырезов с двугранными углами до 180°. Предложены удобные аппроксимационные формулы, позволяющие проводить расчеты на прочность изгибаемых образцов из хрупкого материала с произвольными глубинами и углами острых надрезов без использования расчетов МКЭ.

About the authors

S B Sapozhnikov

South Ural State University

M A Ivanov

South Ural State University

S I Yaroslavtsev

South Ural State University

I A Shcherbakov

South Ural State University

References

Statistics

Views

Abstract - 301

PDF (Russian) - 117