Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear analysis”

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

Critical joints in large composite aircraft structure

A program was conducted at Douglas Aircraft Company to develop the technology for critical structural joints of composite wing structure that meets design requirements for a 1990 commercial transport aircraft. The prime objective of the program was to demonstrate the ability to reliably predict the strength of large bolted composite joints. Ancillary testing of 180 specimens generated data on strength and load-deflection characteristics which provided input to the joint analysis. Load-sharing between fasteners in multirow bolted joints was computed by the nonlinear analysis program A4EJ. This program was used to predict strengths of 20 additional large subcomponents representing strips from a wing root chordwise splice. In most cases, the predictions were accurate to within a few percent of the test results. In some cases, the observed mode of failure was different than anticipated. The highlight of the subcomponent testing was the consistent ability to achieve gross-section failure strains close to 0.005. That represents a considerable improvement over the state of the art.

Nelson, W. D.↗

Statistical foundation of the fluid analogue of the soliton formalism

A fully nonlinear analysis is used to develop a general soliton formalism for the description of the nonlinear evolution of soliton fluctuations in both plasmas and classical fluids. From the Navier-Stokes equations for plasmas and compressible fluids of two scales, two equations for the propagation of density waves are derived. A fast soliton field is spontaneously created by rarefaction, and a slow density wave modulates the field intensity as a ponderomotive force. Constitutive properties are demonstrated using a Lagrangian-kinetic formalism of the fluctuation-dissipation theory.

Tchen, C. M.↗

Efficient Computation Of Behavior Of Aircraft Tires

NASA technical paper discusses challenging application of computational structural mechanics to numerical simulation of responses of aircraft tires during taxing, takeoff, and landing. Presents details of three main elements of computational strategy: use of special three-field, mixed-finite-element models; use of operator splitting; and application of technique reducing substantially number of degrees of freedom. Proposed computational strategy applied to two quasi-symmetric problems: linear analysis of anisotropic tires through use of two-dimensional-shell finite elements and nonlinear analysis of orthotropic tires subjected to unsymmetric loading. Three basic types of symmetry and combinations exhibited by response of tire identified.

Tanner, John A.↗

The development of a mixing layer under the action of weak streamwise vortices

The action of weak, streamwise vortices on a plane, incompressible, steady mixing layer is examined in the large Reynolds number limit. The outer, inviscid region is bounded by a vortex sheet to which the viscous region is confined. It is shown that the local linear analysis becomes invalid at streamwise distances O(epsilon sup -1), where (epsilon much less than 1) is the crossflow amplitude, and a new nonlinear analysis is constructed for this region. Numerical solutions of the nonlinear problem show that the vortex sheet undergoes an O(1) change in position and that the solution is ultimately terminated by a breakdown in the numerical procedure. The corresponding viscous layer shows downstream thickening, but appears to remain well behaved up to the terminal location.

Goldstein, M. E.↗

Modeling and analysis of the space shuttle nose-gear tire with semianalytic finite elements

A computational procedure is presented for the geometrically nonlinear analysis of aircraft tires. The Space Shuttle Orbiter nose gear tire was modeled by using a two-dimensional laminated anisotropic shell theory with the effects of variation in material and geometric parameters included. The four key elements of the procedure are: (1) semianalytic finite elements in which the shell variables are represented by Fourier series in the circumferential direction and piecewise polynominals in the meridional direction; (2) a mixed formulation with the fundamental unknowns consisting of strain parameters, stress-resultant parameters, and generalized displacements; (3) multilevel operator splitting to effect successive simplifications, and to uncouple the equations associated with different Fourier harmonics; and (4) multilevel iterative procedures and reduction techniques to generate the response of the shell. Numerical results of the Space Shuttle Orbiter nose gear tire model are compared with experimental measurements of the tire subjected to inflation loading.

Kim, Kyun O.↗

Optimization of Car Body under Constraints of Noise, Vibration, and Harshness (NVH), and Crash

To be competitive on the today's market, cars have to be as light as possible while meeting the Noise, Vibration, and Harshness (NVH) requirements and conforming to Government-man dated crash survival regulations. The latter are difficult to meet because they involve very compute-intensive, nonlinear analysis, e.g., the code RADIOSS capable of simulation of the dynamics, and the geometrical and material nonlinearities of a thin-walled car structure in crash, would require over 12 days of elapsed time for a single design of a 390K elastic degrees of freedom model, if executed on a single processor of the state-of-the-art SGI Origin2000 computer. Of course, in optimization that crash analysis would have to be invoked many times. Needless to say, that has rendered such optimization intractable until now. The car finite element model is shown. The advent of computers that comprise large numbers of concurrently operating processors has created a new environment wherein the above optimization, and other engineering problems heretofore regarded as intractable may be solved. The procedure, shown, is a piecewise approximation based method and involves using a sensitivity based Taylor series approximation model for NVH and a polynomial response surface model for Crash. In that method the NVH constraints are evaluated using a finite element code (MSC/NASTRAN) that yields the constraint values and their derivatives with respect to design variables. The crash constraints are evaluated using the explicit code RADIOSS on the Origin 2000 operating on 256 processors simultaneously to generate data for a polynomial response surface in the design variable domain. The NVH constraints and their derivatives combined with the response surface for the crash constraints form an approximation to the system analysis (surrogate analysis) that enables a cycle of multidisciplinary optimization within move limits. In the inner loop, the NVH sensitivities are recomputed to update the NVH approximation model while keeping the Crash response surface constant. In every outer loop, the Crash response surface approximation is updated, including a gradual increase in the order of the response surface and the response surface extension in the direction of the search. In this optimization task, the NVH discipline has 30 design variables while the crash discipline has 20 design variables. A subset of these design variables (10) are common to both the NVH and crash disciplines. In order to construct a linear response surface for the Crash discipline constraints, a minimum of 21 design points would have to be analyzed using the RADIOSS code. On a single processor in Origin 2000 that amount of computing would require over 9 months! In this work, these runs were carried out concurrently on the Origin 2000 using multiple processors, ranging from 8 to 16, for each crash (RADIOSS) analysis. Another figure shows the wall time required for a single RADIOSS analysis using varying number of processors, as well as provides a comparison of 2 different common data placement procedures within the allotted memories for each analysis. The initial design is an infeasible design with NVH discipline Static Torsion constraint violations of over 10%. The final optimized design is a feasible design with a weight reduction of 15 kg compared to the initial design. This work demonstrates how advanced methodology for optimization combined with the technology of concurrent processing enables applications that until now were out of reach because of very long time-to-solution.

Kodiyalam, Srinivas↗

Vector algorithms for geometrically nonlinear 3D finite element analysis

Algorithms for geometrically nonlinear finite element analysis are presented which exploit the vector processing capability of the VPS-32, which is closely related to the CYBER 205. By manipulating vectors (which are long lists of numbers) rather than individual numbers, very high processing speeds are obtained. Long vector lengths are obtained without extensive replication or reordering by storage of intermediate results in strategic patterns at all stages of the computations. Comparisons of execution times with those from programs using either scalar or other vector programming techniques indicate that the algorithms presented are quite efficient.

Whitcomb, John D.↗

Failure analysis of thick composite cylinders under external pressure

Failure of thick section composites due to local compression strength and overall structural instability is treated. Effects of material nonlinearity, imperfect fiber architecture, and structural imperfections upon anticipated failure stresses are determined. Comparisons with experimental data for a series of test cylinders are described. Predicting the failure strength of composite structures requires consideration of stability and material strength modes of failure using linear and nonlinear analysis techniques. Material strength prediction requires the accurate definition of the local multiaxial stress state in the material. An elasticity solution for the linear static analysis of thick anisotropic cylinders and rings is used herein to predict the axisymmetric stress state in the cylinders. Asymmetric nonlinear behavior due to initial cylinder out of roundness and the effects of end closure structure are treated using finite element methods. It is assumed that local fiber or ply waviness is an important factor in the initiation of material failure. An analytical model for the prediction of compression failure of fiber composites, which includes the effects of fiber misalignments, matrix inelasticity, and multiaxial applied stresses is used for material strength calculations. Analytical results are compared to experimental data for a series of glass and carbon fiber reinforced epoxy cylinders subjected to external pressure. Recommendations for pretest characterization and other experimental issues are presented. Implications for material and structural design are discussed.

Caiazzo, A.↗

The development of a mixing layer under the action of weak streamwise vortices

The action of weak, streamwise vortices on a plane, incompressible, steady mixing layer is examined in the large Reynolds-number limit. The outer, inviscid region is bounded by a vortex sheet to which the viscous region is confined. It is shown that the local linear analysis becomes invalid at streamwise distances O(epsilon(sup -1)), where epsilon is much less than 1 is the cross flow amplitude, and a new nonlinear analysis is constructed for this region. Numerical solutions of the nonlinear problem show that the vortex sheet undergoes an O(1) change in position and that the solution is ultimately terminated by the appearance of a singularity. The corresponding viscous layer shows downstream thickening, but appears to remain well behaved up to the singular location.

Goldstein, Marvin E.↗

Computer codes for thermal analysis of a solid rocket motor nozzle

A number of computer codes are available for performing thermal analysis of solid rocket motor nozzles. Aerotherm Chemical Equilibrium (ACE) computer program can be used to perform one-dimensional gas expansion to determine the state of the gas at each location of a nozzle. The ACE outputs can be used as input to a computer program called Momentum/Energy Integral Technique (MEIT) for predicting boundary layer development development, shear, and heating on the surface of the nozzle. The output from MEIT can be used as input to another computer program called Aerotherm Charring Material Thermal Response and Ablation Program (CMA). This program is used to calculate oblation or decomposition response of the nozzle material. A code called Failure Analysis Nonlinear Thermal and Structural Integrated Code (FANTASTIC) is also likely to be used for performing thermal analysis of solid rocket motor nozzles after the program is duly verified. A part of the verification work on FANTASTIC was done by using one and two dimension heat transfer examples with known answers. An attempt was made to prepare input for performing thermal analysis of the CCT nozzle using the FANTASTIC computer code. The CCT nozzle problem will first be solved by using ACE, MEIT, and CMA. The same problem will then be solved using FANTASTIC. These results will then be compared for verification of FANTASTIC.

Chauhan, Rajinder Singh↗

A method for nonlinear exponential regression analysis

A computer-oriented technique is presented for performing a nonlinear exponential regression analysis on decay-type experimental data. The technique involves the least squares procedure wherein the nonlinear problem is linearized by expansion in a Taylor series. A linear curve fitting procedure for determining the initial nominal estimates for the unknown exponential model parameters is included as an integral part of the technique. A correction matrix was derived and then applied to the nominal estimate to produce an improved set of model parameters. The solution cycle is repeated until some predetermined criterion is satisfied.

Junkin, B. G.↗

Role of cross-field current instability in substorm onsets and intensifications

A cross field current instability is investigated as a potential mechanism for current reduction/disruption during substorm onsets and intensifications. Linear stability analysis shows that, for sufficiently strong current density, the instability can occur in the plasma sheet with a growth time comparable to the substorm onset time and excite waves with a significant electromagnetic component. Nonlinear analysis shows that the wave growth reaches the nonlinear stage in less than one ion gyroperiod and can reduce the cross tail current by approximately 15 to 28 of its initial value at saturation. The resulting anomalous resistivity is 11 to 12 orders of magnitude above the classical value. For a typical current reduction, the plasma in the disruption region is subjected to an earthward force. The substorm development scenario constructed based on this instability can readily account for a large number of substorm features.

Lui, Anthony T. Y.↗

Three dimensional inelastic finite element analysis of laminated composites

Formulations of the inelastic response of laminated composites to thermal and mechanical loading are used as the basis for development of the computer NALCOM (Nonlinear Analysis of Laminated Composites) computer program which uses a fully three dimensional isoparametric finite element with 24 nodes and 72 degrees of freedom. An incremental solution is performed with nonlinearities introduced as pseudoloads computed for initial strains. Equilibrium iteration may be performed at every step. Elastic and elastic-plastic response of boron/epoxy and graphite/epoxy graphite/epoxy and problems of curing 0/90 sub s Gr/Ep laminates with and without circular holes are analyzed. Mechanical loading of + or - 45sub s Gr/Ep laminates is modeled and symmetry conditions which exist in angle-ply laminates are discussed. Results are compared to experiments and other analytical models when possible. All models are seen to agree reasonably well with experimetnal results for off-axis tensile coupons. The laminate analyses show the three dimensional effects which are present near holes and free corners.

Griffin, O. H., Jr.↗

Application of the reduced basis method to nonlinear transient thermal analysis

An effort to apply the reduced basis method to nonlinear transient thermal analysis is described. The method combines the classical Rayleigh-Ritz and modal superposition techniques with contemporary finite element methods to retain modeling versatility as the degrees of freedom in a problem are reduced. The essence of the method is to use a few thermal modes from eigenvalue analyses as basis vectors to represent the temperature response for a given thermal problem similar to the use of vibration modes to represent displacements in a dynamic response problem. Approximate temperature distributions were obtained using the reduced basis method for a small section of the Shuttle Orbiter lower wing undergoing reentry heating. Good agreement was obtained between the reduced basis method solutions and full system solutions with reductions in the degrees of freedom of up to a factor of four. The good agreement indicates the reduced basis method has the potential for significant reduction in computing effort for thermal analysis; however, considerable work remains to determine techniques for selecting the type and number of basis vectors needed for approximate solutions to more complex transient thermal problems.

Shore, C. P.↗

A new energy transfer model for turbulent free shear flow

A new model for the energy transfer mechanism in the large-scale turbulent kinetic energy equation is proposed. An estimate of the characteristic length scale of the energy containing large structures is obtained from the wavelength associated with the structures predicted by a weakly nonlinear analysis for turbulent free shear flows. With the inclusion of the proposed energy transfer model, the weakly nonlinear wave models for the turbulent large-scale structures are self-contained and are likely to be independent flow geometries. The model is tested against a plane mixing layer. Reasonably good agreement is achieved. Finally, it is shown by using the Liapunov function method, the balance between the production and the drainage of the kinetic energy of the turbulent large-scale structures is asymptotically stable as their amplitude saturates. The saturation of the wave amplitude provides an alternative indicator for flow self-similarity.

Liou, William W.-W.↗

A Numerical Simulator for Three-Dimensional Flows Through Vibrating Blade Rows

The three-dimensional, multi-stage, unsteady, turbomachinery analysis, TURBO, has been extended to predict the aeroelastic and aeroacoustic response behaviors of a single blade row operating within a cylindrical annular duct. In particular, a blade vibration capability has been incorporated so that the TURBO analysis can be applied over a solution domain that deforms with a vibratory blade motion. Also, unsteady far-field conditions have been implemented to render the computational boundaries at inlet and exit transparent to outgoing unsteady disturbances. The modified TURBO analysis is applied herein to predict unsteady subsonic and transonic flows. The intent is to partially validate this nonlinear analysis for blade flutter applications, via numerical results for benchmark unsteady flows, and to demonstrate the analysis for a realistic fan rotor. For these purposes, we have considered unsteady subsonic flows through a 3D version of the 10th Standard Cascade, and unsteady transonic flows through the first stage rotor of the NASA Lewis, Rotor 67, two-stage fan.

Chuang, H. Andrew↗

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.↗

Quantifying Wrinkle Features of Thin Membrane Structures

For future micro-systems utilizing membrane based structures, quantified predictions of wrinkling behavior in terms of amplitude, angle and wavelength are needed to optimize the efficiency and integrity of such structures, as well as their associated control systems. For numerical analyses performed in the past, limitations on the accuracy of membrane distortion simulations have often been related to the assumptions made. This work demonstrates that critical assumptions include: effects of gravity, supposed initial or boundary conditions, and the type of element used to model the membrane. In this work, a 0.2 m x 02 m membrane is treated as a structural material with non-negligible bending stiffness. Finite element modeling is used to simulate wrinkling behavior due to a constant applied in-plane shear load. Membrane thickness, gravity effects, and initial imperfections with respect to flatness were varied in numerous nonlinear analysis cases. Significant findings include notable variations in wrinkle modes for thickness in the range of 50 microns to 1000 microns, which also depend on the presence of an applied gravity field. However, it is revealed that relationships between overall strain energy density and thickness for cases with differing initial conditions are independent of assumed initial conditions. In addition, analysis results indicate that the relationship between wrinkle amplitude scale (W/t) and structural scale (L/t) is independent of the nonlinear relationship between thickness and stiffness.

Jacobson, Mindy B.↗