Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “materials modelling”

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 415 records · Page 23

Representing Matrix Cracks Through Decomposition of the Deformation Gradient Tensor in Continuum Damage Mechanics Methods

A method is presented to represent the large-deformation kinematics of intraply matrix cracks and delaminations in continuum damage mechanics (CDM) constitutive material models. The method involves the additive decomposition of the deformation gradient tensor into 'crack' and 'bulk material' components. The response of the intact bulk material is represented by a reduced deformation gradient tensor, and the opening of an embedded cohesive interface is represented by a normalized cohesive displacement-jump vector. The rotation of the embedded interface is tracked as the material deforms and as the crack opens. The distribution of the total local deformation between the bulk material and the cohesive interface components is determined by minimizing the difference between the cohesive stress and the bulk material stress projected onto the cohesive interface. The improvements to the accuracy of CDM models that incorporate the presented method over existing approaches are demonstrated for a single element subjected to simple shear deformation and for a finite element model of a unidirectional open-hole tension specimen. The material model is implemented as a VUMAT user subroutine for the Abaqus/Explicit finite element software. The presented deformation gradient decomposition method reduces the artificial load transfer across matrix cracks subjected to large shearing deformations, and avoids the spurious secondary failure modes that often occur in analyses based on conventional progressive damage models.

Leone, Frank A., Jr.↗

A combined experimental and numerical approach that eliminates the non-uniqueness associated with the Johnson-Cook parameters obtained using inverse methods

Abstract Johnson-Cook constitutive model is a commonly used material model for machining simulations. The model includes five parameters that capture the initial yield stress, strain-hardening, strain-rate hardening, and thermal softening behavior of the material. These parameters are difficult to determine using experiments since the conditions observed during machining (such as high strain-rates of the order of $$10^5$$ 10 5 /sec - $$10^6$$ 10 6 /sec) are challenging to recreate in the laboratory. To address this problem, several researchers have recently proposed inverse approaches where a combination of experiments and analytical models are used to predict the Johnson-Cook parameters. The errors between the measured cutting forces, chip thicknesses and temperatures and those predicted by analytical models are minimized and the parameters are determined. In this work, it is shown that only two of the five Johnson-Cook parameters can be determined uniquely using inverse approaches. Two different algorithms, namely, Adaptive Memory Programming for Global Optimization (AMPGO) and Particle Swarm Optimization (PSO), are used for this purpose. The extended Oxley’s model is used as the analytical tool for optimization. For determining a parameter’s value, a large range for each parameter is provided as an input to the algorithms. The algorithms converge to several different sets of values for the five Johnson-Cook parameters when all the five parameters are considered as unknown in the optimization algorithm. All of these sets, however, yield the same chip shape and cutting forces in FEM simulations. Further analyses show that only the strain-rate and thermal softening parameters can be determined uniquely and the three parameters present in the strain-hardening term of the Johnson-Cook model cannot be determined uniquely using the inverse method. A combined experimental and numerical approach is proposed to eliminate this determine all parameters uniquely.

42 ENGINEERING↗

Evaluation of advanced materials through experimental mechanics and modelling

Composite materials have been frequently used in aerospace vehicles. Very often defects are inherited during the manufacture and damages are inherited during the construction and services. It becomes critical to understand the mechanical behavior of such composite structure before it can be further used. One good example of these composite structures is the cylindrical bottle of solid rocket motor case with accidental impact damages. Since the replacement of this cylindrical bottle is expensive, it is valuable to know how the damages affects the material, and how it can be repaired. To reach this goal, the damage must be characterized and the stress/strain field must be carefully analyzed. First the damage area, due to impact, is surveyed and identified with a shearography technique which uses the principle of speckle shearing interferometry to measure displacement gradient. Within the damage area of a composite laminate, such as the bottle of solid rocket motor case, all layers are considered to be degraded. Once a lamina being degraded the stiffness as well as strength will be drastically decreased. It becomes a critical area of failure to the whole bottle. And hence the stress/strain field within and around a damage should be accurately evaluated for failure prediction. To investigate the stress/strain field around damages a Hybrid-Numerical method which combines experimental measurement and finite element analysis is used. It is known the stress or strain at the singular point can not be accurately measured by an experimental technique. Nevertheless, if the location is far away from the singular spot, the displacement can be found accurately. Since it reflects the true displacement field locally regardless of the boundary conditions, it is an excellent input data for a finite element analysis to replace the usually assumed boundary conditions. Therefore, the Hybrid-Numerical method is chosen to avoid the difficulty and to take advantage of both experimental technique and finite element analysis. Experimentally, the digital image correlation technique is employed to measure the displacement field. It is done by comparing two digitized images, before and after loading. Numerically, the finite element program, ABAQUS (version 5.2), is used to analyze the stress and strain field. It takes advantage of the high speed and huge memory size of modern supercomputer, CRAY Y-MP, at NASA Marshall Space Flight Center.

Yang, Yii-Ching↗

Thermal Ablation Modeling for Silicate Materials

A thermal ablation model for silicates is proposed. The model includes the mass losses through the balance between evaporation and condensation, and through the moving molten layer driven by surface shear force and pressure gradient. This model can be applied in ablation simulations of the meteoroid or glassy Thermal Protection Systems for spacecraft. Time-dependent axi-symmetric computations are performed by coupling the fluid dynamics code, Data-Parallel Line Relaxation program, with the material response code, Two-dimensional Implicit Thermal Ablation simulation program, to predict the mass lost rates and shape change. For model validation, the surface recession of fused amorphous quartz rod is computed, and the recession predictions reasonably agree with available data. The present parametric studies for two groups of meteoroid earth entry conditions indicate that the mass loss through moving molten layer is negligibly small for heat-flux conditions at around 1 MW/cm(exp. 2).

Ablators↗

Fragmentation Modeling with a Material Point Method [Slides]

Fragmentation modeling showing: development of critical modeling capability for safety; demonstration of crack opening of high explosives; modeling of fragmentation of brittle material; and emphasis on the important of modeling fragments for safety.

42 ENGINEERING↗

Effects of growth substrate on the nucleation of monolayer MoTe 2

The discovery and characterization of two-dimensional (2D) materials beyond graphene has increased dramatically over the past decade with increasingly fine control over the growth dynamics of these materials. MoTe 2 represents a model material for studying phase change properties of 2D materials due to small energy differences between its 2H semiconducting phase and 1T' semimetallic phase. Although some substrates are demonstrated to be better than others for the wafer-scale growth of high quality MoTe 2 films, substrate effects on the nucleation and growth of MoTe 2 are still not well understood. Here, we grow monolayer MoTe 2 by converting MoO x thin films deposited on three different substrates: sapphire Al 2 O 3 (0001), amorphous SiO 2 , and amorphous AlO x , and examine the early stages of the conversion reaction to elucidate the substrate effects on the nucleation of MoTe 2 . We observe that the chemical composition of the substrate is more important than the surface topography and crystallinity of the substrate, with high quality monolayer 2H MoTe 2 formed on both Al 2 O 3 (0001) and AlO x in contrast to mixed phase 2H/1T' MoTe 2 formed on SiO 2 , as determined by Raman spectroscopy, X-ray photoelectron spectroscopy, and atomic force microscopy.

36 MATERIALS SCIENCE↗

On a study of the /Delta T/c and C/asterisk/ integrals for fracture analysis under non-steady creep

Applications of a vector quantity, path-independent integral which has an energy interpretation to the characterization of crack-tip fields in the range from fast to slow crack propagation are examined. The crack tip characterization parameter is defined in terms of a conservation integral for an area around the crack tip in a two-dimensional cracked body. The actual physical interpretation of the parameter is shown to be the difference in crack lengths displayed by two identical bodies which have equal load histories. A steady-state value is obtained for the parameter for cases of steady-state creep and is shown to be related to the standard path-independent integral for macroscopic self-similar crack growth under mode I conditions. A finite element model is developed for viscoplastic material models, using an initial strain approach with steps in a size employed in tangent stiffness methods.

Stonesifer, R. B.↗

A Tabletop Tool for Modeling Life Support Systems

This paper describes the development plan for a comprehensive research and diagnostic tool for aspects of advanced life support systems in space-based laboratories. Specifically it aims to build a high fidelity tabletop model that can be used for the purpose of risk mitigation, failure mode analysis, contamination tracking, and testing reliability. We envision a comprehensive approach involving experimental work coupled with numerical simulation to develop this diagnostic tool. It envisions a 10% scale transparent model of a space platform such as the International Space Station that operates with water or a specific matched index of refraction liquid as the working fluid. This allows the scaling of a 10 ft x 10 ft x 10 ft room with air flow to 1 ft x 1 ft x 1 ft tabletop model with water/liquid flow. Dynamic similitude for this length scale dictates model velocities to be 67% of full-scale and thereby the time scale of the model to represent 15% of the full- scale system; meaning identical processes in the model are completed in 15% of the full- scale time. The use of an index matching fluid (fluid that matches the refractive index of cast acrylic, the model material) allows making the entire model (with complex internal geometry) transparent and hence conducive to non-intrusive optical diagnostics. So using such a system one can test environment control parameters such as core flows (axial flows), cross flows (from registers and diffusers), potential problem areas such as flow short circuits, inadequate oxygen content, build up of other gases beyond desirable levels, test mixing processes within the system at local nodes or compartments and assess the overall system performance. The system allows quantitative measurements of contaminants introduced in the system and allows testing and optimizing the tracking process and removal of contaminants. The envisaged system will be modular and hence flexible for quick configuration change and subsequent testing. The data and inferences from the tests will allow for improvements in the development and design of next generation life support systems and configurations.

Ramachandran, N.↗

Space Laboratory on a Table Top: A Next Generative ECLSS design and diagnostic tool

This paper describes the development plan for a comprehensive research and diagnostic tool for aspects of advanced life support systems in space-based laboratories. Specifically it aims to build a high fidelity tabletop model that can be used for the purpose of risk mitigation, failure mode analysis, contamination tracking, and testing reliability. We envision a comprehensive approach involving experimental work coupled with numerical simulation to develop this diagnostic tool. It envisions a 10% scale transparent model of a space platform such as the International Space Station that operates with water or a specific matched index of refraction liquid as the working fluid. This allows the scaling of a 10 ft x 10 ft x 10 ft room with air flow to 1 ft x 1 ft x 1 ft tabletop model with water/liquid flow. Dynamic similitude for this length scale dictates model velocities to be 67% of full-scale and thereby the time scale of the model to represent 15% of the full- scale system; meaning identical processes in the model are completed in 15% of the full- scale-time. The use of an index matching fluid (fluid that matches the refractive index of cast acrylic, the model material) allows making the entire model (with complex internal geometry) transparent and hence conducive to non-intrusive optical diagnostics. So using such a system one can test environment control parameters such as core flows (axial flows), cross flows (from registers and diffusers), potential problem areas such as flow short circuits, inadequate oxygen content, build up of other gases beyond desirable levels, test mixing processes within the system at local nodes or compartments and assess the overall system performance. The system allows quantitative measurements of contaminants introduced in the system and allows testing and optimizing the tracking process and removal of contaminants. The envisaged system will be modular and hence flexible for quick configuration change and subsequent testing. The data and inferences from the tests will allow for improvements in the development and design of next generation life support systems and configurations. Preliminary experimental and modeling work in this area will be presented. This involves testing of a single inlet-exit model with detailed 3-D flow visualization and quantitative diagnostics and computational modeling of the system.

Ramachandran, N.↗

Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions

Many materials of interest to Sandia transition from fluid to solid or have regions of both phases coexisting simultaneously. Currently there are, unfortunately, no material models that can accurately predict this material response. This is relevant to applications that "birth stress" related to geoscience, nuclear safety, manufacturing, energy production and bioscience. Accurately capturing solidification and residual stress enables fully predictive simulations of the evolving front shape or final product. Accurately resolving flow of proppants or blood could reduce environmental impact or lead to better treatments for heart attacks, thrombosis, or aneurism. We will address a science question in this proposal: When does residual stress develop during the critical transition from liquid to solid and how does it affect material deformation? Our hypothesis is that these early phases of stress development are critical to predictive simulation of material performance, net shape, and aging. In this project, we use advanced constitutive models with yield stress to represent both fluid and solid behavior simultaneously. The report provides an abbreviated description of the results from our LDRD "Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions," since we have written four papers that document the work in detail and which we reference. We give highlights of the work and describe the gravitationally driven flow visualization experiment on a model yield stress fluid, Carbopol, at various concentrations and flow rates. We were able to collapse the data on a single master curve by showing it was self-similar. We also describe the Carbopol rheology and the constitutive equations of interest including the Bingham-Carreau-Yasuda model, the Saramito model, and the HB-Saramito model including parameter estimation for the shear and oscillatory rheology. We present several computational models including the 3D moving mesh simulations of both the Saramito models and Bingham-Carreau-Yasuda (BCY) model. We also show results from the BCY model using a 3D level set method and two different ways of handling reduced order Hele-Shaw modeling for generalized Newtonian fluids. We present some first ever two-dimensional results for the modified Jeffries Kamani-Donley-Rogers constitutive equation developed during this project. We include some recent results with a successful Saramito-level set coupling that allows us to tackle problems with complex geometries like mold filling in a thin gap with an obstacle, without the need for remeshing or remapping. We report on some experiments for curing systems where fluorescent particles are used to track material flow. These experiments were carried out in an oven on Sylgard 184 as a model polymerizing system. We conclude the report with a summary of accomplishments and some thoughts on follow-on work.

36 MATERIALS SCIENCE↗

A Hysteresis Model for Piezoceramic Materials

This paper addresses the modeling of nonlinear constitutive relations and hysteresis inherent to piezoceramic materials at moderate to high drive levels. Such models are, necessary to realize the, full potential of the materials in high performance control applications, and a necessary prerequisite is the development of techniques which permit control implementation. The approach employed here is based on the qualification of reversible and irreversible domain wall motion in response to applied electric fields. A comparison with experimental data illustrates that because the resulting ODE model is physics-based, it can be employed for both characterization and prediction of polarization levels throughout the range of actuator operation. Finally, the ODE formulation is amenable to inversion which facilitates the development of an inverse compensator for linear control design.

Smith, Ralph C.↗

Impact Crater Scaling Laws

The effects of the impactor size and velocity, of the material properties of the impactor and impacted body, and of the gravitational field strength were examined. The dependence of the coupling parameter on the Gruniesen parameter and the us/up slope values that characterize a Tillotson material model were tested. Whether a code calculation can recover the scaling exponents that can be derived theoretically for an idealized perfectly-porous material that crushes at zero strength and subsequently is incompressible was investigated. A model that was close to the idealized case was generated in a form suitable for code calculations and utilized for a series of one dimensional calculations. The code calculation demonstrated the theoretical coupling parameter. Calculations of impacts using a physically real porous material model seem to indicate a coupling parameter of the same form applicable to a nonporous material, distinctly different from that expected from other theoretical and experimental results for porous materials. The one dimensional calculations are being extended to two dimensional impacts.

Holsapple, K. A.↗

Guided Design of Alloys

Using multiscale materials modeling, we are designing stronger alloys for the airspace applications. We focus on improving mechanical properties of materials at operating temperature. We employ basic science and fundamental theory for guided design of the next-generation alloys. Among the examples, we consider precipitated superalloys with a local phase transformation strengthening. We acknowledge funding by NASA’s Aeronautics Research Mission Directorate (ARMD) via Transformational Tools and Technologies (TTT) Project.

multiscale↗

Guided Design of Alloys: Part 2

Using multiscale materials modeling, we are designing stronger alloys for the airspace applications. We focus on improving mechanical properties of materials at operating temperature. We employ basic science and fundamental theory for guided design of the next-generation alloys. Among the examples, we consider precipitated superalloys with a local phase transformation strengthening. We acknowledge funding by NASA’s Aeronautics Research Mission Directorate (ARMD) via Transformational Tools and Technologies (TTT) Project.

Multiscale↗

A Multiscale Nonlocal Progressive Damage Model for Composite Materials

In this paper, the advantages of a nonlocal progressive damage formulation are described and demonstrated. An approximation of the nonlocal formulation was implemented coupled with the MAT162 composite damage model as a User defined material model in the LS DYNA environment. A comparison of the local model and the nonlocal model is simulated for an 8-ply laminate under tension is carried for increasing mesh densities. The results show the regularization achieved by nonlocal models by providing mesh independent results.

Kodagali, Karan↗

Development of a Mesoscale Finite Element Constitutive Model for Fiber Kinking

A mesoscale finite element material model is proposed to analyze structures that fail by the fiber kinking damage mode. To evaluate the assumptions of the mesoscale model, the results were compared with those of a high-fidelity micromechanical model. A direct comparison between the two models shows remarkable correlation, indicating that the key features of the fiber kinking phenomenon are appropriately accounted for in the mesoscale model. The mesoscale model is applied to structural analysis cases to demonstrate the capabilities of the model. A verification study is conducted with an unnotched compression specimen and preliminary validation is demonstrated with a notched compression specimen. The results show that the model is successful at representing the kinematics of fiber kinking while at the same time highlighting the need for further verification and validation.

Bergan, Andrew C.↗

Application of the MAT 213 Composite Impact Model to NASA Problems of Interest

As composite materials are gaining increased use in aircraft components where impact resistance under high-energy impact conditions is important (such as the turbine engine fan case), the need for accurate material models to simulate the deformation, damage, and failure response of polymer matrix composites under impact conditions is becoming more critical.

Polymer Matrix Composites↗

Progress on Quantitative Infrared Thermography at the NASA Langley Aerothermodynamic Laboratory

A quantitative infrared thermography technique is being developed for accurate and high-resolution measurements of surface temperature and heat flux in the hypersonic wind tunnels at NASA Langley Research Center. A hemisphere test campaign was carried out in the Langley Aerothermodynamic Laboratory 20-Inch Mach 6 Air Tunnel to assess the current performance of the technique and to identify remaining challenges. Measurements were obtained for three different model materials over a range of freestream Reynolds numbers. The raw infrared images were converted to temperature via a radiometric calibration of the camera sensor. Corrections were made for losses due to test article emissivity and window transmissivity. A preliminary uncertainty analysis yielded an estimated accuracy within 4 K at moderate viewing angles, similar to that of a standard thermocouple. Surface heat flux was obtained from the measured surface temperatures via a finite-difference solution of the one-dimensional heat equation that modeled thermal variation of the material properties. Comparisons of the measured heating to both computational and theoretical predictions demonstrated agreement within the expected uncertainty of the technique. The heat transfer measurement near the stagnation point was found to be highly sensitive to the test article emissivity, which establishes the need for improved characterization of the wind tunnel model materials for high temperature applications.

infrared thermography↗