JP6179092B2 - 線状構造物の渦励振流体力データ計算プログラム、渦励振を考慮した線状構造物の設計プログラム及び線状構造物の渦励振流体力実験方法 - Google Patents

Description

本発明は、流体中に設けられた線状構造物に発生する渦励振（ＶＩＶ：Vortex-Induced Vibration）に関連した線状構造物の渦励振流体力データ計算プログラム、渦励振を考慮した線状構造物の設計プログラム及び線状構造物の渦励振流体力実験方法に関する。

一様な流体中におかれたリジッドな円柱等の物体の下流側には交互にカルマン渦が発生する。円柱の片側が固定されもう一方がフリーの柱状体や長尺ものの両側がフリーとみなせる円形線状構造物の場合、カルマン渦の周波数と物体の固有振動数がほぼ一致すると、流体力の変化が励振を誘起し、渦励振が物体の流れに直交する方向に発生する。海底での掘削や資源の採取を行う際にライザー等の海中線状構造物を使用するが、このとき、潮流等によって渦励振が発生し、最悪の場合には材質疲労によって線状構造物が破損するおそれがある。また、空中に設けられたポールや電線等の片持ち状態の柱状体には風の影響により柱状体が振動する渦励振が発生する。これに対して、柱状体の渦励振を防止する方法や構造が開示されている（特許文献１及び２）。

特開昭６１−３８２４０号公報 特開２００４−３６７９４号公報

特許文献１や特許文献２の技術は、片側がフリーの柱状体やポール等の構造的な渦励振対策であり、渦励振の予測技術ではない。特に、両側が擬似的にフリーな海中線状構造物であるライザー等は、線状構造物の破損等を防止するために潮流中での線状構造物の挙動予測や疲労寿命評価を適切に行うことが要求されているが、その予測技術は確立されていない。

本発明は、流体中に設けられた線状構造物の挙動予測や疲労寿命評価につながる流体力に関する予測技術を提供することを目的とする。

本発明の線状構造物の渦励振流体力データ計算プログラムは、コンピュータによって、線状構造物である処理対象に働く流体力データを求めるプログラムであって、前記処理対象の構造物条件と、前記処理対象への加振条件と、レイノルズ数を演算条件として取得するステップ１と、線状構造物に対する実験において、加振条件とレイノルズ数との組み合わせに対して得られた実験流体力データを取得するステップ２と、前記実験流体力データを補間して前記演算条件に含まれるレイノルズ数及び加振条件における補間流体力データを得るステップ３と、前記補間流体力データに対し、前記演算条件に含まれる構造物条件に対応する流体力データを計算流体力データとして求めるステップ４と、前記計算流体力データを出力するステップ５を行う。

なお、線状構造物には、海中に設けられるライザーや浮力缶等、また、空中に設けられるポールや電線等、流体中に設けられる線状を成した構造物のすべてを含む。

ここで、前記演算条件に含まれる構造物条件は、前記処理対象に設けられるストレーキの高さ条件とピッチ条件を含み、前記ステップ４は、前記高さ条件に対応する流体力データを求めるステップ４−１と、前記ピッチ条件に対応する流体力データを求めるステップ４−２を有することが好ましい。なお、流体力データを求める順番は、ステップ４−２が先でステップ４−１が後であってもよい。

また、前記演算条件に含まれる構造物条件は、前記処理対象への流れの流入角度条件を含み、前記ステップ４は、前記流入角度条件に対応する流体力データを求めるステップ４−３を有することが好ましい。なお、ステップ４におけるステップ４−３，ステップ４−１，ステップ４−２は任意のステップを選択でき、流体力データを求める順も任意の順が選択できる。

前記演算条件に含まれる構造物条件は、前記処理対象に設けられるストレーキの高さ条件とピッチ条件、前記処理対象への流れの流入角度条件を含み、前記ステップ４は、前記高さ条件に対応する流体力データを求めるステップ４−１と、前記ピッチ条件に対応する流体力データを求めるステップ４−２と、前記流入角度条件に対応する流体力データを求めるステップ４−３を順に行うことが好適である。

前記実験流体力データは、前記実験で用いられた線状構造物に設けられたストレーキの高さとピッチ及び／又は前記実験で用いられた線状構造物への流れの流入角度の組み合わせ毎の離散型データであることが好適である。特に、前記離散型データは、複数のレイノルズ数毎に実験で得られた流体力のデータであることが好適である。

また、前記加振条件は、前記処理対象に対する加振周波数及び振幅であることが好ましい。なお、加振時の波形は、正弦波又は擬似正弦波であることが好ましい。

また、前記実験流体力データは、実験抗力係数、実験揚力方向付加質量係数及び実験揚力方向減衰係数を含み、前記ステップ４は、前記計算流体力データとして計算抗力係数、計算揚力方向付加質量係数及び計算揚力方向減衰係数を求めることが好適である。

また、前記ステップ４−１は、前記離散型データを補間して前記ストレーキの高さ条件に対応する流体力データを求めることが好適である。また、前記ステップ４−２は、前記離散型データを補間して前記ストレーキのピッチ条件に対応する流体力データを求めることが好適である。また、前記ステップ４−３は、前記離散型データを補間して前記流入角度条件に対応する流体力データを求めることが好適である。

前記線状構造物が円形断面を有するストレーキ付き円柱形状である場合、前記実験流体力データは、第１範囲：（Ｈ／Ｄ，Ｐ／Ｄ）＝（０，０），Ｒｅ＝１．５×１０⁵〜５．０×１０⁵，ｆＤ／Ｖ＝０．０１〜０．３，Ａ／Ｄ＝０．１〜１．１、第２範囲：（Ｈ／Ｄ，Ｐ／Ｄ）＝（０．１５，５），Ｒｅ＝１．３×１０⁵〜５．２×１０⁵，ｆＤ／Ｖ＝０．０５〜０．３，Ａ／Ｄ＝０．１〜０．９、第３範囲：（Ｈ／Ｄ，Ｐ／Ｄ）＝（０．１５，１５），Ｒｅ＝１．３×１０⁵〜５．２×１０⁵，ｆＤ／Ｖ＝０．０５〜０．３，Ａ／Ｄ＝０．１〜０．９、第４範囲：（Ｈ／Ｄ，Ｐ／Ｄ）＝（０．２５，１５），Ｒｅ＝１．５×１０⁵〜５．０×１０⁵，ｆＤ／Ｖ＝０．０５〜０．３，Ａ／Ｄ＝０．１〜０．９とし、前記第１範囲から前記第４範囲の少なくとも１つの範囲で取得されることが好ましい。ただし、Ｈ：実験での線状構造物のストレーキの高さ（ｍ）、Ｄ：実験での線状構造物の直径（ｍ）、Ｐ：実験での線状構造物のストレーキの螺旋ピッチ（ｍ）、Ｒｅ：実験時のレイノルズ数、ｆＤ／Ｖ：実験無次元加振周波数、ｆ：実験加振周波数（Ｈｚ）、Ｖ：実験における線状構造物への流れの流速（ｍ／ｓ）、Ａ／Ｄ：実験無次元振幅、Ａ：実験での振動の振幅（ｍ）である。

前記線状構造物が円形断面を有するストレーキ付き円柱形状である場合、前記実験流体力データは、第５範囲：χ＝０〜４５ｄｅｇ，Ｒｅ＝１．５×１０⁵〜５．０×１０⁵，ｆＤ／Ｖ＝０．０１〜０．３，Ａ／Ｄ＝０．１〜１．１とし、前記第５範囲の範囲で取得されることが好適である。ただし、χ：実験における線状構造物への流れの流入角度（線状構造物への流れが直角のときを０とする）、Ｄ：実験での線状構造物の直径（ｍ）、Ｒｅ：実験時のレイノルズ数、ｆＤ／Ｖ：実験無次元加振周波数、ｆ：実験加振周波数（Ｈｚ）、Ｖ：実験における線状構造物への流れの流速（ｍ／ｓ）、Ａ／Ｄ：実験無次元振幅、Ａ：実験での振動の振幅（ｍ）である。

本発明における渦励振を考慮した線状構造物の設計プログラムは、コンピュータによって、上記渦励振流体力データ計算プログラムを用いて得られた前記計算流体力データを含む設定条件を取得するステップ６と、前記設定条件に基づいて前記処理対象の静的な釣り合い条件を求めるステップ７と、前記釣り合い条件及び前記計算流体力データを用いて設定振動周波数における前記処理対象に働く流体力を計算するステップ８とを行い、前記設定振動周波数を変化させつつ前記ステップ８を繰り返して、前記処理対象に働く流体力を計算して前記処理対象が前記釣り合い条件に達したかを判定する。

なお、前記ステップ８で得られる結果には、線状構造物に働く流体力そのもの以外にも、線状構造物の応答や線状構造物の内力等が含まれる。

ここで、前記ステップ８の繰り返し処理において、前記渦励振流体力データ計算プログラムを用いて前記計算流体力データを再計算し、再計算された前記計算流体力データを用いて前記処理対象が前記釣り合い条件に達したかを判定する。

前記ステップ８は、前記処理対象が置かれた潮流条件を用いて前記処理対象に働く流体力を計算することが好適である。なお、潮流条件は、実海域を想定した潮流条件であることが好適である。また、前記ステップ８は、前記処理対象の固有振動数を算出し、前記処理対象に働く流体力を計算することが好適である。

さらに、前記処理対象に働く流体力を用いて前記処理対象の疲労評価を行うステップ９を備えてもよい。この疲労評価を行うステップ９を備えることにより、線状構造物の疲労寿命評価が可能となる。

また、前記ステップ８における流体力の計算は、有限要素法を用いて行われることが好適である。この有限要素法を用いることにより、構造力学分野で広く使われている数値解析手法が利用できる。

本発明における渦励振流体力実験方法は、上記渦励振流体力データ計算プログラムに用いる前記実験流体力データを得るための渦励振流体力実験方法であって、線状構造物の部分模型を強制揺動させて前記実験流体力データを取得する。

ここで、前記部分模型は、実物大模型であることが好適である。実物大模型を用いることにより、より正確に実験流体力データを取得することが可能となる。

本発明の線状構造物の渦励振流体力データ計算プログラムは、コンピュータにより、線状構造物である処理対象に働く流体力データを求めるプログラムであって、前記処理対象の構造物条件と、前記処理対象への加振条件と、レイノルズ数を演算条件として取得するステップ１と、線状構造物に対する実験において、加振条件とレイノルズ数との組み合わせに対して得られた実験流体力データを取得するステップ２と、前記実験流体力データを補間して前記演算条件に含まれるレイノルズ数及び加振条件における補間流体力データを得るステップ３と、前記補間流体力データに対し、前記演算条件に含まれる構造物条件に対応する流体力データを計算流体力データとして求めるステップ４と、前記計算流体力データを出力するステップ５を実行することによって、流体の流れの中での線状構造物に働く流体力を評価するために必要な計算流体力データを求めることができる。このとき、後述する渦励振流体力実験方法によって線状構造物に関する実験流体力データを用いて計算流体力データを求めるので、設計対象となる線状構造物に働く渦励振を考慮した流体力を正確に求めることができる。

ここで、前記演算条件に含まれる構造物条件は、前記処理対象に設けられるストレーキの高さ条件とピッチ条件を含み、前記ステップ４は、前記高さ条件に対応する流体力データを求めるステップ４−１と、前記ピッチ条件に対応する流体力データを求めるステップ４−２を有することによって、ストレーキを備えた線状構造物の計算流体力データを求めることができる。

また、前記演算条件に含まれる構造物条件は、前記処理対象への流れの流入角度条件を含み、前記ステップ４は、前記流入角度条件に対応する流体力データを求めるステップ４−３を有することによって、流体の流れに対して傾きを有する線状構造物の計算流体力データを求めることができる。

また、前記演算条件に含まれる構造物条件は、前記処理対象に設けられるストレーキの高さ条件とピッチ条件、前記処理対象への流れの流入角度条件を含み、前記ステップ４は、前記高さ条件に対応する流体力データを求めるステップ４−１と、前記ピッチ条件に対応する流体力データを求めるステップ４−２と、前記流入角度条件に対応する流体力データを求めるステップ４−３を順に行うことによって、線状構造物への影響が大きい順に流体力データを評価することができ、より正確な計算流体力データを得ることができる。

また、前記実験流体力データは、前記実験で用いられた線状構造物に設けられたストレーキの高さとピッチ及び／又は前記実験で用いられた線状構造物への流れの流入角度の組み合わせ毎の離散型データであることによって、少ない実験条件における実験流体力データを補間処理して計算流体力データを得ることができる。特に、前記離散型データは、複数のレイノルズ数毎に実験で得られた流体力のデータであることによって、レイノルズ数の変化に応じた計算流体力データの評価が可能となる。

前記加振条件は、前記処理対象に対する加振周波数及び振幅であることによって、渦励振流体力実験方法において線状構造物への流体の流れによる渦励振の影響を評価することができる。

前記実験流体力データは、実験抗力係数、実験揚力方向付加質量係数及び実験揚力方向減衰係数を含み、前記ステップ４は、前記計算流体力データとして計算抗力係数、計算揚力方向付加質量係数及び計算揚力方向減衰係数を求めることによって、計算抗力係数、計算揚力方向付加質量係数及び計算揚力方向減衰係数を用いての線状構造物への流体力の影響評価が可能となる。

また、前記ステップ４−１は、前記離散型データを補間して前記ストレーキの高さ条件に対応する流体力データを求めること、又は、前記ステップ４−２は、前記離散型データを補間して前記ストレーキのピッチ条件に対応する流体力データを求めること、又は、前記ステップ４−３は、前記離散型データを補間して前記流入角度条件に対応する流体力データを求めることによって、前記ストレーキの高さ、前記ストレーキのピッチ条件及び前記流入角度条件の少なくとも１つの計算流体力データの評価を行うことができる。

また、前記線状構造物が円形断面を有するストレーキ付き円柱形状である場合、前記実験流体力データは、第１範囲：（Ｈ／Ｄ，Ｐ／Ｄ）＝（０，０），Ｒｅ＝１．５×１０⁵〜５．０×１０⁵，ｆＤ／Ｖ＝０．０１〜０．３，Ａ／Ｄ＝０．１〜１．１、第２範囲：（Ｈ／Ｄ，Ｐ／Ｄ）＝（０．１５，５），Ｒｅ＝１．３×１０⁵〜５．２×１０⁵，ｆＤ／Ｖ＝０．０５〜０．３，Ａ／Ｄ＝０．１〜０．９、第３範囲：（Ｈ／Ｄ，Ｐ／Ｄ）＝（０．１５，１５），Ｒｅ＝１．３×１０⁵〜５．２×１０⁵，ｆＤ／Ｖ＝０．０５〜０．３，Ａ／Ｄ＝０．１〜０．９、第４範囲：（Ｈ／Ｄ，Ｐ／Ｄ）＝（０．２５，１５），Ｒｅ＝１．５×１０⁵〜５．０×１０⁵，ｆＤ／Ｖ＝０．０５〜０．３，Ａ／Ｄ＝０．１〜０．９とし、前記第１範囲から前記第４範囲の少なくとも１つの範囲で取得することによって、例えば、海底での掘削や資源の採取を行うライザーを評価対象とした場合のストレーキの高さやピッチの変化に応じた計算流体力データをより正確に得ることができる。

また、前記線状構造物が円形断面を有するストレーキ付き円柱形状である場合、前記実験流体力データは、第５範囲：χ＝０〜４５ｄｅｇ，Ｒｅ＝１．５×１０⁵〜５．０×１０⁵，ｆＤ／Ｖ＝０．０１〜０．３，Ａ／Ｄ＝０．１〜１．１とし、前記第５範囲の範囲で取得することによって、例えば、海底での掘削や資源の採取を行うライザーを評価対象とした場合の流体の流入角度の変化に応じた計算流体力データをより正確に得ることができる。

本発明における渦励振を考慮した線状構造物の設計プログラムは、コンピュータによって、上記渦励振流体力データ計算プログラムを用いて得られた前記計算流体力データを含む設定条件を取得するステップ６と、前記設定条件に基づいて前記処理対象の静的な釣り合い条件を求めるステップ７と、前記釣り合い条件及び前記計算流体力データを用いて設定振動周波数における前記処理対象に働く流体力を計算するステップ８と、前記設定振動周波数を変化させつつ前記ステップ８を繰り返して、前記処理対象に働く流体力を計算して前記処理対象が前記釣り合い条件に達したかを判定することによって、上記渦励振流体力データ計算プログラムによって得られた計算流体力データを用いて線状構造物へ働く渦励振を考慮した流体力を正確に評価することができる。これにより、渦励振対策を図った線状構造物の設計を容易にすることができる。

ここで、前記ステップ８の繰り返し処理において、前記渦励振流体力データ計算プログラムを用いて前記計算流体力データを再計算し、再計算された前記計算流体力データを用いて前記処理対象が前記釣り合い条件に達したかを判定することによって、上記渦励振流体力データ計算プログラムによって再計算された計算流体力データを用いて線状構造物へ働く流体力をより正確に評価することができる。

さらに、前記処理対象に働く流体力を用いて前記処理対象の疲労評価を行うことによって、線状構造物の流体環境下における疲労寿命評価を行って設計を行うことができる。

本発明における渦励振流体力実験方法は、上記渦励振流体力データ計算プログラムに用いる前記実験流体力データを得るための渦励振流体力実験方法であって、線状構造物の部分模型を強制揺動させて前記実験流体力データを取得することによって、例えば、海底での掘削や資源の採取を行うライザーを評価対象とした場合のストレーキの高さやピッチの変化に応じた実験流体力データを正確に得ることができる。このとき、前記部分模型は実物大模型とすることによって、実物に則した実験流体力データをより正確に得ることができる。

本発明の実施の形態における線形構造物の部分模型を示す斜視図である。 本発明の実施の形態における線形構造物の部分模型の断面図である。 本発明の実施の形態における線形構造物の部分模型への流体の流入角度を説明する図である。 本発明の実施の形態における渦励振流体力データ計算処理及び線状構造物の設計処理を行う装置の構成を示す図である。 本発明の実施の形態における渦励振流体力データ計算処理のフローチャートである。 本発明の実施の形態における渦励振流体力データ計算処理のサブルーチンのフローチャートである。 本発明の実施の形態における線状構造物の設計処理のフローチャートである。

＜渦励振流体力実験方法＞
本実施の形態における線状構造物の渦励振流体力データ計算プログラム、渦励振を考慮した設計プログラムでは、線状構造物の模型に対する実験において得られた実験流体力データを用いて処理対象となる線状構造物に働く流体力等の評価を行う。そこで、まず線状構造物の模型に対する実験による渦励振流体力実験方法について説明する。

本実施の形態では、より正確な実験流体力データを得るために実機と同等の直径を有する線状構造物（ライザー）の部分模型を使用し、曳航中強制動揺実験により実機相当流体力を直接求める方策を採用した。線状構造物は、裸殻状態とＶＩＶを低減させる１つの方法である螺旋ストレーキ（以下，ストレーキという）を取り付けた場合の流体力を取得した。

図１は、実験に使用した線状構造物の部分模型１００の斜視図を示す。部分模型１００は、ライザーの管１２と管１２に取り付けられたストレーキ１０と両端部に設けた端板１４から構成される。曳航中強制動揺実験は、水槽（長さ１５０ｍ、幅７．５ｍ、水深３．５ｍ）で実験を行った。水槽の曳引台車に強制動揺装置を設置し潮流を模擬した前進速度を与えながら、所定の加振条件（加振周波数及び加振振幅）で部分模型１００を上下動させる。部分模型１００の両端部に上下方向と水平方向の力を計測するための２分力計をそれぞれ取り付け、部分模型１００は２分力計を介して端板１４に取り付けることによって計測を行った。端板１４は、強制動揺実験において部分模型１００の端部における３次元的な影響を除くために設けられる。また、部分模型１００の長さは、３次元的な影響を低減するため、直径の５倍以上の長さとした。

部分模型１００は、ストレーキ１０の影響を調査する場合には、裸殻の模型に硬質ゴムスポンジ製のストレーキ１０を取り付けて実験を行う。ストレーキ１０は、図２に示すように、管１２の直径（Ｄ）とストレーキ１０の高さ（Ｈ）との比（Ｈ／Ｄ）、及び、管１２の直径（Ｄ）とストレーキ１０の螺旋ピッチ（Ｐ）との比（Ｐ／Ｄ）を変更して部分模型１００に働く流体力の変化について計測する。なお、実験に使用するストレーキの材料については、硬質ゴムスポンジ以外に、金属、木材、高分子材、セラミック等の各種の材料が選択し得る。

また、図３（ｂ）に示すように、部分模型１００として流入角度の変更対応品を準備し、流れの方向（流入角度χ）を変更して部分模型１００に働く流体力の変化を計測することを行った。ここで、流入角度χとは、図３（ａ）に示すように、部分模型１００の軸方向に垂直な流入角度を０としたときの傾きを意味する。部分模型１００には、部分模型１００の軸方向に垂直な流れの流速成分Ｖ_nが作用することになる。

計測で得られる実験流体力データは、抗力係数（実験抗力係数）、揚力方向付加質量係数（実験揚力方向付加質量係数）、揚力方向減衰係数（実験揚力方向減衰係数）を含んでいる。すなわち、部分模型１００に作用する抗力Ｆ_Dの時系列データ平均値、Ｆ_DMから（数式１）の実験抗力係数Ｃ_Dを算出する。

ここで，ρは水の密度、Ｖ_nは部分模型１００の軸方向に垂直な方向の流れの流速成分である。

また、揚力方向の流体力変動を加振変位と同位相成分の力Ｆ_L同位相（付加質量成分）と、加振変位と９０°位相が異なる成分の力Ｆ_L90°（減衰成分）に分離し、実験揚力方向付加質量係数（Ｃ_LA）と実験減衰係数（Ｃ_LD）を（数式２）、（数式３）から求める．

ここで，Ｍ：ストレーキ１０付きの管１２の質量（空中質量）、Ａ：実験振動振幅、ｆ：実験加振周波数である。

例えば、部分模型１００は、ストレーキ１０を設けていない裸殻状態の管１２において直径（Ｄ）を０．３２０ｍ及び長さ（Ｌ）を２．１６ｍの鋼管とした。なお、ストレーキ１０の管１２に対する巻き数は３とした。また、ストレーキ１０自体の形状は二等辺三角形とし、その底辺幅は断面において５０ｍｍの一定とした。また、管１２の両端には、端部の３次元影響を除くための端板１４を取り付けた。ただし、部分模型１００はこれに限定されるものでなく、実際に使用される線状構造物の一部を模したものとすることが好ましい。なお、ストレーキ１０は、付いたもの、付かないもの、ストレーキ形状も螺旋形以外の例えば直線形のもの、独立型のもの、連続型のもの等、任意の形状が選択し得る。また、管１２も円筒以外の多角管や楕円管等の任意の形状が選択できる。

部分模型１００を用いて実験を行う際に、レイノルズ数Ｒｅ（＝ＶＤ／ν、ν：動粘性係数）は１．０×１０⁵〜５．２×１０⁵の範囲で複数の点を選んだ。また、加振条件として、無次元加振周波数ｆＤ／Ｖ（ｆ：実験加振周波数、Ｖ：流速）は０．０１〜０．３の範囲、無次元振幅Ａ／Ｄ（Ａ：構造物振動振幅）は０．１〜１．１の範囲で複数の条件を選んで実験を行った。また、部分模型１００の構造物条件として、ストレーキ１０は、管１２の直径（Ｄ）とストレーキ１０の高さ（Ｈ）との比（Ｈ／Ｄ）で最大０．２５まで，管１２の直径（Ｄ）とストレーキ１０の螺旋ピッチ（Ｐ）との比（Ｐ／Ｄ）は最大１５までの範囲で複数条件を設定した。なお、（Ｈ／Ｄ）＝０及び（Ｐ／Ｄ）＝０は、部分模型１００にストレーキ１０を設けていないことを示す。

特に、海中線状構造物としてライザーを想定した渦励起流体力データ計算プログラムに用いる実験流体力データを取得する場合、第１範囲：（Ｈ／Ｄ，Ｐ／Ｄ）＝（０，０），Ｒｅ＝１．５×１０⁵〜５．０×１０⁵，ｆＤ／Ｖ＝０．０１〜０．３，Ａ／Ｄ＝０．１〜１．１、第２範囲：（Ｈ／Ｄ，Ｐ／Ｄ）＝（０．１５，５），Ｒｅ＝１．３×１０⁵〜５．２×１０⁵，ｆＤ／Ｖ＝０．０５〜０．３，Ａ／Ｄ＝０．１〜０．９、第３範囲：（Ｈ／Ｄ，Ｐ／Ｄ）＝（０．１５，１５），Ｒｅ＝１．３×１０⁵〜５．２×１０⁵，ｆＤ／Ｖ＝０．０５〜０．３，Ａ／Ｄ＝０．１〜０．９、第４範囲：（Ｈ／Ｄ，Ｐ／Ｄ）＝（０．２５，１５），Ｒｅ＝１．５×１０⁵〜５．０×１０⁵，ｆＤ／Ｖ＝０．０５〜０．３，Ａ／Ｄ＝０．１〜０．９の範囲において条件を設定して実験流体力データを取得することが好ましい。

また、流入角度χの条件については、第５範囲：χ＝０〜４５ｄｅｇ，Ｒｅ＝１．５×１０⁵〜５．０×１０⁵，ｆＤ／Ｖ＝０．０１〜０．３，Ａ／Ｄ＝０．１〜１．１において条件を設定して実験流体力データを取得することが好ましい。このとき、（Ｈ／Ｄ，Ｐ／Ｄ）＝（０，０）、すなわち部分模型１００にストレーキ１０を設けない状態で部分模型１００に対する流れの方向（流入角度χ）を変更して部分模型１００に働く流体力の変化を実験流体力データとすることも可能である。

このような範囲で各パラメータの組み合わせを変更しつつ実験条件を設定し、部分模型１００を強制的に揺動させて各パラメータの組み合わせに対応する離散的な実験流体力データを取得した。

＜渦励振流体力データ計算処理＞
渦励振流体力実験において取得された実験流体力データを用いて渦励振流体力データを求める渦励振流体力データ計算処理について説明する。渦励振流体力データ計算処理では、処理対象となる線状構造物の構造物条件、加振条件及びレイノルズ数に基づいて、処理対象に生ずる渦励振流体力を解析するための計算流体力データを求める。

渦励振流体力データ計算処理は、渦励振流体力データ計算プログラムをコンピュータで実行することによって実現される。図４は、渦励振流体力データ計算プログラムを実行するためのコンピュータ２００の構成を示す図である。コンピュータ２００は、処理部２０、記憶部２２、入力部２４及び出力部２６を含んで構成される。処理部２０は、ＣＰＵ等から構成され、記憶部２２に予め記憶された渦励振流体力データ計算プログラムを実行することによって本実施の形態における渦励振流体力データ計算処理を実現する。記憶部２２は、処理部２０からアクセス可能な半導体メモリ、ハードディスク、光ディスク等の記憶装置を含んで構成される。記憶部２２は、渦励振流体力データ計算プログラム、実験流体力データを含む実験流体力データベース、入力条件や各ステップにおける計算結果等、渦励振流体力データ計算処理に必要な電子情報を記憶する。入力部２４は、キーボード、マウス等の入力装置、インターネット等の電子情報回線からデータを取得するインターフェース等を含んで構成される。入力部２４は、渦励振流体力データ計算処理に必要な電子情報の入力を受けて、記憶部２２に記憶させる。出力部２６は、ディスプレイやプリンタ等の出力装置、インターネット等の電子情報回線へデータを出力するインターフェース等を含んで構成される。出力部２６は、渦励振流体力データ計算処理で得られた演算結果をユーザや外部装置等に出力するために用いられる。

渦励振流体力データ計算処理は、図５のフローチャートに沿って行われる。当該フローチャートにおける各ステップは、記憶部２２に記憶された渦励振流体力データ計算プログラムを処理部２０によって実行することによって実現される。

ステップＳ１０では、処理対象となる線状構造物の構造物条件、当該処理対象への加振条件、及びレイノルズ数を取得する。例えば、処理部２０は、入力部２４からこれらの情報の入力を受けて記憶部２２に記憶させる。処理対象の構造物条件は、処理対象の管の直径（Ｄ）とストレーキの高さ（Ｈ）との比（Ｈ／Ｄ）、処理対象の管の直径（Ｄ）とストレーキの螺旋ピッチ（Ｐ）との比（Ｐ／Ｄ）を含む。また、処理対象への構造物条件は、処理対象への流体の流入角度（χ）を含んでもよい。処理対象への加振条件は、処理対象への無次元加振周波数ｆＤ／Ｖ（ｆ：加振周波数、Ｖ：流速）及び無次元振幅Ａ／Ｄ（Ａ：構造物振動振幅）を含む。

ステップＳ１２では、実験流体力データを取得する。処理部２０は、入力部２４等から渦励振流体力実験方法で得られた実験流体力データを取得して記憶部２２に記憶させる。実験流体力データは、渦励振流体力実験で得られた離散型データであり、予め記憶部２２に記憶させておいてもよい。

ステップＳ１４では、処理対象となる線状構造物に対するレイノルズ数及び加振条件に対応する補間流体力データを求める。処理部２０は、離散型の実験流体力データを補間して、ステップＳ１０で取得したレイノルズ数及び加振条件に対応する流体力データを補間流体力データとして求める。処理部２０は、記憶部２２から実験流体力データを読み出し、ステップＳ１０で取得したレイノルズ数及び加振条件に対応する補間抗力係数、補間揚力方向付加質量係数、補間揚力方向減衰係数を算出する。算出された補間流体力データは、記憶部２２に記憶される。

実験流体力データにおける抗力係数、揚力方向付加質量係数、揚力方向減衰係数は、レイノルズ数及び加振条件以外のパラメータにも依存する。したがって、補間流体力データは、これらのパラメータに依存した抗力係数、揚力方向付加質量係数、揚力方向減衰係数のデータ群となる。

例えば、上記第１範囲〜第４範囲について実験流体力データが取得されている場合、（Ｈ／Ｄ，Ｐ／Ｄ）＝（０，０），（Ｈ／Ｄ，Ｐ／Ｄ）＝（０．１５，５），（Ｈ／Ｄ，Ｐ／Ｄ）＝（０．１５，１５），（Ｈ／Ｄ，Ｐ／Ｄ）＝（０．２５，１５）の４つのデータ群についてステップＳ１０で取得したレイノルズ数及び加振条件に対応する補間抗力係数、補間揚力方向付加質量係数、補間揚力方向減衰係数が算出される。また、実験流体力データの範囲外では、最も近い限界値をレイノルズ数及び加振条件に対応する値として求める。

ステップＳ１６では、処理対象となる線状構造物の構造物条件に対応する計算流体力データを求める。処理部２０は、補間流体力データを補間して、ステップＳ１０で取得した処理対象の構造物条件に対応する流体力データを計算流体力データとして求める。処理部２０は、記憶部２２から補間流体力データを読み出し、ステップＳ１０で取得した構造物条件に対応する計算抗力係数、計算揚力方向付加質量係数、計算揚力方向減衰係数を算出する。算出された計算流体力データは、記憶部２２に記憶される。

ステップＳ１６は、図６のサブルーチンに沿って、ステップＳ１６−１〜Ｓ１６−３の各ステップとして実行される。

ステップＳ１６−１では、処理対象の管の直径（Ｄ）とストレーキの高さ（Ｈ）との比（Ｈ／Ｄ）に応じた抗力係数、揚力方向付加質量係数、揚力方向減衰係数を補間処理により算出する。例えば、上記第１範囲〜第４範囲について補間流体力データが求められている場合、（Ｈ／Ｄ，Ｐ／Ｄ）＝（０，０），（Ｈ／Ｄ，Ｐ／Ｄ）＝（０．１５，１５），（Ｈ／Ｄ，Ｐ／Ｄ）＝（０．２５，１５）の３つのデータ群を用いて、ステップＳ１０で取得した処理対象の（Ｈ／Ｄ）に対応する抗力係数、揚力方向付加質量係数、揚力方向減衰係数が算出される。また、補間流体力データの範囲外では、線形外挿して得られる値を処理対象の（Ｈ／Ｄ）に対応する値として求める。なお、補間（内挿）や外挿は線形的な方法に限られるものでなく、非線形的な方法や多項式補間等の各種の方法が選択し得る。

ステップＳ１６−２では、処理対象の管の直径（Ｄ）とストレーキの螺旋ピッチ（Ｐ）との比（Ｐ／Ｄ）に応じた抗力係数、揚力方向付加質量係数、揚力方向減衰係数に対するピッチ影響補正量（ＲＰ）を求める。例えば、上記第１範囲〜第４範囲について補間流体力データが求められている場合、比（Ｈ／Ｄ）＝０．１５を基準としたデータ群からピッチ影響補正量（ＲＰ）を求める。すなわち、（Ｈ／Ｄ，Ｐ／Ｄ）＝（０．１５，５），（Ｈ／Ｄ，Ｐ／Ｄ）＝（０．１５，１５）の２つのデータ群を補間して、ステップＳ１０で取得した処理対象の（Ｐ／Ｄ）に対応する抗力係数、揚力方向付加質量係数、揚力方向減衰係数に対するピッチ影響補正量（ＲＰ）が算出される。なお、処理対象の比（Ｐ／Ｄ）がデータ群の範囲外である場合には線形外挿してピッチ影響補正量（ＲＰ）を求める。

ステップＳ１６−３では、処理対象への流体の流入角度（χ）に応じた抗力係数、揚力方向付加質量係数、揚力方向減衰係数に対する角度影響補正値（ＲＡ）を算出する。例えば、上記第５範囲について補間流体力データが求められている場合、第５範囲のデータを補間して、処理対象の流入角度（χ）に対応する抗力係数、揚力方向付加質量係数、揚力方向減衰係数に対する角度影響補正値（ＲＡ）が算出される。なお、処理対象の流入角度（χ）がデータ群の範囲外である場合には線形外挿して角度影響補正値（ＲＡ）を求める。

ステップＳ１６−４では、ステップＳ１６−１〜Ｓ１６−３で得られた比（Ｈ／Ｄ）に対する抗力係数、揚力方向付加質量係数、揚力方向減衰係数の補間値、ピッチ影響補正量（ＲＰ）及び角度影響補正値（ＲＡ）から処理対象に対する抗力係数、揚力方向付加質量係数、揚力方向減衰係数の計算流体力データを算出する。すなわち、ステップＳ１６−１で得られた抗力係数、揚力方向付加質量係数、揚力方向減衰係数の各々の補間値に、ステップＳ１６−２で得られたピッチ影響補正量（ＲＰ）及びステップＳ１６−３で得られた角度影響補正値（ＲＡ）を乗算することによって計算流体力データ、すなわち計算抗力係数、計算揚力方向付加質量係数、計算揚力方向減衰係数が算出される。

なお、処理対象への流体の流入角度（χ）を考慮に入れない場合には、ステップＳ１６−３の処理を省略し、抗力係数、揚力方向付加質量係数、揚力方向減衰係数の各々の補間値にピッチ影響補正量（ＲＰ）を乗算することによって計算流体力データを算出することができる。

ステップＳ１８では、ステップＳ１６にて算出された計算流体力データが出力される。処理部２０は、出力部２６から計算流体力データを出力してもよいし、続いて線状構造物の設計処理を行う場合には記憶部２２へ出力し、記憶部２２に計算流体力データを記憶させてもよい。

以上のように、処理対象の管に対する抗力係数、揚力方向付加質量係数、揚力方向減衰係数が計算流体力データとして算出される。

ここで、計算流体力データは、比（Ｈ／Ｄ）、比（Ｐ／Ｄ）及び流入角度（χ）の順に影響を評価して求めることが好適である。これは、抗力係数、揚力方向付加質量係数、揚力方向減衰係数を含む計算流体力データは、線状構造物のストレーキの高さ（Ｈ）とストレーキのピッチ（Ｐ）では、ストレーキの高さ（Ｈ）の影響が大きいので、ストレーキの高さ（Ｈ）の影響評価を優先することが好ましい。また、流入角度（χ）の影響評価は、χ＝０の評価を基準として流入角度の影響を補正するため、最後に行うことが好ましい。補間処理によって計算流体力データを求める場合には、計算流体力データへの影響が大きくパラメータの順に補間処理を行うことによって計算流体力データの精度を高めることができる。

なお、後述する線状構造物の設計処理において有限要素法を適用する場合、線状構造物の要素毎に計算流体力データを求めることが好適である。

＜線状構造物の設計処理＞
渦励振流体力データ計算処理で算出された計算流体力データを用いて渦励振を考慮した線状構造物の設計を行う線状構造物の設計処理について説明する。線状構造物の設計処理では、１本の線状構造物の流れの中での振動を計算する。

線状構造物の設計処理は、有限要素法を用いて、細長い線状構造物を複数の要素に分割して各要素に対する流れの影響を演算する。線状構造物の設計処理は、コンピュータによって線状構造物設計処理プログラムを実行することにより実現される。線状構造物設計処理プログラムを実行するためのコンピュータは、図４に示したコンピュータ２００と同様の構成とすることができる。この場合、記憶部２２に線状構造物設計処理プログラム及び計算流体力データを記憶しておき、処理部２０にて線状構造物設計処理プログラムを実行することによって線状構造物の設計処理を実現する。線状構造物設計処理プログラムは、例えば、有限要素法を扱う周波数領域計算プログラムであるＶＩＶＡＮＡ（ＭＡＲＩＮＴＥＣ製）やＳＨＥＡＲ７（ＭＩＴ製）等を用いて実現することができる。なお、線状構造物の設計処理は、有限要素法以外にもランプドマス法、境界要素法、パネル法等の各種の方法が選択し得るが、有限要素法は構造力学分野で広く使われている数値解析手法であるため利用が容易である。

渦励振を考慮した線状構造物の設計処理は、図７のフローチャートに沿って行われる。当該フローチャートにおける各ステップは、記憶部２２に記憶された線状構造物設計処理プログラムを処理部２０によって実行することによって実現される。

ステップＳ２０では、線状構造物設計に必要な条件の設定が行われる。処理部２０は、入力部２４等から設計対象となる線状構造物の両端の位置、構造物の分割数、剛性等の物理的な条件を取得し、記憶部２２に記憶させる。また、渦励振流体力データ計算処理で得られた計算流体力データを取得して記憶部２２に記憶させる。また、各ステップで得られた計算結果を一時的に記憶してもよい。

ステップＳ２２では、設定された条件下における線状構造物の静的な釣り合い状態を求める。処理部２０は、有限要素法によりステップＳ２０において設定された条件を用いて、流れがない状態における設計対象となる線状構造物の各要素及び要素同士の境界条件について矛盾なく釣り合う条件を算出する。

ステップＳ２４では、線状構造物の固有振動数を算出する。処理部２０は、ステップＳ２２において算出された釣り合い条件に基づいて設計対象となる線状構造物の各要素の固有振動数を求める。

ステップＳ２６では、線状構造物が配置される環境の流体の流れを考慮して線状構造物の各要素に働く流体力を算出する。処理部２０は、線状構造物の各要素に対する計算流体力データ（抗力係数、揚力方向付加質量係数、揚力方向減衰係数）を用いて、設定振動周波数を変更しつつ、線状構造物の各要素に働く流体力を算出し、各要素及び要素同士の境界条件について矛盾なく釣り合う条件を求める。このとき、線状構造物の減衰力等は振動の状況で変化することを考慮することが好適である。また、外力である潮流力も場所により変化していることを考慮することが好適である。

ステップＳ２８では、ステップＳ２６にて算出された流体力が出力される。処理部２０は、出力部２６から流体力を出力してもよいし、記憶部２２に流体力を記憶させてもよい。

このようにして、線状構造物の各要素に働く流体力を計算流体力データに基づいて求めると共に、線状構造物の全体についての釣り合いの条件を探索することによって線状構造物に働く流体力・応答・内力を求めることができる。線状構造物に働く流体力等には、例えば、線状構造物に生ずる振動、周波数、振幅の情報が含まれる。なお、ステップＳ２８の出力として線状構造物に働く曲げモーメントやせん断を得ることもできる。

さらに、得られた流体力等の情報に用いて、線状構造物の疲労評価を行うことができる。疲労評価は、例えば、ステップＳ２８で得られた線状構造物に働く曲げモーメントやせん断力の評価を利用し、材料特性や経時的な条件を考慮して行うことができる。

本実施の形態によれば、線状構造物に生ずる渦励振流体力及びそれによる疲労損傷等への影響をより正確に予測することが可能となる。これによって、線状構造物の設計を容易かつ迅速に行うことを可能とする。例えば、海面から海底まで延びる長いライザーに生ずる渦励振流体力を予測し、用途や利用される海域や要求耐用年数に応じて管径やストレーキ等の仕様（ストレーキ高さ、ピッチ）を最適化して渦励振流体力を低減したライザーを設計することができる。

なお、上記ステップＳ２６では、流体の流れを考慮した線状構造物に働く流体力を算出するが、その算出処理において線状構造物の各要素に対する計算流体力データ（抗力係数、揚力方向付加質量係数、揚力方向減衰係数）に対する条件が変わってしまうことがある。そのような場合、新たな条件を用いて渦励振流体力データ計算処理を再度実行し、線状構造物の各要素に対する計算流体力データを再度算出し、そのデータを用いて線状構造物の設計処理を行うことが好適である。このように、線状構造物の設計処理及び渦励振流体力データ計算処理を繰り返し行うことによって、線状構造物に働く流体力をより正確に求めることができる。

本発明は、流体内に配置される線状構造物に生ずる渦励振流体力の実測及び予測を可能とするものである。本発明の適用範囲は、水中のライザーを含めた線状構造物に限定されるものではなく、空中の屋外設置用ポールや電線等の計算流体力データの算出を通じた疲労破壊の予測や抑制に用いることができる。

１０ ストレーキ、１２ 管、１４ 端板、２０ 処理部、２２ 記憶部、２４ 入力部、２６ 出力部、１００ 部分模型、２００ コンピュータ、Ｓ１０（ステップ１）、Ｓ１２（ステップ２）、Ｓ１４（ステップ３）、Ｓ１６（ステップ４）、Ｓ１８（ステップ５）、Ｓ２０（ステップ６）、Ｓ２２（ステップ７）、Ｓ２６（ステップ８）。

Claims (21)

1. 線状構造物である処理対象に働く流体力データを求めるプログラムであって、
   コンピュータによって、
   前記処理対象の構造物条件と、前記処理対象への加振条件と、レイノルズ数を演算条件として取得するステップ１と、
   線状構造物に対する実験において、加振条件とレイノルズ数との組み合わせに対して得られた実験流体力データを取得するステップ２と、
   前記実験流体力データを補間して前記演算条件に含まれるレイノルズ数及び加振条件における補間流体力データを得るステップ３と、
   前記補間流体力データに対し、前記演算条件に含まれる構造物条件に対応する流体力データを計算流体力データとして求めるステップ４と、
   前記計算流体力データを出力するステップ５を行うことを特徴とする線状構造物の渦励振流体力データ計算プログラム。
2. 請求項１に記載の渦励振流体力データ計算プログラムであって、
   前記演算条件に含まれる構造物条件は、前記処理対象に設けられるストレーキの高さ条件とピッチ条件を含み、
   前記ステップ４は、
   前記高さ条件に対応する流体力データを求めるステップ４−１と、
   前記ピッチ条件に対応する流体力データを求めるステップ４−２を有することを特徴とする渦励振流体力データ計算プログラム。
3. 請求項１又は２に記載の渦励振流体力データ計算プログラムであって、
   前記演算条件に含まれる構造物条件は、前記処理対象への流れの流入角度条件を含み、
   前記ステップ４は、前記流入角度条件に対応する流体力データを求めるステップ４−３を有することを特徴とする渦励振流体力データ計算プログラム。
4. 請求項１に記載の渦励振流体力データ計算プログラムであって、
   前記演算条件に含まれる構造物条件は、前記処理対象に設けられるストレーキの高さ条件とピッチ条件、前記処理対象への流れの流入角度条件を含み、
   前記ステップ４は、
   前記高さ条件に対応する流体力データを求めるステップ４−１と、
   前記ピッチ条件に対応する流体力データを求めるステップ４−２と、
   前記流入角度条件に対応する流体力データを求めるステップ４−３を順に行うことを特徴とする渦励振流体力データ計算プログラム。
5. 請求項１〜４のいずれか１項に記載の渦励振流体力データ計算プログラムであって、
   前記実験流体力データは、前記実験で用いられた線状構造物に設けられたストレーキの高さとピッチ及び／又は前記実験で用いられた線状構造物への流れの流入角度の組み合わせ毎の離散型データであることを特徴とする渦励振流体力データ計算プログラム。
6. 請求項５に記載の渦励振流体力データ計算プログラムであって、
   前記離散型データは、複数のレイノルズ数毎に実験で得られた流体力のデータであることを特徴とする渦励振流体力データ計算プログラム。
7. 請求項１〜５のいずれか１項に記載の渦励振流体力データ計算プログラムであって、
   前記加振条件は、前記処理対象に対する加振周波数及び振幅であることを特徴とする渦励振流体力データ計算プログラム。
8. 請求項１〜７のいずれか１項に記載の渦励振流体力データ計算プログラムであって、
   前記実験流体力データは、実験抗力係数、実験揚力方向付加質量係数及び実験揚力方向減衰係数を含み、
   前記ステップ４は、前記計算流体力データとして計算抗力係数、計算揚力方向付加質量係数及び計算揚力方向減衰係数を求めることを特徴とする渦励振流体力データ計算プログラム。
9. 請求項５又は６に記載の渦励振流体力データ計算プログラムであって、
   前記ステップ４−１は、前記離散型データを補間して前記ストレーキの高さ条件に対応する流体力データを求めることを特徴とする渦励振流体力データ計算プログラム。
10. 請求項５，６及び９のいずれか１項に記載の渦励振流体力データ計算プログラムであって、
    前記ステップ４−２は、前記離散型データを補間して前記ストレーキのピッチ条件に対応する流体力データを求めることを特徴とする渦励振流体力データ計算プログラム。
11. 請求項５，６，９及び１０のいずれか１項に記載の渦励振流体力データ計算プログラムであって、
    前記ステップ４−３は、前記離散型データを補間して前記流入角度条件に対応する流体力データを求めることを特徴とする渦励振流体力データ計算プログラム。
12. 請求項５〜１１のいずれか１項に記載の渦励振流体力データ計算プログラムであって、
    前記線状構造物は、円形断面を有するストレーキ付き円柱形状であり、
    前記実験流体力データは、
    第１範囲：（Ｈ／Ｄ，Ｐ／Ｄ）＝（０，０），Ｒｅ＝１．５×１０⁵〜５．０×１０⁵，ｆＤ／Ｖ＝０．０１〜０．３，Ａ／Ｄ＝０．１〜１．１
    第２範囲：（Ｈ／Ｄ，Ｐ／Ｄ）＝（０．１５，５），Ｒｅ＝１．３×１０⁵〜５．２×１０⁵，ｆＤ／Ｖ＝０．０５〜０．３，Ａ／Ｄ＝０．１〜０．９
    第３範囲：（Ｈ／Ｄ，Ｐ／Ｄ）＝（０．１５，１５），Ｒｅ＝１．３×１０⁵〜５．２×１０⁵，ｆＤ／Ｖ＝０．０５〜０．３，Ａ／Ｄ＝０．１〜０．９
    第４範囲：（Ｈ／Ｄ，Ｐ／Ｄ）＝（０．２５，１５），Ｒｅ＝１．５×１０⁵〜５．０×１０⁵，ｆＤ／Ｖ＝０．０５〜０．３，Ａ／Ｄ＝０．１〜０．９
    ただし、
    Ｈ：実験での線状構造物のストレーキの高さ（ｍ）
    Ｄ：実験での線状構造物の直径（ｍ）
    Ｐ：実験での線状構造物のストレーキの螺旋ピッチ（ｍ）
    Ｒｅ：実験時のレイノルズ数
    ｆＤ／Ｖ：実験無次元加振周波数
    ｆ：実験加振周波数（Ｈｚ）
    Ｖ：実験における線状構造物への流れの流速（ｍ／ｓ）
    Ａ／Ｄ：実験無次元振幅
    Ａ：実験での振動の振幅（ｍ）
    として、前記第１範囲から前記第４範囲の少なくとも１つの範囲で取得されることを特徴とする渦励振流体力データ計算プログラム。
13. 請求項５〜１２のいずれか１項に記載の渦励振流体力データ計算プログラムであって、
    前記線状構造物は、円形断面を有するストレーキ付き円柱形状であり、
    前記実験流体力データは、
    第５範囲：χ＝０〜４５ｄｅｇ，Ｒｅ＝１．５×１０⁵〜５．０×１０⁵，ｆＤ／Ｖ＝０．０１〜０．３，Ａ／Ｄ＝０．１〜１．１
    χ：実験における線状構造物への流れの流入角度（線状構造物への流れが直角のときを０とする）
    Ｄ：実験での線状構造物の直径（ｍ）
    Ｒｅ：実験時のレイノルズ数
    ｆＤ／Ｖ：実験無次元加振周波数
    ｆ：実験加振周波数（Ｈｚ）
    Ｖ：実験における線状構造物への流れの流速（ｍ／ｓ）
    Ａ／Ｄ：実験無次元振幅
    Ａ：実験での振動の振幅（ｍ）
    前記第５範囲の範囲で取得されることを特徴とする渦励振流体力データ計算プログラム。
14. コンピュータによって、
    請求項１〜１３のいずれか１項に記載の渦励振流体力データ計算プログラムを用いて得られた前記計算流体力データを含む設定条件を取得するステップ６と、
    前記設定条件に基づいて前記処理対象の静的な釣り合い条件を求めるステップ７と、
    前記釣り合い条件及び前記計算流体力データを用いて設定振動周波数における前記処理対象に働く流体力を計算するステップ８を行い、
    前記設定振動周波数を変化させつつ前記ステップ８を繰り返して、前記処理対象に働く流体力を計算して前記処理対象が前記釣り合い条件に達したかを判定することを特徴とする渦励振を考慮した線状構造物の設計プログラム。
15. 請求項１４に記載の設計プログラムであって、
    前記ステップ８の繰り返し処理において、前記渦励振流体力データ計算プログラムを用いて前記計算流体力データを再計算し、再計算された前記計算流体力データを用いて前記処理対象のバランス点が前記釣り合い条件に達したかを判定することを特徴とする渦励振を考慮した線状構造物の設計プログラム。
16. 請求項１４又は１５に記載の設計プログラムであって、
    前記ステップ８は、前記処理対象が置かれた潮流条件を用いて前記処理対象に働く流体力を計算することを特徴とする渦励振を考慮した線状構造物の設計プログラム。
17. 請求項１４又は１５に記載の設計プログラムであって、
    前記ステップ８は、前記処理対象の固有振動数を算出し、前記処理対象に働く流体力を計算することを特徴とする渦励振を考慮した線状構造物の設計プログラム。
18. 請求項１４〜１７のいずれか１項に記載の設計プログラムであって、
    前記処理対象に働く流体力を用いて前記処理対象の疲労評価を行うステップ９を備えることを特徴とする渦励振を考慮した線状構造物の設計プログラム。
19. 請求項１４〜１８のいずれか１項に記載の設計プログラムであって、
    前記ステップ８における流体力の計算は、有限要素法を用いて行われることを特徴とする渦励振を考慮した線状構造物の設計プログラム。
20. 請求項１〜１３のいずれか１項に記載の渦励振流体力データ計算プログラムに用いる前記実験流体力データを得るための渦励振流体力実験方法であって、線状構造物の部分模型を強制揺動させて、加振条件とレイノルズ数との組み合わせに対する前記実験流体力データを取得することを特徴とする渦励振流体力実験方法。
21. 請求項２０に記載の渦励振流体力実験方法であって、
    前記部分模型は、実物大模型であることを特徴とする渦励振流体力実験方法。

Images