Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

2257

.pdf
Скачиваний:
0
Добавлен:
16.06.2024
Размер:
6.15 Mб
Скачать

BUILDING MATERIALS AND PRODUCTS

В задаче квадратичного программирования (частный случай задачи нелинейного

программирования (НЛП)) минимизируется сумма линейной и квадратичной форм q x cj xj d j k xj xk

j

j k

при ограничениях вида линейных неравенств

ai k xk bi 0, i

 

 

1, m

k

 

 

 

 

 

и неотрицательности переменных

 

 

 

 

 

 

xi 0, j

 

.

 

1, n

Эффективные вычислительные методы для задачи НЛП, к сожалению, существуют, лишь когда целевая функция имеет единственный оптимум (является и

глобальным). Если квадратичная форма является положительно определённой, то она является выпуклой функцией: при любом 0 1 имеет место:

f x1 1 x2 f x1 1 f x2 .

Линейная форма – также выпуклая функция. Поэтому в рассматриваемом случае целевая функция будет выпуклой. Необходимые условия Куна – Таккера:

L x, λ 0,

если

xj 0,

 

 

 

 

 

j 1, n ,

 

 

 

 

0,

если

xj 0,

 

xj

 

 

 

 

 

 

 

 

 

L x, λ 0,

если

i 0,

 

 

 

 

 

 

i 1, m

 

 

 

 

0,

если

i 0,

 

i

 

 

 

 

 

 

 

 

 

будут и достаточными условиями существования единственного оптимума.

Эти условия можно записать в виде:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

0,

p

 

 

0, x

 

p

 

0,

 

j

 

 

 

,

 

 

 

 

 

 

j

 

j

j

1, n

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

0,

q

 

0,

 

 

q

 

0,

 

 

 

 

 

 

 

 

 

 

 

 

i

 

i

 

i 1, m;

 

 

 

 

 

 

 

 

 

i

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L x, λ q x i fi x = cj xj d j k xj xk i

ai k xk bi

i 1

 

 

 

 

 

j

 

 

 

 

 

 

j

k

 

 

 

 

 

 

 

 

i

 

k

 

 

 

 

 

L x, λ

p

 

,

 

L

q ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xj

 

 

j

 

 

i

 

 

i

 

 

 

 

 

где L x, λ – функция Лагранжа.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Справедливо:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cj 2 d j k

xk i ai j

pj , j

 

;

 

 

 

 

 

1, n

 

 

 

 

 

 

k

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(13)

 

bi aik xk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qi , i 1, m.

 

 

 

 

 

 

 

 

 

 

k

В соответствии с условиями Куна – Таккера требуется найти решение этих уравнений при условиях:

 

 

 

xj

pj 0,

i qi 0,

xj

0;

(14)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i 0,

pj 0, qi 0; i 1, m,

j 1, n.

 

 

 

 

Классический

способ

решения этой задачи состоит в том, что уравнения

fi x 0, i

 

,

m n ,

используются для

исключения из рассмотрения m

1, m

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

q x1, x2 , , xn q1 y1, y2 , , yn m ,

где через y1, y2 , , yn m обозначены неисключённые переменные. Задача сводится к нахождению значений y1, y2 , , yn m , которые обращают в минимум функцию q1 и на

Regional architecture and engineering 2017

№2 41

СТРОИТЕЛЬНЫЕ МАТЕРИАЛЫ И ИЗДЕЛИЯ

которые не наложено никаких ограничений, то есть к задаче на безусловный экстремум.

Если ограничения имеют сложный вид, то исключение с их помощью m переменных из функции q x представляет значительные трудности. В связи с этим

задачи на условный экстремум и сводятся к задаче на безусловный экстремум с использованием функции Лагранжа.

В рассматриваемом случае система (13) содержит n + m уравнений с 2(n + m) переменными xj , i , pj , qi , из которых n + m являются свободными и могут быть

приравнены к нулю. Остальные переменные образуют при этом базисное решение, которое является допустимым, если выполняются условия (14).

Если число переменных в задаче невелико, то можно попытаться угадать допустимое базисное решение, положив произвольные n + m переменных свободными, приравняв их нулю, и, решив систему (14), найти значения базисных переменных. Однако нет никаких гарантий, что полученные значения переменных будут удовлетворять условиям (14). Поэтому попытки угадать допустимое базисное решение приходится проводить многократно.

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

новой базисной переменной необходимо проверять выполнение условий xj pj 0 ,i qi 0 : если в базисе имеется xj или i , то в него не может быть введено pj или qi

соответственно. При большом числе переменных угадать допустимое базисное решение становится чрезвычайно трудоёмким. В этом случае следует обратиться к известным эффективным систематическим методам получения допустимого базисного решения.

Для определения зависимости отклика от количественных факторов x и ошибокнаблюдения отклика используется регрессионная модель:

для u -го наблюдения:

y f x, ;

yu u y f xu , u ;

 

 

 

 

 

 

2

 

 

2

,

E u 0 ,

E u v 0, u v – для равноточных и

E u

 

 

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

функций от факторов, то получится модель регрессионного анализа, линейная по параметрам (линейная модель):

y 1 f1 x1, x2 , , xk 2 f2 x1, x2 , , xk m fm x1, x2 , , xk ,

m

yi fi x ,

i 1

y f т x ;

i – параметры модели (коэффициенты регрессии) i 1,2, , m ; fi x1, x2 , , xk – известные базисные функции переменных x1, x2 , , xk (факторы), не зависящие от параметров модели; f т x f1 x , f2 x , , fm x – вектор-

строка базисных функций (базисная вектор-функция); T 1, 2 , , m – вектор параметров модели.

42 Региональная архитектура и строительство 2017

№2

BUILDING MATERIALS AND PRODUCTS

Частными случаями полиномиальной модели регрессионного анализа (полином по факторам) являются:

– модель регрессионного анализа первого порядка (линейная модель), задавае-

мая полиномом первого порядка:

y 0 1x1 2 x2 k xk

(введя фиктивную переменную x0 1, модель можно представить в виде

k

yi xi );

i 1

модель регрессионного анализа второго порядка (квадратичная модель), зада-

k 1 k 2

ваемая полиномом второго порядка; (в общем случае содержит

 

2

параметров):

 

y 0 1x1 k xk 12 x1x2

k 1,k xk 1xk 11x12 k ,k xk2 .

Если отклик определяется качественными факторами (например, вид материала) и ошибками наблюдений отклика, то рассматривается модель дисперсионного анализа

y 0 1x1 2 x2 k xk ;

xk – дискретные переменные, обычно целочисленные (часто xi либо 0, либо 1). Обыч-

но если количественный фактор принимает небольшое число различных значений в эксперименте, то он рассматривается как качественный. При детерминированных неизвестных параметрах получается модель с постоянными факторами (модель I);

если все i (может быть за исключением одного) являются случайными величинами,

то – модель со случайными факторами (модель II); в промежуточных случаях – смешанная модель. При воздействии на отклик как количественных, так и каче-

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

Ограниченность использования указанных моделей при синтезе композиционных материалов связана с необходимостью установления аналитических зависимостей между параметрами моделей и рецептурно-технологическими факторами: параметры моделей являются лишь интегральными характеристиками [6] такой связи.

Список литературы

1.Garkina, I.A. Methodological principles design of composite materials / I.A. Garkina, A.M. Danilov, E.A. Budylina // Journal of Engineering and Applied Sciences. – 2016. – Vol.11(11). – P. 2524–2527.

2.Garkina, I.A. Materials as complex systems / I.A. Garkina, A.M. Danilov, V.P. Selyaev // Journal of Engineering and Applied Sciences. – 2016. – Vol. 11(11). – P. 2461–2464.

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. Building materials: the modelling, quality functional / I.A. Garkina // J. Ponte. – Mar, 2017. – Vol. 73. – Issue 3. – P.111–116.

5.Garkina, I.A. Parametric identification and optimization of properties of building materials as complex systems / I.A. Garkina, A.M. Danilov / J. Ponte. – Feb. 2017. – Vol. 73. – Issue 2. – P.119–125.

Regional architecture and engineering 2017

№2 43

СТРОИТЕЛЬНЫЕ МАТЕРИАЛЫ И ИЗДЕЛИЯ

6. Сизиков, В.П. Кризис системной методологии: с позиций синтеза систем / В.П. Сизиков, В.И. Разумов // Труды V Международной конференции «Идентификация систем и задачи управления» SICPRO’06. – М.: Институт проблем управления им. В.А.Трапезникова РАН, 2006. – С.1818–1860.

References

1.Garkina, I.A. Methodological principles design of composite materials / I.A. Garkina, A.M. Danilov, E.A. Budylina // Journal of Engineering and Applied Sciences. – 2016. – Vol.11(11). – P. 2524–2527.

2.Garkina, I.A. Materials as complex systems / I.A. Garkina, A.M. Danilov, V.P. Selyaev // Journal of Engineering and Applied Sciences. – 2016. – Vol. 11(11). – P. 2461–2464.

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. Building materials: the modelling, quality functional / I.A. Garkina // J. Ponte. – Mar, 2017. – Vol. 73. – Issue 3. – P.111–116.

5.Garkina, I.A. Parametric identification and optimization of properties of building materials as complex systems / I.A. Garkina, A.M. Danilov / J. Ponte. – Feb. 2017. – Vol. 73. – Issue 2. – P.119–125.

6.Sizikov, V.P. Crisis of system methodology: from positions of systems synthesis / V.P. Sizikov, V.I. Razumov // Proceedings of the V International Conference «System Identification and Control Problems» SICPRO'06. – M.: Institute for Control Sciences. V.A.Trapeznikov of the Russian Academy of Sciences, 2006. – P.1818–1860.

44 Региональная архитектура и строительство 2017

№2

BUILDING MATERIALS AND PRODUCTS

УДК 519.7:691

Пензенский государственный университет архитектуры и строительства

Россия, 440028, г. Пенза, ул. Германа Титова, д.28,

òåë.: (8412) 48-27-37; ôàêñ: (8421) 48-74-77

Данилов Александр Максимович,

доктор технических наук, профессор, советник РААСН, зав. кафедрой «Математика и математическое моделирование»

E-mail: fmatem@pguas.ru

Гарькина Ирина Александровна, доктор технических наук, профессор

кафедры «Математика и математическое моделирование»

E-mail: fmatem@pguas.ru

Penza State University of Architecture and Construction

Russia, 440028, Penza, 28, German Titov St., tel.: (8412) 48-27-37; fax: (8412) 48-74-77

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

Garkina Irina Aleksandrovna,

Doctor of Sciences, Professor of the department «Mathematics and Mathematical Modeling» E-mail: fmatem@pguas.ru

КОМПОЗИТЫ: ГРАФОАНАЛИТИЧЕСКИЙ МЕТОД ПАРАМЕТРИЧЕСКОЙ ИДЕНТИФИКАЦИИ КИНЕТИЧЕСКИХ ПРОЦЕССОВ

А.М. Данилов, И.А. Гарькина

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

Ключевые слова: композиты, формирование свойств, кинетические процессы, параметрическая идентификация, графоаналитический метод

COMPOSITES: GRAPHIC-ANALYTICAL METHOD OF PARAMETRICAL IDENTIFICATION OF KINETIC PROCESSES

A.M. Danilov, I.A. Garkina

The analytical synthesis of composites based on retrospective parametrical identification of kinetic processes with the use of a graphic-analytical method is given.

Keywords: composites, formation of properties, kinetic processes, parametrical identification, graphic-analytical method

Рассмотрим дисперсные системы, для кинетических процессов в которых характерно наличие точки перегиба; описываются обыкновенным дифференциальным уравнением (учитывается скорость и ускорение изменения контролируемого пара-

метра) [1…3]. Поэтому определение абсциссы точки перегиба является важнейшим элементом параметрической идентификации.

В отклонениях от равновесного состояния x xm будем иметь:

z 2 n z 02 z 0

n 0

( W p T1 p 1 k T2 p 1 ) .

Пусть k1,2 1,2 – корни характеристического уравнения k2 2 n k 02 0 .

Regional architecture and engineering 2017

№2 45

СТРОИТЕЛЬНЫЕ МАТЕРИАЛЫ И ИЗДЕЛИЯ

При n2 02 0 имеем

 

 

 

 

 

 

 

z c e 1t c e 2t ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

n n2 2 ,

2

n n2 2 ,

2

0.

1

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

1

 

 

При начальных условиях z 0 xm ,

 

 

z 0 0 будем иметь:

 

z

 

 

2

 

 

x e 1t

 

 

 

 

1

 

 

x e 2t ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

m

 

 

 

 

 

2

 

m

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

xm

 

 

2 e 1 t 1 e 2 t xm .

 

 

 

1 2

 

 

 

Из x 0 следует

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

x e

t

 

 

 

2

 

 

x e

 

t

0 ,

 

 

 

 

2 1

 

 

 

 

1

 

 

 

 

1

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

m

 

 

 

 

 

 

 

2

 

 

m

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

t

n

 

e

 

t

 

 

 

 

 

e

 

t

 

 

 

 

1

2

 

2

 

n ,

 

 

 

 

1

 

1

 

2

 

n .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так что точке перегиба соответствует значение t tn , определяемое из условия

 

 

 

 

 

 

 

 

t

n

 

 

 

 

 

1

 

 

 

ln

1

 

 

 

 

 

 

 

(2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(при t tn вогнутость сменяется на выпуклость).

 

 

Полученное соотношение устанавливает

связьабсциссыточкиперегибаскорнямихарактеристическогополинома(корниразличные).

Приведем методику определения 1

и

2 по экспериментально полученным

значениям

процесса x t . Из

2

1 следует,

 

что

 

составляющая 2

 

xm

 

e 1 t

 

 

 

 

2

процесса x t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

затухает быстрее,

 

 

чем аналогичная составляющая, соответствующая

корню 2.

Поэтому практически значение 2 можно определить по концу экспери-

ментально

полученного

процесса

 

x t .

Без

ограничения общности

рассуждений

можно принять

xm 1

(равносильно масштабированию x t ). При

значениях t ,

превышающих некоторого значения t1 , справедливо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 x t

 

 

 

e

2 t

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

2

 

 

 

1

 

 

 

 

 

 

 

1 e 2 t 2 e 1 t

,

 

e 1 2 t ,

ln

1 2 t .

 

 

 

 

 

 

 

 

 

 

 

Откуда

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln

t , tn t,

 

 

t tn .

 

 

 

 

 

 

 

 

1 2

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 x t

 

 

e

2 t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

 

 

 

 

 

 

 

 

 

при t tn . Введем

y t 1 x t Ae 2 t .

Тогда

y t T Ae 2 t T .

46 Региональная архитектура и строительство 2017

№2

BUILDING MATERIALS AND PRODUCTS

Откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y t

 

Ae 2 t

 

T

 

 

 

 

 

 

 

 

ln

 

 

 

 

 

e 2

 

 

,

2

T ln

, 2

 

T .

y t T

Ae 2 t T

 

 

Далее из x 0 следует

 

 

e 1 tn

 

 

e 2 t ,

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r e r 2 t e 2 t ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

1

1

,

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r e r 1 2 t ,

ln r r 1 2 t,

 

 

1

 

 

 

 

 

 

1

 

 

 

 

ln r

2 t,

ln r

r 1

2 t.

 

r 1

Откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 tn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r r 1 e

.

 

 

 

 

 

 

(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим функцию

Имеем

 

y

1

r 1 ln r

 

r

 

 

y

 

 

 

r 1 2

Отсюда y 0

при r 1.

Справедливо

 

 

 

1

 

 

 

y r

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

r 1

 

 

 

 

 

 

 

 

 

 

ln y

1

 

ln r ,

 

 

 

 

 

 

 

 

 

r 1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

ln r

 

 

 

 

 

 

ln r

 

 

1

r

 

 

1 1

r

 

 

 

 

 

 

 

 

 

 

, y

r

 

 

 

 

 

 

.

 

 

 

 

 

r 1

 

 

 

r 1 2

 

 

r 1 2

 

 

 

 

 

 

 

 

 

 

 

1

 

lim 1

 

1

 

 

 

 

 

lim r r 1

e,

 

r 1

 

 

 

 

0

 

 

 

 

 

 

 

 

r 1.

 

 

 

 

 

 

 

1

 

 

 

 

Так что в интервале (1, )

y r

r 1

не превышает e. Поэтому уравнение (1)

имеет решение r 1 лишь при e 2 tn e .

Откуда следует

2 tn 1,

и 2 должно удовлетворять условию

 

 

 

 

 

 

 

 

t

 

 

 

1

.

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

При этом 2max

n (тогда 1

2 ,

 

 

n 0 ).

Справедливо

 

 

 

 

 

 

1

 

 

 

 

1

 

 

 

 

t

 

 

 

 

ln

,

 

2

n

 

1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

Из

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln r

 

 

 

 

d

 

1 r

 

 

 

 

 

d r

 

 

 

 

r 1 2

 

ln r

.

 

 

r 1

< 0

следует, что 2 tn с ростом r уменьшается.

Regional architecture and engineering 2017

№2 47

СТРОИТЕЛЬНЫЕ МАТЕРИАЛЫ И ИЗДЕЛИЯ

Отметим,

 

1

 

 

 

d y

 

d

 

r r 1

.

d r

d r

 

 

 

 

Вид функции r r , полученный аппроксимацией табличных значений решений

уравнения (1) при различных 2 tn методом наименьших квадратов, приводится на рис.1.

r

45

40

35

r(v) 11376,v 1.5518

30

25

20

15

10

5

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

0,9 v

0

0,1

0,2

0,3

0,4

0,5

0,6

0,7

0,8

Рис. 1. Вид функции r = r( )

Определим зависимость корней 1, 2 (определяют вид кинетического процесса)

от параметров модели 0 и

n (определяют упругие и демпфирующие свойства

материала).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n, n2 02

 

1 n n2 02

2 n,

2 n n2 02

 

следует

 

 

 

 

 

 

 

 

 

 

 

2 n 1 2 n .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При этом 1 2 при

n 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Справедливо

 

 

 

 

 

 

n

 

 

 

 

 

2 1

 

n

 

 

 

 

 

 

1

1

 

 

 

 

 

 

0,

 

 

 

0,

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

n2 02

n

 

 

 

 

n2 02

 

 

1

 

 

 

0

 

 

 

 

0,

 

2

 

 

 

0

 

0,

 

 

 

 

 

 

 

 

 

 

n2 2

 

 

 

0

 

 

 

n2 2

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

0,

2

2

 

 

2

 

 

 

 

0,

 

1

 

 

 

0

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n2

 

 

 

n2 02 3

 

n2

n2 02 3

 

2

 

 

 

 

 

2

 

0,

2

2

 

 

 

2

 

 

 

 

0.

 

1

 

 

 

 

0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

02

 

 

 

n2 02 3

2

 

 

 

 

n2 02 3

 

Введем безразмерный коэффициент демпфирования

n

,

n ,

1. Его

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

величина определяется структурой и физико-химическими свойствами материала.

Имеем

0 ,

2 2 1 0 ; 1 2 2 0

 

1 2 1

2 1;

48 Региональная архитектура и строительство 2017

№2

 

 

 

 

 

 

 

 

 

 

 

 

 

BUILDING MATERIALS AND PRODUCTS

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

2 1

 

 

 

 

2 1

 

 

 

 

 

 

2 1

 

 

 

 

2

 

 

 

 

 

2

 

1

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1 ;

2 1

2 1 2 1

 

 

 

02

 

 

 

 

 

 

 

 

 

 

 

r;

 

 

 

 

 

 

 

 

 

r;

 

 

 

 

 

 

 

1

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

Справедливо

 

 

 

 

 

 

1

 

 

 

0

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r 2

 

 

 

 

 

 

2

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

1

 

 

 

 

 

1

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

2

 

2 1 2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

r ;

 

r

0

 

 

 

 

 

1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 1

 

 

 

 

 

 

 

2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Как видим, с ростом значение

1

 

 

растет.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

r

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

2 r

 

 

 

 

 

 

 

2 2 1 ;

Из

 

 

 

2

 

2

1 r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

2

 

1

 

 

 

 

2

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

2r

 

 

0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 2 1;

 

 

 

3 2 4;

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

2

 

 

 

Так что функция r r

 

имеет перегиб в точке с абсциссой

 

 

 

 

, при этом

 

 

3

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при 1

 

функция выпукла, при

 

 

 

– вогнута.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Легко заметить, что функция

 

 

 

 

 

является решением дифференциального

уравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d r

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

2 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d r

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

,

 

 

 

1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

2

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln

 

 

 

 

 

 

,

 

 

 

ln

 

 

2ln

2 1

ln c,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r c 2 1 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При начальном условии r 0 0

 

 

 

02

 

 

1 2

получим: с = 1, откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r 2 1 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Regional architecture and engineering 2017

№2 49

СТРОИТЕЛЬНЫЕ МАТЕРИАЛЫ И ИЗДЕЛИЯ

Далее рассмотрим связь между характеристиками кинетических процессов и параметрами модели.

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

перегиба и параметрами , n модели (или

n

и ). Справедливо

 

0

 

 

 

 

 

0

 

0

 

 

 

 

 

 

 

 

 

 

tn

 

 

1

ln

n

n2

02

 

 

 

 

 

 

,

 

n2 2

n

n2

 

 

2

 

2

 

 

 

0

 

 

 

 

0

 

откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

n

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

tn

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

ln

 

 

 

 

 

 

 

n

 

n

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

2

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

n2 02

 

 

 

n

 

2

 

 

 

 

 

 

 

 

 

n

n2 02

 

 

 

0

 

0

 

 

 

 

 

 

 

 

 

2

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln r 2 2 1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 n

 

 

 

 

 

 

 

 

 

 

n

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В силу

 

n

2

 

2

отсюда следует,

что знак

 

tn

 

 

определяется знаком функции

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y ln r 2

 

2 1 .

Так как

 

 

y 0

 

 

1,

то значение

t

n

убывает с ростом

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

( 0 ,

n

1 ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Также справедливо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tn

 

 

 

 

 

 

n

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

n

 

n

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 n

 

 

0

 

ln

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

2 n2 02

3

 

 

 

 

 

 

n

 

 

 

 

n

 

 

n

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

1

 

ln r

 

 

 

 

 

 

 

 

 

 

 

 

 

n

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 1

 

 

0

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

2

n2 02

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда следует, что с ростом

 

n

 

значение

 

 

tn

 

уменьшается,

точка

перегиба

смещается влево.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 и

n

Таким

образом,

 

для увеличения

 

 

значения

 

следует

 

уменьшить

(предполагается выполнение условия n 0 ).

Приразработкематериаловследуетучитыватьзависимости 2 иr от , 0 , n .

Имеем

Справедливо

Отсюда с учетом

 

 

r r ,

 

 

 

n

0 , n .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

r

 

d r

 

 

 

 

 

d r

 

 

n

 

 

 

 

 

 

 

 

 

 

.

 

 

 

0

d 0

d

2

 

 

 

 

 

 

 

0

 

d r

0

 

1

следует

 

 

r

0

1.

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

50 Региональная архитектура и строительство 2017

№2

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]