Семинар NAS101 | 2006 | MSC.Software Corporation Постоянное представительство в СНГ Москва Собственные частоты и формы колебаний Раздел 6

2 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Рассмотрим систему с одной степенью свободы без демпфирования, как показано на рисунке: гдеm- масса k- жесткость Уравнение движения при свободных колебаниях (т.е. без внешней нагрузки или демпфирования) выглядит так: или Основные уравнения

3 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Основные уравнения (продолжение) Для системы с несколькими степенями свободы уравнение примет вид: где[K]- Матрица жесткости (такая же как и в статическом анализе). [M]- Матрица масс (представляет инерционные свойства конструкции). Матрицы [K] и [M] должны быть вещественными и симметричными. Запомните: Число степеней свободы равно числу координат, необходимых для описания деформированной формы конструкции в каждый данный момент времени.

4 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Матрица масс представляет инерционные свойства конструкции. MSC Nastran предоставляет пользователю две возможности для определения матрицы масс: Несогласованная матрица масс (используется по умолчанию) Содержит только диагональные члены, связанные с поступательными степенями свободы. Согласованная матрица масс Содержит также вне диагональные члены, связывающие поступательные степени свободы и вращательные степени свободы. (Примечание: для стержневых элементов связываются только степени свободы поступательного движения.) Матрица масс

5 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Матрица масс (продолжение) Пример матрицы масс где - массовая плотность A- площадь сечения Несогласованная матрица масс Согласованная матрица масс

6 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Матрица масс (продолжение) Согласованная и несогласованная матрицы масс Согласованная матрица масс, в общем случае, более точная чем несогласованная матрица масс. Несогласованная матрица масс предпочтительнее с точки зрения скорости расчета при динамическом анализе. Выбираются пользователем согласованные матрицы масс для элементов: PARAM,COUPMASS,1 выбирает согласованную матрицу масс для всех элементов BAR, ROD, и PLATE которые включают в себя изгибную жесткость. По умолчанию выбирается несогласованная матрица масс. Элементы, которые могут иметь как согласованную, так и несогласованную матрицу масс: BAR, BEAM, CONROD, HEXA, PENTA, QUAD4, QUAD8, ROD, TETRA, TRIA3, TRIA6, TRIAX6, TUBE

7 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Матрица масс (продолжение) Элементы имеющие только несогласованную матрицу масс: CONEAX, SHEAR Элементы имеющие только согласованную матрицу масс: BEND, HEX20, TRAPRG, TRIARG Несогласованная матрица масс содержит только диагональные компоненты, связанные с поступательным движением (без вращения). Согласованная матрица масс содержит вне диагональные компоненты, связанные как с поступательным движением, так и с вращательным для элементов типа BAR, BEAM, и BEND.

8 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Рассмотрим уравнение (6-1) Представим гармоническое решение в виде (6-2) (Физический смысл этого в том, что по всем координатам осуществляется синхронное движение и конфигурация системы не меняет свою форму во время движения, меняется только амплитуда.) Из уравнения 6-2 получим (6-3) i Теория

9 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Теория (продолжение) Подставив уравнения (6-2) и (6-3) в уравнение (6-1), получим: которое, после преобразования примет вид: (6-4) Это обобщенная проблема на собственные значения.

10 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Теория (продолжение) Поэтому, существует два варианта решения данного уравнения: 1. Если det ( [ K ] – 2 [ M ] ) 0, то из уравнения (6-4) следует, что есть только одно решение: { } = 0 это, так называемое, тривиальное решение, которое с физической точки зрения интереса не представляет 2. Тогда, чтобы получить нетривиальное решение относительно { }, необходимо, чтобы: det ( [ K ] – 2 [ M ] ) = 0

11 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Теория (продолжение) Проблема собственных значений сводится к решению уравнения det ( [ K ] – 2 [ M ] ) = 0 или det ( [ K ] – [ M ] ) = 0 где = 2

12 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Если конструкция имеет N динамических степеней свободы (степеней свободы с массой), то имеется N чисел, которые будут решением проблемы собственных значений. Эти ( 1, 2,..., N ) и есть собственные частоты, так же известные как резонансные частоты, основные частоты и частотные характеристики. Собственный вектор связанный с собственной частотой и называется собственной формой или формой колебаний. Собственная форма соответствует отклонению модели от состояния покоя. Когда конструкция вибрирует, ее форма в каждый момент времени является линейной комбинацией собственных форм колебаний. Теория (продолжение)

13 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Теория (продолжение) Пример: Свободно опертая балка Форма 1 Форма 2 Форма 3 Форма 2 и т. д.

14 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Для определения динамических характеристик конструкции. Например, если вращающийся механизм помещается в некоторую конструкцию, то для избежания чрезмерных вибраций желательно проверить насколько близки его частоты вращения к собственным частотам конструкции. Для определения динамического усиления нагрузки. Собственные частоты и формы используются для того, чтобы знать как проводить последующий динамический анализ (нестационарные процессы, спектральный анализ), например, если требуется выбрать подходящий временной шаг t для интегрирования уравнения движения в нестационарном анализе. Для чего нужно вычислять собственные частоты и формы колебаний

15 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Для чего нужно вычислять собственные частоты и формы колебаний (продолжение) Для использования собственных частот и форм при последующем динамическом анализе, например, в анализе нестационарных процессов с использованием модальных методов. Для проведения экспериментальных исследований, например, для определения местонахождения акселерометров и т. д.

16 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Важные замечания относительно анализа собственных частот и форм колебаний Если конструкция не полностью закреплена, например, допускается перемещение ее как твердого тела (форма колебаний без напряжений) или механизма, то как минимум одна собственная частота будет равняться нулю. Пример: Следующая не закрепленная конструкция может перемещаться как твердое тело.

17 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Важные замечания относительно анализа собственных частот и форм колебаний (продолжение) Собственные частоты ( 1, 2,...,) имеют размерность радиан/с. Они также могут быть представлены в герцах (цикл/с), с помощью выражения

18 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Важные замечания относительно анализа собственных частот и форм колебаний (продолжение) Масштабирование собственных форм произвольное. Например: представляют одну и ту же "форму колебаний"

19 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Важные замечания относительно анализа собственных частот и форм колебаний (продолжение) Определение собственных частот, т. е., решение уравнения является сложной задачей. Для решения этой задачи должны применяться численные методы.

20 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Методы вычислений MSC Nastran предоставляет пользователю следующие три подхода к нахождению собственных значений: Метод итераций Собственные значения (или собственные частоты) определяются одно за другим методом итераций. Имеются два варианта метода INV и SINV. Эти методы наиболее удобны при определении нескольких собственных частот. В общем SINV более надежен, чем INV. Метод трансформаций Общая проблема собственных значений

21 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Методы вычислений (продолжение) преобразуется в форму: Затем матрица [ A ] трансформируется в трех диагональную матрицу, с использованием методов Гивенса или Хаусхольдера. В конце концов, все собственные значения извлекаются одновременно с использованием QR алгоритма. Имеются два варианта метода Гивенса и два варианта метода Хаусхольдера: GIV, MGIV, HOU и MHOU. Эти методы более эффективны, когда требуется определить большое количество собственных значений

22 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Методы вычислений (продолжение) Метод Ланцоша Это рекомендуемый метод, который является комбинацией методов итераций и трансформаций. Этот метод наиболее эффективен для расчета малого числа собственных значений в больших разреженных задачах (большинство КЭ моделей попадают в эту категорию).

23 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Записи для анализа собственных частот Секция EXECUTIVE CONTROL SOL 103 Секция CASE CONTROL METHOD=x Номер указывающий на записи EIGR или EIGRL, имеющиеся в секции BULK DATA. С дополнительным множителем могут быть использованы сразу несколько SUBCASE Секция BULK DATA EIGR Запись для нахождения собственных значений или EIGRL Запись для метода Ланцоша. Требуются массовые свойства.

24 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Массовые свойства Конструкционная масса Добавляет массу элементов (например, для расчета гравитационного эффекта) Не конструкционная масса Добавляет массу (пример – нагрузка на строительное перекрытие, нагрузка на грузовой корабль) MAT1MIDEGNURHO MAT Плотность в записях MATi, единицы = (масса/объем)

25 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Массовые свойства (продолжение) Масса на единицу размера (масса на единицу площади в данном случае) Сосредоточенная масса Явные массовые свойства в узле (CONM2) (отступ центра тяжести сосредоточенной массы от узловой точки, моменты, инерционные характеристики) PSHELLPIDMID1TMID2121/T 3 MID3TS/TNSM PSHELL

26 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Массовые свойства (продолжение) Единицы массы Программа подразумевает что размерности инерции: PARAM,WTMASS умножает входные данные для получения размерности инерции. Обычно используется для замены весовых размерностей на массовые. Пример: Плотность (RHO) стали задана как lb/ft 3 в MAT1. Используйте PARAM,WTMASS, , который умножит элементы матрицы масс конструкции на 1/g (= 1/ ft/sec 2 ) для того, чтобы перевести плотность в соответствующие единицы инерции. lb-sec 2 /ft (ft-lb-sec система) кг-сек 2 /м

27 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Массовые свойства (продолжение) Генератор весов узловых точек Матрица преобразования из основной системы координат в систему главных осей инерции Главные массы и соответственно центр тяжести Матрица масс жесткого тела (MО) относительно ссылочной точки в основной системе координат. Матрица инерции I(S) в центре тяжести в системе главных осей масс. Матрица инерции I(Q) в центре тяжести в системе главных осей инерции. Матрица преобразования [Q] из S-осей и Q-осей. Примечание: 1. Это стандартный вывод генератора весов узловых точек. Он получается с помощью задания в параметре GRDPNT целого числа, значение которого определяет номер ссылочного узла сетки. Если это число 0 или несуществующий узел, то в качестве ссылочной точки выбирается начало отсчета базовой системы координат. PARAM,GRDPNT,0

28 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Запись SUPORT Запись SUPORT в секции Bulk Data Для вычисления форм колебаний конструкции как твердого тела используются специальные подпрограммы. Esthetics Абсолютные нулевые собственные значения вместо вычисления нулей (для всех методов кроме метода Ланцоша, где программа будет "судить": либо собственное значение должно быть 0.0, либо нет). Cost Отдельная подпрограмма, используемая для расчета форм колебаний конструкции как жесткого тела, может значительно уменьшить время CPU SUPORTIDC C C C SUPORT16125

29 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Запись SUPORT (продолжение) Замечания: Статистически определенный набор закреплений Достаточное число закреплений, поддерживающее все моды твердого тела Метод Ланцоша использует вычисленные собственные вектора

30 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Записи для анализа собственных частот (продолжение) EIGRL Рекомендуемая запись для расчета собственных частот Описывает данные, необходимые для определения собственных частот или анализа устойчивости методом Ланцоша. Поле Содержание SIDИдентификатор (уникальное целое положительное число) EIGRLSIDV1V2NDMSGLVLMAXSETSHFSCLNORM EIGRL

31 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Записи для анализа собственных частот (продолжение) V1, V2Частотный анализ: интересующий диапазон частот Анализ устойчивости: интересующий диапазон (V 1 <V 2, вещественные числа). Если нужны все формы ниже некоторой частоты, задайте V 2 и оставьте поле V 1 пустым. Не рекомендуется задавать 0.0 для V 1, более эффективно использовать маленькое отрицательное число или поле его пустым NDТребуемое число корней (целое положительное число или пробел) MSGLVLУровень диагностики (целое число от 0 до 3 или пробел) MAXSETЧисло векторов в блоке (целое число от 1 до 15 или пробел)

32 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Записи для анализа собственных частот (продолжение) EIGRL Рекомендуемая запись для расчета собственных частот SHFSCLОценка первой из собственных частот (действительное число или пробел) NORMМетод нормирования собственных векторов, одно из BCD значений либо "MASS", либо "MAX". MASSНормирование по обобщенной массе масс (по умолчанию) MAXНормирование по наибольшей компоненте

33 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Записи для анализа собственных частот (продолжение) Основываясь на входных данных, программа будет: Рассчитывать все корни ниже V2 (V1 = пустое, V2 = наибольшее интересующее значение собственной частоты, ND - пробел) Рассчитывать максимум ND корней между V1 и V2 (V1, V1, ND - заданы) Рассчитывать ND корней выше V1 (V1 = самое низкое интересующее значение частоты, V2 = пустое, ND = требуемое число корней) Рассчитывать первые ND корней (V1 и V2 пустые, ND = требуемое число корней). Рассчитывать все корни между V1 и V2 (V1 = самое низкое интересующее значение частоты, V2 = наибольшее интересующее значение собственной частоты, ND = пустое)

Семинар NAS101 | 2006 | MSC.Software Corporation Постоянное представительство в СНГ Москва Пример 8 Анализ собственных частот подкрепленной панели

35 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Анализ собственных частот подкрепленной панели (продолжение)

36 MSC Confidential Семинар NAS101 | 2006 | Раздел 6 | MSC.Software Corporation Анализ собственных частот подкрепленной панели (продолжение) Описание модели Та же самая модель подкрепленной панели, которая была использована в Примере 5. Рассчитайте первые 6 собственных значений частот. Проверьте корректность используемых единиц измерения