Основы проектирования систем виброзащиты с упругими элементами из материала МР Текст научной статьи по специальности «Физика»

УДК 629.7

ОСНОВЫ ПРОЕКТИРОВАНИЯ СИСТЕМ ВИБРОЗАЩИТЫ С УПРУГИМИ ЭЛЕМЕНТАМИ ИЗ МАТЕРИАЛА МР

© 2008 А.М. Уланов, Ю.К. Пономарев

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

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

Виброизоляторы из материала МР (МР -"Металлическая резина" - получается холодным прессованием проволочной спирали) широко используются в виброзащитных системах [1, 2, 4]. Они имеют высокую прочность, высокий коэффициент рассеивания энергии, способны работать в условиях высокой и низкой температуры, в агрессивной среде, вакууме, радиации и так далее. Однако свойства материала МР нелинейно зависят от многих параметров: относительной

р _ рык

плотности материала у ~ (где р - плотность материала проволоки, рык - абсолютная плотность материала МР), амплитуды деформации А, предварительной статической деформации Q. Применяемые модели материала МР (их обзор присутствует в [2]) в большинстве случаев основаны на рассмотрении контактного взаимодействия проволок в МР, очень сложны и малопригодны для практического применения при проектировании виброизоляторов. В настоящей работе предлагается подход к проектированию систем виброзащиты, использующих материал МР, на основе представления о МР как анизотропной сплошной среде. Такой подход также открывает возможность применения конечно-элементных программ при расчете виброизоляторов из материала МР.

Определение параметров МР

Исследовались образцы в виде куба со стороной 1=25 мм, изготовленные из проволоки диаметром 0,1 мм. В принятой системе координат направление прессования - ось Х. Зависимость напряжения прессования от от-

носительной плотности материала можно представить уравнением

а

press

576р17

(1).

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

Коэффициент Пуассона

Коэффициент Уух (описывает влияние силы, приложенной в направлении X, на де-

формацию в направлении Y;

у =--

yx

)

очень мал. В рабочем диапазоне (е < 0.16) Уух < 0.03 . Есть зависимость от относитель-

Ух

ной деформации

^ - 0.2ех . (2)

Аналогично Ух . Поскольку коэффициент Пуассона в данном направлении очень мал, не имеет значения, приклеены или свободны нагрузочные поверхности образца.

При сжатии в направлении Y в направлениях X и 2 наблюдаются значительные петли гистерезиса (рис. 1).

Они приблизительно подобны петле в направлении Y. Если перестроить петли в координаты еу - ех (или, соответственно, еу - ег), результат можно описать приблизительно в виде ех - -аеу ± Ь, где коэффициент а играет роль коэффициента Пуассона, а коэффициент Ь является остаточной деформацией. Знак + относится к процессу нагрузки, знак - к процессу разгрузки.

Коэффициент Пуассона и остаточная деформация при свободных нагрузочных по-

s

Рис. 1. Петли гистерезиса в направлениях X и Z при сжатии в направлении У (приклеенный образец)

верхностях зависят от амплитуды деформации и относительной плотности материала.

« (-118у ±3.]£у)еур . (3) Аналогично V .

г, *(-8.5* + 0.46ev)е„р.

(4)

Аналогично Vxz.

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

Для приклеенных нагрузочных поверхностей вид зависимостей несколько отличается. 8Х *-0.668у ± (-0.015Р + 0.035), (5)

то есть коэффициент Пуассона vxy « 0.66 . Аналогично V . Точно так же

г, * -0.3* ± (-0.033Р + 0.014),

(6)

то есть коэффициент Пуассона vzy « 0.3 . Аналогично vyz.

При растяжении (экспериментально проверено до 8 « -0.05) деформация материала МР в боковых направлениях не изменяется, следовательно, коэффициент Пуассона равен 0.

Модуль сдвига G и сила трения

при сдвиге

Петля гистерезиса при сдвиге (кроме ее концов) близка к линейной (рис. 2), поэтому ее можно приблизительно описать уравнением

г - Gy ±т{, (1)

где г^ - напряжение от сил трения.

Из обработки большого количества измерений (для различной относительной плотности и различного диаметра проволоки) получены соотношения

(8)

1:1.26:1.99; (9)

Gw : G„ : Gy * 1:1.53:2.75

Рис. 2. Деформация образца из МР в одном из направлений сдвига

Gzx - Gyx; ^ - Gxy; Gyz - ^у;

ГШх — ГHyx ; ГHxz — ГHxy ; ГHyz — ГНУ (10)

Первый индекс - направление действия силы, второй - направление, в котором материал закреплен при сдвиге.

Зависимость от относительной плотности

« 6.65р17 МРа, гНух « 0.101Р11 МРа. (11)

Модуль сдвига и напряжение от сил трения в других направлениях можно получить, используя зависимости (8), (9) и (10).

Интересно, что в относительных координатах

о ^ - 6.«?: -

ух

а

press

576р

—1.7

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

Модуль упругости (модуль Юнга) и напряжение трения при сжатии и растяжении

Упругая линия петли гистерезиса на сжатие и растяжение существенно нелинейна. Целесообразно описать петлю как сумму напряжения на упругой линии аь и напряжения трения ан .

а - аь ± ( . (13)

Упругая линия и напряжение трения могут быть определены по экспериментально процессам нагрузки а1 и разгрузки а2 по уравнениям

а + а2

и а

H

соответ-

Т Hvx : Т Hxv : ТH,y

ь 2 Н 2 ств енно в диапазоне деформаций [0.1 Лтп; 0.1 Лтах ] [3]. Для направления Х этот диапазон составил 8Х е [-0.06;0.24], в направлениях Yи 2 8 е [-0.06; 0.16]. Линии ат

Рис. 3. Упругая линия и напряжение трения петли гистерезиса образца материала с р = 0.18 в направлениях X и У

и аН для относительной плотности р = 0.18 в направлениях Х и Y приведены на рис. 3.

Зависимость напряжения трения от £ удовлетворительно описывается полиномом степени 3, зависимость от р в направлении Х сходная с зависимостью при сдвиге (прямая пропорциональность), в направлениях Y и 2 зависимость более сильная (пропорцио-

нально р17

аНх « (0.0181 + 0.231£х + 0.0691£х2 +

+ 7.71£х3)

' Р

0.18

(14)

а

Ну

(0.0597 + 0.819£у + 6.51£у2

-16.8£3)

- Л17

Р

0.18

(15)

Аналогично а

Зависимость упругой линии от £ удовлетворительно описывается полиномом степени 4, зависимость от р в направлении Х сходна с а рге88 (пропорционально р17), в направлениях Y и 2 зависимость более сильная (пропорционально р2).

аЬх « (0.613 + 1.771£х - 20.4£х2 +

+160£_3)£

р ^

0.18

(16)

а

Ьу

(2.69 + 14.2£у +139£

( р \ 2 -632£3)£ р

0.18

(17)

Аналогично а

Lz •

Первый сомножитель в этих произведениях аналогичен модулю упругости Е.

Коэффициент рассеивания энергии в различных направлениях

Уравнения (7) - (11) и (13) - (17) дают возможность построить приблизительно петли гистерезиса (с вертикальными начальными участками нагрузки и разгрузки) и оценить коэффициент рассеивания энергии у . В направлении Y жесткость материала МР существенно больше, чем в направлении Х, а различие силы трения не так велико, поэтому коэффициент рассеивания энергии в направлении Х существенно больше. Например, для амлитуды деформации £ = 0.1 и относительной плотности р = 0.18 получаем

~ 1.4. С увеличением р сила трения в на-

Уу

правлении Y растет несколько быстрее, и различие уменьшается, но все равно при р = 0.3

¥х

получаем — «1.2.

Уу

При сдвиге петля гистерезиса близка к линейной, ее площадь пропорциональна 2ту, а потенциальная энергия деформирования

1 2

пропорциональна ^ УУ . Поэтому коэффициент рассеивания энергии пропорционален — , и из уравнений (8) и (9) следует уух :уху « 1: 0.82 : 0.72 , а из уравнений (11) следует, что у при сдвиге пропорционален (р) 06.

Значения ух и уух сопоставимы. При

£ = у = 0.1 отношение

¥

1.12 , при

£ = у = 0.15 отношение

¥

ух

Ух

0.82

х

Коэффициент трения

Поскольку значительная доля энергии рассеивается в виброизоляторе не при трении внутри материала МР, а при трении материала МР о стальные детали конструкции виброизолятора, необходимо знать коэффициент трения МР по стали/ Для этого образец прижимался известной силой Р1 и смещался силой Р2 в перпендикулярном направлении. Относительная плотность р и направление силы (Р1 в направлении прессования Х, Р2 в перпендикулярном направлении Y, и наоборот) на коэффициент трения не влияют. Обнаружена слабая зависимость / от напряжрния прижатия. При малых значениях о = —1 коэффициент трения/=0.1, далее он слабо растет и после о « 0.75МРа остается постоянным и равным 0,125.

Проектирование

виброзащитной системы

Рассмотрим пример проектирования виброзащитной системы на основе полученных зависимостей для параметров материала МР. Пусть виброизолятор работает на сжатие в направлении Х или Y . Даны: масса защищаемого объекта т, требуемая резонансная частота виброзащитной системы /0, амплитуда действующего ускорения а, допустимое ускорение защищаемого объекта на резонансе .

По этим данным можно, используя приближенную линейную модель, найти необходимую динамическую жесткость виброизолятора

kd = 4я2 /0 т, (18)

необходимый коэффициент передачи на резонансе

Ъ = (19)

амплитуду деформации виброизолятора на резонансе ^

А=^Л (20)

и требуемый коэффициент рассеивания энергии

2п

¥

Ъ

(21)

Для обеспечения большого времени работы виброизолятора желательно, чтобы его относительная деформация не превышала 10%. Примем е0 = 0.1. Тогда высота упругого элемента из МР а

(22)

Н = —

ческая деформация Q определяется весом защищаемого объекта. Однако она может быть увеличена при помощи дополнительного сжатия упругих элементов из материала МР. Получаем

°ь (еч), (23)

Q

где ед = н , д = 9.81т /э2, F - площадь поперечного сечения виброизолятора.

Таким образом, виброизолятор совершает колебания между точками с относительной деформацией

Q + А Q - А

е1 = —--е2 = ---(24)

1 Н ' 2 Н • К }

Тогда динамическая жесткость виброизолятора

К

°Ь (е1) + °Н (е1) - (оь (е2) - °Н (е2) F =

2 А

(АЬ + ь'н )

2А

F

где

АЬ = °Ь (е1) - °Ь (е2) '

(25)

(26)

Минимальная предварительная стати-

ЭН = °Н (е) + °Н (е2) . (27)

Подставив в (25) уравнения (18) и (20), получаем

2^0т = (Аь + Эн )F (28)

Площадь петли гистерезиса (предполагая, что участки первоначальной нагрузки и разгрузки - вертикальные линии) приблизительно равна

АW' « 2А¥(он (ег) + ОН (е2)) = 2А^ .

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

АW « 0.75АW1.5АЕэН .

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

W - 0.25(0 (е1) + ОН (е) - Оь (е0)) +

+ (0Ь (е0) - °Ь (е2) + °Н (е2)))^ = .

= 0.25(А ь + Эн ) AF

Тогда коэффициент рассеивания энергии

¥ =

¥

AW W Из 3Fs

6s г

A L + sff

(25)

A L + sH

2 Akd

F

H

Ak,

С учетом (18) - (21), получаем

FsH = —Г" •

Разделив почленно (28) на (29), получим

3Wo

(30)

A L + Sh

JH

ma

Получаем следующий алгоритм проектирования виброизолятора:

1. Рассчитать значения А и Н по уравнениям (20) и (22).

2. Принять некоторое значение р (для обеспечения минимального веса виброизолятора можно начинать с наименьших технологически достижимых значений, например, р = 0.18. При небольших требованиях к прочности виброизолятора возможно р = 0.15 .

3. Принимая ряд значений Q, получать по уравнениям (24) значения е1 и е2, по уравнениям (14) и (16) (или по уравнениям (15) и (17) для направления, перпендикулярного направлению прессования) значения аН (£1), аН (е2), аь , аь (е2), по уравнениям (26) и (27) значения Аь и sH , добиваться выполнения условия (30).

4. Если условие (30) выполнено, при помощи уравнения (29) найти значение F.

5. Проверить выполнение условия (23). Желательно также, чтобы г2 > 0, иначе необходимо приклеивать материал МР к опорным поверхностям, обеспечивая работу материала на растяжение.

6. Если при принятом значении р условия (30) или (23) не выполняются ни при каком Q, увеличить значение р и повторить Тогда пункты 3, 4, 5. На практике р < 0.35 , и виброизолятор с такой плотностью обычно уже является чрезмерно жестким. Если изменение значения р не дало результата, требования к виброизолятору невыполнимы. (29) Полученный алгоритм легко реализовать

в виде программы для компьютера. В результате получаем высоту H и площадь сечения F виброизолятора, его относительную плотность р и необходимую предварительную деформацию упругого элемента Q . Этого достаточно для изготовления виброизолятора.

Проектирование виброизолятора, работающего на сдвиг, ведется аналогичным образом, с использованием уравнений (7) - (11).

СПИСОК ЛИТЕРАТУРЫ

1. Hongrui Ao, Hongyuan Jiang, Ulanov A.M. Dry Friction Damping Characteristics of a Metallic Rubber Isolator under Two-dimensional Loading Process // Modelling and Simulation in Material Science and Engineering. 2005. №13.

2. Чегодаев Д.Е., Мулюкин О.П., Колтыгин Е.В. Конструирование рабочих органов машин и оборудования из упругопористо-го материала МР. Самара: СГАУ, 1994.

3. Ulanov A.M., Lazutkin G.V. Description of an Arbitrary Multi-axial Load Process for Non-linear Vibration Isolators // Journal of Sound and Vibration. 1997. №203(5).

4. Расчет и конструирование средств виброзащиты сухого трения. Самара: СамГАПС, 2005.

DESIGN OF VIBRATION PROTECTION SYSTEMS WITH ELASTIC ELEMENTS MADE OF MR MATERIAL

© 2008 A.M. Ulanov , Yu.K. Ponomarev

Samara State Aerospace University

MR (metal rubber) is non-linear anisotropic material for vibration protection with high strength and damping. An elastic modulus, shear modulus and Poisson ratio are obtained for this material for different directions. It allows to calculate a deformation of vibration isolator made of MR, to use for it a finite element software. An algorithm for design of MR vibration isolator is given. This algorithm allows to determine size of MR elastic element, its density and preliminary deformation.