Search NASASearch

SEARCH · Search NASA

Results for “Geometrically Nonlinear Analysis”

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 145 records · Page 8

Finite element analysis of instability related delamination growth

Instability-related growth of through-width delaminations in laminated coupons was investigated analytically and experimentally. In the analytical effort, a geometrically nonlinear finite-element stress analysis was used, and the numerical results were compared with exact solutions. In addition, measured lateral deflections of postbuckled through-width delaminations were compared with predicted deflections. Lateral deflections, stress distributions, and strain energy release rates were calculated for various delamination lengths, delamination depths, and applied loads. Delamination growth rates observed for unidirectional graphite/epoxy bonded to aluminum were qualitatively correlated with calculated strain-energy release rates and it was found that the Mode I component of energy release dominates the growth process.

Whitcomb, J. D.

Calculation of skin-stiffener interface stresses in stiffened composite panels

A method for computing the 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 finite-element analysis is standard, while the elasticity solution is based on an eigenvalue expansion of the stress functions. The eigenvalue expansion is assumed to be valid in the local flange termination region and is coupled with the finite-element analysis using collocation of stresses on the local region boundaries. Accuracy and convergence of the local elasticity solution are assessed using a geometrically linear analysis. Using this analysis procedure, the influence of geometric nonlinearities and stiffener parameters on the skin-stiffener interface stresses is evaluated.

Cohen, David

Analysis of nonlinear vibrations of cylinders

In study of geometric nonlinear vibrations, infinite, long, thin-walled cylinder was analyzed under periodic, dynamic loading to demonstrate that some nonlinear phenomena cannot be obtained by straight-forward numerical solution methods. Results demonstrate that nonlinear phenomenon in large amplitude vibration traveling-wave condition can be predicted by analysis.

Chen, J. C.

Buckling Design and Analysis of a Payload Fairing One-Sixth Cylindrical Arc-Segment Panel

Design and analysis results are reported for a panel that is a 16th arc-segment of a full 33-ft diameter cylindrical barrel section of a payload fairing structure. Six such panels could be used to construct the fairing barrel, and, as such, compression buckling testing of a 16th arc-segment panel would serve as a validation test of the buckling analyses used to design the fairing panels. In this report, linear and nonlinear buckling analyses have been performed using finite element software for 16th arc-segment panels composed of aluminum honeycomb core with graphiteepoxy composite facesheets and an alternative fiber reinforced foam (FRF) composite sandwich design. The cross sections of both concepts were sized to represent realistic Space Launch Systems (SLS) Payload Fairing panels. Based on shell-based linear buckling analyses, smaller, more manageable buckling test panel dimensions were determined such that the panel would still be expected to buckle with a circumferential (as opposed to column-like) mode with significant separation between the first and second buckling modes. More detailed nonlinear buckling analyses were then conducted for honeycomb panels of various sizes using both Abaqus and ANSYS finite element codes, and for the smaller size panel, a solid-based finite element analysis was conducted. Finally, for the smaller size FRF panel, nonlinear buckling analysis was performed wherein geometric imperfections measured from an actual manufactured FRF were included. It was found that the measured imperfection did not significantly affect the panel's predicted buckling response

composite structures

The behavior of elastic anisotropic laminated composite flat structures subjected to deterministic and random loadings

Within this research project, the following topics were studied: (1) foundation of the refined theory of flat cross-ply laminated composite flat and curved panels as well as their static and dynamic response analysis; (2) foundation of a geometrically-nonlinear shear-deformable theory of composite laminated flat panels including the effect of initial geometric imperfections and its application in the postbuckling analysis; (3) the study of the dynamic response of shear deformable elastic laminated composite panels to deterministic time-dependent external excitations as the sonic boom and explosive blast type-loadings; (4) the study of the dynamic response of shear deformable elastic laminated composite panels to random excitation as e.g. the one produced by a jet noise or by any time-dependent external excitation whose characteristics are expressed in a statistical sense; and (5) the dynamic stability of fiber-reinforced composite flat panels whose materials (due to e.g. an ambient high temperature field) exhibit a time-dependent physical behavior.

Librescu, Liviu

Advances in Structural Integrity Analysis Methods for Aging Metallic Airframe Structures with Local Damage

Analysis methodologies for predicting fatigue-crack growth from rivet holes in panels subjected to cyclic loads and for predicting the residual strength of aluminum fuselage structures with cracks and subjected to combined internal pressure and mechanical loads are described. The fatigue-crack growth analysis methodology is based on small-crack theory and a plasticity induced crack-closure model, and the effect of a corrosive environment on crack-growth rate is included. The residual strength analysis methodology is based on the critical crack-tip-opening-angle fracture criterion that characterizes the fracture behavior of a material of interest, and a geometric and material nonlinear finite element shell analysis code that performs the structural analysis of the fuselage structure of interest. The methodologies have been verified experimentally for structures ranging from laboratory coupons to full-scale structural components. Analytical and experimental results based on these methodologies are described and compared for laboratory coupons and flat panels, small-scale pressurized shells, and full-scale curved stiffened panels. The residual strength analysis methodology is sufficiently general to include the effects of multiple-site damage on structural behavior.

fatigue-crack growth

Response of nonlinear panels to random loads

Lightweight aircraft structures exposed to a high intensity noise environment can fatigue prematurely if adequate consideration is not given to the problem. Design methods and design criteria for sonic fatigue prevention were developed based on analytical and experimental techniques. Most of the analytical work was based upon small deflection or linear structural theory which did not agree with the experimental results. A large deflection geometrical nonlinearity was incorporated into the analysis methods for determining the structural response to high intensity noise. The Karman-Herrmann large deflection equations with a single mode Galerkin approximation, and the method of equivalent linearization were used to predict mean square amplitude, mean square stresses, and nonlinear frequency at various acoustic loadings for rectangular panels. Both simply supported and clamped support conditions with immovable or movable inplane edges are considered. Comparisons with experimental results are presented.

Mei, C.

Modeling Geometry and Progressive Failure of Material Interfaces in Plain Weave Composites

A procedure combining a geometrically nonlinear, explicit-dynamics contact analysis, computer aided design techniques, and elasticity-based mesh adjustment is proposed to efficiently generate realistic finite element models for meso-mechanical analysis of progressive failure in textile composites. In the procedure, the geometry of fiber tows is obtained by imposing a fictitious expansion on the tows. Meshes resulting from the procedure are conformal with the computed tow-tow and tow-matrix interfaces but are incongruent at the interfaces. The mesh interfaces are treated as cohesive contact surfaces not only to resolve the incongruence but also to simulate progressive failure. The method is employed to simulate debonding at the material interfaces in a ceramic-matrix plain weave composite with matrix porosity and in a polymeric matrix plain weave composite without matrix porosity, both subject to uniaxial cyclic loading. The numerical results indicate progression of the interfacial damage during every loading and reverse loading event in a constant strain amplitude cyclic process. However, the composites show different patterns of damage advancement.

Hsu, Su-Yuen

A triangular thin shell finite element: Nonlinear analysis

Aspects of the formulation of a triangular thin shell finite element which pertain to geometrically nonlinear (small strain, finite displacement) behavior are described. The procedure for solution of the resulting nonlinear algebraic equations combines a one-step incremental (tangent stiffness) approach with one iteration in the Newton-Raphson mode. A method is presented which permits a rational estimation of step size in this procedure. Limit points are calculated by means of a superposition scheme coupled to the incremental side of the solution procedure while bifurcation points are calculated through a process of interpolation of the determinants of the tangent-stiffness matrix. Numerical results are obtained for a flat plate and two curved shell problems and are compared with alternative solutions.

Thomas, G. R.

Three-dimensional analysis of a postbuckled embedded delamination

Delamination growth caused by local buckling of a delaminated group of plies was investigated. Delamination growth was assumed to be governed by the strain energy release rates, G(1), G(2) and G(3). The strain energy release rates were calculated using a geometrically nonlinear three-dimensional finite element analysis. The program is described and several checks of the analysis are discussed. Based on a limited parametric study, the following conclusions were reached: (1) the problem is definitely mixed mode (in some cases G(1) is larger than G(2), for other cases the opposite is true); (2) in general, there is a large gradient in the strain energy release rates along the delamination front; (3) the locations of maximum G(1) and G(2) depend on the delamination shape and the applied strain; (4) the mode 3 component was negligible for all cases considered; and (5) the analysis predicted that parts of the delamination would overlap. The results presented did not impose contact constraints to prevent overlapping. Further work is needed to determine the effects of allowing the overlapping.

Whitcomb, John D.

Three-dimensional analysis of a postbuckled embedded delamination

Delamination growth caused by local buckling of a delaminated group of plies was investigated. Delamination growth was assumed to be governed by the strain energy release rates, G(1), G(2) and G(3). The strain energy release rates were calculated using a geometrically nonlinear three-dimensional finite element analysis. The program is described and several checks of the analysis are discussed. Based on a limited parametric study, the following conclusions were reached: (1) the problem is definitely mixed mode (in some cases G(1) is larger than G(2), for other cases the opposite is true); (2) in general, there is a large gradient in the strain energy release rates along the delamination front; (3) the locations of maximum G(1) and G(2) depend on the delamination shape and the applied strain; (4) the mode 3 component was negligible for all cases considered; and (5) the analysis predicted that parts of the delamination would overlap. The results presented did not impose contact constraints to prevent overlapping. Further work is needed to determine the effects of allowing the overlapping.

Whitcomb, John D.

Nonlinear equations for dynamics of pretwisted beams undergoing small strains and large rotations

Nonlinear beam kinematics are developed and applied to the dynamic analysis of a pretwisted, rotating beam element. The common practice of assuming moderate rotations caused by structural deformation in geometric nonlinear analyses of rotating beams was abandoned in the present analysis. The kinematic relations that described the orientation of the cross section during deformation are simplified by systematically ignoring the extensional strain compared to unity in those relations. Open cross section effects such as warping rigidity and dynamics are ignored, but other influences of warp are retained. The beam cross section is not allowed to deform in its own plane. Various means of implementation are discussed, including a finite element formulation. Numerical results obtained for nonlinear static problems show remarkable agreement with experiment.

Hodges, D. H.

Error detection and control for nonlinear shell analysis

A problem-adaptive solution procedure for improving the reliability of finite element solutions to geometrically nonlinear shell-type problem is presented. The strategy incorporates automatic error detection and control and includes an iterative procedure which utilizes the solution at the same load step on a more refined model. Representative nonlinear shell problem are solved.

Mccleary, Susan L.

Efficiency of unconstrained minimization techniques in nonlinear analysis

Unconstrained minimization algorithms have been critically evaluated for their effectiveness in solving structural problems involving geometric and material nonlinearities. The algorithms have been categorized as being zeroth, first, or second order depending upon the highest derivative of the function required by the algorithm. The sensitivity of these algorithms to the accuracy of derivatives clearly suggests using analytically derived gradients instead of finite difference approximations. The use of analytic gradients results in better control of the number of minimizations required for convergence to the exact solution.

Kamat, M. P.

Progressive failure of structures

A procedure is presented for determining the nonlinear behavior of structures subjected to extreme loading and the possibility of development of potential for progressive failure. The methodology takes into account the effect of both material and geometric nonlinearities. At a given stage of analysis, the individual components of the structure are checked against predetermined failure criteria. Subsequently, the failing components are removed and the modified structure is analyzed for overall failure. Examples, obtained from a computer program based on the proposed procedure, showing the applicability of the method are presented.

Khozeimeh, K.

Helicopter aeroelastic stability and response - Current topics and future trends

This paper presents several current topics in rotary wing aeroelasticity and concludes by attempting to anticipate future trends and developments. These topics are: (1) the role of geometric nonlinearities; (2) structural modeling, and aeroelastic analysis of composite rotor blades; (3) aeroelastic stability and response in forward flight; (4) modeling of coupled rotor/fuselage aeromechanical problems and their active control; and (5) the coupled rotor-fuselage vibration problem and its alleviation by higher harmonic control. Selected results illustrating the fundamental aspects of these topics are presented. Future developments are briefly discussed.

Friedmann, Peretz P.

Optimization of Car Body under Constraints of Noise, Vibration, and Harshness (NVH), and Crash

To be competitive on the today's market, cars have to be as light as possible while meeting the Noise, Vibration, and Harshness (NVH) requirements and conforming to Government-man dated crash survival regulations. The latter are difficult to meet because they involve very compute-intensive, nonlinear analysis, e.g., the code RADIOSS capable of simulation of the dynamics, and the geometrical and material nonlinearities of a thin-walled car structure in crash, would require over 12 days of elapsed time for a single design of a 390K elastic degrees of freedom model, if executed on a single processor of the state-of-the-art SGI Origin2000 computer. Of course, in optimization that crash analysis would have to be invoked many times. Needless to say, that has rendered such optimization intractable until now. The car finite element model is shown. The advent of computers that comprise large numbers of concurrently operating processors has created a new environment wherein the above optimization, and other engineering problems heretofore regarded as intractable may be solved. The procedure, shown, is a piecewise approximation based method and involves using a sensitivity based Taylor series approximation model for NVH and a polynomial response surface model for Crash. In that method the NVH constraints are evaluated using a finite element code (MSC/NASTRAN) that yields the constraint values and their derivatives with respect to design variables. The crash constraints are evaluated using the explicit code RADIOSS on the Origin 2000 operating on 256 processors simultaneously to generate data for a polynomial response surface in the design variable domain. The NVH constraints and their derivatives combined with the response surface for the crash constraints form an approximation to the system analysis (surrogate analysis) that enables a cycle of multidisciplinary optimization within move limits. In the inner loop, the NVH sensitivities are recomputed to update the NVH approximation model while keeping the Crash response surface constant. In every outer loop, the Crash response surface approximation is updated, including a gradual increase in the order of the response surface and the response surface extension in the direction of the search. In this optimization task, the NVH discipline has 30 design variables while the crash discipline has 20 design variables. A subset of these design variables (10) are common to both the NVH and crash disciplines. In order to construct a linear response surface for the Crash discipline constraints, a minimum of 21 design points would have to be analyzed using the RADIOSS code. On a single processor in Origin 2000 that amount of computing would require over 9 months! In this work, these runs were carried out concurrently on the Origin 2000 using multiple processors, ranging from 8 to 16, for each crash (RADIOSS) analysis. Another figure shows the wall time required for a single RADIOSS analysis using varying number of processors, as well as provides a comparison of 2 different common data placement procedures within the allotted memories for each analysis. The initial design is an infeasible design with NVH discipline Static Torsion constraint violations of over 10%. The final optimized design is a feasible design with a weight reduction of 15 kg compared to the initial design. This work demonstrates how advanced methodology for optimization combined with the technology of concurrent processing enables applications that until now were out of reach because of very long time-to-solution.

Kodiyalam, Srinivas

Nonlinear Analysis of Bonded Composite Single-LAP Joints

This study presents a semi-analytical solution method to analyze the geometrically nonlinear response of bonded composite single-lap joints with tapered adherend edges under uniaxial tension. The solution method provides the transverse shear and normal stresses in the adhesive and in-plane stress resultants and bending moments in the adherends. The method utilizes the principle of virtual work in conjunction with von Karman s nonlinear plate theory to model the adherends and the shear lag model to represent the kinematics of the thin adhesive layer between the adherends. Furthermore, the method accounts for the bilinear elastic material behavior of the adhesive while maintaining a linear stress-strain relationship in the adherends. In order to account for the stiffness changes due to thickness variation of the adherends along the tapered edges, their in-plane and bending stiffness matrices are varied as a function of thickness along the tapered region. The combination of these complexities results in a system of nonlinear governing equilibrium equations. This approach represents a computationally efficient alternative to finite element method. Comparisons are made with corresponding results obtained from finite-element analysis. The results confirm the validity of the solution method. The numerical results present the effects of taper angle, adherend overlap length, and the bilinear adhesive material on the stress fields in the adherends, as well as the adhesive, of a single-lap joint

Oterkus, E.