Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite strain”

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 415 records · Page 23

Dynamic Shape Reconstruction of Three-Dimensional Frame Structures Using the Inverse Finite Element Method

A robust and efficient computational method for reconstructing the three-dimensional displacement field of truss, beam, and frame structures, using measured surface-strain data, is presented. Known as shape sensing , this inverse problem has important implications for real-time actuation and control of smart structures, and for monitoring of structural integrity. The present formulation, based on the inverse Finite Element Method (iFEM), uses a least-squares variational principle involving strain measures of Timoshenko theory for stretching, torsion, bending, and transverse shear. Two inverse-frame finite elements are derived using interdependent interpolations whose interior degrees-of-freedom are condensed out at the element level. In addition, relationships between the order of kinematic-element interpolations and the number of required strain gauges are established. As an example problem, a thin-walled, circular cross-section cantilevered beam subjected to harmonic excitations in the presence of structural damping is modeled using iFEM; where, to simulate strain-gauge values and to provide reference displacements, a high-fidelity MSC/NASTRAN shell finite element model is used. Examples of low and high-frequency dynamic motion are analyzed and the solution accuracy examined with respect to various levels of discretization and the number of strain gauges.

Gherlone, Marco↗

Evaluation of advanced materials through experimental mechanics and modelling

Composite materials have been frequently used in aerospace vehicles. Very often defects are inherited during the manufacture and damages are inherited during the construction and services. It becomes critical to understand the mechanical behavior of such composite structure before it can be further used. One good example of these composite structures is the cylindrical bottle of solid rocket motor case with accidental impact damages. Since the replacement of this cylindrical bottle is expensive, it is valuable to know how the damages affects the material, and how it can be repaired. To reach this goal, the damage must be characterized and the stress/strain field must be carefully analyzed. First the damage area, due to impact, is surveyed and identified with a shearography technique which uses the principle of speckle shearing interferometry to measure displacement gradient. Within the damage area of a composite laminate, such as the bottle of solid rocket motor case, all layers are considered to be degraded. Once a lamina being degraded the stiffness as well as strength will be drastically decreased. It becomes a critical area of failure to the whole bottle. And hence the stress/strain field within and around a damage should be accurately evaluated for failure prediction. To investigate the stress/strain field around damages a Hybrid-Numerical method which combines experimental measurement and finite element analysis is used. It is known the stress or strain at the singular point can not be accurately measured by an experimental technique. Nevertheless, if the location is far away from the singular spot, the displacement can be found accurately. Since it reflects the true displacement field locally regardless of the boundary conditions, it is an excellent input data for a finite element analysis to replace the usually assumed boundary conditions. Therefore, the Hybrid-Numerical method is chosen to avoid the difficulty and to take advantage of both experimental technique and finite element analysis. Experimentally, the digital image correlation technique is employed to measure the displacement field. It is done by comparing two digitized images, before and after loading. Numerically, the finite element program, ABAQUS (version 5.2), is used to analyze the stress and strain field. It takes advantage of the high speed and huge memory size of modern supercomputer, CRAY Y-MP, at NASA Marshall Space Flight Center.

Yang, Yii-Ching↗

Engineering Evaluation of International Low Impact Docking System Latch Hooks

The international Low Impact Docking System (iLIDS) provides a structural arrangement that allows for visiting vehicles to dock with the International Space Station (ISS) (Fig 1). The iLIDS docking units are mechanically joined together by a series of active and passive latch hooks. In order to preserve docking capability at the existing Russian docking interfaces, the iLIDS latch hooks are required to conform to the existing Russian design. The latch hooks are classified as being fail-safe. Since the latch hooks are fail-safe, the hooks are not fracture critical and a fatigue based service life assessment will satisfy the structural integrity requirements. Constant amplitude fatigue testing to failure on four sets of active/passive iLIDS latch hooks was performed at load magnitudes of 10, 11, and 12 kips. Failure analysis of the hook fatigue failures identified multi-site fatigue initiation that was effectively centered about the hook mid-plane (consistent with the 3D model results). The fatigue crack initiation distribution implies that the fatigue damage accumulation effectively results in a very low aspect ratio surface crack (which can be simulated as thru-thickness crack). Fatigue damage progression resulted in numerous close proximity fatigue crack initiation sites. It was not possible to determine if fatigue crack coalescence occurs during cyclic loading or as result of the fast fracture response. The presence of multiple fatigue crack initiation sites on different planes will result in the formation of ratchet marks as the cracks coalesce. Once the stable fatigue crack becomes unstable and the fast fracture advances across the remaining ligament and the plane stress condition at a free-surface will result in failure along a 45 deg. shear plane (slant fracture) and the resulting inclined edge is called a shear lip. The hook thickness on the plane of fatigue crack initiation is 0.787". The distance between the shear lips on this plane was on the order of 0.48" and it was effectively centered about the mid-plane of the section. The numerous ratchet marks between the shear lips on the fracture initiation plane are indicative of multiple fatigue initiation sites within this region. The distribution of the fatigue damage about the centerline of the hook is consistent with the analytical results that demonstrate peak stress/strain response at the mid-plane that decreases in the direction of the hook outer surfaces. Scanning electron microscope images of the failed sections detected fatigue crack striations in close proximity to the free surface of the hook radius. These findings were documented at three locations on the fracture surface : 1) adjacent to the left shear lip, 2) adjacent to the right shear lip, and 3) near the centerline of the section. The features of the titanium fracture surface did not allow for a determination of a critical crack size via identification of the region where the fatigue crack propagation became unstable. The fracture based service life projections where benchmarked with strain-life analyses. The strainrange response in the hook radius was defined via the correlated finite element models and the modified method of universal slopes was incorporated to define the strain-life equation for the titanium alloy. The strain-life assessment confirmed that the fracture based projections were reasonable for the loading range of interest. Based upon the analysis and component level fatigue test data a preliminary service life capability for the iLIDS active and passive hooks of 2 lifetimes is projected (includes a scatter factor of 4).

Martinez, J.↗

In situ high energy X-ray diffraction measurement of strain and dislocation density ahead of crack tips grown in hydrogen

The deformation fields near fatigue crack tips grown in hydrogen and in air were measured using high-energy x-ray diffraction. A larger magnitude of elastic strain was observed in the hydrogen case compared to the air case. The magnitude of elastic strain was quantified through an effective crack tip stress intensity factor. The dislocation profile ahead of the crack was probed via x-ray line broadening and electron back-scatter diffraction was used to assess the crack path (intergranular vs. transgranular). Ahead of the crack tip grown in hydrogen, an order of magnitude lower dislocation density, compared to a baseline density far from the crack, was observed. This decrease in dislocation density was not observed in the air case. These differences are discussed in terms of two leading hydrogen embrittlement mechanisms, Hydrogen Enhanced Localized Plasticity (HELP) and Hydrogen Enhanced Decohesion (HEDE). We have observed a decrease in transgranular cohesion (transgranular HEDE), as well as an increase in intergranular fracture. The measurements of dislocation activity support a model of a decrease in intergranular cohesion (intergranular HEDE) which is likely facilitated by the HELP mechanism. Here, this suggests that the increase in fatigue crack growth rate is due to a sum of the two effects of hydrogen, in which the crack grows faster in the transgranular fracture mode and faster due to an increase in a new mode of intergranular fracture. Published by Elsevier Ltd on behalf of Acta Materialia Inc.

36 MATERIALS SCIENCE↗

Yield Behavior of Solution Treated and Aged Ti-6Al-4V

Post yield uniaxial tension-compression tests were run on a solution treated and aged (STA), titanium 6-percent aluminum 4-percent vanadium (Ti-6Al-4V) alloy to determine the yield behavior on load reversal. The material exhibits plastic behavior almost immediately on load reversal implying a strong Bauschinger effect. The resultant stress-strain data was compared to a 1D mechanics model and a finite element model used to design a composite overwrapped pressure vessel (COPV). Although the models and experimental data compare well for the initial loading and unloading in the tensile regime, agreement is lost in the compressive regime due to the Bauschinger effect and the assumption of perfect plasticity. The test data presented here are being used to develop more accurate cyclic hardening constitutive models for future finite element design analysis of COPVs.

Strain↗

Shock-induced chemistry and high strain-rate viscoelastic behavior of a phenolic polymer

We use impact experiments and a finite element model (up to 1.2 GPa), and molecular dynamics simulations (up to 60 GPa), to examine the behavior of a phenolic polymer under shock compression, spanning both nonreactive and reactive regimes. In the nonreactive regime, relaxation following compression at strain rates of ∼105 s−1 can be explained by viscoelasticity observed at ordinary laboratory rates (≲1 s−1) by accounting for the temperature dependence of the phenolic β-transition. Reasonable agreement is found between the measured shock Hugoniot up to 1.2 GPa and molecular dynamics simulation for cross-linked structures of comparable density. We also observed a first-order mechanical transition near 0.36 GPa shock stress and estimated a spall strength of 0.102 GPa and Hugoniot elastic limit of 1–2 GPa. The shock stress is found to vary up to 24% among phenolics made with different resin and/or cure processes. Finally, molecular dynamics simulations are used to identify a reactive regime at shock pressures ≳20 GPa that is characterized by chemically driven, rate-dependent relaxation processes, including dehydrogenation and dehydration reactions that promote the formation of a dense, highly cross-linked carbonaceous solid and the release of light volatiles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Finite element analysis of notch behavior using a state variable constitutive equation

The state variable constitutive equation of Bodner and Partom was used to calculate the load-strain response of Inconel 718 at 649 C in the root of a notch. The constitutive equation was used with the Bodner-Partom evolution equation and with a second evolution equation that was derived from a potential function of the stress and state variable. Data used in determining constants for the constitutive models was from one-dimensional smooth bar tests. The response was calculated for a plane stress condition at the root of the notch with a finite element code using constant strain triangular elements. Results from both evolution equations compared favorably with the observed experimental response. The accuracy and efficiency of the finite element calculations also compared favorably to existing methods.

Dame, L. T.↗

Thermal and structural tests of Rene 41 honeycomb integral-tank concept for future space transportation systems

Two flat 12 by 72 inch Rene 41 honeycomb sandwich panels were tested in a manner to produce combined thermal and mechanical longitudinal stresses that simulated those that would occur in a larger, more complex integral tank and fuselage structure of an earth to orbit vehicle. Elastic strains measured at temperatures below 400 F are compared with calculated values obtained from a linear elastic finite element analysis to verify the analytical model and to establish confidence in the calculated strains. Elastic strain measurement at higher temperatures (between 600 F and 1400 F), where strain measurement is more difficult and less certain, are also compared with calculated strains. Agreement between measured and calculated strains for the lower temperatures is good, but agreement for the higher temperatures is poor because of unreliable strain measurements. Test results indicate that an ascent and entry life cycle of 500 is attainable under high combined thermal and mechanical elastic strains.

Shideler, John L.↗

A microstructure-based modeling approach to predict the mechanical properties of Zr alloy with hydride precipitates

In nuclear reactors, hydrides can form in fuel cladding due to hydrogen absorption in ZIrcaloy and cause embrittlement. This work presents a microstructure-based finite element model to predict the stress-strain response of Zircaloy containing hydrides. Quantitative microstructural details extracted from scanning electron microscopy (SEM) images were used to generate heterogeneous microstructures including the morphology and spatial distribution of hydrides. The constitutive material model for zircaloy in this study is based on crystal plasticity theory which considers the hexagonal close-packed (HCP) atomic structure of Zircaloy material. The hydrides were modeled as brittle material along with a damage model. Hydride formation inside the zircoloy matrix results in residual stress. This phenomenon is also captured in this model. A parametric study has been conducted to understand the effect of volume fraction, orientation, and lamellae thickness of the hydride phase on the mechanical properties of the overall material.

Kulkarni, Shank S.↗

Active interlocking metasurfaces enabled by shape memory alloys

Interlocking metasurfaces (ILMs) are a newly developed joining technology that relies on arrays of interlocking features that transmit force and constrain motion between adjoining bodies in one or more directions. This study explores harnessing the shape memory effect (SME) in Nickel-Titanium shape memory alloys (NiTi SMAs) in structures fabricated using additive manufacturing (AM) to advance the development of active ILMs by creating unit cells that open or close at specific temperatures. The study encompasses designing and fabricating two distinct interlocking array configurations using near-equiatomic NiTi powder and the laser powder bed fusion (L-PBF) AM technique, following a previously developed AM process optimization framework to manufacture defect-free parts. To guide the design process, finite element analysis (FEA) was employed to predict strain values during engage-disengage cycles. The martensitic transformation characteristics of the ILMs were characterized. Thermomechanical testing revealed that the ILMs demonstrate high locking force once engaged, coupled with complete shape recovery and good cyclic stability. Digital image correlation (DIC) was also employed to validate the FEA predictions during the engage-disengage cycles. The results indicate that NiTi SMA-based ILMs can be designed and fabricated into complex shapes using L-PBF. By leveraging the SME, the functionality of an ILM can be improved upon. The combination of computational modeling, additive manufacturing, and thermomechanical and physical property characterization provides a framework for designing future ILMs out of active materials.

Additive manufacturing↗

Stress concentration factors around a circular hole in laminated composites

Stress concentration factors around a circular hole in a composite laminate are determined. The specific case investigated is a four layer (-45/45/45/-45 degs) graphite epoxy laminate. The factors are determined experimentally by means of electrical resistance strain gages, and analytically by using a hybrid finite element analysis.

Ueng, C. E. S.↗

Delamination growth analysis in quasi-isotropic laminates under loads simulating low-velocity impact

A geometrically nonlinear finite-element analysis was developed to calculate the strain energy released by delamination plates during impact loading. Only the first mode of deformation, which is equivalent to static deflection, was treated. Both the impact loading and delamination in the plate were assumed to be axisymmetric. The strain energy release rate in peeling, G sub I, and shear sliding, G sub II, modes were calculated using the fracture mechanics crack closure technique. Energy release rates for various delamination sizes and locations and for various plate configurations and materials were compared. The analysis indicated that shear sliding (G sub II) was the primary mode of delamination growth. The analysis also indicated that the midplane (maximum transverse shear stress plane) delamination was more critical and would grow before any other delamination of the same size near the midplane region. The delamination growth rate was higher (neutrally stable) for a low toughness (brittle) matrix and slower (stable) for high toughness matrix. The energy release rate in the peeling mode, G sub I, for a near-surface delamination can be as high as 0.5G sub II and can contribute significantly to the delamination growth.

Shivakumar, K. N.↗

Ply-level failure analysis of a graphite/epoxy laminate under bearing-bypass loading

A combined experimental and analytical study was conducted to investigate and predict the failure modes of a graphite/epoxy laminate subjected to combined bearing and bypass loading. Tests were conducted in a test machine that allowed the bearing-bypass load ratio to be controlled while a single-fastener coupon was loaded to failure in either tension or compression. Onset and ultimate failure modes and strengths were determined for each test case. The damage-onset modes were studied in detail by sectioning and micrographing the damaged specimens. A two-dimensional, finite-element analysis was conducted to determine lamina strains around the bolt hole. Damage onset consisted of matrix cracks, delamination, and fiber failures. Stiffness loss appeared to be caused by fiber failures rather than by matrix cracking and delamination. An unusual offset-compression mode was observed for compressive bearing-bypass laoding in which the specimen failed across its width along a line offset from the hole. The computed lamina strains in the fiber direction were used in a combined analytical and experimental approach to predict bearing-bypass diagrams for damage onset from a few simple tests.

Naik, R. A.↗

Predicted and observed effects of stacking sequence and delamination size on instability-related delamination growth

A three-dimensional, geometrically nonlinear finite element analysis was used to calculate strain-energy release rate distributions for several laminates which exhibited instability-related delamination growth. Two stacking sequences, (0/90/90/0)6 and (90/0/0/90)6 were considered. The specimens were fabricated with a double layer of 0.013 mm Kapton film located between the fourth and fifth plies, thereby simulating a delamination. Initial delamination sizes were 30, 40, or 60 mm in diameter. The static strain corresponding to delamination growth was determined experimentally. Comparison of the analytical and experimental results indicated that delamination growth was governed by the magnitude of GI.

Whitcomb, John D.↗

Ply-level failure analysis of a graphite/epoxy laminate under bearing-bypass loading

A combined experimental and analytical study was conducted to investigate and predict the failure modes of a graphite/epoxy laminate subjected to combined bearing and bypass loading. Tests were conducted in a test machine that allowed the bearing-bypass load ratio to be controlled while a single-fastener coupon was loaded to failure in either tension or compression. Onset and ultimate failure modes and strengths were determined for each test case. The damage-onset modes were studied in detail by sectioning and micrographing the damaged specimens. A two-dimensional, finite-element analysis was conducted to determine lamina strains around the bolt hole. Damage onset consisted of matrix cracks, delamination, and fiber failures. Stiffness loss appeared to be caused by fiber failures rather than by matrix cracking and delamination. An unusual offset-compression mode was observed for compressive bearing-bypass loading in which the specimen failed across its width along a line offset from the hole. The computed lamina strains in the fiber direction were used in a combined analytical and experimental approach to predict bearing-bypass diagrams for damage onset from a few simple tests.

Naik, R. A.↗

A case study of analysis methods for large deflections of a cantilever beam

A load case study of geometric nonlinear large deflections of a cantilever beam is presented. The bending strain must remain elastic. Closed form solution and finite element methods of analysis are illustrated and compared for three common load cases. A nondimensional nomogram for each case is presented in the summary.

Craig, L. D.↗

Three-dimensional CTOA and constraint effects during stable tearing in a thin-sheet material

A small strain theory, three-dimensional elastic-plastic finite element analysis was used to simulate fracture in thin sheet 2024-T3 aluminum alloy in the T-L orientation. Both straight and tunneled cracks were modeled. The tunneled crack front shapes as a function of applied stress were obtained from the fracture surface of tested specimens. The stable crack growth behavior was measured at the specimen surface as a function of applied stress. The fracture simulation modeled the crack tunneling and extension as a function of applied stress. The results indicated that the global constraint factor, alpha(sub g), initially dropped during stable crack growth. After peak applied stress was achieved, alpha(sub g) began to increase slightly. The effect of crack front shape on alpha(sub g) was small, but the crack front shape did greatly influence the local constraint and through-thickness crack-tip opening angle (CTOA) behavior. The surface values of CTOA for the tunneled crack front model agreed well with experimental measurements, showing the same initial decrease from high values during the initial 3mm of crack growth at the specimen's surface. At the same time, the interior CTOA values increased from low angles. After the initial stable tearing region, the CTOA was constant through the thickness. The three-dimensional analysis appears to confirm the potential of CTOA as a two-dimensional fracture criterion.

Dawicke, D. S.↗