К.т.н. Тихонов В.А.
Институт машиноведения Российской академии наук, Москва,
Россия
Способы расчета слоистых резинометаллических виброизоляторов
Предлагаются варианты расчета
тонкослойных резинометаллических элементов – ТРМЭ, как упругих виброизолирующих
элементов, двумя способами: аналитическим и численным методами решения в разных
задачах. Линейные задачи статики и динамики решаются, чаще всего, полученными
для тонких слоев резины асимптотическими аналитическими решениями. Нелинейные
задачи, такие как устойчивость слоистого пакета при сжатии, косослой, большие
деформации – с применением
модифицированного метода решения краевой задачи. Методики использованы при
проектировании конструкций патрубков трубопроводов высокого давления [1] .
Для анализа виброзащиты машин, оборудования, патрубков чаще всего ставится
общая задача оценки переходных частотных характеристик матрицы жесткости по
соотношению структурной механики между обобщенными силами и смещениями через
матрицу жесткости. При оценках виброизолирующих свойств патрубка наибольший
интерес представляют переходная жесткость и коэффициент передачи.
Задача определения
матрицы жесткости патрубка сводится к расчету наряженного состояния при сжатии,
изгибе и сдвиге, а также действии внутреннего избыточного давления в кольцевом
элементе, который, в свою очередь, после процедуры интегрирования по высоте
тонкого резинового слоя приводятся к краевой задаче Гельмгольца. Задачи
проектирования патрубка и расчета элементов рассмотрены в работах [1-3].
Решения относительно
гидростатического давления s в тонком резиновом слое
находятся из уравнения Гельмгольца на примере кольца:
где
G и K – модули сдвига и
объемного сжатия резины;

Напряжения и перемещения
выражаются через гидростатическое давление. Для линейных задач получены
аналитические решения уравнения в виде функций гидростатического давления. В
реальных условиях приходится решать геометрически и физически нелинейные
задачи, для которых предлагается использовать электронные средства пакета MathCad
, в котором предусмотрено численное решение нелинейных дифференциальных
уравнений. Краевая задача приводится к задаче Коши. Дифференциальное уравнение второго порядка приводится к
системе уравнений первого порядка и строится вектор-функция у:
, компоненты которой определяются так:
Относительно символьной
вектор-функции D (r,y)= {y1 , y2} решается краевая задача, которая с помощью
функций:
приводится
к решению задачи Коши методом Рунге-Кутта:
, где
x1 , x2 – граничные точки интервала
интегрирования – для кольца: ri = Ri / Re и
re = 1; v – вектор начального приближения; с помощью функции score задается нулевое граничное условие y0
– 0 на давление в задачах сжатия и изгиба и равное избыточному давлению y0
– pi
при его наличии, тогда приходится менять порядок краевых точек x2, x1 местами, в том числе и
в функции sbval(v,1,1, r,n,D), которая путем
разделения интервала на n
точек определяет недостающее второе граничное условие на производную ds (в задаче это – вектор) на другом конце
интервала. Решение задачи Коши с помощью функции rkfixed содержится в матрице Z, где первый столбец –
точки разбиения интервала по ширине кольца: rk , k=1,...,n, второй – значения
функции s (r) в этих точках, третий
– значения производной ds/dr, которое используется
для расчета касательных напряжений. При проектировании кольцевого элемента решаются
нелинейные задачи с косослоем и поперечной устойчивостью слоистого пакета.
Неравномерность
резинового слоя по толщине как технологическое несовершенство, выраженное в
непараллельности жестких тарелей и фланцев, задается начальным взаимным
поворотом на угол
тарелей и фланцев –
функцией толщины от полярных координат в виде:
где
– номинальное значение толщины резинового слоя постоянной
толщины (предыдущий вариант),
– полярная координата точки, соответствующей минимальной толщине
слоя при угловой координате изгибной деформации j =0,
при этом решение при изгибе ищется в виде s
= scosj. Напряженное состояние
в резиновом слое определяется решением двухточечной краевой задачи:
Так задача со слоями
постоянной толщины приведены к системе дифференциальных уравнений, решаемых
методом Рунге-Кутта. Данное представление позволяет получить решения нелинейных
задач с учетом зависимостей модулей упругости резины G и K (a = G/K) от интенсивности
деформаций и гидростатического давления. Иллюстрацию некоторых решений
можно найти в [3]. Решена задача потери
поперечной устойчивости при сжатии. Построена нелинейная кривая: нагрузка –
стрела «прогиба» (максимальное смещение сдвига слоя в пакете), аналогичная
нагрузочной кривой потери устойчивости сжатого стержня с начальной погибью.
Литература
1. Кирюхин А.В., Тихонов
В.А., Чистяков А.Г. Вибрационные характеристики углового компенсационного
патрубка с тонкослойными эластомерными элементами // Проблемы машиностроения и
надежности машин. 2007, №1, С. 103-108.
2.Тихонов В.А. Расчет вибрационной жесткости
сферического резинометаллического подшипника // Проблемы машиностроения
и надежности машин. 2004, 6, стр. 9-14.
3.Кирюхин А.В., Тихонов В.А., Чистяков А.Г.
Расчет нелинейных характеристик ТРМЭ при проектировании компенсационных
патрубков. – В сб. тр. XVII Симпозиума: Динамика виброударных (сильно
нелинейных) систем. «DYVIS-2012». Москва-Клин. 2012. С. 96-99.