Open Access. Powered by Scholars. Published by Universities.®

Mechanical Engineering Commons

Open Access. Powered by Scholars. Published by Universities.®

Articles 1 - 30 of 52

Full-Text Articles in Mechanical Engineering

Experimental Investigation Of Active Control Of Bluff Body Vortex Shedding, Ilteris Koc Jan 2008

Experimental Investigation Of Active Control Of Bluff Body Vortex Shedding, Ilteris Koc

Mechanical & Aerospace Engineering Theses & Dissertations

Mean and fluctuating forces acting on a body are strongly related to vortex shedding generated behind it. Therefore, it is possible to obtain substantial reductions of at least the unsteady forces if vortex shedding is controlled or its regularity is reduced. While conventional active flow control methods are mainly concerned with direct interaction with, and alteration of, the mean flow about a body, modern techniques involve altering existing flow instabilities using relatively small inputs to obtain large-scale changes of mean flows. Aerodynamic flow control may be intended to delay or suppress boundary layer separation through creation of a boundary layer …


Micron-Level Actuator For Thermal-Fluid Control In Microchannels, Nurhak Erbas Jul 2006

Micron-Level Actuator For Thermal-Fluid Control In Microchannels, Nurhak Erbas

Mechanical & Aerospace Engineering Theses & Dissertations

Effectiveness of an actuator is investigated for thermal-flow control in microchannels. First, simulations of a single actuator in a quiescent external medium are performed in order to study the parameters characterizing the synthetic jet flow from the actuator. For this purpose, a simplified, two-dimensional configuration is considered. The membrane motion is modeled in a realistic manner as a moving boundary in order to accurately compute the flow inside the actuator cavity. The geometric and actuation parameters of the actuator are investigated to define the effectiveness of the jet flow. The study is done initially at macro scales. Then, the flow …


Boundary Interference Assessment And Correction For Open Jet Wind Tunnels Using Panel Methods, Wael Ahmed Mokhtar Apr 2006

Boundary Interference Assessment And Correction For Open Jet Wind Tunnels Using Panel Methods, Wael Ahmed Mokhtar

Mechanical & Aerospace Engineering Theses & Dissertations

The presence of nearby boundaries in a wind tunnel can lead to aerodynamic measurements on a model in the wind tunnel that differ from those that would be made when the boundaries of the moving fluid were infinitely far away. The differences, referred to as boundary interference or wall interference, can be quite large, such as when testing aircraft models developing high lift forces, or whose wingspan is a large fraction of the wind tunnel width, or high drag models whose frontal area is a large fraction of the tunnel cross section. Correction techniques for closed test section (solid walled) …


Nonlinear Flutter Of Curved Panels Under Yawed Supersonic Flow Using Finite Elements, Mohamed Salim Azzouz Jan 2005

Nonlinear Flutter Of Curved Panels Under Yawed Supersonic Flow Using Finite Elements, Mohamed Salim Azzouz

Mechanical & Aerospace Engineering Theses & Dissertations

In the extensive published literature on panel flutter, a large number of papers are dedicated to investigation of flat plates in the supersonic flow regime. Very few authors have extended their work to flutter of curved panels. The curved geometry generates a pre-flutter behavior, triggering a static deflection due to a static aerodynamic load (SAL) over the panel as well as dynamic characteristics unique to this geometry. The purpose of this dissertation is to provide new insights in the subject of flutter of curved panels. Finite element frequency and time domain methods are developed to predict the pre/post flutter responses …


Computation Of Shock Induced Noise In Imperfectly Expanded Supersonic Jets, Bulent Imamoglu Jan 2004

Computation Of Shock Induced Noise In Imperfectly Expanded Supersonic Jets, Bulent Imamoglu

Mechanical & Aerospace Engineering Theses & Dissertations

Screech noise exists only in imperfectly expanded jets. The exit pressure of imperfectly expanded jets does not match ambient pressure, so expansion or compression waves appear out of the nozzle and generate shock cell patterns. Screech is generated by the interaction of shock cells and instability waves. Many experiments and computations have been done to model screech noise, but it is not yet a very well known subject.

A numerical study is performed to understand screech generation mechanisms and to compare with latest experiments. A supersonic underexpanded jet of 25.4 mm diameter is modeled for cases of Mach numbers of …


Nonlinear Response And Fatigue Estimation Of Aerospace Curved Surface Panels To Acoustic And Thermal Loads, Adam Przekop Jul 2003

Nonlinear Response And Fatigue Estimation Of Aerospace Curved Surface Panels To Acoustic And Thermal Loads, Adam Przekop

Mechanical & Aerospace Engineering Theses & Dissertations

This work presents a finite element modal formulation for large amplitude free vibration of arbitrary laminated composite shallow shells. The system equations of motion are formulated first in the physical structural-node degrees of freedom (DOF). Then, the system is transformed into general Duffing-type modal equations with modal amplitudes of coupled linear bending-inplane modes. The linear bending-inplane coupling is due to the shell curvature as well as unsymmetric lamination stacking. Multiple modes, inplane inertia, and the first-order transverse shear deformation for composites are considered in the formulation. A triangular shallow shell finite element is developed from an extension of the triangular …


New Approaches To Gortler And Tollmien-Schlichting Boundary-Layer Instabilities, Vijay Kalburgi Jan 2003

New Approaches To Gortler And Tollmien-Schlichting Boundary-Layer Instabilities, Vijay Kalburgi

Mechanical & Aerospace Engineering Theses & Dissertations

A detailed understanding of the physical processes that lead to instability and transition in the external viscous boundary layers of airfoils has been the goal of researchers for decades. This study seeks to enhance our understanding of Gortler and Tollmien-Schlichting instabilities through a critical review of prior theoretical methodology and fundamental experiments, the development of dimensionally consistent disturbance equations, and novel interpretation of the resulting flow physics revealed by their solutions to explain the process of transition.

Gortler vortices arise in boundary layers along concave surfaces due to centrifugal effects and these vortices in combination with other instabilities play an …


Numerical Study Of Two-Dimensional Secondary Injection Into A Mach 3.5 Freestream, Stephen C. Coghill Jan 2003

Numerical Study Of Two-Dimensional Secondary Injection Into A Mach 3.5 Freestream, Stephen C. Coghill

Mechanical & Aerospace Engineering Theses & Dissertations

Upstream interaction within a dual-mode scramjet engine is investigated numerically. The upstream interaction is investigated by increasing the freestream-to-injector pressure ratio. The constant area duct and sudden expansion geometries are used to achieve these results. Analysis using the symmetry boundary condition is compared to the full boundary condition for the constant area duct. Numerical analysis of the Mach 3.5 freestream is conducted using normal sonic injection of nitrogen gas to create the upstream interaction. Comparisons, where applicable, are made to experimental results. Due to the high speed of the flowfield, oblique shockwaves are present causing numerical convergence difficult to achieve. …


Modeling Aspects Of Magnetic Actuators And Magnetic Suspension Systems, V. Dale Bloodgood Jr. Apr 2002

Modeling Aspects Of Magnetic Actuators And Magnetic Suspension Systems, V. Dale Bloodgood Jr.

Mechanical & Aerospace Engineering Theses & Dissertations

This dissertation is a study of new modeling techniques developed for magnetic suspension systems. The techniques discussed are modifications of magnetic circuit theory and fundamental eddy current models. The techniques are compared against experimental test results and finite element data. The information gained from the experimental testing is used to provide insight into magnetic bearing design.

A small-gap modeling technique called extended circuit theory is developed that incorporates information about the system gained from finite element data, or experimental data, to be included in the analytic model. The variations between the classical magnetic circuit model and the finite element model …


Supersonic Combustion And Mixing Characteristics Of Hydrocarbon Fuels In Screamjet Engines, Ahmed A. Taha Jan 2002

Supersonic Combustion And Mixing Characteristics Of Hydrocarbon Fuels In Screamjet Engines, Ahmed A. Taha

Mechanical & Aerospace Engineering Theses & Dissertations

The combustion characteristics of gaseous propane in supersonic airflow using the rearward-facing step that is swept inward from both end sides is studied. The effect of sweeping the step on the flow field features of propane combustion is investigated.

The study of the supersonic combustion of ethylene is carried out using different combustor configurations, different main fuel equivalence ratios, and different pilot fuel equivalence ratios.

The swept step shows the ability to hold the propane flame in the supersonic air stream without extinction. It was found that the side sweeping of the combustor exhibits the high temperature and combustion products …


Computational Modeling Of Airborne Noise Demonstrated Via Benchmarks, Supersonic Jet, And Railway Barrier, Moumen Idres Jan 1999

Computational Modeling Of Airborne Noise Demonstrated Via Benchmarks, Supersonic Jet, And Railway Barrier, Moumen Idres

Mechanical & Aerospace Engineering Theses & Dissertations

In the last several years, there has been a growing demand for mobility to cope with the increasing population. All kinds of transportation have responded to this demand by expanding their networks and introducing new ideas. Rail transportation introduced the idea of high-speed trains and air transportation introduced the idea of high-speed civil transport (HSCT). In this expanding world, the noise legislation is felt to inhibit these plans. Accurate computational methods for noise prediction are in great demand.

In the current research, two computational methods are developed to predict noise propagation in air. The first method is based on the …


Continuous Adjoint Sensitivity Analysis For Aerodynamic And Acoustic Optimization, Kaveh Ghayour Jan 1999

Continuous Adjoint Sensitivity Analysis For Aerodynamic And Acoustic Optimization, Kaveh Ghayour

Mechanical & Aerospace Engineering Theses & Dissertations

A gradient-based shape optimization methodology based on continuous adjoint sensitivities has been developed for two-dimensional steady Euler equations on unstructured meshes and the unsteady transonic small disturbance equation. The continuous adjoint sensitivities of the Helmholtz equation for acoustic applications have also been derived and discussed.

The highlights of the developments for the steady two-dimensional Euler equations are the generalization of the airfoil surface boundary condition of the adjoint system to allow a proper closure of the Lagrangian functional associated with a general cost functional and the results for an inverse problem with density as the prescribed target. Furthermore, it has …


Skin-Stiffener Separation In Buckled Composite Plates And Shallow Shells, Mark Edward Robeson Oct 1998

Skin-Stiffener Separation In Buckled Composite Plates And Shallow Shells, Mark Edward Robeson

Mechanical & Aerospace Engineering Theses & Dissertations

The development of an approximate analysis for predicting stresses in the skin-stiffener interface region of composite plates and shallow shells is presented. The analysis determines interlaminar normal stress and transverse shear stresses in the direction of the stiffener axial coordinate and the direction of the inplane coordinate perpendicular to the stiffener axis. The analysis accounts for skin-stiffener interfaces with nonzero thickness and for curvature in the skin and stiffener. The interlaminar normal stress and transverse shear stresses at the interface are applied to the lower face of the skin and the upper face of the stiffener flange as unknown functions, …


Vortex Wake And Exhaust Plume Interaction, Including Ground Effect, Ihab Gaber Adam Jul 1998

Vortex Wake And Exhaust Plume Interaction, Including Ground Effect, Ihab Gaber Adam

Mechanical & Aerospace Engineering Theses & Dissertations

Computational modeling and studies of the near-field wake-vortex turbulent flows, far-field turbulent wake-vortex/exhaust-plume interaction for subsonic and High Speed Civil Transport (HSCT) airplane, and wake-vortex/exhaust-plume interaction with the ground are carried out. The three-dimensional, compressible Reynolds-Averaged Navier-Stokes (RANS) equations are solved using the implicit, upwind, Roe-flux-differencing, finite-volume scheme. The turbulence models of Baldwin and Lomax, one-equation model of Spalart and Allmaras and two-equation shear stress transport model of Menter are implemented with the RANS solver for turbulent-flow modeling.

For the near-field study, computations are carried out on a fine grid for a rectangular wing with a NACA-0012 airfoil section and …


Placement Of Piezoelectric Actuators For Active Control Of Vibration Using Modal Parameters, Xuegeng Zhu Jan 1998

Placement Of Piezoelectric Actuators For Active Control Of Vibration Using Modal Parameters, Xuegeng Zhu

Mechanical & Aerospace Engineering Theses & Dissertations

An equation is derived to model the piezoelectric actuators incorporation with flexible structures. This equation permits the comparison of the performance indices over the entire structure for a piezoelectric actuator with constant area, which is unachievable if the Finite Element Method is used for complicated structures.

An index has been developed for placement of piezoelectric actuator for control of vibration of a flexible structure. This index is derived from the definition of H2norm. Computation of the proposed index requires only the natural frequencies and corresponding mode shapes of the structures of interest. The method is well suited to large …


Three-Dimensional Aerodynamic Design Optimization Using Discrete Sensitivity Analysis And Parallel Computing, Amidu Olawale Oloso Apr 1997

Three-Dimensional Aerodynamic Design Optimization Using Discrete Sensitivity Analysis And Parallel Computing, Amidu Olawale Oloso

Mechanical & Aerospace Engineering Theses & Dissertations

A hybrid automatic differentiation/incremental iterative method was implemented in the general purpose advanced computational fluid dynamics code (CFL3D Version 4.1) to yield a new code (CFL3D.ADII) that is capable of computing consistently discrete first order sensitivity derivatives for complex geometries. With the exception of unsteady problems, the new code retains all the useful features and capabilities of the original CFL3D flow analysis code. The superiority of the new code over a carefully applied method of finite-differences is demonstrated.

A coarse grain, scalable, distributed-memory, parallel version of CFL3D.ADII was developed based on "derivative stripmining". In this data-parallel approach, an identical copy …


Input Design For Systems Under Identification Using Indirect And Direct Methods, Marco P. Schoen Mar 1997

Input Design For Systems Under Identification Using Indirect And Direct Methods, Marco P. Schoen

Mechanical & Aerospace Engineering Theses & Dissertations

The motivation for system identification can be manifold. In this work, the provocation to identify unknown system characteristics is derived from the control engineering point of view. That is, one intends to design a control strategy based on the identified system properties. The used system identification methods are the Open-Loop Kalman filter System Identification method (OKID) and the Closed-Loop System Identification method (CLID). It is shown that the quantitative largest error of the system identification is given by its model representation, that is the attempt to describe a system with model parameters which poses a linear relationship with the input/output …


A Progressive Damage Methodology For Residual Strength Predictions Of Center-Crack Tension Composite Panels, Timothy William Coats Jul 1996

A Progressive Damage Methodology For Residual Strength Predictions Of Center-Crack Tension Composite Panels, Timothy William Coats

Mechanical & Aerospace Engineering Theses & Dissertations

An investigation of translaminate fracture and a progressive damage methodology was conducted to evaluate and develop residual strength prediction capability for laminated composites with through penetration notches. This is relevant to the damage tolerance of an aircraft fuselage that might suffer an in-flight accident such as an uncontained engine failure. An experimental characterization of several composite materials systems revealed an R-curve type of behavior. Fractographic examinations led to the postulate that this crack growth resistance could be due to fiber bridging, defined here as fractured fibers of one ply bridged by intact fibers of an adjacent ply.

The progressive damage …


Numerical Prediction Of Turbulent Diffusion Flames Formed By Cylindrical Tube Injector, Ali S. Kheireddine Jul 1996

Numerical Prediction Of Turbulent Diffusion Flames Formed By Cylindrical Tube Injector, Ali S. Kheireddine

Mechanical & Aerospace Engineering Theses & Dissertations

This work summarizes numerical results for a diffusion flame formed from a cylindrical tube fuel injector, issuing gaseous fuel jet vertically in a quiescent atmosphere. Both pure fuels as well as fuel mixtures are examined. The primary objective is to predict the flame base height as a function of the jet velocity. A finite volume scheme is used to discretize the time-averaged Navier-Stokes equations for the reacting flow, resulting from the turbulent fuel jet motion. The turbulent stresses, and heat and mass fluxes are computed from the Reynolds stress turbulence model. A chemical kinetics model involving a two-step chemical reaction …


Unsteady, Transonic Flow Around Delta Wings Undergoing Coupled And Natural Modes Response: A Multidisciplinary Problem, Margaret Anne Menzies Apr 1996

Unsteady, Transonic Flow Around Delta Wings Undergoing Coupled And Natural Modes Response: A Multidisciplinary Problem, Margaret Anne Menzies

Mechanical & Aerospace Engineering Theses & Dissertations

The unsteady, three-dimensional Navier-Stokes equations coupled with the Euler equations of rigid-body dynamics are sequentially solved to simulate and analyze the aerodynamic response of a high angle of attack delta wing undergoing oscillatory motion. The governing equations of fluid flow and dynamics of the multidisciplinary problem are solved using a time-accurate solution of the laminar, unsteady, compressible, full Navier-Stokes equations with the implicit, upwind, Roe flux-difference splitting, finite-volume scheme and a four-state Runge-Kutta scheme, respectively. The primary model under consideration consists of a 65° swept, sharp-edged, cropped delta wing of zero thickness at 20° angle of attack. In a freestream …


Variational Methods In Sensitivity Analysis And Optimization For Aerodynamic Applications, Adem Ibrahim Hussen Apr 1996

Variational Methods In Sensitivity Analysis And Optimization For Aerodynamic Applications, Adem Ibrahim Hussen

Mechanical & Aerospace Engineering Theses & Dissertations

Variational methods (VM) sensitivity analysis, which is the continuous alternative to the discrete sensitivity analysis, is employed to derive the costate (adjoint) equations, the transversality conditions, and the functional sensitivity derivatives. In the derivation of the sensitivity equations, the variational methods use the generalized calculus of variations, in which the variable boundary is considered as the design function. The converged solution of the state equations together with the converged solution of the costate equations are integrated along the domain boundary to uniquely determine the functional sensitivity derivatives with respect to the design function.

The determination of the sensitivity derivatives of …


Sensitivity Analysis And Optimization Of Aerodynamic Configurations With Blend Surfaces, Almuttil Mathew Thomas Apr 1996

Sensitivity Analysis And Optimization Of Aerodynamic Configurations With Blend Surfaces, Almuttil Mathew Thomas

Mechanical & Aerospace Engineering Theses & Dissertations

A novel (geometrical) parametrization procedure using solutions to a suitably chosen fourth order partial differential equation is used to define a class of airplane configurations. Inclusive in this definition are surface grids, volume grids, and grid sensitivity. The general airplane configuration has wing, fuselage, vertical tail and horizontal tail. The design variables are incorporated into the boundary conditions, and the solution is expressed as a Fourier series. The fuselage has circular cross section, and the radius is an algebraic function of four design parameters and an independent computational variable. Volume grids are obtained through an application of the Control Point …


Numerical Simulation Of Complex, Three-Dimensional, Turbulent-Free Jets, Robert V. Wilson Jan 1996

Numerical Simulation Of Complex, Three-Dimensional, Turbulent-Free Jets, Robert V. Wilson

Mechanical & Aerospace Engineering Theses & Dissertations

Three-dimensional, incompressible turbulent jets with rectangular and elliptical cross-section are simulated with a finite-difference numerical method. The full Navier-Stokes equations are solved at low Reynoids numbers, whereas at high Reynolds numbers filtered forms of the equations are solved along with a sub-grid scale model to approximate the effects of the unresolved scales. A 2-N storage, third-order Runge-Kutta scheme is used for temporal discretization and a fourth-order compact scheme is used for spatial discretization. Although such methods are widely used in the simulation of compressible flows, the lack of an evolution equation for pressure or density presents particular difficulty in incompressible …


An Experimental Investigation Of Wall Cooling Effects On Hypersonic Boundary Layer Stability In A Quiet Wind Tunnel, Alan E. Blanchard Oct 1995

An Experimental Investigation Of Wall Cooling Effects On Hypersonic Boundary Layer Stability In A Quiet Wind Tunnel, Alan E. Blanchard

Mechanical & Aerospace Engineering Theses & Dissertations

One of the primary reasons for developing quiet tunnels is for the investigation of high-speed boundary-layer stability and transition phenomena without the transition-promoting effects of acoustic radiation from tunnel walls. In this experiment, a flared-cone model under adiabatic- and cooled-wall conditions was placed in a calibrated, 'quiet' Mach 6 flow and the stability of the boundary layer was investigated using a prototype constant-voltage anemometer. The results were compared with linear-stability theory predictions and good agreement was found in the prediction of second-mode frequencies and growth. In addition, the same 'N=10' criterion used to predict boundary-layer transition in subsonic, transonic, and …


Buckling Analysis And Optimum Design Of Multidirectionally Stiffened Composite Curved Panel, Navin R. R. Jaunky Oct 1995

Buckling Analysis And Optimum Design Of Multidirectionally Stiffened Composite Curved Panel, Navin R. R. Jaunky

Mechanical & Aerospace Engineering Theses & Dissertations

Continuous filament grid-stiffened structure is a stiffening concept that combines structural efficiency and damage tolerance. However, buckle resistant design optimization of such structures using a finite element method is expensive and time consuming due to the number of design parameters that can be varied. An analytical optimization procedure which is simple, efficient and supports the preliminary design of grid-stiffened structures for application to combined loading cases is needed.

An analytical model for a general grid-stiffened curved panel is developed using an improved smeared theory with a first-order, shear-deformation theory to account for transverse shear flexibilities and local skin-stiffener interaction effects. …


Aerodynamic Design Optimization With Consistently Discrete Sensitivity Derivatives Via The Incremental Iterative Method, Vamshi M. Korivi Apr 1995

Aerodynamic Design Optimization With Consistently Discrete Sensitivity Derivatives Via The Incremental Iterative Method, Vamshi M. Korivi

Mechanical & Aerospace Engineering Theses & Dissertations

In this study which involves advanced fluid-flow codes, an incremental iterative formulation (also known as the "delta" or "correction" form), together with the well-known spatially split approximate-factorization algorithm, is presented for solving the large, sparse systems of linear equations that are associated with aerodynamic sensitivity analysis. For the smaller two dimensional problems, a direct method can be applied to solve these linear equations in either the standard or the incremental form, in which case the two are equivalent. However, iterative methods are needed for larger two-dimensional and three dimensional applications because direct methods require more computer memory than is currently …


Investigation Of Wall Injectors For Supersonic Mixing Enhancement, Yaacov Haimovitch Jul 1994

Investigation Of Wall Injectors For Supersonic Mixing Enhancement, Yaacov Haimovitch

Mechanical & Aerospace Engineering Theses & Dissertations

A comparative study of the interaction between wall mounted swept-ramp injectors and injector nozzle shape has been conducted in a constant area duct to explore techniques to enhance mixing in scramjet combustors. The scramjet combustors are currently being developed for propulsion systems applications on the envisioned hypersonic vehicles. Short combustor residence time, a requirement for fuel injection parallel to the main flow in the combustor, and an overall strong sensitivity of the vehicle performance to the propulsion system motivated the investigation. The swept-ramp injector investigated in this study produces vortex shedding and local separation downstream of the injector's nozzle exit, …


Unsteady Flow Simulations About Moving Boundary Configurations Using Dynamic Domain Decomposition Techniques, Guan-Wei Yen Apr 1994

Unsteady Flow Simulations About Moving Boundary Configurations Using Dynamic Domain Decomposition Techniques, Guan-Wei Yen

Mechanical & Aerospace Engineering Theses & Dissertations

A computational method is developed to solve the coupled governing equations of an unsteady flowfield and those of rigid-body dynamics in six degrees-of-freedom (6-DOF). This method is capable of simulating the unsteady flowfields around multiple component configurations with at least one of the components in relative motion with respect to the others. Two of the important phenomena that such analyses can help us to understand are the unsteady aerodynamic interference and the boundary-induced component of such a flowfield. By hybridizing two dynamic domain decomposition techniques, the grid generation task is simplified, the computer memory requirement is reduced, and the governing …


Flow Simulations About Steady-Complex And Unsteady Moving Configurations Using Structured-Overlapped And Unstructured Grids, James Charles Newman Iii Apr 1994

Flow Simulations About Steady-Complex And Unsteady Moving Configurations Using Structured-Overlapped And Unstructured Grids, James Charles Newman Iii

Mechanical & Aerospace Engineering Theses & Dissertations

The limiting factor in simulating flows past realistic configurations of interest has been the discretization of the physical domain on which the governing equations of fluid flow may be solved. In an attempt to circumvent this problem, many Computational Fluid Dynamic (CFD) methodologies that are based on different grid generation and domain decomposition techniques have been developed. However, due to the costs involved and expertise required, very few comparative studies between these methods have been performed. In the present work, the two CFD methodologies which show the most promise for treating complex three-dimensional configurations as well as unsteady moving boundary …


Finite Element Analysis For Nonlinear Flutter Suppression Of Composite Panels At Elevated Temperatures Using Piezoelectric Materials, Run Chen Zhou Jan 1994

Finite Element Analysis For Nonlinear Flutter Suppression Of Composite Panels At Elevated Temperatures Using Piezoelectric Materials, Run Chen Zhou

Mechanical & Aerospace Engineering Theses & Dissertations

Nonlinear coupled finite element equations of motion are derived for composite panels with embedded piezoelectric layers subjected to aerodynamic, thermal loads and applied electric fields. The nonlinear equations of motion describe the coupling between a structure and an electrical network through the piezoelectric effect. The von Karman large-deflection strain-displacement relations, quasi-steady first-order piston theory aerodynamics, quasi-steady thermal stress theory and linear piezoelectricity theory are used to formulate the nonlinear coupled panel flutter finite element equations of motion in nodal displacements. The governing equations, which are referred to actuator and sensor equations, form a basis for piezoelectric actuation and sensing. Following …