Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Element 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 1,243 records · Page 69

Finite Element Simulations of Two Vertical Drop Tests of F-28 Fuselage Sections

In March 2017, a vertical drop test of a forward fuselage section of a Fokker F-28 MK4000 aircraft was conducted as part of a joint NASA/FAA project to investigate the performance of transport aircraft under realistic crash conditions. In June 2017, a vertical drop test was conducted of a wing-box fuselage section of the same aircraft. Both sections were configured with two rows of aircraft seats, in a triple-double configuration. A total of ten Anthropomorphic Test Devices (ATDs) were secured in seats using standard lap belt restraints. The forward fuselage section was also configured with luggage in the cargo hold. Both sections were outfitted with two hat racks, each with added ballast mass. The drop tests were performed at the Landing and Impact Research facility located at NASA Langley Research Center in Hampton, Virginia. The measured impact velocity for the forward fuselage section was 346.8-in/s onto soil. The wing-box section was dropped with a downward facing pitch angle onto a sloping soil surface in order to create an induced forward acceleration in the airframe. The vertical impact velocity of the wing-box section was 349.2-in/s. A second objective of this project was to assess the capabilities of finite element simulations to predict the test responses. Finite element models of both fuselage sections were developed for execution in LS-DYNA(Registered Trademark), a commercial explicit nonlinear transient dynamic code. The models contained accurate representations of the airframe structure, the hat racks and hat rack masses, the floor and seat tracks, the luggage in the cargo hold for the forward section, and the detailed under-floor structure in the wing-box section. Initially, concentrated masses were used to represent the inertial properties of the seats, restraints, and ATD occupants. However, later simulations were performed that included finite element representations of the seats, restraints, and ATD occupants. These models were developed to more accurately replicate the seat loading of the floor and to enable prediction of occupant impact responses. Models were executed to generate analytical predictions of airframe responses, which were compared with test data to validate the model. Comparisons of predicted and experimental structural deformation and failures were made. Finally, predicted and experimental soil deformation and crater depths were also compared for both drop test configurations.

Jackson, Karen E.↗

Progressive Failure Studies of Stiffened Panels Subjected to Shear Loading

Experimental and analytical results are presented for progressive failure of stiffened composite panels with and without a notch and subjected to in plane shear loading well into their postbuckling regime. Initial geometric imperfections are included in the finite element models. Ply damage modes such as matrix cracking, fiber-matrix shear, and fiber failure are modeled by degrading the material properties. Experimental results from the test include strain field data from video image correlation in three dimensions in addition to other strain and displacement measurements. Results from nonlinear finite element analyses are compared with experimental data. Good agreement between experimental data and numerical results are observed for the stitched stiffened composite panels studied.

Ambur, Damodar R.↗

Calculation of AC Magnetizing Loss of ReBCO Superconducting Tapes Subjected to Applied Distorted Magnetic Fields

In recent years, there are fast-increasing concerns on the utilizations of superconducting rotating electrical machines in different application areas, such as ship propulsion systems, aircraft drivers, and wind turbine generators, since these machines exhibit the merits of high current density, compact design, high power density, light weight, high torque density as well as high efficiency. One of the main limitations in front of the vast use of superconducting tapes in the fabrication of electrical machineries is AC magnetizing loss when tapes are exposed to an external magnetic field, which can decrease the critical current density of wires, as well. In the literature, most of the research works have been done on calculation of the AC magnetizing loss under a pure external magnetic field, while in reality, magnetic flux lines in AC electrical machines are usually distorted with harmonics because of different reasons such as distorted leakage flux, distributed coils of a winding in several slots, cogging fields, mechanical faults, etc. Since these distorted fields contain harmonics, then in this paper, the AC magnetizing loss of superconducting tapes has been electromagnetically modeled and calculated when they are subjected to nonsinusoidally distorted external magnetic fields. The magnetic field dependency of critical current density has been considered in a proposed finite element model. The results have shown that the AC magnetizing loss increases significantly under a distorted applied field compared with a sinusoidal one. In addition, the loss increase depends on the harmonic content which would increases drastically with total harmonic distortion of the applied magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Optimizing the use of pressurized bladders for the assembly of HL-LHC MQXFB magnets

The use of pressurized bladders for stress control of superconducting magnets was firstly proposed at Lawrence Berkeley National Laboratory in the early 2000s. Since then, the so-called ‘bladders and keys’ procedure has become one of the reference techniques for the assembly of high-field accelerator magnets and demonstrators. Exploiting the advantages of this method is today of critical importance for Nb$_{3}$Sn-based accelerator magnets, whose production requires the preservation of tight stress targets in the superconducting coils to limit the effects of the strain sensitivity and brittleness of the conductor. The present manuscript reports on the results of an experimental campaign focused on the optimization of the ‘bladders and keys’ assembly process in the MQXFB quadrupoles. These 7.2 m long magnets shall be among the first Nb$_{3}$Sn cryomagnets to be installed in a particle accelerator as a part of the High Luminosity upgrade of the LHC. One of the main practical implications of the bladders technique, especially important when applied to long magnets like MQXFB, is that to insert the loading keys, the opening of a certain clearance in the support structure is required. The procedure used so far for MQXF magnets involved an overstress in the coils during bladder inflation. The work presented here shows that such an overshoot can be eliminated thanks to additional bladders properly positioned in the structure. This optimized method was validated in a short model magnet and in a full-length mechanical model, becoming the new baseline for the series production at CERN. Furthermore, the results are supported by numerical predictions using finite element models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Predicting microstructurally sensitive fatigue‐crack path in WE43 magnesium using high‐fidelity numerical modeling and three‐dimensional experimental characterization

Abstract Microstructurally small fatigue‐crack growth in polycrystalline materials is highly three‐dimensional due to sensitivity to local microstructural features (e.g., grains). One requirement for modeling microstructurally sensitive crack propagation is establishing the criteria that govern crack evolution, including crack deflection. Here, a high‐fidelity finite‐element modeling framework is used to assess the performance and validity of various crack‐growth criteria, including slip‐based metrics (e.g., fatigue‐indicator parameters), as potential criteria for predicting three‐dimensional crack paths in polycrystalline materials. The modeling framework represents cracks as geometrically explicit discontinuities and involves voxel‐based remeshing, mesh‐gradation control, and a crystal‐plasticity constitutive model. The predictions are compared to experimental measurements of WE43 magnesium samples subject to fatigue loading, for which three‐dimensional grain structures and fatigue‐crack surfaces were measured post‐mortem using near‐field high‐energy x‐ray diffraction microscopy and x‐ray computed tomography. Findings from this work are expected to improve the predictive capabilities of simulations involving microstructurally small fatigue‐crack growth in polycrystalline materials.

Engineering↗

A coupled transport and solid mechanics formulation with improved reaction kinetics parameters for modeling oxidation and decomposition in a uranium hydride bed.

Modeling of reacting flows in porous media has become particularly important with the increased interest in hydrogen solid-storage beds. An advanced type of storage bed has been proposed that utilizes oxidation of uranium hydride to heat and decompose the hydride, releasing the hydrogen. To reduce the cost and time required to develop these systems experimentally, a valid computational model is required that simulates the reaction of uranium hydride and oxygen gas in a hydrogen storage bed using multiphysics finite element modeling. This SAND report discusses the advancements made in FY12 (since our last SAND report SAND2011-6939) to the model developed as a part of an ASC-P&EM project to address the shortcomings of the previous model. The model considers chemical reactions, heat transport, and mass transport within a hydride bed. Previously, the time-varying permeability and porosity were considered uniform. This led to discrepancies between the simulated results and experimental measurements. In this work, the effects of non-uniform changes in permeability and porosity due to phase and thermal expansion are accounted for. These expansions result in mechanical stresses that lead to bed deformation. To describe this, a simplified solid mechanics model for the local variation of permeability and porosity as a function of the local bed deformation is developed. By using this solid mechanics model, the agreement between our reacting bed model and the experimental data is improved. Additionally, more accurate uranium hydride oxidation kinetics parameters are obtained by fitting the experimental results from a pure uranium hydride oxidation measurement to the ones obtained from the coupled transport-solid mechanics model. Finally, the coupled transport-solid mechanics model governing equations and boundary conditions are summarized and recommendations are made for further development of ARIA and other Sandia codes in order for them to sufficiently implement the model.

Salloum, Maher N.↗

Automated design optimization of supersonic airplane wing structures under dynamic constraints.

The problems of the preliminary and first level detail design of supersonic aircraft wings are stated as mathematical programs and solved using automated optimum design techniques. The problem is approached in two phases: the first is a simplified equivalent plate model in which the envelope, plan form and structural parameters are varied to produce a design, the second is a finite element model with fixed configuration in which the material distribution is varied. Constraints include flutter, aeroelastically computed stresses and deflections, natural frequency and a variety of geometric limitations. The Phase I objective is a combination of weight and drag while Phase II is a weight minimization.

Fox, R. L.↗

A structural analysis of an ocean going patrol boat subjected to planning loads

A static structural analysis of an ocean going patrol vessel subjected to hydrodynamic planning loads is discussed. The analysis required the development of a detailed model that included hull plating, five structural bulkheads, longitudinal and transverse stiffners, and a coarse representation of the superstructure. The finite element model was developed from fabrication drawings using the Navy computer aided design system. Various stress and displacement contours are shown for the entire hull. Because several critical areas appeared to be overstressed, these areas were remeshed for detail and are presented for completeness.

Clark, James H.↗

Vibration, performance, flutter and forced response characteristics of a large-scale propfan and its aeroelastic model

An investigation of the vibration, performance, flutter, and forced response of the large-scale propfan, SR7L, and its aeroelastic model, SR7A, has been performed by applying available structural and aeroelastic analytical codes and then correlating measured and calculated results. Finite element models of the blades were used to obtain modal frequencies, displacements, stresses and strains. These values were then used in conjunction with a 3-D, unsteady, lifting surface aerodynamic theory for the subsequent aeroelastic analyses of the blades. The agreement between measured and calculated frequencies and mode shapes for both models is very good. Calculated power coefficients correlate well with those measured for low advance ratios. Flutter results show that both propfans are stable at their respective design points. There is also good agreement between calculated and measured blade vibratory strains due to excitation resulting from yawed flow for the SR7A propfan. The similarity of structural and aeroelastic results show that the SR7A propfan simulates the SR7L characteristics.

August, Richard↗

Vibration, performance, flutter and forced response characteristics of a large-scale propfan and its aeroelastic model

An investigation of the vibration, performance, flutter, and forced response of the large-scale propfan, SR7L, and its aeroelastic model, SR7A, has been performed by applying available structural and aeroelastic analytical codes and then correlating measured and calculated results. Finite element models of the blades were used to obtain modal frequencies, displacements, stresses and strains. These values were then used in conjunction with a 3-D, unsteady, lifting surface aerodynamic theory for the subsequent aeroelastic analyses of the blades. The agreement between measured and calculated frequencies and mode shapes for both models is very good. Calculated power coefficients correlate well with those measured for low advance ratios. Flutter results show that both propfans are stable at their respective design points. There is also good agreement between calculated and measured blade vibratory strains due to excitation resulting from yawed flow for the SR7A propfan. The similarity of structural and aeroelastic results show that the SR7A propfan simulates the SR7L characteristics.

August, Richard↗

Model Of Bearing With Hydrostatic Damper

Improved mathematical model of rotational and vibrational dynamics of bearing package in turbopump incorporates effects of hydrostatic damper. Part of larger finite-element model representing rotational and vibrational dynamics of rotor and housing of pump. Includes representations of deadband and nonlinear stiffness and damping of ball bearings, nonlinear stiffness and damping of hydrostatic film, and stiffness of bearing support. Enables incorporation of effects of hydrostatic damper into overall rotor-dynamic mathematical model without addition of mathematical submodel of major substructure.

Goggin, David G.↗

The effect of friction in the hold down post spherical bearings on hold down post loads

The effect of friction at the connection of the Solid Rocket Booster (SRB) aft skirt and the mobile launch platform (MLP) hold down posts was analyzed. A simplified model of the shuttle response during the Space Shuttle Main Engine (SSME) buildup was constructed. The model included the effect of stick-slip friction for the rotation of the skirt about the spherical bearing. Current finite element models assume the joint is completely frictionless in rotation and therefore no moment is transferred between the skirt and the hold down posts. The model was partially verified against test data and preliminary parameter studies were performed. The parameter studies indicated that the coefficient of friction strongly influenced the moment on the hold down posts. The coefficient of friction had little effect on hold down post vertical loads, however. Further calibration of the model is necessary before the effect of friction on the hold down post horizontal loads can be analyzed.

Richardson, James A.↗

Enhancements to NURBS-Based FEA Airfoil Modeler: SABER

NURBS (Non-Uniform Rational B-Splines) have become a common way for CAD programs to fit a smooth surface to discrete geometric data. This concept has been extended to allow for the fitting of analysis data in a similar manner and "attaching" the analysis data to the geometric definition of the structure. The "attaching" of analysis data to the geometric definition allows for a more seamless sharing of data between analysis disciplines. NURBS have become a useful tool in the modeling of airfoils. The use of NURBS has allowed for the development of software that easily and consistently generates plate finite element models of the midcamber surface of a given airfoil. The resulting displacements can then be applied to the original airfoil surface and the deformed shape calculated.

Saleeb, A. F.↗

Orion Parachute Riser Cutter Development

This paper presents the tests and analytical approach used on the development of a steel riser cutter for the CEV Parachute Assembly System (CPAS) used on the Orion crew module. Figure 1 shows the riser cutter and the steel riser bundle which consists of six individual cables. Due to the highly compressed schedule, initial unavailability of the riser material and the Orion Forward Bay mechanical constraints, JSC primarily relied on a combination of internal ballistics analysis and LS-DYNA simulation for this project. Various one dimensional internal ballistics codes that use standard equation of state and conservation of energy have commonly used in the development of CAD devices for initial first order estimates and as an enhancement to the test program. While these codes are very accurate for propellant performance prediction, they usually lack a fully defined kinematic model for dynamic predictions. A simple piston device can easily and accurately be modeled using an equation of motion. However, the accuracy of analytical models is greatly reduced on more complicated devices with complex external loads, nonlinear trajectories or unique unlocking features. A 3D finite element model of CAD device with all critical features included can vastly improve the analytical ballistic predictions when it is used as a supplement to the ballistic code. During this project, LS-DYNA structural 3D model was used to predict the riser resisting load that was needed for the ballistic code. A Lagrangian model with eroding elements shown in Figure 2 was used for the blade, steel riser and the anvil. The riser material failure strain was fine tuned by matching the dent depth on the anvil with the actual test data. LS-DYNA model was also utilized to optimize the blade tip design for the most efficient cut. In parallel, the propellant type and the amount were determined by using CADPROG internal ballistics code. Initial test results showed a good match with LS-DYNA and CADPROG simulations. Final paper will present a detailed roadmap from initial ballistic modeling and LS-DYNA simulation to the performance testing. Blade shape optimization study will also be presented.

Oguz, Sirri↗

Inverse Force Determination on a Small Scale Launch Vehicle Model Using a Dynamic Balance

A launch vehicle can experience large unsteady aerodynamic forces in the transonic regime that, while usually only lasting for tens of seconds during launch, could be devastating if structural components and electronic hardware are not designed to account for them. These aerodynamic loads are difficult to experimentally measure and even harder to computationally estimate. The current method for estimating buffet loads is through the use of a few hundred unsteady pressure transducers and wind tunnel test. Even with a large number of point measurements, the computed integrated load is not an accurate enough representation of the total load caused by buffeting. This paper discusses an attempt at using a dynamic balance to experimentally determine buffet loads on a generic scale hammer head launch vehicle model tested at NASA Ames Research Center's 11' x 11' transonic wind tunnel. To use a dynamic balance, the structural characteristics of the model needed to be identified so that the natural modal response could be and removed from the aerodynamic forces. A finite element model was created on a simplified version of the model to evaluate the natural modes of the balance flexures, assist in model design, and to compare to experimental data. Several modal tests were conducted on the model in two different configurations to check for non-linearity, and to estimate the dynamic characteristics of the model. The experimental results were used in an inverse force determination technique with a psuedo inverse frequency response function. Due to the non linearity, the model not being axisymmetric, and inconsistent data between the two shake tests from different mounting configuration, it was difficult to create a frequency response matrix that satisfied all input and output conditions for wind tunnel configuration to accurately predict unsteady aerodynamic loads.

testing↗

Aeroelastic Modeling of X-56A Stiff-Wing Configuration Flight Test Data

Aeroelastic stability and control derivatives for the X-56A Multi-Utility Technology Testbed (MUTT), in the stiff-wing configuration, were estimated from flight test data using the output-error method. Practical aspects of the analysis are discussed. The orthogonal phase-optimized multisine inputs provided excellent data information for aeroelastic modeling. Consistent parameter estimates were determined using output error in both the frequency and time domains. The frequency domain analysis converged faster and was less sensitive to starting values for the model parameters, which was useful for determining the aeroelastic model structure and obtaining starting values for the time domain analysis. Including a modal description of the structure from a finite element model reduced the complexity of the estimation problem and improved the modeling results. Effects of reducing the model order on the short period stability and control derivatives were investigated.

Grauer, Jared A.↗

Relevence of 3D Simulations and Sandwich Core Topology for the Modeling of Honeycomb Core Sandwich Panels Undergoing Mode I Delamination

A recently developed cohesive zone traction-separation law, which includes the effects of fiber bridging in a novel way, is extended from 2D to 3D. The proposed cohesive model is applied to low fidelity (i.e. homogenized core representation) and high fidelity (i.e. directly accounting for the core topology) finite element models of a composite panel comprised of carbon fiber reinforced plastic facesheets and a honeycomb sandwich core. This enables the investigation of 2D to 3D parameter transferability, width-dependent effects such as thumbnail-shaped crack growth, and the verification of plane strain / plane stress assumptions. A pronounced curvature of the initial interface-related crack front is observed, while the bridging-related crack front is straight. Furthermore, it is found that the cohesive parameters can easily be transferred from 2D to 3D under plane stress assumptions, but not under plane strain assumptions. The numerical predictions are compared to experimental load-displacement and R-curves.

Composites↗

Railcar Dynamics Model of the ENSA/DOE Multimodal Transportation Campaign Rail Conveyance System

Railcar dynamics models can be applied to calculate loads transmitted to used nuclear fuel (UNF) during normal conditions of transportation (NCT). This work details the development of a railcar dynamics model for the Equipos Nucleares of Spain-U.S. Department of Energy (ENSA/DOE) multimodal transportation campaign. The model was developed using NUCARS and represents the Kasgro 12-axle railcar and the ENSA ENUN-32P cask and cradle used during the campaign. The dynamics of the cask and cradle are captured by modeling them as a single degree-of-freedom (SDOF) system with optimized stiffness and damping parameters. The model was tuned using a statistical design method and manual parameter adjustment, and validated by comparing the predicted response to data collected from captive track tests. The SDOF cask and cradle model simplifies downstream analyses used to predict loads on UNF. Previous methods fed the railcar deck motion into computationally expensive finite element models of the cask and cradle; capturing the dynamics of the cask and cradle in the railcar model bypasses this step. The model provides information useful for understanding the shock and vibration environment of UNF during rail transportation, and the agreement achieved between model and test provides confidence that the methods used in this work could be used for other vehicle configurations such as the ATLAS railcar.

Kadooka, Kevin↗