Abstract

This paper is aimed at studying the effective mechanical property of shape memory polymer composites (SMPC) reinforced with natural short fibers. To this end, a novel modeling scheme was presented. The SMPC was firstly equivalent to the composite laminates, and the natural short fibers are also subtly equivalent to the ellipsoidal inclusions distributed in the matrix materials periodically. Moreover, a represented volume element along laminate thickness can be easily chosen, and its elastic constants are accurately acquired by employing a proper microscopic mechanical model. Herein, the high-fidelity generalized method of cells, which represents a good ability in predicting the effective mechanical behaviors of composites, was used. On this basis, the classic laminate theory was improved to suitable for describing the elastic constants and failure strength the SMPC with respect to ambient temperature. Numerical results show a good consistency to the experimental data. Moreover, a higher ambient temperature tends to sharply decrease their final failure strength. It is also revealed that the presented modeling method shows a great potential in calculating the effectively mechanical property of the natural short fiber-reinforced composites.

1. Introduction

Due to its unique molecular structure and good performances, including light weight, low cost, and high strain recovery rate, shape memory polymer (SMP) exhibits a widespread prospect in medical, energy, and electronic communications [1, 2]. It is reported by Leng et al. [3] that the maximum strain of the SMP reach to 600%. However, its low stiffness and recovery stress seriously limited its application in some extent. To overcome the flaws mentioned above, some reinforced phases, such as carbon fibers, carbon nanotubes, and nanoparticle, were presented to be mixed with the SMP to form the shape memory polymer composites (SMPC) [4]. Due to its intrinsic characteristic of the slow unfolding, the SMPC have been widely prepared as some deployable components in the satellite equipment, such as hinge and truss structure [5]. According to a statistical result by Castet and Saleh [6], the failure of solar wing during the deployment procedure is accounting for 17% of the total failure events, which results in a final failure of the satellite. To maximum guarantee the safe operation of the satellite equipment, it is critical to grasp the effective mechanical property of the SMPC with respect to the ambient temperature.

A large number of experimental tests have been executed on the standard specimen or structure of the SMPC to acquire its mechanical [79], electrical [10, 11], and shape memory performances [12]. The recovery capable of the SMPC under various loads was investigated by Basit et al. [13], and it exhibits a good potential in diverse applications. Le and Goo [14] executed on a folding and deployment test of the SMPC hinge at -10°C to investigate its shape recoverability and revealed that its deployment performance is closely dependent on ambient temperature. To reveal the temperature influences, an infrared camera and thermocouples were employed to observe the temperature distribution. By using a series of experimental methods, Wang et al. [15] studied the thermodynamic and shape memory properties of the SMPC. Moreover, its fracture interface was investigated by scanning electron microscopy. It is indicated that a good mechanical property can be acquired when fiber mass fraction is equal to 8%. Dao et al. [16, 17] measured the moment and blocking force of the SMPC hinge by employing a pulley-mass system. The test method provides a guideline for evaluating its mechanical performances. Annin et al. [18] investigated the deformation and damage of the SMPC during the bending and torsion. Moreover, the fiber off-axis angle influences on the formation process was also revealed. Lelieveld et al. [19] executed on the thermomechanical tests to investigate the actuation characteristics of the SMPC. Kim et al. [20] explored the recovery and physical properties by using the experimental tests on the strength and thermal conductivity.

In recent years, a large number of theoretical methods [21] have been presented to investigate elastic constants, nonlinear deformation, and failure characteristic of the SMPC. With an analogy way with thermomechanical responses at macroscale, Wang et al. [22] presented a chain model to optimize the chain system of the SMPC. It is revealed that the shape memory property is improved by increasing the cross link density. Arvanitakis [23] presented a constitutive level-set model to capture its mechanical behaviors under the thermomechanical cycle loading. Sun et al. and Gu et al. [24, 25] presented a multiscale strategy and established the thermoviscoelastic constitutive equation to investigate the buckling critical stress under the finite deformation, and the numerical results lay a solid foundation for its application and design. Bergman and Yang [26] proposed a macroscopic model combining a nonlinear geometric model with a temperature-dependent constitutive equation to investigate the shape fixation process. Li et al. [27] established a three dimension model by employing the ABAQUS to study the vibration mode and natural frequency of the large spatial deployable structure. The modal test results show a good consistency with the simulation results. Based on the homogenization scheme, Song et al. [28] presented a constitutive model for the SMPC subjected to the thermomechanical load. The experimental tests were executed on the standard specimens for a comparison.

For the SMPC, the mechanical property is closely dependent on its microscopic feature, including inclusion arrangement, interfacial bonding, and inclusion morphology. It is hardly for a macroscopic model to reveal the microscopic characteristic influences. A comprehensive investigation with full consideration of the microscopic factors influence has been executed by some researchers. By choosing a cylindrical representative volume element (RVE), Khalili et al. [29] evaluated the interface property influences on mechanical behaviors of the SMPC. Liu and Jiang [30] presented a novel hierarchical micromechanic method to study the microstructural feature influences on the thermal conducting behavior of the SMPC. However, few studies refer to random inclusions in the matrix materials. The main objective of this paper is to focus on investigating the effective mechanical property of the SMP reinforced with natural short fibers. To this end, an effective modeling scheme based on the microscopic constitutive equation was presented. The outline was summarized as follows: Section 2 presented the modeling scheme and microscopic modeling process of the SMP reinforced with natural short fibers. To validate the effectiveness of the presented method, experimental data and numerical results were both introduced for a comparison in Section 3. On this basis, a series of simulation analysis were executed to fully study the failure strength of the SMPC with respect to ambient temperature in Section 4. Conclusions are shown in Section 5.

2. Modeling Process of the SMPC

2.1. The Equivalent Procedure of the SMPC

For the shape memory polymer composites, short fibers always present obvious random characteristic due to the limitation of the preparation technology as shown in Figure 1(a). In the coordinate system , the short fiber direction can be accurately determined with respect to the fiber angles of and as shown in Figure 1(b).

It is reported by Cai et al. [31, 32] that the random short fibers exhibit an obvious directionality and closely dependent on resin flow direction. To simplify the microscopic modeling procedure, the random short fibers are always considered to be periodically distributed in matrix materials according to the homogenization theory as shown in Figure 2(a). According to the statistical results derived from the microscopic morphology, the random short fibers are equivalently represented by the ellipsoidal model. In addition, the fiber direction and density along thickness direction represent obvious differences. In details, few fibers can be found in the upper and lower edges, and the short fiber density in the center location is much higher. Therefore, the equivalent structure of the SMPC can be considered as the composite laminate (Figure 2(b)), whose upper and lower edges consist of pure polymer, while other regions can be treated as serials of unidirectional (UD) laminas. To sum up, the SMPC are accurately discretized as ten UD laminas, and eight of them are reinforced with ellipsoidal fibers. It should be noted that the elastic modulus of the natural fiber composites along is identical to the equivalent composites.

2.2. Numerical Modeling Process

It is assumed that the SMPC is equivalent to the composite laminates, which consist of serials of the UD laminas mentioned above. The equivalent ellipsoidal inclusions are assumed to be periodically distributed in the SMP. The simplified numerical procedure for the SMPC can be summarized as follows: (1) a proper RVE in each lamina can be easily selected along the thickness direction as shown in Figure 2(c). On this basis, the effective modulus of each lamina can be obtained by employing an effective microscopic mechanical model. (2) The SMPC consists of a series of unidirectional laminas. Therefore, the effective elastic modulus and mechanical property of the SMPC can be acquired by employing the classic laminate theory. (3) The numerical results by the presented method and experimental data are both employed for a comparison. On this basis, the failure strength of the SMPC with respect to ambient temperature is investigated by implanted the failure criterion.

2.2.1. Microscopic Constitutive Relation

By employing a microscopic mechanical theory, such as Eshelby equivalent inclusion theory [33], Mori-Tanaka method [34], self-consistent model, and generalized self-consistent model [35], the elastic constants of the composites can be acquired. However, they limit to analyze some inclusions with special geometric shape. It is hardly for the microscopic theories mentioned above to reveal the relation between geometric shape and inclusion distribution influences on the effective modulus. In recent years, the high-fidelity generalized method of cells (HFGMC) represented a good ability to investigate nonlinear deformation [36] and failure modes [37] of the composites under a coupled multifield environment. Ye et al. [38] presented a multiscale modeling framework by combining with the finite element method (FEM) to reveal the failure mechanism of the composites from microscopic damage to macroscopic fracture. Herein, the HFGMC is employed to calculate the effective modulus of the RVE reinforced with a spheroidal inclusion. During the numerical modeling process, the RVE as shown in Figure 2(c) is discretized as subcells. It should be noted that the parameters , , and indicate the discretized subcell number along , , and , respectively. The subscripts , , and indicate the numbered subcell. The second-order subcell displacement component in the RVE is expressed as the function of the fluctuating average displacements and the higher-order displacement components , that is, where , , and denote the subcell dimension along , , and (Figure 2(c)), respectively.

With full consideration of the displacement continuity condition between the adjacent subcells, the relation between the subcell strain components and average strains can be expressed as follows [36, 39]: where the symbols , , and denote the dimension of the RVE.

Similarly, the normal stress continuity between the adjacent subcells are written as: where indicates subcell stress components, which can be simplified as the function of .

With respect to the axial symmetry characteristic, that is, , the sub-cell shear stress components in the row are expressed:

In the row, the shear stresses can be simplified as follows:

Similarly, the shear stress components and can be simplified as follows:

It is revealed that the matrix materials always represented inelastic deformation. Moreover, the thermal residual stress can be easily discerned due to the mismatch of thermal expansion coefficients between the inclusion and matrix materials. Therefore, the average subcell strain components in the RVE are written as the function of the subcell stress , that is, where and represent the subcell strain and flexibility matrix, respectively. and are the subcell inelastic strain and thermal expansion coefficient. is the ambient temperature variation.

In Equation (16), the subcell thermal residual stress is determined by the temperature variation and thermal expansion coefficients of the constituent materials. Inelastic subcell strain can be calculated according to the loading step size. It is indicated that the subcell average strain can be effectively calculated once the subcell averaged stress is determined. From Equations (8)–(10) and Equations (13)–(15), the subcell normal stress and shear stress components are solved.

Substituting Equations (8)–(10) and Equations (13)–(15) into Equation (16), the subcell average strain components are obtained. Combining with the displacement continuity condition in Equations (2)–(7), the subcell stresses can be written as matrix form, that is, where and indicate the macroscopic average strain and inelastic stress components. denotes the thermal expansion coefficients of the constituent materials. The symbol is the stiffness matrix, which consists of subcell dimension.

Once subcell stress components in Equation (17) is obtained, the macroscopic average stress can be easily solved according to the homogenization theory, given by,

Substituting Equation (17) into Equation (18), the macroscopic average stress is written as the function of macroscopic average strain and inelastic strain , that is, where is the stiffness matrix of the composites.

2.2.2. Microscopic Constitutive Relation

Herein, the SMPC is considered to be composed of some UD laminas. Once the elastic modulus and mechanical behavior of each UD lamina are determined, the overall property of the SMPC can be obtained by employing the classic laminate theory [40, 41]. For each UD lamina, the stress-strain relation in the global coordinate system with respect to inelastic strain and thermal strain is improved as follows: where , , and indicate the stress components in the global coordinate system. , , and are the midplane strain components in the composite laminates. and are the thickness and stiffness matrix of the laminate. , , and are related to the partial differentiation of displacement to the global coordinate.

The total thickness of composite laminates is consider as the sum of each UD lamina, that is, where the symbol is the laminate number.

The internal force components , , and as well as the internal force moments , , and in each lamina can be acquired by executing on an integral operation along the thickness, that is, where and are the bend deflection rates. is the distorting rate in the midplane.

It can be found that the internal strain, midplane bend deflection, and distorting rates are unrelated to the laminate thickness. Equations (22) and (23) can be further simplified. The constitutive relation of the composite laminates with respect to inelastic strain and thermal strain components can be written as: where the symbol indicates the in-plane stiffness matrix. and are the coupling stiffness and bending stiffness matrices. The superscripts and are related to the inelastic and thermal components.

3. Model Validation

The studied SMPC is prepared with T-700SC short carbon fibers and trans-1, 4-polyisoprene (TPI) matrix. The average length and diameter of the carbon fiber are 2 mm and 7 μm, respectively. To verify the proposed method, the next two procedures are executed: (1) the ellipsoid fibers are assumed to be periodically distributed in the cubic matrix to equivalently replace the SMPC reinforced with natural short fibers. On this basis, a RVE with the ellipsoid inclusion can be easily selected, which can be easily discretized by the hexahedral subcells as shown in Figure 2(c). Moreover, its elastic modulus is calculated by the HFGMC and the FEM for a comparison. (2) A novel modeling procedure is proposed to acquire effective property of the SMPC with respect to ambient temperature, and the numerical results are compared with the experimental data. The constituent material parameters of the T700SC inclusion and TPI matrix at room temperature, including elastic modulus and Poisson’s ratio, are shown in Table 1.

3.1. A Comparison between Numerical Methods

To validate effectiveness of the microscopic mechanical theory, a sphere inclusion whose center coincides with the cubic RVE is considered. Figure 3 indicates the elastic modulus of the composites with respect to fiber volume fraction (FVF). For a comparison, numerical results by the FEM are also indicated in the figure when the FVFs are 0.2 and 0.3. It should be pointed out that the identical subcell number are used to discretize the RVE during the numerical modeling. It can be easily found that the numerical results by the 3D HFGMC exhibit a good consistency with the FEM. In addition, it can be easily discerned that the elastic modulus exhibit an exponential growth with the variation of the FVF. In details, elastic modulus is 87.7 MPa when the is considered, while the elastic modulus approximate to 478.8 MPa and 793.2 MPa when the FVF are 0.48 and 0.5, respectively.

To further investigate the inclusion shape influences on the elastic modulus, axial length of the ellipsoid inclusion is investigated. Herein, the RVE dimension of and the are considered in the example. The parameters , , and denote the dimension of the RVE along , , and as shown in Figure 2(c), respectively. It can be easily found from Figure 4 that longitudinal and transverse modulus represented a nonlinear variation with respect to axial length ratio of the ellipsoid inclusion. Moreover, transverse modulus exhibits a sharp decrease with an increase of the axial length ratio. In details, the transverse modulus is approximate to 105.1 MPa when the ellipsoid inclusion dimension is determined. It should be pointed out that the parameters , , and denote the ellipsoid fiber dimension. However, the transverse modulus is equal to 85.1 MPa when the ellipsoid inclusion dimension is considered. However, the longitudinal modulus is 81.3 MPa when the ellipsoid inclusion dimension is considered.

3.2. A Comparison with the Experimental Data

It is reported by Cai et al. [31, 32] (2021, 2022) that the FVF long thickness exhibits an obvious difference. In other words, the fiber distributed in the UD lamina close to the laminate surface is relative fewer. Therefore, it is critical to determine the FVF in each UD lamina, which can be found in Table 2. Herein, the FVF of the SMPC from 5% to 9% are considered in the example. The SMPC is further divided into 10 UD laminas, and the FVF in the surface UD laminas is different. Herein, the symbol “#” indicates the lamina number along thickness direction as shown in Figure 2(b). In addition, the ellipsoid fiber content should also be determined. According to the experimental data at room temperature, the ellipsoid fiber dimension in the 5% SMPC is determined by . Elastic modulus at an elevated temperature is also considered in the example. Herein, it should be noted that the elastic modulus of matrix materials is temperature-dependent, whose material parameters are provided by Zeng et al. [42, 43]. It should be noted that Poisson’s ratio of the constituent materials is considered to be temperature-independent. In other words, this parameter of the constituent materials is a constant without respect to the ambient temperature influence during the numerical calculation. With respect to the raised ambient temperature, a sharp decrease of the elastic modulus for the SMPC is easily discerned as shown in Figure 5. This is attributed to the decreased modulus of matrix materials. Moreover, numerical results represent a good consistency to the experimental data. Overall, numerical results by the presented method exhibit a higher accuracy than the method proposed by Zeng et al. [42]. In details, an approximate results between the present method and the Zeng’s method can be found at 301 K. The maximum error can be found when the ambient temperature is equal to 334 K, and the present method exhibits a much higher accuracy than the method proposed by Zeng et al.

Similarly, procedures are also executed on the 7% and 9% SMPC to determine the ellipsoid fiber dimension at room temperature, and elastic modulus with respect to ambient temperature is also investigated as shown in Figures 67. The ellipsoid inclusion dimension are and , respectively. Combing with the , an obvious difference of the inclusion dimension can be found. A sharp decrease of elastic modulus can be found when a higher ambient temperature is considered. Compared with experimental data, it is revealed that numerical results provided by the present method exhibit a higher accuracy than the theoretical method provided by Zeng et al. [42]. Combined with the numerical results as shown in Figures 46, it can be concluded that the FVF will contribute to an increase of the elastic modulus. In details, the elastic modulus is approximate to 86.99 MPa at the room temperature when the . While the elastic moduli increase to 91.12 MPa and 93.32 MPa when the FVF is 7% and 9%, respectively.

4. Failure Strength Investigations

It is critical for researchers to grasp the failure property of the SMPC in service to maximum their usage life. However, the failure characteristic are closely dependent on many factors, including the FVF, fiber arrangement, and constituent material property. In addition, the service environment also plays an important role in the mechanical property. For instance, the SMPC are exposed to a high temperature condition, which may seriously reduce the failure strength. To fully grasp ambient variation influences, the failure strength at the room temperature and the elevated temperatures 321 K and 341 K are considered. To further valuate the failure strength of the SMPC with respect to the FVF variations, the SMP mixed with 5% and 7% natural short fibers are also investigated. It should be noted that the discretize mode of the composite laminate and inclusion shape in the matrix materials are identical with the study in Section 3. Moreover, the inelastic deformation derived from the matrix materials and thermal residual stress derived from the preparation process are both ignored.

Figure 8 indicates the failure strength of the SMPC with respect to the FVF variations at three ambient temperatures, that is, 301 K, 321 K, and 341 K. Herein, the failure property of the constituent materials at different ambient temperatures are derived from the study by Guo et al. [43]. The fiber longitudinal and transversal failure strengths are 4900 MPa and 80 MPa, respectively. During the prediction of the final failure strength, the maximum stress criterion [44] is employed. It should be pointed out that the failure strength of matrix is temperature-dependent, while the fiber strength always equal to a constant in evaluating the final failure strength at different ambient temperatures. For a further validation, experimental data and numerical results of the failure strength are both indicated. In general, a good consistency between numerical results and experimental data can be found in the figure, and the increased short fibers tend to improve the failure strength of the SMPC in some extent. This is attributed to a higher failure strength of the natural fibers. The maximum error 6.34% can be found under the 301 K when the is considered. In addition, an elevated temperature tends to sharply decrease the final failure strength of the SMPC. In details, the failure strength is approximate to 16.39 MPa when the ambient temperature 301 K is considered. However, the failure strengths decreased to 6.61 MPa and 3.59 MPa when the ambient temperatures 321 K and 341 K are considered. This is attribute to a huge reduction of the elastic modulus with an increase of the ambient temperature. In other words, the lower elastic modulus of the constituent materials sharply decrease the stiffness behavior of the SMPC. Naturally, their failure strength is reduced to some extent.

5. Conclusions

In this study, mechanical behaviors of the SMPC reinforced with natural short fibers were investigated by the presented constitutive model. The SMPC represents a huge deformation in service, and the improved laminate theory is suitable to describing its deformation and inclusion fraction influences on the SMPC. Numerical results of the elastic constants show a good consistency with the experimental data, which demonstrated the improved constitutive model. The conclusions are summarized as follows: (1)The SMPC reinforced with natural short fibers can be considered as the composite laminates, which are composed of the UD laminas with ellipsoid inclusion, whose inclusion dimension can be easily determined according to the experimental data at room temperature(2)The present method exhibits a good accuracy in investigating elastic modulus, which represents a sharp decrease with respect to an elevate temperature. For the 5% SMPC, the ellipsoid fiber dimension represents a good consistency with experimental data. When the SMPC with 7% and 9% short fibers are considered, the ellipsoid inclusion dimensions are and , respectively(3)The final failure strength of the SMPC can be effectively acquired by the proposed ellipsoidal inclusion model, and a higher ambient temperature results in a sharp decrease of their failure strength

The investigations in this article provide a new rapid evaluation method in predicting the mechanical property of natural short fiber composites, while the proposed method limits to a small FVF. A higher FVF gives rise to a negative influence on the failure strength, which may be attributed to the variation of interfacial property.

Data Availability

The raw/processed data required to reproduce these findings cannot be shared at this time due to technical or time limitations.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Authors’ Contributions

Yang Li is responsible for writing—original draft and methodology; Junjie Ye for writing—original draft, methodology, and supervision; Lu Liu for validation, visualization, and software; Baoquan Shi for methodology, writing—review and editing, and supervision; and Yumin He for validation and visualization.

Acknowledgments

This work was supported by the National Natural Science Foundation of China, China (nos. 52175112 and 51675397); the National Natural Science Foundation of Shaanxi Province, China (no. 2018JM5029); and the Fundamental Research Funds for the Central Universities (JB210421).