Thermoelasticity Coupling Transient Response of Composite-Material Direction Pipe

Based on the positive scheme method, the thermal load of the jet flow in an inner composite-material direction pipe is obtained, and the thermoelasticity coupling transient response is investigated. *e positive scheme method with second-order accuracy is extended for solving the axisymmetric Euler equations, and the supersonic axisymmetric jet flow over a missile afterbody containing jet exhaust is simulated. *e correctness of the development for the positive scheme method is verified. With the developed positive scheme method used to simulate the jet flow in the inner direction pipe, the thermal load is obtained. *e thermoelasticity coupling finite element model of the composite-material direction pipe is established, and the stress response under dynamic pressure, unsteady temperature, and coupling state is obtained. Results show that, at the beginning of engine ignition, the effect of dynamic pressure and temperature field on the coupling stress is basically the same, and after that, the contribution of the temperature field to the coupling stress increases, and the thermal stress is the main factor affecting the strength of the composite-material direction pipe.


Introduction
e composite-material direction pipe is an important component of a multiple launch rocket system.During the launch, the transient response of the direction pipe is one of the key factors affecting the initial disturbance and firing intensity of the launch system [1].When the rocket is launched, the inner pipe is subjected to dynamic pressure, an unsteady temperature field, and boundary constraint.Under the thermal load of the jet flow, the transient response of the direction pipe shows a significant thermoelasticity coupling effect.Obtaining the thermal load of the inner jet flow of the direction pipe is the precondition for the study of its thermoelasticity coupling transient response.However, it is difficult to test the obtainment of the thermal load because the time of the formation and development of the rocket gas jet flow is very short.e highresolution numerical scheme has become an important means to calculate the structure of wave systems, pressure fields, and temperature fields [2].Demkowicz and Gopalakrishnan [3] proposed a Runge-Kutta discontinuous Galerkin method that had the same high-order accuracy and high resolution in the spatial and time dimensions.Chang [4][5][6][7][8] carried out a numerical calculation for the unsteady-state underexpanded jet flow with the CE/SE method and captured the complex wave structures such as shear layer, expansion waves, incident shock, and mach shock disk.Liu and Lax [9][10][11][12] proposed a positive scheme method for computational fluid mechanics which was used to solve hyperbolic conservation laws.At present, there are few reports regarding the positive scheme method in the calculation of gas jet flow.For the thermal-elasticity coupling problem of composite material, Yang et al. [13,14] performed a study for the three-dimensional thermalelasticity coupling problem of fiber-reinforced composite material under thermal-mechanical load with the finite element method.Palumbo et al. and Kylili et al. [15,16] analyzed the thermal-mechanical stress for the failure mechanism of composite material and provided an analysis method for the thermal-mechanical stress of glass fiberreinforced composite material.So far, the effects of the pressure field and the temperature field on the coupling state have seldom been studied.e steady-state response of composite materials has been the primary focus, and the study of the transient response at the millisecond level has been even less frequent.
In this study, the thermoelasticity coupling transient response of a composite-material direction pipe under the thermal load is researched.On the basis of the work of Fazio and Jannelli [10], the positive scheme method with secondorder accuracy is extended for solving the axisymmetric Euler equations, the calculation code is programmed and applied to simulate the supersonic axisymmetric jet flow over a missile afterbody containing jet exhaust, and the correctness of the development for positive scheme method is verified.With the developed positive scheme method used in the numerical simulation for the gas jet in the direction pipe, the unsteady process of wave structure, pressure field, and temperature field is revealed.e thermoelasticity finite element model of the composite-material direction pipe is established, and the thermoelasticity coupling transient response of the direction pipe is calculated, which provides reference for the structural design.

Development of Positive Scheme Method
e numerical calculation of the oblique shock regular reflection problem indicates that the positive scheme method is of high calculation accuracy and perfect resolution [9].To apply the positive scheme method to simulate the jet flow of a composite-material direction pipe, the method is developed to solve the axis-symmetric Euler equations.

Governing Equations.
For supersonic jet flow, the governing equations can be expressed as inviscid, axisymmetric Euler equations in the cylindrical coordinate system [9] as where T is the time variable, ρ is the density, u and v are the velocity components in the x and y directions, respectively, E � ρe + ρ(u 2 + v 2 )/2 represents the total energy in unit volume, e is the internal energy in unit mass, p is the pressure for the ideal gas, the state equation is p � (c − 1)ρe, and c is the specific heat capacity.

Space Discretization.
e finite difference equations of the axisymmetric governing Equation (1) discretized using the positive schemes are given in the semidiscrete form as where Δx and Δy, respectively, are the spatial steps in the x and y axes.e construction process of positive scheme numerical flux, taking F i+1/2,j , for example, is briefly described as follows.Other details of the calculation technique can be found in Reference [9].e numerical flux F i+1/2,j in Equation (2) was constructed by mixing a centered difference flux F c i+1/2,j and an upwind flux F up i+1/2,j .Introducing the limiter L 0 � Rdiag(ϕ 0 (θ k ))R −1 , ϕ 0 (θ) satisfies [17,18] e numerical flux could be constructed as where is the absolute value of A, A � (zF(U)/zU), and |Λ| � diag(|λ k |). e numerical flux of the least dissipative scheme is given by is the minmod flux limiter that satisfies [17,19] e numerical flux can be constructed as where ).Similarly, the numerical flux of the more dissipative scheme is expressed as Finally, by combining Equation (5) the least dissipative scheme and Equation (8) the more dissipative scheme, the numerical flux of positive schemes can be written as 2 Shock and Vibration under the condition (Δt/Δx)(α max where Δt is the time step.e constants α and β satisfy 0 ≤ α ≤ 1 and α + β ≥ 1.

Runge-Kutta Time Discretization.
To match the accuracy with the space variables, we constructed the time discretization with the accurate second-order Runge-Kutta method [20] as 2.4.Verification of Supersonic Jet Flow.e positive scheme method was used to solve the axis-symmetric Euler equations, and the supersonic axis-symmetric flow over a missile afterbody was simulated.Half of the calculation domain cut by any meridian plane through the central axis of the jet is shown in Figure 1. e radial direction was 1.5, and the axial direction was 3.5.e calculation domain was dispersed by a uniform grid with Δx � Δy � 1/40, the nozzle surfaces EF and FG used a mirror reflection boundary condition, the upper boundary DE used a supersonic flow condition, the outer boundary CD used a simple wave condition, the central axis of the flow AB used a symmetric condition, the lower boundary BC used an extrapolation boundary condition, and the other calculation domains used the initial condition of supersonic free flow.As shown in Figure 2, the exhaust jet was compressed by the supersonic freestream around the lip of the nozzle exit, also the Mach disc structure almost disappeared near the jet upstream.Meanwhile, two oblique shock waves, a contact discontinuity jet shock wave, and a bottle shock wave were formed.Behind the first oblique shock wave are the expansion waves around the lip of the nozzle exit, and under the second oblique shock wave is a jet shock wave.ese expansion waves propagate across the jet flow and are reflected back into the jet as compression waves from the jet boundary, thus forming an incident shock.Later, a reflected shock wave appears when the incident shock encounters the jet centerline and intersects with jet shock wave.Overall, the several primary flow characteristics between the numerical results and the experimental picture [21] obtained under the same flow conditions as shown in Figure 2(c) agreed reasonably well, indicating that the numerical results obtained by the positive scheme method was reliable.

Positive Scheme Numerical Study for Gas
Jet in a Direction Pipe e nozzle outlet pressure of the rocket engine p c was 5.058 atm, the nozzle outlet density ρ c was 0.81 kg/m 3 , the Mach number of the nozzle exit Ma c was 2.63, the nozzle exit temperature T c was 1,780 K, the nozzle outlet diameter d c was 0.0613 m, and the inner diameter of the direction pipe D was 0.07 m.Setting the characteristic length L 0 � d c , the characteristic density ρ 0 � ρ c , and the characteristic speed U 0 � a c � 881.57m/s as three independent reference quantities, the dimensionless treatment of the above parameters was conducted.
Half of the calculation domain cut by any meridian plane that passes through the nozzle jet is shown in Figure 3, and the nondimensional calculation domain was [0 : 53.8] × [0 : 1.14].e central axis AB of the direction pipe, the inner surface CD of the direction pipe, and the surface EF of the nozzle solid wall all used the mirror reflection boundary condition, the exit cone boundary AF of the nozzle used the cone inflow condition, the lower boundary BC used the extrapolation boundary condition, the upper boundary DE used the static atmosphere condition, and the static atmosphere condition was assigned to the initial jet flow.
e density contour of the inner jet flow at different times (dimensionless time) is shown in Figure 4, which involves 1/4 the length of the pipe.As indicated from Figures 4(a)-4(f ), after the unexpanded gas flowed through the nozzle into the atmospheric environment, the initial shock wave was formed by mutual compression together with the air inside the pipe.Owing to the intersection of the boundary of the jet flow with the direction pipe's inner surface, when the jet flow collided with the direction pipe's wall surface during its propagation process, a back jet flow together with two reflection shock waves were generated.As time developed, the two reflection shock waves intersected at the center axial of the direction pipe, while new reflection shock waves were penetrated and formed on the corresponding directional surface.Along with repeated intersecting and penetrating, the first intersection at the center axial was about 5.5 times the pipe's diameter away from the outer pipe surface (shown in Figure 4(f )). e flow characteristics of the pipe were consistent with the calculation Shock and Vibration results of Reference [22], which demonstrates that the results of the positive scheme method were effective and reliable.e pressure and temperature curves of the inner surface of the direction pipe at t � 7 ms are shown in Figure 4.As shown in Figure 5(a), if the direction pipe is not long enough, the impinging shock intersects and reflects periodically inside the direction pipe, and the distance of the two pressure spikes caused by the impinging shock is approximately the same.At a certain length of the direction pipe, the structure of the impinging shock had nothing to do with the length of the direction pipe.e strength of the impinging shock decreased along the direction pipe, and the variation of pressure curves matched well with those of the relevant literature [23].

Calculation Model of the Direction Pipe.
e structure of the composite-material direction pipe using a cylindrical and spiral design is shown in Figure 6. e front, middle, and back bourrelets were connected with the three clapboards of the launcher.e direction pipe plays a role in fixing the rockets in the state of marching, giving the rocket initial direction and a reasonable muzzle speed at launch time.
ere exists spiral guide groove in the inner surface of the direction pipe; by acting with the directional knob of the rocket, the rocket produces a rotational motion around its axis, which provides a suitable off-pipe rotational speed for the rocket.
e finite element model of the composite-material direction pipe was established using commercial software ABAQUS [24].In the PART module, the geometric model of the direction pipe was established according to the actual structural parameters.In the PROPERTY module, the composite shell section property was assigned to the geometric model; the composite material used the epoxy resin as the matrix and glass fiber as the reinforced material; the composite modes adopted the circumferential stacking, cross stacking, and axial stacking; the ply order was [(90 ∘ / ± 45/0 ∘ ) 2 ]; the single layer thickness was 0.25 mm; the density was ρ � 1.8 × 10 3 kg/m 3 ; the elasticity modulus was In the STEP module, the thermoelasticity coupling algorithm was adopted.In the MESH module, in order to reflect the structure characteristic of the composite-material direction pipe, the 8-node reduced integration brick element was adopted; there were 1762 elements in all.In the LOAD module, the front bourrelet, middle bourrelet, and back bourrelet were constrained.e finite element model of the composite-material direction pipe is shown in Figure 7.

Boundary Condition.
As indicated in Figure 7, three bourrelets were connected with clapboards; therefore, the bourrelets of the direction pipe should meet the displacement boundary conditions where x 0i , y 0i , and z 0i (i � 1, 2, 3) represent, respectively, the three coordinate values of the elements of the bourrelets.e thermal stress effect is generated when the composite material is forced by the temperature field of the gas jet; therefore, the finite element model of the composite-material direction pipe should meet the temperature boundary condition.e environment temperature of the outer surface of the direction pipe was T(x, y, z, t)| S out � 300 K.
e gas jet temperature field of the inner surface of the direction pipe was T(x, y, z, t)| S inner � T(z), where the unsteady temperature e experiment indicated that during the launching process, the pressure curve measured at a certain location was the same as that measured in each location along the inner surface of the direction pipe at a certain moment [23].Based on this conclusion, by transforming the gas jet pressure curve and the temperature distribution curve from the spatial domain to the time domain, the pressure or temperature that acts on the elements or nodes of the direction pipe can be obtained.As pressure is a type of surface load, the pressure should be imposed on the elements of the direction pipe, while the temperature should be imposed on the nodes.

Thermoelasticity Coupling Analysis of
Composite-Material Direction Pipe

Middle bourrelet Front bourrelet
Back bourrelet Guide groove

Shock and Vibration
is the end of the pipe and the right side is the nozzle) under dynamic pressure are shown in Figures 8(a)-8(g).As indicated from Figure 8(a), in the initial period of the launching process, although the engine ignited, the gas jet did not impinge with the direction pipe.e maximum stress occurred near the bourrelets by the action of gravity.From Figures 8(b)-8(g), the maximum stress occurred in the locations where the gas jet impinged with the inner surface of the direction pipe to form an impinging shock.When the rocket moved in the direction pipe, the maximum stress area moved quickly from the end of the direction pipe to the nozzle, and the maximum stress area nearly departed from the direction pipe after 80 ms as shown in Figure 8(g).
Four representative points in the inner surface of the direction pipe were selected along the axial direction, and the stress responses were compared.e distribution of characteristic points A, B, C, and D are shown in Figure 7. e circumferential stress and von Mises stress temporal change curves of the characteristic points are shown in Figures 9(a) and 9(b).As shown in Figure 9(a), the circumferential stress is the tensile stress, along with the movement of the rocket.When the impinging shock reached the point, a phase step occurred for the circumferential stress, and the curve assumed a pulse shape similar to that of the pressure curve.In view of Figure 9(b), a similar distribution law was obtained.As the locations of the characteristic points were different, the peak values of the curves differed.e peak value of the stress of the composite-material direction pipe under dynamic pressure was not large enough (the maximum was about 35 MPa), the stress of the characteristic points increased quickly when the impinging shock reached a certain location, and a pulse shape formed and then decreased quickly.

Response Characteristics under Unsteady Temperature
Field.Under the unsteady temperature field, the von Mises stress nephograms of the inner surface of the composite-material direction pipe at different times are shown in Figures 10(a)-10(g).From Figure 10, when the rocket passed through the direction pipe, a large stress area occurred because of the temperature field of the gas jet, and the stress area moved from the end of the direction pipe to the nozzle.
e distribution laws of stress nephograms differed from those of the stress nephograms under dynamic pressure (Figure 8).In Figure 8, the stress area moved along with the rocket, while in Figure 10, the stress area moved ahead with the temperature field and was located mainly in the middle of the first and second bourrelets and in the middle of the second and third bourrelets.A high-temperature gradient existed between inner and outer surfaces of the direction pipe because the temperature field of more than 1,000 °C acted on the inner surface of the direction pipe.
e direction pipe expanded under the temperature load, but there were restraints from the clapboards in the bourrelets, and the areas near the bourrelets could not expand freely, thus resulting in a high thermal stress.
e circumferential stress and von Mises stress temporal change curves of the characteristic points are shown in Figures 11(a) and 11(b).As shown in Figure 11(a), the circumferential stress of the characteristic points was tensile stress, and the peak values in characteristic points A and D were larger than those of characteristic points B and C. From the von Mises stress curve in Figure 11(b), when the combined influence of the stress components is taken into consideration, the peak values of characteristic points B and C were slightly larger than those of A and D. e stress distribution along the axial direction was inhomogeneous under unsteady temperature with the pulse-shaped stress response curve.
e stress caused by the temperature field was larger than that of the dynamic pressure, indicating that the thermal stress had a greater influence on the strength of the direction pipe.under dynamic pressure, the temperature eld, and the coupling state are shown in Figures 13 and 14.As shown in Figure 13(a), the dynamic pressure stress, the thermal stress, and the coupling stress all increased and decreased quickly.After that, the stress curves did not uctuate wildly, and the peak values of the stress caused by dynamic pressure and the temperature eld were similar, as was the shape of the ascending and the descending segment of the coupling stress curve and the thermal stress curve.As shown in Figure 13(b), the consistent degree was clearer, which indicates that the thermal stress signi cantly a ected the coupling stress.As shown in Figure 14(a), the dynamic stress, thermal stress, and coupling stress quickly increased  Shock and Vibration to the peak value, and the gradient change was large.Slight uctuations occurred and then tended to be steady after the curves decreased.e uctuations were more regular after the rst peak value than those of the corresponding curve of characteristic point A and the amplitude was larger.Similarly, the coupling stress curve and thermal stress curve matched well in the uctuation segment where the rst peak value increased and decreased.
e circumferential stress under the temperature eld and the coupling load was the compressive stress in the initial time when the rst wave crest rose.In a very short time, the stress changed to tensile stress, and the circumferential stress under the dynamic pressure was always the tensile stress.As shown in

Conclusion
(1) e numerical results obtained by the developed positive scheme method matched well with the experimental results.e correctness of the developed positive scheme method was veri ed.(2) e jet ow of the direction pipe numerically calculated by the developed positive scheme method    Shock and Vibration shows that the positive scheme method has a high ability for capturing complex gas jet structures.e ow characteristics and pressure curves of the pipe were consistent with the results of the reference, which further shows that the development of the positive scheme method was successful.
(3) For the stress response of the pipe under dynamic pressure, the stress curve shows a steep peak as the impinging shock wave arrives, then the stress decreases quickly, and the stress of the composite-material direction pipe under dynamic pressure was small.(4) For the stress response of the pipe under the temperature eld, the thermal stress along the axial direction was uniform, the stress curve presented a pulse shape; the peak value of the thermal stress was large.(5) For the stress response of the pipe under dynamic pressure and a temperature eld, at the initial time when the engine ignites, the stress caused by the Shock and Vibration dynamic pressure and temperature field had the same influence on the coupling stress of the direction pipe, and then the contribution to the coupling stress from the temperature field increased.

Figure 2 (
a) shows the density contours, Figure 2(b) the pressure contours, and Figure 2(c) the experimental picture under the same flow condition of a supersonic gas jet field.

Figure 14 (
Figure 14(b), the von Mises stress caused by the temperature eld and the coupling load was almost the same except for the wave crest, and the leading role of thermal stress in the coupling stress is clearly shown.

Figure 14 :Figure 13 :
Figure 14: (a) Circumferential stress and (b) von Mises stress of characteristic point B.