一种基于连续系统仿真验证的剖分有限元方法技术方案

本发明专利技术公开一种基于连续系统仿真验证的剖分有限元方法。该方法针对目前联立法求解动态优化命题的缺陷，即只保证在有限元配置点上满足约束，不保证在非配置点上满足约束，提出了一种基于连续系统仿真验证的剖分有限元方法。该方法采用连续系统仿真验证方法使得有限元非配置点满足约束，采用剖分有限元方法对原有限元序列进行重新配置，能够得到在尽量少增加有限元个数的前提下保证求解精度的有限元配置方案。本发明专利技术针对三维空间下的自由飞行冲突解脱问题采用该方法进行求解，求解结果表明，与原联立法相比，采用这种方法更好的满足了离散化精度和求解精度。

【技术实现步骤摘要】

【技术保护点】
一种基于连续系统仿真验证的剖分有限元方法，其特征在于该方法包括以下步骤：步骤(1).对待测飞行冲突时域tf等分为Nfe段，Nfe≤5，得到Nfe个有限元，则经过等分时域后的有限元序列为每个有限元长度为步骤(2).采用有限元正交配置方法对步骤(1)中的有限元序列进行配置，每个有限元中有3个Radau配置点，每两个配置点之间的点均为非配置点，具体是：将三维空间下的自由飞行冲突解脱问题表述成命题，动态优化命题的一般形式为：minφ(z(tf))dz(t)dt=f(z(t),y(t),u(t),p),z(0)=z0gE(z(t),y(t),u(t),p)=0gI(z(t),y(t),u(y),p)≤0t∈[t0,tf]---(1);]]>其中是φ目标函数，是微分变量，是代数变量，是控制变量，是模型参数，nz、ny、nu、np分别是微分变量、代数变量、控制变量、模型参数的个数，z0是微分变量的初始状态。f表示微分方程，gB表示代数方程，微分代数方程组中微分变量、代数变量和控制变量的边界约束归结到gI中；有限元正交配置方法通过有限元上的正交多项式逼近控制变量和状态变量，定义有限元个数为Nfe，那么t0＜t1＜…＜tNfe＝tf，hi＝ti‑ti‑1；每个有限元hi上通过Lagrange插值多项式对微分变量、控制变量和代数变量进行逼近：z(t)=Σj=0Klj(τ)zi,ju(t)=Σj=1KlJ‾(τ)ui,jy(t)=Σj=1KlJ‾(τ)yi,jt=ti-1+hiτt∈[ti-1,ti],τ∈[0,1]---(2);]]>其中，K为插值的阶次，lj(τ)和分别表示微分变量和控制变量,代数变量的Lagrange插值多项式，可以用如下形式表示：lj(τ)=Πk=0,k≠jKτ-τkτj-τklJ‾(τ)=Πk=1,k≠jKτ-τkτj-τk---(3);]]>Lagrange多项式插值具有如下性质：lj(τk)=1,k=j0,k≠j,lJ‾(τk)=1,k=j0,k≠---(4);]]>即各个变量在配置点上的值正好等于其系数，那么z(ti，j)＝zi，j，y(ti，j)＝yi，j，u(ti，j)＝ui，j    (5)；由于微分变量需要保持状态的连续性，所以在有限元端点上需要通过连接方程来保证微分变量的连续性，而代数变量和控制变量则可以不连续；此外，还必须加上初始和终端条件：zi+1=Σj=0Klj(1)zi,j,z1,0=z(t0),zf=Σj=0Klj(1)zNE,j---(6);]]>将公式(2)和(4)代入(1)并结合连接条件方程(6)，可以得到离散化后NLP问题形式如下：步骤(3).对步骤2得到的离散化后的NLP问题进行求解，得到最优飞行轨迹方案：采用如下的动力学模型进行描述每一架飞行器：dxi(t)dt=vxi(t)dyi(t)dt=vyi(t)dzi(t)dt=vzi(t)---(8);]]>i＝1....N其中，(xi(t)，yi(t)，zi(t))为第i架飞机在t时刻的位置坐标，单位为m；为第i架飞行器在t时刻的速度，单位为m/s，且是操纵飞行器的控制变量；根据ATC标准的飞行安全边界条件可知，任意两架飞行器在同一高度时两者的水平距离Rij不小于R(R＝5nmi)或两者的垂直距离Hij不小于H(H＝1000ft)，可以采用如下的析取表达式进行描述：其中，Rij为第i架飞行器和第j架飞行器之间的水平距离，单位为m；Hij为第i架飞机和第j架飞行器之间的垂直距离，单位为m；假定需在t0到tf这段时间内完成冲突解脱的过程，各架飞行器的初始状态和终止状态为：xi(0)=xi,t0,yi(0)=yi,t0,zi(0)=zi,t0]]>xi(tf)=xi,tf,yi(tf)=yi,tf,zi(tf)=zi,tf---(10);]]>其中，为第i架飞行器的初始位置坐标，单位为m；为第i架飞行器的终止位置坐标，单位为m；t0为开始进行冲突解脱的时刻，单位为s；tf为结束冲突解脱的时刻，单位为s,且tf‑t0＝15min；本方法将飞行冲突解脱过程中的能量消耗作为性能指标，即目标函数J见公式(11)：j=Σi=1Nvi12∫t=t0tf((dxidt)2+(dyidt)2+η2(dzidt)2)2dt---(11...

【技术特征摘要】

【专利技术属性】
技术研发人员：颜丰琳，陈伟锋，邵之江，
申请(专利权)人：浙江大学，
类型：发明
国别省市：浙江;33