Dmitry Nizh kindly worked out the 4P1L triode SPICE model using his great tool:
Here is Dimtry’s model:
** 4P1L_TRIODE ************************************************************
* Created on Mon Jan 07 07:31:48 PST 2013 using tube.model.finder.PaintKIT
* URL: http://www.bartola.co.uk/valves/valve-curves/4p1l/
*————————————————–
.SUBCKT TRIODE_4P1L_TRIODE 1 2 3 ; P G K ;
+ PARAMS: CCG=8P CGP=7P CCP=9P RGI=2000
+ MU=8.232 EX=1.3719 KG1=851.25 KP=108.0 KVB=528.0 VCT=-1.0 ; Vp_MAX=450.0 Ip_MAX=0.08 Vg_step=5.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
*$
