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 415 records · Page 23

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↗

Model reduction in a subset of the original states

A model reduction method is investigated to provide a smaller structural dynamic model for subsequent structural control design. A structural dynamic model is assumed to be derived from finite element analysis. It is first converted into the state space form, and is further reduced by the internal balancing method. Through the co-ordinate transformation derived from the states that are deleted during reduction, the reduced model is finally expressed with the states that are members of the original states. Therefore, the states in the final reduced model represent the degrees of freedom of the nodes that are selected by the designer. The procedure provides a more practical implementation of model reduction for applications in which specific nodes, such as sensor and/or actuator attachment points, are to be retained in the reduced model. Thus, it ensures that the reduced model is under the same input and output condition as the original physical model. The procedure is applied to two simple examples and comparisons are made between the full and reduced order models. The method can be applied to a linear, continuous and time-invariant model of structural dynamics with nonproportional viscous damping.

Yae, K. H.↗

Numerical simulation using vorticity-vector potential formulation

An accurate and efficient computational method is needed for three-dimensional incompressible viscous flows in engineering applications. On solving the turbulent shear flows directly or using the subgrid scale model, it is indispensable to resolve the small scale fluid motions as well as the large scale motions. From this point of view, the pseudo-spectral method is used so far as the computational method. However, the finite difference or the finite element methods are widely applied for computing the flow with practical importance since these methods are easily applied to the flows with complex geometric configurations. However, there exist several problems in applying the finite difference method to direct and large eddy simulations. Accuracy is one of most important problems. This point was already addressed by the present author on the direct simulations on the instability of the plane Poiseuille flow and also on the transition to turbulence. In order to obtain high efficiency, the multi-grid Poisson solver is combined with the higher-order, accurate finite difference method. The formulation method is also one of the most important problems in applying the finite difference method to the incompressible turbulent flows. The three-dimensional Navier-Stokes equations have been solved so far in the primitive variables formulation. One of the major difficulties of this method is the rigorous satisfaction of the equation of continuity. In general, the staggered grid is used for the satisfaction of the solenoidal condition for the velocity field at the wall boundary. However, the velocity field satisfies the equation of continuity automatically in the vorticity-vector potential formulation. From this point of view, the vorticity-vector potential method was extended to the generalized coordinate system. In the present article, we adopt the vorticity-vector potential formulation, the generalized coordinate system, and the 4th-order accurate difference method as the computational method. We present the computational method and apply the present method to computations of flows in a square cavity at large Reynolds number in order to investigate its effectiveness.

Tokunaga, Hiroshi↗

Further generalization of an equivalent plate representation for aircraft structural analysis

Recent developments from a continuing effort to provide an equivalent plate representation for aircraft structural analysis are described. Previous work provided an equivalent plate analysis formulation that is capable of modeling aircraft wing structures with a general planform such as cranked wing boxes. However, the modeling is restricted to representing wing boxes having symmetric cross sections. Further developments, which are described, allow modeling of wing cross sections having asymmetries that can arise from airfoil camber or from thicknesses being different in the upper and lower cover skins. An implementation of thermal loadings, which are described as temperature distributions over the planform of the cover skins, has been included. Spring supports have been added to provide for a more general set of boundary conditions. Numerical results are presented to assess the effect of wing camber on the static and dynamic response of an example wing structure under pressure and thermal loading. These results are compared with results from a finite element analysis program to indicate how well a cambered wing box can be represented with an equivalent plate formulation.

Giles, Gary L.↗

Further generalization of an equivalent plate representation for aircraft structural analysis

Recent developments from a continuing effort to provide an equivalent plate representation for aircraft structural analysis are described. Previous work provided an equivalent plate analysis formulation that is capable of modeling aircraft wing structures with a general planform such as cranked wing boxes. However, the modeling is restricted to representing wing boxes having symmetric cross sections. Further developments, which are described in this paper, allow modeling of wing cross section having asymmetries that can arise from airfoil camber or from thicknesses being different in the upper and lower cover skins. An implementation of thermal loadings, which are described as temperature distributions over the planform of the cover skins, has been included. Spring supports have been added to provide for a more general set of boundary conditions. Numerical results are presented to assess the effect of wing camber on the static and dynamic response of an example wing structure under pressure and thermal loading. These results are compared with results from a finite element analysis program to indicate how well a cambered wing box can be represented with an equivalent plate formulation.

Giles, Gary L.↗

Development of 3D electromagnetic modeling tools for airborne vehicles

The main goal of this project is to develop methodologies for scattering by airborne composite vehicles. Although our primary focus continues to be the development of a general purpose code for analyzing the entire structure as a single unit, a number of other tasks are also pursued in parallel with this effort. These tasks are important in testing the overall approach and in developing suitable models for materials coatings, junctions and, more generally, in assessing the effectiveness of the various parts comprising the final code. Here, we briefly discuss our progress on the five different tasks which were pursued during this period. Our progress on each of these tasks is described in the detailed reports (listed at the end of this report) and the memoranda included. The first task described below is, of course, the core of this project and deals with the development of the overall code. Undoubtedly, it is the outcome of the research which was funded by NASA-Ames and the Navy over the past three years. During this year we developed the first finite element code for scattering by structures of arbitrary shape and composition. The code employs a new absorbing boundary condition which allows termination of the finite element mesh only 0.3 lambda from the outer surface of the target. This leads to a remarkable reduction of the mesh size and is a unique feature of the code. Other unique features of this code include capabilities to model resistive sheets, impedance sheets and anisotropic materials. This last capability is the latest feature of the code and is still under development. The code has been extensively validated for a number of composite geometries and some examples are given. The validation of the code is still in progress for anisotropic and larger non-metallic geometries and cavities. The developed finite element code is based on a Galerkin's formulation and employs edge-based tetrahedral elements for discretizing the dielectric sections and the region between the target and the outer mesh termination boundary (ATB). This boundary is placed in conformity with the target's outer surface, thus resulting in additional reduction of the unknown count.

Volakis, John L.↗