Большая База Данных Рефератов - Оценка параметрической надежности РЭС с использованием моделирования на ЭВМ постепенных отказов

Рус Укр Eng

Оценка параметрической надежности РЭС с использованием моделирования на ЭВМ постепенных отказов

МГТУ ГА

Факультет: компьютерного проектирования

Кафедра: радиоэлектронных средств

Пояснительная записка к курсовому проекту

по предмету: «Теоретические основы конструирования, технологии и надежности»

на тему: «Оценка параметрической надежности РЭС с использованием моделирования на ЭВМ постепенных отказов»

Москва 2002

СОДЕРЖАНИЕ

Введение

1. Постановка задачи

1.1 Анализ исходных данных

1.2 Пояснение решаемой задачи

2. Выбор метода решения поставленной задачи

3. Решение задачи на ЭВМ

3.1 Описание вычислительного алгоритма моделирования температурных и временных изменений параметров.

3.2 Пояснение процедур и функций, используемых в программе

3.3 Обоснование выбора числа реализаций

3.4 Список идентификаторов

4. Описание и анализ полученных результатов

5. Пояснения функциональных частей структурной схемы алгоритма

Заключение и выводы.

Литература.

Приложение 1. Листинг программы.

Приложение 2. Графический материал.

ВВЕДЕНИЕ

В курсовом проекте необходимо произвести оценку параметрической надёжности РЭС, с использованием моделирования на ЭВМ постепенных отказов элементов.

Параметрическая надёжность РЭУ - вероятность отсутствия в изделии постепенных отказов при его работе в заданных условиях эксплуатации в течение времени tзад (в нашем случае tзад = 10000 ч). Параметрическая надёжность связана с понятием постепенных отказов.

Постепенный (параметрический) отказ - отказ, возникающий в результате постепенного изменения значения одного или нескольких параметров изделия.

Основные причины, вызывающие появление постепенных отказов:

1) Производственный разброс выходного параметра, вызываемый действием производственных погрешностей.

2) Отклонение выходного параметра от номинального значения из-за процессов старения.

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

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

Моделируя РЭУ и используя методы математической статистики, проследим как влияют производственный разброс входных параметров, дестабилизирующие факторы и старение на выходной параметр, а следовательно и на параметрическую надежность.

1. ПОСТАНОВКА ЗАДАЧИ

1.1 Анализ исходных данных

Исходные данные к проекту:

1) Схема электрическая принципиальная.

2) Математическая модель для выходного параметра:

(1.1)

3) Сведения о первичных параметрах (параметрах элементов):

а) резисторы R1 = 3 кОм 5% типа ОМЛТ;

б) резисторы R2 = 12 кОм 5% типа ОМЛТ;

в) резисторы R3 = 2,4 кОм 10% типа ОМЛТ;

г) тип микросхемы DA1: 140УД9;

4) Заданное интервал работы РЭС: tзад = 10000 час.

5) Диапазон рабочих температур: Траб = +10…+60 С.

6) Условие параметрической надежности:

Данных, указанных в задании, недостаточно для проведения расчетов и моделирования. Поэтому дополняем необходимые данные из справочников:

7) Согласно [3] температурный коэффициент резисторов типа ОМЛТ:

а) R+ = 710-2 % при Т = +20…+100 С;

б) R- = 1210-2 % при Т = -60… +20 С;

8) Согласно [3] на резисторы типа ОМЛТ величина их сопротивления может измениться на 10% при наработке 25000 часов. Отсюда находим величину коэффициента старения:

СR = = 410-4 % ;

9) Согласно [2] коэффициент усиления Koy и входное сопротивлениеRbx:

Koy35000

Rbx300 кОм

Характеристики первичных параметров представлены в неявной форме, т. е. нет численных значений математического ожидания М(xi) и среднеквадратического отклонения (xi).Вследствие этого необходимо произвести их расчет.

Расчет этих характеристик производят в зависимости от закона распределения первичного параметра. Примем гипотезу о том, что Koy и Rbx распределены по нормальному закону. (Koy)

35000 М(Koy) Koy

Согласно [1] составим систему уравнений:

Koy=5000030%

Аналогично определяем Rbx .Получаем Rbx=430 кОм30%.

Т.о. получили Koy=5000030% Rbx=430 кОм30%

10) На основе данных, приведённых в [2] получили стабильность Koy и Rbx :

а)Температурная : Koy= 2510-2 % при Т = -60…+100 С;

Rbx = 7,510-3 % при Т = -60…+100 С;

б)Временная: С Koy= 310-3%; С Rbx= 510-4 % ;

11) Коэффициент корреляции между Koy и Rbx: r =0.8

1.2 Пояснение решаемой задачи

Оценка параметрической надёжности - определение основных количественных показателей сохранения рабочих функций при возможных постепенных изменениях параметров комплектующих элементов в условиях эксплуатации.

Оценку параметрической надежности будем проводить следующим способом: Подсчитав по формуле (1.1) выходной параметр K (коэффициент передачи) и установив допуск на выходной параметр K, смоделируем n РЭУ. РЭУ будем считать работоспособным, если значение коэффициента передачи лежит в диапазоне установленного допуска, т.е. K K. Таким образом, найдём вероятность отсутствия параметрического отказа (см. раздел 2).

2. ВЫБОР МЕТОДА РЕШЕНИЯ ПОСТАВЛЕННОЙ ЗАДАЧИ

Метод решения задачи состоит в следующем. Определяем выходной параметр по формуле (1.1) по значениям параметров элементов, не учитывая производственные допуска, корреляцию, воздействия температуры и времени. Назовем полученный таким образом коэффициент передачи “идеальным” -- Kи. После чего задаемся допуском на выходной параметр Kи, в пределах которого РЭУ считается исправным.

При помощи ЭВМ моделируем n различных реализаций РЭУ с параметрами элементов, распределенных либо по нормальному закону, либо по равномерному закону. Затем пересчитываем значения параметров элементов при воздействии на них температуры и времени. При этом предполагаем, что температурный коэффициенты R, а также коэффициенты старения СR распределены по нормальному закону, а температура окружающей среды Траб - по равномерному. В связи с тем, что закон распределения температуры окружающей среды был неизвестен, и не было возможности попытаться подобрать закон распределения экспериментально, то была принята гипотеза о том, что температура распределена по равномерному закону, так как эта модель на практике является предельным (наихудшим) случаем разброса параметра. Определяем выходной параметр по формуле (1.1) - этот коэффициент передачи назовем “реальным”(Kр).

По способу, изложенному в подразделе 1.2, вероятность отсутствия параметрического отказа определим следующим образом:

Р (Kн Kр Kв tзад)= , (2.1)

где nисп - число исправных РЭУ на момент времени tзад;

N - общее число смоделированных РЭУ;

Kн - нижнее значение коэффициента передачи Kн = Kи - Kи;

Kв - верхнее значение коэффициента передачи Kв = Kи + Kи.

Определяем математическое ожидание выходного параметра М*(Kр) и его среднеквадратичное отклонение *(Kр) по формулам [1]:

М*(Kр) = , (2.2)

(2.3)

3. РЕШЕНИЕ ЗАДАЧИ НА ЭВМ

3.1 Описание вычислительного алгоритма моделирования температурных и временных изменений параметров

R1, R2, R3 - сопротивления 1-го, 2-го и 3-го резисторов;

Rbx - входное сопротивление, Koy - коэффициент усиления.

1. При помощи стандартной функции Random генерируем равномерно распределённое значение температуры: temp.

Здесь вычислительный алгоритм разделяется на 2 части:

а) Если температура попала в положительную область диапазона рабочих температур т.е 20,

то, используя формулу (3.1) [1] генерируем нормально распределённые значения температурных коэффициентов R+, Rbx : dx1,dx2,dx3,dx 4.

R+ - температурный коэффициент для резисторов в полож-й области температур;

Rbx - температурный коэффициент для входного сопротивления.

dx1, dx2, dx3, dx4 - сгенерированные значения температурных коэффициентов для

1-го, 2-го, 3-го резисторов и входного сопротивления соответственно.

б)Если температура попала в отрицательную область диапазона рабочих температур т.е 20,

R- - температурный коэффициент для резисторов в отриц-й области температур;

Rbx - температурный коэффициент для входного сопротивления.

dx1, dx2, dx3, dx4 - сгенерированные значения температурных коэффициентов для

1-го, 2-го, 3-го резисторов и входного сопротивления соответственно.

x = + m, (3.1)

где x - нормально распределённое случайное число;

m - математическое ожидание;

- среднеквадратичное отклонение;

ri - стандартное равномерно распределенное случайное число в диапазоне 0..1. (ri получаем при помощи стандартной функции Random).

Далее пересчитываем значения первичных параметров (R1,R2, R3, Rbx) с учётом воздействия температуры. Для этого воспользуемся формулами [1]:

(3.2)

где - номинальные значения i-го первичного параметра;

-приращения значений i-го первичного параметра под действием температуры;

Согласно [1] относительное изменение i-го первичного параметра под воздействием температуры (старения) можно выразить следующим образом:

(3.3)

(3.4)

где - температурный коэффициент i-го первичного параметра;

где tср - температура окружающей среды;

сi - коэффициент старения i-го первичного параметра;

- рассматриваемый интервал времени.

В качестве tср для положительной области диапазона рабочих температур примем

наибольшую из возможных температур - Tv, а для отрицательной области примем наименьшую из возможных температур - Tn. С учётом этого и формул (3.3) и (3.4) формула (3.2) примет вид:

для ``+`` -ой области температур:

(3.5)

С учётом этой формулы получаем:

;; ;;

для ``-`` -ой области температур:

(3.6)

С учётом этой формулы получаем:

;;

где Rtemp1, Rtemp2, Rtemp3 - значения сопротивлений 1-го, 2-го и 3-го резисторов соответственно с учётом действия температуры.

RWtemp - значение входного сопротивления под действием температуры.

SR1, SR2, SR3 - номинальные значения 1-го, 2-го и 3-го резисторов соответственно.

SRW - номинальное значение входного сопротивления.

Для получения значений коэффициента усиления (Koy) производим смещение параметров m = m(z) и = (z) его температурного коэффициента ( Koy) с учётом коэффициента парной корреляции , а затем, воспользовавшись подпрограммой формирования случайных нормально распределённых чисел с параметрами m = m(z/x) и = (z/x) генерируем нормально распределённое значение его температурного коэффициента( Koy):dx5.

dx5 - сгенерированное значение температурного коэффициента для коэффициента усиления.

Воспользовавшись формулой (3.5) (для положительной области температур) или (3.6) (для отрицательной области температур) пересчитываем значения коэффициента усиления (Koy) с учётом воздействия температуры:

для ``+`` -ой области температур:

;

для ``-`` -ой области температур: ;

где KOUtemp - значение коэффициента усиления под действием температуры.

SKOU - номинальное значение коэффициента усиления.

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

2. Используя формулу (3.1) генерируем нормально распределённые значения коэффициентов старения СR, С Rbx :dx1,dx2,dx3,dx4.

СR - коэффициент старения для резисторов;

С Rbx - коэффициент старения для входного сопротивления;

dx1, dx2, dx3, dx4 - сгенерированные значения коэффициентов старения для

1-го, 2-го, 3-го резисторов и входного сопротивления соответственно.

Воспользовавшись формулой:

(3.7)

пересчитываем значения первичных параметров (R1,R2, R3, Rbx) с учётом воздействия старения:

;;

где Rtime1, Rtime2, Rtime3 - значения сопротивлений 1-го, 2-го и 3-го резисторов соответственно с учётом действия старения.

RWtime - значение входного сопротивления под действием старения.

SR1, SR2, SR3 - номинальные значения 1-го, 2-го и 3-го резисторов соответственно.

SRW - номинальное значение входного сопротивления.

Для получения значений коэффициента усиления (Koy) производим смещение параметров m = m(z) и = (z) его коэффициента старения(С Koy) с учётом коэффициента парной корреляции , а затем, воспользовавшись подпрограммой формирования случайных нормально распределённых чисел с параметрами m = m(z/x) и = (z/x) генерируем нормально распределённое значение его коэффициента старения(С Koy):dx5.

Воспользовавшись формулой (3.7) пересчитываем значения коэффициента усиления (Koy) с учётом воздействия старения:

;

где KOUtime - значение коэффициента усиления под действием температуры.

SKOU - номинальное значение коэффициента усиления.

По формуле (1.1) определяем значение выходного параметра: коэффициента передачи (Kexit).

3.2 Пояснение процедур и функций, используемых в программе
В написанной программе формула (3.1) реализована через функцию:

Function Generator(m:Real;s:Real):Real;

Label L1;

BEGIN

L1:x:=0;

FOR i:=1 TO 12 DO

BEGIN

k:=Random;

x:=x+k;

END;

x:=x-6;

if (x>3) or (x<-3) then goto L1;

m:=m+s*x;

Generator:=m;

END;

Таким образом, введя Generator(m,s)получим случайное число, распределенное по нормальному закону с параметрами m = m и = s.

В соответствии с [1] формула получения случайных чисел, распределенных по равномерному закону с параметрами a и b следующая:

x = r+ a, (3.8)

где a, b - параметры равномерной модели;

r -стандартное равномерно распределенное случайное число в диапазоне 0..1.

В написанной программе формула (3.8) реализована через функцию:

Function Generator2(m:real;s:real):Real;

BEGIN

k:=Random;

m:=(s-m)*k+m;

Generator2:=m;

end;

Таким образом, введя Generator2(m, s)получим случайное число, распределенное по равномерному закону с параметрами a=m и b = s.

Пусть случайное число x, имеющее нормальное распределение с параметрами m = m(x) и = (x), уже получено. Тогда для получения случайного числа z, имеющего нормальное распределение с параметрами m = m(z) и = (z) и коррелированного с x, необходимо произвести смещение параметров m = m(z) и = (z) с учётом коэффициента парной корреляции, а затем воспользоваться подпрограммой формирования случайных нормально распределённых чисел с параметрами m = m(z/x) и = (z/x):

(3.9)

(3.10)

Определение величины смещения параметров m = M(z) и = (z) с учётом коэффициента парной корреляции в соответствии с формулами (3.9) и (3.10) в программе реализовано следующим образом:

Procedure Corr(x1,mx,mz,sx,sz:real; Var mzx,szx:real);

BEGIN

mzx:=mz+rxz*(sz/sx)*(x1-mx);

szx:=sz*sqrt(1-sqr(rxz));

END;

Страницы: 1, 2

Меню

Оценка параметрической надежности РЭС с использованием моделирования на ЭВМ постепенных отказов

Оценка параметрической надежности РЭС с использованием моделирования на ЭВМ постепенных отказов