Browsing by Author "Heyliger, Paul, committee member"
Now showing 1 - 20 of 34
Results Per Page
Sort Options
Item Open Access A Haversian bone model of fracture healing in a simulated microgravity environment(Colorado State University. Libraries, 2015) Gadomski, Benjamin C., author; Puttlitz, Christian M., advisor; Browning, Raymond, committee member; Donahue, Tammy, committee member; Heyliger, Paul, committee memberGround-based models of weightlessness and microgravity have provided valuable insights into how dynamic physiological systems adapt or react to reduced loading. Almost all of these models have used rodent hindlimb unloading as the means to simulate microgravity on isolated physiological systems. Unfortunately, results derived from rodent studies are significantly limited when one tries to translate them to the human condition due to significant anatomical and physiological differences between the two species. Therefore, it is clear that a novel animal model of ground-based weightlessness that is directly translatable to the human condition must be developed in order for substantial progress to be made in the knowledge of how microgravity affects fracture healing. In light of this, four specific aims are proposed: (1) develop a ground-based, ovine model of skeletal unloading in order to simulate weightlessness, (2) interrogate the effects of the simulated microgravity environment on bone fracture healing using this large animal model, (3) develop a computational model of weightbearing in ovine bone under different experimental conditions in order to characterize the loads experienced by the fracture site, and (4) develop countermeasures that enhance bone fracture healing in the presence of simulated microgravity. Successful completion of this project will substantially elevate the understanding of how fracture site loading affects the subsequent healing cascade in the presence of microgravity and will form the foundation for designing future rehabilitation protocols to facilitate bone healing during long-duration spaceflight.Item Open Access A thermoplastic matrix continuous fiber reinforced composite impregnation method by direct polymer extrusion(Colorado State University. Libraries, 2018) Hedin, Kevin M., author; Radford, Donald W., advisor; Ma, Kaka, committee member; Heyliger, Paul, committee memberDuring component design, continuous fiber reinforced composite material systems are often chosen largely based on their structural efficiency. Their mechanical properties, such as specific strength and specific stiffness, are often cited as significant advantages over the use of other materials. However, composite component production often lacks the capability to provide the local variation necessary to ensure that 1) the reinforcing fibers are best aligned with anticipated loads, and 2) the ideal matrix composition and fiber volume fraction are found throughout the composite part. In practice, these limitations result in composite components that do not demonstrate the maximum possible efficiencies inherent to the fiber-reinforced composite material system. To further increase the flexibility of polymer matrix continuous fiber reinforced composites manufacturing methods, a new thermoplastic impregnation method was developed. This proposed method adds a thermoplastic matrix, which has previously been proven to allow significant variation of local fiber orientation, to the reinforcing fiber just prior to the consolidation of the composite. The increased independence of matrix and fiber addition should allow the local variation of volume and composition of the added matrix, while using less and simpler hardware than previous, similar efforts. In this work, the quality of material deposited from the proposed process is evaluated. The maximum possible quality of the proposed method and also that of a similar process that uses a commercially available material system were determined, primarily using short beam shear (SBS) testing. The material system of both methods consisted of E-glass continuous fiber reinforcement with a PETG matrix. It was found that both manufacturing processes are capable of producing samples with an SBS strength of approximately 53 MPa, and it was concluded that the proposed process has the capability to deposit material of comparable quality to that produced by the baseline method. Subsequent thermal analysis, fiber volume fraction/void content measurement, and metallographic imaging were conducted to investigate the effects of using two different PETG compositions on the SBS strength of composite material produced by the proposed process. It was found that, while using the proposed process, the PETG matrix with a lower glass transition temperature allowed better consolidation of the resulting composite part, ultimately increasing SBS strength. Each process parameter used in the proposed process was evaluated for the practical significance of its effects on SBS strength, which facilitated 1) an understanding of the underlying mechanisms of the process, and 2) a tenable simplification of the process that should reduce operating costs and also demonstrates its robustness via insensitivity to many of the possible process variations. Finally, it was established that the material inputs to the proposed process are relatively inexpensive: Using PETG and continuous E-glass fiber in the proposed process reduces material input cost by at least 52% compared to using commingled PETG and E-glass fibers in the baseline process, on a $/kg basis.Item Open Access Additive manufacture of dissolvable tooling for autoclave processing of fiber reinforced polymer composites(Colorado State University. Libraries, 2022) Morris, Isaac, author; Radford, Donald, advisor; Yourdkhani, Mostafa, committee member; Heyliger, Paul, committee memberAutoclave processing of advanced fiber reinforced polymer composites (AFRPC) uses applied heat and pressure to yield high quality composite components. Geometrically accurate and thermally stable molds or tools are used to maintain the part form until the part cures and rigidizes. For high-volume production runs, molds may be made from materials such as metals, ceramics, or AFRPCs. However, tooling made from these materials can be costly to manufacture and are not suitable for low volume production runs. This is especially true for complex geometries in trapped tooling situations where the cured composite shape prevents tool separation. In this situation, composite manufacturers rely on sacrificial washout tooling materials that are machined or cast to shape to create the tool. However, these sacrificial materials still come with significant challenges. For example, the surfaces of these tools are often porous and require sealing, and their washout can result in corrosive waste that makes disposal challenging. Additionally, these tools are brittle and monolithic in nature, making them fragile to handle and slow to heat up during cure. An alternative may be to use high temperature, dissolvable thermoplastic materials in melt extrusion additive manufacturing to create complex washout tooling. However, there is a lack of information regarding the types of soluble materials and the structural configurations that make this type of tooling successful in autoclave use. To begin to address this, samples made from several materials, and one insoluble model material, were processed in stepwise fashion at increasing autoclave processing temperatures to evaluate the impacts of material and structure on autoclave robustness. Then, mid-sized composite specimens were produced on 3D-printed tooling that evaluated the interaction between the composite and the tool, including surface quality and deformation. Finally, a trapped tooling geometry was used to manufacture several composites at processing conditions of 157°C at 414kPa, well above the use temperature of the tested materials. These trials focused on reducing deformation by adjusting the tool wall thickness and vacuum bagging configuration. It was shown that 3D-printed dissolvable tooling can be used as an alternative to traditional washout tooling for autoclave processing. The materials Stratasys ST-130 and Infinite Material Solutions AquaSys 180 were used to manufacture tools that were processed at autoclave conditions of 121°C at 345kPa with minimal deformation. Surface quality was also found to be acceptable without machining or sealing, eliminating this step from the production of traditional washout tools. Finally, a modified tool design and vacuum bagging technique were demonstrated that significantly reduced the deformation of tooling at processing temperatures that significantly exceed the use temperature of the material.Item Open Access Additive manufacturing of an intervertebral disc repair patch to treat spinal herniation(Colorado State University. Libraries, 2021) Page, Mitchell Ian, author; Puttlitz, Christian, advisor; Heyliger, Paul, committee member; Susan, James, committee member; Kirk, McGilvray, committee memberChronic low back pain is ubiquitous throughout society. The consequences of this disease are extensive and lead to physical, mental, and financial suffering in the affected population. Herniation of the intervertebral disc (IVD) is the primary cause of chronic low back pain due to the essential mechanical role of the IVD in the spinal column. Degenerative changes to the IVD tissues, in particular the annulus fibrosus (AF), lead to a pronounced vulnerability to herniation. Although numerous treatments for intervertebral disc herniation currently exist, these treatments are typically palliative and prone to hernia recurrence. Accordingly, there is a distinct need for an IVD hernia therapy that can provide long-term pain relief and recovery of spinal function. One novel strategy to repair the intervertebral disc is to tissue-engineer a construct that facilitates regeneration of the healthy and functional IVD tissue. Advances in additive manufacturing technology offer the fabrication of complex tissue-engineered structures that augment biological content and biocompatible materials. Therefore, this work sought to engineer an additive manufactured repair patch for IVD herniation towards an improved treatment for chronic low back pain. Specifically, the aims of this work were to leverage experimental and computational methods to: (1) to characterize the mechanics of additive manufactured angle-ply scaffolds, (2) evaluate the tissue response of cell-laden scaffolds cultured with dynamic biaxial mechanical stimuli, and (3) to design and implement an annulus fibrosus repair patch. The mechanics of additive manufactured scaffolds for AF repair were experimentally characterized in a physiologically-relevant, biaxial loading modality. To assess sensitivity of the scaffold mechanics to additive manufacturing parameters, a broad scope of scaffold designs were evaluated with a parameterized finite element model. A custom incubator was developed, cell-laden scaffolds were cultured with a prescribed, multi-axial mechanical loading protocol, and ECM production within the scaffold was evaluated. A finite element model was developed to aid in understanding the relationship between global scaffold loading and the local, inhomogeneous cellular micromechanical environment within the scaffold. The developed TE material was translated into an implant and was implemented in a large animal model. The efficacy of the AF repair strategy was also evaluated in finite element model of the human lumbar spine. This work formed a multi-scale approach to consolidate biological and mechanical efficacy of a novel AF repair strategy. Ultimately, this approach may facilitate regeneration of the AF and represent a revolutionary treatment for chronic low back pain.Item Open Access An analysis of domain decomposition methods using deal.II(Colorado State University. Libraries, 2021) Rigsby, Christina, author; Tavener, Simon, advisor; Bangerth, Wolfgang, committee member; Heyliger, Paul, committee member; Liu, Jiangguo, committee memberIterative solvers have attracted significant attention since the mid-20th century as the computational problems of interest have grown to a size beyond which direct methods are viable. Projection methods, and the two classical iterative schemes, Jacobi and Gauss-Seidel, provide a framework in which many other methods may be understood. Parallel methods or Jacobi-like methods are particularly attractive as Moore's Law and computer architectures transition towards multiple cores on a chip. We implement and explore two such methods, the multiplicative and restricted additive Schwarz algorithms for overlapping domain decomposition. We implement these in deal.II software, which is written in C++ and uses the finite element method. Finally, we point out areas for potential improvement in the implementation and present a possible extension of this work to an agent-based modeling prototype currently being developed by the Air Force Research Laboratory's Autonomy Capability Team (ACT3).Item Open Access Assessment of the effects of ligamentous injury in the human cervical spine(Colorado State University. Libraries, 2012) Leahy, Patrick Devin, author; Puttlitz, Christian, advisor; Heyliger, Paul, committee member; Sakurai, Hiroshi, committee member; Santoni, Brandon, committee memberLigamentous support is critical to constraining motion of the cervical spine. Injuries to the ligamentous structure can allow hypermobility of the spine, which may cause deleterious pressures to be applied to the enveloped neural tissues. These injuries are a common result of head trauma and automobile accidents, particularly those involving whiplash-provoking impacts. The injuries are typically relegated to the facet capsule (FC) and anterior longitudinal (ALL) ligaments following cervical hyperextension trauma, or the flaval (LF) and interspinous (ISL) ligaments following hyperflexion. Impacts sustained with the head turned typically injure the alar ligament. The biomechanical sequelae resulting from each of these specific injuries are currently ill-defined, confounding the treatment process. Furthermore, clinical diagnosis of ligamentous injuries is often accomplished by measuring the range of motion (ROM) of the vertebrae, where current methods have difficulty differentiating between each type of ligamentous injury. Pursuant to enhancing treatment and diagnosis of ligamentous injuries, a finite element (FE) model of the intact human full-cervical (C0-C7) spine was generated from computed tomography (CT) scans of cadaveric human spines. The model enables the quantification of ROM, stresses, and strains, and can be modified to reflect ligamentous injury. In order to validate the model, six human, cadaveric, full-cervical spines were tested under pure ±1.5 Nm moment loadings in the axial rotation, lateral bending, flexion, and extension directions. ROM for each vertebra, facet contact pressures, and cortical strains were experimentally measured. To evaluate injured ligament mechanical properties, a novel methodology was developed where seven alar, fourteen ALL, and twelve LF cadaveric bone-ligament-bone preparations were subjected to a partial-injury inducing, high-speed (50 mm/s) tensile loading. Post-injury stiffnesses and toe region lengths were compared to the pre-injury state for these specimens. These experimental data were incorporated into the FE model to analyze the kinematic and kinetic effects of partial ligamentous injury. For comparison, the model was also adapted to reflect fully injured (transected) ligaments. Injuries simulated at the C5-C6 level included: 1) partial FC injury, 2) full FC injury, 3) partial FC and ALL injury, 4) full FC and ALL injury, 5) partial LF and full ISL jury, 6) full LF and ISL injury, 7) partial FC, ALL, LF, and full ISL injury, and 8) full FC, ALL, LF, and ISL injury. The model was also modified to replicate injury to the right alar ligament. Five cadaveric cervical spines were tested under pure moment conditions with scalpel-sectioning of these ligaments for validation of the full-injury models. Comparisons between the intact and various injury cases were made to determine the biomechanical alterations experienced by the cervical spine due to the specific ligamentous injuries. Variances in ROM and potential impingement on the neural tissues were focused upon. The overarching goals of the study were to identify a unique kinematic response for each specific ligamentous injury to allow for more accurate clinical diagnosis, and to enhance the understanding of the post-injury biomechanical sequelae.Item Open Access Bio-inspired design for engineering applications: empirical and finite element studies of biomechanically adapted porous bone architectures(Colorado State University. Libraries, 2020) Aguirre, Trevor Gabriel, author; Donahue, Seth W., advisor; Ma, Kaka, committee member; Heyliger, Paul, committee member; Simske, Steven, committee memberTrabecular bone is a porous, lightweight material structure found in the bones of mammals, birds, and reptiles. Trabecular bone continually remodels itself to maintain lightweight, mechanical competence, and to repair accumulated damage. The remodeling process can adjust trabecular bone architecture to meet the changing mechanical demands of a bone due to changes in physical activity such as running, walking, etc. It has previously been suggested that bone adapted to extreme mechanical environments, with unique trabecular architectures, could have implications for various bioinspired engineering applications. The present study investigated porous bone architecture for two examples of extreme mechanical loading. Dinosaurs were exceptionally large animals whose body mass placed massive gravitational loads on their skeleton. Previous studies investigated dinosaurian bone strength and biomechanics, but the relationships between dinosaurian trabecular bone architecture and mechanical behavior has not been studied. In this study, trabecular bone samples from the distal femur and proximal tibia of dinosaurs ranging in body mass from 23-8,000 kg were investigated. The trabecular architecture was quantified from micro-computed tomography scans and allometric scaling relationships were used to determine how the trabecular bone architectural indices changed with body mass. Trabecular bone mechanical behavior was investigated by finite element modeling. It was found that dinosaurian trabecular bone volume fraction is positively correlated with body mass like what is observed for extant mammalian species, while trabecular spacing, number, and connectivity density in dinosaurs is negatively correlated with body mass, exhibiting opposite behavior from extant mammals. Furthermore, it was found that trabecular bone apparent modulus is positively correlated with body mass in dinosaurian species, while no correlation was observed for mammalian species. Additionally, trabecular bone tensile and compressive principal strains were not correlated with body mass in mammalian or dinosaurian species. Trabecular bone apparent modulus was positively correlated with trabecular spacing in mammals and positively correlated with connectivity density in dinosaurs, but these differential architectural effects on trabecular bone apparent modulus limit average trabecular bone tissue strains to below 3,000 microstrain for estimated high levels of physiological loading in both mammals and dinosaurs. Rocky Mountain bighorn sheep rams (Ovis canadensis canadensis) routinely conduct intraspecific combat where high energy cranial impacts are experienced. Previous studies have estimated cranial impact forces up to 3,400 N and yet the rams observationally experience no long-term damage. Prior finite element studies of bighorn sheep ramming have shown that the horn reduces brain cavity translational accelerations and the bony horncore stores 3x more strain energy than the horn during impact. These previous findings have yet to be applied to applications where impact force reduction is needed, such as helmets and athletic footwear. In this study, the velar architecture was mimicked and tested to determine suitability as novel material architecture for running shoe midsoles. It was found that velar bone mimics reduce impact force (p < 0.001) and higher energy storage during impact (p < 0.001) and compression (p < 0.001) as compared to traditional midsole architectures. Furthermore, a quadratic relationship (p < 0.001) was discovered between impact force and stiffness in the velar bone mimics. These findings have implications for the design of novel material architectures with optimal stiffness for minimizing impact force.Item Open Access Constitutive modeling of the biaxial mechanics of brain white matter(Colorado State University. Libraries, 2016) Labus, Kevin M., author; Puttlitz, Christian M., advisor; Donahue, Seth, committee member; Heyliger, Paul, committee member; James, Susan, committee memberIt is important to characterize the mechanical behavior of brain tissue to aid in the computational models used for simulated neurosurgery. Due to its anisotropy, it is of particular interest to develop constitutive models of white matter based on experimental data in order to define the material properties in computational models. White matter has been shown to exhibit anisotropic, hyperelastic, and viscoelastic properties. The majority of studies have focused on the shear or compressive properties, while few have tested the tensile properties of the brain. Brain tissue has not previously been tested in a multi-axial loading state, even though in vivo brain tissue is in a constant multi-axial stress state due to fluid pressure, and data from uniaxial experiments do not sufficiently describe multi-axial stresses. The main objective of this project was to characterize the biaxial tensile behavior of brain white matter via experimentation and constitutive modeling. A biaxial experiment was developed specifically for testing brain tissue. Experiments were performed at a quasi-static loading rate, and an Ogden anisotropic hyperelastic model was derived to fit the data. A structural analysis was performed on biaxially tested specimens to relate the structure to the mechanical behavior. The axonal orientation and distribution were measured via histology, and the axon area fraction was measured via transmission electron microscopy. The measured structural parameters were incorporated into the constitutive model. A probabilistic analysis was performed to compare the uncertainty in the stress predictions between models with and without structural parameters. Finally, dynamic biaxial experiments were performed to characterize the anisotropic viscoelastic properties of white matter. Biaxial stress-relaxation experiments were conducted to determine the appropriate form of a viscoelastic model. It was found that the data were accurately modeled by a quasi-linear viscoelastic formulation with an isotropic reduced relaxation tensor and an instantaneous elastic stress defined by an anisotropic Ogden model. Model fits to the stress-relaxation experiments were able to accurately predict the results of dynamic cyclic experiments. The resulting constitutive models from this project build upon previous models of brain white matter mechanics to include biaxial interactions and structural relations, thus improving computational model predictions.Item Open Access Determination of seismic performance factors for cross laminated timber shear wall system based on FEMA P695 methodology(Colorado State University. Libraries, 2018) Amini, Mohammad Omar, author; van de Lindt, John W., advisor; Mahmoud, Hussam, committee member; Heyliger, Paul, committee member; Senior, Bolivar, committee memberCross Laminated Timber (CLT) was initially introduced in Europe and has recently gained popularity in North America where it is seen as a sustainable alternative to steel and concrete in midrise construction. Although most CLT structures to date have been constructed in low seismic regions, recent tests have indicated that CLT based lateral force resisting systems can successfully be utilized in regions of higher seismicity. Despite the many advantages that CLT offers, the lack of a design code and systematic design procedure is one of many challenges inhibiting widespread adoption of CLT in the US. The purpose of this study was to investigate the seismic behavior of CLT based shear wall systems and determine seismic performance factors, namely, the response modification factor (R-factor), the system overstrength factor (Ω), and the deflection amplification factor (Cd), using the FEMA P695 procedure. The methodology is an iterative process that includes establishing design requirements, developing archetypes, performing a series of tests, developing and validating nonlinear models, nonlinear static and dynamic analysis, and evaluating performance; all in conjunction with a peer panel to provide input. Nine index buildings that include, single-family dwellings, multi-family dwellings, and commercial (including mixed-use) mid-rise buildings were developed. Archetypes were then extracted from these index buildings. Testing performed at the component and subassembly levels include connector tests and isolated shear wall tests. A subsequent full-scale shake table test was performed for system level demonstration. A critical aspect of this study is use of generic connectors whose properties are already addressed by a design specification to facilitate building code recognition. Test-based performance for these generic connectors is reported as part of this study to facilitate evaluation of proprietary alternatives for seismic equivalence. Connector tests were performed on angle brackets, used for attachment of the wall to the supporting element, and inter-panel connectors. These tests showed connector thickness to be important in achieving the desired ductile behavior with lesser thickness (12 gauge) being the more favorable. Quasi-static cyclic tests were conducted for a portfolio of CLT shear walls to systematically investigate the effects of various parameters. CLT demonstrated rigid behavior with energy dissipation concentrated in the connectors. Boundary constraints and gravity loading were both found to have a beneficial effect on the wall performance, i.e. higher strength and deformation capacity. Specific gravity also had a significant effect on wall behavior while CLT thickness was less influential. Higher aspect ratio panels (4:1) demonstrated lower stiffness and substantially larger deformation capacity compared to moderate aspect ratio panels (2:1). However, based on the test results there is likely a lower bound for aspect ratio (at 2:1) where it ceases to benefit deformation capacity of the wall. Multi-panel configuration comprised of high aspect ratio panels connected through vertical joint demonstrated considerably larger deformation capacity. Shake table tests showed the proposed system's potential to meet life-safety code requirements and its applicability in US seismic regions. A CLT shear wall design method was developed and refined based on the test results. Phenomenological models were used in modeling CLT shear walls. The archetypes were designed based on the proposed design method and were numerically evaluated by assessing their performance using nonlinear static and dynamic analyses. Based on the rigorous process, an R factor of 3 is proposed for the CLT shear wall systems and an R factor of 4 is proposed for the cases with high aspect ratio panels only. Results from the study will be proposed for implementation in the seismic design codes and standards in the US.Item Open Access Development of novel mechanical diagnostic techniques for early prediction of bone fracture healing outcome(Colorado State University. Libraries, 2021) Wolynski, Jakob G., author; McGilvray, Kirk, advisor; Puttlitz, Christian, advisor; Heyliger, Paul, committee member; James, Susan, committee member; Wang, Zhijie, committee memberTo view the abstract, please see the full text of the document.Item Open Access Dynamic mechanical analysis for quality evaluation of additively manufactured continuous fiber reinforced thermoplastic matrix composites subject to manufacturing defects(Colorado State University. Libraries, 2019) Rodriguez, Patrick A., author; Radford, Donald W., advisor; Ma, Kaka, committee member; Heyliger, Paul, committee memberContinuous fiber reinforced polymers (CFRP) have become integral to modern mechanical design as value-added alternatives to metallic, ceramic and neat polymeric engineering materials. Despite the advantages of CFRP, current methods of preparing laminated continuous fiber reinforced polymers are fundamentally limiting in that reinforcement is typically applied only in the plane of the mold or tool. Additionally, key operations inherent to all CFRP processing approaches require a variety of skilled labor as well as costly net-shape, hard tooling. As such, additive manufacturing has risen to the forefront of manufacturing and processing research and development in the CFRP arena. Additive manufacture of continuous fiber reinforced thermoplastics (CFRTP) exhibits the potential to relieve many of the constraints placed on the current design and manufacturing of continuous fiber reinforced structures. At present, the additive manufacture of CFRTP has been demonstrated successfully to varying extents; however, comprehensive dialogue regarding manufacturing defects and quality of the processed continuous fiber reinforced thermoplastics has been missing from the field. Considering the preliminary nature of additive manufacture of CFRTP, exemplary processed composites are typically subject to various manufacturing defects, namely excessive void content in the thermoplastic matrix. Generally, quality evaluation of processed composites in the literature is limited to test methods that are largely influenced by the properties of the continuous fiber reinforcement, and as such, defects in the thermoplastic matrix are usually less-impactful on the results and overlooked. Hardware to facilitate additive manufacturing of CFRTP was developed and continuous fiber reinforced specimens, with high fiber volume fractions (~ 50 %), were successfully processed. Early efforts at evaluating the processed specimens using defect-sensitive Short-Beam Strength (SBS) analysis exhibited limited sensitivity to void content, coupled with destructive, inelastic failure modes. As a path forward, an expanded study of the effects of void content on the processed specimens was conducted by means of Dynamic Mechanical Analysis (DMA). Utilization of DMA allows for thermomechanical (i.e. highly matrix sensitive) evaluation of the composite specimens, specifically in terms of the measured elastic storage modulus (E'), viscous loss modulus (E"), damping factor (tan δ) and the glass transition temperature (Tg) of the processed composite specimens. The results of this work have shown that DMA exhibits increased sensitivity, as compared to SBS, to the presence of void content in the additively manufactured CFRTP specimens. Within the relevant range of void content, non-destructive specimen evaluation by DMA resulted in a measured, frequency dependent, 5.5 – 5.8 % decrease in elastic storage modulus per 1 % increase in void content by volume. Additionally, quality evaluation by DMA realized a marked decrease in the maximum measured loss modulus in the additively manufactured composites, ranging from 7.0 – 8.2 % per 1 % increase in void content by volume. Effects of void content were also measured in both the damping factor and glass transition temperature, where an approximate 1.6 °C drop in Tg was recorded over the relevant range of void content. The results of this work indicate, firstly, that DMA is a superior evaluation method, as compared to SBS, in terms of sensitivity to void content in additively manufactured CFRTP. Additionally, the results of this work provide a clear expansion of the current state of the literature regarding additive manufacture of CFRTP materials in that the effects of prominent manufacturing defects have been assessed with regard to thermomechanical material performance. Furthermore, and finally, the results of this work establish a direct path forward to characterize long-term effects of manufacturing defects, by means of DMA, on the creep-recovery and stress relaxation behavior of the relevant composite material system.Item Open Access Dynamic structural analysis of ramming in bighorn sheep(Colorado State University. Libraries, 2015) Drake, Aaron Michael, author; Haut Donahue, Tammy L., advisor; Donahue, Seth W., advisor; Stansloski, Mitchell, committee member; Heyliger, Paul, committee memberConcussions are the most common traumatic brain injury and are caused by impulsive loads applied to the skull, resulting in relative motion of the brain within the brain cavity. Despite wearing helmets, athletes involved in full contact sports, such as football, are highly susceptible to concussive injuries. Short term symptoms of concussions include nausea, headache and confusion and there is evidence of more serious, long term effects from repeated concussions. Furthermore, the physical mechanisms of concussions are not well understood, making them difficult to diagnose and treat clinically. Male bighorn sheep sustain massive impact loads to the head during ramming, which is done as a means of determining hierarchy and gaining mating privileges. These large animals thrust themselves, horns first, at one another and collide violently, repeating this ritual for up to several hours until the subdominant male succumbs. After a collision, the animals are stunned momentarily but otherwise appear to suffer no ill effects, based on behavioral observations. This simple fact provided the motivation to examine the dynamic structural behavior of bighorn sheep horns and skulls. For reference, the average translational brain cavity accelerations observed during finite element model impact were found to be 111g (1091 m/s²) and impacts thought to be damaging to human brains occur at around 100g. A dynamic finite element impact model was produced using the geometry, obtained from a CT scan, of a mature male bighorn sheep’s skull and horns. Quantitative and qualitative results of the simulation were examined to determine mechanisms of energy dissipation and stress distribution during an idealized impact event. Video analysis of particularly forceful ramming sequences of wild bighorn sheep was carried out to estimate the dynamics involved with ramming. In order to investigate the relative contributions of the horn curl as well as the internal foamy bone architecture, three separate finite element models were produced. One model had one half of the horn length removed, another had the internal foam-like bone removed and these models were compared to the fully intact model to determine the structural contributions of these features during impact. Removing one half of the horn curl had the effect of increasing the peak brain cavity translational acceleration by 49%. Eliminating the internal foamy bone architecture resulted in a dramatic 442% increase in brain cavity rotational accelerations. The dynamic (vibrational) response of bighorn sheep horns and skulls was investigated using two, related methods: finite element modal analysis and experimental modal analysis. The finite element modal analysis revealed five dominant natural frequencies with values ranging from 118 to 309 Hz. Experimental modal analysis revealed several natural frequencies between 100 and 300 Hz, however, differentiating specific modes was difficult. For both vibrational analyses the dominant vibrational mode shape was side-to-side oscillations of the horn tip. This study hopes to promote and guide further research on the mechanisms of brain trauma prevention in bighorn sheep, with an emphasis on the structural and material characteristics of the horn and skull, to increase our understanding of, and ways to prevent traumatic brain injuries in humans.Item Open Access Effect of bone geometry on stress distribution patterns in the equine metacarpophalangeal joint(Colorado State University. Libraries, 2012) Easton, Katrina L., author; Kawcak, Chris, advisor; McIlwraith, Wayne, committee member; Puttlitz, Christian, committee member; James, Susan, committee member; Shelburne, Kevin, committee member; Heyliger, Paul, committee memberCatastrophic injury of the equine metacarpophalangeal joint is of major concern for both the equine practitioner and the American public. It is one of the major reasons for retirement and sometimes euthanasia of Thoroughbred racehorses. The most common type of catastrophic injury is fracture of the proximal sesamoid bones and lateral condyle of the third metacarpal bone. Many times these injuries are so disastrous that there is no possibility of fixing them. Even in the injuries that are able to be fixed, complications arising from the fracture such as support limb laminitis may ultimately lead to the demise of the horse. Therefore, prevention of these types of injuries is key. In order to decrease the incidence of injury, it is important to understand the risk factors and pathogenesis of disease that leads to them. This project was established to create a finite element model of the equine metacarpophalangeal joint in order to investigate possible risk factors, namely bone geometry, and its effect on the stress distribution pattern in the joint. The first part of the study involved in vitro experiments in order to provide a comprehensive dataset of ligament, tendon, and bone strain and pressure distribution in the joint with which to validate the finite element model. Eight forelimbs from eight different horses were tested on an MTS machine to a load equivalent to that found at the gallop. Beyond providing data for validation, the study was the first to the author's knowledge to measure surface contact pressure between the distal condyles of the third metacarpal bone and the proximal sesamoid bones. A pressure distribution pattern that could lead to an area of high tension in the area of the parasagittal groove was found. This result could help explain the high incidence of lateral condylar fractures that initiate in this location. The second part of this study focused on the development and validation of a finite element model of the metacarpophalangeal joint. A model was created based on computed tomography (CT) data. It included the third metacarpal bone, the proximal phalanx, the proximal sesamoid bones, the suspensory ligament, medial and lateral collateral ligaments, medial and lateral collateral sesamoidean ligaments, medial and lateral oblique sesamoidean ligaments, and the straight sesamoidean ligament. The mesh resolution was varied to create three models to allow for convergence. The converged model was then validated using data from the previous part of the study as well as data from the literature. The result was a finite element model containing 121,533 nodes, 112,633 hexahedral elements, and 10 non-linear springs. The final section of this study used the converged and validated finite element model to study the effect of varying bone geometry. The model was morphed based on CT data from three horses: control, lateral condylar fracture, and contralateral limb to lateral condylar fracture. There was an area similar between all three groups of increased stress in the palmar aspect of the parasagittal grooves where fractures are thought to initiate. Other results showed distinct differences in the stress distribution pattern between the three groups. Further investigation into these differences may help increase the understanding of a horse's predisposition to injury. In conclusion, this study has shown that joint geometry plays a role in the stress distribution patterns found in the equine metacarpophalangeal joint. The differences in these patterns between the three groups may help explain the increased risk of a catastrophic injury for some horses. Further studies are warranted to better define the parameters leading to these changes.Item Open Access Effect of mixed-mode loading on fatigue and fracture assessment of a steel twin box-girder bridge(Colorado State University. Libraries, 2019) Irfaee, Mazin M., author; Mahmoud, Hussam, advisor; Heyliger, Paul, committee member; Atadero, Rebecca, committee member; Stright, Lisa, committee memberSteel twin box-girders are considered an attractive option for the construction of bridges due to their basic design, simple form, and ease of creation. Despite their advantages, they are considered fracture critical and as such there is an additional mandate for these bridges to inspected more in depth. This causes their inspection cost to be approximately two to five times greater than that of bridges with non-fracture critical members. The required additional inspection in the U.S. is mainly driven by rare historical events of bridge collapse for bridges that were not steel twin box girders. In addition, the mandated additional inspection does not reflect the inherent level of redundancy in most bridges. Therefore, it is important to quantify the potential for fracture and the level of redundancy in steel two-girder bridges in general, and twin box girders in particular, to minimize their inspection cost. Recognizing the inherently large scatter in fatigue performance, evaluating crack propagation and potential for fracture should, however, be performed in a probabilistic manner using detailed models that represent accurate behavior of the bridge. In this study, a detailed numerical finite element model of steel twin tub-girder bridge is developed and crack growth analysis, potential for fracture of its main tubs, and its overall redundancy is evaluated. The crack growth analysis is performed using multi-mode linear elastic fracture mechanics while accounting for uncertainties in the random variables associated with crack propagation and fracture. The results of the crack growth analysis are utilized to develop fragility functions that specify inspection intervals versus probability of failure where failure is characterized by dynamic crack growth. The analysis conducted to quantify the potential for fracture show distinct possible failure modes that vary from brittle fracture to ductile fracture. The extreme loading case shows that the bridge overall is not at risk of collapse. It is important to note that this conclusion cannot be generalized for all tub girder bridges since the level of redundancy is expected to vary between bridges depending on many factors such as girders geometries, plate thickness, fabrication, among others. However, the presented approach and the corresponding results provide a systematic way by which fracture critical bridges can be evaluated.Item Open Access Evaluating the bond durability of FRP-concrete systems subjected to environmental exposures(Colorado State University. Libraries, 2012) Mata Carrillo, Oscar Rafael, author; Atadero, Rebecca, advisor; Heyliger, Paul, committee member; Glick, Scott, committee memberThe poor current condition of transportation infrastructure in the U.S. is well documented. With traffic volumes on the rise, as well as limited funding available to maintain and rehabilitate aging bridges, cost effective means of improving the performance and durability of these structures must be employed. Fiber Reinforced Polymers (FRPs) offer one potential solution. Their use has been progressively growing in the field of civil engineering as the material's high strength to weight ratio, non-corrosive nature, and ability to conform to existing geometry make it appealing in the reinforcement of existing reinforced concrete structures. In most applications of FRP to strengthen an existing structure, the FRP-concrete bond is essential. Bond is needed for proper transfer of stresses among interfaces. From a durability standpoint, the long-term bond performance is also a major concern. As a result, a long-term durability study was conducted in the laboratory to evaluate the behavior of the bond between the FRP and concrete. Small concrete specimens were prepared, reinforced with FRP material, and subjected to various environmental scenarios such as wet-dry cycles, freeze-thaw cycles, and constant immersion in water, as well as deicing agents. Direct tension pull-off tests and three-point flexural tests were conducted on these specimens to determine any degradation in bond strength over time. Finally, the pull-off test method was evaluated by means of previous research studies and recommendations about preparation procedures were made.Item Open Access Hydrological assessment of field-scale GeoWaste and waste rock test piles(Colorado State University. Libraries, 2020) Hassanzadeh Gorakhki, Mohammad Reza, author; Bareither, Christopher, advisor; Shackelford, Charles, committee member; Scalia, Joseph, committee member; Heyliger, Paul, committee member; Butters, Greg, committee memberMine waste rock and mine tailings are generated in substantial quantities an d must be managed to protect human health and the environment. Challenges in mine waste management facilities include geotechnical stability, environmental contamination, water management, and post operation (long term) closure. Waste rock and tailings co-disposal is a management technique that can address many of the aforementioned challenges. GeoWaste is a mixture of fast-filtered tailings and waste rock blended to isolate waste rock particles within a tailings-dominated matrix. A field-scale experiment that included a waste rock pile and GeoWaste pile was conducted at a mine in Central America to evaluate if GeoWaste suppresses sulfide oxidation and production of metal-rich acid rock drainage relative to waste rock. The objectives of this study were to (i) evaluate hydrologic performance of the piles, (ii) conduct in situ infiltration tests on the piles, (iii) determine field-scale hydraulic parameters for GeoWaste and waste rock, and (iv) develop numerical models to predict water content and oxygen concentrations within the piles. Water content, temperature, electrical conductivity, and oxygen concentration within the piles were monitored for 26 months. Sealed double ring infiltrometer tests were conducted at the end of the pile experiment and test pile subsequently were excavated to assess the spatial distribution in geotechnical characteristics. Inverse modeling was completed in HYDRUS-2D based on infiltration data to determine hydraulic conductivity and moisture retention parameters for the test piles. Field- and laboratory-scale hydraulic parameters were used in HYDRUS-1D and HYDRUS-2D to develop seepage models to predict moisture movement during the 26-month pile experiment. Oxygen concentration was predicted for the GeoWaste pile in HYDRUS-1D via the solute transport module, Fick's 2nd law, the oxygen consumption rate, oxygen diffusion in gas and water phases, and Henry's constant.Item Open Access Integrated optimization of composite structures(Colorado State University. Libraries, 2022) Lang, Daniel, author; Radford, Donald, advisor; Herber, Daniel, committee member; Chong, Edwin, committee member; Heyliger, Paul, committee memberMany industries are exploring the application of composite materials to structural designs to reduce weight. A common issue that is encountered by these industries, however, is difficulty in developing structural geometries best suited for the materials. Research efforts have begun to develop optimization methodology to help develop structural shapes but have thus far only partially addressed optimization of the geometry. This dissertation provides a literature review of past efforts to develop optimization methodologies. Through that review it is identified that the subprocesses required to fully optimize a composite structure are mold shape optimization, ply draping analysis, kinematic partitioning, connection and joint definition, ply topology optimization and manufacturing simulations. To date, however, these subprocesses have primarily been applied individually and have not been integrated to develop fully optimized designs. In this research, a methodology is proposed to integrate established composite design and subprocesses to develop optimized composite structures. The proposed methodology sequentially and iteratively improves the design through mold shape optimization, ply draping analysis, kinematic partitioning, connection and joint definition, ply topology optimization and manufacturing simulations. Throughout the proposed methodology, checks are also integrated to ensure that the developed design meets design objectives and constraints. To test the methodology a case study is conducted to develop composite rail vehicle structures. As part of this case study, it is hypothesized that a composite structure designed through a fully integrated methodology will demonstrate reduced costs, mass and improved manufacturability compared to a structure where functions have only been partially integrated. When the proposed fully integrated methodology is applied to create a case study design, the hypothesis is validated. The design generated by the fully integrated optimization methodology has a 37% lower mass and a 56% lower cost to manufacture than a design that is developed through a partially integrated methodology. The case study also demonstrates that structures developed through the proposed methodology have improved manufacturability.Item Open Access Interfacial modifications in fiber reinforced geopolymer matrix composites for improved toughness(Colorado State University. Libraries, 2017) Jackson, Patrick R., author; Radford, Donald, advisor; Heyliger, Paul, committee member; Kota, Arun, committee member; Ma, Kaka, committee memberGeopolymers have emerged in the recent decades as a potential matrix material for advanced composites. Geopolymers, or more generically inorganic polymers, extend the use temperature range over more common organic polymers, while retaining relatively low processing temperatures. In fact, some of the techniques used to process geopolymers are very similar to those developed for thermosetting polymers. This allows for processing of near net shape components without many of the complexities that would be associated with conventional ceramics or metals manufacturing technologies. To-date, fiber reinforced geopolymers have seen limited use, primarily in areas that emphasize high temperature resistance and good manufacturability over structural performance. The use of all-oxide geopolymer matrix composites (GMC) for high temperature structural applications remains uncertain due to limited toughness. Attempts to improve toughness in these materials through the use of an interphase material, such as those associated with ceramic matrix composites (CMC), have yielded mixed results. In some cases, an increase in toughness was observed, but at the expense of modulus and sometimes strength. The result was a composite that was less tough, or no tougher, than the composite with the untailored interface condition. Additionally, methods that would indicate differences between the tailored and untailored interface have not been employed leaving uncertainty as to what is providing improved toughness. This research examines the ability of weak interface concepts, often employed in ceramic matrix composites and created using fiber coatings, as means of producing greater toughness in GMCs exposed to elevated temperatures. This was accomplished through examination of composite mechanical properties and interfacial conditions. Geopolymer matrix composites reinforced with the 3M Nextel series of ceramic fibers were fabricated and exposed to elevated temperatures. From the fabricated composites, samples were prepared for testing in flexure, tension, short beam shear, and single fiber push-out. Microscopy techniques were employed to analyze fracture surfaces and results of push-out testing of composites. Both a model coating, carbon and thermally stable oxide coating known as monazite were applied to the surface of fibers and compared against a baseline condition to support the changes observed. The results of the research indicate the importance of ensuring adequate cure time of the geopolymer matrix, which enhances it properties. In GMCs using carbon coated fibers to achieve a weak interfacial condition, low mechanical properties of inadequately cured matrix produced composites with limited shear resistance and limited ability to transfer stresses to fibers. A moderate increase in the mechanical properties of the matrix via extended cure time from 1 to 5 hours resulted in a roughly 50% increase in modulus and 150% increase in strength for GMCs containing the interphase material. In all cases the use of fiber coatings resulted in a reduction of interfacial strength. This was revealed by fiber push-out testing, which constitutes the first known use of this technique in GMCs to directly analyze the strength of the bond between the fiber and the matrix. Analysis of the interfaces in specimens further revealed that simple reductions in bond strength are not sufficient for producing better toughness and mechanical properties in GMCs, but that there is a delicate interplay between the interface properties and improved mechanical behavior. Increased toughness was observed in the specimens containing the carbon coated and monazite coated fiber surfaces except in the instance where coating was degraded by the oxidizing environment. GMCs containing the monazite coated fiber demonstrated the greatest improvement in toughness. The improvement in toughness was the result of increased damage tolerance and also a roughly ~32 to 44% increase in strength as compared to GMCs without coated fiber surfaces. Both limited and elongated elevated temperature exposure did not limit greater toughness from being achieved in monazite coated fiber GMCs as compared to those composites without coatings. In general, the use of fiber coatings did improve the toughness of GMCs as result of weaker interfacial conditions and it was demonstrated that careful tailoring of the interfacial strength can result in retention of mechanical properties.Item Open Access Manufacturing and testing of spline geometry using carbon fiber reinforced composite(Colorado State University. Libraries, 2016) Jambor, Eric, author; Bradley, Thomas, advisor; Radford, Donald, committee member; Heyliger, Paul, committee memberA model and manufacturing process for the design of carbon fiber reinforced composite spline shafts is developed and validated to investigate the feasibility of using composite splines for use in power transmission applications. Composite torque tubes for power transmission have been employed in various industries for over three decades and have shown up to a 50% mass decrease compared to steel shafts designed for the same use. One limiting factor for the amount of weight reduction achievable is the mechanism used to transfer power to and from the composite tube. Most composite shafts use adhesive bonding, fasteners, press fits, or some combination to join a steel or aluminum yolk or spline to the end of the tube. This research will investigate the feasibility of eliminating these mechanisms by replacing them by molding in splines to the composite torque tube. This will additionally reduce part count and manufacturing time as well as eliminating the heavy metal inserts. To achieve this, an analytical model is developed to investigate the strength of composite spline teeth of involute geometry as well as a composite torque tube. Due to the complex nature of designing with composites these models are supplemented by material models using a composite software package and finite element models (FEM). The involute splined shaft was then manufactured using an iterative approach to refine the sample quality and tested in torsion to failure. Although the peak failure torque had a large range over the samples it can be concluded that with improvements in the manufacturing process using molded composite splines is a feasible method of torque transfer. This can be concluded from the failure modes of the splined shaft as they indicate that the splines were able to adequately transfer the load to the torque tube.Item Open Access Methodology and applications for integrating earthquake aftershock risk into performance-based seismic design(Colorado State University. Libraries, 2015) Nazari Khanmiri, Negari, author; van de Lindt, John, advisor; Heyliger, Paul, committee member; Mahmoud, Hussam, committee member; Senior, Bolivar, committee memberAftershocks have the potential to cause severe damage to buildings and contribute to threaten life safety following a major earthquake. However, their effect on seismic hazard is not explicitly accounted for in modern building design codes, nor in emerging methodologies such as performance-based seismic design. In this dissertation a methodology was developed to systematically integrate aftershock seismic hazard into performance-based earthquake engineering (PBEE). This is achieved through a combination of analytical studies with structural degradation models derived from existing publicly available Network for Earthquake Engineering Simulation (NEES) data as well as numerical models. The design adjustments due to aftershock seismic hazard were calculated for the Direct Displacement Design (DDD) approach for a building portfolio. A comprehensive sensitivity analysis was performed to investigate the effect of different factors such as the location and number of stories of the building and magnitudes of mainshocks and aftershocks on the design adjustments needed. The results of this research will have multiple applications such as allowing code developers to investigate different options for change in structural design to account for aftershock hazard. Aftershock consideration can be an option for stake holders in selection of their design criteria to minimize life-cycle cost. Since aftershock hazard is a major consideration when safety tagging a building following an earthquake, the results of this project will provide insight into quantitatively investigating risk for damaged buildings.