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 109 records · Page 6

Fatigue Life Methodology for Tapered Composite Flexbeam Laminates

The viability of a method for determining the fatigue life of composite rotor hub flexbeam laminates using delamination fatigue characterization data and a geometric non-linear finite element (FE) analysis was studied. Combined tension and bending loading was applied to non-linear tapered flexbeam laminates with internal ply drops. These laminates, consisting of coupon specimens cut from a full-size S2/E7T1 glass-epoxy flexbeam were tested in a hydraulic load frame under combined axial-tension and transverse cyclic bending. The magnitude of the axial load remained constant and the direction of the load rotated with the specimen as the cyclic bending load was applied. The first delamination damage observed in the specimens occurred at the area around the tip of the outermost ply-drop group. Subsequently, unstable delamination occurred by complete delamination along the length of the specimen. Continued cycling resulted in multiple delaminations. A 2D finite element model of the flexbeam was developed and a geometrically non-linear analysis was performed. The global responses of the model and test specimens agreed very well in terms of the transverse displacement. The FE model was used to calculate strain energy release rates (G) for delaminations initiating at the tip of the outer ply-drop area and growing toward the thick or thin regions of the flexbeam, as was observed in the specimens. The delamination growth toward the thick region was primarily mode 2, whereas delamination growth toward the thin region was almost completely mode 1. Material characterization data from cyclic double-cantilevered beam tests was used with the peak calculated G values to generate a curve predicting fatigue failure by unstable delamination as a function of the number of loading cycles. The calculated fatigue lives compared well with the test data.

Murri, Gretchen B.↗

A novel post-processing method for progressive failure analysis of brittle composite compression

Finite element analysis of brittle materials in axial compression typically uses element deletion to allow continued global deformation post-element-failure. However, element deletion produces cyclic load-displacement curves that underestimate energy absorption and are not representative of a continuum system. Two key observations support the conclusion that results from an appropriately discretized model can be an adequate representation of a continuum system. Specifically, the frequency of the oscillations in the load-displacement curve is directly dependent upon element length in the loading direction, and the peak amplitudes of oscillations are mesh size independent. A method of post-processing the analysis results, by connecting the peak amplitudes of oscillations, is proposed and applied to a series of continuous carbon fiber composite crush tubes. The load-displacement curve, stable crushing load, and specific energy absorption of the post-processed results compare well to an experimental study of crush tubes with similar layups.

Materials Science↗

A block iterative finite element algorithm for numerical solution of the steady-state, compressible Navier-Stokes equations

An iterative method for numerically solving the time independent Navier-Stokes equations for viscous compressible flows is presented. The method is based upon partial application of the Gauss-Seidel principle in block form to the systems of the nonlinear algebraic equations which arise in construction of finite element (Galerkin) models approximating solutions of fluid dynamic problems. The continuous cubic element on triangles is employed for function approximation. Computational results for a free shear flow at Re = 1000 indicate significant achievement of economy in iterative convergence rate over finite element and finite difference models which employ the customary time dependent equations and symptotic time marching procedure to steady solution. Numerical results are in excellent agreement with those obtained for the same test problem employing time marching finite element and finite difference solution techniques.

Cooke, C. H.↗

Algebraic surface design and finite element meshes

Some of the techniques are summarized which are used in constructing C sup 0 and C sup 1 continuous meshes of low degree, implicitly defined, algebraic surface patches in three dimensional space. These meshes of low degree algebraic surface patches are used to construct accurate computer models of physical objects. These meshes are also used in the finite element simulation of physical phenomena (e.g., heat dissipation, stress/strain distributions, fluid flow characteristics) required in the computer prototyping of both the manufacturability and functionality of the geometric design.

Bajaj, Chandrajit L.↗

An implicit finite element method for high speed flows

A fast algorithm is presented for constructing continuous lines, consisting of element sides, on general unstructured two-dimensional triangular meshes. The lines must pass through each node of the mesh once and only once. The discussion focuses on the use of these lines in a relaxation method for the solution of the equation system arising from an implicit algorithm for the solution of two-dimensional Euler and Navier-Stokes equations on general unstructured grids. It is also shown that the method can be used for solving three-dimensional Navier-Stokes equations on a grid composed of both structured and unstructured regions.

Hassan, O.↗

Bubble Dynamics on a Heated Surface

In this work, we study steady and oscillatory thermocapillary and natural convective flows generated by a bubble on a heated solid surface. The interaction between gas and vapor bubbles with the surrounding fluid is of interest for both space and ground-based processing. A combined numerical-experimental approach is adopted here. The temperature field is visualized using Mach-Zehnder and/or Wollaston Prism Interferometry and the flow field is observed by a laser sheet flow visualization technique. A finite element numerical model is developed which solves the transient two-dimensional continuity, momentum, and energy equations and includes the effects of temperature-dependent surface tension and bubble surface deformation. Below the critical Marangoni number, the steady state low-g and 1-g temperature and velocity fields predicted by the finite element model are in excellent agreement with both the visualization experiments in our laboratory and recently published experimental results in the literature. Above the critical Marangoni number, the model predicts an oscillatory flow which is also closely confirmed by experiments. It is shown that the dynamics of the oscillatory flow are directly controlled by the thermal and hydrodynamic interactions brought about by combined natural and thermocapillary convection. Therefore, as numerical simulations show, there are considerable differences between the 1-g and low-g temperature and flow fields at both low and high Marangoni numbers. This has serious implications for both materials processing and fluid management in space.

Kassemi, M.↗

Review of fatigue and fracture research at NASA Langley Research Center

Most dynamic components in helicopters are designed with a safe-life constant-amplitude testing approach that has not changed in many years. In contrast, the fatigue methodology in other industries has advanced significantly in the last two decades. Recent research at the NASA Langley Research Center and the U.S. Army Aerostructures Directorate at Langley are reviewed relative to fatigue and fracture design methodology for metallic components. Most of the Langley research was directed towards the damage tolerance design approach, but some work was done that is applicable to the safe-life approach. In the areas of testing, damage tolerance concepts are concentrating on the small-crack effect in crack growth and measurement of crack opening stresses. Tests were conducted to determine the effects of a machining scratch on the fatigue life of a high strength steel. In the area of analysis, work was concentrated on developing a crack closure model that will predict fatigue life under spectrum loading for several different metal alloys including a high strength steel that is often used in the dynamic components of helicopters. Work is also continuing in developing a three-dimensional, finite-element stress analysis for cracked and uncracked isotropic and anisotropic structures. A numerical technique for solving simultaneous equations called the multigrid method is being pursued to enhance the solution schemes in both the finite-element analysis and the boundary element analysis. Finally, a fracture mechanics project involving an elastic-plastic finite element analysis of J-resistance curve is also being pursued.

Everett, Richard A., Jr.↗

Full-Scale Direct Numerical Simulation of Two- and Three-Dimensional Instabilities and Rivulet Formulation in Heated Falling Films

A thin film draining on an inclined plate has been studied numerically using finite element method. Three-dimensional governing equations of continuity, momentum and energy with a moving boundary are integrated in an arbitrary Lagrangian Eulerian frame of reference. Kinematic equation is solved to precisely update interface location. Rivulet formation based on instability mechanism has been simulated using full-scale computation. Comparisons with long-wave theory are made to validate the numerical scheme. Detailed analysis of two- and three-dimensional nonlinear wave formation and spontaneous rupture forming rivulets under the influence of combined thermocapillary and surface-wave instabilities is performed.

Krishnamoorthy, S.↗

A Nonlinear Theory for Smart Composite Structures

The paper discusses the following: (1) Development of a completely coupled thermo-piezoelectric-mechanical theory for the analysis of composite shells with segmented and distributed piezoelectric sensor/actuators and shape memory alloys. The higher order displacement theory will be used to capture the transverse shear effects in anisotropic composites. The original theory will be modified to satisfy the stress continuity at ply interfaces. (2) Development of a finite element technique to implement the mathematical model. (3) Investigation of the coupled structures/controls interaction problem to study the complex trade-offs associated with the coupled problem.

Chattopadhyay, Aditi↗

A class of hybrid finite element methods for electromagnetics: A review

Integral equation methods have generally been the workhorse for antenna and scattering computations. In the case of antennas, they continue to be the prominent computational approach, but for scattering applications the requirement for large-scale computations has turned researchers' attention to near neighbor methods such as the finite element method, which has low O(N) storage requirements and is readily adaptable in modeling complex geometrical features and material inhomogeneities. In this paper, we review three hybrid finite element methods for simulating composite scatterers, conformal microstrip antennas, and finite periodic arrays. Specifically, we discuss the finite element method and its application to electromagnetic problems when combined with the boundary integral, absorbing boundary conditions, and artificial absorbers for terminating the mesh. Particular attention is given to large-scale simulations, methods, and solvers for achieving low memory requirements and code performance on parallel computing architectures.

Volakis, J. L.↗

Fluid-Structure Interaction Simulations of the ASPIRE SR03 Supersonic Parachute Flight Test

Research into parachute performance continues to be a source of significant investment from the National Aeronautics and Space Administration to mitigate risks and to enable a variety of exploration missions, including landing on Mars as well as returning to Earth. The cost of flight tests to certify any changes to the current state-of-the-art parachute designs limits the development of next generation parachute systems. Fluid-structure interaction simulations could help accelerate this process once validated. The Launch, Ascent, and Vehicle Aerodynamics team is developing the capability to perform such fluid-structure interaction simulations by coupling a higher-order Cartesian immersed boundary computational fluid dynamics solver with adaptive mesh refinement to a finite element structural dynamics solver in space and time. We continue the effort to validate this tool with the Advanced Supersonic Parachute Inflation Research Experiments SR03 flight test featuring a strengthened parachute akin to the Mars 2020 mission that landed the Perseverance rover on Mars, and a higher freestream dynamic pressure prior to inflation. The effect of the flow conditions’ angle of attack and of the initial parachute shape are quantified. The impact of relaxing modeling assumptions with regards to radial stiffeners on the parachute canopy is also investigated. Results demonstrate improvements in agreement with the pull force recorded during the SR03 flight test as the initial conditions of the flow and parachute are brought closer to those experienced in flight, and further improved when the radial stiffener modeling assumptions are relaxed.

ESM↗

Fluid-Structure Interaction Simulations of the ASPIRE SR03 Supersonic Parachute Flight Test

Research into parachute performance continues to be a source of significant investment from the National Aeronautics and Space Administration to mitigate risks and to enable a variety of exploration missions, including landing on Mars as well as returning to Earth. The cost of flight tests to certify any changes to the current state-of-the-art parachute designs limits the development of next generation parachute systems. Fluid-structure interaction simulations could help accelerate this process once validated. The Launch, Ascent, and Vehicle Aerodynamics team is developing the capability to perform such fluid-structure interaction simulations by coupling a higher-order Cartesian immersed boundary computational fluid dynamics solver with adaptive mesh refinement to a finite element structural dynamics solver in space and time. We continue the effort to validate this tool with the Advanced Supersonic Parachute Inflation Research Experiments SR03 flight test featuring a strengthened parachute akin to the Mars 2020 mission that landed the Perseverance rover on Mars, and a higher freestream dynamic pressure prior to inflation. The effect of the flow conditions’ angle of attack and of the initial parachute shape are quantified. The impact of relaxing modeling assumptions with regards to radial stiffeners on the parachute canopy is also investigated. Results demonstrate improvements in agreement with the pull force recorded during the SR03 flight test as the initial conditions of the flow and parachute are brought closer to those experienced in flight, and further improved when the radial stiffener modeling assumptions are relaxed.

ESM↗

Improved finite-element methods for rotorcraft structures

An overview of the research directed at improving finite-element methods for rotorcraft airframes is presented. The development of a modification to the finite element method which eliminates interelement discontinuities is covered. The following subject areas are discussed: geometric entities, interelement continuity, dependent rotational degrees of freedom, and adaptive numerical integration. This new methodology is being implemented as an anisotropic, curvilinear, p-version, beam, shell, and brick finite element program.

Hinnant, Howard E.↗

Continuously bounds-preserving discontinuous Galerkin methods for hyperbolic conservation laws

For finite element approximations of transport phenomena, it is often necessary to apply a form of limiting to ensure that the discrete solution remains well-behaved and satisfies physical constraints. However, these limiting procedures are typically performed at discrete nodal locations, which is not sufficient to ensure the robustness of the scheme when the solution must be evaluated at arbitrary locations (e.g., for adaptive mesh refinement, remapping in arbitrary Lagrangian–Eulerian solvers, overset meshes, etc.). In this work, a novel limiting approach for discontinuous Galerkin methods is presented which ensures that the solution is continuously bounds-preserving (i.e., across the entire solution polynomial) for any arbitrary choice of basis, approximation order, and mesh element type. Through a modified formulation for the constraint functionals, the proposed approach requires only the solution of a single spatial scalar minimization problem per element for which a highly efficient numerical optimization procedure is presented. Here, the efficacy of this approach is shown in numerical experiments by enforcing continuous constraints in high-order unstructured discontinuous Galerkin discretizations of hyperbolic conservation laws, ranging from scalar transport with maximum principle preserving constraints to compressible gas dynamics with positivity-preserving constraints.

97 MATHEMATICS AND COMPUTING↗

A penalty finite element algorithm for parabolic flow problems

The thin-layer simplification of the two-dimensional Navier-Stokes equations for steady viscous flow are developed using an order of magnitude analysis. A space marching finite element solution algorithm is developed, wherein the first order continuity effects are enforced as a penalty function differential constraint. Numerical results are presented to document accuracy and convergence features of the numerical solution algorithm.

Baker, A. J.↗

Evaluation of on-line pulse control for vibration suppression in flexible spacecraft

On-line pulse control for vibration suppression in a flexible spacecraft was evaluated. A continuous beam vs. a truss was modeled. A linear finite element model was used to determine the truss characteristics. Control issues outlined are ED pulse actuator development, pseudo pulse algorithm development, and large nonlinear simulation problems.

Bekey, G. A.↗

A Domain-Decomposed A-ϕ Formulation Based on Lagrange Multipliers for Low-Frequency Problems

A domain-decomposed A-ϕ formulation based on Lagrange multipliers is proposed to simulate low-frequency elec- tromagnetic problems. This method partitions the computational domain into smaller subdomains, allowing each subdomain to be independently formulated using Lagrange multipliers as Dirichlet boundary conditions, while ensuring continuity of the fields across the interfaces. A mixed finite element method, utilizing both vector and scalar basis functions, is employed to discretize the formulation, resulting in a global system to be solved. The proposed method is validated using TEAM Problem 7 at 50 Hz, demonstrating its effectiveness in handling complex geometries and addressing the low-frequency breakdown issues commonly encountered in traditional finite element methods.

Hossain, Amzad↗

The Design of Mechanically Compatible Fasteners for Human Mandible Reconstruction

Mechanically compatible fasteners for use with thin or weakened bone sections in the human mandible are being developed to help reduce large strain discontinuities across the bone/implant interface. Materials being considered for these fasteners are a polyetherertherketone (PEEK) resin with continuous quartz or carbon fiber for the screw. The screws were designed to have a shear strength equivalent to that of compact/trabecular bone and to be used with a conventional nut, nut plate, or an expandable shank/blind nut made of a ceramic filled polymer. Physical and finite element models of the mandible were developed in order to help select the best material fastener design. The models replicate the softer inner core of trabecular bone and the hard outer shell of compact bone. The inner core of the physical model consisted of an expanding foam and the hard outer shell consisted of ceramic particles in an epoxy matrix. This model has some of the cutting and drilling attributes of bone and may be appropriate as an educational tool for surgeons and medical students. The finite element model was exercised to establish boundary conditions consistent with the stress profiles associated with mandible bite forces and muscle loads. Work is continuing to compare stress/strain profiles of a reconstructed mandible with the results from the finite element model. When optimized, these design and fastening techniques may be applicable, not only to other skeletal structures, but to any composite structure.

Roberts, Jack C.↗