Метод компенсирующих нагрузок для решения задач о циклически симметричном изгибе анизотропных пластин, контактирующих с упругим основанием

Обложка

Цитировать

Полный текст

Аннотация

Цель исследования - построение точных аналитических решений задач статики анизотропных пластин, находящихся под действием циклически симметричных воздействий и лежащих на упругом основании. Для решения поставленных в работе задач используется метод компенсирующих нагрузок, определяются основное и компенсирующие решения. Для нахождения решений используется новый прием, связанный с применением уравнения Нильсена. С помощью метода компенсирующих нагрузок впервые получены точные аналитические решения для задач об анизотропных пластинах, очерченных круговым контуром, лежащих на упругом основании при различных условиях опирания и при действии циклически симметричных нагрузок, распределенных вдоль окружности и по площади кольца. Также рассматривается задача о расчете неограниченной пластины, лежащей на упругом основании и имеющей круговое отверстие. Решения получены в замкнутом виде и выражены в функциях Бесселя.

Полный текст

Введение Пластины постоянной и переменной толщины, а также оболочки вращения находят широкое применение в различных областях: в строительной технике, машиностроении, авиастроении. В частности, фундаментные плиты и междуэтажные перекрытия сооружений, имеющих в плане круговую форму, представляют собой круглую или кольцевую плиту переменной или постоянной толщины. Среди подобных сооружений следует назвать телевизионные и водонапорные башни, дымовые трубы, цилиндрические резервуары. Подобные задачи возникают в машиностроении при расчете турбинных дисков, лопастей турбомашин, дисковых пружин. При защите сооружений от нежелательного уровня вибраций употребляются тарельчатые виброизоляторы, основными конструктивными элементами являются круглые, а иногда и прямоугольные пластины переменной толщины и их системы. Этим вопросам посвящена весьма обширная литература. Назовем, в частности, монографии [1] и [2], в которых отражены многие вопросы теории пластин и оболочек. В [3] рассматривается широкий круг вопросов расчета ортотропных и изотропных пластин различных очертаний при действии сложных нагрузок, применяются точные аналитические методы. В настоящее время существуют мощные программные комплексы, позволяющие изучать работу подобных конструкций. Широко используются численные методы, в частности метод конечных элементов. В [4] исследуются вопросы потери устойчивости ортотропных пластин с различными граничными условиями. В [5] изучаются колебания слоистых пластин, в [6] - вопросы вибрации изотропных и ортотропных прямоугольных пластин линейно-переменной толщины. С помощью численного метода в [7] производится оценка экспериментальных исследований вопросов устойчивости пологих конических оболочек, находящихся под действием внешнего давления. В [8] рассматриваются свободные колебания лопасти переменной толщины при произвольных граничных условиях. Статья [9] посвящена расчету прямоугольной пластины под действием термомеханической нагрузки. С помощью метода конечных элементов в [10] производится оценка работы лопасти ветровой турбины. Оболочки из наноматериала изучаются в [11]. Вибрации слоистых пластин - в [12]. С помощью приближенного аналитического метода декомпозиции уравнений в [13] исследуются некоторые задачи статики, колебаний и устойчивости тонкостенных конструкций. Точные аналитические решения задач расчета комбинированных пластин кусочно-переменной толщины получены в [14]. Однако многие задачи статики и колебаний анизотропных пластин еще не рассмотрены. В настоящей работе изучается весьма актуальная проблема расчета анизотропных круглых пластин, контактирующих с упругим основанием при действии на них циклически симметричных нагрузок. Для решения поставленной задачи используется метод компенсирующих нагрузок (МКН); строится основное решение, удовлетворяющее разрешающему дифференциальному уравнению, и компенсирующее решение, которое в сумме с основным решением удовлетворяет также граничным условиям. Изучается действие разрывных нагрузок и различные условия закрепления. Получено решение задачи об анизотропной неограниченной пластине с круговым отверстием. Разрешающее дифференциальное уравнение задачи о круглой пластине из анизотропного материала, в отличие от случая изотропной пластины, не распадается на два взаимно сопряженных дифференциальных уравнения второго порядка, каждое из которых интегрируется в функциях Бесселя. Ниже используется другой прием, связанный с получением решения в цилиндрических функциях. Вводится в рассмотрение уравнение Нильсена, позволяющее получить интеграл определенного дифференциального уравнения четвертого порядка в цилиндрических функциях. Постановка задачи Рассмотрим задачу о циклически симметричном изгибе круглой анизотропной пластины, лежащей на упругом основании. В этом случае деформации в пластине изменяются по законам или , . Подобный изгиб вызывается самоуравновешенной контурной и поверхностной нагрузкой вида , (1) , (2) где и - интенсивность контурной поперечной силы и контурного изгибающего момента; - интенсивность распределенных по поверхности поперечных сил; , - постоянные величины; - функция только переменной . Циклически симметричный изгиб может возникнуть, когда на фундаментную плиту нагрузка сверху передается через равноотстоящие опоры. Приведем дифференциальное уравнение, описывающее циклически симметричную деформацию круглой анизотропной пластины, лежащей на упругом основании, свойства которого описываются моделью Винклера. Положим, что материал пластины обладает цилиндрической анизотропией и является ортотропным. При внешней нагрузке, определяемой выражением (2), имеем [3] (3) здесь - цилиндрическая жесткость; - коэффициент Пуассона; - коэффициент постели; параметр определяется из соотношений [3], [15]. (4) . Будем разыскивать решение однородного уравнения, соответствующего (3), в виде . (5) Подставляя (5) в (3) и приводя подобные члены, получим (6) В результате ряда преобразований получим следующее разрешающее уравнение: (7) Приведем выражения для изгибающих и крутящих моментов, а также поперечных сил для круглой ортотропной пластины [16]: , (8) , (9) , (10) (11) (12) где . Определим решение дифференциального уравнения (7). Как было указано выше, в отличие от случаев, рассмотренных ранее в литературе, для решения задач о пластинах, сделанных из изотропного материала и лежащих на упругом винклеровском основании при действии неосесимметричных нагрузок [9], разрешающее уравнение не распадается на два сопряженных уравнения второго порядка, каждое из которых интегрируется в функциях Бесселя. В работе используется новый прием: коэффициенты разрешающего дифференциального уравнения (7) будем сопоставлять с коэффициентами уравнения Нильсена, которое имеет вид [16; 17] (13) здесь Как известно, уравнение Нильсена интегрируется в функциях Бесселя [18-20]. В результате сопоставления коэффициентов (7) и (13) имеем (14) отсюда получим следующие значения параметров, при которых решения задачи получаются в цилиндрических функциях: , (15) ; или . Приведем решение однородного уравнения (7), доставляемое функциями Бесселя. В зависимости от значения параметра эти решения можно записать следующим образом: - при (16) - при (17) Далее, используя формулы (8)-(12), определим усилия. Метод компенсирующих нагрузок Основное решение Для решения поставленной задачи применим метод компенсирующих нагрузок. Рассмотрим действие различных циклически симметричных нагрузок и ряд условий закрепления контура. Также изучим задачу о расчете неограниченной анизотропной пластины с круговым отверстием при неосесимметричной нагрузке. Для использования МКН требуется построить основное и компенсирующее решения. Основным назовем такое решение, которое удовлетворяет дифференциальному уравнению и одновременно имеет необходимые особенности, которые соответствуют внешней нагрузке, действующей на пластину. В большинстве случаев основное решение может рассматриваться как решение задачи о неограниченной в плане пластине. Однако основное решение не может удовлетворять граничным условиям, если рассматривается задача о пластине, которая ограничена тем или иным контуром. Компенсирующее решение вводится для того, чтобы вместе с основным решением удовлетворять упомянутым граничным условиям [21-23]. Введем безразмерную координату . При загружении ортотропной пластины сосредоточенной силой основное решение имеет вид . (18) Приведенное выше выражение является основной функцией влияния при . Для отдельных частных задач можно получать основные решения, интегрируя (18). Рассмотрим анизотропную неограниченную пластину, находящуюся под действием сил , равномерно распределенных по окружности, концентрической контуру, с приведенным радиусом . Запишем выражение для прогиба в точке с координатами , пользуясь принципом сложения воздействий: (19) Для вычисления вышеприведенного интеграла (19) воспользуемся формулой сложения функций Бесселя: (20) здесь знак ' символизирует то, что при вводится коэффициент . Указанная формула справедлива при ; если , то в правой части (20) и меняются местами. Положим . Производя интегрирование и отделяя действительную и мнимую части, получим: - при (21) - при (22) При использовании вронскиана уравнения Бесселя можно показать, что полученные выше решения (21), (22) при удовлетворяют условиям сопряжения участков. Решим задачу об анизотропной неограниченной пластине, загруженной по закону , где - заданная функция. Для этого потребуется интегрировать выражения (21) и (22). Пусть на изучаемую пластину действует нагрузка , распределенная по окружности, приведенный радиус которой равен . Разложим указанную нагрузку в ряд Фурье: . Учитывая (21) и (22), получим следующие выражения: - при (23) - при (24) Далее построим основное решение для случая нагрузки, изменяющейся по закону . (25) Для этого будем использовать формулы (23) и (24). Представим нагрузку (25) в виде ряда . (26) В случае, если внешняя нагрузка может быть представлена в виде формулы (27) или , (28) процесс интегрирования заметно упрощается. Приведем формулы дифференцирования функций Бесселя: (29) Положим . Произведем интегрирование; отделим действительную и мнимую части, получим следующие выражения: (30) Рассмотрим действие на неограниченную анизотропную пластину нагрузки, распределенной по площади кругового кольца с внутренним радиусом и наружным , , и изменяющейся по закону (27). Решение можно получить, интегрируя выражения (21) и (22). С этой целью в этих выражениях предварительно заменим нагрузку элементарной нагрузкой . При этом появится множитель , равный . Затем произведем интегрирование указанных выражений в пределах от до . Запишем решение при : (31) - при имеем (32) - при решение принимает вид (33) Приведенные выше формулы (31)-(33) дают выражения прогибов для -го члена ряда (27). Следует отметить имеющий практическое значение случай, когда . Подобное имеет место при расчете фундаментных плит на действие горизонтальных нагрузок. Компенсирующее решение Далее перейдем к рассмотрению анизотропной пластины, лежащей на упругом основании, ограниченной круговым контуром с различными условиями закрепления. Компенсирующее решение представим как результат действия двух компенсирующих нагрузок и . Эти нагрузки действуют по концентрическим окружностям с приведенными радиусами и . Обозначим через приведенный радиус внешнего контура. Компенсирующее решение может интерпретироваться как результат действия некоторой системы сил, которая приложена к пластине бесконечного радиуса. Будем решать задачу для различных условий закрепления контура. Также изучим задачу о неограниченной анизотропной пластине с круговым отверстием. Сначала рассмотрим круглую анизотропную пластину, лежащую на упругом винклеровском основании, защемленную по всему контуру. Представим взятые из основного решения при прогиб и угол наклона нормали к контуру в виде , (34) . (35) Далее определим компенсирующие нагрузки и , которые на контуре пластины в сумме с основным решением должны удовлетворять условиям заделки. Указанное требование приводит к решению системы уравнений (36) (37) Формулы (36) и (37) отражают равенство нулю прогибов и углов поворота на контуре пластины. Для решения этих интегральных уравнений следует разложить функции и в тригонометрические ряды, использовать формулы сложения цилиндрических функций, выполнить интегрирование. Тогда уравнения (36) и (37) перейдут в систему алгебраических уравнений относительно коэффициентов этих рядов. Эта система уравнений здесь не приводится. Компенсирующее решение при можно представить в виде (38) Определим коэффициенты ряда (38), используя выражения (36) и (37); в результате получим следующие системы уравнений: (39) (40) В приведенных выше системах первое и второе уравнения связаны с условиями для прогиба и угла поворота на контуре соответственно. В результате решения указанных систем получим , (41) , (42) , (43) . (44) Для получения компенсирующего решения внесем формулы (41)-(44) в (38). В результате получим (45) Указанным способом можно построить компенсирующие решения для других условий закрепления изучаемых анизотропных пластин на упругом основании. Рассмотрим случай шарнирного опирания по всему контуру. Для того чтобы найти решение поставленной задачи, следует разложить в ряд значения на контуре прогибов и изгибающих моментов, взятых из основного решения. Прогибы будем раскладывать в ряд, пользуясь формулой (34); радиальные изгибающие моменты на контуре могут быть представлены с помощью формулы . (46) Компенсирующее решение в данном случае будет представлено в виде ряда; его коэффициенты определяются из систем уравнений (47) , (48) (49) . (50) Вышеприведенные уравнения получаются в результате приравнивания к нулю суммы прогибов и суммы радиальных изгибающих моментов на контуре для основного и компенсирующего решений соответственно. Введем для краткости следующие обозначения: , (51) , (52) , (53) . (54) Далее внесем полученные значения коэффициентов в ряд (38). Получим следующее выражение для компенсирующего решения: (55) Найдем компенсирующее решение для случая анизотропной круглой пластины со свободным краем. В этом случае, как известно, изгибающий момент и приведенная поперечная сила обращаются на контуре в нуль. Значение приведенной поперечной силы на контуре, взятое из основного решения, разложим в ряд: . (56) Пусть изгибающие моменты даны рядом (46). Уравнения для определения коэффициентов , , , составляются аналогично тому, как это было выполнено в предыдущих примерах. Введем следующие обозначения: , (57) , (58) , (59) . (60) Решая системы уравнений для определения коэффициентов , , , , внесем полученные результаты в ряд (38). Учитывая обозначения (57)-(60), после ряда преобразований получим компенсирующее решение в виде (61) С помощью МКН получим решение задачи об анизотропной неограниченной пластине, лежащей на упругом винклеровском основании, с круговым отверстием при произвольной нагрузке. Искомый результат может быть представлен в виде суммы основного и компенсирующего решений. Основное решение получено выше. Соответствующее компенсирующее решение представим выражением (62) Примем, что начало координат находится в центре отверстия. В компенсирующее решение будут входить только конечные значения функций и и их производных, поскольку точка , в которой названные функции имеют особенности, не входит в изучаемую область. Для определения компенсирующего решения следует в формулах, полученных выше, всюду заменить функции и на и . Приведем компенсирующее решение для случая анизотропной неограниченной пластины со свободным круговым отверстием: (63) Метод компенсирующих нагрузок можно также применить для расчета анизотропной пластины на упругом основании, ограниченной двумя концентрическими окружностями. В этом случае компенсирующее решение следует представить в виде ряда, коэффициенты которого определяются из уравнений, соответствующих существующим граничным условиям. Тогда каждая система будет состоять из четырех уравнений, которые здесь не выписываются. Заключение Впервые получены точные аналитические решения задач об анизотропных пластинах, контактирующих с упругим основанием, при действии на них циклически симметричных нагрузок. Решения выражены в функциях Бесселя. Для интегрирования разрешающего дифференциального уравнения четвертого порядка с переменными коэффициентами применяется новый прием, связанный с использованием уравнения Нильсена. К решению поставленной задачи применяется метод компенсирующих нагрузок. Получено основное решение, удовлетворяющее исходному уравнению. Производится учет действия разрывных нагрузок, распределенных по концентрическим окружностям и по площадям колец. Определяется компенсирующее решение, которое в сумме с основным удовлетворяет граничным условиям. Рассматриваются различные условия закрепления контура. В цилиндрических функциях получено решение об анизотропной неограниченной пластине с круговым отверстием.

×

Об авторах

Елена Борисовна Коренева

Московское высшее общевойсковое командное орденов Ленина и Октябрьской Революции Краснознаменное училище

Автор, ответственный за переписку.
Email: elena.koreneva2010@yandex.ru
SPIN-код: 8804-7930

профессор кафедры общетехнических дисциплин, доктор технических наук

Российская Федерация, 109380, Москва, ул. Головачева, д. 2

Список литературы

  1. Коваленко А.Д. Избранные труды. Киев: Наукова думка, 1976. 703 с.
  2. Коренев Б.Г. Некоторые задачи теории упругости и теплопроводности, решаемые в бесселевых функциях. М.: Физматгиз, 1960. 458 с.
  3. Коренева Е.Б. Аналитические методы расчета пластин переменной толщины и их практические приложения. М.: Изд-во АСВ, 2009. 240 с.
  4. Bank L.C., Yin J. Buckling of orthotropic plates with free and rotationally restrained unloaded edges // Thin-Walled Structures. 1996. Vol. 24. Pp. 83-96.
  5. Chen W.Q., Lüe C.F. 3D free vibration analysis of cross-ply laminated plates with one pair of opposite edges simply supported // Composite Structures. 2005. Vol. 69. Pp. 77-87.
  6. Civalek Ö. Fundamental frequency of isotropic and orthotropic rectangular plates with linearly varying thickness by discrete singular convolution method // Applied Mathematical Modelling. 2009. Vol. 33. Pp. 3825-3835. 10.1016/j.apm.2008.12.019' target='_blank'>https://doi: 10.1016/j.apm.2008.12.019
  7. Karasev Al., Varianychko M., Bessmertnyi Ya., Krasovsky V., Karasev G. Numerical analysis on experimental research on buckling of closed shallow conical shells under external pressure // Journal of Theoretical and Applied Mechanics. 2020. Vol. 58. Issue 1. Pp. 117-126. https://doi.org/10.15632/jtam-pl/115321
  8. Li C., Cheng H. Free vibration analysis of a rotating varying-thickness-twisted blade with arbitrary boundary conditions // Journal of Sound and Vibration. 2020. Vol. 492. Article 115791. https://doi.org/10.1016/j.jsv.2020.115791
  9. Khakpour Komarsofla M., Jedari Salami S., Shakeri M., Khakpour Komarsofla A. Optimization of three-dimensional up to yield bending behaviour using the full layer-wise theory for FGM rectangular plate subjected to thermo-mechanical loads // Compos. Struct. 2021. Vol. 257. 113172. https://doi.org/10.1016/j.compstruct.2020.113172
  10. VanSkike W.P., Hale R.D. Comparative assessment of finite element modelling techniques for wind turbine rotors blades // AIAA Scitech 2020 Forum. Session: Wind Turbine Modeling. 2020. https://doi.org/10.2514/6.2020-0990
  11. Sarafraz A., Sahmani S., Aghdam M.M. Nonlinear primary resonance analysis of nanoshells including vibrational mode interaction based on the surface elasticity theory // Appl. Math. Mech.-Engl. Ed. 2020. Vol. 41. Pp. 233-260. https://doi.org/10.1007/s10483-020-2564-5
  12. Javed S., Al Mukahal F.H.H., El Sayed S.B.A. Geometrical influence on the vibration of layered plates // Hindawi. Shock and Vibration. 2021. Vol. 3. Pp. 1-17. https://doi.org/10.1155/2021/8843358
  13. Koreneva E.B., Grosman V.R. Equation decomposition method for solving of problems of statics, vibration and stability of thin-walled constructions // International Journal for Computational Civil and Structural Engineering. 2020. Vol. 16. No. 2. Pp. 63-70.
  14. Koreneva E.B., Grosman V.R. The problems of computation of combined plates with piecewise variable thickness. Solutions in orthogonal polynomials // International Journal for Computational Civil and Structural Engineering. 2020. Vol. 16. No. 3. Pp. 30-34.
  15. Бурмистров Е.Ф. Симметричная деформация ортотропных оболочек вращения. Саратов.: Изд-во Саратовского университета, 1962. 109 с.
  16. Коренева Е.Б. Метод компенсирующих нагрузок для решения задач об анизотропных средах // International Journal for Computational Civil and Structural Engineering. 2018. Vol. 14. No. 1. Pp. 71-77.
  17. Камке Э. Справочник по обыкновенным дифференциальным уравнениям. М.: Наука, 1965. 703 с.
  18. Abramovitz M., Stigan I.A. Handbook of mathematical functions. 10th ed. National Bureau of Standards, 1972. 820 p.
  19. Коренева Е.Б., Гросман В.Р. Аналитическое решение задачи об изгибе круглой ортотропной пластины переменной толщины, лежащей на упругом основании // Вестник МГСУ. 2011. № 8. С. 156-159.
  20. Гросман В.Р. Собственные колебания круглых ортотропных пластин переменной толщины. Решения в функциях Бесселя // Строительная механика и расчет сооружений. 2012. № 3. С. 52-54.
  21. Коренева Е.Б. Аналитический метод решения задач о колебаниях упругих анизотропных тел // Строительная механика и расчет сооружений. 2018. № 5. С. 47-51.
  22. Коренева Е.Б. Аналитическое моделирование некоторых задач статики и колебаний анизотропных тел // VII International Symposium APCSCE, 1-8 July 2018. Novosibirsk, 2018. P. 478.
  23. Koreneva E.B., Grosman V.R. Forced vibrations of anisotropic elastic solids subjected to an action of complicated loads // International Journal for Computational Civil and Structural Engineering. 2019. Vol. 15. No. 3. Pp. 77-83.

© Коренева Е.Б., 2021

Creative Commons License
Эта статья доступна по лицензии Creative Commons Attribution 4.0 International License.

Данный сайт использует cookie-файлы

Продолжая использовать наш сайт, вы даете согласие на обработку файлов cookie, которые обеспечивают правильную работу сайта.

О куки-файлах