Finite element analysis (FEA) software is used to obtain linear and non-linear solutions to one, two, and three-dimensional (3-D) geometric problems that will see a particular load and constraint case when put into service. Parametric FEA models are commonly used in iterative design processes in order to obtain an optimum model given a set of loads, constraints, objectives, and design parameters to vary. In some instances it is desirable for a designer to obtain some intuition about how changes in design parameters can affect the FEA solution of interest, before simply sending the model through the optimization loop. This could be accomplished by running the FEA on the parametric model for a set of part family members, but this can be very timeconsuming and only gives snapshots of the models real behavior. The purpose of this thesis is to investigate a method of visualizing the FEA solution of the parametric model as design parameters are changed in real-time by approximating the FEA solution using surrogate modeling methods. The tools this research will utilize are parametric FEA modeling, surrogate modeling methods, and visualization methods. A parametric FEA model can be developed that includes mesh morphing algorithms that allow the mesh to change parametrically along with the model geometry. This allows the surrogate models assigned to each individual node to use the nodal solution of multiple finite element analyses as regression points to approximate the FEA solution. The surrogate models can then be mapped to their respective geometric locations in real-time. Solution contours display the results of the FEA calculations and are updated in real-time as the parameters of the design model change.

Measuring the effects of morphology on performance, and performance on fitness, is necessary to gain a full picture of selection, adaptation, ecology, and evolution. The performance of an organism's feeding apparatus, of which teeth are an integral part, has obvious implications for its fitness and survival. Extant shark teeth encompass a wide variety of shapes, and are often ascribed qualitative functions without any biomechanical testing, employing terminology such as gripping, piercing, crushing, cutting, or tearing. Additionally, teeth also comprise the vast majority of the fossil record of sharks. Therefore to understand the evolution of the shark feeding mechanism, we must understand the contribution of all parts of the feeding apparatus, including the teeth. Performance testing of extant and extinct shark teeth, nanoindentation of shark teeth, finite element analysis of tooth morphology, and phylogenetically informed analyses of shark tooth morphology and ecology were employed to elucidate the relationship between performance, ecology, and evolution. Performance testing of teeth in puncture and draw revealed few morphological patterns, indicating that most morphologies are functionally equivalent. Finite element modeling of teeth in puncture, draw, and holding showed that shark teeth are structurally strong and unlikely to fail during feeding. Evolutionary analyses of tooth shape and ecology showed no relationship between morphology, habitat, and diet. These results have significant implications for the shark paleontology, where the shapes of shark teeth are used to make assumptions about ecology and evolution.

The objective of this study was to design and build an after-market suspension for the rear wheels of a manual wheelchair. Suspension for wheelchairs is important because it has been reported that the International Organization for Standards' requirements for vibration loads on wheelchair users (ISO 2631-1), are not meet by today's standard wheelchairs. Today's wheelchairs need to be able to absorb everyday shock loads, thereby minimizing the energy transmitted to the user. The chosen design is based around the concept of adding shock reduction material between the hub of the wheel, and the axel bolt that connects the wheel to the frame of the chair. The approach taken was to design a suspension system that resides between an oversized wheel bearing, and the axle. To do this, ball-race bearings with an inner diameter of 4" were chosen, and polyurethane rubber was used as the shock absorbing material. Pro-Mechanica, a finite element analysis program, was used to analyze the suspension system. Since the most common camber/tilt for wheelchair wheels is three degrees from the vertical, the anticipated loads were applied to the wheel at this angle. A prototype of the suspension system was constructed to verify that the design would work, but no tests were performed on it. This analysis showed that the suspension system should not fail when subjected to 10 times the static load. This load was considered large enough to encompass the forces that a wheelchair chair wheel is typically subjected to. There is room for further work in the area of weight reduction, and in the use of the suspension system on steeper wheel cambers.

The scheduled replacement of the 8th North Bridge, in Salt Lake City, UT, presented a unique opportunity to test a pre-cast concrete deck, steel girder bridge. A live-load test was performed under the directions of Bridge Diagnostic Inc (BDI) and Utah State University. Six different load paths were chosen to be tested. The recorded data was used to calibrate a finite-element model of this superstructure, which was created using solid, shell, and frame elements. A comparison between the measured and finite-element response was performed and it was determined that the finite-element model replicated the measured results within 3.5% of the actual values. This model was later used to obtain theoretical live-load distribution factors, which were compared with the AASHTO LRFD Specifications estimations. The analysis was performed for the actual condition of the bridge and the original case of the bridge, which included sidewalks on both sides. The comparison showed that the code over predicted the behavior of the actual structure by 10%. For the original case, the code's estimation differed by as much as 45% of the theoretical values. Another opportunity was presented to test the behavior of a cast-in-place concrete box girder bridge in Joaquin County, CA. The Walnut Grove Bridge was tested by BDI at the request of Utah State University. The test was performed with six different load paths and the recorded data was used to calibrate a finite-element model of the structure. The bridge was modeled using shell elements and the supports were modeled using solid elements. The model was shown to replicate the actual behavior of the bridge to within 3% of the measured values. The calibrated model was then used to calculate the theoretical live-load distribution factors, which allowed a comparison of the results with the AASHTOO LRFD Specifications equations. This analysis was performed for the real conditions of the bridge and a second case where intermediate diaphragms were not included. It was determined that the code's equations estimated the behavior of the interior girder more accurately for the second model (within 10%) than the real model of the bridge (within 20%).

Underground coal mine explosions pose a significant threat to infrastructure such as mine seals and refuge alternative chambers. After a mine seal failed in the Sago mine disaster, which took the life of 12 miners, design requirements were reexamined and improved. However, most research being completed on the analysis of mine structures during an explosive event focuses solely on peak pressure values, while ignoring the impact of pressure duration. This study investigates the impact pressure duration, waveform shape, and impulse have on structural displacement, while also exploring what pressures and duration can be expected during a mine explosion. Additionally, the use of high explosives to simulation conditions experienced during a mine explosion is examined. Results from this study are produced through experimental testing using a scaled shock tube and theoretical studies using finite element analysis.

In reinforced concrete (RC) structures, adequate bond between the reinforcement and concrete is required to achieve a true composite system, in which reinforcing steel carries tensile stress, once concrete cracks, and concrete and reinforcing steel carry compression. Determining bond strength and required development length for shear transfer between concrete and reinforcement is an ongoing research subject in the field of reinforced concrete with advances in the concrete and reinforcement materials requiring continuous experimental efforts. Finite element analysis (FEA) provides opportunities to explore structural behavior of RC structures beyond the limitations of experimental testing. However, there is a paucity of research studies employing FEA to investigate the reinforcement-concrete bond-zone behavior and related failure mechanism. Instead, most FEA-based research associated with RC bond has centered on developing a bond (or interface) constitutive model for use in FEA that, by itself, can characterize bond-zone behavior, typically represented by the bond stress-slip displacement relationship. This class of bond models is useful for simulating the global behavior of RC structures but is limited in its ability to simulate local bond resistance for geometries and material properties that differ substantially from those used to calibrate the model. To fill this gap in research, this study proposes a finite element (FE) modeling approach that can simulate local bond-zone behavior in reinforced concrete. The proposed FE model is developed in a physics-based way such that it represents the detailed geometry of the bond-zone, including ribs on the deformed reinforcement, and force transfer mechanisms at the concrete-reinforcement interface. The explicit representation of the bond-zone enables simulation of the local concrete compression due to bearing of ribs against concrete and subsequent hoop tension in the concrete. This causes bond failure either due to local concrete crushing (leading to reinforcement pullout) or global concrete splitting. Accordingly, special attention is given to the selection and calibration of a concrete model to reproduce robust nonlinear response. The power of the proposed modeling approach is its ability to predict bond failure and damage patterns, based only on the physical and material properties of the bond area. Thus, the successful implementation and application of this approach enables the use of FEA simulation to support the development of new design specifications for bond zones that include new and improved materials.

Impact strength is one of the most important structural properties for a designer to consider, but is often the most difficult to quantify or measure. A constant concern in the field of composites is the effect of foreign object impact damage because it is often undetectable by visual inspection. An impact can create interlaminar damage that often results in severe reductions in strength and instability of the structure. The main objective of this study is to determine the effectiveness of a damage arrestment device (DAD) on the mechanical behavior of composite sandwiches, following a low-velocity impact. A 7.56-lbf crosshead dropped from a height of 37.5-inches was considered for the low-velocity impact testing. In this study, the experimental and numerical analysis of composite sandwiches were investigated, which included static 4-point bend and vibration testing. Composite sandwiches were constructed utilizing four-plies of Advanced Composites Group LTM45EL/CF1803 bi-directional woven carbon fiber face sheets with a General Plastics Last-A-Foam FR-6710 rigid polyurethane core. Specimens were cured in an autoclave, using the manufacturer’s specified curing cycle. In addition to the experimental and numerical analysis of composite sandwiches, developing and building a data acquisition (DAQ) system for the Dynatup 8250 drop weight impact tester was accomplished. Utilizing National Instruments signal conditioning hardware, in conjunction with LabView and MATLAB, complete testing software was developed and built to provide full data acquisition for an impact test. The testing hardware and software provide complete force vs. time history and crosshead acceleration of the impact event, as well as provide instantaneous impact velocity of the projectile. The testing hardware, software, and procedures were developed and built in the Aerospace Structures/Composites laboratory at Cal Poly for approximately 15% of the cost from the manufacturer. In the first study, static 4-point bend testing was investigated to determine the residual flexural strength of composite sandwich beams following a low-velocity impact. Four different specimen cases were investigated in the 4-point bend test, with and without being impacted: first a control beam with no delamination or DAD, second a control beam with a centrally located 1-inch long initial delamination, third a DAD key beam with two transverse DADs centrally located 1-inch apart, and finally a DAD key beam with a centrally located initial delamination between two transverse DADs. The specimens used followed the ASTM D6272 standard test method. The specimens were 1-inch wide by 11-inch long beams. The experimental results showed that the presence of DAD keys significantly improved both the residual stiffness and ultimate strength of a composite sandwich structure that had been damaged under low-velocity impact loading, even with the presence of an initial face-core delamination. In the second study, vibration testing was investigated as a means to detect a delamination in the structure and the effect of impact damage on the vibrational characteristics, such as damping, on composite sandwich plates. Four different specimen cases were investigated in the vibration test, both with and without being impacted: first a control plate with no delamination or DAD, second three control plates with varying 1-inch initial delamination locations at the 1st, 2nd, and 3rd bending-mode nodes, third a DAD key plate with one DAD running the entire length longitudinally along the center of the plate, and finally three DAD key plates with one DAD running the entire length longitudinally along the center of the plate and varying 1-inch delamination locations at the 1st, 2nd, and 3rd bending mode-nodes. The response accelerometer location was varied at 1-inch increments along the length of the plate. From the experimental results, it was determined that varying the location of the accelerometer had a significant effect on the detection of face-core delamination in a composite sandwich structure. Additionally, it was shown that damping characteristics significantly degraded in control case plates after a low-velocity impact, but they were better retained when a DAD key was added to the structure. Numerical analysis utilizing the finite element method (FEM) was employed to validate experimental testing, as well as provide a means to examine the stress distribution and impact absorption of the structure. The impact event was modeled utilizing the LS-Dyna explicit FE solver, which generated complete force vs. time history of the impact event. Static 4-point bending and vibration analysis were solved utilizing the LS-Dyna implicit solver. Finally a damaged mesh was obtained from the explicit impact solution and subjected to subsequent static 4-point bending and vibration analysis to numerically determine the residual mechanical behavior after impact. All cases showed good agreement between the numerical, analytical, and experimental results.

A primary motivation for this research comes from the need to improve the ability of polymer-matrix composites to withstand lightning strikes. In particular, we are concerned with lightning strike damage in composite wind turbine blades. The direct effects of lightning strike on polymer-matrix composites often include rapid temperature rise, melting or burning at the lightning attachment points, and mechanical damage due to lightning-induced magnetic force and acoustic shock wave. The lightning strike damage accumulation problem is essentially multiphysic. The lightning plasma channel discharges an electric current up to 200 kA, inducing a severe heat flux at the surface of the composite structure, as well as generating Joule heating through the composite structure. The resulting electro-thermo-mechanical response of the composite structure may include matrix degradation and decomposition, delamination, and fiber breakage and sublimation, thus leading to catastrophic failure. The existing studies related to the lightning strike damage in composites ignored the lightning channel radius expansion during the initial lightning discharge and lacked adequate treatment of material phase transitions. These assumptions significantly simplify the mathematical treatment of the problem and affect the predictive capabilities of the models. Another common feature of these limited studies is that they all focused on carbon-fiber-reinforced polymer-matrix (CFRP) composites, which are electrically conductive. In the present thesis, the thermal responses and thermal ablations in a non-conductive glass-fiber-reinforced polymer-matrix (GFRP) composite wind turbine blade and in a conductive CFRP composite wind turbine blade are studied, respectively. In the case of non-conductive GFRP composite wind turbine blade, prior to the thermal response and thermal ablation analysis, a finite element analysis is performed to calculate the electric field due to lightning stepped leader to estimate the dielectric breakdown of the non-conductive composite wind turbine blade. The estimation of dielectric breakdown is used to determine whether Joule heating needs to be included in the problem formulation. To predict the thermal response and thermal ablation in the composite structure due to lightning strike, a physics-based model describing surface interaction between the lightning channel and the composite structure has been developed. The model consists of: (i) spatial and temporal evolution of the lightning channel as a function of the electric current waveform; (ii) temporary and spatially non-uniform heat flux and current density (in the case of electrically conductive CFRP composite or if dielectric breakdown occurs in the case of non-conductive GFRP composite) generated at the composite structure; and (iii) nonlinear transient heat transfer problem formulation for layered anisotropic composites that includes the moving boundary of the expanding lightning channel and the phase transition moving boundary associated with instantaneous material removal due to sublimation. The model has been employed to investigate the thermal responses and thermal ablations in a GFRP composite laminated panel used in a Sandia 100-meter all-glass baseline wind turbine blade (SNL 100-00) and a typical CFRP composite laminated panel subjected to lightning strike. The temperature-dependent directional material properties for both the GFRP and CFRP composites have been determined in this thesis using a micromechanics approach based on the experimental data for fibers and resin. An integrated Matlab-ABAQUS numerical procedure features the aforementioned aspects (i), (ii), and (iii) of the developed model. The obtained results include the evolution of temperature fields in the composite laminated panel and the progressive shape change of the composite laminated panel due to thermal ablation. The predictions of thermal ablation in the CFRP composite laminated panel are validated by reported experimental results.

Aseptic loosening due to wear-induced osteolysis remains a leading cause of failure in total hip arthroplasty (THA), particularly in revision cases beyond the second decade of use. Historically, there have been large amounts of variability of wear within individual THA patient cohorts. Evidence indicates that femoral head damage can be a cause of this variability. While femoral head damage as a result of third body particles and subluxation and dislocation events has been well documented, direct quantifiable linkage between femoral head damage and wear acceleration remains to be established. Due to large ranges of observed retrieval damage, wear testing protocols for simulating third body and other damage effects have been subject to a wide range of variability, making it difficult to know where the clinical reality lies. To study the effect of retrieval femoral head damage on total hip implant wear, a damage-feature-based finite element (FE) formulation which allowed for wear prediction due to individual damage features developed. A multi-scale imaging procedure was also developed to globally map and quantify micron-level damage features appearing on retrieval femoral heads. This allowed for wear simulations of damage patterns observed on specific retrieval femoral heads. Retrieval damage was shown to be highly variable among patients, and capable of producing up to order-of-magnitude wear increases when compared to undamaged head wear rates. Damage following dislocation and subsequent closed reduction maneuvers was particularly detrimental, with average wear rate increases equal to half an order of magnitude. These data were used to develop wear testing protocols for simulating clinically-occurring third body and other damage effects.

The use and implementation of adhesive joints for space structures is necessary for incorporating fiber-reinforced composite materials. Correct modeling and design of cylindrical adhesive joints can increase the dimensional stability of space structures. The few analytical models for cylindrical adhesive joints do not fully describe the displacement or stress field of the joint. A two-dimensional axisymmetric finite element model for the design and analysis of adhesive joints was developed. The model was developed solely for the analysis of cylindrical adhesive joints, but the energy techniques used to develop the model can be applied to other types of joints as well. A numerical program was written to solve the system of equations [K]{d}={R} for the unknown displacements {d}. The displacements found from the program are used to design cylindrical adhesive joints based on dimensional stability. Stresses were calculated from the displacements for comparison with analytical models. The cylindrical joints were assumed to remain within the linear elastic region and no failure criteria was taken into account. The design process for cylindrical joints was developed based on dimensional stability. The nodal displacements found from the finite element model were used in the optimization of geometric parameters of cylindrical joints. The stacking sequence of the composite, the bond length, and the bond thickness were found to have the greatest impact on dimensional stability. Other factors that were found to further reduce the maximum displacements are the implementation of 0° and 90° laminas, the isotropic cylinder thickness, tapering of the isotropic cylinder, and the inside radius of the cylindrical joint. This axisymmetric finite element model is beneficial in that a cylindrical joint can be designed before any testing is performed. The results and cases in this thesis are generalized in order to show how the design process works. The model can be used in conjunction with design requirements for a specific joint to reduce the maximum displacements below any specified operating requirements. The joint is dimensionally stable if the overall displacements meet specific design requirements.