# Abstract

Construction procedure of combined discrete models for the stress state of three-dimensional elastic bodies, complex shape design having an inhomogeneous structure to be analyzed has been suggested. Combined models are composed of first-order homogeneous single grid three dimensional finite elements (FE) of cube shape and inhomogeneous (composite) double grid finite elements (DgFE) of rectangular parallelepiped shape. In the vicinity of the body fixing or complex boundary, single grid FE are used, DgFE are for the rest of the body. Construction procedure of DgFE has been shown. To construct DgFE, two nested grids, both fine and coarse ones are used. Area of the DgFE is a basic (fine) descritization into FE of the first order taking into account its inhomogeneous structure and generating the fine grid. On a basic descritization the total potential energy functional of DgFE (using the approximations constructed on a coarse grid) is given in matrix form projected on a coarse grid. Formulas to calculate the stiffness matrix and the nodal forces vector of DgFE are obtained by minimization condition of the functional obtained by nodal displacement of coarse grid. DgFE characteristics are that the inhomogeneous structure is taken into consideration by using the base fine descritization, discrete models of small dimension are formed and the solutions with a specified error generated. Error of the solution varies with ratio steps of coarse and fine grids of DgFE. Advantages of combined discrete models are that they have a small dimension, take into account a complex shape bodies, inhomogeneous structure and generate solutions with a specified error. The example of calculation has been demonstrated.

# Full Text

Введение Как известно [1], при анализе трехмерных тел неоднородной (композитной) структуры широко используют микро- и макроподходы. В основе макроподхода лежат гипотезы, накладывающие определенные ограничения на поля перемещений, деформаций и напряжений тела, что порождает неустранимую погрешность в решениях. Кроме того, возникают трудности при реализации макроподхода в случае исследования трехмерных неоднородных, композитных тел сложной формы, особенно с малым коэффициентом наполнения неоднородностью. Это связано с тем, что в окрестности границы крепления тела или границы сложной формы гипотезы не выполняются, т.е. в окрестности данных границ тела реализуется трехмерное напряженное состояние. Микроподход дает возможность точно описывать поведение трехмерных композитных (неоднородных) тел сложной формы с любым коэффициентом наполнения неоднородностью. Однако конечноэлементный анализ таких тел с учетом их структуры сводится к построению базовых дискретных моделей очень высокого порядка, что создает проблемы при реализации метода конечных элементов (МКЭ) [2, 3, 4]. Базовые дискретные модели трехмерных тел, имеющих неоднородную структуру, состоят из однородных односеточных трехмерных КЭ первого порядка формы куба [2, 3]. Базовые разбиения учитывают неоднородные структуры и сложную форму тел и при этом порождают решения со сколь угодно малой погрешностью. В данной работе предлагается совместное применение одно- и двухсеточного моделирования трехмерных упругих тел неоднородной структуры и сложной формы. Такое моделирование порождает смешанные дискретные модели, которые состоят из КЭ различной сеточной структуры: односеточных и двухсеточных КЭ. Подобласть тела, которая включает границу сложной формы или крепления тела, представляем базовым (мелким) разбиением, которое состоит из однородных односеточных КЭ первого порядка формы куба и учитывает неоднородную структуру и сложную форму тела, а остальную часть тела покрываем крупным разбиением на ДвКЭ неоднородной (композитной) структуры формы прямоугольного параллелепипеда. Мелкое и крупное разбиения склеиваем с помощью связующих ДвКЭ и в результате получаем смешанную дискретную модель, размерность которой меньше размерности базовой. Для построения ДвКЭ формы прямоугольного параллелепипеда используем две вложенные трехмерные узловые сетки: мелкую и крупную. Мелкая сетка порождена базовым (мелким) разбиением ДвКЭ, которое состоит из однородных односеточных КЭ первого порядка формы куба и учитывает неоднородную структуру данного ДвКЭ. Построение ДвКЭ сводится к исключению параметров МКЭ в узлах мелкой сетки (которые не совпадают с узлами крупной сетки) с помощью аппроксимирующих функций перемещений , построенных на крупной сетке, т.е. узловые перемещения мелкой сетки выражаются через узловые перемещения крупной сетки. Общее число узлов крупной сетки многократно меньше общего числа узлов мелкой, поэтому размерность ДвКЭ многократно меньше размерности его базового разбиения. Существуют два типа ДвКЭ [5]. У ДвКЭ первого типа [5, 6] крупные сетки определяются только на границе его базового разбиения. Крупная сетка ДвКЭ второго типа определяется на всей его области. Изложены процедуры построения ДвКЭ второго типа и связующих ДвКЭ. Недостаток процедуры построения ДвКЭ первого типа состоит в том, что она содержит операцию вычисления обратной матрицы высокой размерности [5, 6]. Приведен пример расчета трехмерного тела сложной формы, имеющего неоднородную структуру. 1. Процедура построения двухсеточных конечных элементов Изложение данной процедуры проведем для трехмерного ДвКЭ (второго типа) неоднородной структуры формы прямоугольного параллелепипеда размерами , для рис. 1 имеем , . Считаем, что между компонентами неоднородной структуры ДвКЭ связи идеальны, а функции перемещений, напряжений и деформаций этих компонентов удовлетворяют закону Гука и соотношениям Коши [7]. Область ДвКЭ представляем базовым разбиением, состоящим из однородных односеточных КЭ первого порядка формы куба со стороной h [2], параметрами МКЭ которых являются значения перемещений в узлах КЭ , – общее число КЭ . На рис. 1 показано базовое (мелкое) разбиение ДвКЭ на КЭ . Базовое разбиение ДвКЭ учитывает его неоднородную (композитную) структуру и порождает мелкую равномерную узловую сетку размерностью с шагом h по осям , , . Для рис. 1 имеем . На мелкой сетке определяем крупную трехмерную узловую сетку размерности с шагами: H1 по оси Ox, H2 по оси Oy и H3 по оси Oz. Сетка вложена в мелкую сетку , при этом имеем: , H2 = k2h, H3 = k3h, где k1, k2, k3 – целые, не менее 2. На рис. 1 узлы крупной сетки отмечены точками, H1 = 3h, H2 = H3 = 2h, n1 = n2 = n3 = 5, k1 = 3, k2 = k3 = 2. Полную потенциальную энергию Пе ДвКЭ представим в матричном виде [8]: (1) где – матрица жесткости; – векторы узловых сил и неизвестных элемента ; – транспонирование. Рис. 1. Мелкая и крупная сетки ДвКЭ С помощью полиномов Лагранжа [8] на крупной сетке определяем аппроксимирующие функции для перемещений u, v, w ДвКЭ, которые соответственно обозначим через и представим в форме , , , (2) где – искомые значения функций в узле (i, j, k) сетки ; i, j, k – координаты целочисленной системы координат ijk, введенной для узлов крупной сетки (см. рис. 1); – базисная функция узла (i, j, k) сетки , , , , , (3) – координаты узла сетки в системе координат Oxyz. Тройке целых чисел i, j, k узла крупной сетки определим целое число и введем следующие обозначения: где . Тогда выражения (2) примут вид . (4) Обозначим через вектор узловых параметров МКЭ крупной сетки , т. е. вектор узловых неизвестных ДвКЭ. Используя (4), компоненты вектора узловых неизвестных КЭ выражаем через компоненты вектора , в результате получим равенство , (5) где – прямоугольная матрица, Подставляя (5) в выражение (1), из условия получаем уравнение , (6) где – матрица жесткости и вектор узловых сил трехмерного ДвКЭ (второго типа) формы прямоугольной призмы, который обозначим через , е – порядковый номер. Замечание 1. Решение, построенное для крупной сетки ДвКЭ, с помощью формулы (5) проецируется на мелкую сетку базового разбиения ДвКЭ, что дает возможность вычислять напряжения в любом КЭ базового разбиения ДвКЭ, следовательно, определять напряжения в любом компоненте неоднородной структуры ДвКЭ. Достоинства двухсеточных конечных элементов · С помощью базового (мелкого) разбиения ДвКЭ учитывается его неоднородная и микронеоднородная структура. · ДвКЭ порождают двухсеточные дискретные модели, общее число узловых неизвестных которых меньше общего числа неизвестных базовых моделей. · ДвКЭ порождают решения, которые отличаются от решений, отвечающих базовым моделям, на заданную величину. · С помощью варьирования соотношений шагов мелкой и крупной вложенных сеток ДвКЭ регулируется погрешность решений, построенных для двухсеточных дискретных моделей. · Напряжения могут быть определены в любом компоненте неоднородной, микронеоднородной структуры ДвКЭ. · Процедура построения ДвКЭ базируется на известных алгоритмах МКЭ и поэтому удобно реализуется на ЭВМ. Реализация МКЭ для двухсеточных дискретных моделей требует меньше ресурсов ЭВМ и временных затрат, чем для базовых моделей. Расчеты показывают [5], что ДвКЭ типа (см. рис. 1) эффективны при анализе трехмерных тел, которые армированы ортогональными регулярными решетками волокон и область которых представляется ДвКЭ формы прямоугольного параллелепипеда. 2. Процедура построения связующих двухсеточных конечных элементов При совместном применении одно- и двухсеточного моделирования тела неоднородной структуры подобласть тела, которая включает границу крепления и границу сложной формы, представляем мелким разбиением, которое состоит из однородных односеточных КЭ первого порядка, а остальную часть тела (область ) – покрываем крупным разбиением на ДвКЭ неоднородной структуры, – область тела. Мелкое разбиение области учитывает неоднородную структуру и сложную форму тела. В результате получаем смешанную дискретную модель тела, состоящую из КЭ различной структуры. В смешанной дискретной модели тела связь мелкого разбиения области с крупным разбиением области осуществляется с помощью связующих ДвКЭ, которые обозначим через ( – порядковый номер КЭ ). Рассмотрим процедуру построения ДвКЭ . Пусть разбиение области состоит из однородных односеточных КЭ первого порядка формы куба со стороной h, т. е. область представляем КЭ базовой модели тела. Связующий ДвКЭ строим на основе ДвКЭ , т. е. ДвКЭ имеет такие же мелкую и крупную сетки, как и ДвКЭ . Пусть связующий ДвКЭ по границе , которая лежит в плоскости yOz (рис. 2), соприкасается с областью . На рис. 2 при показано разбиение тела на ДвКЭ и . Отметим, что на границе сетка разбиения области содержит узлы крупной сетки ДвКЭ . Полную потенциальную энергию ДвКЭ представим в виде (7) где – матрица жесткости базового разбиения ДвКЭ ; – вектор узловых сил и – вектор узловых неизвестных разбиения . Вектор имеет следующую структуру: (8) где – вектор значений перемещений тех узлов мелкой сетки ДвКЭ , которые лежат на границе и не совпадают с узлами крупной сетки ; – вектор значений перемещений остальных узлов сетки , не совпадающих с узлами крупной сетки ; – вектор узловых неизвестных крупной сетки. Используя (4), между векторами и установим связь: (9) где – прямоугольная матрица. Рис. 2. Расчетная схема тела С помощью (8), (9) построим равенство (10) здесь , – булевы матрицы, – вектор узловых неизвестных связующего ДвКЭ , который граничит с областью (см. рис. 2). Применяя (10) в (7), из условия получаем формулы для вычисления матрицы жесткости и вектора узловых сил ДвКЭ . Итак, ДвКЭ связывают мелкое разбиение области с крупным разбиением области . В этом случае параметры МКЭ узлов мелкой сетки границы , не совпадающие с узлами крупной сетки (т. е. компоненты вектора ), не исключаются и, следовательно, являются узловыми неизвестными связующего ДвКЭ . Замечание 2. Обозначим: – вектор функций перемещений базовой дискретной модели трехмерного тела (имеющего неоднородную структуру), которая состоит из однородных односеточных КЭ первого порядка формы куба со стороной h. Пусть , где – точное решение. Пусть , где – вектор функций перемещений смешанной дискретной модели данного тела, т. е. состоящей из однородных односеточных КЭ и двухсеточных КЭ и неоднородной структуры формы прямоугольного параллелепипеда. Тогда имеем , где . Погрешность определяется параметрами базового разбиения тела. Факторы, влияющие на эту погрешность, изучены в теории МКЭ [9, 10]. Расчеты показывают, что погрешность зависит от характерных геометрических размеров подобласти тела, которая представлена односеточными КЭ первого порядка, и от соотношения шагов мелкой и крупной вложенных сеток ДвКЭ . Достоинства смешанных дискретных моделей · Размерности смешанных дискретных моделей трехмерных тел меньше размерностей базовых моделей. · Смешанные дискретные модели учитывают неоднородную структуру и сложную форму трехмерных тел. · Напряжения могут быть определены в любом компоненте неоднородной структуры смешанной дискретной модели трехмерного тела. · Смешанные дискретные модели трехмерных тел порождают решения, отличающиеся от решений базовых дискретных моделей, на заданную величину. · Реализация МКЭ для смешанных дискретных моделей трехмерных тел требует меньше ресурсов ЭВМ и временных затрат, чем для базовых моделей. 3. Результаты численных экспериментов Рассмотрим в декартовой системе координат Oxyz модельную задачу упругости для трехмерного упругого тела неоднородной структуры, которое имеет сложную форму (см. рис. 2). При тело закреплено, т.е. имеем . Граница крепления на рис. 2 показана штриховкой. Тело армировано ортогональной решеткой волокон с поперечным сечением внутри области и на границе. Расстояние между волокнами по оси Ox равно 10h, по осям Oy, Oz – 6h, т.е. волокна с поперечным сечением являются ребрами ДвКЭ и . Максимальное число волокон, параллельных оси Ox, равно 4, параллельных осям Oy, Oz, – 7. Базовая модель тела, состоящая из однородных односеточных КЭ первого порядка формы куба со стороной h [2], учитывает его неоднородную структуру, сложную форму и порождает мелкую сетку . Для узлов сетки введена целочисленная система координат ijk; для рис. 2 имеем , . Тело нагружено силами кг, которые приложены в узлах сетки с целочисленными координатами (i, 25, k), где i =25, 37, 49, 61; k = 1, 9, 17. Модуль Юнга волокон равен кг/см, связующего материала – кг/см, коэффициент Пуассона для всей области тела равен 0,3; h = 0,5 см. Левый торец тела имеет вырезы: вертикальный – размерами и горизонтальный – размерами , которые порождают границу сложной формы. В связи с этим при построении смешанной дискретной модели для данного тела в области (см. рис. 2), т. е. при , используем мелкое (базовое) разбиение, состоящее из однородных односеточных КЭ первого порядка формы куба со стороной h. Остальная часть тела, т.е. при , представлена крупным разбиением на ДвКЭ и связующие ДвКЭ , размерами (см. рис. 1), е = 1,…, 30, . На рис. 2 показано разбиение тела на ДвКЭ и . Результаты расчетов представлены в табл. 1, 2. Анализ результатов показывает, что максимальное значение перемещений смешанной дискретной модели тела отличается от перемещений базовой модели на 1,8 % (в табл. 1 максимальные значения перемещений , выделены жирным шрифтом). Эквивалентные напряжения смешанной дискретной модели тела и – базовой модели вычисляем в центре тяжести КЭ по четвертой теории прочности. Максимальное значение напряжения отличается от напряжения на 0,1 % (в табл. 2 максимальные значения напряжений , выделены жирным шрифтом). Таблица 1 Узловые перемещения (j = 25; см) k \ i 13 25 37 49 61 73 v0, vh 10,055 10,088 23,618 23,037 37,157 36,263 51,212 49,993 65,113 63,584 77,194 75,734 v0 vh 1 9 9,238 9,199 21,397 21,021 34,511 33,857 48,277 47,320 62,040 60,792 74,914 73,547 v0 vh 17 8,177 8,156 19,919 19,537 33,350 32,616 46,893 45,885 60,562 59,263 72,642 71,410 v0 vh Таблица 2 Эквивалентные напряжения (,кг/см) y z \ x 0,5h 1,5h 3,5h 4,5h 6,5h 16,5h σ 10,978 10,958 16,071 16,054 14,423 14,426 13,589 13,602 12,108 12,141 6,413 6,652 σ0 σh 0,5h 0,5h 1,5h 5,543 5,541 1,952 1,949 1,725 1,721 1,561 1,557 1,280 1,276 0,668 0,690 σ0 σh 3,5h 3,872 3,867 2,355 2,351 1,998 1,990 1,793 1,784 1,348 1,338 0,548 0,651 σ0 σh 4,5h 4,811 4,812 2,445 2,442 2,012 2,003 1,903 1,890 1,878 1,855 0,481 0,568 σ0 σh Базовая модель тела содержит 114282 узловых неизвестных, ширина ленты системы уравнений (СУ) МКЭ равна 1956. Смешанная дискретная модель тела имеет 24834 неизвестных, ширина ленты СУ МКЭ равна 3567 и занимает в 2,5 раза меньше объема памяти ЭВМ, чем лента базовой модели. Время реализации МКЭ для смешанной дискретной модели тела в 2,3 раза меньше, чем для базовой модели. Работа выполнена при финансовой поддержке РФФИ (код проекта 11-01-00053).

### Alexsandr Danilovich Matveev

Institute of computational modeling of SB RAS, Krasnoyarsk, Russian Federation

Email: mtv@icm.krasn.ru
50/44, Akademgorodok, 660036, Krasnoyarsk, Russian Federation Ph. D. in Physical and Mathematical Sciences, Associate Professor, Senior Researcher, Federal State Institution of Science ICM SB RAS

# References

1. Фудзии Т., Дзако М. Механика разрушения композиционных материалов. – М.: Мир, 1982. – 232 с.
2. Зенкевич О. Метод конечных элементов в технике. – М.: Мир, 1975. – 541 с.
3. Галлагер Р. Метод конечных элементов. Основы. – М: Мир, 1984. – 430 с.
4. Оден Дж. Конечные элементы в нелинейной механике сплошных сред. – М: Мир, 1976. – 464 с.
5. Матвеев А.Д. Некоторые подходы проектирования упругих многосеточных конечных элементов // Деп. в ВИНИТИ № 2990-В00. Институт вычислительного моделирования СО РАН, Красноярск, 2000. – 30 с.
6. Матвеев А.Д. Многосеточное моделирование композитов нерегулярной структуры с малым коэффициентом наполнения // ПМТФ. – 2004. – № 3. – С. 161–171.
7. Самуль В.И. Основы теории упругости и пластичности. – М.: Высшая школа, 1982. – 264 с.
8. Норри Д., де-Фриз Ж. Введение в метод конечных элементов. – М.: Мир, 1981. – 304 с.
9. Стренг Г., Фикс Дж. Теория метода конечных элементов. – М: Мир, 1977. – 351 с.
10. Ректорис К. Вариационные методы в математической физике и технике. – М.: Мир, 1985. – 591 с.

# Statistics

#### Views

Abstract - 129

PDF (Russian) - 58