Next Article in Journal
Preparation and Comparison of Properties of Three Phase Change Energy Storage Materials with Hollow Fiber Membrane as the Supporting Carrier
Next Article in Special Issue
Finite Element Study on the Impact Resistance of Laminated and Textile Composites
Previous Article in Journal
Synthesis of Highly Thermally Stable Daidzein-Based Main-Chain-Type Benzoxazine Resins
Previous Article in Special Issue
Effect of Porosity on Mechanical Properties of 3D Printed Polymers: Experiments and Micromechanical Modeling Based on X-ray Computed Tomography Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modelling the Mechanical and Strain Recovery Behaviour of Partially Crystalline PLA

IRC in Polymer Science and Technology, Mechanical and Energy Systems Engineering, Faculty of Engineering and Informatics, University of Bradford, Bradford BD7 1DP, UK
*
Author to whom correspondence should be addressed.
Polymers 2019, 11(8), 1342; https://doi.org/10.3390/polym11081342
Submission received: 4 July 2019 / Revised: 7 August 2019 / Accepted: 9 August 2019 / Published: 13 August 2019
(This article belongs to the Special Issue Mechanical Behavior of Polymers)

Abstract

:
This is a study of the modelling and prediction of strain recovery in a polylactide. Strain recovery near the glass transition temperature is the underlying mechanism for the shape memory in an amorphous polymer. The investigation is aimed at modelling such shape memory behaviour. A PLA-based copolymer is subjected to stress–strain, stress relaxation and strain recovery experiments at large strain at 60 °C just below its glass transition temperature. The material is 13% crystalline. Using published data on the mechanical properties of the crystals, finite element modelling was used to determine the effect of the crystal phase on the overall mechanical behaviour of the material, which was found to be significant. The finite element models were also used to relate the stress–strain results to the yield stress of the amorphous phase. This yield stress was found to possess strain rate dependence consistent with an Eyring process. Stress relaxation experiments were also interpreted in terms of the Eyring process, and a two-process Eyring-based model was defined that was capable of modelling strain recovery behaviour. This was essentially a model of the amorphous phase. It was shown to be capable of useful predictions of strain recovery.

Graphical Abstract

1. Introduction

Polylactic acid (PLA) and related polymers and blends are a focus of increasing attention, arising both from their potential for degradability within the natural carbon cycle for disposable products [1,2], and from their bioresorbable properties for prosthetic implants [3,4,5,6]. In the latter category, there are applications for ‘smart’ implants that make use of the polymer’s shape memory properties [7]. As a result, of the load bearing nature of many applications, both potential and realised, it has become important to understand the mechanical behaviour of PLA. This is the motivation of this paper, for which a fundamental study was made of the mechanical behaviour of PLA, with particular emphasis on aspects relevant to shape memory and strain recovery.
The physical and mechanical properties of PLA materials have recently been reviewed by Farah et al. [8] and Garlotta [9], the latter with more emphasis on crystalline properties. Bergström and Hayman [10] have produced a review concentrating more on mechanical properties with relevance to medical applications. PLA materials resemble most other thermoplastics in that they are mechanically nonlinear, time-dependent and have varying levels of crystallinity. At elevated temperatures near the glass transition of around 60 °C, they are capable of being deformed to large strains, acquiring molecular orientation and hence shape memory properties. There have been two notable recent studies of molecular orientation effects arising from biaxial deformations [11,12].
PLA materials are thermoplastics that can be processed using well-established conventional methods including extrusion, injection moulding, injection stretch blow moulding, film or sheet casting and thermoforming [2]. They have also been made into foams to produce superhydrophobic devices for oil–water separation [13]. Additionally, they are a leading polymer in additive manufacturing [14]. To induce shape memory properties into the material, a manufacturing process that creates high levels of molecular orientation is required. Such a method is die drawing, originally developed to produce polymer products with high stiffness and high strength properties [15,16,17]. This process has been successfully applied to PLA-based polymers [18,19] providing a process with potential for shape memory polymer production.
In addition to applications of the pure polymer, PLA has also been used as the principal component of biocompatible blends and composites. Vasile et al. [20] have plasticized PLA with polyethylene glycol (PEG) and used melt-mixing to combine it with biological materials including chitosan in various combinations. The resulting materials showed useful levels of antioxidant and antibacterial activity together with good biocompatibility, with potential for both biomedical and packaging applications. PLA/PEG systems also have a role in the production of polymer-nanoparticles hydrogels [21]. Well-dispersed graphene nanocomposites have been produced using poly(butylene adipate-co-terephthalate) (PBAT) using a PEO/PLA/PBT system [22].
This study is motivated by growing developments in biocompatible shape memory polymers in biomedical applications [23]. Recent examples of specific applications of shape memory PLA systems range from in-vivo devices made from shape memory nanocomposites [24] to melt-blown nonwoven fabrics with potential for load bearing implants, scaffolds and wound dressings [25]. The level of crystallinity is in general an important factor associated with high stiffness and strength, and nucleation techniques have been developed to enhance it in PLA systems [26]. However, in the biomedical context, it has been established that the bioresorbability of poly l-lactide (PLLA) is too low due to its high crystallinity [27], and so less crystalline forms are used. Such a material, a polylactide with a relatively small (~13%) level of crystallinity, is the subject of the experimental study in this paper. The objective is to develop a material model that predicts strain recovery after tensile extension, to form the basis of a model of the shape memory effect. Shape memory is driven by the entropic stresses in the amorphous phase of the polymer. To understand this mechanism on the basis of mechanical experiments on the partially crystalline polymer, it is necessary to assess the effects of the crystalline phase on the measurements. This we do by modelling in a way that shows a significant effect of the crystal phase on the macroscopic elastic properties. Extension of the modelling and comparison with stress–strain experiments allow us to isolate the yield behaviour of the amorphous material and define its strain rate dependence. A constitutive model is then constructed that enables the stress relaxation and recovery behaviour to be simulated and compared with experimental observations of strain recovery. The experiments are a direct measure of shape memory effectiveness and thus an appropriate means of evaluating the model. A testing temperature of 60 °C is chosen as, being slightly below the glass transition temperature, it makes possible high levels of molecular orientation and also enables subsequent strain recovery, replicating shape memory behavior.
The material modelling to assess the effects of the crystalline phase is based on the concept of crystallites as inclusions within an amorphous matrix. This approach to the modelling of semicrystalline polymer was pioneered by Halpin and Kardos [28], who made use of the Halpin–Tsai equations [29]. Halpin and Kardos used an approach that acknowledged the effect of the stiffness of crystals relative to that of the amorphous phase. This gave better predictions than previous models, which had treated crystals as merely sources of more crosslinking sites. Working with crystal volume fractions of up to 30%, they used the Halpin–Tsai composite theory to produce successful models of the stiffness of polyolefins. There is a potential objection to the modelling of crystallites as inclusions isolated from each other, which Halpin and Kardos discussed; the existence of tie-molecules. Partly as a result of the success of their modelling, they describe the contribution of tie-molecules as “at best, of only secondary importance” [28]. It is clear that, at small material strains, tie-molecules will behave no differently than other molecules in the amorphous phase. They will only have an effect of linking crystals to form a network when they are highly extended at large strains, and we are using this form of modelling only up to a strain of 0.025. As a consequence of these arguments, we shall proceed with modelling the material as isolated crystallites within an amorphous matrix.

2. Materials and Methods

The material is a copolymer of l-lactide and dl-Lactide in a 70:30 molar ratio. It was supplied in the form of granules and processed using melt extrusion into hollow tube to be used as specimens for mechanical testing.
Before processing the polymer, granules were dried overnight in a vacuum oven. A single-screw Dr Collin extruder (Maienbeth, Germany) with a 16 mm diameter screw and the final heating zone at 180 °C was used to extrude the tube material. Granules were supplied into the hopper, in which nitrogen was circulated to minimise moisture. The extrudate was pulled through a cooling bath kept at 20 °C by means of a caterpillar. The tube was 3.2 mm in external diameter and 0.8 mm internal diameter.
DSC testing was performed on the extruded material using a TA Instruments (New Castle, DE, USA) Discovery DSC Instrument operating at a heating rate of 10 °C/min. From these results the crystallinity was calculated as 13.1% using the method of Kong and Hay [30], and the glass transition temperature estimated as 64.3 °C. The melt temperature was 165 °C.
Molecular mass of the tube material was determined using gel permeation chromotography (Waters Acquity Advanced Polymer Chromatography System, Milford, MA, USA). The values were measured as M ¯ w = 650,000 Da and M ¯ n = 370,000 Da.
Mechanical tests were carried out in a fan oven at 60 °C using an Instron 5565 testing machine (High Wycombe, UK). Tensile test specimens consisted of lengths of tube with gauge length 75 mm. Specimens were stretched at a range of constant speeds to a strain of 1 (extension ratio 2). For all experiments, tensile stretching was performed at 60 °C at constant speed with strain rates in the range of 2.2 × 10−3–4.0 × 10−2 s−1. Tensile tests were of two kinds: stress relaxation, for which the extension was held constant and the load monitored for a set time interval; and strain recovery, for which one end of the specimen was released from the testing machine on attaining the set extension and strain then allowed to recover under fixed load. For the latter tests, the specimen was released by means of a spring-loaded pin and strain recovery was monitored using video extensometry through the oven window. Images were captured using a PixeLINK model PL-D722MU-T video camera (Rochester, NY, USA), and images were analysed using the Fiji image processing software package based on ImageJ (Eliceiri, University of Wisconsin-Madison, Madison, WI, USA) [31].
Three-point bend tests were carried out on lengths of tube specimen loaded centrally on a span of 25 mm using the same testing machine, oven and temperature using. Testing speeds were such as to provide strain rates close to those in the tensile tests.

3. Results

3.1. Small Strain Structural Modelling

3.1.1. Analytical Approach

The Mori-Tanaka model [32,33], based on the exact solution for an ellipsoidal inclusion within a continuous medium, is a more rigorous approach to composite modelling than the semi-empirical Halpin–Tsai method. The Mori–Tanaka composite is assumed to consist of regularly spaced ellipsoidal inclusions. Each ellipsoid is a spheroid with two equal axes, axisymmetric at the third axis (the symmetry axis). The inclusions are arranged with the symmetry axes aligned parallel to one another. The aspect ratio a is the ratio of the length of the symmetry axis to that of the other two axes; thus, both disc-like (a < 1) and rod-like (a > 1) geometries can be defined. Both matrix and inclusions are of isotropic elastic material. We have adopted the approach to the Mori–Tanaka theory established for ellipsoidal inclusions [34,35].
The amorphous and crystal phases are assigned elastic moduli Ea and Ec respectively, with corresponding Poisson’s ratios νa and νc. As a result of the alignment of the ellipsoids, the composite as a whole is orthotropic, with modulus along the direction of the ellipsoids’ symmetry axes E given by
E E a = A A + ϕ ( A 1 + 2 ν a A 2 )
and the modulus normal to it by
E E a = 2 A 2 A + ϕ ( 2 ν a A 3 + ( 1 ν a ) A 4 + ( 1 + ν a ) A 3 A )
Here ϕ is the volume fraction of inclusions and A, A1, A2, A3, and A4 are functions of Eschelby’s tensor, the elastic constants of both phases, and ϕ [34,35].
As discussed further in Section 3.3.1 below, normal moduli of the α-phase PLLA crystals have been measured using X-ray techniques [36]. Strains in the crystals are deduced from the change in lattice spacing as measured by X-ray diffraction when the material is subject to a known stress. For the a, b and c axes of the orthorhombic crystal, the Young’s moduli are respectively 3.2, 2.8 and 14 GPa, where the chain axis is along c. In a lamellar crystal structure, the chain axis is normal to the plane. In terms of the Mori–Tanaka model, such a structure would be best approximated by disc-like ellipsoids, with the symmetry axis of each ellipsoid identified with the chain axis, which also corresponds to the direction of the largest modulus. However, the theory is limited to isotropic ellipsoids and so we shall explore the extremes Ecrystal = 3 and Ecrystal = 14 GPa modulus in isotropic ellipsoids. In Figure 1, we present predictions of moduli in terms of the ratios of Equations (1) and (2) for disc-like inclusions over a range of aspect ratios a < 1. Here we assume that the amorphous modulus Ea = 0.5 GPa, consistent with the modulus of the bulk material in 3.2.1 below. We take the volume ratio to be ϕ = 0.13 as measured for our material by DSC, and Poisson’s ratios of 0.4 for both materials.
It is clear that the crystals can have a significant effect on the modulus at this concentration. This is consistent with the dynamic mechanical measurements of Tábi et al. [37]. However, in this model the alignment of the ellipsoids ensures that the modulus increase is accompanied by an unrealistic level of anisotropy as defined by E and E . The choice of modulus also has significant effects. This suggests that the issue should be pursued further using a model in which the inclusions are elastically anisotropic and randomly aligned.

3.1.2. Numerical Approach

For a more realistic method, we have created three-dimensional finite element models of the partially crystalline structure. Solid orthorhombic regions representing the crystallites are embedded within a cube of material representing the amorphous polymer, to give a representative volume element (RVE). Model crystallites of identical size and shape are added to the structure at random positions and at random orientations about one of the global axes, and are ensured to be non-intersecting using an acceptance–rejection algorithm [38]. Thus, the orientation of each crystallite is semi-random, in that one set of edges is parallel to the 3 axes but otherwise the shape is angled randomly in the 1–2 plane; see Figure 2. Orthotropic elastic properties based on X-ray measurements on PLA crystals [36] are assigned to the model crystallites. To simulate an infinite body and avoid unrealistic deformations at the RVE boundaries, periodic displacement boundary conditions are applied at the cube surfaces. The overall stress response is evaluated using the ABAQUS finite element package (ABAQUS 6.14-2, Dassault Systèmes, Johnston, RI, USA). Small-strain elastic properties of the composite solid are derived using the stress responses to appropriate boundary conditions applied to the models averaged over repeated realisations. Uniaxial yield properties are also explored at larger strains using the same RVEs with elastic–plastic matrix properties. This enables the overall yield response of the composite solid to be quantitatively related to the yield point of the matrix.
Periodic boundary conditions are enabled by designing the RVE mesh so that there are matching node positions on opposite boundary faces. For an RVE of side length L, for boundary faces normal to the 1 direction X 1 = 0 and X 1 = L the nodal positions are such that for every node at ( 0 , X 2 , X 3 ) there is a corresponding node at ( L , X 2 , X 3 ) and vice versa. Similar conditions apply for the pairs of faces X 2 = 0 ,   X 2 = L and X 3 = 0 ,   X 3 = L .
To evaluate elastic properties, normal moduli are calculated by deforming the RVE along each global direction in turn while restraining boundary planes parallel to this direction of deformation. For an extension Δ L along 1, the displacements u 1 at the boundary planes normal to 1 are related by
u 1 ( L , X 2 , X 3 ) u 1 ( 0 , X 2 , X 3 ) = Δ L
for each pair of corresponding nodes. Plane strain type restraints in the 2 and 3 directions are defined by
u 2 ( X 1 , L , X 3 ) u 2 ( X 1 , 0 , X 3 ) = u 3 ( X 1 , X 2 , L ) u 3 ( X 1 , X 2 , 0 ) = 0
The systems of Equations of (3) and (4) are programmed into ABAQUS. Analogous systems are created for plane strain deformations of the same magnitude along 2 and 3. Each randomly generated realisation of the model is subject to the same three deformations along 1, 2 and 3. The macroscopic stress responses σ 1 ,   σ 2 and σ 3 calculated in the analysis then define a stiffness matrix C for each realisation of the RVE. For the macroscopic strain given by ε = Δ L / L , Hooke’s Law gives for the 1, 2 and 3 directions, respectively,
σ 1 = c 11 ε ,   σ 2 = c 12 ε ,   σ 3 = c 13 ε σ 1 = c 12 ε ,   σ 2 = c 22 ε ,   σ 3 = c 23 ε σ 1 = c 13 ε ,   σ 2 = c 23 ε ,   σ 3 = c 33 ε
thus defining the components cij of C. The stiffness matrix C is inverted to give a compliance matrix S, the diagonal terms of which give the moduli E 1 ,   E 2 and E 3 along global directions 1, 2 and 3, respectively.
For a given set of conditions, several RVEs are generated and solved to give a set of modulus values, allowing a mean and standard deviations to be obtained. A specified crystal volume fraction φ corresponds to a total number N of inclusions. To represent the stochastic nature of the material, the number of inclusions n in each realisation is derived from the expected value N on the basis of a Poisson distribution, such that the probability of n inclusions is given by
p n = N n n ! e N
where the average of n tends to N over a large number of realisations.
The yield properties of the composite system have been explored using uniaxial deformations of selected RVEs. Here the inclusions retain their elastic properties, but the matrix is assumed to be elastic-plastic with a simple von Mises yield criterion. This takes the form
2 σ Y = ( ( σ I σ II ) 2 + ( σ II σ III ) 2 + ( σ III σ I ) 2 ) 1 / 2
for uniaxial yield stress σ Y and principal stresses σ I ,   σ II and σ III . This enables us to explore the yield behaviour of the RVE and how it relates to the underlying yield property of the amorphous polymer.

3.2. Experimental

3.2.1. Stress–strain Results

A full set of stress–strain curves is shown in Figure 3. Each curve is an average of three results. The degree of reproducibility can be gauged from the yield stresses, for which the coefficient of variation is 7% or less, and the stresses at the end of loading, for which the coefficient of variation is 4% or less. There is recognisable yielding and strain hardening at all rates, and a marked strain rate dependence of stress, the latter also noted by Söntjens et al. [39] for a similar material. Zhou et al. [40] produced stress–strain curves of PLA over a range of temperatures at a strain rate of 0.0013 s−1. At 60 °C they report stresses higher than those at the lowest rate in Figure 3, consistent with their PLA material being of higher crystallinity. The initial linear response was detectable and measurable up to very low strain (0.1–0.2%). Since the strain was calculated from the Instron crosshead movement, its accuracy was considered to be insufficient to give reliable Young’s modulus values. Therefore, the tests were supplemented with three-point bend tests as described above, for which larger crosshead motions are required to achieve the strains within the linear range. From these measurements the modulus was found to be 0.53 GPa.

3.2.2. Stress Relaxation

Initial loading was as described in Section 3.2.1 above, with the strain held constant for 1000–1400 s. Stress relaxation curves for the range of strain rates are shown in Figure 4. Each curve is an average of three results. For each set of three results, the coefficient of variation is less than 4% at the start of stress relaxation and has increased to be less than 7% by the end of the experiments.
There is little data on stress relaxation of PLA in the literature. Guedes et al. [41] report stress relaxation on a PLA-PCL (polycaprolactone) blend at room temperature and at strains up to 0.16, and were able to model the results using linear viscoelasticity. Duscunceli et al. [42] have studied stress relaxation of PLA at small strains and up to 50 °C. For conditions closer to those of this study, Sweeney et al. [43] report stress relaxation for a PLA-CaCO3 composite at a strain of 3.0 in the temperature range 60–70 °C, with initial loading at an average strain rate of 0.036 s−1. For these results at 60 °C, the stress at 1000 s has dropped from its initial value by 63%; this compares with the loading rate of 0.04 s−1 in Figure 4 in which the stress drops by 72% at 1000 s. On this basis, the stress relaxation behaviours are roughly equivalent.

3.2.3. Strain Recovery under Load

Initial loading was as described in Section 2 above, with a level of fixed load applied using a 5.9 N weight, corresponding to an engineering stress of 0.89 MPa. Images were captured at 1 s time intervals. Results are shown in Figure 5. Each curve is an average of three tests for which the coefficient of variation of the strain at 1400 s is 6% or less. Recovery levels are similar to those reported for a PLA-PBS blend when stretched to strains of 1.0, 2.0 and 2.5 [44].

3.3. Numerical Analysis of Results

3.3.1. Effects of Inclusions

For this method of modelling, set out in 3.1.2 above, results are not affected by the absolute size of inclusions. However, their aspect ratios can have significant effects. We acknowledge that in practice crystal inclusions will possess a range of aspect ratios, whereas in the current modelling all the inclusions are identical. We have varied the aspect ratio in order to explore its significance on the overall stress–strain behaviour of the RVEs. A typical lamella crystal has a thickness d equal to the chain length along the c direction (see Figure 2), and along the other directions normal to c, (i.e., a and b), there are many chain folds, corresponding to lengths greater than d. We have assumed that the dimensions normal to c are both equal to w. The aspect ratio is defined as d/w and we have adopted three values, d/w = 0.34, 0.17 and 0.086.
Orthotropic elastic properties are assigned to each inclusion based on the X-ray derived moduli of Lee et al. [36] for PLLA and PLDA crystals. Along the chain axis, they derived the value Ec = 14 GPa, which we assign to the c direction of the inclusion. In the transverse directions a and b, they measured 3.2 and 2.8 GPa. Since these values are similar we have assumed the inclusions to be transversely isotropic, with transverse modulus Et in the a–b plane given by Et = 3 GPa. Lee et al. found no significant change in modulus in the temperature range 13–300 K, and with a melt temperature of 165 °C for our material, the room temperature moduli of Lee et al. are appropriate for our 60 °C tests.
For the transversely isotropic inclusions the compliance matrix [45] is given in local a–b–c axes (see Figure 2) as
S = ( s aa s ab s ac s ab s aa s ac s ac s ac s cc )
for normal stresses and strains plus the shear compliance terms s abab , s bcbc which relate the shear stresses and strains:
ε ab = s abab σ ab ε bc = s bcbc σ bc ε ac = s bcbc σ ac
Since ab is an isotropic plane, we have the relation [45]
s abab = 2 ( s aa s ab )
The diagonal terms of Equation (8) are defined by
s aa = 1 E t s cc = 1 E c
and the off-diagonal terms depend on the Poisson’s ratios ν ba , ν ac :
s ab = ν ba s aa = ν ba E t s ac = ν ac s cc = ν ac E c
We have generated values for s ab and s ac using Equations (12) on the assumption that ν ba = ν ac = 0.4 . This then yields a value for s abab using Equation (10). No information is available for s bcbc . Two values for it are used of the same order of magnitude as s abab , but they differ from each other by 50% to investigate the effect of this variation on the RVE stress response.
The components of S generated in this way are
S = ( 0.333 0.133 0.029 0.133 0.333 0.029 0.029 0.029 0.071 ) ( GPa ) 1
for the normal stresses and strains and
s abab = 0.933   ( GPa ) 1 s bcbc = 1.0 ,   0.667   ( GPa ) 1
for the shear components.
For the ABAQUS analysis, stiffness components are required. These were obtained by inverting S and taking reciprocals of s abab and s bcbc .
A typical realisation is shown in Figure 6.
First, the method was used to estimate the modulus at small strains of the amorphous phase. Linear elastic analysis was used. Several realisations were implemented for each inclusion aspect ratio, evaluating a stiffness matrix C for each run using Equations (5). The moduli E1 and E2 of the RVEs were taken as representative of the effect of randomly oriented inclusions, in contrast with E3 which is in the direction along which the inclusions were perfectly aligned. For each inclusion, aspect ratio the average value of E1 and E2 was derived, initially by varying the amorphous modulus Ea on a trial-and-error basis to attain an average modulus consistent with the observed value of 0.53 GPa. As a result, we found that the value Ea = 0.36 GPa gave a set of averaged values as given in Table 1.
From these results we conclude that Ea = 0.36 GPa is an adequate estimate of the amorphous phase modulus at the test temperature of 60 °C. The effect of the variation in s bcbc on the overall modulus is within the range of error.
We have also used the RVE technique to explore the relationship between macroscopic stress–strain behaviour of the partially crystalline material and the yielding of the amorphous phase. As is apparent from Figure 3, deviation from linearity begins at low strains. For this analysis we have selected several RVEs and subjected them to uniaxial extension, having replaced the elastic material model of the amorphous phase with an elastic-plastic model governed by a von Mises criterion (Equation (7)) and allowing nonlinear deformations. We assume that the crystallite inclusions remain elastic as defined above. We assign values of uniaxial yield stress σ Y to the amorphous region of 1.0, 2.5, 5.0 and 10.0 MPa, covering the range of stresses observed in Figure 3 for the various strain rates.
In Figure 7, we show RVE results for the four values of yield stress. The macroscopic stress can be compared with the linear extrapolation of stress from low strains. For yield stresses 1.0, 2.5 and 5 MPa, there is a consistent relationship between the macroscopic stress, the yield stress and the linear extrapolation. For a macroscopic uniaxial stress σ R V E ( ε ) and linear extrapolation σ l i n ( ε ) = E ε , we can calculate a proportional difference p between the two for a given strain:
p = E ε σ R V E ( ε ) σ R V E ( ε )
At the particular strain when σ R V E ( ε ) = σ Y values of p are given in Table 2. In the range 1–5 MPa yield stress, p is essentially constant. This phenomenon provides a means of detecting yield in the amorphous phase from an experimental stress–strain curve: when the offset of the stress from the linear extrapolation is at this constant value of p, the experimental stress is equal to the yield stress. In this region of stress, p is on average 0.138. The deviations from linearity seen in Figure 3 occur within the range of stress 1–5 MPa, so it is appropriate to use this value of p to calculate σ Y for the six strain rates plotted there. The σ Y values thus derived are plotted against logarithm of strain rate in Figure 8. Except for the lowest two strain rates, for which the stress–strain curves are essentially indistinguishable, the yield stresses show a linear dependence on the logarithm of rate. This in contrast to the corresponding plot for the maximum experimental stresses, which shows a distinct upward curvature.

3.3.2. Modelling the Time-Dependent Yield of the Amorphous Phase

The linearity of response shown in Figure 8 suggests that yielding of the amorphous region is governed by an Eyring process [45,46]. In its simplest form, for uniaxial conditions, this process takes the form
ε ˙ p = α sinh ( v k T σ )
where ε ˙ p is the plastic strain rate, α a temperature dependent factor, v the activation volume, k the Boltzmann’s constant and T the absolute temperature. For convenience we write V = v k T , where V is the operational activation volume in units of reciprocal stress so that Equation (16) becomes
ε ˙ p = α sinh ( V σ )
To model yield, it is assumed that the Eyring process acts in series with an elastic element to form a Maxwell-like model. Then, as total strain on this model is increased, stress increases until the plastic strain rate, rising as a nonlinear function of stress, matches the applied strain rate and the elastic element ceases to stretch. Then stress remains constant at a value dependent on the applied strain rate. We may use Equation (17) to model yield by equating the applied strain rate with the plastic strain rate. Also, the sinh function can be approximated with an exponential function on the assumption that the argument of the sinh is large enough. Then Equation (17) becomes
ε ˙ p = 1 2 α exp ( V σ Y )
and can be re-arranged as
σ Y = 1 V ln ( 2 ε ˙ α )
so that the linear slope in Figure 8 can be interpreted. The slope in Figure 8 gives a value V = 0.454 MPa−1, corresponding to an activation volume v = 2.1 nm3. This can be compared with the value 0.72 nm3 obtained by Söntjens et al. [39] for combined results of PLLA, PDLLA and PLDLLA.
Additionally, the value of α can be obtained from the intercept of the fitted line in Figure 8 to give the value α = 0.00752 s−1. This Eyring process is now fully defined.
The analysis so far has only covered the time-dependent behaviour at small strain. The stress relaxation results provides data for large strain, and can be analysed in terms of the Eyring process. Guiu and Pratt [47], again using the model of an Eyring process in series with a spring, showed that the approximation Equation (18) leads to a form of stress relaxation curve for stress at time t after the state of constant strain has been reached at time t = 0:
σ ( t ) σ ( 0 ) = 1 V ln ( 1 + t c )
Here, c is a constant that involves the stiffness of the elastic element, which is assumed to be linear elastic. For large strain Sweeney et al. [43] have generalised the model so that the elastic element operates at large strain and is in the form of a Gaussian model. The Maxwell model is assumed to be stretched at a constant extension ratio λ which is split multiplicatively into elastic and plastic components λ e and λ p respectively, so that
λ = λ e λ p
The plastic strain rate is identified as
ε ˙ p = λ ˙ p λ p
Then, c is given by
c = 1 α VG ( 2 λ e 2 + λ e 1 ) exp ( V σ ( 0 ) )
where G is the strength of the Gaussian process and the stress in uniaxial conditions is given by
σ = G ( λ e 2 λ e 1 )
There is an approximation involved in the derivation of Equation (23), in that the quantity 2 λ e 2 + λ e 1 is assumed to be a constant; this means that the analysis is valid for slowly varying stress, and makes the evaluation of α from the fitted value of c difficult to justify. We have generated optimised fits of Equation (20) to the experimental stress relaxation curves of Figure 4 to give values of V and c. The theory set out above applies to stress relaxation following instantaneous loading and so we base the analysis on the highest strain rate of 0.04 s−1. The fitted and experimental curves are shown in Figure 9. The fitted parameter values are V = 2.2 MPa−1 and c = 3.47 × 10−7 s−1.
The Maxwell model comprising Gaussian and Eyring processes is to be used for the more general case of loading and strain recovery. For this analysis we require the Eyring stress as a function of plastic strain rate, which can be derived from Equation (17) as
σ = 1 V ln [ ε ˙ p α + ( ε ˙ p α ) 2 + 1 ]
This is in equilibrium with the Gaussian element, so we have
1 V ln [ ε ˙ p α + ( ε ˙ p α ) 2 + 1 ] G ( λ e 2 λ e 1 ) = 0
To solve this, we express it in terms of the plastic extension ratio λ p which is related to the plastic strain rate by Equation (22). Also we eliminate λ e by using Equation (21), obtaining
1 V ln [ 1 α λ ˙ p λ p + ( 1 α λ ˙ p λ p ) 2 + 1 ] G ( ( λ λ p ) 2 ( λ λ p ) 1 ) = 0
where λ is the total extension ratio applied to the model. Equation (27) is solved numerically using a time-incremental approach. At each increment, the value of λ p is calculated and the stress is then given by
σ = G ( ( λ λ p ) 2 ( λ λ p ) 1 )
To model the strain recovery experiments, two parallel Maxwell models are required (see Figure 10). The two arms are denoted by using subscripts q and r and each generates a stress, respectively σ q and σ r so that the total stress is given by
σ t o t = σ q + σ r
Each arm is defined by Equations (27) and (28), with parameters V q , V r , α q , α r , G q , G r as appropriate. The total extension ratio λ is common to both arms. When λ is known, such as when the model is being strained at a constant rate (as in the initial loading of a strain recovery experiment), the solution is straightforward. When the total stress σ t o t is known (as in the recovery phase of the experiment), Equations (27) to (29) are solved iteratively to produce the required stress.
The q arm of the model is associated with the Eyring process and stiffness identified in the analysis of loading behaviour. The operational activation volume V q was as obtained from the slope of the plot in Figure 8. The fitted value of α q from the intercept in Figure 8 was found to give an unrealistically fast strain recovery and has been lowered from 0.00745 to 0.004 s−1. The Gaussian parameter G q is taken from the initial modulus for the amorphous phase of E a = 0.36 GPa. The small–strain behaviour of the Gaussian mechanism ensures that G q = E a / 3 . The r arm has a much lower stiffness and is associated with strain recovery; after the initial strain application, it remains in tension and causes recovery. Its stiffness is estimated from the stress and strain at the end of the stress relaxation experiments and assuming that they are related by a Gaussian model of strength G r . The operational activation volume V r is taken from the stress relaxation fit of Figure 9. Due to the approximation involved in Equation (23) mentioned above, α r cannot be evaluated from the fitted value of c and is varied by trial and error to produce realistic predictions.
The parameters used are summarised in Table 3. The predictions of recovery are given in Figure 11.

4. Discussion

Elastic finite element modelling of the partially crystalline material shows that the crystal phase has a significant effect on its stiffness and allows us to estimate the Young’s modulus of the amorphous phase. The analysis has been extended to elastic–plastic behaviour of the amorphous phase, and a method has been developed to identify its yield stress from stress–strain measurements of the material. This method takes the form of an offset criterion, with the yield stress identified by the difference between the applied stress and the stress extrapolated from the linear elastic region. The overall stress in the model material continues to increase after the amorphous yield stress has been exceeded. In Figure 12, we show the inhomogeneous stress field in the amorphous material for the case in which the stress applied to the model is equal to the amorphous yield stress. In some regions the material has yielded, whereas elsewhere the stress is lower due to shielding of the crystal phase. The method has been used to define the strain–rate dependence of the amorphous yield stress, which is found to be consistent with an Eyring process.
We used the above analysis of yield and some stress relaxation data to help define the properties of a two-Eyring process model. This is clearly a simplification in comparison with the real material, which would be expected to include a spectrum of processes equivalent to many different Eyring parameters. The two-arm model comprises the minimum number of arms required to model recovery, when one arm is in tension and is the mechanism for recovery, while the other one resists the recovery and is in compression. The model parameters have in two cases been adjusted to give realistic recovery behaviour. Comparison of Figure 5 and Figure 11 shows that the level of recovery is approximately as observed, and, in common with the observations, increases as the strain rate of the initial stretch is increased. The effect of the initial strain rate is in practice greater than that predicted. The origin of the effect of strain rate is that, during loading, the stiff (q) arm is extended and the time available for the Eyring component to extend is less for higher strain rates. It therefore extends less and a greater proportion of the strain in the arm q is in the elastic element; there is therefore more potential for the model to recover. The other arm r does not contribute to the level of recovery as it resists recovery, being in compression. The dependence of the strain recovery on the initial strain rate could be increased by lowering the stiffness of the elastic part of the q arm, but this quantity is fixed by the observed material stiffness during loading. Introducing a third arm with stiffness lower than the q arm could remedy this problem; however, the extra parameters that would need to be specified would not be directly measurable and the model would have to be fitted by trial and error. Thus, it appears that this deficiency is a consequence of the model’s simplicity.
It should be noted that we have in effect made a material model of the amorphous phase only. The finite element modelling has the effect of isolating the yield behaviour of the amorphous phase to define the q arm, and the r arm is driven by entropic stress that is only a feature of the amorphous phase. This should be an adequate approach for modelling strain recovery behaviour, while more development would be needed to model levels of overall stress.

Author Contributions

Conceptualization, J.S., P.S., K.N. and P.C.; methodology, J.S.; software, P.S. and J.S.; validation, J.S. and P.S.; formal analysis, J.S.; investigation, J.S., P.S. and K.N.; resources, P.C.; data curation, J.S.; writing—original draft preparation, J.S.; writing—review and editing, J.S., P.S., K.N. and P.C.; visualization, J.S. and P.S.; supervision, J.S. and P.C.; project administration, P.C.; funding acquisition, P.C.

Funding

This research was funded by the Engineering and Physical Sciences Research Council, grant number EP/L020572/1.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Avérous, L. Biodegradable Multiphase Systems Based on Plasticized Starch: A Review. J. Macromol. Sci. Polym. Rev. 2004, 44, 231–274. [Google Scholar] [CrossRef]
  2. Castro-Aguirre, E.; Iñiguez-Franco, F.; Samsudin, H.; Fang, X.; Auras, R. Poly(lactic acid)—Mass production, processing, industrial applications, and end of life. Adv. Drug Delivery Rev. 2016, 107, 333–366. [Google Scholar] [CrossRef] [PubMed]
  3. Athanasiou, K.A.; Niederauer, G.G.; Agrawal, C.M. Sterilization, toxicity, biocompatibility and clinical applications of polylactic acid/polyglycolic acid copolymers. Biomaterials 1996, 17, 93–102. [Google Scholar] [CrossRef]
  4. Brkaric, M.; Baker, K.C.; Israel, R.; Harding, T.; Montgomery, D.M.; Herkowit, H.N. Early failure of bioabsorbable anterial cerivcal fusion plates. J. Spinal Disord. Tech. 2007, 20, 248–254. [Google Scholar] [CrossRef] [PubMed]
  5. Smit, T.H.; Engels, T.A.P.; Wuisman, P.I.J.M.; Govaert, L.E. Time-dependent mechanical strength of 70/30 poly(l,dl-lactide): Shedding light on the premature failure of degradable spinal cages. Spine 2008, 33, 14–18. [Google Scholar] [CrossRef]
  6. Rebelo, R.; Fernandes, M.; Fangueiro, R. Biopolymers in medical implants: A brief review. Procedia Eng. 2017, 200, 236–243. [Google Scholar] [CrossRef]
  7. Xu, J.; Song, J. Polylactic acid (PLA)-based shape-memory materials for biomedical applications. In Shape Memory Polymers for Biomedical Applications; Woodhead Publishing: Cambridge, UK, 2015; pp. 197–217. [Google Scholar]
  8. Farah, S.; Anderson, D.G.; Langer, R. Physical and mechanical properties of PLA, and their functions in widespread applications—A comprehensive review. Adv. Drug Delivery Rev. 2016, 107, 367–392. [Google Scholar] [CrossRef]
  9. Garlotta, D.A. Literature review of Poly(Lactic Acid). J. Polym. Environ. 2001, 9, 63–84. [Google Scholar] [CrossRef]
  10. Bergström, J.S. Hayman D An Overview of Mechanical Properties and Material Modeling of Polylactide (PLA) for Medical Applications. Ann. Biomed. Eng. 2016, 44, 330–340. [Google Scholar] [CrossRef]
  11. Al-Itry, R.; Lamnawar, K.; Maazouz, A.; Billon, N.; Combeaud, C. Effect of the simultaneous biaxial stretching on the structural and mechanical properties of PLA, PBAT and their blends at rubbery state. Eur. Polym. J. 2015, 68, 288–301. [Google Scholar] [CrossRef]
  12. Blair, R.W.; Dunne, N.J.; Lennon, A.B.; Menary, G.H. Processing-property relationships of biaxially stretched poly(l-lactic acid) sheet for application in coronary stents. J. Mech. Behav. Biomed. Mater. 2018, 86, 113–121. [Google Scholar] [CrossRef] [PubMed]
  13. Wang, X.; Pan, Y.; Liu, X.; Liu, H.; Li, N.; Liu, C.; Schubert, D.W.; Shen, C. Facile Fabrication of Superhydrophobic and Eco-Friendly Poly(lactic acid) Foam for Oil−Water Separation via Skin Peeling. ACS Appl. Mater. Interfaces 2019, 11, 14362–14367. [Google Scholar] [CrossRef] [PubMed]
  14. Chacón, J.M.; Caminero, M.A.; García-Plaza, M.; Núnez, P.J. Additive manufacturing of PLA structures using fused deposition modelling: Effect of process parameters on mechanical properties and their optimal selection. Mater. Des. 2017, 124, 143–157. [Google Scholar] [CrossRef]
  15. Coates, P.D.; Ward, I.M. Drawing of Polymers through a Conical Die. Polymer 1979, 20, 1553–1560. [Google Scholar] [CrossRef]
  16. Coates, P.D.; Ward, I.M. Die Drawing—Solid-Phase Drawing of Polymers through a Converging Die. Polym. Eng. Sci. 1981, 21, 612–618. [Google Scholar]
  17. Ward, I.M.; Taraiya, A.K.; Coates, P.D. Solid state extrusion and die drawing. In Solid Phase Processing of Polymers; Hanser: Munich, Germany, 2000; pp. 328–365. [Google Scholar]
  18. Li, J.; Li, Z.; Ye, L.; Zhao, X.; Coates, P.D.; Caton-Rose, F. Structure evolution and orientation mechanism of long-chain-branched poly(lactic acid) in the process of solid die drawing. Eur. Polym. J. 2017, 90, 54–65. [Google Scholar] [CrossRef]
  19. Li, J.; Li, Z.; Ye, L.; Zhao, X.; Coates, P.D.; Caton-Rose, F. Structure and biocompatibility improvement mechanism of highly oriented poly(lactic acid) produced by solid die drawing. Eur. Polym. J. 2017, 90, 68–76. [Google Scholar] [CrossRef]
  20. Vasile, C.; Stoleru, E.; Darie-Niţa, R.N.; Dumitriu, R.P.; Pamfil, D.; Tarţau, R. Biocompatible Materials Based on Plasticized Poly(lactic acid), Chitosan and Rosemary Ethanolic Extract I. Effect of Chitosan on the Properties of Plasticized Poly(lactic acid) Materials. Polymers 2019, 11, 941. [Google Scholar] [CrossRef]
  21. Li, Q.; Liu, C.; Wena, J.; Wua, Y.; Shana, Y.; Liao, J. The design, mechanism and biomedical application of self-healing hydrogels. Chin. Chem. Lett. 2017, 28, 1857–1874. [Google Scholar] [CrossRef]
  22. Ye, A.; Wang, S.; Zhao, Q.; Wang, Y.; Liu, C.; Shen, C. Poly(ethylene oxide)-promoted dispersion of graphene nanoplatelets and its effect on the properties of poly(lactic acid)/poly(butylene adipate-co-terephthalate)-based nanocomposites. Mater. Lett. 2019, 253, 34–37. [Google Scholar] [CrossRef]
  23. Peponi, L.; Navarro-Baena, I.; Kenny, J.M. Shape memory polymers: Properties, synthesis and applications. In Smart Polymers and Their Applications; Aguilar, M.R., Román, J.S., Eds.; Woodhead Publishing in Materials, Elsevier: Cambridge, UK, 2014; pp. 204–236. [Google Scholar]
  24. Fillon, T.M.; Xu, J.; Prasad, M.L.; Song, J. In Vivo Tissue Responses to Thermal-responsive Shape Memory Polymer Nanocomposites. Biomaterials 2011, 32, 985–991. [Google Scholar] [CrossRef]
  25. Walzac, J.; Sobota, M.; Chrzanowski, M.; Krucinska, I. Application of the melt-blown technique in the production of shape-memory nonwoven fabrics from a blend of poly(l-lactide) and atactic poly[(R,S)-3-hydroxy butyrate]. Text. Res. J. 2018, 88, 2141–2152. [Google Scholar]
  26. Zhang, H.; Wang, S.; Zhang, S.; Ma, R.; Wang, Y.; Cao, W.; Liu, C.; Shen, C. Crystallization behavior of poly(lactic acid) with a self-assembly aryl amide nucleating agent probed by real-time infrared spectroscopy and X-ray diffraction. Polym. Test. 2017, 64, 12–19. [Google Scholar] [CrossRef]
  27. Peltoniemi, H.; Kontio, R.; Lindqvist, C.; Suuronen, R. The use of bioabsorbable osteofixation devices in craniomaxillofacial surgery. Oral Surg. Oral Med. Oral Pathol. 2002, 84, 5–14. [Google Scholar] [CrossRef]
  28. Halpin, J.C.; Kardos, J.L. Moduli of crystalline polymers employing composite theory. J. Appl. Phys. 1972, 43, 2235–2241. [Google Scholar] [CrossRef]
  29. Halpin, J.C.; Kardos, J.L. The Halpin-Tsai equations: A review. Polym. Eng. Sci. 1976, 16, 344–352. [Google Scholar]
  30. Kong, Y.; Hay, J.N. The measurement of crystallinity of polymers with DSC. Polymer 2002, 3873–3878. [Google Scholar] [CrossRef]
  31. Schindelin, J.; Arganda-Carreras, I.; Frise, E. Fiji: An open-source platform for biological-image analysis. Nat. Methods 2012, 9, 676–682. [Google Scholar] [CrossRef]
  32. Mori, T.; Tanaka, K. Average stress in matrix and average elastic energy of materials with misfitting inclusions. Acta Metall. 1973, 21, 571–574. [Google Scholar] [CrossRef]
  33. Tucker, C.L.; Liang, E. Stiffness predictions for unidirectional short-fiber composites: Review and evaluation. Compos. Sci. Technol. 1999, 59, 655–671. [Google Scholar] [CrossRef]
  34. Spencer, M.W.; Paul, D.R. Modeling the mechanical and thermal expansion behavior of TPO-based nanocomposites. Polymer 2011, 52, 4910–4919. [Google Scholar] [CrossRef]
  35. Mura, T. Micromechanics of Defects in Solids, 2nd ed.; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1991; pp. 77–84. [Google Scholar]
  36. Lee, S.; Kimoto, M.; Tanaka, M.; Tsuji, H.; Nishino, T. Crystal modulus of poly(lactic acid)s, and their stereocomplex. Polymer 2018, 138, 124–131. [Google Scholar] [CrossRef]
  37. Tábi, T.; Sajó, I.E.; Szabó, F.; Luyt, A.S.; Kovács, J.G. Crystalline structure of annealed polylactic acid and its relation to processing. Express Polym. Lett. 2010, 10, 659–668. [Google Scholar] [CrossRef]
  38. Press, W.H.; Flannery, B.P.; Teukolsky, S.A.; Vetterling, W.T. Numerical Recipes, the Art of Scientific Computing. Cambridge University Press: London, UK, 1986. [Google Scholar]
  39. Söntjens, S.H.M.; Engels, T.A.P.; Smit, T.H.; Govaert, L.E. Time-dependent failure of amorphous poly-d,l-lactide: Influence of molecular weight. J. Mech. Behav. Biomed. Mater. 2012, 13, 69–77. [Google Scholar] [CrossRef]
  40. Zhou, C.; Guo, H.; Li, J.; Huang, S.; Li, H.; Meng, Y.; Yu, D.; de Claville Christiansen, J.; Jiang, S. Temperature dependence of poly(lactic acid) mechanical properties. RSC Adv. 2016, 6, 113762–113772. [Google Scholar] [CrossRef]
  41. Guedes, R.M.; Singh, A.; Pinto, V. Viscoelastic Modelling of Creep and Stress Relaxation Behaviour in PLA-PCL Fibres. Fibers Polym. 2017, 18, 2443–2453. [Google Scholar] [CrossRef]
  42. Dusunceli, N.; Drozdov, A.D.; Theilgaard, N. Influence of Temperature on Viscoelastic–Viscoplastic Behavior of Poly(lactic acid) Under Loading–Unloading. Polym. Eng. Sci. 2017, 57, 239–247. [Google Scholar] [CrossRef]
  43. Sweeney, J.; Bonner, M.; Ward, I.M. Modelling of loading, stress relaxation and stress recovery of a shape memory polymer. J. Mech. Behav. Biomed. Mater. 2014, 37, 12–23. [Google Scholar] [CrossRef]
  44. Tcharkhtchi, A.; Abdallah-Elhirsti, S.; Ebrahimi, K.; Fitoussi, J.; Shirinbayan, M.; Farzaneh, S. Some New Concepts of Shape Memory Effect of Polymers. Polymers 2014, 6, 1144–1163. [Google Scholar] [CrossRef] [Green Version]
  45. Ward, I.M.; Sweeney, J. Mechanical Properties of Solid Polymers, 3rd ed.; Wiley: Chichester, UK, 2013. [Google Scholar]
  46. Halsey, G.; White, H.J.; Eyring, H. Mechanical properties of textiles, I. Text. Res. J. 1945, 15, 295–311. [Google Scholar] [CrossRef]
  47. Guiu, F.; Pratt, P.L. Stress relaxation and the plastic deformation of solids. Phys. Status Solidi. 1964, 6, 111–120. [Google Scholar] [CrossRef]
Figure 1. Mori–Tanaka Modulus predictions for (a) Ecrystal = 14 GPa and (b) Ecrystal = 3 GPa.
Figure 1. Mori–Tanaka Modulus predictions for (a) Ecrystal = 14 GPa and (b) Ecrystal = 3 GPa.
Polymers 11 01342 g001
Figure 2. Orientation of model crystallite with respect to global 1, 2, 3 axes coincident with the RVE cube axes. The a crystallite axes are all parallel to the 3 axis and the each crystallite rotates about the a axis in the 1–2 plane to the angle θ. The c axis is the (stiff) chain axis. The crystallite is assumed to be transversely isotropic in the ab plane. Left: view along the 3 axes. Right: general view.
Figure 2. Orientation of model crystallite with respect to global 1, 2, 3 axes coincident with the RVE cube axes. The a crystallite axes are all parallel to the 3 axis and the each crystallite rotates about the a axis in the 1–2 plane to the angle θ. The c axis is the (stiff) chain axis. The crystallite is assumed to be transversely isotropic in the ab plane. Left: view along the 3 axes. Right: general view.
Polymers 11 01342 g002
Figure 3. Stress–strain curves for the range of strain rates shown in s−1 in the caption.
Figure 3. Stress–strain curves for the range of strain rates shown in s−1 in the caption.
Polymers 11 01342 g003
Figure 4. Stress relaxation for a range of initial strain rates shown in s−1 in the caption.
Figure 4. Stress relaxation for a range of initial strain rates shown in s−1 in the caption.
Polymers 11 01342 g004
Figure 5. Strain recovery under load for a range of strain initial rates shown in s−1 in the captions, recovering at a load corresponding to an engineering stress of 0.89 MPa.
Figure 5. Strain recovery under load for a range of strain initial rates shown in s−1 in the captions, recovering at a load corresponding to an engineering stress of 0.89 MPa.
Polymers 11 01342 g005
Figure 6. RVE realisation containing 43 inclusions with d/w = 0.17 and volume fraction φ = 0.13. (a) Model view without inclusions (b) showing inclusions only.
Figure 6. RVE realisation containing 43 inclusions with d/w = 0.17 and volume fraction φ = 0.13. (a) Model view without inclusions (b) showing inclusions only.
Polymers 11 01342 g006
Figure 7. Comparison of RVE stress with linear extrapolation and stress-strain behaviour of the amorphous phase for yield stresses (a) 1.0; (b) 2.5; (c) 5.0; and (d) 10.0 MPa. Aspect ratios of inclusions are d/w = 0.086 except for (b) where aspect ratios of 0.086 and 0.17 have been compared.
Figure 7. Comparison of RVE stress with linear extrapolation and stress-strain behaviour of the amorphous phase for yield stresses (a) 1.0; (b) 2.5; (c) 5.0; and (d) 10.0 MPa. Aspect ratios of inclusions are d/w = 0.086 except for (b) where aspect ratios of 0.086 and 0.17 have been compared.
Polymers 11 01342 g007
Figure 8. σ Y derived from the offset condition p = 0.138 as a function of strain rate.
Figure 8. σ Y derived from the offset condition p = 0.138 as a function of strain rate.
Polymers 11 01342 g008
Figure 9. Stress relaxation for loading at 0.04 s−1 fitted to Equation (20). (a) Complete curve (b) Initial 30 s.
Figure 9. Stress relaxation for loading at 0.04 s−1 fitted to Equation (20). (a) Complete curve (b) Initial 30 s.
Polymers 11 01342 g009
Figure 10. Two-process model.
Figure 10. Two-process model.
Polymers 11 01342 g010
Figure 11. Strain recovery predictions of model in Figure 10 using parameters in Table 3.
Figure 11. Strain recovery predictions of model in Figure 10 using parameters in Table 3.
Polymers 11 01342 g011
Figure 12. Finite element model of material with 89 inclusions at aspect ratio d/w = 0.086. Material yield stress is 2.5 MPa. Stress levels in MPa.
Figure 12. Finite element model of material with 89 inclusions at aspect ratio d/w = 0.086. Material yield stress is 2.5 MPa. Stress levels in MPa.
Polymers 11 01342 g012
Table 1. Average values of E1 and E2 calculated from RVEs. Errors are standard deviations.
Table 1. Average values of E1 and E2 calculated from RVEs. Errors are standard deviations.
Aspect Ratio d/wN target Number of InclusionsNumber of Realisations s bcbc   ( GPa ) 1 Mean Modulus (GPa)
0.08687110.6670.534 ± 0.017
0.08687121.00.523 ± 0.022
0.1743210.6670.512 ± 0.035
0.3422210.6670.507 ± 0.029
Table 2. Values of offset factor p as a function of yield stress.
Table 2. Values of offset factor p as a function of yield stress.
σ Y   ( MPa ) p
1.00.140
2.50.134
5.00.139
10.00.156
Table 3. Values of parameters for the model of Figure 10.
Table 3. Values of parameters for the model of Figure 10.
ArmV (MPa)−1α (s−1)G (MPa)
q0.450.004120
r2.21.0 × 10−60.74

Share and Cite

MDPI and ACS Style

Sweeney, J.; Spencer, P.; Nair, K.; Coates, P. Modelling the Mechanical and Strain Recovery Behaviour of Partially Crystalline PLA. Polymers 2019, 11, 1342. https://doi.org/10.3390/polym11081342

AMA Style

Sweeney J, Spencer P, Nair K, Coates P. Modelling the Mechanical and Strain Recovery Behaviour of Partially Crystalline PLA. Polymers. 2019; 11(8):1342. https://doi.org/10.3390/polym11081342

Chicago/Turabian Style

Sweeney, John, Paul Spencer, Karthik Nair, and Phil Coates. 2019. "Modelling the Mechanical and Strain Recovery Behaviour of Partially Crystalline PLA" Polymers 11, no. 8: 1342. https://doi.org/10.3390/polym11081342

APA Style

Sweeney, J., Spencer, P., Nair, K., & Coates, P. (2019). Modelling the Mechanical and Strain Recovery Behaviour of Partially Crystalline PLA. Polymers, 11(8), 1342. https://doi.org/10.3390/polym11081342

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