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

Реферат

РАСЧЕТ РАБОЧЕГО КОНТУРА С ИСПОЛЬЗОВАНИЕМ ВЫЧИСЛИТЕЛЬНОЙ ТЕХНИКИ

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

Расчет рабочего контура сводится к составлению системы линейных алгебраических уравнений энергетических балансов и решению этой системы уравнений. Однако из-за громоздкости большинства уравнений энергобалансов при решении системы возможны технические ошибки, что может значительно сказаться на точности расчета.

Для решения систем линейных алгебраических уравнений разработаны различные математические методы, для которых составлены и внедрены в математическое обеспечение вычислительной техники достаточно удобные в использовании стандартные программы. Если имеется возможность использовать вычислительную технику с такими программами, то процесс определения расходов сред в ветвях рабочего контура можно существенно упростить. При этом можно значительно расширить количество неизвестных определяемых из системы уравнений, включив в них ряд дополнительных величин. На примере варианта принятого в расчет энергоблока можно показать подготовку системы уравнений к расчету с использованием компьютера. Для расчета системы уравнений необходимо выполнить все подготовительные операции, изложенные в 14 разделе. К таким операциям следует отнести:

• составление перечня неизвестных подлежащих расчету;
• определение значений параметров теплообменивающихся сред. В нашем случае это значения энтальпий, приведенные в таблице (для системы теплофикации) и в таблицах 14.3, 14.4, 14.5 (для рабочего контура);
• составление системы уравнений энергетических балансов в численном виде. Для системы теплофикации уравнения в численном виде приведены в подразделе 14.2, для рабочего контура - в подразделах 14.6, 14.7, 14.8, 14.9.

2. Расчет системы теплофикации

Перечень неизвестных величин подлежащих определению:св - расход сетевой воды, кг/с; гпсвп - расход греющего пара на пиковый подогреватель сетевой воды, кг/с; гпсво2ст - расход греющего пара на основной подогреватель сетевой воды 2-й ступени, кг/с;гпсво1ст - расход греющего пара на основной подогреватель сетевой воды 1-й ступени, кг/с; гпсв - суммарный расход греющего пара на систему теплофикации, кг/с.

Система уравнений тепловых балансов системы теплофикации:

уравнение теплового баланса подогревателя сетевой воды 2-й ступени

,996[Gгпсво2ст (2834,82 - 530,76) + Gгпсвп(659,98 - 530,76)] = Gсв(505,3-378,8);

уравнение теплового баланса подогревателя сетевой воды 1-й ступени

,997[Gгпсво1ст (2710,88-427,71)+( Gгпсвп + Gгпсво2ст)?(530,76-427,71)]= Gсв·(378,8-295);

уравнение суммарного расхода греющего пара

гпсв = Gгпсвп + Gгпсво2ст + Gгпсво1ст.

Результаты расчета системы теплофикации с помощью компьютера и сопоставление их с результатами «ручного» расчета (с помощью калькулятора) представлены в таблице 1. При этом для простоты сравнения результаты компьютерного расчета приведены с точностью до восьми значащих цифр.

Таблица 1 - Результаты расчета системы теплофикации

Расчетная величинаРезультаты «ручного» расчета, АРРезультаты компьютерного расчета, АКПогрешность |(АР - АК)/АК|×100, %Gсв, кг/с314,37314,370500,00016Gгпсвп, кг/с22,5322,5347030,021Gгпсво2ст, кг/с16,0716,0653990,029Gгпсво1ст, кг/с9,839,83097110,010Gгпсв, кг/с48,4348,4310740,002

Сопоставление результатов расчета показывает, что «ручной» расчет, приведенный в примере, выполнен корректно, точность расчета высокая.

. Подготовка системы уравнений рабочего контура. Перечень и обозначения расчетных величин

) G1стЦВД - расход пара через 1 ступень ЦВД, кг/с;

) G2стЦВД - расход пара через 2 ступень ЦВД, кг/с;

) G3стЦВД - расход пара через 3 ступень ЦВД, кг/с;

) G4стЦВД - расход пара через 4 ступень ЦВД, кг/с;

) G5стЦВД - расход пара через 5 ступень ЦВД, кг/с;

) G1стЦНД - расход пара через 1 ступень ЦНД, кг/с;

) G2стЦНД - расход пара через 2 ступень ЦНД, кг/с;

) G3стЦНД - расход пара через 3 ступень ЦНД, кг/с;

) G4стЦНД - расход пара через 4 ступень ЦНД, кг/с;

) G5стЦНД - расход пара через 5 ступень ЦНД, кг/с;

) Gс - расход сепарата из сепаратора, кг/с;

) Gпр - потеря пара в протечках, кг/с;

) Gпп - расход греющего пара на ПП, кг/с;

) Gот1 - расход пара в 1 отборе, кг/с;

) Gот2 - расход пара в 2 отборе, кг/с;

) Gот3 - расход пара в 3 отборе, кг/с;

) Gот4 - расход пара в 4 отборе, кг/с;

) Gот5 - расход пара в 5 отборе, кг/с;

) Gот6 - расход пара в 6 отборе, кг/с;

) Gот7 - расход пара в 7 отборе, кг/с;

) Gот8 - расход пара в 8 отборе, кг/с;

) Gптпн - расход пара на ТПН, кг/с;

) Gт - расход пара на входе в ЦВД, кг/с;

) Gпг - паропроизводительность ПГ, кг/с;

) i пв - энтальпия ПВ на входе в ПГ, кДж/кг;

) Qяр - тепловая мощность ЯР при работе ЯЭУ на номинальной мощности, МВт;

) hЯЭУбрутто - КПД ЯЭУ брутто.

. Система уравнений рабочего контура

Об уравнениях (1) - (11) см. п.14.8.4.

) G1стЦВД = Gт;

) G2стЦВД = G1стЦВД;

) G3стЦВД = G2стЦВД - Gот1;

) G4стЦВД = G3стЦВД - Gот2;

) G5стЦВД = G4стЦВД - Gот3;

) G1стЦНД = G5стЦВД - Gот4 - Gс - Gптпн;

) G2стЦНД = G1стЦНД - Gот5;

) G3стЦНД = G2стЦНД - Gот6;

) G4стЦНД = G3стЦНД - Gот7;

) G5стЦНД = G4стЦНД - Gот8;

) Gс = (G5стЦВД - Gот4) (1 - 0,856);

) Gпр = 0,004 Gт (п.14.8.8);

) 0,99Gпп(2778,8-1207,16) = (0,8199 Gт - 63,9875) 0,856 (2968,29 - 2755,0);

) 0,992 Gот1(2629,24 - 843,68) = (Gт + Gпр) (897,7 - 815,4);

) 0,993[Gот2 (2576,64 - 728,22) + Gот1 (843,68 - 728,22)] = (Gт + Gпр) (815,5 - 702,0);

16) 0,994 {(Gот3 - 41) 2520,99 + 14 ? 546,73 + (Gот1 + Gот2) ? 728,22 + (Gот4 - 22,53) ? 659,98 + (Gт -Gот i) (1 - 0,856) 664,7 + [(Gт -Gот i)? 0,856 + 27 + Gпр + 22,53] 645,51} = (Gт + Gпр) ? 697,1;

) 0,995 (Gот4 - 22,53) ? (2456,02 - 659,98) = [(Gт - Gот i) ? 0,856 + 27 + Gпр + 22,53] (645,51 - 517,16);

18) 0,996(Gот5 - 16,07) (2834,82 - 440,27) = [(Gт - Gот i) 0,856 + 27 + Gпр - 16,07 - 9,83] (517,18- 411,54) + 48,43×(517,18 - 427,71);

) 0,997[ (Gот6 - 9,83) (2710,88 - 349,49) + (Gот5- 16,07)×(440,27 - 349,49)] = [(Gт - Gот i) 0,856 + 27 + Gпр - 9,83 - 16,07]×(411,57 - 321,2);

) 0,998 {Gот7 2586,11 + (Gот5 + Gот6 - 16,07- 9,83)×349,49 + [(Gт - Gот i) 0,856 + 27 + Gпр - Gот i ]×232,67} = [(Gт - Gот i) 0,856 + 27 + Gпр - 16,07 - 9,83]×320,0;

) 0,999{Gот8 2473,33 + [(Gт -Gот i ) 0,856 + 27 + Gпр - Gот i] 150,757} = [(Gт -Gот i) 0,856 + 27 + Gпр - Gот i] 232,67;

22) расход пара на ТПН (п.14.7)

Мощность КГТН

кгтн = pкгтн×Gкгтн×vкгтн / (hкгтн×1000) = 2352×Gпп ×0,0013148 / 0,8.

Расход силовой воды на гидропривод

гп = 1000×N кгтн /(vгп×Dргп×hгп) = 2,352×Gпп ×0,0013148 / (0,8×0,0011015×8,4×0,8);

Массовая подача ТПН

тпн = Gт + Gпр + Gгп = Gт + Gпр + 2,352×Gпп ×0,0013148 / (0,8×0,0011015×8,4×0,8).

Мощность ТПН

тпн = pтпн×Gтпн×vтпн / (hтпн×1000) = 8700×[ Gт + Gпр + 2,352×Gпп ×0,0013148 / (0,8×0,0011015×8,4×0,8)] ´ 0,0011015/0,8.

Расчетное выражение для определения расхода пара на ТПН

птпн = Nтпн /(hмех.тпн×Н i тпн) = 8700×[ Gт + Gпр + 2,352×Gпп ×0,0013148 / (0,8×0,0011015×8,4×0,8)] ´ 0,0011015/(0,8 0,98×573,39);

) уравнение баланса механической энергии комплекса «главная турбина - генератор электрической энергии».

Внутренняя мощность турбины i т = 1041,2327×103 кВт (п.14.8.2).

Внутренние теплоперепады, срабатываемые на ступенях турбины h i ст

Внутренняя мощность ступени турбины

N i ст=G i ст ×h i ст (п.14.8.5).

Расчетное уравнение баланса механической энергии принимает вид:

,58×G1стЦВД + 73,98×G2стЦВД + 52,6×G3стЦВД + 55,65×G4стЦВД + 64,97×G5стЦВД + 133,47 ×G1стЦНД + 123,94 ×G2стЦНД + 124,77 ×G3стЦНД + 112,78 ×G4стЦНД + 139,93 ×G5стЦНД = 1041,2327×103;

) паропроизводительность ПГ(п.14.8.8)

пг = Gт + Gпр + Gпп ;

) энтальпия ПВ на входе в ПГ (п.14.9.)

пв = [(Gт + Gпр ) × 897,7 + Gпп ×1207,16] / Gпг ;

) тепловая мощность ЯР (п.14.9)

Qяр = G пг (2778,8 - i пв) / (0,985 × 103) ;

) КПД ЯЭУ брутто (п.14.9)

hЯЭУбрутто =1000 / Qяр .

5. Результаты расчета. Анализ полученных результатов

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

Таблица 2 - Результаты расчета рабочего контура

Расчетная величинаРезультаты «ручного» расчета, АРРезультаты компьютерного расчета, АКПогрешность |(АР - АК)/АК|×100, %1. G1стЦВД, кг/с1548,9091548,95200,0032. G2стЦВД, кг/с1548,9091548,95200,0033. G3стЦВД, кг/с1476,65241476,69410,0034. G4стЦВД, кг/с1385,03441385,04250,0015. G5стЦВД, кг/с1306,64571306,64590,0006. G1стЦНД, кг/с997,5926997,583890,0017. G2стЦНД, кг/с933,6643933,648700,0028. G3стЦНД, кг/с885,7718885,753500,0029. G4стЦНД, кг/с851,1723851,158530,00210. G5стЦНД, кг/с818,666818,652000,00111. Gс, кг/с173,6608173,659830,00112. Gпр, кг/с6,195646,19580800,00313. Gпп, кг/с141,5105141,511870,00114. Gот1, кг/с72,256672,2579100,002 Gот2, кг/с91,61891,6516040,03716. Gот3, кг/с78,388978,3966010,00117. Gот4, кг/с100,667100,674870,00818. Gот5, кг/с63,928363,9351860,00119. Gот6, кг/с47,892647,8952000,00520. Gот7, кг/с34,599534,5949660,01321. Gот8, кг/с32,506432,5065320,00022. G птпн, кг/с34,726634,7273380,00223. Gт, кг/с1548,9091548,95200,00324. Gпг, кг/с1696,61511696,65970,00325. i пв, кДж/кг923,5113923,510870,00026. Qяр, МВт3195,6453195,73030,00327. hЯЭУбрутто0,31290,312917520,006

Сопоставление результатов расчета показывает, что «ручной» расчет выполнен достаточно корректно, результаты расчета имеют высокую точность.

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

Ниже приведены примеры результатов некоторых исследований.

) Основной вариант установки принят с охладителями дренажа в ПНД-3 и ПНД-4. В исследуемом варианте принята аналогичная установка, но без охладителей дренажа в указанных ПНД. Для расчета показателей такой установки достаточно в исходную систему уравнений внести следующие изменения:

в уравнении теплового баланса ПНД-4 (уравнение 18) значение энтальпии греющей среды на выходе из ПНД-4 вместо i=440,27 кДж/кг (переохлажденный дренаж) следует принять i=530,76 кДж/кг (энтальпия воды на линии насыщения при давлении в ПНД-4 р=0,242 МПа);

в уравнении теплового баланса ПНД-3 (уравнение 19) значение энтальпии греющей среды на выходе из ПНД-3 вместо i=349,49 кДж/кг (переохлажденный дренаж) следует принять i=427,707 кДж/кг (энтальпия воды на линии насыщения при давлении в ПНД-3 р=0,109 МПа). Кроме того, в уравнении 19 значение энтальпии греющей среды на входе - 2, куда каскадно сливается дренаж из ПНД-4, вместо iвх2г=440,27 кДж/кг следует принять iвх2г=530,76 кДж/кг;

в уравнении теплового баланса ПНД-2 (уравнение 20) значение энтальпии среды на входе 2, куда каскадно сливается дренаж из ПНД-3, вместо iвх2г=349,49 кДж/кг следует принять iвх2г=427,707 кДж/кг.

Параметры нагреваемой среды приняты такими же как и в основном варианте.

Как и следовало ожидать, исключение охладителей дренажа в ПНД-3 и ПНД-4 привело к некоторому снижению экономичности установки. Если в основном варианте расчета величина КПД ЯЭУ брутто составляет 0,31291752, то исключение упомянутых охладителей дренажа снизило КПД до 0,31276585. Как следствие исключения охладителей дренажа отмечается некоторое повышение расхода греющего пара на ПНД-3 и ПНД-4: для ПНД-3 Gот6 = 48,876965 кг/с (вместо 47,895199 кг/с в основном варианте); для ПНД-4 Gот5 = 65,839365 кг/с (вместо 63,935186 кг/с в основном варианте). Заметим, что исключение охладителей дренажа в ПНД-3 и ПНД-4 снижает КПД установки весьма незначительно, хотя заметно упрощает ее схему. Можно предположить, что это и побудило проектанта не применять в прототипной установки эти ОД.

) Известно, что при работе установки могут возникнуть какие-либо неисправности в подогревателях высокого давления. В схеме рабочего контура предусмотрена арматура отключения неисправных подогревателей и байпасная линия обвода отключенных ПВД. Для расчета такого режима достаточно в уравнениях тепловых балансов ПВД-6 и ПВД-7 (уравнения 14, 15) а также точки смешения сред на выходе из ПВД-7 (уравнение 25) значения энтальпии нагреваемой среды принять одинаковыми и равными энтальпии рабочего тела на выходе из питательного насоса i = 702 кДж/кг

Результаты расчета рабочего контура в этом режиме показали, что экономичность установки заметно снизилась, КПД ЯЭУ брутто составил 0,30845731 (вместо 0,31291752 в основном варианте). Разумеется расчетные значения отборов пара на отключенные ПВД равны нулю (Gот1 =Gот2 =0). Резко снизилось значение энтальпии на входе в ПГ. Ее значение составило iпв=748,64176 кДж/кг (вместо 923,51087 кДж/кг в основном варианте). Соответственно изменилась и температура питательной воды на входе в ПГ. Она составила tпв = t(pпв =8,12 МПа; iпв=748,64 кДж/кг) = 175,85оС, вместо tпв = 215,19 оС в основном варианте. 3) Рассмотрим вариант отключения ПНД-3. Условно считаем, что у этого ПНД по нагреваемой среде предусмотрена байпасирующая линия. По греющей среде каскадный слив из ПНД-4 переведен на ПНД-2. В результате уменьшилось количество подогревателей в системе регенерации. Однако, последний ПВД (ПВД-7) остался подключенном по отбираемому пару на прежнюю точку проточной части турбины. Следовательно, температура нагреваемой воды будет несколько завышена по сравнению с оптимальной. Кроме того, заметно нарушилась равномерность нагрева воды вдоль тракта конденсатно-питательной системы. Это должно привести к некоторому дополнительному снижению экономичности установки.

Для реализации этого режима в систему уравнениц внесены следующие изменения:

в уравнении теплового баланса ПНД-3 (уравнение 19) энтальпию нагреваемой среды принимают одинаковой и равной энтальпии воды на выходе из КН-2, т. е. iнвх ПНД-3 = iнвых ПНД-3 = 321,2 кДж/кг. Кроме того, расход по каскадному сливу из ПНД-4 принят равным нулю, т.е. изъято слагаемое (Gот5 - 16,7) (440,27 - 349,49);

в уравнении теплового баланса ПНД-4 (уравнение 18) энтальпия нагреваемой среды на входе - 1 принята равной энтальпии в ПНД-3, т. е. iнвх1 ПНД-4 =321,2 кДж/кг (вместо 411,54 кДж/кг в основном варианте);

в уравнении теплового баланса ПНД-2 (уравнение 20) на входе - 2, куда каскадно сливается выход греющей среды из ПНД-4 и ПНД-3, энтальпия среды принята iвх2 ПНД-2 = 440,27 кДж/кг (по значению i гвых ПНД-4 - см. таблицу 14.4) вместо 349,49 кДж/кг в основном режиме.

Анализ полученных результатов показал, что несколько снизилась экономичность установки: КПД установки составил 0,31161317 (вместо 0,31291752 в основном режиме). Кроме того, значительно вырос расход пара на ПНД-4: Gот5 = 103,38895 (вместо 63,935186 в основном режиме). 4) Рассмотрим вариант исследования вклада системы регенерации в КПД АЭС в целом. Для этого решена система уравнений тепловых балансов расчетной схемы с отключением системы регенерации. Результаты решения сопоставлены с результатами основного варианта, в котором все элементы рабочего контура, в том числе и система регенерации, работают по штатной схеме.

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

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

Так как температура питательной воды при отключенной системе регенерации значительно ниже температура воды в штатном режиме, то для упрощения рассмотрения вопроса в уравнениях тепловых балансов всех ВП коэффициенты удержания тепла приняты равными единице. В связи с этим результаты решения задачи сопоставляются не с основным штатным вариантом, а с основным вариантом, в котором система регенерации включена в работу, но коэффициенты удержания тепла также приняты равными единице. Отключение системы регененрации должно привести к тому, что расходы греющего пара должны составить: Gот1=0; Gот2=0; Gот3=Gсн=41; Gот4=Gпсвп=22,53; G от5 =Gпсво2ст=16,07; Gот6=Gпсво1ст=9,83; Gот7=0; Gот8=0. В связи с этим слагаемые уравнений тепловых балансов, в которых значения нулевых расходов умножаются на штатные большие значения энтальпий, можно оставить без изменений. По остальным слагаемым в уравнения тепловым балансам внесены следующие изменения:

а) ур.21 (ПНД1) - iвых=150,757 вместо 232,67;

б) ур. 20 (ПНД2) - iвх3=150,757 вместо 232,67; iвых=150,757 вместо 320,0;

в) ур. 19 (ПНД3) - iнвх=150,757 вместо 321,2; iнвых=150,757 вместо 411,57;

г) ур. 18 (ПНД4).

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

.

рабочий контур теплофикация энергетический

В основное уравнение (18) внесены следующие изменения:

в первом слагаемом правой части iнвых=i1 вместо 517,18; iнвх1=150,757 вместо 411,54;

во втором слагаемом правой части iнвых =i1 вместо 517,18;

д) ур. 17 (ПНД5) - iнвх =i1 вместо 517,16; iнвых =i1 вместо 645,51;

е) ур. 16 (деаэратор).

Для определения энтальпии среды на выходе из деаэратора (обозначена i2) составлено дополнительное уравнение теплового баланса смешения сред, поступающих на деаэратор. С учетом ожидающихся расходов сред в отборах пара из турбины расходы сред на входах в деаэратор составляют:вх1=0; Gвх2=Gсн1=14; Gвх3=0; Gвх4=0; Gвх5=Gс=(Gт-41-22,53)×(1-0,856); Gвх6=(Gт-41-22,53)×0,856+27+Gпр+22,53; Gвых=Gт+Gпр.

Значения энтальпий смешиваемых сред:

iвх2=546,73; iвх5=664,7; iвх6=i1; iвых=i2.

{546,73×14+664,7(Gт -41-22,53)(1- -0,856)+ i1[(Gт - 41-22,53)0,856+27+Gпр+22,53]}/( Gт+Gпр).

В основное уравнение (16) внесены следующие изменения:

в правой части уравнения iвых=i2 вместо 697,1;

в шестом слагаемом левой части уравнения iвх6=i1 вместо 645,51;

ж) ур. 15 (ПВД6).

Дополнительные уравнения: энтальпия нагреваемой среды после прохождения ею ПН i3 = i2 + Diпн. Для расчета приближенно оценено Diпн = 5 кДж/кг. В основное уравнение (15) внесены следующие изменения. iнвх=i3 вместо 702,0; iнвых =i3 вместо 815,5;

з) ур. 14 (ПВД7). iнвх =i3 вместо 815,5; iнвых =i3 вместо 897,7;

и) ур. 25 (точка смешения нагреваемой среды на выходе из ПВД7 и слива конденсата греющего пара из СПП): iнвых ПВД7=i3 вместо 897,7.

Анализ результатов решения:

) подтвердились предположения о значениях расходов в отборах;

) значение энтальпии ПВ на входе в ПГ значительно снизилось и составило 337,19 кДж/кг вместо 923,6 в основном варианте. Это соответствует значительно сниженной величине температуры питательной воды - 79,02оС вместо 215,22 в основном варианте;

) КПД ЯЭУ брутто снизился и составил 0,2826 вместо 0,3137 в основном варианте, т.е. КПД снизился на 3,11%. Как уже отмечалось, умеренное снижение КПД можно объяснить тем, что в его поддержание вносят некоторый вклад сливаемые в поток ПВ горячие сливы от остающихся в штатном режиме работы элементов рабочего контура (СПП, система теплофикации, СН);

) мощность ядерного реактора возросла и составила 3538,72 МВт вместо 3187,8;

) расход пара на турбину снизился и составил 1279,76 кг/с вместо 1544,69 в основном варианте.

Литература

Рожкова Л.Д., Карнеева Л.К., Чиркова Т.В. Электрооборудование электрических станций и подстанций, 2004

Файбисович Д.Л. и др. Справочник по проектированию электрических сетей, 2006

Неклепаев, Б. Н. Электрическая часть электростанций и подстанций, 1989

Бурман А.П., Строев В.А., Основы современной энергетики. Часть 2. Современная электроэнергетика, 2003

Электротехнический справочник. Том 3. Производство, передача и распределение электрической энергии

Ополева Г.Н. Схемы и подстанции электроснабжения, 2006

Крючков И.П., Неклепаев Б.Н., Старшинов В.А. Расчет коротких замыканий и выбор электрооборудования, 2005

Коновалова Л.Л., Рожкова Л.Д. Электроснабжение промышленных предприятий и установок, 2009

Рожкова Л.Д., Козулин В.С. Электрооборудование станций и подстанций, 2007