Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stiffness”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Reliability enhancement of Navier-Stokes codes through convergence acceleration

Methods for enhancing the reliability of Navier-Stokes computer codes through improving convergence characteristics are presented. The improving of these characteristics decreases the likelihood of code unreliability and user interventions in a design environment. The problem referred to as a 'stiffness' in the governing equations for propulsion-related flowfields is investigated, particularly in regard to common sources of equation stiffness that lead to convergence degradation of CFD algorithms. Von Neumann stability theory is employed as a tool to study the convergence difficulties involved. Based on the stability results, improved algorithms are devised to ensure efficient convergence in different situations. A number of test cases are considered to confirm a correlation between stability theory and numerical convergence. The examples of turbulent and reacting flow are presented, and a generalized form of the preconditioning matrix is derived to handle these problems, i.e., the problems involving additional differential equations for describing the transport of turbulent kinetic energy, dissipation rate and chemical species. Algorithms for unsteady computations are considered. The extension of the preconditioning techniques and algorithms derived for Navier-Stokes computations to three-dimensional flow problems is discussed. New methods to accelerate the convergence of iterative schemes for the numerical integration of systems of partial differential equtions are developed, with a special emphasis on the acceleration of convergence on highly clustered grids.

Merkle, Charles L.↗
NASA NTRS1995

Formulation of an improved smeared stiffener theory for buckling analysis of grid-stiffened composite panels

A smeared stiffener theory for stiffened panels is presented that includes skin-stiffener interaction effects. The neutral surface profile of the skin-stiffener combination is developed analytically using the minimum potential energy principle and statics conditions. The skin-stiffener interaction is accounted for by computing the stiffness due to the stiffener and the skin in the skin-stiffener region about the neutral axis at the stiffener. Buckling load results for axially stiffened, orthogrid, and general grid-stiffened panels are obtained using the smeared stiffness combined with a Rayleigh-Ritz method and are compared with results from detailed finite element analyses.

Jaunky, Navin↗
NASA NTRS1994

A survey of the core-congruential formulation for geometrically nonlinear TL finite elements

This article presents a survey of the core-congruential formulation (CCF) for geometrically nonlinear mechanical finite elements based on the total Lagrangian (TL) kinematic description. Although the key ideas behind the CCF can be traced back to Rajasekaran and Murray in 1973, it has not subsequently received serious attention. The CCF is distinguished by a two-phase development of the finite element stiffness equations. The initial phase developed equations for individual particles. These equations are expressed in terms of displacement gradients as degrees of freedom. The second phase involves congruential-type transformations that eventually binds the element particles of an individual element in terms of its node-displacement degrees of freedom. Two versions of the CCF, labeled direct and generalized, are distinguished. The direct CCF (DCCF) is first described in general form and then applied to the derivation of geometrically nonlinear bar, and plane stress elements using the Green-Lagrange strain measure. The more complex generalized CCF (GCCF) is described and applied to the derivation of 2D and 3D Timoshenko beam elements. Several advantages of the CCF, notably the physically clean separation of material and geometric stiffnesses, and its independence with respect to the ultimate choice of shape functions and element degrees of freedom, are noted. Application examples involving very large motions solved with the 3D beam element display the range of applicability of this formulation, which transcends the kinematic limitations commonly attributed to the TL description.

Felippa, Carlos A.↗
NASA NTRS1995

A study of elevated temperature testing techniques for the fatigue behavior of PMCS: Application to T650-35/AMB21

An experimental study was conducted to investigate the mechanical behavior of a T650-35/AMB21 eight-harness satin weave polymer composite system. Emphasis was placed on the development and refinement of techniques used in elevated temperature uniaxial PMC testing. Issues such as specimen design, gripping, strain measurement, and temperature control and measurement were addressed. Quasi-static tensile and fatigue properties (R(sub sigma) = 0.1) were examined at room and elevated temperatures. Stiffness degradation and strain accumulation during fatigue cycling were recorded to monitor damage progression and provide insight for future analytical modeling efforts. Accomplishments included an untabbed dog-bone specimen design which consistently failed in the gage section, accurate temperature control and assessment, and continuous in-situ strain measurement capability during fatigue loading at elevated temperatures. Finally, strain accumulation and stiffness degradation during fatigue cycling appeared to be good indicators of damage progression.

Gyekenyesi, Andrew L.↗
NASA NTRS1992

A three-dimensional nonlinear Timoshenko beam based on the core-congruential formulation

A three-dimensional, geometrically nonlinear two-node Timoshenkoo beam element based on the total Larangrian description is derived. The element behavior is assumed to be linear elastic, but no restrictions are placed on magnitude of finite rotations. The resulting element has twelve degrees of freedom: six translational components and six rotational-vector components. The formulation uses the Green-Lagrange strains and second Piola-Kirchhoff stresses as energy-conjugate variables and accounts for the bending-stretching and bending-torsional coupling effects without special provisions. The core-congruential formulation (CCF) is used to derived the discrete equations in a staged manner. Core equations involving the internal force vector and tangent stiffness matrix are developed at the particle level. A sequence of matrix transformations carries these equations to beam cross-sections and finally to the element nodal degrees of freedom. The choice of finite rotation measure is made in the next-to-last transformation stage, and the choice of over-the-element interpolation in the last one. The tangent stiffness matrix is found to retain symmetry if the rotational vector is chosen to measure finite rotations. An extensive set of numerical examples is presented to test and validate the present element.

Crivelli, Luis A.↗
NASA NTRS1994

Characterization of path effects on the properties of balloon films

The effect of stress and temperature path on the mechanical properties of a balloon film material (Astrofilm E2) have been studied using two techniques. The first technique is based on uniaxial tensile measurements of prestrained strips of the film as a function of temperature. Data are analyzed in terms of two temperature coefficients, defined here as a stiffness coefficient and a strength coefficient; and by comparing the areas under the stress-strain curves. The second technique investigated the effect of stress and temperature path on inflated cylinders of the same material. The changes in the mechanical properties due to preconditioning in the biaxial stress state were evaluated using the ball burst test. Preliminary findings indicate that the material is stress and temperature path dependent. Changes in stiffness coefficient, draw strength and toughness were measured and discussed. The above techniques may be suitable to document path dependent changes, and to discriminate among films of different manufacture.

Said, Magdi A.↗
NASA NTRS1994

Buckling and vibration analysis of laminated panels using VICONOPT

The analysis aspects of the 23,000-line FORTRAN program VIPASA with constraints and optimization (VICONOPT) are described. Overall stiffness matrices assembled from the earlier exact vibration and instability of plate assemblies including shear and anisotropy (VIPASA) flat plate stiffnesses are optionally coupled by Langrangian multipliers to find critical buckling loads, or natural frequencies of undamped vibration, of prismatic assemblies of anisotropic flat plates with arbitrarily located point supports or simple transverse supporting frames. The longitudunal continuity of typical wing and fuselage panels is closely approximated because the solutions are for the infinitely long structure obtained by repeating a bay and its supports longitudinally. Any longitudinally invariant in-plane plate stresses are permitted, and very rapid solutions are guaranteed by numerous refinements, including multilevel substructuring and a method for repetitive cross sections that is exact for regular polygons used to represent cylinders. Modal displacements and stresses in or between plies of laminated plates are calculated and plotted, with values being recovered at all nodes of substructures. Comparison with usual approximate finite-element methods confirms that, for comparably converged solutions, VICONOPT is typically between 100 and 10(exp 4) times faster.

Kennedy, David↗
NASA NTRS1993

Identification of a second-order mechanical system model from a state space realization

Many system identification algorithms produce models in modern state variable form. When a mechanical system is identified one knows that there must be a model in second order form. In structural dynamics one is interested in either a modal model of the system, or in a model including a mass matrix, stiffness matrix and a damping matrix. In this paper, algorithms are developed that convert a modern state space realization into the above two representations of interest in mechanical systems. The algorithm can be used with any identification methods that produce a modern state variable representation such as the Eigensystem Realization Algorithm (ERA), its modified version by data correlation (ERA/DC), and combinations of these with Observer/Kalman Filter Identification (OKID). An algorithm is developed that allows one to identify the damping matrix in a model representation, so that one determines how near to modal the damping is in the system, and can therefore understand the degree to which the modes are coupled by the damping in the system. A second algorithm produced the mass, damping, and stiffness matrices, given the input and output matrices. It is assumed that the number of sensors (or the number of actuators) is greater than or equal to the number of modes in the system, and it is proved that this is a necessary condition to be able to uniquely identify these matrices. Experience with examples suggests that a singular value truncation involved in these algorithms can help one determine the true system order. This truncation has the benefit of knowledge of the correct form for a mechanical system model, which is not present in the modern control identification algorithm.

Tseng, Dong-Huei↗
NASA NTRS1994

Equivalent plate structural modeling for wing shape optimization including transverse shear

A new technique for structural modeling of airplanes wings is presented taking transverse shear effects into account. The kinematic assumptions of first-order shear deformation plate theory in combination with numerical analysis, where simple polynomials are used to define geometry, construction, and displacement approximations, lead to analytical expressions for elements of the stiffness and mass matrices and load vector. Contributions from the cover skins, spar and rib caps, and spar and rib webs are included as well as concentrated springs and concentrated masses. Limitations of wing modeling techniques based on classical plate theory are discussed, and the improved accuracy of the new equivalent plate technique is demonstrated through comparison with finite element analysis and test results. Expressions for analytical derivatives of stiffness, mass, and load terms with respect to wing shape are given. Based on these, it is possible to obtain analytic sensitivities of displacements, stresses, and natural frequencies with respect to planform shape and depth distribution. This makes the new capability an effective structural tool for wing shape optimization.

Livne, Eli↗
NASA NTRS1995

Buffeting tests in a cryogenic windtunnel

Measurements of wing buffeting, using root strain gages, were made in the NASA Langley 0.3 m cryogenic wind tunnel to refine techniques which will be used in larger cryogenic facilities such as the United States National Transonic Facility (NTF) and the European Transonic Wind Tunnel (ETW). The questions addressed included the relative importance variations in frequency parameter and Reynolds number, the choice of model material (considering both stiffness and damping) and the effects of static aeroelastic distortion. The main series of tests was made on three half models of slender 65 deg delta wings with a sharp leading edge. The three delta wings had the same planform but widely differing bending stiffnesses and frequencies (obtained by varying both the material and the thickness of the wings). It was known that the steady flow on this configuration would be insensitive to variations in Reynolds number. On this wing at vortex breakdown the spectrum of the unsteady excitation is unusual, having a sharp peak at particular frequency parameter. Additional tests were made on one unswept half-wing of aspect ratio 1.5 with an NPL 9510 aerofoil section, known to be sensitive to variations in Reynolds number at transonic speeds. The test Mach numbers were M = 0.21 and 0.35 for the delta wings and to M = 0.30 for the unswept wing. On this wing the unsteady excitation spectrum is fairly flat (as on most wings). Hence correct representation of the frequency parameter is not particularly important.

Mabey, D. G.↗
NASA NTRS1995

Failure models for textile composites

The goals of this investigation were to: (1) identify mechanisms of failure and determine how the architecture of reinforcing fibers in 3D woven composites controlled stiffness, strength, strain to failure, work of fracture, notch sensitivity, and fatigue life; and (2) to model composite stiffness, strength, and fatigue life. A total of 11 different angle and orthogonal interlock woven composites were examined. Composite properties depended on the weave architecture, the tow size, and the spatial distributions and strength of geometrical flaws. Simple models were developed for elastic properties, strength, and fatigue life. A more complicated stochastic model, called the 'Binary Model,' was developed for damage tolerance and ultimate failure. These 3D woven composites possessed an extraordinary combination of strength, damage tolerance, and notch insensitivity.

Cox, Brian↗
NASA NTRS1992

Passive damping augmentation of flexible beam-like lattice trusses for large space structures

A Timoshenko beam continuum model is developed to determine the transient response of a beam-like latticed large space structure, subjected to a unit impulse. It is demonstrated that an increase in diagonal stiffness, on account of the stiffness of the vertical girder, leads to a rise in the transverse shear rigidity. This results in higher natural frequencies and reduction in peak displacement. In addition, in an asymmetrical truss configuration, coupling between the extensional and shear modes raises the maximum peak displacement compared with that obtained for a symmetrical truss. The model is modified to investigate the introduction of passive damping in the form of several dynamic vibration absorbers. A unique absorber parameter optimization procedure, based on the classical steady-state criteria, is used to tune the absorbers, having a total allocated mass budget of 10%, to the first two frequencies of the latticed structure. It is found that inclusion of transverse shear rigidity as a design parameter in damping augmentation studies reduces settling time for predetermined maximum peak displacements.

Cohen, K.↗
NASA NTRS1995

A comparison of graphite/epoxy tape laminates and 2-D braided composites mechanical properties

A comparison of the mechanical properties of unidirectional composite tape laminates and of two-dimensional triaxially braided composite was conducted. The tape laminate layups were designed to match the percentage of axial fibers and the angle of the bias tows in the braided composite. The materials system used for the laminates is AS4/3501-6 which was chosen as the closest available match to As4/1895 used for the braids. The strength and stiffness properties measured here include tension, open-hole tension, filled-hole tension, compression and open-hole compression, all of these in both the longitudinal and transverse direction. Results show that the longitudinal modulus of both material forms is quite similar, but that the transverse modulus of the braids is lower. In terms of strength, the longitudinal unnotched strength of the braids is lower than that of the laminates, while the transverse strength is significantly lower. For both strength and stiffness, the crimp in the bias tows of the braid is probably the main cause for reduced properties. On the other hand, a very significant increase in open-hole and filled-hole tension strength was observed for the braids compared to the tape laminates. However, this was not observed in compression where all the braid properties are lower than for the laminates.

Minguet, Pierre J.↗
NASA NTRS1995

TEXCAD: TEXile Composite Analysis for Design

The Textile Composite Analysis for Design (TEXCAD) code provides the materials/design engineer with a user-friendly, desktop computer based tool for the analysis of a wide variety of fabric reinforced woven and braided composites. It can be used to calculate overall thermal and mechanical properties along with engineering estimates of damage progression and strength. TEXCAD also calculates laminate properties for stacked, oriented fabric constructions. It discretely models the yarn centerline paths within the textile repeating unit cell (RUC) by assuming sinusoidal undulations at yarn cross-over points and uses a year discretization scheme (which subdivides each yarn into smaller, piecewise straight yarn slices) together with a 3-D stress averaging procedure to compute overall stiffness properties. In the calculations for strength, it uses a curved beam-on-elastic foundation model for yarn undulating regions together with incremental approach in which stiffness properties for the failed yarn slices are reduced based on the predicted yarn slice failure mode. Nonlinear shear effects and nonlinear geometric effects can be stimulated. Input to TEXCAD consists of: (1) material parameters like impregnated yarn and resin properties such as moduli, Poisson's ratios, coefficients of thermal expansion, nonlinear shear parameters, axial failure strains, and in-plane failure stresses; and (2) fabric parameters like yarn sizes, braid angle, yarn packing density, filament diameter, and overall fiber volume fraction. Output consists of overall thermoelastic constants, yarn slice strains/stresses, yarn slice failure history, in-plane strain response, and ultimate failure strength. Strength can be computed under the combined action of thermal and mechanical loading (tension, compression, and shear). A brief overview of the analytical capabilities, program organization, and modules, input and output parameters, computer platforms, distribution, and modifications/extensions of the TEXCAD code are presented here.

Naik, Rajiv A.↗
NASA NTRS1995

A new merit function for evaluating the flaw tolerance of composite laminates

Advanced composite materials are providing notable savings in numerous aerospace and non-aerospace applications by virtue of their uniquely high strength-to-density ratio and high stiffness-to-density ratio. However, their full potential for high performance has yet to be achieved because of the low tolerance of the 'brittle' material to impact, holes, flaws, and localized discontinuities due to relatively low toughness. In the present paper a new and simple merit function is proposed for maximizing the strength performance of flawed composite laminate for different design requirement. This new merit function takes into account the flaw tolerance of a composite laminated as well as its stiffness. After development of the new merit function, an example is presented of its application to a practical laminate family.

Mikulas, Martin M., Jr.↗
NASA NTRS1995

Angled injection: Hybrid fluid film bearings for cryogenic applications

A computational bulk-flow analysis for prediction of the force coefficients of hybrid fluid film bearings with angled orifice injection is presented. Past measurements on water-lubricated hybrid bearings with angle orifice injection have demonstrated improved rotordynamic performance with virtual elimination of cross-coupled stiffness coefficients and nul or negative whirl frequency ratios. A simple analysis reveals that the fluid momentum exchange at the orifice discharge produces a pressure rise in the recess which retards the shear flow induced by journal rotation, and consequently, reduces cross-coupling forces. The predictions from the model correlate well with experimental measurements from a radial and 45 deg angled orifice injection, five recess water hybrid bearings (C = 125 microns) operating at 10.2, 17.4, and 24.6 krpm and with nominal supply pressures equal to 4, 5.5, and 7 MPa. An application example for a liquid oxygen six recess/pad hybrid journal bearing shows the advantages of tangential orifice injection on the rotordynamic force coefficients and stability indicator for forward whirl motions and without performance degradation on direct stiffness and damping coefficients. The computer program generated, 'hydrojet,' extends and complements previously developed codes.

SanAndres, Luis↗
NASA NTRS1996

Quantification of Gear Tooth Damage by Optimal Tracking of Vibration Signatures

This paper presents a technique for quantifying the wear or damage of gear teeth in a transmission system. The procedure developed in this study can be applied as a part of either an onboard machine health-monitoring system or a health diagnostic system used during regular maintenance. As the developed methodology is based on analysis of gearbox vibration under normal operating conditions, no shutdown or special modification of operating parameters is required during the diagnostic process. The process of quantifying the wear or damage of gear teeth requires a set of measured vibration data and a model of the gear mesh dynamics. An optimization problem is formulated to determine the profile of a time-varying mesh stiffness parameter for which the model output approximates the measured data. The resulting stiffness profile is then related to the level of gear tooth wear or damage. The procedure was applied to a data set generated artificially and to another obtained experimentally from a spiral bevel gear test rig. The results demonstrate the utility of the procedure as part of an overall health-monitoring system.

Choy, F. K.↗
NASA NTRS1995

Interactive flutter analysis and parametric study for conceptual wing design

An interactive computer program was developed for wing flutter analysis in the conceptual design stage. The objective was to estimate the flutter instability boundary of a flexible cantilever wing, when well defined structural and aerodynamic data are not available, and then study the effect of change in Mach number, dynamic pressure, torsional frequency, sweep, mass ratio, aspect ratio, taper ratio, center of gravity, and pitch inertia, to guide the development of the concept. The software was developed on MathCad (trademark) platform for Macintosh, with integrated documentation, graphics, database and symbolic mathematics. The analysis method was based on nondimensional parametric plots of two primary flutter parameters, namely Regier number and Flutter number, with normalization factors based on torsional stiffness, sweep, mass ratio, aspect ratio, center of gravity location and pitch inertia radius of gyration. The plots were compiled in a Vaught Corporation report from a vast database of past experiments and wind tunnel tests. The computer program was utilized for flutter analysis of the outer wing of a Blended Wing Body concept, proposed by McDonnell Douglas Corporation. Using a set of assumed data, preliminary flutter boundary and flutter dynamic pressure variation with altitude, Mach number and torsional stiffness were determined.

Mukhopadhyay, Vivek↗