15
\$\begingroup\$

I'm trying to do find the voltage over R1 in the following circuit, where L1/L2 is an ideal transformer. LTSpice complains that the "Matrix is singular". Why? I've tried to play around with lots of different values in order to see if it's a problem with approximation.

enter image description here

The numbers after "AC" are the max amplitude and phase (in degrees).

\$\endgroup\$
4
  • \$\begingroup\$ The analysis works fine if I remove R1 (replacing it with a gap). \$\endgroup\$ – Karin Oct 12 '11 at 13:33
  • \$\begingroup\$ Can you make L1 and L2 non ideal by a very small amount. Maybe add a milliohm of less. I've had something similar happen a few decades ago :-) \$\endgroup\$ – Russell McMahon Oct 12 '11 at 13:51
  • 1
    \$\begingroup\$ @RussellMcMahon: I tried to set their "parallell resistances" to 1m, if that's what you meant. It resulted in a current over R1 on the order 10^-10A. \$\endgroup\$ – Karin Oct 12 '11 at 14:01
  • \$\begingroup\$ Is it working? Did the 10M between ccts make it work? If you emove the 10M does it work? If you add parallel R does it wiorK. What current are you getting / expecting? \$\endgroup\$ – Russell McMahon Oct 12 '11 at 15:59
18
\$\begingroup\$

You need a DC path between the two circuits. Put a high value resistor between them, say 10M.

I checked that it worked using Pulsonix (actually SIMetrix) SPICE. I got a singular matrix error without the resistor.

\$\endgroup\$
4
  • 9
    \$\begingroup\$ or just ground the secondary \$\endgroup\$ – endolith Oct 12 '11 at 14:49
  • \$\begingroup\$ @endolith: How could grounding both circuits work? There is a difference in potential between all nodes in the circuit as it is designed now. \$\endgroup\$ – Karin Oct 12 '11 at 15:07
  • 1
    \$\begingroup\$ I considered that, but realised that Karin probably wanted the secondary to be isolated. \$\endgroup\$ – Leon Heller Oct 12 '11 at 15:55
  • 8
    \$\begingroup\$ @Karin - the secondary has no reference to anything at all as it is now. As shown, in real; life it could be at -10V average relative to ground, or at + 1,000,000 V relative to ground. The maths may have trouble with this lack of definition. The 10Mohm or the ground boty do the sme job of providing a finite connection between the two. \$\endgroup\$ – Russell McMahon Oct 12 '11 at 16:05
5
\$\begingroup\$

There is a SPICE parameter called RSHUNT which adds shunt resistors to GND on every node. By default it is usually set to ZERO (that meaning no shunt resistors). If you make this a very high value (1e12) then it won't affect the simulation, but it will provide a finite resistance between all nodes, avoiding the singular matrix error.

\$\endgroup\$
2
  • 3
    \$\begingroup\$ This option should be used with care, as every node will be grounded by the gshunt conductance, no matter where, so it can add unwanted poles/zeroes that may ruin the results. E.g.: gshunt=1p with a 1pF capacitor. The same goes for cshunt. \$\endgroup\$ – a concerned citizen Sep 21 '16 at 12:08
  • \$\begingroup\$ I had to use GSHUNT in LTSpice but this solved it for me! To address the above concern, you can try varying the GSHUNT value by an order of magnitude to check if your results are sensitive. This will give an idea of whether it is affecting your results. \$\endgroup\$ – jalalipop May 19 at 19:59
1
\$\begingroup\$

The solution was given by endolith, but if somebody want to have isolated the secondary, he should ground it by a 1Meg or higher resistor, if not, the spice will not able to invert the array to solve the circuit.

The problem did not disappear for me so I restarted the LTspice and the problem was gone. Sometimes LTspice hangs in a simulation and then it works after it is restarted.

\$\endgroup\$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.