Search NASA⌕ Search

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.

At least 667 records · Page 37

Finite element model of SCOLE laboratory configuration

The Spacecraft Control Laboratory Experiment (SCOLE) is defined by element properties: material constants; mast, reflector, rigid links as beam elements; cable as bar element; and space shuttle as very stiff beam. Two boundary conditions are modeled: suspended (6 degrees of freedom for all joints except the top of the cable) and cantilevered cables (shuttle platform fixed in all degrees of freedon). Calculations include stiffness and mass matrices, initial stresses, static displacements and reactions, and eigensolutions.

Lee, Beth↗

Finite element method for non-linear forced vibrations of circular plates

Geometric non-linearities for large amplitude free and forced vibrations of circular plates are investigated. In-plane displacement and in-plane inertia are included in the formulation. The finite element method is used. An harmonic force matrix for non-linear forced vibration analysis is introduced and derived. Various out-of-plane and in-plane boundary conditions are considered. The relations of amplitude and frequency ratio for different boundary conditions and various load conditions are presented.

Mei, Chuh↗

A velocity-pressure integrated, mixed interpolation, Galerkin finite element method for high Reynolds number laminar flows

A velocity-pressure integrated, mixed interpolation, Galerkin finite element method for the Navier-Stokes equations is presented. In the method, the velocity variables were interpolated using complete quadratic shape functions and the pressure was interpolated using linear shape functions. For the two dimensional case, the pressure is defined on a triangular element which is contained inside the complete biquadratic element for velocity variables; and for the three dimensional case, the pressure is defined on a tetrahedral element which is again contained inside the complete tri-quadratic element. Thus the pressure is discontinuous across the element boundaries. Example problems considered include: a cavity flow for Reynolds number of 400 through 10,000; a laminar backward facing step flow; and a laminar flow in a square duct of strong curvature. The computational results compared favorable with those of the finite difference methods as well as experimental data available. A finite elememt computer program for incompressible, laminar flows is presented.

Kim, Sang-Wook↗

Probabilistic structural analysis methods and applications

An advanced algorithm for simulating the probabilistic distribution of structural responses due to statistical uncertainties in loads, geometry, material properties, and boundary conditions is reported. The method effectively combines an advanced algorithm for calculating probability levels for multivariate problems (fast probability integration) together with a general-purpose finite-element code for stress, vibration, and buckling analysis. Application is made to a space propulsion system turbine blade for which the geometry and material properties are treated as random variables.

Cruse, T. A.↗

Compact finite volume methods for the diffusion equation

An approach to treating initial-boundary value problems by finite volume methods is described, in which the parallel between differential and difference arguments is closely maintained. By using intrinsic geometrical properties of the volume elements, it is possible to describe discrete versions of the div, curl, and grad operators which lead, using summation-by-parts techniques, to familiar energy equations as well as the div curl = 0 and curl grad = 0 identities. For the diffusion equation, these operators describe compact schemes whose convergence is assured by the energy equations and which yield both the potential and the flux vector with second order accuracy. A simplified potential form is especially useful for obtaining numerical results by multigrid and alternating direction implicit (ADI) methods. The treatment of general curvilinear coordinates is shown to result from a specialization of these general results.

Rose, Milton E.↗

The effect of coatings and liners on heat transfer in a dry shaft-bush tribosystem

The temperatures due to frictional heating within a solid lubricated or coated journal bearing were analyzed by using a finite element method. A solid model of the shaft-bush tribocontact was generated with an eight-node, three-dimensional, first-order isoparametric heat-transfer element and the Patran solid modeler software. The Patmar (Patran-Marc) translator was used to help develop the Marc-based finite element program for the system; this software was used on the Cray X-MP supercomputer to perform a finite element analysis of the contact. The analysis was performed for various liner materials, for thin, hard, wear-resistant coated bearings, and for different geometries and thermal cooling boundary conditions. The analyses indicated that thermal conductivity of the liner or coating material is the most vital thermal parameter that controls the interface temperature. In addition to design variations, the proximity of the cooling source to the heat-flux-generating interface is critically important to the temperature control in the system.

Ghosh, Mihir K.↗

Analysis of bolt-loaded elliptical holes in laminated composite joints

The feasibility of changing the bolt shape from circular to elliptical in order to increase the joint strength was investigated. An analytical method using a cosine series to represent the bearing stress was derived, and a boundary collocation method was used to determine the unknown coefficients of this cosine series. Stresses at the hole edge predicted by the analytical method agreed very well with finite element solutions. Failure analyses of joints in two different laminates were performed, and each laminate exhibited a different joint failure mode. Results demonstrate that the joint strength of both laminates can be improved substantially by changing the bolt shape to elliptical. The joint that failed in a bearing mode showed a greater strength increase compared to the joint that failed in a shearing mode.

Wang, J. T.↗

Probabilistic assessment of space nuclear propulsion system nozzle

In assessing the reliability of a space nuclear propulsion system (SNPS) nozzle, uncertainties associated with the following design parameters were considered: geometry, boundary conditions, material behavior, and thermal and pressure loads. A preliminary assessment of the reliability was performed using NESSUS (Numerical Evaluation of Stochastic Structures Under Stress), a finite-element computer code developed at the NASA Lewis Research Center. The sensitivity of the nozzle reliability to the uncertainties in the random variables was quantified. With respect to the effective stress, preliminary results showed that the nozzle spatial geometry uncertainties have the most significant effect at low probabilities whereas the inner wall temperature has the most significant effect at higher probabilities.

Shah, Ashwin R.↗

A global/local analysis method for treating details in structural design

A method for analyzing global/local behavior of plate and shell structures is described. In this approach, a detailed finite element model of the local region is incorporated within a coarser global finite element model. The local model need not be nodally compatible (i.e., need not have a one-to-one nodal correspondence) with the global model at their common boundary; therefore, the two models may be constructed independently. The nodal incompatibility of the models is accounted for by introducing appropriate constraint conditions into the potential energy in a hybrid variational formulation. The primary advantage of this method is that the need for transition modeling between global and local models is eliminated. Eliminating transition modeling has two benefits. First, modeling efforts are reduced since tedious and complex transitioning need not be performed. Second, errors due to the mesh distortion, often unavoidable in mesh transitioning, are minimized by avoiding distorted elements beyond what is needed to represent the geometry of the component. The method is applied reduced to a plate loaded in tension and transverse bending. The plate has a central hole, and various hole sixes and shapes are studied. The method is also applied to a composite laminated fuselage panel with a crack emanating from a window in the panel. While this method is applied herein to global/local problems, it is also applicable to the coupled analysis of independently modeled components as well as adaptive refinement.

Aminpour, Mohammad A.↗

Experimental Validation of FE/BEM Dynamic Strain Model Under Diffuse Acoustic Field Loading

Structural finite element (FE) models naturally output displacement or acceleration response data. However, they can also be used to compute stress, internal forces, and strain response. When coupled with a boundary element model (BEM) of the fluid surrounding the structure, a fully coupled analysis can be performed. Modeling a diffuse acoustic field in the BEM fluid provides an excitation like that found when the structure is placed in a reverberation chamber. Fully coupling the structural FE model to the acoustic BEM model provides a means to predict not only the acceleration response of the panel to diffuse field loading, but also the ability to predict the dynamic stress and strain response. This type of model has been available with current predictive tools, but experimental validation of the prediction of dynamic stress or strain is difficult to find. An aluminum panel was instrumented with accelerometers and strain gages and hung in a reverberation room and subjected to a diffuse acoustic field. This paper presents the comparison of the experimental and predicted results.

Tsoi, W. Ben↗

Automatic Incorporation of Flaws into Cad and FEA Models Based On Measured NDE Data

Automatic incorporation of non-destructive evaluation (NDE) data into CAD and finite element (FE) models for structural analysis is essential for accurate residual-life and strength analyses of composite structures. Currently, two significant but independent challenges exist for incorporating NDE data in such analyses: interpretation of the NDE signals to locate and characterize the structural defects and the incorporation of that information into structural analysis tools. The first challenge will be addressed through linear inversion of calculated synthetic responses generated by a simulation engine to evaluate flaw size, location, and other characteristics. The second challenge will be addressed by building on prior and parallel defect modeling efforts to integrate suitable flaw models into the CAD/finite-element representation of a specimen. Ultrasonic NDE information is usually represented using C-Scans of time-of-flight or intensity information. Other modalities such as thermography and computed tomography (CT) could also be used to characterize the structure. This information will first be projected into material coordinates and then aligned with the specimen’s finite element mesh. The NDE data can then be incorporated by inserting flaw geometry into the CAD and FE analysis model and automatic modification of local material properties and boundary conditions. This process will enable the use of structural analysis tools even for practical modalities such as ultrasound and thermography which may not directly give accurate-enough geometric information to fully define the flaw. The full flaw characteristics will be inferred using prior expert and modeling knowledge to categorize flaw signatures. Outputs for this project will include demonstrated methods for creating structural analysis geometric models that incorporate flaw characterization information; these models can then be used directly in popular commercial FEA software. In addition, the primary deliverable for the project will include a software that will be able to automatically incorporate different kinds of flaws into an existing FE models. The models will be generated in standard file formats to enable direct import into the analysis software.

Adarsh Krishnamurthy↗

NASTRAN analysis of Tokamak vacuum vessel using interactive graphics

Isoparametric quadrilateral and triangular elements were used to represent the vacuum vessel shell structure. For toroidally symmetric loadings, MPCs were employed across model boundaries and rigid format 24 was invoked. Nonsymmetric loadings required the use of the cyclic symmetry analysis available with rigid format 49. NASTRAN served as an important analysis tool in the Tokamak design effort by providing a reliable means for assessing structural integrity. Interactive graphics were employed in the finite element model generation and in the post-processing of results. It was felt that model generation and checkout with interactive graphics reduced the modelling effort and debugging man-hours significantly.

Miller, A.↗

Transonic flutter study of a wind-tunnel model of an arrow-wing supersonic transport

A 1/20-size, low-speed flutter model of the SCAT-15F complete airplane was tested on cables to simulate a near free-flying condition. Only the model wing and fuselage were flexible. Flutter boundaries were measured for a nominal configuration and a configuration with wing fins removed at Mach numbers M from 0.76 to 1.2. For both configurations, the transonic dip in the wing flutter dynamic pressure q boundary was relatively small and the minimum flutter q occurred near M = 0.92. Removing the wing fins increased the flutter q about 14 percent and changed the flutter mode from symmetric to antisymmetric. Vibration and flutter analyses were made using a finite-element structural representation and subsonic kernel-function aerodynamics. For the nominal configuration, the analysis (using calculated modal data) predicted the experimental flutter q levels within 10 percent but did not predict the correct flutter mode at the higher M. For the configuration without wing fins, the analysis predicted 16 to 36 percent unconservative (higher than experimental) flutter q levels and showed extreme sensitivity to mass representation details that affected wing tip mode shapes. For high subsonic M, empennage aerodynamics had a significant effect on the predicted flutter boundaries of several symmetric modes.

Ruhlin, C. L.↗

Tectonic motions in California inferred from very long baseline interferometry observations, 1980-1984

Using VLBI, three baselines of 150-300 km length within California, and three baselines of 1500 km length between California and Texas, measured between 1980 and 1984, have been fitted to a steady plane/strain rotation model. The observed rate of displacement across the southern San Andreas fault is in agreement with other contemporary rate estimates near 3 cm/yr. Discrepancies between these rates and indicators of relative plate motion near 6 cm/yr seem to be related to local departures of the San Andreas system from the expected strike of the plate boundary. Results support a regional distribution of elastic strain that is offset in about the same sense and amount as the Big Bend of the San Andreas fault, and finite element simulations suggest that tractions on the base of the crust that accompany mantle downswelling beneath the Transverse Ranges could contribute to these strain inhomogeneities.

Lyzenga, Gregory A.↗

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

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

Methodology for sensitivity analysis, approximate analysis, and design optimization in CFD for multidisciplinary applications

Fundamental equations of aerodynamic sensitivity analysis and approximate analysis for the two dimensional thin layer Navier-Stokes equations are reviewed, and special boundary condition considerations necessary to apply these equations to isolated lifting airfoils on 'C' and 'O' meshes are discussed in detail. An efficient strategy which is based on the finite element method and an elastic membrane representation of the computational domain is successfully tested, which circumvents the costly 'brute force' method of obtaining grid sensitivity derivatives, and is also useful in mesh regeneration. The issue of turbulence modeling is addressed in a preliminary study. Aerodynamic shape sensitivity derivatives are efficiently calculated, and their accuracy is validated on two viscous test problems, including: (1) internal flow through a double throat nozzle, and (2) external flow over a NACA 4-digit airfoil. An automated aerodynamic design optimization strategy is outlined which includes the use of a design optimization program, an aerodynamic flow analysis code, an aerodynamic sensitivity and approximate analysis code, and a mesh regeneration and grid sensitivity analysis code. Application of the optimization methodology to the two test problems in each case resulted in a new design having a significantly improved performance in the aerodynamic response of interest.

Taylor, Arthur C., III↗