Accurate characterization of the mechanical properties of the human brain at both microscopic and macroscopic length scales is a critical requirement for modeling of traumatic brain injury and brain folding. To date, most experimental studies that employ classical tension/compression/shear tests report the mechanical properties of the brain averaged over both the gray and white matter within the macroscopic regions of interest. As a result, there is a missing correlation between the independent mechanical properties of the microscopic constituent elements and the composite bulk macroscopic mechanical properties of the tissue. This microstructural computational study aims to inversely predict the hyperelastic mechanical properties of the axonal fibers and their surrounding extracellular matrix (ECM) from the bulk tissue's mechanical properties. We develop a representative volume element (RVE) model of the bulk tissue consisting of axonal fibers and ECM with the embedded element technique. A multiobjective optimization technique is implemented to calibrate the model and establish the independent mechanical properties of axonal fibers and ECM based on seven previously reported experimental mechanical tests for bulk white matter tissue from the corpus callosum. The result of the study shows that the discrepancy between the reported values for the elastic behavior of white matter in literature stems from the anisotropy of the tissue at the microscale. The shear modulus of the axonal fiber is seven times larger than the ECM, with axonal fibers that also show greater nonlinearity, contrary to the common assumption that both components exhibit identical nonlinear characteristics.
Statement of significance
The reported mechanical properties of white matter microstructure used in traumatic brain injury or brain mechanics studies vary widely, in some cases by up to two orders of magnitude. Currently, the material parameters of the white matter microstructure are identified by a single loading mode or ultimately two modes of the bulk tissue. The presented material models only define the response of the bulk and homogenized white matter at a macroscopic scale and cannot explicitly capture the connection between the material properties of microstructure and bulk structure. To fill this knowledge gap, our study characterizes the hyperelastic material properties of axonal fibers and ECM using microscale computational modeling and multiobjective optimization. The hyperelastic material properties for axonal fibers and ECM presented in this study are more accurate than previously proposed because they have been optimized using seven or six loading modes of the bulk tissue, which were previously limited to only two of the seven possible loading modes. As such, the predicted values with high accuracy could be used in various computational modeling studies. The systematic characterization of the material properties of the human brain tissue at both macro- and microscales will lead to more accurate computational predictions, which will enable a better understanding of injury criteria, and has a positive impact on the improved development of smart protection systems, and more accurate prediction of brain development and disease progression.
more »
« less
An energy-based study of the embedded element method for explicit dynamics
Abstract The embedded finite element technique provides a unique approach for modeling of fiber-reinforced composites. Meshing fibers as distinct bundles represented by truss elements embedded in a matrix material mesh allows for the assignment of more specific material properties for each component rather than homogenization of all of the properties. However, the implementations of the embedded element technique available in commercial software do not replace the material of the matrix elements with the material of the embedded elements. This causes a redundancy in the volume calculation of the overlapping meshes leading to artificially increased stiffness and mass. This paper investigates the consequences in the energy calculations of an explicit dynamic model due to this redundancy. A method for the correction of the edundancy within a finite element code is suggested which removes extra energy and is shown to be effective at correcting the energy calculations for large amounts of redundant volume.
more »
« less
- Award ID(s):
- 1846059
- NSF-PAR ID:
- 10394209
- Date Published:
- Journal Name:
- Advanced Modeling and Simulation in Engineering Sciences
- Volume:
- 9
- Issue:
- 1
- ISSN:
- 2213-7467
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)more » « less
Abstract Composites can be tailored to specific applications by adjusting process variables. These variables include those related to composition, such as volume fraction of the constituents and those associated with processing methods, methods that can affect composite topology. In the case of particle matrix composites, orientation of the inclusions affects the resulting composite properties, particularly so in instances where the particles can be oriented and arranged into structures. In this work, we study the effects of coupled electric and magnetic field processing with externally applied fields on those structures, and consequently on the resulting material properties that arise. The ability to vary these processing conditions with the goal of generating microstructures that yield target material properties adds an additional level of control to the design of composite material properties. Moreover, while analytical models allow for the prediction of resulting composite properties from constituents and composite topology, these models do not build upward from process variables to make these predictions.
This work couples simulation of the formation of microscale architectures, which result from coupled electric and magnetic field processing of particulate filled polymer matrix composites, with finite element analysis of those structures to provide a direct and explicit linkages between process, structure, and properties. This work demonstrates the utility of these method as a tool for determining composite properties from constituent and processing parameters. Initial particle dynamics simulation incorporating electromagnetic responses between particles and between the particles and the applied fields, including dielectrophoresis, are used to stochastically generate representative volume elements for a given set of process variables. Next, these RVEs are analyzed as periodic structures using FEA yielding bulk material properties. The results are shown to converge for simulation size and discretization, validating the RVE as an appropriate representation of the composite volume. Calculated material properties are compared to traditional effective medium theory models. Simulations allow for mapping of composite properties with respect to not only composition, but also fundamentally from processing simulations that yield varying particle configurations, a step not present in traditional or more modern effective medium theories such as the Halpin Tsai or double-inclusion theories.
-
Taylor & Francis (Ed.)n this research, a method is examined by which the behavior of continuous carbon fiber rein-forced additive manufacturing may be simulated using Finite Element Analysis. This technique is used in a simulated tensile test experiment in which the findings are compared to results determined from theoretical calculations according to the Rule of Mixtures method and from existing mechanical testing results. Four different fiber reinforcement configurations are examined with fiber volume fractions ranging from 4% to 32%. It was found that for fiber vol-ume fractions of 11%, the simulation results closely match those predicted theoretically by the Rule of Mixtures as well as the mechanical testing results published in existing research. Lower fiber volume fractions near 4% yield less accurate results, with a 20% error due to the fact that the anisotropic behavior of the polymer matrix is the dominant material trait. Simulation of higher volume fractions near 32% closely approximate theoretical predictions, however neither the theoretical results nor the simulation results accurately reflect real world mechanical testing, indicating that nonideal condition factors such as the effect of micro-voids between the start and end of the fiber reinforcements play a significant role in the overall strength of the material. Thus, for fiber volume fractions near 11%, this simulation method can accurately be used to predict the behavior of end-use components, but more study must be done to increase simulation accuracy in low and high fiber volume fractions.more » « less
-
Abstract Flexible electronics often employ composite inks consisting of conductive flakes embedded in a polymer matrix to transmit electrical signal. Recently, localized necking was identified as a cause of a substantial increase in normalized resistance with applied strain thereby adversely impacting electrical performance. The current study explores two possible contributing factors for the formation of such localization—ink surface roughness and local variations in silver flake volume fraction. Uniaxial tension experiments of a DuPont 5025 type ink are used to inform a constitutive model implemented using finite element method on different substrates. Surface roughness was modeled by sinusoidal variation in ink height, whose amplitude and wavelength are informed by experimental laser profilometry scan data. Local flake fraction variations obtained from experimental measurements before applying any strain, were modeled as local variations in the elastic modulus according to an inverse rule of mixtures between the silver flake and acrylic binder material properties. The study identified that the ink height roughness is the most impactful contributor to the subsequent strain localization. The substrate elastic properties impact the number and magnitude of localization bands, with the stiffer substrate delocalizing strain and averting catastrophic crack formation seen with a more compliant substrate. The model incorporating surface roughness closely matches experimental measurements of local strain across different substrates. The study can inform designers of the adverse impact of ink surface roughness on localization and subsequent detrimental increase of the resistance.more » « less
-
more » « less
Mechanophores (MPs) undergo chemical reactions to become fluorescent in response to a mechanical stimulus that reflects the magnitude and distribution of applied stress. MPs are an emerging technology for self‐reporting damage sensing applications in polymeric materials in the aeronautical, energy generation, and automotive industries. However, quantitative calibration of the MP response to local stresses remains an outstanding challenge. Herein, a method to calibrate the intensity of the MP fluorescent activation (
I ) with local hydrostatic stresses (σ h ) is presented. Uniaxial tension is applied to a simple composite comprised of a rigid sphere (silica) embedded in a MP‐functionalized elastomeric matrix (spiropyran (SPN) functionalized polydimethylsiloxane (PDMS)). By monitoring the fluorescence intensity with a confocal microscope while a quasi‐static deformation is applied, in situ observations of MP activation as a function of applied uniaxial strain are obtained. To calculate the associated stress fields, a finite element analysis (FEA) with cohesive zone elements is employed. By comparingσ h , calculated through FEA with theI of the PDMS/SPN system, a linear relationship betweenI andσ h is directly determined. The technique presented can be employed for many MP‐containing materials systems to calibrateI toσ h .
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1186/s40323-022-00223-x
