Search NASA⌕ Search

SEARCH · Search NASA

Results for “STRESS COMPUTER PROGRAM”

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 577 records · Page 32

HOST structural analysis program overview

Hot-section components of aircraft gas turbine engines are subjected to severe thermal structural loading conditions, especially during the startup and takeoff portions of the engine cycle. The most severe and damaging stresses and strains are those induced by the steep thermal gradients induced during the startup transient. These transient stresses and strains are also the most difficult to predict, in part because the temperature gradients and distributions are not well known or readily predictable and, in part, because the cyclic elastic-viscoplastic behavior of the materials at these extremes of temperature and strain are not well known or readily predictable. A broad spectrum of structures related technology programs is underway to address these deficiencies at the basic as well as the applied level. The three key program elements in the HOST structural analysis program are computations, constitutive modeling, and experiments for each research activity. Also shown are tables summarizing each of the activities.

Thompson, Robert L.↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗

Theoretical Model of the Effect of Crack Tip Blunting on the Ultimate Tensile Strength of Welds in 2219-T87 Aluminum

A theoretical model representing blunting of a crack tip radius through diffusion of vacancies is presented. The model serves as the basis for a computer program which calculates changes, due to successive weld heat passes, in the ultimate tensile strength of 2219-T81 aluminum. In order for the model to yield changes of the same order in the ultimate tensile strength as that observed experimentally, a crack tip radius of the order of .001 microns is required. Such sharp cracks could arise in the fusion zone of a weld from shrinkage cavities or decohered phase boundaries between dendrites and the eutectic phase, or, possibly, from plastic deformation due to thermal stresses encountered during the welding process. Microstructural observations up to X2000 (resolution of about .1 micron) did not, in the fusion zone, show structural details which changed significantly under the influence of a heat pass, with the exception of possible small changes in the configuration of the interdendritic eutectic and in porosity build-up in the remelt zone.

Beil, R. J.↗

Probabilistic/Fracture-Mechanics Model For Service Life

Computer program makes probabilistic estimates of lifetime of engine and components thereof. Developed to fill need for more accurate life-assessment technique that avoids errors in estimated lives and provides for statistical assessment of levels of risk created by engineering decisions in designing system. Implements mathematical model combining techniques of statistics, fatigue, fracture mechanics, nondestructive analysis, life-cycle cost analysis, and management of engine parts. Used to investigate effects of such engine-component life-controlling parameters as return-to-service intervals, stresses, capabilities for nondestructive evaluation, and qualities of materials.

Watkins, T., Jr.↗

Design of ceramic components with the NASA/CARES computer program

The ceramics analysis and reliability evaluation of structures (CARES) computer program is described. The primary function of the code is to calculate the fast-fracture reliability or failure probability of macro-scopically isotropic ceramic components. These components may be subjected to complex thermomechanical loadings, such as those found in heat engine applications. CARES uses results from MSC/NASTRAN or ANSYS finite-element analysis programs to evaluate how inherent surface and/or volume type flaws component reliability. CARES utilizes the Batdorf model and the two-parameter Weibull cumulative distribution function to describe the effects of multiaxial stress states on material strength. The principle of independent action (PIA) and the Weibull normal stress averaging models are also included. Weibull material strength parameters, the Batdorf crack density coefficient, and other related statistical quantities are estimated from four-point bend bar or uniform uniaxial tensile specimen fracture strength data. Parameter estimation can be performed for a single or multiple failure modes by using a least-squares analysis or a maximum likelihood method. Kolmogorov-Smirnov and Anderson-Darling goodness-to-fit-tests, 90 percent confidence intervals on the Weibull parameters, and Kanofsky-Srinivasan 90 percent confidence band values are also provided. Examples are provided to illustrate the various features of CARES.

Nemeth, Noel N.↗

Documentation for a Structural Optimization Procedure Developed Using the Engineering Analysis Language (EAL)

This report describes a structural optimization procedure developed for use with the Engineering Analysis Language (EAL) finite element analysis system. The procedure is written primarily in the EAL command language. Three external processors which are written in FORTRAN generate equivalent stiffnesses and evaluate stress and local buckling constraints for the sections. Several built-up structural sections were coded into the design procedures. These structural sections were selected for use in aircraft design, but are suitable for other applications. Sensitivity calculations use the semi-analytic method, and an extensive effort has been made to increase the execution speed and reduce the storage requirements. There is also an approximate sensitivity update method included which can significantly reduce computational time. The optimization is performed by an implementation of the MINOS V5.4 linear programming routine in a sequential liner programming procedure.

Martin, Carl J., Jr.↗

Automated preliminary design of simplified wing structures to satisfy strength and flutter requirements

A simple structural model of an aircraft wing is used to show the effects of strength (stress) and flutter requirements on the design of minimum-weight aircraft-wing structures. The wing is idealized as an isotropic sandwich plate with a variable cover thickness distribution and a variable depth between covers. Plate theory is used for the structural analysis, and piston theory is used for the unsteady aerodynamics in the flutter analysis. Mathematical programming techniques are used to find the minimum-weight cover thickness distribution which satisfies flutter, strength, and minimum-gage constraints. The method of solution, some sample results, and the computer program used to obtain these results are presented. The results indicate that the cover thickness distribution obtained when designing for the strength requirement alone may be quite different from the cover thickness distribution obtained when designing for either the flutter requirement alone or for both the strength and flutter requirements concurrently. This conclusion emphasizes the need for designing for both flutter and strength from the outset.

Stroud, W. J.↗

ELAS: A general-purpose computer program for the equilibrium problems of linear structures. Volume 2: Documentation of the program

A general purpose digital computer program for the in-core solution of linear equilibrium problems of structural mechanics is documented. The program requires minimum input for the description of the problem. The solution is obtained by means of the displacement method and the finite element technique. Almost any geometry and structure may be handled because of the availability of linear, triangular, quadrilateral, tetrahedral, hexahedral, conical, triangular torus, and quadrilateral torus elements. The assumption of piecewise linear deflection distribution insures monotonic convergence of the deflections from the stiffer side with decreasing mesh size. The stresses are provided by the best-fit strain tensors in the least squares at the mesh points where the deflections are given. The selection of local coordinate systems whenever necessary is automatic. The core memory is used by means of dynamic memory allocation, an optional mesh-point relabelling scheme and imposition of the boundary conditions during the assembly time.

Utku, S.↗

Estimating Aeroheating of a 3D Body Using a 2D Flow Solver

A method for rapidly estimating the aeroheating, shear stress, and other properties of hypersonic flow about a three-dimensional (3D) blunt body has been devised. First, the geometry of the body is specified in Cartesian coordinates. The surface of the body is then described by its derivatives, coordinates, and principal curvatures. Next, previously relatively simple equations are used to find, for each desired combination of angle of attack and meridional angle, a scaling factor and the shape of an equivalent axisymmetric body. These factors and equivalent shapes are entered as inputs into a previously developed computer program that solves the two-dimensional (2D) equations of flow in a non-equilibrium viscous shock layer (VSL) about an axisymmetric body. The coordinates in the output of the VSL code are transformed back to the Cartesian coordinates of the 3D body, so that computed flow quantities can be registered with locations in the 3D flow field of interest. In tests in which the 3D bodies were elliptic paraboloids, the estimates obtained by use of this method were found to agree well with solutions of 3D, finite-rate-chemistry, thin-VSL equations for a catalytic body.

Scott, Carl D.↗

Comparison of predicted and measured elastohydrodynamic film thickness in a 20-millimeter-bore ball bearing

Elastohydrodynamic film thicknesses were measured for a 20-mm bore ball bearing using the capacitance technique. The bearing was thrust loaded to 90, 445, and 778 N (20, 100, and 175 lb). The corresponding maximum contact stress on the inner race was 1.28, 2.09, and 2.45 GPa (185 000, 303,000, and 356, 000 psi). Test speeds ranged from 400 to 15,000 rpm. Measurements were taken with four different lubricants: (1) synthetic paraffinic; (2) synthetic paraffinic with additives; (3) synthetic type II aircraft oil; and (4) synthetic cycloaliphatic hydrocarbon traction fluid. The test bearing was mist lubricated. Test temperatures were 27, 65, and 121 C (80, 150, and 250 F). The measured results for the various test parameters were compared to theoretical predictions from computer programs. Also the data were plotted on dimensionless coordinates and compared to several classical isothermal theories.

Coy, J. J.↗

CODSTRAN - Composite durability structural analysis

CODSTRAN (COmposite Durability STRucture ANalysis) a NASA Lewis Center computer program for the prediction of defect growth and fracture of composite structures when subjected to service loads is presented. Organization, capabilities and present status are discussed. Organizational aspects include executive, input, output, analysis and composite mechanics modules. Capabilities include: durability assessment of large structures and complex structural parts from composites, structural response due to static, cyclic, transient impact and thermal loads, and criteria for static, cyclic, and dynamic fracture. At the present state of development some of CODSTRAN's analysis capabilities include composite mechanics, static failures, and lamination residual stresses. An application in which CODSTRAN is used to predict the defect growth in a flat specimen, with a center through-slit under tension is studied. When completed, CODSTRAN will account for geometry and material nonlinearities, environmental effects as well as static, cyclic and dynamic fracture.

Chamis, C. C.↗

A finite element approach for solution of the 3D Euler equations

Prediction of thermal deformations and stresses has prime importance in the design of the next generation of high speed flight vehicles. Aerothermal load computations for complex three-dimensional shapes necessitate development of procedures to solve the full Navier-Stokes equations. This paper details the development of a three-dimensional inviscid flow approach which can be extended for three-dimensional viscous flows. A finite element formulation, based on a Taylor series expansion in time, is employed to solve the compressible Euler equations. Model generation and results display are done using a commercially available program, PATRAN, and vectorizing strategies are incorporated to ensure computational efficiency. Sample problems are presented to demonstrate the validity of the approach for analyzing high speed compressible flows.

Thornton, E. A.↗

Migration of lymphocytes on fibronectin-coated surfaces: temporal evolution of migratory parameters

Lymphocytes typically interact with implanted biomaterials through adsorbed exogenous proteins. To provide a more complete characterization of these interactions, analysis of lymphocyte migration on adsorbed extracellular matrix proteins must accompany the commonly performed adhesion studies. We report here a comparison of the migratory and adhesion behavior of Jurkat cells (a T lymphoblastoid cell line) on tissue culture treated and untreated polystyrene surfaces coated with various concentrations of fibronectin. The average speed of cell locomotion showed a biphasic response to substrate adhesiveness for cells migrating on untreated polystyrene and a monotonic decrease for cells migrating on tissue culture-treated polystyrene. A modified approach to the persistent random walk model was implemented to determine the time dependence of cell migration parameters. The random motility coefficient showed significant increases with time when cells migrated on tissue culture-treated polystyrene surfaces, while it remained relatively constant for experiments with untreated polystyrene plates. Finally, a cell migration computer model was developed to verify our modified persistent random walk analysis. Simulation results suggest that our experimental data were consistent with temporally increasing random motility coefficients.

NASA Program Fundamental Space Biology↗

Computational Simulation of the High Strain Rate Tensile Response of Polymer Matrix Composites

A research program is underway to develop strain rate dependent deformation and failure models for the analysis of polymer matrix composites subject to high strain rate impact loads. Under these types of loading conditions, the material response can be highly strain rate dependent and nonlinear. State variable constitutive equations based on a viscoplasticity approach have been developed to model the deformation of the polymer matrix. The constitutive equations are then combined with a mechanics of materials based micromechanics model which utilizes fiber substructuring to predict the effective mechanical and thermal response of the composite. To verify the analytical model, tensile stress-strain curves are predicted for a representative composite over strain rates ranging from around 1 x 10(exp -5)/sec to approximately 400/sec. The analytical predictions compare favorably to experimentally obtained values both qualitatively and quantitatively. Effective elastic and thermal constants are predicted for another composite, and compared to finite element results.

Goldberg, Robert K.↗

Development and evaluation of a fault-tolerant multiprocessor (FTMP) computer. Volume 4: FTMP executive summary

The FTMP architecture is a high reliability computer concept modeled after a homogeneous multiprocessor architecture. Elements of the FTMP are operated in tight synchronism with one another and hardware fault-detection and fault-masking is provided which is transparent to the software. Operating system design and user software design is thus greatly simplified. Performance of the FTMP is also comparable to that of a simplex equivalent due to the efficiency of fault handling hardware. The FTMP project constructed an engineering module of the FTMP, programmed the machine and extensively tested the architecture through fault injection and other stress testing. This testing confirmed the soundness of the FTMP concepts.

Smith, T. B., III↗

Description of a computerized method for predicting thermal fatigue life of metals

A computer program is described which can be used to predict the thermal fatigue life of metals and structural components from conventional metal properties. A unique feature of TFLIFE is that it calculates lives according to several different failure criteria for the same input data. These criteria are surface crack initiation, interior crack initiation, and complete fracture of both unnotched and notched fatigue specimens. Sample output tables are shown, together with results for two typical problems: (1) Thermal-mechanical fatigue of bar specimens of the tantalum alloy T-111, and (2) thermal-stress fatigue of wedge specimens of the nickel alloy B-1900.

Spera, D. A.↗

Life prediction and constitutive behavior: Overview

The evolution of programs to investigate high temperature consititutive behavior and develop cyclic life prediction methods is reviewed. Contracts granted for developing and verifying workable engineering methods for the calculation, in advance of service, of the local stress-strain response at the critical life governing location in typical hot section components as well as the resultant cyclic crack initiation and crack growth lifetimes are listed. The Langley fatigue facility is being upgraded to include: (1) a servocontrolled testing machine for high temperature crack growth; (2) three servocontrolled tension/torsion machines for biaxial studies; (3) a HOST/satellite computer for data acquisition, processing, storage, and retrieval; and (4) HCV/LCF machines for cumulative damage studies.

Halford, G. R.↗