LTSpice simulazione circuito LC

[cianfa72]

No attento: quello che ti dà lui è il valore DC, non quello iniziale!

Scusami...cosa intendi per il valore DC della corrente dell'induttanza L1 ? Forse intendiamo cose diverse

Di seguito riporto i primi valori della simulazione (B): come si vede nella quinta colonna al tempo t=0 il valore iniziale per I(L1) e' proprio 10A

Codice: Seleziona tutto

time   V(n001)   I(C1)   I(I1)   I(L1)
0.000000000000000e+000   1.000000e+008   1.000000e-009   1.000000e+001   1.000000e+001
1.828571417458339e-007   9.997600e+007   -1.829629e+006   1.000000e+001   1.829639e+006
3.657142834916679e-007   9.991372e+007   -3.656007e+006   1.000000e+001   3.656017e+006
5.485714252375018e-007   9.981314e+007   -5.479133e+006   1.000000e+001   5.479143e+006
7.314285669833358e-007   9.967428e+007   -7.299009e+006   1.000000e+001   7.299019e+006
9.142857087291697e-007   9.949712e+007   -9.115633e+006   1.000000e+001   9.115643e+006
1.097142850475004e-006   9.928168e+007   -1.092901e+007   1.000000e+001   1.092902e+007
1.279999992220837e-006   9.902794e+007   -1.273913e+007   1.000000e+001   1.273914e+007
2.256562492220837e-006   9.723546e+007   -2.230868e+007   1.000000e+001   2.230869e+007
3.233124992220837e-006   9.459965e+007   -3.167299e+007   1.000000e+001   3.167300e+007
4.209687492220837e-006   9.106410e+007   -4.073507e+007   1.000000e+001   4.073508e+007
5.186249992220837e-006   8.666254e+007   -4.940872e+007   1.000000e+001   4.940873e+007
6.162812492220837e-006   8.143694e+007   -5.761148e+007   1.000000e+001   5.761149e+007
7.139374992220837e-006   7.543710e+007   -6.526534e+007   1.000000e+001   6.526535e+007


[gill90]

Il parametro uic di cui parlavi è rimasto uguale tra la prima simulazione che hai fatto e questa?

[cianfa72]

Il parametro uic di cui parlavi è rimasto uguale tra la prima simulazione che hai fatto e questa?

Si non e' stato mai presente come opzione nella scheda .TRAN

[cianfa72]

Tra l'altro per escludere eventuali problemi ho provato anche con la versione LTSpice IV: stesso comportamento

[gill90]

Secondo me non focalizzarti troppo sull'esempio B) perché non ha senso.
Cioè tu gli stai imponendo una condizione iniziale impropria a cui lui risponde con un risultato sbagliato, quindi potrebbe benissimo essere che LTspice dia quella soluzione e un altro simulatore ne dia una diversa. E nessuno dei due avrebbe comunque ragione.

Nemmeno a me torna quello che hai postato tu come log di simulazione: è come se lui calcolasse il valore a regime (supponendo l'induttanza in corto e il condensatore in aperto) e lo usasse come condizione iniziale in mancanza di altro per poi far evolvere la transient.
Ma questo solo perché il problema iniziale è mal posto.

[PietroBaima]

Ma che bella simulazione di una giunzione Josephson.

[gill90]

Pensa te.. Non si finisce mai di imparare

[cianfa72]

Cerca cerca ho trovato interessanti info non documentate sull'implementazione LTspice

http://ltwiki.org/index.php?title=Undocumented_LTspice

The default value for the parallel damping resistance equals 1e12 times the inductance value (1T x L) and is only applied in the case of a transient analysis.

In effetti con una induttanza di 10uH si ottiene una resistenza parallelo Rpar = 10*10^(-6)*10^(12) = 10*10^6 = 10 MOhm che giustifica il valore della soluzione DC operating point per la tensione del nodo N001 di 100MV = 10A * 10 MOhm.

Tra l'altro forzando a zero il valore Gpar = 1/Rpar per l'induttanza si ritrova che la soluzione per il nuovo DC operating point fornisce una tensione di nodo N001 di 10^18 volt. Questo torna con il fatto che la Gpar di default per un condensatore da 10uF vale Gpar = 10*10^(-6)*10^(-12) = 10*10^(-18) che con una corrente di 10A fornisce una tensione di (10/10)*10^18 = 10^18 volt

simulazione C: nello screenshot il puntatore del mouse si trova sopra il nodo N001

[gill90]

Interessante questa cosa, quindi lui forza automaticamente una resistenza in parallelo per poter simulare.. Non lo sapevo.
A titolo di esempio ti faccio vedere quello che invece restituisce SIMetrix se gli dò in pasto lo stesso circuito: lui si rifiuta proprio!

[cianfa72]

A titolo di esempio ti faccio vedere quello che invece restituisce SIMetrix se gli dò in pasto lo stesso circuito: lui si rifiuta proprio!

Si, se ho capito bene LTspice (SPICE in generale) usa il metodo MNA (Modified Nodal Analysis) per risolvere il circuito. Nel caso in esame, dal punto di vista della soluzione ITS, la condizione iniziale I=0 per l'induttore (al tempo t=0) modella l'induttore come un circuito aperto e pertanto la matrice delle conduttanze coinvolta nel sistema MNA e' assolutamente singolare e di fatto il nodo N001 e' floating