Research on the integral forming process of thin walled and thick mouth seamless gas cylinders
- Select a language for the TTS:
- UK English Female
- UK English Male
- US English Female
- US English Male
- Australian Female
- Australian Male
- Language selected: (auto detect) - EN

Play all audios:

ABSTRACT There is a considerable difference in wall thickness between the mouth and the cavity of thin-walled and thick-mouthed seamless gas cylinders, and the existing manufacturing
processes are unable to effectively meet product requirements. To overcome such issue, a step-by-step boring-necking-spinning solution for gas cylinders was proposed, in which sufficient
wall thickness is reserved for the mouth area of the cylinder blank, followed by necking-spinning to realize the overall forming of thin-walled, thick-mouthed seamless gas cylinders. The
stress–strain distribution and geometric dimensional changes of gas cylinders during the spinning process were investigated by means of finite element simulation, and the effects of
different process parameters on the stress and wall thickness of the bottle mouth were analyzed. Further, multi-objective optimization of the response surface model was performed using the
NSGA-II algorithm to derive a set of optimal process parameters. Finally, the correctness of the simulation and optimization results was verified experimentally, and the expected geometry
and optimal strain state of the gas cylinder were obtained. The newly developed processing solution represents a groundbreaking advancement in the manufacturing of thin-walled and
thick-mouthed gas cylinders. SIMILAR CONTENT BEING VIEWED BY OTHERS OPTIMIZED DESIGN AND KEY PERFORMANCE FACTORS OF A GAS CIRCULATION FILTRATION SYSTEM IN A METAL 3D PRINTER Article Open
access 22 August 2022 IMPACT BEHAVIOUR OF 3D PRINTED CELLULAR STRUCTURES FOR MOUTHGUARD APPLICATIONS Article Open access 07 March 2022 DYNAMIC SEALING SIMULATION AND PERFORMANCE OPTIMIZATION
OF CONICAL RUBBER CORE IN ROTARY BLOWOUT PREVENTER Article Open access 19 March 2025 INTRODUCTION As human civilization advances and science and technology continue to develop, the impact
of energy consumption and global warming has attracted significant attention worldwide. Consequently, research into alternative energy sources, including the transportation and storage of
clean energy sources such as hydrogen and natural gas, has become a priority. Gas cylinders are widely used in the storage and transportation of clean energy. Azeem et al.1 highlighted that
a seamless gas cylinder liner, composed of a metal liner with a high-strength carbon fiber reinforced composite wrapped around the outer layer, is a promising option for high-pressure
storage containers. Such design is advantageous owing to the lighter weight, increased strength, gas-tight properties, and adequate safety features. Due to the limited volume of gas
cylinders, gases are usually stored under high pressure. Such conditions require cylinders to be light in weight and large in volume, while ensuring sufficient strength. Therefore, the
structural design of gas cylinders necessitates a thin wall thickness for the cavity to minimize weight. However, the mouth of the cylinder must have a sufficient thickness to meet the
strength requirements of the thread. The development of the manufacturing process of thin-walled and thick-mounted cylinders, which can largely improve the storage efficiency of hydrogen,
has also attracted considerable attention. At present, a large number of domestic and international scholars have conducted research on gas cylinder-forming technology. The traditional
production process of gas cylinders involves a stamping process to form the upper and lower head parts, respectively, followed by welding. Music et al.2 pointed out that coreless necking and
spinning technology is widely used in the necking and forming process of axisymmetric rotary parts. Wong et al.3 showed that spinning technology, as the main production method for gas
cylinders, fundamentally solves the problems of discontinuity, low strength, brittle cracking and stress concentration of weld seams in the traditional welding production of gas cylinders,
and at the same time, spinning allows the production of parts with high mechanical properties and smooth surface finish. Wang et al.4 explored welded gas cylinders and found that residual
stresses due to weld seams had a significant impact on strength, fatigue life, and gas tightness, thereby posing a considerable safety hazard to the cylinders. Marini et al.5 claimed that
the flow-forming process can be extensively adopted for production of thin-walled, high-precision tubular products. Wong et al.6 used a two-step forming process consisting of
"bending" and flow forming. The method involved the flow of material along a mandrel to shape thin-walled cup-shaped parts, with the process employing two distinct profiles and
axial rolling. Jin et al.7 researched a monolithic manufacturing process for large-diameter seamless steel cylinders. The monolithic forming of cylinders could be achieved by means of cold
spinning of the cylindrical parts obtained by deep drawing. Through such means, the cylinder wall thickness uniformity was greatly improved and the dimensional accuracy was higher. By
investigating the relationship between the microstructure and mechanical properties of gas cylinders, Li et al.8 improved the mechanical properties of gas cylinders based on the
manufacturing process proposed by Jin et al.7. Zoghi and Fallahi Arezoodar9 manufactured pressure vessels using a necking spinning process, which essentially involved forming parts by
applying the required force and displacement to the rotating blank by means of one or more forming rolls. Wang et al.10 claimed that a combination of stamping and deep drawing using steel
plates could produce cylinder blanks as well as cups. Cylinder blanks prepared through such process had the advantages of uniform thickness, high strength, and light oxidation. Music et
al.11 reported that metal spinning can improve the surface finish and mechanical strength of molded parts compared with stamping. M.L. According to A et al.12, the utilization of numerical
modeling via finite element flow formulation can effectively facilitate the comprehension and prediction of different modes of deformation during the end forming of thin-walled tubes. Kuang
et al.13 used Ansys finite element software to establish a three-dimensional finite element model of offset circular tube header spinning. Using the model, during three-dimensional
non-axisymmetric spinning, the effects of metal flow, stress and strain distribution, spinning pressure, and different process parameters on the spinning results were investigated. In
analyzing the evolution of the material shape and thickness as well as the stress and strain distribution generated during the spinning process, Iguchi et al.14 used the dynamic explicit
code DYNA-3D to analyze the spinning manufacturing process of motor vehicle exhaust system components to further understand the failure mechanisms such as fracture and buckling during the
spinning process. The results could provide useful information for failure prediction during the actual spinning process. Xia et al.15 proposed a finite element model for non-axisymmetric
neck spinning using the finite element software MSC to obtain the transient Mises stress distribution in the contact zone between the roll and the billet, as well as the equivalent plastic
strain after spinning. Further, numerical and experimental studies on the thickness distribution of the spun workpiece were conducted. Such studies provided reasonable suggestions for
addressing the occurrence of excessively thin or thick wall thickness in the workpiece. Yao and Makoto16 conducted an experimental study on the near-axis spinning of tube ends and
investigated the effects on thickness strain, torsion angle, spinning force and surface finish of aluminum products with the parameters of spinning pitch and diameter reduction. Through
experiments and 3D finite element simulations, Yt et al.17 investigated the effect of neck length on crack generation during the spinning of SUS409, and determined the spinning conditions
under which cracking would not occur based on the calculated damage values. Huang et al.18 explored the thickness distribution and outer profile of the tube necking spinning process using
the finite element model established by the shell cell. that the research revealed that the thickness of the spun tube had a small amount of thickening with the spinning process, but there
was no explanation given for the deformation in the material thickness. Hamed et al.19 established a three-dimensional finite element model of circular tube spin forming to investigate the
strain distribution in different thickness layers of a tube during the spin-forming process. An increasing trend of strain in the middle layer thickness toward the free end was reported,
which indicated an increase in the tube wall thickness. Biplov Kumar Roy et al.20 investigated the change rule of geometry and thickness during necking and spinning forming of tubes by
combining experimental and numerical simulation. Response surface methodology (RSM), also known as regression design, is a statistical test for optimizing stochastic processes. The purpose
is to analyze the significance of each factor on the response value, to obtain a quantitative and reasonable regression equation, and to guide the engineering application. The central
composite design is a commonly used response surface analysis method, which is proposed on the basis of two-level full-factor partial experimental design. Zhou et al.21 optimized the process
parameters of AL7075 using second-order response surface methodology and experimentally verified the consistency between the predicted and experimental values of the response prediction
model. Asiltirk et al.22 analyzed the effects of feed rate, tool radius of arc, and back engagement on the surface roughness in the turning process as an evaluation metric using
Taguchi's methodology and response surface methodology. Rashmi et al.23,24 systematically investigated the effect of process parameters on milling force, surface roughness and energy
consumption during milling of AA6061 and determined the optimum machining parameters using response surface method and particle swarm optimization algorithm. Based on the aforementioned
studies, the gas cylinder forming process has attained greater maturity, with necking-spinning becoming the dominant production method. During the necking-spinning process, one end of the
billet with uniform wall thickness undergoes a certain degree of thickening in the mouth region. However, in the case of gas cylinders with considerable differences in wall thickness between
the mouth and the cavity, necking-spinning alone cannot achieve the ideal level of thickening. For this reason, the research purpose of this paper aims to solve the manufacturing process
difficulties of thin-walled and thick-mouthed seamless gas cylinders, and puts forward a boring-reducing neck spinning step-by-step scheme for the overall forming of thin-walled and
thick-mouthed seamless gas cylinders, and analyzes the neck-spinning of the cylinders by using finite-element analysis methods, and discusses the stress–strain distributions and geometrical
dimensional changes of the cylinders in the process of spinning; and at the same time, through the central composite experimental method involving a number of groups of experiments and
carrying out Numerical simulation is carried out to study the effects of spindle speed, friction block working angle and friction coefficient on the forming results, and the second-order
response surface method is used to establish the relationship between each process parameter and optimization variables, and then multi-objective optimization is carried out by the NSGA-II
algorithm to obtain the optimal solution of Pareto, and the simulation results are finally compared with the experimental results and analyzed. PROPOSAL FOR INTEGRAL FORMING OF THIN-WALLED
AND THICK-MOUTHED SEAMLESS GAS CYLINDERS ANALYSIS OF FORMING DIFFICULTIES The geometric structure of a thin-walled, thick-mouthed seamless gas cylinder is shown in Fig. 1a. The cylinder was
made of 2Cr13 martensitic stainless steel and had a total length of 300 mm, consisting of the cavity and the mouth. The cylinder had a wall thickness of 5 mm at both ends of the mouth and
0.7 mm in the cavity, resulting in a ratio of over 7 between the two wall thicknesses. The formability of thin-walled, thick-mouthed seamless gas cylinders could be analyzed from both
structural and material aspects. Firstly, from a structural formation perspective of gas cylinders, the ratio of wall thickness between the cylinder's mouth and the cavity exceeded
seven times, making it challenging to thicken the cylinder's mouth solely through a neck reduction spinning process. Secondly, there was a transition zone with gradually increasing wall
thickness between the bottle body and the mouth of the bottle, and the cylinder's cross-section may pose a potential hazard by leading to cracking during the formation process.
Finally, 2Cr13 stainless steel exhibits high yield strength and severe hardening during forming, which also results in higher deformation resistance and increases the risk of cracking. NEW
METHOD Considering the large variability of the wall thickness between the cavity and the mouth of a thin-walled, thick-mouthed seamless gas cylinder, solely using a single processing method
is not possible. As such, a step-by-step boring-necking-spin-forming scheme for gas cylinders was proposed in the present study. The flow chart of the new scheme is shown in Fig. 1b. The
two-step processing scheme was proposed for achieving desired wall thickness characteristics of thin-walled, thick-mouthed seamless gas cylinders. * In the first step, according to the
overall geometry of the gas cylinder, a straight seamless tube is selected as raw material and divided into two areas: the part of the mouth to be necked reduction and the part of the cavity
to be thinned. Boring tools are used to thin the inner wall of the cavity to the target wall thickness, so as to achieve a thin-walled cavity, while retaining the wall thickness in the
mouth area to meet the required thickness after forming. Notably, if there is a section with gradually decreasing wall thickness between the mouth and the cavity of the gas cylinder blanks,
a workpiece with thicker wall thickness at both ends and thinner wall thickness in the middle should be formed (as shown in Fig. 1c). * In the second step, both ends of the workpiece are
processed separately by means of the necking-spinning process to finally form a thin-walled, thick-mouthed seamless gas cylinder. The neck reduction spinning of gas cylinders can be formed
by multi-pass spinning with rolls or by spinning with friction blocks of a certain type of surface. The multi-pass spin-forming process for rolls is more complex and has high production
requirements. When the friction block is used for spinning, the entire surface of the friction block is in contact with the billet, the contact area is large, and the required forming force
is large. Under such conditions, one-time forming can be achieved, and there will be a better thickening effect on the mouth of the gas cylinder, which is suitable for small size workpieces.
Owing to the simplicity and efficiency of the process, the friction block spinning method was adopted for processing of the gas cylinder blanks in the present study. Figure 2a shows the
experimental setup for gas cylinder forming and a photo of the processed finished product. EXPERIMENTAL PROCEDURE The experimental process of neck-spinning of gas cylinders mainly includes
the following steps. * First, the workpiece is clamped in the rotary jig, and the workpiece is exposed to the necked portion using a positioning post at a certain distance. * Secondly, after
the mandrel rotates at high speed, the part of the workpiece to be necked and rotated is heated to the required temperature for forming using a high-frequency induction heater. * Finally,
the position of the friction block is adjusted according to the amount of necking needed to shape the gas cylinder's mouth, and the friction block is not rotated more than 90° around
the axis of rotation, resulting in the final formation of the gas cylinder's mouth. The described steps are repeated for the other end for neck spinning so as to finally form the
finished gas cylinder, as shown in Fig. 2b. FINITE ELEMENT MODELING To investigate the law of neck-spin forming of thin-walled thick nozzle seamless gas cylinders, a 3D finite element model
was developed based on ABAQUS software. Additionally, due to the symmetry of the workpiece itself, a 1/2 finite element model was established, as shown in Fig. 3a. In the model, the slewing
fixture and friction block were defined as a rigid body, the workpiece was set as a deformable body, and the unit body type was C3d8r. Initially, the overall mesh size was 2 mm × 1 mm, and
to ensure the forming accuracy and simulation efficiency, the mesh of the bottle mouth area was further refined to 1 mm × 1 mm, with the total number of meshes being 29,000. To improve the
calculation efficiency, the mass scaling factor was selected as 1000, and the ratio of kinetic energy to internal energy was less than 10% in more than 90% of the spinning process. As such,
the simulation was stable. The material properties are shown in Table 1 and these values were extracted from the study by Ge and Misra et al.25,26. These values were obtained through a
series of thermal compression tests performed on the material. The main process parameters used in the finite element model are given in Table 2, which were set to be consistent with the
actual working conditions.The mandrel speed, friction block working angle, and friction factor were temporarily set to 480r/min, 30° and 0.3, respectively. The schematic diagram of the
friction block working angle is shown in Fig. 3b. To facilitate the investigation into the deformation of the bottle mouth, the bottle mouth was divided into three deformation zones. Zone I
is the portion close to the bottle cavity, Zone III is the free end portion, and Zone II (the current contact area) is the portion between Zone I and III. Moreover, 25 mesh nodes were
selected as feature nodes along the axial direction on each of the three thickness layers. The positions of the feature nodes are shown in Fig. 3c. The comparison of the thickness
distribution obtained through the finite element model with the experimental results is depicted in Fig. 3d. Obviously, there is a good agree-ment between the simulation results and the
experimen-tal ones. ANALYSIS OF THE FORMING PROCESS THE STRESS AND STRAIN DISTRIBUTION DURING THE FORMING PROCESS Figure 4 shows the stress and strain distribution of the gas cylinder mouth
at different stages of the spinning process. The neck spinning process of a gas cylinder can be divided into three stages: Stage 1, the neck of the bottle mouth forming; Stage 2, the
cylinder mouth forming, also the neck shape trimming process; and Stage 3, the cylinder bending forming, but also the neck and mouth finishing process. Figure 4a shows that the equivalent
stress is predominantly distributed within the vicinity of the contact area, with the maximum stress concentration observed at the bottleneck region and the friction block's contact
position. Moreover, as the distance from the contact region increases, the equivalent force diminishes. Such local plastic deformation features of the friction block's neck rotation
spinning govern the observed phenomenon. In the spinning process, the contact area increases first, and then gradually decreases after reaching a certain peak, so the distribution range and
size of the equivalent force also show a trend of increasing and then decreasing. When the free end of the bottle mouth is no longer in contact with the friction block, the free end of the
bottle mouth is less constrained and its stress value becomes small. Figure 4b shows that the equivalent effect variation was distributed hierarchically along the axial direction, and the
size of the equivalent effect variation was equal in the same circumferential direction. With the gradual feeding of the friction block, the equivalent variation of the bottle mouth
increased until the end of spinning. The maximum equivalent effect variation value was located at the neck of the bottle curvature, which indicates that there was a large metal flow, with
the potential occurrence of cracks and buildup (as shown in Fig. S1). GEOMETRIC CHANGES (THICKNESS, ELONGATION, AND PROFILE) Xia et al.27 explored tube neck spinning and reported that the
thickness deviation between the actual dimensions of the workpiece and the nominal dimensions of the part must be satisfied. After the neck-spinning of the cylinder mouth, it is essential to
create threaded holes in its inner wall. However, due to the limited thickness of the mouth, the tolerance of the fit may be impacted. Therefore, it is crucial to ensure that the thickness
of the mouth meets the required wall thickness. Figure 5a shows the thickness and elongation of the bottleneck at different stages of the spinning process. An observation can be made the
elongation of the bottleneck increased from the beginning to the end of the spinning, while the wall thickness exhibited a process of decreasing and then increasing, and the final wall
thickness increased compared with the original wall thickness. Figure 5b shows the internal and external profiles of the bottle mouth during the rotary feeding of the friction block. An
observation can be made that as the spinning proceeded, the mouth of the bottle was gradually compressed to form a bend and the elongation band became more pronounced. The maximum wall
thickness value of the bottle mouth was 5.2 mm, which occurred at the inflection point, as shown in the lead in Fig. 5b. The minimum thickness value appeared at the edge of the free end of
the bottle mouth area, and there was an obvious thinning band at the free end of the bottle mouth. The wall thickness reduction at the edge of the free end became more serious as the
friction block was gradually rotated and fed. THE PLASTIC STRAIN DISTRIBUTION OF DIFFERENT THICKNESS LAYERS With the rotary feed of the friction block, under the high speed operation of the
spindle, with the help of air mold forming, the trajectory of the friction block movement so that it is always parallel to the surface of the workpiece, and keep the flat section and the
workpiece is in complete contact with the mouth of the bottle gradually compressed to form a bending, resulting in a large plastic deformation, so it is necessary to analyze the different
regions of the bottle mouth, Fig. 3c shows the different regions of the bottle neck. To analyze the variation of the bottle mouth thickness and elongation, the radial, circumferential,
axial, and equivalent plastic strains at the gas bottle mouth were examined, as shown in Fig. 6. In manufacturing certain pressure vessels using the tube spinning process, Zoghi and Fallahi
Arezoodar9, highlighted the existence of thickening, thinning, elongation, compression, bending, and shearing in the deformation zone during the neck spinning process of the bottle mouth. To
further investigate the strain state of the bottle mouth, the plastic strain distribution on each thickness layer in Fig. 3c was examined. Figure 7a–d depict the distribution of radial,
circumferential, axial, and equivalent variation for different thickness layers, respectively. From Fig. 6, an observation can be made that the strain distribution of each thickness layer
within Zones I and III was relatively uniform, and the differences between thickness layers were mainly reflected in Zone II (contact area). Figure 6a, b depict the distribution of the
radial strain, which is usually negative due to the necking behavior of the free end of the workpiece. An observation can be made that in Zone III, the strain state was reflected as a
compression strain. Figure 6b shows that positive strain existed in Zone II and the radial positive strain decreased from the outer layer to the inner layer. Combined with Fig. 7a, a
conclusion could be drawn that there was no significant difference in radial strain between the three thickness layers in the bottle mouth region, indicating that the decrease in diameter
during necking rotation was uniform between the thickness layers. Biglari28 reported a similar pattern of radial strain when rotating the plate with a hemispherical mandrel. Figure 6c, d
show the distribution of circumferential strains, which were consistently reflected as positive strains (tensile strains) in different zones between the three thickness layers. As seen in
Fig. 6d, the positive strain in Zone II decreased gradually along the outer layer to the inner layer. While directly subjected to the force of the friction block, the outer layer of the tube
was still undergoing positive circumferential strain (contact area) as a result of the constant volume during plastic deformation. In combination with Fig. 7b, the strain values of the
three thickness layers in Zone II gradually increased, implying an increase in thickness. The strain value of the thickness layer in Zone III gradually decreased, and the circumferential
strain at the edge of Zone III sharply decreased, which indicates a slight thinning in the thickness direction. Figure 6e, f show the distribution of the axial strain. The axial strains in
the three thickness layers in Zones I and II were reflected as negative strains (compressive strains), and the axial strains in Zone III were positive (tensile strains). In combination with
Fig. 7c, a conclusion could be drawn that the negative strain in Zone I gradually increased, which was caused by the local force of the friction block causing the bending of each thickness
layer by axial compression. The negative strain in Zone II gradually decreases and then rapidly changes direction. The axial strain in each thickness layer in Zone III was in a tensile state
after the direction of strain transition, and the strain value gradually increased towards the free end. The reduction of circumferential strain in Zone III was a result of volume
constancy, which compensated to some extent the positive axial strain of volume constancy. The bending due to compressive strain along the axial direction in Zone II and the axial extension
in Zone III together amplified the positive strain. Such findings are consistent with the deformation pattern of Zoghi and Fallahi Arezoodar9 in their research on neck spinning. The
excessive elongation of the outer metal at the rim in Zone III was caused by the unconstrained material at the rim of the bottle mouth, which is consistent with the behavior mentioned by
Hamed et al.19 in their research on tube spinning. According to the principle of volume invariance, due to the necking behavior of the bottle mouth, the metal flowed more along the axial
direction at the free end, then the flow in the circumferential direction had to be reduced. Moreover, the radial flow of the material gradually slowed down along the wall thickness
direction, as shown in Fig. 8. Therefore, the deformation of the gas cylinder mouth was mainly reflected in the axial elongation, while the thickening in the thickness direction was small.
Figure 6g, h show the distribution of the equivalent strain. As expected, there was a larger equivalent strain in Zone II. From Fig. 6h, an observation can be made that the values of the
equivalent strain in the same circumferential direction remained consistent, and the equivalent strain in Zone II gradually decreased from the outer to the inner layers. As can be seen in
Fig. 7d, in the axial direction, the value of the equivalent strain increased and then decreased. The outer layer strain had a larger growth rate. The maximum equivalent strain value was
located at the intersection of Zone III and Zone II (near node 14). Such findings indicate that there was a large amount of metal flow here and the strain gradient was more intense, which
could be attributed to the presence of curvature mutations in Zone II. RESPONSE SURFACE AND PARAMETRIC EQUATIONS FOR SPINNING PROCESS PARAMETER OPTIMIZATION The central composite design is a
commonly used method for response surface analysis, which is based on the two-level full-factor partial test design. It allows a finite set of sample points to be obtained in the design
domain, reducing the number of tests29. The central composite experimental design response surface methodology (RSM) is a combination of mathematical and statistical techniques used to
develop, improve and optimize processes. Using the central composite experimental method, multiple sets of experiments were designed and numerically simulated to analyze the effects of each
process parameter on the spinning results. Response surface methodology (RSM) is an optimization method that integrates experimental design and mathematical modeling, which effectively
reduces the number of experiments and allows for the examination of interactions between influencing factors, and the mapping relationship between each process parameter and the optimization
variables was established by means of the response surface method30. The resulting mapping relations were used to perform multi-objective optimization using the NSGA-II algorithm and obtain
the Pareto optimal solution, so as to optimize the process parameters and improve the quality of gas cylinders31. CENTRAL COMPOSITE EXPERIMENTAL DESIGN The maximum stress, and minimum
thickness were selected as the optimization objectives, and the three parameters of spindle speed, friction block working angle and friction coefficient were taken as the influencing factors
of the central composite experimental design to construct the response surface model. Such means could avoid the stress concentration during the forming process of the gas cylinder mouth
causing damage to the workpiece and further control the wall thickness of the mouth. Table 3 shows the range of values for each parameter. Table S1 shows the central composite experimental
parameters obtained by numerical simulation. The effects of each variable on the maximum stress and minimum wall thickness and their interactions were calculated using Design-Expert13
software. The relationship between the factors (independent variables) and each response (dependent variable) was modeled by fitting a second-order polynomial equation given by the following
equation. The effects of each variable on the maximum stress and minimum wall thickness and their interactions were calculated using Design-Expert13 software. The relationship between the
factors (independent variables) and each response (dependent variable) was modeled by fitting a second-order polynomial equation given by the following equation: $$Y={\alpha }_{0}+{\alpha
}_{1}{X}_{1}+{\alpha }_{2}{X}_{2}+{\alpha }_{3}{X}_{3}+{\alpha }_{12}{X}_{1}{X}_{2}+{\alpha }_{13}{X}_{1}{X}_{3}+{\alpha }_{23}{X}_{2}{X}_{3}+{\alpha }_{11}{{X}_{1}}^{2}+{\alpha
}_{22}{{X}_{2}}^{2}+{\alpha }_{33}{{X}_{3}}^{2}$$ in which \({X}_{1}\), \({X}_{2}\), and \({X}_{3}\) are the independent variables in the response surface model parameters, \({\alpha
}_{0}\), \({\alpha }_{1}\), \({\alpha }_{2}\), \({\alpha }_{3}\), \({\alpha }_{11}\), \({\alpha }_{12}\), \({\alpha }_{13}\), \({\alpha }_{22}\), \({\alpha }_{23}\), and \({\alpha }_{33}\)
are the regression coefficients, and \(Y\) is the response function. RESPONSE SURFACE RESULTS * (1) (1) Mapping relationship between a maximum stress and independent variables and response
surface model. The prediction model of the regression equation for the independent variables \({X}_{1}\), \({X}_{2}\), \({X}_{3}\) and the dependent variable \({Y}_{1}\) (maximum stress) is
shown in Eq. (1). Figure 9a shows the response surface of the maximum stress (\({Y}_{1}\)) with the variation of mandrel speed (\({x}_{1}\)) and friction block working angle (\({x}_{2}\)),
and Fig. 9b shows the response surface of the maximum stress (\({Y}_{1}\)) with the variation of friction block working angle (\({x}_{2}\)) and friction factor (\({x}_{3}\)).
$${Y}_{1}=-190.43150+0.990726{x}_{1}-0.3416{x}_{2}+1034.84125{x}_{3}-0.001684{x}_{1}{x}_{2}-0.151075{x}_{1}{x}_{3}-14.03200{x}_{2}{x}_{3}-0.000969{{x}_{1}}^{2}+0.067795{{x}_{2}}^{2}-1267.28750{{x}_{3}}^{2}$$
(1) Figure 9a reflects that the maximum stress decreased with the gradual increase of the spindle speed, which could be ascribed to the increased mobility of the material in the
circumferential and axial directions with the increase in spindle speed, rendering a decrease in the internal stress. Figure 9b reflects that the maximum stress increased with the gradual
increase in the friction factor, which was due to the increase in the resistance to the flow of the metal at the mouth of the bottle with the increase in the friction factor, and the
formation of larger stress. The effect of the working angle on the maximum stress was smaller. * (2) Mapping relationship between minimum thickness and independent variables and response
surface model The prediction model of the regression equation for the independent variables \({x}_{1}\), \({x}_{2}\), \({x}_{3}\), and the dependent variable \({Y}_{2}\) is shown in Eq. (2).
Figure 9c shows the response surface of the wall thickness value (\({Y}_{2}\)) with the variation of spindle speed (\({x}_{1}\)) and working angle (X2), and Fig. 9d shows the response
surface of the wall thickness value (\({Y}_{2}\)) with the variation of working angle (\({x}_{2}\)) and friction factor (\({x}_{3}\)).
$${Y}_{2}=4.13938+0.004942{x}_{1}-0.003977{x}_{2}-0.245801{x}_{3}-0.000051{x}_{1}{x}_{2}+0.001461{x}_{1}{x}_{3}+0.030487{x}_{2}{x}_{3}-4.62650E-06{{x}_{1}}^{2}+0.00024{{x}_{2}}^{2}-3.19512{{x}_{3}}^{2}$$
(2) Figure 9c, d reflect that the wall thickness value decreased with the increase in spindle speed and friction factor, which could be attributed to the increase in spindle speed and
friction factor, as well as the metal flow at the mouth of the bottle being enhanced and retaining less metal the thickness direction. of the increase in working angle tended to lead to a
decrease in the wall thickness value, which was because the increase in working angle led to the corresponding increase in contact area between the friction block and the mouth of the
bottle, thereby accelerating the material flow along the axial direction leading to a decrease in wall thickness. The normal plots of residuals for Max. Stress and Min. thickness as
demonstrated in Fig. 10a, c are approximately linear, indicating that the residuals are normally distributed and the fittings of regression equations for Max. Stress and Min. thickness are
reasonable. The points cluster near by the diagonal line are observed between the predicted and actual values of Max. Stress and Min. thickness in Fig. 10b, d. It suggests that the
established model is strongly appropriate for quantitatively describing the influence of various factors on the response of Max. Stress and Min. thickness. MULTI-OBJECTIVE OPTIMIZATION
According to the parameter equations discussed in "Response surface results" for the multi-objective optimization, the problem to be investigated could be described in the context
of the actual situation, where 400 ≤ \({x}_{1}\)≤600, 25 ≤ \({x}_{2}\)≤35, and 0.1 ≤ \({x}_{3}\)≤0.3. In order to prevent cylinder failure caused by excessive stress at the bottle mouth and
to ensure adequate thickness in the same Zone, it was necessary to simultaneously achieve the minimum value of the stress variable \({Y}_{1}\) and the maximum value of the thickness variable
\({Y}_{2}\), while keeping the independent variables \({x}_{1}\), \({x}_{2}\), and \({x}_{3}\) within the allowed variation range. Using Design Gateway for multi-objective optimization, the
regression equations were entered into the computer component, the initial values of the design variables \({x}_{1}\), \({x}_{2}\), and \({x}_{3}\) were entered, the upper and lower limits
of the variables were set, and the NSGA-II algorithm was selected to solve for all the Pareto solutions. In the Design Gateway window, the scatter plot between the maximum stress
\({Y}_{1}\), the thickness increment \({Y}_{2}\) and the independent variables can be observed as shown in Fig.S2a. In Fig.S2a, the Pareto solutions (i.e., Pareto fronts) are uniformly
distributed and located at the right fronts of all test points. Fig.S2b shows the distribution of the Pareto optimization solution set on the 3D coordinate axes. By analyzing the
distribution of the Pareto optimized solution set on the 3D coordinate axes, an observation can be made that the optimized solution set was mainly concentrated in a specific range, being a
significant reference value for the optimization and design of the target quantities. The optimal set of solutions provided by Design in the optimized solution set was the spindle speed of
490.4, a working angle of 34.439°, and a friction coefficient of 0.248. VERIFICATION OF FINITE ELEMENT SIMULATION Based on the experimental setup shown in Fig. 2a, spinning experiments were
performed on the workpiece. The process scheme was designed based on the process parameters determined by multi-objective optimization, discussed in "Multi-objective optimization".
Finite element simulations and experiments were performed on the process scheme, and the simulation results were compared with the experimental results, as shown in Fig. S3a, b. Due to the
high temperature of the spinning process, there were difficulties in performing real-time measurements, and it was not possible to investigate the deformation by microstructure due to the
reversion and recrystallization at high temperatures. As such, a decision was made to use the internal and external contours of the gas cylinder after molding, as well as the wall thickness
distribution to verify the simulation. Several researchers32,33 used a similar approach to perform a validation of the reliability of the model. A comparison of the wall thickness
distribution and the cross-sectional profile dimensions is shown in Fig. S3c, d, respectively. The maximum difference between the measured thickness from the workpiece and the predicted
thickness obtained from the simulation was about 8.02%, being within the allowable range. In the experiments and simulations, the minimum wall thicknesses were 5.2 and 5.014 mm,
respectively. In the actual working conditions, with heat dissipation and local heating of the supplemental heat device, the temperature in the deformation area of the billet was not
uniform, while the simulated temperature was constant in the simulation, thereby resulting in errors between the simulation and the actual machining. Such results indicate a good agreement
between the simulation and the experiment. In conclusion, the developed finite element model is reliable and the cross-sectional profile of the part is essentially consistent with the target
shape after spinning. The deviations of both wall thickness and profile are within acceptable limits. Therefore, the proposed scheme can be considered an effective method to achieve
thin-walled and thick-mouthed gas cylinder forming. CONCLUSIONS In the present study, a step-by-step boring-neck spinning process scheme was developed for forming thin-walled and
thick-mouthed seamless gas cylinders. Based on finite element analysis and experimental studies, the following conclusions were drawn: * (1) In the process of neck-spinning of gas cylinders,
there was a large stress in the mouth region and the axial elongation of the mouth was more significant, while there was a smaller thickening effect in the thickness direction. * (2) The
excessive elongation of the metal in the outer layer at the edge of the free end of the bottle neck (zone III) during neck-down spinning is caused by the unconfined material at the edge of
the bottle neck, which is consistent with the behavior mentioned by Hamed et al. researchers in their pipe spinning studies. In the case of constant volume, the necking behavior of the
bottle neck makes the metal flow more in the axial direction and then the flow in the circumferential direction has to be reduced. * (3) As the mandrel speed, friction block working angle,
and friction factor increased, the wall thickness of the bottle mouth decreased, while the maximum stress increased rapidly. Considering the influence of process parameters on the
aforementioned factors, a process parameter scheme was determined as follows: spindle speed of 490 r/min, friction block working angle of 34°, and friction coefficient of 0.25. * (4)
Experimental studies revealed that the desired thin-walled, thick-mouthed seamless gas cylinders can be formed using a step-by-step boring-necking spinning scheme and defined process
parameters. The minimum wall thickness of the gas cylinder mouth in the experiment was 5.2 mm, and the geometric profile accuracy as well as the wall thickness after forming was within the
error range. As such, the proposed step-by-step boring-neck-spinning scheme may become an advanced processing technology for forming thin-walled and thick-mouthed gas cylinder liners. The
accurate prediction of the inner and outer contour molding of the bottle mouth area can be studied in depth in the future to further overcome the molding barriers of thin-walled and
thick-mouthed gas cylinder liners. DATA AVAILABILITY During the current study, the datasets in Table1 Material parameters were available from the corresponding author upon reasonable
request. The rest of the data generated or analyzed are included in this publication and in the supplementary information. REFERENCES * Azeem, M. _et al._ Application of filament winding
technology in composite pressure vessels and challenges: A review. _J. Energy Storage_ 49, 103468 (2022). Article Google Scholar * Music, O., Allwood, J. M. & Kawai, K. A review of the
mechanics of metal spinning. _J. Mater. Process. Tech._ 210, 3–23 (2010). Article CAS Google Scholar * Wong, C., Dean, T. & Lin, J. A review of spinning, shear forming and flow
forming processes. _Int. J. Mach. Tools Manuf._ 43, 1419–1435 (2003). Article Google Scholar * Wang, C., Huang, S. & Xu, S. Optimization of girth welded joint in a high-pressure
hydrogen storage tank based on residual stress considerations. _Int. J. Hydrogen Energy_ 43(33), 16154–16168 (2018). Article CAS Google Scholar * Marini, D., Cunningham, D., Xirouchakis,
P. _et al_. _Flow Forming : A Review of Research Methodologies, Prediction models and their Applications_. Vol. 5 (2016). * Wong, C. C. _et al._ Cold rotary forming of thin-wall component
from flat-disc blank. _J. Mater. Process. Tech._ 208(1–3), 53–62 (2008). Article Google Scholar * Jin, W. _et al._ Reliability analysis of integral hot deep drawing and cold flow forming
process for large-diameter seamless steel gas cylinders. _Int. J. Adv. Manuf. Technol._ 97, 3 (2018). Article Google Scholar * Li, Y. _et al._ Microstructure and mechanical properties of
34CrMo4 steel for gas cylinders formed by hot drawing and flow forming. _Multidiscip. Digit. Publish. Inst._ 12(8), 1351 (2019). CAS Google Scholar * Zoghi, H. & Fallahi, A. A. Finite
element study of stress and strain state during hot tube necking process. _Proc. Inst. Mech. Eng. Part B J. Eng. Manuf._ 227(4), 551–564 (2013). Article Google Scholar * Wang, X., Jin, J.
& Lei, D. Review: State-of-the-art of stamping-forging process with sheet metal blank. _J. Harbin Inst. Technol. English Vers._ 24(3), 16 (2017). Google Scholar * Music, O., Allwood, J.
M. & Kawai, K. A review of the mechanics of metal spinning. _J. Mater. Process. Tech._ 210(1), 3–23 (2010). Article CAS Google Scholar * Alves, M. L. _et al._ End forming of
thin-walled tubes. _J. Mater. Process. Technol._ 177(1–3), 183–187 (2006). Article CAS Google Scholar * Kuang, W., Xia, Q. & Feng, R. Mechanism research for 3D non-axisymmetric
thin-wall tube offset spinning. _Chin. J. Mech. Eng._ 19(002), 217–222 (2006). Article Google Scholar * Iguchi, T., Yoshitake, A., Irie, T. _et al_. _Numerical Simulation and Development
of Tube Spinning Process for Exhaust System Components of Motor Vehicles_ (2004). * Xia, Q. _et al._ Finite element analysis and experimental investigation on deformation mechanism of
non-axisymmetric tube spinning. _Int. J. Adv. Manuf. Technol._ 59, 1–4 (2012). Article Google Scholar * Yao, J. & Makoto, M. An experimental study on paraxial spinning of one tube end.
_J. Mater. Process. Tech._ 128(1–3), 324–329 (2002). CAS Google Scholar * Yt, A. _et al._ Effects of neck length on occurrence of cracking in tube spinning. _Proc. Manuf._ 15, 1200–1206
(2018). Google Scholar * Huang, C. _et al._ Finite element analysis on neck-spinning process of tube at elevated temperature. _Int. J. Adv. Manuf. Technol._ 56, 9–12 (2011). Article Google
Scholar * Hamed, Z., Alireza, F. A. & Mojtaba, S. Enhanced finite element analysis of material deformation and strain distribution in spinning of 42CrMo steel tubes at elevated
temperature. _Mater. Des._ 47, 133 (2012). Google Scholar * Roy, B., Korkolis, Y., Arai, Y. & Araki, W. Experimental and numerical investigation of deformation characteristics during
tube spinning. _Int. J. Adv. Manuf. Technol._ 111, 1851–1867 (2020). Article Google Scholar * Zhou, M., Chen, Y. & Zhang, G. Force prediction and cutting-parameter optimization in
micro-milling Al7075-T6 based on response surface method. _Micromachines_ 11(8), 766 (2020). Article PubMed PubMed Central Google Scholar * Asilürk, I. & Neeli, S. Optimisation of
parameters affecting surface roughness of Co28Cr6Mo medical material during CNC lathe machining by using the Taguchi and RSM methods. _Measurement_ 78, 120–128 (2016). Article ADS Google
Scholar * Lmalghan, R. _et al._ Machining parameters optimization of AA6061 using response surface methodology and particle swarm optimization. _Int. J. Precis. Eng. Manuf._ 19(5), 695–704
(2018). Article Google Scholar * Malghan, R. L. _et al._ Application of particle swarm optimization and response surface methodology for machining parameters optimization of aluminium
matrix composites in milling operation. _J. Braz. Soc. Mech. Sci. Eng._ 39, 3541–3553 (2017). Article CAS Google Scholar * Ge, J. _et al._ Thermal-induced microstructural evolution and
defect distribution of wire-arc additive manufacturing 2Cr13 part: Numerical simulation and experimental characterization. _Appl. Therm. Eng._ 163, 114335 (2019). Article CAS Google
Scholar * Misra, R. D. _et al._ Dynamic recrystallization behavior of AISI 420 stainless steel under hot compression. _High Temp. Mater. Process._ 34(2), 155–161 (2015). ADS Google Scholar
* Xia, Q. X. _et al._ Numerical simulation and experimental research on the multi-pass neck-spinning of non-axisymmetric offset tube. _J. Mater. Process. Tech._ 206(1–3), 500–508 (2008).
Article CAS Google Scholar * Biglari, H. R. B. T. An experimental-numerical investigation of a metal spinning process. _Proc. Inst. Mech. Eng. Part C Mech. Eng. Sci._ 225(3), 1–11 (2011).
Google Scholar * Guo, C. _et al._ Milling optimization of GH4169 nickel-based superalloy under minimal quantity lubrication condition based on multi-objective particle swarm optimization
algorithm. _Int. J. Adv. Manuf. Technol._ 123, 1–12 (2022). Article Google Scholar * Liu, Y. _et al._ Multi-objective optimization of air dehumidification membrane module based on response
surface method and genetic algorithm. _Energy Rep._ 9, 2201–2212 (2023). Article Google Scholar * Deb, K. & Jain, H. An evolutionary many-objective optimization algorithm using
reference-point-based nondominated sorting approach, Part I: Solving problems with box constraints. _IEEE Trans. Evolut. Comput._ 18(4), 577–601 (2014). Article Google Scholar * Lexian, H.
& Dariani, B. M. Effect of roller nose radius and release angle on the forming quality of a hot-spinning process using a non-linear finite element shell analysis. _Proc. Inst. Mech.
Eng. Part B J. Eng. Manuf._ 223(6), 713–722 (2009). Article Google Scholar * Lexian, H. & Dariani, B. M. An analytical contact model for finite element analysis of tube spinning
process. _Proc. Inst. Mech. Eng. Part B J. Eng. Manuf._ 222(11), 1375 (2008). Article Google Scholar Download references ACKNOWLEDGEMENTS This research was funded by the Key Research and
Development Project of Shanxi Province (202202050201017). National Natural Science Foundation of China (52275358). AUTHOR INFORMATION AUTHORS AND AFFILIATIONS * School of Mechanical
Engineering, Taiyuan University of Science and Technology, Taiyuan, 030024, China Chen Wang, Haofei Yu, Wang Tian, Chunjiang Zhao, Lianyun Jiang & Qiaofeng Bai Authors * Chen Wang View
author publications You can also search for this author inPubMed Google Scholar * Haofei Yu View author publications You can also search for this author inPubMed Google Scholar * Wang Tian
View author publications You can also search for this author inPubMed Google Scholar * Chunjiang Zhao View author publications You can also search for this author inPubMed Google Scholar *
Lianyun Jiang View author publications You can also search for this author inPubMed Google Scholar * Qiaofeng Bai View author publications You can also search for this author inPubMed Google
Scholar CONTRIBUTIONS C.W. and H.Y. wrote the main manuscript text , W.T. , C.Z. , L.J.and Q.B. prepared figures 1-12. All authors reviewed the manuscript. CORRESPONDING AUTHOR
Correspondence to Chen Wang. ETHICS DECLARATIONS COMPETING INTERESTS The authors declare no competing interests. ADDITIONAL INFORMATION PUBLISHER'S NOTE Springer Nature remains neutral
with regard to jurisdictional claims in published maps and institutional affiliations. SUPPLEMENTARY INFORMATION SUPPLEMENTARY INFORMATION. RIGHTS AND PERMISSIONS OPEN ACCESS This article is
licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give
appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in
this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's
Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.
To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. Reprints and permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Wang, C., Yu, H., Tian, W. _et al._ Research
on the integral forming process of thin walled and thick mouth seamless gas cylinders. _Sci Rep_ 13, 17021 (2023). https://doi.org/10.1038/s41598-023-44377-z Download citation * Received: 08
August 2023 * Accepted: 07 October 2023 * Published: 09 October 2023 * DOI: https://doi.org/10.1038/s41598-023-44377-z SHARE THIS ARTICLE Anyone you share the following link with will be
able to read this content: Get shareable link Sorry, a shareable link is not currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing
initiative