• WARNING: Tube/Valve amplifiers use potentially LETHAL HIGH VOLTAGES.
    Building, troubleshooting and testing of these amplifiers should only be
    performed by someone who is thoroughly familiar with
    the safety precautions around high voltages.

Koren's spice models

Status
This old topic is closed. If you want to reopen this topic, contact a moderator using the "Report Post" button.
A year ago i spent a lot of time fitting with originlab some precise measurements and found that "standard" tube models didn't give the best fits.

I got an almost perfect fit with some 12EA6 and a complicated function that had nothing to do with mu and other tube parameters, but i dropped it when i noticed that there were big changes when the tube had been warmed for a long time (3 hours or so).

Anything on this?
 
Sorry, fscarpa but this does not work on WinSpice
Code:
* WinSPICE Auto run
*#destroy all
*#run
*#fourier 1000 v(5)
*#linearize 
*#spec 0 20000 50 v(5)
*#let voutdb=db(v(5))
*#plot voutdb
*#setplot op1
*#print v(5)*va2#branch
*#print (v(9)-v(5))*va2#branch
*#print (v(8)-v(3))*va1#branch
*#print (v(3)-v(4))*va1#branch
*#rusage everything

c1  in  g1       3.3uF
x1  a1  g1 k1    k845rca  
rk  k1   0       1

rg  g1 b1  20k


x2  p1 a1  ou 0 LL1620_90mA  


va p1 0        1100v
vp b1 0        -130v

vin in 0 dc 0 ac 100 sin(0 100 1000)
*vin 1 0 dc 0 ac 0.5 sin(0 1.0 10000)
* vin 1 0 dc 0 ac 0
rl ou 0 8
*vout 10 0 dc 0 ac 1

.op 

*.print dc V(3) v(5)
*.print dc I(rk)
*.print dc I(rs)

.ac dec 10 5 1E6 
*.PLOT AC v(2) v(3)

.tran 5000NS 20000US 0uS
*.tran .2n 2000n 0 
*.FOUR 10000 v(3) v(5)


.SUBCKT LL1620_90mA  1 2 3 4
ra1  1   5     308
la1  5   2     40H
ka1  la1 la2   1
la2  6   4     0.0277H 
ra2  3   6     0.2
.ENDS

.SUBCKT k845rca  A G K
+PARAMS: MU=5.355 ERP=1.5
+ KK1=6323  KP=85.64  KVB=65.8 vg0=3 va0=0
+ CGA=13.5E-12 CGK=6E-12 CAK=6.5E-12 RGI=1000;

E6 6 0 VALUE={KP*( (1/MU)+((V(G,K)-vg0)/sqrt(V(A,K)**2+KVB**2)) )}
E8 8 0 VALUE={(V(A,K))/KP*LN(1+EXP(V(6)))}


*.param V_6={KP*( (1/MU)+((V(G,K)-vg0)/sqrt(V(A,K)**2+KVB**2)) )}
* E8    8  0  VALUE={(V(A,K)-va0)/KP*LN(1+EXP(V_6))}

Eam   am  0  VALUE= {2*Pow(V(8),ERP)/KK1}
GA    A  K  VALUE={V(am)}

D3    5 k DX     ; FOR GRID CURRENT
R1    g 5 {RGI}  ; FOR GRID CURRENT
Rak   A  K  1G
Rgk   G  K  1G

C1    G  K  {CGK}
C2    G  A  {CGA}
C3    A  K  {CAK}
.MODEL DX D(IS=1N RS=1)
.ENDS


ERROR: File 'final.cir' line 64:-
B:x1:B_E6 x1:6 0 V=8.5640000000e+01*( (1/5.3550000000e+00)+((v(g1,k1)-3.0000000000e+00)/sqrt(v(a1,k1)**2+6.5800000000
e+01**2)) )
=> Node 'x1:6' has only one component connected

ERROR: File 'final.cir' line 65:-
B:x1:B_E8 x1:8 0 V=(v(a1,k1))/8.5640000000e+01*ln(1+exp(v(x1:6)))
=> Node 'x1:8' has only one component connected
 
Probably I have found the error because now it run, please check
Code:
* TUBE AMP

* WinSPICE Auto run
*#destroy all
*#run
*#plot v(ou)
*#fourier 1000 v(ou)
*#linearize 
*#spec 0 20000 50 v(ou)
*#let voutdb=db(v(ou))
*#plot voutdb
*#setplot op1
*#print va#branch
*#print v(a1)*va#branch
*#rusage everything

c1  in  g1       3.3uF
x1  a1  g1 k1    GM70
rk  k1   0       1

rg  g1 b1  20k

x2  p1 a1  ou 0 LL1620_90mA  

va p1 0        900v
vp b1 0        -80v

vin in 0 dc 0 ac 90 sin(0 90 1000)
rl ou 0 8

.op 
.ac dec 10 5 1E6 
.tran 5000NS 20200US 200uS

.SUBCKT LL1620_90mA  1 2 3 4
ra1  1   5     308
la1  5   2     40H
ka1  la1 la2   1
la2  6   4     0.0277H 
ra2  3   6     0.2
.ENDS

.SUBCKT k845rca  A G K
+PARAMS: MU=5.355 ERP=1.5
+ KK1=6323  KP=85.64  KVB=65.8 vg0=3 va0=0
+ CGA=13.5E-12 CGK=6E-12 CAK=6.5E-12 RGI=1000;

E1 7 0 Value = {(V(A,K)/KP)*LN((1+ EXP((KP/MU) + ((KP*(V(G,K)+.5))/(SQRT(KVB+V(A,K)*V(A,K)))))))}
G1 A K VALUE = {(((V(7))^ERP)/KK1)*(1+SGN(V(7)))}

RE1 7 0 1e12

D3    5 k DX     ; FOR GRID CURRENT
R1    g 5 {RGI}  ; FOR GRID CURRENT
Rak   A  K  1G
Rgk   G  K  1G

C1    G  K  {CGK}
C2    G  A  {CGA}
C3    A  K  {CAK}
.MODEL DX D(IS=1N RS=1)
.ENDS

.SUBCKT GM70  A G K
+PARAMS: MU=8.037 ERP=1.5
+ KK1=4121 KP=182.25 KVB=34 vg0=-5.7
+ CGA=12.p CGK=8.p CAK=4.p RGI=1000

E1 7 0 Value = {(V(A,K)/KP)*LN((1+ EXP((KP/MU) + ((KP*(V(G,K)+.5))/(SQRT(KVB+V(A,K)*V(A,K)))))))}
G1 A K VALUE = {(((V(7))^ERP)/KK1)*(1+SGN(V(7)))}

RE1 7 0 1e12

D3    5 k DX     ; FOR GRID CURRENT
R1    g 5 {RGI}  ; FOR GRID CURRENT
Rak   A  K  1G
Rgk   G  K  1G

C1    G  K  {CGK}
C2    G  A  {CGA}
C3    A  K  {CAK}
.MODEL DX D(IS=1N RS=1)
.ENDS
 

Attachments

  • gm70.gif
    gm70.gif
    12.6 KB · Views: 277
I tried to include the NormanKoren.lib file to my LT Spice on Apple computer. But it crashes.
  1. the picture foir the file is as if it is an executable (so a windows exc). Not very nice to see;
  2. Look in the file with notepad ist is a text file
  3. when I made the spice directive " .inc Koren.lib LTSpice " crashes.
  4. changing the extension to .inc does not solve it.

How can I include it?
Is this an old file incompatib;le with LTSpice

TUBE LIBRARY

* This library was developed by Norman Koren.
*
* For details, refer to the article, "Improved Vacuum-Tube Models
* for SPICE simulations," Glass Audio, Vol. 8, No. 5, 1996,
* available from Audio Amateur Corporation, 305 Union St.,
* PO Box 176, Peterborough, NH 03458 USA. Phone 603-924-9464.
*
* All the usual legal disclaimers apply. The author has made
* every effort to provide correct information, but assumes no
* liabilities for errors, misuse of the models,
* or inevitable changes made by users.
[ ]
*
* Some models are commented out because the evaluation version of
* Pspice has a maximum of twenty.

.SUBCKT 6550 1 2 3 4 ; P G1 C G2
etc etc​
 
Status
This old topic is closed. If you want to reopen this topic, contact a moderator using the "Report Post" button.