Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear elasticity”

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

The structural response of unsymmetrically laminated composite cylinders

The responses of an unsymmetrically laminated fiber-reinforced composite cylinder to an axial compressive load, a torsional load, and the temperature change associated with cooling from the processing temperature to the service temperature are investigated. These problems are considered axisymmetric and the response is studied in the context of linear elastic material behavior and geometrically linear kinematics. Four different laminates are studied: a general unsymmetric laminate; two unsymmetric but more conventional laminates; and a conventional quasi-isotropic symmetric laminate. The responses based on closed-form solutions for different boundary conditions are computed and studied in detail. Particular emphasis is directed at understanding the influence of elastic couplings in the laminates. The influence of coupling decreased from a large effect in the general unsymmetric laminate, to practically no effect in the quasi-isotropic laminate. For example, the torsional loading of the general unsymmetric laminate resulted in a radial displacement. The temperature change also caused a significant radial displacement to occur near the ends of the cylinder. On the other hand, the more conventional unsymmetric laminate and the quasi-isotropic cylinder did not deform radially when subjected to a torsional load. From the results obtained, it is clear the degree of elastic coupling can be controlled and indeed designed into a cylinder, the degree and character of the coupling being dictated by the application.

Butler, T. A.↗

DYCAST: A finite element program for the crash analysis of structures

DYCAST is a nonlinear structural dynamic finite element computer code developed for crash simulation. The element library contains stringers, beams, membrane skin triangles, plate bending triangles and spring elements. Changing stiffnesses in the structure are accounted for by plasticity and very large deflections. Material nonlinearities are accommodated by one of three options: elastic-perfectly plastic, elastic-linear hardening plastic, or elastic-nonlinear hardening plastic of the Ramberg-Osgood type. Geometric nonlinearities are handled in an updated Lagrangian formulation by reforming the structure into its deformed shape after small time increments while accumulating deformations, strains, and forces. The nonlinearities due to combined loadings are maintained, and stiffness variation due to structural failures are computed. Numerical time integrators available are fixed-step central difference, modified Adams, Newmark-beta, and Wilson-theta. The last three have a variable time step capability, which is controlled internally by a solution convergence error measure. Other features include: multiple time-load history tables to subject the structure to time dependent loading; gravity loading; initial pitch, roll, yaw, and translation of the structural model with respect to the global system; a bandwidth optimizer as a pre-processor; and deformed plots and graphics as post-processors.

Pifko, A. B.↗

Recent Enhancements To The FUN3D Flow Solver For Moving-Mesh Applications

An unsteady Reynolds-averaged Navier-Stokes solver for unstructured grids has been extended to handle general mesh movement involving rigid, deforming, and overset meshes. Mesh deformation is achieved through analogy to elastic media by solving the linear elasticity equations. A general method for specifying the motion of moving bodies within the mesh has been implemented that allows for inherited motion through parent-child relationships, enabling simulations involving multiple moving bodies. Several example calculations are shown to illustrate the range of potential applications. For problems in which an isolated body is rotating with a fixed rate, a noninertial reference-frame formulation is available. An example calculation for a tilt-wing rotor is used to demonstrate that the time-dependent moving grid and noninertial formulations produce the same results in the limit of zero time-step size.

Biedron, Robert T,↗

Ultrasonic nondestructive materials characterization

A brief review of ultrasonic wave propagation in solid materials is presented with consideration of the altered behavior in anisotropic and nonlinear elastic materials in comparison with isotropic and linear elastic materials. Some experimental results are described in which ultrasonic velocity and attenuation measurements give insight into materials microstructure and associated mechanical properties. Recent developments with laser beam non-contact generation and detection of ultrasound are presented. The results of several years of experimental measurements using high-power ultrasound are discussed, which provide substantial evidence of the inability of presently accepted theories to fully explain the interaction of ultrasound with solid materials. Finally, a special synchrotron X-ray topographic system is described which affords the possibility of observing direct interaction of ultrasonic waves with the microstructural features of real crystalline solid materials for the first time.

Green, R. E., Jr.↗

Composite Overwrapped Pressure Vessel (COPV) Damage Tolerance Life Analysis Methodology and Test Best Practices

The NASA Engineering and Safety Center (NESC) Deputy Director requested an independent assessment to develop data to understand the limitations of linear elastic fracture mechanics (LEFM) computational methods used to predict fatigue crack growth rate (da/dN) behavior of small detectable cracks in thin metal liners for composite overwrapped pressure vessels (COPVs). The NESC assessment team was also requested to demonstrate a test-based methodology for validating damage tolerance requirements for COPVs with elastically responding metal liners where LEFM methods are not appropriate. This report contains the outcome of the NESC assessment.

Composite Overwrapped Pressure Vessels; Linear Ela↗

Composite Overwrapped Pressure Vessel (COPV) Damage Tolerance Life Analysis Methodology and Test Best Practices: Appendices

The NASA Engineering and Safety Center (NESC) Deputy Director requested an independent assessment to develop data to understand the limitations of linear elastic fracture mechanics (LEFM) computational methods used to predict fatigue crack growth rate (da/dN) behavior of small detectable cracks in thin metal liners for composite overwrapped pressure vessels (COPVs). The NESC assessment team was also requested to demonstrate a test-based methodology for validating damage tolerance requirements for COPVs with elastically responding metal liners where LEFM methods are not appropriate. This document contains the appendices to the main report.

Composite Overwrapped Pressure Vessels; Linear Ela↗

NASA Engineering and Safety Center Technical Bulletin No. 21-04: Evaluating Appropriateness of LEFM Tools for COPV and Metal Pressure Vessel Damage Tolerance Life Verification

Human spaceflight composite overwrapped pressure vessels (COPVs) and metal pressure vessels can use linear elastic fracture mechanics (LEFM) analysis to demonstrate damage tolerance life per ANSI/AIAA-S-081 for COPVs and ANSI/AIAA-S-080 for metal pressure vessels. LEFM analysis assumptions require that the crack tip plastic zone is small relative to the crack size and is completely surrounded by elastically responding material. Test and analysis have shown that LEFM tools (i.e., NASGRO*) can provide unconservative crack growth predictions for cracks in COPV liners that violate LEFM assumptions. COPV and metal pressure vessel designers should evaluate and address the violation of LEFM plasticity assumptions before using LEFM analysis tools for damage tolerance life verification.

COPV↗

Calculating Mode Mixity in Interlaminar Fracture of Composite Laminates

To accurately predict delamination under mixed-mode loading, cohesive elements must partition the energies dissipated in Mode I and Mode II correctly. The ability of cohesive models to predict mode mixity was evaluated by performing analyses of the single leg bending (SLB) test, which exhibits a variable mode mixity with crack length. Accordingly, two numerical integration methods were developed to compute the energy partitioning: (i) a local calculation of the accumulated work of tractions at each point along the crack path, and (ii) a global J-integral “snapshot” calculation along a contour around the process zone. In the case of a composite laminate with properties characterized by short fracture process zone (FPZ), the results of these two methods correlate well with each other and with a linear elastic fracture mechanics solution. However, FPZscan be long due to fiber bridging or due to the intrinsic toughness of the material. Then, the results of mode-mixity calculations obtained with various methods vary widely from each other. It is argued that the global J-integral calculation is the most accurate interpretation of the energy release rate and mode mixity for either short or long process zones. This solution will provide a more accurate characterization of resistance curves and mode mixity in experimental results, and it will also guide the development of more accurate mixed-mode cohesive laws for tougher materials characterized by strong R-curve responses, such as thermo-plastic composites.

cohesive zone analysis↗

Numerical, micro-mechanical prediction of crack growth resistance in a fibre-reinforced/brittle matrix composite

Micromechanics fracture models are incorporated into three distinct fracture process zones which contribute to the crack growth resistance of fibrous composites. The frontal process zone includes microcracking, fiber debonding, and some fiber failure. The elastic process zone is related only to the linear elastic creation of new matrix and fiber fracture surfaces. The wake process zone includes fiber bridging, fiber pullout, and fiber breakage. The R-curve predictions of the model compare well with empirical results for a unidirectional, continuous fiber C/C composite. Separating the contributions of each process zone reveals the wake region to contain the dominant crack growth resistance mechanisms. Fractography showed the effects of the micromechanisms on the macroscopic fracture behavior.

Jenkins, Michael G.↗

Surrogate Modeling of High-Fidelity Fracture Simulations for Real-Time Residual Strength Predictions

A surrogate model methodology is described for predicting, during flight, the residual strength of aircraft structures that sustain discrete-source damage. Starting with design of experiment, an artificial neural network is developed that takes as input discrete-source damage parameters and outputs a prediction of the structural residual strength. Target residual strength values used to train the artificial neural network are derived from 3D finite element-based fracture simulations. Two ductile fracture simulations are presented to show that crack growth and residual strength are determined more accurately in discrete-source damage cases by using an elastic-plastic fracture framework rather than a linear-elastic fracture mechanics-based method. Improving accuracy of the residual strength training data does, in turn, improve accuracy of the surrogate model. When combined, the surrogate model methodology and high fidelity fracture simulation framework provide useful tools for adaptive flight technology.

Spear, Ashley D.↗

Surrogate Modeling of High-Fidelity Fracture Simulations for Real-Time Residual Strength Predictions

A surrogate model methodology is described for predicting in real time the residual strength of flight structures with discrete-source damage. Starting with design of experiment, an artificial neural network is developed that takes as input discrete-source damage parameters and outputs a prediction of the structural residual strength. Target residual strength values used to train the artificial neural network are derived from 3D finite element-based fracture simulations. A residual strength test of a metallic, integrally-stiffened panel is simulated to show that crack growth and residual strength are determined more accurately in discrete-source damage cases by using an elastic-plastic fracture framework rather than a linear-elastic fracture mechanics-based method. Improving accuracy of the residual strength training data would, in turn, improve accuracy of the surrogate model. When combined, the surrogate model methodology and high-fidelity fracture simulation framework provide useful tools for adaptive flight technology.

Spear, Ashley D.↗

Analysis and control of hourglass instabilities in underintegrated linear and nonlinear elasticity

Methods are described to identify and correct a bad finite element approximation of the governing operator obtained when under-integration is used in numerical code for several model problems: the Poisson problem, the linear elasticity problem, and for problems in the nonlinear theory of elasticity. For each of these problems, the reason for the occurrence of instabilities is given, a way to control or eliminate them is presented, and theorems of existence, uniqueness, and convergence for the given methods are established. Finally, numerical results are included which illustrate the theory.

Jacquotte, Olivier P.↗

Generalized characteristics method for elastic wave propagation problems

Characteristic equations are derived in generalized curvilinear coordinates. Linear elastic, isotropic, and homogeneous constitutive equations have been used in the derivation. The generalized characteristic equations readily lend themselves to any requirements of space dimension and geometry. A simple boundary value problem is solved to indicate the applicability of these equations.

Ziv, M.↗

Wave-front singularities for two-dimensional anisotropic elastic waves.

Wavefront singularities for the displacement functions, associated with the radiation of linear elastic waves from a point source embedded in a finitely strained two-dimensional elastic solid, are examined in detail. It is found that generally the singularities are of order d to the -1/2 power, where d measures distance away from the front. However, in certain exceptional cases singularities of order d to the -n power, where n = 1/4, 2/3, 3/4, may be encountered.

Payton, R. G.↗

Rayleigh wave effects in an elastic half-space.

Consideration of Rayleigh wave effects in a homogeneous isotropic linearly elastic half-space subject to an impulsive uniform disk pressure loading. An approximate formula is obtained for the Rayleigh wave effects. It is shown that the Rayleigh waves near the center of loading arise from the portion of the dilatational and shear waves moving toward the axis, after they originate at the edge of the load disk. A study is made of the vertical displacement due to Rayleigh waves at points on the axis near the surface of the elastic half-space.

Aggarwal, H. R.↗

Scaling, elasticity, and CLPT

The first few viewgraphs describe the general solution properties of linear elasticity theory which are given by the following two statements: (1) for stress B.C. on S(sub sigma) and zero displacement B.C. on S(sub u) the altered displacements u(sub i)(*) and the actual stresses tau(sub ij) are elastically dependent on Poisson's ratio nu alone: thus the actual displacements are given by u(sub i) = mu(exp -1)u(sub i)(*); and (2) for zero stress B.C. on S(sub sigma) and displacement B.C. on S(sub u) the actual displacements u(sub i) and the altered stresses tau(sub ij)(*) are elastically dependent on Poisson's ratio nu alone: thus the actual stresses are given by tau(sub ij) = E tau(sub ij)(*). The remaining viewgraphs describe the minimum parameter formulation of the general classical laminate theory plate problem as follows: The general CLT plate problem is expressed as a 3 x 3 system of differential equations in the displacements u, v, and w. The eighteen (six each) A(sub ij), B(sub ij), and D(sub ij) system coefficients are ply-weighted sums of the transformed reduced stiffnesses (bar-Q(sub ij))(sub k); the (bar-Q(sub ij))(sub k) in turn depend on six reduced stiffnesses (Q(sub ij))(sub k) and the material and geometry properties of the k(sup th) layer. This paper develops a method for redefining the system coefficients, the displacement components (u,v,w), and the position components (x,y) such that a minimum parameter formulation is possible. The pivotal steps in this method are (1) the reduction of (bar-Q(sub ij))(sub k) dependencies to just two constants Q(*) = (Q(12) + 2Q(66))/(Q(11)Q(22))(exp 1/2) and F(*) - (Q(22)/Q(11))(exp 1/2) in terms of ply-independent reference values Q(sub ij); (2) the reduction of the remaining portions of the A, B, and D coefficients to nondimensional ply-weighted sums (with 0 to 1 ranges) that are independent of Q(*) and F(*); and (3) the introduction of simple coordinate stretchings for u, v, w and x,y such that the process is neatly completed.

Brunelle, Eugene J.↗

Crack growth in bonded elastic half planes

Two solutions were developed for the two dimensional problem of bonded linearly elastic half-planes. For each solution, numerical results are presented for the stress intensity factors, strain energy release rate, stresses, and displacements. The behavior predicted by the studies was investigated experimentally using polymers for the material pairs. Close agreement was found for the critical stress intensity factor at fracture for the perpendicular crack near the interface. Fracture along the interface proved to be inconclusive due to difficulties in obtaining a brittle bond. Some interesting and predictable behavior regarding the potential for the crack to cross the interface was observed and is discussed.

Goree, J. G.↗