本发明属于高超声速风洞试验技术领域,具体涉及一种连续在轨运动的主体和分离体轨迹捕获实验运动分配方法。
背景技术:
目前,全球各地的飞行器和航天器以及飞行武器研究都离不开风洞实验,风洞实验是分析物体气动特性、研究气流规律、研究相邻飞行器气动力相互作用的有效途径。风洞实验通过缩减飞行器尺寸,营造稳定的高超声速流场,获得最真实、最可靠的气动数据,以验证飞行器的气动特性。
当前的风洞实验已经不满足于仅获得飞行器在固定姿态下的气动力特性,还需要进一步研究飞行器在流场中的运动轨迹。轨迹捕获实验是根据当前气动力解算下一时间步长的飞行器模型理论位姿,并控制飞行器模型运动到理论位姿,逐步捕获飞行器模型运动轨迹的风洞实验。
联动轨迹捕获实验由单一飞行器轨迹捕获实验发展而来,目的在于捕获分离体与主体分离时的相对运动轨迹,以研究分离体与主体之间的相互作用效果。相比原来单自由度的轨迹捕获实验,能够直接将分离体与主体之间的相互作用、激波效应等直接模拟,便于验证理论计算与分析。联动捕获轨迹实验时,主体模型和分离体模型中的测力天平分别测量主体模型和分离体模型当前时刻的当前位姿的气动力数据,根据获得的气动力数据计算下一时间步长的主体模型和分离体模型的理论位姿,并控制主体机构和分离体机构运动至理论位姿。由于主体模型和分离体模型在同一流场中,极易产生同向运动,而主体机构和分离体机构的各自由度行程有限,主体机构或分离体机构将在主体模型或分离体模型的某一自由度到达极限时停止运行,显然,这对实验空间的利用是极低的,会大幅降低捕获的轨迹点位数量。
因此,联动轨迹捕获实验中的重要技术问题之一是联动轨迹捕获实验中主体和分离体解耦机构运动分配,当前亟需发展一种连续在轨运动的主体和分离体轨迹捕获实验运动分配方法。
技术实现要素:
本发明所要解决的技术问题是提供一种连续在轨运动的主体和分离体轨迹捕获实验运动分配方法。
本发明的连续在轨运动的主体和分离体轨迹捕获实验运动分配方法,其特点是,所述的运动分配方法使用的解耦机构包括串联或并联的主体机构和分离体机构;
所述的主体机构具有六个自由度,分别由串联的六部分机构来控制,包括主体x向直线移动机构、主体y向直线移动机构、主体z向直线移动机构、主体α俯仰转动机构、主体β偏航转动机构和主体γ滚转机构;主体x向直线移动机构、主体z向直线移动机构、主体y向直线移动机构均由电机与丝杠配合驱动,主体α俯仰转动机构、主体β偏航转动机构通过电机与直线变圆弧机构实现,主体γ滚转机构由滚转电机配合减速器驱动;主体模型固定在主体机构执行末端的主体支杆上,主体模型的质心位于主体机构的解耦点,主体机构状态为主体解耦机构状态;
所述的分离体机构具有六个自由度,分别由串联的六部分机构来控制,包括分离体x向直线移动机构、分离体y向直线移动机构、分离体z向直线移动机构、分离体α俯仰转动机构、分离体β偏航转动机构和分离体γ滚转机构;分离体x向直线移动机构、分离体z向直线移动机构、分离体y向直线移动机构均由电机与丝杠配合驱动,分离体α俯仰转动机构、分离体β偏航转动机构通过电机与直线变圆弧机构实现,分离体γ滚转机构由滚转电机配合减速器驱动;分离体模型固定在分离体机构执行末端的分离体支杆上,分离体模型的质心位于分离体机构的解耦点,分离体机构状态为分离体解耦机构状态;
所述的运动分配方法包括以下步骤:
a.采集主体模型和分离体模型当前时刻的当前位姿;
主体模型当前位姿azo=(xzo,yzo,zzo,αzo,βzo,θzo),分离体模型当前位姿afo=(xfo,yfo,zfo,αfo,βfo,θfo);
b.通过当前时刻的当前位姿的主体模型和分离体模型的气动力数据,计算主体模型和分离体模型下一时间步长的理论位姿;
主体模型理论位姿azn=(xzn,yzn,zzn,αzn,βzn,θzn),分离体模型理论位姿afn=(xfn,yfn,zfn,αfn,βfn,θfn);
c.从当前位姿到达理论位姿,主体机构的三个角度自由度机构的运动角度分别为αzn-αz0,βzn-βz0,θzn-θz0;分离体机构的三个角度自由度机构的运动角度分别为αfn-αf0,βfn-βf0,θfn-θf0;
d.计算当前位姿和理论位姿状态下,主体模型和分离体模型的x向、y向、z向相对距离和相对距离变动量;
x向,主体模型当前位姿和分离体模型当前位姿相对距离xfo-xzo、主体模型理论位姿和分离体模型理论位姿相对距离xfn-xzn,相对距离变动量sx=xfn-xzn-xfo+xzo;
y向,主体模型当前位姿和分离体模型当前位姿相对距离yfo-yzo、主体模型理论位姿和分离体模型理论位姿相对距离yfn-yzn,相对距离变动量sy=yfn-yzn-yfo+yzo;
z向,主体模型当前位姿和分离体模型当前位姿相对距离zfo-zzo、主体模型理论位姿和分离体模型理论位姿相对距离zfn-zzn,相对距离变动量sz=zfn-zzn-zfo+zzo;
e.计算主体模型当前位姿和分离体模型当前位姿在x向、y向、z向的负端余量和正端余量;
将主体机构的x向、y向、z向的运动行程分别表示为[xzmin,xzmax]、[yzmin,yzmax]、[zzmin,zzmax],则主体机构x向、y向、z向的负端余量分别为xzo-xzmin、yzo-yzmin、zzo-zzmin,正端余量分别为xzmax-xzo、yzmax-yzo、zzmax-zzo;
将分离体机构的x向、y向、z向的运动行程分别表示为[xfmin,xfmax]、[yfmin,yfmax]、[zfmin,zfmax],则分离体机构x向、y向、z向的负端余量分别为xfo-xfmin、yfo-yfmin、zfo-zfmin,正端余量分别为xfmax-xfo、yfmax-yfo、zfmax-zfo;
f.按行程余量比重分配主体机构和分离体机构的运动;
主体机构的x向运动量为sx*(xzo-xzmin)/(xzmin-xzmax),数值为正,则运行方向为主体机构的x向正向,否则为主体机构的x向负向;
主体机构的y向运动量为sy*(yzo-yzmin)/(yzmin-yzmax),数值为正,则运行方向为主体机构的y向正向,否则为主体机构的y向负向;
主体机构的z向运动量为sz*(zzo-zzmin)/(zzmin-zzmax),数值为正,则运行方向为主体机构的z向正向,否则为主体机构的z向负向;
分离体机构的x向运动量为sx*(xfo-xfmin)/(xfmin-xfmax),数值为正,则运行方向为分离体机构的x向正向,否则为分离体机构的x向负向;
分离体机构的y向运动量为sy*(yfo-yfmin)/(yfmin-yfmax),数值为正,则运行方向为分离体机构的y向正向,否则为分离体机构的y向负向;
分离体机构的z向运动量为sz*(zfo-zfmin)/(zfmin-zfmax),数值为正,则运行方向为分离体机构的z向正向,否则为分离体机构的z向负向;
g.控制主体机构和分离体机构分别完成步骤c的角度运动和步骤f的x向、y向、z向运动,到达下一时间步长的主体模型理论位姿和分离体模型理论位姿。
所述的下一时间步长为常值。
所述的步骤a的具体步骤如下:
a1.获取当前时刻主体机构和分离体机构各自由度电机的编码器位置,编码器的值能够准确反映出电机转动的角度,再根据丝杠导程计算丝杠传动的距离:
a2.主体γ滚转机构和分离体γ滚转机构是电机配合减速器直接运动,主体机构和分离体机构的滚转角θγ(t)与滚转电机转角γ(t)关系为:
a3.主体α俯仰转动机构和分离体α俯仰转动机构采用直线变圆弧结构;
a点为丝杠上滑块初始位置,坐标为a(x0,y0),b点为连杆与俯仰机构的转动轴处,转动轴与俯仰中心的连线与x向夹角为
分别计算获得αzo和αfo;
a4.主体β偏航转动机构和分离体β偏航转动机构采用直线变圆弧结构;
c点为丝杠上滑块初始位置,坐标为c(x0,z0),d点为连杆与偏航机构的转动轴处,转动轴与偏航中心的连线与x向夹角为ω0;偏航转轴做弧形运动的圆弧半径od为r′,连杆cd长度为l′;设偏航电机驱动偏航直线滑块c运动位移cc′为s′,则连杆带动偏航部分绕o做旋转运动;偏航转动角度β计算公式如下:
分别计算获得βzo和βfo。
本发明的连续在轨运动的主体和分离体轨迹捕获实验运动分配方法的主体模型和分离体模型同时置于风洞流场中。该方法将主体模型和分离体模型分别安装在主体机构和分离体机构的解耦点,主体机构和分离体机构任一自由度的运动都不会引起其余自由度的变化,主体机构状态的三个角度自由度从(αzo,βzo,θzo)到(αzn,βzn,θzn),分离体机构状态的三个角度自由度从(αfo,βfo,θfo)到(αfn,βfn,θfn),不存在补偿需求且不必考虑相对角度关系,主体机构的三个角度自由度机构直接将主体模型绕质心转动至(αzn,βzn,θzn)即可,分离体机构的三个角度自由度机构直接将分离体模型绕质心转动至(αfn,βfn,θfn)即可。
本发明的连续在轨运动的主体和分离体轨迹捕获实验运动分配方法通过对联动轨迹捕获实验中主体和分离体解耦机构运动分配,避免了主体机构和分离体机构同方向运动带来的短时间内到达行程极限而停止实验的情况。通过按行程余量比重分配主体机构和分离体机构的运动,提高了实验空间的利用率,有效增加了捕获的轨迹点位数量,获得更完整的主体模型和分离体模型运动轨迹。
本发明的连续在轨运动的主体和分离体轨迹捕获实验运动分配方法以保证主体模型和分离体模型的相对位姿正确为前提,通过按行程余量比重分配主体机构和分离体机构的运动,能够在主体机构和分离体机构中合理分配x向、y向、z向的直线运动距离,从而在时间步长固定的情况下,减小主体机构和分离体机构x向、y向、z向电机的最高运动速度,减小对主体机构和分离体机构x向、y向、z向电机的较大冲击,降低主体机构和分离体机构x向、y向、z向电机负载,提高了主体机构和分离体机构的使用寿命。
本发明的连续在轨运动的主体和分离体轨迹捕获实验运动分配方法只需按行程余量比重分配主体机构和分离体机构的运动,方法简洁,不用进行繁琐的计算,高效可靠。
综上所述,本发明的连续在轨运动的主体和分离体轨迹捕获实验运动分配方法能够避免主体机构和分离体机构同向运动,提高实验空间的利用率,获得更长的实验时间和更多的捕获轨迹点位数量,提高主体机构和分离体机构的使用寿命,具有计算简单快速,运行高效可靠的优点。
附图说明
图1为本发明的连续在轨运动的主体和分离体轨迹捕获实验运动分配方法使用的主体机构和分离体机构示意图;
图2为本发明的连续在轨运动的主体和分离体轨迹捕获实验运动分配方法中的α俯仰转动机构直线转圆弧简图;
图3为本发明的连续在轨运动的主体和分离体轨迹捕获实验运动分配方法中的β偏航转动机构直线转圆弧简图;
图4为本发明的连续在轨运动的主体和分离体轨迹捕获实验运动分配方法的按行程余量比重分配主体机构和分离体机构的运动流程图。
具体实施方式
下面结合附图和实施例详细说明本发明。
以下描述用于揭示本发明以使本领域技术人员能够实现本发明。以下描述中的优选实施例只作为举例,本领域技术人员可以想到其他显而易见的变型。在以下描述中界定的本发明的基本原理可以应用于其他实施方案、变形方案、改进方案、等同方案以及没有背离本发明的精神和范围的其他技术方案。
实施例1
如图1所示,本实施例的分离体机构和主体机构以上下结构的形式安装,分离体模型和主体模型在实验时都被置于高超声速风洞流场中。
所述的主体机构具有六个自由度,分别由串联的六部分机构来控制,包括主体x向直线移动机构、主体y向直线移动机构、主体z向直线移动机构、主体α俯仰转动机构、主体β偏航转动机构和主体γ滚转机构;主体x向直线移动机构、主体z向直线移动机构、主体y向直线移动机构均由电机与丝杠配合驱动,主体α俯仰转动机构、主体β偏航转动机构通过电机与直线变圆弧机构实现,主体γ滚转机构由滚转电机配合减速器驱动;主体模型固定在主体机构执行末端的主体支杆上,主体模型的质心位于主体机构的解耦点,主体机构状态为主体解耦机构状态;
所述的分离体机构具有六个自由度,分别由串联的六部分机构来控制,包括分离体x向直线移动机构、分离体y向直线移动机构、分离体z向直线移动机构、分离体α俯仰转动机构、分离体β偏航转动机构和分离体γ滚转机构;分离体x向直线移动机构、分离体z向直线移动机构、分离体y向直线移动机构均由电机与丝杠配合驱动,分离体α俯仰转动机构、分离体β偏航转动机构通过电机与直线变圆弧机构实现,分离体γ滚转机构由滚转电机配合减速器驱动;分离体模型固定在分离体机构执行末端的分离体支杆上,分离体模型的质心位于分离体机构的解耦点,分离体机构状态为分离体解耦机构状态;
所述的运动分配方法包括以下步骤:
a.采集主体模型和分离体模型当前时刻的当前位姿;
主体模型当前位姿azo=(xzo,yzo,zzo,αzo,βzo,λzo),分离体模型当前位姿afo=(xfo,yfo,zfo,αfo,βfo,θfo);
具体步骤如下:
a1.获取当前时刻主体机构和分离体机构各自由度电机的编码器位置,编码器的值能够准确反映出电机转动的角度,再根据丝杠导程计算丝杠传动的距离:
a2.主体γ滚转机构和分离体γ滚转机构是电机配合减速器直接运动,主体机构和分离体机构的滚转角θγ(t)与滚转电机转角γ(t)关系为:
a3.如图2所示,主体α俯仰转动机构和分离体α俯仰转动机构采用直线变圆弧结构;
a点为丝杠上滑块初始位置,坐标为a(x0,y0),b点为连杆与俯仰机构的转动轴处,转动轴与俯仰中心的连线与x向夹角为
分别计算获得αzo和αfo;
a4.如图3所示,主体β偏航转动机构和分离体β偏航转动机构采用直线变圆弧结构;
c点为丝杠上滑块初始位置,坐标为c(x0,z0),d点为连杆与偏航机构的转动轴处,转动轴与偏航中心的连线与x向夹角为ω0;偏航转轴做弧形运动的圆弧半径od为r′,连杆cd长度为l′;设偏航电机驱动偏航直线滑块c运动位移cc′为s′,则连杆带动偏航部分绕o做旋转运动;偏航转动角度β计算公式如下:
分别计算获得βzo和βfo。
b.通过当前时刻的当前位姿的主体模型和分离体模型的气动力数据,计算主体模型和分离体模型下一时间步长的理论位姿;
主体模型理论位姿azn=(xzn,yzn,zzn,αfn,βzn,θzn),分离体模型理论位姿afn=(xfn,yfn,zfn,αfn,βfn,θfn);
c.从当前位姿到达理论位姿,主体机构的三个角度自由度机构的运动角度分别为αzn-αz0,βzn-βz0,θzn-θz0;分离体机构的三个角度自由度机构的运动角度分别为αfn-αf0,βfn-βf0,θfn-θf0;
d.计算当前位姿和理论位姿状态下,主体模型和分离体模型的x向、y向、z向相对距离和相对距离变动量;
x向,主体模型当前位姿和分离体模型当前位姿相对距离xfo-xzo、主体模型理论位姿和分离体模型理论位姿相对距离xfn-xzn,相对距离变动量sx=xfn-xzn-xfo+xzo;
y向,主体模型当前位姿和分离体模型当前位姿相对距离yfo-yzo、主体模型理论位姿和分离体模型理论位姿相对距离yfn-yzn,相对距离变动量sy=yfn-yzn-yfo+yzo;
z向,主体模型当前位姿和分离体模型当前位姿相对距离zfo-zzo、主体模型理论位姿和分离体模型理论位姿相对距离zfn-zzn,相对距离变动量sz=zfn-zzn-zfo+zzo;
e.计算主体模型当前位姿和分离体模型当前位姿在x向、y向、z向的负端余量和正端余量;
将主体机构的x向、y向、z向的运动行程分别表示为[xzmin,xzmax]、[yzmin,yzmax]、[zzmin,zzmax],则主体机构x向、y向、z向的负端余量分别为xzo-xzmin、yzo-yzmin、zzo-zzmin,正端余量分别为xzmax-xzo、yzmax-yzo、zzmax-zzo;
将分离体机构的x向、y向、z向的运动行程分别表示为[xfmin,xfmax]、[yfmin,yfmax]、[zfmin,zfmax],则分离体机构x向、y向、z向的负端余量分别为xfo-xfmin、yfo-yfmin、zfo-zfmin,正端余量分别为xfmax-xfo、yfmax-yfo、zfmax-zfo;
f.如图4所示,按行程余量比重分配主体机构和分离体机构的运动;
主体机构的x向运动量为sx*(xzo-xzmin)/(xzmin-xzmax),数值为正,则运行方向为主体机构的x向正向,否则为主体机构的x向负向;
主体机构的y向运动量为sy*(yzo-yzmin)/(yzmin-yzmax),数值为正,则运行方向为主体机构的y向正向,否则为主体机构的y向负向;
主体机构的z向运动量为sz*(zzo-zzmin)/(zzmin-zzmax),数值为正,则运行方向为主体机构的z向正向,否则为主体机构的z向负向;
分离体机构的x向运动量为sx*(xfo-xfmin)/(xfmin-xfmax),数值为正,则运行方向为分离体机构的x向正向,否则为分离体机构的x向负向;
分离体机构的y向运动量为sy*(yfo-yfmin)/(yfmin-yfmax),数值为正,则运行方向为分离体机构的y向正向,否则为分离体机构的y向负向;
分离体机构的z向运动量为sz*(zfo-zfmin)/(zfmin-zfmax),数值为正,则运行方向为分离体机构的z向正向,否则为分离体机构的z向负向;
g.控制主体机构和分离体机构分别完成步骤c的角度运动和步骤f的x向、y向、z向运动,到达下一时间步长的主体模型理论位姿和分离体模型理论位姿。
所述的下一时间步长为常值。
以上对本发明的一个实施例进行了详细说明,但该内容仅为本发明的较佳实施例,不能被认为用于限定本发明的实施范围。凡依本发明申请范围所作的均等变化与改进等,均应仍归属于本发明的专利涵盖范围之内。