Next Article in Journal
Design and Characterization of the Surface Porous UHMWPE Composite Reinforced by Graphene Oxide
Next Article in Special Issue
Fabrication of High-Quality Straight-Line Polymer Composite Frame with Different Radius Parts Using Fiber Winding Process
Previous Article in Journal
Bio-Packaging Material Impact on Blueberries Quality Attributes under Transport and Marketing Conditions
Previous Article in Special Issue
An Extended Thickness-Dependent Moisture Absorption Model for Unidirectional Carbon/Epoxy Composites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Optimum Fatigue Design of Polymer Composite Compressed Natural Gas Tank Using Hybrid Finite Element-Response Surface Methods

by
Kazem Reza Kashyzadeh
1,2,*,
Seyed Saeid Rahimian Koloor
3,4,
Mostafa Omidi Bidgoli
5,
Michal Petrů
3 and
Alireza Amiri Asfarjani
6
1
Department of Transport, Academy of Engineering, Peoples’ Friendship University of Russia (RUDN University), 6 Miklukho-Maklaya Street, 117198 Moscow, Russia
2
Department of Mechanical and Instrumental Engineering, Academy of Engineering, Peoples’ Friendship University of Russia (RUDN University), 6 Miklukho-Maklaya Street, 117198 Moscow, Russia
3
Institute for Nanomaterials, Advanced Technologies and Innovation, Technical University of Liberec, Studentska 2, 461 17 Liberec, Czech Republic
4
School of Mechanical Engineering, Faculty of Engineering, University Technology Malaysia (UTM), Johor Bahru 81310, Malaysia
5
Department of Mechanical Engineering, Islamic Azad University, Badroud Branch, Badroud, Isfahan 8764136545, Iran
6
Department of Mechanical Engineering, Shahrood University of Technology, Shahrood 3619995161, Iran
*
Author to whom correspondence should be addressed.
Submission received: 8 December 2020 / Revised: 20 January 2021 / Accepted: 28 January 2021 / Published: 3 February 2021
(This article belongs to the Special Issue Damage Mechanics of Polymer Composites)

Abstract

:
The main purpose of this research is to design a high-fatigue performance hoop wrapped compressed natural gas (CNG) composite cylinder. To this end, an optimization algorithm was presented as a combination of finite element simulation (FES) and response surface analysis (RSA). The geometrical model was prepared as a variable wall-thickness following the experimental measurements. Next, transient dynamic analysis was performed subjected to the refueling process, including the minimum and maximum internal pressures of 20 and 200 bar, respectively. The time histories of stress tensor components were extracted in the critical region. Furthermore, RSA was utilized to investigate the interaction effects of various polymer composite shell manufacturing process parameters (thickness and fiber angle) on the fatigue life of polymer composite CNG pressure tank (type-4). In the optimization procedure, four parameters including wall-thickness of the composite shell in three different sections of the CNG tank and fiber angle were considered as input variables. In addition, the maximum principal stress of the component was considered as the objective function. Eventually, the fatigue life of the polymer composite tank was calculated using stress-based failure criterion. The results indicated that the proposed new design (applying optimal parameters) leads to improve the fatigue life of the polymer composite tank with polyethylene liner about 2.4 times in comparison with the initial design.

1. Introduction

Compressed natural gas (CNG) is used as a good alternative to petrol fuel in land transportation (especially, automotive) in many countries. The CNG fuel tank should be manufactured with a high safety factor and pass stringent standards to get the production license because this super-critical component deals directly with human life. On the other hand, hazards from the failure of CNG automotive cylinders can cause serious damage to the vehicle and even the death of the driver or occupants. CNG fuel tanks were first made of metal. Afterward, the usage of composite materials in the industry has led to the discovery of type-3 composite tanks with the aim of weight reduction. This type of tank consists of a metal liner and a polymer composite shell reinforced with carbon or glass fiber. However, due to the bursting of the metal layer due to cyclic internal pressures (fatigue phenomenon), the modified design of the composite tank uses a non-metallic liner such as polyethylene, nylon, polyether ether ketone, polyethylene terephthalate, and polypropylene (finding the type-4 composite tank). One of the most important advantages of this particular type of CNG fuel tank compared to other types is that fatigue failure appears as the leakage of the gas and as a result prevents the tank from bursting suddenly. It is noteworthy that, given the high sensitivity of this component, little research has been done to investigate the behavior of the polymer composite CNG fuel tank under fatigue loads and various working conditions such as in the vicinity of the fire. However, the details of some of these studies are discussed below.
The behavior of a type-3 composite tank considering the low alloy steel liner has been investigated under static loading (maximum internal pressure) [1]. Their research showed that the use of composite tanks compared to the metal tank has more corrosion resistance and the least stress due to fuel filling. In addition, they stated that the weight of the composite tank is much less than the weight of the metal tank. Moreover, finite element simulation was used to study the failure of the composite type-4 CNG tank [2]. In this research, a two-liter tank has been used, which is very small compared to the tanks used in the automotive industry. However, no experimental test was performed and the published results in other studies were used to validate the presented finite element (FE) model and use of the WCM plugin. Finally, various static failure criteria of composite materials including Tsai–Hill, Tsai–Wu, and Hashin have been applied and burst pressure has been calculated. Finally, they reported that the damage rate and breakdown in the spiral layers of the tank are higher than the peripheral layers; also, the regularity of the layer destruction is different.
The cyclic lifetime of the CNG composite tank with the surface defects (notch with the same width and depth, which is considered to be equal to 0.2 wall-thickness) has been predicted by utilizing the strain-life criterion [3]. Kim and Choi have performed a risk analysis of the CNG composite tank to determine the period time of visual inspection [4]. The findings of this research stated that the status of gas leakage, internal pressure, and temperature should be checked every three years.
Nouri et al. have estimated the static strength of the CNG polymer composite cylinder tank subjected to the maximum level of internal pressure [5]. They used different composite failure theories under static loading. By comparing the obtained results with laboratory data, they stated that Hoffmann benchmarking gave the most accurate prediction of static strength. However, due to its complex calculations, it is recommended to use Tsai–Wu, Tsai–Hill, and Christensen criteria to reduce computational costs. Also, the results of this study showed that the shear stress in the XY plane is the main cause of tank failure. Moreover, they reported that the separation between the composite shell and the lower aluminum flange occurs.
Seyedi et al. have analyzed and optimized the type-4 composite fuel tanks (without metal) [6]. In the first step, the stress analysis of a metallic tank (made of carbon-manganese) was performed in finite element software. Then, the prediction of fracture zone was validated by comparison with experimental work. Next, the polymer composite tank was examined. One of the weaknesses of this study was considering the the whole wall of the tank to be a constant thickness. In addition, the geometric modeling was based on the sizes in the standard, which will be different from the actual measurements, due to the different stages of manufacturing. Also, the arrangement of the layers is considered to be 45°. Finally, to optimize the strength and weight of the tank, four variables including the thickness of the liner and the composite shell, and the pressure inside the tank are considered. But the fiber angle is assumed to be constant in the optimization process. Also, to optimize, no special technique has been used, and only for each of the parameters, a range including low and high limits and interval has been selected and the results were compared with each other by trial and error. In this study, they showed that by maintaining the static strength of the tank, they can reduce its weight by more than half.
Reynaldo et al. studied the performance of a type-4 pressure vessel as a fuel tank to store the liquid ammonia [7]. They analyzed different polymer materials that are compatible with this type of fuel (e.g., high-density polyethylene, nylon, polyether ether ketone, polyethylene terephthalate, and polypropylene). Also, two fiber-reinforced polymers including carbon and glass were investigated. They reported that the best structure for this tank is the combination of CFRP and PET as a liner. Of course, this result is obtained only concerning the impact condition. In this research, it is assumed that the designed tank has a fixed wall-thickness and the layout was [90°, θ], in which θ = 10°, 30°, 50°, and 70°. In this area, other studies have been conducted to design, and perform stress analysis, failure analysis, structural stability, and optimization to store hydrogen as a vehicle fuel [8,9]. In addition, the safety analysis of the instantaneous release of CNG from a tank was investigated [10]. Altuwair and Khan reported that the probability of the damage increases with decreasing mechanical energy coefficient and as a consequence, the explosion rate rises.
Tschirschwitz et al. have experimentally studied the hazards from the failure of CNG automotive cylinders in the fire [11]. In this research, various types of polymer composite tanks (type-3 and type-4) have been tested. An assessment of fire safety risks was done for the CNG bus systems [12] and the mean fire fatality risk was obtained to be about 0.23 per 100-million miles. Yersak et al. have studied the thermal gradient and the temperature of the polymer composite tank due to the fast process of refilling [13]. The experimental results showed that for a 18.2 s refill, there is a 23.9 °C m−1 longitudinal thermal gradient and a 125.6 °C gas temperature. Moreover, Dicken and Merida have experimentally evaluated the influence of filling time on temperature distribution [14]. They reported that the final gas temperature decreases by raising the time of the filling process. In this regard, some studies have been performed using CFD simulation [15,16] and experimental setup [17,18].
Despite all the results of the researchers’ studies, the automotive industry is still facing the explosion of CNG fuel tanks, causing personal and financial damage to individuals. One of the most important reasons for the inconsistency of research results with reality is the use of a simplified model (for example, a perfectly symmetrical geometric model, with constant wall-thickness, etc.). Moreover, ignoring the actual performance of the tank during its service life and not performing analysis following reality (using static and quasi-static analysis instead of the dynamic transient with time under the pretext of cost and a long time to solve problems). Therefore, the response of optimization analyzes based on these principles cannot significantly reduce the occurrence of these events in the world. In the present research, the cyclic behavior of a type-4 CNG fuel tank considering variable wall-thickness of the polymer composite shell was investigated using finite element simulation. Moreover, for the first time, the response surface method as a well-known data mining technique was used to optimize the manufacturing process parameters including fiber angle and wall-thickness in different parts of the structure.

2. Methodology

As it is clear from the literature review, several parameters are important and effective in analyzing the stress of a CNG polymer composite tank and examining its strength and efficiency. One of them is the thickness of the composite shell, which varies in the cylindrical fuel tank. Its geometry includes different parts such as the cylinder in the main body, the hemisphere in the head of the tank, and also the connection of these two geometries to each other. In many studies done in the past, this parameter was considered constant and in the present research, in addition to modeling, it varies according to reality; this parameter was optimized in different parts of the CNG polymer composite tank. Another important parameter is the fiber angle and the layout of the composite, which the vast majority of researchers have considered in their studies and tried to increase the strength of the tank by choosing the most appropriate fiber angle. In this study, this factor was considered as a variable in the optimization process. Moreover, one of the advantages of this research compared to the others is the objective function in the optimization process. As in previous research, weight reduction or increase in static strength of the fuel tank has been considered. But, in this study, the service life of the fuel tank under fatigue loading due to the refueling process was considered. The flowchart of the process of getting the longest cyclic lifetime of the polymer composite tank by optimizing various composite shell manufacturing process parameters including fiber angle and wall-thickness is summarized in Figure 1.
In the proposed algorithm, firstly the geometric model of the tank is provided as a parametric model using the empirical measurements [5]. Then, finite element simulation of the problem is done using the geometric model, meshing process, applying the boundary conditions and load history, and defining the material properties [5]. Based on the finite element (FE) results, the critical area is determined using maximum principal stress [19]. Next, the optimization process was performed using the response surface method considering the most effective parameter [20,21,22]. After reaching the optimum design that results in the least principal stress, the fatigue life of the component was calculated by employing one of the well-known stress-based failure criteria for polymer composite materials.

3. Finite Element Simulation

In the present research, a case study was conducted on the 34-L CNG tank for application in light passenger vehicles in an Asian country (Iran). The dimensions and size of the CNG tank are shown in Figure 2. The CATIA software (Dassault Systèmes SE, Vélizy-Villacoublay, France) was used to model the geometry of the tank as having a variable wall-thickness following the experimental measurements. As shown in Figure 3, the liner thickness is uniform through all parts of the tank (7 mm). However, composite shell thickness is variable (minimum thickness is related to the middle part of the tank and it is maximum at the beginning and end of the tank).
This model includes shell (polymer composite) and volumetric (liner and flanges) parts [5]. Hence, two types of elements (shell and solid) were used to mesh the component in the Altair Hypermesh (Figure 4). To increase the accuracy of the solution and reliable response, a 2nd order element with mid-side nodes was used for meshing composite layers and liner. Moreover, the adhesion between different parts (composite shell, liner, and flanges) was considered 100% [5], and the bonded type of contact was defined. However, for meshing aluminum flanges, the first order solid element was used. Also, a valid technique was applied to couple the shell-to-solid elements [23]. The convergence mesh was performed to select an acceptable element size; then further mesh refinement does not significantly affect the resulting maximum principal stress of the tank [19,23,24].
This model of the CNG tank (Type-4) consists of three main parts called: flanges, liner, and polymer composite shell. Therefore, three different materials (aluminum alloy, polyethylene, and glass/epoxy composite, respectively) were defined in the ANSYS WORKBENCH (Ansys Inc., Canonsburg, PA, USA). Mechanical properties of orthotropic and homogeneous materials are presented in Table 1 and Table 2, respectively.
In the manufacturing process of this tank in the industry, the fiber of Glass 2400, Epoxy resin UN3082, and UN 2735 hardener are used, considering a fiber volume fraction of 50%. In addition, the thickness of the layers is constant and equal to 0.55 mm. The static properties required for simulation of each composite layer can be calculated using three-point bending results such as stiffness, Young’s modulus of bending, flexural rigidity, load, and deflection at maximum and breakpoint [25]. Moreover, the fiber angle of ±17° was considered as the tank manufacturer’s recommendation.
In the present FE model, there is no separation between the different parts of the tank during its service life. The internal pressure of 20–200 bar (given equal time for the fill and empty processes) was applied on the full model of the CNG tank and both sides of the tank including flanges that were fixed in all directions. The loading diagram is displayed in Figure 5.

4. Response Surface Method

Various designs of experiment techniques (DOEs) are used to investigate the effect of different parameters on the response when the number of samples is too large and we do not want to examine all of them. In other words, by performing some experiments with predefined settings, the effect of different parameters on the response can be achieved with high accuracy. In this regard, the Taguchi approach represents the minimum number of experiments [20,21,22,24]. But this method can only examine the effect of each parameter on the response separately. However, the response surface method can investigate the effects of changes in two parameters simultaneously on the response. In other words, one of the advantages of using the response surface method (RSM) compared to other data mining methods and design of experiment (DOE) techniques is examining the effect of parameters on the response by considering the interaction between them. Moreover, the output of this method is three-dimensional diagrams including changes of two independent parameters and the response that can be used to optimize the problem. Therefore, this method was used to more accurately investigate the effect of different parameters of the composite tank manufacturing process on the service time of the tank and finally the optimization of the parameters to increase the efficiency of the tank. To this end, four parameters including wall-thickness of polymer composite shell in three different sections of CNG tank and fiber angle were considered as input variables. In addition, five levels were considered for each input variable. The maximum principal stress of the component was considered as an output. A schematic of the used algorithm including the input variables and response is illustrated in Figure 6. Moreover, the variables and different levels considered as inputs for the response surface analysis are given in Table 3. In this study, the linear effects of each parameter and the interaction between them were considered as follows:
A, T1, T2, T3, AT1, AT2, AT3, T1T2, T1T3, and T2T3.
After analyzing the data by RSM [20,21,22] to investigate the effect of parameters (wall-thickness of composite shell in three different sections of the CNG tank and the fiber angle) on the response (maximum principle stress), the optimization tool was used to reduce stress in the critical region of the structure and the optimum values for the parameters were extracted. Since there is a direct relationship between the fatigue life of the component and stress, the optimum tank is expected to have a maximum fatigue life.

5. Results and Discussion

Explicit transient dynamics analysis was performed and the stress tensor components were extracted in terms of time. The small time increments were considered for the convergence of the solution, then, the time history of principal stress was obtained as the FE response [26,27,28] for all cases proposed by RSM. Next, the maximum value in each history was considered as the output in the response surface method (Table 4). To facilitate the analysis by the response surface method, the response values are always considered positive [29]. Also, the principal stress contour in the CNG tank concerning the values of parameters provided by the manufacturer (initial design) is demonstrated in Figure 7. According to the results of the FE analysis, the highest absolute principle stress (S = 206.32 MPa) is related to the joint area of the composite shell to the flange, which is consistent with the results presented by Nouri et al. [5]. Then, the bonding area between the composite layers in the two parts of the main body and the spherical curvature is critical. Therefore, these two regions in the composite CNG tank are identified as critical areas that are apt to fail and future analysis will be carried out to estimate fatigue life only in the critical elements of these two areas.
The design of the CNG tanks is based on high-cycle fatigue (HCF) life. In this type of fatigue, the deformation occurs in the elastic region. Hence, stress-based criteria are considered. Finally, it can be stated that the value of stress created in the component is directly related to its fatigue life and the damage rate. As the stress level decreases, the fatigue life of the component increases. The stress created in the composite shell depends on the wall-thickness. On the other hand, it is clear that by changing the angle of the fibers, the directions of the principal axes also change and the stress created is also affected. Figure 8 presents the effect of each parameter on the stress created in the critical area of the composite shell.
From Figure 8, increasing the fiber angle reduces the stress level in the composite shell. Therefore, it seems better to use the maximum twisting angle of fibers in the manufacturing process. Moreover, simulations showed that increasing the wall-thickness in different parts has different results on the stress level. Increasing the wall-thickness of the middle section leads to a decrease in the level of stress, and conversely, by increasing the wall-thickness of the other parts, the level of stress increases. Therefore, to achieve an optimal design in composite CNG fuel tanks, the thickness of the composite shell should be considered and an intermediate value should be determined for them so that behavioral improvements in one area do not lead to damage to another.
Next, the results of the RSM analysis are illustrated in Figure 9. It should be noted that in these analyses, the constant parameters were assumed to be the average level (L3) [20,21,22].
Moreover, a Pareto chart of the standardized effects was used to determine the most effective parameter in the system response. As shown in Figure 10, the parameter of fiber angle has the greatest effect on the performance of composite CNG tanks; also among the thickness parameters, wall-thickness in the middle section is most important. The least effective parameter on the stress results is the wall-thickness of the cylindrical part of the tank, because the main section has a more uniform structure than the other sections as well as low-stress concentration. These results are completely consistent with the simulation results described above (Figure 8). Eventually, the most important results extracted from the contours of Figure 9 are:
  • The maximum principal stress always decreases by increasing the fiber angle. Therefore, it is recommended to consider the maximum level of fiber angle to achieve the maximum fatigue lifetime of the polymer composite CNG tank.
  • In order to have the minimum principal stress in the polymer composite CNG tank, it is necessary to consider the maximum and minimum level settings for parameters of T2 and T1, respectively. On the other hand, it is better to make the composite shell thickness as uniform as possible in the two curved sections of the CNG tank.
In the following, the optimization process was performed using RSM and the values were compared with the baseline design data in Table 5. The findings of this research reveal that the stress level in the optimum composite CNG tank decreases by 10.12%.
In the last step, the fatigue life analysis in the time domain was carried out by using coupling results of stress analyses under cyclic loading (internal pressure of 20 to 200 Bar) and Fawaz–Ellyin criterion [30] in the NCODE DESIGN LIFE software [28]. The obtained results indicate that the fatigue life of the optimized CNG tank is 218,687 cycles. However, the cyclic lifetime for the initial design was 89,841 cycles. In other words, the service life of the variable wall-thickness composite CNG tank can be extended about 2.4 times by utilizing the proposed optimum design.

6. Conclusions

In the present research, the authors have attempted to extend the service life of the CNG fuel tank. To achieve this purpose, the study was performed on a real tank sample with variable wall-thickness in different parts. A hybrid finite element-response surface method was presented to predict the strength of polymer composite CNG tanks via wall thickness and fiber angle. Eventually, these two parameters were optimized using the RSM, and the fatigue life of the optimal tank was compared with the original tank. The following practical hints may be deduced from the presented results in this study:
  • The FE results indicate that the maximum principal stress created in the initial design of the CNG composite tank (A = 17, T1 = 11 mm, T2 = 4.4 mm, and T3 = 3.3 mm) under static load is about 206 MPa. Also, the fatigue life of the tank under variations of internal pressure (20–200 bar) is 89,841 cycles. It means that the number of refueling times without the risk of tank failure is equal to 89,841.
  • The Pareto chart as a result of response surface analysis reveals that the most and least effective parameters on the stress level of the CNG composite tank are fiber angle and wall-thickness of the main section, respectively. Therefore, the authors strongly recommend to the manufacturers of Type-4 CNG fuel tanks to be more careful about the twist angle of the fibers in the manufacturing process to strengthen this super-critical component.
  • The RSM results show that the maximum principal stress always decreases by increasing fiber angle. Therefore, it is recommended to consider the maximum level of the fiber angle to achieve the maximum fatigue lifetime of the composite CNG tank. Moreover, to have the minimum principal stress in the polymer composite CNG tank, it is necessary to consider the maximum and minimum level settings for parameters T2 and T1, respectively. On the other hand, it is better to make the thickness of the polymer composite shell as uniform as possible in the two curved sections of the CNG tank. In other words, to achieve a longer service life, it is necessary to try to build tanks with invariable thickness.
  • Based on the optimization algorithm presented, the optimal parameters obtained for the tank are A = 30, T1 = 8.8 mm, T2 = 7.7 mm, and T3 = 2.2 mm. Also, the maximum principal stress and fatigue life related to the optimized tank are obtained about 185 MPa and 218,678 cycles, respectively.
  • The findings of this research state that the maximum principal stress decreases by up to 10% and after that the fatigue life of the fuel tank increases about 2.4 times using optimal values of manufacturing process parameters. After confirmation of the results by pre-production tests, it can drastically reduce the financial losses and irreparable damages due to the fuel tank explosion.

Author Contributions

Conceptualization, K.R.K., and M.O.B.; Data curation, K.R.K.; Formal analysis, K.R.K., M.O.B., S.S.R.K., M.P., and A.A.A.; Funding acquisition, S.S.R.K., and M.P.; Investigation, K.R.K., M.O.B., S.S.R.K., M.P., and A.A.A.; Methodology, K.R.K.; Project administration, K.R.K., S.S.R.K., and M.P.; Resources, S.S.R.K., and M.P.; Software, K.R.K., M.O.B., A.A.A., and S.S.R.K.; Supervision, K.R.K., and S.S.R.K.; Validation, K.R.K., M.O.B., S.S.R.K., M.P., and A.A.A.; Visualization, K.R.K., M.O.B., S.S.R.K., M.P., and A.A.A.; Writing—original draft, K.R.K., and A.A.A.; Writing—review & editing, K.R.K., and S.S.R.K. All authors have read and agreed to the published version of the manuscript.

Funding

The research was supported by the Ministry of Education, Youth, and Sports of the Czech Republic and the European Union (European Structural and Investment Funds Operational Program Research, Development, and Education) in the framework of the project “Modular platform for autonomous chassis of specialized electric vehicles for freight and equipment transportation”, Reg. No. CZ.02.1.01/0.0/0.0/16_025/0007293.

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 that they have no conflict of interest.

References

  1. Sawant, T.; Purwant, N.K.; Kulkarni, S.; Karajagi, P. Design & Stress Analysis of a Hoop Wrapped CNG Composite Vessel with an SAE-4135 Low Alloy Steel Liner. Int. J. Eng. Res. 2016, 338–343. [Google Scholar]
  2. Heaidari-Rarani, M.; Ahmadi Jebeli, M. Finite element modeling of failure in IV type composite pressure vessel using WCM plug-in in ABAQUS software. Modares Mech. Eng. 2018, 18, 191–200. [Google Scholar]
  3. Agrawal, A.; Kumar, S. Fatigue life prediction of a Hoop-Wrapped composite CNG cylinder containing surface flow. Int. J. Emerg. Technol. Adv. Eng. 2014, 4, 790–796. [Google Scholar]
  4. Kim, E.S.; Choi, S.K. Risk analysis of CNG composite pressure vessel via computer-aided method and fractography. Eng. Fail. Anal. 2013, 27, 84–98. [Google Scholar] [CrossRef]
  5. Nouri, M.; Ashenai-Ghasemi, F.; Rahimi-Sherbaf, G.; Kashyzadeh, K.R. Experimental and Numerical Study of the Static Performance of a Hoop-Wrapped CNG Composite Cylinder Considering Its Variable Wall Thickness and Polymer Liner. Mech. Compos. Mater. 2020, 56, 339–352. [Google Scholar] [CrossRef]
  6. Seyedi, S.M.; Naddaf Oskouei, A.; Sayah Badkhor, M. Experimental, numerical and Optimization study of Composite Tanks with Non-Metallic Primer (CNG Fourth Type). Modares Mech. Eng. 2020, 20, 1789–1800. [Google Scholar]
  7. Reynaldo, A.; Pramono, H.S.; Santosa, S.P.; Aziz, M. Finite Element Analysis of Liquefied Ammonia Tank for Mobility Vehicles Employing Polymers and Composites. Energies 2020, 13, 5312. [Google Scholar] [CrossRef]
  8. Sapre, S.; Pareek, K.; Vyas, M. Investigation of structural stability of type IV compressed hydrogen storage tank during refueling of fuel cell vehicle. Energy Storage 2020, 2, 150. [Google Scholar] [CrossRef]
  9. Ayakdaş, O.; Aydın, L.; Savran, M.; Küçükdoğan, N.; Öztürk, S. Optimal design of the type III hydrogen storage tank for different carbon/epoxy materials by modified differential evolution method. Res. Eng. Struct. Mat. 2019, 5, 189–201. [Google Scholar]
  10. Altuwair, I.A.; Khan, F. Safety analysis of instantaneous release of compressed natural gas from a cylinder. J. Loss Prev Process Ind. 2020, 68, 104284. [Google Scholar] [CrossRef]
  11. Tschirschwitz, R.; Krentel, D.; Kluge, M.; Askar, E.; Habib, K.; Kohlhoff, H.; Krüger, S.; Neumann, P.P.; Rudolph, M.; Schoppa, A.; et al. Hazards from failure of CNG automotive cylinders in fire. J. Hazard. Mater. 2019, 367, 1–7. [Google Scholar] [CrossRef]
  12. Chamberlain, S.; Modarres, M. Compressed natural gas bus safety: A quantitative risk assessment. Risk Anal. 2005, 25, 377–387. [Google Scholar] [CrossRef] [PubMed]
  13. Yersak, T.A.; Elhamid, M.H.A.; Dailly, A.; Rogers, M.; Jason, P.; Cai, M. Dynamics of a type IV comfortable pressure vessel for natural gas passenger vehicles. Int. J. Press. Vessel. Pip. 2019, 175, 103923. [Google Scholar] [CrossRef]
  14. Dicken, C.J.B.; Merida, W. Measured effects of filling time and initial mass on the temperature distribution within a hydrogen cylinder during refueling. J. Power Sources 2007, 165, 324–336. [Google Scholar] [CrossRef]
  15. Galassi, M.C.; Papanikolaou, E.; Heitsch, M.; Baraldi, D.; Iborra, B.A.; Moretto, P. Assessment of CFD models for hydrogen fast filling simulations. Int. J. Hydrog. Energy 2014, 39, 6252–6260. [Google Scholar] [CrossRef]
  16. Li, Q.; Zhou, J.; Chang, Q.; Xing, W. Effects of geometry and inconstant mass flow rate on temperature within a pressurized hydrogen cylinder during refueling. Int. J. Hydrog. Energy 2012, 37, 6043–6052. [Google Scholar] [CrossRef]
  17. Liu, Y.L.; Zhao, Y.Z.; Zhao, L.; Li, X.; Chen, H.G.; Zhang, L.F.; Zhao, H.; Sheng, R.H.; Xie, T.; Hu, D.H.; et al. Experimental studies on temperature rise within a hydrogen cylinder during refueling. Int. J. Hydrog. Energy 2010, 35, 2627–2632. [Google Scholar] [CrossRef]
  18. Deymi-Dashtebayaz, M.; Farzaneh-Gord, M.; Nooralipoor, N.; Rastgar, S. The full simulation of rapid refueling of a natural gas vehicle on-board cylinder. J. Nat. Gas. Sci. Eng. 2014, 21, 1099–1106. [Google Scholar] [CrossRef]
  19. Reza Kashyzadeh, K.; Kivi, S.A.; Rynkovskaya, M. Fatigue life assessment of unidirectional fibrous composite centrifugal compressor impeller blades based on FEM. Int. J. Emerg. Technol. Adv. Eng. 2017, 7, 6–11. [Google Scholar]
  20. Farrahi, G.H.; Reza Kashyzadeh, K.; Minaei, M.; Sharifpour, A.; Riazi, S. Analysis of Resistance Spot Welding Process Parameters Effect on the Weld Quality of Three-steel Sheets Used in Automotive Industry: Experimental and Finite Element Simulation. Int. J. Eng. Trans. A Basics 2020, 33, 148–157. [Google Scholar]
  21. Reza Kashyzadeh, K.; Ghorbani, S.; Forouzanmehr, M. Effects of Drying Temperature and Aggregate Shape on the Concrete Compressive Strength: Experiments and Data Mining Techniques. Int. J. Eng. 2020, 33, 1780–1791. [Google Scholar]
  22. Ghorbani, S.; Ghorbani, S.; Reza Kashyzadeh, K. Taguchi Approach and Response Surface Analysis for Design of a High-performance Single-walled Carbon Nanotube Bundle Interconnects in a Full Adder. Int. J. Eng. 2020, 33, 1598–1607. [Google Scholar]
  23. Omidi Bidgoli, M.; Reza Kashyzadeh, K.; Rahimian Koloor, S.S.; Petru, M. Estimation of Critical Dimensions for the Crack and Pitting Corrosion Defects in the Oil Storage Tank Using Finite Element Method and Taguchi Approach. Metals 2020, 10, 1372. [Google Scholar] [CrossRef]
  24. Reza Kashyzadeh, K.; Farrahi, G.H.; Shariyat, M.; Ahmadian, M.T. Experimental and finite element studies on free vibration of automotive steering knuckle. Int. J. Eng. 2017, 30, 1776–1783. [Google Scholar]
  25. Vlase, S.; Purcarea, R.; Teodorescu-Draghicescu, H.; Calin, R.; Szava, I.; Mihalcica, M. Behavior of a new Heliopol/Stratimat300 composite laminate. J. Optoelectron. Adv. Mater. 2013, 78, 569–572. [Google Scholar]
  26. Reza Kashyzadeh, K. A new algorithm for fatigue life assessment of automotive safety components based on the probabilistic approach: The case of the steering knuckle. Int. J. Eng. Sci. Technol. 2020, 23, 392–404. [Google Scholar] [CrossRef]
  27. Reza Kashyzadeh, K.; Farrahi, G.H.; Shariyat, M.; Ahmadian, M.T. Experimental accuracy assessment of various high-cycle fatigue criteria for a critical component with a complicated geometry and multi-input random non-proportional 3D stress components. Eng. Fail. Anal. 2018, 90, 534–553. [Google Scholar] [CrossRef]
  28. Reza Kashyzadeh, K. Effects of Axial and Multiaxial Variable Amplitude Loading Conditions on the Fatigue Life Assessment of Automotive Steering Knuckle. J. Fail. Anal. Preven. 2020, 20, 455–463. [Google Scholar] [CrossRef]
  29. Maleki, E.; Unal, O.; Reza Kashyzadeh, K. Efficiency analysis of shot peening parameters on variations of hardness, grain size and residual stress via Taguchi approach. Met. Mater. Int. 2019, 25, 1436–1447. [Google Scholar] [CrossRef]
  30. Fawaz, Z.; Ellyin, F. Fatigue failure model for fiber-reinforced materials under general loading condition. J. Compos. Mater. 1994, 28, 1432–1451. [Google Scholar] [CrossRef]
Figure 1. The proposed algorithm for fatigue life improvement of a composite type-4 compressed natural gas (CNG) tank.
Figure 1. The proposed algorithm for fatigue life improvement of a composite type-4 compressed natural gas (CNG) tank.
Polymers 13 00483 g001
Figure 2. Schematic of the CNG tank with all dimensions in mm [5].
Figure 2. Schematic of the CNG tank with all dimensions in mm [5].
Polymers 13 00483 g002
Figure 3. The naming of different parts of the composite shell based on the wall-thickness in a sheared tank.
Figure 3. The naming of different parts of the composite shell based on the wall-thickness in a sheared tank.
Polymers 13 00483 g003
Figure 4. Finite element model of the CNG fuel tank with 91,622 elements [5].
Figure 4. Finite element model of the CNG fuel tank with 91,622 elements [5].
Polymers 13 00483 g004
Figure 5. Details of the cyclic loading used to predict the fatigue life of the CNG fuel tank.
Figure 5. Details of the cyclic loading used to predict the fatigue life of the CNG fuel tank.
Polymers 13 00483 g005
Figure 6. The schematic of the response surface method (RSM)-based design of experiments used in this study.
Figure 6. The schematic of the response surface method (RSM)-based design of experiments used in this study.
Polymers 13 00483 g006
Figure 7. The contour of principal stress in the initial design of the CNG tank.
Figure 7. The contour of principal stress in the initial design of the CNG tank.
Polymers 13 00483 g007
Figure 8. Maximum principal stress created in the critical area of the composite CNG tank in terms of various parameters including (A) variation of the fiber angle and the assumption T1 = 11 mm, T2 = 4.4 mm, and T3 = 3.3 mm, (B) variation of the wall-thickness of the top and bottom sections and the assumption that the fiber angle = 17, T2 = 4.4 mm, and T3 = 3.3 mm, (C) variation of the wall-thickness of the middle section and the assumption that the fiber angle = 17, T1 = 11 mm, and T3 = 3.3 mm, and (D) variation of the wall-thickness of the main section and the assumption that the fiber angle = 17, T1 = 11 mm, and T2 = 4.4 mm.
Figure 8. Maximum principal stress created in the critical area of the composite CNG tank in terms of various parameters including (A) variation of the fiber angle and the assumption T1 = 11 mm, T2 = 4.4 mm, and T3 = 3.3 mm, (B) variation of the wall-thickness of the top and bottom sections and the assumption that the fiber angle = 17, T2 = 4.4 mm, and T3 = 3.3 mm, (C) variation of the wall-thickness of the middle section and the assumption that the fiber angle = 17, T1 = 11 mm, and T3 = 3.3 mm, and (D) variation of the wall-thickness of the main section and the assumption that the fiber angle = 17, T1 = 11 mm, and T2 = 4.4 mm.
Polymers 13 00483 g008
Figure 9. RSM results for maximum principal stress in a variable wall-thickness composite CNG tank including (A) variations of the fiber angle and the wall-thickness of the top and bottom sections, (B) variations of the fiber angle and the wall-thickness of the middle section, (C) variations of the fiber angle and the wall-thickness of the main section, (D) variations of the wall-thickness of the top and bottom sections and the wall-thickness of the middle section, (E) variations of the wall-thickness of the top and bottom sections and the wall-thickness of the main section, and (F) variations of the wall-thickness of the middle section and the wall-thickness of the main section (all symbols are defined in Table 3).
Figure 9. RSM results for maximum principal stress in a variable wall-thickness composite CNG tank including (A) variations of the fiber angle and the wall-thickness of the top and bottom sections, (B) variations of the fiber angle and the wall-thickness of the middle section, (C) variations of the fiber angle and the wall-thickness of the main section, (D) variations of the wall-thickness of the top and bottom sections and the wall-thickness of the middle section, (E) variations of the wall-thickness of the top and bottom sections and the wall-thickness of the main section, and (F) variations of the wall-thickness of the middle section and the wall-thickness of the main section (all symbols are defined in Table 3).
Polymers 13 00483 g009aPolymers 13 00483 g009b
Figure 10. Pareto chart of the standardized effects.
Figure 10. Pareto chart of the standardized effects.
Polymers 13 00483 g010
Table 1. Mechanical properties of orthotropic materials [5].
Table 1. Mechanical properties of orthotropic materials [5].
ParameterSymbolUnitPolyethyleneValue
Glass/Epoxy Composite
Longitudinal tension modulusE11GPa0.935
Transverse tension modulusE22GPa0.93.5
Shear modulusG12GPa0.3124.6
Poisson’s ratioϑ12 0.450.26
Table 2. Mechanical properties of 7xxx series aluminum alloy (Al 7075-T6) [5].
Table 2. Mechanical properties of 7xxx series aluminum alloy (Al 7075-T6) [5].
ParameterSymbolUnitValue
Tensile ultimate strengthσutMPa572
Tensile yield strengthσytMPa503
Elastic modulusEGPa71.7
Shear modulusGGPa26.9
Poisson’s ratioϑ 0.33
Fatigue limitσfMPa159
Table 3. Variables and their levels are used as input data in the response surface method (RSM).
Table 3. Variables and their levels are used as input data in the response surface method (RSM).
Levels
ParametersSymbolL1L2L3L4L5
Fiber angleA1017202530
Wall-thickness of the top and bottom sectionsT18.89.91112.113.2
Wall-thickness of the middle sectionT23.34.45.56.67.7
Wall-thickness of the main sectionT32.23.34.45.56.6
Table 4. Finite element results for the different settings of wall-thickness parameters and fiber angle based on the RSM.
Table 4. Finite element results for the different settings of wall-thickness parameters and fiber angle based on the RSM.
Simulation No.Level of Input ParametersOutput (MPa)
AT1T2T3
1L5L3L3L3194.84
2L3L3L3L3202.77
3L2L4L4L4203.68
4L4L4L2L2201.27
5L3L1L3L3200.35
6L4L4L4L2196.51
7L4L2L4L4196.78
8L2L4L4L2202.06
9L3L3L3L1199.75
10L3L5L3L3204.23
11L3L3L3L5203.64
12L4L2L2L2199.56
13L4L2L2L4200.97
14L2L2L4L2200.21
15L2L4L2L4208.49
16L3L3L1L3207.85
17L2L2L4L4201.76
18L4L4L4L4198.35
19L4L2L4L2195.03
20L2L4L2L2207.23
21L3L3L5L3198.46
22L2L2L2L2205.17
23L2L2L2L4206.39
24L1L3L3L3207.79
Table 5. A comparison among the parameters of initial and optimum designs of the composite CNG tank.
Table 5. A comparison among the parameters of initial and optimum designs of the composite CNG tank.
DesignFiber AngleWall-Thickness of the Top and Bottom Sections (mm)Wall-Thickness of the Middle Section (mm)Wall-Thickness of the Main Section (mm)Maximum Principal Stress (MPa)
Initial17114.43.3206.32
Optimum308.87.72.2185.45
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kashyzadeh, K.R.; Rahimian Koloor, S.S.; Omidi Bidgoli, M.; Petrů, M.; Amiri Asfarjani, A. An Optimum Fatigue Design of Polymer Composite Compressed Natural Gas Tank Using Hybrid Finite Element-Response Surface Methods. Polymers 2021, 13, 483. https://0-doi-org.brum.beds.ac.uk/10.3390/polym13040483

AMA Style

Kashyzadeh KR, Rahimian Koloor SS, Omidi Bidgoli M, Petrů M, Amiri Asfarjani A. An Optimum Fatigue Design of Polymer Composite Compressed Natural Gas Tank Using Hybrid Finite Element-Response Surface Methods. Polymers. 2021; 13(4):483. https://0-doi-org.brum.beds.ac.uk/10.3390/polym13040483

Chicago/Turabian Style

Kashyzadeh, Kazem Reza, Seyed Saeid Rahimian Koloor, Mostafa Omidi Bidgoli, Michal Petrů, and Alireza Amiri Asfarjani. 2021. "An Optimum Fatigue Design of Polymer Composite Compressed Natural Gas Tank Using Hybrid Finite Element-Response Surface Methods" Polymers 13, no. 4: 483. https://0-doi-org.brum.beds.ac.uk/10.3390/polym13040483

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