К.т.н. Тихонов В.А.

Институт машиноведения Российской академии наук, Москва, Россия

Способы расчета слоистых резинометаллических виброизоляторов

Предлагаются варианты расчета тонкослойных резинометаллических элементов – ТРМЭ, как упругих виброизолирующих элементов, двумя способами: аналитическим и численным методами решения в разных задачах. Линейные задачи статики и динамики решаются, чаще всего, полученными для тонких слоев резины асимптотическими аналитическими решениями. Нелинейные задачи, такие как устойчивость слоистого пакета при сжатии, косослой, большие деформации  – с применением модифицированного метода решения краевой задачи. Методики использованы при проектировании конструкций патрубков трубопроводов высокого давления [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.