Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Strain calculation”

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 361 records · Page 20

Quasi-Static Viscoelastic Finite Element Model of an Aircraft Tire

An elastic large displacement thick-shell mixed finite element is modified to allow for the calculation of viscoelastic stresses. Internal strain variables are introduced at the element's stress nodes and are employed to construct a viscous material model. First order ordinary differential equations relate the internal strain variables to the corresponding elastic strains at the stress nodes. The viscous stresses are computed from the internal strain variables using viscous moduli which are a fraction of the elastic moduli. The energy dissipated by the action of the viscous stresses is included in the mixed variational functional. The nonlinear quasi-static viscous equilibrium equations are then obtained. Previously developed Taylor expansions of the nonlinear elastic equilibrium equations are modified to include the viscous terms. A predictor-corrector time marching solution algorithm is employed to solve the algebraic-differential equations. The viscous shell element is employed to computationally simulate a stair-step loading and unloading of an aircraft tire in contact with a frictionless surface.

Johnson, Arthur R.↗

Modified Mode-I Cracked Sandwich Beam (CSB) Fracture Test

Five composite sandwich panels were fabricated using vacuum assisted resin transfer molding (VARTM). Four of these panels had E-glass/vinylester facesheets and one had carbon/epoxy facesheets. The sandwich panels had different density PVC foam cores. The four E-glass panels had core densities of 80, 100, 130, 200 kg/cu m. The sandwich with carbon/epoxy 3 facesheets had a core with density of 100 kg/cu m. Fracture tests were conducted using a modified Cracked Sandwich Beam (CSB) test configuration. Load displacement curves were obtained for loading and unloading of the specimens during crack growth. Various increments of crack growth were monitored. Critical Strain Energy Release Rates (SERR) were determined from the tests using the area method. The critical values of SERR can be considered the fracture toughness of the sandwich material. The fracture toughness ranged 367 J/sq m to 1350 J/sq m over the range of core densities. These results are compared to the Mode-I fracture toughness of the PVC foam core materials and values obtained for foam-cored sandwiches using the TSD specimen. Finite-element analyses (FEA) were performed for the test configuration and Strain Energy Release Rates were calculated using the Virtual Crack Closure Technique (VCCT). The SERR values determined from the FEA were scaled to the fracture loads, or critical loads, obtained from the modified CSB tests. These critical loads were in close agreement with the test values.

Smith, S. A.↗

A Viscoelastic Hybrid Shell Finite Element

An elastic large displacement thick-shell hybrid finite element is modified to allow for the calculation of viscoelastic stresses. Internal strain variables are introduced at he element's stress nodes and are employed to construct a viscous material model. First order ordinary differential equations relate the internal strain variables to the corresponding elastic strains at the stress nodes. The viscous stresses are computed from the internal strain variables using viscous moduli which are a fraction of the elastic moduli. The energy dissipated by the action of the viscous stresses in included in the mixed variational functional. Nonlinear quasi-static viscous equilibrium equations are then obtained. Previously developed Taylor expansions of the equilibrium equations are modified to include the viscous terms. A predictor-corrector time marching solution algorithm is employed to solve the algebraic-differential equations. The viscous shell element is employed to numerically simulate a stair-step loading and unloading of an aircraft tire in contact with a frictionless surface.

Johnson, Arthur↗

Skin-Stiffener Debond Prediction Based on Computational Fracture Analysis

Interlaminar fracture mechanics has proven useful for characterizing the onset of delaminations in composites and has been used with limited success primarily to investigate onset in fracture toughness specimens and laboratory size coupon type specimens. Future acceptance of the methodology by industry and certification authorities however, requires the successful demonstration of the methodology on structural level. For this purpose a panel was selected that is reinforced with stringers. Shear loading causes the panel to buckle and the resulting out-of-plane deformations initiate skin/stringer separation at the location of an embedded defect. For finite element analysis, the panel and surrounding load fixture were modeled with shell elements. A small section of the stringer foot and the panel in the vicinity of the embedded defect were modeled with a local 3D solid model. Across the width of the stringer foot the mixed-mode strain energy release rates were calculated using the virtual crack closure technique. A failure index was calculated by correlating the results with the mixed-mode failure criterion of the graphite/epoxy material. For small applied loads the failure index is well below one across the entire width. With increasing load the failure index approaches one first near the edge of the stringer foot from which delamination is expected to grow. With increasing delamination lengths the buckling pattern of the panel changes and the failure index increases which suggests that rapid delamination growth from the initial defect is to be expected.

Krueger, Ronald↗

Analysis of Composite Panel-Stiffener Debonding Using a Shell/3D Modeling Technique

Interlaminar fracture mechanics has proven useful for characterizing the onset of delaminations in composites and has been used successfully primarily to investigate onset in fracture toughness specimens and laboratory size coupon type specimens. Future acceptance of the methodology by industry and certification authorities, however, requires the successful demonstration of the methodology on the structural level. For this purpose, a panel was selected that is reinforced with stiffeners. Shear loading causes the panel to buckle, and the resulting out-of-plane deformations initiate skin/stiffener separation at the location of an embedded defect. A small section of the stiffener foot, web and noodle as well as the panel skin in the vicinity of the delamination front were modeled with a local 3D solid model. Across the width of the stiffener foot, the mixedmode strain energy release rates were calculated using the virtual crack closure technique. A failure index was calculated by correlating the results with a mixed-mode failure criterion of the graphite/epoxy material. Computed failure indices were compared to corresponding results where the entire web was modeled with shell elements and only a small section of the stiffener foot and panel were modeled locally with solid elements. Including the stiffener web in the local 3D solid model increased the computed failure index. Further including the noodle and transition radius in the local 3D solid model changed the local distribution across the width. The magnitude of the failure index decreased with increasing transition radius and noodle area. For the transition radii modeled, the material properties used for the noodle area had a negligible effect on the results. The results of this study are intended to be used as a guide for conducting finite element and fracture mechanics analyses of delamination and debonding in complex structures such as integrally stiffened panels.

Krueger, Ronald↗

A Simple Test Method for Large Deformation Bending of Thin High Strain Composite Flexures

A simple test method for large deformation bending of thin composite laminates is investigated using image processing and full-field strain measurements. The assumptions and kinematic equations that represent the test are used to calculate numerically the laminate bending stiffness and strength as well as the curvature and strains at failure. In order to validate the test methodology, a comparison is performed between analytical model predictions and empirical data in terms of computed surface strains versus digital image correlation data and calculated rotation angles of the fixture arms throughout the test versus measured ones. The new test method is then used to calculate the bending stiffness in the D11 and D22 directions as well as failure strains for various thin-ply laminates of interest. These parameters are ultimately compared with predicted values using micromechanics and classical lamination theory analysis. In general, bending stiffness and strain test results and predictions for 0 degree orientation coupons have a maximum difference of 10% and 35%.

Fernandez, Juan M.↗

Real-Time Modulation of Hydrogen Evolution Activity of Graphene Electrodes Using Mechanical Strain

This paper reports the effect of mechanically applied elastic strain on the hydrogen evolution reaction (HER) activity of graphene under acidic conditions. Here, applied tensile strain of 0.2% on a graphene electrode is shown to lead to a 1–3% increase of HER current. Tensile strain increases HER activity while compressive strain decreases it. Density functional theory (DFT) calculations using a periodic graphene slab model show increasing in the adsorption energy of H atom with growing tensile strain, leading to a corresponding enhancement of current density in HER, similar to that observed experimentally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Full Flight Envelope Direct Thrust Measurement on a Supersonic Aircraft

Direct thrust measurement using strain gages offers advantages over analytically-based thrust calculation methods. For flight test applications, the direct measurement method typically uses a simpler sensor arrangement and minimal data processing compared to analytical techniques, which normally require costly engine modeling and multisensor arrangements throughout the engine. Conversely, direct thrust measurement has historically produced less than desirable accuracy because of difficulty in mounting and calibrating the strain gages and the inability to account for secondary forces that influence the thrust reading at the engine mounts. Consequently, the strain-gage technique has normally been used for simple engine arrangements and primarily in the subsonic speed range. This paper presents the results of a strain gage-based direct thrust-measurement technique developed by the NASA Dryden Flight Research Center and successfully applied to the full flight envelope of an F-15 aircraft powered by two F100-PW-229 turbofan engines. Measurements have been obtained at quasi-steady-state operating conditions at maximum non-augmented and maximum augmented power throughout the altitude range of the vehicle and to a maximum speed of Mach 2.0 and are compared against results from two analytically-based thrust calculation methods. The strain-gage installation and calibration processes are also described.

Conners, Timothy R.↗

Charge Separation in Monolayer WSe 2 by Strain Engineering: Implications for Strain-Induced Diode Action

Strain-engineering band structure in transition-metal dichalcogenides (TMDC) is a promising avenue toward capabilities in optoelectronics. For example, controlling the flow of optically generated quasiparticles can be achieved by a localized strain field which reduces the bandgap and generates an energy-band gradient that funnels neutral excitons to the strain apex. It would be even more advantageous to mimic a diode’s internal field, where both conduction and valence bands bend in the same direction, to separate electrons and holes. This can be achieved if the strain in the TMDC layer lowers both the conduction band minimum as well as the valence band maximum during strain-induced band narrowing. Here, we have used density functional theory (DFT) calculations of monolayer WSe 2 electronic structure under biaxial strain to show that WSe 2 has this property. In this work, to test the band bending experimentally, we combined localized strain with electrostatic doping to follow photoluminescence from excitons and positive or negative trions. In unstrained WSe 2 , both positive and negative trion emissions dominate over excitons away from charge neutrality. In contrast, for strained areas, negative trions accumulate, while positive trion emission is near zero away from charge neutrality, indicating a lack of holes. Hence, strain bends both conduction and valence bands down, similarly to the band bending in a PN-diode depletion region, providing an opportunity to separate electrons and holes via localized strain.

2D materials↗

Characterization of oriented cracks with differential strain analysis

Linear strain of a rock sample as a function of hydrostatic pressure can be measured with a precision of 2 millionths. Such high-precision data for three orthogonal directions allow calculation of the distribution function for the porosity due to cracks' closing completely at a given pressure. Such data for at least six independent directions yield the zero-pressure strain tensor due to cracks' closing completely at a given pressure. The principal values and axes of this tensor distribution function provide information about the orientation of cracks as a function of closure pressure. The mathematical basis for the technique is developed, and the technique is illustrated with differential strain data for two samples, the Westerly (Rhode Island) granite and the Twin Sisters (Washington) dunite. Strain-tensor calculations reveal that each of these samples has a different type of anisotropic crack distribution.

Siegfried, R.↗

Elastic-plastic finite-element analyses of thermally cycled double-edge wedge specimens

Elastic-plastic stress-strain analyses were performed for double-edge wedge specimens subjected to thermal cycling in fluidized beds at 316 and 1088 C. Four cases involving different nickel-base alloys (IN 100, Mar M-200, NASA TAZ-8A, and Rene 80) were analyzed by using the MARC nonlinear, finite element computer program. Elastic solutions from MARC showed good agreement with previously reported solutions obtained by using the NASTRAN and ISO3DQ computer programs. Equivalent total strain ranges at the critical locations calculated by elastic analyses agreed within 3 percent with those calculated from elastic-plastic analyses. The elastic analyses always resulted in compressive mean stresses at the critical locations. However, elastic-plastic analyses showed tensile mean stresses for two of the four alloys and an increase in the compressive mean stress for the highest plastic strain case.

Kaufman, A.↗

Strain-energy-release rate analysis of the end-notched flexure specimen using the finite-element method

Two-dimensional finite-element analysis of the end-notched flexure specimen was performed using 8-node isoparametric, parabolic elements to evaluate compliance and mode II strain energy release rates, G sub II. The G sub II values were computed using two different techniques: the virtural crack-closure technique (VCCT) and the rate of change of compliance with crack length (compliance derivative method). The analysis was performed for various crack-length-to-semi-span (a/L) ratios ranging from 0.2 to 0.9. Three material systems representing a wide range of material properties were analyzed. The compliance and strain energy release rates of the specimen calculated with the present finite-element analysis agree very well with beam theory equations including transverse shear. The G sub II values calculated using the compliance derivative method compared extremely well with those calculated using the VCCT. The G sub II values obtained by the compliance derivative method using the top or bottom beam deflections agreed closely with each other. The strain energy release rates from a plane-stress analysis were higher than the plane-strain values by only a small percentage, indicating that either assumption may be used in the analysis. The G sub II values for one material system calculated from the finite-element analysis agreed with one solution in the literature and disagreed with the other solution in the literature.

Salpekar, S. A.↗

Strain energy release rate analysis of the end-notched flexure specimen using the finite-element method

Two-dimensional finite-element analysis of the end-notched flexure specimen was performed using 8-node isoparametric, parabolic elements to evaluate compliance and mode II strain energy release rates, G sub II. The G sub II values were computed using two different techniques: the virtual crack-closure technique (VCCT) and the rate of change of compliance with crack length (compliance derivative method). The analysis was performed for various crack-length-to-semi-span (a/L) ratios ranging from 0.2 to 0.9. Three material systems representing a wide range of material properties were analyzed. The compliance and strain energy release rates of the specimen calculated with the present finite-element analysis agree very well with beam theory equations including transverse shear. The G sub II values calculated using the compliance derivative method compared extremely well with those calculated using the VCCT. The G sub II values obtained by the compliance derivative method using the top or bottom beam deflections agreed closely with each other. The strain energy release rates from a plane-stress analysis were higher than the plane-strain values by only a small percentage, indicating that either assumption may be used in the analysis. The G sub II values for one material system calculated from the finte-element analysis agreed with one solution in the literature and disagreed with the other solution in the literature.

Salpekar, S. A.↗

Analysis of the bonding in tetrahedrane and phosphorus-substituted tetrahedranes

The bonding structures of tetrahedrane, phosphatetrahedrane, diphosphatetrahedrane and triphosphatetrahedrane are studied by employing an intrinsic quasi-atomic orbital analysis. Ethane, cyclopropane and tetrahedral P 4 are employed as reference systems. The orbital analysis is paired with the computation of strain energies via isodesmic reactions. The results show that the increase in geometric strain upon transition from ethane to cyclopropane to tetrahedrane weakens the CC bonds, despite leading to shorter C–C interatomic distances. With the increase in strain, the orbitals centered on C and involved in the bonding of the cage structure are observed to have elevated p-character, and the orbital structure of C deviates from sp 3 hybridization. The systematic substitution of CH groups by P atoms in the cage structure of tetrahedrane leads to stronger CC bonds, larger angles in the cage structures of the resulting phosphatetrahedranes, lower p-character in the orbitals involved in the bonding of the cages, and lower strain energies. It is found that P is more amenable to strained molecular arrangements than is C, and that the propensity of a given atom to hybridize s and p functions, or the lack thereof, has implications in the stability of molecules with strained geometries. The combination of the calculations presented here with the existing literature provides insight into the apparent propensity of tetrahedrane and P 4 to ‘break’ their tetrahedral structures. Trends in the bonding interactions, such as bond strengths, s- and p-orbital characters and charge transfer are identified and related to the strain energy observed in each of the analyzed systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electronic structure and magnetism of pristine, defected, and strained $\mathrm{Ti_2N}$ $\mathrm{MX}$ene

Here, from first principles electronic structure calculations, we unravel the evolution of structural, electronic, and magnetic properties of pristine, defected, and strained titanium nitride MXene with different functional groups (-F, -O, -H, and -OH). The formation and cohesive energies reveal their chemical stability. The dynamical stability of Ti 2 N mono-layer is also confirmed by phonon calculations. The MAX phase and defect free functionalized MXenes are metallic except for oxygen terminated (Ti 2 NO 2 ) one which is 100% spin polarized half-metallic ferromagnet. The spin–orbit coupling significantly influences the bare MXene (Ti 2 N) to exhibit Dirac topology and band inversion near the high symmetry directions. The strain effect sways the Fermi level thereby shifting it towards lower energy state under compression and towards higher energy state under tensile strain in Ti 2 NH 2 . The Ti 2 NO 2 exhibits exotic electronic structure not only in pristine but also in strained and defected structures. Its half-metallic nature changes to semi-metallic under 1% compression and it is completely destroyed under 2% compression. In single vacancy defect, its band structure remarkably transforms from half-metallic to semi-conducting with large band gap in 12.5% Ti, weakly semi-conducting in 5.5% Ti, and semi-metallic in 12.5% O. The 25% N defect changes its half-metallic characteristic to metallic. Further, the 12.5% Co substitution preserves its half-metallic character, whereas Mn substitution allows it to convert half-metallic characteristic into weak semi-metallic characteristic preserving ferromagnetism. However, Cr substitution converts half-metallic ferromagnetic state to half-metallic anti-ferromagnetic state. The understanding made here on collective structural stability, and electronic band structure, and magnetic phenomena in novel 2D Ti 2 N derived MXenes open up their possibility in designing them for synthesis.

36 MATERIALS SCIENCE↗

Calculation of residual principal stresses in CVD boron on carbon filaments

A three-dimensional finite element model of the chemical vapor deposition (CVD) of boron on a carbon substrate (B/C) is developed. The model includes an expansion of the boron after deposition due to atomic rearrangement and includes creep of the boron and carbon. Curves are presented to show how the principal residual stresses and the filament elongation vary as the parameters defining deposition strain and creep are varied. The calculated results are compared with experimental axial residual stress and elongation measurements made on B/C filaments. This comparison requires that for good agreement between calculated and experimental results, the deposited boron must continue to expand after deposition, and that the build-up of residual stresses is limited by significant boron and carbon creep rates.

Behrendt, D. R.↗

Synergistic effect of microstructure and defects on the initiation of fatigue cracks in additively manufactured Inconel 718

Fatigue cracks in additively manufactured (AMed) Inconel 718 (IN-718) in machined surface condition often initiate from persistent slip bands (PSBs) unlike other popular AM alloys such as 17-4 PH stainless steel or Ti-6Al-4V, where fatigue crack initiation is exclusively from volumetric defects; therefore, a competition between PSB- vs. defect- mediated crack initiation clearly exists. To shed light on the factors governing the competition, this study investigates the characteristics of cyclic strain localization, PSB formation, and crack initiation via crystal plasticity (CP) modeling of cyclic loading on polycrystalline aggregates which are then validated by experiments. In this work, a physics-based, free slip distance (FSD) dependent slip strength evolution law is proposed, which is shown to enable the CP model to simulate the heterogeneous strain distribution in IN-718. Implementing a crack initiation criterion based on strain contrast, the locations and lives for crack initiation can also be calculated. It is shown that both FSD and resolve shear stress influence the strain localization and crack initiation behaviors. The distribution of a localization parameter calculated based on the multiplication of FSD and Schmid factor within a grain is found to correlate well with the locations of PSBs. The maximum values of the localization parameters within a microstructure are shown to correlate well with the experimentally obtained crack initiation lives. The presence of volumetric defects in IN-718 generally do not impact strain localization behavior unless their size is large compared to the grain size.

42 ENGINEERING↗

Thermal elastic deformations of the planet Mercury

The variation in solar heating due to the resonance rotation of Mercury produces periodic elastic deformations on the surface of the planet. The thermal stress and strain fields under Mercury's surface are calculated after certain simplifications. It is shown that deformations penetrate to a greater depth than the variation of solar heating, and that the thermal strain on the surface of the planet pulsates with an amplitude of 0.004 and a period of 176 days.

Liu, H.↗