Выбор параметров LC-фильтра инвертора с широтно-импульсной модуляцией и синусоидальным выходным напряжением

№ 1’2011
PDF версия
В данной работе приводится анализ процессов с использованием подхода локального усреднения и с дополнительным учетом некоторых реальных параметров элементов.

В инверторах с широтно-импульсной модуляцией (ШИМ) имеют место сложные процессы, не поддающиеся точному аналитическому описанию. Напряжения и токи ветвей схемы инвертора, кроме составляющих основной частоты, содержат пульсации сложной формы, усложняющие точный расчет. К примеру, на рис. 1 показаны осциллограммы тока дросселя и напряжения конденсатора фильтра при ШИМ, подтверждающие это.

 Осциллограммы токов LC-фильтра при ШИМ

Рис. 1. Осциллограммы токов LC-фильтра при ШИМ

За последнее время инверторы с ШИМ получили широкое практическое применение благодаря хорошим техническим, эксплуатационным и другим характеристикам. Накоплена достаточная база для их разработки и использования. Расчету процессов в таких инверторах посвящена обширная литература [4, 5, 8, 12]. Известны разнообразные методики и подходы к вычислениям, начиная от самых грубых и оценочных и кончая уточненными с интегральным учетом высших гармоник. На современном этапе разработки и исследования этих способов явно недостаточно — нужны более точные подходы и методы.

Расчету параметров LC-фильтра инвертора также посвящена обширная литература [2, 3, 8, 9, 11]. Известны методики вычисления параметров как для немодуляционного случая [8, 9, 11], так и для случая ШИМ [2, 3, 8]. Учет пульсаций при ШИМ ведется по-разному. Так, в [1, 2] при расчете токов фильтра предлагается исходить из максимальной пульсации тока дросселя, а результирующие токи определяются прямым (максимальные значения) или геометрическим (эффективные значения) сложением основных и пульсационных составляющих. В [3] уточненный учет пульсаций проведен путем компьютерных расчетов на модели фильтра, а для выбора параметров L и C предлагается исходить из заданного коэффициента несинусоидальности выходного напряжения ξ и условия минимальности суммарной установленной мощности фильтра.

В работах [4, 5, 7, 10] для уточненного расчета процессов в инверторе с ШИМ, в том числе и выходном фильтре, применен аппарат «локального усреднения». Этот подход получил свое развитие в [4–6] в части распространения на разные структуры и схемы инверторов и более детального и точного расчета величин. В данной работе приводится анализ процессов с использованием подхода «локального усреднения» и с дополнительным учетом некоторых реальных параметров элементов.

 

Принятые обозначения

При аналитическом анализе процессов используются величины:

Uн, Iн — эффективные значения напряжения и тока нагрузки;

Udмин, Udном, Udмакс — минимальное, номинальное и максимальное напряжения источника питания;

Zн, φ — сопротивление и угол активно-индуктивной нагрузки;

f = 1/T, fн = 1/h — выходная и несущая частоты (T, h — их периоды);

Ω = f/fн — относительная выходная частота;

x̄L = 2πfL/ZН, C = 2πfCZН — относительные сопротивление дросселя и проводимость конденсатора фильтра на основной частоте;

α, β, χ = α + β — углы между основными гармониками напряжений и токов фильтра (согласно векторной диаграмме рис. 2б);

n — коэффициент трансформации;

М — максимальная глубина модуляции ШИМ-импульсов;

f(t), (t), (t) — функция, ее основная и пульсационная составляющие;

F(t) — переключающая функция.

Эквивалентная схема LC-фильтра и векторные диаграммы

Рис. 2.
а) Эквивалентная схема LC-фильтраж
б) векторные диаграммы

Относительные значения токов и напряжений обозначаются черточкой сверху, амплитудные — индексом «m», максимальные — индексом «макс», a пульсации — индексом «П». В качестве базовых величин для относительных напряжений и токов выбраны эффективные значения напряжения и тока нагрузки. При анализе частично учитываются прямые падения напряжений на вентилях схемы инвертора (транзистора — ΔUVT, диода — ΔUVD). Остальные элементы схем предполагаются идеальными.

 

Подход «локального усреднения»

Применяемый подход к анализу процессов в ШИМ-инверторах основан на двух широко распространенных методах расчета. Его можно трактовать как конкретизацию метода основной гармоники для данного типа преобразователей [7] и как некую «модернизацию» метода переключающих функций [11].

Временные функции в данном случае искажены пульсацией, но имеют либо явно выраженную синусоидальную основу выходной частоты (токи фильтра), либо почти синусоидальные (выходное напряжение). Анализ фильтра выполняется по следующей общей схеме: определяются напряжения и токи элементов основной частоты; интегрально рассчитываются пульсации токов, а по ним — пульсация и коэффициент несинусоидальности выходного напряжения; «наложением» отдельных составляющих определяются полные токи и напряжения компонентов, а по ним — токи дросселя, конденсатора и проч.

На рис. 2а приведена схема фильтра, возбуждаемого эквивалентным источником ШИМ-напряжения v(t), причем в токах и напряжениях элементов доминируют гармоники основной (выходной) частоты. Это позволяет представить временные функции fi(t) в виде суммы основной и пульсационной составляющих:

Формула

Сущность подхода «локального усреднения» заключается в следующем. Если в пределах периода несущей частоты напряжение нагрузки считать постоянным, то можно определить параметры пульсации токов фильтра, имеющей треугольную форму. Эти параметры зависят от коэффициента заполнения ШИМ-импульсов и меняются от такта к такту. При сравнительно больших отношениях несущей и выходной частот, что обычно имеет место на практике, среднее значение ШИМ-импульса упрощенно представляется синусоидальной функцией выходной частоты и амплитудой, пропорциональной глубине модуляции М. Тем самым пульсации тока дросселя (конденсатора) фильтра «приписывается» модуляционная огибающая, и удается записать конечное и достаточно точное выражение в виде функции от М. В итоге получаются замкнутые аналитические выражения для результирующих токов элементов фильтра. Далее расчет ведется оперированием ими.

 

Расчет основных гармоник

При расчете основных гармоник напряжение и ток нагрузки предполагаются синусоидальными. Из векторных диаграмм (рис. 2б) можно найти напряжения и токи элементов фильтра в относительных единицах:

Формула

Временные функции основных гармоник выражаются в виде:

Формула

где V1m — амплитуда основной гармоники ШИМ-импульсов на входе фильтра, а амплитуды UНm, UL1m, IL1m и IC1m определяются умножением соответствующих эффективных значений на коэффициент √2.

 

Расчет пульсационных составляющих при ШИМ-2

При двухуровневой ШИМ (ШИМ-2) на входе фильтра формируется двуполярное импульсное напряжение v(t), принимающее значения ±Vm. Эти импульсы модулированы по ширине, так что напряжение v(t) можно представить произведением Vm на некую переключающую функцию F(t), т. е.:

Формула

Вид последней показан на рис. 3а. Она, в свою очередь, может представляться суммой основной и пульсационной составляющих:

Формула
Диаграммы переключающих функций F(t)

Рис. 3. Диаграммы переключающих функций F(t):
а) при ШИМ-2;
б) при ШИМ-3

Если в k-м такте UН считать неизменным, то пульсация тока дросселя (рис. 4а) принимает треугольную форму с амплитудой

Формула

где γk = αk/h — коэффициент заполнения импульса в k-м такте.

 Временные диаграммы фильтра

Рис. 4. Временные диаграммы фильтра:
а) при ШИМ-2;
б) при ШИМ-3

Наибольшее значение амплитуды достигается при γ = 0,5:

Формула

или же, в относительных единицах,

Формула

«Текущая» амплитуда связана с максимальной зависимостью

Формула

От такта к такту меняется γk, а вместе с ним и амплитуда ILПmk и уровень выходного напряжения. Последний определяется как среднее начение ШИМ-импульса за такт, т. е.

Формула

В следящей системе с ШИМ от такта к такту UHk меняется по синусоиде, и в итоге формируется основная гармоника напряжения, равная H = UHmsinϑ. Основная гармоника входного напряжения фильтра, согласно рис. 2б, опережает его на угол β, т. е.  = F̂Vm.

Связь между Vm и V1m выражается через максимальную глубину модуляции М:

Формула

так что (ϑ) можно приближенно представить функцией

Формула

Последняя позволяет записать конечные выражения для многих временных переменных в функции от М. Заполнению γ можно также присвоить временную функцию:

Формула

Тогда с учетом (8), (9) и (13) для основной составляющей пульсации тока дросселя получим

Формула

Для относительного эффективного значения этой пульсации легко получить

Формула

На практике М обычно берется в пределах 0,85–0,95 и  ≈ (0,3665–0,3972) LПm.

 

Расчет результирующих величин при ШИМ-2

Для эффективного значения полного тока дросселя легко получить

Формула

Пульсация тока дросселя полностью протекает через конденсатор, поэтому

Формула

Для эффективного значения полного тока конденсатора получается

Формула

Как следует из рис. 5, пульсация напряжения конденсатора изменяется параболически, причем амплитуды обеих парабол в k-м такте определяются следующим образом:

Формула

где UСПm — максимальная амплитуда пульсации, имеющая место при γk = 0,5.

 Диаграммы пульсаций тока и напряжения конденсатора фильтра

Рис. 5. Диаграммы пульсаций тока и напряжения конденсатора фильтра

Если определять эффективное значение параболической функции (19) за k-й такт и геометрически вычитать из последнего псевдопостоянную составляющую (определяется интегрированием пульсации), то для эффективного значения всей пульсации получим:

Формула

В относительных единицах это и есть коэффициент несинусоидальности ξ:

Формула

Теоретически М изменяется в пределах 1≥M≥0, и, как следует из (20), ξ меняется в диапазоне ξ = (0,4346–0,7303) UСПm. Практически М = 0,85–0,95 и ξ = (0,4499–0,4866) ŪСПm. Если рассчитать ξ по максимальной пульсации, как обычно поступают на практике, то мы получили бы 0,7071 ŪСПm, т. е. ошибка составляла бы 31,2–36,4%, что существенно. Учитывая важность ξ для расчета параметров фильтра, формулу (20) можно считать хорошим подспорьем в пользу применяемого здесь подхода к анализу.

Коэффициент трансформации n и максимальная глубина модуляции М зависят от схемы коммутатора инвертора. В частности, для мостовой схемы коммутатора имеем:

Формула

где t̄kмакс = tkмакс/h — максимальное относительное время коммутации.

 

Расчет пульсационных составляющих при ШИМ-3

При трехуровневой ШИМ (ШИМ-3) на входе фильтра формируется напряжение с уровнями ±Vm, 0. Выполнив аналогичные шаги для вывода формул (4)–(12), получим следующие выражения:

Формула

В данном случае коэффициент γ есть модуль переключающей функции F(ϑ), т. е.

Формула

Тогда с учетом (22) для основной составляющей пульсации тока дросселя имеем

Формула

Для эффективного значения этой пульсации можно получить

Формула

Получается, что при тех же М = 0,85–0,95 имеем  ≈ (0,3928–0,4365) LПm, т. е. при одинаковых амплитудах пульсации тока ее эффективное значение при ШИМ-3 примерно на 3–5% больше, чем при ШИМ-2.

 

Расчет результирующих величин при ШИМ-3

Для эффективного значения полного тока дросселя получается

Формула

а для эффективного значения полного тока конденсатора фильтра можно записать

Формула

Пульсация напряжения конденсатора при ШИМ-3 также меняется параболически (рис. 5), причем амплитуды парабол здесь в два раза меньше, т. е.

Формула

Выражение для ξ при ШИМ-3 можно получить аналогично вышеописанному, если при интегрировании положить γ = | (ϑ)| = = М|sin(ϑ+β)|. Тогда имеем

Формула

В относительных единицах последнее и есть ξ:

Формула

Как следует из (29), ξ в данном случае теоретически может изменяться в пределах 0≤x≤0,4599 СПm, а на практике (М = 0,85–0,95) имеем ξ = (0,5992–0,6635)СПm. Если учесть, что СПm в два раза меньше, чем при ШИМ-2, то ξ здесь получается примерно в 1,8 раза меньше, чем при ШИМ-2. Иными словами, при ШИМ-3 требуется фильтр с меньшими параметрами, что и является явным преимуществом этого вида модуляции.

Коэффициент трансформации n и М здесь определяются иначе:

Формула

Выражения (20) и (29) могут быть представлены как произведение числового коэффициента, величины Ω2 и двух функционалов:

Формула

где k = π2/4√2 при ШИМ-2 и k = π2/8√2 при ШИМ-3.

Первый функционал в (30) зависит от параметров фильтра и нагрузки

Формула

a второй — от вида и глубины модуляции

Формула

Такое расчленение формул позволяет выбор параметров фильтра сделать очевидным и легко понятным. Должны быть заранее известны параметры Uн, Iн (либо мощность Sн = UнIн), f, Udмин, Udном, Udмакс, пределы изменения коэффициента мощности нагрузки ((cosφ)мин…(cosφ)макс) и допустимый коэффициент ξдоп.

Расчет ведется в следующей последовательности:

  • Исходя из применяемой в инверторе элементной базы, особенностей реализации системы управления, ожидаемых динамических показателей, ограничения коммутационных потерь и др. выбирается вид модуляции (двухуровневая, трехуровневая) и величина несущей частоты fн. Затем определяется Ω: Ω = f/fн.
  • Глубина модуляции М определяется следующим образом. Исходя из временных параметров используемых ключей коммутатора, конфигурации и конкретных параметров снабберных цепей, а также пределов изменения входных и выходных параметров оценивается максимальная длительность переключения ключей коммутатора tк, а
Формула

Следует учитывать, что в идеальном случае М изменяется в пределах 0,5≤М≤1 (ШИМ-2) и 0≤М≤1 (ШИМ-3), и желательно М выбирать максимально возможной.

  • Для выбранного типа ШИМ из (32) определяется значение Ф2.
  • Если определенные в предыдущих пунктах величины подставлять в (31), то получается неравенство
Формула

где А — постоянная.

Фактически неравенство (34) — единственное для выбора двух параметров фильтра. Причем это неравенство необходимо соблюдать при (cosφ)мин (предполагается также, что –x̄LӯС< 1), т. е.

Формула

Второе условие между параметрами фильтра, в зависимости от конкретной задачи, может быть задано различными способами. Рассмотрим некоторые характерные случаи.

  1. В структуре инвертора с выходным низкочастотным трансформатором последний обычно характеризуется значительными индуктивностями рассеивания обмоток, относительное эквивалентное значение которых может составлять
    Формула

    Особым проектированием трансформатора эту величину без серьезных затруднений можно доводить до 0,1–0,15. Следовательно, можно обойтись без дросселя фильтра, поскольку его роль выполняют индуктивности рассеивания обмоток. Фактически при этом в инверторе отсутствует индуктивный элемент, что существенно упрощает конструкцию. Имея значение L = LS, из (35) сразу определяется ӯС, а затем и необходимая емкость:

    Формула

    Во всяком случае, когда в фильтре используется дроссель, его индуктивность определяется разницей Lд = L–LS.

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

    входное сопротивление фильтра чисто активное и, следовательно, реактивная мощность не потребляется [36]. Если при этом коэффициент мощности нагрузки неизменен (например, когда cosφ = 1), то условия (38) могут быть использованы непосредственно. Если же он изменяется, то в (38) можно подставлять характерное значение sinφ. В любом случае, выражение (38) определяет дополнительную связь между параметрами, и совместным решением (35) и (38) можно определять искомые величины L и ӯС. Следует отметить, что при использовании выражения (38) L получается большим, но улучшаются энергетические показатели устройства.

  3. В качестве искомого второго условия выбора параметров фильтра может задаваться суммарная установленная мощность элементов фильтра, относительное значение которой определяется следующим образом [36]:
    Формула

    Если cosφ изменяется, то максимальная мощность фильтра достигается при (cosφ)макс, так что

    Формула

    Если же, исходя из конструктивных требований, предъявляемых к инвертору, потерь мощности в элементах фильтра или иных требований, ограничивать величину Фмакс, то выражение (39) дает дополнительную связь между параметрами L и ӯС, и совместным решением (34) и (39) определятся искомые параметры фильтра.

  4. Иногда определяющим может быть максимальный ток дросселя, который одновременно является током силовых ключей. При ШИМ этот ток, кроме основной составляющей выходной частоты, содержит также пилообразную высокочастотную составляющую (несущей частоты). Ток дросселя косвенно можно характеризовать коэффициентом амплитуды δi следующим образом [3]:
    – ШИМ-2 (40),
    – ШИМ-3 (41),
    где τ, χ — углы; k — коэффициент.

    Формула

Приведенными соотношениями (40) и (41) задается связь между параметрами фильтра, и их решением совместно с (35) можно определять последние.

 

Результаты моделирования

С целью уточненного расчета параметров фильтра исследуемый инвертор моделировался системой MATLAB (подсистемой Simulink). Исследования проводились для широких пределов изменения коэффициента мощности нагрузки и параметров L и C фильтра (L = 0,05–0,25, C =0,1–0,5), а также при Ω = 0,025 и Ω = 0,0125. Результаты исследований при М = 0,9 приведены на рис. 6 и 7, где сплошные кривые соответствуют случаю ШИМ-2, пунктирные — ШИМ-3. Имея пределы изменения cosφ, желательное значение ξ и выбрав L, по допустимой пульсации тока можно вычислять значение C для каждого вида модуляции. Если значение М отличается от 0,9, то кривые можно аппроксимировать в соответствующем направлении.

 Зависимости коэффициента несинусоидальности от параметров фильтра

Рис. 6. Зависимости коэффициента несинусоидальности от параметров фильтра Ω = 0,025 и коэффициентов мощности нагрузки:
а) 0,6;
б) 0,8;
в) 1,0

Зависимости коэффициента несинусоидальности от параметров фильтра

Рис. 7. Зависимости коэффициента несинусоидальности от параметров фильтра Ω = 0,0125 и коэффициентов мощности нагрузки:
а) 0,6;
б) 0,8;
в) 1,0

В качестве примера рассчитаем параметры выходного LC-фильтра инвертора с однократным преобразованием энергии при входном напряжении 60 В (54–75 В), выходном напряжении 220 В, токе нагрузки 4 А, выходной частоте 50 Гц, коэффициенте мощности нагрузки cosφ = 0,8–1,0 и коэффициенте несинусоидальности выходного напражения ξ≤0,03 (3%). Предположим также, что несущая частота выбрана равной 2 кГц, максимальная длительность коммутации 20 мкс, а тип модуляции — ШИМ-2. Расчет ведется в следующей последовательности.

  1. Формула
  2. Формула
  3. Формула
  4. Формула
  5. Учитывая, что инвертор содержит низкочастотный трансформатор, для второго условия между параметрами фильтра исходим из случая 1. Тогда, к примеру, при L = 0,1 из неравенства (35) сразу получается параметр C:
    Формула

    Из соответствующих графиков находим диапазон изменения C, соответствующий заданным пределам изменения коэффициента мощности нагрузки: C ≈ 0,18–0,21. Следовательно, выбираем C для наихудшего случая — C = 0,21.

  6. По найденным значениям L и C определим требуемые значения параметров фильтра:
    Формула

Следует учесть, что расчетные результаты здесь получены для инвертора с идеальными параметрами. На реальной модели стабилизированного инвертора при номинальном питающем напряжении и в указанных пределах изменения cosφ расчитанный фильтр обеспечивает коэффициент несинусоидальности ξ = 3,84–4,78%. Для обеспечения заданного ξ параметры фильтра следует примерно в 1,6 раза увеличивать, например взять L = 0,1 и C = 0,33.

 

Выводы

  1. Рассмотрение параметров выходного сглаживающего LC-фильтра в совокупности с другими основными характеристиками инвертора и с учетом специфических особенностей ШИМ позволяет оптимизировать выбор фильтра.
  2. Применение аппарата «локального усреднения» в сочетании с методами основной гармоники и переключающих функций позволяет исчерпывающе анализировать процессы в выходном фильтре ШИМ-инверторов, получать компактные и в то же время достаточно точные аналитические выражения для расчета и выбора параметров.
  3. Расчет параметров можно непосредственно выполнять посредством приведенных здесь аналитических соотношений, точность которых у большинства величин не ниже 1%, а в ряде случаев — не ниже 2–3%. В более ответственных случаях можно воспользоваться приведенными результатами.
  4. При прочих равных условиях коэффициент несинусоидальности выходного напряжения инвертора в приемлемом на практике диапазоне глубины модуляции (M = 0,85–0,95) при ШИМ-2 примерно в 1,8 раза больше, чем при ШИМ-3.
Литература
  1. Барегамян Г. В. Исследование транзисторного преобразователя с синусоидальным выходным напряжением // Тезисы доклада на Всесоюзной научно-технич. конф. «Применение в технологических процессах машиностроительного производства полупроводн. преобразоват. частоты». Уфа. 1980.
  2. Барегамян Г. В. Выбор параметров LC-фильтра транзисторного инвертора с ШИМ и синусоидальным выходным напряжением // Тезисы доклада на Межведомственном научно-технич. совещании «Проблемы электромагнитной совместимости силовых полупроводниковых преобразователей». Таллин. 1982.
  3. Барегамян Г. В., Маргарян В. В., Симонян Г. Р. Расчет выходного LC-фильтра инвертора с ШИМ // Деп. в АрмНИИНТИ 08.06.98. № 85-Ар98. Ереван. 1998.
  4. Барегамян Г. В. Автоматизация расчета инвертора с трехуровневой ШИМ // Моделирование, оптимизация, управление: Сборник научных трудов / ГИУА. Ереван. 2001. Вып. 4.
  5. Барегамян Г. В. Автоматизация расчета инверторов с синусоидальным выходным напряжением и двухуровневой широтно-импульсной модуляцией // Изв. НАН РА и ГИУА. Серия ТН. 2001. Т. 55. № 3.
  6. Барегамян Г. В. Основы синтеза стабилизированных транзисторных преобразователей с широтно-импульсной модуляцией: Дис. на соискание д-ра техн. наук. Ереван, 2006.
  7. Ковалев Ф. И., Мустафа Г. М. Использование метода первой гармоники при проектировании стабилизированных инверторов // Электротехника. 1982. № 3.
  8. Моин В. С. Стабилизированные транзисторные преобразователи. М.: Энергоатомиздат. 1986.
  9. Мосткова Г. П., Толкачев А. И. Расчет выходных фильтров автономных инверторов // Электротехн. промышленность. Сер. Преобр. техника. 1971. № 22.
  10. Мустафа Г. М. Расчет импульсно-модулированных инверторов с синусоидальным выходным напряжением // Электротехника. 1988. № 2.
  11. Руденко В. С., Сенько В. И., Чиженко И. М. Основы преобразовательной техники. М.: Высшая школа. 1980.
  12. Адамия Г. Г., Беркович Е. И., Картавых А. С. и др. Статические агрегаты бесперебойного питания. М.: Энергоатомиздат. 1992.

Добавить комментарий

Ваш адрес email не будет опубликован. Обязательные поля помечены *