Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boundary 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 91 records · Page 5

Computational solid mechanics (finite elements and boundary elements) - Present status and future directions

Computational finite-element and boundary-element methods are reviewed, and their application to the mechanics of solids is discussed. Stability conditions for general FEMs are considered in addition to the use of least-order, stable, invariant, or hybrid/mixed isoparametric elements as alternatives to the displacement-based isoparametric elements. The use of symbolic manipulation, adaptive mesh refinement, transient dynamic response, and boundary-element methods for linear elaslticity and finite-strain problems of inelastic materials are also discussed.

Atluri, S. N.↗

Development of an integrated BEM approach for hot fluid structure interaction

A comprehensive boundary element method is presented for transient thermoelastic analysis of hot section Earth-to-Orbit engine components. This time-domain formulation requires discretization of only the surface of the component, and thus provides an attractive alternative to finite element analysis for this class of problems. In addition, steep thermal gradients, which often occur near the surface, can be captured more readily since with a boundary element approach there are no shape functions to constrain the solution in the direction normal to the surface. For example, the circular disc analysis indicates the high level of accuracy that can be obtained. In fact, on the basis of reduced modeling effort and improved accuracy, it appears that the present boundary element method should be the preferred approach for general problems of transient thermoelasticity.

Dargush, G. F.↗

A Hybrid Boundary Element-Finite Volume Method for Unsteady Transonic Airfoil Flows

A hybrid boundary element finite volume method for unsteady transonic flow computation has been developed. In this method, the unsteady Euler equations in a moving frame of reference are solved in a small embedded domain (inner domain) around the airfoil using an implicit finite volume scheme. The unsteady full-potential equation, written in the same frame of reference and in the form of the Poisson equation. is solved in the outer domain using the integral equation boundary element method to provide the boundary conditions for the inner Euler domain. The solution procedure is a time-accurate stepping procedure, where the outer boundary conditions for the inner domain are updated using the integral equation -- boundary element solution over the outer domain. The method is applied to unsteady transonic flows around the NACA0012 airfoil undergoing pitching oscillation and ramp motion. The results are compared with those of an implicit Euler equation solver, which is used throughout a large computational domain, and experimental data.

Hu, Hong↗

Three-dimensional analysis of surface crack-Hertzian stress field interaction

The results are presented of a stress intensity factor analysis of semicircular surface cracks in the inner raceway of an engine bearing. The loading consists of a moving spherical Hertzian contact load and an axial stress due to rotation and shrink fit. A 3-D linear elastic Boundary Element Method code was developed to perform the stress analysis. The element library includes linear and quadratic isoparametric surface elements. Singular quarter point elements were employed to capture the square root displacement variation and the inverse square root stress singularity along the crack front. The program also possesses the capability to separate the whole domain into two subregions. This procedure enables one to solve nonsymmetric fracture mechanics problems without having to separate the crack surfaces a priori. A wide range of configuration parameters was investigated. The ratio of crack depth to bearing thickness was varied from one-sixtieth to one-fifth for several different locations of the Hertzian load. The stress intensity factors for several crack inclinations were also investigated. The results demonstrate the efficiency and accuracy of the Boundary Element Method. Moreover, the results can provide the basis for crack growth calculations and fatigue life prediction.

Ballarini, R.↗

Second-order design sensitivities for linear elastic problems

The subject of this paper is the efficient and accurate determination of second-order design sensitivities in elastic bodies. The approach being carried out here is the direct differentiation of the governing derivative boundary element method (DBEM) formulation of the problem. Second-order sensitivities of boundary stresses are obtained in an elegant manner. A numerical implementation of the method is carried out with isoparametric quadratic boundary elements and numerical results are presented for several sample problems. Considerable savings in computational effort for an optimization procedure is possible through the use of efficiently determined accurate second-order design sensitivities.

Zhang, Qing↗

Conforming versus non-conforming boundary elements in three-dimensional elastostatics

A critical comparison of two basic formulations in three-dimensional elastostatics, using conforming and nonconforming boundary elements, is presented. The basic structure of the boundary element method is developed. The peculiarities that both types of boundary elements present in relation to the numerical implementation are discussed. Through selected examples, key issues such as the computational advantages and disadvantages of both formulations, mesh discretization and accuracy questions, and optimal location of the collocation nodes in the case of nonconforming elements are addressed. It is shown that conforming elements are able to produce more accurate results than nonconforming ones, with substantial economy in the final size of the system equations.

Manolis, G. D.↗

Dynamic response of lined tunnels by boundary elements

The dynamic stress concentration manifested around a lined cylindrical tunnel buried in an infinitely extending linear elastic or viscoelastic medium due to the passage of transient disturbances was investigated. Plane strain is assumed to hold and the transient disturbances can be of any arbitrary time variation. The boundary element method formulated in the Laplace transform domain is employed. Isoparametric boundary elements are used in the discretization of the liner and tunnel surfaces. Viscoelastic material behavior can be readily obtained from the linear elastic case in the Laplace transform domain through the use of the correspondence principle. The transient solution is recovered by numerical inversion of the solution obtained in the transformed domain.

Manolis, G. D.↗

A direct approach to finding unknown boundary conditions in steady heat conduction

The capability of the boundary element method (BEM) in determining thermal boundary conditions on surfaces of a conducting solid where such quantities are unknown was demonstrated. The method uses a non-iterative direct approach in solving what is usually called the inverse heat conduction problem (IHCP). Given any over-specified thermal boundary conditions such as a combination of temperature and heat flux on a surface where such data is readily available, the algorithm computes the temperature field within the object and any unknown thermal boundary conditions on surfaces where thermal boundary values are unavailable. A two-dimensional, steady-state BEM program was developed and was tested on several simple geometries where the analytic solution was known. Results obtained with the BEM were in excellent agreement with the analytic values. The algorithm is highly flexible in treating complex geometries, mixed thermal boundary conditions, and temperature-dependent material properties and is presently being extended to three-dimensional and unsteady heat conduction problems. The accuracy and reliability of this technique was very good but tended to deteriorate when the known surface conditions were only slightly over-specified and far from the inaccessible surface.

Martin, Thomas J.↗

Ares I-X Launch Abort System, Crew Module, and Upper Stage Simulator Vibroacoustic Flight Data Evaluation, Comparison to Predictions, and Recommendations for Adjustments to Prediction Methodology and Assumptions

The National Aeronautics and Space Administration (NASA) Constellation Program (CxP) has identified a series of tests to provide insight into the design and development of the Crew Launch Vehicle (CLV) and Crew Exploration Vehicle (CEV). Ares I-X was selected as the first suborbital development flight test to help meet CxP objectives. The Ares I-X flight test vehicle (FTV) is an early operational model of CLV, with specific emphasis on CLV and ground operation characteristics necessary to meet Ares I-X flight test objectives. The in-flight part of the test includes a trajectory to simulate maximum dynamic pressure during flight and perform a stage separation of the Upper Stage Simulator (USS) from the First Stage (FS). The in-flight test also includes recovery of the FS. The random vibration response from the ARES 1-X flight will be reconstructed for a few specific locations that were instrumented with accelerometers. This recorded data will be helpful in validating and refining vibration prediction tools and methodology. Measured vibroacoustic environments associated with lift off and ascent phases of the Ares I-X mission will be compared with pre-flight vibration predictions. The measured flight data was given as time histories which will be converted into power spectral density plots for comparison with the maximum predicted environments. The maximum predicted environments are documented in the Vibroacoustics and Shock Environment Data Book, AI1-SYS-ACOv4.10 Vibration predictions made using statistical energy analysis (SEA) VAOne computer program will also be incorporated in the comparisons. Ascent and lift off measured acoustics will also be compared to predictions to assess whether any discrepancies between the predicted vibration levels and measured vibration levels are attributable to inaccurate acoustic predictions. These comparisons will also be helpful in assessing whether adjustments to prediction methodologies are needed to improve agreement between the predicted and measured flight data. Future assessment will incorporate hybrid methods in VAOne analysis (i.e., boundary element methods, BEM and finite element methods, FEM). These hybrid methods will enable the ability to import NASTRAN models providing much more detailed modeling of the underlying beams and support structure of the ARES 1-X test vehicle. Measured acoustic data will be incorporated into these analyses to improve correlation for additional post flight analysis.

Smith, Andrew↗

An improved boundary force method for analyzing cracked anisotropic materials

The Boundary Force Method (BFM), a form of indirect boundary element method, is used to analyze composite laminates with cracks. The BFM uses the orthotropic elasticity solution for a concentrated horizontal and vertical force and a moment applied at a point in a cracked, infinite sheet as the fundamental solution. The necessary stress functions for this fundamental solution were formulated using the complex variable theory of orthotropic elasticity. The current method is an improvement over a previous method using only forces and no moment. The improved method was verified by comparing it to accepted solutions for a finite-width, center-crack specimen subjected to uniaxial tension. Four graphite/epoxy laminates were used: (0 + or - 45/90)sub s, (0), (+ or - 45)sub s, and (+ or - 30)sub s. The BFM results agreed well with accepted solutions. Convergence studies showed that with the addition of the moment in the fundamental solution, the number of boundary elements required for a converged solution was significantly reduced. Parametric studies were done for two configurations for which no orthotropic solutions are currently available; a single edge crack and an inclined single edge crack.

Tan, Paul W.↗

An improved boundary force method for analysing cracked anisotropic materials

The Boundary Force Method (BFM), a form of indirect boundary element method, is used to analyze composite laminates with cracks. The BFM uses the orthotropic elasticity solution for a concentrated horizontal and vertical force and a moment applied at a point in a cracked, infinite sheet as the fundamental solution. The necessary stress functions for this fundamental solution were formulated using the complex variables theory of orthotropic elasticity. The current method is an improvement over a previous method using only forces and no moment. The improved method was verified by comparing it to accepted solutions for a finite-width, center-crack specimen subjected to uniaxial tension. Four graphite/epoxy laminates were used: (0 + or - 45/90)sub s, (0), (+ or - 45)sub s, and (+ or - 30)sub s. The BFM results agreed well with accepted solutions. Convergence studies showed that with the addition of the moment in the fundamental solution, the number of boundary elements required for a converged solution was significantly reduced. Parametric studies were done for two configurations for which no orthotropic solutions are currently available; a single edge crack and an inclined single edge crack.

Tan, P. W.↗

Probabilistic Structural Analysis Program

NASA/NESSUS 6.2c is a general-purpose, probabilistic analysis program that computes probability of failure and probabilistic sensitivity measures of engineered systems. Because NASA/NESSUS uses highly computationally efficient and accurate analysis techniques, probabilistic solutions can be obtained even for extremely large and complex models. Once the probabilistic response is quantified, the results can be used to support risk-informed decisions regarding reliability for safety-critical and one-of-a-kind systems, as well as for maintaining a level of quality while reducing manufacturing costs for larger-quantity products. NASA/NESSUS has been successfully applied to a diverse range of problems in aerospace, gas turbine engines, biomechanics, pipelines, defense, weaponry, and infrastructure. This program combines state-of-the-art probabilistic algorithms with general-purpose structural analysis and lifting methods to compute the probabilistic response and reliability of engineered structures. Uncertainties in load, material properties, geometry, boundary conditions, and initial conditions can be simulated. The structural analysis methods include non-linear finite-element methods, heat-transfer analysis, polymer/ceramic matrix composite analysis, monolithic (conventional metallic) materials life-prediction methodologies, boundary element methods, and user-written subroutines. Several probabilistic algorithms are available such as the advanced mean value method and the adaptive importance sampling method. NASA/NESSUS 6.2c is structured in a modular format with 15 elements.

Pai, Shantaram S.↗

Three-dimensional inelastic analysis for hot section components, BEST 3D code

The goal is the development of an alternative stress analysis tool, distinct from the finite element method, applicable to the engineering analysis of gas turbine engine structures. The boundary element method was selected for this development effort on the basis of its already demonstrated applicability to a variety of geometries and problem types characteristic of gas turbine engine components. Major features of the BEST3D computer program are described, and some of the significant developments carried out as part of the Inelastic Methods Contract are outlined.

Wilson, Raymond B.↗

International Congress on Recent Developments in Air- and Structure-Borne Sound and Vibration, 2nd, Auburn Univ., AL, Mar. 4-6, 1992, Proceedings. Vols. 1-3

The papers presented in this volume cover a variety of topics, including sound intensity, structural intensity, modal analysis and synthesis, statistical energy analysis and energy methods, passive and active damping, and boundary element methods. Attention is also given to diagnostics and condition monitoring, material characterization and nondestructive evaluation, active noise and vibration control, sound radiation and scattering, and finite element analysis.

Crocker, Malcolm J.↗

COMGEN-BEM: Boundary element model generation for composite materials micromechanical analysis

Composite Model Generation-Boundary Element Method (COMGEN-BEM) is a program developed in PATRAN command language (PCL) which generates boundary element models of continuous fiber composites at the micromechanical (constituent) scale. Based on the entry of a few simple parameters such as fiber volume fraction and fiber diameter, the model geometry and boundary element model are generated. In addition, various mesh densities, material properties, fiber orientation angles, loads, and boundary conditions can be specified. The generated model can then be translated to a format consistent with a boundary element analysis code such as BEST-CMS.

Goldberg, Robert K.↗

Program Helps Generate Boundary-Element Mathematical Models

Composite Model Generation-Boundary Element Method (COM-GEN-BEM) computer program significantly reduces time and effort needed to construct boundary-element mathematical models of continuous-fiber composite materials at micro-mechanical (constituent) scale. Generates boundary-element models compatible with BEST-CMS boundary-element code for anlaysis of micromechanics of composite material. Written in PATRAN Command Language (PCL).

Goldberg, R. K.↗

Free-vibration analysis of three-dimensional solids by BEM

This paper discusses the calculation of natural frequencies and mode shapes of three-dimensional solids, using the boundary element method. The metod developed is based on the use of particular integrals of the elastic equilibrium equations, and employs only real arithmetic. It represents an extension to three dimensions of work previously carried out for two-dimensional problems by Ahmad and Banerjee (1986) and Nardini and Brebbia (1982). The method has been incorporated, for multiregion analysis, in the BEST3D (Boundary Element Solution Technology, Three-Dimensional) computer program. Results of the boundary element calculations are compared with analytical, finite element, and experimental results.

Wilson, R. B.↗

Development of an integrated BEM approach for hot fluid structure interaction

Significant progress was made toward the goal of developing a general purpose boundary element method for hot fluid-structure interaction. For the solid phase, a boundary-only formulation was developed and implemented for uncoupled transient thermoelasticity in two dimensions. The elimination of volume discretization not only drastically reduces required modeling effort, but also permits unconstrained variation of the through-the-thickness temperature distribution. Meanwhile, for the fluids, fundamental solutions were derived for transient incompressible and compressible flow in the absence of the convective terms. Boundary element formulations were developed and described. For the incompressible case, the necessary kernal functions, under transient and steady-state conditions, were derived and fully implemented into a general purpose, multi-region boundary element code. Several examples were examined to study the suitability and convergence characteristics of the various algorithms.

Dargush, Gary F.↗