Труды сотрудников ИВМ СО РАН

w10=
Найдено документов в текущей БД: 13
   В36
   А654

    Движение бинарной смеси в плоских и цилиндрических областях
[Текст] : монография / В.К. Андреев и др. - Красноярск : Сибирский федеральный университет, 2012. - 187 с. : ил. - Библиогр.: с. 180-187. - ISBN 978-5-7638-2372-1 : Б. ц.
ГРНТИ
ББК В365.5

Аннотация: В монографии представлены результаты исследований конкретных нестационарных движений бинарной смеси с учетом эффекта термодиффузии, возникающих в достаточно длинных плоских и цилиндрических слоях. Рассмотрены свойства инвариантных решений уравнений термодиффузии, когда на границе раздела двух смесей поверхностное натяжение линейно зависит от температуры и концентрации. Для возникающих сопряженных начально-краевых задач получены априорные оценки всех полей, показывающие их экспоненциальную сходимость с ростом времени к стационарным значениям. Приведены результаты численных расчетов поведения скоростей, температур и концентраций в слоях. Дано обобщение решений Остроумова-Бириха на движение смесей в цилиндрической трубе.

Держатели документа:
ИВМ СО РАН : 660036, Красноярск, Академгородок, 50, стр.44

Доп.точки доступа:
Собачкина, Наталья Леонидовна; Sobachkina N.L.; Andreev V.K.
Экземпляры всего: 1
РСФ (1)
Свободны: РСФ (1)
   В19
   З-382

    Об одном классе итерационных схем решения систем линейных уравнений с незнакоопределенной матрицей
[Препринт] : препринт / Ю. Н. Захаров, В. В. Терешкова, Ю. И. Шокин. - Красноярск : ВЦ СО АН СССР, 1990. - 22 с. - (Препринт ВЦ СО АН СССР ; № 14). - Библиогр.: с. 22. - Б. ц.
В надзаг.: Академия наук СССР, Ордена Ленина Сибирское отделение, Вычислительный центр (г. Красноярск)
ББК В192

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

Держатели документа:
ИВМ СО РАН : 660036, Красноярск, Академгородок, 50, стр.44

Доп.точки доступа:
Терешкова, В.В.; Шокин, Юрий Иванович
Экземпляры всего: 1
РСФ (1)
Свободны: РСФ (1)

    Вопросы сходимости конечно-элементных оценок напряженного состояния силовых конструкций с концентраторами напряжений
[Текст] : статья / А. Н. Рогалев, С. В. Доронин // Современные технологии. Системный анализ. Моделирование. - 2015. - № 1. - С. 8-13 . - ISSN 1813-9108
   Перевод заглавия: ISSUES OF CONVERGENCE OF THE FINITE ELEMENT ESTIMATIONS OF STRESS STATE OF LOAD-BEARING STRUCTURES WITH STRESS CONCENTRATORS
Аннотация: Практическая деятельность по решению задач анализа напряженного состояния силовых конструкций с концентраторами напряжений требует развития численных методов моделирования элементов конструкций в области зон высоких градиентов напряжений, обусловленных, в частности, конструктивными концентраторами. Показано, что для геометрически сложных элементов конструкций широко применяемые подходы к реализации метода конечных элементов не гарантируют сходимости численного решения. В статье рассматриваются задачи, связанные с вопросами надежности вычислений с помощью метода конечных элементов. Они зависят как от надежности математической постановки задачи, так и от точности численного решения поставленной задачи. Рассматриваются практические приемы оценки погрешности вычислений. В качестве перспективных путей обеспечения сходимости и достоверности результатов моделирования предлагаются подходы теории технической устойчивости дифференциальных уравнений, специальные вычислительные процедуры, применяемые на стадиях пре- и постпроцессинга.
Practical activity for solving problems of analysis of stress state of load-bearing structures with stress concentrators demands developing numerical methods for simulating high gradient stress zones of structural elements due to structural concentrators. It is shown that widely used approaches to finite element method realization do not make sure convergence for numerical solution for geometrically complex structural elements. The paper is dedicated to problems associated with reliability of calculations by means of finite element method. They depend on validity of mathematical problem definition and adequacy of numerical solution of given problem. Practical techniques for calculating errors assessment are discussed. In the capacity of perspective path to guaranteeing convergence and validity of simulation results approaches of theory for differential equations technical stability and special computational procedures on pre and post processing stages are proposed.

РИНЦ

Держатели документа:
Институт вычислительного моделирования СО РАН
СКТБ «Наука» КНЦ СО РАН

Доп.точки доступа:
Доронин, Сергей Владимирович; Rogalev A.N.
519.63+539.313
А 64

    Анализ конечно-элементных оценок напряженного состояния силовых конструкций с концентраторами напряжений
[Текст] : статья / С. В. Доронин, А. Н. Рогалев, Е. М. Рейзмунт // Современные технологии. Системный анализ. Моделирование. - 2015. - № 2. - С. 26-31 . - ISSN 1813-9108
   Перевод заглавия: ANALYSIS OF THE FINITE ELEMENT ESTIMATIONS OF STRESS STATE OF LOAD-BEARING STRUCTURES WITH STRESS CONCENTRATORS
УДК

Аннотация: Концентрация напряжения в элементах силовых конструкций является одним из основных факторов их прочности, живучести, техногенной безопасности. Несмотря на значительные успехи в области исследования природы и влияния концентрации напряжений на конструкционную прочность, накопленные аналитические, численные, экспериментальные данные не могут охватить все многообразие конструктивных форм и условий нагружения конструкций. Наиболее распространенным в настоящее время является численный подход к анализу концентрации напряжений, реализуемый чаще всего с использованием коммерческих пакетов конечно-элементного анализа. В этом случае принципиальное значение имеет удовлетворение требованиям точности и достоверности результатов моделирования. Опыт численного исследования напряженного состояния в области концентраторов напряжений свидетельствует о многочисленных случаях нарушения этих требований, возникающих вследствие проблем со сходимостью и устойчивостью численных решений. В настоящей работе рассматриваются конечно-элементные решения ряда элементов конструкций с концентраторами напряжений, не удовлетворяющие требованиям сходимости, и анализируются возможные причины возникающих ошибок.
Stress concentration in load-bearing structures elements is one of the main factors of their strength, survivability, safety. In spite of considerable advances of investigations for nature of stress concentration and its influence on structural strength accumulated analytic, numerical, experimental data can’t cover all diversity of structural forms and loading conditions. The numerical approach for analysing stress concentration with the help of commercial finite-element software is the most widely-spread at present. In this case meeting requirements of accuracy and reliability of simulation data is of fundamental importance. The experience of numerical investigation for stress state of structural elements with stress concentrators demonstrates numerous cases of violation of these requirements due to difficulties with convergence and stability of numerical solutions. In this paper finite-element solutions for stress concentrators which are not meeting requirements of convergence are discussed and possible reasons of errors are analysed.

РИНЦ

Держатели документа:
Институт вычислительного моделирования СО РАН
СКТБ «Наука» КНЦ СО РАН

Доп.точки доступа:
Рогалев, Алексей Николаевич; Rogalev A.N.; Рейзмунт, Елена Михайловна; Reizmunt E.M.
539.3
К 65

    Контроль точности решения при анализе напряженно-деформированного состояния высокоответственных технических объектов
[Текст] : статья / А. Н. Рогалев, С. В. Доронин, А. А. Рогалев // Системы. Методы. Технологии. - 2015. - № 3. - С. 32-38 . - ISSN 2077-5415
   Перевод заглавия: Solution accuracy control for analysis of stress-strain state of critical technical objects
УДК

Аннотация: В статье рассматриваются подходы к оценке вычислительной ошибки при решении системы линейных алгебраических уравнений, в качестве матрицы коэффициентов которой рассматривается матрица жесткости конечно-элементной модели технического объекта. Предлагаемый подход предполагает, что уровень вычислительной ошибки определяется структурой и значениями матрицы коэффициентов, и заключается в численном решении системы линейных уравнений с матрицей жесткости и такой специально подобранной правой частью, для которой известно точное решение. Сравнение численного и точного решений позволяет получить оценку вычислительной ошибки, позволяющую судить о приемлемости построенной конечно-элементной модели. Получение указанной оценки является дополнительной процедурой контроля точности численного решения при анализе его сходимости путем последовательного уменьшения шага конечных элементов. Развиваемый подход весьма актуален для конструкций ответственных технических объектов, где цена ошибки при проектных расчетах оказывается неприемлемо высокой. Для реализации предлагаемого подхода организован интерфейс между пакетом конечно-элементного моделирования ANSYS и вычислительным пакетом компьютерной алгебры Wolfram Mathematica. В качестве примера приводится получение оценки вычислительной ошибки численного решения системы линейных алгебраических уравнений с матрицей жесткости силовой конструкции бака высокого давления для перспективных электрореактивных двигателей космических аппаратов. Силовая конструкция представляет собой оболочку давления, подвешенную на системе вантов с регулируемым уровнем натяжения, закрепленных, в свою очередь, на пространственной стержневой системе - силовой структуре корпуса космического аппарата. Для рассматриваемой конструкции найден уровень конечно-элементной дискретизации, обеспечивающий сходимость численного решения.
The paper is devoted to approaches to a problem of numerical error evaluation when solving the system of linear equations. The stiffness matrix of a finite-element model of a technical object is a coefficient matrix of the system of linear equations. The approach proposed supposes that the level of numerical error is determined by a structure and magnitude of coefficient matrix. The approach consists of numerical solving system of linear equations with stiffness matrix and special right-hand member with exact solution known. Comparison of numerical and exact solutions allows evaluating numerical error and making decision on the quality of finite-element model. Evaluation numerical error is a supplementary procedure for checking accuracy of numerical solution within solution convergence analysis by means of cascade reduction mesh spacing. The approach is of great actuality for structures of critical technical objects with great worth of design calculations error. To implement the approach, data interface between the finite-element analysis package ANSYS and computer algebra package Wolfram Mathematica has been created. Evaluated numerical error has been given as an example for numerical solution system of linear equations with stiffness matrix for load-bearing unit of high pressure tank for perspective spacecraft electrojet engines. The load-bearing unit consists of pressure shell suspended by means of cable system with controlled tension. The cable system is attached to spatial bar system - load-bearing frame structure of spacecraft. For the structures considered the level of finite-element discretization has been determined to provide numerical solution convergence.

РИНЦ

Держатели документа:
Институт вычислительного моделирования СО РАН
Институт космических и информационных технологий Сибирского федерального университета
Специальное конструкторско-технологическое бюро «Наука» Красноярского научного центра СО РАН

Доп.точки доступа:
Доронин, С.В.; Doronin S.V.; Рогалев, А.А.; Rogalyov A.A.; Rogalyov A.N.
519.614
О-118

    Ортогонально-степенной метод решения частичной проблемы собственных значений и векторов для симметричной неотрицательно определенной матрицы
[Текст] : статья / И. В. Киреев // Вычислительные методы и программирование: новые вычислительные технологии. - 2016. - Т. 17, № 1. - С. 44-54 . - ISSN 1726-3522
   Перевод заглавия: An orthogonal power method of solving the partial eigenproblem for a symmetric nonnegative definite matrix
УДК

Аннотация: Предложена и обоснована экономичная версия метода сопряженных направлений для построения нетривиального решения однородной системы линейных алгебраических уравнений с вырожденной симметричной неотрицательно определенной квадратной матрицей. Предложено однопараметрическое семейство одношаговых нелинейных итерационных процессов вычисления собственного вектора, отвечающего наибольшему собственному значению симметричной неотрицательно определенной квадратной матрицы. Это семейство включает в себя степенной метод как частный случай. Доказана сходимость возникающих последовательностей векторов к собственному вектору, ассоциированному с наибольшим характеристическим числом матрицы. Предложена двухшаговая процедура ускорения сходимости итераций этих процессов, в основе которой лежит ортогонализация в подпространстве Крылова. Приведены результаты численных экспериментов.
An efficient version of the conjugate direction method to find a nontrivial solution of a homogeneous system of linear algebraic equations with a singular symmetric nonnegative definite square matrix is proposed and substantiated. A one-parameter family of one-step nonlinear iterative processes to determine the eigenvector corresponding to the largest eigenvalue of a symmetric nonnegative definite square matrix is also proposed. This family includes the power method as a special case. The convergence of corresponding vector sequences to the eigenvector associated with the largest eigenvalue of the matrix is proved. A two-step procedure is formulated to accelerate the convergence of iterations for these processes. This procedure is based on the orthogonalization in Krylov subspaces. A number of numerial results are discussed.

РИНЦ

Держатели документа:
Институт вычислительного моделирования СО РАН

Доп.точки доступа:
Kireev I.V.
517.941.1
С117

    Совместное ползущее движение трех вязких жидкостей в плоском слое: априорные оценки и сходимость к стационарному режиму
[Текст] : статья / Виктор Константинович Андреев, Елена Николаевна Черемных // Сибирский журнал индустриальной математики. - 2016. - Т. 19, № 1. - С. 3-17, DOI 10.17377/SIBJIM.2016.19.101 . - ISSN 1560-7518
УДК

Аннотация: Исследовано частично инвариантное решение ранга два и дефекта три уравнения вязкой теплопроводной жидкости. Оно интерпретируется как двумерное движение трех несмешивающихся жидкостей в плоском канале, ограниченном твердыми неподвижными стенками, на которых известно распределение температур. С математической точки зрения возникающая начально-краевая задача является нелинейной и обратной. При некоторых (часто выполняющихся в практических приложениях) предположениях задача заменяется линейной. Для нее получены априорные оценки, найдено стационарное точное решение и доказано, что с ростом времени решение выходит на этот стационарный режим, если со временем температуры на стенках стабилизируются.

РИНЦ

Держатели документа:
Институт вычислительного моделирования СО РАН
Сибирский федеральный университет

Доп.точки доступа:
Черемных, Елена Николаевна; Andreev V.K.
539.3
И118

    Использование критериев обусловленности при численных расчетах напряженного состояния силовых конструкций
[Текст] : статья / А. Н. Рогалев, С. В. Доронин // Системы. Методы. Технологии. - 2016. - № 2. - С. 91-99, DOI 10.18324/2077-5415-2016-2-91-99 . - ISSN 2077-5415
   Перевод заглавия: Conditioning criteria in numerical computation of the stress state of load-bearing structures
УДК

Аннотация: Одним из основных вопросов при конечно-элементном моделировании высокоответственных технических систем является анализ точности полученных результатов. Решение этого вопроса позволяет обосновать надежность характеристик технических объектов. На точность результатов влияют ошибки округления и погрешность приближенных методов линейной алгебры, применяемых в конечно-элементном анализе, а также ошибки, имеющие непосредственное отношение к методу конечных элементов при выборе (построении) сетки конечных элементов, что прямо отражается в свойствах обусловленности матрицы жесткости и величине ошибки решения системы линейных уравнений с матрицей жесткости. Для численных решений, полученных методом конечных элементов, характеристика точности вычислительной ошибки может определяться как величина ошибки решения системы линейных уравнений с матрицей жесткости. Для анализа точности реализован апостериорный анализ ошибок численных вычислений путем двукратного решения системы линейных уравнений с матрицей жесткости, имеющей специально подобранную правую часть. Для полноты анализа результатов конечно-элементного моделирования изучаются числа обусловленности матрицы коэффициентов системы, что позволяет в первом приближении выявлять степень вырожденности матрицы коэффициентов и определять степень чувствительности численных решений к ошибкам. Применение этих процедур позволяет контролировать точность численного решения и анализировать сходимость. Развиваемый подход весьма актуален для конструкций ответственных технических объектов, так как ошибки при проектных расчетах здесь приводят к серьезным последствиям. Оценки вычислительной ошибки численного решения системы линейных алгебраических уравнений с матрицей жесткости применяются в статье для контроля точности в задачах расчета коэффициента концентрации напряжений в круговой пластине с эксцентричным круговым разрезом. Результаты вычислений подтверждают теоретические рассуждения.
One of the main problems in the finite element modeling of highly responsible technical systems is analysis of the accuracy of the results obtained. The solution of this problem allows to prove reliability of technical objects. Rounding errors and errors of approximate methods of linear algebra, used in the finite element analysis, influent the accuracy of the results as well as errors that are directly related to the method of finite elements in the selection (construction) of the finite element grid. This influence is reflected directly in the properties of a stiffness matrix and the size of the error while solving a system of linear equations with the stiffness matrix. For numerical solutions, obtained by finite element method, the characteristics for the accuracy of the computation error can be defined as the value of the error of the solution of linear equations with the stiffness matrix. To analyze the accuracy, a posteriori analysis of numerical computation errors has been realized by solving doubly the systems of linear equations with the stiffness matrix, having specially selected right side. To complete the analysis of the results of the finite element modeling, the condition number of the coefficients matrix of system has been studying which allows to identify the degree of degeneracy of the coefficient matrix at a first approximation and to determine the degree of sensitivity of numerical solutions to the errors. The use of these procedures allows to control the accuracy of the numerical solution and analyze the convergence. The approach proposed is very important for structuring the important technical objects, as errors in the design constructions lead to serious consequences. Estimates of computational error of the numerical solution of a system of linear algebraic equations with the stiffness matrix used in the article for accuracy control in problems of calculation of stress concentration factor in a circular plate with an eccentric circular cutting. The results confirm the theoretical reasoning.

РИНЦ

Держатели документа:
Институт вычислительного моделирования СО РАН
Специальное конструкторско-технологическое бюро «Наука» Института вычислительных технологий СО РАН

Доп.точки доступа:
Доронин, С.В.; Doronin S.V.; Rogalev A.N.
519.6
С 56

    Совместное использование полулагранжева подхода и метода конечных элементов при моделировании течения вязкой несжимаемой жидкости
[Текст] : статья / Е. В. Дементьева, Е. Д. Карепова, В. В. Шайдуров // Решетневские чтения. - 2017. - № 21-2. - С. 16-17 . - ISSN 1990-7702
   Перевод заглавия: A combination of a semi-lagrangian approach and the finite element method for modeling a viscous incompressible fluid flow
УДК

Аннотация: Рассматривается нестационарное течение вязкой несжимаемой жидкости в плоском канале. На границе вытекания ставятся граничные условия «do nothing». Для построения дискретного аналога используется комбинация полулагранжева подхода и метода конечных элементов. Для обеспечения сходимости метода в среднеквадратичной норме используется локальное интегральное балансовое соотношение между соседними временными слоями. Численные эксперименты для модельных задач демонстрируют сходимость предложенного метода.
The two-dimensional time-dependent Navier-Stokes equations are considered for a viscous incompressible fluid in a channel. On the outlet boundary, the modified “do nothing” condition is imposed. To construct a discrete analogue, we use the conforming finite element method in the combination with a semi-Lagrangian approach. To guarantee the energy conservation and the stability in the mean-square norm, we use the discrete analogue of the local integral balance between two neighboring time levels. A numerical experiment shows the convergence of the proposed numerical method.

РИНЦ

Держатели документа:
Институт вычислительного моделирования СО РАН

Доп.точки доступа:
Дементьева, Е.В.; Dementyeva E.V.; Карепова, Е.Д.; Karepova E.D.; Шайдуров, В.В.; Shaidurov V.V.
519.642.2
О-13

    Об аппроксимации интеграла по криволинейному четырехугольнику
[Текст] : статья / А. В. Вяткин, Е. В. Кучунова // Решетневские чтения. - 2017. - № 21-2. - С. 12-13 . - ISSN 1990-7702
   Перевод заглавия: Approximation of integral over quadrangle with curved boundaries
УДК

Аннотация: Представлена асимптотическая оценка точности аппроксимации интеграла по криволинейному четырехугольнику. Представленная оценка может быть использована для теоретического обоснования сходимости полулагранжевого метода.
We present the theoretical estimation of approximation of integral over curvilinear quadrangle. It allows to get a theoretical justification of convergence order for semi-Lagrangian methods.

РИНЦ

Держатели документа:
Институт вычислительного моделирования СО РАН

Доп.точки доступа:
Вяткин, А.В.; Vyatkin A.V.; Кучунова, Е.В.; Kuchunova E.V.

    Оценка точности численного анализа деформированного состояния силовых конструкций технических объектов
[Текст] : статья / Алексей Николаевич Рогалев, Сергей Владимирович Доронин, Владимир Викторович Москвичев // Вычислительные технологии. - 2018. - Т. 23, № 2. - С. 88-101 . - ISSN 1560-7534
   Перевод заглавия: Estimation of the accuracy of numerical analysis for the deformed state of power structures of technical objects
Аннотация: Под силовыми конструкциями понимают технические устройства, составленные из различных частей, воспринимающие комплекс эксплуатационных нагрузок в штатных и аварийных режимах нагружения. При решении прикладных задач исследования напряженно-деформированных состояний силовых конструкций важна оценка степени близости к точному приближенного решения, полученного на вполне определенной сетке конечных элементов с конечной величиной шага сетки. С учетом влияния ошибок округления сходимость метода конечных элементов контролировать сложно: при большом числе конечных элементов решение может расходиться из-за накапливающихся ошибок округления, даже если условия сходимости выполняются. Описанное в статье применение методов обратного анализа ошибок позволяет достаточно точно контролировать точность численных оценок деформированного состояния силовых конструкций, что подтверждают расчеты, выполненные для практических задач.
The solution of applied problems of technogenic safety, survivability, risk and protection is performed for structures which are close to limiting states. These states are characterized by decreasing safety factors down to one. In this case a mistaken estimation for safety factor may cause the situation when the calculated safety factor will be greater than one but the real safety factor will be less than one. Safety factors estimation is performed on the basis of calculation for stress-strain state characteristics. Thus, the issues of accuracy and reliability of determining stresses and deformations are an integral part of the problem of man-made safety. In the numerical analysis of the stress-strain state, the stiffness matrix of the design model is formed, the dimension of which reaches up to tens of millions. A large number of computations for tasks of this dimension is presumably leading to significant rounding errors. Ensuring the grid convergence of results by decreasing the grid spacing is inconsistent with the growth of computational errors due to rounding. For finite element analysis of power structures of technical objects, methods of a posteriori reverse error analysis are proposed that control the effect of rounding errors on the result when solving a solving system of linear algebraic equations. The coefficient matrix of this system is the stiffness matrix of the finite element model. The basic idea is to obtain and solve a system of equations with a known exact solution. Comparison of the results of exact and numerical solutions allows us to estimate the magnitude of the error.

РИНЦ

Держатели документа:
Институт вычислительного моделирования СО РАН
Красноярский филиал Института вычислительных технологий СО РАН

Доп.точки доступа:
Рогалев, Алексей Николаевич; Rogalyov Alexey Nikolaevich; Доронин, Сергей Владимирович; Doronin Sergey Vladimirovich; Москвичев, Владимир Викторович; Moskvichev Vladimir Viktorovich
539.3
M94

    Multigrid finite elements in the calculations of multilayer cylindrical shells
: статья / A. D. Matveev, A. N. Grishanov // Сибирский журнал науки и технологий. - 2018. - Т. 19, № 1. - P27-36 . - ISSN 2587-6066
   Перевод заглавия: Многосеточные конечные элементы в расчетах многослойных цилиндрических оболочек
УДК

Аннотация: An effective numerical method for calculating linearly elastic multilayer cylindrical shells under static loading implemented on the basis of Finite Element Method (FEM) procedures using the multilayer curved Lagrangian multi- grid finite elements (MFE) of the shell type was proposed. Such shells are widely used in rocket-space and aircraft engineering. MFE are developed in local Cartesian coordinate systems based on small (basic) shell partitions that take into account their heterogeneous structure, irregular shape, combined loading and fixing. The stress strained state (SSS) in the MFE was described by the equations of the three-dimensional elasticity problem without using the addi- tional kinematical and static hypotheses, which allow one to use MFE for the shells of various thicknesses to be calcu- lated. The procedure of constructing the Langrage polynomials in local curvilinear coordinate systems used to develop the shell MFE is presented. The displacements in the MFE were approximated by the power and Lagrange polynomials of different orders. When constructing a n -grid finite element (FE), n ≥ 2, n-nested grids were used. The fine grid was generated by the basic partition of the MFE; the other (coarse) grids were used to reduce its dimension. According to the method, the nodes of the coarse MFE grids are located on the common boundaries of the different modular layers of the shell. The proposed law of the expansion in the number of discrete models using MFE with a constant thickness, multiple of the shell thickness, provides a uniform and rapid convergence of approximate solutions, allowing one to frame solutions with a small error. Multigrid discrete models have 10<sup>3</sup>…10<sup>6</sup> times less unknown MFE than the basic ones. The implementation of the MFE for multigrid models requires 10<sup>4</sup>…10<sup>7</sup> times less computer storage space than for the reference models, which allows one using the proposed method to calculate some large shells. An example of calculating a multilayer cylindrical local loading shell of irregular shape was given. In the calculation, three-grid shell - type FE, developed on the basis of the reference models having from 2 million to 3.7 billion of the nodal MFE unknowns were used. To study the approximate solution convergence and error, a well-known numerical method was used.
Предложен эффективный численный метод расчета линейно-упругих многослойных цилиндрических оболо- чек при статическом нагружении с применением многослойных криволинейных лагранжевых многосеточных конечных элементов (МнКЭ) оболочечного типа. Такие оболочки широко используются в ракетно-космической и авиационной технике. МнКЭ проектируются в локальных декартовых системах координат на основе мелких (базовых) разбиений оболочек, которые учитывают их неоднородную структуру, сложную форму, сложное нагружение и закрепление. Напряженное деформированное состояние в МнКЭ описывается уравнениями трехмерной задачи теории упругости без использования дополнительных кинематических и статических гипотез, что позволяет применять МнКЭ для расчета многослойных оболочек различной толщины. Показана процедура построения в локальных криволинейных системах координат полиномов Лагранжа, которые приме- няются при проектировании оболочечных МнКЭ. Перемещения в МнКЭ аппроксимируются степенными и лагранжевыми полиномами различных порядков. При построении n -сеточного конечного элемента (КЭ), n ≥ 2, используют n вложенных сеток. Мелкая сетка порождена базовым разбиением МнКЭ, остальные n - 1 (крупные) сетки применяются для понижения его размерности. В предлагаемом методе узлы крупных сеток МнКЭ расположены на общих границах разномодульных слоев оболочки. Закон измельчения дискретных моде- лей, в которых используются МнКЭ с постоянной толщиной, кратной толщине оболочки, порождает равно- мерную и быструю сходимость приближенных решений, что дает возможность строить решения с малой погрешностью. Многосеточные дискретные модели имеют в 10<sup>3</sup>-10<sup>6</sup> раз меньше узловых неизвестных, чем базовые. Реализация метода конечных элементов (МКЭ) для многосеточных моделей требует в 10<sup>4</sup>-10<sup>7</sup> раз меньше объема памяти ЭВМ, чем для базовых, что позволяет использовать предложенный метод для расчета оболочек больших размеров. В приведенном расчете многослойной цилиндрической оболочки сложной формы, имеющей локальное нагружение, используются оболочечные трехсеточные КЭ, построенные на базовых моде- лях, которые имеют от 2 миллионов до 3,7 миллиарда неизвестных МКЭ. Для анализа сходимости приближен- ных решений используется известный численный метод.

РИНЦ

Держатели документа:
Institute of Computational Modeling
Novosibirsk State Technical University

Доп.точки доступа:
Matveev, A.D.; Матвеев А.Д.; Grishanov, A.N.; Гришанов А.Н.
519.6
E27

    Efficient method of calculating layered conical shells with lagrange multigrid elements use
[Text] : статья / G. I. Rastorguev, A. N. Grishanov, A. D. Matveev // Сибирский журнал науки и технологий. - 2018. - Т. 19, № 3. - P423-431, DOI 10.31772/2587-6066-2018-19-3-423-431 . - ISSN 2587-6066
   Перевод заглавия: Эффективный метод расчета слоистых конических оболочек с применением лагранжевых многосеточных элементов
УДК

Аннотация: The increased requirements for strength calculations of space-rocket and aviation technology designs cause the need for the development and improvement of approximate solutions for elasticity theory tasks with small error algo- rithms. The article considers the numerical method of calculating elastic layered conical shells (LCS) of various thickness under static loading which are widely used in space-rocket technology. The suggested method uses three-dimensional curvilinear Lagrange multigrid finite elements (MGFE). A system of nested grids is used for MGFE constructing. The fine grid is generated by the basic partition that takes into account MGFE heterogeneous structure. The basic partition dimensionality is reduced with the help of large grids which leads to the system of linear algebraic equations of the small dimension finite elements method. Three-dimensional elasticity theory equations use allows to apply MGFE for calculating LCS of any thickness. Displacements in MGFE are approximated by Lagrange polynomials which, in con- trast to power polynomials, gives the opportunity to design big size three-dimensional thin shell elements. Lagrange polynomials nodes coincide in shell thickness with the nodes of MGFE large grids which lie on the shared borders of multi-module layers. The efficiency of the presented method is that the suggested MGFE generate small dimension discrete models that require 10<sup>3</sup>-10<sup>7</sup> times less electronic computing machine (ECM) memory than basic models. The suggested law of dis- crete models grinding generates uniform and fast convergence of numerical solutions which allows to make solutions with the specified (small) error. Examples of LCS calculating (whole ones as well as with holes) under axisymmetric and local loading are given. Comparative analysis of solutions obtained with the help of MGFE, single-grid finite elements and the program com- plex ANSYS has been conducted.
Повышенные требования к прочностным расчетам конструкций ракетно-космической и авиационной техники вызывают необходимость разработки и совершенствования алгоритмов приближенных решений за- дач теории упругости с малой погрешностью. Рассматривается численный метод расчета упругих слоистых конических оболочек (СКО) различной толщины при статическом нагружении, которые широко применяются в ракетно-космической технике. В предлагаемом методе используются трехмерные криволинейные лагранжевые многосеточные конечные элементы (МнКЭ). При построении МнКЭ используется система вложенных сеток. Мелкая сетка порождена базовым разбиением, которое учитывает неоднородную структуру МнКЭ. С помощью крупных сеток пони- жается размерность базового разбиения, что приводит к системе линейных алгебраических уравнений мето- да конечных элементов малой размерности. Использование уравнений трехмерной теории упругости позволя- ет применять МнКЭ для расчета СКО любой толщины. Перемещения в МнКЭ аппроксимируются полиномами Лагранжа, что в отличие от степенных полиномов дает возможность проектировать трехмерные тонкие оболочечные элементы больших размеров. Узлы полиномов Лагранжа по толщине оболочки совпадают с узла- ми крупных сеток МнКЭ, которые лежат на общих границах разномодульных слоев. Эффективность изложенного метода заключается в том, что предлагаемые МнКЭ порождают дискрет- ные модели малой размерности, для которых требуется в 10<sup>3</sup>-10<sup>7</sup> раз меньше объема памяти ЭВМ, чем для базовых моделей. Предложенный закон измельчения дискретных моделей порождает равномерную и быструю сходимость численных решений, что позволяет строить решения с заданной (малой) погрешностью. Приведены примеры расчетов СКО (цельных и с отверстиями) при осесимметричном и локальном нагру- жениях. Выполнен сравнительный анализ решений, полученных с помощью МнКЭ, односеточных конечных элементов и программного комплекса ANSYS

РИНЦ

Держатели документа:
Institute of Computational Modeling SB RAS
Novosibirsk State Technical University

Доп.точки доступа:
Rastorguev, G. I.; Расторгуев Г.И.; Grishanov, A. N.; Гришанов А.Н.; Matveev, A. D.; Матвеев А.Д.