• WARNING: Tube/Valve amplifiers use potentially LETHAL HIGH VOLTAGES.
    Building, troubleshooting and testing of these amplifiers should only be
    performed by someone who is thoroughly familiar with
    the safety precautions around high voltages.

Vacuum Tube SPICE Models

Triode-Connected
Code:
*
* Generic triode model: D3A_PT_AN
* Copyright 2003--2008 by Ayumi Nakabayashi, All rights reserved.
* Version 3.10, Generated on Tue Dec 24 08:47:57 2013
*                 Plate
*                 | Grid
*                 | | Cathode
*                 | | |
.SUBCKT D3A_PT_AN A G K
BGG   GG   0 V=V(G,K)+0.38087015
BM1   M1   0 V=(0.0021582473*(URAMP(V(A,K))+1e-10))**-0.28773641
BM2   M2   0 V=(0.83904987*(URAMP(V(GG)+URAMP(V(A,K))/74.574462)+1e-10))**1.7877364
BP    P    0 V=0.03039942*(URAMP(V(GG)+URAMP(V(A,K))/88.879654)+1e-10)**1.5
BIK   IK   0 V=U(V(GG))*V(P)+(1-U(V(GG)))*0.018902992*V(M1)*V(M2)
BIG   IG   0 V=0.01519971*URAMP(V(G,K))**1.5*(URAMP(V(G,K))/(URAMP(V(A,K))+URAMP(V(G,K)))*1.2+0.4)
BIAK  A    K I=URAMP(V(IK,IG)-URAMP(V(IK,IG)-(0.015714194*URAMP(V(A,K))**1.5)))+1e-10*V(A,K)
BIGK  G    K I=V(IG)
* CAPS
CGA   G    A 2.7p
CGK   G    K 7.3p
CAK   A    K 3.1p
.ENDS
 
Good!
And do not skip this beautiful Ale's triode model.
Btw, my collection of D3A`s.
 

Attachments

  • D3a-triode-SPICE.png
    D3a-triode-SPICE.png
    203.6 KB · Views: 286
  • D3A`s.jpg
    D3A`s.jpg
    59.3 KB · Views: 276
I'll throw my triode connected model into the ring. Made using Curve Captor from Siemens data sheet.
Code:
* D3a triode connected pentode LTspice model
* Modified Koren model (8 parameters): mean fit error 0.499405mA
* Traced by Wayne Clay on 6/24/2015 using Engauge Digitizer and
* Curve Captor v0.9.1 from Siemens datasheet.
.subckt D3a-t  P G K
    Bp  P K  I=(0.2505628718m)*uramp(V(P,K)*ln(1.0+(-0.1791266922)+exp((4.510443556)+(4.510443556)*((84.18566045)+(-3454.657477m)*V(G,K))*V(G,K)/sqrt((29.60063818)**2+(V(P,K)-(-4.621358765))**2)))/(4.510443556))**(1.176920934)
* Capacitances: G2+P, G3+K
    Cgk G K 8.3p ; 1.0p added (7.3p)
    Cpk P K 4.1p ; 1.0p added (3.1p)
    Cgp G P 2.9p ; 0.2p added (2.7p)
    Rpk P K 1G ; to avoid floating nodes
    d3  G K dx1
.model dx1 d(is=1n rs=2k cjo=1pf N=1.5 tt=1n)
.ends D3a-t
 

Attachments

  • D3a-t_CC.png
    D3a-t_CC.png
    32.6 KB · Views: 275
Ale to the rescue!

2P29L – Preamp and driver for 4P1L PSE Amp – Bartola(R) Valves

Pentode:
Code:
**** 2P29L DIRECTLY HEATED PENTODE ** Composite DHT with Advanced Grid Current **************
* Bartola Valves | All about electronic valves and hi-fi
* Plate Data source link: Bartola Valves | All about electronic valves and hi-fi
* Plate Curves image file: 2P29L-triode-model.jpg
*----------------------------------------------------------------------------------
.SUBCKT 2P29L_AM P G K1 K2 ; 1 2 3 4 
+ PARAMS: CCG=6P CGP=8P CCP=4P RFIL=13.5
+ MU=9.1 KG1=2445 KP=140 KVB=21 VCT=-0.46 EX=1.4 
+ VGOFF=-0.6 IGA=0.001 IGB=0.3 IGC=8 IGEX=2 
* Vp_MAX=350 Ip_MAX=140 Vg_step=5 Vg_start=10 Vg_count=9
* Rp=4000 Vg_ac=55 P_max=9 Vg_qui=-48
* X_MIN=77 Y_MIN=89 X_SIZE=492 Y_SIZE=526 FSZ_X=1157 FSZ_Y=724 XYGrid=false
*----------------------------------------------------------------------------------
RFIL_LEFT 3 31 {RFIL/4}
RFIL_RIGHT 4 41 {RFIL/4}
RFIL_MIDDLE1 31 34 {RFIL/4}
RFIL_MIDDLE2 34 41 {RFIL/4}
E11 32 0 VALUE={V(1,31)/KP*LOG(1+EXP(KP*(1/MU+V(2,31)/SQRT(KVB+V(1,31)*V(1,31)))))}
E12 42 0 VALUE={V(1,41)/KP*LOG(1+EXP(KP*(1/MU+V(2,41)/SQRT(KVB+V(1,41)*V(1,41)))))}
RE11 34 0 1G
G11 1 31 VALUE={(PWR(V(32),EX)+PWRS(V(32),EX))/(2*KG1)}
G12 1 41 VALUE={(PWR(V(42),EX)+PWRS(V(42),EX))/(2*KG1)}
RCP1 1 34 1G
C1 2 34 {CCG} ; CATHODE-GRID
C2 2 1 {CGP} ; GRID=PLATE
C3 1 34 {CCP} ; CATHODE-PLATE
RE2 2 0 1G
EGC1 81 0 VALUE={V(2,31)-VGOFF} ; POSITIVE GRID THRESHOLD
GG1 2 31 VALUE={0.5*(IGA+IGB/(IGC+V(1,31)))*(MU/KG1)*(PWR(V(81),IGEX)+PWRS(V(81),IGEX))}
EGC2 82 0 VALUE={V(2,41)-VGOFF} ; POSITIVE GRID THRESHOLD
GG2 2 41 VALUE={0.5*(IGA+IGB/(IGC+V(1,41)))*(MU/KG1)*(PWR(V(82),IGEX)+PWRS(V(82),IGEX))}
.ENDS 
*$

Triode-wired:
Code:
** 2P29L TRIODE ************************************************************
* Created on Fri Jun 15 19:29:17 BST 2012 using tube.model.finder.PaintKIT
* model URL:[url=http://www.bartola.co.uk/valves]Bartola® Valves – All about electronic valves and hi-fi[/url]
*--------------------------------------------------
.SUBCKT 2P29L-T_AM 1 2 3 ; A G K  
+ PARAMS: CCG=6P  CGP=8P CCP=4P RGI=2000
+ MU=9.1 EX=1.4 KG1=2445 KP=140 KVB=21 VCT=-0.46 ; Vp_MAX=200.0 Ip_MAX=0.08 Vg_step=2.0
*--------------------------------------------------
E1 7 0 VALUE={V(1,3)/KP*LOG(1+EXP(KP*(1/MU+(VCT+V(2,3))/SQRT(KVB+V(1,3)*V(1,3)))))} 
RE1 7 0 1G 
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 ; FOR GRID CURRENT 
R1 2 5 {RGI} ; FOR GRID CURRENT 
.MODEL DX D(IS=1N RS=1 CJO=10PF TT=1N) 
.ENDS 
*$

--
 
6267, EF86 equivalent.
Code:
*
* Generic pentode model: 6267
* Copyright 2003--2008 by Ayumi Nakabayashi, All rights reserved.
* Version 3.10, Generated on Sat Mar  8 22:42:16 2008
*               Plate
*               | Screen Grid
*               | |  Control Grid
*               | |  |  Cathode
*               | |  |  |
.SUBCKT 6267_AN A G2 G1 K
BGG   GG   0 V=V(G1,K)+0.59868749
BM1   M1   0 V=(0.010782364*(LIMIT(V(G2,K),0,1e16)+1e-10))^-0.70765893
BM2   M2   0 V=(0.67945278*(LIMIT((V(GG)+LIMIT(V(G2,K),0,1e16)/29.728844),0,1e16)))^2.2076589
BP    P    0 V=0.0013378994*(LIMIT((V(GG)+LIMIT(V(G2,K),0,1e16)/43.754099),0,1e16))^1.5
BIK   IK   0 V=STP(V(GG))*V(P)+(1-STP(V(GG)))*0.00078620809*V(M1)*V(M2)
BIG   IG   0 V=0.00066894969*LIMIT(V(G1,K),0,1e16)^1.5*(LIMIT(V(G1,K),0,1e16)/(LIMIT(V(A,K),0,1e16)+LIMIT(V(G1,K),0,1e16))*1.2+0.4)
BIK2  IK2  0 V=V(IK,IG)*(1-0.4*(EXP(-LIMIT(V(A,K),0,1e16)/LIMIT(V(G2,K),0,1e16)*15)-EXP(-15)))
BIG2T IG2T 0 V=V(IK2)*(0.83966688*(1-LIMIT(V(A,K),0,1e16)/(LIMIT(V(A,K),0,1e16)+10))^1.5+0.16033312)
BIK3  IK3  0 V=V(IK2)*(LIMIT(V(A,K),0,1e16)+7510)/(LIMIT(V(G2,K),0,1e16)+7510)
BIK4  IK4  0 V=V(IK3)-LIMIT(V(IK3)-(0.00071507731*(LIMIT(V(A,K),0,1e16)+LIMIT(LIMIT(V(G2,K),0,1e16)-LIMIT(V(A,K),0,1e16),0,1e16))^1.5),0,1e16)
BIP   IP   0 V=LIMIT(V(IK4,IG2T)-LIMIT(V(IK4,IG2T)-(0.00071507731*LIMIT(V(A,K),0,1e16)^1.5),0,1e16),0,1e16)
BIAK  A    K I=V(IP)+1e-10*V(A,K)
BIG2  G2   K I=LIMIT(V(IK4,IP),0,1e16)
BIGK  G1   K I=V(IG)
* CAPS
CGA   G1  A  0.05p
CGK   G1  K  2.3p
C12   G1  G2 1.5p
CAK   A   K  5.3p
.ENDS

I'm getting really wacky results in a known working circuit. using LTSPICE VII on mac. Duncan's works. Am I doing something wrong using pentode2.asy? Basically the sim shows -54 dB at 100K and sloping down to 20 Hz
 
Try this:

Code:
.SUBCKT NU6P1 1 2 3
+params: MU= 18.10 EX= 4.080 KG1=4270851.9 KP=451.94 KVB= 4.2 VCT= 0.00 RGI=330k CCG=9.1P CGP=2.5P CCP=4.3P
E1 7 0 VALUE= 
+{V(1,3)/KP*LOG(1+EXP(KP*(1/MU+V(2,3)/SQRT(KVB+V(1,3)*V(1,3)))))} 
RE1 7 0 1G 
G1 1 3 VALUE={(PWR(V(7),EX)+PWRS(V(7),EX))/KG1} 
RCP 1 3 1G ; TO AVOID FLOATING NODES IN MU-FOLLOWER 
C1 2 3 {CCG} ; CATHODE-GRID 
C2 2 1 {CGP} ; GRID=PLATE 
C3 1 3 {CCP} ; CATHODE-PLATE 
D3 5 3 DX ; FOR GRID CURRENT 
R1 2 5 {RGI} ; FOR GRID CURRENT 
.MODEL DX D(IS=1N RS=1 CJO=10PF TT=1N) 
.ENDS

Thanks for providing the Nutube model Jazbo.

Anybody know how I can make this connection marked in red in LTSpice?
 

Attachments

  • LT Model.jpg
    LT Model.jpg
    43.6 KB · Views: 291