NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Раздел 20 Анализ корреляции экспериментальных и расчетных результатов

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Раздел 20. Анализ корреляции экспериментальных и расчетных результатов ВВЕДЕНИЕ В АНАЛИЗ КОРРЕЛЯЦИИ ЭКСПЕРИМЕНТАЛЬНЫХ И РАСЧЕТНЫХ РЕЗУЛЬТАТОВ..………………………….………………………………...…… ИНТЕГРАЦИЯ РАСЧЕТОВ И ИСПЫТАНИЙ...…………………………………………..…… ПЛАНИРОВАНИЕ ИСПЫТАНИЙ…………………………………………………………….… АНАЛИЗ РЕЗУЛЬТАТОВ ИСПЫТАНИЙ..…………………………………………………… ВЫВОДЫ.………………………………………………………………………………………… УТОЧНЕНИЕ РАСЧЕТНОЙ МОДЕЛИ……………………………………………………… ЛИТЕРАТУРА.……………………………………………………………………………………

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Введение в анализ корреляции экспериментальных и расчетных результатов q Результаты, полученные с помощью MSC.Nastran, и опытные данные могут не совпадать вследствие недостатков моделирования и проведения испытаний. q Известно выражение: Никто, за исключением инженера- аналитика, не верит результатам расчетов; все, кроме испытателя, верят экспериментальным данным q Источники возникновения трудностей моделирования: äФизическая природа явления äГраничные условия äСвойства материалов äПодатливость соединений äЗадумано конструктором сделано производственником äДемпфирование

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Введение в анализ корреляции экспериментальных и расчетных результатов q Цели анализа корреляции экспериментальных и расчетных результатов: äИсследовать степень корреляции между результатами, полученными с помощью MSC.Nastran, и опытными данными. äУточнить модель для лучшей согласованности данных расчетов и испытаний. q Инженер, исследующий корреляцию, должен понимать данные эксперимента, расчета и знать возможные источники ненадежности как в тех, так и в других.

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Интеграция расчетов и испытаний q Четыре этапа интеграции расчетов и испытаний: Планирование эксперимента (моделирование испытаний) Сбор (регистрация) данных, например, с помощью датчиков ускорений. Обработка и анализ опытных данных (например, определение форм колебаний по экспериментальным данным) Анализ результатов испытаний (оценка уровня согласованности результатов расчетов с помощью MSC.Nastran и опытных данных, уточнение модели для улучшения сходимости данных) q Инженер-аналитик участвует в 1-м и 4-м этапах, этапы 2 и 3 – за инженером-испытателем.

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Планирование испытаний q Создание исходной модели для определения оптимального воздействия на объект и расположения датчиков. Для этого есть два метода: äМоделирование и экспертный анализ äПроверка кросс-ортогональности q Моделирование и экспертный анализ: используя MSC.Nastran, моделируйте испытания и ищите способ воздействия и расположение датчиков для получения максимального отклика. q Проверка кросс-ортогональности: На базе планируемого расположения датчиков создайте A-set. Затем, используя редуцирование Гайана, вычислите собственные формы во всем интересующем частотном диапазоне, нормализуйте их по массе. Назовите эти результаты ( t ) тест-модами. Затем, удалив A-set, повторите вычисление собственных форм, выводя результаты только для тех степеней свободы, которые входили в A-set при предыдущем расчете. Назовите эти результаты ( a ) расчетными модами.

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Планирование испытаний q При третьем запуске MSC.Nastran прочитайте данные и вычислите: q Если планируемое расположение датчиков (A-set) адекватное, диагональные члены полученной матрицы будут равны 1, а вне диагональные – 0. Если вне диагональные члены матрицы не равны 0, то предлагаемое расположение датчиков неадекватное и должно быть пересмотрено. (На практике приемлемо, если вне диагональные члены меньше 0,05). qВ комплекте с MSC.Nastran поставляется процедура на языке DMAP (premaca.vxx) для выполнения планирования испытаний.

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Планирование испытаний q Для лучшего понимания поведения объекта испытаний в MSC.Nastran могут быть также вычислены эффективная модальная масса 1 и модальная кинетическая энергия 2. q Для облегчения планирования испытаний и анализа их результатов расположение датчиков и их ориентация должны полностью соответствовать соответствующим допущениям при проведении расчетного анализа. Для выполнения этих требований в расчетной модели могут использоваться элементы RBAR, межузловые связи (MPC), альтернативные системы координат. q После того, как определятся места расположения датчиков и приложения воздействий, испытания д.б. промоделированы для того, чтобы убедиться, что образец в процессе их проведения не будет перегружен (и, следовательно, поврежден). 1. Для i-ой моды эффективная масса равна где M – матрица масс и Dm – вектор колебаний жесткого тела. 2. Для i-ой моды кинетическая энергия равна / / Ф i ) i 2

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Анализ результатов испытаний q Сравнение результатов расчетов, полученных с помощью MSC.Nastran, и результатов испытаний путем экспертной оценки или проверки кросс-ортогональности. q Экспертная оценка с использованием графических возможностей äГрафики: Нарисуйте опытные и расчетные данные на одном графике и сравните степень их соответствия по величине и частотному составу. Убедитесь, что сравниваемые зависимости соответствуют по расположению и ориентации. äГрафические изображения конструкции: Нарисуйте расчетные о экспериментально полученные формы колебаний и сравните их на предмет соответствия (можно также использовать возможности анимации). q Экспертным оценкам могут быть подвергнуты и неграфические результаты, такие, например, как значения резонансных частот.

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Анализ результатов испытаний q Проверка кросс-ортогональности Величина, подлежащая проверке на кросс-ортогональность при анализе корреляции опытных и расчетных данных, определяется как где t -экспериментально полученные моды M aa -матрица масс, соответствующая A-set a -моды, полученные расчетом для A-set Для соответствия расчетных и опытных данных вне диагональные члены должны быть меньше 10% от диагональных членов. Порядок чередования расчётных и опытных мод может быть нарушен, что приведет к большим значениям внедиагональных членов матрицы корреляции COR. q При несоответствии опытных и расчетных данных может возникнуть необходимость уточнения расчетной модели.

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Анализ результатов испытаний q Пример анализа матрицы COR Рассмотрим плоскую модель балки. Положим, что датчики ускорений в направлениях X и Y располагаются в каждом четном узле.

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Анализ результатов испытаний Частоты и формы колебаний, полученные при эксперименте

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Анализ результатов испытаний Частоты и формы колебаний, полученные при помощи MSC.Nastran Частоты, полученные экспериментально и с помощью MSC.Nastran, близки

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow ID PRETEST, DYNOTES SOL 103 TIME 15 COMPILE MODERS, SOUIN=MSCSOU, NOLIST, NOREF ALTER mxx.*phix$ MATPRN MXX,,,,//$ MATPCH MXX,PHIX,,,// $ CEND TITLE = MSC.Nastran MSC/XL SUBTITLE = MODES CASE CONTROL LABEL = DEFAULT SUBCASE STRUCTURE DISP = ALL SPC=1 METHOD=100 BEGIN BULK PARAM AUTOSPC YES GRDSET,,,,,,,345 BAROR,,,,, 0., 1., 0. GRID GRID GRID GRID GRID GRID GRID GRID GRID GRID GRID CBAR CBAR CBAR CBAR CBAR CBAR CBAR CBAR CBAR CBAR $ REDUCE TO TEST DOF ASET1, 12, 2, 4, 6, 8, 10 SPC PBAR MAT EIGRL ENDDATA Анализ результатов испытаний Замечание: Расположение узлов, входящих в A-set, такое же, как и датчиков при испытаниях.

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Анализ результатов испытаний ID COR, DYN.NOTES TIME 30 SOL 100 COMPILE USERDMAP, SOUIN=MSCSOU, NOLIST, NOREF ALTER 2 $ $ $ DMAP TO COMPUTE CROSS-ORTHOGONALITY $ INPUTS FROM MSC.Nastran RUN: MXX (A-SET MASS) $ PHIX (A-SET MODE SHAPES) $ (PREVIOUS M/N RUN USED MATPCH TO PUNCH DMI ENTRIES) $ INPUT FROM TEST: PHITEST (A-SET MODE SHAPES) $ OUTPUTS: UNITCHK (UNIT MASS CHECK) $ COR (CROSS-ORTHOGONALITY MATRIX) $ $ READ DMI INPUT DMIIN DMI,DMINDX/PHIX,PHITEST,MXX,,,,,,,/ $ $ VERIFY INPUT MATRICES MATPRN PHIX,PHITEST,MXX,,// $ $ MULTIPLY PHIX(TRANS)*MXX = PHITMASS MPYAD PHIX,MXX,/ PHITMASS /1///$ $ MULTIPLY PHITMASS*PHIX = UNITCHK MPYAD PHITMASS,PHIX,/ UNITCHK // $ $ PRINT TITLE AND UNITCHK MATPRN UNITCHK,,,,// $ MESSAGE // CHECK ON UNIT MASS/ $ $ MULTIPLY PHITMASS*PHITEST = COR MPYAD PHITMASS,PHITEST,/ COR // $ $ PRINT TITLE AND COR MATPRN COR,,,,// $ MESSAGE // CROSS-ORTHOGONALITY MATRIX/ $ ENDALTER CEND TITLE = CROSS-ORTHOGONALITY CHECK BEGIN BULK DMI,PHITEST,0,2,1,0,,10,6 DMI,PHITEST,1,1,0.,.0,0.,.15,0.,,.40,0.,.70,0.,1.0 DMI,PHITEST,2,1,.15,0.,.45,0.,.70,,0.,.90,0.,1.0,0... (rest of PHITEST).

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow DMI MXX DMI* MXX E-02 * E-02 DMI* MXX E-01 * E E-02 * E E-05.. (rest of MXX). DMI PHIX DMI* PHIX E-17 * E E E E-17 * E E E E-17 * E+00 DMI* PHIX E-01 * E E E E+00 * E E E E+00 * E-17.. (rest of PHIX). ENDDATA Анализ результатов испытаний

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Анализ результатов испытаний MATRIX COR (GINO NAME 101 ) IS A DB PREC 6 COLUMN X 6 ROW SQUARE MATRIX. COLUMN 1 ROWS 1 THRU ROW 1) D D D D D D-16 COLUMN 2 ROWS 1 THRU ROW 1) D D D D D D-09 COLUMN 3 ROWS 1 THRU ROW 1) D D D D D D-17 COLUMN 4 ROWS 1 THRU ROW 1) D D D D D D-03 COLUMN 5 ROWS 1 THRU ROW 1) D D D D D D-01 COLUMN 6 ROWS 1 THRU ROW 1) D D D D D D-15 THE NUMBER OF NON-ZERO TERMS IN THE DENSEST COLUMN = 6 THE DENSITY OF THIS MATRIX IS PERCENT. ^^^CROSS-ORTHOGONALITY MATRIX

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Выводы 1. Матрица COR указывает на хорошее совпадение форм (вне диагональные члены малы по сравнению с диагональными). 2. Опытная мода 5 - это расчетная мода 6 наоборот. На это указывают большие значения членов (5,6) и (6,5) матрицы COR. 3. Величина диагональных членов не равна 1,0 потому что опытные моды не нормализованы до единичной массы (они нормализованы до единичного максимального компонента). 4. Для подготовки к вводу опытных данных с помощью оператора DMI необходима их некоторая предварительная обработка. 5. Для анализа большого количества мод может быть разработана дополнительная процедура на языке DMAP для упрощения вида результатов расчета кросс-корреляции. 6. Процедура на языке DMAP для выполнения проверки кросс- корреляции (postmaca.vxx) – одна из поставляемых частей MSC.Nastran.

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Уточнение расчетной модели q Существует три пути улучшения расчетной модели для достижения соответствия расчетных и опытных данных: 1. Использование грубой силы 2. Использование матрицы чувствительности 3. Использование методики оптимизации q Все методы подразумевают изменение параметров модели таких как I и A для BARов и t для QUADов. Упругость основания также может рассматриваться в качестве параметра, если она в явном виде представлена в модели (например, с помощью ELASов). q Метод грубой силы: äИзменения модели производятся на базе экспертных оценок и результатов прикидочных расчетов. äПосле изменения модели расчет с помощью MSC.Nastran повторяется и его результаты анализируются на предмет лучшего согласования с опытными данными. äЕсли согласие данных не улучшается, - новые изменения и т.д.

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Уточнение расчетной модели q Матрица чувствительности: Чувствительность – это градиент отклика по параметру модели.Отклики – это результаты вычислений с помощью MSC.Nastran, апараметры – значения свойств. Каждое значение в матрице чувствительности S определяется как где R i - i-ый отклик и P j - j-ый параметр. Чем больше величина S ij, тем больше чувствительность отклика. SOL 200 в MSC.Nastran вычисляет матрицу чувствительности. Матрица S может использоваться для оценки того, какой параметр необходимо изменить для изменения соответствующего отклика.

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Уточнение расчетной модели q Матрица S также может быть получена и использована посредством метода наименьших квадратов для минимизации разницы между результатами испытаний и расчетов: P n = P o + (S T S) -1 S T (R t - R a ) где P n - новое значение параметра, подлежащего изменению P o - значение параметра до изменения S - матрица чувствительности R t - опытный отклик R a - расчетный отклик

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Уточнение расчетной модели q Использование оптимизации: äИнструмент - SOL 200 äЦель: минимизация разницы между опытными и расчетными результатами, но при минимальных изменениях в исходной модели äОператор DEQATN используется для записи целевой функции в видевзвешенной разности между опытными и расчетными результатами. где RT i -i-ый результат испытаний PF j - j-ый параметр модели после уточнения RA i - i-ый расчетный результат PO j - исходное значение j-го параметра WR i - весовой коэффициент для i-го испытания WP j - весовой коэффициент для j-го параметра wt- весовой коэффициент для всех опытных данных wp- весовой коэффициент для всех параметров qSOL 200 минимизирует целевую функцию (E) с учётом ограничений.

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Уточнение расчетной модели q Пример äМодель корпуса дисковода, включающая 1406 узлов, 1354 оболочечных элементов,четыре переменных проектирования (толщины оболочек) и четыре упругих моды, полученные в результате эксперимента äНеобходимо минимизировать разницу между частотами собственных колебаний, вычисленными с помощью MSC.Nastran и полученными экспериментально.

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Уточнение расчетной модели q Пример äВесовой коэффициент для опытных данных в 100 раз больше весового коэффициента для параметров. äУчитываются верхние и нижние ограничения на допустимую толщину оболочечных элементов. äПосле 5 итераций (шагов) достигнута сходимость äРезультаты:

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Уточнение расчетной модели Толщины оболочечных элементов Величина погрешности

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Уточнение расчетной модели qSOL 200 может также использоваться для улучшения корреляции с экспериментальными результатами по определению частотного отклика. q Возможно использование для улучшения корреляции статических результатов. q Для достижения лучших результатов необходимо использовать как можно больше опытных данных (включая вес конструкции). qВ качестве оптимизируемых необходимо рассматривать только те параметры, в значениях которых действительно есть неуверенность.

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow Литература Blakely, K. and Dobbs, M., Integrated System Identification: The Union of Testing and Analysis, Proceedings First International Modal Analysis Conference, November Kientzy, D., Richardson, M. and Blakely, K., Using Finite Element Data to Set Up Modal Tests, Sound and Vibration, June Blakely, K., Updating MSC.Nastran Models to Match Test Data, Proceedings MSC World Users Conference, March Blakely, K. and Rose, T., Cross-Orthogonality Calculations for Pre-Test Planning and Model Verification, Proceedings MSC World Users Conference, May Blakely, K., Matching Frequency Response Test Data with MSC.Nastran, Proceedings MSC World Users Conference, June 1994.

NAS102 Декабрь 2001, Стр MSC Moscow MSC Moscow