×
09.06.2018
218.016.6054

Результат интеллектуальной деятельности: Способ поиска неисправностей в непрерывной динамической системе на основе введения пробных отклонений

Вид РИД

Изобретение

№ охранного документа
0002656923
Дата охранного документа
07.06.2018
Аннотация: Изобретение относится к диагностике систем автоматического управления. В способе поиска неисправностей в непрерывной динамической системе на основе введения пробных отклонений фиксируют неисправности, определяют время контроля и параметр преобразования сигналов. В качестве динамических характеристик используют интегральные оценки сигналов. Фиксируют контрольные точки системы, регистрируют реакцию объекта диагностирования и заведомо исправной системы и определяют оценки сигналов исправной системы. В качестве пробных отклонений рассматривают общее количество комбинаций структурных и топологических неисправностей. Определяют оценки выходных сигналов модели для каждой контрольной точки. Изменяют состояние каждой комбинации структурных блоков динамической системы. Регистрируют, определяют деформации оценок сигналов модели. Определяют нормированные значения деформаций оценок сигналов модели и диагностические признаки. По минимуму диагностического признака определяют комбинацию кратного структурного и топологического дефекта. Расширяются функциональные возможности способа. 1 ил.

Изобретение относится к области контроля и диагностирования систем автоматического управления и их элементов.

Известен способ поиска неисправного блока в динамической системе (Способ поиска неисправного блока в динамической системе: пат.2435189, Рос. Федерация: МПК7 G05B 23/02 (2006.01) / Шалобанов С.В., Шалобанов С.С. - №2009123999/08; заявл. 23.06.2009; опубл. 27.11.2011. Бюл. №33).

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

Наиболее близким техническим решением (прототипом) является способ поиска неисправностей в непрерывной динамической системе на основе введения пробных отклонений (Способ поиска неисправностей в непрерывной динамической системе на основе введения пробных отклонений: пат. 2541857 Рос. Федерация: МПК7 G05B 23/02 (2006.01) / Шалобанов С.С. - №2013149468/08; заявл. 06.11.2013; опубл. 20.02.2015, Бюл. №5).

Недостатком этого способа является то, что он позволяет находить только неисправности в виде изменения топологических связей, то есть обрыва или появления новых межблочных связей.

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

Поставленная задача достигается тем, что регистрируют реакцию заведомо исправной системы ƒj ном(t),j=1…k на интервале t∈[0, ТK] в k контрольных точках, и определяют интегральные оценки выходных сигналов Fjном(α), j=1,…,k, системы, для чего в момент подачи входного сигнала на вход системы с номинальными характеристиками одновременно начинают интегрирование сигналов системы управления в каждой из k контрольных точек с весами e-αt, где путем подачи на первые входы k блоков перемножения сигналов системы управления, на вторые входы блоков перемножения подают экспоненциальный сигнал e-αt, выходные сигналы k блоков перемножения подают на входы k блоков интегрирования, интегрирование завершают в момент времени Тк, полученные в результате интегрирования оценки выходных сигналов Fjном(α), j=1,…,k, регистрируют, определяют интегральные оценки выходных сигналов модели для каждой из k контрольных точек, полученные в результате введения пробных отклонений каждой из m возможных комбинаций структурных и топологических состояний (добавляется дефект/дефекты в структурный блок/блоки и/или удаляется существующая межблочная связь/связи либо вводится новая межблочная связь/связи), для чего поочередно для каждой возможной комбинации динамических блоков и топологических связей системы вводят пробные отклонения параметров структурных блоков и/или состояний топологических связей и находят интегральные оценки выходных сигналов системы для параметра α и входного сигнала x(t), полученные в результате интегрирования оценки выходных сигналов для каждой из k контрольных точек и каждого из m пробных отклонений Pji(α), j=1,…,k; i=1,…,m, регистрируют, определяют деформации интегральных оценок выходных сигналов модели, полученные в результате пробных отклонений параметров соответствующих структурных блоков и/или состояний топологических связей блоков динамической системы ΔРji(α)=Рji(α)-Fjном(α), j=1,…,k; i=1,…,m, определяют нормированные значения деформаций интегральных оценок выходных сигналов модели, полученные в результате пробных отклонений параметров соответствующих структурных блоков и/или состояний топологических связей из соотношения , замещают систему с номинальными характеристиками контролируемой, на вход системы подают аналогичный входной сигнал x(t), определяют интегральные оценки выходных сигналов контролируемой системы для k контрольных точек Fj(α), j=1,…,k для параметра α, определяют отклонения интегральных оценок выходных сигналов контролируемой системы для k контрольных точек от номинальных значений ΔFj(α)=Fj(α)-Fjном(α), j=1,…, k, определяют нормированные значения отклонений интегральных оценок выходных сигналов контролируемой системы из соотношения , определяют диагностические признаки из соотношения , i=1,…, m, по минимуму диагностического признака определяют кратный структурный и/или топологический дефект.

Таким образом, предлагаемый способ поиска комбинации дефектов структурных блоков и неисправных топологических связей блоков сводится к выполнению следующих операций:

1. В качестве динамической системы рассматривают систему, состоящую из произвольно соединенных динамических блоков, с количеством рассматриваемых комбинаций структурных и топологических состояний m.

2. Предварительно определяют время контроля ТК≥ТПП, где ТПП - время переходного процесса системы. Время переходного процесса оценивают для номинальных значений параметров динамической системы.

3. Определяют параметр интегрального преобразования сигналов из соотношения .

4. Фиксируют число контрольных точек k.

5. Предварительно определяют нормированные векторы деформаций интегральных оценок выходных сигналов модели, полученные в результате пробных отклонений параметров структурных блоков и/или состояний топологических связей между блоками каждой из m комбинации динамических блоков и топологических связей системы для номинальных состояний комбинации динамических блоков и топологических связей системы и определенного выше параметра α, для чего выполняют пункты 6-10.

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

7. Регистрируют реакцию системы fj ном(t), j=1,…,k, на интервале t∈[0, ТК] в k контрольных точках и определяют интегральные оценки выходных сигналов Fjном(α), j=1,…,k, системы. Для этого в момент подачи тестового сигнала на вход системы управления с номинальными характеристиками одновременно начинают интегрирование сигналов системы управления в каждой из k контрольных точек с весами е-αt, где , для чего сигналы системы управления подают на первые входы k блоков перемножения, на вторые входы блоков перемножения подают экспоненциальный сигнал e-αt, выходные сигналы k блоков перемножения подают на входы к блоков интегрирования, интегрирование завершают в момент времени Тк, полученные в результате интегрирования оценки выходных сигналов Fjном (α), j=1,…, k регистрируют.

8. Определяют интегральные оценки выходных сигналов модели для каждой из k контрольных точек, полученные в результате каждой из m комбинаций динамических блоков и топологических связей системы, для чего поочередно изменяют состояние каждой комбинации динамических блоков и топологических связей динамической системы (например, изменяют параметры одного или нескольких блоков и/или меняют связь одного или нескольких блоков с соседними из состояния «есть связь» в состояние «нет связи» или наоборот) и выполняют пункты 6 и 7 для одного и того же входного сигнала x(t). Полученные в результате интегрирования оценки выходных сигналов для каждой из k контрольных точек и каждого из m пробных отклонений Pji(α), j=1,…, k; i=1,…,m, регистрируют.

9. Определяют деформации интегральных оценок выходных сигналов модели, полученные в результате пробных отклонений комбинаций динамических блоков и топологических связей динамической системы ΔPji(α)=Pji(α)-Fjном(α), j=1,…, k; i=1,…, m.

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

11. Замещают систему с номинальными характеристиками контролируемой. На вход системы подают аналогичный тестовый сигнал x(t).

12. Определяют интегральные оценки выходных сигналов контролируемой системы для k контрольных точек Fj(α), j=1,…,k, осуществляя операции, описанные в пунктах 6 и 7 применительно к контролируемой системе.

13. Определяют отклонения интегральных оценок выходных сигналов контролируемой системы для k контрольных точек от номинальных значений ΔFj(α)=Fj(α)-Fj ном(α), j=1,…,k.

14. Вычисляют нормированные значения отклонений интегральных оценок выходных сигналов контролируемой системы по формуле .

15. Вычисляют диагностические признаки наличия неисправных структурных блоков и/или топологических связей по формуле

, i=1,…,m.

16. По минимуму значения диагностического признака определяют кратный структурный и/или топологический дефект.

Рассмотрим реализацию предлагаемого способа поиска комбинации структурных и топологических дефектов для системы, схема которой представлена на рисунке (см. чертеж. Структурная схема объекта диагностирования).

Передаточные функции блоков:

; ; ,

номинальные значения параметров: T1=5 с; k1=1; k2=1; Т2=1 с; k3=1; Т3=5 с. При поиске комбинации структурного и топологического дефекта в виде отклонения постоянной времени T1=4 с в первом звене (дефект №1), путем подачи ступенчатого тестового входного сигнала единичной амплитуды и интегрального преобразования сигналов для параметра α=0.5 и Тк=10 с получены значения диагностических признаков на основе пробных отклонений параметров модели при использовании трех контрольных точек, расположенных на выходах блоков: J1=0 (дефект в первом блоке); J2=0.7904 (обрыв связи между первым и вторым блоком); J3=0.8862 (дефект в первом блоке и обрыв связи между первым и вторым блоком); J4=0.08986 (дефект во втором блоке и обрыв связи между третьим и первым блоком); J5=0.0863 (дефект во втором и третьем блоке, а также обрыв связи между вторым и третьим, третьим и первым блоком). Минимальное значение признака J1 однозначно указывает на наличие дефекта в первом блоке, а разность между пятым и первым признаками может количественно характеризовать практическую (апостериорную) различимость этого дефекта.

При поиске комбинации структурного и топологического дефекта в виде обрыва связи между первым и вторым звеньями (дефект №2) для данного объекта диагностирования, путем подачи ступенчатого тестового входного сигнала единичной амплитуды, при том же параметре интегрирования α и при том же времени контроля Тк, получены значения диагностических признаков на основе пробных отклонений состояний топологической связи при использовании трех контрольных точек, расположенных на выходах блоков: J1=0.7904 (дефект в первом блоке); J2=0 (обрыв связи между первым и вторым блоком); J3=0.007189 (дефект в первом блоке и обрыв связи между первым и вторым блоком); J4=0.7086 (дефект во втором блоке и обрыв связи между третьим и первым блоком); J5=0.6724 (дефект во втором и третьем блоке, а также обрыв связи между вторым и третьим, третьим и первым блоком). Минимальное значение признака J2 однозначно указывает на изменение топологической связи между первым и вторым блоком.

Моделирование процессов поиска комбинации структурного (в виде отклонения постоянной времени T1=4 с в первом звене) и топологического дефекта в виде обрыва связи между первым и вторым звеньями (дефект №3) для данного объекта диагностирования, при том же параметре интегрирования α и при единичном ступенчатом входном сигнале дает следующие значения диагностических признаков на основе пробных отклонений параметров модели и состояний топологической связи при использовании трех контрольных точек, расположенных на выходах блоков: J1=0.8862 (дефект в первом блоке); J2=0.01715 (обрыв связи между первым и вторым блоком); J3=0 (дефект в первом блоке и обрыв связи между первым и вторым блоком); J4=0.8194 (дефект во втором блоке и обрыв связи между третьим и первым блоком); J5=0.7836 (дефект во втором и третьем блоке, а также обрыв связи между вторым и третьим, третьим и первым блоком).

Моделирование процессов поиска комбинации структурного (в виде отклонения коэффициента усиления k2=0.8 во втором звене) и топологического дефекта в виде обрыва связи между третьим и первым звеньями (дефект №4) для данного объекта диагностирования, при том же параметре интегрирования α и при единичном ступенчатом входном сигнале дает следующие значения диагностических признаков на основе пробных отклонений параметров модели и состояний топологической связи при использовании трех контрольных точек, расположенных на выходах блоков: J1=0.03641 (дефект в первом блоке); J2=0.617 (обрыв связи между первым и вторым блоком); J3=0.6974 (дефект в первом блоке и обрыв связи между первым и вторым блоком); J4=0 (дефект во втором блоке и обрыв связи между третьим и первым блоком); J5=0.06111 (дефект во втором и третьем блоке, а также обрыв связи между вторым и третьим, третьим и первым блоком).

Моделирование процессов поиска комбинации структурного (в виде отклонения коэффициентов усиления k2=0.8 во втором звене и k3=0.8 в третьем звене) и комбинации топологического дефекта в виде обрывов связей между вторым и третьим, а также третьим и первым звеньями (дефект №5) для данного объекта диагностирования, при том же параметре интегрирования α и при единичном ступенчатом входном сигнале дает следующие значения диагностических признаков на основе пробных отклонений параметров модели и состояний топологической связи при использовании трех контрольных точек, расположенных на выходах блоков: J1=0.1138 (дефект в первом блоке); J2=0.5813 (обрыв связи между первым и вторым блоком); J3=0.6599 (дефект в первом блоке и обрыв связи между первым и вторым блоком); J4=0.07486 (дефект во втором блоке и обрыв связи между третьим и первым блоком); J5=0 (дефект во втором и третьем блоке, а также обрыв связи между вторым и третьим, третьим и первым блоком).

Минимальное значение диагностического признака во всех случаях правильно указывает на наличие топологического дефекта.

Способ поиска неисправностей в непрерывной динамической системе на основе введения пробных отклонений, основанный на том, что фиксируют число возможных неисправностей m, определяют время контроля Т≥Т, где Т - время переходного процесса системы, определяют параметр интегрального преобразования сигналов из соотношения , используют тестовый сигнал на интервале t∈[0, T], в качестве динамических характеристик системы используют интегральные оценки сигналов, полученные для вещественных значений параметра интегрального преобразования α, фиксируют число k контрольных точек системы, регистрируют реакцию объекта диагностирования ƒ(t), j=1, …, k, и реакцию заведомо исправной системы ƒ(t), j=1, …, k, на интервале t∈[0, T] в k контрольных точках, определяют интегральные оценки выходных сигналов F(α), j=1, …, k, исправной системы, для чего в момент подачи тестового сигнала на вход системы с номинальными характеристиками одновременно начинают интегрирование выходных сигналов системы управления в каждой из k контрольных точек с весами e, где , путем подачи на первые входы k блоков перемножения выходных сигналов системы управления, на вторые входы блоков перемножения подают экспоненциальный сигнал e, выходные сигналы k блоков перемножения подают на входы k блоков интегрирования, интегрирование завершают в момент времени Т, полученные в результате интегрирования оценки выходных сигналов F(α), j=1, …, k, регистрируют, определяют интегральные оценки выходных сигналов модели для каждой из k контрольных точек, полученные в результате каждого из m пробных отклонений и находят интегральные оценки выходных сигналов системы для параметра α и входного сигнала x(t), полученные в результате интегрирования оценки выходных сигналов для каждой из k контрольных точек и каждого из m пробных отклонений P(α), j=1, …, k; i=1, …, m, регистрируют, замещают систему с номинальными характеристиками контролируемой, на вход системы подают аналогичный тестовый сигнал x(t), определяют интегральные оценки выходных сигналов контролируемой системы для k контрольных точек F(α), j=1, …, k, для параметра α, определяют отклонения интегральных оценок выходных сигналов контролируемой системы для k контрольных точек от номинальных значений ΔF(α)=F(α)-F(α), j=1, …, k, определяют нормированные значения отклонений интегральных оценок выходных сигналов контролируемой системы из соотношения , отличающийся тем, что в качестве количества пробных отклонений m рассматривают общее количество возможных комбинаций структурных и топологических неисправностей, определяют интегральные оценки выходных сигналов модели для каждой из k контрольных точек, полученные в результате каждого из m пробных отклонений параметров модели и/или состояний топологических связей, для чего поочередно изменяют состояние каждой комбинации структурных блоков и/или топологических связей блоков динамической системы, для которых находят интегральные оценки выходных сигналов системы для параметра α и входного сигнала x(t), полученные в результате интегрирования оценки выходных сигналов для каждой из k контрольных точек и каждого из m пробных отклонений параметров структурных блоков модели и/или состояний топологических связей P(α), j=1, …, k; i=1, …, m, регистрируют, определяют деформации интегральных оценок выходных сигналов модели, полученные в результате пробных отклонений параметров соответствующих структурных блоков и/или состояний топологических связей динамических блоков системы ΔP(α)=P(α)-F(α), j=1, …,k; i=1, …, m, определяют нормированные значения деформаций интегральных оценок выходных сигналов модели, полученные в результате пробных отклонений параметров соответствующих блоков модели и/или состояний топологических связей блоков из соотношения , определяют диагностические признаки из соотношения , i=1, …, m, по минимуму диагностического признака определяют комбинацию кратного структурного и топологического дефекта.
Способ поиска неисправностей в непрерывной динамической системе на основе введения пробных отклонений
Способ поиска неисправностей в непрерывной динамической системе на основе введения пробных отклонений
Способ поиска неисправностей в непрерывной динамической системе на основе введения пробных отклонений
Способ поиска неисправностей в непрерывной динамической системе на основе введения пробных отклонений
Способ поиска неисправностей в непрерывной динамической системе на основе введения пробных отклонений
Способ поиска неисправностей в непрерывной динамической системе на основе введения пробных отклонений
Источник поступления информации: Роспатент

Показаны записи 41-50 из 115.
19.01.2018
№218.016.0ab2

Способ токарной обработки заготовок из углепластика

Изобретение относится к токарной обработке материалов и может быть использовано при механической обработке заготовок из полимерных композиционных материалов, преимущественно из углепластика. Техническим результатом является повышение производительности обработки заготовки и качества...
Тип: Изобретение
Номер охранного документа: 0002632299
Дата охранного документа: 03.10.2017
19.01.2018
№218.016.0ab8

Жидкостно-газовый эжектор

Изобретение относится к области струйной техники, преимущественно к струйным аппаратам для создания вакуума. В эжекторе, содержащем распределительную камеру с соплами, приемную камеру, камеры смешения и сбросную камеру. Каждая камера смешения установлена соосно относительно своего сопла. Сопло...
Тип: Изобретение
Номер охранного документа: 0002632167
Дата охранного документа: 03.10.2017
19.01.2018
№218.016.0b5e

Способ получения материалов на основе алюминида никеля

Изобретение относится к получению материала на основе алюминида никеля. Способ включает приготовление экзотермической шихты путем смешивания порошков алюминия, оксида никеля и по крайней мере одной легирующей добавки и инициирование в экзотермической шихте металлотермической реакции с...
Тип: Изобретение
Номер охранного документа: 0002632341
Дата охранного документа: 04.10.2017
20.01.2018
№218.016.16b5

Способ сортирования волокнистых материалов

Изобретение относится к целлюлозно-бумажному производству и может быть использовано для сортирования волокнистых материалов. Способ сортирования волокнистых материалов, включающий загрузку волокнистого материала через входное центральное отверстие в цапфе во внутреннюю полость сортирующего...
Тип: Изобретение
Номер охранного документа: 0002635792
Дата охранного документа: 15.11.2017
20.01.2018
№218.016.1892

Жидкостно-газовый эжекторный аппарат

Изобретение относится к области струйной техники, преимущественно к струйным аппаратам для создания вакуума. Аппарат содержит распределительную камеру с соплами, приемную камеру, камеры смешения и сбросную камеру, причем каждая камера смешения установлена соосно относительно своего сопла. Сопло...
Тип: Изобретение
Номер охранного документа: 0002635424
Дата охранного документа: 13.11.2017
20.01.2018
№218.016.1a42

Роторно-пульсационный аппарат

Изобретение относится к оборудованию для получения дисперсных систем, преимущественно "жидкость - жидкость", и может быть использовано в химической, пищевой, микробиологической и других отраслях промышленности. Роторно-пульсационный аппарат содержит рабочую камеру с устройством ввода...
Тип: Изобретение
Номер охранного документа: 0002636486
Дата охранного документа: 23.11.2017
04.04.2018
№218.016.3184

Сейсмоакустический преобразователь

Изобретение относится к геофизике, в частности к сейсмоакустическим исследованиям, и может быть использовано для получения прогностических характеристик при контроле трещинообразования в массиве горных пород. Заявлен сейсмоакустический преобразователь, у которого дополнительно в корпусе на...
Тип: Изобретение
Номер охранного документа: 0002645037
Дата охранного документа: 15.02.2018
04.04.2018
№218.016.326d

Смеситель для смешивания жидких и вязких материалов

Изобретение относится к устройству для смешивания жидких и вязких материалов и может найти применение в химической, пищевой, фармацевтической, строительной и других отраслях промышленности. Смеситель содержит цилиндрический неподвижный корпус в виде стакана с размещенным внутри него...
Тип: Изобретение
Номер охранного документа: 0002645537
Дата охранного документа: 21.02.2018
04.04.2018
№218.016.329b

Генератор для получения талой питьевой воды

Изобретение относится к устройствам для доочистки воды. Генератор для получения талой питьевой воды, включающий расположенные последовательно в одном продольном сосуде зону замораживания воды с кольцевой морозильной камерой, зону вытеснения примесей из фронта льда и концентрации примесей в виде...
Тип: Изобретение
Номер охранного документа: 0002645541
Дата охранного документа: 21.02.2018
04.04.2018
№218.016.32a7

Смеситель для жидких и вязких материалов

Изобретение относится к устройству для смешивания жидких и вязких материалов и может найти применение в химической, пищевой, фармацевтической, строительной и других отраслях промышленности. Смеситель для жидких и вязких материалов содержит цилиндрический неподвижный корпус в виде стакана, с...
Тип: Изобретение
Номер охранного документа: 0002645521
Дата охранного документа: 21.02.2018
Показаны записи 21-25 из 25.
02.03.2019
№219.016.d1b9

Способ поиска неисправного блока в непрерывной динамической системе на основе функции чувствительности

Изобретение относится к способу поиска неисправного блока в непрерывной динамической системе на основе функции чувствительности. Для поиска неисправного блока в системе из произвольно соединенных динамических блоков предварительно определяют время контроля, определяют параметр интегрального...
Тип: Изобретение
Номер охранного документа: 0002680928
Дата охранного документа: 28.02.2019
16.01.2020
№220.017.f52d

Способ поиска неисправностей динамического блока в непрерывной системе на основе введения пробных отклонений

Изобретение относится к диагностике систем автоматического управления. В способе поиска неисправностей динамического блока в непрерывной системе на основе введения пробных отклонений определяют время контроля; определяют число контрольных точек; определяют векторы отклонений сигналов модели в...
Тип: Изобретение
Номер охранного документа: 0002710998
Дата охранного документа: 14.01.2020
17.01.2020
№220.017.f6d9

Способ поиска неисправного блока в непрерывной динамической системе на основе введения пробных отклонений и анализа знаков передач

Изобретение относится к диагностике систем автоматического управления. В способе поиска неисправного блока в непрерывной динамической системе на основе введения пробных отклонений и анализа знаков передач, фиксируют число блоков системы, определяют время контроля и параметр интегрального...
Тип: Изобретение
Номер охранного документа: 0002711000
Дата охранного документа: 14.01.2020
24.04.2020
№220.018.184a

Способ поиска неисправного блока в непрерывной динамической системе на основе функции чувствительности

Изобретение относится к диагностике систем автоматического управления. В способе поиска неисправного блока в непрерывной динамической системе на основе функции чувствительности определяют знаки интегральных оценок выходных сигналов модели, полученные в результате структурной функции...
Тип: Изобретение
Номер охранного документа: 0002719747
Дата охранного документа: 22.04.2020
21.05.2020
№220.018.1f63

Способ поиска неисправного блока в непрерывной динамической системе на основе смены позиции входного сигнала и анализа знаков передач

Изобретение относится к диагностике систем автоматического управления. В способе поиска неисправного блока в непрерывной динамической системе на основе смены позиции входного сигнала и анализа знаков передач определяют знаки отклонений интегральных оценок выходных сигналов контролируемой...
Тип: Изобретение
Номер охранного документа: 0002721217
Дата охранного документа: 18.05.2020
+ добавить свой РИД