两自由度体系颤振机理研究论文

纤维增强树脂基复合材料层合结构具有比强度高、比刚度大、阻尼特性好、疲劳寿命长、结构可设计性强等优点,在航空、航天及一些特殊领域中被广泛使用。然而,复合材料的各向异性,非均匀性等特点给复合材料结构的力学分析带来了一系列的挑战。尤其在航空航天领域,飞行器在运行过程中所处的环境和所受的载荷都非常复杂。除了考虑飞行器在这些复杂环境下的自振特性和确定性外载作用下的动力响应外,考虑随机性外载的影响也不容忽视。随机振动理论和方法就是处理这类问题的先进思想和重要手段,但在国内外航空航天领域中还很少实际应用,主要原因之一就是现有随机振动分析方法复杂而且低效,这在很大程度上限制了飞行器设计水平的提高。虚拟激励法是高效精确的随机振动分析方法,迄今已经在大跨度结构抗震、抗风,海洋平台和汽车随机振动等多个工程领域被数以百计的专家针对各工程领域的特点予以发展而取得很多实际成效。但是迄今为止,这一有力的工具却并未在航空航天领域被充分认识和应用,在这些具有战略意义的重要领域中,所应用的随机振动分析方法依然复杂低效,缺乏创新意识。本论文针对这一现状,依据航空航天领域材料和结构的复杂性,以及飞行器所处环境的复杂性,将虚拟激励法作了有针对性的发展,以完全自主版权的DDJ有限元程序系统为开发平台,完成了求解复合材料结构随机振动的高效精确分析程序。本论文中,着重对如下问题进行了研究:1.建立了基于Mindlin一阶剪切变形理论的复合材料层合板有限元分析模型,推导了层合板的有限元列式,在DDJ程序平台上对复合材料层合板的自振频率和模态进行了分析。将虚拟激励法引入到航空航天领域广泛使用的复合材料层合结构的随机振动分析中,针对复杂的复合材料结构有限元模型和非经典阻尼体系,发展了包含全部参振振型和随机激励点之间耦合项的随机振动高效求解方法,比较圆满地解决了传统计算方法精度差、效率低的应用障碍。2.本文推广虚拟激励法于敷设粘弹性阻尼层的复合材料层合结构的平稳和非平稳随机振动分析,建立了高效精确计算方法。尤其是综合考虑了粘弹性阻尼材料的性能参数随频率变化的特点以及复合材料层合结构本身的模态阻尼,建立了组合系统的非经典阻尼表达。为了解决随频率变化的非经典阻尼体系的平稳/非平稳随机响应,本文结合精细积分方法提出了一种直接解法,只需用原系统的实模态对虚拟激励法做出相应的发展,就可精确地求解频变阻尼系统的随机振动。据此对飞机水平尾翼的复合材料安定面结构进行了模拟研究,从精细的计算模型及合理的计算结果可以看出,本文所提出的方法对于这类相当复杂的复合材料结构的随机振动分析十分有效。3.研究飞机对大气紊流响应的主要方法是随机振动功率谱法。用高效、精确的分析方法计算不同飞行环境下飞机的响应,以预测飞机疲劳寿命和可靠度等是航空工程领域研究热点。本文在考虑了二维平面流中简谐振动平板产生的非定常力基础上,又按照虚拟激励法的特点同时考虑了竖向简谐风的影响,进而研究了复合材料二维机翼的大气紊流响应。随机激励谱选用了Dryden紊流频谱模型。结果表明,在处理二维机翼在大气紊流响应的随机问题中,基于简谐响应分析的虚拟激励法不但是精确算法,而且效率非常高,具有很大的实用优势。发展这一方法对于该领域的数值计算是很有价值的。4.计算流体动力学(CFD)是研究流体动力学的有力工具。本文为计算机翼颤振/抖阵分析中的气动参数,首次使用雷诺平均湍流模型对二维翼型截面的颤振导数进行了求解。基于D.K.Sun等最新提出的CFD网格控制算法以及所建立的数值风洞,计算了结构简谐运动下的气动力,并识别了湍流场中NACA0012翼型的颤振导数。将由此得到的颤振导数和气动力应用到大气紊流引起的随机振动计算中,并将计算结果与基于Theodorsen函数得出的响应解析解进行比较,得到了相当满意的一致。本文计算的CFD气动参数充分考虑了气体的分子粘性和紊流粘性,其作用相当于附加阻尼,因此比Theodosen函数方法限制更少、应用范围更广,而且在此基础上还可以考虑三维流和可压缩性。因此本文实施的基于CFD的气动力计算方法具有广阔的应用前景,将成为应用虚拟激励法于航空航天结构时确定气动参数的有力工具。可以说,这一成功的尝试为随机振动方法更广泛地应用于航空航天工程走出了很重要的一步。

机械振动是指物体或质点在其平衡位置附近所作有规律的往复运动。振动的强弱用振动量来衡量，振动量可以是振动体的位移、速度或加速度。

本文从古典耦合颤振理论、分离流颤振模型和三维桥梁额振分析等三个方面简要回顾了空气动力作用下大跨度桥梁风振稳定性研究的历史，比较全面地综述了桥梁额振稳定性理论由简单到复杂，由解析方法到数值方法、由二维颤振到三维额振以及由多模态参与到全模态参与的发展过程。为了便于定量地比较这几种颤振分析理论和方法的适宜条件和精度，以完全流线形的悬臂机翼和钝体截面的上海南浦大桥为例，计算和分析了颤振临界风速的数值结果。关键词：空气动力学 大跨度桥梁 颤振 稳定性 临界风速一、前言浸没在气流中的任一物体，都会受到气流的作用，这种作用通常称为空气力作用。当气流绕过一般为非流线形（钝体）截面的桥梁结构时，会产生涡旋和流动的分离，形成复杂的空气作用力[1].当桥梁结构的刚度较大时，结构保持静止不动，这种空气力的作用只相当于静力作用；当桥梁结构的刚度较小时，结构振动得到激发，这时空气力不仅具有静力作用，而且具有动力作用[2].风的动力作用激发了桥梁风致振动，而振动起来的桥梁结构又反过来影响空气的流场，改变空气作用力，形成了风与结构的相互作用机制。当空气力受结构振动的影响较小时，空气作用力作为一种强迫力，引起结构的强迫振动；当空气力受结构振动的影响较大时，受振动结构反馈制约的空气作用力，主要表现为一种自激力，导致桥梁结构的自激振动。当空气的流动速度影响或改变了不同自由度运动之间的振幅及相位关系，使得桥梁结构能够在流动的气流中不断汲取能量，而该能量又大于结构阻尼所耗散的能量，这种形式的发散性自激振动称为桥梁颤振[3].桥梁颤振物理关系复杂，振动机理深奥，因而桥梁颤振稳定性研究也经历了由古典耦合颤振理论到分离流颤振机理再到三维桥梁颤振分析的发展过程。二、古典耦合颤振理论尽管由气动弹性影响所引起的机翼动力失稳现象早在人类实现空中飞行梦想的最初年代里已经观察到了，但是非定常机翼颤振理论直到20年代初才取得了实质性的进展。1. Theodorsen机翼颤振理论1922年，bimbaum利用Prandtl的约束涡旋理论，提出了第一个简谐振动平板机翼的气动升力解析表达式。此后 Theodorsen，Wagner，Glanert，Kussner，Duncan和 Collar等气动专家对二维振动平板的非定常气动力表达式进行了10多年的深入研究[4，5]，直到 1935年，才由Theodorsen用势能原理第一次求出了这一问题最完整的解答2.Bleich悬索桥颤振分析1940年秋天，美国华盛顿州Tacoma悬索桥风毁失事，人们很自然地将这一风振现象比拟为裹冰状态输电缆的驰振或平板机翼的颤振。Bleich试图用Theodorsen平板机翼颤振理论来解释这一事故，但是他发现居此计算得到的颤振临界风速远高于Tacoma悬索桥破坏当天的实际风速。显然机翼颤振系数不能直接用于气动现象更加复杂的钝体截面中，例如Tacoma悬索桥的桁架加劲梁断面。为此，BIeich又尝试用考虑桥面断面两边涡旋影响的附加升力项来修正Theodorsen气动力表达式，并通过逐次逼近方法计算出了较为合理的悬索桥颤振临界风速，从而建立起了悬索桥古典耦合颤振的分析方法［7］。3. Kloppel/Thiele诺模图1961年，Kltw和Thiele将BIeich悬索桥古典耦合颤振理论的逐次逼近过程编制成计算程序，引入无量纲参数，分别绘制出不同阻尼比条件下颤振方程实部和虚部为零的两条曲线的诺模图，利用诺模图可以直接求出颤振临界风速[8].该方法一直沿用到现在，例如ECCS中的附表[9].4.Van der Put计算公式1976年，van der put在Kloppel和Thiele诺模图的基础上，偏于安全地忽略了阻尼的影响，认为折算风速U/ωB和扭弯频率比ε=ωα/ωh之间具有近似线性关系，从而导出了平板古典耦合颤振临界风速Ucr的实用计算公式。三、分离流颤振机理当气流绕过振动着的非流线性截面时，在迎风面的棱角处气流将发生分离，同时产生涡旋脱落，也可能发生再附，其流态十分复杂，简单地采用Theodorsen表达式已经不能描述气流作用在非流线体上的非定常空气力[11].1.非定常气动力实验测量Theodorsen机翼气动力表达式是建立在有势流沿着翼面流动基础之上的。一旦气流有分离时，这一假定立即失效，而流动分离所引起的失速颤振现象最早是在螺旋浆和航空发动机叶片上观察到的。由于建立在分离流基础之上的非定常气动力表达式无法找到，因此从30年代开始，人们将注意力转向用实验方法来确定非定常气动力，主要通过两种方法来实施。一是直接测量法，即对确定形式振动的物体，采用拾振器、应变计或其他仪器直接测量气动力分量；二是间接测量法，即间接地从振动的物体上计算气动力的大小，这两种方法同样适用于机翼和桥梁断面。1958年，Forshing采用直接测量法测得了各种棱柱体的非定常气动力［12］，而Ukeguchi等人将 Halfman测量机翼非定常气动力的方法，首次用到了桥梁断面非定常气动力的测定中，他们采用机械方法在两个自由度方向用不同频率的简谐波激发刚体桥梁节段模型振动，在模型的支承处测量气动力[13].随后这种强迫振动技术在日本得到了很大的发展，被广泛地用来测定钝体断面的气动力和非线性性能[14-16].近年来，用于高速电子压力扫描阀技术的发展，使得多点同步测量得以实现，这项技术的应用开辟了非定常气动力测量的又一新途径「17」。与直接测量法相反，间接气动力测量方法一般只需要比较简单的实验设备，但是对实验的要求更高，这一方法在桥梁气动力学中的应用是由Scanlan首创的[l1][18]，很快在世界范围得到普及［19］［20］。2.非定常气动力计算模型桥梁结构分离流颤振实验加理论方法的建立与完善是与著名气动力专家R.H.Scanlan的贡献紧密联系在一起的。1967年，Scanlan首先提出对Theodorsen机翼气动力表达式进行修正的建议「11」。Scanlan认为，对于非流线性的钝体截面，不可能从基本的流体力学原理推导出类似于Theodorsen函数的气动函数，但可以通过专门设计的节段模型风洞实验测定小振幅条件下的气动力参数——颤振导数（Flutter Derivatives）来建立线性非定常气动力计算模型［18］1974年，Scanlan利用节段模型风洞实验中实测的颤振导数反算出过渡函数（Indicial Function）[21]，并与Theodorsen函数进行了比较，结果发现两种函数曲线相差很大，从而找到了利用古典耦合颤振理论分析钝体桥梁颤振问题所造成的误差原因。Scanlan还断言，从理论上找到适合于各种非流线型断面的过渡函数是不可能的。3.M维颤振分析方法一旦建立了非定常气动力计算模型，气动失稳临界状态就很容易确定了，其中，最典型的方法就是将所谓阿'片条理？quot；应用于气流与结构相互作用之中，确定出一个垂直于桥轴线方向的二维节段，假定沿着桥轴线方向的任意三维影响都可以忽略不计，由此可得二维颤振方程。与传统的机翼颤振类似，阻力方向的振动影响一般忽略不计。此外，还假定二维节段在h和α两个方向的振动是小振幅的同频简谐振动，这样就可以在传统的颤振分析中采用随折算频率变化的非定常气动力。4.Selberg计算公式在电子计算机诞生之前，颤振分析的数值计算工作一直是一项枯燥繁重的劳动。为了简化这项枯燥的工作，人们提出了许多颤振简化计算方法。对于平板机翼，由于每一种翼型的气动力表达式都有一定的差异，因此需要对一系列的结构参数进行分析，才能有针对性地进行额振计算，Theodorsen和Garrick[24]在这方面作了大量细致的工作，找到了一些实用的机翼颤振计算公式。而在桥梁结构方面，许多研究人员也作出了相同的努力，其中Selberg实用计算公式是被引用得最多的一种二维颤振实用计算公式。四、三线桥梁颤振分析Scanlan提出的非定常气动力计算模型较好地解决了非流线形截面的非定常气动力描述问题，其中二维颤振分析最为简单实用。但是随着桥梁跨径的日益增大，结构刚度急剧下降，特别是侧向刚度的下降，导致了侧弯与扭转振型紧密耦合。此外，结构各阶自振频率的差异很小，两个或两个以上振型参予颤振的可能性逐渐增加，因此，为了提高桥梁颤振分析精度，有必要寻求更精确的三维桥梁额振分析方法。1.时域分析法尽管桥梁颤振分析一般是在频域内进行的，但是也出现了一些时域分析方法。早在70年代初，Scanlan，Beliveau和Budlong采用飞行器设计中的传递函数首先提出了全时域分析方法[21]，Bucher和Lin将这种方法推广到了耦合模态颤振[27].但是，这一方法的主要困难在于寻找与实验所确定的气动导数相对应的合适的过渡函数，特别是当截面为非流线型时，难度更大。近年来，人们之所以投入了大量的精力从事开发有效的非定常气动力的时域表达式，主要是因为这种时域表达式既可与有限元结构计算模型相结合又能包含几乎所有的非线性因素，而这些非线性因素以前是一概忽略的。时域方法的发展是与诸如日本 Akashi桥、丹麦Storebraelt桥和意大利的Messina桥等超大跨度桥梁的Miyata等人清楚地阐明了时域分析方法在桥梁抖振响应估计中的优越性，特别是在采用有限元结构计算模型时的优势[28]，他们在片条假定的前提下，采用了传统的准定常气动力表达式。 Kovacs等人也曾提出过类似的方法[29].但在另一方面，Diana等人应用不同折算风速下的气动力系数等效线性化方法建立了一种所谓精确的准定常理论「30」，这一理论方法除了不能考虑气动力时效影响和升力的展向相关性之外，在许多方面被证明是足够数确的。另一种自激力模型是采用与Laplace变换相对应的有理函数来近似表示非定常气动力。实质上，这种思路与过渡函数是完全类似的。Xie等人将这一思想完善成状态空间法用来分析多模态三维桥梁颤振问题「31，32」。类似的方法还曾经由Lin和Li「31」，M.S.Li「34」，Boonyapinio[35]， Fujino[36]等人提出。2.多模态耦合颤振三维桥梁颤振分析更多地是采用频域分析方法。放弃片条假定后的三维桥梁颤振分析方法的应用还只有很短的历史，这种分析主要通过两种不同的途径来实现：第一条途径是将频率或时域内的非定常气动力直接作用到结构的三维有限元计算模型上，一般称为直接方法；第二条途径是把结构响应看作是分散在各阶模态上的影响，然后将各阶模态所对应的响应叠加起来，称为模态叠加法。直接法是由Miyata和Yamada提出的，他们把直接法归纳为用频域内气动导数所表示的一个复特征值问题「37」，这一方法的基本原理简单，但主要缺陷在于需要大容量的计算机来求解费时的复特征值问题。因此，许多研究者提出了另一种方法，即模态叠加法，现有许多种频域内的多模态参予颤振分析方法。Agar[38，39]和Chen[40，41]采用模态计术来求解线性二次特征值方程。作为机翼颤振分析方法一p-k的推广，Nmini等人「42」和程「43」提出了更加一般性的p-k-F法，通过求解模态方为确定颤振前后的状态。更进一步的还有Lin和Yang[44]，Jones和Scanlan[45]，Tanaka等人[46]，Jain等人[47]直接利用行列式搜索法求解广义特征矩阵的复特征值。几乎所有三维颤振分析都是在频域中进行的，并且基于了模态叠加假定。这一假定认为固有模态之间的动力耦合是通过自激气动力来实现的。然而，值得注意的是，这一假定存在着一些本质上的缺陷。首先多少阶模态和那些模态参与了颤振失稳，特别是在结构跨径很大或在施工过程中结构总体刚度尚未完全达到时，极有可能发生有两个以上的振动模态参与了颤振；其次，这种模态组合仅仅是颤振模态的一种近似表达式，没有任何理由使人们相信这是完全精确的，特别是在参与颤振的模态之间缺乏几何相似性时，颤振模态本身会变得非常复杂。正是考虑到这些因素，有必要建立一种更加综合和精确的方法来分析颤振模态，增进对悬吊体系桥梁气动失稳机理的理解和认识。3.全模态颤振分析全模态颤振分析方法是由本文作者提出的一种适合于大跨度桥梁颤振计算的方法，它是在Scanlan非定常气动力假定基础上建立起来的一种频域内颤振分析的精确方法，是对多模态颤振分析的一种推广[48，49].所谓精确方法，主要体现在两个方面，首先全模态方法不再像多模态方法那样将自激气动力作为外力作用在桥梁结构上，而是把桥梁结构与绕流气体作为一个相互作用的整体系统，建立系统颤振方法。五、算例比较为了比较各种桥梁颤振分析方法的适用性和精确性，现以流线型断面的悬臂机翼结构和钝体截面的上海南浦大桥为例，分析和比较颤振临界风速的计算结果。1.悬臂机翼结构第一个算例涉及到具有流线型断面的悬臂机翼结构，主要考虑到该结构具有精确的非定常气动力表达式，因而可以求得颤振临界风速的解析解——Theodorsen解，表1列出了分别按照六种不同方法计算得到的颤振临界风速及其与Theodorsen解的相对误差，这六种方法包括Theodorsen解，古典耦合颤振的van der Put实用计算公式，分离流二维颤振的Selberg实用计算公式，代表三维颤振时域分析方法的状态空间法，代表三维颤振多模态参与频域方法的p-k-F法，以及本文作者提出的全模态颤振分析法［49］。由于悬臂机翼颤振是一种古典耦合颤振，因此采用Van der Put计算公式精度较高，而采用分离流颤振假定的Selberg计算公式则误差较大，在三种数值计算方法中，全模态颤振分析方法的计算精度最高。2.上海南浦大桥第二个算例为上海南浦大桥，该桥是带有双I字梁钝体截面的结合梁斜拉桥，因此，不存在类似于机翼颤振的精确解。列出了分别按照其余5种方法的计算结果和弯扭两个模态耦合颤振的计算结果及其与全模态分析结果的相对误差。Van der Put计算公式和Selberg计算公式均不能用于钝体截面的桥梁颤振计算，而其余四种数值分析方法的临界风速计算结果均随参与颤振的振型数量的增加而增大。六、结语大跨度桥梁颤振稳定性分析始于Theodorsen的古典耦合颤振理论，Scanlan结合非流线型的桥梁断面提出了分离流非定常气动力表达式及其相应的分离流颤振理论，在此基础上，逐步形成和完善了二维和三维桥梁颤振分析方法。因此，空气动力作用下大跨度桥梁风振稳定性研究经历了由简单到复杂、由解析方法到数值方法、由二维桥梁颤振分析到三维桥梁颤振分析以及由多模态参与颤振到全模态参与颤振的发展过程。更多关于工程/服务/采购类的标书代写制作，提升中标率，您可以点击底部官网客服免费咨询：