Multidisciplinary Sensitivity Analysis and Design Optimization of Flexible Wings Using the Euler Equations on Unstructured Grids

Aeroelasticity is a classical discipline. However, even with recent advancement in computational technology, it still remains a challenging discipline. This is particularly true when aeroelasticity problems are solved in a loosely coupled manner. The advantage of a loosely coupled scheme is that the legacy codes of computational fluid dynamics (CFD) and computational structural mechanics (CSM) can be preserved and used as independent modules in solving the desired aeroelasticity problems. In such a scheme, maintaining proper data transfer (load transfer and deformation tracking between CFD and CSM codes) is crucial in ensuring successful coupled solutions. This can be a challenging task because the interface (the wetted area or the outer mold line) may be discretized differently as required by different levels of computational domains, which leads to mismatch and even gaps. Most published works overcome these drawbacks by transferring the aerodynamic loads and elastic deformation through projection and curve-fitting. The projection is used to find the “host” node, element or Gaussian point from a CSM mesh to the associated CFD mesh, or vice versa. It is then followed by local or global curve-fitting so that the nodal values on the projected surface can be extracted from interpolation. Most of these works can not guarantee a “consistent and conservative” load transfer. Further, they have not adequately demonstrated their availability to support coupled sensitivity analysis. A new remeshing scheme that can guarantee consistent and conservative load transfer and smooth deformation tracking between CFD and CSM is proposed here not only for coupled analysis, but also for coupled sensitivity analysis. The method will introduce an artificial interface structure that is confined with the aerodynamic surface mesh and is supported at the structural surface nodes. This structure is used to redistribute the aerodynamic load as well as the structural deformation. With the help of this artificial structure, the same design parameter that guides the remeshing processes of the structural mesh can be used to guide the remesh processes of the artificial interface structural mesh as well as the aerodynamic interior mesh. This particular feature of the proposed remeshing scheme allows the aerodynamic sensitivity coefficients of a rigid wing to be predetermined and later used for coupled sensitivity analysis that includes only the structural sensitivity code in an iterative routine. A flexible wing with a 3-dimensional Euler flow and a linear finite element model is considered in the present work to demonstrate the proposed scheme for coupled analysis and sensitivity analysis. Preliminary results obtained from the optimization process are presented to substantiate the efficiency of proposed schemes

Predicting Failure Progression and Failure Loads in Composite Open-Hole Tension Coupons

Failure types and failure loads in carbon-epoxy [45n/90n/-45n/0n]ms laminate coupons with central circular holes subjected to tensile load are simulated using progressive failure analysis (PFA) methodology. The progressive failure methodology is implemented using VUMAT subroutine within the ABAQUS(TradeMark)/Explicit nonlinear finite element code. The degradation model adopted in the present PFA methodology uses an instantaneous complete stress reduction (COSTR) approach to simulate damage at a material point when failure occurs. In-plane modeling parameters such as element size and shape are held constant in the finite element models, irrespective of laminate thickness and hole size, to predict failure loads and failure progression. Comparison to published test data indicates that this methodology accurately simulates brittle, pull-out and delamination failure types. The sensitivity of the failure progression and the failure load to analytical loading rates and solvers precision is demonstrated

Test and Analysis Correlation for Sandwich Composite Longitudinal Joint Specimens

The NASA Composite Technology for Exploration (CTE) project is tasked with evaluating methods to analyze and manufacture composite joints for potential use in block upgrades to the Space Launch System (SLS) launch-vehicle structures such as the Payload Attach Fitting (PAF). To perform this task, the CTE project has initiated test and analysis correlation studies for composite joints under various loading conditions. Herein, NASA-developed numerical models are correlated with the experimental results from a series of tension tests. Pretest strain results matched the far-field test data well, but did not capture the nonlinear response in the vicinity of the joint. A refined pretest analytical model was modified to represent progressive failure of the specimens at failure locations observed during the experimental tests. The nonlinear strain response from this progressive failure model predicted the delamination failure load within 15% of the test data, but underpredicted the nonlinearity of the strain response. Further study of composite material models that account for the nonlinear shear response of fabric composites is recommended for the composite joint structures considered in this paper

Influence of Shear Stiffness Degradation on Crack Paths in Uni-Directional Composite Laminates

Influence of shear stiffness degradation in an element, due to damage, on crack paths in uni-directional laminates has been demonstrated. A new shear stiffness degradation approach to improve crack path prediction has been developed and implemented in an ABAQUS/Explicit frame work using VUMAT. Three progressive failure analysis models, built-in ABAQUS (TradeMark), original COmplete STress Reduction (COSTR) and the modified COSTR damage models have been utilized in this study to simulate crack paths in five unidirectional notched laminates, 15deg, 30deg, 45deg, 60deg and 75deg under uniaxial tension load. Results such as crack paths and load vs. edge displacement curves are documented in this report. Modified COSTR damage model shows better accuracy in predicting crack paths in all the uni-directional laminates compared to the ABAQUS (TradeMark) and the original COSTR damage models

Comparison of Damage Path Predictions for Composite Laminates by Explicit and Standard Finite Element Analysis Tools

Splitting, ultimate failure load and the damage path in center notched composite specimens subjected to in-plane tension loading are predicted using progressive failure analysis methodology. A 2-D Hashin-Rotem failure criterion is used in determining intra-laminar fiber and matrix failures. This progressive failure methodology has been implemented in the Abaqus/Explicit and Abaqus/Standard finite element codes through user written subroutines "VUMAT" and "USDFLD" respectively. A 2-D finite element model is used for predicting the intra-laminar damages. Analysis results obtained from the Abaqus/Explicit and Abaqus/Standard code show good agreement with experimental results. The importance of modeling delamination in progressive failure analysis methodology is recognized for future studies. The use of an explicit integration dynamics code for simple specimen geometry and static loading establishes a foundation for future analyses where complex loading and nonlinear dynamic interactions of damage and structure will necessitate it

Testing and Analysis Correlation of Composite Sandwich Longitudinal Bonded Joints for Space Launch Vehicle Structures

The NASA Composite Technology for Exploration (CTE) Project has been developing and demonstrating critical composite technologies with a focus on joints; incorporating materials, design/analysis, manufacturing, and tests that utilize NASA expertise and capabilities. The CTE project has focused on the development of composite longitudinal bonded joint technologies for conical structures such as the SLS Payload Attach Fitting (PAF) due to challenging joint geometries and loads compared to cylindrical jointed structures. The CTE team selected and designed a double-lap composite bonded joint as the most advantageous longitudinal joint to advance for the CTE project. This paper reports on the longitudinal bonded joint sub-element test articles that were fabricated and tested for several loading conditions to test the capability of the bonded joint design. Test and analysis correlation to the sub-element test articles are presented in the paper

Numerical Characterization of a Composite Bonded Wing-Box

The development of composite wing structures has focused on the use of mechanical fasteners to join heavily-loaded areas, while bonded joints have been used only for select locations. The focus of this paper is the examination of the adhesive layer in a generic bonded wing box that represents a "fastenerless" or unitized structure in order to characterize the general behavior and failure mechanisms. A global/local approach was applied to study the response of the adhesive layer using a global shell model and a local shell/solid model. The wing box was analyzed under load to represent a high-g up-bending condition such that the strains in the composite sandwich face sheets are comparable to an expected design allowable. The global/local analysis indicates that at these wing load levels the strains in the adhesive layer are well within the adhesive's elastic region, such that yielding would not be expected in the adhesive layer. The global/local methodology appears to be a promising approach to evaluate the structural integrity of the adhesively bonded structures

Buckling Imperfection Sensitivity of Conical Sandwich Composite Structures for Launch-Vehicles

Structural stability can be an important consideration in the design of large composite shell structures and therefore it is important to understand the buckling response of such structures. It is well known that geometric imperfections can significantly influence the buckling response of such structures by causing the buckling loads to be significantly lower than the theoretical buckling load of a geometrically perfect shell structure. Results are presented of an analytical study on the buckling imperfection sensitivity of large-scale conical sandwich structures for launch vehicles. In particular, representative structures from the Space Launch System launch-vehicle development activities will be considered. The study considered composite sandwich conical structures with multiple sandwich core thicknesses and facesheet layups consisting of tape and fabric composite layups. The results of this analytical study indicate that there is conservatism in the NASA current buckling knockdown factor of 0.33 for conical shell structures. Therefore, it is suggested that the buckling response of composite sandwich cones be further investigated through buckling tests and analytical predictions to potentially revise the buckling design recommendations for conical composite structures

Damage Simulation in Non-Crimp Fabric Composite Plates Subjected to Impact Loads

Progressive failure analysis (PFA) of non-crimp fabric (NCF) composite laminates subjected to low velocity impact loads was performed using the COmplete STress Reduction (COSTR) damage model implemented through VUMAT and UMAT41 user subroutines in the frame works of the commercial finite element programs ABAQUS/Explicit and LS-DYNA, respectively. To validate the model, low velocity experiments were conducted and detailed correlations between the predictions and measurements for both intra-laminar and inter-laminar failures were made. The developed material and damage model predicts the peak impact load and duration very close with the experimental results. Also, the simulation results of delamination damage between the ply interfaces, in-plane matrix damages and fiber damages were all in good agreement with the measurements from the non-destructive evaluation data

Test and Analysis of Sub-Components of Aluminum-Lithium Alloy Cylinders

Integrally machined blade-stiffened panels subjected to an axial compressive load were tested and analyzed to observe the buckling, crippling, and postcrippling response of the panels. The panels were fabricated from aluminum-lithium alloys 2195 and 2050, and both alloys have reduced material properties in the short transverse material direction. The tests were designed to capture a failure mode characterized by the stiffener separating from the panel in the postbuckling range. This failure mode is attributed to the reduced properties in the short transverse direction. Full-field measurements of displacements and strains using three-dimensional digital image correlation systems and local measurements using strain gages were used to capture the deformation of the panel leading up to the failure of the panel for specimens fabricated from 2195. High-speed cameras were used to capture the initiation of the failure. Finite element models were developed using an isotropic strain-hardening material model. Good agreement was observed between the measured and predicted responses for both alloys