Научная статья на тему 'АНАЛИЗ УСТОЙЧИВОСТИ ВЫЧИСЛИТЕЛЬНОГО АЛГОРИТМА К ИЗМЕНЕНИЮ ГЕОМЕТРИЧЕСКИХ ПАРАМЕТРОВ ЦИЛИНДРИЧЕСКИХ ОБОЛОЧЕЧНЫХ КОНСТРУКЦИЙ'

АНАЛИЗ УСТОЙЧИВОСТИ ВЫЧИСЛИТЕЛЬНОГО АЛГОРИТМА К ИЗМЕНЕНИЮ ГЕОМЕТРИЧЕСКИХ ПАРАМЕТРОВ ЦИЛИНДРИЧЕСКИХ ОБОЛОЧЕЧНЫХ КОНСТРУКЦИЙ Текст научной статьи по специальности «Строительство и архитектура»

CC BY
52
11
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ОБОЛОЧЕЧНЫЕ КОНСТРУКЦИИ / ЦИЛИНДРИЧЕСКИЕ ОБОЛОЧКИ / МОДЕЛЬ ТИМОШЕНКО / УСТОЙЧИВОСТЬ ЧИСЛЕННЫХ МЕТОДОВ / УСТОЙЧИВОСТЬ СТРОИТЕЛЬНЫХ КОНСТРУКЦИЙ / ВЫЧИСЛИТЕЛЬНЫЙ ЭКСПЕРИМЕНТ

Аннотация научной статьи по строительству и архитектуре, автор научной работы — Бакусов П. А., Семенов А. А.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

ANALYSIS OF THE STABILITY OF THE COMPUTATIONAL ALGORITHM TO A CHANGE IN THE GEOMETRIC PARAMETERS OF CYLINDRICAL SHELL STRUCTURES

This study deals with testing sustainability of a computational algorithm to a change in geometric parameters of cylindrical shell structures. A change in geometry implies the replacement of one type of a cylindrical shell (elliptic, hyperbolic, parabolic) with another so that the quantitative change (the difference in elevations) in the area under consideration is minimal. On the one hand, this test allows to assessing the correctness of the algorithm itself and is relevant for algorithms that use both numerical methods and symbolic calculations. On the other hand, it allows to evaluating the possibility of simplifying calculations by approximating a complex surface with a simpler one both in understanding the surface definition itself and in expressing its basic characteristics such as Lame coefficients and main curvatures. A mathematical model of deformations of shell structures based on the hypotheses of Timoshenko (Mindlin - Reisner) are used in the work. The model takes into account transverse shifts, geometric nonlinearity and orthotropy of the material, and its written in the form of a functional of the total potential strain energy. The calculation algorithm is built on the basis of the Ritz method to reduce the variational problem of the minimum functional to the solution of a system of nonlinear algebraic equations, and on the method of continuing the solution with the best parameter for its solution. All calculations were carried out in dimensionless parameters. Three types of cylindrical panels are calculated, and critical loads of buckling and deflection fields at subcritical and supercritical moments are obtained. It is shown that for the considered class of problems the previously proposed mathematical model and computational algorithm are resistant to changes in the geometry of the structure.

Текст научной работы на тему «АНАЛИЗ УСТОЙЧИВОСТИ ВЫЧИСЛИТЕЛЬНОГО АЛГОРИТМА К ИЗМЕНЕНИЮ ГЕОМЕТРИЧЕСКИХ ПАРАМЕТРОВ ЦИЛИНДРИЧЕСКИХ ОБОЛОЧЕЧНЫХ КОНСТРУКЦИЙ»

Бакусов П.А., Семенов А.А. Анализ устойчивости вычислительного алгоритма к изменению геометрических параметров цилиндрических оболочечных конструкций // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2021. - № 1. С. 12-21. DOI: 10.15593/peim.mech/2021.1.02

Bakusov P.A., Semenov A.A. Analysis of the stability of the computational algorithm to a change in the geometric parameters of cylindrical shell structures. PNRPU Mechanics Bulletin, 2021, no. 1, pp. 12-21. DOI: 10.15593/perm.mech/2021.1.02

ВЕСТНИК ПНИПУ. МЕХАНИКА № 1,2021 PNRPU MECHANICS BULLETIN

https://ered.pstu.ru/index.php/mechanics/index

Б01: 10.15593/регш.шееЬ/2021.1.02 УДК 539.3

АНАЛИЗ УСТОЙЧИВОСТИ ВЫЧИСЛИТЕЛЬНОГО АЛГОРИТМА К ИЗМЕНЕНИЮ ГЕОМЕТРИЧЕСКИХ ПАРАМЕТРОВ ЦИЛИНДРИЧЕСКИХ ОБОЛОЧЕЧНЫХ КОНСТРУКЦИЙ

П.А. Бакусов, А.А. Семенов

Санкт-Петербургский государственный архитектурно-строительный университет, Санкт-Петербург, Россия

АННОТАЦИЯ

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

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

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

©ПНИПУ

О СТАТЬЕ

Получена: 22 июня 2020 г. Принята: 30 декабря 2020 г. Опубликована: 15 апреля 2021 г.

Ключевые слова: оболочечные конструкции, цилиндрические оболочки, модель Тимошенко, устойчивость численных методов, устойчивость строительных конструкций, вычислительный эксперимент.

© Бакусов Павел Анатольевич - асп., e-mail: [email protected], : 0000-0003-1215-1183 Семенов Алексей Александрович - к.т.н., доц., зав. каф., e-mail: [email protected], ¡D: 0000-0001 -9490-7364

Pavel A. Bakusov - PhD student, e-mail: [email protected], ¡D: 0000-0003-1215-1183 Alexey A. Semenov - CSc in Technical Sciences., Associate Professor, Head of Department, e-mail: [email protected], ¡D: 0000-0001 -9490-7364

Эта статья доступна в соответствии с условиями лицензии Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)

This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)

ANALYSIS OF THE STABILITY OF THE COMPUTATIONAL ALGORITHM TO A CHANGE IN THE GEOMETRIC PARAMETERS OF CYLINDRICAL SHELL STRUCTURES

P.A. Bakusov, A.A. Semenov

Saint Petersburg State University of Architecture and Civil Engineering, Saint Petersburg, Russian Federation

ABSTRACT

This study deals with testing sustainability of a computational algorithm to a change in geometric parameters of cylindrical shell structures. A change in geometry implies the replacement of one type of a cylindrical shell (elliptic, hyperbolic, parabolic) with another so that the quantitative change (the difference in elevations) in the area under consideration is minimal. On the one hand, this test allows to assessing the correctness of the algorithm itself and is relevant for algorithms that use both numerical methods and symbolic calculations. On the other hand, it allows to evaluating the possibility of simplifying calculations by approximating a complex surface with a simpler one both in understanding the surface definition itself and in expressing its basic characteristics such as Lame coefficients and main curvatures.

A mathematical model of deformations of shell structures based on the hypotheses of Timo-shenko (Mindlin - Reisner) are used in the work. The model takes into account transverse shifts, geometric nonlinearity and orthotropy of the material, and its written in the form of a functional of the total potential strain energy.

The calculation algorithm is built on the basis of the Ritz method to reduce the variational problem of the minimum functional to the solution of a system of nonlinear algebraic equations, and on the method of continuing the solution with the best parameter for its solution. All calculations were carried out in dimensionless parameters.

Three types of cylindrical panels are calculated, and critical loads of buckling and deflection fields at subcritical and supercritical moments are obtained. It is shown that for the considered class of problems the previously proposed mathematical model and computational algorithm are resistant to changes in the geometry of the structure.

©PNRPU

ARTICLE INFO

Received: 22 June 2020 Accepted: 30 December 2020 Published: 15 April 2021

Keywords:

shell structures, cylindrical shells, Tymoshenko's model, stability of numerical methods, buckling of building structures, numerical experiment.

Введение

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

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

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

и особенности применения которых также могут оказывать влияние на получаемые результаты [16, 17].

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

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

конструкций покрытий и перекрытий1 (относящееся к СП 52-117-20082), вместо пологой поверхности (считающейся сектором сферической) рассматривается сектор кругового параболоида, т.е. происходит замена (аппроксимация) более сложного выражения геометрии более простым.

Данная работа направлена на оценку устойчивости вычислительного алгоритма по расчету значений критических нагрузок, соответствующих потери локальной и глобальной устойчивости оболочечных конструкций к геометрии самих оболочек (здесь под потерей устойчивости понимается «прохлопывание» конструкции, когда малому изменению нагрузки соответствует большое изменение ее прогиба).

Для оценки влияния вида геометрии на значения критических нагрузок и оценки устойчивости соотношений модели к виду геометрии конструкции предлагается провести следующее исследование. В качестве малого измерения входного параметра рассмотрим три вида цилиндров (цилиндрических поверхностей): эллиптический, гиперболический и параболический. Этот выбор был сделан не случайно. Во-первых, для данных поверхностей достаточно легко подобрать такие коэффициенты, что качественно они будут отличаться мало (в нашем случае наибольшее отклонение сечений составило 1,92 • 10-2 мм). Во-вторых, специфика геометрии (нулевая гауссова кривизна) позволяет несколько упростить и ускорить расчет. В-третьих, цилиндрические поверхности являются одними из самых распространенных, из-за чего исследование их поведения на данный момент крайне актуально [19-29].

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

1. Теория и методы

1.1. Исследуемые поверхности

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

Пусть задана прямоугольная система координат 0xyz, положительные направления осей которой описываются левой тройкой векторов (ось z направлена вверх). Тогда цилиндры будут представлять собой параллельный перенос кривой, заданной в плоскости y0z вдоль оси 0x.

1 Пособие по проектированию железобетонных пространственных конструкций покрытий и перекрытий (к СП 52-117-2008). М.: ОАО «НИЦ "Строительство"», 2010. 159 с.

2 СП 52-117-2008*. Железобетонные пространственные

конструкции покрытий и перекрытий. Методы расчета и кон-

струирование. М.: ОАО «НИЦ "Строительство"», 2010. 143 с.

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

где Щ и Щ - радиусы (или проекции) вдоль осей 0у и 0z соответственно.

Рис. 1. Схемы сечений цилиндров: а - эллиптического; б - гиперболического; в - параболического

Fig. 1. Schemes of cross-section for some types of cylinders: а - elliptic; б - hyperbolic; в - parabolic

Для расчета оболочечных конструкций необходимо знать коэффициенты Ляме и главные кривизны. Для их поиска необходимо обратиться к теории дифференциальной геометрии (например, [30]).

Для эллиптического цилиндра коэффициенты Ляме

A = i, B =

а главные кривизны

r4 + r2 - ri2 y2

- у2 ri2

к = 0, к = -

R1R2

(r4 +(r2 -r2)y2)

3/2 •

Второй рассматриваемой поверхностью является гиперболический цилиндр, поперечное сечение которого изображено на рис. 1, б. Явная форма задания этой поверхности имеет вид

Здесь Щ - это расстояние от фокуса гиперболоида F до асимптоты, а Щ - расстояние от вершины гиперболоида до оси 0у.

Коэффициенты Ляме и главные кривизны гиперболического цилиндра:

A = 1, B =

ri4 + r2 + ri2 y2

(r - y2 ri2

к = 0 , к = -

R14 R2

(ri4 +(r22 + r12) у2 )

32 •

Третья рассматриваемая поверхность - параболический цилиндр. Его поперечное сечение приведено на рис. 1, в.

В явном виде данная поверхность может быть представлена как

z = —

J_ 2R

где ^ - это расстояние от фокуса параболы ^ до директрисы Д или двойное расстояние от фокуса до оси 0у.

Коэффициенты Ляме и главные кривизны этой поверхности имеют следующее выражение:

Л =1, * = +1

к = о,

r''Rf+1

3/2 •

1.2. Математическая модель

Будем рассматривать классическую теорию оболочек нулевой гауссовой кривизны, учитывающую поперечные сдвиги. Срединная поверхность оболочки принимается за координатную. Используемой моделью является модель деформирования оболочки типа Тимошенко (Миндлина - Рейснера), тогда перемещения в слое на расстоянии г (толщина оболочки И) от срединной поверхности составят:

и2 = и+, V2 = V+, W2 = W,

х У

где и = и (х, у) , V = V (х, у), W = W (х, у) представляют собой неизвестные функции перемещений вдоль осей 0х, 0у и 0г соответственно, а х,у),

^ = ^ (х, у) - неизвестные функции углов поворота нормали в плоскостях х0, и у0г .

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

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

1 си 1 за гл7 1 2

е =--+-V--к W + —9,2,

х А дх AB дУ Х 2

1 3v 1 db 1 2

еv =--+-U--к W+ -92,

B ду AB дх у 2 2

1 dV 1 dU 1 TdA 1 dB „ „

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

у ху =--+----U---V — + 9192

ху А дх B ду AB ду AB дх 12

У* = kf (z)[Tx-91 ], уyz = kf (z)[Ty -92] ,

1 1 dW

9 =-|--+ к и

1 ' А дх х

I 1 dW т/

92 =-|--+ кyV

2 . B dy у ,

а функции изменения кривизн и кручения

1 1 йА1Т(

Х1 =--х +--Г

1 дГ у

1 dB,

А дх AB ду

-Г^^-^у > Х2 =

B ду AB дх

Х12

А дх B су AB дх у AB ду х

где гх, еу - деформации удлинения вдоль осей 0х, 0у срединной поверхности; у , уж, у- деформации сдвига в плоскостях х0у, х0,, у0г соответственно; / (г) - функция, характеризующая распределение напряжений и х по толщине обшивки; к = 5/6.

Для слоя обшивки, расположенного на расстоянии г от срединной поверхности, деформации составят:

К =ех + 2Хр , еу =еу + 2X2 , уху =Уху + 22Х12.

Физические соотношения (связь напряжений и деформаций) для ортотропного материала при линейно -упругом деформировании запишутся следующим образом [31]:

СТх = -Т-^1-[ех + ЦыЕу + 2 (Х1 + ^21X2 )] ,

1 "М12^21

СТу -[8у +^128х + 2 (Х2 + М-12Х1 )]

1 — М12М2

Хху = О12 [Уху + 22Х12 ] , Хх, = О13 Ух:, , Ху, = О23 Уу2 .

Здесь Е , Е , Цг , , О12 , О13 , О23 - механические характеристики материала. В случае изотропного материала необходимо принять Е = Е = Е, ц2 = ц21 = ц,

° = ° = °2з = О.

Выражения для усилий и моментов находятся путем интегрирования напряжений по переменной г в пределах от —к/2 до к/2.

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

1

2

1

F = j dz = h ,

-к/2

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

к/ 2

S = j zdz = 0 ,

-h/ 2

момент инерции сечения, приходящийся на единицу длины сечения,

h/2

J = j z2dz = — j 19

12

Поэтому составляющие усилий и моментов будут иметь вид [31]

N = ^ к (в, + у ) , = сг к (Еу + ^Вх ) ;

h3 h3

= Gi— (x + Ц21Х2) > My = G2 — (%2+M12X1);

h3

Ny = G hy ^, M^ = G122 - X12;

Q = ^Gi3h (Y,-01), Q = kG2h (Ту-02),

где

1 -Ml2^21

E

G = - E

1 - ^12^21

Функционал полной потенциальной энергии деформации оболочечной конструкции можно представить в виде [33]

I »'Г 1

Е = -ГГ N в + N в + -(N + N )у +

Í 2JJ XX у у ху ух ) I ху

»0 I-

+МхХ- + МуХ2 + (Мху + Мх )Х-2 +

+Qx (^х -01)+QУ К -02)-] лвахду. (1)

1.3. Алгоритм расчета

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

Согласно этому алгоритму, к функционалу применяется метод Ритца для сведения вариационной задачи к системе нелинейных алгебраических уравнений. Для этого искомые функций представляются в виде

U ■■ = U ( x, y ) = •jn 4n = X TUkiXY, k=1 l=1

V 11 V y 4N 4N = X Y/UX in, k=1 l=1

W ■■ X II 4n 4n =X TWiX in, k=1 l=1

Y x = = Y x(x, y ) •Jn -Jn = TTPSklX4y4 , k=1 l=1

Y y = Yy(x, y ) = ■Jn -Jn =TLpn»xY k=1 l=1

где С/и - Р^ - искомые коэффициенты разложения. Подставив функции (2) в функционал (1), находим производные по неизвестным числовым параметрам ии - Р^. Таким образом, получаем систему нелинейных алгебраических уравнений.

Все расчеты проводятся в безразмерных параметрах, однако все рассуждения приводятся в размерном виде, чтобы не загромождать изложение. Подробно используемые безразмерные параметры и их обоснование приведены в работе [34].

2. Расчеты

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

Зафиксируем два параметра: длину проекции сечения на ось 0у (ширина покрытия Ь) и на ось 0z (высота покрытия d) (рис. 2). Пусть первая Ь = 6 м, а вторая d = 1 м. Таким образом, было зафиксировано три точки: две у основания покрытия и одна на его высоте.

Рис. 2. Схема размеров цилиндров Fig. 2. Scheme of sizes of cylinders

Для начала определим значение параметра R у параболического цилиндра. Как известно, через три точки, не лежащих на одной прямой, можно единственным образом задать параболу. Из этого получаем, что параметр R = 9/2. Используя формулу для расчета длины дуги [30, с. 53]

-А/2

1 - N1+ИJ

находим

l -

IfR*

y2 dy и 6,41870 м.

Теперь, фиксируя значение длины дуги, определим параметры Я и Я2 эллиптического цилиндра. Для этого выразим один параметр через другой, используя зафиксированную высоту

\(

R

1 М

r2

R2 - d ^ R -

bR

2yjd (2R - d) '

Методом подбора был определен параметр Я2 так, чтобы длина дуги эллипса была равна длине дуги параболы с точностью до пятого знака после запятой (до 0,01 мм): Я, = 38883/161, R = 12961,

i-i.

1 +

y2 R2

(R12 - y2) R

-dy -

1+

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

y2129612

388832 1612

y

388832 1612

-dy и 6,41870.

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

' (V 2)21

1 + R

r12

R2 - R - d ^ R -

bR

2^d ( 2R + d ) '

В итоге, параметры R - 51336/185 , R -17112, а длина дуги

i-i

1 +

y2 R22

(R12 + y2) R

-dy -

-i

1 +

y 171122

513362 1852

-+у

513362 1852

dy и 6,41870.

Все полученные геометрические характеристики рассматриваемых оболочек приведены в табл. 1.

Рассматривался вариант, когда граница оболочек закреплена шарнирно-неподвижно. Материал конструкции - сталь (Е = Е = 2,1-105 МПа, |2=|21 = 0,3).

На рис. 3 представлены рассматриваемые оболо-чечные конструкции. Как видно, количественно они практически неотличимы друг от друга.

Таблица 1

Геометрические характеристики оболочек

Table 1

Geometric characteristics of shells

Характеристика Тип цилиндра

эллиптический гиперболический параболический

A 1 1 1

B i R +(R2 - R) y2 R4 +(R2 + R) y2 #+1

(R2 - y2) R12 ^ (R12 + y2) R12

К , м-1 0 0 0

К y, м-1 RR« RlR2 1

(R +(R2 -R)y2f (R + (R + R) y2 f 4 R+1j

R1, м 38 883/161 51 336/185 9/2

R2, м 12 961 17 112 -

h, м 0,03

а, м 6

b, м 6

а1, м 0

а2, м 6

b1, м -3

b2, м 3

в

Рис. 3. Рассматриваемые виды оболочек: а - эллиптический цилиндр; б - гиперболический цилиндр; в - параболический цилиндр

Fig. 3. The types of shells considered: а - elliptic cylinder; б - hyperbolic cylinder; в - parabolic cylinder

Их отличия представлены в виде графиков абсолютного отклонения сечений друг от друга на рис. 4. Как видно, наибольшее отклонение соответствует координатам y и+2,1213.

В ходе вычислительного эксперимента были получены графики нагрузка q - прогиб W, представленные на рис. 5, а также сами значения нагрузки (табл. 2), при которой оболочка теряет устойчивость (на графике этой нагрузке соответствует точка первого экстремума). Поля прогибов, отложенные от поверхности оболочки, представлены на рис. 6-8 (на рис. 7, 8 представлены проекции полей рис. 6 на плоскости x0y, x0z иy0z).

Как видно по графикам рис. 5, характер деформирования у всех трех видов оболочек один и тот же.

2

b

-3

3

-3

Рис. 4. Графики абсолютных отклонений высотных отметок сечений цилиндров различных видов

Fig. 4. Graphs of absolute deviations of elevation of cross-sections for different types of cylinders

Таблица 2

Значения критической нагрузки qcr цилиндрических оболочек

Table 2

Values of the critical load qcr of cylindrical shells

Тип цилиндра Значение нагрузки

Эллиптический цилиндр 0,742655 кПа

Гиперболический цилиндр 0,742323 кПа

Параболический цилиндр 0,742466 кПа

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

Рис. 5. Графики нагрузка q - прогиб W для трех видов цилиндров

Fig. 5. Load q versus deflection Wfor three types of cylinders

ir

f

а б

Рис. 6. Поля прогибов W до (а) и после (б) критической нагрузки qcr, отложенные на поверхности (для наглядности значения прогибов увеличены в 3 раза)

Fig. 6. The fields of deflections W before (а) and after (b) the critical load qcr, deposited on the surface (for clarity, the values of the deflections are tripled)

Рис. 7. Проекции прогибов W до критической нагрузки (для наглядности значения прогибов увеличены в 3 раза)

Fig. 7. The projections of deflections W before the critical load (for clarity, the values of deflections are tripled)

Рис. 8. Проекции прогибов W после критической нагрузки (для наглядности значения прогибов увеличены в 3 раза)

Fig. 8. The projections of deflections W after the critical load (for clarity, the values of the deflections are tripled)

Выводы

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

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

1. Totaro G. Flexural, torsional, and axial global stiffness properties of anisogrid lattice conical shells in composite material // Composite Structures. - 2016. - Vol. 153. - P. 738-745. D01:10.1016/j. compstruct.2016.06.072

2. Efimtsov B.M., Lazarev L.A. Forced vibrations of plates and cylindrical shells with regular orthogonal system of stiffeners // Journal of Sound and Vibration. - 2009. - Vol. 327, № 1-2. -P. 41-54. D0I:10.1016/j.jsv.2009.05.021

3. Sun Y., Qiu Y., Wu Y. Modeling of Wind Pressure Spectra on Spherical Domes // International Journal of Space Structures. - 2013. -Vol. 28, no. 2. - P. 87-100. D0I:10.1260/0266-3511.28.2.87

4. Кривошапко С.Н. О возможностях оболочечных сооружений в современной архитектуре и строительстве // Строительная механика инженерных конструкций и сооружений. - 2013. - № 1 (1). - С. 51-56.

5. Tomás A., Martí P. Shape and size optimisation of concrete shells // Engineering Structures. - 2010. - Vol. 32, no. 6. - P. 16501658. D0I:10.1016/j.engstruct.2010.02.013.

6. Watts G., Singha M.K., Pradyumna S. Nonlinear bending and snap-through instability analyses of conical shell panels using element free Galerkin method // Thin-Walled Structures. - 2018. -Vol. 122. - P. 452-462. D0I:10.1016/j.tws.2017.10.027

7. Frikha A., Dammak F. Geometrically non-linear static analysis of functionally graded material shells with a discrete double directors shell element // Computer Methods in Applied Mechanics and Engineering. - 2017. - Vol. 315. - P. 1-24. D0I:10.1016/j.cma.2016.10.017

8. Timergaliev S.N., Uglov A.N., Kharasova L.S. Solvability of geometrically nonlinear boundary-value problems for shallow shells of Timoshenko type with pivotally supported edges // Russian Mathematics. - 2015. - Vol. 59, no. 5. - P. 41-51. D0I:10.3103/S1066369X15050060

9. Maksimyuk V.A., Storozhuk E.A., Chernyshenko I.S. Variational finite-difference methods in linear and nonlinear problems of the deformation of metallic and composite shells (review) // International Applied Mechanics. - 2012. - Vol. 48, no. 6. - P. 613-687. D0I:10.1007/s10778-012-0544-8

10. Awrejcewicz J., Krysko V.A. 3-D theory versus 2-D approximate theory of free orthotropic (isotropic) plate and shell vibrations, part 2: Numerical algorithms and analysis // Journal of Sound and Vibration. - 1999. - Vol. 226, no. 5. - P. 831-871. D0I: 10.1006/jsvi.1999.2320

11. Timergaliev S.N. A method of integral equations in nonlinear boundary-value problems for flat shells of the Timoshenko type with free edges // Russ Math. - 2017. - Vol. 61, no. 4. - P. 49-64. D0I: 10.3103/S1066369X17040089

12. Savula Y.H., Jarmai K., Mukha I.S. Analysis of shells reinforced by massive stiffening ribs // International Applied

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

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

Mechanics. - 2008. - Vol. 44, no. 11. - P. 1309-1318. DOI: 10.1007/s10778-009-0137-3

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

13. Dung D.V., Chan D.Q. Analytical investigation on mechanical buckling of FGM truncated conical shells reinforced by orthogonal stiffeners based on FSDT // Composite Structures. - 2017. -Vol. 159. - P. 827-841. DOI: 10.1016/j.compstruct.2016.10.006

14. Patel S.N., Datta P.K., Sheikh A.H. Buckling and dynamic instability analysis of stiffened shell panels // Thin-Walled Structures. - 2006. - Vol. 44, no. 3. - P. 321-333. DOI: 10.1016/j.tws.2006.03.004.

15. Buckling of intermediate ring supported cylindrical shells under axial compression / Y. Xiang [et al.] // Thin-Walled Structures. - 2005. - Vol. 43, no. 3. - P. 427-443. DOI: 10.1016/j.tws.2004.07.019

16. Kumar Y. The Rayleigh-Ritz method for linear dynamic, static and buckling behavior of beams, shells and plates: A literature review // Journal of Vibration and Control. - 2017. -P. 107754631769472. DOI:10.1177/1077546317694724

17. Trushin S. Numerical algorithm for solving of nonlinear problems of structural mechanics based on the continuation method in combination with the dynamic relaxation method // MATEC Web of Conferences / ed. V. Andreev. - 2016. - Vol. 86. -P. 01006. DOI:10.1051 /matecconf/20168601006

18. Дикович В.В. Пологие прямоугольные в плане оболочки вращения. - Л.; М.: Госстройиздат, 1960. - 143 с.

19. Тупикова Е.М. Выбор оптимальной оболочки покрытия на квадратном плане в виде поверхности переноса // Строительная механика инженерных конструкций и сооружений. - 2019. - Vol. 15, № 5. - С. 367-373. DOI: 10.22363/1815-5235-2019-15-5-367-373

20. Giloulbe M., Qbaily J. Геометрическое моделирование и линейный статический расчет тонких оболочек в форме цилиндроидов // Строительная механика инженерных конструкций и сооружений. - 2018. - Vol. 14, № 6. - С. 502508. DOI: 10.22363/1815-5235-2018-14-6-502-508

21. Фирсанов В.В. Локальное напряженно-деформированное состояние цилиндрической оболочки на основе трёхмерных уравнений теории упругости // Строительная механика инженерных конструкций и сооружений. - 2013. -№ 2. - С. 10-19.

22. Li D., Qing G., Liu Y. A layerwise/solid-element method for the composite stiffened laminated cylindrical shell structures // Composite Structures. - 2013. - Vol. 98. - P. 215-227. DOI: 10.1016/j.compstruct.2012.11.013

23. Sobhaniaragh B., Nejati M., Mansur W.J. Buckling modelling of ring and stringer stiffened cylindrical shells aggregated by graded CNTs // Composites Part B: Engineering. - 2017. -Vol. 124. - P. 120-133. DOI: 10.1016/j.compositesb.2017.05.045

24. Ng T.Y., Lam K.Y. Dynamic stability analysis of cross-ply laminated cylindrical shells using different thin shell theories // Acta Mechanica. - 1999. - Vol. 134, no. 3-4. - P. 147-167. DOI: 10.1007/BF01312653

25. Martins J.P., Simoes da Silva L., Silvestre N. Energy-based analytical model to predict the elastic critical behaviour of curved panels // Journal of Constructional Steel Research. -2016. - Vol. 127. - P. 165-175. DOI: 10.1016/j.jcsr.2016.07.029

26. Sengupta J., Ghosh A., Chakravorty D. Progressive Failure Analysis of Laminated Composite Cylindrical Shell Roofs // Journal of Failure Analysis and Prevention. - 2015. -Vol. 15, no. 3. - P. 390-400. DOI: 10.1007/s11668-015-9951-6

27. Ahmed M.K. Simplified equations and solutions for the free vibration of an orthotopic oval cylindrical shell with variable thickness // Mathematical Methods in the Applied Sciences. - 2011. - Vol. 34. - P. 1789-1800. DOI: 10.1002/mma.1493

28. The Dynamic Stability of Cylindrical Shells with Variable Geometric and Mechanical Parameters / A.H. Sofiyev, F. Birinci, Al.H. Sofiyev, E. Yusufoglu // Multidiscipline Modeling in Materials and Structures. - 2006. - Vol. 2, no. 3. -P. 345-354. DOI: 10.1163/157361106777641369

References

1. Totaro G. Flexural, torsional, and axial global stiffness properties of anisogrid lattice conical shells in composite material Composite Structures, 2016, vol. 153, pp. 738-745. doi:10.1016/j.compstruct.2016.06.072.

2. Efimtsov B.M., Lazarev L.A. Forced vibrations of plates and cylindrical shells with regular orthogonal system of stiffeners Journal of Sound and Vibration, 2009, vol. 327, no. 1-2, pp. 4154. DOI:10.1016/j.jsv.2009.05.021.

3. Sun Y., Qiu Y., Wu Y. Modeling of Wind Pressure Spectra on Spherical Domes International Journal of Space Structures, 2013, vol. 28, no. 2. pp. 87-100. DOI:10.1260/0266-3511.28.2.87.

4. Krivoshapko S.N. O vozmozhnostjah obolochechnyh sooruzhenij v sovremennoj arhitekture i stroitel'stve [On the possibilities of shell structures in modern architecture and construction]. Stroitel'naja mehanika inzhenernyh konstrukcij i sooruzhenij, 2013, no. 1 (1), pp. 51-56.

5. Tomás A., Martí P. Shape and size optimisation of concrete shells Engineering Structures, 2010, vol. 32, no. 6. pp. 1650-1658. DOI:10.1016/j.engstruct.2010.02.013.

6. Watts G., Singha M.K., Pradyumna S. Nonlinear bending and snap-through instability analyses of conical shell panels using element free Galerkin method Thin-Walled Structures, 2018, vol. 122, pp. 452-462. DOI:10.1016/j.tws.2017.10.027.

7. Frikha A., Dammak F. Geometrically non-linear static analysis of functionally graded material shells with a discrete double directors shell element Computer Methods in Applied Mechanics and Engineering, 2017, vol. 315, pp. 1-24. DOI:10.1016/j.cma.2016.10.017.

8. Timergaliev S.N., Uglov A.N., Kharasova L.S. Solvability of geometrically nonlinear boundary-value problems for shallow shells of Timoshenko type with pivotally supported edges Russian Mathematics, 2015, vol. 59, no. 5. pp. 41-51. DOI:10.3103/S1066369X15050060.

9. Maksimyuk V.A., Storozhuk E.A., Chernyshenko I.S. Variational finite-difference methods in linear and nonlinear problems of the deformation of metallic and composite shells (review) International Applied Mechanics, 2012, vol. 48, no. 6. pp. 613-687. DOI:10.1007/s10778-012-0544-8.

29. Vibration analysis of ring-stiffened conical-cylindrical-spherical shells based on a modified variational approach / Y. Qu, S. Wu, Y. Chen, H. Hua // International Journal of Mechanical Sciences. - 2013. - Vol. 69. - P. 72-84. DOI: 10.1016/j.ijmecsci.2013.01.026

30. Погорелов А.В. Дифференциальная геометрия. - М.: Наука, 1969. - 176 с.

31. Karpov V.V., Semenov A.A. Refined model of stiffened shells // International Journal of Solids and Structures. - 2020. -Vol. 199. - P. 43-56. D0I:10.1016/j.ijsolstr.2020.03.019

32. Шалашилин В.И., Кузнецов Е.Б. Метод продолжения решения по параметру и наилучшая параметризация (в прикладной математике и механике). - М.: Эдиториал УРСС, 1999. - 222 с.

33. Semenov A.A., Leonov S.S. The Continuous Method of Solution Continuation with Respect to the Best Parameter in the Calculation of Shell Structures // Uch. Zap. Kazan. Univ. Ser. Fiz.-Mat. Nauki. - 2019. - Vol. 161, no. 2. - P. 230-249. DOI: 10.26907/2541-7746.2019.2.230-249

34. Karpov V., Semenov A. Dimensionless parameters in the theory of reinforced shells // PNRPU Mechanics Bulletin. -2015. - No. 3. - P. 74-94. DOI: 10.15593/perm.mech/2015.3.07

10. Awrejcewicz J., Krysko V.A. 3-D theory versus 2-D approximate theory of free orthotopic (isotropic) plate and shell vibrations, part 2: Numerical algorithms and analysis Journal of Sound and Vibration, 1999, vol. 226, no. 5. pp. 831-871. DOI: 10.1006/jsvi.1999.2320.

11. Timergaliev S.N. A method of integral equations in nonlinear boundary-value problems for flat shells of the Timoshenko type with free edges Russ Math, 2017, vol. 61, no. 4. pp. 49-64. DOI: 10.3103/S1066369X17040089.

12. Savula Y.H., Jarmai K., Mukha I.S. Analysis of shells reinforced by massive stiffening ribs International Applied Mechanics, 2008, Vol. 44, no. 11, pp. 1309-1318. DOI: 10.1007/s10778-009-0137-3.

13. Dung D.V., Chan D.Q. Analytical investigation on mechanical buckling of FGM truncated conical shells reinforced by orthogonal stiffeners based on FSDT Composite Structures, 2017, vol. 159, pp. 827-841. DOI: 10.1016/j.compstruct.2016.10.006.

14. Patel S.N., Datta P.K., Sheikh A.H. Buckling and dynamic instability analysis of stiffened shell panels Thin-Walled Structures, 2006, vol. 44, no. 3, pp. 321-333. DOI: 10.1016/j.tws.2006.03.004.

15. Xiang Y. et al. Buckling of intermediate ring supported cylindrical shells under axial compression Thin-Walled Structures, 2005, vol. 43, no. 3, pp. 427-443. DOI: 10.1016/j.tws.2004.07.019.

16. Kumar Y. The Rayleigh-Ritz method for linear dynamic, static and buckling behavior of beams, shells and plates: A literature review Journal of Vibration and Control, 2017, pp. 107754631769472. DOI:10.1177/1077546317694724.

17. Trushin S. Numerical algorithm for solving of nonlinear problems of structural mechanics based on the continuation method in combination with the dynamic relaxation method MATEC Web of Conferences / ed. Andreev V, 2016, vol. 86, pp. 01006. DOI:10.1051 /matecconf/20168601006.

18. Dikovich V.V. Pologie prjamougol'nye v plane obolochki vrashhenija [Gentle rectangular in terms of shells of revolution]. Leningrad-Moscow, Gosstrojizdat, 1960. 143 p.

19. Tupikova E.M. Vybor optimal'noj obolochki pokrytija na kvadratnom plane v vide poverhnosti perenosa [The choice of the

optimal coating shell on a square plan in the form of a transfer surface], Stroitel'naja mehanika inzhenernyh konstrukcij i sooruzhenij, 2019, vol, 15, no, 5, pp, 367-373, DOI: 10,22363/1815-5235-201915-5-367-373,

20, Giloulbé M., Qbaily J. Geometricheskoe modelirovanie i linejnyj staticheskij raschet tonkih obolochek v forme cilindroidov Stroitel'naja mehanika inzhenernyh konstrukcij i sooruzhenij, 2018, vol, 14, No, 6, pp, 502-508, DOI: 10,22363/1815-52352018-14-6-502-508,

21, Firsanov V,V, Lokal'noe naprjazhenno-deformirovannoe sostojanie cilindricheskoj obolochki na osnove trjohmernyh uravnenij teorii uprugosti [Local stress-strain state of a cylindrical shell based on three-dimensional equations of the theory of elasticity], Stroitel'naja mehanika inzhenernyh konstrukcij i sooruzhenij, 2013, no, 2, pp, 10-19,

22, Li D,, Qing G,, Liu Y, A layerwise/solid-element method for the composite stiffened laminated cylindrical shell structures Composite Structures, 2013, vol, 98, pp, 215-227, DOI: 10,1016/j,compstruct,2012,11,013,

23, Sobhaniaragh B,, Nejati M,, Mansur W,J, Buckling modelling of ring and stringer stiffened cylindrical shells aggregated by graded CNTs Composites Part B: Engineering, 2017, vol, 124, pp, 120-133, DOI: 10,1016/j,compositesb,2017,05,045,

24, Ng T,Y,, Lam K,Y, Dynamic stability analysis of cross-ply laminated cylindrical shells using different thin shell theories Acta Mechanica, 1999, vol, 134, no, 3-4, pp, 147-167, DOI: 10,1007/BF01312653,

25, Martins J.P., Simoes da Silva L., Silvestre N. Energy-based analytical model to predict the elastic critical behaviour of curved panels Journal of Constructional Steel Research, 2016, vol, 127, pp, 165-175, DOI: 10,1016/j,jcsr,2016,07,029,

26, Sengupta J,, Ghosh A,, Chakravorty D, Progressive Failure Analysis of Laminated Composite Cylindrical Shell Roofs

Journal of Failure Analysis and Prevention, 2015, Vol. 15, No. 3. pp. 390-400. doi: 10.1007/s11668-015-9951-6.

27. Ahmed M.K. Simplified equations and solutions for the free vibration of an orthotopic oval cylindrical shell with variable thickness Mathematical Methods in the Applied Sciences, 2011, vol. 34. pp. 1789-1800. DOI: 10.1002/mma.1493.

28. Sofiyev A.H., Birinci F., Sofiyev Al.H., Yusufoglu E. The Dynamic Stability of Cylindrical Shells with Variable Geometric and Mechanical Parameters Multidiscipline Modeling in Materials and Structures, 2006, vol. 2, no. 3. pp. 345-354. DOI: 10.1163/157361106777641369.

29. Qu Y., Wu S., Chen Y., Hua H. Vibration analysis of ring-stiffened conical-cylindrical-spherical shells based on a modified variational approach International Journal of Mechanical Sciences, 2013, vol. 69, pp. 72-84. DOI: 10.1016/j.ijmecsci.2013.01.026.

30. Pogorelov A.V. Differencial'naja geometrija [Differential geometry].Moscow, Nauka, 1969. 176 p.

31. Karpov V.V., Semenov A.A. Refined model of stiffened shells International Journal of Solids and Structures, 2020, vol. 199. pp. 43-56. D0I:10.1016/j.ijsolstr.2020.03.019.

32. Shalashilin V. I., Kuznetsov E. B. Parameter continuation method and best parameterization (in applied mathematics and mechanics)[ Method of continuation of the solution in a parameter and the best parametrization]. Moscow, Editorial URSS, 1999. 222 p.

33. Semenov A.A., Leonov S.S. The Continuous Method of Solution Continuation with Respect to the Best Parameter in the Calculation of Shell Structures Uch. Zap. Kazan. Univ. Ser. Fiz.-Mat. Nauki, 2019, vol. 161, no. 2. pp. 230-249. DOI: 10.26907/2541-7746.2019.2.230-249.

34. Karpov V., Semenov A. Dimensionless parameters in the theory of reinforced shells PNRPU Mechanics Bulletin, 2015, no. 3. pp. 74-94. DOI: 10.15593/perm.mech/2015.3.07

i Надоели баннеры? Вы всегда можете отключить рекламу.