JP2010170358A - 船体動揺の粘性境界層による減衰比の予測方法 - Google Patents

Abstract

【課題】低コストで高速、高精度な計算が行え、パラメータスタディが容易である船体動揺の粘性境界層による減衰比の予測方法を提供する。
【解決手段】船体がｘｙ平面内で回転動揺している場合を対象とし、船体形状を、船体の静的平衡位置における船体幅方向の対称軸ｙと液面ｘとの交点を原点Ｏとする角座標関数Ｒ＝Ｒ_Ｍ（θ）で表示し、船体から見て境界層の内側の接線方向の流速ν_ｔ．ｉｎと、船体から見て境界層の外側の接線方向の流速ν_{ｔ．ｏｕｔ}を極座標表示で表し、境界層内の、船体に対して相対的な接線方向の流速ν_ｔを、粘性境界層が薄い点に着目して境界条件を満たす式で表し、これから境界層内の流速の極座標成分を求め、液体運動に関する支配方程式系の変分表示を極座標表示で求め、変分直接法に基づき、任意定数に関して変分をとり、各変分の係数を零において、これから固有振動数、固有モードおよび減衰比を求める。
【選択図】図１

Description

本発明は、船体動揺の粘性境界層による減衰比の予測方法に関する。

船体の運動特性を把握するために、フルスケールの実モデルを用いて試験を実施することは通常困難である。そのため、従来から小型模型を用いて試験し、その結果から相似則を用いてフルスケールの実モデルの運動特性を予測することが行われている。
しかし、船体動揺の減衰には、粘性境界層による減衰と、形状抵抗（じゃま板等による剥離渦）による減衰の２つが寄与し、この２つの減衰に関して小型模型と実モデル間の相似則が異なっていて、小型模型では粘性境界層による減衰の寄与が無視できないことが知られている。そのため、粘性境界層による減衰の寄与の正確な評価が重要となる。

そこで船体の横揺れに対する摩擦抵抗及び減衰力について、従来から非特許文献１、２が発表されている。非特許文献１は数値的方法、非特許文献２は実験推定式による方法によるものである。
なお本発明に関連する方法は、非特許文献３および特許文献１に開示されている。

特開２００８−２５７６５４号、「大振幅スロッシング挙動予測方法」

池田良穂、姫野洋司、田中紀男、「横揺れ減衰力について（摩擦成分とビルジキールの直圧力成分）」、関西造船協会誌 第１６１号 昭和５１年６月 加藤弘、「船の横揺れに対する摩擦抵抗について」、造船協会論文集 第１０２号 昭和３２年１１月 林毅、村外志夫、「変分法」、応用数学講座第１３巻、コロナ社、昭和６２年発行

上述したように、粘性境界層による減衰の寄与の正確な評価のために、従来から数値的方法［非特許文献１］や実験推定式による方法［非特許文献２］が適用されているが、これらには計算時間、手間、およびコストが多大であり、パラメータスタディが不便であり、理論的根拠に限界があるという問題がある。

また、ＣＦＤ（流体解析ソフト）等では全減衰を一括して解かざるを得ないため、２つの減衰の寄与の分離評価ができず、小型模型の実験結果から実モデルの減衰を求めるためには適切でないという問題がある。

本発明は、上述した問題点を解決するために創案されたものである。すなわち、本発明の目的は、低コストで高速、高精度な計算が行え、パラメータスタディが容易である船体動揺の粘性境界層による減衰比の予測方法を提供することにある。

本発明によれば、船体がｘｙ平面内で回転動揺している場合を対象とし、
船体形状を、船体の静的平衡位置における船体幅方向の対称軸ｙと液面ｘとの交点を原点Ｏとする角座標関数Ｒ＝Ｒ_Ｍ（θ）（Ｍ：船体と液体の界面の静的平衡位置、θ：ｙ軸からの回転角度）で表示し（後述する式（１））、
船体から見て境界層の内側の接線方向の流速ν_ｔ．ｉｎと、船体から見て境界層の外側の接線方向の流速ν_{ｔ．ｏｕｔ}を極座標表示で表し（後述する式（３１）（３２））、
境界層内の、船体に対して相対的な接線方向の流速ν_ｔを、粘性境界層が薄い点に着目して境界条件を満たす式で表し（後述する式（３９）（４０））、これから境界層内の流速の極座標成分（式（４１））を求め、
液体運動に関する支配方程式系の変分表示を極座標表示で求め、変分直接法に基づき、任意定数に関して変分をとり、各変分の係数を零において、これから固有振動数、固有モードおよび減衰比を求める、ことを特徴とする船体動揺の粘性境界層による減衰比の予測方法が提供される。

本発明の好ましい実施形態によれば、ラグランジュアン密度が液圧に等しいことに基づき、所要の支配方程式系の変分表示を極座標表示で求め（後述する式（１６））、
液体運動の連続性を示すラプラス方程式の解φ（Ｒ，θ，ｔ）（ｔ：時間）を極座標表示で求め（後述する式（１７））、
船体の回転変位の自由振動解θ_Ｇ（ｔ）を極座標表示で表し（後述する式（２２））、
前記φ（Ｒ，θ，ｔ）とθ_Ｇ（ｔ）を前記支配方程式系の変分表示に代入して任意定数に関して変分をとり、各変分の係数を零において、任意定数に関する斉次方程式を標準固有値問題の形に導き（後述する式（２３））、これから固有振動数と固有モードを決定する。

また、前記φ（Ｒ，θ，ｔ）とθ_Ｇ（ｔ）において、ｅ^ｉωｔ→モード座標ｑ（ｔ）の置き換えを行い、その結果（後述する式（２４）（２５））を前記支配方程式系の変分表示（式（１６））に代入してモード座標に関して変分をとり、モード座標の時間微分に関する変分に対して部分積分を行って、モード方程式を導き（後述する式（２６））、
前記境界層内の流速の極座標成分を、Ｒ，θ方向のＮａｖｉｅｒ−Ｓｔｏｋｅｓ方程式から生じる境界層内の仮想仕事に代入して、粘性仮想仕事のモード座標を求め、
前記モード方程式に粘性仮想仕事のモード座標（後述する式（４４））を加算して、粘性を考慮したモード方程式の変分形を得（後述する式（４５））、この式から減衰比ζを計算する（後述する式（４６））。

上述したように、本発明の発明者は、極座標の導入によって船体表面の動径座標を角度の関数（後述する式（１））として与えることにより、任意形状の実船に対して、変分直接法に基づく解析的方法で粘性境界層による減衰が評価できる方法を創案した。
また、変分直接法に基づく解析的方法を用いるために、粘性境界層が薄い点に着目して境界層内流れ場の適切な関数近似（後述する式（４１））を行った。
これにより、低コストで高速、高精度な計算が行え、かつパラメータスタディが容易である効果が得られる。

解析モデルの模式図である。 固有振動数と無次元減衰定数との関係図である。 慣性モーメントと重心位置の関係図である。 慣性モーメントと減衰比の関係図である。 式（４６）中の諸定数の変化を示す図である。 形状ｂ／ａに対する諸パラメータの変化を示す図である。 回転中心ｙ_ｃに対する諸パラメータの変化を示す図である。

以下、本発明の好ましい実施形態を添付図面に基づいて詳細に説明する。なお、各図において共通する部分には同一の符号を付し、重複した説明を省略する。

１． 予測方法
１．１ 問題の設定と座標系
図１のように、船体がｘｙ平面内で回転動揺している場合を考える。
ここで、Ｃ：回転中心、Ｇ：重心、θ_Ｇ（ｔ）：回転変位、Ｍ：船体と液体の界面の静的平衡位置である。
船体は剛体とし、ｘｙ平面に直交するｚ方向に変化のない２次元問題を考える。定式化は、液体領域の境界を有限の位置にある剛体面Ｗ（図には示していない）として行う。

本発明による予測方法は、図１のように、液面（Ｌｉｑｕｉｄ ｓｕｒｆａｃｅ）とｙ軸（船体の静的平衡位置における船体幅方向の対称軸）の交点を原点Ｏとする極座標（Ｒ，θ）（０≦θ≦θ_ｍａｘ＝π／２）を設定することによって、ＭとＷの動径座標を角座標の関数として数１の式（１）（２）の形に表す。
このような極座標の設定により、θが液面上で一定（θ＝±９０°）となるため、液体運動の解の基底となる直交関数系が解析的に決定でき、任意の船体形状に応じて式（１）の関数Ｒ_Ｍ（θ）を与えることにより直交関数展開法による低コストの高速計算が可能となる。図１には一例として、船体形状が式（３）によって与えられる楕円形状の場合を示している。
船体の動揺減衰の模型実験では、固定軸回りに動揺させ粘性摩擦減衰力を計測する。本発明の解析モデルは、このような模型実験のシミュレートを目的とする。

１．２ 変分原理による支配方程式の導出
直交関数展開法を用いるためには、支配方程式を変分原理に基づき重み付きの形に導く必要がある。液体運動に関する変分原理は、ラグランジュアン密度（単位体積当たりのラグランジュアン）が液圧に等しいことに基づき、数２の式（４）のように表される。ここで、Ｌ_ｆ：液体のラグランジュアン、Ｖ： 液体の占める空間、Ｐ_ｌ：液圧である。
液圧Ｐ_ｌは、渦なし非定常流に関する圧力方程式より、速度ポテンシャルφを用いて式（５）のように表される。ここでφ：速度ポテンシャル、ρ_ｆ：液体の密度、ｇ：重力加速度、・付Ｇ（ｔ）：任意の時間関数である。
なお、この出願において、式中に「・」を上に付した記号は、１回微分を示し、式中に「・・」を上に付した記号は、２回微分を示す。また文中では便宜上、それぞれ「・付」「・・付」と記載する。

液体領域が、液面変位と、船体との界面の法線方向変位を介して可変分となることに注意して、式（４）の変分計算を実行する。
次に、船体の回転変位θ_Ｇによるラグランジュアン時間積分の変分式（６）を加算して、系全体のラグランジュアン時間積分の停留条件を導くと、所要の支配方程式系の変分表示が式（７）のように得られる。ここで、θ_Ｇ：船体の回転変位、Ｉ：回転中心回りの慣性モーメント、ｍ：船体の質量である。

速度ポテンシャル、液面変位、船体の変位、任意時間関数の変分が任意かつ独立であることから、支配方程式と境界条件が次のように導かれる。
すなわち、式（７）の第１項より、液体運動の連続性を表すラプラス方程式（８）が導かれる。第２項より、剛体と仮定された境界面で、その法線方向の流速成分が零となる境界条件（９）が導かれる。第３項より、液体Ｆで、その法線方向の液面振動速度と流体粒子速度が等しい条件が導かれる。第４項より、液面で圧力が零となる境界条件が導かれる。第５項より、船体と液体との界面で、その法線方向の船体と液体の速度が等しい条件が導かれる。第６から８項（船体の変位の変分を含む）より、重力と液圧を受ける船体の運動方程式が導かれる。最後の項（第９項）より、液体の非圧縮性の仮定に基づく体積一定条件が導かれる。体積一定条件は、他の運動学的条件から導けるので、以下では考慮しない。

１．３ 変分原理の極座標表示
解析的方法を導入するため、上に導いた変分原理（７）の極座標表示を導く。
まず、式（７）第８項では、運動界面の法線ベクトルＮｓ方向の変位ｕ_ｎとその変分、静圧、運動界面の面積要素の極座標表示を数３の式（１０）〜（１３）のように導いて代入する。
次に、式（７）第５項では、上と異なり、残差（０となるべき括弧中の式）に静的項がないため、運動界面Ｓを静的位置Ｍで近似し、Ｍの面積要素、法線ベクトルを式（１）より式（１４）（１５）のように導いて代入する。
式（７）第２項についても同様に、式（２）より面積要素と法線ベクトルを導いて代入する。
式（７）第３、４項では、次の処理をする：式（５）の液圧の動的項を静的液面上（θ＝θ_ｍａｘ＝π／２）で評価する。
式（５）の液圧の静圧項を運動液面上ｙ＝ηで評価する。
静的液面で、式（１５ａ）となることを用いる。
以上により、式（７）の極座標表示が式（１６）のように導かれる。

１．４ 自由振動解析
式（８）の解は、数４の式（１７）のように表される。ここでａ_ｎは任意の定数で、固有関数と固有値は式（１８）〜式（２０）によって与えられる。
液面の変位ηは、式（１６）第３項から式（２１）によって定められる。
船体の回転変位の自由振動解を式（２２）のように表す。ここで、θ_Ｇ０は任意定数である。
式（１７）、（２２）を式（１６）に代入して任意定数ａ_ｎ、θ_Ｇ０に関して変分をとり、各変分の係数を零におくと、任意定数に関する斉次方程式系が式（２３）の標準固有値問題の形に導かれる。ここで行列Ｍ、Ｋ：質量行列および剛性行列、行列Ｘ：任意定数ａ_ｎ（ｎ≧１）、θ_Ｇ０を集めた列ベクトルである。
式（２３）から固有振動数と固有モードが決定される。

１．５ 粘性減衰比の計算
このモードに関する粘性減衰比を計算するため、モード座標に関する方程式を導く。式（１７）、（２２）で
ｅ^ｉωｔ→ モード座標ｑ（ｔ）
の置き換えることによって解を数５の式（２４）（２５）のように設定する。
式（２４）、（２５）を式（１６）に代入してモード座標に関して変分をとり、モード座標の時間微分に関する変分に対しては部分積分を行うと、モード方程式が式（２６）の形に導かれる。ここでＭ_０、Ｋ_０：積分によって生じる定数であり、それぞれ質量および剛性パラメータである。
式（２６）に、粘性境界層による減衰項を導入する。船体から見て境界層の内側の接線方向の流速ν_ｔ．ｉｎは、船体の接線方向の速度として次式（２７）によって与えられる。ここで、ｔ_Ｍは接線方向の単位ベクトルであり、式（１５）によって与えられるＮ_Ｍとの直交条件より式（２８）のように求められる。
船体から見て境界層の外側の接線方向の流速ν_{ｔ．ｏｕｔ}は、Ｒ，θ方向の流速成分の接線方向成分を加算することにより、式（２９）のように計算できる。
式（２９）の計算に対してｔ_Ｍは、式（２７）と異なり内積計算の便宜のため座標（Ｒ，θ）を用いて表す必要があり、式（１５）の法線ベクトルとの直交性より、式（３０）のようになる。

式（２７）、（２９）を数６の式（３１）（３２）のように記す。ここで、式（３３）（３４）である。
境界層内の、船体に対して相対的な接線方向の流速ν_ｔは、式（３５）〜（３７）の境界条件を満たす。ここで、ξ：境界層に垂直に船体表面から液体に向ってとった座標、ξ_１：液体と剛体との相対運動がそれらの界面の接線方向に、前節で求めた固有振動数で生じる場合の境界層厚さである。
境界層厚さは、式（３８）で表される。
式（３７）は、粘性応力が零になる主流に応力が連続的に接続する条件である。
境界層は非常に薄いので、式（３５）、（３６）、（３７）を満たすν_ｔを式（３９）のように表す。ここで式（４０）である。

Ｎａｖｉｅｒ−Ｓｔｏｋｅｓ方程式の粘性項から粘性応力のなす仮想仕事を計算するために、まず、流速の極座標成分を数７の式（４１）のように表す。
ここでｆ（Ｒ，θ）は、式（４０）の関数を式（４２）の極座標の関数として書いたものである。
Ｒ，θ方向のＮａｖｉｅｒ−Ｓｔｏｋｅｓ方程式から生じる境界層内の仮想仕事は式（４３）によって計算できる。
式（４１）を式（４３）に代入することによって、粘性仮想仕事は式（４４）の形にモード座標で表せる。ここで、Ｃ_０は積分によって生じる定数である。
式（４４）を式（２６）に加算することによって粘性を考慮したモード方程式の変分形が式（４５）のように得られる。
式（４５）より減衰比ζは式（４６）によって計算できる。
ここで、式（４７）のωは前節で求めた固有振動数に等しい。

２． 数値例題
本発明では、船体の形状を式（３）によって表される楕円によって与え、まず、数８の式（４８）に示すパラメータについて数値計算を行った。
最初に、過去の文献の結果との比較による妥当性検証のため、固有振動数（Ｅｉｇｅｎｆｒｅｑｕｅｎｃｙ）と式（４９）の無次元減衰定数（Ｄｉｍｅｎｓｉｏｎｌｅｓｓ ｄａｍｐｉｎｇ ｃｏｎｓｔａｎｔ）との関係を計算した。その結果を図２に示す。
ここで、Δは排水容積、Ｂは船体の幅２ｂであり、Δ＝３５．３４ｍ^３、Ｂ＝０．３ｍである。比較のため、［非特許文献１］の図６に与えられた、［非特許文献２］の方法に基づく値を丸印で示した。この値の妥当性は［非特許文献１］で実施された実験でも確認されている。この実験での重りをのせることによる固有振動数の変化は、図３のような慣性モーメントと重心位置の関係と、ｍ＝Ｉ／（０．１）^２（ｋｇ）によって模擬している。
図２より、本論文の解析と［非特許文献１、２］による計算、実験との一致は、良好であることが確認できる。

図４に、減衰比（Ｄａｍｐｉｎｇ ｒａｔｉｏ）の慣性モーメントＩに対する依存性を示す。図４より、減衰比は慣性モーメントＩが増大すると減少することが分かる。

図５（ａ）、（ｂ）、（ｃ）に、物理的考察のため、減衰比を与える式（４６）に現れる諸定数（減衰定数Ｃ_０、質量パラメータＭ_０、固有振動数ω／２π）の慣性モーメント依存性を示す。まず、図５（ｃ）のように、慣性モーメントＩが増加すると固有振動数ω／２πが高くなることが分かる。これは、図３によって与えられる重心低下に伴う式（７）の重力復元力項の増加を意味する。固有振動数ω／２πが高くなると、式（３８）に基づき境界層が薄くなる。減衰定数Ｃ_０を与える式（４３）において支配的となる∂^２ｆ／∂Ｒ^２は式（４２）に示すように（境界層厚さ）^−２に比例し、Ｒに関する積分範囲は（境界層厚さ）に比例する。したがって減衰定数は（境界層厚さ）^−１に比例し、固有振動数が高くなると増加する。このように、図５（ａ）に示された、慣性モーメントＩの増大による減衰定数Ｃ_０の増加は、固有振動数ω／２πの増加に起因する。このように、減衰定数Ｃ_０は、慣性モーメントＩが増加すると増加する。しかし、式（４６）の分母にある質量パラメータＭ_０と固有振動数ω／２πが、図５（ｂ）、（ｃ）に示すように慣性モーメントＩの増加と共に大きくなるため、減衰比は、図４に示すように慣性モーメントＩが増加すると減少する。

図６（ａ）〜（ｄ）に、形状ｂ／ａの変化に対する、減衰比と式（４６）右辺の諸定数（減衰定数Ｃ_０、質量パラメータＭ_０、固有振動数ω／２π）の変化を示す。慣性モーメントはＩ＝１ｋｇｍ^２とし、ｙ_Ｇ＝−０．０２（ｍ）の場合について計算している。図６（ａ）より、ｂ／ａが１から遠ざかると減衰比が増加することが分かる。図６（ａ）、（ｂ）を比較すると、ｂ／ａが１から増加していく場合には、減衰比が減衰定数Ｃ_０よりかなり顕著に増加していくことが分かる。これは、図６（ｄ）に示すように、ｂ／ａが１を超えると、固有振動数ω／２πのｂ／ａの増加に対する減少率が大きくなり、式（４６）の分母が減少するためである。一方、ｂ／ａが１から減少していく場合には、図６（ｃ）、（ｄ）に示すように質量パラメータＭ_０、固有振動数ω／２πともに増加するので、式（４６）の分母は増加する。このため、減衰比の増加は、減衰定数Ｃ_０の増加に比べて緩和される。

図７（ａ）〜（ｄ）に、回転中心位置ｙ_ｃの零からの下降に対する、減衰比と式（４６）右辺の諸定数（減衰定数Ｃ_０、質量パラメータＭ_０、固有振動数ω／２π）の変化を示す。慣性モーメントはＩ＝１ｋｇｍ^２とし、ｙ_Ｇはｙ_Ｇ＝ｙ_ｃ−０．０２（ｍ）として変化させている。ｙ_ｃが下降すると、図７（ｃ）、（ｄ）のように、質量パラメータＭ_０と固有振動数ω／２πが増加して式（４６）の分母が増加する。このため、図７（ａ）に示されている、ｙ_ｃの下降による減衰比の増加は、図７（ｂ）に示されている減衰定数Ｃ_０の増加に比べて、わずかに緩和される。

上述したように、本発明の発明者は、極座標の導入によって船体表面の動径座標を角度の関数（式（１））として与えることにより、任意形状の実船に対して、変分直接法に基づく解析的方法で粘性境界層による減衰が評価できる方法を創案した。
また、変分直接法に基づく解析的方法を用いるために、粘性境界層が薄い点に着目して境界層内流れ場の適切な関数近似（式（４１））を行った。
これにより、低コストで高速、高精度な計算が行え、かつパラメータスタディが容易である効果が得られる。

すなわち、実施例１に示すように、船体の形状を任意に変化させて、船体の固有振動数、減衰比、減衰定数、質量パラメータ等のパラメータスタディを行い、船体の形状を適正化することができる。
また、この計算を、従来のように、実験に基づく推定式や、多大な計算時間を必要とするＣＦＤ等を必要とせずに、純粋に解析的手法のみで行うことができ、計算時間、手間、およびコストを大幅に低減することができる。

なお、本発明は上述した実施形態に限定されず、特許請求の範囲の記載によって示され、さらに特許請求の範囲の記載と均等の意味および範囲内でのすべての変更を含むものである。

ａ 船体形状を現す楕円の半幅、
ｂ 船体形状を現す楕円の高さの半分、
Ｃ 回転中心、
Ｃ_０ モード方程式の減衰定数、
ｅ_Ｒ，ｅ_θ Ｒ、θ方向の単位ベクトル、
Ｆ 液面、
Ｇ 重心、
・付Ｇ（ｔ）：任意の時間関数、
ｇ 重力加速度、
Ｉ 回転中心のまわりの慣性モーメント、
行列Ｋ 剛性行列、
Ｋ_０ モード方程式の剛性パラメータ
Ｌ_ｆ 液体のラグランジュアン、
Ｍ 船体と液体の界面の静的平衡位置、
Ｍ_０ モード方程式の質量パラメータ、
行列Ｍ 質量行列、
ｍ 船体の質量、
Ｎ_Ｆ Ｆの単位法線ベクトル（液体から見て外向き）、
Ｎ_Ｍ Ｍの単位法線ベクトル（液体から見て内向き）、
Ｎ_Ｓ Ｓの単位法線ベクトル（液体から見て内向き）、
Ｎ_Ｗ Ｗの単位法線ベクトル（液体から見て外向き）、
Ｐ_ｌ 液圧、
Ｒ_Ｍ（θ） ＭのＲ座標、
Ｒ_Ｍθ ｄＲ_Ｍ（θ）／ｄθ、
Ｒ_Ｗ（θ） ＷのＲ座標、
Ｒ_Ｗθ ｄＲ_Ｗ（θ）／ｄθ、
Ｓ 動的変位後の船体と液体との界面、
ｕ_ｎ 船体のＮｓ方向の変位、
ｕ_Ｒ，ｕ_θ Ｒ，θ方向の流体変位成分、
ν_Ｒ，ν_θ Ｒ，θ方向の流体変位成分、
Ｖ 液体の占める空間、
Ｗ 液体領域の，液面と，船体と液体との界面を除く境界、
行列Ｘ 任意定数ａ_ｎ（ｎ≧１）、θ_Ｇ０を集めた列ベクトル
θ_Ｇ（ｔ） 船体の回転変位、
θ_Ｇ０ 任意定数、
φ 速度ポテンシャル、
ξ 船体表面から境界層に垂直に計った座標、
ξ_１ 境界層厚さ、
η 液面変位、
μ 液体の粘性係数、
ρ_ｆ 液体の密度、
ω 固有振動数、
ζ 減衰比、

Claims (3)

1. 船体がｘｙ平面内で回転動揺している場合を対象とし、
   船体形状を、船体の静的平衡位置における船体幅方向の対称軸ｙと液面ｘとの交点を原点Ｏとする角座標関数Ｒ＝Ｒ_Ｍ（θ）（Ｍ：船体と液体の界面の静的平衡位置、θ：ｙ軸からの回転角度）で表示し、
   船体から見て境界層の内側の接線方向の流速ν_ｔ．ｉｎと、船体から見て境界層の外側の接線方向の流速ν_{ｔ．ｏｕｔ}を極座標表示で表し、
   境界層内の、船体に対して相対的な接線方向の流速ν_ｔを、粘性境界層が薄い点に着目して境界条件を満たす式で表し、これから境界層内の流速の極座標成分を求め、
   液体運動に関する支配方程式系の変分表示を極座標表示で求め、変分直接法に基づき、任意定数に関して変分をとり、各変分の係数を零において、これから固有振動数、固有モードおよび減衰比を求める、ことを特徴とする船体動揺の粘性境界層による減衰比の予測方法。
2. ラグランジュアン密度が液圧に等しいことに基づき、所要の支配方程式系の変分表示を極座標表示で求め、
   液体運動の連続性を示すラプラス方程式の解φ（Ｒ，θ，ｔ）（ｔ：時間）を極座標表示で求め、
   船体の回転変位の自由振動解θ_Ｇ（ｔ）を極座標表示で表し、
   前記φ（Ｒ，θ，ｔ）とθ_Ｇ（ｔ）を前記支配方程式系の変分表示に代入して任意定数に関して変分をとり、各変分の係数を零において、任意定数に関する斉次方程式を標準固有値問題の形に導き、これから固有振動数と固有モードを決定する、ことを特徴とする請求項１に記載の予測方法。
3. 前記φ（Ｒ，θ，ｔ）とθ_Ｇ（ｔ）において、ｅ^ｉωｔ→モード座標ｑ（ｔ）の置き換えを行い、その結果を前記支配方程式系の変分表示に代入してモード座標に関して変分をとり、モード座標の時間微分に関する変分に対して部分積分を行って、モード方程式を導き、
   前記境界層内の流速の極座標成分を、Ｒ，θ方向のＮａｖｉｅｒ−Ｓｔｏｋｅｓ方程式から生じる境界層内の仮想仕事に代入して、粘性仮想仕事のモード座標を求め、
   前記モード方程式に粘性仮想仕事のモード座標を加算して、粘性を考慮したモード方程式の変分形を得、この式から減衰比ζを計算する、ことを特徴とする請求項１に記載の予測方法。

Images