Вариационно-разностный метод решения нестационарных волновых задач динамики тонкостенных конструкций тема автореферата и диссертации по механике, 01.02.06 ВАК РФ

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

1.1.Задачи динамики упругопластических сред

1.2.3адачи динамики оболочек. Модель Тимошенко

1.3.Дисперсионные свойства уравнений теории пластин Тимошенко

2. Вариационно-разностный метод

3. Преобразование вариационно-разностных и конечно-элементных 72 численных схем к виду конечно-разностнЫх

3.1.Сеточный аналог формул интегрирования по частям. Конечно- 72 разностное представление вариационно-разностных схем.

3.2.Конечно-разностное представление схем МКЭ

3.3.Примеры преобразования схем МКЭ в конечно- разностные

3.4.Автоматическое построение конечно-разностного представления 97 схем МКЭ. Алгоритм преобразования и программная реализация

3.5.Результаты работы программы построения конечно-разностного 99 представления схем МКЭ

4^. Анализ аппроксимации задач теории пластин типа Тимошенко по 103 пространственным переменным» Проблема малого параметра и крупных сеток

4.1 .Вывод сеточных уравнений вариационно-разностных и конечно- 103 К, элементных схем теории пластин типа Тимошенко

4.2.Анализ численных схем решения одномерных задач теории 113 пластин

4.3.Индексная коммутативность численного дифференцирования 122 4.4.Эквивалентные преобразования разностных схем

4.5.Анализ и тестирование численных схем решения двумерных 133 задач теории пластин и оболочек

4.6.«Ажурные» схемы метода конечного элемента

Т 5. Устойчивость численных схем

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

5.2.Оценки устойчивости разностных схем решения трехмерной за- 157 дачи теории упругости

5.3.Оценки устойчивости одномерных схем теории пластин 164 Тимошенко

5.4.0ценки устойчивости двумерных схем теории пластин 172 Тимошенко

5.5.Неустойчивость типа «песочные часы»

5.6.Граничная неустойчивость численных схем решения задач 177 трехмерной теории упругости —

6. Повышение эффективности численных схем решения задач динамики 183 конструкций. Явно-неявные схемы со стабилизирующим оператором

6.1 .Регуляризация численных схем теории пластин и оболочек

6.2. Разностная схема решения плоской задачи теории упругости с 190 неявным стабилизирующим оператором

6.3.Схемы со стабилизирующим оператором с точки зрения МКЭ

6.4.Регуляризация численных схем решения трехмерной задачи 196 теории упругости

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

7. Численный анализ нелинейных процессов деформирования 211 тонкостенных конструкций при импульсных воздействиях

7.1 .Выпучивание цилиндрической оболочки при неосесимметричном 211 продольном ударе

7.2.0сесимметричное выпучивание сферических куполов под 223 действием импульса давления:

7.3.Выпучивание пологих сферических куполов, квадратных в плане, 229 под действием импульса давления

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

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

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

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

Развитие численных методов решения задач механики сплошных сред тесно связано с прогрессом в вычислительной математике. Основные достижения вычислительной математики в области решения задач математической физики отражены в учебниках и монографиях [19,64,65,74,94,95, 181,187,193-196, 205, 210, 222,223] и других. В настоящее время наиболее актуальной остается проблема эффективности методов, т. е. разработки численных схем, оптимальных по быстродействию. Причина этого - во все возрастающей сложности задач, решаемых численными методами. И хотя рост быстродействия ЭВМ решает многие из проблем, сравнимый с ними эффект дает совершенствование численных методов [65]. Рассмотрим более подробно современное состояние численных методов решения одного из сложных классов задач механики деформируемого твердого тела - нестационарной динамики тонкостенных конструкций.

Характеристика нестационарных динамических процессов в тонкостенных конструкциях. Расчет тонкостенных оболочечных и стержневых конструкций под действием нестационарных динамических нагрузок - задача, обладающая рядом специфических особенностей [11-14,71,72,129,158,162,163]. При решении задач нестационарной динамики конструкций необходимо учитывать следующие факторы: сложность геометрии конструкций, что приводит к необходимости использования неравномерных сеток и криволинейных систем координат, существенно усложняющих системы уравнений; геометрическая нелинейность (поскольку в процессе деформирования конструкций перемещения, с как правило, нельзя считать малыми); физическая нелинейность, обусловленная переходом к нелинейному деформированию материала; вырождение задач, связанное с малостью размеров по одной или двум координатам. Отметим, что нелинейность в задачах динамики конструкций носит несколько иной характер, чем, например, в газовой динамике. Так, в стержнях и оболочках практически невозможно возникновение ударных волн по причине малой толщины. В связи с вышесказанным построение численных схем решения задач данного класса требует серьезного анализа. Неудивительно, что этим вопросам посвящено большое количество работ как в отечественной, так и в зарубежной научной литературе [ 1,2,71,72,89-91,97,105-107,120,139,159,176,182,183,198,202,221]. Отметим также, что многие проблемы построения численных схем метода конечного элемента впервые проявились и в дальнейшем нашли свое решение именно при решении задач, связанных с расчетом тонкостенных конструкций [71,72,97,120,229-233,240,243,245,248,249].

Постановки задач. Математические модели и уравнения. Задачи нестационарной динамики конструкций допускают следующие постановки: начально-краевые задачи для систем дифференциальных уравнений в частных производных; системы законов сохранения в виде интегральных соотношений; вариационные задачи; системы граничных интегральных уравнений. Поскольку последние не получили широкого распространения при расчете тонкостенных конструкций, далее о них упоминать не будем. В ряде случаев численные схемы строятся не на базе континуальных задач, а на основе непосредственного применения физических законов, как это было с первыми схемами метода конечного элемента. При последнем подходе область разбивается на составные части, а далее уравнения неразрывности и движения этих частей составляются на основе физических соображений (см. исторический обзор в [172]). Пример такого подхода приводится, например, в [139,190].

При расчете тонкостенных конструкций используются модели сплошных сред и теории оболочек [10,84,86,92,93,114,133,165,166,175,203]. Учет геометрически нелинейных эффектов производится с использованием нелинейных лагранжевых, эйлеровых или совместных лагранжево-эйлеровых постановок [169]. Учет нелинейного поведения материала при деформировании осуществляется на базе теорий пластичности - деформационных и инкрементальных [115-119,124,133]. Ввиду наличия вырожденной координаты расчеты непосредственно по трехмерной теории зачастую являются неэффективными, что приводит к необходимости выделять процессы по вырожденной координате особо (вводить оболочечные модели), либо использовать сетки, адаптированные к данным задачам (с вытянутыми или сплющенными ячейками). При этом не имеет принципиального значения, вводятся ли гипотезы теории оболочек при выводе уравнений или на этапе построения численной схемы, как это делается в схемах МКЭ [172]. Важен характер принимаемых гипотез. В зависимости от этого приходим к вырожденной или вырождающейся задаче [71,72].

При дифференциальной записи основных уравнений могут использоваться разные формы. В методе конечных разностей наиболее распространена запись в виде уравнений движения замыкаемых соотношениями для деформаций через перемещения и определяющими соотношениями между напряжениями и деформациями, а при необходимости и уравнениями термодинамики. В характеристических схемах используются соотношения на характеристиках [78,127,140-143,149,220}, в схеме Годунова система записывается в дивергентной форме [ 18,95,188].

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

F(4t>Pij>qi.4i>&li,ePij) = Q>Pij =dqt/dxj v где qi - варьируемые параметры (в зависимости от выбранного вариационного принципа это могут быть перемещения, скорости, ускорения, напряжения и т.п. [4,81]). Данная форма записи служит основой для построения широкого класса - вариационно-разностных схем и схем метода конечного элемента {см. [83,156, 181] и др.).

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

Сеточно-характеристические методы развиты в работах В.Н.Кукуджанова, П.З.Лугового и других. Эти методы [76,107,134,140-143,146,149,220] базируются на записи исходных уравнений в виде соотношений на характеристиках. Интегрирование вдоль характеристик обладает следующим преимуществом: близость областей зависимости дифференциальной и сеточной задач. Они нашли достаточно широкое распространение при решении нестационарных волновых задач динамики массивных тел. В то же время при расчете тонкостенных конструкций применимость этих методов ограничена ввиду малости толщины конструкции и вследствие этого - необходимости учета отражения волн (и характеристик) от свободных поверхностей оболочек. В связи с этим возникают очевидные ограничения как на пространственный, так и на временной шаги интегрирования. Поэтому они нашли лишь ограниченное применение при решении задач этого класса [107,131,151]. В этом плане показателен пример монографии [159], авторы которой используют для расчета массивных тел се-точно-характеристическую схему, а при переходе к задачам теории оболочек переключаются на конечно-разностный подход. С другой стороны отметим попытки усовершенствования данного подхода при решении задач теории оболочек на базе модели Тимошенко [107]. Отметим также меньшую алгоритмическую простоту и универсальность данных методов по сравнению с методами конечных разностей и МКЭ (по крайней мере, в многомерном случае). В частности, это касается широкого класса прикладных задач, требующих учета сложной геометрии, неоднородностей материала, геометрической и физической нелинейности и ряда других факторов.

Схема С.К. Годунова. Данная схема [95] использует дивергентную запись уравнений. Следует отметить, что не все уравнения механики сплошных сред допускают запись в дивергентной форме - например, уравнения Прандтля-Рейсса (см. [191]) и теории оболочек. В таких случаях недивергентные слагаемые обычно переносятся в правую часть уравнений [18]. Важным свойством схемы С.К.Годунова является ее монотонность. Однако она обладает большой схемной вязкостью, которая исчезает лишь при числе Куранта к=1 , недостижимом для данной схемы за исключением одномерного случая. Модификации схемы С.К.Годунова, повышающие порядок ее аппроксимации до второго на гладких решениях(схема Лакса-Вендрофа и гибридные схемы) [222] уже не обладают свойством монотонности или имеют смешанный порядок (второй на монотонных решениях и первый - на немонотонных.

Схемы описанных выше классов-оправдывают себя при решении задач на очень коротких промежутках времени и интенсивных переходных процессах (задачи типа пробития, разрушения при ударе и т.п.). Однако эти задачи составляют лишь небольшую и не очень типичную часть задач динамики тонкостенных конструкций. У большинства задач динамики оболочек при характерных временах процесса порядка десяти и более пробегов волн по толщине картина волнового процесса отличается рядом специфических черт» Такие процессы носят волновой характер лишь вдоль срединной поверхности оболочки, а по толщине их можно считать близкими к квазистатическим после "многократного отражения волн от свободных поверхностей. При решении таких задач с помощью характеристических схем и схемы С.К.Годунова эта специфика ' не учитывается. Таким-образом* на задачах этого класса они становятся менее пригодными по сравнению с конечно-разностными и конечноэлементными схемами.

Конечно-разностные методы разрабатывались многими учеными. Большой вклад в их развитие внесли К.И.Бабенко, А.А.Самарский, Г.И.Марчук, Ю И-Шокин, В.С.Рябенький, Р.Рихтмайер, М.Уилкинс и другие ученые. Эти методы являются наиболее универсальными и гибкими при решении любых. задач математической физики. Они имеют наиболее широкое применение при решении задач математической физики [19,65,79,96,130,132,145,149,153,154, 169, 181,187,193-196,202,207, 222,223] и др. Вместе с тем они дают в руки исследователя широчайший произвол при своем построении. Следовательно, важнейшей задачей при построении конечноразностных схем является не расширение, а сужение их класса путем наложения на них различных дополни-. тельных ограничений. Классических критериев качества численных схем (порядок аппроксимации и устойчивость) явно недостаточно. Ниже рассмотрим ряд таких дополнительных требований к разностным схемам. Пока же отметим

N * лишь многолетнюю тенденцию дополнения количественных характеристик (порядок аппроксимации) качественными (консервативность, монотонность и т. д.). Наиболее распространенной для данных схем исходной формой записи является система уравнений движения, дополненная выражениями деформаций через перемещения и определяющими соотношениями. При этом утвердилась .определенная конструкция задания величин в узлах и "ячейках", получившая -■ наибольшее распространение после работ М.Уилкинса [206,207]. В ней имеем дело с разнесенными сетками. По пространству: узлы основной сетки, в которых вычисляются неизвестные перемещения, усилия, скорости и ускорения; "центры ячеек" основной сетки (по сути - узлы дополнительной сетки, смещенной относительно основной на некоторую долю ее шага), в которых вычисляются все величины, являющиеся первыми производными от перемещений (деформации) или связанными с первыми производными функциональными зави- > симостями (напряжения и т.п.). По времени: целые шаги, в которых аппроксимируются перемещения, силы и ускорения; полуцелые шаги, в которых ап-" проксимируются скорости. В итоге получим основную и две смещенные сетки (по пространству и по времени), при этом на основной сетке определены неизвестные и их производные четного порядка, а на смещенных - производные нечетного порядка. Данная конструкция является очень гибкой и удобной по следующим причинам: свойства материала задаются в ячейке совершенно независимо от основной сетки, заменяя "физический блок", легко получить материал с любой реологией; она ориентирована в общем случае на неортогональные и даже на неравномерные сетки, что позволяет применять ее в областях сложной формы; схема получается двухслойной по времени. В случае исключения всех неизвестных, кроме перемещений (это возможно, например, для линейно-упругого изотропного материала и приводит к системе уравнений Ламе) схемы, построенные по данному принципу, преобразуются к виду, аналогичному стандартным схемам для волнового уравнения: "крест", классической неявной схеме. Отметим, что при этом операторы аппроксимации вторых производных получаются автоматически как суперпозиции операторов первых производных.

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

Вариационно-разностные методы развиты в работах С.Г.Михлина, Л.А.Оганесяна, Г.И.Марчука, В.И. Агошкова и других. Применительно к задачам динамики тонкостенных конструкций вариационно-разностный метод развит в работах В.Г.Баженова и его учеников. Эти.методы отличаются от конечно-разностных тем, что в их основе лежит не дифференциальная, а вариационная задача. Близок к данному методу интегро-интерполяционный подход [193]. Под вариационно-разностными часто понимают проекционно-сеточные мето- , ды, т.е. метод конечного элемента. Однако в целом вариационно-разностные методы к ним не сводятся. Вариационно-разностным является любой метод, основанный на сеточной аппроксимации Ъариацйонного уравнения или вариационной задачи для некоторого функционала. При этом построение базисных функций необязательно. Таким образом, не любую вариационно-разностную схему можно считать схемой МКЭ и наоборот. Построение разрешающих соотношений схемы сводится к конечноразностной аппроксимации вариационного уравнения и приравнивания нулю коэффициентов при вариациях узловых перемещений. При этом интеграл представляется в виде суммы интегралов по ячейкам разносной сетки, которые в свою очередь выражаются приближенно через значения перемещений в узлах, принадлежащих ячейке.'Данный класс схем существенно уже конечноразностных. Среди них встречается много удачных схем, что свидетельствует о преимуществах вариационного, подхода; Отметим основные . - достоинства вариационно-разностных методов. 1 )Возможность использования неравномерных и нерегулярных сеток. При этом необходимо определить лишь инцидентность узлов и ячеек. 2)Единообразный расчет внутренних и граничных узлов. 3)Меныиие по сравнению с конечнораз-ностным методом требования к гладкости функций.

Указанные свойства делают вариационно-разностные методы очень удобными для программной реализации. Их можно применять для областей сложной формы. Полученные в результате разностные схемы по форме аналогичны разностной схеме Уилкинса и обладают сходными свойствами. Вместе с тем сам по себе вариационно-разностный метод не является панацеей от всех недостатков. Для высокого качества схем необходимо выполнение ряда дополнительных требований. Еще одним недостатком вариационно-разностных схем является возможность возникновения в них (как и в конечноразностных-. схемах) эффекта неустойчивости "песочные часы" [207} при использовании ячеек с числом узлов больше минимально возможного количества ( трех в двумерных задачах, четырех в трехмерных ).

Метод конечного элемента разрабатывался многими отечественными и зарубежными (в первую очередь - американскими) учеными. Большой вклад в его становление внесли Р.Курант, Дж. Аргирис, Р. Клаф, О.Зенкевич, Р.Галлагер, Дж. Оден, Ж.Деклу, Г.Стренг, Дж. Фикс и другие. Дальнейшее развитие метода связано с именами отечественных ученых: Л.А.Розина, А.С.Сахарова, В.А.Постнова, Р.Б.Рикардса, А.Н.Паутова, С.А.Капустина, А.И.Голованова и других. Метод получил широкое распространение при решении широкого круга задач механики деформирования твердого тела [71,72,77, 81,89-91,94,97,109-111,120-122,125-127,168,170,172,176,181-183,186,189,197, 201,219,224-249]. Он ведет свое происхождение от методов сил и перемещений в строительной механике, а" также от4 методов Рэлея-Ритца и Бубнова-Галеркина. Отличительной чертой МКЭ по сравнению с конечноразностными методами является то, что в нем неизвестные функции определены везде в области определения задачи, а не в дискретном наборе точек. При этом в качестве исходных уравнений для МКЭ могут использоваться как дифференциальные уравнения, так и вариационные принципы в зависимости от метода построения разрешающей системы уравнений (коллокации, взвешенных невязок, Рэлея-Ритца-Галеркина). Однако в большинстве случаев в задачах механики используется вариационный подход, что послужило основанием причислить МКЭ к вариационно-разностным методам. Благодаря своей физической наглядности, высокой алгоритмичности и удобству применения в областях сложной формы МКЭ получил широкое распространение, особенно в задачах статики. На его основе построено большое количество удачных численных схем. При этом многие из них имеют высокий порядок аппроксимации (эрмитовы, сирендипо-вы и т.п. элементы) и по своей сложности многократно превышают известные конечноразностные схемы. И только благодаря специфике подхода МКЭ они достаточно просто реализуются, в том числе и на неравномерных и нерегулярных сетках. Еще одним достоинством МКЭ является простота введения ограничений на вид неизвестных функций, наложения кинематических связей и т.п. (например, введение гипотез теории оболочек) непосредственно в дискретной модели независимо от исходных уравнений. Таким образом, МКЭ является удобным инструментом математического моделирования, а не просто одним из численных методов.

Применительно к задачам динамики конструкций у МКЭ в общем те же проблемы, что и у конечноразностных и вариационно-разностных методов. Это в первую очередь проблема малого параметра или вырождения задачи. Подробно эта проблема обсуждается в [97]. Ряд проблем возникает при решении динамических задач. Последовательный конечноэлементный подход должен относиться ко всей пространственно-временной области переменных. Реально же в абсолютном большинстве случаев используются неизвестные функции, непрерывные по пространству и дискретные по времени. Таким образом, по временной координате имеем дело не с конечноэлементной, а с конечноразно-стной аппроксимацией. Для построения явных конечноэлементных схем приходится при аппроксимации ускорений переходить к несогласованным диагональным матрицам масс, вносящим определенные противоречия в конечноэле-ментный подход. Еще одной проблемой в случае использования явных схем является их устойчивость. Если использовать явные схемы высокого порядка точности, ограничения на временной шаг и гладкость функций будут чересчур жесткими. Поэтому в динамических задачах в сочетании с явными схемами приходится ограничиваться простейшими -линейными и билинейными конечными элементами, особенно для нелинейных задач. Это же замечание справедливо и для неявных схем при использовании итерационных процедур решения полученных систем уравнений [153].

Сделаем еще два общих замечания, касающихся численных схем всех видов.-Первое из них касается порядка точности. Наибольшее распространение* в задачах динамики конструкций (и во многих других) получили схемы второго порядка точности. В газовой динамике распространены схемы первого порядка или гибридные схемы [95,222], обладающие свойством монотонности. Это в первую очередь связано с возникновением в этих задачах ударных врлн и необходимостью сглаживания численных решений. Поскольку возникновение ударных волн в тонкостенных конструкциях невозможно, то решение задач с гладкими начальными и граничными данными остаются гладкими и решения содержат только изначально существующие разрывы. Поэтому проблема монотонности отходит на второй план. Вырождение задач расчета тонкостенных конструкций (см. ниже) требует применения разностных схем не ниже второго порядка аппроксимации, поскольку даже многие схемы второго порядка являются условно аппроксимирующими [55,71,72]. Применение схем более высокого порядка необоснованно по следующим причинам. Повышение порядка схем выше второго по временной координате (а пространство и время в гиперболических системах равноправны) приводит к существенному усложнению схем и проблемам с их устойчивостью. Наиболее часто используемые модели оболочек Кирхгофа-Лява и Тимошенко также имеют второй порядок точности относительно толщины оболочки и при использовании шага сетки порядка толщины нет смысла применять более точную схему их реализации [164]. Отметим также, что класс задач с гладкими решениями довольно узок [153]. Таким образом, построение схем с высоким (выше второго) порядком, аппроксимации хотя и привлекательно, но не нашло до настоящего времени широкого практического применения.

Еще одной общей проблемой для численных методов всех вышеперечисленных классов зачастую является отсутствие возможности их теоретического сравнительного анализа с единой точки зрения. До сих пор основным методом сравнения численных схем остается численный эксперимент [127,135,197,241]. В первую очередь это касается схем, базирующихся на разных постановках задач - вариационной и дифференциальной. Так, существует разрыв между тео- -риями разностных схем и метода конечного элемента. Причины этого иногда объективны, поскольку теория разностных схем в отличие от МКЭ оперирует практически только с задачами на регулярных сетках. Но в случае применения регулярных и особенно равномерных конечноэлементных сеток разрыв между этими теориями может быть преодолен [109,216], и схемы МКР и МКЭ после преобразования могут быть исследованы едиными методами. В целом решение задачи эквивалентных преобразований численных схем позволило бы гораздо более осознанно подходить к их отбору, оценивая аппроксимацию различных форм записи исходных задач [195]. Но ее решение связано с установлением сеточных аналогов ряда формул и теорем математического анализа. На настоящий момент в этом направлении можно отметить только отдельные успехи: разностные аналоги формул интегрирования по частям и Грина [56,145,153,

181], формулы повторного дифференцирования [56,173,174], ряд других соотношений [173,174].

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

Проблема вырождения задачи. Данная проблема [1-3,5-7,71,72,185,200, 227,231,240,248,249] относится как к динамическим, так и к статическим задачам. Задачи теории оболочек относятся к вырожденным (теория Кирхгофа-Лява) или вырождающимся (теория типа Тимошенко и неклассические теории) [71,72]. Первые из них - формально параболического типа (фактически - уравнения типа колебаний), вторые - гиперболического с малым параметром. Если рассмотреть свободные изгибные колебания .пластин типа Тимошенко и Кирхгофа-Лява, то асимптотически при стремлении длины волны к бесконечности (или толщины оболочки к нулю при фиксированной длине волны) решения, полученные по обеим теориям, совпадают. Ту же асимптотику имеет и решение, полученное для тонкого упругого слоя на основе трехмерной теории упругости.

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

-■-"■ 18 ности двух близких величин и при этом имеет порядок h/L, где h -толщина пластины, L - длина волны. Поэтому нужно с особой тщательностью подходить к построению численных схем, чтобы в них погрешности порядка (ЛхЛг)2 ( где Ах - размер ячейки разностной сетки) не привели к фатальным последствиям. Ранние работы по численному решению задач теории оболочек с использованием модели Тимошенко [73,98,177] показывают пример неудачной ко-нечноразностной схемы, дающей приемлемую точность лишь при малом (меньше толщины) шаге пространственной сетки. Другим примером неудачной схемы для данной задачи является схема метода конечного элемента с линейной аппроксимацией и обладающая аналогичным недостатком. Примером удачной численной схемы является схема, построенная вариационно-разностным методом. Отметим, что все приведенные схемы имеют одинаковый второй порядок точности.

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

Этот эффект неоднократно описан под разными названиями в различных работах, особенно по методу конечного элемента. Для устранения указанного недостатка использовались: сокращенное интегрирование [248,249], способ двойной аппроксимации [71,72], другие подходы [1-3,120,121,185,231,240]. Данная проблема в МКЭ подробно обсуждается в [97]. Близко к ней лежит проблема линейного треугольного оболочечного элемента, неудачные опыты с которым по сути привели к отказу от его использования. Подробный анализ причин недостатков схем на треугольных ячейках приведен в [56]. В целом, резюмируя вышесказанное по данной проблеме, отметим следующее. Задачи динамики тонкостенных конструкций являются вырожденными или вырождающимися, т.е. содержащими малый параметр - толщину. В последнем случае низшие и важнейшие частоты длинноволновых процессов получаются как разность близких больших величин. Поэтому малые относительные погрешности этих больших величин могут приводить к большим относительным погрешностям их разности. В связи с этим при построении численных схем следует либо добиваться высокого порядка точности (выше второго по пространственным координатам) и при этом влияние этих малых погрешностей будет действительно мало, либо для схем второго порядка точности очень тщательно следить за согласованностью погрешностей разных членов уравнений. Все рассмотренные выше приемы посвящены последнему подходу. В итоге же неудачные схемы отличаются от удачных отсутствием равномерной сходимости по параметру сеточной задачи AX/h, где АХ - шаг сетки по пространственным координатам, h - толщина оболочки. Реально эти недостатки проявляются при использовании сдвиговой модели оболочки для AX/h >0.3. При использовании трехмерной теории неудачные схемы будут работать плохо на сетках с параметром AXmw/AXmin>5. Проведем сравнение удачных и неудачных схем с точки зрения эффективности. Увеличение шага пространственной сетки в 2 раза позволяет на двумерной задаче уменьшить число пространственных узлов в 4 раза, во столько же раз уменьшается число операций на одном временном шаге. Если реально для широкого класса задач удается без заметной потери точности и трудоемкости вычислений увеличить шаг в 5-10 раз, то эффект очевиден.

Динамические проблемы. Рассмотрим трудности , возникающие исключительно при решении динамических задач. Во-первых, исследователь должен сделать выбор между явными и неявными схемами. Общеизвестно, что у явных схем временной шаг определяется устойчивостью схемы, а у неявных абсолютно устойчивых - точностью. Заметим, что при решении систем итерационными методами в неявных схемах могут возникать ограничения на временной шаг, близкие к условию Куранта-Фридрихса-Леви (К-Ф-Л) [153]. При этом нужно - учитывать, что трудоемкость неявных схем на шаге по времени больше.

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

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

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

Обычно спектр задачи связан со временем процесса. Коротковолновые процессы - это процессы ударные, импульсные, содержащие крутые фронты. Как правило, они и кратковременные, если учитывать диссипацию энергии. В длительных процессах коротковолновая составляющая обычно несущественна. Отсюда можно сделать вывод, что явные схемы предпочтительнее неявных для кратковременных процессов, и конкурентноспособны с ними в некотором промежуточном диапазоне. Для длительных (более периода по основной форме) колебательных процессов с гладкими решениями следует отдать предпочтение неявным схемам. Отметим связь геометрии конструкций с эффективностью явных или неявных схем. Наличие у конструкций различных мелких геометрических особенностей (отверстий, выточек, жестких включений и т. п.) требует более точного описания напряженно-деформированного состояния в их окрестности, а соответственно и более мелкой сетки. В этом случае эффективность явных схем падает в связи с условием К-Ф-Л. Задачи со сложной геометрией являются жесткими [185], поскольку необходимо учитывать разномасштабные по пространству (соответственно и по времени) процессы. В целом, чем проще геометрия^ тем выше эффективность явных схем по сравнению с неявными.

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

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

Способы повышения эффективности численных методов, использующих явные схемы. При решении большинства задач для тонких оболочек достаточная точность достигается при использовании пространственного шага сетки вдоль срединной поверхности порядка толщины оболочки. Если применяется трехмерная теория упругости, выгодно использовать ячейки, сплющенные по толщине. В том и другом случаях условия устойчивости явных схем становятся чересчур обременительными. Этот факт был отмечен давно, и рядом авторов предпринимались попытки тем или иным способом увеличить временной шаг, в целом оставаясь в рамках явного подхода [16,69,85,107,112,221,229-230,237

239]. Все они так-или иначе основаны на противоречии между спецификой волновых.процессов в оболочках и явными схемами. Суть его в том, что первыми численную устойчивость теряют формы колебаний, направленные перпендикулярно срединной поверхности оболочки, т.е. те, для которых волновая динамика фактически отсутствует. Предлагались следующие подходы: с выделением высокочастотных форм - работа [107], в которой рассмотрена сеточно-характеристическая схема с выделением колебаний по толщине оболочки; с подавлением высокочастотных колебаний - работы [51,54].

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

Проблемы геометрической нелинейности [90,120,122,150,155,157,198, 226,233,247]. Эти проблемы связаны с двумя факторами - наличием больших перемещений и углов поворота и использованием криволинейных систем координат. В совокупности они порождают проблему смещения тела как жесткого целого, являющуюся индикатором правильности геометрического описания процесса деформирования. Данные проблемы поднимались в работах [97,198].

Проблема неустойчивости "песочные часы" иди мод нулевой энергии. Данная проблема возникает в конечно-разностных и вариационно-разностных схемах [23,25,33,207] и связана с неполнотой систем разностных операторов, аппроксимирующих функцию и производные в ячейках разностной сетки. Она возникает, когда число узлбв в ячейке превышает число используемых разностных операторов для аппроксимации производных. Типичные случаи - схема на четырехугольных ячейках в плоском случае (три оператора, четыре узла, одна мода нулевой энергии) и шестигранная ячейка в трехмерных задачах (четыре оператора, восемь узлов, четыре моды нулевой энергии). Неполнота систем разностных операторов в ячейках приводит к тому, что сеточный оператор задачи утрачивает свойство положительности. При этом в ячейке возникает как-бы свободно перемещающийся кинематический механизм, закрепление которого происходит только в граничных узлах с заданными условиями в перемещениях. Поэтому в некоторых задачах неустойчивость ''песочные часы" не возникает (это задачи с граничными условиями в перемещениях и сравнительно небольшим числом узлов сетки). Для подавления этой неустойчивости было предложено использовать искусственную вязкость специального типа, повысить порядок аппроксимации по пространственным переменным [242], вводить в расчеты дополнительные треугольные (тетраэдральные) ячейки [32], вычислять моментные составляющих деформаций с помощью аппроксимации линейными функциями[23,25]. В последнем варианте можно получить дискретные соотношения, аналогичные принятым в шестимодальном варианте теории оболочек типа Тимошенко, и при отсутствии локальных импульсных воздействий решать нелинейные задачи динамики тонкостенных конструкций на сетке с одним слоем конечных элементов по толщине. Такой подход позволяет в рамках единого алгоритма исследовать нестационарное деформирование составных конструкций, включающих массивные тела и оболочки. Другие подходы к устранению эффекта «песочных часов» предлагались в работах [99,192,234].

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

Ряд" специфических проблем возникает при решении актуальных в настоящее время задач расчета многослойных конструкций и конструкций из композитных материалов [6,64,67,82].

Большое внимание в литературе уделяется решению задач контактного взаимодействия деформируемых тел и конструкций с твердыми телами, жидкими и грунтовыми средами [24,70,83,100,103,104,113,138,139].

Большая библиография по вопросам, связанным с рассмотренной темой имеется в опубликованных» книгах и обзорах последних (и более ранних) лет [17,22,64,83,101,123,204]. Таким.образом, в целом вопросы математического обоснования и теоретического анализа численных методов решения задач механики конструкций являются недостаточно изученными. Поэтому актуальны вопросы разработки и развития теоретических методов анализа точности, устойчивости и эффективности численных схем решения задач динамики тонкостенных конструкций, а . также применение этих методов к исследованию существующих и разработке новых методов решения задач указанного класса.

Цели работы:

- развитие методов исследования вариационно-разностных и конечно-элементных схем решения задач механики деформируемого твердого тела;

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

- анализ и математическое обоснование точности, устойчивости и методов повышения эффективности численных схем;

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

Научная новизна

1) Разработан метод анализа вариационно-разностных и конечно элементных схем, основанный на преобразовании к конечно-разностному виду, удобному для теоретического исследования. 2) На базе данного подхода проведен анализ точности и сходимости численных, схем решения задач теории упругости и теории оболочек. Впервые получены и теоретически обоснованы условия применения грубых сеток при расчете тонкостенных конструкций. Исследована равномерная сходимость схем по параметру сеточной задачи AX/h (где АХ- диаметр ячейки разностной сетки, h - толщина оболочки). Предложен класс «ажурных» численных схем метода конечного элемента. Проведен анализ устойчивости явных схем типа «крест» решения задач теории упругости и 'теорий оболочек;- Получены точные и* приближенные- оценки устойчивости.

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

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

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

Практическая ценность работы Результаты проведенных в диссертационной работе исследований находят практическое применение в 111111 «Дина-мика-2» и «Динамика-3» решения нестационарных, задач динамики конструкций-и массивных тел, разработанных в НИИ механики при ННГУ. Созданные на их основе методики и рекомендации по их применению внедрены в расчетную практику рядзотраслевых и проектно-конструкторских организаций.

На защиту выносятся:

- метод анализа вариационно-разностных и конечно-элементных схем, основанный на построении эквивалентных конечно-разностных схем;

- результаты исследования точности вариационно-разностных схем теории пластин и оболочек: анализ равномерной сходимости по параметру сеточной задачи AX/h, условия применения грубых сеток в задачах теории оболочек;

- новый класс «ажурных» численных схем метода конечного элемента;

- оценки устойчивости явных схем «крест» решения задач теории упругости и теории оболочек; ■ .———-------

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

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

Апробация работы Результаты диссертационной работы докладывались на следующих конференциях, симпозиумах и семинарах: Республиканская научно-техническая конференция «Механика сплошных сред» (Набережные челны, 1982), Всесоюзная конференция «Численная реализация физико-механических задач прочности» (Горький, 1983), семинар «Прикладные методы расчета физических полей» (Крым, пос. Кацивели, 1984), V Всесоюзная конференция по статике и динамике пространственных конструкций (Киев,

1985), VIII Всесоюзный симпозиум по распространению упругих и упругопла-стических волн (Новосибирск, 1986), II Всесоюзный симпозиум «Устойчивость в механике деформируемого твердого тела» (Калинин, 1986), VI Всесоюзная школа «Теоретические основы и конструирование алгоритмов решения задач математической физики» (Горький, 1986), 1,11,III,IV Всесоюзные школы молодых ученых «Численные методы механики сплошной среды» (Шушенское, 1987, Абакан, 1989, Абрау-Дюрсо, 1991, 1992), XVI,XIX Международные конференции по теории оболочек и пластин (Нижний Новгород, 1993, 1999), XIV Международная конференция «Численные методы решения задач упругости и пластичности» (Волгоград, 1995), X Юбилейная Международная конференция по вычислительной механике и современным программным системам (Переславль-Залесский, 1999), VI Международный симпозиум «Динамические и технологические проблемы механики конструкций и сплошных сред» (Ярополец, 2000), Международная научно-техническая конференция «Испытания материалов и конструкций», посвященная 100-летию со дня рождения проф. М.И.Волского (Нижний Новгород, 2000), .XX Международная конференция «Математическое моделирование в механике сплошных сред. Метод граничных и конечных элементов» (Санкт-Петербург, 2003), на семинаре НИИ механики при Нижегородском государственном университете им. Н.И.Лобачевского.

Публикации По теме диссертации опубликована 41 работа, в том числе 1 монография и 1 учебное пособие. Основные результаты диссертации отражены в работах [4,35,44,45,51 -61,180,212-217].

Структура и объем работы Диссертация состоит из введения, 7 глав, заключения и списка литературы из 249 наименований. Работа содержит 263 страницы основного текста, включая 77 рисунков и 4 таблиц.

Основные результаты и защищаемые положения диссертации опубликованы в работах [4,35,44,45,51-61,180,21 1

239

Заключение

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

1. Разработан метод анализа вариационно-разностных и КЭ численных схем, состоящий в построении эквивалентного им конечно-разностного представления. Он позволяет проводить теоретический анализ методов численного решения широкого класса линейных и нелинейных задач. Его применение позволило провести исследования ряда, свойств дискретных моделей, недоступных при стандартном теоретическом анализе МКЭ.

2. Проведены исследования аппроксимации, устойчивости и сходимости численных схем теории пластин и оболочек Тимошенко. Определены условия равномерной сходимости численных схем по параметру сеточной задачи AXfh. Обосновано применение грубых по сравнению с толщиной оболочки сеток. Предложен новый класс «ажурных» схем МКЭ. 3. Получены оценки устойчивости явных численных схем решения задач теории упругости и теории оболочек. Установлен эффект граничной неустойчивости трехмерных схем на свободных границах, содержащих ребра и углы и проведено исследование частного случая данной неустойчивости.

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

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

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

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

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

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

1. Абдукадыров С.А., Пинчукова Н.И., Степаиенко М.В. Об одном способе численного решения уравнений динамики упругих сред и конструкций//Физико-технические проблемы разработки полезных ископаемых. 1984. №6. С.34-41.

2. Абовский Н.П., Андреев Н.П., Деруга А.П. Вариационные принципы теории упругости и теории оболочек. М.: Наука, 1978. 287 с.

3. Абросимов Н.А., Баженов В.П, Кибец А.И., Садырин А.И., Чекмарев Д.Т. Нелинейные задачи динамики конструкций// Математическое моделирование; 200Q.rT.12.N.6." С:47-50.~- -—

4. Абросимов Н.А., Баженов В.Г. Исследование упруго-пластических процессов деформирования пластин и оболочек вращения при импульсном нагружении в неклассической постановке //Прикх механика. 1985. Т.21. № 1. С. 73-79.

5. Абросимов Н.А., Баженов В.Г.Нелинейные задачи динамики композитных конструкций. Н.Новгород: Издательство ННГУ, 2002. 400 с.

6. Абросимов Н.А., Баженов В.Г. Об одном методе решения нелинейных задач динамики оболочек в уточненной постановке //Прикладные проблемы прочности и пластичности: Всесоюз. межвузю сб. // Горьк. ун-т. Горький, 1975. Вып. 1. С. 58-66.

7. Агамиров В.Л., Вольмир А.С. Об устойчивости цилиндрической оболочки при продольном ударе // Докл. АН СССР, 1964. Т. 157, №. 2. С. 307-308.

8. Адищев В.В., Корнев В.М., Талзи Л.А. Оценки максимальных напряжений в замкнутых цилиндрических сосудах при осесимметричном взрывном нагружении // Новосибирск: Ип-т гидродинамики СО АН СССР. 1983. 75 с. -Деп. в ВИНИТИ 05.12.1983. № 6588-83/

9. Ю.Айнола Л .Я. Вариационные принципы динамики теории оболочек.- ДАН СССР, 1967, т. 172, № 6, с. 1296-1298.

10. Айнола Л.Я. Нелинейная теория типа Тимошенко для упругих оболочек // Изв. АН ЭССР. Сер. физ.-мат. и техн. наук. 1965. Т.14. № 3.C.337-344,

11. Айнола Л.Я. О расчетных моделях упругих пластинок для динамических задач. Изв. АН Эст. ССР, сер. физ.-матем. и техн. наук, 1963, т.12 № 1. С. 31-37.

12. Айнола Л.Я. Уравнения теории типа Тимошенко упругих оболочек в усилиях и моментах,- В сб.: Переходные процессы деформаций оболочек и пластин/ Таллин: 1967.

13. Айнола Л.Я., Нигул У.К. Волновые процессы деформации упругих плит и оболочек // Изв. АН ЭССР: Сер. физ.-мат. и техн. наук. 1965. Т.14. № 1. С.З-63.

14. Алфутов Н.А. Основы расчета на устойчивость упругих систем. М.: Машиностроение, 1978. 311с.

15. Анисимов С.А. Векторное расщепление плоской динамической задачи теории упругости в областях из произвольных четырехугольников // Динамика сплошной среды. Новосибирск: Ин-т гидродинамики СО АН СССР, 1986. Вып. 75. С. 17-26

16. Аннин Б.Д. Механика деформируемого твердого тела в СО РАН в 1988-1997 годы // ПМТФ. 1997. Т. 38. № 4 С 28-45.

17. Афанасьев С.Б., Козлов Е.А. Алгоритм решения двумерных волновых упругопластических задач методом Годунова // Прикладные проблемы прочности и пластичности. Алгоритмизация и автоматизация: Всесоюзн. межвуз. сб. /Горьк. ун-т. Горький, 1987. С.91-100.

18. Бабенко К. И. Основы численного анализа. М.: Наука, 1986. 744 с.

19. Баженов В.Г. и др. Пакет прикладных программ "Динамика-2'7/ Прикладные проблемы прочности и пластичности. Алгоритмизация и автоматизация исследований: Всесоюз. межвуз. сб. / Горьк. ун-т, Горький, 1987. С. 4-13.

20. Баженов В.Г. Численное исследование нестационарных процессов деформации упругопластических оболочек // Проблемы прочности. 1984. №11. С. 51-54.

21. Баженов В.Г. Численное моделирование нестационарных задач динамики упругопластических конструкций // Прикладные проблемы прочности и пластичности: Всесоюз. межвуз. сб. / М. 1995. Вып. 53. С. 17-29.

22. Баженов В.Г., Игоничева Е.В., Кибец А.И., Лаптев П.В., Ломунов В.К. Выпучивание упругих и упругопластических оболочек вращения при осевомударном пагружсшш//Изв.АИП РФ. Волго-Вятское региональное отделение/ НГТУ. Москва-11.Новгород.2001 .С.7-23.

23. Баженов В.Г., Игоничева Е.В.Нелинейные процессы ударного выпучивания упругих элементов конструкций в виде ортотропных оболочек вращения. -Нижний Новгород :ННГУ, 1991.

24. Баженов В.Г., Кибец А.И. КЭ решение трехмерных нестационарных задач упругопластического динамического деформирования элементов конструкций//Сб. научных трудов "Устойчивость, пластичность, ползучесть при сложном нагружении". Вып.2. Тверь,2000. С.55-60.

25. Баженов В.Г., Кибец А.И. Численное моделирование трехмерных задач нестационарного деформирования упругопластических конструкций методом конечных элементов//Изв. РАН. МТТ. 1994. № 1. С. 52-57.

26. Баженов В.Г., Кибец А.И., Цвсткова И.Н. Числсиное моделирование нестационарных процессов ударного взаимодействия деформируемыхэлементов конструкций/ЛПроблемы машиностроения и надежности машин. 1995. №2. С. 20-26.

27. Баженов В.Г., Ломунов В.К. Методика расчета динамического деформирования геометрически изменяемых плоских стержневых систем // Проблемы прочности и пластичности: Межвуз.сб., 2002, Вып. 64, с. 55-63.

28. Баженов В.Г., Ломунов В.К. Экспериментально-теоретическое исследование упругопластического выпучивания цилиндрических оболочек при осевом ударе И Прикладная механика. 1983.Т.19,N 6. С.63-69.

29. Баженов В.Г., Ломунов В.К., Петров М.В. Упругопластическое деформирование цилиндрических оболочек прн магнитно-импульсном нагружении// Прикладные проблемы прочности и пластичности : Всесоюз. Межвуз. сб./Горьк.ун-т. 1979. С. 73-78.

30. Баженов В.Г., Ломунов В.К., Петров М.В., Угодчиков А.Г. Исследование; больших вязкопластических деформаций цилиндрических оболочек с применением магнитно-импульсного способа нагружения. / Машиноведение^ №5, 1983. С.73-80.

31. Баженов В.Г., Ломунов В.К. Исследование упругопластического выпучивания оболочек при ударном нагружении //Прикладные проблемы прочности и пластичности: Всесоюз.межвуз.сб. / Горьк. ун-т. Горький, 1975. Вып.2.С.44-50.

32. Баженов В.Г., Ломунов В.К. Устойчивость и закритическое состояние оболочек вращения при осевом ударе//Прикладная механика. 1986.Т.22, N 9.С.28-33.

33. Баженов В.Г., Мухина А.С., Угодчиков А.Г.Упругопластическое выпучивание цилиндрических оболочек с начальными несовершенствами под действием импульса внешнего давления //Избранные проблемы прикладной механики.- М.:ВИНИТИ,1974.С.73-81.

34. Баженов В.Г., Пирогов С.А., Чекмарев Д.Т. Явная схема со стабилизирующим оператором для решения нестационарных задач динамики конструкций// Изв.РАН.МТТ.2002.№ 5. С. 120-130.

35. Баженов В.Г., Рузанов А.И.,Угодчиков А.Г. О численных методах и результатах решения нестационарных задач теории упругости и пластичности//Численные методы механики сплошной среды. 1985. Т. 16, Ж 4. С. 129-149.

36. Баженов В.Г., Угодчиков А.Г., Шинкаренко А.П.Численный анализ упругопластического деформирования оболочек с криволинейными отверстиями при импульсном нагружении //Прикладная механика. 1979.T.15,N 5.С.48-53.

37. Баженов В.Г., Чекмарев Д.Т. Вариационно-разностные схемы в нестационарных волновых задачах динамики пластин иоболочек //Н.Новгород: Изд-во Нижегород. ун-та, 1992. 159 с.

38. Баженов В.Г., Чекмарев Д.Т. Выпучивание упругопластических сферических куполов под действием импульса давления/ЯТрикладные проблемы прочности и пластичности. Статика и динамика деформируемых:систем Всешюз«Л1ежву^сб^ Л7орьк!.у1^Горьки%

39. Баженов В.Г.,.Чекмарев Д.Т. Вычислительные модели нелинейных задач динамики пространственных конструкций // Актуальные проблемы механики оболочек^ Труды международной конференции. Новое знание* Казань/2000:

40. Баженов В.Г., Чекмарев Д.Т. Об индексной коммутативности численного дифференцирования //Ж. вычисл. математики и мат. физики. 1989. Т. 29. № 5. С. 662-674. .

41. Баженов В.П, Чекмарев Д.Т. Численные методы решения задач динамики тонкостенных конструкций// Изв.РАН.МТТ.2001.№ 5. С.156-173;

42. Баженов В.Г., Чекмарев Д.Т. Решение задач нестационарной динамики пластин и оболочек вариационно-разностным методом .Учебное пособие. Нижний Новгород, изд-во ННГУ, 2000. 118 с.

43. Баженов В.Г., Шинкаренко А.П. Вариационно-разностный метод решения двумерных задач динамики упругопластических оболочек// Прикладные проблемы прочности и пластичности: Всесоюз.межвуз. сб./Горьк.ун-т.Горький, 1976.Вып.З .С.61-69.

44. Баженов В.Г., Шинкаренко А.П. Упругопластическое деформирование составных оболочечных конструкций при импульсных воздействиях // Проблемы прочности. 1981. № 3. С. 25-29. • . : . . . . . .

45. Бакулин В.Н., Образцов И.Ф., Потопахин В.А. Динамические задачи нелинейной теории оболочек. Действие интенсивных термосиловых нагрузок, концентрированных потоков энергии. М.: Наука. Физматлит,1998. 463 с.

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

47. Бермус И.М., Срубщик Л.С. Осесимметричное динамическое прощелкивание упруго-пластических сферических оболочек // Вопросы оптимального проектирования пластин и оболочек /Сарат.ун-т. Саратов, 1981.С. 12-15.

48. Богданович А.Е. Нелинейные задачи динамики цилиндрических композитных оболочек. Рига: Зинатне, 1987. 295 с.

49. Богданович А.Е. Обзор исследований по устойчивости цилиндрических оболочек при осевом динамическом сжатии // Электродинамика и механика ' сплошных сред/Латв.ун-т.Рига, 1980.

50. Богульский И.О. Об одной схеме расщепления решения двумерной задачидинамики твердого тела. Красноярск, .1989. 38 С. Деп в ВИНИТИ 21.03.89. № / 1816-В89.

51. Болдычев В.П. Повышение эффективности метода конечных элементов при решении вырождающихся задач // Вопросы динамики и прочности. Рига: Зинатне, 1983. Вып. 42. С.38-48.

52. Борисенко В.И.,Клокова А.И. Закритическая деформация цилиндрической оболочки при ударе //Прикладная механика. 1966.Т. 11, N 10. С.29-35.

53. Бреббия К., Теплее Ж., Вроубел Л. Методы граничных элементов / Пер. с англ. Л.Г.Корнейчука. М.: Мир, 1987. 524 с.

54. Бригадиров Г.В.Долоконников Л.А.Удар цилиндрической оболочки о жесткую преграду//Изв.АН СССР. MTT.1983,N 3. С.180-181.

55. Бураго Н.Г. Кукуджанов В.Н. Распространение волн в материалах с запаздыванием текучести. // Распространение упругих и упругопластических волн: Матер. V Всесоюз. симп. Алма-Ата, 1973.

56. Бураго Н.Г., Кукуджанов В.Н. Решение упругопластических задач методом конечных элементов. Пакет прикладных задач " Астра". М., 1988. 63 с. (Препринт / Ин-т проблем механики АН СССР; № 326.

57. Бураго Н.Г.,Кукуджанов В.Н.Выпучивание и закритическце деформации упругопластических оболочек в условиях осевой симметрии//Сб. по численным методам в механике деформируемого твердого тела / В1Д АН СССР. М.1978.С.47-66.

58. Валишвили Н.В.Методы расчета оболочек вращения на ЭЦВМ.-М-: Машиностроение, 1976. ,

59. Варвак П.М.,Варвак М.Ш.Приложение двух интегральных представлении -производной в теории решеток//Прикладная механика. 1971.T.7.,N 8. С.77-82.

60. Васидзу К. Вариационные методы в теории упругости и^пластичнек^иг М.г; Мир, 1987. 542 с.

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

62. Вестяк А.В., Горшков А.Г. Тарлаковский Д.В. Нестационарное взаимодействие деформируемых тел с окружающей средой // Итоги науки и техники. Механика деформируемого твердого тела. М.: ВИНИТИ, 1983. Т. 15. С. 69-148.

63. Власов В.З.Общая теория оболочек и ее приложения в технике.М.-Л.:ГИТ-ТЛ,1949.785 с.

64. Вольмир А.С. Нелинейная динамика пластинок и оболочек. М.:Наука. 1972. 432 с.

65. Вольмир А.С.Устойчивость деформируемых систем.-М.Наука, 1967.

66. Ворович И.И.Минакова Н.И. Проблема устойчивости и численные методы в теории сферических оболочек//Итоги науки и техники. Механика деформируемых твердых тел.- М.: ВИНИТИ, 1973 .Т.7.С.2-86.

67. Ву Р.В.Г., Уитмер Е.А.Исследование нелинейных неустановившихся реакций конструкций методом пространственных конечных элементов //Ракетная техника и космонавтика.1973.К8.С.67-70.

68. Ву Р.В.Г., Уитмер Е.А.Расчет неустановившихся больших упруго-пластических деформаций простых конструкций методом конечных:

69. By Р.В.Г., Уитмер Е.А.Устойчивость метода Де Вожела численного интегрирования во времени // Ракетная техника икосмонавтика. 1973.N.10.C.97-100. . .■.:: ^ . . ; ■

70. Галимов К.З. Основы нелинейной теории тонких оболочек. Казань.: Йзд-во Казан.ун-та, 1975.

71. Галимов К.З., Паймушин В.Н. Теория оболочек сложной геометрии. (Геометрические вопросы теории оболочек). Казань: Изд-во Казанск. ун-та, 1985 164 с.

72. Галлагер Р. Метод конечных элементов. Основы. М.: Мир. 1984.

73. Году нов С.К., Забродин А.В. Иванов М.Я. и др. Численное решение многомерных задач газовой динамики. М.: Наука, 1976. 400 с.

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

75. Голованов А.И., Корнишин М.С. Введение в метод конечных элементов статики тонких оболочек. Казань: Физ.-техн. ин-т 1990.269 с.

76. Гордиенко Б.А. О машинном решении задач ударного выпучивания упругих систем методом конечных разностей // Изв. АН СССР. МТТ. 1970. № 3. С. 143-148.

77. Горельский В .А., Хорев И.Е., Югов Н.Т. Динамика трехмерного процесса несимметричного взаимодействия деформируемых тел с жесткой стенкой // ПМТФ. 1985. № 4. С.112-118.

78. Гулидов А.И., Шабалин И.И. Расчет контактных границ с учетом трения при динамическом воздействии деформируемых тел в пространственном случае. // Числ. методы решения задач теории упругости и пластичности: Матер. X Всесоюз. конф., Новосибирск, 1988.

79. Юб.Дресвянников В.И. О численной реализации нелинейных уравнений динамики упруго-пластических оболочек // Прикладные проблемы прочности и пластичности: Всесоюз. межвуз. сб. / Горьк. ун-т. Горький, 1976. Вып: 31 С. 82-90.

80. Ю7.Евсеев Е.Г., Семенов А.Ю. Метод для численного решения уравнений динамики тонкостенных оболочек, основанный на выделении сильноосциллирующих компонент // Докл. АН СССР. 1990. Т.310. № 4. С. 785-788.

81. Ю.Зенкевич O.K. Метод конечных элементов в технике. М.: Мир, 1975. 541 с. .

82. HI.Зенкевич O.K. Метод конечных элементов: от интуиции к общности// Механика. Сб. переводов. 1970 № 6 С. 90-103.

83. Иванов Г.В. Построение схем решения плоском динамической задачи теорий'упругости на основе аппроксимации линейными- ноли«омами. /А^ — Динамика сплошной среды. Новосибирск: Ин-т гидродинамики СО АН СССР, 1978. Вып. 37. С. 63-77. ~

84. И 3 .Иващенко К.Б^ Алгоритм расчета контактных границ при взаимодействии деформируемых твердых тел.//Проблемы прочности. 1989, №2. С.79-82.

85. Ильюшин А.А. Механика сплошной среды. М^ Изд-во МГУ, 1990. 310с.

86. I 5.Ильюшин А.А. Пластичность. М.: Гостехиздат, 1948.

87. Пб.Ильюшин А.А. Пластичность. Основы общей математической теории. М.: Изд-во АН СССР, 1963.

88. Ишлинский А.Ю. Общая теория пластичности с линейным упрочнением // Украинский математический журнал. 1954. № 6. С. 314-325.

89. Кадашевич Ю.И., Новожилов В.В. Теория пластичности, учитывающая остаточные микронапряжения // Прикладная математика и механика. 1958. Т. 22, вып. I. С. 78-89.

90. Кадашевич Ю.И., Новожилов В.В. Теория пластичности, учитывающая- эффект Баушингера //. ДАН CCCPU957.T.U7, ran Jl» С. 586-588

91. Капустин С.А. Метод конечных элементов в задачах механики деформируемых тел.Учебное пособие.Н.Новгород:Изд-во ННГУ, 2002. 180 с

92. Кармишин А.В., Жуков А.И. и др. Методы динамических расчетов и испытаний тонкостенных конструкций. М.: Машиностроение, 1990.

93. Качанов Л.М. Основы теории пластичности. М.: Наука, 1969.420 с.

94. Кийко И.А.Цилиндрическая оболочка под действием.осевой ударной нагрузки //Изв.АН СССР.МТТ. 1969.N 2.С. 135-138.

95. Кильчевский Н.А. Теория нестационарных динамических процессов в оболочках. // Прикладная механика. 1968. Т. 4. Вып. 3. С. 1-18.

96. Киселев А.Б. Развитие метода Уилкинса для решения трехмерных задач соударения деформируемых тел. // Взаимодействие волн в деформируемых средах. М.: МГУ, 1984. С. 93-102.

97. Г.Ковалев A.M. Линейная осесимметричная реакция составной оболочки вращения на ударную нагрузку. Изв. АН СССР. МТТ, 1981, № 1, С. 177-184.

98. Ковеня В.М., Яненко Н.Н. Метод расщепления в задачах газовой динамики. Новосибирск: Наука. Сибирское отд-ние, 1981. 304 с.

99. Коларов Д.,Балтов А.,Бончева Н.Механика пластических сред / пер. с болгарского под ред. Г.С.Шапиро.-М.: Мир, 1979.

100. Кондауров В.И., Кукуджанов В.Н. Соударение жесткого цилиндра со слоистой упругопластической преградой. // Числ. методы решения задач теории упругости и пластичности: Материалы VI Всесоюз. конф. 4.1. -Новосибирск: ИТПМ СО АН СССР, 1980. С.84-91.

101. Корнеев В.Г. Сопоставление метода конечных элементов с вариационно-разностным методом решения задач теории упругости // Изв. Всесоюз. НИИ Гидротехники. 1967. Т. 83. С. 286-307.

102. Зб.Корнишин М.С.Нелинейные задачи теории пластин и пологих оболочек и методы их решения.-М.Наука, 1964.

103. Коротких Ю.Г. Численный метод исследования поведения тел приимпульсных воздействиях // Методы решения задач упругости и . пластичности: Учен. зап. / Горьк. ун-т. 1970. Вып. 122 (3). Сер. механика. С. 51-68.

104. Кошур В.Д., Немировский Ю.В. Континуальные и дискретные модели динамического деформирования элементов конструкций. Новосибирск: Наука, 1990. 199 с.

105. Кукуджанов В.Н. Численные методы решения неодномерных задач динамики упругопластических сред // Численные методы решения задач теории упругости и пластичности: Материалы YI Всесоюз. конф.^Новосибирск, Ш- Г20Г—^

106. Кукуджанов В.Н., Кондауров В.И; Численное решение неоднрмерных задач, динамики твердого деформируемого тела. // Проблемы динамики упругошь

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

108. Ладыженская О.А. Краевые задачи математической физики. М.: Наука, 1973.407 с.

109. Лахе А.Я., Нигул У.К. Анализ нелинейных осесимметричных волновых процессов деформации конических и цилиндрических оболочек методом характеристик// Труды 8 Всесоюз. конф. по теории оболочек и пластин, Ростов-на-Дону, 1971. /М.: Наука, 1973. С. 509-514.

110. Ландау Л.Д., Лившиц Е.М. Теоретическая физика. Т. VII. Теория упругости.- М.: Наука, 1987.248 с.

111. Луговой П.З. Динамика оболочечных конструкций при импульсных нагрузках (обзор) // Прикладная механика. 1990. Т. 26. № 8. С. 3-20.

112. Луговой П.З., Мукоид В.П., Мейш В.Ф. Динамика оболочечныхконструкций при взрывных нагрузках. Киев: Наук, думка, 1991. 278 с.

113. Мазалов В.Н., Немировский Ю.В.Динамика тонкостенных пластических конструкций//Проблемы динамики упруго-пластических сред.-М.:Мир,1975.Вып.5. С. 155-247.

114. Малышев А.П. Исследование переходных процессов в оболочечных конструкциях на основе схемы с минимальной дисперсией // Изв. АН СССР. МТТ. 1981. № 3. С. 66-73:

115. Малышев А.П., Паничкин В.И. Нелинейные волновые процессы в оболочках вращения // Изв АН СССР. МТТ, 1976, № 4, С. 175-178.

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

117. Марчук Г.И. Методы расщепления. М.: Наука, 1988. 263 с.

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

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

120. Морино Л., Лич Дж.В., Уитмер Е.А. Уточненный метод численного расчета нестационарных процессов в упруго-пластических тонких оболочках при больших деформациях//Груды Амер. Об-ва инж. мех:, серия Е, 1971, N 2, ч. 2 С. 157-164.

121. Муштари Х.М., Галимов К.З. Нелинейная теория упругих оболочек./Казань.Таткнигоиздат, 1957.

122. Навал И.К., Пацюк В.И., Римский В.К. Нестационарные волны в деформируемых средах. Кишинев: Штиинца, 1986. 236 с.

123. Нейман Дж., Рихтмайер Р. Метод численного расчета гидродинамических скачков //Механика: Сб. пер. 1951. № Г. V16Г.Нечипорук Г.С. Тен Ен Со. Экспериментальное исследование ударного выпучивания цилиндрических и конических оболочек //Изв. АН

124. CCCPJVITT. 1974.N З.С. 175-J 82^ ;. . . : , ;.Г.

125. Ни гул У.К Сопоставление результатов анализа переходных волновых процессов в оболочках и пластинах по теории упругости и приближенным теориям.//ПММ. 1969, Т. 33. Вып. 2. С. 308-322.

126. Нигул У.К. Асимптотическая теория статики и динамики упругих круговых цилиндрических оболочек. ПММ, 1962, т. 26, вып. 5, с. 923-930;

127. Новожилов В.В. Математические модели и точность инженерных расчетов. Л.: Судостроение. 1979. № 71. С. 5-12.

128. Новожилов В.В. Основы нелинейной теории упругости. М.-Л., Гостехиздат, 1948.

129. Новожилов В.В. Теория упругости Л.; Судпромгиз, 1958.*

130. Новожилов В.В., Слепян Л.И. О принципе Сен-Венана в динамике стержней. Прикладная математика и механика, 1965, т. 29, № 2, с. 261-281.

131. Норри Д. де Фриз Ж. Введение в метод конечных элементов. М.: Мир, 1981.304 с.

132. Паймушин В.Н. Соотношения теории тонких оболочек типа Тимошенко в криволинейных координатах поверхности отсчета // Прикл. математика и механика. 1978. Т. 42. №4. С. 767-772.

133. Паутов А.Н. Треугольный конечный элемент для анализа изгиба пластин с учетом деформации поперечного сдвига // Прикладные проблемы прочности и пластичности. Статика и динамика деформируемых систем: Всесоюз межвуз. сб./Горьк. ун-т, 1983. С. 95-99.

134. Пацюк В.И., Рыбакова ПА., Сабодаш П.Ф. Волновые процессы в цилиндрической оболочке при неосесимметричном продольном ударе // Прикл. механика. 1985. Т. 21. № 1. С. 35-42.

135. Перце»Платонов;^ ллаетин, Л.:,

136. Судостроение. 1987. 317 с.

137. Пирогов С.А., Чекмарев Д.Т. Явно-неявные численные схемы решения двумерных задач теории упругости//Вестник ННГУ. Серия Механика! Н.Новгород: Изд-во ННГУ. 2000. Вып. 2. С.74-81.

138. Победря Б. Е. Численные методы в теории упругости и пластичности. Изд--во МГУ, 1981.343с;

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

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

141. Прочность, устойчивость, колебания: Справочник Т.З. / Под ред.: И.А.Биргера, Я.Г.Пановко. М.: Машиностроение, 1968.

142. Ракитский Ю.В., Устинов С.М., Черноруцкий И.Г. Численные методы решения жестких систем. М.: Наука, 1979.208 с.

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

144. Рихтмайер Р, Мортон К. Разностные методы решения краевых задач. Мл Мир, 1972.418 с.

145. Рождественский Б. Л., Яненко Н.Н. Системы квазилинейных уравнений и их приложение к газовой динамике. 2-е изд., перераб. и доп. М.: Наука, 1978.

146. Розин Л.А. Метод конечных элементов в применении к упругим системам; М.: Стройиздат, 1977. 129 с.

147. КМК, 1999. С: 122-124.::;;; .:;:; ; :. г:.!. „ i.;.:,'

148. Садовский В.М. Гиперболические вариационные неравенства в задачах динамики упругопластических тел IIПММ. 1991. Т. 55; Вып.6. С. 1041-1048.

149. Садырин А.И. Применение треугольных сеток к решению динамических упругопластических задач // Прикладные проблемы прочности и пластичности. Статика и динамика деформируемых систем: Горький: изд-во Горьк. ун-та, 1983. Вып. 24. С. 39-46.

150. Самарский А.А. Теория разностных схем. М. Наука, 1983. 616 с.

151. Самарский А.А., Гулин А.В. Численные методы. М.: Наука, 1989. 430 с.

152. Самарский А.А.,Попов Ю.П.Разностные схемы газовой динамики.-М.:Наука,1980.

153. Самарский А.А;, Тишкин В.Ф., Фаворский А.П., Шашков М.Ю. Операторные-разностные* схемы // Диф; уравнения: 1981/ Т. 17.N 7: С.1317*

154. Сахаров А.С. Кислоокий В.Н., Киричевский В.В. и др. Метод конечных элементов в механике твердых тел. Киев, Вищ. шк., 1982.479 с.

155. Сахаров А.С. Модификация метода Ритца для расчета массивных тел на основе полиномиальных разложений с учетом жестких смещений. // Сопротивление материалов и теория сооружений. Киев, Будивельник, 1974, Вып. 23, С. 62-71.

156. Слепян Л.И. Нестационарные упругие волны. Л.: Судостроение. 1972.

157. Степаненко М.В.Об одном методе расчета нестационарных импульсных процессов деформирования в упругих конструкциях //Физико-технические проблемы разработки полезных ископаемых. 1976. N 2,C.53-57.

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

159. Теоретические основы и конструирование численных алгоритмов задач математической физики. / Под ред. Бабенко К.И. М.г Наука. 1979.

160. Тимошенко С.П. Колебания в инженерном деле. М.: Наука, 1967.

161. Угодчиков А.Г., Баженов В.Г., Рузанов А.И. О численных методах и результатах решения нестационарных задач теории упругости и . пластичности // Численные методы механики сплошной среды./ СО АН СССР. Т. 16. № 4. Новосибирск. 1985. С. 129-149

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

163. Уилкинс М., Френч С., Сорем М. Конечно-разностная схема для решения задач, зависящих от трех пространственных координат й времени// Численные методы в механике жидкостей. Мл Мир, 1973. С. 115-119.

164. Уилкинс МЛ. Расчет упругопластических, течений // Вычислительные методы в гидродинамике / М.: Мир, 1967. С.212-263.

165. Фельдштейн В.А. Упруго-пластические деформации цилиндрической оболочки при продольном ударе//Волны в неупругих средах / Кишинев, 1970.С. 199-204.

166. Фомин В.М., Хакимов Э.И. Численное моделирование волн сжатия иразрежения в металлах. // Журн. прикл. механики и техн. физики, 1979, №5. 2Ю.Харлоу Ф.Х. Численный метод частиц в ячейках для задач гидродинамики

167. Тимошенко//Прикладные проблемы прочности и пластичности. Численное ,М. Товарищ, науч. изданий

168. Численное исследование процессов высокоскоростного деформирования• металлов на основе метода конечных элементов / В.А. Глущенков, И.Е.Гончаренко, Г.З.Исарович, В.Н.Кислоокий // Машиноведение. 1986. №4.

169. Шокин 10.И.,ЯнентН.Н^Методдифференциальнош-приближеш5я^^^».^ Применение к газовой динамике. Новосибирск: Наука, 1985.364 с.

170. Яненко Н.Н. Метод дробных шагов решения многомерных задач-^математической физики. Новосибирск: Наука. Сибирскоеотд^ 1967.195 с. •

171. Ahmand B.M.r Irons ОС., Zienkiewicz О. Analysis of thick shell structures by curved finite elements // Int. J. Num. Methods Eng. 2. 1970. P. 419 451.

172. Anderson D.L., Lindberg H.E. Dynamic pulse buckling of cylindrical shells due to impulsive loading//AIAA Journal,-1968.V.6.No.4.p.589-598.

173. Argyris J.Hi Doltsinis J.St. In the large strain inelastic analysis in natural formulation. Part II. Dynamic problems. // Сотр. Meth. Appl. Mech. Eng. 1980. V. 21.1 I. P. 91-126.

174. Bathe K.J, Dvorkin E.N. A formulation of general shell elements the use of mixed т1ефо1а11оп of tensorial components//International Journal for numerical methods in engineering. 1986. V. 22. P.697-722.

175. Belytschko Ted, Lin Jerry I. Tsay Shyn-Chin Explicit algorithms for the nonlinear dynamics of shells // Comput. Meth. Appl. Mech. and Engng. 1984. V. 42. №2. P.225-251.

176. Dvorkin E.N., Onate E., Oliver J. On non-linear formulation for curved Timoshenko beam elements considering large displacement/rotation increments// IntJ.Numer.Eng.,1988. Vol.26. №7. 1597-1614.

177. Flanagan D.P., Belytschko T. A uniform strain hexahedron and quadrilateral with orthogonal hourglass control//Intern. J. Num. Meth. Eng. l981. V.17. № 5. P.679-706. .

178. Hartzman W., Hutchinson J.R. Nonlinear dynamics of Solids by Finite Element Method//Int. Сотр. and Struct. 1972. v. 2, N1-2. p. 47-77.

179. Herrman W., Bertolf L.D., Thompson S.I. Computational methods for stress wave propagation in nonlinear solid mechanics // Lfect.'Notees Math. 1975. V. 461.• P. 91-127.

180. Huges T.J.R., Liu W.K. Implicit-explicit finite elements in transient analysis: stability theory // J. Appl. Mech. ASME. 45.1978. P. 371-374

181. Huges T.J.R., Liu W.K. Implicit-explicit finite elements in transient analysis: implementation and numerical examples // J. Appl. Mech. ASME. 45. 1978. P^ 375-378 \.

182. Huges T.J.R., Pister K.S., Taylor R.L. Implicit-explicit finite elements in nonlinear transient analysis. // Comput. Meth. Appl, Mech. Eng. 1979. - V. 17-18,4,-P.159-182. :

183. Kikuchi Fumio. An abstract analysis of parameter dependent problems to mixed finite element methods //J.Fac.Univ. Tokyo, 1985.Sec. IA,32,No.3.p.499-538.

184. Mac Neal R., Harder R.L. A proposed standart set of problems to test finite element accuracy// Finite Element Analysis Design. 1985. №1. P.3-20.

185. Petschek A.G., Hanson M.E. Difference equation for two-dimensional elastic flow//J.Comp. Phys. V.3 1968. № 2. P.307-321.

186. Pawsey S.F., Clough R.W. Improved Numerical Integration of Thick Shell Finite Elements // Int. J. Num. Meth. Eng. 1971. V. 3. № 4 P. 575-586.

187. Stephens W.B., Fulton R.E. Axisymmetric static and dynamic buckling of spherical caps due to centrally distributed pressure// ALAA Journal, 1969. V.7.No.6.p.2120-2126.

188. Takemoto H., Cook R.D. Some modifications of an isoparametric shell element // Int. J. Num. Meth. Eng. V. 7. P. 401-405.

189. Timoshenko S.P. On the correction for shear of the differential equation for transverse vibration of prismatic bars. If Phil. Mag. 1921. V. 41. P. 744-746.

190. Witmer E.A., Balmer H.A., Leech J.W., Pian Т.Н. Large dynamic deformation of beams, rings, plates and shells// AIAA Journal 1963. V. 1, № 8. P. 1848-1857

191. Zeng Q.,. Combescure A. A new one-point quadrature, general non-linear quadrilateral shell element with physical stabilization // Int. J. Numer. Meth. Eng. 1998. №42. P. 1307-1338.

192. Zienkiewich O.C., Too J., Taylor R.L. Reduced integration technique in generalanalysis of plates and shells // IntJ. Num.Meth.Engng. 1971 .V.3.No.2.p.275-290.