Theory of elastic shells in the reference state.
Nonlinear field equations of elastic shell in terms of reference state subjected to strain displacement
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.
Nonlinear field equations of elastic shell in terms of reference state subjected to strain displacement
The results of residual strength pressure tests and nonlinear analyses of stringer- and frame-stiffened aluminum fuselage panels with longitudinal cracks are presented. Two types of damage are considered: a longitudinal crack located midway between stringers, and a longitudinal crack adjacent to a stringer and along a row of fasteners in a lap joint that has multiple-site damage (MSD). In both cases, the longitudinal crack is centered on a severed frame. The panels are subjected to internal pressure plus axial tension loads. The axial tension loads are equivalent to a bulkhead pressure load. Nonlinear elastic-plastic residual strength analyses of the fuselage panels are conducted using a finite element program and the crack-tip-opening-angle (CTOA) fracture criterion. Predicted crack growth and residual strength results from nonlinear analyses of the stiffened fuselage panels are compared with experimental measurements and observations. Both the test and analysis results indicate that the presence of MSD affects crack growth stability and reduces the residual strength of stiffened fuselage shells with long cracks.
The results of residual strength pressure tests and nonlinear analyses of stringer- and frame-stiffened aluminum fuselage panels with longitudinal cracks are presented. Two types of damage are considered: a longitudinal crack located midway between stringers, and a longitudinal crack adjacent to a stringer and along a row of fasteners in a lap joint that has multiple-site damage (MSD). In both cases, the longitudinal crack is centered on a severed frame. The panels are subjected to internal pressure plus axial tension loads. The axial tension loads are equivalent to a bulkhead pressure load. Nonlinear elastic-plastic residual strength analyses of the fuselage panels are conducted using a finite element program and the crack-tip-opening-angle (CTOA) fracture criterion. Predicted crack growth and residual strength results from nonlinear analyses of the stiffened fuselage panels are compared with experimental measurements and observations. Both the test and analysis results indicate that the presence of MSD affects crack growth stability and reduces the residual strength of stiffened fuselage shells with long cracks.
Various papers on recent advances in engineering science are presented. Some individual topics addressed include: advances in adaptive methods in computational fluid mechanics, mixtures of two medicomorphic materials, computer tests of rubber elasticity, shear bands in isotropic micropolar elastic materials, nonlinear surface wave and resonator effects in magnetostrictive crystals, simulation of electrically enhanced fibrous filtration, plasticity theory of granular materials, dynamics of viscoelastic media with internal oscillators, postcritical behavior of a cantilever bar, boundary value problems in nonlocal elasticity, stability of flexible structures with random parameters, electromagnetic tornadoes in earth's ionosphere and magnetosphere, helicity fluctuations and the energy cascade in turbulence, mechanics of interfacial zones in bonded materials, propagation of a normal shock in a varying area duct, analytical mechanics of fracture and fatigue.
Fatigue cracks are a major problem in designing structures subjected to cyclic loading. Cracks frequently occur in structures such as aircraft and spacecraft. The inspection intervals of many aircraft structures are based on crack-propagation lives. Therefore, improved prediction of propagation lives under flight-load conditions (variable-amplitude loading) are needed to provide more realistic design criteria for these structures. The main thrust was to develop a three-dimensional, nonlinear, elastic-plastic, finite element program capable of extending a crack and changing boundary conditions for the model under consideration. The finite-element model is composed of 8-noded (linear-strain) isoparametric elements. In the analysis, the material is assumed to be elastic-perfectly plastic. The cycle stress-strain curve for the material is shown Zienkiewicz's initial-stress method, von Mises's yield criterion, and Drucker's normality condition under small-strain assumptions are used to account for plasticity. The three-dimensional analysis is capable of extending the crack and changing boundary conditions under cyclic loading.
The first measurements of the stress induced velocity changes for propagation directions along the direction of applied stress in gr/ep composites have been presented. For propagation and stress direction perpendicular to the fiber direction, the data demonstrated a linear relation between normalized velocity shift and stress. After corrections for the delay line were made, the slope or SAC was determined and compared favorably with the expected value calculated from the previously determined nonlinear coefficients of this material. The ratio of the SAC to the elastic compliance for this direction of loading was evaluated and found to have a value similar to numerous other materials which have very different linear elastic properties. Measurements with stress and propagation along the fibers yielded unusual behavior. The curves were very nonlinear and even shifted direction at higher loads. The large scatter in the data due to bond variations made separation of material effects from bond induced artifacts impossible. Thus the SAC, R, and the remaining two unknown TOEC's could not be determined for this direction of propagation. These measurements further expand the basis of determining nonlinear elastic properties of composite materials. These properties may be useful in developing much needed NDE techniques to determine such important parameters as residual stress after cure and residual strength after impact damage. Additional study is needed to measure the nonlinear behavior in other composite materials including angle ply laminates. Also, other techniques to measure elastic nonlinearity such as harmonic generation should be applied to composites to improve the understanding of these properties and their importance.
The Space Station Freedom photovoltaic solar array blanket assembly is comprised of several layers of materials having dissimilar elastic, thermal, and mechanical properties. The operating temperature of the solar array, which ranges from -75 to +60 C, along with the material incompatibility of the blanket assembly components combine to cause an elastic-plastic stress in the weld points of the assembly. The weld points are secondary structures in nature, merely serving as electrical junctions for gathering the current. The thermal mechanical loading of the blanket assembly operating in low earth orbit continually changes throughout each 90 min orbit, which raises the possibility of fatigue induced failure. A series of structural analyses were performed in an attempt to predict the fatigue life of the solar cell in the Space Station Freedom photovoltaic array blanket. A nonlinear elastic-plastic MSC/NASTRAN analysis followed by a fatigue calculation indicated a fatigue life of 92,000 to 160,000 cycles for the solar cell weld tabs. Additional analyses predict a permanent buckling phenomenon in the copper interconnect after the first loading cycle. This should reduce or eliminate the pulling of the copper interconnect on the joint where it is welded to the silicon solar cell. It is concluded that the actual fatigue life of the solar array blanket assembly should be significantly higher than the calculated 92,000 cycles, and thus the program requirement of 87,500 cycles (orbits) will be met. Another important conclusion that can be drawn from the overall analysis is that, the strain results obtained from the MSC/NASTRAN nonlinear module are accurate to use for low-cycle fatigue analysis, since both thermal cycle testing of solar cells and analysis have shown higher fatigue life than the minimum program requirement of 87,500 cycles.
A rational methodology of modelling both nonlinear and elastic dissipative response of transversely isotropic fibrous composites is developed and illustrated with the aid of the observed response of graphite-polyimide off-axis coupons. The methodology is based on the internal variable formalism employed within the text of classical irreversible thermodynamics and entails extension of Valanis' endochronic theory to transversely isotropic media. Applicability of the theory to prediction of various response characteristics of fibrous composites is illustrated by accurately modelling such often observed phenomena as: stiffening reversible behavior along fiber direction; dissipative response in shear and transverse tension characterized by power-laws with different hardening exponents; permanent strain accumulation; nonlinear unloading and reloading; and stress-interaction effects.
This is the final report for the NASA funded project entitled "Crack Growth Prediction Methodology for Multi-Site Damage." The primary objective of the project was to create a capability to simulate curvilinear fatigue crack growth and ductile tearing in aircraft fuselages subjected to widespread fatigue damage. The second objective was to validate the capability by way of comparisons to experimental results. Both objectives have been achieved and the results are detailed herein. In the first part of the report, the crack tip opening angle (CTOA) fracture criterion, obtained and correlated from coupon tests to predict fracture behavior and residual strength of built-up aircraft fuselages, is discussed. Geometrically nonlinear, elastic-plastic, thin shell finite element analyses are used to simulate stable crack growth and to predict residual strength. Both measured and predicted results of laboratory flat panel tests and full-scale fuselage panel tests show substantial reduction of residual strength due to the occurrence of multi-site damage (MSD). Detailed comparisons of n stable crack growth history, and residual strength between the predicted and experimental results are used to assess the validity of the analysis methodology. In the second part of the report, issues related to crack trajectory prediction in thin shells; an evolving methodology uses the crack turning phenomenon to improve the structural integrity of aircraft structures are discussed, A directional criterion is developed based on the maximum tangential stress theory, but taking into account the effect of T-stress and fracture toughness orthotropy. Possible extensions of the current crack growth directional criterion to handle geometrically and materially nonlinear problems are discussed. The path independent contour integral method for T-stress evaluation is derived and its accuracy is assessed using a p- and hp-version adaptive finite element method. Curvilinear crack growth is simulated in coupon tests and in full-scale fuselage panel tests. Both T-stress and fracture toughness orthotropy are found to be essential to predict the observed crack paths. The analysis methodology and software program (FRANC3D/STAGS) developed herein allows engineers to maintain aging aircraft economically while insuring continuous airworthiness. Consequently, it will improve the technology to support the safe operation of the current aircraft fleet as well as the design of more damage-tolerant aircraft for the next generation fleet.
This paper investigates the thermomechanical response of ceramically coated metal parts in elevated thermal environments. This is made possible through the development of an improved finite element algorithm that enables the efficient and stable solution of the inherently nonlinear elastic-creep (inelastic) type thermomechanical field equations associated with high temperature. Based on the improved algorithm, the results of several numerical experiments are presented. These illustrate the significant influence of inelastic behavior in generating residual stress fields.
A new impact debris propagation code was written to link CTH simulations of space debris shield perforation to the Lagrangian finite element code DYNA3D, for space structure wall impact simulations. This software (DC3D) simulates debris cloud evolution using a nonlinear elastic-plastic deformable particle dynamics model, and renders computationally tractable the supercomputer simulation of oblique impacts on Whipple shield protected structures. Comparison of three dimensional, oblique impact simulations with experimental data shows good agreement over a range of velocities of interest in the design of orbital debris shielding. Source code developed during this research is provided on the enclosed floppy disk. An abstract based on the work described was submitted to the 1994 Hypervelocity Impact Symposium.
The results of an analytical and experimental study of the nonlinear response of thin, unstiffened, aluminum cylindrical shells with a long longitudinal crack are presented. The shells are analyzed with a nonlinear shell analysis code that accurately accounts for global and local structural response phenomena. Results are presented for internal pressure and for axial compression loads. The effect of initial crack length on the initiation of stable crack growth and unstable crack growth in typical shells subjected to internal pressure loads is predicted using geometrically nonlinear elastic-plastic finite element analyses and the crack-tip-opening angle (CTOA) fracture criterion. The results of these analyses and of the experiments indicate that the pressure required to initiate stable crack growth and unstable crack growth in a shell subjected to internal pressure loads decreases as the initial crack length increases. The effects of crack length on the prebuckling, buckling and postbuckling responses of typical shells subjected to axial compression loads are also described. For this loading condition, the crack length was not allowed to increase as the load was increased. The results of the analyses and of the experiments indicate that the initial buckling load and collapse load for a shell subjected to axial compression loads decrease as the initial crack length increases. Initial buckling causes general instability or collapse of a shell for shorter initial crack lengths. Initial buckling is a stable local response mode for longer initial crack lengths. This stable local buckling response is followed by a stable postbuckling response, which is followed by general or overall instability of the shell.
The results of an analytical and experimental study of the nonlinear response of thin, unstiffened, aluminum cylindrical shells with a long longitudinal crack are presented. The shells are analyzed with a nonlinear shell analysis code that accurately accounts for global and local structural response phenomena. Results are presented for internal pressure and for axial compression loads. The effect of initial crack length on the initiation of stable crack growth and unstable crack growth in typical shells subjected to internal pressure loads is predicted using geometrically nonlinear elastic-plastic finite element analyses and the crack-tip-opening angle (CTOA) fracture criterion. The results of these analyses and of the experiments indicate that the pressure required to initiate stable crack growth and unstable crack growth in a shell subjected to internal pressure loads decreases as the initial crack length increases. The effects of crack length on the prebuckling, buckling and postbuckling responses of typical shells subjected to axial compression loads are also described. For this loading condition, the crack length was not allowed to increase as the load was increased. The results of the analyses and of the experiments indicate that the initial buckling load and collapse load for a shell subjected to axial compression loads decrease as the initial crack length increases. Initial buckling causes general instability or collapse of a shell for shorter initial crack lengths. Initial buckling is a stable local response mode for longer initial crack lengths. This stable local buckling response is followed by a stable postbuckling response, which is followed by general or overall instability of the shell.
The objectives were to create a capability to simulate curvilinear crack growth and ductile tearing in aircraft fuselages subjected to widespread fatigue damage and to validate with tests. Analysis methodology and software program (FRANC3D/STAGS) developed herein allows engineers to maintain aging aircraft economically, while insuring continuous airworthiness, and to design more damage-tolerant aircraft for the next generation. Simulations of crack growth in fuselages were described. The crack tip opening angle (CTOA) fracture criterion, obtained from laboratory tests, was used to predict fracture behavior of fuselage panel tests. Geometrically nonlinear, elastic-plastic, thin shell finite element crack growth analyses were conducted. Comparisons of stress distributions, multiple stable crack growth history, and residual strength between measured and predicted results were made to assess the validity of the methodology. Incorporation of residual plastic deformations and tear strap failure was essential for accurate residual strength predictions. Issue related to predicting crack trajectory in fuselages were also discussed. A directional criterion, including T-stress and fracture toughness orthotropy, was developed. Curvilinear crack growth was simulated in coupon and fuselage panel tests. Both T-stress and fracture toughness orthotropy were essential to predict the observed crack paths. Flapping of fuselages were predicted. Measured and predicted results agreed reasonable well.
The results of an analytical and experimental study of the nonlinear response of thin, unstiffened, aluminum cylindrical shells with a long longitudinal crack are presented. The shells are analyzed with a nonlinear shell analysis code that accurately accounts for global and local structural response phenomena. Results are presented for internal pressure and for axial compression loads. The effect of initial crack length on the initiation of stable tearing and unstable crack growth in typical shells subjected to internal pressure loads is predicted using geo- metrically nonlinear elastic-plastic finite element analyses. The results of these analyses and of the experiments indicate that the pressure required to initiate stable tearing and unstable tearing in a shell subjected to internal pressure loads decreases as the crack length increases. The effects of crack length on the prebuckling, buckling and postbuckling responses of typical shells subjected to axial compression are also described. For this loading condition, the crack length is held constant. The results of the analyses illustrate the influence of crack length on shell buckling instabilities. The experimental and analytical results correlate well.
This paper describes our research into the vascular mechanics of the coronary artery and plaque. The three sections describe the determination of arterial mechanical properties using intravascular ultrasound (IVUS), a constitutive relation for the arterial wall, and finite element method (FEM) models of the arterial wall and atheroma. METHODS: Inflation testing of porcine left anterior descending coronary arteries was conducted. The changes in the vessel geometry were monitored using IVUS, and intracoronary pressure was recorded using a pressure transducer. The creep and quasistatic stress/strain responses were determined. A Standard Linear Solid (SLS) was modified to reproduce the non-linear elastic behavior of the arterial wall. This Standard Non-linear Solid (SNS) was implemented into an axisymetric thick-walled cylinder numerical model. Finite element analysis models were created for five age groups and four levels of stenosis using the Pathobiological Determinants of Atherosclerosis Youth (PDAY) database. RESULTS: The arteries exhibited non-linear elastic behavior. The total tissue creep strain was epsilon creep = 0.082 +/- 0.018 mm/mm. The numerical model could reproduce both the non-linearity of the porcine data and time dependent behavior of the arterial wall found in the literature with a correlation coefficient of 0.985. Increasing age had a strong positive correlation with the shoulder stress level, (r = 0.95). The 30% stenosis had the highest shoulder stress due to the combination of a fully formed lipid pool and a thin cap. CONCLUSIONS: Studying the solid mechanics of the arterial wall and the atheroma provide important insights into the mechanisms involved in plaque rupture.
The equations of motion are developed by two complementary methods, Hamilton's principle and the Newtonian method. The resulting equations are valid to second order for long, straight, slender, homogeneous, isotropic beams undergoing moderate displacements. The ordering scheme is based on the restriction that squares of the bending slopes, the torsion deformation, and the chord/radius and thickness/radius ratios are negligible with respect to unity. All remaining nonlinear terms are retained. The equations are valid for beams with mass centroid axis and area centroid (tension) axis offsets from the elastic axis, nonuniform mass and stiffness section properties, variable pretwist, and a small precone angle. The strain-displacement relations are developed from an exact transformation between the deformed and undeformed coordinate systems. These nonlinear relations form an important contribution to the final equations. Several nonlinear structural and inertial terms in the final equations are identified that can substantially influence the aeroelastic stability and response of hingeless helicopter rotor blades.
Effects of inelastic materials models on computed stress-strain solutions for thermally loaded structures were studied by performing nonlinear (elastoplastic creep) and elastic structural analyses on a prismatic, double edge wedge specimen of IN 100 alloy that was subjected to thermal cycling in fluidized beds. Four incremental plasticity creep models (isotropic, kinematic, combined isotropic kinematic, and combined plus transient creep) were exercised for the problem by using the MARC nonlinear, finite element computer program. Maximum total strain ranges computed from the elastic and nonlinear analyses agreed within 5 percent. Mean cyclic stresses, inelastic strain ranges, and inelastic work were significantly affected by the choice of inelastic constitutive model. The computing time per cycle for the nonlinear analyses was more than five times that required for the elastic analysis.