skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Methodological Aspects and Mesh Convergence in Numerical Analysis of Athermal Fiber Network Material Deformation
A balance between model complexity, accuracy, and computational cost is a central concern in numerical simulations. In particular, for stochastic fiber networks, the non-affine deformation of fibers, related non-linear geometric features due to large global deformation, and size effects can significantly affect the accuracy of the computer experiment outputs and increase the computational cost. In this work, we systematically investigate methodological aspects of fiber network simulations with a focus on the output accuracy and computational cost in models with cellular (Voronoi) and fibrous (Mikado) network architecture. We study both p and h-refinement of the discretizations in finite element solution procedure, with uniform and length-based adaptive h-refinement strategies. The analysis is conducted for linear elastic and viscoelastic constitutive behavior of the fibers, as well as for networks with initially straight and crimped fibers. With relative error as the determining criterion, we provide recommendations for mesh refinement, comment on the necessity of multiple realizations, and give an overview of associated computational cost that will serve as guidance toward minimizing the computational cost while maintaining a desired level of solution accuracy.  more » « less
Award ID(s):
2022489
PAR ID:
10509311
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
MDPI
Date Published:
Journal Name:
Fibers
Volume:
12
Issue:
1
ISSN:
2079-6439
Page Range / eLocation ID:
9
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Fiber networks are the primary structural components of many biological structures, including the cell cytoskeleton and the extracellular matrix. These materials exhibit global nonlinearities, such as stiffening in extension and shear, during which the fibers bend and align with the direction of applied loading. Precise details of deformations at the scale of the fibers during strain stiffening are still lacking, however, as prior work has studied fiber alignment primarily from a qualitative perspective, which leaves incomplete the understanding of how the local microstructural evolution leads to the global mechanical behavior. To fill this gap, we studied how axial forces are transmitted inside the fiber network along paths called force chains, which continuously evolve during the course of deformation. We performed numerical simulations on two-dimensional networks of random fibers under uniaxial extension and shear, modeling the fibers using beam elements in finite element software. To quantify the force chains, we identified all chains of connected fibers for which the axial force was larger than a preset threshold and computed the total length of all such chains. To study the evolution of force chains during loading, we computed the derivative of the total length of all force chains with respect to the applied engineering strain. Results showed that the highest rate of evolution of force chains coincided with the global critical strain for strain stiffening of the fiber network. Therefore, force chains are an important factor connecting understanding of the local kinematics and force transmission to the macroscale stiffness of the fiber network. 
    more » « less
  2. We present a study of the mechanical behavior of planar fibrous mats stabilized by inter-fiber adhesion. Fibers of various degrees of tortuosity and of infinite and finite length are considered in separate models. Fibers are randomly distributed, are not cross-linked, and interact through adhesion and friction. The variation of structural parameters such as the mat thickness and the mean segment length between contacts along given fibers with the strength of adhesion is determined. These systems are largely dissipative in that most of the work performed during deformation is dissipated frictionally and only a small fraction is stored as strain energy. The response of the mats to tensile loading has three regimes: a short elastic regime in which no sliding at contacts is observed, a well-defined sliding regime characterized by strain hardening, and a rapid stiffening regime at larger strains. The third regime is due to the formation of stress paths after the fiber tortuosity is pulled out and is absent in mats of finite length fibers. Networks of finite length fibers lose stability during the second regime of deformation. The scaling of the yield stress, which characterizes the transition between the first and the second regimes, and of the second regime's strain hardening modulus, with system parameters such as the strength of adhesion and friction and the degree of fiber tortuosity are determined. The strength of mats of finite length fibers is also determined as a function of network parameters. These results are expected to become useful in the design of electrospun mats and other planar fibrous non-cross-linked networks. 
    more » « less
  3. We study the effect of inter-fiber adhesion on the mechanical behavior of cross-linked ran- dom fiber networks in two dimensions. To this end, we consider networks with connectiv- ity number, z , below, at, and above the isostaticity limit of the structure without adhesion, z c . Fibers store energy in the axial and bending deformation mode and the cross-links are of freely rotating type. Adhesive forces lead to fiber bundling and to a reduction of the total volume of the network. The degree of shrinkage is determined as a function of the strength of adhesion and network parameters. The mechanical response of these struc- tures is further studied in uniaxial tension and compression. The stress-strain curves of networks without inter-fiber adhesion exhibit an initial linear regime, followed by strain stiffening in tension and strain softening and strain localization in compression. In pres- ence of adhesion, the response becomes more complex. The initial linear regime persists, with the effective modulus decreasing and increasing with increasing adhesion in cases with z > z c and z < z c , respectively. The strain range of the linear regime increases signif- icantly with increasing adhesion. Networks with z > z c subjected to tension strain-stiffen at rates that depend on the adhesion strength, but eventually enter a large strain/stress regime in which the response is independent of this parameter. Networks with z < z c are stabilized by adhesion in the unloaded state. Beyond the initial linear regime their tangent modulus gradually decreases, only to increase again at large strains. Adhesive interactions lead to similar effects in compression. Specifically, in the z > z c case, increasing the adhe- sion strength reduces the linear elastic modulus and significantly increases the range of the linear regime, delaying strain localization. This first investigation of the mechanics of cross-linked random networks with inter-fiber adhesion opens the door to the design of soft materials with novel properties. 
    more » « less
  4. Adhesion plays an important role in the mechanics of nanoscale fibers such as various biological filaments, carbon nanotubes and artificial polymeric nanofibers. In this work we study assemblies of non-crosslinked filaments and characterize their adhesion-driven structural evolution and their final stable structure. The key parameters of the problem are the network density, the fiber length, the bending stiffness of fibers and the strength of adhesion. The system of fibers self-organizes in one of three types of structures: locked networks, in which fibers remain in the as-deposited state, cellular networks, in which fibers form bundles and these organize into a larger scale network, and disintegrated networks, in which the network of bundles becomes disconnected. We determine the parametric space corresponding to each of these structures. Further, we identify a triangular structure of bundles, similar to the Plateau triangle occurring in foams, which stabilizes the network of bundles and study in detail the stabilization mechanism. The analysis provides design guidelines and a physical picture of the stability and structure of random fiber networks with adhesion. 
    more » « less
  5. Adhesion plays an important role in the mechanics of nanoscale fibers such as various biological filaments, carbon nanotubes and artificial polymeric nanofibers. In this work we study assemblies of non-crosslinked filaments and characterize their adhesion-driven structural evolution and their final stable structure. The key parameters of the problem are the network density, the fiber length, the bending stiffness of fibers and the strength of adhesion. The system of fibers self-organizes in one of three types of structures: locked networks, in which fibers remain in the as-deposited state, cellular networks, in which fibers form bundles and these organize into a larger scale network, and disintegrated networks, in which the network of bundles becomes disconnected. We determine the parametric space corresponding to each of these structures. Further, we identify a triangular structure of bundles, similar to the Plateau triangle occurring in foams, which stabilizes the network of bundles and study in detail the stabilization mechanism. The analysis provides design guidelines and a physical picture of the stability and structure of random fiber networks with adhesion. 
    more » « less