Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Singularity”

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

Contact and crack problems for an elastic wedge

The contact and the crack problems for an elastic wedge of arbitrary angle are considered. The problem is reduced to a singular integral equation which, in the general case, may have a generalized Cauchy kernel. The singularities under the stamp as well as at the wedge apex were studied, and the relevant stress intensity factors are defined. The problem was solved for various wedge geometries and loading conditions. The results may be applicable to certain foundation problems and to crack problems in symmetrically loaded wedges in which cracks initiate from the apex.

Erdogan, F.↗

Comparison of some biased estimation methods (including ordinary subset regression) in the linear model

Ridge, Marquardt's generalized inverse, shrunken, and principal components estimators are discussed in terms of the objectives of point estimation of parameters, estimation of the predictive regression function, and hypothesis testing. It is found that as the normal equations approach singularity, more consideration must be given to estimable functions of the parameters as opposed to estimation of the full parameter vector; that biased estimators all introduce constraints on the parameter space; that adoption of mean squared error as a criterion of goodness should be independent of the degree of singularity; and that ordinary least-squares subset regression is the best overall method.

Sidik, S. M.↗

Calculations of unsteady turbulent boundary layers with flow reversal

The results are presented of a series of computational experiments aimed at studying the characteristics of time-dependent turbulent boundary layers with embedded reversed-flow regions. A calculation method developed earlier was extended to boundary layers with reversed flows for this purpose. The calculations were performed for an idealized family of external velocity distributions, and covered a range of degrees of unsteadiness. The results confirmed those of previous studies in demonstrating that the point of flow reversal is nonsingular in a time-dependent boundary layer. A singularity was observed to develop downstream of reversal, under certain conditions, accompanied by the breakdown of the boundary-layer approximations. A tentative hypothesis was advanced in an attempt to predict the appearance of the singularity, and is shown to be consistent with the calculated results.

Nash, J. F.↗

Computational aspects of maximum likelihood estimation and reduction in sensitivity function calculations

This paper discusses numerical aspects of computing maximum likelihood estimates for linear dynamical systems in state-vector form. Different gradient-based nonlinear programming methods are discussed in a unified framework and their applicability to maximum likelihood estimation is examined. The problems due to singular Hessian or singular information matrix that are common in practice are discussed in detail and methods for their solution are proposed. New results on the calculation of state sensitivity functions via reduced order models are given. Several methods for speeding convergence and reducing computation time are also discussed.

Gupta, N. K.↗

Small-scale yielding analysis of mixed mode plane-strain crack problems

The small-scale yielding analysis of an elastic-plastic body with a line crack under plane-strain conditions subject to combinations of Mode I and II loadings is examined. The analysis of the near-tip field follows the works of Hutchinson and Rice and Rosengren. Dominant singularity solutions governing the asymptotic behavior of the stresses and strains at the crack tip are obtained for the complete range of loadings between Mode I and II. The results of an accurate finite element technique, which imbeds the dominant singularity solutions, directly relates the near-tip behavior to the elastic stress intensity factors K(I) and K(II). Implications of this study to mixed mode fracture mechanics is also discussed, particularly with respect to the direction of crack initiation and the relation of fracture toughness under mixed modes to that in Mode I. Details of the mixed mode plastic zone sizes and shapes are also given.

Shih, C. F.↗

The problem of edge cracks in an infinite strip

The elastostatic plane problem of an infinite strip containing two symmetrically located internal cracks perpendicular to the boundary is formulated in terms of a singular integral equation with the derivative of the crack surface displacement as the density function. The solution of the problem is obtained for various crack geometries and for uniaxial tension applied to the strip away from the crack region. The limiting case of the edge cracks is then considered in some detail. The fundamental function of the integral equation is obtained and a numerical technique for solving the singular integral equations with this particular type of fundamental function which is characteristic of the edge cracks is described. The stress-intensity factor for the complete range of net ligament-to-width ratio is calculated. The results also include the solution of the edge crack problem in an elastic half plane.

Gupta, G. D.↗

The inclusion problem with a crack crossing the boundary

A solution is given to the problems of a crack in an elastic inclusion with one or both ends approaching and terminating at the interface, of two collinear cracks (one in the inclusion and one in the matrix), and of a crack crossing the interface. The problems are formulated in terms of a system of singular integral equations. In the second case, the dominant parts of the kernels become generalized Cauchy kernels giving rise to stress singularities of powers other than one over the square root. For this case stress intensity factors are defined and some detailed results are presented for various crack-inclusion geometries and material combinations.

Erdogan, F.↗

Pines nonsingular gravitational potential derivation, description and implementation

An engineering interpretation of and some minor corrections to previously completed work on a uniform representation of the gravitational potential and its derivatives (with emphasis on avoiding the usual singularity at the pole) were presented. The physical meaning of the variables was explained, the derivation of results was separated into smaller parts for easier reading, some additional recurrence relations for the derived Legendre polynomials were included and compared, and a computer program implementing this formulation was presented. Numerical experiments have shown that the use of this representation, besides removing the singularity, substantially increases the speed of the computation.

Spencer, J. L.↗

Fracture of composite panels

The fracture problem in panels consisting of periodically arranged load carrying and buffer strips of different materials is considered. The main emphasis is placed on the problem of a crack terminating at and crossing the interfaces and on the stress free end problem. The problem is formulated in terms of a system of singular integral equations, and numerical solutions are obtained for certain material combinations. With the study of possible crack propagation and delamination in mind, certain stress intensity factors are defined and calculated. A main result is that when the crack touches or intersects a bimaterial interface, the stress state has no longer the standard square root singularity, and, to study further propagation of the crack, the conventional fracture models need to be modified, or new models need to be developed.

Erdogan, F.↗

The effect of adhesive layer on crack propagation in laminates

The effect of the adhesive layer on crack propagation in composite materials is investigated. The composite medium consists of parallel load carrying laminates and buffer strips arranged periodically and bonded with thin adhesive layers. The strips, assumed to be isotropic and linearly elastic, contain symmetric cracks of arbitrary lengths located normal to the interfaces. Two problems are considered: (1) thin adhesive layers are approximated by uncoupled tension and shear springs distributed along the interfaces of the strips for which only the case of internal cracks can be treated rigorously; (2) broken laminates and the true singular behavior in the presence of the adhesive layer are studied. The adhesive is then treated as an isotropic, linearly elastic continuum. General expressions for field quantities are obtained in terms of infinite Fourier integrals. These expressions give a system of singular integral equations in terms of the crack surface displacement derivatives. By using appropriate quadrature formulas, the integral equations reduce to a system of linear algebraic equations which are solved numerically.

Gecit, M. R.↗

A lifting surface theory for the analysis of nonplanar lifting systems

A new nonlinear, nonplanar lifting surface theory is presented. The method is regarded as a lifting surface theory in that the effects of wing thickness are neglected, but none of the usual small perturbation assumptions inherent in most other lifting surface theories are made. The method represents nonplanar lifting systems by distributed vorticity, including the leading edge singular behavior characteristic of thin wings. The method is well suited to the computation of induced drag of nonplanar systems because leading edge suction is calculated from the leading edge singularity. The method has been used to compute the induced drag benefit of winglets (vortex diffusers), and the agreement with NASA experimental data is excellent.

Goldhammer, M. I.↗

Fracture and contact problems for an elastic wedge

The paper deals with the plane elastostatic contact problem for an infinite elastic wedge of arbitrary angle. The medium is loaded through a frictionless rigid wedge of a given symmetric profile. Using the Mellin transform formulation the mixed boundary value problem is reduced to a singular integral equation with the contact stress as the unknown function. With the application of the results to the fracture of the medium in mind, the main emphasis in the study has been on the investigation of the singular nature of the stress state around the apex of the wedge and on the determination of the contact pressure.

Erdogan, F.↗

Computation of the inviscid supersonic flow over an external axial corner

A second-order finite-difference procedure is used to evaluate the inviscid supersonic flowfield surrounding an external axial corner composed of swept planar compression surfaces and representing the inlets on existing high-speed aircraft. The governing partial differential equations in conservation-law form are hyperbolic with respect to the axial coordinate and are solved iteratively by means of MacCormack's algorithm. The procedure treats both the peripheral shock wave and vortical singularities as discontinuities. Numerical results are presented for two parametric studies regarding the effects on the flowfield of varying the free-stream Mach number and the leading edge sweep of the horizontal wedge. Results of parametric Mach number study agree with the Mach number independence principle in that as the Mach number increases, such characteristics as shock shape, cross-flow sonic line location, and vortical singularity position approach an asymptote.

Kutler, P.↗

Algorithms for computing the geopotential using a simple density layer

Several algorithms have been developed for computing the potential and attraction of a simple density layer. These are numerical cubature, Taylor series, and a mixed analytic and numerical integration using a singularity-matching technique. A computer program has been written to combine these techniques for computing the disturbing acceleration on an artificial earth satellite. A total of 1640 equal-area, constant surface density blocks on an oblate spheroid are used. The singularity-matching algorithm is used in the subsatellite region, Taylor series in the surrounding zone, and numerical cubature on the rest of the earth.

Morrison, F.↗

Sound propagation in choked ducts

The linearized equations describing the propagation of sound in variable area ducts containing flow are shown to be singular when the duct mean flow is sonic. The singularity is removed when previously ignored nonlinear terms are retained. The results of a numerical study, for the case of plane waves propagating in a one-dimensional converging-diverging duct, show that the sound field is adequately described by the linearized equations only when the axial mean flow Mach number at the duct throat M sub th 0.6. For M sub th 0.6, the numerical results showed that acoustic energy flux was not conserved. An attempt was made to extend the study to include the nonlinear behavior of the sound field. Meaningful results were not obtained due, primarily, to numerical difficulties.

Hersh, A. S.↗

The use of parabolic variations and the direct determination of stress intensity factors using the BIE method

Two advances in the numerical techniques of utilizing the BIE method are presented. The boundary unknowns are represented by parabolas over each interval which are integrated in closed form. These integrals are listed for easy use. For problems involving crack tip singularities, these singularities are included in the boundary integrals so that the stress intensity factor becomes just one more unknown in the set of boundary unknowns thus avoiding the uncertainties of plotting and extrapolating techniques. The method is applied to the problems of a notched beam in tension and bending, with excellent results.

Mendelson, A.↗

Supersonic flow about circular cones at large angles of attack - A floating discontinuity approach

The technique of floating shock fitting is adapted to the computation of the inviscid flowfield about circular cones in a supersonic free stream at angles of attack that exceed the cone half-angle. In those regions in which the governing conical equations are mixed elliptic-hyperbolic, the fully hyperbolic form is obtained by the addition of the temporal derivative. The resulting equations are applicable over the complete range of free-stream Mach numbers, angles of attack and cone half-angles for which the bow shock is attached. An explicit finite-difference algorithm is used to obtain the solution by an unsteady relaxation approach. The bow shock, embedded crossflow shock, and vortical singularity in the leeward symmetry plane are all treated as floating discontinuities in a fixed computational mesh. The method yields excellent results for the bow and embedded shocks, however, the solution in the leeward symmetry plane exhibits viscous-like effects and does not appear to adequately predict the behavior of the vortical singularity.

Daywitt, J.↗

Stress-intensity factors for a wedge-loaded edge crack in a semi-infinite strip

The problem of a semi infinite strip containing an edge crack is considered. It is assumed that the strip is loaded by a frictionless rigid wedge pressed into the crack. The resulting crack contact problem is formulated in terms of a system of singular integral equations. The behavior of the solution near the singular points is studied in detail. A series of numerical examples are given and the results are compared with those obtained by the method of boundary collocation and by the simple beam theory.

Erodogan, F.↗