Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element analysis”

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,081 records · Page 60

Decentralized control experiments on NASA's flexible grid

Methods arising from the area of decentralized control are emerging for analysis and control synthesis for large flexible structures. In this paper the control strategy involves a decentralized model reference adaptive approach using a variable structure control. Local models are formulated based on desired damping and response time in a model-following scheme for various modal configurations. Variable structure controllers are then designed employing co-located angular rate and position feedback. In this scheme local control forces the system to move on a local sliding mode in some local error space. An important feature of this approach is that the local subsystem is made insensitive to dynamical interactions with other subsystems once the sliding surface is reached. Experiments based on the above have been performed for NASA's flexible grid experimental apparatus. The grid is designed to admit appreciable low-frequency structural dynamics, and allows for implementation of distributed computing components, inertial sensors, and actuation devices. A finite-element analysis of the grid provides the model for control system design and simulation; results of several simulations are reported on here, and a discussion of application experiments on the apparatus is presented.

Ozguner, U.↗

Repair Concepts as Design Constraints of a Stiffened Composite PRSEUS Panel

A design and analysis of a repair concept applicable to a stiffened thin-skin composite panel based on the Pultruded Rod Stitched Efficient Unitized Structure is presented. The concept is a bolted repair using metal components, so that it can easily be applied in the operational environment. The damage scenario considered is a midbay-to-midbay saw-cut with a severed stiffener, flange and skin. In a previous study several repair configurations were explored and their feasibility confirmed but refinement was needed. The present study revisits the problem under recently revised design requirements and broadens the suite of loading conditions considered. The repair assembly design is based on the critical tension loading condition and subsequently its robustness is verified for a pressure loading case. High fidelity modeling techniques such as mesh-independent definition of compliant fasteners, elastic-plastic material properties for metal parts and geometrically nonlinear solutions are utilized in the finite element analysis. The best repair design is introduced, its analysis results are presented and factors influencing the design are assessed and discussed.

Przekop, Adam↗

Continuation of tailored composite structures of ordered staple thermoplastic material

The search for the cost effective composite structure has motivated the investigation of several approaches to develop composite structure from innovative material forms. Among the promising approaches is the conversion of a planar sheet to components of complex curvature through sheet forming or stretch forming. In both cases, the potential for material stretch in the fiber direction appears to offer a clear advantage in formability over continuous fiber systems. A framework was established which allows the simulation of the anisotropic mechanisms of deformation of long discontinuous fiber laminates wherein the matrix phase is a viscous fluid. Predictions for the effective viscosities of a hyper-anisotropic medium consisting of collimated, discontinuous fibers suspended in viscous matrix were extended to capture the characteristics of typical polymers including non-Newtonian behavior and temperature dependence. In addition, the influence of fiber misorientation was also modeled by compliance averaging to determine ensemble properties for a given orientation distribution. A design tool is presented for predicting the effect of material heterogeneity on the performance of curved composite beams such as those used in aircraft fuselage structures. Material heterogeneity can be induced during manufacturing processes such as sheet forming and stretch forming of thermoplastic composites. This heterogeneity can be introduced in the form of fiber realignment and spreading during the manufacturing process causing radial and tangential gradients in material properties. Two analysis procedures are used to solve the beam problems. The first method uses separate two-dimensional elasticity solutions for the stresses in the flange and web sections of the beam. The separate solutions are coupled by requiring that forces and displacements match section boundaries. The second method uses an approximate Rayleigh-Ritz technique to find the solutions for more complex beams. Analyses are performed for curved beams of various cross-sections loaded in pure bending and with a uniform distributed load. Preliminary results show that the geometry of the beam dictates the effect of heterogeneity on performance. The role of heterogeneity is larger in beams with a small average radius-to-depth ration, R/t, where R is the average radius of the beam and t is the difference between the inside and outside radii. Results of the anlysis are in the form of stresses and displacements and are compared to both mechanics of materials and numerical solutions obtained using finite element analysis.

Santare, Michael H.↗

Thermal Design, Test and Analysis of PharmaSat, a Small Class D Spacecraft with a Biological Experiment

Small spacecraft have been increasing in popularity because of their low cost, short turnaround and relative efficiency. In the past, small spacecraft have been primarily used for technology demonstrations, but advances in technology have made the miniaturization of space science possible [1,2]. PharmaSat is a low cost, small three cube size spacecraft, with a biological experiment on board, built at NASA (National Aeronautics and Space Administration) Ames Research Center. The thermal design of small spacecraft presents challenges as their smaller surface areas translate into power and thermal constraints. The spacecraft is thermally designed to run colder in the Low Earth Orbit space environment, and heated to reach the temperatures required by the science payload. The limited power supply obtained from the solar panels on small surfaces creates a constraint in the power used to heat the payload to required temperatures. The pressurized payload is isolated with low thermally conductance paths from the large ambient temperature changes. The thermal design consists of different optical properties of section surfaces, Multi Layer Insulation (MLI), low thermal conductance materials, flexible heaters and thermal spreaders. The payload temperature is controlled with temperature sensors and flexible heaters. Finite Element Analysis (FEA) and testing were used to aid the thermal design of the spacecraft. Various tests were conducted to verify the thermal design. An infrared imager was used on the electronic boards to find large heat sources and eliminate any possible temperature runaways. The spacecraft was tested in a thermal vacuum chamber to optimize the thermal and power analysis and qualify the thermal design of the spacecraft for the mission.

Diaz-Aguado, Millan F.↗

Design and Laboratory Validation of a Capacitive Sensor for Measuring the Recession of Thin-Layered Ablator

Flight vehicles are typically instrumented with subsurface thermocouples to estimate heat transfer at the surface using inverse analysis procedures. If the vehicle has an ablating heat shield, however, temperature time histories from subsurface thermocouples no longer provide enough information to estimate heat flux at the surface. In this situation, the geometry changes and thermal energy leaves the surface in the form of ablation products. The ablation rate is required to estimate heat transfer to the surface. A new concept for a capacitive sensor has been developed to measure ablator depth using the ablator's dielectric effect on a capacitor's fringe region. Relying on the capacitor's fringe region enables the gage to be flush mounted in the vehicle's permanent structure and not intrude into the ablative heat shield applied over the gage. This sensor's design allows nonintrusive measurement of the thickness of dielectric materials, in particular, the recession rates of low-temperature ablators applied in thin (0.020 to 0.060 in. (0.05 to 0.15 mm)) layers. Twenty capacitive gages with 13 different sensing element geometries were designed, fabricated, and tested. A two-dimensional finite-element analysis was performed on several candidate geometries. Calibration procedures using ablator-simulating shims are described. A one-to-one correspondence between system output and dielectric material thickness was observed out to a thickness of 0.055 in. (1.4 mm) for a material with a permittivity about three times that of air or vacuum. A novel method of monitoring the change in sensor capacitance was developed. This technical memorandum suggests further improvements in gage design and fabrication techniques.

Gregory K Noffz↗

Dynamic Stability of Uncertain Laminated Beams Under Subtangential Loads

Because of the inherent complexity of fiber-reinforced laminated composites, it can be challenging to manufacture composite structures according to their exact design specifications, resulting in unwanted material and geometric uncertainties. In this research, we focus on the deterministic and probabilistic stability analysis of laminated structures subject to subtangential loading, a combination of conservative and nonconservative tangential loads, using the dynamic criterion. Thus a shear-deformable laminated beam element, including warping effects, is derived to study the deterministic and probabilistic response of laminated beams. This twenty-one degrees of freedom element can be used for solving both static and dynamic problems. In the first-order shear deformable model used here we have employed a more accurate method to obtain the transverse shear correction factor. The dynamic version of the principle of virtual work for laminated composites is expressed in its nondimensional form and the element tangent stiffness and mass matrices are obtained using analytical integration The stability is studied by giving the structure a small disturbance about an equilibrium configuration, and observing if the resulting response remains small. In order to study the dynamic behavior by including uncertainties into the problem, three models were developed: Exact Monte Carlo Simulation, Sensitivity Based Monte Carlo Simulation, and Probabilistic FEA. These methods were integrated into the developed finite element analysis. Also, perturbation and sensitivity analysis have been used to study nonconservative problems, as well as to study the stability analysis, using the dynamic criterion.

Goyal, Vijay K.↗

Implementation of J-A Methodology Elastic-Plastic Crack Instability Analysis Capability into the WARP-3D Code

Characterization of the near crack-tip stress/strain fields is the foundation of fracture mechanics. The description of the near tip stress field and the prediction of when fracture occurs is well established for brittle materials that exhibit linear elastic behavior. However, in ductile materials or conditions that violate linear elastic assumptions (Aluminum alloys, Al 2024-T3, Al 2024- T351 etc.), the elastic-plastic crack-tip stress fields are characterized by the Hutchison-Rice-Rosengren (HRR) field. The J-integral is commonly used to characterize amplitude of the HRR field under elastic-plastic conditions. The J-integral has been demonstrated for crack-tip fields that are under high constraint conditions (i.e., small-scale plasticity where the J-dominance is maintained). However, as the external load increases, yielding changes from small- to largescale plasticity and usually a loss of constraint (i.e., reduction in the triaxial stress field along the crack front). The loss of constraint leads to the deviation of the crack-tip stress fields from that given by the HRR field. Hence, the J-dominance will be gradually lost and additional parameter(s) are required to quantify the crack-tip stress fields and predict fracture behavior. The assessment objectives were to: 1) implement a two-parameter (i.e., J-A) fracture criterion into an elastic-plastic three-dimensional (3D) finite element analysis (FEA), 2) validate the implementation by comparison with the A parameter from literature data, 3) conduct material characterization tests to quantify the material behavior and provide fracture data for validation of the J-A fracture criteria, and (4) perform evaluations to establish if the J-A criteria can be used to predict fracture in a ductile metallic material (e.g., aluminum alloys). The A parameter in these criteria is the second parameter in a three-term elastic-plastic asymptotic expansion of the neartip stress behavior. A series of extensive FEAs were performed using WARP3D software package to obtain solutions for the A parameter for different specimen configurations. The methodology needed for the estimation of the A parameter in the asymptotic expansion was developed and implemented using Matlab®. A user material (UMAT) routine was used to model the material stress-strain response using a Ramberg-Osgood power law with a hardening exponent (n) and a material coefficient (alpha). This UMAT routine was successfully implemented in WARP3D software and validated through comparison with the experimental data. Three configurations were extracted from published results: 1) center cracked plate (CCP), 2) single edge-cracked plate (SECP), and 3) double edge-cracked plate (DECP). These configurations and four other configurations (three-hole tension (THT)), three-point bend (3PTB), three-hole compact tension (3PCT), and compact tension (CT)) were analyzed to verify the methodology that was developed and implemented into WARP3D. Solutions of the A parameter were obtained for remote tension loading conditions that started with small-scale yielding and continued into the large-scale plasticity regime. The results indicate that the methodology developed can be used to calculate the elastic-plastic J-A parameters for test specimens with a range of crack geometries, material strain hardening behaviors, and loading conditions. The J-A parameters were implemented as fracture criteria and used to predict the test results. For comparison, other fracture criteria were used to predict the same test results. Major findings include: The A constraint parameter A varies with specimen type and applied load thus accurate determination is crucial in predicting the failure load, and the A parameter is asymptotic as the failure load is approached, making an accurate determination difficult (i.e., small differences in the A parameter can cause large variations in failure load) for materials exhibiting elastic-plastic behavior. The failure predictions from J-A methodology were more accurate than the traditionally used KC and J methods, and have comparable scatter to that observed when using the crack-tip opening angle (CTOA) method. However, the J-A methodology requires considerable effort (expertise level and labor) to implement and to evaluate the A parameter for different specimen types and materials, or to apply this methodology to part-through crack (e.g., 3D problems) structural applications.

Hamm, Kenneth R., Jr.↗

Experimental and Analytical Characterization of Fluted-Core Sandwich Composite Structures

As composite materials continue to gain wider acceptance for aerospace primary structures, experimental characterization and analysis methods must evolve, especially if new or underutilized composite structural concepts are to be successfully implemented. One such potentially underutilized concept is the fluted-core sandwich composite, which consists of facesheets separated by angled web members, and has a number of potential advantages over traditional-core (honeycomb, foam, etc.) sandwich composites. However, fluted-core sandwich composites also exhibit a behaviors not seen with traditional-core sandwich composites. For example, stable local buckling that can lead to material failure can occur in the facesheets between webs or in the webs themselves. Additionally, in the direction perpendicular to the webs, the effective transverse-shear stiffness is largely determined by the bending stiffnesses of the facesheets and the webs and is generally much lower than in the direction parallel to the webs. These local buckling and transverse-shear behaviors can have a large influence on material failure and global buckling, respectively. Thus, it is important to accurately characterize and model these behaviors. In this report, a series of experiments designed to elicit some of these fluted-core sandwich composite structural responses is discussed and the results are compared with results from a previously developed finite element analysis approach; qualitative and quantitative agreement between the experiment and analysis is shown.

Marc R. Schultz↗

On Calculating the Zero-Gravity Surface Figure of a Mirror

An analysis of the classical method of calculating the zero-gravity surface figure of a mirror from surface-figure measurements in the presence of gravity has led to improved understanding of conditions under which the calculations are valid. In this method, one measures the surface figure in two or more gravity- reversed configurations, then calculates the zero-gravity surface figure as the average of the surface figures determined from these measurements. It is now understood that gravity reversal is not, by itself, sufficient to ensure validity of the calculations: It is also necessary to reverse mounting forces, for which purpose one must ensure that mountingfixture/ mirror contacts are located either at the same places or else sufficiently close to the same places in both gravity-reversed configurations. It is usually not practical to locate the contacts at the same places, raising the question of how close is sufficiently close. The criterion for sufficient closeness is embodied in the St. Venant principle, which, in the present context, translates to a requirement that the distance between corresponding gravity-reversed mounting positions be small in comparison to their distances to the optical surface of the mirror. The necessity of reversing mount forces is apparent in the behavior of the equations familiar from finite element analysis (FEA) that govern deformation of the mirror.

Bloemhof, Eric E.↗

Exploration of Failure and Potential Damage Markers in Ti-6Al-4V

To improve prognostics for predicting component life, coupon tests were conducted on Ti-6Al-4V to investigate the development of damage. Various test types were employed (tension, creep, stress relaxation, and fatigue) under a wide range of temperatures and loading rates to engage various amounts of time-dependent and damage behavior. Both stiffness degradation and Poisson’s ratio were monitored to see if they were useful in indicating the initiation and accumulation of damage and subsequently the onset of failure. Stiffness exhibited large decreases with increasing strain when calculated using engineering stress and strain. However, minimal stiffness change was observed when using the true-stress and true- strain values. Detailed microscopy and image analysis were performed to document the actual physical damage that occurred within the various specimens. Two types of damage were observed: shallow surface cracks and internal pores. The severity of damage was greater in the region of localized necking and diminished further away from the neck. The area fraction of pores was higher in the necked region, yet even there it did not exceed more than a few percent and, in most cases, was significantly lower. It was observed that conditions that induced more time-dependent deformation (e.g., slower strain rate testing) resulted in a larger number of pores and larger pore sizes. The pores were always found to occur in the softer β-phase. The severity of observed damage was inconsistent with the large drop in engineering stiffness, but was more aligned with stiffness changes from the true-stress and -strain values. Poisson’s ratio was monitored on every test and showed promise as a damage indicator. It was found to be more sensitive than either the ultimate tensile strength or the inception of tertiary creep. However, calculation of Poisson’s ratio was found to be more involved than expected because of anisotropy and the presence of strain gradients along the specimen at large applied strains. First, the material had a moderate texture and resulted in anisotropy of the strains. Second, the evolution of Poisson’s ratio as a function of axial strain showed strange and unexpected behaviors, which necessitated more extensive analyses of its development. This led to an investigation on the use of small versus large strain approaches and the importance of making consistent measurements; that is, axial and transverse strain measurements over the same volume of material. Additionally, finite element analysis of a tensile sample indicated that non- uniform strain occurred in the sample prior to attaining the maximum strength. This, along with the location at which transverse strain was measured, had a significant influence on the resulting values of Poisson’s ratio. Based on this, lessons learned are offered for future testing.

Ti-6Al-4V↗

A finite element for thermal stress analysis of shells of revolution

A new finite element is described for performing detailed thermal stress analysis of thin orthotropic shells of revolution. The element provides for temperature loadings which may vary over the surface of the shell as well as through the thickness. In a number of sample calculations, results from the present method are compared with analytical solutions as well as with independent numerical analyses. Such calculations are carried out for two cylinders, a conical frustum, a truncated hemisphere, and an annular plate. Generally, the agreement between the present solution and the other solutions is excellent.

Adelman, H. M.↗

Borehole probes: Summary of previous work

Previous research explored the applicability of some devices to lunar borehold testing and, in particular, examined the use of such instruments to assess the strength of soils and rocks. Special attention was given to borehole jacks, which essentially conduct a plate bearing test across opposing wall areas. Analytical work employed finite element analysis, the theory of plasticity, and the theory of elasticity using the complex variable method. A number of solutions were reached for tractable subproblems in the set of problems posed when a jack is used to expand sectors of a borehole. In addition, an investigation with physical models was made to explore the modes of behavior for varying boundary conditions and materials, including both elastic and plastic media.

Goodman, R. E.↗

Deformation mechanics of deep surface flaw cracks

A combined analytical and experimental program was conducted to determine the deformation characteristics of deep surface cracks in Mode I loading. An approximate plane finite element analysis was performed to make a parameter study on the influence of crack depth, crack geometry, and stress level on plastic zones, crack opening displacement, and back surface dimpling in Fe-3Si steel and 2219-T87 aluminum. Surface replication and profiling techniques were used to examine back surface dimple configurations in 2219-T87 aluminum. Interferometry and holography were used to evaluate the potential of various optical techniques to detect small surface dimples on large surface areas.

Francis, P. H.↗

The transverse Poisson's ratio of composites.

An expression is developed that makes possible the prediction of Poisson's ratio for unidirectional composites with reference to any pair of orthogonal axes that are normal to the direction of the reinforcing fibers. This prediction appears to be a reasonable one in that it follows the trends of the finite element analysis and the bounding estimates, and has the correct limiting value for zero fiber content. It can only be expected to apply to composites containing stiff, circular, isotropic fibers bonded to a soft matrix material.

Foye, R. L.↗