Search NASA⌕ Search

SEARCH · Search NASA

Results for “strain simulation”

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 289 records · Page 16

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials (3MDCP)

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Woven↗

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Discrete Element Method↗

Ground Reaction Forces During Locomotion in Simulated Microgravity

Significant losses in bone density and mineral, primarily in the lower extremities have been reported following exposure to weightlessness. Recent investigations suggest that mechanical influences such as bone deformation and strain rate may be critically important in stimulating new bone formation. It was hypothesized that velocity, cadence and harness design would significantly affect lower limb impact forces during treadmill exercise in simulated zero gravity (0G). A ground-based hypogravity simulator was used to investigate which factors affect limb loading during tethered treadmill exercise. A fractional factorial design was used and 12 subjects were studied. The results showed that running on active and passive treadmills in the simulator with a tethering force close to the maximum comfortable level produced similar magnitudes for the peak ground reaction force. It was also found that these maximum forces were significantly lower than those obtained during overground trials, even when the speeds of locomotion in the simulator were 66 % greater than those in 1 G. Cadence had no effect on any of the response variables. The maximum rate of force application (DFDT-Max) was similar for overground running and exercise in simulated 0G, provided that the "weightless subjects ran on a motorized treadmill. These findings have implications for the use of treadmill exercise as a countermeasure for hypokinetic osteoporosis. As the relationship between mechanical factors and osteogenesis becomes better understood, results from human experiments in 0G simulators will help to design in-flight exercise programs that are more closely targeted to generate appropriate mechanical stimuli.

Davis, B. L.↗

Application of experimental stress analysis techniques to speed ball bearings

The object of the project was to determine experimentally the stress distribution in cylindrically hollow and spherically hollow balls. Ball models several times actual bearing ball size were fabricated. Strain gage rosettes will be mounted on the outside and inside surfaces of the models at specific locations and orientations. The models were subjected to a range of static loads simulating actual ball-race contact loads experienced in high-speed bearings operating at DN values from 2 to 4 million. The strain gages were read, and the principal strains determined for each load setting. Principal stresses were calculated and dimensionless stress coefficients determined to permit the evaluation of stresses in hollow balls of any size and load of geometrically similar configuration.

Nypan, L. J.↗

Testing and Life Prediction for Composite Rotor Hub Flexbeams

A summary of several studies of delamination in tapered composite laminates with internal ply-drops is presented. Initial studies used 2D FE models to calculate interlaminar stresses at the ply-ending locations in linear tapered laminates under tension loading. Strain energy release rates for delamination in these laminates indicated that delamination would likely start at the juncture of the tapered and thin regions and grow unstably in both directions. Tests of glass/epoxy and graphite/epoxy linear tapered laminates under axial tension delaminated as predicted. Nonlinear tapered specimens were cut from a full-size helicopter rotor hub and were tested under combined constant axial tension and cyclic transverse bending loading to simulate the loading experienced by a rotorhub flexbeam in flight. For all the tested specimens, delamination began at the tip of the outermost dropped ply group and grew first toward the tapered region. A 2D FE model was created that duplicated the test flexbeam layup, geometry, and loading. Surface strains calculated by the model agreed very closely with the measured surface strains in the specimens. The delamination patterns observed in the tests were simulated in the model by releasing pairs of MPCs along those interfaces. Strain energy release rates associated with the delamination growth were calculated for several configurations and using two different FE analysis codes. Calculations from the codes agreed very closely. The strain energy release rate results were used with material characterization data to predict fatigue delamination onset lives for nonlinear tapered flexbeams with two different ply-dropping schemes. The predicted curves agreed well with the test data for each case studied.

Murri, Gretchen B.↗

Prediction of fatique crack growth under flight-simulation loading with the modified CORPUS model

The CORPUS (Computation Of Retarded Propagation Under Spectrum loading) crack growth prediction model for variable-amplitude loading, as introduced by De Koning, was based on crack closure. It includes a multiple-overload effect and a transition from plane strain to plane stress. In the modified CORPUS model an underload affected zone (ULZ) is introduced, which is significant for flight-simulation loading in view of the once per flight compressive ground load. The ULZ is associated with reversed plastic deformation induced by the underloads after crack closure has already occurred. Predictions of the crack growth fatigue life are presented for a large variety of flight-simulation test series on 2024-T3 sheet specimens in order to reveal the effects of a number of variables: the design stress level, the gust spectrum severity, the truncation level (clipping), omission of small cycles, and the ground stress level. Tests with different load sequences are also included. The trends of the effects induced by the variables are correctly predicted. The quantitative agreement between the predictions and the test results is also satisfactory.

Padmadinata, U. H.↗

A plane strain analysis of the blunted crack tip using small strain deformation plasticity theory

A deformation plasticity analysis of the tip region of a blunted crack in plane strain is presented. The power hardening material is incompressible both elastically and plastically, in order to simulate behavior of a stress freezing material above critical temperature. Stress and displacement fields surrounding the crack tip are presented. The results indicate that the maximum stress seen at the crack tip is indeed limited and is determined by the tensile properties; however, the scale over which the stresses act is dependent on the loading. Comparisons are good between the forward crack tip displacement and micro-fractographic measurments of stretch zones observed in plane strain fracture toughness tests.

Mcgowan, J. J.↗

Simulation of Hypervelocity Impact on Aluminum-Nextel-Kevlar Orbital Debris Shields

An improved hybrid particle-finite element method has been developed for hypervelocity impact simulation. The method combines the general contact-impact capabilities of particle codes with the true Lagrangian kinematics of large strain finite element formulations. Unlike some alternative schemes which couple Lagrangian finite element models with smooth particle hydrodynamics, the present formulation makes no use of slidelines or penalty forces. The method has been implemented in a parallel, three dimensional computer code. Simulations of three dimensional orbital debris impact problems using this parallel hybrid particle-finite element code, show good agreement with experiment and good speedup in parallel computation. The simulations included single and multi-plate shields as well as aluminum and composite shielding materials. at an impact velocity of eleven kilometers per second.

Fahrenthold, Eric P.↗

Measured Interfacial Residual Strains Produced by In-Flight Ice

The formation of ice on aircraft is a highly dynamic process during which ice will expand and contract upon freezing and undergoing changes in temperature. Finite Element Analysis (FEA) simulations were performed investigating the stress build up in a simplified case with uniform temperature changes between an idealized ice sample and acrylic substrate. These results were used to place strain gages on custom-built acrylic and aluminum specimens; these specimens were then placed in icing conditions such that ice was grown on top of the specimen. Tee rosettes were placed in two configurations adjacent to thermocouple sensors. It was hypothesized that the ice would expand on freezing and contract as the temperature of the interface returned to the equilibrium conditions. While results from the aluminum specimens matched this hypothesis, results from the acrylic specimens show a short period of contraction followed by a much larger expansion at the interface, indicating more complex ice growth thermodynamics than anticipated. Samples were observed to delaminate, and the data suggests that the residual strain is significant to the shedding of ice for in-flight applications.

residual strain↗

A dynamic study of fragmentation and energy loss during high velocity impact

Research conducted under this contract can be divided into two main areas: hypervelocity (in the range up to 7 km/s) and high velocity (less than 1 km/s). Work in the former was performed at NASA-MSFC using the Light Gas Gun Facility. The lower velocity studies were conducted at Auburn University using the ballistic gun. The emphasis of the project was on the hypervelocity phenomenon especially in the characterization of the debris cloud formed by the primary impact events. Special devices were made to determine the angular distributions of momentum and energy of the debris cloud as a function of impact conditions. After several iteration processes, it was decided to concentrate on the momentum effort. Prototype devices were designed, fabricated, and tested. These devices were based on the conservation of momentum. Distributions of the debris cloud formed were measured by determining the amount of momentum transferred from the debris cloud to strategically placed pendulum measurement devices. The motion of the pendula was monitored using itegrated opto-interrupters. The distribution of momentum in the debris cloud was found to be a strong function of the impact condition. Small projectiles at high velocities were observed to produce finely dispersed debris whereas large projectiles generated discrete particles in the debris. Results also show that the momentum in the forward direction was enhanced due to the impact. This phenomenon of momentum multiplication was also observed in other studies and in computer simulations. It was initially planned to determine the energy distribution using deformation energy in a rod with strain gauges. Results from preliminary studies show that this technique is acceptable but too tedious. A new technique was explored based on measuring the heating effect of the debris cloud using an IR camera. The feasibility and sensitivity was established at Auburn University. This type of energy distribution measurement method can easily be adapted to the gas gun facility at MSFC. The objective of the lower velocity studies at Auburn was to simulate the damage produced in advanced materials by the lower energy debris cloud.

Zee, Ralph H.↗

Stress corrosion cracking of titanium alloys

Experiments were conducted to study (1) the basic electrochemical behavior of titanium in acid chloride solutions and (2) the response of the metal to dynamic straining in the same evironment. The aim of this group of experiments was to simulate, as nearly as possible, the actual conditions which exist at the tip of a crack. One of the foremost theories proposed to explain the propagation of stress corrosion cracks is a hydrogen embrittlement theory involving the precipitation of embrittling titanium hydrides inside the metal near the crack tip. An initial survey of the basic electrochemical literature indicated that surface hydrides play a critical role in the electrochemistry of titanium in acid solutions. A comprehensive analysis of the effect of surface films, particularly hydrides, on the electrochemical behavior of titanium in these solution is presented.

May, R. C.↗

Computer simulation of earthquakes

In a computer simulation study of earthquakes a seismically active strike slip fault is represented by coupled mechanical blocks which are driven by a moving plate and which slide on a friction surface. Elastic forces and time independent friction are used to generate main shock events, while viscoelastic forces and time dependent friction add aftershock features. The study reveals that the size, length, and time and place of event occurrence are strongly influenced by the magnitude and degree of homogeneity in the elastic, viscous, and friction parameters of the fault region. For example, periodically reoccurring similar events are observed in simulations with near-homogeneous parameters along the fault, whereas seismic gaps are a common feature of simulations employing large variations in the fault parameters. The study also reveals correlations between strain energy release and fault length and average displacement and between main shock and aftershock displacements.

Cohen, S. C.↗

Experimental evaluation criteria for constitutive models of time dependent cyclic plasticity

Notched members were tested at temperatures far above those recorded till now. Simulation of the notch root stress response was accomplished to establish notch stress-strain behavior. Cyclic stress-strain profiles across the net-section were recorded and on-line direct notch strain control was accomplished. Data are compared to three analysis techniques with good results. The objective of the study is to generate experimental data that can be used to evaluate the accuracy of constitutive models of time dependent cyclic plasticity.

Martin, J. F.↗

Sublaminate analysis of interlaminar fracture in composites

A simple analysis method based upon a transverse shear deformation theory and a sublaminate approach is utilized to analyze a mixed-mode edge delamination specimen. The analysis provides closed form expressions for the interlaminar shear stresses ahead of the crack, the total energy release rate, and the energy release rate components. The parameters controlling the behavior are identified. The effect of specimen stacking sequence and delamination interface on the strain energy release rate components is investigated. Results are compared with a finite element simulation for reference. The simple nature of the method makes it suitable for preliminary design analyses which require a large number of configurations to be evaluated quickly and economically.

Armanios, E. A.↗

Interlaminar analysis of laminated composites using a sublaminate approach

A simple analysis method based upon a transverse shear deformation theory and a sublaminate approach is utilized to analyze a Mixed-Mode edge delamination specimen. The analysis provides a closed form distribution of the interlaminar shear stresses ahead of the crack, and the parameters controlling the behavior are identified. The effect of specimen stacking sequence and delamination interface on the strain energy release rate components is investigated. Results are compared with a finite element simulation for reference. The simple nature of the method makes it suitable for preliminary design analyses which require a large number of configurations to be evaluated quickly and economically.

Armanios, E. A.↗

Computational simulation methods for composite fracture mechanics

Structural integrity, durability, and damage tolerance of advanced composites are assessed by studying damage initiation at various scales (micro, macro, and global) and accumulation and growth leading to global failure, quantitatively and qualitatively. In addition, various fracture toughness parameters associated with a typical damage and its growth must be determined. Computational structural analysis codes to aid the composite design engineer in performing these tasks were developed. CODSTRAN (COmposite Durability STRuctural ANalysis) is used to qualitatively and quantitatively assess the progressive damage occurring in composite structures due to mechanical and environmental loads. Next, methods are covered that are currently being developed and used at Lewis to predict interlaminar fracture toughness and related parameters of fiber composites given a prescribed damage. The general purpose finite element code MSC/NASTRAN was used to simulate the interlaminar fracture and the associated individual as well as mixed-mode strain energy release rates in fiber composites.

Murthy, Pappu L. N.↗

Estimation of distributed acoustic lasers

The purpose of this paper is to present the results of a preliminary investigation of an indirect method which can be used to generate a set of equivalent point loads approximating an applied distributed acoustic load. The method combines the source information with the structure's modal characteristics and forces which give the same response as the applied acoustic load. The technique is demonstrated by predicting a set of point forces acting on a thin aluminum panel which simulate a normally-incident random pressure. The technique is verified by comparing measured strain autospectrums with estimated strain autospectrums generated using the predicting forces. Strain errors within 5 dB were achieved. A limited amount of information about the acoustic load could be deduced using the predicted force autospectrums and phase differences between these forces.

Elliott, K. B.↗

Preferential concentration of particles by turbulence

Direct numerical simulation of isotropic turbulence was used to investigate the effect of turbulence on the concentration fields of heavy particles. The hydrodynamic field was computed using 64-cubed points, and a statistically stationary flow was obtained by forcing the low-wave-number components of the velocity field. The particles used in the simulations were time advanced according to Stokes drag law and were also assumed to be much more dense than the fluid. Properties of the particle cloud were obtained by following the trajectories of 1,000,000 particles through the simulated flow fields. Three values of the ratio of the particle time constant to large-scale turbulence time scale were used in the simulations: 0.075, 0.15, and 0.52. The simulations show that the particles collect preferentially in regions of low vorticity and high strain rate. This preferential collection was most pronounced for particles with the instantaneous number density as much as 25 times the mean value. The fact that dense particles collect in regions of low vorticity and high strain in turn implies that turbulence may actually inhibit rather than enhance mixing of particles.

Squires, Kyle D.↗