Расчет тонкостенных и трехмерных конструкций сложной формы на основе аппроксимирующих функций с конечными носителями

Министерство науки и высшего образования Российской Федерации

Федеральное государственное бюджетное

образовательное учреждение высшего образования

«Казанский национальный исследовательский

технологический университет»

Ф. С. Хайруллин, О. М. Сахбиев

РАСЧЕТ ТОНКОСТЕННЫХ 

И ТРЕХМЕРНЫХ 

КОНСТРУКЦИЙ СЛОЖНОЙ 

ФОРМЫ НА ОСНОВЕ 

АППРОКСИМИРУЮЩИХ 

ФУНКЦИЙ С КОНЕЧНЫМИ 

НОСИТЕЛЯМИ

Монография

Казань

Издательство КНИТУ

2020

УДК 624.07–415 
ББК 38.112

Х12

Печатается по решению редакционно-издательского совета 

Казанского национального исследовательского технологического университета

Рецензенты:

д-р физ.-мат. наук, проф. Р. А. Каюмов

д-р техн. наук, проф. Ю. В. Клочков

Х12

Хайруллин Ф. С.
Расчет тонкостенных и трехмерных конструкций сложной формы на 
основе аппроксимирующих функций с конечными носителями : монография / Ф. С. Хайруллин, О. М. Сахбиев; Минобрнауки России, 
Казан. нац. исслед. технол. ун-т. – Казань : Изд-во КНИТУ, 2020. –
196 с.

ISBN 978-5-7882-2959-1

Представлены результаты исследований в области статического расчета 

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

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

Подготовлена на кафедре теоретической механики и сопротивления ма
териалов.

ISBN 978-5-7882-2959-1
© Хайруллин Ф. С., Сахбиев О. М., 2020
© Казанский национальный исследовательский 

технологический университет, 2020

УДК 624.07–415 
ББК 38.112

О г л а в л е н и е

Предисловие ...........................................................................................................5

1. ОПРЕДЕЛЕНИЕ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО 
СОСТОЯНИЯ ТОНКИХ ОБОЛОЧЕК И СТЕРЖНЕЙ .....................................11

1.1. Основные соотношения теории тонких оболочек..................................11

1.2. Аппроксимирующие функции с конечными носителями для 
четырехугольных подобластей .......................................................................14

1.3. Аппроксимирующие функции с конечными носителями для 
треугольных подобластей................................................................................19

1.4. Вариационный метод расчета тонких оболочек сложной формы 
в плане...............................................................................................................22

1.5. Определяющие уравнения для стержней. ...............................................27

1.6. Вариационный метод расчета стержневых систем ................................32

1.7. Построение матрицы жесткости конструкции .......................................35

1.8. Об особенностях численной реализации задачи ....................................38

2. ПАРАМЕТРИЗАЦИЯ СРЕДИННОЙ ПОВЕРХНОСТИ ОБОЛОЧКИ........44

2.1. Исходные соотношения............................................................................44

2.2. Параметризация граничных линий оболочки.........................................48

2.3. Построение сглаживающей функции двух переменных .......................57

2.4. Численные результаты..............................................................................64

3. МОДЕЛИРОВАНИЕ ДЕФОРМИРОВАННОГО СОСТОЯНИЯ 
СОСТАВНЫХ ТОНКОСТЕННЫХ КОНСТРУКЦИЙ .....................................74

3.1. Определение напряженно-деформированного состояния составных 
оболочек............................................................................................................75

3.2. Основные соотношения для ребер жесткости ........................................79

3.3. Расчет тонких оболочек с ребрами жесткости........................................84

3.4. Определение  напряженно-деформированного состояния оболочечностержневых конструкций ................................................................................89

3.5. Расчет стержневых систем, несущих тонкостенные перекрытия .........93

3.6. Расчет рамной конструкции, имеющей двухстороннюю обшивку ......97

4. РЕЗУЛЬТАТЫ РАСЧЕТОВ ..........................................................................102

4.1. Пластины и оболочки канонической и сложной формы в плане........102

4.2. Составные оболочки ...............................................................................114

4.3. Расчет тонкостенных  конструкций с вмятинами ................................118

4.4. Результаты расчетов стержневых систем..............................................124

4.5. Оболочки с вырождающейся областью.................................................128

4.6. Численный метод определения обобщенных жесткостных 
характеристик сотового поликарбоната.......................................................133

4.7. Оболочечно-стержневые конструкции..................................................139

5. ОПРЕДЕЛЕНИЕ  НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО 
СОСТОЯНИЯ ТРЕХМЕРНЫХ КОНСТРУКЦИЙ ..........................................147

5.1. Основные соотношения теории упругости в криволинейной 
системе координат..........................................................................................147

5.2. Построение аппроксимирующих функций 
для трехмерных подобластей........................................................................149

5.3. Вариационный метод расчета трехмерных конструкций сложной 
формы..............................................................................................................159

5.4. Особенности численной реализации .....................................................162

5.5. Результаты расчетов................................................................................167

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

П р е д и с л о в и е

Тонкостенные конструкции в настоящее время являются одними 

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

Однако реальные объекты, кроме тонкостенных конструкций,

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

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

В последние десятилетия при решении задач механики деформи
руемого твердого тела наибольшее развитие и распространение получил метод конечных элементов (МКЭ), который сочетает универсальность и эффективность с простотой и удобством при численной реализации задачи. Фундаментальным исследованиям по МКЭ и вопросам 
численной реализации метода посвящено большое количество работ, в 
частности [8, 26, 29, 30, 38, 39, 67, 80, 88, 96, 98, 106]. 

Одним из универсальных численных методов расчета некоторых 

видов конструкций является метод конечных разностей (МКР). 

При использовании этого метода исследуемая область разбивается на прямоугольные подобласти, в пределах которых производные от искомых 
функций заменяются разностными отношениями. По данной теме опубликовано довольно много работ. Некоторые проблемы построения и решения разностных схем рассмотрены в работах [6, 9, 22, 34, 47, 55, 63, 105].

В достаточно общей постановке вопросы расчета оболочек слож
ной геометрии исследовались в работах [50, 70–72]. В этих работах для 
оболочек неканонической формы и неканонических очертаний параметризация срединных поверхностей производилась на основе теории 
конечных деформаций поверхностей.

К одним из первых численных методов расчета тонкостенных 

конструкций относятся методы коллокации, в которых неизвестные параметры, определяющие искомые функции, находятся из условия удовлетворения исходных уравнений в заданной системе точек. Начиная с 
первых публикаций [47, 48] методы коллокации успешно использовались при решении задач расчета пластин и оболочек сложной формы, в 
том числе в работах [16, 33, 89, 90, 105].

Другим эффективным методом расчета пластин и оболочек слож
ной формы является метод граничных элементов (МГЭ), в основе которого лежит известный в задачах математической физики метод потенциалов. В отличие от метода конечных элементов, в МГЭ дискретизации подлежат лишь границы рассматриваемых объектов и задача сводится к решению граничных интегральных уравнений. По теоретическим основам метода и вопросам его практического применения имеются многочисленные публикации, в том числе монографии [4, 10, 14, 
15, 20, 45, 54, 56, 104, 111].

Также можно отметить следующие методы расчета оболочек 

сложной формы. В работах [61, 93–94] предлагается постановка задачи 
и метод расчета оболочек с резными срединными поверхностями. Интегрально-проекционный метод для решения задач расчета оболочек 
используется в работах [74, 83, 95]. Возможность использования для 
определения деформаций оболочек соотношений из теории пластин показана в работах [100–101].

Расчет конструкций, составленных из нескольких видов оболо
чек или пластин, т. е. составных конструкций, производится в основном 
вариационными или численными методами, в частности методом конечных элементов или методом суперэлементов. Вопросы постановки 
и численной реализации данных задач рассмотрены, например, в монографиях [17, 60, 62, 81]. 

Методы решения задач механики деформируемого твердого тела, 

основанные на использовании вариационных принципов [1, 19], называются вариационными. К ним относятся метод Ритца, метод Бубнова–
Галеркина, вариационно-разностные методы, метод конечных элементов в вариационной постановке и др. При их использовании возникает 
вопрос выбора аппроксимирующих функций, которые должны обладать определенными свойствами и удовлетворять определенным условиям. Например, при использовании вариационного уравнения Лагранжа, построенного на основе уравнений теории оболочек типа Тимошенко, аппроксимирующие решение функции должны обладать 
гладкостью класса С(0), составлять полную систему функций и удовлетворять геометрическим граничным условиям. Если оболочка имеет 
сложную форму, то выбор таких функций вызывает определенные 
трудности. Одним из методов построения аппроксимирующих функций является метод R-функций В. Л. Рвачева [84–86]. 

Опубликовано довольно много работ, посвященных построению 

аппроксимирующих функций и использованию этих функций для решения задач расчета оболочек сложной формы. В дополнение к тем методам, которые касались этой темы и изложены выше, можно отметить 
работы [58, 59, 87, 143, 144].

В настоящее время разработано достаточно много конечных 

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

В работах [133, 134] представлены трехмерные конечные эле
менты в виде параллелепипеда, в которых в роли аппроксимирующих 
функций выступает Эрмитовый трехмерный кубический сплайн.

Некоторые из объемных конечных элементов, используемых для 

расчета оболочек, как, например, в работах [11, 28, 36, 42–44], представлены в формулировке МКЭ без использования упрощающих гипотез о недеформировании нормали.

При решении задач расчета оболочек сложной геометрии могут 

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

вектора r срединной поверхности оболочки должна производиться с достаточно большой точностью. Как показывают численные эксперименты, возможная осцилляция даже во вторых производных от r может 
привести к большим погрешностям в решении задачи, так как эти производные определяют радиусы кривизны оболочки. Среди основных 
аналитических и численных методов параметризации срединной поверхности оболочек сложной формы можно выделить следующие: метод деформации поверхности отсчета; использование кубических и 
других сплайн-аппроксимаций; использование метода конечных разностей и метода конечных элементов; использование сглаживающей аппроксимации и т. д. Эти и другие методы параметризации поверхностей 
и кривых рассмотрены в монографиях [25, 32, 37, 49, 58, 91, 114, 135].

Анализ приведенных методов расчета показывает, что, хотя и су
ществуют различные методы расчета тонкостенных конструкций сложной формы, однако универсального метода, применимого для любого 
случая, нет. Каждый из этих методов имеет свои положительные и отрицательные стороны, применим для определенных задач. Даже такой 
универсальный метод, как метод конечных элементов, имеет свои недостатки. Как отмечено в монографии А. И. Голованова и соавторов 
[31], несмотря на большое количество работ по методу конечных элементов и большое количество предложенных в этих работах конечных 
элементов, «лишь ограниченное количество их действительно эффективно в расчетах тонких непологих оболочек».  

В работе Эдельмана, Казаринеса, Уолтона [131] исследуется вли
яние порядка аппроксимирующей функции на точность решения. 
На конкретных примерах показывается, что использование высокоточных конечных элементов, использующих полиномы высокого порядка, 
позволяет на малом количестве элементов получить более точные результаты при меньших размерах матрицы жесткости, по сравнению с 
применением большого числа более простых конечных элементов. Однако при использовании функций высокой степени аппроксимации в 
узловых точках требуется задавать производные высоких порядков, 
например, производные второго порядка. Это приводит к усложнению 
формулировки и выполнения граничных условий, а при расчете составных оболочек создает проблемы с выполнением условий сопряжения 
на изломе срединной поверхности оболочки. 

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

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

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

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

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

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

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

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

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

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

В пятой главе дается алгоритм построения аппроксимирующих 

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