Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element 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 343 records · Page 19

On the Representation of Through-The-Thickness Reinforcements in Finite Element Analysis of Stitched, Blade Stiffened Panels

Modern aircraft employ laminated composites for their tailorable in-plane properties, high specific strengths, and high specific stiffnesses. However, laminated composites exhibit relatively poor interlaminar properties without through-the-thickness reinforcements. Quantifying the necessary amount of through-the-thickness reinforcements is necessary to reduce cost and meet damage tolerance certification requirements. In this study, a discrete superposed cohesive element (DSCE) approach is applied to represent the mixed-mode delamination behavior of stitched stiffened panels subjected to seven-point bending. This approach is compared to a one-dimensional embedded spring element (ESE) method. The DSCE approach uses two superposed bilinear traction-separation laws to obtain a representative load-displacement response determined from interlaminar tensile and shear tests. Additionally, several stitch configurations (unstitched, stitched, and overstitched) are evaluated in terms of their load-displacement response and crack-arrestment capability. Results indicate that the DSCE and ESE approaches show good agreement with respect to the predicted load-displacement response, but the ESE method tends to overpredict the crack growth behavior by approximately 13%. Stitches were not observed to fail during skin-stringer separation. Using an overstitched laminate with stitches near the flange edge provides the greatest crack-arrestment capability. Furthermore, the skin retains 92% of its stiffness after skin-stringer separation occurs.

composites↗

On the Representation of Through-The-Thickness Reinforcements in Finite Element Analysis of Stitched, Blade Stiffened Panels

Modern aircraft employ laminated composites for their tailorable in-plane properties, high specific strengths, and high specific stiffnesses. However, laminated composites exhibit relatively poor interlaminar properties without through-the-thickness reinforcements. Quantifying the necessary amount of through-the-thickness reinforcements is necessary to reduce cost and meet damage tolerance certification requirements. In this study, a discrete superposed cohesive element (DSCE) approach is applied to represent the mixed-mode delamination behavior of stitched stiffened panels subjected to seven-point bending. This approach is compared to a one-dimensional embedded spring element (ESE) method. The DSCE approach uses two superposed bilinear traction-separation laws to obtain a representative load-displacement response determined from interlaminar tensile and shear tests. Additionally, several stitch configurations (unstitched, stitched, and overstitched) are evaluated in terms of their load-displacement response and crack-arrestment capability. Results indicate that the DSCE and ESE approaches show good agreement with respect to the predicted load-displacement response, but the ESE method tends to overpredict the crack growth behavior by approximately 13%. Stitches were not observed to fail during skin-stringer separation. Using an overstitched laminate with stitches near the flange edge provides the greatest crack-arrestment capability. Furthermore, the skin retains 92% of its stiffness after skin-stringer separation occurs.

composites↗

An adaptive remeshing method for finite element thermal analysis

A finite element remeshing approach that makes use of quadrilateral and triangular elements is described. The approach uses the solution on a previous mesh to create a new mesh. Meshes are completely unstructured with highly refined elements in regions of steep gradients and larger elements where gradients are smaller. Studies of convergence rates for heat conduction problems with exact solutions show that for problems with highly localized solution variations, the remeshing approach gives smaller solution errors with fewer unknowns than refinement of uniform, structured meshes.

Thornton, Earl A.↗

Finite element methodology for transient conduction/forced-convection thermal analysis

Finite element methodology for steady state thermal analysis of convectively cooled structures has been extended for transient analysis. The finite elements are based on representing the fluid passages by fluid bulk-temperature nodes and fluid-solid interface nodes. The formulation of the finite element equations for a typical flow passage is based on the weighted residual method with upwind weighting functions. Computer implementation of the convective finite element methodology using explicit and implicit time integration algorithms is described. Accuracy and efficiency of the methodology is evaluated by comparisons with analytical solutions and finite-difference lumped-parameter analyses. The comparative analyses demonstrate that finite element conduction/conduction methodology may be used to predict transient temperatures with an accuracy equal or superior to the lumped-parameter finite-difference method.

Thornton, E. A.↗

Fluid-structure finite-element vibrational analysis

A fluid finite element has been developed for a quasi-compressible fluid. Both kinetic and potential energy are expressed as functions of nodal displacements. Thus, the formulation is similar to that used for structural elements, with the only differences being that the fluid can possess gravitational potential, and the constitutive equations for fluid contain no shear coefficients. Using this approach, structural and fluid elements can be used interchangeably in existing efficient sparse-matrix structural computer programs such as SPAR. The theoretical development of the element formulations and the relationships of the local and global coordinates are shown. Solutions of fluid slosh, liquid compressibility, and coupled fluid-shell oscillation problems which were completed using a temporary digital computer program are shown. The frequency correlation of the solutions with classical theory is excellent.

Feng, G. C.↗

Adaptive remeshing method for finite-element thermal analysis

A finite-element remeshing approach that makes use of quadrilateral and triangular elements is described. The approach uses the solution on a previous mesh to create a new mesh. Meshes are completely unstructured with highly refined elements in regions of steep gradients and larger elements where gradients are smaller. Studies of convergence rates for heat conduction problems with exact solutions show that for problems with highly localized solution variations, the remeshing approach gives smaller solution errors with fewer unknowns than refinement of uniform, structured meshes.

Thornton, Earl A.↗

Finite-Element Composite-Analysis Program

Finite Element Composite Analysis Program, FECAP, special-purpose finite-element program for analyzing behavior of composite material with microcomputer. Procedure leads to set of linear simultaneous equations relating unknown nodal displacement to applied loads. Written in HP BASIC 3.0.

Bowles, David E.↗

Simulation of Complex Cracking in Plain Weave C/SiC Composite under Biaxial Loading

Finite element analysis is performed on a mesh, based on computed geometry of a plain weave C/SiC composite with assumed internal stacking, to reveal the pattern of internal damage due to biaxial normal cyclic loading. The simulation encompasses intertow matrix cracking, matrix cracking inside the tows, and separation at the tow-intertow matrix and tow-tow interfaces. All these dissipative behaviors are represented by traction-separation cohesive laws. Not aimed at quantitatively predicting the overall stress-strain relation, the simulation, however, does not take the actual process of fiber debonding into account. The fiber tows are represented by a simple rule-of-mixture model where the reinforcing phase is a hypothetical one-dimensional material. Numerical results indicate that for the plain weave C/SiC composite, 1) matrix-crack initiation sites are primarily determined by large intertow matrix voids and interlayer tow-tow contacts, 2) the pattern of internal damage strongly depends on the loading path and initial stress, 3) compressive loading inflicts virtually no damage evolution. KEY WORDS: ceramic matrix composite, plain weave, cohesive model, brittle failure, smeared crack model, progressive damage, meso-mechanical analysis, finite element.

Cheng, Ron-Bin↗

Two-dimensional finite-element temperature variance analysis

The finite element method is extended to thermal analysis by forming a variance analysis of temperature results so that the sensitivity of predicted temperatures to uncertainties in input variables is determined. The temperature fields within a finite number of elements are described in terms of the temperatures of vertices and the variational principle is used to minimize the integral equation describing thermal potential energy. A computer calculation yields the desired solution matrix of predicted temperatures and provides information about initial thermal parameters and their associated errors. Sample calculations show that all predicted temperatures are most effected by temperature values along fixed boundaries; more accurate specifications of these temperatures reduce errors in thermal calculations.

Heuser, J. S.↗

Optimum element density studies for finite-element thermal analysis of hypersonic aircraft structures

Different finite element models previously set up for thermal analysis of the space shuttle orbiter structure are discussed and their shortcomings identified. Element density criteria are established for the finite element thermal modelings of space shuttle orbiter-type large, hypersonic aircraft structures. These criteria are based on rigorous studies on solution accuracies using different finite element models having different element densities set up for one cell of the orbiter wing. Also, a method for optimization of the transient thermal analysis computer central processing unit (CPU) time is discussed. Based on the newly established element density criteria, the orbiter wing midspan segment was modeled for the examination of thermal analysis solution accuracies and the extent of computation CPU time requirements. The results showed that the distributions of the structural temperatures and the thermal stresses obtained from this wing segment model were satisfactory and the computation CPU time was at the acceptable level. The studies offered the hope that modeling the large, hypersonic aircraft structures using high-density elements for transient thermal analysis is possible if a CPU optimization technique was used.

Ko, William L.↗

An interactive graphics system to facilitate finite element structural analysis

The characteristics of an interactive graphics systems to facilitate the finite element method of structural analysis are described. The finite element model analysis consists of three phases: (1) preprocessing (model generation), (2) problem solution, and (3) postprocessing (interpretation of results). The advantages of interactive graphics to finite element structural analysis are defined.

Burk, R. C.↗

Delamination Modeling of Composites for Improved Crash Analysis

Finite element crash modeling of composite structures is limited by the inability of current commercial crash codes to accurately model delamination growth. Efforts are made to implement and assess delamination modeling techniques using a current finite element crash code, MSC/DYTRAN. Three methods are evaluated, including a straightforward method based on monitoring forces in elements or constraints representing an interface; a cohesive fracture model proposed in the literature; and the virtual crack closure technique commonly used in fracture mechanics. Results are compared with dynamic double cantilever beam test data from the literature. Examples show that it is possible to accurately model delamination propagation in this case. However, the computational demands required for accurate solution are great and reliable property data may not be available to support general crash modeling efforts. Additional examples are modeled including an impact-loaded beam, damage initiation in laminated crushing specimens, and a scaled aircraft subfloor structures in which composite sandwich structures are used as energy-absorbing elements. These examples illustrate some of the difficulties in modeling delamination as part of a finite element crash analysis.

Fleming, David C.↗

Delamination Modeling of Composites for Improved Crash Analysis

Finite element crash modeling of composite structures is limited by the inability of current commercial crash codes to accurately model delamination growth. Efforts are made to implement and assess delamination modeling techniques using a current finite element crash code, MSC/DYTRAN. Three methods are evaluated, including a straightforward method based on monitoring forces in elements or constraints representing an interface; a cohesive fracture model proposed in the literature; and the virtual crack closure technique commonly used in fracture mechanics. Results are compared with dynamic double cantilever beam test data from the literature. Examples show that it is possible to accurately model delamination propagation in this case. However, the computational demands required for accurate solution are great and reliable property data may not be available to support general crash modeling efforts. Additional examples are modeled including an impact-loaded beam, damage initiation in laminated crushing specimens, and a scaled aircraft subfloor structures in which composite sandwich structures are used as energy-absorbing elements. These examples illustrate some of the difficulties in modeling delamination as part of a finite element crash analysis.

Fleming, David C.↗

Coating Development for GRCop-84 Liners for Reusable Launch Vehicles Aided by Modeling Studies

The design of the next generation of reusable launch vehicles calls for using GRCop-84 copper alloy liners based on a composition invented at the NASA Glenn Research Center. Despite its considerable advantage over other copper alloys, it is expected that GRCop-84 will suffer from environmental degradation depending on the type of rocket fuels used and on thermomechanical fatigue. Applying protective coatings on GRCop-84 substrates can minimize or eliminate many of these problems and extend the operational life of the combustion liner. This could increase component reliability, shorten depot maintenance turnaround times, and lower operating costs. Therefore, Glenn is actively pursuing the development of advanced coatings technology for GRCop-84 liners. Technology is being developed in four major areas: (1) new metallic coating compositions, (2) application techniques, (3) test methods, and (4) life prediction design methodology using finite element analysis. The role of finite element analysis in guiding the coating effort is discussed in this report. Thermal analyses were performed at Glenn for different combinations of top- and bondcoat compositions to determine the temperature variation across the coated cross section with the thickness of the top coat. These calculations were conducted for simulated LH2/LO2 booster engine conditions assuming that the bond coat had a constant thickness of 50 m. The preceding graphs show the predicted temperatures at the outer surface of the top coat (hot wall), at the top-coat/bond-coat interface, at the bond-coat/GRCop-84 interface, and at the GRCop-84 cold wall as a function of top-coat thickness for Cu- 26(wt%)Cr top coat (top graph), Ni-17(wt%)Cr-6%Al-0.5%Y top coat and Cu-26%Cr bond coat, and NiAl top coat and Ni bond coat. In all cases, the temperature of the top coat at the hot wall increased with increasing top-coat thickness and with corresponding decreases in the temperatures at the two interfaces and the cold wall. These temperatures are not acutely sensitive to the thermal conductivity of the top coat when it exceeds 25 and 50 W/m/K for low and high heat flux engines. This observation is significant for two reasons. First, several different top-coat compositions can be evaluated as potential protective coatings without loss in the heat-transfer efficiency of the coated system. Second, materials with thermal conductivities less than the critical values of 25 or 50 W/m/K are more likely to act as thermal barrier coatings. The deposition of overlay coatings on GRCop-84 substrates results in the development of residual stresses. The presence of these residual stresses influences the probability of coating spallation, the thermal cycling life, and the fatigue properties of the coated substrate during use. Since it is important to understand how these stresses develop during the vacuum-plasma-spraying coating deposition process, the nature and magnitudes of the cool-down residual stresses were calculated and compared with experimentally determined values across the coated cross section of a disk specimen. The calculations were conducted assuming that the specimen cools down to room temperature from vacuum plasma-spraying temperatures of either 250 or 650 C. The effects of coating the substrate with and without grit blasting were also theoretically examined. The final graph compares the predicted and the experimental results for a GRCop-84 disk coated with about a 50- m-thick Ni bond coat and a 75- to 100- m NiAl top coat, where the curves for NASA-2 assume the presence of a prior residual stress generated by grit blasting under conditions similar to the experimental situation. The predicted cool-down in-plane stresses were compressive in both the NiAl top coat and the Ni bond coat. They were also compressive in the substrate to a depth of about 0.25 mm from the Ni/GRCop-84 interface when the vacuum-plasma-spraying temperature was low. However, using a higher plasma spraying temperaturs likely to leave the substrate under a small tensile stress to counter the compressive stresses in the bond and top coats because of the relaxation of residual stresses generated in the substrate during the grit blasting of its surface prior to spraying. These results suggest that the NiAl and Ni coatings are unlikely to spall after spraying as confirmed by the microstructural observations shown in the following photomicrograph of an as-sprayed specimen. Finally, it is noted that the calculated and experimental results are not in complete agreement, which indicates that both the experimental and modeling techniques need further refinement.

Raj, Sai V.↗

Error analysis of finite element results on plates with nonuniform grids

Discrete Dirac delta functions in two dimensions and the governing plate differential equations are used to produce continuous approximations from discrete finite-element data on nonuniform grids. The resulting continuous approximation when compared to the discrete data provides an error analysis. Moreover, the continuous approximation can be differentiated to compute continuous stresses needed for component failure prediction. Stresses can then be calculated accurately at boundaries and at interfaces where finite-element analysis alone may not provide reliable results. The continuous solution along with the differential equations also provides a residual error which is used in the error analysis. This is a unique feature of this approach. Two application examples of the procedure are provided. One example involves the linear response of a clamped square plate under uniform transverse load. Results from finite-element analysis and the continuous solution from the present analysis are compared to a closed-form elasticity solution. The continuous solution is used in a Zienkiewicz-Zhu error estimator to assess the performance of the present analysis. The continuous solution is observed to dramatically correct the finite-element results. In a second example, results for the nonlinear response of a stiffened panel are compared to the finite-element solution.

Sistla, Rajaram↗

Contact stress analysis of spiral bevel gears using nonlinear finite element static analysis

A procedure is presented for performing three-dimensional stress analysis of spiral bevel gears in mesh using the finite element method. The procedure involves generating a finite element model by solving equations that identify tooth surface coordinates. Coordinate transformations are used to orientate the gear and pinion for gear meshing. Contact boundary conditions are simulated with gap elements. A solution technique for correct orientation of the gap elements is given. Example models and results are presented.

Bibel, G. D.↗

Contact stress analysis of spiral bevel gears using nonlinear finite element static analysis

A procedure is presented for performing three-dimensional stress analysis of spiral bevel gears in mesh using the finite element method. The procedure involves generating a finite element model by solving equations that identify tooth surface coordinates. Coordinate transformations are used to orientate the gear and pinion for gear meshing. Contact boundary conditions are simulated with gap elements. A solution technique for correct orientation of the gap elements is given. Example models and results are presented.

Bibel, G. D.↗

A study of the response of nonlinear springs

The various phases to developing a methodology for studying the response of a spring-reinforced arch subjected to a point load are discussed. The arch is simply supported at its ends with both the spring and the point load assumed to be at midspan. The spring is present to off-set the typical snap through behavior normally associated with arches, and to provide a structure that responds with constant resistance over a finite displacement. The various phases discussed consist of the following: (1) development of the closed-form solution for the shallow arch case; (2) development of a finite difference analysis to study (shallow) arches; and (3) development of a finite element analysis for studying more general shallow and nonshallow arches. The two numerical analyses rely on a continuation scheme to move the solution past limit points, and to move onto bifurcated paths, both characteristics being common to the arch problem. An eigenvalue method is used for a continuation scheme. The finite difference analysis is based on a mixed formulation (force and displacement variables) of the governing equations. The governing equations for the mixed formulation are in first order form, making the finite difference implementation convenient. However, the mixed formulation is not well-suited for the eigenvalue continuation scheme. This provided the motivation for the displacement based finite element analysis. Both the finite difference and the finite element analyses are compared with the closed form shallow arch solution. Agreement is excellent, except for the potential problems with the finite difference analysis and the continuation scheme. Agreement between the finite element analysis and another investigator's numerical analysis for deep arches is also good.

Hyer, M. W.↗