Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mixed finite element”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Advances in reduction techniques for tire contact problems

Some recent developments in reduction techniques, as applied to predicting the tire contact response and evaluating the sensitivity coefficients of the different response quantities, are reviewed. The sensitivity coefficients measure the sensitivity of the contact response to variations in the geometric and material parameters of the tire. The tire is modeled using a two-dimensional laminated anisotropic shell theory with the effects of variation in geometric and material parameters, transverse shear deformation, and geometric nonlinearities included. The contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the stress resultants, the generalized displacements, and the Lagrange multipliers associated with the contact conditions. The elemental arrays are obtained by using a modified two-field, mixed variational principle. For the application of reduction techniques, the tire finite element model is partitioned into two regions. The first region consists of the nodes that are likely to come in contact with the pavement, and the second region includes all the remaining nodes. The reduction technique is used to significantly reduce the degrees of freedom in the second region. The effectiveness of the computational procedure is demonstrated by a numerical example of the frictionless contact response of the space shuttle nose-gear tire, inflated and pressed against a rigid flat surface. Also, the research topics which have high potential for enhancing the effectiveness of reduction techniques are outlined.

Noor, Ahmed K.↗
NASA NTRS1976

Modeling of three-dimensional mixing and reacting ducted flows

A computer code, based upon a finite element solution algorithm, was developed to solve the governing equations for three-dimensional, reacting boundary region, and constant area ducted flow fields. Effective diffusion coefficients are employed to allow analyses of turbulent, transitional or laminar flows. The code was used to investigate mixing and reacting hydrogen jets injected from multiple orifices, transverse and parallel to a supersonic air stream. Computational results provide a three-dimensional description of velocity, temperature, and species-concentration fields downstream of injection. Experimental data for eight cases covering different injection conditions and geometries were modeled using mixing length theory (MLT). These results were used as a baseline for examining the relative merits of other mixing models. Calculations were made using a two-equation turbulence model (k+d) and comparisons were made between experiment and mixing length theory predictions. The k+d model shows only a slight improvement in predictive capability over MLT. Results of an examination of the effect of tensorial transport coefficients on mass and momentum field distribution are also presented. Solutions demonstrating the ability of the code to model ducted flows and parallel strut injection are presented and discussed.

Zelazny, S. W.↗
NASA NTRS1992

Improved accuracy for finite element structural analysis via a new integrated force method

A comparative study was carried out to determine the accuracy of finite element analyses based on the stiffness method, a mixed method, and the new integrated force and dual integrated force methods. The numerical results were obtained with the following software: MSC/NASTRAN and ASKA for the stiffness method; an MHOST implementation method for the mixed method; and GIFT for the integrated force methods. The results indicate that on an overall basis, the stiffness and mixed methods present some limitations. The stiffness method generally requires a large number of elements in the model to achieve acceptable accuracy. The MHOST method tends to achieve a higher degree of accuracy for course models than does the stiffness method implemented by MSC/NASTRAN and ASKA. The two integrated force methods, which bestow simultaneous emphasis on stress equilibrium and strain compatibility, yield accurate solutions with fewer elements in a model. The full potential of these new integrated force methods remains largely unexploited, and they hold the promise of spawning new finite element structural analysis tools.

Patnaik, Surya N.↗
NASA NTRS1992

Improved accuracy for finite element structural analysis via an integrated force method

A comparative study was carried out to determine the accuracy of finite element analyses based on the stiffness method, a mixed method, and the new integrated force and dual integrated force methods. The numerical results were obtained with the following software: MSC/NASTRAN and ASKA for the stiffness method; an MHOST implementation method for the mixed method; and GIFT for the integrated force methods. The results indicate that on an overall basis, the stiffness and mixed methods present some limitations. The stiffness method generally requires a large number of elements in the model to achieve acceptable accuracy. The MHOST method tends to achieve a higher degree of accuracy for course models than does the stiffness method implemented by MSC/NASTRAN and ASKA. The two integrated force methods, which bestow simultaneous emphasis on stress equilibrium and strain compatibility, yield accurate solutions with fewer elements in a model. The full potential of these new integrated force methods remains largely unexploited, and they hold the promise of spawning new finite element structural analysis tools.

Patnaik, S. N.↗
NASA NTRS1995

Reduction technique for tire contact problems

A reduction technique and a computational procedure are presented for predicting the tire contact response and evaluating the sensitivity coefficients of the different response quantities. The sensitivity coefficients measure the sensitivity of the contact response to variations in the geometric and material parameters of the tire. The tire is modeled using a two-dimensional laminated anisotropic shell theory with the effects of variation in geometric and material parameters, transverse shear deformation, and geometric nonlinearities included. The contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the stress resultants, the generalized displacements, and the Lagrange multipliers associated with the contact conditions. The elemental arrays are obtained by using a modified two-field, mixed variational principle. For the application of the reduction technique, the tire finite element model is partitioned into two regions. The first region consists of the nodes that are likely to come in contact with the pavement, and the second region includes all the remaining nodes. The reduction technique is used to significantly reduce the degrees of freedom in the second region. The effectiveness of the computational procedure is demonstrated by a numerical example of the frictionless contact response of the space shuttle nose-gear tire, inflated and pressed against a rigid flat surface.

Noor, Ahmed K.↗
NASA NTRS2008

Mixed Mode Fracture of Plasma Sprayed Thermal Barrier Coatings: Effects of Anisotropy and Heterogeneity

The combined mode I-mode II fracture behavior of anisotropic ZrO2-8wt%Y2O3 thermal barrier coatings was determined in asymmetric flexure loading at both ambient and elevated temperatures. A fracture envelope of KI versus KII was determined for the coating material at ambient and elevated temperatures. Propagation angles of fracture as a function of KI/KII were also determined. The mixed-mode fracture behavior of the microsplat coating material was modeled using Finite Element approach to account for anisotropy and micro cracked structures, and predicted in terms of fracture envelope and propagation angle using mixed-mode fracture theories.

Zhu, Dongming↗
NASA NTRS1989

Finite element solvers for incompressible fluid flows and heat transfer

Two different finite-element solvers for incompressible viscous flow, i.e., the mixed interpolation method and the SIMPLE-type iterative method, are compared and tested with some benchmark problems. The advantages of the SIMPLE-type iterative method are the decoupling of the governing equations and the use of equal-order interpolation functions for both velocity and pressure. Even though there is a significant difference between the two methods in terms of the pressure field, similar solutions are obtained for the velocity field.

Sohn, Jeong L.↗
NASA NTRS2002

Postbuckling and Growth of Delaminations in Composite Plates Subjected to Axial Compression

The postbuckling response and growth of circular delaminations in flat and curved plates are investigated as part of a study to identify the criticality of delamination locations through the laminate thickness. The experimental results from tests on delaminated plates are compared with finite element analysis results generated using shell models. The analytical prediction of delamination growth is obtained by assessing the strain energy release rate results from the finite element model and comparing them to a mixed-mode fracture toughness failure criterion. The analytical results for onset of delamination growth compare well with experimental results generated using a 3-dimensional displacement visualization system. The record of delamination progression measured in this study has resulted in a fully 3-dimensional test case with which progressive failure models can be validated.

Reeder, James R.↗
NASA NTRS2012

Application of Benchmark Examples to Assess the Single and Mixed-Mode Static Delamination Propagation Capabilities in ANSYS

The application of benchmark examples for the assessment of quasi-static delamination propagation capabilities is demonstrated for ANSYS. The examples are independent of the analysis software used and allow the assessment of the automated delamination propagation in commercial finite element codes based on the virtual crack closure technique (VCCT). The examples selected are based on two-dimensional finite element models of Double Cantilever Beam (DCB), End-Notched Flexure (ENF), Mixed-Mode Bending (MMB) and Single Leg Bending (SLB) specimens. First, the quasi-static benchmark examples were recreated for each specimen using the current implementation of VCCT in ANSYS . Second, the delamination was allowed to propagate under quasi-static loading from its initial location using the automated procedure implemented in the finite element software. Third, the load-displacement relationship from a propagation analysis and the benchmark results were compared, and good agreement could be achieved by selecting the appropriate input parameters. The benchmarking procedure proved valuable by highlighting the issues associated with choosing the input parameters of the particular implementation. Overall the results are encouraging, but further assessment for three-dimensional solid models is required.

Krueger, Ronald↗
NASA NTRS2006

On the Development of an Efficient Parallel Hybrid Solver with Application to Acoustically Treated Aero-Engine Nacelles

A finite element solution to the convected Helmholtz equation in a nonuniform flow is used to model the noise field within 3-D acoustically treated aero-engine nacelles. Options to select linear or cubic Hermite polynomial basis functions and isoparametric elements are included. However, the key feature of the method is a domain decomposition procedure that is based upon the inter-mixing of an iterative and a direct solve strategy for solving the discrete finite element equations. This procedure is optimized to take full advantage of sparsity and exploit the increased memory and parallel processing capability of modern computer architectures. Example computations are presented for the Langley Flow Impedance Test facility and a rectangular mapping of a full scale, generic aero-engine nacelle. The accuracy and parallel performance of this new solver are tested on both model problems using a supercomputer that contains hundreds of central processing units. Results show that the method gives extremely accurate attenuation predictions, achieves super-linear speedup over hundreds of CPUs, and solves upward of 25 million complex equations in a quarter of an hour.

Watson, Willie R.↗
NASA NTRS1999

Crack Growth Prediction Methodology for Multi-Site Damage: Layered Analysis and Growth During Plasticity

A finite element program has been developed to perform quasi-static, elastic-plastic crack growth simulations. The model provides a general framework for mixed-mode I/II elastic-plastic fracture analysis using small strain assumptions and plane stress, plane strain, and axisymmetric finite elements. Cracks are modeled explicitly in the mesh. As the cracks propagate, automatic remeshing algorithms delete the mesh local to the crack tip, extend the crack, and build a new mesh around the new tip. State variable mapping algorithms transfer stresses and displacements from the old mesh to the new mesh. The von Mises material model is implemented in the context of a non-linear Newton solution scheme. The fracture criterion is the critical crack tip opening displacement, and crack direction is predicted by the maximum tensile stress criterion at the crack tip. The implementation can accommodate multiple curving and interacting cracks. An additional fracture algorithm based on nodal release can be used to simulate fracture along a horizontal plane of symmetry. A core of plane strain elements can be used with the nodal release algorithm to simulate the triaxial state of stress near the crack tip. Verification and validation studies compare analysis results with experimental data and published three-dimensional analysis results. Fracture predictions using nodal release for compact tension, middle-crack tension, and multi-site damage test specimens produced accurate results for residual strength and link-up loads. Curving crack predictions using remeshing/mapping were compared with experimental data for an Arcan mixed-mode specimen. Loading angles from 0 degrees to 90 degrees were analyzed. The maximum tensile stress criterion was able to predict the crack direction and path for all loading angles in which the material failed in tension. Residual strength was also accurately predicted for these cases.

James, Mark Anthony↗
NASA NTRS1979

Recent developments in finite element analysis for transonic airfoils

The prediction of aerodynamic forces in the transonic regime generally requires a flow field calculation to solve the governing non-linear mixed elliptic-hyperbolic partial differential equations. Finite difference techniques were developed to the point that design and analysis application are routine, and continual improvements are being made by various research groups. The principal limitation in extending finite difference methods to complex three-dimensional geometries is the construction of a suitable mesh system. Finite element techniques are attractive since their application to other problems have permitted irregular mesh elements to be employed. The purpose of this paper is to review the recent developments in the application of finite element methods to transonic flow problems and to report some recent results.

Hafez, M. M.↗
NASA NTRS2008

Guidelines and Parameter Selection for the Simulation of Progressive Delamination

Turon s methodology for determining optimal analysis parameters for the simulation of progressive delamination is reviewed. Recommended procedures for determining analysis parameters for efficient delamination growth predictions using the Abaqus/Standard cohesive element and relatively coarse meshes are provided for single and mixed-mode loading. The Abaqus cohesive element, COH3D8, and a user-defined cohesive element are used to develop finite element models of the double cantilever beam specimen, the end-notched flexure specimen, and the mixed-mode bending specimen to simulate progressive delamination growth in Mode I, Mode II, and mixed-mode fracture, respectively. The predicted responses are compared with their analytical solutions. The results show that for single-mode fracture, the predicted responses obtained with the Abaqus cohesive element correlate well with the analytical solutions. For mixed-mode fracture, it was found that the response predicted using COH3D8 elements depends on the damage evolution criterion that is used. The energy-based criterion overpredicts the peak loads and load-deflection response. The results predicted using a tabulated form of the BK criterion correlate well with the analytical solution and with the results predicted with the user-written element.

Song, Kyongchan↗
NASA NTRS1990

Mixed-Mode Bending Method for Delamination Testing

A mixed mode delamination test procedure was developed combining double cantilever beam (DCB) mode I loading and end-notch fixture (ENF) mode II loading on a split unidirectional laminate. By loading with a lever, a single applied load simultaneously produces mode I and mode II bending loads on the specimen. This mixed-mode bending (MMB) test was analyzed using both finite-element procedures and beam theory to calculate the mode I and mode II components of strain-energy release rate G(sub I) and G(sub II), respectively. A wide range of G(sub I)/G(sub II) ratios can be produced by varying the load position on the lever. As the delamination extended, the G(sub I)/G(sub II) ratios varied by less than 5%. Beam theory equations agreed closely with the finite-element results and provide a basis for selection of G(sub I)/G(sub II) test ratios and a basis for computing the mode I and mode II components of measured delamination toughness. The MMB test was demonstrated using AS4/PEEK (APC2) unidirectional laminates. The MMB test introduced in this paper is rather simple and is believed to offer several advantages over most current mixed-mode test.

Reeder, James R.↗
NASA NTRS1988

A mixed-mode bending apparatus for delamination testing

A mixed-mode delamination test procedure was developed combining double cantilever beam mode I loading and end notch flexure mode II loading on a split unidirectional laminate. By loading the specimen with a lever, a single applied load simultaneously produces mode I and II bending loads on the specimen. This mixed mode bending (MMB) test was analyzed using both finite element procedures and beam theory to calculate the mode I and II components of strain energy release rate, G sub I and G sub II, respectively. The analyses showed that a wide range of G sub I/G sub II ratios could be produced by varying the applied load position on the loading lever. As the delamination extended, the G sub I/G sub II ratios varied by less than 5 percent. The simple beam theory equations were modified to account for the elastic interaction between the two arms of the specimen and to account for shear deformations. The resulting equations agreed closely with the finite element results and provide a basis for selection of G sub I/G sub II test ratios and a basis for computing the mode I and II components of measured delamination toughness. The MMB specimen analysis and test procedures were demonstrated using unidirectional laminates.

Crews, John H., Jr.↗
NASA NTRS2013

Development and Applications of Benchmark Examples for Static Delamination Propagation Predictions

The development and application of benchmark examples for the assessment of quasistatic delamination propagation capabilities was demonstrated for ANSYS (TradeMark) and Abaqus/Standard (TradeMark). The examples selected were based on finite element models of Double Cantilever Beam (DCB) and Mixed-Mode Bending (MMB) specimens. First, quasi-static benchmark results were created based on an approach developed previously. Second, the delamination was allowed to propagate under quasi-static loading from its initial location using the automated procedure implemented in ANSYS (TradeMark) and Abaqus/Standard (TradeMark). Input control parameters were varied to study the effect on the computed delamination propagation. Overall, the benchmarking procedure proved valuable by highlighting the issues associated with choosing the appropriate input parameters for the VCCT implementations in ANSYS® and Abaqus/Standard®. However, further assessment for mixed-mode delamination fatigue onset and growth is required. Additionally studies should include the assessment of the propagation capabilities in more complex specimens and on a structural level.

Krueger, Ronald↗
NASA NTRS1985

Dynamic fracture mechanics

Dynamic fracture and crack propagation concepts for ductile materials are reviewed. The equations for calculating dynamic stress integrity and the dynamic energy release rate in order to study dynamic crack propagation are provided. The stress intensity factor versus crack velocity relation is investigated. The uses of optical experimental techniques and finite element methods for fracture analyses are described. The fracture criteria for a rapidly propagating crack under mixed mode conditions are discussed; crack extension and fracture criteria under combined tension and shear loading are based on maximum circumferential stress or energy criteria such as strain energy density. The development and use of a Dugdale model and finite element models to represent crack and fracture dynamics are examined.

Kobayashi, A. S.↗