Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuous finite element”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32

How Certain Physical Considerations Impact Aerostructural Wing Optimization

Wing design optimization has been studied extensively and is of continued interest as optimization tools are developed and become more accessible. In each of these studies, certain assumptions and simplifications are made to make the design problem tractable. However, it is difficult to find systematic studies in which several considerations are added or removed one at a time to study how much impact they have. In this work, we examine how certain physical considerations (viscous drag, wave drag, thrust loads, and inertial relief from structural, fuel, and engine masses), impact the aerostructural optimization results for three distinct aircraft wings. The goal is to help develop a rough idea of how important these physical considerations are. We do this using gradient-based optimization and a multidisciplinary design optimization framework, OpenMDAO. We use the open-source tool OpenAeroStruct that couples a vortex lattice method to a finite element method. We establish a baseline aerostructural design optimization problem then perform a series of optimizations, each with one physical consideration removed from the baseline case. We find that depending on the size of the aircraft and flight conditions, the importance of some of these physical considerations varies considerably whereas the importance of others do not. Specifically, the optimal designs change radically without proper viscous and wave drag considerations and smaller aircraft with more distributed propulsion are more affected by the inclusion of engine loads.

Aeropropulsion↗

Coherent manipulation of graph states composed of finite-energy Gottesman-Kitaev-Preskill-encoded qubits

Graph states are a central resource in measurement-based quantum information processing. In the photonic qubit architecture based on Gottesman-Kitaev-Preskill (GKP) encoding, the generation of high-fidelity graph states composed of realistic, finite-energy approximate GKP-encoded qubits thus constitutes a key task. We consider the finite-energy approximation of GKP-qubit states given by a coherent superposition of shifted finite-squeezed vacuum states, where the displacements are Gaussian distributed. We present an exact description of graph states composed of such approximate GKP qubits as a coherent superposition of a Gaussian ensemble of randomly displaced ideal GKP-qubit graph states. Using standard Gaussian dynamics, we track the transformation of the covariance matrix and the mean-displacement vector elements of the Gaussian distribution of the ensemble under tools such as GKP-Steane error-correction and fusion operations that can be used to grow large high-fidelity GKP-qubit graph states. The covariance matrix elements capture the noise in the graph state due to the finite-energy approximation of GKP qubits, while the mean displacements relate to the possible absolute shift errors on the individual qubits arising conditionally from the homodyne measurements that are a part of these tools. Our work thus pins down an exact coherent error model for graph states generated from truly finite-energy GKP qubits, which can shed light on their error-correction properties.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Structural Sizing of a 50-m-Tall Thermoplastic Composite Solar Array Truss Tower Structure for the Lunar South Pole

At the lunar south pole, solar arrays elevated by truss towers are exposed to near continuous sunlight since the sun remains near the horizon. Recent studies suggest that tower heights on the order of 50 m are needed to meet the power requirements to support a sustained lunar presence. While eventually towers may be constructed using materials sourced in-situ, initially structural elements will likely be brought from the Earth, deployed, and assembled. Thermoplastic composites offer the benefit of high specific stiffness and strength properties along with welding for assembly joints, and therefore represent a promising material system for this application. The focus of this paper is on the structural sizing of such towers. A set of expressions based on beam theory are presented for preliminary sizing of the truss tower structure for strength and buckling. The loading condition that drives sizing is base excitation resulting from moonquakes. A point design, developed through finite element analysis, verifies the analytical sizing routine and shows the importance of joints. The results show a viable point design and highlight the factors most significant in the structural sizing.

Thermoplastic composites↗

Structural Sizing of a 50-m-Tall Thermoplastic Composite Solar Array Truss Tower Structure for the Lunar South Pole

At the lunar south pole, solar arrays elevated by truss towers are exposed to near continuous sunlight since the sun remains near the horizon. Recent studies suggest that tower heights on the order of 50 m are needed to meet the power requirements to support a sustained lunar presence. While eventually towers may be constructed using materials sourced in-situ, initially structural elements will likely be brought from the Earth, deployed, and assembled. Thermoplastic composites offer the benefit of high specific stiffness and strength properties along with welding for assembly joints, and therefore represent a promising material system for this application. The focus of this paper is on the structural sizing of such towers. A set of expressions based on beam theory are presented for preliminary sizing of the truss tower structure for strength and buckling. The loading condition that drives sizing is base excitation resulting from moonquakes. A point design, developed through finite element analysis, verifies the analytical sizing routine and shows the importance of joints. The results show a viable point design and highlight the factors most significant in the structural sizing.

Thermoplastic composites↗

Structural Sizing of a 50-m-Tall Thermoplastic Composite Truss Tower Structure for the Lunar South Pole

At the lunar south pole, solar arrays elevated by truss towers are exposed to near continuous sunlight since the sun remains near the horizon. Recent studies suggest that tower heights on the order of 50 m are needed to meet the power requirements to support a sustained lunar presence. While eventually towers may be constructed using materials sourced in-situ, initially structural elements will likely be brought from the Earth, deployed, and assembled. Thermoplastic composites offer the benefit of high specific stiffness and strength properties along with welding for assembly joints, and therefore represent a promising material system for this application. The focus of this paper is on the structural sizing of such towers. A set of expressions based on beam theory are presented for preliminary sizing of the truss tower structure for strength and buckling. The loading condition that drives sizing is base excitation resulting from moonquakes. A point design, developed through finite element analysis, verifies the analytical sizing routine and shows the importance of joints. The results show a viable point design and highlight the factors most significant in the structural sizing.

Thermoplastic composites↗

Experimental Verification of Computational Models for Laminated Composites

The objective of the research reported herein is to develop a progressive damage methodology capable of predicting the residual strength of continuous fiber-reinforced, laminated, polymer matrix composites with through-penetration damage. The fracture behavior of center-notch tension panels with thin crack-like slits was studied. Since fibers are the major load-carrying constituent in polymer matrix composites, predicting the residual strength of a laminate requires a criterion for fiber fracture. The effects on fiber strain due to other damage mechanisms such as matrix cracking and delaminations must also be modeled. Therefore, the research herein examines the damage mechanisms involved in translaminate fracture and identifies the toughening mechanisms responsible for damage growth resistance in brittle epoxy matrix systems. The mechanics of matrix cracking and fiber fracture are discussed as is the mathematical framework for the progressive damage model developed by the authors. The progressive damage analysis algorithms have been implemented into a general purpose finite element code developed by NASA, the Computational Structural Mechanics Testbed (COMET). Damage growth is numerically simulated and the analytical residual strength predictions are compared to experimental results for a variety of notched panel configurations and materials systems.

Harris, Charles E.↗

Structural Health Management of Damaged Aircraft Structures Using the Digital Twin Concept

The development of multidisciplinary integrated Structural Health Management (SHM) tools will enable accurate detection, and prognosis of damaged aircraft under normal and adverse conditions during flight. As part of the digital twin concept, methodologies are developed by using integrated multiphysics models, sensor information and input data from an in-service vehicle to mirror and predict the life of its corresponding physical twin. SHM tools are necessary for both damage diagnostics and prognostics for continued safe operation of damaged aircraft structures. The adverse conditions include loss of control caused by environmental factors, actuator and sensor faults or failures, and structural damage conditions. A major concern in these structures is the growth of undetected damage/cracks due to fatigue and low velocity foreign object impact that can reach a critical size during flight, resulting in loss of control of the aircraft. To avoid unstable, catastrophic propagation of damage during a flight, load levels must be maintained that are below a reduced load-carrying capacity for continued safe operation of an aircraft. Hence, a capability is needed for accurate real-time predictions of damage size and safe load carrying capacity for structures with complex damage configurations. In the present work, a procedure is developed that uses guided wave responses to interrogate damage. As the guided wave interacts with damage, the signal attenuates in some directions and reflects in others. This results in a difference in signal magnitude as well as phase shifts between signal responses for damaged and undamaged structures. Accurate estimation of damage size, location, and orientation is made by evaluating the cumulative signal responses at various pre-selected sensor locations using a genetic algorithm (GA) based optimization procedure. The damage size, location, and orientation is obtained by minimizing the difference between the reference responses and the responses obtained by wave propagation finite element analysis of different representative cracks, geometries, and sizes.

Seshadri, Banavara R.↗

Infinite possibilities: Computational structures technology

Computational Fluid Dynamics (or CFD) methods are very familiar to the research community. Even the general public has had some exposure to CFD images, primarily through the news media. However, very little attention has been paid to CST--Computational Structures Technology. Yet, no important design can be completed without it. During the first half of this century, researchers only dreamed of designing and building structures on a computer. Today their dreams have become practical realities as computational methods are used in all phases of design, fabrication and testing of engineering systems. Increasingly complex structures can now be built in even shorter periods of time. Over the past four decades, computer technology has been developing, and early finite element methods have grown from small in-house programs to numerous commercial software programs. When coupled with advanced computing systems, they help engineers make dramatic leaps in designing and testing concepts. The goals of CST include: predicting how a structure will behave under actual operating conditions; designing and complementing other experiments conducted on a structure; investigating microstructural damage or chaotic, unpredictable behavior; helping material developers in improving material systems; and being a useful tool in design systems optimization and sensitivity techniques. Applying CST to a structure problem requires five steps: (1) observe the specific problem; (2) develop a computational model for numerical simulation; (3) develop and assemble software and hardware for running the codes; (4) post-process and interpret the results; and (5) use the model to analyze and design the actual structure. Researchers in both industry and academia continue to make significant contributions to advance this technology with improvements in software, collaborative computing environments and supercomputing systems. As these environments and systems evolve, computational structures technology will evolve. By using CST in the design and operation of future structures systems, engineers will have a better understanding of how a system responds and lasts, more cost-effective methods of designing and testing models, and improved productivity. For informational and educational purposes, a videotape is being produced using both static and dynamic images from research institutions, software and hardware companies, private individuals, and historical photographs and drawings. The extensive number of CST resources indicates its widespread use. Applications run the gamut from simpler university-simulated problems to those requiring solutions on supercomputers. In some cases, an image or an animation will be mapped onto the actual structure to show the relevance of the computer model to the structure. Transferring the digital files to videotape presents a number of problems related to maintaining the quality of the original image, while still producing a broadcast quality videotape. Since researchers normally do not create a computer image using traditional composition theories or video production requirements, often the image loses some of its original digital quality and impact when transferred to videotape. Although many CST images are currently available, those that are edited into the final project must meet two important criteria: they must complement the narration, and they must be broadcast quality when recorded on videotape.

Beam, Sherilee F.↗

Fracture Mechanics Analysis of Stitched Stiffener-Skin Debonding

An analysis based on plate finite elements and the virtual crack closure technique has been implemented to study the effect of stitching on mode I and mode II strain energy release rates for debond configurations. The stitches were modeled as discrete nonlinear fastener elements with a compliance determined by experiment. The axial and shear behavior of the stitches was considered, however, the two compliances and failure loads were assumed to be independent. Both a double cantilever beam (mode I) and a mixed mode skin-stiffener debond configuration were studied. In the double cantilever beam configurations, G(sub I) began to decrease once the debond had grown beyond the first row of stitches and was reduced to zero for long debonds. In the mixed-mode skin-stiffener configurations, G(sub I) showed a similar behavior as in the double cantilever beam configurations, however, G(sub u), continued to increase with increasing debond length.

Glaessgen, E. H.↗

Remote Sensing of Surficial Process Responses to Extreme Meteorological Events

Changes in the frequency and magnitude of extreme meteorological events are associated with changing environmental means. Such events are important in human affairs, and can also be investigated by orbital remote sensing. During the course of this project, we applied ERS-1, ERS-2, Radarsat, and an airborne sensor (AIRSAR-TOPSAR) to measure flood extents, flood water surface profiles, and flood depths. We established a World Wide Web site (the Dartmouth Flood Observatory) for publishing remote sensing-based maps of contemporary floods worldwide; this is also an online "active archive" that presently constitutes the only global compilation of extreme flood events. We prepared an article for EOS concerning SAR imaging of the Mississippi Valley flood; an article for the International Journal of Remote Sensing on measurement of a river flood wave using ERS-2, began work on an article (since completed and published) on the Flood Observatory for a Geoscience Information Society Proceedings volume, and presented lectures at several Geol. Soc. of America Natl. Meetings, an Assoc. of Amer. Geographers Natl. Meeting, and a Binghamton Geomorphology Symposium (all on SAR remote sensing of the Mississippi Valley flood). We expanded in-house modeling capabilities by installing the latest version of the Army Corps of Engineers RMA two-dimensional hydraulics software and BYU Engineering Graphics Lab's Surface Water Modeling System (finite elements based pre- and post-processors for RMA work) and also added watershed modeling software. We are presently comparing the results of the 2-d flow models with SAR image data. The grant also supported several important upgrades of pc-based remote sensing infrastructure at Dartmouth. During work on this grant, we collaborated with several workers at the U.S. Army Corps of Engineers, Remote Sensing/GIS laboratory (for flood inundation mapping and modeling; particularly of the Illinois River using the AIRSAR/TOPSAR/ERS-2 combined data), with Dr. Karen Prestegaard at the University of Maryland (geomorphological responses to the extreme 1993 flood along the Raccoon drainage in central Iowa), and with Mr Tim Scrom of the Albany National Weather Service River Forecast Center (initial planning for the use of Radarsat and ERS-2 for flood warning). The work thus initiated with this proposal is continuing.

Brakenridge, G. Robert↗

Computational procedures for postbuckling of composite shells

A recently developed finite-element capability for general nonlinear shell analysis, featuring the use of three-dimensional constitutive equations within an efficient resultant-oriented framework, is employed to simulate the postbuckling response of an axially compressed composite cylindrical panel with a circular cutout. The problem is a generic example of modern composite aircraft components for which postbuckling strength (i.e., fail-safety) is desired in the presence of local discontinuities such as holes and cracked stiffeners. While the computational software does a reasonable job of predicting both the buckling load and the qualitative aspects of postbuckling (compared both with experiment and another code) there are some discrepancies due to: (1) uncertainties in the nominal layer material properties, (2) structural sensitivity to initial imperfections, and (3) the neglect of dynamic and local material delamination effects in the numerical model. Corresponding refinements are suggested for the realistic continuation of this type of analysis.

Stanley, G. M.↗

A time dependent difference theory for sound propagation in ducts with flow

A time dependent numerical solution of the linearized continuity and momentum equation is developed for sound propagation in a two-dimensional straight hard or soft wall duct with a sheared mean flow. The time dependent governing acoustic-difference equations and boundary conditions are developed along with a numerical determination of the maximum stable time increments. The analysis begins with a harmonic noise source radiating into a quiescent duct. This explicit iteration method then calculates stepwise in real time to obtain the transient as well as the 'steady' state solution of the acoustic field. Example calculations are presented for sound propagation in hard and soft wall ducts, with no flow and with plug flow. Although the problem with sheared flow has been formulated and programmed, sample calculations have not yet been examined. So far, the time dependent finite difference analysis has been found to be superior to the steady state finite difference and finite element techniques because of shorter solution times and the elimination of large matrix storage requirements.

Baumeister, K. J.↗

ZIP2DL: An Elastic-Plastic, Large-Rotation Finite-Element Stress Analysis and Crack-Growth Simulation Program

ZIP2DL is a two-dimensional, elastic-plastic finte element program for stress analysis and crack growth simulations, developed for the NASA Langley Research Center. It has many of the salient features of the ZIP2D program. For example, ZIP2DL contains five material models (linearly elastic, elastic-perfectly plastic, power-law hardening, linear hardening, and multi-linear hardening models), and it can simulate mixed-mode crack growth for prescribed crack growth paths under plane stress, plane strain and mixed state of stress conditions. Further, as an extension of ZIP2D, it also includes a number of new capabilities. The large-deformation kinematics in ZIP2DL will allow it to handle elastic problems with large strains and large rotations, and elastic-plastic problems with small strains and large rotations. Loading conditions in terms of surface traction, concentrated load, and nodal displacement can be applied with a default linear time dependence or they can be programmed according to a user-defined time dependence through a user subroutine. The restart capability of ZIP2DL will make it possible to stop the execution of the program at any time, analyze the results and/or modify execution options and resume and continue the execution of the program. This report includes three sectons: a theoretical manual section, a user manual section, and an example manual secton. In the theoretical secton, the mathematics behind the various aspects of the program are concisely outlined. In the user manual section, a line-by-line explanation of the input data is given. In the example manual secton, three types of examples are presented to demonstrate the accuracy and illustrate the use of this program.

Deng, Xiaomin↗

Finite-Strain Fractional-Order Viscoelastic (FOV) Material Models and Numerical Methods for Solving Them

Fraction-order viscoelastic (FOV) material models have been proposed and studied in 1D since the 1930's, and were extended into three dimensions in the 1970's under the assumption of infinitesimal straining. It was not until 1997 that Drozdov introduced the first finite-strain FOV constitutive equations. In our presentation, we shall continue in this tradition by extending the standard, FOV, fluid and solid, material models introduced in 1971 by Caputo and Mainardi into 3D constitutive formula applicable for finite-strain analyses. To achieve this, we generalize both the convected and co-rotational derivatives of tensor fields to fractional order. This is accomplished by defining them first as body tensor fields and then mapping them into space as objective Cartesian tensor fields. Constitutive equations are constructed using both variants for fractional rate, and their responses are contrasted in simple shear. After five years of research and development, we now possess a basic suite of numerical tools necessary to study finite-strain FOV constitutive equations and their iterative refinement into a mature collection of material models. Numerical methods still need to be developed for efficiently solving fraction al-order integrals, derivatives, and differential equations in a finite element setting where such constitutive formulae would need to be solved at each Gauss point in each element of a finite model, which can number into the millions in today's analysis.

Freed, Alan D.↗

Design and Analysis of Tow-Steered Composite Shells Using Fiber Placement

In this study, a sub-scale advanced composite shell design is evaluated to determine its potential for use on a future aircraft fuselage. Two composite shells with the same nominal 8-ply [+/-45/+/-Theta](sub s) layup are evaluated, where Theta indicates a tow-steered ply. To build this shell, a fiber placement machine would be used to steer unidirectional prepreg tows as they are placed around the circumference of a 17-inch diameter right circular cylinder. The fiber orientation angle varies continuously from 10 degrees (with respect to the shell axis of revolution) at the crown, to 45 degrees on the side, and back to 10 degrees on the keel. All 24 tows are placed at each point on every fiber path in one structure designated as the shell with overlaps. The resulting pattern of tow overlaps causes the laminate thickness to vary between 8 and 16 plies. The second shell without tow overlaps uses the capability of the fiber placement machine to cut and add tows at any point along the fiber paths to fabricate a shell with a nearly uniform 8-ply laminate thickness. Issues encountered during the design and analysis of these shells are presented and discussed. Static stiffness and buckling loads of shells with tow-steered layups are compared with the performance of a baseline quasi-isotropic shell using both finite element analyses and classical strength of materials theory.

Wu, K. Chauncey↗

Trans-Laminar-Reinforced (TLR) Composites

A Trans-Laminar-Reinforced (TLR) composite is defined as composite laminate with up to five percent volume of fibrous reinforcement oriented in a 'trans-laminar' fashion in the through-thickness direction. The TLR can be continuous threads as in 'stitched laminates', or it can be discontinuous rods or pins as in 'Z-Fiber(TM) materials. It has been repeatedly documented in the literature that adding TLR to an otherwise two dimensional laminate results in the following advantages: substantially improved compression-after-impact response; considerably increased fracture toughness in mode 1 (double cantilever beam) and mode 2 (end notch flexure); and severely restricted size and growth of impact damage and edge delamination. TLR has also been used to eliminate catastrophic stiffener disbonding in stiffened structures. TLR directly supports the 'Achilles heel' of laminated composites, that is delamination. As little as one percent volume of TLR significantly alters the mechanical response of laminates. The objective of this work was to characterize the effects of TLR on the in-plane and inter-laminar mechanical response of undamaged composite laminates. Detailed finite element models of 'unit cells', or representative volumes, were used to study the effects of adding TLR on the elastic constants; the in-plane strength; and the initiation of delamination. Parameters investigated included TLR material, TLR volume fraction, TLR diameter, TLR through-thickness angle, ply stacking sequence, and the microstructural features of pure resin regions and curved in-plane fibers. The work was limited to the linear response of undamaged material with at least one ply interface. An inter-laminar dominated problem of practical interest, a flanged skin in bending, was also modeled.

Hinders, Mark↗

Global Dynamic Modeling of Space-Geodetic Data

The proposal had outlined a year for program conversion, a year for testing and debugging, and two years for numerical experiments. We kept to that schedule. In first (partial) year, author designed a finite element for isostatic thin-shell deformation on a sphere, derived all of its algebraic and stiffness properties, and embedded it in a new finite element code which derives its basic solution strategy (and some critical subroutines) from earlier flat-Earth codes. Also designed and programmed a new fault element to represent faults along plate boundaries. Wrote a preliminary version of a spherical graphics program for the display of output. Tested this new code for accuracy on individual model plates. Made estimates of the computer-time/cost efficiency of the code for whole-earth grids, which were reasonable. Finally, converted an interactive graphical grid-designer program from Cartesian to spherical geometry to permit the beginning of serious modeling. For reasons of cost efficiency, models are isostatic, and do not consider the local effects of unsupported loads or bending stresses. The requirements are: (1) ability to represent rigid rotation on a sphere; (2) ability to represent a spatially uniform strain-rate tensor in the limit of small elements; and (3) continuity of velocity across all element boundaries. Author designed a 3-node triangle shell element which has two different sets of basis functions to represent (vector) velocity and all other (scalar) variables. Such elements can be shown to converge to the formulas for plane triangles in the limit of small size, but can also applied to cover any area smaller than a hemisphere. The difficult volume integrals involved in computing the stiffness of such elements are performed numerically using 7 Gauss integration points on the surface of the sphere, beneath each of which a vertical integral is performed using about 100 points.

Bird, Peter↗

Study of Microphonic Effects on the C100 Cryomodule for High Energy Electron Beam Accelerators

The Continuous Electron Beam Accelerator Facility (CEBAF) at Thomas Jefferson National Laboratory (JLab) is a particle accelerator which can accelerate an electron beam to relativistic speeds and apply the beam onto target samples. The C100 superconducting radio frequency (SRF) cavity is the primary accelerating structure of the C100 cryomodule, one of the many cryomodules which compose the CEBAF linear accelerator. SRF cavities are particularly sensitive to internal and external vibrations that can result in a phenomenon called microphonics which degrade the operational stability of a cryomodule. The purpose of this thesis is to investigate the significance of mechanical disturbances on the electromagnetic resonant frequency of a C100 SRF cavity. Knowledge of the mechanical resonance of the cavities and cryomodule sheds light into how these disturbances are most easily realized as deformation which causes radio frequency (RF) detuning. Three studies were conducted: the development and hammer test calibration of a Finite Element Analysis (FEA) model of a C100 cavity, the development and hammer test calibration of an FEA model of a C100 cavity string, and the hammer test of the C100-10R cryomodule at the Cryomodule Test Facility (CMTF). The cavity FEA model was found to accurately predict two modes found in two real cavities in a simply supported configuration. The cavity string FEA model leveraged the calibrated cavity FEA model but was not found to accurately predict the modal behavior of a real cavity string. Even so, the modal behavior of the cavity string inside the C100-10R cryomodule was captured during a hammer test while it was partially assembled. Finally, the C100-10R cryomodule was placed in the CMTF to study RF detuning. The RF detuning spectra during hammer hits and background noise was captured. The results of the hammer testing indicate two strong peaks at low frequencies (9-10 Hz and 22-23 Hz). These two frequencies were found to be nearly coincident to four instances of mechanical resonance found during the hammer testing done on the partially-assembled C100-10R. Because of this, these two modes are believed to contribute to RF detuning of the cryomodule. This test event also included the testing of the effectiveness of a configuration of BNNT canisters designed to act as dampers. While these tests show promising results, the lurking variables render these tests somewhat inconclusive.

Hull, Caleb James↗