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 775 records · Page 43

Full-field quantitative visualization of shock-driven pore collapse and failure modes in PMMA

The dynamic collapse of pores under shock loading is thought to be directly related to hot spot generation and material failure, which is critical to the performance of porous energetic and structural materials. However, the shock compression response of porous materials at the local, individual pore scale is not well understood. This study examines, quantitatively, the collapse phenomenon of a single spherical void in PMMA at shock stresses ranging from 0.4 to 1.0 GPa. Using a newly developed internal digital image correlation technique in conjunction with plate impact experiments, full-field quantitative deformation measurements are conducted in the material surrounding the collapsing pore for the first time. The experimental results reveal two failure mode transitions as shock stress is increased: (i) the first in situ evidence of shear localization via adiabatic shear banding and (ii) dynamic fracture initiation at the pore surface. Numerical simulations using thermo-viscoplastic dynamic finite element analysis provide insights into the formation of adiabatic shear bands (ASBs) and stresses at which failure mode transitions occur. Further numerical and theoretical modeling indicates the dynamic fracture to occur along the weakened material inside an adiabatic shear band. Finally, analysis of the evolution of pore asymmetry and models for ASB spacing elucidate the mechanisms for the shear band initiation sites, and elastostatic theory explains the experimentally observed ASB and fracture paths based on the directions of maximum shear.

42 ENGINEERING↗

PipeSight: A High-Performance Computing Platform for Pipeline Integrity Management

The Phase I feasibility study completed as part of this project has led to a number of innovative technologies being developed and has laid the foundation for a successful Phase II effort to commercialize a platform for managing the integrity of pipelines for the damage mechanisms of the new, hybrid-energy based economy. To ground the development efforts and direction of the project, an extensive market research and customer discovery effort was undertaken early in Phase I. Through this effort, a number of pipeline owners and operators were interviewed, and the following key findings were discovered about the pipeline industry: • Small pipeline operators do not have the central engineering groups necessary to perform their own independent analysis of inspection data, but instead rely on summarized tally sheets provided to them by inspection service providers. • The time it takes to go from an inspection to a completed engineering assessment, even for small segments of pipeline, can take anywhere from 30-120 days. During this delay, critical threats can (and have been known to) cause failures. • Uncertainty is often not accounted for in the assessment of pipeline integrity. The tally sheets provided by third-party service providers are almost always deterministic in nature, identifying threats that present a concern only to the current (not the future) integrity of the pipeline. • It is uncommon to apply the latest technologies to perform advanced assessments of damaged pipelines. There is a desire to use more advanced analysis capabilities to assess threats. Many pipeline operators indicated that they would often excavate a pipeline to perform an inspection and find that the damage was not as bad as they anticipated, thus using limited resources unnecessarily. Companies are not consistent in their use of inspection data to determine corrosion rates, and those that do only calculate deterministic corrosion rates. • The industry has prominently relied on time-based inspections but has recently started to transition to risk-based inspections. However, there appears to be no uniform guidance on how to do so while properly accounting for all sources of uncertainty. • Companies are not storing inspection data in a manner that allows for the ready determination of temporal trends. • Predictive maintenance principles and practices are beginning to be used by early adopters • Some pipelines are being re-purposed to transport different process fluids than they were designed for, e.g., H 2 and CO 2 rich process streams to serve the new hybrid-energy based economy, which are presenting new integrity concerns for the existing pipeline network that crisscrosses the United States. As a result of these discoveries, we were able to target the development efforts in Phase I to best serve the needs of the industry. In Phase I, we developed a way to correlate multiple large-scale scans of the pipeline to determine a probabilistic corrosion rate that accounts for all sources of error and uncertainty in the inspection process. This probabilistic corrosion rate can be used to predict the future thickness distribution of the pipe wall. We demonstrate how this analysis may be performed in an analytical fashion and has been implemented in such a manner that it can be readily distributed using GPU computing through integration of the Kokkos programming model. We also make a very novel extension of the analytical corrosion rate model to Bayesian Networks (an explainable AI technique) that can account for non-parametric distributions of corrosion rates. With the predictions made above for the probabilistic corrosion rate and corresponding future distribution of the pipe wall thickness, we can assess the integrity of the pipeline through the use of a probabilistic engineering assessment. We developed a novel screening data analysis approach that can rapidly identify ‘hotspots’ (local thin areas) where the integrity of the pipeline is a concern. Once more, we implemented this screening approach in C++ to leverage GPU computing via the Kokkos programming model. After the critical hotspots are identified, we developed a program that can automatically generate an advanced finite element model of the damaged regions. Since the number of damaged regions that require advanced analysis can number in the thousands, we integrated an open-source container-native workflow engine for orchestrating parallel jobs on the cloud. Initially, these advanced numerical models were only designed to account for loading due to internal pressure. However, in a slight pivot from the initial Phase I proposal, we developed a complete pipe stress analysis program (called Simflex) which can simulate the complete pipeline and its response to thermal expansion, pressure, thermal bowing, weight, wind, earthquake, support displacement, support friction and external forces. This pipe stress analysis program was written generically, to handle any piping system, but contains the features needed to model long pipelines (i.e., it incorporates a model for soil mechanics and can account for the nonlinear boundary conditions necessary to simulate long underground pipelines). This pipe stress analysis program can simulate any segment of the pipeline (simple or complex) under any set of conditions and loads, to determine the supplemental loads (axial forces and bending moments) at the location of damage. This enables the most accurate state of stress to be accounted for in the pipeline, which can prove critical when evaluating the integrity of a damaged region. In the process of developing the technologies to perform the integrity assessment of the pipeline, we also extended one of the industry standard approaches for performing the assessment of local thin areas that extend more in the circumferential direction than the longitudinal direction of the pipeline. This approach was presented to the API 579-1/AS ME FFS-1 steering committee in November 2021 for consideration in the next edition of the industry standard for Fitness-For-Service (expected to be released in 2023). To help pipeline operators make decisions with the results on any integrity assessment, we developed a new approach to the life-cycle management of pipelines which uses a Bayesian Decision Network. The network is designed to help pipeline operators plan and prioritize inspection activities and ultimately make smarter, more cost-effective decisions. The Bayesian approach accounts for all sources of uncertainty and carries them through to the final optimal decisions, providing a probabilistic framework for optimizing inspection intervals. The proof-of-concept networks developed in the feasibility study are complete, verified, and are focused on a subset of the pipeline. To expand this novel approach to the scale necessary for an entire network of pipelines in Phase II, we will leverage the DOE-funded Bengi solver for industrial-scale decision making with Bayesian Networks [22]. Once implemented, we will be able to provide the pipeline industry with a much-needed tool for optimal inspection planning using truly explainable artificial intelligence (XAI). To handle all of these advanced capabilities into a cloud-based platform, the architecture of the Equity Engineering Cloud (EEC) was extended to include Argo Workflows, a framework capable of distributing and managing a massive number of jobs that consume their own resources, such that thousands of serial finite element simulations can be run in parallel. As part of this substantial undertaking, we also integrated Argo Continuous Delivery (CD) into the EEC, to aid with the rapid prototyping and iterations that will be imperative to the success of the PipeSight platform’s Agile development process in Phase II. As part of the pipe stress analysis program, we also developed a custom visualizer that leverages the DOE-funded VTK visualization library. We added custom contouring capabilities and a means for interacting visually with both the inputs and outputs of the pipe stress analysis program. We also developed routines for automating the post-processing of the finite element simulations to determine if any failure criteria are met and to visualize the deformations, stresses and strains in ParaView using the exodus II file format (a subset of netCDF).

24 POWER TRANSMISSION AND DISTRIBUTION↗

Biaxial load effects in fracture mechanics

It is found that the standard expressions for elastic stress and displacement in the crack-tip region (i.e., the so-called singular solution) cannot be considered to be approximations that are acceptable in a completely general sense. This conclusion is best illustrated by the instance of a biaxially loaded infinite sheet with a flat horizontal central crack, where the effect of load applied parallel to the plane of the crack appears entirely in the second terms of the series representations for local stresses and displacements. An elastoplastic finite-element analysis of the same biaxially loaded finite specimen geometry shows that the global energy release rate, the J-integral, the plastic stress and strain intensity factors (in the sense of Hilton and Hutchinson), and the size of the crack border region plastic yield, all have pronounced biaxial load dependence.

Liebowitz, H.↗

Development of a laminate fatigue analysis

A fatigue analysis technique developed to predict damage growth in notched laminates is described. Features of the analysis include: criteria to relate matrix failure to cyclic stresses in and between plies; correlation of delamination growth with strain-release rate; and in-plane shear modulus change is related to cyclic shear stresses. A simplified finite element model is used to determine stresses in laminates that contain matrix damage. Failure criteria are integrated with the finite element model to form the fatigue analysis.

Roderick, G. L.↗

Elastic-plastic analysis using a triangular ring element in NASTRAN

An elastic plastic triangular ring element is implemented in NASTRAN computer program. The plane strain problem of partially plastic thick walled cylinder under internal pressure is solved and compared with the earlier finite difference solution. A very good agreement has been reached. In order to demonstrate its application to more general problems, an overloaded thread problem for the British Standard Buttress is examined. The maximum axial and principal stresses are located and their values are determined as functions of loadings.

Chen, P. C. T.↗

Characterization of Mode I and mixed-mode delamination growth in T300/5208 graphite/epoxy

The roles played by Mode I and Mode II strain-energy release rates (G-I and G-II, respectively) in inducing delamination growth under static and fatigue loading were investigated, using T300/5208 graphite/epoxy specimens. Double cantilever beam (DCB) specimens and cracked lap shear (CLS) specimens were used for pure Mode I and mixed-mode tests, respectively. Fatigue-induced delamination growth was characterized by constant-amplitude fatigue tests at a minimum to maximum cyclic load ratio of 0.05 and a frequency of 10 Hz. During the tests, the maximum and minimum strain-energy release rates (Gmax, Gmin) and the delamination growth rate (da/dN) were monitored. Static tests on mixed-mode CLS specimens measured the total strain-energy release rate, which was broken into G-I and G-II components using finite-element analysis. A power-law relationship between da/dN and G-Imax, and da/dN and Gmax were obtained from fatigue test results on DCB and CLS specimens, respectively. The power law for a pure Mode II delamination was derived from CLS results by subtracting the contribution due to G-I.

Ramkumar, R. L.↗

A virtual crack-closure technique for calculating stress intensity factors for cracked three dimensional bodies

A three-dimensional virtual crack-closure technique is presented which calculates the strain energy release rates and the stress intensity factors using only nodal forces and displacements from a standard finite element analysis. The technique is an extension of the Rybicki-Kanninen (1977) method, and it assumes that any continuous function can be approximated by a finite number of straight line segments. Results obtained by the method for surface cracked plates with and without notches agree favorably with previous results.

Shivakumar, K. N.↗

Compression Response of a Sandwich Fuselage Keel Panel With and Without Damage

Results are presented from an experimental and analytical study of a sandwich fuselage keel panel with and without damage. The fuselage keel panel is constructed of graphite-epoxy skins bonded to a honeycomb core, and is representative of a highly loaded fuselage keel structure. The face sheets of the panel contain several terminated or dropped plies along the length of the panel. The results presented provide a better understanding of the load distribution in damaged and undamaged thick-face-sheet composite sandwich structure with dropped plies and of the failure mechanisms of such structure in the presence of low-speed impact damage and discrete-source damage. The impact-damage condition studied corresponds to barely visible impact damage (BVID), and the discrete-source damage condition studied is a notch machined through both face sheets. Results are presented from an impact-damage screening study conducted on another panel of the same design to determine the impact energy necessary to inflict BVID on the panel. Results are presented from compression tests of the panel in three conditions: undamaged; BVID in two locations; and BVID in two locations and a notch through both face sheets. Surface strains in the face sheets of the undamaged panel and the notched panel obtained experimentally are compared with finite element analysis results. The experimental and analytical results suggest that for the damage conditions studied, discrete-source damage influences the structural performance more than BVID.

McGowan, David M.↗

Strain Rate Dependent Deformation and Strength Modeling of a Polymer Matrix Composite Utilizing a Micromechanics Approach

Potential gas turbine applications will expose polymer matrix composites to very high strain rate loading conditions, requiring an ability to understand and predict the material behavior under extreme conditions. Specifically, analytical methods designed for these applications must have the capability of properly capturing the strain rate sensitivities and nonlinearities that are present in the material response. The Ramaswamy-Stouffer constitutive equations, originally developed to analyze the viscoplastic deformation of metals, have been modified to simulate the nonlinear deformation response of ductile, crystalline polymers. The constitutive model is characterized and correlated for two representative ductile polymers. Fiberite 977-2 and PEEK, and the computed results correlate well with experimental values. The polymer constitutive equations are implemented in a mechanics of materials based composite micromechanics model to predict the nonlinear, rate dependent deformation response of a composite ply. Uniform stress and uniform strain assumptions are applied to compute the effective stresses of a composite unit cell from the applied strains. The micromechanics equations are successfully verified for two polymer matrix composites. IM7/977-2 and AS4/PEEK. The ultimate strength of a composite ply is predicted with the Hashin failure criteria that were implemented in the composite micromechanics model. The failure stresses of the two composite material systems are accurately predicted for a variety of fiber orientations and strain rates. The composite deformation model is implemented in LS-DYNA, a commercially available transient dynamic explicit finite element code. The matrix constitutive equations are converted into an incremental form, and the model is implemented into LS-DYNA through the use of a user defined material subroutine. The deformation response of a bulk polymer and a polymer matrix composite are predicted by finite element analyses. The results compare reasonably well to experimental values, with some discrepancies. The discrepancies are at least partially caused by the method used to integrate the rate equations in the polymer constitutive model.

Goldberg, Robert K.↗

Transverse Coefficient of Thermal Expansion Measurements of Carbon Fibers Using ESEM at High Temperatures

The transverse coefficient of thermal expansion (CTE) of single IM7, T1000, and P55 carbon fibers are measured at elevated temperatures. The specimens are prepared by press-fitting fiber tows into 0.7mm-diameter cavity in a graphite disk of 5mm in diameter and 3mm high. The specimens are placed on a crucible in an ESEM, and images of the fiber cross section are taken as the fibers are heated up to 800 C. Holding time, heating and cool down cycles are also introduced. The geometrical changes are measured using a graphics tablet. The change in area/perimeter is calculated to determine the strain and transverse CTE for each fiber. In a complimentary computational effort, displacements and stresses are calculated with finite element models.

Ochoa, O.↗

Using Digital Image Correlation for Material Testing at Cryogenic Temperatures

Recent efforts were made to parametrize an aluminum alloy at cryogenic temperatures to calibrate finite element models simulating material deformation. Digital image correlation (DIC) was used to measure material properties using front/back strains of specimens during tensile tests. Data acquired from DIC was compared to physical extensometer measurements. Two camera pairs (front and back) were mounted via aluminum extrusion to the exterior of an environmental chamber of a hydraulic test stand. Images were taken through single-pane optical glass windows of the environment chamber. Positive pressure bags encased the cameras, attaching to the perimeter of the glass to keep out moisture and prevent windows from frosting. Figure 1 shows the environmental setup and the camera setup during cryogenic testing. Some of the challenges during testing at cryogenic temperatures include the gaseous flow of liquid nitrogen in the optical path of the mounted cameras, frost buildup on the windows as seen in Figure 1b, and frost on the test coupon. The temperature gradient present during the introduction of liquid nitrogen into the chamber can increase the projection error of images captured of the specimen. Figure 2 shows the estimated projection error from DIC as a function of temperature during the cool down process. This presentation addresses some of the concerns and observations when using DIC for testing at cryogenic temperatures using methods of isolated camera system and techniques for maintaining the glass panel clear from frost.

Joseph Cochran↗

Strain Rate Dependent Modeling of Polymer Matrix Composites

A research program is in progress to develop strain rate dependent deformation and failure models for the analysis of polymer matrix composites subject to high strain rate impact loads. Strain rate dependent inelastic constitutive equations have been developed to model the polymer matrix, and have been incorporated into a micromechanics approach to analyze polymer matrix composites. The Hashin failure criterion has been implemented within the micromechanics results to predict ply failure strengths. The deformation model has been implemented within LS-DYNA, a commercially available transient dynamic finite element code. The deformation response and ply failure stresses for the representative polymer matrix composite AS4/PEEK have been predicted for a variety of fiber orientations and strain rates. The predicted results compare favorably to experimentally obtained values.

Goldberg, Robert K.↗

Damage development in neutron-irradiated concrete in a test reactor: Hygro-thermal and mechanical simulations

This paper reports the development of a 3D mesoscale hygro-thermal-mechanical simulation approach to predict damage in concrete irradiated in a test reactor. This framework, developed in MOOSE, considers the effects of elevated temperature, moisture content, and high neutron fluence (energy threshold, E > 0.1 MeV) on the mortar and aggregates separately. The first-stage simulation implements hygro-thermal analysis to determine the temperature and RH inside the specimen as a function of imposed radiation energy. These are used as inputs to the second stage, which considers radiation-induced volumetric expansion (RIVE) of aggregates, and creep, shrinkage, and stress-strain response of mortar to predict the expansion, stresses, and damage in specimens made using different coarse aggregates and subjected to different irradiation times. The irradiation time-dependent damage in the mortar is expressed using an isotropic damage parameter. This multi-physics model serves as a predictive tool for damage quantification in concrete due to neutron irradiation.

36 MATERIALS SCIENCE↗

Constitutive modeling of cyclic plasticity and creep, using an internal time concept

Using the concept of an internal time as related to plastic strains, a differential stress-strain relation for elastoplasticity is rederived, such that (1) the concept of a yield-surface is retained; (2) the definitions of elastic and plastic processes are analogous to those in classical plasticity theory; and (3) its computational implementation, via a 'tangent-stiffness' finite element method and a 'generalized-midpoint-radial-return' stress-integration algorithm, is simple and efficient. Also, using the concept of an internal time, as related to both the inelastic strains as well as the Newtonian time, a constitutive model for creep-plasticity interaction, is discussed. The problem of modeling experimental data for plasticity and creep, by the present analytical relations, as accurately as desired, is discussed. Numerical examples which illustrate the validity of the present relations are presented for the cases of cyclic plasticity and creep.

Watanabe, O.↗

Characterization of the Edge Crack Torsion (ECT) Test for Mode III Fracture Toughness Measurement of Laminated Composites

The edge crack torsion (ECT) test is designed to initiate mode III delamination growth in composite laminates. An ECT specimen is a rectangular laminate, containing an edge delamination at the laminate mid-plane. Torsion load is applied to the specimens, resulting in relative transverse shear sliding of the delaminated faces. The test data reduction schemes are intended to yield initiation values of critical mode III strain energy release rate, G(sub IIIc), that are constant with delamination length. The test has undergone several design changes during its development. The objective of this paper was to determine the suitability of the current ECT test design as a mode III fracture test. To this end, ECT tests were conducted on specimens manufactured from IM7/8552 and specimens made from S2/8552 tape laminates. Several specimens, each with different delamination lengths are tested. Detailed, three-dimensional finite element analyses of the specimens were performed. The analysis results were used to calculate the distribution of mode I, mode II, and mode III strain energy release rate along the delamination front. The results indicated that mode III-dominated delamination growth would be initiated from the specimen center. However, in specimens of both material types, the measured values of G(sub IIIc) exhibited significant dependence on delamination length. Furthermore, there was a large amount of scatter in the data. Load-displacement response of the specimens exhibited significant deviation from linearity before specimen failure. X-radiographs of a sample of specimens revealed that damage was initiated in the specimens prior to failure. Further inspection of the failure surfaces is required to identify the damage and determine that mode III delamination is initiated in the specimens.

Ratcliffe, James G.↗

Solar Sail Topology Variations Due to On-Orbit Thermal Effects

The objective of this research was to predict the influence of non-uniform temperature distribution on solar sail topology and the effect of such topology variations on sail performance (thrust, torque). Specifically considered were the thermal effects due to on orbit attitude control maneuvers. Such maneuvers are expected to advance the sail to a position off-normal to the sun by as much as 35 degrees; a solar sail initially deformed by typical pre-tension and solar pressure loads may suffer significant thermally induced strains due to the non-uniform heating caused by these maneuvers. This on-orbit scenario was investigated through development of an automated analytical shape model that iterates many times between sail shape and sail temperature distribution before converging on a final coupled thermal structural affected sail topology. This model utilizes a validated geometrically non-linear finite element model and a thermal radiation subroutine. It was discovered that temperature gradients were deterministic for the off-normal solar angle cases as were thermally induced strains. Performance effects were found to be moderately significant but not as large as initially suspected. A roll torque was detected, and the sail center of pressure shifted by a distance that may influence on-orbit sail control stability.

Banik, Jeremy A.↗

Analysis of Triangular Rollable and Collapsible (TRAC) Composite Booms under the Effects of Gravity and Twist

Finite element analysis is used to investigate Triangular Rollable and Collapsible (TRAC) composite booms. The TRAC booms are subjected to a suite of loading conditions, including in-plane bending, out-of-plane bending, and axial compression. Simulations containing as-designed boom geometry are performed, as are simulations with geometric deviations associated with the as-built conditions of a sample TRAC boom. Namely, a cross-sectional twist along the length of the boom is incorporated into the analysis to study the effect of the twist on the deformation characteristics of the boom. Twisting in the sample boom was observed after short-term storage in a rolled configuration and may be attributed to differential creep deformation between the inner and outer flanges of the TRAC boom. Simulations are conducted with and without the effects of gravity to understand the significance of gravity on testing conditions. Booms of three lengths are evaluated: 3 meters, 7.5 meters and 30 meters. The mechanical responses of the TRAC booms are discussed with reference to the effects of gravity, boom length and cross-sectional twisting. Gravity is found to have a more significant effect on boom mechanical response as boom length is increased. Cross-sectional twisting is found to have a more significant effect on boom mechanical response in the shorter booms.

High-strain composites↗