Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stiffness”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Post-impact fatigue of cross-plied, through-the-thickness reinforced carbon/epoxy composites

An experimental investigation of the post-impact fatigue response of integrally woven carbon/epoxy composites was conducted. Five different through-the-thickness (TTT) reinforcing fibers were used in an experimental textile process that produced an integrally woven (0/90/0/90/0/90/0/90/0)(sub T) ply layup with 21K AS4 carbon tow fiber. The resin was Hercules 3501-6, and the five TTT reinforcing fibers were Kevlar, Toray carbon, AS4 carbon, glass, and IM6 carbon. The purpose of this investigation was to study the post-impact fatigue response of these material systems and to identify the optimum TTT fiber. Samples were impacted with one half inch diameter aluminum balls with an average velocity of 543 ft/sec. Post-impact static compression and constant amplitude tension-compression fatigue tests were conducted. Fatigue tests were conducted with a loading ratio of R=-5, and frequency of 4 Hz. Damage growth was monitored using x-radiographic and sectioning techniques and by examining the stress-strain response (across the impact site) throughout the fatigue tests. The static compressive stress versus far-field strain response was nearly linear for all material groups. All the samples had a transverse shear failure mode. The average compressive modulus (from far-field strain) was about 10 Msi. The average post-impact static compressive strength was about 35.5 Ksi. The IM6 carbon sample had a strength of over 40 Ksi, more than 16 percent stronger than average. There was considerable scatter in the S-N data. However, the IM6 carbon samples clearly had the best fatigue response. The response of the other materials, while worse than IM6 carbon, could not be ranked definitively. The initial damage zones caused by the impact loading and damage growth from fatigue loading were similar for all five TTT reinforcing materials. The initial damage zones were circular and consisted of delaminations, matrix cracks and ply cracks. Post-impact fatigue loading caused delamination growth, ply cracking and fiber bundle failures, typically 45 deg from impact load direction. During the initial 97 percent of fatigue life, delaminations, ply cracks and fiber bundle failures primarily grew at and near the impact site. During the final 3 percent of life, damage grew rapidly transverse to the loading direction as a through-the-thickness transverse shear failure. The stress-strain response was typically linear during the initial 50 percent of life, and stiffness dropped about 20 percent during this period. During the next 47 percent of life, stiffness dropped about 34 percent, and the stress-strain response was no longer linear. The stiffness decreased about 23 percent during the final 3 percent of life. These trends were typical of all the materials tested. Therefore, by monitoring stiffness loss, fatigue failure could be accurately anticipated.

Serdinak, Thomas E.↗
NASA NTRS1995

Aspects of passive magnetic levitation based on high-T(sub c) superconducting YBCO thin films

Passive magnetic levitation systems reported in the past were mostly confined to bulk superconducting materials. Here we present fundamental studies on magnetic levitation employing cylindrical permanent magnets floating above high-T(sub c) superconducting YBCO thin films (thickness about 0.3 mu m). Experiments included free floating rotating magnets as well as well-established flexible beam methods. By means of the latter, we investigated levitation and drag force hysteresis as well as magnetic stiffness properties of the superconductor-magnet arrangement. In the case of vertical motion of the magnet, characteristic high symmetry of repulsive (approaching) and attractive (withdrawing) branches of the pronounced force-displacement hysteresis could be detected. Achievable force levels were low as expected but sufficient for levitation of permanent magnets. With regard to magnetic stiffness, thin films proved to show stiffness-force ratios about one order of magnitude higher than bulk materials. Phenomenological models support the measurements. Regarding the magnetic hysteresis of the superconductor, the Irie-Yamafuji model was used for solving the equation of force balance in cylindrical coordinates allowing for a macroscopic description of the superconductor magnetization. This procedure provided good agreement with experimental levitation force and stiffness data during vertical motion. For the case of (lateral) drag force basic qualitative characteristics could be recovered, too. It is shown that models, based on simple asymmetric magnetization of the superconductor, describe well asymptotic transition of drag forces after the change of the magnet motion direction. Virgin curves (starting from equilibrium, i.e. symmetric magnetization) are approximated by a linear approach already reported in literature only. This paper shows that basic properties of superconducting thin films allow for their application to magnetic levitation or - without need of levitation forces, e.g. microgravity - magnetic damping devices.

Schoenhuber, P.↗
NASA NTRS1995

Parametric studies of stitching effectiveness for preventing substructure disbond

A methodology is desired that will allow a designer to select appropriate amounts of through-thickness reinforcement needed to meet design requirements. The goal is to use a relatively simple analysis to minimize the amount of testing that needs to be performed, and to make test results from simple configurations applicable to more general structures. Using this methodology one should be able to optimize the selection of stitching materials, the weight of the yarn, and the stitching density. The analysis approach is to treat substructure disbond as a crack propagation problem. In this approach, the stitches have little influence until a delamination begins to grow. Once the delamination reaches, or extends beyond a stitch, the stitch serves to reduce the strain-energy-release-rate (G) at the crack tip for a given applied load. The reduced G can then be compared to the unstitched materials toughness to predict the load required to further extend the crack. The current model treats the stitch as a simple spring which responds to displacements in the vertical (through-thickness) direction. In concept, this approach is similar to that proposed by other authors. Test results indicate that the model should be refined to include the shearing stiffness of the stitch. The strain-energy-release-rate calculations are performed using a code which uses interconnected higher-order plates to model built-up composite cross-sections. When plates are stacked vertically, the interfacial tractions between the plates can be computed. The plate differential equations are solved in closed-form. The code, called SUBLAM, was developed as part of this section in one dimension. Because of this limitation, rows of stitches are treated as a two-dimensional sheet. The spring stiffness of a row of stitches can be estimated from the stitch material, weight, and density. As a practical and conservative approach, we can assume that the stitch is bonded until a crack passes the stitch location. After the crack passes, it is fully bonded. A series of tests were performed to exercise this methodology and incorporated an attached flange such that the sudden change in thickness initiated a delamination. The analysis was used to estimate the materials' critical G from that of the unstitched specimens. With this data, a prediction was made for the load required to delaminate the stitched specimens. Using the methodology, design charts have been created for simplified geometries. These charts give stitch force, along with G(sub 1) and G(sub 2) as as function of the stitch spring stiffness. Using the charts, it should be possible to determine the stitch spring stiffness and strength required to reduce the G to a desired level. From these parameters, the actual stitching material, weight, and density can be computed.

Flanagan, Gerry↗
NASA NTRS1996

Passive Magnetic Bearing With Ferrofluid Stabilization

A new class of magnetic bearings is shown to exist analytically and is demonstrated experimentally. The class of magnetic bearings utilize a ferrofluid/solid magnet interaction to stabilize the axial degree of freedom of a permanent magnet radial bearing. Twenty six permanent magnet bearing designs and twenty two ferrofluid stabilizer designs are evaluated. Two types of radial bearing designs are tested to determine their force and stiffness utilizing two methods. The first method is based on the use of frequency measurements to determine stiffness by utilizing an analytical model. The second method consisted of loading the system and measuring displacement in order to measure stiffness. Two ferrofluid stabilizers are tested and force displacement curves are measured. Two experimental test fixtures are designed and constructed in order to conduct the stiffness testing. Polynomial models of the data are generated and used to design the bearing prototype. The prototype was constructed and tested and shown to be stable. Further testing shows the possibility of using this technology for vibration isolation. The project successfully demonstrated the viability of the passive magnetic bearing with ferrofluid stabilization both experimentally and analytically.

Jansen, Ralph↗
NASA NTRS1996

A Feasibility Study on a Parallel Mechanism for Examining the Space Shuttle Orbiter Payload Bay Radiators

The goal of the project was to develop the necessary analysis tools for a feasibility study of a cable suspended robot system for examining the space shuttle orbiter payload bay radiators These tools were developed to address design issues such as workspace size, tension requirements on the cable, the necessary accuracy and resolution requirements and the stiffness and movement requirements of the system. This report describes the mathematical models for studying the inverse kinematics, statics, and stiffness of the robot. Each model is described by a matrix. The manipulator Jacobian was also related to the stiffness matrix, which characterized the stiffness of the system. Analysis tools were then developed based on the singular value decomposition (SVD) of the corresponding matrices. It was demonstrated how the SVD can be used to quantify the robot's performance and to provide insight into different design issues.

Roberts, Rodney G.↗
NASA NTRS1997

Frictionless Contact of Multilayered Composite Half Planes Containing Layers With Complex Eigenvalues

A previously developed local-global stiffness matrix methodology for the response of a composite half plane, arbitrarily layered with isotropic, orthotropic or monoclinic plies, to indentation by a rigid parabolic punch is further extended to accommodate the presence of layers with complex eigenvalues (e.g., honeycomb or piezoelectric layers). First, a generalized plane deformation solution for the displacement field in an orthotropic layer or half plane characterized by complex eigenvalues is obtained using Fourier transforms. A local stiffness matrix in the transform domain is subsequently constructed for this class of layers and half planes, which is then assembled into a global stiffness matrix for the entire multilayered half plane by enforcing continuity conditions along the interfaces. Application of the mixed boundary condition on the top surface of the half plane indented by a rigid punch results in an integral equation for the unknown pressure in the contact region. The integral possesses a divergent kernel which is decomposed into Cauchy-type and regular parts using the asymptotic properties of the local stiffness matrix and a relationship between Fourier and finite Hilbert transform of the contact pressure. The solution of the resulting singular integral equation is obtained using a collocation technique based on the properties of orthogonal polynomials developed by Erdogan and Gupta. Examples are presented that illustrate the important influence of low transverse properties of layers with complex eigenvalues, such as those exhibited by honeycomb, on the load versus contact length response and contact pressure distributions for half planes containing typical composite materials.

Zhang, Wang↗
NASA NTRS1995

Analytic Sensitivity and Approximation of Skin Buckling Constraints in Wing-Shape Synthesis

Explicit expressions for terms of the stiffness and geometric stiffness matrices are derived for the buckling analysis of trapezoidal fiber composite wing skin panels. The formulation is based on Ritz analysis using simple polynomials, and leads to explicit expressions for the analytic sensitivities of the stiffness and geometric stiffness matrices with respect to layer thickness, fiber directions, and panel shape. Integration with wing box analysis using either the equivalent plate approach or the finite element method, makes it possible to obtain sensitivities of panel buckling constraints with respect to wing planform shape or locations of internal ribs and spars. The analytic sensitivities are used to construct approximations of panel buckling constraints for integrated wing/panel design synthesis.

Livne, Eli↗
NASA NTRS2001

Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations

Additive Runge-Kutta (ARK) methods are investigated for application to the spatially discretized one-dimensional convection-diffusion-reaction (CDR) equations. First, accuracy, stability, conservation, and dense output are considered for the general case when N different Runge-Kutta methods are grouped into a single composite method. Then, implicit-explicit, N = 2, additive Runge-Kutta ARK2 methods from third- to fifth-order are presented that allow for integration of stiff terms by an L-stable, stiffly-accurate explicit, singly diagonally implicit Runge-Kutta (ESDIRK) method while the nonstiff terms are integrated with a traditional explicit Runge-Kutta method (ERK). Coupling error terms are of equal order to those of the elemental methods. Derived ARK2 methods have vanishing stability functions for very large values of the stiff scaled eigenvalue, z(exp [I]) goes to infinity, and retain high stability efficiency in the absence of stiffness, z(exp [I]) goes to zero. Extrapolation-type stage-value predictors are provided based on dense-output formulae. Optimized methods minimize both leading order ARK2 error terms and Butcher coefficient magnitudes as well as maximize conservation properties. Numerical tests of the new schemes on a CDR problem show negligible stiffness leakage and near classical order convergence rates. However, tests on three simple singular-perturbation problems reveal generally predictable order reduction. Error control is best managed with a PID-controller. While results for the fifth-order method are disappointing, both the new third- and fourth-order methods are at least as efficient as existing ARK2 methods while offering error control and stage-value predictors.

Kennedy, Christopher A.↗
NASA NTRS2003

Multibody Parachute Flight Simulations for Planetary Entry Trajectories Using "Equilibrium Points"

A method has been developed to reduce numerical stiffness and computer CPU requirements of high fidelity multibody flight simulations involving parachutes for planetary entry trajectories. Typical parachute entry configurations consist of entry bodies suspended from a parachute, connected by flexible lines. To accurately calculate line forces and moments, the simulations need to keep track of the point where the flexible lines meet (confluence point). In previous multibody parachute flight simulations, the confluence point has been modeled as a point mass. Using a point mass for the confluence point tends to make the simulation numerically stiff, because its mass is typically much less that than the main rigid body masses. One solution for stiff differential equations is to use a very small integration time step. However, this results in large computer CPU requirements. In the method described in the paper, the need for using a mass as the confluence point has been eliminated. Instead, the confluence point is modeled using an "equilibrium point". This point is calculated at every integration step as the point at which sum of all line forces is zero (static equilibrium). The use of this "equilibrium point" has the advantage of both reducing the numerical stiffness of the simulations, and eliminating the dynamical equations associated with vibration of a lumped mass on a high-tension string.

Raiszadeh, Ben↗
NASA NTRS2002

Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations

Additive Runge-Kutta (ARK) methods are investigated for application to the spatially discretized one- dimensional convection-diffusion-reaction (CDR) equations. Accuracy, stability, conservation, and dense-output are first considered for the general case when N different Runge-Kutta methods are grouped into a single composite method. Then, implicit-explicit, (N = 2), additive Runge-Kutta (ARK(sub 2)) methods from third- to fifth-order are presented that allow for integration of stiff terms by an L-stable, stiffly-accurate explicit, singly diagonally implicit Runge-Kutta (ESDIRK) method while the nonstiff terms are integrated with a traditional explicit Runge-Kutta method (ERK). Coupling error terms of the partitioned method are of equal order to those of the elemental methods. Derived ARK(sub 2) methods have vanishing stability functions for very large values of the stiff scaled eigenvalue, z['] yields -infinity, and retain high stability efficiency in the absence of stiffness, z['] yield 0. Extrapolation-type stage- value predictors are provided based on dense-output formulae. Optimized methods minimize both leading order ARK(sub 2) error terms and Butcher coefficient magnitudes as well as maximize conservation properties. Numerical tests of the new schemes on a CDR problem show negligible stiffness leakage and near classical order convergence rates. However, tests on three simple singular-perturbation problems reveal generally predictable order reduction. Error control is best managed with a PID-controller. While results for the fifth-order method are disappointing, both the new third- and fourth-order methods are at least as efficient as existing ARK(sub 2) methods.

Kennedy, Christopher A.↗
NASA NTRS2006

Equivalent Plate Analysis of Aircraft Wing with Discrete Source Damage

An equivalent plate procedure is developed to provide a computationally efficient means of matching the stiffness and frequencies of flight vehicle wing structures for prescribed loading conditions. First, the equivalent plate is used to match the stiffness of a stiffened panel without damage and the stiffness of a stiffened panel with damage. For both stiffened panels, the equivalent plate models reproduce the deformation of a corresponding detailed model exactly for the given loading conditions. Once the stiffness was matched, the equivalent plate models were then used to predict the frequencies of the panels. Two analytical procedures using the lumped-mass matrix were used to match the first five frequencies of the corresponding detailed model. In both the procedures, the lumped-mass matrix for the equivalent plate is constructed by multiplying the diagonal terms of the consistent-mass matrix by a proportionality constant. In the first procedure, the proportionality constant is selected such that the total mass of the equivalent plate is equal to that of the detailed model. In the second method, the proportionality constant is selected to minimize the sum of the squares of the errors in a set of pre-selected frequencies between the equivalent plate model and the detailed model. The equivalent plate models reproduced the fundamental first frequency accurately in both the methods. It is observed that changing only the mass distribution in the equivalent plate model did not provide enough flexibility to match all of the frequencies.

Krishnamurthy, T.↗
NASA NTRS1992

Elastomeric load sharing device

An elastomeric load sharing device, interposed in combination between a driven gear and a central drive shaft to facilitate balanced torque distribution in split power transmission systems, includes a cylindrical elastomeric bearing and a plurality of elastomeric bearing pads. The elastomeric bearing and bearing pads comprise one or more layers, each layer including an elastomer having a metal backing strip secured thereto. The elastomeric bearing is configured to have a high radial stiffness and a low torsional stiffness and is operative to radially center the driven gear and to minimize torque transfer through the elastomeric bearing. The bearing pads are configured to have a low radial and torsional stiffness and a high axial stiffness and are operative to compressively transmit torque from the driven gear to the drive shaft. The elastomeric load sharing device has spring rates that compensate for mechanical deviations in the gear train assembly to provide balanced torque distribution between complementary load paths of split power transmission systems.

Isabelle, Charles J.↗
NASA NTRS2010

Frequencies and Flutter Speed Estimation for Damaged Aircraft Wing Using Scaled Equivalent Plate Analysis

Equivalent plate analysis is often used to replace the computationally expensive finite element analysis in initial design stages or in conceptual design of aircraft wing structures. The equivalent plate model can also be used to design a wind tunnel model to match the stiffness characteristics of the wing box of a full-scale aircraft wing model while satisfying strength-based requirements An equivalent plate analysis technique is presented to predict the static and dynamic response of an aircraft wing with or without damage. First, a geometric scale factor and a dynamic pressure scale factor are defined to relate the stiffness, load and deformation of the equivalent plate to the aircraft wing. A procedure using an optimization technique is presented to create scaled equivalent plate models from the full scale aircraft wing using geometric and dynamic pressure scale factors. The scaled models are constructed by matching the stiffness of the scaled equivalent plate with the scaled aircraft wing stiffness. It is demonstrated that the scaled equivalent plate model can be used to predict the deformation of the aircraft wing accurately. Once the full equivalent plate geometry is obtained, any other scaled equivalent plate geometry can be obtained using the geometric scale factor. Next, an average frequency scale factor is defined as the average ratio of the frequencies of the aircraft wing to the frequencies of the full-scaled equivalent plate. The average frequency scale factor combined with the geometric scale factor is used to predict the frequency response of the aircraft wing from the scaled equivalent plate analysis. A procedure is outlined to estimate the frequency response and the flutter speed of an aircraft wing from the equivalent plate analysis using the frequency scale factor and geometric scale factor. The equivalent plate analysis is demonstrated using an aircraft wing without damage and another with damage. Both of the problems show that the scaled equivalent plate analysis can be successfully used to predict the frequencies and flutter speed of a typical aircraft wing.

Krishnamurthy, Thiagarajan↗
NASA NTRS2011

Enhanced Correlation of SMART Active Flap Rotor Loads

This is a follow-on study to a 2010 correlation effort. Measured data from the SMART rotor test in the NASA Ames 40- by 80- Foot Wind Tunnel are compared with CAMRAD II calculations. As background, during the wind tunnel test, unexpectedly high inboard loads were encountered, and it was hypothesized at that time that due to changes in the flexbeams over the years, the flexbeam properties used in the analysis needed updating. Boeing Mesa, recently updated these properties. This correlation study uses the updated flexbeam properties. Compared to earlier studies, the following two enhancements are implemented: i) the inboard loads (pitchcase and flexbeam loads) correlation is included for the first time (reliable prediction of the inboard loads is a prerequisite for any future anticipated flight-testing); ii) the number of blade modes is increased to better capture the flap dynamics and the pitchcase-flexbeam dynamics. Also, aerodynamically, both the rolled-up wake model and the more complex, multiple trailer wake model are used, with the latter slightly improving the blade chordwise moment correlation. This sensitivity to the wake model indicates that CFD is needed. Three high-speed experimental cases, one uncontrolled free flap case and two commanded flap cases, are considered. The two commanded flap cases include a 2o flap deflection at 5P case and a 0o flap deflection case. For the free flap case, selected modifications to the HH-06 section flap airfoil pitching moment table are implemented. For the commanded 2o flap case, the experimental flap variation is approximately matched by increasing the analytical flap hinge stiffness. This increased flap hinge stiffness is retained for the commanded 0o flap case also, which is treated as a free flap case, but with larger flap hinge stiffness. The change in the mid-span and outboard loads correlation due to the updating of the flexbeam properties is not significant. Increasing the number of blade modes results in an effective, commanded flap hinge stiffness of 4X baseline, not 3X as reported earlier. The inboard loads correlation is reasonable, but needs further study. Overall, the free flap case correlation is reasonable, thus confirming the basic correctness of the current semi-empirical modifications; the correlation for the commanded 2o flap at 5P case and the 0o flap case is also reasonable.

Kottapalli, Sesi↗
NASA NTRS2012

Macro Scale Independently Homogenized Subcells for Modeling Braided Composites

An analytical method has been developed to analyze the impact response of triaxially braided carbon fiber composites, including the penetration velocity and impact damage patterns. In the analytical model, the triaxial braid architecture is simulated by using four parallel shell elements, each of which is modeled as a laminated composite. Currently, each shell element is considered to be a smeared homogeneous material. The commercial transient dynamic finite element code LS-DYNA is used to conduct the simulations, and a continuum damage mechanics model internal to LS-DYNA is used as the material constitutive model. To determine the stiffness and strength properties required for the constitutive model, a top-down approach for determining the strength properties is merged with a bottom-up approach for determining the stiffness properties. The top-down portion uses global strengths obtained from macro-scale coupon level testing to characterize the material strengths for each subcell. The bottom-up portion uses micro-scale fiber and matrix stiffness properties to characterize the material stiffness for each subcell. Simulations of quasi-static coupon level tests for several representative composites are conducted along with impact simulations.

Blinzler, Brina J.↗
NASA NTRS2012

Utilization of the Generalized Method of Cells to Analyze the Deformation Response of Laminated Ceramic Matrix Composites

In order to practically utilize ceramic matrix composites in aircraft engine components, robust analysis tools are required that can simulate the material response in a computationally efficient manner. The MAC/GMC software developed at NASA Glenn Research Center, based on the Generalized Method of Cells micromechanics method, has the potential to meet this need. Utilizing MAC/GMC, the effective stiffness properties, proportional limit stress and ultimate strength can be predicted based on the properties and response of the individual constituents. In this paper, the effective stiffness and strength properties for a representative laminated ceramic matrix composite with a large diameter fiber are predicted for a variety of fiber orientation angles and laminate orientations. As part of the analytical study, methods to determine the in-situ stiffness and strength properties of the constituents required to appropriately simulate the effective composite response are developed. The stiffness properties of the representative composite have been adequately predicted for all of the fiber orientations and laminate configurations examined in this study. The proportional limit stresses and strains and ultimate stresses and strains were predicted with varying levels of accuracy, depending on the laminate orientation. However, for the cases where the predictions did not have the desired level of accuracy, the specific issues related to the micromechanics theory were identified which could lead to difficulties that were encountered that could be addressed in future work.

Goldberg, Robert K.↗
NASA NTRS2013

Experimental and Analytical Characterization of the Macromechanical Response for Triaxial Braided Composite Materials

Increasingly, carbon composite structures are being used in aerospace applications. Their highstrength, high-stiffness, and low-weight properties make them good candidates for replacing many aerospace structures currently made of aluminum or steel. Recently, many of the aircraft engine manufacturers have developed new commercial jet engines that will use composite fan cases. Instead of using traditional composite layup techniques, these new fan cases will use a triaxially braided pattern, which improves case performance. The impact characteristics of composite materials for jet engine fan case applications have been an important research topic because Federal regulations require that an engine case be able to contain a blade and blade fragments during an engine blade-out event. Once the impact characteristics of these triaxial braided materials become known, computer models can be developed to simulate a jet engine blade-out event, thus reducing cost and time in the development of these composite jet engine cases. The two main problems that have arisen in this area of research are that the properties for these materials have not been fully determined and computationally efficient computer models, which incorporate much of the microscale deformation and failure mechanisms, are not available. The research reported herein addresses some of the deficiencies present in previous research regarding these triaxial braided composite materials. The current research develops new techniques to accurately quantify the material properties of the triaxial braided composite materials. New test methods are developed for the polymer resin composite constituent and representative composite coupons. These methods expand previous research by using novel specimen designs along with using a noncontact measuring system that is also capable of identifying and quantifying many of the microscale failure mechanisms present in the materials. Finally, using the data gathered, a new hybrid micromacromechanical computer model is created to simulate the behavior of these composite material systems under static and ballistic impact loading using the test data acquired. The model also quantifies the way in which the fiber/matrix interface affects material response under static and impact loading. The results show that the test methods are capable of accurately quantifying the polymer resin under a variety of strain rates and temperature for three loading conditions. The resin strength and stiffness data show a clear rate and temperature dependence. The data also show the hydrostatic stress effects and hysteresis, all of which can be used by researchers developing composite constitutive models for the resins. The results for the composite data reveal noticeable differences in strength, failure strain, and stiffness in the different material systems presented. The investigations into the microscale failure mechanisms provide information about the nature of the different material system behaviors. Finally, the developed computer model predicts composite static strength and stiffness to within 10 percent of the gathered test data and also agrees with composite impact data, where available.

Littell, Justin D.↗
NASA NTRS2016

Effects of Zoledronate and Mechanical Loading during Simulated Weightlessness on Bone Structure and Mechanical Properties

Space flight modulates bone remodeling to favor bone resorption. Current countermeasures include an anti-resorptive drug class, bisphosphonates (BP), and high-force loading regimens. Does the combination of anti-resorptives and high-force exercise during weightlessness have negative effects on the mechanical and structural properties of bone? In this study, we implemented an integrated model to mimic mechanical strain of exercise via cyclical loading (CL) in mice treated with the BP Zoledronate (ZOL) combined with hindlimb unloading (HU). Our working hypothesis is that CL combined with ZOL in the HU model induces additive structural and mechanical changes. Thirty-two C57BL6 mice (male,16 weeks old, n8group) were exposed to 3 weeks of either HU or normal ambulation (NA). Cohorts of mice received one subcutaneous injection of ZOL (45gkg), or saline vehicle, prior to experiment. The right tibia was axially loaded in vivo, 60xday to 9N in compression, repeated 3xweek during HU. During the application of compression, secant stiffness (SEC), a linear estimate of slope of the force displacement curve from rest (0.5N) to max load (9.0N), was calculated for each cycle once per week. Ex vivo CT was conducted on all subjects. For ex vivo mechanical properties, non-CL left femurs underwent 3-point bending. In the proximal tibial metaphysis, HU decreased, CL increased, and ZOL increased the cancellous bone volume to total volume ratio by -26, +21, and +33, respectively. Similar trends held for trabecular thickness and number. Ex vivo left femur mechanical properties revealed HU decreased stiffness (-37),and ZOL mitigated the HU stiffness losses (+78). Data on the ex vivo Ultimate Force followed similar trends. After 3 weeks, HU decreased in vivo SEC (-16). The combination of CL+HU appeared additive in bone structure and mechanical properties. However, when HU + CL + ZOL were combined, ZOL had no additional effect (p0.05) on in vivo SEC. Structural data followed this trend with ZOL not modulating trabecular thickness in CL + NAHU mice. In summary, our integrated model simulates the combination of weightlessness, exercise-induced mechanical strain, and anti-resorptive treatment that astronauts experience during space missions. Based on these results, we conclude that, at the structural and stiffness level, zoledronate treatment during simulated spaceflight does not impede the skeletal response to axial compression. In contrast to our hypothesis, our data show that zoledronate confers no additional mechanical or structural benefit beyond those gained from cyclical loading.

Mechanical Loading↗