Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuous 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 523 records · Page 29

Anisotropic yield models for lattice unit cell structures exploiting orthotropic symmetry

We report numerical homogenization enables efficient computational analysis and design of multiscale structures made of micro-architected materials including lattice unit cells. However, predicting yield is nontrivial because it requires accurate and efficient predictions of the maximum stress inside the homogenized unit cells. To address this challenge, we develop a macroscale anisotropic yield function for micro-architected materials. The yield function depends on the three-dimensional macroscale stress state and the parameters describing a family of micro-architectures, such as the radii of the struts in a lattice unit cell. To ensure accuracy, we determine the maximum stress using three-dimensional continuum finite-element analysis. To ensure efficiency, we construct surrogate models from the aforementioned high-fidelity results for yield prediction. To reduce simulation costs and surrogate modeling complexity, we leverage orthotropic symmetry commonly found in lattice unit cells. In this paper, we provide a thorough presentation of group representations and the systematic procedure to exploit orthotropic domain symmetry in homogenization and surrogate modeling. We illustrate the use of linear homogenization results to predict yield in specific unit cells without further simulations. We furthermore show that surrogate modeling presents a viable option for anisotropic yield prediction in a continuously parametrized family of micro-architectures. More specifically, despite the dimensionality and degree of nonlinearity in the maximum micro von Mises stress within the unit cell, the surrogate models can predict it with less than 5% error at least 90% of the time. Moreover, the largest under-prediction error, which is more critical than the over-prediction error, is typically less than 10%.

42 ENGINEERING↗

The finite element duct eigenvalue problem - An improved formulation with Hermitian elements and no-flow condensation

Hermitian elements are used in a finite element solution for the eigenvalue problem in lined ducts with flow. These elements give significantly greater accuracy for reduced dimensionality when compared with Lagrangian elements. Spurious mode generation associated with the Lagrangian formulation is eliminated. A dramatic improvement in the ratio of the number of reliable eigenvalues to the total number of computed eigenvalues is effected by the use of a condensation scheme based on the no-flow eigenvectors. Results are presented for two dimensional and axisymmetric ducts. In the axisymmetric case good resolution is obtained even for high order, high frequency modes by the use of continuously graded meshes.

Astley, R. J.↗

Buckling and vibration analysis of laminated panels using VICONOPT

The analysis aspects of the 23,000-line FORTRAN program VIPASA with constraints and optimization (VICONOPT) are described. Overall stiffness matrices assembled from the earlier exact vibration and instability of plate assemblies including shear and anisotropy (VIPASA) flat plate stiffnesses are optionally coupled by Langrangian multipliers to find critical buckling loads, or natural frequencies of undamped vibration, of prismatic assemblies of anisotropic flat plates with arbitrarily located point supports or simple transverse supporting frames. The longitudunal continuity of typical wing and fuselage panels is closely approximated because the solutions are for the infinitely long structure obtained by repeating a bay and its supports longitudinally. Any longitudinally invariant in-plane plate stresses are permitted, and very rapid solutions are guaranteed by numerous refinements, including multilevel substructuring and a method for repetitive cross sections that is exact for regular polygons used to represent cylinders. Modal displacements and stresses in or between plies of laminated plates are calculated and plotted, with values being recovered at all nodes of substructures. Comparison with usual approximate finite-element methods confirms that, for comparably converged solutions, VICONOPT is typically between 100 and 10(exp 4) times faster.

Kennedy, David↗

Development of Fuzzy Logic and Soft Computing Methodologies

Our earlier research on computing with words (CW) has led to a new direction in fuzzy logic which points to a major enlargement of the role of natural languages in information processing, decision analysis and control. This direction is based on the methodology of computing with words and embodies a new theory which is referred to as the computational theory of perceptions (CTP). An important feature of this theory is that it can be added to any existing theory - especially to probability theory, decision analysis, and control - and enhance the ability of the theory to deal with real-world problems in which the decision-relevant information is a mixture of measurements and perceptions. The new direction is centered on an old concept - the concept of a perception - a concept which plays a central role in human cognition. The ability to reason with perceptions perceptions of time, distance, force, direction, shape, intent, likelihood, truth and other attributes of physical and mental objects - underlies the remarkable human capability to perform a wide variety of physical and mental tasks without any measurements and any computations. Everyday examples of such tasks are parking a car, driving in city traffic, cooking a meal, playing golf and summarizing a story. Perceptions are intrinsically imprecise. Imprecision of perceptions reflects the finite ability of sensory organs and ultimately, the brain, to resolve detail and store information. More concretely, perceptions are both fuzzy and granular, or, for short, f-granular. Perceptions are f-granular in the sense that: (a) the boundaries of perceived classes are not sharply defined; and (b) the elements of classes are grouped into granules, with a granule being a clump of elements drawn together by indistinguishability, similarity. proximity or functionality. F-granularity of perceptions may be viewed as a human way of achieving data compression. In large measure, scientific progress has been, and continues to be, driven by a quest to progress from perceptions to measurements. Pursuit of this aim has led to brilliant successes. But alongside the successes stand problems whose solutions are not in sight. Representative of such problems is the problem of automation of driving in city traffic. In this case, as in many others, what can be done with ease by humans - without any measurements and a computations - is an intractable task for machines.

Zadeh, L. A.↗

Integrated Force Method for Indeterminate Structures

Two methods of solving indeterminate structural-mechanics problems have been developed as products of research on the theory of strain compatibility. In these methods, stresses are considered to be the primary unknowns (in contrast to strains and displacements being considered as the primary unknowns in some prior methods). One of these methods, denoted the integrated force method (IFM), makes it possible to compute stresses, strains, and displacements with high fidelity by use of modest finite-element models that entail relatively small amounts of computation. The other method, denoted the completed Beltrami Mitchell formulation (CBMF), enables direct determination of stresses in an elastic continuum with general boundary conditions, without the need to first calculate displacements as in traditional methods. The equilibrium equation, the compatibility condition, and the material law are the three fundamental concepts of the theory of structures. For almost 150 years, it has been commonly supposed that the theory is complete. However, until now, the understanding of the compatibility condition remained incomplete, and the compatibility condition was confused with the continuity condition. Furthermore, the compatibility condition as applied to structures in its previous incomplete form was inconsistent with the strain formulation in elasticity.

Hopkins, Dale A.↗

A Multiscale Progressive Failure Modeling Methodology for Composites that Includes Fiber Strength Stochastics

A multiscale modeling methodology was developed for continuous fiber composites that incorporates a statistical distribution of fiber strengths into coupled multiscale micromechanics/finite element (FE) analyses. A modified two-parameter Weibull cumulative distribution function, which accounts for the effect of fiber length on the probability of failure, was used to characterize the statistical distribution of fiber strengths. A parametric study using the NASA Micromechanics Analysis Code with the Generalized Method of Cells (MAC/GMC) was performed to assess the effect of variable fiber strengths on local composite failure within a repeating unit cell (RUC) and subsequent global failure. The NASA code FEAMAC and the ABAQUS finite element solver were used to analyze the progressive failure of a unidirectional SCS-6/TIMETAL 21S metal matrix composite tensile dogbone specimen at 650 degC. Multiscale progressive failure analyses were performed to quantify the effect of spatially varying fiber strengths on the RUC-averaged and global stress-strain responses and failure. The ultimate composite strengths and distribution of failure locations (predominately within the gage section) reasonably matched the experimentally observed failure behavior. The predicted composite failure behavior suggests that use of macroscale models that exploit global geometric symmetries are inappropriate for cases where the actual distribution of local fiber strengths displays no such symmetries. This issue has not received much attention in the literature. Moreover, the model discretization at a specific length scale can have a profound effect on the computational costs associated with multiscale simulations.models that yield accurate yet tractable results.

Generalized Method of Cells↗

Mechanical properties of freestanding few-layer graphene/boron nitride/polymer heterostacks investigated with local and non-local techniques

van der Waals two-dimensional materials and heterostructures combined with polymer films continue to attract research attention to elucidate their functionality and potential applications. This study presents the fabrication and mechanical testing of 2D material heterostacks, consisting of few-layer boron nitride and graphene heterostructures synthesized via chemical vapor deposition, capped with a polymethyl methacrylate layer and suspended across ∼200 μm wide trenches using a combined wet–dry transfer method. The mechanical characterization of the heterostacks was performed using two independent approaches: (a) non-local testing with a custom-built tensile testing platform and (b) local load–displacement testing employing atomic force microscopy probes, complemented by finite element simulations. Both approaches provided new results, which are in good agreement with each other. Overall, our findings offer new insights into a combined load capacity in complex multi-material two-dimensional systems, and can contribute to advancing micro and nano-scale device designs and implementations.

Lespasio, Marcus↗

Calculation of residual principal stresses in CVD boron on carbon filaments

A three-dimensional finite element model of the chemical vapor deposition (CVD) of boron on a carbon substrate (B/C) is developed. The model includes an expansion of the boron after deposition due to atomic rearrangement and includes creep of the boron and carbon. Curves are presented to show how the principal residual stresses and the filament elongation vary as the parameters defining deposition strain and creep are varied. The calculated results are compared with experimental axial residual stress and elongation measurements made on B/C filaments. This comparison requires that for good agreement between calculated and experimental results, the deposited boron must continue to expand after deposition, and that the build-up of residual stresses is limited by significant boron and carbon creep rates.

Behrendt, D. R.↗

High energy modular bootstrap, global symmetries and defects

We derive Cardy-like formulas for the growth of operators in different sectors of unitary 2 dimensional CFT in the presence of topological defect lines by putting an upper and lower bound on the number of states with scaling dimension in the interval [Δ – δ, Δ + δ] for large Δ at fixed δ. Consequently we prove that given any unitary modular invariant 2D CFT symmetric under finite global symmetry G (acting faithfully), all the irreducible representations of G appear in the spectra of the untwisted sector; the growth of states is Cardy like and proportional to the “square” of the dimension of the irrep. In the Schwarzian limit, the result matches onto that of JT gravity with a bulk gauge theory. If the symmetry is non-anomalous, the result applies to any sector twisted by a group element. For c > 1, the statements are true for Virasoro primaries. Furthermore, the results are applicable to large c CFTs. We also extend our results for the continuous U(1) group.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Integrated analysis and design of thick composite structures for optimal passive damping characteristics

The development of novel composite mechanics for the analysis of damping in composite laminates and structures and the more significant results of this effort are summarized. Laminate mechanics based on piecewise continuous in-plane displacement fields are described that can represent both intralaminar stresses and interlaminar shear stresses and the associated effects on the stiffness and damping characteristics of a composite laminate. Among other features, the mechanics can accurately model the static and damped dynamic response of either thin or thick composite laminates, as well as, specialty laminates with embedded compliant damping layers. The discrete laminate damping theory is further incorporated into structural analysis methods. In this context, an exact semi-analytical method for the simulation of the damped dynamic response of composite plates was developed. A finite element based method and a specialty four-node plate element were also developed for the analysis of composite structures of variable shape and boundary conditions. Numerous evaluations and applications demonstrate the quality and superiority of the mechanics in predicting the damped dynamic characteristics of composite structures. Finally, additional development was focused on the development of optimal tailoring methods for the design of thick composite structures based on the developed analytical capability. Applications on composite plates illustrated the influence of composite mechanics in the optimal design of composites and the potential for significant deviations in the resultant designs when more simplified (classical) laminate theories are used.

Saravanos, D. A.↗

3D printed carbon fiber reinforced carbon as an energy efficient alternative to graphite for EFAS tooling

As Electric Field Assisted Sintering (EFAS) gains more industrial acceptance and use, it becomes more important to develop more efficient means to implement this technology. To this aim, 3D printed continuous carbon fiber reinforced carbon (CCC) was manufactured and fabricated into tooling for EFAS systems as an alternative to traditional graphite tooling. The impact of fiber orientation on the thermal and electrical properties of the CCC was characterized. Sample material was sintered in Tokai G535 graphite tooling, under common processing conditions and compared with CCC tooling. There was nearly 50 % energy savings compared to graphite while maintaining equivalent sample density and microstructure plus keeping ram temperatures 39 % cooler. This is due to spatial control of generated heat and thermal diffusivity within the molds, by means of fiber orientation anisotropy. Finite element modeling of the tooling design supported the experimental results as well as displays the effect of optimization of this 3D printed CCC material.

36 MATERIALS SCIENCE↗

Open-Mode Debonding Analysis of Curved Sandwich Panels Subjected to Heating and Cryogenic Cooling on Opposite Faces

Increasing use of curved sandwich panels as aerospace structure components makes it vital to fully understand their thermostructural behavior and identify key factors affecting the open-mode debonding failure. Open-mode debonding analysis is performed on a family of curved honeycomb-core sandwich panels with different radii of curvature. The curved sandwich panels are either simply supported or clamped, and are subjected to uniform heating on the convex side and uniform cryogenic cooling on the concave side. The finite-element method was used to study the effects of panel curvature and boundary condition on the open-mode stress (radial tensile stress) and displacement fields in the curved sandwich panels. The critical stress point, where potential debonding failure could initiate, was found to be at the midspan (or outer span) of the inner bonding interface between the sandwich core and face sheet on the concave side, depending on the boundary condition and panel curvature. Open-mode stress increases with increasing panel curvature, reaching a maximum value at certain high curvature, and then decreases slightly as the panel curvature continues to increase and approach that of quarter circle. Changing the boundary condition from simply supported to clamped reduces the magnitudes of open-mode stresses and the associated sandwich core depth stretching.

Ko, William L.↗

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗

Dynamic fracture toughnesses of reaction-bonded silicon nitride

The room-temperature dynamic fracture response of reaction-bonded silicon nitride is investigated using a hybrid experimental-numerical procedure. In this procedure, experimentally determined crack velocities are utilized to drive a dynamic finite-element code or dynamic finite-difference code in its generation mode in order to extract numerically the dynamic stress intensity factor of the fracturing specimen. Results show that the dynamic fracture toughness vs crack velocity relations of the two reaction-bonded silicon nitrides do not follow the general trend in those relations of brittle polymers and steel. A definite slow crack velocity during the initial phase of dynamic crack propagation is observed in reaction-bonded silicon nitride, which results in a nonunique dynamic fracture toughness vs crack velocity relation. In addition, it is found that a propagating crack will continue to propagate under a static stress intensity factor substantially lower than K(IC).

Kobayashi, A. S.↗

Static shape determination and control for large space structures. I - The flexible beam. II - A large space antenna

A method for determining and controlling the shape of large, continuous space structures by means of discrete or pointwise observations and control devices is presented. The general linear boundary value problem satisfied by a one-dimensional shape function is defined, and the existence of solutions is studied. The static shape control problems for one-dimensional systems with and without rigid body modes and the static shape estimation problem are presented and solved. Eigenfunction expansions are presented which provide approximations to the algorithm terms when the associated Green's function is not known. An integral operator approach is applied to the multidimensional static problem, and the results are illustrated with a finite element model of the disk of a large space antenna which assumes no rigid body modes. It is shown that the shape control algorithm must be modified for systems with rigid body modes.

Weeks, C. J.↗

Progressive Damage Analysis of Post-Buckled Stiffened Panels under Static Compressive Loading

The validation of progressive damage analysis methods requires experimental data that captures the initiation of damage, its evolution with continued loading, and the morphology of damage throughout its evolution. In support of the NASA Advanced Composites Project, several validation specimens of increasing complexity have been tested and characterized in order to provide the data required to assess the capabilities of several progressive damage analysis methods for fiber-reinforced composite materials. The analytical capabilities of interest for this validation program include networks of interacting matrix cracks and delaminations in post-buckled stiffened structures. A series of single stringer compression specimens was tested and characterized at NASA Langley Research Center. The experimental data from these tests are compared with progressive damage finite element analysis results. Test/analysis comparisons are made in terms of load/displacement histories, damage morphology, and damage size versus load metrics.

Frank A Leone↗

Optimization for Load Alleviation of Truss-Braced Wing Aircraft with Variable Camber Continuous Trailing Edge Flap

This paper focuses on load alleviation optimization for a high aspect ratio truss braced wing (TBW) aircraft. The TBW aircraft model is based on the Subsonic Ultra Green Aircraft Research (SUGAR) concept developed by Boeing, with the wing structures of the model modified to include a novel aerodynamic control surface known as the Variable Camber Continuous Trailing Edge Flap (VCCTEF). The purpose of the study is to investigate the effectiveness of a Performance Adaptive Aeroelastic Wing (PAAW) technology, specifically the VCCTEF, for alleviating load on the TBW wing during flight maneuver. The specific flight maneuver under consideration in this study is a 2.5g pull-up maneuver. Constrained gradient-based optimization is conducted to tailor the deflections of the VC- CTEF such that bending moment along the wing is minimized at the 2.5g pull-up flight condition. Aerodynamic modeling for this study is conducted using a vortex-lattice method code called Vorlax. A non-linear finite element analysis (FEA) method is constructed for analyzing the structural deformation and resulting bending moment along the wing of the aircraft with the inclusion of effects from tension-stiffening due to axial loading in the truss. This study is the first phase of several, and involves optimization of a rigid wing aircraft for preliminary analysis. Future studies will incorporate flexible wing structures with aeroelastic interactions and deformations. The results of this first phase positively demonstrate the potential of utilizing the novel control surface on modern aircraft wing designs for shaping control in order to provide load alleviation during flight maneuver.

Lebofsky, Sonia↗

Integrated Modeling Activities for the James Webb Space Telescope (JWST): Structural-Thermal-Optical Analysis

This is a continuation of a series of papers on modeling activities for JWST. The structural-thermal- optical, often referred to as "STOP", analysis process is used to predict the effect of thermal distortion on optical performance. The benchmark STOP analysis for JWST assesses the effect of an observatory slew on wavefront error. The paper begins an overview of multi-disciplinary engineering analysis, or integrated modeling, which is a critical element of the JWST mission. The STOP analysis process is then described. This process consists of the following steps: thermal analysis, structural analysis, and optical analysis. Temperatures predicted using geometric and thermal math models are mapped to the structural finite element model in order to predict thermally-induced deformations. Motions and deformations at optical surfaces are input to optical models and optical performance is predicted using either an optical ray trace or WFE estimation techniques based on prior ray traces or first order optics. Following the discussion of the analysis process, results based on models representing the design at the time of the System Requirements Review. In addition to baseline performance predictions, sensitivity studies are performed to assess modeling uncertainties. Of particular interest is the sensitivity of optical performance to uncertainties in temperature predictions and variations in metal properties. The paper concludes with a discussion of modeling uncertainty as it pertains to STOP analysis.

Johnston, John D.↗