The purpose of this project focuses on analyzing how a typically brittle material, such as PLA, can be manipulated to become deformable, through the development of an origami structure, in this case—the Yoshimuri pattern. The experimental methodology focused on creating a base Solidworks model, with varying hinge depths, and 3D printing these various models. A cylindrical shell was also developed with comparable dimensions to the Yoshimuri dimensions. These samples were then tested through compression testing, with the load-displacement, and thus the stress-strain curves are analyzed. From the results, it was found that generally, the Yoshimuri samples had a higher level of deformation compared to the cylindrical shell. Moreover, the cylindrical shell had a higher stiffness ratio, while the Yoshimuri patterns had strain rates as high as 16%. From this data, it can be concluded that by changing how the structure is created through origami patterns, it is possible to shift the characteristics of a structure even if the material properties are initially quite brittle.

Electrospun nanofibers can be prepared from various kinds of inorganic substances by electro-spinning techniques. They have great potential in many applications including super capacitors, lithium ion batteries, filtration, catalyst and enzyme carriers, and sensors [1]. The traditional way to produce electrospun nanofibers is needle based electro-spinning [1]. However, electrospun nanofibers have not been widely used in practice because of low nanofiber production rates. One way to largely increase the electro-spinning productivity is needleless electro-spinning. In 2005, Jirsak et al. patented a rotating roller fiber generator for the mass production of nanofibers [2]. Elmarco Corporation commercialized this technique to manufacture nanofiber equipment for the production of all sorts of organic and inorganic nanofibers, and named it "NanospiderTM". For this project, my goal is to build a needleless electro-spinner to produce nanofibers as the separator of lithium ion batteries. The model of this project is based on the design of rotating roller fiber generator, and is adapted from a project at North Dakota State University in 2011 [3].

The data and results presented in this paper are part of a continuing effort to innovate and pioneer the future of engineering. The purpose of the following is to demonstrate the mechanical buckling characteristics in stiff thin film and soft substrate systems, and the importance of controlling them. In today's engineering research, wrinkling in systems in beginning to be viewed as a means for engineering innovation rather than failure. This research is important to further progress the possible applications the technology proposes, such as flexible electronics and tunable adhesives. This work utilizes a cost efficient and relatively easy method for generating and analyzing buckled systems. Ultra violate oxidation at ambient temperatures is exploited to create a stiff thin surface on rubber like polydimethylsiloxane, and couple with strain induction wrinkles are generated. Wrinkle characteristics such as amplitude, wavelengths and wetting properties were investigated. In simple cases, trends were confirmed that increased oxidation relates to increased buckle wavelengths, and increase in strain corresponds to a decrease in wavelength. Hierarchical buckles were produced in one-dimensional systems treated with a multi-step method; these were the first to be generated in the ASU labs. Unique topographic changes were produced in two-dimensional systems treated with the same method. Honeycomb or dome like structures were noted to occur, important as they undergo a different energy-reliving configuration compared to traditional parallel buckles. The information provided characterizes many aspects of the buckle phenomena and will allow for further inquiry into specific functions utilizing the technology to continue advancements in engineering.

A novel strain sensing procedure using an optical scanning methodology and diffraction grating is explored. The motivation behind this study is due to uneven thermal strain distribution across semiconductor chips that are composed of varying materials. Due to the unique properties of the materials and the different coefficients of thermal expansion (CTE), one can expect the material that experiences the highest strain to be the most likely failure point of the chip. As such, there is a need for a strain sensing technique that offers a very high strain sensitivity, a high spatial resolution while simultaneously achieving a large field of view. This study goes through the optical setup as well as the evolution of the optical grating in an effort to improve the strain sensitivity of this setup.

In the last several years, there has been a significant growth in research in the field of power harvesting, the process of capturing the energy surrounding a system and converting it into usable electrical energy. This concept has received particular interest in recent years with the ever-increasing production of portable and wireless electronic devices. Many of these devices that are currently in production utilize electrochemical batteries as a power source, which while effective, maintain the drawback of having a finite energy supply, thus requiring periodic replacement. The concept of power harvesting, however, works to solve these issues through electronics that are designed to capture ambient energy surrounding them convert it into usable electronic energy. The use of power harvesting in energy scavenging devices allows for the possible development of devices that are self-powered and do not require their power sources to be replaced. Several models have been developed by Soldano et al [3] and Liao et al [2] that have been proven accurate at predicting the power output of a piezoelectric power harvester in a cantileaver beam configuration. The work in this paper will expand further on the model developed by Liao et al [2], and as its main goal will use a modified form of that model to predict the optimal dimensions for a beam power harvester to achieve the maximum power output possible. The model will be updated b replacing the mode shape function used to approximate the deflection of the beam with the true defletion, which is based on the complex wavenumber that incorporates the complex Young's modulus of the material used. Other changes to account for this replacement will also be presented, along with numerical results of the final model.

Poly(ionic liquid)s (PILs) with an intrinsically conducting pyrrole polymer (ICP) backbone were synthesized and utilized as novel dispersants of carbon nanotubes (CNTs) in various polar and nonpolar solvents. This is due to their highly tunable nature, in which the anions can be easily exchanged to form PILs of varying polarity but with the same polycation. These CNT dispersions were exceedingly stable over many months, and with the addition of hexane, Pickering emulsions with the PIL-stabilized CNTs at the droplet interfaces were formed. Depending on the hydrophobicity of the PIL, hexane-in-water and hexane-in-acetonitrile emulsions were formed, the latter marking the first non-aqueous stabilized-CNT emulsions and corresponding CNT-in-acetonitrile dispersion, further advancing the processability of CNTs. The PIL-stabilized CNT Pickering emulsion droplets generated hollow conductive particles by subsequent drying of the emulsions. With the emulsion templating, the hollow shells can be used as a payload carrier, depending on the solubility of the payload in the droplet phase of the emulsion. This was demonstrated with silicon nanoparticles, which have limited solubility in aqueous environments, but great scientific interest due to their potential electrochemical applications. Overall, this work explored a new class of efficient PIL-ICP hybrid stabilizers with tunable hydrophobicity, offering extended stability of carbon nanotube dispersions with novel applications in hollow particle formation via Pickering emulsion templating and in placing payloads into the shells.

A hybrid molecular dynamics (MD) simulation framework is developed to emulate mechanochemical reaction of mechanophores in epoxy-based nanocomposites. Two different force fields, a classical force field and a bond order based force field are hybridized to mimic the experimental processes from specimen preparation to mechanical loading test. Ultra-violet photodimerization for mechanophore synthesis and epoxy curing for thermoset polymer generation are successfully simulated by developing a numerical covalent bond generation method using the classical force field within the framework. Mechanical loading tests to activate mechanophores are also virtually conducted by deforming the volume of a simulation unit cell. The unit cell deformation leads to covalent bond elongation and subsequent bond breakage, which is captured using the bond order based force field. The outcome of the virtual loading test is used for local work analysis, which enables a quantitative study of mechanophore activation. Through the local work analysis, the onset and evolution of mechanophore activation indicating damage initiation and propagation are estimated; ultimately, the mechanophore sensitivity to external stress is evaluated. The virtual loading tests also provide accurate estimations of mechanical properties such as elastic, shear, bulk modulus, yield strain/strength, and Poisson’s ratio of the system. Experimental studies are performed in conjunction with the simulation work to validate the hybrid MD simulation framework. Less than 2% error in estimations of glass transition temperature (Tg) is observed with experimentally measured Tgs by use of differential scanning calorimetry. Virtual loading tests successfully reproduce the stress-strain curve capturing the effect of mechanophore inclusion on mechanical properties of epoxy polymer; comparable changes in Young’s modulus and yield strength are observed in experiments and simulations. Early damage signal detection, which is identified in experiments by observing increased intensity before the yield strain, is captured in simulations by showing that the critical strain representing the onset of the mechanophore activation occurs before the estimated yield strain. It is anticipated that the experimentally validated hybrid MD framework presented in this dissertation will provide a low-cost alternative to additional experiments that are required for optimizing material design parameters to improve damage sensing capability and mechanical properties.

In addition to the study of mechanochemical reaction analysis, an atomistic model of interphase in carbon fiber reinforced composites is developed. Physical entanglement between semi-crystalline carbon fiber surface and polymer matrix is captured by introducing voids in multiple graphene layers, which allow polymer matrix to intertwine with graphene layers. The hybrid MD framework is used with some modifications to estimate interphase properties that include the effect of the physical entanglement. The results are compared with existing carbon fiber surface models that assume that carbon fiber has a crystalline structure and hence are unable to capture the physical entanglement. Results indicate that the current model shows larger stress gradients across the material interphase. These large stress gradients increase the viscoplasticity and damage effects at the interphase. The results are important for improved prediction of the nonlinear response and damage evolution in composite materials.

The Very High Temperature Reactor (VHTR) is one of six conceptual designs proposed for Generation IV nuclear reactors. Alloy 617, a solid solution strengthened Ni-base superalloy, is currently the primary candidate material for the tubing of the Intermediate Heat Exchanger (IHX) in the VHTR design. Steady-state operation of the nuclear power plant at elevated temperatures leads to creep deformation, whereas loading transients including startup and shutdown generate fatigue. A detailed understanding of the creep-fatigue interaction in Alloy 617 is necessary before it can be considered as a material for nuclear construction in ASME Boiler and Pressure Vessel Code. Current design codes for components undergoing creep-fatigue interaction at elevated temperatures require creep-fatigue testing data covering the entire range from fatigue-dominant to creep-dominant loading. Classical strain-controlled tests, which produce stress relaxation during the hold period, show a saturation in cycle life with increasing hold periods due to the rapid stress-relaxation of Alloy 617 at high temperatures. Therefore, applying longer hold time in these tests cannot generate creep-dominated failure. In this study, uniaxial isothermal creep-fatigue tests with non-traditional loading waveforms were designed and performed at 850 and 950°C, with an objective of generating test data in the creep-dominant regime. The new loading waveforms are hybrid strain-controlled and force-controlled testing which avoid stress relaxation during the creep hold. The experimental data showed varying proportions of creep and fatigue damage, and provided evidence for the inadequacy of the widely-used time fraction rule for estimating creep damage under creep-fatigue conditions. Micro-scale damage features in failed test specimens, such as fatigue cracks and creep voids, were quantified using a Scanning Electron Microscope (SEM) to find a correlation between creep and fatigue damage. Quantitative statistical imaging analysis showed that the microstructural damage features (cracks and voids) are correlated with a new mechanical driving force parameter. The results from this image-based damage analysis were used to develop a phenomenological life-prediction methodology called the effective time fraction approach. Finally, the constitutive creep-fatigue response of the material at 950°C was modeled using a unified viscoplastic model coupled with a damage accumulation model. The simulation results were used to validate an energy-based constitutive life-prediction model, as a mechanistic model for potential component and structure level creep-fatigue analysis.

This investigation develops small-size reduced order models (ROMs) that provide an accurate prediction of the response of only part of a structure, referred to as component-centric ROMs. Four strategies to construct such ROMs are presented, the first two of which are based on the Craig-Bampton Method and start with a set of modes for the component of interest (the β component). The response in the rest of the structure (the α component) induced by these modes is then determined and optimally represented by applying a Proper Orthogonal Decomposition strategy using Singular Value Decomposition. These first two methods are effectively basis reductions techniques of the CB basis. An approach based on the “Global - Local” Method generates the “global” modes by “averaging” the mass property over α and β comp., respectively (to extract a “coarse” model of α and β) and the “local” modes orthogonal to the “global” modes to add back necessary “information” for β. The last approach adopts as basis for the entire structure its linear modes which are dominant in the β component response. Then, the contributions of other modes in this part of the structure are approximated in terms of those of the dominant modes with close natural frequencies and similar mode shapes in the β component. In this manner, the non-dominant modal contributions are “lumped” onto the dominant ones, to reduce the number of modes for a prescribed accuracy. The four approaches are critically assessed on the structural finite element model of a 9-bay panel with the modal lumping-based method leading to the smallest sized ROMs. Therefore, it is extended to the nonlinear geometric situation and first recast as a rotation of the modal basis to achieve unobservable modes. In the linear case, these modes completely disappear from the formulation owing to orthogonality. In the nonlinear case, however, the generalized coordinates of these modes are still present in the nonlinear terms of the observable modes. A closure-type algorithm is then proposed to eliminate the unobserved generalized coordinates. This approach, its accuracy and computational savings, was demonstrated on a simple beam model and the 9-bay panel model.

A method for modelling the interactions of dislocations with inclusions has been developed to analyse toughening mechanisms in alloys. This method is different from the superposition method in that infinite domain solutions and image stress fields are not superimposed. The method is based on the extended finite element method (XFEM) in which the dislocations are modelled according to the Volterra dislocation model. Interior discontinuities are introduced across dislocation glide planes using enrichment functions and the resulting boundary value problem is solved through the standard finite element variational approach. The level set method is used to describe the geometry of the dislocation glide planes without any explicit treatment of the interface geometry which provides a convenient and an appealing means for describing the dislocation. A method for estimating the Peach-Koehler force by the domain form of J-integral is considered. The convergence and accuracy of the method are studied for an edge dislocation interacting with a free surface where analytical solutions are available. The force converges to the exact solution at an optimal rate for linear finite elements. The applicability of the method to dislocation interactions with inclusions is illustrated with a system of Aluminium matrix containing Aluminium-copper precipitates. The effect of size, shape and orientation of the inclusions on an edge dislocation for a difference in stiffness and coefficient of thermal expansion of the inclusions and matrix is considered. The force on the dislocation due to a hard inclusion increased by 8% in approaching the sharp corners of a square inclusion than a circular inclusion of equal area. The dislocation experienced 24% more force in moving towards the edges of a square shaped inclusion than towards its centre. When the areas of the inclusions were halved, 30% less force was exerted on the dislocation. This method was used to analyse interfaces with mismatch strains. Introducing eigenstrains equal to 0.004 to the elastic mismatch increased the force by 15 times for a circular inclusion. The energy needed to move an edge dislocation through a domain filled with circular inclusions is 4% more than that needed for a domain with square shaped inclusions.