Results for “continuous finite element” — Page 13

NASA NTRS2003

High-Fidelity Micromechanics Model Enhanced for Multiphase Particulate Materials

This 3-year effort involves the development of a comprehensive micromechanics model and a related computer code, capable of accurately estimating both the average response and the local stress and strain fields in the individual phases, assuming both elastic and inelastic behavior. During the first year (fiscal year 2001) of the investigation, a version of the model called the High-Fidelity Generalized Method of Cells (HFGMC) was successfully completed for the thermo-inelastic response of continuously reinforced multiphased materials with arbitrary periodic microstructures (refs. 1 and 2). The model s excellent predictive capability for both the macroscopic response and the microlevel stress and strain fields was demonstrated through comparison with exact analytical and finite element solutions. This year, HFGMC was further extended in two technologically significant ways. The first enhancement entailed the incorporation of fiber/matrix debonding capability into the two-dimensional version of HFGMC for modeling the response of unidirectionally reinforced composites such as titanium matrix composites, which exhibit poor fiber/matrix bond. Comparison with experimental data validated the model s predictive capability. The second enhancement entailed further generalization of HFGMC to three dimensions to enable modeling the response of particulate-reinforced (discontinuous) composites in the elastic material behavior domain. Next year, the three-dimensional version will be generalized to encompass inelastic effects due to plasticity, viscoplasticity, and damage, as well as coupled electromagnetothermomechanical (including piezoelectric) effects.

Pindera, Marek-Jerzy↗

NASA NTRS2013

Next Generation Lightweight Mirror Modeling Software

The advances in manufacturing techniques for lightweight mirrors, such as EXELSIS deep core low temperature fusion, Corning's continued improvements in the Frit bonding process and the ability to cast large complex designs, combined with water-jet and conventional diamond machining of glasses and ceramics has created the need for more efficient means of generating finite element models of these structures. Traditional methods of assembling 400,000 + element models can take weeks of effort, severely limiting the range of possible optimization variables. This paper will introduce model generation software developed under NASA sponsorship for the design of both terrestrial and space based mirrors. The software deals with any current mirror manufacturing technique, single substrates, multiple arrays of substrates, as well as the ability to merge submodels into a single large model. The modeler generates both mirror and suspension system elements, suspensions can be created either for each individual petal or the whole mirror. A typical model generation of 250,000 nodes and 450,000 elements only takes 5-10 minutes, much of that time being variable input time. The program can create input decks for ANSYS, ABAQUS and NASTRAN. An archive/retrieval system permits creation of complete trade studies, varying cell size, depth, and petal size, suspension geometry with the ability to recall a particular set of parameters and make small or large changes with ease. The input decks created by the modeler are text files which can be modified by any editor, all the key shell thickness parameters are accessible and comments in deck identify which groups of elements are associated with these parameters. This again makes optimization easier. With ANSYS decks, the nodes representing support attachments are grouped into components; in ABAQUS these are SETS and in NASTRAN as GRIDPOINT SETS, this make integration of these models into large telescope or satellite models possible.

Arnold, William↗

NASA NTRS2013

Next-Generation Lightweight Mirror Modeling Software

The advances in manufacturing techniques for lightweight mirrors, such as EXELSIS deep core low temperature fusion, Corning's continued improvements in the Frit bonding process and the ability to cast large complex designs, combined with water-jet and conventional diamond machining of glasses and ceramics has created the need for more efficient means of generating finite element models of these structures. Traditional methods of assembling 400,000 + element models can take weeks of effort, severely limiting the range of possible optimization variables. This paper will introduce model generation software developed under NASA sponsorship for the design of both terrestrial and space based mirrors. The software deals with any current mirror manufacturing technique, single substrates, multiple arrays of substrates, as well as the ability to merge submodels into a single large model. The modeler generates both mirror and suspension system elements, suspensions can be created either for each individual petal or the whole mirror. A typical model generation of 250,000 nodes and 450,000 elements only takes 5-10 minutes, much of that time being variable input time. The program can create input decks for ANSYS, ABAQUS and NASTRAN. An archive/retrieval system permits creation of complete trade studies, varying cell size, depth, and petal size, suspension geometry with the ability to recall a particular set of parameters and make small or large changes with ease. The input decks created by the modeler are text files which can be modified by any editor, all the key shell thickness parameters are accessible and comments in deck identify which groups of elements are associated with these parameters. This again makes optimization easier. With ANSYS decks, the nodes representing support attachments are grouped into components; in ABAQUS these are SETS and in NASTRAN as GRIDPOINT SETS, this make integration of these models into large telescope or satellite models possible

Arnold, William R., Sr.↗

NASA NTRS2013

Next Generation Lightweight Mirror Modeling Software

The advances in manufacturing techniques for lightweight mirrors, such as EXELSIS deep core low temperature fusion, Corning's continued improvements in the Frit bonding process and the ability to cast large complex designs, combined with water-jet and conventional diamond machining of glasses and ceramics has created the need for more efficient means of generating finite element models of these structures. Traditional methods of assembling 400,000 + element models can take weeks of effort, severely limiting the range of possible optimization variables. This paper will introduce model generation software developed under NASA sponsorship for the design of both terrestrial and space based mirrors. The software deals with any current mirror manufacturing technique, single substrates, multiple arrays of substrates, as well as the ability to merge submodels into a single large model. The modeler generates both mirror and suspension system elements, suspensions can be created either for each individual petal or the whole mirror. A typical model generation of 250,000 nodes and 450,000 elements only takes 5-10 minutes, much of that time being variable input time. The program can create input decks for ANSYS, ABAQUS and NASTRAN. An archive/retrieval system permits creation of complete trade studies, varying cell size, depth, and petal size, suspension geometry with the ability to recall a particular set of parameters and make small or large changes with ease. The input decks created by the modeler are text files which can be modified by any editor, all the key shell thickness parameters are accessible and comments in deck identify which groups of elements are associated with these parameters. This again makes optimization easier. With ANSYS decks, the nodes representing support attachments are grouped into components; in ABAQUS these are SETS and in NASTRAN as GRIDPOINT SETS, this make integration of these models into large telescope or satellite models easier.

Arnold, William R., Sr.↗

NASA NTRS2011

A Deployable 4 Meter 180 to 680 GHz Antenna for the Scanning Microwave Limb Sounder

The Scanning Microwave Limb Sounder (SMLS) is a space-borne heterodyne radiometer which will measure pressure, temperature and atmospheric constituents from thermal emission between 180 and 680 GHz. SMLS, planned for the Global Atmospheric Composition Mission of the NRC Decadal Survey, uses a novel toric Cassegrain antenna to perform both elevation and azimuth scanning. These provide better horizontal and temporal resolution and coverage than were possible with elevation-only scanning at typical Low-Earth orbit spacing in the two previous MLS satellite instruments. Development of the SMLS antenna was the focus of a 2006 Small Business Innovative Research (SBIR) program whose phase II culminated in the fabrication and thermal stability testing of a composite demonstration model of the SMLS primary reflector. This reflector has the full 4m height and 1/3 the width planned for flight. An Instrument Incubator Program (IIP) titled "A deployable 4 Meter 180 to 680 GHz antenna for the Scanning Microwave Limb Sounder" continues development of the SMLS antenna with the study of 5 topics: 1) detailed mathematical modeling of the antenna patterns from which we simulate geophysical parameter retrievals in order to establish FOV performance requirements; 2) thorough correlation of finite element model predictions with measurements made on the SBIR reflector. We will again measure deformations of this reflector, under more flight-like thermal gradients, using higher precision metrology techniques available in a new large-aperture facility at JPL; 3) fabrication of a full-width primary reflector whose asbuilt surface figure will better meet the figure requirements of SMLS than did the SBIR reflector; 4) integration of the primary with other reflectors, and with residual front ends built in a 2007 IIP, in a breadboard antenna; and finally 5) RF testing of the breadboard on a Near Field Range at JPL. We report on significant progress in 3 areas of the current IIP: development of the mathematical model to predict SMLS antenna patterns and their application in a preliminary set of geophysical retrievals; the correlation between surface deformation predicted by finite element models and measurement in the 2009 isothermal stability tests of the SBIR, with implications for the thermal gradient testing to be performed at JPL; and aspects of the conceptual design of the full-width primary reflector to be fabricated and tested in the 2nd and 3rd years of the IIP.

heterodyne radiometer↗

NASA NTRSDate not supplied

Development of a C-ELS Specimen-Based Numerical Benchmark for Mode II Delamination and Assessment of Two VCCT-Based Propagation Strategies

A finite element (FE) benchmark example inspired by the calibrated end-loaded split (C-ELS) specimen is developed and used to assess the performance of delamination propagation capabilities based on linear elastic fracture mechanics (LEFM). The C-ELS specimen has the advantage of a longer region of stable delamination propagation compared to the existing mode II benchmark case. The new benchmark example may therefore provide a better assessment tool by enabling more stable crack growth in regions further away from the boundary conditions or load application. First, a benchmark result is created manually using two-dimensional finite element models of the C-ELS specimen with different delamination lengths. Second, the performance of the virtual crack closure technique (VCCT) delamination propagation capabilities in the Abaqus/Standard®1 FE code and the recently developed Progressive Release eXplicit-VCCT (PRX-VCCT) method are assessed by comparing the results to the benchmark case. Two examples with different starter delamination lengths are studied. A shorter starter length is chosen to create a scenario with unstable delamination propagation. A longer delamination encourages stable delamination propagation. Detailed results from three-dimensional analyses with aligned and misaligned meshes and two levels of mesh refinement are provided. In general, good agreement can be achieved between the results obtained from the quasi-static propagation analysis and the benchmark analysis. Numerical artifacts including anomalous unreleased nodes in the crack wake and zig-zag crack fronts occur for propagation analyses using Abaqus/Standard VCCT. In comparison, continuous, smooth, delamination fronts are observed for PRX-VCCT. The use of the benchmark case to assess different VCCT-based propagation strategies illustrates the value of establishing benchmark cases.

Composite Materials↗

NASA NTRS2002

A Passive Magnetic Bearing Flywheel

A 100 percent passive magnetic bearing flywheel rig employing no active control components was designed, constructed, and tested. The suspension clothe rotor was provided by two sets of radial permanent magnetic bearings operating in the repulsive mode. The axial support was provided by jewel bearings on both ends of the rotor. The rig was successfully operated to speeds of 5500 rpm, which is 65 percent above the first critical speed of 3336 rpm. Operation was not continued beyond this point because of the excessive noise generated by the air impeller and because of inadequate containment in case of failure. Radial and axial stiffnesses of the permanent magnetic bearings were experimentally measured and then compared to finite element results. The natural damping of the rotor was measured and a damping coefficient was calculated.

Siebert, Mark↗

NASA NTRS2004

Thermo-Elastic Triangular Sandwich Element for the Complete Stress Field Based on a Single-Layer Theory

This study presents a new triangular finite element for modeling thick sandwich panels, subjected to thermo-mechanical loading, based on a {3,2}-order single-layer plate theory. A hybrid energy functional is employed in the derivation of the element because of a C interelement continuity requirement. The single-layer theory is based on five weighted-average field variables arising from the cubic and quadratic representations of the in-plane and transverse displacement fields, respectively. The variations of temperature and distributed loading acting on the top and bottom surfaces are non-uniform. The temperature varies linearly through the thickness.

Das, M.↗

NASA NTRS1992

Evolution of attached and detached slabs and their associated mantle dynamics

Over the two years of the NASA grant, this project has produced a significant amount of research results related to the plate subduction process and the surface crustal deformation at convergent boundaries (i.e., above subduction zones). While some research objectives are completely accomplished, other research tasks remain active and continue to be investigated at present. A steady state analytic thermal model for subducting slabs was used to examine the torques acting on a descending slab. It is found that gravitational torque vanishes when a slab is dipping either vertically or horizontally, unlike previous studies indicating that the magnitude of gravitational torque decreases as dip angle increases. Subsequently, a new time-dependent, analytic thermal model for a subducting slab was developed. The new model enables us to study transient phenomena associated with plate subduction analytically. On the basis of this model, the nature of slab dip angles was evaluated. Slab dip angles are found to be transient features. As they penetrate into the mantle and increase their lengths, the associated gravitational torque also increases resulting in a downward pulling of the slab to the steeper dip angle. This is especially true once a slab penetrates the olivine-spinel phase boundary at about 400 km depth. However, if the phase transformation does not follow the equilibrium condition, the gravitational torque may have a different behavior. This problem was investigated. Except for fast descending slabs, non-equilibrium phase transformation can only slow down the transient increase of slab dip angles discussed earlier. Its effect is not sufficiently strong to reverse the downward pulling for most of the slabs. However, when slabs subducting at 10 cm/yr or faster, a sufficient amount of metastable olivine can exist beneath 400 km. Because of its low density compared with the surrounding spinel, an upward buoyancy is produced resulting in an upward bending of the slab and possibly an upward rotation of the slab such that smaller dip angles are formed. Seismic studies of the Japanese Slab seem to support this interpretation. The development of oroclinal geometries at convergent boundaries was also examined to study plate obduction which is an important ingredient to the initiation of plate subduction. Although the study suggests that surface features are better modeled by block models, the large scale deformation can be adequately studied by viscous models. Such a model is now under development to complete our original objective to study the initiation of plate subduction. Finally, a three-dimensional, finite element, spherical convective model is developed to study dynamic plate subductions. The model development is now complete and it is being tested to ensure its proper operation. The model is able to generate convection results with a viscosity contrast of about 100. Our research continues to push the viscosity contrast to a level that is appropriate for a subducting slab.

Hsui, Albert T.↗

NASA NTRS1998

GPS Monitoring of Surface Change During and Following the Fortuitous Occurrence of the M(sub w) = 7.3 Landers Earthquake in our Network

Accomplishments: (1) Continues GPS monitoring of surface change during and following the fortuitous occurrence of the M(sub w) = 7.3 Landers earthquake in our network, in order to characterize earthquake dynamics and accelerated activity of related faults as far as 100's of kilometers along strike. (2) Integrates the geodetic constraints into consistent kinematic descriptions of the deformation field that can in turn be used to characterize the processes that drive geodynamics, including seismic cycle dynamics. In 1991, we installed and occupied a high precision GPS geodetic network to measure transform-related deformation that is partitioned from the Pacific - North America plate boundary northeastward through the Mojave Desert, via the Eastern California shear zone to the Walker Lane. The onset of the M(sub w) = 7.3 June 28, 1992, Landers, California, earthquake sequence within this network poses unique opportunities for continued monitoring of regional surface deformation related to the culmination of a major seismic cycle, characterization of the dynamic behavior of continental lithosphere during the seismic sequence, and post-seismic transient deformation. During the last year, we have reprocessed all three previous epochs for which JPL fiducial free point positioning products available and are queued for the remaining needed products, completed two field campaigns monitoring approx. 20 sites (October 1995 and September 1996), begun modeling by development of a finite element mesh based on network station locations, and developed manuscripts dealing with both the Landers-related transient deformation at the latitude of Lone Pine and the velocity field of the whole experiment. We are currently deploying a 1997 observation campaign (June 1997). We use GPS geodetic studies to characterize deformation in the Mojave Desert region and related structural domains to the north, and geophysical modeling of lithospheric behavior. The modeling is constrained by our existing and continued GPS measurements, which will provide much needed data on far-field strain accumulation across the region and on the deformational response of continental lithosphere during and following a large earthquake, forming the basis for kinematic and dynamic modeling of secular and seismic-cycle deformation. GPS geodesy affords both regional coverage and high precision that uniquely bear on these problems.

Miller, M. Meghan↗

NASA NTRS1989

Fatigue delamination onset prediction in tapered composite laminates

Tapered (0 deg) laminates of S2/CE9000 and S2/SP250 glass/epoxies, and IM6/1827I graphite/epoxy were tested in cyclic tension. The specimens usually showed some initial stable delaminations in the tapered region, but these did not affect the stiffness of the specimens, and loading was continued until the specimens either delaminated unstably, or reached 10(exp 6) to 2 x 10(exp 7) million cycles with no unstable delamination. The final unstable delamination originated at the junction of the thin and tapered regions. A finite-element model was developed for the tapered laminate with and without the initial stable delaminations observed in the tests. The analysis showed that for both cases the most likely place for an opening (Mode 1) delamination to originate is at the junction of the taper and thin regions. For each material type, the models were used to calculate the strain energy release rate, G, associated with delaminations originating at that junction and growing either into the thin region or tapered region. For the materials tested, cyclic G(sub Imax) values from DCB tests were used with the maximum strain energy release rates calculated from the finite-element analysis to predict the onset of unstable delamination at the junction as a function of fatigue cycles. The predictions were compared to experimental values of maximum cyclic load as a function of cycles to unstable delamination from fatigue tests in tapered laminates. For the IM6/1827I and S2/SP250 laminates, the predictions agreed very well with the test data. Predicted values for the S2/CE9000 were conservative compared to the test data.

Murri, Gretchen Bostaph↗

NASA NTRS1991

Fatigue delamination onset prediction in unidirectional tapered laminates

Tapered (0 deg) laminates of S2/CE9000 and S2/SP250 glass/epoxies, and IM6/1827I graphite/epoxy were tested in cyclic tension. The specimens usually showed some initial stable delaminations in the tapered region, but these did not affect the stiffness of the specimens, and loading was continued until the specimens either delaminated unstably, or reached 10(exp 6) to 2 x 10 (exp 7) million cycles with no unstable delamination. The final unstable delamination originated at the junction of the thin and tapered regions. A finite-element model was developed for the tapered laminate with and without the initial stable delaminations observed in the tests. The analysis showed that for both cases the most likely place for an opening (Mode I) delamination to originate is at the junction of the taper and thin regions. For each material type, the models were used to calculate the strain energy release rate, G, associated with delaminations originating at that junction and growing either into the thin region or tapered region. For the materials tested, cyclic G(sub Imax) values from DCB tests were used with the maximum strain energy release rates calculated from the finite-element analysis to predict the onset of unstable delamination at the junction as a function of fatigue cycles. The predictions were compared to experimental values of maximum cyclic load as a function of cycles to unstable delamination from fatigue tests in tapered laminates. For the IM6/1827I and S2/SP250 laminates, the predictions agreed very well with the test data. Predicted values for the S2/CE9000 were conservative compared to the test data.

Murri, Gretchen B.↗

NASA NTRS2009

Using Covariance Analysis to Assess Pointing Performance

A Pointing Covariance Analysis Tool (PCAT) has been developed for evaluating the expected performance of the pointing control system for NASA s Space Interferometry Mission (SIM). The SIM pointing control system is very complex, consisting of multiple feedback and feedforward loops, and operating with multiple latencies and data rates. The SIM pointing problem is particularly challenging due to the effects of thermomechanical drifts in concert with the long camera exposures needed to image dim stars. Other pointing error sources include sensor noises, mechanical vibrations, and errors in the feedforward signals. PCAT models the effects of finite camera exposures and all other error sources using linear system elements. This allows the pointing analysis to be performed using linear covariance analysis. PCAT propagates the error covariance using a Lyapunov equation associated with time-varying discrete and continuous-time system matrices. Unlike Monte Carlo analysis, which could involve thousands of computational runs for a single assessment, the PCAT analysis performs the same assessment in a single run. This capability facilitates the analysis of parametric studies, design trades, and "what-if" scenarios for quickly evaluating and optimizing the control system architecture and design.

Bayard, David↗

NASA NTRS1990

Micromechanics analysis of space simulated thermal deformations and stresses in continuous fiber reinforced composites

Space simulated thermally induced deformations and stresses in continuous fiber reinforced composites were investigated with a micromechanics analysis. The investigation focused on two primary areas. First, available explicit expressions for predicting the effective coefficients of thermal expansion (CTEs) for a composite were compared with each other, and with a finite element (FE) analysis, developed specifically for this study. Analytical comparisons were made for a wide range of fiber/matrix systems, and predicted values were compared with experimental data. The second area of investigation focused on the determination of thermally induced stress fields in the individual constituents. Stresses predicted from the FE analysis were compared to those predicted from a closed-form solution to the composite cylinder (CC) model, for two carbon fiber/epoxy composites. A global-local formulation, combining laminated plate theory and FE analysis, was used to determine the stresses in multidirectional laminates. Thermally induced damage initiation predictions were also made.

Bowles, David E.↗

NASA NTRS1996

First-Order System Least Squares for the Stokes Equations, with Application to Linear Elasticity

Following our earlier work on general second-order scalar equations, here we develop a least-squares functional for the two- and three-dimensional Stokes equations, generalized slightly by allowing a pressure term in the continuity equation. By introducing a velocity flux variable and associated curl and trace equations, we are able to establish ellipticity in an H(exp 1) product norm appropriately weighted by the Reynolds number. This immediately yields optimal discretization error estimates for finite element spaces in this norm and optimal algebraic convergence estimates for multiplicative and additive multigrid methods applied to the resulting discrete systems. Both estimates are uniform in the Reynolds number. Moreover, our pressure-perturbed form of the generalized Stokes equations allows us to develop an analogous result for the Dirichlet problem for linear elasticity with estimates that are uniform in the Lame constants.

Cai, Z.↗

NASA NTRS2003

Simulating Fatigue Crack Growth in Spiral Bevel Pinion

This project investigates computational modeling of fatigue crack growth in spiral bevel gears. Current work is a continuation of the previous efforts made to use the Boundary Element Method (BEM) to simulate tooth-bending fatigue failure in spiral bevel gears. This report summarizes new results predicting crack trajectory and fatigue life for a spiral bevel pinion using the Finite Element Method (FEM). Predicting crack trajectories is important in determining the failure mode of a gear. Cracks propagating through the rim may result in catastrophic failure, whereas the gear may remain intact if one tooth fails and this may allow for early detection of failure. Being able to predict crack trajectories is insightful for the designer. However, predicting growth of three-dimensional arbitrary cracks is complicated due to the difficulty of creating three-dimensional models, the computing power required, and absence of closed- form solutions of the problem. Another focus of this project was performing three-dimensional contact analysis of a spiral bevel gear set incorporating cracks. These analyses were significant in determining the influence of change of tooth flexibility due to crack growth on the magnitude and location of contact loads. This is an important concern since change in contact loads might lead to differences in SIFs and therefore result in alteration of the crack trajectory. Contact analyses performed in this report showed the expected trend of decreasing tooth loads carried by the cracked tooth with increasing crack length. Decrease in tooth loads lead to differences between SIFs extracted from finite element contact analysis and finite element analysis with Hertz contact loads. This effect became more pronounced as the crack grew.

Ural, Ani↗

NASA NTRS1996

Coupled mixed-field laminate theory and finite element for smart piezoelectric composite shell structures

Mechanics for the analysis of laminated composite shells with piezoelectric actuators and sensors are presented. A new mixed-field laminate theory for piezoelectric shells is formulated in curvilinear coordinates which combines single-layer assumptions for the displacements and a layerwise representation for the electric potential. The resultant coupled governing equations for curvilinear piezoelectric laminates are described. Structural mechanics are subsequently developed and an 8-node finite-element is formulated for the static and dynamic analysis of adaptive composite structures of general laminations containing piezoelectric layers. Evaluations of the method and comparisons with reported results are presented for laminated piezoelectric-composite plates, a closed cylindrical shell with a continuous piezoceramic layer and a laminated composite semi-circular cantilever shell with discrete cylindrical piezoelectric actuators and/or sensors.

Saravanos, Dimitris A.↗

NASA NTRS1996

Parallel Newton-Krylov-Schwarz algorithms for the transonic full potential equation

We study parallel two-level overlapping Schwarz algorithms for solving nonlinear finite element problems, in particular, for the full potential equation of aerodynamics discretized in two dimensions with bilinear elements. The overall algorithm, Newton-Krylov-Schwarz (NKS), employs an inexact finite-difference Newton method and a Krylov space iterative method, with a two-level overlapping Schwarz method as a preconditioner. We demonstrate that NKS, combined with a density upwinding continuation strategy for problems with weak shocks, is robust and, economical for this class of mixed elliptic-hyperbolic nonlinear partial differential equations, with proper specification of several parameters. We study upwinding parameters, inner convergence tolerance, coarse grid density, subdomain overlap, and the level of fill-in in the incomplete factorization, and report their effect on numerical convergence rate, overall execution time, and parallel efficiency on a distributed-memory parallel computer.

Cai, Xiao-Chuan↗