JP5151556B2 - Process analysis apparatus, process analysis method, and process analysis program - Google Patents

Process analysis apparatus, process analysis method, and process analysis program Download PDF

Info

Publication number
JP5151556B2
JP5151556B2 JP2008049807A JP2008049807A JP5151556B2 JP 5151556 B2 JP5151556 B2 JP 5151556B2 JP 2008049807 A JP2008049807 A JP 2008049807A JP 2008049807 A JP2008049807 A JP 2008049807A JP 5151556 B2 JP5151556 B2 JP 5151556B2
Authority
JP
Japan
Prior art keywords
signal
signals
process analysis
causal
causal structure
Prior art date
Legal status (The legal status 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 status listed.)
Expired - Fee Related
Application number
JP2008049807A
Other languages
Japanese (ja)
Other versions
JP2009064407A (en
Inventor
真希 遠藤
広典 稲葉
浩輔 鶴田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Omron Corp
Original Assignee
Omron Corp
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 Omron Corp filed Critical Omron Corp
Priority to JP2008049807A priority Critical patent/JP5151556B2/en
Publication of JP2009064407A publication Critical patent/JP2009064407A/en
Application granted granted Critical
Publication of JP5151556B2 publication Critical patent/JP5151556B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
    • 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
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Description

本発明は、生産ライン内の装置間の関係をモデル化する技術に関する。   The present invention relates to a technique for modeling a relationship between apparatuses in a production line.

工場などの生産ラインで異常が発生すると、ライン全体が停止し生産が止まってしまう。そこで、保全員は定期的に生産ライン内を巡回し、異常が発生していないかあるいは異常の予兆がないかを検査する。なお、異常の検知は保全員の五感に頼ることが多かった。   When an abnormality occurs in a production line such as a factory, the entire line stops and production stops. Therefore, maintenance personnel periodically patrol the production line and inspect whether there is an abnormality or a sign of abnormality. The detection of abnormalities often relied on the five senses of maintenance personnel.

何らかの異常あるいはその予兆を検知したときに、異常が検知された装置よりも前の工程の装置に真の異常の原因が存在する場合がある。したがって、真の異常の原因を特定するためには、生産ライン内の装置間の関係を理解している必要がある。生産ラインにおいては、装置の数が膨大(数百以上)となる場合もあるため装置間の関係を把握することは困難であり、熟練した保全員の経験や勘に頼っている。   When an abnormality or a sign of the abnormality is detected, there may be a real cause of the abnormality in the device in the process prior to the device in which the abnormality is detected. Therefore, in order to identify the cause of the true abnormality, it is necessary to understand the relationship between the devices in the production line. In the production line, the number of devices may be enormous (several hundreds or more), so it is difficult to grasp the relationship between the devices, and relies on the experience and intuition of skilled maintenance personnel.

装置間の関係は設計情報によって規定されるが、装置の劣化、ロットの種類、改良保全におけるパラメータ調整などにより日々変化する。したがって、設計情報だけでは非熟練の保全員が装置間の関係を把握しきれない。   The relationship between devices is defined by design information, but changes daily due to device deterioration, lot type, parameter adjustment in improved maintenance, and the like. Therefore, unskilled maintenance personnel cannot grasp the relationship between devices only by design information.

生産ラインに異常が発生した場合の修理時間の大半は、故障原因の究明に費やされる。しかも、この故障原因の究明に要する時間は保全員の熟練度に応じて大幅に異なり、非熟練者の場合は熟練者の数倍程度時間がかかってしまうこともある。このような状況にも拘わらず、最近は非熟練の保全員が一人で保全しなければならない機会も増えている。そこで、熟練者の経験や勘に頼ることなく、コンピュータを利用して保全業務を容易にする技術が研究されている。   Most of the repair time when an abnormality occurs in the production line is spent investigating the cause of the failure. Moreover, the time required for investigating the cause of the failure varies greatly depending on the skill level of maintenance personnel, and in the case of a non-skilled person, it may take several times as long as the skilled person. In spite of this situation, recently, there are increasing opportunities for unskilled maintenance personnel to maintain alone. Therefore, a technique for making maintenance work easier by using a computer without relying on the experience and intuition of an expert is being researched.

例えば、特許文献1には、工程状態データと品質検査データとから、プロセス−品質モデルを作成する技術が開示されている。この技術によれば、工程状態データを監視することで、異常品の発生を検知することができる。   For example, Patent Document 1 discloses a technique for creating a process-quality model from process state data and quality inspection data. According to this technology, it is possible to detect the occurrence of an abnormal product by monitoring the process state data.

また、特許文献2も生産工程と品質の関係を表すモデルを作成する技術を開示する。特許文献2では、製造履歴と製品品質履歴の相関強度を求める。相関強度と製造順序情報を利用して、製造ラインの工程の因果構造モデルを作成する。この因果構造モデルを用いて、製品品質変動の原因となった根本的な変動原因を絞り込むことができる。
特開2005−197323号公報 特開2006−65598号公報
Patent Document 2 also discloses a technique for creating a model representing the relationship between the production process and quality. In Patent Document 2, the correlation strength between the manufacturing history and the product quality history is obtained. A causal structure model of the production line process is created using the correlation strength and the production order information. Using this causal structure model, it is possible to narrow down the fundamental causes of fluctuations that have caused the product quality fluctuations.
JP 2005-197323 A JP 2006-65598 A

しかしながら、上記のような従来技術の場合には、下記のような問題が生じていた。すなわち、生産ラインの異常は、品質不良品が生産されるだけではなく、生産の速度が遅くなるという症状で現れることもある。特許文献1の方法では、品質に現れない程度の装置の異常をモデルで表すことができない。   However, in the case of the prior art as described above, the following problems have occurred. That is, an abnormality in the production line may appear not only as a defective product but also as a symptom that the production speed is slow. In the method of Patent Document 1, it is not possible to represent an abnormality of the apparatus that does not appear in the quality with a model.

特許文献2の方法では、因果構造モデルを作成しており、装置間の関係を把握することが可能となる。しかしながら、この因果構造モデルを作成するためには、製造ラインの設計情報に基づく工程間順序情報をシステムに与える必要がある。したがって、あらかじめ設計情報を取得できない場合や、設計情報と実際の製造ラインが異なる動作をする場合に
は、正しい因果構造モデルを作成することができない。
In the method of Patent Document 2, a causal structure model is created, and the relationship between devices can be grasped. However, in order to create this causal structure model, it is necessary to give order information between processes based on the design information of the production line to the system. Therefore, a correct causal structure model cannot be created when design information cannot be acquired in advance or when the design information and the actual production line operate differently.

現実の製造ラインでは、装置の劣化や、ロットの種類の変化、改良保全などによって、動作順序が変化する。これを毎回取得して設計情報を修正するのは、手間であり非現実的である。   In an actual production line, the operation order changes due to device deterioration, lot type change, improved maintenance, and the like. It is troublesome and unrealistic to correct the design information by acquiring this every time.

本発明は上記実情に鑑みてなされたものであって、その目的とするところは、製造ライン内の装置間の関係を容易にモデル化することにある。   The present invention has been made in view of the above circumstances, and an object of the present invention is to easily model a relationship between apparatuses in a production line.

上記目的を達成するために本発明に係る工程解析装置は、略周期的に入出力される複数の信号間の関係を解析する。解析の対象となる信号は、例えば、複数の製造装置とこれらの製造装置を制御する制御装置から成る生産ラインでの、製造装置と制御装置との間で入出力される制御信号や応答信号などの信号である。このような生産ラインでタクト生産方式が行われる場合、複数の製造装置と制御装置との間で入出力される信号は、略周期的にやりとりされるものである。   In order to achieve the above object, a process analysis apparatus according to the present invention analyzes a relationship between a plurality of signals input and output substantially periodically. Signals to be analyzed include, for example, control signals and response signals that are input and output between the manufacturing apparatus and the control apparatus in a production line that includes a plurality of manufacturing apparatuses and a control apparatus that controls these manufacturing apparatuses. Signal. When the tact production method is performed in such a production line, signals input / output between the plurality of manufacturing apparatuses and the control apparatus are exchanged substantially periodically.

本発明に係る工程解析装置は、信号取得手段と特徴量算出手段と発生順序抽出手段と因果構造作成手段とを備える。信号取得手段は、複数の製造装置と制御装置との間で入出力される信号(制御信号と応答信号)を取得する。特徴量算出手段は、取得された信号から特徴量を算出する。ある信号の特徴量は、任意の基準信号との間の立ち上がり時刻の差として算出することができる。なお、上記のように信号は周期的に入出力されるので、特徴量は各サイクルにおいて算出される。   The process analysis apparatus according to the present invention includes a signal acquisition unit, a feature amount calculation unit, an occurrence order extraction unit, and a causal structure creation unit. The signal acquisition means acquires signals (control signals and response signals) input / output between the plurality of manufacturing apparatuses and the control apparatus. The feature amount calculating means calculates a feature amount from the acquired signal. A feature amount of a certain signal can be calculated as a difference in rising time from an arbitrary reference signal. Note that, as described above, signals are input and output periodically, so that the feature amount is calculated in each cycle.

発生順序抽出手段は、取得された信号の発生順序を抽出する。ここで、発生順序は、各周期における信号の発生順を求め、それを全期間について平均した値に基づいて求めることができる。   The generation order extraction unit extracts the generation order of the acquired signals. Here, the generation order can be determined based on a value obtained by calculating the generation order of signals in each period and averaging them over the entire period.

因果構造作成手段は、特徴量間の相関と信号の発生順序に基づいて、信号間の因果構造を作成する。因果構造作成手段は、各信号について、その信号よりも発生順序が前の信号のうち、その信号との相関が最も高い信号を、その信号の原因として因果構造を作成することが好適である。   The causal structure creation means creates a causal structure between the signals based on the correlation between the feature quantities and the signal generation order. It is preferable that the causal structure creation means creates a causal structure for each signal, using a signal having the highest correlation with the signal among the signals in the order of generation before the signal as the cause of the signal.

このように、信号の発生の先後関係と信号間の相関に着目することで、信号間の因果構造(因果構造ツリー)を作成することができる。因果構造の作成は製造装置と制御装置との間でやりとりされる信号を取得することによって行われるので、装置の劣化や改良保全によるパラメータ修正など動作順序が変化した場合であっても、自動的に因果構造が作成され人の手を煩わすことがない。また、異常が検知された場合に、この因果構造を参照することで、真の異常の原因を探す際に点検すべき製造装置の範囲を絞り込むことが可能となる。   Thus, by paying attention to the relationship between signal generation and correlation, it is possible to create a causal structure (causal structure tree) between signals. Since the causal structure is created by acquiring signals exchanged between the manufacturing equipment and the control equipment, even if the operating order changes, such as equipment deterioration or parameter correction due to improved maintenance, it is automatically performed. In addition, a causal structure is created so that human hands are not bothered. Further, when an abnormality is detected, it is possible to narrow down the range of manufacturing apparatuses to be inspected when searching for the cause of the true abnormality by referring to the causal structure.

また、本発明に係る工程解析装置は、上記のようにして作成した因果構造を表示する表示手段を備えることが好ましい。この表示手段は、信号に対応する名称を設計情報から取得し、その名称を用いて因果構造を表示することがさらに好ましい。   Moreover, it is preferable that the process analysis apparatus according to the present invention includes display means for displaying the causal structure created as described above. More preferably, the display means obtains a name corresponding to the signal from the design information and displays the causal structure using the name.

また、本発明に係る工程解析装置を利用して、異常が発生した場合に、信号間の因果構造のいずれの箇所で異常が発生したかを推定することが好ましい。具体的には、学習用の信号データを利用して、あらかじめ求められた各信号間の因果構造及び因果関係のある信号間の時間差を、学習結果として記憶する記憶手段と、診断用の信号データを利用して、因果関係のある信号間の時間差を算出し、学習結果と比較して診断用データにおける異常
かを推定する異常箇所推定手段と、をさらに有することが好ましい。
Moreover, it is preferable to estimate in which part of the causal structure between signals the abnormality has occurred using the process analysis apparatus according to the present invention. Specifically, using the signal data for learning, the storage means for storing the causal structure between signals obtained in advance and the time difference between signals having a causal relationship as learning results, and the signal data for diagnosis It is preferable to further include an abnormal point estimation unit that calculates a time difference between signals having causal relationships by using and estimates whether there is an abnormality in the diagnostic data in comparison with the learning result.

この場合、異常箇所推定手段は、診断用の信号データにおける因果関係のある信号間の時間差と、学習結果における時間差とのずれが所定のしきい値よりも大きい箇所を異常箇所であると推定することが好適である。   In this case, the abnormal part estimation means estimates that a part where the difference between the time difference between the causal signals in the diagnostic signal data and the time difference in the learning result is larger than a predetermined threshold is an abnormal part. Is preferred.

さらに、信号間の因果構造を、推定された異常箇所を強調しつつ表示する表示部を有することが好適である。   Furthermore, it is preferable to have a display unit that displays the causal structure between signals while emphasizing the estimated abnormal part.

このように、あらかじめ学習のデータを利用して、正常時における信号間の発生時間差を調べておき、異常が発生した場合に信号間の時間差が正常時よりも大きくずれている箇所を異常箇所としてユーザに提示することができる。したがって、異常が検知された場合に、ユーザは迅速に対処・修理することが可能となる。   In this way, using the learning data beforehand, the occurrence time difference between the signals at the normal time is examined, and when the abnormality occurs, the place where the time difference between the signals is greatly deviated from the normal time is regarded as the abnormal place. It can be presented to the user. Therefore, when an abnormality is detected, the user can quickly cope and repair.

なお、本発明は、上記手段の少なくとも一部を有する工程解析装置として捉えることができる。また、本発明は、上記処理の少なくとも一部を含む工程解析方法、または、かかる方法を実現するためのプログラムとして捉えることもできる。上記手段および処理の各々は可能な限り互いに組み合わせて本発明を構成することができる。   In addition, this invention can be grasped | ascertained as a process analysis apparatus which has at least one part of the said means. The present invention can also be understood as a process analysis method including at least a part of the above processing, or a program for realizing the method. Each of the above means and processes can be combined with each other as much as possible to constitute the present invention.

たとえば、本発明の一態様としての工程解析方法は、略周期的に入出力された複数の信号の関係を解析する工程解析方法であって、情報処理装置が、前記複数の信号を取得するステップと、取得された信号から特徴量を算出するステップと、取得された信号の発生順序を抽出するステップと、特徴量間の相関と信号の発生順序とに基づいて、信号間の因果構造を作成するステップと、を実行する   For example, the process analysis method as one aspect of the present invention is a process analysis method for analyzing a relationship between a plurality of signals input and output substantially periodically, and the information processing apparatus acquires the plurality of signals. And the step of calculating the feature quantity from the acquired signal, the step of extracting the generation order of the acquired signal, and the causal structure between the signals based on the correlation between the feature quantity and the signal generation order And step to perform

また、本発明の一態様としての工程解析プログラムは、略周期的に入出力された複数の信号の関係を解析する工程解析プログラムであって、情報処理装置に、前記複数の信号を取得するステップと、取得された信号から特徴量を算出するステップと、取得された信号の発生順序を抽出するステップと、特徴量間の相関と信号の発生順序とに基づいて、信号間の因果構造を作成するステップと、を実行させる。   The process analysis program as one aspect of the present invention is a process analysis program for analyzing a relationship between a plurality of signals input and output substantially periodically, and the information processing apparatus acquires the plurality of signals. And the step of calculating the feature quantity from the acquired signal, the step of extracting the generation order of the acquired signal, and the causal structure between the signals based on the correlation between the feature quantity and the signal generation order And executing a step.

本発明によれば、製造ライン内の装置間の関係を容易にモデル化することが可能となる。   According to the present invention, it is possible to easily model the relationship between devices in a production line.

以下に図面を参照して、この発明の好適な実施の形態を例示的に詳しく説明する。   Exemplary embodiments of the present invention will be described in detail below with reference to the drawings.

(第1の実施形態)
本実施形態に係る工程解析装置1は、製造ラインに係る複数の製造装置2と、これらの製造装置2を制御する制御装置であるプログラマブルロジックコントローラ(PLC)3との間でやりとりされる信号を取得し、これらの信号間の因果構造を作成する。
(First embodiment)
The process analysis apparatus 1 according to the present embodiment transmits signals exchanged between a plurality of manufacturing apparatuses 2 related to a manufacturing line and a programmable logic controller (PLC) 3 that is a control apparatus that controls the manufacturing apparatuses 2. Acquire and create a causal structure between these signals.

図1は、工程解析装置の機能構成を示す図である。工程解析装置1は、信号計測収集部11,特徴量算出部12,発生順位抽出部13,順位付相関係数算出部14,因果構造獲得部15,および表示部17を備える。工程解析装置1はコンピュータ(情報処理装置)によって構成され、メモリに格納されたプログラムがCPU(中央演算装置)によって実行されることで上記の各機能部が実現される。なお、工程解析装置1は、1台のコンピュータによって構成されても良いし、ネットワーク接続された複数のコンピュータによって構成されても良い。   FIG. 1 is a diagram illustrating a functional configuration of the process analysis apparatus. The process analysis apparatus 1 includes a signal measurement collection unit 11, a feature amount calculation unit 12, an occurrence rank extraction unit 13, a ranked correlation coefficient calculation unit 14, a causal structure acquisition unit 15, and a display unit 17. The process analysis apparatus 1 is configured by a computer (information processing apparatus), and each function unit described above is realized by a program stored in a memory being executed by a CPU (central processing unit). In addition, the process analysis apparatus 1 may be comprised by one computer, and may be comprised by the some computer connected by the network.

[信号計測収集部]
PLC3は、センサや応答信号の入力を受け付け、ラダーロジックに従って製造装置2を制御するための制御信号を送る。信号計測収集部11は、製造装置2とPLC3の各接点との間でやりとりされる信号を収集する。信号計測収集部11は、バスを介してPLC3と接続されており、PLC3のI/Oメモリを読み出してファイルとして記憶するユニットである。なお、信号計測収集部11は、バス接続以外にも、イーサネット(登録商標)やシリアル接続などによって、PLC3からデータを収集しても良い。また、信号収集計測部11は、PLC3の信号データを取得可能であれば、PLC3と直接接続されている必要はない。例えば、PLC3の信号データをCD−ROM,DVD−ROM,メモリカードなどの記憶媒体に一旦格納し、この記憶媒体から信号データを取得しても構わない。
[Signal Measurement Collection Unit]
PLC3 receives the input of a sensor and a response signal, and sends the control signal for controlling the manufacturing apparatus 2 according to ladder logic. The signal measurement collection unit 11 collects signals exchanged between the manufacturing apparatus 2 and each contact point of the PLC 3. The signal measurement collection unit 11 is connected to the PLC 3 via a bus, and is a unit that reads the I / O memory of the PLC 3 and stores it as a file. The signal measurement collection unit 11 may collect data from the PLC 3 by Ethernet (registered trademark) or serial connection in addition to the bus connection. Moreover, the signal collection measurement part 11 does not need to be directly connected with PLC3, if the signal data of PLC3 can be acquired. For example, the PLC 3 signal data may be temporarily stored in a storage medium such as a CD-ROM, DVD-ROM, or memory card, and the signal data may be acquired from the storage medium.

図2は、信号計測収集部11の動作の流れを示すフローチャートである。信号計測収集部11は、製造装置2とPLC3との間でやりとりされる信号を計測収集する(S201)。図3は、信号計測収集部11によって収集される信号の例を示す図である。信号計測収集部11は、各時点における信号のオン/オフを、信号情報101としてメモリに記憶し保存する(S202)。図4に、信号情報101の例を示す。信号情報101には、接点ごとに、各時刻tにおける信号がオン(1)であるかオフ(0)であるかが格納される。   FIG. 2 is a flowchart showing an operation flow of the signal measurement collection unit 11. The signal measurement collection unit 11 measures and collects signals exchanged between the manufacturing apparatus 2 and the PLC 3 (S201). FIG. 3 is a diagram illustrating an example of signals collected by the signal measurement collection unit 11. The signal measurement collection unit 11 stores the signal on / off at each time point in the memory as the signal information 101 and stores it (S202). FIG. 4 shows an example of the signal information 101. The signal information 101 stores, for each contact, whether the signal at each time t is on (1) or off (0).

製造装置2からなる製造ラインでは、タクト生産方式によって製品の製造を行っており、各製造装置2の動作は一定の周期(タクトサイクル)に同期されている。したがって、製造装置2とPLC3との間でやりとりされる信号は、このタクトサイクルの周期性を持つ。   In a production line composed of the production apparatuses 2, products are produced by a tact production method, and the operation of each production apparatus 2 is synchronized with a certain cycle (tact cycle). Therefore, the signal exchanged between the manufacturing apparatus 2 and the PLC 3 has the periodicity of this tact cycle.

ここで、任意の信号を基準信号として定義すると、その基準信号の1サイクルの中で動作する他の全ての信号は、そのサイクル内における挙動に関するデータとして捉えることができる。図5は、信号のサイクルについて説明する図である。例えば、図3において一番上の信号を基準信号として定義すると、各サイクルはその基準信号の立ち上がりによって区切られる(図では、i−1サイクル、iサイクル、i+1サイクルの3サイクルが描かれている)。このように定義されたサイクルによって、他の信号のオン/オフがどのサイクルに係るものであるのかが定義される。   Here, when an arbitrary signal is defined as a reference signal, all other signals operating in one cycle of the reference signal can be regarded as data relating to behavior in the cycle. FIG. 5 is a diagram illustrating a signal cycle. For example, when the top signal in FIG. 3 is defined as a reference signal, each cycle is divided by the rising edge of the reference signal (in the figure, three cycles of i−1 cycle, i cycle, and i + 1 cycle are depicted. ). The cycle defined in this way defines which cycle the other signals are turned on / off.

なお、どの信号を基準信号として採用するかに関しては、最初に稼働する装置の信号を基準信号として採用する方法が考えられる。   As to which signal is adopted as the reference signal, a method of adopting the signal of the first operating device as the reference signal can be considered.

[特徴量算出部]
次に、特徴量算出部12について説明する。特徴量算出部12は、信号計測収集部11によって収集された接点ごとの時系列信号の特徴量を算出する。時系列信号の特徴量として採用する値としては種々のものが考えられるが、本工程解析装置1では信号間の関係に着目しているので、基準信号と対象信号との状態変化時刻の差を、対象信号の特徴量とする。より具体的には、図6Aに示すように、基準信号の立ち上がり時刻と、対象信号の立ち上がり時刻の差を特徴量とすることができる。本実施形態では、このような特徴量を採用するが、図6Bに示すように、基準信号の立ち上がり時刻と対象信号の立ち下がり時刻の差(a)、基準信号の立ち下がり時刻と対象信号の立ち上がり時刻の差(b)、または基準信号の立ち下がり時刻と対象信号の立ち下がり時刻の差(c)を特徴量としてもよい。
[Feature amount calculator]
Next, the feature amount calculation unit 12 will be described. The feature amount calculation unit 12 calculates the feature amount of the time series signal for each contact collected by the signal measurement collection unit 11. Although various values can be considered as the feature value of the time series signal, since the present process analysis apparatus 1 focuses on the relationship between the signals, the difference in the state change time between the reference signal and the target signal is calculated. The feature amount of the target signal. More specifically, as shown in FIG. 6A, the difference between the rising time of the reference signal and the rising time of the target signal can be used as the feature amount. In the present embodiment, such a feature amount is adopted. As shown in FIG. 6B, the difference (a) between the rising time of the reference signal and the falling time of the target signal, the falling time of the reference signal and the target signal. The difference between the rise times (b) or the difference between the fall time of the reference signal and the fall time of the target signal (c) may be used as the feature amount.

なお、信号の特徴量は、接点ごとおよびサイクルごとに定義されるものである。   Note that the feature amount of the signal is defined for each contact and each cycle.

図7は、特徴量算出部12が行う特徴量算出処理の流れを示すフローチャートである。特徴量算出部12は、まず、対象信号を表すインデックスであるjに1をセットする(S701)。次に、信号計測収集部11が作成した信号情報101(図4参照)から、信号sの状態変化時刻を抽出して、状態変化時刻情報102を作成する(S702)。すなわち、信号情報101を検査して、信号の状態がオン(1)からオフ(0)、または、オフ(0)からオン(1)に変化する時刻を抽出して状態変化時刻情報102として保存する。状態変化時刻情報102の例を図8に示す。状態変化時刻情報102には、接点ごとに信号のオン/オフ状態が切り替えられる時刻が格納されている。 FIG. 7 is a flowchart showing the flow of the feature amount calculation process performed by the feature amount calculation unit 12. The feature amount calculation unit 12 first sets 1 to j, which is an index representing the target signal (S701). Next, the state change time information 102 is created by extracting the state change time of the signal s j from the signal information 101 (see FIG. 4) created by the signal measurement collection unit 11 (S702). That is, the signal information 101 is inspected, the time when the signal state changes from on (1) to off (0), or from off (0) to on (1) is extracted and stored as state change time information 102. To do. An example of the state change time information 102 is shown in FIG. The state change time information 102 stores the time at which the signal on / off state is switched for each contact.

特徴量算出部12は、N個の信号全てについて状態変化時刻の抽出が完了したか判定し(S703)、処理していない信号がある場合(S703−YES)はjをインクリメントして(S704)、次の信号について状態変化時刻を抽出する。   The feature amount calculation unit 12 determines whether or not the state change time extraction has been completed for all N signals (S703). If there is an unprocessed signal (S703-YES), j is incremented (S704). The state change time is extracted for the next signal.

N個全ての信号について状態変化時刻の抽出が完了したら(S703−NO)、状態変化時刻情報102から特徴量を算出する処理を行う。まず、N個の信号の中から1つの信号を基準信号として設定する(S705)。ここでは、最も早い時刻に状態の変化があった信号を基準信号として選択する。   When the extraction of the state change times for all N signals is completed (S703-NO), processing for calculating the feature amount from the state change time information 102 is performed. First, one signal out of N signals is set as a reference signal (S705). Here, the signal whose state has changed at the earliest time is selected as the reference signal.

そして、対象信号を表すインデックスであるjに1をセットし、対象とするサイクルを表すインデックスであるiに1をセットする(S706)。特徴量算出部12は、状態変化時刻情報102から、サイクルiにおける信号sの特徴量を算出し、特徴量データ103として格納する(S707)。具体的には、サイクルiにおける基準信号の立ち上がり時刻と、信号sの立ち上がり時刻を取得し、その差をサイクルiにおける信号sの特徴量fとして算出する。特徴量データ103の例を図9に示す。特徴量データ103には、信号(接点)ごと、サイクルごとに特徴量の値が格納されている。 Then, 1 is set to j, which is an index representing the target signal, and 1 is set to i, which is an index representing the target cycle (S706). The feature amount calculation unit 12 calculates the feature amount of the signal s j in the cycle i from the state change time information 102 and stores it as the feature amount data 103 (S707). Specifically, the rising time of the reference signal in cycle i and the rising time of signal s j are acquired, and the difference is calculated as feature quantity f j of signal s j in cycle i. An example of the feature amount data 103 is shown in FIG. The feature value data 103 stores the value of the feature value for each signal (contact point) and for each cycle.

次に、N個全ての信号について特徴量の抽出が完了したか判定し(S708)、処理していない信号がある場合は(S708−YES)は、jをインクリメントして(S709)次の信号の特徴量を抽出する。   Next, it is determined whether or not feature amount extraction has been completed for all N signals (S708). If there is an unprocessed signal (S708-YES), j is incremented (S709) and the next signal Extract feature values.

サイクルiにおいてN個全ての信号について特徴量の抽出が完了した場合(S708−NO)は、M個全てのサイクルについて特徴量の抽出が完了したか判定する(S710)。まだ特徴量を抽出していないサイクルがある場合には、サイクルのインデックスiをインクリメントし、信号のインデックスjに1をセットして(S711)、次のサイクルについて特徴量の抽出を行う。   When feature value extraction has been completed for all N signals in cycle i (S708-NO), it is determined whether feature value extraction has been completed for all M cycles (S710). If there is a cycle in which the feature quantity has not yet been extracted, the cycle index i is incremented, 1 is set to the signal index j (S711), and the feature quantity is extracted for the next cycle.

全てのサイクルについて特徴量の抽出が完了した場合(S710−NO)に、特徴量算出部12による特徴量算出処理が終了する。   When the feature amount extraction is completed for all the cycles (S710—NO), the feature amount calculation process by the feature amount calculation unit 12 ends.

[発生順位抽出部]
発生順位抽出部13は、各信号の立ち上がりの発生順位を抽出する。ここで、ある信号の発生順位は、サイクルによって変化することがあるので、各サイクルについて発生順位を求め、その値(平均順位)の大きさに応じて発生順位を決定する。
[Occurrence rank extraction unit]
The generation order extraction unit 13 extracts the generation order of the rise of each signal. Here, since the generation order of a certain signal may change depending on the cycle, the generation order is obtained for each cycle, and the generation order is determined according to the magnitude of the value (average rank).

発生順位抽出部13が行う、発生順位抽出処理の詳細を図10〜12を参照して説明する。図10は、発生順位抽出処理の流れを示すフローチャートである。図11は、各サイクルにおける発生順位を説明する図である。図12は、サイクルごとの発生順位の平均値(重み)と、算出する発生順位の関係を示す図である。   Details of the generation order extraction process performed by the generation order extraction unit 13 will be described with reference to FIGS. FIG. 10 is a flowchart showing the flow of occurrence order extraction processing. FIG. 11 is a diagram for explaining the order of occurrence in each cycle. FIG. 12 is a diagram illustrating the relationship between the average value (weight) of the occurrence order for each cycle and the occurrence order to be calculated.

発生順位抽出部13は、まず、各サイクルにおける、各信号の発生順位を求める。そこで、サイクルを表すインデックスiに1をセットする(S1001)。そして、状態変化時刻情報102(図8)を参照して、サイクルiにおける各信号の立ち上がり時刻を抽出する(S1002)。そして、サイクルiにおける各信号の発生順位を求め、サイクルごとの発生順位104として格納する(S1003)。図11を例に取ると、サイクルi=1において、信号sの発生順位は「3」、信号sの発生順位は「7」、信号sの発生順位は「1」、等という具合になる。 The generation order extraction unit 13 first obtains the generation order of each signal in each cycle. Therefore, 1 is set to the index i representing the cycle (S1001). Then, referring to the state change time information 102 (FIG. 8), the rising time of each signal in cycle i is extracted (S1002). Then, the generation order of each signal in cycle i is obtained and stored as the generation order 104 for each cycle (S1003). Taking FIG. 11 as an example, in cycle i = 1, the generation order of signal s 1 is “3”, the generation order of signal s 2 is “7”, the generation order of signal s 3 is “1”, and so on. become.

M個全てのサイクルについて、サイクルごとの発生順位抽出が終了したか判定し(S1004)、まだ処理していないサイクルがある場合(S1004−YES)は、サイクルのインデックスiをインクリメントして(S1005)、次のサイクルについて処理を行う。   For all M cycles, it is determined whether the generation order extraction for each cycle has been completed (S1004). If there is a cycle that has not yet been processed (S1004-YES), the cycle index i is incremented (S1005). The next cycle is processed.

上述したように、サイクルごとに各信号の発生順にばらつきが見られる場合があり、図11では、サイクルi=1において発生順位が「3」であった信号sは、サイクルi=2においては発生順位が「2」になっている。 As described above, there may be variations in the order of generation of each signal in each cycle. In FIG. 11, the signal s 1 whose generation order is “3” in cycle i = 1 is different in cycle i = 2. The occurrence order is “2”.

全サイクルについて、サイクルごとの発生順位の抽出が完了した場合(S1004−NO)は、それぞれの信号について、全サイクルでの発生順位の平均を算出する(S1006)。この値は、次の式で求められる。なお、ord(i)は、サイクルiにおける信号sの発生順位を示す。 When the extraction of the occurrence rank for each cycle is completed for all cycles (S1004-NO), the average of the occurrence rank for all cycles is calculated for each signal (S1006). This value is obtained by the following formula. Note that ord j (i) indicates the generation order of the signal s j in the cycle i.

Figure 0005151556
Figure 0005151556

図11を例に具体的に説明すると、信号sのサイクルごとの発生順位は「3」、「2」、「3」であるので、平均順位は「2.67」となる。同様に信号sのサイクルごとの発生順位は「7」、「7」、「6」であるので、平均順位は「6.67」になる。その他の信号についても同様である。 Specifically, taking FIG. 11 as an example, the generation order of the signal s 1 for each cycle is “3”, “2”, and “3”, so the average order is “2.67”. Similarly, since the generation rank of the signal s 2 for each cycle is “7”, “7”, “6”, the average rank is “6.67”. The same applies to other signals.

そして、発生順位抽出部13は、サイクルごとの発生順位の平均を昇順にソートして、各信号の発生順位を決定し、発生順位105として記憶する(S1007)。図12は発生順位105の例を示す図であり、サイクルごとの発生順位の平均値(ord)の順に各信号の発生順位(ORD)が決められていることが分かる。 Then, the occurrence rank extraction unit 13 sorts the average of the occurrence ranks for each cycle in ascending order, determines the occurrence rank of each signal, and stores it as the occurrence rank 105 (S1007). FIG. 12 is a diagram showing an example of the generation order 105, and it can be seen that the generation order (ORD j ) of each signal is determined in the order of the average value (ord j ) of the generation order for each cycle.

[順位付相関係数算出部]
次に、順位付相関係数算出部14が行う、信号間の相関を求める処理について説明する。図13は、順位付相関係数算出部14が行う処理の流れを示すフローチャートである。また、図14は、その処理を説明する図である。
[Ranked correlation coefficient calculation unit]
Next, the process for obtaining the correlation between signals performed by the ranked correlation coefficient calculation unit 14 will be described. FIG. 13 is a flowchart showing the flow of processing performed by the ranked correlation coefficient calculation unit 14. FIG. 14 is a diagram for explaining the processing.

順位付相関係数算出部14は、発生順位抽出部13が抽出した発生順位105に基づいて、各信号sの特徴量データ103の並び替えを行う(S1301)。すなわち、次式で表されるように、変数x(i)を、発生順位ORDでラベリングした特徴量f(i)とする。 The ranked correlation coefficient calculation unit 14 rearranges the feature amount data 103 of each signal s j based on the occurrence rank 105 extracted by the occurrence rank extraction unit 13 (S1301). That is, as represented by the following expression, the variable x j (i) is a feature quantity f j (i) labeled with the occurrence order ORD j .

Figure 0005151556
Figure 0005151556

発生順位105が図12のように求められている場合には、図14に示すように、xは発生順位が「1」であるfとなり、xは発生順位が「2」であるfとなる。このように、特徴量データfを発生順に並び替えて、xを得る。発生順位に基づいて信号の並び替えを行っているので、j<kの場合、xはxよりも後に発生する(確率が高い)信号であることになる。 When the occurrence rank 105 is obtained as shown in FIG. 12, as shown in FIG. 14, x 1 is f 3 with the occurrence rank being “1”, and x 2 is the occurrence rank “2”. the f 5. In this way, the feature quantity data f j is rearranged in the order of occurrence to obtain x j . Since the rearrangement of the signals is performed based on the generation order, when j <k, x k is a signal that occurs after x j (having a high probability).

次に、順位付相関係数算出部14は、x(i)間の相関係数行列106を算出する(S1302)。相関係数行列106の各要素ρjkは、以下のようにして算出される。 Next, the ranked correlation coefficient calculation unit 14 calculates the correlation coefficient matrix 106 between x j (i) (S1302). Each element ρ jk of the correlation coefficient matrix 106 is calculated as follows.

Figure 0005151556
Figure 0005151556

算出される相関係数行列106の例を図14に示している。相関係数行列は対称行列となるので、下三角(または上三角)行列のみ求めれば十分である。また、後述するように相関係数行列ρjkにおいて、j≧3のみ求めればよい。 An example of the calculated correlation coefficient matrix 106 is shown in FIG. Since the correlation coefficient matrix is a symmetric matrix, it is sufficient to obtain only the lower triangular (or upper triangular) matrix. Further, as will be described later, it is only necessary to obtain j ≧ 3 in the correlation coefficient matrix ρ jk .

[因果構造獲得部]
因果構造獲得部15は、順位付相関係数算出部14によって算出された相関係数に基づいて、信号間の因果構造を作成する。図15は、因果構造獲得部15が行う因果構造獲得処理の流れを示すフローチャートである。
[Causal structure acquisition unit]
The causal structure acquisition unit 15 creates a causal structure between signals based on the correlation coefficient calculated by the ranked correlation coefficient calculation unit 14. FIG. 15 is a flowchart showing the flow of causal structure acquisition processing performed by the causal structure acquisition unit 15.

まず、因果構造としてx→xを定義する(S1501)。なお、「x→x」は、信号xの発生を原因として信号xが発生するということを意味する。ある信号の原因となる信号は、その信号よりも前に発生する信号であるため、xの原因となる信号は必ずxとなる。 First, x 1 → x 2 is defined as a causal structure (S1501). Note that “x j → x k ” means that the signal x k is generated due to the generation of the signal x j . Signal which causes some signals are the signals generated before the signal, the signal which causes x 2 always becomes x 1.

次に、信号を表すインデックスjに3をセットする(S1502)。そして、相関係数行列106を参照して、ρjk(k<j)が最大となるkを求める。この処理は、k<jを満たすx、すなわち、xよりも発生順序が前である信号のうちxと相関の最も高い信号を探索している。このようにしてkを求めて、因果関係リスト107に、x→xを格納する。 Next, 3 is set to the index j representing the signal (S1502). Then, with reference to the correlation coefficient matrix 106, k that maximizes ρ jk (k <j) is obtained. In this process, x k satisfying k <j, that is, a signal having the highest correlation with x j is searched for among signals whose generation order is earlier than x j . In this way, k is obtained, and x k → x j is stored in the causal relationship list 107.

全ての信号について処理をしたか判定し(S1505),処理していない信号がある場合(S1505−YES)には、jをインクリメントして(S1506)、次の信号についての原因となる信号を求める。   It is determined whether all the signals have been processed (S1505). If there is a signal that has not been processed (S1505-YES), j is incremented (S1506), and the signal that causes the next signal is obtained. .

図16を例に具体的に説明すると、まず最初はx→xが格納されている。そして、j=3として、相関係数ρ31,ρ32のうち最も大きいものを求める。ここでは、ρ31が最も大きいとする。この場合、x→xという因果関係が得られる。 Specifically, taking FIG. 16 as an example, first, x 1 → x 2 is stored. Then, with j = 3, the largest of the correlation coefficients ρ 31 and ρ 32 is obtained. Here, it is assumed that ρ 31 is the largest. In this case, a causal relationship of x 1 → x 3 is obtained.

次に、j=4として、相関係数ρ41、ρ42、ρ43のうち、最も大きいものを求める。ここでは、ρ42が最も大きいとする。この場合、x→xという因果関係が追加される。 Next, assuming j = 4, the largest one of the correlation coefficients ρ 41 , ρ 42 , ρ 43 is obtained. Here, it is assumed that ρ 42 is the largest. In this case, a causal relationship of x 2 → x 4 is added.

さらに、j=5として、相関係数ρ51、ρ52、ρ53、ρ54のうち、最も大きいものを求める。ここでは、ρ52が最も大きいとする。この場合、x→xという因果関係が追加される。 Further, assuming that j = 5, the largest of the correlation coefficients ρ 51 , ρ 52 , ρ 53 , ρ 54 is obtained. Here, it is assumed that ρ 52 is the largest. In this case, a causal relationship of x 2 → x 5 is added.

このようにして、全ての信号について処理をすることで、xを起点とするツリー状の因果構造リストが作成される。上記のように因果構造を作成しているので、ある信号の原因となる信号は1つのみであるのに対して、1つの信号が複数の信号の原因となる場合がある。 In this way, by the processing for all the signals, tree-like causal structure list which starts x 1 is created. Since the causal structure is created as described above, only one signal may cause a certain signal, whereas one signal may cause a plurality of signals.

[表示部]
表示部17は、上記のようにして作成された信号間の因果関係を表示する。この際、設計情報を参照して信号の接点名を取得し、接点名を用いて表示を行う。ここでいう設計情報とは、熟練者などの人の知識、製造ラインの設計図、PLCで保持されている情報(PLCのラダー、信号の名前、説明が記載された文字列、装置の制御に関するパラメータなど)のことをいう。
[Display section]
The display unit 17 displays the causal relationship between the signals created as described above. At this time, the contact name of the signal is acquired by referring to the design information, and the display is performed using the contact name. Design information here refers to knowledge of people such as experts, production line design drawings, information held in the PLC (PLC ladder, signal names, character strings with explanations, device control, etc. Parameter).

図17は、ユーザが因果構造を確認する際の出力画面のイメージである。図に示すように、製造ライン内の信号は、工程解析装置1によってその因果関係が解析され、ディスプレイ装置などの表示部17にその因果構造がグラフとして表示される。   FIG. 17 is an image of the output screen when the user confirms the causal structure. As shown in the figure, the causal relationship of the signals in the production line is analyzed by the process analysis apparatus 1 and the causal structure is displayed as a graph on the display unit 17 such as a display device.

図18は、表示部17が行う表示処理の流れを示すフローチャートである。表示部17は、因果関係リスト107からツリー構造を取得し、設計情報110から変数名に対応する接点名を取得する(S1901)。そして、変数名を信号の接点名で置き換えた因果構造をグラフ表示する(S1902)。なお、変数名を信号の接点名で置き換える場合には、次式にしたがって変数名と信号の接点名の対応付けを行う。   FIG. 18 is a flowchart showing the flow of display processing performed by the display unit 17. The display unit 17 acquires a tree structure from the causal relationship list 107 and acquires a contact name corresponding to the variable name from the design information 110 (S1901). Then, the causal structure in which the variable name is replaced with the signal contact name is displayed in a graph (S1902). When the variable name is replaced with the signal contact name, the variable name is associated with the signal contact name according to the following equation.

Figure 0005151556
Figure 0005151556

表示されるグラフは、図19に示すように、変数名が「シリンダ1出駆動」や「シリンダ1戻センサ」といった信号の接点名で置き換えられるため装置間の関係が分かりやすくなる。   In the displayed graph, as shown in FIG. 19, the variable name is replaced with a contact name of a signal such as “cylinder 1 output drive” or “cylinder 1 return sensor”, so that the relationship between the devices becomes easy to understand.

なお、表示部17は、図20に示すように因果構造の全体を表示しても良く、図21〜23に示すように因果構造上の任意の変数を選択したとき、その変数とxとを単線で結ぶことにより形成される経路を表示しても良い。図21は、変数x(信号の接点名は「シリンダ1出センサ」)から、その原因となる信号を変数xまで辿った経路を表示している。同様に、図22は、変数x(「シリンダ2出センサ」)から原因となる信号を変数xまで辿った経路を示している。図23も同様に、変数x(「シリンダ1戻センサ」)から原因となる信号を変数xまで辿った経路を示している。なお、図20〜23では、接点名の隣に変数名「x」を表示しているが、これは説明のためであり、実際にはこの変数名を表示する必要はない。 The display unit 17 may display the entire causal structure as shown in FIG. 20, when selecting any variable on the causal structure as shown in FIG. 21 to 23, and the variable and x 1 A route formed by connecting the lines with a single line may be displayed. FIG. 21 shows a path obtained by tracing the signal causing the variable x 3 (the contact name of the signal is “cylinder 1 output sensor”) to the variable x 1 . Similarly, FIG. 22 shows a path obtained by tracing the signal causing the variable x 4 (“cylinder 2 output sensor”) to the variable x 1 . Similarly, FIG. 23 shows a path obtained by tracing the signal causing the variable x 5 (“cylinder 1 return sensor”) to the variable x 1 . 20 to 23, the variable name “x j ” is displayed next to the contact name. However, this is for explanation, and it is not actually necessary to display this variable name.

〈実施形態の作用/効果〉
本実施形態に係る工程解析装置1によれば、実際に製造装置2とPLC3の間でやりと
りされる信号の時系列情報から、順序情報(先後関係)を取得し、この順序情報や特徴量の相関に基づいて、信号間の因果関係を作成することができる。この信号間の関係を見ると、装置間の関係を容易に理解できるので、故障やその予兆が検知された場合に、真の原因箇所を特定するために大いに役立つ。特に、図21〜23に示すような一部の因果関係のみを表示する機能によって、異常が見られた信号の原因を辿る経路のみを表示すれば、真の原因箇所が絞り込まれて表示される。保全員は表示された経路に関連する装置のみを調べればよく、保全作業が容易となる。
<Operation / Effect of Embodiment>
According to the process analysis apparatus 1 according to the present embodiment, order information (prior relationship) is acquired from time-series information of signals actually exchanged between the manufacturing apparatus 2 and the PLC 3, and the order information and the feature amount are obtained. Based on the correlation, a causal relationship between the signals can be created. By looking at the relationship between the signals, the relationship between the devices can be easily understood, so that when a failure or a sign thereof is detected, it is very useful for identifying the true cause. In particular, if only the path that follows the cause of the signal in which an abnormality has occurred is displayed by the function of displaying only a part of the causal relationships as shown in FIGS. 21 to 23, the true cause location is narrowed down and displayed. . Maintenance personnel only need to examine the devices related to the displayed route, and maintenance work is facilitated.

また、実際の製造ラインでは、装置の劣化や改良保全などによって動作順序が変化してしまう場合がある。これを毎回取得して設計情報を修正するのは大変な手間がかかるが、本工程解析装置1では、信号をモニタすることで装置間の因果関係を自動的に取得できるので、工数をかけることなく製造ラインの現在の状態を取得することができる。   In an actual production line, the operation order may change due to deterioration of the apparatus or improved maintenance. Although it takes a lot of work to correct the design information by acquiring this every time, the process analysis apparatus 1 can automatically acquire the causal relationship between the apparatuses by monitoring the signal. It is possible to obtain the current state of the production line.

〈変形例〉
本実施形態の説明では、信号を発生順序にしたがって並び替えてから、相関係数を算出し、各信号x(3≦j≦N)について、最も相関の高い信号x(1≦k<j)をxの原因としているが、因果構造を作成するための具体的な方法はこれに限られるものではない。
<Modification>
In the description of the present embodiment, the signals are rearranged according to the generation order, the correlation coefficient is calculated, and for each signal x j (3 ≦ j ≦ N), the signal x k with the highest correlation (1 ≦ k < Although j) is the cause of x j , the specific method for creating the causal structure is not limited to this.

因果構造の作成は、各信号について、その対象信号よりも発生順序が前の信号のうち、対象信号との相関が最も高い信号を、その対象信号の原因として定める方法であればどのような方法によって因果構造の作成が行われても良い。例えば、信号sを並び替えることなく、相関係数を求めて、上記条件を満たす信号を探して対象信号の原因としても良い。 The causal structure can be created by any method as long as the signal having the highest correlation with the target signal is determined as the cause of the target signal among the signals whose generation order is earlier than that of the target signal. The causal structure may be created by. For example, the correlation coefficient may be obtained without rearranging the signal s j , and a signal satisfying the above condition may be searched for as a cause of the target signal.

(第2の実施形態)
上記第1の実施形態では、PLC3と製造装置2との間でやりとりされる各信号間の因果構造を獲得し、ユーザ(保全員)に提供することで製造装置3のメンテナンスに利用している。本実施形態では、獲得した因果構造を活用して製造装置3の異常箇所を特定し、ユーザに通知する。
(Second Embodiment)
In the said 1st Embodiment, the causal structure between each signal exchanged between PLC3 and the manufacturing apparatus 2 is acquired, and it uses for the maintenance of the manufacturing apparatus 3 by providing to a user (maintenance worker). . In the present embodiment, the acquired causal structure is utilized to identify an abnormal part of the manufacturing apparatus 3 and notify the user.

図24は、本実施形態に係る工程解析装置1の機能構成を示す図である。本実施形態にかかる工程解析装置1は、大きく分けて因果構造学習部1aと因果構造活用部1bに分かれる。   FIG. 24 is a diagram showing a functional configuration of the process analysis apparatus 1 according to the present embodiment. The process analysis apparatus 1 according to the present embodiment is broadly divided into a causal structure learning unit 1a and a causal structure utilization unit 1b.

因果構造学習部1aは、製造装置2が正常稼働しているときに、PLC3からの信号(学習用信号データ)に基づいて因果構造をあらかじめ作成・学習しておく。この場合に、製造装置2が正常状態にあるかどうかは人間が判断すれば良く、例えば、装置導入時の試運転終了後、メンテナンス終了後などの信号データや、学習用に別途用意された信号データに基づいて因果構造を獲得すれば良い。獲得された因果構造は学習結果記憶部24に格納される。因果構造活用部1bは、この学習結果を通常運用時の信号(診断用データ)と比較して、生産ラインの稼働状況を診断する。   The causal structure learning unit 1a creates and learns a causal structure in advance based on a signal (learning signal data) from the PLC 3 when the manufacturing apparatus 2 is operating normally. In this case, it is only necessary for a human to determine whether or not the manufacturing apparatus 2 is in a normal state. For example, after completion of a trial run at the time of introduction of the apparatus, after completion of maintenance, signal data prepared separately for learning The causal structure may be acquired based on the above. The acquired causal structure is stored in the learning result storage unit 24. The causal structure utilization unit 1b compares the learning result with a signal (diagnosis data) during normal operation, and diagnoses the operation status of the production line.

なお、学習用の信号データは正常状態のデータであることを想定しているので、以降では「学習用のデータ」等のことを「正常時のデータ」等と表現する場合がある。また、稼働状況の診断は通常運用時に行われることを想定しているので、「診断用のデータ」等のことを「通常時のデータ」「通常運用時のデータ」等と表現する場合がある。   Since it is assumed that the learning signal data is data in a normal state, “learning data” or the like may be hereinafter expressed as “normal data” or the like. In addition, since it is assumed that the diagnosis of the operating status is performed during normal operation, “diagnosis data” may be expressed as “normal data”, “normal operation data”, etc. .

〈因果構造学習部〉
因果構造学習部1aは、図24に示すように、信号計測収集部11、特徴量算出部12
、発生順位抽出部13、順位付相関係数算出部14、因果構造獲得部15、及び時間差算出部20を含む。ここで、時間差算出部20以外は、第1の実施形態で説明したものと同様であるため詳細な説明は省略するが、これら各機能部によって因果関係が作成される。
<Causal Structure Learning Department>
As shown in FIG. 24, the causal structure learning unit 1 a includes a signal measurement collection unit 11 and a feature amount calculation unit 12.
, An occurrence rank extraction unit 13, a ranked correlation coefficient calculation unit 14, a causal structure acquisition unit 15, and a time difference calculation unit 20. Here, except for the time difference calculation unit 20, since it is the same as that described in the first embodiment, a detailed description is omitted, but a causal relationship is created by each of these functional units.

時間差算出部20は、因果関係のある信号間の発生時間の差分(時間差)を算出する。以下、時間差算出部20が行う処理について、図25、26を参照して説明する。図25は、時間差算出部20の処理の流れを示すフローチャートである。図26は、時間差算出部20が時間差を算出する差異のデータ間の関係を表す図である。   The time difference calculation unit 20 calculates a difference (time difference) in occurrence time between signals having a causal relationship. Hereinafter, processing performed by the time difference calculation unit 20 will be described with reference to FIGS. FIG. 25 is a flowchart showing a process flow of the time difference calculation unit 20. FIG. 26 is a diagram illustrating a relationship between difference data for which the time difference calculation unit 20 calculates a time difference.

[時間差算出部]
時間差算出部20は、まず、サイクル数を表すインデックスi、信号数を表すインデックスjのそれぞれを1にセットする(S2501)。そして、因果関係リスト107を参照して、信号xの原因となる信号Xを検索する(S2502)。例えば、図26に示すように、因果関係リスト107に「x→x」という因果関係が格納されている場合は、信号xの原因信号Xは信号xである。また、このときの信号xは原因信号に対する結果信号となる。
[Time difference calculator]
First, the time difference calculation unit 20 sets an index i representing the number of cycles and an index j representing the number of signals to 1 (S2501). Then, the causal relationship list 107 is referenced to search for the signal X j that causes the signal x j (S2502). For example, as shown in FIG. 26, if the causal relationship "x 1x 3" is stored in the causal relation list 107, it causes the signal X 3 of the signal x 3 is a signal x 1. The signal x 3 in this case is the result signal for the cause signal.

次に、時間差算出部20は、学習用の信号データから取得した状態変化時刻情報102とその発生順位105とから、サイクルiにおける信号xとその原因信号Xの発生時間の差(x(i)−X(i))を算出する(S2503)。ここで、信号xは接点sの信号を(数4)の式にしたがって発生順に並び替えたものであるので、信号xと接点sとの対応関係を、発生順位リスト105から取得する。図26に示すように、「x→x」の状況を例にとると、発生順位リスト105から、信号xは接点sの信号に対応し、信号xは接点sの信号に対応することが分かる。信号xと接点sとの対応関係が分かると、状態変化時刻情報102からサイクルiにおける各接点の発生時刻の差を求めることで、信号間の時間差を算出することができる。図26の例では、「x→x」の時間差は、接点sと接点sの信号の発生時間の差(s−s)として求めることができ、サイクルi=1では時間差が10ミリ秒、i=2では時間差が11ミリ秒などと求めることができる。 Next, the time difference calculation unit 20 uses the state change time information 102 acquired from the signal data for learning and the occurrence order 105 to determine the difference between the occurrence times of the signal x j and the cause signal X j in the cycle i (x j (I) −X j (i)) is calculated (S2503). Here, the signal x i are those rearranged in chronological order according to the formula a signal contacts s k (number 4), the correspondence between the signal x and the contact s, obtained from the generator rank list 105. As shown in FIG. 26, taking as an example the situation of "x 1x 3", from the generation rank list 105, the signal x 1 corresponds to the signal contacts s 3, the signal x 3 are contacts s 1 of the signal It can be seen that If the correspondence relationship between the signal x and the contact s is known, the time difference between the signals can be calculated by obtaining the difference in the generation time of each contact in the cycle i from the state change time information 102. In the example of FIG. 26, the time difference “x 1 → x 3 ” can be obtained as a difference in signal generation time between the contact s 3 and the contact s 1 (s 1 −s 3 ). Is 10 milliseconds, and when i = 2, the time difference can be obtained as 11 milliseconds.

信号xについて、全サイクルの時間差を求めたか判断し(S2504)、未処理のサイクルがある場合には、次のサイクルについて処理する(S2505)。全サイクルについて処理をした場合は、全信号について時間差を求めたか判断し(S2506)、未処理の信号がある場合には、次の信号について処理する(S2507)。 It is determined whether or not the time difference of all cycles has been obtained for the signal x j (S2504). If there is an unprocessed cycle, the next cycle is processed (S2505). If processing has been performed for all cycles, it is determined whether time differences have been obtained for all signals (S2506). If there is an unprocessed signal, the next signal is processed (S2507).

このようにして、時間差算出部20は、全ての信号、全てのサイクルについて因果関係のある信号間の発生時間の差を求め、時間差データ111を作成する。作成される時間差データ111の例を図26に示している。この算出結果は、学習用の時間差データ111として、因果関係リスト107と共に学習結果記憶部24に格納される。   In this manner, the time difference calculation unit 20 obtains the difference in generation time between signals that are causal for all signals and all cycles, and creates time difference data 111. An example of the created time difference data 111 is shown in FIG. This calculation result is stored in the learning result storage unit 24 together with the causal relationship list 107 as time difference data 111 for learning.

〈因果構造活用部〉
因果構造活用部1bは、通常運用時などにPLC3からの信号を取得して、正常時の学習データと比較することで、製造装置2が正常に稼働しているかどうか診断する。ここでは、因果構造活用部1bが常時あるいは定期的にPLC3からの信号を取得して、稼働状況を診断することを想定して説明するが、通常は診断せずに保全員が診断を指示したときのみ診断を行っても良い。
<Causal structure utilization department>
The causal structure utilization unit 1b diagnoses whether the manufacturing apparatus 2 is operating normally by acquiring a signal from the PLC 3 during normal operation or the like and comparing it with the learning data during normal operation. Here, explanation will be made on the assumption that the causal structure utilization unit 1b obtains a signal from the PLC 3 constantly or periodically and diagnoses the operation status, but the maintenance staff usually instructs the diagnosis without making a diagnosis. Diagnosis may be made only at times.

因果構造活用部1bは、図24に示すように、診断用データ信号計測収集部21、診断用データ時間差算出部22、及び異常箇所推定部23を含む。なお、図24では、因果構造活用部1bの各機能部を因果構造学習部1aとは異なるものとして示しているが、これ
は説明の便宜のためであり、因果構造学習部1aと同一のプログラムによって各機能部が実現されても良い。
As shown in FIG. 24, the causal structure utilization unit 1 b includes a diagnostic data signal measurement collection unit 21, a diagnostic data time difference calculation unit 22, and an abnormal point estimation unit 23. In FIG. 24, each functional unit of the causal structure utilization unit 1b is shown as being different from the causal structure learning unit 1a. However, this is for convenience of explanation, and the same program as the causal structure learning unit 1a. Each functional unit may be realized by the above.

[診断用データ信号計測収集部]
診断用データ信号計測収集部21は、通常運用時にPLC3から信号情報を取得し、診断用信号情報201を生成する。診断用データ信号計測収集部21が行う処理の内容をフローチャートを図27に示す。この処理の内容は、因果構造学習部1aの信号計測収集部11の処理(図2のフローチャート)と同様であるため、説明は省略する。なお、診断用として取得する信号は、複数サイクルの間にわたって取得しても良く、1サイクルのみを取得しても良い。本実施形態では、診断用の信号を複数サイクルにわたって取得する。
[Data signal measurement and collection unit for diagnosis]
The diagnostic data signal measurement collection unit 21 acquires signal information from the PLC 3 during normal operation, and generates diagnostic signal information 201. FIG. 27 is a flowchart showing the contents of the process performed by the diagnostic data signal measurement / collection unit 21. Since the content of this process is the same as the process (flowchart in FIG. 2) of the signal measurement collection unit 11 of the causal structure learning unit 1a, the description is omitted. The signal acquired for diagnosis may be acquired over a plurality of cycles, or only one cycle may be acquired. In this embodiment, a diagnostic signal is acquired over a plurality of cycles.

[診断用データ時間差算出部]
次に診断用データの時間差算出部22の処理について、図28を参照して説明する。時間差算出部22は、診断用の信号情報201から、各信号の状態変化の時刻を抽出して、診断用の状態変化時刻情報202を作成する(ステップS2801〜S2804)。この処理は、因果構造学習部1aの特徴量算出部12の処理の一部(図7のフローチャートのステップS701〜S704)と同じである。そして、作成された診断用の状態変化時刻情報202から、因果関係のある信号間の時間差を算出し、診断用時間差データ203を作成する(ステップS2805〜S2811)。この処理は上記時間差算出部20の時間差算出処理(図25)と同様である。
[Diagnostic data time difference calculator]
Next, the processing of the time difference calculation unit 22 for diagnostic data will be described with reference to FIG. The time difference calculation unit 22 extracts the state change time information of each signal from the signal information 201 for diagnosis and creates the state change time information 202 for diagnosis (steps S2801 to S2804). This process is the same as part of the process of the feature quantity calculation unit 12 of the causal structure learning unit 1a (steps S701 to S704 in the flowchart of FIG. 7). Then, the time difference between the causal signals is calculated from the created diagnostic state change time information 202, and diagnostic time difference data 203 is created (steps S2805 to S2811). This process is the same as the time difference calculation process (FIG. 25) of the time difference calculation unit 20.

[異常箇所推定部]
異常箇所推定部23は、学習用(正常時)の時間差データ111と、診断用(通常運用時)の時間差データ203を比較して、異常箇所の候補を求める。ここでは、正常時と通常運用時とで、信号間の時間差のずれが大きい箇所を異常のある箇所として推定する。
[Abnormal part estimation part]
The abnormal part estimation unit 23 compares the time difference data 111 for learning (during normal operation) with the time difference data 203 for diagnosis (during normal operation) to obtain abnormal part candidates. Here, a portion where the time difference between the signals is large between normal operation and normal operation is estimated as a portion having an abnormality.

以下、異常箇所推定部23の処理について、図29,30を参照して説明する。図29は異常箇所推定部23が行う異常箇所推定処理の流れを示すフローチャートであり、図30(a)は作成される異常候補リスト204の例、図30(b)は時間差情報付きの因果関係リストの模式図を示す。   Hereinafter, the process of the abnormal part estimation part 23 is demonstrated with reference to FIG. FIG. 29 is a flowchart showing the flow of an abnormal location estimation process performed by the abnormal location estimation unit 23. FIG. 30 (a) is an example of an abnormal candidate list 204 to be created, and FIG. 30 (b) is a causal relationship with time difference information. A schematic diagram of the list is shown.

異常箇所推定部23は、まず、診断用の状態変化時刻情報201から生産サイクルを算出する(S2901)。この生産サイクルは、診断用として取得した信号のうち、マシンサイクルを表す信号のサイクルタイムを計算することによって算出できる。マシンサイクルを表す信号は、因果構造学習部1aの特徴量算出部12で基準信号として用いられた信号であって、診断用データとして再度計測収集されたものを用いる。サイクルタイムは、この基準信号の立ち上がりから次の立ち上がりまでの時間差として表される。   The abnormal part estimation unit 23 first calculates a production cycle from the state change time information 201 for diagnosis (S2901). This production cycle can be calculated by calculating the cycle time of a signal representing a machine cycle among signals acquired for diagnosis. The signal representing the machine cycle is a signal that is used as a reference signal in the feature amount calculation unit 12 of the causal structure learning unit 1a and that is remeasured and collected as diagnostic data. The cycle time is expressed as a time difference from the rising edge of this reference signal to the next rising edge.

次に、異常箇所推定部23は生産サイクルのしきい値を読み込む(S2902)。生産サイクルのしきい値はあらかじめ定められた設定値であり、ユーザが設定可能である。生産サイクルしきい値の設定画面の例を図32に示す。図32に示すようにディスプレイなどの表示装置上のしきい値設定画面で、ユーザはキーボードやマウスなどの入力装置を使って、しきい値の設定を行う。この例では、正常と判断される生産サイクルの上限値を入力する構成となっている。例えば、予定される生産サイクルが30秒の生産ラインで、32秒までは正常であると判断する場合に、図32の設定画面で32秒と入力する。ただし、正常な生産サイクル(例えば、30秒)と許容する遅れの上限値(例えば、2秒)を入力させるようにしても良い。また、表示装置と入力装置を介してしきい値を設定するのではなく、設定ファイルなどを読み込ませてしきい値を設定しても構わない。   Next, the abnormal part estimation part 23 reads the threshold value of a production cycle (S2902). The production cycle threshold value is a predetermined set value and can be set by the user. An example of a production cycle threshold setting screen is shown in FIG. As shown in FIG. 32, on a threshold setting screen on a display device such as a display, the user sets a threshold value using an input device such as a keyboard or a mouse. In this example, an upper limit value of a production cycle that is determined to be normal is input. For example, when it is determined that the production line scheduled for 30 seconds is normal up to 32 seconds, 32 seconds is input on the setting screen of FIG. However, a normal production cycle (for example, 30 seconds) and an allowable upper limit value (for example, 2 seconds) may be input. Further, instead of setting the threshold value via the display device and the input device, the threshold value may be set by reading a setting file or the like.

異常箇所推定部23は、生産サイクルのしきい値を読み込んだ後、ステップS2901
で算出した生産サイクルがしきい値を超えているかを判断する(S2903)。サイクルオーバー(タクト遅れ)が発生していない場合(S2903−NO)は、異常が発生していないと判断して、異常箇所の推定処理は行わない。一方、サイクルオーバーが発生している場合(S2903−YES)は、異常箇所の推定を行う。
The abnormal part estimation unit 23 reads the production cycle threshold value, and then performs step S2901.
It is determined whether the production cycle calculated in (1) exceeds a threshold value (S2903). If cycle over (tact delay) has not occurred (S2903-NO), it is determined that no abnormality has occurred, and the abnormal point estimation process is not performed. On the other hand, when the cycle over has occurred (S2903-YES), the abnormal part is estimated.

ステップS2904で、学習用の時間差データ111と診断用の時間差データ203から、学習用及び診断用のそれぞれについて因果関係のある信号間の時間差の平均を求め、これら2つの平均値の差分を算出する。図26に示すような時間差データには、因果関係のある信号間(x→xなど)の時間差がサイクル毎に格納されている。この平均をとることで、時間差データの平均値を求めることができる。 In step S2904, from the learning time difference data 111 and the diagnostic time difference data 203, an average of the time differences between the causal signals for the learning and the diagnostic is obtained, and the difference between these two average values is calculated. . In the time difference data as shown in FIG. 26, the time difference between causal signals (such as x 1 → x 3 ) is stored for each cycle. By taking this average, the average value of the time difference data can be obtained.

なお、図29のフローチャートでは、学習用データに関して時間差データの平均値を異常箇所推定処理毎に求めているが、学習用データについての平均値は変化しない値であるので、一度だけ求めて記憶しておいても良い。   In the flowchart of FIG. 29, the average value of the time difference data for the learning data is obtained for each abnormal part estimation process. However, since the average value for the learning data is a value that does not change, it is obtained and stored only once. You can keep it.

各信号間の時間差が図30(b)に示すような状態であるとする。この場合、信号xと信号xの時間差の平均は、学習用データでは10ミリ秒、診断用データでは150ミリ秒であり、その差分は140ミリ秒である。同様に信号xと信号xの時間差の平均は、学習用データで30ミリ秒、診断用データで35ミリ秒で、その差分は5ミリ秒である。このように、因果関係のある信号間のそれぞれについて、時間差データ平均値の学習時と診断時とのズレを算出する。算出した時間差データ平均値のズレは、異常候補リスト204に格納される。 Assume that the time difference between the signals is as shown in FIG. In this case, the average time difference between the signals x 1 and the signal x 3 is 10 milliseconds in the learning data, is 150 ms in the diagnostic data, the difference is 140 ms. Similarly the average time difference between the signals x 1 and the signal x 2 is 30 ms in the learning data, at 35 milliseconds diagnosis data, the difference is 5 ms. In this way, for each of the causal signals, the difference between the time difference data average value learning time and the diagnosis time is calculated. The deviation of the calculated time difference data average value is stored in the abnormality candidate list 204.

異常候補リスト204内のデータは、時間差の差分(ズレ)が大きい順にソートされる(S2905)。つまり、異常候補リスト204に、図30(a)に示すように時間差のズレが大きい順に、異常候補順位が付加される。なお、ここでいうソート処理は時間差のズレの大きさの順番を決定する処理を指し、必ずしもメモリ上でデータの移動(入れ替え)を伴う必要はない。   The data in the abnormality candidate list 204 is sorted in descending order of time difference (displacement) (S2905). That is, abnormality candidate ranks are added to the abnormality candidate list 204 in descending order of time difference as shown in FIG. Note that the sort processing here refers to processing for determining the order of time difference deviation, and does not necessarily involve data movement (replacement) on the memory.

〈表示部〉
表示部17は、学習用信号データから作成された因果関係リストを表示する。本実施形態においては、生産ラインに異常があるときは、異常箇所を強調して表示を行う。
<Display section>
The display unit 17 displays a causal relationship list created from the learning signal data. In the present embodiment, when there is an abnormality in the production line, the abnormal part is highlighted and displayed.

表示部17が行う表示処理について、図31のフローチャートを参照して説明する。まず、表示部17は、因果関係リスト107から信号間の因果関係(ツリー構造)を取得し、設計情報110を参照して変数名に各信号に対応する接点名を取得する(S3101)。そして、信号の接点名を利用して、因果構造をグラフ表示する。   The display process performed by the display unit 17 will be described with reference to the flowchart of FIG. First, the display unit 17 acquires a causal relationship (tree structure) between signals from the causal relationship list 107, refers to the design information 110, and acquires a contact name corresponding to each signal as a variable name (S3101). Then, the causal structure is displayed in a graph using the contact name of the signal.

次に、表示部17は異常候補リストが作成されているか、すなわち、生産サイクルに異常があるか判定する(S3103)。異常候補リスト204が作成されていない場合には、これ以上の表示を行わない。したがって、異常候補リスト204が作成されていない場合の、出力画面は図34のようになる。図34は、出力画面と、この出力画面を作成するために使われるデータが示されている。異常がない場合、出力画面には因果構造が設計情報(信号の名称)で置き換えられたツリー構造が表示される。   Next, the display unit 17 determines whether an abnormality candidate list has been created, that is, whether there is an abnormality in the production cycle (S3103). If the abnormality candidate list 204 has not been created, no further display is performed. Accordingly, the output screen when the abnormality candidate list 204 is not created is as shown in FIG. FIG. 34 shows an output screen and data used to create the output screen. When there is no abnormality, a tree structure in which the causal structure is replaced with design information (signal name) is displayed on the output screen.

一方、異常候補リストが作成されている、すなわち、生産サイクルに異常がある(しきい値オーバーしている)場合は(S3103−YES)は、異常候補リスト204、正常時(学習用)時間差データ111、通常時(診断用)時間差データ203、異常候補表示数205、及び異常候補表示しきい値206を参照して、すでに表示している因果構造のうち、異常箇所を強調表示する(S3104)。   On the other hand, if an abnormality candidate list has been created, that is, if there is an abnormality in the production cycle (threshold is exceeded) (S3103—YES), the abnormality candidate list 204, normal (for learning) time difference data 111, with reference to the normal (diagnostic) time difference data 203, the abnormality candidate display number 205, and the abnormality candidate display threshold 206, an abnormal part is highlighted among the already displayed causal structures (S3104). .

異常推定箇所が強調表示された出力画面は図35のようになる。図35は、出力画面と、この出力画面を作成するために使われるデータが示されている。この例では、2箇所が異常箇所として強調表示されている。   FIG. 35 shows an output screen on which the estimated abnormality location is highlighted. FIG. 35 shows an output screen and data used to create the output screen. In this example, two places are highlighted as abnormal places.

ここで、異常候補表示数205と異常候補表示しきい値206の設定に応じて、異常候補のうちのいずれを強調表示するかが決定される。なお、異常候補表示数205及び異常候補表示しきい値206は、図33に示すような、設定画面を通じて入力する。あるいは、設定ファイルを読み込むように構成しても良い。異常候補表示数205は、何カ所を異常候補として表示するかを決定する設定値である。異常候補表示しきい値206は、時間差のずれがどの程度大きい場合を異常として扱うかを決定する設定値である。   Here, depending on the settings of the abnormality candidate display number 205 and the abnormality candidate display threshold 206, it is determined which of the abnormality candidates is highlighted. The abnormality candidate display number 205 and the abnormality candidate display threshold 206 are input through a setting screen as shown in FIG. Or you may comprise so that a setting file may be read. The abnormality candidate display number 205 is a setting value that determines how many places are displayed as abnormality candidates. The abnormality candidate display threshold 206 is a setting value that determines how large a time difference shift is to be treated as an abnormality.

異常候補表示数205と異常候補表示しきい値206は、いずれか一方のみが設定されていればよい。異常候補表示数205のみが設定されている場合には、異常候補リスト204の上位から順に、設定値の数だけ強調表示される。異常候補表示しきい値206のみが設定されている場合には、異常候補リスト204のうち時間差のずれが設定値以上である箇所が、強調表示される。異常候補表示数205と異常候補表示しきい値の206の両方が設定されている場合には、異常候補リスト204のうち、時間差のずれが設定値以上の箇所のうち、ずれが大きいものから設定の数だけ強調表示される。図33の例では異常候補表示数のみが入力されており、この場合は図35に示すように時間差のずれが大きい2箇所が強調表示される。   Only one of the abnormality candidate display number 205 and the abnormality candidate display threshold 206 may be set. When only the abnormality candidate display number 205 is set, the number of set values is highlighted in order from the top of the abnormality candidate list 204. When only the abnormality candidate display threshold 206 is set, a portion of the abnormality candidate list 204 where the time difference shift is equal to or larger than the set value is highlighted. When both of the abnormality candidate display number 205 and the abnormality candidate display threshold 206 are set, from the abnormality candidate list 204, the part having the larger deviation among the positions where the deviation of the time difference is equal to or larger than the set value is set. Will be highlighted. In the example of FIG. 33, only the number of abnormal candidate displays is input, and in this case, as shown in FIG. 35, two locations with a large time difference are highlighted.

図36では、信号x(シリンダ2出駆動)と信号x(シリンダ3出駆動)との間(x→x)で正常時との時間差のズレが大きく、この箇所が異常であると推定されている。そして、ここでは、結果信号の方が強調表示されている。しかしながら、異常箇所の強調表示の方法は、種々の方法を採用可能である。図35に強調表示の例を示す。 In FIG. 36, there is a large time difference between the signal x 2 (cylinder 2 output drive) and the signal x 4 (cylinder 3 output drive) (x 2 → x 4 ) from the normal state, and this point is abnormal. It is estimated that. Here, the result signal is highlighted. However, various methods can be adopted as a method for highlighting abnormal portions. FIG. 35 shows an example of highlighting.

図35(a)は、図34に示すように結果信号の部分を強調表示する方法である。図35(b)は、原因信号の部分を強調表示する方法である。図35(c)は、原因信号と結果信号の遷移部分を強調表示する方法である。図36(d)は、上記の全てを強調表示する方法である。強調表示は、ユーザが異常箇所を出力画面上で認識しやすいような形式であれば、どのような方法によって行っても良い。   FIG. 35A shows a method of highlighting the result signal portion as shown in FIG. FIG. 35B shows a method for highlighting the cause signal portion. FIG. 35C shows a method of highlighting the transition portion between the cause signal and the result signal. FIG. 36D shows a method for highlighting all of the above. The highlighting may be performed by any method as long as the user can easily recognize the abnormal part on the output screen.

〈実施形態の作用・効果〉
本実施形態に係る工程解析装置1によれば、製造装置2とPLC3の間でやりとりされる信号の時系列情報から因果関係を作成するとともに、因果関係のある信号間の発生時間の差を取得できる。そして、通常運用時に、因果関係のある信号間の時間差をモニタリングできる。異常が発生した場合には、この信号間の時間差のずれが大きい箇所を異常箇所として推定して、ユーザの認識しやすい形式で画面表示を行うことができる。
<Operation and effect of the embodiment>
According to the process analysis apparatus 1 according to the present embodiment, a causal relationship is created from time-series information of signals exchanged between the manufacturing apparatus 2 and the PLC 3, and a difference in generation time between signals having a causal relationship is acquired. it can. During normal operation, the time difference between the causal signals can be monitored. When an abnormality occurs, it is possible to estimate a portion where the time difference between the signals is large as an abnormal portion and display the screen in a format that is easy for the user to recognize.

このように本実施形態に係る工程解析装置1では、異常が発生した場合に、異常箇所が自動的に特定されユーザに対して表示されるので、速やかな対処・修理が可能となる。タクト遅れが発生すると、それだけ製造効率が下がってしまう。迅速に対処可能となることで損失を最小限にすることができる。   As described above, in the process analysis apparatus 1 according to the present embodiment, when an abnormality occurs, the abnormal part is automatically identified and displayed to the user, so that prompt handling and repair are possible. If the tact delay occurs, the manufacturing efficiency is lowered accordingly. Loss can be minimized by being able to deal with it quickly.

第1の実施形態に係る工程解析装置の機能ブロックを示す図である。It is a figure which shows the functional block of the process analysis apparatus which concerns on 1st Embodiment. 信号計測収集部11が行う処理の流れを示すフローチャートである。It is a flowchart which shows the flow of the process which the signal measurement collection part 11 performs. 信号計測収集部11によって収集される信号の例を示す図である。It is a figure which shows the example of the signal collected by the signal measurement collection part. 信号計測収集部11によって記憶される信号情報101の例を示す図である。It is a figure which shows the example of the signal information 101 memorize | stored by the signal measurement collection part. 信号のサイクルを説明する図である。It is a figure explaining the cycle of a signal. 特徴量算出部12が算出する信号の特徴量を説明する図である。It is a figure explaining the feature-value of the signal which the feature-value calculation part 12 calculates. 特徴量算出部12が行う処理の流れを示すフローチャートである。It is a flowchart which shows the flow of the process which the feature-value calculation part 12 performs. 特徴量算出部12が作成する状態時刻変化情報102の例を示す図である。It is a figure which shows the example of the state time change information 102 which the feature-value calculation part 12 produces. 特徴量算出部12が作成する特徴量データ103の例を示す図である。It is a figure which shows the example of the feature-value data 103 which the feature-value calculation part 12 produces. 発生順位抽出部13が行う処理の流れを示すフローチャートである。It is a flowchart which shows the flow of the process which the generation | occurrence | production order extraction part 13 performs. 発生順位抽出部13が抽出する各サイクルにおける発生順位を説明する図である。It is a figure explaining the generation | occurrence | production order in each cycle which the generation | occurrence | production order extraction part 13 extracts. 発生順位抽出部13が作成する発生順位105の例を示す図である。It is a figure which shows the example of the generation | occurrence | production order 105 which the generation | occurrence | production order extraction part 13 produces. 順位付相関係数算出部14が行う処理の流れを示すフローチャートである。It is a flowchart which shows the flow of the process which the correlation coefficient calculation part 14 with a ranking performs. 順位付相関係数算出部14が行う相関係数算出処理を説明する図である。It is a figure explaining the correlation coefficient calculation process which the correlation coefficient calculation part 14 with a ranking performs. 因果構造獲得部15が行う処理の流れを示すフローチャートである。It is a flowchart which shows the flow of the process which the causal structure acquisition part 15 performs. 因果構造獲得部15が行う処理を説明する図である。It is a figure explaining the process which the causal structure acquisition part 15 performs. 表示部17に表示される出力画面を表す図である。FIG. 10 is a diagram illustrating an output screen displayed on the display unit 17. 表示部17が行う処理の流れを示すフローチャートである。It is a flowchart which shows the flow of the process which the display part 17 performs. 表示部17が行う処理を説明する図である。It is a figure explaining the process which the display part 17 performs. 表示部17による因果構造のグラフ表示で、因果構造全体を表示した例を示す図である。It is a figure which shows the example which displayed the whole causal structure by the graph display of the causal structure by the display part. 表示部17による因果構造のグラフ表示で、因果構造の一部の経路のみを表示した例を示す図である。It is a figure which shows the example which displayed only the one part path | route of the causal structure by the graph display of the causal structure by the display part 17. FIG. 表示部17による因果構造のグラフ表示で、因果構造の一部の経路のみを表示した例を示す図である。It is a figure which shows the example which displayed only the one part path | route of the causal structure by the graph display of the causal structure by the display part 17. FIG. 表示部17による因果構造のグラフ表示で、因果構造の一部の経路のみを表示した例を示す図である。It is a figure which shows the example which displayed only the one part path | route of the causal structure by the graph display of the causal structure by the display part 17. FIG. 第2の実施形態に係る工程解析装置の機能ブロックを示す図である。It is a figure which shows the functional block of the process analysis apparatus which concerns on 2nd Embodiment. 時間差算出部20が行う処理の流れを示すフローチャートである。It is a flowchart which shows the flow of the process which the time difference calculation part 20 performs. 時間差データ111の作成するためのデータの流れ、および時間差データ111のデータ構造を説明する図である。3 is a diagram for explaining a data flow for creating time difference data 111 and a data structure of time difference data 111. FIG. 診断用データ信号計測収集部21が行う処理の流れを示すフローチャートである。It is a flowchart which shows the flow of the process which the data signal measurement collection part 21 for a diagnosis performs. 診断用データ時間差算出部22が行う処理の流れを示すフローチャートである。It is a flowchart which shows the flow of the process which the diagnostic data time difference calculation part 22 performs. 異常箇所推定部23が行う処理の流れを示すフローチャートである。It is a flowchart which shows the flow of the process which the abnormal location estimation part 23 performs. 異常箇所推定部23が作成する異常候補リスト204の例を示す図である。It is a figure which shows the example of the abnormality candidate list 204 which the abnormal location estimation part 23 produces. 表示部17が行う処理の流れを示すフローチャートである。It is a flowchart which shows the flow of the process which the display part 17 performs. 生産サイクルの異常判定用のしきい値を入力する設定画面の例を示す図である。It is a figure which shows the example of the setting screen which inputs the threshold value for abnormality determination of a production cycle. 異常箇所の表示数および異常とみなすしきい値を入力する設定画面の例を示す図である。It is a figure which shows the example of the setting screen which inputs the display number of an abnormal location, and the threshold value considered as abnormality. 表示部17による因果構造のグラフ表示で、製造装置に異常がなく異常箇所推定が行われていない場合の例を示す図である。It is a figure which shows the example in case there is no abnormality in a manufacturing apparatus and the abnormal location estimation is not performed by the graph display of the causal structure by the display part 17. FIG. 表示部17による因果構造のグラフ表示で、製造装置に異常があり異常推定箇所が強調表示されている場合の例を示す図である。It is a figure which shows the example in case a manufacturing apparatus has abnormality and the abnormal estimated part is highlighted by the graph display of the causal structure by the display part 17. FIG. 表示部17による因果構造のグラフ表示において、異常推定箇所を強調表示する方法の例を示す図である。It is a figure which shows the example of the method of highlighting an abnormality estimated location in the graph display of the causal structure by the display part.

符号の説明Explanation of symbols

1 工程解析装置
2 製造装置
3 プログラマブルロジックコントローラ
11 信号計測収集部
12 特徴量算出部
13 発生順位抽出部
14 順位付相関係数算出部
15 因果構造獲得部
17 表示部
20 時間差算出部
21 診断用データ信号計測収集部
22 診断用データ時間差算出部
23 異常箇所推定部
24 学習結果記憶部
101 信号情報
102 状態変化時刻情報
103 特徴量データ
104 サイクルごとの発生順位
105 発生順位
106 相関係数
107 因果関係リスト
110 設計情報
111 学習用時間差データ
201 診断用信号情報
202 診断用状態変化時刻情報
203 診断用時間差データ
204 異常候補リスト
DESCRIPTION OF SYMBOLS 1 Process analysis apparatus 2 Manufacturing apparatus 3 Programmable logic controller 11 Signal measurement collection part 12 Feature-value calculation part 13 Generation | occurrence | production order extraction part 14 Ranking correlation coefficient calculation part 15 Causal structure acquisition part 17 Display part 20 Time difference calculation part 21 Diagnosis data Signal measurement collection unit 22 Diagnostic data time difference calculation unit 23 Abnormal part estimation unit 24 Learning result storage unit 101 Signal information 102 State change time information 103 Feature data 104 Generation order 105 per cycle Generation order 106 Correlation coefficient 107 Causal relation list 110 Design information 111 Time difference data for learning 201 Signal information for diagnosis 202 State change time information for diagnosis 203 Time difference data for diagnosis 204 Abnormality candidate list

Claims (19)

製造装置と制御装置との間で略周期的に入出力される複数の信号の関係を解析する工程解析装置であって、
前記制御装置から送信される制御信号と、該制御信号に対応する応答信号とを周期ごとに取得する信号取得手段と、
前記信号取得手段によって取得された各々の信号の、前記周期ごとの特徴量を算出する特徴量算出手段と、
前記信号取得手段によって取得された信号の、前記周期ごとの発生順序を抽出する発生順序抽出手段と、
特徴量間の相関と信号の発生順序とに基づいて、信号間の因果構造を作成する因果構造作成手段と、
を備える工程解析装置。
A process analysis device that analyzes the relationship between a plurality of signals that are input and output substantially periodically between a manufacturing device and a control device ,
Signal acquisition means for acquiring a control signal transmitted from the control device and a response signal corresponding to the control signal for each period ;
A feature amount calculating means for calculating a feature amount for each period of each signal acquired by the signal acquiring means;
Generation order extraction means for extracting the generation order for each period of the signals acquired by the signal acquisition means;
A causal structure creation means for creating a causal structure between signals based on a correlation between feature quantities and a signal generation order;
A process analysis apparatus comprising:
前記因果構造作成手段は、各信号について、該信号よりも発生順序が前の信号のうち該信号との相関が最も高い信号を、該信号の原因として因果構造を作成する
ことを特徴とする請求項1に記載の工程解析装置。
The causal structure creating means creates a causal structure for each signal by using a signal having the highest correlation with the signal among the signals whose generation order is earlier than the signal as a cause of the signal. Item 2. The process analysis apparatus according to Item 1.
各信号の発生順序は、各周期における発生順序の平均に基づいて算出される
ことを特徴とする請求項1または2に記載の工程解析装置。
The process analysis apparatus according to claim 1, wherein the generation order of each signal is calculated based on an average of the generation order in each cycle.
各信号の特徴量は、該信号と所定の基準信号との立ち上がり時刻の差として算出される
ことを特徴とする請求項1〜3のいずれかに記載の工程解析装置。
The process analysis apparatus according to claim 1, wherein the feature amount of each signal is calculated as a difference in rise time between the signal and a predetermined reference signal.
作成された因果構造を表示する表示手段を備える
ことを特徴とする請求項1〜4のいずれかに記載の工程解析装置。
The process analysis apparatus according to claim 1, further comprising display means for displaying the created causal structure.
前記表示手段は、信号に対応する名称を設計情報から取得し、該名称を用いて因果構造を表示する
ことを特徴とする請求項5に記載の工程解析装置。
The process analysis apparatus according to claim 5, wherein the display unit acquires a name corresponding to the signal from the design information, and displays the causal structure using the name.
製造装置と制御装置との間で略周期的に入出力される複数の信号の関係を解析する工程解析方法であって、
情報処理装置が、
前記制御装置から送信される制御信号と、該制御信号に対応する応答信号とを周期ごとに取得するステップと、
取得された各々の信号の、前記周期ごとの特徴量を算出するステップと、
取得された信号の、前記周期ごとの発生順序を抽出するステップと、
特徴量間の相関と信号の発生順序とに基づいて、信号間の因果構造を作成するステップと、
を実行する工程解析方法。
A process analysis method for analyzing a relationship between a plurality of signals input and output substantially periodically between a manufacturing apparatus and a control apparatus ,
Information processing device
Obtaining a control signal transmitted from the control device and a response signal corresponding to the control signal for each period ;
Calculating a feature amount for each period of each acquired signal ;
Extracting the order of occurrence of the acquired signals for each period ;
Creating a causal structure between signals based on the correlation between feature quantities and the order of signal generation;
Process analysis method to execute.
前記因果構造を作成するステップでは、各信号について、該信号よりも発生順序が前の信号のうち該信号との相関が最も高い信号を、該信号の原因として因果構造を作成する
ことを特徴とする請求項7に記載の工程解析方法。
In the step of creating the causal structure, for each signal, a causal structure is created by using a signal having the highest correlation with the signal among the signals whose generation order is earlier than the signal as a cause of the signal. The process analysis method according to claim 7.
各信号の発生順序は、各周期における発生順序の平均に基づいて算出される
ことを特徴とする請求項7または8に記載の工程解析方法。
The process analysis method according to claim 7 or 8, wherein the generation order of each signal is calculated based on an average of the generation order in each cycle.
各信号の特徴量は、該信号と所定の基準信号との立ち上がり時刻の差として算出される
ことを特徴とする請求項7〜9のいずれかに記載の工程解析方法。
The process analysis method according to claim 7, wherein the feature amount of each signal is calculated as a difference in rise time between the signal and a predetermined reference signal.
作成された因果構造を表示するステップをさらに実行する
ことを特徴とする請求項7〜10のいずれかに記載の工程解析方法。
The process analysis method according to claim 7, further comprising a step of displaying the created causal structure.
前記因果構造を表示するステップでは、信号に対応する名称を設計情報から取得し、該名称を用いて因果構造を表示する
ことを特徴とする請求項11に記載の工程解析方法。
The process analysis method according to claim 11, wherein in the step of displaying the causal structure, a name corresponding to the signal is acquired from the design information, and the causal structure is displayed using the name.
製造装置と制御装置との間で略周期的に入出力される複数の信号の関係を解析する工程解析プログラムであって、
情報処理装置に、
前記制御装置から送信される制御信号と、該制御信号に対応する応答信号とを周期ごとに取得するステップと、
取得された各々の信号の、前記周期ごとの特徴量を算出するステップと、
取得された信号の、前記周期ごとの発生順序を抽出するステップと、
特徴量間の相関と信号の発生順序とに基づいて、信号間の因果構造を作成するステップと、
を実行させる工程解析プログラム。
A process analysis program for analyzing a relationship between a plurality of signals that are input and output substantially periodically between a manufacturing apparatus and a control apparatus ,
In the information processing device,
Obtaining a control signal transmitted from the control device and a response signal corresponding to the control signal for each period ;
Calculating a feature amount for each period of each acquired signal ;
Extracting the order of occurrence of the acquired signals for each period ;
Creating a causal structure between signals based on the correlation between feature quantities and the order of signal generation;
Process analysis program to execute
学習用の信号データを利用して、あらかじめ求められた各信号間の因果構造及び因果関係のある信号間の時間差を、学習結果として記憶する記憶手段と、
診断用の信号データを利用して、因果関係のある信号間の時間差を算出し、学習結果と比較して診断用データにおける異常箇所を推定する異常箇所推定手段と、を有する
ことを特徴とする請求項1〜6のいずれかに記載の工程解析装置。
A storage means for storing a causal structure between signals obtained in advance and a time difference between signals having a causal relationship as learning results using signal data for learning;
An abnormality location estimating means for calculating a time difference between causal signals using diagnostic signal data and estimating an abnormal location in the diagnostic data by comparing with a learning result. The process analysis apparatus according to claim 1.
前記異常箇所推定手段は、診断用の信号データにおける因果関係のある信号間の時間差と、学習結果における時間差とのずれが所定のしきい値よりも大きい箇所を異常箇所であると推定する
ことを特徴とする請求項14に記載の工程解析装置。
The abnormal location estimating means estimates that a location where a difference between a time difference between causal signals in the diagnostic signal data and a time difference in the learning result is larger than a predetermined threshold is an abnormal location. The process analysis apparatus according to claim 14, wherein the process analysis apparatus is characterized in that:
信号間の因果構造を、推定された異常箇所を強調しつつ表示する表示部を有する
ことを特徴とする請求項14又は15に記載の工程解析装置。
The process analysis apparatus according to claim 14, further comprising: a display unit that displays the causal structure between the signals while highlighting the estimated abnormal part.
前記情報処理装置が、さらに、
学習用の信号データを利用して、信号間の因果構造及び因果関係のある信号間の時間差を求めて記憶するステップと、
診断用の信号データを利用して、因果関係のある信号間の時間差を算出し、学習結果と比較して診断用データにおける異常箇所を推定するステップと、
を実行することを特徴とする請求項7〜11のいずれかに記載の工程解析方法。
The information processing apparatus further includes:
Using the learning signal data, determining and storing a causal structure between signals and a time difference between signals having a causal relationship;
Using the diagnostic signal data, calculating a time difference between the causal signals, and comparing the learning result with estimating the abnormal portion in the diagnostic data;
The process analysis method according to claim 7, wherein the process analysis method is executed.
前記異常箇所を推定するステップでは、診断用の信号データにおける因果関係のある信号間の時間差と、学習結果における時間差とのずれが所定のしきい値よりも大きい箇所を異常箇所であると推定する
ことを特徴とする請求項17に記載の工程解析方法。
In the step of estimating the abnormal location, a location where the difference between the time difference between the causal signals in the signal data for diagnosis and the time difference in the learning result is larger than a predetermined threshold is estimated as an abnormal location. The process analysis method according to claim 17.
前記情報処理装置は、信号間の因果構造を、推定された異常箇所を強調しつつ表示するステップをさらに実行する
ことを特徴とする請求項17又は18に記載の工程解析方法。
The process analysis method according to claim 17 or 18, wherein the information processing apparatus further executes a step of displaying a causal structure between signals while emphasizing an estimated abnormal part.
JP2008049807A 2007-08-10 2008-02-29 Process analysis apparatus, process analysis method, and process analysis program Expired - Fee Related JP5151556B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2008049807A JP5151556B2 (en) 2007-08-10 2008-02-29 Process analysis apparatus, process analysis method, and process analysis program

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2007209872 2007-08-10
JP2007209872 2007-08-10
JP2008049807A JP5151556B2 (en) 2007-08-10 2008-02-29 Process analysis apparatus, process analysis method, and process analysis program

Publications (2)

Publication Number Publication Date
JP2009064407A JP2009064407A (en) 2009-03-26
JP5151556B2 true JP5151556B2 (en) 2013-02-27

Family

ID=40558920

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008049807A Expired - Fee Related JP5151556B2 (en) 2007-08-10 2008-02-29 Process analysis apparatus, process analysis method, and process analysis program

Country Status (1)

Country Link
JP (1) JP5151556B2 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110275489A (en) * 2018-03-13 2019-09-24 发那科株式会社 Data time series analysis device
US11188064B1 (en) 2021-05-04 2021-11-30 Ixden Ltd. Process flow abnormality detection system and method
WO2023233927A1 (en) * 2022-06-03 2023-12-07 オムロン株式会社 Abnormality detection device, abnormality detection method, and program

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5081999B1 (en) * 2011-06-22 2012-11-28 株式会社日立エンジニアリング・アンド・サービス How to display abnormal sign diagnosis results
JP5081998B1 (en) * 2011-06-22 2012-11-28 株式会社日立エンジニアリング・アンド・サービス Abnormal sign diagnostic apparatus and abnormal sign diagnostic method
CN108983710B (en) * 2017-06-02 2021-05-18 欧姆龙株式会社 Process analysis device, process analysis method, and storage medium
JP6551565B2 (en) * 2017-06-02 2019-07-31 オムロン株式会社 Process analysis apparatus, process analysis method, and process analysis program
JP7137810B2 (en) * 2017-11-08 2022-09-15 i Smart Technologies株式会社 Production control equipment and production control system
JP6798968B2 (en) 2017-11-22 2020-12-09 ファナック株式会社 Noise cause estimation device
JP7382882B2 (en) 2020-03-31 2023-11-17 日鉄ソリューションズ株式会社 Information processing device, information processing method and program
JP2022029983A (en) * 2020-08-06 2022-02-18 オムロン株式会社 Display system, display method, and display program
JP2022137914A (en) 2021-03-09 2022-09-22 オムロン株式会社 Information processing device, information processing method and information processing program
JP2022156257A (en) * 2021-03-31 2022-10-14 オムロン株式会社 Information processing device, information processing method and information processing program
JP2023092184A (en) 2021-12-21 2023-07-03 オムロン株式会社 Information processor and information processing method and information processing program

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4495960B2 (en) * 2003-12-26 2010-07-07 キヤノンItソリューションズ株式会社 Model creation device for the relationship between process and quality
JP4314123B2 (en) * 2004-01-30 2009-08-12 株式会社山武 Alarm analysis device, alarm analysis method, and alarm analysis program
JP4239932B2 (en) * 2004-08-27 2009-03-18 株式会社日立製作所 production management system

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110275489A (en) * 2018-03-13 2019-09-24 发那科株式会社 Data time series analysis device
CN110275489B (en) * 2018-03-13 2021-04-20 发那科株式会社 Time series data analysis device
US11188064B1 (en) 2021-05-04 2021-11-30 Ixden Ltd. Process flow abnormality detection system and method
WO2023233927A1 (en) * 2022-06-03 2023-12-07 オムロン株式会社 Abnormality detection device, abnormality detection method, and program

Also Published As

Publication number Publication date
JP2009064407A (en) 2009-03-26

Similar Documents

Publication Publication Date Title
JP5151556B2 (en) Process analysis apparatus, process analysis method, and process analysis program
JP5538597B2 (en) Anomaly detection method and anomaly detection system
JP5945350B2 (en) Equipment condition monitoring method and apparatus
JP5301310B2 (en) Anomaly detection method and anomaly detection system
JP5431235B2 (en) Equipment condition monitoring method and apparatus
JP5048625B2 (en) Anomaly detection method and system
JP5363927B2 (en) Abnormality detection / diagnosis method, abnormality detection / diagnosis system, and abnormality detection / diagnosis program
US9250625B2 (en) System of sequential kernel regression modeling for forecasting and prognostics
US20160110651A1 (en) Method of Sequential Kernel Regression Modeling for Forecasting and Prognostics
US20130024415A1 (en) Monitoring Method Using Kernel Regression Modeling With Pattern Sequences
KR101948604B1 (en) Method and device for equipment health monitoring based on sensor clustering
JP5498540B2 (en) Anomaly detection method and system
Mosallam et al. Component based data-driven prognostics for complex systems: Methodology and applications
WO2021241580A1 (en) Abnormality/irregularity cause identifying apparatus, abnormality/irregularity cause identifying method, and abnormality/irregularity cause identifying program
CN115698882A (en) Abnormal modulation cause identification device, abnormal modulation cause identification method, and abnormal modulation cause identification program
Filios et al. An agnostic data-driven approach to predict stoppages of industrial packing machine in near future
EP4160341A1 (en) Abnormal modulation cause identifying device, abnormal modulation cause identifying method, and abnormal modulation cause identifying program
JP2020071624A (en) Abnormality diagnosing apparatus, abnormality diagnosing method and program
Bond et al. A hybrid learning approach to prognostics and health management applied to military ground vehicles using time-series and maintenance event data
JP6798968B2 (en) Noise cause estimation device
WO2021241579A1 (en) Abnormal modulation cause identifying device, abnormal modulation cause identifying method, and abnormal modulation cause identifying program
KR102469229B1 (en) Fishbone diagram management system based on artifical intelligence and method thereof
Fabricio A review of Condition Based Maintenance approaches and its structure
JP7000266B2 (en) Analysis flow creation system and analysis flow creation method
Juselius Advanced condition monitoring methods in thermal power plants

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20101203

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20120209

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20120214

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20120411

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20121106

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20121119

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20151214

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 5151556

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees