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

• Малинин, Михаил Юрьевич
• кандидат физико-математических наук
• 1984

Оглавление диссертации кандидат физико-математических наук Малинин, Михаил Юрьевич

1. Применение метода конечных элементов к расчету многослойных оболочек ^

1.1. Основные допущения, принятые при расчете многослойных пологих оболочек

1.2. Основные соотношения метода конечных элементов Z

1.3. Использование конечных элементов пластины и пологой оболочки в расчетах оболочек сложной геометрии

1.4. Способы повышения точности расчетов методом конечных элементов для пластин и оболочек типа Тимошенко

2. Матрица жесткости треугольного конечного элемента слоистой пологой оболочки №

2.1. Основные соотношения треугольного конечного элемента пологой оболочки с шестью узлами ^

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

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

2.4. Численное исследование точности расчета оболочек методом конечных элементов 65*

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

3.2. Исследование условий сходимости решения 2i

3.3. Тестовые задачи расчета оболочек с применением четырехугольных элементов ^

4. Сглаживание разрывных функций методом конечных элементов

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

4.2. Аппроксимация с помощью треугольных элементов S

4.3. Аппроксимация с помощью четырехугольных элементов 92.

4.4. Исследование точности решения при использовании аппроксимирующих функций 9Н

4.5. Применение сглаживания в задачах расчета оболочек

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

5.1. Расчет двухслойной пологой оболочки с отверстием 1Q

5.2. Автоматизация подготовки исходных данных и графического вывода результатов расчета пространственных конструкций ^05"

5.3. Расчет тормозного барабана автомобиля КамАЗ-5511 443 • 5.4. Построение итерационного процесса определения толщин по элементам конструкции

5.5. Расчет несущей балки рамы автомобиля КамАЗ-55ПЗ

5.6. Расчет балки передней оси грузовых автомобилей W

5.7. Расчет кронштейна крепления поперечины Л 3 автомобиля КамАЗ

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

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

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

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

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

Теоретические основы расчета многослойных конструкций заложены в работах Амбарцумяна С.А., Александрова А.Я., Болотина В.В., Григолюка Э.И., Григоренко Я.М., Галимова К.З., Королева В.И.* Куршина Л.ГЛ., Муштари Х.М., Новичкова Ю.М., Пелеха Б.Л., Прусакова А.П., Саченкова А.В., Тетерса Г.А., Чулкова А.П. и др. При этом весьма эффективным в практических расчетах показал себя метод гипотез, когда некоторые допущения о напряженно-деформированном состоянии по толщине оболочки или пластины позволяют свести трехмерную задачу к двухмерной.

Можно выделить два направления в применении метода гипотез. Первое направление основано на применении кинематических гипотез для каждого отдельного слоя, когда порядок системы уравнений зависит от числа слоев. Значительное развитие указанное направление в теории слоистых сред получило в работах Болотина В.В. и Новичкова Ю.Н. В работе [ 11] на основе единого подхода предложена теория многослойных конструкций при произвольном числе слоев и широких допущениях о свойствах отдельных слоев. Сравнительный анализ и развитие теорий для многослойных оболочек переменной толщины, которым соответствуют различные предположения об изменении сдвиговых деформаций и напряжений, даны в работах Григоренко Я.М. [■18,19] Л Подобный подход позволяет получить решения для задач со значительным изменением параметров упругости по слоям.

В тех случаях, когда свойства материалов по слоям разнятся незначительно, весьма эффективным оказывается второе направление метода гипотез, когда определенные кинематические допущения принимаются для всего пакета слоев в целом. В этом случае порядок разрешающих уравнений не зависит от числа слоев. Широкое развитие в этом направлении получила модель С.П.Тимошенко. Для пластин неклассическая теория была построена Э.Рейсснером [96] и Р.Миндлиным [S4] . Обобщение модели и ее развитие для оболочек можно найти в работах Галимова К.З., Саченкова А.В. и др. [141 , Пелеха Б.Л. Цурпала И.А.,

Тамурова Н.Г. [6 3] . Развитие модели С.П.Тимошенко в применении к расчету устойчивости и оптимизации многослойных оболочек из композитных материалов приведено в работах Тетерса Г.А.,

Рикардса Р.Б. Г. Расчет реальных многослойных конструкций сложной геометрии позволяют сделать теоретические разработки, предложенные в работах Чулкова П.П. и Паймушина В.Н. [39, <<0,64].

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

За сорок лет (начиная с работ Сига.п( R.C731, Ръадег I Тиъпег М. CtozJ ) ЖЭ получил широкое развитие как за рубежом, так и в отечественной литературе в работах Розина Л.А., Постнова В.А., Шапошникова Н.Н., Григоренко Я.М., Корнишина М.С., 13урмана З.И., Сахарова А.С. и др. Вариационные принципы, на которых базируется метод, могут быть сведены в общую схему, приведенную на рис.1. В верхней части схемы представлены традиционные вариационные принципы.

Принцип минимума потенциальной энергии ' (функционал ХР ) содержит в качестве искомой функции поле перемещений U Перемещения и аппроксимируются пробными функциями А ,

Традиционные вариационные принципы

Принцип Лагранжа u=At J

--[

Совместная модель I

JP(q)

Принцип Ху-Вашицу

Принцип Рейсонера

6"= pfi

Смешанная модель 2

Принцип Кастильяно hrc ссЬ)*0

Равновесная модель I Тс(Р)

Модифицированные вариационные принципы

Принципы Лагранжа cfVr.pi (ч>т а=Сы Т=Мр

Гибридная модель перемещений I

Х.01 (Ы. Р^Х^Ф

Гибридная модель перемещений 2

Принцип Кастильяно

Wu

Гибридная модель .перемещений 3 трз (Ъ)

Гибридная модель напряжений

Жтс (Atf—JiV)

Равновесная модель 2

Смешанный метод

61[Р »T0i%J

Рис.1 Диаграмма вариационных принципов метода конечных элементов где в качестве варьируемых коэффициентов фигурируют узловые обобщенные перемещения ^ . Окончательное расчетное уравнение можно представить в форме :

I) где К - матрица жесткости конструкции

Q - Еектор обобщенных узловых сил.

К уравнению (I) приводит известный метод перемещений [ 16, 2 6, 2 7, 45,50].

Наиболее общий принцип Ху-Вашицу, основанный на варьировании перемещениями, деформациями и напряжениями, позволяет получить принцип Рейсснера (функционал ^ ), базирующийся на вариации перемещений и напряжений [50, 68 ] .В этом случае конечноэлементная формулировка предполагает аппроксимацию пробными функциями как поля перемещений, так и поля напряжений. В качестве искомых величин выступают узловые значения перемещений Ц. и узловые значения напряжений р . Окончательное уравнение смешанного метода, который вытекает из функционала 7ГН , имеет вид : р' U0}

О LQJ где И, б, Q - матрицы, полученные путем суммирования соответствующих матриц конечных элементов. Глобальная матрица в уравнении (2) не является положительно определенной и требует специализированных программ для решения системы уравнений.

Однако принцип Рейсснера можно представить в форме, когда не требуется непрерывность поля напряжений в рассматриваемой области при переходе от элемента к элементу 175] . Тогда имеется возможность аппроксимировать поле напряжений в пределах одного элемента пробной функцией Р с варьируемыми коэффициентами р , которые не несут конкретного физического смысла. Поле перемещений представляется, как и раньше, пробными функциями А и узловыми перемещениями ^ . Вследствие, независимости параметров ft для различных элементов они могут быть выражены через узловые значения перемещений Ц- , что в результате ведет к уравнению метода перемещений (I). Однако,Вебеке показал [75] ,что в большинстве случаев такой подход ведет к тем же самым результатам, что и обычный метод перемещений, основанный на минимизации потенциальной энергии и имеющий более простой алгоритм.

Если аппроксимирующие функции для напряжений 6 удовлетворяют условиям равновесия, то имеем принцип минимума дополнительной энергии. В конечноэлементной формулировке в качестве неизвестных фигурируют узловые значения напряжений р. Такая модель конечного элемента (КЭ) называется равновесной и ведет к разрешающему уравнению метода сил [ 4 J .

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

Например, модифицированный принцип дополнительной энергии, базирующийся на функционале Жпс (см.рис.1), может быть получен путем включения усилий межэлементного взаимодействия, как ограничений, в функциал и соответствующих множителей Лагранжа, которые в данном случае будут иметь физический смысл перемещений границ элемента. Таким образом, этот принцип имеет в своей основе аппроксимацию двух полей: напряжений и перемещений [92- Ю1] . Граничные перемещения и без труда могут быть аппроксимированы через узловые значения перемещений и функций координат L . Поле напряжений аппроксимируется независимо в каждом элементе через параметры fi и функции координат Р . Применив принцип минимума дополнительной энергии на уровне элемента путем варьирования параметрами fi можно выразить последние через узловые перемещения ^ , которые и будут фигурировать в качестве неизвестных при рассмотрении конструкции в целом. Этот метод позволяет получить гибридную модель, основанную на аппроксимации поля напряжений или, просто, гибридную модель напряжений. Этот же модифицированный принцип может быть использован иначе, когда граничное взаимодействие выражается в интегральной форме. Однако и в этом случае узловыми неизвестными остаются перемещения ^ , но это будет уже равновесная модель Я 2 (в отличие от модели № I, описанной выше) [75],

Альтернативный подход возможен при использовании модифицированного принципа минимума потенциальной энергии. В этом случае с помощью параметров об , независимых от элемента к элементу, и функций координат С аппроксимируется поле перемещений в пределах элемента. В качестве множителей Лагранжа будут выступать граничные напряжения. Так как аппроксимация перемещений независима от элемента к элементу, параметры об могут быть выражены через узловые значения напряжений р и исключены на элементном уровне. Результирующие уравнения будут содержать только узловые значения напряжений в качестве неизвестных [7 г J . Это, очевидно, будут уравнения метода сил, а модель конечного элемента называется гибридная модель перемещений I.

Другая модификация принципа минимума потенциальной энергии Л~тРг основывается на независимой аппроксимации поля перемещений и в пределах конечного элемента, граничных перемещений й и введения межэлементного взаимодействия в качестве множителей Лагранжа [fOO]T. Так как Г и и независимы от элемента к элементу, то они могут быть выражены через узловые перемещения ц, , что ведет в итоге к уравнению метода перемещений. Эта модель конечного элемента носит название гибридная модель перемещений }Ь 2 (см.рис.1). Существует еще одна модификация функционала потенциальной энергии J/mp3 ,. которая основывается на независимой аппроксимации поля перемещений внутри элемента U и перемещений границ элемента и через узловые значения перемещений ^ Тогда результирующие уравнения будут уравнениями метода перемещений, а модель будет называться гибридной моделью перемещений В 3

Рассмотрим перечисленные вариационные принципы в приложении к решению задач изгиба пластин и оболочек. Все модели конечных элементов, основанные на гипотезе КирхгоффаДява и базирующиеся на принципе минимума потенциальной энергии (традиционном или модифицированном) требуют непрерывности прогибов и производных по нормали иг)п к границе между элементами. Построить совместную модель, т.е. удовлетворяющую условиям не-"прерыЕНОСти ш и иг)П , для традиционного принципа минимума потенциальной энергии удается только в случае включения в КЭ дополнительных узлов [70, 72] . В большинстве случаев для использования простейших конечных элементов пренебрегают требованием совместности и удовлетворяют лишь условиям непрерывности прогибов иг . В этом случае не требуется введения дополнительных узлов, однако применение подобных КЭ требует большой осторожности, так как во многих случаях результаты расходятся при измельчении сетки конечных элементов, вследствие нарушения условия совместности [ 17, 2B,3Sf Ьк^б, 95].

Простейший конечный элемент пластины Кирхгоффа-Лява, удовлетворяющий условиям совместности, можно построить используя гибридную модель напряжений. В работах[6Ъ69,0% т] представлены матрицы жесткости простейших конечных элементов, основанных на применении модифицированного принципа ffp. Однако количество коэффициентов jb в аппроксимирующих полиномах поля напряжений должно быть достаточно большим. В противном случае можно получить вырожденное решение (движение тела, как жесткого целого [91]). Однако при значительном количестве коэффициентов Jb нет возможности получить аналитическое выражение для матрицы жесткости, т.к. требуется обращение матрицы размером п*л , где л - количество коэффициентов Р (обычно от 18 до 24). Это приводит к необходимости обращения матрицы в процессе вычислений по каждому элементу, что существенно увеличивает время расчета.

Применение принципа представляет определенные трудности, так как требуется аппроксимировать три поля: поле перемещений Енутри элемента, поле граничных перемещений и поле межэлементного взаимодействия. Однако работы [69,100] показывают, что в большинстве случаев, даже если удается получить матрицу жесткости, основанную на этом гринципе, решение оказывается идентичным с обычной совместной моделью. Таким образом предпочтение следует отдать традиционному вариационному принципу ЗГр .

Простейшая гибридная модель перемещений использовалась как для получения матриц жесткости треугольных, так и четырехугольных элементов [&D,$i] . Однако при определенной форме треугольных элементов модель становится крайне нестабильной и дает расходящийся результат . Это является серьезным препятствием при применении указанных элементов.

Применение смешанной модели, основанной на принципе Рейсснера , позволяет достаточно просто построить аппроксимирующие полиномы для прогибов иг и напряжений б, т.к. для и/ требуется непрерывность только самой функции (но не ее производных). Для треугольных и четырехугольных элементов пластины характерна наилучшая точность и сходимость решения при определении поля напряжений по сравнению с другими элементами £6 5 7]. Однако в каждом узле такой конечноэле-ментной схемы содержится практически удвоенное количество неизвестных по сравнению с остальными моделями, что Еедет к существенному увеличению времени решения системы уравнений, увеличению настолько большому, что во многих случаях исключается возможность решения задачи.

Рассмотрим возможность построения модели КЭ пластины типа Тимошенко. В этом случае, применяя принцип минимума потенциальной энергии, можно построить совместную модель простейшего конечного элемента, т.к. независимая аппроксимация прогибов и углов поворота позволяет применить простейшие аппроксимирующие полиномы. Однако некоторые исследования показали, что такие КЭ дают удовлетворительную точность для пластин средней толщины, в случае достаточно тонких пластин модель дает значительную погрешность. Имеются рекомендации по улучшению сходимости путем сокращенного интегрирования £107], двойной аппроксимации угла поворота [6,9] , применения моментной схемы конечных элементов £52] , получившие обобщение в работе В.П.Болдычева LW] . Одним из вариантов улучшения сходимости является введение дополнительных узлов на сторонах конечного элемента [Ьв-^д], однако это увеличивает количество неизвестных в расчетной схеме.

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

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

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

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

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

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

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

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

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

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

Разработанная методика, алгоритмы и программы приняты к использованию в Управлении главного конструктора КамАЗа при создании системы автоматизированного проектирования и доводки (САЛРиД) автомобиля.

Работа прошла апробацию на научно-технических конференциях молодых специалистов КамАЗа (г.Брежнев, 1979,1980,1981г.г.), ЗИЛа (г.Москва, 1981,1983г.г.), научно-технической конференции Казанского инженерно-строительного института (г.Казань, 1982г.), семинаре по строительной механике под руководством профессора Л.А.Розина при Ленинградском политехническом институте (г.Ленинград, 1982г.), семинаре координационного Совета по САПР Министерства автомобильной промышленности СССР (г.Москва,1980г), республиканской научно-технической конференции "Механика сплошных сред" (г.Брежнев, 1982г.), Всесоюзной школе молодых ученых

Актуальные проблемы механики оболочек" (г.Казань, 1983г.), Всесоюзных школах-семинарах по методу конечных элементов (г.Рига, 1981г.; г.Киев, 1983г.).

Основное содержание работы изложено в восьми печатных работах [<2,29-33, £6, S~r] . В работе [t2] автору принадлежит разработка расчетных схем кронштейнов рамы автомобилей КамАЗ, анализ напряженно-деформированного состояния и поиск рациональной конструкции кронштейнов. В работе [3l] автором разработан программный комплекс и методика расчета для балки передней оси грузового автомобиля. В работе [32] автору принадлежит разработка идеологии и программного обеспечения для решения задач по МКЭ на мини-ЭВМ класса СМ-4. В работе [33] автором предложен и внедрен в практику расчетов способ построения матрицы жесткости оболочки произвольной геометрии на основе плоских (или пологих) конечных элементов. В работе [56] автором проведен вывод основных матричных соотношений для решения задачи аппроксимации и внедрение их в практику расчетов. В работе [5г] автору принадлежит выеод матрицы жесткости для треугольного и четырехугольного КЭ, их тестирование и внедрение в практику расчетов.

Заключение диссертации по теме «Механика деформируемого твердого тела», Малинин, Михаил Юрьевич

Основные результаты работы сводятся к следующему :

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

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

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

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

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

6. Решена задача построения итерационного процесса расчета дискретно-равнопрочной конструкции.

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

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

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

10. Разработанные модели и алгоритмы реализованы в комплексе программ на мини-ЭВМ " //OVA -1220" и внедрены для проектирования и доводки автомобилей на Камском комплексе по производству большегрузных автомобилей (КамАЗ).

3 А К Л 10 Ч Е Н И Е

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

Список литературы диссертационного исследования кандидат физико-математических наук Малинин, Михаил Юрьевич, 1984 год

1. Аньютц Ф., Утку 0. Алгоритм автоматической перенумерацииузлов для минимизации ширины ленты в матрице жесткости. -Ракетная техника и космонавтика, 1968, 4, с.212-214.

2. Алфутов Н.А., Попов Б.Г., Быков Е.В. Применение смешанныхфункционалов в численных методах расчета конструкций,- Известия вузов, М., Машиностроение, 1983, 9, с.3-7.

3. Амбарцумян С.А. Общая теория анизотропных оболочек. -М.:1. Наука, 1974.-446с.

4. Аргирис Д. Современные достижения в методах расчета конструкций с применением матриц. -М.: Стройиздат, 1968. -241с.

5. Ашкенази Е.К., Гамов Э.В. Анизотропия конструкционных материалов. -Л.: Машиностроение, 1980. 247с.

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

7. Блох В.И. Теория упругости. -Харьков: Издательство Харьковского государственного университета, 1964. -483с.

8. Болдычев В.П. Двойная аппроксимация угла поворота при расчете пластин средней толщины методом конечных элементов. -Известия ВНИИГ им.Б.Е.Веденеева. Сборник научных трудов, 1979, т.133, с.68-74.

9. Болдычев В.П. О решении статических и динамических задачизгиба пластин средней толщины. -Известия ВНИИГ им.Б.Е.Веденеева. Сборник научных трудов, 1978,т120

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

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

12. Валькова Л.А., Малинин М.Ю. Выбор рациональной конструкциикронштейнов крепления деталей для автомобилей КамАЗ. -Тезисы докладов X научно-технической конференции Управления главного конструктора КамАЗ.: Брежнев, 1983, с.13.

13. Власов В.З. Общая теория оболочек. -М, -Л.: ГТТИ, 1949.-784с.

14. Галимов К.З., Саченков А.В., Артюхин Ю.П., Карасев С.Н.,

15. Драган В.И., Костин В.А., Якушев Н.З.,Теория оболочек с учетом поперечного сдвига. -Казань: Издательство Казанского университета, 1977.- 212с.

16. Гордон JI.A., Розин Л.А. Метод конечных элементов в теориипластин и оболочек. -Известия ШИИТ им.Б. Е.Веденеева. Сборник научных трудов, 1971, т.95, с.85-97.

17. Городецкий А.С., Заворицкий В.И., Лантух-Лященко А.И.,

18. Рассказов А.О. Метод конечных элементов в проектировании транспортных сооружений. -М.: Транспорт, 1981. 143с.

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

20. Григоренко Я.М., Василенко А.Т. О некоторых подходах к построению уточненных моделей теории анизотропных . оболочек переменной толщины. -В сб.: Математические методы и физико-механические поля, Киев, 1978, 7, с.21-25.

21. Григоренко Я.М., Кокошин С.С. Численный анализ напряженного состояния слоистых анизотропных оболочек на базе смешанной модели ЖЭ. -Прикладная механика. К.: Наукова думка, 1982, т.18, 2, с.3-6.

22. Дорофеев С.М. Применение пологих элементов к расчету произвольных оболочек. -Рукопись депонирована в ВШШТИ, 1982, В 1849-82 Деп -8с.

23. Забияка Г.И., Прусаков А.П. О концентрации напряжений вмногослойных оболочках несимметричного строения. -Прикладная механика. К.: Наукова думка, 1974, т.1, 10, с.127-132.

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

25. Камель Х.А., Эйзенштейн Г.К. Автоматическое построениесетки в двух- и трехмерных составных областях. -В сб.: Расчет упругих конструкций с использованием ЭВМ. Л.: Судостроение, 1974, т.2, с. 74-89.

26. Колтунов М.А., Кравчук А.С., Майборода В.П. Прикладная механика деформируемого твердого тела. -М.: Высшая школа, 1983. -349с.

27. Корнеев В.Г. О методе конечных элементов для решения задачупругого равновесия. -В кн.: Строительная механика сооружений. Л.: Издательство ЛПИ, 1971, с.28-45.

28. Корнеев В.Г. Схемы метода конечных элементов высоких порядков точности. -Л.: Издательство Ленинградского университета, 1977. -208с.

29. Малинин М.Ю. Автоматизация расчетов на прочность несущихконструкций автомобиля. -В сб.: Тезисы докладов Республиканской научно-технической конференции "Механика сплошных сред". Набережные Чёлны, 1982, с.148.

30. Малинин М.Ю. Особенности реализации метода конечных элементов на мини-ЭВМ. -Тезисы доклада на семинаре-по строительной механике под руководством проф.Розина Л.А. при ЛПИ. В ж.: Строительная механика и расчет сооружений. М.: Стройиздат, 1983, I.

31. Малинин М.Ю., Воронов В.И., Козлов В.А. Повышение несущейспособности балки передней оси грузовых автомобилей. -Автомобильная промышленность. М.: Машиностроение, 1983, 9, с.16-17.

32. Малинин М.Ю., Козлов В.А. Организация программ МКЭ на мини

33. ЭВМ. -В сб.: Тезисы докладов ХХХ1У научной конференции Казанского инженерно-строительного института. Казань, 1982, с.10.

34. Малинин М.Ю., Снигирев В.Ф. К задаче формирования общейматрицы жесткости оболочки как совокупности плоских конечных элементов. -Депонированные рукописи, 1982, 12, б/о 287. Рук.деп. в ВИНИТИ, В. 3807-82 Деп .-5с.

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

36. Норри Д., де Фриз 1. Введение в метод конечных элементов.-М.: Мир, 1981. -304с.

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

38. Олсон М.Д. Исследование произвольных оболочек с помощьюпологих оболочечных конечных элементов. -В кн.: ! Тонкостенные оболочечные конструкции. М.: Машиностроение, 1980, с.409-437.

39. Паймушин В.Н. К проблеме расчета пластин и оболочек сосложным контуром. -Прикладная механика, 1980, т.16, 4, с.63-70.

40. Паймушин В.Н., Демидов В.Г. Об одном варианте соотношенийтеории среднего изгиба многослойных оболочек сложной геометрии. -В сб.: Статика и динамика оболочек. Казань, 1979, 12, с.53-60.

41. Пелех Б.Л. Теория оболочек с конечной сдвиговой жесткостью.-К.: Наукова думка, 1973. -248с.

42. Пелех Б.Л., Лазько В.А. Слоистые анизотропные пластины иоболочки с концентраторами напряжений. -К.: Наукова думка, 1982. -296с.

43. Пелех Б.Л., Марчук М.В. Метод конечных элементов при решении краевых задач для анизотропных пластин из композитных материалов. -Механика композитных материалов. Рига: Зинатне, 1981, 5, с.815-820.

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

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

46. Расчет упругих конструкций с использованием ЭВМ. -В кн.:

47. Сборник докладов международного симпозиума в Льеже в 1970г. Л.: Судостроение, 1974. -304с.

48. Рикардс Р.Б., Тетере Г.А. Устойчивость оболочек из композитных материалов. -Рига: Зинатне, 1974. -270с.

49. Рикардс Р.Б., Чате А.К. Изопараметрический треугольныйконечный элемент многослойной оболочки по сдвиговой модели Тимошенко. 4.2. Численные примеры. -Механика композитных материалов. Рига: Зинатне, 1981, 5, с.815-820.

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

51. Сакович А.И., Холмянский И.А. .Минимизация ширины лентысистемы уравнений в методе конечных элементов. -Проблемы прочности. К.: Наукова думка, 1981, 4, c.II7-II8.

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

53. Сегерлинд Л. Применение метода конечных элементов. -М.:1. Мир, 1979. -392с.

54. Сливкер В.И. Метод Ритца в задачах теории упругости, основанный на последовательной минимизации двух функционалов. -Известия АН СССР. Механика твердого тела, 1982, 2, с.57-64.

55. Снигирев В.Ф. К расчету слоистых панелей пола рефрижераторного вагона. -Депонированные рукописи, 1981, 12. Рук.деп. в ЦНТИИЭИТШШШ, В 784. -21с.

56. Снигирев В.Ф., Малинин М.Ю. Уточненный конечный элементслоистой пологой оболочки типа Тимошенко. -В сб.: Тезисы докладов Всесоюзной школы молодых ученых "Актуальные проблемы механики оболочек". Казань, 1983, с.199-200.

57. Тетере Г.А., Рикардс Р.В., Нарусберг В.Л. Оптимизация оболочек из слоистых композитов. -Рига: Зинатне,1978. -238с.

58. Тимошенко С.П., Войновский-Кригер С. Пластинки и оболочки.-М.: Физматгиз, 1963. -636с.

59. Тимошенко С.П., Гудьер Дж. Теория упругости. -М.: Наука,1979. -560с.

60. Фесик С.П. Справочник по сопротивлению материалов. -К.:

61. Буд I вельник, 1970. -308с.

62. Фокс А., Пратт М., Вычислительная геометрия. Применениев проектировании и на производстве. -М.: Мир, 1982. -304с.

63. Цурпал И.А., Тамуров Н.Г. Расчет многосвязных слоистыхи нелинейно-упругих пластин и оболочек. -К.: Вища школа, 1977. -22 4с.

64. Чулков П.П., Паймушин В.Н. К проблеме расчета многослойныхоболочек сложной геометрии. -В кн.: Статика сооружений. Киевский инженерно-строительный институт,1978, с.69-72.

65. Шабров Н.Н. Метод конечных элементов в расчетах деталейтепловых двигателей. -Л.: Машиностроение, 1983.- ,212с.

66. Шапошников Н.Н., Тарабасов Н.Д., Петров В.Б., Мяченков В.И.

67. Расчет машиностроительных конструкций на прочностьи жесткость. -М.: Машиностроение, 1981. -333с.

68. Allwood R.J., Cornes G.M.M. A polygonal finite element for plate bending problems using the assumed stress approach, Int. J.Num. Methods Eng., 1969, 1, p.p. 1Э$ Ш9.

69. Argyris J.H. Energy theorems and structural analysis, Butterworths Sci. Publications, London, I960.

70. Atlyri S., Pian T.H.H. Finite element analysis of shell of revolution by two doubly curved quadrilateral elements, J.Struct.Mech., 1973, 1, p.p. 393-U16

71. Bell K. Analysis of thin plates in bending using triangular finite element, Div. of Struct. Mech., Technical Univ. Norway, Trondheim, 1968.

72. Cavendish J.C. Automatic triangulation of arbitrary planar domains for the finite element method, Int. J. Num. Methods Eng., 197k, 8, p.p. 678 699

73. Clough R.'i., Tocher J.L. Finite element stiffness matrices for analysis of plate bending, Proc. Conf. Matrix Methods Struct. Mech., 1-st, AFFDL-TR-66-80, 1966, p.p. 5U6.

74. Fraens de Veubeke B. Displacement and equilibrium models in finite element method, in Stress Analysis, Zienkewicz O.C., Holister G.S.(eds.), London, 1965, p.p. lhS 197.

75. Imafuku I., Kodera Y., Sayawaki M. A generalized automatic mesh generation scheme for finite element method, Int. J. Hum. Methods Eng., 1980, 15, p.p. 713 731 .

76. Jones R.E. A generalization of the direct stiffness method of structural analysis, AIAA J., 196Ц, 2, p.p. 821 826.

77. Katoh H., Tahary H. Pre-processor for finite element analysis and its application to body structures, SAE Tech. Paper in Series 78063, 1973.

78. Kenneth H Huebner. The finite element method for engineers, N.J., 1975»

79. Kikuchi F., Ando Y. A new variational functional for the finite element method and its application to plate and shell problems, Nucl. Eng. Des., 1972, 21, p.p. 95 113.

80. Kikuchi F., Ando Y. Application of scimplified hybrid displacement method to large deflection analysis of elastic-plastic plates and shells, J.Fac. Eng., Univ. Tokyo, 1973, 32, p.p. 117 135.

81. Lynn P.P., Dhillon B.S. Triangular thick plate bending element. Proc. Int. Conf. Struct. Mech, in Reactor Technology, Berlin, 1971, Part M, NS, p.p. 1 - 27

82. Noor A.K., Andersen С.И. Mixed isoparametric finite element models of laminated composite shells, Сотр. Methods Appl. Mech. Eng., 1970, 11, p.p. 255 280.

83. Noor A.K., Hartley S.J. Nonlinear shell analysis via mixed isoparametric element, Сотр. Struct., 1977, 7, p.p. 6l5 626.

84. Noor A.K., Mathers M.D. Shear-flexible finite element models of laminated composite plates and shells, NASA, Т.Н., D 80liU, Dec. 1973.

85. Rorrie D.H., de ?ries G. The. finite element method fundamentals and applications, Academic Press, N.Y., 1973.

86. Pian Т.Н.Н. Element stiffness matrices for boundary compatibility and for prescribed boundary stress, Proc.Conf. Matrix Methods Struct. Mech., AFFDL-TR-66-80, 1966, p.p. U57-U77.

87. Pian T.H.H. Finite element methods by variational principles with relaxed continuity requirements, in variational methods in engineering, Brebbia G.A. and Tottenham (eds.) Southampton Univ. Press, 1973, 1, p.p. 3/1 3/2U

88. Pian T.H.H. Mau S.T. Recent study in assumed stress hybrid model, in Adv. Сотр. Methods Struct. Analysis and Designs, Oden J.T. (её.), . UAM Press, Univ. Alabama, Huntsville, 1972, p.p. 87 106.

89. Pian T.H.H., Tong P. Finite element methods in continuum mechanics, Adv. Appl. Mech. C.S. Yih (ed.), Academic Press, 1972, p.p. 1-58.

90. Przeminiecki J. Theoxy of matrix structural analysis, N.Y., 1969.

91. Reissner E. The effect of transverse shear deformation on the bending of plates, J. Appl. Mech., 19U5, 12, A-69-77.

92. Sovern R.T., Talor T. The finite element method for flexure of slabs when stress distribution are assumed, Proc. Inst. Civ. Eng., 1966, 3h , p.p. 153 170.

93. Steinmuller G. Restriction in the application of automatic mesh generation schemes by isoparametric co-ordinates, Int.J.Num. Methods Eng., 197U, 8, p.p. 289 29b.

94. Synge J.L. The hypercircle in mathematical Physics, Cambridge Univ. Press, 1957.

95. Tong P. New displacement hybrid finite element model for solid continua, Int. J. Num. Methods Eng., 1970, 2, p.p. 78 83.

96. Tong P., Pian T.H.H. A variational principle and the convergence of a finite element method based on assumed stress distribution, Int. J. Solids struct., 1969, 5, p.p. Ц63 H72.

97. Turner K.J., Clough R.W., Martin H.C., Topp L.J. Stiffness and deflection analysis of complex structures, J. Aero Sci., 1956, 23, p.p. 805 -823, p. 8SU.

98. Washizu K. Variational methods in elasticity and plasticity, 2-nd edn., Pergamon Press, Oxford, 1975.loU. V/hetstone W.D. SPAR Structural analysis system reference manual (systemlevel 11), NASA OR 1U5098-1, Feb. 1977.

99. Wolf J.P. Das Flaechentragwerksprogramm von STRIP, Schweiz Bauzeitung, 1972, Jan. 10, 90(3), p.p. Ul 52.

100. Wu SrC., Abel J.F. Representation and dicretization of arbitrary surfaces for finite element shell analysis, Int. J. Num. Methods Eng., 1979, 1U, p.p. 813 836.

101. Zienkievri.cz O.C., Hinton E. Reduced integration, function smoothing and non conformity in finite element analysis (with special reference to thick plates), J. Franklin Inst., 1976, 302, p.p. hk3 U6l.