
- •Содержание
- •1. Основные принципы теории планирования эксперимента
- •1.1 Основные понятия и определения
- •1.2 Классификация планов
- •1.3 Полный факторный эксперимент
- •1.3.1 Принятие решений перед началом эксперимента
- •1.3.2 Подготовка плана и матрицы эксперимента
- •1.3.3 Проведение эксперимента
- •1.3.4 Первичная статическая обработка результатов
- •1.3.5 Выявление грубых промахов
- •1.3.7 Оценка значимости коэффициентов
- •Проверка адекватности модели
- •1.4 Дробный факторный эксперимент
- •1.5 Планирование второго порядка
- •1.5.1 Область применения
- •1.5.2 Модель второго порядка
- •1.5.3 Выбор типа плана
- •1.5.4 Преобразование матрицы цкп к ортогональному виду
- •1.6 Отсеивающие эксперименты
- •1.7 Практическое применение математической модели для оценки свойств объекта
- •1.7.1 Качественный анализ модели
- •Применение модели для количественной оценки свойств объекта
- •2. Статистический анализ данных
- •2.1. Основные понятия и определения
- •2.2. Плотность вероятностей
- •2.3 Функция распределения
- •2.4 Корреляционный анализ
- •2.5 Регрессионный анализ
- •3. Оптимизация инженерных проектов
- •3.1 Методы проектирования и оптимизации
- •3.2 Пример практического применения методов однофакторной оптимизации
- •3.2.1 Описание объекта оптимизации
- •3.2.2 Цель оптимизации
- •3.2.3 Методика оптимизации
- •Выбор методов оптимизации
- •3.3 Пример практического применения методов многофакторной оптимизации
- •3.3.1 Описание объекта оптимизации
- •3.3.1 Цель оптимизации
- •3.3.2 Методика оптимизации
- •3.3.3 Метод Гаусса-Зайделя
- •3.3.4 Метод Бокса-Уилсона (градиентный метод)
- •3.3.5 Последовательный симплексный метод
- •3.4 Пример многокритериальной оптимизации
- •3.4.1 Описание объекта оптимизации
- •3.4.3 Предпосылки оптимизации
- •3.4.4 Выбор методов оптимизации
- •3.4.5 Решение оптимизационной задачи методом Парето
- •Рекомендуемая литература
3.4.5 Решение оптимизационной задачи методом Парето
Пусть первым критерием будет квадрат второй собственной частоты, а вторым – половина динамической массы машины:
К2=m;
Если перемножить критерии, получится уравнение семейства гипербол в пространстве критериев:
К1 · К2 = С1 + 2С0
Крайними
в этом семействе будут гиперболы,
соответствующие крайним значениям
жёсткостей С1
и С0:
и
Эти кривые образуют две границы области допустимых состояний объекта. Другие две границы образованы крайними значениями массы m.
Рассмотрим варианты расположения оптимальной точки.
Может ли быть оптимальной некоторая точка N, расположенная внутри области возможных состояний объекта АВСD?
Если спроектировать её влево до границы AD, получим точку P , которая будет предпочтительнее точки N, поскольку ей соответствует такое же значение критерия К2, но критерий К1 будет лучше.
Если спроектировать её вниз до границы AD, получим точку Q , которая тоже будет предпочтительнее точки N, поскольку ей соответствует такое же значение критерия К1, но критерий К2 будет лучше.
В точке R оба критерия будут лучше, чем в точке N.
Вывод: оптимальные решения лежат не внутри области возможных состояний объекта, а на одной из его границ.
При заданном направлении оптимизации (К1 → min; К2 → min) все оптимальные точки расположены на границе AD.
В каком именно месте этой границы? Ответ на этот вопрос зависит от степени предпочтения (сравнительной важности) критериев:
- если снижение частоты важнее, чем уменьшение массы – оптимальная точка должна располагаться ближе к точке A.
- если уменьшение массы важнее, чем снижение частоты – оптимальная точка должна располагаться ближе к точке D.
Если степень предпочтений задана (например, методом экспертных оценок установлено, что уменьшение частоты в 2 раза важнее, чем уменьшение массы), можно определить точное расположение оптимальной точки следующим образом.
Допустим, мы перемещаемся вдоль границы области возможных состояний объекта из точки 1 в точку 2, и при этом частота уменьшается на 10%. Поскольку критерии противоречивы, это перемещение вызовет увеличение массы, причём в соответствии с принятой степенью важности можно допустить, что оно будет составлять 20%.
Это
условие можно выразить пропорцией:
из которой следует: m 1+2ω1= m 2+2 ω2
Геометрически это соответствует семейству прямых:
m + 2ω = idem. Если мы выберем ту из них, которая является касательной к границе области возможных состояний объекта, мы получим решение, соответствующее принятой степени предпочтения критериев (на рисунке это точка F).
Резюме: практическое применение оптимизирующих методик в процессе проектирования и испытаний новой техники требует от специалистов дополнительных теоретических знаний и практических навыков, однако игра стоит свеч, поскольку оптимизация характеристик способствует повышению качества и конкурентоспособности оборудования.
Рекомендуемая литература
Алёшин Н.В., Ляховицкий А.Г., Царёв Б.А. Методология инженерной и научной деятельности в морской технике. Учебное пособие. Изд. СПбГМТУ, 2000; 294 с.
Половинкин А.И. Основы инженерного творчества. М., «Машиностроение», 1988; 361 стр.
Диксон Дж. Проектирование систем: изобретательство, анализ и принятие решений. М.,«Мир»,1969; 437 стр.
Ляликов А.П. Методы инженерного творчества. Учебно-методическое пособие. Изд.СПбГМТУ, 1996; 212 с.
Альтшуллер Г.С. Творчество как точная наука (теория решения изобретателских задач). М., «Советское радио», 1979, 184 стр.
Адлер Ю.П., Маркова Е.В., Грановский Ю.В. Планирование экспериментов при поиске оптимальных условий. М., «Наука»,1976; 245 стр.
Михайлов В.И., Федосов К.М. Планирование экспериментов в судостроении. Л., «Судостроение»,1978; 159 стр.
Сокальский Г.А. Определение закона распределения случайной величины. Методические разработки. Л., ИПК Судпрома, 1985; 36 стр.
Равин А.А. Экспериментальные методы решения оптимизационных задач. Л., ИПК Судпрома, 1985; 28 стр.
Равин А.А. Лабораторные работы по оптимизации проектных решений с помощью персональной ЭВМ. - Л., ИПК, 1987.
Геминтерн В.И., Каган Б.М. Методы оптимального проектирования. - М., «Энергия», 1980
Хоменюк В.В. Элементы теории многоцелевой оптимизации. – М., Наука, 1988
Равин А.А. Методические указания «Методология инженерной деятельности. Лабораторные работы и тест». СПбГМТУ, 2003, 35 стр.
Приложение 1
Значения G - критерия (Кохрена) для уровня значимости α = 1- Р = 0,05
f2=N |
Степень свободы f1=n -1 |
|||||||||||
l |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
16 |
∞ |
|
2 3 4 5 6 7 8 9 10 12 16 20 24 30 40 ∞ |
0,9985 0,9669 0,9065 0,8412 0,7808 0,7271 0,6798 0,6385 0,6020 0,5410 0,4546 0,3894 0,3434 0,2929 0,2370 0 |
0,9750 0,8709 0,7679 0,6838 0,6161 0,5612 0,5157 0,4775 0,4450 0,3924 0,3218 0,2705 0,2354 0,1980 0,1576 0 |
0,9392 0,7977 0,6841 0,5981 0,5321 0,4800 0,4377 0,4027 0,3733 0,3264 0,2648 0,2205 0,1907 0,1593 0,1259 0 |
0,9057 0,7457 0,6287 0,5441 0,4803 0,4307 0,3919 0,3584 0,3311 0,2880 0,2319 0,1921 0,1656 0,1377 0,1082 0 |
0,8772 0,7071 0,5895 0,5065 0,4447 0,3974 0,3595 0,3286 0,3029 0,2624 0,2103 0,1735 0,1493 0,1237 0,0968 0 |
0,8534 0,6771 0,5598 0,4783 0,4184 0,3726 0,3362 0,3067 0,2823 0,2439 0,1948 0,1602 0,1374 0,1137 0,0887 0 |
0,8332 0,6530 0,5365 0,4564 0,3980 0,3535 0,3185 0,2901 0,2666 0,2299 0,1829 0,1501 0,1286 0,1061 0,0827 0 |
0,8159 0,6333 0,5175 0,4387 0,3817 0,3384 0,3043 0,2768 0,2541 0,2187 0,1736 0,1422 0,1216 0,1002 0,0780 0 |
0,8010 0,6167 0,5017 0,4241 0,3682 0,3259 0,2926 0,2659 0,2439 0,2098 0,1660 0,1357 0,1160 0,0958 0,0745 0 |
0,7880 0,6025 0,4884 0,4118 0,3568 0,3154 0,2829 0,2568 0,2353 0,2020 0,1597 0,1303 0,1113 0,0921 0,0713 0 |
0,734! 0,5466 0,4366 0,3645 0,3135 0,2756 0,2462 0,2226 0,2032 0,1737 0,1365 0,1108 0,0942 0,0771 0,0595 0 |
0,5000 0,3333 0,2500 0,2000 0,1667 0,1429 0,1250 0,1111 0,1000 0,0833 0,0634 0,0500 0,0417 0,0333 0,0250 0 |
Приложение 2
Значения t - критерия (Стьюдента) для доверительной
вероятности Р и числа степеней свободы f = n-1
f |
Доверительная вероятность Р |
||||
0,90 |
0,95 |
0,975 |
0,99 |
0.995 |
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 25 30 40 60 120 ∞ |
6,314 2,920 2,353 2,132 2,015 1,943 1,895 1,860 1,833 1,812 1,796 1.782 1,771 1,761 1,753 1,746 1,740 1,743 1,729 1,725 1,708 1,697 1,684 1,671 1,658 1,645 |
12,706 4,303 3,182 2,776 2,571 2,447 2,365 2,306 2,262 2,228 2,201 2,179 2,160 2,145 2,131 2,120 2,110 2,101 2,093 2,086 2,059 2,042 2,021 2,000 1,980 1,960 |
25,452 6,205 4,176 3,495 3,163 2,969 2,841 2,751 2,685 2,634 2,593 2,560 2,533 2,510 2,490 2,473 2,458 2,445 2,433 2,423 2,385 2,360 2,329 2,299 2,270 2,241 |
63,657 9,925 5,841 4,604 4,032 3,707 3,499 3,355 3,250 3,169 3,106 3,054 3,012 2,977 2,947 2,921 2,898 2,878 2,861 2,845 2,787 2,750 2,705 2,660 2,617 2,576 |
127,32 14,089 7,453 5,598 4,773 4,317 4.029 3,832 3,690 3,581 3,497 3,428 3,372 3,326 3,286 3,252 3,222 3,197 3,174 3,153 3,078 3,030 2,971 2,915 2,860 2,807 |
Дальше: Приложение 3 Значения F-критерия (Фишера) при доверительной вероятности Р=0,95
f2= N(n -1) |
Степень свободы f1 – число незначимых коэффициентов регрессии |
|||||||||||
I |
2 |
3 |
4 |
5 |
6 |
8 |
12 |
16 |
21 |
50 |
∞ |
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 20 30 40 50 100 ∞ |
161,4 19,51 10,13 7,71 6,61 5,99 5,59 5,32 5,12 4,96 4,84 4,75 4,67 4,60 4,54 4,49 4,35 4,17 4,08 4,03 3,94 3,84 |
199,5 19,0 9,55 6,94 5,79 5,14 4,74 4,46 4,26 4,10 3,98 3,88 3,80 3,74 3,68 3,63 3,49 3,32 3,23 3,18 3,09 2,99 |
215,7 19,6 9,28 6,59 5,41 4,76 4,35 4,07 3,86 3,71 3,59 3,49 3,41 3,34 3,29 3,24 3,10 2,92 2,84 2,79 2,60 2,60 |
224,6 19,24 9,12 6,39 5,19 4,53 4,12 3,84 3,63 3,48 3,36 3,26 3,18 3,11 3,06 3,01 2,87 2,69 2,61 2,56 2,46 2,37 |
230,2 19,30 9,01 6,26 5,05 4,39 3,97 3,69 3,48 3,33 3,20 3,11 3,02 2,96 2,90 2,85 2,71 2,53 2,45 2,40 2,30 2,21 |
230,2 19,33 8,94 6,16 4,95 4,28 3,87 3,58 3,37 3,22 3,09 3,00 2,92 2,85 2,79 2,74 2,60 2,42 2,34 2.29 2,19 2,09 |
234,0 19,37 8,84 6,07 4,82 4,15 3,73 3,44 3,23 3,07 2,95 2,85 2,77 2,70 2,64 2,59 2,45 2,27 2,18 2,13 2,03 1,94 |
243,9 19,41 8.74 5,91 4,68 4,00 3,57 3,28 3,07 2,91 2,79 2,69 2,60 2,53 2,48 2,42 2,28 2,09 2,00 1,95 1,85 1,75 |
246,5 19,45 8,69 5,84 4,60 3,92 3,49 3,20 2,98 2,82 2,70 2,60 2,51 2,44 2,39 2,33 2,18 1,99 1,90 1,85 1,75 1,64 |
249,0 19,47 8,64 5.77 4,53 3,84 3,41 3,12 2,90 2,74 2,61 2.50 2,42 2,35 2,29 2,24 2,03 1,89 1,79 1,74 1,63 1,52 |
251,8 19,47 8,58 5,70 4,44 3,75 3,32 3,03 2,80 2,64 2,50 2,40 2,32 2,24 2,18 2,13 1,96 1.76 1,66 1,60 1,48 1,35 |
254,3 19,50 8,53 5,63 4,36 3,67 3,28 2,93 2,71 2,54 2,40 2.30 2,21 2,13 2,07 2,01 1,84 1,62 1,51 1,44 1,28 1,00 |