ЧИСЛЕННЫЙ МЕТОД ОПРЕДЕЛЕНИЯ СОСТОЯНИЯ НЕОДНОРОДНЫХ ОБОЛОЧЕК Старожилова О.В.

Поволжский государственный университет телекоммуникаций и информатики


Номер: 1-1
Год: 2015
Страницы: 10-13
Журнал: Актуальные проблемы гуманитарных и естественных наук

Ключевые слова

Численное моделирование, упруго - пластическое деформирование оболочек, итерационный метод, вторично пластические деформации, зоны пластичности, разгрузки, распределения толщин многослойных оболочек, Iterative method, thin-walled designs, again plastic deformations, zones of plasticity, unloading, digitization of nonlinear operators

Просмотр статьи

⛔️ (обновите страницу, если статья не отобразилась)

Аннотация к статье

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

Текст научной статьи

Широкое применение в современной технике тонкостенных элементов конструкций вызывает необходимость разработки надежных и эффективных методов их расчета. Для решения нелинейных задач используются различные методы последовательных приближений: метод упругих решений, метод переменных параметров упругости, вариационные методы, методы конечных разностей и конечных элементов. Возникающие в различных отраслях прикладные задачи упругопластичности являются существенно нелинейными. Условия эксплуатации тонкостенных конструкций в виде гибких оболочек требуют решения задач в упруго-пластической постановке, использования достаточно эффективных методов линеаризации нелинейных проблем. Математическая модель решения дважды нелинейных задач деформирования гибких неоднородных оболочек разработана на основе теории малых упруго-пластических деформаций А.А.Ильюшина, использует пятимерное девиаторное пространство с последующим применением итерационного процесса. При знакопеременном нагружении используются соотношения В.В.Москвитина. Физические уравнения используемой теории позволяют единообразно представить соотношения различных теорий пластичности. В разрешающей системе уравнений учитывается как физическая, так и геометрическая нелинейность, обусловленная тем, что хотя оболочки рассматриваются при малых деформациях, но прогибах одного порядка с толщиной или значительно превышающие ее. Тонкостенные конструкции оболочечного типа весьма чувствительны к локальным нагрузкам. Исследования несущей способности таких конструкций при локальных силовых воздействиях имеют большое практическое значение. Условия эксплуатации при этом характеризуются внешними воздействиями, которые часто приводят к тому, что материал начинает работать за пределами упругости. Учет нелинейных факторов позволяет более адекватно смоделировать процессы деформирования конструкций. Дискретизация нелинейных операторов осуществлена методом конечных разностей. Схема вычислений построена таким образом, чтобы максимально использовать экстраполяцию решения для получения нулевого приближения не только по параметру нагрузки, но и по другим параметрам задачи (прогиб и др.). Разработанный алгоритм реализован в виде пакета программ, позволяющего единообразно решать широкий класс задач исследования упруго-пластических процессов в пологих оболочках. Для решения упруго-пластических задач введена трехмерная пространственная сетка, в каждом узле которой проводится анализ на принадлежность упругой или упругопластической области деформирования. Внутренние силовые факторы определяются интегрированием напряжений по толщине. Для решения дважды нелинейной краевой задачи используется двухступенчатый итерационный метод. С учетом принятой схемы дискретизации по пространственным переменным система pазностных уравнений после исключения неизвестных в контурных и законтуpных точках примет вид где - вектоp неизвестных pазмеpности , - число узлов сеточной области , - квадpатная матpица, - pазностный аналог нелинейной части уравнений, - известный вектоp. Двухступенчатый метод решения системы имеет вид где - тождественный опеpатоp, - опеpатоp сокращения погpешности за итеpаций в методе пеpеменных напpавлений пpи pешении уpавнения . При реализации алгоритма решения дважды нелинейной краевой задачи на каждом шаге по параметру прослеживания равновесных состояний организовано два итерационных цикла. Внешний цикл, где вычисляются невязки, основан на методе переменных параметров упругости. Уточнение секущего модуля идет до тех пор, пока в каждом узле трехмерной сеточной области интенсивности напряжений и деформаций не будут соответствовать диаграмме деформирования. Во внутреннем цикле решаются интегральные уравнения с использованием метода переменных направлений. В каждом узле пространственной трехмерной сетки проводится анализ на принадлежность его упругой или упруго-пластической области деформирования, по найденным деформациям срединной поверхности находятся усилия и моменты. Итерации выполняются до достижения необходимой степени установления процесса, когда характерные значения параметров напряженно-деформированного состояния оболочки перестают меняться в пределах заданной точности. На каждом шаге по параметру нагружения используется метод экстраполяции для получения нулевых приближений. При этом расчет проводится с учетом физической нелинейности. На рис.1 показано распределение зон пластичности в сечении для жестко заделанной по всему контуру цилиндрической панели ,,, где , - кривизны, -толщина панели, - геометрический размер панели, при различных вариантах нагружения ,,. Интенсивность нагрузки выбиралась таким образом, чтобы главный вектор нагрузки был постоянен для всех схем нагружения и соответствовал равномерно распределенному по всей поверхности давлению интенсивности. Зоны пластичности (рис.1) появляются при нагрузке , с увеличением развиваются более интенсивно, занимают почти всю поверхность и уже при нагрузке появляются зоны разгрузки. Локальное нагружение приводит к существенному изменению характера деформирования пластин и оболочек. Проведенный анализ совместного влияния несимметрии в нагрузке и граничных условиях позволяет сделать вывод, что оболочки, находящиеся под действием равномерно распределенной нагрузки, более чувствительны к изменению граничных условий, чем локально нагруженные оболочки. На основе разработанного алгоритма построен вычислительный комплекс, позволяющий единообразно решать широкий класс задач исследования упруго-пластических процессов в гибких пластинах и оболочках. При разработке вычислительного комплекса учитывались следующие основные требования: наибольшая общность постановок задач из рассматриваемого класса, эффективность в отношении точности решения, экономичность по затратам машинного времени и использованию памяти, компактность и удобство в задании исходной информации, возможность расширения решаемого класса задач. В отдельном постпроцессорном модуле реализована функция визуализации полученных результатов решения упруго-пластических задач: строятся распределения напряжений, деформаций, прогибов, зон пластичности, разгрузки, вторичных пластических деформаций по сечениям оболочки, для характерных точек пластин и оболочек приводятся траектории напряжений и деформаций. Переменность толщины, кривизны, неоднородность материала, нагружение оболочек по малой площадке приводит к появлению областей с большими градиентами напряжений и деформаций. В этом случае необходимо использовать густые сетки, что, в свою очередь, приводит к значительному увеличению размерности систем разностных уравнений. Внешний цикл, где вычисляются невязки, основан на методе переменных параметров упругости. Уточнение секущего модуля идет до тех пор, пока в каждом узле трехмерной сеточной области интенсивности напряжений и деформаций не будут соответствовать диаграмме деформирования. Во внутреннем цикле решаются уравнения с использованием метода переменных направлений. В каждом узле пространственной трехмерной сетки проводится анализ на принадлежность его упругой или упруго-пластической области деформирования, по найденным деформациям срединной поверхности находятся усилия и моменты. Итерации выполняются до достижения необходимой степени установления процесса, когда характерные значения параметров напряженно-деформированного состояния оболочки перестают меняться в пределах заданной точности. На каждом шаге по параметру нагружения используется метод экстраполяции для получения нулевых приближений. Схема вычислений построена таким образом, чтобы максимально использовать экстраполяцию решения для получения нулевого приближения не только по параметру нагрузки, но и по другим параметрам задачи (прогиб и др.). В отдельном постпроцессорном модуле реализована функция визуализации полученных результатов: строятся распределения напряжений, прогибов, зон пластичности, разгрузки, вторичных пластических деформаций по сечениям оболочки, для характерных точек приводятся траектории напряжений и деформаций. Построенная математическая модель учитывает сжимаемость материала и реальный вид диаграммы деформирования. Дано решение широкого класса несимметричных задач упруго-пластического изгиба неоднородных оболочек переменной жесткости. Вычислительные эксперименты показывают, что предложенный метод и алгоритм позволяют учитывать плавное изменение толщины тонкостенных элементов и могут быть применены при их оптимальном проектировании. Реализация предложенного метода и алгоритма компьютерного численного анализа напряженно-деформированного состояния конструкции осуществлена в виде пакета прикладных программ. Программный продукт предназначен для применения в отраслевых САПР и ERP-системах, допускает автономное использование.

Научные конференции

 

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