[发明专利]基于无网格法考虑热效应及剪切效应FGM梁的动力学仿真模型、其建立方法及仿真方法

说明书

本案涉及一种动力学仿真模型、其建立方法及仿真方法，包括：建立考虑热效应及剪切效应下的附加集中质量的FGM梁系统；采用混合坐标法在浮动坐标系中描述附加集中质量的FGM梁系统中的FGM梁上任一点作大范围旋转运动变形的位移场；对FGM梁的大范围旋转变形进行离散；运用第二类Lagrange方程建立附加集中质量的FGM梁系统的一次近似刚柔耦合动力学方程，即得基于无网格法考虑热效应及剪切效应FGM梁的动力学仿真模型；采用Newmark法求解热效应下附加集中质量的FGM梁的动力学仿真模型，得到大范围旋转运动的附加集中质量的FGM梁的变形示意图。本案的仿真模型和仿真方法具有效率高、精度高的优点。

技术领域

本发明涉及系统动力学建模领域，具体涉及一种基于无网格法考虑热效应及剪切效应FGM梁的动力学仿真模型、其建立方法及仿真方法。

背景技术

工程结构中的直升机旋翼叶片、机械臂和涡轮机叶片等正朝着更高的运行速度和精度的方向发展。这些系统由刚性基座和连接在其上的柔性附件组成，是典型的刚柔耦合系统。在实际应用中，系统往往处于高温、高转速等工况，传统的均质材料由于其自身力学性能的局限性已很难满足实际服役环境的要求。因此，在不影响结构强度的前提下，设计出满足服役环境要求的新型复合材料并研究其动力学特性将尤为重要。

在众多新型复合材料中，功能梯度材料(Functionally graded material,FGM)以其耐热、高强度、改善应力集中等优越性而被研究人员广泛关注，并将其应用于航空航天、核工业及生物等领域。

现有技术中对中心刚体-功能梯度材料梁系统进行动力学建模，建模时仅考虑横向弯曲变形量、轴向变形量及横向变形引起的轴向缩短量，忽略了温度对系统的影响。此外，多体系统动力学的相关理论大都采用Euler-Bernoulli梁进行研究，该理论只适用于计算细长梁，对于长细比较小的梁系统，由于Euler-Bernoulli梁理论忽略了剪切效应，会导致一定的系统误差。当采用假设模态法对该系统进行离散时，对于梁发生小变形时，假设模态法可以很好的描述出梁的变形，但当柔性机械臂变形较大时，该方法是存在理论缺陷的。

发明内容

本发明的目的是基于计及剪切效应的Timoshenko梁理论，提供一种新的离散模型，以期解决热效应、剪切效应及附加集中质量对旋转FGM梁动力学特性的影响。

为实现上述目的，本发明所采用的技术方案为：

一种基于无网格法考虑热效应及剪切效应FGM梁的动力学仿真模型的建立方法，其至少包括：

设定FGM梁的几何参数、材料参数，建立考虑热效应及剪切效应下的附加集中质量的FGM梁系统；

采用混合坐标法在浮动坐标系中描述附加集中质量的FGM梁系统中的FGM梁上任一点作大范围旋转运动变形的位移场；

采用无网格点插值法和无网格径向基点插值法对FGM梁的大范围旋转变形进行离散；

运用第二类Lagrange方程建立附加集中质量的FGM梁系统的一次近似刚柔耦合动力学方程，即得基于无网格法考虑热效应及剪切效应FGM梁的动力学仿真模型。

优选的是，所述的建立方法，其中，所述FGM梁的几何参数为：梁的长度为l，宽度为b，集中质量为m_t，集中质量距离固定端的距离为l_t；材料参数为：FGM梁的密度ρ、弹性模量E、热传导系数为K，线膨胀系数为α，表达式分别如下：

其中，角标c表示陶瓷材料，角标m表示金属材料，h为FGM梁的厚度，N为功能梯度指数，y为FGM梁在厚度方向的位置。

优选的是，所述的建立方法，其中：FGM梁上任意一点的变形的矢径r及集中质量m_t的矢径r_mt在惯性坐标系O-XYZ中的表达式为：