Finite-element solution of the incompressible lubrication problem.
Finite element solution of incompressible lubrication problem by minimum principle for transient incompressible Reynolds equation with boundary conditions
Reddi, M. M.↗
SEARCH · Search NASA
Results for “Boundary 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.
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.↗
A two dimensional mesh verification algorithm
A two dimensional mesh verification algorithm is discussed which assumes that a finite element mesh is represented in a program by lists of data, i.e., vertex coordinates, element incidences, boundary data, which describe a collection of triangles. A set of conditions are placed on the triangles (i.e., on the list data) which ensure that the triangles tile a region and which also can be verified by the algorithm.
Simpson, R. B.↗
Implementation of a C-1 triangular element based on the P-version of the finite element method
The implementation of a computer code CONE (for C(1) continuity) based on the p-version of the finite element method is described. A hierarchic family of triangular finite elements of degree p 5 is used. This family enforces C(1)-continuity across interelement boundaries, and the code is applicable to fourth order partial differential equations in two independent variables, in particular to the biharmonic equation. Applications to several benchmark problems in plate bending are presented. Sample results are examined and compared with theoretical predictions. In particular the analysis of the bending of a rhombic plate shows a significant improvement over othr published results.
Wang, D. W.↗
The role of damage-softened material behavior in the fracture of composites and adhesives
Failure mechanisms of materials under very high strains experienced at and ahead of the crack tip such as formation, growth, and interaction of microvoids in ductile materials, microcracks in brittle solids or crazes in polymers and adhesives are represented by one-dimensional, nonlinear stress-strain relations possessing different ways by which the material loses capacity to carry load up to fracture or total separation. A double cantilever beam (DCB) type specimen is considered. The nonlinear material is confined to a thin strip between the two elastic beams loaded by a wedge. The problem is first modeled as a beam on a nonlinear foundation. The pertinent equation is solved numerically as a two-point boundary value problem for both the stationary and the quasi-stationay propagating crack. A finite element model is then used to model the problem in more detail in order to assess the adequacy of the beam model for the reduction of experimental data to determine in-situ properties of the thin interlayer.
Ungsuwarungsri, T.↗
The effect of nonlinearities on the dynamic response of a large shuttle payload
The STS Centaur was designed to be a high energy upper stage for use with the Space Shuttle. Two versions were designed under development when the program was cancelled. The first version, designated G-prime, was for planetary missions. The second version, designated G, was to place spacecraft in geosynchronous orbit. As a part of the STS Centaur finite-element model verification effort, test articles of both versions were subjected to a series of static tests. In addition the Centaur G-prime test article was subjected to a series of dynamic tests including a modal survey. Both the static and dynamic tests showed that nonlinearities existed in the Centaur and its support system. The support system included flight-like latches. The nonlinearities were particularly apparent in tests that loaded the forward support structure of the Centaur. These test results were used to aid in the development of two improved finite-element models. The first was a linear model, while the second contained nonlinear elements at the boundaries. Results from both models were compared with the transient response obtained from a step-relaxation or twang test. The linear model was able to accurately match the low frequency response found in the test data. However, only the nonlinear model was able to match higher frequency response that was present in some of the test data. In addition the nonlinear model was able to predict other nonlinear behavior such as the dynamic jump that occurs in systems with nonlinear stiffness.
Sullivan, Timothy L.↗
The effect of nonlinearities on the dynamic response of a large Shuttle payload
The STS Centaur was designed to be a high energy upper stage for use with the Space Shuttle. Two versions were designed and under development when the program was cancelled. The first version, designated G-prime, was designed for planetary missions. The second version, designated G, was designed to place spacecraft in geosynchronous orbit. As a part of the STS Centaur finite-element model verification effort, test articles of both versions were subjected to a series of static tests. In addition the Centaur G-prime test article was subjected to a series of dynamic tests including a modal survey. Both the static and dynamic tests showed that nonlinearities existed in the Centaur and its support system. The support system included flight-like latches. The nonlinearities were particularly apparent in tests that loaded the forward support structure of the Centaur. These test results were used to aid in the development of two improved finite-element models. The first was a linear model, while the second contained nonlinear elements at the boundaries. Results from both models were compared with the transient response obtained from a step-relaxation or twang test. The linear model was able to accurately match the low frequency response found in the test data. However, only the nonlinear model was able to match higher frequency response that was present in some of the test data. In addition the nonlinear model was able to predict other nonlinear behavior such as the dynamic 'jump' that occurs in systems with nonlinear stiffness.
Sullivan, Timothy L.↗
Development of a model of space station solar array
Space structures, such as the space station solar arrays, must be extremely lightweight, flexible structures. Accurate prediction of the natural frequencies and mode shapes is essential for determining the structural adequacy of components, and designing a control system. The tension preload in the blanket of photovoltaic solar collectors, and the free/free boundary conditions of a structure in space, causes serious reservations on the use of standard finite element techniques of solution. In particular, a phenomena known as grounding, or false stiffening, of the stiffness matrix occurs during rigid body rotation. The grounding phenomena is examined in detail. Numerous stiffness matrices developed by others are examined for rigid body rotation capability, and found lacking. Various techniques are used for developing new stiffness matrices from the rigorous solutions of the differential equations, including the solution of the directed force problem. A new directed force stiffness matrix developed by the author provides all the rigid body capabilities for the beam in space.
Bosela, Paul A.↗
Grounding of space structures
Space structures, such as the Space Station solar arrays, must be extremely light-weight, flexible structures. Accurate prediction of the natural frequencies and mode shapes is essential for determining the structural adequacy of components, and designing a controls system. The tension pre-load in the 'blanket' of photovoltaic solar collectors, and the free/free boundary conditions of a structure in space, causes serious reservations on the use of standard finite element techniques of solution. In particular, a phenomenon known as 'grounding', or false stiffening, of the stiffness matrix occurs during rigid body rotation. This paper examines the grounding phenomenon in detail. Numerous stiffness matrices developed by others are examined for rigid body rotation capability, and found lacking. A force imbalance inherent in the formulations examined is the likely cause of the grounding problem, suggesting the need for a directed force formulation.
Bosela, P. A.↗
A new pre-loaded beam geometric stiffness matrix with full rigid body capabilities
Space structures, such as the Space Station solar arrays, must be extremely light-weight, flexible structures. Accurate prediction of the natural frequencies and mode shapes is essential for determining the structural adequacy of components, and designing a controls system. The tension pre-load in the 'blanket' of photovoltaic solar collectors, and the free/free boundary conditions of a structure in space, causes serious reservations on the use of standard finite element techniques of solution. In particular, a phenomenon known as 'grounding', or false stiffening, of the stiffness matrix occurs during rigid body rotation. The authors have previously shown that the grounding phenomenon is caused by a lack of rigid body rotational capability, and is typical in beam geometric stiffness matrices formulated by others, including those which contain higher order effects. The cause of the problem was identified as the force imbalance inherent in the formulations. In this paper, the authors develop a beam geometric stiffness matrix for a directed force problem, and show that the resultant global stiffness matrix contains complete rigid body mode capabilities, and performs very well in the diagonalization methodology customarily used in dynamic analysis.
Bosela, P. A.↗
A geometrically nonlinear analysis of interlaminar stresses in unsymmetrically laminated plates subjected to uniform thermal loading
An analytical study of interlaminar stresses in unsymmetrically laminated plates is presented. The study examines the linear elastic, large deflection response of square laminated composite plates subjected to uniform thermal loading. Both cross-ply and angle-ply, symmetric and unsymmetric, laminates are examined to evaluate the effects of mismatch between adjacent layers in elastic constants and coefficients of thermal expansion. A geometrically nonlinear kinematic description is used to predict the large out-of-plane (transverse) deflections. The nonlinear, three-dimensional boundary value problems are formulated from elasticity theory and approximate solutions are determined using the finite element method. A global/local analysis procedure is used to obtain improved free edge stress predictions. For the laminates and loading conditions considered, the results indicate that the out-of-plane deflections of the unsymmetric laminates reduce interlaminar shear stresses. These deflections also reduce interlaminar normal stresses in some laminates and increase these stresses for others. The results indicate that as the out-of-plane deflections become large, the differences in interlaminar stress predictions between linear and nonlinear theory can become quite large.
Norwood, D. S.↗
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.↗
Numerical Analysis of Temperature Gradients and Interface Shape During Directional Solidification of Al and Al-Cu Alloy Under Microgravity Conditions
Numerical modeling was undertaken to analyze the influence of radial thermal gradient on solid/liquid (s/1) interface shape and convection patterns during solidification of pure Al and Al-4 wt% Cu alloy. The objective of the numerical task was to predict the influence of convective velocity on an insoluble particle near a s/l interface. These predictions would then be used to define the minimum gravity level (g) required to investigate the fundamental physics of interaction between a particle and a s/I interface. To satisfy this objective, steady state calculations were performed for different gravity levels and orientations with the gravity vector. ne furnace configuration used in this analysis is the proposed International Space Station Furnace, Quench Module Insert (QMI) 1. Results from a thermal model of the furnace core were used as initial boundary conditions for solidification modeling. General model of binary alloy solidification was based on the finite element code FIDAP. It was found that for the worst case orientation of 90 degrees with the gravity vector and a g level of 10(exp -4)g(sub o) (g(sub o) = 9.8 m/s(exp 2)) the dominant forces acting on the particle would be the fundamental drag and interfacial forces.
Bune, Andris V.↗
Quantifying Square Membrane Wrinkle Behavior Using MITC Shell Elements
For future membrane based structures, quantified predictions of membrane wrinkling behavior in terms of amplitude, angle and wavelength are needed to optimize the efficiency and integrity of such structures, as well as their associated control systems. For numerical analyses performed in the past, limitations on the accuracy of membrane distortion simulations have often been related to the assumptions made while using finite elements. Specifically, this work demonstrates that critical assumptions include: effects of gravity. supposed initial or boundary conditions, and the type of element used to model the membrane. In this work, a 0.2 square meter membrane is treated as a structural material with non-negligible bending stiffness. Mixed Interpolation of Tensorial Components (MTTC) shell elements are used to simulate wrinkling behavior due to a constant applied in-plane shear load. Membrane thickness, gravity effects, and initial imperfections with respect to flatness were varied in numerous nonlinear analysis cases. Significant findings include notable variations in wrinkle modes for thickness in the range of 50 microns to 1000 microns, which also depend on the presence of an applied gravity field. However, it is revealed that relationships between overall strain energy density for cases with differing initial conditions are independent of assumed initial con&tions. In addition, analysis results indicate that the relationship between amplitude scale (W/t) and structural scale (L/t) is linear in the presence of a gravity field.
Jacobson, Mindy B.↗
Aeroelastic Simulation of Wind Tunnel Scale Aspect Ratio 13.5 Common Research Model with Trailing Edge Mini-Flaps
This paper presents a computation study of an aspect ratio 13.5 wind tunnel scale Common Research Model (CRM) with trailing edge mini-flaps using FUN3D. The aeroelastic simulations are conducted under free-air and inside wind tunnel conditions. The NASA Langley Research Center Transonic Dynamic Tunnel (TDT) is a closed circuit and slotted-test-section with a 16 ft ×16 ft test section. The structure mode shapes are extracted from the latest updated NASTRAN Finite-Element Model (FEM) model. A dynamic pressure sweep is performed to determine the control reversal boundary of the trailing edge mini-flaps. The wind tunnel wall interference effects on the mini-flap control derivatives are analyzed based on the simulations.
Aeroelastic↗
Additive manufacturing system and method having toolpath analysis
An additive manufacturing (AM) system is disclosed for constructing a three dimensional (3D) part with optimized orthotropy. The system combines an electronic processor which calculates an optimal set of physically achievable toolpaths to meet a given design objective, and a 3D printing direct ink write machine capable of printing inks with reinforcing particles that result in orthotropic materials. The electronic processor may combine a Domain of Interest subsystem that transforms a mathematical description of a desired part orthotropy to a plurality of guidepaths, a toolpath generate subsystem that develops a plurality of physically realizable toolpaths from those guidepaths with a minimum number of starts and stops, a finite element subsystem that computes spatially varying material orthotropy from those toolpaths and then solves a boundary value problem to determine a figure of merit for the design, and an optimization subsystem that uses that figure of merit to update the mathematical description of the part orthotropy to iteratively develop an improved part. The optimization subsystem also includes convergence criteria to indicate when toolpaths have been achieved that yield a sufficiently optimal part has been achieved. The toolpath generate can then output the final toolpaths, which are converted to suitable code that controls the motion of the 3D printer toolhead and allows the optimized 3D part to be manufactured.
King, Michael J.↗
Considerations of crack growth and plasticity in finite element analysis
A finite-element analysis was made of crack growth in a center-cracked specimen subjected to monotonically increasing load until the point of fast fracture. Since part of the specimen experienced unloading, the boundary value problem which was formulated was based upon incremental theory of plasticity. Experimental load and crack size records were utilized. Linear relations between plastic energy and crack growth were observed. Fracture toughness parameters, which were evaluated at the onset of unstable crack propagation from finite-element analysis, were in good agreement with those determined experimentally.
Lee, J. D.↗
A high-order Shifted Interface Method for Lagrangian shock hydrodynamics
Here, we present a new method for two-material Lagrangian hydrodynamics, which combines the Shifted Interface Method (SIM) with a high-order Finite Element Method. Our approach relies on an exact (or sharp) material interface representation, that is, it uses the precise location of the material interface. The interface is represented by the zero level-set of a continuous high-order finite element function that moves with the material velocity. This strategy allows to evolve curved material interfaces inside curved elements. By reformulating the original interface problem over a surrogate (approximate) interface, located in proximity of the true interface, the SIM avoids cut cells and the associated problematic issues regarding implementation, numerical stability, and matrix conditioning. Accuracy is maintained by modifying the original interface conditions using Taylor expansions. We demonstrate the performance of the proposed algorithms on established numerical benchmarks in one, two and three dimensions.
97 MATHEMATICS AND COMPUTING↗