CC BY
257
УДК 539.3
Исследование напряженно-деформированного состояния
ТОНКОСТЕННЫХ ЭЛЕМЕНТОВ КОНСТРУКЦИЙ
О. В. СТАРОЖИЛОВА
ФГОБУ ВПО «Поволжский государственный университет телекоммуникаций и информатики» (ПГУТИ)
Поступила в редакцию 14.06.2013
Аннотация. Статья посвящена численному моделированию тонкостенных элементов конструкций в виде пластин и оболочек, исследованию и построению нелинейных математических моделей напряженно-деформированного состояния. Модели учитывает сжимаемость материала, реальный вид диаграммы и позволяют исследовать деформирование гибких неоднородных оболочек переменной жесткости при поперечном, продольном и комбинированном нагружении.
Ключевые слова: пластины; оболочки; нелинейное деформирование; переменная жесткость.
В современной авиационной технике широко используются тонкостенные элементы конструкций в виде пластин и оболочек, обеспечивающие высокие прочностные показатели при достаточной технологичности.
Тонкостенные конструкции оболочечного типа весьма чувствительны к локальным нагрузкам. Исследования несущей способности таких конструкций при локальных силовых воздействиях имеют большое практическое значение. Условия эксплуатации при этом характеризуются внешними воздействиями, которые часто приводят к тому, что материал начинает работать за пределами упругости. Учет нелинейных факторов позволяет более адекватно смоделировать процессы деформирования конструкций. Реальные условия эксплуатации требуют решения задач в упругопластической постановке, позволяющей определить истинный запас прочности и использовать этот резерв для решения проблемы снижения материалоемкости конструкций. В ряде случаев в сложных условиях эксплуатации конструкций возникает необходимость оценки напряженно-деформированного состояния (НДС), прочности и жесткости конструкции работающих за пределом упругости при значительных перемещениях.
Нелинейные задачи деформирования тонких оболочек принадлежат к числу наиболее сложных из-за геометрически и физически нелиней-
ности задач. Результаты получены для класса задач, ограниченного по характеру распределения нагрузки, граничным условиям, свойствам материала.
Исследован большой класс симметричных задач упругопластического изгиба пластин и оболочек [1], а также апробирован численный алгоритм для класса несимметричных задач упругопластического деформирования, сочетающих простоту реализации с достаточно быстрой сходимостью, позволяющей использовать их в инженерной практике.
1. ПОСТАНОВКА ЗАДАЧИ
Рассматриваются прямоугольные в плане размера 2а и 2Ь неоднородные однослойные и многослойные пластины и оболочки переменной толщины и кривизны: толщина к(х,у), кривизны к , к2, модуль упругости Е, коэффициент Пуассона и, предел текучести <хг являются некоторыми функциями координат (рис. 1).
Математическая модель решения дважды нелинейных задач деформирования гибких неоднородных оболочек использует пятимерное девиаторное пространство с последующим итерационным процессом.
Рассматриваются неоднородные многослойные оболочки переменной толщины и кривиз-
ны, удовлетворяющие условиям текучести Ми-зеса в каждом из слоев модели.
Процесс нагружения реализован компонентами в девиаторных пространствах А. А. Ильюшина: пространства напряжений 7, деформаций Е", деформаций срединной поверхности 7т,с „ „
Е , изменений кривизны срединной поверхности х.
7ц + 7 99 + 7
22
>33
Рис. 1. Пологая оболочка на прямоугольном плане
Исходные уравнения для тонкостенных конструкций оболочечного типа с учетом физической и геометрической нелинейности записываются на основе деформационной теории пластичности и теории гибких пластин и оболочек в предположении малости сдвигов в нормальных сечениях. Применяется гипотеза прямых нормалей Кирхгофа-Лява с учетом деформаций по толщине.
Компоненты деформаций срединной поверхности панелей и кривизны выражаются нелинейно через перемещения соответственно вдоль координатных осей соотношениями по теории Т. Кармана [3].
2. ОСНОВНЫЕ ФИЗИЧЕСКИЕ СООТНОШЕНИЯ
Выбираем в качестве основной неизвестной функции прогиб оболочки w = w( х, у) . В орто-
нормированном базисе (ё, ё,, ё3) напряжения представляются в виде:
7 = 7^^! + 72 ^2 + 73^3
73 = >/2^.
Эц + Э22 I ,
7 (711 722) +(722 733)
7 —(7^ _7з) + 6*(7^ + 7,^ +70
42
2
12 1 ^1^ 23 / ’
7 = — ^7' + 7"
где - компоненты девиатора напряжений, 7ij, 70 - компоненты тензора напряжений и среднее напряжение, 5^- - символ Кронекера, 7М - интенсивность напряжений.
Процесс нагружения в пространстве напряжений однозначно определяется компонентами девиатора напряжений , средним напряжением 70 .
Аналогичные соотношения можно выписать для вектора деформаций.
Е' = Её + Её + Её,
, Е3 "\/2е12 •
=7ij ~5ij70,
Связь между векторами напряжений и деформаций имеет вид:
7 = N • Е’ + д .
Соотношения теории упругости получаются при N = 20, д = 0, соотношения теории малых упругопластических деформаций N = 20^,
д = 0; О = Е /(2(1 + и)) - модуль сдвига, Е -модуль упругости, О =7 /(3^м) - секущий модуль к диаграмме деформирования.
Деформации в эквидистантных слоях в произвольной точке оболочки записываются в виде:
*11 = *п + 2 •Xll,
*22 = *22 + 2 • Х22,
*12 = *12 + 2 • Х12'
Таким образом, деформация оболочки полностью определяется совокупностью шести
функций *, *2С2, *12, Х11, Х22, Х12 • Напряжения, действующие в нормальном сечении, заменяются статически эквивалентной системой усилий и моментов, приложенных к срединнойповерхно-сти.
Расчет упругопластических деформаций в оболочках и пластинах выполняется с использованием последовательных приближений. Усилия и моменты определяются интегрированием
3
напряжений 7 по толщине. Поскольку интегрирование с переменными параметрами упругости затруднительно, то применяется итерационный путь.
3. МЕТОДИКА ИССЛЕДОВАНИЯ
Решение дважды нелинейной краевой задачи проводится на основе сочетания метода переменных параметров упругости и двухступенчатого метода, разработанного автором. Двухступенчатый итерационный метод содержит в основе замену исходного оператора сложной структуры более простым с последующим итерационным процессом.
Дискретизация нелинейных операторов осуществлена методом конечных разностей. Схема вычислений построена таким образом, чтобы максимально использовать экстраполяцию решения для получения нулевого приближения не только по параметру нагрузки, но и по другим параметрам задачи (прогиб и др.). Для получения разрешающей системы уравнений равновесия гибких оболочек в перемещениях и граничных условий используется вариационный принцип Лагранжа для квазистатического процесса деформирования.
Оболочка находится под действием нормальной распределенной по поверхности нагрузки р = р(х, у), по кромкам оболочки реализуются граничные условия жесткого защемления и шарнирного закрепления в любом сочетании.
Разработанный алгоритм позволяет решать несимметричные задачи упругопластического деформирования оболочек при любых сочетаниях граничных условий. Ограничением является лишь то, что вдоль каждой кромки может быть реализовано только одно из них. Физическая нелинейность учитывается по методу переменных параметров упругости, а геометрически нелинейные задачи решаются методом прямого поиска.
Исследуется напряженно-деформированное состояние гибких пластин и оболочек из неоднородного материала. В качестве параметра прослеживания равновесных состояний выбирается параметр нагрузки. Расчет гибких упругопластических слоистых оболочек и пластин основан на гипотезе недеформируемой нормали для всего пакета в целом [1].
Слои обладают соизмеримой жесткостью на сдвиг и деформируются совместно, без проскальзывания (рис. 2).
схема 1
схема 2
схема 3 1, 2 - однослойные 3 - двуслойные 4, 5 - трехслойные
схема 4 □ .
стальной
слой
схема 5 □ дюрал.слой
Рис. 2. Взаимное расположение слоев в многослойных оболочках
Численные результаты представлены для однослойных, двухслойных и трехслойных пластин и оболочек, состоящих из стальных и дюралюминиевых слоев. В качестве материала слоев принималась сталь 30ХГСА с характеристиками:
Е = 194,2гн/м2,О = 74,5гн/м2, /и = 0,3 0к = 14,5 гн / м2 ,стто = 0,95 гн / м2, е = 4,28 • 10~3 гн / м2
то ’
и сплав Д16 с характеристиками:
Е = 74,5гн/м2,0 = 28,6гн/м2,
Л = 0,3 , 0к = 1,55 гн/м2,
<ут = 0,314гн/м2, = 3,65 • 10_3 гн/м2.
Разработанный пакет программ позволяет единообразно проводить расчет гибких однослойных и слоистых пластин и оболочек с учетом упругопластических деформаций, прослеживать развитие зон пластичности, разгрузки, вторичных пластических деформаций. Многослойные оболочки могут иметь симметричную и несимметричную структуру относительно срединной поверхности. В геометрически нелинейной и линейной постановках проведен сравнительный анализ упругопластического поведения однослойных и слоистых пластин и оболочек в зависимости от граничных условий, характера нагружения.
Эпюры для разных типов пластин отличаются количественно и качественно.
Так, наибольший прогиб возникает в пластине с утонченной средней частью. Однако большая интенсивность напряжений возникает в месте изменения толщин в пластине с утолщенной средней частью. Отметим резкое изменение интенсивности напряжений в местах из-
менения толщин. При этом независимо от изменения толщины, напряжения уменьшаются в месте изменения толщин к центру, видно влияние геометрической нелинейности. Большие величины интенсивности напряжений возникают в верхних слоях из-за разномодульности материала, так как модуль упругости на сжатие выше модуля упругости на растяжение.
Рис. 3. Пластина под действием локальной нагрузки: график нагрузка-прогиб и нагрузка-интенсивность напряжений в центре пластины
4. АНАЛИЗ РЕЗУЛЬТАТОВ РАСЧЕТА
Анализ результатов расчета показывает существенное влияние на НДС многослойных пластин геометрической нелинейности и разно-модульности материала. Отметим, что более жесткой и напряженной является пластина с утолщенной средней частью.
Учет пластических деформаций приводит к существенному снижению уровня напряжений 7 на нижней поверхности. Изменение расположения слоев со схемы 4 на схему 5 приводит к перераспределению напряжений. Так, для оболочек с внешним расположением стальных слоев максимальное значение напряжения <7 на нижней поверхности превышает максимальное значение напряжения ох более чем в 2 раза.
Замена внешних слоев на стальные (схема 4) приводит к уменьшению прогиба при нагрузке р = 100 на 30 % по сравнению с панелью из алюминиевого сплава. На внешних поверхностях панели 2 = ±1 происходит перераспределение напряжений по сравнению с однослойной, и максимальное напряжение в трехслойной панели (схема 4) увеличивается в 1,4 раза.
Отметим ряд особенностей, обнаруженных при сравнении результатов расчетов однослойных и двухслойных цилиндрических оболочек с шарнирно закрепленными и заделанными краями. В отличие от шарнирно закрепленных панелей напряжения в панелях с жестко заделанными краями менее чувствительны к изменению характеристик слоев.
Несмотря на различную комбинацию материалов по слоям оболочки, напряжения <7Х в
основном претерпевают изменения лишь в области края, причем для оболочки из сплава Д16 влияние жесткого защемления выражено слабо. В центральной зоне оболочки напряжения 7
мало зависят от комбинации слоев, а влияние жесткого защемления тем значительнее, чем большую толщину имеет стальной слой.
На рис. 4 показано распределение зон пластичности в сечении 2 = —1 для жестко заделанной по всему контуру цилиндрической панели к = 0, к = 20, где к, к2 - кривизны панели при различных вариантах нагружения р = 50, р = 75, р = 100.
Рис. 4. Распределение зон: зона пластичности *, зона разгрузки +
Интенсивность нагрузки выбиралась таким образом, чтобы главный вектор нагрузки был постоянен для всех схем нагружения и соответствовал равномерно распределенному по всей поверхности давлению интенсивности. Зоны пластичности появляются при нагрузке р = 50, с увеличением р развиваются более интенсивно, занимают почти всю поверхность и уже при нагрузке р = 100 появляются зоны разгрузки. Локальное нагружение приводит к существенному изменению характера деформирования пластин и оболочек.
Зоны пластичности и зоны разгрузки в двухслойных панелях с разными толщинами слоев появляются практически при одних и тех же нагрузках, влияние стального слоя проявляется в том, что зоны разгрузки захватывают большую площадь нижнего слоя. При увеличении нагрузки в центральной части нижней лицевой поверхности однослойной оболочки из дюралевого сплава появляется зона вторичных пластических деформаций, а вблизи краев образуется зона разгрузки. В двухслойных оболочках зона вторичных пластических деформаций захватывает большую часть зоны разгрузки, по-
лученной на предыдущей стадии нагружения, а рядом с ней появляются новые точки разгрузки.
Построенная математическая модель учитывает сжимаемость материала и реальный вид диаграммы деформирования. Вычислительные эксперименты показывают, что предложенный метод и алгоритм позволяют учитывать плавное изменение толщины тонкостенных элементов и могут быть применены при их оптимальном проектировании.
Анализ результатов расчетов показывает, что построенный алгоритм позволяет учитывать дискретное подкрепление ребрами и учитывать поперечную нагрузку в виде полосовой.
Рис. 5. Моделирование процесса деформирования гибкой оболочки
Установлено, что неоднородность свойств материала по толщине оболочки может приводить к качественному изменению распределения напряжений. Моделирование задач деформирования нелинейных оболочек показали хорошую сходимость двухступенчатого итерационного метода, разработанного автором, при расчете упругопластического деформирования гибких неоднородных оболочек. Пример моделирования жестко заделанной по всему контуру цилиндрической панели, локально нагруженной, представлен на рис. 5.
Схема вычислений построена таким образом, чтобы максимально использовать экстраполяцию решения для получения нулевого приближения не только по параметру нагрузки, но и по другим параметрам задачи (прогиб и др.).
ВЫВОДЫ
Установлено, что неоднородность свойств материала по толщине оболочки может приводить к качественному изменению распределения напряжений. Результаты получены с учетом геометрической и физической нелинейности.
При увеличении нагрузки в центральной части нижней лицевой поверхности однослойной оболочки из дюралевого сплава появляется
зона вторичных пластических деформаций, а вблизи краев образуется зона разгрузки.
В двухслойных оболочках зона вторичных пластических деформаций захватывает большую часть зоны разгрузки, полученной на предыдущей стадии нагружения, а рядом с ней появляются новые точки разгрузки. Выявлены зоны активного нагружения, разгрузки, вторичных пластических деформаций по объему оболочки, проведен анализ влияния характера распределения нагрузки, граничных условий, распределения толщин на упругопластическое деформирование пластин и оболочек. Исследованы траектории напряжений и деформаций в характерных точках пластин и оболочек. Модель учитывает сжимаемость материала, реальный вид диаграммы и позволяет исследовать деформирование гибких неоднородных оболочек переменной жесткости при поперечном, продольном и комбинированном нагружении.
Реализация предложенного метода и алгоритма компьютерного численного анализа напряженно-деформированного состояния конструкции осуществлена в виде пакета прикладных программ. Реализован экспорт из программного комплекса значений перемещений для расчета напряжений и деформаций программой МКЭ ANSYS. Приводится сравнительный анализ результатов, полученных обоими способами.
Программный продукт предназначен для применения в отраслевых САПР и ЕЯР-системах, допускает автономное использование.
СПИСОК ЛИТЕРАТУРЫ
1. Старожилова О. В. Метод моделирования нелинейных задач деформирования тонких оболочек // Проблемы нелинейной механики деформируемого твердого тела: 2-я Междунар. конф. (8-11 дек. 2009, Казань): тр. конф. Казань: Казанск. гос. ун-т, 2009. С. 367-369.
2. Старожилова О. В. Двухступенчатый итерационный метод в прикладных задачах упруго-пластичности // Обозрение прикладной и промышленной математики. 2002. Т. 9, № 3. С. 655.
3. Малинин Н. Н. Прикладная теория пластичности и ползучести. М.: Машиностроение, 1975. 398 с.
ОБ АВТОРЕ
СТАРОЖИЛОВА Ольга Владимировна, доц. каф. высш. математики, докторант Московск. авиац. ин-та (нац. иссл. ун-та). Дипл. математик-преп. (Самарск. гос. ун-т, 1984). Канд. техн. наук по мех. деформ. тв. тела (СамГТУ,1997). Иссл. в обл. упр.-пласт. деформ. гибких неоднородных оболочек.
METADATA
Title: Investigation of the stress-strain state of thin-walled structural elements.
Authors: O. V. Starozhilova
Affiliation: Povolzhskiy State University of Telecommunications and Informatics (PSUTI), Russia.
Email: [email protected].
Language: Russian.
Source:Vestnik UGATU (scientific journal of Ufa State Aviation Technical University), vol. 17, no. 4 (57), pp. 38-43, 2013. ISSN 2225-2789 (Online), ISSN 1992-6502 (Print).
Abstract: The article is devoted to numerical modeling of thin-walled structural elements in the form of plates and shells, research and the building of nonlinear mathematical models of stress-strain state. The model takes into account compressibility of the material, the real appearance of the chart and allow to investigate the deformation of a flexible inhomogeneous shells of variable rigidity with a transverse, longitudinal and combined loading.
Key words: plate; shell; nonlinear deformation; variable rigidity.
References (English Transliteration):
1. O. V. Starozhilova, "Method of modeling of nonlinear deformation of thin shells," (in Russian), Problems of Nonlinear Mechanics of Deformable Solids: Proc. 2nd Int. Conf. Kazan, 2009. Kazan State Univ., Russia, 2009. pp. 367-369.
2. O. V. Starozhilova, "Two-step iterative method in applied problems elasto-plasticity," (in Russian), Obozrenie Prikladnoy I Promyshlennoy Matematiki (Surveys of Applied and Industrial Mathematics), vol. 9, no. 3, p. 655, 2002.
3. N. Malinin, Applied Theory of Plasticity and Creep, (in Russian). Moscow: Mashinostroyenie, 1975.
About author:
STAROZHILOVA, Olga Vladimirovna, associate professor of the Higher Mathematics Dept., doctoral candidate of the Moscow Aviation Institute (National Research University). Dipl. mathematician (Samara State Univ., 1984). Cand. (PhD) Tech. Sci. (SamGTU, 1997).
