http://pspice.narod.ru

Трансформаторная схема

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

В качестве предмета изучения мы используем трансформаторнную схему (рис. 1) в которой трансформатор имеет одинаковые обмотки а коэффициент трансформации 0.9.
Сопротивление обмоток представлено резисторами R1, R2.

Рис. 1. Трансформаторная схема

[Maple Bitmap]

Составим и решим уравнения Кизгофа для этой схемы.

>    restart: with(MSpice): ESolve(Q,`T1.CKT`):

`Cистема Кирхгофа-Лапласа`

(VOUT-V3)/R2-(V3-V0_1)/s/(L1[1]-K1*sqrt(L1[1]*L2[1])) = 0

(V3-V0_1)/s/(L1[1]-K1*sqrt(L1[1]*L2[1]))+(V2-V0_1)/s/(L2[1]-K1*sqrt(L1[1]*L2[1]))-V0_1/s/K1/sqrt(L1[1]*L2[1]) = 0

-(V2-V1_SIN)/R1-(V2-V0_1)/s/(L2[1]-K1*sqrt(L1[1]*L2[1])) = 0

-VOUT/RL-(VOUT-V3)/R2 = 0

`Решения`

{VOUT, V0_1, V3, V2}

``

MSpice v8.44:   http://pspicelib.narod.ru

Заданы узлы: {VINP} Источники: [V1_SIN]

Решения V_NET: [VOUT, V0_1, V3, V2]

J_NET: [JR1, JR2, JRL, JK1, JV1_SIN, JL1[1], JL2[1]]

``

Найдём операторный коэффициент передачи

>    H:=simplify(VOUT/VINP,'size');

H := -(L1[1]*L2[1])^(1/2)*K1*s*RL/(L1[1]*L2[1]*(K1-1)*(K1+1)*s^2+((-RL-R2)*L2[1]-R1*L1[1])*s-R1*(R2+RL))

Найдём частотный коэффициент передачи

>    K:=simplify(AV(H,f),'size');

K := 2*I*(L1[1]*L2[1])^(1/2)*K1*Pi*f*RL/(4*L1[1]*L2[1]*(K1-1)*(K1+1)*Pi^2*f^2+2*I*((R2+RL)*L2[1]+R1*L1[1])*f*Pi+R1*(R2+RL))

Найдем напряжение на выходе, подав на вход синусоидальный сигнал частотой 50 Гц, амплитудой 1В. Будем использовать преобразование Лапласа.

>    Values(AC,PRLCVI,["RL"]): Digits:=4: H:=simplify(H,'size'):  

>    with(inttrans):
V1_SIN:=sin(2*Pi*f*t);
VOUT:=simplify(invlaplace(laplace(V1_SIN,t,s)*H,s,t),'size');

V1_SIN := sin(2*Pi*f*t)

VOUT := 225./(.2468e7*RL^2*f^2+.8934e8*RL*f^2+.8934e9*f^2+.6250e9*RL^2+.1250e11*RL+.6250e11+.3518e5*f^4)*(15.71*(RL+20.)*RL*f/(25.*RL^2+810.*RL+8100.)^(1/2)*(.2500e5*RL+.2500e6+187.6*f^2)*(-1.*exp(5.26...
VOUT := 225./(.2468e7*RL^2*f^2+.8934e8*RL*f^2+.8934e9*f^2+.6250e9*RL^2+.1250e11*RL+.6250e11+.3518e5*f^4)*(15.71*(RL+20.)*RL*f/(25.*RL^2+810.*RL+8100.)^(1/2)*(.2500e5*RL+.2500e6+187.6*f^2)*(-1.*exp(5.26...
VOUT := 225./(.2468e7*RL^2*f^2+.8934e8*RL*f^2+.8934e9*f^2+.6250e9*RL^2+.1250e11*RL+.6250e11+.3518e5*f^4)*(15.71*(RL+20.)*RL*f/(25.*RL^2+810.*RL+8100.)^(1/2)*(.2500e5*RL+.2500e6+187.6*f^2)*(-1.*exp(5.26...
VOUT := 225./(.2468e7*RL^2*f^2+.8934e8*RL*f^2+.8934e9*f^2+.6250e9*RL^2+.1250e11*RL+.6250e11+.3518e5*f^4)*(15.71*(RL+20.)*RL*f/(25.*RL^2+810.*RL+8100.)^(1/2)*(.2500e5*RL+.2500e6+187.6*f^2)*(-1.*exp(5.26...
VOUT := 225./(.2468e7*RL^2*f^2+.8934e8*RL*f^2+.8934e9*f^2+.6250e9*RL^2+.1250e11*RL+.6250e11+.3518e5*f^4)*(15.71*(RL+20.)*RL*f/(25.*RL^2+810.*RL+8100.)^(1/2)*(.2500e5*RL+.2500e6+187.6*f^2)*(-1.*exp(5.26...

Построим график напряжения на нагрузки. При этом переведём время в Pi-ках,  выполнив замену переменных  t=t[Pi]/(2*f).

>    VOUT:=subs(t=t[Pi]/(2*f),VOUT):

>    f:=50:
subs(t=t[Pi]/(2*f),V1_SIN),eval(VOUT,RL=10),eval(VOUT,RL=500),0:
ploth([%],t[Pi]=0..2,"2) Напряжение на нагрузке f=50 Гц del[V1_SIN,VOUT(RL=10),VOUT(RL=50),V0]");

[Maple Plot]

Интересно посмотреть как ведет себя коэффициент трансформации и фаза.

>    psi:=AVPdeg(H,f); absK:=simplify(AVM(H,f)): abs('K')=absK;
ploth([absK],RL=0..200,"2) Коэффициент трансформации del[absK]");
ploth([psi],RL=0..200,"3) Фазовая нагрузочная характеристика [Фаза(град)]");

psi := 180/Pi*arctan(1.000*RL*(-19.00*Pi^2+100+10*RL),1.000*RL*(10.0*Pi*RL+200.0*Pi))

abs(K) = 25.29*(RL^2*(.5030e6+.4717e5*RL+1359.*RL^2)/(.4026e6+.3774e5*RL+1087.*RL^2)^2)^(1/2)

[Maple Plot]

[Maple Plot]

Hosted by uCoz