Next Article in Journal
Novel Cuboid-like Crystalline Complexes (CLCCs), Photon Emission, Fluorescent Fibers, and Bright Red Fabrics of Eu3+ Complexes Adjusted by Amphiphilic Molecules
Next Article in Special Issue
Modelling of Environmental Ageing of Polymers and Polymer Composites—Durability Prediction Methods
Previous Article in Journal
Effect Mechanism and Simulation of Voids on Hygrothermal Performances of Composites
Previous Article in Special Issue
Influences of Modified Sm2O3 on Thermal Stability, Mechanical and Neutron Shielding Properties of Aminophenol Trifunctional Epoxy Resin
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thickness-Prediction Method Involving Tow Redistribution for the Dome of Composite Hydrogen Storage Vessels

1
Hubei Key Laboratory of Advanced Technology for Automotive Components, Wuhan University of Technology, Wuhan 430070, China
2
Hubei Engineering Research Center for Green & Precision Material Forming, Wuhan 430070, China
3
Hubei Collaborative Innovation Center for Automotive Components Technology, Wuhan 430070, China
*
Author to whom correspondence should be addressed.
Polymers 2022, 14(5), 902; https://doi.org/10.3390/polym14050902
Submission received: 30 December 2021 / Revised: 5 February 2022 / Accepted: 18 February 2022 / Published: 24 February 2022

Abstract

:
Traditional thickness-prediction methods underestimate the actual dome thickness at polar openings, leading to the inaccurate prediction of the load-bearing capacity of composite hydrogen storage vessels. A method of thickness prediction for the dome section of composite hydrogen storage vessels was proposed, which involved fiber slippage and tow redistribution. This method considered the blocking effect of the port on sliding fiber tows and introduced the thickness correlation to predict the dome thickness at polar openings. The arc length corresponding to the parallel circle radius was calculated, and then, the actual radius values corresponding to the bandwidth were obtained by the interpolation method. The predicted thickness values were compared with the actual measured thickness. The maximum relative error of the predicted thickness was 4.19%, and the mean absolute percentage error was 2.04%. The results show that the present method had a higher prediction accuracy. Eventually, this prediction method was used to perform progressive damage analysis on vessels. By comparing with the results of the cubic spline function method, the analysis results of the present method approached the actual case. This showed that the present method improved the accuracy of the design.

Graphical Abstract

1. Introduction

As pollution-free energy, hydrogen energy has become a new energy to replace traditional fossil energy. The burning of conventional fossil fuels produces greenhouse gas and harmful gases, while hydrogen only produces water [1]. Composite hydrogen storage vessels, as the primary equipment of high-pressure hydrogen storage, have the advantages of lightweight and reliable performance [2]. The wide application of hydrogen energy is inseparable from hydrogen storage vessels.
The filament-winding process has been extensively researched and applied [3,4]. The filament-winding process is usually used in the manufacture of hydrogen storage vessels. The geometry and mechanical properties analysis of hydrogen storage vessels, especially in the dome section, are very complicated because of the characteristics of the process and materials. In addition, the thickness and winding angle of composite layers at the dome section of vessels will change significantly, which affects the mechanical properties of composite layers [5,6].
Nowadays, finite-element analysis technology plays an increasingly important role in designing high-performance composite hydrogen storage vessels [7]. The more accurate the hydrogen storage vessel model established, the more accurate the results of the finite-element analysis will be [8]. The changes in angle and thickness at the dome section need to be considered in the accurate performance prediction of a hydrogen storage vessel. It is difficult to predict the dome thickness of vessels because of fiber slippage during the winding process. Therefore, accurate thickness prediction is a prerequisite for the performance prediction of hydrogen storage vessels. The design reliability of vessels can be increased by accurate thickness prediction and finite-element analysis, achieving more effective optimization schemes. The dome thickness of hydrogen storage vessels was predicted to establish the finite-element model to research the performance of the composite hydrogen storage vessel by researchers [9].
However, traditional thickness-prediction methods were only accurate outside the two bandwidths of polar openings. The influence of fiber width was not considered, which caused the predicted thickness value at polar openings to be infinite and did not match actual measurement thickness [10]. Wang et al. [11,12] proposed a thickness-prediction method, which used a cubic spline function and assumed that the amounts of the composite material were constant within the two bandwidths of polar openings. The method provided more accurate thickness-prediction results, but it did not consider the slippage effect of fiber tows during the winding process, which made prediction results at polar openings far lower than actual measurement thickness.
In recent years, many researchers also proposed new thickness-prediction methods. Leh et al. [13] combined existing winding-angle and thickness-prediction methods [8,11,12] to propose a model for predicting the dome geometric characteristics of multi-layer composite pressure vessels, but the problems existing in the existing methods was not improved. Zu et al. [14,15] used a three-dimensional laser scanner to measure the vessel’s profile after winding and obtained the thickness distribution of the dome section. However, the method needed thickness values provided by experiments, which limited the application of this method. Che et al. [16] proposed a prediction method for the thickness of the dome section. This method took the effect of fiber overhead caused by hoop winding, but this method did not consider the influence of bandwidth on the fiber thickness at polar openings. Akhtar et al. [17] presented the analytical modeling involving kinematic constraints based on a geodesic trajectory to predict thickness change of dome section. However, the method did not consider the slippage effect of fiber tows during the winding process, and it needed to divide the dome section into multiple zones for calculation. The number of divided zones affected prediction accuracy.
The predicted thickness values could provide a basis for the establishment of a finite-element model (FEM) of composite hydrogen storage vessels. A progressive damage model has been widely used in progressive damage and failure analysis of composite masteries [18]. Composite material progressive damage analysis method could simulate complex stress conditions and material damage state of composite vessels under internal pressure. Burst pressure and burst mode of composite hydrogen storage vessels could be better predicted by progressive damage analysis [19,20].
The progressive damage analysis of composite hydrogen storage vessels included the establishment of FEM and the analysis of results. ANSYS ACP was used by Kangal et al. [21] to analyze the progressive damage of the glass fiber/carbon fiber hybrid composite pressure vessel. The results of hydrostatic pressure tests were compared with the results of finite-element simulation analysis and showed that numerical simulation results were in good agreement with experimental results. However, the thickness and angle variations were obtained from CADWIND, and the simulation results of dome thickness obtained from CADWIND were inaccurate. A multiscale progressive damage analysis method was proposed by Lin et al. [22], which introduced representative volume elements (RVEs). Puck failure criterion was adopted to predict the burst pressure of composite pressure vessels, but the vessel was modeled by the plug-in WCM. There were some differences between the actual shell composite dimensions and those modeled by WCM. A progressive damage method was used by Alam et al. [23] to investigate the influence of winding angle, the number of layers, and the sequence of layers on the burst strength of composite pressure vessels. However, the method only compared the burst pressure values corresponding to different winding angles, numbers of layers, and sequences of layers without further research on its causes. ANSYS ACP software was used by Sharma et al. [24] to perform progressive failure analysis of composite hydrogen storage vessels. The influence of the interlayer hybridization effect on the burst pressure of composite vessels was analyzed. However, the dome section was divided into multiple areas in that model, resulting in a discontinuous angle and thickness of the dome section, which was different from the actual angle and thickness. A composite material progressive damage model was established by Hu et al. [25] using ABAQUS software. The Hashin failure criterion was used to predict the burst pressure of composite hydrogen storage vessels with different layer-up sequences. The liner used in the article was different from the conventional vessel liner, and the conclusion may not apply to others.
Progressive damage analysis could better simulate the complex stress condition and material damage state of composites [26]. However, previous literature [27,28] mostly used WCM or CADWIND software to model the composite hydrogen storage vessels, and the thickness values at the dome section could not be predicted well by the software, thus affecting the accuracy of the FEM. To accurately predict the burst pressure of composite hydrogen storage vessels, it is beneficial to study a more accurate prediction method of dome thickness.
In this paper, a method of thickness prediction for the dome section of composite hydrogen storage vessels was proposed, which involved the effect of fiber slippage and the redistribution of fiber tows. This method solved the problem that the traditional cubic spline function method underestimated the dome thickness at polar openings. On this basis, this prediction method was used to perform progressive damage analysis on composite hydrogen storage vessels. The influence of the reaming process and the layer-up sequence was studied on load-bearing capacity and damage evolution of hydrogen storage vessels. The analysis results of the present method were compared with the results of the cubic spline function method.

2. Analytical Model for Dome Thickness Prediction

An accurate dome-thickness-prediction model was a key point in the design process of composite hydrogen storage vessels. The results accuracy of finite-element analysis was affected by the prediction accuracy. The cubic spline function prediction method proposed by Wang et al. [12] was a relatively accurate prediction method. However, the thickness-prediction values of the cubic spline function method at polar openings were much lower than the actual measurement. The effect of fiber slippage on dome thickness at polar openings was not considered. The cubic spline function may become undefined when using a large initial winding angle and large bandwidth. A new prediction method is proposed in the following text.

2.1. Dome Thickness at Polar Openings

Figure 1 shows a schematic diagram of fiber slippage during the winding process. Two situations need to be considered when predicting the dome thickness at polar openings: (1) dome thickness at polar openings before reaming; and (2) dome thickness at polar openings after reaming. The dome thickness at polar openings should be calculated separately based on the above two cases.

2.1.1. Dome Thickness before Reaming

The wet winding process was usually used to manufacture hydrogen storage vessels. The fiber tows slid to both sides under the action of winding tension and resin infiltration when fiber tows were wound around polar openings. The redistribution of fiber tows was caused by the effect of fiber slippage. However, the fiber tows cannot slide further because of the blocking effect of the port, which leads to the fiber stacking around polar openings. The thickness of polar openings cannot be calculated according to the theory of the same amounts of fiber tows. The thickness at polar openings and the thickness one bandwidth away from polar openings were affected by factors such as the slipping, realignment, and roving separation of fiber tows. There was a correlation between the thickness of the two places, considering the redistribution of the fiber tows near polar openings. The formula for calculating the dome thickness at polar openings before reaming was expressed as follows:
t ( r ) = m R × n R π × arccos ( r 0 r b ) × t p × γ
where m R stands for the number of fiber tows at the cylindrical section, n R stands for the number of layers in the cylindrical section, r 0 is the radius of the polar opening, r is the radius of the parallel circle, t p is the thickness of a fiber tow, b is the width of a fiber tow, γ is the thickness correlation coefficient, and the value is 0.2~0.3.

2.1.2. Dome Thickness after Reaming

Fiber stacking occurred near polar openings if helical layers were wound at a single helical angle. Winding difficulties and stress concentration in composite layers were caused by fiber stacking, which affected the load-bearing capacity of the vessels. Therefore, the reaming process was used to reduce thickness accumulation and fiber stacking at the dome section. The observation of the fiber slippage effect through actual winding experiments is shown in Figure 1. The fiber tows during the winding process after the reaming process had no longer been obstructed by the port. However, the resin-impregnated fiber tows still slid to both sides under the action of winding tension. As a result, the thickness of the fiber tows around actual polar openings decreased, and the width of fiber tows increased. From the winding measurement, the bandwidth used in the calculation in this paper was 1.2 times the actual bandwidth, and the thickness was divided by 1.2.

2.2. Revised Model of Radius Corresponding to Bandwidth

The radius of the cubic spline function method was calculated by projecting fiber tows trajectory onto the horizontal plane. The formula ( r b = r 0 + b , r 2 b = r 0 + 2 b ) was workable when polar openings radius and bandwidth were small. However, the radius of polar openings of composite layers increased after the reaming process. The projected length of fiber tows bandwidth on the horizontal plane decreased, as shown in Figure 2. The traditional cubic spline function was not suitable to predict the thickness of the dome section. The actual radius values of r b , r 2 b corresponding to the bandwidth were calculated by the actual arc length in this paper. MATLAB R2020b [29] software calculated the arc length corresponding to parallel circle radius. The values used by the cubic spline function prediction method and the present method are shown in Figure 2. Thereafter, the values of r b , r 2 b corresponding to bandwidth were calculated by the interpolation method.

2.3. Revised Model of Fiber Volume within Two Bandwidths

When using the cubic spline function prediction method, the thickness was set as a fixed value during the calculation process for amounts of composite material within two bandwidths of polar openings. In the present method, the values of dome thickness were variable, which varied with the radius of the parallel circle. The amounts of fiber tows within two bandwidths could be expressed as follows:
V c o n s t = r 0 r b 2 π r × m R × n R π × arccos ( r 0 r ) × t p d r +     r b r 2 b 2 π r × m R × n R π × [ arccos ( r 0 r ) arccos ( r 0 + b r ) ] × t p d r
The analytical model for dome thickness prediction could be expressed as follows:
t ( r i ) = m 1 × r i 0 + m 2 × r i 1 + m 3 × r i 2 + m 4 × r i 3 ( r 0 r R )
The unknown coefficients m 1 , m 2 , m 3 , and m 4 could be solved from:
[ m 1 m 2 m 3 m 4 ] = [ 1 r 0 r 0 2 r 0 3 1 r 0 r 2 b 2 r 2 b 3 0 1 2 r 2 b 3 r 2 b 2 π ( r 2 b 2 r 0 2 ) 2 π 3 ( r 2 b 3 r 0 3 ) ( r 2 b 4 r 0 4 ) ( r 2 b 5 r 0 5 ) ] 1 × [ m R × n R π × arccos ( r 0 r b ) × t p × γ m R × n R π × [ arccos ( r 0 r 2 b ) arccos ( r b r 2 b ) ] × t p m R × n R π × ( r 0 r 2 b × r 2 b 2 r 0 2 r b r 2 b × r 2 b 2 r b 2 ) × t p V c o n s t ]
where Vconst could be calculated by:
V c o n s t = r 0 r b 2 π r × m R × n R π × arccos ( r 0 r ) × t p d r + r b r 2 b 2 π r × m R × n R π × [ arccos ( r 0 r ) arccos ( r 0 + b r ) ] × t p d r

2.4. Verification of the Analytical Model

The profiles of three composite vessels in published literature [11,12] were selected to validate the accuracy of the analytical model for dome thickness prediction. The geometric parameters of composite vessels are shown in Table 1. The prediction results are shown in Figure 3. The actual measured thickness values were obtained from the literature [11,12]. The present method had a better agreement with the actual measured thickness values at polar openings and within two bandwidths than the cubic spline function method. Table 2 shows the relative error between predicted values and actual values of type A vessel within two bandwidths.
The double formula method [30] provided a reference for these models. It can be seen in Table 2 that the maximum relative error of the double formula method was 20.23%, the mean absolute percentage error was 6.66%, and the root-mean-square error was 1.95 mm. The maximum relative error of the cubic spline function method was 5.15%, the mean absolute percentage error was 3.22%, and the root-mean-square error was 0.70 mm. The maximum relative error of the present method was 4.19%, the mean absolute percentage error was 2.04%, and the root-mean-square error was 0.48 mm. Compared with the cubic spline function method, the maximum relative error of the present method was reduced by 18.6%, the mean absolute percentage error by 36.6%, and the root-mean-square error by 31.4%. The present method predicted the thickness distribution of the dome section more accurately.

3. Numerical Modeling

3.1. Winding Angle and Thickness

The winding angle and thickness were needed in the finite-element modeling of composite hydrogen storage vessels. The winding angle along the cylindrical section directly followed the design’s winding angle. The winding angle of the helical layers at the dome section changed from the initial winding angle (α0) to 90°. Therefore, the angle of each point of the dome section along the meridian was different. The geodesic trajectories were used for the winding process because of the use of equal domes. The helical winding angle over the dome section was generally assumed to follow geodesic formula [31]:
a = arcsin ( r 0 / r )
The values of dome thickness were obtained with the analytical model from Section 2. The winding angle and dome thickness distribution at the dome of the vessel are shown in Figure 4.
The helical winding angle and the minimum burst pressure of vessels were needed when calculating the thickness of the composite layers of the cylindrical section. The safety factor was set as 2.25 for the 70 MPa hydrogen storage vessel in this paper. Therefore, the minimum burst pressure should not be less than 2.25 times the working pressure.
The geodesic winding method was adopted since a pressure vessel with equal domes was used. The minimum angle of helical winding was calculated by the geodesic formula. The netting theory was commonly used to describe the minimum thickness tαk and tθk of helical and hoop layers, and its expression was:
{ t α k = R p b 2 k s σ f b cos 2 α t θ k = R p b 2 σ f b ( 2 1 k s tan 2 α )
t k = t α k + t θ k = ( 2 k s + 1 ) R p b 2 k s σ f b
where tαk and tθk are the thicknesses of helical layers and hoop layers, tk is the total thickness of composite layers at the cylindrical section, α is the helical winding angle, R is the radius of the cylindrical section, and pb is the minimum burst pressure.
The strength at the dome section was only determined by helical layers. The stress equilibrium factor ks needed to be introduced to increase the thickness of helical layers to strengthen the dome section. The stress equilibrium factor ranged from 0.6 to 0.8 to ensure that the burst point was at the cylindrical section.
The thickness of helical layers increased when the stress equilibrium factor was smaller. The burst point occurred in the cylindrical section of vessels. The burst mode of vessels was safer, but the thickness of composite layers increased. The value of the stress equilibrium factor was 0.7 in this paper. The initial winding angle α = 20°, the fiber tape width b = 6 mm, and the fiber tape thickness tp = 0.125 mm. A total of 26 hoop layers and 24 helical layers were wound.

3.2. Layer-Up Design Schemes

Different layer-up design schemes were used, as shown in Figure 5, to study the influence of the reaming process and layer-up sequence on the load-bearing capacity of vessels. The schemes included different layer-up sequences and helical winding angles. The layer-up schemes were divided into three groups A, B, and C, and each group included five schemes. The thickness of helical layers and hoop layers in these schemes was calculated by the netting theory. All layer-up design schemes had the same number of layers and thicknesses. The scheme A_1 was the initial design scheme. The schemes of group A used single reaming and increased reaming layers. The schemes of group B used double reaming and increased reaming layers. The schemes of group C changed the layer-up sequence of helical layers and hoop layers based on scheme A_4. Each column in Figure 5 represents a complete hoop layer (±90°) or helical layer (±α). The horizontal axis represents the number of layers from inside to outside, and the vertical axis represents the winding angle.

3.3. Finite-Element Model

The FEM of the composite hydrogen storage vessel in this paper was a 70 MPa type III vessel. The liner was made of aluminum alloy 6061-T6, as shown in Figure 6, to meet the requirements of lightweight and hydrogen corrosion resistance of the hydrogen storage vessel. The length of the aluminum liner was 500 mm, the outer diameter was 105 mm, the inner diameter was 101 mm, and the dome section height was 32.7 mm.
As shown in Figure 7, a finite-element model of the type III storage vessel was established to predict the burst pressures of composite hydrogen storage vessels. The cyclic symmetry was applied to improve computer efficiency, and a twelfth model was used to analyze vessels. The internal pressure was applied to the internal surface of the aluminum liner. The axial displacement of the vessel was restricted. Aluminum alloy 6061-T6 was used as the liner of vessels. The liner was modeled using a twenty-node cubic solid element (SOLID186) considering the nonlinearity of aluminum alloy material. The T700 fiber-reinforced composite was used as the composite layers of vessels. The composite layers were modeled using an eight-node structural solid shell element (SOLSH190). The aluminum liner and the composite layers were bonded together.
The liner of the vessel was made of aluminum alloy (6061-T6), and the material model was isotropic and elastoplastic. The bilinear isotropic hardening model was used for the plastic region. The material parameters and stress-strain curve of the aluminum liner obtained from the literature [14] are shown in Figure 8 and Table 3. The composite layers of the vessel were made of T700 fiber-reinforced composites, and the material model was transversely isotropic. The material parameters of the T700 fiber-reinforced composites obtained from the literature [14] are listed in Table 4 and Table 5.

3.4. Progressive Damage Model

A progressive damage model was used to predict structural damage of composite layers and the burst pressure of vessels. The model could accurately predict mechanical responses and the ultimate load-bearing capacity of composite materials. A damage initiation criterion was defined to determine the failure of composite material. The damage evolution law controlled the stiffness reduction of the materials. This progressive damage analysis process is shown in Figure 9. First, a finite-element model of the composite hydrogen storage vessel was established. The internal pressure of the vessel was loaded, and stress calculation was performed. Puck failure criterion was used. A material-stiffness degradation model was adopted according to the corresponding failure mode after the material was damaged. Finally, the calculation was terminated when the hydrogen storage vessel was structurally damaged.
The Puck failure criteria [32] were employed to predict the damage onset of the composite layers. The failure model of composite layers was generally classified as fiber tension and compressive breakage, matrix tensile, and compressive failure.
The Puck failure criterion included fiber failure (FF) and inter-fiber failure (IFF). The fiber failure criterion was written as follows:
f E , F F = 1 ± R t , c [ σ 1 ( ν ν f × m σ f E E f ) ( σ 2 + σ 3 ) ] w i t h                                   R t           f o r   [ ] 0                                         R c           f o r   [ ] 0
where f E , F F is the fiber failure stress exposure of the lamina; R t is the tensile strengths of the materials; R c is the compressive strengths of the materials; σ 1 , σ 2 and σ 3 are the normal stresses in a unidirectional layer; E , E f are the modules of the composite material and fiber direction; ν , ν f are the major Poisson’s ratios of the lamina and fibers; and m σ f is the stress magnification factor for transverse stresses in the fibers and was set as 1.1.
The formulations for the stresses σ n , τ n t and τ n 1 in an arbitrary plane with the inclination angle θ were:
σ n ( θ ) = σ 2 × cos 2 θ + σ 3 × sin 2 θ + 2 × τ 23 × sin θ × cos θ τ n t ( θ ) = ( σ 3 σ 2 ) × sin θ × cos θ + τ 23 × ( cos 2 θ sin 2 θ ) τ n 1 ( θ ) = τ 31 × sin θ + τ 21 × cos θ
The Puck IFF criterion was written as follows:
For σn(θ) ≥ 0:
f E , I F F ( θ ) = [ ( 1 R A t p ψ t R ψ A ) 2 σ n ( θ ) ] + ( τ n t ( θ ) R A ) 2 + ( τ n 1 ( θ ) R ) 2       + p ψ t R ψ A σ n ( θ )
For σn(θ) < 0:
f E , I F F ( θ ) = ( τ n t ( θ ) R A ) 2 + ( τ n 1 ( θ ) R ) 2 + ( p ψ c R ψ A σ n ( θ ) ) 2       + p ψ c R ψ A σ n ( θ )
R A t = R t   ;   R A = R   ;   R A = R c 2 ( 1 + p c ) p ψ t , c R ψ A = p t , c R A cos 2 ψ + p t , c R A sin 2 ψ cos 2 ψ = τ n t 2 τ n t 2 + τ n 1 2   ;   sin 2 ψ = τ n 1 2 τ n t 2 + τ n 1 2 p c = 1 2 cos 2 θ 1   ;   p c R A = p c R cos θ f p = ( τ 21 σ 2 ) 2 ( R A R ) 2 + 1 2 ( 1 + p c )
where R A t and R A are the tensile and shear strengths in the transverse direction of the fiber; R A is the shear strengths in the longitudinal direction of the fiber; R c is the compression strength perpendicular to the fiber direction; p t , c , p t , c are the inclination parameters; and θ f p is the fracture angle of the failure plane. The strength values of the Puck failure criterion are listed in Table 6.
The failure occurred when stress risk factor f E , F F or f E , I F F was greater than 1. Stiffness degradation of the failure elements was performed according to different stiffness-degradation models. Stiffness-degradation criteria [33] are depicted in Table 7.

4. Results and Discussion

The loading history for vessels, as shown in Figure 10, was used in this paper. The working pressure of hydrogen storage vessels was 70 MPa. The safety factor was set as 2.25. The minimum burst pressure was 157.5 MPa. If the hydrogen storage vessels were not structurally damaged at 157.5 MPa, the internal pressure was increased until burst to obtain the ultimate burst pressure of the hydrogen storage vessels.

4.1. Load-Bearing Capacity and Damage Behavior

The explosion simulations were carried out to predict the burst pressures of hydrogen storage vessels. Two burst modes, as shown in Figure 11, may occur; i.e., the safe burst mode and unsafe burst mode. The safe burst mode was attributed to fiber breakage in hoop fibers. The failure appeared in the cylindrical section while the port of the dome section entered the vessel. During this burst, the axial displacement was decreased rapidly after its value reached the maximum displacement point. The unsafe burst mode was attributed to fiber breakage in helical fibers. The failure appeared in the dome section, and the port was ejected at high speed.
Figure 12 shows the burst pressure and the burst mode of the vessels of different layer-up sequences. The results showed that the reaming process reduced fiber thickness near polar openings. The reaming process improved the stress distribution of composite layers, which increased the burst pressure of vessels. The reaming process increased the helical winding angle. The larger the helical winding angle was, the more hoop-bearing capacity could be increased. The schemes A_1 and B_3 showed that the burst pressure of the hydrogen storage vessel with reaming was increased by 9.15%.
The polar openings of the dome section bore not only the internal pressure of its area but also the internal pressure of the port area. A greater bearing capacity was required in composite layers near polar openings than in the other part. The scheme B_5 showed that the load-bearing capacity of the cylindrical section was improved by the reaming process. However, the strength of composite layers at polar openings was reduced by the reaming process, which made the vessel change from the safe mode to the unsafe mode. The load-bearing capacity of the vessel could be improved by using a reasonable reaming process under the premise that the vessel was in safe mode. From the results of groups A and B, it was suggested that the number of reaming layers was less than half of the design number of helical layers.
Figure 13 shows fiber stress of schemes A_4, C_4, and C_1 along the thickness direction of composite layers at the cylindrical section at the test pressure. The helical layers stress was less than the hoop layers stress in the cylindrical section under internal pressure. The fiber strength of helical layers was not fully utilized. The hoop fibers were broken first when the burst appeared in the cylindrical section. The outer portion of composite layers, which consisted of hoop layers, could increase the average stress of the helical layers of the cylindrical section by 5.92% and reduce the average stress of hoop layers by 3.65%.
The ultimate load-bearing capacity of vessels was determined by the maximum stress in the fiber direction. The smaller the stress was, the safer the vessel was. The bearing capacity of vessels could be increased by raising the fiber strength utilization factor of helical layers. The increase in bearing capacity could increase the burst pressure of the hydrogen storage vessels. A gradient change of stress along the thickness direction was caused by the adoption of hoop and helical alternated model, and the stress decreased as the number of layers increased. A large stress difference between the inner and the outer of hoop layers prevented outer hoop layers from completing enough work under the burst pressure, which was not conducive to the utilization of fiber strength. The safest and the highest burst pressure was reached by changing the layer-up sequence as follows: the hoop and helical layers were separated, the winding angles of helical layers were changed from small to large, and the inner portion of the composite layers consisted entirely of helical layers.
Figure 14 shows the axial displacement (axial displacement at polar openings of the dome) and the radial displacement (radial expansion of the cylindrical section) curves of schemes B_3 and B_5 under internal pressure loading. The burst mode of vessels could be predicted by the curves. The axial displacement curve of the scheme B_3 presented the burst point, which showed that it was under the safe burst mode. The axial displacement of scheme B_5 continued increasing, which showed it was under the unsafe burst mode.
The radial and axial displacement curves of the two schemes showed both an obvious linear and nonlinear phase before the vessels were structurally damaged. The initial linearity represented that the material was in the elastic phase. The nonlinear state of curves was caused by elastoplastic deformation of the aluminum alloy inner liner and the damage of composite layers. There were four regions before the vessel was structurally damaged, as shown in the curves of scheme B_3. The first linear curve before 85 MPa denoted the initial elastic behavior of the aluminum liner. The second nonlinear curve before 113 MPa denoted the plastic deformation of the aluminum liner. The third nonlinear curve before 158 MPa indicated resin matrix damage, but the composite vessel still had a great load-bearing capacity. Finally, a sharply rising curve showed that fiber damage occurred at about 158 MPa, and the vessel burst at about 172.9 MPa. The fracture behavior of the fiber determined the load-bearing capacity of vessels, but the matrix damage had little effect on the ultimate load-bearing capacity of vessels.

4.2. Damage-Evolution Laws

As shown in Figure 15a–c, the fiber damaged under the burst pressure of schemes C_3, C_4, and C_5, which were all under the unsafe burst mode. The damage of scheme C_3 appeared at polar openings of the dome section. The damage of scheme C_4 appeared at the connection between cylindrical and dome section, and the damage of scheme C_5 appeared at polar openings of the dome section and the connection between cylindrical and dome section. The fiber overhead occurred when helical angles were changed from large to small, or the inner portion of composite layers consisted entirely of hoop layers. The fiber overhead caused stress concentration and reduced the load-bearing capacity of the dome section, and the unsafe burst mode occurred.
The fiber damage of schemes B_3 and B_5 under the burst pressure is shown in Figure 15d,e. The fiber of hoop layers in the cylindrical section of scheme B_3 was completely damaged under the burst pressure, causing penetration damage. The fiber of helical layers at the dome section was partially damaged without causing structural damage. However, the fiber damage at the dome section of scheme B_5 accumulated continuously to form a penetration trend. The fiber of the cylindrical section was not damaged when structural damage occurred in the hydrogen storage vessel.
Figure 16 shows the damage evolution of fiber and matrix in the scheme B_3. The fiber was mainly subject to tensile stress, and almost no compression failure occurred when the hydrogen storage vessel was under working pressure. Therefore, only fiber tensile damage and matrix tensile and compression damage were shown. A matrix failure occurred at polar openings and the connection between cylindrical and dome sections when internal pressure increased to 113 MPa. The low strength of resin caused the matrix damage. As internal pressure increased, matrix damage developed from the inner to the outside of composite layers. Although the matrix in the dome section had been severely damaged, fiber was not yet broken. The hydrogen storage vessel had no structural damage. The fiber failure occurred at polar openings of the dome section with internal pressure increasing to 158 MPa. This was not enough to cause structural damage to the hydrogen storage vessel because the fiber damage area was small. The fiber tensile failure extended to other areas as the internal pressure increased. Finally, the hoop fibers of the cylindrical section were damaged with an internal pressure increase to 172.9 MPa. The hoop fibers’ damage resulted in stiffness degradation, and the damage rapidly expanded from the inside to the outside of the composite layers. The damage formed a penetrating area and caused structural damage to the hydrogen storage vessel. Research on the damage evolution behavior of the hydrogen storage vessel showed that the ultimate load-bearing capacity of hydrogen storage vessels largely depended on the fiber strength. The matrix had little effect on the burst pressure of hydrogen storage vessels.

4.3. Stress Distribution of the Aluminum Liner

Figure 17 shows the maximum stress of the aluminum liner under different pressures. Each column represents the maximum stress of the aluminum liner. The stress reduction was defined as the change of the maximum stress of the aluminum liner before and after autofrettage at the working pressure. The stress amplitude was defined as half of the difference between the maximum stress of the aluminum liner under the working pressure and the zero pressure.
The maximum stress and the stress amplitudes under the working pressure after autofrettage showed a downward trend as the reaming times and layers increased. The stress distribution of the aluminum liner under the working pressure was improved after the reaming process. The stress reduction in the aluminum liner under the working pressure improved the fatigue performance of vessels. Scheme B_3 showed that the maximum stress of the aluminum liner was reduced by 20.59% after the reaming process, and the stress amplitude was reduced by 9.28%. As the reaming times and layers increased further, the strength of the dome section decreased, and the burst mode changed from the safe mode to the unsafe, but the maximum stress and the stress amplitudes just decreased slightly.
Group C showed that the stress distribution of the aluminum liner was affected by the layer-up sequence. The maximum stress and the stress amplitude of scheme C_1, which adopted hoop and helical layers alternated model, were improved the most, except for the schemes of the unsafe burst mode. The fatigue performance of hydrogen storage vessels was improved by adopting the alternated model. However, the burst pressure of scheme C_1 was only 153.62 MPa. Scheme A_4, which had the same layers and winding angles as scheme C_1, adopted the model with hoop and helical layers separated. The burst pressure of scheme A_4 was 165.73 MPa and was increased by 7.88%. The results showed that the burst pressure of hydrogen storage vessels was improved by adopting the model with hoop and helical layers separated.
Figure 18 shows the stress distribution of the aluminum liner of scheme B_3. The maximum stress of the aluminum liner under working pressure was 291.41 MPa before autofrettage, which exceeded the yield strength of the material. Residual tensile stress was produced in composite layers after autofrettage. The residual compressive stress of the aluminum liner was caused by the restriction of the composite layers. The maximum stress of the aluminum liner under working pressure was reduced to 237.97 MPa, and the fatigue performance of the vessel was improved after autofrettage.

4.4. Comparison with the Cubic Spline Function Method

The comprehensive performance of the scheme B_3 was the best. For comparison, the cubic spline function method was used to analyze the hydrogen storage vessel of scheme B_3. Figure 19 shows the distribution curves of dome thickness predicted by the present method and the cubic spline function method. The cubic spline function method underestimated dome thickness at polar openings, and a sharp thickness peak occurred at about one bandwidth. The curve predicted by the present method was smoother and avoided a sharp thickness peak at about one bandwidth. The present method provided more reasonable and accurate results.
A finite-element model of the vessel was established by using the dome thickness predicted by the cubic spline function method. The burst pressure obtained by the finite-element analysis was 173.3 MPa, with little difference from the 172.9 MPa predicted by the present method. The radial displacement curves of the two methods were similar, as shown in Figure 20. The axial displacement curves of the two methods had no obvious difference before 113 MPa but separated after 113 MPa. The results showed that the pressure of matrix damage was 113 MPa. The axial displacement curve of the cubic spline function method increased more obviously than the curve of the present method. The axial displacement curve of the cubic spline function method did not present a burst point, which showed that the burst mode of the vessel predicted by the cubic spline function method was unsafe.
As shown in Figure 21, the fibers in the dome section were partially damaged under the burst pressure from the results of the present method, but not enough to cause structural damage. The cubic spline function method showed that the fiber damage in the dome section eventually formed a penetration area, resulting in structural damage. The cubic spline function method underestimated the actual thickness of polar openings, leading to a lower strength of composite layers. The finite-element results predicted by the present method were closer to the actual results.

5. Conclusions

In the present research, a method for predicting the thickness of the dome section was proposed, which involved the effect of fiber slippage and the redistribution of fiber tows. This prediction method was used to perform progressive damage analysis on composite hydrogen storage vessels. The analysis results were compared with the results of the cubic spline function method. The following conclusions can be deduced:
(1)
The maximum relative error of the predicted thickness was 4.19%, the mean absolute percentage error was 2.04%, and the root-mean-square error was 0.48 mm. The results verified that the predicted thickness of the present method was in good agreement with the actual measured thickness.
(2)
The results showed that the reaming process increased the burst pressure of the hydrogen storage vessel by 9.15%. The maximum stress of the aluminum liner under the working pressure was reduced by 19.15% after autofrettage, and the stress amplitude was reduced by 9.28%.
(3)
The stress concentration in the composite layers caused by fiber overhead was improved by changing the layer-up sequence as follows: the hoop and helical layers were separated, the winding angles of helical layers were changed from small to large, and the inner portion of the composite layers consisted of helical layers. The burst pressure of the vessel was increased by 7.88% compared with the traditional design sequence.
(4)
The burst pressures predicted by the two methods had little difference. However, the cubic spline function method underestimated the actual thickness of polar openings. The burst mode of the vessel predicted by the cubic spline function method was unsafe. The analysis results of the present method approached the actual case.

Author Contributions

Conceptualization, H.W. and S.F.; methodology, H.W. and S.F.; software, S.F.; validation, H.W., Y.C. and L.H.; formal analysis, Y.C.; investigation, H.W. and Y.C.; resources, L.H.; data curation, S.F.; writing—original draft preparation, H.W. and S.F.; writing—review and editing, Y.C. and L.H.; visualization, S.F.; supervision, L.H.; project administration, H.W.; funding acquisition, H.W. and L.H. All authors have read and agreed to the published version of the manuscript.

Funding

The authors would like to acknowledge the financial support from the National Natural Science Foundation Council of China (Grant No. 51775398, 51805392), the 111 Project (Grant No. B17034), the Program for Innovative Research Team in University of Education Ministry (Grant No. IRT_17R83), the Nature Science Foundation of Hubei Province (Grant No. 2018CFB595), and the Fundamental Research Funds for the Central Universities (Grant No. WUT: 2019IVB022, 2018III074GX, 2018III067GX).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Son, D.S.; Chang, S.H. Evaluation of modeling techniques for a type III hydrogen pressure vessel (70 MPa) made of an aluminum liner and a thick carbon/epoxy composite for fuel cell vehicles. Int. J. Hydrogen Energy 2012, 37, 2353–2369. [Google Scholar] [CrossRef]
  2. Wang, X.; Lin, L.; Lu, S.; Zhang, L.; Li, B.; Zhang, D.; Luo, Y. Evaluation of embedded buckypaper sensors in composite overwrappedped pressure vessels for progressive damage monitoring. Compos. Struct. 2022, 284, 115223. [Google Scholar] [CrossRef]
  3. Ma, L.; He, J.; Gu, Y.; Zhang, Z.; Yu, Z.; Zhou, A.; Tam, L.H.; Wu, C. Structure design of gfrp composite leaf spring: An experimental and finite element analysis. Polymers 2021, 13, 1193. [Google Scholar] [CrossRef] [PubMed]
  4. Eggers, F.; Almeida, J.H.S.; Lisbôa, T.V.; Amico, S.C. Creep and residual properties of filament-wound composite rings under radial compression in harsh environments. Polymers 2021, 23, 33. [Google Scholar] [CrossRef]
  5. Mishra, R.; Behera, B.K.; Mukherjee, S.; Petru, M.; Muller, M. Axial and radial compression behavior of composite rocket launcher developed by robotized filament winding: Simulation and experimental validation. Polymers 2021, 13, 517. [Google Scholar] [CrossRef] [PubMed]
  6. Boon, Y.D.; Joshi, S.C.; Bhudolia, S.K. Review: Filament winding and automated fiber placement with in situ consolidation for fiber reinforced thermoplastic polymer composites. Polymers 2021, 13, 1951. [Google Scholar] [CrossRef] [PubMed]
  7. Kashyzadeh, K.R.; Koloor, S.S.R.; Bidgoli, M.O.; Petrů, M.; Asfarjani, A.A. An optimum fatigue design of polymer composite compressed natural gas tank using hybrid finite element-response surface methods. Polymers 2021, 13, 483. [Google Scholar] [CrossRef]
  8. Jiao, W.; Wang, R.; Liu, W.; Yang, F. Dome thickness prediction of composite pressure vessels. Fuhe Cailiao Xuebao (Acta Mater. Compos. Sin.) 2010, 27, 116–121. [Google Scholar]
  9. Krikanov, A.A. Refined thickness of filament wound shells. Sci. Eng. Compos. Mater. 2002, 10, 241–248. [Google Scholar] [CrossRef]
  10. Hartung, R.F. Planar-wound filamentary pressure vessels. AIAA J. 1963, 1, 2842–2844. [Google Scholar] [CrossRef]
  11. Wang, R.; Jiao, W.; Liu, W.; Yang, F. A new method for predicting dome thickness of composite pressure vessels. J. Reinf. Plast. Compos. 2010, 29, 3345–3352. [Google Scholar] [CrossRef]
  12. Wang, R.; Jiao, W.; Liu, W.; Yang, F. Dome thickness prediction of composite pressure vessels by a cubic spline function and finite element analysis. Polym. Polym. Compos. 2011, 19, 227–234. [Google Scholar] [CrossRef]
  13. Leh, D.; Saffré, P.; Francescato, P.; Arrieux, R. Multi-sequence dome lay-up simulations for hydrogen hyper-bar composite pressure vessels. Compos. Part A Appl. Sci. Manuf. 2013, 52, 106–117. [Google Scholar] [CrossRef]
  14. Zu, L.; Xu, H.; Jia, X.; Zhang, Q.; Wang, H.; Zhang, B. Winding path design based on mandrel profile updates of composite pressure vessels. Compos. Struct. 2020, 235, 111766. [Google Scholar] [CrossRef]
  15. Zhang, Q.; Xu, H.; Jia, X.; Zu, L.; Cheng, S.; Wang, H. Design of a 70 MPa type IV hydrogen storage vessel using accurate modeling techniques for dome thickness prediction. Compos. Struct. 2020, 236, 1–10. [Google Scholar] [CrossRef]
  16. Che, J.-l.; Han, M.G.; Chang, S.H. Prediction of composite layer thickness for Type III hydrogen pressure vessel at the dome part. Compos. Struct. 2021, 271, 114177. [Google Scholar] [CrossRef]
  17. Akhtar, M.K.; Qureshi, W.A.; Jamshed, R.; Raza, M. Analytical modeling and numerical analysis of the effect of mosaic patterns on composite pressure vessels with dome. J. Reinf. Plast. Compos. 2021, 40, 621–631. [Google Scholar] [CrossRef]
  18. Han, Z.; Li, H.; Xu, X.; Wang, H.; Li, H.; Tao, J. Crushing characteristics of aluminum/CFRP/aluminum hybrid tubes prepared by spinning forming. Compos. Struct. 2020, 249, 112551. [Google Scholar] [CrossRef]
  19. Hu, Z.; Chen, M.; Pan, B. Simulation and burst validation of 70 MPa type IV hydrogen storage vessel with dome reinforcement. Int. J. Hydrogen Energy 2021, 46, 23779–23794. [Google Scholar] [CrossRef]
  20. Liao, B.-B.; Wang, D.-L.; Jia, L.-Y.; Zheng, J.-Y.; Gu, C.-H. Continuum damage modeling and progressive failure analysis of a Type III composite vessel by considering the effect of autofrettage. J. Zhejiang Univ. Sci. A 2019, 20, 36–49. [Google Scholar] [CrossRef]
  21. Kangal, S.; Kartav, O.; Tanoğlu, M.; Aktaş, E.; Artem, H.S. Investigation of interlayer hybridization effect on burst pressure performance of composite overwrapped pressure vessels with load-sharing metallic liner. J. Compos. Mater. 2020, 54, 961–980. [Google Scholar] [CrossRef]
  22. Lin, S.; Yang, L.; Xu, H.; Jia, X.; Yang, X.; Zu, L. Progressive damage analysis for multiscale modelling of composite pressure vessels based on Puck failure criterion. Compos. Struct. 2021, 255. [Google Scholar] [CrossRef]
  23. Alam, S.; Yandek, G.R.; Lee, R.C.; Mabry, J.M. Design and development of a filament wound composite overwrapped pressure vessel. Compos. Part C Open Access 2020, 2, 100045. [Google Scholar] [CrossRef]
  24. Sharma, P.; Bera, T.; Semwal, K.; Badhe, R.M.; Sharma, A.; Ramakumar, S.S.V.; Neogi, S. Theoretical analysis of design of filament wound type 3 composite cylinder for the storage of compressed hydrogen gas. Int. J. Hydrogen Energy 2020, 45, 25386–25397. [Google Scholar] [CrossRef]
  25. Hu, Z.; Chen, M.H.; Zu, L.; Jia, X.; Shen, A.; Yang, Q.; Xu, K. Investigation on failure behaviors of 70 MPa Type IV carbon fiber overwound hydrogen storage vessels. Compos. Struct. 2021, 259, 1–11. [Google Scholar] [CrossRef]
  26. Weerts, R.A.J.; Cousigné, O.; Kunze, K.; Geers, M.G.D.; Remmers, J.J.C. Novel material model to predict the residual strength of a composite overwrapped pressure vessel after impact. Int. J. Impact Eng. 2022, 160. [Google Scholar] [CrossRef]
  27. Nguyen, B.N.; Roh, H.S.; Merkel, D.R.; Simmons, K.L. A predictive modeling tool for damage analysis and design of hydrogen storage composite pressure vessels. Int. J. Hydrogen Energy 2021, 46, 20573–20585. [Google Scholar] [CrossRef]
  28. Ahmadi Jebeli, M.; Heidari-Rarani, M. Development of Abaqus WCM plugin for progressive failure analysis of type IV composite pressure vessels based on Puck failure criterion. Eng. Fail. Anal. 2022, 131, 105851. [Google Scholar] [CrossRef]
  29. MathWorks. Available online: https://ww2.mathworks.cn/en/ (accessed on 1 December 2021).
  30. Gramoll, K.; Onoda, J.; Namiki, F. Dome thcikness of filament wound pressure vessels. Trans. Japan Soc. Aeronaut. Space Sci. 1990, 33, 69–79. [Google Scholar]
  31. Peters, S.T. Composite Filament Winding; ASM International: Cleveland, OH, USA, 2011; pp. 19–28. ISBN 978-1-61503-722-3. [Google Scholar]
  32. Puck, A.; Schürmann, H. Failure analysis of FRP laminates by means of physically based phenomenological models. Compos. Sci. Technol. 2002, 62, 1633–1662. [Google Scholar] [CrossRef]
  33. Camanho, P.P.; Matthews, F.L. A progressive damage model for mechanically fastened joints in composite laminates. J. Compos. Mater. 1999, 33, 2248–2280. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of fiber slippage at different winding angles.
Figure 1. Schematic diagram of fiber slippage at different winding angles.
Polymers 14 00902 g001
Figure 2. Comparison of analytical results of the arc length corresponding to the radius of parallel circle between two methods.
Figure 2. Comparison of analytical results of the arc length corresponding to the radius of parallel circle between two methods.
Polymers 14 00902 g002
Figure 3. Comparison of analytical and measurement results for dome thickness distribution within two bandwidths: (a) Type A vessel, (b) Type B vessel, and (c) Type C vessel.
Figure 3. Comparison of analytical and measurement results for dome thickness distribution within two bandwidths: (a) Type A vessel, (b) Type B vessel, and (c) Type C vessel.
Polymers 14 00902 g003
Figure 4. Analytical results of the variations of winding angle and dome thickness.
Figure 4. Analytical results of the variations of winding angle and dome thickness.
Polymers 14 00902 g004
Figure 5. Layer-up design schemes used in the finite-element modeling: (a) Group A, (b) Group B, and (c) Group C.
Figure 5. Layer-up design schemes used in the finite-element modeling: (a) Group A, (b) Group B, and (c) Group C.
Polymers 14 00902 g005
Figure 6. Cross-section of lining structure.
Figure 6. Cross-section of lining structure.
Polymers 14 00902 g006
Figure 7. Finite-element model of the composite hydrogen storage vessel.
Figure 7. Finite-element model of the composite hydrogen storage vessel.
Polymers 14 00902 g007
Figure 8. Stress-strain curve of the aluminum alloy.
Figure 8. Stress-strain curve of the aluminum alloy.
Polymers 14 00902 g008
Figure 9. Flow chart of progressive damage process.
Figure 9. Flow chart of progressive damage process.
Polymers 14 00902 g009
Figure 10. Loading history for FEM of hydrogen storage composite vessels.
Figure 10. Loading history for FEM of hydrogen storage composite vessels.
Polymers 14 00902 g010
Figure 11. Comparison of numerical results of two burst modes: (a) safe mode, (b) unsafe mode.
Figure 11. Comparison of numerical results of two burst modes: (a) safe mode, (b) unsafe mode.
Polymers 14 00902 g011
Figure 12. Numerical results of the burst pressure and the burst mode.
Figure 12. Numerical results of the burst pressure and the burst mode.
Polymers 14 00902 g012
Figure 13. Comparison of numerical results of the stress in fiber direction at the cylindrical section: (a) Schemes A_4 and C_4, (b) Scheme C_1.
Figure 13. Comparison of numerical results of the stress in fiber direction at the cylindrical section: (a) Schemes A_4 and C_4, (b) Scheme C_1.
Polymers 14 00902 g013
Figure 14. Comparison of numerical results of load-displacement curve between schemes B_3 and B_5.
Figure 14. Comparison of numerical results of load-displacement curve between schemes B_3 and B_5.
Polymers 14 00902 g014
Figure 15. Comparison of numerical results of fiber damage at the burst pressure: (a) Scheme C_3, (b) Scheme C_4, (c) Scheme C_5, (d) Scheme B_3, and (e) Scheme B_5.
Figure 15. Comparison of numerical results of fiber damage at the burst pressure: (a) Scheme C_3, (b) Scheme C_4, (c) Scheme C_5, (d) Scheme B_3, and (e) Scheme B_5.
Polymers 14 00902 g015
Figure 16. Numerical results of progressive damage of the composite layers in the composite hydrogen storage vessels: (a) fiber tensile, (b) matrix tensile, and (c) matrix compression.
Figure 16. Numerical results of progressive damage of the composite layers in the composite hydrogen storage vessels: (a) fiber tensile, (b) matrix tensile, and (c) matrix compression.
Polymers 14 00902 g016
Figure 17. Numerical results of the maximum stress of the aluminum liner: (a) under different pressure; (b) stress reduction and stress amplitude.
Figure 17. Numerical results of the maximum stress of the aluminum liner: (a) under different pressure; (b) stress reduction and stress amplitude.
Polymers 14 00902 g017
Figure 18. Numerical results of loading history of the aluminum liner at (a) service pressure before autofrettage, (b) zero pressure after autofrettage, and (c) service pressure after autofrettage.
Figure 18. Numerical results of loading history of the aluminum liner at (a) service pressure before autofrettage, (b) zero pressure after autofrettage, and (c) service pressure after autofrettage.
Polymers 14 00902 g018
Figure 19. Comparison of analytical results of dome thickness between two methods.
Figure 19. Comparison of analytical results of dome thickness between two methods.
Polymers 14 00902 g019
Figure 20. Comparison of numerical results of load-displacement curves between two methods.
Figure 20. Comparison of numerical results of load-displacement curves between two methods.
Polymers 14 00902 g020
Figure 21. Comparison of numerical result of the fiber damage at the burst pressure: (a) present method, (b) the cubic spline function method.
Figure 21. Comparison of numerical result of the fiber damage at the burst pressure: (a) present method, (b) the cubic spline function method.
Polymers 14 00902 g021
Table 1. The parameters of three composite vessels.
Table 1. The parameters of three composite vessels.
Vessel’s TypeDome Length h (mm)Fiber Band b (mm)Thickness of Fiber Tows tp (mm)Polar Radius r0 (mm)Radius of Cylinder Section R (mm)
a805.390.7659.5132
b2135.390.76560132
c21312.50.48110372.5
Table 2. Relative error between the analytical results and the measurement results.
Table 2. Relative error between the analytical results and the measurement results.
The Present MethodThe Double Formula MethodThe Cubic
Spline Function Method
r (mm)Actual Thickness (mm)Predicted Value (mm)Relative
Error (%)
Predicted Value (mm)Relative
Error (%)
Predicted Value (mm)Relative
Error (%)
9.508.508.25−2.940-5.31-
11.2022.0021.12−4.0020.69−5.9521.77−1.05
11.6023.5022.96−2.2922.67−3.5324.052.34
12.5026.0025.69−1.1226.271.0427.345.15
13.5027.3026.77−1.9429.427.7728.424.10
14.5026.6026.19−1.5731.9820.2327.282.56
15.4024.8024.63−0.6924.23−2.3025.101.21
16.3021.5022.404.1919.73−8.2322.132.93
17.5018.7018.931.2316.54−11.5518.00−3.74
18.5615.8015.961.0114.69−7.0315.00−5.06
19.5013.8013.810.0713.45−2.5413.38−3.04
20.3013.0012.63−2.8512.59−3.1512.45−4.23
Table 3. Material properties of aluminum alloy 6061-T6.
Table 3. Material properties of aluminum alloy 6061-T6.
PropertyE (GPa)μσs (MPa)σb (MPa)δ (%)ρ (kg/m3)
6061-T674.12 GPa0.2834028111.572700
Table 4. Material properties of the T700/epoxy composite.
Table 4. Material properties of the T700/epoxy composite.
PropertyE11 *
(GPa)
E22 *
(GPa)
E33 *
(GPa)
G12
(MPa)
G23
(MPa)
G13
(MPa)
μ12μ23μ13
T700/epoxy
composite
1347.427.423710479037100.280.30.28
* Direction 1—fiber direction, Direction 2—transverse direction, Direction 3—normal to layer direction.
Table 5. Material properties of the T700/epoxy composite.
Table 5. Material properties of the T700/epoxy composite.
PropertyXT
(MPa)
XC
(MPa)
YT
(MPa)
YC
(MPa)
S12 (MPa)S23 (MPa)S13 (MPa)ρ
(kg/m3)
T700/epoxy
composite
23001250741805090501680
Direction 1—fiber direction, Direction 2—transverse direction, Direction 3—normal to layer direction.
Table 6. Parameters of Puck criterion.
Table 6. Parameters of Puck criterion.
Parameter R t   ( MPa ) R c   ( MPa ) E   ( GPa ) E f   ( GPa ) ν ν f  
Valve230012502301340.280.2
R A t (MPa) R A (MPa) R c (MPa) p t , c p t p c
74501800.30.30.35
Table 7. Stiffness-degradation criteria of composite materials.
Table 7. Stiffness-degradation criteria of composite materials.
Failure ModesStiffness-Degradation Criteria
Fiber tensile failureE11 = 0.07E11, E22 = E33 = 0.07E22, υ12 = υ13 = 0.07υ12,
υ23 = 0.07υ23, G12 = G13 = 0.07G12, G23 = 0.07G23
Fiber compression failureE11 = 0.14E11, E22 = E33 = 0.14E22, υ12 = υ13 = 0.14υ12,
υ23 = 0.14υ23, G12 = G13 = 0.14G12, G23 = 0.14G23
Matrix tensile failureE22 = 0.2E22, G12 = 0.2G12, G23 = 0.2G23
Matrix compression failureE22 = 0.4E22, G12 = 0.4G12, G23 = 0.4G23
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, H.; Fu, S.; Chen, Y.; Hua, L. Thickness-Prediction Method Involving Tow Redistribution for the Dome of Composite Hydrogen Storage Vessels. Polymers 2022, 14, 902. https://doi.org/10.3390/polym14050902

AMA Style

Wang H, Fu S, Chen Y, Hua L. Thickness-Prediction Method Involving Tow Redistribution for the Dome of Composite Hydrogen Storage Vessels. Polymers. 2022; 14(5):902. https://doi.org/10.3390/polym14050902

Chicago/Turabian Style

Wang, Hui, Shuang Fu, Yizhe Chen, and Lin Hua. 2022. "Thickness-Prediction Method Involving Tow Redistribution for the Dome of Composite Hydrogen Storage Vessels" Polymers 14, no. 5: 902. https://doi.org/10.3390/polym14050902

APA Style

Wang, H., Fu, S., Chen, Y., & Hua, L. (2022). Thickness-Prediction Method Involving Tow Redistribution for the Dome of Composite Hydrogen Storage Vessels. Polymers, 14(5), 902. https://doi.org/10.3390/polym14050902

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop