弯曲波对含多圆孔薄板的散射与动应力集中

采用复变函数法和多极坐标方法,研究了弯曲波对含有多圆孔薄板的散射问题。通过板的弯曲波动方程和内力方程的推导,求出在入射弯曲波条件下该问题的一般解的函数逼近序列和边界条件的表达式。用展开正交函数的方法将待解的问题归结为对一组无穷代数方程组的求解。最后,给出了含3圆孔薄板的孔边动应力集中系数的结果,并分析了孔间距和波数对动应力分布的影响。     

板中正交静应力与任意方向兰姆波波速的关系

运用部分波分析法(或子波分析)理论推导了在正交静应力下,板中任意方向的Lamb波的频散方程,给出了任意方向Lamb波波速、正交静应力和频率之间的关系. 进行了数值计算并讨论了对称方向的Lamb波波速对单向静应力的依赖关系,为将Lamb波用于应力测量提供了理论依据.      

串联电气设备支架隔震体系地震响应半解析法

通过并联橡胶隔震支座,建立串联高架电气设备支架隔震体系及力学模型,应用分布参数梁振动理论, 通过边界条件引入集中参数,推导其频率方程,并用数值方法求得频率及振型. 应用Betti定律,推导具有集中分布参数柔性节点的多节电瓷型高压电气设备的正交条件,得到广义质量及广义载荷. 将该串联隔震体系的非比例阻尼分解为比例阻尼部分和非比例阻尼部分,应用Hamilton原理推导出非比例阻尼部分等效振型阻尼比,实现串联电气设备支架隔震体系振动方程的解耦,然后通过振型叠加法求得结构的地震响应. 最后应用该半解析法与有限元数值积分法对一330kV电压互感器地震响应进行分析. 结果表明:该隔震体系能显著减小电气设备的地震响应,同时该半解析法求解的地震响应与有限元数值积分求解的结果相一致,说明该隔震体系的有效性与该半解析方法求解非比例阻尼串联电气设备支架隔震体系地震响应的正确性.      

高压电气设备支架串联体系地震响应半解析法

针对具有集中分布参数柔性节点的多节电瓷型高压电气设备支架串联体系,通过应用分布参数梁振动理论与集中参数通过边界条件的引入,推导出其频率方程,然后通过数值方法求得频率及振型。通过应用Betti定律,推导出具有集中分布参数柔性节点的多节电瓷型高压电气设备支架串联体系振型的正交条件,应用此正交条件对集中参数与分布参数的振动方程解耦,从而得到广义质量及广义荷载,然后通过振型叠加法求解结构的地震响应。最后应用该半解析法与有限元法数值积分法对一330 KV电压互感器支架串联体系地震响应进行分析,分析结果表明,该半解析法求解的地震响应与有限元直接积分法求解的地震响应相一致,从而为具有集中分布参数柔性节点的多节电瓷型高压电气设备支架串联体系的地震响应分析和抗震设计提供了一种新的途径。     

计算结构动力方程三次精度的直接积分法

以往计算结构动力方程的无条件稳定积分格式的证明,是在阻尼矩阵满足振型正交条件下得到的,文本给出的三次精度无条件稳定积分格式的证明,可不要求阻尼矩阵满足振型正交条件。此外本文提出的高精度方法和以往的高精度方法相比具有存贮空间小,计算量少的优点,本文方法还具有没超越现象的优点。     

采用精确化理论求解厚板任意形开孔动应力集中

摘要：本文基于复变函数与保角映射法,采用平板弯曲振动精确化方程[9],对含任意形开孔平板中弹性波散射与动应力集中问题进行了研究。利用正交函数展开的方法将待解的问题归结为对一组无穷代数方程组的求解。作为算例,计算了自由边界条件下圆孔和椭圆孔的动弯矩集中系数的数值结果,并对板厚与孔径比对动弯矩分布的影响做了分析研究。结果表明：入射波数、平板厚度和椭圆偏心率等参数对动弯矩的分布都有很大的影响。在较低频率和平板较薄的情况下,基于文献[9]的方程与基于Mindlin板的动弯矩结果在数值分布上是基本一致的;在较高频率和平板较厚的情况下,基于文献[9]的方程与基于Mindlin板的动弯矩结果在数值分布计算结果相差较大。由于文献[9]给出的平板振动精确化方程是在没有任何工程假设条件下得到的,因此本文的分析计算结果更精确一些。     

考虑可控性的压电作动器拓扑优化设计

压电作动器可以把电能转换成机械能,在结构主动振动控制中具有应用背景. 由于压电作动器的布局对振动控制效果影响很大,因此作动器布局优化一直是结构控制研究的关键之一. 为了提高压电结构控制能量的利用效率,本文提出了以提高结构可控性为目标的压电作动器的拓扑优化方法. 基于经典层合板理论对压电结构进行了有限元建模,并采用模态叠加法将动力控制方程映射到模态空间,推导了基于控制矩阵奇异值的可控性指标. 优化模型中,选取可控性指标指数形式为目标函数,将设计变量定义为作动器单元的相对密度,并基于人工密度惩罚模型构造了压电系数惩罚模型,给出了基于控制矩阵奇异值的可控性指标关于设计变量的灵敏度分析方法. 优化问题采用基于梯度的数学规划法求解. 数值算例验证了灵敏度分析方法和优化模型的有效性,并讨论了主要因素对优化结果的影响.      

四边固支蜂窝夹层板1:3内共振非线性动力学分析

基于经典叠层板理论和几何大变形理论,将铝基蜂窝芯层等效为一正交异性层,等效弹性参数由修正后的Gibson公式得出,对四边固支蜂窝夹层板非线性动力学特性进行了分析。考虑横向阻尼的影响,建立了四边固支蜂窝夹层板受横向激振力作用的受迫振动微分方程,通过振型正交化将蜂窝夹层板受迫振动微分方程简化成双模态下的动力学控制方程,利用同伦分析方法对双模态下蜂窝夹层板的动力学控制方程进行研究,得到了1:3内共振下的幅频特性曲线,研究了不同结构尺寸对动力学特性的影响以及蜂窝夹层板作稳态运动时的稳定性问题。本文得到的结果为蜂窝夹层板的设计和实际应用提供了理论依据和数值参考。     

考虑可控性的压电作动器拓扑优化设计

压电作动器可以把电能转换成机械能,在结构主动振动控制中具有应用背景.由于压电作动器的布局对振动控制效果影响很大,因此作动器布局优化一直是结构控制研究的关键之一.为了提高压电结构控制能量的利用效率,本文提出了以提高结构可控性为目标的压电作动器的拓扑优化方法.基于经典层合板理论对压电结构进行了有限元建模,并采用模态叠加法将动力控制方程映射到模态空间,推导了基于控制矩阵奇异值的可控性指标.优化模型中,选取可控性指标指数形式为目标函数,将设计变量定义为作动器单元的相对密度,并基于人工密度惩罚模型构造了压电系数惩罚模型,给出了基于控制矩阵奇异值的可控性指标关于设计变量的灵敏度分析方法.优化问题采用基于梯度的数学规划法求解.数值算例验证了灵敏度分析方法和优化模型的有效性,并讨论了主要因素对优化结果的影响.     

基于平板振动精确化方程求解动应力集中问题

基于文献[8]给出的平板弯曲振动精确化方程,对含圆孔平板中弹性波散射与动应力集中问题进行了研究.文中给出了分别基于Mindlin板与精确化板方程在不同参数下圆孔动弯矩集中系数的数值结果,并对结果进行了对比分析和讨论.结果表明:在较低频率和薄板情况下,基于文献[8]的方程与基于Mindlin板理论得到的动弯矩结果是基本一致的;在较高频率和厚板情况下,基于文献[8]的方程与基于Mindlin板理论的动弯矩结果相差较大,最大值超出可达16%.由于文献[8]给出的平板振动精确化方程是在没有任何工程假设条件下得到的,因此其分析计算结果更精确一些.      

LINEAR ACTIVE STRUCTURES AND MODES (Ⅱ)--DISCRETE SYSTEMS AND BEAMS

The basic concepts about the active structures and some attributes of the modes were presented in paper "Liner Active Structures and Modes ( Ⅰ ) ". The characteristics of the active discrete systems and active beams were discussed, especially, the stability of the active structures and the orthogonality of the eigenvectors. The notes about modes were portrayed by a model of a seven-storeyed building with sensors and actuators. The concept of the adjoint active structure was extended from the discrete systems to the beams that were the representations of the continuous structures. Two types of beams with different placements of the measuring and actuating systems were discussed in detail. One is the beam with the discrete sensors and actuators, and the other is the beam with distributed sensor and actuator function. The orthogonality conditions were derived with the modal shapes of the active beam and its adjoint active beam. An example shows that the variation of eigenvalues with feedback amplitude for the homo-configuration and non-homo-configuration active structures.     

LINEAR ACTIVE STRUCTURES AND MODES （ Ⅱ ） — DISCRETE SYSTEMS AND BEAMS

The basic concepts about the active structures and some attributes of the modes were presented in paper “Liner Active Structures and Modes ( Ⅰ ) “. The characteristics of the active discrete systems and active beams were discussed, especially, the stability of the active structures and the orthogonality of the eigenvectors. The notes about modes were portrayed by a model of a seven-storeyed building with sensors and actuators. The concept of the adjoint active structure was extended from the discrete systems to the beams that were the representations of the continuous structures. Two types of beams with different placements of the measuring and actuating systems were discussed in detail. One is the beam with the discrete sensors and actuators, and the other is the beam with distributed sensor and actuator function. The orthogonality conditions were derived with the modal shapes of the active beam and its adjoint active beam. An example shows that the variation of eigenvalues with feedback amplitude for the homo-configuration and non-homo-configuration active structures.     

Remarks on the numerical solution of the adjoint quasi‐one‐dimensional Euler equations

We examine the numerical solution of the adjoint quasi‐one‐dimensional Euler equations with a central‐difference finite volume scheme with Jameson‐Schmidt‐Turkel (JST) dissipation, for both the continuous and discrete approaches. First, the complete formulations and discretization of the quasi‐one‐dimensional Euler equations and the continuous adjoint equation and its counterpart, the discrete adjoint equation, are reviewed. The differences between the continuous and discrete boundary conditions are also explored. Second, numerical testing is carried out on a symmetric converging–diverging duct under subsonic flow conditions. This analysis reveals that the discrete adjoint scheme, while being manifestly less accurate than the continuous approach, gives nevertheless more accurate flow sensitivities. Copyright © 2011 John Wiley & Sons, Ltd.     

流固耦合的周期加强板的振动及声辐射研究

研究了流体负载下的无穷大双周期加强板, 在周期谐振力作用下的振动响应和声辐射,并提出了一种基于有限元和空间波数法的半解析半数值方法. 首先利用有限元的方法对周期结构进行单元离散, 并将结构对薄板的作用力等效为节点力的作用. 然后通过周期结构的振动方程, 结合薄板与结构的位移边界条件, 建立了节点力与薄板节点位移的函数方程. 最后应用空间波数法和傅里叶变换, 并采用数值计算的方法求解出薄板的节点位移, 得到了周期加强板关于离散节点位移的振动和辐射声压方程. 在数值算例中, 对该方法的正确性进行了验证, 并且分析了周期结构对薄板的振动和声辐射的影响.     

Frequency equations for the in-plane vibration of orthotropic circular annular plate

Orthotropic circular annular plates have a lot of applications in engineering such as space structures and rotary machines. In this paper, frequency equations for the in-plane vibration of the orthotropic circular annular plate for general boundary conditions were derived. To obtain the frequency equation, first the equation of motion for the circular annular plate in the cylindrical coordinate is derived by using the stress-strain- displacement expressions. Helmholtz decomposition is used to uncouple the equations of motion. The wave equation is obtained by assumption a harmonic solution for the uncoupled equations. Using the separation of the variables leads to the general wave equation solution and the in-plane displacements in the r and θ directions. Finally, boundary conditions are exerted and the natural frequency is derived for general boundary conditions. The obtained results are validated by comparing with the previously reported and those from finite element analysis.     

Design optimization of unsteady airfoils with continuous adjoint method

In this paper, a new unsteady aerodynamic design method is presented based on the Navier-Stokes equations and a continuous adjoint approach. A basic framework of time-accurate unsteady airfoil optimization which adopts time-averaged aerodynamic coefficients as objective functions is presented. The time-accurate continuous adjoint equation and its boundary conditions are derived. The flow field and the adjoint equation are simulated numerically by the finite volume method (FVM). Feasibility and accuracy of the approach are perfectly validated by the design optimization results of the plunging NACA0012 airfoil.     

Algorithm for transient growth of perturbations in channel Poiseuille flow

This study develops a direct optimal growth algorithm for three-dimensional transient growth analysis of perturbations in channel flows which are globally stable but locally unstable. Different from traditional non-modal methods based on the OrrSommerfeld and Squire(OSS) equations that assume simple base flows, this algorithm can be applied to arbitrarily complex base flows. In the proposed algorithm, a reorthogonalization Arnoldi method is used to improve orthogonality of the orthogonal basis of the Krylov subspace generated by solving the linearized forward and adjoint Navier-Stokes(N-S) equations. The linearized adjoint N-S equations with the specific boundary conditions for the channel are derived, and a new convergence criterion is proposed. The algorithm is then applied to a one-dimensional base flow(the plane Poiseuille flow) and a two-dimensional base flow(the plane Poiseuille flow with a low-speed streak)in a channel. For one-dimensional cases, the effects of the spanwise width of the channel and the Reynolds number on the transient growth of perturbations are studied. For two-dimensional cases, the effect of strength of initial low-speed streak is discussed. The presence of the streak in the plane Poiseuille flow leads to a larger and quicker growth of the perturbations than that in the one-dimensional case. For both cases, the results show that an optimal flow field leading to the largest growth of perturbations is characterized by high-and low-speed streaks and the corresponding streamwise vortical structures.The lift-up mechanism that induces the transient growth of perturbations is discussed.The performance of the re-orthogonalization Arnoldi technique in the algorithm for both one-and two-dimensional base flows is demonstrated, and the algorithm is validated by comparing the results with those obtained from the OSS equations method and the crosscheck method.     

First-order perturbation solution to the complex eigenvalues

The matrix perturbation method is extended to discrete linear nonconservative system with unsymmetrical matrices in this article. By introducing the concept of the adjoint complex eigenvector and by making use of the orthogonality relationship in the complex mode theory, the first-order perturbation solution to the complex eigenvalues is derived. Numerical example shows that this method is efficient and practicable.     

A Discrete Adjoint Approach for the Optimization of Unsteady Turbulent Flows

In this paper we present a discrete adjoint approach for the optimization of unsteady, turbulent flows. While discrete adjoint methods usually rely on the use of the reverse mode of Automatic Differentiation (AD), which is difficult to apply to complex unsteady problems, our approach is based on the discrete adjoint equation directly and can be implemented efficiently with the use of a sparse forward mode of AD. We demonstrate the approach on the basis of a parallel, multigrid flow solver that incorporates various turbulence models. Due to grid deformation routines also shape optimization problems can be handled. We consider the relevant aspects, in particular the efficient generation of the discrete adjoint equation and the parallel implementation of a multigrid method for the adjoint, which is derived from the multigrid scheme of the flow solver. Numerical results show the efficiency of the approach for a shape optimization problem involving a three dimensional Large Eddy Simulation (LES).     

Fracture analysis of mode III crack problems for the piezoelectric bimorph

In this paper, a symplectic method based on the Hamiltonian system is proposed to analyze the interfacial fracture in the piezoelectric bimorph under anti-plane deformation. A set of Hamiltonian governing equations is derived from the Hamiltonian function by introducing dual variables of generalized displacements and stresses which can be expanded in series in terms of the symplectic eigensolutions. With the aid of the adjoint symplectic orthogonality, coefficients of the series are determined by the boundary conditions along the crack faces and along the external geometry. The stress\electric displacement intensity factors and energy release rates (G) directly relate to the first few terms of the nonzero eigenvalue solutions. The two ideal crack boundary conditions, namely the electrically impermeable and permeable crack assumptions, are considered. Numerical examples including the complex mixed boundary conditions are considered to show fracture behaviors of the interface crack and discuss the influencing factors.