Введение 4
1. Постановка задачи и численное исследование 6
2. Сравнение с полуаналитическим результатом Коллатца 10
3. Сравнение с прямоугольной пластинкой 11
4. Возможности аналитического решения 12
Заключение 13
Список литературы 13
Приложение 1
Работа посвящена исследованию колебаний тонких изотропных эллиптических пластин. Для описания пластин используется линейная теория тонкостенных пластин и оболочек Кирхгофа – Лява. При решении задачи нахождения собственных частот формулируется краевая задача для уравнения Жермен – Лагранжа. Это уравнение можно решить приближенно или аналитически с использованием собственных функций для двойного лапласиана в эллиптической области (функций Матьё).
Элементы в форме эллиптических пластин широко распространены в строительстве, приборостроении и автомобилестроении: в накопителях на жестком диске, токарных инструментах, датчиках давления, как компоненты инженерных структур [1, 4], в многослойных термовязкоупругих пластинах [17]. Колеблющиеся эллиптические пластины используются в акустических и музыкальных инструментах: сабвуферы, усилители, проигрыватели. Эллиптические пластины одни из немногих, для которых можно построить аналитическое решение уравнения колебаний.
Первая статья, посвященная сравнению колебаний круглой и эллиптической пластин, была опубликована в 1880 году (Francis E. Cabot [2]). Francis E. Cabot сделал предположение о характере колебаний эллиптической пластины и о колебании диска с круглым вырезом в центре. Однако до начала прошлого века не уделялось достаточное внимание аналитическому описанию колебаний эллиптических пластин. В исследованиях XIX века использовались лишь приближенные формулы и численные методы.
Поперечные колебания неоднородных эллиптических пластин рассмотрены во многих работах. Например, в 1992 году B. Singh, S. Chakraverty [13] описали решение задачи о нахождении частот поперечных колебаний эллиптической и круглой неоднородных пластин методом Рэлея-Ритца, используя в качестве функции формы полиномиальную функцию, обращающуюся в нуль на границе эллиптической пластины и обладающую нулевой производной на этой границе. В статье [1] N. Bhardwaj, A. P. Guptaand K. K. Choong исследовали первые шесть собственных частот для пластины, имеющей форму четверти эллипса, разрезанного по главным осям, которая покоится на упругом основании Винклера. Они использовали ортогональные полиномы, которые были получены при использовании процесса Грама-Шмидта, в качестве функций формы для метода Рэлея-Ритца. В 2005 году S. Chakraverty, RaginiJindal, V. K. Agarwal [4] изучали поведение собственных частот эллиптической пластины в зависимости от неоднородности пластины с помощью метода Рэлея-Ритца, используя в качестве функций формы ортогональные полиномы, которые были получены при использовании процесса Грама-Шмидта.
Аналитическое решение для задачи о поперечных колебаниях может быть построено только для пластинок простых форм. Колебания пластин круглых, квадратных и прямоугольных форм хорошо изучены и представлены в различных литературных источниках, например, в справочнике Arthur W. Leisa [8].
Самые ранние работы по нахождению колебаний эллиптической пластинки с помощью функций Матьё были проведены YoshioShibaoka [11], N. W. McLachlan [9]. Статьи этих авторов были опубликованы в середине прошлого века. Shibaoka утверждал, что N. W. McLachlan допустил серьезную ошибку, используя при решении уравнения колебания для эллиптической пластины тот же способ, что и для круглой пластины. N. W. McLachlanпредположил, что для вычисления любой моды колебания эллиптической пластины достаточно получить одно частное решение дифференциального уравнения колебания пластины. В случае поперечных колебаний эллиптической пластины для нахождения каждой частоты необходимо рассмотреть линейную композицию всех допустимых частных решений. Работа YoshioShibaoka была использована KenzoSato [6] для нахождения частот эллиптической пластины с упруго защемленным краем.
Исследование колебаний эллиптических пластин аналитическим методом с использованием функций Матьё было проведено в 2014 году S. M. Hasheminejad, AliGhaheri [5]. Они изучили изменение собственных частот эллиптической и круглой пластин с круглым и эллиптическим вырезом в зависимости от расположения выреза и угла поворота для эллиптического выреза.
Вычисление колебаний аналитическим способом для большинства форм пластинок является затруднительным, поэтому часто используются численные методы. Такие методы описаны в книге Л. Д. Акуленко и С. В. Нестерова [14]. Самыми распространенными приближенными методами являются: метод Рэлея-Ритца [1, 3, 4, 10, 12, 13] и метод конечных элементов.
Весьма важной является идея использования форм колебаний круглой пластины в качестве в качестве функций форм для метода Рэлея-Ритца при изучении пластинок других форм. Многие авторы используют этот метод для нахождения частот колебания эллиптических пластин. Так, например, в 1993 C. Ragalingham, R. B. Bhat и G. D. Xistris в статье [10] изучили зависимость спектра собственных частот жестко закрепленной эллиптической пластины от соотношения полуосей.
Метод конечных элементов является наиболее распространенным способом решения задачи на нахождение собственных колебаний пластинок разных форм. T. LakshmiReddy, P.V. PawanKumar, AkshayPrajapati в статье [7] провели численный анализ частот эллиптической пластины с жестко закрепленной границей в зависимости от соотношения полуосей, используя программный пакет ANSYS, и сравнили с результатами других авторов использовавших метод Рэлея-Ритца [3, 12]. А.Л. Смирнов в статье [18] получил спектр частот колебания свободной эллиптической пластины и исследовал его изменения в зависимости от соотношения большой и малой полуосей. В 2018 Е. А. Долгова в своей дипломной работе [16] изучала поведение собственных частот колебаний эллиптических пластин в зависимости от эксцентриситета и от размера центрального выреза, используя программный пакет Ansys.
При сохранении общего объема эллиптической пластины увеличение ее эксцентриситета приводит к росту всех собственных частот, связанному с увеличением жесткости конструкции при ее неизменной массе. При этом наиболее быстрый рост наблюдается у некратных частот, которые в случае нулевого эксцентриситета соответствуют осесимметричным колебаниям круглой пластины (формы c волновыми числами (0,n)). Частоты, соответствующие формам (m,n) круглой пластины, кратные для m > 0 в силу поворотной симметрии, расщепляются и растут с разной скоростью при увеличении эксцентриситета, причем быстрее растут частоты форм, вытянутых вдоль малой оси. Если в качестве параметра эллиптичности использовать отношение малой и большой осей k, то зависимость нижних частот колебаний при малых значенияхk близка к линейной.
[1] N. Bhardwaj A. P. Gupta, Choong K. K. Vibration of Rectangular Orthotropic Quarter Elliptic Plates with Simply-Supported Curved Boundary and Other Complicated Effects // Tamkang Journal of Science and Engineering. 2006. Noveber. Vol. 9, no. 1. P. 1-18.
[2] E. F. Cabot. Vibrations of Circular and Elliptical Plates // Proceedings
of the American Academy of Arts and Sciences. 1879. May. Vol. 15.
P. 219-222.
[3] Chakraverty S, M. Petyt, Natural frequencies for free vibration of non homogeneous elliptic and circular plates using two dimensional orthogonal polynomials. AppliedMathematicalModelling 07/1997: 21(7-21): 399-417.
[4] S. Chakraverty, Ragini Jindal, V. K. Agarwal. Flexural vibrations of non-homogeneous elliptic plates // Indian Journal of Engineering and Materials Sciences. 2005. December. Vol. 12. P. 521-528.
[5] S. M. Hasheminejad, Ali Ghaheri. Exact solution for free vibration analysis of an eccentric elliptical plate // Archive of Applied Mechanics. 2014. April. Vol. 84. No. 4. P. 543-552
[6] Kenzo Sato. Free flexural vibrations of an elliptical plate with edge restrained elastically // JSME. 1976. March. Vol. 19. No 129. P. 260-264.
[7] T. Lakshmi Reddy, P.V. Pawan Kumar, AkshayPrajapati. Modal analysis of an elliptical plate clamped along its boundary // International Research Journal of Engineering and Technology. 2015. Dec. Vol. 10. Issue. 9. P. 2030-2034.
[8] A. W. Leissa. Vibration of plates. Office of Technology utilization. 1969
[9] McLachlan N. W. Vibrational problems in elliptical coordinates // University of London. 1947. Vol. 5, no. 3. P. 289-297.
[10] C. Ragalingham, R. B. Bhat, G. D. Xistris. Vibration of clamped elliptical plates using exact circular plate modes as shape functions in Rayleigh-Ritz method // Int. J. Mech. Sci. 1994. Vol. 36. No. 3 P. 231-246.
[11] Shibaoka Y. On the Transverse Vibration of an Elliptic Plate with Clamped Edge // Journal of the Physical Society of Japan. 1956. July. Vol. 11, no. 7.
[12] Singh B, Chakraverty S, Use of characteristic orthogonal polynomials in two dimensions for transverse vibration of elliptic and circular plates with variable thickness. Journal of Sound and Vibration 1994: 173 (3): 289-99.
[13] B. Singh S. Chakraverty. Trannseverse vibration of circular elliptic plates with quadratically varying thickness // Appl. Math. 1992. May. Vol. 16. P. 269-274.
[14] Leonid D. Akulenko, Sergei V. Nesterov. High-Precision methods in eigenvalue problems and their applications. //2005
[15] И.А. Биргер, Я.Г. Пановко, Прочность. Устойчивость. Колебания. Том. 3 // М.: Машиностроение. 1963. Стр. 395-397
[16] Е. А. Долгова, “Свободные колебания тонких эллиптических пластин”, СПбГУ, 2018
[17] В. А. Ковалев. Динамика многослойных термовязкоупругих пластин // Изв. Сарат. ун-та. 2009. Том 9. Вып. 4. Ч. 1. Стр. 61-78.
[18] A. L. Smirnov. Free vibrations of annular circular and elliptic plates // EccomasProceedia COMPDYN. 2019. June. P. 3547-3555.