Расчетно-экспериментальная методика достижения нормативного каскадного перепада характеристик вибрации амортизированного судового оборудования

ПРОБЛЕМЫ КОРАБЛЕСТРОЕНИЯ И ОКЕАНОТЕХНИКИ

УДК 624. 042:629.5.024

1 2 В.И. Сутырин , Е.И. Короткая

РАСЧЕТНО-ЭКСПЕРИМЕНТАЛЬНАЯ МЕТОДИКА ДОСТИЖЕНИЯ НОРМАТИВНОГО КАСКАДНОГО ПЕРЕПАДА ХАРАКТЕРИСТИК ВИБРАЦИИ АМОРТИЗИРОВАННОГО СУДОВОГО ОБОРУДОВАНИЯ

Балтийский федеральный университет им. И. Канта, г. Калининград1, Балтийская государственная академия рыбопромыслового флота, г. Калининград2

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

Ключевые слова: Анализ колебаний, корпусные конструкции, метод конечных элементов, квазистатическая модель, прием конденсации.

Судовое оборудование и механизмы (двигатели, генераторы, и др.) крепятся с целью снижения вибрации к корпусу судна через амортизирующие элементы-амортизаторы различных конструкций (рис. 1). В случае применения двухкаскадной амортизации амортизируемый объект 2 через амортизаторы первого каскада опирается на промежуточную опорную раму 3, которая, в свою очередь, через амортизаторы второго каскада - на палубу 1, платформу или днище судна.

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

Поскольку наибольший интерес представляют, как правило, вертикальные колебания системы, решение задачи в первом приближении обычно строится на основе трехмассовой модели (рис. 2), объединяющей М2 - массу опорной платформы; М1 - массу оборудования; С1, С2 - каскадные жесткости (их значения могут оказаться неизвестными); М3 и С3 - массу и жесткость, приведенные к месту закрепления оборудования (А) со стороны корпуса. Учитывается, что первая масса находится под действием гармонической вынуждающей силы с амплитудой Г и частотой изменения ш. Для максимального упрощения математического описания колебательного процесса в настоящей статье рассматривается вариант опирания обо-

© Сутырин В.И., Короткая Е.И., 2013.

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

Рис. 1. Двухкаскадная амортизация

Рис. 2. Трехмассовая модель

При подсоединении к оборудованию трубопровода 5 в расчет вводятся дополнительные обобщенные коэффициенты массы МТ и жесткости СТ первого каскада, получаемые путем предварительного приведения массы и жесткости [1] трубопровода к точке его соединения с оборудованием (В) в соответствии с равенствами [2, 3]:

СТ rB \

-ко}-II^qq|| -{Rcb}; Mt = mB +{GlG}T;

{G} = -{Lf -1M оо 1Г2; \L}= Hl RanV - ROb },

где гв, тв - коэффициенты жесткости и массы точки В в исходных матрицах трубопровода;

и |Млп1! - блок-матрицы коэффициентов жесткости и массы исключаемых узлов (об-

т

ласти модели О) трубопровода; {^ВО}и {Япв}= [ЯВо] - векторы коэффициентов жесткости взаимного влияния узлов.

На первом этапе решения задачи определяются, если это необходимо, жесткости каскадов системы существующей амортизации. Система разрешающих уравнений имеет вид:

где X={ Х1,Х2 ,Х3} - вектор амплитуд колебаний масс.

Переход к алгебраическому виду и последующие преобразования позволяют получить зависимости для определения каскадных жесткостей С1 и С2:

С1 ^ Ст - С1 0 М1 + Мт 0 0 X1 f

С1 + С2 - С2 -ш2 М 2 X2 = 0

симм. С2 + С3 симм . м 3 X3 0

С = {f + X! - [ш2(М1 + Mt ) - Ст ]}• X - X2 )-1;

С2 =[х2 • (ш2 • М2 -(Х2 -Х3)-1.

Для определения числовых значений С1 и С2 необходимо и достаточно экспериментально измерить амплитуды колебаний масс Х1, Х2 и основания Х3. Если колебаниями основания можно пренебречь, то амплитуду Х3 следует положить в (3) равной нулю.

На втором этапе осуществляется подбор новых каскадных жесткостей для максимизации отношения амплитуд (Х1 / Х3) вплоть до нормативного значения и выше. В случае, когда колебания основания (корпуса) могут оказать существенное влияние на результат решения системы (2), необходимо дополнительно определить обобщенные коэффициенты С3 и М3.

В зависимости от конкретных условий определение обобщенной массы и обобщенной жесткости возможно расчетно-экспериментальным способом. Он заключается в использовании конечно-элементной модели конструкции корпуса либо ее фрагмента. Фрагмент выбирается в таких границах, чтобы дальнейшее их расширение не сказывалось существенным образом на получаемом результате, т.е. не изменяло результат за пределы допустимой погрешности. Затем выполняется прямой расчет конструкции на действие вертикальной единичной силы, приложенной в месте опирания оборудования А (рис. 1).

Обобщенная жесткость С3 может быть получена как обратная величина податливости конструкции в точке приложения силы. Экспериментальный путь определения С3 предполагает непосредственное измерение податливости относительно выбранной (несмещаемой) базы экспериментальными средствами. Аналогичный экспериментальный способ можно применить и для определения Ст.. В этом случае, устанавливая вибратор массой Мv в точку В отсоединенного трубопровода и определяя амплитуду Хв, можно определить и Мт:

Мт =

Из (2) также следует:

Мз = [Хз •

Найденные значения С2, С3 и М3 позволяют определить жесткости и произвести подбор амортизаторов при условии (Х[/Х3) = / (С, М, Г, ш-) ^ тах > (Х1/Х3)доп. для заданной области допустимых значений перемещений и жесткостей Сгт1П < С^ < Сгтах, (/=1,2); шт1П < ш < ш тах, где (Х1/Х3)доп - нормативное (допускаемое) значение каскадного перепада амплитуд колебаний. Вариация частоты ш при решении системы (2) с новыми каскадными жесткостями обеспечивает оценку и контроль резонансных частот системы.

В общем случае объекты колебательной системы могут иметь более сложную структуру, что предопределяет необходимость достаточно сложного в реализации расчетного анализа пространственных деформаций конструктивных элементов оборудования (механизмов), промежуточных рам и корпусных конструкций судна (рис.5). Пространственные расчетные модели приводят к существенному усложнению расчетного анализа. Поэтому в указанном случае целесообразно построить преобразованную расчетную схему колебательной системы, объединяющую ряд узлов приведения, назначаемых в местах опирания, концентрации масс и приложения изменяющейся внешней нагрузки. В зависимости от характера решаемой задачи в узлах приведения необходимо сохранять от 1 до 6 степеней свободы.

Матричные равенства (1) позволяют распространить изложенную расчетную методику на указанный случай. С этой целью их можно представить в следующем виде [3-4]:

\Срр II = I\Крр II- II • II 1II ;

\М л М»1+1и -I и

И 414 -IМ

1 2 •

- Я

Яр

где ||СРР|| и |МРР|| - матрицы коэффициентов жесткости узлов приведения; ||^рр|| , |МРР|| -аналогичные матрицы коэффициентов в исходных матрицах колебательной системы; ||^ЯЯ||, |Мяп|| - блок - матрицы коэффициентов жесткости и массы исключаемых узлов (области модели Я) системы; ||^рЯ|| и ||^Яр|| = ||КрЯ|| - блок- матрицы коэффициентов жесткости взаимного влияния узлов.

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

В случае гармонического внешнего возбуждения средством структуризации системы могут служить матрично-векторные равенства, объединяющие динамические матрицы жесткости конструкции:

Спр -ш2 -Мпр

Срр-2 -МРР

Сри -ш2 -МрП

С„„ -ш2 -М„

-1

Срй -ш2 -МрП

т

т

где ||Спр||, ||Мпр|| - матрицы приведенных коэффициентов жесткости и массы; р, п обозначают соответственно сохраняемые и исключаемые степени свободы (р << п); Т - знак транспонирования. Применение матричных равенств типа (6)-(8) значительно расширяет спектр решаемых задач, поскольку позволяет учитывать жесткостные и инерционные свойства всех компонент колебательной системы, включая корпус судна.

На рис. 3 показана двухкаскадная система виброизоляции электродвигателя, массой 50 кг, состоящая из основной и промежуточной опорных плит (массой 80 кг каждая) и восьми амортизаторов (жесткостью 316667 н/м).

На рис.4 приведена расчетная схема системы, установленной на консольно закрепленной платформе, состоящей из 80 стержней, длиной 0,1м с моментами инерции относительно главных центральных осей 12 = 3,33Е-9 м4, 1Х = 5,33Е-8 м4, 1у = 5,33Е-8 м4. Материал -сталь с модулем Юнга Е = 2 1011 Па, коэффициентом Пуассона ц= 0,3. В расчетах рассматривалось гармоническое (синусное) силовое воздействие на систему (рис. 4). Расчет выполнялся с применением специализированной программы «Парус-2000».

Рис. 3. Виброизоляция электродвигателя

Рис. 4. Расчетная схема перекрытия

Преобразованная расчетная схема объединяла вертикальные перемещения крайних узлов I, II и III (рис. 4) конструкции. На рис. 5 представлены графики гармонического отклика системы в указанных узлах в диапазоне частот 10-40 Гц.

Жесткость амортизатора 316667 н/м

Жесткость амортизатора 158334 н/м

п

е з

6 , 00E-03

5 , 00E-03

| 4 , 00E-03

s

х

3 3, 00E-03 ш I ш

щ 2, 00E-03 с

1, 00E-03 0 , 00E+00

10

LJL

15

20 25 30 Частота, Гц

35

40

1,00E-02

9,00E-03

8,00E-03

s 7,00E-03

1 6,00E-03 a>

| 5,00E-03 г

g, 4,00E-03 0>

c 3,00E-03 2,00E-03 1,00E-03 0,00E+00

.. j 1. 1. 1  / ,

20 30

Частота, Гц

П

е з

2,00E-02

1,80E-02

1,60E-02

E 1,40E-02

I 1,20E-02 z

=T 1,00E-02 o>

1 8,00E-03 a.

1= 6,00E-03 4,00E-03 2,00E-03 0,00E+00

1

1 1 1

20 25 30 35 Частота, Гц

3,50E-02 3,00E-02 2,50E-02 2,00E-02 1,50E-02 1,00E-02 5,00E-03 0,00E+00

1, . 1 ^ - »

10

20 30

Частота, Гц

40

е з

1,40E-02

1,20E-02

s 1,00E-02 к

I 8,00E-03 ш

I 6,00E-03 ш а ш

c 4,00E-03 2,00E-03 0,00E+00

/ ч 1 »

10

15

20 25 30 35 Частота, Гц

40

1,40E-02 1,20E-02 2 1,00E-02 8,00E-03

Ei

| 6,00E-03

а

0)

c 4,00E-03 2,00E-03 0,00E+00

1 V Л / f

10

20 30

Частота, Гц

40

40

10

40

Рис. 5. Графики гармонических откликов конструкции в узлах I, II, III при действии единичной гармонической вынуждающей силы от электродвигателя с эксцентриком

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

XJXn С

XX

XX

чо

ЧО ЧО

H3.00E+01 * 2,00E+01 1,00E+01

1

1 1 1

10

20 25 30

Частота, Гц

5,00E+01 4,00E+01 н 3,00E+01 ^ 2,00E+01 1,00E+01 0,00E+00

к \ \

20 25 30

Частота, Гц

m m oo in

1,20E+02 1,00E+02 8,00E+01 * 6,00E+01 4,00E+01 2,00E+01 0,00E+00

t> Ii f

20 25 30

Частота, Гц

7,00E+01 6,00E+01 5,00E+01 0 4,00E+01 * 3,00E+01 2,00E+01 1,00E+01 0,00E+00

f 1 i i j

20 25 30

Частота, Гц

5,00E+01

4.00E+01

0,00E+00

35

40

5

35

40

10

15

35

40

0

15

35

40

Рис. 6. Каскадные перепады амплитуд (X) вынужденных колебаний конструкции

в узлах I, II и I, III

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

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

1. Шиманский, Ю.А. Динамический расчет судовых конструкций / Ю.А. Шиманский. - Л.: Судпрогиз, 1948. - 407 с.

2. Сутырин, В.И. Варианты конденсации в практике расчетов динамики конструкций методом конечных элементов // Известия вузов. Машиностроение. 1990. №8. С. 7-11.

3. Сутырин, В.И. Применение способа фронтальной конденсации в расчетах динамики судовых корпусных конструкций // Судостроение. 1996. №7. С. 10-14.

4. Сутырин, В.И. Метод приведения в практике расчетов корпусных конструкций // Судостроение. 2008. №4. С. 17-22.

5. Сутырин, В.И. Методология конечно-элементного анализа судовых корпусных конструкций при их проектировании // Труды НГТУ им. Р.Е. Алексеева. - Нижний Новгород, 2010. №2 (81). С. 171-174.

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

V.I. Sutyrin1, E.I. Korotkaya2

THE SETTLEMENT-EXPERIMENTAL TECHNIQUE OF ACHIEVE STANDARD CASCADE DIFFERENCE OF VIBRATION INSULATION FOR THE AMORTIZED SHIP EQUIPMENT

Immanuel Kant Baltic Federal University, Kaliningrad1, Baltic State Academy of Fishing Fleet, Kaliningrad2

The settlement-experimental technique of achieve standard cascade difference of vibration insulation for the amortized ship equipment is offered and realized in given paper. Main goal consists of finding of cascade rigidity which returns system of the amortized equipment in limits of standard cascade change of vibration parameters. Solution of such task is possible by settlement and experimental way which is based on a structural quasistatic model. This method consists of construction of the transformed settlement scheme and unites limited number of checkouts of initial finalelement splitting in analyzed object. Also results of the solution test tasks are presented and analyzed in paper. The made researches confirmed possibilities of preservation of high precision received result.

Key words: Fluctuation analysis, ship hull structures, finite elements method, quasistatic model, condensation reception.