Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite element model”

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 865 records · Page 48

Incorporation of a progressive failure analysis method in the CSM testbed software system

Analysis of the postbuckling behavior of composite shell structures pose many difficult and challenging problems in the field of structural mechanics. Current analysis methods perform well for most cases in predicting the postbuckling response of undamaged components. To predict component behavior accurately at higher load levels, the analysis must include the effects of local material failures. The CSM testbed software system is a highly modular structural analysis system currently under development at Langley Research Center. One of the primary goals of the CSM testbed is to provide a software environment for the development of advanced structural analysis methods and modern numerical methods which will exploit advanced computer architecture such as parallel-vector processors. Development of a progressive failure analysis method consists of the design and implementation of a processor which will perform the ply-level progressive failure analysis and the development of a geometrically nonlinear analysis procedure which incorporates the progressive failure processor. Regarding the development of the progressive failure processor, two components are required: failure criteria and a degradation model. For the initial implementation, the failure criteria of Hashin will be used. For a matrix failure which typically indicates the development of transverse matrix cracks, the ply properties will be degraded. Work to date includes the design of the progressive failure analysis processor and initial plans for the controlling geometrically nonlinear analysis procedure. The implementation of the progressive failure analysis has begun. Access to the model database and the Hashin failure criteria are completed. Work is in progress on the input/output operations for the processor related data and the finite element model updating procedures. In total the progressive failure processor is approximately one-third complete.

Arenburg, Robert T.↗

Time simulation of flutter with large stiffness changes

Time simulation of flutter, involving large local structural changes, is formulated with a state-space model that is based on a relatively small number of generalized coordinates. Free-free vibration modes are first calculated for a nominal finite-element model with relatively large fictitious masses located at the area of structural changes. A low-frequency subset of these modes is then transformed into a set of structural modal coordinates with which the entire simulation is performed. These generalized coordinates and the associated oscillatory aerodynamic force coefficient matrices are used to construct an efficient time-domain, state-space model for a basic aeroelastic case. The time simulation can then be performed by simply changing the mass, stiffness, and damping coupling terms when structural changes occur. It is shown that the size of the aeroelastic model required for time simulation with large structural changes at a few apriori known locations is similar to that required for direct analysis of a single structural case. The method is applied to the simulation of an aeroelastic wind-tunnel model. The diverging oscillations are followed by the activation of a tip-ballast decoupling mechanism that stabilizes the system but may cause significant transient overshoots.

Karpel, Mordechay↗

Time simulation of flutter with large stiffness changes

Time simulation of flutter, involving large local structural changes, is formulated with a state-space model that is based on a relatively small number of generalized coordinates. Free-free vibration modes are first calculated for a nominal finite-element model with relatively large fictitious masses located at the area of structural changes. A low-frequency subset of these modes is then transformed into a set of structural modal coordinates with which the entire simulation is performed. These generalized coordinates and the associated oscillatory aerodynamic force coefficient matrices are used to construct an efficient time-domain, state-space model for basic aeroelastic case. The time simulation can then be performed by simply changing the mass, stiffness and damping coupling terms when structural changes occur. It is shown that the size of the aeroelastic model required for time simulation with large structural changes at a few a priori known locations is similar to that required for direct analysis of a single structural case. The method is applied to the simulation of an aeroelastic wind-tunnel model. The diverging oscillations are followed by the activation of a tip-ballast decoupling mechanism that stabilizes the system but may cause significant transient overshoots.

Karpel, M.↗

Development of Reliability Based Life Prediction Methods for Thermal and Environmental Barrier Coatings in Ceramic Matrix Composites

Literature survey related to the EBC/TBC (environmental barrier coating/thermal barrier coating) fife models, failure mechanisms in EBC/TBC and the initial work plan for the proposed EBC/TBC life prediction methods development was developed as well as the finite element model for the thermal/stress analysis of the GRC-developed EBC system was prepared. Technical report for these activities is given in the subsequent sections.

Shah, Ashwin↗

Component Characterization of an eVTOL Reference Model for Crashworthiness Studies

Researchers at the National Aeronautics and Space Administration (NASA) Langley Research Center (LaRC) have conducted a series of structural component and seat level tests to improve finite element model (FEM) characterization of a representative vertical take-off and landing (eVTOL) test article developed by NASA. A full-scale dynamic test was conducted on the representative eVTOL test article in November of 2022. The test article represented a high wing, six passenger eVTOL design concept and is referred to as the lift plus cruise (LPC) test article. The full-scale test identified limitations in the analytical models used to predict aircraft structural response, in particular the composite material models did not effectively capture brittle failure of the structure which were measured during dynamic loading. To better understand the mechanism behind the composite material failure mechanisms observed and to improve the FEM, intact sample specimens of the composite airframe structure were recovered from the test article post-test and used in material characterization testing. In addition, the seat configurations used in the LPC test article were further studied using isolated seat and anthropomorphic test device (ATD) drop tower testing. Dynamic compression tests and three-point bend tests, conducted at varied impact speeds, were performed on the recovered frame section specimens. Additional testing was conducted to characterize the material properties of the forming foam, which remained in the frames after fabrication. These tests were used to improve characterization of the damage and failure parameters of the composite material model used in the FE model of the LPC test article. Seat level tests were conducted on the seats used in the LPC test article using acceleration pulses inclusive of current general aviation and rotorcraft certification load levels as well as conditions representative of those measured at the seat base during the LPC test. The structural material models and seat environment models of the LPC test article FEM were calibrated using the generated component test data. The updates made to these models were then integrated into the LPC FEM and simulated in the full-scale test condition. Results demonstrated the effectiveness of component testing to improve predictive capability of composite aerospace structural models within the crash and dynamic loading environments. Demonstration of the LPC FEM response across an accumulation of coupon, component, seat environment, and full-scale test levels provides confidence in the predictive capability of this model for future use in the study of occupant safety within eVTOL relevant crash environments.

Craswhorthiness↗

Component Characterization of an eVTOL Reference Model for Crashworthiness Studies

Researchers at the National Aeronautics and Space Administration (NASA) Langley Research Center (LaRC) have conducted a series of structural component and seat level tests to improve finite element model (FEM) characterization of a representative vertical take-off and landing (eVTOL) test article developed by NASA. A full-scale dynamic test was conducted on the representative eVTOL test article in November of 2022. The test article represented a high wing, six passenger eVTOL design concept and is referred to as the lift plus cruise (LPC) test article. The full-scale test identified limitations in the analytical models used to predict aircraft structural response, in particular the composite material models did not effectively capture brittle failure of the structure which were measured during dynamic loading. To better understand the mechanism behind the composite material failure mechanisms observed and to improve the FEM, intact sample specimens of the composite airframe structure were recovered from the test article post-test and used in material characterization testing. In addition, the seat configurations used in the LPC test article were further studied using isolated seat and anthropomorphic test device (ATD) drop tower testing. Dynamic compression tests and three-point bend tests, conducted at varied impact speeds, were performed on the recovered frame section specimens. Additional testing was conducted to characterize the material properties of the forming foam, which remained in the frames after fabrication. These tests were used to improve characterization of the damage and failure parameters of the composite material model used in the FE model of the LPC test article. Seat level tests were conducted on the seats used in the LPC test article using acceleration pulses inclusive of current general aviation and rotorcraft certification load levels as well as conditions representative of those measured at the seat base during the LPC test. The structural material models and seat environment models of the LPC test article FEM were calibrated using the generated component test data. The updates made to these models were then integrated into the LPC FEM and simulated in the full-scale test condition. Results demonstrated the effectiveness of component testing to improve predictive capability of composite aerospace structural models within the crash and dynamic loading environments. Demonstration of the LPC FEM response across an accumulation of coupon, component, seat environment, and full-scale test levels provides confidence in the predictive capability of this model for future use in the study of occupant safety within eVTOL relevant crash environments.

Craswhorthiness↗

Simulation of Delamination-Migration and Core Crushing in a CFRP Sandwich Structure

Following the onset of damage caused by an impact load on a composite laminate structure, delaminations often form propagating outwards from the point of impact and in some cases can migrate via matrix cracks between plies as they grow. The goal of the present study is to develop an accurate finite element modeling technique for simulation of the delamination-migration phenomena in laminate impact damage processes. An experiment was devised where, under a quasi-static indentation load, an embedded delamination in the facesheet of a laminate sandwich specimen migrates via a transverse matrix crack and then continues to grow on a new ply interface. The quasistatic nature of the indentation results in structural behavior equivalent to that seen in low-velocity impact and also allows for highly detailed real time damage characterization. Several finite element damage simulation methods were investigated. Comparing the experimental results with those of the different models reveals certain modeling features that are important to include in a numerical simulation of delamination-migration and some that may be neglected.

McElroy, M.↗

Fiber shape effects on metal matrix composite behavior

The effects of different fiber shapes on the behavior of a SiC/Ti-15 metal matrix composite is computationally simulated. A three-dimensional finite element model consisting of a group of nine unidirectional fibers is used in the analysis. The model is employed to represent five different fiber shapes: a circle, an ellipse, a kidney, and two different cross shapes. The distribution of microstresses and the composite material properties, such as moduli, coefficients of thermal expansion, and Poisson's ratios, are obtained from the finite element analysis for the various fiber shapes. Comparisons of these results are used to determine the sensitivity of the composite behavior to the different fiber shapes and assess their potential benefits. No clear benefits result from different fiber shapes though there are some increases/decreases in isolated properties.

Brown, H. C.↗

Design/Analysis of Metal/Composite Bonded Joints for Survivability at Cryogenic Temperatures

A major design and analysis challenge for the JWST ISM structure is the metal/composite bonded joints that will be required to survive down to an operational ultra-low temperature of 30K (-405 F). The initial and current baseline design for the plug-type joint consists of a titanium thin walled fitting (1-3mm thick) bonded to the interior surface of an M555/954-6 composite truss square tube with an axially stiff biased lay-up. Metallic fittings are required at various nodes of the truss structure to accommodate instrument and lift-point bolted interfaces. Analytical experience and design work done on metal/composite bonded joints at temperatures below liquid nitrogen are limited and important analysis tools, material properties, and failure criteria for composites at cryogenic temperatures are virtually nonexistent. Increasing the challenge is the difficulty in testing for these required tools and parameters at 30K. A preliminary finite element analysis shows that failure due to CTE mismatch between the biased composite and titanium or aluminum is likely. Failure is less likely with Invar, however an initial mass estimate of Invar fittings demonstrates that Invar is not an automatic alternative. In order to gain confidence in analyzing and designing the ISM joints, a comprehensive joint development testing program has been planned and is currently running. The test program is designed for the correlation of the analysis methodology, including tuning finite element model parameters, and developing a composite failure criterion for the effect of multi-axial composite stresses on the strength of a bonded joint at 30K. The testing program will also consider stress mitigation using compliant composite layers and potential strength degradation due to multiple thermal cycles. Not only will the finite element analysis be correlated to the test data, but the FEA will be used to guide the design of the test. The first phase of the test program has been completed and the preliminary analysis has been revisited based on the test data In this work, we present an overview of the test plan, results today, and resulting design improvements.

Bartoszyk, Andrew E.↗

Pressure distribution under flexible polishing tools. II - Cylindrical (conical) optics

A previously developed eigenvalue model is extended to determine polishing pressure distribution by rectangular tools with unequal stiffness in two directions on cylindrical optics. Tool misfit is divided into two simplified one-dimensional problems and one simplified two-dimensional problem. Tools with nonuniform cross-sections are treated with a new one-dimensional eigenvalue algorithm, permitting evaluation of tool designs where the edge is more flexible than the interior. This maintains edge pressure variations within acceptable parameters. Finite element modeling is employed to resolve upper bounds, which handle pressure changes in the two-dimensional misfit element. Paraboloids and hyperboloids from the NASA AXAF system are treated with the AXAFPOD software for this method, and are verified with NASTRAN finite element analyses. The maximum deviation from the one-dimensional azimuthal pressure variation is predicted to be 10 percent and 20 percent for paraboloids and hyperboloids, respectively.

Mehta, Pravin K.↗

Soft Soil Impact Testing and Simulation of Aerospace Structures

In June 2007, a 38-ft/s vertical drop test of a 5-ft-diameter, 5-ft-long composite fuselage section that was retrofitted with a novel composite honeycomb Deployable Energy Absorber (DEA) was conducted onto unpacked sand. This test was one of a series of tests to evaluate the multi-terrain capabilities of the DEA and to generate test data for model validation. During the test, the DEA crushed approximately 6-in. and left craters in the sand of depths ranging from 7.5- to 9-in. A finite element model of the fuselage section with DEA was developed for execution in LS-DYNA, a commercial nonlinear explicit transient dynamic code. Pre-test predictions were generated in which the sand was represented initially as a crushable foam material MAT_CRUSHABLE_FOAM (Mat 63). Following the drop test, a series of hemispherical penetrometer tests were conducted to assist in soil characterization. The penetrometer weighed 20-lb and was instrumented with a tri-axial accelerometer. Drop tests were performed at 16-ft/s and crater depths were measured. The penetrometer drop tests were simulated as a means for developing a more representative soil model based on a soil and foam material definition MAT_SOIL_AND FOAM (Mat 5) in LS-DYNA. The model of the fuselage with DEA was reexecuted using the updated soil model and test-analysis correlations are presented.

Fasanella, Edwin L.↗

Harnessing Collaborative Learning Automata to Guide Multi-objective Optimization based Inverse Analysis for Structural Damage Identification

Structural damage identification based on physical models is often transformed into an optimization problem that minimizes the difference between measurement information of structure being monitored and the model prediction in the parametric space. However, the objective function in this context often exhibits multimodality, involving high-dimensional variables due to the reliance on finite element models for damage identification. These features pose challenges to optimization algorithms, where entrapment in local solutions can lead to false positives and false negatives in damage identification. In this research, we propose a reinforcement learning based multi-swarm optimizer to tackle such challenges in pursuit of a small yet diverse solution set that can capture the true damage scenario as one of the solutions. The proposed method leverages the flexibility of the particle swarm optimizer and incorporates novel strategies of metaheuristics to realize targeted improvement. To enable the particle swarm to adaptively select the appropriate search strategy based on the current environment, we adopt the learning automata technique, which sidesteps the need for reward strategy selection that is usually ad hoc at each step of the search. The integration harnesses the automatic learning and self-adaptation capabilities of learning automata, enabling the particles to navigate based on environmental signals. This leads to accumulated probabilities tied to advantageous movements, fostering an adaptive exploration of particles in the search space. The proposed approach is first validated through implementing into benchmark test cases with comparisons. It is then applied to structural damage identification with piezoelectric admittance experimental signals. `The results highlight the capability of the algorithm to identify a small solution set with high accuracy to match the actual damage scenario.

Yang Zhang↗

Shunted Piezoelectric Vibration Damping Analysis Including Centrifugal Loading Effects

Excessive vibration of turbomachinery blades causes high cycle fatigue problems which require damping treatments to mitigate vibration levels. One method is the use of piezoelectric materials as passive or active dampers. Based on the technical challenges and requirements learned from previous turbomachinery rotor blades research, an effort has been made to investigate the effectiveness of a shunted piezoelectric for the turbomachinery rotor blades vibration control, specifically for a condition with centrifugal rotation. While ample research has been performed on the use of a piezoelectric material with electric circuits to attempt to control the structural vibration damping, very little study has been done regarding rotational effects. The present study attempts to fill this void. Specifically, the objectives of this study are: (a) to create and analyze finite element models for harmonic forced response vibration analysis coupled with shunted piezoelectric circuits for engine blade operational conditions, (b) to validate the experimental test approaches with numerical results and vice versa, and (c) to establish a numerical modeling capability for vibration control using shunted piezoelectric circuits under rotation. Study has focused on a resonant damping control using shunted piezoelectric patches on plate specimens. Tests and analyses were performed for both non-spinning and spinning conditions. The finite element (FE) shunted piezoelectric circuit damping simulations were performed using the ANSYS Multiphysics code for the resistive and inductive circuit piezoelectric simulations of both conditions. The FE results showed a good correlation with experimental test results. Tests and analyses of shunted piezoelectric damping control, demonstrating with plate specimens, show a great potential to reduce blade vibrations under centrifugal loading.

Min, James B.↗

Adaptation of two-dimensional transonic analysis code, TSFOIL, for use in modeling adaptive wall test sections

The effects of the elastic deformation of the wind tunnel wall held to a streamline-like shape were simulated numerically. The wall itself is simulated by a finite element model and is allowed to deform under the pressure loading developing in the wind tunnel with an airfoil model present. A modified version of the transonic analysis program, TSFOIL, is then used to determine the resulting flow field with the effects of the deformed wall included, the shape of the wall and the flow field. Once a result from a particular operating condition was obtained, the pressure distribution on the airfoil in the wind tunnel model is compared with solutions generated by TSFOIL in the free air mode.

Webb, J.↗

Polar Plate Theory for Orthogonal Anisotropy

Laminated fiber-reinforced (or filamentary) composites are used today for their high strength-to-weight and stiffness-to-weight ratios. However, because of the anisotropic behavior of composites, determining the response on a macroscopic scale is challenging. This is particularly evident in the evaluation of the governing differential equations of a circular disk with the fibers of the lamina oriented with rectilinear orthogonality. This includes any situation involving a composite plate of circular geometry in which out-of-plane displacements due to load are desired, such as fastener pull through loading of a composite plate. Current analysis techniques use numerical methods with rectilinear coordinate systems to solve problems with circular geometry. These analyses over predict plate stiffness by 20% and underpredict failure by 70%. Consequently, there is a need to transform classical composite plate theory to a polar coordinate system. In order to better analyze structures with circular geometries the classical composite plate equations are transformed into the plate equations for a rectilinearly anisotropic composite in polar coordinates. A composite plate is typically a laminate of fibers in rectilinear directions. Subsequent to the lay-tip the necessary geometry is cut out of a rectangular plate. In a similar manner, the derivation of the plate equation starts with the fundamental definitions of strain, displacement and curvature and incorporates the material property angular dependence into the equilibrium equations for a differential polar element. In the transformed state, the stiffness coefficients are no longer constant, adding to the complexity of the governing differential equations. This paper discusses the new derivation and evaluation of the plate equations for a circular composite disk with orthogonal rectilinear anisotropy. The resultant new three partial differential equations, which describe the circular anisotropic plate, can be used to evaluate out-of-plane displacements for given load conditions in design of composite Structures. Without the formation of these unique plate equations, design of such structures is not a precise engineering accomplishment due to the lack of precise design tools. Such structures are over designed to compensate and costly tests need to be performed after production to validate safety. Though impossible to present the mathematical solution within the limits of this paper, the application of the newly derived plate equation in its orthotropic form for a balanced symmetric laminate compared favorably with test results and surpassed the results of standard numerical methods. Within the limitations of the problem foundation the predictive model provides a mechanism for establishing out-of-plane deflection levels for a circular composite plate. Future work in this area should include incorporating the circular composite plate equations in Finite Element Models. The new anisotropic polar plate equations call be utilized as a design tool so that the design more accurately meets the requirements, thus taking full advantage of the weight savings of composites. These accomplishments have not previously been presented to the engineering community nor ever proposed to the IAF. Should the proposed paper be accepted for presentation, attendance is assured.

Bailey, Michelle D.↗

Static shape determination and control for a large space antenna

An integral operator approach is used to derive solutions to static shape determination and control problems associated with large space structures. Problem assumptions include a linear self-adjoint system model, observations and control forces at discrete points, and quadratic performance criteria for the comparison of estimates or control forces. Results are illustrated by simulations with a finite element model of a large space antenna. Modal expansions for terms in the solution algorithms are presented, using modes from the static or associated dynamic model. These expansions provide approximate solutions in the event that a closed form analytical solution to the system boundary value problem is not available.

Weeks, C.↗

Influence of local fiber undulations on the global buckling behavior of filament-wound cylinders

The predicted compressive stiffness and buckling strength of filament-wound cylinders using classical lamination theory is significantly higher than those observed experimentally. This discrepancy is partially influenced by the variation of mechanical properties in the region of fiber undulations. These regions are localized geometric defects intrinsic to the filament-winding, weaving, and braiding processes. In the present work, the average mechanical properties of the fiber undulation region are quantified using modified models of woven-fabric composites to account for the 3-dimensional effects. The mechanical properties thus determined can be incorporated as local element properties into global finite-element models. Preliminary results from large-displacement analyses of filament-wound cylinders are relatively more accurate when fiber undulations are accounted for.

Jensen, D. W.↗

A software tool to design thermal barrier coatings

This paper summarizes work completed for a NASA Phase 1 SBIR program which demonstrated the feasibility of developing a software tool to aid in the design of thermal barrier coating (TBC) systems. Toward this goal, three tasks were undertaken and completed. Task 1 involved the development of a database containing the pertinent thermal and mechanical property data for the top coat, bond coat and substrate materials that comprise a TBC system. Task 2 involved the development of an automated set-up program for generating two dimensional (2D) finite element models of TBC systems. Most importantly, Task 3 involved the generation of a rule base to aid in the design of a TBC system. These rules were based on a factorial design of experiments involving FEM results, and were generated using a Yates analysis. A previous study has indicated the suitability and benefit of applying finite element analysis to perform computer based experiments to decrease but not eliminate physical experiments on TBC's. This program proved feasibility by expanding on these findings by developing a larger knowledge base and developing a procedure to extract rules to aid in TBC design.

Petrus, G.↗