Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

On the application of subcell resolution to conservation laws with stiff source terms

LeVeque and Yee recently investigated a one-dimensional scalar conservation law with stiff source terms modeling the reacting flow problems and discovered that for the very stiff case most of the current finite difference methods developed for non-reacting flows would produce wrong solutions when there is a propagating discontinuity. A numerical scheme, essentially nonoscillatory/subcell resolution - characteristic direction (ENO/SRCD), is proposed for solving conservation laws with stiff source terms. This scheme is a modification of Harten's ENO scheme with subcell resolution, ENO/SR. The locations of the discontinuities and the characteristic directions are essential in the design. Strang's time-splitting method is used and time evolutions are done by advancing along the characteristics. Numerical experiment using this scheme shows excellent results on the model problem of LeVeque and Yee. Comparisons of the results of ENO, ENO/SR, and ENO/SRCD are also presented.

Chang, Shih-Hung↗
NASA NTRS1992

Composite laminates with spatially varying fiber orientations - 'Variable stiffness panel concept'

A solution has been obtained to the plane elasticity problem for a symmetrically laminated composite panel with spatially varying fiber orientations. Since variation of the fiber angles along the length of a composite laminate results in stiffness properties that change as a function of location, the laminates are called variable stiffness panels. An analysis of the stiffness variation and its effect on the elastic response of the panel is presented here. A numerical solution has been obtained using an iterative collocation technique. Corresponding closed-form solutions are given for three different sets of boundary conditions. Two of the cases considered have exact solutions and thus serve to validate the numerical model.

Guerdal, Zafer↗
NASA NTRS1993

X-Ray And Acoustic Measurements Yield Stiffnesses

Analytical technique combines ultrasonic scanning measurements of local velocity of sound in specimen of material with x-ray computed tomographic measurements of local mass density to compute local stiffness of material. Stiffnesses at various locations in specimen then used in finite-element mathematical model of elastic behavior of specimen to compute local stresses, local strains, and overall deformations. Technique enhances value of quantitative nondestructive x-ray and ultrasonic measurements. Especially useful in characterization of carbon/carbon composites and other advanced materials not homogeneous and for which customary simplifying assumption of constant density and/or constant stiffness not valid.

Madaras, Eric Irvine↗
NASA NTRS1992

Stiffness of magnetic bearings subjected to combined static and dynamic loads

This paper investigates the stiffness of a magnetic bearing that is subjected to the combined action of static and dynamic loads. Since their sum cannot exceed the saturation load, a large static load will imply that the bearing can carry only a small dynamic load. This smaller dynamic load together with the practical vibration amplitude define a practical upper bound to the dynamic stiffness. This paper also presents approximate design formulas and curves for this stiffness capacity as a function of the ratio of dynamic and static loads. In addition, it indicates that vibrations larger than a certain gap fraction can destabilize the magnetic bearing. This gap fraction, called the critical gap fraction, depends on the dynamic and static load ratio. For example, if the dynamic load is half of the static load, the use of more than 25 percent of gap can destabilize the bearing.

Rao, D. K.↗
NASA NTRS1991

Improved compression buckling for rectangular composite plates by stiffness tailoring

Buckling resistance is often a controlling criterion in the design of structural elements. Design concepts that lead to increased buckling loads (or strains) can directly lower the structural cost and/or weight by a number of means. This study quantifies the improvements that can be achieved in compression buckling loads of rectangular composite plates by using a simple stiffness tailoring concept. The approach is to position the unidirectional lamina through the thickness and over the planform of the plate so that the buckling load is increased with no loss in in-plane stiffness or increase in weight. Finite element analyses have been used to determine the effects of tailoring on the buckling load of plates with various boundary conditions, aspect ratios, thicknesses, and membrane stiffnesses. Increases in buckling loads (or strains) of nearly 200 percent over the uniform plate buckling loads are shown possible with this tailoring concept.

Biggers, Sherrill B.↗
NASA NTRS1993

An inverse method for computation of structural stiffness distributions of aeroelastically optimized wings

An inverse method has been developed to compute the structural stiffness properties of wings given a specified wing loading and aeroelastic twist distribution. The method directly solves for the bending and torsional stiffness distribution of the wing using a modal representation of these properties. An aeroelastic design problem involving the use of a computational aerodynamics method to optimize the aeroelastic twist distribution of a tighter wing operating at maneuver flight conditions is used to demonstrate the application of the method. This exercise verifies the ability of the inverse scheme to accurately compute the structural stiffness distribution required to generate a specific aeroelastic twist under a specified aeroelastic load.

Schuster, David M.↗
NASA NTRS1993

A unit-cell model of textile composite beams for predicting stiffness properties

Flexural stiffness properties of a textile composite beam are obtained from a finite-element model of the unit cell. Three linearly independent deformations, namely, pure extension, pure bending and pure shear, are applied to the unit cell. The top and bottom surfaces of the beam are assumed to be traction free. Periodic boundary conditions on the lateral boundaries of the unit cell are enforced by multi-point constraint elements. From the forces acting on the unit cell, the flexural stiffness coefficients of the composite beam are obtained. The difficulties in determining the transverse shear stiffness are discussed, and a modified approach is presented. The methods are first verified by applying them to isotropic and bimaterial beams for which the results are known, and then illustrated for a simple plain-weave textile composite.

Sankar, Bhavani V.↗
NASA NTRS1996

Effects of Core Thickness and Fiber Orientation on Composite Beam Stiffness

This laboratory experiment is intended for students in an introductory polymer materials and processes course or engineering materials course. It can be conducted as an introduction to the hand lay-up process, with additional observations regarding the stiffness of the completed composite beams based on core thickness and fiber orientation. Students gain hands-on experience with the hand lay-up process by constructing glass/epoxy composite panels. Each lab group produces a panel with different core thickness or fiber orientation. The panels are then cut into strips and tested for flexural stiffness in a three-point bending fixture. Students plot deflection versus load data for composite beams with two different fiber orientations, two core thicknesses and one beam with laminate plies only (no core). The deflection plots highlight the effects of core thickness and fiber orientation on composite beam stiffness.

Hillsman, Vernon S.↗
NASA NTRS2001

Dynamic Stiffness and Damping Characteristics of a High-Temperature Air Foil Journal Bearing

Using a high-temperature optically based displacement measurement system, a foil air bearing's stiffness and damping characteristics were experimentally determined. Results were obtained over a range of modified Sommerfeld Number from 1.5E6 to 1.5E7, and at temperatures from 25 to 538 C. An Experimental procedure was developed comparing the error in two curve fitting functions to reveal different modes of physical behavior throughout the operating domain. The maximum change in dimensionless stiffness was 3.0E-2 to 6.5E-2 over the Sommerfeld Number range tested. Stiffness decreased with temperature by as much as a factor of two from 25 to 538 C. Dimensionless damping was a stronger function of Sommerfeld Number ranging from 20 to 300. The temperature effect on damping being more qualitative, showed the damping mechanism shifted from viscous type damping to frictional type as temperature increased.

Howard, Samuel A.↗
NASA NTRS2001

Estimation of left ventricular operating stiffness from Doppler early filling deceleration time in humans

Shortened early transmitral deceleration times (E(DT)) have been qualitatively associated with increased filling pressure and reduced survival in patients with cardiac disease and increased left ventricular operating stiffness (K(LV)). An equation relating K(LV) quantitatively to E(DT) has previously been described in a canine model but not in humans. During several varying hemodynamic conditions, we studied 18 patients undergoing open-heart surgery. Transesophageal echocardiographic two-dimensional volumes and Doppler flows were combined with high-fidelity left atrial (LA) and left ventricular (LV) pressures to determine K(LV). From digitized Doppler recordings, E(DT) was measured and compared against changes in LV and LA diastolic volumes and pressures. E(DT) (180 +/- 39 ms) was inversely associated with LV end-diastolic pressures (r = -0.56, P = 0.004) and net atrioventricular stiffness (r = -0.55, P = 0.006) but had its strongest association with K(LV) (r = -0.81, P < 0.001). K(LV) was predicted assuming a nonrestrictive orifice (K(nonrest)) from E(DT) as K(nonrest) = (0.07/E(DT))(2) with K(LV) = 1.01 K(nonrest) - 0.02; r = 0.86, P < 0.001, DeltaK (K(nonrest) - K(LV)) = 0.02 +/- 0.06 mm Hg/ml. In adults with cardiac disease, E(DT) provides an accurate estimate of LV operating stiffness and supports its application as a practical noninvasive index in the evaluation of diastolic function.

Non-NASA Center↗
NASA NTRS2004

Overcoming Geometry-Induced Stiffness with IMplicit-Explicit (IMEX) Runge-Kutta Algorithms on Unstructured Grids with Applications to CEM, CFD, and CAA

My goal is to develop and implement efficient, accurate, and robust Implicit-Explicit Runge-Kutta (IMEX RK) methods [9] for overcoming geometry-induced stiffness with applications to computational electromagnetics (CEM), computational fluid dynamics (CFD) and computational aeroacoustics (CAA). IMEX algorithms solve the non-stiff portions of the domain using explicit methods, and isolate and solve the more expensive stiff portions using implicit methods. Current algorithms in CEM can only simulate purely harmonic (up to lOGHz plane wave) EM scattering by fighter aircraft, which are assumed to be pure metallic shells, and cannot handle the inclusion of coatings, penetration into and radiation out of the aircraft. Efficient MEX RK methods could potentially increase current CEM capabilities by 1-2 orders of magnitude, allowing scientists and engineers to attack more challenging and realistic problems.

Kanevsky, Alex↗
NASA NTRS2007

Low-Friction, High-Stiffness Joint for Uniaxial Load Cell

A universal-joint assembly has been devised for transferring axial tension or compression to a load cell. To maximize measurement accuracy, the assembly is required to minimize any moments and non-axial forces on the load cell and to exhibit little or no hysteresis. The requirement to minimize hysteresis translates to a requirement to maximize axial stiffness (including minimizing backlash) and a simultaneous requirement to minimize friction. In practice, these are competing requirements, encountered repeatedly in efforts to design universal joints. Often, universal-joint designs represent compromises between these requirements. The improved universal-joint assembly contains two universal joints, each containing two adjustable pairs of angular-contact ball bearings. One might be tempted to ask why one could not use simple ball-and-socket joints rather than something as complex as universal joints containing adjustable pairs of angularcontact ball bearings. The answer is that ball-and-socket joints do not offer sufficient latitude to trade stiffness versus friction: the inevitable result of an attempt to make such a trade in a ball-and-socket joint is either too much backlash or too much friction. The universal joints are located at opposite ends of an axial subassembly that contains the load cell. The axial subassembly includes an axial shaft, an axial housing, and a fifth adjustable pair of angular-contact ball bearings that allows rotation of the axial housing relative to the shaft. The preload on each pair of angular-contact ball bearings can be adjusted to obtain the required stiffness with minimal friction, tailored for a specific application. The universal joint at each end affords two degrees of freedom, allowing only axial force to reach the load cell regardless of application of moments and non-axial forces. The rotational joint on the axial subassembly affords a fifth degree of freedom, preventing application of a torsion load to the load cell.

Lewis, James L.↗
NASA NTRS1992

Composite Stiffened Panels: Part I - Stiffness, Thermal Expansion, and Thermal Bending Formulation for Finite Element Analysis

A method is presented for formulating stiffness terms and thermal coefficients of fiber-reinforced composite, stiffened panels for input to finite element analysis (FEA). The method is robust enough to handle panels with general cross sectional shapes, including those which are unsymmetric and/or unbalanced. New thermal coefficients are introduced to quantify panel response from through-the-thickness temperature gradients. Equations are defined for stiffness, thermal expansion, and thermal bending that consider the full complement of membrane, bending, membrane-bending coupling, and transverse shear behavior. A technique of implementing this capability with a single plane of shell finite elements using the MSC/NASTRAN(TradeMark) FEA program is revealed. Finally, an example of a composite, hat-stiffened panel is included to demonstrate errors that occur when an unsymmetric panel is symmetrically formulated as traditionally done. These erroneous values and the correct ones produced from the presented method are listed for the panel's stiffness matrices, thermal expansion and bending vectors, and the forces and moments that are thermally induced from in-plane and through-the-thickness temperature gradients.

Craig S Collier↗
DOE OSTI2024

Variable stiffness magnetic spring

This disclosure presents a new type of variable stiffness magnetic spring, which can have a highly linear translational force characteristic. The variable stiffness is achieved through the rotation of a central magnet. Both positive and negative spring constants can be created. Using an analytic-based field analysis modelling technique, the operating principle and linearity characteristics of the adjustable magnetic spring are studied. The use of a magnetic spring with an adjustable negative spring constant could enable an ocean generator to continuously operate in a resonant state, thereby greatly increasing its power generation capability. The described variable stiffness spring could also be useful in other energy harvesting applications, robotic actuator applications, and/or other applications.

Bird, Jonathan↗
DOE OSTI2025

Data-Driven Closures and Assimilation for Stiff Multiscale Random Dynamics

Here, we introduce a data-driven and physics-informed framework for propagating uncertainty in stiff, multiscale random ordinary differential equations (RODEs) driven by correlated (colored) noise. Unlike systems subjected to Gaussian white noise, a deterministic equation for the joint probability density function (PDF) of RODE state variables does not exist in closed form. Moreover, such an equation would require as many phase-space variables as there are states in the RODE system. To alleviate this curse of dimensionality, we instead derive exact, albeit unclosed, reduced-order PDF (RoPDF) equations for low-dimensional observables/quantities of interest. The unclosed terms take the form of state-dependent conditional expectations, which are directly estimated from data at sparse observation times. However, for systems exhibiting stiff, multiscale dynamics, data sparsity introduces regression discrepancies that compound during RoPDF evolution. This is overcome by introducing a kinetic-like defect term to the RoPDF equation, which is learned by assimilating in sparse, low-fidelity RoPDF estimates. Two assimilation methods are considered, namely nudging and deep neural networks, which are successfully tested against Monte Carlo simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Manufacturing and stiffness constraints for topology optimized periodic structures

Topology optimization (TO) is commonly applied to design the unit cells of periodic structures. For example, metamaterials, lattice structures, phononic crystals (PhC), and photonic crystals (PC) have all been previously designed via TO. Unfortunately, the optimal structures for certain design objectives, e.g., bandgaps, are often impossible to manufacture as they have disconnected regions or “islands” of solid material (ISM) that are not self-supporting. Additionally, designs with enclosed void space (EVS) are problematic for additive manufacturing (AM) since support material or pre-sintered powder cannot be removed after manufacturing. We present a series of constraints that may be incorporated into any TO framework to ensure structures are self-supporting without enclosed voids. Additionally, we employ homogenization-based constraints that allow the designer to tune the elastic stiffness and isotropy of the optimized design. The proposed constraints are evaluated on example microstructures and utilized in a simple optimization test problem to highlight their abilities and limitations so that guidelines for appropriate combinations of constraints may be proposed. Effective constraint combinations are demonstrated on the design of 3D photonic crystals for maximum bandgap subject to manufacturing and stiffness constraints.

42 ENGINEERING↗
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)↗
DOE OSTI2022

Multiscale mechanical design of the lightweight, stiff, and damage-tolerant cuttlebone: A computational study

Cuttlebone, the endoskeleton of cuttlefish, offers an intriguing biological structural model for designing low-density cellular ceramics with high stiffness and damage tolerance. Cuttlebone is highly porous (porosity ~93%) and lightweight (density less than 20% of seawater), constructed mainly by brittle aragonite (95 wt%), but capable of sustaining hydrostatic water pressures over 20 atmospheres and exhibits energy absorption capability under compression comparable to many metallic foams (~4.4 kJ/kg). Here, in this work, we computationally investigate how such remarkable mechanical efficiency is enabled by the multiscale structure of cuttlebone. Using the common cuttlefish, Sepia Officinalis, as a model system, we first conducted high-resolution synchrotron micro-computed tomography (µ-CT) and quantified the cuttlebone's multiscale geometry, including the 3D asymmetric shape of individual walls, the wall assembly patterns, and the long-range structural gradient of walls across the entire cuttlebone (ca. 38 chambers). The acquired 3D structural information enables systematic finite-element simulations, which further reveal the multiscale mechanical design of cuttlebone: at the wall level, wall asymmetry provides optimized energy absorption while maintaining high structural stiffness; at the chamber level, variation of walls (number, pattern, and waviness amplitude) contributes to progressive damage; at the entire skeletal level, the gradient of chamber heights tailors the local mechanical anisotropy of the cuttlebone for reduced stress concentration. Our results provide integrated insights into understanding the cuttlebone's multiscale mechanical design and provide useful knowledge for the designs of lightweight cellular ceramics.

36 MATERIALS SCIENCE↗