Tissue engineering involves three-dimensional scaffolds to support cell culture activities and provide mechanical support. One of the potential scaffolds used in tissue engineering is an electrospun scaffold consisting fibres ranging from nano- to micrometer scales deposited on layer stack. The finite element models have been used to study the in-plane deformation of two-dimensional single layer fibrous networks and without considers out-of-plane deformation. While the existing study focuses on two-dimensional study, the out-of-plane deformation of layer structured of electrospun scaffolds through the scaffolds thickness has not been studied. In this study, three-dimensional finite element model was constructed to investigate the fracture of multilayer fibrous networks. The three-dimensional results were compared with the fracture on two-dimensional single layer fibrous network. The result shows that these two models had identical fracture behaviour and similar deformation at the crack-tip region, where the fibres are rearranged and reoriented with similar stress distribution. The work here concludes that two-dimensional single layer fibrous network model is a simple yet effective model for the study of homogeneous fibrous networks

Chung, S. M

Theng, Koh Ching

English

Journals of Universiti Tun Hussein Onn Malaysia (UTHM)

Finite Element Modelling for Fracture of Multilayer Fibrous Networks

International Journal of Integrated Engineering

penerbit.uthm.edu.my/ojs/index.php/ijie  Abstract: Tissue engineering involves three-dimensional scaffolds to support cell culture activities and provide mechanical support. One of the potential scaffolds used in tissue engineering is an electrospun scaffold consisting fibres ranging from nano- to micrometer scales deposited on layer stack. The finite element models have been used to study the in-plane deformation of two-dimensional single layer fibrous networks. The two-dimensional models do not consider the out-of-plane deformation of layer structured of electrospun scaffolds through the scaffolds thickness. In this study, three-dimensional finite element model was constructed to investigate the fracture of multilayer fibrous networks. The three-dimensional results were compared with the fracture on two-dimensional single layer fibrous network. The result shows that these two models had identical fracture behaviour and similar deformation at the crack-tip region, where the fibres are rearranged and reoriented with similar stress distribution. The work here concludes that two-dimensional single layer fibrous network model is a simple yet effective model for the study of homogeneous fibrous networks. Keywords: Fibrous networks, multilayer, fracture INTERNATIONAL JOURNAL OF INTEGRATED ENGINEERING VOL. 12 NO. 8 (2020) 117-122    © Universiti Tun Hussein Onn Malaysia Publisher’s Office  IJIE Journal homepage: http://penerbit.uthm.edu.my/ojs/index.php/ijie ISSN : 2229-838X e-ISSN : 2600-7916 International Journal of Integrated Engineering   Finite Element Modelling for Fracture of Multilayer Fibrous Networks  S. M. Chung1, C. T. Koh1* 1Faculty of Mechanical and Manufacturing Engineering,  Universiti Tun Hussein Onn Malaysia, 86400 Parit Raja, Batu Pahat, Johor, MALAYSIA  *Corresponding Author  DOI: https://doi.org/10.30880/ijie.2020.12.08.011 Received 2 May 2019; Accepted 8 June 2020; Available online 30 August 2020    1. Introduction Electrospinning technique is a simple technique to produce fibrous materials with fibers ranging from nano- to micrometer scales that deposited on layer stack [1-2]. The simplicity and diversity of electrospinning technique attract the use of electropun scaffolds in various applications. One of the applications is to act as scaffold in tissue engineered construct. The microstructure architecture of fibrous scaffolds is designed to meet a certain requirement of tissue engineering. Sufficient mechanical property is required in tissue engineering to support loading applied on the tissue engineered constructs [3]. In recent years, the approach to understand the mechanical behavior of fibrous networks has been intensive. These approaches are done by either experimental samples or computational models [4–13]. The current approaches to generate computational model for fibrous networks are based on multiscale simulation combining microscopic and macroscopic models [6], post processing of SEM images [14, 15], fiber deposition method and dynamic analysis of fibrous networks [16, 17]. The priority for a computational model is to develop a simple and effective model to study the mechanical properties of fibrous networks under different loading condition. Hence, a certain assumption is usually found in computational models. One of the assumptions is treating a fibrous network as a two-dimensional single layer fibrous network without considering out-of-plane deformation [18]. The two-dimensional models are often modelled with large thickness and assume that the networks deform in in-plane condition. This leaves an interesting question on the effect of this projected thickness on two-dimensional model. Further, it remains unclear on how significant the out-of-plane   *Corresponding author: ctkoh@uthm.edu.my 2020 UTHM Publisher. All rights reserved. 117 118  S. M. Chung et al., Int. J. of Integrated Engineering Vol. 12 No. 8 (2020) p. 117-122   deformation affects the mechanics of thin scaffolds. The study on out-of-plane deformation require a three-dimensional multilayer fibrous networks model that consider other parameters like torsional stiffness of fibers. The objective of this work is to compare the fracture behavior of two-dimensional (2D) single layer fibrous network and three-dimensional (3D) multilayer fibrous networks. Both 2D and 3D consist identical network architecture of an electrospun scaffold. Further, a study on the effect of projected thickness on 2D single layer fibrous network was constructed. Such findings suggest a simple and effective computational model for the study of homogeneous fibrous networks.  2. Methodology 2.1 Modelling of Fibrous Networks MATLAB code from previous work [16] was used to construct 2D single layer fibrous networks. Fibrous networks were generated in MATLAB 2017b software by placing random line at random positions with random angles. These random lines were then extended to model boundary and were randomly bonded at the same node depending on user- defined cross-link percentage (i.e. the ratio between the sum of bonding points and the sum of fiber intersection points). The fiber density (i.e. the sum of fibers per unit area) and cross-link density (i.e. the sum of bonding points per unit area) of fibrous networks were calculated. In this work, 2D fibrous networks were constructed with fiber density of 6403.6 mm-1 and were cross-linked at cross-link percentage of 15 % and 50 %. The fiber density is defined as the sum of fiber length per unit area. The fibers were modelled with rectangular profile with the width of 300 nm and thicknesses of 1 m [18], 1 mm, 10 μm and 1 μm. The thickness was defined as model thickness in this paper. The 3D multilayer fibrous networks were generated by using MATLAB 2017b software. The networks were constructed layer by layer with 2D model construction method. Each layer was constructed by sequence and connected to each other. The networks were cross-linked within layer by intralayer cross-link and bonded between layer to layer by interlayer cross-link. Fig. 1 shows multilayer fibrous networks that consists of 3 layers of fibrous networks, which are bonded between layers at the thickness of a. In this study, 3D multilayer fibrous networks were constructed in 3 and 5 layers and with constant fiber density of 6558.01 ± 124.03 mm-1 for each layer. The networks were cross-linked within layer and between layers at cross-link percentage of 15 %. The total cross-link density of multilayer fibrous networks are made of the sum of intralayer and interlayer bonding points for the first network layer. The fibres were modelled with circular profile with diameter of 300 nm.  Fig. 1 - Schematic illustration of the finite element model of a multilayer fibrous network  All the networks models were imported into finite element software, ABAQUS 2017 and were analyzed by using nonlinear finite element analysis, which considers large strain and rotation. The fibers were defined with Young’s modulus, E of 100 MPa and fracture strength, σf of 30 ± 4 MPa [18].  2.2 Fracture Analysis The 2D and 3D fibrous networks were constructed in a circular unit with 75 μm radius and with a 75 μm notch length. The boundary conditions were subjected to the displacement field associated with the macroscopic crack tip field for a homogeneous and isotropic materials. The boundary conditions were calculated from the equations obtained in Kanninen and Popelar [19]:  u = 1 r K1 κ − 1 + 2sin2 θ cos θ     (1) 1 2 2π G  2  2 119  1 r   K1 2 2π  G S. M. Chung et al., Int. J. of Integrated Engineering Vol. 12 No. 8 (2020) p. 117-122   u  = κ + 1 − 2 cos2 θ sin θ    (2) 2  2  2  Where (u1, u2) are the displacements for x and y directions at the model boundary, r is the distance of the node from the crack-tip, θ is the angle between node and crack-tip and K1 is the stress intensity factor for mode I. The fibrous network models were expected to have plane stress condition, where κ = (3 − v) (1+ v) with Poisson’s ratio, v of 0.3. The shear modulus of fibrous network, G was assumed to be 4 MPa [16]. The fracture toughness was characterized by measuring stress intensity factor (SIF) once failure criterion met. The network was assumed to fracture once the first fibre exceeded its fracture strength.  3. Results 3.1 Fracture of 2D Single Layer Fibrous Network The 2D fibrous networks were constructed to study the effect of model thickness on fracture behavior of fibrous networks. Fig. 2 shows the stress intensity factors for 2D single layer fibrous networks with various model thicknesses. The result shows that the networks have consistent stress intensity factor for model thickness ranging from 1 µm to 1 m. Moreover, the deformation of networks at the crack-tip is similar for all the networks. The critical failure location of these networks are found at the same fiber and location. The result indicates that 2D fibrous network with model thickness ranging from 1 µm to 1 m has similar fracture behavior. This indicates that the 1 µm thickness is large enough to allow the networks deform in in-plane manner.  Fig. 2 - Stress intensity factor of 2D single layer fibrous networks constructed with various model thicknesses  3.2 Fracture of 3D Multilayer Fibrous Networks The 3D multilayer fibrous networks were constructed to study the effect of out-of-plane deformation on thin scaffolds. Fig. 3 shows the comparison of the stress intensity factor for 2D single layer fibrous networks modelled with model thickness of 1 m and 3D multilayer fibrous networks modelled with circular cross-section. The result shows that 3D multilayer fibrous networks constructed with 3 and 5 layers have comparable stress intensity factor to 2D fibrous network with cross-link percentage of 50 % instead of 15 %. This is attributed to the additional interlayer cross-link from 3D fibrous networks that act towards each network layer. The total cross-link density of multilayer fibrous networks were now close to the cross-link density of 2D fibrous networks that constructed with cross-link percentage of 50 %. Further, the deformation of these networks at the crack-tip is found to be similar, where the fibres are rearranged and reoriented with similar stress distribution (Fig. 4). 120  S. M. Chung et al., Int. J. of Integrated Engineering Vol. 12 No. 8 (2020) p. 117-122    Fig. 3 - Comparison of the stress intensity factor for 2D single layer fibrous networks and 3D multilayer fibrous networks   Fig. 4 - The stress distribution for (a) 2D single layer fibrous network with cross-link percentage of 50 % and (b-d) 3D multilayer fibrous networks 121  S. M. Chung et al., Int. J. of Integrated Engineering Vol. 12 No. 8 (2020) p. 117-122     Fig. 5 - The out-of-plane deformation for (a) first layer, (b) second layer and (c) third layer of 3D multilayer fibrous networks  Fig. 5 shows the out-of-plane deformation for each network layer of 3D multilayer fibrous networks. The finite element result shows that 3D fibrous networks deform along the thickness of networks, where the first and the last layer of networks deformed towards the mid-plane of networks (i.e. second layer of networks). The out-of-plane deformation of each network layer are relatively small as compared to the thickness of networks. The network layers deform with a maximum out-of-plane displacement of 80 nm as compared to the networks thickness of 600 nm.  4. Discussion Multilayer fibrous networks are often simplified into 2D fibrous networks that consider only in-plane deformation. The out-of-plane deformation of multilayer fibrous networks considers other parameters like torsional stiffness and bending stiffness and could possibly affect the fracture of fibrous networks. In this work, 2D and 3D models having homogeneous fibrous networks are shown to have identical stress intensity factor. To explain these results, analysis on the out-of-plane deformation, deformation mechanism and stress distribution of multilayer fibrous networks were considered. The relatively small deformation of each network layer along the thickness indicates that multilayer fibrous networks deform mainly in in-plane condition with minimum influence from out-of-plane deformation. The deformation mechanism of multilayer fibrous networks at the crack-tip region are similar to 2D fibrous network, which involves of fiber rearrangement and reorientation. This shows that the toughening mechanism of both models are similar. The 2D studies had highlighted the importance of fiber rearrangement and reorientation on the fracture of fibrous networks [16, 20, 21]. Further, the stress distribution of multilayer fibrous networks at the crack-tip region are identical to 2D fibrous network. Both 2D and 3D models had highest stress at the crack-tip and the stress are progressively distributed to the far-field. The result suggests that 3D multilayer fibrous networks can be reduced to 2D single layer fibrous network without any significant impact in simulating the fracture of homogeneous fibrous networks. The assumption of in plane deformation however relies on the thickness and consistent distribution of fibrous network structures. The homogeneous fibrous networks investigated here had uniform random fibrous networks throughout the thickness. The 2D models require further validation for the fibrous materials having non-uniform microstructures [24-25]. In general, finite element analysis is an effective tool to study mechanical behaviors  of medical devices under different loading conditions [26].  5. Conclusion The computational models for fibrous networks are often treated as a 2D single layer fibrous network with projected thickness. This work demonstrated that such model is an appropriate fracture model, where the fracture behaviour and deformation at the crack-tip region are similar to 3D multilayer fibrous networks. In addition, the work here suggests that 2D single layer fibrous networks modelled with projected thickness ranging from 1 µm to 1 m predict similar fracture behaviour. The work here demonstrates the suitability of 2D finite element models in simulating the fracture of homogeneous fibrous networks.  Acknowledgement The authors would like to express their gratitude to University of Tun Hussein Onn Malaysia [GPPS/U807], Ministry of Science, Technology and Innovation Malaysia [03-01-13-SF0112-S026] and the Ministry of Education Malaysia [FRGS/1/2014/TK01/UTHM/02/2-1462] for funding the research work. 122  S. M. Chung et al., Int. J. of Integrated Engineering Vol. 12 No. 8 (2020) p. 117-122   References [1] W. Khoo, C.T. Koh, S.C. Lim (2017). Synthetic and Natural Fibrous Scaffolds for Soft Tissue Engineering Applications. Journal of Mechanical Engineering 4(3), 223-233. [2] S.C. Lim, W. Khoo, C.T. Koh. (2017). Electrospinning of Polycaprolactone (PCL) and Gelatin Polymeric Fibers. Journal of Mechanical Engineering 4(3), 246-254. [3] P. Lu and B. Ding, (2008), Applications of electrospun fibers, Recent Pat. Nanotechnol. 2(3), pp. 169–82. [4] W.J. Li, C.T. Laurencin, E.J. Caterson, R.S. Tuan, and F.K. Ko (2002), Electrospun nanofibrous structure: a novel scaffold for tissue engineering, J. Biomed. Mater. Res. 60 (4), pp. 613–621. [5] D.W. Hutmacher (2000), Scaffolds in tissue engineering bone and cartilage, Biomaterials 21, pp. 2529–2543. [6] Y. Yin and J. Xiong (2018), Finite element analysis of electrospun nanofibrous mats under biaxial tension, Nanomaterials 8(5), p. 348. [7] Y. Yin, Z. Pan, and J. Xiong (2018) A tensile constitutive relationship and a finite element model of electrospun nanofibrous mats, Nanomaterials 8(1), p. 29. [8] T. Stylianopoulos, C. A. Bashur, A. S. Goldstein, S. A. Guelcher, and V. H. Barocas (2008), Computational predictions of the tensile properties of electrospun fibre meshes: effect of fibre diameter and fibre orientation, J. Mech. Behav. Biomed. Mater. (4), pp. 326–335. [9] F. Croisier, A. S. Duwez, C. Jérôme, A. F. Léonard, K. O. Van der Werf, P. J. Dijkstra, and M. L. Bennink (2012), Mechanical testing of electrospun PCL fibers, Acta Biomaterlia 8(1), pp. 218–224. [10] C. Xiang and M. W. Frey (2016), Increasing mechanical properties of 2-D-structured electrospun nylon 6 non- woven fiber mats, Materials (Basel) 9(4). [11] X. Mao (2018), Mechanics of disordered fiber networks, Gels and Other Soft Amorphous Solids, pp. 199–210. [12] S. C. Wong, A. Baji, and S. Leng (2008), Effect of fiber diameter on tensile properties of electrospun poly(ε- caprolactone), Polymer (Guildf) 49 (21), pp. 4713–4722. [13] C. T. Koh, C. Y. Low, and Y. Yusri (2015), Structure-property relationship of bio-inspired fibrous materials, Procedia Comput. Sci. 76, pp. 411–416. [14] K. Mubyana, R. A. Koppes, K. L. Lee, J. A. Cooper, and D. T. Corr (2016), The influence of specimen thickness and alignment on the material and failure properties of electrospun polycaprolactone nanofiber mats, J. Biomed. Mater. Res. A, pp. 1–7. [15] T. P. Driscoll, N. L. Nerurkar, N. T. Jacobs, D. M. Elliott, and R. L. Mauck (2011), Fiber angle and aspect ratio influence the shear mechanics of oriented electrospun nanofibrous scaffolds, J. Mech. Behav. Biomed. Mater 4(8), pp. 1627–1636. [16] J. B. Carleton, A. D’Amore, K. R. Feaver, G. J. Rodin, and M. S. Sacks (2015), Geometric characterization and simulation of planar layered elastomeric fibrous biomaterials, Acta Biomater. 12(1), pp. 93–101. [17] A. D’Amore, N. Amoroso, R. Gottardi, C. Hobson, C. Carruthers, S. Watkins, W. R. Wagner, and M. S. Sacks (2014), From single fiber to macro-level mechanics: a structural finite-element model for elastomeric fibrous biomaterials, J. Mech. Behav. Biomed. Mater. 39, pp. 146–161. [18] C. T. Koh and M. L. Oyen (2012), Branching toughens fibrous networks, J. Mech. Behav. Biomed. Mater., vol. 12, pp. 74–82. [19] C. T. Koh and M. L. Oyen (2015), Toughening in electrospun fibrous scaffolds, APL Mater., vol. 3, p. 14908. [20] T. van Dillen, P. R. Onck, and E. Van der Giessen (2008), Models for stiffening in cross-linked biopolymer networks: a comparative study, J. Mech. Phys. Solids, 56 (6), pp. 2240–2264. [21] M. F. Kanninen and C. H. Popelar (1985), Advanced fracture mechanics. [22] U. Stachewicz, I. Peker, W. Tu, and A. H. Barber (2011), Stress delocalization in crack tolerant electrospun nanofiber networks, ACS Appl. Mater. Interfaces 3 (6), pp. 1991–1996. [23] Theng, K. C. and Chung, S. (2020) “Effects of Microstructure Architecture on the Fracture of Fibrous Materials”, International Journal of Integrated Engineering, 12(1), pp. 138-143. [24] W. Khoo, S.M. Chung, S.C. Lim, C. Y. Low, J. M. Shapiro, C. T. Koh (2019) “Fracture behavior of multilayer fibrous scaffolds featuring microstructural gradients”, Materials & Design 184, 108184. [25] W. Khoo, S. M. Chung, S. C. Lim, C. Y. Low, J. M. Shapiro, C. T. Koh (2019) “Dataset on microstructural characteristics and mechanical performance of homogeneous and functionally graded fibrous scaffolds”, Data in Brief 27, 104718. [26] A. Kareem, O. Fakhri, A. Ismail, I. Taib, I. Ali, M.N. Nemah, (2020). Comparative Computational Study of Mechanical Behavior in Self-Expanding Femoropopliteal Stents. International Journal of Integrated Engineering, 12(1), 210-224 

https://publisher.uthm.edu.my/ojs/index.php/ijie/article/download/4283/3749

Finite Element Modelling for Fracture of Multilayer Fibrous Networks

Abstract

Similar works

Full text

Available Versions

Journals of Universiti Tun Hussein Onn Malaysia (UTHM)

International Journal of Integrated Engineering