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