A NASTRAN implementation of the doubly asymptotic approximation for underwater shock response

A detailed description is given of how the decoupling approximation known as the doubly asymptotic approximation is implemented with NASTRAN to solve shock problems for submerged structures. The general approach involves locating the nonsymmetric terms (which couple structural and fluid variables) on the right hand side of the equations. This approach results in coefficient matrices of acceptable bandwidth but degrades numerical stability, requiring a smaller time step size than would otherwise be used. It is also shown how the structure's added (virtual) mass matrix, is calculated with NASTRAN

Recent improvements to BANDIT

The NASTRAN preprocessor BANDIT, which improves NASTRAN's computer efficiency by resequencing grid point labels for reduced matrix bandwidth is described. The addition of (1) the Gibbs-Poole-Stockmeyer (GPS) algorithm, and (2) the user option to reduce matrix profile rather than matrix bandwidth is also described. It is shown that, compared to the Cuthill-McKee algorithm on which BANDIT was originally based, GPS is faster and achieves similar results. Current capabilities and options of BANDIT are summarized

Reduction of matrix wavefront for NASTRAN

The three grid point resequencing algorithms most often run by NASTRAN users are compared for their ability to reduce matrix root-mean-square (rms) wavefront, which is the most critical parameter in determining matrix decomposition time in NASTRAN. The three algorithms are Cuthill-McKee (CM), Gibbs-Poole-Stockmeyer (GPS), and Levy. The first two (CM and GPS) are in the BANDIT program, and the Levy algorithm is in WAVEFRONT. Results are presented for a diversified collection of 30 test problems ranging in size from 59 to 2680 nodes. It is concluded that GPS is exceptionally fast and, for the conditions under which the test was made, the algorithm best able to reduce rms wavefront consistently well

Comparison of finite element analysis of a piping tee using NASTRAN and CORTES/SA

A comparison of finite element analyses of a piping tee was made using NASTRAN and CORTES/SA, a modified version of SAP3 having a special purpose input processor for generating geometries for a wide variety of tee joints. Four finite element models were subjected in force, moment, and pressure loadings. Flexibility factors and principal stresses were computed for each model and compared with results obtained experimentally by Combustion Engineering, Inc. Results from the NASTRAN analyses were in good agreement with experimental results for all loadings except internal pressure. The CORTES/SA analyses gave good results for the internal pressure loading, but poorer results for out of plane bending moments or forces resulting in out of plane bending. Two of the basic load cases in CORTES/SA were found to contain errors that could not be easily corrected. COST COMPARison of NASTRAN and CORTES/SA showed NASTRAN to be less expensive to two than CORTES/SA for identical meshes

The dynamic analysis of submerged structures

Methods are described by which the dynamic interaction of structures with surrounding fluids can be computed by using finite element techniques. In all cases, the fluid is assumed to behave as an acoustic medium and is initially stationary. Such problems are solved either by explicitly modeling the fluid (using pressure or displacement as the basic fluid unknown) or by using decoupling approximations which take account of the fluid effects without actually modeling the fluid

Finite element analysis of fluid-filled elastic piping systems

Two finite element procedures are described for predicting the dynamic response of general 3-D fluid-filled elastic piping systems. The first approach, a low frequency procedure, models each straight pipe or elbow as a sequence of beams. The contained fluid is modeled as a separate coincident sequence axial members (rods) which are tied to the pipe in the lateral direction. The model includes the pipe hoop strain correction to the fluid sound speed and the flexibility factor correction to the elbow flexibility. The second modeling approach, an intermediate frequency procedure, follows generally the original Zienkiewicz-Newton scheme for coupled fluid-structure problems except that the velocity potential is used as the fundamental fluid unknown to symmetrize the coefficient matrices. From comparisons of the beam model predictions to both experimental data and the 3-D model, the beam model is validated for frequencies up to about two-thirds of the lowest fluid-filled labor pipe mode. Accurate elbow flexibility factors are seen to be crucial for effective beam modeling of piping systems

Finite Element Prediction of Loss Factors for Structures with Frequency-dependent Damping Treatments

A finite element procedure is described for calculating the loss factors for elastic structures to which frequency-dependent viscoelastic damping treatments were applied. The frequency dependence of the viscoelastic damping material is treated by approximating its shear modulus with a second-order polynomial so that the stiffnesses associated with the constant, linear, and quadratic terms can be combined, respectively, with the stiffness, damping, and mass matrices assembled for the rest of the structure. A single complex eigenvalue analysis is then performed in which the eigenvalues are purely imaginary. The loss factor is computed by the modal strain energy (MSE) approach. In the the MSE approach, the loss factor of a composite structure vibrating in one of its natural modes may be visualized as a weighted average of the loss factors of the component parts, with the relative stored energies as weighting constants. The finite element procedure, which can treat very general geometries, is illustrated for the case of a vibrating constrained-layer damped plate

A general low frequency acoustic radiation capability for NASTRAN

A new capability called NASHUA is described for calculating the radiated acoustic sound pressure field exterior to a harmonically-excited arbitrary submerged 3-D elastic structure. The surface fluid pressures and velocities are first calculated by coupling a NASTRAN finite element model of the structure with a discretized form of the Helmholtz surface integral equation for the exterior fluid. After the fluid impedance is calculated, most of the required matrix operations are performed using the general matrix manipulation package (DMAP) available in NASTRAN. Far field radiated pressures are then calculated from the surface solution using the Helmholtz exterior integral equation. Other output quantities include the maximum sound pressure levels in each of the three coordinate planes, the rms and average surface pressures and normal velocities, the total radiated power and the radiation efficiency. The overall approach is illustrated and validated using known analytic solutions for submerged spherical shells subjected to both uniform and nonuniform applied loads