METHOD AND SYSTEM FOR OPTIMIZING DESIGN OF HIGH-PRESSURE HYDROGEN STORAGE COMPOSITE WOUND GAS CYLINDER

Abstract

Provided are a method and a system for optimizing a design of a high-pressure hydrogen storage composite wound gas cylinder. The method includes: S1, determining a preliminary design result according to a film theory and a grid theory; S2, establishing a three-dimensional solid refined model of the preliminary design result based on a finite element analysis software; and S3, applying the three-dimensional solid refined model for a simulation analysis, and applying a response surface or a genetic algorithm for an optimization design.

Claims

1. A method for optimizing a design of a high-pressure hydrogen storage composite wound gas cylinder, comprising following steps: S1, determining a preliminary design result according to a film theory and a grid theory, wherein the preliminary design result comprises a thickness of an inner container of a high-pressure composite wound gas cylinder and a thickness of a composite reinforced layer; S2, establishing a three-dimensional solid refined model of the preliminary design result based on a finite element analysis method; and S3, carrying out a multi-objective optimization design of the high-pressure composite wound gas cylinder by applying a response surface or a genetic algorithm based on the three-dimensional solid refined model.

2. The method for optimizing the design of the high-pressure hydrogen storage composite wound gas cylinder according to claim 1, wherein a method for calculating the thickness of the inner container of the high-pressure composite material wound gas cylinder comprises following steps: $t_i=\frac{P_{i}D}{2R_{\mathrm{mi}}+P_i}$ wherein t.sub.i represents the thickness of the inner container, P.sub.i represents a minimum design bursting pressure of the inner container, D represents an outer diameter of the inner container, and R.sub.mi represents a minimum guaranteed value of a tensile strength of an inner container material after heat treatment; a minimum thickness of a neck diameter of the inner container is calculated according to a following formula: $t_{i0\min }=\frac{P_b{d}_{i0}}{2R_{\mathrm{mi}}+P_b}$ wherein t.sub.i0min represents the minimum thickness of the neck diameter of the inner container, P.sub.b represents an actual design bursting pressure of the gas cylinder, and d.sub.i0 represents an outer diameter of the inner container bottle mouth.

3. The method for optimizing the design of the high-pressure hydrogen storage composite wound gas cylinder according to claim 2, wherein a method for calculating the thickness of the composite reinforced layer comprises: $t_{\mathrm{os}}=\frac{P_{b}D-4R_{\mathrm{mi}}{t}_i}{4k{\cos }^2{R}_{\mathrm{mo}}}$ wherein t.sub.os represents a layer thickness of a spiral layer of the composite reinforced layer, R.sub.mo represents a minimum guaranteed value of a tensile strength of the composite, represents a winding angle and k represents a balance coefficient; $t_{\mathrm{oh}}=\frac{P_{b}D}{4R_{\mathrm{mo}}}\left(2-\frac{1}{k}{\tan }^2\right)-\frac{R_{\mathrm{mi}}{t}_i}{R_{\mathrm{mo}}}$ wherein t.sub.oh represents a thickness of a circumferential layer of the composite reinforced layer, and a circumferential reinforcement structure consists of two parts, one part is provided by the circumferential layer, and an other part is provided by the spiral layer; and a calculation formula of the thickness t.sub.o of the composite reinforced layer is: $t_o=t_{\mathrm{os}}+t_{\mathrm{oh}}.$

4. The method for optimizing the design of the high-pressure hydrogen storage composite wound gas cylinder according to claim 1, wherein in the S2, a process of establishing the three-dimensional solid refined model comprises an establishment of the inner container model and an establishment of the composite reinforced layer.

5. The method for optimizing the design of the high-pressure hydrogen storage composite wound gas cylinder according to claim 4, wherein a method for establishing the composite reinforced layer comprises following steps: the reinforced layer of composite material may be divided into a circumferential layer modeling and a spiral layer modeling; the circumferential layer is performed with a ply stacking modeling by using 90 fibers, the spiral layer is performed with a ply stacking modeling according to a preliminarily designed ply angle, and the ply stacking modeling method with variable angle and thickness is used at a head; the modified cubic spline thickness formula is used for a data calculation and then the ply stacking modeling is carried out.

6. The method for optimizing the design of the high-pressure hydrogen storage composite wound gas cylinder according to claim 4, wherein a method for carrying out a ply stacking modeling after the data calculation by using the modified cubic spline thickness formula comprises: a modified cubic spline function for calculating a thickness of the head is as follows: $t\left(r_i\right)=m_1{r}_i^0+m_2{r}_i^1+m_3{r}_i^2+m_4{r}_i^3$ wherein m.sub.1, m.sub.2, m.sub.3, m.sub.4 represent undetermined coefficients and r.sub.i represents a distance from a central axis; a formula for determining m.sub.1, m.sub.2, m.sub.3, m.sub.4 is as follows: $\left[\begin{array}{c}{m}_1\\ m_2\\ m_3\\ m_4\end{array}\right]={\left[\begin{array}{cccc}1& r_0& r_0^2& r_0^3\\ 1& r_{2b}& r_{2b}^2& r_{2b}^3\\ 0& 1& 2r_{2b}& 3r_{2b}^2\\ 1& r_{2\text{.2}b}& r_{2.2b}^2& r_{2.2b}^3\end{array}\right]}^{-1}\left[\begin{array}{c}{t}_R.Math.R.Math.\cos {}_0/\left(m_0.Math.b\right)\\ \frac{m_R.Math.n_R}{}.Math.\left[\mathrm{arc}\cos \left(\frac{r_0}{r_{2b}}\right)-\mathrm{arc}\cos \left(\frac{r_0+b}{r_{2b}}\right)\right].Math.t_p\\ \frac{m_R.Math.n_R}{}.Math.\left(\frac{r_0}{r_{2b}\sqrt{r_{2b}^2-r_0^2}}-\frac{r_b}{r_{2b}\sqrt{r_{2b}^2-r_b^2}}\right).Math.t_p\\ \frac{m_R.Math.n_R}{}.Math.\left[\mathrm{arc}\cos \left(\frac{r_0}{r_{2.2b}}\right)-\mathrm{arc}\cos \left(\frac{r_0+b}{r_{2.2b}}\right)\right].Math.t_p\end{array}\right]$ wherein m.sub.R represents a number of yarn bands in each single layer of a cylinder body, n.sub.R represents a total number of single layers, t.sub.p represents a thickness of single layer yarn bands, m.sub.0 represents a number of yarn pieces at a polar hole, m.sub.R.Math.n.sub.R represents a total number of yarn bands, and m.sub.R.Math.n.sub.R.Math.t.sub.p represents a total thickness of yarn bands in the cylinder body; formulas for calculating a winding angle at the polar hole of the head are as follows: a spherical head formula: $=\mathrm{arc}\sin \left(\frac{r_0}{R}\right)\frac{180}{}$ an elliptical head formula: $=\mathrm{arc}\sin \left(\frac{r_0^2}{r_0^2+\left(R^2-r_0^2\right)\frac{h^2}{R^2}}\right)\frac{180}{}$ wherein r.sub.0 represents the polar hole radius, and R represents an outside radius of the cylinder body; formulas for calculating a winding angle of the head are as follows: a spherical head formula: $=\mathrm{arc}\sin \left(\frac{r_0}{r_i}\right)\frac{180}{}$ an elliptical head formula: $=\mathrm{arc}\cos \left(\cos .Math.\cos \sqrt{\frac{R^2-\left(R^2-h^2\right).Math.\frac{{\sin }^2}{{\sin }^2+{\cos }^2.Math.{\sin }^2}}{R^2-\left(R^2-h^2\right).Math.\frac{{\cos }^2.Math.{\cos }^2.Math.{\sin }^2}{{\sin }^2+{\cos }^2.Math.{\sin }^2}}}\right)\frac{180}{}$ wherein r.sub.i represents a radius from point i at the head to a rotating shaft, represents an included angle between a plane where a geodesic line is located and an XZ plane, and represents an included angle between a plane of meridian and the XZ plane, and a method for calculating and is as follows: ${\sin }^2=\frac{r_0^2}{R^2}$ ${\sin }^2=\frac{R^2-r^2}{r^2}.Math.\frac{r_0^2}{R^2-r_0^2}.$

7. The method for optimizing the design of the high-pressure hydrogen storage composite wound gas cylinder according to claim 1, wherein an application and a solution of boundary conditions are further comprised in the S2: one end of the gas cylinder is fixed or other fixed constraints are selected; a load is pressurized on an inner surface of the inner container, wherein a self-tightening pressure of a metal inner container needs to be calculated, and CAE analysis and calculation are carried out after loading; a pressure loading process is as follows: process 1: an internal pressure of the gas cylinder rises from a zero pressure to the self-tightening pressure; process 2: the internal pressure of the gas cylinder declines from the self-tightening pressure to the zero pressure; process 3: the internal pressure of the gas cylinder rises from the zero pressure to a nominal working pressure; process 4: the internal pressure of the gas cylinder rises from the nominal working pressure to a hydraulic fatigue test pressure; process 5: the internal pressure of the gas cylinder rises from the hydraulic fatigue test pressure to a hydraulic test pressure of the gas cylinder; and process 6: the internal pressure of the gas cylinder rises from the hydraulic test pressure of the gas cylinder to a hydraulic bursting pressure of the gas cylinder; wherein the self-tightening pressure is estimated by several groups of data for trial calculation, and then the self-tightening pressure is determined according to following restrictions: a, under a working pressure, a fiber stress of the winding layer cannot exceed 30% of a strength limit; b, under the working pressure, a stress of the metal inner container cannot exceed 60% of a yield limit; c, under the zero pressure, a compressive stress of the metal inner container cannot exceed 95% of the yield limit, but is above 60% of the yield limit; and d, the maximum stress at the head of the winding layer of the gas cylinder is always less than the maximum stress at the cylinder body.

8. A system for optimizing a design of a high-pressure hydrogen storage composite wound gas cylinder, comprising a preliminary design result acquisition module, a model establishment module and a multi-objective optimization design module; the preliminary design result acquisition module is used for determining a preliminary design result according to a film theory and a grid theory, and the preliminary design result comprises a thickness of an inner container of a high-pressure composite wound gas cylinder and a thickness of a composite reinforced layer; the model establishment module is used for establishing a three-dimensional solid refined model of the preliminary design result based on a finite element analysis method; and the multi-objective optimization design module is used for carrying out a multi-objective optimization design of the high-pressure composite wound gas cylinder by applying a response surface or a genetic algorithm based on the three-dimensional solid refined model.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0051] In order to explain the technical scheme of this disclosure more clearly, the drawings needed in the embodiments are briefly introduced below. Obviously, the drawings in the following description are only some embodiments of this disclosure. For ordinary technicians in this field, other drawings may be obtained according to these drawings without paying creative labor.

[0052] FIG. 1 is a method step diagram of a method and a system for optimizing a design of a high-pressure hydrogen storage composite wound gas cylinder according to an embodiment of the present disclosure.

[0053] FIG. 2 is a method step diagram of carrying out a multi-objective optimization design of the high-pressure composite wound gas cylinder based on a response surface or a genetic algorithm in the method and the system for optimizing the design of the high-pressure hydrogen storage composite wound gas cylinder according to an embodiment of the present disclosure.

[0054] FIG. 3 is a schematic diagram of a geodesic line and a meridian line according to an embodiment of the present disclosure.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0055] In the following, technical schemes in the embodiments of the disclosure is clearly and completely described with reference to the attached drawings in the embodiments of the disclosure. Obviously, the described embodiments are only a part of embodiments of the disclosure, but not all embodiments. Based on the embodiments in this disclosure, all other embodiments obtained by ordinary technicians in this field without creative work belong to the protection scope of this disclosure.

[0056] In order to make the above objects, features and advantages of this disclosure more obvious and easier to understand, the disclosure is further described in detail with the attached drawings and specific embodiments.

Embodiment 1

[0057] In this embodiment, as shown in FIG. 1, a method for optimizing a design of a high-pressure hydrogen storage composite wound gas cylinder includes: [0058] S1, determining a preliminary design result according to a film theory and a grid theory, where the preliminary design result includes a thickness of an inner container of a high-pressure composite wound gas cylinder and a thickness of a composite reinforced layer; [0059] in the S1, a process of determining the preliminary design result according to the thin film theory and the grid theory includes: [0060] calculating the thickness of the inner container and the thickness of the reinforced layer of the high-pressure composite wound gas cylinder.

[0061] A method for calculating the thickness of the inner container of the high-pressure composite material wound gas cylinder includes following steps: [0062] (1) a formula of the thickness of the inner container is as follows:

[00013]$t_i=\frac{P_{i}D}{2R_{\mathrm{mi}}+P_i}$ [0063] where t.sub.i represents the thickness of the inner container, P.sub.i represents a minimum design bursting pressure of the inner container, D represents an outer diameter of the inner container, and R.sub.mi represents a minimum guaranteed value of a tensile strength of an inner container material after heat treatment; [0064] a minimum thickness of a neck diameter of the inner container is calculated according to a following formula:

[00014]$t_{i0\min }=\frac{P_b{d}_{i0}}{2R_{\mathrm{mi}}+P_b}$ [0065] where t.sub.i0min represents the minimum thickness of the neck diameter of the inner container, P.sub.b represents an actual design bursting pressure of the gas cylinder, and d.sub.i0 represents an outer diameter of the inner container bottle mouth. [0066] (2) a formula of a composite reinforced layer is as follows:

[00015]$t_{\mathrm{os}}=\frac{P_{b}D-4R_{\mathrm{mi}}{t}_i}{4k{\cos }^2{R}_{\mathrm{mo}}}$ [0067] where t.sub.os represents a layer thickness of a spiral layer of the composite reinforced layer, R.sub.mo represents a minimum guaranteed value of a tensile strength of the composite, represents a winding angle, k represents a balance coefficient, and the value range is 0.7-0.8 according to experience;

[00016]$t_{\mathrm{oh}}=\frac{P_{b}D}{4R_{\mathrm{mo}}}\left(2-\frac{1}{k}{\tan }^2\right)-\frac{R_{\mathrm{mi}}{t}_i}{R_{\mathrm{mo}}}$ [0068] where t.sub.oh represents a thickness of a circumferential layer of the composite reinforced layer, and a circumferential reinforcement structure consists of two parts, one part is provided by the circumferential layer, and the other part is provided by the spiral layer.

[0069] A calculation formula of a minimum winding angle .sub.0 is as follows:

[00017]${}_0={\sin }^{-1}\frac{d_{i0}}{D}$

[0070] When different winding angles are adopted to avoid pole hole accumulation and meet process requirements, a cosine or sine square ratio may be used for conversion. [0071] a calculation formula of the thickness t.sub.o of the composite reinforced layer is as follows:

[00018]$t_o=t_{\mathrm{os}}+t_{\mathrm{oh}}.$

[0072] Finally, through the calculation results and experience, as well as the thickness of single-layer composite, the ply stacking sequence, ply stacking angle and ply stacking thickness are preliminarily determined, and then the refined model is analyzed and designed according to a preliminary design result.

[0073] S2, establishing a three-dimensional solid refined model of the preliminary design result based on a finite element analysis method:

(1) Establishing a Three-Dimensional Solid Refined Model

[0074] The refined model establishing is divided into two parts, one part is the establishment of an inner container model, and the other part is the establishment of a composite reinforced layer. The inner container is made of metal material or non-metal material, the metal inner container is a revolving body, and the non-metal inner container is a revolving assembly (the cylinder is non-metal, and the polar hole is metal); the establishment of the composite reinforced layer is complicated, and is divided into a circumferential layer modeling and a spiral layer modeling. The circumferential layer is performed with a ply stacking modeling by using 90 fibers, the spiral layer is performed with a ply stacking modeling according to a preliminarily designed ply angle, and the ply stacking modeling method with variable angle and thickness is used at the head. The modified cubic spline thickness formula is used for a data calculation and then the ply stacking modeling is carried out.

[0075] A process of establishing the three-dimensional solid refined model includes: the establishment of the inner container model, and the establishment of the composite reinforced layer.

[0076] A method for establishing the composite reinforced layer includes following steps: [0077] the reinforced layer of composite material may be divided into a circumferential layer modeling and a spiral layer modeling. The circumferential layer is performed with a ply stacking modeling by using 90 fibers, the spiral layer is performed with a ply stacking modeling according to a preliminarily designed ply angle, and the ply stacking modeling method with variable angle and thickness is used at the head. The modified cubic spline thickness formula is used for a data calculation and then the ply stacking modeling is carried out. [0078] a method for carrying out a ply stacking modeling after the data calculation by using the modified cubic spline thickness formula includes: [0079] a modified cubic spline function for calculating a thickness of the head is as follows:

[00019]$t\left(r_i\right)=m_1{r}_i^0+m_2{r}_i^1+m_3{r}_i^2+m_4{r}_i^3$ [0080] where m.sub.1, m.sub.2, m.sub.3, m.sub.4 represent undetermined coefficients and r.sub.i represents a distance from a central axis;

[0081] The modified cubic spline function uses the calculation formula of r.sub.2.2b instead of the conventional calculation formula of fiber volume in two bandwidth ranges, the calculation formula of r.sub.2.2b is more convenient for calculation and also avoids a peak of fiber accumulation within one-time bandwidth, a formula for determining m.sub.1, m.sub.2, m.sub.3, m.sub.4 is as follows:

[00020]$\left[\begin{array}{c}{m}_1\\ m_2\\ m_3\\ m_4\end{array}\right]={\left[\begin{array}{cccc}1& r_0& r_0^2& r_0^3\\ 1& r_{2b}& r_{2b}^2& r_{2b}^3\\ 0& 1& 2r_{2b}& 3r_{2b}^2\\ 1& r_{2.2b}& r_{2.2b}^2& r_{2.2b}^3\end{array}\right]}^{-1}\left[\begin{array}{c}{t}_R.Math.R.Math.\cos {}_0/\left(m_0.Math.b\right)\\ \frac{m_R.Math.n_R}{}.Math.\left[\mathrm{arc}\cos \left(\frac{r_0}{r_{2b}}\right)-\mathrm{arc}\cos \left(\frac{r_0+b}{r_{2b}}\right)\right].Math.t_p\\ \frac{m_R.Math.n_R}{}.Math.\left(\frac{r_0}{r_{2b}\sqrt{r_{2b}^2-r_0^2}}-\frac{r_b}{r_{2b}\sqrt{r_{2b}^2-r_0^2}}\right).Math.t_p\\ \frac{m_R.Math.n_R}{}.Math.\left[\mathrm{arc}\cos \left(\frac{r_0}{r_{2.2b}}\right)-\mathrm{arc}\cos \left(\frac{r_0+b}{r_{2.2b}}\right)\right].Math.t_p\end{array}\right]$ [0082] where m.sub.R represents a number of yarn bands in each single layer of a cylinder body, n.sub.R represents a total number of single layers, t.sub.p represents a thickness of single layer yarn bands, m.sub.0 represents a number of yarn pieces at a polar hole, m.sub.R.Math.n.sub.R represents a total number of yarn bands, and m.sub.R.Math.n.sub.R.Math.t.sub.p represents a total thickness of yarn bands in the cylinder body.

[0083] Formulas for calculating a winding angle at the polar hole of the head are as follows: [0084] a spherical head formula:

[00021]$=\mathrm{arc}\sin \left(\frac{r_0}{R}\right)\frac{180}{}$ [0085] an elliptical head formula:

[00022]$=\mathrm{arc}\sin \left(\frac{r_0^2}{r_0^2+\left(R^2-r_0^2\right)\frac{h^2}{R^2}}\right)\frac{180}{}$ [0086] where r.sub.0 represents the polar hole radius, and R represents an outside radius of the cylinder body; [0087] formulas for calculating a winding angle of the head are as follows: [0088] a spherical head formula:

[00023]$=\mathrm{arc}\sin \left(\frac{r_0}{r_i}\right)\frac{180}{}$ [0089] an elliptical head formula:

[00024]$=\mathrm{arc}\cos \left(\cos .Math.\cos \sqrt{\frac{R^2-\left(R^2-h^2\right).Math.\frac{{\sin }^2}{{\sin }^2+{\cos }^2.Math.{\sin }^2}}{R^2-\left(R^2-h^2\right).Math.\frac{{\cos }^2.Math.{\cos }^2.Math.{\sin }^2}{{\sin }^2+{\cos }^2.Math.{\sin }^2}}}\right)\frac{180}{}$ [0090] where r.sub.i represents a radius from point i at the head to a rotating shaft, represents an included angle between a plane where a geodesic line is located and an XZ plane, and represents an included angle between a plane of meridian and the XZ plane, as shown in FIG. 3, a method for calculating and is as follows:

[00025]${\sin }^2=\frac{r_0^2}{R^2}$${\sin }^2=\frac{R^2-r^2}{r^2}.Math.\frac{r_0^2}{R^2-r_0^2}.$

(2) Material Attribute Setting

[0091] Failure criteria need to be added to the material property setting, such as three-dimensional Hashin and LaRC04 failure criteria, or progressive failure model to improve the calculation accuracy, but the calculation time becomes longer.

(3) Application and Solution of Boundary Conditions

[0092] One end of the gas cylinder is fixed or other fixed constraints are selected; the load is pressurized on an inner surface of the inner container, where a self-tightening pressure of a metal inner container needs to be calculated, and a CAE analysis and calculation are carried out after loading. The pressure loading process is as follows: [0093] process 1: the internal pressure of the gas cylinder rises from a zero pressure to the self-tightening pressure; [0094] process 2: the internal pressure of the gas cylinder declines from the self-tightening pressure to the zero pressure; [0095] process 3: the internal pressure of the gas cylinder rises from the zero pressure to a nominal working pressure; [0096] process 4: the internal pressure of the gas cylinder rises from the nominal working pressure to a hydraulic fatigue test pressure (1.3 times the working pressure); [0097] process 5: the internal pressure of the gas cylinder rises from the hydraulic fatigue test pressure to a hydraulic test pressure of the gas cylinder (1.5 times the working pressure); and [0098] process 6: the internal pressure of the gas cylinder rises from the hydraulic test pressure of the gas cylinder to a hydraulic bursting pressure of the gas cylinder (2.35 times the working pressure); [0099] where the self-tightening pressure is estimated by several groups of data for trial calculation, and then the self-tightening pressure is determined according to following restrictions: [0100] a, under the working pressure, the fiber stress of the winding layer cannot exceed 30% of a strength limit; [0101] b, under the working pressure, the stress of the metal inner container cannot exceed 60% of a yield limit; [0102] c, under the zero pressure, a compressive stress of the metal inner container cannot exceed 95% of the yield limit, but is above 60% of the yield limit; and [0103] d, the maximum stress at the head of the winding layer of the gas cylinder is always less than the maximum stress at the cylinder body.

[0104] In the process of CAE solution, nonlinear and large deformation options need to be turned on to obtain more accurate results.

(4) Post-Processing of Results

[0105] The overall stress evaluation of gas cylinder inner container may be based on Von Mises stress less than tensile strength, and the local stress may be evaluated according to the stress analysis method provided by analysis and design standards, that is, the primary bending stress is less than 1.5 times the yield strength, the local film stress is less than 1.5 times the yield strength, and the primary stress plus the secondary stress is less than 3 times the yield strength.

[0106] The stress evaluation of gas cylinder winding layer is evaluated according to LaRC04 criterion, which may be used to judge whether the gas cylinder winding layer is invalid under the blasting pressure, and this method may also be used to calculate the actual blasting pressure.

[0107] S3, carrying out a multi-objective optimization design of the high-pressure composite wound gas cylinder by applying a response surface or a genetic algorithm based on the three-dimensional solid refined model.

[0108] Based on the response surface or the genetic algorithm, the multi-objective optimization design of high-pressure composite wound gas cylinder is carried out. The multi-objective optimization design of high-pressure composite wound gas cylinder takes ply stacking angle, ply stacking sequence and ply stacking number as design variables, takes weight of gas cylinder winding layer, fiber direction stress and fiber tangential stress as objective functions and takes the blasting pressure as a constraint condition.

[0109] Optionally, in S3, a process of applying the response surface or the genetic algorithm to carry out the optimization design includes: [0110] carrying out a simulation analysis according to the built refined simulation model; and [0111] applying the response surface or the genetic algorithm to carry out the optimization design.

[0112] A mathematical model is built as follows:

[00026]$\min F\left(x\right)=\left(f_1\left(x\right),f_2\left(x\right),f_3\left(x\right)\right)$ [0113] where F(x) represents the objective function, .sub.1(x) represents the mass of winding layer; .sub.2(x) represents the fiber direction stress, .sub.3(x) represents the fiber tangential stress, and P(x)P.sub.b represents the constraint condition.

[0114] The mathematical model is an optimized design objective function model.

[0115] As shown in FIG. 2, firstly, the simulation model is established, and then the refined model is used for simulation analysis. Finally, the response surface or the genetic algorithm is used for optimization design. According to the computer capacity, the optimization method, the sample number and the iteration times are selected. Generally, the number of candidate points is three, and the minimum value of the objective function meeting the constraint condition is the optimization result. Otherwise, the parameters are modified and the optimization program is continued.

Embodiment 2

[0116] The disclosure relates to a system for optimizing a design of a high-pressure hydrogen storage composite wound gas cylinder, the system includes a preliminary design result acquisition module, a model establishment module and a multi-objective optimization design module; [0117] the preliminary design result acquisition module is used for determining a preliminary design result according to a film theory and a grid theory, and the preliminary design result includes a thickness of an inner container of a high-pressure composite wound gas cylinder and a thickness of a composite reinforced layer; [0118] the model establishment module is used for establishing a three-dimensional solid refined model of the preliminary design result based on a finite element analysis method; and [0119] the multi-objective optimization design module is used for carrying out a multi-objective optimization design of the high-pressure composite wound gas cylinder by applying a response surface or a genetic algorithm based on the three-dimensional solid refined model.

[0120] The above-mentioned embodiments are only a description of the preferred mode of this disclosure, not a limitation on the scope of this disclosure. Without departing from the design spirit of this disclosure, various modifications and improvements made by ordinary technicians in this field to the technical scheme of this disclosure shall fall within the protection scope determined by the claims of the disclosure.