Введение к работе
Актуальность проблемы. Условия эксплуатации современных машин, аппаратов и сооружений характеризуются высоким уровнем нагрузок, температур и других типов физических полей, действув-щих в широкой диапазоне скоростей их изменения, сложным характером взаимодействия этих полей во времени и пространстве. Среди всего многообразия такого типа воэдействий, для практических целей оценки прочности конструкций, особый интерес представляют режимы так называемых квазистатических нагруяений, при которых времена приложения и действия нагрузок, температур и перемещений соизмеримы и больше характерных времен, определяемых минимальной частотой собственных колебаний конструкции. К такому классу нагружензй относятся как медленно меняющиеся и статически прикладываемые воздействия, так и динамические-воздействия, при которых волновые аффекты не играют судаствеиного влияния на процесс деформирования.
Постоянно растущие требования к снижению веса и металлоёмкости современных конструкций при увеличении их нагруженности допускают возможность работы отдельных их элементов за пределом упругости и даіе в условиях локального нарушения прочности, если последние не приводят к выходу из строя конструкции в целом. В связи с этим возникает естественная необходимость исследования поведения конструкций с учетом явлений пластичности, ползучести и развивающейся повреаденности материала, приводящей к разрушению конструкции. Теоретическое исследование поведения конструкций с учетом названных эффектов необратимого деформирования и разрушения представляет собой сложную проблему, успешное решение которой возможно лишь при условии решения комплекса самостоятельных- вопросов различного масштаба, и их число входят, вопросы создания моделей материала, описывающих процессы развития пластических деформация, ползучести и поврелдаяностн в точке тела ігри заданных законах изменения температуры, напряжений или деформаций, построение математической модели исследуемого объекта, описывающей поведение конструкции с учетом выбранных моделей материала, развитием методов и алгоритмов решения сложных нелинейных краевых задач к которым, в конечном счете, сводится исследование поведения конструкций. Каждый из названных вопросов представляет собой самостоятельную научную
_ 4 -
проблему, однако, практическое использование результатов решения каждой из них для оценки прочности конструкций возможно только при создании достаточно эффективной вычислительной модели, связывающей воедино все, используемые для анализа конструкций, частные модели, схемы и алгоритмы и обеспечивающей их оптимальное сочетание в составе единого вычислительного процесса.
Цель работы состоит в разработке комплексной методики и единого вычислительного алгоритма исследования процессов необратимого деформирования и разрушения конструкций, состоящих из элементов пластин, плит, оболочек различной толщины, а также массивных конструктивных образований при квазистатических силовых и температурных нагрукениях на основе современных моделей, описывающих поведение конструкционных материалов и численных методов решения нелинейных задач механики деформируемого твердого тела.
Научная новизна. В диссертации разработан единый подход к численному решении задач исследования процессов деформирования и разрушения оболочечшх конструкций при квазистатических силовых и температурных воздействиях и создан соответствупщий комплекс методических и программно-алгоритмических средств.
На основе изучения , анализа и теоретического обобщения работ отечественных и зарубежных ученых, предложен вариант модели поврежденного материала, описывающей процессы деформирования и разрушения начально изотропных конструкционных материалов с учетом взаимного влияния эффектов пластичности, ползучести и накопления повреждений. Модель устанавливает состав, структуру и функциональные связи между переменными, определяющими развитие названных эффектов как формально независимых процессов, и позволяет использовать для описания последних различные варианты конкретных математических моделей пластичности, ползучести и накопления повреждений. Получена система уравнений, описывающих процессы деформирования и разрушения конструкций, состоящих из элементов пластин, плит, оболочек различной толщины, а также фрагментов массивных образований при квазистатических силовых и температурных воздействиях с учетом эффектов пластичности, ползучести, развивающейся повреждвнности в геометрически линейной и нелинейной постановках. Разработана схема решения сформулированной нелинейной краевой задачи, заключающаяся в пошаговом интегрировании определяющих уравнений с использова-
ниєм предложенных в работе двухуровневой шаговой схемы и схем экстраполяции решения, применяемых для ускорения сходимости итерационных процессов, анализа поведения решения и его продолжения по параметрам нагружешзя.
Дано развитие и обоснование эффективности применения моделей МКЭ для анализа плит, пластан и оболочек различной толщины с высокой точностью и малой трудоёмкостью. Предложена комбинированная схема МКЭ, позволяющая с единых позиций исследовать НДС деформируемых систем, состоящих из тонкостенных и массивных образований на основе соотношений, наиболее соответствующих реализуемым в каждом из них видом НДС.
На основе предложенных моделей и численных схем разработана комплексная методика и эффективные алгоритмы численного исследования нелинейных процессов деформирования и разрушения оболочечкых конструкций, бешены новые задачи исследования процессов деформирования и разрушения конструкций с учетом эффектов пластичности, ползучести, накопления повреждений, геометрической нелинейности, а также связности процессов названных эффектов.
Достоверность полученных результатов подтверждается математическим обоснованием ряда принимаемых положений, решением большого числа тестовых задач с использованием различных альтернативных методик и программных средств,' сравнением получаемых решений с известными теоретическими и экспериментальными результатами, проведением исследования сходимости решений путем последовательного увеличения числа конечных элементов, шагов интегрирования уравнений или точности удовлетворения определяющих уравнений.
Практическую ценность работы составляют комплекс методических и программно-алгоритмических.средств для исследования процессов деформирования и разрушения оболочечкых конструкций на основе МКЭ, реализованных на различных типах ЭВМ, а также, имеющие самостоятельное значение отдельные его фрагменты, касавшиеся реализации моделей поведения материалов и моделей МКЭ, схем численного решения нелинейных задач, средств информационного обеспечения. Применение разработанных средств для моделирования процессов образования остаточных напряжений ш деформаций при сварке. Результаты численного исследования погрешности применения соотношений теории тонких оболочек для
некоторых типовых случаев локального распределения напряжений и деформаций, проведенные о целью получения практических рекомендаций при составлении расчетных схем конструкций. Результаты численного решения ряда прикладных задач. Перечисленные методические, и программные средства,а также некоторые из полученных в работе численных результатов внедрены в расчетную практику ряда научных и проектно-вонструкторских организаций.
На защиту выносятся следующие основные результаты диссертации :
вариант модели поврежденного материала, описывающей процессы деформирования и разрушения начально изотропных конструкционных материалов с учетом эффектов пластичности, ползучести и накопления повреждений;
формулировка нелинейных краевых задач исследования процессов деформирования и разрушения оболочечных конструкций при квазистатических термосилових воздействиях и схемы их численного решения;
комбинированная схема МКЭ и эффективные модели КЗ, позволяющие с единых позиций исследовать ВДС деформируемых систем, состоящих из тонкостенных и массивных образований на основе соотношений, наиболее соответствующих реализуемым в каждом из них видом ВД>,
комплексная методика, алгоритмы и, созданные на их основе, программные средства для численного исследования процессов деформирования и разрушения оболочечных конструкций;
результаты численного решения конкретных задач.
Апробация работы. Изложенные в работе результаты доложены на 71, УИЛ УШ Всесоюзных конференциях по прочности и пластичности (Москва, 1975,. Горький, 1978, Дермь, 1983), IX, XI, ХП, ХУ Всесоюзных конференциях по теории оболочек и пластин (Ленинград, 1973, Харьков, 1977, Ереван, 1980, Казань, 1990), Ш, ІУ, У Всесоюзных симпозиумах по малоцикловой усталости элементов конструкций (Вильнюс, 1979, Краснодар, 1983, Волгоград, І987), Всесоюзном симпозиуме "Ползучесть в конструкциях" (Днепропетровск, 1982, Новосибирск, 1984), 1,П Всесоюзной конференции "Численная реализация физико-механических задач прочности" (Горький, 1983; 1987), УІ.УП.УШ,IX Всесоюзных семинарах по комплексам программ математической физики (Днепропетровск,
1979, Горький, 1981, Ташкент, 1983, Шушенское, 1986), Всесоюзном симпозиуме по математическим методам МДТТ (Москва, 1984), Всесоюзной конференции по статике и динамике пространственных конструкций (Киев, 1985), Всесоюзной конференции по смешанным задачам механики деформируемого тела (Харьков, 1985), Всесоюзной иноле-севднаре "Математическое моделирование в науке и технике" (Пермь, 1986), Всесоюзном семинаре "Актуальные проблемы механики оболочек" (Казань, 1988), П.У.УШ школах по МКЭ (Горький, 1975, Рига, 1981, Усть-Нарва, 1987), П Всесоюзном симпозиуме "Прочность материалов и элементов конструкций при сложном НДС (Киев, 1984), П Всесоюзной конференции "Проблемы снижения материалоёмкости силовых конструкций (Горький, 1989), Всесоюзных конференциях "Теоретические основа конструирования численных алгоритмов решения задач математической физики"(Горький, 1986, Москва, 1990), Ш, XI, ХП Всесоюзных конференциях по численним методам решения задач теории упругости « пластичности (Новосибирск, 1974, Волгоград, 1989, Тверь, 1991).
Публикации. Основные результаты работа отражены в 41 публикации и 6 отчетах о ШР. В автореферате приведен список 20 основных работ.
Структура и объём работы. Диссертация состоит из введения, пяти глав, заключения, списка литературы и приложения. Работа содержит 297' страниц основного текста, 108 страниц иллюстраций (134 рис., 13 табл.). Список литературы включает 259 наименований. В приложении приведены акты о внедрении результатов работ.
Во введении приведен краткий обзор литературных источников, посвященных развитию рассматриваемых в работе основных направлений, показана актуальность намеченных исследований, сформулированы цель работы,её научная новизна и практическая ценность.
Изучению общих закономерностей развития процессов упруго-пластического деформирования и ползучести при термосилових яа-гружениях и развитию основных направлений в области создания математических моделей посвящэны работы А.А.Ильюшина, В.В.Нз-вожилойа, Ю.Н.Работнова, А.Ю.Иіалинского, В.С.Ленского, И.А.Бир-
гера, В.В.Москвитина, В.Д.Клюшникова.Н.Н.малинина, Прагера, Бейли, Мруза, Ольшака, Пежины, Надай, Ходжа и других.
Среди наиболее важных направлений в создании математических моделей пластичности следует выделить теории малых упруго-пластических деформаций, развитую в работах А.А.Ильюшина, Й.А.Биргера, В.В.Москвитина, Ю.Н.Шевченко, Генки, Надай, теорию упругопластических процессов, разработанную и развитую в работах А.А.Ильюшина, В.О.Ленского, Р.А.Васина, В.Г.Зубчанинова, В.И.Малого, Ю.Н.Шевчвнко, А.С.Кравчука, Н.Н.Столярова, группу структурных моделей, рассмотренных Мазингом, Бесселингом, Д.А.Гохфельдом и О.С.Садаковцм, Н.Н.Афанасьевым, В.С .Зарубиным, Ю.Н.Шевченко, а также болыцутэ группу теорий течения, различающихся между ообой различными способами процесса упругоплаотл-ческого упрочнения. Теория пластичности с изотропным упрочнением рассмотрена в работах Прагера, Й.А.Биргера; теории с кинематическим упрочнением - в работах А.Ю.Ишинского, Прагера, Ю.И.Кадашевича и В.В.Новожилова. Различные варианты теорий с комбинированным упрочнением рассмотрены в работах В.В.Новожилова и Ю.И.Кадашевича, Г.Б.Талнпова, А.А.Вакулешсо, Й.А.Биргера, И.В.Демьянушко и Ю.Н.Темиса, Ю.Г.Коротких, B.C.Бондаря и других.
Развитие этих моделей для учета явлений ползучести производилось по двум основным направлениям. Согласно первому -необратимые деформации разделялись на мгновенную пластическую составляющую и зависящую от времени составляющую деформаций ползучести, модели этого направления рассмотрены в работах Й.А.Биргера, И.В.Демьянушко, Ю.Н.Темиса, Ю.Г.Коротких, Ю.Н.Шевченко, Мруза и других. Согласно второму направлению - необратимая деформация не разделяется на мгновенную и временную составляющие и определяется на основе- обобщенной теории неупругости. Появление таких теорий и их дальнейшее развитие связано с именами Бейли, Орована, В.В.Новожилова и Ю.И.Кадашевича, И.З.Паллейі В.0.Бондаря, Миллера.
Следует иметь в виду, что наличие одних моделей пластичности и ползучести еще не позволяет получить прямой ответ на вопрос о текущем состоянии поврежденности исследуемых конструкций. При атом оценка их работоспособности решается приближенно на основе простейших критериев прочности. В то же время, современные исследования в области прочности материалов позво-
ляют представить разрушение как многостадийный процесс возникновения и развития необратимых дефектов, определяемый всей кинетикой НДС конструкций в процессе их нагруления.
Разработке принципиальных вопросов, основных направлений построения кинетических уравнений и конкретных моделей накопления повреждений посвящены работы А.А.Ильюшина, Ю.Н.Работнова, Л.М,Качанова, В.В.Нзвожилова, В.Л.Колмогорова, Б.Е.Прбедри,. Г.С.Писаренко, В.В.Моснвитина, В.В.Болотина, С.В.Серенсена, Р.М.Шнейдеровича, А.Н.Ромагава, Н.А.Махутова, А.П.Гусенкова, В.П.Когаева, С.А.Шестерикова, Ю.Г.Коротких,'В.С.Бондаря, А.А.Мовчана, О.Г.Рыбакиной, Р.А.Дульнева, П.И.Котова, П.А.Павлова, Мэнсона, Мруза, Боднера, Кремпла, Мартина, Мураками, Ша-боши, Леметра и других. Следует иметь в виду,- что в большей части современных работ, посвященных моделировании накопления повреждений,предполагается, что развивавшаяся в материале пов-режденность не сказывается непосредственно на характеристиках процесса деформирования. Однако такое предположение справедливо лишь на ранних стадиях процесса разрушения, а о развитием поврежденное ти она начинает все сильнее влиять на механические характеристики материала. Это приводит к необходимости учета при формулировке уравнений состояния материалов взаимного влияния эффектов деформирования и повренденности, т.е. применения соотношений механики поврежденной среды.
Разработке основных положений, подходов и конкретных моделей механики поврежденной среды, а также применению их для решения конкретных задач, посвящены работы Ю.Н.Работнова, Л.М.Качанова, Мураками, Шабоши, Хулта, Ю.Г.Корртких, В.С.Бондаря, Ю.Н.Шевченко і ь.Н.Мазура, работы автора и других.
Исследование поведения конструкций с учетом явлений пластичности, ползучеста и тем более в рамках соотношений механики поврежденной среды приводит обычно к сложной системе существенно нелинейных уравнений, решение которых, в большинстве случаев, осуществляется путем предварительной линеаризации задачи и сведении её к последовательности линейных краевых задач..Построению различных методов и численных схем решения нелйнёйгагх задач механики посвящены работы И.А.Биргера, В.В.Петрова, Э.Н. Іриголпка и В.И.Шалаашлина, Д.Ф.Давядеяно, М.С.Коршшйна* Б.Е.Победри, В.И.Мяченкова, А.Е.Фролова, В.А.Поотнова, Н.Н.Шапошникова, И.В.Демьянушко, Ю.В.Темиса, B.C.Бондаря, Е.М.Моро-
зова, В.Н.Кукуджанова, А.С.Сахарова, В.И.Гуляева, В.Н.Паймуши-на, Н.Н.Столярова, А.С,Городецкого, Ю.В.Липозцева, Зенкевича, Одена, Аргириса, Стриклина, Хейслера, Ривеманна, Маркала и других. Для решения линейных краевых вадач статического анализа конструкций в настоящее время наибольшее распространение находит МКЭ. При этом отмечается, что широкие возможности МКЭ для решения многих задач механики деформируемого тела, применение его для анализа оболочечнах конструкций, состоящих из фрагментов пластин, плит, оболочек различной толщины и массивных фрагментов, наталкивается на ряд трудностей .которые не получили должного разрешения.
Анализ современного состояния исследований в рамках рассмотренных в работе проблем показал, что наряду со значительными успехами в математическом моделировании процессов необратимого деформирования и разрушения конструкционных материалов, созданием методов и численных схем решения нелинейных и линеаризованных краевых задач, в настоящее время в области практического их использования для исследования реальных конструкций наблюдается заметное отставание. Число публикаций, посвященных исследованию процессов разрушения элементов и узлов конкретных * конструкций весьма невелико, особенно при исследования этих процессС'З в рамках соотношении механики поврежденной среды. Остаются актуальными разработка надежных моделей, позволяющих описать с единых методологических позиций эффекты пластичности, ползучести и развивающейся поврежденности с учетом их взаимного влияния в общем случае термосилового нагружения, построение экономичных и устойчивых схем решения нелинейных задач и создание на их обнове комплексной методики, связывающей воедино используемые для анализа конструкций частные модели, схемы, алгоритмы и обеспечивающей их оптимальное сочетание в составе единого вычислительного процесса.
Первая глава посвящена построению моделей, используемых в работе для описания процессов деформирования и разрушения конструкционных материалов. На основе анализа литературных данных рассмотрены основные закономерности пластического деформирования и ползучести начально изотропных конструкционных материалов при кваз иста тиче ских термосиловых воздействиях, зависимость их от условий нагружения и возможности описания этих
- II -
процессов в рамках феноменологических моделей. Рассмотрены основные закономерности развития повреаденности материалов в условиях пластического деформирования и ползучести, влияния на них температуры, вида ВДС, шаровой составляющей тензора напряжений и истории нагрукения, взаимодействия между собой различных зядов поврэдденяоста. Отмечено, что заметное влияние повреаденности на процесс деформирования сказывается лишь на завершавшей его стадии и быстро становится определяющим фактором этого процесса, проявляясь в образовании пластической нестабильности материала, при возрастании скорости деформирования на третьем участке кривой ползучести. На основе анализа рассмотренных моделей сформулированы перечень и структура соотношений, описывающих процессы пластичности, ползучести, накопления повреждений, проведен анализ и классификация входящих в них величин как формальных параметров, определяющих содержание соответствующих формальных моделей.
Предложен вариант составной иерархической модели поврежденного материала, описывающей процессы деформирования и разрушения начально изотропных конструкционных материалов с учетом взаимного влияния эффектов пластичности, ползучести и накопления повреждений. Модель устанавливает состав, структуру и функциональные связи между перемен:!*!!!, определяющими развитие названных эффектов как формально независимых процессов, описываемых, соответствующими частными моделями. Взаимодействие и взаимное влияние таких частных моделей при описании реальных процессов осуществляется в общей составной модели путем согласования их функций и общих параметров. При этом частные модели могут быть определены в общей модели достаточно формально, без детального описания способов реализация выполняемых функций и конкретных величин, используемых внутри частных моделей. Такой подход позволяет, с одной сторона, представить сложный процесс развития взаимосвязанных эффектов необратимого деформирования ' и поврежденности в виде совокупности более простых и формально несвязанных процессов, а с другой - использовать для описания этих процессов широкий набор имеющихся и разрабатываемых альтернативных средств без изменения формулировки общей модели поврежденного материала.
В частности, исходя из представлений современных теорий
пластичности, использующих понятие поверхности текучести, элементарное изменение пластических деформаций дбу № малом интервале изменения внешних воздействий, характеризуемых изменением температуры J Т и деформаций А в и (или напряжений й 6а) может быть однозначно определено, если известны исходный уровень действующих напряжений 61; температуры Т и набора скалярных и тензорных параметров %# { & = 1,2... п ), характеризуючих упрочнение материала и другие проявления истории упругопластического деформирования. К названному перечню должна быть добавлены необходимые материальные функции и комплекс уравнений, описывающих начальное и'текущее условия текучести, закономерности изменения пластических деформаций двft
Р
я параметров t . Обобщая функции таких уравнений можно ввести формальную модель' термопластичнооти, записываемую в виде:
В рамках модели (I) может быть рассмотрен широкий класс различных конкретных моделей термопластичности, отличающихся числом и конкретным содержанием параметров 1^ . Определяющие соотношения модели (I) записываются для неповрежденного материала и формально не включают в себя какой-либо прямой зависимости от текущей поврежденкости (эта зависимость проявляется лишь через уровень действующих напряжений &[ , фигурирующих в модели-в качестве входных параметров).
По аналогии с моделью пластичности формулируются функции формальной модели термоползучести: определить в точке тела изменение деформаций ползучести Аб"; и, текущих значений параметров ї( q. " 1,2... m ), характеризующих вязкое упрочнение материала на малом интервале изменения внешних воздействий, определяемом изменением деформаций ДЄ;; и температур л Т , происшедшем за промежуток вршени jt s t ~І » "PS. известном нсходном_состоянші, заданным уровнем напряжений 6 ц , температуры 7" и sкачением параметров 1^ :
АЄ*:деІ(еу,7,4Є^,і,Ц); (2)
- ІЗ -
При построении формальных моделей накопления повреждений предполагается, что в материале, в процессе его деформирования, могут независимо развиваться несколько различных типов повреж-денности, характеризуемых соответствующими функциями ^.( к -1,2... I ), Влияние каждого из них на характеристики процесса деформирования описывается о помощью скалярной функции сО (функции целостности), представляющей собой меру уменьшения эффективных площадок действия напряжении по отношению к их начальному неповрежденному значении ( u) = 0+1).
Предполагается также, что изменение повреждений каждого вида А У* , в общем случае, определяется уровнем действующих напряжений 6ц , температуры Т. . изменением необратимых деформаций Ав*: за рассматриваемый интервал нагружйния и значением некоторых параметров I^*( S = 1,2... р ) (в случае пластических повреждений ДЄК - Дёц , 1съ), а также значением накопленной поврежденное ти данного вида и константами материала W* :
аг-^(6^ ,7Л, ^,1^^), (3)
Вклад поврежденности каждого вида в изменение функции целостности представляется в виде:
где сО - накопленное значение функции сО , О. » ( Ji - 1,2... ) - некоторые константы материала.
Полное значение ей , соответствующее хонду рассматриваемого шага нагруже ния определяется на оонове используемого в конкретной модели алгоритма суммирования прирааззниЯ А и)к от каждого вида поврежденности. С учетом рассмотреяной выше сруктуры уравнений накопления повреждений й входящих в них величин, формальную модель накопления повреждений можно представить в виде:
В рамках формальной модели (5) может быть рассмотрен ши
рокий класс кинетических уравнений, описывающих накопление по
вреждений. При атом для получения рабочей модели достаточно
конкретизировать соотношения (3) и (4), алгоритм суммирования
повреждений, а также установить число и конкоетный вид пара-
метров г? . '
Рассмотренные формальные модели используются как составные части обшей модели поврежденного материала, функциональная совместимость которых может быть обеспечена при условии согласования состава и структуры входящих в них внешних параметров. Ввиду того, что в (I) + (5) все величины, кроме параметров Z#, ta U If определены однозначно для любых моделей, рассматриваемых в рамках формального подхода, условия согласованности сводятся к тому, чтобы входящие в модель накопления повреждений параметры 1Н. ' содержались в числе параметров . %* и t в, конкретных моделей пластичности и ползучести.
При построении общей (составной) модели поврежденного материала, связывающей изменения напряжений и деформаций на интервале изменения внешних воздействий, в рассмотрение введены две характеристики напряжений: аффективные напряжения 6г > действующие на поврежденных площадках и приведенные 6г , статически эквивалентные первым, но отнесенные к неповрежденным площадкам. Первые фигурируют во всех частных моделях,определяющих состояние материала в точке тела, вторые используются на уровне описания исследуемого объекта при формулировке уравнений равновесия и статических граничных условий. Поэтому определяющие соотношения модели поврежденного материала записываются для изменений приведенных напряжений:
aG*=G*-G*\ G*=(1-b)G; G*=(i-u»Gi
АК*=К*-ІЇ*; К*=(1-їд)К; К*:(і-й)К;
№ к = іс(Т); к=К{7)\ 3=G-(T)',GsG(D-um7M 0бъёи'
ной и сдвиговой деформации неповрежденного материала, отнесен,-ные к уровню температуры в исходном и текущем состояниях, МеСГу изменение тепловой деформации на интервале нагружения.
Рассмотрены вопросы конкретизации соотношений предложенных формальных моделей с использованием известных моделей пластичности, ползучести и накопления повреждений. В частности, получена полная система уравнений модели поврежденного материала и дана конкретизация всех входящих в неё параметров на основе соотношений модели термовязкопластичности с комбинированным упрочнением и уравнений накопления повреждений, развиваемых в работах Ю.Г.Коротких. При описании поврежденности материала предусмотрена возможность его разрушения за счет пластического накопления повреждений и хрупкого разрушения. С зтой целью в первом случае использовались кинетические уравнения, определяющие изменение энергии пластического разрыхления, типа:
dty = dW/W*; dWyt]del i7)
где f:j - текущее значение тензора остаточных микронапряжений, {(П) - Функция вида НДС ( Qif(fl)
- tl
of«), =у-и) t'dty , С 8 )
где Ц, - константа материала.
Для описания изменения поврежденности при хрупком разрушении использован критерий максимального нормального напряжешш:
*
Ad)z = О , если 6Х<6*0 ( g )
Л и>: s(f-iO), если 6j}6t
где 60=6 (Т)~ разрушавшее напряжение неповрежденного материала при одноосном'раотяжении. Определение текущего значения функции сд строилось на основе линейного правила суммирования повреждений, определяемых соотношениями (7) и (8). Приведенная конкретная модель использовалась в качестве базовой при проведении рассмотренных в работе численных исследований процессов дефор-
мирования и разрушения конкретных конструкций.
Вторая глава посвяиена формированию и решению систем уравнений, описывающих поведение конструкций, состоящих из элементов пластин, плит, оболочек различной толщины, а такие фрагментов массивных образований при, квазистатических термосиловых на-груяениях на основе полученных в первой главе физических уравнений в геометрически линейной и нелинейной постановках. Ввиду существенной нелинейности этих уравнений, а также в связи с тем, что значительная часть входящих в них величин является функционалами, зависящими от предшествующей истории процесса деформирования и может быть определена лишь путем интегрирования для конкретных траекторий нагрукения, при построении разрешающей системы использована инкрементальная формулировка исходных уравнений.
С целью удобства изложения в первых разделах главы рассмотрены вопросы исследования поведения конструкций в геометрически линейкой постановке. Получена система уравнений, описывающих поведение конструкций для конечного интервала изменения внешних воздействий (этапа нагружения), характеризуемого значением температуры 7* и объёмных сил F- ( і =1,3), заданных в объёме конструкции V . поверхностных сил Р; , заданных на части поверхности , , а также сметаниями граничных_поверх-ностей uf на части поверхности Z-z в исходном {t=t,Т, F^ Р[,й^) и текущем (tst,T, ?i, Рс,й* ) состояниях при известных в исходном состоянии начальной конфигурации конструкции ее . , действующих напряжений 6п и параметров, характеризующих историю необратимого деформирования материала ( Х#, Z\ , йо ) В состав этой системы включены полученные в первой главе физические соотношения, линейные геометрические уравнения, связывающие изменения деформаций Аг- и перемещений Ли і , уравнения равновесия,,статические и геометрические граничные условия, записанные для изменений напряжений дбу и граничных смещений 4U? на этапе нагружения и две дополнительные группы уравнений. Необходимость введения первой группы вызвана возможностью нарушения непрерывности перемещений и условий равновесия вдоль некоторых поверхностей ,2Гд і являющихся границей раздела смежных конструктивных фрагментов, в которых поля перемещений определяются на основе различных допущений. При этом для каждой из границ J ' ' разделяющих -такие фрагменты „о<" и „ д" записываются
дополнительные уравнения
au"-au{ =0 , Щ.-Є; +Аб:. і =0 ; ( Ю )
где с и с; -направляющие косинусы фрагментов „ос и „/3 к границе сопряжения.
Вторая группа уравнений вводится для дополнительной коррекции глобальных условий равновесия, которые могут нарушаться в исходном состоянии за счет накопления погрешностей при пошаговом решении задачи
6LJiJ+F.=o SV, 6lflrP.=0 на Z, (II)
Дальнейшее применение перечисленных уравнений связано с их линеаризацией в рамках метода начальных напряжений, позволяющего свести решение нелинейной задачи на шаге нагружения к последовательности линейных задач с уточняемыми в процессе приближений величинами, определяемыми эффектами поврежденности и необратимого деформирования материала (коэффициентами Ad;: в уравнениях (6). Разрешающая система линеаризованных уравнений, описывающая поведение конструкций в текущем приближении (при фиксированных значениях коэффициентов ДСІц) получена из условия стационарности функционала:
n*?U6ii иец-Мрс/У-ЬЪщ dV-. dZ+
( 12 )
V J J J v Z1
+ v
Jce.jAe.-F^u^dV-j^alZ^Xiu^Z;
V J J *i 2x
В представленном функционале изменения напряжений дб;; и деформаций AS с; связаны перечисленными выше физическими и геометрическими уравнениями с изменением перемещений AUi, удовлетворяющих заданным кинематическим граничным условиям на части поверхности Z2 , {иіУ'Рх/іи"-A U. - скачок функций перемещений на внутренних границах 2Jt'^ некоторых конструктивных образований вдоль которых эти функции могут претерпевать разрыв, Л; - соответствующие множители Лагранжа, имеющие смысл поверхностных сил, распределенных вдоль поверхностей 5Г\.
Проведен краткий анализ существующих шаговых методов решения задач исследования неупругого поведения конструкций при квазистаткческих нпгрухсчялх, рассмотрена преимущества и недо-
статки различных конкретних формулировок, отмечены основные источники накопления погрешностей, возникающих при численном решении задач на основе шаговых методов в возможные пути снижения этих погрешностей. Предложена двухуровневая шаговая схема решения нелинейных уравнений, описывающих процессы термовязко-пластического деформирования конструкций, обеспечивающая высокую точность и устойчивость вычислений при существенном снижении их трудоемкости. Шаги каждого уровня в этой схеме несут различную функциональную нагрузку. На шагах верхнего уровня (этапах нагружения) осуществляется внешняя линеаризация задачи. Реальная траектория нагружения представляется'в виде совокупности прямолинейных участков. Величина шага здесь определяется только условиями удовлетворительной аппроксимации траектории нагружения и для монотонных процессов может быть достаточно большой, вплоть до конечного изменения внешних воздействий. Решение нелинейных задач на каждом этапе нагружения осуществляется в форме метода начальных напряжений путем итерационного уточнения равновесного состояния системы на основе схемы промежуточных экстраполяции, позволившей значительно ускорить сходимость процесса последовательных приближений. Для определения изменений необратимых деформаций внутри этапа нагружения строится внутренняя шаговая схема, позволяющая с заданной точностью вычислить значения скоростей их изменения (или их малых приращений) для внутренних точек траекторий деформирования и проинтегрировать в пределах текущего этапа. На внутренних шагах все вычисления строятся независимо для каждой точки тела, где ожидается изменение необратимых, деформаций без коррекции уравнений равновесия, поэтому увеличение числа внутренних шагов не сказывается существенно на увеличение трудоемкости задачи. Размеры шагов внутренней схемы выбираются автоматически путем согласования прогнозируемых величин изменения напряжений на текущем шаге с некоторым заданным допуском.
При решении задач термопластичности в каждой точке тела, где ожидаются изменения необратимых деформаций, текущий этап нагружения разбивается на Д тагов протяженностью otz (2lSij'l). размеры которых выбираются из условий, чтобы весь участок упругого деформирования (если таковой имеет место)' полностью укладывался в один первый шаг О ti , а последующие шаги не приводили к превышению прогнозируемых величин изменения
напряжений некоторой заданной доли текущей поверхности текучести. Для каждого текущего шага вычисляются соответствующие изменения полных деформаций (.5ві;)і -&ij6tx и температуры 7^г/г-?+ aTqTx , а затем на основе выбранной модели пластичности, вычисляются приращения пластических деформаций (<$вс;)і . текущих значений напряжений и параметров, характеризующих упрочшние материала. Изменение пластических деформаций на этапе получается в результате суммирования их приращений на внутренних шагах, при решении задач ползучести так же, как и для термспластичнос-ти, текущий этап, характеризуемый временной протяженностью A t . подразделяется наряд шагов, длительностьюSt!!(At:iZu{8t1) из условия, чтобы изменение напряжений в рассматриваемой точке за текущий шаг не превышала заданной доли радиуса поверхности текучести. Далее на основе выбрангой модели ползучести вычисляются текущие скорости ползучести (P('j)j и их приращения (О6у )г (интегрирование скоростей в пределах этапа осуществляется на основе метода Эйлера с модификацией типа "предиктор-корректор"), после чего корректируется текущий уровень действующих напряжений.и параметров, характеризующих упрочнение материала. При решении задач термовязкопластичности для каждого текущего шага производится независимое вычисление протяженное те й шага из условий пластичности (8tP)l ползучести {StH)x и выбирается наименьшая из протяженностей 6ti $(Sth)i , (dt^i Далее, с учетом выбранной протяженности шага, производится вычисление приращений деформаций ползучести ОЄц , корректируется уровень текущих напряжений и осуществляется аналогичный шаг по пластичности. После этого производится вычисление протяженности и приращений необратимых деформаций для следующего шага.
Применение описанной шаговой схемы позволило сократить трудоёмкость решения задач термовязкопластичности за счет возможности значительного увеличения размера шагов верхнего уровня. С целью дальнейшего снижения трудоёмкости в работе предложены специальные приемы ускорения сходимости итерационных процессов решения задач на этапе нагружения. Первый из них основан на том, что для прост'ейших задач пластичности, решаемых на основе метода начальных напряжений, сходимость процесса последовательных приближений происходит по закону, близкому к закону геометрической прогрессии. Это обстоятельство позволило применить для ускорения сходимости подобных итерационных процессов преобразо-
вание Эйткина для определения предельного значения последовательностей, меняющихся по законам, слизким к геометрической прогресом или показательной функции вида:
5"*3 + А(9)\ (13>
где 5 - постоянная часть функции 5-*5 при п~>оо . предел последовательности S такой функции может быть найден по трем различным значениям 5 :
5 - "' ; ( 14 }
Применение преобразования (14) к последовательности векторов узловых смещений конструкции, определяемых в процессе итерационного уточнения равновесного состояния, позволяет найти приближенное значение этого вектора по значениям его компонент на трех смежных итерациях.
Наряду с возможностью значительного снижения трудоёмкости решаемых задач за счет сокращения общего числа итераций рассмотренная схема может бать эффективно использована для' экстраполяции решения по параметру нагруяения на основе результатов решения задачи для трех смежных точек. В этом случае сначала находятся значения коэффициентов 5 . А и Q, в последовательности (ІЗ), а затаї находится искомое значение SN для заданной величины параметра нагруяеняя N .
Как показал опыт решения большого числа конкретных задач, эффективность применения рассмотренной схемы существенно зависит от того, на сколько установилось в процессе последовательных приближений монотонное изменение всех,компонент решения. Это обстоятельство ограничивает область практического примзшша схемы, т.к. в процессе упругопластического деформирования и ползучести реальных конструкций возможно значительное перераспределение напряжений и деформаций. С целкю преодоления отмеченных недостатков в работе предложена ещё одна схема ускорения сходимости, названная схемой промежуточных экстраполяции. Для её построения используется разложение уравнения/ :Ни-Х = 0 , описывающего нелинейную часть перемещений и на текущем этапе нагружения (i:l(U) , U*U+U . где \JC- линейная часть решения, определяемая упругим приближением задачи) в ряд Тей-
пора в окрестности точки и :
fzfh + (tt)^U-.HAU-5z"-Jlh=0, ( is)
где Н - линейный оператор, определяемый упругими характеристиками материала в текущем состоянии,St = -Hu"t lhCU") вектор невязки уравнения в ґі -ом прибліяешіи,4гп=('Л/^а)'^^. Пренебрежение в (15) величиной AZ* приводит к простейшей схеме последовательшх приближений в рамках метода начальных напряжений:
Ди*6и»:Н'Ъ" иПН=и"+Ди (16)
Ввиду сложности и значительной- трудоемкости непосредственного вычисления вектора а X л для его определения делается предположение о коллинеарности векторов Л" и^Е*1 или о возможности представления Д$п- &t*+8l"-m-'8%* , где т -некоторый скалярный множитель. Значение его определяется как среднеквадратичная величина по значениям /г>; (со специально подобранными весами), вычисленным на основе преобразования Эйт-кина для каждой из компонент решения. Прн этом для каждого шага предполагаемой экстраполяции сразу же после вычисления значений /ті- осуществляется проверка выполнений некоторых условий монотонности для каадой из компонент решения. В случае невыполнения какого-либо из этих условий вычисленное значение пі* корректируется соответствующим весом и текущая компонента считается "дефектной". Если отношение числа таких "дефектных" компонент О. к размерности вектора решения /V превысит заранее установленную норму $ (^//V Ъ Q ), то процесс приближений повторяется по схеме (16) до выполнения указанного условия. Такая процедура экстраполяции осуществляется параллельно с процессом приближений по схеме (16) до достижения заданной точности, причем каждый последующий шаг производится не ранее, чем через три итерации после предыдущего. Благодаря своей простоте, возможности однократного формирования и обращения матриц жесткости определяющих систем, а также значительному сокращению числа итераций, применение рассмотренной схемы оказалось чрезвычайно эффективным для решения широкого класса нелинейных задач.
На основе изчоженяых методических положений предложены
алгоритмы построения вычислительного процесса для типовых постановок задач с использованием схем экстраполяции решения по параметру нагружения, дополнительной коррекции глобальных условий равновесия и других приемов, позволивших, в приемлемое для инженерных расчетов время, исследовать процессы деформирования и разрушения реальных конструкций.
В заключительном, шестом параграфе главы сделано обобщение рассмотренных методик и алгоритмов для исследования поведения конструкций с учетом конечности деформаций, перемещений 7 углов поворота. Для характеристики состояния конструкций в процессе их нагружения, по аналогии с геометрически линейной постановкой, использованы три отсчетшх конфигурации, определяемые координатами точек, составляющих конструкцию: начальная xf = X-L (t) . исходная xL = X-Jt) и текущая # = Хі (і) В процессе пошагового решения задачи исходная конфигурация 5С; используется в качестве начальной Лагракжевой системы на шаге нагружения для определения текущего состояния и текущей конфигурации Хі в конце шага. В качестве меры изменения деформаций на шаге использован модифицированный тензор приращений деформаций Грина, определенный в метрике исходной конфигурации. Для характеристики напряженного состояния конструкции использован модифицированный тензор напряжений Кирхгофа Q-- = 6с;+Дбіі'< где бц - тензор напряжений Эйлера в исходном состоянии перед приложением очередного шага нагрузки, a6;j- приращение тензора напряжений Кирхгофа на шаге нагружения, отнесенного к исходному состоянии. В конце каждого шага тензор ф.. преобразуется в тензор напряжений Эйлера, который используется в качестве тензора начальных напряжений для следующего шага:
fy"" ш~: э1* it ty, 3tJ zd&i/axj (I?}
Разрешающая система линеаризованных уравнений, описывающих геометрически нелинейное поведение конструкций на шаге нагружения, по аналогии с линейной постановкой, получена из условия стационарности функционала типа (12) с добавлением в последний членов, учитывающих влияние начальных напряжений. Численная реализация методики осуществлена на основе алгоритмов и в рамках общей схемы организации вычислительного процесса, предложенных для геометрически линейной постановки, дополненных
ограничениями размера шага верхнего уровня, корректировкой на каждом из них текущей геометрии конструкции и глобальных условий равновесия её текущей конфигурации.
Третья глава посвящена построению решений линеаризованных задач расчета конструкций, состоящих из фрагментов пластин, плит и оболочек различной толщины, а также массивных, образований на основе МКЭ.
Рассмотрены вопросы применения МКЭ для анализа тонкостенных плит и оболочек, составлявших наиболее распространенный тип конструктивных фрагментов в составе конструкций рассматриваемого класса. Отмечены основные трудности формулировки моделей КЭ для исследования пластин и оболочек, связанные о выполнением условий межэлементной совместности, аппроксимацией искривленных поверхностей, смещений элементов как жесткого целого, сдвигового и мембранного запирания. Рассмотрены различные подходы и конкретные формулировки МКЭ на основе классических соотношений теории тонких оболочек, соотношений, учитывающих деформации поперечного сдвига, а также общих соотношений трехмерной теории. Детально исследованы вопросы применения для анализа плит и оболочек пространственных изопараметрических элементов при расположении последних в один ряд по толщине. Остановлены причины низкой точности и плохой сходимости конечно-элементных решений, получаемых на основе простейших элементов этого типа, заключающиеся в нарушении условий полноты аппроксимирующих функций из-за несогласованности их распределения при использовании неполных полиномов. Показано, что аппроксимация функций в элементе неполными полиномами приводит к появлению ошибок в представлении сдвиговых компонент, которые существенно снижают точность и сходимость КЭ решений при относительном уменьшении одного из размеров элемента. Показано также, что повышение порядка аппроксимирующих функций позволяет последовательно улучшать показатели точности КЭ решений, однако при этом резко возрастают трудности вычислительного характера. В настоящее время существует много других подходов, позволяющих повысить эффективности элементов без существенного повышения их сложности: способы доведения аппроксимирующих функций до полных полиномов путем включения дополнительных узловых илі{ внеузловых неизвестных, способ двойной аппрокси-
мании, моментная схема конечных олементов и различше варианты метода сокращенного интегрирования. Метод сокращенного интегрирования, предложенный вначале на основе интуитивынх соображений, получил в настоящее время достаточную физическую и математическую интерпретацию и рассматривается как некоторая техно- , логия корректировки свойств КЭ без изменения принятых координатных функций. Для преодоления произвола в выборе порядка и точек численного интегрирования при формулировке конкретных моделей КЭ в работе предложен подход, позволяющий получать высокоэффективные модели без нарушения основных условий, обеспечивающих сходимость и общепринятых технологий построения матриц жесткости. Согласно этому подходу интегрирование должно осуществляться в точках элемента, где достигается наивысшая точность вычисления напряжений и деформаций и обеспечивается точное вычисление объёма и всех производных от полных полиномов, входящих в аппроксимирующий полином. При этом не должно появляться сингулярности разрешающей системы алгебраических уравнений в связи с появлением в элементах форм перемещений, соответствующих нулевой энергии деформации. Для КЭ, функции которых аппроксимируютмя полными полиномами, точки наилучшего определения деформаций совпадают с точками, обеспечивающими точнее интегрирование всех членов энергии деформации. В КЭ сирендипо-ва типа, функции перемещений которых аппроксимируются неполными полиномами, такие точки совпадают с'корнями полиномов 1е-жандра, степень которых определяется степенью полного полинома, используемого для аппроксимации функций в элементе. Иначе говоря, эти точки совпадают с точками .квадратур Іаусса, по которым точно интегрируется часть энергии деформации, определяемая полными полиномами, входящими в выражения аппроксимации функций в элементе. При этом сохраняется максимально возможный для рассматриваемых элементов порядок сходимости. Особенно благоприятным применение такого подхода оказывается для квадратичных элементов сирендипова типа, т.к. возникающая при этом сингулярность не проявляется в системах, состоящих более, чем из одного элемента и в отдельном элементе исчезает при задании кинематических граничных условий. Поэтому в работе, в качестве базовых моделей, используемых для анализа массивных толстостенных и даже ряда тонкостенных фрагментов оболочечных конструкций, предложено использовать семейство изопарвметричес-
ких квадратичных КЭ о сирендігаовой аппроксимацией перемещений а рассмотренной схемой численного интегрирования. Проведено исследование точности, сходимости и условий применимости таких моделей на основе энергетических оценок, многочисленних результатов численных экспериментов и оценки трудоёмкости вичислеїшй, показавшее их высокую эффективность при решении рассматриваемых классов задач.
Для анализа тонкостенных фрагментов в составе оболочечннх конструкций в работе предложен изопараметрический сдвиговой конечный элемент оболочки общего вида с квадратичнім законом распределения по элементу смещений срединной поверхности и углов поворота сечений. Элемент имеет сорок степеней свободы и представляет собой участок гладкой оболочки переменной толщины в виде топологического четырехугольника с четырьмя узлами в углах и четырьмя узлами на срединах ребер. В отличие от аналогичных элементов, предложенных ранее Ьенкевичем,- Ахмадом, Pu\VSSy и другими, построенных на основе трехмерных уравнений, предлагаемый элемент базируется на соотношениях теории оболочек, учитывающей деформации поперечного сдвига. При этом наряду с некоторым снижением общности модели (неучета изменения метрики по толщине элемента) удалось значительно сократить трудоемкость вычислений (прежде всего за счет независимого интегрирования по толщине жесткостей оболочки, а также дополнительных сил и моментов, обусловленных изменением необратимых деформаций) и повысить его эксплуатационные возможности при анализе составных конструкций. Приведенные в работе результаты многочисленных исследований, связанных с возможностями геометрической аппроксимации, возможностями описания краевых эффектов и диапазона относительных толщин, а также решением различных тестовых задач, позволяют судить о высокой вычислительной эффективности предложенного элемента.
Наряду с моделями КЭ, предназначенных для анализа типовых фрагментов конструкций, в работе предложена комбинированная схема, позволяющая в рамках единой расчетной схемы исследовать ВДС конструкций, состоящих из различных тонкостенных и массивных образований с учетом эффектов концентрации напряжений в районах геометрических особенностей и зонах с высокими градиентам паї »зок и температур. Исследуемая деформітруепая система расчленяется на ряд подгсоксгрукций - конструктивных фраг-
ментов достаточно простой геометрической форма с характерным для каждой из них материалом и видом .НДС. Каждая подконструк-ция, в свою очередь, аппроксимируется совокупностью КЭ определенного типа, наиболее соответствующего рассматриваемому виду ВДС. Для описания сложных видов НДС в массивных деталях, переходных узлах, отдельных участках тонкостенных деталей, содержащих зоны концентрации напряжений, используются пространственные КЭ. При описании ВДС в тонкостенных участках конструкций, где без большого ущерба для точности могут быть использованы' соотношения прикладных теорий, используются КЭ оболочки. Необходимые условия сопряжения подконструкций, набранных из КЭ различных типов, реализуются в вариационной форме путем введения в исходный функционал соответствующих множителей Лагранжа, которые представляются в виде интерполяционных функций вдоль границ сопрягаемых элементов. Такая формулировка условий сопряжения позволяет избавиться от появления в районе сопряжения-возмущений, связанных с различиями локальных распределений перемещений в сопрягаемых подобластях, заменив локальные условия совместности интегральными. При этом вдоль каждой из сопрягаемых границ приравниваются не сами перемещения, а их интегралы по поверхности границы, взятые с некоторыми весовыми коэффициентами, определяемыми принятым законом распределения множителей Лагранжа. Такая схема позволяет легко осуществлять сопряжение самых разнородных систем и обеспечить гибкое выполнение условий совместности: от полностью нестесненного состояния до наложения жестких связей. С целью сохранения хорошей узколенточной структуры разрешающей системы алгебраических уравнений, параметры интерполяции и часть компонент узловых перемещений, принадлежащих сопрягаемым подконструкциям, исключаются из числа варьируемых параметров. В силу локальности распределения функций множителей Лагранжа такое исключение удается осуществить аналитически для каждой типовой поверхности сопряжения и получить конкретные соотношения, связывающие узловые параметры сопрягаемых КЭ. В результате реализация условий сопряжения подконструкций, набранных из различных типов КЭ в комбинированной схеме МКЭ, сводится к преобразованию векторов узловых перемещений, узловых сил и матриц жесткости отдельных элементов, примыкающих к границе сопряжения, которое выполняется на уровне сборки матриц жесткости подконструкций. На основе та-
кого подхода получены конкретные соотношения, используемые для реализации условий сопряжения предложенных в работе моделей КЭ и проведены численные исследования, подтвердившие высокую эффективность применения комбинированной схемы для решения конкретных задач.
Применение предлагаемой схемы для анализа НДС реальных конструкций принципиально позволяет получать решения задач с нужной степенью точности, однако для её реального использования необходимо иметь определенные представления о том, в каких случаях и к каким ошибкам могут привести соотношения прикладных теорий, а также о зависимости этих ошибок,от характерных параметров геометрии и нагружения исследуемых конструкций. Несмотря на большое число работ, посвященных вопросам обоснования и применения соотношений теории оболочек, конкретные оценки погрешностей получены лишь для небольшого числа частных задач. Поэтому в работе, не претендуя на полноту решения этой проблемы в общем плане, получены качественные и количественные оценки погрешности теории оболочек для некоторых типовых случаев локального распределения напряжений и деформаций в процессе упругого и упрутопластического нагружения оболочечных конструкций. В частности, рассмотрены погрешности определения напряжений в оболочках при действии локальных нагрузок (мягкие локальные нагружения) в зависимости от размера области приложения нагрузки и относительной толщины оболочки, а также при различных условиях локальных ограничений (жесткие локальные ограничения).
В результате проведенных исследований было установлено, что при мягких локальных нагружэниях по участку, характерный размер которого 0.)(1.5 r2.0)h (где h - толщина оболочки) на основе прикладных теорий можно получать достоверные результаты даже для сравнительно толстых оболочек. Погрешность определения максимальных значений компонент напряжений О при этом не превышает относительной толщины оболочки о ~t . При уменьшении размера области приложения нагрузки О. <1,5h разница в результатах, полученных на основе теории оболочек по сравнению с точным решением начинает резко возрастать и становиться особенно заметной в процессе появления и развития пластических деформаций. В районах жестких локальных нагруже-ний закон расггре.пэления напряжений и деформаций носит очень
сложный характер и не может быть смоделирован в рамках соотношений прикладных теорий. Размер зоны затухания возмущений X в обоих случаях локальных нагружений не превышает величины Z s (1.0 + 1.5) А . Как показали результаты исследования конкретных конструкций, все другие случаи локальных возмущений, возникающих при анализе оболочечных конструкций в районах высоких градиентов температур, геометрических концентраторов, подкрепляющих элементов, угловых соединений и др. можно рассматривать как промежутоные варианты мягкого и жесткого локальных воздействий.
Четвертая глава посвящена вопросам создания вычислительного комплекса, реализующего предложенные в работе методические положения для численного моделирования процессов деформирования и разрушения оболочечных конструкций на современных ЭВМ. Рассмотрена общая архитектура комплекса, представляющего собой интегрированный пакет, состоящий-из нескольких, формально независимых компонент: блока формирования информационной модели исследуемых конструкций (препроцессор), блока решения краевых задач (процессор), блока, обеспечиааютцего вывод и визуализацию информации об исследуемой конструкции на различных этапах решения задачи (постпроцессор), а также средств информационного обеспечения, поддерживающего связь всех фрагментов комплекса в едином вычислительном процессе, реализующего хранение, накопление и доступ к данным, используемым в процессе решения задачи. Рассмотрены функциональное назначение каждой из компонент комплексен состав входящих в них программных средств.
Проведен анализ состава и специфики преобразований информации, порождаемой решением рассматриваемых в рамках настоящей работы задач и на основе его предложены принципы управления данными, определяющими состояние исследуемых конструкций в процессе их нагружения.и данными по физикомеханическим характеристикам конструкционных материалов, необходимыми для реализации используемых моделей материала.
Пятая глава посвящена описанию" некоторых результатов численного моделирования процессов деформирования и разрушения оболочечных конструкций при квазистатических термосилових вдгружениях, полученных на основе рассмотренных в работе методических и программных средств.
Приведеш результаты исследования закономерностей упруго-пластического я упруговязкопластического деформирования материала в простейших элементах оболочечных конструкций для ряда типовых нагружений (цилиндрическая оболочка при действии кольцевой нагрузки, меняющейся по закону симметричного цикла, при различных историях термосилового нагружения, ребристая оболочка при действии давления и движущегося температурного поля, трубчатый образец при жестком термосиловом циклическом нагружений с выдержками), позволившие выявить возможности описания поведения конструкций в реальных условиях термосилових нагружений на основе разработанных моделей, методик и численных схем.
Рассмотрены задачи исследования процессов упругопласти-ческого и упруговязкопластического деформирования конструкций, состоящих из оболочек различной толщины и массивных фрагментов при осесиммегричных термосиллвых нагружениях (замкнутый сосуд с полостью между нагружной и внутренней стенками', используемой в качестве рубашки подогрева, осесиммегричная конструкция, состоящая из набора подкрепленных ребрами тонких оболочек и массивного шпангоута, корпус теплообменного аппарата, представляющий собой сочлененную систему цилиндрических, конических, то-рообразных оболочек и толстостенной плиты). При дискретизации таких конструкций в их составе выделялись зоны, в которых соотношения прикладных теорий могли привести к погрешности в описании ВДС и для их анализа использовались КЭ осесимметричного тела; оставшиеся части аппроксимировались КЭ оболочек. В результате проведенных исследований установлены характерные закономерности распределения напряжений и деформаций в переходных узлах конструкций, а также особенности развития необратимых деформаций в процессе их нагружения. В качестве приложений к решению задач термовязкопластичности рассмотрены задачи численного моделирования процессов образования остаточных напряжений при сварке осесимметричных конструкций. При этом предполагалось, что сварные швы в конструкциях формируются в результате нескольких кольцевых проходов с помощью мощного быстродвижущегося источника температуры, в результате чего распределение температурного поля вдоль окружной координаты можно считать постоянным, а соответствующую задачу термовязко-пласгичносгл - осе симметричной. Па основе этих предположений проведено исследование сварочных деформаций и напряжений в
- зо -
цилиндрической оболочке, сваренной из двух обечаек, выполненных из алюминиевого сплава, нержавеющей стали, а также сферической крышки свариваемой с плоским фланцем. Сравнение результатов, полученных на основе численного моделирования с известными экспериментальными данными показало их хорошее качественное и количественное согласование. В частности, подтвержден известный экспериментальный факт о выпучивании наружу сварного шва оболочки, сваренной из алюминиевых обечаек, а из стальных - внутрь.
Приведены некоторые результаты исследования НДС оболочеч-ных конструкций в пространственной постановке: конической оболочки, подкрепленной системой кольцевых ребер, сферического купола переменной толщины с нецентрально врезанным коническим патрубком, цилиндрического сосуда со сферической крышкой, имеющей четыре нецентральных круговых выреза.
Проведено исследование процессов упругопластического де-формированияи разрушения цилиндрического образца, представляющего собой сочлененную систему цилиндрических оболочек различной толщины с массивным фланцем в условиях осевого растяжения по схе?.іам мягкого и жесткого нагружений. Моделирование процесса осуществлялось в геометрически лилейной и нелинейной постановках до полного исчерпания образцом его несущей способности. В результате проведенных исследований отмечена значительная неравномерность развития пластических деформаций и поврежденнос-ти материала образца, особенно заметная на заключительной стадии процесса и вызванная постоянным перераспределением внутренних усилий, несмотря на монотонный характер, нагружения. Установлено, что непосредственное влияние поврежденности на процесс деформирования и разрушения начинает заметно сказываться лишь при достижении функцией целостности сО величин порядка сО * (0.005 + 0.01), а потеря несущей способности образца происходит при значениях этой функции в ряде точек значительно меньших своего предельного значения. Максимальное и) и среднее по толщине сечения иЭер значения функции целостности в зоне локализации поврежденности к моменту разрушения составили и)тая 0.053, vdcfl~ 0.046. Установлено также, что результаты численного моделирования и разрушения образца, полученные с использованием геометрически нелинейных соотношений, хорошо согласуются и имеющимися экспериментальными данными
- ЗІ -
как по уровню деформаций в процессе нагр ужения, так и по з на- * чению разрушающей нагрузки. Нзучет геометрической нелинейности позволил получить качественно приемлемые результаты, однако количественная погрешность при этом в определении предельных перемещений составила Suz 1455, а предельного значения осевой силы і^г9.$.
Аналогичные исследования проведены для сферического купола, находящегося под действием поперечной нагрузки, распределенной по малому участку. Материал купола - чугун, характерной особенностью которого является низкий предел прочности при растяжении, поэтому в зонах растяжения купола элементарные акты разрушения носили хрупкий характер. В результате численного моделирования получена подробная картина развития процессов упругопластического деформирования и разрушения купола до исчерпания им несущей способности. При атом отмечено, что результаты численного решения задачи на основе линейных соотношений ещё хуже согласуются с результатами, полученными на основе более точной нелинейной постановки. В частности, разница в определении предельной нагрузки на основе этих двух постановок превысила 69$.
В приложении приведены акты о внедрении.