Vacuum Tube SPICE Models

OldHector · 2021-02-06 10:58 am

Probably a daft question - If one has a model of a pentode. and it is triode wired, is the performance of the tube still model-able if the connection is made explicitly in LTSpice, from g2 to anode? I suppose I am trying to understand the rationale for seperate models for triode-connected pentodes.

I derived a triode model for EF12 using Dmitry Nizhegorodov's curve mapping tool. I could only manage a close match, after spending a long time on finding a correlation. I like this empirical approach to validating designs, and I am interested in how people build confidence in the models they are using.

Adrian Immler · 2021-02-07 5:53 am

Hi OldHector,

if your pentode model delivers correct results for ANY Vg2, you don't need an extra triode-connected model.

regards, Adrian

OldHector · 2021-02-08 2:23 pm

Adrian Immler said:
if your pentode model delivers correct results for ANY Vg2, you don't need an extra triode-connected model.

Hi Adrian, How do I confirm that the model incorporates Vg2? Should I try to emulate curves where Vg2 is changed, and look for a correlation?
Cheers, Richard

tikiroo · 2021-02-08 6:36 pm

The pentode EF12 model I posted in #2810 incorporates g2. However the datasheet only shows curves for Vg2 = 100. I did check the triode connected curves and they looked reasonable to me, but this is a case where a specific triode model may be better.

Adrian Immler · 2021-02-08 8:09 pm

Hi Richard, hi tikiroo

The Telefunken datasheet on frank's site also showes triode mode curves. So one can easy check how accurate your triode-connected EF12 pentode model is.v

cheers, Adrian

tikiroo · 2021-02-08 10:56 pm

Yes, I did a quick visual check against the Telefunken curves but didn't go to the trouble of overlaying the model on the published curves. Richard could do that if he wanted.

OldHector · 2021-02-08 11:17 pm

That sounds like a good exercise in improving my understanding! I'll give it a bash tomorrow. Thanks for the pointers!

Koonw · 2021-02-09 2:39 pm

Here is 807/1625 RCA model, 807 heater is 6.3V, 1625 is 12.6V and pin-out is different. Also attached is Extract V3.000 1625 model but the positive range is somewhat out. These models are similar to 6L6 but more extended Vg1 range.

Code:

 *** 807 ************************************G4******
* Created on 02/09/2021 22:16 using paint_kip.jar
* [URL="http://www.dmitrynizh.com/tubeparams_image.htm"]www.dmitrynizh.com/tubeparams_image.htm[/URL]
* Plate Curves image file: 807.png
* Data source link: <plate curves URL>
*----------------------------------------------------------------------------------
.SUBCKT 807 P G2 G K ; LTSpice tetrode.asy pinout
* .SUBCKT 807 P G K G2 ; Koren Pentode Pspice pinout
+ PARAMS: MU=8.952 KG1=1178.02 KP=42.4 KVB=1811.65 VCT=-2.18 EX=1.277 KG2=1106.5 KNEE=8.776 KVC=1.809
+ KLAM=3.125E-9 KLAMG=2.935E-4  KD=45.73 KC=10.12 KR1=0.001196 KR2=0.0256 KVBG=0.024 KB1=5.6 KB2=2.12 KB3=1.24 KB4=0.285 KVBGI=6.104E-9 KNK=0.27 KNG=0.02253 KNPL=0.007808 KNSL=0.00325 KNPR=0.379 KNSR=29.18
+ CCG=12P CGP=0.2P CCP=7P VGOFF=-2.25 IGA=0.003648 IGB=0.039 IGC=0.1 IGEX=1.88
* Vp_MAX=800 Ip_MAX=500 Vg_step=5 Vg_start=30 Vg_count=17
* X_MIN=53 Y_MIN=66 X_SIZE=717 Y_SIZE=561 FSZ_X=1296 FSZ_Y=736 XYGrid=true
* Rp=1400 Vg_ac=20 P_max=25 Vg_qui=-10 Vp_qui=300 
* showLoadLine=n showIp=y isDHP=n isPP=n isAsymPP=n isUL=n showDissipLimit=y 
* showIg1=y isInputSnapped=y addLocalNFB=n
* XYProjections=n harmonicPlot=y dissipPlot=n 
* UL=0.43 EG2=250 gridLevel2=y addKink=y isTanhKnee=n advSigmoid=y 
*----------------------------------------------------------------------------------
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
E4   8 0  VALUE={V(P,K)/KNEE/(KVBGI+V(6)*KVBG)}
E5  81 0  VALUE={PWR(V(8),KB1)}
E6  82 0  VALUE={PWR(V(8),KB2)}
E7  83 0  VALUE={PWR(V(8),KB3)}
E8   9 0  VALUE={PWR(1-EXP(-V(81)*(KC+KR1*V(82))/(KD+KR2*V(83))),KB4)*1.5708}
RE4  8 0  1
RE5 81 0  1
RE6 82 0  1
RE7 83 0  1
RE8  9 0  1
RE21 21 0 1
E21  21 0 VALUE={V(6)/KG1*V(9)} ; Ip with knee but no slope and no kink
RE22 22 0 1 ; E22: kink curr deviation for plate
E22  22 0 VALUE={V(21)*LIMIT(KNK-V(G,K)*KNG,0,0.3)*(-ATAN((V(P,K)-KNPL)/KNSL)+ATAN((V(P,K)-KNPR)/KNSR))} 
G1   P K  VALUE={V(21)*(1+KLAMG*V(P,K))+KLAM*V(P,K) + V(22)}
G2   G2 K  VALUE={V(6)/KG2*(KVC-V(9))/(1+KLAMG*V(P,K)) - V(22)}
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
RE23 G 0 1G 
GG G K VALUE={(IGA+IGB/(IGC+V(P,K)))*(MU/KG1)*
+(PWR(V(G,K)-VGOFF,IGEX)+PWRS(V(G,K)-VGOFF,IGEX))}
.ENDS
*$

Code:

**************************************************** 
.SUBCKT 1625 1 2 3 4 ; A G2 G1 C; 
* Extract V3.000 * Model created: 29-Apr-2020 * 
X1 1 2 3 4 BTetrodeDE MU= 9.64 EX=1.372 kG1= 519.7 KP= 42.1 kVB = 4096.6 kG2= 7370.3 
+Sc=.26E-01 ap= .018 w= 20. nu= 3.12 lam= 6.2 
+ Ookg1mOokG2=.179E-02 Aokg1=.11E-05 alkg1palskg2=.179E-02 be= .057 als= 6.85 RGI=2000 
+ CCG1=11P CCG2 = 0.0p CPG1 = 0.2p CG1G2 = 0.0p CCP=7P ; 
.ENDS 
**************************************************** 
.SUBCKT BTetrodeDE 1 2 3 4; A G2 G1 C 
* * NOTE: LOG(x) is base e LOG or natural logarithm. 
* For some Spice versions, e.g. MicroCap, this has to be changed to LN(x). 
* 
RE1 7 0 1MEG ; DUMMY SO NODE 7 HAS 2 CONNECTIONS 
E1 7 0 VALUE= 
+{V(2,4)/KP*LOG(1+EXP(KP*(1/MU+V(3,4)/SQRT(KVB+V(2,4)*V(2,4)))))} 
E2 8 0 VALUE = {Ookg1mOokG2 + Aokg1*V(1,4) - alkg1palskg2*Exp(-be*V(1,4)*SQRT(be*V(1,4)))} 
E3 9 0 VALUE = {Sc/kG2*V(1,4)*(1+tanh(-ap*(V(1,4)-V(2,4)/lam+w+nu*V(3,4))))} 
G1 1 4 VALUE = {0.5*(PWR(V(7),EX)+PWRS(V(7),EX))*(V(8)-V(9))} 
G2 2 4 VALUE = {0.5*(PWR(V(7),EX)+PWRS(V(7),EX))/KG2 *(1+als*Exp(-be*V(1,4) * SQRT(be*V(1,4))))} 
RCP 1 4 1G ; FOR CONVERGENCE A - C 
C1 3 4 {CCG1} ; CATHODE-GRID 1 C - G1 
C4 2 4 {CCG2} ; CATHODE-GRID 2 C - G2 
C5 2 3 {CG1G2} ; GRID 1 -GRID 2 G1 - G2 
C2 1 3 {CPG1} ; GRID 1-PLATE G1 - A 
C3 1 4 {CCP} ; CATHODE-PLATE A - C 
R1 3 5 {RGI} ; FOR GRID CURRENT G1 - 5 
D3 5 4 DX ; FOR GRID CURRENT 5 - C 
.MODEL DX D(IS=1N RS=1 CJO=10PF TT=1N) 
.ENDS BTetrodeDE

Koonw · 2021-02-10 9:35 am

This 807 model is more robust than the previous model as it can be ran in more schematics. I disable some of the options like making KB1, KB3 =0 in Advance Knee Params. The makes the model less complex and easier to be handled by Spice software. Also I recommend using the following option in LTSpice:

.options GminSteps=0
.options NoOpIter
.options SrcSteps=380000

This will make the source steps very small and allow Spice to have sufficient step to complete simulation without premature quitting. It'll take a little longer but it will be completed.

Code:

**** 807 ******************************************
* Created on 02/10/2021 17:34 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: 807.png
* Data source link: <plate curves URL>
*----------------------------------------------------------------------------------
.SUBCKT 807 P G2 G K ; LTSpice tetrode.asy pinout
* .SUBCKT 807 P G K G2 ; Koren Pentode Pspice pinout
+ PARAMS: MU=8.952 KG1=471.94 KP=42.36 KVB=1630.48 VCT=-7.89 EX=1.08 KG2=353.2 KNEE=0.3401 KVC=1.809
+ KLAM=3.047E-10 KLAMG=2.398E-4  KD=1211.72 KC=359999.82 KR1=0.1866 KR2=4.125E6 KVBG=4.259E-5 KB1=0 KB2=1.681 KB3=0 KB4=2.433 KVBGI=5.341E-7 KNK=0.06215 KNG=0.004617 KNPL=0.001662 KNSL=0.7774 KNPR=56.92 KNSR=44.23
+ CCG=12P CGP=0.2P CCP=7P VGOFF=-2.25 IGA=0.001788 IGB=0.039 IGC=0.1 IGEX=1.88
* Vp_MAX=800 Ip_MAX=500 Vg_step=5 Vg_start=30 Vg_count=17
* X_MIN=48 Y_MIN=10 X_SIZE=717 Y_SIZE=561 FSZ_X=1288 FSZ_Y=723 XYGrid=true
* Rp=1400 Vg_ac=20 P_max=25 Vg_qui=-10 Vp_qui=300 
* showLoadLine=n showIp=y isDHP=n isPP=n isAsymPP=n isUL=n showDissipLimit=y 
* showIg1=y isInputSnapped=y addLocalNFB=n
* XYProjections=n harmonicPlot=y dissipPlot=n 
* UL=0.43 EG2=250 gridLevel2=y addKink=y isTanhKnee=n advSigmoid=y 
*----------------------------------------------------------------------------------
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
E4   8 0  VALUE={V(P,K)/KNEE/(KVBGI+V(6)*KVBG)}
E5  81 0  VALUE={PWR(V(8),KB1)}
E6  82 0  VALUE={PWR(V(8),KB2)}
E7  83 0  VALUE={PWR(V(8),KB3)}
E8   9 0  VALUE={PWR(1-EXP(-V(81)*(KC+KR1*V(82))/(KD+KR2*V(83))),KB4)*1.5708}
RE4  8 0  1
RE5 81 0  1
RE6 82 0  1
RE7 83 0  1
RE8  9 0  1
RE21 21 0 1
E21  21 0 VALUE={V(6)/KG1*V(9)} ; Ip with knee but no slope and no kink
RE22 22 0 1 ; E22: kink curr deviation for plate
E22  22 0 VALUE={V(21)*LIMIT(KNK-V(G,K)*KNG,0,0.3)*(-ATAN((V(P,K)-KNPL)/KNSL)+ATAN((V(P,K)-KNPR)/KNSR))} 
G1   P K  VALUE={V(21)*(1+KLAMG*V(P,K))+KLAM*V(P,K) + V(22)}
G2   G2 K  VALUE={V(6)/KG2*(KVC-V(9))/(1+KLAMG*V(P,K)) - V(22)}
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
RE23 G 0 1G 
GG G K VALUE={(IGA+IGB/(IGC+V(P,K)))*(MU/KG1)*
+(PWR(V(G,K)-VGOFF,IGEX)+PWRS(V(G,K)-VGOFF,IGEX))}
.ENDS
*$

Koonw · 2021-02-11 6:31 am

This is Mullard EL34 model, the curve is smooth so no Knee/Slope or Advance Knee/Slope is used, so no Option is required and it runs faster in simulation.

Code:

*** EL34_MU-2 ******************************6************
* Created on 02/11/2021 13:57 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: el34-mu.png
* Data source link: <plate curves URL>
*----------------------------------------------------------------------------------
.SUBCKT EL34 P G2 G K ; LTSpice tetrode.asy pinout
* .SUBCKT EL34 P G K G2 ; Koren Pentode Pspice pinout
+ PARAMS: MU=12.04 KG1=1305.54 KP=62.44 KVB=5480.79 VCT=1.73 EX=1.494 KG2=1134.84 KNEE=14.78 KVC=1.799
+ KLAM=2E-6 KLAMG=5.51E-4   KNK=5.296E-4 KNG=0.01341 KNPL=0.2079 KNSL=90.05 KNPR=67.88 KNSR=27.4
+ CCG=3P CGP=1.4P CCP=1.9P VGOFF=-1.5 IGA=7.283E-4 IGB=0.237 IGC=7.12 IGEX=1.6
* Vp_MAX=500 Ip_MAX=350 Vg_step=5 Vg_start=0 Vg_count=8
* X_MIN=64 Y_MIN=28 X_SIZE=758 Y_SIZE=570 FSZ_X=1296 FSZ_Y=736 XYGrid=false
* Rp=1400 Vg_ac=20 P_max=25 Vg_qui=-17.5 Vp_qui=300
* showLoadLine=n showIp=y isDHP=n isPP=n isAsymPP=n isUL=n showDissipLimit=y
* showIg1=y isInputSnapped=y addLocalNFB=n
* XYProjections=n harmonicPlot=y dissipPlot=n
* UL=0.43 EG2=250 gridLevel2=y addKink=y 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
RE21 21 0 1
E21  21 0 VALUE={V(6)/KG1*ATAN(V(P,K)/KNEE)} ; Ip with knee but no slope and no kink
RE22 22 0 1 ; E22: kink curr deviation for plate
E22  22 0 VALUE={V(21)*LIMIT(KNK-V(G,K)*KNG,0,0.3)*(-ATAN((V(P,K)-KNPL)/KNSL)+ATAN((V(P,K)-KNPR)/KNSR))}
G1   P  K  VALUE={V(21)*(1+KLAMG*V(P,K))+KLAM*V(P,K) + V(22)}
* Alexander Gurskii screen current, see audioXpress 2/2011, with slope and kink added
RE43 43 K 1G ; Dummy
E43  43 G2 VALUE={0} ; Dummy
G2   43 K  VALUE={V(6)/KG2*(KVC-ATAN(V(P,K)/KNEE))/(1+KLAMG*V(P,K))-V(22)}
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
RE23 G 0 1G
GG G K VALUE={(IGA+IGB/(IGC+V(P,K)))*(MU/KG1)*
+(PWR(V(G,K)-VGOFF,IGEX)+PWRS(V(G,K)-VGOFF,IGEX))}
.ENDS
*$

Koonw · 2021-02-11 6:20 pm

This is 6L6GC RCA model being made less complex by setting Knee/Slope KB1 and KB3 to 0.

It works with the following options:
.options GminSteps=0
.options NoOpIter
.options SrcSteps=20

Code:

**** 6L6GC_RCA *************************L6*****************
* Created on 02/12/2021 01:44 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: 6L6GC-RCA.png
* Data source link: <plate curves URL>
*----------------------------------------------------------------------------------
.SUBCKT 6L6GC P G2 G K ; LTSpice tetrode.asy pinout
* .SUBCKT 6L6GC P G K G2 ; Koren Pentode Pspice pinout
+ PARAMS: MU=9.37 KG1=1057.35 KP=43.88 KVB=21148.73 VCT=-6.207 EX=1.311 KG2=1521.16 KNEE=0.05078 KVC=1.759
+ KLAM=7.263E-8 KLAMG=4.379E-4  KD=2.28E8 KC=9765.62 KR1=1.808 KR2=0.06272 KVBG=6.026E-4 KB1=0 KB2=2.23 KB3=0 KB4=1.384 KVBGI=0.03725 KNK=0.03138 KNG=0.00513 KNPL=4.198E-10 KNSL=15.93 KNPR=53.36 KNSR=65.29
+ CCG=10P CGP=0.6P CCP=6.5P VGOFF=-1.5 IGA=0.001 IGB=0.3 IGC=8 IGEX=1.54
* Vp_MAX=700 Ip_MAX=300 Vg_step=10 Vg_start=0 Vg_count=8
* X_MIN=55 Y_MIN=1 X_SIZE=778 Y_SIZE=627 FSZ_X=1296 FSZ_Y=736 XYGrid=true
* Rp=1400 Vg_ac=20 P_max=30 Vg_qui=-35 Vp_qui=300 
* showLoadLine=n showIp=y isDHP=n isPP=n isAsymPP=n isUL=n showDissipLimit=y 
* showIg1=y isInputSnapped=y addLocalNFB=n
* XYProjections=n harmonicPlot=y dissipPlot=n 
* UL=0.43 EG2=300 gridLevel2=y addKink=y isTanhKnee=n advSigmoid=y 
*----------------------------------------------------------------------------------
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
E4   8 0  VALUE={V(P,K)/KNEE/(KVBGI+V(6)*KVBG)}
E5  81 0  VALUE={PWR(V(8),KB1)}
E6  82 0  VALUE={PWR(V(8),KB2)}
E7  83 0  VALUE={PWR(V(8),KB3)}
E8   9 0  VALUE={PWR(1-EXP(-V(81)*(KC+KR1*V(82))/(KD+KR2*V(83))),KB4)*1.5708}
RE4  8 0  1
RE5 81 0  1
RE6 82 0  1
RE7 83 0  1
RE8  9 0  1
RE21 21 0 1
E21  21 0 VALUE={V(6)/KG1*V(9)} ; Ip with knee but no slope and no kink
RE22 22 0 1 ; E22: kink curr deviation for plate
E22  22 0 VALUE={V(21)*LIMIT(KNK-V(G,K)*KNG,0,0.3)*(-ATAN((V(P,K)-KNPL)/KNSL)+ATAN((V(P,K)-KNPR)/KNSR))} 
G1   P K  VALUE={V(21)*(1+KLAMG*V(P,K))+KLAM*V(P,K) + V(22)}
G2   G2 K  VALUE={V(6)/KG2*(KVC-V(9))/(1+KLAMG*V(P,K)) - V(22)}
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
RE23 G 0 1G 
GG G K VALUE={(IGA+IGB/(IGC+V(P,K)))*(MU/KG1)*
+(PWR(V(G,K)-VGOFF,IGEX)+PWRS(V(G,K)-VGOFF,IGEX))}
.ENDS
*$

Koonw · 2021-02-13 12:44 pm

Refer to Vacuum Tube SPICE Models
This PL519 Philips has been remodeled again, it runs with following options:
.options GminSteps=0
.options NoOpIter
.options SrcSteps=20

Code:

**** PL519_P ************************L1-3******************
* Created on 02/13/2021 20:13 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: pl519_p.png
* Data source link: <plate curves URL>
*----------------------------------------------------------------------------------
.SUBCKT PL519 P G2 G K ; LTSpice tetrode.asy pinout
* .SUBCKT PL519 P G K G2 ; Koren Pentode Pspice pinout
+ PARAMS: MU=3.486 KG1=26089.53 KP=22.1 KVB=0 VCT=7.685 EX=2.3 KG2=34084.71 KNEE=6.86 KVC=1.644
+ KLAMG=5.702E-4  KD=0.02666 KC=9.92 KR1=0.001573 KR2=0.1656 KVBG=0.005306 KB1=1.027 KB2=1 KB3=0.0756 KB4=2.387 KVBGI=11.84 KNK=7.629E-7 KNG=0.002511 KNPL=8.043E-4 KNSL=8.257 KNPR=181.28 KNSR=55.94
+ CCG=0.2P CGP=3P CCP=1P VGOFF=-1 IGA=0.001 IGB=0.3 IGC=7.44 IGEX=1.3
* Vp_MAX=500 Ip_MAX=2000 Vg_step=10 Vg_start=0 Vg_count=22
* X_MIN=68 Y_MIN=51 X_SIZE=736 Y_SIZE=598 FSZ_X=1296 FSZ_Y=736 XYGrid=false
* Rp=1400 Vg_ac=20 P_max=45 Vg_qui=-105 Vp_qui=300 
* showLoadLine=n showIp=y isDHP=n isPP=n isAsymPP=n isUL=n showDissipLimit=y 
* showIg1=y isInputSnapped=y addLocalNFB=n
* XYProjections=n harmonicPlot=y dissipPlot=n 
* UL=0.43 EG2=190 gridLevel2=y addKink=y isTanhKnee=n advSigmoid=y 
*----------------------------------------------------------------------------------
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
E4   8 0  VALUE={V(P,K)/KNEE/(KVBGI+V(6)*KVBG)}
E5  81 0  VALUE={PWR(V(8),KB1)}
E6  82 0  VALUE={PWR(V(8),KB2)}
E7  83 0  VALUE={PWR(V(8),KB3)}
E8   9 0  VALUE={PWR(1-EXP(-V(81)*(KC+KR1*V(82))/(KD+KR2*V(83))),KB4)*1.5708}
RE4  8 0  1
RE5 81 0  1
RE6 82 0  1
RE7 83 0  1
RE8  9 0  1
RE21 21 0 1
E21  21 0 VALUE={V(6)/KG1*V(9)} ; Ip with knee but no slope and no kink
RE22 22 0 1 ; E22: kink curr deviation for plate
E22  22 0 VALUE={V(21)*LIMIT(KNK-V(G,K)*KNG,0,0.3)*(-ATAN((V(P,K)-KNPL)/KNSL)+ATAN((V(P,K)-KNPR)/KNSR))} 
G1   P K  VALUE={V(21)*(1+KLAMG*V(P,K)) + V(22)}
G2   G2 K  VALUE={V(6)/KG2*(KVC-V(9))/(1+KLAMG*V(P,K)) - V(22)}
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
RE23 G 0 1G 
GG G K VALUE={(IGA+IGB/(IGC+V(P,K)))*(MU/KG1)*
+(PWR(V(G,K)-VGOFF,IGEX)+PWRS(V(G,K)-VGOFF,IGEX))}
.ENDS
*$

Koonw · 2021-02-13 12:52 pm

This EL509 JJ model, it runs with the following options:
.options GminSteps=0
.options NoOpIter
.options SrcSteps=20

Code:

**** EL509_JJ ************************L1-2******************
* Created on 02/13/2021 17:55 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: el509_jj.png
* Data source link: <plate curves URL>
*----------------------------------------------------------------------------------
.SUBCKT EL509 P G2 G K ; LTSpice tetrode.asy pinout
* .SUBCKT EL509 P G K G2 ; Koren Pentode Pspice pinout
+ PARAMS: MU=3.591 KG1=3431.47 KP=27.42 KVB=131.34 VCT=4.523 EX=1.874 KG2=5324.45 KNEE=2.478 KVC=1.587
+ KLAM=2.57E-13 KLAMG=4.199E-4 KNEE2=27.32 KNEX=0.596  KNK=0.02886 KNG=0.01649 KNPL=20.67 KNSL=31.95 KNPR=50.92 KNSR=145.97
+ CCG=25P CGP=2.5P CCP=18.5P VGOFF=-1 IGA=0.001 IGB=0.3 IGC=8 IGEX=1.54
* Vp_MAX=400 Ip_MAX=1600 Vg_step=10 Vg_start=0 Vg_count=15
* X_MIN=61 Y_MIN=78 X_SIZE=696 Y_SIZE=523 FSZ_X=1296 FSZ_Y=736 XYGrid=false
* Rp=1400 Vg_ac=20 P_max=35 Vg_qui=-70 Vp_qui=300 
* showLoadLine=n showIp=y isDHP=n isPP=n isAsymPP=n isUL=n showDissipLimit=y 
* showIg1=y isInputSnapped=y addLocalNFB=n
* XYProjections=n harmonicPlot=y dissipPlot=n 
* UL=0.43 EG2=160 gridLevel2=y addKink=y isTanhKnee=y 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
RE21 21 0 1
E21  21 0 VALUE={V(6)/KG1*ATAN((V(P,K)+KNEX)/KNEE)*TANH(V(P,K)/KNEE2)} ; Ip with knee but no slope and no kink
RE22 22 0 1 ; E22: kink curr deviation for plate
E22  22 0 VALUE={V(21)*LIMIT(KNK-V(G,K)*KNG,0,0.3)*(-ATAN((V(P,K)-KNPL)/KNSL)+ATAN((V(P,K)-KNPR)/KNSR))} 
G1   P  K  VALUE={V(21)*(1+KLAMG*V(P,K))+KLAM*V(P,K) + V(22)}
* Alexander Gurskii screen current, see audioXpress 2/2011, with slope and kink added
RE43 43 K 1G ; Dummy
E43  43 G2 VALUE={0} ; Dummy
G2   43 K  VALUE={V(6)/KG2*(KVC-ATAN((V(P,K)+KNEX)/KNEE)*TANH(V(P,K)/KNEE2))/(1+KLAMG*V(P,K))-V(22)}
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
RE23 G 0 1G 
GG G K VALUE={(IGA+IGB/(IGC+V(P,K)))*(MU/KG1)*
+(PWR(V(G,K)-VGOFF,IGEX)+PWRS(V(G,K)-VGOFF,IGEX))}
.ENDS
*$

Adrian Immler · 2021-02-14 7:16 am

Adrian Immler said:
You 're welcome, rongon!
BTW, you asked for a 6N6P model in the post #2738. Some samples just arrived! (...)

... and here it is, the freshly baked 6N6P model!

Based on the most representative triode out of 4 tubes, all burnt-in for 100 hours.

next tube will be the 6N15P, but it will take a bit more time as I have first to buy some 7pin sockets for the burnin.

cheers, Adrian

Code:

*6N6P LTspice model based on the generic triode model from Adrian Immler, version i4
*A version log is at the end of this file
*100h BurnIn of 4 Novosibirsk factory tubes, sample selection and measurements done in Febr. 2021
*Params fitted to the measured values by Adrian Immler, Febr. 2021
*The high fit quality is presented at adrianimmler.simplesite.com
*History's best of tube decribing art (plus some new ideas) is merged to this new approach.
*@ neg. Vg, Ia accuracy is similar to Koren models.
*@ small neg. Vg, the "Anlauf" current is considered.
*@ pos. Vg, Ig and Ia accuracy is on a unrivaled level.
*This offers new simulation possibilities like bias point setting with MOhm grid resistor,
*Audion radio circuits, low voltage amps, guitar distortion stages or pulsed stages.
*               anode (plate)
*               | grid
*               | | cathode
*               | | |
.subckt 6N6P.i4 A G K
+ params:
*Parameters for the space charge current @ Vg <= 0
+ mu   = 20.3 ;Determines the voltage gain @ constant Ia
+ rad  = 1k1  ;Differential anode resistance, set @ Iad and Vg=0V
+ Vct  = -0.85;Offsets the Ia-traces on the Va axis. Electrode material's contact potential
+ kp   = 90   ;Mimics the island effect
+ xs   = 1.60 ;Determines the curve of the Ia traces. Typically between 1.2 and 1.8
*
*Parameters for assigning the space charge current to Ia and Ig @ Vg > 0
+ kB   = 0.4  ;Describes how fast Ia drops to zero when Va approaches zero.
+ radl = 10   ;Differential resistance for the Ia emission limit @ very small Va and Vg > 0
+ tsh  = 5    ;Ia transmission sharpness from 1th to 2nd Ia area. Keep between 3 and 20. Start with 20.
+ xl   = 1.2  ;Exponent for the emission limit
*
*Parameters of the grid-cathode vacuum diode
+ kvdg  = 27  ;virtual vacuumdiode. Causes an Ia reduction @ Ig > 0.
+ kg = 590    ;Inverse scaling factor for the Va independent part of Ig (caution - interacts with xg!)
+ Vctg = -0.75;Offsets the log Ig-traces on the Vg axis. Electrode material's contact potential
+ xg   = 1.4  ;Determines the curve of the Ig slope versus (positive) Vg and Va >> 0
+ VT   = 0.11 ;Log(Ig) slope @ Vg<0. VT=k/q*Tk (cathodes absolute temp, typically 1150K)
+ kVT=0.1     ;Va dependant koeff. of VT
+ Vft2 = 0.0 gft2 = 20 ;finetunes the gridcurrent @ low Va and Vg near zero
*
*Parameters for the caps
+ cag  = 3p5  ;From datasheet
+ cak  = 1p85 ;From datasheet
+ cgk  = 4p4  ;From datasheet
*
*special purpose parameters
+ os = 1      ;Overall scaling factor, if a user wishes to simulate manufacturing tolerances
*
*Calculated parameters
+ Iad = {100/rad} ;Ia where the anode a.c. resistance is set according to rad.
+ ks = {pow(mu/(rad*xs*Iad**(1-1/xs)),-xs)} ;Reduces the unwished xs influence to the Ia slope
+ ksnom = {pow(mu/(rad*1.5*Iad**(1-1/1.5)),-1.5)} ;Sub-equation for calculating Vg0
+ Vg0 = {Vct + (Iad*ks)**(1/xs) - (Iad*ksnom)**(2/3)} ;Reduces the xs influence to Vct.
+ kl = {pow(1/(radl*xl*Ild**(1-1/xl)),-xl)} ;Reduces the xl influence to the Ia slope @ small Va
+ Ild = {sqrt(radl)*1m} ;Current where the Il a.c. resistance is set according to radl.
*
*Space charge current model
Bggi GGi 0 V=v(Gi,K)+Vg0 ;Effective internal grid voltage.
Bahc Ahc 0 V=uramp(v(A,K)) ;Anode voltage, hard cut to zero @ neg. value
Bst   St 0 V=uramp(max(v(GGi)+v(A,K)/(mu), v(A,K)/kp*ln(1+exp(kp*(1/mu+v(GGi)/(1+v(Ahc)))))));Steering volt.
Bs    Ai K I=os/ks*pow(v(St),xs) ;Langmuir-Childs law for the space charge current Is
*
*Anode current limit @ small Va
.func smin(z,y,k) {pow(pow(z+1f, -k)+pow(y+1f, -k), -1/k)} ;Min-function with smooth trans.
Ra  A Ai 1
Bgl Gi A I=min(i(Ra)-smin(1/kl*pow(v(Ahc),xl),i(Ra),tsh),i(Bgvd)*exp(4*v(G,K))) ;Ia emission limit
*
*Grid model
Bvdg G Gi I=1/kvdg*pwrs(v(G,Gi),1.5) ;Reduces the internal effective grid voltage when Ig rises
Rgip G Gi 1G ;avoids some warnings
Cvdg G Gi 0p1;this small cap improves convergence
.func fVT() {VT*exp(-kVT*sqrt(v(A,K)))}
.func Ivd(Vvd, kvd, xvd, VTvd)  {if(Vvd < 3, 1/kvd*pow(VTvd*xvd*ln(1+exp(Vvd/VTvd/xvd)),xvd), 1/kvd*pow(Vvd, xvd))} ;Vacuum diode function
Bgvd Gi K I=Ivd(v(G,K) + Vctg, kg/os, xg, fVT())
.func ft2() {gft2*(1-tanh(3*(v(G,K)+Vft2)))} ;Finetuning-func. improves ig-fit @ Vg near -0.5V, low Va.
Bgr Gi Ai I=ivd(v(GGi),ks/os, xs, 0.8*VT)/(1+ft2()+kB*v(Ahc));Is reflection to grid when Va approaches zero
Bs0 Ai K  I=ivd(v(GGi),ks/os, xs, 0.8*VT)/(1+ft2()) - os/ks*pow(v(GGi),xs) ;Compensates neg Ia @ small Va and Vg near zero
*
*Caps
C1 A G {cag}
C2 A K {cak}
C3 G K {cgk}
.ends
*
*Version log
*i1 :Initial version
*i2 :Pin order changed to the more common order âA G Kâ (Thanks to Markus Gyger for his tip)
*i3 :bugfix of the Ivd-function: now also usable for larger Vvd
*i4: Rgi replaced by a virtual vacuum diode (better convergence). ft1 deleted (no longer needed)
;2 new prarams for Ig finetuning @ Va and Vg near zero. New emission skaling factor ke for aging etc.

rongon · 2021-02-15 3:40 pm

Fabulous! Thank you for this and for all the work you've put in. I'll be giving this model a go later today.

Koonw · 2021-02-16 5:29 am

Here is Svetlana 4X150A/7034 model and sch if you need it.
options GminSteps=0
.options NoOpIter
.options SrcSteps=200
The Knee/Slope is sharp and is sensitive to accuracy of the sch in sim.

Code:

**** 4X150A ***********************L1-3*******************
* Created on 02/16/2021 13:01 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: 4x150a.png
* Data source link: <plate curves URL>
*----------------------------------------------------------------------------------
.SUBCKT 4X150A P G2 G K ; LTSpice tetrode.asy pinout
* .SUBCKT 4X150A P G K G2 ; Koren Pentode Pspice pinout
+ PARAMS: MU=8.28 KG1=1059.3 KP=21.48 KVB=5038.08 VCT=1.704 EX=1.4 KG2=41583.36 KNEE=20.88 KVC=3.752
+ KLAM=6.25E-9 KLAMG=8.91E-5  KD=1.231E6 KC=0.01918 KR1=0 KR2=0.1693 KVBG=0.00104 KB1=9.32 KB2=0 KB3=0 KB4=0.081 KVBGI=0.003125 KNK=0.0938 KNG=0.002097 KNPL=15.05 KNSL=37.38 KNPR=275.88 KNSR=84.22
+ CCG=15.7P CGP=0.04P CCP=4.5P VGOFF=-7.5 IGA=7.1E-4 IGB=0.084 IGC=15.52 IGEX=1.84
* Vp_MAX=1500 Ip_MAX=1400 Vg_step=20 Vg_start=20 Vg_count=18
* X_MIN=40 Y_MIN=29 X_SIZE=801 Y_SIZE=557 FSZ_X=1296 FSZ_Y=736 XYGrid=false
* Rp=1400 Vg_ac=20 P_max=250 Vg_qui=-150 Vp_qui=300 
* showLoadLine=n showIp=y isDHP=n isPP=n isAsymPP=n isUL=n showDissipLimit=y 
* showIg1=y isInputSnapped=y addLocalNFB=n
* XYProjections=n harmonicPlot=y dissipPlot=n 
* UL=0.43 EG2=300 gridLevel2=y addKink=y isTanhKnee=n advSigmoid=y 
*----------------------------------------------------------------------------------
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
E4   8 0  VALUE={V(P,K)/KNEE/(KVBGI+V(6)*KVBG)}
E5  81 0  VALUE={PWR(V(8),KB1)}
E6  82 0  VALUE={PWR(V(8),KB2)}
E7  83 0  VALUE={PWR(V(8),KB3)}
E8   9 0  VALUE={PWR(1-EXP(-V(81)*(KC+KR1*V(82))/(KD+KR2*V(83))),KB4)*1.5708}
RE4  8 0  1
RE5 81 0  1
RE6 82 0  1
RE7 83 0  1
RE8  9 0  1
RE21 21 0 1
E21  21 0 VALUE={V(6)/KG1*V(9)} ; Ip with knee but no slope and no kink
RE22 22 0 1 ; E22: kink curr deviation for plate
E22  22 0 VALUE={V(21)*LIMIT(KNK-V(G,K)*KNG,0,0.3)*(-ATAN((V(P,K)-KNPL)/KNSL)+ATAN((V(P,K)-KNPR)/KNSR))} 
G1   P K  VALUE={V(21)*(1+KLAMG*V(P,K))+KLAM*V(P,K) + V(22)}
G2   G2 K  VALUE={V(6)/KG2*(KVC-V(9))/(1+KLAMG*V(P,K)) - V(22)}
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
RE23 G 0 1G 
GG G K VALUE={(IGA+IGB/(IGC+V(P,K)))*(MU/KG1)*
+(PWR(V(G,K)-VGOFF,IGEX)+PWRS(V(G,K)-VGOFF,IGEX))}
.ENDS
*$

Koonw · 2021-02-16 6:35 am

Ia in left side scale of the above model is wrong, now corrected and it's Eimac 4-150A model.

Code:

**** 4X150A *********************L1-3*********************
* Created on 02/16/2021 14:25 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: 4x150a.png
* Data source link: <plate curves URL>
*----------------------------------------------------------------------------------
.SUBCKT 4X150A P G2 G K ; LTSpice tetrode.asy pinout
* .SUBCKT 4X150A P G K G2 ; Koren Pentode Pspice pinout
+ PARAMS: MU=8.28 KG1=741.51 KP=21.48 KVB=5038.08 VCT=1.704 EX=1.4 KG2=23761.92 KNEE=20.88 KVC=3.752
+ KLAM=6.25E-9 KLAMG=8.91E-5  KD=1.589E-4 KC=0.01676 KR1=0.002223 KR2=0.0576 KVBG=0.006022 KB1=5.84 KB2=1.34 KB3=0.73 KB4=0.162 KVBGI=0.00118 KNK=0.091 KNG=0.00233 KNPL=3.225 KNSL=46.25 KNPR=278.16 KNSR=86.77
+ CCG=15.7P CGP=0.04P CCP=4.5P VGOFF=-7.5 IGA=7.1E-4 IGB=0.084 IGC=15.52 IGEX=1.84
* Vp_MAX=1500 Ip_MAX=1400 Vg_step=20 Vg_start=20 Vg_count=18
* X_MIN=40 Y_MIN=196 X_SIZE=801 Y_SIZE=390 FSZ_X=1296 FSZ_Y=736 XYGrid=false
* Rp=1400 Vg_ac=20 P_max=250 Vg_qui=-150 Vp_qui=300 
* showLoadLine=n showIp=y isDHP=n isPP=n isAsymPP=n isUL=n showDissipLimit=y 
* showIg1=y isInputSnapped=y addLocalNFB=n
* XYProjections=n harmonicPlot=y dissipPlot=n 
* UL=0.43 EG2=300 gridLevel2=y addKink=y isTanhKnee=n advSigmoid=y 
*----------------------------------------------------------------------------------
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
E4   8 0  VALUE={V(P,K)/KNEE/(KVBGI+V(6)*KVBG)}
E5  81 0  VALUE={PWR(V(8),KB1)}
E6  82 0  VALUE={PWR(V(8),KB2)}
E7  83 0  VALUE={PWR(V(8),KB3)}
E8   9 0  VALUE={PWR(1-EXP(-V(81)*(KC+KR1*V(82))/(KD+KR2*V(83))),KB4)*1.5708}
RE4  8 0  1
RE5 81 0  1
RE6 82 0  1
RE7 83 0  1
RE8  9 0  1
RE21 21 0 1
E21  21 0 VALUE={V(6)/KG1*V(9)} ; Ip with knee but no slope and no kink
RE22 22 0 1 ; E22: kink curr deviation for plate
E22  22 0 VALUE={V(21)*LIMIT(KNK-V(G,K)*KNG,0,0.3)*(-ATAN((V(P,K)-KNPL)/KNSL)+ATAN((V(P,K)-KNPR)/KNSR))} 
G1   P K  VALUE={V(21)*(1+KLAMG*V(P,K))+KLAM*V(P,K) + V(22)}
G2   G2 K  VALUE={V(6)/KG2*(KVC-V(9))/(1+KLAMG*V(P,K)) - V(22)}
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
RE23 G 0 1G 
GG G K VALUE={(IGA+IGB/(IGC+V(P,K)))*(MU/KG1)*
+(PWR(V(G,K)-VGOFF,IGEX)+PWRS(V(G,K)-VGOFF,IGEX))}
.ENDS
*$

Koonw · 2021-02-17 11:40 am

This 4X150A RCA model, work without any options:

Code:

**** 4X150A *******************RCA-7***********************
* Created on 02/17/2021 15:16 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: 4x150a.png
* Data source link: <plate curves URL>
*----------------------------------------------------------------------------------
.SUBCKT 4X150A P G2 G K ; LTSpice tetrode.asy pinout
* .SUBCKT 4X150A P G K G2 ; Koren Pentode Pspice pinout
+ PARAMS: MU=5.913 KG1=2283.64 KP=19.38 KVB=5038.08 VCT=-5.253 EX=1.669 KG2=6032.41 KNEE=1881.54 KVC=1.954
+ KLAM=4.043E-12 KLAMG=6.77E-5 KNEE2=50.56 KNEX=4.795E7  KNK=0.1336 KNG=0.001604 KNPL=0.04744 KNSL=176.92 KNPR=243.07 KNSR=91.92
+ CCG=16P CGP=0.03P CCP=4.4P VGOFF=-7.5 IGA=7.1E-4 IGB=0.084 IGC=15.52 IGEX=1.84
* Vp_MAX=2000 Ip_MAX=1600 Vg_step=10 Vg_start=30 Vg_count=17
* X_MIN=60 Y_MIN=27 X_SIZE=718 Y_SIZE=574 FSZ_X=1296 FSZ_Y=736 XYGrid=false
* Rp=1400 Vg_ac=20 P_max=250 Vg_qui=-50 Vp_qui=300 
* showLoadLine=n showIp=y isDHP=n isPP=n isAsymPP=n isUL=n showDissipLimit=y 
* showIg1=y isInputSnapped=y addLocalNFB=n
* XYProjections=n harmonicPlot=y dissipPlot=n 
* UL=0.43 EG2=250 gridLevel2=y addKink=y isTanhKnee=y 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
RE21 21 0 1
E21  21 0 VALUE={V(6)/KG1*ATAN((V(P,K)+KNEX)/KNEE)*TANH(V(P,K)/KNEE2)} ; Ip with knee but no slope and no kink
RE22 22 0 1 ; E22: kink curr deviation for plate
E22  22 0 VALUE={V(21)*LIMIT(KNK-V(G,K)*KNG,0,0.3)*(-ATAN((V(P,K)-KNPL)/KNSL)+ATAN((V(P,K)-KNPR)/KNSR))} 
G1   P  K  VALUE={V(21)*(1+KLAMG*V(P,K))+KLAM*V(P,K) + V(22)}
* Alexander Gurskii screen current, see audioXpress 2/2011, with slope and kink added
RE43 43 K 1G ; Dummy
E43  43 G2 VALUE={0} ; Dummy
G2   43 K  VALUE={V(6)/KG2*(KVC-ATAN((V(P,K)+KNEX)/KNEE)*TANH(V(P,K)/KNEE2))/(1+KLAMG*V(P,K))-V(22)}
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
RE23 G 0 1G 
GG G K VALUE={(IGA+IGB/(IGC+V(P,K)))*(MU/KG1)*
+(PWR(V(G,K)-VGOFF,IGEX)+PWRS(V(G,K)-VGOFF,IGEX))}
.ENDS
*$

woodo · 2021-02-17 9:59 pm

Posted a query in the other spice thread. I’m just trying to model a prospective 6AS7 amp designed by @suncalc. I couldn’t find a proper audio potentiometer model so consulted the wiki and YouTube. I came up with the attached and by luck it seems to work. I don’t know how an actual audio pot is supposed to behave and would appreciate any advice to improve it and fix errors. Thanks. https://www.diyaudio.com/forums/sof...ltxvii-beginner-advanced-259.html#post6528950

rongon · 2021-02-18 12:37 am

Are the best ECC88 - 6DJ8 models around the Ayumi one and the 6DJ8 low current model by cognscogs?

Vacuum Tube SPICE Models

Member

Member

Member

Member

Member

Member

Member

Disabled Account

Attachments

Disabled Account

Attachments

Disabled Account

Attachments

Disabled Account

Attachments

Disabled Account

Attachments

Disabled Account

Attachments

Member

Attachments

Member

Disabled Account

Attachments

Disabled Account

Attachments

Disabled Account

Attachments

Member

Member