Symulacja SPICE zweryfikowana ze zmontowanym układem na tym samy egzemplarzu lampy dała zgodność widma THD w dokładnością lepszą od 1.5%.
Kod: Zaznacz cały
**** 6J45B_V_30 ******************************************
* Created on 06/06/2026 09:33 using paint_kip.jar
* www.dmitrynizh.com/tubeparams_image.htm
* Plate Curves image file: 6J45B-V_30.bmp
* Data source link: <plate curves URL>
*----------------------------------------------------------------------------------
.SUBCKT 6J45B_V_30 P G2 G K ; LTSpice tetrode.asy pinout
* .SUBCKT 6J45B_V_30 P G K G2 ; Koren Pentode Pspice pinout
+ PARAMS: MU=25.6 KG1=616 KP=100.61 KVB=0.04172 VCT=5.391E-4 EX=1.288 KG2=609 KNEE=1.35 KVC=1.799
+ KLAM=4.88E-6 KLAMG=2.88E-4 KD=1.136 KC=0.002538 KR1=6.825E-4 KR2=0.1869 KVBG=0.0381 KB1=1.94 KB2=2.06 KB3=2.52 KB4=2.28 KVBGI=6.1E-4 KNK=-0.044 KNG=0.006 KNPL=50 KNSL=11 KNPR=120 KNSR=29
+ CCG=6.1P CGP=0.05P CCP=2.1P RGI=2000.0
* Vp_MAX=50 Ip_MAX=8 Vg_step=0.25 Vg_start=0 Vg_count=6
* X_MIN=191 Y_MIN=133 X_SIZE=947 Y_SIZE=663 FSZ_X=1903 FSZ_Y=1024 XYGrid=false
* Rp=1400 Vg_ac=20 P_max=25 Vg_qui=-0.625 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=30 gridLevel2=n 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
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
*$