Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 37 records · Page 2

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

Impact damage resistance of flat glass/cellular glass solar mirror experimental panels

An attempt is made to characterize the hail-impact behavior of a solar concentrator structural material which consists of thin, second-surface silvered glass sheets bonded to structural support panels of cellular glass. The resulting glass/cellular glass panel can be described as a thin, brittle plate supported by an elastic foundation, with the adhesive bonding the two as a third component. The plate will bend during impact. Attention is experimentally given to how variations in the characteristics or geometry of the flat glass, adhesive and cellular glass affect hail-impact damage resistance. The critical energy needed to cause impact damage is calculated from the critical velocity of the iceball by means of the kinetic energy formula. The complexity of this structural system and the variability of the specimens render it impossible to completely develop an analytical treatment for the prediction of impact damage resistance.

Varner, J. R.↗

Delamination failure in a unidirectional curved composite laminate

Delamination failure in a unidirectional curved composite laminate was investigated. The curved laminate failed unstably by delaminations developing around the curved region of the laminate at different depths through the thickness until virtually all bending stiffness was lost. Delamination was assumed to initiate at the location of the highest radial stress in the curved region. A closed form curved beam elasticity solution and a 2D finite element analysis (FEA) were conducted to determine this location. The variation in the strain energy release rate, G, with delamination growth was then determined using the FEA. A strength-based failure criteria adequately predicted the interlaminar tension failure which caused initial delamination onset. Using the G analysis the delamination was predicted to extend into the arm and leg of the laminate, predominantly in mode I. As the initial delamination grew around the curved region, the maximum radial stress in the newly formed inner sublaminate increased to a level sufficient to cause a new delamination to initiate in the sublaminate with no increase in applied load. This failure progression was observed experimentally.

Martin, Roderick H.↗

BUCKLING AND POST-BUCKLING OF ELASTIC SHELLS

Some theoretical investigations of buckling of elastic shells are surveyed in this report. Only geometrically perfect shells are considered; initial dents and out-of-roundness are not taken into account. Several questions raised by the studies are: (a) Under what conditions is the infinitesimal theory of buckling of shells adequate? (b) How does the energy theory of buckling of shells correlate with the method based on equilibrium equations for bending moments, tensions, and shears in a buckled configuration? (c) How important are nonlinear terms in the tangential displacements u, v in the strain-displacement relations for buckling and post-buckling studies? (d) How important are the boundary conditions for u, v in affecting stability? (e) If a condition of snap-through is approached, how much external work is required to push the shell "over the hump" into the buckled configuration? (f) How reliable are mathematical approximations used previously in the infinitesimal theory of buckling of shells? Tentative and incomplete answers to some of these questions are suggested.

ELASTIC BUCKLING↗

Delamination failure in a unidirectional curved composite laminate

Delamination failure in a unidirectional curved composite laminate was investigated. The curved laminate failed unstably by delaminations developing around the curved region of the laminate at different depths through the thickness until virtually all bending stiffness was lost. Delamination was assumed to initiate at the location of the highest radial stress in the curved region. A closed form curved beam elasticity solution and a 2-D finite element analysis (FEA) were conducted to determine this location. The variation in the strain energy release rate, G, with delamination growth was then determined using the FEA. A strength-based failure criteria adequately predicted the interlaminar tension failure which caused initial delamination onset. Using the G analysis the delamination was predicted to extend into the arm and leg of the laminate, predominantly in mode I. As the initial delamination grew arould the curved region, the maximum radial stress in the newly formed inner sublaminate increased to a level sufficient to cause a new delamination to initiate in the sublaminate with no increase in applied load. This failure progression was observed experimentally.

Martin, Roderick H.↗

Two–Photon Printing of Shape–Memory Microstructures and Metasurfaces via Radical–Mediated Thiol–Vinyl Hydrothiolation

Shape-memory resists capable of high-resolution curing into arbitrarily designed structures are increasingly demanded for soft robotics, optical sensors, microscale manufacturing, and biomedicine. Amorphous, shape-memory thiol-vinyl networks were printed using two-photon polymerization (2PP) curing of a simple resin formulated with commercially available reagents. The ability to print high-resolution feature sizes down to 200 nm is attributed to the use of radicalmediated, thiol-vinyl step-growth polymerization that quickly crosslinks the resin, limiting diffusive transport. The thermomechanical behavior of the 2PP-cured material analyzed in compression, tension and three-point bending is similar to the behavior of the UV-polymerized samples. To demonstrate the ability to design, field, and test 4D responsive microstructures, an array of nine springs with coil diameters of 330 μm was printed. Following compressive shapefixing, printed arrays can release 11 μJ of stored elastic strain energy when reheated. Further, a new concept of dichroic-memory of a metamaterial device is demonstrated by printing a twisted woodpile structure with circular dichroism as characterized by Mueller Matrix ellipsometry. Furthermore, our results demonstrate how combining high-resolution 2PP curing with stimuli-responsive molecular architectures can further the engineering of responsive microstructures and metamaterials.

36 MATERIALS SCIENCE↗

Nonlinear elasticity and damping govern ultrafast dynamics in click beetles

Many small animals use springs and latches to overcome the mechanical power output limitations of their muscles. Click beetles use springs and latches to bend their bodies at the thoracic hinge and then unbend extremely quickly, resulting in a clicking motion. When unconstrained, this quick clicking motion results in a jump. While the jumping motion has been studied in depth, the physical mechanisms enabling fast unbending have not. Here, we first identify and quantify the phases of the clicking motion: latching, loading, and energy release. We detail the motion kinematics and investigate the governing dynamics (forces) of the energy release. We use high-speed synchrotron X-ray imaging to observe and analyze the motion of the hinge’s internal structures of four Elater abruptus specimens. We show evidence that soft cuticle in the hinge contributes to the spring mechanism through rapid recoil. Using spectral analysis and nonlinear system identification, we determine the equation of motion and model the beetle as a nonlinear single-degree-of-freedom oscillator. Quadratic damping and snap-through buckling are identified to be the dominant damping and elastic forces, respectively, driving the angular position during the energy release phase. The methods used in this study provide experimental and analytical guidelines for the analysis of extreme motion, starting from motion observation to identifying the forces causing the movement. Here, the tools demonstrated here can be applied to other organisms to enhance our understanding of the energy storage and release strategies small animals use to achieve extreme accelerations repeatedly.

59 BASIC BIOLOGICAL SCIENCES↗

A three-dimensional nonlinear Timoshenko beam based on the core-congruential formulation

A three-dimensional, geometrically nonlinear two-node Timoshenkoo beam element based on the total Larangrian description is derived. The element behavior is assumed to be linear elastic, but no restrictions are placed on magnitude of finite rotations. The resulting element has twelve degrees of freedom: six translational components and six rotational-vector components. The formulation uses the Green-Lagrange strains and second Piola-Kirchhoff stresses as energy-conjugate variables and accounts for the bending-stretching and bending-torsional coupling effects without special provisions. The core-congruential formulation (CCF) is used to derived the discrete equations in a staged manner. Core equations involving the internal force vector and tangent stiffness matrix are developed at the particle level. A sequence of matrix transformations carries these equations to beam cross-sections and finally to the element nodal degrees of freedom. The choice of finite rotation measure is made in the next-to-last transformation stage, and the choice of over-the-element interpolation in the last one. The tangent stiffness matrix is found to retain symmetry if the rotational vector is chosen to measure finite rotations. An extensive set of numerical examples is presented to test and validate the present element.

Crivelli, Luis A.↗

Scalability in SLA lattice through lattice orientation and hybrid frame and plate architectures

Lightweighting has been a key goal for engineers and designers, with lattice structures widely explored as the building blocks for structural components. Cellular structure-inspired lattice truss frame designs made of struts, have been extensively studied. All plate-based lattice designs have superior mechanical performance. Limits in scalability occur from formation of closed pockets limiting uv curing in processes like stereolithography (SLA). We examine hybrid frame and plate body-centered cubic-simply cubic (BCC-SC) lattices under compression. Unit cells and scaled up lattice exhibit an increase in yield stress and modulus with the addition of plates. The loading direction on the hybrid frame and plate unit cells affected the magnitude of improvement. Simulations and measurements indicated that the optimal lightweight lattice was determined to be when two plates were placed opposite each other with plates buttressing the struts, inhibiting buckling of the struts aligned with the loading direction. This lattice resulted in 63% improvement in specific modulus, a 137% improvement in specific yield point, and a 360% improvement in specific energy absorption (SEA) and the scaled up 4 × 4 × 4 scaled-up structures, showed a 107%, 148%, and 297% in specific modulus, specific yield point, and SEA, respectively A combined stretching-bending behavior was identified in optimal orientations reflecting the delayed buckling mechanism paired to a rising stress-strain curve past the elastic yield indicating bending resistance. The mass moment of inertia was found to be a key parameter correlating optimum orientation for the same number of plates added to the BCC-SC frame.

Mahan Ghosh↗

OpenFAST Modeling of the T-Omega Wind Floating Offshore Wind Turbine System

To reduce the cost of offshore wind energy through a more efficient design of the floating support structure, T-Omega Wind developed a novel lightweight and shallow-draft platform concept that aims to achieve a wave-following behavior without resonance amplification. The rotor and generator are carried by four tower legs with each leg supported by a shallow-draft float at the base. In preparation for a more detailed analysis, a coupled aero-hydro-elastic model of the proposed concept is developed in OpenFAST. The detailed modeling approach is presented, including a novel application of the SubDyn substructure dynamics module of OpenFAST to approximate the loads on the axle tube at the tower top from the rotor hub bearings. Preliminary results obtained with the OpenFAST model by rigidifying the substructure and blades indicate that the original sizing of the design can lead to large hub acceleration in the axial/surge and pitch directions due to platform-pitch motion. The high axial acceleration also potentially leads to large bending moments in the four tower legs. To address these issues, an updated design is developed with, among other changes, increased float spacing at the tower base and improved shear-transmitting geometry of the tower legs. The new design suggests significantly reduced extreme hub accelerations under the same conditions even with structural flexibility and will be further analyzed in the future through full aero-hydro-servo-elastic simulations.

ENGINEERING,WIND ENERGY↗

OpenFAST Modeling of the T-Omega Wind Floating Offshore Wind Turbine System: Preprint

To reduce the cost of offshore wind energy through a more efficient design of the floating support structure, T-Omega Wind developed a novel lightweight and shallow-draft platform concept that aims to achieve a wave-following behavior without resonance amplification. The rotor and generator are carried by four tower legs with each leg supported by a shallow-draft float at the base. In preparation for a more detailed analysis, a coupled aero-hydro-elastic model of the proposed concept is developed in OpenFAST. The detailed modeling approach is presented, including a novel application of the SubDyn substructure dynamics module of OpenFAST to approximate the loads on the axle tube at the tower top from the rotor hub bearings. Preliminary results obtained with the OpenFAST model by rigidifying the substructure and blades indicate that the original sizing of the design can lead to large hub acceleration in the axial/surge direction and in the pitch direction due to platform-pitch motion. The high hub acceleration also potentially leads to large bending moments in the four tower legs. To address the issue identified, an updated design is developed with, among other changes, increased float spacing at the tower base and improved geometry of the tower legs. The new design suggests significantly reduced extreme platform pitch angles under the same conditions and will be further analyzed in the future through full aero-hydro-servo-elastic simulations.

FOWT↗

OpenFAST Modeling of the T-Omega Wind Floating Offshore Wind Turbine System

To reduce the cost of offshore wind energy through a more efficient design of the floating support structure, T-Omega Wind developed a novel lightweight and shallow-draft platform concept that aims to achieve a wave-following behavior without resonance amplification. The rotor and generator are carried by four tower legs with each leg supported by a shallow-draft float at the base. In preparation for a more detailed analysis, a coupled aero-hydro-elastic model of the proposed concept is developed in OpenFAST. The detailed modeling approach is presented, including a novel application of the SubDyn substructure dynamics module of OpenFAST to approximate the loads on the axle tube at the tower top from the rotor hub bearings. Preliminary results obtained with the OpenFAST model by rigidifying the substructure and blades indicate that the original sizing of the design can lead to large hub acceleration in the axial/surge and pitch directions due to platform-pitch motion. The high axial acceleration also potentially leads to large bending moments in the four tower legs. To address these issues, an updated design is developed with, among other changes, increased float spacing at the tower base and improved shear-transmitting geometry of the tower legs. The new design suggests significantly reduced extreme hub accelerations under the same conditions even with structural flexibility and will be further analyzed in the future through full aero-hydro-servo-elastic simulations.

FOWT↗

Implementation of Creep-Damage Model for Concrete Fracture in MOOSE

Numerical implementation of an isotropic creep-damage model for concrete in multiphysics object-oriented simulation environment (MOOSE) finite element framework is presented in this paper. The constitutive model considers the combined effect of instantaneous and delayed strains on damage propagation. The implementation considers creep using generalized Maxwell or Kelvin-Voigt models. Using strain splitting assumptions, the total mechanical strains are split into elastic and creep components. Damage is considered to evolve as a function of the elastic and creep strains. This work considers damage as a function of fracture energy using the characteristic length of each finite element. This approach preserves the energy release rate of each element and avoids vanishing energy dissipation as the mesh is refined. Further, a creep-damage parameter is used to quantify the effect of creep strain on damage. The model is tested against published results on notched three-point bending specimens involving non-linear creep and predicts that about a third of the creep strain contributes towards damage for the experiments simulated. Results show that the proposed framework has predictive capabilities, and the model can be extended for more complex systems.

36 MATERIALS SCIENCE↗