Engineering topics

Kapania, Rakesh K.

Publications and source records attributed to Kapania, Rakesh K..

At least 19 records

NASA NTRS2019

Interference Drag Associated with Engine Locations for Multidisciplinary Design Optimization

This research aims to quantify the interference drag for various engine locations on a traditional “tube and wing,” 150-passenger commercial aircraft flying at 35,000 ft and Mach 0.8. Engine locations are varied in the chord wise, span wise, and vertical directions near the wing, both under and above the wing, as well as along the fuselage. Euler simulations are performed with representative powered modern engines. The results are intended to supplement empirical drag estimates suitable for multidisciplinary design environments. Large interference drag increases, as compared to the isolated air frame and engine geometry, are found to occur when the engine is placed directly above or below the wing. Interference effects are significantly reduced, and in some instances result in benefits compared to the isolated bodies, when the engines are placed fore or aft of the wing. Interference drag increases are partially explained by flow channels leading to choked flow and shock interactions between bodies.

Blaesser, Nathaniel J.↗

NASA NTRS2018

Multidisciplinary Design Optimization and Cruise Mach Number Study of Truss-Braced Wing Aircraft

The Subsonic Ultra Green Aircraft Research (SUGAR) Phase III was led by Dr. Rakesh K. Kapania and Dr. Joseph A. Schetz at the Multidisciplinary Analysis and Design Center for Advanced Vehicles, Department of Aerospace and Ocean Engineering, Virginia Tech, Blacksburg VA. The research was performed from December 2014 to December 2015. Three major areas were investigated: Multidisciplinary Design Optimization (MDO) studies of truss braced wing (TBW) and strut braced wing (SBW) vehicles at cruise Mach numbers of 0.7 and 0.8 for a flight mission similar to current market single aisle configurations. The performance and the characteristics of the optimized vehicles were compared to the SUGAR Phase II TBW vehicle. These results were obtained without applying any of the extended transonic aerodynamic and aeroelastic tools that will be discussed later. It was found that the cruise Mach number has a large effect on the best “truss” configuration. At Mach 0.7, an SBW has a better fuel consumption and better take-off gross weight (TOGW). However, at Mach 0.8, the TBW is superior because the jury strut aids in satisfying the flutter constraint; Two-dimensional, steady, transonic aerodynamic analysis of the Boeing Airfoil J (BACJ) airfoil was performed for a range of thickness ratios, Mach numbers and lift coefficients. Reynolds-averaged Navier-Stokes (RANS) equations were solved to obtain the lift-curve slope, wave drag coefficient, the location of the center of pressure and to predict the separation at the trailing edge, which may lead to buffeting. One of the goals was to develop a database of lift-curve slope and the location of center of pressure, which could be used in a transonic aeroelastic analysis. Another goal was to compare the wave drag coefficients to those predicted by Lock’s fourth-power law and also to compare the transonic effects obtained from RANS simulations to those predicted by the Korn equations. A third goal was to develop a buffet boundary, which can be integrated into the MDO framework to prevent the optimized designs from probable buffeting; A state-space transonic aeroelastic analysis tool was developed, which can incorporate the nonlinear transonic effects in the unsteady aerodynamics but is yet computationally cheap when used within the VT MDO framework. The aeroelastic analysis uses Leishman- Beddoes (LB) indicial functions, which generated a state-space representation of the aeroelastic system. The indicial functions allow the incorporation of data for steady lift-curve slope and location of the center of pressure. Thus, the steady transonic effects are included, and the unsteady aerodynamic responses are a linearization about the steady results. The aeroelastic approach discretizes the wing into numerous strips, which results in a large eigenvalue problem as each strip has eight augmented aerodynamic states as per the LB theory. Thus, to reduce the computation expense, a reduced order model (ROM) was developed. The approach was validated using a few examples.

Kapania, Rakesh K.↗

NASA NTRS2017

SpaRibs Geometry Parameterization for Wings with Multiple Sections using Single Design

The SpaRibs topology of an aircraft wing has a significant effect on its structural behavior and stability as well as the flutter performance. The development of additive manufacturing techniques like Electron Beam Free Form Fabrication (EBF3) has made it feasible to manufacture aircraft wings with curvilinear spars, ribs (SpaRibs) and stiffeners. In this article a new global-local optimization framework for wing with multiple sections using curvilinear SpaRibs is described. A single design space is used to parameterize the SpaRibs geometry. This method has been implemented using MSC-PATRAN to create a broad range of SpaRibs topologies using limited number of parameters. It ensures C0 and C1 continuities in SpaRibs geometry at the junction of two wing sections with airfoil thickness gradient discontinuity as well as mesh continuity between all structural components. This method is advantageous in complex multi-disciplinary optimization due to its potential to reduce the number of design variables. For the global-local optimization the local panels are generated by an algorithm which is totally based on a set algebra on the connectivity matrix data. The great advantage of this method is that it is completely independent of the coordinates of the nodes of the finite element model. It is also independent of the order in which the elements are distributed in the FEM. The code is verified by optimizing of the CRM Baseline model at trim condition at Mach number equal to 0.85 for five different angle of attack (-2deg, 0deg,2deg,4deg and 6deg). The final weight of the wing is 19,090.61 lb. This value is comparable to that obtained by Qiang et al. 6 (19,269 lb).

wing sections↗

NASA NTRS2015

Novel Control Effectors for Truss Braced Wing

At cruise flight conditions very high aspect ratio/low sweep truss braced wings (TBW) may be subject to design requirements that distinguish them from more highly swept cantilevered wings. High aspect ratio, short chord length and relative thinness of the airfoil sections all contribute to relatively low wing torsional stiffness. This may lead to aeroelastic issues such as aileron reversal and low flutter margins. In order to counteract these issues, high aspect ratio/low sweep wings may need to carry additional high speed control effectors to operate when outboard ailerons are in reversal and/or must carry additional structural weight to enhance torsional stiffness. The novel control effector evaluated in this study is a variable sweep raked wing tip with an aileron control surface. Forward sweep of the tip allows the aileron to align closely with the torsional axis of the wing and operate in a conventional fashion. Aft sweep of the tip creates a large moment arm from the aileron to the wing torsional axis greatly enhancing aileron reversal. The novelty comes from using this enhanced and controllable aileron reversal effect to provide roll control authority by acting as a servo tab and providing roll control through intentional twist of the wing. In this case the reduced torsional stiffness of the wing becomes an advantage to be exploited. The study results show that the novel control effector concept does provide roll control as described, but only for a restricted class of TBW aircraft configurations. For the configuration studied (long range, dual aisle, Mach 0.85 cruise) the novel control effector provides significant benefits including up to 12% reduction in fuel burn.

White, Edward V.↗

NASA NTRS2014

Artificial Intelligence Based Control Power Optimization on Tailless Aircraft

Traditional methods of control allocation optimization have shown difficulties in exploiting the full potential of controlling large arrays of control devices on innovative air vehicles. Artificial neutral networks are inspired by biological nervous systems and neurocomputing has successfully been applied to a variety of complex optimization problems. This project investigates the potential of applying neurocomputing to the control allocation optimization problem of Hybrid Wing Body (HWB) aircraft concepts to minimize control power, hinge moments, and actuator forces, while keeping system weights within acceptable limits. The main objective of this project is to develop a proof-of-concept process suitable to demonstrate the potential of using neurocomputing for optimizing actuation power for aircraft featuring multiple independently actuated control surfaces. A Nastran aeroservoelastic finite element model is used to generate a learning database of hinge moment and actuation power characteristics for an array of flight conditions and control surface deflections. An artificial neural network incorporating a genetic algorithm then uses this training data to perform control allocation optimization for the investigated aircraft configuration. The phase I project showed that optimization results for the sum of required hinge moments are improved by more than 12% over the best Nastran solution by using the neural network optimization process.

Gern, Frank↗

NASA NTRS2014

Supersonic Wing Optimization Using SpaRibs

This research investigates the advantages of using curvilinear spars and ribs, termed SpaRibs, to design a supersonic aircraft wing-box in comparison to the use of classic design concepts that employ straight spars and ribs. The objective is to achieve a more efficient load-bearing mechanism and to passively control the deformation of the structure under the flight loads. Moreover, the use of SpaRibs broadens the design space and allows for natural frequencies and natural mode shape tailoring. The SpaRibs concept is implemented in a new optimization MATLAB-based framework referred to as EBF3SSWingOpt. This optimization scheme performs both the sizing and the shaping of the internal structural elements, connecting the optimizer with the analysis software. The shape of the SpaRibs is parametrically defined using the so called Linked Shape method. Each set of SpaRibs is placed in a one by one square domain of the natural space. The set of curves is subsequently transformed in the physical space for creating the wing structure geometry layout. The shape of each curve of each set is unique; however, mathematical relations link the curvature in an effort to reduce the number of design variables. The internal structure of a High Speed Commercial Transport aircraft concept developed by Boeing is optimized subjected to stress, subsonic flutter and supersonic flutter constraints. The results show that the use of the SpaRibs allows for the reduction of the aircraft's primary structure weight without violating the constraints. A weight reduction of about 15 percent is observed.

Locatelli, David↗

NASA NTRS2012

Design, Optimization, and Evaluation of A1-2139 Compression Panel with Integral T-Stiffeners

A T-stiffened panel was designed and optimized for minimum mass subjected to constraints on buckling load, yielding, and crippling or local stiffener failure using a new analysis and design tool named EBF3PanelOpt. The panel was designed for a compression loading configuration, a realistic load case for a typical aircraft skin-stiffened panel. The panel was integrally machined from 2139 aluminum alloy plate and was tested in compression. The panel was loaded beyond buckling and strains and out-of-plane displacements were extracted from 36 strain gages and one linear variable displacement transducer. A digital photogrammetric system was used to obtain full field displacements and strains on the smooth (unstiffened) side of the panel. The experimental data were compared with the strains and out-of-plane deflections from a high-fidelity nonlinear finite element analysis.

Mulani, Sameer B.↗

NASA NTRS2012

Compression After Impact on Honeycomb Core Sandwich Panels With Thin Facesheets: Experiments - Part 1

McQuigg, Thomas D.↗

NASA NTRS2012

Compression After Impact on Honeycomb Core Sandwich Panels with Thin Facesheets, Part 2: Analysis

A two part research study has been completed on the topic of compression after impact (CAI) of thin facesheet honeycomb core sandwich panels. The research has focused on both experiments and analysis in an effort to establish and validate a new understanding of the damage tolerance of these materials. Part 2, the subject of the current paper, is focused on the analysis, which corresponds to the CAI testings described in Part 1. Of interest, are sandwich panels, with aerospace applications, which consist of very thin, woven S2-fiberglass (with MTM45-1 epoxy) facesheets adhered to a Nomex honeycomb core. Two sets of materials, which were identical with the exception of the density of the honeycomb core, were tested in Part 1. The results highlighted the need for analysis methods which taken into account multiple failure modes. A finite element model (FEM) is developed here, in Part 2. A commercial implementation of the Multicontinuum Failure Theory (MCT) for progressive failure analysis (PFA) in composite laminates, Helius:MCT, is included in this model. The inclusion of PFA in the present model provided a new, unique ability to account for multiple failure modes. In addition, significant impact damage detail is included in the model. A sensitivity study, used to assess the effect of each damage parameter on overall analysis results, is included in an appendix. Analysis results are compared to the experimental results for each of the 32 CAI sandwich panel specimens tested to failure. The failure of each specimen is predicted using the high-fidelity, physicsbased analysis model developed here, and the results highlight key improvements in the understanding of honeycomb core sandwich panel CAI failure. Finally, a parametric study highlights the strength benefits compared to mass penalty for various core densities.

Mcquigg, Thomas D.↗

NASA NTRS2011

Design, Optimization, and Evaluation of Integrally-Stiffened Al-2139 Panel with Curved Stiffeners

A curvilinear stiffened panel was designed, manufactured, and tested in the Combined Load Test Fixture at NASA Langley Research Center. The panel is representative of a large wing engine pylon rib and was optimized for minimum mass subjected to three combined load cases. The optimization included constraints on web buckling, material yielding, crippling or local stiffener failure, and damage tolerance using a new analysis tool named EBF3PanelOpt. Testing was performed for the critical combined compression-shear loading configuration. The panel was loaded beyond initial buckling, and strains and out-of-plane displacements were extracted from a total of 20 strain gages and 6 linear variable displacement transducers. The VIC-3D system was utilized to obtain full field displacements/strains in the stiffened side of the panel. The experimental data were compared with the strains and out-of-plane deflections from a high fidelity nonlinear finite element analysis. The experimental data were also compared with linear elastic finite element results of the panel/test-fixture assembly. Overall, the panel buckled very near to the predicted load in the web regions.

Havens, David↗

NASA NTRS2011

Design, Optimization and Evaluation of Integrally Stiffened Al 7050 Panel with Curved Stiffeners

A curvilinear stiffened panel was designed, manufactured, and tested in the Combined Load Test Fixture at NASA Langley Research Center. The panel was optimized for minimum mass subjected to constraints on buckling load, yielding, and crippling or local stiffener failure using a new analysis tool named EBF3PanelOpt. The panel was designed for a combined compression-shear loading configuration that is a realistic load case for a typical aircraft wing panel. The panel was loaded beyond buckling and strains and out-of-plane displacements were measured. The experimental data were compared with the strains and out-of-plane deflections from a high fidelity nonlinear finite element analysis and linear elastic finite element analysis of the panel/test-fixture assembly. The numerical results indicated that the panel buckled at the linearly elastic buckling eigenvalue predicted for the panel/test-fixture assembly. The experimental strains prior to buckling compared well with both the linear and nonlinear finite element model.

Slemp, Wesley C. H.↗

NASA NTRS2011

Response of Honeycomb Core Sandwich Panel with Minimum Gage GFRP Face-Sheets to Compression Loading After Impact

A compression after impact study has been conducted to determine the residual strength of three sandwich panel constructions with two types of thin glass fiber reinforced polymer face-sheets and two hexagonal honeycomb Nomex core densities. Impact testing is conducted to first determine the characteristics of damage resulting from various impact energy levels. Two modes of failure are found during compression after impact tests with the density of the core precipitating the failure mode present for a given specimen. A finite element analysis is presented for prediction of the residual compressive strength of the impacted specimens. The analysis includes progressive damage modeling in the face-sheets. Preliminary analysis results were similar to the experimental results; however, a higher fidelity core material model is expected to improve the correlation.

McQuigg, Thomas D.↗

NASA NTRS2010

Interlaminar Stresses by Refined Beam Theories and the Sinc Method Based on Interpolation of Highest Derivative

Computation of interlaminar stresses from the higher-order shear and normal deformable beam theory and the refined zigzag theory was performed using the Sinc method based on Interpolation of Highest Derivative. The Sinc method based on Interpolation of Highest Derivative was proposed as an efficient method for determining through-the-thickness variations of interlaminar stresses from one- and two-dimensional analysis by integration of the equilibrium equations of three-dimensional elasticity. However, the use of traditional equivalent single layer theories often results in inaccuracies near the boundaries and when the lamina have extremely large differences in material properties. Interlaminar stresses in symmetric cross-ply laminated beams were obtained by solving the higher-order shear and normal deformable beam theory and the refined zigzag theory with the Sinc method based on Interpolation of Highest Derivative. Interlaminar stresses and bending stresses from the present approach were compared with a detailed finite element solution obtained by ABAQUS/Standard. The results illustrate the ease with which the Sinc method based on Interpolation of Highest Derivative can be used to obtain the through-the-thickness distributions of interlaminar stresses from the beam theories. Moreover, the results indicate that the refined zigzag theory is a substantial improvement over the Timoshenko beam theory due to the piecewise continuous displacement field which more accurately represents interlaminar discontinuities in the strain field. The higher-order shear and normal deformable beam theory more accurately captures the interlaminar stresses at the ends of the beam because it allows transverse normal strain. However, the continuous nature of the displacement field requires a large number of monomial terms before the interlaminar stresses are computed as accurately as the refined zigzag theory.

Slemp, Wesley C. H.↗

NASA NTRS2003

Control Law Design in a Computational Aeroelasticity Environment

A methodology for designing active control laws in a computational aeroelasticity environment is given. The methodology involves employing a systems identification technique to develop an explicit state-space model for control law design from the output of a computational aeroelasticity code. The particular computational aeroelasticity code employed in this paper solves the transonic small disturbance aerodynamic equation using a time-accurate, finite-difference scheme. Linear structural dynamics equations are integrated simultaneously with the computational fluid dynamics equations to determine the time responses of the structure. These structural responses are employed as the input to a modern systems identification technique that determines the Markov parameters of an "equivalent linear system". The Eigensystem Realization Algorithm is then employed to develop an explicit state-space model of the equivalent linear system. The Linear Quadratic Guassian control law design technique is employed to design a control law. The computational aeroelasticity code is modified to accept control laws and perform closed-loop simulations. Flutter control of a rectangular wing model is chosen to demonstrate the methodology. Various cases are used to illustrate the usefulness of the methodology as the nonlinearity of the aeroelastic system is increased through increased angle-of-attack changes.

Newsom, Jerry R.↗

NASA NTRS2002

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.↗

NASA NTRS2001

Free Vibration of Uncertain Unsymmetrically Laminated Beams

Monte Carlo Simulation and Stochastic FEA are used to predict randomness in the free vibration response of thin unsymmetrically laminated beams. For the present study, it is assumed that randomness in the response is only caused by uncertainties in the ply orientations. The ply orientations may become random or uncertain during the manufacturing process. A new 16-dof beam element, based on the first-order shear deformation beam theory, is used to study the stochastic nature of the natural frequencies. Using variational principles, the element stiffness matrix and mass matrix are obtained through analytical integration. Using a random sequence a large data set is generated, containing possible random ply-orientations. This data is assumed to be symmetric. The stochastic-based finite element model for free vibrations predicts the relation between the randomness in fundamental natural frequencies and the randomness in ply-orientation. The sensitivity derivatives are calculated numerically through an exact formulation. The squared fundamental natural frequencies are expressed in terms of deterministic and probabilistic quantities, allowing to determine how sensitive they are to variations in ply angles. The predicted mean-valued fundamental natural frequency squared and the variance of the present model are in good agreement with Monte Carlo Simulation. Results, also, show that variations between plus or minus 5 degrees in ply-angles can affect free vibration response of unsymmetrically and symmetrically laminated beams.

Kapania, Rakesh K.↗

NASA NTRS2000

Efficient Simulation of Wing Modal Response: Application of 2nd Order Shape Sensitivities and Neural Networks

At the preliminary design stage of a wing structure, an efficient simulation, one needing little computation but yielding adequately accurate results for various response quantities, is essential in the search of optimal design in a vast design space. In the present paper, methods of using sensitivities up to 2nd order, and direct application of neural networks are explored. The example problem is how to decide the natural frequencies of a wing given the shape variables of the structure. It is shown that when sensitivities cannot be obtained analytically, the finite difference approach is usually more reliable than a semi-analytical approach provided an appropriate step size is used. The use of second order sensitivities is proved of being able to yield much better results than the case where only the first order sensitivities are used. When neural networks are trained to relate the wing natural frequencies to the shape variables, a negligible computation effort is needed to accurately determine the natural frequencies of a new design.

Kapania, Rakesh K.↗

NASA NTRS2000

Modal Response of Trapezoidal Wing Structures Using Second Order Shape Sensitivities

The modal response of wing structures is very important for assessing their dynamic response including dynamic aeroelastic instabilities. Moreover, in a recent study an efficient structural optimization approach was developed using structural modes to represent the static aeroelastic wing response (both displacement and stress). In this paper, the modal response of general trapezoidal wing structures is approximated using shape sensitivities up to the 2nd order. Also different approaches of computing the derivatives are investigated.

Liu, Youhua↗