Search NASA⌕ Search

SEARCH · Search NASA

Results for “Elastic bending energy”

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 19 records

A novel technique for minimizing energy functional using neural networks

An energy functional describes the equilibrium state of a system. In this work, we present a novel technique, Functional Optimization using Neural Networks (FONN), for minimizing the system’s energy. FONN utilizes neural networks to process information at discrete grid points, considering their interactions with neighboring grid points, to update the state of the system. The training process involves formulating a loss function based on the system’s energy, and with the help of multiple fine-tuning steps, the method employs a progressive energy reduction technique that decreases the energy in multiple steps. FONN’s effectiveness is demonstrated across various problems, including the minimization of the heat and Lyapunov energy. Furthermore, the paper explores the minimization of the elastic bending energy with an area constraint.

97 MATHEMATICS AND COMPUTING↗

Patterning of multicomponent elastic shells by gaussian curvature

Recent findings suggest that shell protein distribution and the morphology of bacterial microcompartments regulate the chemical fluxes facilitating reactions which dictate their biological function. Here, we explore how the morphology and component patterning are coupled through the competition of mean and gaussian bending energies in multicomponent elastic shells that form three-component irregular polyhedra. We observe two softer components with lower bending rigidities allocated on the edges and vertices while the harder component occupies the faces. When subjected to a nonzero interfacial line tension, the two softer components further separate and pattern into subdomains that are mediated by the gaussian curvature. We find that this degree of fractionation is maximized when there is a weaker line tension and when the ratio of bending rigidities between the two softer domains ≈2. Our results reveal a patterning mechanism in multicomponent shells that can capture the observed morphologies of bacterial microcompartments, and moreover, can be realized in synthetic vesicles.

Monte Carlo methods↗

Photomechanical vibration of thin crystals of polar semiconductors

It was found that thin crystals of polar (non-centrosymmetric) semiconductors constitute a new type of photosensitive system in which incident illumination is converted into mechanical energy: thus, illumination-induced elastic deformation (bending) was observed on thin (00.1) CdS and (111) GaAs crystals; furthermore, by employing chopped light the crystals were excited to their resonant vibration (photomechanical vibration); the dependence of the amplitude of this vibration on the energy of the incident radiation was found to be similar to the dependence of the surface photovoltage on the energy of the incident radiation (surface photovoltage spectrum). The present findings are consistent with a model based on light-induced modulation of the piezoelectric surface stresses.

Lagowski, J.↗

Sculpting 2D Crystals via Membrane Contractions before and during Solidification

When phospholipids crystallize within the otherwise fluid membranes of giant unilamellar vesicles, the resulting molecularly thin “2D” solids exhibit great variety in their morphology evolution. For example, within membranes containing moderate amounts of the crystallizing component, crystals grow with a fixed morphology depending on vesicle size. Conversely for membranes containing large amounts of the crystallizing species, we find small compact crystals on vesicles of all sizes. However, on large vesicles, growing crystals sprout flower petals that lengthen progressively. These behaviors result from two combined mechanisms: first, like other 2D solids, the shear rigidity of phospholipid crystals renders them intolerant to morphologies with nonzero Gaussian curvature. As a result and especially at elevated membrane tension, the cost of bending elasticity is reduced at the expense of line energy by the formation of flowers as opposed to compact crystals. Second, the composition-dependent tension rise during cooling relaxes via water permeation of the membrane with a time constant scaling as R2. The amount of crystal formed for a small decrease in temperature determines this composition-dependent increase in stress from thermal contractions versus solidification. Surface Evolver computations were motivated using the predicted tension evolution to develop a processing space that maps to experimental observations for initial and growing crystal morphology. Important variable groups are identified, including a scaled ratio of bending to line energy, a vesicle-size-independent group for membrane contractions, and a time constant for stress relaxation. Though processing stresses ultimately relax, the crystal morphology persists well beyond the processing window.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Upper Atmosphere Research Satellite In-Flight Dynamics

Upper Atmosphere Research Satellite flight data from the first 737 days after launch (September 1991) was used to investigate spacecraft disturbances and responses. The investigation included two in-flight dynamics experiments (approximately three orbits each). Orbital and configuration influences on spacecraft dynamic response were also examined. Orbital influences were due to temperature variation from crossing the Earth's terminator and variation of the solar incident energy as the orbit precessed. During the terminator crossing, the rapid ambient temperature change caused the spacecraft's two flexible appendages to experience thermal elastic bending (thermal snap). The resulting response was dependent upon the orientation of the solar array and the solar incident energy. Orbital influences were also caused by on-board and environmental disturbances and spacecraft configuration changes resulting in dynamic responses which were repeated each orbit. Configuration influences were due to solar array rotation changing spacecraft modal properties. The investigation quantified the spacecraft dynamic response produced by the solar array and high gain antenna harmonic drive disturbances. The solar array's harmonic drive output resonated two solar array modes. Friction in the solar array gear drive provided sufficient energy dissipation which prevented the solar panels from resonating catastrophically; however, the solar array vibration amplitude was excessively large. The resulting vibration had a latitude-specific pattern.

Woodard, Stanley E.↗

Range and strength of mechanical interactions of force dipoles in elastic fiber networks

Mechanical forces generated by myosin II molecular motors drive diverse cellular processes, most notably shape change, division and locomotion. These forces may be transmitted over long range through the cytoskeletal medium – a disordered, viscoelastic network of biopolymers. The resulting cell size scale force chains can in principle mediate mechanical interactions between distant actomyosin units, leading to self-organized structural order in the cell cytoskeleton. Inspired by such force transmission through elastic structures in the cytoskeleton, we consider a percolated fiber lattice network, where fibers are represented as linear elastic elements that can both bend and stretch, and the contractile activity of myosin motors is represented by force dipoles. Then, by using a variety of metrics, we show how two such contractile force dipoles interact with each other through their mutual mechanical deformations of the elastic fiber network. As a prelude to two-dipole interactions, we quantify how forces propagate through the network from a single anisotropic force dipole by analyzing clusters of nodes connected by highly strained bonds, as well as through the decay rate of strain energy with distance from a force dipole. We show that predominant fiber bending screens out force propagation, resulting in reduced and strongly network configuration-dependent dipole interactions. On the other hand, stretching-dominated networks support longer-ranged inter-dipole interactions that recapitulate the predictions of linear elasticity theory. By characterizing the differences between tensile and compressive force propagation in the fiber network, we show how inter-dipole interaction depends on the dipoles’ mutual separation and orientation. The resulting elastic interaction energy may mediate a force between multiple distant dipoles, leading to their self-organization into ordered configurations. In conclusion, this provides a potential pathway for active mechanical force-driven structural order in elastic biopolymer networks.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Sensitivity of inelastic response to numerical integration of strain energy

The exact solution to the quasi-static, inelastic response of a cantilever beam of rectangular cross section subjected to a bending moment at the tip is obtained. The material of the beam is assumed to be linearly elastic-linearly strain-hardening. This solution is then compared with three different numerical solutions of the same problem obtained by minimizing the total potential energy using Gaussian quadratures of two different orders and a Newton-Cotes scheme for integrating the strain energy of deformation. Significant differences between the exact dissipative strain energy and its numerical counterpart are emphasized. The consequence of this on the nonlinear transient responses of a beam with solid cross section and that of a thin-walled beam on elastic supports under impulsive loads are examined.

Kamat, M. P.↗

New O+CO2 Scattering Cross Sections and Their Impact on Photochemical Atmospheric Escape of Oxygen from Mars

Dissociative recombination of O2+ and CO2+ molecular ions in the upper atmosphere of Mars is a major source of suprathermal (hot) oxygen atoms. A significant fraction of hot O atoms are sufficiently fast to overcome Mars' gravitational potential and escape to space. This mechanism is known as photochemical escape (PE) and found to be one of the major escape mechanisms presently active on Mars. The Mars Atmosphere and Volatile Evolution Mission (MAVEN) can constrain the PE rate from measurements of electron temperature and density and ion temperature and density. Escape probabilities of hot O can then be constructed from in-situ neutral densities and escaping O flux loss due to collisions with atmospheric background gases. Due to CO2 abundance in the upper atmosphere, O+CO2 scattering is the most important sink of the escaping hot O flux that strongly affects estimates of total atomic oxygen escape rate obtained from MAVEN data, as outlined above. New O+CO2 elastic and inelastic cross sections were constructed from first principles for collision energies from 0.01-5 eV. O(3P)+CO2 collisions were described using new 4D electronic potential energy surfaces (PESs) for restricted geometries. Three lowest triplet states of CO2 were considered for symmetric, bending, and asymmetric vibrational modes. The cross sections were determined from close-coupling calculations performed on individual surfaces using vibrating rotor approximations for CO2 molecule. Differential and momentum transfer cross sections, of use in energy transfer atmospheric models, were also constructed. The obtained cross sections are in excellent agreement with published results. Our new elastic cross sections are significantly smaller than assumed in earlier atmospheric escape models, resulting in larger photochemical escape rate of O from Mars, as estimated from a 1D escape model. Our results agree well with estimated hot O flux (about 9 x 10(sup 25) s(sup -1)) derived from altitude profiles of 130.4 nm O emission measured up to ~4,000 km by the Imaging Ultraviolet Spectrograph on MAVEN orbiter.

Atmospheric↗

Contribution to the problem of buckling of orthotropic plates, with special reference to plywood

Planar stress-strain relations and bending stress-strain relations are presented for elastic orthotropic plates and specialized to plywood. These relations are used to derive the differential equation and energy expression for the buckling of orthotropic rectangular plates whose principal stiffness directions are not parallel to the plate edges. Buckling analyses are made for the case of pure compression and pure shear of a long plate-strip.

WOODS↗

Effect of load introduction in compression testing of composite laminates

Compression testing of composite materials is affected by the manner in which the compressive load is introduced. Two such effects are studied in this paper: (a) the constrained edge effect, in which transverse expansion of the edges is prevented while the axial load is introduced, and (b) nonuniform gripping, as manifested by inplane bending of the test specimen. The principle of minimum complementary energy is used to develop an analytical model that quantifies these two effects upon the measured elastic properties of laminated composites. Numerical results are presented for selected high-strength graphite/epoxy composites.

Reiss, R.↗

Combining four-point bending and corrosion to map stress-dependent corrosion susceptibility of 316H stainless steel in FLiNaK

Understanding the effect of stress on the corrosion of steels in molten salts is important for material screening and safety assessment of molten salt reactors. We present a method that combines four-point bending with corrosion testing, enabling the mapping of corrosion susceptibility as a function of local stress from a single specimen. Guided by finite element analysis, a four-point bending setup was used to bend a 316H stainless-steel bar under controlled elastic stress during a 100-hour exposure in FLiNaK at 700 °C. Post-exposure characterization using scanning electron microscopy and energy-dispersive X-ray spectroscopy revealed a pronounced correlation between stress and corrosion. Regions under tensile stress exhibited significantly deeper attack depths and greater chromium depletion along grain boundaries, whereas compressive zones showed shallower corrosion cracks. The stress effect is attributed to stress concentration under tensile loading, which promotes chromium dissolution and crack propagation. This technique can be readily applied to other structural alloys, enabling quantitative measurement of corrosion susceptibility as a function of local stress.

316H↗

Transverse shear stiffness of laminated anisotropic shells

Equations are derived for the transverse shear stiffness of laminated anisotropic shells. Without making assumptions for thickness distribution for either transverse shear stresses or strains, constitutive equations for the transverse shear deformation theory of anisotropic heterogeneous shells are found. The equations are based on Taylor series expansions about a generic point for stress resultants and couples, identically satisfying plate equilibrium equations. These equations are used to find statically correct expressions for in-surface stresses, transverse shear stresses, and the area density of transverse shear strain energy, in terms of transverse shear stress resultants and redundants. The application of Castigliano's theorem of least work minimizes shear strain energy with respect to the redundants. Examples are presented for several laminated walls. Good agreement is found between the results and those of exact three-dimensional elasticity solutions for the cylindrical bending of a plate.

Cohen, G. A.↗

Stress Analysis of Composite Cylindrical Shells With an Elliptical Cutout

A special-purpose, semi-analytical solution method for determining the stress and deformation fields in a thin laminated-composite cylindrical shell with an elliptical cutout is presented. The analysis includes the effects of cutout size, shape, and orientation; nonuniform wall thickness; oval-cross-section eccentricity; and loading conditions. The loading conditions include uniform tension, uniform torsion, and pure bending. The analysis approach is based on the principle of stationary potential energy and uses Lagrange multipliers to relax the kinematic admissibility requirements on the displacement representations through the use of idealized elastic edge restraints. Specifying appropriate stiffness values for the elastic extensional and rotational edge restraints (springs) allows the imposition of the kinematic boundary conditions in an indirect manner, which enables the use of a broader set of functions for representing the displacement fields. Selected results of parametric studies are presented for several geometric parameters that demonstrate that analysis approach is a powerful means for developing design criteria for laminated-composite shells.

Nemeth, M. P.↗

Stress Analysis of Composite Cylindrical Shells with an Elliptical Cutout

A special-purpose, semi-analytical solution method for determining the stress and deformation fields in a thin laminated-composite cylindrical shell with an elliptical cutout is presented. The analysis includes the effects of cutout size, shape, and orientation; non-uniform wall thickness; oval-cross-section eccentricity; and loading conditions. The loading conditions include uniform tension, uniform torsion, and pure bending. The analysis approach is based on the principle of stationary potential energy and uses Lagrange multipliers to relax the kinematic admissibility requirements on the displacement representations through the use of idealized elastic edge restraints. Specifying appropriate stiffness values for the elastic extensional and rotational edge restraints (springs) allows the imposition of the kinematic boundary conditions in an indirect manner, which enables the use of a broader set of functions for representing the displacement fields. Selected results of parametric studies are presented for several geometric parameters that demonstrate that analysis approach is a powerful means for developing design criteria for laminated-composite shells.

Oterkus, E.↗

A Nonlinear Theory of Bending and Buckling of Thin Elastic Shallow Spherical Shells

The problem of the finite displacement and buckling, of a shallow spherical dome is investigated both theoretically and experimentally. Experimental results seem to indicate that the classical criterion of buckling is applicable to very shallow spherical domes for which the theoretical calculation was made. A transition to energy criterion for higher domes is also indicated.

Kaplan, A↗

A simple and efficient shear-flexible plate bending element

A shear-flexible triangular element formulation, which utilizes an assumed quadratic displacement potential energy approach and is numerically integrated using Gauss quadrature, is presented. The Reissner/Mindlin hypothesis of constant cross-sectional warping is directly applied to the three-dimensional elasticity theory to obtain a moderately thick-plate theory or constant shear-angle theory (CST), wherein the middle surface is no longer considered to be the reference surface and the two rotations are replaced by the two in-plane displacements as nodal variables. The resulting finite-element possesses 18 degrees of freedom (DOF). Numerical results are obtained for two different numerical integration schemes and a wide range of meshes and span-to-thickness ratios. These, when compared with available exact, series or finite-element solutions, demonstrate accuracy and rapid convergence characteristics of the present element. This is especially true in the case of thin to very thin plates, when the present element, used in conjunction with the reduced integration scheme, outperforms its counterpart, based on discrete Kirchhoff constraint theory (DKT).

Chaudhuri, Reaz A.↗

Extension-torsion coupling behavior of advanced composite tilt-rotor blades

An analytic model was developed to study the extension-bend-twist coupling behavior of an advanced composite helicopter or tilt-rotor blade. The outer surface of the blade is defined by rotating an arbitrary cross section about an initial twist axis. The cross section can be nonhomogeneous and composed of generally anisotropic materials. The model is developed based upon a three dimensional elasticity approach that is recast as a coupled two-dimensional boundary value problem defined in a curvilinear coordinate system. Displacement solutions are written in terms of known functions that represent extension, bending, and twisting and unknown functions for local cross section deformations. The unknown local deformation functions are determined by applying the principle of minimum potential energy to the discretized two-dimensional cross section. This is an application of the Ritz method, where the trial function family is the displacement field associated with a finite element (8-node isoparametric quadrilaterals) representation of the section. A computer program was written where the cross section is discretized into 8-node quadrilateral subregions. Initially the program was verified using previously published results (both three-dimensional elasticity and technical beam theory) for pretwisted isotropic bars with an elliptical cross section. In addition, solid and thin-wall multi-cell NACA-0012 airfoil sections were analyzed to illustrate the pronounced effects that pretwist, initial twist axis location, and spar location has on coupled behavior. Currently, a series of advanced composite airfoils are being modeled in order to assess how the use of laminated composite materials interacts with pretwist to alter the coupling behavior of the blade. These studies will investigate the use of different ply angle orientations and the use of symmetric versus unsymmetric laminates.

Kosmatka, J. B.↗

A theory for the fracture of thin plates subjected to bending and twisting moments

Stress fields near the tip of a through crack in an elastic plate under bending and twisting moments are reviewed assuming both Kirchhoff and Reissner plate theories. The crack tip displacement and rotation fields based on the Reissner theory are calculated. These results are used to calculate the J-integral (energy release rate) for both Kirchhoff and Reissner plate theories. Invoking Simmonds and Duva's (1981) result that the value of the J-integral based on either theory is the same for thin plates, a universal relationship between the Kirchhoff theory stress intensity factors and the Reissner theory stress intensity factors is obtained for thin plates. Calculation of Kirchhoff theory stress intensity factors from finite elements based on energy release rate is illustrated. It is proposed that, for thin plates, fracture toughness and crack growth rates be correlated with the Kirchhoff theory stress intensity factors.

Hui, C. Y.↗