ALBERT

Description: In many structures the final stress states are dependent on the sequence of construction or the stress states at various stages of construction are of interest. Such problems can be analyzed using finite element programs that have the capability of adding (birthing) elements to simulate the progress of construction. However, the usual procedure of assembling elements may lead to numerical instabilities or stress states that are unrealistic. Both problems are demonstrated in the analysis of a structure using the program ADINA. A technique which combines application of a preload with element birthing to overcome these problems is described and illustrated.

Description: A technique utilizing finite element analysis, liquid impact kinematics, and momentum theory is described and compared to single-drop impact test data performed on various configurations of coated ceramic material. The method correlates well with test data and is useful in predicting the single-drop impact damage velocity threshold for low-density, coated ceramic materials.

Description: An algorithm is outlined for simulating the contact surface erosion for impact problems. The algorithm dynamically relocates the contact surface as projectile and target materials exceed their failure criterion. Example computations of axisymmetric and oblique impacts are compared with experimental data.

Description: A study of the dynamic characteristics of a coupled translational-rotational system is given. The formulation of the problem considers the soil-structure interaction effects by utilizing the impedance functions at the foundation of a structure. Due to the fact that the coefficient matrix in the characteristic equation is frequency dependent in nature, iterations have to be performed to find the nature frequencies of the system. Examples and discussions are presented. Comparisons of the analytical results from various approaches are also given.

Description: Two analytical techniques applicable to large deflection dynamic response calculations for pressure loaded composite sandwich panels are demonstrated. One technique utilizes finite element modeling with a single equivalent layer representing the face sheets and core. The other technique utilizes the modal analysis computer code DEPROP which was recently modified to include transverse shear deformation in a core layer. The example problem consists of a simply supported rectangular sandwich panel. Included are comparisons of linear and nonlinear static response calculations, in addition to dynamic response calculations.

Description: The design of thin shell structures with respect to elastoplastic buckling requires an extended analysis of the influence of initial imperfections. For conservative design, the most critical defect should be assumed with the maximum allowable magnitude. This defect is closely related to the initial postbuckling behavior. An algorithm is given for the quasi-static analysis of the postbuckling behavior of structures that exhibit multiple buckling points. the algorithm based upon an energy criterion allows the computation of the critical perturbation which will be employed for the definition of the critical defect. For computational efficiency, the algorithm uses the reduced basis technique with automatic update of the modal basis. The method is applied to the axisymmetric buckling of cylindrical shells under axial compression, and conclusions are given for future research.

Description: A general approach is required for describing matter of behavior when the failure is likely to involve growth and coalescence of a large number of fractures. Failures of this kind appear frequently in rapid dynamic processes, particularly in the formation of spall fragments. An approach to formulating constitutive relations that accounts for the opening, shear and growth of an ensemble of cracks is discussed. The approach accounts for plastic flow accompanying fragmentation. The resulting constitutive relations were incorporated into a Lagrangian computer program. A theoretical approach to coalescence is described. The simplest formulation uses a linear Liouville equation, with crack growth limited by the mean free path of cracks, assumed constant. This approach allows for an anisotropic distribution of cracks. An alternative approach in which the decrease of the mean free path with increasing crack size is accounted for, but the crack distribution is assumed isotropic is described. A reduction of the governing Liouville equation to an ordinary differential equation of third order is possible, and the result can be used to determine how mean free path decreases with increasing crack size.

Description: The design of a composite panel requires some way of finding the minimum thickness laminate which will withstand the load requirements without failure. The mathematical complexity of this problem dictates the use of nonlinear optimization techniques. Specialized laminate optimization programs were developed which are compact and efficient enough to run on microcomputers. Only stresses at a point and inplane loads and deflections are considered. The programs are simple to use and require no knowledge of optimization. Techniques are developed which find minimum thickness laminates with either ply ratios or ply angles as design variables. A method is presented for finding the optimum orientation for the axis of symmetry of an orthotropic laminate. The orthotropic laminate program uses an approximate failure theory, which speed up computations dramatically.

Description: The evaluation of the total probability of a plastic collapse failure P sub f for a highly redundant structure of random interdependent plastic moments acted on by random interdepedent loads is a difficult and computationally very costly process. The evaluation of reasonable bounds to this probability requires the use of second moment algebra which involves man statistical parameters. A computer program which selects the best strategy for minimizing the interval between upper and lower bounds of P sub f is now in its final stage of development. The relative importance of various uncertainties involved in the computational process on the resulting bounds of P sub f, sensitivity is analyzed. Response sensitivities for both mode and system reliability of an ideal plastic portal frame are shown.

Description: A finite element modeling technique which utilizes a triangular element called TSHEL with 45 degrees of freedom and seven point integration was tested for analysis of thin plate and shell structures. The element formulation is based on the degenerate solid shell concept and the mixed formulation with assumed independent inplane and transverse shear strains. The effectiveness of the present modeling technique which features combined use of elements with kinematic modes and those without kinematic modes to eliminate both locking and spurious kinematic modes at the global structural levels are shown.

Description: The stress derivative technique for control of keystoning deformation in under-integrated finte elements is based on expansion of the stress in a Taylor series about the element center and retention of additional terms beyond the constant stress term. It has the advantage over other control techniques that keystoning resistance is provided by actual rather than artificial material properties. Application of this technique to the quadrilateral ring elements used for modelling solids of revolution subjected to axisymmetric loads is described. In a cylindrical coordinate system additional terms appear in the formulation which must be dealt with in arriving at a workable keystoning control scheme.

Description: An existing program is currently being adapted to perform finite element analysis by distributing substructures over a network of four Apple IIe microcomputers connected to a shared disk. In this network, one microcomputer controls the entire process while the others perform the analysis on each substructure in parallel. This substructure analysis is used in an iterative, fully stressed, structural resizing procedure. This procedure allows experimentatation with resizing in which all analyses are not completed during a single iteration. This research gives some insight on how to configure multidiscriplinary analysis and optimization procedures for decomposable engineering systems using either high performance engineering workstations or a parallel processor supercomputer. In addition, the operational experience gained facilitates the implementation of analysis programs on these new computers when they become available in an engineering environment.

Description: Described are the experiences gained from solving for the dynamic response of two simple structures on an experimental Multiple Instruction Multiple Data (MIMD) computer called the finite element machine. Introduced are MIMD computing concepts, describing how the concurrent algorithmic techniques implemented and giving results for the two example problems. The results show computational speedups of up to 7.83 using eight of the finite element machine processors and indicate that significant computational speedups are possible for large order structural computations.

Description: The FLEX/32 Multicomputer is a generic environment for cooperating multiple processors. The FLEX/32 supports a number of different processors, making it heterogeneous in terms of the instruction sets it supports, and homogeneous in its ability to provide consistent storage and input/output facilities to its differing processors. These facilities are accessed through standard 32-bit VMEbus connections. The FLEX/32 supports the full UNIX System V Operating System and languages associated with it, plus the extended ConCurrent C and Concurrent FORTRAN 77 languages that allow programming of concurrent software at a high level. Direct programming support at all levels is provided by the environment hardware for concurrent software execution and optimization, including hardware support for shared resource access arbitration, conditional critical region arbitration, and interprocessor messages.

Description: Research aimed at faster, more cost effective parallel machines and algorithms for improving designer productivity with finite element computations is discussed. A set of 8 boards, containing 4 nearest neighbor connected arrays of commercially available floating point chips and substantial memory, are inserted into a commercially available machine. One-tenth Mflop (64 bit operation) processors provide an 89% efficiency when solving the equations arising in a finite element problem for a single variable regular grid of size 40 by 40 by 40. This is approximately 15 to 20 times faster than a much more expensive machine such as a VAX 11/780 used in double precision. The efficiency falls off as faster or more processors are envisaged because communication times become dominant. A novel successive overrelaxation algorithm which uses cyclic reduction in order to permit data transfer and computation to overlap in time is proposed.

Description: The design, manufacture, and testing of an isolation system placed under a bench which simulates a space platform during testing in a hard vacuum are described. A low natural frequency isolation system is used to decouple ground vibrations from the bench. The materials used are vacuum compatible and do not introduce unwanted matter which would contaminate the payload optics and/or the testing environment. The system accommodates payloads of varying weights and envelopes.

Description: The potential cost and performance advantages of welding was understood but ignored by solar panel manufacturers in the U.S. Although NASA, DOD and COMSAT have supported welding development efforts, soldering remains the only U.S. space qualified method for interconnecting solar cells. The reason is that no U.S. satellite prime contractor found it necessary, due to mission requirements, to abandon the space proven soldering process. It appears that the proposed NASA space station program will provide an array requirement, a 10 year operation in a low Earth orbital environment, that mandates welding. The status of welding technology in the U.S. is assessed.

Description: Researchers define damage-tolerant structural systems as those systems which not only have adeqate intact strength to withstand initial failure but also adequate residual strength to minimize the possibility of, and hence the consequences of, further failure. The incorporation of damage tolerance cannot be done in total isolation of the function being required of the system and the costs associated with obtaining improved damage tolerance. The approach, therefore, is to formulate multiple-objective, multi-level decision support problems (DSP), the solutions of which represent a compromise between higher costs and higher damage tolerance. Mulitple-objective decision support problems are easily solved in the linear domain. These formulations, however, include both linear and nonlinear constraints and goals, which in the past, have not been considered due to the resulting complexity. Here, researchers: (1) present a complete discussion and description of decision support problems; (2) identify what further research needs to be done in order to obtain information that is required but not known for solving problems using these models; and (3) identify what needs to be done to implement this prototype method in practice.

Description: Ocean-going vehicles and platforms are among the largest structures in the world and are subjected to relatively harsh conditions of motions and loads. Some of them, such as semi-submersible platforms, are a relatively new type of structure and hence there is no formal, well evolved and established structural design code as there is for more traditional structures. More recently, efforts have also been made to develop a design method of this type for ships and other ocean structures. One of the many advantages of a rationally based design method is versatility; it can be used for structures that have widely differing purposes, measures of merit, shapes and sizes. The purpose is to describe a rationally based design method that has been developed within the field of ocean structures, in order that persons dealing with other types of structure can judge whether and to what extent its various features may be useful for those other types. Also, even though some features may not be applicable they might stimulate some useful ideas.

Description: The use of microcomputers in the design of a cable catenary large space antenna system is discussed. The development of a system design capability, data base utilization, systems integration, program structure and logic, and integrated graphics output are discussed.

Description: An exceptionally elegant method for structural optimization with constraints on the static response presented by Shield and Prager is discussed. Their derivation of the optimality condition was facilitated by a reformulation of the structural elasticity equations in terms of what was then a new variational principle, the principle of stationary mutual potential energy. Their optimality condition relates the design variable to an appropriately defined mutual strain energy. An alternative but related approach, based upon the principle of stationary mutual complementary energy, presented by N. C. Haung, is also discussed. The simplicity of these principles lies in the facts that the energy functionals are stationary at the solution to the field equations and that their stationary value is proportional to the quantity to be optimized.

Description: The effects of space environments on damping materials and damping designs on flexible structures were investigated. The following items were examined: damping of flexible spacecraft appendages; composite loss factor (n sub s) vs. time in high vacuum for damped test beams and damping of flexible structures. The STEP experiments show inherent damping of flexible structures in space effective possible damping design configurations for space structures, effects of passively damped components on the system loss factor of flexible structures and the effect of space environment on properties of damping materials.

Description: Interplanetary dust particles (micrometeoroids) are expected to form well-defined craters upon impacting exposed material in space. Studying the frequency and features of these craters will provide data on the mass-flux distribution of micrometeoroids and, to a lesser extent, on the velocity magnitude and direction. This experiment will study impact craters produced by micrometeoroids on selected materials (metals and glasses in the form of thick targets) to obtain valuable technological and scientific data. Specifically, the studies will focus on determining micrometeoroid composition and mass-flux distribution. Analyses will also be made on the distribution of impact velocity vectors.

Description: In-flight calibration for the solr and Earth flux channels was examined. Earth Radiation on Budget (ERB) channel components were exposed to the space environment and then retrieved and resubmitted to radiometric calibration after exposure. It is suggested that corrections may be applied to ERB results and information will be obtained to aid in the selection of components for future operational solar and Earth radiation budget experiments. To assure that these high accuracy devices are measuring real variations and are not responding to changes induced by the space environment, it is desirable to test such devices radiometrically after exposure to the best approximation of the orbital environment.

Description: The extent and thermal stratification of the region of convective overshoot underneath the convection zone of the sun are investigated. The phenomenon of convective overshoot in general is discussed, and some of the modal and model approaches to studying it are briefly reviewed. A detailed theoretical description of the motion of plumes in a stably stratified medium is given, leading to a 'derivation' of the plume equations from the hydrodynamic equations. Entrainment is discussed, and it is shown how the plume equations can be used to compute convective overshoot in the sun. The limitations of the plume model are addressed, arguing that a thin boundary layer must exist which separates convective and radiative regions. The results of numerical integrations of the plume equations, as applied to the region of convective overshoot underneath the solar convective zone, are discussed.

Description: It is demonstrated that the common assumption made in solar flare beam transport theory that the beam-accompanied return current is purely electrostatically driven is incorrect, and that the return current is both electrostatically and inductively driven, in accordance with Lenz's law, with the inductive effects dominating for times greater than a few plasma periods. In addition, it is shown that a beam can only exist in a solar plasma for a finite time which is much smaller than the inductive return current dissipation time. The importance of accounting for the role of the acceleration mechanism in forming the beam is discussed. In addition, the role of return current driven anomalous resistivity and its subsequent anomalous Joule heating during the flare process is elucidated.

Description: The way in which the initial development of solar filament radiative cooling and the magnetic reconnection of a solar flare can occur in the center of a field-shear layer is demonstrated. Since the present treatment unites these two mechanisms, it indicates the common as well as the disparate features they possess. Unstable radiation serves to increase the Coulomb resistivity at the X-point, so that the reconnection is not self-quenching. The surprising dominance of the magnetic component of the perturbation in the midwavelength range indicates the need to examine the nonlinear saturation of the energy transport of the radiative mode, taking the accompanying magnetic reconnection and potential-energy release into account, for comparison with observations of filaments as well as for clues to the character of the preflare state.

Description: Using NASA's Tracking and Data Relay Satellite as a communications link, astronomers are able to receive scans from the Solar Maximum Mission (SMM) satellite immediately and regularly at the Goddard Space Flight Center. This major operational improvement permits the examination of SMM imagery and spectra as they arrive, as well as the formulation of future observational sequences on the basis of the solar activity in progress. Attention is given to aspects of the sun that change in the course of the 11-year sunspot cycle's movement from maximum to minimum. Proof has been obtained by means of SMM for the near-simultaneity of X-ray and UV bursts at flare onset.

Description: A geometrically nonlinear finite element analysis of cohesive failure in typical joints is presented. Cracked-lap-shear joints were chosen for analysis. Results obtained from linear and nonlinear analysis show that nonlinear effects, due to large rotations, significantly affect the calculated mode 1, crack opening, and mode 2, inplane shear, strain-energy-release rates. The ratio of the mode 1 to mode 2 strain-energy-release rates (G1/G2) was found to be strongly affected by the adhesive modulus and the adherend thickness. The ratios between 0.2 and 0.8 can be obtained by varying adherend thickness and using either a single or double cracked-lap-shear specimen configuration. Debond growth rate data, together with the analysis, indicate that mode 1 strain-energy-release rate governs debond growth. Results from the present analysis agree well with experimentally measured joint opening displacements. Previously announced in STAR as N83-13497

Description: Recent investigations using measurements at 1 AU have discovered three types of long term variation in the interplanetary magnetic field: solar minimum decreases, solar maximum enhancements, and small decreases around solar reversal. In this study the 1972-1982 Helios 1, 2, ISEE-3, and Pioneer 10, 11 observations between 0.3 and 12 AU are examined to further investigate these changes. It was found that all three IMF solar cycle effects are also present in the Helios and Pioneer measurements, confirming that these variations occur throughout the low latitude heliosphere. In addition, the comparison of measurements by identical magnetometers on ISEE-3, Pioneer 10 and Pioneer 11 has revealed a more rapid decrease in IMF intensity than predicted by classical Parker theory. Causes and ramifications of both the long term variations and steeper-than-expected radial gradients in the interplanetary magnetic field are discussed.

Description: A numerical investigation is conducted into the way in which a solar wind model initially satisfying both steady state and energy balance conditions is disturbed and deformed, under the assumption of heating that correspoonds to the energy release of solar flares of an importance value of approximately 1 which occur in radial open field regions. Flare-associated solar wind transient behavior is modeled for 1-8 solar radii. The coronal temperature around the heat source region rises, and a large thermal conductive flux flows inward to the chromosphere and outward to interplanetary space along field lines. The speed of the front of expanding chromospheric material generated by the impingement of the conduction front on the upper chromosphere exceeds the local sound velocity in a few minutes and eventually exceeds 100 million cm/sec.

Description: The numbers and spectra of the accelerated protons and nuclei that produce the neutrons and gamma-rays observed in solar flares are derived, and the results are compared with interplanetary observations of flare protons. The two most widely studied flare acceleration mechanisms, stochastic and diffusive shock acceleration, are discussed, and the arguments favoring the thick-target interaction model for neutron and gamma-ray production at the sun are briefly reviewed. The pertinent results of the theory of neutron and gamma-ray production are presented. The number and spectrum of the accelerated particles are derived from observations of nuclear deexcitation lines and the 2.223 MeV line from several flares. The June 21, 1980 and June 3, 1982 flares, from which a wealth of neutron, gamma-ray and energetic-particle data has recently become available, are discussed.

Description: Time sequences of recurrent mass ejections have been observed during a coordinated SMY program (Sept. 1, 1980 - Sept. 23, 1980 - Oct. 2, 1980). Comparison of the temporal evolution of H-alpha and CIV brightnesses shows a weak phase lag between H-alpha and CIV maxima, in the case of homologous flares, with CIV brightness maxima preceding H-alpha maxima. The analysis of the variation of the ejection velocities is expected to lead to the determination of an energy balance. Such recurrent ejections could be due to periodic energy storage and periodic reorganization of magnetic field as envisaged to occur for flares, but at lower energy levels.

Description: In NOAA Active Region 2372 (April 1980), 4 x 10 to the 20th maxwells of magnetic flux concentrated in an area 30 arcsec across disappeared overnight. Vector magnetograms show that all components of the magnetic field weakened together. If the field had weakened through diffusion or fluid flow, 90 percent of the original flux would still have been detected by the magnetograph within a suitably enlarged area. In fact there was a threefold decrease in detected flux. Evidently, magnetic field was removed from the photosphere. Since the disappearing flux was located in a region of low magnetic shear and low activity in H-alpha and Ly-alpha, it is unlikely that the field dissipated through reconnection. It is argued that the most likely possibility is that flux submerged. The observations suggest that even during the growth phase of active regions, submergence is a strong process comparable in magnitude to emergence.

Description: Coronal bullets are small ejecta of cool, dense plasma observed to accelerate through the solar atmosphere from 20 to 450 km/s. The NRL Dynamic Flux Tube Model has been used to simulate the evolving physical properties of these dynamic events. The present calculations utilize an adaptive-gridding technique to resolve the fine structure within and around the bullets. In this work, an identification was made of a component of shocked plasma which piles up ahead of the bullet and eventually dominates both the dynamics and heating of the original bullet mass. The observational consequences of this shocked component are discussed in terms of the available HRTS EUV data, and suggestions are made for optimizing future observations of this phenomenon. An investigation has also been conducted of the structure of the bullet material visible in EUV spectral lines and the observable characteristics of the EUV-emitting plasma. Finally, the most likely mechanisms for accelerating the bullets, as well as favorable sites of origin are evaluated.

Description: Studies have been conducted to develop efficient techniques to simulate crack extension and to examine various local and global fracture criteria. Of the considered criteria, the crack-tip-opening angle (CTOA) or displacement (CTOD) at a specified distance from the crack tip was shown to be most suited for modeling stable crack growth and instability during the fracture process. The results obtained in a number of studies show the necessity for studying different crack configurations when assessing the validity of any fracture criteria. One of the objectives of the present investigation is related to a critical evaluation of the CTOD growth criterion using an elastic-plastic finite element analysis under monotonic loading to failure. The analysis was found to predict three stages of crack growth behavior under monotonic loading to failure. Calculated CTOD values agreed well with experimental values for crack growth initiation.

Description: Low-noise (S/N greater than 100), high spectral resolution observations of two pure rotation transitions of OH from the solar photosphere are used to make inferences concerning the thermal structure and inhomogeneity of the upper photosphere. It is found that the v = O R22(24.5)e line strengthens at the solar limb, in contradiction to the predictions of current one-dimensional photospheric models. The results for this line support a two-dimensional model in which horizontal thermal fluctuations in the upper photosphere are of the order plus or minus 800 K. This thermal bifurcation may be maintained by the presence of magnetic flux tubes and may be related to the solar limb extensions observed in the 30-200-micron region.

Description: Solar irradiance measurements from the ACRIM experiment show a clear response to the rotation periods of g-mode oscillations (l = 1, 2, and 3) and their first harmonics. Peaks in the ACRIM spectrum at 16.6, 18.3, 20.7, 36.5, and about 71 days all lie within about 1 percent of periods arising from g-mode rotation. This means that the g-modes are a fundamental cause of irradiance fluctuations. On time scales of months and less they modulate the irradiance by means of transient flows of global scale which they stimulate in the sun's convective envelope. Dimensional arguments indicate that the flows carry up heat at an average rate of about 0.001 solar luminosities, which is not in conflict with observed changes in the irradiance. Five additional tests for g-modes and large-scale convection are given. An instability is described which undermines diffusion models of sunspot energy storage.

Description: The result of a study on the application of an improved statistical prediction method for estimating the intermediate-term (months) and long-term (years) behavior of solar flux is discussed. The study indicates that better predictions, in a chi square sense, are possible by selecting sets of the solar flux data such that each set (cycle) starts and ends at the maxima (or minima) for the data base and initialization point of the procedure. Then one applies a Lagrangian least-squares statistical technique. Evidence is also presented to support the existence of an aperiodic variation in the periods as well as the amplitudes.

Description: Attenuation of the load transmitted to a vehicle occupant by the structure, either by modifying structural assembly, changing geometry of its elements, or adding specific load limiting devices to help dissipate kinetic energy was studied. General aviation aircraft have fuselage subfloors of a built up structure which are generally very stiff perpendicular to the floor of the cabin. The subfloor structure, designed to crush at an appropriate force level, can be used to advantage in a crash by dissipating energy through plastic buckling of the floor beams and frames. Simple closed form solutions to predict the mean crushing force levels of subfloor designs is useful in engineering practice. The crushing process of thin walled, plate formed, open structures with particular emphasis on L and cruciform shapes was analyzed. Lower and upper bound solutions for the mean crushing strength of cruciforms are obtained by considering modes of deformation which account for both bending and extensional deformation. It is the importance of extensional deformation to the energy absorption process, representing at least one third of the dissipated energy is shown. Previously announced in STAR as N82-33765

Description: This work relates to a series of collaborative investigations involving the application of a computational model for the determination of the detailed plasma and magnetic field properties associated with the global interaction of the solar wind with various planetary obstacles throughout the solar system. The theoretical method is based on an established single fluid, steady, dissipationless, magnetohydrodynamic continuum model, and is appropriate for the calculation of supersonic, super-Alfvenic solar wind flow past planetary obstacles. The investigations undertaken relate to studies of various solar wind interaction phenomena with Venus, Earth, Jupiter, and Saturn.

Description: The methods for predicting the strength of composite laminates with fastener holes are reviewed for the cases of unloaded as well as loaded holes. Three categories of methods are discussed: (1) the hole-boundary stress (HBS) methods which compare the peak stresses and the laminate strength; (2) the linear elastic fracture (LEFM) methods which assume cracks emanating from the hole and predict failure when the stress intensity factor for a crack equals the laminate fracture toughness; and (3) two different two-parameter (TP) methods: the average-stress (ASTP) method, which predicts failure when the average stress over a distance from the hole equals the laminate strength; and the point-stress (PSTP) Whitney-Nuismer (1974) approach, which predicts failure when the stress near the hole equals the material strength. The PSTP is relatively accurate and is by far the most widely used of all the prediction methods. However, for large ranges of variables and cases, the current strength predictions may not always be accurate.

Description: Optimality criteria methods take advantage of some concepts as those of statically determinate or indeterminate structures, and certain variational principles of structural dynamics, to develop efficient algorithms for the sizing of structures that are subjected to stiffness-related constraints. Some of the methods and iterative strategies developed over the last decade for calculations of the Lagrange multipliers in stressand displacement-limited problems, as well as for satisfying the appropriate optimality criterion, are discussed. The application of these methods are illustrated by solving problems with stress and displacement constraints.

Description: The stiffened composite structural panel analysis and sizing code designated 'PASCO' encompasses both the generality required for the exploitation of composite materials' design flexibility and an accurate buckling analysis for the detection of complex buckling modes. PASCO can accordingly design for buckling, frequency, material strength, and panel stiffness requirements. Attention is given to an additional thermal loading design capability. Design studies illustrate the importance of the multiple load condition capability when thermal loads are present.

Description: The results of observations of solar flares, obtained by means of the UV Spectrophotometer and Polarimeter (UVSP) instrument on board the Solar Maximum Mission satellite are summarized and discussed. The results are grouped into three main topics: (1) plasma diagnostics in the flare transition zone plasmas, (2) spatial and temporal evolutions of the UV and hard X-ray bursts, and (3) energy release processes in the impulsive phase. The methods of spectral UV analysis, comparison with the hard X-ray burst results, and the interpretation of the results are summarized. It is concluded that the energy release processes in the flare phase can be best interpreted in terms of multiple large and small interloops interacting with each other either mechanically or inductively. Furthermore, the majority of impulsive UV and hard X-ray bursts occurs in small compact loops with high densities and transition-zone temperatures. The impulsive hard X-ray and UV bursts are emitted by nonthermal particles, accelerated by the loop interactions and impinging on the footpoints of these loops.

Description: A brief review is given of non-flare investigations using data obtained by the Ultraviolet Spectrometer and Polarimeter on the Solar Maximum Mission. The major topics described are sunspot research including magnetic field measurements, oscillations, and models; mass motions in quiet and active regions including steady flows and acoustic waves; and prominence research including physical conditions, dynamics, and mass motions around prominences. Also discussed are studies of UV bursts, the formation of the Cl I line at 1351 A, ozone in the terrestrial atmosphere, and active regions using correlated observations from other instruments on the spacecraft or on the ground.

Description: A detailed comparison is presented of the predicted eigenfrequencies of twisted rotating plates as obtained by using two different shape functions. Primarily, rotating twisted plates of two different aspect ratios and two different thickness ratios are considered. The effects of rotation are included by using a 'stress smoothing' technique when calculating the augmented stiffness matrix. In addition, the effects of Coriolis acceleration, contributions from membrane behavior, setting angle and sweep angle are considered. The effects of geometric nonlinearity are briefly discussed. Finally, results of a brief study of cambered plates are presented.

Description: Delamination growth in compressively loaded composite laminates was studied analytically and experimentally. The configuration used was a laminate with an across-the-width delamination. An approximate super-position stress analysis was developed to quantify the effects of various geometric, material, and load parameters on mode 1 and mode 2 strain energy release rates G sub 1 and G sub 2, respectively. Calculated values of G sub 1 and G sub 2 were then compared with measured cyclic delamination growth rates to determine the relative importance of G sub 1 and G sub 2. High growth rates were observed only when G sub 1 was large. However, slow growth was observed even when G sub 1 was negligibly small. This growth was apparently due to a large value of G sub 2.

Description: New results obtained with the Clark Lake multifrequency radioheliograph at meter-decameter wavelengths and from satellite multifrequency directive observations at hectometer and kilometer wavelengths are reviewed. Evidence is presented that type III electrons propagate in dense coronal streamers and that frequently observed microbursts (presumably type III) at meter-decameter wavelengths are due to plasma radiation. Observations of hectometer and kilometer type III radio storms which reveal information about active region structures, the interplanetary magnetic field configuration, and solar wind acceleration are discussed. Kilometer type II bursts and interactions between type III electrons and interplanetary shocks are examined, and some new results on shock-associated events are presented.

Description: An attempt is made to develop a self-consistent model which accounts for the line and continuum data generated by the three X-ray imaging instruments on the SMM satellite. The intensities measured covered the 4-500 kV energy range. The model is based on a differential emission measure and electron beam parameters and is used to predict absolute signals detected by the 15 channels of the SMM sensors. Consideration is given to the thermal contribution, instrumental characteristics, thin target excitation and thick target bremsstrahlung. In comparison with data from a flare event on June 29, 1980, model predictions provide a good fit, including the identification of hard electrons with a 5.3 index during the impulsive phase.

Description: Solar Maximum Mission observations have been used to study the origin and amount of energy, mechanism of storage and release, and conditions for the occurrence of solar flares, and some results of these studies as they pertain to homologous flares are briefly discussed. It was found that every set of flares produced 'rafales' of homologous flares, i.e., two, three, four, or more flares separated in time by an hour or less. No great changes in macroscopic photospheric patterns were observed during these flaring periods. A quantitative brightness parameter of the relation between homologous flares is defined. Scale changes detected in the dynamic spectrum of flare sites are in good agreement with a theoretical suggestion by Sturrock. Statistical results for different homologous flare active regions show the existence in homologous flaring areas of a 'pivot' of previous filaments interpreted as a signature of an anomaly in the solar rotation.

Description: The observational difficulties of obtaining the magnetic field distribution in the chromosphere and corona of the sun has led to methods of extending photospheric magnetic mesurements into the solar atmosphere by mathematical procedures. A new approach to this problem presented here is that a constant alpha force-free field can be uniquely determined from the tangential components of the measured photospheric flux alone. The vector magnetographs now provide measurements of both the solar photospheric tangential and the longitudinal magnetic field. This paper presents derivations for the computation of the solar magnetic field from these type of measurements. The fields considered are assumed to be a constant alpha force-free fields or equivalent, producing vanishing Lorentz forces. Consequently, magnetic field lines and currents are related by a constant and hence show an identical distribution. The magnetic field above simple solar regions are described from the solution of the field equations.

Description: The Solar Maximum Mission (SMM) spacecraft has provided high time resolution observational data regarding the soft X-ray emission from solar-flare plasma during 1980. The present investigation is concerned with the characteristics of a soft X-ray flare and the energetics of the impulsive phase on the basis of the data collected with the aid of two of the instruments on board the SMM, taking into account the Hard X-ray Burst Spectrometer (HXRBS) and the Bent Crystal Spectrometer (BCS). Attention is given to an analysis of soft X-ray flare spectra, the relative motion of the soft X-ray sources, the phenomenology of the soft X-ray flare, energy and mass transport during the impulsive phase, and energy deposition in the chromosphere during evaporation.

Description: The mass ejections of 1 September, 1980 are studied from observations obtained with the MSDP spectrograph and with the Ultraviolet Spectrometer and Polarimeter aboard the Solar Maximum Mission satellite. The analysis is focused on observations in the chromospheric H-alpha line and the transition region C IV 1548 A line. It is noted that cold and hot material had the same projection, although the upward C IV velocity structure was more extended than the H-alpha one. It is shown that the observed contrast of the H-alpha absorbing structure can be interpreted in terms of a dynamic cloud model overlying the chromosphere. Radial velocities of 25-30 km/s and -40 km/s are estimated for the first and second phases of ejection, respectively.

Description: A special three-dimensional element based on the total Lagrangian description of the motion of a layered anisotropic composite medium is developed, validated and employed to analyze laminated anisotropic composite shells. The element contains the following features: geometric nonlinearity, dynamic (transient) behavior and arbitrary lamination scheme and lamina properties. Numerical results of nonlinear bending, natural vibration, and transient response are presented to illustrate the capabilities of the element.

Description: Based on the principal component analysis technique and evidence for a 22-yr double-sunspot cycle periodicity. The time series of sunspot numbers is represented as a sum of mutually orthogonal eigenvectors in the time domain. It is shown that the first two eigenvectors account for about 90 percent of the cumulative 'signal power,' and that this is sufficient for reconstruction of the raw data curve. It is also noted that the second eigenvector behaves as the time derivative of the first, and that a phase-plane plot of these eigenvectors (i.e. a plot of a variable vs. its rate of change) suggests that the sun's sunspot cycle is driven by an oscillator; the implication is that, embedded within the sun, a chronometer is at work (e.g. Dicke, 1979).

Description: Analytical and observational data are presented to show that the lower transition zone, a 100 km thick region at 10,000-200,000 K between the solar chromosphere and corona, is heated by local electric currents. The study was spurred by correlations between the enhanced atmospheric heating and magnetospheric flux in the chromospheric network and active regions. Field aligned current heated flux loops are asserted to mainly reside in and make up most of the transition region. It is shown that thermal conduction from the sides of hot gas columns generated by the current dissipation is the source of the observed temperature distribution in the transition regions.

Description: Selected plasma parameters observed by Pioneer 10 and Pioneer 11 between launch (1972 and 1973) and the end of 1979 are used to find the large-scale radial structure of the solar wind. Comparison of data from the two spacecraft is used to separate temporal from spatial variations. The average bulk speed is found to remain constant at about 430 km/s, with stream structure still evident, though of diminished amplitude, at 20.5 AU (Pioneer 10's distance by the end of 1979). Proton density, flux, pressure, and kinetic energy flux are found to have radial profiles consistent with 1/R-squared. Proton temperatures decrease as R to the -0.6 power, too slowly for an adiabatic expansion.

Description: Aluminum alloys of types 1100, 3003, 5052, 6061, and 2024 were used to study the sensitivity of the acousto-elastic constant to changes in the microstructure. Results show that there is a strong relationship between the acousto-elastic constants and the yield strength and hardness. This relationship depends on whether the alloy is strain hardened or precipitation hardened. In strain hardened alloys, the constants increase as the amount of solid solution is decreased, while the behavior is the opposite in precipitation hardened alloys.

Description: Observations of oscillations in filaments, which are cold magnetic structures suspended in the corona, are used to test theories of coronal heating by acoustic waves high in the atmosphere. Long lived stationary motions observed in filaments suggest exchanges of material between them and the surrounding corona. A coupling between upflows in filaments and subphotospheric convective motions, dragging the field lines of the prominence magnetic support, is possible.

Description: Using the 'dynamo theory' method to predict solar activity, a value for the smoothed sunspot number of 109 + or - 20 is obtained for solar cycle 22. The predicted cycle is expected to peak near December, 1990 + or - 1 year. Concommitantly, F(10.7) radio flux is expected to reach a smoothed value of 158 + or - 18 flux units. Global mean exospheric temperature is expected to reach 1060 + or - 50 K and global total average total thermospheric density at 400 km is expected to reach 4.3 x 10 to the -15th gm/cu cm + or - 25 percent.

Description: Considering the solar atmosphere in its entirety, that large-scale nonadiabatic processes (MHD waves, and so on) effectively tend to yield an 'equalization of temperature', such that the atmospheric temperature is limited to the base temperature associated with its heat source. This conjecture suggests that (1) the chromospheric temperature is limited by the granulation base temperature (10,000 K), (2) the spicule temperatures are limited by the base temperature (100,000 K) where the supergranular cells form, and (3) the quiet coronal temperature is less than or equal to the convection zone base temperature (2,000,000 K). Thermodynamical arguments are provided which may serve to augment the detailed heating models wherein large-scale mechanical energy is transported into the solar atmosphere involving MHD waves, current dissipation, and other nonthermal processes.

Description: The two-dimensional, linear hydrodynamics of quiet solar and umbral model atmospheres in a plane-parallel, adiabatic approximation are investigated. The 5.5-8.5 mHz oscillations observed in umbral chromospheres and transition regions are interpreted as acoustic waves propagating parallel, or nearly parallel, to the temperature gradient. These waves are not totally internally reflected by the steep temperature gradient and, thus, are not trapped. Partial reflections, however, are effective in modulating the transmission as a function of frequency. The resonant transmission mechanism of Zugzda, Locans, and Staude (1983) is found to produce a spectrum of resonances in the transmission of acoustic waves in any atmosphere with a temperature minimum. Since the observed umbral oscillations display power in only a narrow range of frequencies, characteristics of the umbral models, wave propagation, and observations that would tend to suppress the higher frequency resonances are examined.

Description: The possibility of the existence in soft X-ray flare plasmas of conditions that result in a steady state departure of ion abundances from ionization equilibrium values is considered. The observed flare plasma is assumed to be a result of many small 'elementary bursts' that occur on time scales comparable to the ionization and recombination times of highly ionized atoms of iron and calcium. Specific models are adopted, the time-dependent equations for ion abundances are solved numerically, and X-ray line intensities and line ratios are computed and averaged over the effective time of a single burst. The computed results are compared to observed variations for a number of different line ratios. Although the behavior of certain line ratios can be explained in the context of the burst models considered in this paper, the behavior of the set of all the available line ratios cannot be explained in this manner. The observed departures of line ratios from equilibrium values that can be explained in terms of a burst scenario can also be accounted for by uncertainties in the atomic physics.

Description: A cylindrical axisymmetric tearing mode model for solar flares is investigated numerically. Large magnetic energy release only occurs when there are at least two mode rational surfaces in the current-carrying plasma.

Description: The effect of nonlinearity on the higher harmonic control (HHC) of helicopter vibration is investigated using a nonlinear aeroelastic simulation. A nonlinear solution is proposed which relates the HHC inputs to vibration outputs on the basis of a Volterra functional series. The Volterra series solution is shown to reduce to a vector polynomial equation relating HHC inputs to vibration outputs at any harmonic frequency. The nonlinear transfer relationship is identified from a nonlinear vibration analysis computer program, and the identification model is examined in detail. Improvements to current HHC algorithms are presented, and several Kalman filter divergence problems are quantified.

Description: Prograde and retrograde sectoral oscillations of the sun have been observed so as to determine frequency differences produced by rotation. Oscillations in the frequency range 2.1-3.7 mHz and with spherical harmonic degrees from 1 to 100 have been identified. Average frequency shifts due to rotation in a sidereal reference frame are found to range from a high of about 660 nHz at degree 1 to a low of about 423 nHz at degree 6, rising to about 471 nHz at degree 100. These results indicate that most of the sun's volume rotates at a rate close to that of the surface, but also that the energy-generating core may rotate more rapidly than the surface.

Description: The frequency difference between prograde and retrograde sectoral solar oscillations is analyzed to determine the rotation rate of the solar interior, assuming no latitudinal dependence. Much of the solar interior rotates slightly less rapidly than the surface, while the innermost part apparently rotates more rapidly. The resulting solar gravitational quadrupole moment is J2 = (1.7 + or - 0.4) x 10 to the -7th and provides a negligible contribution to current planetary tests of Einstein's theory of general relativity.

Description: A statistical analysis of intensities and radial velocities of several solar plage filaments (i.e. prominences seen on the solar disk) observed at disk center is presented. Intensity and radial (= vertical) velocity maps were derived from simultaneous 2D measurements of the H-alpha chromospheric line with the Multichannel Subtractive Double Pass spectrograph operating on the Meudon solar tower, and also obtained from 2D recordings of the C IV transition-zone line (1548 A) with the UV Spectrometer and Polarimeter aboard the Solar Maximum Mission satellite. A good correlation (around 0.5) is found between intensities in both lines, as well as between velocities. Persistent upflows are measured in both lines at the filament location. The mean vertical velocities are respectively 0.5 km/s in H-alpha and 5.6 km/s in C IV. The analysis of mass fluxes suggests that C IV upflows occur in the transition region around prominences rather than below, in the chromosphere-corona transition zone.

Description: The Active Cavity Radiometer Irradiance Monitor (ACRIM) of the Solar Maximum Mission satellite measures the radiant power emitted by the sun in the direction of the earth and has worked flawlessly since 1980. The main motivation for ACRIM's use to measure the solar constant is the determination of the extent to which this quantity's variations affect earth weather and climate. Data from the solar minimum of 1986-1987 is eagerly anticipated, with a view to the possible presence of a solar cycle variation in addition to that caused directly by sunspots.

Description: The heating of minor ions in solar flares by wave-wave-particle interaction with Langmuir waves, or ion acoustic waves, can be described by a diffusion equation in velocity-space for the particle distribution function. The dependence of the heating on the ion charge and mass, and on the composition of the plasma, is examined in detail. It is found that the heat mechanisms proposed by Ibragimov and Kocharov cannot account for the enhanced abundances of heavy elements in the solar cosmic rays.

Description: The paper presents a broad range of complementary observations (SMM and ground-based) of the onset and impulsive phase of the fairly large (1B, M1.2) but simple two-ribbon flare which occurred at 19:15 UT on November 1, 1980 in the northern part of the active region Boulder No. AR2776. It is found that the overall magnetic field configuration in which the flare occurred was a fairly simple, closed arch containing nonpotential substructure; the flare occurred spontaneously within the arch (it was not triggered by emerging magnetic flux). The two major spikes of the impulsive energy release are examined, and the three immediate products of this energy release are discussed.

Description: A general crack opening stress equation is presented which may be used to correlate crack growth rate data for various materials and thicknesses, under constant amplitude loading, once the proper constraint factor has been determined. The constraint factor, alpha, is a constraint on tensile yielding; the material yields when the stress is equal to the product of alpha and sigma. Delta-K (LEFM) is plotted against rate for 2024-T3 aluminum alloy specimens 2.3 mm thick at various stress ratios. Delta-K sub eff was plotted against rate for the same data with alpha = 1.8; the rates correlate well within a factor of two.

Description: The transient response of a 0.4 by 0.6 cm rectangular bore rail accelerator was analyzed by a three dimensional finite element code. The copper rail deflected to a peak value of 0.08 mm in compression and then oscillated at an amplitude of 0.02 mm. Simultaneously the insulating side wall of glass fabric base, epoxy resin laminate (G-10) was compressed to a peak value of 0.13 mm and rebounded to a steady state in extension. Projectile pinch or blowby due to the rail extension or compression, respectively, can be identified by examining the time history of the rail displacement. The effect of blowby was most significant at the side wall characterized by mm size displacement in compression. Dynamic stress calculations indicate that the G-10 supporting material behind the rail is subjected to over 21 MPa at which the G-10 could fail if the laminate was not carefully oriented. Results for a polycarbonate resin (Lexan) side wall show much larger displacements and stresses than for G-10. The tradeoff between the transparency of Lexan and the mechanical strength of G-10 for sidewall material is obvious. Displacement calculations from the modal method are smaller than the results from the direct integration method by almost an order of magnitude, because the high frequency effect is neglected. Previously announced in STAR as N83-35412

Description: Transient behavior of flare-associated solar wind in the nonradial open field region is numerically investigated, taking into account the thermal and dynamical coupling between the chromosphere and the corona. A realistic steady solar wind is constructed which passes through the inner X-type critical point in the rapidly diverging region. The wind speed shows a local maximum at the middle, O-type, critical point. The wind's density and pressure distributions decrease abruptly in the rapidly diverging region of the flow tube. The transient behavior of the wind following flare energy deposition includes ascending and descending conduction fronts. Thermal instability occurs in the lower corona, and ascending material flows out through the throat after the flare energy input ceases. A local density distribution peak is generated at the shock front due to the pressure deficit just behind the shock front.

Description: The unified constitutive theories for application to typical isotropic cast nickel base supperalloys used for air-cooled turbine blades were evaluated. The specific modeling aspects evaluated were: uniaxial, monotonic, cyclic, creep, relaxation, multiaxial, notch, and thermomechanical behavior. Further development of the constitutive theories to model thermal history effects, refinement of the material test procedures, evaluation of coating effects, and verification of the models in an alternate material will be accomplished in a follow-on for this base program.

Description: A state-of-the-art review of applicable constitutive models with selection of two for detailed comparison with a wide range of experimental tests was conducted. The experimental matrix contained uniaxial and biaxial tensile, creep, stress relaxation, and cyclic fatigue tests at temperatures to 1093 C and strain rates from .0000001 to .001/sec. Some nonisothermal cycles will also be run. The constitutive models will be incorporated into the MARC finite element structural analysis program with a demonstration computation made for advanced turbine blade configuration. In the code development work, particular emphasis is being placed on developing efficient integration algorithms for the highly nonlinear and stiff constitutive equations. Another area of emphasis is the appropriate and efficient methodology for determing constitutive constants from a minimum extent of experimental data.

Description: About 1% of all sudden ionospheric disturbances (SIDs) observed at the Panska Ves Observatory (Czechoslovakia), were found to be not of solar-XUV origin. Among them, the very rare SWF events (observed at L = 2.4) of corpuscular origin are the most interesting. The IMF sector structure effects in the midlatitude lower ionosphere are minor in comparison with effects of solar flares, geomagnetic storms, etc. There are two basic types of effects. The first type is a disturbance, best developed in geomagnetic activity, and observed in the night-time ionosphere. It can be interpreted as a response to sector structure related changes of geomagnetic (= magnetospheric) activity. The other type is best developed in the tropospheric vorticity area index and is also observed in the day-time ionosphere in winter. This effect is quietening in the ionosphere as well as troposphere. While the occurrence of the former type is persistent in time, the latter is severely diminished in some periods. All the stratosphere, the 10-mb level temperature and height above Berlin-Tempelhof do not display any observable IMF section structure effect.

Description: Radio wave absorption data on 1539 kHz for the summer period of 1978 to 1980 are considered in relation to variations of solar X-ray and L-alpha radiation. It is shown that under non-flare conditions L-alpha dominates in controlling absorption and that X-rays contribute about 10% to the total absorption. Optimum regression equations show that absorption is proportional to the m-th power of ionizing flux where m 1. The role of correcting L-alpha values, measured by the AE-E satellite, is discussed.

Description: Static aerothermoelastic design/analysis of axial-flow compressors, modal flutter analysis of axial-flow turbomachines, forced vibration analysis of rotating cyclic structures and modal flutter analysis of advanced turbopropellers with highly swept blades are four new capabilities developed and implemented in NASTRAN Level 17.7. The contents, applicability and usefulness of these capabilities which were developed and documented under the sponsorship of NASA's Lewis Research Center are discussed. Overall flowcharts and selected examples are presented.

Description: COSMIC/NASTRAN was converted to the CRAY computer systems. The CRAY version is currently available and provides users with access to all of the machine independent source code of COSMIC/NASTRAN. Future releases of COSMIC/NASTRAN will be made available on the CRAY soon after they are released by COSMIC.

Description: Computer aided engineering systems are invaluable tools in performing NASTRAN finite element analysis. These techniques are implemented in both the pre-processing and post-processing phases of the NASTRAN analysis. The finite element model development, or pre-processing phase, was automated with a computer aided modeling program called Supertabl, and the review and interpretation of the results of the NASTRAN analysis, or post-processing phase, was automated with a computer aided plotting program called Output Display. An intermediate program, Nasplot, which was developed in-house, has also helped to cut down on the model checkout time and reduce errors in the model. An interface has been established between the finite element computer aided engineering system and the Learjet computer aided design system whereby data can be transferred back and forth between the two. These systems have significantly improved productivity and the ability to perform NASTRAN analysis in response to product development requests.

Description: A methodology that supports forced transient response dynamic solutions when both static and kinetic friction effects are included in a structural system model is described. Modifications that support this type of nonlinear transient response solution are summarized for the transient response dynamics (TRD) NASTRAN module. An overview of specific modifications for the NASTRAN processing subroutines, INITL, TRD1C, and TRD1D, are described with further details regarding inspection of nonlinear input definitions to define the type of nonlinear solution required, along with additional initialization requirements and specific calculation subroutines to successfully solve the transient response problem. The extension of the basic NASTRAN nonlinear methodology is presented through several stages of development to the point where constraint equations and residual flexibility effects are introduced into the finite difference Newmark-Beta recurrsion formulas. Particular emphasis is placed on cost effective solutions for large finite element models such as the Space Shuttle with friction degrees of freedom between the orbiter and payloads mounted in the cargo bay. An alteration to the dynamic finite difference equations of motion is discussed, which allows one to include friction effects at reasonable cost for large structural systems such as the Space Shuttle. Data are presented to indicate the possible impact of transient friction loads to the payload designer for the Space Shuttle. Transient response solution data are also included, which compare solutions without friction forces and those with friction forces for payloads mounted in the Space Shuttle cargo bay. These data indicate that payload components can be sensitive to friction induced loads.

Description: The results of a study to determine the accuracy of the three dimensional solid elements available in NASTRAN for predicting displacements is presented. Of particular interest in the study is determining how to effectively use solid elements in analyzing thick optical mirrors, as might exist in a large telescope. Surface deformations due to thermal and gravity loading can be significant contributors to the determination of the overall optical quality of a telescope. The study investigates most of the solid elements currently available in either COSMIC or MSC NASTRAN. Error bounds as a function of mesh refinement and element aspect ratios are addressed. It is shown that the MSC solid elements are, in general, more accurate than their COSMIC NASTRAN counterparts due to the specialized numerical integration used. In addition, the MSC elements appear to be more economical to use on the DEC VAX 11/780 computer.

Description: Some innovative techniques applicable to sensitivity analysis of discretized structural systems are reviewed. These techniques include a finite-difference step-size selection algorithm, a method for derivatives of iterative solutions, a Green's function technique for derivatives of transient response, a simultaneous calculation of temperatures and their derivatives, derivatives with respect to shape, and derivatives of optimum designs with respect to problem parameters. Computerized implementations of sensitivity analysis and applications of sensitivity derivatives are also discussed. Finally, some of the critical needs in the structural sensitivity area are indicated along with Langley plans for dealing with some of these needs.

Description: In developing the constitutive equations an interdisciplinary approach is being pursued. Specifically, both metallurgical and continuum mechanics considerations are recognized in the formulation. Experiments will be utilized to both explore general qualitative features of the material behavior that needs to be modeled and to provide a means of assessing the validity of the equations being developed. The model under development explicitly recognizes crystallographic slip on the individual slip systems. This makes possible direct representation of specific slip system phenomena. The present constitutive formulation takes the anisotropic creep theory and incorporates two state variables into the model to account for the effect of prior inelastic deformation history on the current rate-dependent response of the material.

Description: A brief overview of the status of the fatigue and fracture programs is given. The programs involve the development of appropriate analytic material behavior models for cyclic stress-strain-temperature-time/cyclic crack initiation, and cyclic crack propagation. The underlying thrust of these programs is the development and verification of workable engineering methods for the calculation, in advance of service, of the local cyclic stress-strain response at the critical life governing location in hot section compounds, and the resultant crack initiation and crack growth lifetimes.

Description: For economic reasons, the steam drive for a large compressor was replaced by a large synchronous electric motor. Due to the resulting large increase in mass and because the unit was mounted on a steel frame approximately 18 feet above ground level, it was deemed necessary to determine if a steady state or transient vibration problem existed. There was a definite possibility that a resonant or near resonant condition could be encountered. The ensuing analysis, which led to some structural changes as the analysis proceeded, did not reveal any major steady state vibration problems. However, the analysis did indicate that the system would go through several natural frequencies of the support structure during start-up and shutdown. This led to the development of special start-up and shutdown procedures to minimize the possibility of exciting any of the major structural modes. A coast-down could result in significant support structure and/or equipment damage, especially under certain circumstances. In any event, dynamic field tests verified the major analytical results. The unit has now been operating for over three years without any major vibration problems.

Description: Seismic analysis is a technique which pertains to loading described in terms of boundary accelerations. Earthquake shocks to buildings is the type of excitation which usually comes to mind when one hears the word seismic, but this technique also applied to a broad class of acceleration excitations which are applied at the base of a structure such as vibration shaker testing or shocks to machinery foundations. Four different solution paths are available in NASTRAN for seismic analysis. They are: Direct Seismic Frequency Response, Direct Seismic Transient Response, Modal Seismic Frequency Response, and Modal Seismic Transient Response. This capability, at present, is invoked not as separate rigid formats, but as pre-packaged ALTER packets to existing RIGID Formats 8, 9, 11, and 12. These ALTER packets are included with the delivery of the NASTRAN program and are stored on the computer as a library of callable utilities. The user calls one of these utilities and merges it into the Executive Control Section of the data deck to perform any of the four options are invoked by setting parameter values in the bulk data.

Description: A user's guide to the COSMIC/NASTRAN free field input for the Bulk Data section of the NASTRAN program is proposed. The free field input is designed to be user friendly and the user is not forced out of the computer system due to input errors. It is easy to use, with only a few simple rules to follow. A stand alone version of the COSMIC/NASTRAN free field input is also available. The use of free field input is illustrated by a number of examples.

Description: The annual static loading tests were made on the assembled wing/fuselage structure of the RP-1 sailplane. Simple wing bending and wing bending/torsion tests up to about four G's were made, using sand bag loading on the inverted aircraft. Bending a torsion deflections were measured and compared to those obtained in previous years' tests. It appears that the structure is retaining its original strength and stiffness very well. This aircraft has been disassembled and stored on the Jonsson Engineering Center under ambient, indoor conditions of temperature and humidity.

Description: The appearance of solar magnetic flux is discussed; both the mechanisms and phenomena associated with flux emergence are considered. The dynamics of solar surface magnetic flux are addressed, including both the transport and structure of magnetic-flux-carrying elements as seen in the chromosphere and corona. The disappearance of magnetic flux from the surface of the Sun is also discussed. Standard solar models representing generally accepted views are considered along with observations which seem to fall outside the scope of these models.

Description: Progress has been made on the design fabrication and installation of a tension carry-through linkage between the aft parts of the port and starboard wings to resist forward bending deflection on the RP-2 sailplane. In addition, towhook/release mechanisms, landing gear and wheel brake assembly, compound curved sandwich composite fuselage shells and an enclosed trailer to house and transport the complete disassembled aircraft, have all undergone various stages of design, fabrication, installation and testing.

Description: Mean frequencies, amplitudes, and linewidths for the solar 5 min p mode oscillations of degree 0, 1, and 2 have been obtained from approx. 280 days of SMM-ACRIM total irradiance data. The frequencies are in good agreement with measurements obtained from velocity data. The amplitudes of the modes lie along a well defined envelope of power vs. frequency, which peaks at 3.1 mHz and has a width of 0.7 mHz (FWHM). The r.m.s. amplitude of the highest peak in the spectrum (n=21, l=1) is approx. 3 ppm of the total flux. The linewidths of the narrowest l=O modes are approx. 1 micro Hz (FWHM). A broad continuum of power caused both by solar surface granulation and by instrumental noise interferes with the analysis of 5 min modes. The continuum spectral power in a 1 micro Hz band near 3 mHz corresponds to an apparent r.m.s. variation of approx. 0.5 parts per million of the mean solar flux.

Description: Solar UV irradiance variations with solar activity are examined using a three component model of the CaII K chromospheric emission. This model, developed from ground based observations of the location, area and relative intensity of CaII K plage, in conjunction with measurements throughout solar cycle 21 of the full disc CaII K emission, includes the contributions to the ultraviolet flux from both plage and active network emission. The model successfully replicates changes in the Lyman alpha flux related to the 27 day rotation of solar plage, outbreaks (or rounds) of activity over periods of a year or more, and the growth and accumulation of active regions over the eleven year solar activity cycles. Estimates of the magnitude of the solar cycle variability of the UV emission between 200 and 300 nm are presented but cannot currently be verified by available observations.

Description: The principal uncertainties in the equation of state involve the treatment of pressure ionization, the Debye-Huckel coulomb corrections, and the treatment of many-particle interaction effects. It is found that, for the lowest degree modes (l between 0 and 3), the terms and procedures used in the equation of state which deal with these uncertainties introduce changes in the frequencies which are less than 4 micro Hz. Recently, Shibahashi, Noels and Gabriel (1983) published solar eigenfrequencies using a theory with an equation of state improved with respect to the theory used earlier by Shibahashi and Osaki (1981). Their comparison between the two sets of results suggested that uncertainties in the frequencies as large as 10 micro Hz could be caused by the equation of state. It is felt that since the entire effect of the uncertain terms is only 4 micro Hz and since the uncertainties are only a fraction of each term, the 10 micro Hz changes found by Shibahashi et al must be a consequence of differences between the earlier and later calculations in areas other than the equation of state.

Description: The nature of the fine structure of high order, low degree five minute period solar oscillations following from various postulated forms of spherical rotation is predicted. The first and second order effects of rotation are included.

Description: The overstability of acoustic modes trapped in the solar convection zone is studied with mechanical and thermal effects of turbulence included, in an approximate manner, through the eddy transport coefficients. Many of these acoustic modes are found to be overstable with the most rapidly growing modes occupying a region centered around 3.2 mHz and spread over a wide range of length-scales. The numerical results are in reasonable accord with the observed power-spectrum of the five-minute oscillations of intermediate and high degree. The oscillations are probably driven by a simultaneous operation of the kappa-mechanism and the turbulent conduction (convective Cowling) mechanism, the dominant contribution to the generation of self-excited acoustic waves arising from the convective Cowling mechanism.

Description: The decay rate for many of the low degree p modes observed as 5 minute oscillations of the Sun is discussed and decay rates for modes ranging from radial (l=0) to the nonradial ones with l=5 for overtones 10 through 28. Parameters needed for he solar model are given and the hydrogen mass fraction composition structure is shown as well as the structure given by Christensen-Dalsgaard (1982) for an evolved solar model. The special equation of state and opacity table with X=0.74 for the hydrogen mass fraction in the outer 0.40 of the mass needs slightly more hydrogen in the central regions than obtained by Christensen-Dalsgaard in order to give a complete and consistent model. The difference in helium production between these two models is about 10%, meaning that the total energy radiated by the Sun during its lifetime thus far agrees satisfactorily with accurately calculated evolution sequences.