Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “static condensation”

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 19 records

An easily implemented static condensation method for structural sensitivity analysis

A black-box approach to static condensation for sensitivity analysis is presented with illustrative examples of a cube and a car structure. The sensitivity of the structural response with respect to joint stiffness parameter is calculated using the direct method, forward-difference, and central-difference schemes. The efficiency of the various methods for identifying joint stiffness parameters from measured static deflections of these structures is compared. The results indicate that the use of static condensation can reduce computation times significantly and the black-box approach is only slightly less efficient than the standard implementation of static condensation. The ease of implementation of the black-box approach recommends it for use with general-purpose finite element codes that do not have a built-in facility for static condensation.

Gangadharan, S. N.↗

An exact zooming method for finite element analyses

An exact zooming technique which employs static condensation and exact structural reanalysis methods was developed. Successive application of static condensation reduces the system to one that is only associated with the degrees of freedom (dof) of the original model. Application of an exact static reanalysis technique permits the displacements at the dof of the original model that are contained in the zoomed portion of the structure to be obtained first. The response external to the zoom, as well as the response of additional dof within various levels of zooming, is computed. With the triangular factor of the stiffness matrix of the original system available, this approach involves only the solution of a system of equations of small order.

Hirai, I.↗

An independent refinement and integration procedure in multiregion finite element analysis

An independent refinement and integration procedure is developed to couple together independently modeled (global and local) regions in a single analysis. The models can have different levels of refinement and along the interface between them the finite element nodes need not coincide with one another. In the local model all the nodes except the nodes at the interface are statically condensed and the reduced stiffness matrix is obtained. For this static condensation a modified frontal solution technique is employed. A spline interpolation function that satisfies the linear isotropic plate bending differential equation is used to relate the local model interface nodal displacements to the global model interface displacements. The proposed independent refinement and integration procedure is evaluated by applying it to two- and three-dimensional cases involving inplane and out-of-plane deformation. The procedure yielded very accurate results for all the examples studied.

Krishnamurthy, T.↗

Space station static and dynamic analyses using parallel methods

Algorithms for high-performance parallel computers are applied to perform static analyses of large-scale Space Station finite-element models (FEMs). Several parallel-vector algorithms under development at NASA Langley are assessed. Sparse matrix solvers were found to be more efficient than banded symmetric or iterative solvers for the static analysis of large-scale applications. In addition, new sparse and 'out-of-core' solvers were found superior to substructure (superelement) techniques which require significant additional cost and time to perform static condensation during global FEM matrix generation as well as the subsequent recovery and expansion. A method to extend the fast parallel static solution techniques to reduce the computation time for dynamic analysis is also described. The resulting static and dynamic algorithms offer design economy for preliminary multidisciplinary design optimization and FEM validation against test modes. The algorithms are being optimized for parallel computers to solve one-million degrees-of-freedom (DOF) FEMs. The high-performance computers at NASA afforded effective software development, testing, efficient and accurate solution with timely system response and graphical interpretation of results rarely found in industry. Based on the author's experience, similar cooperation between industry and government should be encouraged for similar large-scale projects in the future.

Gupta, V.↗

Rapid iterative reanalysis for automated design

A method for iterative reanalysis in automated structural design is presented for a finite-element analysis using the direct stiffness approach. A basic feature of the method is that the generalized stiffness and inertia matrices are expressed as functions of structural design parameters, and these generalized matrices are expanded in Taylor series about the initial design. Only the linear terms are retained in the expansions. The method is approximate because it uses static condensation, modal reduction, and the linear Taylor series expansions. The exact linear representation of the expansions of the generalized matrices is also described and a basis for the present method is established. Results of applications of the present method to the recalculation of the natural frequencies of two simple platelike structural models are presented and compared with results obtained by using a commonly applied analysis procedure used as a reference. In general, the results are in good agreement. A comparison of the computer times required for the use of the present method and the reference method indicated that the present method required substantially less time for reanalysis. Although the results presented are for relatively small-order problems, the present method will become more efficient relative to the reference method as the problem size increases. An extension of the present method to static reanalysis is described, ana a basis for unifying the static and dynamic reanalysis procedures is presented.

Bhatia, K. G.↗

Computational procedures for evaluating the sensitivity derivatives of vibration frequencies and Eigenmodes of framed structures

Computational procedures are presented for evaluating the sensitivity derivatives of the vibration frequencies and eigenmodes of framed structures. Both a displacement and a mixed formulation are used. The two key elements of the computational procedure are: (a) Use of dynamic reduction techniques to substantially reduce the number of degrees of freedom; and (b) Application of iterative techniques to improve the accuracy of the derivatives of the eigenmodes. The two reduction techniques considered are the static condensation and a generalized dynamic reduction technique. Error norms are introduced to assess the accuracy of the eigenvalue and eigenvector derivatives obtained by the reduction techniques. The effectiveness of the methods presented is demonstrated by three numerical examples.

Fetterman, Timothy L.↗

Dynamic analysis of pretwisted elastically-coupled rotor blades

This paper addresses the accuracy of using a one-dimensional analysis to predict frequencies of elastically-coupled highly-twisted rotor blades. Degrees of freedom associated with shear deformation are statically condensed from the formulation, so the analysis uses only those degrees of freedom associated with classical beam theory. The effects of cross section deformation (warping) are considered, and are shown to become significant for some types of elastic coupling. Improved results are demonstrated for highly-coupled blade structures through account of warping in a local cross section analysis, without explicit inclusion of these effects in the beam analysis. A convergence study is also provided which investigates the potential for improving efficiency of elastically-coupled beam analysis through implementation of a p-version beam finite element.

Nixon, Mark W.↗

Damage detection using experimentally measured mass and stiffness matrices

A method is presented for locating physical damage or change in a structure using experimentally measured mass and stiffness matrices. The approach uses a recently developed algorithm for transforming a state-space realization into a second order structural model with physical displacements as the generalized coordinates. This is accomplished by first rotating a state-space model of the identified structural dynamics into modal coordinates and approximating the mass normalized modal vectors for the output measurement set. Next, the physical mass, damping and stiffness matrices are synthesized directly from the measured modal parameters. This yields experimental mass and stiffness matrices for the structure without the use of a finite element model or a numerical search. The computed mass and stiffness are asymptotically equivalent to a static condensation of the global physical coordinate model. Techniques for solving the inverse connectivity problem are then developed whereby it is possible to assess the stiffness in a region of the structure bounded by several sensors. Applications to both simulated data and experimental data are used to discuss the effectiveness of the approach.

Peterson, L. D.↗

An alternative to Guyan reduction of finite-element models

Structural modeling is a key part of structural system identification for large space structures. Finite-element structural models are commonly used in practice because of their general applicability and availability. The initial models generated by using a standard computer program such as NASTRAN, ANSYS, SUPERB, STARDYNE, STRUDL, etc., generally contain tens of thousands of degrees of freedom. The models must be reduced for purposes of identification. Not only does the magnitude of the identification effort grow exponentially as a function of the number of degrees of freedom, but numerical procedures may also break down because of accumulated round-off errors. Guyan reduction is usually applied after a static condensation. Misapplication of Guyan reduction can lead to serious modeling errors. It is quite unfortunate and disappointing, since the accuracy of the original detailed finite-element model one tries very hard to achieve is lost by the reduction. First, why and how Guyan reduction always causes loss of accuracy is examined. An alternative approach is then introduced. The alternative can be thought of as an improvement of Guyan reduction, the Rayleigh-Ritz method, and in particular the recent algorithm of Wilson, Yuan, and Dickens. Unlike Guyan reduction, the use of the alternative does not need any special insight, experience, or skill for partitioning the structural degrees of freedom. In addition to model condensation, this alternative approach can also be used for predicting analytically, quickly, and economically, what are those structural modes that are excitable by a force actuator at a given trial location. That is, in the excitation of the structural modes for identification, it can be used for guiding the placement of the force actuators.

Lin, Jiguan Gene↗

Dynamic analysis of pretwisted elastically-coupled rotor blades

The accuracy of using a one-dimensional analysis to predict frequencies of elastically-coupled highly-twisted rotor blades is addressed. Degrees of freedom associated with shear deformation are statically condensed from the formulation, so the analysis uses only those degrees of freedom associated with classical beam theory. The effects of cross section deformation (warping) are considered, and are shown to become significant for some types of elastic coupling. Improved results are demonstrated for highly-coupled blade structures through account of warping in a local cross section analysis, without explicit inclusion of these effects in the beam analysis. A convergence study is also provided which investigates the potential for improving efficiency of elastically-coupled beam analysis through implementation of a p-version beam finite element.

Nixon, Mark W.↗

Static aeroelastic analysis of wings using Euler/Navier-Stokes equations coupled with improved wing-box finite element structures

Strong interactions between flow about an aircraft wing and the wing structure can result in aeroelastic phenomena which significantly impact aircraft performance. Time-accurate methods for solving the unsteady Navier-Stokes equations have matured to the point where reliable results can be obtained with reasonable computational costs for complex non-linear flows with shock waves, vortices and separations. The ability to combine such a flow solver with a general finite element structural model is key to an aeroelastic analysis in these flows. Earlier work involved time-accurate integration of modal structural models based on plate elements. A finite element model was developed to handle three-dimensional wing boxes, and incorporated into the flow solver without the need for modal analysis. Static condensation is performed on the structural model to reduce the structural degrees of freedom for the aeroelastic analysis. Direct incorporation of the finite element wing-box structural model with the flow solver requires finding adequate methods for transferring aerodynamic pressures to the structural grid and returning deflections to the aerodynamic grid. Several schemes were explored for handling the grid-to-grid transfer of information. The complex, built-up nature of the wing-box complicated this transfer. Aeroelastic calculations for a sample wing in transonic flow comparing various simple transfer schemes are presented and discussed.

Guruswamy, Guru P.↗

Advanced Modeling Strategies for the Analysis of Tile-Reinforced Composite Armor

A detailed investigation of the deformation mechanisms in tile-reinforced armored components was conducted to develop the most efficient modeling strategies for the structural analysis of large components of the Composite Armored Vehicle. The limitations of conventional finite elements with respect to the analysis of tile-reinforced structures were examined, and two complementary optimal modeling strategies were developed. These strategies are element layering and the use of a tile-adhesive superelement. Element layering is a technique that uses stacks of shear deformable shell elements to obtain the proper transverse shear distributions through the thickness of the laminate. The tile-adhesive superelement consists of a statically condensed substructure model designed to take advantage of periodicity in tile placement patterns to eliminate numerical redundancies in the analysis. Both approaches can be used simultaneously to create unusually efficient models that accurately predict the global response by incorporating the correct local deformation mechanisms.

Davila, Carlos G.↗

Discontinuous dual-primal mixed finite elements for elliptic problems

We propose a novel discontinuous mixed finite element formulation for the solution of second-order elliptic problems. Fully discontinuous piecewise polynomial finite element spaces are used for the trial and test functions. The discontinuous nature of the test functions at the element interfaces allows to introduce new boundary unknowns that, on the one hand enforce the weak continuity of the trial functions, and on the other avoid the need to define a priori algorithmic fluxes as in standard discontinuous Galerkin methods. Static condensation is performed at the element level, leading to a solution procedure based on the sole interface unknowns. The resulting family of discontinuous dual-primal mixed finite element methods is presented in the one and two-dimensional cases. In the one-dimensional case, we show the equivalence of the method with implicit Runge-Kutta schemes of the collocation type exhibiting optimal behavior. Numerical experiments in one and two dimensions demonstrate the order accuracy of the new method, confirming the results of the analysis.

Bottasso, Carlo L.↗

Methodology for Life Testing of Refractory Metal/Sodium Heat Pipes

The focus of this work was to establish an approach to generate carefully controlled data that can conclusively establish heat pipe operating life with material-fluid combinations capable of extended operation. To accomplish this goal acceleration is required to compress 10 years of operational life into 3 years of laboratory testing through a combination of increased temperature and mass fluence. Specific test series have been identi3ed, based on American Society for Testing and Materials (ASTM) specifications, to investigate long term corrosion rates. The refractory metal selected for demonstration purposes is a Molybdenum-44.5%Rhenium alloy formed by powder metallurgy. The heat pipe makes use of an annular crescent wick design formed by hot isostatic pressing of Molybdenum-Rhenium wire mesh. The heat pipes are filled using vacuum distillation and purity sampling is considered. Testing of these units is round-the-clock with 6-month destructive and non-destructive inspection intervals to identify the onset and level of corrosion. Non-contact techniques are employed for providing power to the evaporator (radio frequency induction heating at I to 5 kW per unit) and calorimetry at the condenser (static gas gap coupled water cooled calorimeter). The planned operating temperature range would extend from 1123 to 1323 K. Accomplishments prior to project cancellation included successful demonstration of the heat pipe wick fabrication technique, establishment of all engineering designs, baselined operational test requirements and procurement/assembly of supporting test hardware systems.

Martin, James J.↗

Accuracy of an approximate static structural analysis technique based on stiffness matrix eigenmodes

Use of the stiffness matrix eigenmodes, instead of the vibration eigenmodes, as generalized coordinates is proposed for condensation of static load deflection equations in finite element stiffness method. The modes are selected by strain energy criteria and the resulting fast, approximate analysis technique is evaluated by applications to idealized built-up wings and a fuselage segment. The best results obtained are a two-order of magnitude reduction of the number of degrees of freedom in a high aspect ratio wing associated with less than one percent error in prediction of the largest displacement.

Sobieszczanski-Sobieski, J.↗

GC/MS Method Development for Separating Lunar Volatile Ice Simulant Headspace Gases

Various investigators propose the lunar surface contains widely distributed volatiles, especially water- like species, i.e. OH and H2O. Surface volatiles are theorized to exist as a hydrated regolith layer, concentrated in extremely cold polar permanently shadowed regions (PSR), and/or solar wind implantation reservoirs in lunar glasses. The proposed sources of lunar surface volatiles range from cometary impacts, solar wind, or a supply present during moon formation. Future Artemis missions aim to collect and return the samples containing volatiles collected near lunar polar craters or PSRs. We, as advanced curation scientists, are responsible for developing techniques and methodologies for preserving returned sample integrity as much as possible. Pristine volatile-bearing samples are invaluable to the scientific community seeking to unravel the history of the solar system. Realistically, a sample will experience alteration during collection, transportation back to earth, and storage. The Planetary Exploration and Astromaterials Research Lab (PEARL) seeks to understand temperature and pressure effects on high-fidelity volatile-containing regolith simulants, the foundation for the future of cold curation. This abstract outlines the separation, identification, and quantification of headspace gases over volatile ice feed stock material using gas chromatography/mass spectrometry (GC/MS). Preliminary objectives concentrated on sample handling, reproducibility, and understanding the elution characteristics for each analyte. Initial GC/MS method development experiments utilized diluted static headspace sample preparation. Diluted samples were used because sampling headspace gases directly from a vial containing liquid analyte resulted in overloading of the column and detector. Overloading is evident based on chromatogram peak shapes and instrument contamination, or carry over, between experiments. A mixture of three alcohols were used for a majority of the sample handling and reproducibility studies. Reproducibility was tested via multiple users, calibration curves, and check standards. Stock solutions of condensed lunar volatile analytes included methanol, ammonia in methanol, hydrogen sulfide in water, and an equal volume mixture of methanol, ethanol, and isopropanol. Current samples use room air as the headspace sample matrix, however future experiments will incorporate an inert purge gas, such as argon or nitrogen. Three mL of each analyte solution were capped in separate 20 mL crimp top GC vials. Dilutions were carried out by removing an aliquot of headspace gases with a calibrated 1 mL gastight syringe and immediately transferring to a 20 mL capped crimp top vial. The GC/MS is a Thermo Fisher Trace 1310/ISQ 7000 with a TriPlus RSH autosampler and split/splitless injector module. The experiments outlined in this abstract use the following hardware: a 2.5 mL gastight headspace syringe tool, 1 mm ID x 78.5 mm length ultra-inert straight injection liner, and a TG-BondQ 30 m × 0.32 mm × 10 μm column. Various parameters, such as hardware selection and the temperature, pressure, and split ratio set points, continue to evolve as the overall experiment is refined. Diluted headspace chromatograms were collected for the individual stock solutions. Retention times, peak shapes, and mass spectra were evaluated and added to the data processing method for each molecule of interest. Figure 1 shows the total ion chromatograms for the three major lunar volatile simulant stock solutions: methanol, 7 N ammonia in methanol, and 0.4% hydrogen sulfide in water. Tailing peak shapes for ammonia (2.98 min rt) and water (4.06 min rt) indicate the molecules are not properly eluting from the selected column with the current separation method. Additionally, hydrogen sulfide and ammonia have overlapping peak windows, which could impact quantification. Ongoing experiments aim to address the peak shape and overlapping via the separation method and hardware selection. Sample preparation reproducibility experiments used stock solution containing equal volumes of a non- interactive mixture of methanol, ethanol, and isopropanol. Mass spectrum ion traces were used to identify and quantify all three alcohols. Peaks were automatically detected, identified, and integrated through the mass spectra detection and processing parameters. Calibration response curves and check standards were used to evaluate the validity of the sample preparation procedure. Figure 2 shows the methanol chromatogram peak area versus total headspace dilution volume transferred from the alcohol mixture vial. The calibration response curves and check standards validate sample preparation procedure. Continuing data analysis efforts are working towards correlating the peak area and instrument response factor to the headspace analyte concentration and condensed phase composition. Static headspace gas chromatography theory relies on Dalton’s law, Raoult’s law, Henry’s Law, and the Kolb and Ettre equation to associate peak area to the analyte composition in a non-ideal solution. Equation 1 is a simplified expression derived from the aforementioned theories. Future experiments involve liquid injections of the individual stock solutions, liquid and headspace analysis of various stock solution combinations, and the addition of regolith simulants to the mixtures. Temperature is another variable expected to affect reaction rates and will be explored.

Cecilia L. Amick↗