428 / 275 / 52
Регистрация: 24.09.2017
Сообщений: 1,623
Записей в блоге: 2
1

ГОСТ Р 8.585-2001 ГСИ. Термопары

02.04.2021, 10:21. Показов 1254. Ответов 30
Метки нет (Все метки)

В ГОСТ Р 8.585-2001 ГСИ отсутствуют обратные полиномы для термопар типа "К" от минус 270 до минус 200 град. и для термопар типа "В" от нуля до плюс 250 град. Так должно и быть, или это ошибка в ГОСТ?
0
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
02.04.2021, 10:21
Ответы с готовыми решениями:

Обмен с электросчетчиком CE301, СЕ303. ГОСТ Р МЭК 61107-2001
Здравствуйте! Примеров работы по Modbus полно, а вот по МЭК 61107-2001 нигде не встречал ни одного...

Не работает шифрование по ГОСТ 28147-89, ГОСТ Р 34.10-2001,ГОСТ Р 34.11-94 (ASP.NET)
Создавал веб-приложение(ASP.Net) для шифрования по данным гостам, нашёл по данной теме код, но он...

ГОСТ р 34.10-2001
Кто нибудь реализовывал алгоритм ЦП по алгоритму из этого госта?

Реализация ГОСТ Р 34.10-2001
Нужно реализовать стандар Электроной цифровой подписи ГОСТ Р 34.10-2001 на C# (windows form...

30
Модератор
Эксперт по электронике
7651 / 3815 / 1483
Регистрация: 01.02.2015
Сообщений: 11,813
Записей в блоге: 2
02.04.2021, 18:19 2
Смотрите внимательно - всё есть.
Найдите другой скан ГОСТ

ГОСТ Р 8.585-2001 ГСИ. Термопары


ГОСТ Р 8.585-2001 ГСИ. Термопары
0
428 / 275 / 52
Регистрация: 24.09.2017
Сообщений: 1,623
Записей в блоге: 2
02.04.2021, 20:57  [ТС] 3
Цитата Сообщение от ФедосеевПавел Посмотреть сообщение
Смотрите внимательно - всё есть.
Цитата Сообщение от Argus19 Посмотреть сообщение
отсутствуют обратные полиномы

Обратный полином для термопар типа "М" во всех источниках указан от 0 до 5, а везде только от 0 до 3, поэтому считает неправильно.
0
Модератор
Эксперт по электронике
7651 / 3815 / 1483
Регистрация: 01.02.2015
Сообщений: 11,813
Записей в блоге: 2
02.04.2021, 22:12 4
Посмотрите и сравните другие источники - предыдущий ГОСТ, аналог ГОСТ - ICO, т.к. термопары не только в России применяются.
0
428 / 275 / 52
Регистрация: 24.09.2017
Сообщений: 1,623
Записей в блоге: 2
02.04.2021, 23:00  [ТС] 5
ФедосеевПавел, Я нашёл только старый ГОСТ 3044-84. Там только прямые полиномы. Поисковик не находит ISO.
0
Модератор
Эксперт по электронике
7651 / 3815 / 1483
Регистрация: 01.02.2015
Сообщений: 11,813
Записей в блоге: 2
03.04.2021, 12:31 6
В ГОСТ Р 8.585-2001
Предисловие
......................
3. ВЗАМЕН ГОСТ Р 50431-92. МИ 2559-99
В https://ru.wikipedia.org/wiki/... 1%80%D0%B0
Точный состав сплава термоэлектродов для термопар из неблагородных металлов в МЭК 60584-1 не приводится. НСХ для хромель-копелевых термопар ТХК и вольфрам-рениевых термопар определены только в ГОСТ Р 8.585-2001. В стандарте МЭК данные термопары отсутствуют. По этой причине характеристики импортных датчиков из этих металлов могут существенно отличаться от отечественных, например импортный Тип L и отечественный ТХК не взаимозаменяемы. При этом, как правило, импортное оборудование не рассчитано на отечественный стандарт.

В настоящее время стандарт МЭК 60584 пересматривается. Планируется введение в стандарт вольфрам-рениевых термопар типа А-1, НСХ для которых будет соответствовать российскому стандарту, и типа С по стандарту АСТМ[7].

В 2008 г. МЭК ввел два новых типа термопар: золото-платиновые и платино-палладиевые. Новый стандарт МЭК 62460 устанавливает стандартные таблицы для этих термопар из чистых металлов. Аналогичный Российский стандарт пока отсутствует.
Плохой из вас шпион - не узнаете из открытых источников особо секретные данные.
0
428 / 275 / 52
Регистрация: 24.09.2017
Сообщений: 1,623
Записей в блоге: 2
03.04.2021, 18:41  [ТС] 7
Цитата Сообщение от ФедосеевПавел Посмотреть сообщение
Плохой из вас шпион - не узнаете из открытых источников особо секретные данные.
Конечно плохой.
https://srdata.nist.gov/its90/... lcoeff.tab
Там нет того, что мне нужно.
Для термопар типа "М" совсем ничего.
Придётся спрашивать знакомых киповцев. Может, у них что-то есть.
0
Модератор
Эксперт по электронике
7651 / 3815 / 1483
Регистрация: 01.02.2015
Сообщений: 11,813
Записей в блоге: 2
03.04.2021, 20:36 8
Не могу объяснить причину отсутствия некоторых диапазонов в обратных функциях, при наличии этих диапазонов в "прямых" функциях.
Но если составители ГОСТ и IEC посчитали возможным исключить (пропустить) диапазоны, значит тому была причина.

Раз их нет, значит решение не в поиске чужих сомнительных аппроксимаций. Я бы тогда попробовал решить уравнение любым методом - хорд, деления пополам, Ньютона, Золотого сечения или что придумается.

Т.е. для термопары типа В имеются функции:
- E(t) для 2 диапазонов
1) 0.000 ... 630.615 (°C)
2) 630.615 ... 1820.000 (°C)
- t(E) для 2 диапазонов
1) 250. ... 700. (°C)
2) 700. ... 1800. (°C)
Т.е. для обратной функции пропущены диапазоны
1) 0. ... 250. (°C)
2) 1800. ... 1820. (°C)

Для определения температуры в диапазоне 0. ... 250. (°C) от термоЭДС в диапазоне 0.000 ... 0.291 мВ нужно решить уравнение E(t)=Eизмеренное
где E(t) - "прямая" функция для диапазона 0.000 ... 630.615 (°C)

Оценивая значения термоЭДС для отсутствующего диапазона 0...250 (°C) - 0.000 ... 0.291 мВ, вижу, что там очень ничтожные величины, которые на практике в промышленности уже измерять не будут, т.к. они затеряются в шумах и помехах.
Другими словами, уравнение можно решить, но практического значения оно уже иметь не будет.

Добавлено через 20 минут
Для обратного аппроксимирующего полинома термопары типа M над знаком сигмы число 5, а приведено 4 коэффициента с индексами от 0 до 3. Вероятно, это опечатка - над сигмой должно быть число 3.

Добавлено через 42 минуты
Для термопары типа M (ТМК) в скане ГОСТ опечатка.
Тестовая программа
Pascal
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
{
Вычисление термоЭДС для термопары типа M (ТМК)
по полиному из
ГОСТ Р 8.585-2001 ГСИ.
Термопары. Номинальные статические характеристики преобразования
}
program type_M;
 
const
  E20  = 0.873; //термоЭДС типа M (ТМК) при 20 градусах
  E100 = 4.722; //термоЭДС типа M (ТМК) при 100 градусах
 
type
  TArrayTE = array[0..3] of real;
  TArrayET = array[0..3] of real;
const
  //для диапазона температур от минус 200 °С до плюс 100 °С
  A: TArrayTE = (2.4455560e-6, 4.2638917e-2, 5.0348392e-5, -4.4974485e-8);
  //для диапазона температур от минус 200 °С до плюс 100 °С
  C: TArrayET = (0.4548090, 2.2657698e-2, 7.7935652e-7, 1.1786931e-10);
  //зависимость термоЭДС от температуры
  function E_typeM(t: real): real;
  var
    i: integer;
    r: real;
  begin
    r := 0;
    for i := high(A) downto low(A) do
      r := r * t + A[i];
    E_typeM := r;
  end;
  //обратная зависимость температуры от термоЭДС
  function T_typeM(E: real): real;
  var
    i: integer;
    r: real;
  begin
    r := 0;
    for i := high(C) downto low(C) do
      r := r * E + C[i];
    T_typeM := r;
  end;
 
begin
  writeln(E_typeM(20): 0: 3, ' = ', E20: 0: 3);
  writeln(E_typeM(100): 0: 3, ' = ', E100: 0: 3);
  writeln(T_typeM(E20): 0: 3, ' = ', 20.0: 0: 3);
  writeln(T_typeM(E100): 0: 3, ' = ', 100.0: 0: 3);
end.
Прошу прощения, что выводятся числа без пояснительных сообщений, но это тестовая программа для проверки "прямого" и "обратного" полиномов в двух точках.
Прогон0.873 = 0.873
4.722 = 4.722
0.475 = 20.000
0.562 = 100.000
Как видно из прогона - совпадают расчётные и табличные значения.
1
428 / 275 / 52
Регистрация: 24.09.2017
Сообщений: 1,623
Записей в блоге: 2
03.04.2021, 22:41  [ТС] 9
Цитата Сообщение от ФедосеевПавел Посмотреть сообщение
Для термопары типа M (ТМК) в скане ГОСТ опечатка.
У меня почему-то не получается. Температура в ТЭДС считается нормально, а ТЭДС в температуру даёт почему-то значения не более 0, 56 град. Странно.

Добавлено через 42 минуты
Проверил обратный полином на VB 6.0 тоже ерунда получается.
0
Модератор
Эксперт по электронике
7651 / 3815 / 1483
Регистрация: 01.02.2015
Сообщений: 11,813
Записей в блоге: 2
03.04.2021, 23:53 10
Видимо - ГОСТ небезупречен.
У меня обратный полином для типа M тоже даёт значения, не соответствующие таблице.
0
428 / 275 / 52
Регистрация: 24.09.2017
Сообщений: 1,623
Записей в блоге: 2
04.04.2021, 00:55  [ТС] 11
Цитата Сообщение от ФедосеевПавел Посмотреть сообщение
У меня обратный полином для типа M тоже даёт значения, не соответствующие таблице.
В чём и проблема.
В зарубежном аналоге тоже не хватает обратных полиномов. Ладно. Пусть, так и задумано, но с типом "М" явная ошибка и в зарубежном этого типа совсем нет. Я задал вопрос ещё в двух местах. Но ответ жду не ранее вечера понедельника.
0
78 / 81 / 6
Регистрация: 30.10.2014
Сообщений: 247
04.04.2021, 14:16 12
Эксель умеет находить полином по графику.
Как то баловался с полиномом для В.
1
Вложения
Тип файла: xlsx полином.xlsx (18.2 Кб, 7 просмотров)
428 / 275 / 52
Регистрация: 24.09.2017
Сообщений: 1,623
Записей в блоге: 2
04.04.2021, 14:29  [ТС] 13
Цитата Сообщение от Олег_ Посмотреть сообщение
Как то баловался с полиномом для В
Я не силён в VBA.
Это прямой полином температура в ТЭДС?
0
78 / 81 / 6
Регистрация: 30.10.2014
Сообщений: 247
04.04.2021, 14:35 14
Там нет VBA, в свойствах графика можно добавить линию с полиномом
0
Миниатюры
ГОСТ Р 8.585-2001 ГСИ. Термопары  
428 / 275 / 52
Регистрация: 24.09.2017
Сообщений: 1,623
Записей в блоге: 2
05.04.2021, 00:13  [ТС] 15
Цитата Сообщение от Олег_ Посмотреть сообщение
свойствах графика можно добавить линию с полиномом
Ввёл данные из ГОСТ на тип "М". Получил полиномы 3 и 5 степеней. Оба привирают.
0
78 / 81 / 6
Регистрация: 30.10.2014
Сообщений: 247
05.04.2021, 08:07 16
Не обязательно расчёт вести через полином, можно по табличным данным линейно кусочную интерполяцияю сделать.
Громоздко, но зато никаких расхождений с гостовскими таблицами.
0
Модератор
Эксперт по электронике
7651 / 3815 / 1483
Регистрация: 01.02.2015
Сообщений: 11,813
Записей в блоге: 2
05.04.2021, 08:59 17
Или же, т.к. один из полиномов для типа M корректно вычисляет - обратную для него функцию не выводить в виде готовой формулы, а значение находить численным методом, решая нелинейное уравнение методом дихотомии (или Ньютона) - там буквально 10-15 итераций, т.к. диапазон измерений для типа M составляет (-200...+100).

Зачем-то же нас учили чисметам?!. Пришла пора применять.

Добавлено через 36 минут
QBasic/QuickBASIC
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
'--------------------------------------------------------------------------
'Вычисление термоЭДС для термопары типа M (ТМК)
'по полиному из
'ГОСТ Р 8.585-2001 ГСИ.
'Термопары. Номинальные статические характеристики преобразования
'
'С учётом поправки ИУС 5-2003.
'Вычисления разбиты на диапазоны:
'1) аппроксимирующий полином E(t) - зависимость термоЭДС от температуры
'-200...+100 °С
'2) обратный аппроксимирующий полином t(E) - зависимость температуры от термоЭДС
'-200...+100 °С
'--------------------------------------------------------------------------
 
FUNCTION Polynomial( x AS DOUBLE, a() AS DOUBLE) AS DOUBLE
        DIM AS DOUBLE r
        r=0
        FOR i AS INTEGER = UBOUND(a) TO LBOUND(a) STEP -1
                r = r*x+a(i)
        NEXT i
        RETURN r
END FUNCTION
 
'зависимость термоЭДС от температуры
FUNCTION E_typeM (T90 AS DOUBLE) AS DOUBLE
        'коэффициенты полинома для вычисления значения термоЭДС по значению
        'температуры рабочего спая
        'для диапазона температур -200...+100 °С
        DIM A(3) AS DOUBLE={ 2.4455560e-6, 4.2638917e-2, 5.0348392e-5, -4.4974485e-8}
        RETURN Polynomial(T90, A())
END FUNCTION
'обратная зависимость температуры от термоЭДС
FUNCTION T_typeM (EMF AS DOUBLE) AS DOUBLE
        'коэффициенты полинома для вычисления температуры рабочего спая по значению
        'термоЭДС
        'для диапазона температур -200...+100 °С
        DIM C(3) AS DOUBLE ={0.4548090, 2.2657698e-2, 7.7935652e-7, 1.1786931e-10}
        RETURN Polynomial(EMF, C())
END FUNCTION
'--------------------------------------------------------------------------
'Т.к. функция из ГОСТ Р 8.585-2001 для термопары типа M вероятно опубликована
'с ошибкой, то измеренную температуру по ТЭДС находим численным методом
'деления отрезка пополам
FUNCTION T_typeM_nm (EMF AS DOUBLE) AS DOUBLE
        DIM AS DOUBLE left_bound, right_bound, middle
        CONST   eps=0.001
        left_bound=-200 'левая граница измеряемого термопарой типа M диапазона
        right_bound=100 'правая граница измеряемого термопарой типа M диапазона
        DO
                middle=left_bound+(right_bound-left_bound)/2
                IF E_typeM(middle)>EMF _
                        THEN right_bound=middle _
                        ELSE left_bound=middle
        LOOP WHILE (right_bound-left_bound)>=eps
        RETURN left_bound+(right_bound-left_bound)/2
END FUNCTION
'--------------------------------------------------------------------------
'проверка вычисления термоЭДС и сравнение с известными константами
SUB test(T90 AS DOUBLE, EMF AS DOUBLE)
        PRINT USING "test t=####, EMF=###.###";T90;EMF
        PRINT "         calc    exact"
        'проверка вычисления термоЭДС и сравнение с известными константами
        PRINT  USING " E(t):  ####.### = ####.###";E_typeM(T90);EMF
        'проверка вычисления температуры и сравнение с известными константами
        PRINT  USING " t(E):  ####.### = ####.###";T_typeM(EMF);T90
        PRINT  USING " t_nm(E):####.### = ####.###";T_typeM_nm(EMF);T90
 
END SUB
'--------------------------------------------------------------------------
'проверка вычисления термоЭДС и сравнение с известными константами
CONST E0  = 0.000 'термоЭДС типа M (ТМК) при 0 °С
CONST E20  = 0.873 'термоЭДС типа M (ТМК) при 20 °С
CONST E100 = 4.722 'термоЭДС типа M (ТМК) при 100 °С
 
'проверка вычисления температуры и сравнение с известными константами
test(0.0, E0)
test(20.0, E20)
test(100.0, E100)
 
'ожидание нажатия любой клавиши
'do
'loop while inkey=""
 
END
Прогонtest t= 0 EMF= 0.000
calc exact
E(t): 0.000 = 0.000
t(E): 0.455 = 0.000
tnm(E): 0.000 = 0.000
test t= 20 EMF= 0.873
calc exact
E(t): 0.873 = 0.873
t(E): 0.475 = 20.000
tnm(E): 20.010 = 20.000
test t= 100 EMF= 4.722
calc exact
E(t): 4.722 = 4.722
t(E): 0.562 = 100.000
tnm(E): 99.992 = 100.000
Теперь уже точно видно, что температура по ТЭДС численным методом даёт приемлемый результат.
Тем более, что на практике температуру точнее 0,1°С даже не отображают, а реально измеренное значение интересно вообще до целых разрядов.
0
428 / 275 / 52
Регистрация: 24.09.2017
Сообщений: 1,623
Записей в блоге: 2
05.04.2021, 09:49  [ТС] 18
Пока можно сказать, что ошибка в ГОСТ в степенях, если в строке №37 записать так:
QBasic/QuickBASIC
1
DIM C(3) AS DOUBLE ={0.4548090, 2.2657698e1, 7.7935652e-1, 1.1786931e-1}
, то результат будет ближе.
Excell показывает так:
C(0) = 0.4407
C(1) = 22.647
C(2) = -0.77622
C(3) = 0.1189
0
Модератор
Эксперт по электронике
7651 / 3815 / 1483
Регистрация: 01.02.2015
Сообщений: 11,813
Записей в блоге: 2
05.04.2021, 13:02 19
Проверьте график - я как-то аппроксимировал или интерполировал полиномом Ньютона расходную характеристику диафрагмы (корнеизвлекающую) - хотел не кусочно-линейной функцией реализовывать, а простенькой функцией. В узлах было точное совпадение, а между узлами - "синусы" с огромными полуволнами.
Поэтому думаю, что составление этих полиномов - ручная работа, или близкая к ручной.

Добавлено через 1 минуту
Хотя для педантичности - искал бы численными методами. Так и никакой отсебятины и ГОСТ привлечён.
0
428 / 275 / 52
Регистрация: 24.09.2017
Сообщений: 1,623
Записей в блоге: 2
05.04.2021, 13:17  [ТС] 20
Цитата Сообщение от ФедосеевПавел Посмотреть сообщение
Проверьте график
Вот:
0
Вложения
Тип файла: xlsx Книга1.xlsx (20.4 Кб, 5 просмотров)
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
05.04.2021, 13:17

Помощь в написании контрольных, курсовых и дипломных работ здесь.

Реализация ГОСТ Р34.11-2001
Всем привет, пытаюсь реализовать ГОСТ р34.11-2001, дошел до проверки ЭЦП, P = u1G + u2Ha, u1, u2...

Как реализовать ГОСТ Р34.11-2001?
Всем привет, пытаюсь реализовать ГОСТ р34.11-2001, дошел до проверки ЭЦП, P = u1G + u2Ha, u1, u2...

Реализация алгоритма ЭЦП ГОСТ 34.10-2001 (или 34.10-1994) на C++
Где можно взять реализацию алгоритма ЭЦП ГОСТ 34.10-2001 (или 34.10-1994) на C++ или может есть у...

Моделирование скоростной термопары при помощи одной малоинерционной термопары
Доброго дня, уважаемые специалисты. У меня возникла необходимость программно получить сигнал от...

Проблема с select * from contest Where id_contest='6' and day_contest between #29/07/2001# and #04/08/2001#
Подскажите в чём ошибка! Есть таблица с конкурсными фотографиями из неё должны демонстрироваться...

Lenovo z 585 не включается
Здраствуйте прошу помощи сил уже нет. После синего экрана которой видавал ошибку 0x000000a5 ....


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
20
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.