Vacuum Tube SPICE Models

euro21 · 2021-08-22 5:11 pm

Koonw said:
8185 model, try it:

Thanks!

n33d2learn · 2021-08-24 12:28 pm

Anyone have a model for the 12DT5?

Koonw · 2021-08-24 8:10 pm

12DT5 Tung Sol model:

Code:

* Created on 08/25/2021 03:00 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: 12dt5.png
* Data source link: <plate curves URL>
*----------------------------------------------------------------------------------
.SUBCKT 12DT5 P G2 G K ; LTSpice tetrode.asy pinout
* .SUBCKT 12DT5 P G K G2 ; Koren Pentode Pspice pinout
+ PARAMS: MU=11.88 KG1=1313.26 KP=41.34 KVB=12 VCT=0.2 EX=1.471 KG2=3682.33 KNEE=12.58 KVC=1.81
+ KLAM=3.125E-10 KLAMG=4.433E-7 KNEE2=16.44 KNEX=0.08206  KNK=-0.2 KNG=0.01861 KNPL=6.546 KNSL=0.06755 KNPR=124.44 KNSR=53.97
+ CCG=12.5P CGP=0.57P CCP=4.9P VGOFF=-0.6 IGA=0.001 IGB=0.3 IGC=8 IGEX=2
* Vp_MAX=500 Ip_MAX=250 Vg_step=4 Vg_start=0 Vg_count=11
* X_MIN=107 Y_MIN=15 X_SIZE=620 Y_SIZE=613 FSZ_X=1296 FSZ_Y=736 XYGrid=false
* Rp=1400 Vg_ac=20 P_max=9 Vg_qui=-20 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
*$

n33d2learn · 2021-08-24 8:45 pm

Awesome. Thank You! I couldn't get the curves to line up right, the included screen cap will be a great help if I need to model something else in the future.

OldHector · 2021-08-27 11:42 am

I am trying to model the Tubelab SP for operation with a PL84/EL86. I had to remodel the power supply so that the driver/splitter B+ is closer to the original design, so hence the two B+ supplies in the schematic (275v + 216v).

I have done something stupid and I just cannot get the model to work, at least I do not see any output at the loudspeaker.

If someone could take a quick look and point me in the right direction, I'd appreciate that!

Koonw · 2021-08-27 2:22 pm

Got an error "Couldn't find symbol(s) Penthode", try to use Tetrode symbol instead for 6cw5, it works.

OldHector · 2021-08-27 4:08 pm

Great. Thanks!

Koonw · 2021-08-30 3:00 pm

vt4c tube datasheet - Google Search

211-VT4C model:

Code:

* [url=http://www.dmitrynizh.com/tubeparams_image.htm]Model Paint Tools: Trace Tube Parameters over Plate Curves, Interactively[/url]
* Plate Curves image file: 211-GE-VT4C.png
* Data source link: 
*----------------------------------------------------------------------------------
.SUBCKT 211_GE 1 2 3 ; Plate Grid Cathode
+ PARAMS: CCG=6P  CGP=14P CCP=5.5P 
+ MU=12 KG1=3928.74 KP=451.14 KVB=249.72 VCT=-2.52 EX=1.458 
+ VGOFF=-2.62 IGA=8.906E-6 IGB=0.04944 IGC=2.4 IGEX=2 
* Vp_MAX=2200 Ip_MAX=160 Vg_step=20 Vg_start=40 Vg_count=15
* Rp=4000 Vg_ac=55 P_max=75 Vg_qui=-48 Vp_qui=300
* X_MIN=90 Y_MIN=40 X_SIZE=765 Y_SIZE=563 FSZ_X=1296 FSZ_Y=736 XYGrid=false
* showLoadLine=n showIp=y isDHT=n isPP=n isAsymPP=n showDissipLimit=y 
* showIg1=y gridLevel2=y 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 
RE2 2 0 1G 
EGC 8 0 VALUE={V(2,3)-VGOFF} ; POSITIVE GRID THRESHOLD
GG 2 3 VALUE={(IGA+IGB/(IGC+V(1,3)))*(MU/KG1)*(PWR(V(8),IGEX)+PWRS(V(8),IGEX))}
.ENDS 
*$

Koonw · 2021-08-30 9:42 pm

http://www.china-hifi-audio.com/images/PSVANE-tube/PSVANE-211-TII-vacuum-tube-007.jpg
Psvane 211-HiFi Series
Characteristics are quite close to 211 RCA curve.

Code:

PSVANE 211-TII model:
**** 211_TII ** Advanced Grid Current **********************************
* Created on 08/31/2021 04:28 using paint_kit.jar 3.1
* [url=http://www.dmitrynizh.com/tubeparams_image.htm]Model Paint Tools: Trace Tube Parameters over Plate Curves, Interactively[/url]
* Plate Curves image file: 211-TII.png
* Data source link: 
*----------------------------------------------------------------------------------
.SUBCKT 211_TII 1 2 3 ; Plate Grid Cathode
+ PARAMS: CCG=6P  CGP=14.5P CCP=5.5P 
+ MU=11.64 KG1=2208.15 KP=589.39 KVB=889.01 VCT=-3.98 EX=1.346 
+ VGOFF=-4 IGA=8.906E-6 IGB=0.04944 IGC=2.4 IGEX=2 
* Vp_MAX=1800 Ip_MAX=400 Vg_step=20 Vg_start=20 Vg_count=15
* Rp=4000 Vg_ac=55 P_max=75 Vg_qui=-48 Vp_qui=300
* X_MIN=56 Y_MIN=16 X_SIZE=787 Y_SIZE=475 FSZ_X=1296 FSZ_Y=736 XYGrid=false
* showLoadLine=n showIp=y isDHT=n isPP=n isAsymPP=n showDissipLimit=y 
* showIg1=y gridLevel2=y 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 
RE2 2 0 1G 
EGC 8 0 VALUE={V(2,3)-VGOFF} ; POSITIVE GRID THRESHOLD
GG 2 3 VALUE={(IGA+IGB/(IGC+V(1,3)))*(MU/KG1)*(PWR(V(8),IGEX)+PWRS(V(8),IGEX))}
.ENDS 
*$

Koonw · 2021-08-30 9:44 pm

211 RCA model:

Code:

* Created on 08/31/2021 04:50 using paint_kit.jar 3.1
* [URL="http://www.dmitrynizh.com/tubeparams_image.htm"]Model Paint Tools: Trace Tube Parameters over Plate Curves, Interactively[/URL]
* Plate Curves image file: 211-RCA.png
* Data source link: 
*----------------------------------------------------------------------------------
.SUBCKT 211_RCA 1 2 3 ; Plate Grid Cathode
+ PARAMS: CCG=5.4P  CGP=14P CCP=4.8P 
+ MU=11.64 KG1=1397.72 KP=3536.33 KVB=4994.4 VCT=-7.995 EX=1.27 
+ VGOFF=-12 IGA=8.906E-6 IGB=0.04944 IGC=2.4 IGEX=2 
* Vp_MAX=2000 Ip_MAX=1600 Vg_step=25 Vg_start=200 Vg_count=15
* Rp=4000 Vg_ac=55 P_max=75 Vg_qui=-48 Vp_qui=300
* X_MIN=38 Y_MIN=8 X_SIZE=802 Y_SIZE=647 FSZ_X=1296 FSZ_Y=736 XYGrid=false
* showLoadLine=n showIp=y isDHT=n isPP=n isAsymPP=n showDissipLimit=y 
* showIg1=y gridLevel2=y 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 
RE2 2 0 1G 
EGC 8 0 VALUE={V(2,3)-VGOFF} ; POSITIVE GRID THRESHOLD
GG 2 3 VALUE={(IGA+IGB/(IGC+V(1,3)))*(MU/KG1)*(PWR(V(8),IGEX)+PWRS(V(8),IGEX))}
.ENDS 
*$

prazza1 · 2021-08-31 6:47 am

Hello,

Does the Spice model for the DHT Triode 31 exist ?

Thanks.

Koonw · 2021-08-31 1:57 pm

31 RCA model, try it:

Code:

* Created on 08/31/2021 20:43 using paint_kit.jar 3.1
* [url=http://www.dmitrynizh.com/tubeparams_image.htm]Model Paint Tools: Trace Tube Parameters over Plate Curves, Interactively[/url]
* Plate Curves image file: 31.png
* Data source link: 
*----------------------------------------------------------------------------------
.SUBCKT 31 1 2 3 ; Plate Grid Cathode
+ PARAMS: CCG=3.5P  CGP=5.7P CCP=2.7P 
+ MU=3.8 KG1=9180 KP=51.2 KVB=210 VCT=-0.2175 EX=1.414 
+ VGOFF=-0.6 IGA=0.001 IGB=0.3 IGC=8 IGEX=1.6 
* Vp_MAX=250 Ip_MAX=30 Vg_step=5 Vg_start=0 Vg_count=20
* Rp=4000 Vg_ac=55 P_max=3 Vg_qui=-48 Vp_qui=300
* X_MIN=49 Y_MIN=35 X_SIZE=763 Y_SIZE=571 FSZ_X=1296 FSZ_Y=736 XYGrid=false
* showLoadLine=n showIp=y isDHT=n isPP=n isAsymPP=n showDissipLimit=y 
* showIg1=y gridLevel2=y 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 
RE2 2 0 1G 
EGC 8 0 VALUE={V(2,3)-VGOFF} ; POSITIVE GRID THRESHOLD
GG 2 3 VALUE={(IGA+IGB/(IGC+V(1,3)))*(MU/KG1)*(PWR(V(8),IGEX)+PWRS(V(8),IGEX))}
.ENDS 
*$

prazza1 · 2021-08-31 6:20 pm

Thank you koonw !

Koonw · 2021-09-01 4:19 am

Paint Tool can generate DHT model for A1 or A2. 31 has 2V and 130mA filament so Rfil=2/0.13=15.4 is used. Please try if any better:

Code:

* Created on 09/01/2021 11:05 using paint_kit.jar 3.1
* [url=http://www.dmitrynizh.com/tubeparams_image.htm]Model Paint Tools: Trace Tube Parameters over Plate Curves, Interactively[/url]
* Plate Curves image file: 31.png
* Data source link: 
*----------------------------------------------------------------------------------
.SUBCKT DHT_31_A2 1 2 3 4 ; P G K1 K2
+ PARAMS: CCG=3.5P  CGP=5.7P CCP=2.7P RFIL=15.4
+ MU=3.8 KG1=9180 KP=51.2 KVB=210 VCT=-0.2175 EX=1.414 
+ VGOFF=-0.6 IGA=0.001 IGB=0.3 IGC=8 IGEX=1.6 
* Vp_MAX=250 Ip_MAX=30 Vg_step=5 Vg_start=20 Vg_count=20
* Rp=4000 Vg_ac=55 P_max=3 Vg_qui=-48 Vp_qui=300
* X_MIN=49 Y_MIN=35 X_SIZE=763 Y_SIZE=571 FSZ_X=1296 FSZ_Y=736 XYGrid=false
* showLoadLine=n showIp=y isDHT=y isPP=n isAsymPP=n showDissipLimit=y 
* showIg1=y gridLevel2=y isInputSnapped=n 
* XYProjections=n harmonicPlot=n dissipPlot=n 
*----------------------------------------------------------------------------------
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 
*$

Koonw · 2021-09-05 12:29 pm

6E5P-I (6Э5П-И) Russian military pulse tetrode tubes model:

**** 6E5P ******************************************

Code:

* Created on 09/05/2021 19:20 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: 6e5p.png
* Data source link: <plate curves URL>
*----------------------------------------------------------------------------------
.SUBCKT 6E5P P G2 G K ; LTSpice tetrode.asy pinout
* .SUBCKT 6E5P P G K G2 ; Koren Pentode Pspice pinout
+ PARAMS: MU=34.61 KG1=417.43 KP=2365.08 KVB=43.39 VCT=0.001396 EX=1.777 KG2=1270.08 KNEE=13.09 KVC=2.792
+ KLAM=2E-6 KLAMG=8.015E-4  KD=173704.09 KC=110716.95 KR1=1.33 KR2=0.24 KVBG=0.2139 KB1=7.179 KB2=0 KB3=0 KB4=0.1965 KVBGI=-0.46 KNK=0.00178 KNG=0.03673 KNPL=0.6205 KNSL=0.002871 KNPR=47.52 KNSR=25.06
+ CCG=15P CGP=0.065P CCP=2.55P VGOFF=-0.6 IGA=0.001 IGB=0.3 IGC=7.6 IGEX=1.6
* Vp_MAX=240 Ip_MAX=130 Vg_step=0.5 Vg_start=0 Vg_count=9
* X_MIN=38 Y_MIN=25 X_SIZE=829 Y_SIZE=570 FSZ_X=1296 FSZ_Y=736 XYGrid=false
* Rp=1400 Vg_ac=20 P_max=8.3 Vg_qui=-2 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=150 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
*$

PIL · 2021-09-06 9:14 am

Hello, I am bringing old tubes back to life
I couldn't find a model for the 6B7 tube
Can you help me?

Koonw · 2021-09-06 11:45 am

6B7 RCA pentode section model:
Screen current is derived from transconductance curve of 2mA max.

Code:

**** 6B7 ******************************************
* Created on 09/06/2021 18:36 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: 6b7.png
* Data source link: <plate curves URL>
*----------------------------------------------------------------------------------
.SUBCKT 6B7 P G2 G K ; LTSpice tetrode.asy pinout
* .SUBCKT 6B7 P G K G2 ; Koren Pentode Pspice pinout
+ PARAMS: MU=10 KG1=8962.8 KP=46.82 KVB=15.24 VCT=0 EX=1.4 KG2=53760 KNEE=14.79 KVC=1.799
+ KLAM=1.588E-8 KLAMG=3.186E-4  KD=0.44 KC=0.14 KR1=0.0013 KR2=0.064 KVBG=0.015 KNK=-0.044 KNG=0.006 KNPL=50 KNSL=11 KNPR=120 KNSR=29
+ CCG=3.5P CGP=0.007P CCP=9.5P VGOFF=-0.6 IGA=0.001 IGB=0.3 IGC=8 IGEX=1.54
* Vp_MAX=400 Ip_MAX=10 Vg_step=3 Vg_start=0 Vg_count=20
* X_MIN=30 Y_MIN=14 X_SIZE=800 Y_SIZE=502 FSZ_X=1296 FSZ_Y=736 XYGrid=false
* Rp=1400 Vg_ac=20 P_max=2.25 Vg_qui=-28.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=100 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={PWR(V(1,3)/KNEE/(V(6)*KVBG),2)}
E5   9 0  VALUE={SQRT(1-EXP(-V(8)*(KC+KR1*V(8))/(KD+KR2*V(8))))*1.5708}
RE4  8 0  1
RE5  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-09-06 4:03 pm

EL84 Mullard model:

Code:

**** EL84_ML ******************************************
* Created on 09/06/2021 22:58 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: el84-ML.png
* Data source link: <plate curves URL>
*----------------------------------------------------------------------------------
.SUBCKT EL84_ML P G2 G K ; LTSpice tetrode.asy pinout
* .SUBCKT EL84_ML P G K G2 ; Koren Pentode Pspice pinout
+ PARAMS: MU=20.34 KG1=548.3 KP=141.5 KVB=153.6 VCT=0.2 EX=1.316 KG2=1731.12 KNEE=36.88 KVC=1.814
+ KLAM=5.06E-9 KLAMG=3.815E-7 KNEE2=1.438 KNEX=3.534  KNK=-0.01791 KNG=0.01267 KNPL=0.00988 KNSL=0.6384 KNPR=14.52 KNSR=113.68
+ CCG=10.8P CGP=0.5P CCP=6.5P VGOFF=-0.6 IGA=0.001 IGB=0.3 IGC=8 IGEX=1.32
* Vp_MAX=400 Ip_MAX=175 Vg_step=2 Vg_start=0 Vg_count=10
* X_MIN=74 Y_MIN=29 X_SIZE=620 Y_SIZE=544 FSZ_X=1296 FSZ_Y=736 XYGrid=false
* Rp=1400 Vg_ac=20 P_max=12 Vg_qui=-9 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=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
*$

PIL · 2021-09-06 4:23 pm

Thank you very much Koonw, I will test the model today

PIL

PIL · 2021-09-06 4:52 pm

Hello Koonw,
I get the following error message during the curve reading test: "Fatal Error: u1:e4: Unknown circuit node: "u1:1" requested in behavioral source"

Vacuum Tube SPICE Models

Member

Member

Disabled Account

Attachments

Member

Member

Attachments

Disabled Account

Member

Disabled Account

Attachments

Disabled Account

Attachments

Disabled Account

Attachments

Member

Disabled Account

Attachments

Member

Disabled Account

Attachments

Disabled Account

Attachments

Member

Attachments

Disabled Account

Attachments

Disabled Account

Attachments

Member

Member