Purdue University Graduate School
PhD_Thesis_Orzuri_Rique.pdf (11.56 MB)

Multiscale thermoviscoelastic modeling of composite materials

Download (11.56 MB)
posted on 2021-05-05, 23:47 authored by Orzuri Rique GaraizarOrzuri Rique Garaizar
Polymer matrices present in composite materials are prone to have time-dependent behavior very sensitive to changes in temperature. The modeling of thermoviscoelasticity is fundamental for capturing the performance of anisotropic viscoelastic materials subjected to both mechanical and thermal loads, or for manufacturing simulation of composites. In addition, improved plate/shell and beam models are required to efficiently design and simulate large anisotropic composite structures. Numerical models have been extensively used to capture the linear viscoelasticity in composites, which can be generalized in integral or differential forms. The hereditary integral constitutive form has been adopted by many researchers to be implemented into finite element codes, but its formulation is complex and time consuming as it is function of the time history. The differential formulation provides faster computation times, but its applicability has been limited to capture the behavior of three-dimensional thermoviscoelastic orthotropic materials.

This work extends mechanics of structure genome (MSG) to construct linear thermoviscoelastic solid, plate/shell and beam models for multiscale constitutive modeling of three-dimensional heterogeneous materials made of time and temperature dependent constituents. The formulation derives the transient strain energy based on integral formulation for thermorheologically simple materials subject to finite temperature changes. The reduced time parameter is introduced to relate the time-temperature dependency of the anisotropic material by means of master curves at reference conditions. The thermal expansion creep is treated as inherent material behavior. Exact three-dimensional thermoviscoelastic homogenization solutions are also formulated for laminates modeled as an equivalent, homogeneous, anisotropic solid. The new model is implemented in SwiftComp, a general-purpose multiscale constitutive modeling code based on MSG, to handle real heterogeneous materials with arbitrary microstructures, mesostructures or cross-sectional shapes.

Three-dimensional representative volume element (RVE) analyses and direct numerical simulations using a commercial finite element software are conducted to verify the accuracy of the MSG-based constitutive modeling. Additionally, MSG-based plate/shell results are validated against thin-ply high-strain composites experimental data showing good agreement. Numerical cases with uniform and nonuniform cross-sectional temperature distributions are studied. The results showed that unlike MSG, the RVE method exhibits limitations to properly capture the long-term behavior of effective coefficients of thermal expansion (CTEs) when time-dependent constituent CTEs are considered. The analyses of the homogenized properties also revealed that despite the heterogeneous nature of the composite material, from a multiscale analysis perspective, the temperature dependencies of the effective stiffness and thermal stress properties are governed by the same shift factor as the polymer matrix. This conclusion remains the same for MSG-based solid, plate/shell and beam models with uniform temperature distributions.


Degree Type

  • Doctor of Philosophy


  • Aeronautics and Astronautics

Campus location

  • West Lafayette

Advisor/Supervisor/Committee Chair

Wenbin Yu

Additional Committee Member 2

R. Byron Pipes

Additional Committee Member 3

Weinong Chen

Additional Committee Member 4

Johnathan E. Goodsell