Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite element model”

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 505 records · Page 28

Nucleation and growth of rolling contact failure of 440C bearing steel

A 'two-body' elasto-plastic finite element model of 2-dimensional rolling and rolling-plus-sliding was developed to treat the effect of surface irregularities. The model consists of a smooth cylinder in contact with a semi-infinite half-space that is either smooth or fitted with one of 0.4 microns deep or 7 microns deep groove, or a 0.4 microns high ridge-like asperity. The model incorporates elastic-linear-kinematic hardening-plastic (ELKP) and non-linear-kinematic hardening-plastic (NLKP) material constitutive relations appropriate for hardened bearing steel and the 440C grade. The calculated contact pressure distribution is Hertzian for smooth body contact, and it displays intense, stationary, pressure spikes superposed on the Hertzian pressure for contact with the grooved and ridged surface. The results obtained for the 0.4 microns deep groove compare well with those reported by Elsharkawy and Hamrock for an EHD lubricated contact. The effect of translating the counterface on the half space as opposed to indenting the half space with the counter face with no translation is studied. The stress and strain values near the surface are found to be similar for the two cases, whereas they are significantly different in the subsurface. It is seen that when tiny shoulders are introduced at the edge of the groove in the finite element model, the incremental plasticity and residual stresses are significantly higher in the vicinity of the right shoulder (rolling direction is from left to right) than at the left shoulder. This may explain the experimental observation that the spall nucleation occurs at the exit end of the artificially planted indents. Pure rolling calculations are compared with rolling + sliding calculations. For a coefficient of friction, mu = 0.1, the effect of friction is found to be small. Efforts were made to identify the material constitutive relations which best describe the deformation characteristics of the bearing steels in the initial few cycles. Elastic-linear-kinematic hardening-plastic (ELKP) material constitutive relations produce less net plastic deformation in the initial stages for a given stress, than seen in experiments. A new set of constitutive relations: non-linear-kinematic hardening-plastic (NLKP) was used. This material model produces more plasticity than the ELKP model and shows promise for treating the net distortions in the early stages. Techniques for performing experimental measurements that can be compared with the finite element calculations were devised. The measurements are being performed on 9mm-diameter, 440C steel cylindrical rolling elements in contact with 12.5 mm-diameter, 52100 steel balls in a 3-ball-rod fatigue test machine operating at 3600 RPM. Artificial, 7 microns deep, indents were inserted on the running track of the cylindrical rolling elements and profilometer measurements of these indents made, before and after the rolling. These preliminary measurements show that the indents are substantially deformed plastically in the process of rolling. The deformations of the groove calculated with the finite element model are comparable to those measured experimentally.

Gupta, V.↗

Low and High Frequency Models of Response Statistics of a Cylindrical Orthogrid Vehicle Panel to Acoustic Excitation

This presentation further develops the orthogrid vehicle panel work. Employed Hybrid Module capabilities to assess both low/mid frequency and high frequency models in the VA One simulation environment. The response estimates from three modeling approaches are compared to ground test measurements. Detailed Finite Element Model of the Test Article -Expect to capture both the global panel modes and the local pocket mode response, but at a considerable analysis expense (time & resources). A Composite Layered Construction equivalent global stiffness approximation using SEA -Expect to capture response of the global panel modes only. An SEA approximation using the Periodic Subsystem Formulation. A finite element model of a single periodic cell is used to derive the vibroacoustic properties of the entire periodic structure (modal density, radiation efficiency, etc. Expect to capture response at various locations on the panel (on the skin and on the ribs) with less analysis expense

Smith, Andrew↗

Experimental demonstration of static shape control

Results of a microprocessor-controlled implementation of static shape control using a specially constructed flexible beam facility are presented. The discussion covers the development of shape control algorithms, adaptation of the algorithms for use with finite element models, construction of a flexible beam, characterization and calibration of the facility, development of a finite element model for the beam, and the development of computer hardware and software. It is shown that feedback control yields better results than open-loop control, and that the use of more than two sensors in the control loop has little effect on the system performance.

Eldred, D.↗

Control of large flexible structures - An experiment on the NASA Mini-Mast facility

The output variance constraint controller design procedure is integrated with model reduction by modal cost analysis. A procedure is given for tuning MIMO controller designs to find the maximal rms performance of the actual system. Controller designs based on a finite-element model of the system are compared with controller designs based on an identified model (obtained using the Q-Markov Cover algorithm). The identified model and the finite-element model led to similar closed-loop performance, when tested in the Mini-Mast facility at NASA Langley.

Hsieh, Chen↗

Static Testing of an Inflatable/Rigidizable Hexapod Structure

Scientific application of large space antennas and telescopes is fostering the study of ultralightweight and inflatable structures. These so-called gossamer structures are tightly packaged for launch and subsequent deployment. Of particular interest is a hexapod configuration often found in vehicle ride simulators, ground-based telescopes, and antennas. One such configuration was modeled and tested at NASA Langley. Discrepancies between dynamic test data and predictions from a finite element model prompted additional static testing. This paper discusses static tests that were conducted to update stiffness parameters in the finite element model.

Berger, K. T.↗

The NASA/Industry Design Analysis Methods for Vibrations (DAMVIBS) Program - A government overview

An overview is presented of government contributions to the program called Design Analysis Methods for Vibrations (DAMV) which attempted to develop finite-element-based analyses of rotorcraft vibrations. NASA initiated the program with a finite-element modeling program for the CH-47D tandem-rotor helicopter. The DAMV program emphasized four areas including: airframe finite-element modeling, difficult components studies, coupled rotor-airframe vibrations, and airframe structural optimization. Key accomplishments of the program include industrywide standards for modeling metal and composite airframes, improved industrial designs for vibrations, and the identification of critical structural contributors to airframe vibratory responses. The program also demonstrated the value of incorporating secondary modeling details to improving correlation, and the findings provide the basis for an improved finite-element-based dynamics design-analysis capability.

Kvaternik, Raymond G.↗

Effects of finite element discretization and model simplification on calculations of ductile failure initiation

The finite element method is a scheme to discretize the infinite number of degrees of freedom in continuum-level problems down to a finite number of degrees of freedom. This discretization is done in conjunction with methods that also reduce the field differential equations to sets of algebraic ones that can be solved by arithmetical operations. Therefore, solutions attained by finite element models are approximations to the exact solutions of the field equations.

97 MATHEMATICS AND COMPUTING↗

Nonlinear fracture of concrete and ceramics

The nonlinear fracture process zones in an impacted unnotched concrete bend specimen, a prenotched ceramic bend specimen, and an unnotched ceramic/ceramic composite bend specimen were estimated through hybrid experimental numerical analysis. Aggregate bridging in concrete, particulate bridging in ceramics, and fiber bridging in ceramic/ceramic composite are modeled by Barenblatt-type cohesive zones which are incorporated into the finite-element models of the bend specimens. Both generation and propagation analyses are used to estimate the distribution of crack closure stresses in the nonlinear fracture process zones. The finite-element models are then used to simulate fracture tests consisting of rapid crack propagation in an impacted concrete bend specimen, and stable crack growth and strain softening in a ceramic and ceramic/ceramic composite bend specimens.

Kobayashi, Albert S.↗

Parameter identification of material constants in a composite shell structure

One of the basic requirements in engineering analysis is the development of a mathematical model describing the system. Frequently comparisons with test data are used as a measurement of the adequacy of the model. An attempt is typically made to update or improve the model to provide a test verified analysis tool. System identification provides a systematic procedure for accomplishing this task. The terms system identification, parameter estimation, and model correlation all refer to techniques that use test information to update or verify mathematical models. The goal of system identification is to improve the correlation of model predictions with measured test data, and produce accurate, predictive models. For nonmetallic structures the modeling task is often difficult due to uncertainties in the elastic constants. A finite element model of the shell was created, which included uncertain orthotropic elastic constants. A modal survey test was then performed on the shell. The resulting modal data, along with the finite element model of the shell, were used in a Bayes estimation algorithm. This permitted the use of covariance matrices to weight the confidence in the initial parameter values as well as confidence in the measured test data. The estimation procedure also employed the concept of successive linearization to obtain an approximate solution to the original nonlinear estimation problem.

Martinez, David R.↗

Modal Correlation of Complex Aerospace Joints Using Automated Variable Substitution

A critical task involved with being able to predict flight loads accurately in aerospace finite element models (FEMs) is the prior verification of the FEMs by conducting modal survey testing (MST). Experience comparing dynamic response of initial FEMs to MST data tends to demonstrate that FEMs can have unacceptable accuracy even when best modeling practices are followed. One inherent source of inaccuracy in linear dynamic FEMs is the modeling of nonlinear joints with mechanisms such as spherical bearings. These joints are usually designed to freely translate or rotate under the high levels of loading experienced in flight. Engineers who create linear FEMs conventionally model these joints without any stiffness in the mechanism degrees of freedom to meet this design intent. However, inaccuracy is observed during test validation of these FEMs, which usually relies on low-level force excitation orders of magnitude below flight load levels. This low-level modal test rarely overcomes the joint friction that is present, and thus the mechanism joints are able to react loads. This divide between the test results and the FEM creates a significant challenge to the engineer who is performing the correlation, in that the engineer has no basis for what stiffness value should be used to make the FEM match the test results. A compounding challenge is that complicated built-up aerospace structures commonly have multiple joints through a load path where each joint will “stick” and “slip” at different levels of force input. Explicitly matching the dynamics of a system containing these nonlinear mechanisms would require a nonlinear FEM, which is prohibitively costly for dynamic simulations of most aerospace systems. The objective of this paper is to present a workflow that can efficiently cycle through many iterations of a FEM, allowing a Monte Carlo style examination of the design space to identify candidate stiffness values for nonlinear mechanism joints. The outlined approach is specific to MSC Nastran and utilizes MSC Nastran’s symbolic substitution capabilities, coupled with the IMAT™ and Attune™ software packages developed by ATA Engineering, Inc. (ATA). The workflow is demonstrated with a case study from the correlation effort for The Boeing Company’s Crew Space Transportation (CST)-100 Starliner FEM. Keywords: Correlation, MSC Nastran, IMAT™, Attune™, Finite Element Model, Modal Testing

Correlation↗

Reduced Uncertainties in the Flutter Analysis of the Aerostructures Test Wing

Tuning the finite element model using measured data to minimize the model uncertainties is a challenging task in the area of structural dynamics. A test validated finite element model can provide a reliable flutter analysis to define the flutter placard speed to which the aircraft can be flown prior to flight flutter testing. Minimizing the difference between numerical and experimental results is a type of optimization problem. Through the use of the National Aeronautics and Space Administration Dryden Flight Research Center s (Edwards, California, USA) multidisciplinary design, analysis, and optimization tool to optimize the objective function and constraints; the mass properties, the natural frequencies, and the mode shapes are matched to the target data and the mass matrix orthogonality is retained. The approach in this study has been applied to minimize the model uncertainties for the structural dynamic model of the aerostructures test wing, which was designed, built, and tested at the National Aeronautics and Space Administration Dryden Flight Research Center. A 25-percent change in flutter speed has been shown after reducing the uncertainties

Pak, Chan-gi↗

Reduced Uncertainties in the Flutter Analysis of the Aerostructures Test Wing

Tuning the finite element model using measured data to minimize the model uncertainties is a challenging task in the area of structural dynamics. A test validated finite element model can provide a reliable flutter analysis to define the flutter placard speed to which the aircraft can be flown prior to flight flutter testing. Minimizing the difference between numerical and experimental results is a type of optimization problem. Through the use of the National Aeronautics and Space Administration Dryden Flight Research Center's (Edwards, California) multidisciplinary design, analysis, and optimization tool to optimize the objective function and constraints; the mass properties, the natural frequencies, and the mode shapes are matched to the target data, and the mass matrix orthogonality is retained. The approach in this study has been applied to minimize the model uncertainties for the structural dynamic model of the aerostructures test wing, which was designed, built, and tested at the National Aeronautics and Space Administration Dryden Flight Research Center. A 25 percent change in flutter speed has been shown after reducing the uncertainties.

Pak, Chan-Gi↗

Implementation of a Smeared Crack Band Model in a Micromechanics Framework

The smeared crack band theory is implemented within the generalized method of cells and high-fidelity generalized method of cells micromechanics models to capture progressive failure within the constituents of a composite material while retaining objectivity with respect to the size of the discretization elements used in the model. An repeating unit cell containing 13 randomly arranged fibers is modeled and subjected to a combination of transverse tension/compression and transverse shear loading. The implementation is verified against experimental data (where available), and an equivalent finite element model utilizing the same implementation of the crack band theory. To evaluate the performance of the crack band theory within a repeating unit cell that is more amenable to a multiscale implementation, a single fiber is modeled with generalized method of cells and high-fidelity generalized method of cells using a relatively coarse subcell mesh which is subjected to the same loading scenarios as the multiple fiber repeating unit cell. The generalized method of cells and high-fidelity generalized method of cells models are validated against a very refined finite element model.

Pineda, Evan J.↗

Progressive Failure of a Unidirectional Fiber-Reinforced Composite Using the Method of Cells: Discretization Objective Computational Results

The smeared crack band theory is implemented within the generalized method of cells and high-fidelity generalized method of cells micromechanics models to capture progressive failure within the constituents of a composite material while retaining objectivity with respect to the size of the discretization elements used in the model. An repeating unit cell containing 13 randomly arranged fibers is modeled and subjected to a combination of transverse tension/compression and transverse shear loading. The implementation is verified against experimental data (where available), and an equivalent finite element model utilizing the same implementation of the crack band theory. To evaluate the performance of the crack band theory within a repeating unit cell that is more amenable to a multiscale implementation, a single fiber is modeled with generalized method of cells and high-fidelity generalized method of cells using a relatively coarse subcell mesh which is subjected to the same loading scenarios as the multiple fiber repeating unit cell. The generalized method of cells and high-fidelity generalized method of cells models are validated against a very refined finite element model.

Pineda, Evan J.↗

A robust interface finite element formulation for modeling brittle material failure problems

Abstract Failure of many brittle materials and structures can be modeled using interface‐oriented finite elements combined with intrinsic cohesive zone models. The discontinuous Galerkin (DG) finite element method provides an innovative framework for modeling brittle crack propagation with zero‐thickness interface elements, which can accommodate extrinsic cohesive laws to avoid the artificial compliance required in intrinsic cohesive models. However, robust formulations and implementations of DG methods are critical in alleviating the well‐known convergence issues for both crack nucleation and propagation with reduced instability. This paper presents a robust interface element formulation by modifying the incomplete interior penalty Galerkin (IIPG) method, which successfully avoids the initial element interface penetration across elements that occurs prior to crack nucleation, and thereby greatly reduces the instability issue as cracks open. We further verified and validated our implementation by using a bar tension test and a beam fracturing benchmark. The robustness of our proposed interface element method was demonstrated by a micromechanics fiber/matrix debonding problem with 64 fibers embedded in a bulk matrix.

Liu, Ruijie↗

Optimum element density studies for finite-element thermal analysis of hypersonic aircraft structures

Different finite element models previously set up for thermal analysis of the space shuttle orbiter structure are discussed and their shortcomings identified. Element density criteria are established for the finite element thermal modelings of space shuttle orbiter-type large, hypersonic aircraft structures. These criteria are based on rigorous studies on solution accuracies using different finite element models having different element densities set up for one cell of the orbiter wing. Also, a method for optimization of the transient thermal analysis computer central processing unit (CPU) time is discussed. Based on the newly established element density criteria, the orbiter wing midspan segment was modeled for the examination of thermal analysis solution accuracies and the extent of computation CPU time requirements. The results showed that the distributions of the structural temperatures and the thermal stresses obtained from this wing segment model were satisfactory and the computation CPU time was at the acceptable level. The studies offered the hope that modeling the large, hypersonic aircraft structures using high-density elements for transient thermal analysis is possible if a CPU optimization technique was used.

Ko, William L.↗

Development of a laminate fatigue analysis

A fatigue analysis technique developed to predict damage growth in notched laminates is described. Features of the analysis include: criteria to relate matrix failure to cyclic stresses in and between plies; correlation of delamination growth with strain-release rate; and in-plane shear modulus change is related to cyclic shear stresses. A simplified finite element model is used to determine stresses in laminates that contain matrix damage. Failure criteria are integrated with the finite element model to form the fatigue analysis.

Roderick, G. L.↗

Vertical solidification of dendritic binary alloys

Three numerical techniques are employed to analyze the influence of thermosolutal convection on defect formation in directionally solidified (DS) alloys. The finite-element models are based on the Boussinesq approximation and include the plane-front model and two plane-front models incorporating special dendritic regions. In the second model the dendritic region has a time-independent volume fraction of liquid, and in the last model the dendritic region evolves as local conditions dictate. The finite-element models permit the description of nonlinear thermosolutal convection by treating the dendritic regions as porous media with variable porosities. The models are applied to lead-tin alloys including DS alloys, and severe segregation phenomena such as freckles and channels are found to develop in the DS alloys. The present calculations and the permeability functions selected are shown to predict behavior in the dendritic regions that qualitatively matches that observed experimentally.

Heinrich, J. C.↗