Пример эксперимента по дробному плану
12.3. Пример эксперимента по дробному плану
В последовательном экспериментировании, если для достижения желаемого уровня точности не нужно использовать всё число опытов полного или повторяемого дробного факторного плана, то лучше всего начать с эксперимента по дробному факторному плану. Такие планы всегда могут быть позднее увеличены, если в этом появится необходимость. В этом случае эксперимент по первому дробному плану покажет, как нужно выбрать следующий дробный план для соединения в больший план таким образом, чтобы ответить на вопросы и разрешить неопределённости.
Рассмотрим эксперимента по плану 2III7–4 [Боков (1988), Bokov (2002)]. Этот эксперимент поставлен с целью изучения времени срабатывания пневматического преобразователя линейных перемещений, схема которого показана на Рис.12.3.1. В преобразователе сжатый воздух через регулятор давления поступает в два канала. В нижнем канале воздух через входной дроссель подаётся в камеру противодавления и истекает из неё через регулируемый дроссель в атмосферу. В верхнем канале сжатый воздух поступает в рабочую камеру через переменный входной дроссель с коническим клапаном и вытекает из неё через дроссель сопло-заслонка. Разделяет камеры очень эластичная резиновая мембрана с жёстким центром. Жёсткий центр мембраны через шток связан с индикатором перемещений и на штоке установлен конический клапан, служащий для автоматического изменения проходного сечения входного дросселя рабочей камеры.
Рис.12.3.1. Схема пневматического преобразователя линейных перемещений
Для установки преобразователя в рабочее положение между соплом и заслонкой устанавливается рабочий зазор, создающий в рабочей камере измерительное давление. Выходной дроссель камеры противодавления при этом регулируется так, чтобы давления в камерах были приблизительно равны и конический клапан находился в рабочем положении. В этом положении создаётся проходное сечение переменного дросселя равное половине наибольшего значения, и индикатор показывает нуль. Перемещение заслонки ведёт к изменению зазора между ней и соплом и давления в рабочей камере, что приводит к изменению положения конического клапана и показания индикатора. С целью устранения сухого трения шток устанавливается в аэростатических направляющих.
Когда заслонка отсутствует, то сжатый воздух свободно истекает из сопла и конический клапан находится в верхнем положении. При этом расход воздуха через преобразователь максимальный. При появлении заслонки перед соплом давление в рабочей камере увеличивается, что вызывает перемещение вниз жёсткого центра мембраны и штока с коническим клапаном. Это перемещение будет до тех пор, пока не установится новое равновесие сил между разностью давлений в камерах, весом подвижных деталей и усилием пружины сжатия. Отрезок времени между моментом появления перед соплом заслонки и моментом установки подвижных деталей в новое положение является временем срабатывания преобразователя. Так как в преобразователе используется очень эластичная мембрана и давление в камере противодавления почти не изменяется, то масса воздуха в рабочей камере сильно не меняется и поэтому время его срабатывания короткое.
Исследование было начато постановкой эксперимента из восьми опытов по дробному факторному плану для изучения влияния семи факторов на время срабатывания преобразователя. Опыты проводились в случайной последовательности. Используемые в эксперименте семь факторов показаны в таблице 12.3.1. Эти семь факторов определяют различные состояния элементов преобразователя и условия его работы.
Для создания дробного плана эксперимента использовался полный план 23. На его основе с использованием произведений Адамара были получены вектор-столбцы уровней взаимодействий исходных факторов, как показано в таблице 12.2.1. Все воздействия взаимодействий двух и трёх факторов были совмещены с воздействиями новых факторов х4, х5, х6 и х7 из таблицы 12.3.1. Полученный в результате дробный план эксперимента со столбцами уровней всех факторов показан в таблице 12.3.2. Этот факторный план 27–4 является одной шестнадцатой частью полного факторного плана 27 и обладает проективностью Р=2.
Рекомендуемые материалы
Таблица 12.3.1. Факторы эксперимента и обозначения их воздействий на отклик
Факторы | Нормир. факторы | Воздействия факторов |
Эффективная площадь мембраны, ξ1 | х1 | β1 |
Площадь проходного сечение дросселя сопло-заслонка, ξ2 | х2 | β2 |
Площадь проходного сечения входного дросселя, ξ3 | х3 | β3 |
Объём рабочей камеры, ξ4 | х4 | β4 |
Жёсткость пружины, ξ5 | х5 | β5 |
Объём камеры противодавления, ξ6 | х6 | β6 |
Давление сжатого воздуха на входе в преобразователь, ξ7 | х7 | β7 |
Создающими план 27–4 на основе плана 23 генерирующими равенствами являются х4=х1◦х2, х5=х1◦х3, х6=х2◦х3, х7=х1◦х2◦х3 или равносильно их можно записать в виде определяющего отношения 1=х1◦х2◦х4=х1◦х3◦х5=х2◦х3◦х6=х1◦х2◦х3◦х7, где х1, х2, ..., х7 – вектор-столбцы уровней факторов х1, х2, ..., х7 в таблице 12.3.2. Однако это определяющее отношение является неполным, так как на его основе невозможно найти все структуры совмещённых воздействий. Для нахождения полного определяющего отношения надо ещё добавить все результаты произведений Адамара, получающиеся в результате всех возможных перемножений четырёх членов исходного определяющего отношения. Например, умножение первых двух членов даёт х1◦х2◦х4◦х1◦х3◦х5 =х1◦х1◦х2◦х3◦х4◦х5 =х2◦х3◦х4◦х5, так как х1◦х1=1. В результате, полным определяющим отношением является 1=х1◦х2◦х4=х1◦х3◦х5=х2◦х3◦х6 =х1◦х2◦х3◦х7 =х2◦х3◦х4◦х5 =х1◦х3◦х4◦х6 =х3◦х4◦х7 =х1◦х2◦х5◦х6 =х2◦х5◦х7 =х1◦х6◦х7=х4◦х5◦х6 =х1◦х4◦х5◦х7 =х2◦х4◦х6◦х7 =х3◦х5◦х6◦х7 =х1◦х2◦х3◦х4◦х5◦х6◦х7.
Таблица 12.3.2. План из восьми опытов для изучения влияний семи факторов
Опыты | х1 | х2 | х3 | х4 | х5 | х6 | х7 | у1 | у2 | у3 | |
1 | –1 | –1 | –1 | +1 | +1 | +1 | –1 | 144 | 136 | 148 | 143 |
2 | +1 | –1 | –1 | –1 | –1 | +1 | +1 | 88 | 84 | 91 | 88 |
3 | –1 | +1 | –1 | –1 | +1 | –1 | +1 | 140 | 148 | 152 | 147 |
4 | +1 | +1 | –1 | +1 | –1 | –1 | –1 | 136 | 137 | 137 | 137 |
5 | –1 | –1 | +1 | +1 | –1 | –1 | +1 | 99 | 98 | 96 | 98 |
6 | +1 | –1 | +1 | –1 | +1 | –1 | –1 | 89 | 89 | 86 | 88 |
7 | –1 | +1 | +1 | –1 | –1 | +1 | –1 | 142 | 135 | 144 | 140 |
8 | +1 | +1 | +1 | +1 | +1 | +1 | +1 | 110 | 105 | 112 | 109 |
В этом определяющем отношении есть члены, имеющие наименьшее число сомножителей равное трём. Поэтому данный план имеет разрешающую способность III и записывается в виде 2III7–4. Произведения Адамара членов определяющего отношения на х1 преобразуют его к виду х1=х2◦х4=х3◦х5=х1◦х2◦х3◦х6 =х2◦х3◦х7 =х1◦х2◦х3◦х4◦х5 =х3◦х4◦х6 =х1◦х3◦х4◦х7 =х2◦х5◦х6 =х1◦х2◦х5◦х7 =х6◦х7=х1◦х4◦х5◦х6 =х4◦х5◦х7 =х1◦х2◦х4◦х6◦х7=х1◦х3◦х5◦х6◦х7 =х1◦х2◦х3◦х4◦х5◦х6◦х7. В полученном, если считать, что все взаимодействия трёх и большего числа факторов дают пренебрежимо малые воздействия, то соответствующие им члены определяющего отношения можно сократить и на основе оставшихся х1=х2◦х4=х3◦х5 =х6◦х7 получается сокращённая структура β1+β24+β35+β67 совместных воздействий фактора х1 и двухфакторных взаимодействий остальных факторов. Сокращённые структуры совместных воздействий всех остальных факторов с двухфакторными взаимодействиями приведены в таблице 12.3.3.
Таблица 12.3.3. Совместные воздействия факторов и двухфакторных взаимодействий и результаты их оценки для данных эксперимента по плану 2III7–4
Факторы | Совместные воздействия | Результаты оценки |
1. Эффективная площадь мембраны, х1 | β1+β24+β35+β67 | b1=–13,25 |
2. Проходное сечение дросселя сопло-заслонка, х2 | β2+β14+β36+β57 | b2=14,583 |
3. Проходное сечение входного дросселя, х3 | β3+β15+β26+β47 | b3=–9,833 |
4. Объём рабочей камеры, х4 | β4+β12+β37+β56 | b4=2,917 |
5. Жёсткость пружины, х5 | β5+β13+β27+β46 | b5=3,0 |
6. Объём камеры противодавления, х6 | β6+β17+β23+β45 | b6=1,333 |
7. Давление на входе в преобразователь, х7 | β7+β16+β25+β34 | b7=–8,333 |
Среднее переменных отклика в опытах эксперимента | b0=118,583 |
Эти структуры совместных воздействий можно найти и с использованием матрицы совмещения для воздействий факторов эксперимента по плану 2III7–4 с воздействиями двухфакторных взаимодействий. Функция модели этого эксперимента с учётом всех двухфакторных взаимодействий записывается в виде Е(у)=Х1β1+Х2β2, где вектор параметров β1Т=[β0, β1, β2, β3, β4, β5, β6, β7], вектор параметров β2Т=[β12, β13, β14, β15, β16, β17, β23, β24, β25, β26, β27, β34, β35, β36, β37, β45, β46, β47, β56, β57, β67], Х1 – матрица модели эксперимента без учёта двухфакторных взаимодействий и Х2 – матрица уровней двухфакторных взаимодействий факторов из таблицы 12.3.1. Как и в предыдущем разделе, матрица совмещения воздействий находится по формуле А=(X1TX1)–1X1TХ2. Результат её вычисления представлен в таблице 12.3.4.
Таблица 12.3.4. Матрица коэффициентов совмещения воздействий по плану 2III7–4
β12 | β13 | β14 | β15 | β16 | β17 | β23 | β24 | β25 | β26 | β27 | β34 | β35 | β36 | β37 | β45 | β46 | β47 | β56 | β57 | β67 | |
β0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
β1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | +1 | 0 | 0 | 0 | 0 | +1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | +1 |
β2 | 0 | 0 | +1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | +1 | 0 | 0 | 0 | 0 | 0 | +1 | 0 |
β3 | 0 | 0 | 0 | +1 | 0 | 0 | 0 | 0 | 0 | +1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | +1 | 0 | 0 | 0 |
β4 | +1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | +1 | 0 | 0 | 0 | +1 | 0 | 0 |
β5 | 0 | +1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | +1 | 0 | 0 | 0 | 0 | 0 | +1 | 0 | 0 | 0 | 0 |
β6 | 0 | 0 | 0 | 0 | 0 | +1 | +1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | +1 | 0 | 0 | 0 | 0 | 0 |
β7 | 0 | 0 | 0 | 0 | +1 | 0 | 0 | 0 | +1 | 0 | 0 | +1 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
Вместе с этой лекцией читают "14.3 Русско-японская война 1904-1905 гг".
Все опыты эксперимента повторялись трижды и их результаты в миллисекундах представлены в таблице 12.3.2 в столбцах у1, у2 и у3. На основе этих данных можно оценить дисперсию ошибок опытов следующим образом. По формуле (6.4.1) находится сумма квадратов чистых ошибок, а по формуле (6.4.2) - число её степеней свободы. В разделе 11.3 показано, что среднеквадратичное суммы квадратов чистых ошибок является несмещённым результатом оценки дисперсии ошибок. Отсюда получаем результат s2=16 оценки дисперсии. Результат оценки дисперсии оценочных воздействий получается 0,125s2=2, а равная корню квадратному из результата оценки дисперсии их стандартная ошибка равна ±1,414.
Оценка совмещённых воздействий делается методом наименьших квадратов с использованием постулируемой модели у=β0+β1х1+β2х2+β3х3+β4х4+β5х5+β6х6+β7х7+ε и вектор-столбца усреднённых значений переменных отклика из таблицы 12.3.2. Результаты оценки представлены в таблице 12.3.3. Сравнение результатов оценки со стандартной ошибкой показывает, что выделенные жирным шрифтом значения b1, b2, b3, b4, b5 и b7 явно выделяются из шума, так как отличаются от стандартной ошибки больше чем в 2 раза.
Простейшее объяснение результатов оценки совместных воздействий состоит в том, что из семи факторов только четыре: х1 - эффективная площадь мембраны, х2 - проходное сечение дросселя сопло-заслонка, х3 - проходное сечение входного дросселя камеры противодавления и х7 - давление сжатого воздуха на входе в преобразователь, оказывают существенное влияние и их воздействия совмещены с воздействиями двухфакторных взаимодействий. Однако, при четырёх активных факторах, показанный в таблице 12.2.1 план 24–1 состоит также из восьми опытов. Поэтому, если считать, что эксперимент проводился по плану 24–1 из таблицы 12.2.1 с четырьмя факторами, то получается, что b1, b2, b3 и b7 являются результатами оценки для структур совмещения воздействий самих этих факторов с воздействиями трёхфакторных взаимодействий.
Эксперимент по плану 24–1 представлен графически на Рис.12.3.2. Из него видно, что изменение факторов х1, х3 и х7 с нижнего на верхний уровень ведёт к уменьшению времени срабатывания преобразователя, а изменение фактора х2 с нижнего на верхний уровень ведёт к увеличению его времени срабатывания. Это следует также из анализа результатов оценки совмещённых воздействий в таблице 12.3.3. Результат b7 оценки совместных воздействий фактора х7 и двухфакторных взаимодействий выделяется из шума и поэтому фактором х7 нельзя пренебрегать.
Рис.12.3.2. Графическое изображение эксперимента по плану 24–1