The derivation of a thick and thin plate formulation without ad hoc assumptions

For the plate formulation considered in this paper, appropriate three-dimensional elasticity solution representations for isotropic materials are constructed. No a priori assumptions for stress or displacement distributions over the thickness of the plate are made. The strategy used in the derivation is to separate functions of the thickness variable z from functions of the coordinates x and y lying in the midplane of the plate. Real and complex 3-dimensional elasticity solution representations are used to obtain three types of functions of the coordinates x, y and the corresponding differential equations. The separation of the functions of the thickness coordinate can be done by separately considering homogeneous and nonhomogeneous boundary conditions on the upper and lower faces of the plate. One type of the plate solutions derived involves polynomials of the thickness coordinate z. The other two solution forms contain trigonometric and hyperbolic functions of z, respectively. Both bending and stretching (or in-plane) solutions are included in the derivation.     

Comparison of stability and control characteristics of two twin-lift helicopter configurations

The stability and control characteristics of two twin-lift helicopter configurations are analyzed in this paper. In order to address the issue of configuration selection from a handling qualities viewpoint, their open-and closed-loop characteristics are compared. The two twin-lift configurations considered are the twin-lift with spreader bar and twin-lift without spreader bar. The nonlinear models describing the dynamics of these two configurations in the lateral/vertical plane are derived. The open-loop characteristics of the two systems are compared by linearizing the nonlinear models about a symmetric hovering equilibrium condition. The closed-loop characteristics of the two systems are compared using nonlinear controllers based on feedback linearization schemes. The performance of the resulting closed-loop systems in carrying out a typical twin-lift mission is evaluated through nonlinear simulation. Also, the effects of helicopter performance degradation and measurement errors on the overall system performance are discussed.[B] Matrix multiplying the control vector in the nonlinear model[B₁] Matrix multiplying the control vector in the linear model[C] Matrix defining vector of variables to be controlled[C₁] Damping matrixC_ijElement of the damping matrix e Parameter used in the linear model = M ₁ h ₁/I ₁=M ₂ h ₂/I ₂,/ft{f} Vector independent of controls in the nonlinear model g Acceleration due to gravity, ft/sec² h₁, h₂Distance of tether attachment point to the center of gravity for helicopters 1 and 2, ft h Parameter used in the linear model, =h ₁=h ₂, ft h Distance between rotor hub and the helicopter center of gravity, ft h h/l H Distance of the load from the spreader bar c.g., ftH₁, H₂Length of tethers 1 and 2, ftI_RMass moment of inertia of spreader bar, slug-ft² I₁, I₂Roll moments of inertia of helicopters 1 and 2, slug-ft² k Non-dimensional hub control moment coefficientK_DDerivative gainsK_IIntegral gainsK_PProportional gains[K_i] Stiffness matrixK_ijElement of the stiffness matrix l Parameter used in the linear model, =H ₁=H ₂, ft L Spreader bar length, ftNomenclature     

Secondary cells and separation in developing laminar curved-pipe flows

Laminar flows through 180° curved bends of circular cross section are investigated numerically. For small curvature ratio, , defined as pipe radius over mean bend radius, the governing equations could be parabolized. The equations are solved for an range of from 0.04 to 0.143, a Dean number (De) range of from 277.5 to 1360, and for a uniform flow, a potential vortex, and a parabolic flow inlet condition. In all these studies a zero cross-stream flow at the inlet is assumed. A detailed study of the effects of , De, and inlet condition on the secondary flow pattern is carried out. Within the range of parameters investigated, up to three secondary cells are found in the cross-stream half-plane of a curved pipe. They are the Dean-type secondary cell, a secondary separation cell near the inner bend (closest to the center of curvature of the bend), and a third cell near the pipe center. The number of secondary cells in the cross-stream half-plane is greatly influenced by the inlet flow, and to a much lesser extent by and De. For example, only the Dean cell is found in a curved-pipe flow where and De are small and the inlet flow is either uniform or a potential vortex. When the inlet condition of the same case is changed to a parabolic flow, a three-cell structure results. Furthermore, as De increases to 1180, incipient axial flow separation begins at around 23° downstream of the curved-pipe entrance. The formation and extent of the separation and third cells are investigated together with their dependence on the parameters studied. This investigation further shows that, within the range of parameters examined, there is no secondary cell occurring near the outer bend, contrary to some earlier findings on fully developed curved-pipe flows.This work was supported by the Office of Naval Research under Grant No. N0014-81-K-0428 and by DTRC, Annapolis, Maryland, under Contract No. N00167-86-K-0075. Also, support in the form of an IPA awarded to RMCS during his sabbatical leave at DTRC, Annapolis, Maryland, in the spring of 1990 is gratefully acknowledged.     

Reduced‐order modelling of an adaptive mesh ocean model

A novel proper orthogonal decomposition (POD) model has been developed for use with an advanced unstructured mesh finite‐element ocean model, the Imperial College Ocean Model (ICOM, described in detail below), which includes many recent developments in ocean modelling and numerical analysis. The advantages of the POD model developed here over existing POD approaches are the ability:

• 1. To increase accuracy when representing geostrophic balance (the balance between the Coriolis terms and the pressure gradient). This is achieved through the use of two sets of geostrophic basis functions where each one is calculated by basis functions for velocities u and v.
• 2. To speed up the POD simulation. To achieve this a new numerical technique is introduced, whereby a time‐dependent matrix in the discretized equation is rapidly constructed from a series of time‐independent matrices. This development imparts considerable efficiency gains over the often‐used alternative of calculating each finite element over the computational domain at each time level.
• 3. To use dynamically adaptive meshes in the above POD model.

Copyright © 2008 John Wiley & Sons, Ltd.     

Nonlinear Strong Shock Interactions: A Shock-Fitted Approach

This paper addresses nonlinear effects which result from the interaction of shock waves with vortices. A series of experiments are carried out, which involve the interaction of a strong shock wave with a single plane vorticity wave and a randomly distributed wave system. These experiments are first conducted in the linear regime to obtain a mutual verification of theory and computation. They are subsequently extended into the nonlinear regime. A systematic study of the interaction of a plane shock wave and a single vortex is then conducted. Specifically, we investigate the conditions under which nonlinear effects become important, both as a function of shock Mach number, M ₁, and incident vortex strength (characterized by its circulation Γ). The shock Mach number is varied from 2 to 8, while the circulation of the vortex is varied from infinitesimally small values (linear theory) to unity. Budgets of vorticity, dilatation, and pressure are obtained. They indicate that nonlinear effects become more significant as both the shock Mach number and the circulation increase. For Mach numbers equal to 5 and above, the dilatation in the vortex core grows quadratically with circulation. An acoustic wave propagates radially outward from the vortex center. As circulation increases, its upstream-facing front steepens at low Mach numbers, and its downstream-facing front steepens at high Mach numbers. A high Mach number asymptotic expansion of the Rankine--Hugoniot conditions reveals that nonlinear effects dominate both the shock motion and the downstream flow for ΓM ₁ > 1. Received 28 June 1997 and accepted 25 November 1997     

Multiple Internal Resonance in Suspended Cables under Random In-Plane Loading

The random excitation of a suspended cable with simultaneous internal resonances is considered. The internal resonances can take place among the first in-plane and the first two out-of-plane modes. The external loading is represented by a wide-band random process. The response statistics are estimated using the Fokker-Planck-Kolmogorov (FPK) equation, together with Gaussian and non-Gaussian closures. Monte Carlo simulation is also used for numerical verification. The unimodal in-plane motion exists in regions away from the internal resonance condition. The mixed mode interaction is manifested within a limited range of internal detuning parameters, depending on the excitation power spectrum density and damping ratios. The Gaussian closure scheme failed to predict bounded solutions of mixed mode interaction. The non-Gaussian closure results are in good agreement with the Monte Carlo simulation. The on-off intermittency of the autoparametrically excited modes is observed in the Monte Carlo simulation over a small range of excitation levels. The influence of the cable parameters, such as damping ratios, sag-to-span ratio, internal detuning parameters, and excitation level on the autoparametric interaction, is studied. It is found that the internal detuning and excitation level are the two main parameters which affect the autoparametric interaction among the three modes. Due to the system's nonlinearity, the response of the three modes is strongly non-Gaussian and the coupled modes experience irregular modulation.     

Nonlinear Nonplanar Dynamics of Parametrically Excited Cantilever Beams

The nonlinear nonplanar response of cantilever inextensional metallic beams to a principal parametric excitation of two of its flexural modes, one in each plane, is investigated. The lowest torsional frequencies of the beams considered are much larger than the frequencies of the excited modes so that the torsional inertia can be neglected. Using this condition as well as the inextensionality condition, we develop a Lagrangian whose variation leads to two integro-partial-differential equations governing the motions of the beams. The method of time-averaged Lagrangian is used to derive four first-order nonlinear ordinary-differential equations governing the modulation of the amplitudes and phases of the two interacting modes. These modulation equations exhibit symmetry properties. A pseudo arclength scheme is used to trace the branches of the equilibrium solutions and an investigation of the eigenvalues of the Jacobian matrix is used to assess their stability. The equilibrium solutions experience pitchfork, saddle-node, Hopf, and codimension-2 bifurcations. A detailed bifurcation analysis of the dynamic solutions of the modulation equations is presented. Five branches of dynamic (periodic and chaotic) solutions were found. Two of these branches emerge from two Hopf bifurcations and the other three are isolated. The limit cycles undergo symmetry-breaking, cyclic-fold, and period-doubling bifurcations, whereas the chaotic attractors undergo attractor-merging and boundary crises.     

The Influence of Moduli Slope of a Linearly Graded Matrix on the Bulk Moduli of Some Particle- and Fiber-Reinforced Composites

The stored energy functional of a homogeneous isotropic elastic body is invariant with respect to translation and rotation of a reference configuration. One can use Noether's Theorem to derive the conservation laws corresponding to these invariant transformations. These conservation laws provide an alternative way of formulating the system of equations governing equilibrium of a homogeneous isotropic body. The resulting system is mathematically identical to the system of equilibrium equations and constitutive relations, generally, of another material. This implies that each solution of the system of equilibrium equations gives rise to another solution, which describes the reciprocal deformation and solves the system of equilibrium equations of another material. In this paper we derive conservation laws and prove the theorem on conjugate solutions for two models of elastic homogeneous isotropic bodies – the model of a simple material and the model of a material with couple stress (Cosserat continuum). This revised version was published online in August 2006 with corrections to the Cover Date.     

Flow mechanisms in dense silica-metal oxide-water suspensions

The rheological properties of dense silica in water suspensions (approx. 50% solids by volume) containing additions of metal oxides were examined. Metal oxides used were ferric, zinc and stannic. To prevent settling, testing was performed in a rheometer which was modified to provide for continual stirring of the materials. Relatively small oxide additions had the effect of thickening the mixtures and making them non-Newtonian. Different rate-limiting steps for flow were identified depending on the particular mixture, testing temperature and shear strain rate. Flow could be described using empirical equations which are identical to those often used to describe plastic flow in solid crystalline materials.     

On the equivalent inclusion method and impotent eigenstrains

The equivalency equations and the nature of the solution are investigated when an inhomogeneity under applied stresses is simulated by an inclusion with eigenstrains. The equivalency equations by which the equivalent eigenstrain is obtained becomes singular when the inhomogeneity is void and the applied stress has a form of polynomials of coordinates of degree one. The solutions of the system of the equivalency equations are not uniquely determined. Explicit expressions are given for the impotent eigenstrains which do not generate any stress field throughout a material.This research was supported by the U.S. Army Research Office under Grant No. DAAG 29-77-G-0042 to Northwestern University.On leave of absence of Meiji University, Tokyo.