Решение задач механики деформирования оболочек методом функционального нормирования тема автореферата и диссертации по механике, 01.02.04 ВАК РФ

Институт прикладной механики Российской академии наук

На правах рукописи

РЕШЕНИЕ ЗАДАЧ МЕХАНИКИ ДЕФОРМИРОВАНИЯ ОБОЛОЧЕК МЕТОДОМ ФУНКЦИОНАЛЬНОГО НОРМИРОВАНИЯ

01.02.04 - Механика деформируемого твердого тела

Диссертация на соискание ученой степени кандидата физико-математических наук

Научный руководитель:

доктор физико-математических наук

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

КАЗАНЬ - 1999

СОДЕРЖАНИЕ

ВВЕДЕНИЕ 6

ГЛАВ А 1. МЕТОД ФУНКЦИОНАЛЬНОГО НОРМИРОВАНИЯ РЕШЕНИЙ ОБЫКНОВЕННЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ КРАЕВЫХ ЗАДАЧ 22

ГЛАВ А 2. ЦИЛИНДРИЧЕСКИЕ СЛОИСТЫЕ ОРТОТРОПНЫЕ ОБОЛОЧКИ

2.1. Постановка краевой задачи 29

2.2. Разделение переменных. Аналитическое решение обыкновенных дифференциальных уравнений 34

2.3. Исследование влияния обобщенных жесткостей слоистого материала на изменяемость решений обыкновенных дифференциальных уравнений 40

2.4. Функциональное нормирование фундаментальной системы решений обыкновенных дифференциальных уравнений краевой задачи 52

2.5. Алгоритм решения краевых задач 55

ГЛАВ А 3. КОНИЧЕСКИЕ СЛОИСТЫЕ ОРТОТРОПНЫЕ ОБОЛОЧКИ

3.1. Постановка краевой задачи 57

3.2. Разделение переменных. Комплексное преобразование системы обыкновенных дифференциальных уравнений ^

3.3. Фундаментальная система решений дифференциального уравнения в окрестности начала координат 66

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

3.4.1. Степенные ряды 67

3.4.2. Асимптотические разложения 69 3.4.3 Численное интегрирование 74

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

3.5.1. Степенные ряды 77

3.5.2. Асимптотические разложения 79

3.6. Фундаментальная система решений дифференциального уравнения в окрестности бесконечно удаленной точки комплексной плоскости 81

3.7. Универсальное асимптотическое решение дифференциального уравнения при больших значениях независимого переменного и

параметров 84

3.8. Выражения для физических величин задачи 90

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

3.10. Выражения для производных и интегралов от решений дифференциального уравнения 101

3.11. Осесимметричное нагружение конической оболочки

3.11.1 Разрешающее уравнение 105

3.11.2. Решение дифференциального уравнения 107

3.11.3. Осесимметричное кручение оболочки 111

3.11.4. Функциональное нормирование фундаментальной системы решений дифференциального уравнения краевой

задачи \ \ 3

3.12. Антисимметричное нагружение конической оболочки

3.12.1. Разрешающее уравнение U4

3.12.2. Построение фундаментальной системы решений однородного дифференциального уравнения

3.12.3. Частные решения неоднородного дифференциального уравнения и смещения оболочки как жесткого целого Ц8

3.12.4. Функциональное нормирование фундаментальной системы решений дифференциального уравнения краевой

задачи 119

3.13. Алгоритм решения краевых задач 120

ГЛАВ А 4. СФЕРИЧЕСКИЕ СЛОИСТЫЕ ОРТОТРОПНЫЕ ОБОЛОЧКИ

4.1. Постановка краевой задачи 121

4.2. Разделение переменных. Комплексное преобразование обыкновенных дифференциальных уравнений 123

4.3. Решение дифференциальных уравнений 130

4.4. Функциональное нормирование фундаментальной системы решений обыкновенных дифференциальных уравнений краевой задачи 138

4.5. Определение перемещений 141

4.6. Особенности осесимметричного нагружения 141

4.7. Особенности антисимметричного нагружения 146

4.8. Алгоритм решения краевых задач 150

ГЛАВ А 5. ВЫЧИСЛИТЕЛЬНЫЕ

ЭКСПЕРИМЕНТЫ 151

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

ЛИТЕРАТУРА

183 186

ВВЕДЕНИЕ

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

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

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

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

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

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

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

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

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

Научная новизна диссертации состоит в следующем:

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

2. Предложен метод функционального нормирования решений обыкновенных дифференциальных уравнений.

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

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

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

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

Практическая ценность работы состоит в следующем: - обеспечена возможность применения на ЭВМ многочисленных аналитических решений дифференциальных уравнений при решении краевых задач;

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

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

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

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

и

СС «-"51 СС 55

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

Для сферических оболочек уравнения содержат переменные коэффициенты и сводятся к гипергеометрическому уравнению, которое имеет второй порядок. Интересный анализ свойств различных решений этого уравнения применительно ко сферической оболочке приведен в работе В.И.Авдеева [2], где подробно рассмотрен вопрос о точности результатов, получаемых при разных формах решения, и даны рекомендации по выбору ФСР разрешающего уравнения. Однако даже после выбора удачной ФСР необходимо нормирование полученных решений, без которого для длинных оболочек невозможно получение хорошо обусловленной матрицы системы алгебраических уравнений для определения постоянных из граничных условий. Поэтому естественным продолжением этих исследований, так же как и для цилиндрических оболочек, должно стать обоснование возможности расчета по единой методике оболочек произвольной длины с помощью не только правильного выбора ФСР, но и нормирования входящих в нее функций.

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

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

Из работ, посвященных этой проблеме, выделим книгу А.Д.Коваленко, Я.М.Григоренко и Л.А.Ильина [60]. В ней получены результаты, послужившие основой для дальнейшего развития аналитического метода расчета конических оболочек.

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

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

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

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