斜浪规则波中船舶运动与载荷的CFD-FEM模拟分析
CFD-FEM simulation of ship motions and loads in oblique regular waves
-
摘要:目的 旨在建立基于CFD-FEM双向流固耦合的船舶水弹性方法,预报船舶在斜浪规则波中的运动与波浪载荷,并充分考虑不对称载荷与结构响应及非线性载荷因素的影响。方法 首先,基于CFD软件建立斜浪规则波计算域,在流场域中求解三维N-S方程,实时求解非线性波浪载荷;然后,建立船体梁与船体外表面耦合的有限元模型,通过在船体梁上均匀施加集中质量和横摇惯性矩,实现船模总质量和横摇转动惯量与实船的相似;最后,通过船体外表面实时进行信息传递,以在2个求解器之间进行双向迭代耦合,进而针对非对称波浪载荷及船体结构响应开展模拟预报。结果 与试验结果的比较表明,该数值方法能够较为准确地预报纵摇、垂向弯矩和扭矩的趋势与幅值;在高海况下,波频的水平弯矩和扭矩会随波高线性增加,而在高频信号下则表现出明显的非线性增大。结论 采用所建立的CFD-FEM方法能够准确预报船舶在不对称波浪下的运动与载荷响应,可为复杂海况下船舶的水动力与结构载荷性能预报评估提供一种新的途径,且在恶劣海况下,船体水弹性效应对垂向弯矩、水平弯矩和扭矩均具有显著影响。Abstract:Objectives With the continuous development of modern ships with increasing size and speed, and the application of high-strength steel, the problem of ship hydroelastic vibration has become more prominent due to the reduction in hull natural frequencies. Springing and whipping loads significantly contribute to the wave loads of large ships. However, existing research on ship wave loads and hydroelasticity mainly focuses on the symmetric responses of ships sailing in head regular waves. In actual sea conditions, ships encounter waves from various directions. Due to the complexity of analyzing ship asymmetric responses, research on wave loads and hydroelasticity of ships in oblique waves remains limited. It is of great significance to accurately predict the asymmetric wave loads and structural responses of ships under extreme wave conditions. This study aims to develop a hydroelasticity analysis method for ships using two-way CFD-FEM fluid-structure coupling to predict the motions and wave loads of ships in oblique regular waves, accounting for asymmetric loads, structural responses, and nonlinear load effects. This method can provide a new approach for evaluating the hydrodynamic and structural load performance of ships in oblique sea conditions, enhancing the understanding of hydroelastic effects on ships.Methods The methodology employed in this research involves several steps. First, a computational domain for oblique regular waves is established using CFD software. The three-dimensional N-S equations are solved within the fluid domain to calculate the nonlinear wave loads on the model-scale ship. The computational domain of the numerical wave tank consists of a background region and an overset region. The Euler Overlay method is used to generate fifth-order Stokes waves. Second, a finite element model integrating the hull beam and the ship hull is developed. The massless hull surface is represented using shell elements, while the backbone beam is modeled with 3D uniform beam elements. The total mass and roll moment of inertia of the ship model are similar to those of the full-scale ship. To ensure the longitudinal weight distribution and roll moment of inertia in the FE model to be consistent with the experimental model, concentrated mass and moment of inertia are added to the reference points in two stages. Finally, a two-way fluid-structure coupling analysis is performed. Both the motions and structural deformations of flexible structure derived from the FEA are fed back to the CFD solver to update the hydrodynamic grid data. The fluid loads on the deformed structure, calculated using CFD with a morphing grid technique, are then applied to the structural FE model for the subsequent FEA.Results The results show that the numerical method is effective, as confirmed by the CFD grid and time step sensitivity analysis and the comparison with free roll decay test results. It can accurately predict the trends and amplitudes of pitch motion, vertical bending moment (VBM), and torsional moment (TM) compared to experimental results. For instance, in the motion response analysis of a ship in oblique waves, the numerical simulation results of heave, pitch, and roll are generally consistent with the experimental results, with a pitch peak error within 8%. In terms of load analysis, the VBM and TM calculated by the numerical method align well with the experimental trends, though there are some differences. Under extreme sea conditions, the high-frequency components caused by slamming primarily contribute to the total bending moments. The wave-frequency horizontal bending moment (HBM) and TM increase linearly with the wave height, while the high-frequency components exhibit significant nonlinear growth. Under typical extreme sea conditions, the HBM at typical sections is comparable to the VBM loads.Conclusions In conclusion, the established CFD-FEM method can accurately predict the motion and load responses of ships under asymmetric waves. It provides a novel approach for evaluating the hydrodynamic and structural load performance of ships in oblique sea conditions. The research also reveals the significant impact of hydroelastic effects on VBM, HBM, and TM under severe sea conditions. This study provides valuable insights into ship design and performance evaluation in complex sea environments, promoting the development of ship hydroelasticity research.