Need help to find the mistake.

Hello everybody!

I need a help to understand where is the mistake in my SPICE (programm: SPICE 2G.6) netlist. I have drawn an average-complexity (as for me) AC a circuit to test myself, and calculated everything what I have learned (except power), every property (voltages, currents, and their angles), but those values seems to be incorrect. However, when I check my circuit with EWB (programm: Electronics Workbench Version 5.12), it says me that everything is okay. But SPICE tells me that I am completely wrong!

Some images are too big for this forum, so I will paste links to them.

Here is my drawn circuit with netlist:

http://i42.tinypic.com/2i1cc29.jpg

Here are my attempts to solve the circuit (most calculation were made on scientific calculator, so no need simple formulas):

http://i41.tinypic.com/33a3slu.jpg

http://i40.tinypic.com/33uaekn.jpg

Here is the netlist:

circuit
v1 1 0 ac 24 sin
vit 1 2 ac 0
vilr 3 4 ac 0
vicr 3 6 ac 0
l1 2 3 560m
l2 4 5 22m
r1 5 0 470
c1 6 7 .33u
r2 7 0 910
.ac lin 1 50 50
.print ac v(2,3) vp(2,3) i(vit) ip(vit)
.print ac v(4,5) vp(4,5) i(vilr) ip(vilr)
.print ac v(5,0) vp(5,0) i(vilr) ip(vilr)
.print ac v(6,7) vp(6,7) i(vicr) ip(vicr)
.print ac v(7,0) vp(7,0) i(vicr) ip(vicr)
.print ac i(vit)
.end

Click to expand...

Here is what SPICE tells me:

***********************************************
V(2,3) VP(2,3) I(VIT) IP(VIT)
8.544E+00 7.107E+01 4.857E-02 -1.893E+01
***********************************************
V(4,5) VP(4,5) I(VILR) IP(VILR)
3.340E-01 6.831E+01 4.832E-02 -2.169E+01
***********************************************
V(5) VP(5) I(VILR) IP(VILR)
2.271E+01 -2.169E+01 4.832E-02 -2.169E+01
***********************************************
V(6,7) VP(6,7) I(VICR) IP(VICR)
2.261E+01 -2.623E+01 2.345E-03 6.377E+01
***********************************************
V(7) VP(7) I(VICR) IP(VICR)
2.134E+00 6.377E+01 2.345E-03 6.377E+01
***********************************************
I(VIT) <------------------------| total current
4.857E-02
***********************************************

Click to expand...

And here are the EWB results for currents only, and that is enough for me for thinking that I am on the right way:



What do you think?

I haven't done AC analyses in any SPICE for many years, so I could be wrong on this, but doesn't the .AC analysis generate an impulse response of the circuit in the frequency domain? Your EWB results look like steady state results in the time domain.

Hopefully that at least gives you someplace to explore.

Two comments from my side:
* AC analysis of a circuit using the SPICE code is performed NOT via an impulse response. It is based on classical calculations in the frequency domain.
* There is no ground reference in the shown circuit diagram.

Note that I'm only saying what it generates, not how it generates it. I believe that what it generates is the impulse response in the frequency domain, often called the transfer function.

Impulse response? Transfer function? Frequency domain? Time domain? May God has mercy!!! Will these words appear in the future chapters of e-book?!

* There is no ground reference in the shown circuit diagram.

Click to expand...

Look, I have just added ground to the circuit — nothing have changed.



I realised that I had done few analyses on the such AC circuits before, and hand-writed results always matched with results of two programms. Now only with one.

In your netlist I see a part called L1 that does not appear in the diagram.

LvW said:

L1 does not appear in the diagram.

Click to expand...

No, I think it is in place.

Sorry, my mistake. Im was blind.

I do a I quick calculation check of your and this is what I find

Ztot = XL1 + ((Xc1 + R1)*(XL2 + R2))/((Xc1 + R1) + (XL2 + R2)) = 899.226 + j98.4753

Ztot = 904.602Ω

and

Itot = 24V/904.602Ω = 0.026531A = 26.531mA

And
Ic1 = 0.00247957A = 2.47957mA
IL2 = 0.0263132A = 26.3132A

I think some mystery hiding inside this program when analyzing my circuit. I am very concerned that it does not give me the true results. Now as before I should always rely on my own calculations.
I solved the circuit successfully — isn't that a holiday?!
Thanks for all for taking a part in this discussion!