Engineering topics
Publications and source records attributed to Thurston, Gaylen A..
Numerical algorithms derived and compared with classical analytical methods. In method, expansions replaced with integrals evaluated numerically. Resulting numerical solutions retain linear independence, main advantage of asymptotic solutions.
Explore the source record for details and available documents.
EAC computer program designed for analysis of errors in results of STAGS computer program (COSMIC Program HQN-10967). Requires input data for geometry of plate, properties of material, and set of boundary conditions. These input data come from STAGS code. (The specific link between input and output data from STAGS and input data for EAC is POSTP, postprocessor program in STAGS processors.) EAC computes continuous solution from discrete results of STAGS in order to estimate error of results of STAGS. Written in FORTRAN 77.
The accurate calculation of stresses at boundaries and interfaces where FEM analysis may be unreliable is presently undertaken by an error analysis that derives a continuous approximation to discrete finite-element data, which can be differentiated to compute continuous stresses for component-failure predictions. An evaluation is conducted of this approximation in the context of the nonlinear PDEs. A novel interpolation formula which is a simple modification of the double Fourier sine series is used to reduce truncation errors near the rectangular plate boundary by means of an 'extended grid'. Results are presented from a FEM solution, a conventional double-Fourier series' continuous approximation, and a solution applying interpolation on the extended grid, which yields superior convergence properties near the plate boundaries.
Explore the source record for details and available documents.
This paper is one of a series on error analysis and correction of finite element solutions for plates and shells. The error analysis in the earlier papers used half-range double Fourier sine series for numerical harmonic analysis. The half-range formulas are simple to apply, but they can be inaccurate near the ends of the ranges of the independent variables. The Gibbs' phenomenon exhibited by half-range sine series in one independent variable has a two-dimensional analog; a classic example is the Navier solution in a double half-range sine series for the simply supported plate under a uniform load. A simple change of variables is introduced in the paper to improve the accuracy of the double Fourier sine series without adding complexity to the numerical analysis. The change of variables is applied to the problem of approximating a transverse load that is tabulated on a rectangular grid. A solution based on the change of variables is compared with results from the Navier solution for the simply supported plate problem and finite element results for the same problem.
Plates can have more than one buckled solution for a fixed set of boundary conditions. The theory for the identification and the computation of multiple solutions in buckled plates is examined. The theory predicts modal interaction (which is also called change in buckle pattern or secondary buckling) in experiments on certain plates with multiple theoretical solutions. A set of coordinate functions is defined for Galerkin's method so that the von Karman plate equations are reduced to a coupled set of cubic equations in generalized coordinates that are uncoupled in the linear terms. An iterative procedure for solving modal interaction problems is suggested based on this cubic form.
The postbuckling response of circular rings (or long cylinders) is examined. The rings are subjected to four types of external pressure loadings; each type of pressure is defined by its magnitude and direction at points on the buckled ring. Newton's method is applied to the nonlinear differential equations of the exact inextensional theory for the ring problem. A zeroth approximation for the solution of the nonlinear equations, based on the mode shape corresponding to the first buckling pressure, is derived in closed form for each of the four types of pressure. The zeroth approximation is used to start the iteration cycle in Newton's method to compute numerical solutions of the nonlinear equations. The zeroth approximations for the postbuckling pressure-deflection curves are compared with the converged solutions from Newton's method and with similar results reported in the literature.
Classical asymptotic analysis of ordinary differential equations derives approximate solutions that are numerically stable. However, the analysis also leads to tedious expansions in powers of the relevant parameter for a particular problem. The expansions are replaced with integrals that can be evaluated by numerical integration. The resulting numerical solutions retain the linear independence that is the main advantage of asymptotic solutions. Examples, including the Falkner-Skan equation from laminar boundary layer theory, illustrate the method of asymptotic analysis with numerical integration.
A computer code is now available for estimating the error in results from the STAGS finite element code for a shell unit consisting of a rectangular orthotropic plate. This memorandum contains basic information about the computer code EAC (Error Analysis and Correction) and describes the connection between the input data for the STAGS shell units and the input data necessary to run the error analysis code. The STAGS code returns a set of nodal displacements and a discrete set of stress resultants; the EAC code returns a continuous solution for displacements and stress resultants. The continuous solution is defined by a set of generalized coordinates computed in EAC. The theory and the assumptions that determine the continuous solution are also outlined in this memorandum. An example of application of the code is presented and instructions on its usage on the Cyber and the VAX machines have been provided.
A general method of error analysis and correction is investigated for the discrete finite-element results for cylindrical shell structures. The method for error analysis is an adaptation of the method of successive approximation. When applied to the equilibrium equations of shell theory, successive approximations derive an approximate continuous solution from the discrete finite-element results. The advantage of this continuous solution is that it contains continuous partial derivatives of an order higher than the basis functions of the finite-element solution. Preliminary numerical results are presented in this paper for the error analysis of finite-element results for a postbuckled stiffened cylindrical panel modeled by a general purpose shell code. Numerical results from the method have previously been reported for postbuckled stiffened plates. A procedure for correcting the continuous approximate solution by Newton's method is outlined.
Explore the source record for details and available documents.
An error analysis of results from finite-element solutions of problems in shell structures is further developed, incorporating the results of an additional numerical analysis by which oscillatory behavior is eliminated. The theory is extended to plates with initial geometric imperfections, and this novel analysis is programmed as a postprocessor for a general-purpose finite-element code. Numerical results are given for the case of a stiffened panel in compression and a plate loaded in shear by a 'picture-frame' test fixture.
A completely parallel algorithm for the symmetric eigenproblem AX = Lambda BX is outlined. The algorithm is parallel in the sense that the numerical operations do not occur in a fixed sequence. Therefore, a large number of operations can be programmed to be performed concurrently on a computer with multiple central processing units. The standard symmetric eigenvalue problem AX = Lambda X has the property that the n eigenvalues of the principal submatrix of A of order n are separated by the (n-1) eignvalues of the principal submatrix of order (n-1). The separation property delineated n intervals containing one eigenvalue. Each eigenvalue and corresponding eigenvector can be computed independently. The n eigenproblem calculations can be divided among multiple processing units.
A direct algorithm is presented for computing the Euler load of a column from experimental data. The method is based on exact inextensional theory for imperfect columns, which predicts two distinct deflected shapes at loads near the Euler load. The bending stiffness of the column appears in the expression for the Euler load along with the column length, therefore the experimental data allows a direct computation of bending stiffness. Experiments on graphite-epoxy columns of rectangular cross-section are reported in the paper. The bending stiffness of each composite column computed from experiment is compared with predictions from laminated plate theory.