"ГОСТ Р 57700.43-2023. Национальный стандарт Российской Федерации. Численное моделирование полимерных композиционных материалов. Верификация программного обеспечения"
(утв. и введен в действие Приказом Росстандарта от 17.11.2023 N 1426-ст)

1 РАЗРАБОТАН Федеральным государственным бюджетным образовательным учреждением высшего профессионального образования "Московский государственный технический университет имени Н.Э. Баумана" (ФГБОУ ВПО "МГТУ им. Н.Э. Баумана") и Федеральным государственным бюджетным учреждением "Российский институт стандартизации" (ФГБУ "Институт стандартизации") при участии Объединения юридических лиц "Союз производителей композитов" (Союзкомпозит)

2 ВНЕСЕН Техническим комитетом по стандартизации ТК 700 "Математическое моделирование и высокопроизводительные вычислительные технологии"

3 УТВЕРЖДЕН И ВВЕДЕН В ДЕЙСТВИЕ Приказом Федерального агентства по техническому регулированию и метрологии от 17 ноября 2023 г. N 1426-ст

4 ВВЕДЕН ВПЕРВЫЕ

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

Разработка программных комплексов класса CAE для решения задач прогнозирования свойств полимерных композиционных материалов на основе численного 3D-моделирования требует создания набора верификационных задач, которые позволили бы проводить верификацию создаваемых программных комплексов и тем самым обеспечивать высокую точность получаемых результатов при использовании программных комплексов в практической деятельности.

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

Включенные в данный ГОСТ верификационные задачи основаны на теории асимптотического осреднения (homogenization method) периодических структур, которая в настоящее время хорошо математически обоснована в ряде классических трудов отечественных и зарубежных основоположников этого метода (см. [1] - [6]). Аналитические решения верификационных задач, представленные в настоящем стандарте, получены в [3].

В настоящем стандарте применяется терминология и классификации, введенные в ГОСТ Р 57700.42.

Настоящий стандарт устанавливает общие требования к верификации программного обеспечения, предназначенного для численного моделирования свойств полимерных композиционных материалов (ПО МСКМ), на основе метода конечных элементов (МКЭ).

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

Термины и определения, относящиеся к композиционным материалам (КМ) и их свойствам, соответствуют терминам, приведенным в ГОСТ Р 57700.42.

Настоящий стандарт должен применяться совместно с вышеуказанным стандартом.

Данный стандарт применим для численного моделирования эффективных характеристик полимерных композиционных материалов (ПКМ) методом МКЭ, с использованием только моделей линейно-упругих и теплофизических свойств композитов.

Полимерные композиционные материалы представляют собой КМ на основе полимерной матрицы (см. ГОСТ Р 57700.42, ГОСТ 32794, ГОСТ 33742).

В настоящем стандарте использованы нормативные ссылки на следующие стандарты:

ГОСТ 32794 Композиты полимерные. Термины и определения

ГОСТ 33742 Композиты полимерные. Классификация

ГОСТ Р 57700.10 Численное моделирование физических процессов. Определение напряженно-деформированного состояния. Верификация и валидация численных моделей сложных элементов конструкций в упругой области

ГОСТ Р 57700.14 Верификация получаемых сеточными методами численных решений задач механики сплошной среды

ГОСТ Р 57700.42 Численное моделирование композиционных материалов. Общие положения

Примечание - При пользовании настоящим стандартом целесообразно проверить действие ссылочных стандартов в информационной системе общего пользования - на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет или по ежегодному информационному указателю "Национальные стандарты", который опубликован по состоянию на 1 января текущего года, и по выпускам ежемесячного информационного указателя "Национальные стандарты" за текущий год. Если заменен ссылочный стандарт, на который дана недатированная ссылка, то рекомендуется использовать действующую версию этого стандарта с учетом всех внесенных в данную версию изменений. Если заменен ссылочный стандарт, на который дана датированная ссылка, то рекомендуется использовать версию этого стандарта с указанным выше годом утверждения (принятия). Если после утверждения настоящего стандарта в ссылочный стандарт, на который дана датированная ссылка, внесено изменение, затрагивающее положение, на которое дана ссылка, то это положение рекомендуется применять без учета данного изменения. Если ссылочный стандарт отменен без замены, то положение, в котором дана ссылка на него, рекомендуется применять в части, не затрагивающей эту ссылку.

3.1.1

Примечания

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

2 Условие существования большого числа компонентов является существенным: например, материал, образованный соединением двух или трех различных материалов, не рассматривается в настоящем стандарте в качестве композиционного материала.

3.1.2

3.1.3

3.1.4

3.1.5

3.1.6

3.1.7

3.1.8

3.2.1

3.2.2 эффективный тензор модулей упругости КМ: Тензор модулей упругости однородного материала, эквивалентного в некотором смысле исходному КМ.

Примечание - Как правило, эффективный тензор модулей упругости КМ связывает тензоры напряжений и деформации, осредненные по ЯП КМ.

3.2.3 технические упругие константы: Константы модели линейно-упругой среды, которые обычно определяются непосредственно из экспериментов (для КМ они вычисляются): модули упругости, коэффициенты Пуассона, модули сдвига, модуль объемного сжатия и др.

3.2.4 эффективный тензор теплопроводности КМ: Тензор теплопроводности однородного материала, эквивалентного в некотором смысле исходному КМ.

3.2.5 эффективный тензор теплового расширения КМ: Тензор теплового расширения однородного материала, эквивалентного в некотором смысле исходному КМ.

В настоящем стандарте применены следующие сокращения:

ДАКМ - дисперсно-армированный композиционный материал;

КЭ - конечный элемент;

КЭМ - конечно-элементная модель;

ТКМ - тканевый композиционный материал;

1D КМ - однонаправленно-армированный композиционный материал.

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

5.1.2 Задачами верификации ПО для численного моделирования полимерных композиционных материалов являются:

- подтверждение достоверности результатов численного моделирования упругих свойств КМ, получаемых с помощью верифицируемого ПО;

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

5.1.3 Настоящий ГОСТ необходимо применять вместе с ГОСТ Р 57700.42. Все термины и определения, применяемые в настоящем стандарте, относящиеся как к классификации КМ и их свойствам, так и к моделям свойств КМ, сформулированы в ГОСТ Р 57700.42.

5.1.4 При проведении верификации ПО для численного моделирования полимерных композиционных материалов необходимо выполнять требования ГОСТ Р 57700.10 и ГОСТ Р 57700.14.

5.2.1.1 Задача для расчета компонент эффективного тензора модулей упругости и технических упругих констант СКМ.

5.2.1.2 Задача для расчета эффективного тензора модулей упругости и технических упругих констант СВКМ.

5.2.1.3 Задача для расчета компонент эффективного тензора модулей упругости и технических упругих констант ДАКМ.

5.2.1.4 Задача для расчета компонент эффективного тензора модулей упругости и технических упругих констант 1D КМ.

5.2.1.5 Задача для расчета компонент эффективного тензора модулей упругости и технических упругих констант ТКМ.

5.2.2.1 Задача для расчета компонент эффективного тензора теплопроводности СКМ.

5.2.2.2 Задача для расчета компонент эффективного тензора теплопроводности СВКМ.

5.2.2.3 Задача для расчета компонент эффективного тензора теплопроводности тканевого композиционного материала ТКМ.

5.2.2.4 Задача для расчета компонент эффективного тензора теплового расширения СКМ.

5.2.2.5 Задача для расчета компонент эффективного тензора теплового расширения СВКМ.

5.2.2.6 Задача для расчета компонент эффективного тензора теплового расширения ТКМ.

5.3.1 Верификация ПО проводится для оценки адекватности реализации заложенных в ПО математических моделей для расчета упругих и теплофизических свойств КМ, а также для оценки точности конечно-элементных алгоритмов, реализованных в тестируемом ПО.

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

5.3.3 Верификация осуществляется на основе решения тестовых задач с входными данными (рекомендованными значениями характеристик компонентов КМ).

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

5.3.4 Оценка относительного отклонения полученных результатов и результатов, приведенных в тестовых задачах, проводится для каждого значения эффективной характеристики.

5.3.5 Рекомендуемые примерные размеры КЭМ (количество КЭ и количество узлов в КЭМ) для проведения расчетов представлены в каждой тестовой задаче. Размер КЭМ может меняться в зависимости от математического алгоритма решения задачи, применяемого моделирования эффективных характеристик КМ.

Обязательным требованием является применение именно МКЭ для расчета эффективных характеристик КМ.

А.1.1.1 Описание и постановка задачи

Рассматривается СКМ, ЯП которого состоит из четырех различных компонентов (слоев) (см. рисунок А.1).

Каждый слой является линейно-упругим изотропным материалом, технические упругие константы которых приведены в таблице А.1.

ЯП имеет форму куба с длиной ребра, равной 1 (безразмерные координаты). Безразмерная толщина каждого слоя равна 0,25.

Ортонормированная система координат Oe_i, i = 1, 2, 3, ориентирована следующим образом: векторы e₁, e₂ коллинеарны базовой плоскости слоев, а вектор e₃ ортогонален к этой плоскости.

Задача заключается в том, чтобы вычислить:

- компоненты эффективного тензора модулей упругости (см. [1]) композита ,

- эффективные технические упругие константы композита (см. [7]):

E_i - три модуля упругости в направлениях Oe_i,

G_ij - три модуля сдвига в плоскостях Oe_ie_j, i /= j,

- три коэффициента Пуассона, i /= j,

используя для этого верифицируемое ПО, основанное на решении данной задачи с помощью МКЭ.

Для решения задачи следует использовать постановку трехмерной квазистатической задачи линейной теории упругости (см. [9 - 12]) с граничными условиями, которые зависят от применяемой модели вычисления компоненты эффективного тензора модулей упругости и численного метода реализации этого метода.

Для решения задачи должен быть использован МКЭ.

Пример математической модели и постановки математической задачи для вычисления компонент эффективного тензора модулей упругости КМ на основе теории асимптотического осреднения приведен в [13] и [14].

Для СВКМ пример конечно-элементной тетраэдральной сетки размером 206254 элементов показан на рисунке А.2.

А.1.1.2 Результаты

Результаты решения данной задачи, полученные с помощью верифицируемого ПО, сравнивают со следующими вариантами расчетов:

а) результатами аналитического решения задачи, полученными с помощью теории асимптотического осреднения для СКМ [см. [1], формула (1.12)]

где C_ijkl - компоненты тензора модулей упругости компонент композита, представленные в виде разрывной функции координат;

<C_ijkl> - осреднение по ЯП;

- обратная матрица,

эффективные технические упругие константы композита в этом случае вычисляют с помощью обратного тензора (см. [5]) (тензора упругих податливостей) и следующих формул (см. [4, 5]):

Суммирование по i и j отсутствует.

Замечание: индексацию коэффициентов Пуассона вводят таким образом, чтобы удовлетворялись соотношения взаимности (см. [4], [5]).

Суммирование по i и j отсутствует.

б) результатами предельных вариантов расчетов, в которых характеристики всех слоев выбираются одинаковыми (например, могут быть выбраны характеристики слоя 1 или 2 из таблицы А.1), в этом случае:

1) эффективные модули упругости композита совпадают с модулями упругости C_ijkl выбранного слоя композита (например, 1 или 2)

2) эффективные технические упругие константы композита совпадают с техническими упругими константами выбранного слоя

(все значения , , должны быть получены численно, с помощью конечно-элементного решения задачи);

в) результатами численного решения данной задачи с помощью другого ПО, ранее прошедшего верификацию на данных задачах.

Пример сравнения расчетов по варианту 3 приведен в таблице А.2.

А.1.2.1 Описание и постановка задачи

Рассматривается СВКМ, ЯП которого состоит из четырех различных компонентов (слоев) (см. рисунок А.1).

Каждый слой СВКМ представляет собой 1D КМ с эффективными характеристиками и имеет угол армирования (см. рисунок А.3), значения которого указаны в таблице А.3.

1D КМ в каждом слое является одним и тем же однородным линейно-упругим трансверсально-изотропным материалом, отличается только углом армирования, технические упругие константы 1D КМ приведены в таблице А.4.

Единая для всех слоев СВКМ ортонормированная система координат Oe_i, i = 1, 2, 3 ориентирована следующим образом: векторы e₁, e₂ коллинеарны базовой плоскости слоев, а вектор e₃ ортогонален к этой плоскости.

Собственная система координат каждого слоя ориентирована следующим образом: вектор для всех слоев совпадает, векторы направлены под углами к вектору , где s - номер слоя в ЯП.

ЯП имеет форму куба с длиной ребра, равной 1 (безразмерные координаты). Толщина каждого слоя равна 0,25.

Задача заключается в том, чтобы вычислить:

- компоненты эффективного тензора модулей упругости композита СВКМ;

- эффективные технические упругие константы композита СВКМ:

E_i - три модуля упругости в направлениях Oe_i,

G_ij - три модуля сдвига в плоскостях Oe_ie_j, i /= j,

- три коэффициента Пуассона, i /= j,

используя для этого верифицируемое ПО, основанное на решении данной задачи с помощью МКЭ.

Для решения задачи следует использовать постановку трехмерной квазистатической задачи линейной теории упругости (см. [9] - [12]) с граничными условиями, которые зависят от применяемой модели вычисления компонент эффективного тензора модулей упругости и численного метода реализации этого метода.

Для решения задачи должен быть использован МКЭ.

Пример математической модели и постановки математической задачи для вычисления компонент эффективного тензора модулей упругости КМ на основе теории асимптотического осреднения приведен в [13] и [14].

Пример КЭМ размером 206 254 элемента с тетраэдальным конечным элементом показан на рисунке А.2.

А.1.2.2 Результаты

Результаты решения данной задачи с помощью верифицируемого ПО сравнивают со следующими вариантами расчетов:

а) результатами аналитического решения задачи, полученными с помощью теории асимптотического осреднения для СКМ [см. [1], формула (1.12)].

1) эффективные модули упругости СВКМ

где - компоненты тензора модулей упругости s-го слоя, вычисляемые по формуле

где - осреднение по ЯП,

- обратная матрица,

- компоненты тензора модулей упругости 1D КМ в собственной системе координат ,

- матрица поворота системы координат слоя.

2) эффективные технические упругие константы СВКМ, которые вычисляют с помощью обратного тензора

(тензора упругих податливостей) и формул (А.2) - (А.6);

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

Пример сравнения расчетов по методу 1 приведен в таблице А.5.

А.1.3.1 Описание и постановка задачи

Рассматривается ДАКМ, ЯП которого состоит из двух различных компонентов: наполнителя шаровой формы и матрицы (см. рисунок А.4).

Каждый компонент является упругим изотропным материалом, технические упругие константы которых приведены в таблице А.6.

ЯП имеет форму куба с длиной ребра, равной 1 (безразмерные координаты). Концентрация наполнителя равна 0,6.

Ортонормированная система координат Oe_i, i = 1, 2, 3, ориентирована по ребрам ЯП.

Задача заключается в том, чтобы вычислить:

- компоненты эффективного тензора модулей упругости композита ,

- эффективные технические упругие константы композита:

E_i - три модуля упругости в направлениях Oe_i,

G_ij - три модуля сдвига в плоскостях Oe_ie_j, i /= j,

- три коэффициента Пуассона, i /= j,

используя для этого верифицируемое ПО, основанное на МКЭ.

Для решения задачи следует использовать постановку трехмерной квазистатической задачи линейной теории упругости (см. [9] - [12]) с граничными условиями, которые зависят от применяемой модели вычисления компонент эффективного тензора модулей упругости и численной реализации этого метода.

Для решения задачи должен быть использован МКЭ.

Пример математической модели и постановки математической задачи для вычисления компонент эффективного тензора модулей упругости КМ на основе метода асимптотического осреднения приведен в [13] и [14].

Пример КЭМ на основе тетраэдрального конечного элемента размером 150 979 элементов, построенной для 1/8 ЯП ДАКМ, показан на рисунке А.5. Полная ЯП строится отражением 1/8 ЯП относительно трех плоскостей Oe₁e₃, Oe₂e₃ и Oe₁e₂.

А.1.3.2 Результаты

Результаты решения данной задачи с помощью верифицируемого ПО сравнивают со следующими вариантами расчетов:

а) результатами предельных вариантов расчетов, в которых характеристики матрицы и наполнителя выбирают одинаковыми, в этом случае:

1) эффективные модули упругости композита совпадают с модулями упругости C_ijkl компонентов композита (формула (А.7)),

2) эффективные технические упругие константы композита совпадают с техническими упругими константами выбранного компонента [формулы (А.8) и (А.9)];

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

Пример сравнения расчетов по методу 2 приведен в таблице А.7.

А.1.4.1 Описание и постановка задачи

Рассматривается 1D КМ, ЯП которого состоит из двух компонентов: матрицы и волокна цилиндрической формы (см. рисунок А.6).

Каждый компонент является линейно-упругим изотропным материалом, технические упругие константы которых приведены в таблице А.8.

ЯП имеет форму куба с длиной ребра, равной 1 (безразмерные координаты). Диаметр цилиндрических волокон равен 0,74252.

Объемная концентрация волокон варьируется от 0,4 до 0,78 (для варианта 1) и (для варианта 2), при этом диаметр цилиндрических волокон равен 0,74252.

Ортонормированная система координат Oe_i, i = 1, 2, 3 ориентирована по ребрам ЯП. Направление Oe₃ ориентировано по оси цилиндрических волокон.

Задача заключается в том, чтобы вычислить:

- компоненты эффективного тензора модулей упругости композита ,

- эффективные технические упругие константы композита:

E_i - три модуля упругости в направлениях Oe_i,

G_ij - три модуля сдвига в плоскостях Oe_ie_j, i /= j,

- три коэффициента Пуассона, i /= j,

используя для этого верифицируемое ПО, основанное на конечно-элементном методе решения данной задачи.

Для решения задачи следует использовать постановку трехмерной квазистатической задачи линейной теории упругости (см. [9] - [12]) с граничными условиями, которые зависят от применяемой модели вычисления компонент эффективного тензора модулей упругости и численной реализации этого метода.

Для решения задачи должен быть использован МКЭ.

Пример математической модели и постановки математической задачи для вычисления компонент эффективного тензора модулей упругости КМ на основе теории асимптотического осреднения приведен в [13] и [14].

Пример КЭМ размером 17 513 элементов, построенной для 1/4 ЯП, показан на рисунке А.7. Полная ЯП строится отражением 1/4 ЯП относительно плоскостей Oe₁e₃ и Oe₂e₃.

А.1.4.2 Результаты

Результаты решения данной задачи с помощью верифицируемого ПО сравнивают со следующими вариантами расчетов:

а) результатами аналитического решения задачи, полученными с помощью теории асимптотического осреднения и метода функций комплексного переменного для волокнистых композитов (таблица А.9) (см. [1], таблицы 2.2 и 3.2); сравниваются эффективные модули упругости 1D КМ ;

б) результатами предельных вариантов расчетов, в которых характеристики всех компонентов (матрицы и волокна) выбирают одинаковыми, например могут быть выбраны совпадающими с характеристиками наполнителя (или матрицы) из таблицы А.8, в этом случае имеют место формулы (А.7) - (А.9);

в) результатами численного решения данной задачи с помощью другого ПО, ранее прошедшего верификацию на данных задачах.

Пример сравнения расчетов по методу 3 приведен в таблице А.10, характеристики компонентов при этом соответствуют варианту 2 из таблицы А.8.

А.1.5.1 Описание и постановка задачи

Рассматривается ТКМ, ЯП которого состоит из пяти различных компонентов: четырех нитей и матрицы (см. рисунок А.8). Форма нитей близка к синусоидальной, две нити ориентированы в направлении Oe₁, а две другие - в направлении Oe₂. Система нитей образует взаимное переплетение без касания.

Каждый компонент ТКМ является линейно-упругим изотропным материалом, технические упругие константы которых приведены в таблице А.11.

ЯП имеет форму прямоугольного параллелепипеда, длины двух ребер которого равны 1, а длина третьего ребра равна 0,1 (безразмерные величины). Суммарная объемная концентрация нитей в ЯП равна 0,6.

Ортонормированная система координат Oe_i, i = 1, 2, 3, ориентирована по ребрам ЯП.

Задача заключается в том, чтобы вычислить:

- компоненты эффективного тензора модулей упругости композита ,

- эффективные технические упругие константы композита:

E_i - три модуля упругости в направлениях Oe_i,

G_ij - три модуля сдвига в плоскостях Oe_ie_j, i /= j,

- три коэффициента Пуассона, i /= j,

используя для этого верифицируемое ПО, основанное на решении данной задачи с помощью МКЭ.

Для решения задачи следует использовать постановку трехмерной квазистатической задачи теории линейной упругости (см. [9] - [12] с граничными условиями, которые зависят от применяемой модели вычисления компонент эффективного тензора модулей упругости и численной реализации этого метода.

Для решения задачи должен быть использован МКЭ.

Пример математической модели и постановки математической задачи для вычисления компонент эффективного тензора модулей упругости КМ на основе теории асимптотического осреднения приведен в [13] и [14].

Пример КЭМ размером 375 707 элементов, на основе тетраэдрального конечного элемента, построенной для 1/4 ЯП ТКМ с полотняным плетением нитей, показан на рисунке А.9. Полная ЯП строится отражением 1/4 ЯП относительно плоскостей Oe₁e₃ и Oe₂e₃.

А.1.5.2 Результаты

Результаты решения данной задачи с помощью верифицируемого ПО сравнивают со следующими вариантами расчетов:

а) результатами предельных вариантов расчетов, в которых характеристики матрицы и нитей выбирают одинаковыми, в этом случае:

1) эффективные модули упругости композита совпадают с модулями упругости C_ijkl компонентов композита [формула (А.7)];

2) эффективные технические упругие константы композита совпадают с техническими упругими константами выбранного компонента КМ [формулы (А.8) и (А.9)];

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

Пример сравнения расчетов по методу 2 приведен в таблице А.12.

А.2.1.1 Описание и постановка задачи

Рассматривается СКМ, ЯП которого состоит из четырех различных компонентов (слоев) (рисунок А.1).

Каждый слой является изотропным материалом, коэффициенты теплопроводности которых приведены в таблице А.13.

ЯП имеет форму куба с длиной ребра равной 1 (безразмерные координаты). Толщина каждого слоя равна 0,25.

Ортонормированная система координат Oe_i, i = 1, 2, 3 ориентирована следующим образом: векторы e₁, e₂ коллинеарны базовой плоскости, а вектор e₃ ортогонален к этой плоскости.

Задача заключается в том, чтобы вычислить компоненты эффективного тензора теплопроводности композита в направлениях Oe₁, используя для этого верифицируемое ПО, основанное на конечно-элементном методе решения данной задачи.

Для решения задачи следует использовать постановку трехмерной стационарной задачи линейной теории теплопроводности (см. [9] - [12]) с граничными условиями, которые зависят от применяемой модели вычисления компонент эффективного тензора модулей упругости и численной реализации этого метода.

Для решения задачи необходимо использовать МКЭ.

Пример математической модели и постановки математической задачи для вычисления компонент эффективного тензора теплопроводности КМ на основе теории асимптотического осреднения приведен в [15].

Пример КЭМ размером 206 254 элемента показан на рисунке А.2. КЭМ построена на основе тетраэдальных КЭ.

А.2.1.2 Результаты

Результаты решения данной задачи с помощью верифицируемого ПО сравнивают со следующими вариантами расчетов:

а) результатами аналитического решения задачи, полученными с помощью теории асимптотического осреднения для СКМ см. [[1], формула (3.11)]

где - компоненты тензора модулей упругости компонентов композита, представленные в виде разрывной функции координат;

- осреднение по ЯП;

б) результатами предельных вариантов расчетов, в которых характеристики всех слоев выбираются одинаковыми (например, могут быть выбраны характеристики слоя 1 или 2 из таблицы А.1), в этом случае компоненты эффективного тензора теплопроводности композита совпадают с компонентами тензора теплопроводности выбранного слоя (например, 1 или 2, 3, 4)

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

в) результатами численного решения данной задачи с помощью другого ПО, ранее прошедшего верификацию на данных задачах.

Пример сравнения расчетов по методу 2 приведен в таблице А.14.

А.2.2.1 Описание и постановка задачи

Рассматривается СВКМ, ЯП которого состоит из четырех различных компонентов (слоев) (см. рисунок А.1).

Каждый слой СВКМ представляет собой 1D КМ (см. рисунок А.3), с эффективными характеристиками. Каждый слой имеет угол армирования , значения которого указаны в таблице А.3.

1D КМ в каждом слое является одним и тем же трансверсально-изотропным материалом, отличается только углом армирования, коэффициенты теплопроводности 1D КМ приведены в таблице А.15.

Единая для всех слоев СВКМ ортонормированная система координат Oe_i, i = 1, 2, 3 ориентирована следующим образом: векторы e₁, e₂ коллинеарны базовой плоскости слоев, а вектор e₃ ортогонален к этой плоскости.

Собственная система координат каждого слоя ориентирована следующим образом: вектор для всех слоев совпадает с e₃, вектор направлен под углом к вектору e₁, где s - номер слоя в ЯП (см. рисунок А.3).

ЯП имеет форму куба с длиной ребра, равной 1 (безразмерные координаты). Толщина каждого слоя равна 0,25.

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

Для решения задачи следует использовать постановку трехмерной стационарной задачи линейной теории теплопроводности (см. [9] - [12]) с граничными условиями, которые зависят от применяемой модели вычисления компонент эффективного тензора модулей упругости и численной реализации этого метода.

Для решения задачи должен быть использован МКЭ.

Пример математической модели и постановки математической задачи для вычисления компонент эффективного тензора теплопроводности КМ на основе теории асимптотического осреднения приведен в [15].

Пример КЭМ размером 206 254 элементов показан на рисунке А.2.

А.2.2.2 Результаты

Результаты решения данной задачи с помощью верифицируемого ПО сравнивают со следующими вариантами расчетов:

а) результатами аналитического решения задачи для компоненты эффективного тензора теплопроводности СВКМ, полученными с помощью теории асимптотического осреднения для СКМ [см. [1], формула (3.11)].

где - компоненты тензора теплопроводности s-го слоя, вычисляемые по формуле

где - осреднение по ЯП;

- компоненты тензора теплопроводности 1D КМ в собственной системе координат , - матрица поворота системы координат слоя, определяемые формулами (А.12), (А.13);

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

Пример сравнения расчетов по методу 2 приведен в таблице А.16.

А.2.3.1 Описание и постановка задачи

Рассматривается ТКМ, ЯП которого состоит из пяти различных компонентов: четырех нитей и матрицы (см. рисунок А.8). Форма нитей близка к синусоидальной, две нити ориентированы в направлении Oe₁, а две другие - в направлении Oe₂. Система нитей образует взаимное переплетение без касания.

Каждый компонент является изотропным материалом, константы теплопроводности которых приведены в таблице А.17.

ЯП имеет форму куба с длиной ребра, равной 1 (безразмерные координаты). Суммарная объемная концентрация нитей в ЯП равна 0,6.

Ортонормированная система координат Oe_i, i = 1, 2, 3 ориентирована по ребрам ЯП.

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

Для решения задачи следует использовать постановку трехмерной стационарной задачи линейной теории теплопроводности (см. [9] - [12]) с граничными условиями, которые зависят от применяемой модели вычисления компоненты эффективного тензора модулей упругости и численной реализации этого метода.

Для решения задачи должен быть использован МКЭ.

Пример математической модели и постановки математической задачи для вычисления компонент эффективного тензора теплопроводности КМ на основе теории асимптотического осреднения приведен в [15].

Пример КЭМ размером 375 707 элементов показан на рисунке А.9.

А.2.3.2 Результаты

Результаты решения данной задачи с помощью верифицируемого ПО сравнивают со следующими вариантами расчетов:

а) результатами предельных вариантов расчетов, в которых характеристики всех слоев выбираются одинаковыми (например, могут быть выбраны характеристики слоя 1 или 2 из таблицы А.1), в этом случае компоненты эффективного тензора теплопроводности композита совпадают с компонентами тензора теплопроводности выбранного слоя (например, 1 или 2, 3, 4) [см. формулу (А.15)];

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

Пример сравнения расчетов по методу 2 приведен в таблице А.18.

А.2.4.1 Описание и постановка задачи

Рассматривается СКМ, ячейка периодичности которого состоит из четырех различных компонентов (слоев) (см. рисунок А.1).

Каждый слой является изотропным материалом, технические упругие константы и коэффициенты теплового расширения которых приведены в таблице А.19.

ЯП имеет форму куба с длиной ребра, равной 1 (безразмерные координаты). Толщина каждого слоя равна 0,25.

Ортонормированная система координат Oe_i, i = 1, 2, 3 ориентирована следующим образом: векторы e₁, e₂ коллинеарны базовой плоскости, а вектор e₃ - ортогонален к этой плоскости.

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

Для решения задачи следует использовать постановку трехмерной квазистатической задачи линейной теории термоупругости (см. [9] - [12]) с граничными условиями, которые зависят от применяемой модели вычисления компонент эффективного тензора теплового расширения и численной реализации этого метода.

Для решения задачи должен быть использован МКЭ.

Пример математической модели и постановки математической задачи для вычисления компонент эффективного тензора теплового расширения КМ на основе теории асимптотического осреднения приведен в [16].

Пример КЭМ размером 206 254 элемента на основе тетраэдального КЭ показан на рисунке А.2.

А.2.4.2 Результаты

Результаты решения данной задачи с помощью верифицируемого ПО сравнивают со следующими вариантами расчетов:

а) результатами аналитического решения задачи, полученными с помощью теории асимптотического осреднения для СКМ [см. [1], формула (3.11)]

где - компоненты эффективного тензора теплового расширения композита;

- компоненты тензора теплового расширения слоев композита, представленные в виде разрывной функции координат в рамках ЯП;

C_ijkl - компоненты тензора модулей упругости слоев композита, также представленные в виде разрывной функции координат в рамках ЯП;

- компоненты эффективного тензора упругих податливостей композита;

- обратная матрица;

- осреднение по ЯП;

б) результатами предельных вариантов расчетов, в которых характеристики всех слоев выбираются одинаковыми (например, могут быть выбраны характеристики слоя 1 или 2 из таблицы А.1), в этом случае компоненты эффективного тензора теплового расширения композита совпадают с компонентами тензора теплового расширения выбранного слоя (например, 1 или 2, 3, 4)

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

в) результатами численного решения данной задачи с помощью другого ПО, ранее прошедшего верификацию на данных задачах.

Пример сравнения расчетов по способу 2 приведен в таблице А.20.

А.2.5.1 Описание и постановка задачи

Рассматривается СВКМ, ЯП которого состоит из четырех различных компонентов (слоев) (см. рисунки А.1 и А.3).

Каждый слой СВКМ представляет собой 1D КМ с эффективными характеристиками. Каждый слой имеет угол армирования , значения которого указаны в таблице А.3.

1D КМ в каждом слое является одним и тем же линейно-термо-упругим трансверсально-изотропным материалом, отличается только углом армирования, технические упругие константы 1D КМ приведены в таблице А.4, а коэффициенты теплового расширения 1D КМ - в таблице А.21.

Единая для всех слоев СВКМ ортонормированная система координат Oe_i, i = 1, 2, 3 ориентирована следующим образом: векторы e₁, e₂ коллинеарны базовой плоскости, а вектор e₃ ортогонален к этой плоскости. Собственная система координат каждого слоя ориентирована следующим образом: вектор для всех слоев совпадает, векторы направлены под углами к вектору e₁, где s - номер слоя в ЯП.

ЯП имеет форму куба с длиной ребра, равной 1 (безразмерные координаты). Толщина каждого слоя равна 0,25.

Задача заключается в том, чтобы вычислить компоненты эффективного тензора теплового расширения композита в направлениях Oe_i, используя для этого верифицируемое ПО, основанное на МКЭ.

Для решения задачи следует использовать постановку трехмерной квазистатической задачи линейной теории термоупругости (см. [9] - [12]) с граничными условиями, которые зависят от применяемой модели вычисления компонент эффективного тензора теплового расширения и численной реализации этого метода.

Для решения задачи должен быть использован МКЭ.

Пример математической модели и постановки математической задачи для вычисления компонент эффективного тензора теплового расширения КМ на основе теории асимптотического осреднения приведен в [16].

Пример КЭМ размером 206 254 элемента показан на рисунке А.2.

А.2.5.2 Результаты

Результаты решения данной задачи с помощью верифицируемого ПО сравнивают со следующими вариантами расчетов:

а) результатами аналитического решения задачи, полученными с помощью теории асимптотического осреднения для СКМ [см. [1], формула (3.15)].

где - компоненты тензора модулей упругости;

- компоненты тензора теплового расширения s-го слоя в единой системе координат Oe_i вычисляемые по формулам

где - компоненты тензора модулей упругости;

- компоненты тензора теплового расширения 1D КМ в собственной системе координат 

- матрица поворота системы координат слоя, определяемая формулами (А.12), (А.13).

В формулах (А.23) - (А.26) также обозначены: - осреднение по ЯП, - обратная матрица, - компоненты эффективного тензора теплового расширения композита, - компоненты тензора теплового расширения слоев композита, представленные в виде разрывной функции координат в рамках ЯП, C_ijkl - компоненты тензора модулей упругости слоев композита, представленные в виде разрывной функции координат в рамках ЯП, - компоненты эффективного тензора упругих податливостей композита, - обратная матрица, - осреднение по ЯП;

б) результатами предельных вариантов расчетов, в которых характеристики всех слоев выбираются одинаковыми (например, могут быть выбраны характеристики слоя 1 или 2 из таблицы А.1), в этом случае компоненты эффективного тензора теплового расширения композита совпадают с компонентами тензора теплового расширения выбранного слоя (например, 1 или 2, 3, 4) [см. формулу (А.22)];

в) результатами численного решения данной задачи с помощью другого ПО, ранее прошедшего верификацию на данных задачах.

Пример сравнения расчетов по методу 3 приведен в таблице А.22.

А.2.6.1 Описание и постановка задачи

Рассматривается ТКМ, ЯП которого состоит из пяти различных компонентов: четырех нитей и матрицы (см. рисунок А.8). Форма нитей близка к синусоидальной, две нити ориентированы в направлении Oe₁, а две другие - в направлении Oe₂. Система нитей образует взаимное переплетение без касания.

Каждый компонент является упругим изотропным материалом, технические упругие константы и коэффициенты теплового расширения которых приведены в таблице А.23.

ЯП имеет форму куба с длиной ребра, равной 1 (безразмерные координаты). Толщина каждого слоя равна 0,25.

Ортонормированная система координат Oe_i, i = 1, 2, 3 ориентирована следующим образом: векторы e₁, e₂ коллинеарны базовой плоскости, а вектор e₃ ортогонален к этой плоскости.

Задача заключается в том, чтобы вычислить компонент эффективного тензора теплового расширения композита в направлениях Oe_i, используя для этого верифицируемое ПО, основанное на МКЭ.

Для решения задачи следует использовать постановку трехмерной квазистатической задачи линейной теории термоупругости (см. [9] - [12]) с граничными условиями, которые зависят от применяемой модели вычисления компонент эффективного тензора теплового расширения и численной реализации этого метода.

Для решения задачи должен быть использован МКЭ.

Пример математической модели и постановки математической задачи для вычисления компонент эффективного тензора теплового расширения КМ на основе теории асимптотического осреднения приведен в [16].

Пример конечно-элементной тетраэдральной сетки размером 206 254 элемента показан на рисунке А.9.

А.2.6.2 Результаты

Результаты решения данной задачи с помощью верифицируемого ПО сравнивают со следующими вариантами расчетов:

а) результатами предельных вариантов расчетов, в которых характеристики всех слоев выбираются одинаковыми (например, могут быть выбраны характеристики слоя 1 или 2 из таблицы А.23), в этом случае компоненты эффективного тензора теплового расширения композита совпадают с компонентами тензора теплового расширения выбранного слоя (например, 1 или 2, 3, 4) [см. формулу (А.22)];

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

Пример сравнения расчетов по методу 2 приведен в таблице А.24.