Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear 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 325 records · Page 18

Multiscale analysis in solids with unseparated scales: fine-scale recovery, error estimation, and coarse-scale adaptivity

There are several engineering applications in which the assumptions of homogenization and scale separation may be violated, in particular, for metallic structures constructed through additive manufacturing. Instead of resorting to direct numerical simulation of the macroscale system with an embedded fine scale, an alternative approach is to use an approximate macroscale constitutive model, but then estimate the model-form error using a posteriori error estimation techniques and subsequently adapt the macroscale model to reduce the error for a given boundary value problem and quantity of interest. Here, we investigate this approach to multiscale analysis in solids with unseparated scales using the example of an additively manufactured metallic structure consisting of a polycrystalline microstructure that is neither periodic nor statistically homogeneous. As a first step to the general nonlinear case, we focus here on linear elasticity in which each grain within the polycrystal is linear elastic but anisotropic.

42 ENGINEERING↗

Dynamic characteristics of a variable-mass flexible missile

The general motion of a variable mass flexible missile with internal flow and aerodynamic forces is considered. The resulting formulation comprises six ordinary differential equations for rigid body motion and three partial differential equations for elastic motion. The simultaneous differential equations are nonlinear and possess time-dependent coefficients. The differential equations are solved by a semi-analytical method leading to a set of purely ordinary differential equations which are then solved numerically. A computer program was developed for the numerical solution and results are presented for a given set of initial conditions.

Meirovitch, L.↗

Dynamic response analysis of geometrically non-linear structures subjected to high impact.

Description of an efficient digital computer method for the determination of the propagation of elastic stresses and deformations in certain geometrically nonlinear structures subjected to high impact loading. The finite-element matrix displacement approach utilizing curved quadrilateral shell elements in conjunction with a nodewise predictor-corrector method employing Runge-Kutta extrapolation techniques has been adopted for the present solution. The related computer program written in FORTRAN V for the UNIVAC 1108 computer has proved to be effective for the solution of a range of practical problems including rectangular and cylindrical panels. Numerical results are presented for a relevant structure, the cell container, and the negative electrode of an impact-resistant battery subjected to high impact, simulating its free landing on a planetary surface.

Gupta, K. K.↗

The effect of surface contamination on thermocapillary flow in a two-dimensional slot

In an analysis of the effect of insoluble surfactants on the steady capillary flow in a shallow slot, the problem is formulated as a two-dimensional channel heated differentially from the side, with a free surface on which an insoluble surface species is adsorbed. The vertical boundaries are considered to be isothermal and impenetrable and the lower horizontal boundary to be insulated and impenetrable. Lubrication theory and matched asymptotic expansions are used to reduce the problem to a single nonlinear integral-algebraic equation. For small elasticity E, the surface velocity is nearly constant, with the result that, for small Pe, the surface concentration is nearly uniform while, for large Pe, it develops a boundary layer near the stagnation point at the cold wall, resulting in high local surface concentrations. Condensed surfactant layers will form only with E much less than one and Pe much greater than one.

Homsy, G. M.↗

On finite element implementation and computational techniques for constitutive modeling of high temperature composites

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.

Saleeb, A. F.↗

A computational procedure for multibody systems including flexible beam dynamics

A computational procedure suitable for the solution of equations of motions for flexible multibody systems has been developed. A fully nonlinear continuum approach capable of accounting for both finite rotations and large deformations has been used to model a flexible beam component. The beam kinematics are referred directly to an inertial reference frame such that the degrees of freedom embody both the rigid and flexible deformation motions. As such, the beam inertia expression is identical to that of rigid body dynamics. The nonlinear coupling between gross body motion and elastic deformation is contained in the internal force expression. Numerical solution procedures for the integration of spatial kinematic systems can be directily applied to the generalized coordinates of both the rigid and flexible components. An accurate computation of the internal force term which is invariant to rigid motions is incorporated into the general solution procedure.

Downer, J. D.↗

A variable-order laminated plate theory based on the variational-asymptotical method

The variational-asymptotical method is a mathematical technique by which the three-dimensional analysis of laminated plate deformation can be split into a linear, one-dimensional, through-the-thickness analysis and a nonlinear, two-dimensional, plate analysis. The elastic constants used in the plate analysis are obtained from the through-the-thickness analysis, along with approximate, closed-form three-dimensional distributions of displacement, strain, and stress. In this paper, a theory based on this technique is developed which is capable of approximating three-dimensional elasticity to any accuracy desired. The asymptotical method allows for the approximation of the through-the-thickness behavior in terms of the eigenfunctions of a certain Sturm-Liouville problem associated with the thickness coordinate. These eigenfunctions contain all the necessary information about the nonhomogeneities along the thickness coordinate of the plate and thus possess the appropriate discontinuities in the derivatives of displacement. The theory is presented in this paper along with numerical results for the eigenfunctions of various laminated plates.

Lee, Bok W.↗

Structural Integrity of a Wind Tunnel Balance

The National Aeronautics and Space Administration (NASA) Langley Research Center (LaRC) has been designing strain-gage balances for utilization in wind tunnels since its inception. The utilization of balances span over a wide variety of aerodynamic tests. A force balance is an inherently critically stressed component due to the requirements of measurement sensitivity. Research and analyses are done in order to investigate the structural integrity of the balances as well as developing an understanding of their performance in order to enhance their capability. Maximum loading occurs when all 6 components of the loads are applied simultaneously with their maximum value allowed (limit load). This circumstance normally does not occur in the wind tunnel. However, if it occurs, is the balance capable of handling the loads with an acceptable factor of safety? LaRC Balance 1621 was modeled and meshed in PATRAN for analysis in NASTRAN. For a complete analysis, it is necessary to consider all the load cases as well as use dense mesh near all the edges. Because of computer limitations, it is not possible to have one model with the dense mesh near all edges. In the present study, a dense mesh is limited to the surface corners where the cage and axial sections meet. Four different load combinations are used for the current analysis. Linear analysis is performed for each load case. In the case where the stress value is above linear elastic region, it is necessary to perform nonlinear analysis. It is also important to investigate the variables limiting the structural integrity of the balances. In order to investigate the possibility of modifying the existing balances to enhance the structural integrity, some modifications are done on this balance. The structural integrity of the balance after modification is investigated.

Karkehabadi, R.↗

Parametric Study on the Response of Compression-Loaded Composite Shells With Geometric and Material Imperfections

The results of a parametric study of the effects of initial imperfections on the buckling and postbuckling response of three unstiffened thinwalled compression-loaded graphite-epoxy cylindrical shells with different orthotropic and quasi-isotropic shell-wall laminates are presented. The imperfections considered include initial geometric shell-wall midsurface imperfections, shell-wall thickness variations, local shell-wall ply-gaps associated with the fabrication process, shell-end geometric imperfections, nonuniform applied end loads, and variations in the boundary conditions including the effects of elastic boundary conditions. A high-fidelity nonlinear shell analysis procedure that accurately accounts for the effects of these imperfections on the nonlinear responses and buckling loads of the shells is described. The analysis procedure includes a nonlinear static analysis that predicts stable response characteristics of the shells and a nonlinear transient analysis that predicts unstable response characteristics.

Hilburger, Mark W.↗

Design and Analysis of a Stiffened Composite Structure Repair Concept

A design and analysis of a repair concept applicable to a stiffened thin-skin composite panel based on the Pultruded Rod Stitched Efficient Unitized Structure is presented. Since the repair concept is a bolted repair using metal components, it can easily be applied in the operational environment. Initial analyses are aimed at validating the finite element modeling approach by comparing with available test data. Once confidence in the analysis approach is established several repair configurations are explored and the most efficient one presented. Repairs involving damage to the top of the stiffener alone are considered in addition to repairs involving a damaged stiffener, flange and underlying skin. High fidelity finite element modeling techniques such as mesh-independent definition of compliant fasteners, elastic-plastic metallic material properties and geometrically nonlinear analysis are utilized in the effort. The results of the analysis are presented and factors influencing the design are assessed and discussed.

Przekop, Adam↗

Interlaminar R-Curve Effects in Thermoplastic Composites: Experiments and Model Correlations

The interlaminar fracture behavior of thermoplastic composites under Mode-I and Mode-II static loading is investigated in this research. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK)†. Characterization experiments were conducted to measure the load-displacement response, peak load, and propagation of delamination from a pre-implanted Kapton® insert. The Mode-I response was characterized using a double cantilever beam (DCB) test according to the ASTM D5528-22 standard. It was observed that the propagation exhibited significant resistance to delamination growth, i.e., an R-Curve effect. The test data indicated that the steady-state fracture toughness, G SS I , was approximately twice the mean initiation fracture toughness, G Ic . The Mode-II response of an end notch flexure (ENF) specimen tested according to the ASTM D7905-19 standard exhibited an unexpected stable behavior upon attaining peak load. Prior to attaining peak load, the Mode-II response exhibited some non-linearity, which may be attributed to nonlinear bulk material response, e.g., hyper-elasticity, viscoelasticity, plasticity, etc., and potentially friction at the load/support rollers and R-Curve effects. Finite element models based on the Cohesive Zone Method (CZM) and the Virtual Crack Closure Technique (VCCT) were used to predict the DCB and ENF specimen responses. DCB analysis results obtained using G SS I predicted propagation consistent with experimental data. However, results obtained with the as characterized G Ic were inaccurate in predicting delamination growth onset and peak load. ENF analysis results obtained with the as-characterized Mode-II fracture toughness, G IIc , were inaccurate in predicting the peak load and propagation. Hence, R-Curve effects were included in both DCB and ENF simulations, which enabled improved predictions. The experimentally determined Mode-I R-Curve was used for the DCB simulations, whereas an R-Curve was derived for Mode-II assuming that the change in compliance of the specimen was associated with crack growth only.

Prabhakar M Rao↗

Interlaminar R-Curve Effects in Thermoplastic Composites: Experiments and Model Correlations

The interlaminar fracture behavior of thermoplastic composites under Mode-I and Mode-II static loading is investigated in this research. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK)†. Characterization experiments were conducted to measure the load-displacement response, peak load, and propagation of delamination from a pre-implanted Kapton® insert. The Mode-I response was characterized using a double cantilever beam (DCB) test according to the ASTM D5528-22 standard. It was observed that the propagation exhibited significant resistance to delamination growth, i.e., an R-Curve effect. The test data indicated that the steady-state fracture toughness, G SS I , was approximately twice the mean initiation fracture toughness, G Ic . The Mode-II response of an end notch flexure (ENF) specimen tested according to the ASTM D7905-19 standard exhibited an unexpected stable behavior upon attaining peak load. Prior to attaining peak load, the Mode-II response exhibited some non-linearity, which may be attributed to nonlinear bulk material response, e.g., hyper-elasticity, viscoelasticity, plasticity, etc., and potentially friction at the load/support rollers and R-Curve effects. Finite element models based on the Cohesive Zone Method (CZM) and the Virtual Crack Closure Technique (VCCT) were used to predict the DCB and ENF specimen responses. DCB analysis results obtained using G SS I predicted propagation consistent with experimental data. However, results obtained with the as characterized G Ic were inaccurate in predicting delamination growth onset and peak load. ENF analysis results obtained with the as-characterized Mode-II fracture toughness, G IIc , were inaccurate in predicting the peak load and propagation. Hence, R-Curve effects were included in both DCB and ENF simulations, which enabled improved predictions. The experimentally determined Mode-I R-Curve was used for the DCB simulations, whereas an R-Curve was derived for Mode-II assuming that the change in compliance of the specimen was associated with crack growth only.

Prabhakar Rao↗

Compression failure of angle-ply laminates

The present work deals with modes and mechanisms of failure in compression of angle-ply laminates. Experimental results were obtained from 42 angle-ply IM7/8551-7a specimens with a lay-up of ((plus or minus theta)/(plus or minus theta)) sub 6s where theta, the off-axis angle, ranged from 0 degrees to 90 degrees. The results showed four failure modes, these modes being a function of off-axis angle. Failure modes include fiber compression, inplane transverse tension, inplane shear, and inplane transverse compression. Excessive interlaminar shear strain was also considered as an important mode of failure. At low off-axis angles, experimentally observed values were considerably lower than published strengths. It was determined that laminate imperfections in the form of layer waviness could be a major factor in reducing compression strength. Previously developed linear buckling and geometrically nonlinear theories were used, with modifications and enhancements, to examine the influence of layer waviness on compression response. The wavy layer is described by a wave amplitude and a wave length. Linear elastic stress-strain response is assumed. The geometrically nonlinear theory, in conjunction with the maximum stress failure criterion, was used to predict compression failure and failure modes for the angle-ply laminates. A range of wave length and amplitudes were used. It was found that for 0 less than or equal to theta less than or equal to 15 degrees failure was most likely due to fiber compression. For 15 degrees less than theta less than or equal to 35 degrees, failure was most likely due to inplane transverse tension. For 35 degrees less than theta less than or equal to 70 degrees, failure was most likely due to inplane shear. For theta less than 70 degrees, failure was most likely due to inplane transverse compression. The fiber compression and transverse tension failure modes depended more heavily on wave length than on wave amplitude. Thus using a single parameter, such as a ratio of wave amplitude to wave length, to describe waviness in a laminate would be inaccurate. Throughout, results for AS4/3502, studied previously, are included for comparison. At low off-axis angles, the AS4/3502 material system was found to be less sensitive to layer waviness than IM7/8551-7a. Analytical predictions were also obtained for laminates with waviness in only some of the layers. For this type of waviness, laminate compression strength could also be considered a function of which layers in the laminate were wavy, and where those wavy layers were. Overall, the geometrically nonlinear model correlates well with experimental results.

Peel, Larry D.↗

Towards an Aero-Propulso-Servo-Elasticity Analysis of a Commercial Supersonic Transport

This paper covers the development of an aero-propulso-servo-elastic (APSE) model using computational fluid dynamics (CFD) and linear structural deformations. The APSE model provides the integration of the following two previously developed nonlinear dynamic simulations: a variable cycle turbofan engine and an elastic supersonic commercial transport vehicle. The primary focus of this study is to provide a means to include relevant dynamics of a turbomachinery propulsion system into the aeroelastic studies conducted during a vehicle design, which have historically neglected propulsion effects. A high fidelity CFD tool is used here for the integration platform. The elastic vehicle neglecting the propulsion system serves as a comparison of traditional approaches to the APSE results. An overview of the methodology is presented for integrating the propulsion system and elastic vehicle. Static aeroelastic analysis comparisons between the traditional and developed APSE models for a wing tip detection indicate that the propulsion system impact on the vehicle elastic response could increase the detection by approximately ten percent.

aeroservoelasticity↗

Nonlinear Aeroelastic Equations of Motion of Twisted, Nonuniform, Flexible Horizontal-Axis Wind Turbine Blades

The second-degree nonlinear equations of motion for a flexible, twisted, nonuniform, horizontal axis wind turbine blade were developed using Hamilton's principle. A mathematical ordering scheme which was consistent with the assumption of a slender beam was used to discard some higher-order elastic and inertial terms in the second-degree nonlinear equations. The blade aerodynamic loading which was employed accounted for both wind shear and tower shadow and was obtained from strip theory based on a quasi-steady approximation of two-dimensional, incompressible, unsteady, airfoil theory. The resulting equations had periodic coefficients and were suitable for determining the aeroelastic stability and response of large horizontal-axis wind turbine blades.

Kaza, K. R. V.↗

Viscoelastic flow modeling in the extrusion of a dough-like fluid

This work attempts to investigate the effect of viscoelasticity and three-dimensional geometry in screw channels. The Phan-Thien Tanner (PTT) constitutive equation with simplified model parameters was solved in conjunction with the flow equations. Polyflow, a commercially available finite element code was used to solve the resulting nonlinear partial differential equations. The PTT model predicted one log scale lower pressure buildup compared to the equivalent Newtonian results. However, the velocity profile did not show significant changes for the chosen PTT model parameters. Past Researchers neglected viscoelastic effects and also the three dimensional nature of the flow in extruder channels. The results of this paper provide a starting point for further simulations using more realistic model parameters, which may enable the food engineer to more accurately scale-up and design extrusion processes.

NASA Discipline Life Support Systems↗

Influence of geometric nonlinearities on skin-stiffener interface stresses

A method for computing skin-stiffener interface stresses in stiffened composite panels is developed. Both geometrically linear and nonlinear analyses are considered. Particular attention is given to the flange termination region where stresses are expected to exhibit unbounded characteristics. The method is based on a finite-element analysis and an elasticity solution. The results indicate that the inclusion of geometric nonlinearities is very important for an accurate determination of the interface stresses. Membrane flattening of the panel tends to reduce the tendency of the stiffener to separate.

Cohen, D.↗

Measurement of third-order elastic constants using thermal modulation of ultrasonic waves

Third-order elastic constants (TOEC) play an important role in nonlinear material characterization, but measurements of TOEC are laborious with large error margins. This Letter presents the equations of wave velocity changes caused by homogeneous temperature variation and uniaxial stress in isotropic media and the expression of TOEC in terms of thermally induced velocity change and thermal strain. TOEC of an aluminum sample were experimentally determined by measuring ultrasonic wave velocity changes in the uniaxial loading test and the thermal modulation test. Experimental results showed good agreement between the two test methods. Owing to the simple test setup and high measurement sensitivity, the thermal modulation test is a potential experimental method to determine TOEC and absolute acoustic nonlinearity parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗