Презентация "Анализ отклика на случайное воздействие в MSC" по информатике – проект, доклад

Раздел 14

Анализ отклика на случайное воздействие

Раздел 14. Анализ отклика на случайное воздействие

ТИПЫ ДИНАМИЧЕСКИХ ПРОЦЕССОВ………………………………… 14 - 4 АНАЛИЗ ОТКЛИКА НА СЛУЧАЙНОЕ ВОЗДЕЙСТВИЕ..……...…….. 14 - 5 ТЕОРИЯ СЛУЧАЙНЫХ ПРОЦЕССОВ......……………………………… 14 - 6 ПРИМЕРЫ СЛУЧАЙНЫХ ПРОЦЕССОВ………………..……………… 14 - 7 ПРИМЕР АНСАМБЛЯ ЭРГОДИЧЕСКИХ СЛУЧАЙНЫХ ПРОЦЕССОВ……………………………………..……….. 14 - 8 АВТОКОРРЕЛЯЦИЯ И АВТОСПЕКТР………..………………………… 14 - 9 ВЫЧИСЛЕНИЕ ОТКЛИКА ЛИНЕЙНОЙ СИСТЕМЫ НА ЭРГОДИЧЕСКОЕ СЛУЧАЙНОЕ ВОЗДЕЙСТВИЕ.……………………. 14 - 11 СПЕКТРАЛЬНЫЕ СООТНОШЕНИЯ ДЛЯ ЛИНЕЙНОЙ СИСТЕМЫ С МНОГИМИ ВХОДАМИ И ВЫХОДАМИ..……………… 14 - 12 АНАЛИЗ ОТКЛИКА НА СЛУЧАЙНОЕ ВОЗДЕЙСТВИЕ В MSC.Nastran…………………………………………………..…………. 14 - 16 ОПЕРАТОР RANDPS……………………………………………………… 14 - 17 ОПЕРАТОР TABRND….………………………………………………….. 14 - 19

Анализ отклика на случайное воздействие (продолж.)

УПРАВЛЕНИЕ РЕШЕНИЕМ ПРИ АНАЛИЗЕ ОТКЛИКА НА СЛУЧАЙНОЕ ВОЗДЕЙСТВИЕ………………….. 14 - 21 РЕКОМЕНДАЦИИ ПО ВЫПОЛНЕНИЮ РАСЧЕТА..………………...… 14 - 22 ПРИМЕР №10 – ОТКЛИК ПРИ СЛУЧАЙНОМ ВОЗДЕЙСТВИИ ПО ОДНОМУ ВХОДУ………………………………..…. 14 - 23 ВХОДНОЙ ФАЙЛ ДЛЯ ПРИМЕРА №10...………………………………. 14 - 26 РЕЗУЛЬТАТЫ РЕШЕНИЯ ПРИМЕРА №10…………………………….. 14 - 27 ПРИМЕР №11 - ОТКЛИК ПРИ СЛУЧАЙНОМ ВОЗДЕЙСТВИИ ПО НЕСКОЛЬКИМ ВХОДАМ..……………………….. 14 - 35 ВХОДНОЙ ФАЙЛ ДЛЯ ПРИМЕРА №11...……………………………… 14 - 37 РЕЗУЛЬТАТЫ РЕШЕНИЯ ПРИМЕРА №11……………………………. 14 - 40

Типы динамических процессов

Случайные колебания – это колебания, которые могут быть описаны только статистически. Их истинное значение в конкретный момент времени неизвестно, чего нельзя сказать о вероятности превышения определенного уровня. Примеры: перемещение почвы при землетрясении, флюктуации давления вокруг самолетов и высоких зданий, акустическое воздействие шума ракетных и авиадвигателей. MSC.Nastran выполняет анализ отклика на случайное воздействие как постпроцессинг анализа частотного отклика. Исходной информацией являются результаты анализа частотного отклика и заданные пользователем нагрузки в форме авто- и взаимных спектральных плотностей. Результаты расчета – спектральные плотности, автокорреляционные функции, значения частоты положительных пересечений нулевого уровня и среднеквадратические отклонения откликов. Теорию см.: Random Vibration in Mechanical Systems, by S. H. Crandall and W. D. Mark, Academic Press, 1963. Дополнительная информация - в MSC.Nastran Advanced Dynamics User’s Guide.

Теория случайных процессов

Существуют различные формулировки характеристик случайных процессов. Будьте внимательны при использовании MSC.Nastran (см. MSC.Nastran Advanced Dynamics User’s Guide и книгу Бендата и Пирсола [13] для понимания соответствующих формулировок). При анализе отклика на случайное воздействие с помощью MSC.Nastran предполагается эргодичность процессов. Должна быть ясность в понятиях автокорреляции, автоспектра (спектра мощности), взаимной корреляции и кросспектра. Среднеквадратичное отклонение (значение) и кажущаяся частота – основные статистические характеристики процесса.

Примеры случайных процессов

Пример ансамбля эргодических случайных процессов

Автокорреляция и автоспектр

Автокорреляционная функция: Примечание: Rj(o) - дисперсия uj(t). Функция автоспектра (спектр, спектральная плотность): Преобразование Фурье: Дисперсия:

Кажущаяся частота (частота пересечения нулевого уровня) N0:

Вычисление отклика линейной системы на эргодическое случайное воздействие

Из анализа частотного отклика где Hja(w) -передаточная функция от входа Fa к выходу uj. При наличии нескольких входов

Спектральные соотношения для линейной системы с многими входами и выходами

В матричной форме имеем Fa (w) uj(w) = [ Hja(w) Hjb(w) … ] Fb (w) . . Автоспектр на выходе равен Fa (w) H*ja Sujuj = [ Hja Hjb … ] Fb (w) [F*a(w)F*b(w)…] H*jb . . . .

Спектральные характеристики на входе

Спектральные характеристики при многих входах и выходах:

Входная спектральная матрица Saa(w) Sab(w) … [ S ]in = Sba(w) Sbb(w) … . . . . и ее специальные свойства

Обычно используемые частные случаи Анализ с одним входом (полностью коррелированный вход) Несколько некоррелированных воздействий

Анализ отклика на случайное воздействие в MSC.Nastran

Предполагается, что результат анализа частотного отклика - Hja(w). Вычисление не производится. В случае необходимости вычисления Hja(w) задавайте F(w) = 1,0.

Управление решением при анализе отклика на случайное воздействие

Executive Control Section SOL (выбирается в соответствие с методом анализа частотного отклика) Case Control Section RANDOM (инициализирует операторы RANDPS, RANDT в Bulk Data Section и должен быть впереди операторов SUBCASE) Bulk Data Section RANDPS (задает значения спектральной плотности мощности) RANDT1 (задает временные сдвиги для вычисления автокорреляционной функции)

Рекомендации по выполнению расчета

В большинстве случаев спектр задают в логарифмических координатах. Используйте соответствующие возможности оператора TABRND1. При возможности всегда вычисляйте входную спектральную плотность (СП). Стройте графики выходной СП. Не пользуйтесь интегральными результатами вслепую. Задавайте несколько частот вблизи каждой моды. При модальном методе для этого хорошо подходят операторы FREQ1 (или FREQ2) и FREQ4. На низких частотах (

Пример №10

Отклик при случайном воздействии по одному входу

Пример №10. Отклик при случайном воздействии по одному входу

Приложите к модели пластины случайное воздействие в направлении Z с заданной СП. Левое ребро соедините с помощью элемента RBE2 с узлом 9999 и приложите вынужденное перемещение к этому узлу Используйте модальное решение с остаточным вектором

Задайте постоянное модальное демпфирование на уровне 0,03 от критического. Используйте задание СП на входе в двойных логарифмических координатах. Определите СП ускорения в точке воздействия, а также в угле и в середине свободного ребра (узлы 33 и 55)

Входной файл для Примера №10

ID SEMINAR, PROB10 SOL 111 CEND TITLE= RANDOM ANALYSIS - BASE EXCITATION SUBTITLE= USING THE MODAL METHOD WITH LANCZOS ECHO= UNSORTED SPC= 101 SET 111= 33, 55, 9999 ACCELERATION(SORT2, PHASE)= 111 METHOD= 100 FREQUENCY= 100 SDAMPING= 100 RANDOM= 100 DLOAD= 100 $ OUTPUT(XYPLOT) XTGRID= YES YTGRID= YES XBGRID= YES YBGRID= YES YTLOG= YES XTITLE= FREQUENCY YTTITLE= ACCEL RESPONSE BASE, MAGNITUDE YBTITLE= ACCEL RESPONSE AT BASE, PHASE XYPLOT ACCEL RESPONSE / 9999 (T3RM, T3IP) YTTITLE= ACCEL RESPONSE AT TIP CENTER, MAGNITUDE YBTITLE= ACCEL RESPONSE AT TIP CENTER, PHASE XYPLOT ACCEL RESPONSE / 33 (T3RM, T3IP) YTTITLE= ACCEL RESPONSE AT OPPOSITE CORNER, MAGNITUDE YBTITLE= ACCEL RESPONSE AT OPPOSETE CORNER, PHASE XYPLOT ACCEL RESPONSE / 55 (T3RM, T3IP) $ $ PLOT OUTPUT IS ONLY MEANS OF VIEWING PSD DATA $ XGRID= YES YGRID= YES XLOG= YES YLOG= YES YTITLE= ACCEL P S D AT LOADED CORNER XYPLOT ACCEL PSDF / 9999(T3) YTITLE= ACCEL P S D AT TIP CENTER XYPLOT ACCEL PSDF / 33(T3) YTITLE= ACCEL P S D AT OPPOSITE CORNER XYPLOT ACCEL PSDF / 55(T3) $

BEGIN BULK param,resvec,yes PARAM,COUPMASS,1 PARAM,WTMASS,0.00259 $ INCLUDE 'plate.bdf' $ GRID, 9999, , 0., 1., 0. $ RBE2, 101, 9999, 12345, 1, 12, 23, 34, 45 $ SPC1, 101, 12456, 9999 $ $ EIGENVALUE EXTRACTION PARAMETERS $ EIGRL, 100 , , 2000. $ $ SPECIFY MODAL DAMPING $ TABDMP1, 100, CRIT, +, 0., .03, 10., .03, ENDT $ $ POINT LOADING AT TIP CENTER $ RLOAD2, 100, 600, , , 310,,A spcd,600,9999,3,1.0 spc1,101,3,9999 $ TABLED1, 310, +, 10., 1., 1000., 1., ENDT $ $ SPECIFY FREQUENCY STEPS $ FREQ,100,30. FREQ1,100,20.,20.,50 FREQ4,100,20.,1000.,.03,5 $ $ SPECIFY SPECTRAL DENSITY $ RANDPS, 100, 1, 1, 1., 0., 111 $ TABRND1, 111,LOG,LOG +, 20., 0.1, 30., 1., 100., 1., 500., .1, +, 1000., .1, ENDT $ ENDDATA

Результаты решения Примера №10

Пример №11

Отклик при случайном воздействии по нескольким входам

Пример №11. Отклик при случайном воздействии по нескольким входам

Используя модальный метод, определите спектр перемещения в центральной точке пластины под действием случайного давления и сосредоточенных сил. Используйте комплексное представление взаимной спектральной плотности.

Входной файл для Примера №11

ID SEMINAR, PROB11 SOL 111 TIME 30 CEND TITLE= FREQUENCY RESPONSE WITH PRESSURE AND POINT LOADS SUBTITLE= USING THE MODAL METHOD WITH LANCZOS ECHO= UNSORTED SPC= 1 SET 111= 11, 33, 55 DISPLACEMENT(PLOT, PHASE)= 111 METHOD= 100 FREQUENCY= 100 SDAMPING= 100 RANDOM= 100 SUBCASE 1 LABEL= PRESSURE LOAD DLOAD= 100 LOADSET= 100 SUBCASE 2 LABEL CORNER LOAD DLOAD= 200 LOADSET= 100 $

$ OUTPUT (XYPLOT) $ XTGRID= YES YTGRID= YES XBGRID= YES YBGRID= YES YTLOG= YES YBLOG= NO XTITLE= FREQUENCY (HZ) YTTITLE= DISPLACEMENT RESPONSE AT LOADED CORNER, MAGNITUDE YBTITLE= DISPLACEMENT RESPONSE AT LOADED CORNER, PHASE XYPLOT DISP RESPONSE / 11 (T3RM, T3IP) YTTITLE= DISPLACEMENT RESPONSE AT TIP CENTER, MAGNITUDE YBTITLE= DISPLACEMENT RESPONSE AT TIP CENTER, PHASE XYPLOT DISP RESPONSE / 33 (T3RM, T3IP) YTTITLE= DISPLACEMENT RESPONSE AT OPPOSITE CORNER, MAGNITUDE YBTITLE= DISPLACEMENT RESPONSE AT OPPOSITE CORNER, PHASE XYPLOT DISP RESPONSE / 55 (T3RM, T3IP) $ $ PLOT OUTPUT IS ONLY MEANS OF VIEWING PSD DATA $ XGRID= YES YGRID= YES XLOG= YES YLOG= YES YTITLE= DISP P S D AT LOADED CORNER XYPLOT DISP PSDF / 11(T3) YTITLE= DISP P S D AT TIP CENTER XYPLOT DISP PSDF / 33(T3) YTITLE= DISP P S D AT OPPOSITE CORNER XYPLOT DISP PSDF / 55(T3)

BEGIN BULK PARAM,COUPMASS,1 PARAM,WTMASS,0.00259 $ $ MODEL DESCRIBED IN NORMAL MODES EXAMPLE $ INCLUDE ’plate.bdf’ $ $ EIGENVALUE EXTRACTION PARAMETERS $ EIGRL, 100, 10., 2000. $ $ SPECIFY MODAL DAMPING $ TABDMP1, 100, CRIT, +, 0., .03, 10., .03, ENDT $ $ FIRST LOADING $ RLOAD2, 100, 300, , , 310 $ TABLED1, 310, +, 10., 1., 1000., 1., ENDT $ $ UNIT PRESSURE LOAD TO PLATE $ LSEQ, 100, 300, 400 $ PLOAD2, 400, 1., 1, THRU, 40 $ $ SECOND LOADING $ RLOAD2, 200, 600, , , 310 $ $ POINT LOAD AT TIP CENTER $ DAREA, 600, 11, 3, 1. $

$ SPECIFY FREQUENCY STEPS $ FREQ1, 100, 20., 20., 49 $ $ SPECIFY SPECTRAL DENSITY $ RANDPS, 100, 1, 1, 1., 0., 100 RANDPS, 100, 2, 2, 1., 0., 200 RANDPS, 100, 1, 2, 1., 0., 300 RANDPS, 100, 1, 2, 0., 1., 400 $ TABRND1, 100, +, 20., 0.1, 30., 1., 100., 1., 500., .1, +, 1000., .1, ENDT $ TABRND1, 200, +, 20., 0.5, 30., 2.5, 500., 2.5, 1000., 0., +, ENDT $ TABRND1, 300, +, 20., -.099619, 100., -.498097, 500., .070711, 1000., 0., +, ENDT $ TABRND1, 400, +, 20., .0078158, 100., .0435791, 500., -.70711, 1000., 0., +, ENDT $ ENDDATA

Результаты решения Примера №11

R E A L E I G E N V A L U E S MODE EXTRACTION EIGENVALUE RADIANS CYCLES GENERALIZED GENERALIZED NO. ORDER MASS STIFFNESS 1 1 7.056994E+05 8.400591E+02 1.336996E+02 1.000000E+00 7.056994E+05 2 2 1.878432E+07 4.334088E+03 6.897916E+02 1.000000E+00 1.878432E+07 3 3 2.811467E+07 5.302327E+03 8.438915E+02 1.000000E+00 2.811467E+07 1 FREQUENCY RESPONSE WITH PRESSURE AND POINT LOADS APRIL 8, 1998 MSC.Nastran 4/ 6/98 PAGE 13 USING THE MODAL METHOD WITH LANCZOS 0 *** USER INFORMATION MESSAGE 5222 (SQFREQ) UNCOUPLED SOLUTION ALGORITHM USED. 1 FREQUENCY RESPONSE WITH PRESSURE AND POINT LOADS APRIL 8, 1998 MSC.Nastran 4/ 6/98 PAGE 14 USING THE MODAL METHOD WITH LANCZOS 0 1 FREQUENCY RESPONSE WITH PRESSURE AND POINT LOADS APRIL 8, 1998 MSC.Nastran 4/ 6/98 PAGE 15 USING THE MODAL METHOD WITH LANCZOS 0 PRESSURE LOAD SUBCASE 1 0 X Y - O U T P U T S U M M A R Y ( R E S P O N S E ) 0 SUBCASE CURVE FRAME XMIN-FRAME/ XMAX-FRAME/ YMIN-FRAME/ X FOR YMAX-FRAME/ X FOR ID TYPE NO. CURVE ID. ALL DATA ALL DATA ALL DATA YMIN ALL DATA YMAX 0 1 DISP 1 11( 5,--) 2.000000E+01 1.000000E+03 4.931745E-04 1.000000E+03 2.648644E-01 1.400000E+02 2.000000E+01 1.000000E+03 4.931745E-04 1.000000E+03 2.648644E-01 1.400000E+02 0 1 DISP 1 11(--, 11) 2.000000E+01 1.000000E+03 2.067484E+01 1.000000E+03 3.594679E+02 2.000000E+01 2.000000E+01 1.000000E+03 2.067484E+01 1.000000E+03 3.594679E+02 2.000000E+01 0 2 DISP 2 11( 5,--) 2.000000E+01 1.000000E+03 7.310375E-05 4.000000E+02 6.715895E-02 1.400000E+02 2.000000E+01 1.000000E+03 7.310375E-05 4.000000E+02 6.715895E-02 1.400000E+02 0 2 DISP 2 11(--, 11) 2.000000E+01 1.000000E+03 1.839137E+02 2.600000E+02 3.595085E+02 2.000000E+01 2.000000E+01 1.000000E+03 1.839137E+02 2.600000E+02 3.595085E+02 2.000000E+01 0 1 DISP 3 33( 5,--) 2.000000E+01 1.000000E+03 5.224569E-04 1.000000E+03 2.652415E-01 1.400000E+02 2.000000E+01 1.000000E+03 5.224569E-04 1.000000E+03 2.652415E-01 1.400000E+02 0 1 DISP 3 33(--, 11) 2.000000E+01 1.000000E+03 2.008288E+01 1.000000E+03 3.594677E+02 2.000000E+01 2.000000E+01 1.000000E+03 2.008288E+01 1.000000E+03 3.594677E+02 2.000000E+01 0 2 DISP 4 33( 5,--) 2.000000E+01 1.000000E+03 4.193472E-05 6.000000E+02 6.768133E-02 1.400000E+02 2.000000E+01 1.000000E+03 4.193472E-05 6.000000E+02 6.768133E-02 1.400000E+02 0 2 DISP 4 33(--, 11) 2.000000E+01 1.000000E+03 1.822376E+02 3.400000E+02 3.594847E+02 2.000000E+01 2.000000E+01 1.000000E+03 1.822376E+02 3.400000E+02 3.594847E+02 2.000000E+01 0 1 DISP 5 55( 5,--) 2.000000E+01 1.000000E+03 4.931745E-04 1.000000E+03 2.648644E-01 1.400000E+02 2.000000E+01 1.000000E+03 4.931745E-04 1.000000E+03 2.648644E-01 1.400000E+02 0 1 DISP 5 55(--, 11) 2.000000E+01 1.000000E+03 2.067484E+01 1.000000E+03 3.594679E+02 2.000000E+01 2.000000E+01 1.000000E+03 2.067484E+01 1.000000E+03 3.594679E+02 2.000000E+01 0 2 DISP 6 55( 5,--) 2.000000E+01 1.000000E+03 1.551823E-04 1.000000E+03 6.802692E-02 1.400000E+02 2.000000E+01 1.000000E+03 1.551823E-04 1.000000E+03 6.802692E-02 1.400000E+02 0 2 DISP 6 55(--, 11) 2.000000E+01 1.000000E+03 7.881378E+00 7.599999E+02 3.594572E+02 2.000000E+01 2.000000E+01 1.000000E+03 7.881378E+00 7.599999E+02 3.594572E+02 2.000000E+01

1 FREQUENCY RESPONSE WITH PRESSURE AND POINT LOADS APRIL 8, 1998 MSC.Nastran 4/ 6/98 PAGE 16 USING THE MODAL METHOD WITH LANCZOS 0 0 X Y - O U T P U T S U M M A R Y ( A U T O O R P S D F ) 0 PLOT CURVE FRAME RMS NO. POSITIVE XMIN FOR XMAX FOR YMIN FOR X FOR YMAX FOR X FOR* TYPE TYPE NO. CURVE ID. VALUE CROSSINGS ALL DATA ALL DATA ALL DATA YMIN ALL DATA YMAX 0 PSDF DISP 7 11( 5) 1.377514E+00 1.401104E+02 2.000E+01 1.000E+03 2.432E-08 1.000E+03 5.941E-02 1.400E+02 0 PSDF DISP 8 33( 5) 1.377778E+00 1.360495E+02 2.000E+01 1.000E+03 2.730E-08 1.000E+03 5.963E-02 1.400E+02 0 PSDF DISP 9 55( 5) 1.379390E+00 1.432761E+02 2.000E+01 1.000E+03 2.432E-08 1.000E+03 5.951E-02 1.400E+02