Damage assessment and residual useful life estimation (RULE) are essential for aerospace, civil and naval structures. Structural Health Monitoring (SHM) attempts to automate the process of damage detection and identification. Multiscale modeling is a key element in SHM. It not only provides important information on the physics of failure, such as damage initiation and growth, the output can be used as "virtual sensing" data for detection and prognosis. The current research is part of an ongoing multidisciplinary effort to develop an integrated SHM framework for metallic aerospace components. In this thesis a multiscale model has been developed by bridging the relevant length scales, micro, meso and macro (or structural scale). Micro structural representations obtained from material characterization studies are used to define the length scales and to capture the size and orientation of the grains at the micro level. Parametric studies are conducted to estimate material parameters used in this constitutive model. Numerical and experimental simulations are performed to investigate the effects of Representative Volume Element (RVE) size, defect area fraction and distribution. A multiscale damage criterion accounting for crystal orientation effect is developed. This criterion is applied for fatigue crack initial stage prediction. A damage evolution rule based on strain energy density is modified to incorporate crystal plasticity at the microscale (local). Optimization approaches are used to calculate global damage index which is used for the RVE failure prediciton. Potential cracking directions are provided from the damage criterion simultaneously. A wave propagation model is incorporated with the damage model to detect changes in sensing signals due to plastic deformation and damage growth.

Structural health management (SHM) is emerging as a vital methodology to help engineers improve the safety and maintainability of critical structures. SHM systems are designed to reliably monitor and test the health and performance of structures in aerospace, civil, and mechanical engineering applications. SHM combines multidisciplinary technologies including sensing, signal processing, pattern recognition, data mining, high fidelity probabilistic progressive damage models, physics based damage models, and regression analysis. Due to the wide application of carbon fiber reinforced composites and their multiscale failure mechanisms, it is necessary to emphasize the research of SHM on composite structures. This research develops a comprehensive framework for the damage detection, localization, quantification, and prediction of the remaining useful life of complex composite structures. To interrogate a composite structure, guided wave propagation is applied to thin structures such as beams and plates. Piezoelectric transducers are selected because of their versatility, which allows them to be used as sensors and actuators. Feature extraction from guided wave signals is critical to demonstrate the presence of damage and estimate the damage locations. Advanced signal processing techniques are employed to extract robust features and information. To provide a better estimate of the damage for accurate life estimation, probabilistic regression analysis is used to obtain a prediction model for the prognosis of complex structures subject to fatigue loading. Special efforts have been applied to the extension of SHM techniques on aerospace and spacecraft structures, such as UAV composite wings and deployable composite boom structures. Necessary modifications of the developed SHM techniques were conducted to meet the unique requirements of the aerospace structures. The developed SHM algorithms are able to accurately detect and quantify impact damages as well as matrix cracking introduced.

Composite materials are increasingly being used in aircraft, automobiles, and other applications due to their high strength to weight and stiffness to weight ratios. However, the presence of damage, such as delamination or matrix cracks, can significantly compromise the performance of these materials and result in premature failure. Structural components are often manually inspected to detect the presence of damage. This technique, known as schedule based maintenance, however, is expensive, time-consuming, and often limited to easily accessible structural elements. Therefore, there is an increased demand for robust and efficient Structural Health Monitoring (SHM) techniques that can be used for Condition Based Monitoring, which is the method in which structural components are inspected based upon damage metrics as opposed to flight hours. SHM relies on in situ frameworks for detecting early signs of damage in exposed and unexposed structural elements, offering not only reduced number of schedule based inspections, but also providing better useful life estimates. SHM frameworks require the development of different sensing technologies, algorithms, and procedures to detect, localize, quantify, characterize, as well as assess overall damage in aerospace structures so that strong estimations in the remaining useful life can be determined. The use of piezoelectric transducers along with guided Lamb waves is a method that has received considerable attention due to the weight, cost, and function of the systems based on these elements. The research in this thesis investigates the ability of Lamb waves to detect damage in feature dense anisotropic composite panels. Most current research negates the effects of experimental variability by performing tests on structurally simple isotropic plates that are used as a baseline and damaged specimen. However, in actual applications, variability cannot be negated, and therefore there is a need to research the effects of complex sample geometries, environmental operating conditions, and the effects of variability in material properties. This research is based on experiments conducted on a single blade-stiffened anisotropic composite panel that localizes delamination damage caused by impact. The overall goal was to utilize a correlative approach that used only the damage feature produced by the delamination as the damage index. This approach was adopted because it offered a simplistic way to determine the existence and location of damage without having to conduct a more complex wave propagation analysis or having to take into account the geometric complexities of the test specimen. Results showed that even in a complex structure, if the damage feature can be extracted and measured, then an appropriate damage index can be associated to it and the location of the damage can be inferred using a dense sensor array. The second experiment presented in this research studies the effects of temperature on damage detection when using one test specimen for a benchmark data set and another for damage data collection. This expands the previous experiment into exploring not only the effects of variable temperature, but also the effects of high experimental variability. Results from this work show that the damage feature in the data is not only extractable at higher temperatures, but that the data from one panel at one temperature can be directly compared to another panel at another temperature for baseline comparison due to linearity of the collected data.

Recent studies of the occurrence of post-flutter limit cycle oscillations (LCO) of the F-16 have provided good support to the long-standing hypothesis that this phenomenon involves a nonlinear structural damping. A potential mechanism for the appearance of nonlinearity in the damping are the nonlinear geometric effects that arise when the deformations become large enough to exceed the linear regime. In this light, the focus of this investigation is first on extending nonlinear reduced order modeling (ROM) methods to include viscoelasticity which is introduced here through a linear Kelvin-Voigt model in the undeformed configuration. Proceeding with a Galerkin approach, the ROM governing equations of motion are obtained and are found to be of a generalized van der Pol-Duffing form with parameters depending on the structure and the chosen basis functions. An identification approach of the nonlinear damping parameters is next proposed which is applicable to structures modeled within commercial finite element software.

The effects of this nonlinear damping mechanism on the post-flutter response is next analyzed on the Goland wing through time-marching of the aeroelastic equations comprising a rational fraction approximation of the linear aerodynamic forces. It is indeed found that the nonlinearity in the damping can stabilize the unstable aerodynamics and lead to finite amplitude limit cycle oscillations even when the stiffness related nonlinear geometric effects are neglected. The incorporation of these latter effects in the model is found to further decrease the amplitude of LCO even though the dominant bending motions do not seem to stiffen as the level of displacements is increased in static analyses.

Advanced aerospace materials, including fiber reinforced polymer and ceramic matrix composites, are increasingly being used in critical and demanding applications, challenging the current damage prediction, detection, and quantification methodologies. Multiscale computational models offer key advantages over traditional analysis techniques and can provide the necessary capabilities for the development of a comprehensive virtual structural health monitoring (SHM) framework. Virtual SHM has the potential to drastically improve the design and analysis of aerospace components through coupling the complementary capabilities of models able to predict the initiation and propagation of damage under a wide range of loading and environmental scenarios, simulate interrogation methods for damage detection and quantification, and assess the health of a structure. A major component of the virtual SHM framework involves having micromechanics-based multiscale composite models that can provide the elastic, inelastic, and damage behavior of composite material systems under mechanical and thermal loading conditions and in the presence of microstructural complexity and variability. Quantification of the role geometric and architectural variability in the composite microstructure plays in the local and global composite behavior is essential to the development of appropriate scale-dependent unit cells and boundary conditions for the multiscale model. Once the composite behavior is predicted and variability effects assessed, wave-based SHM simulation models serve to provide knowledge on the probability of detection and characterization accuracy of damage present in the composite. The research presented in this dissertation provides the foundation for a comprehensive SHM framework for advanced aerospace materials. The developed models enhance the prediction of damage formation as a result of ceramic matrix composite processing, improve the understanding of the effects of architectural and geometric variability in polymer matrix composites, and provide an accurate and computational efficient modeling scheme for simulating guided wave excitation, propagation, interaction with damage, and sensing in a range of materials. The methodologies presented in this research represent substantial progress toward the development of an accurate and generalized virtual SHM framework.

Damage detection in heterogeneous material systems is a complex problem and requires an in-depth understanding of the material characteristics and response under varying load and environmental conditions. A significant amount of research has been conducted in this field to enhance the fidelity of damage assessment methodologies, using a wide range of sensors and detection techniques, for both metallic materials and composites. However, detecting damage at the microscale is not possible with commercially available sensors. A probable way to approach this problem is through accurate and efficient multiscale modeling techniques, which are capable of tracking damage initiation at the microscale and propagation across the length scales. The output from these models will provide an improved understanding of damage initiation; the knowledge can be used in conjunction with information from physical sensors to improve the size of detectable damage. In this research, effort has been dedicated to develop multiscale modeling approaches and associated damage criteria for the estimation of damage evolution across the relevant length scales. Important issues such as length and time scales, anisotropy and variability in material properties at the microscale, and response under mechanical and thermal loading are addressed. Two different material systems have been studied: metallic material and a novel stress-sensitive epoxy polymer.

For metallic material (Al 2024-T351), the methodology initiates at the microscale where extensive material characterization is conducted to capture the microstructural variability. A statistical volume element (SVE) model is constructed to represent the material properties. Geometric and crystallographic features including grain orientation, misorientation, size, shape, principal axis direction and aspect ratio are captured. This SVE model provides a computationally efficient alternative to traditional techniques using representative volume element (RVE) models while maintaining statistical accuracy. A physics based multiscale damage criterion is developed to simulate the fatigue crack initiation. The crack growth rate and probable directions are estimated simultaneously.

Mechanically sensitive materials that exhibit specific chemical reactions upon external loading are currently being investigated for self-sensing applications. The "smart" polymer modeled in this research consists of epoxy resin, hardener, and a stress-sensitive material called mechanophore The mechanophore activation is based on covalent bond-breaking induced by external stimuli; this feature can be used for material-level damage detections. In this work Tris-(Cinnamoyl oxymethyl)-Ethane (TCE) is used as the cyclobutane-based mechanophore (stress-sensitive) material in the polymer matrix. The TCE embedded polymers have shown promising results in early damage detection through mechanically induced fluorescence. A spring-bead based network model, which bridges nanoscale information to higher length scales, has been developed to model this material system. The material is partitioned into discrete mass beads which are linked using linear springs at the microscale. A series of MD simulations were performed to define the spring stiffness in the statistical network model. By integrating multiple spring-bead models a network model has been developed to represent the material properties at the mesoscale. The model captures the statistical distribution of crosslinking degree of the polymer to represent the heterogeneous material properties at the microscale. The developed multiscale methodology is computationally efficient and provides a possible means to bridge multiple length scales (from 10 nm in MD simulation to 10 mm in FE model) without significant loss of accuracy. Parametric studies have been conducted to investigate the influence of the crosslinking degree on the material behavior. The developed methodology has been used to evaluate damage evolution in the self-sensing polymer.

All structures suffer wear and tear because of impact, excessive load, fatigue, corrosion, etc. in addition to inherent defects during their manufacturing processes and their exposure to various environmental effects. These structural degradations are often imperceptible, but they can severely affect the structural performance of a component, thereby severely decreasing its service life. Although previous studies of Structural Health Monitoring (SHM) have revealed extensive prior knowledge on the parts of SHM processes, such as the operational evaluation, data processing, and feature extraction, few studies have been conducted from a systematical perspective, the statistical model development.

The first part of this dissertation, the characteristics of inverse scattering problems, such as ill-posedness and nonlinearity, reviews ultrasonic guided wave-based structural health monitoring problems. The distinctive features and the selection of the domain analysis are investigated by analytically searching the conditions of the uniqueness solutions for ill-posedness and are validated experimentally.

Based on the distinctive features, a novel wave packet tracing (WPT) method for damage localization and size quantification is presented. This method involves creating time-space representations of the guided Lamb waves (GLWs), collected at a series of locations, with a spatially dense distribution along paths at pre-selected angles with respect to the direction, normal to the direction of wave propagation. The fringe patterns due to wave dispersion, which depends on the phase velocity, are selected as the primary features that carry information, regarding the wave propagation and scattering.

The following part of this dissertation presents a novel damage-localization framework, using a fully automated process. In order to construct the statistical model for autonomous damage localization deep-learning techniques, such as restricted Boltzmann machine and deep belief network, are trained and utilized to interpret nonlinear far-field wave patterns.

Next, a novel bridge scour estimation approach that comprises advantages of both empirical and data-driven models is developed. Two field datasets from the literature are used, and a Support Vector Machine (SVM), a machine-learning algorithm, is used to fuse the field data samples and classify the data with physical phenomena. The Fast Non-dominated Sorting Genetic Algorithm (NSGA-II) is evaluated on the model performance objective functions to search for Pareto optimal fronts.

There are many applications for polymer matrix composite materials in a variety of different industries, but designing and modeling with these materials remains a challenge due to the intricate architecture and damage modes. Multiscale modeling techniques of composite structures subjected to complex loadings are needed in order to address the scale-dependent behavior and failure. The rate dependency and nonlinearity of polymer matrix composite materials further complicates the modeling. Additionally, variability in the material constituents plays an important role in the material behavior and damage. The systematic consideration of uncertainties is as important as having the appropriate structural model, especially during model validation where the total error between physical observation and model prediction must be characterized. It is necessary to quantify the effects of uncertainties at every length scale in order to fully understand their impact on the structural response. Material variability may include variations in fiber volume fraction, fiber dimensions, fiber waviness, pure resin pockets, and void distributions. Therefore, a stochastic modeling framework with scale dependent constitutive laws and an appropriate failure theory is required to simulate the behavior and failure of polymer matrix composite structures subjected to complex loadings. Additionally, the variations in environmental conditions for aerospace applications and the effect of these conditions on the polymer matrix composite material need to be considered. The research presented in this dissertation provides the framework for stochastic multiscale modeling of composites and the characterization data needed to determine the effect of different environmental conditions on the material properties. The developed models extend sectional micromechanics techniques by incorporating 3D progressive damage theories and multiscale failure criteria. The mechanical testing of composites under various environmental conditions demonstrates the degrading effect these conditions have on the elastic and failure properties of the material. The methodologies presented in this research represent substantial progress toward understanding the failure and effect of variability for complex polymer matrix composites.

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.

Multiaxial mechanical fatigue of heterogeneous materials has been a significant cause of concern in the aerospace, civil and automobile industries for decades, limiting the service life of structural components while increasing time and costs associated with inspection and maintenance. Fiber reinforced composites and light-weight aluminum alloys are widely used in aerospace structures that require high specific strength and fatigue resistance. However, studying the fundamental crack growth behavior at the micro- and macroscale as a function of loading history is essential to accurately predict the residual fatigue life of components and achieve damage tolerant designs. The issue of mechanical fatigue can be tackled by developing reliable in-situ damage quantification methodologies and by comprehensively understanding fatigue damage mechanisms under a variety of complex loading conditions. Although a multitude of uniaxial fatigue loading studies have been conducted on light-weight metallic materials and composites, many service failures occur from components being subjected to variable amplitude, mixed-mode multiaxial fatigue loadings. In this research, a systematic approach is undertaken to address the issue of fatigue damage evolution in aerospace materials by:

(i) Comprehensive investigation of micro- and macroscale crack growth behavior in aerospace grade Al 7075 T651 alloy under complex biaxial fatigue loading conditions. The effects of variable amplitude biaxial loading on crack growth characteristics such as crack acceleration and retardation were studied in detail by exclusively analyzing the influence of individual mode-I, mixed-mode and mode-II overload and underload fatigue cycles in an otherwise constant amplitude mode-I baseline load spectrum. The micromechanisms governing crack growth behavior under the complex biaxial loading conditions were identified and correlated with the crack growth behavior and fracture surface morphology through quantitative fractography.

(ii) Development of novel multifunctional nanocomposite materials with improved fatigue resistance and in-situ fatigue damage detection and quantification capabilities. A state-of-the-art processing method was developed for producing sizable carbon nanotube (CNT) membranes for multifunctional composites. The CNT membranes were embedded in glass fiber laminates and in-situ strain sensing and damage quantification was achieved by exploiting the piezoresistive property of the CNT membrane. In addition, improved resistance to fatigue crack growth was observed due to the embedded CNT membrane.