Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable stiffness”

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 253 records · Page 14

Configuration optimization of space structures

The objective is to develop a computer aid for the conceptual/initial design of aerospace structures, allowing configurations and shape to be apriori design variables. The topics are presented in viewgraph form and include the following: Kikuchi's homogenization method; a classical shape design problem; homogenization method steps; a 3D mechanical component design example; forming a homogenized finite element; a 2D optimization problem; treatment of volume inequality constraint; algorithms for the volume inequality constraint; object function derivatives--taking advantage of design locality; stiffness variations; variations of potential; and schematics of the optimization problem.

Felippa, Carlos↗

Modeling the Nonlinear, Strain Rate Dependent Deformation of Woven Ceramic Matrix Composites With Hydrostatic Stress Effects Included

An analysis method based on a deformation (as opposed to damage) approach has been developed to model the strain rate dependent, nonlinear deformation of woven ceramic matrix composites with a plain weave fiber architecture. In the developed model, the differences in the tension and compression response have also been considered. State variable based viscoplastic equations originally developed for metals have been modified to analyze the ceramic matrix composites. To account for the tension/compression asymmetry in the material, the effective stress and effective inelastic strain definitions have been modified. The equations have also been modified to account for the fact that in an orthotropic composite the in-plane shear stiffness is independent of the stiffness in the normal directions. The developed equations have been implemented into a commercially available transient dynamic finite element code, LS-DYNA, through the use of user defined subroutines (UMATs). The tensile, compressive, and shear deformation of a representative plain weave woven ceramic matrix composite are computed and compared to experimental results. The computed values correlate well to the experimental data, demonstrating the ability of the model to accurately compute the deformation response of woven ceramic matrix composites.

Goldberg, Robert K.↗

Design of structures for optimum geometry

A method is presented for configuration optimization of finite element structures, given a reasonable initial geometry. The objective is to minimize weight or cost. Design variables include geometric as well as member sizing parameters. The number of elements and joints, and the element-joint relationships are prescribed and are not changed during the optimization process. However, the joint locations are changed. The structure is assumed to be linearly elastic and may be statically indeterminate. Multiple loading conditions are allowed. Constraints include limits on stiffness as well as strength. The method is demonstrated with application to truss design, subject to minimum size, strength, buckling, and displacement constraints. Major design improvements are achieved through configurations changes.

Vanderplaats, G. N.↗

Stacking-sequence optimization for buckling of laminated plates by integer programming

Integer-programming formulations for the design of symmetric and balanced laminated plates under biaxial compression are presented. Both maximization of buckling load for given total thickness, and the minimization of total thickness subject to a buckling constraint are formulated. The design variables that define the stacking sequence of the laminate are zero-one integers. It is shown that the formulation results in a linear optimization problem that can be solved on readily aviable software. Constraints on the stacking sequence such as a limit on the number of contiguous plies of the same orientation and limits on in-plane stiffnesses are easily accommodated. Examples are presented for graphite-epoxy plates under uniaxial and biaxial compression using a commercial software package based on the branch-and-bound algorithm.

Haftka, Raphael T.↗

An analytical/numerical correlation study of the multiple concentric cylinder model for the thermoplastic response of metal matrix composites

The utility of a recently developed analytical micromechanics model for the response of metal matrix composites under thermal loading is illustrated by comparison with the results generated using the finite-element approach. The model is based on the concentric cylinder assemblage consisting of an arbitrary number of elastic or elastoplastic sublayers with isotropic or orthotropic, temperature-dependent properties. The elastoplastic boundary-value problem of an arbitrarily layered concentric cylinder is solved using the local/global stiffness matrix formulation (originally developed for elastic layered media) and Mendelson's iterative technique of successive elastic solutions. These features of the model facilitate efficient investigation of the effects of various microstructural details, such as functionally graded architectures of interfacial layers, on the evolution of residual stresses during cool down. The available closed-form expressions for the field variables can readily be incorporated into an optimization algorithm in order to efficiently identify optimal configurations of graded interfaces for given applications. Comparison of residual stress distributions after cool down generated using finite-element analysis and the present micromechanics model for four composite systems with substantially different temperature-dependent elastic, plastic, and thermal properties illustrates the efficacy of the developed analytical scheme.

Pindera, Marek-Jerzy↗

Design of laminated plates for maximum buckling load

The buckling load of laminated plates having midplane symmetry is maximized for a given total thickness. The thicknesses of the layers are taken as the design variables. Buckling analysis is carried out using the finite element method. The optimality equations are solved by a homotopy method which permits tracing optima as a function of total thickness. It is shown that for any design with a given stacking sequence of ply orientations, there exists a design associated with any other stacking sequence which possesses the same bending stiffness matrix and same total thickness. Hence, from the optimum design for a given stacking sequence, one can directly determine the optimum design for any rearrangement of the ply orientations, and the optimum buckling load is independent of the stacking sequence.

Shin, Yung S.↗

Scaled Aeroelastic Wind Tunnel Model Design with Topology Optimization

The goal of this work is to use numerical aeroelastic topology optimization to design flutter models to be tested in the Transonic Dynamics Tunnel, as part of an envisioned validation experiment for transonic flutter. The topological design develops the structural details of a variable-thickness aluminum plate located at the camber line of a scaled transport aircraft model, in order to satisfy the various conflicting design requirements, including a prescribed flutter boundary, flutter frequency, and static stress levels. Four designs are demonstrated here, where it is shown that the simultaneous design of a cylindrical tip store can help ease the conflict between low stiffness (to allow a wing to flutter in the tunnel envelope) and high strength (to ensure safe structural response). The tip store complicates the geometrical features of the wing, however, and is generally undesirable from the vantage point of a validation experiment.

Aeroelasticity↗

Solution of the subsynchronous whirl problem in the high-pressure hydrogen turbomachinery of the Space Shuttle Main Engine

Subsynchronous whirl of the high-pressure fuel turbopump limited operation of the Space Shuttle Main Engine for some months in early 1976. The means by which this problem was successfully attacked is of particular interest to the rotor-dynamics community, not only because this machinery was the highest power-to-weight ratio known (77,000 hp...760 pounds), but because of the multiple forcing functions involved and the means, both analytical and experimental, which were utilized in separating variables, pointing toward successful solutions, and evaluating results. The general means of identifying fundamental characteristics, analyzing data, and conducting computer investigations are delineated. The results of analysis and testing are discussed. Since whirl inception occurs at a shaft speed greater than twice the first system critical, this was increased by stiffening the shaft and bearing supports; the system damping and system stiffness was additionally increased by proper interstage seal design to the extent that the instability threshold is now beyond the operating range.

Ek, M. C.↗

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

Modeling the Nonlinear, Strain Rate Dependent Deformation of Shuttle Leading Edge Materials with Hydrostatic Stress Effects Included

An analysis method based on a deformation (as opposed to damage) approach has been developed to model the strain rate dependent, nonlinear deformation of woven ceramic matrix composites, such as the Reinforced Carbon Carbon (RCC) material used on the leading edges of the Space Shuttle. In the developed model, the differences in the tension and compression deformation behaviors have also been accounted for. State variable viscoplastic equations originally developed for metals have been modified to analyze the ceramic matrix composites. To account for the tension/compression asymmetry in the material, the effective stress and effective inelastic strain definitions have been modified. The equations have also been modified to account for the fact that in an orthotropic composite the in-plane shear response is independent of the stiffness in the normal directions. The developed equations have been implemented into LS-DYNA through the use of user defined subroutines (UMATs). Several sample qualitative calculations have been conducted, which demonstrate the ability of the model to qualitatively capture the features of the deformation response present in woven ceramic matrix composites.

Goldberg, Robert K.↗

Implicit Formulation of Muscle Dynamics in OpenSim

Astronauts lose bone and muscle mass during spaceflight. Exercise countermeasure is the primary method for counteracting bone and muscle mass loss in space. New spacecraft exercise device concepts are currently being developed for the NASAs new crew exploration vehicle. The NASA Digital Astronaut Project (DAP) uses computational modeling to help determine if the new exercise devices will be effective as countermeasures. The NASA Digital Astronaut Project is developing the ability to utilize predictive simulation to provide insight into the change in kinematics and kinetics with a change in device and gravitational environment (1-g versus 0-g). For example, in space exercise the subject's body weight is applied in addition to the loads prescribed for musculoskeletal maintenance. How and where these loads are applied obviously directly impacts bone and tissue loads. Additionally, due to space vehicle structural requirements, exercise devices are often placed on vibration isolation systems. This changes the apparent impedance or stiffness of the device as seen by the user. Data collection under these conditions is often impractical and limited. Predictive modeling provides a means to have a virtual subject to test hypotheses. Predictive simulation provides a virtual subject for which we are able to perform studies such as sensitivity to device loading and vibration isolation without the need for laboratory kinematic or kinetic test data.Direct Collocation optimization provides an efficient means to perform task based optimization and predictive modeling. It is relatively straight forward to structure a physical exercise task in a Direct Collocation mathematical formulation: perform a motion such that you start at an initial pose, achieve a given amount of deflection i.e a squat, return to the initial pose, and minimize muscle activation cost. Direct Collocation is advantageous in that it does not require numerical integration to evaluate the objective function. Instead, the system dynamics are transformed to discrete time and the optimizer is constrained such that the solution is not considered to be a valid unless the dynamic equations are satisfied at all time points. The simulation and optimization are effectively done simultaneously. Due to the implicit integration, time steps can be more coarse than in a differential equation solver. In a gait scenario this means that that the model constraints and cost function are evaluated at 100 nodes in the gait cycle versus 10,000 integration steps in a variable-step forward dynamic simulation. Furthermore, no time is wasted on accurate simulations of movements that are far from the optimum. Constrained optimization algorithms require a Jacobian matrix that contains the partial derivatives of each of the dynamic constraints with respect to of each of the state and control variables at all time points. This is a large but sparse matrix. An implicit dynamics formulation requires computation of the dynamic residuals f as a function of the states x and their derivatives, and controls u:f(x, dxdt, u) 0If the dynamics of musculoskeletal system are formulated implicitly, the Jacobian elements are often available analytically, eliminating the need for numerical differentiation; this is obviously computationally advantageous. Additionally, implicit formulation of musculoskeletal dynamics do not suffer from singularities from low mass bodies, zero muscle activation, or other stiff system or

physical exercise↗

Load redistribution considerations in the fracture of ceramic matrix composites

Using a macroscopic viewpoint, composite laminae are homogeneous orthotropic solids whose directional strengths are random variables. Incorporation of these random variable strengths into failure models, either interactive or noninteractive, allows for the evaluation of the lamina reliability under a given stress state. Using a noninteractive criterion for demonstration purposes, laminate reliabilities are calculated assuming previously established load sharing rules for the redistribution of load as the failure of laminae occur. The matrix cracking predicted by ACK theory is modeled to allow a loss of stiffness in the fiber direction. The subsequent failure in the fiber direction is controlled by a modified bundle theory. Results are compared with previous models which did not permit separate consideration of matrix cracking, as well as to results obtained from experimental data. The effects of variations from the ideal physical geometry which is normally used to depict the matrix cracking are also studied.

Thomas, David J.↗

Self-Deploying Trusses Containing Shape-Memory Polymers

Composite truss structures are being developed that can be compacted for stowage and later deploy themselves to full size and shape. In the target applications, these smart structures will precisely self-deploy and support a large, lightweight space-based antenna. Self-deploying trusses offer a simple, light, and affordable alternative to articulated mechanisms or inflatable structures. The trusses may also be useful in such terrestrial applications as variable-geometry aircraft components or shelters that can be compacted, transported, and deployed quickly in hostile environments. The truss technology uses high-performance shape-memory-polymer (SMP) thermoset resin reinforced with fibers to form a helical composite structure. At normal operating temperatures, the truss material has the structural properties of a conventional composite. This enables truss designs with required torsion, bending, and compression stiffness. However, when heated to its designed glass transition temperature (Tg), the SMP matrix acquires the flexibility of an elastomer. In this state, the truss can be compressed telescopically to a configuration encompassing a fraction of its original volume. When cooled below Tg, the SMP reverts to a rigid state and holds the truss in the stowed configuration without external constraint. Heating the materials above Tg activates truss deployment as the composite material releases strain energy, driving the truss to its original memorized configuration without the need for further actuation. Laboratory prototype trusses have demonstrated repeatable self-deployment cycles following linear compaction exceeding an 11:1 ratio (see figure).

Schueler, Robert M.↗

Rotor Design Options for Improving V-22 Whirl-Mode Stability

With current tilt-rotor technology, very stiff, thick wings of limited aspect ratio are essential to meet aeroelastic stability requirements, which severely limits cruise efficiency and maximum speed. Larger and more efficient tiltrotors; will need longer and lighter wings, for which whirl-mode flutter is a serious design issue. Numerous approaches to improving the whirl-mode airspeed boundary have been investigated, typically focusing on wing design, active stability augmentation, and variable geometry rotors. The research to be reported applies the much simpler approaches of sweeping the outboard blade sections and extending balance weights at the tip. Improving proprotor whirl-mode stability margins is an ongoing research activity at NASA Ames Research Center. Previous efforts focussed on the XV-15 rotor. This is the first application of these design principles to a full V-22 rotor model in CAMRAD II (a previous publication presented results for a preliminary analytical model). The analysis to be discussed in the proposed paper includes a dual load-path blade model, a compressible-flow aerodynamic model, a representation of swashplate/rotor modal coupling, and other new features. The paper will include discussions of the mechanisms by which blade sweep and tip-mass offsets improve whirl-flutter margins. The V-22 was modeled with CAMRAD II; the rotor model is shown. It is also shows an example swept blade derived from the V-22 rotor system. For this rotor, the primary significance of sweep is the improved whirl-flutter boundary, not the reduced Mach-number effects. An offset tip mass is also shown.

Acree, C. W., Jr.↗

Sliding control of pointing and tracking with operator spline estimation

It is shown how a variable structure control technique could be implemented to achieve precise pointing and good tracking of a deformable structure subject to fast slewing maneuvers. The correction torque that has to be applied to the structure is based on estimates of upper bounds on the model errors. For a rapid rotation of the deformable structure, the elastic response can be modeled by oscillators driven by angular acceleration, and where stiffness and damping coefficients are also angular velocity and acceleration dependent. By transforming this slew-driven elastic dynamics into bilinear form (be regarding the vector made up of the angular velocity, squared angular velocity and angular acceleration components, which appear in the coefficients as the input to the deformation dynamics), an operator spline can be constructed, that gives a low order estimate of the induced disturbance. Moreover, a worst case error bound between the estimated deformation and the unknown exact deformation is also generated, which can be used where required in the sliding control correction.

Dwyer, Thomas A. W., III↗

Model reduction in the physical coordinate system

In the dynamics modeling of a flexible structure, finite element analysis employs reduction techniques, such as Guyan's reduction, to remove some of the insignificant physical coordinates, thus producing a dynamics model that has smaller mass and stiffness matrices. But this reduction is limited in the sense that it removes certain degrees of freedom at a node points themselves in the model. From the standpoint of linear control design, the resultant model is still too large despite the reduction. Thus, some form of the model reduction is frequently used in control design by approximating a large dynamical system with a fewer number of state variables. However, a problem arises from the placement of sensors and actuators in the reduced model, because a model usually undergoes, before being reduced, some form of coordinate transformations that do not preserve the physical meanings of the states. To correct such a problem, a method is developed that expresses a reduced model in terms of a subset of the original states. The proposed method starts with a dynamic model that is originated and reduced in finite element analysis. Then the model is converted to the state space form, and reduced again by the internal balancing method. At this point, being in the balanced coordinate system, the states in the reduced model have no apparent resemblance to those of the original model. Through another coordinate transformation that is developed, however, this reduced model is expressed by a subset of the original states.

Yae, K. Harold↗

Optimization of Residual Stresses in MMC's Using Compensating/Compliant Interfacial Layers. Part 2: OPTCOMP User's Guide

A user's guide for the computer program OPTCOMP is presented in this report. This program provides a capability to optimize the fabrication or service-induced residual stresses in uni-directional metal matrix composites subjected to combined thermo-mechanical axisymmetric loading using compensating or compliant layers at the fiber/matrix interface. The user specifies the architecture and the initial material parameters of the interfacial region, which can be either elastic or elastoplastic, and defines the design variables, together with the objective function, the associated constraints and the loading history through a user-friendly data input interface. The optimization procedure is based on an efficient solution methodology for the elastoplastic response of an arbitrarily layered multiple concentric cylinder model that is coupled to the commercial optimization package DOT. The solution methodology for the arbitrarily layered cylinder is based on the local-global stiffness matrix formulation and Mendelson's iterative technique of successive elastic solutions developed for elastoplastic boundary-value problems. The optimization algorithm employed in DOT is based on the method of feasible directions.

Pindera, Marek-Jerzy↗

Reliability analysis of ceramic matrix composite laminates

At a macroscopic level, a composite lamina may be considered as a homogeneous orthotropic solid whose directional strengths are random variables. Incorporation of these random variable strengths into failure models, either interactive or non-interactive, allows for the evaluation of the lamina reliability under a given stress state. Using a non-interactive criterion for demonstration purposes, laminate reliabilities are calculated assuming previously established load sharing rules for the redistribution of load as the failure of laminae occur. The matrix cracking predicted by ACK theory is modeled to allow a loss of stiffness in the fiber direction. The subsequent failure in the fiber direction is controlled by a modified bundle theory. Results using this modified bundle model are compared with previous models which did not permit separate consideration of matrix cracking, as well as to results obtained from experimental data.

Thomas, David J.↗