Раздел 20 Анализ корреляции экспериментальных и расчетных результатов
Раздел 20. Анализ корреляции экспериментальных и расчетных результатов ВВЕДЕНИЕ В АНАЛИЗ КОРРЕЛЯЦИИ ЭКСПЕРИМЕНТАЛЬНЫХ И РАСЧЕТНЫХ РЕЗУЛЬТАТОВ..………………………….………………………………...…… 20 - 3 ИНТЕГРАЦИЯ РАСЧЕТОВ И ИСПЫТАНИЙ...…………………………………………..…… 20 - 5 ПЛАНИРОВАНИЕ ИСПЫТАНИЙ…………………………………………………………….… 20 - 6 АНАЛИЗ РЕЗУЛЬТАТОВ ИСПЫТАНИЙ..…………………………………………………… 20 - 9 ВЫВОДЫ.………………………………………………………………………………………… 20 - 18 УТОЧНЕНИЕ РАСЧЕТНОЙ МОДЕЛИ……………………………………………………….. 20 - 19 ЛИТЕРАТУРА.…………………………………………………………………………………… 20 - 27
Введение в анализ корреляции экспериментальных и расчетных результатов Результаты, полученные с помощью MSC.Nastran, и опытные данные могут не совпадать вследствие недостатков моделирования и проведения испытаний. Известно выражение: “Никто, за исключением инженера-аналитика, не верит результатам расчетов; все, кроме испытателя, верят экспериментальным данным” Источники возникновения трудностей моделирования: Физическая природа явления Граничные условия Свойства материалов Податливость соединений Задумано конструктором сделано производственником Демпфирование
Введение в анализ корреляции экспериментальных и расчетных результатов Цели анализа корреляции экспериментальных и расчетных результатов: Исследовать степень корреляции между результатами, полученными с помощью MSC.Nastran, и опытными данными. Уточнить модель для лучшей согласованности данных расчетов и испытаний. Инженер, исследующий корреляцию, должен понимать данные эксперимента, расчета и знать возможные источники ненадежности как в тех, так и в других.
Интеграция расчетов и испытаний Четыре этапа интеграции расчетов и испытаний: Планирование эксперимента (моделирование испытаний) Сбор (регистрация) данных, например, с помощью датчиков ускорений. Обработка и анализ опытных данных (например, определение форм колебаний по экспериментальным данным) Анализ результатов испытаний (оценка уровня согласованности результатов расчетов с помощью MSC.Nastran и опытных данных, уточнение модели для улучшения сходимости данных) Инженер-аналитик участвует в 1-м и 4-м этапах, этапы 2 и 3 – за инженером-испытателем.
Планирование испытаний Создание исходной модели для определения оптимального воздействия на объект и расположения датчиков. Для этого есть два метода: Моделирование и экспертный анализ Проверка кросс-ортогональности Моделирование и экспертный анализ: используя MSC.Nastran, моделируйте испытания и ищите способ воздействия и расположение датчиков для получения максимального отклика. Проверка кросс-ортогональности: На базе планируемого расположения датчиков создайте A-set. Затем, используя редуцирование Гайана, вычислите собственные формы во всем интересующем частотном диапазоне, нормализуйте их по массе. Назовите эти результаты ( t) “тест”-модами. Затем, удалив A-set, повторите вычисление собственных форм, выводя результаты только для тех степеней свободы, которые входили в A-set при предыдущем расчете. Назовите эти результаты ( a) “расчетными” модами.
Планирование испытаний При третьем запуске MSC.Nastran “прочитайте” данные и вычислите: Если планируемое расположение датчиков (A-set) адекватное, диагональные члены полученной матрицы будут равны 1, а внедиагональные – 0. Если внедиагональные члены матрицы не равны 0, то предлагаемое расположение датчиков неадекватное и должно быть пересмотрено. (На практике приемлемо, если внедиагональные члены меньше 0,05). В комплекте с MSC.Nastran поставляется процедура на языке DMAP (premaca.vxx) для выполнения планирования испытаний.
Планирование испытаний Для лучшего понимания поведения объекта испытаний в MSC.Nastran могут быть также вычислены эффективная модальная масса1 и модальная кинетическая энергия2. Для облегчения планирования испытаний и анализа их результатов расположение датчиков и их ориентация должны полностью соответствовать соответствующим допущениям при проведении расчетного анализа. Для выполнения этих требований в расчетной модели могут использоваться элементы RBAR, межузловые связи (MPC), альтернативные системы координат. После того, как определятся места расположения датчиков и приложения воздействий, испытания д.б. промоделированы для того, чтобы убедиться, что образец в процессе их проведения не будет перегружен (и, следовательно, поврежден). 1.Для i-ой моды эффективная масса равна где M – матрица масс и Dm – вектор колебаний жесткого тела. 2.Для i-ой моды кинетическая энергия равна
Анализ результатов испытаний Сравнение результатов расчетов, полученных с помощью MSC.Nastran, и результатов испытаний путем экспертной оценки или проверки кросс-ортогональности. Экспертная оценка с использованием графических возможностей Графики: Нарисуйте опытные и расчетные данные на одном графике и сравните степень их соответствия по величине и частотному составу. Убедитесь, что сравниваемые зависимости соответствуют по расположению и ориентации. Графические изображения конструкции: Нарисуйте расчетные о экспериментально полученные формы колебаний и сравните их на предмет соответствия (можно также использовать возможности анимации). Экспертным оценкам могут быть подвергнуты и “неграфические” результаты, такие, например, как значения резонансных частот.
Анализ результатов испытаний Проверка кросс-ортогональности Величина, подлежащая проверке на кросс-ортогональность при анализе корреляции опытных и расчетных данных, определяется как где t - экспериментально полученные моды Maa - матрица масс, соответствующая A-set a - моды, полученные расчетом для A-set Для соответствия расчетных и опытных данных внедиагональные члены должны быть меньше 10% от диагональных членов. Порядок чередования “расчётных” и “опытных” мод может быть нарушен, что приведет к большим значениям внедиагональных членов матрицы корреляции COR. При несоответствии опытных и расчетных данных может возникнуть необходимость уточнения расчетной модели.
Анализ результатов испытаний Пример анализа матрицы COR Рассмотрим плоскую модель балки. Положим, что датчики ускорений в направлениях X и Y располагаются в каждом четном узле.
Анализ результатов испытаний Частоты и формы колебаний, полученные при эксперименте
Анализ результатов испытаний Частоты и формы колебаний, полученные при помощи MSC.Nastran Частоты, полученные экспериментально и с помощью MSC.Nastran, близки
Анализ результатов испытаний Замечание: Расположение узлов, входящих в A-set, такое же, как и датчиков при испытаниях.
Анализ результатов испытаний
Анализ результатов испытаний
Анализ результатов испытаний
Выводы 1. Матрица COR указывает на хорошее совпадение форм (внедиагональные члены малы по сравнению с диагональными). 2. “Опытная” мода 5 - это “расчетная” мода 6 наоборот. На это указывают большие значения членов (5,6) и (6,5) матрицы COR. 3. Величина диагональных членов не равна 1,0 потому что “опытные” моды не нормализованы до единичной массы (они нормализованы до единичного максимального компонента). 4. Для подготовки к вводу опытных данных с помощью оператора DMI необходима их некоторая предварительная обработка. 5. Для анализа большого количества мод может быть разработана дополнительная процедура на языке DMAP для упрощения вида результатов расчета кросс-корреляции. 6. Процедура на языке DMAP для выполнения проверки кросс- корреляции (postmaca.vxx) – одна из поставляемых частей MSC.Nastran.
Уточнение расчетной модели Существует три пути улучшения расчетной модели для достижения соответствия расчетных и опытных данных: 1. Использование “грубой” силы 2. Использование матрицы чувствительности 3. Использование методики оптимизации Все методы подразумевают изменение параметров модели таких как I и A для BAR’ов и t для QUAD’ов. Упругость основания также может рассматриваться в качестве параметра, если она в явном виде представлена в модели (например, с помощью ELAS’ов). Метод “грубой” силы: Изменения модели производятся на базе экспертных оценок и результатов “прикидочных” расчетов. После изменения модели расчет с помощью MSC.Nastran повторяется и его результаты анализируются на предмет лучшего согласования с опытными данными. Если согласие данных не улучшается, - новые изменения и т.д.
Уточнение расчетной модели Матрица чувствительности: Чувствительность – это градиент отклика по параметру модели. “Отклики” – это результаты вычислений с помощью MSC.Nastran, а “параметры” – значения свойств. Каждое значение в матрице чувствительности S определяется как где Ri - i-ый отклик и Pj - j-ый параметр. Чем больше величина Sij, тем больше чувствительность отклика. SOL 200 в MSC.Nastran вычисляет матрицу чувствительности. Матрица S может использоваться для оценки того, какой параметр необходимо изменить для изменения соответствующего отклика.
Уточнение расчетной модели Матрица S также может быть получена и использована посредством метода наименьших квадратов для минимизации разницы между результатами испытаний и расчетов: Pn = Po + (STS)-1ST(Rt - Ra) где Pn - новое значение параметра, подлежащего изменению Po - значение параметра до изменения S - матрица чувствительности Rt - “опытный” отклик Ra - “расчетный” отклик
Уточнение расчетной модели Использование оптимизации: “Инструмент” - SOL 200 Цель: минимизация разницы между опытными и расчетными результатами, но при минимальных изменениях в исходной модели Оператор DEQATN используется для записи целевой функции в виде “взвешенной” разности между опытными и расчетными результатами. где RTi - i-ый результат испытаний PFj - j-ый параметр модели после уточнения RAi - i-ый расчетный результат POj - исходное значение j-го параметра WRi - весовой коэффициент для i-го испытания WPj - весовой коэффициент для j-го параметра wt - весовой коэффициент для всех опытных данных wp - весовой коэффициент для всех параметров SOL 200 минимизирует целевую функцию (E) с учётом ограничений.
Уточнение расчетной модели Пример Модель корпуса дисковода, включающая 1406 узлов, 1354 оболочечных элементов,четыре переменных проектирования (толщины оболочек) и четыре упругих моды, полученные в результате эксперимента Необходимо минимизировать разницу между частотами собственных колебаний, вычисленными с помощью MSC.Nastran и полученными экспериментально.
Уточнение расчетной модели Пример Весовой коэффициент для опытных данных в 100 раз больше весового коэффициента для параметров. Учитываются верхние и нижние ограничения на допустимую толщину оболочечных элементов. После 5 итераций (шагов) достигнута сходимость Результаты:
Уточнение расчетной модели
Уточнение расчетной модели SOL 200 может также использоваться для улучшения корреляции с экспериментальными результатами по определению частотного отклика. Возможно использование для улучшения корреляции статических результатов. Для достижения лучших результатов необходимо использовать как можно больше опытных данных (включая вес конструкции). В качестве оптимизируемых необходимо рассматривать только те параметры, в значениях которых действительно есть неуверенность.
Литература Blakely, K. and Dobbs, M., “Integrated System Identification: The Union of Testing and Analysis,” Proceedings First International Modal Analysis Conference, November 1982. Kientzy, D., Richardson, M. and Blakely, K., “Using Finite Element Data to Set Up Modal Tests,” Sound and Vibration, June 1989. Blakely, K., “Updating MSC.Nastran Models to Match Test Data,” Proceedings MSC World Users’ Conference, March 1991. Blakely, K. and Rose, T., “Cross-Orthogonality Calculations for Pre-Test Planning and Model Verification,” Proceedings MSC World Users’ Conference, May 1993. Blakely, K., “Matching Frequency Response Test Data with MSC.Nastran,” Proceedings MSC World Users’ Conference, June 1994.