KR20080072855A - 가상 온도 측정 포인트 - Google Patents

가상 온도 측정 포인트 Download PDF

Info

Publication number
KR20080072855A
KR20080072855A KR1020087011939A KR20087011939A KR20080072855A KR 20080072855 A KR20080072855 A KR 20080072855A KR 1020087011939 A KR1020087011939 A KR 1020087011939A KR 20087011939 A KR20087011939 A KR 20087011939A KR 20080072855 A KR20080072855 A KR 20080072855A
Authority
KR
South Korea
Prior art keywords
temperature
layer
average
wall
determination method
Prior art date
Application number
KR1020087011939A
Other languages
English (en)
Other versions
KR101334885B1 (ko
Inventor
올트리히 차비스카
Original Assignee
지멘스 악티엔게젤샤프트
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 지멘스 악티엔게젤샤프트 filed Critical 지멘스 악티엔게젤샤프트
Publication of KR20080072855A publication Critical patent/KR20080072855A/ko
Application granted granted Critical
Publication of KR101334885B1 publication Critical patent/KR101334885B1/ko

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K19/00—Testing or calibrating calorimeters
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K7/00—Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements
    • G01K7/42—Circuits effecting compensation of thermal inertia; Circuits for predicting the stationary value of a temperature
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F27—FURNACES; KILNS; OVENS; RETORTS
    • F27D—DETAILS OR ACCESSORIES OF FURNACES, KILNS, OVENS, OR RETORTS, IN SO FAR AS THEY ARE OF KINDS OCCURRING IN MORE THAN ONE KIND OF FURNACE
    • F27D21/00—Arrangements of monitoring devices; Arrangements of safety devices
    • F27D21/0014—Devices for monitoring temperature
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K17/00—Measuring quantity of heat
    • G01K17/06—Measuring quantity of heat conveyed by flowing media, e.g. in heating systems e.g. the quantity of heat in a transporting medium, delivered to or consumed in an expenditure device
    • G01K17/08—Measuring quantity of heat conveyed by flowing media, e.g. in heating systems e.g. the quantity of heat in a transporting medium, delivered to or consumed in an expenditure device based upon measurement of temperature difference or of a temperature
    • G01K17/20—Measuring quantity of heat conveyed by flowing media, e.g. in heating systems e.g. the quantity of heat in a transporting medium, delivered to or consumed in an expenditure device based upon measurement of temperature difference or of a temperature across a radiating surface, combined with ascertainment of the heat transmission coefficient

Abstract

본 발명은 두꺼운 벽 또는 샤프트내 온도 프로파일 및 적분 평균 온도 및/또는 축 온도를 결정하기 위한 방법에 관한 것이다. 다층 모델에서 가열 또는 냉각 처리 동안 평균 적분 벽 온도를 결정하기 위하여, 평균 적분 벽 온도는 각각의 층의 평균 온도로부터 계산된다. 다층 모델은 가열 또는 냉각 처리들 동안 평균 적분 벽 온도를 결정하기 위하여 사용되고 각각의 층의 평균 온도를 끌어낸다.

Description

가상 온도 측정 포인트{VIRTUAL TEMPERATURE MEASURING POINT}
본 발명은 예를들어 스팀 수집기들, 스팀 라인들, 밸브 하우징들, 터빈 하우징들 또는 샤프트들 또는 등등 같은 벽들 또는 두꺼운 벽 구성요소들의 샤프트들 내 온도 프로파일 및 평균 적분 온도 및/또는 축 온도를 결정하기 위한 방법에 관한 것이다.
구성요소 벽들, 예를들어 스팀 터빈, 밸브 하우징 또는 스팀 라인에서 발생하는 바와 같은 가열 및 냉각 처리들 동안, 특히 동작 모드를 변경할 때, 온도 기울기들은 이들 구성요소들의 두꺼운 벽들에서 형성되고 상당한 재료 스트레스들을 유도할 수 있다. 이들 재료 스트레스들은 형태를 파괴하는 크기까지 너무 이른 재료 손상을 유도할 수 있다.
스팀 전력 플랜트들의 애플리케이션들에서 특히 상기 온도 기울기들을 모니터하기 위하여, 이전에 적어도 하나 또는 그 이상의 온도 측정 포인트들은 구성요소 벽에 통합되었다. 벽의 온도 및 작업 매체의 온도에 대해 결정된 측정 값들은 구성요소 벽 내에서 온도차들을 평가하고 특히 할당된 평균 적분 벽 온도를 결정하기 위하여 사용될 수 있다. 평균 적분 온도와 허용 가능한 제한 값들의 비교는 허용 가능한 제한들 내에서 열적 재료 스트레스들을 유지할 수 있게 한다. 그러나, 이 방법은 비교적 값이 비싸고 에러가 발생하기 쉽다.
선택적으로, 평균 적분 벽 온도는 측정 포인트(예를들어 터빈 샤프트)가 제공될 수 없는 벽내 또는 구성요소들의 경우 통합된 측정 포인트들이 값비싼 비용 및 에러 발생 경향 없이 계산될 수 있다. 하나의 가능한 방법은 수학적 대체 모델, 특히 금속 로드의 열 도전성에 대한 베셀(Bessel) 방정식을 바탕으로 이 온도를 계산하는 것이다. 그러나, 예를들어 스팀 전력 플랜트들의 튜브들 같은 산업 플랜트들의 기계 및 제어를 바탕으로 이전에 구현된 시스템들은 상기 방식으로 얻어진 온도 값들의 신뢰성 있는 평가를 제한하는 작업 매체의 온도 변화 기간에 따른 변동들을 겪는 경향을 가진다.
그러므로 본 발명의 목적은 온도 프로파일의 특히 정확한 픽쳐(picture)를 형성하고, 동시에 관련 벽에서 온도 측정 포인트들의 사용 없이 특히 견고하고 본질적으로 안정한 평균 적분 벽 온도/축 온도를 결정하기 위한 방법을 제공하는 것이다.
이 목적은 가열 또는 냉각 처리들 동안 평균 적분 벽 온도를 결정하기 위하여 각각의 층의 평균 온도를 바탕으로 다층 모델을 사용함으로써 본 발명에 따라 달성된다.
다층 모델이 상기 방식으로 사용될 때, 구성요소 벽은 표면에 평행하게 놓인 다수의 층들로 부가적으로 분할되고, 상기 층들의 수는 벽 두께에 종속한다. 각각의 층에 사용된 재료 데이터(열 용량, 열 도전성)는 층 기하구조에 무관하다. 진입하고 나가는 열 흐름의 일시적인 밸런싱은 각각의 층에서 발생한다. 얻어진 일시적인 열 밸런스는 대응 평균 층 온도들을 결정하기 위하여 사용된다.
다층 모델은 바람직하게 측정 값들로서 스팀 온도(
Figure 112008035304085-PCT00001
) 및 스팀 질량(
Figure 112008035304085-PCT00002
)뿐 아니라 밸런스된 처음 상태에서 처음 벽 온도(
Figure 112008035304085-PCT00003
)에 의해 표현될 수 있는 벽의 처음 온도 프로파일만을 사용한다. 만약 스팀 질량 흐름 측정이 없다면, 스팀 처리량은 압력(PAM) 및 밸브 위치(HAV) 또는 자유 흐름 단면을 바탕으로 대체 모델에 의해 계산된다. 이들 처리 변수들은 쉽게 얻어지고 임의의 경우 기술적 플랜트의 기계 및 제어에 일반적으로 이용할 수 있다. 특히, 관련된 벽에서 통합되어야 하는 부가적인 측정 포인트들은 요구되지 않는다.
본 발명은 벽의 온도 프로파일을 계산하고, 결과적으로 벽에 통합되는 비싸고 에러가 잘 발생하는 측정 포인트들 없이 및 직접적인 온도 측정이 가능하지 않은 경우들에서 가열 및 냉각 처리들 동안 다층 모델에 의해 충분히 정확하고 안정하게 평균 적분 벽 온도를 계산하는 것이 가능하다는 생각을 바탕으로 한다. 이런 목적을 위하여 일시적 열 흐름 밸런스의 함수로서 순간 벽 온도 프로파일의 결정은 예상된다. 본래, 두꺼운 벽 구성요소의 내부 및 외부 표면 온도들 또는 심지어 작업 매체의 온도 및 주위 또는 절연 온도 또는 심지어 표면 온도(예를들어 샤프트의 경우)를 사용하여 작업하는 것은 가능하다.
그러나, 두꺼운 벽들을 다수의 층들로 분할하는 것은 특히 바람직한 것으로 증명되었다. 결과적인 장점은, 두꺼운 벽내의 일시적인 온도 프로파일이 강한 비선형성을 가지기 때문에, 벽 온도 프로파일의 우수한 결정, 및 결과적으로 평균 적분 벽 온도의 보다 우수한 계산이 가능하다. 그 이유는 특히, 재료의 열 도전성 및 재료의 특정 열 용량 자체가 온도 종속적이기 때문이다. 다층 모델의 다른 장점은 만약 벽이 충분히 미세하게 다수의 층들로 분할되면, 순방향 계산 구조가 온도 종속 열 도전성 및 특정 열 용량을 위하여 사용될 수 있고, 즉 현재 층 대신 이전 층의 평균 온도가 사용될 수 있고, 이에 따라 양의 부호를 가질 수 있는 피드백을 피하고, 이에 따른 계산 회로는 보다 강건한 작용을 가진다.
열 전달 계수(α)의 계산은 바람직하게 스팀 응축, 습식 스팀 및 과열 스팀에 대한 허용 오차와 함께 이루어진다. 이로 인해, 작업 매체의 상태 검출은 하나의 모듈에서 이루어진다. 대응 스팀 엘리먼트 및 물 엘리먼트 모두의 가능한 응축, 및 과열 스팀 상태는 검출된다. 만약 과열된 스팀만 작업 매체로서 제공되면, 작업 매체로부터 벽의 제 1 층으로 열 흐름의 전달을 위한 열 전달 계수(αAM)는 바람직하게 스팀 처리량(
Figure 112008035304085-PCT00004
)의 함수로서 형성된다.
다른 한편, 만약 스팀 응축이 발생하면, 전달 계수(α)는 바람직하게 작업 매체의 응축 엘리먼트, 소위 응축 성분을 위하여 사용된 일정한 열 전달 계수(αW)에 의해 계산된다. 이런 응축 성분을 결정하기 위하여, 포화 온도(TS)는 압력(PAM), 작업 매체의 온도(TAM) 및 가열/냉각 표면(T1)의 온도(제 1 층의 평균 온도)의 함수로서 사용된다.
두꺼운 벽 구성요소(T1)의 제 1 층의 온도는 두 개의 값들 중 보다 큰 것으로부터 감산되고, 그 결과는 설정될 수 있는 상수(K)와 비교된다. 이들 두 개의 값들 중 보다 큰 것은, 피제수로서 작업 매체의 온도 및 포화 온도(TAM-TS) 사이의 차 및 포화 온도 및 두꺼운 벽 구성요소(TS-T1)의 제 1 층의 온도 사이의 차를 가지는 두 개의 몫 중 제수를 형성한다. 제 1 제수는 만약 양이면 과열된 스팀의 열 전달 계수(αAM)에 의해 곱셈되고, 제 2 제수는 만약 양이면 응축될 수 있도록 물에 대한 열 전달 계수(αW)에 의해 곱셈된다. 두 개의 곱셈들의 합은 과열된 스팀의 열 전달 계수(αAM)와 비교된다. 두 개의 값들 중 큰 것은 결과적인 열 전달 계수(α)이다.
평균 적분 벽 온도(
Figure 112008035304085-PCT00005
)의 계산은 n개의 개별 층들에서 진입하고 떠나는 열 흐름들의 일시적 밸런싱으로부터 특히 바람직하게 얻어진다. 이것은 소위 n개의 층 모듈에서 발생한다.
제 1 층 모듈에서, 제 1 층(
Figure 112008035304085-PCT00006
)으로 작업 매체의 열적 흐름 및 제 1 층에서 제 2 층(
Figure 112008035304085-PCT00007
)으로 열 흐름은 열 전달 계수(αl), 작업 매체의 온도(TAM) 및 평균 온도의 도움으로 계산된다. 관련된 층의 처음 온도(TAnf)에서, 제 1 층의 평균 온도(T1)는 제 1 층으로 작업 매체의 열 흐름 및 제 1 층에서 제 2 층으로 작업 매체의 열 흐름(
Figure 112008035304085-PCT00008
) 사이의 일시적 차이의 시간에 따른 적분에 의해 얻어진다.
k번째 층 모듈에서, k번째 층(Tk)의 평균 온도는 (k-1)번째 층(
Figure 112008035304085-PCT00009
)의 일시적 열 흐름 밸런스 및 k번째 층에서 (k+1)번째 층(
Figure 112008035304085-PCT00010
)으로 일시적 열 흐름의 도움으로 계산된다. k번째 층의 처음 온도(TAnf _k)에서, k 번째 층의 평균 온도(Tk)는 k번째 층으로 및 k번째 층으로부터 열 흐름들(
Figure 112008035304085-PCT00011
) 사이의 일시적 차이의 시간에 따른 적분에 의해 얻어진다.
마지막으로 최종 층 모듈에서, 최종(n번째) 층의 평균 온도(Tn)는 최종에서 하나 전 (n-1) 번째 층에서 최종(n번째) 으로 및 최종 층에서 열 절연부(
Figure 112008035304085-PCT00012
)로 일시적 열 전달 밸런스로부터 계산된다.
k번째 층에서 열 도전성(λk)의 온도 종속성 및 특정 열 용량(ck)은 바람직하게 지수적으로 2차의 다항식으로 근사화되거나, 대응 함수들에 의해 특정된다.
최종적으로, 평균 적분 벽 온도(
Figure 112008035304085-PCT00013
)는 층 재료의 질량 및 벽의 등가 부분 재료의 질량에 대한 허용오차와 함께 개별 층들의 평균 온도들(Tk)의 질량에 의해 특히 바람직한 방식으로 하나의 모듈에서 결정된다.
전체 방법은 특히 강화된 데이터 처리 시스템, 바람직하게 스팀 전력 플랜트의 기계 및 제어 시스템에서 수행된다.
본 발명으로 달성되는 장점들은 특히 두꺼운 벽 구성요소들의 벽 온도 프로파일들 및 평균 적분 벽 온도가 작업 매체의 질량 흐름 온도의 처리 파라미터들뿐 아니라 벽의 처음 온도 분배의 처리 파라미터들로부터 단독으로 신뢰성 있고 안정하게 지정될 수 있고, 만약 온도 처리량의 직접적인 측정이 이용될 수 없거나 가능하지 않으면, 부가적으로 압력 및 밸브 위치 또는 자유 흐름 단면으로 신뢰성 있고 안정하게 지정될 수 있고, 이에 따라 구성요소 벽들에 통합되는 측정 포인트들이 요구되지 않는다. 여기서 선택된 층들의 수가 많을수록, 평균 적분 벽 온도/축 온도의 크기가 보다 정확하게 된다.
3개의 층의 사용 및 절연 허용(제 4 층)을 가진 본 발명의 예시적인 실시예는 도면을 기초로 보다 상세히 설명된다.
도 1은 3개의 층들로 분할된 두꺼운 벽의 예로서 스팀 튜브를 통한 단면을 도시한다.
도 2는 열 전달 계수의 계산을 위한 모듈의 블록도를 도시한다.
도 3은 제 1 층의 평균 온도의 계산을 위한 모듈의 블록도를 도시한다.
도 4는 제 2 층의 평균 온도의 계산을 위한 모듈의 블록도를 도시한다.
도 5는 제 3 층의 평균 온도의 계산을 위한 모듈의 블록도를 도시한다.
도 6은 평균 적분 벽 온도의 계산을 위한 모듈의 블록도를 도시한다.
동일한 부분들은 모든 도면들에서 동일한 부호들이 제공된다.
도 1은 두꺼운 벽의 예로서 튜브 부분(1)을 단면으로 도시한다. 스팀 튜브 의 내부 공간(2)은 작업 매체(스트림)가 흐르고, 여기에서 열 전달은 제 1 층(4)으로 전달된다. 이런 흐름은 제 2 층(6) 및 제 3 층(8)으로 이어진다. 튜브 부분(1)은 절연부(10)에 의해 밀봉된다.
도 2에 따라, 스트림 처리량(
Figure 112008035304085-PCT00014
)의 측정된 값은 함수 발생기(32)에 입력 신호로서 공급되고, 상기 함수 발생기는 상기 입력 신호로부터 스팀의 경우 작업 매체(
Figure 112008035304085-PCT00015
)의 함수, 즉
Figure 112008035304085-PCT00016
)으로서 열 전달 계수(αAM)를 계산한다. 이 함수는 다수의 보간 포인트들에 의해 제공되고, 중간 값들은 적당한 보간 방법들에 의해 형성된다.
부분 응축의 경우에서, 작업 매체(PAM)의 압력은 함수 발생기(34)의 입력으로 통과되고, 상기 함수 발생기는 포화 함수(TS=f(PAM))를 복사하고, 결과적으로 각각의 압력에 대한 포화 온도(TS)를 그 출력에 공급한다. 이 함수는 보간 포인트들(스팀 테이블들의 압력 및 온도들)에 의해 제공되고, 중간 값들은 적당한 보간 방법들에 의해 결정된다.
작업 매체(TAM)의 온도는 최대치 발생기(36)에 의해 포화 온도(TS)와 비교된다. 제 1 층(T1)의 평균 온도는 감산기(38)에 의해 보다 큰 값으로부터 감산된다. 그 차이는 최대치 발생기(40)에 의해 설정될 수 있는 상수(K)와 비교된다. 결과적으로, 신호
Figure 112008035304085-PCT00017
는 최대치 발생기(40)의 출력에 제공된다. 상기 신호는 두 개의 분할기들(42 및 44)의 제수 입력들로 통과된다.
분할기(42)는 피제수 입력에서 감산기(46)에 의해 형성된 차(TAM-TS)를 수신한다. 함수 발생기(48)는 만약 양이면 곱셈기(50)의 하나의 입력에 신호
Figure 112008035304085-PCT00018
만을 통과시킨다. 신호는 증발된 작업 매체의 퍼센트, 소위 스팀 성분을 가리킨다. 만약 TAM-TS가 음이면, 즉 작업 매체의 온도가 포화 온도보다 낮으면, 신호 "영"은 곱셈기(50)의 대응 입력에 제공된다.
곱셈기(50)의 다른 입력에서, 스팀에 대한 열 전달 계수(αAM)는 제공된다. 그러므로, 스팀 성분이 웨이팅된 열 전달 계수(αPD)는 가산기(58)의 하나의 입력에 전달된다.
분할기(44)는 감산기(52)에 의해 형성된 차(TS-T1)를 그 피제수 입력에서 수신한다. 함수 발생기(54)는 만약 양이면 곱셈기(56)의 하나의 입력에 신호
Figure 112008035304085-PCT00019
만을 통과시킨다. 신호는 응결 성분에 의해 형성된 퍼센트를 나타낸다. 만약 차(TS-T1)가 음이면, 즉 제 1 층의 평균 온도가 포화 온도보다 높으면, 신호 "영"은 곱셈기(56)의 대응 입력에 제공된다.
물에 대한 열 전달 계수(αW)는 곱셈기(56)의 다른 입력에 제공된다. 그러므로, 응결 성분이 웨이팅된 열 전달 계수(αpW)는 가산기(59)의 제 2 입력으로 통과된다.
최대치 발생기(60)에서, 스팀의 경우에 대한 열 전달 계수는 하나의 입력에 제공된다. 가산기(58)에 의해 형성된 부분 응결의 경우에 대한 열 전달 계수(
Figure 112008035304085-PCT00020
)는 제 2 입력에 제공된다. 두 개의 값들 중 큰 것은 현재 열 전달 계수(α)이다.
만약 스트림 질량 흐름 측정치가 없으면, 스팀 질량 흐름은 예를들어 다음 계산 회로의 도움으로 계산된다. 함수 발생기(12)에서, 밸브 위치(HAV)의 실제 값은 자유 흐름 영역(AAV)으로 전환된다. 자유 흐름 영역은 곱셈기들(14 및 16)에 의해 적당한 전환 인자들(KU1 및 KU2)이 제공되고 추가 곱셈기(18)로 통과된다. 작업 매체(PAM)의 압력은 - 적당한 전환 상수(KU4)를 가진 곱셈기(20)에 의해 - 곱셈기(18)의 제 2 입력으로 통과되고, 그 결과는 곱셈기(22)의 입력으로 통과된다. 곱셈기(24)에 의해 적당한 전환 인자(KU5)가 제공된 작업 매체의 온도(TAM)는 분할기(26)의 분모 입력으로 통과되고, 분자 입력은 일을 수신한다. 역수 값은 출력에 제공된다. 역수 값의 루트는 루트 추출기(28)에 의해 곱셈기(22)의 제 2 입력으로 통과된다. 곱셈기(30)에 의해 적당한 전환 인자(KU3)가 제공된 곱셈기(22)의 출력 신호는 스팀 처리량(
Figure 112008035304085-PCT00021
)을 나타낸다. 대략, 다음이 스팀 처리량의 계산 동안 얻어진다:
Figure 112008035304085-PCT00022
.
도 3에 따른 제 1 층에 대한 모듈은 일시적 열 전달 밸런스로부터 제 1 층(T1)의 평균 온도를 결정한다. 이런 목적을 위하여, 온도 차(TAM-T1)는 우선 감산기(62)에 의해 우선 작업 매체(TAM)의 온도 및 제 1 층(T1)의 평균 온도로부터 형성되고 곱셈기(64)에 의해 열 전달 계수(α)에 의해 곱셈된다. 곱셈기(66)는 작업 매체로부터 구성요소 벽으로 열 전달을 위하여 설정할 수 있는 적당한 계수(KAL)를 가지며 등가의 제 1 표면을 나타내는 신호를 제공한다. 곱셈기(66)의 출력에는 감산기(68)의 피감수 입력으로 통과되는 작업 매체로부터 제 1 층으로의 열 전달 신호
Figure 112008035304085-PCT00023
이다.
예시적인 실시예에서, 제 1 층의 열 도전성(λ1) 및 특정 열 용량(C1)의 온도 종속성은 이차 다항식으로 근사화되고, 상기 다항식은 계수들 WO1, W11 및 W21뿐 아니라 C01, C11 및 C21에 의해 표현된다. 이 실시예에서 사용된 다항식들을 다음 형태를 가진다:
Figure 112008035304085-PCT00024
.
이것은 3개의 곱셈기들(70,72 및 74)의 입력들로 통과되는 작업 매체(TAM)의 온도에 의해 회로 측면에서 복사된다. 가능한 양의 피드백(재료의 특성들에 따라) 및 결과적으로 시스템 안정성의 증가를 위하여, 순방향 구조는 사용된다, 즉 작업 매체(TAM)의 온도는 제 1 층(T1)의 평균 온도 대신 사용된다.
열 도전성을 계산하기 위하여, 다항식 상수(W11)는 곱셈기(70)의 제 2 입력에 제공된다. 출력은 가산기(76)의 입력에 접속된다.
스쿼러(squarer)로서 접속된 곱셈기(72)의 출력에는 작업 매체(
Figure 112008035304085-PCT00025
)의 온도의 제곱 신호이다. 이것은 곱셈기(78)에 의해 다항식 상수(W21)에 의해 곱셈되고 추후 가산기(76)의 제 2 입력으로 통과된다.
다항식 상수(W01)는 가산기(76)의 제 3 입력으로 스위칭된다. 출력에는 상기 표현에 의해 제공된 온도 종속 열 도전성(λ1)이 제공된다.
특정 열 용량을 계산하기 위하여, 다항식 상수(C11)는 곱셈기(74)의 제 2 입력에 제공된다. 곱셈기(74)의 출력은 가산기(80)의 입력에 놓인다. 가산기(80)의 제 2 입력에는 다항식 상수(C01)가 제공된다. 곱셈기(72)의 출력에 제공된 작업 매체(
Figure 112008035304085-PCT00026
)의 온도 제곱은 곱셈기(82)에 의해 다항식 계수(C21)가 제공되고 추후 가산기(80)의 제 3 입력으로 통과된다. 그 출력에는 상기 표현에 의해 제공된 온도 종속 특정 열 용량(C1)이 제공된다.
감산기(84)는 제 1 층의 평균 온도들과 다음 층(T1-T2)으로부터의 온도 차를 형성한다. 이것은 곱셈기(86)에 의해 가산기(76)의 출력으로부터의 온도 종속 열 도전성(λ1)에 의해 곱셈되고 곱셈기(88)에 의해 층 두께 및 등가 표면에 따른 상수(KW1)에 의해 곱셈된다. 곱셈기(88)의 출력에는 제 1 층으로부터 제 2 층으로 열 흐름에 대한 신호인
Figure 112008035304085-PCT00027
가 제공된다.
이 신호는 감산기(68)의 감수 입력으로 통과된다. 그 출력에는 층 재료의 질량에 따른 제 1 층의 온도 변화율을 허용하는 계수(KT1)가 곱셈기(90)에 의해 제공되는 열 흐름 차(
Figure 112008035304085-PCT00028
)에 대한 신호가 제공된다.
결과적인 신호는 분할기(92)에 의해 온도 종속 특정 열 용량(C1)에 대해 가산기(80)의 출력에 제공되는 신호에 의해 분할된다.
내부 층의 평균 온도는 시간(t)에 따른 열 흐름 차의 적분에 의해 얻어진다.
Figure 112008035304085-PCT00029
분할기(92)의 출력에 제공된 신호는 적분기(94)에 제공되고, 상기 적분기는 처음 조건으로서 제 1 층(TAnf1)의 처음 온도를 사용한다.
도 4에 따른 제 2 층에 대한 모듈은 일시적 열 흐름 밸런스로부터 제 2 층(T2)의 평균 온도를 결정한다. 이런 목적을 위하여, 차(T2-T3)는 처음에 감산 기(96)에 의해 제 3 층(T3)의 온도 및 제 2 층(T2)의 평균 온도로부터 형성되고 곱셈기(98)에 의해 제 2 층의 온도 종속 열 도전성(λ2)에 의해 곱셈된다. 곱셈기(100)는 설정될 수 있고 층 두께 및 표면상 열 도전성에 의존하는 적당한 계수(KW2)를 가진 신호를 제공한다. 곱셈기(100)의 출력에는 감산기(102)의 감수 입력으로 통과되는 제 2 층으로부터 제 3 층으로 열 흐름에 대한 신호
Figure 112008035304085-PCT00030
가 제공된다.
감산기(102)의 피감수 입력에는 제 1 층으로부터 제 2 층으로 열 흐름(
Figure 112008035304085-PCT00031
)에 대한 신호가 제공된다. 입력은 열 흐름 차(
Figure 112008035304085-PCT00032
Figure 112008035304085-PCT00033
)를 공급한다. 곱셈기(104)는 층 재료의 질량에 따라 제 2 층의 온도 변화율을 허용하고 설정할 수 있는 계수(KT2)를 가진 신호를 제공한다. 추후, 상기 신호는 분할기(106)에 의해 제 2 층의 온도 종속 특정 열 용량(C2)에 의해 분할되고 그 다음 적분기(108)의 입력으로 통과된다. 적분기(108)는 처음 조건으로서 제 2 층의 처음 온도(TAnf2)를 사용한다. 그 출력에는 제 2 층의 평균 온도
Figure 112008035304085-PCT00034
가 제공된다.
열 도전성(λ2)의 온도 종속성 및 제 2 층의 특정 열 용량(C2)은 다시 계수들 W02, W12 및 W22뿐 아니라 C02, C12 및 C22를 가진 다항식들에 의해 근사화된다. 상 기 다항식들은 다음과 같다:
Figure 112008035304085-PCT00035
이것은 3개의 곱셈기들(110,112 및 114)의 입력들로 통과되는 제 1 층의 평균 온도(T1)에 의해 회로 측면에서 복사된다. 가능한 양의 피드백(재료의 특성들에 따라)을 피하고 및 결과적으로 시스템의 안정성 증가를 위하여, 순방향 구조는 사용되고, 즉 제 1 층(T1)의 평균 온도는 제 2 층(T2)의 평균 온도 대신 사용된다.
열 도전성을 계산하기 위하여, 다항식 상수(W12)는 곱셈기(110)의 제 2 입력에 제공된다. 출력은 가산기(116)의 입력에 접속된다.
스쿼러로서 접속된 곱셈기(112)의 출력에는 제 1 층(
Figure 112008035304085-PCT00036
)의 평균 온도의 제곱에 대한 신호이다. 이것은 곱셈기(118)에 의해 다항식 상수(W22)에 의해 곱셈되고 추후 가산기(116)의 제 2 입력으로 통과된다. 다항식 상수(W02)는 가산기(116)의 제 3 입력으로 스위칭된다. 출력에는 상기 표현에 의해 제공된 온도 종속 열 도전성(λ2)이 제공된다.
온도 종속 특정 열 용량을 계산하기 위하여, 다항식 상수(C12)는 곱셈기(114)의 제 2 입력에 제공된다. 곱셈기(114)의 출력은 가산기(120)의 입력에 놓인다. 가산기(120)의 제 2 입력에는 다항식 상수(C02)가 제공된다. 곱셈기(112)의 출력에 제공된 제 1 층의 평균 온도의 제곱(
Figure 112008035304085-PCT00037
)은 곱셈기(112)에 의해 다항식 계수(C22)가 제공되고 추후 가산기(120)의 제 3 입력으로 통과된다. 그 출력에는 상기 표현에 의해 제공된 온도 종속 특정 열 용량(C2)이 제공된다.
도 5에 따른 제 2 층에 대한 모듈은 열 흐름 밸런스로부터 제 3 층의 평균 온도(T3)를 결정한다. 이런 목적을 위하여, 온도 차(T3-TISOL)는 우선 감산기(124)에 의해 절연부의 온도(TISOL) 및 제 3 층(T3)의 온도로부터 형성되고 곱셈기(126)에 의해 설정될 수 있고 절연부의 열 손실 크기를 나타내는 적당한 상수(KISOL)에 의해 곱셈된다. 곱셈기(126)의 출력에는 감산기(128)의 감수로 통과되는 제 3 층으로부터 절연부로 열 흐름 신호가 제공된다(여기서 직접적으로 지정된 절연부의 열 손실 가능성이 있다).
Figure 112008035304085-PCT00038
감산기(128)의 피감수에는 제 2 층으로부터 제 3 층으로 열 흐름에 대한 신호(
Figure 112008035304085-PCT00039
)가 제공된다. 입력은 열 흐름 차(
Figure 112008035304085-PCT00040
)를 공급한다. 곱셈기(130)는 설정할 수 있고 층의 재료 질량에 따른 제 3 층 온도의 변화율을 허용하는 계수(KT3)를 가진 신호를 제공한다. 추후, 신호는 분할기(132)에 의해 제 3 층의 온도 종속 특정 열 용량(C3)에 의해 분할되고 적분기(134)의 입력으로 통과된다. 적분기(134)는 처음 조건으로서 제 3 층(TAnf3)의 처음 온도를 사용한다. 출력에는 제 3 층의 평균 적분 온도가 제공된다.
Figure 112008035304085-PCT00041
제 3 층의 특정 열 용량(C3)의 온도 종속성은 계수들(C03, C13 및 C23)을 가진 다항식에 의해 근사화된다.
다항식은 다음과 같다:
Figure 112008035304085-PCT00042
가능한 양의 피드백을 피하고(재료의 특성들에 따라) 및 결과적으로 시스템의 안정성을 증가시키기 위하여, 순방향 구조는 여기서 사용된다, 즉 제 2 층(T2)의 평균 온도는 제 3 층(T3)의 평균 온도 대신 사용된다.
이것은 두 개의 곱셈기들(136 및 138)의 입력들로 통과되는 제 2 층의 평균 온도(T2)에 의해 회로 측면에서 복사된다. 계수(C13)는 곱셈기(136)의 제 2 입력에 제공된다. 곱셈기(136)의 출력은 가산기(140)의 입력에 놓인다. 가산기(140)의 제 2 입력에는 다항식 상수(C03)이 제공된다. 곱셈기(138)의 출력에 제공된 제 2 층의 평균 온도의 제곱(
Figure 112008035304085-PCT00043
)은 곱셈기(142)에 의해 다항식 계수(C23)가 제공되고 추후 가산기(140)의 제 3 입력으로 통과된다. 출력에는 상기 표현에 의해 제공된 온도 종속 특정 열 용량(C3)이 제공된다.
도 6에 따라, 평균 적분 벽 온도(
Figure 112008035304085-PCT00044
)는 개별 층들의 평균 온도(T1,T2 및 T3)로부터 결정된다. 3개의 곱셈기들(144,146 및 148)은 층의 재료 질량에 대응하는 방식으로 개별 층들의 평균 온도들을 웨이팅하는 적당한 웨이팅 인자들(KG1,KG2 및 KG3)을 가진 온도 신호들을 제공한다. 웨이팅된 온도 신호들을 가산기(150)로 통과한다. 출력 신호는 곱셈기(152)에 의해 벽의 등가 부분 재료의 전체 웨이트의 영향을 허용하는 계수(KG)가 제공된다. 곱셈기(152)의 출력에는 평균 적분 벽 온도(
Figure 112008035304085-PCT00045
)에 대한 신호가 제공된다.

Claims (9)

  1. 두꺼운 벽 또는 샤프트에서 온도 프로파일 및 적분 평균 온도 및/또는 축 온도를 결정하기 위한 방법으로서,
    평균 적분 벽 온도는 가열 또는 냉각 처리들 동안 평균 적분 벽 온도를 결정하기 위하여 다층 모델의 각각의 층의 평균 온도로부터 계산되는,
    온도 결정 방법.
  2. 제 1 항에 있어서, 작업 매체의 온도 및 질량 흐름의 변수뿐 아니라 두꺼운 벽의 처음 온도는 계산에 사용되는,
    온도 결정 방법.
  3. 제 1 항 또는 제 2 항에 있어서, 스팀 처리량은 작업 매체의 압력 및 온도뿐 아니라 밸브 단면으로부터 결정되는,
    온도 결정 방법.
  4. 제 1 항 내지 제 3 항 중 어느 한 항에 있어서, 스팀만 과열된 경우 구성요소 벽에 작업 매체의 열 전달을 위한 열 전달 계수는 스팀 처리량의 함수로서 결정되는,
    온도 결정 방법.
  5. 제 1 항 내지 제 3 항 중 어느 한 항에 있어서, 응결의 존재시 두꺼운 벽 구성요소에 작업 매체의 열 전달을 위한 열 전달 계수는 응결 성분에 따라 결정되는,
    온도 결정 방법.
  6. 제 1 항 내지 제 5 항 중 어느 한 항에 있어서, 각각의 층의 평균 온도는 이 층에 대한 일시적 열 흐름 밸런스로부터 계산되는,
    온도 결정 방법.
  7. 제 6 항에 있어서, 상기 각각의 층의 열 도전성의 온도 종속성은 다항식, 특히 2차 다항식으로 근사화되는,
    온도 결정 방법.
  8. 제 6 항에 있어서, 상기 각각의 층의 특정 온도 용량의 온도 종속성은 다항식, 특히 2차 다항식으로 근사화되는,
    온도 결정 방법.
  9. 제 1 항 내지 제 8 항 중 어느 한 항에 있어서, 상기 평균 적분 벽 온도는 개별 층들의 재료의 질량 및 벽의 등가 부분의 전체 질량에 대한 허용오차를 가진 층들의 평균 온도들로부터 결정되는,
    온도 결정 방법.
KR1020087011939A 2005-10-19 2006-09-25 적분 평균 온도를 결정하기 위한 방법 KR101334885B1 (ko)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP05022820A EP1777505A1 (de) 2005-10-19 2005-10-19 Virtuelle Temperaturmessstelle
EP05022820.4 2005-10-19
PCT/EP2006/066704 WO2007045546A1 (de) 2005-10-19 2006-09-25 Virtuelle temperaturmessstelle

Publications (2)

Publication Number Publication Date
KR20080072855A true KR20080072855A (ko) 2008-08-07
KR101334885B1 KR101334885B1 (ko) 2013-12-02

Family

ID=36354090

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020087011939A KR101334885B1 (ko) 2005-10-19 2006-09-25 적분 평균 온도를 결정하기 위한 방법

Country Status (7)

Country Link
US (1) US7909506B2 (ko)
EP (2) EP1777505A1 (ko)
JP (1) JP4694625B2 (ko)
KR (1) KR101334885B1 (ko)
CN (1) CN101292141B (ko)
IL (1) IL190449A (ko)
WO (1) WO2007045546A1 (ko)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101520759B1 (ko) * 2013-09-27 2015-05-19 한국전력공사 보일러 튜브의 온도분포를 측정하기 위한 시스템 및 방법
JP6482970B2 (ja) * 2015-06-30 2019-03-13 大和ハウス工業株式会社 熱伝導率算出装置、温度予測装置、コンピュータプログラム、熱伝導率算出方法及び温度予測方法
DE102015226076A1 (de) * 2015-12-18 2017-06-22 Robert Bosch Gmbh Verfahren und Vorrichtung zum Bestimmen einer Temperatur einer Komponente in einem mechatronischem System

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE2504679A1 (de) * 1974-02-21 1975-08-28 Koppers Co Inc Einrichtung zum bestimmen des temperaturprofils einer heizzugwand bei einer koksofenbatterie
JPS61184432A (ja) * 1985-02-12 1986-08-18 Babcock Hitachi Kk 熱応力監視方法
KR900005989B1 (ko) * 1985-02-27 1990-08-18 미쓰비시전기 주식회사 가열로의 로온도 설정 및 가열제어방법
JPH0722129B2 (ja) * 1985-10-28 1995-03-08 株式会社日立製作所 酸化・拡散装置内ウエハ温度制御システム
JPS63217621A (ja) * 1987-03-06 1988-09-09 Hitachi Ltd 酸化・拡散装置内ウエハ間温度履歴均一化方法
JPH02291950A (ja) * 1989-05-02 1990-12-03 Ishikawajima Harima Heavy Ind Co Ltd 熱伝導率測定方法
DE4026731A1 (de) * 1990-08-24 1992-02-27 Peter Prof Dipl Ing Wippenbeck Verf. u. vorr. zur beherrschung der isochoren prozessphase beim spritzgiessen thermoplastischer kunststoffe
US5890805A (en) * 1997-09-26 1999-04-06 Usx Corporation Method for monitoring the wear and extending the life of blast furnace refractory lining
US6331075B1 (en) 1998-05-01 2001-12-18 Administrator, National Aeronautics And Space Administration Device and method for measuring thermal conductivity of thin films
US6487866B1 (en) 2000-07-10 2002-12-03 The United States Of America As Represented By The National Aeronautics & Space Administration Multipurpose thermal insulation test apparatus
CN1688948B (zh) * 2002-07-19 2010-05-26 布鲁克斯器具有限公司 在质量流动控制器中用于压力补偿的方法和装置
US6824305B1 (en) * 2002-08-16 2004-11-30 The Texas A & M University System Local wall heat flux/temperature meter for convective flow and method of utilizing same

Also Published As

Publication number Publication date
US7909506B2 (en) 2011-03-22
KR101334885B1 (ko) 2013-12-02
JP2009512842A (ja) 2009-03-26
JP4694625B2 (ja) 2011-06-08
EP1938067A1 (de) 2008-07-02
EP1777505A1 (de) 2007-04-25
WO2007045546A1 (de) 2007-04-26
US20090132196A1 (en) 2009-05-21
IL190449A (en) 2013-04-30
CN101292141B (zh) 2013-04-17
IL190449A0 (en) 2008-11-03
CN101292141A (zh) 2008-10-22

Similar Documents

Publication Publication Date Title
Ma et al. Distributed fusion filter for networked stochastic uncertain systems with transmission delays and packet dropouts
RU2509991C2 (ru) Способ и система для корректировки сигнала измерения температуры
JP2017215864A (ja) 監視装置、対象装置の監視方法、およびプログラム
JPS5814718B2 (ja) 応答遅れ補償回路
KR101334885B1 (ko) 적분 평균 온도를 결정하기 위한 방법
Dev et al. Prediction-based super twisting sliding mode load frequency control for multi-area interconnected power systems with state and input time delays using disturbance observer
JP2017072125A (ja) 蒸気タービンロータ応力を管理するための制御システム及び使用方法
WO2018181120A1 (ja) 情報処理装置、情報処理方法およびプログラム
JP2010507159A (ja) 遅延を伴った制御対象を同定するための方法および装置、制御装置およびコンピュータプログラム製品
Iannelli An implicit Galerkin finite element Runge–Kutta algorithm for shock-structure investigations
CN113031435B (zh) 数字式双频扩张状态观测器及扰动观测方法
JP4664842B2 (ja) エネルギープラントの最適運用システムと方法、およびプログラム
KR101397566B1 (ko) 3차원 그린함수와 가중계수를 사용한 열응력 계산방법
Gao et al. High-gain observer-based parameter identification with application in a gas turbine engine
Mukherjee et al. Order reduction of linear time invariant continuous systems using mixed approach
CN115630875B (zh) 一种同步电网的转动惯量评估方法及系统
JP6937787B2 (ja) 湿球温度の演算装置、湿球温度の演算方法
SU928041A1 (ru) Устройство дл контрол прогрева ротора турбины
Nalepa Correction of Dynamic Error by the „Blind” Method. A Differential Algorithm Simulation Study
SU1157271A1 (ru) Устройство дл эксплуатационного контрол осевых зазоров в проточной части цилиндра паровой турбины
JPS5843602Y2 (ja) タ−ビンケ−シングの熱応力監視装置
EP2837777A1 (en) Method for calculating thermal stress in a thick-walled component of a thermal system
Sunder Raj Multi-Pressure Surface Condenser Performance Evaluation: A Case Study
JPH02123944A (ja) 巻線温度予測方法
Naveros et al. Dynamic Heat Transfer In Walls: Limitations Of Heat Flux Meters

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
LAPS Lapse due to unpaid annual fee