Construction of high-precision low-dimensional MgFE using local approximations and generating FE

Full Text

Введение

В современной авиационной и ракетно-космической технике широко применяются упругие конструкции, пластины, балки, оболочки, имеющие неоднородную (волокнистую) структуру. При анализе напряженно-деформированного состояния (НДС) упругих конструкций (тел) активно используется метод конечных элементов (МКЭ) [1–7]. Базовые дискретные модели (БМ) тел, которые учитывают их неоднородную, микронеоднородную структуру в рамках микроподхода [8] с помощью конечных элементов (КЭ) 1-го порядка, имеют очень высокую размерность.

Для понижения размерностей дискретных моделей очень результативно используется метод многосеточных конечных элементов (ММКЭ) [9–11]. Особенно эффективно ММКЭ используется при решении задач теории упругости [12–16] для композитных тел (КТ). ММКЭ реализуется на основе функционала Лагранжа (т. е. в перемещениях) [17] с применением многосеточных конечных элементов (МнКЭ) [18–21]. Поскольку при построении n-сеточного КЭ используется не одна, а  вложенных сеток (n ≥ 2), то ММКЭ можно считать обобщением МКЭ, т. е. МКЭ – частный случай ММКЭ. Отсюда следует, что если в расчетах тел по МКЭ применяются МнКЭ, то в этом случае, по сути, реализуется ММКЭ. Неоднородные, микронеоднородные структуры КТ учитываются в рамках микроподхода с помощью МнКЭ. Существующие подходы и численные методы решения задач упругости для КТ имеют сложные формулировки и труднореализуемы [22–30]. Для тел с неоднородной регулярной структурой с различными коэффициентами наполнения, которые с позиций макроподхода рассматриваются как изотропные однородные тела, НДС определяется с помощью фиктивных модулей упругости. Процедуры нахождения фиктивных модулей упругости показаны в работе [31] для двумерных композитов, в [32] – для трехмерных. Основные достоинства ММКЭ состоят в том, что он порождает для КТ дискретные модели малой размерности и реализуется на базе процедур МКЭ.

При построении по известным процедурам многосеточного конечного элемента (кратко – стандартного МнКЭ) используются мелкая сетка и крупные, вложенные в мелкую. Мелкая (базовая) сетка порождена базовым разбиением, которое учитывает неоднородную, микронеоднородную  структуру стандартного МнКЭ, крупные сетки применяются для понижения его размерности. Для стандартных МнКЭ характерно следующее. Всякая крупная сетка стандартного МнКЭ и отвечающие ей аппроксимации перемещений определяются на всей его области. Это приводит к увеличению размерности стандартного МнКЭ при повышении его порядка точности, так как в этом случае на крупных сетках используются аппроксимаци перемещений высокого порядка. Для уменьшения погрешности приближенных решений используются (стандартные) высокоточные МнКЭ, т. е. МнКЭ высокого порядка точности, которые имеют большую размерность. Однако применение высокоточных МнКЭ затруднительно, так как они порождают дискретные модели тел высокой размерности.

Основной недостаток МКЭ и ММКЭ состоит в том, что при построении по МКЭ, ММКЭ решений с малой погрешностью для КТ с помощью известных односеточных КЭ, МнКЭ необходимо применять дискретные модели тел высокой размерности. Для решения данной проблемы здесь предлагается использовать высокоточные МнКЭ малой размерности, которые позволяют для крупных дискретных моделей КТ определять напряжения с малой погрешностью.

В данной работе предлагается метод локальных аппроксимаций (МЛА) для построения высокоточных МнКЭ малой размерности (кратко – малоразмерных МнКЭ), которые применяются при расчете НДС по ММКЭ упругих однородных и КТ. МЛА применяется для КТ, которые удовлетворяют следующим основным положениям.

Положение 1. КТ состоят из разномодульных изотропных однородных тел, связи между которыми идеальны, т. е. на общих границах разномодульных изотропных однородных тел функции перемещений и напряжений являются непрерывными.

Положение 2. Перемещения, деформации и напряжения разномодульных изотропных однородных тел отвечают соотношениям Коши и закону Гука трехмерной задачи линейной теории упругости [16].

Следует отметить, что МЛА, по сути, порождает новые подходы построения МнКЭ. Рассмотрены два подхода, т. е. два типа малоразмерных МнКЭ. Малоразмерные МнКЭ 1-го типа проектируются на базе стандартных с применением локальных аппроксимаций перемещений, которые определяются на подобластях стандартных МнКЭ, 2-го типа – с применением образующих КЭ. Краткая суть построения малоразмерных МнКЭ 1-го типа состоит в следующем. Согласно МЛА в центральной части области V₀ стандартного МнКЭ применяются локальные аппроксимации перемещений высокого порядка (построенные на крупных сетках), в окрестности границы области V₀ – малого порядка, что позволяет проектировать высокоточный МнКЭ V_p малой размерности, т. е. малоразмерный МнКЭ 1-го типа, размерность которого равна размерности стандартного. Следует отметить следующее: При увеличении порядка точности МнКЭ V_p размерность его не меняется, т. е. не увеличивается, и поэтому он называется высокоточным МнКЭ малой размерности. Подробно изложена процедура построения малоразмерных МнКЭ 1-го типа. Как известно, расчет на статическую прочность упругих конструкций сводится к нахождению для них максимальных эквивалентных напряжений, определение которых с малой погрешностью для КТ является в настоящее время актуальной проблемой.

1. Расчеты показывают, что малоразмерные МнКЭ 1-го типа порождают в КТ максимальные эквивалентные напряжения, погрешности которых в 25–50 раз меньше погрешностей аналогичных напряжений, полученных с применением стандартных МнКЭ, которые имеют такие же формы, характерные размеры, размерности, неоднородную структуру и образуют такие же дискретные модели, как малоразмерные МнКЭ. Приведенный пример расчета КТ по ММКЭ с применением малоразмерных МнКЭ 1-го типа показывает их высокую эффективность.
2. Основные достоинства малоразмерных МнКЭ 1-го типа состоят в следующем. Во-первых,
3. малоразмерные МнКЭ 1-го типа порождают в КТ максимальные эквивалентные напряжения с меньшей погрешностью, чем стандартные, и образуют дискретные модели КТ малой размерности. Во-вторых, применение в расчетах по ММКЭ малоразмерных МнКЭ 1-го типа позволяет для крупных дискретных моделей КТ определять максимальные эквивалентные напряжения с малой погрешностью. Отсюда следует, что малоразмерные МнКЭ 1-го типа более эффективны, чем стандартные (на базе которых построены малоразмерные).
4. Кратко показана процедура построения малоразмерных МнКЭ 2-го типа на базе стандартных высокоточных МнКЭ с применением образующих КЭ. Малоразмерные МнКЭ 2-го типа имеют такой же порядок точности, как стандартные, но размерность их меньше размерности стандартных. Поэтому малоразмерные МнКЭ 2-го типа порождают дискретные модели КТ меньшей размерности, чем стандартные.

1. Процедура построения стандартных многосеточных конечных элементов

Основные положения построения стандартных МнКЭ на основе функционала Лагранжа, не теряя общности суждений, покажем на примере лагранжевого двухсеточного конечного элемента  (2сКЭ) $V_d^{\left(2\right)}$ размерами 8h × 8h × 8h (рис. 1), h – задано. Здесь и далее МнКЭ на рисунках показаны в локальной декартовой системе координат Oxyz. Рассматриваемый 2сКЭ $V_d^{\left(2\right)}$ армирован непрерывными волокнами сечением h × h, которые параллельны оси Oy, расстояние между которыми равно h. На рис. 2 показано сечение 2сКЭ в плоскости , сечения волокон закрашены. Считаем, что между компонентами неоднородной структуры 2сКЭ связи идеальны, а функции перемещений, напряжений и деформаций этих компонентов удовлетворяют закону Гука и соотношениям Коши, отвечающие трехмерной задачи линейной теории упругости [12–16], т. е. во всей области 2сКЭ $V_d^{\left(2\right)}$ реализуется трехмерное НДС [16]. Область 2сКЭ $V_d^{\left(2\right)}$ представляем БМ  R_d, состоящей из однородных односеточных КЭ (1сКЭ) $V_j^h$ первого порядка формы куба со стороной h [6; 7], j = 1,..., M; M – общее число 1сКЭ $V_j^h$, для рис. 1 имеем M = 512. На рис. 1 показана БМ R_d 2сКЭ $V_d^{\left(2\right)}$, которая учитывает неоднородную структуру 2сКЭ $V_d^{\left(2\right)}$ и порождает мелкую (базовую) узловую равномерную сетку h_d размерности 9 × 9 × 9 с шагом h. 

 

Рис. 1. Сетки 2сКЭ ${\mathit{V}}_{\mathbf{d}}^{\mathbf{\left(}\mathbf{2}\mathbf{\right)}}$ Fig. 1. Grids 2gFE ${\mathit{V}}_{\mathbf{d}}^{\mathbf{\left(}\mathbf{2}\mathbf{\right)}}$	Рис. 2. Сечение 2сКЭ ${\mathit{V}}_{\mathbf{d}}^{\mathbf{\left(}\mathbf{2}\mathbf{\right)}}$ Fig. 2. Section 2gFE ${\mathit{V}}_{\mathbf{d}}^{\mathbf{\left(}\mathbf{2}\mathbf{\right)}}$

На мелкой сетке h_d определяем крупную равномерную сетку $H_d^{\left(2\right)}$ размерности 5 × 5 × 5 с шагом 2h. На рис. 1 узлы крупной сетки $H_d^{\left(2\right)}$ отмечены точками – 125 узлов. Полную потенциальную энергию П_d базового разбиения R_d 2сКЭ $V_d^{\left(2\right)}$ (т. е. функционал Лагранжа [17]) представим в матричной форме [6; 7]:

${П}_d=\sum _{j=1}^M(\frac{1}{2}{\mathbf{q}}_j^T\left[K_j^h\right]{\mathbf{q}}_j-{\mathbf{q}}_j^{\mathit{T}}{\mathbf{P}}_j),$                                                                   (1)

где $\left[K_j^h\right]$ – матрица жесткости; P_j q_j– векторы узловых сил и неизвестных 1сКЭ $V_j^h$; T – транспонирование; М – общее число 1сКЭ $V_j^h$.

С помощью полиномов Лагранжа [6] на сетке $H_d^{\left(2\right)}$ определяем аппроксимирующие функции u₂, v₂, w₂ для перемещений u, v, w 2сКЭ $V_d^{\left(2\right)}$, которые запишем в форме

$u_2=\sum _{i=1}^5\sum _{j=1}^5\sum _{k=1}^5{N}_{ijk}{u}_{ijk}$,   $v_2=\sum _{i=1}^5\sum _{j=1}^5\sum _{k=1}^5{N}_{ijk}{v}_{ijk}$,   $w_2=\sum _{i=1}^5\sum _{j=1}^5\sum _{k=1}^5{N}_{ijk}{w}_{ijk}$,         (2)

где $u_{ijk},v_{ijk},w_{ijk}$ – искомые значения функций u₂, v₂, w₂ в узле i, j, k сетки $H_d^{\left(2\right)}$; i, j, k – координаты целочисленной системы координат ijk, которая введена для узлов крупной сетки  (рис. 1); N_ijk = N_ijk (x, y, z) – базисная функция узла i, j, k сетки $H_d^{\left(2\right)}$, i, j, k =1,.....,5, где

$N_{ijk}=L_i\left(x\right)L_j\left(y\right)L_k\left(z\right)$,                                                                             (3)

где $L_i\left(x\right)=\prod _{\alpha =\mathrm{1,}\alpha \ne i}^5\frac{x-x_{\alpha }}{x_i-x_{\alpha }}$;  $L_j\left(y\right)=\prod _{\alpha =\mathrm{1,}\alpha \ne j}^5\frac{y-y_{\alpha }}{y_j-y_{\alpha }}$;  $L_k\left(z\right)=\prod _{\alpha =\mathrm{1,}\alpha \ne k}^5\frac{z-z_{\alpha }}{z_k-z_{\alpha }}$,

x_i, y_j, z_k – координаты узла i, j, k сетки $H_d^{\left(2\right)}$ в системе координат Oxyz.

Обозначим: N_b = N_ijk , u_β = u_ijk , v_β = v_ijk , w_β = w_ijk , где i , j, k = 1,...,5 , т. е. β = 1,...,125.  Тогда выражения (2) принимают вид

$u_2=\sum _{\beta =1}^{125}{N}_{\beta }{u}_{\beta }$,  $v_2=\sum _{\beta =1}^{125}{N}_{\beta }{v}_{\beta }^{}$,  $w_2=\sum _{\beta =1}^{125}{N}_{\beta }{w}_{\beta }^{}$.                                        (4)

Обозначим: q_d = {u₁,...,u₁₂₅, v₁,...,v₁₂₅, w₁,...,w₁₂₅}^T – вектор узловых перемещений крупной сетки $H_d^{\left(2\right)}$, т. е. вектор неизвестных 2сКЭ $V_d^{\left(2\right)}$. Используя (4), компоненты вектора q_j узловых неизвестных 1сКЭ $V_j^h$ выражаем через компоненты вектора q_d, в результате получим

${\mathbf{q}}_j=\mathit{\left[}{A}_{\mathit{j}}^{\mathit{d}}\right]{\mathbf{\text{d}}}_d$,                                                                                               (5)

где $\left[A_j^d\right]$ – прямоугольная матрица; j = 1,.....M.

Подставляя (5) в выражение (1), из условия ∂П_d /∂q_d = 0  получаем матричное соотношение вида [K_d ] q_d = F_d , где

$\text{[}{K}_d]=\sum _{j=1}^M{\left[A_j^d\right]}^T\left[K_j^h\right]\left[A_j^d\right]$, ${\mathbf{F}}_d=\sum _{j=1}^M{\left[A_j^d\right]}^T{\mathbf{P}}_{\mathbf{j}}$ ,                                                     (6)

где [K_d], F_d – матрица жесткости (размерности 375 × 375) и вектор узловых сил (размерности 375) стандартного 2сКЭ $V_d^{\left(2\right)}$.

Особенность стандартных МнКЭ заключается в том, что всякая крупная сетка стандартного МнКЭ и отвечающие ей аппроксимации перемещений определяются на всей его области.

Стандартный 2сКЭ формы куба, имеющий равномерную крупную сетку размерности (n + 1)×(n + 1)×(n + 1), называется 2сКЭ n-го порядка. Поскольку равномерная крупная сетка $H_d^{\left(2\right)}$ 2сКЭ $V_d^{\left(2\right)}$ имеет размерность 5 × 5 × 5 , то 2сКЭ $V_d^{\left(2\right)}$ называется 2сКЭ 4-го порядка. Отметим, что размерность лагранжевого 2сКЭ n-го порядка формы куба при использовании крупной равномерной сетки равна 3(n+1)³, т. е. с увеличением порядка n размерность 2сКЭ резко возрастает.

Замечание 1. Решение, построенное для крупной сетки $H_d^{\left(2\right)}$ 2сКЭ $V_d^{\left(2\right)}$, с помощью формулы (5) проецируется на мелкую сетку h_d базового разбиения R_d 2сКЭ, что дает возможность вычислять напряжения в любом 1сКЭ $V_j^h$ разбиения R_d, следовательно, можно определять напряжения в любом компоненте неоднородной структуры 2сКЭ $V_d^{\left(2\right)}$.

Замечание 2. В силу (5) размерность вектора q_d (размерность 2сКЭ $V_d^{\left(2\right)}$ ) не зависит от числа М, т. е. от размерности разбиения R_d. Следовательно, для учета в 2сКЭ $V_d^{\left(2\right)}$ сложной неоднородной (микронеоднородной) структуры, можно использовать сколь угодно мелкие базовые разбиения R_d, состоящие из 1сКЭ $V_j^h$. В этом случае в 1сКЭ $V_j^h$ сколь угодно точно описывается трехмерное напряженное состояние (без упрощающих гипотез).

Замечание 3. Отметим случай, когда 2сКЭ $V_d^{\left(2\right)}$ имеет сложную форму и его крупная сетка имеет внешние узлы, которые совпадают с узлами крупных сеток соседних с ним 2сКЭ. В этом случае при построении 2сКЭ $V_d^{\left(2\right)}$ во всех узлах его мелкой сетки искомые перемещения u, v, w  выражаются через узловые перемещения крупной сетки $H_d^{\left(2\right)}$ 2сКЭ $V_d^{\left(2\right)}$, кроме тех узлов мелкой сетки, которые совпадают с узлами (со стыковочными узлами) крупных сеток соседних 2сКЭ и сетки $H_d^{\left(2\right)}$, что обеспечивает стыковку 2сКЭ $V_d^{\left(2\right)}$ с соседними с ним 2сКЭ.

Рис. 3. Крупная сетка H_d⁽³⁾

Fig. 3. Large grid H_d⁽³⁾

 

При построении (трехсеточного) 3сКЭ используем конечное число 2сКЭ, их крупные сетки образуют мелкую сетку 3сКЭ, на которой определяем крупную сетку для 3сКЭ. Процедура построения 3сКЭ изложена в работах [9; 19]. Рассмотрим частный случай, когда построение 3сКЭ $V_d^{\left(3\right)}$ проводится на базе одного 2сКЭ $V_d^{\left(3\right)}$ . Крупную сетку $H_d^{\left(2\right)}$ 2сКЭ $V_d^{\left(2\right)}$ считаем мелкой сеткой 3сКЭ. На сетке $H_d^{\left(3\right)}$ для 3сКЭ $V_d^{\left(3\right)}$ определяем крупную сетку $H_d^{\left(3\right)}$. На рис. 3 узлы равномерной

сетки $H_d^{\left(3\right)}$ с шагом 4h отмечены точками – 27 узлов.

Затем по выше описанной процедуре определяем матрицу жесткости и вектор узловых сил стандартного 3сКЭ . Построение стандартного 3сКЭ (МнКЭ) $V_d^{\left(3\right)}$ на базе одного стандартного 2сКЭ целесообразно применять в случае, когда стандартный 2сКЭ имеет большие геометрические характерные размеры.

2. Построение малоразмерных МнКЭ с применением локальных аппроксимаций

Малоразмерные МнКЭ 1-го типа проектируются на базе стандартных с применением локальных аппроксимаций перемещений, которые определяются на подобластях стандартных МнКЭ. Рассмотрим процедуру построения малоразмерных МнКЭ 1-го типа.

Предлагаемая процедура реализуется на базе более мелкой сетки, чем мелкая базовая сетка стандартного МнКЭ. Основные положения рассматриваемой процедуры кратко заключаются в следующем. Согласно МЛА на области  стандартного МнКЭ определяем более мелкую сетку H, чем базовая сетка H₀, H₀ ⊂ H. Область V₀ представляем n тонкостенными граничными и m внутренними областями. Граничные (внутренние) области имеют общую границу (не имеют общей границы) с областью V₀ (общие границы областей не вырождаются в точку). На граничных (внутренних) областях определяем n(m) крупных сеток, которые вложены в мелкую сетку H и порождают аппроксимации перемещений малого (высокого) порядка. На границе области V₀ число узлов крупных сеток мало. На граничных и внутренних областях, используя их мелкие и крупные сетки, строим n граничных и m внутренних 2сКЭ, которые образуют высокоточный р-сеточный КЭ (pсКЭ), где p = 1 + n + m. Отметим, что при построении pсКЭ используются одна мелкая сетка H и n + m (в общем случае различных) крупных сеток. Выражая в pсКЭ перемещения внутренних узлов крупных сеток (с помощью метода конденсации [6]) через перемещения граничных узлов, получаем высокоточный pсКЭ малой размерности, т. е. малоразмерный МнКЭ 1-го типа. Следует отметить, что математические операции метода конденсации являются математическими тождественными преобразованиями, т. е. они не влияют на погрешность решения.

Для простоты изложения, рассмотрим стандартный МнКЭ, мелкая базовая сетка H которого имеет большую размерность. На мелкой сетке H определяем крупную сетку на всей области МнКЭ, которая также имеет высокую размерность. Тогда в центральной части области МнКЭ на крупной сетке определяем аппроксимации перемещений высокого порядка, в окрестности границы области – малого порядка, что позволяет с помощью различных локальных аппроксимаций варьировать размерность и порядок точности малоразмерного МнКЭ. Основные положения процедуры построения малоразмерных МнКЭ 1-го типа показаны в примере 1.

Пример 1. Рассмотрим модельную задачу определения НДС по ММКЭ для тела $V_0^1$ с неоднородной (волокнистой) структурой размерами 16h × 64h × 16h, которое лежит в прямоугольной декартовой системе координат Oxyz (рис. 4), где h – задано. В расчетах тела $V_0^1$  используем его БМ R₀ и дискретные модели R₁ и R₂, состоящие соответственно из стандартных и малоразмерных МнКЭ 1-го типа. Тело $V_0^1$ армировано непрерывными волокнами сечением h × h, параллельными оси Oy, расстояние между волокнами равно h при y = 0: u, v, w = 0, т. е. при y =0 тело жестко закреплено. Для модельной задачи имеем следующие исходные данные:

h = 0,5;  1 E_c = ,  E_v =10 ,  v_c = v_v = 0,3 ,                                          (7)

Рис. 4. Размеры КТ V₀¹

Fig. 4. Sizes CB V₀¹

 

где E_c,  E_v(v_c, v_v) – модули Юнга (коэффициенты Пуассона) соответственно связующего материала и волокна; в точках тела $V_0^1$ с координатами

x_i , y_j , z , где z =16h , 8 ( 1) xi = h_i - , i = 1,2,3 , y_j = 8h_j , j =1,...,8 , действует нагрузка 0,35  →  Fz = 0.35 (рис. 4).

Сечение тела $V_0^1$ в плоскости Oxz показано на рис. 5, сечения волокон закрашены.

БМ R0 КТ  состоит из односеточных конечных элементов (1сКЭ) $V_j^h$ 1-го порядка формы куба со стороной h (в которых реализуется трехмерное НДС [16]). БМ R₀ учитывает в рамках микроподхода неоднородную структуру КТ $V_0^1$ и порождает равномерную (базовую) сетку с шагом h размерности 17 × 65 × 17.

Дискретную модель R₁ КТ $V_0^1$ образуют четыре одинаковых стандартных 2сКЭ $V_q^{\left(2\right)}$ 4-го порядка формы куба, построенных на области $V_q$ размерами 16h × 16h ×16h (рис. 6). На рис. 5, 6 показаны сечение 2сКЭ $V_q^{\left(2\right)}$ и его крупная равномерная сетка с шагом 4h , узлы крупной сетки отмечены точками – 125 узлов. Матрицу жесткости (размерности 375 × 375) и вектор узловых сил (размерности 375) 2сКЭ $V_q^{\left(2\right)}$ определяем по процедуре п. 1.

Рис. 5. Сечение КТ ${\mathit{V}}_{\mathbf{0}}^{\mathbf{1}}$ (МнКЭ ${\mathit{V}}_{\mathbf{p}}^{\mathbf{1}}$)

Fig. 5. Section CB ${\mathit{V}}_{\mathbf{0}}^{\mathbf{1}}$ (MgFE ${\mathit{V}}_{\mathbf{p}}^{\mathbf{1}}$)

 

Рис. 6. Область V_q, 2сКЭ ${\mathit{V}}_{\mathbf{q}}^{\mathbf{\left(}\mathbf{2}\mathbf{\right)}}$, МнКЭ ${\mathit{V}}_{\mathbf{p}}^{\mathbf{1}}$

Fig. 6. Region V_q, 2gFE ${\mathit{V}}_{\mathbf{q}}^{\mathbf{\left(}\mathbf{2}\mathbf{\right)}}$, MgFE ${\mathit{V}}_{\mathbf{p}}^{\mathbf{1}}$

 

Рис. 7. Внутренняя область V₁

Fig. 7. Inner region V₁

 

Дискретная модель R2 КТ $V_0^1$ состоит из 4-х одинаковых малоразмерных МнКЭ $V_0^1$ (1-го типа) размерами 16h × 16h ×16h, сечения которых показаны на рис. 5. Согласно МЛА, построение МнКЭ $V_p^1$ 1-го типа (рис. 6) на базе стандартного 2сКЭ $V_q^{\left(2\right)}$, т. е. на области $V_{\mathrm{q}}$ размерами  16h × 16h ×16h, сводится к следующему. В центральной части области $V_{\mathrm{q}}$ выделяем внутреннюю область $V_1$ размерами 12h × 12h ×12h (рис. 7), которая состоит из восьми одинаковых внутренних 2сКЭ $V_q^{\left(2\right)}$ 3-го порядка размерами 6h × 6h ×6h (рис. 8).

Рис. 8. Сетки внутреннего 2сКЭ

Fig. 8. Grids of the inner 2gFE

 

Рис. 9. Сечение 2сКЭ

 Fig. 9. Section 2gFE

 

Крупная и мелкая сетки и сечение 2сКЭ $V_1^{\left(2\right)}$ показаны на рис. 8, 9, сечения волокон (параллельных оси Oy) размерами h × h закрашены, узлы крупной сетки отмечены точками, 64 узла. Матрицу жесткости (размерности 192 × 192) и вектор узловых сил (размерности 192) 2сКЭ $V_1^{\left(2\right)}$ находим по процедуре п. 1. Внутреннюю область V₁ окружают восемь одинаковых по форме и характерным размерам  граничных тонкостенных областей V₂ сложной формы, толщиной 2h (рис. 10). На области V₂ строим граничный 2сКЭ $V_1^{\left(2\right)}$, используя крупную (равномерную) сетку $H_2^{\left(2\right)}$ размерами 8h × 8h × 8h с шагом 4h, т. е. малой размерности 3 × 3 × 3 (рис. 11), узлы крупной сетки $H_2^{\left(2\right)}$ отмечены точками – 27 узлов. На рис. 10 представлена мелкая равномерная сетка 2сКЭ $V_1^{\left(2\right)}$ с шагом h. Отметим, что 8 узлов крупной сетки $H_2^{\left(2\right)}$ лежат вне области V₂, но совпадают с узлами крупной сетки 2сКЭ $V_1^{\left(2\right)}$.

При построении 2сКЭ $V_1^{\left(2\right)}$ во всех узлах его мелкой сетки искомые перемещения u, v, w, выражаем через узловые перемещения крупной сетки $H_2^{\left(2\right)}$, кроме тех узлов мелкой сетки, которые совпадают с узлами сетки $H_2^{\left(2\right)}$ (27 узлов) и с граничными узлами крупной сетки 2сКЭ $V_1^{\left(2\right)}$ (37 узлов), которые обеспечивают стыковку 2сКЭ $V_1^{\left(2\right)}$ и $V_1^{\left(2\right)}$. На рис. 10 эти стыковочные узлы 2сКЭ $V_1^{\left(2\right)}$ отмечены точками (см. замечание 3 п. 1). Итак, 2сКЭ $V_1^{\left(2\right)}$ имеет 64 узла, в которых определяются перемещения, т. е. 2сКЭ $V_1^{\left(2\right)}$  имеет матрицу жесткости размерности 192 × 192 и вектор узловых сил размерности 192. В области V_q восемь одинаковых внутренних 2сКЭ $V_1^{\left(2\right)}$ и восемь граничных 2сКЭ типа $V_1^{\left(2\right)}$ образуют высокоточный МнКЭ ${\text{R}}_{\text{p}}^{\text{1}}$.

Рис. 10. Граничная область V₂ (2сКЭ${\mathit{V}}_{\mathbf{2}}^{\mathbf{\left(}\mathbf{2}\mathbf{\right)}}$)

Fig. 10. Boundary region V₂ (2gFE${\mathit{V}}_{\mathbf{2}}^{\mathbf{\left(}\mathbf{2}\mathbf{\right)}}$)

 

Рис. 11. Сетка ${\mathit{H}}_{\mathbf{2}}^{\mathbf{\left(}\mathbf{2}\mathbf{\right)}}$ 2сКЭ ${\mathit{V}}_{\mathbf{2}}^{\mathbf{\left(}\mathbf{2}\mathbf{\right)}}$

Fig. 11. Grid ${\mathit{H}}_{\mathbf{2}}^{\mathbf{\left(}\mathbf{2}\mathbf{\right)}}$ of the 2gFE ${\mathit{V}}_{\mathbf{2}}^{\mathbf{\left(}\mathbf{2}\mathbf{\right)}}$

 

Выражая в МнКЭ ${\text{R}}_{\text{p}}^{\text{1}}$ перемещения внутренних узлов (с помощью метода конденсации [6]) через перемещения граничных узлов МнКЭ ${\text{R}}_{\text{p}}^{\text{1}}$ , получаем малоразмерный МнКЭ ${\mathit{\text{V}}}_{\text{p}}^{\text{1}}$ 1-го типа. Отметим, что МнКЭ ${\mathit{\text{V}}}_{\text{p}}^{\text{1}}$ 1-го типа имеет такую же размерность, как стандартный МнКЭ ${\mathit{\text{V}}}_{\text{q}}^{\text{(2)}}$ (рис. 6), в котором с помощью метода конденсации исключены перемещения внутренних узлов.

Результаты расчетов КТ ${\mathit{\text{V}}}_{\text{0}}^{\text{1}}$ по ММКЭ с применением дискретных моделей R₀, R₁ и R₂  и  даны в таблице, где σ_n – максимальное эквивалентное напряжение модели R_n (напряжение σ_n определяется по 4-й теории прочности [33]), $N_n^o$ и  $b_n^o$– размерность и ширина ленты системы уравнений ММКЭ модели R_n, n = 0,1,2. Считаем, что БМ R₀ КТ ${\mathit{\text{V}}}_{\text{0}}^{\text{1}}$ порождает точное решение, т. е. напряжение σ_o отвечает точному решению. Тогда относительная погрешность δ_n(%) для напряжения σ_n, отвечающего модели R_n, где n = 1.2,  определяется по формуле

${\delta }_n(\%)=100\%\times |{\sigma }_0-{\sigma }_n|/{\sigma }_0$.                                                     (8)

Анализ результатов таблице показывает, что погрешность напряжения σ₂, которое отвечает модели R₂, состоящей из малоразмерных МнКЭ ${\mathit{\text{V}}}_{\text{0}}^{\text{1}}$ 1-го типа, в k₁ = δ₁(%) /δ₂(%) = 25,122 раз меньше погрешности напряжения σ₁, отвечающего модели R₁, состоящей из стандартных 2сКЭ ${\mathit{\text{V}}}_{\text{0}}^{\text{1}}$ . Для приближенных решений используются скорректированные условия прочности, которые учитывают погрешность напряжений и представлены в следующей теореме.

Результаты расчетов моделей R₀, R₁ и R₂

n	Rn	σn	δn(%)	Nn	bno
0	R0	4,999	–	55488	924
1	R1	4,374	12,511	1200	375
2	R2	5,0246	0,498	1200	375

 

Теорема. Пусть для коэффициента запаса n₀ упругого тела V₀ заданы условия прочности

$n_1\le n_0\le n_2$,                                                                   (9)

где n₁, n₂ – заданы, n₁> 1, n_o = σ_T/σ₀, σ_T– предельное напряжение тела V₀, σ₀ – максимальное эквивалентное напряжение тела V₀, которое отвечает точному решению задачи теории упругости, построенному для тела V₀.

Пусть коэффициент запаса n_b тела V₀, отвечающий приближенному решению задачи теории упругости, удовлетворяет скорректированным условиям прочности вида

$\frac{n_1}{1-{\delta }_{\alpha }}\le n_b\le \frac{n_2}{1+{\delta }_{\alpha }}$.                                                      (10)

Тогда коэффициент запаса n₀ тела V₀, отвечающий точному решению задачи теории упругости, удовлетворяет заданным условиям прочности (9), где , n_b = σ_T/σ_b, σ_b – максимальное эквивалентное напряжение тела V₀, отвечающее приближенному решению задачи теории упругости, построенному для тела V₀, и найденное с такой погрешностью δ_b, что

$\left|{\delta }_b\right|\le {\delta }_{\alpha }<C_{\alpha }=\frac{n_2-n_1}{n_1+n_2}<1$,                                           (11)

где δ_α – верхняя оценка относительной погрешности δ_b, δ_α – задано, погрешность δ_b для напряжения σ_b определяется по формуле δ_b = (σ₀ - σ_b)/σ₀, 0 < δ_α <1.

Теорема доказана в работе [34].

Следует отметить, что с увеличением погрешности δ_b, т. е. с увеличением оценки δ_α, диапазон для коэффициента запаса n_b в скорректированных условиях прочности (10) уменьшается, становится меньше диапазона заданных условий прочности (9). Например, согласно (11), при δ_α = Cα диапазон для скорректированных условий прочности вырождается в точку, что трудно выполнить на практике. Поэтому для практики важно применять приближенные решения с малой погрешностью. Напряжение σ₂, отвечающее модели R₂, отличается от напряжения σ₀ базовой модели R₀ (которое считаем точным) на0,498 % (см. таблицу). Для малых значений погрешности максимальных эквивалентных напряжений, меньше одного процента, оценку δ_α можно принять δ_α(%) = 1%, т. е. δ_α = 0,01. При δ_α =0,01 диапазон Δ₁ скорректированных условий прочности (10) мало отличается от диапазона Δ₂ заданных условий прочности (9), т. е.

${\mathrm{\Delta }}_1\approx {\mathrm{\Delta }}_2, \mathrm{где} {\Delta }_1=n_2/(1+{\delta }_{\alpha })-n_1/(1-{\delta }_{\alpha }), {\Delta }_2=n_2-n_1$.

Дискретная модель R₂ (и модель R₁) требует в k₂ =(N₀× b₀)/(N₁⁰ ×b₁⁰) = 113,94 раз меньше объема памяти ЭВМ, т. е. почти в 114 раз меньше, чем БМ R₀ КТ $V_0^1$ (см. таблицу). Анализ результатов расчетов показывает, что предлагаемые малоразмерные МнКЭ $V_p^1$  (1-го типа) более эффективны, чем стандартные 2сКЭ $V_q^{\left(2\right)}$, которые имеют такие же форму, характерные размеры, размерность и такую же неоднородную структуру, как малоразмерные.

Важно отметить следующее. Согласно замечанию 2 (см. п. 1), при построении МнКЭ можно использовать сколь угодно мелкую базовую сетку. Тогда крупные сетки малоразмерного МнКЭ (1-го типа) могут иметь сколь угодно высокую размерность и, следовательно, порождать на внутренних областях МнКЭ локальные аппроксимации перемещений сколь угодно высокого порядка и число внутренних областей в МнКЭ возрастает. При этом порядок локальных аппроксимаций перемещений на граничных областях и их число не меняется. Это приводит к увеличению порядка точности малоразмерных МнКЭ при постоянстве их размерности. Однако порядок точности малоразмерных МнКЭ не может быть сколь угодно большим, так как реализация метода конденсации, связанная с матрицами высокого порядка, в этом случае требует большого объема памяти ЭВМ, который ограничен.

3. Построение малоразмерных МнКЭ с применением образующих КЭ

Основные положения процедуры построения малоразмерных МнКЭ 2-го типа с применением образующих КЭ показаны в следующем примере.

Пример 2. Основные положения построения малоразмерных МнКЭ 2-го типа с применением образующих КЭ покажем на примере высокоточного лагранжевого 2сКЭ $V_3^{\left(2\right)}$ размерами 9h×10h × 9h  (рис. 12), где h – задано, который расположен в декартовой системе координат Oxyz. Рассматриваемый 2сКЭ $V_3^{\left(2\right)}$ армирован непрерывными волокнами сечением h × h, которые параллельны оси Oy. На рис. 13 показано сечение 2сКЭ $V_3^{\left(2\right)}$, сечения волокон закрашены, расстояние между волокнами равно h. Считаем, что между компонентами неоднородной структуры 2сКЭ связи идеальны, а функции перемещений, напряжений и деформаций этих компонентов удовлетворяют закону Гука и соотношениям Коши, отвечающие трехмерной задачи линейной теории упругости [16], т. е. в области 2сКЭ $V_3^{\left(2\right)}$ реализуется трехмерное НДС.

 

Рис. 12. Сетки 2сКЭ ${\mathit{V}}_{\mathbf{3}}^{\mathbf{\left(}\mathbf{2}\mathbf{\right)}}$ 

Fig. 12. Grids of the 2gFE ${\mathit{V}}_{\mathbf{3}}^{\mathbf{\left(}\mathbf{2}\mathbf{\right)}}$

 

Область 2сКЭ $V_3^{\left(2\right)}$ представляем БМ $R_d^3$, состоящей из однородных 1сКЭ $V_j^h$ первого порядка формы куба со стороной h [6; 7], в которых реализуется трехмерное НДС [16]. На рис. 12 показано базовое разбиение $R_d^3$ 2сКЭ $V_3^{\left(2\right)}$, которое учитывает неоднородную структуру 2сКЭ $V_3^{\left(2\right)}$ и порождает мелкую (базовую) узловую равномерную сетку $h_d^3$ размерности 10 × 11 × 10 с шагом h. На мелкой сетке $h_d^3$ определяем крупную равномерную сетку $H_d^3$ размерности 4 × 6 × 4 с шагом 3h по осям Ox, Oz и с шагом 2h по оси Oy. На рис. 12 узлы крупной сетки $H_d^3$ отмечены точками – 96 узлов. В данном случае с помощью полиномов Лагранжа [6] на крупной сетке $H_d^3$ аппроксимирующие функции u₂, v₂, w₂ для перемещений u, v, w  2сКЭ $V_3^{\left(2\right)}$ представим в форме

 

Рис. 13. Сечение 2сКЭ ${\mathit{V}}_{\mathbf{3}}^{\mathbf{\left(}\mathbf{2}\mathbf{\right)}}$

Fig. 13. Section of the 2gFE  ${\mathit{V}}_{\mathbf{3}}^{\mathbf{\left(}\mathbf{2}\mathbf{\right)}}$

 

$\begin{array}{c}{u}_2=\sum _{i=1}^4\sum _{j=1}^6\sum _{k=1}^4{N}_{ijk}{u}_{ijk},\\ v_2=\sum _{i=1}^4\sum _{j=1}^6\sum _{k=1}^4{N}_{ijk}{v}_{ijk},\\ w_2=\sum _{i=1}^4\sum _{j=1}^6\sum _{k=1}^4{N}_{ijk}{w}_{ijk}.\end{array}$                                         (12)

где $u_{ijk},v_{ijk},w_{ijk}$ – искомые значения функций $u_2,v_2,w_2$ в узле $i,j,k$ крупной сетки $H_d^3$; $i,j,k$ – координаты целочисленной системы координат $ijk$, которая введена для узлов крупной сетки $H_d^3$ (рис. 12); $N_{ijk}$ – базисная функция узла $i,j,k$ сетки $H_d^3$, i, k = 1,... 4, j = 1,..., 6, где

$N_{ijk}=L_i\left(x\right)L_j\left(y\right)L_k\left(z\right)$,                                                        (13)

где $L_i\left(x\right)=\prod _{\alpha =\mathrm{1,}\alpha \ne i}^4\frac{x-x_{\alpha }}{x_i-x_{\alpha }}$, $L_j\left(y\right)=\prod _{\alpha =\mathrm{1,}\alpha \ne j}^6\frac{y-y_{\alpha }}{y_j-y_{\alpha }}$, $L_k\left(z\right)=\prod _{\alpha =\mathrm{1,}\alpha \ne k}^4\frac{z-z_{\alpha }}{z_k-z_{\alpha }}$;

$x_i,y_j,z_k$ – координаты узла $i,j,k$ крупной сетки $H_d^3$ в системе координат Oxyz.

Введем обозначения: $N_{\beta }=N_{ijk}$, $u_{\beta }=u_{ijk}$, $v_{\beta }=v_{ijk}$, $w_{\beta }=w_{ijk}$,где i, k = 1,..., 4, j = 1 т. е. имеем β = 1,.... 96. Тогда выражения (12) принимают вид

$u_2=\sum _{\beta =1}^{96}{N}_{\beta }{u}_{\beta }$,    $v_2=\sum _{\beta =1}^{96}{N}_{\beta }{v}_{\beta }^{}$,   $w_2=\sum _{\beta =1}^{96}{N}_{\beta }{w}_{\beta }^{}$.             (14)

Используя (12)–(14), по процедуре п. 1 для 2сКЭ $V_3^{\left(2\right)}$ определяем матрицу жесткости (размерности 288 × 288) и вектор узловых сил (размерности 288), т. е. размерность 2сКЭ $V_3^{\left(2\right)}$ равна 288. Для понижения размерности 2сКЭ $V_3^{\left(2\right)}$ применяем метод конденсации [6]. Выражаем перемещения внутренних узлов крупной сетки  через  перемещения узлов сетки $H_d^3$, которые лежат на границе 2сКЭ $V_3^{\left(2\right)}$. В результате получаем 2сКЭ $V_3^{\left(2\right)}$, который имеет 240 узловых неизвестных. Итак, размерность 2сКЭ $V_3^{\left(2\right)}$ равна 240. Отметим, что 2сКЭ $V_3^{\left(2\right)}$ является высокоточным, так как он имеет 3-й порядок аппроксимации перемещений по осям Ox, Oz и 5-й порядок – по оси Oy.

Кратко рассмотрим процедуру построения малоразмерного 2сКЭ $V_3^{\left(2\right)}$ 2-го типа на основе высокоточного стандартного лагранжевого 2сКЭ $V_3^{\left(2\right)}$ (рис. 12), т. е. высокого порядка точности, с применением образующего КЭ V_L размерами 9h × 9h. На рис. 14 показаны мелкая сетка, отвечающая 2сКЭ $V_3^{\left(2\right)}$, и сетка КЭ V_L, узлы которой отмечены точками – 12 узлов.

Согласно методу образующих КЭ [21], область 2сКЭ $V_4^{\left(2\right)}$ получается путем параллельного перемещения образующего КЭ V_L   квадратной формы вдоль оси Oy на заданное расстояние d = 10h (рис. 15). Сетка образующего КЭ V_L образует крупную сетку  2сКЭ $V_4^{\left(2\right)}$. Отметим, что неоднородные структуры 2сКЭ $V_3^{\left(2\right)}$ и $V_4^{\left(2\right)}$ одинаковы. Поэтому БМ $R_d^4$ 2сКЭ $V_4^{\left(2\right)}$, как и БМ 2сКЭ $V_3^{\left(2\right)}$ , состоит из однородных 1сКЭ $V_j^h$ 1-го порядка формы куба со стороной h, где , где j = 1,....M – общее число 1сКЭ $V_j^h$, M = 810. БМ $R_d^4$ учитывает неоднородную структуру и сложную форму 2сКЭ $V_4^{\left(2\right)}$ и порождает мелкую сетку $h_d^4$.

Рис. 14. Сетки образующего КЭ V_L

Fig. 14. Grids of the generating FE V_L

 

Согласно методу образующих КЭ, общее число узлов N_d крупной сетки $H_d^4$ 2сКЭ $V_4^{\left(2\right)}$, которая вложена в мелкую сетку $h_d^4$, определяется по формуле N_d = n_xzn_y, где n_xz – число узлов образующего КЭ V_L, для которого n_xz = 12, n_y – число узлов крупной сетки $H_d^4$, лежащих на оси Oy (на оси j), 2сКЭ $V_4^{\left(2\right)}$ имеет 6 узлов на оси Oy, т. е. n_y = 6 , тогда N_d = 12 × 6 = 72. Узлы крупной сетки $H_d^4$ равномерно расположены по оси Oy с шагом 2h. На рис. 15 узлы сетки $H_d^4$ отмечены точками – 72 узла. Для 2сКЭ $V_4^{\left(2\right)}$ введены две локальные системы координат: декартовая Oxyz и для узлов крупной сетки $H_d^4$ – целочисленная $ijk$, где $i,j,k=\mathrm{1,...,6}$  (рис. 14, 15). В области 1сКЭ $V_j^h$ НДС описывается уравнениями трехмерной задачи линейной теории упругости [16]. Следовательно, в области 2сКЭ $V_4^{\left(2\right)}$ реализуется трехмерное НДС.

Рис. 15. Сетки 2сКЭ ${\mathit{V}}_{\mathbf{4}}^{\mathbf{\left(}\mathbf{2}\mathbf{\right)}}$

Fig. 15. Grids of the 2gFE ${\mathit{V}}_{\mathbf{4}}^{\mathbf{\left(}\mathbf{2}\mathbf{\right)}}$

 

Поскольку сетка образующего КЭ V_L имеет 12 узлов (построенная на мелкой сетке 2сКЭ $V_3^{\left(2\right)}$, т. е. на мелкой сетке образующего КЭ V_L), которые на рис. 14 отмечены точками, то для описания перемещений в образующем КЭ V_L используем полином P(x, z) 3-го порядка, который в локальной декартовой системе координат Oxz (рис. 14) имеет вид [6]:

$\begin{array}{c}P(x,z)=a_1+a_{2}x+a_{3}z+a_{4}xz+a_5{x}^2+a_6{z}^2+\\ + a_7{x}^{2}z+a_{8}x{z}^2++a_{9}x{z}^3+a_{10}{x}^{3}z+a_{11}{x}^3+a_{12}{z}^3\end{array}$,                    (15)

где a_i – постоянные, i = 1,...., 12.

Базисную функцию $N_{ijk}$ для узла $i,j,k$ крупной сетки $H_d^4$ согласно методу образующих КЭ [20; 21] определяем в виде

$N_{ijk}(x,y,z)=N_{ik}(x,z)L_j\left(y\right)$,                                                   (16)

где $N_{ik}$ – базисная функция узла i, k образующего КЭ V_L, отвечающая полиному P(x, z) вида (15), i, k = 1,....6, $L_j\left(y\right)$  – полином Лагранжа 5-го порядка:

$L_j\left(y\right)=\prod _{p=\mathrm{1,}p\ne j}^6\frac{y-y_p}{y_j-y_p}$,                                                        (17)

где y_p – координата р-го узла сетки $H_d^4$, лежащего на оси J, j // J; , J = $\overline{\mathrm{1,6}}$ (рис. 15).

Используя (15) – (17), на крупной сетке $H_d^4$ аппроксимирующие функции u₂, v₂, w₂ для перемещений u, v, w  2сКЭ $V_4^{\left(2\right)}$ формально представим

$\begin{array}{c}{u}_2=\sum _{i=1}^6\sum _{j=1}^6\sum _{k=1}^6{N}_{ijk}{u}_{ijk},\\ v_2=\sum _{i=1}^6\sum _{j=1}^6\sum _{k=1}^6{N}_{ijk}{v}_{ijk},\\ w_2=\sum _{i=1}^6\sum _{j=1}^6\sum _{k=1}^6{N}_{ijk}{w}_{ijk},\end{array}$                                                      (18)

где $u_{ijk},v_{ijk},w_{ijk}$ – искомые значения функций $u_2,v_2,w_2$ в узле $i,j, k$ крупной сетки $H_d^4$; $i,j, k$ – координаты целочисленной системы координат ijk, которая введена для узлов крупной сетки $H_d^4$ (рис. 15); $N_{ijk}$ – базисная функция узла $i,j,k$ сетки $H_d^4$, i, k, j = 1,....,6.

Обозначим $N_{\beta }=N_{ijk}, u_{\beta }=u_{ijk}, v_{\beta }=v_{ijk}, w_{\beta }=w_{ijk},$где $u_{ijk},v_{ijk},w_{ijk}$– искомые перемещения в узле $i,j, k$ сетки $H_d^4$; i, k = 1,6; i = 3,4 , k = 2,5 ; i = 2,5 , k = 3,4 ; j = 1,...,6, рис. 14, 15, т. е. β = 1,....,75. Тогда функции перемещений (18) для сетки $H_d^4$ принимают вид

$\begin{array}{c}{u}_2=\sum _{\beta =1}^{72}{N}_{\beta }{u}_{\beta },\\ v_2=\sum _{\beta =1}^{72}{N}_{\beta }{v}_{\beta },\\ w_2=\sum _{\beta =1}^{72}{N}_{\beta }{w}_{\beta },\end{array}$                                                                (19)

где $N_{\beta }$, $u_{\beta }$, $v_{\beta }$, $w_{\beta }$ – базисная функция и перемещения β-го узла сетки $H_d^4$, β = $\overline{\mathrm{1,72}}$.

Используя (15)–(19), по процедуре п. 1 для 2сКЭ $V_4^{\left(2\right)}$ определяем матрицу жесткости (размерности 216 × 216) и вектор узловых сил (размерности 216). Для понижения размерности 2сКЭ  применяем метод конденсации [6], т. е. выражаем перемещения внутренних узлов крупной сетки  через перемещения узлов сетки $H_d^4$, которые лежат на границе 2сКЭ $V_4^{\left(2\right)}$.В результате получаем 2сКЭ $V_4^{\left(2\right)}$ (2-го типа), имеющий 120 узловых неизвестных. Итак, размерность 2сКЭ $V_4^{\left(2\right)}$ равна 120. В силу (15)–(17) 2сКЭ $V_4^{\left(2\right)}$ имеет 3-й порядок аппроксимации перемещений по осям Ox, Oz и 5-й порядок – по оси Oy, т. е. является высокоточным. Отметим, что 2сКЭ $V_4^{\left(2\right)}$ имеет такой же порядок аппроксимаций перемещений по осям Ox, Oy, Oz, такие же размеры 9h × 10h × 9h  и неоднородную структуру, как стандартный 2сКЭ $V_3^{\left(2\right)}$, но размерность 2сКЭ $V_4^{\left(2\right)}$, равная 120, в 2 раза меньше размерности 2сКЭ $V_3^{\left(2\right)}$, размерность которого равна 240. Итак, малоразмерные 2сКЭ $V_4^{\left(2\right)}$ 2-го типа порождают дискретные модели КТ меньшей размерности, чем стандартные 2сКЭ $V_3^{\left(2\right)}$. Следует отметить следующее. Если крупная и мелкая сетки стандартного высокоточного лагранжевого МнКЭ имеют большую размерность, то при построении аппроксимирующих функций перемещений для образующего КЭ целесообразно использовать локальные аппроксимации перемещений, построенные на его мелкой сетке. Особенности малоразмерных МнКЭ 1-го и 2-го типов состоят в следующем. Малоразмерные МнКЭ 1-го типа имеют более высокий порядок точности, чем стандартные, что позволяет проектировать дискретные модели КТ малой размерности, которые порождают напряжения с малой погрешностью. Малоразмерные МнКЭ 2-го типа имеют такой же порядок точности, как стандартные высокоточные МнКЭ, но образуют дискретные модели меньшей размерности, чем стандартные.

Заключение

В данной работе предлагается метод локальных аппроксимаций (МЛА) для построения высокоточных МнКЭ малой размерности, кратко – малоразмерных МнКЭ, которые проектируются на базе известных (стандартных) МнКЭ. Рассмотрены малоразмерные МнКЭ двух типов. Построение малоразмерные МнКЭ 1-го типа проводится с применением локальных аппроксимаций перемещений, определяемых на подобластях МнКЭ, 2-го типа – с применением образующих конечных элементов. Расчеты композитных тел (КТ) показывают, что малоразмерные МнКЭ 1-го типа порождают максимальные эквивалентные напряжения, погрешности которых в 25–50 раз меньше погрешностей аналогичных напряжений, полученных с помощью стандартных МнКЭ, которые имеют такие же размерности, формы, размеры и неоднородные структуры, как малоразмерные. Малоразмерные МнКЭ 1-го типа для крупных дискретных моделей КТ порождают максимальные эквивалентные напряжения с малой погрешностью. Малоразмерные МнКЭ 2-го типа образуют дискретные модели КТ меньшей размерности, чем стандартные.

About the authors

Alexander D. Matveev

Institute of Computational Modeling SB RAS

Author for correspondence.
Email: mtv241@mail.ru

Cand. Sc., associate Professor, senior researcher

Russian Federation, 50/44, Akademgorodok, Krasnoyarsk, 660036

References

Supplementary files