Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,369 records · Page 76

Frequency-Domain Analysis of Diffusion-Cooled Hot-Electron Bolometer Mixers

A new theoretical model is introduced to describe heterodyne mixer conversion efficiency and noise (from thermal fluctuation effects) in diffusion-cooled superconducting hot-electron bolometers. The model takes into account the non-uniform internal electron temperature distribution generated by Wiedemann-Franz heat conduction, and accepts for input an arbitrary (analytical or experimental) superconducting resistance-versus- temperature curve. A non-linear large-signal solution is solved iteratively to calculate the temperature distribution, and a linear frequency-domain small-signal formulation is used to calculate conversion efficiency and noise. In the small-signal solution the device is discretized into segments, and matrix algebra is used to relate the heating modulation in the segments to temperature and resistance modulations. Matrix expressions are derived that allow single-sideband mixer conversion efficiency and coupled noise power to be directly calculated. The model accounts for self-heating and electrothermal feedback from the surrounding bias circuit.

Skalare, A.↗

Fracture Mechanisms For SiC Fibers And SiC/SiC Composites Under Stress-Rupture Conditions at High Temperatures

The successful application of SiC/SiC ceramic matrix composites as high-temperature structural materials depends strongly on maximizing the fracture or rupture life of the load-bearing fiber and matrix constituents. Using high-temperature data measured under stress-rupture test conditions, this study examines in a mechanistic manner the effects of various intrinsic and extrinsic factors on the creep and fracture behavior of a variety of SiC fiber types. It is shown that although some fiber types fracture during a large primary creep stage, the fiber creep rate just prior to fracture plays a key role in determining fiber rupture time (Monkman-Grant theory). If it is assumed that SiC matrices rupture in a similar manner as fibers with the same microstructures, one can develop simple mechanistic models to analyze and optimize the stress-rupture behavior of SiC/SiC composites for applied stresses that are initially below matrix cracking.

DiCarlo, James A.↗

Thermal Dispersion Within a Porous Medium Near a Solid Wall

The regenerator is a key component to Stirling cycle machine efficiency. Typical regenerators are of sintered fine wires or layers of fine-wire screens. Such porous materials are contained within solid-waH casings. Thermal energy exchange between the regenerator and the casing is important to cycle performance for the matrix and casing would not have the same axial temperature profile in an actual machine. Exchange from one to the other may allow shunting of thermal energy, reducing cycle efficiency. In this paper, temperature profiles within the near-wall region of the matrix are measured and thermal energy transport, termed thermal dispersion, is inferred. The data show how the wall affects thermal transport. Transport normal to the mean flow direction is by conduction within the solid and fluid and by advective transport within the matrix. In the near-wall region, both may be interrupted from their normal in-core pattern. Solid conduction paths are broken and scales of advective transport are damped. An equation is presented which describes this change for a wire screen mesh. The near-wall layer typically acts as an insulating layer. This should be considered in design or analysis. Effective thermal conductivity within the core is uniform. In-core transverse thermal effective conductivity values are compared to direct and indirect measurements reported elsewhere and to 3D numerical simulation results, computed previously and reported elsewhere. The 3-D CFD model is composed of six cylinders in cross flow, staggered in arrangement to match the dimensions and porosity of the matrix used in the experiments. The commercial code FLUENT is used to obtain the flow and thermal fields. The thermal dispersion and effective thermal conductivities for the matrix are computed from the results.

Simon, T.↗

Swift/BAT Calibration and Spectral Response

The Burst Alert Telescope (BAT) aboard NASA#s Swift Gamma-Ray Burst Explorer is a large coded aperture gamma-ray telescope consisting of a 2.4 m (8#) x 1.2 m (4#) coded aperture mask supported 1 meter above a 5200 square cm area detector plane containing 32,768 individual 4 mm x 4 mm x 2 mm CZT detectors. The BAT is now completely assembled and integrated with the Swift spacecraft in anticipation of an October 2004 launch. Extensive ground calibration measurements using a variety of radioactive sources have resulted in a moderately high fidelity model for the BAT spectral and photometric response. This paper describes these ground calibration measurements as well as related computer simulations used to study the efficiency and individual detector properties of the BAT detector array. The creation of a single spectral response model representative of the fully integrated BAT posed an interesting challenge and is at the heart of the public analysis tool #batdrmgen# which computes a response matrix for any given sky position within the BAT FOV. This paper will describe the batdrmgen response generator tool and conclude with a description of the on-orbit calibration plans as well as plans for the future improvements needed to produce the more detailed spectral response model that is required for the construction of an all-sky hard x-ray survey.

Parsons, A.↗

Dynamic Analysis and Control of Lightweight Manipulators with Flexible Parallel Link Mechanisms

The objective is the theoretical analysis and the experimental verification of dynamics and control of a two link flexible manipulator with a flexible parallel link mechanism. Nonlinear equations of motion of the lightweight manipulator are derived by the Lagrangian method in symbolic form to better understand the structure of the dynamic model. The resulting equation of motion have a structure which is useful to reduce the number of terms calculated, to check correctness, or to extend the model to higher order. A manipulator with a flexible parallel link mechanism is a constrained dynamic system whose equations are sensitive to numerical integration error. This constrained system is solved using singular value decomposition of the constraint Jacobian matrix. Elastic motion is expressed by the assumed mode method. Mode shape functions of each link are chosen using the load interfaced component mode synthesis. The discrepancies between the analytical model and the experiment are explained using a simplified and a detailed finite element model.

Lee, Jeh Won↗

Historical (1850‐2014) Aerosol Evolution and Role on Climate Forcing Using the GISS Modele2.1 Contribution to CMIP6

The Earth’s climate is rapidly changing. Over the past centuries, aerosols, via their ability to absorb or scatter solar radiation and alter clouds, played an important role in counterbalancing some of the greenhouse gas (GHG) caused global warming. The multi-century anthropogenic aerosol cooling effect prevented present-day climate from reaching even higher surface air temperatures and subsequent more dramatic climate impacts. Trends in aerosol concentrations and optical depth show that in many polluted regions such as Europe and the United States of America, aerosol precursor emissions decreased back to levels of the 1950s. More recent polluting countries such as China may have reached a turning point in recent years as well, while India still follows an upward trend. Here we study aerosol trends in the CMIP6 simulations of the GISS ModelE2.1 climate model using a fully coupled atmosphere composition onfiguration, including interactive gas-phase chemistry, and either an aerosol microphysical (MATRIX) or a mass-based (OMA) aerosol module. Results show that whether global aerosol radiative forcing is already declining depends on the aerosol scheme used. Using the aerosol microphysical scheme, where the aerosol system reacts more strongly to the trend in sulfur dioxide (SO2) emissions, global peak direct aerosol forcing was reached in the 1980’s, whereas the mass-based scheme simulates peak direct aerosol forcing around 2010.

Aerosol forcing↗

Analysis of thermal stresses in polymer matrix composites

In the present micromechanics study of the thermally-induced stress field which emerges at very low temperatures in the matrix and fibers of continuous fiber-reinforced polymer composites, the results obtained by a specialized FEM analysis are compared to an analytical solution of the composite cylinder model for several composite materials. Attention is given to the influence of microstructural geometry and constituent properties; it is established that the distributions and magnitudes of the thermally-induced stresses are affected by the assumed microstructural geometry of the sample. While matrix stresses are not a strong function of fiber properties, the temperature dependence of the matrix properties significantly affects the magnitudes of thermally-induced matrix stresses.

Bowles, David E.↗

Fracture and crack growth in orthotropic laminates. Part 2: Experimental determination of internal damage growth in unidirectional boron/aluminum composite laminates

The fracture behavior of unidirectional boron/aluminum composite laminates is investigated in order to verify the results of mathematical models. These models predict the stresses and displacements of fibers and the amount of damage growth in a center-notched lamina as a function of the applied remote stress and the matrix and fiber material properties. The damage may take the form of longitudinal yielding and splitting in the matrix as well as stable transverse damage consisting of broken fibers and matrix yielding ahead of the notch. A brittle lacquer coating is used to detect the yielding in the matrix while X-ray techniques are used to detemine the number of broken fibers in the laminate. The notched strengths and the amounts of damage found in the experimental specimens agree well with those predicted by the mathematical model.

Goree, J. G.↗

A proposed method for enhanced eigen-pair extraction using finite element methods: Theory and application

The paper covers two distinct parts: theory and application. The goal of this work was the reduction of model size with an increase in eigenvalue/vector accuracy. This method is ideal for the condensation of large truss- or beam-type structures. The theoretical approach involves the conversion of a continuum transfer matrix beam element into an 'Exact' dynamic stiffness element. This formulation is implemented in a finite element environment. This results in the need to solve a transcendental eigenvalue problem. Once the eigenvalue is determined the eigenvectors can be reconstructed with any desired spatial precision. No discretization limitations are imposed on the reconstruction. The results of such a combined finite element and transfer matrix formulation is a much smaller FEM eigenvalue problem. This formulation has the ability to extract higher eigenvalues as easily and as accurately as lower eigenvalues. Moreover, one can extract many more eigenvalues/vectors from the model than the number of degrees of freedom in the FEM formulation. Typically, the number of eigenvalues accurately extractable via the 'Exact' element method are at least 8 times the number of degrees of freedom. In contrast, the FEM usually extracts one accurate (within 5 percent) eigenvalue for each 3-4 degrees of freedom. The 'Exact' element results in a 20-30 improvement in the number of accurately extractable eigenvalues and eigenvectors.

Jara-Almonte, J.↗

Motion models in attitude estimation

Attitude estimator use observations from different times to reduce the effects of noise. If the vehicle is rotating, the attitude at one time needs to be propagated to that at another time. If the vehicle measures its angular velocity, attitude propagating entails integrating a rotational kinematics equation only. If a measured angular velocity is not available, torques can be computed and an additional rotational dynamics equation integrated to give the angular velocity. Initial conditions for either of these integrations come from the estimation process. Sometimes additional quantities, such as gyro and torque parameters, are also solved for. Although the partial derivatives of attitude with respect to initial attitude and gyro parameters are well known, the corresponding partial derivatives with respect to initial angular velocity and torque parameters are less familiar. They can be derived and computed numerically in a way that is analogous to that used for the initial attitude and gyro parameters. Previous papers have demonstrated the feasibility of using dynamics models for attitude estimation but have not provided details of how each angular velocity and torque parameters can be estimated. This tutorial paper provides some of that detail, notably how to compute the state transition matrix when closed form expressions are not available. It also attempts to put dynamics estimation in perspective by showing the progression from constant to gyro-propagated to dynamics-propagated attitude motion models. Readers not already familiar with attitude estimation will find this paper an introduction to the subject, and attitude specialists may appreciate the collection of heretofore scattered results brought together in a single place.

Chu, D.↗

A spectral element approach to wave motion in layered solids

A matrix methodology similar to that of the finite element method is developed for the analysis of stress waves in layered solids. Because the mass distribution is modeled exactly, the approach gives the exact frequency response of each layer. The fast Fourier transform and Fourier series are used for inversion to the time/space domain. The impact of a structured medium with multiple layers is used to demonstrate the method. Comparison with existing propagator and direct global matrix methods show the present approach to be computationally more efficient.

Rizzi, S. A.↗

Stiffness and Fatigue Life Estimator for Polymer Composite Laminates Using Machine Learning

Machine learning (ML) models are increasingly being used in many engineering fields due to the advancements in ML algorithms and availability of high-speed computing power. One of the most popular ML class of models is artificial neural networks (ANN). ML is increasingly being used in the design and analysis of composite materials and structures, specifically in the constitutive modeling of composite materials with the focus on greatly accelerating multiscale analyses of composite materials and structures through development of surrogate models. Towards that end, Python-based neural nets have been developed to predict initial stiffness and fatigue life of an eight-ply symmetric polymer matrix composite laminate. Two types of neural networks, a Multilayer Perceptron (MLP) and a Recurrent Neural Network (RNN), have been established. Results show that both neural net type algorithms can provide an excellent estimate of initial laminate stiffness as well as fatigue life of eight-ply symmetric polymer matrix composite laminates (PMCs). RNNs are better able to capture the shape of the fatigue curve of a laminate. The resulting tool and GUI can be very useful for system level studies to obtain an estimate of desired properties and life of PMC composite laminates. Further, the associated surrogate models can also be used in composite multiscale analyses to replace the actual physics-based calculations at lower scales and thereby significantly increase the computational efficiency of such analyses and thus make micromechanics-based multiscale analyses a viable industrial tool for large scale structural problems.

multiscale analysis↗

Stiffness and Fatigue Life Estimator for Polymer Composite Laminates Using Machine Learning

Machine learning (ML) models are increasingly being used in many engineering fields due to the advancements in ML algorithms and availability of high-speed computing power. One of the most popular ML class of models is artificial neural networks (ANN). ML is increasingly being used in the design and analysis of composite materials and structures, specifically in the constitutive modeling of composite materials with the focus on greatly accelerating multiscale analyses of composite materials and structures through development of surrogate models. Towards that end, Python-based neural nets have been developed to predict initial stiffness and fatigue life of an eight-ply symmetric polymer matrix composite laminate. Two types of neural networks, a Multilayer Perceptron (MLP) and a Recurrent Neural Network (RNN), have been established. Results show that both neural net type algorithms can provide an excellent estimate of initial laminate stiffness as well as fatigue life of eight-ply symmetric polymer matrix composite laminates (PMCs). RNNs are better able to capture the shape of the fatigue curve of a laminate. The resulting tool and GUI can be very useful for system level studies to obtain an estimate of desired properties and life of PMC composite laminates. Further, the associated surrogate models can also be used in composite multiscale analyses to replace the actual physics-based calculations at lower scales and thereby significantly increase the computational efficiency of such analyses and thus make micromechanics-based multiscale analyses a viable industrial tool for large scale structural problems.

multiscale analysis↗

Cost analysis of composite fan blade manufacturing processes

The relative manufacturing costs were estimated for large high technology fan blades prepared by advanced composite fabrication methods using seven candidate materials/process systems. These systems were identified as laminated resin matrix composite, filament wound resin matrix composite, superhybrid solid laminate, superhybrid spar/shell, metal matrix composite, metal matrix composite with a spar and shell, and hollow titanium. The costs were calculated utilizing analytical process models and all cost data are presented as normalized relative values where 100 was the cost of a conventionally forged solid titanium fan blade whose geometry corresponded to a size typical of 42 blades per disc. Four costs were calculated for each of the seven candidate systems to relate the variation of cost on blade size. Geometries typical of blade designs at 24, 30, 36 and 42 blades per disc were used. The impact of individual process yield factors on costs was also assessed as well as effects of process parameters, raw materials, labor rates and consumable items.

Stelson, T. S.↗

Design, fabrication and test of graphite/polyimide composite joints and attachments for advanced aerospace vehicles

Bonded and bolted designs are presented for each of four major attachment types. Prepreg processing problems are discussed and quality control data are given for lots 2W4604, 2W4632 and 2W4643. Preliminary design allowables test results for tension tests and compression tests of laminates are included. The final small specimen test matrix is defined and the configuration of symmetric step-lap joint specimens are shown. Finite element modeling studies of a double lap joint were performed to evaluate the number of elements required through the adhesive thickness to assess effects of various joint parameters on stress distributions. Results of finite element analyses assessing the effect of an adhesive fillet on the stress distribution in a double lap joint are examined.

Skoumal, D. E.↗

Application of a systematic finite-element model modification technique to dynamic analysis of structures

A systematic finite-element model modification technique has been applied to two small problems and a model of the main wing box of a research drone aircraft. The procedure determines the sensitivity of the eigenvalues and eigenvector components to specific structural changes, calculates the required changes and modifies the finite-element model. Good results were obtained where large stiffness modifications were required to satisfy large eigenvalue changes. Sensitivity matrix conditioning problems required the development of techniques to insure existence of a solution and accelerate its convergence. A method is proposed to assist the analyst in selecting stiffness parameters for modification.

Robinson, J. C.↗

Extended aeroelastic analysis for helicopter rotors with prescribed hub motion and blade appended penduluum vibration absorbers

The mathematical development for the expanded capabilities of the G400 rotor aeroelastic analysis was examined. The G400PA expanded analysis simulates the dynamics of all conventional rotors, blade pendulum vibration absorbers, and the higher harmonic excitations resulting from prescribed vibratory hub motions and higher harmonic blade pitch control. The methodology for modeling the unsteady stalled airloads of two dimensional airfoils is discussed. Formulations for calculating the rotor impedance matrix appropriate to the higher harmonic blade excitations are outlined. This impedance matrix, and the associated vibratory hub loads, are the rotor dynamic characteristic elements for use in the simplified coupled rotor/fuselage vibration analysis (SIMVIB). Updates to the development of the original G400 theory, program documentation, user instructions and information are presented.

Bielawa, R. L.↗