Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shape function”

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 73 records · Page 4

Wave drag as the objective function in transonic fighter wing optimization

The original computational method for determining wave drag in a three dimensional transonic analysis method was replaced by a wave drag formula based on the loss in momentum across an isentropic shock. This formula was used as the objective function in a numerical optimization procedure to reduce the wave drag of a fighter wing at transonic maneuver conditions. The optimization procedure minimized wave drag through modifications to the wing section contours defined by a wing profile shape function. A significant reduction in wave drag was achieved while maintaining a high lift coefficient. Comparisons of the pressure distributions for the initial and optimized wing geometries showed significant reductions in the leading-edge peaks and shock strength across the span.

Phillips, P. S.↗

Curvilinear Displacement Transfer Functions for Deformed Shape Predictions of Curved Structures Using Distributed Surface Strains

Curvilinear Displacement Transfer Functions were formulated for deformed shape predictions of different curved structures using surface strains. The embedded curved beam (depth-wise cross section of a curved structure along a surface strain-sensing line) was discretized into multiple small domains, with domain junctures matching the strain-sensing stations. Thus, the surface strain distribution can be described with a piecewise linear or a piecewise nonlinear function. The discrete approach enabled piecewise integrations of a curvature-strain differential equation for the embedded curved beam to yield closed-form Curvilinear Displacement Transfer Functions, which are written in terms of embedded curved-beam geometrical parameters and surface strains. By inputting the surface strain data, the Curvilinear Displacement Transfer Functions can transform surface strains into deflections along each embedded curved beam for mapping out the overall structural deformed shapes. The finite-element method was used to analytically generate the surface strains of the curved beams. The deformed shape prediction accuracies were then determined by comparing the theoretical deflections with the finite-element-generated deflections, which were used as yardsticks. By introducing the correction factors in simple mathematical forms, the Curvilinear Displacement Transfer Functions can be quite accurate for shape predictions of different curved-beam structures ranging from limit case of straight beam up to semicircular curved beam.

Ko, William L.↗

Thermal analysis of a functionally graded material subject to a thermal gradient using the boundary element method

The boundary element method is utilized in this study to conduct thermal analysis of functionally graded composites, materials in which the internal microstructure or properties are explicitly tailored in order to obtain an optimal response, on the micromechanical (constituent) scale. A unique feature of the boundary element formulations used here is the use of circular shape functions to convert the two-dimensional integrations of the composite fibers to one dimensional integrations. Using the computer code BEST-CMS, the through the thickness temperature profiles are computed for a representative material with varying numbers of fibers and fiber spacing in the thickness direction. The computed temperature profiles are compared to those obtained using an alternate analytical theory which explicitly couples the heterogeneous microstructure to the global analysis. The boundary element results compared favorably to the analytical calculations, with discrepancies that are explainable based on the boundary element formulation. The results serve both to demonstrate the ability of the boundary element method to analyze these types of materials, and to verify the accuracy of the analytical theory.

Goldberg, Robert K.↗

A high-precision finite element method for shock-tube calculations

A two-pass explicit scheme is developed in order to exploit some of the capabilities of finite difference modeling (FDM) for finite element modeling (FEM), which offers the opportunity to account for any type of geometry in fluid flow modeling. Features of the first-order upwind and the Lax-Wendroff high precision explicit finite difference algorithms are reviewed. A flux limiter is developed for FEM to serve as an analog for the single limiter function which has been defined for the various FDMs. It is shown that an antidiffusive limiter must be introduced into the weighting function which normally multiplies the time-derivative term in the variational equation. The two-pass scheme which results is demonstrated to be the equivalent of FDMs with five-point support. However, the present scheme is valid only for one-dimensional calculations and linear shape functions for shock tube flow phenomena. Further work is required for its use with nonlinear hyperbolic systems.

Hughes, T. J. R.↗

Application of numerical optimization to the design of supercritical airfoils without drag-creep

Recent applications of numerical optimization to the design of advanced airfoils for transonic aircraft have shown that low-drag sections can be developed for a given design Mach number without an accompanying drag increase at lower Mach numbers. This is achieved by imposing a constraint on the drag coefficient at an off-design Mach number while the drag at the design Mach number is the objective function. Such a procedure doubles the computation time over that for single design-point problems, but the final result is worth the increased cost of computation. The ability to treat such multiple design-point problems by numerical optimization has been enhanced by the development of improved airfoil shape functions. Such functions permit a considerable increase in the range of profiles attainable during the optimization process.

Hicks, R. M.↗

Three-dimensional canard-wing shape optimization in aircraft cruise and maneuver environments

This paper demonstrates a numerical technique for canard-wing shape optimization at two operating conditions. For purposes of simplicity, a mean surface wing paneling code is employed for the aerodynamic calculations. The optimization procedures are based on the method of feasible directions. The shape functions for describing the thickness, camber, and twist are based on polynomial representations. The primary design requirements imposed restrictions on the canard and wing volumes and on the lift coefficients at the operating conditions. Results indicate that significant improvements in minimum drag and lift-to-drag ratio are possible with reasonable aircraft geometries. Calculations were done for supersonic speeds with Mach numbers ranging from 1 to 6. Planforms were mainly of a delta shape with aspect ratio of 1.

De Silva, B. M. E.↗

First-and Second-Order Displacement Transfer Functions for Structural Shape Calculations Using Analytically Predicted Surface Strains

New first- and second-order displacement transfer functions have been developed for deformed shape calculations of nonuniform cross-sectional beam structures such as aircraft wings. The displacement transfer functions are expressed explicitly in terms of beam geometrical parameters and surface strains (uniaxial bending strains) obtained at equally spaced strain stations along the surface of the beam structure. By inputting the measured or analytically calculated surface strains into the displacement transfer functions, one could calculate local slopes, deflections, and cross-sectional twist angles of the nonuniform beam structure for mapping the overall structural deformed shapes for visual display. The accuracy of deformed shape calculations by the first- and second-order displacement transfer functions are determined by comparing these values to the analytically predicted values obtained from finite element analyses. This comparison shows that the new displacement transfer functions could quite accurately calculate the deformed shapes of tapered cantilever tubular beams with different tapered angles. The accuracy of the present displacement transfer functions also are compared to those of the previously developed displacement transfer functions.

Ko, William L.↗

Multibody dynamics: Modeling component flexibility with fixed, free, loaded, constraint, and residual modes

The assumed-modes method in multibody dynamics allows the elastic deformation of each component in the system to be approximated by a sum of products of spatial and temporal functions commonly known as modes and modal coordinates respectively. The choice of component modes used to model articulating and non-articulating flexible multibody systems is examined. Attention is directed toward three classical Component Mode Synthesis (CMS) methods whereby component normal modes are generated by treating the component interface (I/F) as either fixed, free, or loaded with mass and stiffness contributions from the remaining components. The fixed and free I/F normal modes are augmented by static shape functions termed constraint and residual modes respectively. A mode selection procedure is outlined whereby component modes are selected from the Craig-Bampton (fixed I/F plus constraint), MacNeal-Rubin (free I/F plus residual), or Benfield-Hruda (loaded I/F) mode sets in accordance with a modal ordering scheme derived from balance realization theory. The success of the approach is judged by comparing the actuator-to-sensor frequency response of the reduced order system with that of the full order system over the frequency range of interest. A finite element model of the Galileo spacecraft serves as an example in demonstrating the effectiveness of the proposed mode selection method.

Spanos, John T.↗

Functionally Graded Adhesives for Composite Joints

Adhesives with functionally graded material properties are being considered for use in adhesively bonded joints to reduce the peel stress concentrations located near adherend discontinuities. Several practical concerns impede the actual use of such adhesives. These include increased manufacturing complications, alterations to the grading due to adhesive flow during manufacturing, and whether changing the loading conditions significantly impact the effectiveness of the grading. An analytical study is conducted to address these three concerns. An enhanced joint finite element, which uses an analytical formulation to obtain exact shape functions, is used to model the joint. Furthermore, proof of concept testing is conducted to show the potential advantages of functionally graded adhesives. In this study, grading is achieved by strategically placing glass beads within the adhesive layer at different densities along the joint.

Stapleton, Scott E.↗

PROFILE: Airfoil Geometry Manipulation and Display

This report provides user information for program PROFILE, an aerodynamics design utility for plotting, tabulating, and manipulating airfoil profiles. A dozen main functions are available. The theory and implementation details for two of the more complex options are also presented. These are the REFINE option, for smoothing curvature in selected regions while retaining or seeking some specified thickness ratio, and the OPTIMIZE option, which seeks a specified curvature distribution. Use of programs QPLOT and BPLOT is also described, since all of the plots provided by PROFILE (airfoil coordinates, curvature distributions, pressure distributions)) are achieved via the general-purpose QPLOT utility. BPLOT illustrates (again, via QPLOT) the shape functions used by two of PROFILE's options. These three utilities should be distributed as one package. They were designed and implemented for the Applied Aerodynamics Branch at NASA Ames Research Center, Moffett Field, California. They are all written in FORTRAN 77 and run on DEC and SGI systems under OpenVMS and IRIX.

Collins, Leslie↗

Modified Displacement Transfer Functions for Deformed Shape Predictions of Slender Curved Structures with Varying Curvatives

To eliminate the need to use finite-element modeling for structure shape predictions, a new method was invented. This method is to use the Displacement Transfer Functions to transform the measured surface strains into deflections for mapping out overall structural deformed shapes. The Displacement Transfer Functions are expressed in terms of rectilinearly distributed surface strains, and contain no material properties. This report is to apply the patented method to the shape predictions of non-symmetrically loaded slender curved structures with different curvatures up to a full circle. Because the measured surface strains are not available, finite-element analysis had to be used to analytically generate the surface strains. Previously formulated straight-beam Displacement Transfer Functions were modified by introducing the curvature-effect correction terms. Through single-point or dual-point collocations with finite-elementgenerated deflection curves, functional forms of the curvature-effect correction terms were empirically established. The resulting modified Displacement Transfer Functions can then provide quite accurate shape predictions. Also, the uniform straight-beam Displacement Transfer Function was applied to the shape predictions of a section-cut of a generic capsule (GC) outer curved sandwich wall. The resulting GC shape predictions are quite accurate in partial regions where the radius of curvature does not change sharply.

displacement transfer functions↗

Image Correlation Pattern Optimization for Micro-Scale In-Situ Strain Measurements

The accuracy and precision of digital image correlation (DIC) is a function of three primary ingredients: image acquisition, image analysis, and the subject of the image. Development of the first two (i.e. image acquisition techniques and image correlation algorithms) has led to widespread use of DIC; however, fewer developments have been focused on the third ingredient. Typically, subjects of DIC images are mechanical specimens with either a natural surface pattern or a pattern applied to the surface. Research in the area of DIC patterns has primarily been aimed at identifying which surface patterns are best suited for DIC, by comparing patterns to each other. Because the easiest and most widespread methods of applying patterns have a high degree of randomness associated with them (e.g., airbrush, spray paint, particle decoration, etc.), less effort has been spent on exact construction of ideal patterns. With the development of patterning techniques such as microstamping and lithography, patterns can be applied to a specimen pixel by pixel from a patterned image. In these cases, especially because the patterns are reused many times, an optimal pattern is sought such that error introduced into DIC from the pattern is minimized. DIC consists of tracking the motion of an array of nodes from a reference image to a deformed image. Every pixel in the images has an associated intensity (grayscale) value, with discretization depending on the bit depth of the image. Because individual pixel matching by intensity value yields a non-unique scale-dependent problem, subsets around each node are used for identification. A correlation criteria is used to find the best match of a particular subset of a reference image within a deformed image. The reader is referred to references for enumerations of typical correlation criteria. As illustrated by Schreier and Sutton and Lu and Cary systematic errors can be introduced by representing the underlying deformation with under-matched shape functions. An important implication, as discussed by Sutton et al., is that in the presence of highly localized deformations (e.g., crack fronts), error can be reduced by minimizing the subset size. In other words, smaller subsets allow the more accurate resolution of localized deformations. Contrarily, the choice of optimal subset size has been widely studied and a general consensus is that larger subsets with more information content are less prone to random error. Thus, an optimal subset size balances the systematic error from under matched deformations with random error from measurement noise. The alternative approach pursued in the current work is to choose a small subset size and optimize the information content within (i.e., optimizing an applied DIC pattern), rather than finding an optimal subset size. In the literature, many pattern quality metrics have been proposed, e.g., sum of square intensity gradient (SSSIG), mean subset fluctuation, gray level co-occurrence, autocorrelation-based metrics, and speckle-based metrics. The majority of these metrics were developed to quantify the quality of common pseudo-random patterns after they have been applied, and were not created with the intent of pattern generation. As such, it is found that none of the metrics examined in this study are fit to be the objective function of a pattern generation optimization. In some cases, such as with speckle-based metrics, application to pixel by pixel patterns is ill-conditioned and requires somewhat arbitrary extensions. In other cases, such as with the SSSIG, it is shown that trivial solutions exist for the optimum of the metric which are ill-suited for DIC (such as a checkerboard pattern). In the current work, a multi-metric optimization method is proposed whereby quality is viewed as a combination of individual quality metrics. Specifically, SSSIG and two auto-correlation metrics are used which have generally competitive objectives. Thus, each metric could be viewed as a constraint imposed upon the others, thereby precluding the achievement of their trivial solutions. In this way, optimization produces a pattern which balances the benefits of multiple quality metrics. The resulting pattern, along with randomly generated patterns, is subjected to numerical deformations and analyzed with DIC software. The optimal pattern is shown to outperform randomly generated patterns.

Bomarito, G. F.↗

In-orbit flexible spacecraft dynamics program

A continuous flexible body nonlinear dynamics computer program is used for simulating the spinning mode performance of a spacecraft under applied control torques. The program takes into account the continuous flexible nature of the antennas by representing deflections in terms of shape functions and integrated the spatial dependence in the formulation of the equations of motion. Comparison of RAE flight data on roll, pitch, and yaw with predictions of the computer program in the gravity gradient mode show fair to good agreement in roll and pitch and excellent agreement in the yaw angle.

Fedor, J. V.↗

A computerized symbolic integration technique for development of triangular and quadrilateral composite shallow-shell finite elements

Computerized symbolic integration was used in conjunction with group-theoretic techniques to obtain analytic expressions for the stiffness, geometric stiffness, consistent mass, and consistent load matrices of composite shallow shell structural elements. The elements are shear flexible and have variable curvature. A stiffness (displacement) formulation was used with the fundamental unknowns consisting of both the displacement and rotation components of the reference surface of the shell. The triangular elements have six and ten nodes; the quadrilateral elements have four and eight nodes and can have internal degrees of freedom associated with displacement modes which vanish along the edges of the element (bubble modes). The stiffness, geometric stiffness, consistent mass, and consistent load coefficients are expressed as linear combinations of integrals (over the element domain) whose integrands are products of shape functions and their derivatives. The evaluation of the elemental matrices is divided into two separate problems - determination of the coefficients in the linear combination and evaluation of the integrals. The integrals are performed symbolically by using the symbolic-and-algebraic-manipulation language MACSYMA. The efficiency of using symbolic integration in the element development is demonstrated by comparing the number of floating-point arithmetic operations required in this approach with those required by a commonly used numerical quadrature technique.

Anderson, C. M.↗

A general algorithm using finite element method for aerodynamic configurations at low speeds

A finite element algorithm for numerical simulation of two-dimensional, incompressible, viscous flows was developed. The Navier-Stokes equations are suitably modelled to facilitate direct solution for the essential flow parameters. A leap-frog time differencing and Galerkin minimization of these model equations yields the finite element algorithm. The finite elements are triangular with bicubic shape functions approximating the solution space. The finite element matrices are unsymmetrically banded to facilitate savings in storage. An unsymmetric L-U decomposition is performed on the finite element matrices to obtain the solution for the boundary value problem.

Balasubramanian, R.↗

A finite element simulation of sound attenuation in a finite duct with a peripherally variable liner

Using multimodal analysis, a variational finite element method is presented for analyzing sound attenuation in a three-dimensional finite duct with a peripherally variable liner in the absence of flow. A rectangular element, with cubic shaped functions, is employed. Once a small portion of a peripheral liner is removed, the attenuation rate near the frequency where maximum attenuation occurs drops significantly. The positioning of the liner segments affects the attenuation characteristics of the liner. Effects of the duct termination are important in the low frequency ranges. The main effect of peripheral variation of the liner is a broadening of the attenuation characteristics in the midfrequency range. Because of matrix size limitations of the presently available computer program, the eigenvalue equations should be solved out of core in order to handle realistic sources.

Watson, W. R.↗

Reproducibility of structural strength and stiffness for graphite-epoxy aircraft spoilers

Structural strength reproducibility of graphite epoxy composite spoilers for the Boeing 737 aircraft was evaluated by statically loading fifteen spoilers to failure at conditions simulating aerodynamic loads. Spoiler strength and stiffness data were statistically modeled using a two parameter Weibull distribution function. Shape parameter values calculated for the composite spoiler strength and stiffness were within the range of corresponding shape parameter values calculated for material property data of composite laminates. This agreement showed that reproducibility of full scale component structural properties was within the reproducibility range of data from material property tests.

Howell, W. E.↗

Proton-driven electromagnetic instabilities in high-speed solar wind streams

Electromagnetic instabilities of the field-aligned, right-hand circularly polarized magnetosonic wave and the left-hand circularly polarized Alfven wave driven by two drifted proton components are analyzed for model parameters determined from Imp 7 solar wind proton data measured during high-speed flow conditions. Growth rates calculated using bi-Lorentzian forms for the main and beam proton as well as core and halo electron velocity distributions do not differ significantly from those calculated using bi-Maxwellian forms. Using distribution parameters determined from 17 measured proton spectra, we show that considering the uncertainties the magnetosonic wave may be linearly stable and the Alfven wave is linearly unstable. Because proton velocity distribution function shapes are observed to persist for times long compared to the proton gyroperiod, the latter result suggests that linear stability theory fails for proton-driven ion cyclotron waves in the high-speed solar wind.

Abraham-Shrauner, B.↗