Abstract
Under the service conditions, steel pipelines coated with the thermal protection system are subjected to cyclic loadings of axial tension and hydrostatic pressure. The finite element method generally used to simulate the behavior of composite structures under these loadings allows us to estimate the stresses generated in the system and to conclude on several origins of damage. However, for the framework of displacement or deformation analyses in such multilayer systems, these calculations do not allow a better prediction of their behavior. The methods used do not take sufficiently into account the characteristics of the different coating materials to predict their response in service conditions under cyclic loading. In this paper, we consider the viscosity of the thermoplastic materials used for the five layers coating system. Finite element calculations allow us to observe the areas of highest stress concentration at the interface with the steel pipe. Simulations allowed us to observe that the applied loads lead to increases in residual deformation in the thermoplastic matrix composite material. Cyclic tensile loading causes cracks in the matrix of the syntactic foam material. The study carried out here makes it possible to justify the origin of the failure mechanism in the composite material at the time of the installation of the pipelines which could limit the duration of their use in an offshore environment. The tensile failure of the syntactic foam considered as the polypropylene matrix composite material on which cyclic loads have been applied, is due to the stress level at a given temperature.
Introduction
The use of offshore pipelines requires a certain level of performance in the face of mechanical loads that may arise from installation processes and during the fluid transport phase at sea surface level. 1 High mechanical loads and strong thermal gradients can cause a situation of asphalt formation and paraffin deposit which may limit for this purpose the flow rate of the fluid along the steel pipe. In order to maintain the oil production, thermal insulation systems have been proposed. However the type of system consisting of five layers has been widely used by the companies. 1 This system consists of adhesive Epoxy, adhesive polypropylene, solid polypropylene and syntactic foam with a polypropylene matrix (GSPP).1–3 This system is recognized as bringing together the different advantages that a single material could not have had. 2 During the manufacturing of thermally insulated pipelines, the steel pipe coating process is carried out by the heating, and then water-cooling of the whole structure.2,3 At the end of this manufacturing, residual stresses are developed in the thickness of materials polymers. 3 Commissioning is a process during which the pipelines held in any position must support the various loads applied.4,5 under service conditions, pipelines must be subjected to hydrostatic pressure type loads.3,6–9 The experimental analysis carried out in Rizzi et al. 10 made it possible to underline the plastic behavior of the syntactic epoxy matrix foam under various loads. A nonlinear Drucker-Prager model was then proposed for the numerical simulation of the mechanical behavior of the material. 10 Attention is paid to the stress states developed in the pipes in installations. The use of long life pipes requires the inclusion of a time-dependent overall response for the multilayer thermal insulation system. The objective of our work is therefore to improve the tools for optimizing mechanical performance during the design and sizing of passive pipeline protection systems used in deep water. This study consists in numerically evaluating the time-dependent mechanical behavior of the materials constituting the multilayered pipe system under the influence of the various loads applied to the assembly. Within the framework of this modeling, the loading conditions were represented by the application of a tension and external pressure at the external borders. The characteristics of the constituent materials derived from the experimental analysis are obtained in the literature.1,3,11–14
Materials and geometry of the structure
The coating system used to provide the thermal insulation function of the steel pipeline in the offshore environment can be either three-layer or five-layer.
15
We are interested in the five-layer system which also has the advantage of good resistance to hydrostatic pressure by inserting a thickness of syntactic foam material.2,3 Figure 1 shows the geometry and constituent materials of the five-layer surfacing system. According to the geometry of this structure, the five-layer system is based on the first three materials of the three-layer system that has been in use for several years. Tables 1 and 2 show the thicknesses and the elastic properties of the materials used in both types of systems.

The five and three layers systems of thermal insulation.
Materials thickness. 2
Polymers and steel pipe properties at the temperature of 20°C. 2
During the numerical calculations two types of loads are used in the static and time domain. Tensile loading conditions consist in applying a force to the outer boundaries of the structure using a triangular mesh and the displacements are blocked in the direction perpendicular to the loading direction. The pressure loading conditions remain the same but the loads are applied to all the external boundaries of the structure. Figure 2(a) to (d) shows the different constituents of pipe and loading conditions used in the finite element code of the Comsol Multiphysics Software.

(a) 3D geometry of the five layers systems under cyclic tensile and pressure loadings; (b) 2D axisymmetric geometry of the three layers system under pressure loadings; (c) 2D axisymmetric geometry of the assembly of the polypropylene materials and glass syntactic polypropylene foam under tension force, (d) Pressure loadings conditions applied to the three layers system
Methods
Overview
The glass syntactic polypropylene foam is a thermoplastic composite material used in the lining of Offshore pipeline to provide adequate performance under severe environmental conditions. The characterization of the behavior of this material is therefore carried out to optimize the tensile and pressure strength of thermal insulation system in installation and immersion conditions. Creep loading tests are therefore well suited to study the behavior of this material. The temperature and the stress must be analyzed simultaneously in the permanent deformation of the material.
Constitutive equations of the models
The overall resistance of a thermal insulation system will be analyzed in this manuscript by studying the individual behavior of the constituent materials. The metal steel tube is the tube in which the fluid circulates. This tube is protected by a system of five layers of materials mainly made of polypropylene material (5LPP).2,3,12 The time-dependent behavior patterns selected for each polymer material are shown below. We assume that each material has a response that takes the following general form:
The steel material has an elastic behavior (
The first layer of the coating system is epoxy and acts as an adhesive between the metal tube and the thermal insulation materials. This layer is only a few millimeters thick enough to fulfill its role.12,15 The viscous strain of this material will not also be considered here owing to the fact that its mechanical response consists mainly of the instantaneous strain under the level of applied creep stress. 15
Then, three layers of polypropylene materials are deposited on the steel tube. 2 The first layer acts as an adhesive to improve the adhesion of materials with the steel tube. The second layer acts as a solid polymer to protect the primary insulation material Glass Syntactic polypropylene from impact loading. The third layer protects all of the constituent materials against external corrosion.1,3,12 The viscoelastic character of the polymer material will be simulated with the generalized Maxwell model.3,15
The model connects the different deviatoric parts of the strain and the stress.
11
The deviatoric stress
The layer of composite material, which acts as the main thermal insulation material, has a greater thickness.2,3,12 Its resistance is at the center of studies of the overall system due to the particular properties of the material, namely the reduction in the weight of the structure, the maintenance of the temperature of the hot fluid, the guarantee of the longevity of the system and its buoyancy at the sea.
2
Recently a spectral model has been presented for polymeric and composite materials.16,17 The model is used for the simulation of the syntactic foam material response.
3
This spectral model gives the viscous elementary strain
We can operate this model using the details presented by Treasurer et al.
19
By solving the differential equation (3), we will introduce the total linear viscoelastic deformation
for
for
The normalization of the weighting coefficients is given by:
The experimental analyses of the creep behavior of composite materials were carried out in the literature.16,18 They made it possible to identify the irreversible nature of the latter. The model proposed to take into account the irreversible mechanisms in the deformation of polymer and composite materials is that of Zapas and Crissmann. 21 This model will permit to link the response of the Glass syntactic polypropylene material to the three variables time, temperature and stress by Eqn. (9):
where
Results and discussions
In this part, we perform analysis of the constituent materials response of the thermal insulation coating system. The objective is to determine their behavior in service when the axial tension and the external pressure are applied.22–24
Two organic coatings under tensile loads
The installation and immersion phases of the pipes can be carried out by the J-shaped installation method. 25 This method has the advantage of reducing the effects of certain installation constraints.3,22,24,25 During the installation process, axial tension loads can occur in the thicknesses of the different materials constituting these structures.4,25 We consider a five-layer system of Glass Syntactic Polypropylene (see Figure 1)2,3,12 and a three layers system of polypropylene (see Figure 1)26,27 under the effects of loadings in axial tensile load. The experiments details of quasi-static tensile tests are presented by Vestrum et al.13,14 The response of the pipeline samples was obtained using a deformation hydraulic piston. An Instron general purpose machine with the ability to increase the force was used to perform all tests. The force was measured by a load cell and the displacement was recorded by two sets.13,14
The geometric characteristics of the two types of thermal insulation system considered here have been presented in Aris-Brossou et al.,
12
Joliff et al.
15
and Bouchonneau et al.1,2 During the axial tension experiments performed on two pipe sections, the thermal insulation system showed an elastic behavior even under the two forces of 300 kN and 600 kN. The axial tension experiments performed on the materials use machines whose forces does not usually exceed 300 kN. In order to note the consequences of such a critical force on the structure we have applied it in the numerical study. The geometry of the system is chosen in 3D dimension (see Figure 2(a)) to be able to study the mechanical strength under the application of the force of 300 kN.
4
Axial tension forces are considered only in the longitudinal direction

Displacements distribution in the five and three layers system of thermal insulation.
Three layers coating under external pressure
We study the effect of external pressure applied to the polypropylene thermal insulation structure. We choose three pressure levels (7, 15 and 30 MPa) representing any depth of immersion.1,28,29 The distributions of stresses and displacements will be obtained in the thickness of the thermal insulation system in Polypropylene PP material. The experimental details of the hydrostatic compression tests are presented by Bouchonneau et al.1,2 A high-pressure vessel was used for the prototype testing. The pressure and water temperature were regulated and monitored inside the pressure tank. Figure 2(b) presents the 2D axisymmetric geometry and the constituent materials of the pipe structure type. The mesh used for the computation by finite elements is of triangular type. 1 The loading conditions are represented at the external borders of the structure for all the levels of external pressure (see Figure 2(d)).
Stress distribution
The quantification of the stresses is carried out in the pipe structure in order to evaluate the resistance of the materials in their elastic domain during the application of the external pressure. 30 The Von Mises criterion is used to identify the failure mechanisms of these materials.
Figure 4 shows the maximum stress distributions of Von Mises in the polypropylene (PP) coating system. The maximum Von Mises stresses are obtained for different levels of external pressure applied to the pipe structure. We can observe that for the pressure of 7 MPa, the maximum level of the Von Mises stress is 46,243 MPa, then for the pressure of 15 MPa, the maximum level is 99,914 MPa and finally for the pressure of 30 MPa this maximum is 197.05 MPa. Observations on the various Von Mises maximum stresses show that these are slightly above the breaking stress of the polypropylene material at low pressure (7 MPa) and are rather very high compared to the limit stress of this material for higher pressures (15 and 30 MPa). The response of the polypropylene material under the external pressure loading may involve a domain of plasticity just as in the case already observed on the polyethylene material 15 and the steel material. 28 We can say that the presence of plastic deformations of the material can be more or less marked according to the level of external pressure.1,3,6,8 Within the framework of a study of the viscoelastic and plastic properties of the primary epoxy and polyethylene materials, the influence of the coating thickness was observed during the manufacturing process of the passive protection system. 15 This justifies in the context of this study the dependence of the response of the polypropylene material on the level of external pressure. The viscous properties of the polypropylene material are furthermore temperature sensitive. 12

Von Mises stress (MPa) distribution in the polypropylene system for the pressure of 7, 15 and 30 MPa.
Total displacement distribution
We can also observe the displacements in the three-layer polypropylene system. Figure 5 shows the maps of the deformities during the application of the different pressure levels on the polypropylene thermal protection system. We can then note for the pressure of 7 MPa a maximum displacement of 0.054 mm, then for the pressure of 15 MPa the maximum displacement is 0.1158 mm and finally for the pressure of 30 MPa this maximum displacement is 0.2521 mm. The calculations show that the viscous response of the material consists largely of a viscoelastic strain. 29 The results allow to note that the time-dependent response of the material can be sufficiently described under a low pressure level by a linear viscoelastic behavior model.3,12

Distributions of the total displacement (mm) in the polypropylene system for the pressure of 7, 15 and 30 MPa.
Comparison under hydrostatic pressure
We now compare the resistance of syntactic foam coating systems to those of polypropylene. A study had been carried out by Grosjean et al.
6
on the influence of the hydrostatic pressure applied to the pipe structure coated with syntactic foam. The experiment details of the confined Compression test is found in Grosjean et al.
6
A confined compression set-up allowed the application of a non-isotropic solicitation with a strain rate of
Figure 6 shows the different distributions of the main stresses in the thickness of the Glass Syntactic Polypropylene material. The results show higher inward Von Mises stresses. The Von Mises criterion is moreover based on the presence of deviatoric stresses. This will justify the damage of the Glass syntactic polypropylene material which seems to increase with the presence of these deviatoric stresses directed toward the assembly with the steel material. 1 The results of the numerical calculations carried out in this work show deviatoric stresses in the thickness of the syntactic foam admitting a non-isotropic pressure loading.1,6 While these stress levels are lower in the thicknesses of polymer materials (see Figure 7). This difference in distribution is also due to the particular properties of the composite with respect to those of the polymers.1,2,13,14 In literature, Grber et al. 9 studied the buckling performance of thick composite cylinders loaded in compression. The results showed that the buckling failure would occurs before the strength failure under the hydrostatic pressure. The calculations carried out within the framework of this study (see Figures 6 and 7) show that the behavior of polymer structures can be modeled with the finite element calculation code and can also be extended to cyclic loading to include the fatigue of the constituent materials. 29 From the numerical modeling, we can see that damage could occur in the thickness of the syntactic foam under hydrostatic pressure as predicted by the literature.6,9

Stress repartitions in the GSPP thickness (GPa).

Stress repartitions in the polypropylene thickness.
Time-dependent mechanical responses
Viscoelasticity of the polypropylene material
The analysis of the response of the polymer and the composite materials is carried out here under the uniaxial tensile stress. We consider the creep phenomena which will translate the application of a load in the thickness of materials during a rather significant period. 3
Experimental studies of the assembly behavior were carried out by Cognard et al. 31 The analysis therefore made it possible to improve the results with the Arcan-modified assembly.32,33 Recently, a model was developed for an epoxy adhesive layer by the same experimental technique. 17 This model is based on the Gaussian form spectral approach. 17 Our objective is to develop a database for the analysis of the long-term behavior of a multilayer system using numerical analysis. This approach seems advantageous to us because the authors implemented the model in the finite element calculation code of the Comsol software. A few years ago, the use of the polypropylene material for the pipeline coating was recognized as a good candidate for a three-layer assembly due to its toughness in the marine environment.26,27 This type of material still lends itself well to use as a coating in the multilayer system. 12 Under creep-type loading, the viscoelastic behavior of a polypropylene material can be modeled using two branches of the generalized Maxwell model presented in section (3.2). 3 The installation phase always imposes axial stresses in the thicknesses of the pipeline structural materials. Until now, the numerical calculations of the stresses and strains generated by this type of loading have only concerned the steel tube.4,5 Significant shifts in the thicknesses of the materials studied have been observed in the previous sections. Furthermore, a study in the quasi-static field carried out by Vestrum et al.13,14 made it possible to note the contribution of the polymer coating in the overall resistance of the multilayer structure. A numerical simulation of this installation phase for very short-term creep phenomena under the different types of loading is necessary to predict the level of stress developed. This would in fact make it possible to predict damage and obtain the optimal performance of multilayer structures. 3 Thus, to carry out the numerical calculations for such phenomena, we considered the characteristics of the materials of polypropylene and of polypropylene matrix syntactic foam already provided by Phan et al. 3
The numerical simulation of the response of the assembly is carried out in time regime and for an elastic model chosen for the Glass Syntactic Polypopylene by considering the temperature level of 20°C and a force of 300 kN applied to the external borders. We can represent in 2D axisymmetric (see Figure 2(c)), the geometry of the assembly of the Glass Syntactic Polypropylene material from the inside to the outside with the solid polypropylene.2,12,13 The viscoelastic parameter of the polypropylene materials are given in Table 3. It is thus possible to observe in Figure 8 the stress distributions of Von Mises and those of the total displacement in the thicknesses of the materials of the thermal insulation system for the applied force. The numerical calculations were carried out in the time-dependent framework and we can observe on the mapping the zones for levels of the maximum Von Mises stress in Figure 8(a) and of the maximum displacement in Figure 8(b). These areas of high stress concentration are found at the interfaces of polypropylene PP and Glass syntactic polypropylene materials. Once more the numerical calculations make it possible to observe the zones of high stresses in the interface part. 30 Interface areas are critical areas for the type of assembly under the axial tensile load applied in the thicknesses of the materials during the installation process. The modeling of a bonded assembly of polypropylene and syntactic foam was carried out in the study by Phan et al. 3 to analyze the process of making a field-joint which is the zone made up of two sections of pipe. The results obtained during this study made it possible to observe the fracture envelope of this assembly of the two materials. The levels of displacements obtained by our numerical calculations are comparable to the levels of displacements admitted during the experimental characterization carried out by Phan et al. 3
Viscoelastic parameters of the polypropylene material at different temperature. 11

(a) Von Mises stress repartition (MPa) and (b) distribution of the total displacement (mm).
Viscoelastic-viscoplastic behaviors of the GSPP
The mechanical behavior of polypropylene matrix syntactic foam has been studied in numerous works.1,3,10 The linear viscoelasticity hypothesis of syntactic foam has been best presented in the work of Phan et al.3,11 Indeed, the material describes a complex behavior including the constituents (matrix and hollow glass microspheres). We use the approaches of this viscoelasticity to resume the reversible part of the viscous deformation of the syntactic foam material. Table 4 provides us with the viscoelastic properties obtained during the characterization tests in uniaxial creep tension at different temperatures (T denotes temperature). In Figure 9, we can observe the results of the correlation of the simulation with the experiment of Phan et al. 3 The response of the material in creep under uniaxial tension include other time-dependent phenomena related to the nature of the matrix. In order to validate the numerical modeling, time-dependent tests have been conducted on a block of syntactic foam. Usually, A 100 KN capacity load cell with a temperature control system depending on the loadings conditions allows to perform the Creep-recovery tests. 18 The mechanical behavior of this material in creep consists of elastic, viscoelastic deformations. 3 Hook’s law and the spectral models make it possible to predict the volumetric deformation of the material under the loading in hydrostatic compression. 3 Within the framework of a cyclic loading, the model of Zapas-Crissmann will make it possible to describe the irreversibility of the strain in creep of the material under uniaxial tension. First, we consider that when the material is subjected to creep at a stress of 0.5 MPa, its total strain is then mostly reversible. We assume that the stress levels do not generate any damage to the material. Temperature dramatically alters the properties of the syntactic foam material. 3 In order to take this effect into account, the temperature levels chosen are 25° and 80° for the situations of pipeline storage and transport of the fluid in the marine environment. Figure 9 shows the material response obtained at 25°C. During the cyclic loading, a stress of 0.5 MPa is applied on the polypropylene syntactic foam. 3
Viscoelastic parameters of the GSPP material at different temperature. 3

Correlation between numerical simulation and experiment 3 at 25°C.
The evolution of the applied load in the band of values generally admissible in the real context had not been addressed. The numerical calculations using the viscoelastic parameters of the material for another load allowed to conclude on the observation that the temperature could modify the mechanical behavior of the composite materials for each load as predicted by Hamel et al. 34
Figure 10(a) and (b) presents the evolutions of the reversible creep strains for the two temperatures considered. We thus show an influence of temperature on the creep strain of the composite material. This influence is noted by observing the significant difference between the two strains of the material at stresses of 0.5 MPa and 1 MPa. The results show an increase in the reversible strain which occurs with temperature for the same level of creep stress. The response of the material is sufficiently described here with the spectral model proposed by Phan et al. 3

Axial strain with temperature and time: (a) for
We are now interested in the contribution of the irreversible deformation of the material during the loading and unloading cycles. We consider the stress levels low. Through the experimental analysis and by noting the level of contribution of the residual part in the recovery phase, we can observe a part dependent on the stress and the time for the temperature of 25°C (see Figure 11). 3 Figure 11(a) presents the evolution of the residual part with time through the Zapas-Crissmann model then its contribution to the axial creep strain. 21 The stress level chosen here is 0.5MPa. The number of load cycles used in the creep-recovery tests is five. The applied load induces cracks in the matrix of the composite material without causing microspheres failure. This allows to neglect the effect of the latter in the residual deformation. Thus all the residual deformation during tensile creep becomes mainly viscoplastic deformation due to the irreversibility of the matrix. We observe that despite this evolution over time, the contribution of this residual part compared to the instantaneous part is negligible. 29 We show for this level of applied stress that the viscous part can be assimilated to the linear viscoelasticity of the matrix. 3 Three longitudinal stress levels of 1.5, 2.5 and 3.5 MPa are chosen for the numerical simulation using the Table 4. We use Table 5 for a comparison with the additional residual strain. These longitudinal stress levels were chosen to determine the effect of residual strains in the viscous strain of the material. We can analyze and better define the residual part of the axial deformations (see Figure 11(b)). The contribution of the breakage of glass microspheres has been mentioned in the case of hydrostatic pressure loadings of 30 MPa in Phan et al. 1 This hydrostatic pressure level is however located below the collapse or crush pressure level of the material obtained experimentally at nearly 60 MPa. 3

(a) Residual strain evolution with time at 25°C. (b) Effect of residual strain on mechanical response by varying stress level at 25°C.
Viscoplastic parameters for the GSPP material.
We then show the changes over time of the axial strain as a function of different longitudinal stresses applied in uniaxial traction in Figure 11(b). We can observe the differences in axial strain between the reversible and irreversible process which become much more significant with increasing stress. Under conditions of use in an offshore environment, the thermal insulation systems of the pipes will have an overall behavior depending on the weather, temperature and loading level. During their deformation mechanism, the response of the materials will involve their viscoelastic and viscoplastic character and mainly when the load is not initially damaging. Similar conclusions were retained in the studies carried out in the literature.16,18,19 The use of the power-law offer a potential tool for predicting the time and stress of long-term creep rupture of the syntactic foam. The conclusion is the same as in Hamel et al. 34 So, the failure strain at multiple stress levels is observed and the tensile failure of the Glass Syntactic polypropylene foam is due to the stress level.
Conclusion
In this paper, we analyzed the mechanical behavior of polypropylene and syntactic foam materials for thermal insulation lining of offshore pipelines. We studied the effects of mechanical loadings applied to three-layer and five-layer systems under suspension conditions in offshore environments. Results showed a contribution of polymer and composite coatings in the overall strength of the structure when axial tensile forces are applied. Using Maxwell’s model, the results showed quite significant admissible displacements in material assemblies. A characterization of the syntactic foam showed the viscoelastic and viscoplastic behaviors. The residual part can be neglected at low stresses. This will require to consider the behavior of the material as linear viscoelastic. While for a rather significant contribution of this residual deformation, the temperature variation may allow a significant increase in the deformation of the syntactic foam material. Under this effect of temperature, the damage of a polypropylene matrix composite material occurs before its strength failure during cyclic loadings, highlighting the criticality of polymer-composite materials assemblies used in deep sea where polymers are used as adhesives. In perspective, an analysis of the assembly behavior with syntactic foam and polymeric materials should be performed to allow a more robust modeling of the thermal insulation material response under the conditions of mechanical loading and temperature.
Supplemental Material
Supplemental Material, sj-pdf-1-jtc-10.1177_08927057211051417 - Time-dependent response of thermoplastic matrix composite material under the cyclic loadings
Supplemental Material, sj-pdf-1-jtc-10.1177_08927057211051417 for Time-dependent response of thermoplastic matrix composite material under the cyclic loadings by Ismael Figapka Pagore, Guy Richard Kol and Jean Gambo Betchewe in Journal of Thermoplastic Composite Materials
Footnotes
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
