ホーム > 特許ランキング > 国立大学法人横浜国立大学
知財求人 - 知財ポータルサイト「IP Force」
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2025098528
(43)【公開日】2025-07-02
(54)【発明の名称】情報処理システム、情報処理方法、プログラム
(51)【国際特許分類】
G01M 7/02 20060101AFI20250625BHJP
G01M 99/00 20110101ALI20250625BHJP
【ＦＩ】
G01M7/02 H
G01M99/00 Z
【審査請求】未請求
【請求項の数】9
【出願形態】ＯＬ
(21)【出願番号】P 2023214722
(22)【出願日】2023-12-20
【新規性喪失の例外の表示】特許法第３０条第２項適用申請有り 〔１〕開催日（公開日） 令和５年９月１８日 集会名、開催場所 防災ＫＯＫＵＤＡＩ（国立大学法人横浜国立大学 主催）（会期：令和５年９月１７日～１８日） 先端科学高等研究院 ＹＮＵ Ｄｉａｌｏｇｕｅ＆サイエンスカフェ 「デジタル社会における防災」 国立大学法人横浜国立大学 経営講義棟１号棟１Ｆ １０８室 ＜資 料＞ 防災ＫＯＫＵＤＡＩ 開催概要 ＜資 料＞ 研究発表 スライド資料
【国等の委託研究の成果に係る記載事項】（出願人による申告）令和５年度、消防庁「ＩＣＴ・ＩｏＴ技術を活用した石油コンビナート災害対応システムの開発と社会実装」委託研究、産業技術力強化法第１７条の適用を受ける特許出願
(71)【出願人】
【識別番号】504182255
【氏名又は名称】国立大学法人横浜国立大学
(74)【代理人】
【識別番号】110001634
【氏名又は名称】弁理士法人志賀国際特許事務所
(72)【発明者】
【氏名】吉田 聖一
【テーマコード（参考）】
2G024
【Ｆターム（参考）】
2G024AD38
2G024BA12
2G024CA04
2G024DA12
2G024FA06
2G024FA14
2G024FA15
(57)【要約】
【課題】異なる複数のタンクの損傷に関する指標をより簡易な準備で容易に算出することができる情報処理システムを提供する。
【解決手段】側板と底板とを溶接して形成されたタンクの振動による前記溶接部近傍の隅角部の浮き上がり変位波形のピークが示す値と、前記タンクに応じた係数とを隅角部ひずみ算出式に入力し、前記溶接部近傍における隅角部ひずみを算出する。
【選択図】図１
【解決手段】側板と底板とを溶接して形成されたタンクの振動による前記溶接部近傍の隅角部の浮き上がり変位波形のピークが示す値と、前記タンクに応じた係数とを隅角部ひずみ算出式に入力し、前記溶接部近傍における隅角部ひずみを算出する。
【選択図】図１
【特許請求の範囲】
【請求項1】
振動による、タンクを構成する側板と底板との接合部近傍を示す隅角部の浮き上がり変位波形のピークが示す値と、前記タンクに応じた係数とを隅角部ひずみ算出式に入力し、前記隅角部における隅角部ひずみを算出する隅角部ひずみ算出部を備える情報処理システム。
【請求項2】
前記隅角部ひずみ算出部は、隅角部ひずみをε、前記底板を構成する鋼板の種別に応じた係数をＣ、前記タンクにおける底部の液体圧力をｐ、前記浮き上がり変位のピーク値をδ、前記底板の板厚をｔａ、前記底板の降伏応力をσｙとして、式（１）により表される前記隅角部ひずみ算出式を用いて、前記隅角部ひずみを算出する
請求項１に記載の情報処理システム。
【数1】
【請求項3】
前記隅角部ひずみ算出部は、前記振動を示す波形に基づいて前記タンクに生じた１回目の前記浮き上がり変位のピーク値に対応する第一の前記隅角部ひずみと、前記波形に基づいて前記タンクに生じた２回目以降の前記浮き上がり変位のピーク値それぞれに対応する第二の前記隅角部ひずみと、を算出し、
前記第一の前記隅角部ひずみの最大と最小との差を示すひずみ振幅を式（２）で示す第一ひずみ振幅算出式を用いて算出し、
前記第二の前記隅角部ひずみの最大と最小との差を示すひずみ振幅を式（３）で示す第二ひずみ振幅算出式を用いて算出する、
ひずみ振幅算出部とを備える請求項２に記載の情報処理システム。
前記第一の前記隅角部ひずみの最大と最小との差を示すひずみ振幅を式（２）で示す第一ひずみ振幅算出式を用いて算出し、
【数2】
【数3】
【請求項4】
前記波形に基づいて前記タンクに生じた前記浮き上がり変位のピークの繰り返し発生に応じた当該繰り返しの回数における前記隅角部の累積損傷度を、累積損傷度算出式を用いて算出する累積損傷度算出部と、
を備える請求項３に記載の情報処理システム。
を備える請求項３に記載の情報処理システム。
【請求項5】
前記隅角部ひずみ算出部は、対象のタンクに応じた前記係数であって、前記底板を構成する鋼板の種別に応じた係数をＣ、前記タンクにおける底部の液体圧力をｐ、前記浮き上がり変位のピーク値をδ、前記底板の板厚をｔａ、前記底板の降伏応力をσｙと、前記隅角部ひずみ算出式を用いて、前記隅角部ひずみを算出する
請求項２に記載の情報処理システム。
請求項２に記載の情報処理システム。
【請求項6】
発生した前記振動の波形と、前記タンクに生じた前記浮き上がり変位を１質点系非線形水平ばねモデルと、当該モデルから得られた水平方向の振動ピーク値と浮き上がり変位算出式とを用いて、前記浮き上がり変位のピーク値を算出する浮き上がり変位算出部と、
を備える請求項１から請求項５の何れか一項に記載の情報処理システム。
を備える請求項１から請求項５の何れか一項に記載の情報処理システム。
【請求項7】
前記振動の波形として地震動の加速度波形を前記地震動が発生している際に取得する取得部と、
を備える請求項６に記載の情報処理システム。
を備える請求項６に記載の情報処理システム。
【請求項8】
振動による、タンクを構成する側板と底板との接合部近傍を示す隅角部の浮き上がり変位波形のピークが示す値と、前記タンクに応じた係数とを隅角部ひずみ算出式に入力し、前記隅角部における隅角部ひずみを算出する
情報処理方法。
情報処理方法。
【請求項9】
情報処理システムのコンピュータを、
振動による、タンクを構成する側板と底板との接合部近傍を示す隅角部の浮き上がり変位波形のピークが示す値と、前記タンクに応じた係数とを隅角部ひずみ算出式に入力し、前記隅角部における隅角部ひずみを算出する隅角部ひずみ算出手段、
として機能させるプログラム。
振動による、タンクを構成する側板と底板との接合部近傍を示す隅角部の浮き上がり変位波形のピークが示す値と、前記タンクに応じた係数とを隅角部ひずみ算出式に入力し、前記隅角部における隅角部ひずみを算出する隅角部ひずみ算出手段、
として機能させるプログラム。
【発明の詳細な説明】
【技術分野】
【0001】
この開示は、情報処理システム、情報処理方法、プログラムに関する。
【背景技術】
【0002】
石油などの対象を保管するタンクは地震動などの振動の影響により変形する場合がある。例えば短周期地震動（高周波振動）によるタンク側板主体の振動であるバルジングや、長周期地震動（低周波数振動）による液面の揺動であるスロッシングが発生する。また短周期地震動（高周波振動）により水平方向の慣性力でタンクが横方向に振動するロッキング振動が起き、タンクを構成する側板と底板との接合部近傍を示す隅角部が地面から浮き上がり、これが地震発生などの振動発生時におけるタンクの主要な損傷原因の一つであることが知られている。
【0003】
このようなタンクの損傷を推定するものとして、石油タンクの健全性評価システム（ＳＵＳＴＡＩＮＥＲ）が開発されている。また関連する技術が特許文献１に開示されている。
【先行技術文献】
【特許文献】
【0004】
【特許文献1】特開２００６－２９９２５号公報
【発明の概要】
【発明が解決しようとする課題】
【0005】
ところで上述のようなタンクは、その大きさや内部に保管する対象がそれぞれ異なる場合がある。これら異なる複数のタンクの損傷に関する指標をより簡易な準備で容易に算出することのできる技術が求められている。
【0006】
この開示は、上記の課題を解決する情報処理システム、情報処理方法、プログラムを提供することを目的としている。
【課題を解決するための手段】
【0007】
この開示の第１の態様によれば、情報処理システムは、振動による、タンクを構成する側板と底板との接合部近傍を示す隅角部の浮き上がり変位波形のピークが示す値と、前記タンクに応じた係数とを隅角部ひずみ算出式に入力し、前記隅角部における隅角部ひずみを算出する隅角部ひずみ算出部を備える。
【0008】
この開示の第２の態様によれば、情報処理方法は、振動による、タンクを構成する側板と底板との接合部近傍を示す隅角部の浮き上がり変位波形のピークが示す値と、前記タンクに応じた係数とを隅角部ひずみ算出式に入力し、前記隅角部における隅角部ひずみを算出する隅角部ひずみ算出部を備える。
【0009】
この開示の第３の態様によれば、プログラムは、情報処理システムのコンピュータを、振動による、タンクを構成する側板と底板との接合部近傍を示す隅角部の浮き上がり変位波形のピークが示す値と、前記タンクに応じた係数とを隅角部ひずみ算出式に入力し、前記隅角部における隅角部ひずみを算出する隅角部ひずみ算出手段、として機能させる。
【発明の効果】
【0010】
異なる複数のタンクの損傷に関する指標をより簡易な準備で容易に算出することができる情報処理システムを提供することができる。
【図面の簡単な説明】
【0011】
【図1】この開示の一実施形態による情報処理システムの概略構成図である。
【図2】この開示の一実施形態による情報処理装置のハードウェア構成図である。
【図3】この開示の一実施形態による情報処理装置の機能ブロック図である。
【図4】この開示の一実施形態による浮き上がり変位の算出処理の概要を示す図である。
【図5】この開示の一実施形態による復元力特性ＱＲを示すグラフである。
【図6】この開示の一実施形態による情報処理装置の処理フローを示す図である。
【図7】この開示の一実施形態による地震動の地震動変位に応じた隅角部の浮き上がり変位δの波形を示す。
【図8】この開示の一実施形態による隅角部ひずみと浮き上がり変位との関係の比較例を示す第一の図である。
【図9】この開示の一実施形態による隅角部ひずみと浮き上がり変位との関係の比較例を示す第二の図である。
【図10】この開示の一実施形態による情報処理装置の他の構成を示す図である。
【図11】この開示の一実施形態による他の構成による情報処理装置の処理フローを示す図である。
【発明を実施するための形態】
【0012】
図１はこの開示の一実施形態による情報処理システムの概略構成図である。
図１に示すようにこの開示の情報処理システム１００は情報処理装置１０を備える。情報処理システム１００は、情報処理装置１０以外の情報処理装置も備えており、他の情報処理装置では、地震計やその他のセンサから得られた情報を入力したデータを用いて、地震動観測記録処理、タンク毎の地震動算出処理、タンクの側板や底板の損傷評価、スロッシング影響評価、波高、タンク内包対象の溢流、浮き屋根損傷度算出などの処理を行っている。これらの処理結果は、タンク算出結果ＤＢ３０が記憶する。タンク毎地震動ＤＢ２０には、他の情報処理装置が算出したタンク毎の地震動算出結果が記録されている。
図１に示すようにこの開示の情報処理システム１００は情報処理装置１０を備える。情報処理システム１００は、情報処理装置１０以外の情報処理装置も備えており、他の情報処理装置では、地震計やその他のセンサから得られた情報を入力したデータを用いて、地震動観測記録処理、タンク毎の地震動算出処理、タンクの側板や底板の損傷評価、スロッシング影響評価、波高、タンク内包対象の溢流、浮き屋根損傷度算出などの処理を行っている。これらの処理結果は、タンク算出結果ＤＢ３０が記憶する。タンク毎地震動ＤＢ２０には、他の情報処理装置が算出したタンク毎の地震動算出結果が記録されている。
【0013】
情報処理装置１０は、タンク毎地震動ＤＢ２０より、対象のタンクの地震動算出結果である地震動波形を取得する。また情報処理装置１０は、対象のタンクに関する係数を所定の記憶部から取得する。情報処理装置１０は、振動による、タンクを構成する側板と底板との接合部近傍を示す隅角部の浮き上がり変位波形のピークが示す値と、タンクに応じた係数とを隅角部ひずみ算出式に入力し、隅角部における隅角部ひずみを算出する。
【0014】
隅角部ひずみを算出する場合、有限要素法の手法を用いてタンク毎に隅角ひずみを算出する必要があった。しかしながら本開示の情報処理装置１０の手法により、どのような大きさや保管対象のタンクの隅角部ひずみも、有限要素法を行う場合の事前準備や、膨大な有限要素法の処理を行うことなく、簡易に算出することができる。
【0015】
図２は情報処理装置のハードウェア構成図である。
図２に示すように情報処理装置１０は、ＣＰＵ（Central Processing Unit）１０１、ＲＯＭ（Read Only Memory）１０２、ＲＡＭ（Random Access Memory）１０３、ＨＤＤ（Hard Disk Drive）１０４、通信モジュール１０５等の各ハードウェアを備えたコンピュータである。情報処理システム１００は、複数の情報処理装置１０を用いて構成されたものであってよい。つまり他の情報処理装置も、情報処理装置１０と同様のコンピュータである。
図２に示すように情報処理装置１０は、ＣＰＵ（Central Processing Unit）１０１、ＲＯＭ（Read Only Memory）１０２、ＲＡＭ（Random Access Memory）１０３、ＨＤＤ（Hard Disk Drive）１０４、通信モジュール１０５等の各ハードウェアを備えたコンピュータである。情報処理システム１００は、複数の情報処理装置１０を用いて構成されたものであってよい。つまり他の情報処理装置も、情報処理装置１０と同様のコンピュータである。
【0016】
図３は情報処理装置の機能ブロック図である。
情報処理装置１０のＣＰＵ１０１がタンク管理プログラムを実行する。これにより、情報処理装置１０は、取得部１１、浮き上がり変位算出部１２、隅角部ひずみ算出部１３、ひずみ振幅算出部１４、累積損傷度算出部１５、出力部１６の各機能を発揮する。
情報処理装置１０のＣＰＵ１０１がタンク管理プログラムを実行する。これにより、情報処理装置１０は、取得部１１、浮き上がり変位算出部１２、隅角部ひずみ算出部１３、ひずみ振幅算出部１４、累積損傷度算出部１５、出力部１６の各機能を発揮する。
【0017】
取得部１１は、情報処理装置１０が算出処理に利用する各種情報を外部のデータベース（ＤＢ）や他の情報処理装置から取得する。
浮き上がり変位算出部１２は、発生した地震動の加速度波形に基づいてタンクに生じた浮き上がり変位を、１質点系非線形水平ばねモデルによる地震応答解析で当該モデルから得られた水平方向の振動値と浮き上がり変位算出式とを用いて算出する。
浮き上がり変位算出部１２は、発生した地震動の加速度波形に基づいてタンクに生じた浮き上がり変位を、１質点系非線形水平ばねモデルによる地震応答解析で当該モデルから得られた水平方向の振動値と浮き上がり変位算出式とを用いて算出する。
【0018】
隅角部ひずみ算出部１３は、振動によるタンクの隅角部の浮き上がり変位波形のピークが示す値と、タンクに応じた係数とを隅角部ひずみ算出式に入力し、隅角部ひずみを算出する。
【0019】
ひずみ振幅算出部１４は、振動に応じて底板と側板との隅角部（底板と側板との接合部近傍）の浮き上がり変位から求まる隅角部ひずみの最大と最小との差を示すひずみ振幅を算出する。隅角部ひずみの最大と最小との差は浮き上がり変位が零からピークまでのひずみの変動の値である。
【0020】
累積損傷度算出部１５は、タンクに生じた浮き上がり変位のピークの繰り返し発生に応じた隅角部ひずみ振幅の当該繰り返しの回数における隅角部の累積損傷度を算出する。
【0021】
出力部１６は、情報処理装置１０が算出した浮き上がり変位量や、隅角部ひずみ、ひずみ振幅、累積損傷度などの情報を記憶部に記憶する、または表示部に表示する処理を行う。
【0022】
図４は浮き上がり変位の算出処理の概要を示す図である。
図４（４ａ）は、地震動によるタンクの隅角部の浮き上がりの状態を示している。隅角部の浮き上がりは、短周期地震動（高周波振動）により水平方向の慣性力でタンクやタンク内部の保管対象（石油など）が横方向に振動するロッキング振動が主要な原因となる。したがって本開示において情報処理装置１０は、このロッキング振動を、１質点系非線形水平ばねモデルによりモデル化し、保管対象を含むタンク全体の質点の応答変位Δを求める。応答変位Δは保管対象を含むタンク全体の質点の横方向の移動量である。
図４（４ａ）は、地震動によるタンクの隅角部の浮き上がりの状態を示している。隅角部の浮き上がりは、短周期地震動（高周波振動）により水平方向の慣性力でタンクやタンク内部の保管対象（石油など）が横方向に振動するロッキング振動が主要な原因となる。したがって本開示において情報処理装置１０は、このロッキング振動を、１質点系非線形水平ばねモデルによりモデル化し、保管対象を含むタンク全体の質点の応答変位Δを求める。応答変位Δは保管対象を含むタンク全体の質点の横方向の移動量である。
【0023】
【数1】
【0024】
式（１）が示す１質点系非線形水平ばねモデルにおいて、Ｍ１はタンクとその保管対象全体における有効質量（応答加速度に影響を与えるマスの質量）、Ｃｅは質点の１回の応答速度（１回の質点の横方向の移動）に対応する減衰、ＱＲはタンクの横方向に加わった力と隅角部の浮き上がり変位に応じた質点の復元力特性、Δは質点の応答変位、ｕｇは地震動変位（地震動波形）を示す。また式（１）においてΔやｕの上に付された［・・］は二階時間微分を示し、またΔの上に付された［・］は一階時間微分を示す。タンクとその保管対象全体における有効質量Ｍ１、質点の１回の応答速度に対応する減衰Ｃｅ、タンクの横方向に加わった力と隅角部の浮き上がり変位に応じた質点の復元力特性ＱＲは、タンクに特有の係数であり、タンク毎に異なってよい。
【0025】
【0026】
情報処理装置１０は式（１）が示す１質点系非線形水平ばねモデルに基づいて地震動の地震波形（地震動変位ｕｇの二階時間微分である地震動加速度（ｕｇのｕの上に［・・］付与））に応じた、質点の応答変位Δの波形を算出する。その質点の応答変位Δの波形に対応する、タンク隅角部の浮き上がり変位δを浮き上がり変位算出式（２）を用いて算出する。なお、以下の文章において数式中の地震動加速度（ｕｇのｕの上に［・・］付与）は「ｕｇ’」と別文字で表して説明する。
【0027】
【数2】
【0028】
浮き上がり変位算出式（２）において、Ｈ１はタンク全体の質点の高さ、Ｄは円筒形のタンクの直径、Ｋｂは１質点系非線形水平ばねモデルにおける初期浮き上がり段階のばね定数である。タンク全体の質点の高さＨ１、タンクの直径Ｄ、ばね定数Ｋｂは、タンクに特有の係数であり、タンク毎に異なってよい。そして、情報処理装置１０は、隅角部の浮き上がり変位波形のピークが示す値δと、タンクに応じた係数とを隅角部ひずみ算出式（３）に入力し、タンクの隅角部における隅角部ひずみεを算出する。
【0029】
【数3】
【0030】
浮き上がり変位算出式において、εは隅角部ひずみ、Ｃは底板（アニュラ板）を構成する鋼板の種別に応じた係数、ｐはタンクにおける底部の液体圧力（Ｎ（ニュートン）／ｍｍ２）、δは浮き上がり変位のピーク値（ｍｍ（ミリメートル））、ｔａは底板（アニュラ板）の板厚（ｍｍ）、σｙは底板（アニュラ板）の降伏応力（Ｎ／ｍｍ２）である。これら底板を構成する鋼板の種別に応じた係数Ｃ、タンクにおける底部の液体圧力ｐ、浮き上がり変位のピーク値δ、底板（アニュラ板）の板厚ｔａ、底板（アニュラ板）の降伏応力σｙはタンクに特有の係数であり、タンク毎に異なってよい。
【0031】
図６は情報処理装置の処理フローを示す図である。
図７は地震動の地震動加速度に応じた隅角部の浮き上がり変位δの波形を示す。
次に情報処理装置の処理フローについて順を追って説明する。
まず情報処理装置１０の取得部１１は情報取得タイミングを検知する（ステップＳ１０１）。情報取得タイミングは例えば地震発生時などであってよい。または情報取得タイミングは地震発生時以外のタイミングであってもよい。取得部１１は算出対象のタンクの識別情報をＨＤＤ１０４等の記憶部から取得する。算出対象のタンクは、たとえば地震発生エリアの全てのタンクであってもよい。取得部１１は、算出対象を決定するために地震発生の震源地の情報などを外部の情報提供サーバから取得してもよい。算出対象のタンクの決定はどのような手法であってもよい。算出対象のタンクは予め決められていてもよい。取得部１１は算出対象の１つまたは複数のタンクの識別情報（ＩＤなど）に紐づいて記録されるタンク毎の地震波形を示す地震動加速度ｕｇ’（ｕｇのｕの上に［・・］を付した数式中文字に相当）をタンク毎地震動ＤＢ２０から取得する（ステップＳ１０２）。取得部１１が取得する地震波形（地震動加速度ｕｇ’）は、地震計が記録した地震動あるいは模擬地震動から算出対象として決定したタンクの位置における地盤条件に応じて算出した地震波形でもよい。取得部１１は、タンクの識別情報とそのタンクの位置における地震波形（地震動加速度ｕｇ’）との関係を示す地震波形データをタンク毎地震ＤＢ２０から取得する。
図７は地震動の地震動加速度に応じた隅角部の浮き上がり変位δの波形を示す。
次に情報処理装置の処理フローについて順を追って説明する。
まず情報処理装置１０の取得部１１は情報取得タイミングを検知する（ステップＳ１０１）。情報取得タイミングは例えば地震発生時などであってよい。または情報取得タイミングは地震発生時以外のタイミングであってもよい。取得部１１は算出対象のタンクの識別情報をＨＤＤ１０４等の記憶部から取得する。算出対象のタンクは、たとえば地震発生エリアの全てのタンクであってもよい。取得部１１は、算出対象を決定するために地震発生の震源地の情報などを外部の情報提供サーバから取得してもよい。算出対象のタンクの決定はどのような手法であってもよい。算出対象のタンクは予め決められていてもよい。取得部１１は算出対象の１つまたは複数のタンクの識別情報（ＩＤなど）に紐づいて記録されるタンク毎の地震波形を示す地震動加速度ｕｇ’（ｕｇのｕの上に［・・］を付した数式中文字に相当）をタンク毎地震動ＤＢ２０から取得する（ステップＳ１０２）。取得部１１が取得する地震波形（地震動加速度ｕｇ’）は、地震計が記録した地震動あるいは模擬地震動から算出対象として決定したタンクの位置における地盤条件に応じて算出した地震波形でもよい。取得部１１は、タンクの識別情報とそのタンクの位置における地震波形（地震動加速度ｕｇ’）との関係を示す地震波形データをタンク毎地震ＤＢ２０から取得する。
【0032】
また取得部１１は、算出対象として決定したタンクごとの隅角部の浮き上がり変位の算出に用いる係数を所定の記憶部から取得する（ステップＳ１０３）。所定の記憶部は外部のデータベースでもよいし、自装置のＨＤＤ１０４であってもよい。隅角部の浮き上がり変位の算出に用いる係数は、上述したように、タンクとその保管対象全体における有効質量Ｍ１、質点の応答速度に対応する減衰Ｃｅ、タンクの横方向に加わった力と隅角部の浮き上がり変位に応じた質点の復元力特性ＱＲ、タンク全体の質点の高さＨ１、タンクの直径Ｄ、ばね定数Ｋｂ、底板を構成する鋼板の種別に応じた係数Ｃ、タンクにおける底部の液体圧力ｐ、底板（アニュラ板）の板厚ｔａ、底板（アニュラ板）の降伏応力σｙである。
【0033】
そして浮き上がり変位算出部１２は、１つまたは複数のタンクのうち、算出対象のタンクとその保管対象全体における有効質量Ｍ１、質点の応答速度に対応する減衰Ｃｅ、タンクの横方向に加わった力と隅角部の浮き上がり変位に応じた質点の復元力特性ＱＲと、ステップＳ１０２で取得した当該タンクに関する地震動加速度ｕｇ’を１質点系非線形水平ばねモデル（１）入力する。これにより浮き上がり変位算出部１２は、地震動加速度ｕｇ’に応じた質点の応答変位Δを算出する。浮き上がり変位算出部１２は、地震動加速度ｕｇ’が示す波形の所定間隔の地震動加速度ｕｇ’の値毎に、応答変位Δの算出を繰り返す。これにより浮き上がり変位算出部１２は、地震動加速度ｕｇ’に応じた質点の応答変位Δの波形を算出する（ステップＳ１０４）。１質点系非線形水平ばねモデル（１）を用いた地震動加速度ｕｇ’に応じた質点の応答変位Δの波形は他の情報処理装置が算出し、これを情報処理装置１０が取得してもよい。
【0034】
また浮き上がり変位算出部１２は、地震動加速度ｕｇ’に応じた質点の応答変位Δの値と、タンク全体の質点の高さＨ１、タンクの直径Ｄ、ばね定数Ｋｂを浮き上がり変位算出式（２）に入力して浮き上がり変位δを算出する。浮き上がり変位算出部１２は、地震動加速度ｕｇ’に応じた質点の応答変位Δが示す波形の所定間隔の応答変位Δの値毎に、浮き上がり変位δの算出を繰り返し、浮き上がり変位δの時刻歴応答波形（図７）を算出する（ステップＳ１０５）。
【0035】
図７が示す地震動の地震動加速度ｕｇ’に応じた隅角部の浮き上がり変位δの波形から、浮き上がり変位δのピークがδ１，δ２，・・・，δ１０と１０回発生していることが分かる。浮き上がり変位算出部１２は、算出した浮き上がり変位δの波形のデータから、ピーク値を特定し、それら浮き上がり変位δのピーク値（図７の場合は、δ１，δ２，・・・，δ１０の１０回発生の各ピーク値）を隅角部ひずみ算出部１３へ出力する。
【0036】
なお、浮き上がり変位δの時刻歴応答波形の算出や、浮き上がり変位δのピーク値の特定も、他の情報処理装置が算出し、これらの結果を情報処理装置１０が取得してもよい。つまりこの場合、情報処理装置１０は、以降の説明の隅角部ひずみの算出からの処理を行うようにしてもよい。
【0037】
隅角部ひずみ算出部１３は地震動が発生している間の隅角部の浮き上がり変位δのピーク値（図７の場合は、δ１，δ２，・・・，δ１０の１０回発生の各ピーク値）を取得する。隅角部ひずみ算出部１３は１回目の発生の浮き上がり変位のピーク値δ１と、算出対象のタンク特有の係数である、底板を構成する鋼板の種別に応じた係数Ｃ、タンクにおける底部の液体圧力ｐ、浮き上がり変位のピーク値δ、底板（アニュラ板）の板厚ｔａ、底板（アニュラ板）の降伏応力σｙを、隅角部ひずみ算出式（３）に入力する。これにより隅角部ひずみ算出部１３は、１回目の発生の浮き上がり変位のピーク値δ１に対応する、タンクの隅角部における隅角部ひずみε１を算出する。
【0038】
同様に、隅角部ひずみ算出部１３は、２回目発生の浮き上がり変位のピーク値δ２と、算出対象のタンク特有の係数である、底板を構成する鋼板の種別に応じた係数Ｃ、タンクにおける底部の液体圧力ｐ、浮き上がり変位のピーク値δ、底板（アニュラ板）の板厚ｔａ、底板（アニュラ板）の降伏応力σｙを、隅角部ひずみ算出式（３）に入力する。これにより隅角部ひずみ算出部１３は、２回目の発生の浮き上がり変位のピーク値δ２に対応する、タンクの隅角部における隅角部ひずみε２を算出する。これにより、隅角部ひずみ算出部１３は、各ピークが発生した隅角部の浮き上がり変位δ（δ１，δ２，・・・，δ１０）毎の隅角部ひずみε（ε１，ε２，・・・，ε１０）を算出する（ステップＳ１０６）。隅角部ひずみ算出部１３は、各ピークが発生した隅角部の浮き上がり変位δ毎の隅角部ひずみεを、ひずみ振幅算出部１４へ出力する。
【0039】
なお情報処理装置１０は少なくとも、隅角部ひずみεを算出するものであってもよい。関連する技術においては隅角部ひずみεを算出する場合には、有限要素法を用いて算出する必要があった。有限要素法を用いて隅角部ひずみεを算出する場合、事前準備に多くの時間とコスト、労力がかかり、地震発生直後などの所望のタイミングで、異なる大きさや、異なる対象を保管するタンクそれぞれの隅角部ひずみεを直ちに算出することができなかった。本開示によれば、有限要素法を用いずに隅角部ひずみ算出式（３）の簡易な式で、事前準備のための多くの時間やコスト、労力をかけずに直ちに異なる大きさや、異なる対象を保管するタンクそれぞれの隅角部ひずみεを直ちに算出することができる。これにより、異なる大きさや、異なる対象を保管するタンクそれぞれに関する振動発生による隅角部ひずみεに基づいて、損傷の程度を迅速に把握できる。
【0040】
ひずみ振幅算出部１４は各ピークが発生した隅角部の浮き上がり変位δ毎の隅角部ひずみεを取得する。ひずみ振幅算出部１４は、１回目発生の浮き上がり変位のピーク値δ１に対応するタンクの隅角部における隅角部ひずみε１を、第一のひずみ振幅算出式（４）に入力し、１回目発生の浮き上がり変位のピーク値δ１に対応するひずみ振幅Δεを算出する（ステップＳ１０７）。
【0041】
【数4】
【0042】
ひずみ振幅算出部１４は、２回目以降発生の浮き上がり変位のピーク値δ２～δ１０に対応するタンクの隅角部における隅角部ひずみε２～ε１０を、第二のひずみ振幅算出式（５）に入力し、２回目以降発生の浮き上がり変位のピーク値δそれぞれに対応するひずみ振幅Δεを算出する（ステップＳ１０８）。
【0043】
【数5】
【0044】
なお第二のひずみ振幅算出式（５）における係数０．３の値はタンクの大きさなどの種別によって異なってよい。
【0045】
以上の処理により、隅角部ひずみ算出部１３は、振動を示す波形に基づいてタンクに生じた１回目の浮き上がり変位のピーク値に対応する第一の前記隅角部ひずみと、波形に基づいてタンクに生じた２回目以降の前記浮き上がり変位のピーク値それぞれに対応する第二の前記隅角部ひずみと、を算出している。また 以上の処理によりひずみ振幅算出部１４は、波形が示す振動に応じて底板の隅角部の浮き上がり変位から求まる隅角部ひずみの最大と最小との差を示すひずみ振幅であって、第一の隅角部ひずみが発生した際のひずみ振幅を式（４）で示す第一ひずみ振幅算出式を用いて算出し、第二の隅角部ひずみが発生した際のひずみ振幅を式（５）で示す第二ひずみ振幅算出式を用いて算出している。
【0046】
ひずみ振幅算出部１４は、１回目発生の浮き上がり変位のピーク値δ１に対応するひずみ振幅Δεや、２回目以降発生の浮き上がり変位のピーク値δそれぞれに対応するひずみ振幅Δεを、累積損傷度算出部１５へ出力する。
【0047】
累積損傷度算出部１５は、１回目発生の浮き上がり変位のピーク値δ１に対応するひずみ振幅Δεや、２回目以降発生の浮き上がり変位のピーク値δそれぞれに対応するひずみ振幅Δεを取得する。今、浮き上がり変位δに対応する隅角部ひずみεが式（６）で表され、
【0048】
【数6】
【0049】
浮き上がり変位のピーク値が発生する回数（サイクル）ごとの係数をｃとすると、ひずみ振幅Δεは、式（７）で表すことができる。なおひずみ振幅とは、浮き上がり変位が零からピークまでのひずみの変動の値である。
【0050】
【数7】
【0051】
ここで地震発生時の浮き上がりに対する隅角部の強度評価項目は、低サイクル疲労であるとする。隅角部が繰り返し浮き上がる現象は低サイクル疲労である。低サイクル疲労による損傷度は、飯田の最適疲労曲線に当てはめることで得られる。累積損傷度算出部１５は、飯田の最適疲労曲線をＳ－Ｎ曲線として、疲労寿命Ｎｃを式（８）で求める。
【0052】
【数8】
【0053】
なお、式（８）においてｆεは、ひずみ振幅に対する安全率でその値を１．１とする。なお安全率とは、基準疲労寿命に対する許容疲労寿命の比率であり、基準疲労寿命÷許容疲労寿命である。
【0054】
累積損傷度算出部１５は、マイナー則を適用し、ｊ回目の地震に対する損傷度Ｄｊを、当該地震の地震波形（地震動加速度ｕｇ’）において発生する浮き上がり変位のピークの回数をｉとして、下記の累積損傷度算出式（９）により算出する（ステップＳ１０９）。
【0055】
【数9】
【0056】
ここで、Δεiは、浮き上がり変位のピークの回数に対応する浮き上がりのひずみ振幅を示す。またＮｃ（Δεi）は、浮き上がり変位のピークの回数に対応する浮き上がりのひずみ振幅Δεiの疲労寿命を示す。ｆｃは繰り返し回数に対する安全率を示す。
【0057】
ここで、ｊ回目の地震に対する損傷度Ｄｊの複数の地震の合計である累積損傷度Ｄが１を超えた時が寿命になる。したがって、累積損傷度算出部１５は、ｊ回の地震に対する累積損傷度を、各１回の地震に対する損傷度Ｄｊを合計して算出する（ステップＳ１１０）。また累積損傷度算出部１５は、式（１０）が示す判定式を満たさない場合、つまり、累積損傷度Ｄが１を超えた場合、タンクの寿命であると判定し、警告情報などを所定の出力先に出力するようにしてよい。
【0058】
【数10】
【0059】
情報処理装置１は全ての算出対象のタンクについて処理を終了したかを判定する（ステップＳ１１１）。情報処理装置１は全ての算出対象のタンクについて処理を終了していない場合には、ステップＳ１０３からの処理を繰り返す。
【0060】
以上の処理によれば、情報処理装置１０は、事前準備のための多くの時間やコスト、労力をかけずに算出した所定の１つまたは複数のタンクにおいて、隅角部の浮き上がり変位δのピーク毎の隅角部ひずみεを用いて、１回の地震における当該タンクの損傷度Ｄｊや、複数の地震発生による累積損傷度を算出することができる。これにより異なる大きさや、異なる対象を保管するタンクそれぞれの１回の地震における当該タンクの損傷度Ｄｊや、複数の地震発生による累積損傷度を、多くの時間やコスト、労力をかけずに算出することができる。
【0061】
情報処理装置１０の出力部１６は、上述した隅角部ひずみや、損傷度Ｄｊや、累積損傷度Ｄの数値を取得する。そして出力部１６は、そのタンク毎の隅角部ひずみや、損傷度Ｄｊや、累積損傷度Ｄの数値や、その数値に応じた表示態様の警告情報をタンク毎に示す表示情報を生成して、所定の出力装置に出力してよい。これにより、異なる大きさや、異なる対象を保管するタンクのうちの、どのタンクで、どの程度の隅角部ひずみが発生しているか、１度の地震でどの程度の損傷度Ｄｊが発生したか、過去の複数回の地震でどの程度の累積損傷度Ｄが発生しているかを確認することが可能となる。そしてこのような表示情報を、地震が終了した数分後に確認することで、タンクが受けた被害に対して迅速に復旧作業などの対処を行うことができる。
【0062】
図８は隅角部ひずみと浮き上がり変位との関係の比較例を示す第一の図である。
図９は隅角部ひずみと浮き上がり変位との関係の比較例を示す第二の図である。
図８と図９とを用いて異なるタンクについて、本開示による手法（隅角部ひずみ算出式（３）：簡易式）と、有限要素法（ＦＥＭ：Finite Element Method）のそれぞれで算出した隅角部ひずみと隅角部の浮き上がり変位との関係を比較する。図８で示す例はタンクＡ（ＴａｎｋＡ）について本開示による手法と、有限要素法のそれぞれで算出した隅角部ひずみと、隅角部の浮き上がり変位との関係を示す。図８で示すように本開示による手法と、有限要素法で算出した、隅角部ひずみと隅角部の浮き上がり変位との関係は、いずれもほぼ同じような結果となっている。
図９は隅角部ひずみと浮き上がり変位との関係の比較例を示す第二の図である。
図８と図９とを用いて異なるタンクについて、本開示による手法（隅角部ひずみ算出式（３）：簡易式）と、有限要素法（ＦＥＭ：Finite Element Method）のそれぞれで算出した隅角部ひずみと隅角部の浮き上がり変位との関係を比較する。図８で示す例はタンクＡ（ＴａｎｋＡ）について本開示による手法と、有限要素法のそれぞれで算出した隅角部ひずみと、隅角部の浮き上がり変位との関係を示す。図８で示すように本開示による手法と、有限要素法で算出した、隅角部ひずみと隅角部の浮き上がり変位との関係は、いずれもほぼ同じような結果となっている。
【0063】
図9で示す例はタンクＢ（ＴａｎｋＢ）について本開示による手法と、有限要素法のそれぞれで算出した隅角部ひずみと隅角部の浮き上がり変位との関係を示す。図９で示すように本開示による手法と、有限要素法で算出した、隅角部ひずみと隅角部の浮き上がり変位との関係は、いずれもほぼ同じような結果となっている。つまり、本開示の隅角部ひずみを算出するにあたり、隅角部ひずみ算出式（３）の簡易な式で有限要素法と同様の精度の高い隅角部ひずみを算出することが可能となる。
【0064】
図１０は情報処理装置の他の構成を示す図である。
図１１は他の構成による情報処理装置の処理フローを示す図である。
本開示の情報処理装置１０は少なくとも隅角部ひずみ算出部１３を備えればよい。
隅角部ひずみ算出部１３は、振動による、タンクを構成する側板と底板との接合部近傍を示す隅角部の浮き上がり変位波形のピークが示す値と、タンクに応じた係数とを隅角部ひずみ算出式に入力し、隅角部における隅角部ひずみを算出する（ステップＳ１１１）。
図１１は他の構成による情報処理装置の処理フローを示す図である。
本開示の情報処理装置１０は少なくとも隅角部ひずみ算出部１３を備えればよい。
隅角部ひずみ算出部１３は、振動による、タンクを構成する側板と底板との接合部近傍を示す隅角部の浮き上がり変位波形のピークが示す値と、タンクに応じた係数とを隅角部ひずみ算出式に入力し、隅角部における隅角部ひずみを算出する（ステップＳ１１１）。
【0065】
なお上述の情報処理装置１０の例では、タンクに与える振動の例が地震動である場合を例に説明した。しかしながら他の例では、タンクに与える振動の例が地震動以外であってもよい。例えば、風によって振動するタンクの隅角ひずみや、ひずみ振幅、損傷度、累積損傷度を算出するようにしてもよい。
【0066】
なお、上記の情報処理装置は、行う処理の全部または一部を実行するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより各部の処理を行う。ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ（Read Only Memory）、ＣＤ－ＲＯＭ（Compact Disc Read Only Memory）等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。また上記プログラムは、前述した機能の一部を実現するためのものであってもよく、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであってもよい。
【0067】
なお、上記の実施形態は下記のように定義されてもよい。
【0068】
（付記１）
振動による、タンクを構成する側板と底板との接合部近傍を示す隅角部の浮き上がり変位波形のピークが示す値と、前記タンクに応じた係数とを隅角部ひずみ算出式に入力し、前記隅角部における隅角部ひずみを算出する隅角部ひずみ算出部を備える情報処理システム。
振動による、タンクを構成する側板と底板との接合部近傍を示す隅角部の浮き上がり変位波形のピークが示す値と、前記タンクに応じた係数とを隅角部ひずみ算出式に入力し、前記隅角部における隅角部ひずみを算出する隅角部ひずみ算出部を備える情報処理システム。
【0069】
（付記２）
前記隅角部ひずみ算出部は、隅角部ひずみをε、前記底板を構成する鋼板の種別に応じた係数をＣ、前記タンクにおける底部の液体圧力をｐ、前記浮き上がり変位のピーク値をδ、前記底板の板厚をｔａ、前記底板の降伏応力をσｙとして、式（１）により表される前記隅角部ひずみ算出式を用いて、前記隅角部ひずみを算出する
付記１に記載の情報処理システム。
前記隅角部ひずみ算出部は、隅角部ひずみをε、前記底板を構成する鋼板の種別に応じた係数をＣ、前記タンクにおける底部の液体圧力をｐ、前記浮き上がり変位のピーク値をδ、前記底板の板厚をｔａ、前記底板の降伏応力をσｙとして、式（１）により表される前記隅角部ひずみ算出式を用いて、前記隅角部ひずみを算出する
【数11】
【0070】
（付記３）
前記隅角部ひずみ算出部は、前記振動を示す波形に基づいて前記タンクに生じた１回目の前記浮き上がり変位のピーク値に対応する第一の前記隅角部ひずみと、前記波形に基づいて前記タンクに生じた２回目以降の前記浮き上がり変位のピーク値それぞれに対応する第二の前記隅角部ひずみと、を算出し、
前記第一の前記隅角部ひずみの最大と最小との差を示すひずみ振幅を式（２）で示す第一ひずみ振幅算出式を用いて算出し、
前記第二の前記隅角部ひずみの最大と最小との差を示すひずみ振幅を式（３）で示す第二ひずみ振幅算出式を用いて算出する、
ひずみ振幅算出部とを備える付記２に記載の情報処理システム。
前記隅角部ひずみ算出部は、前記振動を示す波形に基づいて前記タンクに生じた１回目の前記浮き上がり変位のピーク値に対応する第一の前記隅角部ひずみと、前記波形に基づいて前記タンクに生じた２回目以降の前記浮き上がり変位のピーク値それぞれに対応する第二の前記隅角部ひずみと、を算出し、
前記第一の前記隅角部ひずみの最大と最小との差を示すひずみ振幅を式（２）で示す第一ひずみ振幅算出式を用いて算出し、
【数12】
【数13】
【0071】
（付記４）
前記波形に基づいて前記タンクに生じた前記浮き上がり変位のピークの繰り返し発生に応じた当該繰り返しの回数における前記隅角部の累積損傷度を、累積損傷度算出式を用いて算出する累積損傷度算出部と、
を備える付記３に記載の情報処理システム。
前記波形に基づいて前記タンクに生じた前記浮き上がり変位のピークの繰り返し発生に応じた当該繰り返しの回数における前記隅角部の累積損傷度を、累積損傷度算出式を用いて算出する累積損傷度算出部と、
を備える付記３に記載の情報処理システム。
【0072】
（付記５）
前記隅角部ひずみ算出部は、対象のタンクに応じた前記係数であって、前記底板を構成する鋼板の種別に応じた係数をＣ、前記タンクにおける底部の液体圧力をｐ、前記浮き上がり変位のピーク値をδ、前記底板の板厚をｔａ、前記底板の降伏応力をσｙと、前記隅角部ひずみ算出式を用いて、前記隅角部ひずみを算出する
付記２に記載の情報処理システム。
前記隅角部ひずみ算出部は、対象のタンクに応じた前記係数であって、前記底板を構成する鋼板の種別に応じた係数をＣ、前記タンクにおける底部の液体圧力をｐ、前記浮き上がり変位のピーク値をδ、前記底板の板厚をｔａ、前記底板の降伏応力をσｙと、前記隅角部ひずみ算出式を用いて、前記隅角部ひずみを算出する
付記２に記載の情報処理システム。
【0073】
（付記６）
発生した前記振動の波形と、前記タンクに生じた前記浮き上がり変位を１質点系非線形水平ばねモデルと、当該モデルから得られた水平方向の振動ピーク値と浮き上がり変位算出式とを用いて、前記浮き上がり変位のピーク値を算出する浮き上がり変位算出部と、
を備える付記１から付記５の何れか一つに記載の情報処理システム。
発生した前記振動の波形と、前記タンクに生じた前記浮き上がり変位を１質点系非線形水平ばねモデルと、当該モデルから得られた水平方向の振動ピーク値と浮き上がり変位算出式とを用いて、前記浮き上がり変位のピーク値を算出する浮き上がり変位算出部と、
を備える付記１から付記５の何れか一つに記載の情報処理システム。
【0074】
（付記７）
前記振動の波形として地震動の加速度波形を前記地震動が発生している際に取得する取得部と、
を備える付記６に記載の情報処理システム。
前記振動の波形として地震動の加速度波形を前記地震動が発生している際に取得する取得部と、
を備える付記６に記載の情報処理システム。
【0075】
（付記８）
振動による、タンクを構成する側板と底板との接合部近傍を示す隅角部の浮き上がり変位波形のピークが示す値と、前記タンクに応じた係数とを隅角部ひずみ算出式に入力し、前記隅角部における隅角部ひずみを算出する
情報処理方法。
振動による、タンクを構成する側板と底板との接合部近傍を示す隅角部の浮き上がり変位波形のピークが示す値と、前記タンクに応じた係数とを隅角部ひずみ算出式に入力し、前記隅角部における隅角部ひずみを算出する
情報処理方法。
【0076】
（付記９）
前記隅角部ひずみをε、前記底板を構成する鋼板の種別に応じた係数をＣ、前記タンクにおける底部の液体圧力をｐ、前記浮き上がり変位のピーク値をδ、前記底板の板厚をｔａ、前記底板の降伏応力をσｙとして、式（１）により表される前記隅角部ひずみ算出式を用いて、前記隅角部ひずみを算出する
付記８に記載の情報処理方法。
前記隅角部ひずみをε、前記底板を構成する鋼板の種別に応じた係数をＣ、前記タンクにおける底部の液体圧力をｐ、前記浮き上がり変位のピーク値をδ、前記底板の板厚をｔａ、前記底板の降伏応力をσｙとして、式（１）により表される前記隅角部ひずみ算出式を用いて、前記隅角部ひずみを算出する
【数14】
【0077】
（付記１０）
前記隅角部ひずみ算出部は、前記振動を示す波形に基づいて前記タンクに生じた１回目の前記浮き上がり変位のピーク値に対応する第一の前記隅角部ひずみと、前記波形に基づいて前記タンクに生じた２回目以降の前記浮き上がり変位のピーク値それぞれに対応する第二の前記隅角部ひずみと、を算出し、
前記第一の前記隅角部ひずみの最大と最小との差を示すひずみ振幅を式（２）で示す第一ひずみ振幅算出式を用いて算出し、
前記第二の前記隅角部ひずみの最大と最小との差を示すひずみ振幅を式（３）で示す第二ひずみ振幅算出式を用いて算出する、
付記９に記載の情報処理方法。
前記隅角部ひずみ算出部は、前記振動を示す波形に基づいて前記タンクに生じた１回目の前記浮き上がり変位のピーク値に対応する第一の前記隅角部ひずみと、前記波形に基づいて前記タンクに生じた２回目以降の前記浮き上がり変位のピーク値それぞれに対応する第二の前記隅角部ひずみと、を算出し、
前記第一の前記隅角部ひずみの最大と最小との差を示すひずみ振幅を式（２）で示す第一ひずみ振幅算出式を用いて算出し、
【数15】
【数16】
【0078】
（付記１１）
前記波形に基づいて前記タンクに生じた前記浮き上がり変位のピークの繰り返し発生に応じた当該繰り返しの回数における前記隅角部の累積損傷度を、累積損傷度算出式を用いて算出する
を備える付記１０に記載の情報処理方法。
前記波形に基づいて前記タンクに生じた前記浮き上がり変位のピークの繰り返し発生に応じた当該繰り返しの回数における前記隅角部の累積損傷度を、累積損傷度算出式を用いて算出する
を備える付記１０に記載の情報処理方法。
【0079】
（付記１２）
対象のタンクに応じた前記係数であって、前記底板を構成する鋼板の種別に応じた係数をＣ、前記タンクにおける底部の液体圧力をｐ、前記浮き上がり変位のピーク値をδ、前記底板の板厚をｔａ、前記底板の降伏応力をσｙと、前記隅角部ひずみ算出式を用いて、前記隅角部ひずみを算出する
付記９に記載の情報処理方法。
対象のタンクに応じた前記係数であって、前記底板を構成する鋼板の種別に応じた係数をＣ、前記タンクにおける底部の液体圧力をｐ、前記浮き上がり変位のピーク値をδ、前記底板の板厚をｔａ、前記底板の降伏応力をσｙと、前記隅角部ひずみ算出式を用いて、前記隅角部ひずみを算出する
付記９に記載の情報処理方法。
【0080】
（付記１３）
発生した前記振動の波形と、前記タンクに生じた前記浮き上がり変位を１質点系非線形水平ばねモデルと、当該モデルから得られた水平方向の振動ピーク値と浮き上がり変位算出式とを用いて、前記浮き上がり変位のピーク値を算出する
付記８から付記１２の何れか一つに記載の情報処理方法。
発生した前記振動の波形と、前記タンクに生じた前記浮き上がり変位を１質点系非線形水平ばねモデルと、当該モデルから得られた水平方向の振動ピーク値と浮き上がり変位算出式とを用いて、前記浮き上がり変位のピーク値を算出する
付記８から付記１２の何れか一つに記載の情報処理方法。
【0081】
（付記１４）
前記振動の波形として地震動の加速度波形を前記地震動が発生している際に取得する
付記１３に記載の情報処理方法。
前記振動の波形として地震動の加速度波形を前記地震動が発生している際に取得する
付記１３に記載の情報処理方法。
【0082】
（付記１５）
情報処理システムのコンピュータを、
振動による、タンクを構成する側板と底板との接合部近傍を示す隅角部の浮き上がり変位波形のピークが示す値と、前記タンクに応じた係数とを隅角部ひずみ算出式に入力し、前記隅角部における隅角部ひずみを算出する隅角部ひずみ算出手段、
として機能させるプログラム。
情報処理システムのコンピュータを、
振動による、タンクを構成する側板と底板との接合部近傍を示す隅角部の浮き上がり変位波形のピークが示す値と、前記タンクに応じた係数とを隅角部ひずみ算出式に入力し、前記隅角部における隅角部ひずみを算出する隅角部ひずみ算出手段、
として機能させるプログラム。
【0083】
（付記１６）
前記隅角部ひずみ算出手段は、隅角部ひずみをε、前記底板を構成する鋼板の種別に応じた係数をＣ、前記タンクにおける底部の液体圧力をｐ、前記浮き上がり変位のピーク値をδ、前記底板の板厚をｔａ、前記底板の降伏応力をσｙとして、式（１）により表される前記隅角部ひずみ算出式を用いて、前記隅角部ひずみを算出する
付記１５に記載のプログラム。
前記隅角部ひずみ算出手段は、隅角部ひずみをε、前記底板を構成する鋼板の種別に応じた係数をＣ、前記タンクにおける底部の液体圧力をｐ、前記浮き上がり変位のピーク値をδ、前記底板の板厚をｔａ、前記底板の降伏応力をσｙとして、式（１）により表される前記隅角部ひずみ算出式を用いて、前記隅角部ひずみを算出する
【数17】
【0084】
（付記１７）
前記隅角部ひずみ算出手段は、前記振動を示す波形に基づいて前記タンクに生じた１回目の前記浮き上がり変位のピーク値に対応する第一の前記隅角部ひずみと、前記波形に基づいて前記タンクに生じた２回目以降の前記浮き上がり変位のピーク値それぞれに対応する第二の前記隅角部ひずみと、を算出し、
前記第一の前記隅角部ひずみの最大と最小との差を示すひずみ振幅を式（２）で示す第一ひずみ振幅算出式を用いて算出し、
前記第二の前記隅角部ひずみの最大と最小との差を示すひずみ振幅を式（３）で示す第二ひずみ振幅算出式を用いて算出する、
ひずみ振幅算出手段、
として機能させる付記１６に記載のプログラム。
前記隅角部ひずみ算出手段は、前記振動を示す波形に基づいて前記タンクに生じた１回目の前記浮き上がり変位のピーク値に対応する第一の前記隅角部ひずみと、前記波形に基づいて前記タンクに生じた２回目以降の前記浮き上がり変位のピーク値それぞれに対応する第二の前記隅角部ひずみと、を算出し、
前記第一の前記隅角部ひずみの最大と最小との差を示すひずみ振幅を式（２）で示す第一ひずみ振幅算出式を用いて算出し、
【数18】
【数19】
として機能させる付記１６に記載のプログラム。
【0085】
（付記１８）
前記波形に基づいて前記タンクに生じた前記浮き上がり変位のピークの繰り返し発生に応じた当該繰り返しの回数における前記隅角部の累積損傷度を、累積損傷度算出式を用いて算出する累積損傷度算出手段、
として機能させる付記１７に記載のプログラム。
前記波形に基づいて前記タンクに生じた前記浮き上がり変位のピークの繰り返し発生に応じた当該繰り返しの回数における前記隅角部の累積損傷度を、累積損傷度算出式を用いて算出する累積損傷度算出手段、
として機能させる付記１７に記載のプログラム。
【0086】
（付記１９）
前記隅角部ひずみ算出手段は、対象のタンクに応じた前記係数であって、前記底板を構成する鋼板の種別に応じた係数をＣ、前記タンクにおける底部の液体圧力をｐ、前記浮き上がり変位のピーク値をδ、前記底板の板厚をｔａ、前記底板の降伏応力をσｙと、前記隅角部ひずみ算出式を用いて、前記隅角部ひずみを算出する
付記１６に記載のプログラム。
前記隅角部ひずみ算出手段は、対象のタンクに応じた前記係数であって、前記底板を構成する鋼板の種別に応じた係数をＣ、前記タンクにおける底部の液体圧力をｐ、前記浮き上がり変位のピーク値をδ、前記底板の板厚をｔａ、前記底板の降伏応力をσｙと、前記隅角部ひずみ算出式を用いて、前記隅角部ひずみを算出する
付記１６に記載のプログラム。
【0087】
（付記２０）
発生した前記振動の波形と、前記タンクに生じた前記浮き上がり変位を１質点系非線形水平ばねモデルと、当該モデルから得られた水平方向の振動ピーク値と浮き上がり変位算出式とを用いて、前記浮き上がり変位のピーク値を算出する浮き上がり変位算出手段、
として機能させる付記１５から付記１９の何れか一つに記載のプログラム。
発生した前記振動の波形と、前記タンクに生じた前記浮き上がり変位を１質点系非線形水平ばねモデルと、当該モデルから得られた水平方向の振動ピーク値と浮き上がり変位算出式とを用いて、前記浮き上がり変位のピーク値を算出する浮き上がり変位算出手段、
として機能させる付記１５から付記１９の何れか一つに記載のプログラム。
【0088】
（付記２１）
前記振動の波形として地震動の加速度波形を前記地震動が発生している際に取得する取得手段、
として機能させる付記２０に記載のプログラム。
前記振動の波形として地震動の加速度波形を前記地震動が発生している際に取得する取得手段、
として機能させる付記２０に記載のプログラム。
【符号の説明】
【0089】
１０・・・情報処理装置
２０・・・タンク毎地震動ＤＢ
３０・・・タンク算出結果ＤＢ
１１・・・取得部
１２・・・浮き上がり変位算出部
１３・・・隅角ひずみ算出部
１４・・・ひずみ振幅算出部
１５・・・累積損傷度算出部
１６・・・出力部
２０・・・タンク毎地震動ＤＢ
３０・・・タンク算出結果ＤＢ
１１・・・取得部
１２・・・浮き上がり変位算出部
１３・・・隅角ひずみ算出部
１４・・・ひずみ振幅算出部
１５・・・累積損傷度算出部
１６・・・出力部
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】