.SUBCKT NTCLE101E3C90173__THERMISTOR__1 1 2 Tin
R1 Tin 0 1T
Vsense 1 N002 0V
B3 N002 2 V= I(Vsense)*Rntc(V(Ttot),beta)
B1 N003 0 V= V(Tin)
B2 Ttot N003 V= +abs(V(vcth))
R2 Ttot 0 1T
I2 Beta 0 -1
R4 Vcth N001 1
C1 0 Vcth {Cth/Gth}
B4 N001 0 V= +abs(V(1,2)*I(Vsense))/Gth
.FUNC Rntc(x,y)
+ {if(x<25,+Rref*exp((y-Beta)*(1/(x+273.15)-1/(273.15+Tref)))*exp(B*(pwr(1/(273.15+x),1)-pwr(1/(273.15+Tref),1))+C*(pwr(1/(273.15+x),2)-pwr(1/(273.15+Tref),2))+D*(pwr(1/(273.15+x),3)-pwr(1/(273.15+Tref),3))),
+ +Rref*exp((y-Beta)*(1/(x+273.15)-1/(273.15+Tref)))*exp(B1*(pwr(1/(273.15+x),1)-pwr(1/(273.15+Tref),1))+C1*(pwr(1/(273.15+x),2)-pwr(1/(273.15+Tref),2))+D1*(pwr(1/(273.15+x),3)-pwr(1/(273.15+Tref),3))))}
.PARAM A=-12.89228328
.PARAM B=4245.148
.PARAM C=-87493
.PARAM D=-9588114
.PARAM a1=-12.89228328
.PARAM b1=4245.148
.PARAM c1=-87493
.PARAM d1=-9588114
.PARAM Gth=0.0068
.PARAM Rref=10000
.PARAM tolB=0.01
.PARAM tolRref=0.01
.PARAM Tref=25
.PARAM Cth=0.0476
.PARAM beta=3435
.IC V(vcth)=0
.ENDS NTCLE101E3C90173
*====================================================
Comments