Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stress equilibrium”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Gradient Driven Fluctuations

We have worked with our collaborators at the University of Milan (Professor Marzio Giglio and his group-supported by ASI) to define the science required to measure gradient driven fluctuations in the microgravity environment. Such a study would provide an accurate test of the extent to which the theory of fluctuating hydrodynamics can be used to predict the properties of fluids maintained in a stressed, non-equilibrium state. As mentioned above, the results should also provide direct visual insight into the behavior of a variety of fluid systems containing gradients or interfaces, when placed in the microgravity environment. With support from the current grant, we have identified three key systems for detailed investigation. These three systems are: 1) A single-component fluid to be studied in the presence of a temperature gradient; 2) A mixture of two organic liquids to be studied both in the presence of a temperature gradient, which induces a steady-state concentration gradient, and with the temperature gradient removed, but while the concentration gradient is dying by means of diffusion; 3) Various pairs of liquids undergoing free diffusion, including a proteidbuffer solution and pairs of mixtures having different concentrations, to allow us to vary the differences in fluid properties in a controlled manner.

Cannell, David↗
NASA NTRS2007

A Semi-Analytical Method for Determining the Energy Release Rate of Cracks in Adhesively-Bonded Single-Lap Composite Joints

A semi-analytical method for determining the strain energy release rate due to a prescribed interface crack in an adhesively-bonded, single-lap composite joint subjected to axial tension is presented. The field equations in terms of displacements within the joint are formulated by using first-order shear deformable, laminated plate theory together with kinematic relations and force equilibrium conditions. The stress distributions for the adherends and adhesive are determined after the appropriate boundary and loading conditions are applied and the equations for the field displacements are solved. Based on the adhesive stress distributions, the forces at the crack tip are obtained and the strain energy release rate of the crack is determined by using the virtual crack closure technique (VCCT). Additionally, the test specimen geometry from both the ASTM D3165 and D1002 test standards are utilized during the derivation of the field equations in order to correlate analytical models with future test results. The system of second-order differential field equations is solved to provide the adherend and adhesive stress response using the symbolic computation tool, Maple 9. Finite element analyses using J-integral as well as VCCT were performed to verify the developed analytical model. The finite element analyses were conducted using the commercial finite element analysis software ABAQUS. The results determined using the analytical method correlated well with the results from the finite element analyses.

Yang, Charles↗
NASA NTRS2012

Adhesive-Bonded Composite Joint Analysis with Delaminated Surface Ply Using Strain-Energy Release Rate

This paper presents an analytical model to determine the strain energy release rate due to an interlaminar crack of the surface ply in adhesively bonded composite joints subjected to axial tension. Single-lap shear-joint standard test specimen geometry with thick bondline is followed for model development. The field equations are formulated by using the first-order shear-deformation theory in laminated plates together with kinematics relations and force equilibrium conditions. The stress distributions for the adherends and adhesive are determined after the appropriate boundary and loading conditions are applied and the equations for the field displacements are solved. The system of second-order differential equations is solved to using the symbolic computation tool Maple 9.52 to provide displacements fields. The equivalent forces at the tip of the prescribed interlaminar crack are obtained based on interlaminar stress distributions. The strain energy release rate of the crack is then determined by using the crack closure method. Finite element analyses using the J integral as well as the crack closure method are performed to verify the developed analytical model. It has been shown that the results using the analytical method correlate well with the results from the finite element analyses. An attempt is made to predict the failure loads of the joints based on limited test data from the literature. The effectiveness of the inclusion of bondline thickness is justified when compared with the results obtained from the previous model in which a thin bondline and uniform adhesive stresses through the bondline thickness are assumed.

Chadegani, Alireza↗
DOE OSTI2024

Stress-hybrid virtual element method on six-noded triangular meshes for compressible and nearly-incompressible linear elasticity

In this paper, we present a first-order Stress-Hybrid Virtual Element Method (SH-VEM) on six-noded triangular meshes for linear plane elasticity. Here, we adopt the Hellinger–Reissner variational principle to construct a weak equilibrium condition and a stress based projection operator. In each element, the stress projection operator is expressed in terms of the nodal displacements, which leads to a displacement based formulation. This stress-hybrid approach assumes a globally continuous displacement field while the stress field is discontinuous across each element. The stress field is initially represented by divergence-free tensor polynomials based on Airy stress functions, but we also present a formulation that uses a penalty term to enforce the element equilibrium conditions, referred to as the Penalty Stress-Hybrid Virtual Element Method (PSH-VEM). Numerical results are presented for PSH-VEM and SH-VEM, and we compare their convergence to the composite triangle FEM and B-bar VEM on benchmark problems in linear elasticity. The SH-VEM converges optimally in the L 2 norm of the displacement, energy seminorm, and the L 2 norm of hydrostatic stress. Furthermore, the results reveal that PSH-VEM converges in most cases at a faster rate than the expected optimal rate, but it requires the selection of a suitably chosen penalty parameter.

42 ENGINEERING↗
NASA NTRS1992

The behavior of a liquid drop levitated and drastically flattened by an intense sound field

The deformation and break-up are studied of a liquid drop in levitation through the radiation pressure. Using high-speed photography ripples are observed on the central membrane of the drop, atomization of the membrane by emission of satellite drops from its unstable ripples, and shattering of the drop after upward buckling like an umbrella, or after horizontal expansion like a sheet. These effects are captured on video. The ripples are theorized to be capillary waves generated by the Faraday instability excited by the sound vibration. Atomization occurs whenever the membrane becomes so thin that the vibration is sufficiently intense. The vibration leads to a destabilizing Bernoulli correction in the static pressure. Buckling occurs when an existent equilibrium is unstable to a radial (i.e., tangential) motion of the membrane because of the Bernoulli effect. Besides, the radiation stress at the rim of the drop is a suction stress which can make equilibrium impossible, leading to the horizontal expansion and the subsequent break-up.

Lee, C. P.↗
DOE OSTI2022

Disclination-dislocation based model for grain boundary stress field evolution due to slip transmission history and influence on subsequent dislocation transmission

This work demonstrates how the structure of a grain boundary (GB) and its evolution due to slip transmission history influences subsequent dislocation transmission. First, a model for the evolution of stress fields within grain boundaries that accounts for the effects of coherent dislocation transmission is introduced. Starting with a disclination-based construct of GBs at minimum energy (equilibrium), the model describes the evolution of the GB stress field to a state characteristic of excess energy (non-equilibrium) due to the incorporation of residual Burgers vector content following sequential slip transmission events. Several essential features of this model are verified via molecular dynamics simulations of lattice dislocation absorption. Second, this model is implemented into a discrete dislocation dynamics (DDD) code and simulations are performed to understand the influence of Non-equilibrium GB stress fields, conditioned by the slip transmission history, on subsequent dislocation transmission. DDD simulations reveal that the critical resolved stress necessary for slip propagation can be reduced by with continued absorption of residual dislocation content. Moreover, DDD simulations prove that a comprehensive consideration of both the binding and driving stresses, and the evolution of the transmission configuration is necessary to quantify the influence of the mechanical state of the GB on slip propagation. In general, this work provides important insights into the role of GB structure evolution, conditioned by prior deformation history, on intergranular plasticity.

42 ENGINEERING↗
DOE OSTI2024

Iterative Stress Reconstruction Algorithm to Estimate Three-Dimensional Residual Stress Fields in Manufactured Components

Residual stress (RS) significantly impacts the mechanical performance of components. Measurement of RS often provides incomplete data in terms of components of stress and spatial density. Employing such fields in finite element simulations results in significant modification of the field to achieve equilibrium and compatibility among strains. To overcome this, an iterative stress reconstruction algorithm (ISRA) is developed to estimate 3D RS fields that satisfy equilibrium, are stress component-wise complete, and represent the characterized data sampled. An Al 7075-T651 plate and an additively manufactured (AM) A36 steel wall are considered for RS reconstruction using measurement data from the literature. A maximum variation of ~2.5 MPa in the Al plate, and ~10 MPa in the steel wall are observed between the reconstructed and measured stresses. Furthermore, unknown stress components emerge and reach significant magnitudes (upto ~2.3 MPa in the Al plate and ~45 MPa in the AM wall) during ISRA. Indeed, it is found that minor errors in measurement or data processing are eliminated through the physical requirements during ISRA. Employing a reconstructed RS field is hence not just more accurate given its compatibility, but it additionally corrects for minor errors in measurement. Furthermore, it is found that spatially dense measurement data result in convergence with fewer iterations. Finally, although ISRA yields a nonunique solution dependent on boundary conditions, measurement errors, fitting errors, and mesh density, it accommodates for uncertainties and inaccuracies in measurement, as opposed to failing to reach a physically realistic converged solution.

42 ENGINEERING↗
DOE OSTI2026

Measurements of the inert Hugoniot and observation of mechanical ignition in Ni(V) + Al reactive multilayers via x-ray diffraction

This study investigates the inert Hugoniot response, mechanical ignition, and reaction dynamics of Ni(V)+Al multilayers during longitudinal, laser-driven shock compression experiments. Ni(V)+Al multilayers, known for their self-propagating exothermic reactions, were subjected to longitudinal stresses exceeding 50 GPa using the laser shock facility within the Dynamic Compression Sector (DCS) at the Advanced Photon Source (APS). In situ x-ray diffraction (XRD) revealed that Ni(V) and Al were not in equilibrium during compression, with stress discrepancies attributed to twinning, grain structure effects, and/or dislocation density. However, the measured inert Hugoniot closely matched prior experimental and computational studies, confirming the utility of XRD for measuring the equation of state of thin, complex materials. Additionally, reaction was observed at significantly higher stresses than reported previously using laser-launched flyers. This discrepancy suggests a strong influence of externally imposed shear stress on reaction thresholds, which likely arose from deviations in flyer planarity during past experiments. Full reaction of the multilayer occurred within 40 ns after shock-wave passage, evidenced by complete melting of the constituents. Eulerian hydrocode simulations replicated experimental conditions, providing insights into equilibrium dynamics and experimental artifacts. The results highlight how even small shear forces facilitate ignition in Ni(V)+Al multilayers at lower stresses.

Crystallography↗
NASA NTRS1984

Hybrid Semiloof elements for plates and shells based upon a modified Hu-Washizu principle

Hybrid SemiLoof elements for plates and shells are developed based upon modified Hu-Washizu principle. In the new version of the assumed stress hybrid formulation the equilibrium equations are satisfied through the introduction of internal displacement parameters as Lagrange multipliers. The inversion of the resulting H-matrices is simplified particularly when the stresses are expressed in terms of natural coordinates. A 24-DOF triangular element and a 32-DOF quadrilateral element based on shallow shell theory are derived and evaluated.

Pian, T. H. H.↗
DOE OSTI2024

Initial Testing of an In Situ Load Retention Aging Vessel

A thermal aging vessel instrumented with load cells was fabricated. The primary function of the vessel is to continuously monitor the in situ load retention of up to three compressed polymer coupons undergoing thermally accelerated aging under nitrogen. A secondary function is to enable gas sampling of the vessel headspace during thermal aging. Heating of the vessel is achieved using a custom heater jacket. To improve upon our conventional aging study methods which require periodic interruption of aging to perform load testing in an Instron machine at room temperature, this technology aims to automate/facilitate data acquisition/analysis, improve data quality, and enable uninterrupted compression of the polymer which represents the service condition. As an example case to assess functionality of the in situ vessel, the load retention of a siloxane elastomer material additively manufactured by direct-ink-writing (DIW) was measured at three different isothermal aging temperatures for ~1 month. Initial compression of the coupons while near the aging temperature was achieved by temporarily opening the heated vessel to access the interior chamber and manually tightening four nuts to drive the heated compression plate down onto the heated coupons. Initial testing demonstrated achievement of the primary load retention monitoring function. Unfortunately, the vessel leaked which prevented gas sampling; an active purge was used to maintain a nitrogen atmosphere. Welded or otherwise sealed joints, which could be implemented in a future design, would likely eliminate leak paths. To apply time-temperature superposition (TTS), a technique used to provide long-term prediction of the load retention from short-term isothermal data, the load retention needed to be calculated relative to the load at an estimated “equilibrium” time, after most of the transient viscoelastic physical relaxation occurred. The peak load immediately after compression could not be used as the load retention basis for two reasons: (1) age-related changes must be isolated from non-age-related physical relaxation before applying TTS and (2) the manual mechanism used to compress the specimens at the aging temperature was neither smooth nor repeatable which affected the peak load value. To better understand the effect of the mode of initial compression on the measured load, and possibly better estimate “equilibrium” physical relaxation times, systematic stress relaxation experiments were performed using an Instron machine with a thermal chamber. At a given temperature, the DIW polymer was compressed to a fixed strain in either a stepped or continuous manner at two different rates, then held at that strain for 24 hrs. The results indicated that, at a given temperature, the different stress relaxation curves appeared to converge to the same curve at some “equilibrium” time when the non-age-related physical relaxation was mostly complete. Though this observation suggests that the discontinuous manual compression employed by the vessel is feasible, a compression mechanism that is rapid, smooth, and repeatable would enhance its use.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Static shape and instability of an acoustically levitated liquid drop

Dynamical aspects of a drop drastically flattened by acoustic radiation stress are considered. Its static equilibrium has been studied, starting with a dislike shape and modeling the sound field and the associated radiation stress according to this geometry. It is suggested that, at low viscosity, the ripples are capillary waves generated by the parametric instability excited by the membrane vibration, which is driven by the sound pressure. Atomization occurs whenever the membrane becomes so thin that the vibration is sufficiently intense. Buckling occurs when an existent equilibrium is unstable to a radial oscillation of the membrane because of the Bernoulli effect. The radiation stress at the rim of the flattened drop is also destabilizing and leads to horizontal expansion and subsequent breakup.

Lee, C. P.↗
NASA NTRS2000

Measurement of the Surface Dilatational Viscosity of an Insoluble Surfactant Monolayer at the Air/Water Interface Using a Pendant Drop Apparatus

When a fluid interface with surfactants is at rest, the interfacial stress is isotropic (as given by the equilibrium interfacial tension), and is described by the equation of state which relates the surface tension to the surfactant surface concentration. When surfactants are subjected to shear and dilatational flows, flow induced interaction of the surfactants; can create interfacial stresses apart from the equilibrium surface tension. The simplest relationship between surface strain rate and surface stress is the Boussinesq-Scriven constitutive equation completely characterized by three coefficients: equilibrium interfacial tension, surface shear viscosity, and surface dilatational viscosity Equilibrium interfacial tension and surface shear viscosity measurements are very well established. On the other hand, surface dilatational viscosity measurements are difficult because a flow which change the surface area also changes the surfactant surface concentration creating changes in the equilibrium interfacial tension that must be also taken into account. Surface dilatational viscosity measurements of existing techniques differ by five orders of magnitude and use spatially damped surface waves and rapidly expanding bubbles. In this presentation we introduce a new technique for measuring the surface dilatational viscosity by contracting an aqueous pendant drop attached to a needle tip and having and insoluble surfactant monolayer at the air-water interface. The isotropic total tension on the surface consists of the equilibrium surface tension and the tension due to the dilation. Compression rates are undertaken slow enough so that bulk hydrodynamic stresses are small compared to the surface tension force. Under these conditions we show that the total tension is uniform along the surface and that the Young-Laplace equation governs the drop shape with the equilibrium surface tension replaced by the constant surface isotropic stress. We illustrate this technique using DPPC as the insoluble surfacant monolayer and measured for it a surface dilatational viscosity in the LE phase that is 20 surface poise.

Lorenzo, Jose↗
NASA NTRS1995

Fully-Explicit and Self-Consistent Algebraic Reynolds Stress Models

A fully-explicit, self-consistent algebraic expression for the Reynolds stress, which is the exact solution to the Reynolds stress transport equation in the 'weak equilibrium' limit for two-dimensional mean flows for all linear and some quasi-linear pressure-strain models, is derived. Current explicit algebraic Reynolds stress models derived by employing the 'weak equilibrium' assumption treat the production-to-dissipation (P/epsilon) ratio implicitly, resulting in an effective viscosity that can be singular away from the equilibrium limit. In the present paper, the set of simultaneous algebraic Reynolds stress equations are solved in the full non-linear form and the eddy viscosity is found to be non-singular. Preliminary tests indicate that the model performs adequately, even for three dimensional mean flow cases. Due to the explicit and non-singular nature of the effective viscosity, this model should mitigate many of the difficulties encountered in computing complex turbulent flows with the algebraic Reynolds stress models.

Girimaji, Sharath S.↗
NASA NTRS1979

Analysis of a unidirectional composite containing broken fibers and matrix damage

An analytical solution is developed for the determination of the stresses and displacements in a unidirectional fiber-reinforced composite containing an arbitrary number of broken fibers as well as longitudinal yielding and splitting of the matrix. The solution is developed using a materials-modeling approach which is based on a shear-lag stress transfer mechanism. The equilibrium equation in the axial direction gives a pair of integral equations which are solved numerically. Excellent agreement is shown to exist between the solution and experimental results for notched unidirectional boron/aluminum laminates without splitting. For brittle matrix composites (i.e. epoxy) equally good results are indicated for both matrix yielding and splitting. For yielding without splitting the fracture strength depends on crack length while for large splitting it is crack length independent.

George, J. G.↗
NASA NTRS1980

Analysis of a unidirectional composite containing broken fibers and matrix damage

An analytical solution is developed for the determination of the stresses and displacements in a unidirectional fiber-reinforced composite containing an arbitrary number of broken fibers as well as longitudinal yielding and splitting of the matrix. The solution is developed using a 'materials-modeling' approach which is based on a shear-lag stress transfer mechanism. The equilibrium equation in the axial direction gives a pair of integral equations which are solved numerically. Excellent agreement is shown to exist between the solution and experimental results for notched unidirectional boron/aluminum laminates without splitting. For brittle matrix composites (i.e., epoxy) equally good results are indicated for both matrix yielding and splitting. For yielding without splitting the fracture strength is found to depend on crack length while for large splitting it is crack length independent.

Goree, J. G.↗
NASA NTRS1992

Bending response of thin-walled laminated composite cylinders

Closed-form linear response solutions are presented for the bending response of thin-walled laminated composite cylinders. Displacements and intralaminar stresses are computed using the Donnell theory, and interlaminar stresses are computed using 3D equilibrium equations of elasticity. Laminates with stacking sequencies and an aluminum cylinder are used to illustrate the displacement response, intralaminar, and interlaminar stresses for two length-to-radius ratios of 1 and 5. Data obtained indicate that if only fiber breakage causes ultimate failure, material failure will occur on the compression side of the cylinders, at the ends, due to fiber compression. The interlaminar normal stress is an order of magnitude smaller than the interlaminar shear stresses and it is not likely to initiate failure.

Fuchs, H. P.↗
NASA NTRS1979

Large Eddy Simulation of turbulent shear flows

The conceptual foundation underlying Large Eddy Simulation (LES) is summarized, and the numerical methods developed for simulation of the time-developing turbulent mixing layer and turbulent plane Poiseuille flow are discussed. Computational results show that the average Reynolds stress profile nearly attains the equilibrium shape which balances the downstream mean pressure gradient in the regions away from the walls. In the vicinity of the walls, viscous stresses are shown to be significant; together with the Reynolds stresses, these stresses balance the mean pressure gradient. It is stressed that the subgrid scale contribution to the total Reynolds stress is significant only in the vicinity of the walls. The continued development of LES is urged.

Moin, P.↗