本发明属于多体系统动力学领域，具体地说，是一种中心刚体-fgm楔形梁系统末端动力学响应计算方法。
背景技术：
针对做大范围运动的中心刚体-柔性梁系统的动力学响应问题，选择合适的坐标系建立精确的动力学模型，并在完整模型的基础上进行简化，在满足计算精度的情况下，得到计算效率较高的动力学模型，成为求解这类动力学问题的关键所在。
librescu首次建立了中心刚体-fgm梁模型，并在此基础上对这一模型进行了振动分析。2005年，librescu在以往建模方法的基础上，对圆柱薄壁梁在大范围运动下的动力学特性做出了研究。2012年，张伟将高阶剪切理论运用在建模过程中，考虑到离心力的作用，建立了旋转运动下fgm板的动力学方程。黎亮首次提出了倾角坐标，基于细长梁假设，对fgm梁系统在大范围运动下的动力学问题进行了研究。目前已有的工作主要是运用传统的直角坐标系计算柔性梁的动力学问题，计算效率较为低下，选用合理的建模方法，建立满足计算精度，计算效率较高的模型，成为这类问题研究中的重点。
技术实现要素：
本发明的目的是，针对大范围旋转运动下中心刚体-fgm楔形梁系统末端动力学响应问题，提供一种数值仿真的计算方法，将fgm楔形梁几何参数、功能梯度参数、材料组成分别进行设置，得到fgm楔形梁的末端横向变形与轴向变形。
实现本发明目的地技术解决方案为：一种中心刚体-fgm楔形梁系统末端动力学响应计算方法，包括以下步骤：
(1)设定中心刚体-fgm楔形梁系统相关参数：中心刚体转动惯量、楔形梁几何尺寸、fgm梁组成材料组成、功能梯度指数，并给出大范围运动角速度规律；
(2)选用弧长坐标中心刚体-fgm楔形梁系统进行建模，运用几何关系描述中心刚体-fgm楔形梁系统的变形场，得出柔性梁末端位移表达式；
(3)取中心刚体-fgm楔形梁系统的一段微元进行分析，写出柔性梁系统在大范围转动下的动能和势能表达式；
(4)运用假设模态法对每段微元的横向弯曲角、纵向拉伸量及剪切角进行离散，并将动能与势能带入第二类lagrange方程，并将方程中二次以上项舍去，得到中心刚体-柔性梁系统的刚柔耦合动力学方程；
(5)针对中心刚体-fgm楔形梁系统，运用梁高比rh，梁宽比rb描述楔形梁几何形状；运用梁悬臂端与自由端材料参数及功能梯度参数描述fgm梁材料组成；
(6)根据步骤(4)中动力学方程和步骤(5)给定的参数，得出fgm楔形梁末端横向变形及轴向变形随时间变化规律数据。
步骤(1)中，大范围运动角速度规律为：
式中，ω为转动角速度，ω0为初始转动角速度，t为大范围转动计算时长；
步骤(2)中，柔性梁末端位移表达式为：
式中，u(t)为柔性梁末端轴向变形，v(t)为柔性梁末端横向变形，α(s,t)为弧长坐标s处横截面的弯曲角度，ε(s,t)为弧长坐标s处轴向拉伸量，l为柔性梁长度。
步骤(3)中，柔性梁系统的动能表达式为：
式中，joh为中心刚体转动惯量，θ0为中心刚体角位移，ρ(s)为柔性梁沿轴向密度函数，a(s)为柔性梁沿轴向横截面积函数，x0、y0为梁轴线上一点处坐标分量，γ(s,t)为弧长坐标s处横截面的剪切角。
柔性梁势能表达式为：
式中，e(s)为柔性梁沿轴向弹性模量函数，g(s)为柔性梁沿轴向剪切模量函数，k为剪切修正系数。
步骤(4)中，采用假设模态法对柔性梁的变形进行描述，将倾角、纵向拉伸量、剪切角进行离散处理：
其中，φi(s)为一端固支一端自由杆试函数行向量，a(t)、b(t)、c(t)为与时间相关项列向量。将上式带入第二类lagrange方程，并舍去部分高阶项，得到非惯性系下中心刚体-柔性梁系统的动力学方程：
式中各项分别为：
步骤(5)中梁高比rh，梁宽比rb的取值范围分别为0≤rh≤1，0≤rb≤1，材料参数分别需设定梁固定端及自由端材料密度与弹性模量。
步骤(6)中中心刚体-fgm楔形梁系统需设置的参数分别为：柔性梁长度、悬臂端截面积及惯性矩、梁高比及梁宽比、中心刚体转动惯量、悬臂端及固定端材料密度与弹性模量、功能梯度指数。
本发明与现有技术相比，其显著优点为：
(1)在建模过程中，采用倾角坐标，可通过描述梁中轴线上点的位置，对柔性梁的变形问题进行描述，在建模过程中较传统的坐标系更加方便简洁；
(2)在公式推导中，基于timoshenko梁假设，考虑到了横向弯曲、轴向拉伸及剪切效应；计算精度较高；
(3)在计算末端响应时，在满足精度的前提下，省略了部分高阶项，计算效率较高；
(4)fgm楔形梁参数设置时，可对中心刚体转动惯量、几何尺寸、材料组成、功能梯度参数等进行设置，可运用于多种形式的fgm楔形梁。
附图说明
图1是柔性梁变形示意图。
图2是fgm楔形梁几何尺寸示意图。
图3是文件“0.txt”数据。
图4是c++程序运行过程。
图5是计算完成输出文件“v.txt”及“u.txt”。
图6是本发明方法的实现流程图。
具体实施方式
下面结合实施例及附图对本发明进行进一步介绍。
实例：一种中心刚体-fgm楔形梁系统末端动力学响应计算方法，如图6所示，该方法包括以下步骤：
(1)设定旋转速度规律：
(2)设定中心刚体-fgm梁系统的相关参数，悬臂端材料为铝，自由端材料为陶瓷。具体数值在表1中给出。
表1中心刚体-fgm梁系统的几何参数
表2中心刚体-fgm梁系统的材料参数
(3)在文件“0.txt”中输入数据，如图3所示；
(4)编写c++程序求解动力学方程(4)，运行程序并读入文件“0.txt”中数据，计算过程如图4所示；
(5)程序运行完成后，输出fgm楔形梁系统在大范围运动下横向变形与轴向变形随时间变化规律文件“v.txt”及“u.txt”(图5)，文件中数据可用于生成曲线来研究fgm楔形梁末端变形随时间变化规律。