講義予定 信頼性工学 第10回 2016.12.13 千葉大学 工学部 都市環境システム学科 山崎 文雄 http://ares.tu.chiba-u.jp/ 1 (1) 2016年10月 4日（火） 信頼性と信頼性工学（イントロダクション） (2) 2016年10月11日（火） 信頼性解析の基礎数理1（確率論の基礎） (3) 2016年10月18日（火） 信頼性解析の基礎数理2（信頼性の基本量） (4) 2016年10月25日（火） 信頼性解析の基礎数理3（故障率と確率分布） (5) 2016年11月 1日（火） データの統計解析1（統計データ処理） (6) 2016年11月 8日（火） データの統計解析2（最尤法と確率紙） (7) 2016年11月15日（火）データの統計解析3／システムの信頼性1 (8) 2016年11月22日（火） 中間試験 (9) 2016年12月 6日（火） システムの信頼性2（一般システムと信頼性設計） (10) 2016年12月13日（火） 故障モードの同定（FMEA, FTA, ETA） (11) 2016年12月20日（火） 構造物の信頼性工学1（破壊確率と信頼性指標） (12) 2017年 1月10日（火） 構造物の信頼性工学2（信頼性解析モデル） (13) 2017年 1月17日（火） モンテカルロ法 (14) 2017年 1月24日（火） 確率過程のモデル化 2 (15) 2017年 1月 31日（火） 期末試験 5.1 FMEA/FMECA 5．故障（破損）モードの同定 システム信頼性を評価するためには，カットセットに対 応して，故障（破損）モードをまず定義する必要がある． どのような故障モードが起こり，どの故障モードが致命 的であるか，系統的に同定する方法が必要となる． FMEA(故障モードと影響解析，failure mode and effect analysis) サブシステムの破損モードの生起によって，システムに どのような影響が現れるかを組織的に評価する方法． FMECA(故障モードと影響および致命度解析，failure mode, effect and criticality analysis) FMEA をさらに進めて，破損モードの起こり易さやその 影響の度合いをより詳細に評価する方法． （致命度の評価をFMEAに加えたもの） 5.1 FMEA/FMECA 5.2 FTA/ETA 3 4 FMEA解析の例 FMEA/FMECA の実施手順の例 大型舶用スターリングエンジン http://www.nmri.go.jp/eng/khirata/design/ch10/ch10_01.html 表2 影響度の評価点数 1. システムの機能（達成すべき使命）を明らかにする． 2. 故障の定義を明らかにする． 3. システムの構造，すなわち，システムを構成するサブ システム，コンポーネント，パーツを明確にする． 4. 各構成要素の故障モードとその原因を列挙する． 5. 故障モードの生起によるシステムへの影響を吟味す る． 6. 故障の生起確率や影響の度合いを検討する． 7. 故障防止，あるいは軽減対策を立てる． 表3 発生頻度の評価点数 エンジン全体を構成品に分割し，起こりうる故障とその原因を考える（表1）． その故障が構成品およびエンジン全体へ与える影響を考え，影響度の点数（表2）を記 入． 発生頻度（表3）を考え，故障等級（＝影響度×発生頻度）を計算． 5 表１ FMEA解析シートの例 大型舶用スターリングエンジン 故障等級（＝影響度×発生頻度） 7 このような方法で解析すると，故障等級が高い構成品（部品）ほど，信頼性が低いことを 表している。例えば，故障等級30以上を「要注意（再設計）」とすれば，影響度を減らす 6 か，あるいは発生頻度を減らすような設計が必要になることとなる。 FMECAの実施例（航空機の部品） 8 FMECAの各項目の記入方法(2) FMECAの各項目の記入方法(1) 6. 致命度：各故障の致命度を分類記号によって記入する． 1. 部品名，部品の機能：対象システム名，部品名，部品 の機能を記入する． 2. 故障または誤操作モード：考えられるハードウエアの 故障モードおよび誤操作のモードを別々に記入する． 3. 原因：各モードについて考えられる原因を列挙する． 4. 運用状態：システムがどのような運用状態の場合に当 該故障が起こり得るかを記入する． 5. システムへの影響：各故障によって生じる影響を詳細 に記入する． 9 5.2 FTA/ETA クラス１：安全(safe) – 当該故障がシステムの主要機能を低下させたり，人命に影 響を及ぼすには至らない． クラス２：限界(marginal) – 当該故障がシステム機能の低下をもたらすが，重大な 損害または人体への損傷に結びつくことなくコントロールすることが可能であ る． クラス３：致命(critical) – システムに重大な損害もしくは人命への危険につながる ので直ちに処置が必要である． クラス４：破局(catastrophic) – システム機能の停止もしくは人命の損失を伴う． 7. 発進の可否：当該故障が発生した場合，そのままの状 態でシステム（この例の場合は航空機）が発進可能か否 かをYes, No で記入する． 8. 識別およぶ処置：故障の表示や処置に関して表題部に 書かれた各事項について記入する． 9.故障の波及：故障が更に他の故障に波及するおそれが 10 ないかを記入する． 例1) 電気回路のFT FTA(故障の木解析，fault tree analysis) top event 故障（破損）のメカニズムを階層的に表現する方法． 対象とするシステムの破損を頂上事象として，それに関 する下位事象を“AND” あるいは“OR” の論理で結合して いく．(トップダウン解析手法 ) 下位事象について同様な操作を繰り返し，それ以上展 開することができない基本事象に到るまで続けていく． 基本事象の確率値を定めると，頂上事象の確率が求ま る． basic event 11 12 FTAとシステム構造関数 FTAに用いられる記号 要素の状態指標変数ai a 1 基本事象iが発生している i 0 基本事象iが発生していない 要素状態ベクトル a (a1 , a2 , , an )T システムの状態指標変数 a (a) (4.12) (4.13) S システム構造関数 1 頂上事象が発生している aS 0 頂上事象が発生していない FTAでは故障確率を計算するので，基本事象は故障となる． 指標変数aiを0または1の値をとる確率変数Aiで表すと Fi PAi 1 システムの不信頼度は FS PAS 1 要素の不信頼度は 13 FTAとシステム構造関数(2) FS E (A) (4.14’) 信頼度 ORゲート（直列システム） 0.9 0.95 0.8 電球 電池 電線 0.05 0.2 不信頼度 0.1 電球が点灯しない (4.16’) 0.323 電線 OR 0.99 0.01 電球 電池 － ＋ T AND 14 スイッチ スイッチ FS E (A ) T (4.20’) 直列ブロック図とFT図の関係 システム構造関数はシステムの内部構造を表す． ANDゲート（並列システム) (4.19’) 直列システムの信頼度： 電球 切れ OR 0.1 電流が 流れない Rs = 0.2476 構造関数を用いたORゲートの不信頼度： OR Fs1= A1 A2 ... An A1 A2 ... An 断線 0.2 15 電池 切れ 0.05 スイッチ 開放 0.01 Fs= 信頼度は Rs= 1- Fs= 16 セミナー開催不能となるFTA 並列ブロック図とFT図の関係 信頼度 0.95 0.9 0.05 0.95 電球 電池2 0.1 不信頼度 電線は正常 電池1 電球 電池１ － ＋ 0.05 電池２ － ＋ 電球が点灯しない 0.898 直並列システムの信頼度： OR Rs = 電球 切れ 電池切れ 0.1 0.0025 ANDゲートの不信頼度： Fs1= AND ORゲートの不信頼度： 電池1 切れ 0.05 電池2 切れ 0.05 Fs= 信頼度は Rs= 1- Fs= 17 １日/年間稼働日数＝１日/260日＝4×10-3 18 http://www.geocities.jp/takaro_u/fta.html 電気回路のETと生起確率 ETA(事象の木解析，event tree analysis) 損失の発端となる事象（初期事象）を出発点として，それ が拡大していく過程を様々な現象の発生有無によって枝 分かれしたイベントツリー(ET)を作成して解析するもの FTA: 頂上事象（結果）複数の基本事象（原因） ETA: １つの基本事象（原因） 頂上事象（結果） スイッチ スイッチ 電球の点灯 電球 電池 電池 － ＋ 電球 電線は正常 正常 0.9 低頻度高インパクトの事故について，シナリオ分析する場 合に用いる． ツリーに初期事象の発生頻度と事象の分岐確率を与える ことにより，中間あるいは末端に現われる種々の損失事 象がどの程度の頻度で起こりうるかを算出. 19 正常 正常 0.99 0.95 故障 正常 故障 電球をつけ ようとする 故障 正常 正常 故障 Ｎ段階の分岐に対し，2N通り の経路がある 故障 故障 正常 故障 P=0.9x0.95x0.99 Ｘ Ｘ Ｘ Ｘ Ｘ Ｘ Ｘ 20 ガス爆発事故のETの一例 0.1 不注意で再 度ガス供給 新型インフルエンザの感染シナリオ 0.8 ガス爆発 p1=0.1x0.5x0.9x0.1x0.8 着火源が存在 0.2 =0.0036 事故発生せず 0.9 0.9 ガス供給停 止成功 yes 0.5 0.1 異臭を検出 no 0.5 0.1 ガス漏れ ガス漏れの 発生確率p 換気不十分 0.9 事故発生せず 0.8 ガス爆発 p2=0.1x0.5x0.1x0.8 着火源が存在 0.2 0.8 事故発生せず ガス爆発 p3= 着火源が存在 0.2 =0.004 事故発生せず =0.04 事故発生せず http://www.tg-inet.co.jp/financial/products/opr.htm ガス爆発の発生確率 = =0.0476p 21 タイタニック号事故の概要 東田ほか， 2010地域安全学会論文集 22 タイタニック号は1912年4月10日 予定より約1ヶ月遅れて処女航海に向けて多数の著名人を乗せて 出航した。1ヶ月前の3月中旬に出航していれば、流氷の量は少なく氷山との衝突事故の発生確率は少 であったと言われている。次に、出航日においても予定より出航時刻が1時間遅れた。もし出来事の発 生が全て1時間早まっていたならば、衝突後救難無線を発した時刻が夜中の12時前になり、すみやか に周辺に存在していた船舶が救助に駆けつけることができたと言われている。 出航後、流氷原が行く手の海域に存在するという警告を他船から無線で受信していたが、船長はそ れほど深刻には受け止めていなかった。航海中も流氷原があるという警告を更に受けるが、速力を減 速することなく高速（20.5ノット）で航行を続けた。 やがて午後11:40に氷山に衝突するが、氷山発見は衝突の約450ｍ手前であった。この夜は月齢 26.1の新月に近い暗闇であった事と、珍しいほどの無風で鏡の様な海面で氷山の周囲に白波が全く 立っていなかったことも発見が遅れた要因であった。更に、見張り員が双眼鏡無しで見張り台に立って いた。ただし、この夜の条件では双眼鏡があったとしても、より早期の氷山の発見は無理であったとの 意見もある。 氷山の発見後、回避行動をとるが船腹をなでる様に氷山をかすめて通り、そのため却って船体の多数 区画の損傷によって多量の浸水をもたらし沈没に至ってしまった。 http://www.nmri.go.jp/sed/psa/titanic/ 独立行政法人 海上技術安全研究所 23 氷山との衝突の後、速やかには救難無線を出さず午前0時14分になって初めて救難無線を発してい る。しかし、タイタニック号から19海里の距離にいたカリフォルニア号の無線は午前零時を過ぎたその時 には切られていた。58海里の距離にいたカルパチア号が救難無線を受信し救助に駆けつけた。 午前零時44分には信号灯を打ち上げ、カリフォルニア号の乗組員がこの信号灯を視認したと言われ ているが、信号灯の意味するところを理解せず救助には向かわなかった。 午前2時20分についにタイタニック号は沈没してしまった。最初の救助船であるカルパチア号が到着し たのは午前4時10分頃であった。 24 ETA解析による遭難者数と累積発生確率 タイタニック号事故の イベントツリー 設定したイベント・ツリーの各分岐確率を用いて各シーケンスの発生頻度の点推定値が得 られる。発生したタイタニック号事故とまったく同一経過をたどるシーケンスの発生頻度は 1.07×10-8/航海と得られた。 しかし、異なった経過をたどっても結果的に多数の死者の出るシーケンスは多数存在する。 それらを集計した結果を図８に示す。遭難者数と累積発生確率の関係から タイタニック号 事故と同程度（死者1490人）あるいはそれ以上の事故が発生する確率は2.05×10-2/航海 となっている。 25 図８：解析結果 2016年中間試験 成績 受験者 最低点 最高点 平均点 51 20 100 64.14 中間 2016年中間試験 ヒストグラム 14 自主課題 中間試験があまり芳しくなかったと思われる人は，自主的に以下の課題のレポート （３－4枚程度）を提出することにより，配慮します． 以下のいずれかの課題を選択し，回答しなさい． 提出期限： 2017年１月17日（火） ４０点台 （１）HPに掲載したExcelシートは，「気象庁ホームページ」→「気象統計情報」 →「過去の気象 データ検索」から作成した東京都心における1876年から2015年までの月間降水量データであ る． http://www.data.jma.go.jp/obd/stats/etrn/index.php まず，6月，9月，および年間の降水量の年変化の様子をグラフで示し，平均値，標準偏差を 計算しなさい．これより，どのような傾向が読み取れるか考察しなさい．次に，Excelの「データ 分析」の「ヒストグラム」機能を用いて，年間降水量の140個のデータのヒストグラムを作成せ よ．またこの分布が正規分布で表現できるかどうかについて，カイ自乗検定を行いなさい． 120.0% 12 100.0% 頻度 累… 10 80.0% 8 60.0% 6 40.0% 4 100 参考： Excelによるヒストグラム作成とカイ自乗検定のやり方は，「環境プラニング演習II」講義ノート （http://ares.tu.chiba-u.jp/~lecturenote/index.htm）のLecture 10,11に掲載． 60点以下 次の級 90‐99 80‐89 70‐79 60‐69 50‐59 40‐49 30‐39 0.0% 20‐29 0 10‐19 20.0% 0 2 1‐9 頻度 ４０点未満 26 点数区間 27 （２）最近，注目を集めた事故（トラブル）について，ETAまたはFTAを用いて，原因と なる要因と結果の関係を構造的に示しなさい．さらに，分岐または基本事象の確率を 適宜与えて，事故（トラブル）の発生確率を試算しなさい． 対象とするテーマの例： •福島第一原発事故，巨大地震の発生見逃がし，タイ洪水による製造業への影響 •JR西福知山線脱線事故，2012中央道トンネル崩落，高齢者による交通事故多発 28 •その他，火災，医療事故，航空事故，交通事故，食の安全など
© Copyright 2024 Paperzz
