• 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

Hi all:

Does anyone have a model of the 6ES8-ecc189 triode? I can't seem to find one anywhere.

Thanks,
Code:
* ============================================================== 
* 6ES8 LTSpice model 
* Koren model (5 parameters): mean fit error 0.405269mA 
* Traced by Wayne Clay on 01/01/2019 using Engauge Digitizer 5.2 
* and Curve Captor v0.9.1 from Sylvania data sheet. 
* ============================================================== 
.subckt 6ES8  P G K 
    Bp  P K  I= 
+ (0.001667389859m)*uramp(V(P,K)*ln(1.0+exp((0.02294225145)+ 
+ (0.02294225145)*(3057.880457)*V(G,K)/sqrt((0.1158340263k)+ 
+ (V(P,K))**2)))/(0.02294225145))**(1.268488137) 
    Cgp G P 2.55p ; 0.7p added (1.85) 
    Cgk G K 4.0p ; 0.7p added (3.3p) 
    Cpk P K 0.38p ; 0.2p added (0.18p) 
    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 6ES8
Code:
* ============================================================== 
* 
* Generic triode model: ECC189_AN 
* Copyright 2003--2008 by Ayumi Nakabayashi, All rights reserved. 
* Version 3.10, Generated on Tue Dec 25 13:38:23 2018 
*                 Plate 
*                 | Grid 
*                 | | Cathode 
*                 | | | 
.SUBCKT ECC189_AN A G K 
BGG   GG   0 V=V(G,K)+0.32201735 
BM1   M1   0 V=(0.35943878*(URAMP(V(A,K))+1e-10))**-21.749142 
BM2   M2   0 V=(0.06451851*(URAMP(V(GG)+URAMP(V(A,K))/2.602617)+1e-10))**23.249142 
BP    P    0 V=0.0091825815*(URAMP(V(GG)+URAMP(V(A,K))/40.339075)+1e-10)**1.5 
BIK   IK   0 V=U(V(GG))*V(P)+(1-U(V(GG)))*1.6647528e+23*V(M1)*V(M2) 
BIG   IG   0 V=0.0045912908*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.0049348506*URAMP(V(A,K))**1.5)))+1e-10*V(A,K) 
BIGK  G    K I=V(IG) 
* CAPS 
CGA   G    A 1.9p 
CGK   G    K 3.5p 
CAK   A    K 1.7p 
.ENDS
 
Hi stijena1973
FYI:
My triode model is able to mimic the small initial velocity current. But in electrometer tubes, there are also gridcurrents in the sub-nA-range of interest, where further effects play a role. This area is still uncovered by the known spice model approaches...
All the best, Adrian
There is not much on the topic in the internet. The only mention of the component currents I saw is in "Electrometer tubes, part 1" by Dolezalek. But if the curve is known from the datasheet, maybe some kind of a look-up table is possible?
 
6N21B model

Hi guys! In case anyone needs model for Russian 6N21B (6Н21Б) triode, here it is, created by myself using model paint tools:


Code:
.SUBCKT 6N21B 1 2 3 ; Plate Grid Cathode
+ PARAMS: CCG=2.7P  CGP=0.6P CCP=1.4P RGI=2000
 + MU=113.16 KG1=585 KP=276 KVB=225 VCT=0.912 EX=1.526 
* Vp_MAX=320 Ip_MAX=14 Vg_step=0.5 Vg_start=0 Vg_count=8
* Rp=4000 Vg_ac=55 P_max=1 Vg_qui=-48 Vp_qui=300
* X_MIN=70 Y_MIN=11 X_SIZE=667 Y_SIZE=593 FSZ_X=2578 FSZ_Y=1458 XYGrid=true
* showLoadLine=y showIp=y isDHT=n isPP=n isAsymPP=n showDissipLimit=y 
* showIg1=n gridLevel2=n isInputSnapped=n  
* XYProjections=n harmonicPlot=n dissipPlot=n 
*----------------------------------------------------------------------------------
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  ; 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 
*$