X-Authentication-Warning: delorie.com: mail set sender to geda-user-bounces using -f
X-Recipient: geda-user AT delorie DOT com
X-Injected-Via-Gmane: http://gmane.org/
To: geda-user AT delorie DOT com
From: Petro <x DOT piter AT gmail DOT com>
Subject: [geda-user] ngspice, strange voltage
Date: Sat, 02 Mar 2013 10:17:26 +0100
Lines: 68
Message-ID: <87wqtq2mah.fsf@cica.cica>
Mime-Version: 1.0
Content-Type: text/plain; charset=us-ascii
X-Complaints-To: usenet AT ger DOT gmane DOT org
X-Gmane-NNTP-Posting-Host: 194.228.32.161
User-Agent: Gnus/5.13 (Gnus v5.13) Emacs/23.2 (gnu/linux)
Cancel-Lock: sha1:LV/+hhQTNgpSdIGYDpKC6t1rqs8=
Reply-To: geda-user AT delorie DOT com

Hi all.

I try to model a tube ampllifier using ngspice.
Here is my netlist. Some how I get high voltage on the second grid (node
5).
Whos mistake is it, ngspice or mine?
Thanks
-------------------------------------------
* Spice netlister for gnetlist
* Cph 7 4 4.7n
* Ra2 9 4 2k
Ra2 7 4 20k
Cc2 0 8 10u
Rc2 8 0 10
* La 7 9 2H
X2 7 5 8 TRIODE_6S51N
V1 6 0 SIN(0 0.01 1kHz) 
V2 4 0 DC 50
Cgi 3 6 150n
Cout 5 1 150n
Rgi 3 0 1M
Rout 5 0 1M
Ra 4 1 10k
Rc 2 0 180
Cc 0 2 10u
X1 1 3 2 TRIODE_6S51N2
* Subcircuits
*------------------------------------------------------------
.func PWR(x,a) {abs(x) ** a}
.func PWRS(x,a) {sgn(x) * PWR(x,a)}
*____________________________________________________________
* triode             P G C
.SUBCKT TRIODE_6S51N 1 2 3 CCG=4.2P  CGP=2.2P CCP=1.6P RGI=2000  MU=51.8522 EX=1.3719 KG1=118.77  KP=110.41 KVB=383.7899 
*--------------------------------------------------
E1 7 0 VALUE={V(1,3)/KP*LN(1+EXP(KP*(1/MU+(V(2,3))/SQRT(KVB+V(1,3)*V(1,3)))))} 
RE1 7 0 1G   ; TO AVOID FLOATING NODES
G1 1 3 VALUE={(PWR(V(7),EX)+PWRS(V(7),EX))/KG1} 
RCP 1 3 1G   ; TO AVOID FLOATING NODES
C1 2 3 {CCG} ; CATHODE-GRID 
C2 2 1 {CGP} ; GRID=PLATE 
C3 1 3 {CCP} ; CATHODE-PLATE 
D3 5 3 DX    ; POSITIVE GRID CURRENT 
R1 2 5 {RGI} ; POSITIVE GRID CURRENT 
.MODEL DX D(IS=1N RS=1 CJO=10PF TT=1N) 
.ENDS 
.SUBCKT TRIODE_6S51N2 1 2 3 CCG=4.2P  CGP=2.2P CCP=1.6P RGI=2000  MU=51.8522 EX=1.3719 KG1=118.77  KP=110.41 KVB=383.7899 
*--------------------------------------------------
E1 7 0 VALUE={V(1,3)/KP*LN(1+EXP(KP*(1/MU+(V(2,3))/SQRT(KVB+V(1,3)*V(1,3)))))} 
RE1 7 0 1G   ; TO AVOID FLOATING NODES
G1 1 3 VALUE={(PWR(V(7),EX)+PWRS(V(7),EX))/KG1} 
RCP 1 3 1G   ; TO AVOID FLOATING NODES
C1 2 3 {CCG} ; CATHODE-GRID 
C2 2 1 {CGP} ; GRID=PLATE 
C3 1 3 {CCP} ; CATHODE-PLATE 
D3 5 3 DX    ; POSITIVE GRID CURRENT 
R1 2 5 {RGI} ; POSITIVE GRID CURRENT 
.MODEL DX D(IS=1N RS=1 CJO=10PF TT=1N) 
.ENDS 
*____________________________________________________________

* .DC Vp 0 100 .1 Vg 0.0 5.0 0.5
.TRAN 0.01ms 5ms 
*.op
.control
run
plot V(5) V(1)
.END