SEARCH · Engineering Papers

Results for “stiff problem”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1986

On the symbolic manipulation and code generation for elasto-plastic material matrices

A computerized procedure for symbolic manipulations and FORTRAN code generation of elastoplastic material matrix for finite element applications is presented. Special emphasis is placed on expression simplifications during intermediate derivations, optimal code generation, and interface with the main program. A systematic procedure is outlined to avoid redundant algebraic manipulations. Symbolic expressions of the derived material stiffness matrix are automatically converted to RATFOR code which is then translated into FORTRAN statements through a preprocessor. To minimize the interface problem with the main program, a template file is prepared so that the translated FORTRAN statements can be merged into the file to form a subroutine (or a submodule). Three constitutive models; namely, von Mises plasticity, the Drucker-Prager model, and a concrete plasticity model, are used as illustrative examples.

Chang, T. Y.↗

NASA NTRS1991

On the symbolic manipulation and code generation for elasto-plastic material matrices

A computerized procedure for symbolic manipulations and FORTRAN code generation of an elasto-plastic material matrix for finite element applications is presented. Special emphasis is placed on expression simplifications during intermediate derivations, optimal code generation, and interface with the main program. A systematic procedure is outlined to avoid redundant algebraic manipulations. Symbolic expressions of the derived material stiffness matrix are automatically converted to RATFOR code which is then translated into FORTRAN statements through a preprocessor. To minimize the interface problem with the main program, a template file is prepared so that the translated FORTRAN statements can be merged into the file to form a subroutine (or a submodule). Three constitutive models; namely, von Mises plasticity, Drucker-Prager model, and a concrete plasticity model, are used as illustrative examples.

Chang, T. Y.↗

NASA NTRS1991

Controls of interaction dynamics of orbital assembly

Building structures and spacecraft in orbit will require technologies for positioning, docking/berthing, and joining orbital structures. A fundamental problem underlying the operation of docking and berthing is that of controlling the contact dynamics of mechanical structures actuated by active mechanisms such as robotic devices. Control systems must be designed to control these active mechanisms so that both the free space motions and contact motions are stable and satisfy specifications on position accuracy and bounds on contact forces. For the large orbital structures of the future, the problem of interactive dynamics and control is fundamentally different in several ways than it was for spacecraft docking in the past. First, future space structures must be treated as flexible structures - the operations of docking, berthing, and assembly will need to respect the vibrations of the structures. Second, the assembly of these structures will require multiple-point contact, rather than the essentially single-point positioning of conventional spacecraft docking. Third, some assembly operations require the subassemblies to be brought and held in contact so that successful joining can be accomplished. A preliminary study of contact stability and compliance control design has resulted in the development of an analytical method and a design method to analyze stability. The analytical method analyzes the problem of stability when an actively-controlled structure contacts a passive structure. This method makes it possible to accurately estimate the stiffness of the passive structures with which the contact motion will become unstable. The analytic results suggest that passivity is neither achievable in practice, nor necessary as a design concept. A contact control system need only be passive up to a certain frequency; beyond that frequency the system can be stabilized with sufficiently small gains. With this concept the Center developed a design methodology for achieving desired compliant contact motions. This design method is based on H-infinity norm optimization, which makes it possible to consider both driving point mechanical impedance and systems robustness to modeling uncertainty. A laboratory facility was set up to verify experimentally the analytical and design theory.

Su, Renjeng↗

NASA NTRS2000

Adjoint Techniques for Topology Optimization of Structures Under Damage Conditions

The objective of this cooperative agreement was to seek computationally efficient ways to optimize aerospace structures subject to damage tolerance criteria. Optimization was to involve sizing as well as topology optimization. The work was done in collaboration with Steve Scotti, Chauncey Wu and Joanne Walsh at the NASA Langley Research Center. Computation of constraint sensitivity is normally the most time-consuming step of an optimization procedure. The cooperative work first focused on this issue and implemented the adjoint method of sensitivity computation (Haftka and Gurdal, 1992) in an optimization code (runstream) written in Engineering Analysis Language (EAL). The method was implemented both for bar and plate elements including buckling sensitivity for the latter. Lumping of constraints was investigated as a means to reduce the computational cost. Adjoint sensitivity computation was developed and implemented for lumped stress and buckling constraints. Cost of the direct method and the adjoint method was compared for various structures with and without lumping. The results were reported in two papers (Akgun et al., 1998a and 1999). It is desirable to optimize topology of an aerospace structure subject to a large number of damage scenarios so that a damage tolerant structure is obtained. Including damage scenarios in the design procedure is critical in order to avoid large mass penalties at later stages (Haftka et al., 1983). A common method for topology optimization is that of compliance minimization (Bendsoe, 1995) which has not been used for damage tolerant design. In the present work, topology optimization is treated as a conventional problem aiming to minimize the weight subject to stress constraints. Multiple damage configurations (scenarios) are considered. Each configuration has its own structural stiffness matrix and, normally, requires factoring of the matrix and solution of the system of equations. Damage that is expected to be tolerated is local and represents a small change in the stiffness matrix compared to the baseline (undamaged) structure. The exact solution to a slightly modified set of equations can be obtained from the baseline solution economically without actually solving the modified system.. Shennan-Morrison-Woodbury (SMW) formulas are matrix update formulas that allow this (Akgun et al., 1998b). SMW formulas were therefore used here to compute adjoint displacements for sensitivity computation and structural displacements in damaged configurations.

Akgun, Mehmet A.↗

NASA NTRS1999

Optimization of Aerospace Structure Subject to Damage Tolerance Criteria

The objective of this cooperative agreement was to seek computationally efficient ways to optimize aerospace structures subject to damage tolerance criteria. Optimization was to involve sizing as well as topology optimization. The work was done in collaboration with Steve Scotti, Chauncey Wu and Joanne Walsh at the NASA Langley Research Center. Computation of constraint sensitivity is normally the most time-consuming step of an optimization procedure. The cooperative work first focused on this issue and implemented the adjoint method of sensitivity computation in an optimization code (runstream) written in Engineering Analysis Language (EAL). The method was implemented both for bar and plate elements including buckling sensitivity for the latter. Lumping of constraints was investigated as a means to reduce the computational cost. Adjoint sensitivity computation was developed and implemented for lumped stress and buckling constraints. Cost of the direct method and the adjoint method was compared for various structures with and without lumping. The results were reported in two papers. It is desirable to optimize topology of an aerospace structure subject to a large number of damage scenarios so that a damage tolerant structure is obtained. Including damage scenarios in the design procedure is critical in order to avoid large mass penalties at later stages. A common method for topology optimization is that of compliance minimization which has not been used for damage tolerant design. In the present work, topology optimization is treated as a conventional problem aiming to minimize the weight subject to stress constraints. Multiple damage configurations (scenarios) are considered. Each configuration has its own structural stiffness matrix and, normally, requires factoring of the matrix and solution of the system of equations. Damage that is expected to be tolerated is local and represents a small change in the stiffness matrix compared to the baseline (undamaged) structure. The exact solution to a slightly modified set of equations can be obtained from the baseline solution economically without actually solving the modified system. Sherrnan-Morrison-Woodbury (SMW) formulas are matrix update formulas that allow this. SMW formulas were therefore used here to compute adjoint displacements for sensitivity computation and structural displacements in damaged configurations.

Akgun, Mehmet A.↗

NASA NTRS1974

Analysis of longitudinal thickness-shear stiffness of a monolayer filamentary composite

The longitudinal thickness-shear behavior of a rectangular cross-section prism containing circular cross-section fibers embedded in a matrix is analyzed. The problem is formulated as a Saint-Venant flexure problem modeled by a tip-loaded cantilever composite beam and solution is obtained by the boundary-point least-squares method. It is shown that the longitudinal thickness-shear modulus is much greater than the in-plane shear modulus for a monofilament composite and therefore an independent analysis is necessary to predict the longitudinal thickness-shear behavior of the composite. Numerical results are presented in graphical form for various values of fiber volume fraction and material parameters typical of modern advanced composites.

Chang, S.↗

NASA NTRS2002

Dynamic Stability of Uncertain Laminated Beams Under Subtangential Loads

Because of the inherent complexity of fiber-reinforced laminated composites, it can be challenging to manufacture composite structures according to their exact design specifications, resulting in unwanted material and geometric uncertainties. In this research, we focus on the deterministic and probabilistic stability analysis of laminated structures subject to subtangential loading, a combination of conservative and nonconservative tangential loads, using the dynamic criterion. Thus a shear-deformable laminated beam element, including warping effects, is derived to study the deterministic and probabilistic response of laminated beams. This twenty-one degrees of freedom element can be used for solving both static and dynamic problems. In the first-order shear deformable model used here we have employed a more accurate method to obtain the transverse shear correction factor. The dynamic version of the principle of virtual work for laminated composites is expressed in its nondimensional form and the element tangent stiffness and mass matrices are obtained using analytical integration The stability is studied by giving the structure a small disturbance about an equilibrium configuration, and observing if the resulting response remains small. In order to study the dynamic behavior by including uncertainties into the problem, three models were developed: Exact Monte Carlo Simulation, Sensitivity Based Monte Carlo Simulation, and Probabilistic FEA. These methods were integrated into the developed finite element analysis. Also, perturbation and sensitivity analysis have been used to study nonconservative problems, as well as to study the stability analysis, using the dynamic criterion.

Goyal, Vijay K.↗

DOE OSTI2021

Plastic work constrained elastoplastic topology optimization

An elastoplastic topology optimization framework for limiting plastic work generation while maximizing stiffness is presented. The kinematics and constitutive model are based on finite strain linear isotropic hardening plasticity, and the balance laws are solved using a total Lagrangian finite element formulation. Aggregation of the specific plastic work combined with an adaptive normalization scheme efficiently constrains the maximum specific plastic work. The optimization problem is regularized using an augmented partial differential equation filter, and is solved by the method of moving asymptotes where path-dependent sensitivities are derived using the adjoint method. The numerical examples show a clear dependence on the optimized maximum stiffness structures for different levels of constrained specific plastic work. It is also shown that due to the history dependency of the plasticity, the load path significantly influences the structural performance and optimized topology.

42 ENGINEERING↗

NASA NTRS1983

The design and development of a mounting and jettison assembly for the shuttle orbiter advanced gimbal system

This paper describes the requirements, design development, and qualification of the mounting and jettison assembly (MJA) which serves as the base structure for the advanced gimbal system (AGS) developed for NASA, Marshall Space Flight Center, for use during shuttle missions. An engineering model of the MJA has been built and subjected to the following testing: stiffness and modal characterization, sine and random vibration, and a jettison function and energy release. A qualitative summary of the results and the problems encountered during testing, together with the design solutions, is presented.

Korzeniowski, E. S.↗

NASA NTRS1985

Alternative approximation concepts for space frame synthesis

A method for space frame synthesis based on the application of a full gamut of approximation concepts is presented. It is found that with the thoughtful selection of design space, objective function approximation, constraint approximation and mathematical programming problem formulation options it is possible to obtain near minimum mass designs for a significant class of space frame structural systems while requiring fewer than 10 structural analyses. Example problems are presented which demonstrate the effectiveness of the method for frame structures subjected to multiple static loading conditions with limits on structural stiffness and strength.

Lust, R. V.↗

NASA NTRS1983

Optimal trajectories for aeroassisted orbital transfer

Consideration is given to classical and minimax problems involved in aeroassisted transfer from high earth orbit (HEO) to low earth orbit (LEO). The transfer is restricted to coplanar operation, with trajectory control effected by means of lift modulation. The performance of the maneuver is indexed to the energy expenditure or, alternatively, the time integral of the heating rate. Firist-order optimality conditions are defined for the classical approach, as are a sequential gradient-restoration algorithm and a combined gradient-restoration algorithm. Minimization techniques are presented for the aeroassisted transfer energy consumption and time-delay integral of the heating rate, as well as minimization of the pressure. It is shown that the eigenvalues of the Jacobian matrix of the differential system is both stiff and unstable, implying that the sequential gradient restoration algorithm in its present version is unsuitable. A new method, involving a multipoint approach to the two-poing boundary value problem, is recommended.

Miele, A.↗

NASA NTRS1972

A study of modal coupling procedures for the space shuttle

Various procedures for analytically coupling two or more substructures to obtain modal data for assembly, using data obtained from mode surveys of the individual components are discussed. Particular attention was paid to the applicability of the techniques investigated to the space shuttle, with the possible replacement of full scale mode surveys of the shuttle assembly by component mode surveys along with analytical coupling. The synthesis procedures formulated in this report can all handle redundant connections, although the presence of redundancies does add complication to one of the techniques. All of the procedures use the test data directly, without going through an intermediate analytical model, and no analytical stiffness data is required to supplement the mode survey data. The coupling precedures were initially verified and evaluated by applying them to analytical check problems.

Goldenberg, S.↗

NASA NTRS1978

The ATLAS integrated structural analysis and design software system

The ATLAS system provides an extensive set of integrated technical computer-program modules for the analysis and design of general structural configurations, as well as capabilities that are particularly suited for the aeroelastic design of flight vehicles. The system is based on the stiffness formulation of the finite element structural analysis method and can be executed in batch and interactive computing environments on CDC 6600/CYBER computers. Problem-definition input data are written in an engineering-oriented language using a free field format. Input-data default values, generation options, and data quality checks provided by the preprocessors minimize the amount of data and flowtime for problem definition/verfication. Postprocessors allow selected input and calculated data to be extracted, manipulated, and displayed via on-line and off-line prints or plots for monitoring and verifying problem solutions. The sequence and mode of execution of selected program modules are controlled by a common user-oriented language.

Dreisbach, R. L.↗

NASA NTRS1979

Evaluation of an operating MOD-OA 200 kW wind turbine blade

Operating loads and structural damage were monitored during operation of the MOD-OA electric generating system. The turbine was damaged locally between stations 48 and 125 after 2.8 million rotations. Loads due to degraded yaw stiffness and fretting at rib station 48 were identified as primary to this distress. The repaired blades operated an additional 4.8 million rotations without problems.

Donham, R. E.↗

NASA NTRS1982

The mechanics of fatigue in composite laminates

A status report is presented for the development of the mechanics of fatigue in composite materials. A concensus is noted on the need to construct mechanics treatments on the basis of precisely defined damage states which are both characteristic of the material and independent of as many extensive variables, such as load history, as possible. This approach is in keeping with the use of a single crack as the characteristic damage state in homogeneous materials, which has provided a well defined and precise boundary value problem that can be used in the description and prediction of damaged material response. The utility of stress redistribution concepts in determining residual strength, stiffness, and life in composite laminates is also recognized.

Reifsnider, K. L.↗

NASA NTRS1992

Acoustic intensity calculations for axisymmetrically modeled fluid regions

An algorithm for calculating acoustic intensities from a time harmonic pressure field in an axisymmetric fluid region is presented. Acoustic pressures are computed in a mesh of NASTRAN triangular finite elements of revolution (TRIAAX) using an analogy between the scalar wave equation and elasticity equations. Acoustic intensities are then calculated from pressures and pressure derivatives taken over the mesh of TRIAAX elements. Intensities are displayed as vectors indicating the directions and magnitudes of energy flow at all mesh points in the acoustic field. A prolate spheroidal shell is modeled with axisymmetric shell elements (CONEAX) and submerged in a fluid region of TRIAAX elements. The model is analyzed to illustrate the acoustic intensity method and the usefulness of energy flow paths in the understanding of the response of fluid-structure interaction problems. The structural-acoustic analogy used is summarized for completeness. This study uncovered a NASTRAN limitation involving numerical precision issues in the CONEAX stiffness calculation causing large errors in the system matrices for nearly cylindrical cones.

Hambric, Stephen A.↗

DOE OSTI2025

An iterative dynamic chemical stiffness removal method for reacting flow simulations

Abstract An iterative dynamic chemical stiffness removal method (IDCSR) based on quasi-steady-state approximation (QSSA) is proposed. The IDCSR method is built on a previously developed non-iterative method which has proved to work well for small timestep sizes. A novel iterative procedure is designed in IDCSR to enable explicit time integration of stiff chemistry at relatively large timestep sizes relevant to practical reacting flow simulations. The effectiveness of the iterative procedure is first demonstrated with a toy problem and homogeneous auto-ignition with fixed integration step sizes, showing that larger timestep sizes can be allowed for explicit time integration using IDCSR compared with the previous non-iterative method. IDCSR is then compared with existing explicit chemistry solvers for simulations of homogeneous auto-ignition and shows similar or lower computational cost but significantly higher accuracy across a wide range of timestep sizes. IDCSR is further combined with an automatic adaptive time-stepping scheme for simulations of 0-D homogeneous auto-ignition and a 2-D laminar lifted n -dodecane jet flame. For the 0-D auto-ignition simulations, IDCSR is shown to reduce both the error (by 43%–90%) and computational cost (by 6–15 times) compared with existing explicit solvers, while achieving speed-up factors of up to 400 compared with VODE for a wide range of timestep sizes and reaction mechanisms. For the 2-D jet flame simulations, speed-up factors of 15 and 31 for chemistry integration, and 5 and 9 for overall simulation, are achieved by IDCSR compared with CVODE with and without analytic Jacobian, respectively.

Xu, Chao (ORCID:0000000153074159)↗

NASA NTRS1974

Analysis of twisting stiffness for a multifiber composite layer

The twisting stiffness of a rectangular cross section consisting of a single row of solid circular cross-section fibers embedded in a matrix is analyzed. The problem is formulated as a Dirichlet torsion problem of a multielement region and solved by the boundary-point least-squares method. Numerical results for a single-fiber square cross section compare favorably with previous relaxation-method results. New numerical results for three and five-fiber composites suggest that the torsional rigidity of a multifiber composite can be approximated from the torsional rigidities of single and three-fiber models.

Bert, C. W.↗