Влияние жесткости шпангоута на механику деформирования цилиндрической оболочки Текст научной статьи по специальности «Механика и машиностроение»

ВИНОГРАДОВ Юрий Иванович

(МГТУ им. Н.Э. Баумана)

VINOGRADOV Yuriy Ivanovich

(Moscow, Russian Federation, Bauman Moscow State Technical University)

УДК 539.3

Влияние жесткости шпангоута на механику деформирования цилиндрической оболочки

Ю.И. Виноградов

При воздействии нагрузки на оболочку через шпангоут она распределяется по ее кольцевой поверхности. В зависимости от жесткостных характеристик шпангоута напряженно-деформированное состояние оболочки изменяется. Для оценки прочности подкрепленной шпангоутом оболочки важно знать это изменение. Для этого необходимо построить математическую модель механики совместного деформирования оболочки и шпангоута и решить численно задачу прочности подкрепленной оболочки. Решение данной задачи предполагает, что шпангоут идеализируется упругой кольцевой линией, наделенной естественными жесткостными характеристиками, которая совмещается со срединной поверхностью оболочки. Математическая модель механики деформирования шпангоута (система уравнений равновесия элемента кольцевой линии, геометрические и физические соотношения) описывается системой из четырех линейных обыкновенных дифференциальных уравнений для определения составляющих перемещения ее точек и угла закручивания. Построенная математическая модель позволила дать количественную оценку имеющимся моделям, которые использовались при анализе конструкторских решений. Расчеты в широком диапазоне изменения параметров оболочки и шпангоута показали, что использование упрощенных математических моделей шпангоута приводит к недопустимому, даже качественному, изменению результатов, которые определяют прочность подкрепленной шпангоутом оболочки.

Ключевые слова: цилиндрическая оболочка, шпангоут, математическая модель механики деформирования шпангоута, жесткость шпангоута.

The influence of the frame stiffness on the deformation of a cylindrical shell

Y.I. Vinogradov

The load acting on a shell through its frame is distributed circumferentially over the shell surface. The stress-strain state of this shell depends on the frame stiffness. To assess the strength of the stiffened shell, it is necessary to develop a mathematical model of joint deformation of the shell and frame and solve this problem numerically. In this paper, the frame is represented by an elastic circle with appropriate stiffness characteristics, which lies in the middle surface of the shell. The mathematical model of deformation of the frame includes the set of equilibrium equations for a circular element and geometrical and physical relationships. In fact, it is a system of four linear ordinary differential equations

МАШИНОСТРОЕНИ

solvable for the displacement components of a frame point and the angle of twist. The proposed mathematical model made it possible to estimate available models used for the analysis of design solutions. The calculations in a wide range of parameters ofshells and frames showed that simplified mathematical models of a frame do not provide an adequate accuracy of the results, which is unacceptable for evaluating strength of shells stiffened by frames.

Keywords: cylindrical shell, frame, mathematical model of deformation of a frame, frame stiffness.

Избежать локальной передачи сил и моментов на оболочку в составе конструкции, например, рамного соединения ступеней ракет, крепления рулей летательных и подводных аппаратов и т. п., невозможно. Для снижения напряжений в оболочке нагрузку от сил и моментов распределяют по оболочке. Типовым конструкторским решением такой задачи является подкрепление оболочки шпангоутом, на который действуют силы и моменты. При этом силы и моменты передаются распределенными. Для рационального подкрепления оболочки необходимо подобрать жесткость шпангоута, отвечающую противоречивым требованиям: с одной стороны, обеспечить прочность оболочки, с другой, — требование по допустимой величине веса оболочки, подкрепленной шпангоутом.

Известно, что распределение сил и моментов, приходящих на оболочку, зависит от его жесткости. Чем больше жесткость шпангоута, тем равномернее передается нагрузка на оболочку. Однако при этом растет вес подкрепленной оболочки, иногда до недопустимых значений.

Таким образом, в конструкторской практике актуальной является задача выбора жестко-стных характеристик шпангоута для обеспечения весового совершенства конструкции. Эту задачу обычно решают вариантными расчетами на прочность оболочки и подкрепляющего ее шпангоута. При этом, как правило, параметры оболочки задаются, а место размещения и же-

сткостные характеристики шпангоута подбираются вариантными вычислениями. Проблема решается путем математического моделирования механики деформирования оболочки и шпангоута, а также моделирования подкрепления оболочки шпангоутом и количественного анализа на ЭВМ.

В предлагаемой работе выбрана математическая модель механики деформирования замкнутой в окружном направлении цилиндрической оболочки постоянной толщины из изотропного материала [1, 2], построена математическая модель механики деформирования шпангоута, исследовано влияние распространенных в расчетной практике допущений на жесткостные характеристики шпангоута: жесткость на изгиб при деформировании в его плоскости, из его плоскости, жесткости на растяжение и кручение.

Математическая модель механики деформирования шпангоута. Шпангоут моделируется упругой кольцевой линий, наделенной его же-сткостными характеристиками.

Положительные направления осей координат, перемещений, усилий и моментов для упругой кольцевой линии шпангоута показаны на рис. 1.

Уравнения равновесия элемента шпангоута имеют следующий вид:

Рис. 1. Силовые факторы, действующие в сечении шпангоута, интенсивности внешних усилий и моментов

= 0; +Qz - 4у*=0; йМх

О.

йф

^у -9-Я = ^-ая=0; (1)

йМу йМг

_ф- + ^ - т>Л = 0; + ОхЯ - Му = 0,

где Qy, Qz, Qx — соответственно усилия растяжения кольцевой линии и перерезывающие усилия, направленные по радиусу и перпендикулярно плоскости кольцевой линии; Му, М 7, Мх — моменты в сечении шпангоута, векторы которых направлены соответственно по касательной у к кольцевой линии, по нормали I и бинормали х; 9у, , 9Х — интенсивности внешних усилий, направленные соответственно по касательной у, нормали I и бинормали х к упругой кольцевой линии шпангоута; ту, т7, тх — интенсивности внешних моментов, изгибающих упругую кольцевую линию шпангоута в ее плоскости, из ее плоскости и закручивающий ее момент.

Исключая перерезывающие усилия Qx, Qz из уравнений равновесия (1), получаем

1 йМу й2 М йОу йМ л п

у г - ^ =0; -йфг + ях- = 0;

Я йф Яйф2

й2 М йМу

—Г ^ - = 0; ^ + М - туЯ = 0. (2)

Яйф2

Здесь Я — радиус упругой кольцевой линии.

Физические соотношения имеют следующий вид:

ау = ЕЕв, Мх = Е1Х Х1; Му = 0/ кр 0, М 7 = X 2

(3)

где ЕЕ, Е1Х, Е17, 01 кр — жесткости шпанго-

ута на растяжение (сжатие), на изгиб в плоскости кольцевой линии и из ее плоскости и на закручивание соответственно; в, х1, X2, © — линейная, изгибные и угловая деформации.

Геометрические соотношения имеют вид

в=

1 Я

Я2

йv \ 1

йф + 7 Х1 = Я2

й 2ы \

й 2w + йv йф

\

2 - Яр ; 0 = —2 йрф И) Я2

1

йрф

/

йи я йв | йф йф)

(4)

Здесь в — угол поворота сечения шпангоута при закручивании.

Исключая деформации в физических соотношениях (3) с помощью геометрических соотношений (4), подставляя полученные выражения для усилия растяжения (сжатия) и моментов в уравнения равновесия, получаем систему линейных обыкновенных дифференциальных уравнений в перемещениях (математическую модель механики деформирования шпангоута):

01 к

Я3 01

й2и й2в\

+ Я

йф

/л 2

йф2

Я3

й4и й2в\

- Я

йф йф

= яЛ

Я3

й2и й2в\

+ Я

уйф2 йф2 )

ш.

+

й2и

Я3

\

йф2

- Яв

ту

ЕТХ ( А4 й ™ й 3 V] ЕЕ

Я3 1йф4 йфъ) Я {

йv \

ЕЕ

Я3

й V й2 VI - +

йф йф )

) ЕЕ 'йм> й 2 V)

+- — + йф2)

) я \йф

*+йф)=(5)

= 9, Я,

где и, V, ^ — составляющие перемещения точек упругой кольцевой линии; 1Х, I у, / , Е — моменты инерции и площадь сечения шпангоута; Е — модуль упругости материала шпангоута.

Как и для оболочек, которые подкрепляют шпангоуты, решение системы уравнений (5) ищем в форме тригонометрических рядов:

и= 2и» ^Ы; V=2 V» ^(пф);

п=1 »=1

™=2™» со8(»Ф); в=2в» со§(»Ф). (6)

п =1

п=1

При этом внешняя нагрузка, например, в виде сосредоточенных сил и моментов, раскладывается в ряды:

Р = 0 +29- ^М;

п=1

Му = ту0 + 2ту» со8(»ф) и т. д. (7)

п =1

В итоге, система дифференциальных уравнений (5) приводится к системе алгебраических уравнений:

X

2

01

я 01,

кр

3 (ип + яв„ )п

2 Е] (п \ + Явп )п2 = дшЯ;

я

я

Е]

^ (ип + Явп )п2 - ^ (п \ + Явп) = туп;

Е1Г 3 ЕЕ

—Г (п^п + Уп )п3 - я ^ + пУп) = д^Я; (8)

Е1Х 2 ЕЕ

- ~яг(™п + V)п - я + п^„)п = ду„я-

Разложения (6) и (7) справедливы, если внешняя нагрузка на шпангоут симметрична относительно начала координаты ф.

В результате решения системы алгебраических уравнений для каждого номера п гармоники определяем параметры ип, vn, , вп тригонометрических рядов (6) искомых величин и, V, м, р.

Усилия и моменты в шпангоуте находим с помощью тригонометрических рядов:

Оу = ^уп ео8(пф); Мх = ^ИХп ео8(пф);

п=1 п=1

МУ =^Муп втМ; =2 Мп ео8(пф). (9)

п =1 п =1

Значения функциональных коэффициентов О-уп, Мхп, Муп, Мп этих рядов определяем по формулам:

ЕЕ

Е1х

Оуп = ^т К + "Vп); мп = —2^ (шп + Vn )п;

~Уп я Х " п ' '' ' п /7 хп я2 0/

Муп =

я

Е1

кр

2 (ип- яРп)п;

(10)

Мп (п \ + явп),

которые получены с помощью физических соотношений (3), геометрических соотношений (4) и представления перемещений в виде рядов (6).

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

Взаимодействие оболочки и шпангоута моделируем действием сил и моментов, которые возникают в деформированной оболочке, на упругую кольцевую линию. Эти силы и моменты имеют смысл внешней нагрузки.

Систему алгебраических уравнений (8) представим в матричной форме:

1-V2

(11 0 0 (14 ип

0 ( 22 ( 23 0 Vn

0 ( 32 ( 33 0 п

(41 0 0 (44 яр,

ЛЕ

1п

ЛS,

ЕН

1-V2

ло

ЕН

12(1- V2)

1п

ЛМ

ЕН

12(1- V2)

1п

ЕН

(11)

В сокращенной записи

Gp = f, (12)

где f — столбец параметров сил и моментов, действующих на упругую кольцевую линию со стороны сопряженных с ней оболочек; p — столбец параметров перемещений и угла закручивания; 0 — матрица жесткости упругой кольцевой линии, элементами которой являются следующие параметры:

(11 =-п2(а1п2 + Р1); (14 =-п2(а1 + Р1);

(41 =■ ( 32 =■ ( 22 = "

-п2 (а2 + р2); (44 =-(а2 + п 2р2);

-п(а3п + Р3); (33 =-(а3п + Р3); -п2 (а4 + р4); (23 = -п(а4 п2 + Р4);

•2 Е1Х 1- V2

я3 ЕН

я3 ЕН

Р1 =

0] кр 1

V

2

а

я3 ЕН

а

2

-; Р4 =

ЕЕ 1-V2 я ЕН

а

2

= П. в = РЪ в = Ёъ с 2 =.

2 ; в2 2 ; Из 2; с

Н2

0~ 12я

Влияние жесткостных характеристик шпангоута на состояние оболочки. Для определения влияния жесткостных характеристик шпангоута на состояние оболочки выполнены параметрические исследования для консольно закрепленной цилиндрической оболочки, усиленной шпангоутом на свободном краю. При этом использовалась математическая модель механики деформирования цилиндрической оболочки В.З. Власова [1, 2]. На шпангоут действует сосредоточенная радиальная сила (рис. 2). Задачу решали методом сопряжения

п

а1 =

а4 =

а 2 =

интервалов [3—5]. Основное внимание уделено влиянию жесткостных характеристик шпангоута на внутренние силовые факторы, которые возникают в оболочке в окрестности сопряжения ее со шпангоутом. Дана количественная оценка широко распространенных допущений об отсутствии жесткости шпангоута при изгибе из его плоскости и об отсутствии жесткости при его закручивании.

Рис. 2. Расчетная схема, внутренние силовые факторы в сечениях оболочки и шпангоута

Исследования выполняли для оболочки

1 -х R ™

с параметрами — = 3 и — = 50 и шпангоута R h

с квадратным поперечным сечением 0^ х 0,5h; 5h х5h; 20h х 20h; 30h х30h; 50h х50h при следующих условиях:

шпангоут наделен всеми жесткостными характеристиками (рис. 4, кривая а);

отсутствует жесткость шпангоута на изгиб из его плоскости, EJZ — 0 (рис. 4, кривая б);

отсутствует жесткость при закручивании шпангоута, GJy — 0 (рис. 4, кривая в);

шпангоут обладает изгибной жесткостью только в своей плоскости, GJy — 0, EJZ — 0 (рис. 4, кривая г).

Моменты инерции Jx, Jy, JZ для шпангоута

Ь 4

определяли по формулам Jx — J7 —

12;

Jy — кхЬ4, где Ь — 0,5^ 5h и т. д., k — 0,1406 [6].

Изменения искомых величин в безразмер-

— пR — п

ном виде 71,2 — Ти2 — и Ml2 — Ми — вдоль об-

разующей оболочки, проходящей через точку приложения силы Р1 в малой ее окрестности Ь

0 < £ < —, и шпангоута сечением 30h х 30h по-5R

казаны на рис. 3—6 сплошными линиями.

Рис. 3. Внутренние силовые факторы

Рис. 4. Изменение усилий в оболочке вдоль образующей при ф — 0

Результаты вычислений показали, что значения искомых величин для оболочки, подкрепленной шпангоутом с размерами 0^ х0,5^ совпадают (что и следовало ожидать) с соответствующими для неподкрепленной оболочки. Они приведены на рис. 5 и 6 пунктирными и штрихпунктирными линиями соответственно.

Из параметрического исследования, некоторые результаты которого представлены на рис. 3—6, следует, что необоснованные часто используемые допущения при построении математических моделей механики деформирования шпангоутов для расчета подкрепленных

МАШИНОСТРОЕНИ

Ml-10

Mi

6

/Q fa

1 ( s*

I

Mr 1 0

1

Рис. 5. Изменение продольного момента в оболочке вдоль образующей при ф = 0

Рис. 6. Изменение окружного момента в оболочке вдоль образующей при ф = 0

оболочек могут приводить к результатам даже качественно отличающимся от тех, которые находят без каких-либо допущений (см. рис. 5).

Выводы

Полученная математическая модель механики деформирования шпангоута позволила количественно оценить используемые для конструкторских решений упрощенные математиче-с к и е м о д е л и м е х а н и к и д е ф о р м и р о в а н и я шпангоута. Показано, что упрощенные математические модели механики деформирования шпангоута приводят к результатам, характери-

зующим прочность оболочки с недопустимой погрешностью.

Литература

1. Власов В.З. Общая теория оболочек и ее приложение в технике. М.: Гостехиздат, 1949. 784 с.

2. Власов В.З. Избранные труды. М.: Изд-во АН СССР, 1962. Т. 1. 528 с.

3. Виноградов Ю.И, Клюев Ю.И. Тонкостенные осесим-метричные конструкции из композиционных материалов (метод решения задач статики и динамики)//Механика композиционных материалов и конструкций. ИПриМ РАН, 1984. Т. 4. № 1. С. 57-72.

4. Vinogradov Yu, Bakulin V Numerical method investigation of shells under the impact of concentrated load// Matematicke metody v technickych vedach mathematical in engineering/ Karlovy Vary.: Czechoslovak Scince and Technikal Socity Charles University Prague — Faculty of Mathematies and Phaysics Slovak society for Mechanics subsidiary of Slovak Academy of Sciences, Central Research Institute SKODA Concern, Concern Enterprise Pzen, 1986. Vol. 1. P. 411—412.

5. Виноградов Ю.И. Методы вычислений и построение алгоритмов решения краевых задач строительной механики // Докл. АН СССР, 1988. Т. 298. № 2. С. 308—313.

6. Варвак П.М, Рябов А.Ф. Справочник по теории упругости. Киев: Будевельник, 1971. 419 с.

References

1. Vlasov V.Z. Obshchaia teoriia obolochek i ee prilozhenie v tekhnike [General theory of shells and its application in engineering]. Moscow, Gostekhizdat publ., 1949, 784 p.

2. Vlasov V.Z. Izbrannye trudy [Selected Works]. Moscow, AN SSSRpubl., 1962, vol. 1. 528 p.

3. Vinogradov Iu.I., Kliuev Iu.I. Tonkostennye osesimmetrichnye konstruktsii iz kompozitsionnykh materialov (metod resheniia zadach statiki i dinamiki) [Thin-walled axisymmetric structures made of composite materials (method of solving problems of statics and dynamics)]. Mekhanika kompozitsionnykh materialov i konstruktsii [Journal on Composite Mechanics and Design]. Institute of Applied Mechanics of Russian Academy of Sciences publ., 1984, vol. 4, no. 1, pp. 57—72.

4. Vinogradov Yu., Bakulin V. Numerical method investigation of shells under the impact of concentrated load. Matematicke metody v technickych vedach mathematical in engineering. Karlovy Vary, Czechoslovak Scince and Technikal Socity Charles University Prague — Faculty of Mathematies and Phaysics Slovak society for Mechanics subsidiary of Slovak Academy of Sciences, Central Research Institute SKODA Concern, Concern Enterprise Pzen, 1986, vol. 1, pp. 411—412.

5. Vinogradov Iu.I. Metody vychislenii ipostroenie algoritmov resheniia kraevykh zadach stroitel'noi mekhaniki [Calculation methods and the construction of algorithms for solving boundary value problems of structural mechanics]. Report of the Academy of Sciences of the USSR publ., 1988, vol. 298, no. 2, pp. 308—313.

6. Varvak P.M., Riabov A.F. Spravochnikpo teorii uprugosti [Handbook of the theory of elasticity]. Kiev, Budevel'nik publ., 1971. 419 p.

Статья поступила в редакцию 26.04.2013

Информация об авторе

ВИНОГРАДОВ Юрий Иванович (Москва) — доктор физико-математических наук, профессор кафедры «Аэрокосмические системы». МГТУ им. Н.Э. Баумана (105005, Москва, Российская Федерация, 2-я Бауманская ул., д. 5, стр. 1, e-mail: vino.yuri@rambler.ru).

Information about the author

VINOGRADOV Yuriy Ivanovich (Moscow) — Dr. Sc. (Phys. Math.), Professor of «Aerospace Systems» Department. Bauman Moscow State Technical University (BMSTU, building 1, 2-nd Baumanskaya str., 5, 105005, Moscow, Russian Federation, e-mail: vino.yuri@rambler.ru).