Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element analysis”

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 613 records · Page 34

Ground shake test of the UH-60A helicopter airframe and comparison with NASTRAN finite element model predictions

Sikorsky Aircraft, together with the other major helicopter airframe manufacturers, is engaged in a study to improve the use of finite element analysis to predict the dynamic behavior of helicopter airframes, under a rotorcraft structural dynamics program called DAMVIBS (Design Analysis Methods for VIBrationS), sponsored by the NASA-Langley. The test plan and test results are presented for a shake test of the UH-60A BLACK HAWK helicopter. A comparison is also presented of test results with results obtained from analysis using a NASTRAN finite element model.

Howland, G. R.↗

Effect of microcracks on the thermal expansion of composite laminates

A finite element analysis was used to quantitatively predict the effect of matrix microcracks in the 90 deg plies of graphite/epoxy composites on the coefficient of thermal expansion in the 0 deg direction, alpha (y) (perpendicular to the cracks). Results were generated for (0m/90n)s, (0/+ or 45/90)s and (0/90/ + or - 45)s graphite/epoxy laminate configurations. Analytical predictions were compared with experimental results for the two quasi-isotropic laminate configurations. Both analytical and experimental results showed that microcracks reduced the effective stiffness of the 90 deg plies, thus causing the laminates, thermal response to be more like that of a (0) laminate. The change in alpha(y) was a function of lamina material properties, microcrack density, fiber orientation, and stacking sequence. A combination of classical lamination theory and finite element analysis was used to predict the effect of microcracks in both the 90 deg and 0 deg plies. Analytical results showed that the addition of microcracks in the 0 deg plies do affect alpha(y), but to a lesser extent than those in the 90 deg plies.

Bowles, D. E.↗

Effect of Ex Vivo Ionizing Radiation on Static and Fatigue Properties of Mouse Vertebral Bodies

For a variety of medical and scientific reasons, human bones can be exposed to a wide range of ionizing radiation levels. In vivo radiation therapy (0.05 kGy) is used in cancer treatment, and ex vivo irradiation (25-35 kGy) is used to sterilize bone allografts. Ionizing radiation in these applications has been shown to increase risk of fracture, decrease bone quality and degrade collagen integrity. Past studies have investigated the deleterious effects of radiation on cortical or trabecular bone specimens individually, but to date no studies have examined whole bones containing both cortical and trabecular tissue. Furthermore, a clear relationship between the dose and the mechanical and biochemical response of bone's extracellular matrix has yet to be established for doses ranging from cancer therapy to allograft sterilization (0.05-35 kGy). To gain insight into these issues, we conducted an ex vivo radiation study to investigate non-cellular (i.e. matrix) effects of ionizing radiation dose on vertebral whole bone mechanical properties, over a range of radiation doses (0.05-35 kGy), with a focus on any radiation-induced changes in collagen. With underlying mechanisms of action in mind, we hypothesized that any induced reductions in mechanical properties would be associated with changes in collagen integrity. METHODS: 20-week old female mice were euthanized and the lumbar spine was dissected using IACUC approved protocols. The lumbar vertebrae (L1- S1) were extracted from the spine via cuts through adjacent intervertebral discs, and the endplates, posterior processes, surrounding musculature, and soft tissues were removed (approx. 1.5mm diameter, approx. 2mm height). Specimens were randomly assigned to one of five groups for ex vivo radiation exposure: x-ray irradiation at 0.05, 1, 17, or 35 kGy, or a 0 kGy control. Following irradiation, the vertebrae were imaged using microcomputed tomography (micro-CT) and then subjected to either monotonic compressive loading to failure or uniform cyclic compressive loading. During cyclic testing, samples were loaded in force control to a force level that corresponded to a strain of 0.46%, as determined in advance by a linearly elastic micro-CT-based finite element analysis for each specimen. Tests were stopped at imminent fracture, defined as a rapid increase in strain. The main outcome for the monotonic test was the strength (maximum force); for cyclic testing it was the fatigue life (log of the number of cycles of loading at imminent failure). A fluorometric assay was used on the S1 vertebrae to measure the number of non-enzymatic collagen crosslinks[4]. A one-way ANOVA was performed on mechanical properties and collagen crosslinks; means were compared with controls using Dunnett's method, with a Tukey-Kramer post-hoc analysis when significance was found (p < or = 0.05). RESULTS: Compared to the unirradiated control group, the concentration of non-enzymatic collagen crosslinks was significantly increased for all irradiated groups (p < 0.0001), and being higher by at least 50% (Figure 1a). By contrast, the radiation effects on the collagen were only evident at the higher doses. For irradiation exposures of 17 kGy or more, strength decreased substantially as the radiation level was increased, but no effect was evident below 17 kGy (Figure 1b). There was no significant change in the stiffness or maximum displacement for any radiation dose (p>0.05). The finite element analysis prescribed force level for cyclic loading exceeded the measured (monotonic) strength of the 17 and 35 kGy irradiated groups (mean +/- SD, 20.6 +/- 5.6 N; 13.2 +/- 3.7 N, respectively) and therefore these groups were eliminated from the fatigue study. The fatigue life for the 0.05 and 1 kGy groups were similar to each other and were not statistically significantly different from the control group (Figure 1c).

radiation↗

A methodology to predict damage initiation, damage growth and residual strength in titanium matrix composites

In this research, a methodology to predict damage initiation, damage growth, fatigue life, and residual strength in titanium matrix composites (TMC) is outlined. Emphasis was placed on micromechanics-based engineering approaches. Damage initiation was predicted using a local effective strain approach. A finite element analysis verified the prevailing assumptions made in the formulation of this model. Damage growth, namely, fiber-bridged matrix crack growth, was evaluated using a fiber bridging (FB) model which accounts for thermal residual stresses. This model combines continuum fracture mechanics and micromechanics analyses yielding stress-intensity factor solutions for fiber-bridged matrix cracks. It is assumed in the FB model that fibers in the wake of the matrix crack are idealized as a closure pressure, and an unknown constant frictional shear stress is assumed to act along the debond length of the bridging fibers. This frictional shear stress was used as a curve fitting parameter to the available experimental data. Fatigue life and post-fatigue residual strength were predicted based on the axial stress in the first intact 0 degree fiber calculated using the FB model and a three-dimensional finite element analysis.

Bakuckas, J. G., Jr.↗

Intelligent resolution: Integrating Cryo-EM with AI-driven multi-resolution simulations to observe the severe acute respiratory syndrome coronavirus-2 replication-transcription machinery in action

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) replication transcription complex (RTC) is a multi-domain protein responsible for replicating and transcribing the viral mRNA inside a human cell. Attacking RTC function with pharmaceutical compounds is a pathway to treating COVID-19. Conventional tools, e.g., cryo-electron microscopy and all-atom molecular dynamics (AAMD), do not provide sufficiently high resolution or timescale to capture important dynamics of this molecular machine. Consequently, we develop an innovative workflow that bridges the gap between these resolutions, using mesoscale fluctuating finite element analysis (FFEA) continuum simulations and a hierarchy of AI-methods that continually learn and infer features for maintaining consistency between AAMD and FFEA simulations. We leverage a multi-site distributed workflow manager to orchestrate AI, FFEA, and AAMD jobs, providing optimal resource utilization across HPC centers. Our study provides unprecedented access to study the SARS-CoV-2 RTC machinery, while providing general capability for AI-enabled multi-resolution simulations at scale.

Trifan, Anda↗

Model reduction in the physical coordinate system

In the dynamics modeling of a flexible structure, finite element analysis employs reduction techniques, such as Guyan's reduction, to remove some of the insignificant physical coordinates, thus producing a dynamics model that has smaller mass and stiffness matrices. But this reduction is limited in the sense that it removes certain degrees of freedom at a node points themselves in the model. From the standpoint of linear control design, the resultant model is still too large despite the reduction. Thus, some form of the model reduction is frequently used in control design by approximating a large dynamical system with a fewer number of state variables. However, a problem arises from the placement of sensors and actuators in the reduced model, because a model usually undergoes, before being reduced, some form of coordinate transformations that do not preserve the physical meanings of the states. To correct such a problem, a method is developed that expresses a reduced model in terms of a subset of the original states. The proposed method starts with a dynamic model that is originated and reduced in finite element analysis. Then the model is converted to the state space form, and reduced again by the internal balancing method. At this point, being in the balanced coordinate system, the states in the reduced model have no apparent resemblance to those of the original model. Through another coordinate transformation that is developed, however, this reduced model is expressed by a subset of the original states.

Yae, K. Harold↗

New Representation of Bearings in LS-DYNA

Non-linear, dynamic, finite element analysis is used in various engineering disciplines to evaluate high-speed, dynamic impact and vibration events. Some of these applications require connecting rotating to stationary components. For example, bird impacts on rotating aircraft engine fan blades are a common analysis performed using this type of analysis tool. Traditionally, rotating machines utilize some type of bearing to allow rotation in one degree of freedom while offering constraints in the other degrees of freedom. Most times, bearings are modeled simply as linear springs with rotation. This is a simplification that is not necessarily accurate under the conditions of high-velocity, high-energy, dynamic events such as impact problems. For this reason, it is desirable to utilize a more realistic non-linear force-deflection characteristic of real bearings to model the interaction between rotating and non-rotating components during dynamic events. The present work describes a rolling element bearing model developed for use in non-linear, dynamic finite element analysis. This rolling element bearing model has been implemented in LS-DYNA as a new element, *ELEMENT_BEARING.

bearings↗

Contact pressure analysis to allow improved design for the clearing plate in a biomass comminution system

For biomass processing systems, feedstocks with small but equal sizes can improve flowability and thermochemical outputs. The Crumbler rotary shear system was designed to achieve such ideal feedstocks to reduce costs and energy throughout the biomass supply chain. Improving the wear resistance of major components in this system is important to decrease the cost. Accelerating the feedstock flow through the rotary shear machine plays an equivalent role of improving its efficiency. An analytical analysis and a finite-element analysis of the stationary clearing plate are presented in this work. Here, these tools are used to investigate the effect of contact surface curvature on the impact pressure. The optimized new clearing plate design is expected to improve the component’s lifetime and wood chips flowability. The field tests confirm that the trend of the new clearing plate design guided by FEA. The optimization method, model verification, and validation experience applied in this work can also be applied to other static components with similar simple contact wear problems.

09 BIOMASS FUELS↗

NASTRAN data generation of helicopter fuselages using interactive graphics

The development and implementation of a preprocessor system for the finite element analysis of helicopter fuselages is described. The system utilizes interactive graphics for the generation, display, and editing of NASTRAN data for fuselage models. It is operated from an IBM 2250 cathode ray tube (CRT) console driven by an IBM 370/145 computer. Real time interaction plus automatic data generation reduces the nominal 6 to 10 week time for manual generation and checking of data to a few days. The interactive graphics system consists of a series of satellite programs operated from a central NASTRAN Systems Monitor. Fuselage structural models including the outer shell and internal structure may be rapidly generated. All numbering systems are automatically assigned. Hard copy plots of the model labeled with GRID or elements ID's are also available. General purpose programs for displaying and editing NASTRAN data are included in the system. Utilization of the NASTRAN interactive graphics system has made possible the multiple finite element analysis of complex helicopter fuselage structures within design schedules.

Sainsbury-Carter, J. B.↗

Implementation of design sensitivity analysis with existing finite element codes

A numerical method is presented to implement structural design sensitivity analysis theory, using the versatility and convenience of existing finite element structural analysis programs. Design variables such as thickness and cross-sectional areas of components of individual members and built-up structures are considered. Structural performance functionals considered include displacement and stress. The method is also applicable for eigenvalue problem design sensitivity analysis. It is shown that calculations can be carried out outside existing finite element codes, using postprocessing data only. Thus, design sensitivity analysis software does not have to be embedded in an existing finite element code. Feasability of the method is shown through analysis of several problems, including a built-up structure. Accurate design sensitivity results are obtained without the uncertainty of numerical accuracy associated with selection of finite difference perturbations.

Choi, Kyung K.↗

Finite element engine blade structural optimization

The Structural Tailoring of Engine Blades (STAEBL) computer program was developed to perform engine fan blade numerical optimizations. These blade optimizations seek a minimum weight or cost design that satisfies realistic blade design constraints, by tuning one to twenty design variables. The STAEBL system has been generalized to include both fan and compressor blade numerical optimizations. The system analyses have been significantly improved through the inclusion of an efficient plate finite element analysis for blade stress and frequency determinations. Additionally, a finite element based approximate severe foreign object damage (FOD) analysis has been included. The new FOD analysis gives very accurate estimates of the full nonlinear bird ingestion solution. Optimizations of fan and compressor blades have been performed using the system, showing significant cost and weight reductions, while comparing very favorably with refined design validation procedures.

Brown, K. W.↗

Charges on a suspended silicon nitride membrane under a high-energy electron beam

Thin silicon nitride (SiN x ) membranes are widely used in gas and liquid phase transmission electron microscopy (TEM) and as phase plates to enhance imaging contrast. SiN x contains trap sites for both positive and negative charges, which can be manipulated by high-energy electron irradiation, external potential biasing, or light exposure. Charge accumulation on the membrane can significantly affect in situ TEM processes, including chemical and electrochemical reactions, nanoparticle dynamics, and catalytic activity, or introduce unwanted phase shifts when used as a phase plate. Here, in this study, charge accumulation on suspended SiN x membranes was investigated using off-axis electron holography combined with model-free charge analysis, supported by custom finite element analysis (FEA) simulations. An average residual positive charge density of approximately 2.8 × 10 −4 C m −2 was measured. Localized and stable regions of both positive and negative charges were identified on the membrane. The global positive and localized positive/negative charges give rise to strong electric fields and electroosmotic slip velocities at the membrane surface, which are sufficient to induce non-Brownian particle behavior and directional fluid flow, offering a physical explanation for previously observed anomalies in particle dynamics, nucleation, and growth during gas and liquid phase TEM experiments. These results provide a benchmark for understanding charge behavior at SiNx interfaces in gas and liquid phase TEM. Furthermore, the FEA simulations establish a framework for future investigations into charge distribution, electrostatic potentials, and electrical double layers at solid–liquid interfaces, particularly in complex geometries and chemically dynamic environments.

Suspended silicon↗

Sierra/SD - Its2Sierra - User's Manual - 5.20

The Integrated Tiger Series (ITS) generates a database containing energy deposition data. This data, when stored on an Exodus file, is not typically suitable for analysis within SierraMechanics for finite element analysis. The its2sierra tool maps data from the ITS database to the Sierra database. This document provides information on the usage of its2sierra.

97 MATHEMATICS AND COMPUTING↗

Finite element stress analysis of polymers at high strains

A numerical analysis is presented for the problem of a flat rectangular rubber membrane with a circular rigid inclusion undergoing high strains due to the action of an axial load. The neo-hookean constitutive equations are introduced into the general purpose TITUS program by means of equivalent hookean constants and initial strains. The convergence is achieved after a few iterations. The method is not limited to any specific program. The results are in good agreement with those of a company sponsored photoelastic stress analysis. The theoretical and experimental deformed shapes also agree very closely with one another. For high strains it is demonstrated that using the conventional HOOKE law the stress concentration factor obtained is unreliable in the case of rubberlike material.

Durand, M.↗

Nonlinear crack analysis with finite elements

The application of finite element techniques to the analytic representation of the nonlinear behavior of arbitrary two-dimensional bodies containing cracks is discussed. Specific methods are proposed using which it should be possible to obtain information concerning: the description of the maximum, minimum, and residual near-tip stress and strain fields; the effects of crack closure on the near-tip behavior of stress and strain fields during cyclic loading into the plastic range; the stress-strain and displacement field behavior associated with a nonstationary crack; and the effects of large rotation near the crack tip.

Armen, H., Jr.↗

A triangular thin shell finite element: Linear analysis

The formulation of the linear stiffness matrix for a doubly-curved triangular thin shell element, using a modified potential energy principle, is described. The strain energy component of the potential energy is expressed in terms of displacements and displacement gradients by use of consistent Koiter strain-displacement equations. The element inplane and normal displacement fields are approximated by complete cubic polynomials. The interelement displacement admissibility conditions are met in the global representation by imposition of constraint conditions on the interelement boundaries; the constraints represent the modification of the potential energy. Errors due to the nonzero strains under rigid body motion are shown to be of small importance for practical grid refinements through performance of extensive comparison analyses.

Thomas, G. R.↗

A numerically efficient finite element hydroelastic analysis. Volume 2: Implementation in NASTRAN, part 1

Details are presented of the implementation of the new formulation into NASTRAN including descriptions of the DMAP statements required for conversion of the program and details pertaining to problem definition and bulk data considerations. Details of the current 1/8-scale space shuttle external tank mathematical model, numerical results and analysis/test comparisons are also presented. The appendices include a description and listing of a FORTRAN program used to develop harmonic transformation bulk data (multipoint constraint statements) and sample bulk data information for a number of hydroelastic problems.

Coppolino, R. N.↗