Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “polymer theory”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Molecular Dynamics Simulation and Theoretical Analysis of Structural Relaxation, Bond Exchange Dynamics, and Glass Transition in Vitrimers

Vitrimers are a class of polymer networks featuring dynamic covalent cross-links that can undergo associative bond exchange. There has been recent interest in these materials due to their promise as recyclable thermosets or self-healing polymers because of the ability of vitrimer networks to rearrange at the molecular level and undergo macroscopic flow. However, the practical use of these materials often occurs in the supercooled regime or glassy state, where the implications of dynamic bonds are complicated by the interplay between slow activated segmental dynamics, cross-link (i.e., bond-exchange) kinetics, and ultimately material properties. Here, in this paper, we combine coarse-grained molecular dynamics simulation and microscopic statistical mechanical theory to understand how cross-linking kinetics affect material dynamics and how this couples to segmental relaxation of the polymeric network strands across a spectrum of length and time scales, especially in the supercooled regime. We characterize the Kuhn segmental alpha relaxation time and bond exchange time for vitrimer systems across various cross-link densities, temperatures, and bond exchange rates. Simulation and theory both exhibit a bending-up behavior for bond exchange time upon cooling, suggesting a coupling between bond exchange dynamics and segmental relaxation that intensifies with faster bond exchange kinetics. We also found bond exchange dynamics have an impact on Kuhn segment alpha relaxation time, which is most significant at higher cross-link densities. Both these effects are most prominent when the bond exchange time is similar to the Kuhn segment alpha relaxation time, and the resulting coupling of these two relaxation processes is tied to both the probability of a free end to find a bonded pair and the time scale of the constraints imposed by the dynamic cross-links. This relationship is reflected by a cross-link dependence of a theoretical parameter which represents the quantitative degree of coupling between bond exchange and segmental dynamics. Overall, the combination of simulation and theory clarifies the intricate interaction between bond kinetics and segmental relaxation and demonstrates the ability to provide molecular-level insights into vitrimer dynamics over a wide temperature range.

dynamic relaxation↗
DOE OSTI2021

Self-Assembly of Flexible Linear–Semiflexible Bottlebrush–Flexible Linear Triblock Copolymers

Block copolymer (BCP) self-assembly is a fundamental process in which incompatible blocks spontaneously form organized microstructures with broad practical applications. Classical understanding is that the domain spacing is limited by the contour length of the polymer backbone. Here, using a combination of molecular design, chemical synthesis, small-/wide-angle X-ray scattering, transmission electron microscopy, and electron tomography, we discover that this molecular picture does not hold for architecturally semiflexible BCPs. For strongly segregated linear–semiflexible bottlebrush–linear triblock copolymers, the size of the bottlebrush domain can be twice the bottlebrush backbone contour length. The mechanism of such anomalous self-assembly is likely that the interfacial repulsion between the incompatible blocks is large enough to pull a part of the linear end blocks into the bottlebrush domain. This effectively increases the bottlebrush domain size. Moreover, the semiflexible bottlebrush widens the regime for the cylinder morphology that is associated with the volume fraction of the end blocks fC SFB ∈ (0.10, >0.41). This window is much wider than that for flexible linear BCPs, fC F ∈ (0.14, 0.35), and that predicted by the recent self-consistent field theory for linear-bottlebrush BCPs of the same chemistry and molecular architecture. Finally, our experimental findings reveal previously unrecognized mechanisms for the self-assembly of architecturally complex BCPs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Applying Improved Optical Recognition with Machine Learning on Sorting Cu Impurities in Steel Scrap

Annually, 20–55 million tons of electronic waste (e-waste) is produced worldwide (5% of all municipal solid waste). Although e-waste embodies only 2% of America’s municipal waste, it accounts for a significantly larger proportion of the heavy metals and flame retardants present in the waste stream. Currently, less than 20% of all e-waste is recycled in the United States because the heterogeneity of the feedstock limits the opportunity for reuse in high value products and processing of the waste itself is often too costly to justify handling. To address this concern, the aims of this study are: 1) to identify the major plastic and metal compositions within electronic shredder residue (ESR), 2) to formulate solvents and processing conditions to separate 90% of the plastics targeted from consumer shred ESR, and 3) to develop a process design model to estimate the cost and energy efficiency of the proposed solvent-based processing. In this study, a pre-sorted heterogeneous ESR feedstock (one where aluminum, magnetic components, and hazardous battery materials removed by an e-waste recycling facility) was used, with the major compositions of the ESR characterized. It was found that 25 wt.% of the feedstock was composed of plastics, 6 wt.% rubber, 27 wt.% printed circuit boards, 23% wire, and the remainder metals and capacitors. Within the plastic portion, polystyrene (PS, 40 wt.%), acrylonitrile butadiene styrene (ABS, 25 wt.%), and styrene-acrylonitrile (SAN, 9wt.%) were identified to compose the majority of the screened plastics using Fourier transform infrared spectroscopy (FTIR). Next, selective solvents were screened using Hansen Solubility Parameter Theory (HSP) for dissolving PS and ABS. The pre-screening results show that methylene chloride (dichloromethane, DCM) and tetrahydrofuran (THF) are capable of dissolving the most PS and ABS, while methanol (MeOH) and ethylene glycol (EG) are capable of precipitating the most PS and ABS. These solvents were subsequently used to recover polymers and remove. flame retardants within the ESR feedstock. By optimizing the dissolution time and the solvents used, the highest polymer dissolution yield (99 wt.%%) was achieved using DCM for 48 hr. Both pre-screened anti-solvents (MeOH and EG) showed the highest polymer precipitation yield (71 wt.%). In terms of flame retardant removal rate, EG was found to have a high phosphorus-containing flame retardant removal rate (up to 98%). Characterization shows that the proposed solvent-based processing can preserve a high molecular weight fraction of the polymers and effectively remove flame retardants. Cost analysis indicates that the amount of the solvent/anti-solvent recovered after the reaction would play a critical role in reducing the operating costs. The energy analysis shows that the proposed solvent-based processes can save up to 60% of the embodied energy used to manufacture plastics used in electronics (PS and ABS were used for calculations). The results from this project prove the potential of solvent-based processing to produce secondary materials (plastics and metals) from e-waste for cross-industry reuse.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Theoretical Pathway toward Improved Reverse Osmosis Membrane Selectivity for Neutral Solutes: Inspiration from Gas Separations

We report highly rigid membrane materials with tailored structures have exhibited permeabilities and selectivities that exceed the polymer upper bound in gaseous and organic solvent separations for challenging mixtures containing species that are similar in size and shape. One potential question is whether such membrane materials can maintain meaningful guest diffusivities in situations where the microporous spaces are essentially full of guest molecules. Here, we use a simplified transition state theory approach to estimate the diffusivity of water and small organics within a microporous membrane. The transition state theory model is parameterized using experimental values from zeolites and carbon molecular sieve materials found in the literature. We demonstrate the differences in transport and Maxwell–Stefan diffusivities based on guest species loading with different isotherm behaviors. These calculations theorize a path forward for highly selective reverse osmosis membranes for aqueous phase separations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

A description of structured waves in shock compressed particulate composites

Dynamic compression of composite materials is of scientific interest because the mechanical mismatch between internal phases challenges continuum theories. Typical assumptions about steady wave propagation and quasi-instantaneous state changes require reexamination along with the need for time-dependent models. To that end, data and models are presented here for the shock compression of an idealized particulate composite. To serve as a generic representative of this material class, a polymer matrix was filled with tungsten particles, ranging from 1 to 50 vol. %. This creates a simple microstructure containing randomly scattered particles with an extreme impedance mismatch to the binding matrix. These materials were parallel plate impact loaded by Al flyers traveling at 1.8–5.0 km/s. Velocimetry provided records of the equilibrium state and the compression wave structure for each case with trends quantified by an empirical fit. The same quantities were also studied as a function of the wave's propagation distance. A homogenized viscoelastic model then made it possible to progress from cataloging trends to making predictions. Starting from a Mie–Gruneisen equation of state, additional time varying terms were added to capture the transient response. After calibration, accurate predictions of the steady wave structure were possible.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Temperature dependence of nuclear quantum effects on liquid water via artificial neural network model based on SCAN meta-GGA functional

We investigate the temperature dependence of nuclear quantum effects (NQEs) on structural and dynamic properties of liquid water by training a neural network force field using first-principles molecular dynamics (FPMD) based on the strongly constrained and appropriately normed meta-generalized gradient approximation exchange-correlation approximation. The FPMD simulation based on density functional theory has become a powerful computational approach for studying a wide range of condensed phase systems. However, its large computational cost makes it difficult to incorporate NQEs in the simulation and investigate temperature dependence of various properties. To circumvent this difficulty, we use an artificial neural network model and employ the thermostatted ring polymer MD approach for studying the temperature dependence of NQEs on various properties. The NQEs generally bring the radial distribution functions closer to the experimental measurements. Translational diffusivity and rotational dynamics of water molecules are both slowed down by the NQEs. The competing inter-molecular and intra-molecular quantum effects on hydrogen bonds, as discussed by Habershon, Markland, and Manolopoulos [J. Chem. Phys. 131(2), 024501 (2019)], can explain the observed temperature dependence of the NQEs on the dynamical properties in our simulation.

Yao, Yi (ORCID:0000000299576979)↗
DOE OSTI2021

Chemical Recycling of Mixed Plastics and Valuable Metals in the Electronic Waste Using Solvent-Based Processing

Annually, 20-55 million tons of electronic waste (e-waste) is produced worldwide (5% of all municipal solid waste). Although e-waste embodies only 2% of America’s municipal waste, it accounts for a significantly larger proportion of the heavy metals and flame retardants present in the waste stream. Currently, less than 20% of all e-waste is recycled in the United States because the heterogeneity of the feedstock limits the opportunity for reuse in high value products and processing of the waste itself is often too costly to justify handling. To address this concern, the aims of this study are: 1) to identify the major plastic and metal compositions within electronic shredder residue (ESR), 2) to formulate solvents and processing conditions to separate 90% of the plastics targeted from consumer shred ESR, and 3) to develop a process design model to estimate the cost and energy efficiency of the proposed solvent-based processing. In this study, a pre-sorted heterogeneous ESR feedstock (one where aluminum, magnetic components, and hazardous battery materials removed by an e-waste recycling facility) was used, with the major compositions of the ESR characterized. It was found that 25 wt.% of the feedstock was composed of plastics, 6 wt.% rubber, 27 wt.% printed circuit boards, 23% wire, and the remainder metals and capacitors. Within the plastic portion, polystyrene (PS, 40 wt.%), acrylonitrile butadiene styrene (ABS, 25 wt.%), and styrene-acrylonitrile (SAN, 9wt.%) were identified to compose the majority of the screened plastics using Fourier transform infrared spectroscopy (FTIR). Next, selective solvents were screened using Hansen Solubility Parameter Theory (HSP) for dissolving PS and ABS. The pre-screening results show that methylene chloride (dichloromethane, DCM) and tetrahydrofuran (THF) are capable of dissolving the most PS and ABS, while methanol (MeOH) and ethylene glycol (EG) are capable of precipitating the most PS and ABS. These solvents were subsequently used to recover polymers and remove. flame retardants within the ESR feedstock. By optimizing the dissolution time and the solvents used, the highest polymer dissolution yield (99 wt.%%) was achieved using DCM for 48 hr. Both pre-screened anti-solvents (MeOH and EG) showed the highest polymer precipitation yield (71 wt.%). In terms of flame retardant removal rate, EG was found to have a high phosphorus-containing flame retardant removal rate (up to 98%). Characterization shows that the proposed solvent-based processing can preserve a high molecular weight fraction of the polymers and effectively remove flame retardants. Cost analysis indicates that the amount of the solvent/anti-solvent recovered after the reaction would play a critical role in reducing the operating costs. The energy analysis shows that the proposed solvent-based processes can save up to 60% of the embodied energy used to manufacture plastics used in electronics (PS and ABS were used for calculations). The results from this project prove the potential of solvent-based processing to produce secondary materials (plastics and metals) from e-waste for cross-industry reuse.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Rational sub-nanometer manipulation of polymer morphology for efficient chemical separations

The ability to control polymer morphology on the sub-nanometer length scale has broad implications for chemical separations. To achieve such control on easily processable systems, this proposal focuses on the synthesis and characterization of polymers containing appended labile moieties that are easily detached by thermolysis or UV irradiation deep within the glassy state. Once liberated, these moieties can diffuse from the polymer matrix as gaseous products, leaving behind templated pathways for selective diffusion and sorption of small molecules. With a specific target of creating polymeric membrane materials with unprecedented diffusion and (ad)sorption characteristics for chemical separations, synthesis of new materials will be complemented with advanced metrologies, simulations, and evaluation of thermodynamic and transport theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Bending and shear improvements in 3D-printed core sandwich composites through modification of resin uptake in the skin/core interphase region

Fused deposition modeling (FDM) printed polymers are rarely used as a structural material due to anisotropic and low mechanical properties compared with conventional composites. In recent years, greater need has been expressed for recycling of materials, such as recyclable FDM, at the end of service life to reduce environmental pollution and manufacture cost. However, how the amount of resin uptake in the skin and skin/core interphase affects the bending and shear performance of the sandwich composites when replacing the low strength and ductile core (conventional core) with a high strength and brittle core (FDM printed PLA (polylactic acid) core) still remains unclear. A new manufacturing routine is needed to improve the incorporation of FDM printed polymers in composite structures. In this work, FDM printed PLA was used as core material and sandwiched between two unidirectional glass fiber reinforced polymer (GFRP) skins to form a sandwich composite by compression-molding (CM) process, which provides a good manufacturing strategy for skin/core interphase modification. The significance of the CM process is proved by investigating the effect of resin uptake on bending and in-plane/out-of-plane shear performances. Current first order shear deformation (FSDT) theory lacks a direct connection between the in-plane shear stress and out-of-shear stress in the core region of sandwich composites. With the help of DIC, a connection between the in-plane shear and the out-of-plane shear strain was built and in-plane shear properties can acquire through out-of-plane shear properties, hence reducing the redundancy of sample preparation or the need for simulation. A significant improvement was found compared with the optimized resin uptake (Optimized resin uptake range: 20.43%–22.86 wt%) 3D-printed PLA core sandwich composite and lowest performance sandwich composite (Improvement: in-plane shear strength (~34%)/modulus (~29%), out-of-plane shear strength (~25%)/modulus (~31%), specific peak bending load (~19%)). Finally, compared with balsa core sandwich composites, the 3D-printed cores are suitable for use in composite sandwich structures in many applications with a satisfactory strength-to-weight ratio.

3D-printed core sandwich composite↗
DOE OSTI2025

Backbone Stitching in Bottlebrush Copolymer Mesodomains and the Impact of Side Chain Crystallization

We synthesized bottlebrush statistical copolymers (BSCPs) having poly­(ethylene oxide) (PEO) and poly­(dimethylsiloxane) (PDMS) side chains attached to a polynorbornene backbone. Small-angle X-ray scattering analysis showed that for densely grafted BSCPs, the scattering length density gradually transitions between the PEO and PDMS domains. For loosely grafted BSCPs, the polymer backbone formed a distinct mesodomain, with a lower electron and mass density than both the PEO and PDMS domains. The bottlebrush backbone essentially “stitches” the PEO and PDMS side chains, looping back and forth from the PEO to PDMS domains with the backbone segments oriented normal to the domain interfaces. Self-consistent field theory (SCFT) calculations validated the stitching of the backbone driven by the microphase separation of PEO and PDMS, along with a strong segmental order of the side chains in the melt. The reduced birefringence upon PEO crystallization suggests the disruption of the strong segmental order by the crystallization. Both the static intrinsic and the form birefringences of the BSCPs decreased upon PEO crystallization. Solid-state NMR confirmed the rigidity of PEO crystallites and the bottlebrush backbone. Self-assembly of BSCPs containing polyhedral oligomeric silsesquioxane (POSS) pendent groups was also evaluated by X-ray scattering, showing the formation of lamellar microdomains that inhibited POSS crystallization.

Hu, Mingqiu↗
DOE OSTI2025

Photoinduced electronic excitations drive polymerization of carbon monoxide: A first-principles study

Under pressure, carbon monoxide (CO) transforms into a polymer that can be recovered to ambient conditions. While this transformation can occur without additional stimuli, experimental observations have shown that laser irradiation can induce a similar transformation at reduced pressure. The resulting polymeric phase, which is metastable under ambient conditions, releases energy through decomposition into more stable configurations. Using time-dependent density functional theory and Born–Oppenheimer molecular dynamics simulations, we investigate the mechanism by which electronic excitation facilitates CO polymerization. Our calculations reveal that electronic excitation enhances carbon–carbon bonding, enabling polymerization at pressures significantly lower than those required by conventional compression methods. In conclusion, these findings suggest that a photo-assisted approach could be employed to synthesize novel, potentially energetic materials under less demanding pressure conditions.

Born-Oppenheimer molecular dynamics↗
DOE OSTI2022

Path Integrals for Nonadiabatic Dynamics: Multistate Ring Polymer Molecular Dynamics

This review focuses on a recent class of path-integral-based methods for the simulation of nonadiabatic dynamics in the condensed phase using only classical molecular dynamics trajectories in an extended phase space. Specifically, a semiclassical mapping protocol is used to derive an exact, continuous, Cartesian variable path-integral representation for the canonical partition function of a system in which multiple electronic states are coupled to nuclear degrees of freedom. Building on this exact statistical foundation, multistate ring polymer molecular dynamics methods are developed for the approximate calculation of real-time thermal correlation functions. As a result, the remarkable promise of these multistate ring polymer methods, their successful applications, and their limitations are discussed in detail.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1980

Extensional flow of dilute polymer solution

Simultaneous effects of temperature and concentration on tensile viscosity are investigated using a semi-empirical scaling law established to relate the concentration, temperature and deformation rate by a reduced deformation rate variable. Results show the onset of high viscosity occurring in the vicinity of the onset of non-Newtonian viscosity observed in a simple shear flow experiment. High viscosity is also compared with steady and dynamic shear data. Finally, the asymptotic solution of Hassager's kinetic theory of a bead-rod model for macromolecular solutions is compared with the data, indicating that tensile viscosity depends strongly on the molecular weight, the higher the molecular weight, the more sensitive the tensile viscosity in terms of the deformation rate.

Peng, S. T. J.↗
DOE OSTI2025

Machine learning-accelerated path integral molecular dynamics simulations of reactive organic electrolytes

Hydrogen bonded electrolytes that exhibit accelerated proton transport via sequential reactive hops have drawn interest for their promise in clean energy applications. Molecular dynamics simulations of these electrolytes offer the opportunity to uncover microscopic mechanistic details that could be used to design and tune the properties of candidate electrolyte technologies. However, accurately modeling the proton transfer reactions and transport properties that give rise to high charge conductivites in these electrolytes proves computationally challenging because of the need to perform lengthy condensed phase simulations, treating both the electronic and nuclear degrees of freedom quantum mechanically. In this paper, we demonstrate that such a modeling task can be efficiently achieved with the use of density functional theory (DFT)-trained machine learning potentials (MLP) to accelerate path integral molecular dynamics (PIMD) simulations. We highlight the practical utility of this approach by using it to benchmark how closely PIMD simulations employing different DFT exchange–correlation functionals reproduce the composition-dependent densities, diffusion coefficients, and electrical conductivities of mixtures consisting of imidazole and levulinic acid. Even with the speedup afforded by our MLPs, PIMD simulations remain quite expensive. Furthermore, in order to render PIMD more computationally tractable, we introduce and benchmark the accuracy of a ring polymer contraction approach that leverages a computationally efficient short-range MLP to accelerate our PIMD simulations by an additional factor of four.

Chemical bonding↗
DOE OSTI2026

Supramolecular Arene‐Perfluoroarene Assembly Enhances Photoiniferter Polymerization Kinetics

Ultra-high-molecular-weight (UHMW, >106 Da) polymers have unique thermomechanical properties due to their large number of physical entanglements. Producing UHMW styrenic polymers is challenging because of their prohibitively slow homopolymerization kinetics. In this report, we show that alternating photoiniferter copolymerization between methoxy-functionalized styrenic monomers and pentafluorostyrene can circumvent this limitation under conditions where arene–perfluoroarene (AP) interactions contribute to favorable monomer association and accelerate copolymerization. Increasing methoxy substitution on the styrene arene ring leads to pronounced rate enhancements that correlate with electronic complementarity between the styrenic comonomers. Density functional theory calculations reveal increasingly favorable AP interaction energies across the series, consistent with monomer association contributing to faster propagation. The result is up to 2000% enhancement in propagation rate in copolymerizations compared to styrene homopolymerization. Despite these rate enhancements, the resulting materials retain broadly similar glass-transition temperatures and network-like thermomechanical behavior, with only modest softening across the series. Here, we then leverage the approach to produce UHMW trialkoxystyrene copolymers bearing sterically encumbered pendants, which are ultrasoft (Young's modulus of 8 kPa) and highly extensible (560%). These results establish that designed supramolecular association between styrenics is a powerful tool for controlling reactivity in radical polymerization and generating otherwise difficult-to-obtain materials.

Marquez, Joshua D. [University of Florida, Gainesv↗
NASA NTRS1988

Demixing of aqueous polymer two-phase systems in low gravity

Experiments on the effect of gravity on phase emulsion demixing were performed in order to study the potential of two-phase partitioning in space. The behavior of phase systems possessing phase density differences and systems with essentially identical phase densities was investigated during one-g and low-g parabolic aircraft maneuvers. It is found that demixing can occur rather rapidly in space, although more slowly than on earth. The Ostwald ripening theory predicts demixing rates many orders of magnitude lower than the observed values.

Bamberger, S.↗
NASA NTRS2026

Micromechanical Design of Carbon Nanotube Ribbon Reinforced Polymer Composite Materials

Lightweight materials are an important component of the design of aerospace structures. Carbon nanotube materials have been considered for this purpose due to the strength and stiffness of individual nanotubes, and the commercial availability of bulk formats such as fibers. These fibers can have a ribbon cross section which results in a different design space for their composites relative to traditional reinforcements which have a round cross section. This work applies brick-and-mortar micromechanical models and classical lamination theory with an inverse approach to investigate the design space of these composites. Using this approach, the influences of fiber geometry and axial and transverse mechanical properties are mapped. Finally, a sensitivity study is performed and the relative impacts of ±10% variations in the constituent material and geometric properties are ranked. Lamina axial moduli were found to range from a maximum of 3x to 1x minimum relative to a target quasi-isotropic laminate modulus depending on the anisotropy and shear modulus of the lamina. The fiber targets depended strongly on the fiber volume fraction in the lamina and the fiber axial modulus target was found to range from 4.8x to 3.2x the quasi-isotropic laminate target. The sensitivity analysis found that the largest driver of performance was the volume fraction, followed by the fiber axial modulus. While bio-based brick-and-mortar composites, such as nacre, can benefit from reinforcement aspect ratios above 10, for carbon nanotube ribbon (or carbon fiber)/polymer composites the sensitivity study indicated that the optimal cross-sectional aspect ratio was relatively smaller, potentially less than three.

Micromechanics↗
NASA NTRS2016

Tool for Generation of MAC/GMC Representative Unit Cell for CMC/PMC Analysis

This document describes a recently developed analysis tool that enhances the resident capabilities of the Micromechanics Analysis Code with the Generalized Method of Cells (MAC/GMC) 4.0. This tool is especially useful in analyzing ceramic matrix composites (CMCs), where higher fidelity with improved accuracy of local response is needed. The tool, however, can be used for analyzing polymer matrix composites (PMCs) as well. MAC/GMC 4.0 is a composite material and laminate analysis software developed at NASA Glenn Research Center. The software package has been built around the concept of the generalized method of cells (GMC). The computer code is developed with a user friendly framework, along with a library of local inelastic, damage, and failure models. Further, application of simulated thermomechanical loading, generation of output results, and selection of architectures to represent the composite material have been automated to increase the user friendliness, as well as to make it more robust in terms of input preparation and code execution. Finally, classical lamination theory has been implemented within the software, wherein GMC is used to model the composite material response of each ply. Thus, the full range of GMC composite material capabilities is available for analysis of arbitrary laminate configurations as well. The primary focus of the current effort is to provide a graphical user interface (GUI) capability that generates a number of different user-defined repeating unit cells (RUCs). In addition, the code has provisions for generation of a MAC/GMC-compatible input text file that can be merged with any MAC/GMC input file tailored to analyze composite materials. Although the primary intention was to address the three different constituents and phases that are usually present in CMCs-namely, fibers, matrix, and interphase-it can be easily modified to address two-phase polymer matrix composite (PMC) materials where an interphase is absent. Currently, the tool capability includes generation of RUCs for square packing, hexagonal packing, and random fiber packing as well as RUCs based on actual composite micrographs. All these options have the fibers modeled as having a circular cross-sectional area. In addition, a simplified version of RUC is provided where the fibers are treated as having a square cross section and are distributed randomly. This RUC facilitates a speedy analysis using the higher fidelity version of GMC known as HFGMC. The first four mentioned options above support uniform subcell discretization. The last one has variable subcell sizes due to the primary intention of keeping the RUC size to a minimum to gain the speed ups using the higher fidelity version of MAC. The code is implemented within the MATLAB (The Mathworks, Inc., Natick, MA) developmental framework; however, a standalone application that does not need a priori MATLAB installation is also created with the aid of the MATLAB compiler.

Materials Engineering↗