SEARCH · Engineering Papers
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.
Crustal motions in Alaska are computed by calculating the finite element solution for an elastic spherical shell problem. Very long baseline interferometry and geologic data are used and seem to be in rough agreement.
The purpose of this investigation was to determine to what extent ordinary deflection formulas, which neglect shear deformations, are in error when applied to beams of various sections, and to develop reasonably accurate yet comparatively simple formulas which take into account such deformations. A great many tests were made to determine the amount of shear deformation for beams of various sections tested over many different spans. As the span over which the beam is tested is increased the error introduced by neglecting shear deformations becomes less, and the values obtained by substituting measured deflections in the ordinary formulas approach more nearly the modulus of elasticity in tension and compression. For short spans the error is considerable and increases rapidly as the span is reduced. Two formulas were developed for estimating the magnitude of shear deformations, both of which have been verified by tests. The first assumes the parabolic distribution of shear on a cross section of a beam and, starting with a differential volume, the distortion due to shear is determined by the ordinary methods of summarizing the work. The second assumes that the deflections due to shear in any two beams of the same length, height, and moment of inertia, which are similarly loaded, are proportional to the summations of the shear stresses on their respective vertical sections. Both formulas check experimental results very closely when the calculations are made with great refinement.
Herein we discuss a laboratory study of inelastic deformation and localization was conducted to evaluate microseismicity and fracture in fluid-saturated rock. Plane strain compression experiments were performed on oil- and water-saturated Berea sandstone under dry, drained, undrained, and quasi-unjacketed conditions. Parameters associated with the elastic and inelastic material response were accurately measured and used in a constitutive model to predict localization of deformation under the drained and undrained conditions. Acoustic emission (AE) activity was recorded and the onset of inelastic response in fluid-saturated rock, similar to dry rock, coincided with an increase in AE rate. Fracture initiation and propagation were successfully monitored by AE locations. Clustering of AE events and the start of inhomogeneous deformation inferred from displacement measurements are shown to be well correlated with the model prediction of the onset of localized deformation under various rock-fluid interaction conditions.
Probe-induced soft sample damage in atomic force microscopy (AFM), as well as the resulting alteration of local mechanical and electrical properties of the material are explored, specifically comparing contact-mode and intermittent-contact-mode imaging methods. In our experiments, performed on conductive polymer films, induced changes are present in contact-mode imaging while they are negligible or absent in tapping-mode imaging. Here, to understand this result, a viscoelastic parameter extraction is performed, which suggests that permanent sample deformation can readily occur for tip-sample interactions with a duration on the timescale of contact-mode interactions. Using the extracted viscoelastic parameters, a dynamic AFM simulation is conducted, which suggests that the material responds more elastically with reduced or absent sample damage in tapping-mode AFM, due to the higher rate of mechanical deformation and shorter timescales.
The combination of thermal and mechanical loading expected in practice means that constitutive equations of metal matrix composites must be developed which deal with time-independent and time-dependent irreversible deformation. Also, the internal state of composites is extremely complicated which underlines the need to formulate macroscopic constitutive equations with a limited number of state variables which represent the internal state at the micro level. One available method for calculating the macro properties of composites in terms of the distribution and properties of the constituent materials is the method of homogenization whose formulation is based on the periodicity of the substructure of the composite. A homogenization procedure was developed which lends itself to the use of the finite element procedure. The efficiency of these procedures, to determine the macroscopic properties of a composite system from its constituent properties, was demonstrated utilizing an aluminum plate perforated by directionally oriented slits. The selection of this problem is based on the fact that, extensive experimental results exist, the macroscopic response is highly anisotropic, and that the slits provide very high stress gradients which severely test the effectiveness of the computational procedures. Furthermore, both elastic and plastic properties were investigated so that the application to practical systems with inelastic deformation should be able to proceed without difficulty. The effectiveness of the procedures was rigorously checked against experimental results and with the predictions of approximate calculations. Using the computational results it is illustrated how macroscopic constitutive equations can be expressed in forms of the elastic and limit load behavior.
The use of nanoglass (NG) microstructures has been shown as an effective strategy to improve the ductility of rather brittle metallic glass alloys. To evaluate the effects of grain size and gradient design on the mechanical properties of NGs, we perform molecular dynamics simulations of nanoindentation on Cu 64 Zr 36 NGs. We consider samples with uniform 3 and 7 nm grain sizes as well as gradient microstructures with gradient grain sizes varying from 3 to 7 nm. Here, the results show that the deformation mechanism in NGs with small grain sizes is dominated by the activation and evolution of multiple shear transformation zones. Increasing grain sizes enhances the elastic modulus and hardness at the cost of reduced abrasion resistance. Though the average grain size at the indentation surface plays a crucial role in the deformation behavior of gradient NG models, their strain localization and plastic deformation states are affected by the grain sizes far from the indenter. Desired combinations of mechanical properties can be realized by different NG microstructure designs.
Potential applications involving high-temperature shape memory alloys have been growing in recent years. Even in those cases where promising new alloys have been identified, the knowledge base for such materials contains gaps crucial to their maturation and implementation in actuator and other applications. We begin to address this issue by characterizing the mechanical behavior of a Ni19.5Pd30Ti50.5 high-temperature shape memory alloy in both uniaxial tension and compression at various temperatures. Differences in the isothermal uniaxial deformation behavior were most notable at test temperatures below the martensite finish temperature. The elastic modulus of the material was very dependent on strain level; therefore, dynamic Young#s Modulus was determined as a function of temperature by an impulse excitation technique. More importantly, the performance of a thermally activated actuator material is dependent on the work output of the alloy. Consequently, the strain-temperature response of the Ni19.5Pd30Ti50.5 alloy under various loads was determined in both tension and compression and the specific work output calculated and compared in both loading conditions. It was found that the transformation strain and thus, the specific work output were similar regardless of the loading condition. Also, in both tension and compression, the strain-temperature loops determined under constant load conditions did not close due to the fact that the transformation strain during cooling was always larger than the transformation strain during heating. This was apparently the result of permanent plastic deformation of the martensite phase with each cycle. Consequently, before this alloy can be used under cyclic actuation conditions, modification of the microstructure or composition would be required to increase the resistance of the alloy to plastic deformation by slip.
Ambidirectionality, which is the ability of structural elements to move beyond a reference state in two opposite directions, is common in nature. However, conventional soft materials are typically limited to a single, unidirectional deformation unless complex hybrid constructs are used. We exploited the combination of mesogen self-assembly, polymer chain elasticity, and polymerization-induced stress to design liquid crystalline elastomers that exhibit two mesophases: chevron smectic C (cSmC) and smectic A (SmA). Inducing the cSmC-SmA–isotropic phase transition led to an unusual inversion of the strain field in the microstructure, resulting in opposite deformation modes (e.g., consecutive shrinkage or expansion and right-handed or left-handed twisting and tilting in opposite directions) and high-frequency nonmonotonic oscillations. Further, this ambidirectional movement is scalable and can be used to generate Gaussian transformations at the macroscale.
Abstract Many energy systems demand heat transfer at high temperatures to keep up with high demand for power, so high‐temperature material that can perform and last under these harsh conditions is needed for heat exchangers. The engineering requirements for these high‐temperature heat exchanger material call for high thermal conductivity, high resistance to fracture, high resistance to creep deformation, environmental stability in environments associated with the application, and high modulus of elasticity while maintaining low cost to make and maintain. Naturally, ceramics are a good solution for this endeavor. In the past, high‐temperature heat exchangers made from ceramics have been used. We provide examples of ceramics in relevant heat exchange applications and provide motivation where additive manufacturing (AM) can improve efficiency. AM for the relevant material is under development, and we provide insight on the AM of ceramic materials and examples of AM heat exchangers keeping cost in mind. The motivation of the review paper is to provide a framework for material and manufacturing selection for high‐temperature heat exchangers for AM to keep up with the demand for better efficiency, better material, better manufacturing, and cost moving forward with AM technology in high‐temperature ceramic heat exchangers.
High-entropy alloys (HEAs) are materials that consist of equimolar or near-equimolar multiple principal components but tend to form single phases, which is a new research topic in the field of metallurgy, and have attracted extensive attention in the past decade. The HEAs families contain the face-centered-cubic (fcc), body-centered-cubic (bcc), and hexagonal-close-packed (hcp)-structured HEAs. On the one hand, mechanical properties, e.g., hardness, strength, ductility, fatigue, and elastic moduli, are essential for practical applications of HEAs. Scientists have explored in this direction since the advent of HEAs. On the other hand, the pursuit of high strength and good plasticity is the critical research issue of materials. Hence, strengthening of HEAs is a crucial issue.In this chapter, we reviewed the recent work on the room-temperature elastic properties and mechanical behavior of HEAs, including the mechanisms behind the plastic deformation of HEAs at both low and high temperatures. Furthermore, the present work examined the strengthening strategies of HEAs, e.g., strain hardening, grain-boundary strengthening, solid-solution strengthening, and particle strengthening. The fatigue, creep, and fracture properties were briefly introduced. Lastly, the future scientific issues and challenges of HEAs were discussed.
This work examines the frame-invariance (and the lack thereof) exhibited in simulated anisotropic elasto-plastic responses generated from supervised machine learning of classical multi-layer and informed-graph-based neural networks, and proposes different remedies to fix this drawback. The inherent hierarchical relations among physical quantities and state variables in an elasto-plasticity model are first represented as informed, directed graphs, where three variations of the graph are tested. While feed-forward neural networks are used to train path-independent constitutive relations (e.g., elasticity), recurrent neural networks are used to replicate responses that depends on the deformation history, i.e. or path dependent. In dealing with the objectivity deficiency, we use the spectral form to represent tensors and, subsequently, three metrics, the Euclidean distance between the Euler Angles, the distance from the identity matrix, and geodesic on the unit sphere in Lie algebra, can be employed to constitute objective functions for the supervised machine learning. In this, the aim is to minimize the measured distance between the true and the predicted 3D rotation entities. Following this, we conduct numerical experiments on how these metrics, which are theoretically equivalent, may lead to differences in the efficiency of the supervised machine learning as well as the accuracy and robustness of the resultant models. Neural network models trained with tensors represented in component form for a given Cartesian coordinate system are used as a benchmark. Our numerical tests show that, even given the same amount of information and data, the quality of the anisotropic elasto-plasticity model is highly sensitive to the way tensors are represented and measured. The results reveal that using a loss function based on geodesic on the unit sphere in Lie algebra together with an informed, directed graph yield significantly more accurate rotation prediction than the other tested approaches.
𝛾 decays of the isovector giant dipole resonance (IVGDR) of the deformed nucleus 154 Sm were measured using 2$^{+}_{1}$-Smekal-Raman and elastic scattering of linearly polarized, quasimonochromatic photon beams. The two scattering processes were disentangled through their distinct angular distributions. Their branching ratio and cross sections were determined at six excitation energies covering the 154 Sm IVGDR. Both agree with the predictions of the geometrical model for the IVGDR and confirm 𝛾 decay as an observable sensitive to the structure of the resonance. Consequently, the data place strong constraints on the nuclear shape, including the degree of triaxiality. The derived 154 Sm shape parameters 𝛽 = 0.2925(25) and 𝛾 = 5.0(15)° agree well with other measurements and recent Monte Carlo shell-model calculations.
Programs for the analyses of the free or forced, undamped vibrations of one or two elastically-coupled lumped parameter teams are presented. Bearing nonlinearities, casing and rotor distributed mass and elasticity, rotor imbalance, forcing functions, gyroscopic moments, rotary inertia, and shear and flexural deformations are all included in the system dynamics analysis. All bearings have nonlinear load displacement characteristics, the solution is achieved by iteration. Rotor imbalances allowed by such considerations as pilot tolerances and runouts as well as bearing clearances (allowing concail or cylindrical whirl) determine the forcing function magnitudes. The computer programs first obtain a solution wherein the bearings are treated as linear springs of given spring rates. Then, based upon the computed bearing reactions, new spring rates are predicted and another solution of the modified system is made. The iteration is continued until the changes to bearing spring rates and bearing reactions become negligibly small.
The nonlinear curvature expressions for a twisted rotor blade or a beam undergoing transverse bending in two planes, torsion, and extension were developed. The curvature expressions were obtained using simple geometric considerations. The expressions were first developed in a general manner using the geometrical nonlinear theory of elasticity. These general nonlinear expressions were then systematically reduced to four levels of approximation by imposing various simplifying assumptions, and in each of these levels the second degree nonlinear expressions were given. The assumptions were carefully stated and their implications with respect to the nonlinear theory of elasticity as applied to beams were pointed out. The transformation matrices between the deformed and undeformed blade-fixed coordinates, which were needed in the development of the curvature expressions, were also given for three of the levels of approximation. The present curvature expressions and transformation matrices were compared with corresponding expressions existing in the literature.
The development of high-precision instrumentation for monitoring benchmark stability and evaluating coseismic strain and tilt signals is reviewed. Laser strainmeter and tilt observations are presented. Examples of coseismic deformation in several geographic locations are given. Evidence suggests that the Earth undergoes elastic response to abrupt faulting.
The research work performed during the past year on finite element implementation and computational techniques pertaining to high temperature composites is outlined. In the present research, two main issues are addressed: efficient geometric modeling of composite structures and expedient numerical integration techniques dealing with constitutive rate equations. In the first issue, mixed finite elements for modeling laminated plates and shells were examined in terms of numerical accuracy, locking property and computational efficiency. Element applications include (currently available) linearly elastic analysis and future extension to material nonlinearity for damage predictions and large deformations. On the material level, various integration methods to integrate nonlinear constitutive rate equations for finite element implementation were studied. These include explicit, implicit and automatic subincrementing schemes. In all cases, examples are included to illustrate the numerical characteristics of various methods that were considered.
The thermal stresses on a cryogenic storage tank contribute strongly to the state of stress of the tank material and its ability to withstand operational stresses. These thermal stresses also affect the growth of any surface damage that might occur in the tank walls. These stresses are particularly of concern during the initial cooldown period for a new tank placed into service, and during any subsequent thermal cycles. A previous preliminary elastic analysis showed that the thermal stress on the inner wall would reach approximately 1,000MPa (145,000 psi). This stress far exceeds the ASTM specified room temperature values for both yield (170MPa) and ultimate (485 MPa) strength for 304L stainless steel. The present analysis determines the thermal stresses using an elastic-plastic model. The commercial software application ANSYS was used to determine the transient spatial temperature profile and the associated spatial thermal stress profiles in a segment of a thick-walled vessel during a typical cooldown process. A strictly elastic analysis using standard material properties for 304L stainless steel showed that the maximum thermal stress on the inner and outer walls was approximately 960 MPa (tensile) and - 270 MPa (compressive) respectively. These values occurred early in the cooldown process, but at different times, An elastic-plastic analysis showed significantly reducing stress, as expected due to the plastic deformation of the material. The maximum stress for the inner wall was approximately 225 MPa (tensile), while the maximum stress for the outer wall was approximately - 130 MPa (compressive).
Maps of post-seismic surface displacement after the 1992, Landers, California earthquake, generated by interferometric processing of ERS-1 Synthetic Aperture Radar (SAR) images, reveal effects of various deformation processes near the 1992 surface rupture.