Advanced
Hydro-elastic analysis of marine propellers based on a BEM-FEM coupled FSI algorithm
Hydro-elastic analysis of marine propellers based on a BEM-FEM coupled FSI algorithm
International Journal of Naval Architecture and Ocean Engineering. 2014. Sep, 6(3): 562-577
Copyright © 2014, The Society of Naval Architects of Korea
This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Published : September 30, 2014
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Hyoungsuk Lee
Maritime Research Institute, Hyundai Heavy Industries, Co., Ltd., Ulsan, Korea
Min-Churl Song
Industrial Technology Institute, Hyundai Heavy Industries, Co., Ltd., Ulsan, Korea
Jung-Chun Suh
Dept. of Naval Architecture and Ocean Engineering & Research Institute of Marine Systems Engineering, Seoul National University, Seoul, Korea
Bong-Jun Chang
Maritime Research Institute, Hyundai Heavy Industries, Co., Ltd., Ulsan, Korea
Abstract
A reliable steady/transient hydro-elastic analysis is developed for flexible (composite) marine propeller blade design which deforms according to its environmental load (ship speed, revolution speed, wake distribution, etc.) Hydro-elastic analysis based on CFD and FEM has been widely used in the engineering field because of its accurate results however it takes large computation time to apply early propeller design stage. Therefore the analysis based on a boundary element method-Finite Element Method (BEM-FEM) Fluid-Structure Interaction (FSI) is introduced for computational efficiency and accuracy. The steady FSI analysis, and its application to reverse engineering, is designed for use regarding optimum geometry and ply stack design. A time domain two-way coupled transient FSI analysis is developed by considering the hydrodynamic damping ffects of added mass due to fluid around the propeller blade. The analysis makes possible to evaluate blade strength and also enable to do risk assessment by estimating the change in performance and the deformation depending on blade position in the ship’s wake. To validate this hydro-elastic analysis methodology, published model test results of P5479 and P5475 are applied to verify the steady and the transient FSI analysis, respectively. As the results, the proposed steady and unsteady analysis methodology gives sufficient accuracy to apply flexible marine propeller design .
Keywords
INTRODUCTION
- Research background
A marine propeller operates in the complex flow field called the ship’s wake and the transition of flow angle of attack leads to rapid changes in pressure distribution on the propeller blade. As the trend of high power engines follows increases in vessel size, a marine propeller is being exposed to various risks including cavitation erosion, pressure fluctuation-induced vibration, underwater radiation noise, etc.
Compared with propellers made of metal, those made of composite materials - a combination of two or more physically / chemically heterogeneous materials - have superior properties for managing these risks. Composite materials are light weight, and have high stiffness and better sound absorption ability ( ITTC, 2005 ). Despite their higher costs, the excellent manufacturing and maintenance ability of an additive production process is expanding the application of composite materials to all industrial areas, including aviation, environment/energy, information, electricity, automobile and household supplies.
Attempts are being made to use composite materials in the marine propeller blades of warships, submarines, and torpedoes for two reasons. First, composites can improve chance of survival because their good sound absorption ability reduces noncevitating propeller noise. Second, their blade deformation characteristics of the wake distribution delays cavitation inception and reduces cavitation volume. Their superior properties have led to their expanded application in yachts, tugs, and fishing boats. Composite materials also make it possible to reduce A/S costs for preventing heat damages in shaft systems, and to simplify the vessel stern structures due to their light weight - less than half the weight of Nickel-Aluminum-Bronze (NAB) alloy ( ITTC, 2005 ).
The amount and geometry of a composite material’s deformation from the acting load can be controlled with different combinations of laminated material, stacking order, and ply angle. Thus, if its deformable characteristics can be properly considered at the design stage, a composite propeller can reduce fuel oil consumption in its operating profile compared with an NAB propeller which does not change shape ( Young, 2007a ; 2007b ; 2008 ).
While recognizing these virtues, the difficulty is that it is not easy to accurately and practically simulate the hydro-elastic behavior of the composite marine propellers. These simulation difficulties act as a barrier to the accurate prediction of propulsion performance, initial geometry decisions, and strength evaluations at the propeller design stage.
- Previous studies
Inviscid numerical methodologies for calculating the flow field around propeller blades include the lifting line method, the lifting surface method (or vortex-lattice method, VLM), and the boundary element method (BEM). These inviscid methods are widely used in propeller performance prediction because of their short computing times and general usability. Computational fluid dynamics (CFD) analyses (RANS, LES, etc.) that consider the viscous effect tend to use a propeller flow analysis, which requires increased computing power. Even though the viscous flow analyses give a more accurate solution, BEM-based FSI analysis is performed in this study due to the considerations of computing time for the evaluation of propeller design applications.
The lifting line theory, proposed by Prandtl (1921) , and the VLM, which was applied to the propeller by Kerwin et al. (1978) , are useful methods for initial propeller design due to their ease of calculation. They have limitations, however, when considering three-dimensional blade geometry, BEM is introduced to achieve a more accurate FSI analysis. BEM, also called the panel method, was adapted to the propeller by dividing the potential basis and velocity basis methods of Hess and Valarezo (1985) , Lee (1987) , Kerwin et al. (1987) , Hoshino (1989) , etc. After Lee (1987) showed the excellent performance of the potential basis method, this method has been mainly used in the field of propeller analysis. Hsin (1990) and Suh et al. (1992) proposed the Kutta condition and Greeley and Kerwin (1982) and Lee (1987) introduced the flow-adapted grid to include a propeller wake roll-up model for improved convergence in the BEM analysis. These methods apply to the BEM code, KPA14, used in this study.
Various other numerical and experimental studies have also been conducted to analyze composite material (or elastomer) propellers. Lin and Lin (1996) simulated the elastic behavior of marine propellers in noncavitating conditions with a VLMFEM interaction. Young (2007a ; 2008) showed a hydro-elastic response of a composite propeller in steady/unsteady conditions by using BEM-FEM FSI and validated the results with model test results. In particular, Young’s study conducted BEM analysis to predict the hydrodynamic damping effect for unsteady behavior such as blade “fluttering” induced by a change in inflow conditions at the hull behind the wake field. Chen et al. (2006) evaluated the behavior of the composite marine propeller in four-cycle inflow conditions by performing model tests and a numerical simulation. Motley et al. (2009) and Blasques et al. (2010) studied optimum material stacking to improve propulsion efficiency and Jang et al. (2012) proposed a reverse engineering procedure for the initial design of an elastomer propeller considering blade deformation. Lee et al. (2009) , Lee et al. (2012) and Lee et al. (2013) have carried out the study of the relationship between manufacturing and performance of flexible (composite) propeller.
- Objective
The hydro-elastic analysis of a flexible marine propeller, which is a simplified orthotropic-homogeneous material model of a composite propeller, has been developed to accurately and efficiently predict the propeller’s propulsion performance and evaluate its strength. For practical use in the propeller design stage, the inviscid potential-based BEM code, KPA14, is used to calculate the loads on the blades. A commercial Finite Element Analysis (FEA) solver (Abaqus Standard) is introduced for predicting blade deformation and evaluating propeller strength. An interface between these two solvers is developed to analyze the hydro-elastic behavior of the composite-material propeller, and the accuracy and usability of the developed algorithm are confirmed by comparing our results with the model test results in previously published literature.
FORMULATION AND NUMERICAL METHODS
The flow field around the propeller blade, the structural deformation, and their relationship at the blade surface are numerically formulated for the FSI analysis to predict the performance of the marine propeller.
- Governing equation of the flow field around the propeller blades
As mentioned in the previous section, the potential basis analysis, from the studies of Hsin (1990) , Suh et al. (1992) , Greeley and Kerwin (1982) and Lee (1987) , are applied to analyze the steady or unsteady flow fields. We assume that the flow field consists of the three-dimensional closed boundary S and the unit normal vector 𝐧 toward the flow to formulate the flow around the propeller blades. The boundary S is composed of the blade and hub SB , the wake SW , and the outer boundary S . The propeller is rotating in the incompressible, inviscid, and irrotational flow with inflow velocity U . The general notation of object in fluid for the application of Green’s theorem is shown in Fig. 1 .
PPT Slide
Lager Image
General notation of object in fluid for the application of green’s theorem.
The perturbed velocity potential φ satisfies the Laplace equation.
PPT Slide
Lager Image
When applying the boundary condition to the boundaries, the boundary value problem can be formulated and the kinematic boundary condition must be satisfied on the propeller blade surface S as follows:
PPT Slide
Lager Image
Assume that wake surface SW has no thickness and there is no normal velocity jump, then
PPT Slide
Lager Image
The Kutta condition of Eq. (4) is applied to the unique circulation at the trailing edge of the blades.
PPT Slide
Lager Image
Namely, the flow velocity is finite at the trailing edge and at the outer boundary, the perturbed velocity induced by the propeller is vanishing far away from the blade surface as follows:
PPT Slide
Lager Image
The disturbed velocity potential φ can be expressed as Eq. (6) according to the governing equation because it should satisfy Green’s 2 nd identity for all points on the blade surface.
PPT Slide
Lager Image
where
  • S∞:Sblade+Shub
  • p: point induced velocity potential is calculated
  • q: point positioned singularity
  • R( p; q): distance betweenpandq
The geometry of the propeller blades, hub, and wake could be discretized as a panel to numerically solve the governing equation in Eq. (7).
PPT Slide
Lager Image
The effective coefficients ai,j and bi,j are the velocity potentials at the i th panel induced by the unit strength dipole and the source distributed at the j th panel, respectively. The wake effective coefficient Wi,m is the velocity potential at the i th panel induced by the unit strength dipole at the m th strip at the wake panel. (Δ φ ) m is the difference in potential between the upper and lower side of the trailing edge and σj is the source strength on the blades. The diagram of discritized surface panel of blade, hub and wake and coordinate system is shown in Fig. 2 .
PPT Slide
Lager Image
Discretized surface panel of blade, hub and wake and coordinate system.
To solve the unsteady transient flow problem defined by the time domain, a steady analysis for the key blade is performed for every time step and the singularity of the other blades is updated for every step shown in Eq. (8).
PPT Slide
Lager Image
where
PPT Slide
Lager Image
and variable n means time step, superscript k means the indices of blades, i.e. k =1 means the key blade.
The pressure acting on the blade surface can be obtained by using the velocity potential derived from the discretized integral equation in Eqs. (7) or (8) as follows:
PPT Slide
Lager Image
where the operator ∂ / ∂ t in the inertial frame into ( ∂′/ ∂t − V F ⋅∇) in the moving frame.
- Structural equations of motion for blade deformation
The equation of motion for the structural deformation corresponding to the propeller blade fixed coordinate is as Eq. (10),
PPT Slide
Lager Image
where M s is the mass matrix, C s is the damping matrix, K s is the stiffness matrix which represents the material mass, damping and stiffness of the composite propeller blade. The values of mass, damping and stiffness in the validation have use of published experimental results.
PPT Slide
Lager Image
,
PPT Slide
Lager Image
, d are acceleration, velocity, and displacement for the computing point, respectively.
F ST in Eq. (11) represents the total of all loads acting on the blade structure. These consist of the hydrodynamic pressure load F hp and the body force due to rigid body rotation, the Coriolis force F cori and the centrifugal force F cent , and the coupling component force F fs acting on the interface between the fluid and the structure. The summation of the structural loads except F fs can be written as F s .
PPT Slide
Lager Image
The structural analysis of the hydro-elastic response is conducted with the commercial FEA solver Abaqus 6.12. Since the implicit method is adapted for both the steady and unsteady analyses, the static option of the Abaqus Standard is used for the steady state analysis and the dynamic option is used for the transient analysis.
The Newmark integration method of Eqs. (12) and (13) is applied to solve the nonlinear dynamic response problem and the model response at time t + Δ t .
PPT Slide
Lager Image
PPT Slide
Lager Image
where parameters β and γ are the value-related stability and the accuracy of the convergence of analysis, respectively. 0.25 and 0.5 are applied for each parameter, providing second-order analysis accuracy.
The 8-node continuum shell element is used in the structural analysis. This element can reflect the three-dimensional geometry of the propeller blade and has the merit of giving a shell element to the plate-like structure, as well as having automatic mesh generation and ease in applying the material properties according to the ply.
- Fluid-structure interaction
Assuming that deformation of the blade for the inflow variation is relatively large and that the change in the hydrodynamic characteristics due to geometrical alteration is not small, the FSI effect must be applied for the estimation on each hydrodynamic force and structural deformation. We propose using a two-way coupled FSI analysis that includes structural nonlinear analysis because the amount of blade deformation is not small to assume linear according to Young (2008) . To analyze strictly for the elastically vibrating response due to the changing flow field in the ship’s wake, the “wet” condition (considering not only hydrodynamic pressure distribution on blade but including hydrodynamic added mass/damping effect) must be considered. The acoustic fluid medium finite element model is adapted to consider this effect. This acoustic medium model was applied to ship hull vibration in contact with water by Paik (2010) .
Let the fluid field around propeller blade be compressible and inviscid then the equation of equilibrium is as follows:
PPT Slide
Lager Image
where p is the hydrodynamic pressure, ρ is the density of fluid, γ is the volumetric constant (volumetric drag) for dissipation, and
PPT Slide
Lager Image
and
PPT Slide
Lager Image
represent velocity and acceleration of fluid particles, respectively.
Fluid volumetric strain can be represented as in Eq. (15) by the definition of the relationship between the bulk modulus and pressure.
PPT Slide
Lager Image
The equation of motion for hydro pressure can be derived from Eqs. (14) and (15).
PPT Slide
Lager Image
The equation of motion that considers the coupling of the boundary between the fluid and the structure can be expressed as the vector form.
PPT Slide
Lager Image
where Mf , Cf , Kf and S fs are 1 / Kf , γ / ( ρ Kf ), −1 / ρ 2 and coupling matrix respectively. And p represents a pressure vector in the fluid domain, drepresents a unit vector on the fluid-structure boundary. Eq. (18) is derived from Eqs. (10) and (11) for structural deformation and then the system equation of FSI (Eq. 19) is obtained by substitution with Eq. (17).
PPT Slide
Lager Image
PPT Slide
Lager Image
The tetra mesh (AC3D4 known as an element type of Abaqus Standard) of the acoustic medium domain is composed to perform structural analysis with consideration to the hydrodynamic damping due to added mass are shown in Fig. 3 .
PPT Slide
Lager Image
Grid system for structural analysis in acoustic medium (‘wet’ condition).
METHODOLOGY FOR HYDRO-ELASTICTY ANALYSIS OF A PROPELLER BLADE
Based on the formulated mathematical model, we implement a steady and unsteady transient FSI algorithm to predict the performance and evaluate the structure of the flexible propeller.
- Basic concept of FSI analysis
A basic interface structure of the two-way coupled FSI, shown in Fig. 4 , consists of a BEM panel code (KPA14) for calculating the hydrodynamic force acting on the blade surface and a nonlinear FEM solver (Abaqus Standard) for calculating the structural deforming response, rotational body force (Coriolis force and centrifugal force), and added mass damping of the unsteady transient flow condition.
PPT Slide
Lager Image
Schematic diagram for propeller blade FSI.
In steady-state FSI analysis, it is not necessary to consider the change in pressure distribution and the blade fluttering effect according to the blade position, so ★ items are excluded. The convergence is determined by the progress of comparison of geometrical differences from the former step (steady) and the former revolution (transient). Thus it is considered as converge when the maximum difference of nodal point is less than 0.01% of propeller diameter in the iterative FSI analysis.
- Process for initial geometry prediction
The pattern of blade deformation depends on the properties of the materials used in manufacture and the ply stacking condition such as the main orientation, as well as the loads acting on the blade. This means that the process of predicting the initial blade geometry, which is dependent on the materials and ply stacking conditions, requires that the design geometry be determined at given propeller execution conditions.
Assuming that the optimum geometry for the propeller design condition is the “A” blade and the flexible propeller is manufactured with this geometry and is in operation, its shape would be deformed to “B”. The “B” blade has completely different geometry than the optimal design in Fig. 5 . Thus, the initial design process should be adapted to the geometry design of the marine propeller with the shape of the “C” blade, which is estimated by the reverse engineering method presented by Jang (2012) . Then, the blade would be changed to the target geometry for the existing driving conditions. The reverse engineering algorithm for blade initial design shown below Fig. 6 .
PPT Slide
Lager Image
Blade deformation aspects of design load condition and initial design considering deformation.
PPT Slide
Lager Image
Reverse engineering algorithm for blade initial design.
- Unsteady transient FSI procedure
Time domain-based FSI analysis of the propeller’s blade position in the wake field is required to more accurately predict the performance and evaluate the structural strength. Rotating blades in the nonuniform inflow flow field have a “fluttering” vibration behavior caused by the change in pressure distribution and blade deformation from the response of the pressure load on the surface. Unlike the aviation propeller and wind turbine blade which act in low density air, the marine propeller cannot ignore the damping effect of the added mass on the blade surface and vibrates in the ship’s wake. Thus, in this study, we considered the hydrodynamic damping effect of the added mass by applying the acoustic medium around the propeller blades.
To perform unsteady transient analysis (dynamic analysis of Abaqus Standard), it is impossible to carry out the prediction of deformed shape at T + dt priory in the transient analysis. Therefore, considering the amount of infinitesimally small time step dt , and assuming pressure varies with linearly at the blade surface, the transient analysis is performed by applying blade shape of the time T + dt with the shape at time T (final shape of the deformation in the previous time step) in the way that considering only changes in pressure variation of blade position in wake ( Fig. 7 ).
PPT Slide
Lager Image
Transient FSI algorithm for composite propeller.
Steady analysis for the initial convergence is performed prior to the unsteady analysis in order to reduce computing time since the transient analysis can start with the result. The convergence check of the total unsteady analysis is determined by comparing the geometrical difference between the present and the former revolution positions.
VALIDATIONS
A comparison and review of experimental and computational results in the literature were conducted to verify the BEM-FEM-based steady and transient FSI analysis presented in this study. Since published references that contain both measured experimental data as well as the geometric/material properties of propellers is limited, the propellers validated in each analysis are different.
- Steady state analysis
Chen et al. (2006) and Young (2008) present results for their experiments and FSI analysis of a six-bladed 610 mm diameter carbon fiber reinforced plastic propeller, manufactured at A.I.R. Fertigung-Technologie, GmbH in Hohen Luckow, Germany, and tested in the Naval Surface Warfare Center, Carderock Division (NSWCCD). In particular, Young (2008) conducted the validation of the proposed analytical method using the BEM-FEM FSI algorithm to compare the thrust, torque, and efficiency of a full-scaled open water model test for a “rigid” propeller P5471 and a “flexible” propeller P5479.
In this study, we validate this analytical method by comparing our results with Young’s (2008) experimental and analytical results for propeller P5479. The geometry of propeller P5479 is shown in Fig. 8 . This is a useful exercise for our verification because this propeller has significant hydro-elastic behavior. Isotropic homogeneous material (not a real ply stack model) is applied in the analysis. The material properties extracted by Liu et al. (2007) are E 1 = 138.9 GPa , E 2 = 10.2 GPa , G 12 = 5.7 GPa , and ν 12 = 0.3 and the main fiber orientation is applied at + 32° (see the definition of fiber orientation in Young (2008) ).
PPT Slide
Lager Image
Geometry of propeller P5479.
The mesh size of the BEM analysis is 40 × 20 for each chord (from face to back side) and span direction. We apply a cosine spacing for chord direction and a half cosine spacing for span direction to dense toward tip to faithfully reflect the leading/ trailing edge and tip geometry of the blades. A total of 6400 elements (40 × 40 × 4 for chord, span, and thickness direction) are composed for the structural finite element analysis. We introduce an inverse distance weighting algorithm for the three-dimensional interpolation to share the different information at the calculating points.
FSI analyses are conducted under two operating conditions, RPM = 909 and 454, and they are compared to the NSWCCD test and to Young’s (2008) calculations at the advance ratio ( J ) 0.5-0.9. Fig. 9 shows that our study results and Young’s (2008) analysis with the full-scale model test. It seems that the results of this study are quite good for operation advance ratio range. But, from the high advance ratio region, it contains a relatively large degree of error. It can be estimated that the reason of these errors is due to the feature of KPA14, the used BEM code, does not fully considering the viscous characteristic of Kaplan type propeller which has a finite tip chord length. Thus, it is expected to be improved by tuning the viscous correction factor in the BEM code with more tested data. The accuracy of the proposed FSI algorithm at the design advance ratio, J = 0.66 is sufficient for designing a composite propeller.
PPT Slide
Lager Image
Comparison of model test and BEM-FEM FSI prediction of propeller P5479.
- Unsteady transient analysis
A marine propeller operates in the complex flow field condition known as the ship’s wake. The composite material marine propeller, in particular, must undergo hydro-elastic response analysis for unsteady conditions because its shape changes according to the blade position in the wake. Therefore, consideration of its structural safety, such as the risk of delamination due to fatigue and accurately predicting its performance, cannot be ignored.
The damping induced by the hydrodynamic added mass significantly affects the elastic response of the propeller blade in water. Thus hydro-elastic analysis that considers the so called “wet” condition effect should be conducted. In the research of Chen et al. (2006) , Young (2007a ; 2008) , and Motley et al. (2009) , this elastic vibration effect was estimated using the potential- based BEM and the results were applied directly to the structural mass and damping coefficients.
In this study, however, the added mass effect is considered by adapting the acoustic medium field and applying it as a pressure component on the blade surface directly. We compare the results from our proposed method with the measured and calculated results of Chen et al. (2006) .
The propeller used in the analysis is a P5475 ( Fig. 10 ), which deforms in actual conditions significantly enough to validate the analysis algorithm and the Laser Doppler Velocimetry (LDV) measured result of the pitch deformation published in Chen et al. (2006) . Material properties of the P5475 are E 1 = 139.0 GPa , E 2 = 8.06 GPa , G 12 = 4.88 GPa , and ν 12 = 0.28 ( Young, 2007b ).
PPT Slide
Lager Image
Geometry of propeller P5475.
A wake screen was adapted to simulate the 4-cycle wake flow for measuring the hydro-elastic behavior in the unsteady flow field at NSWCCD, and the 4-cycle wake flow condition is applied to the FSI analysis. The mean velocity of inflow in the tunnel was 5.03 m/s and velocity distribution measured by LDV is the dotted line in Fig. 11 . The wake applied in this study is the solid line in Fig. 11 and it has minor differences from the measured velocity.
PPT Slide
Lager Image
Inflow measured by LDV (NSWCCD) and applied to FSI analysis (present).
Pitch variations were measured by video cameras at two points, 10° and 55°. Propeller rotation was 750 RPM and the measured difference of pitch angle at the blade tip between the deformed and undeformed shape according blade position is represented by the squared symbol in Fig. 12 . The calculation results from Chen et al. (2006) using the VLM-FEM FSI method and from this study’s BEM-FEM FSI method (solid line) under the same conditions is plotted also. To examine the effect of hydrodynamic added mass damping, the “dry” condition only material damping is displayed in the dotted line.
PPT Slide
Lager Image
Pitch change for P5475 (measured, predicted by Chen et al. and present).
The proposed unsteady transient algorithm is more accurate than the VLM-FEM FSI method regarding mean deformation. The difference between an analysis result that considers added mass due to the vibrating blade (“wet” condition) and a result that considers only material/structural damping (“dry” condition) is not significant. But considering that the purpose of the unsteady FSI analysis of a marine propeller in a wake field is for accurate performance prediction and structural safety evaluation, we believe that the proposed unsteady transient analysis under “wet” conditions can give more accurate and realistic prediction results for deformed geometry and variations.
The calculated variations of the thrust and torque coefficients for the key blade in the four-cycle wake flow and the measured data from the composite propeller P5475 at J = 0.66/RPM = 780 for open water conditions ( Chen et al., 2006 ) is shown in Fig. 13 . Even though differences exist in the inflow conditions, four-cycle and uniform, the average values of thrust and torque variations from the unsteady FSI analysis and the measured values from uniform flow at the same advance ratio seem very close. These results show that the proposed practical methodology is sufficiently accurate for use in composite marine propeller design with the validation results from steady state analysis.
PPT Slide
Lager Image
Thrust coefficient and torque coefficient for key blade of P5475 (present).
- Natural frequency analysis for dynamic stability evaluation
Natural frequency analyses of propeller P5479 rotating at 909RPM for “dry” and “wet” conditions are conducted and the analyzed natural frequencies of each mode are compared with the measured (“wet” condition, mode 1) and calculated results of Young (2008) in Fig. 14 . And the blade shapes under “wet” conditions for each mode are shown in Fig. 15 .
PPT Slide
Lager Image
Natural frequency prediction.
PPT Slide
Lager Image
Mode shapes for P5479 at ‘wet’ condition.
CONCLUSIONS
Since a flexible (or composite material) propeller is affected not only by the pressure changes acting on the blade but also the body forces from the rotational motion, such as the centrifugal force, the Coriolis force, and the “fluttering” hydro-elastic behavior in wake field; it is not easy to predict the propulsion performance, the initial geometry with respect to deformation, and the structural strength at the design stage.
The BEM-FEM-based practical hydro-elastic algorithm for a flexible marine propeller is developed for blade and ply orientation design to address performance prediction and structural evaluation. The study contents are as follows:
  • 1) BEM code does not consider viscous effect and effective wake. However, it still gives sufficient accuracy and short computing time to apply propeller design stage. Therefore BEM-FEM FSI methodology is proposed for practical application to flexible marine propeller design instead of the more elaborate CFD-FEM FSI analysis.
  • 2) Panel code (KPA14) is used for hydrodynamic analysis to predict pressure loads on the propeller blade and a commercial FEA tool (Abaqus 6.12 Standard, Static and Dynamic option) is used for nonlinear analysis to calculate blade deformation and evaluate strength. An interface code is developed to connect these two analysis tools and perform FSI analysis.
  • 3) A three-dimensional continuum shell element is introduced for the FEA, which has merits in its easy modeling of composite material ply and its ability to reflect the structural behavior of objects having three-dimensional geometry.
  • 4) The hydrodynamic damping effect of added mass, which is induced by the elastic behavior of the propeller blade, are applied to the FEM for acoustic fluid mediums and as composing system equations of motion for pressure on the surface and the structural deformation of blade geometry, unlike existing studies.
  • 5) A validation study for the proposed steady and unsteady transient algorithm for composite propeller FSI analysis is conducted by comparing our study results with published experimental and analytical results on flexible propellers composed of orthotropic-homogeneous materials. As a result, the usefulness of the proposed practical methodology was confirmed.
Further study of the application of composite material ply structure in the FSI methodology will be conducted for further analysis of the composite material propeller. And future research will address the optimum stack algorithm that satisfies both structural strength and performance optimization for off-design operation conditions. In addition, CFD-FEM FSI methodology will be developed to improve hydro-elastic analysis accuracy.
Acknowledgements
This research was supported partly by the Industrial Convergence Strategic Technology Development Program (No. 1004- 4499) funded by the Ministry of Trade, Industry and Energy (Korea).
References
Blasques J.P. , Berggreen C. , Anderson P. 2010 Hydro-elastic analysis and optimization of a composite marine propeller Marine Structures 23 22 - 38
Chen B.Y.H. , Neely S.K. , Michael T.J. , Gowing S. , Szwerc R.P. , Buchler D. , Schult R. 2006 Design, fabrication and testing of pitch-adapting (flexible) composite propellers SNAME Propellers/Shafting Symposium Williamsburg,VA 9 January 2006 8/1 - 11
Greeley D.S. , Kerwin J.E. 1982 Numerical methods for propeller design and analysis in steady flow Transactions of Society of Naval Architects and Marine Engineers 90 415 - 453
Hess J.L. , Valarezo W.O. 1985 Calculation of steady flow about propellers by means of a surface panel method Journal of Propulsion and Power 1 (6) 470 - 476
Hoshino T. 1989 Hydrodynamic analysis of propellers in steady flow using a surface panel method Journal of the Society of Naval Architects of Japan 165 55 - 70
Hsin C.Y. 1990 Development and analysis of panel method for propellers in unsteady flow. PhD thesis MIT, Dept. of Ocean Engineering
2005 The Propulsion Committee Edinburg Final report and recommendations to the 24th ITTC 4-10 September 73 - 136
Jang H.G. , Noh I.S. , Hong C.H. , Lee C.S. 2012 Design algorithm of flexible propeller by fluid-structure interactive analysis Journal of the Society of Naval Architects of Korea 49 528 - 533
Jang H.G. , Noh I.S. , Hong C.H. , Lee C.S. 2013 Design algorithm of flexible propeller by fluid-structure interactive analysis Proceeding of the 12th International Symposium on Practical Design of Ships and Other Floating Structures Changwon City 20-25 October 2013 1158 - 1164
Kerwin J.E. , Lee C.S. 1978 Prediction of steady and unsteady marine propeller performance by numerical lifting surface theory Transactions of Society of Naval Architects and Marine Engineers 86 218 - 253
Kerwin J.E. , Kinnas S.A. , Lee J.T. , Shih W.Z. 1987 A surface panel method for the hydrodynamic analysis of ducted propellers Transactions of Society of Naval Architects and Marine Engineers 95 93 - 122
Lee J.T. 1987 A potential based panel method for the analysis of marine propellers in steady flow. PhD thesis MIT, Dept. of Ocean Engineering
Lee G.H. , Jang H.G. , Lee C.S. , Nho I.S. , Lee S.G. , Hyun B.S. 2013 Effective of manufacturing accuracy of flexible propeller on the open water performance Journal of the Society of Naval Architects of Korea 50 349 - 354
Lee S.G. , Byun J.H. , Paik B.G. , Hyun B.S. 2009 Production & performance assessment of composite materialflexible propeller Journal of the Society of Naval Architects of Korea 46 667 - 674
Lee S.G. , Nam J.H. , Hyun B.S. , Paik B.G. , Lee C.S. , Jang H.G. , Nho I.S. 2012 Enhanced manufacturing and performance analysis of flexible composite propeller Journal of the Society of Naval Architects of Korea 49 521 - 527
Lin H. , Lin J. 1996 Nonlinear hydroelastic behavior of propellers using a finite element method and lifting surface theory Journal of Marine Science and Technology 1 (2) 114 - 124
Liu Z. , Young Y.L. 2007 Utilization of deformation coupling in self-twisting composite propellers Proceedings of 16th International Conference on Composite Materials Kyoto, Japan 8-13 July 2007
Motley M.R. , Liu Z. , Young Y.L. 2009 Utilizing fluid-structure interactions to improve energy efficiency of composite marine propellers in spatially varying wake Composite Structure 90 304 - 313
Paik K.J. 2010 Simulation of fluid-structure interaction for surface ships with linear/nonlinear deformations. PhD thesis University of Iowa
Prandtl L. 1921 Application of modern hydrodynamics to aeronautics, National Advisory Committee for Aeronautics Annual Report 7th NASA USA
Suh J.C. , Lee J.T. , Suh S.B. 1992 A bilinear source and doublet distribution over a planar panel, and its application to surface panel method Proceedings of 19th Symposium on Naval Hydrodynamics Seoul 23-28 August 1992
Young Y.L. 2007 Time-dependent hydroelastic analysis of cavitating propulsors Journal of Fluids and Structures 23 269 - 295
Young Y.L. , Liu Z. 2007 Hydroelastic tailoring of composite naval propulsors Proceedings of the 26th International Conference on Offshore Mechanics and Arctic Engineering San Diego, California, USA 10-15 June 2007 777 - 787
Young Y.L. 2008 Fluid-structure interaction analysis of flexible composite marine propellers Journal of Fluids and Structures 24 799 - 818