Построение графика переходного процесса

Построение графика переходного процесса (рис.14) в системе по средствам программной среды MathCAD. Для этого требуется решить систему (40). Для решения преобразуем структуру системы к виду, показанному на рис.13 [1].

 

 

Рис.15. Преобразованная структура системы

Новой структуре соответствует система уравнений [1]:

                                              (40)

Первое уравнение является дифференциальным, второе – алгебраическим, так как содержит производные, находимые из первого уравнения.

Для построения графика переходного процесса (рис. 14) необходимо решить дифференциальное уравнение:

                                                          (41)

При численном решении дифференциального уравнения n-го порядка преобразуем в систему из n уравнений первого порядка [1].

                           z(t)=z₁(t), z^’(t)=z₂(t).                                       (42)

В результате подстановки выражения (42) в (40) можно записать систему уравнений следующего вида:       

                              (43)

Решив систему (43) определим:

                                            y(t) = 3.58z₁(t)                       (44)

Для решение системы в MathCAD необходимо указать следующие параметры:

Исходное уравнение:

                             (45)

Решение дифференциального уравнения методом Рунге-Кутта:

                                               (46)
где: - начальные условия.

n – 2000 число точек решения; tk = 200 конечное время решения; D – исходное уравнение.

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

                                                             (47)

Представим найденные значение в таблице 3.

Расчет точек для графика переходного процесса представим в таблице 4.

 

Таблица 3 Расчет точек переходного процесса

Таблица 4 Расчет точек для графика переходного процесса              

Рис. 16. График переходного процесса системы

 2.2.3 Оценка качества исследуемой системы

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

1. Вид переходного процесса

2. Длительность переходного процесса

3. Величина перерегулирования

4. Запас устойчивости системы по фазе

5. Запас устойчивости системы по амплитуде

Согласно графику переходного процесса (рис.14) переходный процесс имеет плавный апериодический характер.

Длительность переходного процесса определим как время, прошедшие от начала переходного процесса (t=0) до момента установления величины выходного сигнала, отличающейся не более чем на 5% от установившегося значения. По графику (рис.14)  переходного процесса определяем его длительность: t_пп = 0,72 с.

Величина перерегулирования определяется по графику переходного процесса и так как в системе отсутствуют колебания, то величина перерегулирования будет равна нулю.

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

Запас устойчивости системы по фазе определяется через фазовый угол системы на частоте среза ω_с:

 

                                                                (48)

где: - фазовый сдвиг, определяем по графику на чертеже КР-2068.998-26-06-00.00.000.Д1:

φ_з =180⁰ – 140⁰ = 40⁰                            (49)

Что бы система обладала достаточным качеством, запас устойчивости по фазе должен лежать в пределах 20˚- 50˚.

Запас устойчивости по амплитуде определяется по графику на чертеже КР-2068.998-26-06-00.00.000.Д1 как ордината ЛАХ на частоте фазового угла, равного π:

                                          (50)

                                                (51)

Для обеспечения перерегулирования в системе не более 20 % запас по амплитуде должен быть не менее -15дБ. Даная система удовлетворяет этому требованию.

Оценка точности системы

Оцениваются статическая и вынужденная ошибки системы. Для этого необходимо определить передаточную функцию замкнутой системы по ошибке [1]:

 Ф_х(р) = 1- Ф(р) =   =                                         (52)

Для уравнения (52) в общем виде:

                    Ф_х(р) = ,                                    (53)

где:  = 0;  = 1; = T_об = 0,05;  =0 ;  = k = 3,58;  = 1;  = T_об = 0,05;.

Представим Ф_х(р) в виде бесконечного степенного ряда [1]:

                              Ф_х(р) = S₀+ S₁·p+S₂·p²+…,                        (54)

где коэффициенты ошибок S_n определяются по формуле [1]:

                              S_n =  ( - )                        (55)

Определим коэффициенты ошибок:

                     S₀ =  = 0;                                         (56)
                           S₁ =  ( - S₀· )                            (57)

                     S₁=  (1- 0) = 0,05                     (58)

                              S₂ =  ( - (S₁· + S₀· ))                                   (59)

                                      S₂= (0,05- (0,05·1+0))= 0                              (60)

Находим выражение вынужденной ошибки системы:

                                  ∆_x(t) = S₀·x_вх(t) + S₁                            (61)

где: S₀-коэффициент статической ошибки,

 S₁-коэффициент скоростной ошибки,

 x_вх(t)-заданное значение выходного параметра системы.

                               ∆_x(t) = 0,05                      (62)

При x_вх(t) = 1 Dx(t) = 0 – статическая ошибка системы отсутствует (т.к. система астатическая)[1];

При x_вх(t) = t Dx(t) =0,05– скоростная ошибка;

⇐ Предыдущая123456Следующая ⇒

Поделиться: