Re: Peripetie s termočlánky (Vianocna list..

Daniel Valuch balu na k-net.fr
Úterý Prosinec 24 16:13:57 CET 2024


poslane sukromne...

Ano, ja som to tiez minuly rok studoval spracovanie signalov od RTD a 
termoclankov cele vianoce. Zaujimava problematika. Dokumenty priamo z 
BIPM na temu ITS-90 a tabulky koeficientov pre RTD aj termoclanky od 
NISTu su velmi dobry zaklad a su velmi informativne.

Prikladam vycuc z kodu

   if (MeasureTC1) {
       // measure cold junction temperature first
       RColdJunction = (float)ADCavrgCJ / (float)ADCavrgRef * R0; // 
scaling of cold junction and reference cancels out
       TemperatureCJ = 
RTDpolynomial(WCorrection(RColdJunction/R0nominal)); // calculate cold 
junction temperature

       VCJ = ThermoCoupleK_TtoE(TemperatureCJ - 273.16); // equivalent 
Cold Junction voltage
       VTC1raw = (float)ADCavrgTC1 * QueueBitShift * LSBmV * 
TC1InvGain;  // convert measured ADC number to voltage in mV
       VTC1raw = VTC1raw - 
TC1offsetVoltage;                              // suppress analogue 
offset (e.g. from the op-amp)
       VTC1comp = VTC1raw + 
VCJ;                                          // cold junction voltage 
added to measured voltage
       TemperatureTC1 = ThermoCoupleK_EtoT(VTC1comp);
       // 0=Celsius, 1=Kelvin, 2=Farenheit, 3=Ohm, 4=microVolt
       if (unit == 0) {
         PrepareText(TemperatureTC1, 0);
       } else if (unit == 1) {
         PrepareText(TemperatureTC1 + 273.15, 1);
       } else if (unit == 2) {
         PrepareText(TemperatureTC1 * 1.8 + 32.0, 2);
       } else if (unit == 3) {
         PrepareText(VTC1raw, 4);
       }
     }


float ThermoCoupleK_EtoT(float e) {
   // type K ITS-90 polynomial
   // https://srdata.nist.gov/its90/type_k/kcoefficients_inverse.html
   const float d0m = 0.0e0;
   const float d1m = 2.5173462e1;
   const float d2m = -1.1662878e0;
   const float d3m = -1.0833638e0;
   const float d4m = -8.9773540e-1;
   const float d5m = -3.7342377e-1;
   const float d6m = -8.6632643e-2;
   const float d7m = -1.0450598e-2;
   const float d8m = -5.1920577e-4;

   const float d0p = 0.000000e0;
   const float d1p = 2.508355e1;
   const float d2p = 7.860106e-2;
   const float d3p = -2.503131e-1;
   const float d4p = 8.315270e-2;
   const float d5p = -1.228034e-2;
   const float d6p = 9.804036e-4;
   const float d7p = -4.413030e-5;
   const float d8p = 1.057734e-6;
   const float d9p = -1.052755e-8;

   float E = 0.0;
   float t90 = 0.0;

   E = e;

   if (e < 0.00) {  // negative temperature
     t90 = t90 + d1m * E;
     E = E * e;
     t90 = t90 + d2m * E;
     E = E * e;
     t90 = t90 + d3m * E;
     E = E * e;
     t90 = t90 + d4m * E;
     E = E * e;
     t90 = t90 + d5m * E;
     E = E * e;
     t90 = t90 + d6m * E;
     E = E * e;
     t90 = t90 + d7m * E;
     E = E * e;
     t90 = t90 + d8m * E;
   } else if (e >= 0 && e < 20.644) {  // positive temperature <500 C
     t90 = t90 + d1p * E;
     E = E * e;
     t90 = t90 + d2p * E;
     E = E * e;
     t90 = t90 + d3p * E;
     E = E * e;
     t90 = t90 + d4p * E;
     E = E * e;
     t90 = t90 + d5p * E;
     E = E * e;
     t90 = t90 + d6p * E;
     E = E * e;
     t90 = t90 + d7p * E;
     E = E * e;
     t90 = t90 + d8p * E;
     E = E * e;
     t90 = t90 + d9p * E;
   } else {
     t90 = -273.0;
   }

   return t90;
}


On 24/12/2024 15:26, Pavel Hudeček wrote:
> No já si vždycky napřed zkusím nastudovat problematiku sám a až pak se 
> možná ptát.
>
> Předpokládám, že ta implementace je:
>>
>>   * Software do mikrokontroléra - kontaktujte nás
>>
>
> PH
>
>
> Dne 24.12.2024 v 12:47 Daniel Valuch napsal(a):
>> robili sme to minuly rok na ZENITe 
>> https://dvaluch.web.cern.ch/zenit/2024.html
>>
>> aj s implementaciou ITS polynomov do Arduina. Nie je to zlozite, 
>> stacilo sa opytat :-)
>>
>> b.
>>
>>
>> On 24/12/2024 12:08, Pavel Hudeček wrote:
>>> Jo to je pravda, to byla blbost:-)
>>>
>>> Teď, když jsem se začal zabývat tím, jak TČ použít ve svém SW, jsem 
>>> asi zjistil, jakej je skutečnej důvod, proč se ty měřáky u těch -50, 
>>> nebo -80 tolik lišej a stejně tak to očekávám na +300 a víc:
>>>
>>> Tady je komplet všechno hezky vysvětlený, chybí ale tabulka 
>>> koeficientů pro jeden směr:
>>> https://www.ti.com/lit/an/sbaa274a/sbaa274a.pdf?ts=1734992321711&ref_url=https%253A%252F%252Fwww.google.com%252F 
>>>
>>>
>>> Tady jsou obě* a taky pro všechny možný další TČ:
>>> https://www.omega.co.uk/temperature/z/pdf/z198-201.pdf
>>>
>>> *Tím se dostáváme k podstatě:
>>> Nějak jsem si myslel, že vezmu polynom, k němu nějaký koeficienty, 
>>> jedny na U->T, pak druhý na T->U pro kompenzaci druhýho konce a 
>>> hotovo. Jenže těch sad koeficientů je víc pro různý rozsahy a ještě 
>>> se ty rozsahy lišej pro jeden a druhej směr. Celkem třeba 5 ks pro K 
>>> v tom druhým PDF (u jedný chybí - v popisku). No a v případě K je 
>>> tam v jednom směru ještě další, tentokrát exponenciální člen.
>>>
>>> Tohle asi 90 % měřáků nemá komplet implementováno a tak se jejich 
>>> zjednodušený verze lišej:-)
>>> A i kdyby, tak kdoví, jestli koeficienty nevzali odněkud, kde 
>>> nechybělo jen - v popisku, ale třeba byl překlep v tabulce.
>>>
>>> Tady si teda někdo dal tu práci:
>>> https://github.com/IowaScaledEngineering/arduino-thermocouple/blob/master/src/Thermocouple.cpp 
>>>
>>> Ale víme zda zas někde není překlep?
>>>
>>> Tadyten si to přechroustal na tabulky intů:
>>> https://github.com/Miceuz/k-thermocouple-lib/blob/master/thermocouple.c
>>> To bude určitě do MCU rychlejší a menší.
>>> Ale šance zkontrolovat jeho tabulky se blíží 0.
>>>
>>> PH
>>>
>>> Dne 22.12.2024 v 8:23 Jaroslav Buchta napsal(a):
>>>> Noooo, na tom snad nezalezi, z termoclanku tusim lezou pidivolti 
>>>> umerne rozdilu teplot mezi konci, jedno ktery je teply a ktery 
>>>> studeny, predpokladam...
>>>> Dne 21.12.2024 v 20:42 Pavel Hudeček napsal(a):
>>>>> U suchýho ledu se mi taky teploměry docela rozcházely. Oni pak z 
>>>>> toho termočlánku lezou dost pidivolti tak se nedivim.
>>>>>
>>>>> S tím že je to míň než má var nejnižší složka, nemám problém.
>>>>> AZEotropní směs, no:-)
>>>>>
>>>>> PH
>>>>>
>>>>> Dne 21.12.2024 v 20:15 ajtservis napsal(a):
>>>>>> https://www.rajce.idnes.cz/tomik-k/album/zamrazac
>>>>>> v ramci sobotni nepresnosti, kdy jsem musel odlozit drink a 
>>>>>> polonahej jit pro ten cajk do nevytapeny dilny :-) proste zvedavost.
>>>>>> pokus uz v teple kanclu. uz tam toho v tom spreji moc neni.
>>>>>> pdf. https://termopasty.com/en/products/freeze/
>>>>>> chemii nerozumim, koupil jsem si radeji nehorlavej = 4x drazsi. 
>>>>>> jezdil jsem leta na lpg a ten smrad varuje :-)
>>>>>> t.
>>>>>>
>>>>>> Dne 21.12.2024 v 19:29 AZe napsal(a):
>>>>>>> Bod varu dimethyeteru při 1 baru je -24°C. V tom spreji je ještě 
>>>>>>> něco málo (pod 10 proc.) propanu, bod varu -42°C.
>>>>>>> Takže udávaná teplota -55°C v titulku je blbost a v textu 
>>>>>>> dokonce  -67°C je kravínium.
>>>>>>> Navíc je to hořlavé, postříkat s tím cokoli žhavého může být 
>>>>>>> zajímavý zážitek.
>>>>>>> aze
>>>>>>>
>>>>>>> -----Původní zpráva----- From: Pavel Hudeček
>>>>>>> Sent: Saturday, December 21, 2024 5:25 PM
>>>>>>> To: hw-list na list.hw.cz
>>>>>>> Subject: Re: Vianocna listovica bratislava piatok 20.12.
>>>>>>>
>>>>>>> Mě tam zaujal sprej Zamrazacz. Změřili jste kolik umí?
>>>>>>>
>>>>>>> PH
>>>>>>>
>>>>>>> Dne 21.12.2024 v 16:58 ajtservis napsal(a):
>>>>>>>> fotky jsem vystavil sem:
>>>>>>>> https://www.rajce.idnes.cz/tomik-k/album/2x-max6350cpa-mereni-3x- 
>>>>>>>> dmm/1662580062
>>>>>>>> https://www.rajce.idnes.cz/tomik-k/album/2x-max6350cpa-mereni-3x- 
>>>>>>>> dmm/1662580061
>>>>>>>> https://www.rajce.idnes.cz/tomik-k/album/2x-max6350cpa-mereni-3x- 
>>>>>>>> dmm/1662580066
>>>>>>>> je jeste jedna s lidma, tu nevim jestli se dneska smi jen tak 
>>>>>>>> zverejnit ? :-)
>>>>>>>> t.
>>>>>>>>
>>>>>>>>
>>>>>>>> Dne 21.12.2024 v 15:05 Daniel Valuch napsal(a):
>>>>>>>>> poslal som sukromne, velky subor
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> On 21/12/2024 15:02, Jindrich Fucik wrote:
>>>>>>>>>> ... a kam došla fotka slavnostního předávání #1?
>>>>>>>>>>
>>>>>>>>>> Dne 20.12.2024 v 19:36 balu napsal(a):
>>>>>>>>>>> Okalibrovane v krcme. Kam posleme fakturu???
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>> On 13 Dec 2024, at 23:08, ajtservis <info na ajtservis.cz> wrote:
>>>>>>>>>>>>
>>>>>>>>>>>> diky !
>>>>>>>>>>>> vysledky mereni jsou paradni ! a teda zatim i ten muj 
>>>>>>>>>>>> keithley(vsak ma novej vfd).
>>>>>>>>>>>> samozrejme si zelam "aj s protokolom na podpivniku" :-)
>>>>>>>>>>>> patronskemu labaku zdar !
>>>>>>>>>>>> t.
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> Dne 13.12.2024 v 22:45 Daniel Valuch napsal(a):
>>>>>>>>>>>>> Dnes prisiel balik s putovnym etalonom do podatelne na 
>>>>>>>>>>>>> fakulte. Hned som ho vyzdvihol a nechal odmerat mnohomerom 
>>>>>>>>>>>>> Fluke 8588A, kalibracia 19.1.2024.
>>>>>>>>>>>>> Vysledky v texte a prilohe.
>>>>>>>>>>>>> Ak to chcete okalibravovane na zasadani Patronskeho 
>>>>>>>>>>>>> elektrotechnickeho spolku, metrologicka sekcia tak rada 
>>>>>>>>>>>>> urobi aj s protokolom na podpivniku. Ku meraniu sa ale 
>>>>>>>>>>>>> prida este neistota typu C a bude to velmi drahe. Ak 
>>>>>>>>>>>>> klient suhlasi, buduci piatok na listovici urobime. V 
>>>>>>>>>>>>> opacnom pripade zakazku uzatvarame a etalon posielame 
>>>>>>>>>>>>> naspat ku zakaznikovi.
>>>>>>>>>>>>> S pozdravom,
>>>>>>>>>>>>> Kolektiv DMI
>>>>>>>>>>>>>> On 10/12/2024 17:36, ajtservis wrote:
>>>>>>>>>>>>>> hezky den.
>>>>>>>>>>>>>> ja bych poslal do "labaku" 2x ref. 5V MAX6350EPA+.
>>>>>>>>>>>>>> https://www.rajce.idnes.cz/tomik-k/album/2x-max6350cpa- 
>>>>>>>>>>>>>> mereni-3x- dmm/1661677708
>>>>>>>>>>>>>> https://www.rajce.idnes.cz/tomik-k/album/2x-max6350cpa- 
>>>>>>>>>>>>>> mereni-3x- dmm/1661677712
>>>>>>>>>>>>>> https://www.rajce.idnes.cz/tomik-k/album/2x-max6350cpa- 
>>>>>>>>>>>>>> mereni-3x- dmm/1661677711
>>>>>>>>>>>>>> standart 4mm zdirky. napajeni baterka 9V. kdyby neprezila 
>>>>>>>>>>>>>> cestu(neni lithiovka, omylem), na vymenu nutny sroubovak 
>>>>>>>>>>>>>> ph1.
>>>>>>>>>>>>>> v napajeni ochrana pri pripojovani baterky shottky dioda, 
>>>>>>>>>>>>>> lze dip_switchem vyradit/zkratovat. zlousel jsem a hejba 
>>>>>>>>>>>>>> se posledni cislo na dmm o 0,00002V, nevim proc :-)
>>>>>>>>>>>>>> ted jen kam a jakou zvolit spolehlivou sluzbu cz/sk ?
>>>>>>>>>>>>>> t.
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Dne 10.12.2024 v 16:39 Daniel Valuch napsal(a):
>>>>>>>>>>>>>>> co presne budeme kalibrovat? Nech sa mozeme riadne 
>>>>>>>>>>>>>>> pripravit.
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> On 10/12/2024 12:47, Jindrich Fucik wrote:
>>>>>>>>>>>>>>>> A ještě jedna praktická - přivezu ten čínský napěťový 
>>>>>>>>>>>>>>>> "standard". Doufám, že budou k dispozici dostatečné 
>>>>>>>>>>>>>>>> měřící kapacity. Je v oblasti dostatek podpivníkov, aby 
>>>>>>>>>>>>>>>> bylo na co vypracovat protokol? Nebo si mám nějaké dovézt?
>
> _______________________________________________
> HW-list mailing list  -  sponsored bywww.HW.cz
> Hw-list na list.hw.cz
> http://list.hw.cz/mailman/listinfo/hw-list
------------- další část ---------------
HTML příloha byla odstraněna...
URL: <http://list.hw.cz/pipermail/hw-list/attachments/20241224/85672bd7/attachment.htm>


Další informace o konferenci Hw-list