Эффективные методы расчета колебаний тонкостенных конструкций, содержащих жидкие и вязкоупругие среды тема диссертации и автореферата по ВАК РФ 01.02.06, доктор технических наук Мокеев, Владимир Викторович

• Мокеев, Владимир Викторович
• доктор технических наук
• 1999

Оглавление диссертации доктор технических наук Мокеев, Владимир Викторович

Введение

1 Обзор методов исследования динамики упругих конструкций с жидкими и вязкоупругими средами

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

1.2 Обзор методов исследования динамики конструкций с жидкостью

1.3 Обзор методов исследования динамики конструкций с вязкоупругими средами

1.4 Цель и задачи работы

2 Метод решения проблемы собственных значений больших матричных систем

2.1 Теоретические основы метода частотной конденсации

2.2 Алгоритм метода.

2.3 Сегментирование фронтальных матриц.

2.4 Выбор основных и вспомогательных степеней свободы

2.5 Надёжность метода.

2.6 Точность метода.

2.6.1 Критерии для оценки точности.

2.6.2 Исследование точности метода на примерах решения задач собственных колебаний.

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

2.7 Быстродействие метода.

Выводы

Метод исследования динамики конструкций с жидкостью

3.1 Основные уравнения и модели описания свойств жидкости

3.1.1 Основные соотношения метода конечных элементов для задач динамики конструкций с жидкостью

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

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

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

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

3.3 Нахождение собственных значений и собственных векторов

3.4 Сравнение результатов с известными решениями.

3.5 Схема приближенного учёта сжимаемости жидкости

3.6 Об эффективности вычислений собственных частот и форм конструкций с жидкостью.

3.7 О точности схемы приближенного учёта сжимаемости жидкости

3.8 Выводы.

Метод исследования динамики конструкций с вязкоупругими средами

4.1 Основные уравнения и модели вязкоупругих материалов

4.2 Идентификация вязкоупругой модели.

4.3 О влиянии температуры на вязкоупругие свойства материала

4.4 Собственные колебания конструкций с вязкоупругими средами

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

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

Особое место в инженерных расчетах занимают методы исследования вибронагруженности конструкций при различного рода динамических воздействиях. Эксплуатация конструкций при интенсивных ударных и вибрационных воздействиях может привести к их разрушению либо нарушению заданных режимов работы. Динамические нагрузки, действующие на ракету при разделении ступеней импульсом избыточного давления, являются не только опасными для приборов системы управления, но и определяют прочность корпуса приборного отсека, корпуса и узлов крепления боевых блоков, топливного заряда и т.п. Для оборудования нефтегазопро-мыслов и нефтегазопереработки серьезную опасность могут представлять различного рода вибрационные и ударные нагрузки, такие как ветровые, сейсмические, воздействие ударной волны при взрывах, падение тяжёлых предметов (элементов бурового оборудования массой до 20 т. и т.д.) при проведении ремонтных работ и др. Подробная и оперативная информация о режимах динамического нагружения требуется на всех этапах процесса проектирования, так как по ней определяются геометрические размеры и материал несущих элементов конструкции. Актуальность темы возрастает ещё и в связи с тем, что интенсивная и продолжительная вибрация оболочек топливных емкостей, цистерн для перевозки токсичных и горючих жидкостей может приводить к авариям и экологическим катастрофам. Большие объёмы перевозимых "экологически грязных жидкостей" предъявляют особые требования к прочности и надёжности транспортных средств, поскольку последствия загрязнения окружающей среды становятся все более необратимыми, а расходы на их устранение и восстановление естественных природных условий всё обременительнее.

Задачами исследования динамики тонкостенных конструкций с жидкими и вязкоупругими средами занимались многие видные российские и зарубежные ученые. Большой вклад в разработку теоретических основ современных методов расчета тонкостенных конструкций внесли работы Н.А Алфутова, H.A. Амбарцумяна, Л.И. Балабуха, В.В. Болотина, А.Н. Гузя, И.И. Гольденблата, В.Т. Лизина, И.Ф. Образцова, В.И. Усюкина и многих других. Успешное исследования динамики конструкций возможно, если решены две задачи. Первая задача заключается в создании модели конструкции, вторая - в решении уравнений движений, полученных для разработанной модели. Модель конструкции должна учитывать геометрическую форму всей конструкции и отдельных ее элементов, реальные изменения механических и массовых свойств, позволять увеличивать степень её детализации и сложности в процессе проектирования, позволять находить не только кинематические параметры, но и характеристики напряженно-деформированного состояния элементов конструкции. Наиболее полно этим требованиям удовлетворяют модели, построенные на основе метода конечных элементов (МКЭ). Широкому распространению, теоретическому обоснованию и развитию МКЭ способствовали работы российских и зарубежных ученых Постнова В.А. [72], Розина JI.A. [74], Корнеева В.Г.[31], Зенкевича О. [28, 29], Стренга Г., Фикса Дж. [76], Галагера Р. [19] и др. В указанных работах установлена связь МКЭ с другими численными методами, разработана общая методология построения конечно-элементных моделей. Использование конечно-элементных моделей, учитывающих особенности динамического взаимодействия тонкостенных конструкций с жидкими и вязкоупругими средами, позволяет повысить надежность проектируемых конструкций, снизить их материалоемкость и повысить качество всего процесса проектирования. Расчет динамического нагруженности конструкций является составной частью таких сложных проектных задач как: многовариантный анализ, оптимизация, структурный синтез. В этих случаях расчет динамической нагруженности конструкции приходится выполнять многократно. При этом высокая эффективность используемых методов исследования динамики конструкций есть необходимая предпосылка успешного выполнения многовариантного проектного анализа и оптимизации сложных технических объектов. Анализ наиболее распространенных конечно-элементных программ, используемых при проектировании конструкций, показывает, что в большинстве из них возможности расчета колебаний конструкций с жидкими и вязкоупругими средами либо отсутствуют либо ограничены и не обладают достаточной надежностью и эффективностью, что объективно отражает состояние развития методов исследования в данной области.

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

Задачи моделирования динамического поведения топливных баков ракетных конструкций в полете и при наземной эксплуатации, ёмкостей и цистерн при транспортировке больших объёмов жидкостей, оборудования нефтегазопромыслов приводят к необходимости исследовать колебания системы: "упругая конструкция - жидкость - поверхностные волны". Решение задач динамики конструкций с жидкостью рассматривалось в работах как российских учёных И.Б. Богоряда, A.C. Вольмира, А.Г. Горшкова, М.А. Ильгамова, К.С. Колесникова, Г.Н. Микишева, H.H.

Моисеева, Ю.С. Павлюка, Б.И. Рабиновича, В.В. Румянцева, К.В. Фролова, Ф.Н. Шклярчука, В.П. Шмакова, так и зарубежных - H.N. Abramson, H.F. Bauer, D.D. Капа, О.С. Zienkiewicz и др.

Одна из основных проблем, возникающих при построении моделей твердотопливных ракет, оболочек с демпфирующими покрытиями, конструкций из полимеров и пластмасс и т.п., связана с тем, что элементы этих конструкций обладают вязкоупругими свойствами. В настоящее время теория вязкоупругости достаточно хорошо разработана в трудах российских и зарубежных ученых Работнова Ю.Н. [73], Москвитина В.В. [65], Кристенсена P.M.[39], Бленда Д. [14], Колтунова М.А. [41] и др., в которых для описания вязкоупругих свойств используется широкий набор моделей, базирующихся на интегральном и дифференциальном исчислении. Решение некоторых задач расчета конструкций из вязкоупругих материалов приведено в работах К.Д. Джонсона, Д.А. Кинхольца, В.П. Майбороды, М.И. Мирсаидова, В.Г. Пальчиковского, Я.С. Садикова, И.Е. Трояновского и др.

Важным элементом динамического анализа является решение уравнений, описывающих в рамках выбранной модели колебания конструкции. Основные трудности здесь обусловлены высоким порядком матриц, полученных при конечно-элементной дискретизации конструкции. Из существующих в настоящее время методов решения матричных уравнений наиболее часто применяется модальный метод, использующий разложение вектора перемещений конструкции в ряд по собственным формам. При этом задача сводится к проблеме нахождения собственных значений матричных уравнений большого порядка. Решению этой проблемы посвящены работы Кандидова В.П. [32], Кислоокого В.Н. [36, 37], Банах Л.Я.[9, 10], Бате К.Д. и Вилсона E.JI.[11], Гайана Р [18], Андерсона JI.P. и Халлоуэра B.JI. [2] и др. В отличие от упругой задачи матрицы, получаемые при исследовании конструкций с вязкоупругими и жидкими средами, являются неположительно определенными, что предъявляет к надежности метода нахождения собственных значений довольно жесткие требования.

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

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

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

- создание эффективного метода исследования динамики конструкций с жидкостью;

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

- создание программ, реализующих разработанные методы;

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

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

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

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

Новыми результатами в диссертации, которые автор выносит на защиту, являются: ;

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

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

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

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

- при оценки динамической нагруженности твердотопливных и жидкостных ракет в Государственном ракетном центре "КБ им. академика В.П.Макеева,"

- для определения динамических характеристик топливных баков ракето-носителя "Энергия" в Волжском филиале НПО "Энергия" (г. Самара),

- при проектировании энергонапряженных машиностроительных конструкций в КБ ПО "Полет" (г. Омск), КБ "Арсенал" (Санкт-Петербург) и других предприятиях. I

Данная фа/бота поддержана грантами Министерства общего и профессионального образования РФ по фундаментальным исследованиям в области транспортных наук (96-10-4.1-15), машиностроения (97-24-12.3-1 ), в области фундаментальных проблем охраны окружающей среды и экологии человека (проект 1-21).

Заключение диссертации по теме «Динамика, прочность машин, приборов и аппаратуры», Мокеев, Владимир Викторович

4.9 Выводы

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

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

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

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

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

6. Проведено изучение влияния на амплитуды вынужденных колеба

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

Глава 5

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

5.1 Исследование продольных колебаний топливных баков жидкостных ракет

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

Рассматриваются собственные колебания полусферической бака, заполненного водой и шарнирно закрепленного по верхнему краю А (рис.5.1). Параметры оболочки: радиус 0,133м, толщина оболочки 0,0007ж, модуль упругости 4,016 • 109Да, плотность материала 1180кг/ж3, коэффициент Пуассона 0,3. Подобная оболочка была аналитически и экспериментально исследована в работе [85]. При проведении эксперимента оболочка закреплялась по контуру кольцами к стакану, который крепился к подвиж

Конечно-элементная модель полусферической оболочки с водой А

Рис. 5.1. ному столу электродинамического вибратора. Для определения частот на поверхности полусферы в 2-х плоскостях наклеено четыре проволочных датчика сопротивления. Усиление сигналов датчика проводилось универсальным тензометрической установкой УТС1-ВТ-12. Колебания записывались шлейфовым осциллографом К-9-21 [85]. В табл. 5.1 представлены значения первых трёх частот осесимметричных колебаний оболочки, полученные в аналитическим и экспериментальным методами в работе [85].

Заключение

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

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

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

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

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

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

Список литературы диссертационного исследования доктор технических наук Мокеев, Владимир Викторович, 1999 год

1. Абросимов H.A., Баженов В.Г. К расчету стеклопластиковых цилиндрических оболочек при ударном нагружении.//Механика композитных материалов. - Рига, 1983.- N 5. - С. 820-823.

2. Андерсон JI.P., Халлоуэр В.Л. Метод понижения порядка в задачах динамики конструкций, основанный на итерационном решении уравнения Риккати. //Ракетная техника и космонавтика. 1981, 19, N 8. -С. 147-153.

3. Анисимов A.M. Осесимметричные колебания сферической оболочки, заполненной жидкостью . //Изв. Вузов, Авиационная техника. 1963, N 2.

4. Ахенбах Д. Динамическое поведение длинного скрепленного с корпусом вязкоупругого цилиндра. //Ракетная техника и космонавтика. 1965, - N 4. - С. 124-129.

5. Айронс Б.М. Задачи о собственных значениях: исключение лишних переменных. //Ракетная техника и космонавтика. 1965, 3, N 5. -С. 207-209.

6. Балабух JIM., Ганичев А.И., Молчанов А.Г. Две задачи о собственных колебаниях упругих систем с жидким наполнением. В сб.: Расчеты на прочность. М.: Машиностроение, 1966, вып. 12.

7. Балабух Л.И., Молчанов А.Г. Осесимметричные колебания сферической оболочки, частично заполненной жидкостью . //Инж. ж. Механика твердого тела. 1967, N 5.

8. Балакирев Ю.Г., Шмаков В.П. Осесимметричные колебания цилиндрической оболочки с полусферическим днищем //Колебание упругихконструкций с жидкостью: Сборник научных докладов 3 симпозиума. ЦНТИ "Волна", 1974. С. 28-31.

9. Банах Л.Я. Упрощение расчетных динамических систем. //Колебания и динамическая прочность элементов машин. М: Наука, 1976. - С. 36-46.

10. Банах Л.Я. Уменьшение числа степеней свободы при исследование многомерных систем. //Машиноведение. 1979, N 1. - С. 21-26.

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

12. Бегли Р.Л., Торвик П.Д. Дифференциальное исчисление, основанное на производных дробного порядка новый подход к расчету конструкций с вязкоупругим демпфированием. //Аэрокосмическая техника. - 1984, 2, N 2. - С. 84-93.

13. Белов И.В., Мокеев В.В. Исследование динамики вязкоупругих композиционных материалов с использованием метода конечных элементов. //Техника, экономика, информация. Сер. Конструкции из композитных материалов. М., 1982, вып. 1. - С. 32-35.

14. Бленд Д.Р. Теория линейной вязкоупругости. М.: Мир, 1965. - 197 с.

15. Богданов Ф.И., Халмагпов Г.А. Алгоритмы и программы диспетчера для решения геометрически нелинейных задач теории вязкоупругости оболочек и пластин. //Алгоритмы и программы. Ташкент, 1974, вып. 16. - С. 9-32.

16. Болотин В.В. Колебания и устойчивость упругой цилиндрической оболочки в потоке сжимаемой жидкости. Инж. сб., 1956, т.24.

17. Волъмир A.C., Михнев В.Ф., Тихомиров А.Е., Хорхордин В.И. Многоуровневые методы в динамике роторов авиационных двигателей. //Прикладная механика. 1984, 20, N 12. - С. 58-63.

18. Гайан Р.Д. Приведение матриц жесткости и масс. //Ракетная техника и космонавтика. 1965, N 2. С. 287-289.

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

20. Графтон П.Е. Строум Д.Р. Расчет осесимметричных оболочек методом прямого определения жесткости. //Ракетная техника и космонавтика. 1963, 1, N 10. С. 129-130.

21. Григорьев В.Г. Расчет динамических характеристик сложных оболо-чечных конструкций с жидкостью. //Колебание упругих конструкций с жидкостью: Сборник научных докладов 4 симпозиума. ЦНТИ "Волна", 1976. С. 386-396.

22. Гриненко Н.И., Мокеев В.В. Использование частотнодинамичес-кой конденсации при исследовании колебаний сложных конструкций методом конечных элементов. //Колебания упругих конструкций с жидкостью: Сб. науч. докладов 5 симпозиума. М., 1984. - С. 100102.

23. Гриненко Н.И., Мокеев В.В. Применение теории нормальных форм при решении задач динамики конструкций, обладающих вязкоуп-ругими свойствами. //Техника, экономика, информация. Сер. Конструкции из композитных материалов. М., 1985, вып. 2. - С. 10-15.

24. Гриненко Н.И., Мокеев В.В. Анализ динамики неоднородных вяз-коупругих конструкций методом конечных элементов. // Механика композиционных материалов. 1989, N 3. С. 475-481.

25. Гриненко Н.И., Мокеев В.В. О задачах исследований колебаний конструкций методом конечных элементов. //Прикладная механика. -1985, 21, N 3. С. 25-30.

26. Дауне К. Точная редукция матрицы масс и жесткости при анализе колебаний и рациональные основы выбора основных степеней свободы. //Труды Амер. общества инженеров-механиков. Конструирование и технология машиностроения. 1980, 102, N 2. - С. 132-138.

27. Джонсон К.Д., Кинхолъц Д.А. Расчет демпфированных колебаний в конструкциях, содержащих вязкоупругие слои, методом конечных элементов. //Аэрокосмическая техника. 1983, 1, N 4. - С. 120-125.

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

29. Зенкевич О.С. Метод конечных элементов в технике. М.: Мир, 1975.- 541 с.

30. Ицкович М.А., Нигул У.К. Эволюция нестационарной волны в наследственных упругих средах. //Механика твердого тела.- 1984, N 1.- С. 72-77.

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

32. Кандидов В.П. Расчет устойчивости изгибно-крутильных колебаний крыла в дозвуковом потоке методом конечных элементов. //Уч. записки ЦАГИ. 1972. 3, N 1. - С. 111-117.

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

34. Кауфман СХолл Д.В. Приведение матриц масс и нагрузки. //Ракетная техника и космонавтика. 1968, 6, N 3. - С. 235-237.

35. Киддер Р.Л. Приведение частотных уравнений для конструкций. //Ракетная техника и космонавтика. 1976, 16, N 6. - С. 156-157.

36. Кислоокий В.Н., Гранат С.Я. Исследование колебаний пластин и комбинированных систем методом базисных векторов на основе разностной дискретизации. //Сопротивленние материалов и теория сооружений. Киев: Будевильник, 1979, вып. 16. - С. 190-194.

37. Кислоокий В.H., Легостаев А.Д. Реализация МКЭ в задачах исследования колебаний оболочек и пластин. //Сопротивление материалов и теория сооружений. Киев: Будевильник, 1979, вып. 24. - С. 25-32.

38. Кифлинг Л., Фен Г. Расчет колебаний конструкций с жидкостью методом конечных элементов. //Ракетная техника и космонавтика. 1964, N 7. С. 89 - 94.

39. Кристенсен P.M. Введение в теорию вязкоупругости. М.: Мир, 1974. - 328 с.

40. Колтунов М.А. Ползучесть и релаксация. М.: Высшая школа, 1976. -297.

41. Колтунов М.А., Майборода В.П., Зубчанинов В.Г. Прочностные расчеты из полимерных материалов. М.: Машиностроение, 1983. - 420 с.

42. Колтунов М.А., Мирсаидов М.И., Трояновский И.Е. Установившиеся колебания осесимметричных вязкоупругих оболочек. //Механика полимеров. 1978, N 2. - С. 290-295.

43. Кокошвили С.М., Калнин П.П., Янсон Ю.О. Квазистатические методы исследования жестких полимеров. //Механика полимеров. Рига, 1974, N 4. - С. 683-688.

44. Кузнецов Г.В., Матвиенко В.П., Трояновский И.Е. Вынужденные колебания короткого вязкоупругого цилиндра. //Вопросы теории упругости и вязкоупругости.- Свердловск, 1978. С. 34-38.

45. Кухар Е.Д., Сталь C.B. Метод динамического преобразования в синтезе форм колебаний. //Ракетная техника и космонавтика. 1974, 12, N 5. - С. 120-129.

46. Лампер P.E. Осесимметричные колебания цилиндрического бака. //Изв. Вузов, Авиационная техника. 1967, N 4.

47. Лизин В. Т., Пяткин В.А. Проектирование тонкостенных конструкций. М.: Машиностроение, 1976. - 408 с.

48. Майборода В.П., Трояновский И.Е. Собственные колебания неоднородных вязкоупругих тел. //Механика твердого тела. 1983, N 2. -С. 117-123.

49. Малмейстер А.К. Статистическая интерпретация реологических уравнений. //Механика полимеров. Рига, 1966, N 2.- С. 197-213.

50. Микишев Г.Н. Рабинович Б.И. Динамика тонкостенных конструкций с отсеками, содержащими жидкость. М.: Машиностроение, 1971. -463 с.

51. Мнев E.H. Колебания круговой цилиндрической оболочки погруженной в замкнутую полость, заполненную сжимаемой жидкостью. В сб.: Теория пластин и оболочек. - Киев, АН УССР, 1962.

52. Мокеев В.В. О решении задач динамики вязкоупругих тел с использованием метода конечных элементов. //Динамика и прочность конструкций: тематич. сб. науч. тр. Челябинск: ЧПИ, 1982. - С. 4752.

53. Мокеев В.В. Применение теории нормальных форм для анализа вынужденных колебаний неоднородных вязкоупругих конструкций. //Динамика машин и рабочих процессов: тематич. сб. науч. тр. -Челябинск: ЧПИ, 1983. С. 87-92.

54. Мокеев В.В. Исследование установившихся колебаний неоднородных вязкоупругих конструкций с помощью метода конечных элементов. //Вопросы динамики и прочности механических систем: тематич. сб. науч. тр. Челябинск: ЧПИ, 1984. - С. 13-24.

55. Мокеев B.B. О задаче нахождения собственных значений и векторов больших матричных систем. //Журнал Вычислительной Математики и Математической Физики. 1992, 32, N 10. - С. 1652-1657.

56. Мокеев В.В., Павлюк Ю.С. Эффективная процедура решения задач о собственных значениях при исследованиях взаимодействия конструкция-жидкость но основе конечноэлементных моделей. //Механика твёрдого тела. 1992, N4. - С. 178-182.

57. Мокеев В. В., Фот Е.Я. Метод частотной конденсации эффективный подход в решении задач динамики конструкций //Динамика, прочность и износостойкость машин". 1997, N 3. - С. 63 - 67.

58. Мокеев В. В. Колебания упругих систем с акустическими средами // Динамика, прочность и износостойкость машин". 1997, N 3. -С. 58 - 62.

59. Мокеев В. В., Фот Е.Я. Модальный анализ гидроупругого взаимодействия конструкций с жидкостью и газом //Проблемы проектирования неоднородных конструкций. Труды 15 Российской школы. 1996. 93-97 с.

60. Mokeyev V. V. A frequency condensation method for the eigenvalue problem //Communication in Numerical Method in Engineering. 1998, v.14. — P. 1-8.

61. Мокеев В. В. Моделирование вязкоупругого поведения конструкций из композиционных материалов при вибрациях и ударах. //Механикакомпозиционных материалов и конструкций. 1997, т.З, N 4. - С. 321.

62. Мокеев В. В. Моделирование динамического поведения конструкций из вязкоупругих материалов с помощью метода конечных элементов. // Динамика, прочность и износостойкость машин. 1998, N 4. -С. 44-50 .

63. Москвитин В.В. Сопротивление вязкоупругих тел. М.: Наука, 1974. - 328 с.

64. Моргунов Б.И. О колебаниях длиной вязкоупругой пластинки. //Труды Моск. энерг. института. Механика. -1974, вып. 39. С. 121129.

65. Нашиф А., Джоунс Д., Хендерсон Дж. Демпфирование колебаний М.: Мир, 1988.

66. Оялво А.У., Нъюмак М. Формы колебаний больших систем, полученные методом автоматической редукции матриц. //Ракетная техника и космонавтика. 1970, 8, N 7. - С. 62-68.

67. Парлетт Б. Симметричная проблема собственных значений. М.: Мир, 1983.

68. Пожалостин A.A. К расчету собственных колебаний пологой сферической оболочки. //Изв. Вузов, Машиностроение, 1963, N 10.

69. Пожалостин A.A. Осесимметричные колебания упругого бака с жидкостью. //Изв. Вузов, Аввиационная техника, 1964, N 3. С. 148-151.

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

71. Работное Ю.Н. Элементы наследственной механики твердых тел. -М.: Наука, 1977. 384 с.

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

73. Самойлов Е.А., Павлов Б.С. Колебаний полусферической оболочки, заполненной жидкостью . //Изв. Вузов, Авиационная техника. 1964, N 3.

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

75. Сорокин Е. С. К теории внутреннего трения при колебании упругих систем. М.: Госстройиздат, 1960. - 450 с.

76. Сорокин Е. С. Еще раз об учете внутреннего сопротивления при исследовании случайных колебаний конструкций. //Строительная механика и расчет сооружений. 1971, N 1. - С. 10-15.

77. Сорокин Е.С. Уравнения динамической теории упругости с учетом внутреннего трения. //Труды МИИТ. 1971, вып. 343 - С. 121-132.

78. Скучик Е.А. Простые и сложные колебания систем. М.: Мир, 1971. - 557 с.

79. Системы автоматизированного проектирования: В 9-ти кн. Кн. 1. Принципы построения и структура. Учебное пособие для втузов. / И.П. Норенков. М.: Высшая школа, 1986. - 127 с.

80. Садиков Я.С., Палъчиковский В.Г., Булъбович Р.В., Жерноклеева Г. И. Исследование зависимости динамического модуля упругости некоторых смесевых твердых топлив от частоты и деформации. //Оборонная техника. 1975, N 9. - С. 12-14.

81. Садиков Я.С., Палъчиковский В.Г., Булъбович Р.В., Жерноклеева Г. И. Определение основной резонансной частоты двухслойной конструкции, содержащей вязкоупругий заполнитель. //Оборонная техника. 1975, N 11. - С. 12-14.

82. Сахаров И.Е. Экспериментальные исследования и расчет собственных колебаний оболочек вращения с набором кольцевых ребер. //Теория оболочек и пластин. М.: Наука, 1975. - С. 560-565.

83. Самойлов Е.А., Павлов Б.С. Колебания полусферической оболочки заполненной с жидкостью. //Известия ВУЗов. 1964, - N 3. - С. 1015.

84. Сиделъников Р.В., Ямчук В.В. Расчет колебаний осесимметричных конструкций с жидкостью методом конечных элементов. //В сб.: Динамика и прочность конструкций. Сборник научных трудов N 227. -Челябинск: ЧПИ, 1979. С. 24 - 29.

85. Сорокин Е.С. О погрешностях общеизвестного метода теории колебаний диссипативных систем в применении к неоднородному демпфированию. //Строительная механика и расчет сооружений. 1984, N 2. - С. 29-34.

86. Coy ер с Д. Уплотнение матрицы масс свободного тела с помощью коэффициентов податливости. //Ракетная техника и космонавтика. 1976, 16, N 3. - С. 101-102.

87. Уржумцев Ю.С., Янсон Ю.О. О паспортизации вязкоупругих характеристик полимерных материалов. //Механика композиционных материалов 1979, N 5. - С. 900-907.

88. Фрид А. Конденсация в задачах о собственных значениях, решаемых методом конечных элементов. //Ракетная техника и космонавтика. -1972, 10, N 11. С. 185-186.

89. Филатов A.M. Асимптотические методы в теории дифференциальных уравнений. Ташкент: ФАИ, 1974. 215 с.

90. Фролов К.В., Антонов В.Н. Колебания оболочек в жидкости. М.: Наука, 1983. 143 с.

91. Хант Д. Идеализация жидкости структурными элементами. //Ракетная техника и космонавтика. 1971, N 3. 136 142 с.

92. Харрис С. М., Крид Ч.И. Справочник по ударным нагрузкам. JL: Судостроение, 1980 360 с.

93. Хессиг Д. Исключение степеней свободы с помощью структурного сплининга. //Ракетная техника и космонавтика. 1973, 11, N 3. -С. 178-179.

94. Шкллрчук Ф.Н. Осесимметричные колебания жидкости внутри упругой цилиндрической оболочки с упругим днищем. //Изв. Вузов, Авиационная техника, 1965, N 4.

95. Шмаков В.П. Об уравнениях осесимметричных колебаний цилиндрической оболочки с жидким заполнением. //Изв. АН СССР, Механика и машиностроение. 1964, N 1.

96. Шмаков В.П. Об одном методе интегрирования уравнений осесимметричных колебаний сферической оболочки . //Инж. ж. Механика твердого тела. 1966, N 2.

97. Шкллрчук Ф.Н. Колебания упругой оболочки, содержащей тяжелую сжимаемую жидкость. //Сборник научных докладов 3 симпозиума. ЦНТИ "Волна", 1976. С. 386-396

98. Юй Подпространственное итерирование в задачах о собственных значениях при исследовании гидроупругих взаимодействий. //Теоретические основы инженерных расчетов. 1988, N 1. С. 176-182.

99. Archer J.S. Consistent mass matrix for distributed mass system. //Proc. of the Amer. Society of Civil Eng.; J of Struct. Duv. 1963, 89, N 4. -P. 17-20.

100. Bleich H.H., Baron M.L. Free and forced vibration of infinitely long cylindrical shell an infinite acoustical medium. //J. Appl. Mech. 1954, 21, N 2. - P. 167-177.

101. Crandel S. H. Dynamic responses of system with structural damping. //Air. Space and Instruments. New-York: Mc.Craw-Hill Book Co. -1963. - P. 183-193.

102. Chopra A.K., Chakrabarti P. Eartquake analysis of concrete gravity dams including dam-water -foundation rock interaction. //Earthquake Eng. Struct. Dynam. 1982, 10. - P. 363-383.

103. Dubois J.J., de Rouvray A.L. An improved fluid superelement for the coupled solid-fluid-surface wave dynamic interaction problem. //Earthquake Eng. Struct. Dynam. 1978, 6. - P. 235-245.

104. Daniel W.J. T. Peformance of reduction methods for fluid- structure and acoustic eigenvalue problems. //Int. J. Num. Meth. Engng. 1980, 15. - P. 1585-1594.

105. Daniel W.J.T. Modal mehtods in finite element fluid-structure eigenvalue problems. //Int. J. Num. Meth. Engng. 1980, 15. - P. 15611575.

106. Gegodin L., Gallagher K. A frontal-based solver for frequency analysis. //Int. J. Num. Meth. Engng. 1978, 12, N 11. - P. 1659-1666.

107. Gilroy L.E. Coupled finite element analyses of sea shest vibration. //Noise Contr. Eng. 1994, 42, N 1. - P. 19-24.

108. Gladwell G.M., Zimmermann G. On energy and complementary energy formulations of acoustic and structural vibration problems. //J. Sound and Vib. 1966, 3. - P. 233-241.

109. Hall J.F., Chopra A.K. Two-dimensional dynamic behavior of concrete and embankment dams including hydrodynamical effects. //Earthquake Eng. Struct. Dynam. 1982, 9. - P. 305-332.

110. Irons B.M. A frontal solution program finite element analysis. //Int. J. for Num. Meth. in Eng. 1970, N 2. - C. 5-32.

111. Jennigs A. Mass condensation and simultaneous iteration for vibration problem. //Int. J. for Num. Meth. and Eng. 1973, N 6. - P. 543-552.

112. Kana D., Gormley J. Longitudinal vibration of model vehicle propelant tank. //J. Spacecraft and Rockets. 1967, 4, N 12. - P. 50-56.

113. Kawabe H, Kuwahare K. Some analyses for dynamic characteristics of the materials of relaxation type in term of the 3-lumped parameter vibration model. //Trans. Jap. Inst. Metals. 1979, 20, N 5. - C. 209-218.

114. Muller W.C. Simplified analysis of linear fluid-structure interaction. //Int. J. Num. Meth. Engng. 1981, 17. P. 113-121.

115. Nagaga K. Vibration of viscoelastic plate having a circular ounter boundary for varios edge condition. //J. of Sound and Vibration. 1979, 63, N 1. - P. 73-85.

116. O'Callahan J.C. A procedure for an improved reduced system (IRS) model. //Procceeding of the 15th International Modal Analysis Conference, (Las Vegas, NV), Union College, Schenectady, NY, 1989. P. 17-21.

117. Oden J.T., Zienkiewicz O.C., Gallager R.H., Taylor C. Finite elements in fluids. Vol I and II, Wiley, New York. 1975.

118. Papdopoulos M., Garcia E. Improvement in model reduction schemes using the system equivalent reduction expansion process. //AIAA Journal. 1996, v 34, N 3. - P. 2217-2219.

119. Shah V.N., Raymund M. Analytical selection of masters for the reduced eigenvalue problem // Int. J. Numer. Methods. Eng. 1982, 18. - P. 8998.

120. Zienkiewicz O.C., Bettes P. Fluid-structure dynamic interaction and wave forces, an introduction to numerical treatment. //Int. J. Num. Meth. Engng. 1978, 13. - P. 1-16.

121. Wilson E.L. Finite elements for the dynamic analysis of fluid-solid system. //Int. J. Num. Meth. Engng. 1983, 19. - P. 1657-1668.7 / /? ¡р Г)

122. ДИНАМИЧЕСКИМ АНАЛИЗ КОНСТРУКЦИИ КОМПЛЕКС ПРОГРАММ ДИАС

123. Описание комплекса программо г т>оупти С ¿С и у1. Челябинск 1999 г.

124. ПЛ. Структура комплекса программ ДИАС

125. П. 1.1. Назначение комплекса программ ДИАС

126. При создании комплекса ДИАС основными направлениями разработки являлись:

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

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

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

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

131. П.1.2. Процессоры комплекса

132. Протокол работы процессора включает сообщения о начале и завершении работы процессора, все сообщения об ошибках.

133. П.1.3. План вычислительного процесса (ПВП).

134. Обращение к процедуре запуска ПВП имеет вид

135. Имя ПВП путь.имя-данных [/HELP] [/RE:number] [/NOL] [/Fil] [имя-листинга]1. Здесь:

136. RE:mimber — режим рестарта; number номер процессора, с которого надо начать выполнение;

137. NOL отказ от вывода листинга на диск;

138. F11 вывод таблиц осуществляется в 11-ом формате (по умолчанию таблицы выводятся в 12-ом формате);путь.имя-данных имя файла данных;имя-листинга. имя файла с листингом; если имя не указано, то совпадает с именем файла данных;

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

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

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

142. П.1.4. Служба рабочих файлов и диагностика ошибок

143. П. 1.5. Вспомогательная служба

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

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

146. П.2. Подсистемы комплекса ДИАС.

147. Комплекс программ ДИАС состоит из подсистем, в которых сгруппированы программы, ориентированные на решение однотипных задач динамического анализа. В табл. П.1 дано перечисление этих подсистем с указанием какие программы в них входят.