JP4838829B2 - プロペラによる船体表面の変動圧力推定装置及び方法並びにプログラム - Google Patents
プロペラによる船体表面の変動圧力推定装置及び方法並びにプログラム Download PDFInfo
- Publication number
- JP4838829B2 JP4838829B2 JP2008198257A JP2008198257A JP4838829B2 JP 4838829 B2 JP4838829 B2 JP 4838829B2 JP 2008198257 A JP2008198257 A JP 2008198257A JP 2008198257 A JP2008198257 A JP 2008198257A JP 4838829 B2 JP4838829 B2 JP 4838829B2
- Authority
- JP
- Japan
- Prior art keywords
- propeller
- analysis
- hull
- shape
- cavitation
- 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
Links
- 238000000034 method Methods 0.000 title claims description 43
- 238000004458 analytical method Methods 0.000 claims description 142
- 238000004088 simulation Methods 0.000 claims description 43
- 238000004364 calculation method Methods 0.000 claims description 41
- 230000000694 effects Effects 0.000 claims description 10
- 239000012530 fluid Substances 0.000 claims description 10
- 230000002123 temporal effect Effects 0.000 claims description 4
- 238000009434 installation Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000004904 shortening Methods 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
本発明は、船体に設けられたプロペラに発生するキャビテーションにより、船体表面に誘起される変動圧力を推定するプロペラによる船体表面の変動圧力推定装置及び方法並びにプログラムに関するものである。
船舶のプロペラに発生するキャビテーションにより船体表面に誘起される変動圧力を推定する方法として、例えば、数値計算を用いた流体シミュレーションによる解法が知られている。流体シミュレーションによる解法としては、粘性を考慮しないポテンシャル計算によるものや、粘性影響も考慮するＲＡＮＳ（Raynolds Averaged Navier-Stokes）を基礎においた流体解析技術（ＣＦＤ：Computational Fluids Dynamics）によるもの等が知られている。
特開平１０−１１４２０号公報
特開平９−３２９１０９号公報
しかしながら、上記ポテンシャル計算による解法は、比較的短時間で計算できるが、プロペラに発生するキャビテーションの計算精度が悪く、全体として十分な精度が得られないという問題があった。
また、熱流体解析技術による解法では、プロペラに発生するキャビテーションの計算精度はポテンシャル計算を用いる場合よりも高いが、計算に多大の時間を要するという問題があった。
また、熱流体解析技術による解法では、プロペラに発生するキャビテーションの計算精度はポテンシャル計算を用いる場合よりも高いが、計算に多大の時間を要するという問題があった。
本発明は、上記問題を解決するためになされたもので、所定の精度を確保しながら、処理の低減及び時間短縮を図ることのできるプロペラによる船体表面の変動圧力推定装置及び方法並びにプログラムを提供することを目的とする。
上記課題を解決するために、本発明は以下の手段を採用する。
本発明は、船体に取り付けられたプロペラの形状を入力情報として取得し、取得した該プロペラを包含する解析空間を設定する解析空間設定手段と、境界条件としてプロペラへの流入速度を与え、粘性影響を考慮する第１のシミュレーション計算モデルを用いて、該解析空間における該プロペラのキャビテーション形状を算出する第１の解析手段と、該船体の形状、該プロペラの形状、及び該第１の解析手段によって解析された該キャビテーション形状を入力情報として取得し、これらの情報を用いて、粘性影響を考慮しない第２のシミュレーション計算を行うことで、該プロペラの取り付け位置周辺における該船体表面の変動圧力を算出する第２の解析手段とを具備するプロペラによる船体表面の変動圧力推定装置を提供する。
本発明は、船体に取り付けられたプロペラの形状を入力情報として取得し、取得した該プロペラを包含する解析空間を設定する解析空間設定手段と、境界条件としてプロペラへの流入速度を与え、粘性影響を考慮する第１のシミュレーション計算モデルを用いて、該解析空間における該プロペラのキャビテーション形状を算出する第１の解析手段と、該船体の形状、該プロペラの形状、及び該第１の解析手段によって解析された該キャビテーション形状を入力情報として取得し、これらの情報を用いて、粘性影響を考慮しない第２のシミュレーション計算を行うことで、該プロペラの取り付け位置周辺における該船体表面の変動圧力を算出する第２の解析手段とを具備するプロペラによる船体表面の変動圧力推定装置を提供する。
上記構成によれば、キャビテーションの計算を粘性影響を考慮する第１のシミュレーション計算モデルを用いて行い、このシミュレーション結果であるキャビテーション形状を用いて粘性影響を考慮しない第２のシミュレーション計算を行うことで、船体表面の変動圧力を算出する。
ここで、粘性影響を考慮する第１のシミュレーション計算モデルは、粘性影響を考慮しない第２のシミュレーション計算モデルに比べて精度が高く、一方、第２のシミュレーション計算モデルは、第１のシミュレーション計算モデルよりも処理が少なく、処理時間の短縮を図ることができる。このように、本発明では、船体表面の変動圧力を算出するまでの全ての工程を同じシミュレーション方法を用いて計算するのではなく、船体表面の変動圧力を算出するまでの工程を大きく２つに分け、一方の工程に精度の高い第１のシミュレーションモデルを採用するとともに、他方の工程に時間短縮を図ることのできる第２のシミュレーションモデルを採用したので、両シミュレーション技術の特徴を効果的に用いて船体表面の変動圧力を算出することが可能となる。また、精度の高い第１のシミュレーション計算モデルを、船体表面の変動圧力を算出する工程ではなく、敢えて、キャビテーションの算出に用いることにより、計算精度を効果的に向上させることができる。
ここで、粘性影響を考慮する第１のシミュレーション計算モデルは、粘性影響を考慮しない第２のシミュレーション計算モデルに比べて精度が高く、一方、第２のシミュレーション計算モデルは、第１のシミュレーション計算モデルよりも処理が少なく、処理時間の短縮を図ることができる。このように、本発明では、船体表面の変動圧力を算出するまでの全ての工程を同じシミュレーション方法を用いて計算するのではなく、船体表面の変動圧力を算出するまでの工程を大きく２つに分け、一方の工程に精度の高い第１のシミュレーションモデルを採用するとともに、他方の工程に時間短縮を図ることのできる第２のシミュレーションモデルを採用したので、両シミュレーション技術の特徴を効果的に用いて船体表面の変動圧力を算出することが可能となる。また、精度の高い第１のシミュレーション計算モデルを、船体表面の変動圧力を算出する工程ではなく、敢えて、キャビテーションの算出に用いることにより、計算精度を効果的に向上させることができる。
上記プロペラによる船体表面の変動圧力推定装置において、前記第１の解析手段は、数値流体力学のシミュレーション計算モデルを用いて、該プロペラのキャビテーション形状を算出することとしてもよい。
上記プロペラによる船体表面の変動圧力推定装置において、前記第２の解析手段は、境界要素法を用いた音場解析により、該プロペラの取り付け位置周辺における該船体表面の変動圧力を算出することとしてもよい。
上記プロペラによる船体表面の変動圧力推定装置において、前記第１の解析手段は、前記キャビテーション形状を推定した後に、プロペラの表面及び翼断面に複数の解析点を設定し、各解析点におけるキャビテーションの厚さを算出し、各該解析点におけるキャビテーションの厚さ情報を前記第２の解析手段に出力することとしてもよい。
キャビテーションの厚さ情報を第２の解析手段に出力することで、第２の解析手段ではキャビテーションの厚さ情報をそのまま用いてシミュレーション計算を行うことができる。これにより、第２の解析手段における処理負担を軽減させることができる。
上記プロペラによる船体表面の変動圧力推定装置において、前記第１の解析手段は、前記キャビテーション形状を推定した後に、該プロペラを包含する円筒面を設定し、該円筒面における速度の時間変動を算出し、算出した該円筒面における速度の時間変動を前記第２の解析手段に出力することとしてもよい。
円筒面における速度の時間変動を第２の解析手段に出力することで、第２の解析手段では円筒面における速度の時間変動をそのまま用いてシミュレーション計算を行うことができる。これにより、第２の解析手段における処理負担を軽減させることができる。
上記プロペラによる船体表面の変動圧力推定装置は、前記第２の解析手段によって算出された該船体表面の変動圧力を用いて、船体における振動分析を行う第３の解析手段を備えていてもよい。
このように、第３の解析手段を備えることで、船体における振動分析をも算出することが可能となる。
本発明は、船体に取り付けられたプロペラの形状を入力情報として取得し、取得した該プロペラを包含する解析空間を設定する第１過程と、境界条件としてプロペラへの流入速度を与え、粘性影響を考慮する第１のシミュレーション計算モデルを用いて、該解析空間における該プロペラのキャビテーション形状を算出する第２過程と、該船体の形状、該プロペラの形状、及び該第２過程において解析された該キャビテーション形状を入力情報として取得し、これらの情報を用いて、粘性影響を考慮しない第２のシミュレーション計算を行うことで、該プロペラの取り付け位置周辺における該船体表面の変動圧力を算出する第３過程とを有するプロペラによる船体表面の変動圧力推定方法を提供する。
本発明は、船体に取り付けられたプロペラの形状を入力情報として取得し、取得した該プロペラを包含する解析空間を設定する第１処理と、粘性影響を考慮する第１のシミュレーション計算モデルを用いて、該解析空間における該プロペラのキャビテーション形状を算出する第２処理と、プロペラへの流入速度場等の境界条件、該船体の形状、該プロペラの形状、及び該第２処理において解析された該キャビテーション形状を入力情報として取得し、これらの情報を用いて、粘性影響を考慮しない第２のシミュレーション計算を行うことで、該プロペラの取り付け位置周辺における該船体表面の変動圧力を算出する第３処理とをコンピュータに実行させるためのプロペラによる船体表面の変動圧力推定プログラムを提供する。
本発明によれば、所定の精度を確保しながら、処理の低減及び時間短縮を図ることができるという効果を奏する。
以下に、本発明に係るプロペラによる船体表面の変動圧力推定装置及び方法並びにプログラムの一実施形態について、図面を参照して説明する。
図１は、本発明の一実施形態に係るプロペラによる船体表面の変動圧力推定装置（以下「変動圧力推定装置」という。）の概略構成を示したブロック図である。図１に示すように、本実施形態に係る変動圧力推定装置１０は、コンピュータシステム（計算機システム）であり、ＣＰＵ（中央演算処理装置）１１、ＲＡＭ（Random Access Memory）等の主記憶装置１２、ＨＤＤ（Hard Disk Drive）等の補助記憶装置１３、キーボードやマウスなどの入力装置１４、及びモニタやプリンタなどの出力装置１５、外部の機器と通信を行うことにより情報の授受を行う通信装置１６などで構成されている。
補助記憶装置１３には、各種プログラム（例えば、プロペラによる船体表面の変動圧力推定プログラム等）が格納されており、ＣＰＵ１１が補助記憶装置１３から主記憶装置１２にプログラムを読み出し、実行することにより種々の処理を実現させる。
図１は、本発明の一実施形態に係るプロペラによる船体表面の変動圧力推定装置（以下「変動圧力推定装置」という。）の概略構成を示したブロック図である。図１に示すように、本実施形態に係る変動圧力推定装置１０は、コンピュータシステム（計算機システム）であり、ＣＰＵ（中央演算処理装置）１１、ＲＡＭ（Random Access Memory）等の主記憶装置１２、ＨＤＤ（Hard Disk Drive）等の補助記憶装置１３、キーボードやマウスなどの入力装置１４、及びモニタやプリンタなどの出力装置１５、外部の機器と通信を行うことにより情報の授受を行う通信装置１６などで構成されている。
補助記憶装置１３には、各種プログラム（例えば、プロペラによる船体表面の変動圧力推定プログラム等）が格納されており、ＣＰＵ１１が補助記憶装置１３から主記憶装置１２にプログラムを読み出し、実行することにより種々の処理を実現させる。
図２は、変動圧力推定装置１０が備える機能を展開して示した機能ブロック図である。図２に示されるように、変動圧力推定装置１０は、解析空間設定部（解析空間設定手段）２１と、第１の解析部（第１の解析手段）２２と、第２の解析部（第２の解析手段）２３とを備えている。
解析空間設定部２１には、船体に取り付けられたプロペラの形状、船体におけるプロペラの取り付け位置等、シミュレーションを行うのに必要な情報が入力されるようになっている。解析空間設定部２１は、船体に取り付けられたプロペラの形状を入力情報として取得すると、図３に示すように、取得した該プロペラ３０の表面を分割して、プロペラ３０の表面上に複数のグリッドを作成する。更に、解析空間設定部２１は、図４に示すように、該プロペラを包含する解析空間を設定し、この解析空間を分割することにより、複数のグリッドを作成する。このとき、プロペラ近傍のグリッドは細かく、プロペラから離れたところのグリッドは粗く設定するとよい。このようにすることで、グリッド数を抑えることができる。このときのグリッド数は、例えば、１００万個のオーダーである。
解析空間設定部２１には、船体に取り付けられたプロペラの形状、船体におけるプロペラの取り付け位置等、シミュレーションを行うのに必要な情報が入力されるようになっている。解析空間設定部２１は、船体に取り付けられたプロペラの形状を入力情報として取得すると、図３に示すように、取得した該プロペラ３０の表面を分割して、プロペラ３０の表面上に複数のグリッドを作成する。更に、解析空間設定部２１は、図４に示すように、該プロペラを包含する解析空間を設定し、この解析空間を分割することにより、複数のグリッドを作成する。このとき、プロペラ近傍のグリッドは細かく、プロペラから離れたところのグリッドは粗く設定するとよい。このようにすることで、グリッド数を抑えることができる。このときのグリッド数は、例えば、１００万個のオーダーである。
第１の解析部２２は、粘性影響を考慮する第１のシミュレーション計算モデルを用いて、解析空間設定部２１によって設定された該解析空間におけるプロペラ３０のキャビテーション形状を算出する。例えば、第１の解析部２２は、数値流体力学（ＣＦＤ）のシミュレーション計算モデルを用いて、該プロペラのキャビテーション形状を算出する。このように、粘性影響を考慮するＣＦＤ等のシミュレーション計算モデルを用いることにより、精度を高いシミュレーション結果を得ることが可能となる。
具体的には、第１の解析部２２は、解析空間設定部２１によって設定されたグリッドを用いて、以下に示す式を解くことで、該解析空間におけるプロペラ３０の動的なキャビテーション形状を算出する。本実施形態では、プロペラ周囲の粘性流場解析の基礎式として、以下の（１）式に示す非定常非圧縮性流体ＲＡＮＳ（Reynolds-Averaged Navier-Stokes）方程式と、（２）式に示す連続の条件式とを用いる。なお、以下の（１）式、（２）式は、デカルト座標テンソル表記にて示されている。また、定常解析の場合は、非定常項が省略される。
乱流に対しては、ｕi、ｐ、τij等にアンサンブル平均値（時間平均値）が当てはめられ、そのとき、τijは以下の（３）〜（４）式のように表される。
第１の解析部２２は、上記（１）式及び（２）式を解くことにより、時々刻々と変化するプロペラ周辺に発生するキャビテーション形状を時間毎に推定する。
図５及び図６は、第１の解析部２２による解析結果の一例を示している。図５において、範囲Ｒは、キャビテーションが発生する範囲を示している。また、図６は、図５に示すＡ−Ａの半径位置でのキャビテーションの発生状況を示した図である。
図５及び図６は、第１の解析部２２による解析結果の一例を示している。図５において、範囲Ｒは、キャビテーションが発生する範囲を示している。また、図６は、図５に示すＡ−Ａの半径位置でのキャビテーションの発生状況を示した図である。
続いて、第１の解析部２２は、第２の解析部２３がシミュレーションを行いやすいように、上記解析結果から任意の解析点におけるキャビテーションの厚み情報を算出する。以下、具体的に説明する。
まず、第１の解析部２２は、図６に示されるプロペラの翼断面及び翼表面に複数の解析点（座標）を設定し、この解析点における法線ベクトル及びキャビテーションの厚みを求める。解析点は、例えば、プロペラの半径方向に１０点、翼弦方向に１０点、回転方向に３６点（１０度間隔）設定される。図７は、ある解析点Ａにおける法線ベクトル及びキャビテーションの厚みの一例を示した図である。
第１の解析部２２は、各解析点におけるキャビテーションの厚みを求めると、このキャビテーションの厚み情報とキャビテーションの推定結果とを第２の解析部２３に出力する。
まず、第１の解析部２２は、図６に示されるプロペラの翼断面及び翼表面に複数の解析点（座標）を設定し、この解析点における法線ベクトル及びキャビテーションの厚みを求める。解析点は、例えば、プロペラの半径方向に１０点、翼弦方向に１０点、回転方向に３６点（１０度間隔）設定される。図７は、ある解析点Ａにおける法線ベクトル及びキャビテーションの厚みの一例を示した図である。
第１の解析部２２は、各解析点におけるキャビテーションの厚みを求めると、このキャビテーションの厚み情報とキャビテーションの推定結果とを第２の解析部２３に出力する。
第２の解析部２３は、各解析点におけるキャビテーションの厚み情報を第１の解析部２２から受け取ると、このキャビテーションの厚み情報、プロペラの形状、該プロペラが取り付けられている船体の形状、及び該船体におけるプロペラの取り付け位置等を入力情報として用いて、粘性影響を考慮しない第２のシミュレーション計算を行うことで、該プロペラの取り付け位置周辺における該船体表面の変動圧力を算出する。
具体的には、第２の解析部２３は、図８に示すように、プロペラの取り付け位置に近い船体部分、例えば、船体後尾部分の表面に複数のグリッドを設定する。このとき、プロペラ近傍のグリッドは細かく、プロペラから離れたところのグリッドは粗く設定するとよい。このようにすることで、グリッド数を抑えることができ、処理の軽減及び処理時間の短縮を図ることができる。
続いて、第２の解析部２３は、上記キャビテーションの厚みを音源としてみなし、境界要素法（ＢＥＭ：Boundary Element Method）を用いた音場解析を行うことで、該船体表面の変動圧力を算出する。
ここで、シミュレーションに必要となる上記各種条件は、入力装置１４（図１参照）から適宜入力されるものとする。
第２の解析部２３は、各時間におけるプロペラによる船体表面の変動圧力を予測すると、その変動圧力分布を解析結果として出力する。図９は、第２の解析部２３の解析結果である変動圧力分布の一例を示した図である。
続いて、第２の解析部２３は、上記キャビテーションの厚みを音源としてみなし、境界要素法（ＢＥＭ：Boundary Element Method）を用いた音場解析を行うことで、該船体表面の変動圧力を算出する。
ここで、シミュレーションに必要となる上記各種条件は、入力装置１４（図１参照）から適宜入力されるものとする。
第２の解析部２３は、各時間におけるプロペラによる船体表面の変動圧力を予測すると、その変動圧力分布を解析結果として出力する。図９は、第２の解析部２３の解析結果である変動圧力分布の一例を示した図である。
次に、上述した変動圧力推定装置１０の作用について図８を参照して説明する。なお、図１０に示された各処理については、ＣＰＵ１１が補助記憶装置１３に記憶されている各種プログラムを主記憶装置１２に読み出して実行することにより実現されるものである。
まず、プロペラ形状等の入力情報が解析空間設定部２１に入力されると、解析空間設定部２１によりプロペラの表面に複数のグリッドが設定されるとともに、該プロペラを包含する解析空間が設定され、更に、この解析空間が分割されて複数のグリッドが設定される（ステップＳＡ１）。続いて、該解析空間におけるプロペラのキャビテーション形状を予測するための数値流体力学を用いたシミュレーション解析が第１の解析部２２により行われる。続いて、第１の解析部２２により、推定したプロペラのキャビテーション形状に各解析点が設定され、各解析点におけるキャビテーションの厚みが算出される。このキャビテーションの厚み情報は、プロペラのキャビテーション形状とともに、第２の解析部２３に出力される（ステップＳＡ２）。
第２の解析部２３においては、第１の解析部２２から入力されたキャビテーションの厚み情報等を用いて第２のシミュレーション計算が行われることにより、プロペラの取り付け位置周辺における該船体表面の変動圧力が算出される（ステップＳＡ３）。
第２の解析部２３においては、第１の解析部２２から入力されたキャビテーションの厚み情報等を用いて第２のシミュレーション計算が行われることにより、プロペラの取り付け位置周辺における該船体表面の変動圧力が算出される（ステップＳＡ３）。
以上、説明してきたように、本実施形態に係る変動圧力推定装置及び方法並びにプログラムによれば、キャビテーションの計算をＣＦＤ等のシミュレーションモデルで行い、このシミュレーション結果であるキャビテーション形状を用いてＢＥＭによる音場解析を行うことにより、船体表面の変動圧力を予測する。
ここで、ＣＦＤは音場解析等のポテンシャル計算に比べて精度が高く、また、音場解析等のポテンシャル計算はＣＦＤに比べて処理時間が短いという長所を有する。このように、本実施形態によれば、船体表面の変動圧力を算出するまでの全ての工程を同じシミュレーション方法を用いて計算するのではなく、船体表面の変動圧力を算出するまでの工程を大きく２つに分け、一方の工程に精度の高いＣＦＤを用いた解析技術を採用し、他方の工程に時間短縮を図ることのできるポテンシャル計算を用いた解析技術を採用したので、両解析技術の特徴を効果的に用いて船体表面の変動圧力を算出することが可能となる。これにより、所定の精度を確保しながら、処理の低減及び時間短縮を図ることができる。
更に、解析空間設定部２１並びに第２の解析部２３が、プロペラ近傍のグリッドを細かく設定し、プロペラから離れた位置におけるグリッドを比較的大きめに設定するので、計算時間や計算機のメモリ容量を更に削減することが可能となる。更に、第２の解析部２３が、船体後尾部分の表面にグリッドを設定するので、計算時間や計算機のメモリの容量を更に削減することができる。
なお、上述した本実施形態においては、第１の解析部２２が、推定したキャビテーションの形状に基づいて各解析点におけるキャビテーションの厚さを算出し、これらの情報を第２の解析部２３に出力していたが、第２の解析部２３に出力する情報は、キャビテーションの厚さに限られない。
例えば、キャビテーション形状を推定した後に、図１１及び図１２に示されるように、プロペラを包含する円筒面を設定し、この円筒面における速度の時間変動を求め、この情報を第２の解析部２３に出力することとしてもよい。このとき、円筒面は、例えば、プロペラよりも１割程度大きいものを設定するとよい。
例えば、キャビテーション形状を推定した後に、図１１及び図１２に示されるように、プロペラを包含する円筒面を設定し、この円筒面における速度の時間変動を求め、この情報を第２の解析部２３に出力することとしてもよい。このとき、円筒面は、例えば、プロペラよりも１割程度大きいものを設定するとよい。
円筒面における速度の時間変動の情報が入力された第２の解析部２３は、円筒面における速度の時間変動を音源とみなして、上記境界要素法を用いた音場解析を行うことで、該船体表面の変動圧力を算出する。
このように、キャビテーションの厚さ情報に代えて、プロペラを包含する円筒面における速度の時間変動を第２の解析部２３に出力することによっても、上述した実施形態と同様の効果を得ることができる。
このように、キャビテーションの厚さ情報に代えて、プロペラを包含する円筒面における速度の時間変動を第２の解析部２３に出力することによっても、上述した実施形態と同様の効果を得ることができる。
また、本実施形態においては、該船体表面の変動圧力の分布を最終的な出力情報として出力することとしていたが、この船体表面の変動圧力に基づいて更に船体における振動分析を行う第３の解析部（第３の解析手段：図示略）を設け、船体における振動分析結果を最終的な出力とすることとしてもよい。
１０ プロペラによる船体表面の変動圧力推定装置
１１ ＣＰＵ
１２ 主記憶装置
１３ 補助記憶装置
１４ 入力装置
１５ 出力装置
１６ 通信装置
２１ 解析空間設定部
２２ 第１の解析部
２３ 第２の解析部
３０ プロペラ
１１ ＣＰＵ
１２ 主記憶装置
１３ 補助記憶装置
１４ 入力装置
１５ 出力装置
１６ 通信装置
２１ 解析空間設定部
２２ 第１の解析部
２３ 第２の解析部
３０ プロペラ
Claims (8)
- 船体に取り付けられたプロペラの形状を入力情報として取得し、取得した該プロペラを包含する解析空間を設定する解析空間設定手段と、
境界条件としてプロペラへの流入速度を与え、粘性影響を考慮する第１のシミュレーション計算モデルを用いて、該解析空間における該プロペラのキャビテーション形状を算出する第１の解析手段と、
該船体の形状、該プロペラの形状、及び該第１の解析手段によって解析された該キャビテーション形状を入力情報として取得し、これらの情報を用いて、粘性影響を考慮しない第２のシミュレーション計算を行うことで、該プロペラの取り付け位置周辺における該船体表面の変動圧力を算出する第２の解析手段と
を具備するプロペラによる船体表面の変動圧力推定装置。 - 前記第１の解析手段は、数値流体力学のシミュレーション計算モデルを用いて、該プロペラのキャビテーション形状を算出する請求項１に記載のプロペラによる船体表面の変動圧力推定装置。
- 前記第２の解析手段は、境界要素法を用いた音場解析により、該プロペラの取り付け位置周辺における該船体表面の変動圧力を算出する請求項１または請求項２に記載のプロペラによる船体表面の変動圧力推定装置。
- 前記第１の解析手段は、前記キャビテーション形状を推定した後に、プロペラの表面及び翼断面に複数の解析点を設定し、各解析点におけるキャビテーションの厚さを算出し、各該解析点におけるキャビテーションの厚さ情報を前記第２の解析手段に出力する請求項１から請求項３のいずれかに記載のプロペラによる船体表面の変動圧力推定装置。
- 前記第１の解析手段は、前記キャビテーション形状を推定した後に、該プロペラを包含する円筒面を設定し、該円筒面における速度の時間変動を算出し、算出した該円筒面における速度の時間変動を前記第２の解析手段に出力する請求項１から請求項３のいずれかに記載のプロペラによる船体表面の変動圧力推定装置。
- 前記第２の解析手段によって算出された該船体表面の変動圧力を用いて、船体における振動分析を行う第３の解析手段を備える請求項１から請求項５のいずれかに記載のプロペラによる船体表面の変動圧力推定装置。
- 船体に取り付けられたプロペラの形状を入力情報として取得し、取得した該プロペラを包含する解析空間を設定する第１過程と、
境界条件としてプロペラへの流入速度を与え、粘性影響を考慮する第１のシミュレーション計算モデルを用いて、該解析空間における該プロペラのキャビテーション形状を算出する第２過程と、
該船体の形状、該プロペラの形状、及び該第２過程において解析された該キャビテーション形状を入力情報として取得し、これらの情報を用いて、粘性影響を考慮しない第２のシミュレーション計算を行うことで、該プロペラの取り付け位置周辺における該船体表面の変動圧力を算出する第３過程と
を有するプロペラによる船体表面の変動圧力推定方法。 - 船体に取り付けられたプロペラの形状を入力情報として取得し、取得した該プロペラを包含する解析空間を設定する第１処理と、
粘性影響を考慮する第１のシミュレーション計算モデルを用いて、該解析空間における該プロペラのキャビテーション形状を算出する第２処理と、
プロペラへの流入速度場等の境界条件、該船体の形状、該プロペラの形状、及び該第２処理において解析された該キャビテーション形状を入力情報として取得し、これらの情報を用いて、粘性影響を考慮しない第２のシミュレーション計算を行うことで、該プロペラの取り付け位置周辺における該船体表面の変動圧力を算出する第３処理と
をコンピュータに実行させるためのプロペラによる船体表面の変動圧力推定プログラム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2008198257A JP4838829B2 (ja) | 2008-07-31 | 2008-07-31 | プロペラによる船体表面の変動圧力推定装置及び方法並びにプログラム |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2008198257A JP4838829B2 (ja) | 2008-07-31 | 2008-07-31 | プロペラによる船体表面の変動圧力推定装置及び方法並びにプログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2010036600A JP2010036600A (ja) | 2010-02-18 |
| JP4838829B2 true JP4838829B2 (ja) | 2011-12-14 |
Family
ID=42009618
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2008198257A Expired - Fee Related JP4838829B2 (ja) | 2008-07-31 | 2008-07-31 | プロペラによる船体表面の変動圧力推定装置及び方法並びにプログラム |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4838829B2 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010038556A (ja) * | 2008-07-31 | 2010-02-18 | Mitsubishi Heavy Ind Ltd | プロペラの圧力変動推定装置及び方法並びにプログラム |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101879515B1 (ko) * | 2016-12-19 | 2018-07-18 | 한국해양과학기술원 | 쌍축선의 실시간 진동 정보와 프로펠러 회전각 조절을 통한 변동압력 저감 방법 |
| KR101884534B1 (ko) * | 2016-12-19 | 2018-08-01 | 한국해양과학기술원 | 쌍축선의 프로펠러 회전각 조절을 통한 선체 변동압력 저감 방법 |
| CN110489902B (zh) * | 2019-08-26 | 2022-07-29 | 安徽工业大学 | 一种螺旋桨空化尾流精细特征多元统计建模方法 |
| CN114611214B (zh) * | 2022-03-15 | 2024-05-07 | 中国舰船研究设计中心 | 基于螺旋桨高效建模的敞水性能仿真方法 |
| CN115169261A (zh) * | 2022-07-14 | 2022-10-11 | 哈尔滨工业大学 | 一种复合材料螺旋桨的空化斗计算方法、存储介质及设备 |
| CN118243343B (zh) * | 2024-05-28 | 2024-08-06 | 中国人民解放军海军工程大学 | 一种舰船水压场模拟发生装置和模拟试验方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH1011420A (ja) * | 1996-06-20 | 1998-01-16 | Toshio Miyauchi | 物体周りの流れ推定方法 |
| JP2003200890A (ja) * | 2001-12-28 | 2003-07-15 | Hyundai Heavy Industries Co Ltd | チップ部分に曲線型レーキを形成した船舶用プロペラ |
| JP2003320989A (ja) * | 2002-05-08 | 2003-11-11 | Mitsubishi Heavy Ind Ltd | 流場計算方法、及び流場計算用プログラム |
| JP4184401B2 (ja) * | 2006-10-20 | 2008-11-19 | 川崎重工業株式会社 | ジェット騒音の音源モデル化方法、ジェット騒音の解析方法および航空機の設計方法 |
-
2008
- 2008-07-31 JP JP2008198257A patent/JP4838829B2/ja not_active Expired - Fee Related
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010038556A (ja) * | 2008-07-31 | 2010-02-18 | Mitsubishi Heavy Ind Ltd | プロペラの圧力変動推定装置及び方法並びにプログラム |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2010036600A (ja) | 2010-02-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4838829B2 (ja) | プロペラによる船体表面の変動圧力推定装置及び方法並びにプログラム | |
| Offermans et al. | On the strong scaling of the spectral element solver Nek5000 on petascale systems | |
| Biancolini et al. | Static aeroelastic analysis of an aircraft wind-tunnel model by means of modal RBF mesh updating | |
| EP3851996A1 (en) | Automotive body vibration performance optimization analysis method and device | |
| Kim et al. | Development of a numerical simulation tool for efficient and robust prediction of ship resistance | |
| Chan et al. | Prediction of vibration response levels of mistuned integral bladed disks (blisks): robustness studies | |
| Desai et al. | Uncertainty quantification of a nonlinear aeroelastic system using polynomial chaos expansion with constant phase interpolation | |
| Denimal et al. | Topological optimization of under-platform dampers with moving morphable components and global optimization algorithm for nonlinear frequency response | |
| Yıldız et al. | URANS prediction of roll damping for a ship hull section at shallow draft | |
| He et al. | Two-scale methodology for URANS/large eddy simulation solutions of unsteady turbomachinery flows | |
| Woolliscroft et al. | A fast-running CFD formulation for unsteady ship maneuvering performance prediction | |
| JP4973296B2 (ja) | 振動対策支援システム、振動対策支援方法、振動対策支援プログラム | |
| CN105989205A (zh) | 飞行器表面脉动压力的确定方法 | |
| JP6308117B2 (ja) | 振動評価方法 | |
| JP4744542B2 (ja) | 設計支援装置 | |
| CN107817300A (zh) | 使用基于声学特征的技术来确定结构内的裂纹扩展长度的系统和方法 | |
| Bakica et al. | Hydroelastic analysis of pre-swirl stator | |
| CN102270249B (zh) | 识别零部件特征频率的方法 | |
| Owen et al. | Recent ship hydrodynamics developments in the parallel two-fluid flow solver Alya | |
| Lu et al. | A rapid response calculation method for symmetrical floating structures based on state–space model solving in hybrid time-Laplace domain | |
| JP5573446B2 (ja) | 液体と接触する構造体の応答量推定方法と装置 | |
| Sachdev et al. | A mesh adjustment scheme for embedded boundaries | |
| CN113836773A (zh) | 声固耦合系统结构振动响应预示方法、系统及可存储介质 | |
| JP2022136412A (ja) | 性能予測装置、性能予測方法及び性能予測プログラム | |
| Totounferoush et al. | A data-based inter-code load balancing method for partitioned solvers |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20100212 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20110816 |
|
| 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: 20110913 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20110930 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20141007 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20141007 Year of fee payment: 3 |
|
| LAPS | Cancellation because of no payment of annual fees |