Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stiffness”

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 595 records · Page 33

Sensitivity of lag-damping correlations to structural and aerodynamic approximations of isolated experimental rotors in forward flight

The predictions of regressive lag-mode damping levels are correlated with the database of an isolated, soft-inplane, three-blade rotor operated untrimmed. The database was generated at the Army Aeroflightdynamics Directorate at Ames. The correlation covers a broad range of data, from near-zero thrust conditions in hover to high-thrust and highly stalled conditions in forward flight with advance ratio as high as 0.55 and shaft angle as high as 20 degrees. In the experimental rotor, the airfoil or blade portion has essentially uniform mass and stiffness distributions, but the root flexure has highly nonuniform mass and stiffness distributions. Accordingly, the structural approximations refer to four models of root-flexure-blade assembly. They range from a rigid flap-lag model to three elastic flap-lag-torsion models, which differ in modeling the root flexure. The three models of root-flexure are: three root springs in which the bending-torsion couplings are fully accounted for; a finite-length beam element with some average mass and stiffness distributions such that the fundamental frequencies match those of the experimental model; and accurate modal representation in which the actual mass and stiffness distributions of the experimental root-flexure-blade assembly are used in calculating the nonrotating mode shapes. The four models of root-flexure-blade assembly are referred to as the rigid flap-lag model, spring model, modified model and modal model. For each of these four models of the root-flexure-blade assembly, the predictions are based on the following five aerodynamic theories: ear theory, which accounts for large angle-of-attack and reverse-flow effects on lift, and has constant drag and pitching moment; quasisteady stall theory, which includes quasisteady stall lift, drag and pitching moment characteristics of the airfoil section, dynamics stall theory, which uses the ONERA dynamic stall models of lift, drag and pitching moment; dynamic wake theory, which is based on a finite-state three-dimensional wake model and includes all wake effects including both shed and trailing vorticity; and dynamics and wake theory, which combines both dynamic stall theory and dynamic wake theory and is a relatively complete aerodynamic representation.

Gaonkar, G. H.↗

Distributed parameter modeling of repeated truss structures

A new approach to find homogeneous models for beam-like repeated flexible structures is proposed which conceptually involves two steps. The first step involves the approximation of 3-D non-homogeneous model by a 1-D periodic beam model. The structure is modeled as a 3-D non-homogeneous continuum. The displacement field is approximated by Taylor series expansion. Then, the cross sectional mass and stiffness matrices are obtained by energy equivalence using their additive properties. Due to the repeated nature of the flexible bodies, the mass, and stiffness matrices are also periodic. This procedure is systematic and requires less dynamics detail. The first step involves the homogenization from a 1-D periodic beam model to a 1-D homogeneous beam model. The periodic beam model is homogenized into an equivalent homogeneous beam model using the additive property of compliance along the generic axis. The major departure from previous approaches in literature is using compliance instead of stiffness in homogenization. An obvious justification is that the stiffness is additive at each cross section but not along the generic axis. The homogenized model preserves many properties of the original periodic model.

Wang, Han-Ching↗

Dynamics of a split torque helicopter transmission

Split torque designs, proposed as alternatives to traditional planetary designs for helicopter main rotor transmissions, can save weight and be more reliable than traditional designs. This report presents the results of an analytical study of the system dynamics and performance of a split torque gearbox that uses a balance beam mechanism for load sharing. The Lagrange method was applied to develop a system of equations of motion. The mathematical model includes time-varying gear mesh stiffness, friction, and manufacturing errors. Cornell's method for calculating the stiffness of spur gear teeth was extended and applied to helical gears. The phenomenon of sidebands spaced at shaft frequencies about gear mesh fundamental frequencies was simulated by modeling total composite gear errors as sinusoid functions. Although the gearbox has symmetric geometry, the loads and motions of the two power paths differ. Friction must be considered to properly evaluate the balance beam mechanism. For the design studied, the balance beam is not an effective device for load sharing unless the coefficient of friction is less than 0.003. The complete system stiffness as represented by the stiffness matrix used in this analysis must be considered to precisely determine the optimal tooth indexing position.

Krantz, Timothy L.↗

Thin tailored composite wing for civil tiltrotor

The tiltrotor aircraft is a flight vehicle which combines the efficient low speed (i.e., take-off, landing, and hover) characteristics of a helicopter with the efficient cruise speed of a turboprop airplane. A well-known example of such vehicle is the Bell-Boeing V-22 Osprey. The high cruise speed and range constraints placed on the civil tiltrotor require a relatively thin wing to increase the drag-divergence Mach number which translates into lower compressibility drag. It is required to reduce the wing maximum thickness-to-chord ratio t/c from 23% (i.e., V-22 wing) to 18%. While a reduction in wing thickness results in improved aerodynamic efficiency, it has an adverse effect on the wing structure and it tends to reduce structural stiffness. If ignored, the reduction in wing stiffness leads to susceptibility to aeroelastic and dynamic instabilities which may consequently cause a catastrophic failure. By taking advantage of the directional stiffness characteristics of composite materials the wing structure may be tailored to have the necessary stiffness, at a lower thickness, while keeping the weight low. The goal of this study is to design a wing structure for minimum weight subject to structural, dynamic and aeroelastic constraints. The structural constraints are in terms of strength and buckling allowables. The dynamic constraints are in terms of wing natural frequencies in vertical and horizontal bending and torsion. The aeroelastic constraints are in terms of frequency placement of the wing structure relative to those of the rotor system. The wing-rotor-pylon aeroelastic and dynamic interactions are limited in this design study by holding the cruise speed, rotor-pylon system, and wing geometric attributes fixed. To assure that the wing-rotor stability margins are maintained a more rigorous analysis based on a detailed model of the rotor system will need to ensue following the design study. The skin-stringer-rib type architecture is used for the wing-box structure. The design variables include upper and lower skin ply thicknesses and orientation angles, spar and rib web thicknesses and cap areas, and stringer cross-sectional areas. These design variables will allow the maximum tailoring of the structure to meet the design requirements most efficiently. Initial dynamic analysis has been conducted using MSC/NASTRAN to determine the baseline wing's frequencies and mode shapes. For the design study we intend to use the finite-element based code called WIDOWAC (Wing Design Optimization With Aeroeastic Constraints) that was developed at NASA Langley in early 1970's for airplane wing structural analysis and preliminary design. Currently, the focus is on modification and validation of this code which will be used for the civil tiltrotor design efforts.

Rais-Rohani, Masoud↗

A new look at the simultaneous analysis and design of structures

The minimum weight optimization of structural systems, subject to strength and displacement constraints as well as size side constraints, was investigated by the Simultaneous ANalysis and Design (SAND) approach. As an optimizer, the code NPSOL was used which is based on a sequential quadratic programming (SQP) algorithm. The structures were modeled by the finite element method. The finite element related input to NPSOL was automatically generated from the input decks of such standard FEM/optimization codes as NASTRAN or ASTROS, with the stiffness matrices, at present, extracted from the FEM code ANALYZE. In order to avoid ill-conditioned matrices that can be encountered when the global stiffness equations are used as additional nonlinear equality constraints in the SAND approach (with the displacements as additional variables), the matrix displacement method was applied. In this approach, the element stiffness equations are used as constraints instead of the global stiffness equations, in conjunction with the nodal force equilibrium equations. This approach adds the element forces as variables to the system. Since, for complex structures and the associated large and very sparce matrices, the execution times of the optimization code became excessive due to the large number of required constraint gradient evaluations, the Kreisselmeier-Steinhauser function approach was used to decrease the computational effort by reducing the nonlinear equality constraint system to essentially a single combined constraint equation. As the linear equality and inequality constraints require much less computational effort to evaluate, they were kept in their previous form to limit the complexity of the KS function evaluation. To date, the standard three-bar, ten-bar, and 72-bar trusses have been tested. For the standard SAND approach, correct results were obtained for all three trusses although convergence became slower for the 72-bar truss. When the matrix displacement method was used, correct results were still obtained, but the execution times became excessive due to the large number of constraint gradient evaluations required. Using the KS function, the computational effort dropped, but the optimization seemed to become less robust. The investigation of this phenomenon is continuing. As an alternate approach, the code MINOS for the optimization of sparse matrices can be applied to the problem in lieu of the Kreisselmeier-Steinhauser function. This investigation is underway.

Striz, Alfred G.↗

Post-impact fatigue of cross-plied, through-the-thickness reinforced carbon/epoxy composites

An experimental investigation of the post-impact fatigue response of integrally woven carbon/epoxy composites was conducted. Five different through-the-thickness (TTT) reinforcing fibers were used in an experimental textile process that produced an integrally woven (0/90/0/90/0/90/0/90/0)(sub T) ply layup with 21K AS4 carbon tow fiber. The resin was Hercules 3501-6, and the five TTT reinforcing fibers were Kevlar, Toray carbon, AS4 carbon, glass, and IM6 carbon. The purpose of this investigation was to study the post-impact fatigue response of these material systems and to identify the optimum TTT fiber. Samples were impacted with one half inch diameter aluminum balls with an average velocity of 543 ft/sec. Post-impact static compression and constant amplitude tension-compression fatigue tests were conducted. Fatigue tests were conducted with a loading ratio of R=-5, and frequency of 4 Hz. Damage growth was monitored using x-radiographic and sectioning techniques and by examining the stress-strain response (across the impact site) throughout the fatigue tests. The static compressive stress versus far-field strain response was nearly linear for all material groups. All the samples had a transverse shear failure mode. The average compressive modulus (from far-field strain) was about 10 Msi. The average post-impact static compressive strength was about 35.5 Ksi. The IM6 carbon sample had a strength of over 40 Ksi, more than 16 percent stronger than average. There was considerable scatter in the S-N data. However, the IM6 carbon samples clearly had the best fatigue response. The response of the other materials, while worse than IM6 carbon, could not be ranked definitively. The initial damage zones caused by the impact loading and damage growth from fatigue loading were similar for all five TTT reinforcing materials. The initial damage zones were circular and consisted of delaminations, matrix cracks and ply cracks. Post-impact fatigue loading caused delamination growth, ply cracking and fiber bundle failures, typically 45 deg from impact load direction. During the initial 97 percent of fatigue life, delaminations, ply cracks and fiber bundle failures primarily grew at and near the impact site. During the final 3 percent of life, damage grew rapidly transverse to the loading direction as a through-the-thickness transverse shear failure. The stress-strain response was typically linear during the initial 50 percent of life, and stiffness dropped about 20 percent during this period. During the next 47 percent of life, stiffness dropped about 34 percent, and the stress-strain response was no longer linear. The stiffness decreased about 23 percent during the final 3 percent of life. These trends were typical of all the materials tested. Therefore, by monitoring stiffness loss, fatigue failure could be accurately anticipated.

Serdinak, Thomas E.↗

Aspects of passive magnetic levitation based on high-T(sub c) superconducting YBCO thin films

Passive magnetic levitation systems reported in the past were mostly confined to bulk superconducting materials. Here we present fundamental studies on magnetic levitation employing cylindrical permanent magnets floating above high-T(sub c) superconducting YBCO thin films (thickness about 0.3 mu m). Experiments included free floating rotating magnets as well as well-established flexible beam methods. By means of the latter, we investigated levitation and drag force hysteresis as well as magnetic stiffness properties of the superconductor-magnet arrangement. In the case of vertical motion of the magnet, characteristic high symmetry of repulsive (approaching) and attractive (withdrawing) branches of the pronounced force-displacement hysteresis could be detected. Achievable force levels were low as expected but sufficient for levitation of permanent magnets. With regard to magnetic stiffness, thin films proved to show stiffness-force ratios about one order of magnitude higher than bulk materials. Phenomenological models support the measurements. Regarding the magnetic hysteresis of the superconductor, the Irie-Yamafuji model was used for solving the equation of force balance in cylindrical coordinates allowing for a macroscopic description of the superconductor magnetization. This procedure provided good agreement with experimental levitation force and stiffness data during vertical motion. For the case of (lateral) drag force basic qualitative characteristics could be recovered, too. It is shown that models, based on simple asymmetric magnetization of the superconductor, describe well asymptotic transition of drag forces after the change of the magnet motion direction. Virgin curves (starting from equilibrium, i.e. symmetric magnetization) are approximated by a linear approach already reported in literature only. This paper shows that basic properties of superconducting thin films allow for their application to magnetic levitation or - without need of levitation forces, e.g. microgravity - magnetic damping devices.

Schoenhuber, P.↗

Parametric studies of stitching effectiveness for preventing substructure disbond

A methodology is desired that will allow a designer to select appropriate amounts of through-thickness reinforcement needed to meet design requirements. The goal is to use a relatively simple analysis to minimize the amount of testing that needs to be performed, and to make test results from simple configurations applicable to more general structures. Using this methodology one should be able to optimize the selection of stitching materials, the weight of the yarn, and the stitching density. The analysis approach is to treat substructure disbond as a crack propagation problem. In this approach, the stitches have little influence until a delamination begins to grow. Once the delamination reaches, or extends beyond a stitch, the stitch serves to reduce the strain-energy-release-rate (G) at the crack tip for a given applied load. The reduced G can then be compared to the unstitched materials toughness to predict the load required to further extend the crack. The current model treats the stitch as a simple spring which responds to displacements in the vertical (through-thickness) direction. In concept, this approach is similar to that proposed by other authors. Test results indicate that the model should be refined to include the shearing stiffness of the stitch. The strain-energy-release-rate calculations are performed using a code which uses interconnected higher-order plates to model built-up composite cross-sections. When plates are stacked vertically, the interfacial tractions between the plates can be computed. The plate differential equations are solved in closed-form. The code, called SUBLAM, was developed as part of this section in one dimension. Because of this limitation, rows of stitches are treated as a two-dimensional sheet. The spring stiffness of a row of stitches can be estimated from the stitch material, weight, and density. As a practical and conservative approach, we can assume that the stitch is bonded until a crack passes the stitch location. After the crack passes, it is fully bonded. A series of tests were performed to exercise this methodology and incorporated an attached flange such that the sudden change in thickness initiated a delamination. The analysis was used to estimate the materials' critical G from that of the unstitched specimens. With this data, a prediction was made for the load required to delaminate the stitched specimens. Using the methodology, design charts have been created for simplified geometries. These charts give stitch force, along with G(sub 1) and G(sub 2) as as function of the stitch spring stiffness. Using the charts, it should be possible to determine the stitch spring stiffness and strength required to reduce the G to a desired level. From these parameters, the actual stitching material, weight, and density can be computed.

Flanagan, Gerry↗

Passive Magnetic Bearing With Ferrofluid Stabilization

A new class of magnetic bearings is shown to exist analytically and is demonstrated experimentally. The class of magnetic bearings utilize a ferrofluid/solid magnet interaction to stabilize the axial degree of freedom of a permanent magnet radial bearing. Twenty six permanent magnet bearing designs and twenty two ferrofluid stabilizer designs are evaluated. Two types of radial bearing designs are tested to determine their force and stiffness utilizing two methods. The first method is based on the use of frequency measurements to determine stiffness by utilizing an analytical model. The second method consisted of loading the system and measuring displacement in order to measure stiffness. Two ferrofluid stabilizers are tested and force displacement curves are measured. Two experimental test fixtures are designed and constructed in order to conduct the stiffness testing. Polynomial models of the data are generated and used to design the bearing prototype. The prototype was constructed and tested and shown to be stable. Further testing shows the possibility of using this technology for vibration isolation. The project successfully demonstrated the viability of the passive magnetic bearing with ferrofluid stabilization both experimentally and analytically.

Jansen, Ralph↗

A Feasibility Study on a Parallel Mechanism for Examining the Space Shuttle Orbiter Payload Bay Radiators

The goal of the project was to develop the necessary analysis tools for a feasibility study of a cable suspended robot system for examining the space shuttle orbiter payload bay radiators These tools were developed to address design issues such as workspace size, tension requirements on the cable, the necessary accuracy and resolution requirements and the stiffness and movement requirements of the system. This report describes the mathematical models for studying the inverse kinematics, statics, and stiffness of the robot. Each model is described by a matrix. The manipulator Jacobian was also related to the stiffness matrix, which characterized the stiffness of the system. Analysis tools were then developed based on the singular value decomposition (SVD) of the corresponding matrices. It was demonstrated how the SVD can be used to quantify the robot's performance and to provide insight into different design issues.

Roberts, Rodney G.↗

Frictionless Contact of Multilayered Composite Half Planes Containing Layers With Complex Eigenvalues

A previously developed local-global stiffness matrix methodology for the response of a composite half plane, arbitrarily layered with isotropic, orthotropic or monoclinic plies, to indentation by a rigid parabolic punch is further extended to accommodate the presence of layers with complex eigenvalues (e.g., honeycomb or piezoelectric layers). First, a generalized plane deformation solution for the displacement field in an orthotropic layer or half plane characterized by complex eigenvalues is obtained using Fourier transforms. A local stiffness matrix in the transform domain is subsequently constructed for this class of layers and half planes, which is then assembled into a global stiffness matrix for the entire multilayered half plane by enforcing continuity conditions along the interfaces. Application of the mixed boundary condition on the top surface of the half plane indented by a rigid punch results in an integral equation for the unknown pressure in the contact region. The integral possesses a divergent kernel which is decomposed into Cauchy-type and regular parts using the asymptotic properties of the local stiffness matrix and a relationship between Fourier and finite Hilbert transform of the contact pressure. The solution of the resulting singular integral equation is obtained using a collocation technique based on the properties of orthogonal polynomials developed by Erdogan and Gupta. Examples are presented that illustrate the important influence of low transverse properties of layers with complex eigenvalues, such as those exhibited by honeycomb, on the load versus contact length response and contact pressure distributions for half planes containing typical composite materials.

Zhang, Wang↗

Analytic Sensitivity and Approximation of Skin Buckling Constraints in Wing-Shape Synthesis

Explicit expressions for terms of the stiffness and geometric stiffness matrices are derived for the buckling analysis of trapezoidal fiber composite wing skin panels. The formulation is based on Ritz analysis using simple polynomials, and leads to explicit expressions for the analytic sensitivities of the stiffness and geometric stiffness matrices with respect to layer thickness, fiber directions, and panel shape. Integration with wing box analysis using either the equivalent plate approach or the finite element method, makes it possible to obtain sensitivities of panel buckling constraints with respect to wing planform shape or locations of internal ribs and spars. The analytic sensitivities are used to construct approximations of panel buckling constraints for integrated wing/panel design synthesis.

Livne, Eli↗

Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations

Additive Runge-Kutta (ARK) methods are investigated for application to the spatially discretized one-dimensional convection-diffusion-reaction (CDR) equations. First, accuracy, stability, conservation, and dense output are considered for the general case when N different Runge-Kutta methods are grouped into a single composite method. Then, implicit-explicit, N = 2, additive Runge-Kutta ARK2 methods from third- to fifth-order are presented that allow for integration of stiff terms by an L-stable, stiffly-accurate explicit, singly diagonally implicit Runge-Kutta (ESDIRK) method while the nonstiff terms are integrated with a traditional explicit Runge-Kutta method (ERK). Coupling error terms are of equal order to those of the elemental methods. Derived ARK2 methods have vanishing stability functions for very large values of the stiff scaled eigenvalue, z(exp [I]) goes to infinity, and retain high stability efficiency in the absence of stiffness, z(exp [I]) goes to zero. Extrapolation-type stage-value predictors are provided based on dense-output formulae. Optimized methods minimize both leading order ARK2 error terms and Butcher coefficient magnitudes as well as maximize conservation properties. Numerical tests of the new schemes on a CDR problem show negligible stiffness leakage and near classical order convergence rates. However, tests on three simple singular-perturbation problems reveal generally predictable order reduction. Error control is best managed with a PID-controller. While results for the fifth-order method are disappointing, both the new third- and fourth-order methods are at least as efficient as existing ARK2 methods while offering error control and stage-value predictors.

Kennedy, Christopher A.↗

Multibody Parachute Flight Simulations for Planetary Entry Trajectories Using "Equilibrium Points"

A method has been developed to reduce numerical stiffness and computer CPU requirements of high fidelity multibody flight simulations involving parachutes for planetary entry trajectories. Typical parachute entry configurations consist of entry bodies suspended from a parachute, connected by flexible lines. To accurately calculate line forces and moments, the simulations need to keep track of the point where the flexible lines meet (confluence point). In previous multibody parachute flight simulations, the confluence point has been modeled as a point mass. Using a point mass for the confluence point tends to make the simulation numerically stiff, because its mass is typically much less that than the main rigid body masses. One solution for stiff differential equations is to use a very small integration time step. However, this results in large computer CPU requirements. In the method described in the paper, the need for using a mass as the confluence point has been eliminated. Instead, the confluence point is modeled using an "equilibrium point". This point is calculated at every integration step as the point at which sum of all line forces is zero (static equilibrium). The use of this "equilibrium point" has the advantage of both reducing the numerical stiffness of the simulations, and eliminating the dynamical equations associated with vibration of a lumped mass on a high-tension string.

Raiszadeh, Ben↗

Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations

Additive Runge-Kutta (ARK) methods are investigated for application to the spatially discretized one- dimensional convection-diffusion-reaction (CDR) equations. Accuracy, stability, conservation, and dense-output are first considered for the general case when N different Runge-Kutta methods are grouped into a single composite method. Then, implicit-explicit, (N = 2), additive Runge-Kutta (ARK(sub 2)) methods from third- to fifth-order are presented that allow for integration of stiff terms by an L-stable, stiffly-accurate explicit, singly diagonally implicit Runge-Kutta (ESDIRK) method while the nonstiff terms are integrated with a traditional explicit Runge-Kutta method (ERK). Coupling error terms of the partitioned method are of equal order to those of the elemental methods. Derived ARK(sub 2) methods have vanishing stability functions for very large values of the stiff scaled eigenvalue, z['] yields -infinity, and retain high stability efficiency in the absence of stiffness, z['] yield 0. Extrapolation-type stage- value predictors are provided based on dense-output formulae. Optimized methods minimize both leading order ARK(sub 2) error terms and Butcher coefficient magnitudes as well as maximize conservation properties. Numerical tests of the new schemes on a CDR problem show negligible stiffness leakage and near classical order convergence rates. However, tests on three simple singular-perturbation problems reveal generally predictable order reduction. Error control is best managed with a PID-controller. While results for the fifth-order method are disappointing, both the new third- and fourth-order methods are at least as efficient as existing ARK(sub 2) methods.

Kennedy, Christopher A.↗

Equivalent Plate Analysis of Aircraft Wing with Discrete Source Damage

An equivalent plate procedure is developed to provide a computationally efficient means of matching the stiffness and frequencies of flight vehicle wing structures for prescribed loading conditions. First, the equivalent plate is used to match the stiffness of a stiffened panel without damage and the stiffness of a stiffened panel with damage. For both stiffened panels, the equivalent plate models reproduce the deformation of a corresponding detailed model exactly for the given loading conditions. Once the stiffness was matched, the equivalent plate models were then used to predict the frequencies of the panels. Two analytical procedures using the lumped-mass matrix were used to match the first five frequencies of the corresponding detailed model. In both the procedures, the lumped-mass matrix for the equivalent plate is constructed by multiplying the diagonal terms of the consistent-mass matrix by a proportionality constant. In the first procedure, the proportionality constant is selected such that the total mass of the equivalent plate is equal to that of the detailed model. In the second method, the proportionality constant is selected to minimize the sum of the squares of the errors in a set of pre-selected frequencies between the equivalent plate model and the detailed model. The equivalent plate models reproduced the fundamental first frequency accurately in both the methods. It is observed that changing only the mass distribution in the equivalent plate model did not provide enough flexibility to match all of the frequencies.

Krishnamurthy, T.↗

Elastomeric load sharing device

An elastomeric load sharing device, interposed in combination between a driven gear and a central drive shaft to facilitate balanced torque distribution in split power transmission systems, includes a cylindrical elastomeric bearing and a plurality of elastomeric bearing pads. The elastomeric bearing and bearing pads comprise one or more layers, each layer including an elastomer having a metal backing strip secured thereto. The elastomeric bearing is configured to have a high radial stiffness and a low torsional stiffness and is operative to radially center the driven gear and to minimize torque transfer through the elastomeric bearing. The bearing pads are configured to have a low radial and torsional stiffness and a high axial stiffness and are operative to compressively transmit torque from the driven gear to the drive shaft. The elastomeric load sharing device has spring rates that compensate for mechanical deviations in the gear train assembly to provide balanced torque distribution between complementary load paths of split power transmission systems.

Isabelle, Charles J.↗

Frequencies and Flutter Speed Estimation for Damaged Aircraft Wing Using Scaled Equivalent Plate Analysis

Equivalent plate analysis is often used to replace the computationally expensive finite element analysis in initial design stages or in conceptual design of aircraft wing structures. The equivalent plate model can also be used to design a wind tunnel model to match the stiffness characteristics of the wing box of a full-scale aircraft wing model while satisfying strength-based requirements An equivalent plate analysis technique is presented to predict the static and dynamic response of an aircraft wing with or without damage. First, a geometric scale factor and a dynamic pressure scale factor are defined to relate the stiffness, load and deformation of the equivalent plate to the aircraft wing. A procedure using an optimization technique is presented to create scaled equivalent plate models from the full scale aircraft wing using geometric and dynamic pressure scale factors. The scaled models are constructed by matching the stiffness of the scaled equivalent plate with the scaled aircraft wing stiffness. It is demonstrated that the scaled equivalent plate model can be used to predict the deformation of the aircraft wing accurately. Once the full equivalent plate geometry is obtained, any other scaled equivalent plate geometry can be obtained using the geometric scale factor. Next, an average frequency scale factor is defined as the average ratio of the frequencies of the aircraft wing to the frequencies of the full-scaled equivalent plate. The average frequency scale factor combined with the geometric scale factor is used to predict the frequency response of the aircraft wing from the scaled equivalent plate analysis. A procedure is outlined to estimate the frequency response and the flutter speed of an aircraft wing from the equivalent plate analysis using the frequency scale factor and geometric scale factor. The equivalent plate analysis is demonstrated using an aircraft wing without damage and another with damage. Both of the problems show that the scaled equivalent plate analysis can be successfully used to predict the frequencies and flutter speed of a typical aircraft wing.

Krishnamurthy, Thiagarajan↗