Reliability Assessment of Water Hydraulic-Drive Wave-Energy Converters

: The main objective of the current paper is to determine fatigue critical details of wave-energy converters (WECs) and to evaluate their reliability. For this purpose, a new model of a water-hydraulic drive WEC is designed based on the IPS (company Interproject Service) buoy working principle, and the hydrodynamic frequency domain analysis method is carried out accordingly. A reliability assessment method, which combines the stress-number (SN) curves approach and linear damage cumulative theory, has been developed speciﬁcally for the proposed WEC. In order to obtain the stochastic wave load spectrum, 22 years of wave data of the intended deployment site (South China Sea), covering the span of 1988–2009, were analyzed, discussed and processed. A ﬁnite element analysis model is also established to calculate the maximum Von Mises stress of the fatigue critical component. The reliability indices and conﬁguration parameters of the water hydraulic-drive WEC are obtained by applying ﬁrst order reliability method, considering the uncertainty of the model parameters and the fatigue limit state. Optimization of conﬁguration parameters is also developed based on the reliability method for the sake of approaching the expected service life. The results show that the fatigue life of the water turbine can reach the design service life of 20 years as the conﬁguration parameters meet the speciﬁc condition. The cumulative reliability index and annual reliability index are recommended as 2.1 and 3.5, respectively.


Introduction
With the upgrading of energy consumption globally and rising industrialization, renewable energy development has become one of the hottest issues in the world [1]. Many countries and regions have drawn up sustainable development plans and have committed to the development of marine renewable energy resources [2]. Wave energy is regarded as having the greatest potential for development [3,4], and the technologies harvesting energy from waves have received more consideration in the last decade [5]. Ocean wave energy is clean and free from pollution, and the world's largest estimated reserves reach 80,000 TW h/year [6]. Wave energy has the characteristics of good predictability and high energy density, which makes wave energy utilization favored on coastal regions.
Wave energy is produced by wind action and has only a small energy loss even if it travels thousands of kilometers. The magnitude of wave energy is often represented in terms of power per unit wave-front (kW/m), and the main parameters representing wave-energy potential are the significant wave height and average energy period [7]. The evaluation of wave-energy resources is the premise for the exploitation and utilization and the design of relevant devices [8]. Research on out accordingly. The purpose of this work is to perform a reliability assessment method for a water hydraulic-drive WEC and calibrate a preliminary scheme and system reliability index.
The paper is organized in six sections: Section 1 introduces the background, methods and research purpose of the problem. The structural design and hydrodynamic analysis of the water hydraulic drive WEC are illustrated in Section 2. Section 3 presents the reliability model of fatigue failure. Wave characteristics of the target sea area and fatigue-critical details of the present WEC are carried out in Section 4. The results are developed and discussed in Section 5, followed by conclusions in Section 6.

Structural Design
The proposed water hydraulic drive WEC mainly includes a floater and a vertical tube suspended in sea water, which is similar to the IPS buoy [29], oscillating up and down under the action of waves. The tube is rigidly connected to the floater through the filter bracket which is formed by welding the support structure arranged in intervals. An axial flow turbine, as a PTO element, is located in the middle of the tube and moves with the floater at the same amplitude. When the floater oscillates in heave, the water contained inside the tube drives the turbine under the action of inertia force. The energy is obtained from ocean waves when the water in the tube moves relative to the floater. The basic working principle and structural parameters are represented in Figure 1. The floater consists of a cylindrical buoy and an inverted cone, both of radius . The submerged depth of the cylindrical part equal to radius below the surface of calm water. The base angle of the conical part of the floater equal to . The tube, with radius and length , is submerged vertically in seawater of a certain depth (relatively large) and open at both ends. Considering the barrier of the tube wall, the water column contained in the tube can be considered independent of the direct action of wave load.
The water column contained in the tube will be resisted by external seawater when it flows out or into the tube, which can be equivalent to the added mass of the water column. The added mass of the water in the tube, assuming the radius is and the thickness of the two ends of the wall is negligible, is [30]. Here, = 0.6133 is an approximated distance of the water column that can be considered as independent of wave frequency . The floater consists of a cylindrical buoy and an inverted cone, both of radius R. The submerged depth of the cylindrical part equal to radius R below the surface of calm water. The base angle of the conical part of the floater equal to γ. The tube, with radius r and length L, is submerged vertically in seawater of a certain depth (relatively large) and open at both ends. Considering the barrier of the tube wall, the water column contained in the tube can be considered independent of the direct action of wave load.
The water column contained in the tube will be resisted by external seawater when it flows out or into the tube, which can be equivalent to the added mass of the water column. The added mass of the water in the tube, assuming the radius is r and the thickness of the two ends of the wall is negligible, is ρlπr 2 [30]. Here, l = 0.6133r is an approximated distance of the water column that can be considered as independent of wave frequency ω. The turbine, which is capable of self-rectifying, has a rotor with fan blades and the rotor blades are installed in the hub through a pivot mounted near the blade root; hence, they are able to swing back and forth within two set angles, which are restricted by pre-installed mechanical stops. Both action planes of the mechanical stops located in the rotor hub come into contact with the inner wall of the hub when the rotor blades reach the limit position (see Figure 2). The oscillation range of the rotor blade is from the top dead center (TDC) to the bottom dead center (BDC) (see Figure 3). The rotor blades are mounted on a spherical hub to keep the hub clearance constant while the blade rotates. The outer casing is also spherical to ensure a stable tip clearance. The turbine, which is capable of self-rectifying, has a rotor with fan blades and the rotor blades are installed in the hub through a pivot mounted near the blade root; hence, they are able to swing back and forth within two set angles, which are restricted by pre-installed mechanical stops. Both action planes of the mechanical stops located in the rotor hub come into contact with the inner wall of the hub when the rotor blades reach the limit position (see Figure 2). The oscillation range of the rotor blade is from the top dead center (TDC) to the bottom dead center (BDC) (see Figure 3). The rotor blades are mounted on a spherical hub to keep the hub clearance constant while the blade rotates. The outer casing is also spherical to ensure a stable tip clearance.  The guide vanes sets are installed on each sides of the hub. Alternately, the rotor blades swing in hydrodynamic action of water column in the tube, and automatically select the correct positions for efficient operation whenever the water changes direction. When the water is moving in a positive direction, the rotor blades rotate from TDC to BDC under the action of hydrodynamic force. On the contrary, the rotor blades rotate from BDC to TDC when the water moves in the opposite direction, as shown in Figure 1. Therefore, the rotation direction of the rotor will not change when the water direction varies.
In present work, a simplified computational model is adopted, which assumes that the motion of the floater is restricted to the direction of heave and the effect of radiation from both ends of the tube on the wave field is negligible. If the turbine part is replaced by a compact piston, the water column contained in the tube should move with the same amplitude as the floater. Therefore, the force driving the movement of the water column in the tube should be the same as that of the turbine according to the law of conservation of energy.

Hydrodynamic Analysis
The heaving motions of the floater (obviously includes the turbine and tube) is denoted by ( ) the coordinate, with = 0 in the equilibrium position (calm water surface) and increasing upwards. Let ( ) be the oscillation displacement of the water column in the tube relative to the  The turbine, which is capable of self-rectifying, has a rotor with fan blades and the rotor blades are installed in the hub through a pivot mounted near the blade root; hence, they are able to swing back and forth within two set angles, which are restricted by pre-installed mechanical stops. Both action planes of the mechanical stops located in the rotor hub come into contact with the inner wall of the hub when the rotor blades reach the limit position (see Figure 2). The oscillation range of the rotor blade is from the top dead center (TDC) to the bottom dead center (BDC) (see Figure 3). The rotor blades are mounted on a spherical hub to keep the hub clearance constant while the blade rotates. The outer casing is also spherical to ensure a stable tip clearance.  The guide vanes sets are installed on each sides of the hub. Alternately, the rotor blades swing in hydrodynamic action of water column in the tube, and automatically select the correct positions for efficient operation whenever the water changes direction. When the water is moving in a positive direction, the rotor blades rotate from TDC to BDC under the action of hydrodynamic force. On the contrary, the rotor blades rotate from BDC to TDC when the water moves in the opposite direction, as shown in Figure 1. Therefore, the rotation direction of the rotor will not change when the water direction varies.
In present work, a simplified computational model is adopted, which assumes that the motion of the floater is restricted to the direction of heave and the effect of radiation from both ends of the tube on the wave field is negligible. If the turbine part is replaced by a compact piston, the water column contained in the tube should move with the same amplitude as the floater. Therefore, the force driving the movement of the water column in the tube should be the same as that of the turbine according to the law of conservation of energy.

Hydrodynamic Analysis
The heaving motions of the floater (obviously includes the turbine and tube) is denoted by ( ) the coordinate, with = 0 in the equilibrium position (calm water surface) and increasing upwards. Let ( ) be the oscillation displacement of the water column in the tube relative to the The guide vanes sets are installed on each sides of the hub. Alternately, the rotor blades swing in hydrodynamic action of water column in the tube, and automatically select the correct positions for efficient operation whenever the water changes direction. When the water is moving in a positive direction, the rotor blades rotate from TDC to BDC under the action of hydrodynamic force. On the contrary, the rotor blades rotate from BDC to TDC when the water moves in the opposite direction, as shown in Figure 1. Therefore, the rotation direction of the rotor will not change when the water direction varies.
In present work, a simplified computational model is adopted, which assumes that the motion of the floater is restricted to the direction of heave and the effect of radiation from both ends of the tube on the wave field is negligible. If the turbine part is replaced by a compact piston, the water column contained in the tube should move with the same amplitude as the floater. Therefore, the force driving the movement of the water column in the tube should be the same as that of the turbine according to the law of conservation of energy.

Hydrodynamic Analysis
The heaving motions of the floater (obviously includes the turbine and tube) is denoted by x(t) the coordinate, with x = 0 in the equilibrium position (calm water surface) and x increasing upwards. Let y(t) be the oscillation displacement of the water column in the tube relative to the floater. The energy conversion characteristics of the water turbine are mainly determined by the flow rate and rotation velocity, while the flow rate through the runner is related to the relative speed of the water column in the tube. In order to describe the performance of the turbine, a linearized PTO is assumed that can be represented as follows: where the constant B PTO is the equivalent PTO damping and the time-dependent power obtained by the turbine is P = f PTO .
y. Since the wave scale and the scale of floater motions is small relative to the wave length, the hydrodynamic analysis could be performed by using the linearized water wave theory. In fact, hydrodynamic analysis under linear wave conditions has been studied in other similar WEC devices [29,31,32]. The motion equations of the floater under the condition that the harmonic waves frequency is ω can be obtained as follows: ..
where m 1 and µ 1 are the mass and added mass of the floater, respectively; m 2 is the sum of mass and added mass of the water column in the tube, B is radiation damping for the floater, ρ is seawater density, 1025 kg/m 3 ; g is acceleration of gravity, 9.81 m/s 2 ; and S is the cross-sectional area of the cylinder tangent to the calm water surface, f e represents the action load of random waves, f PTO is the vertical load of the turbine that drive water column contained in the tube to oscillate relative to the tube. The hydrodynamic coefficients µ 1 and B and the wave action load f e are relative to the shape of the floater and the wave frequency ω. Considering the assumption of the linearized water wave theory, the WECs can be seen as a normal sinusoidal wave load with frequency ω acting on a linear system [33]. The time-dependent periodic wave excitation force, the motion of the floater and the relative motion of the water column in the tube can be described similarly in Equation (6): where F e , X and Y are complex amplitudes and F e (ω) = A w Γ(ω), A w is the incoming wave amplitude and Γ(ω) is the coefficient of action load. The magnitude of Γ(ω) can be involved radiation damping coefficient B(ω) by the Haskind relation [34]: By introducing the representation of complex amplitude, the following frequency-domain equations can be easily obtained: For incoming wave frequency ω and the magnitude of action load F e , the absolute value of the complex amplitudes X and Y can be yielded by solving the both of linearized frequency domain Equations (8) and (9): Through simple analysis, it is not difficult to obtain that the displacement of floater and the relative displacement of the water in the tube are only related to wave frequency ω and PTO damping B PTO . Based on the assumption of a linear time-harmonic system, the time average power obtained from the water turbine (PTO), P, is half of the maximum power, P max . The relationship between the two can be described as follows [35].

Torque Characteristics of the Power Take-Off (PTO)
Having obtained the maximum absorb power of PTO, the torque characteristics of the turbine blades can be further analyzed. As is known, the power obtained by the hydraulic turbine can be expressed by the mechanical torque T m and rotational angular velocity ω r .
Since the linear speed of wave motion is relatively low, the hydraulic turbine shows the characteristics of low rotational velocity and high torque after absorbing wave energy when acceleration devices are not present. One of the main purposes of this paper is to study the reliability of the water hydraulic drive WEC, and the hydraulic turbine is the torque transmission component of the WEC. Therefore, it is necessary to consider the maximum stress characteristics of the critical parts (e.g., blades) of the hydraulic turbine. Material fatigue is the most common fault mode of a hydraulic turbine [36]. Duo to the continuous propagation of the cracks, the turbine parts will gradually fail under the action of alternating and cyclic load. For the purpose of developing the stress characteristics of fatigue-critical components, the maximum power that the turbine can absorb during each wave period should be considered when reliability analysis is carried out.
It is assumed that the water hydraulic drive WEC is installed with λ 1 turbines, and each turbine carries λ 2 blades which bear an alternating maximum stress during a wave period. The maximum torque, T max , sustained by a single blade at the wave frequency ω can be expressed as follows: The position and magnitude of the maximum alternating stress of a single blade under each wave condition can be determined by finite element analysis. As the constant speed control strategy is relatively easy to implement in the robust operation of the turbine and the technology is relatively mature [37], the study in this paper is based on the assumption that the rotational velocity of the turbine remains constant. In order to achieve the specified service life, the maximum alternating stress in the critical parts can be adjusted through the number of turbines, λ 1 , the number of blades, λ 2 , as well as the PTO damping, B PTO .

Reliability Modelling of Fatigue Failure
The methods, of which the reliability model of fatigue failure can be performed, are the objective of this chapter. The SN-curves approach in combination with the PM-rules for considering sequence-independent linear damage accumulation, is an appropriate strategy when the condition monitoring is not present. SN curves indicate the relationship between the maximum cycles number and the failure stress of fatigue critical location. The SN-curves method has been widely used in ocean engineering design [38], which is more embodied in the reliability assessment of offshore wind turbines (OWTs) and WECs in recent years.

Fatigue Assessments
The parametric equations considering the bilinear SN curves can be expressed as follows [21]: where N represents the cycles number causing fatigue failure, S is the corresponding stress amplitude, K 1 , K 2 , m 1 as well as m 2 are material parameters (structure intensity and crack growth) about the SN curves.
Considering the stress range ∆σ can be performed based on action load ∆T m (e.g., torque) through the FEA method: where τ is the stress coefficient corresponding to the torque load of ∆T m . In addition, the stress intervals of fatigue critical components, which are divided into n groups, are assumed. The intervals number of failure stress falling to the group i each year is denoted by n i . (∆T m , n i ) is obtained by analyzing the statistical wave data of previous years. The fatigue failure state calculation equation based on signs using P-M rules can be expressed as: where K c 1 and K c 2 are the characteristic values of K 1 and K 2 (logK c i is the first order origin moment of logK i minus two standard deviation of logK i ), n ijk is the cycles number each year of stress intervals S ijk under a wave condition with H s i and T av j , P H s i , T av j is the probability of a particular wave condition determined from H s i and T av j .
The fatigue design factor (FDF) value, which is used to evaluate the reliability of fatigue critical details, can be expressed as: where T F is the fatigue life for fatigue design; T L is the service life of the offshore facility, typically 20 years for wave energy converter. The reliability of the device and its safety indices could be performed using the load value ∆T m obtained by finite element analysis, and the limit state equation is written [20]: where t indicates the service time (0 ≤ t ≤ T L ), K 1 and K 2 are stochastic variables determined by the SN curve, ∆ represents the modelling uncertainty based on P-M rules, S ijk = X M X SCF ∆σ represents the stress range given H s i and T av j . The modelling uncertainty X M , which is approximately follows LN (logarithmic normal) distributed, is composed of wave condition uncertainty and action load uncertainty. The uncertainty of SCF (stress concentration factor), X SCF , which is considered to follow LN distributed, contains the geometry uncertainty and fatigue critical location uncertainty. According to the particularity of working environment of wave energy converters, the representative SN curves of offshore installment with cathodic protection are considered [38]. Table 1 describes partial recommended parameters of stochastic variables applied to fatigue estimation. The values of the COV (coefficient of variation) are used to reduce the inaccuracy caused by the uncertainty of the stochastic variables.  [21,38] It is worth mentioning that the above stochastic variable models are taken into account the assumption that the annual variation of wave load can be ignored. However, developing appropriate control strategies can greatly reduce such variations.

Determination of Reliability Assessment Indices
Generally, material properties of the installment, wave loads as well as environmental features are uncertain. However, these factors are not considered in the traditional deterministic analysis. It is a feasible scheme to establish parameterized stochastic model by using a probability approach, which can consider the variables' uncertainty. The typical uncertainty of system behavior mainly includes epistemic uncertainties and aleatory uncertainties. The first may be avoided by acquiring more real operational data and more insight into the mechanism. However, the second is inevitable due to the physical properties of the material [20].
The probability of failure, P F , can be expressed as the following equation: Having acquired the probability expression of fatigue failure for critical components, the FORM can be applied to assess the corresponding reliability indices. FORM makes the stochastic parameters follow the standard normal distribution independent of each other by a special transformation approach. The basic principle of the approach is that stochastic variables can be approximately represented by their first-and second-order origin moments. The reliability indices could be evaluated by the moment of stochastic variables based on limit state equation. The probability function of failure can be expressed as: where Φ( ) represents the function follows standard normal distribution; β is the reliability index of the fatigue critical components, which is determined by the shortest failure distance in the state space. Mathematically, it can be expressed as [39]: where µ and σ are mean value and standard deviation, respectively; The subscripts R and L represent the strength and load effect, i.e., the first term and all the following terms in Equation (2). The typical values of β and the corresponding probability of failure for standard normal distribution, P F , as shown in Table 2. For time-dependent failure probabilities, P F (t), the annual probability of failure ∆P F (t), which is the probability of surviving reliably to time t, can be developed as follows [40]: where ∆t is equal to 1 year. Obviously, the annual reliability index ∆β = −Φ −1 (∆P F (t)).

Target Reliability Level for Fatigue Failure
The target reliability level needs to be defined when performing reliability-based analysis. Target reliability level is significantly determined by the impact of the failure mode and the implementation location, such as casualty risk and property damage. If the designed reliability indices are less than the target value, the fatigue critical components under consideration will not be able to withstand the applied alternating stress. The wave-energy converters have the characteristics of low environmental impact, no fire hazard and few casualties since they are operating unmanned. This indicates that the failure of WECs results mainly in economic influences. The acceptable annual reliability index (see Reference [20]) considering various cost of condition monitoring strategies on basis of the consequences for certain failure mode, as shown in Table 3. Because the WECs are operating unattended, the corresponding costs for condition monitoring are higher than for offshore oil and gas platforms. Hence, the acceptable target annual reliability values of 3.1-3.7 are appropriate, which is consistent with the target reliability evaluation range of OWTs. In addition to considering the annual reliability of the failure mode, focusing on the cumulative reliability over the entire life cycle is also an important indicator. The magnitude of β, typically from 2.3 to 3.1 for offshore structures, is largely determined by the implementation of condition monitoring based on the limit state equation [20]. The minimum values range of β can be mapped from OWTs to WECs because of the similar failure consequences (almost no casualties and environmental pollution). Table 4 shows the cumulative and annual reliability indices for OWTs in different corrosive environments [21], which can be used as a useful reference for fatigue design of WECs. The principles of fatigue design are determined by the extent of condition monitoring, the installation location of critical components and the consequences of fatigue failure. The floating WECs are far away from the land, and have been operated in the vast ocean for a long time. The fatigue critical components (e.g., turbine) are in an inaccessible area, in which is difficult to take measurements for condition monitoring during service. Therefore, the FDF values are recommended to be 3 according to Reference [41].
From what has been discussed above, the minimum accumulated reliability index of 2.4-3.1 and the minimum annual reliability index of 3.1-3.7 are appropriable for required values of fatigue design factor. These values can reduce the probability of failure and compensate the uncertainties of the device's operation.

Wave Characteristics and Fatigue-Critical Details
In this section, the primary purpose is to acquire the characteristic period of the planned deployment location by analyzing the wave climate of the target sea area (South China Sea), so as to provide theoretical support for the selection of PTO damping of the water hydraulic-drive WEC. Secondly, a whole FEA model of turbine blade is developed using the ANSYS 19.0 software package to obtain the fatigue-critical location and the maximum stress amplitude under alternating loads.

Analysis of Wave Climate
Ocean waves are generally considered to consist of variety stochastic wave conditions, which includes a series of wave sequences with different frequencies as well as amplitudes. In addition, the spectral moments function could be utilized to determine significant wave height, H s , and average energy period, T av , on the basis of the nth spectral moment [12].
where E(ω, θ) is the variance density spectrum, ω is the wave frequency; both of them are depended on the relationship of DSD (Doppler-Shifted-Dispersion). Therefore, the H s and T av of ocean waves can be estimated from the frequency domain as: and For deep water waves, the included power potential per unit wave-front can be calculated as follows [42]: where ρ is the seawater density, g is gravitational acceleration. The 22-year wave data in South China Sea, with the span of 1988-2009, were performed according to the WAVEWATCH-III (WW3) model [11]. Table 5 shows that the 3-hourly wave energy statistics data (times) and probability of occurrence (%) of H s and T av in South China Sea. The background color scale represents the occurrence probability of the wave condition.
As can be seen from Table 5, the most probable average significant wave height values are below 3.5 m and that most of values are in the 0.5-2.5 m range. The wave periods are basically scattered at the interval of 4.5-8.5 s.
According to the contents of Table 5 and Equation (30), the energy distribution probabilities, P e , of each wave condition can be obtained by the following Equation.
where P i and N i are the wave power and occurrences of the wave condition in group i, respectively; n is the total number of groups. The wave-power levels (kW/m) and energy distribution probabilities (%) at each wave condition can be estimated as shown in Table 6. The background color scale illustrates the energy distribution probability. The maximum wave power is distributed in H s of 1.5-3.5 m and T av of 6.5-8.5 s. It also can be seen that in wave condition with H s = 2.5 m and T av = 6.5 s, the largest energy contribution, which account for about 15.24% of the total wave energy, is obtained. Table 5. Three-hourly wave energy statistics data (upper data, times) and probability of occurrence (underlying data, %).   Figure 4 illustrates that the fitting curve of energy distribution probabilities under different wave frequencies. The energy contained in wave conditions with a frequency of 0.837 rad/s (average energy period T av = 7.5 s) accounts for about 25% of the total wave energy in South China Sea. Therefore, the proposed WEC should achieve maximum efficiency in this wave condition.

Finite-Element Analysis (FEA) Model of Turbine Blade
The single blade system of the turbine is composed of a fan blade and a limit switch with mechanical stop function, as depicted in Figure 5. The basic design parameters of the blade system are provided in Table 7.
Based on the purpose of illustrating the application scope of the proposed method, the turbine blade and limit switch are considered to be made of structural steel (S355) widely used in offshore implementations, in this study. Table 8 shows that the main performance parameters of S355 steel.

Finite-Element Analysis (FEA) Model of Turbine Blade
The single blade system of the turbine is composed of a fan blade and a limit switch with mechanical stop function, as depicted in Figure 5. The basic design parameters of the blade system are provided in Table 7.
Based on the purpose of illustrating the application scope of the proposed method, the turbine blade and limit switch are considered to be made of structural steel (S355) widely used in offshore implementations, in this study. Table 8 shows that the main performance parameters of S355 steel. The single blade system of the turbine is composed of a fan blade and a limit switch with mechanical stop function, as depicted in Figure 5. The basic design parameters of the blade system are provided in Table 7.
Based on the purpose of illustrating the application scope of the proposed method, the turbine blade and limit switch are considered to be made of structural steel (S355) widely used in offshore implementations, in this study. Table 8 shows that the main performance parameters of S355 steel.  In order to obtain the maximum stress position of the turbine blade, a numerical simulation was carried out using the Mechanical APDL solver in the software of ANSYS Workbench. Two different boundary conditions are considered in present paper, that is fixed support and moment loads. Both of the action planes of the limit switch are supported with fixed support to restrict its free rotation  In order to obtain the maximum stress position of the turbine blade, a numerical simulation was carried out using the Mechanical APDL solver in the software of ANSYS Workbench. Two different boundary conditions are considered in present paper, that is fixed support and moment loads. Both of the action planes of the limit switch are supported with fixed support to restrict its free rotation (see Figure 2). The moment load boundary conditions, whose values are obtained by analyzing the 42 wave conditions mentioned in Section 4.1, are applied to the blade system.
Mesh generation is an important part in the establishment of finite element model, and an unstructured mesh technology was used for generating the volume meshes. Generally, the precision of numerical calculations is improved by increasing the mesh number. However, a quite dense mesh results in burdensome calculations and seriously effects the efficiency of simulation. Therefore, it is reasonable to generate an appropriate number of meshes with acceptable computational accuracy.
A mesh independence validation was performed to find the appropriately converged solutions. The whole blade system was discretized into mesh cells of 22,700; 35,500; 55,100; 90,400; and 180,100 approximately by changing the base size of the mesh element. In this section, the turbine blade was analyzed under wave fatigue loads generated in the South China Sea, which has a significant wave height of 1.5 m and average energy period of 7.5 s. The estimated applied torque load is calculated from the analytical model with λ 1 = 1 and λ 2 = 6 (see Sections 2.2 and 2.3). Figure 6 shows the relationship of the maximum Von Mises stress and the number of elements. Specific mesh information as shown in Table 9. approximately by changing the base size of the mesh element. In this section, the turbine blade was analyzed under wave fatigue loads generated in the South China Sea, which has a significant wave height of 1.5 m and average energy period of 7.5 s. The estimated applied torque load is calculated from the analytical model with = 1 and = 6 (see Sections 2.2 and 2.3). Figure 6 shows the relationship of the maximum Von Mises stress and the number of elements. Specific mesh information as shown in Table 9.  The maximum Von Mises stress tends to be stable when the element size is 0.8 mm, which has a relative difference (0.31%) when compared to a denser meshing, as shown in Figure 6 and Table 8. Therefore, we choose to use a moderate number of mesh cells (90,400) to save computing resources and improve work efficiency.
Simulation analysis shows that the blade root bears the maximum Von Mises stress (consider stress concentration) and the blade tip obtains the maximum deformation when the blade system is  The maximum Von Mises stress tends to be stable when the element size is 0.8 mm, which has a relative difference (0.31%) when compared to a denser meshing, as shown in Figure 6 and Table 8. Therefore, we choose to use a moderate number of mesh cells (90,400) to save computing resources and improve work efficiency.
Simulation analysis shows that the blade root bears the maximum Von Mises stress (consider stress concentration) and the blade tip obtains the maximum deformation when the blade system is subjected to moment load, as shown in Figure 7. Therefore, the blade root is the fatigue critical detail of the blade system, and increasing its thickness can effectively improve the fatigue strength of the blade. subjected to moment load, as shown in Figure 7. Therefore, the blade root is the fatigue critical detail of the blade system, and increasing its thickness can effectively improve the fatigue strength of the blade.

Results and Discussion
The reliability assessment of water hydraulic drive WEC is carried out in this section. The reliability index of the proposed WEC is evaluated by means of frequency domain hydrodynamic analysis, finite element analysis, wave climate analysis and reliability analysis. The flowchart of the reliability analysis is illustrated in Figure 8.
In order to make the calculated results more specific, it is assumed that the structural parameters

Results and Discussion
The reliability assessment of water hydraulic drive WEC is carried out in this section. The reliability index of the proposed WEC is evaluated by means of frequency domain hydrodynamic analysis, finite element analysis, wave climate analysis and reliability analysis. The flowchart of the reliability analysis is illustrated in Figure 8.

Parameters Values Parameters Values
3.0 m 30° 1.2 m 12 m Figure 9. Plots of added mass and radiation damping versus energy period .

Optimized PTO Damping
Since the proposed water hydraulic drive WEC is planned to be installed in the South China Sea, the design parameters of the device must be adapted to the characteristics of the sea conditions in the designated sea area. According to the analysis in Section 4.1, the wave energy contained in the wave condition with a period of 7.5 s is the largest. Therefore, one of the focuses in this section is to make the WEC obtains the maximum power under this wave condition. Figure 10 illustrates the curves of the maximum power, , absorbed by the turbine and the relative displacement of the water column contained in the tube, , under different dimensionless coefficients * = / . Obviously, the maximum power absorbed by the water hydraulic drive WEC is obtained in = 5.4 , which can be achieved through a reasonable design of water turbine. In order to make the calculated results more specific, it is assumed that the structural parameters of the proposed floater of water hydraulic drive WEC are described in Table 10. The added mass µ 1 and the radiation damping coefficient B were calculated by the software ANSYS AQWA for different wave frequencies ω and the condition of infinite depth. The plots of µ 1 and B versus wave period T = 2πω is shown in Figure 9. It is easy to see that the radiation damping coefficient of the floater tends to rise first and then fall with the increase of the wave period, and the maximum value is obtained at T = 5.3 s. However, the added mass shows an opposite trend, and the minimum value is obtained at a wave period of 3.5 s.

Parameters Values Parameters Values
3.0 m 30° 1.2 m 12 m Figure 9. Plots of added mass and radiation damping versus energy period .

Optimized PTO Damping
Since the proposed water hydraulic drive WEC is planned to be installed in the South China Sea, the design parameters of the device must be adapted to the characteristics of the sea conditions in the designated sea area. According to the analysis in Section 4.1, the wave energy contained in the wave

Optimized PTO Damping
Since the proposed water hydraulic drive WEC is planned to be installed in the South China Sea, the design parameters of the device must be adapted to the characteristics of the sea conditions in the designated sea area. According to the analysis in Section 4.1, the wave energy contained in the wave condition with a period of 7.5 s is the largest. Therefore, one of the focuses in this section is to make the WEC obtains the maximum power under this wave condition. Figure 10 illustrates the curves of the maximum power, P max , absorbed by the turbine and the relative displacement of the water column contained in the tube, Y, under different dimensionless coefficients B * = B PTO /B. Obviously, the maximum power absorbed by the water hydraulic drive WEC is obtained in B PTO = 5.4B, which can be achieved through a reasonable design of water turbine. The equivalent PTO damping is invariable as the rotational angular velocity of the designed turbine is constant. Therefore, the optimized PTO damping obtained under the wave condition with the period of 7.5 s should be applied to other situations if the rotational speed control strategy is adopted. The relationship between the maximum absorbed power of the turbine and relative displacement of oscillating water column and the wave periods is plotted in Figure 11. It is shown that the optimized PTO damping can enable the turbine to obtain the maximum power when the wave period is 7.5 s. Thus, it is realized that the maximum absorption of wave energy under the most energetic wave conditions.

Fatigue Reliability
The reliability index can represent the reliability level of the devices, which is a probabilistic assurance of component reliability. In fatigue design, the reliability index is closely related to time because the cumulated cycles number of fatigue loads increased with the service time. Without The equivalent PTO damping is invariable as the rotational angular velocity of the designed turbine is constant. Therefore, the optimized PTO damping obtained under the wave condition with the period of 7.5 s should be applied to other situations if the rotational speed control strategy is adopted. The relationship between the maximum absorbed power of the turbine and relative displacement of oscillating water column and the wave periods is plotted in Figure 11. It is shown that the optimized PTO damping can enable the turbine to obtain the maximum power when the wave period is 7.5 s. Thus, it is realized that the maximum absorption of wave energy under the most energetic wave conditions. The equivalent PTO damping is invariable as the rotational angular velocity of the designed turbine is constant. Therefore, the optimized PTO damping obtained under the wave condition with the period of 7.5 s should be applied to other situations if the rotational speed control strategy is adopted. The relationship between the maximum absorbed power of the turbine and relative displacement of oscillating water column and the wave periods is plotted in Figure 11. It is shown that the optimized PTO damping can enable the turbine to obtain the maximum power when the wave period is 7.5 s. Thus, it is realized that the maximum absorption of wave energy under the most energetic wave conditions.

Fatigue Reliability
The reliability index can represent the reliability level of the devices, which is a probabilistic assurance of component reliability. In fatigue design, the reliability index is closely related to time because the cumulated cycles number of fatigue loads increased with the service time. Without considering the probability of occurrence of the wave condition with period , the cycle number,

Fatigue Reliability
The reliability index can represent the reliability level of the devices, which is a probabilistic assurance of component reliability. In fatigue design, the reliability index is closely related to time because the cumulated cycles number of fatigue loads increased with the service time. Without considering the probability of occurrence of the wave condition with period T av j , the cycle number, n ijk , of fatigue loads per year can be calculated by the following Equation: For structural steels, the stress and strain show a linear relationship when the stress in the components is less than Young's modulus. Therefore, the maximum stress of the turbine blade under the action of moment load is proportional to the applied moment. Based on the simulation analysis of the finite element model established in Section 4.2, the maximum Von Mises stress of the blade system under the action of unit moment load (Nm) is τ = 3.731MPa. A turbine performance analysis with variable configuration parameters, λ 1 λ 2 ω r , is carried out, which can be used to study the fatigue life of the blade system in combination with the SN-curve approach and limit state equation. By solving Equation (20), it can be obtained that the service life of the water turbine blade can exceed 20 years under the condition of λ 1 λ 2 ω r ≥ 839 (FDF = 3), as shown in Figure 12. The application value of this conclusion is that the target service life of the WEC can be realized through reasonable adjustment of the turbine configuration parameters on basis of the reliability approach. of the finite element model established in Section 4.2, the maximum Von Mises stress of the blade system under the action of unit moment load (Nm) is = 3.731MPa. A turbine performance analysis with variable configuration parameters, λ , is carried out, which can be used to study the fatigue life of the blade system in combination with the SN-curve approach and limit state equation. By solving Equation (20), it can be obtained that the service life of the water turbine blade can exceed 20 years under the condition of λ ≥ 839 (FDF = 3), as shown in Figure 12. The application value of this conclusion is that the target service life of the WEC can be realized through reasonable adjustment of the turbine configuration parameters on basis of the reliability approach. In view of the lower linear speed of wave motion, the rotational velocity of the direct drive water turbine is generally lower than other WECs, e.g., an oil hydraulic system. Therefore, the optimal configuration parameters of the water turbine considered in this paper are shown in Table 11. Table 11. Optimal configuration parameters of the water turbine.

Parameters Values Parameters Values Parameters
Values 3 2   6 47.1 rad/s The calculated cumulative reliability index for the WEC in stochastic wave condition is described in Figure 13. As predicted, the reliability index declined monotonously with service life. However, the blade system does not achieve the target reliability level, = 2.4, recommended for other marine projects (see Section 3.3). Therefore, a calibrated minimum target cumulative reliability index of 2.1 for moving components of water hydraulic drive WEC is recommended in this paper (FDF = 3). In view of the lower linear speed of wave motion, the rotational velocity of the direct drive water turbine is generally lower than other WECs, e.g., an oil hydraulic system. Therefore, the optimal configuration parameters of the water turbine considered in this paper are shown in Table 11. The calculated cumulative reliability index for the WEC in stochastic wave condition is described in Figure 13. As predicted, the reliability index declined monotonously with service life. However, the blade system does not achieve the target reliability level, β = 2.4, recommended for other marine projects (see Section 3.3). Therefore, a calibrated minimum target cumulative reliability index of 2.1 for moving components of water hydraulic drive WEC is recommended in this paper (FDF = 3).
The calculated cumulative reliability index for the WEC in stochastic wave condition is described in Figure 13. As predicted, the reliability index declined monotonously with service life. However, the blade system does not achieve the target reliability level, = 2.4, recommended for other marine projects (see Section 3.3). Therefore, a calibrated minimum target cumulative reliability index of 2.1 for moving components of water hydraulic drive WEC is recommended in this paper (FDF = 3).   Figure 14 shows that the annual reliability index trend of the WEC under the specific fatigue load. It can be observed that the turbine blade achieves the 20-year design life with the annual reliability index ∆β = 3.5. This is consistent with the range of annual reliability index, ∆β ∈ (3.1, 3.7), proposed in literature [20].
Energies 2019, 12,4189 18 of 21 Figure 14 shows that the annual reliability index trend of the WEC under the specific fatigue load. It can be observed that the turbine blade achieves the 20-year design life with the annual reliability index ∆ = 3.5. This is consistent with the range of annual reliability index, ∆ ∈ (3.1,3.7), proposed in literature [20]. In the present section, reliability assessment based on probabilistic method has been carried out for different configuration parameters of the water turbine, which can be used to improve the reliability design of the proposed WEC under the condition of meeting the target service life. It must be pointed out that the coefficients of variation of stochastic variables are derived from the recommendation of previous references, some of which assumed values for the proposed water hydraulic drive WEC when more realistic data was not available. In addition, the efficiency of the water turbine is not considered when the torque loads of the blade system are calculated in this work, and the reliability analysis results are relatively conservative. The detailed reliability data can be obtained through experimental methods and real-time monitoring systems, which should be conducive to establish more accurate reliability models of stochastic variables.

Conclusions
In this paper, a structure design and frequency domain analysis of the water hydraulic-drive WEC has been developed and used to assess the fatigue reliability of system-critical details. A reliability model based on the SN-curves and linear damage cumulative theory is proposed and solved, considering the uncertainty of the model parameters. The FORM is developed to assess the reliability indices, and the configuration parameters of the turbine blade system are obtained through reliability assessment. The following conclusions can be obtained based on previous analysis:

•
The frequency domain analysis method can effectively evaluate both the average and maximum power absorbed by the water hydraulic drive WEC, which provides a theoretical basis for the improvement of the system performance under random wave conditions. In the present section, reliability assessment based on probabilistic method has been carried out for different configuration parameters of the water turbine, which can be used to improve the reliability design of the proposed WEC under the condition of meeting the target service life. It must be pointed out that the coefficients of variation of stochastic variables are derived from the recommendation of previous references, some of which assumed values for the proposed water hydraulic drive WEC when more realistic data was not available. In addition, the efficiency of the water turbine is not considered when the torque loads of the blade system are calculated in this work, and the reliability analysis results are relatively conservative. The detailed reliability data can be obtained through experimental methods and real-time monitoring systems, which should be conducive to establish more accurate reliability models of stochastic variables.

Conclusions
In this paper, a structure design and frequency domain analysis of the water hydraulic-drive WEC has been developed and used to assess the fatigue reliability of system-critical details. A reliability model based on the SN-curves and linear damage cumulative theory is proposed and solved, considering the uncertainty of the model parameters. The FORM is developed to assess the reliability indices, and the configuration parameters of the turbine blade system are obtained through reliability assessment. The following conclusions can be obtained based on previous analysis:

•
The frequency domain analysis method can effectively evaluate both the average and maximum power absorbed by the water hydraulic drive WEC, which provides a theoretical basis for the improvement of the system performance under random wave conditions. • The wave condition with H s = 2.5 m and T av = 6.5 s contains the largest energy contribution, which accounts for about 15.24% of the total wave energy in the South China Sea. The energy contained in the wave condition with a period of 7.5 s accounts for about 25% of the total energy.

•
The blade root is the fatigue critical detail bearing the maximum Von Mises stress, and the blade tip obtains the maximum deformation under the moment load. The reliability assessment method employed and analyzed in the present paper is able to effectively evaluate the reliability of the water hydraulic-drive WEC and calibrate the preliminary structure parameters based on probabilistic methods. The proposed method can be extended to relevant engineering design of marine-energy development and utilization.
Author Contributions: H.L. and S.T. were responsible for conceptualization, date curation, writing and formal analysis. W.W. was responsible for funding acquisition, investigation and validation. L.M. and H.M. were responsible for project administration and visualization. G.Z. and J.L. was responsible for methodology and software.