Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element simulation”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Experimental Evidence of Mechanical Isotropy in Porcine Lung Parenchyma

Pulmonary injuries are a major source of morbidity and mortality associated with trauma. Trauma includes injuries associated with accidents and falls as well as blast injuries caused by explosives. The prevalence and mortality of these injuries has made research of pulmonary injury a major priority. Lungs have a complex structure, with multiple types of tissues necessary to allow successful respiration. The soft, porous parenchyma is the component of the lung which contains the alveoli responsible for gas exchange. Parenchyma is also the portion which is most susceptible to traumatic injury. Finite element simulations are an important tool for studying traumatic injury to the human body. These simulations rely on material properties to accurately recreate real world mechanical behaviors. Previous studies have explored the mechanical properties of lung tissues, specifically parenchyma. These studies have assumed material isotropy but, to our knowledge, no study has thoroughly tested and quantified this assumption. This study presents a novel methodology for assessing isotropy in a tissue, and applies these methods to porcine lung parenchyma. Briefly, lung parenchyma samples were dissected so as to be aligned with one of the three anatomical planes, sagittal, frontal, and transverse, and then subjected to compressive mechanical testing. Stress-strain curves from these tests were statistically compared by a novel method for differences in stresses and strains at percentages of the curve. Histological samples aligned with the anatomical planes were also examined by qualitative and quantitative methods to determine any differences in the microstructural morphology. Our study showed significant evidence to support the hypothesis that lung parenchyma behaves isotropically.

Materials Science↗
DOE OSTI2022

Study of the Printability, Microstructures, and Mechanical Performances of Laser Powder Bed Fusion Built Haynes 230

The nickel-based superalloy, Haynes 230 (H230), is widely used in high-temperature applications, e.g., heat exchangers, because of its excellent high-temperature mechanical properties and corrosion resistance. As of today, H230 is not yet in common use for 3D printing, i.e., metal additive manufacturing (AM), primarily because of its hot cracking tendency under fast solidification. The ability to additively fabricate components in H230 attracts many applications that require the additional advantages leveraged by adopting AM, e.g., higher design complexity and faster prototyping. In this study, we fabricated nearly fully dense H230 in a laser powder bed fusion (L-PBF) process through parameter optimization. The efforts revealed the optimal process space which could guide future fabrication of H230 in various metal powder bed fusion processes. The metallurgical analysis identified the cracking problem, which was resolved by increasing the pre-heat temperature from 80 °C to 200 °C. A finite element simulation suggested that the pre-heat temperature has limited impacts on the maximum stress experienced by each location during solidification. Additionally, the crack morphology and the microstructural features imply that solidification and liquation cracking are the more probable mechanisms. Both the room temperature tensile test and the creep tests under two conditions, (a) 760 °C and 100 MPa and (b) 816 °C and 121 MPa, confirmed that the AM H230 has properties comparable to its wrought counterpart. The fractography showed that the heat treatment (anneal at 1200 °C for 2 h, followed by water quench) balances the strength and the ductility, while the printing defects did not appreciably accelerate part failure.

14 SOLAR ENERGY↗
DOE OSTI2017

Crystal Plasticity Modeling and Experimental Validation with an Orientation Distribution Function for Ti-7Al Alloy

An orientation distribution function based model is used for micromechanical modeling of the titanium-aluminum alloys, Ti-0 wt % Al and Ti-7 wt % Al, which are in demand for many aerospace applications. This probability descriptor based modeling approach is different than crystal plasticity finite element techniques since it computes the averaged material properties using upper bound averaging. A rate-independent single-crystal plasticity model is implemented to compute the effect of macroscopic strain on the polycrystal. An optimization problem is defined for calibrating the basal, prismatic, pyramidal slip system and twin parameters using the available tension and compression experimental data. The crystal plasticity parameters of Ti-7 wt % Al are not studied extensively in literature, and therefore the optimization results for the crystal plasticity model realization produce unique data, which will be beneficial to future studies in the field. The sensitivities of the slip and twin parameters to the design objectives are also investigated to identify the most critical slip system parameters. Using the optimum design parameters, the microstructural textures, during the tension test, are predicted by the crystal plasticity finite element simulations, and compared to the available experimental texture and scanning electron microscope—digital image correlation data.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Investigation of Goubau Wave Propagation on Large Pipes for Sensing Applications

We examine the application of guided waves on a single conductor (Goubau waves) for sensing. In particular, the use of such waves to remotely interrogate surface acoustic wave (SAW) sensors mounted on large-radius conductors (pipes) is considered. Experimental results using a small-radius (0.0032 m) conductor at 435 MHz are reported. The applicability of published theory to conductors of large radius is examined. Finite element simulations are then used to study the propagation and launching of Goubau waves on steel conductors up to 0.254 m in radius. Simulations show that waves can be launched and received, although energy loss into radiating waves is a problem with current launcher designs.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Investigating the Electromechanical Sensitivity of Carbon-Nanotube-Coated Microfibers

The piezoresistance of carbon nanotube (CNT)-coated microfibers is examined using diametric compression. Diverse CNT forest morphologies were studied by changing the CNT length, diameter, and areal density via synthesis time and fiber surface treatment prior to CNT synthesis. Large-diameter (30–60 nm) and relatively low-density CNTs were synthesized on as-received glass fibers. Small-diameter (5–30 nm) and-high density CNTs were synthesized on glass fibers coated with 10 nm of alumina. The CNT length was controlled by adjusting synthesis time. Electromechanical compression was performed by measuring the electrical resistance in the axial direction during diametric compression. Gauge factors exceeding three were measured for small-diameter (<25 μm) coated fibers, corresponding to as much as 35% resistance change per micrometer of compression. The gauge factor for high-density, small-diameter CNT forests was generally greater than those for low-density, large-diameter forests. A finite element simulation shows that the piezoresistive response originates from both the contact resistance and intrinsic resistance of the forest itself. The change in contact and intrinsic resistance are balanced for relatively short CNT forests, while the response is dominated by CNT electrode contact resistance for taller CNT forests. These results are expected to guide the design of piezoresistive flow and tactile sensors.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Limiting Rotation Rate of Neutron Stars from Crust Breaking and Gravitational Waves

Abstract Neutron stars are not observed to spin faster than about half their breakup rate. This limiting rotational frequency may be related to the strength of their crusts. As a star spins up from accretion, centrifugal forces stress the crust. We perform finite-element simulations of rotating neutron stars and find that the crust fails at rotation rates about half the breakup rate. Given uncertainties in microphysics, we have not determined the crust configuration after this failure. Instead, we argue that the crust may fail in an asymmetric way and could produce a configuration with a significant ellipticity (fractional difference in moments of inertia). If the ellipticity is large, a rotating star will radiate gravitational waves that may limit further spin-up. These stars may be promising sources for LIGO/VIRGO and next-generation gravitational-wave detectors.

Morales, J. A.↗
DOE OSTI2017

The Role of Second Phase Hard Particles on Hole Stretchability of two AA6xxx Alloys

The hole stretchability of two Aluminum Alloys (AA6111 and AA6022) are studied by using a two stages integrated finite element framework where the edge geometry and edge damages from the hole piercing processes were considered in the subsequent hole expansion processes. Experimentally it has been found that AA6022 has higher hole expansion ratios than those of AA6111. This observation has been nicely captured by finite element simulations. The main cause of differences have been identified to the volume fractions of the random distributed second phase hard particles which play a critical role in determining the fracture strains of the materials.

Hu, Xiaohua↗
DOE OSTI2022

Local Thermomechanical Processing for Improving Formability of High Strength Aluminum Sheets

Limited room temperature formability hinders the wide-spread use of high strength aluminum alloys in structural body-in-white parts. Stamping or extrusion at warm temperatures or from softer tempers are the current solutions. In this work, our approach is to start with age-hardened sheets from 7xxx and 6xxx family of alloys and improve their formability using local thermomechanical processing only in the regions demanding highest ductility in the forming processes. We achieved local formability improvements with friction stir processing and introduce another process named roller-bending as a concept and showed its feasibility through finite element simulations. Initial results from FSP indicated significant deformation in the processed zones with minimal sheet distortion. FSP also resulted in dynamically recrystallized, fine grained (d < 5 µm) microstructures in the processed regions with textures significantly different from the base material. This resulted in formability improvements of > 30% in 7085-T76 and > 50% in 6111-T6 alloys as measured by minimum bend radius (r/t) ratios established with V-bending tests performed at 90° and 180° bend angles. Based on the encouraging initial results, we envision applicability of these processes in both pre- and post-stamping stages through integration of the process tools as end-effectors to robotic arms.

Nasim, Wahaz↗
DOE OSTI2022

pyvisco [SWR-22-30]

pyvisco is a Python library that supports the identification of Prony series parameters for linear viscoelastic materials described by a Generalized Maxwell model. The necessary material model parameters are identified by fitting a Prony series to the experimental measurement data. pyvisco allows for the identification of Prony series parameters from experimental data measured in either the frequency-domain (via Dynamic Mechanical Thermal Analysis) or time-domain (via relaxation measurements). The experimental data can be provided as raw measurement sets at different temperatures or as pre-processed master curves. An optional minimization routine is included to reduce the number of Prony elements. This routine is helpful in Finite Element simulations where reducing the computational complexity of the linear viscoelastic material models can shorten the simulation time. See also, https://pypi.org/project/pyvisco/

Springer, Martin↗
NASA NTRS1976

Finite-element analysis of dynamic fracture

Applications of the finite element method to the two dimensional elastodynamics of cracked structures are presented. Stress intensity factors are computed for two problems involving stationary cracks. The first serves as a vehicle for discussing lumped-mass and consistent-mass characterizations of inertia. In the second problem, the behavior of a photoelastic dynamic tear test specimen is determined for the time prior to crack propagation. Some results of a finite element simulation of rapid crack propagation in an infinite body are discussed.

Aberson, J. A.↗
NASA NTRS1980

Self-adaptive incremental Newton-Raphson algorithms

Multilevel self-adaptive Newton-Raphson type strategies are developed to improve the solution efficiency of nonlinear finite element simulations of statically loaded structures. The overall strategy involves three basic levels. The first level involves preliminary solution tunneling via primative operators. Secondly, the solution is constantly monitored via quality/convergence/nonlinearity tests. Lastly, the third level involves self-adaptive algorithmic update procedures aimed at improving the convergence characteristics of the Newton-Raphson strategy. Numerical experiments are included to illustrate the results of the procedure.

Padovan, J.↗
NASA NTRS1984

On the development of hierarchical solution strategies for nonlinear finite element formulations

This paper develops a hierarchical type solution scheme which can handle the field equations associated with nonlinear finite element simulations. The overall procedure possesses various levels of application namely degree of freedom, nodal, elemental, substructural as well as global. In particular iteration, updating, assembly and solution control occurs at the various hierarchical levels. Due to the manner of formulation, the degree of matrix inversion depends on the size of the various hierarchical partitioned groups. In this context, degree of freedom partitioning requires no inversion. To benchmark the overall scheme, the results of several numerical examples are presented.

Padovan, J.↗
NASA NTRS1985

Orbit on demand - Structural analysis finds vertical launchers weigh less

Structural considerations arising from favored design concepts for the next generation on-demand launch vehicles are explored. The two emerging concepts are a two stage fully reusable vertical take-off vehicle (V-2) and a horizontal take-off, two stage subsonic boost launch vehicle (H-2-Sub). Both designs have an 1100 n. mi. cross-range capability, with the V-2 orbiter having small wings with winglets for hypersonic trim and the H-2-Sub requiring larger, swept wings. The rockets would be cryogenic, while airbreathing initial boosters would be either turbofans, turbojets and/or ramjets. Dynamic loading is lower in the launch of a V-2. The TPS is a critical factor due to thinner leading edges than on the Shuttle and may require heat-pipe cooling. Airframe structures made of metal matrix composites have passed finite element simulations of projected loads and can now undergo proof-of-concept tests, although whisker-reinforced materials may be superior once long-whisker technology is developed.

Taylor, A. H.↗
NASA NTRS1986

Sublaminate analysis of interlaminar fracture in composites

A simple analysis method based upon a transverse shear deformation theory and a sublaminate approach is utilized to analyze a mixed-mode edge delamination specimen. The analysis provides closed form expressions for the interlaminar shear stresses ahead of the crack, the total energy release rate, and the energy release rate components. The parameters controlling the behavior are identified. The effect of specimen stacking sequence and delamination interface on the strain energy release rate components is investigated. Results are compared with a finite element simulation for reference. The simple nature of the method makes it suitable for preliminary design analyses which require a large number of configurations to be evaluated quickly and economically.

Armanios, E. A.↗
NASA NTRS1986

Interlaminar analysis of laminated composites using a sublaminate approach

A simple analysis method based upon a transverse shear deformation theory and a sublaminate approach is utilized to analyze a Mixed-Mode edge delamination specimen. The analysis provides a closed form distribution of the interlaminar shear stresses ahead of the crack, and the parameters controlling the behavior are identified. The effect of specimen stacking sequence and delamination interface on the strain energy release rate components is investigated. Results are compared with a finite element simulation for reference. The simple nature of the method makes it suitable for preliminary design analyses which require a large number of configurations to be evaluated quickly and economically.

Armanios, E. A.↗
NASA NTRS1986

Constrained hierarchical least square nonlinear equation solvers

The current paper develops a constrained hierarchical least square nonlinear equation solver. The procedure can handle the response behavior of systems which possess indefinite tangent stiffness characteristics. Due to the generality of the scheme, this can be achieved at various hierarchical application levels. For instance, in the case of finite element simulations, various combinations of either degree of freedom, nodal, elemental, substructural, and global level iterations are possible. Overall, this enables a solution methodology which is highly stable and storage efficient. To demonstrate the capability of the constrained hierarchical least square methodology, benchmarking examples are presented which treat structure exhibiting highly nonlinear pre- and postbuckling behavior wherein several indefinite stiffness transitions occur.

Padovan, J.↗
NASA NTRS1987

Structural modeling for multicell composite rotor blades

Composite material systems are currently good candidates for aerospace structures, primarily for the design flexibility they offer, i.e., it is possible to tailor the material and manufacturing approach to the application. A working definition of elastic or structural tailoring is the use of structural concept, fiber orientation, ply stacking sequence, and a blend of materials to achieve specific performance goals. In the design process, choices of materials and dimensions are made which produce specific response characteristics, and which permit the selected goals to be achieved. Common choices for tailoring goals are preventing instabilities or vibration resonances or enhancing damage tolerance. An essential, enabling factor in the design of tailored composite structures is structural modeling that accurately, but simply, characterizes response. The objective of this paper is to present a new multicell beam model for composite rotor blades and to validate predictions based on the new model by comparison with a finite element simulation in three benchmark static load cases.

Rehfield, Lawrence W.↗