Ch-ki i modele nietypowych lamp.

[STUDI_bis]

Podrzucam model lampy 6J45B-V estymowany ręcznie za pomocą paint_kip dla napięcia Ug2 - 30V.
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
*$

6J45-V_30V_spice_model_fitting.png

[STUDI_bis]

6Ż52P. (Doskonale znana lampa szczególnie tym co serwisowali Rubiny )

Ug2=30V

6J52P_30.png

Ug2=60V

6J52P_60.png

Ug2=90V

6J52P_90.png

Ug2=120V

6J52P_120.png

Ug2=150V

6J52P_150.png

Ug2=180V

6J52P_180.png

Spakowane razem pliki danych pomiarowych dla powyższych zestawów charakterystyk:

6J52P_utd-files.ZIP

I jeszcze region, który minie interesuje, czyli niskonapięciowy. Ug2 = 20V, Ua <= 80V.
Dopasowanie modelu do charakterystyk za pomica Paint_KIP.jar

ModelFitting.png

Zestaw danych z jednej serii pomiarowej

6J52P_20-utd.zip

Estymowany model:

Kod: Zaznacz cały

**** 6J52P_20 *******************************************
* Created on 06/13/2026 16:15 using paint_kip.jar
* www.dmitrynizh.com/tubeparams_image.htm
* Plate Curves image file: 6J52P_20.bmp
* Data source link: <plate curves URL>
*----------------------------------------------------------------------------------
.SUBCKT 6J52P_20 P G2 G K ; LTSpice tetrode.asy pinout
* .SUBCKT 6J52P_20B P G K G2 ; Koren Pentode Pspice pinout
+ PARAMS: MU=17.54 KG1=616 KP=141.15 KVB=0.4575 VCT=0.00525 EX=2.576 KG2=945 KNEE=0.7087 KVC=1.773
+ KLAM=2.187E-8 KLAMG=1.552E-4  KD=2.32E8 KC=0.0588 KR1=0.001287 KR2=0.04224 KVBG=0.0153 KB1=2.18 KB2=1.8 KB3=2.96 KB4=0.98 KVBGI=1.65E-4 KNK=-0.044 KNG=0.006 
+ CCG=13.5P CGP=0.05P CCP=1.8P RGI=2000.0
* Vp_MAX=80 Ip_MAX=7 Vg_step=0.15 Vg_start=0 Vg_count=6
* X_MIN=192 Y_MIN=133 X_SIZE=945 Y_SIZE=661 FSZ_X=1691 FSZ_Y=1022 XYGrid=false
* Rp=1400 Vg_ac=20 P_max=25 Vg_qui=-0.375 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=20 gridLevel2=n addKink=n 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
G1   P K  VALUE={V(6)/KG1*V(9)*(1+KLAMG*V(P,K))+KLAM*V(P,K)}
G2   G2 K  VALUE={V(6)/KG2*(KVC-V(9))/(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))

**** 6J52P_20B ******************************************
* Created on 06/13/2026 16:15 using paint_kit.jar 4.7 
* www.dmitrynizh.com/tubeparams_image.htm
* Plate Curves image file: 6J52P_20.bmp
* Data source link: <plate curves URL>
*----------------------------------------------------------------------------------
.SUBCKT TRIODE_6J52P_20 1 2 3 ; Plate Grid Cathode
+ PARAMS: CCG=13.5P  CGP=2.4P CCP=1.8P RGI=2000
+ MU=17.54 KG1=1368,62 KP=141.15 KVB=0.4575 VCT=0.00525 EX=2.576 
* Vp_MAX=80 Ip_MAX=7 Vg_step=0.15 Vg_start=0 Vg_count=6
* Rp=1400 Vg_ac=20 P_max=25 Vg_qui=-0.375 Vp_qui=300
* X_MIN=192 Y_MIN=133 X_SIZE=945 Y_SIZE=661 FSZ_X=1691 FSZ_Y=1022 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 
*$

Uwaga: Dla drugiego modelu, lampy połączonej jak trioda brak jest podanej wartości pojemności zwrotnej (Cag1), założyłem że jest ona o 50% większa od pojemności wyjściowej (Cak) równej 1.8pF. Dla modelu pentody wartości pojemności pobrane z karty katalogowej https://frank.pocnet.net/sheets/112/6/6Z52P.pdf