Практические методы и руководство по вычислению параметров трещиностойкости конструкций с применением систем MSC.Software презентация

систем MSC.Software

Внедрение прогрессивных методов ресурсного проектирования деталей авиационных ГТД с использованием продуктов MSC Часть 1

Слезкин Д.В., представительство MSC.Software в СНГ
IX Российская конференция пользователей систем MSC

разрушения различных конструкций в США составляют порядка 4% от ВВП!!!

П. И. Баранова разработаны новые нормы прочности, требующие проводить прочностные расчеты с учетом начальных дефектов в деталях авиационных ГТД
НТЦ им. Люльки – необходимость при расчетах учитывать возможное наличие начальных дефектов элементов ГТД
Требуются ответы на следующие вопросы:
какие методы вычисления Коэффициентов Интенсивности Напряжений (КИН) использовать?
как вычислять параметры трещиностойкости с помощью программного обеспечения MSC.Software (пошаговая инструкция)?
какие требования к параметрам КЭ сетки?
Обязательное требование:
ТЕСТИРОВАНИЕ НА ЗАДАЧАХ, ИМЕЮЩИХ ИЗВЕСТНОЕ РЕШЕНИЕ

к трем типам: расчет плоских элементов конструкции с трещинами, расчет осесимметричных элементов конструкции с трещинами, а также расчет сложных пространственных элементов конструкции с эллиптическими трещинами

моделей с трещинами
Пластина с трещиной. Изучение методов вычисления Коэффициентов Интенсивности Напряжений (КИН) с помощью Nastran. Определение необходимого размера КЭ сетки в области вершины трещины на стандартном образце (пластина с краевой трещиной).
Осесимметричная задача. Изучение разных методов для вычисления КИН с помощью Nastran, Marc
Трехмерная задача. Создание подробной методики вычисления J-интеграла с помощью Patran/Marc для краевых и эллиптических трещин

расчет (расчеты велись для плоского образца с краевой трещиной)
Вычисление КИН 1 по стандартной справочной формуле 1 – использован соответствующий справочник [3], линейный расчет удаленного поля напряжений
Вычисление КИН 1 по формуле 2 - с учетом поправки на пластичность по Ирвину, линейный расчет (расчеты цилиндрического образца с круговой трещиной велись в осесимметричной постановке)
Вычисление КИН 1 по перемещениям раскрытия берегов трещины (формула 3) –линейный расчет (расчеты цилиндрического образца с круговой трещиной велись в осесимметричной постановке)
Вычисление КИН 1 по раскрытию вблизи вершины (формула 4) – нелинейно-упругий расчет (расчеты цилиндрического образца с круговой трещиной велись в осесимметричной постановке)
Вычисление КИН 1 с помощью J-интеграла – формула (5) (расчеты велись для объемного образца с эллиптической трещиной)

вычисления сингулярных полей напряжений/деформаций (CRAC2D/CRAC3D) и вычисления Коэффициентов Интенсивности Напряжений (КИН) для произвольной геометрии (и нагружения) тела с трещиной;

CRAC2D

CRAC3D

справочное значение KI = 2.35 кгс/мм3/2 [3]
Более подробно – см. [7],

ВЫЧИСЛЕНИЕ КИН С ПОМОЩЬЮ NASTRAN

2.35 кгс/мм3/2 [3]
Более подробно – см. [7]

ВЫЧИСЛЕНИЕ КИН С ПОМОЩЬЮ NASTRAN

Важно: применение CRAC элементов ограничено наличием малых зон пластики, при развитых пластических зонах следует применять другие подходы

Изменение К1 в зависимости от размера элемента

разных методов для цилиндрического образца

D=10 мм., d=5 мм., c=2.5 мм., от σ=10 кгс/мм2
Характеристики материала при t=20° C: Е=20500 кгс/мм2, μ=0.3, σ0.2=68 кгс/мм2 , К1с=340 кгс/мм3/2

разрушения

Анализ области применимости линейной механики разрушения в зависимости от интенсивности внешней нагрузки (характеризующей размеры пластической зоны вблизи вершины трещины)

K1lin(r) – КИН1, вычисленный по формуле (2) с учетом поправки по Ирвину, K1(J) – КИН1, вычисленный с помощью J-интеграла, K1(delta) – КИН1, вычисленный по раскрытию вблизи вершины трещины, K1(linear) – КИН1, вычисленный по справочной формуле, K1(Uy) – КИН1, вычисленный по перемещениям раскрытия берегов трещины

K1c

мм., t = 20 мм., b = 80 мм., H = 120 мм., 0 < L/a < 1.0, 0 < L/t =< 0.4. Характеристики материала при t=20° C: Е=20500 кгс/мм2, μ=0.3, σ0.2=68 кгс/мм2 , К1с=340 кгс/мм3/2.Нгаружение σ=50 кгс/мм2

Справочное значение KI = 186.7 кгс/мм3/2 [3]

значение J-интеграла: J=1.45
Соответственно, КИН1:

Конечноэлементная модель - 13167 степеней свободы, размер элемента вдоль фронта трещины = 1.22 мм.

Проведены расчеты для КЭ сеток разного размера:
Размер элемента вдоль фронта трещины 0.61 мм., К1=183.8, 37359 ст.св.
Размер элемента вдоль фронта трещины 0.245 мм., К1=184.4, 267000 ст.св.

разрушения

НЕЛИНЕЙНАЯ МЕХАНИКА РАЗРУШЕНИЯ

J-интеграл

Вычисление J-интеграла с помощью метода Де-Лоренци
Линейный и нелинейный материалы
Учет больших деформаций
Температурные нагрузки
Возможность учета контакта между берегами трещины, включая трение
Автоматически распознаются граничные условия симметрии
Нагружение свободных поверхностей трещины
Вычисление динамических КИН

разрушения

ЛИНЕЙНАЯ МЕХАНИКА РАЗРУШЕНИЯ

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

трещин и J-интеграла с помощью систем MSC
На основе полученных решений разработаны и опробованы специалистами НТЦ им. Люльки пошаговые руководства по вычислению КИН и J-интеграла на примере расчета тестовых образцов (плоская пластина с трещиной, цилиндрический образец с круговой трещиной, массивный образец с эллиптической трещиной)
Описанные в работе подходы применимы не только при расчете элементов конструкции ГТД, но и при расчете деталей других конструкций из различных материалов
Накоплен важный опыт решения задач трещиностойкости, которым могут воспользоваться специалисты других предприятий по соответствующему запросу руководителей (необходимо быть официальным пользователем систем MSC.Software)
Выявлены области применения каждого из использованных в работе методов определения КИН
На основе анализа результатов решений полученных методом конечного элемента (МКЭ) можно рекомендовать, как достаточно оптимальный, диапазон размеров конечных элементов в вершине трещины от 0.1 мм до 1 мм. для вычисления КИН и J-интеграла

Parekh, Rocky R. Arnold, Peter J. Woytowitz, Anamet Laboratories Inc., San Carlos, California 94070; Presented at the MSC/NASTRAN User’s conference, Universal City, California 91608, March 20-21, 1986;
Hughes T. J. R., Akin J. E., “Techniques for Developing Special Finite Element Shape Functions with Particular Reference to Singularities”, International Journal for Numerical Methods in Engineering, Vol. 15, pp. 733-751, 1980.
Справочник по коэффициентам интенсивности напряжений: в 2-х томах. Ито Ю., Мураками Ю., Хасебэ Н., Юуки Р., Тоя М., Того К., Мията Х., Терада Х., Миядзаки Н., Аоки С. Перевод с английского. – М.: Мир, 1990, 448 стр.
О. Зенкевич, Метод конечных элементов в технике. Перевод на русский язык, М.: Мир, 1975, 440 стр.
В. М. Пестриков, Е. М. Морозов Механика разрушения твердых тел: курс лекций. – СПб.: Профессия, 2002. – 320 стр.
Хеллан К. Введение в механику разрушения. Пер. с англ. – М.:Мир, 1998. – 364 стр.
Слезкин Д. В. Методика вычисления Коэффициентов Интенсивности Напряжений с применением MSC.Nastran, 2006 год. Методическое пособие MSC.Software Moscow

СПИСОК ЛИТЕРАТУРЫ