Coupled LC circuits - contradictory equations

I consider two identical LC circuits, inductively coupled with the coupling coefficient k=1. Let y1 be the charge on the capacitor of the first circuit, y2 - charge on the capacitor of the second circuit. Let y1(0) = q > 0, y2(0) = 0. The equations that I get are y1'' + y2'' + ω^2 * y1 =0, y1'' + y2'' + ω^2 * y2 =0, where ω^2 = 1/(LC).
I subtract the second equation from the first and I obtain ω^2 * (y1-y2)=0, which contradicts the initial conditions y1(0)=q>0, y2(0)=0.
I do not understand why this problem appears and what to do.

k=1 means both inductances respectively both capacitors are connected in parallel with zero impedance inbetween. It's physically
impossible to have different voltage at the capacitors in this case, your assumed initial conditions can't exist, all calculations based on
it are meaningless.

LC circuits are INDUCTIVELY coupled, so my initial conditions are clearly possible.

I don't think so. In my opinion, it's similar to a problem previously discussed at edaboard: shorting two capacitors with a zero impedance switch.

It's the "same" behavior when you connect two caps in parallel with differents voltages.


Solve with k<>1 and see what happen when k->1.
The solution is the sum of two sinewaves, one of w = 1/sqrt(LC) and other of w->infinity.

In real cases, the second sinewaves vanish quickly because you always have losses.

No. You have inductance which prevents appearance of any much higher frequency components. It may be similar to connecting two capacitors by means of an inductive connector.

Perhaps you "see" the obvious facts from an equivalent circuit:

Please, just solve with k<>1 and see what happen when k is close to 1

Perhaps you are right. The inductance gets cancelled, because the secondary coil is initially shorted by the uncharged capacitor. So what happens, say, for k=.99?

Including the coupling coefficient k in the equations

y1'' + k*y2'' + ω^2 * y1 =0
k*y1'' + y2'' + ω^2 * y2 =0


Adding and substracting both equations

(1+k)*(y1+y2)'' + ω^2 *(y1+y2) = (1+k)*U1'' + ω^2 *U1 = 0
(1-k)*(y1-y2)'' + ω^2 *(y1-y2) = (1-k)*U2'' + ω^2 *U2 = 0


Then, the solutions are trivials.

U1 = y1+y2 : general sinewave of frequency w/sqrt(1+k)
U2 = y1-y2 : general sinewave of frequency w/sqrt(1-k)

--> y1 = (U1+U2)/2 and y2 = (U1-U2)/2

Is easy to see what when k-->1 or k-->-1 , one natural frequency -->w/sqrt(2) and the other goes to infinity.


In real circuits, the high frequency vanish as k goes close to +/-1 due the losses.

The image is a simulation with 0.1 ohm of windings resistance for k = 0.995 and k=0.9999



At k=+/-1 with no losses you have a singularity: Infinite current at infinite frequency.

Hi Eduardo,

I have a different opinion on the contradiction showed:

In fact you have analyzed the circuit in frequency domain, These one approach is usable only in stady state condition (sinusoidal in this case) then without initial condition.

Initial condition are usable only in transient analisys.

For me the contradiction don't exist.

Ninux

In fact you have analyzed the circuit in frequency domain

Click to expand...

He didn't. Look sharp, it's the differential equation solution in time domain.

Hi experts,

Has Anybody done some work on coupled resonators modelling in MATLAB. Any example would be might helpful for me. Thanks!!! Waiting for your reply.