This research summarizes the validation testing completed for the material model MAT213, currently implemented in the LS-DYNA finite element program. Testing was carried out using a carbon fiber composite material, T800-F3900. Stacked-ply tension and compression tests were performed for open-hole and full coupons. Comparisons of experimental and simulation results showed a good agreement between the two for metrics including, stress-strain response and displacements. Strains and displacements in the direction of loading were better predicted by the simulations than for that of the transverse direction.

Double cantilever beam and end notched flexure tests were performed experimentally and through simulations to determine the delamination properties of the material at the interlaminar layers. Experimental results gave the mode I critical energy release rate as having a range of 2.18 – 3.26 psi-in and the mode II critical energy release rate as 10.50 psi-in, both for the pre-cracked condition. Simulations were performed to calibrate other cohesive zone parameters required for modeling.

Samples of tested T800/F3900 coupons were processed and examined with scanning electron microscopy to determine and understand the underlying structure of the material. Tested coupons revealed damage and failure occurring at the micro scale for the composite material.

The exceptional mechanical properties of polymers with heterogeneous structure, such as the high toughness of polyethylene and the excellent blast-protection capability of polyurea, are strongly related to their morphology and nanoscale structure. Different polymer microstructures, such as semicrystalline morphology and segregated nanophases, lead to coordinated molecular motions during deformation in order to preserve compatibility between the different material phases. To study molecular relaxation in polyethylene, a coarse-grained model of polyethylene was calibrated to match the local structural variable distributions sampled from supercooled atomistic melts. The coarse-grained model accurately reproduces structural properties, e.g., the local structure of both the amorphous and crystalline phases, and thermal properties, e.g., glass transition and melt temperatures, and dynamic properties: including the vastly different relaxation time scales of the amorphous and crystalline phases. A hybrid Monte Carlo routine was developed to generate realistic semicrystalline configurations of polyethylene. The generated systems accurately predict the activation energy of the alpha relaxation process within the crystalline phase. Furthermore, the models show that connectivity to long chain segments in the amorphous phase increases the energy barrier for chain slip within crystalline phase. This prediction can guide the development of tougher semicrystalline polymers by providing a fundamental understanding of how nanoscale morphology contributes to chain mobility. In a different study, the macroscopic shock response of polyurea, a phase segregated copolymer, was analyzed using density functional theory (DFT) molecular dynamics (MD) simulations and classical MD simulations. The two models predict the shock response consistently up to shock pressures of 15 GPa, beyond which the DFT-based simulations predict a softer response. From the DFT simulations, an analysis of bond scission was performed as a first step in developing a more fundamental understanding of how shock induced material transformations effect the shock response and pressure dependent strength of polyurea subjected to extreme shocks.

Hydrodynamic phenomena such as the Rayleigh-Taylor (RT) and Richtmyer-Meshkov (RM) instabilities can be described by exponential/linear growth of surface perturbations at a bimaterial interface when subjected to constant/impulsive acceleration. A challenge in designing systems to mitigate or exploit these effects is the lack of accurate material models at large dynamic strain rates and pressures. In particular, little stress-strain constitutive information at large strain rates and pressures is available for transient material phases formed at high pressures, and the continuum effect the phase transformation process has on the instability evolution. In this work, a phase-aware isotropic strength model is developed and partially validated with a novel RM-based instability experiment in addition to existing data from the literature. With the validated material model additional simulations are performed to provide insight into to the role that robust material constitutive behavior (e.g., pressure, temperature, rate dependence) has on RM instability and how RM instability experiments can be used to characterize and validated expected material behavior. For phase aware materials, particularly iron in this work, the simulations predict a strong dependence on the Atwood number that single phase materials do not have. At Atwood numbers close to unity, and pressures in the high pressure stability region, the high pressure phase dominates the RM evolution. However, at Atwood numbers close to negative one, the RM evolution is only weakly affected by the high-pressure phase even for shocks well above the phase transformation threshold. In addition to RM evolution this work looks at the closely related shock front perturbation evolution. Existing analytical models for isentropic processes in gases and liquids are modified for metal equation of states and plastic behavior for the first time. It is found that the presence of a volume collapsing phase transformation with increased pressure causes shock front perturbations to decay sooner, while plastic strength has the opposite effect which is significantly different from the effect viscosity has. These results suggest additional experimental setups to validate material models, or relevant material parameters that can be optimized for system design objectives, e.g., minimize feed through perturbations in inertial confinement fusion capsules.

This study explores the possibility of two matrices containing metallic particulates to act as smart materials by sensing of strain due to the presence of the conducting particles in the matrix. The first matrix is a regular Portland cement-based one while the second is a novel iron-based, carbonated binder developed at ASU. Four different iron replacement percentages by volume (10%, 20%, 30% and 40%) in a Portland cement matrix were selected, whereas the best performing iron carbonate matrix developed was used. Electrical impedance spectroscopy was used to obtain the characteristic Nyquist plot before and after application of flexural load. Electrical circuit models were used to extract the changes in electrical properties under application of load. Strain sensing behavior was evaluated with respect to application of different stress levels and varying replacement levels of the inclusion. A similar approach was used to study the strain sensing capabilities of novel iron carbonate binder. It was observed that the strain sensing efficiency increased with increasing iron percentage and the resistivity increased with increase in load (or applied stress) for both the matrices. It is also found that the iron carbonate binder is more efficient in strain sensing as it had a higher gage factor when compared to the OPC matrix containing metallic inclusions.

Analytical equations (Maxwell) were used to extract frequency dependent electrical conductivity and permittivity of the cement paste (or the host matrix), interface, inclusion (iron) and voids to develop a generic electro-mechanical coupling model to for the strain sensing behavior. COMSOL Multiphysics 5.2a was used as finite element analysis software to develop the model. A MATLAB formulation was used to generate the microstructure with different volume fractions of inclusions. Material properties were assigned (the frequency dependent electrical parameters) and the coupled structural and electrical physics interface in COMSOL was used to model the strain sensing response. The experimental change in resistance matched well with the simulated values, indicating the applicability of the model to predict the strain sensing response of particulate composite systems.

Pultrusion manufacturing technique stands at the forefront for efficient production of continuous, uniform concrete composites for use in large scale structural applications. High volume and low labor, among other benefits such as improved impregnation and better sample consistency, stand as some of the crucial advances found in automated pultrusion. These advantages introduce textile reinforced concrete (TRC) composites as a potential surrogate for wood, light gauge steel, and other common structural materials into an ever changing and broadening market of industrial grade structural sections. With the potential modifications of textile geometry, textile type, section geometry, and connection type, the options presented by TRC sections seem nearly boundless. Automated pultrusion presents the ability to manufacture many different TRC composite types in at a quickened rate opening up a new field of study of structural materials.

The objective of this study centered on two studies including the development of an automated pultrusion system for the manufacturing of TRC composites and ultimately the assessment of composites created with the pultrusion technique and their viability as a relevant structural construction material. Upon planning, fabrication, and continued use of an automated pultrusion system in Arizona State University’s Structures Lab, an initial, comparative study of polypropylene microfiber composites was conducted to assess fiber reinforced concrete composites, manufactured with Filament Winding Technique, and textile reinforced concrete composites, manufactured with Automated Pultrusion Technique, in tensile and flexural mechanical response at similar reinforcement dosages. A secondary study was then conducted to measure the mechanical behavior of carbon, polypropylene, and alkali-resistant glass TRC composites and explore the response of full scale TRC structural shapes, including angle and channel sections. Finally, a study was conducted on the connection type for large scale TRC composite structural sections in tension and compression testing.

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.

Sustainable materials and methods have achieved a pivotal role in the research plethora of the new age due to global warming. Cement production is responsible in contributing to 5% of global CO2 emissions. Complete replacement of cement by alkaline activation of aluminosilicate waste materials such as slag and fly ash is a major advancement towards reducing the adverse impacts of cement production. Comprehensive research has been done, to understand the optimized composition and hydration products. The focus of this dissertation is to understand the multiscale behavior ranging from early age properties, fundamental material structure, fracture and crack resistance properties, durability responses and alternative activation methods to existing process.

The utilization of these materials has relied primarily on the dual benefits of reduced presence in landfills and cost. These have also proven to yield a higher service life as opposed to conventional ordinary portland cement (OPC) concrete due to an enhanced microstructure. The use of such materials however has not been readily acceptable due to detrimental early age behavior. The influence of design factors is studied to understand the reaction mechanism. Silicon polymerization at the molecular level is studied to understand the aluminosilicate interactions which are responsible for prevention of any leaching of ions. A comparative study between fly ash and slag binders is carried out to evaluate the stable states of sodium, aluminum and silicon in both these binders, since the likelihood of the sodium ions leaching out is high.

Compressive and flexural strength have been reported in previous literature, but the impact of crack resistance was unevaluated from an approach of characterizing the fracture process zone. Alternative routes of activation are explored with an intent to reduce the high alkalinity by use of neutral salts such as sodium sulfate. High volume OPC replacement by both class C and F fly ash is performed to evaluate the differences in hydration phase formation responsible for its pore refinement and strength. Spectroscopic studies have also allowed to study the fundamental material structure. Durability studies are also performed on these samples to understand the probability external sulfate attacks as opposed to OPC mixes.

Finite element simulations modeling the hydrodynamic impact loads subjected to an elastomeric coating were performed to develop an understanding of the performance and failure mechanisms of protective coatings for cavitating environments.

In this work, two major accomplishments were achieved: 1) scaling laws were developed from hydrodynamic principles and numerical simulations to allow conversion of measured distributions of pressure peaks in a cavitating flow to distributions of microscopic impact loadings modeling individual bubble collapse events, and 2) a finite strain, thermo-mechanical material model for polyurea-based elastomers was developed using a logarithmic rate formulation and implemented into an explicit finite element code.

Combining the distribution of microscopic impact loads and finite element modeling, a semi-quantitative predictive framework is created to calculate the energy dissipation within the coating which can further the understanding of temperature induced coating failures.

The influence of coating thickness and elastomer rheology on the dissipation of impact energies experienced in cavitating flows has also been explored.

The logarithmic formulation has many desired features for the polyurea constitutive model, such as objectivity, integrability, and additive decomposition compatibility.

A review and discussion on the kinematics in large deformation, including a comparison between Lagrangian and Eulerian descriptions, are presented to explain the issues in building rate-dependent constitutive models in finite strains.

When comparing the logarithmic rate with other conventional rates in test examples, the logarithmic rate shows a better conservation of objectivity and integrability.

The modeling framework was validated by comparing predictions against temperatures measured within coatings subjected to a cavitating jet.

Both the experiments and models show that the temperatures generated, even under mild flow conditions, raise the coating temperature by a significant amount, suggesting that the failure of these coatings under more aggressive flows is thermally induced.

The models show that thin polyurea coatings synthesized with shorter molecular weight soft segments dissipate significantly less energy per impact and conduct heat more efficiently.

This work represents an important step toward understanding thermally induced failure in elastomers subjected to cavitating flows, which provides a foundation for design and optimization of coatings with enhanced erosion resistance.

Tall building developments are spreading across the globe at an ever-increasing rate (www.ctbuh.org). In 1982, the number of ‘tall buildings’ in North America was merely 1,701. This number rose to 26,053, in 2006. The global number of buildings, 200m or more in height, has risen from 286 to 602 in the last decade alone. This dissertation concentrates on design optimization of such, about-to-be modular, structures by implementing AISC 2010 design requirements. Along with a discussion on and classification of lateral load resisting systems, a few design optimization cases are also being studied. The design optimization results of full scale three dimensional buildings subject to multiple design criteria including stress, serviceability and dynamic response are discussed. The tool being used for optimization is GS-USA Frame3D© (henceforth referred to as Frame3D). Types of analyses being verified against a strong baseline of Abaqus 6.11-1, are stress analysis, modal analysis and buckling analysis.

The provisions in AISC 2010 allows us to bypass the limit state of flexural buckling in compression checks with a satisfactory buckling analysis. This grants us relief from the long and tedious effective length factor computations. Besides all the AISC design checks, an empirical equation to check beams with high shear and flexure is also being enforced.

In this study, we present the details of a tool that can be useful in design optimization - finite element modeling, translating AISC 2010 design code requirements into components of the FE and design optimization models. A comparative study of designs based on AISC 2010 and fixed allowable stresses, (regardless of the shape of cross section) is also being carried out.

Composite materials are widely used in various structural applications, including within the automotive and aerospace industries. Unidirectional composite layups have replaced other materials such as metals due to composites’ high strength-to-weight ratio and durability. Finite-element (FE) models are actively being developed to model response of composite systems subjected to a variety of loads including impact loads. These FE models rely on an array of measured material properties as input for accuracy. This work focuses on an orthotropic plasticity constitutive model that has three components – deformation, damage and failure. The model relies on the material properties of the composite such as Young’s modulus, Poisson’s ratio, stress-strain curves in the principal and off-axis material directions, etc. This thesis focuses on two areas important to the development of the FE model – tabbing of the test specimens and data processing of the tests used to generate the required stress-strain curves. A comparative study has been performed on three candidate adhesives using double lap-shear testing to determine their effectiveness in composite specimen tabbing. These tests determined the 3M DP460 epoxy performs best in shear. The Loctite Superglue with 80% the ultimate stress of the 3M DP460 epoxy is acceptable when test specimens have to be ready for testing within a few hours. JB KwikWeld is not suitable for tabbing. In addition, the Experimental Data Processing (EDP) program has been improved for use in post-processing raw data from composites test. EDP has improved to allow for complete processing with the implementation of new weighted least squares smoothing options, curve averaging techniques, and new functionality for data manipulation.