Engineering topics

Belytschko, Ted

Publications and source records attributed to Belytschko, Ted.

At least 19 records

NASA NTRS1993

Elastic interactions of a fatigue crack with a micro-defect by the mixed boundary integral equation method

In this paper, the mixed boundary integral equation method is developed to study the elastic interactions of a fatigue crack and a micro-defect such as a void, a rigid inclusion or a transformation inclusion. The method of pseudo-tractions is employed to study the effect of a transformation inclusion. An enriched element which incorporates the mixed-mode stress intensity factors is applied to characterize the singularity at a moving crack tip. In order to evaluate the accuracy of the numerical procedure, the analysis of a crack emanating from a circular hole in a finite plate is performed and the results are compared with the available numerical solution. The effects of various micro-defects on the crack path and fatigue life are investigated. The results agree with the experimental observations.

Lua, Yuan J.↗

NASA NTRS1993

Adaptive methods for nonlinear structural dynamics and crashworthiness analysis

The objective is to describe three research thrusts in crashworthiness analysis: adaptivity; mixed time integration, or subcycling, in which different timesteps are used for different parts of the mesh in explicit methods; and methods for contact-impact which are highly vectorizable. The techniques are being developed to improve the accuracy of calculations, ease-of-use of crashworthiness programs, and the speed of calculations. The latter is still of importance because crashworthiness calculations are often made with models of 20,000 to 50,000 elements using explicit time integration and require on the order of 20 to 100 hours on current supercomputers. The methodologies are briefly reviewed and then some example calculations employing these methods are described. The methods are also of value to other nonlinear transient computations.

Belytschko, Ted↗

NASA NTRS1993

Probabilistic finite elements for fatigue and fracture analysis

An overview of the probabilistic finite element method (PFEM) developed by the authors and their colleagues in recent years is presented. The primary focus is placed on the development of PFEM for both structural mechanics problems and fracture mechanics problems. The perturbation techniques are used as major tools for the analytical derivation. The following topics are covered: (1) representation and discretization of random fields; (2) development of PFEM for the general linear transient problem and nonlinear elasticity using Hu-Washizu variational principle; (3) computational aspects; (4) discussions of the application of PFEM to the reliability analysis of both brittle fracture and fatigue; and (5) a stochastic computational tool based on stochastic boundary element (SBEM). Results are obtained for the reliability index and corresponding probability of failure for: (1) fatigue crack growth; (2) defect geometry; (3) fatigue parameters; and (4) applied loads. These results show that initial defect is a critical parameter.

Belytschko, Ted↗

NASA NTRS1992

A stochastic damage model for the rupture prediction of a multi-phase solid. I - Parametric studies. II - Statistical approach

A stochastic damage model for predicting the rupture of a brittle multiphase material is developed, based on the microcrack-macrocrack interaction. The model, which incorporates uncertainties in locations, orientations, and numbers of microcracks, characterizes damage by microcracking and fracture by macrocracking. A parametric study is carried out to investigate the change of the stress intensity at the macrocrack tip by the configuration of microcracks. The inherent statistical distribution of the fracture toughness arising from the intrinsic random nature of microcracks is explored using a statistical approach. For this purpose, a computer simulation model is introduced, which incorporates a statistical characterization of geometrical parameters of a random microcrack array.

Lua, Yuan J.↗

NASA NTRS1992

Probabilistic finite elements for fatigue and fracture analysis

Attenuation is focused on the development of Probabilistic Finite Element Method (PFEM), which combines the finite element method with statistics and reliability methods, and its application to linear, nonlinear structural mechanics problems and fracture mechanics problems. The computational tool based on the Stochastic Boundary Element Method is also given for the reliability analysis of a curvilinear fatigue crack growth. The existing PFEM's have been applied to solve for two types of problems: (1) determination of the response uncertainty in terms of the means, variance and correlation coefficients; and (2) determination the probability of failure associated with prescribed limit states.

Belytschko, Ted↗

NASA NTRS1991

Fatigue crack growth reliability by probabilistic finite elements

Fusion of the probabilistic finite-element method and reliability analysis for probabilistic fatigue-crack growth is presented. A comprehensive method for determining the probability of fatigue failure for mixed-mode cyclic loading is also presented. The loading is mixed-mode with randomness in the initial and final crack lengths, initial crack angle and position, material properties, crack-growth law, crack-direction law, and loading. The methodology consists of calculating the reliability index via an optimization procedure which is used to calculate the probability of fatigue failure. Performance of the methodology presented is demonstrated on a classical mode-I fatigue problem.

Besterfield, Glen H.↗

NASA NTRS1990

Explicit time integration of finite element models on a vectorized, concurrent computer with shared memory

The implementation of a nonlinear explicit program on a vectorized, concurrent computer with shared memory is described and studied. The conflict between vectorization and concurrency is described and some guidelines are given for optimal block sizes. Several example problems are summarized to illustrate the types of speed-ups which can be achieved by reprogramming as compared to compiler optimization.

Gilbertsen, Noreen D.↗

NASA NTRS1990

Parallel processors and nonlinear structural dynamics algorithms and software

Techniques are discussed for the implementation and improvement of vectorization and concurrency in nonlinear explicit structural finite element codes. In explicit integration methods, the computation of the element internal force vector consumes the bulk of the computer time. The program can be efficiently vectorized by subdividing the elements into blocks and executing all computations in vector mode. The structuring of elements into blocks also provides a convenient way to implement concurrency by creating tasks which can be assigned to available processors for evaluation. The techniques were implemented in a 3-D nonlinear program with one-point quadrature shell elements. Concurrency and vectorization were first implemented in a single time step version of the program. Techniques were developed to minimize processor idle time and to select the optimal vector length. A comparison of run times between the program executed in scalar, serial mode and the fully vectorized code executed concurrently using eight processors shows speed-ups of over 25. Conjugate gradient methods for solving nonlinear algebraic equations are also readily adapted to a parallel environment. A new technique for improving convergence properties of conjugate gradients in nonlinear problems is developed in conjunction with other techniques such as diagonal scaling. A significant reduction in the number of iterations required for convergence is shown for a statically loaded rigid bar suspended by three equally spaced springs.

Belytschko, Ted↗

NASA NTRS1990

Measurement and exploitation of mesh structure on the Connection Machine

A four-node shell element has been implemented on the Connection Machine using the direct vector method and an exchange algorithm previously used with a simpler element. It is shown that the method retains a large measure of its speed in the more complicated setting of a shell element. A hybrid interprocessor communication scheme for the exchange algorithm which uses NEWS within structured subdomains is described.

Plaskacz, Edward J.↗

NASA NTRS1989

Parallel processors and nonlinear structural dynamics algorithms and software

The adaptation of a finite element program with explicit time integration to a massively parallel SIMD (single instruction multiple data) computer, the CONNECTION Machine is described. The adaptation required the development of a new algorithm, called the exchange algorithm, in which all nodal variables are allocated to the element with an exchange of nodal forces at each time step. The architectural and C* programming language features of the CONNECTION Machine are also summarized. Various alternate data structures and associated algorithms for nonlinear finite element analysis are discussed and compared. Results are presented which demonstrate that the CONNECTION Machine is capable of outperforming the CRAY XMP/14.

Belytschko, Ted↗

NASA NTRS1989

Parallel processors and nonlinear structural dynamics algorithms and software

A nonlinear structural dynamics finite element program was developed to run on a shared memory multiprocessor with pipeline processors. The program, WHAMS, was used as a framework for this work. The program employs explicit time integration and has the capability to handle both the nonlinear material behavior and large displacement response of 3-D structures. The elasto-plastic material model uses an isotropic strain hardening law which is input as a piecewise linear function. Geometric nonlinearities are handled by a corotational formulation in which a coordinate system is embedded at the integration point of each element. Currently, the program has an element library consisting of a beam element based on Euler-Bernoulli theory and trianglar and quadrilateral plate element based on Mindlin theory.

Belytschko, Ted↗

NASA NTRS1989

Probalistic Finite Elements (PFEM) structural dynamics and fracture mechanics

The purpose of this work is to develop computationally efficient methodologies for assessing the effects of randomness in loads, material properties, and other aspects of a problem by a finite element analysis. The resulting group of methods is called probabilistic finite elements (PFEM). The overall objective of this work is to develop methodologies whereby the lifetime of a component can be predicted, accounting for the variability in the material and geometry of the component, the loads, and other aspects of the environment; and the range of response expected in a particular scenario can be presented to the analyst in addition to the response itself. Emphasis has been placed on methods which are not statistical in character; that is, they do not involve Monte Carlo simulations. The reason for this choice of direction is that Monte Carlo simulations of complex nonlinear response require a tremendous amount of computation. The focus of efforts so far has been on nonlinear structural dynamics. However, in the continuation of this project, emphasis will be shifted to probabilistic fracture mechanics so that the effect of randomness in crack geometry and material properties can be studied interactively with the effect of random load and environment.

Liu, Wing-Kam↗

NASA NTRS1989

Algorithms and software for nonlinear structural dynamics

The objective of this research is to develop efficient methods for explicit time integration in nonlinear structural dynamics for computers which utilize both concurrency and vectorization. As a framework for these studies, the program WHAMS, which is described in Explicit Algorithms for the Nonlinear Dynamics of Shells (T. Belytschko, J. I. Lin, and C.-S. Tsay, Computer Methods in Applied Mechanics and Engineering, Vol. 42, 1984, pp 225 to 251), is used. There are two factors which make the development of efficient concurrent explicit time integration programs a challenge in a structural dynamics program: (1) the need for a variety of element types, which complicates the scheduling-allocation problem; and (2) the need for different time steps in different parts of the mesh, which is here called mixed delta t integration, so that a few stiff elements do not reduce the time steps throughout the mesh.

Belytschko, Ted↗

NASA NTRS1989

Explicit-explicit subcycling with non-integer time step ratios for structural dynamic systems

An explicit-explicit subcycling procedure for the finite element analysis of structural dynamics is developed. This procedure has relaxed the usual constraint of requiring integer time step ratios for adjacent nodal groups. This allows for greater advantage to be taken of local stability criteria, and thus improves the efficiency of the explicit time integrator. Example problems are included to demonstrate the accuracy and stability of the method.

Neal, Mark O.↗

NASA NTRS1989

An expert system for setting time steps in dynamic finite element programs

An expert system, ETUDES - Expert Time integration control Using Deep and Surface Knowledge System, which addresses the determination of the timestep for time integration of linear structural dynamic equations is described. This timestep may also be applicable for a moderately nonlinear simulation of the same structure. The program also determines whether an explicit or implicit method is most efficient for the particular simulation. A production rule programming system written in OPS5 is used for the implementation of this prototype expert system. Issues relating to the expert system architecture for this application, such as knowledge representation and structure, as well as domain knowledge are discussed. The prototype is evaluated by measuring its performance in various benchmark model problems.

Ramirez, Martin R.↗

NASA NTRS1989

Analytical and computational models of shells; Proceedings of the Symposium, ASME Winter Annual Meeting, San Francisco, CA, Dec. 10-15, 1989

Topics presented include asymptotic analysis and computation for shells, the edge effects in the Reissner-Mindlin plate theory, nonlinear problems of geometrically exact shell theories, and developments in variational methods for high performance plate and shell elements. Also presented are an assumed strain solid element model for geometrically nonlinear shell analysis, shell finite elements with six degrees of freedom per node, hierarchic plate and shell models based on p-extension, and a simple shell element formulation for large-scale elastoplastic analysis. Also discussed are the assessment of computational models for multilayered composite cylinders, shell models for impact analysis, analysis of shell structures subjected to contact-impacts, and the application of shell theory to cardiac mechanics.

Noor, Ahmed K.↗

NASA NTRS1988

Use of adjoint methods in the probabilistic finite element approach to fracture mechanics

The adjoint method approach to probabilistic finite element methods (PFEM) is presented. When the number of objective functions is small compared to the number of random variables, the adjoint method is far superior to the direct method in evaluating the objective function derivatives with respect to the random variables. The PFEM is extended to probabilistic fracture mechanics (PFM) using an element which has the near crack-tip singular strain field embedded. Since only two objective functions (i.e., mode I and II stress intensity factors) are needed for PFM, the adjoint method is well suited.

Liu, Wing Kam↗

NASA NTRS1987

Finite element methods in probabilistic mechanics

Probabilistic methods, synthesizing the power of finite element methods with second-order perturbation techniques, are formulated for linear and nonlinear problems. Random material, geometric properties and loads can be incorporated in these methods, in terms of their fundamental statistics. By construction, these methods are applicable when the scale of randomness is not too large and when the probabilistic density functions have decaying tails. By incorporating certain computational techniques, these methods are shown to be capable of handling large systems with many sources of uncertainties. Applications showing the effects of combined random fields and cyclic loading/stress reversal are studied and compared with Monte Carlo simulation results.

Liu, Wing Kam↗