SEARCH · Engineering Papers
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.
Finite element method analysis of axisymmetrically loaded thin shells of revolution, considering nonlinearity due to material properties and shell geometry changes
A computational strategy is presented for the nonlinear static and postbuckling analyses of large complex structures on massively parallel computers. The strategy is designed for distributed-memory, message-passing parallel computer systems. The key elements of the proposed strategy are: (1) a multiple-parameter reduced basis technique; (2) a nested dissection (or multilevel substructuring) ordering scheme; (3) parallel assembly of global matrices; and (4) a parallel sparse equation solver. The effectiveness of the strategy is assessed by applying it to thermo-mechanical postbuckling analyses of stiffened composite panels with cutouts, and nonlinear large-deflection analyses of HSCT models on Intel Paragon XP/S computers. The numerical studies presented demonstrate the advantages of nested dissection-based solvers over traditional skyline-based solvers on distributed memory machines.
This ASTM Round Robin was conducted to evaluate the state of the art in stress analysis of adhesively bonded joint specimens. Specifically, the participants were asked to calculate the strain-energy-release rate for two different geometry cracked lap shear (CLS) specimens at four different debond lengths. The various analytical techniques consisted of 2- and 3-dimensional finite element analysis, beam theory, plate theory, and a combination of beam theory and finite element analysis. The results were examined in terms of the total strain-energy-release rate and the mode I to mode II ratio as a function of debond length for each specimen geometry. These results basically clustered into two groups: geometric linear or geometric nonlinear analysis. The geometric nonlinear analysis is required to properly analyze the CLS specimens. The 3-D finite element analysis gave indications of edge closure plus some mode III loading. Each participant described his analytical technique and results. Nine laboratories participated.
This ASTM Round Robin was conducted to evaluate the state of the art in stress analysis of adhesively bonded joint specimens. Specifically, the participants were asked to calculate the strain-energy-release rate for two different geometry cracked lap shear (CLS) specimens at four different debond lengths. The various analytical techniques consisted of 2- and 3-dimensional finite element analysis, beam theory, plate theory, and a combination of beam theory and finite element analysis. The results were examined in terms of the total strain-energy-release rate and the mode I to mode II ratio as a function of debond length for each specimen geometry. These results basically clustered into two groups: geometric linear or geometric nonlinear analysis. The geometric nonlinear analysis is required to properly analyze the CLS specimens. The 3-D finite element analysis gave indications of edge closure plus some mode III loading. Each participant described his analytical technique and results. Nine laboratories participated.
A nonlinear device modeling algorithm (NOLID) has been developed for use in spacecraft power system analysis. This algorithm is designed to explore the effects of nonlinear devices and loads on a spacecraft power system. Application of this harmonic modeling algorithm in spacecraft power system management programs such as harmonic power flow analysis packages is discussed. It is shown that the NOLID algorithm can be applied in conjunction with a harmonic power flow to give a more accurate description of system state.
Evaluating methods of numerically integrating equations of motion for nonlinear dynamic analyses of shells of revolution by matrix displacement method
A computational capability for predicting the nonlinear thermomechanical structural response of fiber-reinforced superalloy (FRS) turbine blades is described. This capability is embedded in a special purpose computer code (COBSTRAN) developed at the NASA Lewis Research Center. Special features of this computational capability include accounting for: fiber/matrix reaction, nonlinear and anisotropic material behavior, complex stress distribution due to local and global heterogeneity, and residual stresses due to initial fabrication and/or inelastic behavior during subsequent missions. Numerical results are presented from analyses of a hypothetical FRS turbine blade subjected to a fabrication process and subsequent mission cycle. The results demonstrate the capabilities of this computational tool to; predict local stress/strain response and capture trends of local nonlinear and anisotropic material behavior, relate the effects of this local behavior to the global response of a multilayered fiber-composite turbine blade, and trace material history from fabrication through successive missions.
In this report, we study the sensitivity of discrete-time dynamic programs with nonlinear dynamics and objective to perturbations in the initial conditions and reference parameters. Under uniform controllability and boundedness assumptions for the problem data, we prove that the directional derivative of the optimal state and control at time $k$, $x^*_k$, and $u^*_k$, with respect to the reference signal at time $i$, $d_i$, will have exponential decay in terms of $|k-i|$ with a decay rate $\rho$ independent of the temporal horizon length. The key technical step is to prove that a version of the convexification approach proposed by Verschueren et al. can be applied to the KKT conditions and results in a convex quadratic program with uniformly bounded data. In turn, Riccati techniques can be further employed to obtain the sensitivity result, borne from the observation that the directional derivatives are solutions of quadratic programs with structure similar to the KKT conditions themselves. We validate our findings with numerical experiments on a small nonlinear, nonconvex, dynamic program.
A finite element approach is presented for determining the nonlinear flutter characteristics of composite panels using unsteady, third-order piston theory aerodynamics. Both nonlinear structural (large-amplitude) and nonlinear aerodynamics terms are considered in the finite element formulation. Solution procedures are presented to solve the nonlinear panel flutter and the large-amplitude free vibration finite element equations. Nonlinear aerodynamic and linear structural finite element flutter results for composite panels are also presented. An efficient, vector-version generalized unsymmetric eigenequation solver is developed for large-amplitude vibration and nonlinear panel flutter analyses on high-performance computers.
A simplified inelastic analysis computer program was developed for predicting the stress-strain history of a thermomechanically cycled structure from an elastic solution. The program uses an iterative and incremental procedure to estimate the plastic strains from the material stress-strain properties and a simulated plasticity hardening model. The simplified method was exercised on a number of problems involving uniaxial and multiaxial loading, isothermal and nonisothermal conditions, and different materials and plasticity models. Good agreement was found between these analytical results and nonlinear finite element solutions for these problems. The simplified analysis program used less than 1 percent of the CPU time required for a nonlinear finite element analysis.
Results are presented of applying percolation analysis to several two-dimensional N-body models which simulate the formation of large-scale structure. Three parameters are estimated: total area (a(c)), total mass (M(C)), and percolation density (rho(c)) of the percolating structure at the percolation threshold for both unsmoothed and smoothed (with different scales L(s)) nonlinear with filamentary structures, confirming early speculations that this type of model has several features of filamentary-type distributions. Also, it is shown that, by properly applying smoothing techniques, many problems previously considered detrimental can be dealt with and overcome. Possible difficulties and prospects with the use of this method are discussed, specifically relating to techniques and methods already applied to CfA deep sky surveys. The success of this test in two dimensions and the potential for extrapolation to three dimensions is also discussed.
The application of finite element techniques to the analytic representation of the nonlinear behavior of arbitrary two-dimensional bodies containing cracks is discussed. Specific methods are proposed using which it should be possible to obtain information concerning: the description of the maximum, minimum, and residual near-tip stress and strain fields; the effects of crack closure on the near-tip behavior of stress and strain fields during cyclic loading into the plastic range; the stress-strain and displacement field behavior associated with a nonstationary crack; and the effects of large rotation near the crack tip.
Two new equivalent linearization implementations for geometrically nonlinear random vibrations are presented. Both implementations are based upon a novel approach for evaluating the nonlinear stiffness within commercial finite element codes and are suitable for use with any finite element code having geometrically nonlinear static analysis capabilities. The formulation includes a traditional force-error minimization approach and a relatively new version of a potential energy-error minimization approach, which has been generalized for multiple degree-of-freedom systems. Results for a simply supported plate under random acoustic excitation are presented and comparisons of the displacement root-mean-square values and power spectral densities are made with results from a nonlinear time domain numerical simulation.
A multiphase selfadaptive predictor corrector type algorithm was developed. This algorithm enables the solution of highly nonlinear structural responses including kinematic, kinetic and material effects as well as pro/post buckling behavior. The strategy involves three main phases: (1) the use of a warpable hyperelliptic constraint surface which serves to upperbound dependent iterate excursions during successive incremental Newton Ramphson (INR) type iterations; (20 uses an energy constraint to scale the generation of successive iterates so as to maintain the appropriate form of local convergence behavior; (3) the use of quality of convergence checks which enable various self adaptive modifications of the algorithmic structure when necessary. The restructuring is achieved by tightening various conditioning parameters as well as switch to different algorithmic levels to improve the convergence process. The capabilities of the procedure to handle various types of static nonlinear structural behavior are illustrated.
A problem-adaptive solution procedure for improving the reliability of finite element solutions to geometrically nonlinear shell-type problem is presented. The strategy incorporates automatic error detection and control and includes an iterative procedure which utilizes the solution at the same load step on a more refined model. Representative nonlinear shell problem are solved.
Automated fault management is at the forefront of next-generation optical communication networks. The increase in complexity of modern networks has triggered the need for programmable and software-driven architectures to support the operation of agile and self-managed systems. In these scenarios, the European Telecommunications Standards Institute zero-touch network and service management approach is imperative. The need for machine learning algorithms to process the large volume of telemetry data brings safety concerns as distributed cloud-computing solutions become the preferred approach for deploying reliable communication network automation. This paper’s contribution is twofold. First, we propose a simple yet effective method to guarantee the confidentiality of the telemetry data based on feature scrambling. The method allows the operation of third-party computational services without direct access to the full content of the collected data. Additionally, the effectiveness of four unsupervised machine learning algorithms for soft-failure detection is evaluated when applied to the scrambled telemetry data. The methods are based on factor analysis, principal component analysis, nonlinear principal component analysis, and singular value decomposition. Most dimensionality reduction algorithms have the common property that they can maintain similar levels of fault classification performance while hiding the data structure from unauthorized access. Evaluations of the proposed algorithms demonstrate this capability.
Preliminary analysis showed that the doubly curved triangular shell elements were too stiff for these shell structures. The doubly curved quadrilateral shell elements were found to give much improved results. A total of six load cases were analyzed in this study. The load cases were either those resulting from a static test using reaction straps to simulate the drop conditions or under assumed hydrodynamic conditions resulting from a drop test. The latter hydrodynamic conditions were obtained through an emperical fit of available data. Results obtained from a linear analysis were found to be consistent with results obtained elsewhere with NASTRAN and BOSOR. The nonlinear analysis showed that the originally assumed loads would result in failure of the shell structures. The nonlinear analysis also showed that it was useful to apply internal pressure as a stabilizing influence on collapse. A final analysis with an updated estimate of load conditions resulted in linear behavior up to full load.
The elastic instability of an internally-pressurized cylindrical tank with a torispherical head is investigated using a triangular, doubly curved, thin-shell finite element. The formulation of the finite element, which is based upon cubic displacement functions and a modified principle of potential energy, is first described. Then, the element is verified by comparing numerical results for the linear, stable analysis to alternative solutions for the same problem. The subsequent instability investigation includes the solution of the linearized problem of equilibrium bifurcation, that is, of the classical eigenvalue problem, and a general nonlinear analysis, based on tracing the nonlinear load-displacement path. The critical pressure, obtained with use of the general nonlinear analysis, agrees closely with an experimental result as well as with a numerical solution stemming from an axisymmetric formulation.