4. NOWA METODA KOMPENSACJI BŁĘDÓW TEMPERATUROWYCH

s. 38

4. 1. Wprowadzenie

Zaproponowana niżej nowa metoda kompensacji błędów temperaturowych [81, 82] eliminuje niedogodności znanych metod, umożliwiając jednocześnie tożsamościowe sprowadzenie błędu temperaturowego przetwornika do zera.

Dla ujednoznacznienia przedstawionych niżej rozważań zostaną podane ograniczenia dla klasy przetworników, których kompensacja będzie analizowana: a/ przedmiotem rozważań będą przetworniki o liniowej charakterystyce przetwarzania ciśnienia na sygnał elektryczny, b/ zakłada się, że jedynym źródłem błędów temperaturowych przetwornika jest czujnik. Jeśli będzie wymagane uwzględnienie dryftów temperaturowych innych członów przetwornika znajdujących się w tej samej temperaturze co czujnik /np. elektronicznego układu przetwarzającego/ powinny one być sprowadzone do poziomu czujnika, a kompensacji powinien podlegać zastępczy błąd czujnika /w rozumieniu pkt. 2 niniejszej pracy/, c/ przyjmuje się addytywność wartości mierzonego parametru i wartości błędu temperaturowego, tzn. zakłada się, że wartość rozpatrywanego parametru w dowolnej temperaturze jest opisana sumą wartości parametru w temperaturze odniesienia oraz błędu temperaturowego, d/ przedmiotem rozważań będą błędy w stanach ustalonych temperatury; nie będzie natomiast analizowany wpływ temperatury w stanie nieustalonym, ponieważ błąd ten jest związany nierozerwalnie z konstrukcją przetwornika.

s. 39

4. 2. Warunki wystarczające kompensacji temperaturowej charakterystyki przetwarzania przetwornika

Niech charakterystyka przetwarzania (2) przetwornika będzie opisana zależnością liniową /rys. 14/:

U=Up+α1X(3)U = U_{p} + \alpha_{1} X \qquad (3)

gdzie: UpU_{p} - wartość początkowa sygnału wyjściowego przetwornika przy wartości początkowej XpX_{p} wielkości nieelektrycznej α1\alpha_{1} - czułość przetwarzania

Zmiany sygnału wyjściowego pod wpływem temperatury określa się równaniem różniczkowym cząstkowym

∂U∂T=∂Up∂T+X∂α1∂T(4)\frac{\partial U}{\partial T} = \frac{\partial U_{p}}{\partial T} + X \frac{\partial \alpha_{1}}{\partial T} \qquad (4)

W celu spełnienia warunku pełnej kompensacji temperaturowej sygnału wyjściowego

∂U∂T=0(5)\frac{\partial U}{\partial T} = 0 \qquad (5)

otrzymuje się dwa rozwiązania równania (4)

a.

∂Up∂T=−X∂α1∂T(6)\frac{\partial U_{p}}{\partial T} = - X \frac{\partial \alpha_{1}}{\partial T} \qquad (6)

b.

{∂Up∂T=0∂α1∂T=0(7)\left\{ \begin{aligned} \frac{\partial U_{p}}{\partial T} &= 0 \\ \frac{\partial \alpha_{1}}{\partial T} &= 0 \end{aligned} \right. \qquad (7)
s. 40
Rys. 14. Charakterystyka przetwarzania przetwornika
Rys. 14. Charakterystyka przetwarzania przetwornika

Warunek a. nie ma znaczenia praktycznego ze względu na zmiany w funkcji wartości mierzonej ciśnienia.

Kompensacja temperaturowa większości przetworników pomiarowych jest oparta na warunku b., który możemy zdefiniować następująco:

Warunkiem wystarczającym do skompensowania temperaturowego charakterystyki przetwarzania (3) jest uniezależnienie od temperatury wartości początkowej UpU_{p} oraz czułości przetwarzania α1\alpha_{1}. Warunek ten możemy zmodyfikować: do kompensacji temperaturowej charakterystyki przetwarzania (3) wystarczające jest uniezależnienie od temperatury wartości początkowej UpU_{p} i wartości końcowej UkU_{k} zakresu pomiarowego przetwornika. Jest to oczywiste, gdyż wartość końcowa UkU_{k} jest funkcją wartości początkowej i czułości:

Uk=Up+α1(Xk−Xp)(8)U_{k} = U_{p} + \alpha_{1} (X_{k} - X_{p}) \qquad (8)

W rzeczywistości charakterystyka przetwarzania wykazuje również pewną nieliniowość 1/

U=Up+α1(X−Xp)+α2(X−Xp)2+α3(X−Xp)3+…αn(X−Xp)n(9)U = U_{p} + \alpha_{1} (X - X_{p}) + \alpha_{2} (X - X_{p})^{2} + \alpha_{3} (X - X_{p})^{3} + \ldots \alpha_{n} (X - X_{p})^{n} \qquad (9)

1. Twierdzenie Weierstrassa/udowodnione w 1885r [83] /: dla każdej funkcji P(X) ciągłej w przedziale XpX_{p}, XkX_{k} i dowolnego ε>0\varepsilon > 0 istnieje wielomian

U(X)=αo+α1X+…+αnXnU(X) = \alpha_{o} + \alpha_{1} X + \ldots + \alpha_{n} X^{n}
s. 41

w związku z czym kompensacja dwuparametrowa nie zapewnia idealnej kompensacji wpływu temperatury na sygnał wyjściowy przetwornika. Ponieważ jednak zarówno wartość dopuszczalnej nieliniowości, jak i błędu temperaturowego są ze sobą związane 2/, wartość członów nieliniowości jest mała

∣Up+α1⋅(X−Xp)∣≫∣α2⋅(X−Xp)2+…+αn⋅(X−Xp)n∣(10)\left| U_{p} + \alpha_{1} \cdot (X - X_{p}) \right| \gg \left| \alpha_{2} \cdot (X - X_{p})^{2} + \ldots + \alpha_{n} \cdot (X - X_{p})^{n} \right| \qquad (10)

i wpływ zmian współczynników α2,α3,…,αn\alpha_{2}, \alpha_{3}, \ldots, \alpha_{n} spowodowanych działaniem temperatury jest pomijalnie mały /na zasadzie małej wyższych rzędów/.

Rozważmy układ blokowy przetwornika z punktu widzenia potrzeb kompensacji /rys. 14/. Czujnik C wytwarza sygnał u uzależniony od wartości mierzonej wielkości nieelektrycznej X i temperatury T zależnością (1). Układ kompensacji D wytwarza sygnał kompensujący k uzależniony od temperatury za pomocą funkcji H(T):

k=H(T)(11)k = H(T) \qquad (11)

Sygnał wyjściowy przetwornika U powstaje w układzie normalizującym N, sumującym sygnał czujnika i układu kompensacji z odpowiednimi współczynnikami proporcjonalności /wzmocnienia/ b1b_{1}, b2b_{2} oraz przesunięcia b3b_{3} wartości początkowej sygnału wyjściowego


stopnia n(ε\varepsilon) taki, że:

max⁡∣U(X)−P(X)∣<ε\max \left| U(X) - P(X) \right| < \varepsilon

a więc dowolną charakterystykę czujników rzeczywistych można dowolnie dokładnie opisać zależnością (9).

2/ Wg PN-73/M-42024 maksymalna dopuszczalna wartość nieliniowości charakterystyki przetwarzania jest o jedną pozycję w ciągu dokładności większa od maksymalnej dopuszczalnej wartości błędu temperaturowego, wyrażonego w procentach zmian wartości zakresu pomiarowego na 10K. Przykładowo: dla przetwornika klasy 0,4% dopuszcza się nieliniowość 0,4% oraz dryft temperaturowy zera i zakresu 0,25%/10K.

s. 42
Rys. 15. Układ blokowy przetwornika z punktu widzenia potrzeb kompensacji
Rys. 15. Układ blokowy przetwornika z punktu widzenia potrzeb kompensacji
U=b1⋅f(X,T)+b2⋅H(T)+b3(12)U = b_{1} \cdot f(X, T) + b_{2} \cdot H(T) + b_{3} \qquad (12)

gdzie: b1b_{1}, b2b_{2}, b3b_{3} są niezależne od temperatury /ich dryft temperaturowy jest włączony do zastępczego błędu temperaturowego czujnika/ i wartości mierzonej.

Różniczkując i przekształcając (3) otrzymujemy wyrażenie na wartość początkową UpU_{p} i czułość przetwarzania α1\alpha_{1}:

{Up=U−X∂U∂Xα1=∂U∂X(13)\left\{ \begin{aligned} U_{p} &= U - X \frac{\partial U}{\partial X} \\ \alpha_{1} &= \frac{\partial U}{\partial X} \end{aligned} \right. \qquad (13)

Wstawiając do (13) wyrażenie (12) otrzymuje się równania:

{Up=b1[f(X,T)−X∂f(X,T)∂X]+b2⋅H(T)+b3α1=b1∂f(X,T)∂X(14)\left\{ \begin{aligned} U_{p} &= b_{1} \left[ f(X, T) - X \frac{\partial f(X, T)}{\partial X} \right] + b_{2} \cdot H(T) + b_{3} \\ \alpha_{1} &= b_{1} \frac{\partial f(X, T)}{\partial X} \end{aligned} \right. \qquad (14)

uzależniające wartość początkową UpU_{p} i czułość przetwarzania α1\alpha_{1} od funkcji f(X, T) i H(T). W celu określenia ograniczeń dla tych funkcji i w związku z przyjęciem liniowego modelu charakterystyki

s. 43

przetwarzania określonego zależnością (3), sprawdzimy dla (14) warunki kompensacji temperaturowej (7). Warunki te doprowadzają do dwóch równań różniczkowych cząstkowych stopnia drugiego:

{∂H(T)∂T=b1b2⋅[X⋅∂2f(X,T)∂X ∂T−∂f(X,T)∂T]∂2f(X,T)∂X⋅∂T=0(15)\left\{ \begin{aligned} \frac{\partial H(T)}{\partial T} &= \frac{b_{1}}{b_{2}} \cdot \left[ X \cdot \frac{\partial^{2} f(X, T)}{\partial X \, \partial T} - \frac{\partial f(X, T)}{\partial T} \right] \\ \frac{\partial^{2} f(X, T)}{\partial X \cdot \partial T} &= 0 \end{aligned} \right. \qquad (15)

po rozwiązaniu których mamy poszukiwane warunki dla funkcji H(T)H(T) oraz f(X,T)f(X, T):

H(T)=b1b2⋅[X⋅∂f(X,T)∂X−f(X,T)]+f1(X)+C1(16)H(T) = \frac{b_{1}}{b_{2}} \cdot \left[ X \cdot \frac{\partial f(X, T)}{\partial X} - f(X, T) \right] + f_{1}(X) + C_{1} \qquad (16)
f(X,T)=f2(X)+f3(T)+C2(17)f(X, T) = f_{2}(X) + f_{3}(T) + C_{2} \qquad (17)

gdzie f1f_{1}, f2f_{2}, f3f_{3} - dowolne funkcje jednej zmiennej C1C_{1}, C2C_{2} - dowolne stałe

Tak więc funkcja H(T)H(T) musi mieć przebieg odwrotny do błędu temperaturowego czujnika, określonego pierwszym członem zależności (16) oraz może być powiększona o dowolną funkcję wielkości mierzonej f1(X)f_{1}(X) a także o dowolną stałą C1C_{1}. Funkcja f(X,T)f(X, T) może natomiast być jedynie sumą niezależnych funkcji ciśnienia f2(X)f_{2}(X) i temperatury f3(T)f_{3}(T), powiększonych o dowolną stałą.

Zależność (17) wyraża więc warunek addytywności wartości mierzonego parametru i błędu temperaturowego, sformułowany w p. 4. 1. c.

s. 44

4. 3. Istota pełnej kompensacji temperaturowej

4. 3. 1. Synteza metody

Wykorzystując opisaną uprzednio możliwość kompensacji temperaturowej charakterystyki przetwarzania przez kompensację wartości początkowej i końcowej zakresu pomiarowego, przedstawimy nową me metodę kompensacji każdego z tych parametrów.

Określmy bezwzględny błąd temperaturowy czujnika jako różnicę pomiędzy wartością sygnału wyjściowego czujnika przy rozpatrywanym ciśnieniu XX w danej temperaturze TT i temperaturze odniesienia ToT_{o}, oznaczając go przez δ(T)\delta(T):

δ(T)=f(X,T)−f(X,To)(18)\delta(T) = f(X, T) - f(X, T_{o}) \qquad (18)

Oznaczmy moduł wartości maksymalnej błędu 18 przez [?:δm\delta_{m} — only the subscript "m" is typed, the δ is missing]

δm=∣δ(T)∣max⁡(19)\delta_{m} = \left| \delta(T) \right|_{\max} \qquad (19)

W analogiczny sposób zdefiniujemy bezwzględny błąd temperaturowy przetwornika Δ(T)\Delta(T)

Δ(T)=g(X,T)−g(X,To)(20)\Delta(T) = g(X, T) - g(X, T_{o}) \qquad (20)

oraz moduł wartości maksymalnej błędu (20) przez Δm\Delta_{m}

Δm=∣Δ(T)∣max⁡(21)\Delta_{m} = \left| \Delta(T) \right|_{\max} \qquad (21)

Jeśli oznaczymy przez UpU_{p} i UkU_{k} odpowiednio wartości sygnału wyjściowego czujnika odpowiadającą początkowej i końcowej wartości zakresu mierzonego ciśnienia XX, to efektem działania układu kompensacji będzie zmniejszenie maksymalnego względnego błędu pomiarowego przetwornika, w porównaniu z maksymalnym błędem względem czujnika /rys. 16/:

δm∣Uk−Up∣>Δm∣Uk−Up∣(22)\frac{\delta_{m}}{\left| U_{k} - U_{p} \right|} > \frac{\Delta_{m}}{\left| U_{k} - U_{p} \right|} \qquad (22)
s. 45
Rys. 16. Ilustracja działania układu kompensacji
Rys. 16. Ilustracja działania układu kompensacji

Wstawiając do zależności (20) wyrażenie (12) z uwzględnieniem (18) , otrzymujemy określenie błędu temperaturowego przetwornika w postaci

Δ(T)=b1⋅δ(T)+b2⋅H(T)(23)\Delta(T) = b_{1} \cdot \delta(T) + b_{2} \cdot H(T) \qquad (23)

Niech układ kompensacji D /rys. 15/ zbudowany będzie z /n+1/ generatorów HoH_{o}, H1H_{1}, ....., HnH_{n} /rys. 17/. Każdy z generatorów wytwarza sygnały kompensujące odpowiednio ho(T)h_{o}(T),..., hn(T)h_{n}(T). Sygnał wyjściowy H(T) układu kompensacji stanowi sumę tych sygnałów:

H(T)=∑i=0nhi(T)(24)H(T) = \sum_{i=0}^{n} h_{i}(T) \qquad (24)

Załóżmy, że sygnały kompensujące ho(T)h_{o}(T),..., h1(T)h_{1}(T)[?:likely hn(T)] są wzajemnie ze sobą powiązane za pomocą cząstkowych funkcji

s. 46

fo(T)f_{o}(T),....,fn(T)f_{n}(T), uzależnionych od temperatury, w sposób następujący

ho(T)=woh1(T)=w1⋅fo(T)hi(T)=wi⋅∏j=0i=1fj(T)hn(T)=wn⋅∏j=0n−1fj(T)}(25)\left. \begin{aligned} h_{o}(T) &= w_{o} \\ h_{1}(T) &= w_{1} \cdot f_{o}(T) \\ h_{i}(T) &= w_{i} \cdot \prod_{j=0}^{i=1} f_{j}(T) \\ h_{n}(T) &= w_{n} \cdot \prod_{j=0}^{n-1} f_{j}(T) \end{aligned} \right\} \qquad (25)

[?:upper limit of the product in the third line is typed "i=1", likely i-1]

Rys. 17. Układ blokowy układu kompensacji
Rys. 17. Układ blokowy układu kompensacji

gdzie: wow_{o}, w1w_{1},...wnw_{n} - stałe ii, jj, nn - liczby naturalne

Po uwzględnieniu (25) i (24) , wyrażenie opisujące błąd temperaturowy przetwornika (23) przybiera postać

Δ(T)=b1⋅δ(T)+b2⋅wo+b2⋅∑i=1n[wi⋅∏j=0i−1fj(T)](26)\Delta(T) = b_{1} \cdot \delta(T) + b_{2} \cdot w_{o} + b_{2} \cdot \sum_{i=1}^{n} \left[ w_{i} \cdot \prod_{j=0}^{i-1} f_{j}(T) \right] \qquad (26)

Takie uzależnienie sygnałów kompensujących charakteryzuje się następującą właściwością:

TWIERDZENIE 1

Jeśli w pewnej temperaturze TxT_{x} wartość funkcji cząstkowej fk(Tx)=0f_{k}(T_{x})=0,

s. 47

gdzie kk - liczba naturalna k<nk < n, oraz wartość stałej wk≠0w_{k} \neq 0, to funkcje kompensujące hm(Tx)=0h_{m}(T_{x}) = 0 dla naturalnych mm z przedziału k+1⩽m⩽nk + 1 \leqslant m \leqslant n, gdzie nn jest liczbą funkcji cząstkowych występujących w (25) .

D o w ó d: zapiszmy na podstawie (25)

hm(Tx)hk(Tx)=wmwk⋅fk(Tx)⋅∏j=k+1mfj(Tx)\frac{h_{m}(T_{x})}{h_{k}(T_{x})} = \frac{w_{m}}{w_{k}} \cdot f_{k}(T_{x}) \cdot \prod_{j=k+1}^{m} f_{j}(T_{x})

skąd

hm(Tx)=wmwk⋅fk(Tx)⋅hk(Tx)⋅∏j=k+1mfj(Tx)h_{m}(T_{x}) = \frac{w_{m}}{w_{k}} \cdot f_{k}(T_{x}) \cdot h_{k}(T_{x}) \cdot \prod_{j=k+1}^{m} f_{j}(T_{x})

dla k+1⩽m⩽nk + 1 \leqslant m \leqslant n,

więc gdy fk(Tx)=0f_{k}(T_{x}) = 0, to również hm(Tx)=0h_{m}(T_{x}) = 0

c. b. d. d.

Właściwość ta zostanie wykorzystana do jednokierunkowego uniezależnienia od siebie nastaw układu kompensacji w różnych temperaturach. W tym celu uzależnijmy funkcje cząstkowe fj(T)f_{j}(T) od temperaratury odniesienia ToT_{o} i temperatur kompensacji T1T_{1},..., Tn−1T_{n-1} następująco:

fj(T)=aj(T−Tj)(27)f_{j}(T) = a_{j}(T - T_{j}) \qquad (27)

gdzie aja_{j} - stałe j=0,1,...,n−1j = 0, 1, ..., n-1

Przyjęta forma (27) funkcji cząstkowych fj(T)f_{j}(T) ma właściwość tożsamościowego przyjmowania wartości równej zeru w temperaturach TjT_{j}

fj(Tj)≡0(28)f_{j}(T_{j}) \equiv 0 \qquad (28)
s. 48

Po uwzględnieniu (27), wyrażenie opisujące błąd temperaturowy przetwornika (26) przyjmuje postać

Δ(T)=b1⋅δ(T)+b2wo+b2⋅∑i=1nwi⋅∏j=0i−1[aj(T−Tj)](29)\Delta(T) = b_{1} \cdot \delta(T) + b_{2} w_{o} + b_{2} \cdot \sum_{i=1}^{n} w_{i} \cdot \prod_{j=0}^{i-1} \left[ a_{j}(T - T_{j}) \right] \qquad (29)

4. 3. 2. Procedura kompensacji

Przeprowadźmy obecnie zerowanie błędu temperaturowego przetwornika w temperaturach T0T_{0},..., TnT_{n} za pomocą zmiany wartości współczynników w1w_{1} [?:likely wiw_{i}], przy czym dla uproszczenia rozważań posługiwać się będziemy zależnością (26) opisującą błąd temperaturowy przetwornika. Otrzymujemy kolejno: T=ToT = T_{o}; temperatura odniesienia

Wykorzystując (28) stwierdzimy

fo(To)=0f_{o}(T_{o}) = 0

i po uwzględnieniu (26) otrzymamy wartość błędu temperaturowego

Δ(To)=b1⋅δ(To)+b2⋅wo\Delta(T_{o}) = b_{1} \cdot \delta(T_{o}) + b_{2} \cdot w_{o}

Ponieważ z założenia δ(To)=0\delta(T_{o}) = 0, więc zerowanie się błędu Δ(To)\Delta(T_{o}) nastąpi dla

wo=0(30)w_{o} = 0 \qquad (30)

Współczynnik wow_{o} mogłby więc być pominięty w dalszej analizie. Ponieważ jednak w temperaturze odniesienia powinno się przeprowadzić normalizację rozważanych punktów charakterystyki przetwarzania, np. zerowanie lub ustawianie wymaganej wartości sygnału zakresowego, w dalszych rozważaniach pozostawiono współczynnik wow_{o} jako symbol sprowadzania błędu przetwarzania do zera /a nie błędu temperaturowego/.

s. 49

T=T1T = T_{1}: temperatura pierwszej kompensacji

W tej temperaturze

f1(T1)=0f_{1}(T_{1}) = 0

w wyniku czego błąd temperaturowy (26), po uwzględnieniu warunku dla wow_{o} otrzymanego przy T=ToT=T_{o}, zostaje określony:

Δ(T1)=b1⋅δ(T1)+b2⋅w1⋅fo(T1)\Delta(T_{1}) = b_{1} \cdot \delta(T_{1}) + b_{2} \cdot w_{1} \cdot f_{o}(T_{1})

Zerowanie się błędu Δ(T1)\Delta(T_{1}) uzyskujemy dla wartości w1w_{1}:

w1=−b1  (T1)b2  fo(T1)(31)w_{1} = - \frac{b_{1} \; (T_{1})}{b_{2} \; f_{o}(T_{1})} \qquad (31)

[?:in (31) the numerator is typed "b1 /T1/" with the δ missing, likely δ(T1)\delta(T_{1})]

Wymagana wartość w1w_{1} uzależniona jest jedynie od stosunku wartości błędu czujnika δ(T)\delta(T) do wartości funkcji cząstkowej fo(T)f_{o}(T) w tej temperaturze, zmodyfikowanych przez współczynniki układu normalizującego b1b_{1}, b2b_{2}

T=T2T = T_{2}: temperatura drugiej kompensacji

Analogicznie jak poprzednio otrzymujemy

f2(T2)=0f_{2}(T_{2}) = 0

i po uwzględnieniu (26) otrzymujemy wyrażenie na błąd temperaturowy

Δ(T2)=b1⋅δ(T2)+b2⋅[wo+w1⋅fo(T2)+w2⋅fo(T2)⋅f1(T2)]\Delta(T_{2}) = b_{1} \cdot \delta(T_{2}) + b_{2} \cdot \left[ w_{o} + w_{1} \cdot f_{o}(T_{2}) + w_{2} \cdot f_{o}(T_{2}) \cdot f_{1}(T_{2}) \right]

skąd po przyrównaniu do zera i podstawieniu warunków (30) i (31) dla stałych wow_{o} i w1w_{1} uzyskuje się warunek dla stałej w2w_{2}

w2=−b1b2⋅δ(T2)−δT1fo(T1)⋅fo(T2)fo(T2)⋅f1(T2)(32)w_{2} = - \frac{b_{1}}{b_{2}} \cdot \frac{\delta(T_{2}) - \frac{\delta T_{1}}{f_{o}(T_{1})} \cdot f_{o}(T_{2})}{f_{o}(T_{2}) \cdot f_{1}(T_{2})} \qquad (32)

[?:in (32) the inner numerator is typed "δ T1" without the slashes, likely δ(T1)\delta(T_{1})]

s. 50

Tak więc wartość stałej w2w_{2} wymagana do wyzerowania błędu tempe- [?:typed "tempera", the "ra" overstruck with a hyphen] raturowego przetwornika w temperaturze T2T_{2}, uzależniona jest od wartości błędu czujnika w pierwszej i drugiej temperaturze kompensacji /odpowiednio δ(T1)\delta(T_{1}) i δ(T2)\delta(T_{2}) oraz od wartości funkcji cząstkowych fo(T)f_{o}(T), f1(T)f_{1}(T) w obu temperaturach. Nie jest natomiast uzależniona od wartości błędu czujnika i wartości funkcji cząstkowych w następnych temperaturach kompensacji. Jednocześnie należy podkreślić, że ustalenie wartości stałej w2w_{2} w temperaturze T2T_{2} nie wpływa na wartość stałej w1w_{1}, ustalonej w temperaturze pierwszej kompensacji T1T_{1}.

T=T3T = T_{3}.....TkT_{k}....Tn−1T_{n-1}: temperatury kolejnych kompensacji

Postępując w sposób analogiczny do wyżej przedstawionego, uzyskuje się w każdej z temperatur kompensacji TkT_{k} warunek dla stałej wkw_{k} /k = 3, 4,....., n-1/ w postaci ogólnej

wk=−b1⋅δ(Tk)+b2⋅wo+b2⋅∑i=1k−1wi⋅∏j=0i−1fj(Tk)b2⋅∏j=0k−2fj(Tk)(33)w_{k} = - \frac{b_{1} \cdot \delta(T_{k}) + b_{2} \cdot w_{o} + b_{2} \cdot \sum_{i=1}^{k-1} w_{i} \cdot \prod_{j=0}^{i-1} f_{j}(T_{k})}{b_{2} \cdot \prod_{j=0}^{k-2} f_{j}(T_{k})} \qquad (33)

Podobnie jak poprzednio, wartość stałej wkw_{k} niezbędna dla uzyskania zerowej wartości błędu temperaturowego przetwornika w temperaturze TkT_{k} jest uzależniona od wartości stałych wow_{o},...wk−1w_{k-1} ustalonych w poprzednich temperaturach i od wartości funkcji cząstkowych fo(T)f_{o}(T),......, fk−1(T)f_{k-1}(T) w temperaturze TkT_{k} oraz od wartości współczynników układu normalizującego b1b_{1}, b2b_{2}. Wartość stałej wkw_{k} nie zależy jednak od wartości stałych wk+1w_{k+1},....., wn−1w_{n-1} i funkcji cząstkowych fk(T)f_{k}(T),....., fn−2(T)f_{n-2}(T) w następnych temperaturach kompensacji Tk+1T_{k+1},.....Tn−1T_{n-1}.

Wykorzystując tę właściwość do praktycznego zastosowania przedstawionej metody stwierdzimy, że istnieje jednokierunkowa niezależność nastaw elementów ustalających wartość stałych wow_{o},...., wn−1w_{n-1}.

s. 51

T=TnT = T_{n}: temperatura ostatniej kompensacji

W poprzednich temperaturach zostały wyczerpane własności (28) zerowania się wszystkich istniejących funkcji cząstkowych fj(T)f_{j}(T), /j=0, 1,..., n-1/, przy ustalaniu niezbędnych dla skompensowania błędu wartości stałych wiw_{i}, /i = 1, 2,...,n-1/. Dobierzmy obecnie stałą wnw_{n} tak, aby po raz ostatni uzyskać zerową wartość błędu Δ(Tn)\Delta(T_{n})

Δ(Tn)=b1⋅δ(Tn)+b2⋅wo+b2⋅∑i=1n−1wi⋅∏j=0i−1fj(Tn)+b2⋅wn⋅∏j=0n−2fj(Tn)\Delta(T_{n}) = b_{1} \cdot \delta(T_{n}) + b_{2} \cdot w_{o} + b_{2} \cdot \sum_{i=1}^{n-1} w_{i} \cdot \prod_{j=0}^{i-1} f_{j}(T_{n}) + b_{2} \cdot w_{n} \cdot \prod_{j=0}^{n-2} f_{j}(T_{n})

skąd

wn=−b1⋅δ(Tn)+b2⋅wo+b2⋅∑i=1n−1wi⋅∏j=0i−1fj(Tn)b2⋅∏j=0n−2fj(Tn)(34)w_{n} = - \frac{b_{1} \cdot \delta(T_{n}) + b_{2} \cdot w_{o} + b_{2} \cdot \sum_{i=1}^{n-1} w_{i} \cdot \prod_{j=0}^{i-1} f_{j}(T_{n})}{b_{2} \cdot \prod_{j=0}^{n-2} f_{j}(T_{n})} \qquad (34)

Warunek (34)ma analogiczną postać i właściwości, jak wyrażenie (33) z tą tylko różnicą, że bez następnych temperatur kompensacji. W wyniku zastosowania przedstawionej metody otrzymuje się zerową wartość błędu temperaturowego rozważanego punktu charakterystyki przetwornika pomiarowego w n temperaturach T1T_{1}, T2T_{2},.... ..., TnT_{n}, co oznacza, że w tych temperaturach wartość rozpatrywanego punktu charakterystyki równa jest jego wartości w temperaturze odniesienia ToT_{o}. W interpretacji graficznej działania przedstawionej metody kompensacji wystąpi przecinanie się wykresu odwzorowującego zależność składowej błędu czujnika b1b_{1} (T) [?:δ\delta missing] w sygnale wyjściowym przetwornika z wykresem składowej sygnału kompensującego b2b_{2} H(T)H(T) o zmienionym znaku, w /n+1/ punktach temperatury /rys. 18/.

s. 52
Rys. 18. Interpretacja graficzna działania kompensacji n-punktowej
Rys. 18. Interpretacja graficzna działania kompensacji n-punktowej

4. 3. 3. Warunki wystarczające kompensacji temperaturowej pełnej

Warunek zerowego błędu temperaturowego przetwornika w n temperaturach kompensacji nie jest wystarczający do określenia jego maksymalnej wartości. Niezbędna jest bowiem znajomość wartości błędu w całym rozpatrywanym zakresie temperatur pracy. Obecnie rozważy się warunki wystarczające kompensacji temperaturowej pełnej w sensie minimalizacji wartości błędu temperaturowego. Wprowadźmy następujący opis błędu temperaturowego czujnika δ(T)\delta(T)

δ(T)=co+c1T+c...+cs  Ts(35)\delta(T) = c_{o} + c_{1} T + c ... + c_{s} \; T^{s} \qquad (35)

gdzie coc_{o},..., csc_{s} - stałe

Na mocy twierdzenia Weierstrassa [83] , /patrz pkt. 4. 2. / wyrażenie (35) przybliża dostatecznie dokładnie dowolną krzywą ciągłą, przy czym rząd s wielomianu aproksymującego jest funkcją dokładności odwzorowania. Załóżmy więc, że opis (35) błędu temperaturowego czujnika jest dostatecznie dokładnym odwzorowaniem rzeczywistego błędu (18) czujnika. Porównując wyrażenia (29) i (23) możemy określić sygnał kompensujący następująco

s. 53
H(T)=wo+∑i=1nwi⋅∏j=0i−1aj⋅(T−Tj)(36)H(T) = w_{o} + \sum_{i=1}^{n} w_{i} \cdot \prod_{j=0}^{i-1} a_{j} \cdot (T - T_{j}) \qquad (36)

Wyrażenie (36) jest wielomianem stopnia n. Tak więc błąd temperaturowy (23) przetwornika jest sumą dwóch wielomianów (35),(36), zmodyfikowanych przez wprowadzenie współczynników b1b_{1}, b2b_{2} układu normalizującego. Minimalną wartością, jaką może przyjąć suma dwóch wielomianów jest zero. Jeśli dwa wielomiany są tego samego stopnia i współczynniki przy zmiennej również tego samego stopnia mają jednakową wartość modułu o przeciwnych znakach, to otrzyma się zerową wartość sumy tych wielomianów niezależnie od wartości zmiennej. Możemy więc sformułować wniosek, że interesującym nas warunkiem optymalnym kompensacji temperaturowej jest uzyskanie tożsamościowej równości

b1⋅δ(T)+b2⋅H(T)≡0(37)b_{1} \cdot \delta(T) + b_{2} \cdot H(T) \equiv 0 \qquad (37)

Warunkiem koniecznym dla osiągnięcia (37) jest równość rzędów wielomianów (35) i (36):

n=s(38)n = s \qquad (38)

co jest równoznaczne ze zrównaniem liczby n temperatur kompensacji z rzędem s wielomianu aproksymującego błąd temperaturowy czujnika. W wyniku przeprowadzenia kompensacji zgodnie z pkt. 4. 3. 2. uzyskuje się /n+1/ zer równania:

Δ(T)=0\Delta(T) = 0

Ponieważ przez /n+1/ punktów można poprowadzić tylko jedną krzywą stopnia n, oznacza to, że procedura kompensacji doprowadzi w sposób jednoznaczny do równości tożsamościowej (37). Jest to równoważne z tożsamościowym zerowaniem się błędu temperaturowego przetwornika w całym zakresie temperatur pracy, niezależnie od wyboru temperatur kompensacji. W przeprowadzonej powyżej analizie założono, że generatory sygnałów kompensujących H1H_{1}...HnH_{n} wytwarzają sygnały uzależnione liniowo od temperatury. Wypadkowy sygnał kompensujący jest wytwarzany

s. 54

przez ustalenie odpowiedniej kombinacji iloczynów tych sygnałów, dając w wyniku tego wielomian stopnia nn. Załóżmy obecnie, że chcemy uzyskać wielomian stopnia nn przy mniejszej liczbie temperatur kompensacji. W tym celu uzależnijmy funkcje cząstkowe fj(T)f_{j}(T) od temperatury w sposób następujący:

fj(T)=(T−Tj)⋅∑p=0m−1ajp⋅(T−Tj)p(39)f_{j}(T) = (T - T_{j}) \cdot \sum_{p=0}^{m-1} a_{jp} \cdot (T - T_{j})^{p} \qquad (39)

gdzie ajpa_{jp} - stałe j=0,1,...,n−1j = 0, 1, ..., n-1
mm - liczba naturalna

Wyrażenie (39) jest wielomianem stopnia mm i na mocy twierdzenia Weierstrassa[83] /patrz pkt. 4. 2. /, może ono opisać każdą funkcję ciągłą, przy czym rząd wielomianu mm jest uzależniony od wymaganej dokładności odwzorowania tej funkcji. Przyjęta forma (39) funkcji cząstkowych fj(T)f_{j}(T) ma również właściwość (28) tożsamościowego zerowania się w temperaturach TjT_{j}. Załóżmy z kolei, że znana jest wstępnie pewna ilość RR własności zależności (35) aproksymującej błąd temperaturowy czujnika. Niech rząd mm funkcji cząstkowej (39) określony będzie następująco

m=R+r(40)m = R + r \qquad (40)

gdzie rr jest liczbą naturalną

Podwyższając rząd mm uzyskuje się wielomian kompensujący H(T)H(T) stopnia n(R+r)n(R+r). W celu uzyskania tożsamości (37) musi być spełniona zależność

s=n(R+r)(41)s = n(R + r) \qquad (41)

skąd uzyskuje się wymaganą liczbę nn temperatur kompensacji w przypadku przyjęcia r=0r = 0, a więc w przypadku maksymalnego wykorzystania RR znanych własności:

n=sR(42)n = \frac{s}{R} \qquad (42)
s. 55

55

Wnioski: optymalne warunki kompensacji temperaturowej pełnej można osiągnąć bądź przez wytworzenie liniowych w funkcji temperatury funkcji cząstkowych fj(T)f_j(T) i przeprowadzanie kompensacji w liczbie temperatur równej rzędowi wielomianu opisującego błąd temperaturowy czujnika /warunek 38/, bądź też przez podwyższenie stopnia wielomianu opisującego funkcję cząstkową fj(T)f_j(T) opierając się na znanych wstępnie własnościach błędu temperaturowego czujnika i zmniejszeniu liczby temperatur kompensacji /warunek 42/.

Przykładem wstępnej informacji o charakterze błędu temperaturowego czujnika jest informacja o stosunku wartości błędów w pewnych temperaturach. I tak w zintegrowanych czujnikach ciśnienia [40] dryft temperaturowy zakresu aproksymowany jest krzywą czwartego rzędu i iloraz wartości błędu temperaturowego w temperaturze 352K i 241K jest funkcją rezystywności materiału podłoża. W takim przypadku wystarczające może się okazać wytworzenie funkcji fj(T)f_j(T) kompensującej błąd temperaturowy jednego czujnika, co umożliwi następnie kompensację wszystkich czujników za pomocą mnożenia /dzielnikowania/ wartości sygnału kompensującego przez odpowiedni współczynnik proporcjonalności.

s. 56

56

Metoda zostanie obecnie wyjaśniona na przykładach.

Przykład 1

Niech błąd czujnika będzie określony zależnością kwadratową

δ(T)=B⋅(T−To)⋅(T+Tx)\delta(T) = B \cdot (T - T_{o}) \cdot (T + T_{x})

gdzie TxT_{x} - stała oraz niech b1=b2=1b_{1} = b_{2} = 1

Zrealizujemy dwa warianty układu kompensacji dla m⋅n=2m \cdot n = 2;poddając je następnie kompensacji przeprowadzimy dla nich procedury kompensacji

wariant a/ n=2n = 2 m=1m = 1, więc

fo=a1⋅(T−To)f_{o} = a_{1} \cdot (T - T_{o})
f1=a2⋅(T−T1)f_{1} = a_{2} \cdot (T - T_{1})

a po wstawieniu do( 26) otrzymamy

Δ(T)=B⋅(T−To)⋅(T+Tx)+wo+w1⋅a1⋅(T−To)−w2⋅a1⋅a2⋅(T−To)⋅(T−T1)\Delta(T) = B \cdot (T - T_{o}) \cdot (T + T_{x}) + w_{o} + w_{1} \cdot a_{1} \cdot (T - T_{o}) - w_{2} \cdot a_{1} \cdot a_{2} \cdot (T - T_{o}) \cdot (T - T_{1})

Procedura kompensacji:

T=To;Δ(To)=woT = T_{o}; \quad \Delta(T_{o}) = w_{o}

a stąd uzyskujemy warunek zerowania się błędu Δ(To)\Delta(T_{o})

Wo=0W_{o} = 0

[?:likely wow_{o}; capital W as printed]

T=T1;Δ(T1)=B⋅(T1−To)⋅(T1+Tx)+w1⋅a1⋅(T−To)T = T_{1}; \quad \Delta(T_{1}) = B \cdot (T_{1} - T_{o}) \cdot (T_{1} + T_{x}) + w_{1} \cdot a_{1} \cdot (T - T_{o})

oraz warunek zerowania się błędu Δ(T1)\Delta(T_{1})

w1=−B⋅(T1+Txa1w_{1} = - \frac{B \cdot (T_{1} + T_{x}}{a_{1}}

[?:closing slash after TxT_{x} missing in print; the "+" is overtyped]

T=T2;T = T_{2};
Δ(T2)=B⋅(T2−To)⋅(T2+Tx)+B⋅(T1+Tx)⋅(T2−To)−a1⋅a2⋅w2⋅(T2−To)⋅(T2−T1)\Delta(T_{2}) = B \cdot (T_{2} - T_{o}) \cdot (T_{2} + T_{x}) + B \cdot (T_{1} + T_{x}) \cdot (T_{2} - T_{o}) - a_{1} \cdot a_{2} \cdot w_{2} \cdot (T_{2} - T_{o}) \cdot (T_{2} - T_{1})
s. 57

57

i wreszcie warunek zerowania się błędu Δ(T2)\Delta(T_{2})

w2=−Ba1⋅a2w_{2} = - \frac{B}{a_{1} \cdot a_{2}}

Podstawiając uzyskane w ten sposób wartości wow_{o}, w1w_{1}, w2w_{2} do wyrażenia (26) uzyskuje się tożsamość

Δ(T)≡0\Delta(T) \equiv 0

wariant b/ n=1n = 1, m=2m = 2

fo(T)=(T−To)⋅(ao1+ao2T)f_{o}(T) = (T - T_{o}) \cdot (a_{o1} + a_{o2} T)

Po wstawieniu do (26) uzyskamy

Δ(T)=B⋅(T−To)⋅(T+Tx)+wo+w1⋅(T−To)⋅(ao1+ao2⋅T)\Delta(T) = B \cdot (T - T_{o}) \cdot (T + T_{x}) + w_{o} + w_{1} \cdot (T - T_{o}) \cdot (a_{o1} + a_{o2} \cdot T)

Procedura kompensacji

T=To;Δ(To)=woskądT = T_{o}; \quad \Delta(T_{o}) = w_{o} \qquad \text{skąd}
wo=0w_{o} = 0
T=T1;Δ(T1)=B⋅(T1+To)⋅(T1+Tx)+w1⋅(T1−To)⋅(ao1+ao2⋅T1)T = T_{1}; \quad \Delta(T_{1}) = B \cdot (T_{1} + T_{o}) \cdot (T_{1} + T_{x}) + w_{1} \cdot (T_{1} - T_{o}) \cdot (a_{o1} + a_{o2} \cdot T_{1})

[?:likely T1−ToT_{1} - T_{o} in the first factor; "+" as printed]

skąd

w1=−B⋅(T1+Txao1+a02⋅T1w_{1} = - \frac{B \cdot (T_{1} + T_{x}}{a_{o1} + a_{02} \cdot T_{1}}

[?:closing slash after TxT_{x} missing in print; denominator subscript typed "02" (digit zero) as printed, likely ao2a_{o2}]

Po uwzględnieniu otrzymanych wartości wow_{o} i w1w_{1} otrzymamy

Δ(T)=B⋅(T−To)⋅[T+Tx−(T1+Tx)⋅(T+ao1ao2)T1+ao1ao2]\Delta(T) = B \cdot (T - T_{o}) \cdot \left[ T + T_{x} - \frac{(T_{1} + T_{x}) \cdot (T + \frac{a_{o1}}{a_{o2}})}{T_{1} + \frac{a_{o1}}{a_{o2}}} \right]

Otrzymana wartość Δ(T)\Delta(T) jest tożsamościowo równa zeru dla warunku

s. 58

58

dodatkowego, uzależniają jego charakterystykę błędu czujnika z charakterystyką sygnału kompensującego w postaci

Tx=ao1ao2T_{x} = \frac{a_{o1}}{a_{o2}}

Jest to oczywiste, gdyż pełna kompensacja błędu opisanego funkcją kwadratową wymaga trzech warunków swobody /zrównoważenia trzech współczynników ToT_{o}, TxT_{x}, B, podczas gdy przy procedurze kompensacji w wariancie b mamy jedynie dwa warunki swobody, a mianowicie współczynniki wow_{o} i w1w_{1}. Procedura przedstawiona zgodnie z wariantem a jest oparta na trzech warunkach swobody - współczynnikami wow_{o}, w1w_{1} i w2w_{2} - i w związku z tym nie był wymagany żaden dodatkowy warunek.

4. 4. Istota niepełnej kompensacji temperaturowej

Rozważymy obecnie wartość błędu (26) przetwornika przy założeniu niepełnej kompensacji. Pod nazwą niepełna kompensacja rozumiane będzie ograniczenie liczby temperatur kompensacji w porównaniu z wartościami uzyskiwanymi z zależności (38) lub (42).

m⋅n<s(43)m \cdot n < s \qquad (43)

Dla uproszczenia kompensacja ta zostanie pokazana w przypadku rzędu s wielomianu aproksymującego błąd temperaturowy czujnika równego dwa i trzy, co praktycznie wyczerpuje prawie wszystkie przypadki spotykane w praktyce. W przypadku s=1s = 1 błąd temperaturowy sprowadzany jest tożsamościowo do zera układem kompensacji o parametrach m=1m = 1, n=1n = 1.

4. 4. 1. Przypadek: s=2s = 2, m=1m = 1, n=1n = 1

W tym przypadku błąd czujnika opisany jest zależnością

δ(T)=B⋅(T−To)⋅(T+Tx)(44)\delta(T) = B \cdot (T - T_{o}) \cdot (T + T_{x}) \qquad (44)

gdzie ToT_{o} - temperatura odniesienia B, TxT_{x} - stałe

s. 59

59

Na podstawie (27) mamy fo=ao⋅(T−To)f_{o} = a_{o} \cdot (T - T_{o}) i wykorzystując (26) zapiszemy błąd Δ(T)\Delta(T) przetwornika w postaci

Δ(T)=b1⋅B⋅T2+[b1⋅B⋅(Tx−To)+b2⋅w1⋅ao]⋅T−b1⋅B⋅To⋅Tx−b2⋅w1⋅ao⋅To\Delta(T) = b_{1} \cdot B \cdot T^{2} + \left[ b_{1} \cdot B \cdot (T_{x} - T_{o}) + b_{2} \cdot w_{1} \cdot a_{o} \right] \cdot T - b_{1} \cdot B \cdot T_{o} \cdot T_{x} - b_{2} \cdot w_{1} \cdot a_{o} \cdot T_{o}

Procedura kompensacji jest następująca:

a. T=ToT = T_{o}; Δ(To)=0\Delta(T_{o}) = 0

b. T=T1T = T_{1};

Δ(T1)=b1⋅B⋅[T12+(Tx−To)⋅T1−To⋅Tx]−w1⋅b2⋅ao⋅(To−T1)\Delta(T_{1}) = b_{1} \cdot B \cdot \left[ T_{1}^{2} + (T_{x} - T_{o}) \cdot T_{1} - T_{o} \cdot T_{x} \right] - w_{1} \cdot b_{2} \cdot a_{o} \cdot (T_{o} - T_{1})

W celu sprowadzenia błędu Δ(T1)\Delta(T_{1}) do zera niezbędna jest wartość w1w_{1}:

w1=b1⋅B⋅[T12+(Tx−To)⋅T1−To⋅Tx]b2⋅ao⋅(To−T1)(45)w_{1} = \frac{b_{1} \cdot B \cdot \left[ T_{1}^{2} + (T_{x} - T_{o}) \cdot T_{1} - T_{o} \cdot T_{x} \right]}{b_{2} \cdot a_{o} \cdot (T_{o} - T_{1})} \qquad (45)

Wprowadzając otrzymaną wartość w1w_{1} do (45), otrzymujemy zależność opisującą błąd temperaturowy skompensowanego przetwornika

Δ(T)=b1⋅B⋅(T−To)⋅(T−T1)(46)\Delta(T) = b_{1} \cdot B \cdot (T - T_{o}) \cdot (T - T_{1}) \qquad (46)

Narzućmy warunek dodatkowy żądając, aby wartości błędu (73) na krańcach TaT_{a}, TbT_{b} zakresu temperatur pracy były sobie równe oraz aby były równe - z przeciwnym znakiem - wartości błędu w temperaturze TmT_{m}, w której występuje ekstrmum funkcji /rys. 19a./:

Δ(Ta)=Δ(Tb)=−Δ(Tm)(47)\Delta(T_{a}) = \Delta(T_{b}) = -\Delta(T_{m}) \qquad (47)

łatwo więc można wykazać, że

Tm=To+T12(48)T_{m} = \frac{T_{o} + T_{1}}{2} \qquad (48)

oraz

Δ(Tm)=−b1⋅B⋅(To−T1)24(49)\Delta(T_{m}) = - \frac{b_{1} \cdot B \cdot (T_{o} - T_{1})^{2}}{4} \qquad (49)

Zakładając Ta<TbT_{a} < T_{b}, To<T1T_{o} < T_{1}, wartość temperatury odniesienia ToT_{o} będzie określona zależnością

s. 60

60

To=Ta+2−122⋅(Tb−Ta)(50)T_{o} = T_{a} + \frac{\sqrt{2} - 1}{2\sqrt{2}} \cdot (T_{b} - T_{a}) \qquad (50)

natomiast temperatura pierwszej kompensacji T1T_{1}:

T1=Tb−2−122⋅(Tb−Ta)(51)T_{1} = T_{b} - \frac{\sqrt{2} - 1}{2\sqrt{2}} \cdot (T_{b} - T_{a}) \qquad (51)

Współczynnik η\eta zmniejszenia błędu temperaturowego czujnika będzie określony następująco:

η=∣δmΔm∣=∣Tb+Txb1⋅(Tb−T1)∣(52)\eta = \left| \frac{\delta m}{\Delta m} \right| = \left| \frac{T_{b} + T_{x}}{b_{1} \cdot (T_{b} - T_{1})} \right| \qquad (52)

Jak widać warunek (47) jednoznacznie uzależnia wartości temperatury odniesienia i pierwszej kompensacji od krańcowych temperatur pracy, niezależnie od parametrów czujnika i układu kompensacji /które powinny jedynie doprowadzać do warunku 47 /.

Wartość maksymalnego błędu temperaturowego przetwornika jest również w takim przypadku niezależna od parametrów układu kompensacji. Wstawiając bowiem (50) i (51) do (49) uzyskujemy wartość ekstremalnego błędu przetwornika Δ(Tm)\Delta(T_{m}):

Δ(Tm)=−b1⋅B8⋅(Tb−Ta)2(53)\Delta(T_{m}) = - \frac{b_{1} \cdot B}{8} \cdot (T_{b} - T_{a})^{2} \qquad (53)

a więc wartość ta zależy jedynie od szerokości zakresu temperatur pracy i nieskompensowanej stałej B charakterystyki błędu czujnika, oraz stałej b1b_{1} układu normalizującego.

Przykład 2

Przetwornik powinien pracować w zakresie temperatur 253K ÷ 333K. Na podstawie (51) i (52) obliczamy: ToT_{o} = 264, 7K, T1T_{1} = 321, 3K. Błąd czujnika odniesiony do temperatury ToT_{o} opiszemy zależnością /wyrażo-

s. 61
Rys. 19. Błąd temperaturowy czujnika δ(T) i przetwornika Δ(T) w przypadku s = 2
Rys. 19. Błąd temperaturowy czujnika δ(T) i przetwornika Δ(T) w przypadku s = 2
s. 62

62

ną w procentach żakresu pomiarowego/:

Δ(T)=−6 10−4⋅(T−264,7)⋅(T−423)\Delta(T) = - 6\ 10^{-4} \cdot (T - 264, 7) \cdot (T - 423)

i dla uproszczenia przyjmiemy parametry układu normalizującego b1=b2=1b_{1} = b_{2} = 1. Wykorzystując (46) otrzymamy zależność dla błędu temperaturowego przetwornika /wyrażonego w procentach zakresu pomiarowego/:

Δ(T)=−6 10−4⋅(T−264,7)⋅(T−321,3)\Delta(T) = -6\ 10^{-4} \cdot (T - 264, 7) \cdot (T - 321, 3)

Ilustrację graficzną obu błędów przedstawiono na rys. 19b. Współczynnik η\eta zmniejszenia błędu ma wartość 7, 7.

Odnosząc błąd do temperatury ToT_{o} stwierdzamy, że błąd przetwornika zawarty jest między prostymi o nachyleniu +0, 41%/10K i 0, 07%/10K /linie przerywane przechodzące przez punkt (283K, 0) na rys. 19b/. Regulując przetwornik tak, aby uzyskać wartość odniesienia rozpatrywanego parametru w temperaturze TmT_{m} = 293K bez zmiany parametrów kompensacji, otrzyma się ograniczenie błędu prostymi o nachyleniu ± 0, 24%/10K /linie przerywane przecinają się w 293K na rys. 19b/. Maksymalne nachylenie błędu czujnika δ(T)\delta(T) /pomiędzy temperaturami 253K i 293 K/ wynosi 1, 05%/10K.

4. 4. 2. Przypadek: s=3s = 3, m=1m = 1, n=1n = 1

Dla tego przypadku błąd czujnika opisany jest zależnością:

δ(T)=B⋅(T−To)⋅(T+Tx)⋅(T+Ty)(54)\delta(T) = B \cdot (T - T_{o}) \cdot (T + T_{x}) \cdot (T + T_{y}) \qquad (54)

gdzie ToT_{o} - temperatura odniesienia B, TxT_{x}, TyT_{y} - stałe

Na podstawie (27) zapisujemy więc

fo=ao⋅(T−To)f_{o} = a_{o} \cdot (T - T_{o})

i uwzględniając wo=0w_{o} = 0 otrzymujemy błąd przetwornika w postaci

Δ(T)=b1⋅B⋅(T−To)⋅(T+Tx)⋅(T+Ty)+b2⋅w1⋅ao⋅(T−To)\Delta(T) = b_{1} \cdot B \cdot (T - T_{o}) \cdot (T + T_{x}) \cdot (T + T_{y}) + b_{2} \cdot w_{1} \cdot a_{o} \cdot (T - T_{o})
s. 63

63

Procedura kompensacji temperaturowej jest następująca:

T=To;Δ(To)=0T = T_{o}; \quad \Delta(T_{o}) = 0
T=T1;Δ(T1)=b1⋅B⋅(T1−To)⋅(T1+Ty)+b2⋅w1⋅ao⋅(T1−To)T = T_{1}; \quad \Delta(T_{1}) = b_{1} \cdot B \cdot (T_{1} - T_{o}) \cdot (T_{1} + T_{y}) + b_{2} \cdot w_{1} \cdot a_{o} \cdot (T_{1} - T_{o})

W celu sprowadzenia Δ(T1)\Delta(T_{1}) do zera wymagana jest wartość w1w_{1}

w1=−b1⋅B⋅(T1+Tx)⋅(T1+Ty)b2⋅ao(55)w_{1} = - \frac{b_{1} \cdot B \cdot (T_{1} + T_{x}) \cdot (T_{1} + T_{y})}{b_{2} \cdot a_{o}} \qquad (55)

Po uwzględnieniu otrzymanej wartości (55) uzyskuje się dla w1w_{1} zależność na błąd skompensowanego przetwornika

Δ(T)=b1⋅B⋅(T−To)⋅[T2+(Tx+Ty)⋅T−T12⋅(Tx+Ty)](56)\Delta(T) = b_{1} \cdot B \cdot (T - T_{o}) \cdot \left[ T^{2} + (T_{x} + T_{y}) \cdot T - T_{1}^{2} \cdot (T_{x} + T_{y}) \right] \qquad (56)

Można założyć, że z punktu widzenia potrzeb kompensacji korzystne będzie, a by zależność (56) charakteryzowała się następującymi cechami: a/ aby miała 3 pierwiastki rzeczywiste; zostanie to spełnione, gdy wyróżnik D równania kwadratowego w nawiasie prostokątnym zależności (56)

Δ=(Tx+Ty+2T1)2\Delta = (T_{x} + T_{y} + 2T_{1})^{2}

będzie większy od zera, a więc dla każdego T1T_{1} spełniającego warunek:

T1≠−Tx+Ty2T_{1} \neq - \frac{T_{x} + T_{y}}{2}

b/ aby oba pierwiastki równania w nawiasie prostokątnym zależności (56) leżały w zakresie temperatur pracy: warunku tego nie można jednak zagwarantować przy jednym warunku swobody w postaci stałej T1T_{1}.

Minimalizację błędu (56) przy założonych rozrzutach parametrów TxT_{x} i TyT_{y} najszybciej można przeprowadzić za pomocą maszynowych obliczeń wariantów na podstawie wyników pomiaru błędu czujnika.

Przykład 3

Zakres temperatury pracy od 253K do 333K. Błąd czujnika

s. 64
Rys. 20. Błąd temperaturowy czujnika δ(T) oraz przetwornika Δ(T) w przypadku s = 3, m = 1, n = 1.
Rys. 20. Błąd temperaturowy czujnika δ(T) oraz przetwornika Δ(T) w przypadku s = 3, m = 1, n = 1.
s. 65
δ(T)=−1,8⋅10−5⋅(T−263)(T−238,3)⋅(T−382,7)\delta(T) = -1, 8 \cdot 10^{-5} \cdot (T - 263)(T - 238, 3) \cdot (T - 382, 7)

Stałe b1=b2=1b_{1} = b_{2} = 1

Ponieważ TxT_{x} = 238, 3K, TyT_{y} = 382, 7K więc na podstawie (56) zapiszemy

Δ(T)=−1,8⋅10−5⋅(T−263)⋅(T2−621⋅T+621⋅T12)\Delta(T) = - 1, 8 \cdot 10^{-5} \cdot (T - 263) \cdot (T^{2} - 621 \cdot T + 621 \cdot T_{1}^{2})

Na rys. 20a przedstawiono zależność błędu czujnika od temperatury, natomiast na rys. 20b - zależność błędu przetwornika od temperatury dla trzech wartości T1T_{1}: 288K, 303K, 318K.

Należy zaznaczyć, że istnieje możliwość zmiany wartości ToT_{o} /tj. przesuwania punktu odniesienia po krzywej 20a/ bez zmiany kształtu krzywej, tak aby minimalizować błąd przetwornika.

Z przedstawionego przykładu wynika, że jakość kompensacji jest niska.

4. 4. 3. Przypadek: s = 3, m = 1, n = 2

Błąd temperaturowy czujnika opisany jest zależnością (54). Na podstawie (27) zapiszemy

fo=ao⋅(T−To)f_{o} = a_{o} \cdot (T - T_{o})
f1=a1⋅(T−T1)f_{1} = a_{1} \cdot (T - T_{1})

i po uwzględnieniu, że wo=0w_{o} = 0 błąd przetwornika przyjmie postać

Δ(T)=b1⋅B(T−To)⋅(T+Tx)⋅(T+Ty)+b2⋅w1⋅ao⋅(T−To)+\Delta(T) = b_{1} \cdot B (T - T_{o}) \cdot (T + T_{x}) \cdot (T + T_{y}) + b_{2} \cdot w_{1} \cdot a_{o} \cdot (T - T_{o}) +
+b2⋅w2⋅ao⋅a1(T−To)⋅(T−T1)+ b_{2} \cdot w_{2} \cdot a_{o} \cdot a_{1} (T - T_{o}) \cdot (T - T_{1})

Procedura kompensacji temperaturowej jest następująca:

T=To;Δ(To)=0,wo=0T = T_{o}; \quad \Delta(T_{o}) = 0, \qquad w_{o} = 0
T=T1;Δ(T1)=b1⋅B(T1−To)⋅(T1+Tx)⋅(T1+Ty)+b2⋅w1⋅ao⋅(T1−To)+T = T_{1}; \quad \Delta(T_{1}) = b_{1} \cdot B (T_{1} - T_{o}) \cdot (T_{1} + T_{x}) \cdot (T_{1} + T_{y}) + b_{2} \cdot w_{1} \cdot a_{o} \cdot (T_{1} - T_{o}) +
+b2⋅w2⋅ao⋅a1⋅(T1−To)+ b_{2} \cdot w_{2} \cdot a_{o} \cdot a_{1} \cdot (T_{1} - T_{o})
s. 66

Wartość w1w_{1} sprowadzającą błąd Δ(T1)\Delta(T_{1}) do zera zapiszemy

w1=−b1⋅B(T1+Tx)⋅(T1+Ty)b2⋅ao(57)w_{1} = - \frac{b_{1} \cdot B (T_{1} + T_{x}) \cdot (T_{1} + T_{y})}{b_{2} \cdot a_{o}} \qquad (57)

Po uwzględnieniu znalezionej wartości w1w_{1} wyrażenie na błąd przetwornika przyjmie formę

Δ(T)=b1⋅B(T−To)⋅(T+Tx)⋅(T+Ty)−b1⋅B(T1+Tx)⋅(T1+Ty)⋅(T−To)+\Delta(T) = b_{1} \cdot B (T - T_{o}) \cdot (T + T_{x}) \cdot (T + T_{y}) - b_{1} \cdot B (T_{1} + T_{x}) \cdot (T_{1} + T_{y}) \cdot (T - T_{o}) +
+b2⋅w2⋅ao⋅a1⋅(T−To)⋅(T−T1)+ b_{2} \cdot w_{2} \cdot a_{o} \cdot a_{1} \cdot (T - T_{o}) \cdot (T - T_{1})

T=T2T = T_{2}; określając błąd analogicznie jak wyżej, otrzymujemy zerującą go wartość w2w_{2}:

w2=−b1⋅B(T2+Tx)⋅(T2+Ty)−b1⋅B(T1+Tx)⋅(T1+Ty)b2⋅ao⋅a1⋅(T2−T1)(58)w_{2} = - \frac{b_{1} \cdot B (T_{2} + T_{x}) \cdot (T_{2} + T_{y}) - b_{1} \cdot B (T_{1} + T_{x}) \cdot (T_{1} + T_{y})}{b_{2} \cdot a_{o} \cdot a_{1} \cdot (T_{2} - T_{1})} \qquad (58)

Po wstawieniu znalezionych wartości w1w_{1} i w2w_{2} otrzymamy zależność opisującą błąd temperaturowy skompensowanego przetwornika:

Δ(T)=b1⋅B(T−To)⋅(T+T1)⋅(T+T2)(59)\Delta(T) = b_{1} \cdot B (T - T_{o}) \cdot (T + T_{1}) \cdot (T + T_{2}) \qquad (59)

[?:likely (T−T1)⋅(T−T2)(T - T_{1}) \cdot (T - T_{2}); "+" signs as printed]

Błąd ten jest uzależniony od wartości nieskompensowanej stałej B błędu czujnika, od wartości temperatur kompensacji T1T_{1} i T2T_{2} oraz odniesienia ToT_{o} i jednocześnie zmodyfikowany przez współczynnik b1b_{1} układu normalizującego. Zależność (59) gwarantuje uzyskanie trzech zer w zakresie temperatur pracy /rys. 21/. Narzucając warunek równych pod względem modułu wartości błędu Δ(T)\Delta(T) w ekstremach występujących w temperaturach TcT_{c}, TdT_{d}, uzyska się powiązanie temperatury odniesienia z temperaturami kompensacji:

To=T1+T22(60)T_{o} = \frac{T_{1} + T_{2}}{2} \qquad (60)

Wartość maksymalnego błędu można kształtować przez zmianę różnic między temperaturami kompensacji i zmianę wartości temperatury odniesienia.

s. 67
Rys. 21. Pożądany przebieg błędu przetwornika w przypadku s = 3, m = 1, n = 2
Rys. 21. Pożądany przebieg błędu przetwornika w przypadku s = 3, m = 1, n = 2

Przykład 4

Zakres temperatur pracy od 253K do 333K Błąd czujnika opisano krzywą /rys. 22/

δ(T)=−1,8⋅10−5⋅(T−293)⋅(T−373)⋅(T−218)\delta(T) = - 1,8 \cdot 10^{-5} \cdot (T - 293) \cdot (T - 373) \cdot (T - 218)

stałe b1=b2=1b_{1} = b_{2} = 1

A więc, Tx=−373KT_{x} = - 373K, Ty=−218KT_{y} = - 218K, B=−1,8⋅10−5B = - 1,8 \cdot 10^{-5}

Zależność błędu /59/ skompensowanego przetwornika od temperatury jest uzależniona od doboru temperatur kompensacji T1T_{1} i T2T_{2}:

Δ(T)=−1,8⋅10−5⋅(T−293)(T+T1)⋅(T+T2)\Delta(T) = - 1,8 \cdot 10^{-5} \cdot (T - 293)(T + T_{1}) \cdot (T + T_{2})

Wykresy tego błędu dla różnych par T1T_{1}, T2T_{2} spełniających (60) oraz jednej pary nie spełniającej (60) przedstawiono na rys. 22.b.

s. 68
Rys. 22. Błąd temperaturowy czujnika: δ(T) oraz przetwornika Δ(T) przypadku s = 3, m = 1, n = 2
Rys. 22. Błąd temperaturowy czujnika: δ(T) oraz przetwornika Δ(T) przypadku s = 3, m = 1, n = 2
s. 69

4. 5. Podsumowanie

Zaproponowana metoda kompensacji błędu temperaturowego wartości początkowej i końcowej zakresu pomiarowego przetwornika ciśnienia polega na sumowaniu wartości rozpatrywanego parametru z wartością sygnału generatora sygnałów kompensujących, tak aby w n temperaturach uzyskać zrównanie się wartości tego parametru z jego wartością w temperaturze odniesienia. Liczba n temperatur kompensacji jest uzależniona od wartości błędu czujnika i wartości dopuszczalnego błędu przetwornika.

Jeśli liczba n temperatur kompensacji będzie równa stopniowi wielomianu opisującego błąd temperaturowy danego parametru charakterystyki czujnika, a sygnały generatorów kompensujących będą uzależnione liniowo od temperatury, wtedy przedstawiona procedura kompensacji jednoznacznie doprowadzi do tożsamościowego zerowania się błędu przetwornika w całym rozpatrywanym zakresie temperatur pracy, niezależnie od doboru wartości temperatur kompensacji.

Jeśli znana jest wstępnie pewna liczba własności charakterystyki błędu temperaturowego czujnika, fakt ten można wykorzystać do zmniejszenia liczby temperatur niezbędnych do tożsamościowego sprowadzenia błędu przetwornika do zera.

Jeśli nie jest wymagana pełna kompensacja błędu, proponuje się stosowanie metody z liczbą n temperatur kompensacji mniejszą o jedność od rzędu wielomianu opisującego błąd temperaturowy czujnika, przy liniowym uzależnieniu od temperatury sygnałów generatorów kompensujących.

Procedura kompensacji danego parametru charakteryzuje się niezależnością jednokierunkową, co oznacza, że ustalanie warunków kompensacji tego parametru w bieżącej temperaturze kompensacji nie zmienia warunków w poprzednich temperaturach kompensacji. Dzięki temu eliminuje się konieczność iteracyjnego wyznaczania co wymaga kilku cykli temperaturowych lub obliczania /najczęściej za pomocą maszyny cyfrowej/ wartości elementów układu kompensacji niezbędnych do skompensowania błędu temperaturowego rozpatrywanego para-

s. 70

metru charakterystyki przetwornika.

Porównując przedstawioną metodę z innymi znanymi metodami [1, 32 ÷ 76], należy stwierdzić, że jedynie przedstawiona w pracy metoda umożliwia tożsamościowe sprowadzenie błędu temperaturowego przetwornika do zera oraz jednokierunkową niezależność ustalania warunków kompensacji.

W przedstawionej metodzie zmiany charaketrystyki temperaturowej sygnału kompensującego dokonuje się za pomocą zmiany współczynników proporcjonalności niezależnych od temperatury. W zastosowaniu praktycznym metoda kompensacji może polegać na odpowiednim ustawieniu dzielników potencjometrycznych sygnałów kompensujących. Ponieważ na dzielnikach takich występowałby wówczas tylko sygnał rzędu błędów temperaturowych czujnika, a więc mały procent sygnału zakresowego, można by złagodzić wymagania metrologiczne stawiane tym elementom.

Zastosowanie wspomnianych dzielników potencjometrycznych eliminowałoby konieczność magazynowania dużej liczby wartości elementów kompensujących /które potencjalnie mogą być użyte/ lub wykonywanie ich po uzyskaniu informacji o błędach czujników.

Prezentowana metoda może być stosowana we wszelkiego typu przetwornikach pomiarowych wielkości nieelektrycznych spełniających warunki wymienione w p. 4. 1.

Na rys. 23 przedstawiono schematycznie przebieg błędu temperaturowego skompensowanego przetwornika w przypadku kompensacji pełnej i kompensacji niepełnej dla najczęściej spotykanych w praktyce przepadków.

s. 71
Rys. 23. Schematyczny przebieg błędu temperaturowego Δ(T) skompensowanego przetwornika w przypadku kompensacji pełnej i niepełnej; a. błąd czujnika δ(T) opisany funkcją liniową; b. błąd czujnika δ(T) opisany funkcją kwadratową; c. błąd czujnika δ(T) opisany funkcją sześcienną
Rys. 23. Schematyczny przebieg błędu temperaturowego Δ(T) skompensowanego przetwornika w przypadku kompensacji pełnej i niepełnej; a. błąd czujnika δ(T) opisany funkcją liniową; b. błąd czujnika δ(T) opisany funkcją kwadratową; c. błąd czujnika δ(T) opisany funkcją sześcienną