2255
.pdf
BUILDING MATERIALS AND PRODUCTS
Традиционно под оптимальной системой понимается наилучшая в принятом смысле система. Оценка качества сложной системы в общем случае является многокритериальной и противоречивой. В простейшем случае (однокритериальная система) среди всех возможных вариантов системы требуется найти оптимальный (наилучший при принятом критерии, однозначно характеризующий возможный вариант реализации системы). При однокритериальной оптимизации каждому варианту исполнения системы может быть поставлено в соответствие некоторое число. Цель управления – достижение экстремума критерия оптимальности. Для многокрите-
риальной системы в общем случае выбранное решение всегда будет компромиссным (противоречивость частных критериев: значения одних критериев улучшаются, а других – ухудшаются).
Обычно в строительном материаловедении аналитические зависимости частных критериев (отдельные свойства) аппроксимируются полиномиальными моделями
|
|
k |
|
k |
|
|
b0 bi xi bi j xi xj bii xi2 , |
||
Y |
||||
|
|
i 1 |
i j |
i 1Б |
где b0 ,bi ,bij ,bii определяются с использованием статистических методов планирова-
ния эксперимента; – случайная величина, учитывающая совокупность ошибок эксперимента. Решаются задачи однокритериальной оптимизации свойств. К сожа-
лению, очень редко производится многокритериальная оптимизация; анализируется зависимость лишь отдельных свойств, а не обобщенного критерия качества материала.
Свойства материала определяются рецептурно-технологическими параметрами (факторами) [1]; их прогнозирование требует решения двух задач:
–оптимизация параметров принятой модели для выбранного вида критерия качества (параметрическая идентификация целевой функции);
–составление интерполяционной модели коэффициентов целевой функции от рецептурно-технологических параметров.
При решении этих задач обычно используются методы математического планирования эксперимента: строятся модели для каждого из частных критериев, определяются интерполяционные полиномы для линий равного уровня (эллипс, гипербола, парабола в двухфакторном пространстве). В частности, зависимость частного
критерия Z от рецептурно-технологических параметров X1, , Xm можно описать при их малых изменениях линейной моделью вида
Z a0 a1 X1 a2 X2 am Xm . |
( 1 ) |
||||||||||
Используя серию измерений величин Z , |
X j , |
можно определить параметры ai . По r |
|||||||||
|
|
|
|
и |
|
|
|
– средние значения Z , |
X j для ука- |
||
совокупностям измерений вычисляются Z |
X |
j |
|||||||||
|
|
|
|
|
|
. Тогдасоотношение(1) приметвид |
|||||
занныхсерийизмерений. Введём z Z Z |
|
, x X X |
|||||||||
z a1x1 a2 x2 am xm , |
|
||||||||||
или в векторной форме |
|
|
|
|
|
|
|
|
|||
z xт a , |
|
|
|
|
|||||||
где x, a – вектор-столбцы с элементами xj, aj соответственно. |
|
||||||||||
Количествоr последовательныхизмеренийудовлетворяетсоотношениям |
|
||||||||||
z 1 |
x 1 Ta, |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
z μ |
x μ Ta, |
(2) |
|||||||||
z r x r Ta,
xт 1, r .
Regional architecture and engineering 2018 |
№2 51 |
СТРОИТЕЛЬНЫЕ МАТЕРИАЛЫ И ИЗДЕЛИЯ
Введём вектор и матрицу U следующим образом:
z 1 z z r т ;
|
x 1 T |
|
|
x11 |
xj1 |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
T |
|
|
|
|
|||
U |
|
|||||||
x |
|
|
x1 |
xj |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
xjr |
|
|
x r T |
|
x1r |
|||||
|
|
|
|
|
|
|
|
|
|
xm1 |
|
|
|
|
xm .
xmr
Тогда система уравнений (2) может быть записана в векторной форме:
Ua . |
(3) |
Предполагая, что компоненты вектора a в уравнении (3) |
являются оценками aˆ |
истинного вектора a , можно с помощью уравнения (3) получить такие оценки χˆ
вектора , что
χˆ Uaˆ .
Легко показать, что наилучшая в смысле наименьших квадратов оценка aˆ вектора a удовлетворяет уравнению
UтUaˆ Uт ,
так что
aˆ UT U 1 UT χ aˆ1, aˆ2 , , aˆm T ,
что и позволяет построить процедуру идентификации вектора a на основе линейной рагрессии и метода наименьших квадратов. Отметим, что матрица UT U 1
существует только тогда, когда матрица U не является особенной.
Число измерений r должно быть больше числа идентифицируемых параметровr m . Если r m , то в оценке шум измерений не будет сглажен. Поэтому для
адекватной идентификации требуется, по крайней мере, m 1 измерение, причём в течение этого периода система предполагается стационарной.
|
Так, указанным подходом [2, 3] определялась зависимость предела прочности при |
||||||||||||||||||||||||||
сжатии Rcж |
от твёрдости T и модуля деформации Е эпоксидных композитов по |
||||||||||||||||||||||||||
данным эксперимента (табл. 1). |
|
|
|
|
|
|
|
|
|
Т а б л и ц а |
1 |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Z |
117 |
|
|
100 |
|
|
120 |
57 |
|
|
99 |
102 |
79 |
64 |
74 |
87 |
|
||||||||||
X1 |
6,01 |
|
|
5,05 |
|
|
6,15 |
2,94 |
|
|
5,06 |
5,23 |
4,05 |
3,40 |
3,79 |
4,44 |
|
||||||||||
X2 |
3,62 104 |
4,71 104 |
3,51 104 |
1,06 104 |
4,71 104 |
4,48 104 |
6,57 104 |
8,52 104 |
7,25 104 |
5,71 104 |
|||||||||||||||||
|
Принято: |
|
X1 T, X2 E, Z Rсж . |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
1 = 4,61; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
Имеем: X |
X |
2 = 5,01 104; Z |
= 89,9. |
|
|
|
|
|
|
|
||||||||||||||||
|
Для центрированных значений переменных данные эксперимента приводятся в |
||||||||||||||||||||||||||
табл. 2. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Т а б л и ц а |
2 |
|
Z |
|
27,1 |
|
|
10,1 |
|
|
30,1 |
|
–32,9 |
|
9,1 |
|
12,1 |
|
–10,9 |
|
–25,9 |
–15,9 |
–2,9 |
|
||||||
X1 |
|
1,40 |
|
|
0,44 |
|
|
1,54 |
|
–1,67 |
|
0,45 |
|
0,62 |
|
–0,56 |
|
–1,21 |
–0,82 |
–0,17 |
|||||||
X2 |
|
–1,39 104 |
|
–0,3 104 |
|
–1,5 104 |
|
–3,95 104 |
|
–0,3 104 |
|
–0,53 104 |
|
1,56 104 |
|
3,51 104 |
2,24 104 |
0,7 104 |
|||||||||
52 Региональная архитектура и строительство 2018 |
№2 |
BUILDING MATERIALS AND PRODUCTS
В соответствии с предыдущим значение z при -измерении
z a1x1 a2 x2 ,
и для параметров линейной модели будем иметь:
aˆ aˆ1, aˆ2 = (19,1; – 1,26 10 –5).
Окончательно получим
Rcж = 2,48 + 19,1 T – 1,26 10 -5 E.
Как видим, относительная ошибка вычисления по модели Rcж не превышает
приблизительно 1 %. Одновременно отметим, что исходная таблица измерений обладает избыточностью (достаточно знания одного из двух параметров Rcж или T).
При этом Rcж и T практически от E не зависят.
Заметим, предлагаемая методика в основном ориентирована на изучение свойств материалов в локальных областях факторного пространства. Метод легко обобщается и на случай систем со многими входами и выходами, а также на случай идентифи-
кации динамических процессов, описываемых нормальной системой обыкновенных дифференциальных уравнений.
Все вычислительные процедуры сводятся к простому вычислению неизвестных коэффициентов с доступностью математических выводов. Недостатки определяются чувствительностью оценок метода наименьших квадратов к резким выбросам, встречающимся в экспериментальных данных. Естественно, всегда предполагается, что используемая в ходе регрессионного анализа форма связи наилучшим образом отражает зависимость выходных переменных Z от факторной переменной X.
В заключение приведем пример многокритериального синтеза эпоксидных композитов повышенной плотности.
Использование математических методов планирования эксперимента позволило получить следующие аналитические зависимости средней плотности , кг/м3, и
предела прочности Rсж, Мпа, на сжатие образцов 20 20 20 см:
X1, X2 3642,8 129,1X1 668,5X2 53, 2X1 X2 513,7X22 ,
Rсж X1, X2 118,5 19,5X1 20,9X2 3, 2X22 ,
где X1, X2 – кодированные значения соответственно концентрации x1 пластификатора (в % от массы смолы) и степени наполнения x2 (П:H по массе); x10 = 25, I1 = 25; x20 =1/7,
I2Н = 5. Линиями уровня X1, X2 const и Rсж X1, X2 const являются соот-
ветственно ветви гиперболы и параболы (квадратичные модели целевых функций). В качестве оптимальных значений факторов будут: x1 = 2,5; x2 = 10,2 (им соответствуют значения плотности = 3955 кг/м3 и предела прочности Rсж = 145 МПа).
Аналогично решается задача и в случае большего числа переменных и частных критериев качества [4…9].
Список литературы
1.Гарькина, И.А. Эволюция представлений о композиционных материалах с позиций смены парадигм / И.А. Гарькина, А.М. Данилов, Е.В. Королев // Строи-
тельные материалы. – 2018. – №1–2. – С.60–63.
2.Системный анализ в строительном материаловедении: монография / Ю.М. Баженов, И.А. Гарькина, А.М. Данилов, Е.В. Королев. – М.: МГСУ: Библиотека научных разработок и проектов, 2012. – 432 с.
3.Данилов, А.М. Математическое и компьютерное моделирование сложных систем / А.М. Данилов, И.А. Гарькина, Э.Р. Домке. – Пенза: ПГУАС, 2011. – 296 с.
4.Данилов, А.М. Методология проектирования сложных систем при разработке материалов специального назначения / А.М. Данилов, И.А. Гарькина // Известия высших учебных заведений. Строительство. – 2011. – № 1. – С. 80–85.
Regional architecture and engineering 2018 |
№2 53 |
СТРОИТЕЛЬНЫЕ МАТЕРИАЛЫ И ИЗДЕЛИЯ
5.Гарькина, И.А. Опыт разработки композиционных материалов: некоторые аспекты математического моделирования / И.А. Гарькина, А.М. Данилов // Известия высших учебных заведений. Строительство. – 2013. – № 8 (656). – С. 28–33.
6.Гарькина, И.А. Формализация оценки структуры и свойств композиционных материалов специального назначения / И.А. Гарькина // Строительные материалы. – 2007. – № 1. – С. 70–72.
7.Garkina, I.A. Modeling of kinetic processes in composite materials / I.A. Garkina // Contemporary Engineering Sciences. – 2015. – Vol. 8. – № 9. – P. 421–425.
8.Garkina, I.A. Modeling of building materials as complex systems / I.A. Garkina, A.M. Danilov, Y.P. Skachkov // Key Engineering Materials. – 2017. – Vol. 730. – P. 412.
9.Данилов, А.М. Математическое моделирование сложных систем: состояние, перспективы, пример реализации / А.М. Данилов, И.А. Гарькина // Вестник гражданских инженеров. – 2012. – № 2. – С. 333–337.
References
1.Garkina, I.A. Evolution of the concept of composite materials from the standpoint of changing paradigms / I.A. Garkina, A.M. Danilov, E.V. Korolev // Construction materials. – 2018. – №1–2. – P. 60–62.
2.The system analysis in construction materials science: monograph / Yu.M. Bazhenov, I.A. Garkina, A.M. Danilov, E.V. Korolev. – M.: MGSU: Library of scientific developments and projects, 2012. – 432 p.
3.Danilov, A.M. Mathematical and computer modeling of difficult systems / A.M. Danilov, I.A. Garkina, E.R. Domke. – Penza: PGUAS, 2011. – 296 p.
4.Danilov, A.M. Methodology of design of difficult systems when developing materials of a special purpose / A.M. Danilov, I.A. Garkina // News of higher educational institutions. Construction. – 2011. – №1. – P. 80–85.
5.Garkina, I.A. Experience of development of composite materials: some aspects of mathematical modeling / I.A. Garkina, A.M. Danilov // News of higher educational institutions. Construction. – 2013. – №. 8 (656). – P. 28–33.
6.Garkina, I.A. Formalization of assessment of structure and properties of composite materials of a special purpose / I.A. Garkina // Construction materials. – 2007. – №. 1. – P. 70–72.
7.Garkina, I.A. Modeling of kinetic processes in composite materials / I.A. Garkina // Contemporary Engineering Sciences. – 2015. – Vol. 8. – № 9. – P. 421–425.
8.Garkina, I.A. Modeling of building materials as complex systems / I.A. Garkina, A.M. Danilov, Y.P. Skachkov // Key Engineering Materials. – 2017. – Vol. 730. – P. 412.
9.Danilov, A.M. Mathematical modeling of difficult systems: state, prospects, example realization / A.M. Danilov, I.A. Garkina // Herald of Civil Engineers – 2012. – №. 2. – P. 333–337.
54 Региональная архитектура и строительство 2018 |
№2 |
BUILDING MATERIALS AND PRODUCTS
УДК 519.7: 691
Пензенский государственный университет |
Penza State University of Architecture |
архитектуры и строительства |
and Construction |
Россия, 440028, г. Пенза, ул. Германа Титова, д.28,
òåë.: (8412) 48-27-37; ôàêñ: (8421) 48-74-77
Гарькина Ирина Александровна, доктор технических наук, профессор
кафедры «Математика и математическое моделирование»
E-mail: fmatem@pguas.ru
Данилов Александр Максимович,
доктор технических наук, профессор, советник РААСН, зав. кафедрой «Математика и математическое моделирование»
E-mail: fmatem@pguas.ru
Russia, 440028, Penza, 28, German Titov St., tel.: (8412) 48-27-37; fax: (8412) 48-74-77
Garkina Irina Aleksandrovna,
Doctor of Sciences, Professor of the department «Mathematics and Mathematical Modeling» E-mail: fmatem@pguas.ru
Danilov Alexander Maksimovich, Doctor of Sciences, Professor, Adviser of the Russian Academy of Architectural and
Construction Sciences, Head of the department «Mathematics and Mathematical Modeling» E-mail: fmatem@pguas.ru
ФОРМАЛИЗАЦИЯ ПРОЦЕССОВ ФОРМИРОВАНИЯ СВОЙСТВ КОМПОЗИТОВ
И.А. Гарькина, А.М. Данилов
Предлагаются эффективные методы аппроксимации кинетических процессов формирования физико-механических характеристик композиционных материалов для решения задач как однокритериальной, так и многокритериальной оптимизации их структуры и свойств. Подробно рассматривается аппроксимация кинетических процессов на заданном интервале с определением частичных интервалов аппроксимации максимальной длины. Приводится алгоритм аппроксимации кубическими сплайнами. Для определения передаточных функций композитов как динамических систем со многими входами и выходами рекомендуется использовать полиномы Лагерра. Не исключается возможность и когнитивного подхода при непосредственном подборе вида аппроксимирующей функции.
Ключевые слова: композиционные материалы, свойства, кинетические процессы, параметры, аппроксимация, методы
FORMALIZATION THE PROCESSES OF COMPOSITES
PROPERTIES FORMATION
I.A. Garkina, A.M. Danilov
Effective methods are proposed for approximating the kinetic processes formation of physical and mechanical characteristics of composite materials for solving the problems of both single-criterion and multi-criteria optimization of their structure and properties. The approximation of kinetic processes on a given interval with definition of partial intervals of the maximum length approximation is considered in detail. An approximation algorithm on cubic splines is given. To determine the transfer functions of composites, as dynamic systems with many inputs and outputs, it is recommended to use Laguerre polynomials. The possibility of a cognitive approach is not excluded with direct selection of the form of the approximating function.
Keywords: composite materials, properties, kinetic processes, parameters, approximation, methods
Можно считать общепринятым подход к оптимизации структуры и свойств композиционных материалов по аналитическому описанию кинетических процессов формирования их физико-механических характеристик. Среди целого ряда методов аппроксимации этих процессов особое место в силу его простоты занимает кусочнолинейная аппроксимация, как правило, с равномерным шагом. Однако во многих
Regional architecture and engineering 2018 |
№2 55 |
СТРОИТЕЛЬНЫЕ МАТЕРИАЛЫ И ИЗДЕЛИЯ
случаях целесообразна разбивка кинетических процессов на отдельные участки,
исходя из обеспечения максимальной длины участков аппроксимации. В качестве параметров кинетических процессов могут использоваться различные характеристики, интегрально определяющие процессы формирования свойств [1…4]. Часто выход параметров на эксплуатационное значение носит экспоненциальный апериодический характер (рис. 1).
x |
|
|
Mm |
|
|
|
Mn |
|
xm |
|
|
x0 |
|
|
|
|
|
tm |
tn |
t |
Рис. 1. Апериодический кинетический процесс |
|
|
Ниже предлагается алгоритм кусочно-линейной аппроксимации кинетического процесса f x с интервалами максимальной длины на отрезке a,b .
Ввод данных. Задается интервал a,b , относительная погрешность аппрокси-
мации 0 в процентах. Функция f x |
реализуется в виде программы вычисления ее |
|||||||||||||||||||||
значений в любой точке x a,b |
или хотя бы в точках xj , j |
|
; x0 a, xN |
b , |
||||||||||||||||||
0, N |
||||||||||||||||||||||
таких, что |
fср ; fср 1 f j , |
f j f xj . |
|
|
||||||||||||||||||
|
|
f j f j 1 2 0 |
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N |
|
|
|
|
|
|
|
|
|
|
|
1000 |
|
|
|
|
|
|
N |
|
j 1 |
|
|
|
|
xj , f j N |
|
|||
|
|
|
|
|
|
|
|
|
|
|||||||||||||
Реализация алгоритма осуществлялась по табличным значениям |
при |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
j 0 |
|
|
|
|
|
|
|
|
0 |
|
1 |
N |
|
|
|
|
|
|
|||||||
абсолютной погрешности |
|
|
|
f j |
|
|
при |
плотности |
табулирования |
|||||||||||||
|
|
|
||||||||||||||||||||
100 N |
||||||||||||||||||||||
|
|
|
|
|
|
j 1 |
|
|
|
|
|
|
|
|
|
|
||||||
f j f j 1 0,2 .
Табулирование. Здесь реализована программа получения информации. Использовался нижеприведённый алгоритм.
Для h 0,04 b a , N b a вычисляются: |
|
|
|
|
|
|
|
|
|
|
|||||||
xj a jh,; |
|
h |
fср |
1 f j , |
|
0 fср . |
|||||||||||
f j f xj , j 0, N ; h max f j f j 1 ; |
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
N |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 j N |
|
|
|
|
|
N |
j 1 |
|
|
|
|
100 |
|
|
|
|
|
|
|
|
|
|
||||||||||
При h 0,2 в два раза уменьшалось значение h . |
|
|
|
|
|
|
|
|
|
|
|||||||
56 Региональная архитектура и строительство 2018 |
№2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
BUILDING MATERIALS AND PRODUCTS |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
Кусочно-линейная аппроксимация. Аппроксимирующая |
функция |
f x |
пред- |
||||||||||||||||||||||||||||||||||||||||||||||||
ставлялась в виде |
f x Ar Kr x zr , zr x zr 1, r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0, Na 1, |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
причем z0 |
x0 , zNa |
xN , |
f zr f zr ; |
|
f x f x |
|
|
|
для x a,b . |
|
|
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||
|
|
Функция |
f x определяется таблицей |
zr , Ar , Kr rNa0 ; |
|
zr – |
узлы, |
Na |
– |
число |
|||||||||||||||||||||||||||||||||||||||||||
узлов, |
|
Hr |
zr 1 zr |
– |
интервалы аппроксимации, |
|
|
|
Kr |
– |
|
угловые коэффициенты, |
|||||||||||||||||||||||||||||||||||||||||
z0 x0 , A0 |
f0 , Ar 1 Ar |
Kr Hr (непрерывность |
f x ). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
|
|
Таким образом, |
кусочно-линейная аппроксимация определяется параметрами Kr , Hr , |
||||||||||||||||||||||||||||||||||||||||||||||||||
r 0 . |
Максимальность |
Hr |
определяется выбором точки |
|
zr 1, r 0 , |
максимально |
|||||||||||||||||||||||||||||||||||||||||||||||
удаленной от zr ; |
|
|
f x f x |
|
, x zr , zr 1 ; |
f x f x , f zr 1 f zr 1 . |
|||||||||||||||||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
Вычисление zr , Ar , Kr . Принимается z0 |
x0 , A0 f0 . По вычисленным значениям |
||||||||||||||||||||||||||||||||||||||||||||||||||
zr , Ar |
|
находятся |
zr 1 |
, |
Hr |
zr 1 |
zr , |
|
Ar 1 f zr 1 , |
|
Kr |
|
Ar 1 Ar |
|
(точка zr 1 |
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Hr |
|
|
|
|
|
|
является одной из точек табулирования xj ). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
Алгоритм вычисления zr 1 |
|
сводится к выбору |
|
точки табулирования |
xj |
zr с |
|||||||||||||||||||||||||||||||||||||||||||||
учетом условия |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f j Ar |
xi |
zr |
|
|
; |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
max |
fi Ar |
|
|
|
|
|
|
|
|
|
|
|
|
(1) |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
jr i j |
|
|
|
|
|
xj zr |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
jr |
– номер точки табулирования xj r , соответствующей zr |
zr |
xjr |
. |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||
|
|
Как только при некотором |
|
j условие (1) нарушится, то запоминаем |
xj |
как zr0 1 |
|||||||||||||||||||||||||||||||||||||||||||||||
(номер запоминается как |
j0 |
|
) |
|
и принимается K |
0 |
|
f j |
Ar |
, A0 |
f |
j |
. Далее прове- |
||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r 1 |
|
|
|
|
|
|
|
|
|
|
|
r |
|
|
|
xj |
zr |
|
|
r 1 |
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
ряется условие (1) для xj |
zr0 1 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Если |
|
|
|
xj zr 1 |
, xN |
|
|
|
условие |
|
|
(1) |
не |
|
|
|
|
|
выполнено, |
|
|
принимается |
|||||||||||||||||||||||||||||
z |
r |
1 |
z0 |
, |
A |
|
A0 |
, K |
r |
K |
0 |
, |
|
j |
j0 |
и переходим к вычислению z |
|
|
. |
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
r 1 |
|
|
r 1 |
|
|
|
r 1 |
|
|
r |
|
|
r 1 |
|
r 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r |
2 |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
Если для некоторого |
xj zr 1 |
, xN |
выполняется условие (1), то осуществляется |
||||||||||||||||||||||||||||||||||||||||||||||||
переход к предыдущему этапу ( r не увеличивается). |
|
|
|
|
|
|
|
|
N 1 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
Так |
будут определены |
|
все |
триады |
|
|
|
|
|
|
|
|
|
|
|
|
. |
В |
последней |
триаде |
|||||||||||||||||||||||||||||||
|
|
|
|
zr , Ar , Kr r a0 |
|||||||||||||||||||||||||||||||||||||||||||||||||
zNa , ANa , KNa достаточно вычисления zNa |
zr 1 |
xN , Kr |
|
KNa 1 . |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||
|
|
Заметим, |
если |
f x |
удовлетворяет |
|
условию |
|
f |
x |
|
x a,b , |
то |
вместо |
|||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||
условия |
|
(1) |
|
|
|
(используется |
усредненная |
|
|
|
|
относительная |
|
|
погрешность |
||||||||||||||||||||||||||||||||||||||
|
|
|
f |
f |
|
100% ) можно воспользоваться и условием |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
o |
|
|
|
|
fср |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
max |
fi |
|
Ar |
f j Ar |
xi zr |
|
|
|
0 |
|
fi |
|
. |
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
jr i j |
|
|
|
|
|
xj |
zr |
|
|
|
|
|
|
|
|
|
100 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Regional architecture and engineering 2018 |
№2 57 |
СТРОИТЕЛЬНЫЕ МАТЕРИАЛЫ И ИЗДЕЛИЯ
Если при выполнения условия (1) |
xj xN , то |
|
|
|
|||||||
z |
z |
r 1 |
x |
N |
, K |
r |
K |
Na 1 |
|
fN Ar |
; |
|
|||||||||||
|
Na |
|
|
|
|
xN zr |
|||||
|
|
|
|
|
|
|
|
|
|
||
вычисления завершаются.
Приведем и другой возможный вариант аппроксимации кинетического процесса,
основанный на использовании интерполяции многочленами на отдельных участках, а
не на основе построения глобального интерполяционного многочлена на всём промежутке. Получаемые при этом гладкие кусочно-многочленные функции (составленные из многочленов одной и той же степени) являются сплайнами. Ограничимся куби-
ческими сплайнами (составлен из многочленов третьей степени). Ниже предлагается алгоритм построения кубических сплайнов класса C2 (дважды непрерывно дифференцируемые функции).
Здесь |
предполагается, |
что |
на |
отрезке |
a,b |
в |
узлах |
сетки |
||||||
a x1 |
x2 |
, , xN b заданы значения функции |
f x : fi |
f xi ,i |
|
. По |
||||||||
1, N |
||||||||||||||
предыдущему |
интерполяционный кубический сплайн |
S x должен |
удовлетворять |
|||||||||||
условиям: |
S xi |
fi , S xi f xi , S xi f xi . |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
(2) |
|||||||
Кубический сплайн S x |
на каждом из отрезков xi , xi 1 определяется четырьмя |
|||||||||||||
коэффициентами; для его построения на всём промежутке a,b |
требуется определить |
|||||||||||||
4 N коэффициента. Из условия принадлежности сплайна классу C2 предполагается |
||||||||||||||
непрерывность во всех внутренних узлах интерполяции |
xi , i |
|
|
только |
||||||||||
1, N 1 не |
||||||||||||||
сплайна |
S x , но и его производных S x и S x . |
Из этих условий получим |
||||||||||||
3 N 1 |
уравнений для определения неизвестных коэффициентов сплайна. Добавив |
|||||||||||||
N 1 уравнение из (2), будем иметь |
4N 2 уравнения; два недостающих уравнения |
|||||||||||||
получим из краевых условий (ограничения на значения сплайна и его производных на концах промежутка a,b ). Наиболее употребительны следующие краевые условия:
I. S a f a , S b f b . II. S a f a , S b f b . III. S k a S k b , k 1, 2 .
IV. S xp 0 S xp 0 , p 1, N 1.
(так называемые “естественные” условия имеют вид: S a 0, S b 0 ). На каждом из промежутков xi , xi 1 сплайн представим не в общем виде
S x ai bi x ci x2 di x3 ,
а в некотором специальном, позволяющем уменьшить число неизвестных коэффициентов сплайна. Для этого введём обозначения:
S xi mi ,i |
|
; hi xi 1 |
xi , t |
x xi |
. |
|
0, N |
||||||
|
||||||
|
|
|
|
hi |
||
На отрезке xi , xi 1 кубический сплайн можно записать в виде
S x fi 1 t 2 1 2t fi 1t2 3 2t mi hit 1 t 2 mi 1t2 1 t hi .
На каждом из промежутков xi , xi 1 он непрерывен вместе со своей первой производной всюду на a,b . Выберем параметры таким образом, чтобы и вторая
58 Региональная архитектура и строительство 2018 |
№2 |
BUILDING MATERIALS AND PRODUCTS
производная была непрерывна во всех внутренних узлах. Получим следующую систему уравнений:
|
i mi 1 2mi i mi 1 |
|
|
i |
fi 1 fi |
|
i |
fi fi 1 |
|
|
|
||||||||||
|
3 |
|
|
, |
(3) |
||||||||||||||||
|
|
|
hi 1 |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
hi |
|
|
|
|
|
|
|
|||
где |
|
hi 1 |
|
|
|
|
|
|
|
hi |
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
, i 1 i |
|
|
|
, |
|
i |
|
. |
|
|
|
||||||
|
|
|
|
0, N |
|
|
|
||||||||||||||
|
|
|
hi 1 hi |
|
|
|
|||||||||||||||
|
|
hi 1 hi |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
К уравнениям (3) следует добавить уравнения |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
2m m |
3 |
f1 f0 |
, |
m |
N 1 |
2m |
N |
3 |
fN fN 1 |
, |
|
|
|
||||||||
|
|
|
|
|
|||||||||||||||||
0 |
|
1 |
|
|
h0 |
|
|
|
|
|
hN 1 |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
полученные из граничных условий.
Задача построения кубического сплайна свелась к решению линейной системы для неизвестных коэффициентов mi:
2m0 m1 3 f1 f0 ,
h0
i mi 1 2mi i mi 1 |
3 |
|
i |
fi 1 fi |
i |
fi fi 1 |
|
, |
(4) |
|
|
|
|||||||||
|
hi 1 |
|||||||||
|
|
|
|
hi |
|
|
|
|||
mN 1 2mN 3 fN hN f1N 1 .
Нами при решении практических задач система (4) решалась методом прогонки (трехдиагональность матрицы системы).
При решении целого ряда задач по изучению композиционных материалов как нестационарных динамических систем возникает необходимость аппроксимации не только входных и выходных сигналов («черный ящик»), но и определение вида частных передаточных функций (параметры связи выходных сигналов с входными). Нами для этого использовалось разложение кинетического процесса по различным ортонормированным системам специальных функций, предпочтение отдавалось полиномам Лаггера [5]:
|
|
|
|
n |
|
|
n! t |
|
|
|
|
|
t e t dt 0 при |
m n |
|
|||||
L |
t |
|
|
|
|
|
|
|
; L |
t |
|
L |
, |
|||||||
|
|
|
2 |
|
|
|
|
|
||||||||||||
n |
|
|
|
! |
n ! |
n |
|
m |
|
|
m n |
|
||||||||
|
|
|
|
0 |
|
0 |
|
|
|
|
1 при |
|
||||||||
так как в этом случае передаточная функция W p |
легко определяется как дробно- |
|||||||||||||||||||
рациональная функция |
|
Y p |
с простыми правилами вычисления коэффициентов. |
|||||||||||||||||
|
X p |
|
||||||||||||||||||
При решении практических задач важную роль играет и человеческий фактор: выбор вида аппроксимирующей функции часто определяется интуицией экспериментатора. Так, по данным таблицы нами определялась аналитическая зависимость вязкости B эпоксидного композита от температуры t, oC, и процентного содержания x, %, специальной добавки.
t |
50 |
55 |
60 |
65 |
70 |
75 |
80 |
85 |
90 |
95 |
100 |
x |
|
|
|
|
|
|
|
|
|
|
|
0 |
172 |
119 |
84 |
65 |
50 |
40 |
30 |
25 |
22 |
20 |
19 |
1 |
210 |
145 |
105 |
80 |
65 |
52 |
43 |
35 |
30 |
27 |
25 |
5 |
250 |
160 |
115 |
90 |
75 |
60 |
48 |
40 |
33 |
30 |
28 |
10 |
97 |
65 |
45 |
33 |
22 |
17 |
13 |
10 |
11 |
9 |
8 |
15 |
132 |
91 |
62 |
47 |
36 |
30 |
21 |
18 |
17 |
15 |
12 |
Regional architecture and engineering 2018 |
№2 59 |
СТРОИТЕЛЬНЫЕ МАТЕРИАЛЫ И ИЗДЕЛИЯ
Исходя из когнитивного анализа, аппроксимирующая функция задавалась в виде B B x,t . Функция B B x,t const представлялась в виде полинома третьей
степени, а B B x const,t – функциями вида y ae kt , k k x при различных
степенях полинома. По этим результатам было получено представление
B c0 c1x c2 x2 c3 x3 er x t 50 , 50 t 100.
Предложенные методы прошли апробацию при разработке ряда композиционных материалов специального назначения [6, 7].
Список литературы
1.Garkina, I.A. Modeling of kinetic processes in composite materials / I.A. Garkina // Contemporary Engineering Sciences. – 2015. – Vol. 8. – № 9. – P. 421–425.
2.Garkina, I.A. Mathematical Methods of System Analysis in Construction Materials / I.A. Garkina, A.M. Danilov // IOP Conf. Series: Materials Science and Engineering. – 2017. – №245. – Р.062014.
3.Garkina, I.A. Analytical design of building materials / I.A. Garkina, A.M. Danilov // Journal of Basic and Applied Research International. – 2016. – Vol.18. – Issue 2. – P. 95–99.
4.Гарькина, И.А. Опыт разработки композиционных материалов: некоторые аспекты математического моделирования / И.А. Гарькина, А.М. Данилов // Известия высших учебных заведений. Строительство. – 2013. – № 8 (656). – С. 28–33.
5.Бейтмен, Г. Высшие трансцендентные функции / Г.Бейтмен, А.Эрдейи. – М.:
Наука, 1973. – Т. 2. – 297 с.
6.Данилов, А.М. Методология проектирования сложных систем при разработке материалов специального назначения / А.М. Данилов, И.А. Гарькина // Известия высших учебных заведений. Строительство. – 2011. – № 1. – С. 80–85.
7.Гарькина, И.А. Опыт разработки композиционных материалов: некоторые аспекты математического моделирования / И.А. Гарькина, А.М. Данилов // Известия высших учебных заведений. Строительство. – 2013. – № 8 (656). – С. 28–33.
References
1.Garkina, I.A. Modeling of kinetic processes in composite materials / I.A. Garkina // Contemporary Engineering Sciences. – 2015. – Vol. 8. – № 9. – P. 421–425.
2.Garkina, I.A. Mathematical Methods of System Analysis in Construction Materials / I.A. Garkina, A.M. Danilov // IOP Conf. Series: Materials Science and Engineering. – 2017. – №245. – Р.062014.
3.Garkina, I.A. Analytical design of building materials / I.A. Garkina, A.M. Danilov // Journal of Basic and Applied Research International. – 2016. – Vol.18. – Issue 2. – P. 95–99.
4.Garkina, I.A. Experience of development of composite materials: some aspects of mathematical modeling / I.A. Garkina, A.M. Danilov // News of higher educational institutions. Construction. – 2013. – No. 8 (656). – P.28–33.
5.Beytmen, G. The highest transcendental functions / G. Beytmen, A. Erdeyi. – M.: Science, 1973. – Vol. 2. – 297 p.
6.Danilov, A.M. Methodology of design of difficult systems when developing materials of a special purpose / A.M. Danilov, I.A. Garkina // News of higher educational institutions. Construction. – 2011. – №. 1. – P. 80–85.
7.Garkina, I.A. Experience of development of composite materials: some aspects of mathematical modeling / I.A., Garkina A.M. Danilov // News of higher educational institutions. Construction. – 2013. – №. 8 (656). – P. 28–33.
60 Региональная архитектура и строительство 2018 |
№2 |
