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 235 records · Page 13

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.↗

Elasticity theory in general relativity

Elasticity theory in general relativity formulated from classical nonlinear three dimensional theory, discussing thermodynamics and weak field limit

Hernandez, W. C., Jr.↗

Theoretical study of stress concentrations at circular holes and inclusions in strain hardening materials.

Nonlinear boundary value problems of an infinite elastic-plastic plate with a circular hole subjected to pure tension and pure shear at infinity are solved by a method involving Fourier series and finite difference. On the basis of these solutions, the validity of Neuber's relationship between the stress and strain concentration factors for the plane stress problems is examined and a generalized Stowell formula for the stress concentration factor is proposed for problems in which the applied loading may be pure shear as well as pure tension and, furthermore, other stress states. By the same method of solution, the stress distributions around a rigid circular cylindrical inclusion embedded in an infinite rigid-plastic matrix subjected to uniform transverse pure shear and tension are obtained.

Huang, W. C.↗

An analysis of rotor blade twist variables associated with different Euler sequences and pretwist treatments

A nonlinear analysis which is necessary to adequately model elastic helicopter rotor blades experiencing moderately large deformations was examined. The analysis must be based on an appropriate description of the blade's deformation geometry including elastic bending and twist. Built-in pretwist angles complicate the deformation process ant its definition. Relationships between the twist variables associated with different rotation sequences and corresponding forms of the transformation matrix are lasted. Relationships between the twist variables associated with first, the pretwist combined with the deformation twist are included. Many of the corresponding forms of the transformation matrix for the two cases are listed. It is shown that twist variables connected with the combined twist treatment are related to those where the pretwist is applied initially. A method to determine the relationships and some results are outlined. A procedure to evaluate the transformation matrix that eliminates the Eulerlike sequence altogether is demonstrated. The resulting form of the transformation matrix is unaffected by rotation sequence or pretwist treatment.

Alkire, K.↗

BEST3D user's manual: Boundary Element Solution Technology, 3-Dimensional Version 3.0

The theoretical basis and programming strategy utilized in the construction of the computer program BEST3D (boundary element solution technology - three dimensional) and detailed input instructions are provided for the use of the program. An extensive set of test cases and sample problems is included in the manual and is also available for distribution with the program. The BEST3D program was developed under the 3-D Inelastic Analysis Methods for Hot Section Components contract (NAS3-23697). The overall objective of this program was the development of new computer programs allowing more accurate and efficient three-dimensional thermal and stress analysis of hot section components, i.e., combustor liners, turbine blades, and turbine vanes. The BEST3D program allows both linear and nonlinear analysis of static and quasi-static elastic problems and transient dynamic analysis for elastic problems. Calculation of elastic natural frequencies and mode shapes is also provided.

Source record↗