Приведение краевых задач к начальным и исследование концентрации напряжений в тонкостенных конструкциях мультипликативным методом тема диссертации и автореферата по ВАК РФ 05.13.18, кандидат физико-математических наук Петров, Виталий Игоревич

• Петров, Виталий Игоревич
• кандидат физико-математических наук
• 2008

Оглавление диссертации кандидат физико-математических наук Петров, Виталий Игоревич

Введение и краткий обзор литературы.

01. Физическая постановка задачи.

02. Математические модели.

03. Методы анализа.

Глава 1. Математические модели механики деформирования оболочек.

1.1 Математическая модель механики деформирования слоистых оболочек Я.М.Григоренко и А.Т.Василенко.

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

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

1.2 Математическая модель механики деформирования изотропных оболочек В.З.Власова.

1.2.1 Разрешающая система уравнений общей моментной технической теории В.З.Власова.

1.2.2 Матричная форма уравнений механики деформирования цилиндрической, канонических и сферических оболочек.

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

Глава 2. Решение линейных обыкновенных дифференциальных уравнений в матричной форме.

2.1 Матричный ряд Тейлора.

2.2 Мультипликативный интеграл Вольтерра и матричный бином Ньютона.

2.3 Частное решение.

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

3.1 Приведение краевых задач к начальным до счета.

3.2 Приведение краевых задач к начальным с помощью ЭВМ.

3.3 Математическое моделирование локальной нагрузки.

3.4 Идеализация локальной нагрузки сосредоточенной.

3.5 Вычислительные эксперименты.

3.5.1 Определение относительной критической длины оболочки.

3.5.2 Удвоение относительной критической длины оболочки.

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

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

4.1 Теорема.

4.2 Алгоритм формирования начальных условий.

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

5.1 Теоретические основы алгоритма.Т.

5.2 Матричная форма уравнений механики деформирования шпангоута.

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

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

5.5 Вычислительный эксперимент.

Глава 6. Исследование концентрации напряжений в транспортно-пусковых стаканах (ТПС) летательных аппаратов.

6.1 Локальное воздействие на днище ТПС по круглым площадкам.

-46.2 Локальное воздействие на днище ТПС по площадкам, очерченным линиями главных кривизн.

6.3 ТПС под внутренним давлением.

6.4 ТПС при запуске летательных аппаратов.

О погрешностях.

Введение диссертации (часть автореферата) на тему «Приведение краевых задач к начальным и исследование концентрации напряжений в тонкостенных конструкциях мультипликативным методом»

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

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

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

0.1. Физическая постановка задачи.

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

A.M. Доценко, В.А. Коргопольцев, М.Н. Коган, А.Б. Корнилов, A.A. Орлов, В.Н. Семенов, JI.JI. Теперин, М.В. Устинов экспериментально в ЦАГИ им. проф. Н.Е. Жуковского показали, что использование нанотехноло-гий в местах концентрации напряжений для предотвращения появления и развития трещин увеличивает прочностной ресурс тонкостенных конструкций до 3-х раз.

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

0.2. Математические модели.

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

Фундаментальные основы теории оболочек приводятся в работах В.З. Власова [84-87], С.П. Тимошенко [239], А.Л. Гольденвейзера [96, 97], В.В. Новожилова [180-182], А.Н.Лурье, И. Геккелера и получили развитие в работах многих других авторов, приведенных в списке литературы. Изучение напряженно-деформированного состояния оболочек при силовых и^ температурных воздействий получило развитие в работах И.А. Биргера [23], В.Л. Бидермана [19-22], В.В. Болотина [24], И.Н. Векуа [41], Я.М. Григорен-ко и А.Т. Василенко [107-112], А.Д. Коваленко, В.В. Мяченкова, [175, 176], К.Ф. Черных и многих других авторов, также приведенных в списке литературы.

В качестве математической1 модели механики деформирования оболочек нами принимаются дифференциальные уравнения, полученные-для анизотропных слоистых оболочек Я.М. Григоренко и А.Т. Василенко [108]. Они обладают в сравнении с многими другими известными тем свойством, что позволяют исследовать концентрацию напряжений в анизотропных, орто-тропных и изотропных оболочках классических форм. Важным является и то, что эти уравнения уточнены повторным выводом Ю.А. Гусевым [114]. Им обнаружено и исправлено более 10 существенных опечаток.

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

Я8 + а316 + а2Я4 + а{Х2 +а0= О разрешающих обыкновенных дифференциальных уравнений, полученных методом Фурье разделения переменных для цилиндрических оболочек

- В.З. Власовым а0 = п4{п2 -1)2, а{ = -4п6 + 8и4 - 2цп4 - 2п2, а2 = 6п2(п2 -!) + (! +4к4), а3=-4п2 + /л\

- Б.Г. Галеркиным а0 = п6(п2 -1), ах = -4 л6 + (5 - 3/ф4 - (1 - 2/л)п2, а2 = 6п4-(4-/и)п2+ 4к4, а3 = -4п2 +2//;

- А.Л. Гольденвейзером а0=п4(п2-1)2, а1=-4п2(п2 -I)2, а2=6п4-Ъп2+2п2/л2+4к4, а3=-4п2;

- В.М. Даревским а0 = п4(п2-1)2, а^-Лп2^2-!)2, а2=6п4 -1п2 +п21л2 +4к4, а3=-4п2;

- В. Флюгге

2п4 - (4 - /и)п2 + (2 - //)!, а2 = 6и2(и2-1) + 4А:4, аъ=-2{2п2-ц), где за )л —, // - коэффициент Пуассона, п - номер гармоники в методе Фурье разделения переменных. я0 = гс4(и2-1)2, -2п2

На различие математических моделей и, следовательно, коэффициентов в соответствующих характеристических уравнениях обратил внимание, вероятно впервые, в 1964 году Д.А. Рояк [220]. Однако анализ различных математических моделей с целью установления наиболее совершенной, адекватной реальному объекту, нам представляется непродуктивным. Более результативным и менее трудоемким представляется сравнительный количественный вычислительный эксперимент. Ожидается, что количественные параметрические исследования, обязательно с контролируемой погрешностью, дадут количественно мало отличающие результаты. Отличия, как показали наши некоторые решения для цилиндрической оболочки, не превышали 5%, что и предсказывал В.З. Власов [84]. Очевидно, что вычислительный эксперимент сталкивается с проблемой контроля погрешностей в результатах.

Требуется ответ и на другие вопросы, связанные с математическими моделями. Один из вопросов связан с правомерностью использования теории оболочек для исследования в них концентрации напряжений. Такой вопрос возникает при локальном воздействии на оболочки. Ответ дали В.И. Леонов и Х.С. Хазанов [153]. Они показали, что математические модели теории оболочек могут использоваться при исследовании в них напряжений, если радиус круглой площадки внешнего воздействия более половины толщины оболочки. Если площадка воздействия меньше, то концентрацию напряжений под ней и в малой её окрестности, не превышающей толщины оболочки, необходимо исследовать, решая трехмерную задачу.

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

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

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

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

- концентрироваться - собираться, скапливаться в каком-нибудь месте;

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

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

0.3. Методы анализа.

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

Метод Ритца-Тшюшенко. Предложен В. Ритцем и распространён С.П. Тимошенко на задачи строительной механики. Он позволяет получить приближенное решение в перемещениях на основе вариационного принципа Ла-гранжа. Идея метода, определяющая его приближенный характер, связана с тем, что искомые перемещения отыскиваются в классе заранее заданных функций, которые выбираются на основе опыта, интуитивно или на основе решений более простых задач. Для построения приближенного решения перемещения" представляются в виде рядов заданных функций, коэффициенты которых определяются при решении системы алгебраических уравнений, полученной по принципу Лагранжа. При этом необходимо помнить требования к задаваемым функциям. Они должны быть дифференцируемыми и удовлетворять геометрическим граничным условиям задач. Возникает и чисто математический вопрос, связанный с полнотой заданных функций и сходимостью рядов к точному решению. Метод Ритца излагается в работе [215]. Свое развитие он получил в работе [239], посвященной изучению устойчивости упругих систем. Прочностной расчет многослойных пластин с помощью метода Ритца выполнялся в работе В.В. Болотина, Ю.Н. Новичкова [24]. Здесь же определяются частоты собственных колебаний многослойных плит.

Метод Бубнова-Галеркина. Метод был построен И.Г. Бубновым и широко использовался при решении различных задач строительной механики Б. Г. Галеркиным. Метод основан на принципе возможных перемещений и сведении задач к решению систем алгебраических уравнений. Искомые перемещения, как и в методе Ритца-Тимошенко, ищутся в форме рядов в классе заданных функции. Заданные функции должны удовлетворять геометрическим граничным условиям и быть непрерывными. Если приравнять нулю возможную работу на каждом из заданных перемещений, то получаем метод, принятый называть обобщенным методом Бубнова-Галеркина. Если же при выборе заданных функций потребовать, чтобы они удовлетворяли и статическим граничным условиям, то полученный таким образом метод называется методом Бубнова-Галеркина. При этом коэффициенты рядов решений могут определяться в соответствии с методом Ритца-Тимошенко или в соответствии с методом Бубнова-Галеркина. Теоретические основы и определение области применения метода Бубнова-Галеркина изложены в монографиях [143, 211, 232, 248]. Метод получил широкое применение для расчета оболочек и тонкостенных конструкций на прочность, устойчивость и колебания. Этим методом решались краевые задачи в работах [10, 143]. Необходимо отметить, что в задачах определения напряженно-деформированного состояния при нестационарном нагружении методы Ритца-Тимошенко и Бубнова-Галеркина, базирующиеся на глобальных базисных функцях, обладают существенными' недостатками. Решения, полученные с помощью этих методов, сильно зависят от конкретного вида базисных функций, аппроксимирующих перемещения.

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

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

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

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

Дифференциально-разностный метод (метод прямых или полос). Вероятно впервые метод сформулирован М.Г. Слободянским [234] и развивался С.Г. Михлиным [172, 173], И.С. Березиным, Н.П. Жидковым [17, 18], А.П. Филиным [245, 246] и др [110]. По методу по одной из независимых переменных в двумерных задачах и по двум независимым переменным в трехмерных задачах производные заменяются разностными выражениями. Такая процедура обеспечивает замену краевой задачи для дифференциальных уравнений в частных производных краевой задачей для обыкновенных дифференциальных уравнений. Обыкновенные дифференциальные уравнения интегрируются пошаговыми методами [6], чаще всего методом- Рунге-Кутта четвертого порядка. Краевые условия по направлениям аппроксимации производных конечными разностями реализуются с использованием разностных аппроксимаций производных в краевых условиях. По направлению .численного интегрирования используются алгоритмы переноса краевых условий методом Коши, например, в произвольную точку краевого интервала встречной прогонкой. При этом преодолевается проблема устойчивости счета, например ортонормированием [94].

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

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

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

Метод "прогонки " Гелъфанда - Локуциевского. Впервые изложен в виде дополнения к книге [95]. Метод относится к методам переноса краевых условий, в основе которых лежат численные пошаговые методы интегрирования линейных обыкновенных дифференциальных уравнений. И. Бабушка, Э. Витасек, М. Прагер [6] отмечают, что при интегрировании жестких дифференциальных уравнений, к которым относятся и уравнения механики деформирования тонкостенных конструкций, возникают вычислительные трудности, связанные с неустойчивостью счета. Метод не всегда обеспечивает устойчивость счета. Особенно это проявляется при решении задач устойчивости и колебаний оболочек [21]. Метод Гельфанда-Локуциевского развивается в работах [47]: строятся его модификации, а также дискретный аналог, по существу алгоритма - новый метод.

Метод Абрамова переноса граничных условий. Метод предложен A.A. Абрамовым [2]. Всегда обеспечивает устойчивость счета при решении краевых задач механики деформирования.оболочек. Недостатком метода является увеличение в два раза порядка разрешающей системы обыкновенных дифференциальных уравнений с появлением сложных математических функций в преобразованных правых частях уравнений и, следовательно, неоправданный рост времени счета и необходимой оперативной памяти ЭВМ. При этом также теряется простой и понятный физический, смысл входящих в разрешаюч i » щую систему уравнений неизвестных. В работе [49] показано, что> устойчивость счета обеспечивается пошаговым ортонормированием в процессе интегрирования дифференциальных уравнений. -Метод Годунова. Метод численного решения краевых задач для систем линейных обыкновенных дифференциальных уравнений. Впервые предложен С.К. Годуновым [94]. Он получил развитие применительно к решению задач прочности, устойчивости и колебаний пластин, оболочек и тонкостенных конструкций в работах В.Л. Бидермана [19-22], В.В Болотина, Ю.Н. Но-вичкова [24], Я.М. Григоренко, А.Т. Василенко, Е.И. Беспаловой [111], В.И. Мяченкова, И.В. Григорьева, В.П. Мальцева [175, 176] и др. На основании метода созданы и опубликованы пакеты прикладных программ для решения краевых задач теории оболочек и механики деформирования тонкостенных конструкций. Устойчивость счета при решении краевых задач обеспечивается предложенной С.К. Годуновым процедурой ортогонализации и нормиро ; • . .■ -17- ^ . : вания решений в. процессе интегрирования систем линейных обыкновенных дифференциальных уравнений методом Рунге-Кутта.

Метод Виноградова-Клюева. Метод численного и аналитического решения краевых задач механики деформирования тонкостенных конструкций сопряжением участков устойчивого счета. Метод предложен Ю.И. Виноградовым и развит Ю.И. Клюевым для решения задач устойчивости и колебаний оболочек и тонкостенных конструкций [57-65, 67, 68]. Метод можно считать численным, если решения» линейных обыкновенных дифференциальных уравнений с постоянными коэффициентами определяются» с помощью матричного бинома Ньютона, а с переменными коэффициентами - с помощью^ мультипликативного интеграла Вольтерра, аналога- интеграла как суммы* бесконечно малых [92]. Метод можно считать аналитическим., если решение:линейных: обыкновенных дифференциальных уравнению определяется? с помощью матричного рядя: Тейлора. Сопряжение участков устойчивого счета' приводит к системе алгебраических уравнений. При решении'системы алгебраических уравнений4учитываютсязаданные краевые условия и нагрузка на оболочку или; тонкостенную конструкцию. Таким образом; вариантные крае. ? Л.вые задачи при различных внешних воздействиях решаются при повторном решении алгебраических уравнений; в результате которых определяются искомые величины задач с априори назначенной погрешностью.

Мультипликативный метод Виноградова. Метод предложен Ю.И. Виноградовым [51]. Эффективные алгоритмы строились в работах А.Ю. Виноградова [45, 54, 55] и работах Ю.А. Гусева [114]. Идея метода состоит в переносе краевых условий; в произвольную точку краевого интервала с помощью решений линейных обыкновенных дифференциальных уравнений механики деформирования оболочек, записанных в канонической форме. Решения дифференциальных уравнений определяются, с помощью матричного бинома Ньютона, определения мультипликативного интеграла Вольтерра или матричного ряда Тейлора. При этом определяются значения функций Коши

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

Метод Виноградова А.Ю. и Виноградова Ю.И. Метод разработан Ю.И. Виноградовым и А.Ю. Виноградовым для решения краевых задач на основе линейных обыкновенных дифференциальных уравнений, содержащих только чётные производные [42]. Многие прикладные задачи изгиба балок,"изотропных, ортотропных, трехслойных пластин, круговых цилиндрических оболочек с уравнениями механики их деформирования Власова, Гольденвейзера, Новожилова, Флюгге, Бейларда, Даревского, Тимошенко-Лява, Доннелла-Власова-Лурье и также варианты полубезмоментных теорий, задачи механики деформирования оболочек на основе моментной технической теории; общей моментной теории, теории пологих оболочек математически моделируются после разделения переменных линейными обыкновенными дифференциальными уравнениями, содержащими только четные производные [21, 177, 184]. Метод решения таких задач построен на идее сведения разрешающего дифференциального уравнения, содержащего только четные производные, к системе дифференциальных уравнений, каждое из которых второго порядка, и записи её в матричной форме. Решение такого матричного однородного уравнения определяется по известным [92] формулам через матричные синусы и косинусы, если дифференциальное уравнение имеет постоянные коэффициенты. Если коэффициенты переменные, то они осредняются на участках и строятся соответствующие алгоритмы решения краевых задач. Построено несколько модификаций метода, предлагаются рекуррентные формулы для сокращения затрат машинного времени на построение общего решения системы обыкновенных дифференциальных уравнений.

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

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

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

Однако метод Гельфанда-Локуциевского и Годунова объединяет тот факт, что для начала вычислений они требуют начальные значения прого-ночной матрицы и вектора. Для метода Абрамова этого не требуется.

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

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

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

Методы Виноградова-Клюева, Виноградова и Виноградовых принципиально отличаются от других тем, что для дифференциальных уравнений решения определяются по формулам без учета краевых условий. Общим существенным признаком этих формул является то, что с их помощью определяются значения функций Коши-Крылова, удовлетворяющих произвольным начальным условиям задач. Если в основе алгоритмов используется матричный бином Ньютона и интеграл Вольтерра для определения решений обыкновенных дифференциальных уравнений, то методы будут относиться к численным, так как при этом нет возможностей определения погрешностей методов. Если в основе тех же алгоритмов используются матричный ряд Тейлора или матричные сходящиеся ряды для определения матричных синусов и косинусов для определения решений обыкновенных дифференциальных уравнений, то методы будут относиться к аналитическим, т.к. при этом решения дифференциальных уравнений определяются с помощью сходящихся рядов с контролируемой погрешностью. Эффективностьf этих методов в сравнении с методами Гельфанда-Локуциевского, Абрамова и Годунова по программной реализации очевидно проще, затраты, машинного времени сокращаются при решении одних и тех же задач на один-два порядка, значительно сокращаются требованиям оперативной памяти ЭВМ и, что принципиально важно, эти методы позволяют решать задачи механики деформирования тонкостенных конструкций с априори заданной погрешностью. Однако их алгоритмы построены для определения^ напряженно-деформированного состояния во всей области краевого интервала.

Метод граничных элементов (МГЭ). Широко используется при решении задач теории упругости. Основой МГЭ является интегральное представление решения системы дифференциальных уравнений, описывающих поведение исследуемого объекта. Впервые такое представление в теории потенциала было получено А.Г.Грином [103]. Применительно к задачам теории упругости В.Д .Купрадзе [150] ввел некоторые интегральные уравнения на основе упругих потенциалов простого и двойного слоя, доказал существование решений этих уравнений и предложил приближенный метод решения статических задач для однородных упругих тел. Развитие метод получил в работах В.З. Партона, П.И. Перлина на основе упругих потенциалов. Прикладное значение решения задач статики теории упругости определено в работах Cruse Т.А., Rizzo F.Y. [147, 148]. В монографях [16, 28, 144, 168] суммируется опыт применения МГЭ в самых различных разделах механики. Достоинством МГЭ являются понижение размерности решаемой задачи за счет поиска неизвестных только на границе исследуемой области и возможность определения параметров напряженно-деформированного состояния во внутренних точках тела, а не в местах их концентрации.

Метод конечных разностей (МКР) и метод конечных элементов (МКЭ). Универсальные численные методы. Используются при создании программных комплексов для решения многомерных задач параметрического анализа. В основе МКР - разностная аппроксимация производных с образованием разностных сеток при решении краевых и начально-краевых заач. В основе МКЭ, широко распространенном при исследовании механики деформирования оболочек и тонкостенных конструкций, - аппроксимация искомых величин с помощью конечного числа выбранных локализованных функций, чаще всего полиномов. Теоретическое обоснование МКР дано в работах С.К. Годунова, B.C. Рябенького [95], JI.B. Канторовича, В.Г. Крылова [135], А.А Самарского [226-231]. Этим методом решались краевые задачи статики, устойчивости и колебания пластин, оболочек и тонкостенных конструкций в работах Э.И. Григолюка, В.М. Толкачева [106], В.А. Иванова, Б.В. Гулина

132], В.П. Баженова и др. [7-9, 32, 36; 216, 217]. Фундаментальные результа ты теории МКЭ получены в работах О.Зенкевича [127-129], С.Г.Михлйна [172, 173], К.Ректориса [213], J.T. Oden [188, 189], Стренг, Фикс Дж. [235]. Этим методом решались задачи механики деформирования оболочек и тонкостенных конструкций в работах А.И. Голованова [98-100], М.С. Корниши-на [141], Бате К., Вилсон Е. [11], Постнов В.А.[202-206], Р.Б. Рикардса [214] и др. [195, 219, 237]. Список таких работ велик и продолжает быстро увеличиваться. Развитие универсальных МКР и МКЭ идет по пути решения многомерных задач в автоматическом режиме путем распараллеливания алгоритма и организации вычислительного процесса в диалоговом режиме с использованием многопроцессорной техники, с обработкой многомерных массивов данных и визуализацией результатов счета. Программные комплексы совершенствуются в направлении расширения их возможностей и упрощения использования. В то же время мало внимания уделяется проблеме контроля за погрешностями счета [25, 26, 34, 102, 118, 121, 130; 131, 139, 159, 190, 221, 222, 253]. Универсальные МКР и МКЭ, как и выше описанные, которые можно отнести к специализированным, реализованные в программных комплексах, аппроксимируют оболочки или тонкостенные конструкции во всех их элементах, что не требуется при исследовании напряжений в местах их концентрации. Они особенно не надежны при повышенных требованиях к точности исследования концентрации напряжений* и определении их максимальных значений.

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

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

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

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

Группа методов Гельфанда-Локуциевского, Абрамова, Годунова, Виноградова и Виноградовых, алгоритмы решения краевых задач которых связаны с переносом краевых условий в произвольную точку краевого интервала, по своей сути основаны на замене краевой задачи для одномерного процесса решением ряда задач с начальными условиями, то есть, решением ряда задач методом Коши. Идея решения, краевых задач методом начальных параметров весьма простая: известные краевые условия'или условия, дополненные в символическом виде неизвестными начальными параметрами, переносятся в произвольную точку краевого интервала. Совокупность начальных условий позволяет определить напряженно-деформированное состояние в произвольной- точке. Метод начальных параметров используется в строительной механике тонкостенных конструкций в виде различных модификаций, предложенных в разное время Коши, Клебшем, П.Ф. Папковичем [192], Ш.Е.Микеладзе [169], В.З. Власовым [84-87] и* другими. Метод начальных параметров при замене краевой задачи совокупностью задач Коши не обеспечивает устойчивости счета для1 жестких дифференциальных уравнений механики деформирования тонкостенных элементов конструкции. При этом добиться удовлетворительной точности численных расчетов становится весьма затруднительным.

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

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

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

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

Научную новизну, работы.составляют:

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

- перенос краевых условий: к местам концентрации напряжений мультипликативным методом, (доказана,теорема); формирование: соответствующих начальных условий и исследование концентрации напряжений в элементах тонкостенных, конструкций- мультипликативным? методом* с: априори? контролируемой погрешностью,

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

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

- исследование концентрации напряжений в тонкостенных элементах транс-портно пусковых стаканов летательных аппаратов:

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

Прикладная ценность работы состоит в:

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

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

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

- возможности априори определять погрешности исследования концентрации напряжений,

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

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

Основные положения, выносимые на защиту.

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

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

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

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

Апробация работы и публикации. Основные результаты диссертации докладывались и обсуждались на 3-й, 4-й, 5-й и 6-й Международных конференциях "Авиация и космонавтика - 2004, 2005, 2006 и 2007" (Москва, 2004,

2005, 2006, 2007 г.г.); Первой международной научно-технической конфе ренции "Аэрокосмические технологии" (Москва - Реутов, 2004 г.); Международной конференции "DYNAMICAL SYSTEM MODELLING AND STABILITY INVESTIGATION" (Киев, 2005 г.); XIV и XV Международных конференциях по вычислительной механике и современным прикладным системам (Алушта, 2005 и 2007 г.г.); Международном симпозиуме "Образование через науку" (Москва 2005 г.); XII и XIII Международных симпозиумах "Динамические и технологические проблемы механики конструкций и сплошных сред" (Москва, 2006 и 2007 г.г.); 3-й Международной конференции "Ракетно-космическая техника: фундаментальные и прикладные проблемы" (Москва, 2007), кафедре "Аэрокосмические системы" МГТУ им. Н.Э. Баумана.

По теме диссертации опубликовано 17 работ.

Объем и структура диссертации. Диссертация состоит из введения и краткого обзора литературы, 6-ти глав, заключения, списка литературы, включающего 263 наименования. Объем работы составляет 167 страниц машинописного текста и 38 страниц рисунков.

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Петров, Виталий Игоревич

ЗАКЛЮЧЕНИЕ

1. В диссертации обсуждаются математические модели механики деформирования оболочек. Дается обоснованный выбор математической модели Я.М. Григоренко и А.Т. Василенко механики деформирования слоистых анизотропных оболочек, записанных в канонической форме, и технической моментной теории В.З. Власова механики деформирования изотропных оболочек как простейшей.

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

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

2. Обсуждаются математические основы универсальных численных методов конечных разностей и конечных элементов, методов Ритца-Тимошенко, Бубнова-Галеркина, Власова-Канторовича, Папковича-Треффца, дифференциально-разностных методов Слободянского и Виноградова, вариационно-разностного метода, метода прогонки Гельфанда-Локуциевского, Абрамова и Годунова, метода Виноградова-Клюева, мультипликативного метода Виноградова, метода Виноградова А.Ю. — Виноградова Ю.И. и метода граничных элементов.

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

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

Уравнения в частных производных механики деформирования оболочек методом Фурье разделения переменных приводятся к обыкновенным, для которых решения определяются по формулам.

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

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

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

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

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

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

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

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

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

Крылова - решения матричного дифференциального уравнения механики деформирования оболочек. Краевые условия переносятся к местам концентрации напряжений мультипликативным методом Виноградова через участх ки, например, для цилиндрической оболочки — «1. При этом устойчивость Л счета обеспечивается построчным ортонормированием матриц без искажения значений искомых величин и, следовательно, без необходимости запоминать ортонормированные матрицы. С помощью перенесенных краевых условий у мест концентрации напряжений формируются соответствующие начальные условия. Напряжения исследуются мультипликативным методом.

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

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

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

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

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

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

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

Материал транспортно пускового стакана:

- титан;

- цилиндрическая оболочка - композиционный материал, сферическая оболочка - титан.

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

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

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

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

Список литературы диссертационного исследования кандидат физико-математических наук Петров, Виталий Игоревич, 2008 год

1. Абовский Н.П., Андреев Н.П., Деруга А.П., Савченков В.И. Численные методы в теории упругости и теории оболочек. -Красноярск: Изд-во Красноярского университета, 1986. - 383 с.

2. Абрамов A.A. О переносе граничных условий для систем линейных дифференциальных уравнений (вариант метода прогонки) // Журнал вычислительной математики и математической физики, 1961. -Т.1. -№3. С. 542-545.

3. Андреев JI.B., Зюзин В.А., Муляр Ю.М. Применение метода начальных параметров к расчету сложных оболочечных конструкций // Гидроаэромеханика и теория упругости. -Днепропетровск: Изд-во Днепропетровского университета, 1971. -С.17-23.

4. Артюхин Ю.П., Грибов А.П. Исследования изгиба пластин сложной формы под действием температурного поля и нормального давления методом граничных элементов // Прикладные задачи напряженного состояния упругих тел, Саратов, 1987. -С.50-54.

5. Babuska I., Aziz A.K. Survey lectures on the mathematical foundation of the finite element method with application to partial differential equations // Edt. A.K.Aziz new York: Academic Press, 1972. -P.359.

6. Бабушка И., Витасек Э., Прагер М. Численные процессы решения дифференциальных уравнений. -М.: Мир, 1969. -368 с.

7. Баженов В.Г., Чекмарев Д.Т. О конечно-разностном решении волновых уравнений теории оболочек типа Тимошенко // Прикладные проблемы прочности и пластичности. Всесоюз. межвуз. сб. / Горьковский университет, Горький, 1981. -С.41-50.

8. Баженов В.Г., Шинкаренко А.П. Вариационно-разностный метод решения двумерных задач динамики упругопластических оболочек // Прикладные проблемы прочности и пластичности, Горький, ГГУ, 1976. -Вып.№. -С.14-21.

9. Баженов В.Г., Журавлев Е.А. Вариационно-разностный метод решения нелинейных осесимметричных задач динамики слоистых оболочек // Прикладные проблемы прочности и пластичности, 1979. -Вып.13. -С.36-45.

10. Бастатский Б.И. Модификация метода Бубнова-Галеркина в задачах теории пологих оболочек//Строительная механика и расчет сооружений. -М., 1985.-№5.-С. 4-8.

11. Бате К., Вилсон Е. Численные методы анализа и метод конечных элементов. -М.: Стройиздат, 1982. -494 с.

12. Бахвалов Н.С. Численные методы. Анализ, алгебра, обыкновенные дифференциальные уравнения. -М.: Мир, 1984. -496 с.

13. Бахвалов Н.С. Осреднение дифференциальных уравнений с частными производными с быстро осциллирующими коэффициентами // ДАН СССР, 1975. -Т.221. -№3. -С.516-519.

14. Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы. -М.: Наука, 1987. -600 с.

15. Беллман Р. Введение в теорию матриц. -М.: Наука, 1976. -352 с.

16. Бенерджи П., Баттерфилд Р. Методы граничных элементов в прикладных науках. -М.: Мир, 1984. -494 с.

17. Березин И.С., Жидков Н.П. Методы вычислений. -М.: Физматгиз, 1962.-Т.1.-464 с.

18. Березин И.С., Жидков Н.П. Методы вычислений. -М.: Физматгиз, 1962. -Т.2. -635 с.

19. Бидерман В Л. Применение метода прогонки для численного решения задач строительной механики // Изв. АН СССР Механика твердого тела. -М., 1967.-№5.-С. 62-66.

20. Бидерман В.Л. Некоторые вычислительные методы решения задач строительной механики, приводимых к обыкновенным дифференциальнымуравнениям. -В кн.: Расчеты на прочность. -М.: Машиностроение, 1976.-С. 8-36.

21. Бидерман В.Л. Механика тонкостенных конструкций. -М.; Машиностроение, 1977. -488 с.

22. Бидерман В.Л. Прикладная теория механических колебаний. -М.: Высшая школа, 1972. -416 с.

23. Биргер И.А. Стержни, пластинки и оболочки. -М.: Наука, 1992. -392 с.

24. Болотин В.В., Новичков Ю.Н. Механика многослойных конструкций. -М.: Машиностроение, 1980. -375 с.

25. Бреббия К., Теллес Ж., Броубел Л. Методы граничных элементов. > -М.: Мир, 1987. -524 с.

26. Бреббия К., Уокер С. Применение метода граничных элементов в технике. -М.: Мир, 1982." -248 с.

27. Бублик Б.Н. Численное решение задач динамики пластин и оболочек. -К.: Наукова думка, 1969. -147 с.

28. Вазов В.Р., Форсайт Г.Е Разностные методы решения дифференциальных уравнений. -М.: ИЛ, 1969.

29. Васильев В.В. Механика конструкций из композиционных материалов. -М.: Машиностроение, 1988. -272 с.

30. Ван Тассел Д. Стиль, разработка, эффективность, отладка и испытание программ. -М.: Мир, 1985. -332 с.

31. Вайнберг Д.В., Ждан В.З. Матричные алгоритмы в теории оболочек. -К.: КГУ, 1967. -165 с.

32. Варвак П.М., Варвак Л.П. Метод сеток в задачах расчета строительных конструкций. -М.: Стройиздат, 1977. -159 с.

33. Вайнберг Д.В., Синявский А.Л. Расчет оболочек. -К.: Изд-во Гос-стройиздат УССР, 1961. -119 с.

34. Вахитов М.Б., Сафарнев М.С., Халилулин В.И. Расчет консольных пластин методом прямых // Труды КАИ,-Вып.166. -Казань, 1974. -С.52-61.

35. Вахитов М.Б. Интегрирующие матрицы аппарат численного решения дифференциальных уравнений строительной механики // Изв. ВУЗов. Авиационная техника, 1966. -№3. -С. 50-61.

36. Вахитов М.Б., Фирсов В.А. Численные методы решения одномерных задач строительной механики летательных аппаратов. -Казань: КАИ, 1985. -66 с.

37. Векуа И.Н. Некоторые общие методы построения различных вариантов теории оболочек. -М.: Наука, 1982. -286 с.

38. Виноградов А.Ю., Виноградов Ю.И. Метод решения краевых задач механики деформирования пластин и оболочек для дифференциальных уравнений только с четными производными // Докл. РАН, 1993. -Т.330. -№1. -С. 41-42.

39. Виноградов А.Ю. Численное решение задач механики деформирования пластин и оболочек на основе уравнений с четными производными // Труды XVI Международной конференции по теории оболочек и пластин. -Т.З. Н.Новгород: Изд-во НГУ, 1994. -С. 53-57.

40. Виноградов А.Ю., Виноградов Ю.И. Совершенствование метода прогонки Годунова для задач строительной механики // Изв. РАН Механика твердого тела, 1994. -№4. -С. 187-191.

41. Виноградов А.Ю. Вычисление начальных векторов для численного решения краевых задач // Журнал вычислительной математики и математической физики, 1995. -Т.З5. -№1. -С. 156-159.

42. Виноградов А.Ю. Дискретный аналог метода Гельфанда-Локуциевского для краевых задач теории оболочек и пластин // Труды XVI Международной конференции по теории оболочек и пластин, Казань: Изд-во КГТУ, 1995.

43. Виноградов А.Ю. Модификация метода Годунова // Труды Международной научно-технической конференции "Современные проблемы машиноведения", Гомель: ГПИ им. П.О.Сухого, 1996. -С. 39-41.

44. Виноградов. А.Ю. Методы переноса краевых условий для дифференциальных уравнений теории оболочек // Диссертация на соискание учетной степени к.ф. -м.н. -М., 1996. -163 с.

45. Виноградов А.Ю., Виноградов Ю.И. Функции Коши-Крылова и алгоритмы решения краевых задач теории оболочек // ДАН РФ, 2000. -Т.375.-№3.-С. 331-333.

46. Виноградов А.Ю. Перенос краевых условий для дифференциальных уравнений теории оболочек только с четными производными // Изв. РАН МТТ, 2001.

47. Виноградов А.Ю., Виноградов Ю.И., Гусев Ю.А, Клюев Ю.И. Перенос краевых условий функциями Коши-Крылова и его свойства // Изв. РАН МТТ, 2001.

48. Виноградов А.Ю. Дифференциально-разностный метод и его возможности // Тез. докладов 4-й Всероссийской научной конференции "Краевые задачи и математическое моделирование", Новокузнецк, 2001.

49. Виноградов Ю.И. Напряженно-деформированное состояние цилиндрической оболочки при сосредоточенном нагружении // Изв. ВУЗов Машиностроение. -М., 1973. -№11. -С. 5-9.

50. Виноградов Ю.И., Клюева Г.П. Матричный алгоритм численного метода решения задач о локальном и сосредоточенном нагружении оболочеквращения. -В кн.: Расчет тонкостенных элементов конструкций. Труды МВТУ. -1976. -№206. -С. 67-77.

51. Виноградов Ю.И. Жесткость цилиндрических оболочек при локальном нагружении. -В кн.: Жесткость машиностроительных конструкций: Тез. докладов Всесоюзной научно-технической конференции. -Брянск, 1976. -С. 47-51.

52. Виноградов Ю.И. Изгиб оболочки и методы его вычисления при сосредоточенном нагружении. -В кн.: Шестой Всесоюзный съезд по теоретической и прикладной механике: Тез. докладов. -Ташкент, 1986. -С. 167.

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

54. Виноградов Ю.И., Клюев Ю.И. Устойчивость и динамика конструкций при начальных напряжениях // Механика композиционных материалов и конструкций. ИПРИМ РАН. -1998. -Т.4. -№2. -С.49-60.

55. Виноградов Ю.И., Гусев Ю.А., Золотухин B.C. Специальный численный метод исследования концентрации напряжений в оболочках с контролируемой погрешностью // Тез. Докладов 3-й международной конференции "Авиация и космонавтика 2004", Москва, 2004. -С.54.

56. Виноградов Ю.И., Золотухин B.C. Метод исследования концентрации напряжений в тонкостенных элементах аэрокосмических систем // Тез. Докладов Международного симпозиума "Образование через науку", Москва: МГТУ им. Н.Э. Баумана, 2005. -С. 516.

57. Виноградов Ю.И., Петров В.И. Исследование локальной прочности тонкостенных конструкций методом Коши // Тез. Докладов 4-й международной конференции "Авиация и космонавтика 2005", Москва, 2005. -С. 161.

58. Виноградов Ю.И., Гусев Ю.А., Золотухин B.C. Методы исследования концентрации напряжений в оболочках // Вестник Московского авиационного института, 2005, т. 12, № 3. -С. 61-65.

59. Виноградов Ю.И. Петров В.И. Методы исследования концентрации напряжений в тонкостенных элементах аэрокосмических систем // Сб. докладов Международного симпозиума "Образование через науку", Москва: МГТУ им. Н.Э. Баумана, 2006. -С. 132-140.

60. Виноградов Ю.И. Гусев Ю.А., Петров В.И. Мультипликативный метод Коши исследования концентрации напряжений в оболочках // Тез.докладов 5-й международной конференции "Авиация и космонавтика 2006", Москва, 2006.-С. 110-111.

61. Виноградов Ю.И., Петров В.И. Простейшие мультипликативные алгоритмы исследования напряжений в оболочках // Тез. докладов XV Международной конференции по вычислительной механики и современным прикладным системам, Алушта 2007. -С 122.

62. Виноградов Ю.И., Петров В.И. Параметрические исследования концентрации напряжений в оболочках с контролируемой погрешностью // Тез. докладов 6-й международной конференции "Авиация и космонавтика -2007", Москва, 2007, -С 36.

63. Виноградов Ю.И., Петров В.И. Исследование концентрации напряжений в транспортно пусковом стакане // Труды 3-й Международной конференции "Ракетно-космическая техника: фундаментальные и прикладные проблемы", Москва, МГТУ им. Н.Э. Баумана, 2007. -С.

64. Виноградов Ю.И. Гусев. Ю.А. Никитенко В.В. Петров В.И. Локальная прочность тонкостенных конструкций // Известия РАН МТТ, 2007, №. -С 158-169.

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

66. Власов В.З. Строительная механика тонкостенных пространственных систем. -М.: Госстройиздат, 1949.

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

68. Власов В.З. Тонкостенные пространственные системы. -М.:Гостехиздат, 1958. -472 с.

69. Владимиров B.C. Уравнения математической физики. -М.: Наука, Гл. ред. физ.-мат. лит., 1988. -512 с.

70. Воеводин В.В., Кузнецов Ю.А. Матрицы и вычисления. -М.: Наука, 1984. -320 с.90."Гавурин М.К. Лекции по методам вычислений. -М.: Наука, 1971. -248 с.

71. Галлагер Р. Метод конечных элементов. Основы. -М.; Мир, 1984.428 с.

72. Гантмахер Ф.Р. Теория матриц. -М.: Наука, 1988. -548 с.

73. Гельфанд И.М., Локуциевский О.В. Метод "прогонки". Дополнение к книге С.К.Годунова, В.С.Рябенького Введение в теорию разностных схем. -М., Физматгиз, 1962. -С. 283-309.

74. Годунов С.К. О численном решении краевых задач для систем линейных обыкновенных дифференциальных уравнений // Успехи математических наук, 1961. -Т.16. -Вып.3.(99). -С.171-174.

75. Годунов С.К., Рябенький B.C. Разностные схемы. -М.: Наука, 1973. -400 с.

76. Гольденвейзер А.Л. Теория упругих тонких оболочек. -М.: Наука, 1976. -510 с.

77. Гольденвейзер А.Л. К вопросу о расчете оболочек на сосредоточенные силы // Прикладная математика и механика. -Т8. -№2, 1954. -С. 181-186.

78. Голованов А.И. Конечные элементы тонких не пологих оболочек. Способы построения//Прикладные проблемы прочности и пластичности. -Н.Новгород, 1991. -С.58-65.

79. Голованов А.И. Расчет составных оболочек произвольной геометрии // Проблемы механики оболочек. -Калинин,1988. -С.33-40.

80. Голованов А.И., Корнишин М.С. Введение в метод конечных элементов статики тонких оболочек. -Казань: Казанский физико-технический институт, 1989. -269 с.

81. Green G. An assay on the application of mathematical anlusis to the theory of electricity and magnetism/ZNottingham, 1828.4

82. Григолюк Э.И., Кабанов B.B. Устойчивость оболочек. -М.: Наука, 1978. -360 с.

83. Григолюк Э.И., Шалашилин В.И. Проблемы нелинейного деформирования: Метод продолжения решения по параметру в нелинейных задачах механики твердого деформируемого тела. -М.: Наука, 1988. -231 с.

84. Григолюк Э.И., Толкачев В.М. О расчете цилиндрических оболочек, загруженных по линиям//Прикладная математика и механика, 1967. -Т.31. -Вып.6. -С.1141-1146.

85. Григоренко Я.М. Изотропные и анизотропные слоистые оболочки вращения переменной жесткости. -К.: Наукова^ думка, 1973. -288 с.

86. Григоренко Я.М., Василенко А.Т. Задачи статики анизотропных неоднородных оболочек. -М.: Наука, 1992. -336 с.

87. Григоренко Я.М, Василенко А.Т., Голуб Г.П. Статика анизотропных оболочек с конечной сдвиговой жесткостью. -К.: Нукова думка, 1987. -216 с.

88. Григоренко Я.М., Крюков H.H. Решение линейных и нелинейных краевых задач теории пластин и оболочек на основе метода линий// Прикладная механика, 1993. -№4. -С.3-11.

89. Григоренко Я.М., Василенко А.Т., Беспалова Е.И. и др. Численное решение краевых задач статики ортотропных оболочек с переменными параметрами. -К.: Наукова думка, 1975. -183 с.

90. Гусев. Ю.А. Мультипликативные алгоритмы переноса краевых условий в задачах механики деформирования оболочек // Диссертация на соискание ученой степени к.ф.-м.н. -М., 2003. -176 с.

91. Гурьянов Н.Г. Замкнутая цилиндрическая оболочка под действием сосредоточенной силы // Исследования по теории пластин и оболочек. -Казань: Изд-во Казанского государственного университета, 1970. -№6. -С.29-36.

92. Dabrowski Otton, Szmigielski Roman. Solution of shallow shells by boundary element method. Problem of corners. -Mrch. teor.: Stosow. -1988, 26, №4. -P.603-610.

93. Даревский B.M. Определение перемещений и напряжений в цилиндрической оболочке при локальных нагрузках. -В сб. Прочность и динамика авиационных двигателей. -М.: Машиностроение, 1964. -Вып.1. -С.23-83.

94. Демидович Б.П., Марон И.А. Основы вычислительной математики. -М.: Физматгиз, 1966. -664 с.

95. Демьянович Ю.С. К вопросу об изгибе цилиндрической оболочки сосредоточенной силой. -В сб. Исследование по упругости и пластичности. -Л.: Изд-во ЛГУ, 1968. -Вып.2. -С.121-131.

96. Джордж А., Лю Дж. Численное решение больших разреженных систем уравнений. -М.: Мир, 1984. -309 с.

97. Du I., Huí D. Frequency-load interaction of imperfect angle-ply cylindrical panels under compression and pressure// AIAA Journal, 1977. -25. -13. -P.484-491.

98. Elsbemd G.F., Leissa A.W. The vibration of non-circular cylindrical shells with initial stresses//J. Sound and Vibr., 1973/ -29. №3. -P.309-329.

99. Жигалко Ю:П. Расчет тонких упругих цилиндрических оболочек на локальные нагрузки (обзор литературы, метод и результаты). -В сб. Исследования по теории пластин и оболочек. -Казань: Изд-во Казанского университета, 1966. -Вып.4. -С.3-41.

100. Зенкевич О., Чанг И. Метод конечных элементов в теории сооружений и механики сплошных сред. -М.: Недра, 1974. -240 с.

101. Зенкевич О. Метод конечных элементов в технике: Перевод с английского. -М.: Мир, 1975. -544 с.

102. Зенкевич О., Морган К. Конечные элементы и аппроксимация. -М.: Мир, 1986. -318 с.

103. Зибаров А.В., Бабаев Д.В., Миронов А.А. Система визуализации Mirage // Тезисы докладов XI международной конференции по вычислительной механике и современны прикладным программным системам, Москва-Истра, 2001. С. 199-200.

104. Зибаров А.В., Бабаев Д.Б., Миронов А.А., Синютин А.Н., Лохов

105. B.C. Технологии распараллеливания и распределения вычислений//Тезисы докладов XI международной конференции по вычислительной механике и современным прикладным программным системам, Москва-Истра, 2001.1. C.201-202.

106. Ильгамов М.А., Иванов В.А., Гулин Б.В. Прочность, устойчивость и динамика оболочек с упругим заполнителем. -М.: Наука, 1977. -351 с.

107. Калиткин Н.Н. Численные методы. -М.: Наука, 1978. -512 с.

108. Kamenskii G.A., Myshkis A.D. On the mixed type functional-differential equations // Nonlinear Analysis. TMA. -1887/ -V.30-, No.5, pp.25772584.

109. Канторович Л.В., Крылов В.И. Приближенные методы высшего анализа. -М., -Л.: Физматгиз, 1962. -695 с.

110. Кармишин А.В., Лясковец В .А., Мяченков В.И., Фролов А.Н. Статика- и динамика тонкостенных обол очечных конструкций. -М.: Машиностроение, 1975. -375 с.

111. Conte S.D. The numerical solution of linear boundary vaiue problems. -Siam Review, 1966. -Vol.8.

112. Коноплев Ю.Г. Устойчивость цилиндрических панелей. -Казань: Изд-во Казанского государственного университета, 1976. -20 с.

113. Корн Г., Корн Т. Справочник по математике. Для научных работников и инженеров. -М.: Наука, 1984. -831 с.

114. Котов Е.А. Построение базисных функций для решения осесим-метричных задач нелинейной теории оболочек вращения методом Бубнова-Галеркина. -В кн.: Прочность материалов и элементов судовых конструкций. -Л.: Наука, 1985.-С. 50-57.

115. Крауч С., Старфилд С. Методы граничных элементов в механике твердого тела. -М.: Мир, 1987. -328 с.

116. Крылов А.Н. О расчете балок, лежащих на упругом основании. -Л.: АН СССР, 1931.-154 с.

117. Крылов В.И., Бобков В.В., Монастырский П.И. Вычислительные методы. -М.: Наука, 1976. -Т.1. -304 с.

118. Cruse Т.A., Rizzo F.J. A direct formulation and numerical solution of the general transient elastodynamic problem. -Int. J. Math. Anal. Appics, 1968, 22. -P.244-259.

119. Cruse T.A., Rizzo F.J. (eds.) Boundary integral equations method: computational applications in applied mechanics. ADM 11. -New York: Am Soc. Mech/Engrs, 1975.

120. Кулаков B.M., Толкачев B.M. Изгиб пластины произвольного очертания // ДАН СССР, 1976. -Т.230, -Ч. -С.56-59.

121. Купрадзе В.Д. Методы потенциала в теории упругости, -М.: Физматгиз, 1963. -472 с.

122. Куранов Б.А., Гусев С.С. Применение метода суперэлементов для расчета сложных машиностроительных конструкций/ТРасчеты на прочность. -М.: Машиностроение, 1983. -Вып.26. -С.174-182.

123. Курант Р., Фридрихе, Леви Г. О разностных уравнениях математической физики//Успехи математических наук, 1940. -Вып.8. -С. 112-125.

124. Леонов В.И., Хазанов Х.С. Напряженное состояние цилиндрической оболочки, нагруженной через круглую площадку нормальными поверхностными силами // Труды КУАИ, вып. 60, 1973.

125. Ли C.B., Пиан Т.Х. Усовершенствование метода расчета конечных элементов для пластинок и оболочек с помощью смешанного подхода//РТК, 1978. -Т.16. -№1. -С.38-45.

126. Лукасевич С. Локальные нагрузки в пластинах и оболочках. -М.: Мир, 1982. -542 с.

127. Loo Wen-da, Gao Shi-qiao. The effect local geometric imperfections of rotational shell on its natural frequencies and models// Appl. Math. And Mech., 1987.-Vol.8.-№11. -РЛ013-1018.

128. Мак-Кракен Д., Дорн У. Численные методы и программирование на ФОРТРАНе. -М. : Мир, 1977.-584с.

129. Мальцев А.И. Основы линейной алгебры. -М.: Наука, 1975. -400 с.

130. Мартыненко С.И. Универсальная многосеточная технология для численного решения дифференциальных уравнений в частных производных на структурированных сетках// Вычислительные методы и- программирование, 2000. -Т.1, раздел 1. -С.85-104.

131. Мартыненко С.И. Программное обеспечение для универсальной многосеточной технологии: строительные блоки и диагностические инструменты// Вычислительные методы и программирование, 2001. -Т.2, раздел 4. -С.1-6.

132. Марчук Г.И. Методы вычислительной математики. -М.: Наука, 1973. -352 с.

133. Марчук Г.И. Методы вычислительной математики. -М.: Наука, 1989. -608 с.

134. Марчук Г.И., Агошков В.И. Введение в проекционно-сеточные методы. -М.: Наука, 1981. -416 с.

135. Masataka Т. Introduction to boundary element method/ /J. Jap. Soc. Simulat. Technol. -1987. -6. -l2. -P.77-83.

136. Masinda J. Application of the boundary element method to elasticity and thermoelasticity problems. "Monogr. and Mem. Nat. Res. Inst. Mach. Des/ Praha-Bechavice", 1986, '36. -P.58.

137. Метод конечных элементов в механике твердых тел/Под общей редакцией А.С.Сахарова и И.Альтенбаха. -Киев: Вища школа, головное изд-во, 1982.

138. Методы расчета оболочек. -Т.4. Теория оболочек переменной жесткости/ Под редакцией Григоренко Я.М., Василенко А.Т. -К.: Наукова думка, 1981.-544 с.

139. Метод граничных интегральных уравнений. Вычислительные аспекты и приложения в механике. -Серия: Механика. Новое в зарубежной науке. -М.: Мир, 1978. -210 с.

140. Микеладзе Ш.Е. Новые методы интегрирования дифференциальных уравнений. -М.: ГТИИ, 1951.-291 с.

141. Митчелл Э., Уэйт Р. Метод конечных элементов для уравнений с частными производными. -М.: Мир, 1981. -216 с.

142. Mitchell S.A., Vavasis S.A. Quality Mesh Generation in Higher Dimension. 12th ACM Symposium on Computational Geometry, 1996, pp.48-57.

143. Михлин С.Г. Численная реализация вариационных методов. -М.: Наука, 1966.-432 с.

144. Михлин С.Г. Вариационные методы в математической физике. -М.: Наука, 1970. -512 с.

145. Моссаковский В.И., Макаренков А.Г., Никитин П.И. и др. Прочность решетных конструкций. -М.: Высшая школа, 1990. -359 с.

146. Мяченков В.И., Григорьев И.В. Расчет составных оболочечных конструкций на ЭВМ: Справочник. -М.: Машиностроение, 1981. -216 с.

147. Мяченков В.И., Мальцев В.П. Методы и алгоритмы расчета пространственных конструкций на ЭВМ ЕС. -М.: Машиностроение, 1984. -280 с.

148. Нерубайло Б.В. Локальные задачи прочности цилиндрических оболочек. -М.: Машиностроение, 1983. -248 с.

149. Нерубайло Б.В. Применение вариационного метода В.З.Власова к расчету круговой цилиндрической оболочки на местную нагрузку//Изв. ВУЗов Авиационная техника, 1967. -№1. -С.45-50.

150. Нерубайло Б.В., Сибиряков В.А. Напряжения в круговой цилиндрической оболочке при радиальном локальном воздействии. -В сб. Прочность и устойчивость авиационных конструкций.//Труды МАИ, 1971. -Вып. 180. -С.25-44.

151. Новожилов В.В. Теория тонких оболочек. -Л.: Судпромгиз, 1962. -432 с.

152. Новожилов В.В. Вопросы механики сплошной среды. -Л.: Судостроение, 1989. -400 с.

153. Новожилов В.В., Черных К.Ф., Михайловский Е.И. Линейная теория тонких оболочек. -Л.: Политехника, 1991. -656 с.

154. Noor А.К., Anderssen С.М. Mixed models and reduced selective Integration displacement models for nonlinear shell analysis// Proc/ ASME winter annual meeting/ -Washington: AMD, 1981. -Vol.48. -P.l 19-146.

155. Образцов И.Ф., Булычев Л.А., Васильев B.B., Елпатьевский А.Н., Жеков К.А., Иванов Ю.И., Коновалов Б.А., Матюшев Ю.С., Шклярчук Ф.Н. Строительная механика летательных аппаратов. -М.: Машиностроение, 1986. -536 с.

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

157. Оден Д. Конечные элементы в нелинейной механике сплошных сред. -М.: Мир, 1976. -464 с.

158. Oden J.T., Carey G.F. Finite elements. Mathematical aspects. -Englwood Cliffs, N.Y.Prentice-Hall, -1983. vol.4. -195 p.

159. Ортега Дж. Введение в параллельные и векторные методы решения линейных систем. -М.: Мир, 1991.-364 с.

160. Ортега Дж., Пул У. Введение в численные методы решения дифференциальных уравнений. -М.: Наука, 1986. -288 с.

161. Папкович П.Ф. Строительная механика корабля, ч.21 -Л.: Судпром-гиз, 1941.-960 с.

162. Паймушин В.Н. К вариационным методам решения нелинейных пространственных задач сопряжения деформируемых тел//ДАН СССР, 1983. -Т.273, -№5. -С.1083-1086.

163. Паймушин В.Н., Андреев C.B. К нелинейной теории трехслойных оболочек со слоями переменной толщины и сложной геометрии. -Казань: Изд-во КГУ, 1981. -Вып. 16. -С.29-36."

164. Паймушин В.Н., Фирсов В.А. Оболочки из стекла. Расчет напряженно-деформированного состояния. -М.: Машиностроение, 1993. -208 с.

165. Паймушин В.Н., Сидоров И.Н. Вариант метода граничных интегральных уравнений для решения задач статики изотропных оболочек произвольной геометрии//Изв. АН СССР Механика твердого тела, 1991. -№11. С.160-169.

166. Паймушин В.Н., Сидоров И.Н., Сулейманов И.М. Построение граничных интегральных уравнений различных вариантов теории оболочексложной геометрии//Изв. РАН Механика твердого тела, 1997. -№4 -С.133-143.

167. Петров.В.И. Исследование концентрации напряжений в оболочках с априори контролируемой погрешностью // Материалы XII Международного симпозиума "Динамические и технологические проблемы механики конструкций и сплошных сред", Москва, 2006. -С. 256-258.

168. Петров В.И. Концентрация напряжений в тонкостенных конструкциях // Материалы XIII Международного симпозиума "Динамические и технологические проблемы механики конструкций и сплошных сред", Москва 2007. -С.

169. Победря Б.Е. Численные методы в теории упругости и пластичности. -М.: Изд-во Московского университета, 1981. -344 с.

170. Попов Г.Я. Метод ортогональных многочленов в контактных задачах теории упругости//ПММ, 1969. -Т.ЗЗ. -Вып.З. -С.543-547.

171. Постнов В.А., Хархурин И.Я. Метод конечных элементов в расчетах судовых конструкций. -Л.: Судостроение, 1974. -342 с.

172. Постнов В.А. Численные методы расчета судовых конструкций. -Л.: Судостроение, 1977. -280 с.

173. Постнов В.А., Родионов A.A. Метод суперэлементов в расчете инженерных сооружений. -Л.: Судостроение, 1979. -288 с.

174. Постнов В.А., Суслов В.П. Строительная механика корабля и теория упругости. В 2-х т. -Л.: Судостроение, 1987. -Т.1: Теория упругости и численные методы решения задач строительной механики корабля. -288 с.

175. Постнов В.А., Ростовцев Д.М., Суслов В.П., Кочанов Ю.П. в 2-х т. -Л.: Судостроение, 1987. -Т.2: Изгиб и устойчивость стержней, стержневых систем, пластин и оболочек. -416 с.

176. Прокопович М.Е., Слезингер И.Н., Штейнберг М.В. Расчет цилиндрических оболочек. -К.: Будщивильник, 1967. -240 с.

177. Пшеничнов Г.И. Метод декомпозиций решения уравнений и краевых задач//ДАН СССР, 1985. -Т.282. -№4. -С.792-794.

178. Работнов Ю.Н. Изгиб цилиндрической оболочки сосредоточенной силой// Докл. АН СССР, 1946. -Т.52. -№4. -С.110-114.

179. Расчеты машиностроительных конструкций методом конечных элементов: Справочник/ В.И.Мяченков, В.П.Мальцев, В.П.Майборода и др.; Под общей редакцией В.И.Мяченкова. -М.: Машиновтроение, 1989. -520 с.

180. Рвачев B.JI. Об удовлетворении краевым условиям в методе Буб-нова-Галеркина с помощью Я-функций//Сб. статей Успехи механики деформируемых сред. -М.; Наука, 1975. -С.488-501.

181. Reddy T.N., Chao W.C. Large-deflection and large-amplitude free vibrations of laminated composite-material plates//Comput. and struct., 1981. -43. -P.341-342.

182. Ректорис К. Вариационные методы в математической физике и технике. -М.: Мир, 1985. -589 с.

183. Рикардс Р.Б. Метод конечных элементов в теории оболочек и пластин. -Рига: Зинатне, 1988. -284 с.

184. Ritz W.H. J. Reine angew. Math., 1908. -T.135; Ann Physik, IV, 28, 737, 1909; Gesamm Werke, 1941. -S.228.

185. Рихтмайер P.Д. Разностные методы решения краевых задач. -М.: ИЛ, 1960.

186. Рихтмайер Р., Мортон К. Разностные методы решения краевых задач. -М.: Мир, 1972.-418 с.

187. Розин Л.А. О связи метода конечных элементов с методами Бубнова- Галеркина и Ритца// Строительная механика сооружений. -Л., 1971. -С.6-27.

188. Розин Л.А. Метод конечных элементов в применении к упругим системам. М.: Стройиздат, 1977. -129 с.- 189220. Рояк Д.А. Расчет весьма длинной цилиндрической оболочки по методу наложения // Инженерный журнал, Т.4, вып.1, 1964.

189. Сайтов И.Х. Метод декомпозиции оболочечных систем применительно к интегральной формулировке решения фрагмента//Теория и методы исследования пластин и оболочек сложной формы: Межвуз. сб. -Казань: КАИ, 1987. -С.69-75.

190. Сайтов И.Х., Ларионов Н.Г. Матричная форма интегрально-разностного метода в решении двумерных краевых задач теории пологих оболочек // Вопросы прочности тонкостенных авиационных конструкций. Межвуз сб. -Казань, КАИ, 1987. -С. 109-115.

191. Самарский A.A. Введение в теорию разностных схем. -М.: Наука, 1971.-550 с.

192. Самарский A.A., Попов Ю.П. Разностные схемы газовой динамики. -М.: Наука, 1975.- 190228. Самарский A.A., Николаев Е.С. Методы решений сеточных уравнений. -М.: Наука, 1978. -592 с.

193. Самарский A.A., Гулин A.B. Численные методы. -М.: Наука, 1989. -432 с.

194. Самарский A.A. Теория разностных схем -М.: Наука, 1989. -616 с.

195. Самарский A.A. Всесоюзная школа молодых ученых. Теория и прикладные проблемы вычислительной математики и математической физики // Ж.ВМиМФ, 1983. -Т.23. -№1. -С.246-247.

196. Свирский И.В. Методы типа Бубнова-Галеркина и последовательных приближений. -М.: Наука, 1968. -200 с.

197. Skubachevskii A.L. Elliptic functional differential equations and applications. Basel- Boston-Berlin: Birkhauser, 1997.

198. Слободянский М.Г. Способ приближенного интегрирования уравнений с частными производными и его применение к задачам теории упругости. -ПММ, 1939, т.З, вып. 1. -С. 75-82.

199. Стренг Г., Фикс Дж. Теория метода конечных элементов. -М.:Мир,1877. -349 с.

200. Стеклов В.А. Основы теории интегрирования обыкновенных дифференциальных уравнений. -М., -JL: Гос.издательство, 1927. -419 с.

201. Суперэлементный расчет подкрепленных оболочек/3.И.Бурман, О.М.Аксенов, В.И. Лукашенко, М.Т.Тимофеев. -М.: Машиностроение, 1982. -256 с.

202. Терегулов И.Г. К построению уточненных теорий пластин и оболочек // Прикладная математика и механика, 1962. -Т.26. -№2. -С.346-358.

203. Тимошенко С.П. Устойчивость стержней, пластин и оболочек. М.: Наука, 1971.-807 с.

204. Ting L., Yuan S.W. On radial deflection of a cylinder of finite length with various end conditions. Journal of the Aeronautical Sciences, 1958. -Vol.25. -4. -P.230-234.

205. Толкачев В.М. Метод компенсирующих нагрузок в теории изгиба пластин // Изв. АН СССР Механика твердого тела, 1988. -№1. -С.155-160.

206. Угодчиков А.Г., Хуторянский Н.М. Метод граничных элементов в механике деформированного твердого тела. -Казань: Изд-во Казанского университета, 1986. -295 с.

207. Федорюк М.В. Асимптотические методы для линейных обыкновенных дифференциальных уравнений. -М.: Наука, 1983. -352 с.

208. Федосьев В.И. Сопротивление материалов. -М.: Наука, 1974. -559 с.

209. Филин А.П. Элементы теории оболочек. -JL: Стройиздат, 1975. -256 с.

210. Филин А.П. Прикладная механика твердого деформируемого тела, Т.П. -М.: Наука, 1978. -616 с.

211. Флюгге В. Статика и динамика оболочек. -М.: Госстройиздат, 1963.-306 с.

212. Флетчер К. Численные методы на основе метода Галеркина. -М.: Мир, 1988. -352 с.

213. Хижняк В.К., Шевченко В.П. Действие сосредоточенных сил на анизотропные оболочки// Изв. АН СССР Механика твердого тела, 1972. -№4. -С.123-128.

214. Халилулин В.И., Рахматкулов Н.У. К расчету монолитных крыльев с -произвольно расположенными подкрепляющими элементами//Вопросы прочности и долговечности элементов авиационных конструкций: Межвуз. сб., Куйбышев, 1983. -С.3-10.

215. Holt P.J., Webber J.P. Finite element for honeycomb sandwich plates and shells//Aeron. J., 1980. -Vol.84. -№831. -P.l 13-123; №832. -P.157-167.

216. Хорн P., Джонсон Ч. Матричный анализ. -M.: Мир, 1989. -655 с.

217. Хьюз Д., Мичтом Д. Структурный подход к программированию. -М.: Мир, 1980. -278 с.

218. Чернышев Г.Н. О действии сосредоточенных сил и моментов на упругую тонкую оболочку произвольного очертания/ЛТрикладная математика и механика, 1963. -Т.27. -Вып.1. -С. 126-134.

219. Чернышев Г.Н. Расчет сферических оболочек на действие сосредоточенных сил//Изв. АН СССР Механика, 1963. -№1. -С.99-108.

220. Чернышев Г.Н. Асимптотические методы в теории оболочек (Сосредоточенные нагрузки)//Труды VI Всесоюзной конференции по теории оболочек и пластинок, Москва. -С.799-810.

221. Чернышев Г.Н. Прогиб под сосредоточенной силой в оболочках положительной кривизны//Прикладная математика и механика, 1967. -Т.31. -№5, -С.883-886.

222. Чернина B.C. Статика тонкостенных оболочек вращения. -М.: Наука, 1968. -455 с.

223. Шевляков Ю.А., Шевченко В.П. Пологая сферическая оболочка под действием сосредоточенных сил и моментов//Прикладная механика,1965. -Вып.2. -С.74-77.

224. Шевляков Ю.А., Шевченко В.П. К вопросу о действии сосредоточенных воздействий на пологие оболочки. -В сб. Концентрация напряжений. -К.: Наукова думка, 1965. -Вып.1. -С.326-337.

225. Шевляков Ю.А., Шевченко В.П. О действии сосредоточенных сил и изгибающих моментов на пологую цилиндрическую оболочку//Прикладная механика, 1966, -Т.2. -Вып.1. -С.120-123.

226. Шевченко В.П. К решению граничных задач круговых цилиндрических оболочек при локальных нагрузках. -В кн. Теория оболочек и пластин. -М.: Наука, 1973. -С.763-767.

227. Шевченко В.П. К вопросу о действии сосредоточенных моментов на упругие тонкие оболочки произвольной кривизны. -В сб. Теоретическая и прикладная механика. -Киев-Донецк: Вища школа, 1977. -Вып.8. -С.42-47.