Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Virtual Element Method”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

A General Interface Method for Aeroelastic Analysis of Aircraft

The aeroelastic analysis of an aircraft requires an accurate and efficient procedure to couple aerodynamics and structures. The procedure needs an interface method to bridge the gap between the aerodynamic and structural models in order to transform loads and displacements. Such an interface method is described in this report. This interface method transforms loads computed by any aerodynamic code to a structural finite element (FE) model and converts the displacements from the FE model to the aerodynamic model. The approach is based on FE technology in which virtual work is employed to transform the aerodynamic pressures into FE nodal forces. The displacements at the FE nodes are then converted back to aerodynamic grid points on the aircraft surface through the reciprocal theorem in structural engineering. The method allows both high and crude fidelities of both models and does not require an intermediate modeling. In addition, the method performs the conversion of loads and displacements directly between individual aerodynamic grid point and its corresponding structural finite element and, hence, is very efficient for large aircraft models. This report also describes the application of this aero-structure interface method to a simple wing and an MD-90 wing. The results show that the aeroelastic effect is very important. For the simple wing, both linear and nonlinear approaches are used. In the linear approach, the deformation of the structural model is considered small, and the loads from the deformed aerodynamic model are applied to the original geometry of the structure. In the nonlinear approach, the geometry of the structure and its stiffness matrix are updated in every iteration and the increments of loads from the previous iteration are applied to the new structural geometry in order to compute the displacement increments. Additional studies to apply the aero-structure interaction procedure to more complicated geometry will be conducted in the second phase of the present contract.

Tzong, T.↗
NASA NTRS1987

Physically based sensitivity derivatives for finite-element programs

Attention is given to a virtual-work approach to the definition of sensitivity in terms of physical quantities; sensitivity fields are treated as a solution to loadings that consist of body forces, initial strains, and initial stresses. The method can accordingly be implemented in a finite element package, with only minimal knowledge of its programming detail and without access to its source code. The application of the approach to one- and two-dimensional problems is demonstrated. The sensitivity of the structural response for truss structures was computed with respect to the cross-sectional areas of the members.

Barthelemy, Bruno↗
NASA NTRS1993

Non-oscillatory and non-diffusive solution of convection problems by the iteratively reweighted least-squares finite element method

A comparative description is presented for the least-squares FEM (LSFEM) for 2D steady-state pure convection problems. In addition to exhibiting better control of the streamline derivative than the streamline upwinding Petrov-Galerkin method, numerical convergence rates are obtained which show the LSFEM to be virtually optimal. The LSFEM is used as a framework for an iteratively reweighted LSFEM yielding nonoscillatory and nondiffusive solutions for problems with contact discontinuities; this method is shown to convect contact discontinuities without error when using triangular and bilinear elements.

Jiang, Bo-Nan↗
NASA NTRS2002

Formation metrology and control for large separated optics space telescopes

In this paper we present formation flying performance analysis initial results for a representative large space telescope composed of separated optical elements [Mett 02]. A virtual-structure construct (an equivalent rigid body) is created by unique metrology and control that combines both centralized and decentralized methods. The formation may be in orbit at GEO for super-resolution Earth observation, as in the case of Figure 1, or it may be in an Earth-trailing orbit for astrophysics, Figure 2. Extended applications are envisioned for exo-solar planet interferometric imaging by a formation of very large separated optics telescopes, Figure 3. Space telescopes, with such large apertures and f/10 to f/100 optics, are not feasible if connected by massive metering structures. Instead, the new virtual-structure paradigm of information and control connectivity between the formation elements provides the necessary spatial rigidity and alignment precision for the telescope.

formation flying virtual-structure telescopes↗
NASA NTRS1999

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

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

Mesh Convergence Requirements for Composite Damage Models

The ability of the finite element method to accurately represent the response of objects with intricate geometry and loading renders the finite element method as an extremely versatile analysis technique for structural analysis. Finite element analysis is routinely used in industry to calculate deflections, stress concentrations, natural frequencies, buckling loads, and much more. The method works by discretizing complex problems into smaller, simpler approximations that are valid over small uniform domains. For common analyses, the maximum size of the elements that can be used is often be determined by experience. However, to verify the quality of a solution, analyses with several levels of mesh refinement should be performed to ensure that the solution has converged. In recent years, the finite element method has been used to calculate the resistance of structures, and in particular that of composite structures. A number of techniques such as cohesive zone modeling, the virtual crack closure technique, and continuum damage modeling have emerged that can be used to predict cracking, delaminations, fiber failure, and other composite damage modes that lead to structural collapse. However, damage models present mesh refinement requirements that are not well understood. In this presentation, we examine different mesh refinement issues related to the representation of damage in composite materials. Damage process zone sizes and their corresponding mesh requirements will be discussed. The difficulties of modeling discontinuities and the associated need for regularization techniques will be illustrated, and some unexpected element size constraints will be presented. Finally, some of the difficulties in constructing models of composite structures capable of predicting transverse matrix cracking will be discussed. It will be shown that to predict the initiation and propagation of transverse matrix cracks, their density, and their saturation may require models that are significantly more refined than those that have been contemplated in the past.

Davila, Carlos G.↗
NASA NTRS2001

Numerical Simulation of Delamination Growth in Composite Materials

The use of decohesion elements for the simulation of delamination in composite materials is reviewed. The test methods available to measure the interfacial fracture toughness used in the formulation of decohesion elements are described initially. After a brief presentation of the virtual crack closure technique, the technique most widely used to simulate delamination growth, the formulation of interfacial decohesion elements is described. Problems related with decohesion element constitutive equations, mixed-mode crack growth, element numerical integration and solution procedures are discussed. Based on these investigations, it is concluded that the use of interfacial decohesion elements is a promising technique that avoids the need for a pre-existing crack and pre-defined crack paths, and that these elements can be used to simulate both delamination onset and growth.

Camanho, P. P.↗
NASA NTRS1990

Evaluation of finite-element models and stress-intensity factors for surface cracks emanating from stress concentrations

A re-evaluation of the 3-D finite-element models and methods used to analyze surface crack at stress concentrations is presented. Previous finite-element models used by Raju and Newman for surface and corner cracks at holes were shown to have ill-shaped elements at the intersection of the hole and crack boundaries. Improved models, without these ill-shaped elements, were developed for a surface crack at a circular hole and at a semi-circular edge notch. Stress-intensity factors were calculated by both the nodal-force and virtual-crack-closure methods. Comparisons made between the previously developed stress-intensity factor equations and the results from the improved models agreed well except for configurations with large notch-radii-to-plate-thickness ratios. Stress-intensity factors for a semi-elliptical surface crack located at the center of a semi-circular edge notch in a plate subjected to remote tensile loadings were calculated using the improved models.

Tan, P. W.↗
NASA NTRS2013

The Virtual Observatory Service TheoSSA: Establishing a Database of Synthetic Stellar Flux Standards I. NLTE Spectral Analysis of the DA-Type White Dwarf G191-B2B *,**,***,****

Hydrogen-rich, DA-type white dwarfs are particularly suited as primary standard stars for flux calibration. State-of-the-art NLTE models consider opacities of species up to trans-iron elements and provide reliable synthetic stellar-atmosphere spectra to compare with observations. Aims. We will establish a database of theoretical spectra of stellar flux standards that are easily accessible via a web interface. Methods. In the framework of the Virtual Observatory, the German Astrophysical Virtual Observatory developed the registered service TheoSSA. It provides easy access to stellar spectral energy distributions (SEDs) and is intended to ingest SEDs calculated by any model-atmosphere code. In case of the DA white dwarf G191−B2B, we demonstrate that the model reproduces not only its overall continuum shape but also the numerous metal lines exhibited in its ultraviolet spectrum. Results. TheoSSA is in operation and contains presently a variety of SEDs for DA-type white dwarfs. It will be extended in the near future and can host SEDs of all primary and secondary flux standards. The spectral analysis of G191−B2B has shown that our hydrostatic models reproduce the observations best at Teff =60 000 +/- 2000K and log g=7.60 +/- 0.05.We newly identified Fe vi, Ni vi, and Zn iv lines. For the first time, we determined the photospheric zinc abundance with a logarithmic mass fraction of −4.89 (7.5 × solar). The abundances of He (upper limit), C, N, O, Al, Si, O, P, S, Fe, Ni, Ge, and Sn were precisely determined. Upper abundance limits of about 10% solar were derived for Ti, Cr, Mn, and Co. Conclusions. The TheoSSA database of theoretical SEDs of stellar flux standards guarantees that the flux calibration of all astronomical data and cross-calibration between different instruments can be based on the same models and SEDs calculated with different model-atmosphere codes and are easy to compare.

Observatory↗
NASA NTRS1988

A re-evaluation of finite-element models and stress-intensity factors for surface cracks emanating from stress concentrations

A re-evaluation of the 3-D finite-element models and methods used to analyze surface crack at stress concentrations is presented. Previous finite-element models used by Raju and Newman for surface and corner cracks at holes were shown to have ill-shaped elements at the intersection of the hole and crack boundaries. These ill-shaped elements tended to make the model too stiff and, hence, gave lower stress-intensity factors near the hole-crack intersection than models without these elements. Improved models, without these ill-shaped elements, were developed for a surface crack at a circular hole and at a semi-circular edge notch. Stress-intensity factors were calculated by both the nodal-force and virtual-crack-closure methods. Both methods and different models gave essentially the same results. Comparisons made between the previously developed stress-intensity factor equations and the results from the improved models agreed well except for configurations with large notch-radii-to-plate-thickness ratios. Stress-intensity factors for a semi-elliptical surface crack located at the center of a semi-circular edge notch in a plate subjected to remote tensile loadings were calculated using the improved models. The ratio of crack depth to crack length ranged form 0.4 to 2; the ratio of crack depth to plate thickness ranged from 0.2 to 0.8; and the ratio of notch radius to the plate thickness ranged from 1 to 3. The models had about 15,000 degrees-of-freedom. Stress-intensity factors were calculated by using the nodal-force method.

Tan, P. W.↗
NASA NTRS1993

Contact-impact analysis of deployable space systems

This paper presents simulation results of a flexible, deployable space system with multiple contact-impact events, using three different modeling approaches. In the first two approaches, the equations of motion of the system are written in terms of generalized coordinates of each component, with modal coordinates used to account for flexibility in the components. Two different methods are used to model contact-impact events in the systems. In the momentum balance method, as contact-impact occurs, momentum equations of the systems are solved for the jump discontinuity in velocity of generalized coordinates. Integration of the system is restarted with the new initial conditions from solutions of momentum equations. In the logical function method, a virtual spring connecting the two contact points is assumed. The virtual spring is active only when impact occurs. No special treatment is needed in the integration algorithm using this method. In the third approach, instead of using modal coordinates, a finite element based formulation for the transient dynamics of mechanical system is used to model flexibility of the components. The logical function method is used to model contact-impact in the third approach.

Wu, Shih-Chin↗
NASA NTRS1988

Providing structural modules with self-integrity monitoring

With the advent of complex space structures (i.e., U.S. Space Station), the need for methods for remotely detecting structural damage will become greater. Some of these structures will have hundreds of individual structural elements (i.e., strut members). Should some of them become damaged, it could be virtually impossible to detect it using visual or similar inspection techniques. The damage of only a few individual members may or may not be a serious problem. However, should a significant number of the members be damaged, a significant problem could be created. The implementation of an appropriate remote damage detection scheme would greatly reduce the likelihood of a serious problem related to structural damage ever occurring. This report presents the results of the research conducted on remote structural damage detection approaches and the related mathematical algorithms. The research was conducted for the Small Business Innovation and Research (SBIR) Phase 2 National Aeronautics and Space Administration (NASA) Contract NAS7-961.

Walton, W. B.↗
NASA NTRS1993

surf3d: A 3-D finite-element program for the analysis of surface and corner cracks in solids subjected to mode-1 loadings

A computer program, surf3d, that uses the 3D finite-element method to calculate the stress-intensity factors for surface, corner, and embedded cracks in finite-thickness plates with and without circular holes, was developed. The cracks are assumed to be either elliptic or part eliptic in shape. The computer program uses eight-noded hexahedral elements to model the solid. The program uses a skyline storage and solver. The stress-intensity factors are evaluated using the force method, the crack-opening displacement method, and the 3-D virtual crack closure methods. In the manual the input to and the output of the surf3d program are described. This manual also demonstrates the use of the program and describes the calculation of the stress-intensity factors. Several examples with sample data files are included with the manual. To facilitate modeling of the user's crack configuration and loading, a companion program (a preprocessor program) that generates the data for the surf3d called gensurf was also developed. The gensurf program is a three dimensional mesh generator program that requires minimal input and that builds a complete data file for surf3d. The program surf3d is operational on Unix machines such as CRAY Y-MP, CRAY-2, and Convex C-220.

Raju, I. S.↗
NASA NTRS1998

Global Qualitative Flow-Path Modeling for Local State Determination in Simulation and Analysis

For qualitative modeling and analysis, a general qualitative abstraction of power transmission variables (flow and effort) for elements of flow paths includes information on resistance, net flow, permissible directions of flow, and qualitative potential is discussed. Each type of component model has flow-related variables and an associated internal flow map, connected into an overall flow network of the system. For storage devices, the implicit power transfer to the environment is represented by "virtual" circuits that include an environmental junction. A heterogeneous aggregation method simplifies the path structure. A method determines global flow-path changes during dynamic simulation and analysis, and identifies corresponding local flow state changes that are effects of global configuration changes. Flow-path determination is triggered by any change in a flow-related device variable in a simulation or analysis. Components (path elements) that may be affected are identified, and flow-related attributes favoring flow in the two possible directions are collected for each of them. Next, flow-related attributes are determined for each affected path element, based on possibly conflicting indications of flow direction. Spurious qualitative ambiguities are minimized by using relative magnitudes and permissible directions of flow, and by favoring flow sources over effort sources when comparing flow tendencies. The results are output to local flow states of affected components.

Malin, Jane T.↗
NASA NTRS1998

Superelement Analysis of Tile-Reinforced Composite Armor

Super-elements can greatly improve the computational efficiency of analyses of tile-reinforced structures such as the hull of the Composite Armored Vehicle. By taking advantage of the periodicity in this type of construction, super-elements can be used to simplify the task of modeling, to virtually eliminate the time required to assemble the stiffness matrices, and to reduce significantly the analysis solution time. Furthermore, super-elements are fully transferable between analyses and analysts, so that they provide a consistent method to share information and reduce duplication. This paper describes a methodology that was developed to model and analyze large upper hull components of the Composite Armored Vehicle. The analyses are based on two types of superelement models. The first type is based on element-layering, which consists of modeling a laminate by using several layers of shell elements constrained together with compatibility equations. Element layering is used to ensure the proper transverse shear deformation in the laminate rubber layer. The second type of model uses three-dimensional elements. Since no graphical pre-processor currently supports super-elements, a special technique based on master-elements was developed. Master-elements are representations of super-elements that are used in conjunction with a custom translator to write the superelement connectivities as input decks for ABAQUS.

Davila, Carlos G.↗
NASA NTRS2018

Simulating the Clamped Tapered Beam Specimen Under Quasi-Static and Fatigue Loading Using Floating Node Method

As part of the NASA Advanced Composites Project (ACP), a sub-element has been designed to provide validation data for progressive damage analysis models. The clamped tapered beam is a cross-ply laminated composite specimen designed to validate the simulation of the onset of matrix cracks and their interaction with delaminations, including delamination migration. A tapered geometry was used to localize the first damage occurrence in the tapered region, without prescribing an initial crack. The boundary and loading conditions were chosen to favor delamination growth and subsequent migration after the first damage occurrence. The typical sequence of events consists of a matrix crack located at the tapered region, leading to delamination onset, followed by delamination growth and subsequent delamination migration to a different interface via a dominant matrix crack. The Clamped Tapered Beam (CTB) was tested in both quasi-static and fatigue regimes. The results obtained are used in this study to assess and validate a methodology based on the Floating Node Method (FNM) implemented as an Extended Interface Element. In this methodology, quasi-static and fatigue damage formation and development are modeled by combining FNM to represent crack networks, with Directional Cohesive Zone Elements (DCZE) and Virtual Crack Closure Technique (VCCT), respectively. Qualitatively, the methodology is capable of predicting the sequence of events and overall failure morphology. Quantitatively, the simulation results generally bound the experimental data, based on the range of the characterization data used. In this paper, the results from quasi static and fatigue simulations are compared and correlated with experimental data.

Seshadri, B. R.↗