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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

A waved journal bearing concept with improved steady-state and dynamic performance

Analysis of the waved journal bearing concept featuring a waved inner bearing diameter for use with a compressible lubricant (gas) is presented. A three wave, waved journal bearing geometry is used to show the geometry of this concept. The performance of generic waved bearings having either three, four, six, or eight waves is predicted for air lubricated bearings. Steady-state performance is discussed in terms of bearing load capacity, while the dynamic performance is discussed in terms of dynamic coefficients and fluid film stability. It was found that the bearing wave amplitude has an important influence on both steady-state and dynamic performance of the waved journal bearing. For a fixed eccentricity ratio, the bearing steady-state load capacity and direct dynamic stiffness coefficient increase as the wave amplitude increases. Also, the waved bearing becomes more stable as the wave amplitude increases. In addition, increasing the number of waves reduces the waved bearing's sensitivity to the direction of the applied load relative to the wave. However, the range in which the bearing performance can be varied decreases as the number of waves increases. Therefore, both the number and the amplitude of the waves must be properly selected to optimize the waved bearing design for a specific application. It is concluded that the stiffness of an air bearing, due to the hydrodynamic effect, could be doubled and made to run stably by using a six or eight wave geometry with a wave amplitude approximately half of the bearing radial clearance.

Dimofte, Florin↗
NASA NTRS1994

Isothermal Damage and Fatigue Behavior of SCS-6/Timetal 21S [0/90](Sub S) Composite at 650 Deg C

The isothermal fatigue damage and life behaviors of SCS-6/Timetal 21S (0/90)s were investigated at 650 C. Strain ratcheting and degradation of the composite's static elastic modulus were carefully monitored as functions of cycles to indicate damage progression. Extensive fractographic and metallographic analyses were conducted to determine damage/failure mechanisms. Resulting fatigue lives show considerable reductions in comparison to (0) reinforced titanium matrix composites subjected to comparable conditions. Notable stiffness degradations were found to occur after the first cycle of loading, even at relatively low maximum stress levels, where cyclic lives are greater than 25,000 cycles. This was attributed to the extremely weak fiber/matrix bond which fails under relatively low transverse loads. Stiffness degradations incurred on first cycle loadings and degradations thereafter were found to increase with increasing maximum stress. Environmental effects associated with oxidation of the (90) fiber interfaces clearly played a role in the damage mechanisms as fracture surfaces revealed environment assisted matrix cracking along the (90) fibers. Metallographic analysis indicated that all observable matrix fatigue cracks initiated at the (90) fiber/matrix interfaces. Global de-bonding in the loading direction was found along the (90) fibers. No surface initiated cracks were evident and minimal if any (0) fiber cracking was visible.

Castelli, Michael G.↗
NASA NTRS1994

Scaling, elasticity, and CLPT

The first few viewgraphs describe the general solution properties of linear elasticity theory which are given by the following two statements: (1) for stress B.C. on S(sub sigma) and zero displacement B.C. on S(sub u) the altered displacements u(sub i)(*) and the actual stresses tau(sub ij) are elastically dependent on Poisson's ratio nu alone: thus the actual displacements are given by u(sub i) = mu(exp -1)u(sub i)(*); and (2) for zero stress B.C. on S(sub sigma) and displacement B.C. on S(sub u) the actual displacements u(sub i) and the altered stresses tau(sub ij)(*) are elastically dependent on Poisson's ratio nu alone: thus the actual stresses are given by tau(sub ij) = E tau(sub ij)(*). The remaining viewgraphs describe the minimum parameter formulation of the general classical laminate theory plate problem as follows: The general CLT plate problem is expressed as a 3 x 3 system of differential equations in the displacements u, v, and w. The eighteen (six each) A(sub ij), B(sub ij), and D(sub ij) system coefficients are ply-weighted sums of the transformed reduced stiffnesses (bar-Q(sub ij))(sub k); the (bar-Q(sub ij))(sub k) in turn depend on six reduced stiffnesses (Q(sub ij))(sub k) and the material and geometry properties of the k(sup th) layer. This paper develops a method for redefining the system coefficients, the displacement components (u,v,w), and the position components (x,y) such that a minimum parameter formulation is possible. The pivotal steps in this method are (1) the reduction of (bar-Q(sub ij))(sub k) dependencies to just two constants Q(*) = (Q(12) + 2Q(66))/(Q(11)Q(22))(exp 1/2) and F(*) - (Q(22)/Q(11))(exp 1/2) in terms of ply-independent reference values Q(sub ij); (2) the reduction of the remaining portions of the A, B, and D coefficients to nondimensional ply-weighted sums (with 0 to 1 ranges) that are independent of Q(*) and F(*); and (3) the introduction of simple coordinate stretchings for u, v, w and x,y such that the process is neatly completed.

Brunelle, Eugene J.↗
NASA NTRS1992

On a finite-state inflow application to flap-lag-torsion damping in hover

An aerodynamic model with a coupled set of generalized dynamic wake equations and hybrid equations of motion for an elastic blade are applied here to a two-blade untwisted stiff in-plane hingeless small-scale model rotor with torsionally soft blades. Blade root offset, precone, blade droop, pitch control stiffness, and blade pitch angle are included in the model rotor. Numerical results show that 3D tip relief effects within the nonuniform steady-state inflow are significant to predict steady-state aerodynamic loads and blade deflections. Eigenvalue results confirm the importance of unsteady 3D aerodynamics in predicting lead-lag damping and frequency. Eigenvector analysis correlations reinforced qualitative and quantitative shortcomings associated with quasi-steady 2D aerodynamic theory for aeroelastic applications in hover.

De Andrade, Donizeti↗
NASA NTRS1994

Mechanical property characterization and impact resistance of selected graphite/PEEK composite materials

To use graphite polyetheretherketone (PEEK) material on highly curved surfaces requires that the material be drapable and easily conformable to the surface. This paper presents the mechanical property characterization and impact resistance results for laminates made from two types of graphite/PEEK materials that will conform to a curved surface. These laminates were made from two different material forms. These forms are: (1) a fabric where each yarn is a co-mingled Celion G30-500 3K graphite fiber and PEEK thermoplastic fiber; and (2) an interleaved material of Celion G30-500 3K graphite fabric interleaved with PEEK thermoplastic film. The experimental results from the fabric laminates are compared with results for laminates made from AS4/PEEK unidirectional tape. The results indicate that the tension and compression moduli for quasi-isotropic and orthotropic laminates made from fabric materials are at least 79 percent of the modulus of equivalent laminates made from tape material. The strength of fabric material laminates is at least 80 percent of laminates made from tape material. The evaluation of fabric material for shear stiffness indicates that a tape material laminate could be replaced by a fabric material laminate and still maintain 89 percent of the shear stiffness of the tape material laminate. The notched quasi-isotropic compression panel failure strength is 42 to 46 percent of the unnotched quasi-isotropic laminate strength. Damage area after impact with 20 ft-lbs of impact energy is larger for the co-mingled panels than for the interleaved panels. The inerleaved panels have less damage than panels made from tape material. Residual compression strength of quasi-isotropic panels after impact of 20 ft-lbs of energy varies between 33 percent of the undamaged quasi-isotropic material strength for the tape material and 38 percent of the undamaged quasi-isotropic material strength for the co-mingled fabric material.

Baker, Donald J.↗
NASA NTRS1994

Micromechanics of fatigue in woven and stitched composites

The goals of this research program were to: (1) determine how microstructural factors, especially the architecture of reinforcing fibers, control stiffness, strength, and fatigue life in 3D woven composites; (2) identify mechanisms of failure; (3) model composite stiffness; (4) model notched and unnotched strength; and (5) model fatigue life. We have examined a total of eleven different angle and orthogonal interlock woven composites. Extensive testing has revealed that these 3D woven composites possess an extraordinary combination of strength, damage tolerance, and notch insensitivity in compression and tension and in monotonic and cyclic loading. In many important regards, 3D woven composites far outstrip conventional 2D laminates or stitched laminates. Detailed microscopic analysis of damage has led to a comprehensive picture of the essential mechanisms of failure and how they are related to the reinforcement geometry. The critical characteristics of the weave architecture that promote favorable properties have been identified. Key parameters are tow size and the distributions in space and strength of geometrical flaws. The geometrical flaws should be regarded as controllable characteristics of the weave in design and manufacture. In addressing our goals, the simplest possible models of properties were always sought, in a blend of old and new modeling concepts. Nevertheless, certain properties, especially regarding damage tolerance, ultimate failure, and the detailed effects of weave architecture, require computationally intensive stochastic modeling. We have developed a new model, the 'binary model,' to carry out such tasks in the most efficient manner and with faithful representation of crucial mechanisms. This is the final report for contract NAS1-18840. It covers all work from April 1989 up to the conclusion of the program in January 1993.

Cox, B. N.↗
NASA NTRS1994

Analysis of the Effects of Surface Pitting and Wear on the Vibrations of a Gear Transmission System

A comprehensive procedure to simulate and analyze the vibrations in a gear transmission system with surface pitting, 'wear' and partial tooth fracture of the gear teeth is presented. An analytical model was developed where the effects of surface pitting and wear of the gear tooth were simulated by phase and magnitude changes in the gear mesh stiffness. Changes in the gear mesh stiffness were incorporated into each gear-shaft model during the global dynamic simulation of the system. The overall dynamics of the system were evaluated by solving for the transient dynamics of each shaft system simultaneously with the vibration of the gearbox structure. In order to reduce the number of degrees-of-freedom in the system, a modal synthesis procedure was used in the global transient dynamic analysis of the overall transmission system. An FFT procedure was used to transform the averaged time signal into the frequency domain for signature analysis. In addition, the Wigner-Ville distribution was also introduced to examine the gear vibration in the joint time frequency domain for vibration pattern recognition. Experimental results obtained from a gear fatigue test rig at NASA Lewis Research Center were used to evaluate the analytical model.

Choy, F. K.↗
NASA NTRS1994

Beta testing of MTI seal codes

An evaluation of MTI seal codes is made by comparing cylindrical air and water seals. Results are presented in viewgraph format and show that: ICYL and GCYL geometry variations are desirable; load and direct stiffness calculations are good; damping and cross-coupled stiffness predictions are poor; added mass coefficients should be calculated; and variation in inlet tangential velocity is critical to design.

Scharrer, Joseph↗
NASA NTRS1994

EG and G and NASA face seal codes comparison

This viewgraph presentation presents the following results for the example comparison: EG&G code with face deformations suppressed and SPIRALG agree well with each other as well as with the experimental data; 0 rpm stiffness data calculated by EG&G code are about 70-100 percent lower than that by SPIRALG; there is no appreciable difference between 0 rpm and 16,000 rpm stiffness and damping coefficients calculated by SPIRALG; and the film damping above 500 psig calculated by SPIRALG is much higher than the O-Ring secondary seal damping (e.g. 50 lbf.s/in).

Basu, Prit↗
NASA NTRS1994

AXAF SIM focus mechanism study

The design requirements and initial design concept for the AXAF-I Science Instrument Module (SIM) were reviewed at Ball on September 29, 1993. The concept design SIM focus mechanism utilizes a planetary gearset, with redundant motors, to drive a large ring (called 'main housing bearing') via a spur gearset. This large drive ring actuates three tangent bar links (called 'push rods'), which in turn actuate three levers (called 'pin levers'). Each of the three pin levers rotates an 'eccentric pin,' which in turn moves the base of a bipod flexure in both the radial (normal to optical axis) and axial (focus along optical axis) directions. Three bipod flexures are employed, equally spaced at 120 degrees apart, the base of each being translated in the two directions as described above. A focus adjustment is made by rotating the drive ring, which drives the push rods and therefore the pin levers, which in turn rotate the eccentric pins, finally imparting the two motions to the base of each of the bipod flexures. The axial translation (focus adjustment) of the focused structure is the sum of the direct axial motion plus axial motion which comes from uniformly squeezing the three bipod bases radially inward. SAO documented the following concerns regarding the focus mechanism in memo WAP-FY94-001, dated October 7, 1993: (1) The focus adjustment depends, in large part, on the structural properties (stiffnesses and end fixities) of the bipod flexures, push rods, pin levers and eccentric pins. If these properties are not matched very well, then lateral translations as well as unwanted rotations of the focussed structure will accompany focus motion. In addition, the stackup of linkage tolerances and any nonuniform wear in the linkages will result in the same unwanted motions. Thermal gradients will also affect these motions. At the review Ball did not present supporting analyses to support their choice of this design concept. (2) The proposed 'primary' method of measuring focus is by counting motor steps. The 'backup' method is by a pot mounted on the drive ring. Neither method provides for a direct measurement of the quantity desired (focus position). This is of concern because of the long and indirect relationship between focus and the sensed quantity (drive ring rotation). There are three sinusoidal relationships and structural stiffness in the path, and the resulting calibration is likely to be highly nonlinear. These methods would require an accurate ground calibration. (3) Ground calibration (and verification) of focus vs. drive position must be done in 1-g on the ground. This calibration will be complicated by both the structural characteristics of the bipods and the fact that the CG of the translating portion of the SIM is not on the optical axis (thereby causing unwated rotations and changing the focus position vs. motor step and pot readout relationships). The SIM translating weight could be offloaded, but the calibration then becomes sensitive to any errors in offloading (both magnitude and direction). There are concerns as to whether a calibration to the required accuracy can be accomplished on the ground. (4) The choice of a potentiometer as the focus position sensor is questionable in terms of reliability for a five year mission. The results of SAO's study of items 1, 2 and 3 described above are presented in this report.

Tananbaum, H. D.↗
NASA NTRS1994

The buckling response of symmetrically laminated composite plates having a trapezoidal planform area

The focus of this work is the buckling response of symmetrically laminated composite plates having a planform area in the shape of an isosceles trapezoid. The loading is assumed to be inplane and applied perpendicular to the parallel ends of the plate. The tapered edges of the plate are assumed to have simply supported boundary conditions, while the parallel ends are assumed to have either simply supported or clamped boundary conditions. A semi-analytic closed-form solution based on energy principles and the Trefftz stability criterion is derived and solutions are obtained using the Rayleigh-Ritz method. Intrinsic in this solution is a simplified prebuckling analysis which approximates the inplane force resultant distributions by the forms Nx=P/W(x) and Ny=Nxy=0, where P is the applied load and W(x) is the plate width which, for the trapezoidal planform, varies linearly with the lengthwise coordinate x. The out-of-plane displacement is approximated by a double trigonometric series. This analysis is posed in terms of four nondimensional parameters representing orthotropic and anisotropic material properties, and two nondimensional parameters representing geometric properties. For comparison purposes, a number of specific plate geometry, ply orientation, and stacking sequence combinations are investigated using the general purpose finite element code ABAQUS. Comparison of buckling coefficients calculated using the semi-analytical model and the finite element model show agreement within 5 percent, in general, and within 15 percent for the worst cases. In order to verify both the finite element and semi-analytical analyses, buckling loads are measured for graphite/epoxy plates having a wide range of plate geometries and stacking sequences. Test fixtures, instrumentation system, and experimental technique are described. Experimental results for the buckling load, the buckled mode shape, and the prebuckling plate stiffness are presented and show good agreement with the analytical results regarding the buckling load and the prebuckling plate stiffness. However, the experimental results show that for some cases the analysis underpredicts the number of halfwaves in the buckled mode shape. In the context of the definitions of taper ratio and aspect ratio used in this study, it is concluded that the buckling load always increases as taper ratio increases for a given aspect ratio for plates having simply supported boundary conditions on the parallel ends. There are combinations of plate geometry and ply stackling sequences, however, that reverse this trend for plates having clamped boundary conditions on the parallel ends such that an increase in the taper ratio causes a decrease in the buckling load. The clamped boundary conditions on the parallel ends of the plate are shown to increase the buckling load compared to simply supported boundary conditions. Also, anisotropy (the D16 and D26 terms) is shown to decrease the buckling load and skew the buckled mode shape for both the simply supported and clamped boundary conditions.

Radloff, H. D., II↗
NASA NTRS1994

Higher harmonic control analysis for vibration reduction of helicopter rotor systems

An advanced higher harmonic control (HHC) analysis has been developed and applied to investigate its effect on vibration reduction levels, blade and control system fatigue loads, rotor performance, and power requirements of servo-actuators. The analysis is based on a finite element method in space and time. A nonlinear time domain unsteady aerodynamic model, based on the indicial response formulation, is used to calculate the airloads. The rotor induced inflow is computed using a free wake model. The vehicle trim controls and blade steady responses are solved as one coupled solution using a modified Newton method. A linear frequency-domain quasi-steady transfer matrix is used to relate the harmonics of the vibratory hub loads to the harmonics of the HHC inputs. Optimal HHC is calculated from the minimization of the vibratory hub loads expressed in term of a quadratic performance index. Predicted vibratory hub shears are correlated with wind tunnel data. The fixed-gain HHC controller suppresses completely the vibratory hub shears for most of steady or quasi-steady flight conditions. HHC actuator amplitudes and power increase significantly at high forward speeds (above 100 knots). Due to the applied HHC, the blade torsional stresses and control loads are increased substantially. For flight conditions where the blades are stalled considerably, the HHC input-output model is quite nonlinear. For such cases, the adaptive-gain controller is effective in suppressing vibratory hub loads, even though HHC may actually increase stall areas on the rotor disk. The fixed-gain controller performs poorly for such flight conditions. Comparison study of different rotor systems indicates that a soft-inplane hingeless rotor requires less actuator power at high speeds (above 130 knots) than an articulated rotor, and a stiff-inplane hingeless rotor generally requires more actuator power than an articulated or a soft-inplane hingeless rotor. Parametric studies for a hingeless rotor operating in a transition flight regime and for an articulated rotor operating at the level-flight boundary (high speed and high thrust conditions) indicate that blade parameters including flap, lag, torsion stiffness distributions, linear pretwist, chordwise offset of center-of-mass from elastic axis and chordwise offset of elastic axis from aerodynamic center can be selected to minimize the actuator power requirements for HHC.

Nguyen, Khanh Q.↗
NASA NTRS1994

TEXCAD: Textile Composite Analysis for Design. Version 1.0: User's manual

The Textile Composite Analysis for Design (TEXCAD) code provides the materials/design engineer with a user-friendly desktop computer (IBM PC compatible or Apple Macintosh) 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 yarn discretization scheme (which subdivides each yarn not 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 an 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 simulated. Input to TEXCAD consists of: (1) materials parameters like impregnated yarn and resin properties such moduli, Poisson's ratios, coefficients of thermal expansion, nonlinear 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 stress-strain response and ultimate failure strength. Strength can be computed under the combined action of thermal and mechanical loading (tension, compression and shear).

Naik, Rajiv A.↗
NASA NTRS1994

Experimental verification of a progressive damage model for composite laminates based on continuum damage mechanics

Progressive failure is a crucial concern when using laminated composites in structural design. Therefore the ability to model damage and predict the life of laminated composites is vital. The purpose of this research was to experimentally verify the application of the continuum damage model, a progressive failure theory utilizing continuum damage mechanics, to a toughened material system. Damage due to tension-tension fatigue was documented for the IM7/5260 composite laminates. Crack density and delamination surface area were used to calculate matrix cracking and delamination internal state variables, respectively, to predict stiffness loss. A damage dependent finite element code qualitatively predicted trends in transverse matrix cracking, axial splits and local stress-strain distributions for notched quasi-isotropic laminates. The predictions were similar to the experimental data and it was concluded that the continuum damage model provided a good prediction of stiffness loss while qualitatively predicting damage growth in notched laminates.

Coats, Timothy William↗
NASA NTRS1994

Deflections of anisotropic sandwich beams with variable face sheets and core thicknesses

A sandwich construction consists of a low-density core material with high strength face sheets bounded to the top and bottom surfaces. The construction has been widely used in the aerospace and marine industries due to its outstanding characteristics such as noise absorption, weight minimization, heat insulation, and better bending stiffness. In sandwich structures used in high-performance aircraft, the face sheets are often made of fiber-reinforced composite materials and the core is made of honeycomb. The structures may also have variable thickness so as to satisfy aerodynamic requirements. In the stress analysis, the constant-thickness face sheets are usually considered as membrane and the core is assumed to be inextensible but deformable in the thickness direction. The static behavior of variable-thickness, isotropic and homogeneous sandwich beams was successfully studied by employing a constant-thickness theory but allowing stiffnesses to vary in accordance with local thickness variations. It has been recently found in a refined theory that the analyses based on the constant thickness theory locally can lead to significant errors in structural responses if the sandwich beam is thickness-tapered and the cores are deformable in transverse shear. The errors arise mainly from two factors: (1) the transverse shear components of the membrane forces in the face sheets alter the transverse shears carried by the core; and (2) the face-sheet membrane strains arise from transverse shear deformation of the core. In practice the variable thickness may not only exist in core but also in face sheets. The thickness-variations may even be a type of step function. In this case the transverse shear stress in the face sheets and bending stress in the core should be taken into account in the refined theory mentioned. In the present study, energy principles are employed in deriving governing equations for general bending of anisotropic sandwich beams with variable thickness in both face sheets and cores. Solutions to these equations are based on a finite difference scheme. As an example in application, a simply supported thickness-tapered sandwich beam subject to a concentrated load at its center is considered. Let W' be the maximum deflection of the beam in which face sheets are considered as membrane, while W'' is that based on using the modified refined theory. It is found that W' is always larger than W'', however, the magnitude of (W'- W'') appears to be insensitive to the change of the taper of the beam.

Lu, Chu-Ho↗
NASA NTRS1993

Towards a generalized computational fluid dynamics technique for all Mach numbers

Currently there exists no single unified approach for efficiently and accurately solving computational fluid dynamics (CFD) problems across the Mach number regime, from truly low speed incompressible flows to hypersonic speeds. There are several CFD codes that have evolved into sophisticated prediction tools with a wide variety of features including multiblock capabilities, generalized chemistry and thermodynamics models among other features. However, as these codes evolve, the demand placed on the end user also increases simply because of the myriad of features that are incorporated into these codes. In order for a user to be able to solve a wide range of problems, several codes may be needed requiring the user to be familiar with the intricacies of each code and their rather complicated input files. Moreover, the cost of training users and maintaining several codes becomes prohibitive. The objective of the current work is to extend the compressible, characteristic-based, thermochemical nonequilibrium Navier-Stokes code GASP to very low speed flows and simultaneously improve convergence at all speeds. Before this work began, the practical speed range of GASP was Mach numbers on the order of 0.1 and higher. In addition, a number of new techniques have been developed for more accurate physical and numerical modeling. The primary focus has been on the development of optimal preconditioning techniques for the Euler and the Navier-Stokes equations with general finite-rate chemistry models and both equilibrium and nonequilibrium thermodynamics models. We began with the work of Van Leer, Lee, and Roe for inviscid, one-dimensional perfect gases and extended their approach to include three-dimensional reacting flows. The basic steps required to accomplish this task were a transformation to stream-aligned coordinates, the formulation of the preconditioning matrix, incorporation into both explicit and implicit temporal integration schemes, and modification of the numerical flux formulae. In addition, we improved the convergence rate of the implicit time integration schemes in GASP through the use of inner iteration strategies and the use of the GMRES (General Minimized Resisual) which belongs to the class of algorithms referred to as Krylov subspace iteration. Finally, we significantly improved the practical utility of GASP through the addition of mesh sequencing, a technique in which computations begin on a coarse grid and get interpolated onto successively finer grids. The fluid dynamic problems of interest to the propulsion community involve complex flow physics spanning different velocity regimes and possibly involving chemical reactions. This class of problems results in widely disparate time scales causing numerical stiffness. Even in the absence of chemical reactions, eigenvalue stiffness manifests itself at transonic and very low speed flows which can be quantified by the large condition number of the system and evidenced by slow convergence rates. This results in the need for thorough numerical analysis and subsequent implementation of sophisticated numerical techniques for these difficult yet practical problems. As a result of this work, we have been able to extend the range of applicability of compressible codes to very low speed inviscid flows (M = .001) and reacting flows.

Walters, R. W.↗
NASA NTRS1995

Analysis and Development of Finite Element Methods for the Study of Nonlinear Thermomechanical Behavior of Structural Components

Underintegrated methods are investigated with respect to their stability and convergence properties. The focus was on identifying regions where they work and regions where techniques such as hourglass viscosity and hourglass control can be used. Results obtained show that underintegrated methods typically lead to finite element stiffness with spurious modes in the solution. However, problems exist (scalar elliptic boundary value problems) where underintegrated with hourglass control yield convergent solutions. Also, stress averaging in underintegrated stiffness calculations does not necessarily lead to stable or convergent stress states.

Oden, J. Tinsley↗
NASA NTRS1993

Test results from large wing and fuselage panels

This paper presents the first results in an assessment of the strength, stiffness, and damage tolerance of stiffened wing and fuselage subcomponents. Under this NASA funded program, 10 large wing and fuselage panels, variously fabricated by automated tow placement and dry-stitched preform/resin transfer molding, are to be tested. The first test of an automated tow placement six-longeron fuselage panel under shear load was completed successfully. Using NASTRAN finite-element analysis the stiffness of the panel in the linear range prior to buckling was predicted within 3.5 percent. A nonlinear analysis predicted the buckling load within 10 percent and final failure load within 6 percent. The first test of a resin transfer molding six-stringer wing panel under compression was also completed. The panel failed unexpectedly in buckling because of inadequate supporting structure. The average strain was 0.43 percent with a line load of 20.3 kips per inch of width. This strain still exceeds the design allowable strains. Also, the stringers did not debond before failure, which is in contrast to the general behavior of unstitched panels.

Madan, Ram C.↗