A technique was developed for analyzing and developing mechanical properties for porous geomaterials subjected to the high pressures encountered in penetration and blast-type loadings. A finite element (FE) code was developed to verify laboratory test results or to predict unavailable laboratory test data for porous media loaded to multikilobar stress levels. This FE program eliminates a deficiency in the process of analyzing and developing mechanical properties for porous geomaterials by furnishing an advanced analysis tool to the engineer providing properties to material modelers or ground shock calculators. The FE code simulates quasi-static, axisymmetric, laboratory mechanical property tests, i.e., the laboratory tests are analyzed as boundary value problems. The code calculates strains, total and effective stresses, and pore fluid pressures for fully- and partially-saturated porous media. The time dependent flow of the pore fluid is also calculated. An elastic-plastic strain-hardening cap model calculates the time-independent skeletal responses of the porous solids. This enables the code to model nonlinear irreversible stress-strain behavior and shear-induced volume changes. Fluid and solid compressibilities were incorporated into the code, and partially-saturated materials were simulated with a "homogenized" compressible pore fluid. Solutions for several verification problems are given as proof that the program works correctly, and numerical simulations of limestone behavior under drained and undrained boundary conditions are also presented. / Ph. D.

Sizings, a thin polymer coating applied to the surface of the carbon fiber before impregnation with the matrix, have been shown to affect the mechanical performance of the composite. These sizings affect the processability of the carbon fiber that translates into a composite with less fiber breakage and improved fiber/matrix adhesion. In addition, the interdiffusion of the sizing and the bulk matrix results in the formation of an interphase. This interphase can alter damage initiation and propagation that can ultimately affect composite performance. The overall objective of the work detailed in this thesis is to ascertain the effects that thermoplastic sizing agents have on composite performance and determine the phenomenological events associated with the effects. All of the thermoplastic sizings had improved processability over the traditional G' sizing. These improvements in processability translated into a composite with less fiber damage and improved surface quality. In addition, all of the thermoplastic sizings outperformed the industrial benchmark sizing G' by at least 25% in static tensile strength, 11% in longitudinal flexure strength, and 30% in short beam shear strength. All moduli were found to be unaffected by the addition of a sizing. The interphase formed in K-90 PVP sized carbon fiber composites was fundamentally predicted from the constitutive properties of K-90 PVP/Derakane™ interdiffusion and fundamental mass transport equations. The K-90 PVP sizing material interdiffusing with the Derakane™ matrix was found to be dissolution controlled. The dissolution diffusion coefficient had an exponential concentration dependence. Fundamental mass transport models were utilized to predict the interphase profile. The predicted K-90 PVP interphase concentration profile displayed steep gradients at the fiber/matrix interface but essentially no gradients at points distant from the fiber surface. The predicted mechanical property profile was essentially flat for the modulus but did show a steep gradient in the strain-to-failure and shrinkage properties. However, the K-90 PVP interphase compared to the unsized/pure Derakane™ interphase showed improvements in strength and strain-to-failure and a reduction in cure shrinkage without significantly affecting the interphase tensile or shear moduli. / Ph. D.

A two part study is performed to investigate the application of reliability-based design optimization (RBDO) approach to design elastic-plastic stiffener beams with Tsection. The objectives of this study are to evaluate the benefits of reliability-based optimization over deterministic optimization, and to illustrate through a practical design example some of the difficulties that a design engineer may encounter while performing reliability-based optimization. Other objectives are to search for a computationally economic RBDO method and to utilize that method to perform RBDO to design an elastic-plastic T-stiffener under combined loads and with flexural-torsional buckling and local buckling failure modes. First, a nonlinear elastic-plastic T-beam was modeled using a simple 6 degree-of-freedom non-linear beam element. To address the problems of RBDO, such as the high non-linearity and derivative discontinuity of the reliability function, and to illustrate a situation where RBDO fails to produce a significant improvement over the deterministic optimization, a graphical method was developed. The method started by obtaining a deterministic optimum design that has the lowest possible weight for a prescribed safety factor (SF), and based on that design, the method obtains an improved optimum design that has either a higher reliability or a lower weight or cost for the same level of reliability as the deterministic design. Three failure modes were considered for an elastic-plastic beam of T cross-section under combined axial and bending loads. The failure modes are based on the total plastic failure in a beam section, buckling, and maximum allowable deflection. The results of the first part show that it is possible to get improved optimum designs (more reliable or lighter weight) using reliability-based optimization as compared to the design given by deterministic optimization. Also, the results show that the reliability function can be highly non-linear with respect to the design variables and with discontinuous derivatives. Subsequently, a more elaborate 14-degrees-of-freedom beam element was developed and used to model the global failure modes, which include the flexural-torsional and the out-of-plane buckling modes, in addition to local buckling modes. For this subsequent study, four failure modes were specified for an elasticplastic beam of T-cross-section under combined axial, bending, torsional and shear loads. These failure modes were based on the maximum allowable in-plane, out-ofplane and axial rotational deflections, in addition, to the web-tripping local buckling. Finally, the beam was optimized using the sequential optimization with reliabilitybased factors of safety (SORFS) RBDO technique, which was computationally very economic as compared to the widely used nested optimization loop techniques. At the same time, the SOPSF was successful in obtaining superior designs than the deterministic optimum designs (either up to12% weight savings for the same level of safety, or up to six digits improvement in the reliability for the same weight for a design with Safety Factor 2.50). / Ph. D.

The vestibular system consists of several organs that contribute to ones sense of balance. One set of organs, otoconial organs, have been shown to respond to linear acceleration (1949). Hair bundles (and hair cells), which are the mechano-electric transducers found within otoconial organs, respond to displacement of the overlying otoconial membrane (OM). Structure, position and orientation of the OM within the head may influence the stimulus of hair bundles by changing the deformation characteristics of the OM. Therefore, studying the deformation characteristics of the OM with finite element models presents a unique advantage: the ability to study how different variables may influence the deformation of the OM. Previous OM models have ignored complicated OM geometry in favor of single degree of freedom (De Vries 1951)or distributed parameter models (Grant et al. 1984; Grant and Cotton 1990; Grant et al. 1994). Additionally, OMs have been modeled considering three dimensional geometry (Benser et al. 1993; Kondrachuk 2000; 2001a), however OM layer thicknesses were assumed to be constant. Further, little research has investigated the effect of position and orientation of otoconial organs on the deformation of the OM (Curthoys et al. 1999), due to natural movement of the head. The effect of structure, position and orientation of the utricle of a red ear slider turtle on the stimulation of hair bundles in the OM is investigated here. Using confocal images, a finite element model of the utricle OM is constructed considering its full 3D geometry and varying OM layer thickness. How specific geometric variables, which are missing from other OM models, effect the deformation of the utricle OM is studied. Next, since hair bundles are part of the structure of the OM, their contribution to the deformation of the utricular OM is quantified. Then, using computed tomography of a turtle head and high speed video of turtle feeding strikes, acceleration at the utricle during natural motion is estimated. Finally, the effects of orientation of the utricle in the head on the stimulus of hair bundles within the organ is investigated. In summary, a model and methods are developed through which deformation of the turtle utricle OM through natural movements of the head may be studied. Variables that may contribute to utricle OM deformation are investigated. Utricle OM geometry, hair bundles, position and orientation all play a role in utricle OM deflection and therefore hair bundle stimulus. Their effects are quantified and their roles are discussed in this dissertation. / Ph. D.

Polymeric foam cores are widely used as core materials in sandwich panels subject to blast loads, where high strain rates of the order of 4000 /s are observed. Unlike metallic foams polymeric foams exhibit transversely isotropic response when tested in a laboratory setting. More specifically, they exhibit different hardening along the foam thickness than that in a direction transverse to the thickness. Furthermore, polymeric foams harden differently in tension and compression. In this thesis we adopt ideas from the constitutive model developed by Hoo Fatt et al. cite{hoofatt2}, which captures strain hardening, transverse isotropy and distinguishes the response in tension and in compression, to include isotropic strain rate hardening in our constitutive model. A one dimensional prototype of the model is used to aid in the physical explanation of various variables, and the model is generalized to three dimensions. The material model is implemented as a VUMAT (user defined) subroutine in the commercial finite element software ABAQUS Explicit. We show that the model works robustly in uniaxial deformations as well as in sandwich problems using the test data available in the literature. We provide values of the 39 material parameters for H45, H60, H80, H100, H130 and H200 foams. The constitutive relation is utilized in an optimization problem in which the surrogate optimizer is utilized to minimize the backface deflection of a blast loaded clamped sandwich plate of a fixed mass. The core in the optimized sandwich structure has a stratified configuration (not functionally graded) and has 24% less maximum back face deflection as compared to that in which the six core layers vary from highest density to lowest density or vice a versa. For a sandwich panel subject to a blast load, when the strain rate hardening effect are neglected, we observed a 12% reduction in the predicted peak deflection from that when strain rate effects are considered. It is counter intuitive and needs further investigation. / Master of Science / Sandwich panels are widely used in high performance structures requiring high stiffness, low weight and the ability to withstand blasts. Sandwich panels consist of several layers, and it is possible to vary the material and thickness of each layer to arrive at a sandwich panel design which performs optimally. In this thesis, we numerically find an optimum sandwich panel design so that it deflects the least when exposed to a given blast. The problem is studied using ABAQUS Explicit and the Surrogate Optimization solver built into MATLAB. The outer layers of the sandwich panel are made of a highly stiff material and their thicknesses are fixed. The remaining six inner layers are allowed to be any of six different H45, H60, H80, H100, H130, H200 Divinycell polymeric foams and are allowed to vary in thickness. In order to draw a fair comparison between the designs, we constrain the total mass of the sandwich panel to be 1 kg. In our quest to find the best sandwich panel design, we develop and implement, in ABAQUS Explicit, a custom mathematical model which captures the complex behavior of the polymeric foams. Experimental data in the literature and other techniques were utilized to check that this mathematical model accurately predicts the physical response of polymeric foams in different scenarios. The reader is given all of the theory and physical constants needed to use this mathematical model for the six foams. The optimal sandwich panel deflects 24% less than a baseline design, and it is found that the material properties of the six foams do not vary gradually as they did in the baseline designs.

Surface charging is a fundamental interaction process in space plasma engineering. A three-dimensional Immersed Finite Element Particle-In-Cell (IFE-PIC) method is developed to model surface charging involving complex boundary conditions. This method extends the previous IFE-PIC algorithm to explicitly include charge deposition on a dielectric surface for charging calculations. Three simulation studies are carried out using the new algorithm to model current collection and charging in both the orbital motion limited (OML) and space charge limited regime. The first one is a full particle simulation of the charging process of single small sphere and clusters of multiple small spheres in plasma. We find that while single sphere charging agrees well with the predictions of the OML theory, the charging of a sphere in a cluster is significantly, indicating that the often used OML charging model is not an accurate one to model charging in dusty plasma. The second one concerns a secondary electron emission experiment. The simulation includes detailed experimental setup in a vacuum chamber and the results are compared against experimental data. The simulation is used to determine the facility error in experiments. The third one is a full particle simulation of charging on lunar surface. The simulation concerns both flat and non-flat surface, and spacecraft on lunar surface, in the lunar polar region. The surface sees a mesothermal solar wind plasma flow and the emission of photoelectrons and secondary electrons. At a small sun elevation angle, the surface landscape generates a complex plasma flow field and local differential charging on surface. The results will be useful for further study of charging and levitation of lunar dust. / Ph. D.

Inflatable dams are flexible, cylindrical structures anchored to a foundation. They are used for a variety of purposes, e.g. diverting water for irrigation or groundwater recharging, impounding water for recreational purposes, and raising the height of existing dams or spillways. The vibration behavior of such dams is analyzed. Single-anchor inflatable dams with fins are considered. First, a static analysis is performed which yields the equilibrium shapes of the dam, both in the presence and absence of water. Then, a dynamic analysis is undertaken which analyzes the small vibrations of the inflatable dam about the equilibrium configuration, both in the presence of water (hydrostatic water as well as parallel flowing water) and absence of water. The dam is modeled as an elastic shell. It is assumed to be air-inflated and resting on a rigid foundation. The cross-sectional perimeter, material thickness, modulus of elasticity, and Poisson's ratio are given. The analysis is performed for different values of internal pressure and external water heads. Initially, the dam is assumed to lie flat. The internal pressure is then increased slowly until it reaches the desired value. Then the external water is applied and the equilibrium configuration is obtained. Small vibrations about this configuration are considered. The water is assumed to be inviscid and incompressible, and potential theory is used. The infinite-frequency limit is assumed on the free surface. A boundary element technique is utilized to determine the behavior of the water, and the finite element program ABAQUS is used to analyze the structural behavior. Both the cases of fluid at rest and flowing parallel to the dam are considered. The vibration frequencies and mode shapes are computed. The effect of the internal pressure of the dam is investigated, and the results are compared to those for the dam in the absence of external water. / Master of Science

An important part of an aerodynamic design process is optimizing designs to maximize quantities such as lift and the lift-to-drag ratio, in a process known as shape optimization. It is the goal of this thesis to develop and apply understanding of mixed finite element method and sensitivity analysis in a way that sets the foundation for shape optimization. The open-source Incompressible Flow Iterative Solution Software (IFISS) mixed finite element method toolbox for MATLAB developed by Silvester, Elman, and Ramage is used. Meshes are produced for a backward-facing step problem, using built-in tools from IFISS as well as the mesh generation software Gmsh, and grid convergence studies are performed for both sets of meshes along a sampled data line to ensure that the simulations converge asymptotically with increasing mesh resolution. As a preliminary study of sensitivity analysis, analytic sensitivities of velocity components along the backward-facing step data line to inflow velocity parameters are determined and verified using finite difference and complex step sensitivity values. The method is then applied to pressure drag calculated by integrating the pressure over the surface of a circular cylinder in a freestream flow, and verified and validated using published simulation data and experimental data. The sensitivity analysis study is extended to shape optimization, wherein the shape of a circular cylinder is altered and the sensitivities of the pressure drag coefficient to the changes in the cylinder shape are determined and verified. / Master of Science

Spacecraft structures rely on honeycomb panels to provide a light weight means to support the vehicle. Honeycomb panels can carry significant load but are most vulnerable to structural failure at their joints where panels connect. This research shows that predicting sandwich panel joint capability using finite element analysis (FEA) is possible. This allows for the potential elimination of coupon testing early in a spacecraft design program to determine joint capability. Linear finite element analysis (FEA) in NX Nastran was used to show that adhesive failure can be predicted with reasonable accuracy by including a fillet model on the edge of the fitting. Predicting the ultimate failure of a joint using linear FEA requires that engineering judgment be used to determine whether failure of certain bonds in a fitting will lead to ultimate joint failure or if other bonds will continue to carry the joint's load. The linear FEA model is also able to predict when the initiation of core failure will begin. This has the limitation that the joint will still be able to continue to carry significantly more load prior to joint ultimate failure even after the core has begun to buckle. A nonlinear analysis is performed using modified Riks' method in Abaqus FEA to show that this failure mode is predictable. The modified Riks' analysis showed that nonlinear post-buckling analysis of a honeycomb coupon can predict ultimate core failure with good accuracy. This solution requires a very high quality mesh in order to continue to run after buckling has begun and requires imperfections based on linear buckling mode shapes and thickness tolerance on the honeycomb core to be applied. / Master of Science

Finite Element Analysis (FEA) is proposed to simulate the connection response of practical in-service conditions and test the performance of Oil Country Tubular Goods (OCTG) premium threaded connections. A plug-in is developed in Abaqus/CAE for creating the 360-degree full 3D parametric finite element model with helical threads as an effective design and analysis tool. All size, position and material data of the model are parameterized. The premium connection plug-in accepts input from the Graphical User Interface (GUI) for further modification. Each premium connection component is programed as a collection of single-purpose independent functions organized as an independent module in order to allow users to modify specific function behavior conveniently. A main program is designed as an Abaqus kernel plug-in to achieve all the functions by calling these independent functions, making the plug-in flexible. Each single script file is not too long to jeopardize readability. The GUI of the plug-in is designed with proper layout arrangement and illustrations to make the plug-in user-friendly and easy to use. The premium connection FE model is used in a virtual test to validate the model against the ISO 13679 test protocol and is used to develop the seal metrics for points on the ISO 13679 sealability envelope. The plug-in can be used to develop and evaluate the design envelope of the premium connection. / Master of Science