International Journal of Rotating Machinery

Volume 2012 (2012), Article ID 632856, 11 pages

doi:10.1155/2012/632856

Research Article

Scaling of the Transient Hydroelastic Response and Failure Mechanisms of Self-Adaptive Composite Marine Propellers

Michael R. Motley and Yin L. Young

Department of Naval Architecture and Marine Engineering, University of Michigan, Ann Arbor, MI 48109-2145, USA

Received 11 February 2012; Revised 18 June 2012; Accepted 20 June 2012

Academic Editor: Moustafa Abdel-Maksoud

Copyright © 2012 Michael R. Motley and Yin L. Young. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The load dependent deformation responses and complex failure mechanisms of self-adaptive composite propeller blades make the design, analysis, and scaling of these structures nontrivial. The objective of this work is to investigate and verify the dynamic similarity relationships for the hydroelastic response and potential failure mechanisms of self-adaptive composite marine propellers. A fully coupled, three-dimensional boundary element method-finite element method is used to compare the model and full-scale responses of a self-adaptive composite propeller. The effects of spatially varying inflow, transient sheet cavitation, and load-dependent blade deformation are considered. Three types of scaling are discussed: Reynolds scale, Froude scale, and Mach scale. The results show that Mach scaling, which requires the model inflow speed to be the same as the full scale, will lead to discrepancies in the spatial load distributions at low speeds due to differences in Froude number, but the differences between model and full-scale results become negligible at high speeds. Thus, Mach scaling is recommended for a composite marine propeller because it allows the same material and layering scheme to be used between the model and the full scale, leading to similar 3D stress distributions, and hence similar failure mechanisms, between the model and the full scale.

1. Introduction

In recent years, advanced composite materials have become an increasingly popular alternative to traditional metallic alloys for aerospace and marine applications, including rotors such as propellers and turbines. In addition to having the benefits of higher specific strength and stiffness, composites can provide improved performance over metallic alloys through exploitation of the intrinsic bend-twist coupling characteristics. The anisotropic properties of composites can be used to elastically tailor the rotor blades to achieve improved performance through passive pitch adaptation. However, the load-dependent deformation responses and complex failure mechanisms of composite blades make the design, analysis, and scaling of these structures nontrivial.

Over the last two decades, much research on composite rotors focused on the utilization of fluid-structure interactions (FSI) to improve the performance of aerospace structures, notably helicopter, aircraft, and wind turbine blades [1–7]. More recently, the use of advanced composites to improve the performance of marine rotors has been demonstrated experimentally [8–10] and numerically [11–23]. It has been shown that self-adaptive composite rotors can help to delay cavitaton, increase energy efficiency, and decrease fuel consumption when compared to rigid metallic rotors in spatially varying flows and in off-design conditions. To the knowledge of the authors, nearly all published systematic experimental studies of composite marine rotors in the open literature have been conducted in model-scale cavitation tunnel and towing tank facilities. In order to predict the full-scale, load-dependent deformation response and potential failure mechanisms of self-adaptive composite marine rotors, appropriate hydroelastic scaling laws are needed.

While hydrodynamic similarity relationships are well-defined for traditional rigid, metallic marine rotors, there exist very few works that discuss the hydroelastic scaling of self-adaptive marine rotors. Hydroelastic scaling of wave-structure interaction problems can be found in [24, 25], though these and other similar wave-structure interaction studies typically do not involve rotating components or cavitation. Hydroelastic scaling of surface-piercing propellers have been discussed in [26], but it applies only to isotropic metallic blades. Young [27] derived and validated dynamic hydroelastic similarity relations for self-adaptive composite rotors and demonstrated the importance of material scaling to ensure similar load-deformation characteristics. For flexible composite rotors, scaling of the material is highly nontrivial, especially for the prediction of material failure. The effects of specimen size, material properties, stacking sequence, number of plies, and fiber orientation, among other characteristics, have been shown to have a significant effect on the failure strength of composite materials, as well as the failure mode of the test specimen [28–33].

The objectives of this work are to investigate and to verify the dynamic similarity relationships for the hydroelastic response and potential failure mechanisms of self-adaptive composite marine propellers. A fully coupled, three-dimensional (3D), boundary element method-finite element method (BEM-FEM) is used to compare the model and full-scale responses of a self-adaptive composite propeller designed for a naval combatant. The 3D BEM-FEM solver is summarized in Section 2, the propeller characteristics are described in Section 3, the scaling results are shown in Section 4, and the major findings are reported in Section 5.

2. Methodology

2.1. 3-D BEM-FEM Solver

A 3D coupled BEM-FEM solver is applied herein for the analysis of propeller performance. The solver considers the effects of nonlinear geometric coupling due to thickness and 3D effects, spatially varying inflow, transient sheet cavitation, load dependent FSI response, and potential strength and stability issues. The fluid behavior is assumed to be governed by the incompressible Euler equations in a blade-fixed rotating coordinate system as follows:

where is the total velocity, is the physical time, is the hydrodynamic pressure, is the water density, is the gravitational acceleration, is the propeller rotational speed vector, and is the non-inertial blade-fixed coordinates vector that rotates with the reference blade.

The total velocity () is expressed as the summation of the inflow velocity () and the perturbation potential velocity ():

where . The effective wake velocity contains both the nominal inflow velocity (i.e., in the absence of the propeller) and the vortical interactions between the propeller and the inflow [35]. This allows the perturbation flow field to be treated as incompressible, inviscid, and irrotational and it is governed by the Laplace equation: . Readers should refer to [14, 15, 36–41] for more details.

To consider FSI effects, the perturbation potential is further decomposed into components due to rigid blade rotation, , and elastic blade deformation, , which are solved using the 3D BEM with proper kinematic and dynamic boundary conditions. Similarly, the total hydrodynamic pressure is decomposed into components due to rigid blade rotation, , and elastic blade deformation, . Integration of over the wetted blade area can be expressed as the added mass matrix () times the nodal acceleration vector () and the added damping matrix () times the nodal velocity vector (), resulting in the equation of motion defined with respect to the rotating blade-fixed coordinate system in the time domain:

where is the structural nodal displacement vector; , , and are the structural mass, damping, and stiffness matrices, respectively; , , and are the centrifugal force, the Coriolis force, and the hydrodynamic force (due to rigid blade rotation) vectors, respectively. Detailed formulation of these vectors and matrices can be found in [14, 15].

The dynamic equation of motion is solved using the commercial FEM solver, ABAQUS/Standard [42], where , , and are obtained from the BEM solver. The direct cyclic algorithm in ABAQUS/Standard, which combines a modified Newton method with a Fourier representation of the solution and residual vectors, is used to calculate the dynamic blade response in unsteady flows and the effects of noncacheable large blade deformations are considered by iterating between the fluid and solid solvers until the solution converges, which generally occurs within 6-7 iterations. The propeller blades are discretized using 3D, reduced integration quadratic continuum solid elements. The composite material is modeled using orthotropic material properties for each of the elements along with a primary axis to represent the orientation angle of the fibers. Material failure initiation is modeled using the Hashin failure initiation criteria [43]. Further details of the formulation, including numerical implementation validation studies, can be found in [13–15, 39–41, 44–46].

3. Propeller Characteristics

In previous works [20–22], the authors designed and analyzed a pair of self-adaptive composite propellers, shown in Figure 1 and modeled after the classic propeller, DTMB 4383, details of which can be found in [47]. The propeller was assumed to be made of carbon fiber reinforced polymer (CFRP) and has a diameter of 5.18 m. The inflow wake in the propeller plane, shown in Figure 1, is based on data presented in [34] and is asymmetric because of the upstream hull form of the twin-shafted naval combatant craft.

(a)

(b)

Figure 1: (a) Undeformed propeller geometry for the modified DTMB 4383 made of CFRP. (b) Assumed inflow wake at the propeller plane behind a twin-shafted naval combatant [34].

4. Comparison of the Model and Full-Scale Responses and Potential Failure Mechanisms of a Self-Adaptive Composite Marine Propeller

A self-adaptive composite propeller is designed to depitch under normal forward loading [21]. For a self-adaptive propeller, the blade deformation and the resulting hydrodynamic performance depend on the total dimensional load corresponding to a specific advance speed, , and rotational frequency, . As increases, the depitching action increases, which in turn requires a greater increase in in order to meet the vessel thrust requirement. In order to validate the predicted propeller performance and to investigate potential failure mechanisms of self-adaptive composite propellers, proper scaling relationships must be developed that capture the load-deformation response of the propeller blades.

4.1. Dynamic Hydroelastic Similarity Relations

In a previous work by Young [27], hydroelastic similarity conditions and scaling factors for self-adaptive composite marine propellers were presented in detail and are summarized here for clarity. Scaling factors are defined as the ratio of model-scale to full-scale parameters and are denoted by with appropriate subscripts representing the parameter of interest. Assuming geometric similarity, the characteristic length scale is , where is the diameter of the model-scale propeller and is the diameter of the full-scale propeller. For a typical towing tank or cavitation tunnel test facility with water as the fluid medium, the scaling factors for gravity (), fluid density (), viscosity (), and speed of sound () are approximately 1. To simulate the same operating conditions between the model and the prototype, the scaling ratios for the advance coefficient (, where ) and cavitation number (, where , is the absolute hydrostatic pressure at the propeller shaft axis, and is the saturated vapor pressure) must also be 1, which can be achieved by controlling the inflow velocity, propeller rotational frequency, and pressure inside a cavitation tunnel. Additionally, to achieve similar FSI response between the model and full-scale propellers, the scaling ratios for solid density () and Poisson’s ratio () must equal 1 and the effective structural stiffness ratios must be scaled as follows:

where is the Young’s modulus in the -direction and is the shear modulus in the -orientation. Assuming that the above similarity conditions are met, only three critical nondimensional parameters remain: Reynolds number (), Froude number (), and Mach number (), where the parameters and are respectively the fluid kinematic viscosity and speed of sound. Reynolds number, Froude number, and Mach number similarity cannot be satisfied simultaneously. Hence, the selection of the appropriate scaling depends on the objective of the model-scale study.

Reynolds number similarity should be applied when viscous effects and the influence of large-scale vortices are critical, but it is difficult to achieve in typical cavitation tunnels because it requires the model-scale velocity to be faster than the full-scale velocity, or . For typical forward operating conditions, the Reynolds number should be large enough such that viscous forces should be small compared to inertial and gravitational forces, and hence, Reynolds number similarity is not considered herein.

Froude number similarity is critical for flow conditions where gravitational forces are important (i.e., at small Froude numbers). There are four principal forces that must be considered for a flexible composite propeller: solid elastic restoring force, gravitational force, hydrodynamic inertial force, and rotor inertial force. When gravitational forces are significant, for example, at lower speeds, Froude scaling is the only model that maintains the same ratios between all four of the dominant forces [27]. Further, Froude scaling allows smaller , which is easier to achieve in cavitation tunnel and towing tank studies.

Mach number similarity is typically required when flow compressibility is an issue, which should not be a concern for marine propellers at both the model and full-scale because the Mach number is typically significantly less than 1. However, Mach scaling allows the same material and layering scheme to be used between the model and full-scale propellers [27]. In order to absorb the high hydrodynamic loads, self-adaptive composite marine propellers typically require complex 3D geometry with small aspect ratio and solid material layup. Hence, it is extremely difficult, if not impossible, to find the proper fiber and matrix properties and layering scheme, such that the 3D distribution of the structural density, bending rigidity, torsional rigidity, and coupled bending-torsional rigidity are the same between the model and the full-scale propellers unless the same material and layering scheme are used.

The objectives of this work are to verify the dynamic similarity relationships presented in [27] by comparing the model-scale and full-scale responses as well as potential failure mechanisms of the self-adaptive composite marine propeller shown in Figure 1. For convenience, the relevant scaling parameters for Froude and Mach scaling as derived in [27] are summarized in Table 1.

Table 1: Scaling ratios for critical parameters assuming geometric similarity and the following conditions to be true [27]: .

Parameter
Expression
Froude scale
Mach scale

Froude number

1

Reynolds number

Mach number

1

Axial velocity

1

Angular velocity

Elastic modulus

1

Shear modulus

1

Poisson's ratio

1
1

Elastic force

Gravitational force

Hydrodynamic inertial force

Rotor inertial force

Natural frequencies

4.2. Model and Full-Scale Parameters

The validity of the Mach-scale and Froude-scale similarity relationships shown in Table 1 is demonstrated using the 3D BEM-FEM solver described in Section 2. The 5.18 m (17 ft) full-scale propeller is described in Section 3, and shown in Figure 1. The model-scale propeller is assumed to be a geometrically similar 1/17-scale model with a diameter of m (1 ft). The relevant parameters for the full-scale, Mach-scale, and Froude-scale composite propellers are shown in Table 2. Note that the longitudinal (), transverse (), and shear () strength components are included. It should be noted that both the Mach-scale and Froude-scale studies are assumed to be conducted in cavitation tunnels such that the cavitation number, , can be made to be the same as the full-scale propeller by controlling tunnel pressure. Additionally, no values are listed for the Froude-scale material strength parameters in the last six rows Table 2 because they cannot be derived theoretically since, as noted in the previous section, it is difficult to find the fiber and matrix combinations that will satisfy all the effective structural density and moduli scaling requirements.

Table 2: Geometric and material parameters of the full-scale, Mach-scale, and Froude-scale composite propellers.

Parameter
Full scale
Mach scale
Froude scale

1
1/17
1/17

D (m)
5.18161
0.3048
0.3048

(knots)
0–20
0–20
0–4.8

(rpm)
0–110
0–1680
0–400

(kg/m³)
1025
1025
1025

(kg/m³)
2150
2150
2150

(GPa)
80.0
80.0
4.71

(GPa)
10.0
10.0
0.59

(GPa)
3.30
3.30
0.19

(GPa)
3.90
3.90
0.23

0.32
0.32
0.32

0.45
0.45
0.45

(MPa)
1950
1950
—

(MPa)
1480
1480
—

(MPa)
48
48
—

(MPa)
200
200
—

(MPa)
79
79
—

(MPa)
50
50
—

4.3. Steady-State Response in Uniform Inflow

A comparison of the load-dependent deformation responses and resulting performance of the full-scale and model-scale propellers under steady, uniform inflow conditions is shown in Figure 2. By applying the scaling relationships shown in Table 1, the variations of the normalized deformation responses (change in pitch angle, , change in skew, , and normalized change in rake, ) and hydrodynamic load coefficients (thrust coefficient, , torque coefficient, , and efficiency, ) for the deformed geometry, with for both the Froude-scale and Mach-scale propellers agree well with the full-scale propeller. Note that while is the same between the Mach-scale and the full-scale propellers, the advance speed for the Froude-scale propeller should be reduced by . For ease of comparison, the full-scale, Mach-scale, and Froude-scale results are shown in the same graph as a function of the full-scale .

(a)

(b)

Figure 2: Comparison of propeller performance curves (a) and deformation responses (b) of the full-scale, Froude-scale, and Mach-scale propellers as a function of the full-scale advance speed .

Figure 3 shows the first five wetted resonant frequencies of the full-scale, Mach-scale, and Froude-scale propellers normalized by the appropriate propeller rotational frequency corresponding to a full-scale advance speed of knots. As expected, because , the normalized frequencies are the same.

Figure 3: First five wetted excitation frequencies, , normalized by the propeller rotational frequency, , for the full-scale, Froude-scale, and Mach-scale propellers. The results correspond to the full-scale propeller operating at knots.

The results demonstrate that by following the similarity relations shown in Table 2, both the Froude-scale and Mach-scale composite propellers are able to correctly predict the load-dependent deformation response, hydrodynamic performance, and susceptibility to transient and/or resonant vibration of the full-scale composite propeller.

4.4. Transient Response in Spatially Varying Wake Flow

To further verify the hydroelastic similarity relations, results are shown in this section for the full-scale, Mach-scale, and Froude-scale composite propellers operating in the spatially varying wake shown in Figure 1.

Comparisons of the circumferentially averaged values of the hydrodynamic coefficients and blade deformations corresponding to full-scale advance speeds of knots and knots are shown in Table 3. For both speeds, the Reynolds numbers are high enough for the full-scale and Mach-scale propellers such that viscous forces should be negligible compared to inertial forces. It should be noted that transition may occur on the Froude-scale propeller, and hence, special treatment may be needed at the blade leading edge to ensure fully turbulent flow. However, since the 3D BEM-FEM model assumes inviscid flow, viscous effects will not be discussed herein. Although the Froude number is rather low at 10 knots, the mean axial force coefficient, , and axial moment coefficient, , where and are the axial force and moment, respectively, as well as the changes in blade tip pitch angle, , and the normalized blade tip deflections, for the full-scale, Mach-scale, and Froude-scale propellers are in good agreement with each other.

Table 3: Circumferentially averaged values of the hydrodynamic load coefficients and blade deformations in wake fow for the full-scale, Mach-scale, and Froude-scale propellers ( kg/m³,  m²/s,  m/s²,  m/s).

Full scale
Mach scale
Froude scale

(knots)
10
20
10
20
2.425
4.851

Fr
0.340
1.435
5.780
24.406
0.340
1.435

21.45
44.05
1.26
2.59
0.305
0.628

6.41
12.8
6.41
12.8
1.55
3.11

0.0440
0.0409
0.0440
0.0408
0.0439
0.0403

0.0118
0.0106
0.0117
0.0106
0.0118
0.0105

(deg)
0.569
2.092
0.569
2.075
0.574
2.112

/D
0.00363
0.0135
0.00363
0.135
0.00366
0.0135

Comparisons of the time histories of the hydrodynamic responses, deformations, and cavitation volumes for the full-scale, Mach-scale and Froude-scale propellers are shown in Figures 4, 5 and 6. Good agreement is observed between the full-scale and Froude-scale propellers because Froude scaling has the benefit of preserving the ratios between the four dominant forces, as noted in Table 1. However, as shown in Figures 4, 5 and 6, some discrepancies could be observed between the full-scale and Mach-scale propellers, particularly at the lower full-scale speed of knots. As shown in Table 3, the Froude number is higher for the Mach-scale propeller, which will lead to under-prediction of the gravitational force relative to the other forces, which is consistent with the relations shown in Table 1. Consequently, the time histories of the per-blade axial force and moment coefficients of the Mach-scale propeller are slightly different from the full-scale propeller, although the mean values are approximately the same (see Table 3). Moreover, the Mach-scale propeller tends to under-predict the cavitation volumes, as shown in Figure 6. As the speed increases, the relative importance of gravitational forces decreases, which leads to better agreement between the Mach-scale and full-scale propellers at knots.

(a)

(b)

Figure 4: Unsteady axial force coefficient, , and axial moment coefficient , of the reference blade as a function of blade angle at a full-scale advance speed of knots (a) and knots (b) for the full-scale, Froude-scale, and Mach-scale propellers.

(a)

(b)

Figure 5: Unsteady changes in blade tip pitch angle (a), , and normalized tip deflection (b), , as a function of blade angle at a full-scale advance speed of knots and knots for the full-scale, Froude-scale, and Mach-scale propellers.

(a)

(b)

Figure 6: Unsteady normalized backside cavitation volume as a function of blade angle at a full-scale advance speed of knots (a) and knots (b) for the full-scale, Froude-scale, and Mach-scale propellers.

To better illustrate the difference between the Mach-scale and full-scale propellers at knots, comparisons of the pressure coefficient contours are shown in Figure 7. The pressure coefficient is defined as , where is the absolute total pressure and is the absolute hydrostatic pressure at the propeller shaft axis. The darker regions indicate where cavitation develops, that is, . As shown in Figure 7, the pressure distribution and cavitation coverage between the Mach-scale and full-scale propellers are very similar. The results demonstrate that although Froude number effects do influence the spatial variations of the dynamic blade loads and cavitation volumes, the effects are very limited at high speeds, where cavitation, material and/or stability failure are potential concerns.

(a)

(b)

Figure 7: Unsteady backside pressure coefficient, , contours at knots for the full-scale (a) and Mach-scale (b) propellers.

Depending on the length scale, , constructing a geometrically similar model-scale propeller following Froude similarity with the required 3D distribution of the structural density, as well as bending, torsional, and bending-torsional rigidity can be very difficult, especially for a self-adaptive composite propeller made of anisotropic laminates stacked in complex, 3D configurations. As shown in Table 2, the difference in the bending and shear moduli can be significant. It is difficult to ensure that all of the moduli ( and ) and Poisson’s ratios () scale according to the values given in Table 2 while keeping the effective solid density () the same between the model and the full-scale propellers. Changing even one of the material parameters can have a significant effect on the 3D load-deformation characteristics. Additionally, assessing structural integrity of the blades is not feasible for Froude similarity because of the need to scale the structural strength parameters.

As shown in Table 2, the same material and layering scheme could be used for the Mach-scale propeller, which significantly simplifies the material scaling. To accommodate the reduced size of the Mach-scale model, ply-level scaling could be used to scale the composite layup using the same type of laminates to ensure similar elastic response of the composite blade [29]. However, the failure strengths may differ between the model and the full-scale propellers even with Mach-scaling of the geometry and operating conditions, and ply-level scaling of the composite with the same material properties, because of size effects attributed to material uncertainty with size [28–30]. Nevertheless, estimates of the susceptibility for first-ply or initial material failure can be achieved. The dominant failure modes for a composite rotor blade constructed of solid multilayered laminates in flexure are matrix tensile failure and delamination [22]. While there are many different models for the prediction of composite failure initiation, the commonly used Hashin failure initiation models [43] are applied herein and are defined as,

where is the normal stress in the -direction, is the shear stress in the -direction, and , , , , and are material strength parameters as defined in Table 2. Initial material failure is assumed to occur when or .

Comparisons of the matrix tensile and delamination failure initiation indicators at the blade root between the Mach-scale and full-scale propellers operating in steady, uniform inflow at knots are shown in Figure 8. It should be noted that the trends shown for steady, uniform inflow are similar to those for the propeller operating in a spatially varying wake. For this propeller, failure initiation occurs at the blade root in the trailing edge region because of the high skew and the assumed fixed boundary at the root, which tends to overestimate the stress concentrations. Nevertheless, the failure indicator contours are very similar between the Mach-scale and full-scale propellers when operating in both steady and unsteady flow conditions. The results demonstrate the critical advantage of Mach scaling—the ability to preserve similar 3D stress distribution by allowing the same material and layering scheme to be used, and hence allow investigation of potential failure mechanisms of the full-scale self-adaptive composite material propeller when conducting model scale testing.

Figure 8: Comparisons of the matrix tensile and delamination composite failure initiation indicators at the blade root for the full-scale (top) and Mach-scale (bottom) propellers operating in steady, uniform inflow at knots.

5. Conclusions

A previously validated 3D BEM-FEM solver is used to compare the model (1/17-scale) and full-scale hydroelastic responses and potential failure mechanisms of a self-adaptive composite propeller designed for a naval combatant. The effects of spatially varying inflow, transient sheet cavitation, and load-dependent blade deformations are considered.

The critical scaling ratios are shown in Table 1. The results show that Froude scaling has the benefit of being able to maintain the ratios of the solid and fluid inertial forces, gravitational force, and elastic restoring force. However, it will be very difficult to properly scale the solid density and all the elastic material properties of a Froude-scale self-adaptive composite propeller, and will be nearly impossible to properly scale the material failure strengths. Mach scaling, on the other hand, has the benefit of allowing the same material and layering scheme to be used between the model and full-scale propellers, which helps to preserve the 3D stress distributions and potential failure mechanisms. It should be emphasized that flow compressibility effects are typically negligible for marine propellers and are ignored in the current analysis. Mach number similarity simply implies that the relative inflow velocity should be the same between the model and full-scale propellers, which can be achieved in cavitation tunnel studies. The results show that Mach scaling will under-predict the gravitational force compared to the other three dominant forces because of the higher Froude number at model scale, particularly in the lower-speed range. Nevertheless, the results show that both the Mach-scale and Froude-scale propellers are able to reproduce the average hydrodynamic load coefficients, load-dependent deformations, and susceptibility to resonant vibrations of the full-scale self-adaptive composite propeller. However, some differences between the Mach-scale and full-scale propellers could be observed in the spatial variation of the dynamic blade loads and cavitation volumes when operating in a spatially varying wake, particularly at the low speed range. At the high speed range where cavitation, material and/or stability failure are potential concerns, very limited differences are observed between the dynamic response of the Mach-scale and full-scale propellers. Moreover, the results demonstrate that the Mach-scale propeller is able to emulate the 3D distribution of the material failure initiation indicators, which is critical to assessing the structural integrity and safe operating envelope of the full-scale self-adaptive composite marine propeller.

It should be noted that the current work assumes ply-level scaling of the composite layup, and that the material properties and manufacturing processes are the same between the Mach-scale and full-scale composite propellers. However, it is well known that the failure strengths of CFRPs typically decrease with increasing size due to the increasing material and geometry uncertainties, for example, misalignment or kinking of fibers, existence of voids, uneven distribution of fiber or matrix volumes, and so forth. Moreover, different manufacturing techniques may be required between the model and the full-scale because of challenges with the solid material layup and complex 3D geometry of marine propellers. Hence, additional research is needed to address composite scaling issues, particularly related to the failure strengths, influence of residual stresses, load-sequence effects, and fatigue strengths when subject to long term salt water immersion and potential large temperature variations.

The analyses shown in this paper are limited to inviscid flow assumptions, and hence, Reynolds effects are not considered. Although viscous forces should be negligible compared to inertial forces for most speeds of interest in normal forward operating modes, transition and viscous effects on the blade tip may be of concern at the model-scale because of the reduced Reynolds number. Moreover, viscous effects may dominate for extreme off-design conditions such as crashback, where the flow is dominated by large-scale flow separations and transient ring vortices. Hence, additional research is also needed to investigate viscous effects on the dynamic hydroelastic response and potential failure mechanisms of self-adaptive composite marine propellers.

Acknowledgments

The authors are grateful to the Office of Naval Research (ONR) and Dr. Ki-Han Kim (program manager) through Grant nos. N00014-09-1-1204 and N00014-10-1-0170 for their financial support. This work was also supported in part by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MEST) through the GCRC-SOP Grant no. 2012-0004783.

References

R. Ganguli and I. Chopra, “Aeroelastic optimization of a composite helicopter rotor,” in Proceedings of the 4th AIAA Symposium on Multidisciplinary Analysis and Optimization, pp. 21–23, Cleveland, Ohio, USA, 1992.

R. Ganguli and I. Chopra, “Aeroelastic tailoring of composite couplings and blade geometry of a helicopter rotor using optimization methods,” Journal of the American Helicopter Society, vol. 42, no. 3, pp. 218–228, 1997.

B. Glaz, P. P. Friedmann, and L. Liu, “Helicopter vibration reduction throughout the entire flight envelope using surrogate-based optimization,” Journal of the American Helicopter Society, vol. 54, no. 1, Article ID 012007, 2009.

O. Soykasap and D. H. Hodges, “Performance enhancement of a composite tilt-rotor using aeroelastic tailoring,” Journal of Aircraft, vol. 37, no. 5, pp. 850–858, 2000.

A. T. Lee and R. G. J. Flay, “Compliant blades for passive power control of wind turbines,” Wind Engineering, vol. 24, no. 1, pp. 3–11, 2000.

D. W. Lobitz and P. S. Veers, “Load mitigation with bending/twist-coupled blades on rotors using modern control strategies,” Wind Energy, vol. 6, no. 2, pp. 105–117, 2003.

A. Maheri and A. T. Isikveren, “Performance prediction of wind turbines utilizing passive smart blades: approaches and evaluation,” Wind Energy, vol. 13, no. 2-3, pp. 255–265, 2010.

S. Gowing, P. Coffin, and C. Dai, “Hydrofoil cavitation improvements with elastically coupled composite materials,” in Proceedings of the 25th American Towing Tank Conference, Iowa City, Iowa, USA, 1998.

B. Y. H. Chen, S. K. Neely, T. J. Michael et al., “Design, fabrication and testing of pitch-adapting (Flexible) composite propellers,” in Proceedings of the Society of Naval Architects and Marine Engineers Propellers/Shafting 2006 Symposium, Virginia Beach, Va, USA, September 2006.

C. C. Lin, Y. J. Lee, and C. S. Hung, “Optimization and experiment of composite marine propellers,” Composite Structures, vol. 89, no. 2, pp. 206–215, 2009.

Y. J. Lee and C. C. Lin, “Optimized design of composite propeller,” Mechanics of Advanced Materials and Structures, vol. 11, no. 1, pp. 17–30, 2004.

C. C. Lin and Y. J. Lee, “Stacking sequence optimization of laminated composite structures using genetic algorithm with local improvement,” Composite Structures, vol. 63, no. 3-4, pp. 339–345, 2004.

Y. L. Young, T. J. Michael, M. Seaver, and S. T. Trickey, “Numerical and experimental investigations of composite marine propellers,” in Proceedings of the 26th Symposium on Naval Hydrodynamics, Rome, Italy, September 2006.

Y. L. Young, “Time-dependent hydroelastic analysis of cavitating propulsors,” Journal of Fluids and Structures, vol. 23, no. 2, pp. 269–295, 2007.

Y. L. Young, “Fluid-structure interaction analysis of flexible composite marine propellers,” Journal of Fluids and Structures, vol. 24, no. 6, pp. 799–818, 2008.

Z. Liu and Y. L. Young, “Utilization of bend-twist coupling for performance enhancement of composite marine propellers,” Journal of Fluids and Structures, vol. 25, no. 6, pp. 1102–1116, 2009.

Z. Liu and Y. L. Young, “Static divergence of self-twisting composite rotors,” Journal of Fluids and Structures, vol. 26, no. 5, pp. 841–847, 2010.

M. R. Motley, Z. Liu, and Y. L. Young, “Utilizing fluid-structure interactions to improve energy efficiency of composite marine propellers in spatially varying wake,” Composite Structures, vol. 90, no. 3, pp. 304–313, 2009.

M. Motley and Y. Young, “Reliability-based global design of selfadaptive marine rotors,” in Proceedings of the 3rd ASME Joint USEuropean Fluids Engineering Summer Meeting and 8th International Conference on Nanochannels, Microchannels, and Minichannels, Montreal, Canada, 2010.

M. Motley and Y. Young, “Performance-based design of adaptive composite marine propellers,” in Proceedings of the 28th Symposium on Naval Hydrodynamics, Pasadena, Calif, USA, 2010.

M. R. Motley and Y. L. Young, “Performance-based design and analysis of flexible composite propulsors,” Journal of Fluids and Structures, vol. 27, no. 8, pp. 1310–1325, 2011.

M. R. Motley and Y. L. Young, “Influence of design tolerance on the hydrodynamic response of self-adaptive marine rotors,” Composite Structures, vol. 94, no. 1, pp. 114–120, 2011.

M. R. Motley, M. Nelson, and Y. L. Young, “Integrated probabilistic design of marine propulsors to minimize lifetime fuel consumption,” Ocean Engineering, vol. 45, pp. 1–8, 2012.

C. Shumin, A. S. J. Swamidas, and J. J. Sharp, “Similarity method for modeling hydroelastic offshore platforms,” Ocean Engineering, vol. 23, no. 7, pp. 575–595, 1996.

O. Faltinsen, Hydrodynamics of High-Speed Marine Vehicles, Cambridge University Press, New York, NY, USA, 2005.

N. Olofsson, Force and flow characteristics of a partially submerged propeller [Ph.D. thesis], Chalmers University of Technology, Goteborg, Sweden, 1996.

Y. L. Young, “Dynamic hydroelastic scaling of self-adaptive composite marine rotors,” Composite Structures, vol. 92, no. 1, pp. 97–106, 2010.

S. Kellas and J. Morton, “Strength scaling in fiber composites,” AIAA Journal, vol. 30, no. 4, pp. 1074–1080, 1992.

K. E. Jackson, “Scaling effects in the flexural response and failure of composite beams,” AIAA Journal, vol. 30, no. 8, pp. 2099–2105, 1992.

K. E. Jackson, S. Kellas, and J. Morton, “Scale effects in the response and failure of fiber reinforced composite laminates loaded in tension and in flexure,” Journal of Composite Materials, vol. 26, no. 18, pp. 2674–2705, 1992.

S. R. Swanson, “Scaling of impact damage in fiber composites from laboratory specimens to structures,” Composite Structures, vol. 25, no. 1–4, pp. 249–255, 1993.

A. P. Christoforou and A. S. Yigit, “Scaling of low-velocity impact response in composite structures,” Composite Structures, vol. 91, no. 3, pp. 358–365, 2009.

A. Rehan and K. H. Grote, “Composite strength scaling effect using progressive degradation model,” in Proceedings of the WRI World Congress on Computer Science and Information Engineering (CSIE'09), pp. 1–5, Los Angeles, Calif, USA, April 2009.

M. Hugel, An evaluation of propulsors for several navy ships [M.S. thesis], Massachusetts Institute of Technology, Cambridge, Mass, USA, 1992.

J. K. Choi and S. A. Kinnas, “Prediction of non-axisymmetric effective wake by a three-dimensional euler solver,” Journal of Ship Research, vol. 45, no. 1, pp. 13–33, 2001.

J. E. Kerwin, S. A. Kinnas, J. T. Lee, and W. Z. Shih, “A surface panel method for the hydrodynamic analysis of ducted propellers,” Transactions of Society of Naval Architects and Marine Engineers, vol. 95, 1987.

S. A. Kinnas and N. E. Fine, “A numerical nonlinear analysis of the flow around two-and three-dimensional partially cavitating hydrofoils,” Journal of Fluid Mechanics, vol. 254, pp. 151–181, 1993.

Y. Young and S. Kinnas, “A BEM for the prediction of unsteady midchord face and/or back propeller cavitation,” Journal of Fluids Engineering, vol. 123, pp. 311–319, 2001.

Y. L. Young and S. A. Kinnas, “Application of BEM in the modeling of supercavitating and surface-piercing propeller flows,” Journal of Computational Mechanics, vol. 32, no. 4–6, pp. 269–280, 2003.

Y. L. Young and S. A. Kinnas, “Performance prediction of surface-piercing propellers,” Journal of Ship Research, vol. 48, no. 4, pp. 288–304, 2004.

Y. L. Young and Y. T. Shen, “A numerical tool for the design/analysis of dual-cavitating propellers,” Journal of Fluids Engineering, vol. 129, no. 6, pp. 720–730, 2007.

ABAQUS, ABAQUS Version 6.5 Documentation, ABAQUS, Inc., 1080 Main Street, Pawtucket, RI 02860, 2005

Z. Hashin, “Failure criteria for unidirectional composites,” Journal of Applied Mechanics, vol. 47, no. 2, pp. 329–334, 1980.

Y. L. Young, “Hydroelastic response of composite marine propellers,” in Proceedings of the Society of Naval Architects and Marine Engineers Propellers/Shafting 2006 Symposium, Williamsburg, Va, USA, September 2006.

Y. L. Young, Z. Liu, and M. Motley, “Influence of material anisotropy on the hydroelastic behaviors of composite marine propellers,” in Proceedings of the 27th Symposium on Naval Hydrodynamics, Seoul, Korea, 2008.

Y. L. Young and B. R. Savander, “Design, analysis, and challenges of largescale surface-piercing propellers,” Ocean Engineering, vol. 38, no. 13, pp. 1368–1381, 2011.

R. Boswell, “Design, cavitation performance and open-water performance of a series of research skewed propellers,” Tech. Rep. 3339, DTNSRDC, 1971.