Содержание к диссертации
Введение
Применение метода конечных элементов к расчету многослойных оболочек
1.1. Основные допущения, принятые при расчете многослойных пологих оболочек 20
1.2. Основные соотношения метода конечных элементов Z5
1.3. Использование конечных элементов пластины и пологой оболочки в расчетах оболочек сложной геометрии 32
1.4. Способы повышения точности расчетов методом конечных элементов для пластин и оболочек типа Тимошенко 36
2. Матрица жесткости треугольного конечного элемента слоистой пологой оболочки
2.1. Основные соотношения треугольного конечного элемента пологой оболочки с шестью узлами
2.2. Основные соотношения треугольного конечного элемента с учетом эффекта постоянства перерезывающих усилий
2.3. Исследование условий сходимости решения с использованием треугольных конечных элементов о1
2.4. Численное исследование точности расчета оболочек методом конечных элементов 65
3. Применение четырехугольных конечных элементов для расчета слоистой пологой оболочки 70
3.1. Матрица жесткости и вектор нагрузок четырехугольного элемента оболочки с учетом эффекта постоянства перерезывающих усилий 7Q
3.2. Исследование условий сходимости решения 2І
3.3. Тестовые задачи расчета оболочек с применением четырехугольных элементов *
4. Сглаживание разрывных функций методом конечных элементов
4.1. Применение метода конечных элементов к задаче аппроксимации напряжений
4.2. Аппроксимация с помощью треугольных элементов S8
4.3. Аппроксимация с помощью четырехугольных элементов 91
4.4. Исследование точности решения при использовании аппроксимирующих функций 94
4.5. Применение сглаживания в задачах расчета оболочек 96
5. Анализ напряженно-деформированного состояния пространственных тонкостенных конструкций "100
5.1. Расчет двухслойной пологой оболочки с отверстием
5.2. Автоматизация подготовки исходных данных и графического вывода результатов расчета пространственных конструкций 405*
5.3. Расчет тормозного барабана автомобиля КамАЗ-
5.4. Построение итерационного процесса определения толщин по элементам конструкции
5.5. Расчет несущей балки рамы автомобиля КамАЗ-55113 124
5.6. Расчет балки передней оси грузовых автомобилей во
5.7. Расчет кронштейна крепления поперечины В 3 автомобиля КамАЗ-5320
Заключение 9
Литература
- Основные допущения, принятые при расчете многослойных пологих оболочек
- Основные соотношения треугольного конечного элемента пологой оболочки с шестью узлами
- Матрица жесткости и вектор нагрузок четырехугольного элемента оболочки с учетом эффекта постоянства перерезывающих усилий
- Применение метода конечных элементов к задаче аппроксимации напряжений
- Расчет двухслойной пологой оболочки с отверстием
Введение к работе
В современном машиностроении все большее применение находят пространственные тонкостенные конструкции,?которые позволяют в значительной мере повысить прочность различных деталей и узлов при минимальных расходах материала. Наличие всевозможных концентраторов и подкрепляющих элементов требует проведения комплекса исследований напряженно-деформированного состояния при создании подобных конструкций.
Увеличение нагрузок и задача снижения материалоемкости конструкций предъявляют высокие требования к их прочностным качествам и ведут к существенному усложнению геометрии.
Реальные машиностроительные конструкции состоят из набора пластин, цилиндрических, конических и пологих оболочечных участков, имеющих всевозможные вырезы и отверстия. Многие детали представляют собой двухслойные и многослойные пластины и оболочки.
Решить поставленные задачи проектирования в большой мере позволяет развитие современных мєтодое расчета тонкостенных конструкций.
Теоретические основы расчета многослойных конструкций заложены е работах Амбарцумяна С.А., Александрова А.Я., Болотина В.В., Григолюка Э.И., Григоренко Я.ЇЛ., Галимова К.З., Королева B.PI..,. Куршина Л.М., Муштари Х.М., Новичкова Ю.М., Пелеха Б.Л., Прусакова А.П., Саченкова А.В., Тетерса Г.А., Чулкова А.П. и др. При этом весьма эффективным в практических расчетах показал себя метод гипотез, когда некоторые допущения о напряженно-деформированном состоянии по толщине оболочіш или пластины позволяют свести трехмерную задачу к двухмерной.
Можно выделить два направления в применении метода гипотез. Первое направление основано на применении кинематических гипотез для каждого отдельного слоя, когда порядок системы уравнений зависит от числа слоев. Значительное развитие указанное направление в теории слоистых сред получило в работах Болотина В.В. и Новичкова Ю.Н. В работе [ -и j на основе единого подхода предложена теория многослойных конструкций при произвольном числе слоев и широких допущениях о свойствах отдельных слоев. Сравнительный анализ и развитие теорий для многослойных оболочек переменной толщины, которым соответствуют различные предположения об изменении сдвиговых деформаций и напряжений, даны в работах Григоренко Я.М. [-18,19] Л Подобный подход позволяет получить решения для задач со значительным изменением параметров упругости по слоям.
В тех случаях, когда свойства материалов по слоям разнятся незначительно, весьма эффективным оказывается второе направление метода гипотез, когда определенные кинематические допущения принимаются для всего пакета слоев в целом. В этом случае порядок разрешающих уравнений не зависит от числа слоев. Широкое развитие в этом направлении получила модель С.П.Тимошенко. Для пластин неклассическая теория была построена Э.Рейсснером [96] и Р.Миндлиным l&hj . Обобщение модели и ее развитие для оболочек можно найти в работах Галимова К.З., Саченкова А.В. и др. [14] , Пелеха Б.Л. [М,42], Цурпала И.А., Тамурова Н.Г. [6 3] . Развитие модели С.П.Тимошенко в применении к расчету устойчивости и оптимизации многослойных оболочек из композитных материалов приведено в работах Тетерса Г.А.,
Рикардса Р.Б. [Ь?,8] . Расчет реальных многослойных конструкций сложной геометрии позволяют сделать теоретические разработки, предложенные в работах Чулкога П.П. и Паймушина В.Н. [39,
Однако даже при наличии разработанного математического аппарата исследования многослойных конструкций для эффективного решения задачи рационального проектирования машиностроительных конструкций при всем их многообразии необходимо применение численных методов расчета. Метод конечных элементов (МКЭ) является в этом отношении наиболее универсальным, так как позволяет рассчитывать конструкции произвольной геометрии при сложном характере нагружения. Основное достоинство МКЭ заключается в том, что он объединяет преимущества численных методов решения Еариационных задач (Ритца, БубноЕа-Галеркина и т.д.) и сеточных методов.
За сорок лет (начиная с работ Сига.л R.C?3]^ Ръадег V/ЛЩ
Тиьпег М. CtozJ ) МКЭ получил широкое развитие как за рубежом, так и в отечественной литературе в работах Розина Л.А., Постнова В.А., Шапошникова Н.Н., Григоренко Я.М., Корнишина М.С., Еіурмана З.И., Сахарова А.С. и др. Вариационные принципы, на которых базируется метод, могут быть сведены в общую схему, приведенную на рис.І. В верхней части схемы представлены традиционные вариационные принципы.
Принцип минимума потенциальной энергии ' (функционал ХР ) содержит в качестве искомой функции поле перемещений U Перемещения и аппроксимируются пробными функциями А , - 7 —
Традиционные вариационные принципы
Принцип Ху-Вашицу
Принцип Рейсснера
Принцип Кастильяно forcC Совместная модель I Смешанная модель 2 Смешанная модель Равновесная модель I ТГс(Р) Модифщированньге вариационные принципы Принципы Лагранжа а=сы Гибридная модель перемещений I Смешанный метод ЧГ-Н6ЇР] Го Л^(«д,">^ АКргС"*")^ T~A/f Гибридная модель перемещений 2 Гибридная модель .перемещений 3 Метод сил С/р-Р Принціш Кастильяно Гибридная модель напряжений г^Ц I СГ= Pfi | [Равновесная —модель 2 *н*.(АЛ)-*-7іп«Ф Рис.1 Диаграмма вариационных принципов метода конечных элементов — 8 — где в качестве варьируемых коэффициентов фигурируют узловые обобщенные перемещения у . Окончательное расчетное уравнение можно представить в форме : K-y=Q, (і) где К - матрица жесткости конструкции U - Еектор обобщенных узловых сил. К уравнению (I) приводит известный метод перемещений [ 16, 2 6,2 7, 44, 45,50]. Наиболее общий принцип Ху-Вашицу, основанный на варьировании перемещениями, деформациями и напряжениями, позволяет получить принцип Рейсснера (функционал ^ ), базирующийся на вариации перемещений и напряжений [50,68] .В этом случае конечноэлементная формулировка предполагает аппроксимацию пробными функциями как поля перемещений, так и поля напряжений. В качестве искомых величин выступают узловые значения перемещений Ц. и узловые значения напряжений р . Окончательное уравнение смешанного метода, который вытекает из функционала 7fR , имеет вид : -Н &1 fal (0] (2) где И, G, Q - матрицы, полученные путем суммирования соответствующих матриц конечных элементов. Глобальная матрица в уравнении (2) не является положительно определенной и требует специализированных программ для решения системы уравнений. Однако принцип Рейсснера можно представить в форме, когда не требуется непрерывность поля напряжений в рассматриваемой области при переходе от элемента к элементу [ 75] . Тогда имеется возможность аппроксимировать поле напряжений в пределах одного элемента пробной функцией Р с варьируемыми коэффициентами р , которые не несут конкретного физического смысла. Поле перемещений представляется, как и раньше, пробными функциями А и узловыми перемещениями cj, . Вследствие, независимости параметров J5 для различных элементов они могут быть выражены через узловые значения перемещений Ц- , что в результате ведет к уравнению метода перемещений (I). Однако,Вебеке показал [75] ,что в большинстве случаев такой подход ведет к тем же самым результатам, что и обычный метод перемещений, основанный на минимизации потенциальной энергии и имеющий более простой алгоритм. Если аппроксимирующие функции для напряжений 6 удовлетворяют условиям равновесия, то имеем принцип минимума дополнительной энергии. В конечноэлементной формулировке в качестве неизвестных фигурируют узловые значения напряжений р. Такая модель конечного элемента (КЭ) называется равновесной и ведет к разрешающему уравнению метода сил [b] . При применении вышеизложенных традиционных принципов предполагалось, что искомые поля перемещений и напряжений аппроксимировались непрерывными функциями в пределах рассматриваемого объема деформируемого тела. Однако эти условия можно смягчить и применять разрывную аппроксимацию искомых функции, БКЛЮчрів при этом в минимизируемый функционал условие непрерыв- --/0- ности межэлементного взаимодействия, как ограничение, с соответствующим неопределенным множителем Лагранжа. Такой подход (модификация вариационных принципов) носит название гибридной формулировки.' Очевидно, что в этом случае параметры, описывающие искомую функцию, независимы от элемента к элементу и могут быть выражены через параметры, описывающие межэлементное взаимодействие путем применения вариационного принципа на урогне одного конечного элемента. Например, модифицированный принцип дополнительной энергии, базирующийся на функционале Жпе (см.рис.1), может быть получен путем включения усилий межэлементного взаимодействия, как ограничений, в функциал и соответствующих множителей Лагранжа, которые в данном случае будут иметь физический смысл перемещений границ элемента. Таким образом, этот принцип имеет в своей основе аппроксимацию двух полей: напряжений и перемещений [92- - 101] . Граничные перемещения и без труда могут быть аппроксимированы через узловые значения перемещений Ц> и функций координат L . Поле напряжений аппроксимируется независимо в каждом элементе через параметры fi и функции координат Р . Применив принцип минимума дополнительной энергии на уровне элемента путем варьирования параметрами fi можно выразить последние через узловые перемещения Cj, , которые и будут фигурировать в качестве неизвестных при рассмотрении конструкции в целом. Этот метод позволяет получить гибридную модель, основанную на аппроксимации поля напряжений или, просто, гибридную модель напряжений. Этот же модифицированный принцип может быть использован иначе, когда граничное взаимо- -н- дейстЕие выражается в интегральной форме. Однако и в этом случае узловыми неизвестными остаются перемещения О- , но это будет уже равновесная модель & 2 (в отличие от модели В I, описанной выше) [75]. Альтернативный подход возможен при использовании модифицированного принципа минимума потенциальной энергии. В этом случае с помощью параметров об , независимых от элемента к элементу, и функций координат С аппроксимируется поле перемещений в пределах элемента. В качестве множителей Лагранжа будут выступать граничные напряжения. Так как аппроксимация перемещений независима от элемента к элементу, параметры об могут быть выражены через узловые значения напряжений р и исключены на элементном уровне. Результирующие уравнения будут содержать только узловые значения напряжений в качестве неизвестных [7 г J . Это, очевидно, будут уравнения метода сил, а модель конечного элемента называется гибридная модель перемещений \1 I. Другая модификация принципа минимума потенциальной энергии Л~тРг основывается на независимой аппроксимации поля перемещений и в пределах конечного элемента, граничных перемещений и и введения межэлементного взаимодействия в качестве множителей Лагранжа [ЮО}1. Так как Т ш U независимы от элемента к элементу, то они могут быть выражены через узловые перемещения Cj, , что ведет в итоге к уравнению метода перемещений. Эта модель конечного элемента носит название гибридная модель перемещений }Ь 2 (см.рис.1). Существует еще одна модификация функционала потенциальной энергии ЗГтрг ,. которая основывается на независимой аппроксимации поля перемещений внутри элемента U и перемещений границ элемента и через узловые значения перемещений а Тогда результирующие уравнения будут уравнениями метода перемещений, а модель будет называться гибридной моделью перемещений В 3 Рассмотрим перечисленные вариационные принципы в приложении к решению задач изгиба пластин и оболочек. Все модели конечных элементов, основанные на гипотезе Кирхгоффа/Лява и базирующиеся на принципе минимума потенциальной энергии (традиционном или модифицированном) требуют непрерывности прогибов и производных по нормали иг)п к границе между элементами. Построить совместную модель, т.е. удовлетворяющую условиям непрерывности иг и иг)П , для традиционного принципа минимума потенциальной энергии удается только в случае включения в КЭ дополнительных узлов [70, 72] . В большинстве случаев для использования простейших конечных элементов пренебрегают требованием совместности и удовлетворяют лишь условиям непрерывности прогибов иг . В этом случае не требуется введения дополнительных узлов, однако применение подобных КЭ требует большой осторожности, так как во многих случаях результаты расходятся при измельчении сетки конечных элементов, вследствие нарушения условия совместности [17,2ЬІЗВ) ^^}бб, 95]. Простейший конечный элемент пластины Кирхгоффа-Лява, удовлетворяющий условиям совместности, можно построить используя гибридную модель напряжений. В работах 6^ ^^ 10*1 представлены матрицы жесткости простейших конечных элементов, — р. Однако количество коэффициентов & в аппроксимирующих полиномах поля напряжений должно быть достаточно большим. В противном случае можно получить вырожденное решение (движение тела, как жесткого целого [91]), Однако при значительном количестве коэффициентов р нет возможности получить аналитическое выражение для матрицы жесткости, т.к. требуется обращение матрицы размером п*п , где л - количество коэффициентов Р (обычно от 18 до 24). Это приводит к необходимости обращения матрицы в процессе вычислений по каждому элементу, что существенно увеличивает время расчета. Применение принципа ^»Рг представляет определенные трудности, так как требуется аппроксимировать три поля: поле перемещений Енутри элемента, поле граничных перемещений и поле межэлементного взаимодействия. Однако работы [69,100} показывают, что в большинстве случаев, даже если удается получить матрицу жесткости, основанную на этом іринципе, решение оказывается идентичным с обычной совместной моделью. Таким образом предпочтение следует отдать традиционному вариационному принципу VTp . Простейшая гибридная модель перемещений ^рз использовалась как для получения матриц жесткости треугольных, так и четырехугольных элементов [&D,$i] . Однако при определенной форме треугольных элементов модель становится крайне нестабильной и дает расходящийся результат /"3 3J . Это является серьезным препятствием при применении указанных элементов. Применение смешанной модели, основанной на принципе Рейсснера 3^ , позволяет достаточно просто построить ап- проксимирующие полиномы для прогибов их и напряжений б, т.к. для и/ требуется непрерывность только самой функции (но не ее производных). Для треугольных и четырехугольных элементов пластины характерна наилучшая точность и сходимость решения при определении поля напряжений по сравнению с другими элементами 35-аг]. Однако в каждом узле такой конечноэле-ментной схемы содержится практически удвоенное количество неизвестных по сравнению с остальными моделями, что еєдєт к существенному увеличению времени решения системы уравнений, увеличению настолько большому, что во многих случаях исключается возможность решения задачи. Рассмотрим возможность построения модели КЭ пластины типа Тимошенко. В этом случае, применяя принцип минимума потенциальной энергии, можно построить совместную модель простейшего конечного элемента, т.к.. независимая аппроксимация прогибов и углов поворота позволяет применить простейшие аппроксимирующие полиномы. Однако некоторые исследования показали, что такие КЭ дают удовлетворительную точность для пластин средней толщины, в случае достаточно тонких пластин модель дает значительную погрешность. Имеются рекомендации по улучшению сходимости путем сокращенного интегрирования [10г], двойной аппроксимации угла поворота [6,9] , применения моментной схемы конечных элементов [52] , получившие обобщение е работе В.П.Болдычева 1ю] . Одним из вариантов улучшения сходимости является введение дополнительных узлов на сторонах конечного элемента Ц&~4У], однако это увеличивает количество неизвестных в расчетной схеме. В связи с этим представляет интерес разработка конечного — 45 — элемента на основе модели СП.Тимошенко, имеющего минимальное количество узлов, обладающего высокой сходимостью и позволяющего выполнять расчеты многослойных пластин и оболочек в широком диапазоне толщин (разработки по КЭ многослойных пластин и оболочек на основе других моделей можно найти в работах [гяЛз]. При разработке КЭ требуется уделить особое внимание вопросам сглаживания разрывного поля напряжений, получаемого при расчетах по принципу Лагранжа, автоматизации подготовки исходных данных и обработки результатов расчета при проектировании равнопрочных конструкций. Настоящая работа посвящена исследованию напряженно-деформированного состояния однослойных и многослойных оболочек сложной геометрии методом конечных элементов. На защиту еыносятся следующие положения: а) получение матриц жесткости конечных элементов много слойной пологой оболочки, обладающих повышенной сходимостью решения при минимальном количестве узлов; б) метод аппроксимации разрывной функции непрерывной на основе решения вариационного уравнения методом конечных элемен тов,* в) алгоритм автоматизированного построения сетки конечных элементов для пространственных тонкостенных конструкций и гра фического отображения результатов расчета на периферийных уст ройствах; г) методика построения дискретно-равнопрочной конструкции путем варьирования толщиной конечных элементов в расчетной схеме; д) исследование напряженно-деформированного состояния тонкостенных узлов автомобиля КамАЗ. Теоретические предпосылки, методы и результаты решения поставленных задач изложены в пяти главах диссертации. В первой главе приведены основные допущения, принятые при расчетах многослойных пологих оболочек. Рассмотрены пути расчета тонкостенных конструкций методом конечных элементов, предложен способ построения расчетной схемы оболочки произвольной геометрии на основе плоских или пологих конечных элементов. Исследованы известные способы построения аппроксимирующих функций для компонент перемещений с целью повышения точности расчетов для пластин и оболочек типа Тимошенко. Во второй главе получены основные соотношения для треугольного шестиузлового и трехузлового элементов многослойной пологой оболочки несимметричного строения. При Еыводе соотношений шестиузлового элемента применялась квадратичная аппроксимация поперечного прогиба иґ путем введения соответствующих степеней свободы в узлах, лежащих на серединах сторон КЭ. Для аппроксимации тангенциальных смещений и и if применялся линейный полином. Квадратичная аппроксимация поперечного прогиба для трехузлового элемента строилась из условия постоянства перерезывающих усилий вдоль сторон элемента, что позволяло определить дополнительные коэффициенты в аппроксимации иґ через коэффициенты линейных полиномов смещений и и if . Исследована сходимость результатов расчета на основе разработанных КЭ при решении тестовых задач. В третьей главе приведен вывод основных соотношений для произвольного четырехугольного элемента многослойной пологой оболочки. Аппроксимация поперечного прогиба строилась на основе гипотезы постоянства перерезывающих усилий по граням КЭ. Это позволило повысить степень аппроксимирующего полинома для иґ , причем дополнительные коэффициенты, так же как и для трехузло-вого элемента, определялись через коэффициенты билинейной аппроксимации смещений и и if . Все зависимости строились в локальной системе координат с последующим отображением в глобальную систему с помощью Якобиана преобразования координат. В четвертой главе разработана методика аппроксимации разрывного поля напряжений, полученного при решении задач методом перемещений. Методика основана на решении вариационного уравнения задачи аппроксимации методом конечных элементов. Для решения задачи аппроксимации использовалась та же сетка КЭ, которая применялась при решении основной задачи определения напряженно-деформированного состояния оболочки. Полученные матрицы представлены в яеном виде и позволяют распространить методику для определения значений любой разрывной функции в узлах конечно-элементной сетки. Непрерывная аппроксимирующая функция позволяет автоматизировать процесс обработки результатов расчета по В пятой главе произведен расчет реальных тонкостенных конструкций, как однослойных, так и многослойных. На примере расчета пологой сферической оболочки показана эффективность МКЭ. Для облегчения этапа подготовки данных к численному расчету и существенного ускорения процесса определения рациональной геометрии конструкции разработаны алгоритмы автоматизированной генерации конечноэлементной сетки, графического отображения ре- -48 — зультатоЕ расчета и определения толщин,по элементам, с целью удовлетворения условиям равнопрочности. На основе разработанных алгоритмов произведен расчет узлоЕ автомобилей КамАЗ, даны рекомендации по их рациональному проектированию, что позволило получить значительную экономию металла. Результаты выполненной работы могут быть использованы в машиностроении при проектировании узлов и деталей средней толщины, в судостроении при расчетах корпусных конструкций, в строительстве при анализе напряженно-деформированного состояния несущих панелей и куполоЕ, в автомобилестроении при проектировании и доводке рам, состоящих из тонкостенных элементов, кронштейнов сложной геометрии с отверстиями и вырезами; кузовных конструкций как металлических, так и выполненных из многослойных композиционных материалов. Разработанная методика, алгоритмы и программы приняты к использованию в Управлении главного конструктора КамАЗа при создании системы автоматизированного проектирования и доводки (САПРиД) автомобиля. Работа прошла апробацию на научно-технических конференциях молодых специалистов КамАЗа (г.Брежнев, І979,І980,І98Іг.г.), ЗИЛа (г.Москва, 1981,1983г.г.), научно-технической конференции Казанского инженерно-строительного института (г.Казань, 1982г.), семинаре по строительной механике под руководством профессора Л.А.Розина при Ленинградском политехническом институте (г.Ленинград, 1982г.), семинаре координационного Совета по САПР Министерства автомобильной промышленности СССР (г.Москва, 1980г), республиканской научно-технической конференции "Механика сплошных сред" (г.Брежнев, 1982г.), Всесоюзной школе молодых ученых — A9 — "Актуальные проблемы механики оболочек" (г.Казань, 1983г.), Всесоюзных школах-семинарах по методу конечных элементов (г.Рига, 1981г.; г.Киев, 1983г.). Основное содержание работы изложено в восьми печатных работах [f2,29 -33, 56, 5Г] . В работе [t2] автору принадлежит разработка расчетных схем кронштейнов рамы автомобилей КамАЗ, анализ напряженно-деформированного состояния и поиск рациональной конструкции кронштейнов. В работе [зі] автором разработан программный комплекс и методика расчета для балки передней оси грузового автомобиля. В работе [32] автору принадлежит разработка идеологии и программного обеспечения для решения задач по МКЭ на мини-ЭВМ класса СМ-4. В работе [33] автором предложен и внедрен в практику расчетов способ построения матрицы жесткости оболочки произвольной геометрии на осноеє плоских (или пологих) конечных элементов. В работе [S6J автором проведен вывод основных матричных соотношений для решения задачи аппроксимации и внедрение их в практику расчетов. В работе [572 автору принадлежит выеод матрицы жесткости для треугольного и четырехугольного КЭ, их тестирование и внедрение в практику расчетов. — 20 — Г І А В A I Применение метода конечных элементов к расчету многослойных оболочек Рассмотрим многослойную пологую оболочку средней толщины, собранную из П упругих трансверсально анизотропных слоев. Слои имеют различные толщины и различные упругие характеристики. Координатная поверхность эквидистантна внешним поверхностям оболочки. Положение произвольной точки К определяется ортогональными координатами /, . По толщине оболочки положение точки определяется координатой % , отсчитываемой от поверхности о Для упрощения решения задачи об изгибе оболочки необходимо ввести допущения относительно напряженно-деформированного состояния, позволяющие сделать переход от трехмерной задачи к двухмерной. Согласно сдвиговой модели оболочки типа Тимошенко принимаем, что прямолинейный элемент нормали к срединной поверхности всего пакета слоев остается прямолинейным, не изменяет своей длины и поворачивается на некоторый угол относительно поверхности S при деформации оболочки. Рассмотрим вывод матрицы жесткости элемента, не являющегося простейшим. Однако такой элемент имеет достаточно высокую точность и позволит провести ряд сравнительных расчетов, сопоставить трудоемкость и точность вычислений при применении других элементов. Форму элемента и порядок аппроксимации принимаем согласно рис.1.5 исходя из соображений, изложенных в п.1.4. Тогда аппроксимирующие функции для перемещени. Использование треугольных КЭ позволяет решить весьма широкий круг задач"при сложной геометрии конструкции. Однако достаточно часто в машиностроительных узлах и деталях встречаются топологически регулярные участки четырехугольной (или прямоугольной) формы. Разбиение таких участков на треугольные I приводит к завышенному количеству требуемых элементов и узлов. Гораздо эффективнее, в подобных случаях применять четырехугольные КЭ, образуемые пересечением двух семейств линий (рис.1), т.к. при том же количестве узлов понадобится меньшее количество КЭ, а при одинаковой точности - меньшее количество и узлов, и элементов по сравнению с треугольной сеткой. Это объясняется наличием в аппроксимирующих функциях четырехугольного элемента ела-, гаемых более высокого порядка, чем у треугольного, т.е. сходимость имеет порядок на единицу выше. Рассмотрим четырехугольный КЭ с минимальным количеством узлов : четыре узла в вершинах элемента (рис.3.1). В каждом узле имеется пять степеней свободы. По аналогии с треугольным элементом следует ожидать, что эффект постоянства перерезывающих сил позволит повысить степень полішомиальной аппроксимации прогиба it и, тем самым, избежать завышения сдвиговой жесткости. Принадлежность пространству Lz означает, что напряжения могут претерпевать разрыв при переходе через границу элемента. Построение непрерывной функции напряжений позволяет во многом облегчить обработку результатов расчета и построение графиков, эпюр, полей напряжений и деформаций. Как правило, непрерывная аппроксимация позволяет уточнить полученное разрывное решение для напряжений. В настоящее время известны несколько способов аппроксимации напряжений. Все они базируются на вычислении значений напряжений в узлах конечноэлеглентнои сетки и аппроксимации их в пределах элемента с помощью функций формы, обеспечивающих непрерывность при переходе от элемента к элементу. В работе /"2 3/ напряжения в каждом узле вычисляются путем их осреднения по значениям в элементах, окружающих узел. Этот метод позволяет получить приемлемые результаты только во внутренних узлах сетки при незначительных градиентах напряжений. Предварительное вычисление напряжений в узлах каждого конечного элемента с последующим осреднением значений в узлепредлагается в работе [52J . Очевидно, подобный подход не на много отличается от первого, а в случае элемента с постоянным значением напряжений по области элемента оба подхода дают идентичный результат, Решение задачи расчета оболочки с отверстием проводилось для определения точности вычислений с использованием пологих конечных элементов, предложенных в п.2.2 и 3.1, и оценки эффективности применения ЫКЭ в задачах с большими несимметричными отверстиями, когда аналитические методы не позволяют получить решение. Рассмотрим оболочку, очерченную по части сферическом поверхности и имеющую подъем, не превышающий 1/5 диаметра основания оболочки, которую по В.З.Власову \13\ можно считать пологой. Оболочка имеет два слоя различной толщины, ослаблена центральным отверстием радиуса Ч-о , защемлена по периметру и нагружена равномерным внутренним давлением р . Решение для такой задачи известно [2"Й , причем так же, как и в [22], будем считать, что отверстие закрыто крышкой, которая передает только поперечные нагрузки по краю отверстия. Расчеты будем производить для безразмерной величины коэффициента концентрации Кв, определяемой соотношением.Основные допущения, принятые при расчете многослойных пологих оболочек
Основные соотношения треугольного конечного элемента пологой оболочки с шестью узлами
Матрица жесткости и вектор нагрузок четырехугольного элемента оболочки с учетом эффекта постоянства перерезывающих усилий
Применение метода конечных элементов к задаче аппроксимации напряжений
Расчет двухслойной пологой оболочки с отверстием
Похожие диссертации на Исследование напряженно-деформированного состояния пространственных тонкостенных конструкций сложной геометрии методом конечных элементов
