0
Research Papers

Modeling of the Long-Term Behavior of Glassy Polymers

[+] Author and Article Information
Mathias Wallin

Division of Solid Mechanics,
Lund University,
P. O. Box 118,
SE-221 00 Lund, Sweden

Contributed by the Materials Division of ASME for publication in the Journal of Engineering Materials and Technology. Manuscript received March 23, 2012; final manuscript received August 17, 2012; published online December 21, 2012. Assoc. Editor: Yoshihiro Tomita.

J. Eng. Mater. Technol 135(1), 011001 (Dec 21, 2012) (11 pages) Paper No: MATS-12-1055; doi: 10.1115/1.4007499 History: Received March 23, 2012; Revised August 17, 2012

The constitutive model for glassy polymers proposed by Arruda and Boyce (BPA model) is reviewed and compared to experimental data for long-term loading. The BPA model has previously been shown to capture monotonic loading accurately, but for unloading and long-term behavior, the response of the BPA model is found to deviate from experimental data. In the present paper, we suggest an efficient extension that significantly improves the predictive capability of the BPA model during unloading and long-term recovery. The new, extended BPA model (EBPA model) is calibrated to experimental data of polycarbonate (PC) in various loading–unloading situations and deformation states. The numerical treatment of the BPA model associated with the finite element analysis is also discussed. As a consequence of the anisotropic hardening, the plastic spin enters the model. In order to handle the plastic spin in a finite element formulation, an algorithmic plastic spin is introduced. In conjunction with the backward Euler integration scheme use of the algorithmic plastic spin leads to a set of algebraic equations that provides the updated state. Numerical examples reveal that the proposed numerical algorithm is robust and well suited for finite element simulations.

© 2013 by ASME
Your Session has timed out. Please sign back in to continue.

References

Figures

Grahic Jump Location
Fig. 1

Rheological illustration of the BPA model. The model is governed by the elements: (a) elastic spring, (b) viscoplastic dashpot, and (c) nonlinear Langevin spring.

Grahic Jump Location
Fig. 2

The idealized chain structure, according to the 8-chain model. The dimension of the 8-chain cube is a0 and φ0 denotes the initial orientation angle between a chain and the principal axes. The unit vectors eα, α = 1,2,3, align with the principal directions of V¯p.

Grahic Jump Location
Fig. 4

The stress–strain response of the BPA model during uniaxial compression for repeated unloadings to π = 1.2 MPa. The experimental results are taken from Ref. [10].

Grahic Jump Location
Fig. 3

The stress–strain response of the BPA model during uniaxial compression. The experimental results are taken from Ref. [10].

Grahic Jump Location
Fig. 5

Rheological illustration of the EBPA model. The model is governed by the elements: (a) elastic spring, (b) two viscoplastic dashpots, and (c) nonlinear Langevin spring.

Grahic Jump Location
Fig. 6

Stress–strain curves for uniaxial compression of PC. Experimental data is taken from Ref. [7].

Grahic Jump Location
Fig. 7

Stress–strain curves for plane strain compression of PC. Experimental data is taken from Ref. [7].

Grahic Jump Location
Fig. 8

Stress–strain curves for simple shear of PC. Experimental data is taken from Ref. [5].

Grahic Jump Location
Fig. 9

True stress versus true strain for uniaxial compression of bisphenol A polycarbonate. Experimental data is taken from Ref. [10].

Grahic Jump Location
Fig. 13

The force–displacement response of the cylindrical bar. The deformation rate, u·/L = 0.001 1/s, was used in the simulation. s0 is a constitutive parameter given in Table 3.

Grahic Jump Location
Fig. 14

(a) The geometry and finite element mesh of the cylindrical bar. The deformed mesh and the plastic stretch, λ¯p, at (b) u/L = 0.08, (c) u/L = 0.10, and (d) u/L = 0.18.

Grahic Jump Location
Fig. 10

Athermal shear strength, s, and its components s1 and s2 versus true strain, according to the EBPA model. The backstresses, βBPAdev and βEBPAdev, in the direction of the applied load.

Grahic Jump Location
Fig. 11

Uniaxial compression responses for bisphenol A polycarbonate, according to (a) the BPA and (b) the EBPA model. The repeated unloadings are performed to π = 59 MPa. Experimental data is taken from Ref. [10].

Grahic Jump Location
Fig. 12

Uniaxial compression responses for bisphenol A polycarbonate, according to (a) the BPA and (b) the EBPA model. The repeated unloadings are performed to π = 1.2 MPa. Experimental data is taken from Ref. [10].

Grahic Jump Location
Fig. 15

(a) The geometry of the specimen subjected to combined tension and bending. Visualization of the plastic stretch, λ¯p, in the deformed meshes at the end of (b) loading (by EBPA), (c) unloading (by EBPA), and (d) dwell period of two months (by both models).

Grahic Jump Location
Fig. 16

The displacement component, ux versus time during unloading and the dwell period of two months. Solid and dashed line present the EBPA and BPA response, respectively. The markers, ▪ and ♦, indicate the final position.

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In