KR102416995B1 - Torque Coefficient Calculation Method - Google Patents

Torque Coefficient Calculation Method Download PDF

Info

Publication number
KR102416995B1
KR102416995B1 KR1020200155393A KR20200155393A KR102416995B1 KR 102416995 B1 KR102416995 B1 KR 102416995B1 KR 1020200155393 A KR1020200155393 A KR 1020200155393A KR 20200155393 A KR20200155393 A KR 20200155393A KR 102416995 B1 KR102416995 B1 KR 102416995B1
Authority
KR
South Korea
Prior art keywords
torque
vibration component
torque coefficient
generator
power system
Prior art date
Application number
KR1020200155393A
Other languages
Korean (ko)
Other versions
KR20220068521A (en
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 연세대학교 산학협력단
Priority to KR1020200155393A priority Critical patent/KR102416995B1/en
Publication of KR20220068521A publication Critical patent/KR20220068521A/en
Application granted granted Critical
Publication of KR102416995B1 publication Critical patent/KR102416995B1/en

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H1/00—Measuring characteristics of vibrations in solids by using direct conduction to the detector
    • G01H1/003—Measuring characteristics of vibrations in solids by using direct conduction to the detector of rotating machines
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L3/00—Measuring torque, work, mechanical power, or mechanical efficiency, in general
    • G01L3/02—Rotary-transmission dynamometers
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L5/00—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
    • G01L5/16—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes for measuring several components of force
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

본 발명은 진동성분 검출을 통해 신속하게 토크 계수를 산정할 수 있는 방법에 관한 것으로, 고속 푸리에 변환 없이도 특정 주파수 영역을 필터링할 수 있는 제어기를 통해 발전기의 각속도 및 토크를 변화시키는 진동 성분만을 도출하여 토크 계수를 신속하게 산정 가능한 방법을 제공하고, 데이터 저장 개수 및 출력 데이터 용량이 적게 소모되므로 컴퓨터 성능과 저장 용량에 영향을 적게 받으며, 전력계통상 발전기와 직렬보상장치, HVDC 등과 같은 특수설비 사이에 발생되는 상호 뒤틀림 현상의 발생 여부를 특수설비가 전력계통에 설치되기 전 미리 검증할 수 있어 전력계통 시스템의 안정성을 향상시킬 수 있다.The present invention relates to a method for quickly calculating a torque coefficient by detecting vibration components, and by deriving only vibration components that change the angular speed and torque of a generator through a controller that can filter a specific frequency region without a fast Fourier transform It provides a method that can quickly calculate the torque coefficient, and consumes less data storage and output data capacity, so it is less affected by computer performance and storage capacity. It is possible to verify the occurrence of mutual distortion before special equipment is installed in the power system, thereby improving the stability of the power system system.

Description

토크 계수 산정 방법{Torque Coefficient Calculation Method}Torque Coefficient Calculation Method

본 발명은 발전기 및 발전기 회전축의 진동 성분 검출을 통해 신속하게 토크 계수를 산정할 수 있는 방법 및 연산장치에 관한 것이다.The present invention relates to a method and an arithmetic device capable of quickly calculating a torque coefficient through detection of a vibration component of a generator and a generator rotating shaft.

기존에는 복소 토크 계수를 산정하기 위한 방법은 EMT 모의실험 환경에서 터빈 발전기와 특수설비가 포함된 계통을 모델링하고, 상기 터빈 발전기에 특정 주파수에 해당하는 임의의 진동을 발생시켜 발전기의 시간별 각속도, 전기적 토크, 주파수 등의 데이터를 저장한 후, EMT 모의실험으로 얻은 데이터(각속도, 전기적 토크, 주파수 등)를 고속 푸리에 변환하는 방법이었다. 하지만 상당한 양의 데이터 저장과 저장된 데이터의 고속 푸리에 변환 처리에 필요한 툴을 통해 이산적으로 토크 계수를 도출해야 하므로 컴퓨터 성능과 추가 연산 툴이 필요하다는 문제가 있다. 따라서 본 발명은 상기와 같은 문제점을 해결하기 위해 안출한 것으로서, 복소 토크 계수 산정에 데이터 저장과 연산 단계를 줄이는 동시에 EMT 모의실험 상에서 결과를 도출하는 것에 그 목적이 있다.Conventionally, a method for calculating the complex torque coefficient is to model a system including a turbine generator and special equipment in an EMT simulation environment, and generate random vibrations corresponding to a specific frequency in the turbine generator to determine the angular velocity of the generator, electrical After storing data such as torque and frequency, it was a fast Fourier transform method of data (angular velocity, electrical torque, frequency, etc.) obtained from EMT simulation. However, there is a problem in that computer performance and additional calculation tools are required because torque coefficients must be discretely derived through a tool necessary for storing a significant amount of data and processing the stored data for a fast Fourier transform. Therefore, the present invention has been devised to solve the above problems, and it is an object of the present invention to reduce the data storage and calculation steps in calculating the complex torque coefficient and to derive the results on the EMT simulation.

공개특허공보 제10-2018-0086628호(2018.08.01.)Laid-open Patent Publication No. 10-2018-0086628 (2018.08.01.)

이에, 본 발명은 고속 푸리에 변환 없이도 특정 주파수 영역을 필터링할 수 있는 제어기를 통해 발전기 회전축의 각속도 및 발전기의 전기적 토크를 변화시키는 진동 성분만을 도출하여 토크 계수를 신속하게 산정 가능한 방법을 제공하는 것을 목적으로 한다.Accordingly, the present invention aims to provide a method capable of quickly calculating a torque coefficient by deriving only vibration components that change the angular velocity of the generator rotation shaft and the electrical torque of the generator through a controller capable of filtering a specific frequency region without fast Fourier transform do it with

상기 목적을 달성하기 위해, (S1) 각속도 및 토크를 전기적 신호로 출력하도록 모델링 된 전력계통상의 발전기 모델에서 상기 발전기에 포함된 회전축에 상기 각속도 및 발전기의 전기적 토크 중 어느 하나 이상을 변화시키는 제1 신호를 주입하는 단계; (S2) 상기 (S1) 단계에 의해 변화된 각속도 및 전기적 토크 중 어느 하나 이상을 출력하는 단계; (S3) 상기 (S2) 단계에서 상기 출력된 값을 제1 필터에 통과시켜 상기 제1 신호에 의한 진동 성분을 검출하는 단계; 및 (S4) 상기 (S3) 단계에서 검출된 상기 진동 성분을 이용하여 토크계수를 연산하는 단계;를 포함하고, 상기 연산된 토크계수로부터 상기 제1 신호에 대한 상기 전력계통의 안정도를 판별하는 토크계수 산정방법을 제공한다.In order to achieve the above object, (S1) a first to change any one or more of the angular velocity and the electrical torque of the generator to the rotating shaft included in the generator in the generator model of the power system modeled to output the angular velocity and torque as an electrical signal injecting a signal; (S2) outputting any one or more of the angular velocity and the electrical torque changed by the step (S1); (S3) passing the value output in step (S2) through a first filter to detect a vibration component by the first signal; and (S4) calculating a torque coefficient using the vibration component detected in step (S3). A method for calculating coefficients is provided.

또한, 모델링된 전력계통의 토크계수를 연산하는 장치로써, 상기 장치는 입력부, 연산부 및 출력부를 포함하고, 상기 연산부는 (S1) 각속도 및 토크를 전기적 신호로 출력하도록 모델링 된 전력계통상 발전기의 모델에서 상기 발전기에 포함된 회전축에 상기 각속도 및 발전기의 전기적 토크 중 어느 하나 이상을 변화시키는 제1 신호를 주입하는 단계; (S2) 상기 (S1) 단계에 의해 변화된 각속도 및 전기적 토크 중 어느 하나 이상을 출력하는 단계; (S3) 상기 (S2) 단계에서 상기 출력된 값을 제1 필터에 통과시켜 상기 제1 신호에 의한 진동 성분을 검출하는 단계; 및 (S4) 상기 (S3) 단계에서 검출된 상기 진동 성분을 이용하여 토크계수를 연산하는 단계;를 포함하여 수행하고, 상기 연산부는 상기 연산된 토크계수로부터 상기 전력계통의 안정도를 판단하는 토크계수 연산장치를 제공한다.In addition, as a device for calculating a torque coefficient of a modeled power system, the device includes an input unit, a calculation unit, and an output unit, and the calculation unit (S1) In the model of the generator on the power system modeled to output the angular velocity and torque as an electrical signal injecting a first signal for changing at least one of the angular velocity and the electrical torque of the generator to the rotation shaft included in the generator; (S2) outputting any one or more of the angular velocity and the electrical torque changed by the step (S1); (S3) passing the value output in step (S2) through a first filter to detect a vibration component by the first signal; and (S4) calculating a torque coefficient by using the vibration component detected in the step (S3), wherein the calculating unit determines the stability of the power system from the calculated torque coefficient. It provides an arithmetic unit.

본 발명에 따른 토크계수 산정방법은 발전기의 출력값 중에서 각속도 및 전기적 토크를 진동 및 발산하게 하는 진동 성분을 도출하여 산정하므로 토크계수를 산정하기 위한 연산 단계가 간소화되어 신속하게 도출할 수 있다.The torque coefficient calculation method according to the present invention is calculated by deriving and calculating the vibration component that causes the angular velocity and electric torque to vibrate and diverge from the output value of the generator, so the calculation step for calculating the torque coefficient is simplified and can be quickly derived.

또한, 기존의 고속 푸리에 변환 방법에 비해 데이터 저장 개수 및 출력 데이터 용량이 적게 소모되므로 컴퓨터 성능과 저장 용량에 영향을 적게 받는다.In addition, since the number of data storage and output data capacity are less consumed compared to the conventional fast Fourier transform method, the computer performance and storage capacity are less affected.

또한, 전력계통상 발전기와 직렬보상장치, HVDC 등과 같은 특수설비 사이에 발생되는 상호 뒤틀림 현상의 발생 여부를 특수설비가 전력계통에 설치되기 전 미리 검증할 수 있어 전력계통 시스템의 안정성을 향상시킬 수 있다.In addition, it is possible to verify in advance whether or not the occurrence of mutual distortion between the generator and the special equipment such as the series compensator and HVDC in the electric power system before the special equipment is installed in the electric power system, the stability of the electric power system system can be improved. .

도 1은 본 발명의 토크계수 산정을 위해 전력계통을 모델링한 것을 나타낸 것이다.
도 2은 본 발명의 토크계수 산정 방법을 위한 순서도를 나타낸 것이다.
도 3은 제1 신호가 주입되는 터빈과 발전기 사이의 회전축을 나타낸 것이다.
도 4는 제1 신호가 주입된 발전기와 터빈사이의 회전축의 각속도 변화를 나타내내 그래프이다.
도 5는 제1 신호가 발전기와 터빈 사이의 회전축에 주입된 후 발전기의 전기적 토크변화를 나타낸 그래프이다.
도 6은 발전기의 출력으로부터 진동 성분을 검출하기 위한 제1 필터를 나타낸 것이다.
도 7은 검출된 진동 성분을 저주파 통과 필터를 통해 위상을 변화시키는 방법을 도시한 것이다..
도 8은 검출된 진동 성분의 위상을 변화시켜 복소평면상에 나타낸 것이다.
도 9은 제1 신호의 주파수 변화에 따라 토크계수의 값을 나타낸 것이다.
도 10는 고속 푸리에 변환 방법론과 본 발명의 토크계수 산정방법의 데이터 저장 개수 및 출력 데이터 용량을 비교한 표이다.
도 11은 토크계수 연산장치에 관한 구성을 나타낸 것이다.1 shows the modeling of the power system for calculating the torque coefficient of the present invention.
Figure 2 shows a flowchart for the torque coefficient calculation method of the present invention.
Figure 3 shows the axis of rotation between the turbine and the generator to which the first signal is injected.
4 is a graph showing the change in the angular velocity of the rotation shaft between the generator and the turbine to which the first signal is injected.
5 is a graph showing a change in electrical torque of the generator after the first signal is injected into the rotation shaft between the generator and the turbine.
6 shows a first filter for detecting a vibration component from the output of the generator.
7 shows a method of changing a phase of a detected vibration component through a low-pass filter.
8 is a graph showing the phase of the detected vibration component is changed on a complex plane.
9 shows the values of the torque coefficient according to the change in the frequency of the first signal.
10 is a table comparing the data storage number and output data capacity of the fast Fourier transform methodology and the torque coefficient calculation method of the present invention.
11 shows the configuration of the torque coefficient calculating device.

이하 설명하는 발명은 다양한 변경을 가할 수 있고 여러 가지 실시예를 가질 수 있는 바, 특정 실시예들을 도면에 예시하고 상세하게 설명하고자 한다. 그러나, 이는 이하 설명하는 발명을 특정한 실시 형태에 대해 한정하려는 것이 아니며, 이하 설명하는 기술의 사상 및 기술 범위에 포함되는 모든 변경, 균등물 내지 대체물을 포함하는 것으로 이해되어야 한다.The invention to be described below may have various changes and may have various embodiments, and specific embodiments will be illustrated in the drawings and described in detail. However, this is not intended to limit the invention described below to specific embodiments, and it should be understood to include all modifications, equivalents, and substitutes included in the spirit and scope of the technology described below.

도 1 내지 도 10에는 본 발명의 복소 토크계수 산정방법에 관하여 설명하기 위한 자료들이 나타나 있다. 이하, 본 발명의 이해를 돕기 위해 첨부된 상기 도면들을 참고하여 본 발명을 보다 상세히 설명한다. 그러나 하기의 실시예는 본 발명을 보다 쉽게 이해하기 위해 제공되는 것일 뿐, 하기 실시예에 의해 본 발명의 내용이 한정되는 것은 아니다.1 to 10 show data for explaining the method for calculating the complex torque coefficient of the present invention. Hereinafter, the present invention will be described in more detail with reference to the accompanying drawings to help the understanding of the present invention. However, the following examples are only provided for easier understanding of the present invention, and the contents of the present invention are not limited by the following examples.

도 1은 본 발명에 따른 토크 계수 산정 방법에서 발전기가 포함된 전력계통을 모델링한 것을 나타낸 것이다. 도 1을 참조하여 설명하면, 차동기 공진 현상 중 하나인 차동기의 상호 비틀림 현상은 전력계통에 직렬보상장치, HVDC 등과 같은 특수설비와 터빈 발전기 사이에 발생하는 문제를 의미한다. 여기서 차동기 공진 현상은 전력계통과 발전기 간의 상호작용으로 특정 운전 조건에서 터빈과 발전기 사이의 축의 진동이 증폭되는 현상을 말한다. 상기 특정 운전 조건은 계통의 외란에 따라 발전기가 접속된 전원 계통이 상용주파수보다 작은 값의 공진주파수 성분을 함께 포함하는 조건을 의미한다. 특히, 발전기가 직렬커패시터에 의해 보상된 장거리 송전선로에 접속되는 경우 차동기 공진 현상에 취약하며, 이때 터빈과 발전기 사이의 축에 유기되는 진동 토크의 증폭으로 상기 축의 피로 및 손상을 가져올 수 있다.1 shows a modeling of a power system including a generator in the torque coefficient calculation method according to the present invention. Referring to FIG. 1 , the mutual torsion of the differential, which is one of the resonance phenomena of the differential, means a problem occurring between the turbine generator and special equipment such as a series compensator and HVDC in the power system. Here, the differential resonance phenomenon refers to a phenomenon in which the vibration of the shaft between the turbine and the generator is amplified under specific operating conditions due to the interaction between the power system and the generator. The specific operating condition refers to a condition in which the power supply system to which the generator is connected according to system disturbance includes a resonance frequency component of a value smaller than the commercial frequency. In particular, when the generator is connected to a long-distance power transmission line compensated by a series capacitor, it is vulnerable to a differential resonance phenomenon, and at this time, the amplification of the vibration torque induced in the shaft between the turbine and the generator may cause fatigue and damage of the shaft.

상기 문제점의 발생 여부에 관해 직렬커패시터, FACTS, HVDC 등과 같은 특수설비가 전력계통에 설치되지 전 계획단계에서 발전기와 상기 특수설비 사이에 상기 문제점이 발생할 수 있는지를 쉽게 검증할 수 있도록 전력계통 시스템을 모델링 하여 발전기와 직렬보상설비 연계를 통한 계통 구성과 계통 운전 상황별 모의실험 환경을 구성하였다. 상기 모의실험 환경은 전력계통 해석프로그램인 PSCAD/EMTDC를 이용하여 모델링을 하였으며, 이를 통해 시뮬레이션을 수행하였다.With regard to the occurrence of the problem, the power system system is designed to easily verify whether the problem can occur between the generator and the special equipment in the planning stage before special equipment such as series capacitors, FACTS, HVDC, etc. are installed in the electric power system. By modeling, the system configuration by linking the generator and the series compensation facility and the simulation environment for each system operation situation were constructed. The simulation environment was modeled using PSCAD/EMTDC, a power system analysis program, and simulation was performed through this.

도 2는 본 발명에 따른 토크계수 산정방법을 순서도로 나타낸 것이다. 도 2를 참조하여 설명하면, 본 발명에 따른 토크계수 산정방법은 (S1) 각속도 및 토크를 전기적 신호로 출력하도록 모델링 된 전력계통상의 발전기 모델에서 상기 발전기에 포함된 회전축에 상기 각속도 및 발전기의 전기적 토크 중 어느 하나 이상을 변화시키는 제1 신호를 주입하는 단계;, (S2) 상기 (S1) 단계에 의해 변화된 각속도 및 전기적 토크 중 어느 하나 이상을 출력하는 단계;, (S3) 상기 (S2) 단계에서 상기 출력된 값을 제1 필터에 통과시켜 상기 제1 신호에 의한 진동 성분을 검출하는 단계; 및 (S4) 상기 (S3) 단계에서 검출된 상기 진동 성분을 이용하여 토크계수를 연산하는 단계;를 포함한다.2 is a flowchart illustrating a method for calculating a torque coefficient according to the present invention. Referring to FIG. 2, the torque coefficient calculation method according to the present invention (S1) is the angular velocity and the electrical power of the generator on the rotating shaft included in the generator in the generator model of the power system modeled to output the angular velocity and torque as electrical signals. injecting a first signal for changing any one or more of the torques; (S2) outputting any one or more of the angular velocity and the electrical torque changed by the (S1) step; (S3) the (S2) step detecting a vibration component by the first signal by passing the output value through a first filter; and (S4) calculating a torque coefficient using the vibration component detected in step (S3).

제1 신호가 주입되면 발전기로부터 출력되는 전기적 토크 편차(

Figure 112020124266603-pat00001
)는 하기의 수학식(1)으로 표현될 수 있다.When the first signal is injected, the electric torque deviation (
Figure 112020124266603-pat00001
) can be expressed by the following Equation (1).

Figure 112020124266603-pat00002
수학식(1)
Figure 112020124266603-pat00002
Equation (1)

도 3은 제1 신호가 주입되는 위치를 나타낸 것이다. 도 3을 참조하면, 상기 제1 신호는 발전기의 회전축 또는 터빈과 발전기 사이의 회전축에 주입되어 상기 회전축의 각속도를 변화시 수 있다.3 shows a position where the first signal is injected. Referring to FIG. 3 , the first signal may be injected into the rotation shaft of the generator or the rotation shaft between the turbine and the generator to change the angular velocity of the rotation shaft.

도 4는 제1 신호가 터빈과 발전기 사이의 회전축에 주입된 후 상기 회전축의 시간에 대한 각속도 그래프를 나타낸 것이다. 도 4를 참조하여 설명하면, 상기 제1 신호는 각속도 및 토크를 전기적 신호로 출력하는 발전기와 터빈 사이의 회전축에 주입할 수 있다. 상기 발전기와 터빈 사이의 회전축은 제1 신호가 주입되기 전에는 발전기에서 출력되는 전기적 토크가 일정하게 유지될 수 있으나 상기 제1 신호가 상기 발전기와 터빈 사이의 회전축에 주입되면 상기 회전축에 진동이 발생하게 되어 상기 발전기의 전기적 토크에 변화가 발생하게 될 수 있다. 상기 제1 신호는 발전기의 전기적 토크에 변화를 주어 교란하는 섭동 신호로써, 각속도의 진폭이 시간에 따라 점점 증가하는 발산형의 형태로 상기 회전축에 주입될 수 있다.4 is a graph showing the angular velocity versus time of the rotation shaft after the first signal is injected into the rotation shaft between the turbine and the generator. Referring to FIG. 4 , the first signal may be injected into a rotating shaft between a generator and a turbine that outputs angular velocity and torque as electrical signals. The rotation shaft between the generator and the turbine may maintain a constant electrical torque output from the generator before the first signal is injected, but when the first signal is injected into the rotation shaft between the generator and the turbine, vibration occurs in the rotation shaft This may cause a change in the electrical torque of the generator. The first signal is a perturbation signal that is disturbed by changing the electrical torque of the generator, and may be injected into the rotation shaft in a divergent type in which the amplitude of the angular velocity gradually increases with time.

상기 제1 신호는 각도에 대한 미소 신호(

Figure 112020124266603-pat00003
)를 만든 후 미분기를 통과시켜 각속도에 대한 미소 신호(
Figure 112020124266603-pat00004
)를 만들어 상기 미소 신호가 포함된 새로운 각속도(
Figure 112020124266603-pat00005
)로 주입된다. 상기 제1 신호는 발전기와 터빈 사이의 회전축에 주입되어 발전기의 전기적 토크 편차 및 각속도 편차가 발생시킬 수 있다.The first signal is a minute signal for an angle (
Figure 112020124266603-pat00003
) and pass it through a differentiator to obtain a microsignal (
Figure 112020124266603-pat00004
) to create a new angular velocity (
Figure 112020124266603-pat00005
) is injected. The first signal may be injected into the rotation shaft between the generator and the turbine to generate an electrical torque deviation and an angular velocity deviation of the generator.

도 5는 터빈에서 출력되는 전기적 토크의 시간에 대한 전기적 토크를 그래프로 나타낸 것이다. 도 5를 참조하여 설명하면, 상기 제1 신호가 상기 발전기와 터빈 사이의 회전축에 주입되면 발전기에서는 출력되는 토크를 신호화한 전기적 토크에 변화가 발생 된다. 상기 발전기의 전기적 토크는 신호가 주입되기 전에는 일정하게 유지되고 있으나 상기 제1 신호가 주입된 후 전기적 토크의 진폭이 점점 증가하는 발산형의 형태로 그래프가 나타날 수 있다.5 is a graph showing the electrical torque versus time of the electrical torque output from the turbine. Referring to FIG. 5 , when the first signal is injected into the rotation shaft between the generator and the turbine, a change occurs in the electrical torque signaled by the torque output from the generator. Although the electrical torque of the generator is maintained constant before the signal is injected, the graph may appear in the form of a divergent type in which the amplitude of the electrical torque gradually increases after the first signal is injected.

상기 발전기에서 출력되는 전기적 토크에는 상기 제1 신호에 의한 축의 진동으로 발생 되는 진동 성분이 포함되어 있으므로 상기 진동 성분만을 검출할 필요가 있다. 상기 회전축의 각속도 변화로 상기 발전기의 전기적 토크(

Figure 112020124266603-pat00006
)의 변화가 발생하게 되며, 상기 전기적 토크는 토크 평균(
Figure 112020124266603-pat00007
)과 진동 성분(
Figure 112020124266603-pat00008
)으로 나타낼 수 있다. 전기적 토크 편차(
Figure 112020124266603-pat00009
)는 전기적 토크(
Figure 112020124266603-pat00010
)와 토크 평균(
Figure 112020124266603-pat00011
)의 차에 해당된다. 상기 전기적 토크 편차(
Figure 112020124266603-pat00012
)는 하기의 수학식과 같이 표현될 수 있다. Since the electrical torque output from the generator includes a vibration component generated by the vibration of the shaft by the first signal, it is necessary to detect only the vibration component. Electrical torque of the generator (
Figure 112020124266603-pat00006
) occurs, and the electrical torque is the torque average (
Figure 112020124266603-pat00007
) and the vibration component (
Figure 112020124266603-pat00008
) can be expressed as electrical torque deviation (
Figure 112020124266603-pat00009
) is the electrical torque (
Figure 112020124266603-pat00010
) and the torque average (
Figure 112020124266603-pat00011
) corresponds to the car of The electrical torque deviation (
Figure 112020124266603-pat00012
) can be expressed as the following equation.

Figure 112020124266603-pat00013
수학식(2)
Figure 112020124266603-pat00013
Equation (2)

Figure 112020124266603-pat00014
수학식(3)
Figure 112020124266603-pat00014
Equation (3)

상기 전기적 토크 편차(

Figure 112020124266603-pat00015
)는 복소평면에 하기의 수학식과 같이
Figure 112020124266603-pat00016
와
Figure 112020124266603-pat00017
로 표현될 수 있다.The electrical torque deviation (
Figure 112020124266603-pat00015
) is in the complex plane as in the following equation
Figure 112020124266603-pat00016
Wow
Figure 112020124266603-pat00017
can be expressed as

Figure 112020124266603-pat00018
수학식(4)
Figure 112020124266603-pat00018
Equation (4)

Figure 112020124266603-pat00019
수학식(5)
Figure 112020124266603-pat00019
Equation (5)

도 6은 발전기에서 출력되는 전기적 토크와 각속도를 제1 필터에 통과시켜 진동 성분을 검출하는 과정을 나타낸 것이다. 도 6을 참조하면, 상기 발전기에서 출력되는 각속도를 적분하여 각도 형태로 만들어 제1 필터에 통과시켜 각도 편차를 구할 수 있다.6 illustrates a process of detecting a vibration component by passing an electrical torque and an angular velocity output from the generator through a first filter. Referring to FIG. 6 , the angular velocity output from the generator may be integrated to form an angle, and the angular deviation may be obtained by passing it through the first filter.

또한, 상기 발전기에서 출력되는 전기적 토크에서 일정한 범위 이상의 주파수를 갖는 진동 성분을 검출하기 위해 발전기에서 출력되는 각속도 및 전기적 토크를 고역 통과 필터(high pass filter)와 대역 통과 필터(band pass filter)를 직렬로 연결한 제1 필터로 통과시킬 수 있다.In addition, in order to detect a vibration component having a frequency of a certain range or more in the electrical torque output from the generator, the angular velocity and the electrical torque output from the generator are combined with a high pass filter and a band pass filter in series. It can pass through the first filter connected to

상기 고역 통과 필터는 1차 또는 2차 필터일 수 있고, 밴드 패스 필터도 1차 또는 2차 필터일 수 있다.The high-pass filter may be a first-order or second-order filter, and the band-pass filter may also be a first-order or second-order filter.

상기 고역 통과 필터와 대역 통과 필터는 서로 직렬로 연결함으로써 진동 성분 검출의 정확성이 더욱 높아질 수 있다. 고역 통과 필터와 대역 통과 필터 중 어느 하나만 사용하게 되면 검출되는 성분에 진동 성분(

Figure 112020124266603-pat00020
) 뿐만 아니라 토크 평균(
Figure 112020124266603-pat00021
) 성분도 일부 포함될 수 있다.By connecting the high-pass filter and the band-pass filter in series with each other, the accuracy of detecting a vibration component may be further increased. If only one of the high-pass filter and the band-pass filter is used, the detected component is added to the vibration component (
Figure 112020124266603-pat00020
) as well as the torque average (
Figure 112020124266603-pat00021
) components may also be included.

본 발명은 모델링된 전력 계통상에서 발전기와 직렬보상장치, HVDC 등과 같은 특수설비 사이에 발생하는 차동기 공진 현상 중 하나인 차동기 상호 비틀림 현상의 발생여부를 판별하기 위한 것이다. 차동기 공진 현상은 50Hz 혹은 60Hz와 같은 규정주파수 이하의 단일 혹은 복수 주파수에서 임의의 한 발전기와 이를 제외한 나머지 전력계통이 에너지를 주고받는 상태를 의미하므로, 60Hz 이하의 주파수 영역에서 모델링된 전력계통의 발전기와 특수설비 사이에 차동기 상호 비틀림 현상이 발생하는지를 판별할 수 있다. 따라서 이러한 상기 비틀림 현상이 발생하는지는 알고자 하는 60Hz 이하 규정주파수를 상기 고역 통과 필터 및 대역 통과 필터에 입력하여 토크계수를 산정하여 판별할 수 있다.An object of the present invention is to determine whether a differential mutual torsion phenomenon, which is one of differential resonance phenomena occurring between a generator, a series compensator, and special equipment such as HVDC, occurs in a modeled power system. Since the differential resonance phenomenon refers to a state in which an arbitrary generator and the rest of the power system exchange energy at single or multiple frequencies below a specified frequency such as 50Hz or 60Hz, the power system modeled in the frequency domain below 60Hz It can be determined whether a differential mutual torsion phenomenon occurs between the generator and special equipment. Therefore, it is possible to determine whether the twisting phenomenon occurs by inputting a specified frequency of 60 Hz or less to the high-pass filter and the band-pass filter to calculate the torque coefficient.

도 7은 검출된 진동 성분을 제2 필터 저역 통과 필터에 통과시켜 위상을90°변화시키는 방법을 나타낸 것이고, 도 8은 위상이 변화된 진동 성분과 위상이 변화전 진동 성분을 복소평면상에 나타내어 진동성분의 크기 및 위상을 구하는 방법을 나타낸 것이다. 도 7 및 도 8을 참조하여 설명하면, 상기 진동 성분은 상기 제1 신호에 의해 변화하는 상기 회전축의 각도 편차 및 상기 제1 신호에 의해 변화하는 상기 발전기의 전기적 토크 편차를 포함하고, 상기 검출된 진동 성분을 복소 평면상에서 크기 및 위상을 구하기 위해 발전기에서 출력되는 전기적 토크 편차 및 각도 편차를 제2 필터인 저역 통과 필터(low pass filter)에 통과시켜 위상을 90°변화 시킬 수 있다((S5) 단계). 상기 로우 패스 필터를 통과한 상기 전기적 토크 편차 및 각도 편차는 복소평면상 허수축에 나타낼 수 있고, 로우 패스 필터를 통과하기 전의 전기적 토크 편차 및 각도 편차는 복소평면상 실수축에 나타낼 수 있어 전기적 토크 편차 및 각도 편차의 위상 및 크기를 구할 수 있게 된다((S6) 단계). 7 shows a method of changing the phase by 90° by passing the detected vibration component through the second filter low-pass filter, and FIG. A method for obtaining the magnitude and phase of a component is shown. 7 and 8, the vibration component includes an angular deviation of the rotation shaft that is changed by the first signal and an electrical torque deviation of the generator that is changed by the first signal, and the detected In order to obtain the magnitude and phase of the vibration component on a complex plane, the electric torque deviation and angle deviation output from the generator are passed through a low pass filter, which is a second filter, to change the phase by 90° ((S5) step). The electric torque deviation and angular deviation passing through the low-pass filter may be represented on an imaginary axis on a complex plane, and the electrical torque deviation and angular deviation before passing through the low-pass filter may be indicated on a real axis on the complex plane, so that the electrical torque It is possible to obtain the phase and magnitude of the deviation and the angle deviation (step S6).

상기 저역 통과 필터에 입력하는 주파수도 차동기 상호 비틀림 현상이 발생하는지는 알고자 하는 60HZ 이하의 주파수를 입력한다. 상기 입력하는 주파수는 고역 통과 필터, 대역 통과 필터 및 저역 통과 필터 모두 동일한 값으로 입력한다.The frequency input to the low-pass filter is also input to a frequency of 60 Hz or less, which is to be determined whether the mutual distortion phenomenon of the differential occurs. The input frequency is input as the same value for the high-pass filter, the band-pass filter, and the low-pass filter.

Figure 112020124266603-pat00022
와
Figure 112020124266603-pat00023
는 전기적 토크 편차의 위상 및 크기를 통해 하기의 수학식으로 산정할 수 있다.
Figure 112020124266603-pat00022
Wow
Figure 112020124266603-pat00023
can be calculated by the following equation through the phase and magnitude of the electrical torque deviation.

Figure 112020124266603-pat00024
수학식(6)
Figure 112020124266603-pat00024
Equation (6)

Figure 112020124266603-pat00025
수학식(7)
Figure 112020124266603-pat00025
Equation (7)

Figure 112020124266603-pat00026
수학식(8)
Figure 112020124266603-pat00026
Equation (8)

산정한

Figure 112020124266603-pat00027
와
Figure 112020124266603-pat00028
를 통해 댐핑 토크 계수(
Figure 112020124266603-pat00029
) 및 동기 토크 계수(
Figure 112020124266603-pat00030
)를 연산할 수 있다. 상기 댐핑 토크 계수 및 동기 토크 계수는 하기의 수학식과 같이 표현될 수 있다.calculated
Figure 112020124266603-pat00027
Wow
Figure 112020124266603-pat00028
damping torque coefficient (
Figure 112020124266603-pat00029
) and synchronous torque coefficient (
Figure 112020124266603-pat00030
) can be calculated. The damping torque coefficient and the synchronous torque coefficient may be expressed as follows.

Figure 112020124266603-pat00031
수학식(9)
Figure 112020124266603-pat00031
Equation (9)

Figure 112020124266603-pat00032
수학식(10)
Figure 112020124266603-pat00032
Equation (10)

Figure 112020124266603-pat00033
수학식(11)
Figure 112020124266603-pat00033
Equation (11)

각도 편차 및 각속도 편차를 수학식(9)에 대입하면 하기의 수학식(12)과 같이 표현될 수 있다. Substituting the angular deviation and the angular velocity deviation into Equation (9), it can be expressed as Equation (12) below.

Figure 112020124266603-pat00034
수학식(12)
Figure 112020124266603-pat00034
Equation (12)

Figure 112020124266603-pat00035
수학식(13)
Figure 112020124266603-pat00035
Equation (13)

댐핑 토크 계수(

Figure 112020124266603-pat00036
) 및 동기 토크 계수(
Figure 112020124266603-pat00037
)는 하기의 수학식(15) 및 수학식(16)으로 산정할 수 있다((S7) 단계).damping torque coefficient (
Figure 112020124266603-pat00036
) and synchronous torque coefficient (
Figure 112020124266603-pat00037
) can be calculated by the following Equations (15) and (16) (step S7).

Figure 112020124266603-pat00038
수학식(14)
Figure 112020124266603-pat00038
Equation (14)

Figure 112020124266603-pat00039
수학식(15)
Figure 112020124266603-pat00039
Equation (15)

Figure 112020124266603-pat00040
수학식(16)
Figure 112020124266603-pat00040
Equation (16)

Figure 112020124266603-pat00041
수학식(17)
Figure 112020124266603-pat00041
Equation (17)

전력계통의 안정도란 주어진 초기 운전조건에서 외란발생시 평형 운전상태로 복귀하기 위한 전력계통의 능력을 의미한다. 이러한 전력계통의 안정도는 토크계수를 산정하여 그 값으로 판단할 수 있다. 상기 토크계수는 외란에 의해 평형운전상태를 벗어난 발전기의 전기적 토크의 진동을 평형운전상태 내의 전기적 토크 범위로 되돌릴 수 있는 능력이 있는지 판단할 수 있다. 상기 토크계수에는 댐핑 토크 계수(

Figure 112020124266603-pat00042
) 및 동기 토크 계수(
Figure 112020124266603-pat00043
)가 있고, 이 중 상기 댐핑 토크 계수(
Figure 112020124266603-pat00044
)는 평형상태를 벗어나 진동하는 전기적 토크가 발산하는 것을 감쇠하여 다시 평형상태로 할 수 있는 능력이 있는지 즉, 안정도가 확보된 시스템인지 아닌지를 판단할 수 있도록 하는 계수이다. The stability of the power system refers to the power system's ability to return to the balanced operating state in the event of disturbance in a given initial operating condition. The stability of such a power system can be determined by calculating the torque coefficient. The torque coefficient can determine whether there is an ability to return the vibration of the electric torque of the generator out of the balanced operation state due to the disturbance to the electric torque range within the balanced operation state. The torque coefficient includes the damping torque coefficient (
Figure 112020124266603-pat00042
) and synchronous torque coefficient (
Figure 112020124266603-pat00043
), of which the damping torque coefficient (
Figure 112020124266603-pat00044
) is a coefficient that allows it to be judged whether the system has the ability to return to the equilibrium state by damping the dissipation of the oscillating electrical torque out of the equilibrium state, that is, whether the system is stable or not.

도 9는 제1 신호의 연속적인 주파수 변화에 따른 댐핑 토크 계수 및 동기 토크 계수 값의 변화를 그래프로 나타낸 것이다. 연산된 댐핑 토크 계수(

Figure 112020124266603-pat00045
) 의 부호가 양수 또는 0인 경우, 발전기 회전축에 발생 되는 진동 성분을 감쇠시켜 초기의 평형상태로 되돌릴 수 있다는 것을 의미하고, 음수인 경우 상기 발전기 회전축에 발생 되는 진동 성분을 감쇠시킬 수 있는 능력이 없어 발산시킬 수 있는 구성이다. 도 9를 참조하면, 주파수의 변화에 따라 댐핑 토크 계수(
Figure 112020124266603-pat00046
) 및 동기 토크 계수(
Figure 112020124266603-pat00047
)의 값이 변화하는데 특정한 주파수의 범위에서 상기 댐핑 토크 계수(
Figure 112020124266603-pat00048
) 및 동기 토크 계수(
Figure 112020124266603-pat00049
) 부호가 음수에 해당되는 언더 슈트가 발생하였으므로, 상기 특정 주파수 범위에서 발전기와 특수설비 사이에 안정도는 불안정하여 문제가 발생할 수 있다는 것을 나타낸다. 상기 특정 주파수 범위 내의 진동 성분이 유입된다면 발전기의 전기적 토크는 상기의 진동 성분을 감쇠하여 다시 평형상태로 되돌릴 수 없고 전기적 토크가 발산하게 하므로 결국 상기 발전기와 특수설비 사이에 문제가 발생될 수 있다는 것을 예측할 수 있다. 9 is a graph showing changes in values of a damping torque coefficient and a synchronous torque coefficient according to a continuous frequency change of the first signal. Calculated damping torque coefficient (
Figure 112020124266603-pat00045
), if the sign is positive or 0, it means that the vibration component generated on the generator rotation shaft can be attenuated and returned to the initial equilibrium state, and when it is negative, the ability to attenuate the vibration component generated on the generator rotation shaft is It is a composition that can be diffused without it. Referring to FIG. 9 , the damping torque coefficient (
Figure 112020124266603-pat00046
) and synchronous torque coefficient (
Figure 112020124266603-pat00047
) changes in the damping torque coefficient (
Figure 112020124266603-pat00048
) and synchronous torque coefficient (
Figure 112020124266603-pat00049
) sign corresponds to a negative number, so the stability between the generator and the special equipment in the specific frequency range is unstable, indicating that a problem may occur. If the vibration component within the specific frequency range is introduced, the electrical torque of the generator attenuates the vibration component and cannot be returned to the equilibrium state again, and the electrical torque is emitted. predictable.

PSCAD/EMTDC 전력계통 해석프로그램을 활용한 모의실험에서는 이산적/연속적 결과 모두 얻어낼 수 있기에 컴퓨터 성능과 저장 용량에 크게 영향을 받지 않는다. 도 10은 고속 푸리에 변환을 사용하기 위해 얻은 데이터와 본 발명을 통해 최종적으로 도출한 결과의 데이터 용량을 비교한 표이다. 고속 푸리에 변환 방법론과 본 발명의 토크계수 산정방법은 데이터의 개수 및 용량에서 확연한 차이를 확인할 수 있다. 이러한 차이는 본 발명은 토크 계수 산정 과정에서 데이터 저장과 연산 단계를 줄이며, 모의실험 상에서 그 결과를 즉시 도출할 수 있어 고속 푸리에 변환 방법을 이용하여 토크계수를 산정하는 경우보다 더 쉽고 신속하게 결과를 도출할 수 있다.In simulation using PSCAD/EMTDC power system analysis program, both discrete and continuous results can be obtained, so it is not significantly affected by computer performance and storage capacity. 10 is a table comparing the data capacity of data obtained to use the fast Fourier transform and the result finally derived through the present invention. The fast Fourier transform methodology and the torque coefficient calculation method of the present invention can confirm a clear difference in the number and capacity of data. The present invention reduces the data storage and calculation steps in the torque coefficient calculation process, and the results can be obtained immediately in the simulation, so that the results are easier and faster than the case of calculating the torque coefficient using the fast Fourier transform method. can be derived

본 발명에 따른 토크계수 산정방법을 활용하면 발전기와 직렬보상장치, HVDC, 직렬캐패시터, TCSC(Thyristor Controlled Series Capacitor) 및 STATCOM(Static Synchronous Compensator) 등과 같은 특수설비 사이에 발생하는 뒤틀림 문제의 발생 가능 여부를 상기 특수설비가 전력계통에 설치되기 전에 토크계수를 산정하여 미리 판별할 수 있다.If the torque coefficient calculation method according to the present invention is used, the distortion problem occurring between the generator and the series compensator, HVDC, series capacitor, TCSC (Thyristor Controlled Series Capacitor) and STATCOM (Static Synchronous Compensator), etc. is possible can be determined in advance by calculating the torque coefficient before the special equipment is installed in the power system.

도 11 은 본 발명에 따른 토크계수 연산장치에 관한 것이다. 11 relates to a torque coefficient calculating device according to the present invention.

상기 연산장치는 입력부, 연산부 및 출력부를 포함할 수 있다. 상기 연산부는 상기 (S1), (S2), (S3) 및 (S4) 단계를 포함하여 수행하고, 상기 연산부는 (S1) 각속도 및 토크를 전기적 신호로 출력하도록 모델링 된 전력계통상 발전기의 모델에서 상기 발전기에 포함된 회전축에 상기 각속도 및 발전기의 전기적 토크 중 어느 하나 이상을 변화시키는 제1 신호를 주입하는 단계;, (S2) 상기 (S1) 단계에 의해 변화된 각속도 및 전기적 토크 중 어느 하나 이상을 출력하는 단계;, (S3) 상기 (S2) 단계에서 상기 출력된 값을 제1 필터에 통과시켜 상기 제1 신호에 의한 진동 성분을 검출하는 단계; 및 (S4) 상기 (S3) 단계에서 검출된 상기 진동 성분을 이용하여 토크계수를 연산하는 단계;를 포함하여 수행하고, 상기 연산된 토크계수로부터 상기 전력계통의 안정도를 판단을 수행할 수 있다.The arithmetic unit may include an input unit, an arithmetic unit, and an output unit. The calculation unit is performed including the steps (S1), (S2), (S3) and (S4), and the calculation unit (S1) is modeled to output angular velocity and torque as electrical signals. injecting a first signal for changing any one or more of the angular velocity and the electrical torque of the generator to the rotating shaft included in the generator; (S2) outputting any one or more of the angular velocity and the electrical torque changed by the (S1) step (S3) passing the value output in the step (S2) through a first filter to detect a vibration component by the first signal; and (S4) calculating a torque coefficient using the vibration component detected in the step (S3); and determining the stability of the power system from the calculated torque coefficient.

또한, 상기 연산부는 상기 (S4) 단계에 (S5) 검출된 진동 성분을 저역 통과 필터(low pass filter)에 통과시켜 위상을 변화시키는 단계;, (S6) 상기 (S3) 단계에서 검출된 진동 성분을 실수축으로 하고, 상기(S5) 단계에서 위상을 변화시킨 진동 성분을 허수축으로 하여 크기 및 위상을 구하는 단계;, (S7) 상기 (S6) 단계에서 상기 진동 성분의 크기 및 위상을 이용하여 토크계수를 산정하는 단계;를 포함하여 수행할 수 있다.In addition, the calculating unit (S4) step (S5) passing the detected vibration component through a low pass filter (low pass filter) to change the phase; (S6) the vibration component detected in the step (S3); is the real axis, and the imaginary axis is the vibration component whose phase is changed in step (S5) to obtain the magnitude and phase; (S7) using the magnitude and phase of the vibration component in step (S6) Calculating the torque coefficient; can be performed including.

상기 연산부는 상기 연산된 토크계수로 전력계통상 발전기와 특수설비 사이의 안정도를 판단할 수 있고, 구체적으로 산정된 댐핑토크계수의 부호가 음수이면 상기 안정도는 불안정한 시스템으로 판단하고, 부호가 0 또는 양수이면 안정한 시스템으로 판단할 수 있다. 상기 안정도는 발전기가 출력하는 전기적 토크에 포함된 진동성분을 감쇠시켜 평형상태로 되돌아갈 수 있는지를 판단할 수 있는 지표에 해당할 수 있다.The calculation unit may determine the stability between the generator and the special equipment in the power system with the calculated torque coefficient, and if the sign of the specifically calculated damping torque coefficient is negative, the stability is determined as an unstable system, and the sign is 0 or positive This can be judged as a stable system. The stability may correspond to an index capable of determining whether it is possible to return to the equilibrium state by attenuating the vibration component included in the electrical torque output by the generator.

상기 출력부는 상기 판단부에서 판단한 결과를 사용자가 알기 쉽게 디스플레이 화면상에 표시할 수 있다.The output unit may display a result determined by the determination unit on a display screen for a user to easily understand.

이상 실시예를 통해 본 기술을 설명하였으나, 본 기술은 이에 제한되는 것은 아니다. 상기 실시예는 본 기술의 취지 및 범위를 벗어나지 않고 수정되거나 변경될 수 있으며, 본 기술분야의 통상의 기술자는 이러한 수정과 변경도 본 기술에 속하는 것임을 알 수 있을 것이다.Although the present technology has been described through the above embodiments, the present technology is not limited thereto. The above embodiments may be modified or changed without departing from the spirit and scope of the present technology, and those skilled in the art will recognize that such modifications and changes also belong to the present technology.

Claims (12)

(S1) 각속도 및 토크를 전기적 신호로 출력하도록 모델링 된 전력계통상의 발전기 모델에서 상기 발전기에 포함된 회전축에 상기 각속도 및 발전기의 전기적 토크 중 어느 하나 이상을 변화시키는 제1 신호를 주입하는 단계;
(S2) 상기 (S1) 단계에 의해 변화된 각속도 및 전기적 토크 중 어느 하나 이상을 출력하는 단계;
(S3) 상기 (S2) 단계에서 출력된 값을 제1 필터에 통과시켜 상기 제1 신호에 의한 진동 성분을 검출하는 단계; 및
(S4) 상기 (S3) 단계에서 검출된 상기 진동 성분을 이용하여 토크계수를 연산하는 단계;를 포함하고,
상기 연산된 토크계수로부터 상기 제1 신호에 대한 상기 전력계통의 안정도를 판별하며,
상기 (S4) 단계는,
(S5) 검출된 진동 성분을 저역 통과 필터(low pass filter)에 통과시켜 위상을 변화시키는 단계;
(S6) 상기 (S3) 단계에서 검출된 진동 성분을 실수축으로 하고, 상기(S5) 단계에서 위상을 변화시킨 진동 성분을 허수축으로 하여 크기 및 위상을 구하는 단계; 및
(S7) 상기 (S6) 단계에서 상기 진동 성분의 크기 및 위상을 이용하여 토크계수를 산정하는 단계를 포함하는 토크계수 산정방법.(S1) injecting a first signal for changing any one or more of the angular velocity and the electrical torque of the generator to the rotating shaft included in the generator in the generator model of the power system modeled to output the angular velocity and torque as an electrical signal;
(S2) outputting any one or more of the angular velocity and the electrical torque changed by the step (S1);
(S3) passing the value output in step (S2) through a first filter to detect a vibration component by the first signal; and
(S4) calculating a torque coefficient using the vibration component detected in the step (S3);
Determining the stability of the power system with respect to the first signal from the calculated torque coefficient,
The step (S4) is,
(S5) changing the phase by passing the detected vibration component through a low pass filter;
(S6) taking the vibration component detected in step (S3) as a real axis and the vibration component having the phase changed in step (S5) as an imaginary axis to obtain magnitudes and phases; and
(S7) A torque coefficient calculation method comprising the step of calculating a torque coefficient using the magnitude and phase of the vibration component in the step (S6). 제 1 항에 있어서,
상기 전력계통의 안정도를 판별하는 방법으로써
상기 연산된 토크계수의 부호가 음수이면 상기 전력계통의 안정도는 불안정한 것으로 판단하고, 0 또는 양수이면 안정한 것으로 판단하는 토크계수 산정방법.The method of claim 1,
As a method of determining the stability of the power system
If the sign of the calculated torque coefficient is negative, it is determined that the stability of the power system is unstable, and if it is 0 or a positive number, it is determined that the stability of the power system is stable. 제 1 항에 있어서,
상기 진동 성분은 상기 제1 신호에 의해 변화하는 상기 회전축의 각도 편차로부터 발생하는 진동 성분 및 상기 제1 신호에 의해 변화하는 상기 발전기의 전기적 토크 편차로부터 발생하는 진동 성분 중 어느 하나 이상을 포함하는 토크계수 산정방법.The method of claim 1,
The vibration component is a torque including any one or more of a vibration component generated from an angular deviation of the rotation shaft that is changed by the first signal and a vibration component generated from an electrical torque deviation of the generator that is changed by the first signal Coefficient calculation method. 제 1 항에 있어서,
상기 제1 필터는 고역 통과 필터(high pass filter) 및 대역 통과 필터(band pass filter)인 토크계수 산정방법.The method of claim 1,
The first filter is a high pass filter (high pass filter) and a band pass filter (band pass filter) is a torque coefficient calculating method. 삭제delete 제 1 항에 있어서,
상기 (S7) 단계의 토크계수는 댐핑 토크 계수 및 동기 토크 계수를 포함하며, 상기 댐핑 토크 계수 및 상기 동기 토크 계수는 하기 수학식으로 연산되는 토크계수 산정방법.
Figure 112022024911094-pat00050

(
Figure 112022024911094-pat00051
: 전기적 토크 편차,
Figure 112022024911094-pat00052
: 댐핑 토크 계수,
Figure 112022024911094-pat00053
: 동기 토크 계수,
Figure 112022024911094-pat00054
: 각도 편차 크기
Figure 112022024911094-pat00055
: 각속도 편차 크기,
Figure 112022024911094-pat00056
: 초기 각속도,
Figure 112022024911094-pat00057
: 진동성분의 실수축 크기,
Figure 112022024911094-pat00058
: 진동성분의 허수축 크기,
Figure 112022024911094-pat00059
: 변화된 주파수)The method of claim 1,
The torque coefficient in step (S7) includes a damping torque coefficient and a synchronous torque coefficient, and the damping torque coefficient and the synchronous torque coefficient are calculated by the following equation.
Figure 112022024911094-pat00050

(
Figure 112022024911094-pat00051
: electrical torque deviation,
Figure 112022024911094-pat00052
: damping torque coefficient,
Figure 112022024911094-pat00053
: synchronous torque coefficient,
Figure 112022024911094-pat00054
: Angular deviation size
Figure 112022024911094-pat00055
: angular velocity deviation magnitude,
Figure 112022024911094-pat00056
: initial angular velocity,
Figure 112022024911094-pat00057
: Real axis size of vibration component,
Figure 112022024911094-pat00058
: The size of the imaginary contraction of the vibration component,
Figure 112022024911094-pat00059
: changed frequency) 모델링된 전력계통의 토크계수를 연산하는 장치로써,
상기 장치는 입력부, 연산부 및 출력부를 포함하고,
상기 연산부는
(S1) 각속도 및 토크를 전기적 신호로 출력하도록 모델링 된 전력계통상 발전기의 모델에서 상기 발전기에 포함된 회전축에 상기 각속도 및 발전기의 전기적 토크 중 어느 하나 이상을 변화시키는 제1 신호를 주입하는 단계;
(S2) 상기 (S1) 단계에 의해 변화된 각속도 및 전기적 토크 중 어느 하나 이상을 출력하는 단계;
(S3) 상기 (S2) 단계에서 출력된 값을 제1 필터에 통과시켜 상기 제1 신호에 의한 진동 성분을 검출하는 단계; 및
(S4) 상기 (S3) 단계에서 검출된 상기 진동 성분을 이용하여 토크계수를 연산하는 단계;를 포함하여 수행하고,
상기 연산부는 상기 연산된 토크계수로부터 상기 전력계통의 안정도를 판단하며,
상기 (S4) 단계는,
(S5) 검출된 진동 성분을 저주파 통과 필터(low pass filter)에 통과시켜 위상을 변화시키는 단계;
(S6) 상기 (S3) 단계에서 검출된 진동 성분을 실수축으로 하고, 상기(S5) 단계에서 위상을 변화시킨 진동 성분을 허수축으로 하여 크기 및 위상을 구하는 단계; 및
(S7) 상기 (S6) 단계에서 상기 진동 성분의 크기 및 위상을 이용하여 토크계수를 산정하는 단계를 포함하는 토크계수 연산장치.As a device for calculating the torque coefficient of the modeled power system,
The device includes an input unit, a calculation unit and an output unit,
the calculation unit
(S1) injecting a first signal for changing any one or more of the angular velocity and the electrical torque of the generator to the rotation shaft included in the generator in the model of the generator on the power system modeled to output the angular velocity and torque as an electrical signal;
(S2) outputting any one or more of the angular velocity and the electrical torque changed by the step (S1);
(S3) passing the value output in step (S2) through a first filter to detect a vibration component by the first signal; and
(S4) calculating a torque coefficient using the vibration component detected in step (S3);
The calculation unit determines the stability of the power system from the calculated torque coefficient,
The step (S4) is,
(S5) passing the detected vibration component through a low pass filter to change the phase;
(S6) taking the vibration component detected in step (S3) as a real axis and the vibration component having the phase changed in step (S5) as an imaginary axis to obtain magnitudes and phases; and
(S7) A torque coefficient calculating device comprising the step of calculating a torque coefficient by using the magnitude and phase of the vibration component in the step (S6). 제 7 항에 있어서,
상기 전력계통의 안정도를 판별하는 방법으로써
상기 연산된 토크계수의 부호가 음수이면 상기 전력계통의 안정도는 불안정한 것으로 판단하고, 0 또는 양수이면 안정한 것으로 판단하는 토크계수 연산장치.8. The method of claim 7,
As a method of determining the stability of the power system
If the sign of the calculated torque coefficient is negative, it is determined that the stability of the power system is unstable, and if it is 0 or a positive number, it is determined that the stability of the power system is stable. 제 7 항에 있어서,
상기 진동 성분은 상기 제1 신호에 의해 변화하는 상기 회전축의 각도 편차로부터 발생하는 진동 성분 및 상기 제1 신호에 의해 변화하는 상기 발전기의 전기적 토크 편차로부터 발생하는 진동 성분 중 어느 하나 이상을 포함하는 토크계수 연산장치.8. The method of claim 7,
The vibration component is a torque including any one or more of a vibration component generated from an angular deviation of the rotation shaft that is changed by the first signal and a vibration component generated from an electrical torque deviation of the generator that is changed by the first signal counting unit. 제 7 항에 있어서,
상기 제1 필터는 고역 통과 필터(high pass filter) 및 대역 통과 필터(band pass filter)인 토크계수 연산장치.8. The method of claim 7,
The first filter is a high pass filter (high pass filter) and a band pass filter (band pass filter) is a torque coefficient calculating device. 삭제delete 제 7 항에 있어서,
상기 (S7) 단계의 토크계수는 댐핑 토크 계수 및 동기 토크 계수를 포함하며, 상기 댐핑 토크 계수 및 상기 동기 토크 계수는 하기 수학식으로 연산되는 토크계수 연산장치.
Figure 112022024911094-pat00060

(
Figure 112022024911094-pat00061
: 전기적 토크 편차,
Figure 112022024911094-pat00062
: 댐핑 토크 계수,
Figure 112022024911094-pat00063
: 동기 토크 계수,
Figure 112022024911094-pat00064
: 각도 편차 크기
Figure 112022024911094-pat00065
: 각속도 편차 크기,
Figure 112022024911094-pat00066
: 초기 각속도,
Figure 112022024911094-pat00067
: 진동성분의 실수축 크기,
Figure 112022024911094-pat00068
: 진동성분의 허수축 크기,
Figure 112022024911094-pat00069
: 변화된 주파수)8. The method of claim 7,
The torque coefficient of step (S7) includes a damping torque coefficient and a synchronous torque coefficient, and the damping torque coefficient and the synchronous torque coefficient are calculated by the following equation.
Figure 112022024911094-pat00060

(
Figure 112022024911094-pat00061
: electrical torque deviation,
Figure 112022024911094-pat00062
: damping torque coefficient,
Figure 112022024911094-pat00063
: synchronous torque coefficient,
Figure 112022024911094-pat00064
: Angular deviation size
Figure 112022024911094-pat00065
: angular velocity deviation magnitude,
Figure 112022024911094-pat00066
: initial angular velocity,
Figure 112022024911094-pat00067
: Real axis size of vibration component,
Figure 112022024911094-pat00068
: The size of the imaginary contraction of the vibration component,
Figure 112022024911094-pat00069
: changed frequency)
KR1020200155393A 2020-11-19 2020-11-19 Torque Coefficient Calculation Method KR102416995B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020200155393A KR102416995B1 (en) 2020-11-19 2020-11-19 Torque Coefficient Calculation Method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020200155393A KR102416995B1 (en) 2020-11-19 2020-11-19 Torque Coefficient Calculation Method

Publications (2)

Publication Number Publication Date
KR20220068521A KR20220068521A (en) 2022-05-26
KR102416995B1 true KR102416995B1 (en) 2022-07-05

Family

ID=81808991

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020200155393A KR102416995B1 (en) 2020-11-19 2020-11-19 Torque Coefficient Calculation Method

Country Status (1)

Country Link
KR (1) KR102416995B1 (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005180924A (en) 2003-12-16 2005-07-07 Matsushita Electric Ind Co Ltd Gear damage detector and gear damage detection method
KR101432786B1 (en) * 2013-11-14 2014-09-23 엠앤디테크놀로지 주식회사 Fault diagnosis method of motor and system thereof
KR101901127B1 (en) 2017-09-29 2018-09-27 주식회사 만도 Apparatus and method for detecting a faulty motor

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0751999A (en) * 1993-08-06 1995-02-28 Fanuc Ltd Tool breakage detecting method
KR20180086628A (en) 2017-01-23 2018-08-01 현대엘리베이터주식회사 Remote diagnosis method for elevator using by real time frequency analysis of motor current and speed

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005180924A (en) 2003-12-16 2005-07-07 Matsushita Electric Ind Co Ltd Gear damage detector and gear damage detection method
KR101432786B1 (en) * 2013-11-14 2014-09-23 엠앤디테크놀로지 주식회사 Fault diagnosis method of motor and system thereof
KR101901127B1 (en) 2017-09-29 2018-09-27 주식회사 만도 Apparatus and method for detecting a faulty motor

Also Published As

Publication number Publication date
KR20220068521A (en) 2022-05-26

Similar Documents

Publication Publication Date Title
Kejian et al. A uniform control method for imbalance compensation and automation balancing in active magnetic bearing-rotor systems
US9006933B2 (en) Power system stabilization
US10892695B2 (en) Motor drive system
Xu et al. Characterization of subsynchronous oscillation with wind farms using describing function and generalized Nyquist criterion
Zilli et al. Nonlinear dynamics of a simplified model of an overhung rotor subjected to intermittent annular rubs
CN111512134B (en) Method for estimating internal effective torque of torque generator
Tchokogoué et al. The effects of gravity on the response of centrifugal pendulum vibration absorbers
Xie et al. Numerical tracking of limit points for direct parametric analysis in nonlinear rotordynamics
KR102416995B1 (en) Torque Coefficient Calculation Method
Riemann et al. Real gyroscopic uncertainties in robust control of flexible rotors
Hollenbach et al. A fluid–structure interaction tool using a van der pol-based reduced-order model for buffet and nonsynchronous vibrations
Alzarooni et al. Breathing crack model effect on rotor's postresonance backward whirl
CN103117692A (en) Control method of mechanical elastic energy storing permanent magnet motor group under various external disturbances
JP5358668B2 (en) Power system stable equilibrium point calculation device
Sinha Vibration absorbers for a mistuned bladed disk
Chen et al. Active vibration control based on linear matrix inequality for rotor system under seismic excitation
Villa et al. The invariant manifold approach applied to nonlinear dynamics of a rotor-bearing system
CN101515719B (en) Method for establishing generator shafting multimass block damping variable model and device thereof
CN104931129B (en) A kind of electric generating set rotor torsional vibration measurement system and measuring method
Miranda et al. Finite element method applied to the eigenvalue analysis of flexible rotors supported by journal bearings
Peng et al. Synchronous vibration control for a class of cross‐coupled antisymmetric MSR systems
Wramner Numerical Algorithms for Simulation of One-Dimensional Mechanical Systems With Clearance-Type Nonlinearities
CN114709873B (en) Method for judging stability of small interference of grid-connected voltage source type converter with phase-locked loop
CN118094784B (en) Determination method of joint response in double-degree-of-freedom accessory system under joint excitation
Behrens et al. Modelling of technical systems: application to hydrodynamic torque converters and couplings

Legal Events

Date Code Title Description
E701 Decision to grant or registration of patent right
GRNT Written decision to grant