УДК 681.2:53.088
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРИ АНАЛИЗЕ МЕТРОЛОГИЧЕСКОЙ НАДЕЖНОСТИ АНАЛОГОВЫХ БЛОКОВ ИНФОРМАЦИОННО-ИЗМЕРИТЕЛЬНЫХ СИСТЕМ
Т. И. Чернышова, В. В. Третьяков
Кафедра «Конструирование радиоэлектронных и микропроцессорных систем», ФГБОУ ВПО «ТГТУ»; crems@crems.jesby.tstu.ru
Ключевые слова и фразы: аналоговый блок; внешний влияющий фактор; информационно-измерительная система; математическая модель; метрологическая надежность; метрологическая характеристика; метрологический ресурс.
Аннотация: Разработана математическая модель изменения во времени метрологической характеристики аналоговых блоков информационно-измерительных систем с учетом внешних влияющих факторов, а также обобщенная математическая модель параметров элементной базы, используемая в процедуре статистического моделирования метрологических характеристик исследуемых измерительных средств, позволяющая анализировать уровень их метрологической надежности на разных этапах эксплуатации. В качестве основного показателя метрологической надежности аналоговых блоков информационно-измерительных систем рассмотрен метрологический ресурс. Приведены результаты исследования метрологической надежности типового блока, входящего в состав измерительного канала информационно-измерительных систем, получены значения метрологического ресурса при различных значениях дестабилизирующих факторов и в различных временных сечениях. Основным преимуществом предлагаемого подхода к математическому моделированию метрологических характеристик является возможность определения уровня метрологической надежности измерительного средства на любом этапе эксплуатации с учетом совместного действия изменяющихся внешних факторов.
Информационно-измерительные системы (ИИС) используются во всех отраслях промышленности, машиностроения, науки и производства. Сложность и ответственность выполняемых ими функций обуславливает необходимость обеспечения их высокого уровня метрологической надежности (МН). Как известно, метрологическая надежность определяется характером и темпом изменения нормируемых метрологических характеристик (МХ) аналоговых блоков (АБ) ИИС [1]. Основным показателем МН является метрологический ресурс. Оценка метрологического ресурса на этапе проектирования ИИС позволяет потребителю делать выводы о состоянии МН ИИС в произвольные моменты времени предстоящей эксплуатации, а также определить необходимые сроки поверок и технического обслуживания ИИС.
Решение задачи оценки МН при проектировании АБ ИИС формируется с применением метода аналитико-вероятностного прогнозирования, основанного на анализе нестационарных случайных процессов изменения во времени метрологических характеристик АБ ИИС с использованием их математических моделей (ММ) [1].
Такие ММ строятся на основе статистического моделирования значений МХ по данным об изменении параметров элементной базы в процессе предстоящей эксплуатации. Как показали исследования, внешние факторы, влияющие на процесс изменения во времени МХ проектируемых АБ, существенно ускоряют процесс старения параметров АБ и способствуют уменьшению МН ИИС.
Основными дестабилизирующими факторами, оказывающими влияние на изменение метрологической характеристики АБ ИИС, являются температура, влажность и давление окружающей среды, а также радиационный фон.
Для конкретного АБ имеется ряд ММ изменения во времени параметров элементов:
41 (<м ) = /1 (<М, ^01);
42 (Ф1 ) = /2 (^ 402 );
Ч (ф ) = / (Ф л 40/);
4п (Ф^) = /п (,t, 40п ).
где ф = {,ф2,...,фI} - вектор влияющих факторов; I - число влияющих факторов; п - число элементов АБ, входящих в математическую модель метрологической характеристики; t - время эксплуатации; - реальное значение параметра /-го элемента, / = 1, ..., п; ^ - номинальное значение параметра /-го элемента.
С учетом указанных внешних влияющих факторов математическая модель МХ АБ, построенная с применением анализа структурной и принципиальной схем блока, запишется в виде
£ = Е (х, §, ф) = Е [ х, §, (Т, Е, Р, Е, t ), (2)
где х - входной сигнал; § = {§1, §2,-., - вектор параметров элементов; Т - температура, °С; Е - влажность, %; Р - давление, атм; Е - радиационный фон, мкР/ч.
Построенная таким образом модель вида (2) используется далее в процедуре статистического моделирования [1 - 3]. Алгоритм моделирования состоит из последовательного расчета характеристик закона распределения значений параметров комплектующих элементов исследуемых блоков и моделирования реализаций МХ в различных временных сечениях области контроля t/, / = 1, ..., К, где К - число временных сечений.
Результатом статистического моделирования является совокупность значений математического ожидания МХ АБ в различные моменты времени эксплуатации т,,^ 1), ..., т,(К) и значений среднеквадратического отклонения ..., о,(К) при вариации внешних факторов.
На основе результатов статистического моделирования строится математическая модель изменения во времени метрологической характеристики [1 - 3], представляющая собой совокупность аналитических зависимостей, полученных для функций изменения во времени математического ожидания т((), и функций, характеризующих изменение во времени границ отклонений возможных значений МХ от ее математического ожидания:
т,(ф, t);
(3)
V± о(ф, t) = т, (ф, t) + с о, (ф, t),
где с - коэффициент, выбираемый в зависимости от уровня доверительной вероятности р, с = 3 при р = 0,9973.
Построенная ММ (3) является основой для оценки метрологического ресурса путем экстраполяции зависимостей (3) на область будущих значений времени эксплуатации. Очевидно, что точность полученных результатов прогнозирования будет зависеть от адекватности применяемых ММ (1) для элементной базы АБ.
Определение математического описания процесса старения элементной базы АБ ИИС под действием внешних факторов является важным условием для оценки влияния внешних воздействий на процесс изменения во времени метрологической характеристики АБ ИИС. Применение этого описания в априорной модели процесса старения элементной базы АБ ИИС позволяет корректировать априорные данные, и тем самым оценку метрологической надежности с учетом реальных условий эксплуатации. В качестве внешних дестабилизирующих факторов, как было указано выше, принимаются температура Т, влажность F, давление P и радиационный фон E окружающий среды. Проведенные исследования показали, что температура Т и влажность F оказывают влияние на все элементы АБ ИИС, тогда как давление P существенно влияет только на объемно пористые конденсаторы, которые могут применяться далеко не во всех ИИС. Радиационный фон E оказывает влияние на элементы микроэлектроники (микросхемы, большие и сверхбольшие интегральные схемы), а также на активные элементы (диоды, транзисторы и др.). Таким образом, при построении ММ МХ необходимо учитывать типы элементов, применяемых в АБ ИИС, так как при отсутствии тех или иных типов элементов нет необходимости вводить все упомянутые факторы в модель (1).
Известно, что деградация элементной базы происходит как на параметрическом, так и на структурном уровнях [2]. На параметрическом уровне изменения происходят из-за негативного влияния окружающей среды на параметры элементов, такие как проводимость, емкость, диэлектрическая постоянная и др.
Скорость деградации элементов в нормальных условиях Vo определяется максимальным изменением параметра е^ за время т (срок сохраняемости)
(4)
В условиях влияния повышенной влажности и температуры окружающей среды, а также в случае применения в АБ ИИС средств микроэлектроники, активных элементов и объемно пористых конденсаторов, учитывается влияние повышенного давления и радиационного фона, и старение ускоряется.
Количественно старение характеризуется коэффициентом ускорения и описывается уравнением [1, 2]
( v\n (
Af = — =
F_
vF0,
:exp
exp
Ea
V k V
1 _ J_
T T0
W
ln (P _ Po)
P _ P0
(5)
4M1M 2 (1 + M 2)
2 NEn _ Eo
где Ар - ускоряющий фактор; V - скорость деградации элементов в реальных условиях; Eа - энергия активации деградационного процесса, эВ; k - постоянная Больцмана; п - степенной показатель влажности; T, ^ - температура в реальных и нормальных условиях эксплуатации соответственно, °С; F, Fo - влажность в реальных и нормальных условиях эксплуатации соответственно, %; P, Po - давление в реальных и нормальных условиях эксплуатации соответственно, атм;
0
I
г - степенной показатель давления; М, Му, М2, Еп, Ео - справочные значения параметров окружающей среды, определяемые на молекулярном уровне.
Согласно выражениям (4) и (5) ММ изменения во времени параметров элементной базы, в которой учитываются воздействия температуры, влажности, давления и радиации, имеет вид
К Е, Т, Р, Е ) (1 + ^с А^ )(1 + а(Т - То ))(1 + р( Е - Е0 ))х х(1 + р(Р - Ро))(1 + а(Е - Ео )) =
^ ¡п (Р - Р0 )г
= ^0
1 + ^ т
/ f е* Г F л
V F0 J
exp
Г e Г
^a
Tf
0 J
J
P - P0
(6)
с exp
4M1M 2
-NE„-E
(М1 + М 2 )2
х(1 + а(Т - То))(1 + р(Е - Ео ))(1 + р(Р - Ро))(1 + а(Е - Ео)),
где а, р, р, с - температурный, влажностный, манометрический и радиационный коэффициенты параметра элемента соответственно.
Математическая модель вида (6) является обобщенным аналитическим описанием процесса изменения параметров элементной базы АБ ИИС при комплексном воздействии внешних факторов. Такая структура ММ применяется в процедуре статистического моделирования МХ АБ ИИС.
Изложенный подход рассмотрен на примере исследования МН преобразователя напряжение-частота (ПНЧ). Данный АБ входит в состав измерительного канала ИИС. Электрическая принципиальная схема этого блока представлена на рис. 1.
Метрологической характеристикой блока является основная относительная погрешность ПНЧ 5. Условие сохранения метрологической исправности записывается в виде: |5(/)| < |5доп|, 5доп = ± 5 %. Математическая модель исследуемой МХ данного преобразователя представлена в виде
5 = ■
Кпр (t,T,F,E)- Кпр.н
кп
пр. н
Кпр (t TT F E ) = -R
1
1
TT 2R6 + U
U"R~+R~+UБЭ c^in r7
Г R v
1 + Ri R3
1 + -
U
VD1
U
ОУ„
(7)
где Кпр(/, Т, Е, Е) - дрейф коэффициента преобразования ПНЧ; Кпр н - номинальный коэффициент преобразования; ип - напряжение нагрузки, В; иБэ - напряжение база-эмиттер, В;
- значение падения напряжения, В; ^оУтах - максимальное выходное напряжение операционного усилителя, В; С1, Л2, •••, Л -емкостные и резистивные параметры элементной базы соответственно.
Рис. 1. Схема электрическая принципиальная ПНЧ
Необходимо отметить, что предложенная обобщенная структура ММ изменения во времени параметров элементной базы АБ с учетом внешних влияющих факторов (6) может быть упрощена, учитывая воздействие температуры Т, влажности F и радиации Е на соответствующие параметры элементов, входящих в (7). Давление Р исключается из общей структуры ММ (6), так как в рассматриваемой схеме ПНЧ отсутствуют элементы, на которые указанный фактор оказывает существенное влияние.
Тогда выражение (6) примет вид
%(t,F, T,P,E) = ^ (1 + «0AFt)(1 + a(T _ T))(l + ß(F _ F0 ))(l + a (E _ E0 )) =
= ^0
1 + -
i( F
F0
nf
exp
Ea ( 1
k V T T0
exp
>
4M1M 2 ( + M 2 )2
-NE„ _ E.
(8)
<(1 + a(T _T0))( + ß(F _F0))(1 + a(E _E0)).
Следующим этапом является статистическое моделирование изменения во времени основной относительной погрешности блока ПНЧ в соответствии с выбранной структурой ММ его элементной базы (8). Процедура статистического моделирования производится с учетом указанных внешних дестабилизирующих воздействий. По результатам статистического моделирования построена многофакторная ММ изменения во времени исследуемой МХ блока, учитывающая условия эксплуатации Т, F, Е:
у-0 (,Т^,Е) = 45,93 -1,77 • 10-3г - 0,773Т - 0,66^ + 1,81Е +
+5,023 • 10-6 Т + 3,967 • 10-6 ^ + 3,123 • 10-6 гЕ + 3,769 • 10-5 Т¥Е -
-3,737 • 10-10 г2 +1,12-10-3Т 2 +1,153 •Ю-4 F2 +1,149 • 10-4 Е2;
ж, (г,Т^,Е) = 45,9303 -1,788 •Ю-3г -0,77Т -0,65F +1,807Е +
+5,0243 • 10_6 tT + 3,9683 • 10_6 tF + 3,1243 • 10_6 tE + 3,7703 • 10_5 TFE _ _3,7363 •10_1012 +1,095 ^10_3T2 +1,1604-10_4 F2 +1,1507-10_4 E2;
V+a (t,T,F,E) = 45,9306 _ 1,791 • 10_31 _ 0,766T _ 0,636F +1,804E + +5,025 • 10_6 tT + 3,969 • 10_6 tF + 3,125 • 10_6 tE + 3,771 • 10_5 TFE _ _3,735 •10_1012 +1,05-10_3T2 +1,17-10_4 F2 +1,1521 -10_4 E2.
(9)
Математическая модель (9) позволяет получить массив значений метрологического ресурса, как основного показателя МН при различных сочетаниях выделенных внешних влияющих факторов. Для отслеживания процесса изменения МХ ИИС может быть построено множество графических интерпретаций для каждого из возможных комбинаций значений внешних факторов. В данной работе проиллюстрированы две графические зависимости: для наиболее благоприятных условий эксплуатации (рис. 2, а) и наиболее экстремальных условий (рис. 2, б). График на рисунке 2, а соответствует следующим условиям окружающей среды: Т = 20 °С; F = 50 % и Е = 50 мкР/ч. Метрологический ресурс при этом составляет гр = 85 500 ч. График на рисунке 2, б соответствует: Т = 70 °С; F = 90 % и Е = 150 мкр/ч. Метрологический ресурс при этом гр = 28 000 ч. На приведенных графиках видно, как существенно изменяется значение метрологического ресурса при различных значениях дестабилизирующих факторов.
+ 5дОП 4 3 2 1 О - 1 - 2 - 3 -4
"бдйп
5, % +5
О-
у/ 1
-у/ .-"-"I 1
1 t 1 2 tP 3 1 4 1 1 5 6 1 7 1 8
ix 10"
2
t х10-
б)
Рис. 2. Изменение во времени МХ ИИС:
Т = 20 °С, Е = 50 %, Е = 50 мкр/ч; б - Т = 70 °С, Е = 90 %, Е = 150 мкр/ч;
7 - Т, Е, Е); 2 - гаД Т, Г, Е); 3 - у-а(/, Т, Е, Е)
ч
3
ч
а
Таким образом, учет в ММ параметров элементной базы при оценке МН проектируемых АБ ИИС в предполагаемых условиях эксплуатации показывает их существенное влияние на метрологический ресурс как основной показатель МН. Анализ полученных результатов по оценке МН при проектировании АБ ИИС позволяет сформировать базу знаний, содержащую априорные значения дестабилизирующих факторов, а также законы изменения МХ в зависимости от изменения внешних влияющих воздействий.
Список литературы
1. Мищенко, С. В. Метрологическая надежность измерительных средств / С. В. Мищенко, Э. И. Цветков, Т. И. Чернышова. - М. : Машиностроение, 2001. -218 с.
2. Чернышова, Т. И. Оценка и повышение метрологической надежности измерительных средств с учетом климатических воздействий / Т. И. Чернышова, М. И. Нистратов // Вестн. Тамб. гос. техн. ун-та. - 2011. - Т. 17, № 1. - С. 24 - 31.
3. Чернышова, Т. И. Оценка метрологического ресурса процессорных средств теплофизических измерений с учетом температурного режима эксплуатации / Т. И. Чернышова, Д. В. Игнатов // Вестн. Тамб. гос. техн. ун-та. - 2005. -Т. 11, № 1Б. - С. 241 - 244.
Mathematical Modeling in the Analysis of Metrological Reliability of Analog Parts of Information-Measuring Systems
T. I. Chernyshova, V. V. Tretyakov
Department "Designing of Electronic and Microprocessor Systems", TSTU; crems@crems.jesby.tstu.ru
Key words and phrases: analog part; outward influencing factor; information-measuring system; mathematical model; metrological reliability; metrological feature; metrological resource.
Abstract: Mathematic model of metrological characteristic changing of analog blocks of information measuring systems was developed with regard to the outward influencing factors. Also, general mathematic model of elements parameters was developed. It can be used in statistic modeling procedure of metrological characteristics to examine measuring devices that allow analyzing metrological reliability level in different operation periods. The main characteristic of metrological reliability of analog blocks of information measuring systems is metrological resource. There are research results of metrological reliability of analog block in the article, that enter into the composition of measuring tract of information measuring systems. Also, there are metrological resource values according to different destabilizing factors values and different time periods. The main advantage of this method of mathematic modeling of metrological characteristics is the possibility to determine the metrological reliability level of measuring device according to any operation period given the combined outward influencing factors.
References
1. Mishchenko S.V., Tsvetkov E.I., Chernyshova T.I.. Metrologicheskaya nadezhnost' izmeritel'nykh sredstv (Metrological reliability of measuring devices), Moscow: Mashinostroenie, 2001, 218 p.
2. Chernyshova T.I. Nistratov M.I. Transactions of the Tambov State Technical University, 2011, vol. 17, no. 1, pp. 24-31.
3. Chernyshova T.I. Ignatov D.V. Transactions of the Tambov State Technical University, 2005, vol. 11, no. 1B, pp. 241-244.
Mathematische Modellierung bei der Analyse der metrologischen Zuverlässigkeit der Analogblöcke der Informationsmeßsysteme
Zusammenfassung: Es ist das mathematische Modell der Veränderung in der Zeit der metrologischen Charakteristiken der Analogblöcke Informationsmeßsysteme unter Berücksichtigung der äußerlichen beeinflussenden Faktoren, sowie das verallgemeinerte mathematische Modell der Parameter der Elementbasis, verwendet in der Prozedur der statistischen Modellierung der metrologischen Charakteristiken der untersuchten Messmittel, zulassend entwickelt, ihr Niveau der metrologischen Zuverlässigkeit in verschiedenen Etappen des Betriebes zu analysieren. Als Hauptkennziffer der metrologischen Zuverlässigkeit der Analogblöcke der Informationsmeßsysteme ist die metrologische Ressource betrachtet. Es sind die Ergebnisse der Forschung der metrologischen Zuverlässigkeit des typisierten Blocks,
der den Messkanal der Informationsmeßsysteme bildet gebracht, es sind die Bedeutungen der metrologischen Ressource bei verschiedenen Bedeutungen der destabilisierenden Faktoren und in verschiedenen vorübergehenden Schnitten bekommen. Ein Hauptvorteil des angebotenen Herangehens an die mathematische Modellierung der metrologischen Charakteristiken ist die Möglichkeit der Bestimmung des Niveaus der metrologischen Zuverlässigkeit des Messmittels in einer beliebigen Etappe des Betriebes unter Berücksichtigung der gemeinsamen Handlung der sich ändernden äußerlichen Faktoren.
Modélage mathématique lors de l'analyse de la sécurité métrologique des unités analogiques des systèmes d'information et de mesure
Résumé: Est élaboré le modèle mathématique du changement dans le temps de la caractéristique des unités analogiques des systèmes d'information et de mesure compte tenu des facteurs extérieurs ainsi que le modèle mathématique généralisé des paramètres de la base élémentaire utilisée dans la procédure du modélage statistique des caractéristiques métrologiques des moyens de mesure étudiés permettant d'analyser le niveau de leur sécurité métrologique aux différentes étapes de l'exploitation. En qualité de l'indice essentielle de la sécurité métrologique des unités analogiques des systèmes d'information et de mesure est examinée une ressource métrologique. Sont cités les résultats de l'étude de la sécurité métrologique de l'unité type, sont obtenues les valeurs de la ressource métrologique. L'avantage principale de l'approche proposée envers le modélage mathématique des caractéristiques métrologiques est la possibilité de la définition du niveau de la sécurité métrologique du moyen de mesure à n'importe quelle étape de l'exploitation compte tenu l'action commune des facteurs extérieurs.
Авторы: Чернышова Татьяна Ивановна - доктор технических наук, профессор кафедры «Конструирование радиоэлектронных и микропроцессорных систем», директор института энергетики, приборостроения и радиоэлектроники; Третьяков Владимир Владиславович - аспирант кафедры «Конструирование радиоэлектронных и микропроцессорных систем», ФГБОУ ВПО «ТГТУ».
Рецензент: Пудовкин Анатолий Петрович - доктор технических наук, профессор, заведующий кафедрой «Радиотехника», ФГБОУ ВПО «ТГТУ».