Active vibration control of piezoelectric bonded smart structures using PID algorithm

Thin-walled structures are sensitive to vibrate under even very small disturbances. In order to design a suitable controller for vibration suppression of thin-walled smart structures, an electro-mechanically coupled finite element(FE) model of smart structures is developed based on first-order shear deformation(FOSD) hypothesis. Considering the vibrations generated by various disturbances, which include free and forced vibrations, a PID control is implemented to damp both the free and forced vibrations. Additionally, an LQR optimal control is applied for comparison.The implemented control strategies are validated by a piezoelectric layered smart plate under various excitations.     

Multiphysics simulation method of Lamb wavepropagation with piezoelectric transducers under load condition

Lamb Wave(LW) simulation under time-varying conditions is an effective and low cost way to study the problem of the low reliability of the structural health monitoring methods based on the LW and Piezoelectric Transducer(PT). In this paper, a multiphysics simulation method of the LW propagation with the PTs under load condition is proposed. With this method, two key mechanisms of the load influence on the LW propagation are considered and coupled with each other. The first mechanism is the acoustoelastic effect which is the main reason of the LW velocity change. The second key mechanism is the load influence on piezoelectric materials, which results in a change of the amplitude. Based on the computational platform of the COMSOL Multiphysics, a multiphysics simulation model of the LW propagation with the PTs under load condition is established. The simulation model includes two physical phenomena. The first one is called solid mechanics, which is used to simulate the acoustoelastic effect being combined with the hyperelastic material properties of the structure in which the LW propagates. The second one is called electromechanical coupling, which considers the simulation of the piezoelectric effect of the PTs for the LW excitation and sensing. To simulate the load influence on piezoelectric materials, a non-linear numerical model of the relationship between the load and the piezoelectric coefficient d31 is established based on an experiment of the load influence on the LW. The simulation results under uniaxial tensile load condition are obtained and are compared with the data obtained from the experiment. It shows that the variations of the phase velocity and amplitude of the LW obtained from the simulation model match the experimental results well.     

Gust load alleviation wind tunnel tests of a large-aspect-ratio flexible wing with piezoelectric control

An active control technique utilizing piezoelectric actuators to alleviate gust-response loads of a large-aspect-ratio flexible wing is investigated. Piezoelectric materials have been exten-sively used for active vibration control of engineering structures. In this paper, piezoelectric mate-rials further attempt to suppress the vibration of the aeroelastic wing caused by gust. The motion equation of the flexible wing with piezoelectric patches is obtained by Hamilton's principle with the modal approach, and then numerical gust responses are analyzed, based on which a gust load alle-viation (GLA) control system is proposed. The gust load alleviation system employs classic propor tional-integral-derivative (PID) controllers which treat piezoelectric patches as control actuators and acceleration as the feedback signal. By a numerical method, the control mechanism that piezo-electric actuators can be used to alleviate gust-response loads is also analyzed qualitatively. Further-more, through low-speed wind tunnel tests, the effectiveness of the gust load alleviation active control technology is validated. The test results agree well with the numerical results. Test results show that at a certain frequency range, the control scheme can effectively alleviate the z and x wing-tip accelerations and the root bending moment of the wing to a certain extent. The control system gives satisfying gust load alleviation efficacy with the reduction rate being generally over 20%.     

压电驱动器的气动弹性应用

 随着压电智能材料与结构的发展,压电驱动器在气动弹性控制领域占据重要地位。使用压电驱动器控制翼面变形,利用而不是抵抗气动弹性效应可以控制升力、力矩以及它们的分布。采用基本相同的智能结构翼面控制系统,根据不同的控制目标需求,使用压电智能材料驱动器可以达到多种目的,包括静态的形状控制与动态的颤振抑制、抖振控制与阵风响应控制。静态控制方面例如改变翼面形状获得附加空气动力以增加升力、提供横滚力矩、改变升力分布以减小诱导阻力或减小翼根弯矩等;动态控制例如利用改变翼面形状产生的附加空气动力作为控制载荷,改变气动弹性系统的耦合程度,根据控制效果要求可作为气动阻尼、气动刚度或气动质量。这种控制方法可以减轻结构重量,提高操纵效率,扩大飞行包线,提高材料利用率,已成为可变形飞行器的重要研究内容。本文主要阐述压电驱动器气动弹性应用的动机与机理、发展与成就以及问题与展望。     

压电智能反射面静态形状控制与作动器位置优化

为了提高反射面的精度,建立了压电智能反射面的形状控制模型,对该结构的力学建模方法、形状控制方法和作动器优化配置算法进行了研究。首先,将蜂窝夹层结构的压电智能反射面等效为多层复合板,根据虚功原理推导了结构的有限元方程。然后,根据建立的有限元方程,推导了反射面变形的均方根误差与作动器控制电压的关系式,以均方根误差最小为优化目标,建立了形状控制优化模型,将作动器控制电压的优化转化为约束优化问题的求解。最后,采用模拟退火算法对压电作动器进行了优化配置。为了验证形状控制的可行性及优化算法的有效性,以300mm口径的平面压电智能反射面为例进行仿真,分析结果表明,通过压电作动器的控制,可以使反射面的重力变形误差减小97%以上,对于给定数目的作动器,通过模拟退火算法优化,可使其布置在最佳的位置。     

面向结构振动控制的压电作动器优化配置研究

 为了研究压电主动结构振动控制当中作动器的位置优化问题,从系统的状态空间方程出发,在系统可控性Gram矩阵特征值的基础上来描述性能指标,以输入的能量吸收率为优化目标函数,提出了一种新的位置优化配置准则。利用有限元分析(FEA)方法分析作动器的配置,并与遗传算法(GA)结合进行优化计算,计算过程中对作动器的位置采用二进制编码加以描述。通过对一压电简支板结构的仿真计算对该方法进行了验证,并与其他几种不同的配置方法进行对比,从而证明了新方法的优越性。     

含脱胶压电传感器/驱动器的智能结构的有限元分析

将压电传感器、驱动器粘贴在复合材料结构的表面可实现结构的振动主动控制 ,但若传感器、驱动器部分脱开会对结构的静、动态特性产生显著影响。建立了一个新的加强假定应变压电固体单元 ,用于压电自适应层合结构振动主动控制的模拟仿真。与现有的压电固体单元相比 ,所建单元性能更优越 ,精度和计算效率更高 ,并能用于壳体结构的分析。采用相同坐标值但不同的节点号的方法模拟脱层 ,利用该单元分析了传感器、驱动器脱开对结构动力学特性的影响。     

1-3型压电纤维主动叠层板扭转特性的研究

对4类典型铺层的1-3型压电纤维复合材料主动板的变形模式进行了理论分析。在此基础上对具有扭转变形模式的主动叠层板进行了振动特性分析,得到了该类型主动叠层板在外界电场激励下的稳态解。最后对该类主动叠层板的扭转振动特性进行了深入的研究,包括作用电场强度的大小、频率对扭转特性的影响、压电纤维铺层角和压电纤维体积比对主动叠层板扭转特性的影响等。算例表明本文计算模型和有限元得到的结果符合良好,同时算例结果表明含1-3型压电纤维复合材料主动叠层板具有良好的主动扭转特性。     

Optimal positioning of piezoelectric actuators on a smart fin using bio-inspired algorithms

In this paper a novel approach is developed for optimization of piezoelectric actuators in vibration suppression. A scaled model of a vertical tail of F/A-18 is developed in which piezoelectric actuators are bounded to the surface. The frequency response function (FRF) of the system is then recorded and maximization of the FRF peaks is considered as the objective function of the optimization algorithm to enhance the actuator authority on the mode, which assigns the optimal placement of the pair of piezoelectric actuators on the smart fin. Six multi-layer perceptron neural networks are employed to perform surface fitting to the discrete data generated by the finite element method (FEM). Invasive weed optimization (IWO), a novel numerical stochastic optimization algorithm, is then employed to maximize the FRF peak which in due reduces the vibration of the smart fin. Results indicated an accurate surface fitting for the FRF peak data as well as the optimal placement of the piezoelectric actuators for vibration suppression.     

Nonlinear modal electromechanical coupling factor for piezoelectric structures containing nonlinearities

Within the linear framework, the Modal Electromechanical Coupling Factor (MEMCF) is an important indicator to quantify the dynamic conversion of mechanical energy and electrical energy of piezoelectric structures. It is also an important tool to guide the piezoelectric damping design of linear structures. Advanced aircraft often fly in maneuvers, and the variable working conditions induce drastic changes in the load level on structures. Geometric and contact nonlinearities of thin-walled structures and joint structures are often activated. To achieve a good vibration reduction effect covering all working conditions, one cannot directly use linear electromechanical coupling theory to instruct the piezoelectric damping design for nonlinear structures. Therefore, this paper defines the Nonlinear Modal Electromechanical Coupling Factor (NMEMCF) and proposes the corresponding numerical method for the first time to quantitatively evaluate the electromechanical coupling capability of nonlinear piezoelectric structures. Three candidate definitions of the NMEMCF are given, including two frequency definitions and one energy definition. The energy definition is the most promising one. It is not only applicable to both conservative and dissipative nonlinear structures but also compatible with the linear MEMCF. In addition, based on the energy formula, the NMEMCF can be obtained by only performing one nonlinear modal analysis in the open-circuit state. The analytical findings and the numerical tool are validated against two piezoelectric structures with different types of nonlinearities. A strong correlation among the NMEMCF, geometric parameters, and energy dissipation is observed. The results confirm that the proposed NMEMCF captures the physics of the electromechanical coupling phenomenon associated with nonlinear piezoelectric structures and can be used as an essential design indicator of piezoelectric damping, especially for variable working conditions.     

Efficient aerodynamic shape optimization using variable-fidelity surrogate models and multilevel computational grids

A variable-fidelity method can remarkably improve the efficiency of a design optimization based on a high-fidelity and expensive numerical simulation, with assistance of lower-fidelity and cheaper simulation(s). However, most existing works only incorporate “two” levels of fidelity, and thus efficiency improvement is very limited. In order to reduce the number of high-fidelity simulations as many as possible, there is a strong need to extend it to three or more fidelities. This article proposes a novel variable-fidelity optimization approach with application to aerodynamic design. Its key ingredient is the theory and algorithm of a Multi-level Hierarchical Kriging (MHK), which is referred to as a surrogate model that can incorporate simulation data with arbitrary levels of fidelity. The high-fidelity model is defined as a CFD simulation using a fine grid and the lower-fidelity models are defined as the same CFD model but with coarser grids, which are determined through a grid convergence study. First, sampling shapes are selected for each level of fidelity via technique of Design of Experiments (DoE). Then, CFD simulations are conducted and the output data of varying fidelity is used to build initial MHK models for objective (e.g. C_D) and constraint (e.g. C_L, C_m) functions. Next, new samples are selected through infill-sampling criteria and the surrogate models are repetitively updated until a global optimum is found. The proposed method is validated by analytical test cases and applied to aerodynamic shape optimization of a NACA0012 airfoil and an ONERA M6 wing in transonic flows. The results confirm that the proposed method can significantly improve the optimization efficiency and apparently outperforms the existing single-fidelity or two-level-fidelity method.     

Design of aircraft structures against threat of bird strikes

In this paper, a method to design bird-strike-resistant aircraft structures is presented and illustrated through examples. The focus is on bird strike experiments and simulations. The explicit finite element software PAM-CRASH is employed to conduct bird strike simulations, and a coupled Smooth Particles Hydrodynamic (SPH) and Finite Element (FE) method is used to simulate the interaction between a bird and a target structure. The SPH method is explained, and an SPH bird model is established. Constitutive models for various structural materials, such as aluminum alloys, composite materials, honeycomb, and foam materials that are used in aircraft structures, are presented, and model parameters are identified by conducting various material tests. Good agreements between simulation results and experimental data suggest that the numerical model is capable of predicting the dynamic responses of various aircraft structures under a bird strike, and numerical simulation can be used as a tool to design bird-strike-resistant aircraft structures.     

蜂窝夹芯旋转壳的屈曲分析

蜂窝夹芯复合材料旋转壳是航天器中的重要部件,且常常与其它部件互相联结,联结区局部高应力往往诱发复杂的局部屈曲模态,为此发展了一套有限元求解方案。针对蜂窝夹芯层合壳体构造了一种基于相对自由度技术的32节点三层壳元,这种单元避免了传统壳元的转动自由度,易与三维实体单元连接,使变厚度、带有补强的复合材料层合壳体等复杂结构得以正确建模。同时运用旋转周期结构有限元技术对大规模的空间蜂窝夹芯层合结构成功实施了屈曲分析。数值算例表明了计算策略的有效性和优越性。     

An alternative numerical approach for full scale characterisation for riveted joint design

This paper describes a pure numerical methodology (FE database) to improve the representativeness of joint equivalent models for airframe crashworthiness. This method is based on material constitutive models and failure criterion in accurate 3D FE simulations. The interest of FE simulations is to define the dynamic strength of many types of riveted joints with a reduced cost compared to a pure experimental way. The FE database method is carried out on elementary riveted joints to predict and to analyse: first the post-riveting initial strain and stress state, and second several kinds of joint failure (e.g., crack propagation, rivet shearing or pull-out). The aim of the first step is to start the mechanical strength simulations with a correct deformed shape and post-riveting state. The responses of 3D riveted joint simulations can then be considered as reference and be used to optimise the mechanical properties of equivalent joint elements. A new equivalent joint element is developed to improve the representativeness of an airframe crash simulation.     

失谐周期压电复合材料结构中的波动局部化研究

考虑力电耦合效应的影响,研究了层状失谐周期压电复合材料结构中的波动局部化问题。根据界面上力电连续条件,推导了结构中相邻单胞间的传递矩阵。以力场和电场变量为状态向量,给出了结构中局部化因子的表达式。作为算例,计算了结构中的波动局部化因子。计算结果表明,压电陶瓷的压电效应对周期压电复合材料的波动局部化特性有显著影响,压电常数越大局部化因子值越大,结构的局部化程度越强;结构的失谐度越大,频率通带区间内的局部化因子值越大,局部化程度越强。分析结果对于周期压电复合材料结构的优化设计和振动控制具有理论参考价值。     

Coupling fluid (CFD) and structural (FE) models using finite interpolation elements

In the field of aeroelasticity, interactions between elastic structures and fluid flow are investigated. Recently, numerical aeroelastic models have been built composing those of the combining fluid dynamics (CFD) and the computational structural dynamics (CSD) domains. Since the fluid and the structural models differ in their formulation and discretization, an interface model has to be introduced that represents the connectivity and physical interaction between the two single domain models. In the following, a scheme for coupling fluid (CFD) and structural models (FE) in space is presented which is based on finite interpolation elements. It is applied to static aeroelastic problems, in order to predict the equilibrium of elastic wing models in transonic fluid flow. The structure is represented by finite elements and the related equations are solved using commercial FE analysis codes. The transonic fluid flow is described by the three-dimensional Euler equations, solved by an upwind scheme procedure. The resulting coupled field problem containing the fluid and the structural state equations, is solved by applying a partitioned solution procedure. In each solution step the interface and boundary conditions are exchanged and updated. Here, a fixed-point iteration scheme is used for the coupled aeroelastic equations.     

脉冲子结构与有限元刚-弹混合连接的子结构方法

航天器结构的日益复杂庞大对系统级的动力学建模仿真以及进一步的结构优化提出了巨大挑战。为提高动力学求解效率,通常引用动态子结构方法。本文利用适于处理瞬态冲击问题的脉冲子结构(IBS)方法,并对其进行改进,将基于脉冲响应函数(IRF)的子结构与有限元建立的子结构综合,同时考虑刚性、弹性以及刚-弹混合连接情形下的子结构综合格式。通过3个数值算例,验证了方法的正确性。最后将刚-弹混合连接下的子结构方法应用到月球探测器软着陆的动态响应预测,结果表明该方法适于对月球探测器软着陆动态响应进行高精度快速预测,并且可以应用于月球探测器的局部动力学结构优化。     

A novel surrogate modeling strategy of the mechanical properties of 3D braided composites

In this paper, a surrogate-based modeling methodology is developed and presented to predict the elastic properties of three dimensional (3D) four-directional braided composites. Using this approach, the prediction process becomes feasible with only a limited number of training points. The surrogate models constructed using Finite Element (FE) method and Diffuse Approximation, reduce the computational time and cost for preparing experimental samples. In the FE model, multiscale method is applied to couple the computations of elastic properties at microscale and mesoscale. Subsequently, a parametric study is performed to analyze the effects of the three design parameters on the elastic properties. Satisfactory results are obtained via the surrogate-based modeling predictions, which are compared with the experimental measurements. Moreover, the predictions obtained from surrogate models concur well with the FE predictions. This study orients a new direction for predicting the mechanical properties based on surrogate models which can effectively reduce the sample preparation cost and computational efforts.     

Experimental,analytical and numerical investigation on tensile behavior of twisted fiber yarns

Stitched composite materials are emerging as a promising material due to their high interlaminar strength, combined performance and light weight. The mechanical properties of stitch yarns are very essential for stitched composite structures. In this study, the tensile behaviors of the twisted fiber yarn in stitched composites were investigated experimentally, analytically and numerically. Two kinds of cross-sectional area of twisted yarn are proposed and discussed. The paper presents an intersecting circle model to describe the cross-section of twisted fiber yarns, and a physics-based theoretical model to predict the effective tensile moduli. The numerical models take into account the cross-sectional characteristic and the twist architecture. The investigation shows that: the sum of each fiber area should be used for experimental analysis; and the cross-sectional area surrounded by the yarn profile should be used for theoretical predictions and finite element (FE) simulations. The relative errors of the prediction method and the FE simulation are less than 2% and 1%, respectively. The friction between the fibers is derived, and the effect of friction on mechanical properties is discussed. The investigation method will serve as a fundamental component of twisted fiber bundle/yarn analysis.     

基于拓扑优化的叶盘压电阻尼器性能研究

针对压电分支阻尼技术在航空发动机叶盘结构振动抑制问题中的应用,提出了一种拓扑优化方法,可给出限定用量的压电材料在轮盘上的最佳位置,以提升压电分支阻尼的上限。在轮盘上布置压电材料还可防止对叶片通道内流场的影响,避免降低流体效率。首先,论述了该拓扑优化方法的原理,推导了模态机电耦合系数这一核心参数的计算公式及其与最佳阻尼比、模态应变场的关系。其次,建立了基于模态应变场的压电材料分布拓扑优化方法,可用于任意有限元模型。最后,将该优化方法应用于一个航空发动机压气机叶盘结构模型,分别针对单一和多阶模态进行了压电材料在轮盘上分布的拓扑优化,研究这种铺设方式对各典型模态(轮盘主导、叶片主导、耦合振动)的振动抑制效果。结果表明,在仅采用占轮盘质量5%的压电材料的情况下,优化后的压电阻尼器最多可以为轮盘振动主导模态及叶片-轮盘耦合振动模态提供约13%的模态阻尼比,为叶片主导模态提供的模态阻尼比集中在2%～4%。