任意形状功能梯度板的自由振动分析

OpenJournalofAcousticsandVibration声学与振动,2019,7(1),28-40PublishedOnlineMarch2019inHans.://doi.org/10.12677/ojav.2019.71004文章引用:玉云艳,朱翔,李天匀,郭文杰.任意形状功能梯度板的自由振动分析[J].声学与振动,2019,7(1):28-40.DOI:10.12677/ojav.2019.71004FreeVibrationAnalysisofFunctionalGradientMindlinPlateofArbitraryShapeYunyanYu1,2,3*,XiangZhu1,2,3,TianyunLi1,2,3,WenjieGuo1,4#1SchoolofNavalArchitectureandOceanEngineering,HuazhongUniversityofScienceandTechnology,WuhanHubei2HubeiKeyLaboratoryofNavalArchitectureandOceanEngineeringHydrodynamics,WuhanHubei3CollaborativeInnovationCenterforAdvancedShipandDeep-SeaExploration,Shanghai4EngineeringResearchCenterofRailwayEnvironmentVibrationandNoiseMinistryEducation,EastChinaJiaotongUniversity,NanchangJiangxiReceived:Feb.16th,2019;accepted:Mar.6th,2019;published:Mar.13th,2019AbstractInthispaper,thefreevibrationcharacteristicsoffunctionallygradedplatesinarbitraryshapewerestudiedbasedontheimprovedRayleighRitzmethod.Assumingthatthematerialofthefunctionallygradedplatechangesexponentiallyalongthethicknessdirection,thevibrationoftheplateisdescribedbyMindlinplatetheoryandtheimprovedFourierseriesisusedasthedis-placementtolerancefunction.Theenergyfunctionalexpressionsofthefunctionallygradedplatearederivedandthenaturalfrequencyisobtainedbysolvingthem.Theconvergenceandaccuracyoftheproposedmethodareverifiedbycomparedwiththeexistingliterature.Thentheproposedmethodwasappliedtothetriangularandcircularplatesrespectively,inwhichtheversatilityoftheproposedmethodwasshowing.Finallytheinfluencesofgradientindex,boundaryconditionandthicknessofthefunctionallygradedplateonthefreevibrationarediscussed.KeywordsArbitraryShape,FunctionallyGraded,TypicalBoundaryConditions,Rayleigh-RitzMethod,ImprovedFourierSeriesMethod任意形状功能梯度板的自由振动分析玉云艳1,2,3*，朱翔1,2,3，李天匀1,2,3，郭文杰1,4#1华中科技大学船舶与海洋工程学院，湖北武汉2船舶与海洋水动力湖北省重点实验室，湖北武汉3高新船舶与深海开发装备协同创新中心，上海*第一作者。#通讯作者。玉云艳等DOI:10.12677/ojav.2019.7100429声学与振动4华东交通大学，铁路环境振动与噪声教育部工程研究中心，江西南昌收稿日期：2019年2月16日；录用日期：2019年3月6日；发布日期：2019年3月13日摘要本文基于改进的Rayleigh-Ritz法对任意形状功能梯度材料板的自由振动特性进行了分析。假设功能梯度板的材料沿厚度方向指数变化，用Mindlin板理论描述板的振动。引入改进的Rayleigh-Ritz法，利用改进的傅里叶级数作为位移容许函数，结合Mindlin板理论推导得到功能梯度板的能量泛函表达式，根据能量泛函变分原理得到了振动系统的特征方程，求解得到功能梯度板的振动固有频率。通过与文献中功能梯度板的固有频率值对比，验证了本文方法的收敛性与准确性。然后分别对三角形板和圆形板的固有频率进行了求解，表明本文方法的通用性。然后通过算例分析讨论了功能梯度板的梯度指数、边界条件及板厚等参数对固有频率的影响。关键词任意形状，功能梯度，典型边界条件，Rayleigh-Ritz法，改进的傅里叶级数法Copyright©2019byauthor(s)andHansPublishersInc.ThisworkislicensedundertheCreativeCommonsAttributionInternationalLicense(CCBY).引言日本材料学家新野正之，平井敏雄和渡边龙三等人在20世纪80年代中后期提出了功能梯度材料的概念[1][2]。功能梯度材料选择两种不同性能的材料，根据使用要求采用材料复合技术，使中间的组成呈现连续的梯度变化而内部不存在明显的界面，材料的性质和功能沿厚度方向呈梯度变化。这样就可以从材料的本身改善材料的各种性能。由于功能梯度材料具有良好的力学性能，日益成为研究者研究的重点[3]。针对功能梯度材料板壳的振动问题，已经有大量文献开展了研究。王明禄[4]考虑了由材料的非均匀性引起的中面应变，求解了固有频率。曹志远[5]利用特殊类型各项异性及薄膜弯曲耦联的二维常系数方程研究功能梯度材料圆柱壳的固有频率的解析解。HenryKhov[6]拓展了选用的位移函数，使其能够精准的满足边界点的控制微分方程。徐坤[7]假设剪应力沿厚度方向呈抛物线分布，利用剪切余能与剪切应变能相等原理与虚位移原理及Navier方法研究平板的振动。梁斌[8]利用Rayleigh-Ritz方法建立了功能梯度圆柱壳自由振动固有频率的特征方程，其中材料参数为指数型体积分数，推导了一端固定一端自由和两端简支两种基本边界情况下的圆柱壳的自振特性。李伟[9]基于哈密顿原理，采用了动力变分方程，并把任意四边边界条件转化为统一的振型参数，对正交各向异性功能梯度开口圆柱壳的固有频率与模态进行求解。陈淑萍[10]基于一阶剪切变形理论，求解了两端简支Timoshenko梁自由振动的固有频率，并研究了材料梯度指数对梁动力响应的影响。尹硕辉[9]基于一阶剪切变形理论，应用等几何有限元法对功能梯度板进行研究。此外还有大量的数值方法已经得到应用[12][13][14]。以上研究多是应用解析方法对单一形状的功能梯度板壳的研究，或者是应用数值方法，对于应用解析方法对任意形状功能梯度板壳问题的研究成果还非常少。OpenAccess玉云艳等DOI:10.12677/ojav.2019.7100430声学与振动本文采用改进的Rayleigh-Ritz法，基于一阶剪切理论对功能梯度材料板的自由振动进行研究。考虑材料沿厚度方向梯度指数变化的功能梯度板，引入改进的Rayleigh-Ritz法，利用改进的傅里叶级数和Mindlin板理论推导得到功能梯度板的能量泛函表达式，求解得到功能梯度板的振动固有频率。通过数值算例的计算并和文献对比证明本方法的收敛性与准确性，然后对其他形状板的自由振动进行了求解。2.理论模型2.1.功能梯度材料对于功能梯度材料板结构，通常研究其材料属性沿厚度方向的变化。采用空间连续函数模型，则由两种材料组成的功能梯度材料沿厚度方向材料组分的分布规律：()()()mmccXzXVzXVz=+(1)其中，()Xz为功能梯度材料的物理属性如密度、杨氏模量等，mX和cX为两种组分材料的物理属性，()mVz和()cVz为两种材料的体积分数，且满足：()()1mcVzVz+=(2)采用幂函数形式，则体积分数表示为：()202kmzhVkh+=≥(3)式中k为梯度指数。代入式(1)可得()()22kmcczhXzXXXh+=−+(4)2.2.改进的Rayleigh-Ritz法改进的Rayleigh-Ritz法与传统的Rayleigh-Ritz法相比的改进之处在于对边界条件的处理以及位移容许函数的选取方法。2.2.1.函数域的拓展对Rayleigh-Ritz法中的函数域进行扩展，模型如图1。Figure1.Calculationmodelofarbitraryshapedplate图1.任意形状板的计算模型玉云艳等DOI:10.12677/ojav.2019.7100431声学与振动S为薄板的求解域，Sa为试函数的积分域，长度分别为Lx和Ly。当求解域S为任意的不规则形状时，可以假想将结构域扩展为一个略大于真实域的矩形域Sa，矩形域要涵盖待求取结构的真实求解域。一般来说，取求解域距x、y轴最大距离作为矩形域的边界。当采用改进的傅里叶级数作为结构的试函数时，使用矩形域能更便捷地构造出合理的试函数。而在求板的应变能、动能时，只需对真实求解域进行离散并积分，然后求取泛函极值得到振动方程。2.2.2.曲边边界的虚拟弹簧处理方法经典的边界条件，比如固支和简支等往往不能准确模拟结构的实际边界，对于实际边界条件，可采用弹簧模型，通过在结构域的边界上采用平动弹簧和转动弹簧来模拟弹性边界。以矩形板为例施加弹簧边界，如图2所示。Figure2.Modelofrectangularplate图2.矩形板物理模型假设平动弹簧和转动弹簧刚度值分别为kij(N/m)和Kij(Nm/rad)，各种复杂边界条件可通过设置两种弹簧的刚度值模拟。经典边界条件下的刚度值如表1所示。模拟弹性边界条件时则将刚度系数取为相应的值即可。Table1.Classicalboundaryconditionscorrespondingspringvalues.表1.经典边界条件下的边界弹簧取值固支(C)简支(S)自由(F)kl∞∞0K∞00通过弹性边界设置，这将任意边界对结构振动的影响转化为边界弹簧的弹性势能对结构总刚度的影响。3.公式推导功能梯度板的振型函数可以表示为：()()()()()()()()()000000,,e,,e,,eMNitmnmnmnMNitxmnmnmnMNitymnmnmnwxytAfxgyxytBfxgyxytCfxgyωωωββ=========∑∑∑∑∑∑(4)其中Amn，Bmn，Cmn为未知系数，()mfx，()ngy分别为与x、y有关的函数，eitω为时间因子，在下文中忽略。本文方法选用改进的傅里叶级数作为试函数[11]，具体形式为：玉云艳等DOI:10.12677/ojav.2019.7100432声学与振动()()()sin055πcos5mmmfxxmLxmfxxmLxπ=−=≥(5a)()()()πsi