Search NASA⌕ Search

SEARCH · Search NASA

Results for “shear displacement”

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 271 records · Page 15

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗

Jet transport performance in thunderstorm wind shear conditions

Several hours of three dimensional wind data were collected in the thunderstorm approach-to-landing environment, using an instrumented Queen Air airplane. These data were used as input to a numerical simulation of aircraft response, concentrating on fixed-stick assumptions, while the aircraft simulated an instrument landing systems approach. Output included airspeed, vertical displacement, pitch angle, and a special approach deterioration parameter. Theory and the results of approximately 1000 simulations indicated that about 20 percent of the cases contained serious wind shear conditions capable of causing a critical deterioration of the approach. In particular, the presence of high energy at the airplane's phugoid frequency was found to have a deleterious effect on approach quality. Oscillations of the horizontal wind at the phugoid frequency were found to have a more serious effect than vertical wind. A simulation of Eastern flight 66, which crashed at JFK in 1975, served to illustrate the points of the research. A concept of a real-time wind shear detector was outlined utilizing these results.

Mccarthy, J.↗

Shear deformation plate continua of large double layered space structures

A simple method is presented to model large rigid-jointed lattice structures as continuous elastic media with couple stresses using energy equivalence. In the analysis, the transition from the discrete system to the continuous media is achieved by expanding the displacements and the rotations of the nodal points in a Taylor series about a suitable chosen origin. The strain energy of the continuous media with couple stresses is then specialized to obtain shear deformation plate continua. Equivalent continua for single layered grids, double layered grids, and three-dimensional lattices are then obtained.

Hefzy, Mohamed Samir↗

Vibration computer programs E13101, E13102, E13104, and E13112 and application to the NERVA program. Project 187: Methodology documentation

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.

Mironenko, G.↗

An efficient quadrilateral element for plate bending analysis

A simple, shear flexible, quadrilateral plate element is developed based on the Hellinger/Reissner mixed variational principle with independently assumed displacement and stress fields. The crucial point of the selection of appropriate stress parameters is emphasized in the formulation. For this purpose, a set of guidelines is formulated based on the follwoing considerations: (1) suppression of all kinematic deformation modes, (2) the element has a favorable value for the constraint index in the thin plate limit, and (3) element properties are frame-invariant. For computer implementation the components of the element stiffness matrix are evaluated analytically using the symbolic manipulation package MACSYMA. The effectiveness and practical usefulness of the proposed element are demonstrated by the numerical results of a variety of problems involving thin and moderately thick plates under different loading and support conditions.

Saleeb, A. F.↗

Mixed finite element models for free vibrations of thin-walled beams

Simple, mixed finite element models are developed for the free vibration analysis of curved thin-walled beams with arbitrary open cross section. The analytical formulation is based on a Vlasov's type thin-walled beam theory with the effects of flexural-torsional coupling, transverse shear deformation and rotary inertia included. The fundamental unknowns consist of seven internal forces and seven generalized displacements of the beam. The element characteristic arrays are obtained by using a perturbed Lagrangian-mixed variational principle. Only C(sup o) continuity is required for the generalized displacements. The internal forces and the Lagrange multiplier are allowed to be discontinuous at interelement boundaries. Numerical results are presented to demonstrate the high accuracy and effectiveness of the elements developed. The standard of comparison is taken to be the solutions obtained by using 2-D plate/shell models for the beams.

Noor, Ahmed K.↗

Mixed finite element models for free vibrations of thin-walled beams

Simple, mixed finite element models are developed for the free vibration analysis of curved thin-walled beams with arbitrary open cross section. The analytical formulation is based on a Vlasov's type thin-walled beam theory with the effects of flexural-torsional coupling, transverse shear deformation and rotary inertia included. The fundamental unknowns consist of seven internal forces and seven generalized displacements of the beam. The element characteristic arrays are obtained by using a perturbed Lagrangian-mixed variational principle. Only C(sup o) continuity is required for the generalized displacements. The internal forces and the Lagrange multiplier are allowed to be discontinuous at interelement boundaries. Numerical results are presented to demonstrate the high accuracy and effectiveness of the elements developed. The standard of comparison is taken to be the solutions obtained by using 2-D plate/shell models for the beams.

Noor, Ahmed K.↗

Research aircraft observations of the mesoscale and microscale structure of a cold front over the eastern Pacific Ocean

The structure of an oceanic cold front is described on the basis of research aircraft observations taken during the Ocean Storms field experiment. Synoptic and mesoscale analyses compare the structure of an upper-level jet-front system observed slightly downstream from the wind speed maximum to its structure in the upstream entrance region. Stratospheric potential vorticity and ozone were found within the frontal zone down to about 800 mb. Microscale analyses of the front near the sea surface were carried out for a portion of the front having the signature of a 'rope' cloud in satellite imagery. A narrow (less than 1 km) zone of upward motion (about 4 m/s) and of horizontal shear (about 0.01/s) characterized the front near the surface. Significant alongfront variability was found, including lateral displacements in the frontal zone where there were weaker updrafts.

Bond, Nicholas A.↗

An assumed-stress hybrid 4-node shell element with drilling degrees of freedom

An assumed-stress hybrid/mixed 4-node quadrilateral shell element is introduced that alleviates most of the deficiencies associated with such elements. The formulation of the element is based on the assumed-stress hybrid/mixed method using the Hellinger-Reissner variational principle. The membrane part of the element has 12 degrees of freedom including rotational or 'drilling' degrees of freedom at the nodes. The bending part of the element also has 12 degrees of freedom. The bending part of the element uses the Reissner-Mindlin plate theory which takes into account the transverse shear contributions. The element formulation is derived from an 8-node isoparametric element by expressing the midside displacement degrees of freedom in terms of displacement and rotational degrees of freedom at corner nodes. The element passes the patch test, is nearly insensitive to mesh distortion, does not 'lock', possesses the desirable invariance properties, has no hidden spurious modes, and for the majority of test cases used in this paper produces more accurate results than the other elements employed herein for comparison.

Aminpour, M. A.↗

Numerical simulation of the nonlinear response of composite plates under combined thermal and acoustic loading

A time-domain study of the random response of a laminated plate subjected to combined acoustic and thermal loads is carried out. The features of this problem also include given uniform static inplane forces. The formulation takes into consideration a possible initial imperfection in the flatness of the plate. High decibel sound pressure levels along with high thermal gradients across thickness drive the plate response into nonlinear regimes. This calls for the analysis to use von Karman large deflection strain-displacement relationships. A finite element model that combines the von Karman strains with the first-order shear deformation plate theory is developed. The development of the analytical model can accommodate an anisotropic composite laminate built up of uniformly thick layers of orthotropic, linearly elastic laminae. The global system of finite element equations is then reduced to a modal system of equations. Numerical simulation using a single-step algorithm in the time-domain is then carried out to solve for the modal coordinates. Nonlinear algebraic equations within each time-step are solved by the Newton-Raphson method. The random gaussian filtered white noise load is generated using Monte Carlo simulation. The acoustic pressure distribution over the plate is capable of accounting for a grazing incidence wavefront. Numerical results are presented to study a variety of cases.

Mei, Chuh↗

Design, Optimization, and Evaluation of Integrally-Stiffened Al-2139 Panel with Curved Stiffeners

A curvilinear stiffened panel was designed, manufactured, and tested in the Combined Load Test Fixture at NASA Langley Research Center. The panel is representative of a large wing engine pylon rib and was optimized for minimum mass subjected to three combined load cases. The optimization included constraints on web buckling, material yielding, crippling or local stiffener failure, and damage tolerance using a new analysis tool named EBF3PanelOpt. Testing was performed for the critical combined compression-shear loading configuration. The panel was loaded beyond initial buckling, and strains and out-of-plane displacements were extracted from a total of 20 strain gages and 6 linear variable displacement transducers. The VIC-3D system was utilized to obtain full field displacements/strains in the stiffened side of the panel. The experimental data were compared with the strains and out-of-plane deflections from a high fidelity nonlinear finite element analysis. The experimental data were also compared with linear elastic finite element results of the panel/test-fixture assembly. Overall, the panel buckled very near to the predicted load in the web regions.

Havens, David↗

Interfacial shear behavior of sapphire-reinforced NiAl composites

The interfacial shear behavior in near-equiatomic NiAl reinforced by sapphire filaments has been examined at room temperature using a fiber pushout test technique. The load-displacement data indicate a large variability in the initial interface failure stress, although reverse push behavior indicates a comparatively constant interfacial sliding friction stress. The observed behavior suggests that the presence of asperities on the fiber surfaces and nonuniformities in fiber diameter require constrained plastic flow within the NiAl matrix in order for interfacial shear to occur. The location, shape, severity, and distribution of fiber asperities as well as the uniformity of fiber diameter are critical to the interfacial shear process.

Moose, C. A.↗

Fatigue Crack Growth and Crack Bridging in SCS-6/Ti-24-11

Interfacial damage induced by relative fiber/matrix sliding was found to occur in the bridged zone of unidirectional SCS-6/Ti-24Al-11Nb intermetallic matrix composite specimens subjected to fatigue crack growth conditions. The degree of interfacial damage was not uniform along the bridged crack wake. Higher damage zones were observed near the machined notch in comparison to the crack tip. The interfacial friction shear strength tau(sub f) measured in the crack wake using pushout testing revealed lower values than the as-received interface. Interfacial wear also reduced the strength of the bridging fibers. The reduction in fiber strength is thought to be a function of the magnitude of relative fiber/matrix displacements ind the degree of interfacial damage. Furthermore, two different fiber bridging models were used to predict the influence of bridging on the fatigue crack driving force. The shear lag model required a variable tau(sub f) in the crack wake (reflecting the degradation of the interface) before its predictions agreed with trends exhibited by the experimental data. The fiber pressure model did an excellent job in predicting both the FCG data and the DeltaCOD in the bridged zone even though it does not require a knowledge of tau(sub f).

Ghosn, Louis J.↗

Two bonded half planes with a crack going through the interface

The plane problem of two bonded elastic half planes containing a finite crack perpendicular to and going through the interface is considered. The problem is formulated as a system of singular integral equations with generalized Cauchy kernels. Even though the system has three irregular points, it is shown that the unknown functions are algebraically related at the irregular point on the interface and the integral equations can be solved by a method developed previously. The system of integral equations is shown to yield the same characteristic equation as that for two bonded quarter planes in the general case of the through crack, and the characteristic equation for a crack tip terminating at the interface in the special case. The numerical results given in the paper include the stress intensity factors at the crack tips, the normal and shear components of the stress intensity factors at the singular point on the interface, and the crack surface displacements.

Erdogan, F.↗

Evidence of a major fault zone along the California-Nevada state line 35 deg 30 min to 36 deg 30 min north latitude

The author has identified the following significant results. Geologic reconnaissance guided by analysis of ERTS-1 and Apollo-9 satellite imagery and intermediate scale photography from X-15 and U-2 aircraft has confirmed the presence of a major fault zone along the California-Nevada state line, between 35 deg 30 min and 36 deg 30 min north latitude. The name Pahrump Fault Zone has been suggested for this feature after the valley in which it is best exposed. Field reconnaissance has indicated the existence of previously unreported faults cutting bedrock along range fronts, and displacing Tertiary and Quaternary basin sediments. Gravity data support the interpretation of regional structural discontinuity along this zone. Individual fault traces within the Pahrump Fault Zone form generally left-stepping en echelon patterns. These fault patterns, the apparent offset of a Laramide age thrust fault, and possible drag folding along a major fault break suggest a component of right lateral displacement. The trend and postulated movement of the Pahrump Fault Zone are similar to the adjacent Las Vegas Shear Zone and Death Valley-Furnace Creek Faults, which are parts of a regional strike slip system in the southern Basin-Range Province.

Liggett, M. A.↗

A plate with a crack, stiffened by a partially debonded stringer

The effect of a partially debonded infinite stringer on the stress intensity factor at the crack tip is investigated. It is assumed that the stringer is bonded to an isotropic plate through an adhesive, along a line perpendicular to the crack. Due to the high stress concentration around the crack tip and on the adhesive, partial debonding is assumed to develop between the stringer and the plate. The integral equation, obtained from the continuity of displacements along the line where the stringer and the plate are bonded, is solved numerically to give the shear stress distribution between the plate and the stringer and the stress intensity factors at both crack tips.

Arin, K.↗

Recent developments in equivalent plate modeling for wing shape optimization

A new technique for structural modeling of airplane wings is presented taking transverse shear effects into account. The kinematic assumptions of first order shear deformation plate theory in combination with numerical analysis based on simple polynomials which define geometry, construction and displacement approximations lead to analytical expressions for elements of the stiffness and mass matrices and load vector. Contributions from the cover skins, spar and rib caps and spar and rib webs are included as well as concentrated springs and concentrated masses. Limitations of current equivalent plate wing modeling techniques based on classical plate theory are discussed, and the improved accuracy of the new equivalent plate technique is demonstrated through comparison to finite element analysis and test results. Analytical derivatives of stiffness, mass and load terms with respect to wing shape lead to analytic sensitivities of displacements, stresses and natural modes with respect to planform shape and depth distribution. This makes the new capability an effective structural tool for wing shape optimization.

Livne, Eli↗

Forced axisymmetric response of fluid filled spherical shells.

A general solution for the forced, linear axisymmetric response of a fluid filled spherical shell is derived with the effects of shell transverse shear and rotational inertia included. The fluid is assumed to be inviscid and compressible. Expressions for the shell radial displacement are computed for the problem of a uniform radial Heaviside load on a cap of the shell surface. In addition, the time history of the shell inner fiber stresses at the impact pole is obtained. Some comparisons with the in vacuo shell theory are also given.

Lee, Y. C.↗