Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite strains”

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 181 records · Page 10

Mechanical characteristics of filament-wound pressure vessel (burst pressure)

The finite element method is used to analyze the mechanical characteristics of a pressurized filament-wound (FW) pressure vessel, and to predict its burst pressure. The analysis takes into account the bending moment, the stretch-bend coupling effect, nonlinear stress-strain relations, and finite deflection. The analysis is based on two initial failure criteria for laminae, and two ultimate fracture criteria for laminated structures. The numerical results, obtained by applying the load incremental method to the isotensoid CFRP pressure vessel used in the launching of the Zikiken satellite, are in good agreement with the experimental burst pressure and fracture behaviors.

Iida, H.↗

A microstructural lattice model for strain oriented problems: A combined Monte Carlo finite element technique

A specialized, microstructural lattice model, termed MCFET for combined Monte Carlo Finite Element Technique, was developed which simulates microstructural evolution in material systems where modulated phases occur and the directionality of the modulation is influenced by internal and external stresses. In this approach, the microstructure is discretized onto a fine lattice. Each element in the lattice is labelled in accordance with its microstructural identity. Diffusion of material at elevated temperatures is simulated by allowing exchanges of neighboring elements if the exchange lowers the total energy of the system. A Monte Carlo approach is used to select the exchange site while the change in energy associated with stress fields is computed using a finite element technique. The MCFET analysis was validated by comparing this approach with a closed form, analytical method for stress assisted, shape changes of a single particle in an infinite matrix. Sample MCFET analytical for multiparticle problems were also run and in general the resulting microstructural changes associated with the application of an external stress are similar to that observed in Ni-Al-Cr alloys at elevated temperature.

Gayda, J.↗

Analytical and finite element solutions of some problems using a viscoplastic model

The feasibility of using a viscoplastic model developed by Robinson to perform a nonlinear structural analysis was investigated. The paper presents analytical solutions for three classical problems: (1) a pressurized thick-walled cylinder, (2) a thin rotating disk, and (3) a pressurized thick-walled sphere. The analytical expressions derived for the stress and the strain rates for these components are general in nature as they consider both the mechanical and thermal loadings to be time-dependent. A computer program VISTAN (VIscoplastic STructural ANalyzer) was developed to obtain the stress and strain distributions. The finite element solutions for these problems are also presented. The numerical results pertaining to isothermal loading conditions are provided. The results obtained demonstrate the feasibility of using the viscoplastic model developed by Robinson to perform nonlinear structural analyses.

Arya, V. K.↗

Fracture Mechanics Analyses of Reinforced Carbon-Carbon Wing-Leading-Edge Panels

Fracture mechanics analyses of subsurface defects within the joggle regions of the Space Shuttle wing-leading-edge RCC panels are performed. A 2D plane strain idealized joggle finite element model is developed to study the fracture behavior of the panels for three distinct loading conditions - lift-off and ascent, on-orbit, and entry. For lift-off and ascent, an estimated bounding aerodynamic pressure load is used for the analyses, while for on-orbit and entry, thermo-mechanical analyses are performed using the extreme cold and hot temperatures experienced by the panels. In addition, a best estimate for the material stress-free temperature is used in the thermo-mechanical analyses. In the finite element models, the substrate and coating are modeled separately as two distinct materials. Subsurface defects are introduced at the coating-substrate interface and within the substrate. The objective of the fracture mechanics analyses is to evaluate the defect driving forces, which are characterized by the strain energy release rates, and determine if defects can become unstable for each of the loading conditions.

Raju, Ivatury S.↗

The premixed flame in uniform straining flow

Characteristics of the premixed flame in uniform straining flow are investigated by the technique of activation-energy asymptotics. An inverse method is used, which avoids some of the restrictions of previous analyses. It is shown that this method recovers known results for adiabatic flames. New results for flames with heat loss are obtained, and it is shown that, in the presence of finite heat loss, straining can extinguish flames. A stability analysis shows that straining can suppress the cellular instability of flames with Lewis number less than unity. Strain can produce instability of flames with Lewis number greater than unity. A comparison shows quite good agreement between theoretical deductions and experimental observations of Ishizuka, Miyasaka & Law (1981).

Durbin, P. A.↗

Blasim: A computational tool to assess ice impact damage on engine blades

A portable computer called BLASIM was developed at NASA LeRC to assess ice impact damage on aircraft engine blades. In addition to ice impact analyses, the code also contains static, dynamic, resonance margin, and supersonic flutter analysis capabilities. Solid, hollow, superhybrid, and composite blades are supported. An optional preprocessor (input generator) was also developed to interactively generate input for BLASIM. The blade geometry can be defined using a series of airfoils at discrete input stations or by a finite element grid. The code employs a coarse, fixed finite element mesh containing triangular plate finite elements to minimize program execution time. Ice piece is modeled using an equivalent spherical objective that has a high velocity opposite that of the aircraft and parallel to the engine axis. For local impact damage assessment, the impact load is considered as a distributed force acting over a region around the impact point. The average radial strain of the finite elements along the leading edge is used as a measure of the local damage. To estimate damage at the blade root, the impact is treated as an impulse and a combined stress failure criteria is employed. Parametric studies of local and root ice impact damage, and post-impact dynamics are discussed for solid and composite blades.

Reddy, E. S.↗

Mechanical Performance of the First Two Prototype 4.5 m Long Nb$_3$Sn Low-$β$ Quadrupole Magnets for the Hi-Lumi LHC Upgrade

The U.S. High-Luminosity LHC Accelerator Upgrade Project (HL-LHC AUP) team is collaborating with CERN in the design and fabrication of the first 4.5 m long MQXFA magnets, a 150 mm aperture high-field Nb$_3$Sn quadrupole magnet that uses the aluminum shell-based bladder-and-key technology. The first two prototype magnets, MQXFAP1 and MQXFAP2, were assem-bled and tested while the first pre-series structure (MQXFA03) was in fabrication. This paper summarizes the mechanical perfor-mance of these prototype structures based on the comparison of measured strain gauge data with finite element model analyses from all load steps to powering. The MQXFAP1 magnet almost reached ultimate current before a short to ground was detected and the test was stopped. The MQXFAP2 magnet experienced a low training performance due to a fractured aluminum shell. MQXFAP1b was rebuilt with a new replacement coil, but an old coil limited the magnet from achieving the ultimate current. The mitigations and analyses of these prototype magnets are discussed in the context of the transition to pre-series production.

43 PARTICLE ACCELERATORS↗

Strain energy release rate determination of stress intensity factors by finite element methods

The stiffness derivative finite element technique is used to determine the Mode I stress intensity factors for three-crack configurations. The geometries examined include the double edge notch, single edge notch, and the center crack. The results indicate that when the specified guidelines of the Stiffness Derivative Method are used, a high degree of accuracy can be achieved with an optimized, relatively coarse finite element mesh composed of standard, four-node, plane strain, quadrilateral elements. The numerically generated solutions, when compared with analytical ones, yield results within 0.001 percent of each other for the double edge crack, 0.858 percent for the single edge crack, and 2.021 percent for the center crack.

Walsh, R. M., Jr.↗

BLASIM - A computational tool to assess ice impact damage on engine blades

A portable computer code called BLASIM is developed at NASA LeRC to assess the ice impact damage on aircraft engine blades. In addition to the ice impact analyses, the code is also capable of carrying out static, dynamic, resonance margin and flutter analyses. The blade can be solid, hollow, superhybrid or composite material. An optional preprocessor (input generator) is also developed to generate input to the code through interactive process. The blade geometry can be defined either by a series of airfoils at discrete input stations or by a finite element grid. The code employs a coarse fixed finite element mesh with triangular plate finite elements and has quick turnaround time. The ice piece is modeled as an equivalent spherical object and has the velocity opposite to that of the aircraft with direction parallel to the engine axis. For the local impact damage assessment, the impact force is considered as a distributed load acting over a region around the impact point and the average radial strain of the finite elements along the leading edge is taken as a measure of the local damage. To estimate the damage at the blade root, the impact is considered to be an impulse and a combined stress failure criteria is employed. Parametric studies for local and root ice impact damage, and post-impact dynamics are discussed for solid and composite blades.

Reddy, E. S.↗

A finite viscoelastic constitutive model for low to high strain rate response of elastomers with application of strain rate-induced glass transition

Amorphous elastomers exhibit significant rate-stiffening and unique viscous flow characteristics across a wide range of strain rates, often undergoing glass transition above a strain rate threshold. We have developed a thermodynamically-consistent and micromechanically-inspired constitutive model for soft elastomeric materials to capture the rate-dependent stress-strain behavior and hysteresis when subjected to low to high strain rates. Here, our proposed constitutive model encapsulates the viscous flow of materials through molecular motion at low strain rates and local rearrangement and alignment of the molecules trying to overcome the intermolecular resistance at high strain rates, essentially covering the glass transition. We applied our constitutive model to uniaxial compression experiments performed at low and high strain rates for polyborosiloxane (PBS) to identify the material parameters, and subsequently, performed numerical simulations of single and multi-cycle compression, stress relaxation, and small amplitude oscillatory tension-compression. Our analyses indicate that the model predicts higher total energy dissipation with increasing strain rate; however, dissipation associated with molecular relaxation decreases (forming a cusp) because, beyond a crossover strain rate, intermolecular rearrangement and alignment become dominant, which is consistent with the onset of the glass transition. For cyclic loading-unloading, we observed that dissipation over a cycle remains constant at low strain rates but decreases non-monotonically at high strain rates before becoming constant, with the peak stress over the cycle becoming higher, which can be interpreted as more loading being carried elastically by the polymer network as the intermolecular rearrangement process occurs. Additionally, our model was able to predict the qualitative nature of the storage modulus and loss modulus in the limit of small strain over a wide range of frequency sweeps.

36 MATERIALS SCIENCE↗

Finite-element formulations for problems of large elastic-plastic deformation

An Eulerian finite element formulation is presented for problems of large elastic-plastic flow. The method is based on Hill's variational principle for incremental deformations, and is ideally suited to isotropically hardening Prandtl-Reuss materials. Further, the formulation is given in a manner which allows any conventional finite element program, for 'small strain' elastic-plastic analysis, to be simply and rigorously adapted to problems involving arbitrary amounts of deformation and arbitrary levels of stress in comparison to plastic deformation moduli. The method is applied to a necking bifurcation analysis of a bar in plane-strain tension. The paper closes with a unified general formulation of finite element equations, both Lagrangian and Eulerian, for large deformations, with arbitrary choice of the conjugate stress and strain measures. Further, a discussion is given of other proposed formulations for elastic-plastic finite element analysis at large strain, and the inadequacies of some of these are commented upon.

Mcmeeking, R. M.↗

Corotational velocity strain formulations for nonlinear analysis of beams and axisymmetric shells

Finite element formulations for large strain, large displacement problems are formulated using a kinematic description based on the corotational components of the velocity strain. The corotational components are defined in terms of a system that rotates with each element and approximates the rotation of the material. To account for rotations of the material relative to this element system, extra terms are introduced in the velocity strain equations. Although this formulation is incremental, in explicitly integrated transient problems it compares very well with formulations that are not.

Belytschko, T.↗

Flux-Closure Domain Structures in Ferroelectric K 0.5 Na 0.5 NbO 3 Thin Films

Topological domain structures in ferroelectric materials have garnered increasing attention due to their intriguing physical properties and promising applications. While most existing topological structures in ferroelectric perovskite oxides originate from tetragonal or rhombohedral bulk phases, much less is understood about their counterparts in orthorhombic ferroelectrics. Here, in this work, we employ ferroelectric K 0.5 Na 0.5 NbO 3 (KNN) thin films as a model system and leverage phase-field simulations to theoretically predict the static structures and dynamic behaviors of three types of flux-closure domain configurations: in-plane (Type-I), out-of-plane (Type-II), and superdomain (Type-III) flux-closure structures. We systematically investigate the effects of finite size, misfit strains, and electrical boundary conditions on the formation and switching of these topological structures. For the Type-I structure, size reduction or small misfit strain facilitates a transition of the flux-closure pattern to polar vortices. Type-II structures emerge under open-circuit electrical boundary conditions of the film, forming at the junctions of specific domain walls with the film surface or the film–substrate interface. The formation mechanisms of these two flux-closure structures are rationalized from an energy perspective. We demonstrated switching capabilities of Type-II and Type-III structures by obtaining polarization–electric field hysteresis loops. Our simulations also reveal a reversible electric-field-induced transition between the orthorhombic and rhombohedral ferroelectric phases with a checkerboard domain pattern, during which the integrity of the flux-closure structure is preserved. These findings provide theoretical insights and practical guidance for identifying and manipulating topological structures in low-symmetry ferroelectrics, paving the way for developing energy-efficient microelectronic devices based on topological structures.

P-E loop↗

Accumulation of collagen molecular unfolding is the mechanism of cyclic fatigue damage and failure in collagenous tissues

Overuse injuries to dense collagenous tissues are common, but their etiology is poorly understood. The predominant hypothesis that micro-damage accumulation exceeds the rate of biological repair is missing a mechanistic explanation. Here, we used collagen hybridizing peptides to measure collagen molecular damage during tendon cyclic fatigue loading and computational simulations to identify potential explanations for our findings. Our results revealed that triple-helical collagen denaturation accumulates with increasing cycles of fatigue loading, and damage is correlated with creep strain independent of the cyclic strain rate. Finite-element simulations demonstrated that biphasic fluid flow is a possible fascicle-level mechanism to explain the rate dependence of the number of cycles and time to failure. Molecular dynamics simulations demonstrated that triple-helical unfolding is rate dependent, revealing rate-dependent mechanisms at multiple length scales in the tissue. The accumulation of collagen molecular denaturation during cyclic loading provides a long-sought “micro-damage” mechanism for the development of overuse injuries.

59 BASIC BIOLOGICAL SCIENCES↗

Low-cycle fatigue analysis of a cooled copper combustion chamber

A three-dimensional finite element elastoplastic strain analysis was performed for the throat section of regeneratively cooled rocket engine combustion chamber. The analysis included thermal and pressure loads, and the effects of temperature dependent material properties, to determine the strain range corresponding to the engine operating cycle. The strain range was used in conjunction with OFHC copper isothermal fatigue test data to predict engine low-cycle fatigue life. The analysis was performed for chamber configuration and operating conditions corresponding to a hydrogen-oxygen chamber which was fatigue tested to failure at the NASA Lewis Research Center.

Miller, R. W.↗

A triangular element based on generalized potential energy concepts

Stiffness equations are formulated for a doubly-curved triangular thin shell finite element. The strain energy component of the potential energy is first expressed in terms of displacements and displacement gradients with the aid of consistent deep shell strain-displacement equations. The element in-plane and normal displacement fields are approximated by complete cubic polynomials. These functions do not satisfy the interelement displacement admissibility conditions. Satisfaction is forced by the imposition of constraint conditions on the interelement boundaries; the constraints represent the modification of the potential energy. Some numerical results for a pinched cylinder, a cylindrical sphere, and a pinched sphere are examined.

Thomas, G. R.↗