Multiscale computational homogenisation of 3D textile-based fiber reinforced polymer composites

Zahur Ullah, Lukasz Kaczmarczyk, Edward Archer, AT McIlhagger, Eileen Harkin-Jones

Research output: Contribution to conferenceOther

Abstract

This paper presents a multiscale computational homogenisation approach for the calculation of homogenised
structural level mechanical properties of 3D textile/woven based fiber reinforced polymer (FRP) composites. Textile
or woven composites, in which interlaced fibres are used as reinforcement, are a class of FRP composites which
provide flexibility of design and functionality and are used in many engineering applications, including ships,
aircrafts, automobiles, civil structures and prosthetics [1]. The more recently developed 3D-textile composites,
consisting of 3D arrangements of yarns in a polymer matrix, allow weaving of near-net-shape and complex structures
as compared to the traditional 2D-textile composites. In addition, these 3D-textile composites provide high throughthickness
mechanical properties, lower manufacturing cost and improved impact and delamination resistance. The
macro or structural level mechanical properties of these composites are rooted in their underlying complicated and
heterogeneous micro structures. The heterogeneous microstructure of these composites requires a detailed multiscale
computational homogenisation, which results in the macroscopic constitutive behaviour based on their
microscopically heterogeneous representative volume elements (RVE). Elliptical cross sections and cubic splines are
used respectively to model the cross sections and paths of the yarns within these RVEs. The RVE geometry along
with other input parameters, e.g. material properties and boundary conditions, are modelled in CUBIT/Trelis using a
parameterised Python script.

The multiscale computational homogenisation scheme, with a unified imposition of RVE boundary conditions, is
implemented in MoFEM (Mesh Oriented Finite Element Method) [2], which allows convenient switching between
linear displacement, uniform traction and periodic boundary conditions. MoFEM utilises hierarchic basis functions
[3], which permits the use of arbitrary order of approximation leading to accurate results for relatively coarse meshes.
The matrix and yarns within the RVEs are modelled by considering isotropic and transversely isotropic materials
models respectively. The principal direction of the yarns required for the transversely isotropic material model is
calculated using a computationally inexpensive potential flow analysis along these yarns. Furthermore, the
computational framework is designed to take advantage of distributed memory high-performance computing. The
implementation and performance of the computational tool is demonstrated with a variety of 2.5D and 3D woven
based FRP composites including 3D orthogonal interlock, 3D orthogonal layer-to-layer interlock, 3D orthogonal
through-the-thickness angle interlock, 2.5D layer-to-layer angle interlock and 2.5D layer-layer angle interlock [4].


REFERENCES:
[1] Z. Ullah, Ł. Kaczmarczyk, S. A. Grammatikos, M. C. Evernden and C. J. Pearce (2016). Multi-scale
computational homogenisation to predict the long-term durability of composite structures. Computers and Structures, 181 . pp. 21-31.
[3] Ł. Kaczmarczyk, Z. Ullah, K. Lewandowski, X. Meng, X. -Y. Zhou, I. Athanasiadis, I and C. J. Pearce (2017). MoFEM-v0.6.20. Zenodo. http://doi.org/10.5281/zenodo.1053811
[3] M. Ainsworth and J. Coyle (2003). Hierarchic finite element bases on unstructured tetrahedral meshes.
International Journal for Numerical Methods in Engineering, 58 (14): 2103–2130,
[4] Y. Rahali, M. Assidi, I. Goda, A. Zghal, and J. F. Ganghoffer (2016). Computation of the effective mechanical properties including nonclassical moduli of 2.5 D and 3D interlocks by micromechanical approaches. Composites Part B: Engineering, 98, 194-212.

Conference

Conference6th European Conference on Computational Mechanics (ECCM 6)
7th European Conference on Computational Fluid Dynamics (ECFD 7)
Abbreviated titleECCM 6 and ECFD 7
CountryUnited Kingdom
CityGlasgow
Period11/06/1815/06/18
Internet address

Fingerprint

Textiles
Fibers
Composite materials
Polymers
Yarn
Boundary conditions
Finite element method
Mechanical properties
Microstructure
Potential flow
Composite structures
Prosthetics
Polymer matrix
Delamination
Splines
Automobiles
Numerical methods
Materials properties
Reinforcement
Ships

Keywords

  • Fiber reinforced polymer composites
  • 3D textile/woven composites
  • Finite element analysis
  • Multiscale computational homogenization

Cite this

Ullah, Z., Kaczmarczyk, L., Archer, E., McIlhagger, AT., & Harkin-Jones, E. (2018). Multiscale computational homogenisation of 3D textile-based fiber reinforced polymer composites. 6th European Conference on Computational Mechanics (ECCM 6)
7th European Conference on Computational Fluid Dynamics (ECFD 7) , Glasgow, United Kingdom.
Ullah, Zahur ; Kaczmarczyk, Lukasz ; Archer, Edward ; McIlhagger, AT ; Harkin-Jones, Eileen. / Multiscale computational homogenisation of 3D textile-based fiber reinforced polymer composites. 6th European Conference on Computational Mechanics (ECCM 6)
7th European Conference on Computational Fluid Dynamics (ECFD 7) , Glasgow, United Kingdom.
@conference{e33853f2db834c2b9bbebbdd6c147b95,
title = "Multiscale computational homogenisation of 3D textile-based fiber reinforced polymer composites",
abstract = "This paper presents a multiscale computational homogenisation approach for the calculation of homogenisedstructural level mechanical properties of 3D textile/woven based fiber reinforced polymer (FRP) composites. Textileor woven composites, in which interlaced fibres are used as reinforcement, are a class of FRP composites whichprovide flexibility of design and functionality and are used in many engineering applications, including ships,aircrafts, automobiles, civil structures and prosthetics [1]. The more recently developed 3D-textile composites,consisting of 3D arrangements of yarns in a polymer matrix, allow weaving of near-net-shape and complex structuresas compared to the traditional 2D-textile composites. In addition, these 3D-textile composites provide high throughthicknessmechanical properties, lower manufacturing cost and improved impact and delamination resistance. Themacro or structural level mechanical properties of these composites are rooted in their underlying complicated andheterogeneous micro structures. The heterogeneous microstructure of these composites requires a detailed multiscalecomputational homogenisation, which results in the macroscopic constitutive behaviour based on theirmicroscopically heterogeneous representative volume elements (RVE). Elliptical cross sections and cubic splines areused respectively to model the cross sections and paths of the yarns within these RVEs. The RVE geometry alongwith other input parameters, e.g. material properties and boundary conditions, are modelled in CUBIT/Trelis using aparameterised Python script.The multiscale computational homogenisation scheme, with a unified imposition of RVE boundary conditions, isimplemented in MoFEM (Mesh Oriented Finite Element Method) [2], which allows convenient switching betweenlinear displacement, uniform traction and periodic boundary conditions. MoFEM utilises hierarchic basis functions[3], which permits the use of arbitrary order of approximation leading to accurate results for relatively coarse meshes.The matrix and yarns within the RVEs are modelled by considering isotropic and transversely isotropic materialsmodels respectively. The principal direction of the yarns required for the transversely isotropic material model iscalculated using a computationally inexpensive potential flow analysis along these yarns. Furthermore, thecomputational framework is designed to take advantage of distributed memory high-performance computing. Theimplementation and performance of the computational tool is demonstrated with a variety of 2.5D and 3D wovenbased FRP composites including 3D orthogonal interlock, 3D orthogonal layer-to-layer interlock, 3D orthogonalthrough-the-thickness angle interlock, 2.5D layer-to-layer angle interlock and 2.5D layer-layer angle interlock [4].REFERENCES:[1] Z. Ullah, Ł. Kaczmarczyk, S. A. Grammatikos, M. C. Evernden and C. J. Pearce (2016). Multi-scalecomputational homogenisation to predict the long-term durability of composite structures. Computers and Structures, 181 . pp. 21-31.[3] Ł. Kaczmarczyk, Z. Ullah, K. Lewandowski, X. Meng, X. -Y. Zhou, I. Athanasiadis, I and C. J. Pearce (2017). MoFEM-v0.6.20. Zenodo. http://doi.org/10.5281/zenodo.1053811[3] M. Ainsworth and J. Coyle (2003). Hierarchic finite element bases on unstructured tetrahedral meshes.International Journal for Numerical Methods in Engineering, 58 (14): 2103–2130,[4] Y. Rahali, M. Assidi, I. Goda, A. Zghal, and J. F. Ganghoffer (2016). Computation of the effective mechanical properties including nonclassical moduli of 2.5 D and 3D interlocks by micromechanical approaches. Composites Part B: Engineering, 98, 194-212.",
keywords = "Fiber reinforced polymer composites, 3D textile/woven composites, Finite element analysis, Multiscale computational homogenization",
author = "Zahur Ullah and Lukasz Kaczmarczyk and Edward Archer and AT McIlhagger and Eileen Harkin-Jones",
year = "2018",
month = "6",
day = "11",
language = "English",
note = "6th European Conference on Computational Mechanics (ECCM 6)<br/>7th European Conference on Computational Fluid Dynamics (ECFD 7) , ECCM 6 and ECFD 7 ; Conference date: 11-06-2018 Through 15-06-2018",
url = "http://www.eccm-ecfd2018.org/frontal/default.asp",

}

Ullah, Z, Kaczmarczyk, L, Archer, E, McIlhagger, AT & Harkin-Jones, E 2018, 'Multiscale computational homogenisation of 3D textile-based fiber reinforced polymer composites' 6th European Conference on Computational Mechanics (ECCM 6)
7th European Conference on Computational Fluid Dynamics (ECFD 7) , Glasgow, United Kingdom, 11/06/18 - 15/06/18, .

Multiscale computational homogenisation of 3D textile-based fiber reinforced polymer composites. / Ullah, Zahur; Kaczmarczyk, Lukasz; Archer, Edward; McIlhagger, AT; Harkin-Jones, Eileen.

2018. 6th European Conference on Computational Mechanics (ECCM 6)
7th European Conference on Computational Fluid Dynamics (ECFD 7) , Glasgow, United Kingdom.

Research output: Contribution to conferenceOther

TY - CONF

T1 - Multiscale computational homogenisation of 3D textile-based fiber reinforced polymer composites

AU - Ullah, Zahur

AU - Kaczmarczyk, Lukasz

AU - Archer, Edward

AU - McIlhagger, AT

AU - Harkin-Jones, Eileen

PY - 2018/6/11

Y1 - 2018/6/11

N2 - This paper presents a multiscale computational homogenisation approach for the calculation of homogenisedstructural level mechanical properties of 3D textile/woven based fiber reinforced polymer (FRP) composites. Textileor woven composites, in which interlaced fibres are used as reinforcement, are a class of FRP composites whichprovide flexibility of design and functionality and are used in many engineering applications, including ships,aircrafts, automobiles, civil structures and prosthetics [1]. The more recently developed 3D-textile composites,consisting of 3D arrangements of yarns in a polymer matrix, allow weaving of near-net-shape and complex structuresas compared to the traditional 2D-textile composites. In addition, these 3D-textile composites provide high throughthicknessmechanical properties, lower manufacturing cost and improved impact and delamination resistance. Themacro or structural level mechanical properties of these composites are rooted in their underlying complicated andheterogeneous micro structures. The heterogeneous microstructure of these composites requires a detailed multiscalecomputational homogenisation, which results in the macroscopic constitutive behaviour based on theirmicroscopically heterogeneous representative volume elements (RVE). Elliptical cross sections and cubic splines areused respectively to model the cross sections and paths of the yarns within these RVEs. The RVE geometry alongwith other input parameters, e.g. material properties and boundary conditions, are modelled in CUBIT/Trelis using aparameterised Python script.The multiscale computational homogenisation scheme, with a unified imposition of RVE boundary conditions, isimplemented in MoFEM (Mesh Oriented Finite Element Method) [2], which allows convenient switching betweenlinear displacement, uniform traction and periodic boundary conditions. MoFEM utilises hierarchic basis functions[3], which permits the use of arbitrary order of approximation leading to accurate results for relatively coarse meshes.The matrix and yarns within the RVEs are modelled by considering isotropic and transversely isotropic materialsmodels respectively. The principal direction of the yarns required for the transversely isotropic material model iscalculated using a computationally inexpensive potential flow analysis along these yarns. Furthermore, thecomputational framework is designed to take advantage of distributed memory high-performance computing. Theimplementation and performance of the computational tool is demonstrated with a variety of 2.5D and 3D wovenbased FRP composites including 3D orthogonal interlock, 3D orthogonal layer-to-layer interlock, 3D orthogonalthrough-the-thickness angle interlock, 2.5D layer-to-layer angle interlock and 2.5D layer-layer angle interlock [4].REFERENCES:[1] Z. Ullah, Ł. Kaczmarczyk, S. A. Grammatikos, M. C. Evernden and C. J. Pearce (2016). Multi-scalecomputational homogenisation to predict the long-term durability of composite structures. Computers and Structures, 181 . pp. 21-31.[3] Ł. Kaczmarczyk, Z. Ullah, K. Lewandowski, X. Meng, X. -Y. Zhou, I. Athanasiadis, I and C. J. Pearce (2017). MoFEM-v0.6.20. Zenodo. http://doi.org/10.5281/zenodo.1053811[3] M. Ainsworth and J. Coyle (2003). Hierarchic finite element bases on unstructured tetrahedral meshes.International Journal for Numerical Methods in Engineering, 58 (14): 2103–2130,[4] Y. Rahali, M. Assidi, I. Goda, A. Zghal, and J. F. Ganghoffer (2016). Computation of the effective mechanical properties including nonclassical moduli of 2.5 D and 3D interlocks by micromechanical approaches. Composites Part B: Engineering, 98, 194-212.

AB - This paper presents a multiscale computational homogenisation approach for the calculation of homogenisedstructural level mechanical properties of 3D textile/woven based fiber reinforced polymer (FRP) composites. Textileor woven composites, in which interlaced fibres are used as reinforcement, are a class of FRP composites whichprovide flexibility of design and functionality and are used in many engineering applications, including ships,aircrafts, automobiles, civil structures and prosthetics [1]. The more recently developed 3D-textile composites,consisting of 3D arrangements of yarns in a polymer matrix, allow weaving of near-net-shape and complex structuresas compared to the traditional 2D-textile composites. In addition, these 3D-textile composites provide high throughthicknessmechanical properties, lower manufacturing cost and improved impact and delamination resistance. Themacro or structural level mechanical properties of these composites are rooted in their underlying complicated andheterogeneous micro structures. The heterogeneous microstructure of these composites requires a detailed multiscalecomputational homogenisation, which results in the macroscopic constitutive behaviour based on theirmicroscopically heterogeneous representative volume elements (RVE). Elliptical cross sections and cubic splines areused respectively to model the cross sections and paths of the yarns within these RVEs. The RVE geometry alongwith other input parameters, e.g. material properties and boundary conditions, are modelled in CUBIT/Trelis using aparameterised Python script.The multiscale computational homogenisation scheme, with a unified imposition of RVE boundary conditions, isimplemented in MoFEM (Mesh Oriented Finite Element Method) [2], which allows convenient switching betweenlinear displacement, uniform traction and periodic boundary conditions. MoFEM utilises hierarchic basis functions[3], which permits the use of arbitrary order of approximation leading to accurate results for relatively coarse meshes.The matrix and yarns within the RVEs are modelled by considering isotropic and transversely isotropic materialsmodels respectively. The principal direction of the yarns required for the transversely isotropic material model iscalculated using a computationally inexpensive potential flow analysis along these yarns. Furthermore, thecomputational framework is designed to take advantage of distributed memory high-performance computing. Theimplementation and performance of the computational tool is demonstrated with a variety of 2.5D and 3D wovenbased FRP composites including 3D orthogonal interlock, 3D orthogonal layer-to-layer interlock, 3D orthogonalthrough-the-thickness angle interlock, 2.5D layer-to-layer angle interlock and 2.5D layer-layer angle interlock [4].REFERENCES:[1] Z. Ullah, Ł. Kaczmarczyk, S. A. Grammatikos, M. C. Evernden and C. J. Pearce (2016). Multi-scalecomputational homogenisation to predict the long-term durability of composite structures. Computers and Structures, 181 . pp. 21-31.[3] Ł. Kaczmarczyk, Z. Ullah, K. Lewandowski, X. Meng, X. -Y. Zhou, I. Athanasiadis, I and C. J. Pearce (2017). MoFEM-v0.6.20. Zenodo. http://doi.org/10.5281/zenodo.1053811[3] M. Ainsworth and J. Coyle (2003). Hierarchic finite element bases on unstructured tetrahedral meshes.International Journal for Numerical Methods in Engineering, 58 (14): 2103–2130,[4] Y. Rahali, M. Assidi, I. Goda, A. Zghal, and J. F. Ganghoffer (2016). Computation of the effective mechanical properties including nonclassical moduli of 2.5 D and 3D interlocks by micromechanical approaches. Composites Part B: Engineering, 98, 194-212.

KW - Fiber reinforced polymer composites

KW - 3D textile/woven composites

KW - Finite element analysis

KW - Multiscale computational homogenization

M3 - Other

ER -

Ullah Z, Kaczmarczyk L, Archer E, McIlhagger AT, Harkin-Jones E. Multiscale computational homogenisation of 3D textile-based fiber reinforced polymer composites. 2018. 6th European Conference on Computational Mechanics (ECCM 6)
7th European Conference on Computational Fluid Dynamics (ECFD 7) , Glasgow, United Kingdom.