Разрешающие уравнения расчета тонкостенных многострингерных кессонов с учетом геометрической нелинейности Текст научной статьи по специальности «Строительство и архитектура»

УДК 624.157

РАЗРЕШАЮЩИЕ УРАВНЕНИЯ РАСЧЁТА ТОНКОСТЕННЫХ МНОГОСТРИНГЕРНЫХ КЕССОНОВ С УЧЕТОМ ГЕОМЕТРИЧЕСКОЙ НЕЛИНЕЙНОСТИ

© 2008 О. Р. Кузнецов

Саратовский государственный технический университет

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

Призматические оболочки, геометрическая нелинейность, краевая задача, численное интегрирование

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

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

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

поперечного сечения ¥{ (і - номер узла поперечного сечения). Под узлом понимается точка излома контура или точка расположения стрингера. Продольные перемещения точек контура определяются продольными перемещениями ее узлов. Продольные элементы (стрингеры) воспринимают только нормальные напряжения. Касательные и нормальные напряжения по толщине оболочки постоянные.

Статико-геометрическая модель (рис. 1). В соответствии с методом В. З. Вла-

А У

сова расчёта призматических оболочек средней длины компоненты вектора перемещений точек контура в продольном (в направлении оси z) и вдоль контура оболочки (в направлении оси ^) задаются соответственно в виде

U(2 = В^г (Z)фг (S) = Щ (Z)Фг (S); < І > ,

і=1

(1)

V(г, я) = В £ Гк (г)ун (я) = ВГк (г)ун ^); < к >,

к=0

(2)

где ф (я), у к (я), (г = 1,..., п; к = 0,..., т) - задаваемые безразмерные аппроксимирующие функции; и (2), Vh (2) - неизвестные безразмерные функции, подлежащие определению. В дальнейшем для сокращения записи символ “X” опускается, вместо него используется символ < >, который означает суммирование по соответствующему индексу. К системам аппроксимирующих функций

Ф (я), Ук (я) предъявляются следующие требования:

- линейная независимость, желательно ортогональность, то есть выполнение условий:

§ Фг ^ Ф і ^ )* = 0 (і * І ),

к ^у* (sй =0 (к * *);

(3)

= йі/В - относительная длина і-й стороны контура;

8г = 8І/8 - относительная толщина і-й пластинки;

2 = z|l; 5 = ^В - безразмерные координаты точек контура.

Здесь й, В - характерные поперечные размеры контура; 8 - характерная толщина пластинок, составляющих оболочку.

Отметим, что во всех нижеприведённых выражениях величины приведены в безразмерном виде для безразмерных 2 и 5 . Чтобы не загромождать запись, верхняя черточка, которая символизирует безразмерный вид, опускается. В частности, под z и 5 подразумевается 2 и 5 .

В соответствии с выбранной расчётной моделью оболочки отличными от нуля компонентами деформации являются продольные деформации в направлении оси z и деформации сдвига срединной поверхности пластинок, составляющих оболочку. Зависимости между деформациями и перемещениями принимаются в виде

1 Эи(2, 5) ГІі Э

+

+ -

2т7і2

Эи(2, 5) Э2

2

+

ЭУ(2, 5) Э2

2

+

ЭЖ(2, 5) Э2

2

- полнота;

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

Геометрические размеры оболочки задаются следующими безразмерными параметрами:

щ = В/ё - относительная толщина профиля;

щ = I/В - относительное удлинение оболочки;

ц2 = В/8 - параметр, характеризующий степень тонкостенности конструкции;

"12

1

+---

Эи 1 ЭУ

Є 52 = ~Э~ +--------------Э“ +

Э5 Ч1 Э2

Эи Эи ЭУ ЭУ ЭЖ ЭЖ

-------------1--------------1-------------

Э2 Э5 Э2 Э5 Э2 Э5

(4)

Отличными от нуля компонентами тензора напряжений являются (2,5) и

<72 (2,5) . По толщине пластинок, составляющих оболочку, напряжения распределены равномерно. Продольные элементы (стрингеры) воспринимают только нормальные напряжения. Зависимость между напряжениями и деформациями полагаем линейной:

1

®2 ®г (2, 5) Ее2; ^ 12 ^52 (2, 5) 11Ее 12;

1

О = 11Е; 11

2(1 + V)

ЭУ(2, 5)

Э5

= 0, Ж(2,5) = Ук (2)с(s), < к >,

, У,(:)Щх;м

Э2 У к (5 )

можно пренебречь по сравнению с единицей. Тогда имеем

1 Эи 1

Є 2 =-----— +

12

Л1 Э2 2ЛГ

Эи +1_ ЭУ

Э5 Ч1 Э2

ЭУ

Эг

2

+

ЭЖ

Э2

2

Лс2 = ЕЛе 2, Лс12 = 11Е4е

мк = У/(УкУй+с с,); < * >,

или с учётом (3)

(12)

(5)

На этом этапе возможны следующие упрощения.

В выражении (4) учтем, что в случае жёсткого контура поперечного сечения

(6)

где Ж(г, я) - безразмерное перемещение точек контура по нормали к нему; %к (я) - соответствующая аппроксимирующая функция. Затем полагаем, что в выражениях (4), записанных с учетом (6), величинами

М = 1^1+х1). (13)

Напряжённо-деформированное состояние оболочки полностью определяется приращениями обобщённых перемещений на текущем этапе нагружения

иг (2X и'(2X ик (2X

и'к(2) (г = 1,..., п, к = 0,...,т) (14)

и их накопленными на ж-м этапе нагружения значениями

Щ (2),(и'(2))-, у; (2),(УЙ'(2))-.

(15)

(7)

Вывод разрешающих уравнений. Уравнения равновесия оболочки для определения обобщённых перемещений получаются на основе принципа возможных перемещений Лагранжа. На этапе нагружения имеем

8ЛТ - 8ЛЖ = 0.

(16)

, (8)

(9)

Здесь 8ЛТ - работа внешних сил, которая записывается в виде

8ЛТ = В21 § (Я( 2)8ик + Ру (2)8иг й+В2 (Я*8ьк + Р*8иу)

2=1 2, (17)

где

Для линеаризации нелинейных зависимостей используется метод последовательных нагружений [1]. Имеем для приращений величин следующие выражения:

Ле2 = — ифг + -1икМ,; < і >, < к >,

Ч ч2

Л , 1 , , (10) Ле12 = т+—чУк; <г > < к >,

Кк(2) = § Р1( 2, к(^,

Р (2) = § Р2(2, 5Ф (^,

К( 2) = § р*( 2 у к(^,

К(2) = § р* (2Ф(s)ds, (2 =0;1),

(18)

(11)

где иг, ик - приращения искомых обобщенных перемещений и использовано обозначение

р1( 2, я), р2( 2, я) - распределённая по поверхности оболочки нагрузка, действующая соответственно в направлении осей я и 2;

р*(2), р*(2), (2 = 0; 1) - нагрузка, приложенная в торцах оболочки, отнесённая к единице длины контура и действующая соответственно в направлении осей я и 2.

Расчётная модель оболочки задаётся следующим выражением:

1

8ЛЖ = В1§§(Ж2Л8е2 + Л£12Л8е12)2&, (19)

Л (У-Ч2)) • Лу- (2) = (2),

(2 = 0; 1),

(21)

Л(Г --1( 2 )}:

л11 л12 л13 л14 л15

л21 л22 л23 л24 л25

, (22)

л (У ) =

-а‘ 07

„ 7 7

1 І I 1

-фр'Мк) 1йк „ г* + -гз^МкМ")

41 41 I *1/1

<і >, < * >; (у, ,і = 1,...,п) (к,* = 0,...,т); (х = 0;1).

(23)

Элементы матрицы А У (2) имеют вид

г—-1.

Л11 „2 ау, Л12 0у,

41

1 I*

Л13 = -Щ, Л14 = „з \М Фу/’ „1

где

Ш2 = Е8Ле2, Л£12 = (1ХЕ8) Ле12.

После выполнения преобразований получаем краевую задачу для системы линейных обыкновенных дифференциальных уравнений 2(п + т + 1) порядка с переменными коэффициентами относительно неизвестных (14).

Уравнения равновесия записываются в

виде

Л(¥ Ч 2)) • Луя (2) = Л0 (2). (20)

Статические граничные условия имеют вид

Л = -—

Л15 -

„1

Л = -1ФМк), 42 = ^ „і „1

г‘ + £< М‘М*

11

Л23 0к, Л24 „2

„1

Л=5 = Т1 ((М*М^ „1

(24)

Отметим, что элементы матриц ^ и сами

являются матрицами. Индексы у к задают номер строки; і, * - номер столбца матрицы.

Векторы-столбцы искомых приращений обобщенных перемещений на ж-м этапе нагружения имеют вид

где Л(Уг-1(2)), Л *(У г-1(2)) - матрицы коэффициентов, вычисленные с учетом накопленных на ( ж -1) этапе нагружения обобщённых перемещений:

ЛУ-(2) =( <(2X и\(2), иі(2),

и*(2X и*(2))- (і = 1,..., п * = 0,...,т);

ЛУ-( 2) = ( Щ( 2 X иі(2), и *(2))

(25)

(і = 1,..., п, * = 0,...,т); 2 = 0;1.

Векторы-столбцы правых частей имеют вид

Л@ (2 ) = ~Е Р (х), я, (г))Т

Е

(у = 1,..., п, к = 0,...,т);

Л0 *(2) = е18 (Р*(х), (х ))Т

(у = 1,..., п, * = 0,...,т); (2 = 0;1).

(26)

Метод решения. Полученная краевая задача может быть решена численным интег-

0

Ь

С

рированием путем сведения к серии задач Коши. Численное интегрирование полученной краевой задачи сопровождается явлениями численной неустойчивости. Для преодоления этих проблем использован метод С. К. Годунова [2].

Рассмотрим порядок получения решения для нагрузки

0 (2 ) = N • Л0( 2 );

(0*(2) = N• Л0*(2), (2 = 0; 1)),

(27)

где ^-количество ступеней нагружения. На первом этапе нагружения (при ж =1) в коэффициентах (23), (24) полагаем Мк = М'к = 0,

т. е. рассматривается геометрически линейная задача, а также полагаются равными нулю обобщённые перемещения (15). В результате получается линейная краевая задача (20), (21) с постоянными коэффициентами. Вид этих разрешающих уравнений точно совпадает с разрешающими уравнениями В. З. Власова для оболочек такого типа. В результате решения краевой задачи получаются значения приращений обобщённых перемещений (14). Затем обобщённые перемещения (15) полагаются равными полученным приращениям обобщённых перемещений (14). Вычисляются коэффициенты (23), (24) с учетом полученных значений обобщённых перемещений (15) и, соответственно, отличными от

нуля значениями Мк, М'к. Эти коэффициенты становятся переменными, зависящими от 2. Затем решается линейная краевая задача (20), (21), но теперь с переменными коэффициентами. Находятся приращения обобщённых перемещений (14) на втором этапе нагружения ( ж =2). К значениям обобщённых перемещений (15) прибавляются полученные ( ж =2) приращения обобщённых перемещений (14). Полученные таким образом значения обобщённых перемещений (15) дают решение задачи для нагрузки:

0 (2 ) = 2 • Л0 (2);

(0 *(2) = 2Л0 *(2), (2 = 0; 1))

(28)

ной Q(z), Q*(z) (выражения (27)). Основные трудности вычислительного плана возникают при пересчёте на каждом этапе нагружения коэффициентов, связанных с нелинейностью. Они на этом этапе таковы, что ставят под сомнение возможность получения правильных численных результатов. Поэтому важным этапом при решении задач такого типа является обоснованное упрощение выражений для коэффициентов, связанных с нелинейностью. Первый этап упрощений реализован выражениями (6), (7). Второй этап проведём на основе анализа выражений для коэффициентов (23), (24). Этот анализ показывает, что с учетом вида выражения (13) для Мк можно пренебречь членами, содержащими их квадраты, и производной по 2 от этого

квадрата, а также произведениями Мк, Мк на продольные перемещения и производные по ъ от продольных перемещений как членами более высокой малости по сравнению с другими. Тогда имеем для коэффициентов (22) матрицы Л:

Л21 0, Л22 йк,

„1

Л23 = а Л24 = ~~2 ГІ , Л25 = 0

„1

(30)

и для матрицы Л* граничных условий (23):

Л (У --1( 2)) =

0к

¥к -г*

„1

, (31)

где коэффициенты матрицы вычисляются следующим образом:

а

с* (2) =< у>, йк(2) =< уф; > , К(2) =< >,

(32)

Описанная процедура выполняется N раз до достижения значений нагрузки, рав-

< М*ф7 >= §Мф7 8 йк,

< (М* )фг >= § М'*ф7 8 .

Для записи коэффициентов символ < > означает интеграл по контуру. Напомним, что

р О) = ЭР(5)/ Э5.

Таким образом, имеем краевую задачу (20), (21) с матрицами коэффициентов (22), (24), (31) с учётом (30) и выражениями (32), (33).

Продолжая анализ выражений для коэффициентов краевой задачи, отметим, что нелинейность учитывается коэффициентами (33). В случае жёсткого контура поперечного сечения эти коэффициенты точно равны нулю

для всех р7 (5), которые задают продольные

перемещения точек контура как жёсткого диска (при неподвижном полюсе) плюс депла-национные перемещения. Эти коэффициенты отличны от нуля только для функции

р0(5), которая задает поступательное перемещение точек контура в продольном направлении. Поясним сказанное на следующем примере. Рассмотрим прямой однозамкнутый кессон с прямоугольным контуром поперечного сечения, который закреплён, например, как показано на рис. 2. Кессон нагружен некоторой изгибающей и крутящей нагрузкой.

Система аппроксимирующих функций в этом случае имеет вид [3], приведённый на

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

Функции у0 (5) , С0 (5) - перемещение точек контура в его плоскости соответственно при

кручении относительно оси ъ; у1( 5), %х( 5) -при изгибе в плоскости ^. Вид функций р1, р3 может уточняться, но их характер симметрии относительно осей Xи У остаётся неизменным.

С учетом выражений (13) коэффициенты (33) записываются в виде

< Мкр7 >= V'(< р7у1 > + < р7х1 >) = ^ (34)

так как

< ру2 >= IрУ8ё5=0 < Руй >= ^ (35)

как интеграл от произведения кососимметричных функций р7 на симметричные

у 1; . Поэтому для возможности учёта гео-

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

Функция р0( 5) задает поступательные

перемещения точек контура в продольном направлении, т. е. растяжимость продольной

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

Рис. 2 133

X

д)

г)

I У

1 1 1 Г+1

1 1 м |У1(5)| 1+ X —► С1(5) 1 1 - ) 1

X --►

е)

Рис. 3

ж)

носительно этих осей нагрузки приводит к и(г, 5) = Ви0(г)у0(5), (36)

несимметричному относительно осей х, у

распределению напряжённо-деформирован- где р1(5), у0(5), р0(5) приведены на рис. 3,а; ного состояния точек контура оболочки.

3,б и 4,а, соответственно. Уравнения равно-

Пример расчета. Рассмотрим кессон, г г г г ■> весия (20) принимают вид

изображённый на рис. 2, который находится под действием сосредоточенного крутящего

момента (рис. 5). Геометрические размеры — а!Х + —2а\щ- 1Ь%0 -Щщ + —3 < М°р > и"0-

оболочки задаются следующими параметрами:

В = ё2; ё = й{; Л = 2; 4 = 4; В = ё2 ё 2= 1; 1.

- к ' V

1

с7- л < Мр >

7и, 0 7,

Л

7 >

и'0 = 0 (7 = 0; 1),

/

ё1 = -; 81 = 82

Л

(37)

Л. ёV +-^ ё Щ + -Л г 0и' = 0 Приращения перемещений в соответ- л°° Ч°' Ч2°°

ствии с (1), (2) задаются в виде щ(^ 5) = В(щ0(г) • Р0(5) + щ1(р(5)),

Статические граничные условия (21) при г = 1 принимают вид

аV + — а7 Щ + < М0р > и0 = 0 (7 = 0; 1),

Ч 7 Л1 7 7

7 /0 7/1 11 0 / 1 тл*

11ё0Щ0 + 11ё0Щ1 + ' и0 = ^0 .

Ч Е8

(38)

Геометрические граничные условия имеют вид

и0 = щ = и0 = 0 (г = 0) .

(39)

Для заданной схемы нагружения правая часть (26) уравнений равновесия (20) и граничных условий (21) в соответствии с (18) имеет вид

Р>( г) = Р( г) = К«(г) = #(1) = Р*(1) = 0; Я0(1) = Мр.

(40)

Получена краевая задача (37), (38), (39) для системы линейных дифференциальных уравнений с переменными коэффициентами шестого порядка относительно шести неизвестных:

и0(гХ Щ0(^Х и'(гХ Щ1(гХи0(гХ и0(г). (41)

М

Рис. 5

Коэффициенты этой краевой задачи вычисляются с помощью выражений (32),

(33). Алгоритм решения описан в разделе «метод решения».

Значение правой части в статических граничных условиях (38) на этапе нагружения принимается равным

2Л1(1 + П) Я„* = 0,0028.

Е8 0

Это соответствует значению максимального депланационого напряжения в угловой точке контура

Лаг = — Ей', • 0,1 = 42,85 МПа.

Л1

Нагружение ведётся до четырнадцатой ступени нагружения # = 14, что соответствует

а1 = 600 МПа.

На рис. 6 приведены графики изменения вдоль оболочки:

- приращений продольных перемещений, связанных с депланацией контура поперечного сечения щ1(г) • р1 (т.1) =щ1(г) 0,1 -кривая 1 (линейная составляющая). Запись п = 1-14 на рис. 6 означает, что рассматриваемая кривая одинакова на всех ступенях нагружения;

- приращений продольных перемещений, связанных с растяжимостью продольной

оси оболочки и0(г) • р0(5) = и0(г) - кривая 2, для седьмой (п = 7) и четырнадцатой (п = 14) ступени нагружения (нелинейная составляющая).

I

На рис. 7 приведены графики изменения вдоль оболочки:

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

щ(г) •р1(т.1) = щ(г) • 0,1 - кривая 1 (линейная составляющая);

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

и0 (г) • (р0 (^) = и0 (г) - кривая 2 для четырнадцатой (п = 14) ступени нагружения (нелинейная составляющая).

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

Анализ графиков на рис. 6, 7 показывает, что учёт геометрической нелинейности приводит к увеличению продольных перемещений (рис. 6, кривая 2) и появлению нормальных напряжений (рис. 7, кривая 2) в сечении (г = 1), где они отсутствовали. Перемещения в сечении г = 1 увеличиваются на 30 %, дополнительные напряжения в сечении

2 = 1 равны 15 МПа.

Поэтому в данном случае следует конструктивно предусматривать возможность для поступательного перемещения сечений оболочки на величину, равную до 30 % от значений максимальных депланационных перемещений, и учитывать появление нормальных напряжений в сечениях (в данном случае г = 1).

Начиная с уровня напряжений, равных примерно 500 МПа, перемещения и0(г) и

и'0(г), связанные с учетом геометрической нелинейности, начинают резко расти. Например, приращения этих величин на одном четырнадцатом ( ж =14) этапе нагружения равны сумме приращений на первых пяти этапах.

Выводы

1. Учет геометрической нелинейности эквивалентен отказу от гипотезы о нерастя-жимости продольной оси оболочки.

2. При внешней нагрузке, которая вызывает в оболочке нормальные напряжения порядка 200-250 МПа, гипотеза о нерастяжи-мости продольной оси выполняется с высокой степенью точности. Вклад геометрической нелинейности для нормальных напряжений составляет до 1,2 % от максимальных; для перемещений - от 1 % до 8 % от максимальных.

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

Рис. 6 Рис. 7

обходимо учитывать в практике проектирования соответствующих конструкций.

4. При внешней нагрузке, которая вызывает в оболочке напряжения, равные примерно 500 МПа, перемещения и напряжения, связанные с учётом геометрической нелинейности, т.е. и0(г) и и0(г), начинают резко расти.

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

Библиографический список

1. Петров В. В. Метод последовательных нагружений в нелинейной теории пластин и оболочек. - Саратов: Изд-во Сарат. гос. ун-та, 1975.

2. Власов В. З. Избранные труды. - М.: Наука, 1964.

3. Образцов И. Ф. Вариационные методы расчёта тонкостенных авиационных конструкций. - М.: Машиностроение, 1966.

4. Образцов И. Ф., Онанов Г. Г. Строительная механика скошенных тонкостенных систем. - М.: Машиностроение, 1973.

5. Петров В. В., Кузнецов О. Р. Некоторые нелинейные задачи расчёта прямых замкнутых призматических оболочек. (Моногра-фия)/Сарат. гос.техн. ун-т. - Саратов, 1996. -85 с. - Деп. в ВИНИТИ, 31.10.96, № 3197-В96 ДЕП.

References

1. V. V. Petrov. Successive loading method in the non-linear theory of plates and shells. -Saratov: Saratov State University Publishing House, 1975.

2. V. Z. Vlasov. Selected works. - Moscow: Nauka (Science), 1964.

3. I. F. Obraztsov. Variational methods of thin-walled aircraft structure design. - Moscow: Mashinostroyeniye (Machine building), 1966.

4. I. F. Obraztsov, G. G. Onanov. Structural

mechanics of skewed thin-walled systems. -Moscow: Mashinostroyeniye (Machine

building), 1973.

5. V. V. Petrov, O. R. Kuznetsov. Some non-linear problems of straight closed prismatic shell design. (Monograph)/Saratov State Technical University - Saratov, 1996. - 85 pp. -Dep. at All-Russian Institute of Scientific and Technical Information, 31.10.96, No. 3197-B 96 Dep.

SOLVING EQUATIONS OF THIN-WALLED MULTISTRINGER TORSION BOX DESIGN IN VIEW OF GEOMETRICAL NON-LINEARITY

© 2008 O. R. Kuznetsov

Saratov State Technical University

A procedure of thin-walled multistringer torsion box design in view of geometrical non-linearity is proposed. Initial non-linear relations are linearized by the successive loading method. Solving equations are obtained by using the Lagrange principle; these are presented in the form of a boundary problem for a system of ordinary linear differential equations.

Prismatic shells, geometrical non-linearity, boundary problem, numerical integration

Информация об авторе Кузнецов Олег Рафаэльевич, доцент, кандидат технических наук, кафедра «Механика деформируемого твёрдого тела» Саратовского государственного технического университета. Область научных интересов: расчёт прямых замкнутых призматических оболочек с учётом нелинейных факторов.

Kuznetsov, Oleg Raphaelievitch, Associate Professor, Candidate of Technical Science, Department of deformable solid body mechanics, Saratov State Technical University, e-mail [email protected]. Area of research: Design of straight closed prismatic shells in view of non-linear factors.