Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear problems”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Stability, accuracy, and efficiency of some underintegrated methods in finite element computations

In an attempt to increase computational efficiency in the numerical solution of highly nonlinear problems in solid and fluid mechanics, underintegrated finite element methods have been employed by many analysts. Underintegration refers to the use of a rule of an order lower than that required to integrate polynomial integrands exactly. The main drawback of this technique is related to the production of rank-deficient stiffness matrices, or equivalently an expanded kernel of the governing linear momentum operators. Such a development can introduce numerical instabilities. In order to overcome this difficulty, artificial stiffness or viscosity methods, or other stabilization methods have been proposed. One approach involves the elimination of spurious modes in a postprocessing operation. The present study is concerned with this a posteriori elimination method, taking into account the results which can be expected from it, and some of its possible extensions.

Jacquotte, O.-P.↗
NASA NTRS1985

Crystal Melts Under Reduced Gravity

To date, fluid motion due to thermocapillary convection and centrifugal buoyancy has been studied theoretically and numerically, by adapting specialized programs for rotating viscous fluids. It has been shown that by rotating the crystal melt, thermocapillary flow can be confined to a thin layer at the melt-gas interface even when the interfacial velocity is large enough for non-linear effects to be important. Theoretical analysis of the basic processes proceeds, but a much larger effort is now being devoted to the numerical simulation of the actual motion in a rotating crystal melt. A numerical program is under development for the study of the fully nonlinear problem. Although the flow is viewed, to start, as one of rapid rotation (in order that certain boundary layer simplifications can be invoked), a more general approach will allow consideration of the situations and conditions of real practical interest.

Greenspan, H. P.↗
NASA NTRS1985

Analysis of a finite composite plate with smooth rigid pin

An elastic continuum method of analysis developed for a smooth rigid pin in a finite composite plate by using complex potentials that satisfy the governing orthotropic differential equation is presented. The plate is subjected to uniaxial loading and the nonlinear problems of partial contact in the pin-plate are studied. Three types of pin fits are examined: interference, push, and clearance. The problems of stresses due to pure interference, load contact behavior, and stresses variations are addressed. The effects of the rectangular geometry and the orthotropic property of the plate are demonstrated through numerical results. It is found that the use of large hole diameter/plate width ratios have deleterious effects and the use of interference has beneficial effects. The solutions obtained will help in the development of special finite elements for fastener joints. Finally, the boundary conditions on the finite outer boundaries are found to be satisfied by a successive integration technique.

Naidu, A. C. B.↗
NASA NTRS1986

Nonlinear response - A time domain approach

The present paper reviews the basic concepts of nonlinear response of panels to surface flow and acoustic pressures, simulation of random processes, time domain solutions and the Monte Carlo Method. Applications of this procedure to the orbit-on-demand space vehicles, acoustic fatigue and composite materials are discussed. Numerical examples are included for a variety of nonlinear problems to illustrate the applicability of this method.

Vaicaitis, R.↗
NASA NTRS1987

A creep-rupture model of filament-wound spherical pressure vessels

The creep-rupture model is that of a quasi-isotropic filament-wound spherical pressure vessel, subjected to internal pressurization. The matrix material of the composite system is assumed to be linearly viscoelastic. Internal damage resulting from the relaxation of the matrix and the corresponding increase in microcracks is represented by a functional relationship between circumferential strain and transverse modulus. The numerical solution to this nonlinear problem is an iterative technique, whereby the elastic-viscoelastic correspondence principle is employed. In the Laplace domain, the associated elastic solution is obtained and this solution is inverted by the multidata method to yield the time-dependent solution.

Dozier, Jan D.↗
NASA NTRS1987

Acoustic fatigue - A Monte Carlo approach

The present paper reviews the basic concepts of nonlinear structural response to surface flow and acoustic pressure inputs. A time domain solution and the Monte Carlo method are utilized to estimate the required statistics of the nonlinear response. Simulation procedures of multi-dimensional and multi-variate random processes are discussed. Application of these procedures to acoustic fatigue of flight structures are considered. Numerical examples are included to illustrate the applicability of this method for the solution of nonlinear problems.

Vaicaitis, R.↗
NASA NTRS1987

Viscoelastic/damage modeling of filament-wound spherical pressure vessels

A model of the viscoelastic/damage response of a filament-wound spherical vessel used for long-term pressure containment is developed. The matrix material of the composite system is assumed to be linearly viscoelastic. Internal accumulated damage based upon a quadratic relationship between transverse modulus and maximum circumferential strain is postulated. The resulting nonlinear problem is solved by an iterative routine. The elastic-viscoelastic correspondence is employed to produce, in the Laplace domain, the associated elastic solution for the maximum circumferential strain which is inverted by the method of collocation to yield the time-dependent solution. Results obtained with the model are compared to experimental observations.

Hackett, Robert M.↗
NASA NTRS1988

An approach to CMG steering using feedback linearization

This paper presents an approach for controlling spacecraft equipped with control moment gyroscopes. A technique from feedback linearization theory is used to transform the original nonlinear problem to an equivalent linear form without approximating assumptions. In this form, the spacecraft dynamics appear linearly, and are decoupled from redundancy in the system of gyroscopes. A general approach to distributing control effort among the available actuators is described which includes provisions for redistribution of rotors, explicit bounds in gimbal rates, and guaranteed operation at or near singular configurations. A particular algorithm is developed for systems of double-gimbal devices, and demonstrated in two examples for which existing approaches fail to give adequate performance.

Dzielski, John↗
NASA NTRS1987

Finite element methods in probabilistic mechanics

Probabilistic methods, synthesizing the power of finite element methods with second-order perturbation techniques, are formulated for linear and nonlinear problems. Random material, geometric properties and loads can be incorporated in these methods, in terms of their fundamental statistics. By construction, these methods are applicable when the scale of randomness is not too large and when the probabilistic density functions have decaying tails. By incorporating certain computational techniques, these methods are shown to be capable of handling large systems with many sources of uncertainties. Applications showing the effects of combined random fields and cyclic loading/stress reversal are studied and compared with Monte Carlo simulation results.

Liu, Wing Kam↗
NASA NTRS1989

Domain decomposition methods for systems of conservation laws: Spectral collocation approximations

Hyperbolic systems of conversation laws are considered which are discretized in space by spectral collocation methods and advanced in time by finite difference schemes. At any time-level a domain deposition method based on an iteration by subdomain procedure was introduced yielding at each step a sequence of independent subproblems (one for each subdomain) that can be solved simultaneously. The method is set for a general nonlinear problem in several space variables. The convergence analysis, however, is carried out only for a linear one-dimensional system with continuous solutions. A precise form of the error reduction factor at each iteration is derived. Although the method is applied here to the case of spectral collocation approximation only, the idea is fairly general and can be used in a different context as well. For instance, its application to space discretization by finite differences is straight forward.

Quarteroni, Alfio↗
NASA NTRS1989

Unconditionally stable concurrent procedures for transient finite-element analysis

A family of algorithms was outlined which would appear to be particularly well-suited for implementation in a parallel environment. This is due to the fact that for any partition of the mesh each subdomain in the partition can be processed over a time step simultaneously and independently of the rest. The method eliminates the need for assembling and factorizing large global arrays while retaining the unconditional stability properties of the algorithms used at the local level. To critically appraise the proposed methodology, two limiting cases were considered: element-by-element mesh partitions, and coarse mesh partitions. It was concluded that while the proposed methodology can be useful in sequential machines, it would appear to be promising as it bears on computation. It should also be emphasized that extensions of the method to nonlinear problems are possible.

Ortiz, Michael↗
NASA NTRS1990

Concurrent Algorithm For Particle-In-Cell Simulations

Separate decompositions used for particle-motion and field calculations. General Concurrent Particle-in-Cell (GCPIC) algorithm used to implement motions of individual plasma particles (ions and electrons) under influence of particle-in-cell (PIC) computer codes on concurrent processors. Simulates motions of individual plasma particles under influence of electromagnetic fields generated by particles themselves. Performed to study variety of nonlinear problems in plasma physics, including magnetic and inertial fusion, plasmas in outer space, propagation of electron and ion beams, free-electron lasers, and particle accelerators.

Liewer, Paulett C.↗
NASA NTRS1990

Global nonlinear optimization of spacecraft protective structures design

The global optimization of protective structural designs for spacecraft subject to hypervelocity meteoroid and space debris impacts is presented. This nonlinear problem is first formulated for weight minimization of the space station core module configuration using the Nysmith impact predictor. Next, the equivalence and uniqueness of local and global optima is shown using properties of convexity. This analysis results in a new feasibility condition for this problem. The solution existence is then shown, followed by a comparison of optimization techniques. Finally, a sensitivity analysis is presented to determine the effects of variations in the systemic parameters on optimal design. The results show that global optimization of this problem is unique and may be achieved by a number of methods, provided the feasibility condition is satisfied. Furthermore, module structural design thicknesses and weight increase with increasing projectile velocity and diameter and decrease with increasing separation between bumper and wall for the Nysmith predictor.

Mog, R. A.↗
NASA NTRS1990

Modeling of vortex-induced oscillations based on indicial response approach

The indicial response approach, a modeling approach originally used for studying nonlinear problems in flight dynamics is applied to the study of vortex-induced oscillations phenomena. The indicial response of the velocity field is derived for the problem studied with emphasis on physical postulates involved. A full account of fluid dynamics effects is taken by considering the incompressible Navier-Stokes equations. The theory is applied to the particular case of flow past a cylinder in periodically forced motion to derive some salient features. The indicial response approach is shown to be equivalent to a currently popular approach based on the use of the amplitude equation. Jump and hysteresis phenomena that experiments indicate occur within the lock-in regime (where the velocity field oscillates at the forcing frequency) are captured by the approach.

Truong, Khiem Van↗
NASA NTRS1990

Dynamical approach study of spurious steady-state numerical solutions of nonlinear differential equations. Part 1: The ODE connection and its implications for algorithm development in computational fluid dynamics

Spurious stable as well as unstable steady state numerical solutions, spurious asymptotic numerical solutions of higher period, and even stable chaotic behavior can occur when finite difference methods are used to solve nonlinear differential equations (DE) numerically. The occurrence of spurious asymptotes is independent of whether the DE possesses a unique steady state or has additional periodic solutions and/or exhibits chaotic phenomena. The form of the nonlinear DEs and the type of numerical schemes are the determining factor. In addition, the occurrence of spurious steady states is not restricted to the time steps that are beyond the linearized stability limit of the scheme. In many instances, it can occur below the linearized stability limit. Therefore, it is essential for practitioners in computational sciences to be knowledgeable about the dynamical behavior of finite difference methods for nonlinear scalar DEs before the actual application of these methods to practical computations. It is also important to change the traditional way of thinking and practices when dealing with genuinely nonlinear problems. In the past, spurious asymptotes were observed in numerical computations but tended to be ignored because they all were assumed to lie beyond the linearized stability limits of the time step parameter delta t. As can be seen from the study, bifurcations to and from spurious asymptotic solutions and transitions to computational instability not only are highly scheme dependent and problem dependent, but also initial data and boundary condition dependent, and not limited to time steps that are beyond the linearized stability limit.

Yee, H. C.↗
NASA NTRS1990

Geometric programming prediction of design trends for OMV protective structures

The global optimization trends of protective honeycomb structural designs for spacecraft subject to hypervelocity meteroid and space debris are presented. This nonlinear problem is first formulated for weight minimization of the orbital maneuvering vehicle (OMV) using a generic monomial predictor. Five problem formulations are considered, each dependent on the selection of independent design variables. Each case is optimized by considering the dual geometric programming problem. The dual variables are solved for in terms of the generic estimated exponents of the monomial predictor. The primal variables are then solved for by conversion. Finally, parametric design trends are developed for ranges of the estimated regression parameters. Results specify nonmonotonic relationships for the optimal first and second sheet mass per unit areas in terms of the estimated exponents.

Mog, R. A.↗
NASA NTRS1990

Innovative design of composite structures: Axisymmetric deformations of unsymmetrically laminated cylinders loaded in axial compression

The study focuses on the axisymmetric deformation response of unsymmetrically laminate cylinders loaded in axial compression by known loads. A geometrically nonlinear analysis is used. Though buckling is not studied, the deformations can be considered to be the prebuckling response. Attention is directed at three 16 layer laminates: a (90 sub 8/0 sub 8) sub T; a (0 sub 8/90 sub 8) sub T and a (0/90) sub 4s. The symmetric laminate is used as a basis for comparison, while the two unsymmetric laminates were chosen because they have equal but opposite bending-stretching effects. Particular attention is given to the influence of the thermally-induced preloading deformations that accompany the cool-down of any unsymmetric laminate from the consolidation temperature. Simple support and clamped boundary conditions are considered. It is concluded that: (1) The radial deformations of an unsymmetric laminate are significantly larger than the radial deformations of a symmetric laminate, although for both symmetric and unsymmetric laminates the large deformations are confined to a boundary layer near the ends of the cylinder; (2) For this nonlinear problem the length of the boundary layer is a function of the applied load; (3) The sign of the radial deformations near the supported end of the cylinder depends strongly on the sense (sign) of the laminate asymmetry; (4) For unsymmetric laminates, ignoring the thermally-induced preloading deformations that accompany cool-down results in load-induced deformations that are under predicted; and (5) The support conditions strongly influence the response but the influence of the sense of asymmetry and the influence of the thermally-induced preloading deformations are independent of the support conditions.

Hyer, M. W.↗
NASA NTRS1990

Parallel processors and nonlinear structural dynamics algorithms and software

Techniques are discussed for the implementation and improvement of vectorization and concurrency in nonlinear explicit structural finite element codes. In explicit integration methods, the computation of the element internal force vector consumes the bulk of the computer time. The program can be efficiently vectorized by subdividing the elements into blocks and executing all computations in vector mode. The structuring of elements into blocks also provides a convenient way to implement concurrency by creating tasks which can be assigned to available processors for evaluation. The techniques were implemented in a 3-D nonlinear program with one-point quadrature shell elements. Concurrency and vectorization were first implemented in a single time step version of the program. Techniques were developed to minimize processor idle time and to select the optimal vector length. A comparison of run times between the program executed in scalar, serial mode and the fully vectorized code executed concurrently using eight processors shows speed-ups of over 25. Conjugate gradient methods for solving nonlinear algebraic equations are also readily adapted to a parallel environment. A new technique for improving convergence properties of conjugate gradients in nonlinear problems is developed in conjunction with other techniques such as diagonal scaling. A significant reduction in the number of iterations required for convergence is shown for a statically loaded rigid bar suspended by three equally spaced springs.

Belytschko, Ted↗