Оптимизация при наличии ограничений.

 

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

Ограничение диапазона варьирования переменных связано, прежде всего, с тем, что в условиях реального технологического процесса (или объекта) наряду с целевой функцией имеется ряд контролируемых выходных параметров (где L = II,III,…, G), также зависящих от k -мерного вектора управляющих воздействий . Обычно аналитическая форма этих зависимостей, как и характер целевой функции , априори неизвестны и значения и в каждой точке факторного пространства могут быть определены только экспериментально. На значения контролируемых входных параметров обычно накладываются ограничения

; L = II,III,…, G, (1)

; L = II,III,…, G.

Входные параметры x_i также имеют ограничения

, i = 1,2,…, k. (2)

Ограничения (1) и (2) выделяют в факторном пространстве область допустимых значений вектора управляющих параметров . На рисунке 27 приведен пример формирования области допустимых значений для двухфакторного процесса.

 

Dk

yw=const

yIV=const

yIII=const

yII=const

yI=const

x 2

x 1

x 2max

x 2min

x 1max

x 1min

 

Рис.27. Формирование области допустимых значений

для двухфакторного процесса.

 

Концентрические кривые показывают линии равных значений целевой функции, т.е. . Так как ограничения (1) задаются в виде неравенства , то в факторном пространстве им тоже соответствуют поверхности равных значений ; L = II,III,…, G (кривые со штриховкой на рис. 27). Поверхности ограничения, пересекаясь, выделяют в факторном пространстве область D_k допустимых значений со сложной конфигурацией.

Экстремальная точка может лежать не только внутри, но и за пределами допустимой области. В этом случае задача сводится к отысканию точки условного экстремума функции у_I на границе области допустимых значений D_k.

Задача оптимизации при наличии ограничений формулируется следующим образом: необходимо найти координаты максимума целевой функции на подмножестве всех , принадлежащих факторному пространству переменных процесса и удовлетворяющих условиям

; ; L = II,III,…, G (3)

 

т. е. найти .

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

Рассмотрим процедуру эксперимента для нахождения точки в факторном пространстве двух переменных х₁ и х₂, отвечающую условиям

(4)

(5)

(6)

(7)

Ограничения (5) и (6) наложены непосредственно на значения варьируемых параметров, а ограничения (7) – на значения выходного показателя , которые определяются только из опыта.

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

Статистическая обработка результатов ПФЭ также выполняется по стандартной схеме, т. е. проверяется воспроизводимость опытов, значимость коэффициентов b ₀, b ₁, b ₂ уравнения регрессии, адекватность линейной модели

(8)

Если при осуществлении факторного эксперимента ни в одной из пробных точек ограничения (5), (6) и (7) не были нарушены, то из центра планирования проводится крутое восхождение в направлении .

Для того чтобы величина градиента не зависела от скорости возрастания функции в точке , его составляющие следует пронормировать

; (9)

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

Крутое восхождение прекращается в следующих случаях

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

2. Если нарушается ограничение типа (7), тогда для корректировки направления движения ставится новый факторный эксперимент с центром в точке .

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

, (10)

где – вектор, нормальный к эквипотенциальной поверхности (поверхности равного уровня) в точке и ведет в область допустимых значений (рис. 28).

 

 

Рис.28. Выбор компромисного направления.

 

Для формирования вектора наряду с гиперплоскостью (8) находится также уравнение линейной аппроксимации

(11)

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

(12)

Если в окрестностях точки нарушается ограничение

,

то вектор направлен внутрь допустимой области, если же нарушено условие

,

то он имеет направление из допустимой области.

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

если если если

 

Следовательно, координаты искомого вектора будут

.

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

Если после отбрасывания одной из составляющих и движения вдоль границы точка выходит на ограничение типа (7) или попадает в область, где одновременно выполняются условия (5) и (6), то снова проводится корректировка направления с помощью ПФЭ.

В процессе экспериментальной оптимизации точка не должна подходить к границам типа (5) и (6) ближе величины шага варьирования , т.е. практически вместо неравенств типа (5) и (6) следует использовать неравенства типа

, i = 1,2,…, k

Результатом применения описанной методики могут быть четыре случая, служащих сигналом окончания поиска:

1) В одном из очередных экспериментов , т. е. точка находится в окрестностях явного экстремума целевой функции.

2) и вектор . Этот случай соответствует условному экстремуму y_i на границе типа (7)

3) В некоторой точке нарушаются оба условия (5) и (6) и вектор ведет из допустимой области (условный экстремум находится в углу допустимой области)

4) После корректировки направления и точек, где одновременно нарушались ограничения (5), (6) и (7) уровень выхода y_I начинает уменьшаться (или увеличиваться). Условный экстремум находится в углу допустимой области Эти четыре случая указывают на то, что экстремальное значение целевой функции достигнуто. Однако в связи с тем, что при описанном выше алгоритме метода возможен выход на локальный экстремум, рекомендуется повторить поиске новой исходной точки.

 

Литература.

 

 

1. Алабужев Л.М. Теория подобия и размерностей / Л.М.Алабужев, В.Б.Геронибус, Л.М.Минкевич и др. – М.: Высшая школа, 1985.

2. Годлевский В.А. Введение в анализ экспериментальных данных. Учебное пособие – Иваново, 1993.

3. Гречишников В.А. Основы научных исследований в области проектирования и эксплуатации режущего инструмента /В.А.Гречишников, Н.В.Колесов, Е.В.Козлов.– М.:МОССТАНКИН, 1990.

4. Захарцев С.Н. Математическая статистика и планирование эксперимента в технологии машиностроения, ч.I,II – М. Издат. МГТУ им. Баумана, 1991.

5. Кане М.М. Основы научных исследований в технологии машиностроения, – Минск, Высшая школа, - 1987.

6. Кацев П.Г. Статистические методы исследования режущего инструмента. – М.,Машиностроение, - 1974.

7. Круг Г.К. Теоретические основы планирования экспериментальных исследований. Лабораторный практикум – М. МЭИ, 1969.

8. Крутов В.И. Основы научных исследований /В.И.Крутов, И.М.Грушко, В.В.Попов и др. – М.: Высшая школа, 1989.

9. Матевосов Л.М. Охрана промышленной собственности – М.: ИНИЦ Роспатента, 2003.

10. Налимов В.В., Статистические методы планирования экспериментов / В.В.Налимов, Н.А.Дернова – М.: Наука, 1995.

11. Селиванов С.Г. Инноватика /С.Г.Селиванов, М.Б.Гузаиров, А.А.Кутин. – М.: Машиностр., 2008. – 721с.

12. Силин С.С. Теория подобия в приложении к технологии машиностроения. – Ярославль, ЯПИ, 1989.

13. Спирин Н.А. Методы планирования и обработки результатов инженерного эксперимента /Н.А.Спирин, В.В.Лавров, А.Н.Бондин и др. – Екатеринбург, 2003.

14. Шаламов В.Г. Математическое моделирование при резании металлов. – Челябинск, Челяб.ПИ, 1995.

15. Ящерицын П.И. Планирование эксперимента в машиностроении /П.И.Ящерицын, Е.И.Мехаринский, Справочное пособие. – Минск, Высшая школа, 1985.

 

Приложение 1

 

Информационно-поисковые системы и БД
российских и зарубежных патентных ведомств

 

1. Россия – данные WEB-сайт является частью WEB-сайта Роспатента и предназначен для обеспечения доступа к российским базам данных по изобретениям, полезным моделям и товарным знакам. Открыт свободный доступ к БД «Рефераты российских заявок и патентов на изобретения на русском языке 1994-2000гг.» (RUABEN), «Рефераты российских патентов на изобретения на английском языке 1994-2000гг.» (RUABEN). Полнотекстовые БД российских патентов на изобретения (RUPAT), БД полезных моделей (RUABU1) и БД товарных знаков (RUTM и WRUTM) предоставляются в доступ за плату. Для осуществления поиска используются два интерфейса: полный и сокращенный.

2. PCT – Patent Cooperation Treaty (PCT) System – содержит патенты с 1997г. по настоящее время, позволяет проводить бесплатный поиск и просмотр титульной страницы. Сайт принадлежит Всемирной организации интеллектуальной собственности (WIPO) http://petgazette.wipo.int/ и содержит библиографические данные, рефераты и рисунки патентов, опубликованных с 1 января 1997г.

3. GB and European patents – сайт http://gb.espacenet.com/ принадлежит Патентному бюро Великобритании и содержит библиографические данные и рисунки патентов Великобритании, опубликованные Патентным бюро Великобритании с 5 июня 1996г. База данных содержит также европейские патенты, мировые патенты (база данных PCT), японские патенты и патенты других стран. Присмотр рисунков возможен с использованием Adobe Acrobat 3.x.

4. Europ’s Network Patent Databases – esp@cent ­– сайт Европейского Патентного Бюро, содержит полные тексты европейских патентов с октября 1997г., а также национальные патенты европейских стран, мировые патенты, опубликованные Всемирной организации интеллектуальной собственности (WIPO) с ноября 1997г., японские патенты с 1980г. и приложения патентов, опубликованных в более чем 50 странах мира. Сайт http://ep.espacenet/com/.

5. IBM Patent Search – поддерживается IBM – International Business Machines Corporation, содержит патенты США с 1971г. и неполные данные по патентам предыдущих годов, патенты Европейского Патентного Бюро (EPO) с 1979г. и патенты Patent Cooperation Treaty (PCT) с 1997г. Сайт http://www.patents.ibm.com/, на нем доступны заявки и титульные страницы.

6. USPTO Web Patent Database – содержит полные тексты и рисунки патентов США с 1976г. по настоящее время. Бесплатный поиск и доступ к полному тексту может быть проведен по номеру патента, ключевым словам названия и реферата, фамилии и другим полям. Сайт http://www.uspto.gov/patft/index.html поддерживается USPTO – Бюро Патентов и Торговых марок США.

7. US Patent and Trademark Office Official Gazette – официальный бюллетень Патентов и Торговых марок США, публикуется еженедельно с 1964г. по настоящее время для представления патентов, которые будут опубликованы, и торговых марок, которые будут зарегистрированы. Сайт http://www.uspto.gov/web/offices/com/sol/og/ позволяет проводить бесплатный поиск и просмотр.

8. PatentMiner – предоставляет возможность бесплатного поиска (регистрация обязательна) патентной библиографической информации, просмотра PDF-рисунков патентов США, поубликованных с 1970г. Эта БД поддерживает связь с рисунками патентов БД IBM Patent imag. и позволяет получить PDF-изображения патентов, включая рисунки (услуга платная). Сайт http://www.patentminer.com/LIVE/cgi-bin/pm.cgi управляется Manning and Napier Information Services.

9. STO’s Internet Patent Search System – база данных с 1970г. по настоящее время, поддерживается системой Source Translation and Optimization (STO), которая может быть использована для поиска по БД патентов и торговых марок США по всем годам с использованием классификационных классов USPTO или по номеру патента. Возможен просмотр рефератов патентов, опубликованных в период 1981-1989гг. Сайт http://metalab.unc.edu/patents/intropat.html.

10. Micro Patent WEB – он-лайновый бюллетень США (официальный бюллетень USPTO), который является частью сайта Micro Patent WEB, содержит полную суммарную информацию, включая рефераты, для каждого из 2500 новых патентов США, регистрируемых каждую неделю. Сайт http://www.micropat.com/og/front.html (платная база данных, но доступ к он-лайновому бюллетеню бесплатный).

11. БД патентов Великобритании – содержит библиографические данные и иллюстрации патентов Великобритании, опубликованных Патентным Бюро с 5 июня 1996г. В БД включены также национальные патенты государств – членов Европейского Патентного Бюро, европейские патенты, японские патенты, мировые патенты (БД PCT).

 

 

Приложение 2

 

2.1. Задание

 

По результатам полного факторного эксперимента 2³ при количестве параллельных опытов n₀ = 3 построить математическую модель для методов и условий обработки, приведенных ниже. Исследовать математическую модель на однородность дисперсий, статистическую значимость коэффициентов матмодели (уравнения регрессии), адекватность результатов эксперимента и матмодели при достоверности P_D = 0,95.

1. Эмпирическую формулу для расчета скорости резания при сверлении стали 09X18 спиральными сверлами диаметром D =5...15 мм на глубину l =3D; 5D; 7D; 10D с подачами S = 0,05...0,15мм/об, обеспечивающей стойкость Т=20...40мин. Матмодель представить в следующем виде:

 

.

 

Результаты экспериментов для четырех вариантов обработки (различных «l») приведены в таблице планирования (таблица 1).

 

Таблица 1

 

 

 

 

№ опыта	Факторы	Значение V, м/мин
х 1 (D)	х2 (Т)	х3 (S)		1) l = 3D	2) l = 5D	3) l = 7D	4) l = 10D
V1	V2	V3	V1	V2	V3	V1	V2	V3	V1	V2	V3
	+	+	+	6,8	7,1	7,4	4,4	4,7	4,8	4,1	3,7	4,3	3,6	2,9	3,3
	–	+	+	4,1	4,8	4,6	3,4	2,7	3,1	2,6	2,9	2,2	1,8	2,2	2,5
	+	–	+	8,6	9,1	7,7	5,1	5,6	6,3	4,7	5,3	4,3	3,5	4,0	4,4
	–	–	+	5,5	5,2	4,9	3,1	3,3	3,8	3,0	2,6	3,2	2,1	2,5	2,8
	+	+	–	16,8	19,3	21,2	10,9	13,9	12,6	9,4	11,0	11,9	7,7	8,7	9,7
	–	+	–	11,0	12,3	11,5	7,7	7,0	8,2	6,1	6,5	7,0	5,8	5,0	5,3
	+	–	–	24,2	22,5	20,2	14,8	12,9	15,9	12,5	13,3	11,3	9,2	11,0	9,6
	–	–	–	13,0	14,1	13,6	8,9	9,4	8,5	7,1	8,0	7,6	5,9	6,3	6,6

 

2. Эмпирическую формулу для расчета скорости резания при продольном точении стали 10Х12НВМА резцами с главным углом в плане φ = 30°; 45°; 60°; 90° с глубиной резания t = 1...3 мм; подачей S = 0,1...0,33 мм/об, стойкостью Т = 30...90 мин. Матмодель представить в следующем виде:

.

 

Результаты экспериментов приведены в матрице планирования (таблица 2).

 

 

Таблица 2

 

 

 

 

№ опыта

Факторы

Значение V, м/мин

х 1 (T)

х2 (S)

х3 (t)

1) φ = 30°

2) φ = 45°

3) φ = 60°

4) φ = 90°

V1

V2

V3

V1

V2

V3

V1

V2

V3

V1

V2

V3

+

+

+

96,0

–

+

+

+

–

+

–

–

+

+

+

–

82.

–

+

–

+

–

–

–

–

–

 

3. Эмпирическую формулу для расчета скорости резания при продольном точении никелевого сплава ХН78Т резцами с главным углом в плане φ = 45°; 60°; 75°; 90° с глубиной резания t = 1...3 мм; подачей S = 0,1…0,3 мм/об, стойкостью Т = 30...90мин. Матмодель представить в виде формулы (2). Результаты экспериментов приведены в матрице планирования (таблица 3).

 

Таблица 3

 

 

 

 

 

№ опыта	Факторы	Значение V, м/мин *
х 1 (T)	х2 (S)	х3 (t)	1) φ = 45°	2) φ = 60°	3) φ = 75°	4) φ = 90°
V1	V2	V3	V1	V2	V3	V1	V2	V3	V1	V2	V3
	+	+	+	9,2	8,1	10,4	8,3	7,3	9,4	7,8	6,9	8,9	7,3	6,4	8,5
	–	+	+	9,3	10,6	12,0	8,4	9,7	11,0	7,9	8,4	10,3	7,4	8,5	9,6
	+	–	+	13,1	10,2	11,6	12,1	9,1	10,6	10,2	8,2	10,7	10,3	8,3	9,3
	–	–	+	11,8	15,1	13,4	10,8	13,7	12,2	10,1	11,5	12,8	9,4	11,9	10,7
	+	+	–	10,4	8,9	9,9	11,3	9,2	8,1	10,1	8,6	8,1	8,9	8,1	7,1
	–	+	–	13,4	10,7	10,6	10,8	12,2	9,5	10,2	8,9	11,4	9,5	10,7	8,4
	+	–	–	12,9	11,6	13,6	12,0	10,5	11,9	11,9	11,2	9,9	11,8	10,5	9,3
	–	–	–	17,0	15,4	12,2	15,6	14,0	10,7	14,6	13,2	10,7	10,8	13,7	12,3

 

4. Эмпирическую формулу для расчета шероховатости поверхности R_a при чистовом точении стали 45 (НВ=2000Мпа) резцом с пластинкой Т15К6; φ =45...75°, φ₁ = 25°, 30°, 35°; r = 0,5...1,5 мм; t = 1мм; подачей S = 0,08...0,18 мм/об; V = 50, 100, 150 м/мин. Для указанных значений φ₁ и V матмодель представить в следующем виде:

, (мкм)

Результаты экспериментов при разных значениях φ₁ и V приведены в таблицах 4, 5.

 

 

Таблица 4

 

 

 

 

 

№ опыта	Факторы	Значение Ra, мкм при V = 50 м/мин
х 1 (S)	х2 (φ)	х3 (r)	1) φ1 = 25°	2) φ1 = 30°	3) φ1 = 35°
Ra 1	Ra 2	Ra 3	Ra 1	Ra 2	Ra 3	Ra 1	Ra 2	Ra 3
	+		+	2,41	2,20	2,59	2,61	2,38	2,81	2,78	2,53	2,98
	–	+	+	1,34	1,52	1,67	1,45	1,66	1,82	1,55	1,76	1,93
	+	–	+	2,14	1,72	1,93	2,36	1,87	2,11	2,46	1,99	2,23
	–	–	+	1,10	1,21	1,33	1,19	1,32	1,45	1,27	1,40	1,54
	+	+	–	4,52	4,07	5,12	4,91	4,42	5,54	5,16	4,65	5,86
	–	+	–	3,12	2,81	2,51	3,38	3,08	2,79	3,58	3,24	2,91
	+	–	–	3,47	4,01	3,65	3,76	4,35	3,96	3,98	4,61	4,20
	–	–	–	2,53	2,31	2,10	2,75	2,52	2,31	2,91	2,68	2,42

 

Таблица 5

 

 

 

 

№ опыта	Факторы	Значение Ra, мкм при V = 100 м/мин
х 1 (S)	х2 (φ)	х3 (r)	1) φ1 = 25°	2) φ1 = 30°	3) φ1 = 35°
Ra 1	Ra 2	Ra 3	Ra 1	Ra 2	Ra 3	Ra 1	Ra 2	Ra 3
	+	+	+	2,30	2,11	2,47	2,49	2,27	2,68	2,65	2,41	2,83
	–	+	+	1,28	1,45	1,59	1,53	1,58	1,73	1,48	1,68	1,76
	+	–	+	2,04	1,64	1,84	2,25	1,78	2,01	2,34	1,91	2,12
	–	–	+	1,05	1,15	1,27	1,15	1,26	1,38	1,21	1,33	1,42
	+	+	–	4,31	3,88	4,87	4,68	4,21	5,28	4,91	4,48	5,61
	–	+	–	2,97	2,68	2,39	3,22	2,93	2,66	3,40	3,06	2,81
	+	–	–	3,31	3,84	3,48	3,58	4,14	3,77	3,79	4,39	4,01
	–	–	–	2,41	2,20	2,00	2,62	2,40	2,30	2,77	2,55	2,38

 

 

Таблица 6

 

 

 

 

 

№ опыта	Факторы	Значение Ra, мкм при V = 150 м/мин
х 1 (S)	х2 (φ)	х3 (r)	1) φ1 = 25°	2) φ1 = 30°	3) φ1 = 35°
Ra 1	Ra 2	Ra 3	Ra 1	Ra 2	Ra 3	Ra 1	Ra 2	Ra 3
	+	+	+	2,21	2,07	2,42	2,44	2,23	2,63	2,15	2,36	2,72
	–	+	+	1,25	1,41	1,55	1,50	1,55	1,69	1,43	1,65	1,72
	+	–	+	2,00	1,61	1,81	2,20	1,74	1,97	2,29	1,88	2,08
	–	–	+	1,03	1,13	1,24	1,13	1,23	1,35	1,19	1,30	1,39
	+	+	–	4,23	3,80	4,82	4,59	4,13	5,18	4,82	4,39	5,50
	–	+	–	2,91	2,63	2,34	3,16	2,87	2,61	3,34	3,01	2,76
	+	–	–	3,26	3,77	3,41	3,51	4,06	3,70	3,72	4,31	3,73
	–	–	–	2,36	2,16	1,96	2,57	2,35	2,25	2,72	2,49	2,33

 

5. Эмпирическую формулу для расчета шероховатости поверхности R_a при точении стали 45 (HRC 50) алмазным резцом с r = 0,5...2 мм, γ = +4°...–10°, S = 0,05...0,2 мм/об. Для скоростей резания V = 50; 100; 150 м/мин матмодель представить в следующем виде

, где 90­ – γ = δ – угол резания.

Результаты экспериментов при разных значениях V приведены в таблице 7.

Таблица 7

 

 

 

 

№ опыта	Факторы	Значение Ra, мкм
х 1 (S)	х2 (90–γ)	х3 (r)	1) V = 50 м/мин	2) V = 100 м/мин	3) V = 150 м/мин
Ra 1	Ra 2	Ra 3	Ra 1	Ra 2	Ra 3	Ra 1	Ra 2	Ra 3
	+	4-	+	0,513	0,528	0,491	0,451	0,463	0,431	0,418	0,428	0,401
	-	+	+	0,214	0,225	0,235	0,188	0,197	0,206	0,174	0,182	0,191
	+	-	+	0,321	0,352	0,336	0,282	0,309	0,295	0,261	0,286	0,273
	-	-	+	0,141	0,148	0,155	0,124	0,130	0,136	0,115	0,121	0,126
	+	+	-	0,767	0,798	0,736	0,672	0,700	0,648	0,622	0,648	0,600
	-	+	-	0,322	0,336	0,351	0,282	0,295	0,308	0,261	0,273	0,285
	+	-	-	0,483	0,528	0,503	0,424	0,463	0,441	0,393	0,429	0,408
	-	-	-	0,233	0,221	0,210	0,204	0,194	0,184	0,189	0,180	0,171

 

 

Приложение 2 (продолж.)

 

2.2. Критерий грубых ошибок Груббса β _max

 

Таблица 8

 

 

N изм.	β max при PD	N изм.	β max при PD
0,9	0,95	0,99	0,9	0,95	0,99
	1,41	1,41	1,41		2,33	2,49	2,81
	1,64	1,69	1,72		2,35	2,52	2,84
	1,79	1,87	1,96		2,38	2,55	2,87
	1,89	2,00	2,13		2,40	2,58	2,90
	1,97	2,09	2,26		2,43	2,60	2,93
	2,04	2,17	2;37		2,45	2,62	2,96
	2,10	2,24	2,46		2,54	2,72	3,07

⇐ Предыдущая123456

Поделиться:

Воспользуйтесь поиском по сайту: