Research on coupling solution of soft-rigid multibody dynamics and numerical fluid dynamics

软刚多体动力学与数值流体动力学耦合求解研究

• 批准号: 22K04572
• 负责人: 荒木 元輝
• 金额: $ 2.66万
• 依托单位: National Institute of Maritime, Port and Aviation Technology
• 依托单位国家: 日本
• 项目类别: Grant-in-Aid for Scientific Research (C)
• 财政年份: 2022
• 资助国家: 日本
• 起止时间: 2022-04-01 至 2025-03-31
• 项目状态: 未结题
• 来源: https://kaken.nii.ac.jp/grant/KAKENHI-PROJECT-22K04572/
• 关键词: マルチボディダイナミクス; 数値流体力学; マリンオペレーション

本研究は粘性を考慮可能な数値流体力学モデル（CFD）とマルチボディダイナミクスと絶対節点座標法（ANCF）を基にした剛体、線状構造体、3次元構造体の柔剛混在多体系動力学モデルを開発し、様々なマリンイペレーションに応用可能な解析手法を確立することを目標としている。そこで初年度は線状構造体用のANCFコードを開発した。ANCFはShabanaら（1996）により提案された非線形有限要素法の一種で節点座標として節点の絶対位置および勾配で表現することにより、従来の有限要素法に比べ少ない少ない要素数で大変形を精確に表現することが可能である。本研究においては基本的にBerzeri & Shabana(2000)のモデルを用いることとした。ANCFでは仮想仕事の原理を用いてラグランジュ方程式を導出し、各ケーブル要素の運動方程式が得られるが、その中の剛性行列においては、計算負荷及び精度を考慮し、曲げ剛性行列に関しては軸ひずみが小さい場合を仮定してT1モデル、軸剛性行列にかんしては一般的なひずみ状態を扱えるL2モデルを用いることとした。また今後のマルチボディダイナミクスモデルとの連成を考慮し、微分代数方程式（DAE）となるようラグランジュ未定乗数法を適用し、運動方程式の解法としては一般化アルファ法を用いた。開発した線状構造体用のANCFコードの検証のために、ANCFコードのベンチマークとしてよく用いられる、原点でピン支持された柔軟な梁を自由落下させた時の変形を計算し、従来の研究結果とほぼ一致することを確認した。

This study は viscous を may consider な the numerical fluid mechanics モ デ ル (CFD) と マ ル チ ボ デ ィ ダ イ ナ ミ ク ス と unique node coordinates method (ANCF) を seaborne に し た rigid body, linear structure, three dimensional structure の soft just mixed in system dynamics モ デ ル を open 発 し, others 々 な マ リ ン イ ペ レ ー シ ョ ン に 応 with analytical technique may な を establish す る Youdaoplaceholder3 とを objective と て る る る. Youdaoplaceholder0 で で in the first year, <s:1> ANCFコ ドを ドを are used to develop <s:1> た. ANCF は Shabana ら (1996) に よ り proposal さ れ た nonlinear finite element method is の a で node coordinates と し て node の unique position お seaborne よ び hook with performance で す る こ と に よ り, 従 の finite elements method に less than べ な い less な で い elements number big - shaped を precise に performance す る こ と が may で あ る. The basic にBerzeri & Shabana(2000) モデ モデ を を is used in にお る とと とと た た た in this study. ANCF で は 仮 thinking about its principle of の を using い て ラ グ ラ ン ジ ュ equations derived を し, various ケ ー ブ ル elements の movement equations が ら れ る が, そ の の rigidity in the ranks of に お い て は, and び を consider し accuracy, computational load げ rigid ranks に masato し て は shaft ひ ず み が small さ い occasions を 仮 set し て T1 モ デ ル ranks, axial rigidity に か ん し て The general なひずみ state を handle えるL2モデ モデ を を is represented by なひずみ る とと とと た た た た. ま た future の マ ル チ ボ デ ィ ダ イ ナ ミ ク ス モ デ ル と の connect を し, differential algebraic equations (DAE will work) と な る よ う ラ グ ラ ン ジ ュ undetermined number of 乗 を applicable し, motion equation is の hydrolysis と し て は generalization ア ル フ を ァ method with い た. Open 発 し た linear structure body with の ANCF コ ー ド の 検 card の た め に, ANCF コ ー ド の ベ ン チ マ ー ク と し て よ く with い ら れ る, origin で ピ ン support さ れ た soft な beam を free falling さ せ た の when calculating the shape を - し 従 to の research results と ほ ぼ consistent す る こ と を confirm し た.