Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonlinear structural 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 487 records · Page 27

Multi-point Jig Twist Optimization of Mach 0.745 Transonic Truss-Braced Wing Aircraft and High-Fidelity CFD Validation

This paper presents a multi-point jig twist optimization study of the Mach 0.745 Transonic Truss-Braced Wing (TTBW) aircraft using an in-house developed aero-structural analysis solver VSPAERO coupled to BEAM3D. A vortex-lattice model of the TTBW model is developed, and a transonic small disturbance /integral boundary layer correction method is implemented in the VSPAERO model to account for transonic and viscous flow effects. A correction method for the wing-strut interference aerodynamics is developed and applied to the VSPAERO solver. Also, a structural dynamic finite-element model of the TTBW aircraft is developed. This finite-element model includes the geometric nonlinear effect due to the tension in the struts which causes a deflection-dependent nonlinear stiffness. The VSPAERO model is coupled to the corresponding finite-element model to provide a rapid aero-structural analysis. The VSPAERO model has been validated against wind tunnel test data and FUN3D CFD simulation data. Three different flight conditions corresponding to Mach 0.745 at 40,000 ft are selected for the TTBW aircraft jig twist optimization to reduce the drag coefficients at the design and off-design cruise lift coefficients. The preliminary optimization results show that the TTBW aircraft with the optimized jig twist achieves a drag reduction of about 9 counts, 7 counts, and8 counts corresponding to the lift coefficients 0.68, 0.73, and 0.78, respectively. A high-fidelity CFD solver FUN3D is used to validate the jig twist optimization.

Multi-point↗

Non-isothermal elastoviscoplastic analysis of planar curved beams

The development of a general mathematical model and solution methodologies, to examine the behavior of thin structural elements such as beams, rings, and arches, subjected to large nonisothermal elastoviscoplastic deformations is presented. Thus, geometric as well as material type nonlinearities of higher order are present in the analysis. For this purpose a complete true abinito rate theory of kinematics and kinetics for thin bodies, without any restriction on the magnitude of the transformation is presented. A previously formulated elasto-thermo-viscoplastic material constitutive law is employed in the analysis. The methodology is demonstrated through three different straight and curved beams problems.

Simitses, G. J.↗

Group implicit concurrent algorithms in nonlinear structural dynamics

During the 70's and 80's, considerable effort was devoted to developing efficient and reliable time stepping procedures for transient structural analysis. Mathematically, the equations governing this type of problems are generally stiff, i.e., they exhibit a wide spectrum in the linear range. The algorithms best suited to this type of applications are those which accurately integrate the low frequency content of the response without necessitating the resolution of the high frequency modes. This means that the algorithms must be unconditionally stable, which in turn rules out explicit integration. The most exciting possibility in the algorithms development area in recent years has been the advent of parallel computers with multiprocessing capabilities. So, this work is mainly concerned with the development of parallel algorithms in the area of structural dynamics. A primary objective is to devise unconditionally stable and accurate time stepping procedures which lend themselves to an efficient implementation in concurrent machines. Some features of the new computer architecture are summarized. A brief survey of current efforts in the area is presented. A new class of concurrent procedures, or Group Implicit algorithms is introduced and analyzed. The numerical simulation shows that GI algorithms hold considerable promise for application in coarse grain as well as medium grain parallel computers.

Ortiz, M.↗

Probability Bounds Analysis Applied to Multi-Purpose Crew Vehicle Nonlinearity

The Multi-Purpose Crew Vehicle (MPCV) Program Orion vehicle finite element model (FEM) was updated based on a modal test performed by Lockheed Martin. Due to nonlinearity observed in the test results, linear low force level (LL) and high force level (HL) FEMs were developed for use during various Space Launch System (SLS) flight regimes depending on expected forcing levels. Uncertainty models were derived for the combined MPCV and MPCV Stage Adaptor LL and HL Hurty/Craig-Bampton (HCB) components based on the MPCV structural test article Configuration 4 modal test-analysis correlation results. Subsequently, system-level uncertainty quantification analyses were performed using both models for various SLS flight configurations to determine the impact of the nonlinearity on important system metrics. The system metrics included both transfer functions associated with attitude control and dynamic loads associated with aerodynamic buffeting during ascent. In each case, an independent Monte Carlo (MC) analysis was performed, and no attempt was made to combine the results. The Hybrid Parametric Variation (HPV) method was used to develop the LL and HL MPCV HCB uncertainty models. The HPV method provides both parametric and non-parametric components of uncertainty. The non-parametric uncertainty accounts for the difference in model-form between the linearized analytical model and the corresponding linearized component test results in the form of mode shapes and frequencies at that force level. This linear model-form uncertainty is implemented in the HPV method using random matrix theory. However, the HPV uncertainty models developed for the linear LL and HL MPCV components do not account for the nonlinearity in the MPCV. With respect to the linearized models, this nonlinearity is also an uncertainty in model form, but in this case, it must be treated independently as an epistemic uncertainty. It represents a lack of knowledge, in contrast to an aleatory uncertainty due to the randomness of a variable. In the case of an epistemic variable, the true value is unknown, only the interval within which it lies is known. Epistemic uncertainty can be reduced with increased knowledge, while in general, aleatory uncertainty cannot. This work combines the epistemic uncertainty due to the MPCV nonlinearity with the parametric and non-parametric uncertainty within the HPV method using a second order propagation approach. The LL and HL test data is augmented with surrogate test data derived from a nonlinear MPCV representation. The impact of the MPCV nonlinearity on system response statistics is determined using a series of cumulative distribution functions in the form of a horsetail plot, or p-box. This results in an interval of probabilities for a specific response value, or an interval of response values at a specific probability.

Daniel C Kammer↗

Theoretical hypervelocity ballistic limit for single or double plates using nonlinear modal analysis

The original research on the use on nonlinear vibration technique to solve for the hypervelocity ballistic limit for double plates is examined. Such structure is commonly found in typical Space Station design where the incoming space or man-made debris would be fragmented upon hitting the outer plate (shield) and the subsequent impact on the main wall would result in a much reduced damge of the space station or spacecraft. The existing few theoretical impact equations do not agree well with each other. The existing computer code bumper used at NASA-Johnson Space Center appears to predict an unconservative ballistic limit when compared with experimental data where the velocity ranges from 3 km/s to 8 km/s. Such unconservative prediction is unacceptable from a practical safe design point of view. The bumper code is based on Wilkinson's (1968) paper and his equations have not been improved nor modified even though they are viewed with suspicion due to lack of agreement with experiments. To make matters worse, there is not another theory which is better than Wilkinson's equation and the designers are forced to use purely empirical Nysmith (1969) or semiempirical equations developed by Cour-Palais in 1969. The Cour-Palais equations were later modified empirically in 1989. The purpose of the present investigation is to examine the many assumptions of the Wilkinson equation. An attempt is made to present design charts based on the modified-Wilkinson equation so that the designer can get a feel of the ranges of the parameters which are of interest and discard a huge range of parameters, thus, significantly reducing the number of test shots required. The analysis is based on a solution of the governing nonlinear differential equations for a plate, assuming axisymmetric behavior using polar coordinates.

Hui, David↗

Engine dynamic analysis with general nonlinear finite element codes. Part 2: Bearing element implementation overall numerical characteristics and benchmaking

Finite element codes are used in modelling rotor-bearing-stator structure common to the turbine industry. Engine dynamic simulation is used by developing strategies which enable the use of available finite element codes. benchmarking the elements developed are benchmarked by incorporation into a general purpose code (ADINA); the numerical characteristics of finite element type rotor-bearing-stator simulations are evaluated through the use of various types of explicit/implicit numerical integration operators. Improving the overall numerical efficiency of the procedure is improved.

Padovan, J.↗

Mapping domain structures near a grain boundary in a lead zirconate titanate ferroelectric film using X-ray nanodiffraction

The effect of an electric field on local domain structure near a 24° tilt grain boundary in a 200 nm-thick Pb(Zr 0.2 Ti 0.8 )O 3 bi-crystal ferroelectric film was probed using synchrotron nanodiffraction. The bi-crystal film was grown epitaxially on SrRuO 3 -coated (001) SrTiO 3 24° tilt bi-crystal substrates. From the nanodiffraction data, real-space maps of the ferroelectric domain structure around the grain boundary prior to and during application of a 200 kV cm −1 electric field were reconstructed. In the vicinity of the tilt grain boundary, the distributions of densities of c -type tetragonal domains with the c axis aligned with the film normal were calculated on the basis of diffracted intensity ratios of c - and a -type domains and reference powder diffraction data. Diffracted intensity was averaged along the grain boundary, and it was shown that the density of c -type tetragonal domains dropped to ∼50% of that of the bulk of the film over a range ±150 nm from the grain boundary. This work complements previous results acquired by band excitation piezoresponse force microscopy, suggesting that reduced nonlinear piezoelectric response around grain boundaries may be related to the change in domain structure, as well as to the possibility of increased pinning of domain wall motion. The implications of the results and analysis in terms of understanding the role of grain boundaries in affecting the nonlinear piezoelectric and dielectric responses of ferroelectric materials are discussed.

36 MATERIALS SCIENCE↗

International Space Station 2A Array Modal Analysis

On December 9th 2009, the International Space Station (ISS) 2A solar array mast experienced prolonged longeron shadowing during a Soyuz undocking. Analytical reconstruction of induced thermal and dynamic structural loads showed an exceedance of the mast buckling limit. Possible structural damage to the solar array mast could have occurred during this event. A Low fidelity video survey of the 2A mast showed no obvious damage of the mast longerons or battens. The decision was made to conduct an on-orbit dynamic test of the 2A array on December 18th, 2009. The test included thruster pluming on the array while photogrammetry data was recorded. The test was similar to other Dedicated Thruster Firings (DTFs) that were performed to measure structural frequency and damping of a solar array. Results of the DTF indicated lower frequency mast modes than model predictions, thus leading to speculation of mast damage. A detailed nonlinear analysis was performed on the 2A array model to assess possible solutions to modal differences. The setup of the parametric nonlinear trade study included the use of a detailed array model and the reduced mass and stiffness matrices of the entire ISS being applied to the array interface. The study revealed that the array attachment structure is nonlinear and thus was the source of error in the model prediction of mast modes. In addition, a detailed study was performed to determine mast mode sensitivity to mast longeron damage. This sensitivity study was performed to assess if the ISS program has sufficient instrumentation for mast damage detection.

Laible, Michael↗

Computational methods for structural load and resistance modeling

An automated capability for computing structural reliability considering uncertainties in both load and resistance variables is presented. The computations are carried out using an automated Advanced Mean Value iteration algorithm (AMV +) with performance functions involving load and resistance variables obtained by both explicit and implicit methods. A complete description of the procedures used is given as well as several illustrative examples, verified by Monte Carlo Analysis. In particular, the computational methods described in the paper are shown to be quite accurate and efficient for a material nonlinear structure considering material damage as a function of several primitive random variables. The results show clearly the effectiveness of the algorithms for computing the reliability of large-scale structural systems with a maximum number of resolutions.

Thacker, B. H.↗

Aeroelastic analysis for helicopter rotor blades with time-variable, non-linear structural twist and multiple structural redundancy: Mathematical derivation and program user's manual

The differential equations of motion for the lateral and torsional deformations of a nonlinearly twisted rotor blade in steady flight conditions together with those additional aeroelastic features germane to composite bearingless rotors are derived. The differential equations are formulated in terms of uncoupled (zero pitch and twist) vibratory modes with exact coupling effects due to finite, time variable blade pitch and, to second order, twist. Also presented are derivations of the fully coupled inertia and aerodynamic load distributions, automatic pitch change coupling effects, structural redundancy characteristics of the composite bearingless rotor flexbeam - torque tube system in bending and torsion, and a description of the linearized equations appropriate for eigensolution analyses. Three appendixes are included presenting material appropriate to the digital computer program implementation of the analysis, program G400.

Bielawa, R. L.↗

Effects of nonlinear aerodynamics and static aeroelasticity on mission performance calculations for a fighter aircraft

During conceptual design studies of advanced aircraft, the usual practice is to use linear theory to calculate the aerodynamic characteristics of candidate rigid (nonflexible) geometric external shapes. Recent developments and improvements in computational methods, especially computational fluid dynamics (CFD), provide significantly improved capability to generate detailed analysis data for the use of all disciplines involved in the evaluation of a proposed aircraft design. A multidisciplinary application of such analysis methods to calculate the effects of nonlinear aerodynamics and static aeroelasticity on the mission performance of a fighter aircraft concept is described. The aircraft configuration selected for study was defined in a previous study using linear aerodynamics and rigid geometry. The results from the previous study are used as a basis of comparison for the data generated herein. Aerodynamic characteristics are calculated using two different nonlinear theories, potential flow and rotational (Euler) flow. The aerodynamic calculations are performed in an iterative procedure with an equivalent plate structural analysis method to obtain lift and drag data for a flexible (nonrigid) aircraft. These static aeroelastic data are then used in calculating the combat and mission performance characteristics of the aircraft.

Giles, Gary L.↗

TAP 1: A Finite Element Program for Steady-State Thermal Analysis of Convectively Cooled Structures

The program has a finite element library of six elements: two conduction/convection elements to model heat transfer in a solid, two convection elements to model heat transfer in a fluid, and two integrated conduction/convection elements to represent combined heat transfer in tubular and plate/fin fluid passages. Nonlinear thermal analysis due to temperature dependent thermal parameters is performed using the Newton-Raphson iteration method. Program output includes nodal temperatures and element heat fluxes. Pressure drops in fluid passages may be computed as an option. A companion plotting program for displaying the finite element model and predicted temperature distributions is presented. User instructions and sample problems are presented in appendixes.

Thornton, E. A.↗

Component mode synthesis and large deflection vibrations of complex structures

The accuracy of the NASTRAN modal synthesis analysis was assessed by comparing it with full structure NASTRAN and nine other modal synthesis results using a nine-bay truss. A NASTRAN component mode transient response analysis was also performed on the free-free truss structure. A finite element method was developed for nonlinear vibration of beam structures subjected to harmonic excitation. Longitudinal deformation and inertia are both included in the formula. Tables show the finite element free vibration results with and without considering the effects of longitudinal deformation and inertia as well as the frequency ratios for a simply supported and a clamped beam subjected to a uniform harmonic force.

Mei, C.↗

Impact of land-surface moisture variability on local shallow convective cumulus and precipitation in large-scale models

Numerical experiments using state-of-the-art high-resolution mesoscale cloud model showed that land-surface moisture significantly affects the timing of onset of clouds and the intensity and distribution of precipitation. In general, landscape discontinuity enhances shallow convective precipitation. Two mechanisms that are strongly modulated by land-surface moisture-namely, random turbulent thermal cells and organized sea-breeze-like mesoscale circulations-also determine the horizontal distribution of maximum precipitation. However, interactions between shallow cumulus and land-surface moisture are highly nonlinear and complicated by different factors, such as atmospheric thermodynamic structure and large-scale background wind. This analysis also showed that land-surface moisture discontinuities seem to play a more important role in a relatively dry atmsophere, and that the strongest precipitation is produced by a wavelength of land-surface forcing equivalent to the local Rossby radius of deformation. A general trend between the maximum precipitation and the normalized maximum latent heat flux was identified. In general, large values of mesoscale latent heat flux imply strongly developed mesoscale circulations and intense cloud activity, accompanied by large surface latent heat fluxes that transport more water vapor into the atmosphere.

Chen, Fei↗

Scalability of Cohesive Fatigue Analyses Using Explicit Solvers

A cohesive fatigue law has been integrated into a constitutive material model compatible with an explicit finite element solver. The cohesive fatigue model response is based on engineering approximations of the endurance limit and the Goodman diagram. This approach can predict stress-life diagrams for crack initiation, the Paris law regime, and transient effects of crack initiation and stable tearing. Simplified cyclic loading is utilized so that the applied load(or displacement) corresponds to the peak load of a fatigue cycle. Loads are held constant during fatigue while damage develops with increasing solution increments. An automatically-calculated ratio of fatigue cycles per solution increment controls the rate of damage growth, ensuring that damage growth is modeled with a sufficient minimum number of increments and damage growth advances to a minimum desired extent within the explicit analysis step time. The compatibility with an explicit finite element solver enables the analysis of structures that are computationally intractable for implicit finite element solvers. Scalability studies are conducted for geometrically nonlinear problems involving fiber-reinforced composite structures that exhibit fatigue damage growth of interacting matrix cracks and delaminations

Frank A Leone↗

Effects of joints in truss structures

The response of truss-type structures for future space applications, such as Large Deployable Reflector (LDR), will be directly affected by joint performance. Some of the objectives of research at BAC were to characterize structural joints, establish analytical approaches that incorporate joint characteristics, and experimentally establish the validity of the analytical approaches. The test approach to characterize joints for both erectable and deployable-type structures was based upon a Force State Mapping Technique. The approach pictorially shows how the nonlinear joint results can be used for equivalent linear analysis. Testing of the Space Station joints developed at LaRC (a hinged joint at 2 Hz and a clevis joint at 2 Hz) successfully revealed the nonlinear characteristics of the joints. The Space Station joints were effectively linear when loaded to plus or minus 500 pounds with a corresponding displacement of about plus or minus 0.0015 inch. It was indicated that good linear joints exist which are compatible with errected structures, but that difficulty may be encountered if nonlinear-type joints are incorporated in the structure.

Ikegami, R.↗

Nonlinear probabilistic finite element models of laminated composite shells

A probabilistic finite element analysis procedure for laminated composite shells has been developed. A total Lagrangian finite element formulation, employing a degenerated 3-D laminated composite shell with the full Green-Lagrange strains and first-order shear deformable kinematics, forms the modeling foundation. The first-order second-moment technique for probabilistic finite element analysis of random fields is employed and results are presented in the form of mean and variance of the structural response. The effects of material nonlinearity are included through the use of a rate-independent anisotropic plasticity formulation with the macroscopic point of view. Both ply-level and micromechanics-level random variables can be selected, the latter by means of the Aboudi micromechanics model. A number of sample problems are solved to verify the accuracy of the procedures developed and to quantify the variability of certain material type/structure combinations. Experimental data is compared in many cases, and the Monte Carlo simulation method is used to check the probabilistic results. In general, the procedure is quite effective in modeling the mean and variance response of the linear and nonlinear behavior of laminated composite shells.

Engelstad, S. P.↗

High temperature stress-strain analysis

The objectives of the high-temperature structures program are threefold: to assist in the development of analytical tools needed to improve design analyses and procedures for the efficient and accurate prediction of the nonlinear structural response of hot-section components; to aid in the calibration, validation, and evaluation of the analytical tools by comparing predictions with experimental data; and to evaluate existing as well as advanced temperature and strain measurement instrumentation. As the analytical tools, test methods, tests, instrumentations, as well as data acquisition, management, and analysis methods are developed and evaluated, a proven, integrated analysis and experiment method will result in a more accurate prediction of the cyclic life of hot section components.

Thompson, Robert L.↗