经典理论下的功能梯度材料板自振分析.docVIP

经典理论下的功能梯度材料板自振分析 　　【摘要】基于经典板理论，研究了物理中面下功能梯度材料（FGM）板的自由振动响应。应用哈密顿原理推导了经典理论下任意形状任意边界条件的有限元形式的无阻尼自由振动方程，并与均匀板自振方程比较两者之间的相似性，由此得到了两者之间固有频率的相似转换关系。最后给出计算实例来验证了经典理论下FGM板与均匀板之间固有频率的相似转换关系，给出有限元数值模拟结果与其他文献结果的比较。 　　【关键词】物理中面；FGM；自由振动；转换关系 　　前言 　　功能梯度材料在1984年经日本学者平井敏雄提出后，因其优于传统层合式复合材料的力学性能以及良好的可设计性，由早期主要用作航空航天设备的涂层，目前已经扩展到各种材料体系和结构构件中，在许多工程应用领域都有广阔的应用前景[1-5]。He等[6]用有限元法，对上下表面粘贴压电层的功能梯度板的振动问题和稳定问题进行了研究，其选取的组分材料性质是温度相关的，分析了体积分数对振动响应的影响；Yang和Shen[7]通过微分求积法和模态叠加法对施加了横向冲击荷载的FGM板进行了研究，板被置于弹性地基上，研究了梯度指数、地基刚度、板的形状、冲击荷载的持续时间和面内初应力等对动态响应的影响，发现材料性质介于两种基材之间的FGM板，其固有频率和响应也介于两者之间；Abrate[8]通过分析大量文献中的数值结果，将文献中的数据进行曲线拟合后发现FGM板的弯曲位移、固有频率和临界荷载与相应的均匀板之间存在近似的比例关系，并通过分析控制方程中的刚度系数给出消除拉弯耦合项的条件，这一方法与Zhang等[9]中物理中面的定义方式相同；Li等[10]进一步讨论了功能梯度板与对应的均匀板固有频率之间的相似转换关系，用解析方法推导出相似转换系数的解析表达式，并通过对比文献中的数值得到这一频率之间的相似转换关系对于考虑剪切变形的FGM板也近似成立的结论。 　　同文献[6-10]，本文基于物理中面和经典板理论，通过哈密顿原理推导有限元形式的振动方程，推导了单元刚度矩阵和单元质量矩阵，得到了与文献[10]相同的相似转换关系的解析表达式，从而可将功能梯度板自由振动问题转化为求解相同尺寸、荷载及边界条件的均匀板的计算问题。借助高精度的功能梯度板有限元数值解和其他文献的结果，来证明转换关系的准确性。 　　1 问题的基本方程 　　由方程（14）可知，功能梯度材料板的自由振动频率解可以由相同边界条件和相同尺寸的均匀板的自由振动频率解转换得到，它们之间成精确比例关系，其中参数 可以理解为材料性质的非均匀性对功能梯度材料板自振频率的影响。Abrate[8]通过分析大量文献中的数值结果，也得出了这一结论；李世荣等[10]则通过消除功能梯度材料板的自由振动控制微分方程中的面内位移，得到了相似转换关系的解析表达式，并得出对于考虑横向剪切变形的中厚度功能梯度板并不精确满足这个比例关系的结论。 　　3 数值结果比较 　　经典理论下，功能梯度薄板自由振动有限元方程的求解采用科学计算软件MATLAB来计算。考虑长度为 ，宽度为 ，厚度为 的功能梯度矩形薄板，下表面为金属，上表面为陶瓷，组分材料参数由表1给出。 　　已知四边简支矩形均匀板的自由振动频率的解析解为 。考虑 ， 的四边简支矩形氧化铝板，借助FGM薄板自由振动有限元MATLAB程序（令 ），分别求得在不同网格划分下的自由振动频率参数 。他们与解析解的对比列于表2。 　　由上表可以看出，方板网格划分为100×100时，有限元计算所得结果已经与解析解的结果非常接近，但是求解所花费的时间也呈几何级数上升；当划分为50×50时，前5阶频率的最大相对误差为0.13%，最后决定采用50×50的网格划分。 　　表3给出由转换公式所得Al/Al2O3方板前三阶固有频率与文献[11]中的精确解以及有限元（FEM）所得结果的比较。计算中取无量纲频率参数为 ，给定板的厚宽比为 。 　　结论 　　基于物理中面和经典板理论，假设FGM板材料参数随厚度方向连续变化，推导出有限元形式的振动方程，根据物理中面的定义，对比均匀板自振方程，得到相同条件下FGM板与均匀板的固有频率之间的相似转换关系。通过有限元数值计算结果与相似转换结果以及其他文献结果进行对比，证明了转换关系的正确性，便于工程应用。 　　参考文献 　　[1]沈惠申. 功能梯度复合材料板壳结构的弯曲、屈曲和振动[J]. 力学进展， 2004（01）：53-60. 　　[2]仲政， 吴林志， 陈伟球. 功能梯度材料与结构的若干力学问题研究进展[J]. 力学进展， 2010（05）：528-541. 　　[3]李华东， 朱锡， 梅志远， 等. 功能梯度板壳的力学研究进展[J]. 材料导报， 2012（0