Испытания на основе математического моделирования объектов испытаний
Методы испытаний ОИ на основе их математических моделей по решаемым задачам можно разделить на две группы.
Первую группу составляют исследовательские доводочные испытания, когда необходимо оценить достигнутые при разработке показатели надежности, прочности или устойчивости ОИ к внешним воздействиям, а также предварительно оценить техническую совместимость узлов (блоков), в один из которых были внесены определенные конструктивные или технологические изменения.
В эту группу входят расчетные (функциональные и статистические) методы: максимума-минимума, квадратичного сложения, вероятностный и статистических испытаний. Они позволяют аналитическим путем получить количественную оценку изменения выходного параметра ОИ при наличии его математической модели и данных по изменению ее параметров с течением времени (в том числе статистических данных по разбросу параметров) или вследствие воздействия внешних факторов (климатических, механических и др.) и входных параметров ОИ. Результаты расчетов позволяют определить оценку надежности, прочности и устойчивости ОИ к внешним воздействиям.
Вторую группу составляют четыре вида методов, позволяющих имитировать изменение параметров математической модели по заданному закону:
1. Метод физического аналога состоит в том, что в качестве модели ОИ используется физическое устройство, аналогичное ОИ по составу и структуре, искусственное изменение параметров которого отражает их изменения в реальных условиях. Такой физический аналог может строиться с использованием реальных блоков или в виде специальных вычислительных алгоритмов. Если модель ОИ представлена физическим устройством, то испытания сводятся к изменению параметров его блоков и контролю за соответствующим изменением выходных параметров модели ОИ. Если модель ОИ представлена в виде специального вычислительного алгоритма, то изменяются уравнения, описывающие работу соответствующих блоков, и фиксируется изменение выходных параметров модели ОИ. Рассматриваемый вид испытаний используется главным образом на этапе разработки ОИ в качестве исследовательских, а также для отработки действий в нештатных ситуациях или оценки схемной защищенности от экстремальных нагрузок.
2. Метод граничных испытаний преследует цель определения зависимостей между предельно допустимыми значениями параметров ОИ и режимом эксплуатации и имеет те же области применения, что и первый метод. Если модель ОИ представлена физическим устройством, то в процессе моделирования его параметры изменяются так, чтобы отразить воздействия, соответствующие экстремальным режимам эксплуатации. При этих режимах проводятся измерения значений параметров устройства. Результаты измерений используются как контрольные точки, по которым отлаживается математическое моделирование зависимости параметров модели от режима эксплуатации. Если модель ОИ реализуется в виде программы для ЭВМ, то указанные воздействия отражаются путем изменения исходных данных или корректировки программы. Отличие второго метода от первого состоит в том, что на его основе оценивается только предельная нагрузка, соответствующая границе работоспособности ОИ.
3. Метод матричных испытаний является дальнейшим развитием метода граничных испытаний. Сущность его заключается в математическом моделировании рабочей области объекта испытаний и сопоставлении ее с областью работоспособности ОИ. Это моделирование осуществляется следующим образом. Диапазон возможных значений каждого входного параметра 0у min—0у max разбивается на интервалы (кванты), «представителями» которых являются средние значения параметров в середине кванта. Выбираются определенные сочетания значений квантов, их представители вводятся в модель, расчет по которой дает значение выходного параметра. Перебором сочетаний значений квантов различных входных параметров обеспечивают реализацию процесса массовых испытаний.
Если реализация оказывается работоспособной или неработоспособной при значении параметра qi, то эту реализацию относят ко всем значениям параметра, лежащим в данном кванте.
Все реализации записывают в виде матрицы, которая учитывает возможные несовместные комбинации представителей квантов. Общее число реализаций будет определяться числом возможных сочетаний величин матрицы
011 |
012 — |
01/ — |
0і/ |
021 |
022 — |
02/ — |
02/ |
• • • 0«1 |
0/»2 •” |
• • • Япі ••• |
• t « *nf |
где qn,…,qjf — кванты допуска на/-параметры.
Упорядоченную матричным построением последовательность возможных ситуаций обозначим аі, а2,…,(Ху,…,а^. Перебор ситуаций {oty} в модели позволяет учесть влияние всех факторов и определить некоторую упорядоченную последовательность Q отказов, обозначаемую Pi, Р2 > * • •»Рц > * * * > Рр •
При обработке результатов достаточно знать положение ситуации в матрице, чтобы определить, какими квантами каждого параметра образуется данная ситуация.
В процессе матричных испытаний (они обычно выполняются с использованием специальных автоматов) исследуется переходная характеристика каждой реализации модели ОИ и в соответствии с выбранным критерием отказа делается заключение о работоспособности или неработоспособности данной реализации. По результатам испытаний строится л-мерная область работоспособности. По изменению формы и положения этой области можно давать качественную оценку изменения надежности, прочности или устойчивости к внешним воздействиям модели ОИ. Количественная оценка определяется как вероятность нахождения выходного параметра модели ОИ в границах области работоспособности.
Метод матричных испытаний используется для предварительной оценки работоспособности с учетом законов распределения производственных погрешностей значений параметров элементов ОИ и при учете влияния дестабилизирующих факторов (температуры, влажности, давления, радиации и др.), значений погрешностей параметров элементов и процессов их старения; показателей надежности, прочности и устойчивости ОИ, влияния технологических допусков (погрешностей) на работоспособность ОИ. Кроме того, по результатам матричных испытаний могут определяться допустимые значения погрешности измерения параметров элементов, узлов и блоков, при которых не нарушается их техническая совместимость.
4. Метод статистических испытаний (метод Монте-Карло) состоит в построении искусственного случайного процесса, реализуемого с использованием ЭВМ. Искусственный случайный процесс в виде последовательности случайных чисел имитирует изменение соответствующих внутренних параметров математической модели ОИ или внешнего воздействия, что позволяет оценить статистические характеристики процесса изменения выходных параметров ОИ, т. е. оценить надежность, прочность и устойчивость к внешним воздействиям.
Способы реализации метода статистических испытаний весьма разнообразны, поэтому рассмотрим его общий принцип. Датчик случайных чисел генерирует случайные числа, равномерно распределенные в йнтервале 0,…,1. Эти случайные числа поступают в преобразователь, на выходе которого получаются случайные числа, принадлежащие к генеральной совокупности {Л}, имеющей заданную функцию плотности распределения /(л). Преобразование выполняется путем решения относительно величины Л уравнения
(4.1)
Модель ОИ представляется в виде алгоритма, описываемого в общем случае зависимостями вида
Nj =/(?1……….
j = 1 ,k; і = 1,т; Aj < Nj < Bj,
где N — значения выходных параметров ОИ; qt — значения параметров входящих в ОИ элементов; Aj, Bj — предельно допустимые по ТУ значения выходных параметров.
Для нахождения выходных параметров в уравнение (4.2) подставляются значения Яіся параметров элементов > опреде
ленные случайным образом для любого значения нагрузки (воздействия) при оценке прочности и устойчивости или времени при оценке надежности:
ЯіаііЛ = Яю + Д*(ф) + /[МфН], (4.3)
где (qio+Aqj) — математическое ожидание параметра qі при нагрузке (времени) Ф; /Га9,(Ф)т|, ] — случайная величина, зависящая от вида закона распределения q / , среднеквадратического отклонения <^/(Ф) и случайного числа Т|/, имеющего плотность распределения /(rj). Для нормального закона распределения Г|/ выражение (4.3) приобретает вид
Яі сл (ф) = тд, (Ф) + (Ф), (4.4)
где ntqi(Ф) — математическое ожидание qi в зависимости от Ф.
Сравнивая значения N-, полученные при подстановке (4.3) и (4.4) в (4.2), с предельно допустимыми значениями А-, В= по числу отказов (здесь Ф = Т — время наблюдения), получим одну из характеристик надежности — вероятность безотказной работы:
Р{Т) = 1-1/L,
где / — число реализаций, при которых наблюдается выход хотя бы одного из параметров N, за нормы ТУ; L — заданное число реализаций, выбираемое исходя из требуемой точности расчета.
Аналогично решением уравнения (4.5) относительно Ф™ могут быть определены значения предельной нагрузки Ф_, которые соответствуют отказу с заданной вероятностью 1 — ^(Фцр) •
Подробное изложение метода статистических испытаний приведено в [99].