Next Article in Journal
Improved Measurements of the Physical Properties of Oriental Lacquers Using Atomic Force Microscopy and a Nanoindenter
Previous Article in Journal
Fabrication of Polycaprolactone/Nano Hydroxyapatite (PCL/nHA) 3D Scaffold with Enhanced In Vitro Cell Response via Design for Additive Manufacturing (DfAM)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A General Temperature-Dependent Stress–Strain Constitutive Model for Polymer-Bonded Composite Materials

1
Graduate School of China Academy of Engineering Physics, 10 Xibeiwang E. Rd., Beijing 100193, China
2
Institute of Chemical Materials, China Academy of Engineering Physics, 64 Mianshan Rd., Mianyang 621900, China
3
School of Reliability and Systems Engineering, Beihang University, 37 Xueyuan Rd., Beijing 100191, China
*
Author to whom correspondence should be addressed.
Polymers 2021, 13(9), 1393; https://doi.org/10.3390/polym13091393
Submission received: 1 April 2021 / Revised: 17 April 2021 / Accepted: 21 April 2021 / Published: 25 April 2021
(This article belongs to the Section Polymer Processing and Engineering)

Abstract

:
This study develops a general temperature-dependent stress–strain constitutive model for polymer-bonded composite materials, allowing for the prediction of deformation behaviors under tension and compression in the testing temperature range. Laboratory testing of the material specimens in uniaxial tension and compression at multiple temperatures ranging from −40 C to 75 C is performed. The testing data reveal that the stress–strain response can be divided into two general regimes, namely, a short elastic part followed by the plastic part; therefore, the Ramberg–Osgood relationship is proposed to build the stress–strain constitutive model at a single temperature. By correlating the model parameters with the corresponding temperature using a response surface, a general temperature-dependent stress–strain constitutive model is established. The effectiveness and accuracy of the proposed model are validated using several independent sets of testing data and third-party data. The performance of the proposed model is compared with an existing reference model. The validation and comparison results show that the proposed model has a lower number of parameters and yields smaller relative errors. The proposed constitutive model is further implemented as a user material routine in a finite element package. A simple structural example using the developed user material is presented and its accuracy is verified.

1. Introduction

Polymer-bonded composites materials (PBMs) have drawn increasing attention in recent years due to their advantage in complex shaping, bi-injection molding, and ease of machining. By compounding different filler materials with polymer, new functional materials can be made, providing a viable technology for on-demand material design and engineering [1]. Applications of PBMs include, but are not limited to, polymer-bonded magnetics [2,3], energetics [4,5], aluminum [6], fiber [7,8], and many others [9,10,11]. To utilize those functional materials in critical applications, understanding their mechanical properties and deformation behaviors under different loading and temperature conditions is of great importance in meeting the basic design criteria.
The mechanical properties such as compliance, modulus, toughness, and strength of polymer materials depend largely on strain (rate) and temperature [12,13]. The stress–strain behavior of PBMs is usually that of viscoelastic plastic, which can be approximated as an elastic (linear) regime followed by a plastic (nonlinear) regime [14] in phenomenological view. The deformation is homogeneous in the elastic regime and can be related to the applying stress by Hooke’s law [15]. Once the stress of the material reaches a certain stress level, the material will produce irreversible strain and plastic deformation [16]. The maximum stress below which the resulting strain remains proportional to the applying stress is the elastic limit of the polymer materials. However, the filler materials can drastically alter the elastic limit and change the stress–strain behavior in both elastic and plastic regimes [17,18,19]. In addition, the plastic regime of the PBM can also be vastly different to that of the regular polymer materials due to the use of different filler materials.
The need to incorporate temperature effects into the deformation behavior of PBMs has been identified in several experimental studies on a variety of materials, for example, EDC37 in the temperature range of −65 C to 60 C by Williamson et al. [20], Rowanex between −60 C and 60 C by Walley et al. [21], HTPB at low temperatures by Chen et al. [22], short carbon-fiber-reinforced polyether-ether-ketone (SCFR/PEEK) composites between 20 C and 235 C by Zheng et al. [23], epoxy mortar at temperature between 80 C to 210 C by Vecchio et al. [24], nanoparticle/epoxy nanocomposites between 23 C to 53 C by Unger et al. [25], and many others [26,27]. Several models have been reported to describe the mechanical responses of PBM under different temperatures, including the elastic–plastic fracture-mechanism-based model using the glass-transition temperature ( T g ) as a piecewise knot point with a total number of 39 parameters [28], temperature-dependent mesoscale interface-based model [29], and the temperature-dependent visco-hyperelastic model [30]. Models based on the concept of temperature correction factor and its variants have been seen, for example, in [31,32]. Experimental evidence indicates that the mechanical behaviors of PBMs under tension and compression are quite different, and the stress–strain behavior under tension and compression may be described separately [20,21,28]. The aforementioned models in [29,30,31] were established for a single mode, and their applicability to the other mode is unknown. Other models such as the one reported in [28] and its variants requires a dozen parameters, demanding much more testing data for parameter calibration, which may not be realistic in all cases.
The goal of this study is to develop a general temperature-dependent stress–strain constitutive model for PBMs suitable for both tension and compression loading modes with a minimal set of parameters. Compression and tension testing in the temperature range from −40 C to 75 C is performed to obtain stress–strain response data. Based on the observation of the stress–strain curves at individual temperatures, a nonlinear model based on the Ramberg–Osgood (RO) relationship is proposed to describe the stress–strain under a given temperature. By correlating the model parameter with temperature, a general temperature-dependent constitutive model can be formulated. The reminder of the study is organized as follows. First, the experimental work is presented in detail and the data in compression and tension are acquired. Next, the constitutive model is developed based on the testing data. Following that, the model validation using both an independent set of testing data and third-party data is performed. The effectiveness of the model is demonstrated, and the performance of the model is compared with a reference model having a comparable number of parameters. The developed constitutive model is further implemented as a user material in a finite element package for structural applications. Finally, conclusions are drawn based on the current results.

2. Experimental Testing

Figure 1 presents the overall process of experimental testing, model development, validation, and comparison. In the experiment, test specimens were prepared according to the national standard [33] and were arbitrarily divided into two groups, e.g., 80% for model development and 20% for validation. The testing is performed at various temperatures in the range of −40 C to 75 C under tension and compression. In the model development part, a model based on the RO relationship is proposed to describe the stress–strain data under each of the temperatures. The temperature dependence on model parameters is established using the response surface model and the final temperature-dependent constitutive model is formulated. In model validation, the validation data (20%) and third party data are used to verify the effectiveness of the proposed model. In model comparison, the performance of the proposed model is compared with an existing reference model having a comparable number of parameters.

2.1. Specimens Preparation and Test Conditions

The composition of the PBM used in this study is barium sulfate of 94 wt.% embedded in 6 wt.% fluororubber binder. The conversion volume percentages are 87.6% and 12.4%, respectively. The size of the filler particles is in the range of 0.5 mm and 3 mm. The test specimens were obtained from a vendor and the actual image of the PBM sample is shown in Figure 2a. The PBM was inspected using cone-beam computed tomography (CT), as shown in Figure 2b, to ensure the integrity before testing. In the shown CT image, taken from a piece of the material sample, the light-colored particles are filler materials and the dark-colored portion are the binder materials. After CT inspection, the cylinder specimens for compression ( Φ 20 mm × 20 mm) and dog-bone specimens for tension ( Φ 15 mm × 50 mm) were fabricated according to the standard [33], with the dimensions shown in Figure 3a,b, respectively.
A uniform scheme of design of experiments was used to determine the design points of the temperature. In this study, 20 test conditions are designed in the temperature range of −40 C to 75 C and the temperature interval is less than 10 C, totaling a number of 191 test specimens. The temperatures and the number of specimens at the given temperature are summarized in Table 1. Moreover, one validation specimen is selected for every 4 points, excluding the boundary temperature points.

2.2. Test Results

Uniaxial tension and compression tests were performed on a universal testing machine in an environmental chamber at the specified temperature. The loading rate is 0.5 mm/min, representing a quasi-static condition. Figure 4 and Figure 5 present the average of the tensile and compressive stress–strain curves at different temperatures from initial state to rupture. It should be noted that the stress and strain under tension and compression are expressed by positive numbers in this paper.
It is observed from Figure 4 that the ultimate strength decreases significantly as the temperature increases, e.g., from 11.55 MPa at −40 C to 2.69 MPa at 75 C. The compression results shown in Figure 5 can be divided into two regimes. As the strain increases, the stress increases monotonically until the ultimate strength is reached, following by the postultimate stage before fracture. With the further increase of the strain, the stress decreases gradually. For compression results, the ultimate strength decreases from 42.41 MPa at −40 C to 12.64 MPa at 75 C. The maximum resistance of a material to external forces is called the stress at break [34]. The ultimate strength in tension and compression indicates the failure. Therefore, the data ranging from the initial state to the stress at break are used for the model development. The main transition temperature for the used binder material is T g = 80 C, which is consistent with that of a similar material [35]. The stress–strain behavior at a higher temperature that is close to T g may be influenced by this transition temperature. This can be observed from Figure 4 and Figure 5, where at higher temperatures, the stress–strain curves start to deviate much more than those at lower temperatures. As the proposed constitutive model is empirical in nature, the developed model does not deal with the mechanism related to the transition temperature.

3. Constitutive Model Development

The stress–strain model under single temperature is proposed first, and the temperature dependence of the model parameters is established using the tension and compression data. After that, the general temperature-dependent stress–strain constitutive model is obtained.

3.1. Stress–Strain Model under Single Temperature

The stress–strain curves presented in Figure 4 and Figure 5 indicate that the total response consists of a short linear part and a nonlinear part. Therefore, the total strain ε can be written as the linear accumulation of the two parts such as
ε = ε e + ε p ,
where ε e is linear (elastic) strain and ε p is nonlinear (plastic) strain. For the linear strain, it can be described by the stress σ and Young’s modulus E as ( σ / E ) . For the nonlinear part, a general exponential form as in the Ramberg–Osgood relationship [36] can be used. Substituting the linear and nonlinear parts, for monotonic loading,
ε = ε e + ε p = σ E + K ( σ E ) n ,
to determine the adequacy of Equation (2) in describing the testing data at a single temperature. The testing data of T = −40 C are arbitrarily chosen to fit the model parameters, and the fitted model parameters are presented in Table 2 for both tension and compression cases. Using the fitting parameters, the results of the mean fit and the testing data are presented in Figure 6 for comparison purposes. It can be observed that the model can capture the linear and nonlinear portions of the testing data fairly well in the given strain range for both tension and compression cases, indicating the effectiveness of Equation (2) for correlating the stress–strain behavior of the PBM under investigation.

3.2. Tension

3.2.1. Model Parameters under Different Temperatures

The nonlinear least squares estimator is used to obtain the parameters ( E , K , n ) in Equation (2) for tension testing data of each of the specimens at each of the temperatures. For convenience, the parameter ln K instead of K is used. The resulting parameters vs. temperatures are presented in the Figure 7. Every solid dot in Figure 7 is the fitted parameter associated with one specimen. From the results of Figure 7, it can be found that the variation of each of the three parameters with temperatures can loosely be divided into two linear regimes. For parameters ln K and n, the change points of the two regimes are about −20 C, indicated by the vertical dash lines in Figure 7a,b. The two parameters decrease linearly from −40 C to the change point, and increase linearly from that point to the end. For parameter E, the change point is also −20 C but the variation in the two regimes is reversed compared to that of the other two parameters. The scattering in the fitted parameters shown in Figure 7 is caused by the inherent uncertainty in material properties and stochastic nature of the filler particle distribution.

3.2.2. Temperature Dependence Modeling

Based on the above discussion, it can be observed that the temperature dependence can be described using the piecewise linear relationship; therefore, a piecewise linear model is used to correlate the variation of the three parameters with temperatures.
{   ln K = α 1 + α 2 · T + α 3 · ( T T k 1 ) · H ( T T k 1 ) n = β 1 + β 2 · T + β 3 · ( T T k 2 ) · H ( T T k 2 ) E = γ 1 + γ 2 · T + γ 3 · ( T T k 3 ) · H ( T T k 3 ) ,
where α i , β i , and γ i ( i = 1 , , 3 ) are fitting coefficients without unit; T is the temperature in degrees centigrade; and T k i ( i = 1 , , 3 ) are the knot points of the piecewise linear function. The term H ( x ) is the Heaviside step function defined as
H ( x ) = {   1 x 0 0 otherwise .
According to Figure 7, T k i ( i = 1 , , 3 ) are determined as −20 C. Using the parameter-fitting results presented in Figure 7 and Equation (3), the temperature dependence fitting coefficients α , β , and γ are obtained using the ordinary least square estimator as
Tension :     {   α = [ α 1 , α 2 , α 3 ] = [ 3.522 , 0.06351 , 0.3691 ] β = [ β 1 , β 2 , β 3 ] = [ 1.516 , 0.01018 , 0.05640 ] γ = [ γ 1 , γ 2 , γ 3 ] = [ 13 , 374 , 70.93 , 166.2 ] .
Figure 8 presents the temperature dependence fitting results, where the mean fit (solid lines) is shown. It can be seen that the mean fit of the piecewise linear model in Equation (3) can roughly describe the temperature dependence of the model parameters ln K , n, and E.
Substituting Equation (3) into Equation (2) with parameters in Equation (5), the final temperature-dependent stress–strain constitutive model under tension is obtained as
Tension : ε = σ E ( T ) + K ( T )   [ σ E ( T ) ] n ( T ) ,
where the temperature terms are
Tension : { ln K ( T ) = 3.522 0.06351 · T + 0.3691 · ( T + 20 ) · H ( T + 20 ) n ( T ) = 1.516 0.01018 · T + 0.05640 · ( T + 20 ) · H ( T + 20 ) E ( T ) = 13 , 374 + 70.93 · T 166.2 · ( T + 20 ) · H ( T + 20 ) .
The temperature-dependent stress–strain model (Equation (6)) is used to predict the strain given the temperature and stress. The overall residuals on all the tension data are evaluated and presented in Figure 9. The value corresponding to the sum of squared errors (SSE) is 3.0686 × 10 4 using Equation (8). The SSE is given as
SSE = i = 1 N ( y i y ^ i ) 2 ,
where y i is the actual value, y ^ i is the model prediction, and i = 1 , , N represents the index of a total number of N data points. The standard deviation of the residuals is estimated as 5.133 × 10 5 .

3.3. Compression

3.3.1. Model Parameters under Different Temperatures

The same fitting and processing method in Section 3.2.1 is applied to the compression data, and the resulting parameters in Equation (2) under different temperatures are presented in Figure 10.
Similar to the tension results, the variation of parameters with the temperatures under compression are also piecewise linear. Both ln K and n are linear proportional to the temperature in the range of −40–40 C. In the range of 40–75 C, the two parameters are inversely proportional to the temperature. The parameter E shows the same behavior but with a different change point of −5 C. Therefore, for compression data, the same equation, Equation (3), is used with T k 1 = T k 2 = 40 C and T k 3 = −5 C.

3.3.2. Temperature Dependence Modeling

Using the parameter-fitting results presented in Figure 10 and Equation (3), the temperature dependence fitting coefficients α , β , and γ are obtained using the ordinary least squares estimator as
Compression : {   α = [ α 1 , α 2 , α 3 ] = [ 13.74 , 0.1125 , 0.4976 ] β = [ β 1 , β 2 , β 3 ] = [ 6.797 , 0.04343 , 0.1937 ] γ = [ γ 1 , γ 2 , γ 3 ] = [ 7665 , 0.7296 , 69.69 ] .
Substituting Equation (3) into Equation (2) with parameters in Equation (9), the final temperature-dependent stress–strain constitutive model under compression is obtained as
Compression : ε = σ E ( T ) + K ( T )   [ σ E ( T ) ] n ( T ) ,
where the temperature terms are
Compression : { ln K ( T ) = 13.74 + 0.1125 · T 0.4976 · ( T 40 ) · H ( T 40 ) n ( T ) = 6.797 + 0.04343 · T 0.1937 · ( T 40 ) · H ( T 40 ) E ( T ) = 7665 + 0.7269 · T 69.69 · ( T + 5 ) · H ( T + 5 ) .
Figure 11 presents the temperature dependence fitting results, where the mean fits (solid lines) are shown. It can be observed that the temperature dependence can be described using the piecewise linear functions. Figure 12 presents the residuals of the model fitting for compression data. The SSE in this case is 0.07323 and the standard deviation of the residuals is 0.001089.

4. Model Validations and Comparisons

Independent sets of testing data are used to validate the performance of the proposed model. Third-party data published for a PBM are further used to validate the generality of the model. To compare the proposed model with existing models, a model with a comparable number of parameters is adopted and results from the two models are compared. In the validation and comparison, the following two general measures, namely, the mean absolute error (MAE) and the root mean square error (RMSE), are used.
MAE = i = 1 n | y i y ^ i | N ,
RMSE = SSE N ,
where y i , y ^ i , i, and N are defined as in Equation (8).

4.1. Validation Using Independent Sets of Testing Data

The data from the validation specimens (labeled as ‘V’ in Table 1) are not used for model development and are used for validation. It is noted that both tension validation specimens and compression validation specimens are available. We utilize data associated with the validation specimens. Prediction results of the developed model and the actual testing values are compared. In addition, 2 groups of modeling data near the piecewise point are added to compare with the model fitting results. Figure 13 presents the model prediction results for tension specimens compared with testing results. Figure 14 shows the model prediction results for compression specimens compared with testing results. In both modes, a close agreement between the prediction and testing data can be observed.
The maximum MAE and RMSE are 6.200 × 10 5 and 8.625 × 10 5 at 25 C, respectively, for tension data. For compression data, the maximum MAE and RMSE are 0.001741 and 0.002460 at 10 C, respectively. The prediction performance for each of the temperatures are detailed in Table 3.

4.2. Validation Using Third-Party Data

To further verify the generality of the proposed model, third-party testing data on PBX 9502 from [28] consisted of 5 sets of tension data and 5 sets of compression data (at temperatures −52 C, −20 C, 20 C, 50 C, and 72 C) are used to validate the generality of the proposed model.
It should be noted that the results of model parameters in Equations (7) and (11) cannot be directly used for prediction purposes as those parameters are material specific. The reported data on PBX 9502 are used to fit the developed model using Equation (2) for tension and compression data at each of the temperatures, and the temperature dependence for both modes are obtained using Equation (3). The same least squares fitting scheme is used. The prediction results of the proposed model and the actual testing data at the five temperatures are compared and shown in Figure 15. Under tension, the maximum MAE and RMSE are 5.4472 × 10 5 and 7.8292 × 10 5 at 50 C, respectively. For compression, the maximum MAE is 0.001040 at 74 C and the maximum RMSE is 0.001184 at 20 C. For both cases, the model can capture the behaviors of the stress–strain variations fairly well.

4.3. Model Comparisons

To further investigate the performance of the proposed model, a stress–strain constitutive model reported by Chang et al. [37] is adopted for comparison purposes. This model is formed by appending a quadratic term and a cubic term of the strain to the regular linear model, as in Equation (14):
σ = E 0 ε + C 1 ε 2 + C 2 ε 3 ,
where the parameters E 0 , C 1 , and C 2 are fitting parameters. It is worth mentioning that other reported models can also be adopted for comparisons; however, this model is chosen due to the fact that the model has a comparative number of parameters and it requires only stress–strain data for fitting. Other models reviewed in the introduction either contain too many parameters or require special testing data.

4.3.1. Tension

The same amount of tension data is used to fit model parameters using Equation (14) for each of the temperatures. The fitted model parameters ( E 0 , C 1 , C 2 ) are further correlated with temperatures to obtain the following temperature dependence equations in Equation (15):
{   E 0 = 9691 56.78 · T 0.9012 · T 2 + 0.01122 · T 3 C 1 = 2.311 × 10 6 + 74 , 213 · T 878.7 · T 2 32.35 · T 3 + 0.4033 · T 4 C 2 = 1.061 × 10 8 3.471 × 10 7 · T + 5.458 × 10 5 · T 2 + 16 , 787 · T 3 222.5 · T 4 .
The two models with established model parameters are used to predict the validation data (those labeled as ‘V’ in Table 1), and the results are shown in Figure 16. In general, both models can describe the data satisfactorily; however, for data at 50 C, the proposed model outperforms the reference model. The errors in terms of MAE and RMSE of the two models are presented in Table 4. From Table 4, it can be seen that the proposed model yields smaller errors for all data sets except the one associated with 0 C.

4.3.2. Compression

The same amount of data under compression are used to fit the reference model, resulting in the following temperature dependence functions in Equation (4).
{   E 0 = 8335 80.90 · T 1.069 · T 2 + 0.02098 · T 3 C 1 = 7.200 × 10 5 + 9184 · T + 106.7 · T 2 3.504 · T 3 + 0.01548 · T 4 1.420 × 10 4 · T 5 C 2 = 2.201 × 10 9 3.700 × 10 7 · T 5363 · T 2 + 172.5 · T 3 .
The comparisons for compression data between the two models are made following the same procedure described above and are shown in Figure 17, with the error results compared in Table 5. In general, the proposed model can yield improved prediction results, particularly for the tail regions of the stress–strain data. The detailed error results in Table 5 indicate that the proposed model outperforms the reference model at all temperatures except 10 C. It should be noted that other probabilistic methods such as Bayesian model assessment [38,39,40] can also be employed to compare the performance of the models.

5. A User Material Implementation of the Model

To use the proposed temperature-dependent stress–strain model for structural applications, a user material subroutine for the proposed model is developed for finite element analysis. The dimensional of the structural model is the same as the actual specimens. The structural models for the tension and compression specimens are created and meshed using quadratic hexahedron elements with an average size of 2 mm.
For the tension specimen, a constant stress of 2 MPa is applied to the top face with the bottom face fixed in all degrees of freedom, representing a tension test condition. For the compression specimen, a constant stress of −20 MPa is applied on the top face, and the same boundary conditions are used for the bottom face. Figure 18 demonstrates the resulting normal stresses after applying the load at −40 C, which are shown for the tension and compression specimens in Figure 18a,b, respectively. Under tension, it can be found that the stress on the top face is about 2 MPa, which is the same as the applied load, and the maximum stress occurs in the middle section of the structural part. Under compression, the maximum stress appears at the fixed position. From the top face, the stress gradually decreases along the UX direction. The stress on the top face is the same as the applied load of 20 MPa. For both cases, the resulting stresses are consistent with the actually applied loads.
Figure 19 presents the FEA results using the specimen structural models in Figure 18 with the developed user material routine for temperatures of −40 C, 20 C, and 75 C and compares them with the theoretical results of the proposed temperature-dependent stress–strain model. Both tension and compression results are compared. The maximum relative error (RE) is 0.01239, occurring at −40 C for the tension case, and the maximum RE is 0.01537 at 75 C for the compression case, indicating the effectiveness of the developed user material routine for structural applications.

6. Conclusions

A general temperature-dependent stress–strain constitutive model for polymer-bonded composite materials is proposed, developed, validated, and compared. Experimental testing data were acquired in the temperature range of −40–75 C for both tension and compression. The model under a given temperature is proposed based the linear accumulation of an elastic term and a plastic term, resembling the Ramberg–Osgood relationship. The model parameters are obtained by fitting the testing data at each of the temperatures, and the temperature dependence of the model parameters can be established using the response surface method. The proposed model is validated using independent sets of testing data and third-party data, and is further compared with an existing model with a comparable number of parameters. A user material routine of the proposed model is implemented and verified for structural applications. Based on the current results, the following conclusions are drawn.
  • The proposed method can unify the temperature effect by modeling the temperature dependence of the constitutive model parameters, thus providing an alternative to the existing temperature correction factor-based methods. Furthermore, such a treatment allows a minimal number of fitting parameters compared to the existing models.
  • The proposed method provides a general approach to model the stress–strain behaviors of PBMs. The basic assumption that the stress–strain response is an accumulation of an elastic term and a power-law plastic term is observed in many other materials such as various metals. The effectiveness and generality of the model are validated and demonstrated by two independent preponderate data sets in the temperature range −40–75 C.
It is worth mentioning that the developed model is a semiempirical data-driven model, and it does not deal with the microscopic mechanism. The dynamic mechanical analysis on the material can be carried out to enhance the understanding of the mechanism, and SEM imaging may be used to reveal the interface interaction between the filler particles and the binder materials. The constitutive model established in this study does not incorporate the strain rate effect. The strain-rate-induced hardening phenomenon has been reported in several studies, e.g., [41,42]. The incorporation of the strain rate effect to the proposed model will be studied in future work. The proposed model considers the PBM as a homogeneous material in macroscopic perspective. The local variations of the material properties due to the microscopic effects of filler materials and polymer binder interface [43,44,45] may be explored in the future. In addition, the effectiveness of the model is verified using the presented experimental data acquired in this study and a third-party published data. Due to the limit on currently accessible data, the verification of the developed model using other PBMs will be investigated in the future. The current modeling adopts a deterministic approach, and a probabilistic approach can also be employed for reliability analysis [46,47,48].

Author Contributions

Conceptualization, X.D. and X.G.; methodology, X.D., X.G., and J.H.; software, X.D.; validation, X.D., X.G., and J.H.; formal analysis, X.D.; investigation, X.D.; resources, H.Y. and W.T.; data curation, H.Y. and W.T.; writing—original draft preparation, X.D.; writing—review and editing, X.G.; visualization, X.D. and X.G.; project administration, X.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by The National Natural Science Foundation of China, Nos. 51975546, U1930403, and the research fund from the Institute of Chemical Materials, CAEP. The support is gratefully acknowledged.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

This study did not involve humans.

Data Availability Statement

All data are included within the text.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
PBMsPolymer-bonded composites materials
RORamberg–Osgood
CTComputed tomography
SSESum of squared errors
MAEMean absolute error
RMSERoot mean square error
FEAFinite element analysis

References

  1. Vavlekas, D.; Ansari, M.; Hao, H.; Fremy, F.; McCoy, J.L.; Hatzikiriakos, S.G. Zero Poisson’s ratio PTFE in uniaxial extension. Polym. Test. 2016, 55, 143–151. [Google Scholar] [CrossRef]
  2. Dutoit, B.; Besse, P.A.; Blanchard, H.; Guérin, L.; Popovic, R. High performance micromachined Sm2Co17 polymer bonded magnets. Sens. Actuators A Phys. 1999, 77, 178–182. [Google Scholar] [CrossRef]
  3. Wang, B.; Kari, L. Constitutive model of isotropic magneto-sensitive rubber with amplitude, frequency, magnetic and temperature dependence under a continuum mechanics basis. Polymers 2021, 13, 472. [Google Scholar] [CrossRef]
  4. Palmer, S.; Field, J.E.; Huntley, J. Deformation, strengths and strains to failure of polymer bonded explosives. Proc. R. Soc. Lond. Ser. A Math. Phys. Sci. 1993, 440, 399–419. [Google Scholar]
  5. Siviour, C.; Gifford, M.; Walley, S.; Proud, W.; Field, J. Particle size effects on the mechanical properties of a polymer bonded explosive. J. Mater. Sci. 2004, 39, 1255–1258. [Google Scholar] [CrossRef]
  6. Sang, J.; Hirahara, H.; Aisawa, S.; Kang, Z.; Mori, K. Direct Bonding of Polypropylene to Aluminum Using Molecular Connection and the Interface Nanoscale Properties. ACS Appl. Polym. Mater. 2019, 1, 2450–2459. [Google Scholar] [CrossRef]
  7. Zhang, J.; Duan, Y.; Wang, B.; Zhang, X. Interfacial enhancement for carbon fibre reinforced electron beam cured polymer composite by microwave irradiation. Polymer 2020, 192, 122327. [Google Scholar] [CrossRef]
  8. Li, Y.F.; Sio, W.K.; Yang, T.H.; Tsai, Y.K. A Constitutive Model of High-Early-Strength Cement with Perlite Powder as a Thermal-Insulating Material Confined by Caron Fiber Reinforced Plastics at Elevated Temperatures. Polymers 2020, 12, 2369. [Google Scholar] [CrossRef]
  9. Skokov, K.; Karpenkov, D.Y.; Kuz’Min, M.; Radulov, I.; Gottschall, T.; Kaeswurm, B.; Fries, M.; Gutfleisch, O. Heat exchangers made of polymer-bonded La (Fe, Si) 13. J. Appl. Phys. 2014, 115, 17A941. [Google Scholar] [CrossRef]
  10. Berezvai, S.; Kossa, A. Performance of a parallel viscoelastic-viscoplastic model for a microcellular thermoplastic foam on wide temperature range. Polym. Test. 2020, 84, 106395. [Google Scholar] [CrossRef]
  11. Yang, T.C. Effect of extrusion temperature on the physico-mechanical properties of unidirectional wood fiber-reinforced polylactic acid composite (WFRPC) components using fused deposition modeling. Polymers 2018, 10, 976. [Google Scholar] [CrossRef] [Green Version]
  12. Guo, H.; Guo, W.; Amirkhizi, A.V.; Zou, R.; Yuan, K. Experimental investigation and modeling of mechanical behaviors of polyurea over wide ranges of strain rates and temperatures. Polym. Test. 2016, 53, 234–244. [Google Scholar] [CrossRef] [Green Version]
  13. Abdel-Wahab, A.A.; Ataya, S.; Silberschmidt, V.V. Temperature-dependent mechanical behaviour of PMMA: Experimental analysis and modelling. Polym. Test. 2017, 58, 86–95. [Google Scholar] [CrossRef] [Green Version]
  14. Senden, D.J.; Krop, S.; van Dommelen, J.; Govaert, L. Rate-and temperature-dependent strain hardening of polycarbonate. J. Polym. Sci. Part B Polym. Phys. 2012, 50, 1680–1693. [Google Scholar] [CrossRef]
  15. Sakai, T.; Kurakazu, M.; Akagi, Y.; Shibayama, M.; Chung, U.I. Effect of swelling and deswelling on the elasticity of polymer networks in the dilute to semi-dilute region. Soft Matter 2012, 8, 2730–2736. [Google Scholar] [CrossRef]
  16. Van Loock, F.; Fleck, N.A. Deformation and failure maps for PMMA in uniaxial tension. Polymer 2018, 148, 259–268. [Google Scholar] [CrossRef]
  17. Li, Y.; Liu, Z. A novel constitutive model of shape memory polymers combining phase transition and viscoelasticity. Polymer 2018, 143, 298–308. [Google Scholar] [CrossRef]
  18. Johnsen, J.; Clausen, A.H.; Grytten, F.; Benallal, A.; Hopperstad, O.S. A thermo-elasto-viscoplastic constitutive model for polymers. J. Mech. Phys. Solids 2019, 124, 681–701. [Google Scholar] [CrossRef]
  19. Somarathna, H.; Raman, S.N.; Mohotti, D.; Mutalib, A.A.; Badri, K.H. Hyper-viscoelastic constitutive models for predicting the material behavior of polyurethane under varying strain rates and uniaxial tensile loading. Constr. Build. Mater. 2020, 236, 117417. [Google Scholar] [CrossRef]
  20. Williamson, D.; Siviour, C.; Proud, W.; Palmer, S.; Govier, R.; Ellis, K.; Blackwell, P.; Leppard, C. Temperature–time response of a polymer bonded explosive in compression (EDC37). J. Phys. D Appl. Phys. 2008, 41, 085404. [Google Scholar] [CrossRef]
  21. Walley, S.; Taylor, N.; Williamson, D. Temperature and strain rate effects on the mechanical properties of a polymer-bonded explosive. Eur. Phys. J. Spec. Top. 2018, 227, 127–141. [Google Scholar] [CrossRef]
  22. Chen, X.D.; Chang, X.L.; Zhang, Y.H.; Wang, B.; Zhang, Q.; Zhang, X. Tensile Mechanical Properties of HTPB Propellant at Low Temperature. In Key Engineering Materials; Trans Tech Publications Ltd.: Bäch, Switzerland, 2018; Volume 765, pp. 54–59. [Google Scholar]
  23. Zheng, B.; Deng, T.; Li, M.; Huang, Z.; Zhou, H.; Li, D. Flexural behavior and fracture mechanisms of short carbon fiber reinforced polyether-ether-ketone composites at various ambient temperatures. Polymers 2019, 11, 18. [Google Scholar] [CrossRef] [Green Version]
  24. Vipulanandan, C.; Mebarkia, S. Effect of strain rate and temperature on the performance of epoxy mortar. Polym. Eng. Sci. 1990, 30, 734–740. [Google Scholar] [CrossRef]
  25. Unger, R.; Arash, B.; Exner, W.; Rolfes, R. Effect of temperature on the viscoelastic damage behaviour of nanoparticle/epoxy nanocomposites: Constitutive modelling and experimental validation. Polymer 2020, 191, 122265. [Google Scholar] [CrossRef]
  26. Wang, P.; Yu, H.; Ma, R.; Wang, Y.; Liu, C.; Shen, C. Temperature-dependent orientation of poly (ether ether ketone) under uniaxial tensile and its correlation with mechanical properties. J. Therm. Anal. Calorim. 2020, 141, 1361–1369. [Google Scholar] [CrossRef]
  27. Pan, Z.; Liu, Z. A novel fractional viscoelastic constitutive model for shape memory polymers. J. Polym. Sci. Part B Polym. Phys. 2018, 56, 1125–1134. [Google Scholar] [CrossRef]
  28. Zubelewicz, A.; Thompson, D.G.; Ostoja-Starzewski, M.; Ionita, A.; Shunk, D.; Lewis, M.W.; Lawson, J.C.; Kale, S.; Koric, S. Fracture model for cemented aggregates. AIP Adv. 2013, 3, 012119. [Google Scholar] [CrossRef] [Green Version]
  29. Ge, S.; Zhang, W.; Sang, J.; Yuan, S.; Lo, G.V.; Dou, Y. Mesoscale Simulation to Study Constitutive Properties of TATB/F2314 PBX. Materials 2019, 12, 3767. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Li-Mayer, J.; Williamson, D.; Lewis, D.; Connors, S.; Iqbal, M.; Charalambides, M. Characterization and Micromechanical Modelling of a Temperature Dependent Hyper-Viscoelastic Polymer Bonded Explosive. Cranfield Online Res. Data (CORD) 2017. [Google Scholar] [CrossRef]
  31. Pulla, S.; Karaca, H.; Lu, Y. Numerical design of shape memory polymer composites with temperature-responsive SMA fillers. Compos. Part B Eng. 2016, 96, 287–294. [Google Scholar] [CrossRef] [Green Version]
  32. Tschoegl, N.W.; Knauss, W.G.; Emri, I. The effect of temperature and pressure on the mechanical properties of thermo-and/or piezorheologically simple polymeric materials in thermodynamic equilibrium—A critical review. Mech. Time-Depend. Mater. 2002, 6, 53–99. [Google Scholar] [CrossRef]
  33. GJB 772A-1997. Explosive Test Method; China Ordnance Industrial Standardization Research Institute: Beijing, China, 1997. [Google Scholar]
  34. ISO 527-1. Plastics-Determination of Tensile Properties. Part 1: General Principles; International Organization for Standardization: Geneva, Switzerland, 1993. [Google Scholar]
  35. Khanna, Y.P.; Kumar, R. The glass transition and molecular motions of poly (chlorotrifluoroethylene). Polymer 1991, 32, 2010–2013. [Google Scholar] [CrossRef]
  36. Ramberg, W.; Osgood, W.R. Description of Stress-Strain Curves by Three Parameters; NASA Technical Note NO.902; NASA: Washington, DC, USA, 1943. [Google Scholar]
  37. Chang, B.; Wang, X.; Long, Z.; Li, Z.; Gu, J.; Ruan, S.; Shen, C. Constitutive modeling for the accurate characterization of the tension behavior of PEEK under small strain. Polym. Test. 2018, 69, 514–521. [Google Scholar] [CrossRef]
  38. Guan, X.; Jha, R.; Liu, Y. Model selection, updating, and averaging for probabilistic fatigue damage prognosis. Struct. Saf. 2011, 33, 242–249. [Google Scholar] [CrossRef]
  39. Huo, H.; He, J.; Guan, X. A Bayesian fusion method for composite damage identification using Lamb wave. Struct. Health Monit. 2020. [Google Scholar] [CrossRef]
  40. He, J.; Chen, J.; Guan, X. Lifetime distribution selection for complete and censored multi-level testing data and its influence on probability of failure estimates. Struct. Multidiscip. Optim. 2020, 62, 1–17. [Google Scholar] [CrossRef]
  41. Okereke, M.I.; Akpoyomare, A.I. Two-process constitutive model for semicrystalline polymers across a wide range of strain rates. Polymer 2019, 183, 121818. [Google Scholar] [CrossRef]
  42. Ou, Y.; Zhu, D.; Zhang, H.; Huang, L.; Yao, Y.; Li, G.; Mobasher, B. Mechanical characterization of the tensile properties of glass fiber and its reinforced polymer (GFRP) composite under varying strain rates and temperatures. Polymers 2016, 8, 196. [Google Scholar] [CrossRef] [Green Version]
  43. Yang, L.; Yan, Y.; Liu, Y.; Ran, Z. Microscopic failure mechanisms of fiber-reinforced polymer composites under transverse tension and compression. Compos. Sci. Technol. 2012, 72, 1818–1825. [Google Scholar] [CrossRef]
  44. Su, C.; Wei, Y.; Anand, L. An elastic–plastic interface constitutive model: Application to adhesive joints. Int. J. Plast. 2004, 20, 2063–2081. [Google Scholar] [CrossRef]
  45. Koloor, S.S.R.; Rahimian-Koloor, S.; Karimzadeh, A.; Hamdi, M.; Petr, M.; Tamin, M. Nano-level damage characterization of graphene/polymer cohesive interface under tensile separation. Polymers 2019, 11, 1435. [Google Scholar] [CrossRef] [Green Version]
  46. He, J.; Huo, H.; Guan, X.; Yang, J. A Lamb wave quantification model for inclined cracks with experimental validation. Chin. J. Aeronaut. 2021, 34, 601–611. [Google Scholar] [CrossRef]
  47. Du, Y.M.; Ma, Y.H.; Wei, Y.F.; Guan, X.; Sun, C. Maximum entropy approach to reliability. Phys. Rev. E 2020, 101, 012106. [Google Scholar] [CrossRef] [PubMed]
  48. Du, Y.M.; Chen, J.F.; Guan, X.; Sun, C. Maximum Entropy Approach to Reliability of Multi-Component Systems with Non-Repairable or Repairable Components. Entropy 2021, 23, 348. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Flow diagram of the study process. There are 5 main parts: experiment, model development, validation, comparison, and FEM usermat.
Figure 1. Flow diagram of the study process. There are 5 main parts: experiment, model development, validation, comparison, and FEM usermat.
Polymers 13 01393 g001
Figure 2. (a) A sample of the PBM investigated in this study, and (b) schematic diagram of the testing process.
Figure 2. (a) A sample of the PBM investigated in this study, and (b) schematic diagram of the testing process.
Polymers 13 01393 g002
Figure 3. Geometry and dimension of the specimen. (a) The compression specimen, and (b) the tension specimen.
Figure 3. Geometry and dimension of the specimen. (a) The compression specimen, and (b) the tension specimen.
Polymers 13 01393 g003
Figure 4. Stress–strain curves under tension.
Figure 4. Stress–strain curves under tension.
Polymers 13 01393 g004
Figure 5. Stress–strain curves under compression.
Figure 5. Stress–strain curves under compression.
Polymers 13 01393 g005
Figure 6. Results of model mean fit using Equation (2) for testing data at T = −40 C. (a) Tension, and (b) compression.
Figure 6. Results of model mean fit using Equation (2) for testing data at T = −40 C. (a) Tension, and (b) compression.
Polymers 13 01393 g006
Figure 7. Temperature dependence of parameters under tension. (a) ln K , (b) n, and (c) E.
Figure 7. Temperature dependence of parameters under tension. (a) ln K , (b) n, and (c) E.
Polymers 13 01393 g007
Figure 8. Mean fit of the temperature dependence data under tension. (a) ln K , (b) n, and (c) E.
Figure 8. Mean fit of the temperature dependence data under tension. (a) ln K , (b) n, and (c) E.
Polymers 13 01393 g008
Figure 9. Histogram of residuals of the temperature-dependent stress–strain model (Equation (6)) on tension data.
Figure 9. Histogram of residuals of the temperature-dependent stress–strain model (Equation (6)) on tension data.
Polymers 13 01393 g009
Figure 10. Temperature dependence of parameters under compression. (a) ln K , (b) n, and (c) E.
Figure 10. Temperature dependence of parameters under compression. (a) ln K , (b) n, and (c) E.
Polymers 13 01393 g010
Figure 11. Mean fit of the temperature dependence data under compression. (a) ln K , (b) n, and (c) E.
Figure 11. Mean fit of the temperature dependence data under compression. (a) ln K , (b) n, and (c) E.
Polymers 13 01393 g011
Figure 12. Histogram of residuals of the temperature-dependent stress–strain model (Equation (10)) on compression data.
Figure 12. Histogram of residuals of the temperature-dependent stress–strain model (Equation (10)) on compression data.
Polymers 13 01393 g012
Figure 13. Prediction results of the test data under tension. (a,b) modeling data; (cf) validation data.
Figure 13. Prediction results of the test data under tension. (a,b) modeling data; (cf) validation data.
Polymers 13 01393 g013aPolymers 13 01393 g013b
Figure 14. Prediction results of the test data under compression. (a,b) modeling data; (cf) validation data.
Figure 14. Prediction results of the test data under compression. (a,b) modeling data; (cf) validation data.
Polymers 13 01393 g014aPolymers 13 01393 g014b
Figure 15. Validation results on third-party data in [28]. Solid lines are the model prediction results and discrete markers are testing data. (a) Tension, and (b) compression.
Figure 15. Validation results on third-party data in [28]. Solid lines are the model prediction results and discrete markers are testing data. (a) Tension, and (b) compression.
Polymers 13 01393 g015
Figure 16. Comparison of the viscoelastic model and the proposed model at validation temperatures under tension. The discrete markers represent testing data. (a) −35 C, (b) 0 C, (c) 25 C, and (d) 50 C.
Figure 16. Comparison of the viscoelastic model and the proposed model at validation temperatures under tension. The discrete markers represent testing data. (a) −35 C, (b) 0 C, (c) 25 C, and (d) 50 C.
Polymers 13 01393 g016
Figure 17. Comparison of the viscoelastic model and the proposed model at validation temperatures under compression. The discrete markers represent testing data. (a) −20 C, (b) 10 C, (c) 30 C, and (d) 50 C.
Figure 17. Comparison of the viscoelastic model and the proposed model at validation temperatures under compression. The discrete markers represent testing data. (a) −20 C, (b) 10 C, (c) 30 C, and (d) 50 C.
Polymers 13 01393 g017
Figure 18. The stress distribution results calculated by the finite element model at −40 C. (a) Tension, and (b) compression.
Figure 18. The stress distribution results calculated by the finite element model at −40 C. (a) Tension, and (b) compression.
Polymers 13 01393 g018
Figure 19. Comparison of FEA results and the proposed model at several temperatures. Solid lines represent ground truth by the proposed model, and discrete markers are numerical results extracted from FEA (user material). (a) Tension, and (b) compression.
Figure 19. Comparison of FEA results and the proposed model at several temperatures. Solid lines represent ground truth by the proposed model, and discrete markers are numerical results extracted from FEA (user material). (a) Tension, and (b) compression.
Polymers 13 01393 g019
Table 1. Test temperature, the number of repeated specimens, and the usage of data. The symbols ‘M’ and ‘V’ denote modeling and validation, respectively. Striking lines denote no testing performed at that temperature.
Table 1. Test temperature, the number of repeated specimens, and the usage of data. The symbols ‘M’ and ‘V’ denote modeling and validation, respectively. Striking lines denote no testing performed at that temperature.
TensionCompression
T ( C)Num.UsageNum.Usage
−405M5M
−355V
−306M6M
−25 5M
−205M5V
−155M5M
−105M
−55M5M
05V
55M5M
10 5V
153M5M
204M4M
258V6M
30 5V
355M5M
405M5M
455M6M
504V6V
553M6M
603M1M
653M2M
755M5M
Table 2. Fitting parameters using Equation (2) and testing data of T = −40 C.
Table 2. Fitting parameters using Equation (2) and testing data of T = −40 C.
ParametersTensionCompression
E (MPa) 1.120 × 10 4 7103
K214.5 2.664 × 10 13
n1.7796.955
Table 3. Model prediction performance on test data.
Table 3. Model prediction performance on test data.
ModeUsageTemperature ( C)MAERMSE
TensionM−30 4.592 × 10 5 3.311 × 10 5
M−20 2.225 × 10 5 3.188 × 10 5
V−35 3.008 × 10 5 4.547 × 10 5
V0 2.139 × 10 5 3.328 × 10 5
V25 6.200 × 10 5 8.625 × 10 5
V50 4.191 × 10 5 5.904 × 10 5
CompressionM−25 5.585 × 10 4 8.580 × 10 4
M−15 3.483 × 10 4 4.195 × 10 4
V−20 3.926 × 10 4 6.233 × 10 4
V100.0017410.002460
V300.0011000.001778
V500.0011950.001688
Table 4. Comparisons of the performance between the proposed model and the reference model for tension data.
Table 4. Comparisons of the performance between the proposed model and the reference model for tension data.
ModelTemperature ( C)MAERMSE
Proposed Model−35 3.008 × 10 5 4.547 × 10 5
0 2.139 × 10 5 3.328 × 10 5
25 6.200 × 10 5 8.625 × 10 5
50 4.191 × 10 5 5.904 × 10 5
Reference Model [37]−35 3.267 × 10 5 5.594 × 10 5
0 1.681 × 10 5 2.381 × 10 5
25 6.232 × 10 5 8.786 × 10 5
50 6.754 × 10 5 1.203 × 10 4
Table 5. Comparisons of the performance between the proposed model and reference model for compression data.
Table 5. Comparisons of the performance between the proposed model and reference model for compression data.
ModelTemperature ( C)MAERMSE
Proposed Model−20 3.926 × 10 4 6.233 × 10 4
100.0017410.002460
300.0011000.001778
500.0011950.001688
Reference Model [37]−20 4.722 × 10 4 7.916 × 10 4
10 7.109 × 10 4 9.527 × 10 4
300.0012620.002824
500.0021110.003094
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Duan, X.; Yuan, H.; Tang, W.; He, J.; Guan, X. A General Temperature-Dependent Stress–Strain Constitutive Model for Polymer-Bonded Composite Materials. Polymers 2021, 13, 1393. https://doi.org/10.3390/polym13091393

AMA Style

Duan X, Yuan H, Tang W, He J, Guan X. A General Temperature-Dependent Stress–Strain Constitutive Model for Polymer-Bonded Composite Materials. Polymers. 2021; 13(9):1393. https://doi.org/10.3390/polym13091393

Chicago/Turabian Style

Duan, Xiaochang, Hongwei Yuan, Wei Tang, Jingjing He, and Xuefei Guan. 2021. "A General Temperature-Dependent Stress–Strain Constitutive Model for Polymer-Bonded Composite Materials" Polymers 13, no. 9: 1393. https://doi.org/10.3390/polym13091393

APA Style

Duan, X., Yuan, H., Tang, W., He, J., & Guan, X. (2021). A General Temperature-Dependent Stress–Strain Constitutive Model for Polymer-Bonded Composite Materials. Polymers, 13(9), 1393. https://doi.org/10.3390/polym13091393

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop