Vacuum Tube SPICE Models

Adrian Immler · 2018-10-28 2:38 pm

Set of new approaches with accurate grid current

Hello

God news for all looking for tube spice models with accurate grid current!
I'm proud to introduce a set of new approaches which represents the merge of the very best ideas found in the web to this topic. It's the result of a two-year-project of mine.

Up to now, I have fitted 6 triodes: 6021, 6111, 6948, 7963, ECC81 and ECC86

You can find the models and also diagrams showing the fit quality here:
http://adrianimmler.simplesite.com/440951057

Any feedback, your own fits of tubes (using my approach) or any improvement proposals are highly welcome!

kind regards, Adrian Immler

bea · 2018-10-28 9:36 pm

Really Nice!

tikiroo said:
ECF802 Triode
....
ECF802 Pentode

Thanks.
Are these self contained?
Will these run in LTSpice, or will they need edits?
I am not familiar enough with LTSPICE - (how) can i use the models with the already supplied symbols for triode and pentode?

tikiroo · 2018-10-29 5:42 am

Yes they are self contained and formatted for LTSpice. Add a new Spice Directive, copy and paste the model text into it. Add a standard triode symbol, right click on it and on the "SpiceModel" line, type the name of the model e.g ECF802T. Similar with tetrode symbol.

jazbo8 · 2018-10-29 6:18 am

Adrian Immler said:
Hello

God news for all looking for tube spice models with accurate grid current!
I'm proud to introduce a set of new approaches which represents the merge of the very best ideas found in the web to this topic. It's the result of a two-year-project of mine.

Up to now, I have fitted 6 triodes: 6021, 6111, 6948, 7963, ECC81 and ECC86

You can find the models and also diagrams showing the fit quality here:
http://adrianimmler.simplesite.com/440951057

Any feedback, your own fits of tubes (using my approach) or any improvement proposals are highly welcome!

kind regards, Adrian Immler

Thanks for sharing your innovative work and the improved SPICE models! As a suggestion, perhaps you can add a reference section at the end with links to the articles/papers/books that you sited in your paper. Also you may be interested in Dr. Reefman's paper: http://www.dos4ever.com/uTracer3/Theory.pdf

Adrian Immler · 2018-10-29 6:49 am

Hello jazbo8,
yes a reference section would be a good idea. It will need some time because my priorities have changed now.
I also studied Derk Reefmans great paper, it is a "must" when it comes to tetrodes/pentodes, especially his way to mimic second emission effects. Regarding Korens triode model, he was too fast satisfied in my opinion ;-).

kind regards, Adrian Immler

tikiroo · 2018-10-29 7:27 am

tikiroo said:
Add a standard triode symbol, right click on it and on the "SpiceModel" line, type the name of the model e.g ECF802T. Similar with tetrode symbol.

Model name should actually go in "Value" line, not "SpiceModel" line, sorry

.

See example attached.

bea · 2018-10-29 12:09 pm

thanks.

rufus74pz · 2018-10-30 7:43 pm

6HE7

6HE7 in triode

.SUBCKT 6HE7 1 2 3 ; Plate Grid Cathode
+ PARAMS: CCG=19P CGP=0.38P CCP=8P RGI=2000
+ MU=5.347 KG1=442.5 KP=18.56 KVB=0.1945 VCT=5.064E-4 EX=1.457
* Vp_MAX=400 Ip_MAX=200 Vg_step=7 Vg_start=0 Vg_count=15
* Rp=4000 Vg_ac=55 P_max=40 Vg_qui=-48 Vp_qui=300
* X_MIN=35 Y_MIN=10 X_SIZE=532 Y_SIZE=654 FSZ_X=1085 FSZ_Y=738 XYGrid=false
* showLoadLine=n 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
*$

rufus74pz · 2018-10-30 9:51 pm

I'm looking for a model for 6GE5 in triode?

choique · 2018-10-30 11:39 pm

Hi,
I'm looking for LTspice Model of
6BM8/ECL82 Pentode Section's Triode Connection.
Thanks

RajkoM said:
Means Triode section. And here is Pentode section:
(Of course, pay attention to the Ayumi 6BM8 models)

Many Thanks for Yr Supports,

But I'm Looking for 6BM8/ECL82 Pentode Section's Triode connection Data.

Hi tikiroo,
U R Hero..!!
Thanks,

rongon · 2018-10-31 12:02 am

rufus74pz said:
I'm looking for a model for 6GE5 in triode?

There's this:
http://www.bartola.co.uk/valves/2015/09/19/6ge5-pentode-model/

According to this,
"6JN6 is equivalent to 6GE5, 6FW5, 6DQ6B and 6GV5. So I would check to see if something is available for one of them."

If that's true, then these LTspice Ayumi models should work...

6DQ6 Pentode:

Code:

*
* Generic pentode model: 6DQ6
* Copyright 2003--2008 by Ayumi Nakabayashi, All rights reserved.
* Version 3.10, Generated on Sat Mar  8 22:39:50 2008
*            Plate
*            | Screen Grid
*            | | Control Grid
*            | | | Cathode
*            | | | |
.SUBCKT 6DQ6 A G2 G1 K
BGG   GG   0 V=V(G1,K)+-1
BM1   M1   0 V=(0.17812317*(URAMP(V(G2,K))+1e-10))**-1.6443126
BM2   M2   0 V=(0.47705181*(URAMP(V(GG)+URAMP(V(G2,K))/2.9358797)))**3.1443126
BP    P    0 V=0.0034562529*(URAMP(V(GG)+URAMP(V(G2,K))/6.1542156))**1.5
BIK   IK   0 V=U(V(GG))*V(P)+(1-U(V(GG)))*0.0040197798*V(M1)*V(M2)
BIG   IG   0 V=0.0017281264*URAMP(V(G1,K))**1.5*(URAMP(V(G1,K))/(URAMP(V(A,K))+URAMP(V(G1,K)))*1.2+0.4)
BIK2  IK2  0 V=V(IK,IG)*(1-0.4*(EXP(-URAMP(V(A,K))/URAMP(V(G2,K))*15)-EXP(-15)))
BIG2T IG2T 0 V=V(IK2)*(0.980449351*(1-URAMP(V(A,K))/(URAMP(V(A,K))+10))**1.5+0.019550649)
BIK3  IK3  0 V=V(IK2)*(URAMP(V(A,K))+2287.5)/(URAMP(V(G2,K))+2287.5)
BIK4  IK4  0 V=V(IK3)-URAMP(V(IK3)-(0.0026038839*(URAMP(V(A,K))+URAMP(URAMP(V(G2,K))-URAMP(V(A,K))))**1.5))
BIP   IP   0 V=URAMP(V(IK4,IG2T)-URAMP(V(IK4,IG2T)-(0.0026038839*URAMP(V(A,K))**1.5)))
BIAK  A    K I=V(IP)+1e-10*V(A,K)
BIG2  G2   K I=URAMP(V(IK4,IP))
BIGK  G1   K I=V(IG)
* CAPS
CGA   G1  A  0.5p
CGK   G1  K  8.7p
C12   G1  G2 5.8p
CAK   A   K  6.5p
.ENDS

6DQ6 Triode:

Code:

*
* Generic triode model: 6DQ6T
* Copyright 2003--2008 by Ayumi Nakabayashi, All rights reserved.
* Version 3.10, Generated on Sat Mar  8 22:39:50 2008
*             Plate
*             | Grid
*             | | Cathode
*             | | |
.SUBCKT 6DQ6T A G K
BGG   GG   0 V=V(G,K)+-1
BM1   M1   0 V=(0.17812317*(URAMP(V(A,K))+1e-10))**-1.6443126
BM2   M2   0 V=(0.47705181*(URAMP(V(GG)+URAMP(V(A,K))/2.9358797)+1e-10))**3.1443126
BP    P    0 V=0.0034562529*(URAMP(V(GG)+URAMP(V(A,K))/6.1542156)+1e-10)**1.5
BIK   IK   0 V=U(V(GG))*V(P)+(1-U(V(GG)))*0.0040197798*V(M1)*V(M2)
BIG   IG   0 V=0.0017281264*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.0026038839*URAMP(V(A,K))**1.5)))+1e-10*V(A,K)
BIGK  G    K I=V(IG)
* CAPS
CGA   G    A 6.3p
CGK   G    K 8.7p
CAK   A    K 6.5p
.ENDS

rufus74pz · 2018-10-31 10:24 am

Thanks.

Adrian Immler · 2018-10-31 11:42 am

jazbo8 said:
Thanks for sharing your innovative work and the improved SPICE models! As a suggestion, perhaps you can add a reference section at the end with links to the articles/papers/books that you sited in your paper. Also you may be interested in Dr. Reefman's paper: http://www.dos4ever.com/uTracer3/Theory.pdf

Hi jazbo8
I have added now your proposed reference section at the end of my paper

:
http://adrianimmler.simplesite.com/440956786

kind regards, Adrian

jazbo8 · 2018-10-31 11:43 am

Awesome!

rufus74pz · 2018-11-02 8:20 am

EF55 in triode

EF55 in triode

.SUBCKT EF55T 1 2 3 ; Plate Grid Cathode
+ PARAMS: CCG=15P CGP=12P CCP=0.15P RGI=2000
+ MU=24.86 KG1=372.4 KP=148 KVB=729 VCT=0.316 EX=1.414
* Vp_MAX=400 Ip_MAX=70 Vg_step=2 Vg_start=0 Vg_count=10
* Rp=4000 Vg_ac=55 P_max=12 Vg_qui=-48 Vp_qui=300
* X_MIN=51 Y_MIN=20 X_SIZE=1057 Y_SIZE=687 FSZ_X=1638 FSZ_Y=762 XYGrid=false
* showLoadLine=n 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
*$

bluerooster · 2018-11-05 9:21 pm

jazbo8 said:
The following models are generated from a curve traced sample of 813, so they do not match up well with the GK-71 nor the 813 datasheets, use them at your risk...

The triode one seems to work fine. But the "pentode" only has 4 elements, where's the supressor grid? Also, need a symbol to go with it, as the 813 has no seperate cathode, but a filament instead, which is also the cathode.

jazbo8 · 2018-11-07 10:49 am

The suppressor grid is usually not modeled as it is typically tied to the cathode. You can also make your own schematic symbol for the filament/cathode if that's what you want to see, it does not matter to the model.

Adrian Immler · 2018-11-07 11:16 am

Adrian Immler said:
Hello

God news for all looking for tube spice models with accurate grid current!
I'm proud to introduce a set of new approaches which represents the merge of the very best ideas found in the web to this topic. It's the result of a two-year-project of mine.

Up to now, I have fitted 6 triodes: 6021, 6111, 6948, 7963, ECC81 and ECC86

You can find the models and also diagrams showing the fit quality here:
http://adrianimmler.simplesite.com/440951057

Any feedback, your own fits of tubes (using my approach) or any improvement proposals are highly welcome!

kind regards, Adrian Immler

Hello
I just put all my triode models in a zip file, so it is now more convenient to download them. Thank you Retrovert for your tip!
http://adrianimmler.simplesite.com/440951057

kind regards, Adrian Immler

artosalo · 2018-11-07 4:02 pm

Could any of you gurus create the LT Spice model from these datas ?
These, 12E1 and S11E12, are identical, but S11E12 is a Hi-rel type and therefore has lower Pd_max limit.

Koonw · 2018-11-07 8:05 pm

artosalo said:
Could any of you gurus create the LT Spice model from these datas ?
These, 12E1 and S11E12, are identical, but S11E12 is a Hi-rel type and therefore has lower Pd_max limit.

Code:

**** S11E12_P ******************************************
* Created on 11/08/2018 03:53 using paint_kip.jar
* [url=http://www.dmitrynizh.com/tubeparams_image.htm]Model Paint Tools: Trace Tube Parameters over Plate Curves, Interactively[/url]
* Plate Curves image file: s11e12-p.png
* Data source link: <plate curves URL>
*----------------------------------------------------------------------------------
.SUBCKT S11E12_P P G2 G K ; LTSpice tetrode.asy pinout
* .SUBCKT S11E12_P P G K G2 ; Koren Pentode Pspice pinout
+ PARAMS: MU=5.5 KG1=1235.52 KP=34.63 KVB=0 VCT=2.687 EX=1.426 KG2=1538.46 KNEE=10.54 KVC=1.635
+ KLAM=7.424E-5 KLAMG=8.311E-5   KNK=-0.044 KNG=0.006 
+ CCG=19.5P CGP=1.8P CCP=16.5P RGI=2000.0
* Vp_MAX=300 Ip_MAX=350 Vg_step=4 Vg_start=0 Vg_count=13
* X_MIN=101 Y_MIN=102 X_SIZE=556 Y_SIZE=648 FSZ_X=1550 FSZ_Y=878 XYGrid=false
* Rp=1400 Vg_ac=20 P_max=28 Vg_qui=-24 Vp_qui=300 
* showLoadLine=n showIp=y isDHP=n isPP=n isAsymPP=n isUL=n showDissipLimit=y 
* showIg1=n isInputSnapped=y addLocalNFB=n
* XYProjections=n harmonicPlot=y dissipPlot=n 
* UL=0.43 EG2=150 gridLevel2=n addKink=n isTanhKnee=n advSigmoid=n 
*----------------------------------------------------------------------------------
RE1  7 0  1G    ; DUMMY SO NODE 7 HAS 2 CONNECTIONS
E1   7 0  VALUE=  ; E1 BREAKS UP LONG EQUATION FOR G1.
+{V(G2,K)/KP*LOG(1+EXP((1/MU+(VCT+V(G,K))/SQRT(KVB+V(G2,K)*V(G2,K)))*KP))}
RE2  6 0  1G    ; DUMMY SO NODE 6 HAS 2 CONNECTIONS
E2  6 0  VALUE={(PWR(V(7),EX)+PWRS(V(7),EX))} ; Kg1 times KIT current
G1   P K  VALUE={V(6)/KG1*ATAN(V(P,K)/KNEE)*(1+KLAMG*V(P,K))+KLAM*V(P,K)}
* Alexander Gurskii screen current, see audioXpress 2/2011
RE4K 4K K 1G ; Dummy, per Alex request
E4K 4K  4 VALUE={0} ; Dummy, per Alex request
G4K 4K K VALUE={V(6)/KG2*(KVC-ATAN(V(P,K)/KNEE))/(1+KLAMG*V(P,K))}
RCP  P K  1G     ; FOR CONVERGENCE
C1   K G  {CCG}  ; CATHODE-GRID 1
C2   G P  {CGP}  ; GRID 1-PLATE
C3   K P  {CCP}  ; CATHODE-PLATE
R1   G 5  {RGI}  ; FOR GRID CURRENT   
D3   5 K  DX     ; FOR GRID CURRENT } 
.MODEL DX D(IS=1N RS=1 CJO=10PF TT=1N)
.ENDS
*$

 
* The following triode model is derived from pentode model, see above.
* In the triode model, all spice parameters come directly from the pentode model, except for Kg1,
* which for triode-strapped pentodes is derived from pentode's Kg1, Kg2 and Kvc as
*
*    4Kg1Kg2 / ((2Kvc-Pi)(2Kg1+PiKg2))

**** S11E12_P ******************************************
* Created on 11/08/2018 03:53 using paint_kit.jar 4.7 
* [url=http://www.dmitrynizh.com/tubeparams_image.htm]Model Paint Tools: Trace Tube Parameters over Plate Curves, Interactively[/url]
* Plate Curves image file: s11e12-p.png
* Data source link: <plate curves URL>
*----------------------------------------------------------------------------------
.SUBCKT TRIODE_S11E12_P 1 2 3 ; Plate Grid Cathode
+ PARAMS: CCG=3P  CGP=1.4P CCP=1.9P RGI=2000
+ MU=5.5 KG1=8106.89 KP=34.63 KVB=0 VCT=2.687 EX=1.426 
* Vp_MAX=300 Ip_MAX=350 Vg_step=4 Vg_start=0 Vg_count=13
* Rp=1400 Vg_ac=20 P_max=28 Vg_qui=-24 Vp_qui=300
* X_MIN=101 Y_MIN=102 X_SIZE=556 Y_SIZE=648 FSZ_X=1550 FSZ_Y=878 XYGrid=false
* showLoadLine=n showIp=y isDHT=n isPP=n isAsymPP=n showDissipLimit=y 
* showIg1=n gridLevel2=n isInputSnapped=y  
* XYProjections=n harmonicPlot=y 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 
*$

Vacuum Tube SPICE Models

Member

Member

Member

Ex-Moderator

Member

Member

Attachments

Member

Member

Member

Member

Member

Member

Member

Ex-Moderator

Member

Member

Ex-Moderator

Member

Member

Attachments

Disabled Account

Attachments