alainstas
Published © GPL3+

LTspice for PID control of a Thermoelectric Cooling Element

A LTspice XVII simulation of a thermoelectric Peltier cooling element with 100% analogue optimized PID circuit and a Vishay NTC thermistor

IntermediateProtip2 hours3,029
LTspice for PID control of a Thermoelectric Cooling Element

Things used in this project

Hardware components

Vishay NTCALUG03AE3103H
×1

Story

Read more

Custom parts and enclosures

Thermoelectric Cooling Element controlled by PID and NTC Thermistor

LT spice XVII simulation - please download the free LTspice XVII on the Analog Devices Website

Schematics

LTspice simulation of Thermoelectric Element PID control with NTC

download the freeware LTspice XVII on the Analog Devices website

Code

LTspice model of VISHAY NTCALUG03A

Plain text
import it in LTspcie - also to be downloaded on Vishay website ltspice_volt_as_temp_ntcalug03-39.zip
* Voltage_driven_thermistor NTCALUG03A/39A
* Copyright:
* Vishay Intertechnology, Inc.
* customer : general
* Revision:   1.0  $
* Author:   Alain Stas  $
* Date:   14 April 2017 $
* Revision History :correction <> introduced
*---------------------------------------------------------------
* data sheet :http://www.vishay.com/docs/29114/ntcalug3.pdf
* the different R25 and B are as follows : NTCALUGxxA103   = 10Kohms B25/85=3984K
*                                                    103_L = 10Kohms       =3435K
*                                                    123   = 12Kohms       =3740K
*                                                    473   = 47Kohms       =4090K
* The models of this file must be used with the transfer function model transferLUG03-39 in order to reproduce the thermal response when screwed on a surface
* The order of connection : voltage source(surface temperature) connected to pin IN of transferLG03-39 whose pin OUT is connected to Ti
*                electrical pins (to be connected to electrical circuit)
*$                |  |
*                 |  | temperature pin ( for simulation only ,to be connected to pin OUT of transferLUG03-39)
*                 |  | |
.SUBCKT NTC_BASE RN Rp Ti Params: W=1 X=1 Y=1 Z=1 GTH=1 GTH1=1 CTH=1 A=1 R25=1 B=1 C=1 D=1 T0=273.15 TRmc=5 TBmc=1
R_RTR    0  TR  {mc(1,TR-1)}
I_ITR    0 TR 1
R_RTB    0 TB  {mc(1,TB-1)}
I_ITB    0 TB  1
G_G1         AOUT 0 VALUE { if(V(Ti)>25,V(AOUT,
+  0)/(R25*V(TR)*exp(((D*V(TB)/(T0+abs(V(H))+V(Ti))+C*V(TB))/(T0+abs(V(H))+V(Ti))+B*V(TB))/(T0+V(Ti)+abs(V(H)))+A*V(TB))),0)
+  }
G_G2         AOUT 0 VALUE { IF(V(Ti)>25,0,V(AOUT,
+  0)/(R25*V(TR)*exp(((Z*V(TB)/(T0+abs(V(H))+V(Ti))+Y*V(TB))/(T0+abs(V(H))+V(Ti))+X*V(TB))/(T0+abs(V(H))+V(Ti))+W*V(TB))))
+  }
G_G3         H 0 VALUE {
+  if(V(Ti)>25,-V(RP,RN)*V(RP,RN)/(R25*V(TR)*exp(((D*V(TB)/(T0+abs(V(H))+V(Ti))+C*V(TB))/(T0+abs(V(H))+V(Ti))+B*V(TB))/
+  (T0+V(Ti)+abs(V(H)))+A*V(TB))),0)}
G_G4         H 0 VALUE {
+  if(V(Ti)>25,0,-V(RP,RN)*V(RP,RN)/(R25*V(TR)*exp(((Z*TB/(T0+abs(V(H))+V(Ti))+Y*V(TB))/(T0+abs(V(H))+V(Ti))+X*V(TB))/
+  (T0+V(Ti)+abs(V(H)))+W*V(TB))))}
G_G5         RP RN VALUE { V(RP, RN)/V(AOUT) }
G_G6         H 0 VALUE { V(H)*(Gth + Gth1*(V(Ti)-25)) }
I_I1         0 AOUT DC 1Adc
R_R1         0 AOUT  1T TC=0,0
R_R2         0 H  1T TC=0,0
C_C1         0 H  {Cth}
R_R99        0 Ti 1T
.IC V(H) =0
.ENDS
*$
.SUBCKT NTCALUG03A123 RN Rp Ti  PARAMS: TOLR=0 TOLB=0
X43 RN Rp Ti NTC_BASE Params:
+ w=-13.8973229
+ x=4557.725
+ y=-98275
+ z=-7522357
+ gth=0.010 gth1 = 0.00000
+ cth=0.018
+ a=-13.89732294
+ r25=12000
+ b=4557.725
+ c=-98275
+ d=-7522357
+ T0=273.15
+ TR={1+TOLR/100}
+ TB={1+TOLB/100}
.ENDS
*$
.SUBCKT NTCALUG03A103_L RN Rp Ti  PARAMS: TOLR=0 TOLB=0
X45 RN Rp Ti NTC_BASE Params:
+ w=-12.8922832
+ x=4245.148
+ y=-87493
+ z=-9588114
+ gth=0.010 gth1 = 0.00000
+ cth=0.018
+ a=-12.89228328
+ r25=10000
+ b=4245.148
+ c=-87493
+ d=-9588114
+ T0=273.15
+ TR={1+TOLR/100}
+ TB={1+TOLB/100}
.ENDS
*$
.SUBCKT NTCALUG03A103 RN Rp Ti  PARAMS: TOLR=0 TOLB=0
X46 RN Rp Ti NTC_BASE Params:
+ w=-14.6571976
+ x=4798.842
+ y=-115334
+ z=-3730535
+ gth=0.010 gth1 = 0.00000
+ cth=0.018
+ a=-14.65719769
+ r25=10000
+ b=4798.842
+ c=-115334
+ d=-3730535
+ T0=273.15
+ TR={1+TOLR/100}
+ TB={1+TOLB/100}
.ENDS
*$
.SUBCKT NTCALUG03A473 RN Rp Ti  PARAMS: TOLR=0 TOLB=0
X87 RN Rp Ti NTC_BASE Params:
+ w=-13.8973229
+ x=4557.725
+ y=-98275
+ z=-7522357
+ gth=0.010 gth1 = 0.0000
+ cth=0.018
+ a=-13.89732294
+ r25=47000
+ b=4557.725
+ c=-98275
+ d=-7522357
+ T0=273.15
+ TR={1+TOLR/100}
+ TB={1+TOLB/100}
.ENDS
*$
.SUBCKT NTCALUG39A123 RN Rp Ti  PARAMS: TOLR=0 TOLB=0
X43 RN Rp Ti NTC_BASE Params:
+ w=-13.8973229
+ x=4557.725
+ y=-98275
+ z=-7522357
+ gth=0.010 gth1 = 0.00000
+ cth=0.018
+ a=-13.89732294
+ r25=12000
+ b=4557.725
+ c=-98275
+ d=-7522357
+ T0=273.15
+ TR={1+TOLR/100}
+ TB={1+TOLB/100}
.ENDS
*$
.SUBCKT NTCALUG39A103_L RN Rp Ti  PARAMS: TOLR=0 TOLB=0
X45 RN Rp Ti NTC_BASE Params:
+ w=-12.8922832
+ x=4245.148
+ y=-87493
+ z=-9588114
+ gth=0.010 gth1 = 0.00000
+ cth=0.018
+ a=-12.89228328
+ r25=10000
+ b=4245.148
+ c=-87493
+ d=-9588114
+ T0=273.15
+ TR={1+TOLR/100}
+ TB={1+TOLB/100}
.ENDS
*$
.SUBCKT NTCALUG39A103 RN Rp Ti  PARAMS: TOLR=0 TOLB=0
X46 RN Rp Ti NTC_BASE Params:
+ w=-14.6571976
+ x=4798.842
+ y=-115334
+ z=-3730535
+ gth=0.010 gth1 = 0.00000
+ cth=0.018
+ a=-14.65719769
+ r25=10000
+ b=4798.842
+ c=-115334
+ d=-3730535
+ T0=273.15
+ TR={1+TOLR/100}
+ TB={1+TOLB/100}
.ENDS
*$
.SUBCKT NTCALUG39A473 RN Rp Ti  PARAMS: TOLR=0 TOLB=0
X87 RN Rp Ti NTC_BASE Params:
+ w=-13.8973229
+ x=4557.725
+ y=-98275
+ z=-7522357
+ gth=0.010 gth1 = 0.0000
+ cth=0.018
+ a=-13.89732294
+ r25=47000
+ b=4557.725
+ c=-98275
+ d=-7522357
+ T0=273.15
+ TR={1+TOLR/100}
+ TB={1+TOLB/100}
.ENDS
*$

Credits

alainstas
67 projects • 35 followers
product marketing engineer at Vishay. began to simulate in spice programs in 2014.

Comments