TW201734870A - Operation support device and wind power system - Google Patents

Operation support device and wind power system Download PDF

Info

Publication number
TW201734870A
TW201734870A TW106104632A TW106104632A TW201734870A TW 201734870 A TW201734870 A TW 201734870A TW 106104632 A TW106104632 A TW 106104632A TW 106104632 A TW106104632 A TW 106104632A TW 201734870 A TW201734870 A TW 201734870A
Authority
TW
Taiwan
Prior art keywords
failure
data
risk
maintenance
probability
Prior art date
Application number
TW106104632A
Other languages
Chinese (zh)
Other versions
TWI607328B (en
Inventor
Norio Takeda
Hiroshi Shintani
Kazuo Muto
Tomoaki Yamashita
Original Assignee
Hitachi Ltd
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 Hitachi Ltd filed Critical Hitachi Ltd
Publication of TW201734870A publication Critical patent/TW201734870A/en
Application granted granted Critical
Publication of TWI607328B publication Critical patent/TWI607328B/en

Links

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03DWIND MOTORS
    • F03D17/00Monitoring or testing of wind motors, e.g. diagnostics
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03DWIND MOTORS
    • F03D80/00Details, components or accessories not provided for in groups F03D1/00 - F03D17/00
    • F03D80/50Maintenance or repair
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/72Wind turbines with rotation axis in wind direction

Abstract

The purpose of the present invention is to provide an operation support system which is capable of referring to high-precision evaluation values of reliability of a plurality of components which configure a product, and of carrying out a management and maintenance plan of the product. Provided is an operation support device 100 of an arbitrary product, comprising a fault risk evaluation unit 2 which: derives a fault risk of a plurality of components which configure the product, said fault risk being computed on the basis of information which includes at least one of environmental data, working data, design data, and/or materials data, of the product from the past to the present; and derives fault risk estimation values of the plurality of components which fluctuate if a current management and maintenance plan of the product is varied. The operation support device 100 further comprises a maintenance/management scenario creation unit 3 which is for referring to each of the computed fault risk estimation values of the plurality of components, and assigning the management of the product and maintenance times of the plurality of components.

Description

運作輔助裝置及風力發電系統 Operation aids and wind power systems

本發明,特別有關任意製品的運作輔助裝置及風力發電系統。 The invention relates in particular to an operation aid for an arbitrary product and a wind power generation system.

在廠房等的製品中組裝有多數的感測器,藉由感測器計測出的資訊,係被利用在基於穩定運作目的之控制、製品的構成要素之可靠性分析(故障徵兆診斷、剩餘壽命診斷等)的維護計劃之擬定等。例如,發電廠中係藉由多種多樣的感測器來實施主要3種類的計測。該3種類的計測,分別被稱為控制計測(SCADA,Supervisory Control And Data Aquisition)、狀態監視(CMS、Condition Monitoring System)及構造物監視(SHM、Structural Health Monitoring)。風力發電廠的情形下,控制計測(SCADA)中,目的在於掌握風車的環境條件或運作狀態來適當地控制風車,而會計測風速、風向、發電量、發電機的轉數、溫度等各式各樣的物理量。狀態監視(CMS)中,目的在於偵測風車故障的徵兆,將故障所造成的被害抑制在最小限度,而進行計測。此外構造物監視(SHM) 中,目的在於評估風車的葉片等之健全性,而計測構造物的形變等。一般的風車中,會實施這樣的控制計測、狀態監視及構造物監視的全部或部分,來做風車的控制同時評估可靠性,實現風車的穩定運轉。 A large number of sensors are assembled in products such as factories, and the information measured by the sensors is used for reliability analysis based on the control of stable operation and the components of the products (diagnosis of failure signs, remaining life) Diagnosis, etc.) The maintenance plan is prepared. For example, in a power plant, the main three types of measurements are implemented by a variety of sensors. The three types of measurement are referred to as SCADA (Supervisory Control And Data Aquisition), Condition Monitoring (CMS, Condition Monitoring System), and Structure Monitoring (SHM, Structural Health Monitoring). In the case of wind power plants, in the control measurement (SCADA), the purpose is to grasp the environmental conditions or operating conditions of the windmill to properly control the windmill, and to measure the wind speed, wind direction, power generation, number of revolutions of the generator, temperature, etc. Various physical quantities. In the state monitoring (CMS), the purpose is to detect the symptoms of the windmill failure, and to minimize the damage caused by the failure, and to perform the measurement. In addition to structure monitoring (SHM) In order to evaluate the soundness of the blade of the windmill, etc., the deformation of the structure and the like are measured. In a general windmill, all or part of such control measurement, state monitoring, and structure monitoring is performed to control the windmill and evaluate the reliability to achieve stable operation of the windmill.

專利文獻1中,揭示一種運轉控制程式及風車,係「一種風車的運轉控制裝置,具備:疲勞劣化排程,讓風車的累積運轉時間與前述風車的最佳疲勞劣化度建立對應;及疲勞劣化演算手段,演算前述風車的目前的疲勞劣化度;及運轉控制手段,因應藉由前述疲勞劣化演算手段而演算出的前述風速的疲勞劣化度與從前述疲勞劣化排程取得的目前的最佳疲勞劣化度之關係,來控制前述風車的運轉」(請求項1)。 Patent Document 1 discloses an operation control program and a windmill, which is a type of operation control device for a windmill, which includes: a fatigue deterioration schedule, and an accumulated operation time of the wind turbine is associated with an optimum fatigue deterioration degree of the wind turbine; and fatigue degradation The calculation means calculates the current fatigue deterioration degree of the windmill; and the operation control means, the fatigue deterioration degree of the wind speed calculated by the fatigue deterioration calculation means, and the current optimum fatigue obtained from the fatigue deterioration schedule The operation of the windmill is controlled by the relationship of the degree of deterioration (request 1).

專利文獻2中,揭示「一種風力電廠的運轉控制系統,係具備複數個風車之風力電廠的運轉控制系統,其特徵為,具備:剩餘壽命預測部,針對各風車預測零件的剩餘壽命;及售電收入預測部,針對各風車預測複數個輸出限制條件下的售電收入;及維護成本預測部,針對各風車依據前述零件的剩餘壽命預測各輸出限制條件下的維護成本;及輸出限制條件選擇部,針對各風車依據對每一前述輸出限制條件預測出的售電收入及維護成本,來針對各風車選擇會讓從風力電廠獲得的收益成為最大之輸出限制條件;及運轉司令部,依據選擇出的輸出限制條件將運轉指令送至各風車」(請求項5)。 Patent Document 2 discloses "an operation control system for a wind power plant, which is an operation control system for a wind power plant having a plurality of windmills, and is characterized in that: a remaining life prediction unit is provided, and the remaining life of the components is predicted for each windmill; The electric revenue prediction unit predicts the sales revenue under the plurality of output restriction conditions for each windmill; and the maintenance cost prediction unit predicts the maintenance cost under each output restriction condition for each windmill based on the remaining life of the component; and the output restriction condition selection For each windmill, based on the sales revenue and maintenance cost predicted for each of the aforementioned output constraints, the choice of each windmill will maximize the output gains from the wind power plant; and the operation command, based on the selection The output restriction condition is sent to each windmill" (request item 5).

專利文獻3中,揭示一種廠房機器的運用診斷裝置, 「具備:將構成廠房的機器及各個構件的目前運轉環境下之對於各個破壞現象的破壞機率,乘上對每一該構件事先決定好的每一損傷形態的權重係數而得之最大值,訂為廠房風險推定值之手段;及計算運用條件而求出廠房運轉限制值之手段,該運用條件係讓從將前述各構件的設想運轉條件下之對於各個破壞現象的破壞機率,乘上每一前述損傷形態的權重係數而得之數值中求出最大值而成之廠房風險值,不超出設定值;及從前述構件的目前時間點的剩餘壽命評估資訊,依據運用計劃計算今後的壽命消費之手段;及基於目前時間點的剩餘壽命資訊、及將依據運用計劃而計算出的將來的剩餘壽命預測為基準,評估每一運用計劃的破壞機率的演變,將前述各構件的每一損傷形態的權重係數乘上前述破壞機率的演變資料,算出廠房運用風險推定值之手段」(摘要)。 Patent Document 3 discloses an operation diagnostic device for a plant machine. "With: the maximum probability of damage to each damage phenomenon in the current operating environment of the machine and each component that constitutes the plant, multiplied by the weight coefficient of each damage form determined in advance for each component, Means for estimating the value of the plant risk; and means for calculating the operation limit value of the plant by calculating the operating conditions, which are used to multiply the probability of destruction of each failure phenomenon under the assumed operating conditions of the aforementioned components The plant risk value obtained by determining the maximum value of the damage form of the damage form does not exceed the set value; and the remaining life evaluation information from the current time point of the member is calculated, and the future life consumption is calculated according to the operation plan. Means; and based on the remaining life information at the current time point and the future remaining life prediction calculated based on the operational plan, the evolution of the damage probability of each application plan is evaluated, and each damage form of each of the aforementioned components is Multiply the weighting factor by the evolution data of the above-mentioned damage probability, and calculate the means by which the plant uses the risk estimation value. Yes).

〔先前技術文獻〕 [Previous Technical Literature] 〔專利文獻〕 [Patent Document]

〔專利文獻1〕日本特開2006-241981號公報 [Patent Document 1] Japanese Patent Laid-Open Publication No. 2006-241981

〔專利文獻2〕日本特開2013-170507號公報 [Patent Document 2] Japanese Patent Laid-Open Publication No. 2013-170507

〔專利文獻3〕日本特開2002-73155號公報 [Patent Document 3] Japanese Patent Laid-Open Publication No. 2002-73155

上述專利文獻1、2中,是以疲勞損傷率及剩餘壽命來作為可靠性的評估基準,但該些文獻中,並未揭示當藉 由複數個零件來獲得疲勞損傷率或剩餘壽命的情形下,該依據哪個零件的可靠性(疲勞損傷率或剩餘壽命)來實施廠房的運轉控制或維護計劃。另一方面,專利文獻3中,是對複數個構件的各者,計算將破壞機率乘上權重係數而得之數值,並將其最大值訂為廠房風險推定值來運用廠房,因此揭示了該著眼於哪一零件的可靠性來運用廠房,並計劃維護。但,廠房的維護計劃中,不僅是風險值成為最大值之零件,還必須考量應維護之其他零件的風險值來擬定維護計劃。但,上述任一專利文獻中,均未揭示或教示考量複數個零件的可靠性評估值(風險值等)來擬定運用、維護計劃,以使廠房高效率且穩定地運作之手段。又,專利文獻3,是設想蒸氣渦輪這類火力發電廠,可以想見並未將曝露於戶外嚴苛環境變化的風車、建設機械等中使用之物視為對象。 In the above Patent Documents 1 and 2, the fatigue damage rate and the remaining life are used as the evaluation criteria for reliability, but in these documents, it is not disclosed. In the case where the fatigue damage rate or the remaining life is obtained from a plurality of parts, the operation control or maintenance plan of the plant is implemented based on which part reliability (fatigue damage rate or remaining life). On the other hand, in Patent Document 3, for each of a plurality of members, a numerical value obtained by multiplying the probability of destruction by a weighting coefficient is calculated, and the maximum value is set as a factory risk estimation value to operate the plant, thereby revealing the Focus on the reliability of which parts to use the plant and plan maintenance. However, in the maintenance plan of the plant, not only the parts whose risk value becomes the maximum value, but also the risk value of other parts to be maintained must be considered to formulate the maintenance plan. However, none of the above-mentioned patent documents discloses or teaches a reliability evaluation value (risk value, etc.) of a plurality of parts to formulate an operation and maintenance plan to enable the plant to operate efficiently and stably. Further, Patent Document 3 assumes a thermal power plant such as a steam turbine, and it is conceivable that objects used in windmills, construction machines, and the like that are not exposed to harsh outdoor environments are not considered.

此外,由疲勞損傷率、剩餘壽命、破壞機率、風險值等所得的可靠性評估,有著如下課題,即,特別是當對象製品為新開發之製品的情形下,因缺乏故障資料等,可靠性的評估精度未必良好。上述任一專利文獻,亦揭示了從對象製品的資料來評估可靠性,並將其利用在維護或運用,但相對於此,並未揭示或教示例如使用和對象製品同型之製品、類似之製品的可靠性的評估值,來更新對象製品的可靠性評估值,並依據高精度化後的評估值來運用製品、計劃維護。 In addition, the reliability evaluation obtained from the fatigue damage rate, the remaining life, the probability of failure, the risk value, and the like has a problem that, in particular, when the target product is a newly developed product, reliability is lacking due to failure data, etc. The accuracy of the evaluation may not be good. Any of the above patent documents also discloses the evaluation of reliability from the material of the object product and uses it for maintenance or application, but in contrast, it does not disclose or teach examples such as the use of products of the same type as the object product, similar articles. The reliability evaluation value is used to update the reliability evaluation value of the object product, and the product is planned and maintained according to the high-accuracy evaluation value.

本發明有鑑於以上問題點,目的在於達成可考量構成 製品之複數零件的高精度的可靠性評估值,來實施製品的運用、維護計劃。 The present invention has the above problems in view of the above problems, and aims to achieve a consideration. The high-precision reliability evaluation value of the multiple parts of the product is used to implement the operation and maintenance plan of the product.

按照本發明之第1解決手段,一種運作輔助裝置,具備:故障風險評估部;及維護/運用腳本制定部;前述故障風險評估部,利用從對象製品的複數個感測器輸入之環境資料及運轉資料、與事先訂定好之設計資料及材料資料,來演算視為對象之零件p的時刻t1下之破壞機率F(t1),於時刻t1,藉由零件p的時刻t1下的破壞機率F(t1)與當零件p損壞了的情形下之事先訂定好之每一零件p的影響度C(p)之乘積,來計算製品中包含之零件p的故障風險RS(t1,p),前述維護/運用腳本制定部,依據從前述故障風險評估部送來的時刻t1下之零件p的故障風險RS(t1、p)、及已從前述故障風險評估部送來而記憶於故障風險資料庫之從過去至時刻t1為止的複數個故障風險,而將事先由對象製品輸入的環境資料及運轉資料選擇出之對故障風險造成影響的物理量x、以及時間t訂為變數,生成故障風險的趨勢曲線,依據故障風險的趨勢曲線,求出從目前時間點至往前 了事先訂定好之時間的故障風險的預測值,依據故障風險的預測值,遵照依維護/運用者所做的來自輸入部之輸入而手動地設定出之、或藉由事先訂定好之處理而自動地設定出之製品的維護/運用腳本,而生成以時刻t、物理量x、以及由維護資料及/或運轉資料選擇出之對故障風險造成影響的物理量y作為變數之故障風險預測模型,來預測零件p的將來的故障風險,將每一零件的將來的故障風險,依故障風險的預測值較高之順序予以整理並顯示於顯示部及/或記憶於記憶部,並以事先訂定好之複數個閾值予以劃分群組,藉由依維護/運用者所做的來自輸入部之輸入而手動地、或藉由事先訂定好之處理而自動地對每一群組設定包含各零件的維護時期在內之維護運用內容,並將維護運用內容記憶於記憶部及/或顯示於顯示部。 According to a first aspect of the present invention, a service assisting device includes: a failure risk assessment unit; and a maintenance/operation scenario creation unit; and the failure risk assessment unit uses environmental information input from a plurality of sensors of the target product and The operation data, the design data and the material data set in advance, the failure probability F(t1) at the time t1 of the part p regarded as the object, and the probability of destruction at the time t1 of the part p at the time t1 (t1) Calculate the failure risk RS(t1, p) of the part p contained in the product by multiplying the influence degree C(p) of each part p which is determined in advance in the case where the part p is damaged. The maintenance/operation scenario development unit stores the failure risk data (t1, p) of the component p at the time t1 sent from the failure risk assessment unit, and the failure risk data that has been sent from the failure risk assessment unit. The plurality of failure risks from the past to the time t1, and the physical quantity x and the time t that affect the failure risk selected from the environmental data and the operation data input from the target product are set as variables, and the failure is generated. The trend curve of risk, based on the trend curve of the fault risk, is obtained from the current time point to the previous The predicted value of the failure risk set in advance is determined according to the predicted value of the failure risk, manually set according to the input from the input unit by the maintenance/operator, or by the predetermined processing. Automatically setting the maintenance/application script of the product, and generating a fault risk prediction model with the physical quantity y that is affected by the failure risk by the time t, the physical quantity x, and the maintenance data and/or the operation data as a variable. Predict the future failure risk of the part p, and sort out the future failure risk of each part according to the higher predicted value of the failure risk and display it on the display part and/or in the memory part, and set it in advance. The plurality of thresholds are divided into groups, and the maintenance period including each part is automatically set for each group manually by the input from the input unit by the maintenance/operator or manually by the predetermined processing. The maintenance content is used internally, and the maintenance application content is memorized in the memory unit and/or displayed on the display unit.

按照本發明之第2解決手段,提供一種風力發電系統,具備:如上述般之運作輔助裝置;及風力發電機,為具有複數個感測器之對象製品。按照本發明之第3解決手段,提供一種運作輔助裝置,具備:故障風險評估/更新部;及維護/運用腳本制定部;故障資料庫,蓄積構成對象製品與其同型機及/或類似機之複數個零件的故障資料; 前述故障風險評估/更新部,基於對象零件p達破壞之壽命的機率密度函數、及由前述故障資料庫中包含之故障資料計算出的概度,活用貝氏定理,求出考量了故障資料之更新後的壽命的機率密度函數,藉由更新後的壽命的機率密度函數,利用從對象製品的複數個感測器輸入之環境資料及運轉資料、與事先訂定好之設計資料及材料資料,來演算視為對象之零件p的時刻t1下之更新後的破壞機率F(t1)’,於時刻t1,藉由零件p的時刻t1下的更新後的破壞機率F(t1)’與當零件p損壞了的情形下之事先訂定好之每一零件p的影響度C(p)之乘積,來計算製品中包含之零件p的更新後的故障風險RS(t1,p)’,前述維護/運用腳本制定部,依據從前述故障風險評估/更新部送來的時刻t1下之零件p的更新後的故障風險RS(t1、p)’、及已從前述故障風險評估/更新部送來而記憶於故障風險資料庫之從過去至時刻t1為止的複數個更新後的故障風險,而將事先由對象製品輸入的環境資料及運轉資料選擇出之對故障風險造成影響的物理量x、以及時間t訂為變數,生成故障風險的趨勢曲線,依據故障風險的趨勢曲線,求出從目前時間點至往前之事先訂定好之時間的故障風險的預測值,依據故障風險的預測值,遵照依維護/運用者所做的 來自輸入部之輸入而手動地設定出之、或藉由事先訂定好之處理而自動地設定出之製品的維護/運用腳本,而生成以時刻t、物理量x、以及由維護資料及/或運轉資料選擇出之對故障風險造成影響的物理量y作為變數之故障風險預測模型,來預測零件p的將來的故障風險,並將預測值和趨勢曲線一起記憶於記憶部及/或顯示於顯示部。 According to a second solution of the present invention, there is provided a wind power generation system comprising: the operation aid device as described above; and the wind power generator, which is a target product having a plurality of sensors. According to a third aspect of the present invention, there is provided a operation assisting device comprising: a failure risk assessment/update unit; and a maintenance/operation script development unit; a fault data base for accumulating a plurality of constituent products and their identical machines and/or the like Fault data of parts; The fault risk assessment/update unit described above calculates the fault data based on the probability density function of the life of the target component p and the fault data calculated from the fault data contained in the fault database. The probability density function of the updated life is obtained by using the probability density function of the updated life, using environmental data and operation data input from a plurality of sensors of the target product, and design data and material materials set in advance. The calculation is regarded as the updated failure probability F(t1)' at time t1 of the part p of the object, and at time t1, the updated failure probability F(t1)' at time t1 of the part p and the part p In the case of damage, the product of the influence degree C(p) of each part p which is determined in advance is calculated to calculate the updated failure risk RS(t1, p)' of the part p contained in the product, the aforementioned maintenance/ The script creation unit is based on the updated failure risk RS(t1, p)' of the part p at the time t1 sent from the failure risk assessment/update unit, and has been sent from the failure risk assessment/update unit. Memory in faulty wind The plurality of updated failure risks from the past to the time t1 of the database, and the physical quantity x and the time t that are selected to affect the failure risk by the environmental data and the operation data input in advance by the target product are set as variables. The trend curve of the fault risk is generated, and the predicted value of the fault risk from the current time point to the previous predetermined time is obtained according to the trend curve of the fault risk, and according to the predicted value of the fault risk, according to the maintenance/operator made The maintenance/application script of the product that is manually set by input from the input unit or automatically set by the predetermined processing, and generated at time t, physical quantity x, and maintenance data and/or operation The data selects the physical quantity y that affects the risk of failure as a variable risk prediction model to predict the future failure risk of the part p, and memorizes the predicted value and the trend curve together in the memory and/or on the display.

按照本發明之第4解決手段,提供一種風力發電系統,具備:如上述般之運作輔助裝置;及第1風力發電機,為具有複數個感測器之對象製品;第2風力發電機,具有複數個感測器,和前述第1風力發電機為同型機或類似機。 According to a fourth aspect of the present invention, there is provided a wind power generation system comprising: the operation assisting device as described above; and the first wind power generator, which is a target product having a plurality of sensors; and the second wind power generator having A plurality of sensors are the same type of machine as the first wind turbine or the like.

按照本發明,能夠達成可考量構成製品之複數零件的高精度的可靠性評估值,來實施製品的運用、維護計劃。 According to the present invention, it is possible to realize a high-accuracy reliability evaluation value of a plurality of components constituting a product, and to implement an operation and maintenance plan of the product.

1‧‧‧對象製品 1‧‧‧ object products

2‧‧‧故障風險評估部 2‧‧‧Fault Risk Assessment Department

3‧‧‧維護/運用腳本制定部 3‧‧‧Maintenance/Using Script Development Department

4‧‧‧故障風險預測部 4‧‧‧Fault Risk Forecasting Department

5‧‧‧設計/材料資料庫 5‧‧‧Design/Materials Database

6‧‧‧破壞機率計算 6‧‧‧Destruction probability calculation

7‧‧‧風險計算 7‧‧‧ Risk calculation

8‧‧‧影響度資料庫 8‧‧‧Impact database

10‧‧‧故障風險趨勢分析 10‧‧‧ Failure risk trend analysis

11‧‧‧故障風險預測 11‧‧‧Fault risk prediction

12‧‧‧維護/運用腳本制定 12‧‧‧Maintenance/application scripting

13‧‧‧同型機、類似機 13‧‧‧The same type of machine, similar machine

14‧‧‧故障風險評估/更新部 14‧‧‧Fault Risk Assessment/Update Department

15‧‧‧故障資料庫 15‧‧‧ Fault Database

16‧‧‧等價應力振幅下之壽命的機率密度函數 16‧‧‧The probability density function of the lifetime under the equivalent stress amplitude

17‧‧‧外部資料庫 17‧‧‧External database

18‧‧‧破壞機率評估部 18‧‧‧Destruction of the probability assessment department

19‧‧‧損傷度評估部 19‧‧‧Damage Assessment Department

20‧‧‧疲勞壽命的機率密度函數 20‧‧‧ probability density function of fatigue life

21‧‧‧應力頻率分布 21‧‧‧stress frequency distribution

22‧‧‧破壞機率P%的疲勞壽命曲線 22‧‧‧Fatigue life curve with probability of failure P%

30‧‧‧故障風險資料庫 30‧‧‧Fault Risk Database

31‧‧‧破壞機率資料庫 31‧‧‧Destroy probability database

32‧‧‧損傷度資料庫 32‧‧‧Damage database

100‧‧‧運作輔助系統 100‧‧‧Operational support system

〔圖1〕依本發明第1實施形態之運作輔助系統的主要構成要素、及提供運作輔助系統中利用之資料的製品及資料庫、及它們的關係概略示意方塊圖。 Fig. 1 is a schematic block diagram showing the main components of the operation support system according to the first embodiment of the present invention, and the products and databases for providing the data used in the operation support system, and their relationships.

〔圖2〕依本發明第1實施形態之運作輔助系統當中,故障風險評估部的主要機能展開示意方塊圖。 [Fig. 2] A schematic block diagram showing the main functions of the failure risk assessment unit in the operation support system according to the first embodiment of the present invention.

〔圖3〕依本發明第1實施形態之運作輔助系統當 中,由對象製品中包含之零件的應力歷程而以剩餘壽命評估來演算破壞機率所必須之P-S-N線圖。 [Fig. 3] The operation assisting system according to the first embodiment of the present invention In the case of the stress history of the parts contained in the object product, the P-S-N line diagram necessary for the probability of destruction is calculated by the remaining life evaluation.

〔圖4〕依本發明第1實施形態之運作輔助系統當中,由在對象製品中包含之零件發生的應力頻率分布與P-S-N線圖而演算損傷度之方法示意圖。 [Fig. 4] A schematic diagram showing a method of calculating the degree of damage by the stress frequency distribution generated in the component included in the target product and the P-S-N line diagram in the operation assisting system according to the first embodiment of the present invention.

〔圖5〕依本發明第1實施形態之運作輔助系統當中,利用破壞壽命的機率密度函數來演算破壞機率之方法示意圖。 [Fig. 5] A schematic diagram of a method for calculating the probability of destruction using the probability density function of the damage life in the operation support system according to the first embodiment of the present invention.

〔圖6〕依本發明第1實施形態之運作輔助系統當中,由製品的計測資料來選擇和破壞有關連之物理量,並將其適當地變換而求出壽命的機率密度函數之方法示意圖。 [Fig. 6] A schematic diagram of a method for determining the probability density function of the life by selecting and destroying the physical quantity associated with the measurement data of the product according to the measurement data of the first embodiment of the present invention.

〔圖7〕依本發明第1實施形態之運作輔助系統當中,維護/運用腳本制定部的主要機能展開示意方塊圖。 [Fig. 7] A schematic block diagram showing the main functions of the maintenance/operation scenario creation unit in the operation support system according to the first embodiment of the present invention.

〔圖8〕依本發明第1實施形態之運作輔助系統當中,藉由故障風險預測部的故障風險趨勢分析而作成之故障風險的趨勢曲線的一例示意圖。 [Fig. 8] Fig. 8 is a schematic diagram showing an example of a trend curve of a failure risk generated by a failure risk trend analysis of a failure risk prediction unit in the operation support system according to the first embodiment of the present invention.

〔圖9〕依本發明第1實施形態之運作輔助系統當中,當更改藉由維護/運用腳本制定部執行之維護/運用腳本的情形下之風險預測處理的一例示意圖。 [Fig. 9] Fig. 9 is a diagram showing an example of risk prediction processing in the case where the maintenance/operation script executed by the maintenance/application scenario creation unit is changed in the operation support system according to the first embodiment of the present invention.

〔圖10〕依本發明第1實施形態之運作輔助系統當中,將製品的主要構成要素解釋成零件,依故障風險的預測值較高的順序予以整理顯示,並劃分群組的例子示意圖。 [Fig. 10] In the operation support system according to the first embodiment of the present invention, the main components of the product are explained as parts, and the example in which the predicted value of the failure risk is displayed in a high order is shown and the group is divided.

〔圖11〕依本發明第1實施形態之運作輔助系統當中,對隸屬製品的主要構成要素之每一零件演算故障風險的預測值,並依該預測值較高的順序整理零件而顯示之圖。 [Fig. 11] In the operation support system according to the first embodiment of the present invention, the predicted value of the failure risk is calculated for each of the main components of the subsidiary product, and the parts are displayed in accordance with the order in which the predicted values are high. Figure.

〔圖12〕依本發明第1實施形態之運作輔助系統當中,從破壞機率演算至維護/運用腳本制定為止之手續示意圖。 [Fig. 12] A schematic diagram of a procedure from the destruction probability calculation to the maintenance/application script development in the operation support system according to the first embodiment of the present invention.

〔圖13〕依本發明第2實施形態之運作輔助系統的主要構成要素、及對運作輔助系統提供資料的製品、同型機、類似機及資料庫、及它們的關係概略示意方塊圖。 Fig. 13 is a schematic block diagram showing the main components of the operation support system according to the second embodiment of the present invention, and the products, the same type machines, the similar machines and the data banks which provide information to the operation support system, and their relationships.

〔圖14〕依本發明第2實施形態之運作輔助系統當中,當由製品中包含之零件的應力歷程而以剩餘壽命評估來演算破壞機率的情形下,以等價應力振幅描繪壽命的機率密度函數,並依據貝氏統計將其更新之例子示意圖。 [Fig. 14] In the operation assisting system according to the second embodiment of the present invention, when the probability of failure is calculated by the residual life evaluation of the stress history of the components included in the product, the probability density of the life is plotted with the equivalent stress amplitude. A schematic diagram of the function and updating it based on Bayesian statistics.

〔圖15〕依本發明第2實施形態之運作輔助系統當中,當由製品中包含之零件的應力歷程而藉由剩餘壽命評估來演算破壞機率的情形下,從等價應力振幅演算至更新破斷壽命的機率密度函數為止之手續示意圖。 [Fig. 15] In the operation assisting system according to the second embodiment of the present invention, when the probability of failure is calculated by the residual life evaluation from the stress history of the components included in the product, the calculation from the equivalent stress amplitude to the update is broken. A schematic diagram of the procedure for the probability density function of the broken life.

〔圖16〕依本發明第2實施形態之運作輔助系統當中,當依據製品中包含之零件的破壞壽命的機率密度函數來演算破壞機率的情形下,依據貝氏統計來更新事前設定之壽命的機率密度函數之例子示意圖。 [FIG. 16] In the operation assisting system according to the second embodiment of the present invention, when the probability of failure is calculated based on the probability density function of the breaking life of the component included in the product, the life of the pre-set is updated according to the Bayesian statistics. A schematic diagram of an example of a probability density function.

〔圖17〕依本發明第2實施形態之運作輔助系統當中,當以包含時間之多變量的機率密度函數來表現破壞壽 命,而演算破壞機率的情形下,依據貝氏統計來更新事前設定之機率密度函數之例子示意圖。 [Fig. 17] In the operation assisting system according to the second embodiment of the present invention, the failure rate is expressed by a probability density function including a variable amount of time. In the case of a life-destroying probability, the example of the probability density function set in advance is updated according to the Bayesian statistics.

〔圖18〕依本發明第3實施形態之運作輔助系統的主要構成要素、及提供運作輔助系統中利用之資料的製品及資料庫、及它們的關係概略示意方塊圖。 Fig. 18 is a schematic block diagram showing the main components of the operation support system according to the third embodiment of the present invention, and the products and databases for providing the data used in the operation support system, and their relationships.

〔圖19〕依本發明第3實施形態之運作輔助系統的維護/運用腳本制定部的主要構成要素,與在要素間往來之資料的流向概略示意方塊圖。 Fig. 19 is a schematic block diagram showing main components of a maintenance/operation scenario creation unit of the operation support system according to the third embodiment of the present invention, and flow of information between elements.

〔圖20〕依本發明第4實施形態之運作輔助系統的主要構成要素、及提供運作輔助系統中利用之資料的製品及資料庫、及它們的關係概略示意方塊圖。 Fig. 20 is a schematic block diagram showing the main components of the operation support system according to the fourth embodiment of the present invention, and the products and databases for providing the data used in the operation support system, and their relationships.

〔圖21〕依本發明第4實施形態之運作輔助系統的維護/運用腳本制定部的主要構成要素,與在要素間往來之資料的流向概略示意方塊圖。 Fig. 21 is a schematic block diagram showing main components of the maintenance/operation scenario creation unit of the operation support system according to the fourth embodiment of the present invention, and the flow of information between the elements.

〔圖22〕依本發明第4實施形態之運作輔助系統當中,當作成破壞機率的趨勢曲線的情形下之維護/運用腳本制定部的主要構成要素,與在要素間往來之資料的流向概略示意方塊圖。 [Fig. 22] In the operation support system according to the fourth embodiment of the present invention, the main components of the maintenance/operation scenario creation unit in the case of the tendency curve of the destruction probability, and the flow of the data between the elements are schematically illustrated. Block diagram.

〔圖23〕依本發明第5實施形態之運作輔助系統的主要構成要素、及提供運作輔助系統中利用之資料的製品及資料庫、及它們的關係概略示意方塊圖。 Fig. 23 is a schematic block diagram showing the main components of the operation support system according to the fifth embodiment of the present invention, and the products and databases for providing the data used in the operation support system, and their relationships.

〔圖24〕依本發明第5實施形態之運作輔助系統當中,當作成損傷度的趨勢曲線的情形下之維護/運用腳本制定部的主要構成要素,與在要素間往來之資料的流向概 略示意方塊圖。 [Fig. 24] In the operation support system according to the fifth embodiment of the present invention, the main components of the maintenance/operation scenario creation unit in the case of the trend curve of the damage degree, and the flow of information between the elements. Slightly schematic block diagram.

A.概要 A. Summary

本實施形態包含複數個解決上述問題之手段,但若舉其一例,則能夠為一種任意製品的運作輔助系統, This embodiment includes a plurality of means for solving the above problems, but as an example, it can be an operation assisting system for an arbitrary product.

其特徵為:除了構成前述製品之複數個零件的故障風險,亦即包含從過去至現在之前述製品的環境資料、運轉資料、設計資料、材料資料的至少一者之資訊外,還依據前述製品的同型機、類似機的故障資料之資訊來執行故障風險的評估與更新,具備求出當於現在更改前述製品的維護/運用計劃的情形下會變動之複數個前述零件的故障風險推定值之演算裝置,具備參照演算出之複數個前述零件的故障風險推定值各者而分派前述製品的運用或複數個前述零件的維護時期之手段。 It is characterized in that, in addition to the risk of failure of a plurality of parts constituting the aforementioned product, that is, information including at least one of environmental materials, operation materials, design materials, and material materials of the aforementioned products from the past to the present, The information on the fault data of the same type machine and the similar machine is used to perform the evaluation and update of the fault risk, and the fault risk estimation value of the plurality of the aforementioned parts that may change when the maintenance/application plan of the product is changed now is determined. The calculation device includes means for assigning the use of the product or a plurality of maintenance periods of the plurality of parts to each of the estimated risk risk estimates of the plurality of parts.

按照本實施形態,可提供一種能夠制定高可靠性製品的運用、維護計劃,並進行穩定的運作之製品的運作輔助系統。 According to the present embodiment, it is possible to provide an operation support system for a product capable of setting up an operation and maintenance plan for a highly reliable product and performing stable operation.

B.運作輔助裝置及風力發電系統 B. Operation aids and wind power systems

以下,利用圖面說明本發明之實施形態。 Hereinafter, embodiments of the present invention will be described using the drawings.

〔第1實施形態〕 [First Embodiment]

以下,利用圖1~10,說明依本發明第1實施形態之 運作輔助系統。本實施形態中,作為製品1係舉出風力發電廠,但本發明及/或本實施形態的運用並非限定於風力發電廠。 Hereinafter, the first embodiment of the present invention will be described with reference to Figs. Operational assistance system. In the present embodiment, the product 1 is a wind power plant, but the operation of the present invention and/or the present embodiment is not limited to a wind power plant.

按照第1實施形態,是將構成製品之複數零件予以群組化,藉此,能夠達成可考量構成製品之複數零件的高精度的可靠性評估值,來實施製品的運用、維護計劃。 According to the first embodiment, the plurality of components constituting the product are grouped, whereby a highly accurate reliability evaluation value of a plurality of components constituting the product can be achieved, and the operation and maintenance plan of the product can be implemented.

圖1為本實施形態之製品的運作輔助系統的主要構成要素、及提供運作輔助系統中利用之資料的製品及資料庫、及它們的關係概略示意方塊圖。圖1所示之運作輔助系統100,具備故障風險評估部2、維護/運用腳本制定部3。維護/運用腳本制定部3中,包含故障風險預測部4,其將製品中包含之複數個零件訂為對象,並預測當更改維護/運用計劃的情形下變動之該零件的故障風險。圖中雖省略,但運作輔助系統100,具備輸入部、顯示部、對其他裝置的輸出部。 BRIEF DESCRIPTION OF THE DRAWINGS Fig. 1 is a schematic block diagram showing the main components of the operation support system of the product of the present embodiment, and the products and databases for providing the data used in the operation support system, and their relationships. The operation support system 100 shown in FIG. 1 includes a failure risk assessment unit 2 and a maintenance/operation scenario creation unit 3. The maintenance/application scenario creation unit 3 includes a failure risk prediction unit 4 that targets a plurality of parts included in the product and predicts the risk of failure of the part that changes when the maintenance/operation plan is changed. Although omitted in the figure, the operation support system 100 includes an input unit, a display unit, and an output unit for other devices.

製品1中組裝有多數個用來計測使用環境或運轉狀態之感測器。藉由該些感測器計測出的環境資料、運轉資料,會被送至運作輔助系統100的故障風險評估部2,而被利用於評估。此處所謂環境資料,為包含和製品所曝露之環境有關的資料,風力發電廠的情形下,例如風車的風速、風向等風況資料會包含在環境資料中。當為設置於海洋上之風力發電廠的情形下,除風況資料外,波長或波高等海況資料亦為環境資料的範疇。此外所謂運轉資料,為速度、加速度、旋轉速度、旋轉角等和製品之運作狀態關 連的資料。風力發電廠的情形下,風車的發電量、發電機的旋轉速度、方位角(azimuth)、機艙(nacelle)角等為運轉資料的範疇。廠房中,環境資料、運轉資料多半被計測作為控制計測(SCADA)。但,針對包含狀態監視(CMS)或構造物監視(SHM)之廠房,若藉由它們計測之資料中亦為和製品的使用環境或運作狀態關連的資料,則包含在本實施形態之環境資料或運轉資料。 A plurality of sensors for measuring the use environment or operating state are assembled in the product 1. The environmental data and the operational data measured by the sensors are sent to the failure risk assessment unit 2 of the operation support system 100, and are used for evaluation. The environmental information referred to here is information related to the environment exposed by the product. In the case of a wind power plant, wind conditions such as wind speed and wind direction of the windmill are included in the environmental data. In the case of a wind power plant installed on the ocean, in addition to the wind condition data, sea state data such as wavelength or wave height is also the category of environmental data. In addition, the so-called operating data is the speed, acceleration, rotation speed, rotation angle, etc. Even the information. In the case of a wind power plant, the amount of power generated by the windmill, the rotational speed of the generator, the azimuth angle, and the nacelle angle are the categories of operational data. In the plant, environmental data and operational data are mostly measured as control measurement (SCADA). However, for the plant containing state monitoring (CMS) or structure monitoring (SHM), if the data measured by them is also related to the use environment or operational status of the product, the environmental data included in the embodiment is included. Or operational data.

除環境資料或運轉資料外,設計資料或材料資料亦在故障風險評估部2中受到利用。此處設計資料,例如包含製品的圖面等和製品形狀有關的資料。此外,材料資料,包含構成製品之材料的特性、或螺栓緊固或熔接配管等之構造的特性。 In addition to environmental data or operational data, design data or material data are also used in the failure risk assessment department 2. Here, the design information, such as the drawing containing the product, and the shape of the product. Further, the material information includes characteristics of materials constituting the product, or characteristics of a structure such as a bolt fastening or a welding pipe.

圖2為故障風險評估部2中執行之處理概略示意方塊圖。 FIG. 2 is a schematic block diagram showing the processing executed in the failure risk assessment unit 2.

故障風險評估部2,利用環境資料、運轉資料、設計資料及材料資料的至少其中一者,演算故障風險。於時刻t1,製品中包含之某一零件p的故障風險RS(t1,p),能夠由對象零件p的時刻t1下之破壞機率F(t1)、及當對象零件p損壞了的情形下之影響度C(p),而依下式計算。另,每一零件p的影響度C(p)的資料,是事先蓄積於影響度資料庫8。 The failure risk assessment unit 2 calculates the failure risk using at least one of environmental data, operational data, design data, and material data. At time t1, the failure risk RS(t1, p) of a certain part p included in the product can be caused by the destruction probability F(t1) at the time t1 of the target part p, and when the target part p is damaged. The degree of influence C(p) is calculated according to the following formula. In addition, the information of the degree of influence C(p) of each part p is accumulated in advance in the influence degree database 8.

RS(t1,p)=C(p)×F(t1)‧‧‧(1) RS(t1,p)=C(p)×F(t1)‧‧‧(1)

將式(1)的故障風險RS(t1,p)訂為可靠性之指標,藉此,便能考量當某一零件故障了的情形下波及之影 響大小,來計劃零件之維護/運用。 The fault risk RS(t1,p) of equation (1) is set as an indicator of reliability, so that it can be considered when a certain part fails. Respond to the size, to plan the maintenance / use of parts.

破壞機率F(t1),能夠藉由剩餘壽命評估或徵兆診斷來計算。 The probability of destruction F(t1) can be calculated by residual life assessment or symptom diagnosis.

首先,說明剩餘壽命評估。 First, the remaining life assessment is explained.

將負荷反覆作用於製品而發生之疲勞現象視為對象,若欲藉由剩餘壽命評估來求出破壞機率F(t1),首先,利用環境資料、運轉資料及設計資料,計算從過去至目前時間點t1為止在零件發生的應力歷程。接著,對此應力歷程,運用雨流計數法(rainflow)等的頻率分析法,作成應力頻率分布,該應力頻率分布是將某一大小的應力是以怎樣的頻率發生予以整理而得。然後,利用應力頻率分布與和對象零件關連的材料資料,求出破壞機率F(t1)。 The fatigue phenomenon that occurs when the load is applied to the product is considered as the object. If the damage probability F(t1) is to be obtained by the remaining life assessment, first, the environmental data, the operational data, and the design data are used to calculate the time from the past to the present. The stress history of the part at point t1. Next, in this stress history, a frequency analysis method such as a rainflow method is used to create a stress frequency distribution obtained by sorting out a certain magnitude of stress. Then, using the stress frequency distribution and the material data associated with the target part, the probability of failure F(t1) is obtained.

圖3為由對象製品中包含之零件的應力歷程而以剩餘壽命評估來演算破壞機率所必須之P-S-N線圖。 Figure 3 is a P-S-N diagram necessary to calculate the probability of failure by the residual life assessment from the stress history of the parts contained in the target article.

此處利用的材料資料,理想是被稱為P-S-N線圖之如圖3的疲勞壽命曲線,本實施例中係蓄積於設計/材料資料庫5。P-S-N線圖,為由依縱軸的各應力振幅實施疲勞試驗而得到之疲勞壽命的機率密度函數20來求出達破壞機率P%之反覆次數,並將它們予以連結圖示而成者(圖3)。 The material data used herein is ideally the fatigue life curve of Fig. 3, which is called a P-S-N line diagram, and is accumulated in the design/material database 5 in this embodiment. The PSN line diagram is a probability density function 20 of the fatigue life obtained by performing a fatigue test on each of the stress amplitudes of the vertical axis, and the number of times of the failure probability P% is obtained, and these are connected as shown in the figure (Fig. 3). ).

圖4為由在對象製品中包含之零件發生的應力頻率分布與P-S-N線圖而演算損傷度之方法示意圖。 Fig. 4 is a schematic view showing the method of calculating the degree of damage by the stress frequency distribution and the P-S-N line diagram of the parts included in the object product.

如圖3般,當能夠假定各應力振幅下之疲勞壽命的機 率密度函數20全部為同一的情形下,可依下述手續求出破壞機率F(t1)。首先,如圖4所示,能夠由將頻率分析法運用至從過去至目前時間點t1為止發生的應力歷程而得到之應力頻率分布21、及事先訂定好之破壞機率P%的疲勞壽命曲線22,而依下式計算相對於疲勞破壞P%而言之損傷度D(t1)。 As shown in Figure 3, when it is possible to assume the fatigue life of each stress amplitude When the rate density function 20 is all the same, the probability of destruction F(t1) can be obtained by the following procedure. First, as shown in FIG. 4, the stress frequency distribution 21 obtained by applying the frequency analysis method to the stress history occurring from the past to the current time point t1, and the fatigue life curve 22 having the predetermined breaking probability P% can be set. And the damage degree D(t1) with respect to the fatigue damage P% is calculated according to the following formula.

D(t1)=(n1/N1)+(n2/N2)+…+(nm/Nm) ‧‧‧(2) D(t1)=(n 1 /N 1 )+(n 2 /N 2 )+...+(n m /N m ) ‧‧‧(2)

此處,n1、n2、nm,分別為藉由應力歷程的頻率分析而得到之應力振幅S1、S2,Sm的反覆次數(m為整數)。此外,N1,N2,Nm,分別為當反覆負荷應力振幅S1、S2,Sm的情形下,在破壞機率P%下發生疲勞破壞之破斷反覆次數。當假定疲勞壽命的機率密度函數20和應力振幅無關而是全部同一的情形下,導致損傷度D(t1)發生之反覆次數N(t1)可依下式求出。 Here, n 1 , n 2 , and n m are the repetition times (m is an integer) of the stress amplitudes S 1 , S 2 , and S m obtained by frequency analysis of the stress history, respectively. Further, N 1 , N 2 , and N m are the number of times of breaking and the occurrence of fatigue fracture at the failure probability P% when the load stress amplitudes S 1 , S 2 , and S m are repeated. When it is assumed that the probability density function 20 of the fatigue life is not the same as the stress amplitude, the number of times of the occurrence of the damage degree D(t1) N(t1) can be obtained by the following equation.

N(t1)=D(t1)×Np‧‧‧(3) N(t1)=D(t1)×Np‧‧‧(3)

Np為事先訂定好之破壞機率P%的反覆次數,藉由輸入部(未圖示)將應力振幅訂定為適當的Si(例如,平均值、中間值、接近該些的值等),藉此如圖4中所示,能夠由N(t1)與疲勞壽命的機率密度函數20,藉由圖4的式子,求出破壞機率F(t1)。也就是說,若將從製品運轉開始至目前為止之經過時間訂為t1,則破壞機率F(t1)為將密度函數f(N)從0至N(t1)為止積分而得到的值。 Np is the number of times of the failure probability P% set in advance, and the stress amplitude is set to an appropriate Si (for example, an average value, an intermediate value, a value close to the value, etc.) by an input unit (not shown). As shown in FIG. 4, the probability of failure F(t1) can be obtained from the probability density function 20 of N(t1) and fatigue life by the equation of FIG. In other words, if the elapsed time from the start of the product operation to the present is set to t1, the probability of destruction F(t1) is a value obtained by integrating the density function f(N) from 0 to N(t1).

接著,說明剩餘壽命評估的其他例子。 Next, another example of the remaining life evaluation will be described.

圖5為利用破壞壽命的機率密度函數來演算破壞機率之其他方法示意圖。 Figure 5 is a schematic diagram of another method for calculating the probability of failure using the probability density function of the damage life.

此外,達破壞為止的破壞機率F(t1),亦能直接定義如圖5般之壽命的機率密度函數f(t1)而求出。也就是說,若將從製品運轉開始至目前為止之經過時間訂為t1,則破壞機率F(t1)為將密度函數f(t)從0至t1為止積分而得到的值。此處,若將達破壞機率50%之時間定義為壽命T,則能將目前時間點t1與壽命時間T的差想成是剩餘壽命,因此像這樣直接定義壽命的機率密度函數之方法亦能解釋為剩餘壽命評估。另,壽命的機率密度函數f(t1)之資料,例如能夠對每一對象製品事先記憶於故障風險評估部2的內部記憶體等適當的記憶體。此外,壽命的機率密度函數f(t1)之資料,能夠使用環境資料、運轉資料等而事先訂定。 In addition, the probability of failure F(t1) up to the damage can be obtained by directly defining the probability density function f(t1) of the life as shown in FIG. In other words, if the elapsed time from the start of the product operation to the present is set to t1, the probability of destruction F(t1) is a value obtained by integrating the density function f(t) from 0 to t1. Here, if the time to reach the failure rate of 50% is defined as the life T, the difference between the current time point t1 and the life time T can be regarded as the remaining life. Therefore, the method of directly defining the probability density function of the life can also be used. Interpreted as residual life assessment. Further, the information of the probability density function f(t1) of the life can be stored in an appropriate memory such as the internal memory of the failure risk assessment unit 2 in advance for each target product. In addition, the data of the probability density function f(t1) of the life can be determined in advance using environmental data, operation data, and the like.

接著,說明徵兆診斷。 Next, the symptom diagnosis will be explained.

圖6為由製品的計測資料來選擇和破壞有關連之物理量,並將其適當地變換而求出壽命的機率密度函數之方法示意圖。 Fig. 6 is a schematic view showing a method for selecting and destroying the physical quantity associated with the measured data of the product and appropriately converting the probability density function of the life.

在一種方式的徵兆診斷中,事先整理對象零件正常時與異常時之運轉資料,並監視目前時間點的運轉資料來判定對象零件的故障。若將異常時的運轉資料群如圖6般理解成破壞的機率密度函數,則藉由將目前時間點的運轉資料的位置描繪於圖6,便能求出目前時間點t1下之破壞機 率F(t1)。有由異常時的運轉資料來事先選擇和著眼的零件的破壞有關連之物理量,並作成以選擇出的物理量作為機率變數之機率密度函數的情形,或亦有如圖6般,將運轉資料中包含之物理量變換成和著眼的零件的破壞有關連之形式,並以該變換物理量x1’、x2’作為機率變數來作成機率密度函數之情形。此處,所謂變換物理量,例如可舉出將運轉資料中包含之加速度資料做高速傅立葉變換而求出的加速度頻譜當中,某一特定的頻率的頻譜值。另,機率密度函數,例如能夠對每一對象製品事先記憶於故障風險評估部2的內部記憶體或設計/材料資料庫5等適當的記憶體。此外,機率密度函數,能夠使用環境資料、運轉資料等而事先訂定。 In the symptom diagnosis of one type, the operation data of the normal and abnormal time of the target part is sorted in advance, and the operation data at the current time point is monitored to determine the failure of the target part. If the operational data group at the time of the abnormality is understood as the probability density function of the failure as shown in FIG. 6, by describing the position of the operation data at the current time point in FIG. 6, the destruction machine at the current time point t1 can be obtained. Rate F(t1). There is a physical quantity related to the destruction of the part selected and focused by the operation data at the time of the abnormality, and the physical quantity is selected as the probability density function of the probability variable, or as shown in Fig. 6, the operation data is included. The physical quantity is converted into a form related to the destruction of the component under consideration, and the converted physical quantity x1', x2' is used as the probability variable to create a probability density function. Here, the physical quantity to be converted includes, for example, a spectral value of a specific frequency among acceleration spectra obtained by performing fast Fourier transform on acceleration data included in the operation data. Further, the probability density function can be stored, for example, in an appropriate memory such as the internal memory of the failure risk assessment unit 2 or the design/material database 5 for each target product. In addition, the probability density function can be determined in advance using environmental data, operation data, and the like.

接著,說明C(p)。 Next, C(p) will be described.

影響度資料庫8中蓄積著式(1)中包含之每一零件p的影響度C(p)的資料。作為影響度C(p),若採用當實際上某一零件p損壞了的情形下所需的全部或事先訂定好之範圍等的費用,則故障風險RS(t1,p),能夠想成是於目前時間點t1因故障而肇生的全部或事先訂定好之範圍等的損失費用的期望值。某一零件損壞了的情形下所需的費用中,包含新零件本體的費用、零件的更換費用、零件的搬運費用、製品的運作停止所致之發電機會的損失費用等。另一方面,亦可不是具體的費用,而是因應當零件損壞了的情形下所肇生之影響的大小來分派一整數,並將該整數採用作為影響度C(p)。這樣的情形 下,能夠在零件間相對比較故障風險RS(t1,p),而判斷應留意哪一零件的可靠性。 The influence degree database 8 stores information on the degree of influence C(p) of each part p included in the equation (1). As the degree of influence C(p), if the cost of all or a predetermined range required in the case where a certain part p is actually damaged is used, the risk of failure RS(t1, p) can be considered as It is an expected value of the loss cost of all or a predetermined range which is generated by the failure at the current time point t1. The cost of a certain part is damaged, including the cost of the new part body, the replacement cost of the part, the handling cost of the part, and the loss of the generator due to the operation stoppage of the product. On the other hand, instead of a specific cost, an integer may be assigned due to the influence of the influence of the damage of the part, and the integer is adopted as the influence degree C(p). Such a situation In the next case, it is possible to compare the failure risk RS(t1, p) between the parts and determine which part of the reliability should be noted.

如此演算出的故障風險RS(t1,p),會和維護/運用資料、環境資料、運轉資料、設計資料、材料資料一起被輸入至維護/運用腳本制定部3。 The fault risk RS (t1, p) calculated in this way is input to the maintenance/application script development unit 3 together with the maintenance/application data, the environmental data, the operation data, the design data, and the material data.

圖7為維護/運用腳本制定部中執行之處理概略示意方塊圖。維護/運用腳本制定部3,具有故障風險預測部4,故障風險預測部4中,是由從故障風險評估部2送來的時刻t1下之故障風險RS(t1、p)、及已從故障風險評估部2送來而被事先記憶於故障風險資料庫30之從過去至時刻t1為止的各時刻tn(n=0、-1、-2、-3、‧‧‧)的複數個故障風險RS(tn、p)、及環境資料、運轉資料,來分析故障風險的變動趨勢。 Fig. 7 is a schematic block diagram showing the processing executed in the maintenance/operation scenario creation unit. The maintenance/application scenario creation unit 3 includes a failure risk prediction unit 4, and the failure risk prediction unit 4 is the failure risk RS(t1, p) at the time t1 sent from the failure risk assessment unit 2, and the failure fault. The plurality of failure risks sent by the risk assessment unit 2 and memorized in advance in the failure risk database 30 from the past to the time t1 (n=0, -1, -2, -3, ‧‧) RS (tn, p), and environmental data, operational data, to analyze the trend of failure risk.

圖8為藉由故障風險預測部的故障風險趨勢分析而作成之故障風險的趨勢曲線的一例示意圖。 FIG. 8 is a diagram showing an example of a trend curve of a failure risk generated by a failure risk trend analysis of a failure risk prediction unit.

例如,位於日本的風力發電廠的情形下,冬季有風較強的傾向,因此依零件不同,如圖8般能夠描繪將橫軸訂為時間、縱軸訂為故障風險之趨勢曲線。但,即使有季節變動,依破壞機率的演算中使用之P-S-N線或機率密度函數的分布形狀不同,無法得到如圖8般有規則性的趨勢曲線之情形亦不少。又,即使是構成同一風力發電廠的零件,短期性的風的紊亂程度所致之故障風險RS(t1、p)的變動,仍有和季節變動同等或是比季節變動還大之情形,因此必須對每一零件由環境資料或運轉資料來選擇參 與趨勢變動之物理量,除時間外還將選擇出的物理量訂為變數來描繪故障風險的趨勢曲線。也就是說,故障風險預測部4,若將時間訂為t、選擇出的物理量訂為x,則因應製品或構成製品之零件,選擇下式(4)、(5)等來決定故障風險的趨勢曲線RS。 For example, in the case of a wind power plant in Japan, there is a tendency for winds to be strong in winter. Therefore, depending on the components, as shown in Fig. 8, a trend curve in which the horizontal axis is set as time and the vertical axis is set as a risk of failure can be drawn. However, even if there is a seasonal change, the distribution shape of the P-S-N line or the probability density function used in the calculation of the probability of destruction is different, and there are many cases where a regular trend curve as shown in Fig. 8 cannot be obtained. In addition, even if it is a component that constitutes the same wind power plant, the fluctuation risk RS (t1, p) due to the short-term wind disturbance is still the same as the seasonal change or larger than the seasonal change. Each part must be selected from environmental data or operational data. The physical quantity with the trend change, in addition to the time, the selected physical quantity is set as a variable to describe the trend curve of the failure risk. In other words, when the time limit is set to t and the selected physical quantity is set to x, the failure risk prediction unit 4 selects the following items (4), (5), etc. to determine the risk of failure in response to the product or the component constituting the product. Trend curve RS.

RS(t,p)=g1(t,p)‧‧‧(4) RS(t,p)=g1(t,p)‧‧‧(4)

RS(t,x,p)=g2(t,x,p)‧‧‧(5) RS(t,x,p)=g2(t,x,p)‧‧‧(5)

式(4)(5),是以時刻t或物理量x為變數之故障風險的趨勢曲線,但亦能採用含有從遡及某一時間點之過去至現在為止的物理量或誤差項之自迴歸移動平均模型,或以時刻、環境資料、運轉資料作為輸入而學習而成之神經網路(neural network)來作為趨勢曲線。 Equation (4)(5) is a trend curve of the risk of failure with time t or physical quantity x as a variable, but it can also adopt an autoregressive moving average containing physical quantities or error terms from the past and the present time to a certain point in time. A model, or a neural network that is learned by inputting time, environmental data, and operational data as a trend curve.

另,作為簡單的數式,例如有RS(t,x,p)=α(p)‧t+β(p)‧x+γ(p)(α,β,γ為常數,惟和零件有關)等,但不限於此。 In addition, as a simple formula, for example, RS(t, x, p) = α(p) ‧ t + β (p) ‧ x + γ (p) (α, β, γ is a constant, but related to the part ), etc., but not limited to this.

故障風險預測部4,依據故障風險的趨勢曲線、維護/運用資料、維護/運用腳本制定12中訂定之維護/運用腳本,來實施故障風險預測11。所謂維護/運用資料,為過去之資料,例如包含製品的定期檢驗之資訊、用以更改運用之控制變更資訊、因故障而做的檢驗實施之資訊。所謂維護/運用腳本,例如指何年何月更換哪一零件、以怎樣方式運用等這類計劃。維護/運用腳本制定12,可藉由事先訂定好之手法自動地作成維護/運用腳本,亦可從輸入部手動地輸入維護/運用腳本。 The failure risk prediction unit 4 implements the failure risk prediction 11 based on the trend curve of the failure risk, the maintenance/operation data, and the maintenance/application script specified in the maintenance/application script. The so-called maintenance/use information is information about the past, such as information on periodic inspections of products, information on changes in control used to change operations, and inspections performed on failures. The so-called maintenance/application scripts, for example, refer to plans such as which parts are replaced in what month and how they are used. The maintenance/application script development 12 can be automatically created as a maintenance/application script by a predetermined method, or a maintenance/application script can be manually input from the input unit.

圖9為當更改藉由維護/運用腳本制定部執行之維護/運用腳本的情形下之風險預測處理的一例示意圖。 FIG. 9 is a diagram showing an example of risk prediction processing in the case where the maintenance/operation script executed by the maintenance/application scenario development unit is changed.

當於將來亦會繼續目前採用的定期檢驗、運用方法(製品的控制方法等)的情形下,從目前時間點至往前之事先訂定好之時間△的將來的故障風險(預測風險)a,會如圖9中的預測風險a般,成為沿著過往的趨勢曲線之物。 In the case where the current periodic inspection and application method (product control method, etc.) will be continued in the future, the future failure risk (predicted risk) a from the current time point to the prior time is set in advance. It will become the trend curve along the past as shown by the predicted risk a in Figure 9.

接著,如圖7所示,維護/運用腳本制定部3,使用預測風險進行維護/運用腳本制定12。也就是說,依據預測風險的值,藉由維護運用者等所做的來自輸入部之輸入而手動地、或藉由或事先訂定好之處理而自動地檢討製品之維護或運用方法之變更。例如,當預測風險比事先訂定好之閾值還高的情形下,維護/運用腳本制定12會提議變更成更加平穩的製品之運用方法,當預測風險比事先訂定好之閾值還低的情形下,會提議更加激烈的製品之運用方法或頻繁的維護檢驗。提議出的維護/運用腳本再度被送至故障風險預測11,因應該維護/運用腳本,如圖9中般預測將來風險b或c。為了預測當更改維護或運用的情形下之將來風險b或c,除時刻t、對故障風險造成影響之物理量x(由環境資料、運轉資料選擇)外,還從維護/運轉資料選擇對故障風險造成影響之物理量y,而能夠利用以時刻t、物理量x及y作為變數之下式的故障風險預測模型g3或g4。 Next, as shown in FIG. 7, the maintenance/operation scenario creation unit 3 performs maintenance/operation script development 12 using the predicted risk. That is to say, based on the value of the predicted risk, the maintenance or the change of the method of the product is automatically reviewed by the maintenance of the input from the input unit by the operator or the like manually or by a predetermined process. For example, when the predicted risk is higher than the pre-defined threshold, the maintenance/application scripting 12 will propose a change to a more stable method of use, when the predicted risk is lower than the pre-defined threshold. More intense use of the product or frequent maintenance inspections will be proposed. The proposed maintenance/application script is again sent to the failure risk prediction 11 as the future risk b or c is predicted as in Figure 9. In order to predict the future risk b or c in the case of changing maintenance or operation, in addition to the physical quantity x (selected by environmental data and operational data) that affects the risk of failure, the risk of failure is also selected from the maintenance/operation data. The physical quantity y of the influence is caused, and the failure risk prediction model g3 or g4 having the time t, the physical quantity x, and y as the variable expression can be utilized.

RS(t+△T,y,p)=g3(t+△T,y,p)‧‧‧(6) RS(t+△T,y,p)=g3(t+△T,y,p)‧‧‧(6)

RS(t+△T,x,y,p)=g4(t+△T,x,y,p)‧‧‧(7) RS(t+△T,x,y,p)=g4(t+△T,x,y,p)‧‧‧(7)

式(6)(7),是以時刻t+△T或將來的物理量x、y為變數之趨勢曲線,但亦能採用含有從遡及某一時間點之過去至時刻t+△T為止的物理量或誤差項之自迴歸移動平均模型、或以時刻、環境資料、運轉資料、維護/運用資料作為輸入而學習而成之神經網路來作為預測模型。 Equations (6) and (7) are trend curves with time t + ΔT or future physical quantities x and y as variables, but physical quantities or errors including from 遡 and a certain time point to time t + ΔT can also be used. The autoregressive moving average model of the item or the neural network learned by inputting time, environmental data, operational data, and maintenance/utilization data as a predictive model.

接著,圖10為將製品的主要構成要素解釋成零件,依故障風險的預測值較高的順序予以整理顯示,並劃分群組的例子示意圖。考量群組化來進行圖7的維護/運用腳本制定12。例如,制定哪一群組於幾年後維護等計劃。 Next, FIG. 10 is a schematic diagram showing an example in which the main components of the product are interpreted as parts, and the predicted values of the failure risk are ranked in a high order, and the group is divided. Consider the grouping to perform the maintenance/application script formulation 12 of FIG. For example, plan which group to maintain and other plans after a few years.

圖9的故障風險預測,可將構成製品的全部零件視為對象來實施。但,僅將故障風險(或,預測風險或將來風險)較高的零件視為對象來預測故障風險,藉此能夠削減預測所必須之工程數,能夠有效率地計劃維護或運用。維護/運用腳本制定部3,針對故障風險的預測值,例如如圖10般,是將預測值成為最大的零件、和其最大值一起,從故障風險(或,預測風險或將來風險)的預測值較高者予以依序整理,顯示於顯示部。針對群組化的具體處理,例如,能夠藉由每一群組的上限及下限之閾值,來分派群組及維護時期,以使全部的零件的故障風險不超過某一閾值。例如,藉由基因演算法(genetic algorithm),可達成這樣的分派。維護/運用腳本制定部3的維護/運用腳本制定12中,參照圖10,藉由維護運用者等所做的 來自輸入部之輸入而手動地、或藉由事先訂定好之處理而自動地,對每一群組計劃何時維護哪一零件等維護運用內容。或是,亦可從輸入部手動地輸入維護運用內容等維護/運用腳本。維護/運用腳本制定部3的維護/運用腳本制定12中,會將計劃出的維護運用內容記憶於適當的記憶部及/或顯示於顯示部。若故障風險的預測值(預測風險或將來風險等)是因應大小來整理,則例如如圖10所示,能夠藉由事先訂定好之複數個閾值,將零件群予以群組化成群組A、群組B、群組C等來計劃維護。藉由這樣的群組化所致之維護計劃,能夠配合事先計劃好的維護排程,進行故障風險的預測值(預測風險或將來風險)較高的零件之更換等,例如,在已排定距當前最近的維護時期實施隸屬群組A之零件的更換等。用以分派這樣的維護時期之群組化,能夠利用基因演算法或分支界限法(branch and bound)等的組合最佳化法來實施。 The failure risk prediction of Fig. 9 can be implemented by considering all the components constituting the product as objects. However, only parts with high risk of failure (or predicted risk or future risk) are considered as objects to predict the risk of failure, thereby reducing the number of projects necessary for prediction and efficiently planning maintenance or operation. The maintenance/application scenario development unit 3 predicts the risk of failure, for example, as a component having the largest predicted value, together with the maximum value thereof, as predicted from the failure risk (or predicted risk or future risk). Those with higher values are sorted in order and displayed on the display. For specific processing of grouping, for example, the group and maintenance period can be assigned by the threshold of the upper and lower limits of each group, so that the risk of failure of all parts does not exceed a certain threshold. Such assignments can be achieved, for example, by a genetic algorithm. The maintenance/application script creation 12 of the maintenance/application scenario creation unit 3 is performed by the maintenance operator, etc. with reference to FIG. From the input of the input unit, the maintenance operation content is automatically scheduled for each group, such as which part is scheduled, manually or by a predetermined process. Alternatively, a maintenance/application script such as maintenance operation contents may be manually input from the input unit. In the maintenance/operation scenario creation 12 of the maintenance/application scenario creation unit 3, the planned maintenance operation content is stored in an appropriate storage unit and/or displayed on the display unit. If the predicted value of the risk of failure (predicted risk or future risk, etc.) is sorted according to the size, for example, as shown in FIG. 10, the group of parts can be grouped into group A by setting a plurality of thresholds in advance. Group B, group C, etc. plan maintenance. With such a maintenance plan due to grouping, it is possible to match the planned maintenance schedule with the pre-planned maintenance schedule, and to replace the parts with higher predicted risk (predicted risk or future risk), for example, in the scheduled Replacement of parts belonging to group A is performed from the current maintenance period. The grouping of such maintenance periods can be implemented by a combination optimization method such as a gene algorithm or a branch and bound method.

此外,圖11為對隸屬製品的主要構成要素之每一零件演算故障風險的預測值,並依該預測值較高的順序整理零件而顯示之圖。 In addition, FIG. 11 is a diagram showing the predicted values of the risk of failure of each component of the main components of the belonging product, and sorting the parts in the order in which the predicted values are higher.

上述圖10,是將製品中包含之組件解釋成零件,並將零件與故障風險的預測值之關係予以整理而成之物。例如,將訂定哪一組件中包含哪一零件之檔案事先記憶於適當的記憶部,維護/運用腳本制定部3,能夠參照它而將組件、零件、故障風險顯示於顯示部。像這樣將組件視為對象來整理故障風險的預測值(預測風險或將來風險), 藉此能夠將對製品造成影響之組件明確化。另一方面,如圖11般,亦能以組件中包含之零件為單位來整理故障風險的預測值。像這樣以零件為單位來整理故障風險的預測值,藉此能夠使做維護所應準備之零件變得明確。為了獲得這樣的效果,本實施形態之運作輔助系統,如圖10、11般具備整理並顯示故障風險的預測值之裝置。 Figure 10 above is a figure in which the components included in the product are interpreted as parts, and the relationship between the parts and the predicted value of the failure risk is sorted. For example, the file of which component is included in the predetermined memory is stored in the appropriate memory unit, and the maintenance/application scripting unit 3 can display the component, the component, and the risk of failure on the display unit with reference to it. Treat components as objects to sort out the predicted value of the risk of failure (predicted risk or future risk), This makes it possible to clarify the components that affect the product. On the other hand, as shown in Fig. 11, the predicted value of the failure risk can also be sorted in units of parts included in the assembly. In this way, the predicted value of the failure risk is sorted in units of parts, whereby the parts to be prepared for maintenance can be made clear. In order to obtain such an effect, the operation assisting system of the present embodiment has means for sorting and displaying the predicted value of the risk of failure as shown in Figs.

運作輔助系統100的故障風險評估部2中的故障風險評估、維護/運用腳本制定部3中的故障風險預測、維護/運用腳本制定,是以某一時間間隔執行。此時間間隔可和計測環境資料或運轉資料之時間間隔同一,亦可相異。風力發電廠中採用之控制計測(SCADA)中,環境資料或運轉資料的統計值(最大值、最小值、平均值等),例如是以10分間隔演算,該統計值蓄積於由PC等所構成之伺服器。例如,訂定配合此10分間隔來執行故障風險評估、故障風險預測,將故障風險預測的△T訂為例如數個月,則藉由本實施形態之運作輔助系統,能夠充分地預測風力發電廠中包含之零件的故障,同時使風車穩定運作。 The failure risk assessment in the failure risk assessment unit 2 of the operation support system 100, the failure risk prediction in the maintenance/operation scenario creation unit 3, and the maintenance/operation script creation are performed at certain time intervals. This time interval can be the same as or different from the time interval between the measurement environment data and the operation data. In the control measurement (SCADA) used in a wind power plant, the statistical values (maximum value, minimum value, average value, etc.) of the environmental data or the operational data are calculated, for example, at intervals of 10 minutes, and the statistical values are accumulated in the PC or the like. The server that constitutes it. For example, if the fault risk assessment and the fault risk prediction are performed in accordance with the 10-point interval, and the fault risk prediction ΔT is set to, for example, several months, the wind power plant can be sufficiently predicted by the operation assisting system of the present embodiment. The failure of the parts contained in it also stabilizes the windmill.

圖12揭示依第1實施形態之運作輔助系統的流程圖。首先,故障風險評估部2,演算視為對象之零件p的時刻t1下之破壞機率F(t1)及故障風險RS(t1,p)(S11、S12)。接著,故障風險預測部4,演算故障風險的趨勢曲線RS(t,x,p)及將來的故障風險RS(t+△T,x,y,p)(S13、S14)。然後,依據計算出的將來的故障風險,藉由與事先訂定好之閾值比較等而自動地、 或藉由以維護/運用者等的手動而來自輸入部的輸入,來判斷是否變更維護/運用。然後,若判斷為有變更的必要,則藉由事先訂定好之處理而自動地、或藉由維護/運用者等所做的來自輸入部之輸入而手動地設定/輸入對維護/運用腳本制定12變更之物理量/條件的參數等,維護/運用腳本制定部3遵照該設定/輸入來制定維護/運用腳本(S15、S16、S14)。這樣的手續,如前述般例如以10分間隔反覆,來輔助製品的穩定運作。 Fig. 12 is a flow chart showing the operation assisting system according to the first embodiment. First, the failure risk assessment unit 2 calculates the failure probability F(t1) and the failure risk RS(t1, p) at the time t1 of the component p which is the target (S11, S12). Next, the failure risk prediction unit 4 calculates a trend curve RS(t, x, p) of the failure risk and a future failure risk RS (t + ΔT, x, y, p) (S13, S14). Then, based on the calculated future risk of failure, by automatically comparing with a predetermined threshold, etc. Alternatively, it is determined whether or not the maintenance/operation is changed by input from the input unit by manual operation such as maintenance/operator. Then, if it is determined that there is a need for the change, the maintenance/application script is manually set/inputted automatically by the predetermined processing or by the input from the input unit by the maintenance/operator or the like. The maintenance/application scenario creation unit 3 creates a maintenance/operation script in accordance with the setting/input (the parameters of the physical quantity/condition of the change) (S15, S16, and S14). Such procedures are repeated at intervals of, for example, 10 minutes as described above to assist in the stable operation of the article.

〔第2實施形態〕 [Second Embodiment]

接著,利用圖13~19,說明依本實施形態第2實施形態之運作輔助系統。 Next, an operation assisting system according to a second embodiment of the present embodiment will be described with reference to Figs.

按照第2實施形態,能夠達成使用和對象製品同型的製品、類似的製品之可靠性評估值,並考量構成製品之複數零件的高精度的可靠性評估值,來實施製品的運用、維護計劃。 According to the second embodiment, it is possible to achieve a reliability evaluation value of a product similar to the target product and a similar product, and to evaluate the high-accuracy reliability evaluation value of the plurality of components constituting the product, and to implement the operation and maintenance plan of the product.

本實施形態,如圖13所示,除了含有從過去至現在之製品的環境資料、運轉資料、設計資料、材料資料的至少任一者之資訊外,還利用製品1與其同型機、類似機13的故障資料之資訊,而藉由故障風險評估/更新部14執行故障風險的評估及更新。構成製品1與其同型機、類似機13之複數個零件的故障資料,係蓄積於故障資料庫15。此處所謂故障資料,例如包含從運轉開始至故障為止的運作時間、或從運轉開始至故障時為止的環境資料、運 轉資料、維護/運用資料。此外所謂同型機,係包含和製品1為相同型式的製品惟運作的場所相異之製品,所謂類似機,係包含和製品1為不同型式之製品。 In the present embodiment, as shown in FIG. 13, in addition to the information of at least one of the environmental data, the operation data, the design data, and the material data of the products from the past to the present, the product 1 and its like machine and the similar machine 13 are used. The fault information is updated, and the fault risk assessment/update unit 14 performs the evaluation and update of the fault risk. The failure data of the plurality of components constituting the product 1 and its homogenizer and similar machine 13 is stored in the fault database 15. The fault data here includes, for example, the operation time from the start of the operation to the failure, or the environmental data and operation from the start of the operation to the time of the failure. Transfer data, maintain/use data. In addition, the so-called homogenizer is a product which is different from the place where the product of the product 1 is of the same type, and the so-called similar machine comprises a product of a different type from the product 1.

圖14為當由製品中包含之零件的應力歷程而以剩餘壽命評估來演算破壞機率的情形下,以等價應力振幅描繪壽命的機率密度函數,並依據貝氏統計將其更新之例子示意圖。 Fig. 14 is a diagram showing an example of a probability density function for depicting the life with an equivalent stress amplitude when the probability of failure is calculated by the residual life evaluation of the component included in the article, and updating it according to the Bayesian statistics.

故障風險評估/更新部14中,是將在製品反覆承受負荷而發生之疲勞現象視為對象,當藉由剩餘壽命評估來求出破壞機率F(t1)的情形下,將雨流計數法等頻率分析法運用於從過去至目前時間點t1為止在零件發生的應力歷程,求出應力頻率分布(參照圖4)。實際上,應力歷程中包含各式各樣大小的應力,因此會利用應力頻率分布21來表現它,但當假定僅發生了某一大小的應力的情形下,等價應力振幅Seq例如能夠使用圖14中的式子而由應力頻率分布來計算。只要獲得等價應力振幅Seq,便能描繪等價應力振幅Seq中零件達破壞之壽命的機率密度函數16。故障風險之更新,能夠藉由此破壞壽命的機率密度函數之更新而實現。也就是說,能夠基於由故障資料庫15中包含之故障資料而計算出的概度、及更新前已利用之事前的壽命的密度函數,活用下式之貝氏定理,來獲得考量了故障資料之更新後的壽命的密度函數。 In the failure risk assessment/update unit 14, the fatigue phenomenon that occurs when the product is repeatedly subjected to the load is taken as a target, and when the failure probability F(t1) is obtained by the remaining life evaluation, the rain flow counting method or the like is used. The frequency analysis method is applied to the stress history occurring in the part from the past to the current time point t1, and the stress frequency distribution is obtained (see Fig. 4). In fact, the stress history contains various kinds of stresses, so the stress frequency distribution 21 is used to express it, but when it is assumed that only a certain magnitude of stress occurs, the equivalent stress amplitude Seq can be used, for example. The equation in 14 is calculated from the stress frequency distribution. As long as the equivalent stress amplitude Seq is obtained, the probability density function 16 of the life of the component in the equivalent stress amplitude Seq can be depicted. The update of the risk of failure can be achieved by updating the probability density function that destroys the lifetime. That is to say, based on the probability calculated by the fault data included in the fault database 15 and the density function of the pre-renewed life before the update, the Bayes' theorem of the following formula can be used to obtain the fault data. The density function of the updated lifetime.

更新後的密度函數=概度×事前的密度函數‧‧‧(8) Updated density function = probability × density function beforehand ‧ ‧ (8)

圖15揭示這樣的破斷壽命的密度函數之更新流程圖。首先,故障風險評估/更新部14,由對從對象零件p的過去至目前時間點t1為止之應力歷程做頻率分析而獲得之應力頻率分布,利用圖14中的式子計算等價應力振幅Seq(p)(S21)。然後,故障風險評估/更新部14,參照事先記憶於設計/材料資料庫5等之對象零件p的材料資料(P-S-N線圖),求出等價應力振幅Seq(p)下之破斷壽命的機率密度函數f(N)(S22)。此處,假設故障資料庫15中存在k個和對象零件p為相同零件,而裝載於同型機、類似機之零件pj的故障資料(j=1~k)。故障風險評估/更新部14,由從該些零件的運轉開始至故障時為止之環境資料、運轉資料、維護/運用資料及設計/材料資料,求出故障時為止之應力歷程及應力頻率分布,並遵照圖15中的右側的式子演算等價應力振幅Seq(pj)下之破斷壽命Nf j(S23)。故障風險評估/更新部14,能夠由如此求出的k個破斷壽命Nf j,遵照圖15中的左側的式子計算概度L(S24)。接著,故障風險評估/更新部14,能夠由事前的機率密度函數f(N)及概度L,遵照式(8)獲得更新後的破斷壽命的機率密度函數f(N)’(S25)。當處理破壞壽命的情形下,一般而言,事前的機率密度函數、概度的機率分布形狀會服從韋伯分布或對數常態分布。若參照更新後的破斷壽命的機率密度函數f(N)’,便能描繪更新後的P-S-N線圖。 Figure 15 reveals an updated flow chart of the density function for such breaking life. First, the failure risk assessment/update unit 14 calculates the stress frequency distribution obtained by performing frequency analysis on the stress history from the past to the current time point t1 of the target part p, and calculates the equivalent stress amplitude Seq using the equation in FIG. (p) (S21). Then, the failure risk assessment/update unit 14 refers to the material data (PSN line diagram) of the target part p previously stored in the design/material database 5, etc., and finds the breaking life under the equivalent stress amplitude Seq(p). Probability density function f(N) (S22). Here, it is assumed that there are k fault data in the fault database 15 and the target part p is the same part, and the fault data (j=1 to k) of the part p j mounted on the same machine or the like. The failure risk assessment/update unit 14 obtains the stress history and the stress frequency distribution from the time of the failure, the environmental data, the operation data, the maintenance/operation data, and the design/material data from the start of the operation of the components to the failure time. The breaking life N f j (S23) at the equivalent stress amplitude Seq(p j ) is calculated in accordance with the equation on the right side in Fig. 15 . The failure risk assessment/update unit 14 can calculate the probability L from the k breaking lifes N f j thus obtained in accordance with the equation on the left side in FIG. 15 (S24). Next, the failure risk assessment/update unit 14 can obtain the probability density function f(N)' of the updated breaking life in accordance with the equation (8) from the prior probability density function f(N) and the degree L (S25). . In the case of dealing with a destructive life, in general, the probability density function of the prior probability density and the probability distribution shape of the probation will obey the Weber distribution or the lognormal distribution. Referring to the updated probability density function f(N)' of the breaking life, the updated PSN line graph can be drawn.

接著,如圖4所示,能夠由將頻率分析法運用至從過去至目前時間點t1為止發生的應力歷程而得到之應力頻率分布、及事先訂定好之破壞機率P%的疲勞壽命曲線,而依下式計算相對於疲勞破壞P%而言之損傷度D(t1)。 Next, as shown in FIG. 4, it is possible to apply the frequency analysis method to the stress frequency distribution obtained from the stress history occurring from the past to the current time point t1, and the fatigue life curve obtained by setting the damage probability P% in advance. The damage degree D(t1) with respect to the fatigue damage P% is calculated according to the following formula.

D(t1)=(n1/N1)+(n2/N2)+…+(nm/Nm)‧‧‧(2) D(t1)=(n 1 /N 1 )+(n 2 /N 2 )+...+(n m /N m )‧‧‧(2)

此處,n1、n2、nm,分別為藉由應力歷程的頻率分析而得到之應力振幅S1、S2,Sm的反覆次數(m為整數)。此外,N1,N2,Nm,分別為當反覆負荷應力振幅S1、S2,Sm的情形下,在破壞機率P%下發生疲勞破壞之破斷反覆次數。此外,導致損傷度D(t1)發生之反覆次數N(t1)可依下式求出。 Here, n 1 , n 2 , and n m are the repetition times (m is an integer) of the stress amplitudes S 1 , S 2 , and S m obtained by frequency analysis of the stress history, respectively. Further, N 1 , N 2 , and N m are the number of times of breaking and the occurrence of fatigue fracture at the failure probability P% when the load stress amplitudes S 1 , S 2 , and S m are repeated. Further, the number of times of repetition N(t1) which causes the degree of damage D(t1) to occur can be obtained by the following equation.

N(t1)=D(t1)×Np‧‧‧(3) N(t1)=D(t1)×Np‧‧‧(3)

Np為事先訂定好之破壞機率P%的反覆次數,藉由將應力振幅訂為事先訂定好之應力振幅Si或Seq(p),如圖4中所示,能夠由N(t1)與疲勞壽命的機率密度函數,藉由圖4的式子,求出破壞機率F(t1)’。也就是說,若將從製品運轉開始至目前為止之經過時間訂為t1,則破壞機率F(t1)’為將密度函數f(N)’從0至N(t1)為止積分而得到的值。 Np is the number of times of the failure probability P% set in advance, and the stress amplitude is set to a predetermined stress amplitude Si or Seq(p), as shown in FIG. 4, which can be obtained by N(t1) and fatigue life. The probability density function, by the equation of Fig. 4, finds the probability of failure F(t1)'. In other words, if the elapsed time from the start of the product operation to the present is set to t1, the probability of destruction F(t1)' is a value obtained by integrating the density function f(N)' from 0 to N(t1). .

利用此更新後的P-S-N線圖來更新破壞機率,並將更新後的破壞機率F(t1)’與影響度C(p)相乘,藉此便能遵照式(1)演算更新後的故障風險RS(t1,p)’。像 這樣基於應力頻率分布、P-S-N線圖、故障資料之故障風險的更新方法係為一例,例如亦可不透過等價應力振幅而演算更新後的破壞機率。此外,亦可因應和視為對象之製品1為同型機、類似機13的使用環境、運轉狀況、構造的類似度等,來變化基於貝氏定理之破壞壽命的密度函數的更新式。 The updated PSN line graph is used to update the probability of destruction, and multiply the updated failure probability F(t1)' by the influence degree C(p), thereby being able to follow the updated (1) algorithm to calculate the risk of failure after updating. RS(t1,p)'. image Such an update method based on the stress frequency distribution, the P-S-N line graph, and the fault risk of the fault data is an example. For example, the probability of the damage after the update may be calculated without transmitting the equivalent stress amplitude. In addition, the update formula of the density function based on the damage life of the Bayes' theorem can be changed in response to the use of the product 1 as the object, the use environment of the similar machine 13, the operation state, the similarity of the structure, and the like.

此外,圖16為當依據製品中包含之零件的破壞壽命的機率密度函數來演算破壞機率的情形下,依據貝氏統計來更新事前設定之壽命的機率密度函數之例子示意圖。 In addition, FIG. 16 is a diagram showing an example of a probability density function for updating the life set in advance according to Bayesian statistics in the case where the probability of failure is calculated according to the probability density function of the breaking life of the components included in the product.

故障風險評估/更新部14中,當不參照應力頻率分布或P-S-N線圖,而是直接利用圖5的壽命的密度函數來演算破壞機率的情形下,也能利用式(8)所代表之貝氏定理。也就是說,能夠將壽命的密度函數與故障資料代入式(8)來更新壽命的密度函數,並藉由此更新後的機率密度函數來求出破壞機率F(t1)’(圖16),並將其乘上影響度來更新故障風險。但,若和製品1為同型機、類似機13而當使用環境或運轉狀況大幅相異的情形下,可以想見即使運作時間相同,和製品1的零件為同型機、類似機13之零件其破壞壽命仍會大幅相異。這樣的情形下,除時間外,還會由環境資料、運轉資料選擇幾個對破壞壽命造成影響之物理量來作為變數,而作成多變量的壽命的機率密度函數。 In the failure risk assessment/update unit 14, when the probability of failure is directly calculated by using the density function of the life of FIG. 5 without referring to the stress frequency distribution or the PSN line diagram, the shell represented by the equation (8) can also be utilized. Theorem. That is to say, the density function of the lifetime and the fault data can be substituted into the equation (8) to update the density function of the lifetime, and the probability of failure F(t1)' can be obtained by the updated probability density function (Fig. 16). And multiply it by the impact to update the risk of failure. However, in the case where the product 1 is the same type machine or the like machine 13 and the use environment or the operation state is greatly different, it is conceivable that even if the operation time is the same, the parts of the product 1 are the same type machine, the parts of the similar machine 13 Destructive life will still vary significantly. In such a case, in addition to the time, several physical quantities affecting the damage life are selected as the variables from the environmental data and the operational data, and the probability density function of the multivariate life is created.

另,壽命的機率密度函數f(t1)之資料,例如能夠對每一對象製品事先記憶於故障風險評估部2的內部記憶 體或設計/材料資料庫5等適當的記憶體。此外,壽命的機率密度函數f(t1)之資料,能夠使用環境資料、運轉資料等而事先訂定。 In addition, the data of the probability density function f(t1) of the lifetime can be stored, for example, in the internal memory of the failure risk assessment unit 2 for each object product. Appropriate memory such as body or design/material database. In addition, the data of the probability density function f(t1) of the life can be determined in advance using environmental data, operation data, and the like.

圖17為當以包含時間之多變量的機率密度函數來表現破壞壽命,而演算破壞機率的情形下,依據貝氏統計來更新事前設定之機率密度函數之例子示意圖。 Fig. 17 is a diagram showing an example of updating the probability density function set in advance according to Bayesian statistics in the case where the failure life is expressed by the probability density function including the variable of time.

圖17為2變量的壽命的機率密度函數之例子,本例中是將把從環境資料或運轉資料選擇出的物理量x予以變換而成之值x’,和時間一起訂為機率變數。此變換中,包含對於隨時間變化的物理量的統計量(某一時間間隔下之平均值、最大值等)之變換、或對於依據頻率分析的等價物理量之變換等。即使壽命的密度函數為多變量,如圖17所示,如同1變量的情形般,仍能遵照貝氏定理來更新壽命的密度函數。此外,亦能作成包含對破壞壽命造成大幅影響的物理量或把其變換而成的物量x’以及時間t之函數r(t、x’),並將該函數解釋成機率變數來作成破壞壽命的機率密度函數。此外,故障風險評估/更新部14中,當藉由徵兆診斷來演算破壞機率的情形下,同樣地利用式(8)所代表之貝氏定理,來更新破壞的機率密度函數即可。由此更新後的機率密度求出破壞機率,並將其乘上影響度來更新故障風險即可。 Fig. 17 is an example of the probability density function of the life of the two variables. In this example, the value x' obtained by converting the physical quantity x selected from the environmental data or the operation data is set as the probability variable together with the time. This transformation includes a statistic for a physical quantity that changes with time (a mean value, a maximum value, etc. at a certain time interval), or a transformation of an equivalent physical quantity according to frequency analysis, and the like. Even if the density function of the lifetime is multivariate, as shown in Fig. 17, as in the case of the 1 variable, the density function of the lifetime can be updated in accordance with Bayes' theorem. In addition, it is also possible to create a function r(t, x') which includes a physical quantity which greatly affects the damage life or converts it into a quantity x' and a time t, and interprets the function as a probability variable to create a failure life. Probability density function. Further, in the failure risk assessment/update unit 14, when the probability of destruction is calculated by the symptom diagnosis, the probability density function of the failure can be updated by using the Bayes' theorem represented by the equation (8) in the same manner. The probability of damage is determined by the updated probability density, and the probability of failure is multiplied to increase the risk of failure.

另,2變量的壽命的機率密度函數,例如能夠對每一對象製品事先記憶於故障風險評估部2的內部記憶體或設計/材料資料庫5等適當的記憶體。此外,2變量的壽命 的機率密度函數,能夠使用環境資料、運轉資料等而事先訂定。 Further, the probability density function of the life of the two variables can be stored, for example, in an appropriate memory such as the internal memory of the failure risk assessment unit 2 or the design/material database 5 for each target product. In addition, the life of 2 variables The probability density function can be determined in advance using environmental data, operational data, and the like.

其後,故障風險評估/更新部14及維護/運用腳本制定部3,運用更新後的故障風險RS(t1,p),如同圖12所示之第1實施形態般,執行故障風險趨勢分析(S13)、將來風險預測(S14)、維護/運用變更(S15)、維護/運用腳本制定(S16)之各處理。 Then, the failure risk assessment/update unit 14 and the maintenance/operation scenario creation unit 3 perform the failure risk trend analysis using the updated failure risk RS(t1, p) as in the first embodiment shown in FIG. S13), future risk prediction (S14), maintenance/operation change (S15), and maintenance/application script development (S16).

藉由本實施形態中利用了製品1及同型機、類似機13的故障資料來做故障風險之更新,能夠使破壞機率的演算高精度化,故可更高精度地實施製品1中包含之零件的故障風險之評估、預測。故障風險之更新的時間間隔,可和故障風險的評估間隔相同,亦可相異。藉由增多故障風險之更新的頻率,便能更高精度評估製品中包含之零件的故障風險。 By using the failure data of the product 1 and the similar machine and the similar machine 13 in the present embodiment to update the failure risk, the calculation of the destruction probability can be made highly accurate, so that the components included in the product 1 can be implemented with higher precision. Assessment and prediction of failure risk. The time interval for updating the risk of failure may be the same as the interval for assessing the risk of failure, or it may be different. By increasing the frequency of the update of the risk of failure, the risk of failure of the parts contained in the product can be assessed with greater precision.

〔第3實施形態〕 [Third embodiment]

接著,利用圖18、19,說明依本實施形態第3實施形態之運作輔助系統。 Next, an operation assisting system according to a third embodiment of the present embodiment will be described with reference to Figs.

圖18為依本發明第3實施形態之運作輔助系統的主要構成要素、及提供運作輔助系統中利用之資料的製品及資料庫、及它們的關係概略示意方塊圖。 Fig. 18 is a schematic block diagram showing the main components of the operation support system according to the third embodiment of the present invention, and the products and databases for providing the data used in the operation support system, and their relationships.

本實施形態之維護/運用腳本制定部3中,除了藉由組裝於製品1之感測器而計測的環境資料、運轉資料,或針對製品1的設計資料、材料資料、維護/運用資料外, 還利用和製品1無關之來自外部資料庫17的資訊。所謂外部資料庫17中包含之外部資料,例如包含藉由大型計算機計算出的氣象、海象的將來預測資料、或資源的供給預測資料、資源的埋藏預測資料等。這樣的外部資料,不會影響製品1的運轉狀態,是故外部資料並非和製品1有關之物。 In the maintenance/operation scenario creation unit 3 of the present embodiment, in addition to the environmental data and the operation data measured by the sensor incorporated in the product 1, or the design data, the material data, and the maintenance/application data for the product 1, Information from the external database 17 that is not related to the article 1 is also utilized. The external data included in the external database 17 includes, for example, weather forecasting data of a seahorse calculated by a large computer, prediction data of resources, and burial prediction data of resources. Such external data does not affect the operation state of the product 1, and therefore the external data is not related to the product 1.

圖19為依本發明第3實施形態之運作輔助系統的維護/運用腳本制定部的主要構成要素,與在要素間往來之資料的流向概略示意方塊圖。 Fig. 19 is a schematic block diagram showing the main components of the maintenance/operation scenario creation unit of the operation support system according to the third embodiment of the present invention, and the flow of information between the elements.

如圖19般,維護/運用腳本制定部3的故障風險預測11,是利用故障風險的趨勢曲線、維護/運用腳本及外部資料,演算將來的故障風險。也就是說,能夠由環境資料或運轉資料中包含之物理量x、維護/運用資料中包含之物理量y、外部資料中包含之物理量z,依下式演算故障風險的預測值。 As shown in Fig. 19, the failure risk prediction 11 of the maintenance/operation scenario creation unit 3 calculates the risk of future failure by using the trend curve of the failure risk, the maintenance/operation script, and the external data. That is to say, the physical quantity x contained in the environmental data or the operational data, the physical quantity y included in the maintenance/application data, and the physical quantity z included in the external data can be calculated according to the following formula.

RS(t+△T,y,z,p)=g5(t+△T,y,z,p)‧‧‧(9) RS(t+△T, y, z, p)=g5(t+△T, y, z, p)‧‧‧(9)

RS(t+△T,x,y,z,p)=g6(t+△T,x,y,z,p)‧‧‧(10) RS(t+△T,x,y,z,p)=g6(t+△T,x,y,z,p)‧‧‧(10)

藉由外部資料之利用,能夠更廣泛地考量製品1周遭之環境,因此能夠更精度良好地預測將來的故障風險。 By utilizing the use of external data, it is possible to more widely consider the environment in which the product is exposed for one week, and thus it is possible to predict the risk of future failure more accurately.

〔第4實施形態〕 [Fourth embodiment]

接著,利用圖20、21,說明依本實施形態第4實施形態之運作輔助系統。 Next, an operation assisting system according to a fourth embodiment of the present embodiment will be described with reference to Figs.

圖20為依本發明第4實施形態之運作輔助系統的主 要構成要素、及提供運作輔助系統中利用之資料的製品及資料庫、及它們的關係概略示意方塊圖。 Figure 20 is a diagram showing the main function of the operation assisting system according to the fourth embodiment of the present invention. A schematic block diagram of the components and the products and databases that provide the elements used in the operational support system and their relationships.

此外,圖21為依本發明第4實施形態之運作輔助系統的維護/運用腳本制定部的主要構成要素,與在要素間往來之資料的流向概略示意方塊圖。 In addition, FIG. 21 is a schematic block diagram showing the main components of the maintenance/operation scenario creation unit of the operation support system according to the fourth embodiment of the present invention, and the flow of information between the elements.

本實施形態中將在製品1計測出的環境資料、運轉資料輸入至破壞機率評估部18,並在該處演算製品1中包含之零件的破壞機率F。演算出的破壞機率F,會成為維護/運用腳本制定部3的輸入,並在該處將破壞機率與影響度相乘而作成故障風險的趨勢曲線。 In the present embodiment, the environmental data and the operation data measured by the product 1 are input to the destruction probability evaluation unit 18, and the destruction probability F of the components included in the product 1 is calculated therein. The calculated failure probability F becomes an input to the maintenance/operation scenario creation unit 3, where the probability of failure is multiplied by the degree of influence to create a trend curve of the risk of failure.

圖22為當作成破壞機率的趨勢曲線的情形下之維護/運用腳本制定部的主要構成要素,與在要素間往來之資料的流向概略示意方塊圖。 Fig. 22 is a schematic block diagram showing the flow of the main components of the maintenance/operation scenario creation unit in the case of the tendency curve of the destruction probability and the flow of information between the elements.

是故,本實施形態中如圖21般,影響度資料庫8配置於維護/運用腳本制定部3。另一方面,如圖22所示,亦可設想不作成故障風險的趨勢曲線,而是作成破壞機率的趨勢曲線之形態。這樣的情形下,能夠將零件p視為對象,由時間t與對破壞機率或破壞機率造成影響之物理量x,作成如下式般之破壞機率F的趨勢曲線。 Therefore, in the present embodiment, as shown in Fig. 21, the influence degree database 8 is placed in the maintenance/operation scenario creation unit 3. On the other hand, as shown in FIG. 22, it is also conceivable to form a trend curve which does not cause a risk of failure, but to form a trend curve which destroys the probability. In such a case, the part p can be regarded as an object, and the time t and the physical quantity x which affects the probability of destruction or the probability of destruction are formed as a trend curve of the destruction probability F as follows.

F(t,p)=h1(t,p)‧‧‧(11) F(t,p)=h1(t,p)‧‧‧(11)

F(t,x,p)=h2(t,x,p)‧‧‧(12) F(t,x,p)=h2(t,x,p)‧‧‧(12)

式(11)為僅時間t會對破壞機率造成影響之情形,式(12)為時間t與計測出的其他物理量x會對破壞機率或破壞機率造成影響之情形。故障風險預測11中,依下式計算當更改維護/運用的情形下之將來的破壞機率F。 Equation (11) is a case where only the time t affects the probability of destruction, and the equation (12) is a case where the time t and the measured other physical quantity x have an influence on the probability of destruction or the probability of destruction. In the failure risk prediction 11, the future failure probability F in the case of changing maintenance/utilization is calculated according to the following formula.

F(t+△T,y,p)=h3(t+△T,y,p)‧‧‧(13) F(t+△T,y,p)=h3(t+△T,y,p)‧‧‧(13)

F(t+△T,x,y,p)=h4(t+△T,x,y,p)‧‧‧(14) F(t+△T,x,y,p)=h4(t+△T,x,y,p)‧‧‧(14)

此處,y為維護/運轉資料當中,對故障風險RS造成影響之物理量。將此將來的破壞機率或破壞機率乘上影響度C(p),便能如下式般決定將來的故障風險RS。 Here, y is the physical quantity that affects the failure risk RS in the maintenance/operation data. By multiplying the future probability of destruction or the probability of failure by the influence degree C(p), the future failure risk RS can be determined as follows.

RS(t+△T,y,p)=C(p)×h3(t+△T,y,p)‧‧‧(15) RS(t+△T,y,p)=C(p)×h3(t+△T,y,p)‧‧‧(15)

RS(t+△T,x,y,p)=C(p)×h4(t+△T,x,y,p)‧‧‧(16) RS(t+△T,x,y,p)=C(p)×h4(t+△T,x,y,p)‧‧‧(16)

〔第5實施形態〕 [Fifth Embodiment]

圖23為依本發明第5實施形態之運作輔助系統的主要構成要素、及提供運作輔助系統中利用之資料的製品及資料庫、及它們的關係概略示意方塊圖。 Fig. 23 is a schematic block diagram showing the main components of the operation support system according to the fifth embodiment of the present invention, and the products and databases for providing the data used in the operation support system, and their relationships.

以下,利用圖23,說明依本實施形態第5實施形態之運作輔助系統。本實施形態中將在製品1計測出的環境資料、運轉資料輸入至損傷度評估部19,並在該處演算製品1中包含之零件的損傷度。損傷度如圖4所示,當由應力歷程而以剩餘壽命評估求出破壞機率的情形下,圖3的P-S-N線圖中,能夠選擇某一破壞機率的疲勞壽命曲線來求出。演算出的損傷度,會成為維護/運用腳本制定部3的輸入,並在該處參照材料資料來計算破壞機率,並將其乘上影響度而作成故障風險的趨勢曲線。是故,本實施形態中如圖21般,影響度資料庫8配置於維護/運用腳本制定部3。另一方面,此處亦如圖22所示,亦可設想不作成故障風險的趨勢曲線,而是作成破壞機率的趨勢曲 線之形態。 Hereinafter, an operation assisting system according to a fifth embodiment of the present embodiment will be described with reference to Fig. 23 . In the present embodiment, the environmental data and the operation data measured by the product 1 are input to the damage degree evaluation unit 19, and the damage degree of the components included in the product 1 is calculated there. As shown in FIG. 4, when the damage probability is obtained from the stress history and the remaining life evaluation, the P-S-N diagram of FIG. 3 can be obtained by selecting a fatigue life curve of a certain probability of destruction. The calculated degree of damage is input to the maintenance/application scripting unit 3, and the material probability is calculated by referring to the material data, and multiplied by the influence degree to create a trend curve of the risk of failure. Therefore, in the present embodiment, as shown in Fig. 21, the influence degree database 8 is placed in the maintenance/operation scenario creation unit 3. On the other hand, as shown in Fig. 22, it is also possible to assume a trend curve that does not pose a risk of failure, but to create a trend of destruction probability. The shape of the line.

又,圖24為當作成損傷度的趨勢曲線的情形下之維護/運用腳本制定部的主要構成要素,與在要素間往來之資料的流向概略示意方塊圖。如圖24所示,亦可設想作成損傷度的趨勢曲線之形態。也就是說,亦能作成將來的損傷度d3、d4,與影響度C(p)相乘,而如下式般決定故障風險的趨勢曲線。 In addition, FIG. 24 is a schematic block diagram showing the main components of the maintenance/operation scenario creation unit in the case of the trend curve of the damage degree, and the flow of the data between the elements. As shown in Fig. 24, it is also conceivable to form a trend curve of the degree of damage. In other words, the future damage degree d3, d4 can be multiplied by the influence degree C(p), and the trend curve of the failure risk is determined as follows.

F(t+△T,y,p)=K(p)×d3(t+△T,y,p)‧‧‧(17) F(t+△T,y,p)=K(p)×d3(t+△T,y,p)‧‧‧(17)

F(t+△T,x,y,p)=K(p)×d4(t+△T,x,y,p)‧‧‧(18) F(t+△T,x,y,p)=K(p)×d4(t+△T,x,y,p)‧‧‧(18)

RS(t+△T,y,p)=C(p)×K×d3(t+△T,y,p)‧‧‧(19) RS(t+△T,y,p)=C(p)×K×d3(t+△T,y,p)‧‧‧(19)

RS(t+△T,x,y,p)=C(p)×K×d4(t+△T,x,y,p)‧‧‧(20) RS(t+△T,x,y,p)=C(p)×K×d4(t+△T,x,y,p)‧‧‧(20)

此處K(p)是為了由損傷度計算破壞機率所必須之換算常數。K(p)亦可蓄積於損傷度資料庫、影響度資料庫或是其他資料庫。本實施形態中,故障機率或損傷度是暫且獨立地被演算,因此有故障機率或損傷度的可視化(以顯示部顯示)變得容易之優點。 Here K(p) is the conversion constant necessary to calculate the probability of damage from the damage degree. K(p) can also be accumulated in the damage database, impact database or other database. In the present embodiment, since the failure probability or the damage degree is temporarily calculated independently, there is an advantage that the failure probability or the damage degree is visualized (displayed by the display unit).

遵照式(1)演算的故障風險,為將從製品的運轉開始時至時刻t1為止之期間中,於製品累積的損傷予以全部納入考量而得的故障風險。另一方面,亦能計算僅將從目前時間點t1至將來的某一時間點t1+△T為止之期間中,於製品累積的損傷予以納入考量而得的故障風險,並以其作為故障風險的預測值。這樣的情形下,首先依下式計算從目前時間點t1至將來的某一時間點t1+△T之期間 中,製品中包含之某一零件損壞的機率P。 The risk of failure according to the formula (1) is a risk of failure that is included in the damage accumulated in the product from the start of the operation of the product to the time t1. On the other hand, it is also possible to calculate the risk of failure due to the damage accumulated in the product from the current time point t1 to a certain time point t1 + ΔT in the future, and use it as a risk of failure. Predictive value. In such a case, first calculate the period from the current time point t1 to a certain time point t1 + ΔT in the future according to the following equation Medium, the probability of damage to a part contained in the product.

P(t1,t1+△T)=(F(t1+△T)-F(t1))/(1-F(t1)) P(t1, t1+△T)=(F(t1+△T)-F(t1))/(1-F(t1))

‧‧‧(21) ‧‧‧(twenty one)

此處,F(t1)、F(t1+△T)分別為時刻t1及t1+△T下之破壞機率。然後,從目前時間點t1至將來的某一時間點t1+△T為止之期間的故障風險,能夠由對象零件的影響度C(p)、破壞機率P(t1,t1+△T)而依下式計算。 Here, F(t1) and F(t1+ΔT) are the probability of destruction at time t1 and t1+ΔT, respectively. Then, the risk of failure from the current time point t1 to a certain time point t1 + ΔT in the future can be determined by the influence degree C(p) of the target component and the probability of destruction P(t1, t1 + ΔT). Calculation.

RS(t1+△T)=C(p)×P(t1,t1+△T)‧‧‧(22) RS(t1+△T)=C(p)×P(t1,t1+△T)‧‧‧(22)

式(21)、(22)的破壞機率、故障風險,會因應環境資料、運轉資料、維護/運用資料而變動,因此破壞機率、故障風險並非僅是時間的函數,但式(21)、(22)中為求簡單而記述成時間的函數。 The failure rate and failure risk of equations (21) and (22) will change according to environmental data, operation data, and maintenance/application data. Therefore, the probability of failure and the risk of failure are not only a function of time, but equation (21), 22) A function described as time for simplicity.

C.附記 C. Notes

另,本發明並非限定於上述實施例,還包含各式各樣的變形例。例如,上述實施例是為了以易於理解的方式說明本發明而詳細說明,未必限定於具備所說明的全部構成。此外,可將某一實施例的構成的一部分置換成其他實施例的構成,此外,亦可在某一實施例的構成中加入其他實施例的構成。此外,針對各實施例的構成的一部分,可追加/刪除/置換其他構成。 Further, the present invention is not limited to the above embodiment, and various modifications are also included. For example, the above-described embodiments are described in detail to explain the present invention in an easy-to-understand manner, and are not necessarily limited to having all of the configurations described. Further, a part of the configuration of a certain embodiment may be replaced with a configuration of another embodiment, and a configuration of another embodiment may be added to the configuration of a certain embodiment. Further, other components may be added, deleted, or replaced for a part of the configuration of each embodiment.

此外,上述各構成、機能、處理部、處理手段等,亦可將它們的一部分或全部例如藉由以積體電路來設計等而 以硬體實現。此外,上述各構成、機能等,亦可藉由解譯、執行以處理器實現各個機能之程式而以軟體實現。實現各機能之程式、表格、檔案等資訊,能夠置放於記憶體、或硬碟、SSD(Solid State Drive)等記錄裝置,或IC卡、SD卡、DVD等記錄媒體。 Further, each of the above-described configurations, functions, processing units, processing means, and the like may be designed by an integrated circuit or the like, for example, by a part or all of them. Implemented in hardware. Further, each of the above-described configurations, functions, and the like can be realized by software by interpreting and executing a program that realizes each function by a processor. Information such as programs, tables, and files for each function can be placed in a memory, or a recording device such as a hard disk or an SSD (Solid State Drive), or a recording medium such as an IC card, an SD card, or a DVD.

此外,控制線或資訊線係揭示說明上認有必要者,就製品上未必揭示了全部的控制線或資訊線。實際上亦可想成幾乎全部的構成均相互連接。 In addition, if the control line or information line reveals that it is necessary, it may not reveal all the control lines or information lines on the product. In fact, it is also conceivable that almost all of the components are connected to each other.

1‧‧‧對象製品 1‧‧‧ object products

2‧‧‧故障風險評估部 2‧‧‧Fault Risk Assessment Department

3‧‧‧維護/運用腳本制定部 3‧‧‧Maintenance/Using Script Development Department

4‧‧‧故障風險預測部 4‧‧‧Fault Risk Forecasting Department

5‧‧‧設計/材料資料庫 5‧‧‧Design/Materials Database

100‧‧‧運作輔助系統 100‧‧‧Operational support system

Claims (15)

一種運作輔助裝置,具備:故障風險評估部;及維護/運用腳本制定部;前述故障風險評估部,利用從對象製品的複數個感測器輸入之環境資料及運轉資料、與事先訂定好之設計資料及材料資料,來演算視為對象之零件p的時刻t1下之破壞機率F(t1),於時刻t1,藉由零件p的時刻t1下的破壞機率F(t1)與當零件p損壞了的情形下之事先訂定好之每一零件p的影響度C(p)之乘積,來計算製品中包含之零件p的故障風險RS(t1,p),前述維護/運用腳本制定部,依據從前述故障風險評估部送來的時刻t1下之零件p的故障風險RS(t1、p)、及已從前述故障風險評估部送來而記憶於故障風險資料庫之從過去至時刻t1為止的複數個故障風險,而將從事先由對象製品輸入的環境資料及運轉資料選擇出之對故障風險造成影響的物理量x、以及時間t訂為變數,生成故障風險的趨勢曲線,依據故障風險的趨勢曲線,求出從目前時間點至往前之事先訂定好之時間的故障風險的預測值,依據故障風險的預測值,遵照依維護/運用者所做的來自輸入部之輸入而手動地設定出之、或藉由事先訂定好 之處理而自動地設定出之製品的維護/運用腳本,而生成以時刻t、物理量x、以及由維護資料及/或運轉資料選擇出之對故障風險造成影響的物理量y作為變數之故障風險預測模型,來預測零件p的將來的故障風險,將每一零件的將來的故障風險,依故障風險的預測值較高之順序予以整理並顯示於顯示部及/或記憶於記憶部,並以事先訂定好之複數個閾值予以劃分群組,藉由依維護/運用者所做的來自輸入部之輸入而手動地、或藉由事先訂定好之處理而自動地對每一群組設定包含各零件的維護時期在內之維護運用內容,並將維護運用內容記憶於記憶部及/或顯示於顯示部。 An operation assisting device comprising: a failure risk assessment unit; and a maintenance/application script development unit; the failure risk assessment unit uses environmental data and operation data input from a plurality of sensors of the target product, and a predetermined design The data and material data are used to calculate the probability of failure F(t1) at time t1 of the part p regarded as the object. At time t1, the probability of failure F(t1) at time t1 of part p is damaged when part p is In the case of the predetermined product C (p) of each part p, the risk of failure of the part p contained in the product is calculated (t1, p), and the maintenance/application scripting unit is based on The failure risk RS (t1, p) of the part p at the time t1 sent from the failure risk assessment unit, and the memory that has been sent from the failure risk assessment unit and memorized in the failure risk database from the past to the time t1 A plurality of failure risks, and the physical quantity x and the time t which are selected from the environmental data and the operation data input in advance by the target product are selected as variables, and a trend curve of the failure risk is generated, The trend curve of the risk, and the predicted value of the risk of failure from the current time point to the predetermined time in advance is obtained, and the predicted value of the risk is determined according to the input from the input unit according to the maintenance/operator. Set it up, or set it in advance The process of automatically setting and maintaining the maintenance/application script of the product, and generating a fault risk prediction with the physical quantity y that affects the failure risk selected by the maintenance data and/or the operation data as the variable at time t, the physical quantity x, and the maintenance data and/or the operation data Model to predict the future failure risk of part p, and to sort out the future failure risk of each part according to the higher predicted value of the failure risk and display it on the display part and/or in the memory part, and A plurality of thresholds are set in advance to divide the group, and each part is automatically set for each group by manual input by the maintenance/operator from the input unit or by a predetermined process. The maintenance content is maintained during the maintenance period, and the maintenance application content is memorized in the memory unit and/or displayed on the display unit. 如申請專利範圍第1項所述之運作輔助裝置,其中,前述故障風險評估部,利用環境資料、運轉資料及設計資料,計算從過去至目前時間點t1為止在零件發生的應力歷程,對應力歷程,運用雨流計數法或其他頻率分析法,作成應力頻率分布,由應力頻率分布、及事先訂定好之疲勞壽命曲線,計算損傷度D(t1),求出導致損傷度D(t1)發生之反覆次數N(t1),將作為材料資料之事先訂定好之疲勞壽命的機率密度函數從0至N(t1)為止積分,求出破壞機率F(t1)。 The operation assisting device according to claim 1, wherein the fault risk assessment unit calculates the stress history of the component from the past to the current time point t1 using the environmental data, the operation data, and the design data. In the course of the process, the rain flow counting method or other frequency analysis method is used to create the stress frequency distribution, and the stress frequency distribution and the fatigue life curve determined in advance are calculated, and the damage degree D(t1) is calculated, and the damage degree D(t1) is determined. The number of times of repetition N (t1) is integrated from 0 to N (t1) as a function of the probability density of the fatigue life defined in advance as the material data, and the probability of failure F(t1) is obtained. 如申請專利範圍第1或2項所述之運作輔助裝 置,其中,前述維護/運用腳本制定部,具備:記憶體,事先記憶有將製品予以分類而成的各組件中包含之一或複數個零件,參照前述記憶體,再將對於組件之各零件的將來的故障風險的預測值予以整理並顯示於顯示部及/或記憶於記憶部。 Operation aids as described in claim 1 or 2 In the above-described maintenance/application scenario development unit, the memory includes a memory, and one or a plurality of components included in each component for classifying the product are stored in advance, and the components of the component are referred to by referring to the memory. The predicted value of the future failure risk is sorted and displayed on the display unit and/or stored in the memory unit. 如申請專利範圍第1項所述之運作輔助裝置,其中,具備:外部資料庫,事先記憶有和製品無關之外部資料,前述維護/運用腳本制定部,利用故障風險的趨勢曲線、及維護/運用腳本、及外部資料,演算將來的故障風險。 The operation support device according to the first aspect of the invention, comprising: an external database, in which external data irrelevant to the product is stored in advance, and the maintenance/application script development unit uses a trend curve of failure risk and maintenance/ Use scripts and external data to calculate future failure risks. 如申請專利範圍第1或2項所述之運作輔助裝置,其中,前述維護/運用腳本制定部,作成破壞機率F的趨勢曲線、及/或損傷度的趨勢曲線。 The operation assisting device according to claim 1 or 2, wherein the maintenance/operation scenario creating unit creates a trend curve for destroying the probability F and/or a trend curve for the damage degree. 一種風力發電系統,具備:如申請專利範圍第1或2項所述之運作輔助裝置;及風力發電機,為具有複數個感測器之對象製品。 A wind power generation system comprising: the operation auxiliary device according to claim 1 or 2; and the wind power generator, which is an object product having a plurality of sensors. 一種運作輔助裝置,具備:故障風險評估/更新部;及 維護/運用腳本制定部;故障資料庫,蓄積構成對象製品與其同型機及/或類似機之複數個零件的故障資料;前述故障風險評估/更新部,基於對象零件p達破壞之壽命的機率密度函數、及由前述故障資料庫中包含之故障資料計算出的概度,活用貝氏定理,求出考量了故障資料之更新後的壽命的機率密度函數,藉由更新後的壽命的機率密度函數,利用從對象製品的複數個感測器輸入之環境資料及運轉資料、與事先訂定好之設計資料及材料資料,來演算視為對象之零件p的時刻t1下之更新後的破壞機率F(t1)’,於時刻t1,藉由零件p的時刻t1下的更新後的破壞機率F(t1)’與當零件p損壞了的情形下之事先訂定好之每一零件p的影響度C(p)之乘積,來計算製品中包含之零件p的更新後的故障風險RS(t1,p)’,前述維護/運用腳本制定部,依據從前述故障風險評估/更新部送來的時刻t1下之零件p的更新後的故障風險RS(t1、p)’、及已從前述故障風險評估/更新部送來而記憶於故障風險資料庫之從過去至時刻t1為止的複數個更新後的故障風險,而將事先由對象製品輸入的環境資料及運轉資料選擇出之對故障風險造成影響的物理量x、以及時間t訂為變數,生成故障風險的趨勢曲線, 依據故障風險的趨勢曲線,求出從目前時間點至往前了事先訂定好之時間的故障風險的預測值,依據故障風險的預測值,遵照依維護/運用者所做的來自輸入部之輸入而手動地設定出之、或藉由事先訂定好之處理而自動地設定出之製品的維護/運用腳本,而生成以時刻t、物理量x、以及由維護資料及/或運轉資料選擇出之對故障風險造成影響的物理量y作為變數之故障風險預測模型,來預測零件p的將來的故障風險,並將預測和趨勢曲線一起記憶於記憶部及/或顯示於顯示部。 An operation auxiliary device having: a failure risk assessment/update unit; and Maintenance/application scripting unit; fault database, accumulating fault data of a plurality of parts constituting the object product and its same type machine and/or the like; the foregoing fault risk assessment/update unit, based on the probability density of the life of the target part p The function and the probability calculated from the fault data contained in the fault database described above, using the Bayes' theorem, the probability density function that takes into account the updated lifetime of the fault data, and the probability density function of the updated lifetime Using the environmental data and operation data input from a plurality of sensors of the target product, and the design data and material materials set in advance, the updated probability of destruction F at the time t1 of the part p regarded as the object is calculated ( T1)', at time t1, the degree of influence of each part p which is determined in advance in the case where the part p is damaged by the updated failure probability F(t1)' at the time t1 of the part p (p) the product to calculate the updated failure risk RS(t1, p)' of the part p included in the product, and the maintenance/application scripting unit is sent from the failure risk assessment/update unit The updated failure risk RS(t1, p)' of the part p at the time t1, and a plurality of pieces from the past to the time t1 that have been received from the failure risk assessment/update unit and memorized in the failure risk database The updated risk of failure, and the physical quantity x and the time t that affect the risk of failure selected by the environmental data and the operational data input from the target product are set as variables, and a trend curve of the risk of failure is generated. According to the trend curve of the fault risk, the predicted value of the fault risk from the current time point to the time set in advance is obtained, and the input from the input unit according to the maintenance/operator is followed according to the predicted value of the fault risk. And the maintenance/application script of the product that is manually set or automatically set by the predetermined processing to generate the pair selected by the time t, the physical quantity x, and the maintenance data and/or the operation data The physical quantity y affected by the failure risk is used as a variable risk prediction model to predict the future failure risk of the part p, and the prediction and trend curves are memorized together in the memory and/or displayed on the display. 如申請專利範圍第7項所述之運作輔助裝置,其中,前述故障風險評估/更新部,由對從零件p的過去至目前時間點t1為止之應力歷程做頻率分析而獲得之應力頻率分布,利用以下式子計算等價應力振幅Seq(p), n i :應力振幅S i 的頻率m:疲勞壽命曲線的斜率參照零件p的材料資料,求出等價應力振幅Seq(p)下之破斷壽命的機率密度函數f(N),若前述故障資料庫中,存在k個和對象零件p為相同 零件,而裝載於同型機、類似機之零件pj的故障資料(j=1~k)、則由從該些零件的運轉開始至故障時為止之環境資料、運轉資料、維護/運用資料及設計/材料資料,求出故障時為止之應力歷程及應力頻率分布,並遵照下式演算等價應力振幅Seq(p)下之破斷壽命 n i :應力振幅S i 的頻率m:疲勞壽命曲線的斜率由求出的k個破斷壽命,遵照下式計算概度L, 由事前的機率密度函數f(N)及概度L,遵照下式獲得更新後的破斷壽命的機率密度函數f(N)’,更新後的機率密度函數f(N)’=概度×事前的機率密度函數f(N)藉由更新後的破斷壽命的機率密度函數f(N)’計算更新後的破壞機率F(t1)’,將更新後的破壞機率F(t1)’與事先訂定好之每一零件p的影響度C(p)相乘,藉此演算更新後的故障風險RS(t1、p)’。 The operation assisting device according to claim 7, wherein the failure risk assessment/update unit obtains a stress frequency distribution obtained by performing frequency analysis on a stress history from a past part to a current time point t1 of the part p, Calculate the equivalent stress amplitude Seq(p) using the following equation, n i : frequency m of stress amplitude S i : slope of fatigue life curve Referring to the material data of part p, the probability density function f(N) of the breaking life under the equivalent stress amplitude Seq(p) is obtained, if the above failure In the database, there are k and the object parts p are the same parts, and the fault data (j=1~k) of the parts p j loaded on the same machine or the like, from the start of the operation of the parts to the failure The environmental data, operation data, maintenance/application data, and design/material data so far, and the stress history and stress frequency distribution up to the fault are obtained, and the breaking life under the equivalent stress amplitude Seq(p) is calculated according to the following formula , n i : frequency m of the stress amplitude S i : slope of the fatigue life curve from the obtained k breaking lives , calculate the probability L according to the following formula, From the probability density function f(N) and the probability L beforehand, the probability density function f(N)' of the updated breaking life is obtained according to the following equation, and the updated probability density function f(N)'=probability× The ex ante probability density function f(N) calculates the updated failure probability F(t1)' by the updated probability density function f(N)' of the breaking life, and the updated failure probability F(t1)' The influence degree C(p) of each part p is multiplied in advance, thereby calculating the updated failure risk RS(t1, p)'. 如申請專利範圍第7或8項所述之運作輔助裝置,其中,前述故障風險評估/更新部,利用環境資料、運轉資料及設計資料,計算從過去至目前時間點t1為止在零件發生的應力歷程,對應力歷程,運用雨流計數法或其他頻率分析法,作成應力頻率分布,由應力頻率分布、及事先訂定好之破壞機率P%的疲勞壽命曲線,計算相對於疲勞破壞P%而言之損傷度D(t1),求出導致損傷度D(t1)發生之反覆次數N(t1),將獲得之更新後的疲勞壽命的機率密度函數f(N)’從0至N(t1)為止積分,求出更新後的破壞機率F(t1)’。 The operation assisting device according to claim 7 or 8, wherein the fault risk assessment/update unit uses the environmental data, the operation data, and the design data to calculate the stress occurring in the part from the past to the current time point t1. History, the stress history, using the rain flow counting method or other frequency analysis method, the stress frequency distribution, the stress frequency distribution, and the fatigue life curve of the pre-determined failure probability P%, calculated relative to the fatigue damage P% The damage degree D(t1) is obtained by determining the number of times of overshoot N(t1) which causes the damage degree D(t1) to occur, and the probability density function f(N)' of the obtained fatigue life is obtained from 0 to N(t1) The integral is obtained, and the probability of destruction F(t1)' after the update is obtained. 如申請專利範圍第7或8項所述之運作輔助裝置,其中,前述故障風險評估/更新部,藉由事先記憶好之壽命的機率密度函數f(t1)及前述概度,利用貝氏定理,更新壽命的密度函數,將更新後的機率密度函數f(t)’從0至t1積分而求出更新後的破壞機率F(t1)’,並將其乘上影響度C(p)來更新故障風險。 The operation assisting device according to claim 7 or 8, wherein the failure risk assessment/update unit utilizes the Bayesian theorem by using the probability density function f(t1) of the lifetime that is memorized in advance and the aforementioned generality. , update the density function of the life, and integrate the updated probability density function f(t)' from 0 to t1 to obtain the updated probability of failure F(t1)', and multiply it by the influence degree C(p) Update the risk of failure. 如申請專利範圍第7或8項所述之運作輔助裝置,其中, 前述故障風險評估/更新部,依據事先由製品的計測資料選擇出之和破壞有關連的複數個物理量,藉由依異常時的運轉資料群而得之破壞的機率密度函數與前述概度,利用貝氏定理,更新破壞的機率密度函數,對更新後的機率密度函數,描繪目前時間點的運轉資料的位置,藉此求出目前時間點t1下之更新後的破壞機率F(t1)’。 The operation auxiliary device as described in claim 7 or 8, wherein The foregoing failure risk assessment/update unit selects and destroys a plurality of physical quantities related to each other based on the measurement data of the product in advance, and uses the probability density function of the damage according to the operation data group at the time of the abnormality and the aforementioned generality, The theorem, update the probability density function of the damage, and describe the position of the running data at the current time point for the updated probability density function, thereby obtaining the updated probability of failure F(t1)' at the current time point t1. 如申請專利範圍第7項所述之運作輔助裝置;及具備:外部資料庫,事先記憶有和製品無關之外部資料,前述維護/運用腳本制定部,利用故障風險的趨勢曲線、及維護/運用腳本、及外部資料,演算將來的故障風險。 For example, the operation auxiliary device described in item 7 of the patent application; and having an external database, which previously stores external data irrelevant to the product, the maintenance/application scripting department, using the trend curve of the risk of failure, and maintenance/operation Scripts, and external materials, calculate the risk of future failures. 如申請專利範圍第4或12項所述之運作輔助裝置,其中,前述外部資料,包含事先計算出之氣象及/或海象的將來預測資料、資源的供給預測資料、資源的埋藏預測資料的任一或複數。 The operation auxiliary device as described in claim 4 or 12, wherein the external data includes pre-calculated meteorological and/or walrus future prediction data, resource supply prediction data, and resource burial prediction data. One or plural. 如申請專利範圍第7項所述之運作輔助裝置;及前述維護/運用腳本制定部,作成破壞機率的趨勢曲線、及/或損傷度的趨勢曲線。 The operation assisting device described in claim 7 of the patent application; and the maintenance/application scripting unit, the trend curve for destroying the probability and/or the trend curve for the damage degree. 一種風力發電系統,具備:如申請專利範圍第7項所述之運作輔助裝置;及 第1風力發電機,為具有複數個感測器之對象製品;第2風力發電機,具有複數個感測器,和前述第1風力發電機為同型機或類似機。 A wind power generation system comprising: the operation auxiliary device as described in claim 7; and The first wind power generator is a target product having a plurality of sensors; the second wind power generator has a plurality of sensors, and the first wind power generator is a same type machine or the like.
TW106104632A 2016-03-25 2017-02-13 Operational auxiliary device and wind power generation system TWI607328B (en)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016061725A JP2019113883A (en) 2016-03-25 2016-03-25 Utilization assisting apparatus and wind power generation system

Publications (2)

Publication Number Publication Date
TW201734870A true TW201734870A (en) 2017-10-01
TWI607328B TWI607328B (en) 2017-12-01

Family

ID=59899924

Family Applications (1)

Application Number Title Priority Date Filing Date
TW106104632A TWI607328B (en) 2016-03-25 2017-02-13 Operational auxiliary device and wind power generation system

Country Status (3)

Country Link
JP (1) JP2019113883A (en)
TW (1) TWI607328B (en)
WO (1) WO2017163561A1 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI688846B (en) * 2017-12-12 2020-03-21 日商日立製作所股份有限公司 Maintenance plan generation system
TWI709092B (en) * 2018-03-16 2020-11-01 日商日立製作所股份有限公司 Failure probability evaluation system and method
TWI760791B (en) * 2019-07-19 2022-04-11 日商三菱動力股份有限公司 Operation comparison device, operation comparison method and program

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019175984A1 (en) * 2018-03-13 2019-09-19 三菱電機株式会社 Failure level calculation device and failure level calculation system
WO2020030235A1 (en) 2018-08-07 2020-02-13 Vestas Wind Systems A/S Operating a wind turbine during service
JP6753442B2 (en) * 2018-09-10 2020-09-09 日本電気株式会社 Model generator, model generator, and program
JP7248877B2 (en) * 2018-12-21 2023-03-30 日立建機株式会社 Selection system for maintenance parts for working machines
CN109858689A (en) * 2019-01-21 2019-06-07 北京航空航天大学 A kind of Product Assembly system health risk analysis method of reliability guiding
CN110084465B (en) * 2019-03-11 2022-03-18 杭州电子科技大学 Wind power generation system cost/power supply reliability evaluation method based on energy storage
JP7145821B2 (en) * 2019-07-05 2022-10-03 株式会社日立製作所 Failure probability evaluation system and method
JP7429619B2 (en) * 2020-08-17 2024-02-08 株式会社日立製作所 Crack growth evaluation device and crack growth evaluation program
WO2023026668A1 (en) * 2021-08-26 2023-03-02 株式会社カネカ Monitoring system
CN114183314A (en) * 2021-11-25 2022-03-15 吉林省电力科学研究院有限公司 Wind turbine generator opportunity maintenance method based on reliability
CN115130702B (en) * 2022-09-02 2022-12-02 山东汇泓纺织科技有限公司 Textile machine fault prediction system based on big data analysis
JP2024035858A (en) * 2022-09-05 2024-03-15 株式会社日立製作所 Risk data management device and risk data management method
KR102655555B1 (en) * 2023-11-20 2024-04-11 한국전력기술 주식회사 System and method for maintenance planning in offshore wind farm using deep reinforcement learning

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4528533B2 (en) * 2004-02-03 2010-08-18 株式会社東芝 Nuclear power plant equipment inspection method and nuclear power plant equipment inspection system
JP4977064B2 (en) * 2008-03-12 2012-07-18 株式会社東芝 Maintenance plan support system
US8326577B2 (en) * 2011-09-20 2012-12-04 General Electric Company System and method for predicting wind turbine component failures
CN102637289B (en) * 2012-04-16 2016-02-24 天津大学 Containing the security value appraisal procedure of the Power System Planning scheme of large-scale wind power
US10861090B2 (en) * 2013-11-27 2020-12-08 Apple Inc. Provisioning of credentials on an electronic device using passwords communicated over verified channels
CN103914791A (en) * 2014-03-17 2014-07-09 中国南方电网有限责任公司超高压输电公司检修试验中心 Electrical equipment state maintenance system

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI688846B (en) * 2017-12-12 2020-03-21 日商日立製作所股份有限公司 Maintenance plan generation system
TWI709092B (en) * 2018-03-16 2020-11-01 日商日立製作所股份有限公司 Failure probability evaluation system and method
TWI760791B (en) * 2019-07-19 2022-04-11 日商三菱動力股份有限公司 Operation comparison device, operation comparison method and program
TWI802281B (en) * 2019-07-19 2023-05-11 日商三菱動力股份有限公司 Operation comparison device, operation comparison method and program

Also Published As

Publication number Publication date
JP2019113883A (en) 2019-07-11
WO2017163561A1 (en) 2017-09-28
TWI607328B (en) 2017-12-01

Similar Documents

Publication Publication Date Title
TWI607328B (en) Operational auxiliary device and wind power generation system
US11835030B2 (en) Methods and systems for predicting risk of observable damage in wind turbine gearbox components
JP6865189B2 (en) Failure probability evaluation system and method
Nielsen et al. On risk-based operation and maintenance of offshore wind turbine components
JP7320368B2 (en) FAILURE PREDICTION DEVICE, FAILURE PREDICTION METHOD AND COMPUTER PROGRAM
US7895016B2 (en) System and method for wind turbine health management
Besnard et al. Reliability-Centred Asset Maintenance—A step towards enhanced reliability, availability, and profitability of wind power plants
Butler et al. Exploiting SCADA system data for wind turbine performance monitoring
Tiddens et al. The adoption of prognostic technologies in maintenance decision making: a multiple case study
US11796992B2 (en) Condition-based method for malfunction prediction
TWI546762B (en) Wind power generation equipment of the stress estimation device and wind power equipment, the stress estimation method, wind power generation system
He et al. Condition-based maintenance optimization for multi-component systems considering prognostic information and degraded working efficiency
Zhang et al. Condition based maintenance and operation of wind turbines
Halme et al. Fault tree analysis for maintenance needs
Bangalore Load and risk based maintenance management of wind turbines
Igba et al. Through-life engineering services: a wind turbine perspective
Gray et al. A practical approach to the use of scada data for optimized wind turbine condition based maintenance
Geiss et al. Reliability-centered asset management of wind turbines—A holistic approach for a sustainable and cost-optimal maintenance strategy
Freitag et al. A concept for the dynamic adjustment of maintenance intervals by analysing heterogeneous data
Li et al. Advances in Intelligent Reliability and Maintainability of Energy Infrastructure Assets
CN112668794A (en) Maintenance optimization method for aircraft engine under uncertain life prediction information
Sayar et al. Real-time prediction of electricity distribution network status using artificial neural network model: A case study in Salihli (Manisa, Turkey)
Ambühl et al. Different reliability assessment approaches for wave energy converters
WO2023286659A1 (en) Failure predicting device, failure predicting method, and program
WO2024042733A1 (en) Failure probability evaluation system