close

Вход

Забыли?

вход по аккаунту

?

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

код для вставкиСкачать
Механика и машиностроение
УДК 621.44.533.697
СОБСТВЕННЫЕ ЧАСТОТЫ ТРЕХСЛОЙНОЙ ПАНЕЛИ С ЛЕГКИМ ЗАПОЛНИТЕЛЕМ
© Н.В. Осадчий1, В.А. Малышев2, В.Т. Шепель3
1,3
ОАО «Научно-производственное объединение “Сатурн”»,
152903, Россия, Ярославская обл., г. Рыбинск, пр. Ленина, 163.
2
Рыбинский государственный авиационный технический университет им. П.А. Соловьёва,
152934, Россия, Ярославская обл., г. Рыбинск, ул. Пушкина, 53.
Расчет собственных частот панелей звукопоглощающих конструкций с сотовым заполнителем сведен к модели
трехслойной панели, толщина легкого заполнителя которой существенно больше толщины обшивок. Предложены численные методы расчета собственных частот, в частности метод сеток и метод конечного элемента. Погрешность оценки собственных частот колебаний, полученных методом конечных элементов и методом сеток, не
превышает 5%. Получено также аналитическое выражение для оценки главной собственной частоты квадратной
трехслойной панели, позволившее протестировать указанные численные методы.
Ил. 4. Табл. 1. Библиогр. 2 назв.
Ключевые слова: панель звукопоглощающей конструкции; собственная частота; система дифференциальных уравнений.
EIGEN FREQUENCIES OF SANDWICH PANEL WITH LIGHT FILLER
N.V. Osadchii, V.A. Malyshev, V.T. Shepel
NPO Saturn JSC,
163 Lenin pr., Rybinsk, Yaroslavl region, 152903, Russia.
Solovyov Rybinsk State Aviation Technical University,
53 Pushkin St., Rybinsk, Yaroslavl region, 152934, Russia.
Calculation of eigen frequencies of liner panels with honeycomb core is reduced to the model of a sandwich panel, the
light filler thickness of which significantly exceeds the thickness of paneling. The article proposes numerical methods to
calculate normal frequencies, in particular, a net method and a finite element method. Estimation error for eigen oscill ation frequency obtained by the finite element method and the net method does not exceed 5%. An analytic expression for
the estimation of the main eigen frequency of a square sandwich panel has been also derived. It allowed to test the spe cified numerical methods.
3 figures. 1 table. 2 sources.
Key words: liner panel; eigen frequency; system of differential equations.
Постановка задачи
Некоторые элементы авиационных конструкций, например панели звукопоглощающих конструкций (ЗПК) с
сотовым заполнителем, могут быть представлены в виде трехслойных панелей, толщина легкого заполнителя
которых существенно больше толщины обшивок (рис. 1). Легкий заполнитель – это заполнитель, который не сопротивляется растяжению/сжатию вдоль продольной оси образца.
а)
б)
Рис. 1. Панель в плане (а) и ее поперечное сечение (б)
___________________________
1
Осадчий Николай Васильевич, кандидат технических наук, эксперт конструкторского отдела прочности, тел.: 89206522794.
Osadchii Nikolai, Candidate of technical sciences, Expert of the Construction Department of Durability, tel.: 89206522794.
2
Малышев Владимир Александрович, доктор физико-математических наук, профессор кафедры математики, тел.:
89051378949, e-mail: malysheva@mail.ru
Malyshev Vladimir, Doctor of Physical and Mathematical sciences, Professor of the Mathematics Department, tel.: 89051378949, email: malysheva@mail.ru
3
Шепель Вячеслав Тимофеевич, доктор технических наук, профессор, начальник конструкторского отдела сертификации,
тел.: 89605386407, e-mail: sshepel@yandex.ru
Shepel Vyacheslav, Doctor of technical sciences, Professor, Head of the Construction Department of Certification, tel.: 89605386407,
e-mail: sshepel@yandex.ru
ВЕСТНИК ИрГТУ №9 (92) 2014
45
Механика и машиностроение
На рис. 1  – толщина обшивок трехслойной панели; h – расстояние между срединными поверхностями
обшивок.
Знание параметров собственных частот ЗПК является важным моментом их динамического исследования,
поскольку с их помощью можно определить степень опасности возможных резонансных режимов. Поэтому в
данной статье ставится задача аналитической и численной оценки собственных частот трехслойной панели с
легким заполнителем.
В соответствии с работой [1] уравнения динамики трехслойной панели могут быть записаны в виде
M y M xy
M x M xy
Qx Qy
2w

 Qx  0;

 Qy  0;

F  2 ,
x
y
y
x
x
y
t
(1)
где: Mx, My, Mxy – изгибающие моменты; Qx, Qy – поперечные силы; F – интенсивность внешней нагрузки;
w – прогиб;  – плотность.
Выражения для изгибающих моментов и поперечных сил имеют вид:
 

Mx  D


x
y


 
 

; M y  D 
; M

y
x 


xy
   
 D(1   ) 


 y x 

w 
w 

Qx  Gh   
,
 ; Qy  Gh   
x 
y 


(2)
(3)
где  и  – относительные смещения верхнего и нижнего слоев по осям x и y ; E – модуль упругости; 0    1 –
коэффициент Пуассона; G – модуль сдвига в плоскости панели; D – цилиндрическая жесткость, вычисляемая по
формуле
D  Eh2 / 2(1   2 ) .
Подставляя выражения (2) и (3) в уравнение (1), получаем систему дифференциальных уравнений:
 2
 2  2  Gh
Gh w

(1


)



 0;
2
2
x
y
xy D
D x
2
 2   2 Gh
Gh w

(1


)



 0;
y 2
x 2 xy D
D y
(4)
   2 w  2 w F
 2w

 2  2

.
x y x
y
Gh Gh t 2
Решение системы уравнений (4) строится в прямоугольнике
0  х  а; 0  у  b.
Условия для прогиба w(x, y, t) и смещений (x, y, t) и  (x, y, t) задаются в виде краевых условий Дирихле:
w(0, y, t )  w(a, y, t )  w( x,0, t )  w( x, b, t )  0;
 (0, y, t )   (a, y, t )   ( x,0, t )   ( x, b, t )  0;
 (0, y, t )   (a, y, t )   ( x,0, t )   ( x, b, t )  0.
(5)
Задача определения собственных частот панели решается при нулевой нагрузке F(x, y, t) = 0 и следующих
начальных условиях:
w( x, y,0)  w0 ( x, y );
 ( x, y,0)   0 ( x, y);
46

w( x, y,0)  0;
t
 ( x, y,0)   0 ( x, y).
ВЕСТНИК ИрГТУ №9 (92) 2014
(6)
Механика и машиностроение
Начальный прогиб w0(x, y) и смещения 0(x, y) и 0(x, y) определяются из решения задачи статики для
начальной нагрузки F(x, y) = F0(x, y). В рамках краевых условий Дирихле решение задачи статики сводится к решению системы уравнений (4), в которой формально полагается  = 0.
Собственные частоты трехслойной панели
В случае краевых условий Дирихле (5) решение системы дифференциальных уравнений (4) с начальными
условиями (6) выражается через базисные функции W(x, y), A(x, y), B(x, y) системы дифференциальных уравнений:
2 A
 2 A  2 B Gh
Gh W

(1


)


A
0
2
2
x
y
xy D
D x
2B
 2 B  2 A Gh
Gh W
 (1   ) 2 

B
0
2
y
x
xy D
D y
(7)
A B  2W  2W


 2  W  0
x y x 2
y
Величины , для которых система дифференциальных уравнений (7) имеет ненулевые решения W(x, y), A(x,
y), B(x, y), называются собственными числами. Проверим, все ли собственные числа рассматриваемой краевой
задачи Дирихле строго положительны:   0.
Зададим функции
P
2 A
2 A 2 B
2 B
2 B 2 A



,
Q




,
*
*
x 2
y 2 xy
y 2
x 2 xy
где *  1   . Запишем третье соотношение системы (7) в виде
P Q
Gh


W.
x y
D
Отсюда
 P
Q

  x W  y W    
Gh
W2 ,
D 
где интегрирование выполняется по прямоугольнику 0  х  а; 0  у  b. Интегрируя по частям с учетом краевых
условий Дирихле, получаем:
 P
Q

 W
  x W  y W      P x
Q
W
y

  D
D

 D

P  A  Q
Q  B   
( P 2  Q 2 )   ( PA  QB).
    P 
Gh 

 Gh

  Gh

Далее
2
2
 A 2
  2 A

2 A 2B 
 A  B A  B  

   *   
 2  * 2 
 A 
 
y
xy 
y  y x  y  
 x
 x 



 ( PA  QB)      2 B  2 B  2 A     

2

B

A

B



 .
   2  * 2 



B
 *  x  x y

x
xy  
  y
2
2
 A B  2
 A 
 B  
   

  *    * 
 
 x  
 x y 
 y 
В результате
 W 2 
D2
G 2 h2
2
2
 (P  Q ) 
2
2
2
 A 
D  A B 
 B  









*
*
 .
Gh   x y 
 x  
 y 
ВЕСТНИК ИрГТУ №9 (92) 2014
(8)
47
Механика и машиностроение
Если A(x, y) = 0 и B(x, y) = 0, то из первых двух уравнений системы (7) следует, что W(x, y) = 0. Поэтому данный случай не рассматривается. Если А  0 или В  0, то из условий Дирихле и равенства (8) следует, что W  0.
Поэтому, действительно, все собственные числа  краевой задачи Дирихле трехслойной панели строго положительны, то есть   0.
Возьмем произвольное собственное число  краевой задачи Дирихле трехслойной панели и соответствующую данному числу тройку собственных функций W(x, y), A(x, y), B(x, y). Зададим тройку функций:
w( x, y, t )  W ( x, y )cos(t );
 ( x, y, t )  A( x, y )cos(t );
 ( x, y, t )  B( x, y )cos(t ),
где
 2  Gh /  .
Непосредственно проверяется, что тройка функций w(x, y, t), (x, y, t),  (x, y, t) удовлетворяет системе дифференциальных уравнений (4) при нулевой нагрузке F(x, y, t) = 0 и начальных условиях:
w0(x, y) = W(x, y); 0 = A(x, y); 0(x, y) = B(x, y).
Данная тройка функций w(x, y, t), (x, y, t),  (x, y, t) задает колебание трехслойной панели с собственной частотой
   / 2 .
Возьмем минимальное собственное число min краевой задачи Дирихле трехслойной панели и соответствующую данному числу величину
2
min
 Ghmin /  .
Тогда величина
 min  min / 2
является главной собственной частотой трехслойной панели.
Процедура вычисления собственных частот
Процедура вычисления собственных частот трехслойной панели основана на идеологии метода сеток [2].
Осуществим замену дифференциальных операторов в системе уравнений (7) на разностные операторы по следующим правилам:
f
 fi 1, j f
f
f
f
 i 1, j
;  i , j 1 i , j 1 ;
x
2x
y
2y
2
fi 1, j  2 fi , j  fi 1, j  2 f
f
f
f
f
 f

;
 i 1, j 1 i 1, j 1 i 1, j 1 i 1, j 1 ;
2
2
x
x
xy
4xy
2
f
 2 fi , j  fi , j 1
 f
 i , j 1
,
2
y
y 2
где 0  i; j  N + 1. Отметим, что в силу краевых условий Дирихле возникающие системы линейных уравнений
решаются в предположении, что
W0, j  WN 1, j  Wi ,0  Wi , N 1  0;
A0, j  AN 1, j  Ai ,0  Ai , N 1  0;
B0, j  BN 1, j  Bi ,0  Bi , N 1  0.
Используя переменные Wi, j, где 1 i, j  N, как параметры, решим систему линейных уравнений, соответствующую первым двум соотношениям системы дифференциальных уравнений (7), относительно переменных Ai , j и
Bi , j , где 1  i, j  N. После подстановки выражений линейного характера для переменных Ai , j и Bi , j через пере-
менные Wi , j в систему линейных уравнений, соответствующую третьему соотношению системы дифференци-
48
ВЕСТНИК ИрГТУ №9 (92) 2014
Механика и машиностроение
альных уравнений (7), получаем систему линейных уравнений вида
LW = W
относительно вектора W размерности N с переменными Wi, j, где 1  i, j  N.
Расположим собственные числа матрицы L в порядке возрастания
2
0  1  2  ...  N 2
.
Тогда при достаточно большом значении N величины
0   1   2  ...   N 2
,
где
k 
1
2
Gh

k ,
2
могут служить приближениями к первым N собственным частотам трехслойной панели. В частности, приближением к главной частоте  min трехслойной панели может служить величина  1 .
Полезно отметить, что в частном случае при N = 2 квадратной панели a  b удается получить явное выражение для величины
1 
1
2
Gh 54 39 D  12 D(1   )  a 2Gh
.
 a 2 117 D  36 D(1   )  4a 2Gh
(9)
Полученное выражение может служить полезной рабочей формулой для оценки главной частоты квадратной
трехслойной панели.
Численный эксперимент
Проиллюстрируем полученные результаты на примере трехслойной панели с параметрами:
4
2
2
-5
3
а = b = 100 мм; h = 10 мм; τ = 1 мм; Е = 2·10 кгс/мм ; G = 100 кгс/мм ;  = 0,3;  = 1,1·10 кг/мм .
В случае N = 7 первые 49 = 77 собственные частоты данной трехслойной панели расположены по строкам
следующей 77-матрицы:
 63
161

206

 241
255

273
 301

100 100 128 139 140 161 
176 176 188 193 193 206 
216 216 221 221 229 229 

214 241 242 242 243 243 
255 261 261 264 264 273 

282 282 285 293 293 301
312 312 317 327 327 337  .
Значения собственных частот округлены до целых значений Гц. Отметим, что точное значение главной частоты  1 равно 63,357 Гц. При этом приближенное значение главной частоты  1 , вычисленное по рабочей формуле (9), равно 62,719 Гц. В данном примере, в частности, рабочая формула (9) дает 1% относительной ошибки
при вычислении главной частоты. Кроме того, среди собственных частот могут быть кратные. Например, имеют
место равенства 2 = 3; 7 = 8; 9 = 10.
На рис. 2 изображены базисные формы колебаний Wk(x, y), соответствующие собственным частотам
1   2   3   4 .
С целью проверки точности определения главной собственной частоты было осуществлено моделирование
колебания трехслойной панели с указанными выше параметрами, при которых начальный прогиб w0(x, y) и смещения 0(x, y) и 0(x, y) соответствовали единичной внешней нагрузке F0(x, y) = 1. Вид начального прогиба изображен на рис. 3.
ВЕСТНИК ИрГТУ №9 (92) 2014
49
Механика и машиностроение
W1(x, y)
W3(x, y)
W2(x, y)
W4(x, y)
Рис. 2. Базисные формы колебаний
В процессе численного эксперимента решение системы дифференциальных уравнений (4) проводилось с
нулевой внешней нагрузкой F(x, y, t) = 0. При этом на квадрате 100100 была введена сетка, образованная из
400 (20·20) подквадратов размером 55, а шаг интегрирования по времени t был принят равным 0,0001 c.
На рис. 3,б изображен график первого периода колебания функции w(50, 50, t), отображающий зависимость
отклонения центра панели от нейтрального положения при начальном прогибе (рис. 3,а). По графику функции
w(50, 50, t) период колебаний T = 0,0161 с, и, следовательно, частота колебания v =1/T = 62,111 Гц. Как видим, в
данном примере частота колебания v совпадает с главной частотой v1 с точностью относительной ошибки 2%.
50
ВЕСТНИК ИрГТУ №9 (92) 2014
Механика и машиностроение
а)
б)
Рис. 3. Начальный прогиб (а) и график колебания центра панели (б)
Метод конечных элементов
Расчет собственных частот методом конечных элементов проводился также с помощью программы ANSYS
на объемной модели, которая состоит из следующих элементов:
– Solid 45 – для заполнителя;
– Shell 181 – для обшивок.
Рис. 4. Первые четыре формы собственных колебаний, определенные методом конечных элементов
ВЕСТНИК ИрГТУ №9 (92) 2014
51
Механика и машиностроение
Закрепление модели соответствовало краевым условиям (5). На рис. 4 представлены первые четыре формы
собственных колебаний, найденные методом конечных элементов. По характеру и числовым значениям они совпадают с базисными функциями, найденными методом сеток. В таблице представлены значения первых десяти
собственных частот, значения которых получены как методом сеток, так и методом конечных элементов. Из таблицы видно, что значения собственных частот, определенные методом конечных элементов, подчиняются условию (9), при этом погрешность определения не превышает 5%.
Значения собственных частот
Значения собственных частот, Гц, рассчитанные
Номер собственной
частоты
методом конечных элементов
методом сеток
Относительная
погрешность, %
1
60
63
5,08
2
95
100
4,92
3
95
100
4,92
4
123
128
4,22
5
138
139
0,42
6
139
140
0,57
7
159
161
1,09
8
159
161
1,09
9
184
176
4,20
10
184
176
4,20
Итак, в данной работе представлены результаты аналитической оценки и верифицированные на ее основе
вычислительные модели определения собственных частот трехслойных панелей с легким заполнителем. Визуальные картины собственных колебаний по первым четырем формам, определенным методом сеток и методом
конечных элементов, представлены на рис. 3, 4. Значения собственных частот, а также их относительная погрешность приведены в таблице. Погрешность в оценке собственных частот, определенных обоими методами, не
превышает 5%. Полученные результаты, по нашему мнению, могут быть использованы при конструировании и
доводке панелей звукопоглощающих конструкций, а также при их отстройке от возможных резонансных режимов.
Статья поступила 22.07.2014 г.
Библиографический список
1. Вольмир М.А. Гибкие пластины и оболочки. М.: Физматгиз, 1970. 312 с.
2. Самарский А.А. Теория разностных схем. М.: Наука, 1989. 616 с.
52
ВЕСТНИК ИрГТУ №9 (92) 2014
Документ
Категория
Без категории
Просмотров
8
Размер файла
3 515 Кб
Теги
заполнителей, частоты, трехслойного, панели, легких, собственных
1/--страниц
Пожаловаться на содержимое документа