保守系统动力学方程的符号推导方法及系统、计算机程序技术方案

本发明专利技术属于动力学分析领域，公开了一种保守系统动力学方程的符号推导方法及系统、计算机程序，巧妙将LpY(拉格朗日函数对广义速度的导数)对时间的导数分成了两个部分之和(对广义位移求导乘上广义速度之和；对广义速度求导乘上广义加速度之和)，最终得到了保守系统动力学符号方程。本发明专利技术实现了保守系统动力学方程的自动符号推导，克服了手工推导繁琐、易错、效率低下的缺陷；本发明专利技术符号推导动力学方程的过程比较简单，且容易用计算机程序实现；本发明专利技术的符号推导的动力学方程，代入其初始条件，即可进行数值求解和仿真。该方法能够应用于保守系统的动力学方程的推导，为系统的动力学仿真、分析、优化和控制等奠定了基础。

【技术实现步骤摘要】
保守系统动力学方程的符号推导方法及系统、计算机程序
本专利技术属于动力学分析领域，尤其涉及一种保守系统动力学方程的符号推导方法及系统、计算机程序。
技术介绍
目前，业内常用的现有技术是这样的：动力学方程的准确获得是进行机械系统动力学分析和设计的前提。目前，机械系统的动力学方程的推导主要有三种方法：一是采用矢量动力学的方法，这种方法需要将机械系统隔离开来分析，具有分析过程复杂、自由度数目多以及求解规律大的缺点；二是采用拉格朗日第一类方程来完成，这种方法具有程式化的优点，但方法具有自由度数目多、求解规律大的缺点；三是采用拉格朗日第二类方程手工推导来完成，得到的动力学方程具有方程数目最少，求解规律最少的优点。然而，手工推导的过程十分复杂且容易出错，即便是自由度数目较少的机械系统，利用拉格朗日第二类方程来手工推导其动力学方程也是相当困难的。尽管拉格朗日第二类方程手工推导过程繁琐，但多数人仍然愿意采用这种方法，因为同一系统的动力学模型的动力学方程是相同的，具有“一次推导，终身受用”的特点。随着计算机技术的发展，符号推导演算功能日益强大，成为了复杂机械系统动力学推导强有力的工具。利用计算机技术进行符号推导能够直接导出复杂机械系统动力学微分方程，可避免大量的累积误差，提高计算效率和计算精度，减轻人的劳动并实现快速分析。保守系统是指系统在运动以及变化的过程中，机械能始终不向外流失，动能、势能之和为一恒定值。保守系统在机械系统中是一种典型的系统，利用拉格朗日第二类方程推导其动力学方程具有极强的程式化特点，这种系统适用于计算机符号推导来得到系统的动力学方程。因此，迫切需要利用计算机技术研究一种动力学方程的符号推导方法。综上所述，现有技术存在的问题是：(1)现有技术中利用矢量动力学方法和拉格朗日第一类方程方法具有程式化特征，但得到的方程规模通过较大，不利于方程的求解。(2)现有技术多通过手工推导系统的动力学方程，而手工推导繁琐易错，不利于大规模复杂问题的动力学方程的推导。(3)现有技术没有公开利用拉格朗日第二类方程进行系统的动力学方程的符号推导方法，同时也没有公开利用拉格朗日第二类方程进行系统动力学方程符号推导的过程。解决上述技术问题的难度和意义：本专利技术实现了动力学方程的自动推导，提高计算效率、减轻人的劳动，使人们从繁琐又容易出错的方程推导中解脱出来，将主要的精力投入到更具有创造性的工作中去；目前文献没有公开利用拉格朗日第二类方程进行系统符号推导的方法，因此没有改进的保守系统动力学方程的符号推导方法及系统。
技术实现思路
针对现有技术存在的问题，本专利技术提供了一种保守系统动力学方程的符号推导方法及系统、计算机程序。本专利技术是这样实现的，一种保守系统动力学方程的符号推导方法，为：在求拉格朗日函数对广义速度的导数LpY对时间的导数时，分成对广义位移求导乘上广义速度之和、对广义速度求导乘上广义加速度之和两个部分；最终得到保守系统动力学符号方程。进一步，所述保守系统动力学方程的符号推导方法，包括：定义系统的常量参数；定义系统的符号变量参数；根据系统的自由度数目，定义系统的运动学符号变量参数，分别为广义位移列向量X、广义速度列向量Y和广义加速度列向量Z；利用定义的常量和变量参数写出系统的动能T和势能U，并构造拉格朗日函数L＝T-U；拉格朗日函数L对广义位移求偏导，记为LpX：LpX(i)＝diff(L,X(i))(i＝1,2...,n)(1)拉格朗日函数L对广义速度的偏导，记为LpY：LpY(i)＝diff(L,Y(i))(i＝1,2...,n)(2)式(1)和(2)中：diff(L,X(i))表示L对第i个广义位移的导数，diff(L,Y(i))表示L对第i个广义速度的导数，X(i)表示第i个广义位移，Y(i)表示第i个广义速度，n表示广义位移的数目；对LpY再取时间的导数，记为dLpY；考虑到LpY中既含有广义位移X又含有广义速度Y，故将dLpY分解成两部分：对LpY广义位移求导乘上广义速度之和、对LpY广义速度求导乘上广义加速度之和；第一部分记为：第二部分记为：故dLpY为式(3)和(4)之和，即：dLpY(i)＝dLpYX(i)+dLpYY(i)(i＝1,2...,n)(5)式(3)、(4)和(5)中：diff(LpY(i),X(j))和diff(LpY(i),Y(j))分别表示dLpY的第i个表达式对X的第j个广义位移的导数和对Y的第j个广义速度的导数，X(j)、Y(j)和Z(j)分别表示第j个广义位移、广义速度和广义加速度；利用拉格朗日第二类方程，建立系统的动力学方程为：EQ(i)＝dLpY(i)-LpX(i)(i＝1,2...,n)(6)式(6)中：EQ(i)表示第i个动力学方程。本专利技术的另一目的在于提供一种实现所述保守系统动力学方程的符号推导方法的计算机程序。本专利技术的另一目的在于提供一种一种实现所述保守系统动力学方程的符号推导方法的信息数据处理终端。本专利技术的另一目的在于提供一种计算机可读存储介质，包括指令，当其在计算机上运行时，使得计算机执行所述的保守系统动力学方程的符号推导方法。本专利技术的另一目的在于提供一种四自由度系统动力学模型包括：1)定义四自由度系统的常量参数；常量参数有弹簧的刚度k1、k2，刚体的质量m1、m2、m3和m4，刚体的转动惯量J1、J2、J3和J4,刚体B3和B4的长度分别为L1、L2，重力加速度g；2)定义四自由度系统的符号变量参数；所述弹簧有4个自由度，定义为广义位移列向量X＝[x1,x2,x3,x4]T、广义速度列向量Y＝[y1,y2,y3,y4]T和广义加速度列向量Z＝[z1,z2,z3,z4]T；其中：x1和x2分别表示为刚体B1和刚体B2沿x方向的位移，x3和x4分别表示B3绕C1点的转角和B4绕C2点的转角，y1和y2分别表示为刚体B1和B2沿x方向的速度，y3和y4分别表示B3绕C1点的角速度和B4绕着C2点的角速度，z1和z2分别表示为刚体B1和B2沿x方向的加速度，z3和z4分别表示B3绕C1点的角加速度和B4绕着C2点的角加速度；3)利用定义的四自由度系统的常量和符号变量参数写出系统的动能T和势能U：在此基础上构造拉格朗日函数L＝T-U；拉格朗日函数L对广义位移求偏导，记为LpX：LpX(i)＝diff(L,X(i))(i＝1,2...,n)(1)拉格朗日函数L对广义速度的偏导，记为LpY：LpY(i)＝diff(L,Y(i))(i＝1,2...,n)(2)式(1)和(2)中：diff(L,X(i))表示L对第i个广义位移的导数，diff(L,Y(i))表示L对第i个广义速度的导数，X(i)表示第i个广义位移，Y(i)表示第i个广义速度，n表示广义位移的数目，其中n＝4；对LpY再取时间的导数，记为dLpY；考虑到LpY中既含有广义位移X又含有广义速度Y，故将dLpY分解成两部分：对LpY广义位移求导乘上广义速度之和、对LpY广义速度求导乘上广义加速度之和；第一部分记为：第二部分记为：故dLpY为式(3)和(4)之和，即：dLpY(i)＝dLpYX(i)+dLpYY(i)(i＝1,2...,n)(5)式(3)、(4)和(5)中：diff(LpY(i),X(j))和diff(本文档来自技高网...

【技术保护点】
1.一种保守系统动力学方程的符号推导方法，其特征在于，所述保守系统动力学方程的符号推导方法为：在求拉格朗日函数对广义速度的导数LpY对时间的导数时，分成对广义位移求导乘上广义速度之和、对广义速度求导乘上广义加速度之和两个部分；最终得到保守系统动力学符号方程。

【技术特征摘要】
1.一种保守系统动力学方程的符号推导方法，其特征在于，所述保守系统动力学方程的符号推导方法为：在求拉格朗日函数对广义速度的导数LpY对时间的导数时，分成对广义位移求导乘上广义速度之和、对广义速度求导乘上广义加速度之和两个部分；最终得到保守系统动力学符号方程。2.如权利要求1所述的保守系统动力学方程的符号推导方法，其特征在于，所述保守系统动力学方程的符号推导方法，包括：定义系统的常量参数；定义系统的符号变量参数；根据系统的自由度数目，定义系统的运动学符号变量参数，分别为广义位移列向量X、广义速度列向量Y和广义加速度列向量Z；利用定义的常量和变量参数写出系统的动能T和势能U，并构造拉格朗日函数L＝T-U；拉格朗日函数L对广义位移求偏导，记为LpX：LpX(i)＝diff(L,X(i))(i＝1,2...,n)(1)拉格朗日函数L对广义速度的偏导，记为LpY：LpY(i)＝diff(L,Y(i))(i＝1,2...,n)(2)式(1)和(2)中：diff(L,X(i))表示L对第i个广义位移的导数，diff(L,Y(i))表示L对第i个广义速度的导数，X(i)表示第i个广义位移，Y(i)表示第i个广义速度，n表示广义位移的数目；对LpY再取时间的导数，记为dLpY；考虑到LpY中既含有广义位移X又含有广义速度Y，故将dLpY分解成两部分：对LpY广义位移求导乘上广义速度之和、对LpY广义速度求导乘上广义加速度之和；第一部分记为：第二部分记为：故dLpY为式(3)和(4)之和，即：dLpY(i)＝dLpYX(i)+dLpYY(i)(i＝1,2...,n)(5)式(3)、(4)和(5)中：diff(LpY(i),X(j))和diff(LpY(i),Y(j))分别表示dLpY的第i个表达式对X的第j个广义位移的导数和对Y的第j个广义速度的导数，X(j)、Y(j)和Z(j)分别表示第j个广义位移、广义速度和广义加速度；利用拉格朗日第二类方程，建立系统的动力学方程为：EQ(i)＝dLpY(i)-LpX(i)(i＝1,2...,n)(6)式(6)中：EQ(i)表示第i个动力学方程。3.一种实现权利要求1～2任意一项所述保守系统动力学方程的符号推导方法的计算机程序。4.一种实现权利要求1～2任意一项所述保守系统动力学方程的符号推导方法的信息数据处理终端。5.一种计算机可读存储介质，包括指令，当其在计算机上运行时，使得计算机执行如权利要求1～2任意一项所述的保守系统动力学方程的符号推导方法。6.一种如权利要求1所述保守系统动力学方程的符号推导方法的四自由度系统动力学模型，其特征在于，所述四自由度系统动力学模型包括：1)定义四自由度系统的常量参数；常量参数有弹簧的刚度k1、k2，刚体的质量m1、m2、m3和m4，刚体的转动惯量J1、J2、J3和J4,刚体B3和B4的长度分别为L1、L2，重力加速度g；2)定义四自由度系统的符号变量参数；所述弹簧有4个自由度，定义为广义位移列向量X＝[x1,x2,x3,x4]T、广义速度列向量Y＝[y1,y2,y3,y4]T和广义加速度列向量Z＝[z1,z2,z3,z4]T；其中：x1和x2分别表示为刚体B1和刚体B2沿x方向的位移，x3和x4分别表示B3绕C1点的转角和B4绕C2点的转角，y1和y2分别表示为刚体B1和B2沿x方向的速度，y3和y4分别表示B3绕C1点的角速度和B4绕着C2点的角速度，z1和z2分别表示为刚体B1和B2沿x方向的加速度，z3和z4分别表示B3绕C1点的角加速度和B4绕着C2点的角加速度；3)利用定义的四自由度系统的常量和符号变量参数写出系统的动能T和势能U：在此基础上构造拉格朗日函数L＝T-U；拉格朗日函数L对广义位移求偏导，记为LpX：LpX(i)＝diff(L,X(i))(i＝1,2...,n)(1)拉格朗日函数L对广义速度的偏导，记为LpY...

【专利技术属性】
技术研发人员：伍建伟，刘夫云，赵亮亮，唐荣江，何水龙，杨孟杰，
申请(专利权)人：桂林电子科技大学，
类型：发明
国别省市：广西,45