Применение инновационной техники Structural Intensity Analysis для повышения вибрационной надежности и усовершенствования акустического дизайна рулевой колонки автомобиля Текст научной статьи по специальности «Физика»

УДК 623113

Н.А. Кузьмин, В.А. Шапкин, Ю.В. Шапкина

ПРИМЕНЕНИЕ ИННОВАЦИОННОМ ТЕХНИКИ STRUCTURAL INTENSITY ANALYSIS ДЛЯ ПОВЫШЕНИЯ ВИБРАЦИОННОЙ НАДЕЖНОСТИ И УСОВЕРШЕНСТВОВАНИЯ АКУСТИЧЕСКОГО ДИЗАЙНА РУЛЕВОЙ КОЛОНКИ АВТОМОБИЛЯ

Нижегородский государственный технический университет им. Р.Е. Алексеева

В данной статье представлена система численного моделирования и визуализации Structural Intensity Analysis, базирующаяся на методе конечно элементных моделей (FEM-models). Приведена математическая формулировка структурной интенсивности, система численного моделирования осуществлена в программном пакете ANSYS, техника визуализации векторных полей потока энергии осуществлена в программном пакете HyperView. Рассчитаны и показаны в структуре shell-элементов: векторы интенсивности, карта потоков и распределение механической энергии колебаний. Рассмотрены случаи использования различных возбудителей гармонических колебаний (точечной силы и давления на несколько элементов) и применения различных коэффициентов поглощения колебаний структурой. Используя разработанную систему, был осуществлен Structural Intensity Analysis рулевой колонки автомобиля Volkswagen, визуализованы поля структурной интенсивности на глобальной поверхности и исследованы доминантные пути потока энергии. Structural Intensity Analysis, вместе с техникой его визуализации, представляет собой мощный инструмент для контроля над вибрацией и шумом.

Ключевые слова: Structural Intensity Analysis, рулевая колонка автомобиля, виброакустическое поведение, FEM-моделирование, шум и вибрация, надежность систем автомобиля.

Введение

В связи с все более ужесточающимися требованиями к качеству автомобилей и возрастающими ожиданиями клиентов в отношении комфорта, исследование и контроль вибрации и акустического поведения продуктов автомобильной промышленности стало важным критерием в процессе их разработки. Для дальнейшего усовершенствования в конструкции транспортного средства и его составляющих, с учетом акустических характеристик, а также повышения надежности автомобиля при воздействии вибрационных нагрузок, необходимо использовать инновационные технологии на этапе его расчета и FEM -анализ модели [1].

Рис. 1. Активная Структурная Интенсивность при второй резонансной частоте. Масленый поддон двигателя, Дармштадский технический университет, 2012

© Кузьмин Н.А., Шапкин В.А., Шапкина Ю.В., 2013

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

Так как виброакустическое поведение продукта в основном определяется поведением вовлеченных механических структур, мы должны сосредоточиться на численном прогнозировании этого поведения. Реализация разработок в вибрационном поведении автомобиля и его частей требует детального знания динамических процессов возбуждения и излучения звука в структуре. Существует целый ряд численных методов, используемых для структурно-динамического анализа (Structural Dynamic Analysis) автомобильных структур экспериментально или в технике FEM-моделирования. Но ни один из них не может снабдить нас достаточно точной информацией о волновом потоке в структуре, который производит шум и вибрацию.

Structural Intensity Analysis (Анализ Структурной Интенсивности) - это новый шаг в анализе виброакустического поведения моделей, наиболее эффективная методика определения и контроля вибрации. Данный уникальный метод позволяет детально рассмотреть поток от точки возбуждения вибрации в структуре до точки ее поглощения. Анализ Структурной Интенсивности позволяет рассчитать поток колебательной энергии в пределах структуры, основываясь на методе конечных элементов, и визуализировать его в векторном виде. Составляющие потока энергии могут быть разбиты на отдельные компоненты вибрации, такие как поперечные и продольные волны [2].

а) б)

Рис. 2. Активная Структурная Интенсивность на модели с закреплением некоторых элементов (а) и свободной (б). Пол автомобиля, BMW-групп и Дармштадский технический университет, AAC конференция 2011

Знаний, полученные при использовании Structural Intensity Analysis, могут быть использованы для снижения или перенаправления потока энергии в целевых районах с высокой интенсивностью для обеспечения оптимального вибрационного поведения продукции.

Теоретические основы Структурной Интенсивности

Термин Структурная Интенсивность и ее математическое определение были введены впервые Нойзеуксом [3] и Павиком [4] в 70-х годах прошлого столетияс целью расширения описания вибрационных явлений в структуре с применением векторных полей.

Структурная Интенсивность (Structure Intensity) - это усредненный во времени поток чистой механической энергии через единичную площадку, перпендикулярную направлению потока. Интерес к исследованию Структурной Интенсивности имеет под собой практическое основание - поля интенсивности показывают величину и направление потока колебательной

энергии в любой точке структуры. Техника вибрационной интенсивности не зависит от граничных условий структуры и дает возможность исследовать краевой эффект от передачи вибрационной мощности по структуре [5, 6]. Поля Структурной Интенсивности могут быть визуализированы с помощью векторной карты и предоставляют точную информацию о потоке энергии, доминирующих путях его передачи, положении источников и поглотителей механической энергии. Применение данной техники снабжает нас преимуществом при изучении влияния от добавления диссипативных элементов, изменения механического состава или при активном контроле вибрации.

Структурная Интенсивность выражается как произведение напряжения на поверхностную скорость, и во временной доминанте будет иметь вид:

I =(!г )) = (-*у (*>у (*)), (1)

4" хх^х + Хху^у + Х хЛ

!у =- Х ух^х + yyVy + X у^ , (2)

к ] _ХЛ + Хгу*у + ^ _

где I - это интенсивность в направлении /; V - скорость частиц в направлении /; стг7 (^) -

компонент нормального напряжения.

Колебательная энергия в структурах передается при помощи поперечных сил и Qy), сил в плоскости (Мх, Ыу и Ыху), крутящих моментов (Мху и Мух). Таким образом, компонент вибрационной интенсивности 1х в направлении х записывается в виде:

, д% ,, д2Т лг д£, ЛГ 5л

4 = Ох—- мх——+мху——+ых—+ыху—

х х дл дхд1 ху дудг х дг у дг

(3)

Эти силы и моменты связаны с поперечным смещением следующим образом:

О =- в

гд% д3Т Л +

дх3 дхду

2

Мх = -£

Гд\ д2^ -Т+Цдт

)

N =■

ЕИ

дТ дл . —+ ц—1 1 - ц2 ^ дх дуу

\

ЕИ

дх'

Мху = - £>(1 -ц)

ду2

N = Е^ дТ +дЛ ху 2(1 + ц) ^ду дх)

дхду'

С

\

в = ■

ЕИ3

12(1 -ц2)

Формулировка Структурной Интенсивности для пластины

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

Рис. 3. Элемент пластины с определенными силами и перемещениями

Так же из данной теории следует, что все структурные точки нормали к средней плоскости остаются на прямой после деформации. 1Х в (2) может быть выражена в глобальной системе координат следующим образом:

Ix =-

i -и2

fe , ,,3Ф* , ,Лм .....^z, fe , ,x , dU , ,z 1

+у—+ y—

сХ dx dz dz

x

V

Ф x V

dx

• Vy + G •

?M ^ | ^м ^.^z dz dz dx dx

V, (4)

(4)

где E - модуль упругости Юнга; ц - коэффициент Пуассона; G = E / (1 - 2ц) - модуль сдвига; фх и ф - вращательные перемещения в направлениях х и у, vx , vy и vz - скорости в х, у, z

направлениях. На рис. 1 показана положительная ориентация внутренних сил и перемещений.

Гаврич и Павик [6] разработали эквивалентные алгоритмы структурной интенсивности вибрации, которые могут быть реализованы во временной и частотной доменах. Формулировки структурной интенсивности по х и у для плоской пластины в частотной домене могут быть выражены следующим образом:

Ix = -(ш / 2) Im[/~x~ * + КуГ+ QxW + Mx 9у* + Mxy~;l (5)

Iy = -(ш / 2)1т[/~у~* + Nyx~* + Q,W* + My Qy* + Myx0y*\, (6)

где Nx, Ny и N^ = N комплексные мембранные силы на единицу ширины пластины; Mx My и M = M- комплексные изгибающие и крутящие моменты на единицу ширины пластины; Q и Q - комплексные поперечные силы сдвига на единицу ширины пластины; u , V и w комплексно сопряженные поступательные перемещения в направлениях х, у и z; 9 * и 9y комплексно сопряженные вращательные перемещения в направлениях х и у.

Система моделирования Structural Intensity Analysis

Для эффективного анализа полей структурной интенсивности и их визуализации была разработана система численного моделирования. Система состоит из программы ANSYS для FEM-моделирования и процедуры анализа, и программы HyperView для визуализации полей структурной интенсивности на модели.

С использованием данной системы был произведен Анализ Структурной Интенсивности простой пластины, для верификации процедуры. Затем был произведен анализ рулевой колонки автомобиля, испытывающей вибрационные нагрузки, предающиеся на нее от подвески автомобиля.

Результаты и визуализация Structural Intensity Analysis для конечно-элементной модели пластины

Для пластины были взяты геометрические размеры, граничные условия, частоты возбуждения, свойства материала и условия нагружения, соответствующие параметрам в работе S. Kuhl «Strukturintensitätsanalyse als Werkzeug der Maschinenakustik», Дармштадский технический университет [7], затем была осуществлена симуляция по разработанной методе. Сравнение Анализа Структурной Интенсивности с предыдущими результатами подтвердило правильность использованного алгоритма.

Пластина смоделирована с помощью стандартных shell-элементов. Затухание представлено как коэффициент демпфирования. Пластина закреплена по краям. Пластина находится под гармоническим возбуждением. Результаты моделирования потока интенсивности в векторной форме показаны на рис. 4.

Рис. 4. Сравнение результатов данного исследования (а, в, д) с результатам из работы S. Kuhl "Gezielte Leitung von Körperschall unter Zuhilfenahme der Strukturintensitätsrechnung" (б, г, е). Векторные поля активной (а, б ) и реактивной (в, г) Структурной Интенсивности, форма изгиба (д, е) пластины при третьей натуральной частоте

Точка приложения возбуждающей силы выступает в качестве источника энергии, точка затухания колебаний - в качестве поглотителя энергии, и основной поток энергии проходит путь от точки приложения силы до точки демпфирования.

Рис. 5. Распределение структурной интенсивности при возбуждении колебаний в пластине при помощи точечной силы (а) и давления (б). Коэффициент демпфирования 0,6

Из представленных результатов видно, что поток энергии от точки возбуждения (источника энергии) к точке демпфирования (концентратора энергии) может быть четко идентифицирован векторами Структурной Интенсивности.

Было произведено сравнение картин распределения Интенсивностей, вызванных воздействием точечной силы и давлением на несколько элементов пластины, как это представлено на рис. 5. В результате было установлено, что форма векторной картины Структурной Интенсивности не зависит от вида возбуждающей силы, при одинаковом расположении точки различие возбуждения выражается только в величине векторов.

Изучены случаи применения различных демпфирующих коэффициентов. Визуализация результирующих векторных потоков интенсивности энергии приведена в табл. 1. При использовании более высокого коэффициента затухания наблюдается потеря определенных точек демпфирования и поглощение волн всей структурой равномерно.

Таблица 1

Сравнение картин векторных полей Структурной Интенсивности в пластине при различных коэффициентах демпфирования

0,1

0,3

0,5

0,6

, Time/Freq 305 000000

Contour Plot Actlve(Z, Max) Global System Advanced Average - I 077E-03 B—9.576E-04 —8.380E-04 —7.183E-04

-5.986E-04

M— L

Max = 1.077E-03 №de 462 Min = 9 145E-08

Substep 1, Time/Freq 305.0000

/ —\ \ if f , ^ . ~ , , у____. \ ^ 1 I t / ■

ii/ , I I f

I i i v . , t t t t > It , -\ \ \ ' 'tt> f^i ■

--""'л '

• - - "^T7/

в

г)

Рис. 6. Структурная интенсивность (а) на пластине; активная (б); реактивная (в);

напряжения по фон Мизу (г)

Structural Intensity Analysis для рулевой колонки автомобиля

Техника Структурной Интенсивности успешно реализована в области вибрации элементарной пластины для определения основных путей передачи энергии. Однако на данный момент существует очень мало исследований с использованием подхода структурной интенсивности для оценки вибрации реальных структур.

В данной работе проведено численное моделирование системы для анализа и визуализации структурной интенсивности сложных структур, таких как структура частей автомобиля. Анализ Структурной Интенсивности для рулевой колонки и ее части был проведен в их разработанных проектных состояниях. Целью являлось исследование характеристик передачи энергии колебаний для определения доминирующего пути потока энергии.

Приведена визуализации результатов воздействия силы и момента возбуждения на основную структурную деталь рулевой колонки - шпиндель, данная деталь является основным передатчиком вибрации на рулевую колонку. По результатам анализа, оценены доминирующие пути потока вибрационной энергии от источника возбуждения к деталям рулевой колонки.

Деталь шпиндель закреплена с одного конца и к другому концу приложена возбуждающая гармоническая сила. Результаты картины структурной интенсивности наблюдались при второй натуральной частоте. Из рис. 7 видно, что поля распространения вибрационной энергии, визуализированные при помощи техники структурной интенсивность, не возможно наблюдать при других видах анализа - сравнение с нижними картинами напряжений и деформации детали.

Рис. 7. Шпиндель рулевой колонки. Активная Структурная Интенсивность (а), Реактивная Структурная Интенсивность (б), напряжение в направлении оси z (в), перемещения (г).

Коэффициент затухание 0,0001, /=193 Гц

Процедура анализа проведена для различных значений поглощения колебаний структурой. На рис. 8 можно наблюдать различие в картинах распространения Активной Структурной Интенсивности (меры переносимой энергии) при различных коэффициентах затухания.

а)

б)

Рис. 8. Поля Активной Структурной Интенсивности на шпинделе рулевой колонки при различных коэффициентах затухания 0.0001 (а) и 0.1 (б) колебаний в материале, /=193 Гц

В работе проведен Structural Intensity Analysis рулевой колонки автомобиля Volkswagen. Смоделированы условия закрепления модели и произведено возбуждение колебаний при помощи гармонической силы. На рис. 9 и рис. 10 показаны картины распространения полей Активной Структурной Интенсивности при первой и второй натуральных частотах. Как известно, возбуждающие колебания, совпадающие с собственными колебаниями структуры, являются наиболее опасными с точки зрения ее разрушения, поэтому исследование распространения вибрационной энергии в данных случаях является наиболее важным.

Contour Plot VW PQ35SP, fullout/down, 04 07.2012 : Substep 3, Tlme/Freq 59.825268

г '-'Г '' "-i-i i energy Mid) ,-„.., В ' Angle 315.00GG00

Simple Average

Рис. 9. Распространение вибрационной энергии и Активная Структурная Интенсивность при первой натуральной частоте. Рулевая колонка Volkswagen

Рис. 10. Распространение вибрационной энергии и Активная Структурная Интенсивность

при второй натуральной частоте. Рулевая колонка Volkswagen

Заключение

В данном исследовании разработана система численного моделирования для Structural Intensity Analysis и его визуализации, приведены его теоретические основы и осуществлено моделирование на реальных компонентах. Из Structural Intensity Analysis рулевой колонки Volkswagen удостоверено, что система моделирования может эффективно применяться для структурного анализа интенсивности и визуализации результатов на глобальной структуре автомобиля или его составных частей. Более того, установлено, что доминирующим путем передачи вибрационной энергии от источников возбуждения является шпиндель рулевой колонки, он же является деталью, воспринимающей наибольшие вибрационные нагрузки. Получена детальная информация о путях передачи вибрационной энергии по структуре рулевой колонки и установлены наиболее уязвимые с акустической точки зрения места. Данная информация не может быть получена через стандартные методы анализ структуры под нагрузкой.

Structural Intensity Analysis может быть эффективно применен для принятия надлежащих контрмер против проблемной вибрации и шума в конструкции рулевой колонки. Этот метод позволяет проводить детальное исследование в районах высокой концентрации вибрационного потока энергии, которые, следовательно, подвергаются риску повреждения или распространяют звуковые волны в окружающую среду. Он может также рассматриваться как идентификатор регионов для применения дополнительного затухания в целях снижения уровня вибрации и шума. Было показано, что Structural Intensity Analysis может выступать в качестве нового критерия надежности конструкции.

Библиографический список

1. Stoewer, T. Efficiency improvements in the acoustical design of vehicle structures by means of energy flow calculation / T. Stoewer, H. Schneeweiß [BMW Group], T. Hering, Dr.-Ing. // Bös , Prof. DrIng. H. Hanselka [TU Darmstadt]: AAC conference, 2011.

2. Hering ,T. Strukturintensitätsanalyseals Werkzeug der Maschinenakustik / T. Hering. - TU Darmstadt, 2012.

3. Noiseux, D.U. Measurement of power flow in uniform beams and plates // Acoustical Society of America. 1970. 47. Р. 238-247.

4. Pavic, G. Measurement of structure borne wave intensity - Part I: formulation of the methods // Journal of Sound and Vibration. 1976. 49(2). Р. 221-230.

5. Gavric, L. A finite element method for computation of vibration intensity by the normal mode approach / L. Gavric, G. Pavic // Journal of Sound and Vibration. Vol. 164. No. 1. Р. 29-43 (1993).

6. Hambric, S. A. Comparison of experimental and finite element structure-borne flexural wave power measurements for straight beam / S. Hambric, A. P. D. Taylor // Journal of Sound and Vibration. 1994. Vol. 175. No. 5. Р. 595-605.

7. Kuhl, S. Gezielte Leitung von Körperschall unter Zuhilfenahme der Strukturintensitätsrechnung / S. Kuhl. - TU Darmstadt, 2010.

8. Кузьмин, Н.А. Процессы и причины изменения работоспособности автомобилей: учеб. пособие / Н.А Кузьмин; НГТУ. - Н. Новгород, 2005. - 160 с.

9. Шапкина, Ю.В. Влияние конструктивных параметров большегрузных автомобилей на характеристики управляемости / Ю.В. Шапкина, В.Б. Платонов // Проблемы транспортных и технологических комплексов: сб. научн. ст. III Международной научно-технической конференции, посвященной 40-летию кафедры «Строительные и дорожные машины»; НГТУ. - Н. Новгород, 2012. С. 129-131.

Дата поступления в редакцию 05.02.2013

N. A. Kuzmin, J. V. Shapkina, V.A. Shapkin

THE STRUCTURAL INTENSITY ANALYSIS FOR PREDICTION OF VIBRATION ENERGY FLOW AND IMPROVEMENTS IN THE ACOUSTICAL DESIGN OF STEERING COLUMN PART

Nizhny Novgorod state technical university n.a. R.E. Alexeev

Purpose: to develop the new system for Structural Intensity Analysis and visualization to apply for shell-element automobiles structures based on the FEM- method.

Design/methodology/approach: The mathematical formulation of the Structural Intensity has introduced, the numerical simulation system has implemented in ANSYS program software, the visualization technique of the vector fields of energy flow has realized in HyperView software.

Findings: The intensity vectors, streamline map and energy distribution in shell-elements plate under excitation force are calculated and visualized to predict mechanical energy transmission. The energy sources, sinks and transmission paths are clearly shown with visualization technique. Using the system, the Structural Intensity Analysis for the part of Steering Column was carried out to visualize structural intensity fields on the global shell-elements structure. Research limitations/implications: The present study provides the new powerful tool for vibration and noise control. The application of the Structural Intensity technique together with visualization methods has improved the quality of structure-borne noise diagnostics. This method enables the investigation in the regions of high concentration of vibration energy flow which consequently is exposed to the risk of damage or is propagating the sound waves to the environment. It can be also considered as the identification of the regions for application of additional damping in purpose of lowering of vibration level and resulting noise radiation. It has been shown that the Structural Intensity Analysis can act as new criteria of reliability for structural design.

Originality/value: For the present moment it is quite innovative study because the structural intensity approach had been implemented for noise and vibration of real automobile structures.

Key words: Structural Intensity Analysis, Steering Column, vibro-acoustic behavior, FEM-modeling, noise and vibration.