Implementation of elastic-plastic structural analysis into NASTRAN

Elastic-plastic analytic capabilities were incorporated into the NASTRAN program. The present implementation includes a general rigid format and additional bulk data cards as well as to two new modules. The modules are specialized to include only perfect plasticity of the CTRMEN and CROD elements but can easily be expanded to include other plasticity theories and elements. The practical problem of an elastic-plastic analysis of a ship's bracket connection is demonstrated and compared to an equivalent analysis using Grumman's PLANS program. The present work demonstrates the feasibility of incorporating general elastic-plastic capabilities into NASTRAN

DYCAST: A finite element program for the crash analysis of structures

DYCAST is a nonlinear structural dynamic finite element computer code developed for crash simulation. The element library contains stringers, beams, membrane skin triangles, plate bending triangles and spring elements. Changing stiffnesses in the structure are accounted for by plasticity and very large deflections. Material nonlinearities are accommodated by one of three options: elastic-perfectly plastic, elastic-linear hardening plastic, or elastic-nonlinear hardening plastic of the Ramberg-Osgood type. Geometric nonlinearities are handled in an updated Lagrangian formulation by reforming the structure into its deformed shape after small time increments while accumulating deformations, strains, and forces. The nonlinearities due to combined loadings are maintained, and stiffness variation due to structural failures are computed. Numerical time integrators available are fixed-step central difference, modified Adams, Newmark-beta, and Wilson-theta. The last three have a variable time step capability, which is controlled internally by a solution convergence error measure. Other features include: multiple time-load history tables to subject the structure to time dependent loading; gravity loading; initial pitch, roll, yaw, and translation of the structural model with respect to the global system; a bandwidth optimizer as a pre-processor; and deformed plots and graphics as post-processors

Groundbreaking Ceremony Address

Address by Governor Ogilvie from the groundbreaking ceremony for Governors State University on June 12, 1971

Magnetorotational-type instability in Couette-Taylor flow of a viscoelastic polymer liquid

We describe an instability of viscoelastic Couette-Taylor flow that is directly analogous to the magnetorotational instability (MRI) in astrophysical magnetohydrodynamics, with polymer molecules playing the role of magnetic field lines. By determining the conditions required for the onset of instability and the properties of the preferred modes, we distinguish it from the centrifugal and elastic instabilities studied previously. Experimental demonstration and investigation should be much easier for the viscoelastic instability than for the MRI in a liquid metal. The analogy holds with the case of a predominantly toroidal magnetic field such as is expected in an accretion disk and it may be possible to access a turbulent regime in which many modes are unstable.Comment: 4 pages, 4 figures, to be published in Physical Review Letter

Effective sigma models and lattice Ward identities

We perform a lattice analysis of the Faddeev-Niemi effective action conjectured to describe the low-energy sector of SU(2) Yang-Mills theory. To this end we generate an ensemble of unit vector fields ("color spins") n from the Wilson action. The ensemble does not show long-range order but exhibits a mass gap of the order of 1 GeV. From the distribution of color spins we reconstruct approximate effective actions by means of exact lattice Schwinger-Dyson and Ward identities ("inverse Monte Carlo"). We show that the generated ensemble cannot be recovered from a Faddeev-Niemi action, modified in a minimal way by adding an explicit symmetry-breaking term to avoid the appearance of Goldstone modes.Comment: 25 pages, 17 figures, JHEP styl

Global axisymmetric Magnetorotational Instability with density gradients

We examine global incompressible axisymmetric perturbations of a differentially rotating MHD plasma with radial density gradients. It is shown that the standard magnetorotational instability, (MRI) criterion drawn from the local dispersion relation is often misleading. If the equilibrium magnetic field is either purely axial or purely toroidal, the problem reduces to finding the global radial eigenvalues of an effective potential. The standard Keplerian profile including the origin is mathematically ill-posed, and thus any solution will depend strongly on the inner boundary. We find a class of unstable modes localized by the form of the rotation and density profiles, with reduced dependence on boundary conditions.Comment: 22 pages, 5 figure

Global m=1 modes and migration of protoplanetary cores in eccentric protoplanetary discs

We calculate global $m=1$ modes with low pattern speed corresponding to introducing a finite eccentricity into a protoplanetary disc. We consider disc models which are either isolated or contain one or two protoplanets orbiting in an inner cavity. Global modes that are strongly coupled to inner protoplanets are found to have disc orbits which tend to have apsidal lines antialigned with respect to those of the inner protoplanets. Other modes corresponding to free disc modes may be global over a large range of length scales and accordingly be long lived. We consider the motion of a protoplanet in the earth mass range embedded in an eccentric disc and determine the equilibrium orbits which maintain fixed apsidal alignment with respect to the disc gas orbits. Equilibrium eccentricities are found to be comparable or possibly exceed the disc eccentricity. We then approximately calculate the tidal interaction with the disc in order to estimate the orbital migration rate. Results are found to deviate from the case of axisymmetric disc with near circular protoplanet orbit once eccentricities of protoplanet and disc orbits become comparable to the disc aspect ratio in magnitude. Aligned protoplanet orbits with very similar eccentricity to that of the gas disc are found to undergo litle eccentricity change while undergoing inward migration in general. However, for significantly larger orbital eccentricities, migration may be significantly reduced or even reverse from inwards to outwards. Thus the existence of global non circular motions in discs with radial excursions comparable to the semi-thickness may have important consequences for the migration and survival of protoplanetary cores in the earth mass range.Comment: Accepted for publication by A &

On the tilting of protostellar disks by resonant tidal effects

We consider the dynamics of a protostellar disk surrounding a star in a circular-orbit binary system. Our aim is to determine whether, if the disk is initially tilted with respect to the plane of the binary orbit, the inclination of the system will increase or decrease with time. The problem is formulated in the binary frame in which the tidal potential of the companion star is static. We consider a steady, flat disk that is aligned with the binary plane and investigate its linear stability with respect to tilting or warping perturbations. The dynamics is controlled by the competing effects of the m=0 and m=2 azimuthal Fourier components of the tidal potential. In the presence of dissipation, the m=0 component causes alignment of the system, while the m=2 component has the opposite tendency. We find that disks that are sufficiently large, in particular those that extend to their tidal truncation radii, are generally stable and will therefore tend to alignment with the binary plane on a time-scale comparable to that found in previous studies. However, the effect of the m=2 component is enhanced in the vicinity of resonances where the outer radius of the disk is such that the natural frequency of a global bending mode of the disk is equal to twice the binary orbital frequency. Under such circumstances, the disk can be unstable to tilting and acquire a warped shape, even in the absence of dissipation. The outer radius corresponding to the primary resonance is always smaller than the tidal truncation radius. For disks smaller than the primary resonance, the m=2 component may be able to cause a very slow growth of inclination through the effect of a near resonance that occurs close to the disk center. We discuss these results in the light of recent observations of protostellar disks in binary systems.Comment: 21 pages, 7 figures, to be published in the Astrophysical Journa

Carbenicillin-Induced Coagulopathy

In high dosage, carbenicillin may interfere with the conversion of fibrinogen to fibrin and result in a haemorrhagic diathesis. The effect is dose-dependent and requires a high concentration of carbenicillin in plasma. Such a level may be attained in renal failure unless the dose of the drug is appropriately reduced. In such situations the screening coagulation tests may be prolonged and this should alert one to the development of a haemorrhagic diathesis due to the drug.S. Afr. Med. J., 48, 457 (1974