CN113761653B - 子母弹分离安全性评价方法及系统及子母弹分离仿真平台 - Google Patents

子母弹分离安全性评价方法及系统及子母弹分离仿真平台 Download PDF

Info

Publication number
CN113761653B
CN113761653B CN202110950777.4A CN202110950777A CN113761653B CN 113761653 B CN113761653 B CN 113761653B CN 202110950777 A CN202110950777 A CN 202110950777A CN 113761653 B CN113761653 B CN 113761653B
Authority
CN
China
Prior art keywords
separation
shrapnel
interpolation
aerodynamic
derivative
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110950777.4A
Other languages
English (en)
Other versions
CN113761653A (zh
Inventor
王波兰
刘广
许斌
胡珊
葛志闪
俞刘建
刘陆广
张保刚
李微微
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Institute of Electromechanical Engineering
Original Assignee
Shanghai Institute of Electromechanical Engineering
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Institute of Electromechanical Engineering filed Critical Shanghai Institute of Electromechanical Engineering
Priority to CN202110950777.4A priority Critical patent/CN113761653B/zh
Publication of CN113761653A publication Critical patent/CN113761653A/zh
Application granted granted Critical
Publication of CN113761653B publication Critical patent/CN113761653B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/13—Differential equations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Operations Research (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供了一种子母弹分离安全性评价方法及系统及子母弹分离仿真平台，包括：计算子母弹分离各个特征点的气动力系数；根据气动力系数构建子母弹分离变构型气动力响应面模型；构建子母弹分离过程中变构型变质量组合体多体动力学模型；通过多体动力学模型与气动力响应面模型耦合模拟方法，构建子母弹分离流固耦合动力学模型；对子母弹分离流固耦合动力学模型不同的初始边界条件进行仿真计算，获得子母弹分离过程的安全包络，对子母弹分离安全性进行动态评估，获得最优的子母弹分离初始边界条件。本发明解决了本领域对子母弹高速分离安全性评价和过程控制的技术难题；为制定有利于提高子母弹高速分离安全性的技术措施提供量化的设计依据。

Description

子母弹分离安全性评价方法及系统及子母弹分离仿真平台
技术领域
本发明涉及飞行器分离装置总体设计技术领域，具体地，涉及一种子母弹分离安全性评价方法及系统及子母弹分离仿真平台。
背景技术
子母弹用于对付集群目标，作为大纵深、大面积的压制兵器，已广泛地用于炮弹、航弹、火箭弹和导弹上。子母弹分离气动特性是子母弹研制中的关键技术之一，分离时刻母弹的飞行状态(飞行马赫数、攻角、侧滑角等)，以及子弹的抛射速度、排列及抛射方式等，都对子弹群落点分布的形状、抛撒面积大小及其均匀程度存在影响，并进而影响子母弹毁伤目标的效果。
采用从尾到头的逆向子母弹高速分离技术时，子弹在弹射器弹射力作用下，沿着与母弹固连的分离装置上的导轨运动，后、中、前三吊挂先后与分离装置上的导轨分离。当后、中吊挂脱离导轨后，仅剩前吊挂在导轨内，此时子弹的运动轨迹和姿态受到母弹的约束和牵连，处于半约束状态，分离动力学模拟困难。其特点是：(1)多学科耦合。空气动力学、飞行力学和多体动力学之间高度耦合。(2)接触非连续。子弹与母弹构成的两体或多体系统之间存在非连续接触(如母体的滑轨与分离体的滑块之间的接触)。(3)变构型变质量。分离过程中子弹与母弹接触时组合系统的外形和质量特性(质心、转动惯量等)存在时变现象。(4)分离瞬时扰动大。母弹与子弹质量比小(相比飞机导弹质量比小一个量级)，分离瞬时子弹对母弹存在很大冲击扰动。(5)离轨条件难确定。分离过程中母弹无动力、无驾驶员在回路进行稳定控制，高度、速度和姿态等难以保持，而子弹与母弹做牵连运动，导致子弹的离轨条件散布大，相比飞机投弹存在更大不确定性。因此子母弹高速分离过程非常复杂，具有极大的安全隐患。分离过程中可能会出现后吊挂卡滞、子弹与母弹外包络碰撞、子弹分离后初始扰动过大等问题。为了为子母弹高速分离系统设计提供重要的理论支撑依据，确保子母弹安全分离，需构建子母弹高速分离安全性评价系统及方法对子母弹高速分离过程进行评价。
子母弹分离性安全性评价的专利目前没有类似专利可以查询到。
发明内容
针对现有技术中的缺陷，本发明的目的是提供一种子母弹分离安全性评价方法及系统及子母弹分离仿真平台。
根据本发明提供的一种子母弹分离安全性评价方法，包括：
步骤S1：计算子母弹分离各个特征点的气动力系数；
步骤S2：根据气动力系数构建子母弹分离变构型气动力响应面模型；
步骤S3：构建子母弹分离过程中变构型变质量组合体多体动力学模型；
步骤S4：通过多体动力学模型与气动力响应面模型耦合模拟，构建子母弹分离流固耦合动力学模型；
步骤S5：对子母弹分离流固耦合动力学模型不同的初始边界条件进行仿真计算，获得子母弹分离过程的安全包络，对子母弹分离安全性进行动态评估，获得最优的子母弹分离初始边界条件。
优选地：
所述的子母弹分离过程中各个特征点气动力系数，应用计算流体力学软件Fluent计算，揭示变构型子母弹组合体的非定常气动特征；
所述的子母弹分离时，其周围的流场为三维粘性非定常的湍流流场，湍流流场采用分离涡模拟方法模拟，其控制方程的运输方程形式为：
Figure GDA0003590407020000021
其中，ρ为空气密度，
Figure GDA0003590407020000022
为流场通量，t为时间，div()为散度，u为速度矢量，div(grad())表示梯度的散度，Γ为扩散系数，S为源项；
所述的分离涡模拟方法，在薄边界层中分离涡模拟使用RANS湍流模型，RANS为雷诺时均模拟；在非薄边界层区域使用LES方法，分离涡模拟方法在不同区域将两种模型各自的优点结合起来，LES为大涡模拟方法。
优选地，所述步骤S2包括：根据计算获得的气动力系数应用三次样条插值方法构建子母弹分离变构型气动力响应面模型，为子母弹分离流固耦合动力学分析奠定基础；
三次样条插值方法具体步骤如下：
已知气动力响应面模型S(t)满足插值条件yi＝f(xi,ti)，i＝1,2,3...,n,并且有S(ti)＝yi，i＝1,2,...,n及其相应的边界条件为：
(1)第一类边界条件：S'(t1)＝f1'，S'(tn)＝f'n；
(2)第二类边界条件：S”(t1)＝f1”，S”(tn)＝f”n；
yi表示第i个插值节点处的插值函数序列值；
xi表示第i个插值节点处的自变量序列值；
ti表示第i个插值节点；
n表示插值节点i的最大值；
S(ti)表示第i个插值节点处的气动力序列值；
S'(t1)表示第1个插值节点处的一阶导数；
f1'表示插值函数第1个插值节点处的一阶导数；
S'(tn)表示第n个插值节点处的一阶导数；
f'n表示插值函数第n个插值节点处的一阶导数；
S”(t1)表示第1个插值节点处的二阶导数；
f1”表示插值函数第1个插值节点处的二阶导数；
S”(tn)表示第n个插值节点处的二阶导数；
f”n表示插值函数第n个插值节点处的二阶导数。
利用三次样条函数的二阶导数值S”(ti)＝mi，i＝1,2,3...,n表达样条函数S(t)，由于在区间[ti,ti+1]上是次数不超过三次的多项式，因此在区间上是线性的，具体表示为
Figure GDA0003590407020000031
其中，ti为插值节点，hi＝ti+1-ti；
S”(ti)表示第i个插值节点处的二阶导数；
mi表示插值节点t处的二阶导数；
S”(t)表示插值节点t处的二阶导数；
对上式进行两次积分并利用插值条件，得积分常数，于是得到三次样条函数为：
Figure GDA0003590407020000041
采用第一类边界条件和第二类边界条件求得mi，i＝1,2,3...,n，由于在区间[ti-1,ti]和[ti,ti+1]上的表达式不同，有要保证在节点ti处的光滑性，即ti处左右导数相等，Si'(ti-0)＝S'i+1(ti+0)，S'i(ti-0)表示插值节点i处的一阶左导数，S'i+1(ti+0)表示插值节点i处的一阶右导数；由上述条件得到关于mi的三对角方程组，由于方程组的系数矩阵是对角占优矩阵，所以所求方程有唯一解，用追赶法求出mi，代入上式即得到所求的三次样条插值函数S(t)，即所需的气动力响应面模型。
优选地，所述步骤S3包括：
考虑非线性气动载荷和非线性非连续接触碰撞的因素，采用牛顿欧拉矢量力学的方法构建分离过程中变构型、变质量组合体多体动力学模型；
所述的牛顿欧拉矢量力学的方法建模过程为：以变构型变质量子母弹分离机构中每一个物体为单元，建立固结在刚体上的坐标系，刚体的位置相对于一个公共参考基进行定义，其位置坐标统一为刚体坐标系基点的笛卡尔坐标与坐标系的方位坐标，方位坐标选用欧拉角或欧拉参数；单个物体位置坐标在三维系统中为6个；对于由N个刚体组成的变构型变质量子母弹分离机构，位置坐标阵q中的坐标个数为6N，由于铰约束的存在，这些位置坐标不独立；
完整约束条件下变构型变质量子母弹分离的多体动力学微分-代数方程表示为：
Figure GDA0003590407020000042
其中，q、
Figure GDA0003590407020000043
分别为多体系统的广义坐标、广义坐标对时间的一阶倒数、广义坐标对时间的二阶导数，M(q,t)为多体系统的广义质量矩阵，Φ(q,t)和Φq(q,t)为多体系统的约束函数向量及其对广义坐标q的雅克比矩阵，λ为约束的拉格朗日乘子，
Figure GDA0003590407020000044
为广义外力向量，包含燃气弹射力Fg、法向接触力FN、切向接触力FT和气动力
Figure GDA0003590407020000045
的广义力。
优选地，所述步骤S4包括：
建立含非连续接触碰撞的多体动力学模型与气动力响应面模型耦合模拟方法，构建多体非连续接触动态分离安全性评估的时空数值仿真平台，为分离安全性研究奠定基础；
所述的子母弹分离流固耦合动力学模型的构建基本思路是：利用接口程序连接子母弹分离变构型气动力响应面模型求解器和子母弹分离变构型变质量多体动力学模型求解器，气动力响应面模型求解器实时提供子母弹在当前姿态下的气动载荷，多体动力学模型求解器实时提供子母弹的分离姿态数据，构建多体非连续接触动态分离安全性评估的时空数值仿真平台。
优选地，所述步骤S5包括：
采用蒙特卡洛方法赋予子母弹分离过程流固耦合动力学模型不同的初始边界条件，利用构建的时空数值仿真平台中的接口程序进行仿真计算；接口程序将当前气动载荷传递给考虑非线性非连续接触的多体动力学模型求解器，多体动力学模型求解器计算出当前气动载荷作用下子母弹的姿态数据；接口程序再将子母弹姿态数据传递给气动力响应面模型求解器用于子母弹气动力计算；如此反复迭代完成子母弹分离过程的流固耦合仿真计算；
评估内容包含分离过程中是否会发生卡滞、子弹与母弹外包络是否发生碰撞、子弹分离后的姿态角和姿态角速度是否满足控制系统要求。
根据本发明提供的一种子母弹分离安全性评价系统，包括：
模块M1：计算子母弹分离各个特征点的气动力系数；
模块M2：根据气动力系数构建子母弹分离变构型气动力响应面模型；
模块M3：构建子母弹分离过程中变构型变质量组合体多体动力学模型；
模块M4：通过多体动力学模型与气动力响应面模型耦合模拟，构建子母弹分离流固耦合动力学模型；
模块M5：对子母弹分离流固耦合动力学模型不同的初始边界条件进行仿真计算，获得子母弹分离过程的安全包络，对子母弹分离安全性进行动态评估，获得最优的子母弹分离初始边界条件。
优选地，所述的子母弹分离过程中各个特征点气动力系数，应用计算流体力学软件Fluent计算，揭示变构型子母弹组合体的非定常气动特征；
所述的子母弹分离时，其周围的流场为三维粘性非定常的湍流流场，湍流流场采用分离涡模拟方法模拟，其控制方程的运输方程形式为：
Figure GDA0003590407020000051
其中，ρ为空气密度，
Figure GDA0003590407020000052
为流场通量，t为时间，div()为散度，u为速度矢量，div(grad())表示梯度的散度，Γ为扩散系数，S为源项；
所述的分离涡模拟方法，在薄边界层中分离涡模拟使用RANS湍流模型，RANS为雷诺时均模拟；在非薄边界层区域使用LES方法，分离涡模拟方法在不同区域将两种模型各自的优点结合起来，LES为大涡模拟方法。
优选地，所述模块M2包括：
根据计算获得的气动力系数应用三次样条插值方法构建子母弹分离变构型气动力响应面模型，为子母弹分离流固耦合动力学分析奠定基础；
三次样条插值方法具体步骤如下：
已知气动力响应面模型S(t)满足插值条件yi＝f(xi,ti)，i＝1,2,3...,n,并且有S(ti)＝yi，i＝1,2,...,n及其相应的边界条件为：
(1)第一类边界条件：S'(t1)＝f1'，S'(tn)＝f'n；
(2)第二类边界条件：S”(t1)＝f1”，S”(tn)＝f”n；
yi表示第i个插值节点处的插值函数序列值；
xi表示第i个插值节点处的自变量序列值；
ti表示第i个插值节点；
n表示插值节点i的最大值；
S(ti)表示第i个插值节点处的气动力序列值；
S'(t1)表示第1个插值节点处的一阶导数；
f1'表示插值函数第1个插值节点处的一阶导数；
S'(tn)表示第n个插值节点处的一阶导数；
f'n表示插值函数第n个插值节点处的一阶导数；
S”(t1)表示第1个插值节点处的二阶导数；
f1”表示插值函数第1个插值节点处的二阶导数；
S”(tn)表示第n个插值节点处的二阶导数；
f”n表示插值函数第n个插值节点处的二阶导数。
利用三次样条函数的二阶导数值S”(ti)＝mi，i＝1,2,3...,n表达样条函数S(t)，由于在区间[ti,ti+1]上是次数不超过三次的多项式，因此在区间上是线性的，具体表示为
Figure GDA0003590407020000061
其中，ti为插值节点，hi＝ti+1-ti；
S”(ti)表示第i个插值节点处的二阶导数；
mi表示插值节点t处的二阶导数；
S”(t)表示插值节点t处的二阶导数；
对上式进行两次积分并利用插值条件，得积分常数，于是得到三次样条函数为：
Figure GDA0003590407020000071
采用第一类边界条件和第二类边界条件求得mi，i＝1,2,3...,n，由于在区间[ti-1,ti]和[ti,ti+1]上的表达式不同，有要保证在节点ti处的光滑性，即ti处左右导数相等，S'i(ti-0)＝S'i+1(ti+0)，S'i(ti-0)表示插值节点i处的一阶左导数，S'i+1(ti+0)表示插值节点i处的一阶右导数；由上述条件得到关于mi的三对角方程组，由于方程组的系数矩阵是对角占优矩阵，所以所求方程有唯一解，用追赶法求出mi，代入上式即得到所求的三次样条插值函数S(t)，即所需的气动力响应面模型；
优选地，在所述模块M3中：：
考虑非线性气动载荷和非线性非连续接触碰撞的因素，采用牛顿欧拉矢量力学的方法构建分离过程中变构型、变质量组合体多体动力学模型；
所述的牛顿欧拉矢量力学的方法建模过程为：以变构型变质量子母弹分离机构中每一个物体为单元，建立固结在刚体上的坐标系，刚体的位置相对于一个公共参考基进行定义，其位置坐标统一为刚体坐标系基点的笛卡尔坐标与坐标系的方位坐标，方位坐标选用欧拉角或欧拉参数；单个物体位置坐标在三维系统中为6个；对于由N个刚体组成的变构型变质量子母弹分离机构，位置坐标阵q中的坐标个数为6N，由于铰约束的存在，这些位置坐标不独立；
完整约束条件下变构型变质量子母弹分离的多体动力学微分-代数方程表示为：
Figure GDA0003590407020000072
其中，q、
Figure GDA0003590407020000073
分别为多体系统的广义坐标、广义坐标对时间的一阶倒数、广义坐标对时间的二阶导数，M(q,t)为多体系统的广义质量矩阵，Φ(q,t)和Φq(q,t)为多体系统的约束函数向量及其对广义坐标q的雅克比矩阵，λ为约束的拉格朗日乘子，
Figure GDA0003590407020000074
为广义外力向量，包含燃气弹射力Fg、法向接触力FN、切向接触力FT和气动力
Figure GDA0003590407020000075
的广义力；
优选地，在所述模块M4中：
建立含非连续接触碰撞的多体动力学模型与气动力响应面模型耦合模拟方法，构建多体非连续接触动态分离安全性评估的时空数值仿真平台，为分离安全性研究奠定基础；
所述的子母弹分离流固耦合动力学模型的构建基本思路是：利用接口程序连接子母弹分离变构型气动力响应面模型求解器和子母弹分离变构型变质量多体动力学模型求解器，气动力响应面模型求解器实时提供子母弹在当前姿态下的气动载荷，多体动力学模型求解器实时提供子母弹的分离姿态数据，构建多体非连续接触动态分离安全性评估的时空数值仿真平台；
优选地，在所述模块M5中：
采用蒙特卡洛方法赋予子母弹分离过程流固耦合动力学模型不同的初始边界条件，利用构建的时空数值仿真平台中的接口程序进行仿真计算；接口程序将当前气动载荷传递给考虑非线性非连续接触的多体动力学模型求解器，多体动力学模型求解器计算出当前气动载荷作用下子母弹的姿态数据；接口程序再将子母弹姿态数据传递给气动力响应面模型求解器用于子母弹气动力计算；如此反复迭代完成子母弹分离过程的流固耦合仿真计算；
评估内容包含分离过程中是否会发生卡滞、子弹与母弹外包络是否发生碰撞、子弹分离后的姿态角和姿态角速度是否满足控制系统要求。
根据本发明提供的一种子母弹分离仿真平台，采用所述的子母弹分离安全性评价方法，或者所述的子母弹分离安全性评价系统，对子母弹高速分离进行仿真，得到安全性的仿真结果并展示给用户。
与现有技术相比，本发明具有如下的有益效果：
1、本发明解决了本领域对子母弹高速分离安全性评价和过程控制的技术难题；
2、本发明为制定有利于提高子母弹高速分离安全性的技术措施提供量化的设计依据；
3、本发明可为子母弹高速分离安全性设计提供理论依据，具有重要的工程价值。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
图1为本发明的子母弹高速分离安全性评价系统组成图；
图2为本发明的子母弹高速分离安全性评价方法步骤图；
图3为本发明实施例中在Matlab/Simulink中封装好的三次样条插值子弹气动力响应面模型；
图4为本发明实施例中在Matlab/Simulink中封装好的三次样条插值母弹气动力响应面模型；
图5为本发明实施例中在Matlab/Simulink中封装好的子母弹高速分离变构型变质量多体动力学模型图；
图6为本发明实施例中在Matlab/Simulink中封装好的子母弹高速分离流固耦合模型；
图7为本发明实施例中的子母弹高速分离流固耦合仿真模型的输入输出关系；
图8为本发明实施例中的气动力插值模型和子母弹高速分离多体动力学模型间实时耦合计算过程图。
图9为本发明实施例中的多工况子母弹高速分离过程子弹俯仰角速度、偏航角速度和滚转角速度曲线。
图中：
pitch angular velocity为俯仰角速度；
yaw angular velocity为偏航角速度；
roll angular velocity为滚转角速度。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变化和改进。这些都属于本发明的保护范围。
实施例1
根据本发明提供的一种子母弹分离安全性评价方法，包括：
步骤S1：计算子母弹分离各个特征点的气动力系数；
步骤S2：根据气动力系数构建子母弹分离变构型气动力响应面模型；
步骤S3：构建子母弹分离过程中变构型变质量组合体多体动力学模型；
步骤S4：通过多体动力学模型与气动力响应面模型耦合模拟，构建子母弹分离流固耦合动力学模型；
步骤S5：对子母弹分离流固耦合动力学模型不同的初始边界条件进行仿真计算，获得子母弹分离过程的安全包络，对子母弹分离安全性进行动态评估，获得最优的子母弹分离初始边界条件。
具体地：
所述的子母弹分离过程中各个特征点气动力系数，应用计算流体力学软件Fluent计算，揭示变构型子母弹组合体的非定常气动特征；
所述的子母弹分离时，其周围的流场为三维粘性非定常的湍流流场，湍流流场采用分离涡模拟方法模拟，其控制方程的运输方程形式为：
Figure GDA0003590407020000101
其中，ρ为空气密度，
Figure GDA0003590407020000102
为流场通量，t为时间，div()为散度，u为速度矢量，div(grad())表示梯度的散度，Γ为扩散系数，S为源项；
所述的分离涡模拟方法，在薄边界层中分离涡模拟使用RANS湍流模型，RANS为雷诺时均模拟；在非薄边界层区域使用LES方法，分离涡模拟方法在不同区域将两种模型各自的优点结合起来，LES为大涡模拟方法。
具体地，所述步骤S2包括：根据计算获得的气动力系数应用三次样条插值方法构建子母弹分离变构型气动力响应面模型，为子母弹分离流固耦合动力学分析奠定基础；
三次样条插值方法具体步骤如下：
已知气动力响应面模型S(t)满足插值条件yi＝f(xi,ti)，i＝1,2,3...,n,并且有S(ti)＝yi，i＝1,2,...,n及其相应的边界条件为：
(1)第一类边界条件：S'(t1)＝f1'，S'(tn)＝f'n；
(2)第二类边界条件：S”(t1)＝f1”，S”(tn)＝f”n；
yi表示第i个插值节点处的插值函数序列值；
xi表示第i个插值节点处的自变量序列值；
ti表示第i个插值节点；
n表示插值节点i的最大值；
S(ti)表示第i个插值节点处的气动力序列值；
S'(t1)表示第1个插值节点处的一阶导数；
f1'表示插值函数第1个插值节点处的一阶导数；
S'(tn)表示第n个插值节点处的一阶导数；
f'n表示插值函数第n个插值节点处的一阶导数；
S”(t1)表示第1个插值节点处的二阶导数；
f1”表示插值函数第1个插值节点处的二阶导数；
S”(tn)表示第n个插值节点处的二阶导数；
f”n表示插值函数第n个插值节点处的二阶导数。
利用三次样条函数的二阶导数值S”(ti)＝mi，i＝1,2,3...,n表达样条函数S(t)，由于在区间[ti,ti+1]上是次数不超过三次的多项式，因此在区间上是线性的，具体表示为
Figure GDA0003590407020000111
其中，ti为插值节点，hi＝ti+1-ti；
S”(ti)表示第i个插值节点处的二阶导数；
mi表示插值节点t处的二阶导数；
S”(t)表示插值节点t处的二阶导数；
对上式进行两次积分并利用插值条件，得积分常数，于是得到三次样条函数为：
Figure GDA0003590407020000112
采用第一类边界条件和第二类边界条件求得mi，i＝1,2,3...,n，由于在区间[ti-1,ti]和[ti,ti+1]上的表达式不同，有要保证在节点ti处的光滑性，即ti处左右导数相等，S'i(ti-0)＝S'i+1(ti+0)，S'i(ti-0)表示插值节点i处的一阶左导数，S'i+1(ti+0)表示插值节点i处的一阶右导数；由上述条件得到关于mi的三对角方程组，由于方程组的系数矩阵是对角占优矩阵，所以所求方程有唯一解，用追赶法求出mi，代入上式即得到所求的三次样条插值函数S(t)，即所需的气动力响应面模型。
具体地，所述步骤S3包括：
考虑非线性气动载荷和非线性非连续接触碰撞的因素，采用牛顿欧拉矢量力学的方法构建分离过程中变构型、变质量组合体多体动力学模型；
所述的牛顿欧拉矢量力学的方法建模过程为：以变构型变质量子母弹分离机构中每一个物体为单元，建立固结在刚体上的坐标系，刚体的位置相对于一个公共参考基进行定义，其位置坐标统一为刚体坐标系基点的笛卡尔坐标与坐标系的方位坐标，方位坐标选用欧拉角或欧拉参数；单个物体位置坐标在三维系统中为6个；对于由N个刚体组成的变构型变质量子母弹分离机构，位置坐标阵q中的坐标个数为6N，由于铰约束的存在，这些位置坐标不独立；
完整约束条件下变构型变质量子母弹分离的多体动力学微分-代数方程表示为：
Figure GDA0003590407020000121
其中，q、
Figure GDA0003590407020000122
分别为多体系统的广义坐标、广义坐标对时间的一阶倒数、广义坐标对时间的二阶导数，M(q,t)为多体系统的广义质量矩阵，Φ(q,t)和Φq(q,t)为多体系统的约束函数向量及其对广义坐标q的雅克比矩阵，λ为约束的拉格朗日乘子，
Figure GDA0003590407020000123
为广义外力向量，包含燃气弹射力Fg、法向接触力FN、切向接触力FT和气动力
Figure GDA0003590407020000124
的广义力。
具体地，所述步骤S4包括：
建立含非连续接触碰撞的多体动力学模型与气动力响应面模型耦合模拟方法，构建多体非连续接触动态分离安全性评估的时空数值仿真平台，为分离安全性研究奠定基础；
所述的子母弹分离流固耦合动力学模型的构建基本思路是：利用接口程序连接子母弹分离变构型气动力响应面模型求解器和子母弹分离变构型变质量多体动力学模型求解器，气动力响应面模型求解器实时提供子母弹在当前姿态下的气动载荷，多体动力学模型求解器实时提供子母弹的分离姿态数据，构建多体非连续接触动态分离安全性评估的时空数值仿真平台。
具体地，所述步骤S5包括：
采用蒙特卡洛方法赋予子母弹分离过程流固耦合动力学模型不同的初始边界条件，利用构建的时空数值仿真平台中的接口程序进行仿真计算；接口程序将当前气动载荷传递给考虑非线性非连续接触的多体动力学模型求解器，多体动力学模型求解器计算出当前气动载荷作用下子母弹的姿态数据；接口程序再将子母弹姿态数据传递给气动力响应面模型求解器用于子母弹气动力计算；如此反复迭代完成子母弹分离过程的流固耦合仿真计算；
评估内容包含分离过程中是否会发生卡滞、子弹与母弹外包络是否发生碰撞、子弹分离后的姿态角和姿态角速度是否满足控制系统要求。
根据本发明提供的一种子母弹分离安全性评价系统，其特征在于，包括：
模块M1：计算子母弹分离各个特征点的气动力系数；
模块M2：根据气动力系数构建子母弹分离变构型气动力响应面模型；
模块M3：构建子母弹分离过程中变构型变质量组合体多体动力学模型；
模块M4：通过多体动力学模型与气动力响应面模型耦合模拟，构建子母弹分离流固耦合动力学模型；
模块M5：对子母弹分离流固耦合动力学模型不同的初始边界条件进行仿真计算，获得子母弹分离过程的安全包络，对子母弹分离安全性进行动态评估，获得最优的子母弹分离初始边界条件。
具体地，所述的子母弹分离过程中各个特征点气动力系数，应用计算流体力学软件Fluent计算，揭示变构型子母弹组合体的非定常气动特征；
所述的子母弹分离时，其周围的流场为三维粘性非定常的湍流流场，湍流流场采用分离涡模拟方法模拟，其控制方程的运输方程形式为：
Figure GDA0003590407020000131
其中，ρ为空气密度，
Figure GDA0003590407020000132
为流场通量，t为时间，div()为散度，u为速度矢量，div(grad())表示梯度的散度，Γ为扩散系数，S为源项；
所述的分离涡模拟方法，在薄边界层中分离涡模拟使用RANS湍流模型，RANS为雷诺时均模拟；在非薄边界层区域使用LES方法，分离涡模拟方法在不同区域将两种模型各自的优点结合起来，LES为大涡模拟方法。
具体地，所述模块M2包括：
根据计算获得的气动力系数应用三次样条插值方法构建子母弹分离变构型气动力响应面模型，为子母弹分离流固耦合动力学分析奠定基础；
三次样条插值方法具体步骤如下：
已知气动力响应面模型S(t)满足插值条件yi＝f(xi,ti)，i＝1,2,3...,n,并且有S(ti)＝yi，i＝1,2,...,n及其相应的边界条件为：
(1)第一类边界条件：S'(t1)＝f'1，S'(tn)＝f'n；
(2)第二类边界条件：S”(t1)＝f1”，S”(tn)＝f”n；
yi表示第i个插值节点处的插值函数序列值；
xi表示第i个插值节点处的自变量序列值；
ti表示第i个插值节点；
n表示插值节点i的最大值；
S(ti)表示第i个插值节点处的气动力序列值；
S'(t1)表示第1个插值节点处的一阶导数；
f1'表示插值函数第1个插值节点处的一阶导数；
S'(tn)表示第n个插值节点处的一阶导数；
f'n表示插值函数第n个插值节点处的一阶导数；
S”(t1)表示第1个插值节点处的二阶导数；
f1”表示插值函数第1个插值节点处的二阶导数；
S”(tn)表示第n个插值节点处的二阶导数；
f”n表示插值函数第n个插值节点处的二阶导数。
利用三次样条函数的二阶导数值S”(ti)＝mi，i＝1,2,3...,n表达样条函数S(t)，由于在区间[ti,ti+1]上是次数不超过三次的多项式，因此在区间上是线性的，具体表示为
Figure GDA0003590407020000141
其中，ti为插值节点，hi＝ti+1-ti；
S”(ti)表示第i个插值节点处的二阶导数；
mi表示插值节点t处的二阶导数；
S”(t)表示插值节点t处的二阶导数；
对上式进行两次积分并利用插值条件，得积分常数，于是得到三次样条函数为：
Figure GDA0003590407020000142
采用第一类边界条件和第二类边界条件求得mi，i＝1,2,3...,n，由于在区间[ti-1,ti]和[ti,ti+1]上的表达式不同，有要保证在节点ti处的光滑性，即ti处左右导数相等，S'i(ti-0)＝S'i+1(ti+0)，S'i(ti-0)表示插值节点i处的一阶左导数，S'i+1(ti+0)表示插值节点i处的一阶右导数；由上述条件得到关于mi的三对角方程组，由于方程组的系数矩阵是对角占优矩阵，所以所求方程有唯一解，用追赶法求出mi，代入上式即得到所求的三次样条插值函数S(t)，即所需的气动力响应面模型；
具体地，在所述模块M3中：：
考虑非线性气动载荷和非线性非连续接触碰撞的因素，采用牛顿欧拉矢量力学的方法构建分离过程中变构型、变质量组合体多体动力学模型；
所述的牛顿欧拉矢量力学的方法建模过程为：以变构型变质量子母弹分离机构中每一个物体为单元，建立固结在刚体上的坐标系，刚体的位置相对于一个公共参考基进行定义，其位置坐标统一为刚体坐标系基点的笛卡尔坐标与坐标系的方位坐标，方位坐标选用欧拉角或欧拉参数；单个物体位置坐标在三维系统中为6个；对于由N个刚体组成的变构型变质量子母弹分离机构，位置坐标阵q中的坐标个数为6N，由于铰约束的存在，这些位置坐标不独立；
完整约束条件下变构型变质量子母弹分离的多体动力学微分-代数方程表示为：
Figure GDA0003590407020000151
其中，q、
Figure GDA0003590407020000152
分别为多体系统的广义坐标、广义坐标对时间的一阶倒数、广义坐标对时间的二阶导数，M(q,t)为多体系统的广义质量矩阵，Φ(q,t)和Φq(q,t)为多体系统的约束函数向量及其对广义坐标q的雅克比矩阵，λ为约束的拉格朗日乘子，
Figure GDA0003590407020000153
为广义外力向量，包含燃气弹射力Fg、法向接触力FN、切向接触力FT和气动力
Figure GDA0003590407020000154
的广义力；
具体地，在所述模块M4中：
建立含非连续接触碰撞的多体动力学模型与气动力响应面模型耦合模拟方法，构建多体非连续接触动态分离安全性评估的时空数值仿真平台，为分离安全性研究奠定基础；
所述的子母弹分离流固耦合动力学模型的构建基本思路是：利用接口程序连接子母弹分离变构型气动力响应面模型求解器和子母弹分离变构型变质量多体动力学模型求解器，气动力响应面模型求解器实时提供子母弹在当前姿态下的气动载荷，多体动力学模型求解器实时提供子母弹的分离姿态数据，构建多体非连续接触动态分离安全性评估的时空数值仿真平台；
具体地，在所述模块M5中：
采用蒙特卡洛方法赋予子母弹分离过程流固耦合动力学模型不同的初始边界条件，利用构建的时空数值仿真平台中的接口程序进行仿真计算；接口程序将当前气动载荷传递给考虑非线性非连续接触的多体动力学模型求解器，多体动力学模型求解器计算出当前气动载荷作用下子母弹的姿态数据；接口程序再将子母弹姿态数据传递给气动力响应面模型求解器用于子母弹气动力计算；如此反复迭代完成子母弹分离过程的流固耦合仿真计算；
评估内容包含分离过程中是否会发生卡滞、子弹与母弹外包络是否发生碰撞、子弹分离后的姿态角和姿态角速度是否满足控制系统要求。
根据本发明提供的一种子母弹分离仿真平台，采用上述所述的子母弹分离安全性评价方法，或者上述所述的子母弹分离安全性评价系统，对子母弹高速分离进行仿真，得到安全性的仿真结果并展示给用户。
实施例2
本领域技术人员可以将本实施例理解为实施例1的优选例。
本发明首先通过建立子母弹高速分离变构型非定常气动力模型并通过计算获得子母弹高速分离过程中的子弹和母弹上的气动力数据；其次，采用三次样条插值方法构建子母弹高速分离变构型气动力插值响应面模型；第三，采用牛顿-欧拉方法构建子母弹高速分离变构型变质量多体动力学模型；第四，综合变构型气动力插值响应面模型和变构型变质量多体动力学模型构建子母弹高速分离过程流固耦合动力学模型；第五，采用蒙特卡洛方法赋予子母弹高速分离过程流固耦合动力学模型不同的初始边界条件进行仿真计算，获得子母弹高速分离过程的安全包络，对子母弹高速分离安全性进行动态评估，最终获得最优的子母弹高速分离初始边界条件。
如图1所示，本发明实施例提供了一种子母弹高速分离安全性评价系统，包括子母弹高速分离变构型非定常气动力计算单元、子母弹高速分离变构型气动力响应面模型构建单元、子母弹高速分离变构型变质量多体动力学模型构建单元、子母弹高速分离流固耦合动力学模型构建单元和子母弹高速分离安全性评估单元。
所述子母弹高速分离变构型非定常气动力计算单元应用计算流体力学软件计算子母弹高速分离各个特征点的气动力系数，揭示变构型子母弹组合体的非定常气动特征；
所述子母弹高速分离变构型气动力响应面模型构建单元根据计算获得的气动力系数应用三次样条插值方法构建子母弹高速分离变构型气动力响应面模型，为子母弹高速分离流固耦合动力学分析奠定基础；
所述子母弹高速分离变构型变质量多体动力学模型构建单元考虑非线性高精度气动载荷和非线性非连续接触碰撞等因素，构建分离过程中变构型、变质量组合体多体动力学模型，解决后退式分离过程中复杂牵连运动问题，准确表征小质量比母弹和子弹非连续接触变构型变质量后退拔出过程。
所述子母弹高速分离流固耦合动力学模型构建单元建立含非连续接触碰撞的多体动力学模型与气动力响应面模型耦合模拟方法，构建多体非连续接触动态分离安全性评估的时空高精度数值仿真平台，为分离安全性研究奠定基础。
所述子母弹高速分离安全性评估单元采用蒙特卡洛方法赋予子母弹高速分离过程流固耦合动力学模型不同的初始边界条件进行仿真计算，获得子母弹高速分离过程的安全包络，对子母弹高速分离安全性进行动态评估，最终获得最优的子母弹高速分离初始边界条件。评估内容包含分离过程中是否会发生卡滞、子弹与母弹外包络是否发生碰撞；子弹分离后的姿态角和姿态角速度是否满足控制系统要求等。
如图2所示，本发明还提供了一种对子母弹高速分离过程进行安全性评价的方法，其包括以下步骤：
S1、子母弹高速分离变构型非定常气动力计算：应用计算流体力学软件计算子母弹高速分离各个特征点的气动力系数，揭示变构型子母弹组合体的非定常气动特征；
S2、子母弹高速分离变构型气动力响应面模型构建：根据计算获得的气动力系数应用三次样条插值方法构建子母弹高速分离变构型气动力响应面模型，为子母弹高速分离流固耦合动力学分析奠定基础；
S3、子母弹高速分离变构型变质量多体动力学模型构建：考虑非线性高精度气动载荷和非线性非连续接触碰撞等因素，构建分离过程中变构型、变质量组合体多体动力学模型，解决后退式分离过程中复杂牵连运动问题，准确表征小质量比母弹和子弹非连续接触变构型变质量后退拔出过程；
S4、子母弹高速分离流固耦合动力学模型构建；建立含非连续接触碰撞的多体动力学模型与气动力响应面模型耦合模拟方法，构建多体非连续接触动态分离安全性评估的时空高精度数值仿真平台，为分离安全性研究奠定基础；
S5，子母弹高速分离安全性评估：采用蒙特卡洛方法赋予子母弹高速分离过程流固耦合动力学模型不同的初始边界条件进行仿真计算，获得子母弹高速分离过程的安全包络，对子母弹高速分离安全性进行动态评估，最终获得最优的子母弹高速分离初始边界条件。评估内容包含分离过程中是否会发生卡滞、子弹与母弹外包络是否发生碰撞；子弹分离后的姿态角和姿态角速度是否满足控制系统要求等。
下面结合图2以及具体实施例对本发明的方法的工作流程做进一步解释：
子母弹高速分离过程中各个特征点气动力系数采用计算流体力学软件Fluent来计算获得。
子母弹高速分离时，其周围的流场为三维粘性非定常的湍流流场，其控制方程的运输方程形式为：
Figure GDA0003590407020000181
其中，t为时间，ρ为空气密度，u为速度矢量，
Figure GDA0003590407020000182
为流场通量，Γ为扩散系数，S为源项。
子母弹高速分离时湍流流场采用分离涡模拟(Detached-Eddy Simulation，DES)方法模拟。
分离涡模拟(Detached-Eddy Simulation，DES)方法，能够捕捉到丰富的流动现象和非定常特性。在薄边界层中，DES使用RANS湍流模型；而在其它区域，使用LES方法。DES方法在不同区域将两种模型各自的优点结合起来。
子母弹高速分离变构型气动力响应面模型采用三次样条插值方法实现，具体的，
已知气动力响应面模型S(t)满足插值条件yi＝f(xi,ti)，(i＝1,2,3...,n)并且有S(ti)＝yi，i＝1,2,...,n及其相应的边界条件为：(1)第一类边界条件：S'(t1)＝f1'，S'(tn)＝f'n；(2)第二类边界条件：S”(t1)＝f1”，S”(tn)＝f”n。利用三次样条函数的二阶导数值S”(ti)＝mi，i＝1,2,3...,n表达样条函数S(t)，由于在区间[ti,ti+1]上是次数不超过三次的多项式，可知在区间上是线性的，具体表示为
Figure GDA0003590407020000183
对上式进行两次积分并利用插值条件，可得积分常数。于是得到三次样条函数为：
Figure GDA0003590407020000184
采用第一类边界条件和第二类边界条件求得mi(i＝1,2,3...,n)。由于在区间[ti-1,ti]和[ti,ti+1]上的表达式不同，有要保证在节点ti处的光滑性，即ti处左右导数相等，可知S'i(ti-0)＝S'i+1(ti+0)。由上述条件可得关于mi的三对角方程组，由于方程组的系数矩阵是对角占优矩阵，所以所求方程有唯一解，可用追赶法求出mi，代入上式即得到所求的三次样条插值函数S(t)，即所需的气动力响应面模型，图3为本发明实施例中在Matlab/Simulink中封装好的三次样条插值子弹气动力响应面模型；图4为本发明实施例中在Matlab/Simulink中封装好的三次样条插值母弹气动力响应面模型。
子母弹高速分离变构型变质量多体动力学模型采用牛顿-欧拉矢量力学的方法来构建。图5为本发明实施例中在Matlab/Simulink中封装好的子母弹高速分离变构型变质量多体动力学模型图。牛顿-欧拉矢量力学的方法建模过程为：以变构型变质量子母弹分离机构中每一个物体为单元，建立固结在刚体上的坐标系，刚体的位置相对于一个公共参考基进行定义，其位置坐标(也可称为广义坐标)统一为刚体坐标系基点的笛卡尔坐标与坐标系的方位坐标，方位坐标可以选用欧拉角或欧拉参数。单个物体位置坐标在三维系统中为6个。对于由N个刚体组成的变构型变质量子母弹分离机构，位置坐标阵q中的坐标个数为6N，由于铰约束的存在，这些位置坐标不独立。
完整约束条件下变构型变质量子母弹分离的多体动力学微分-代数方程可表示为：
Figure GDA0003590407020000195
其中，q、
Figure GDA0003590407020000191
和
Figure GDA0003590407020000192
分别为多体系统的广义坐标及其广义坐标对时间的一阶、二阶导数，M(q,t)为多体系统的广义质量矩阵，Φ(q,t)和Φq(q,t)为多体系统的约束函数向量及其对广义坐标q的雅克比矩阵，λ为约束的拉格朗日乘子，
Figure GDA0003590407020000193
为广义外力向量，包含燃气弹射力Fg、法向接触力FN、切向接触力FT和气动力
Figure GDA0003590407020000194
的广义力。
子母弹高速分离流固耦合动力学模型的构建基本思路是：利用接口程序联接子母弹高速分离变构型气动力响应面模型求解器和子母弹高速分离变构型变质量多体动力学模型求解器，气动力响应面模型求解器实时提供子母弹在当前姿态下的气动载荷，多体动力学模型求解器实时提供子母弹的分离姿态数据，构建多体非连续接触动态分离安全性评估的时空高精度数值仿真平台。图6为本发明实施例中在Matlab/Simulink中封装好的子母弹高速分离流固耦合模型；图7为本发明实施例中的子母弹高速分离流固耦合仿真模型的输入输出关系。
子母弹高速分离安全性评估基本过程为：采用蒙特卡洛方法赋予子母弹高速分离过程流固耦合动力学模型不同的初始边界条件，利用所构建的时空高精度数值仿真平台中的接口程序进行子母弹高速分离流固耦合仿真计算。
子母弹高速分离流固耦合仿真计算过程为：接口程序将当前气动载荷传递给考虑非线性非连续接触的多体动力学模型求解器，多体动力学模型求解器计算出当前气动载荷作用下子母弹的姿态等数据；接口程序再将子母弹姿态等数据传递给气动力响应面模型求解器用于子母弹气动力计算；如此反复迭代完成子母弹分离过程的流固耦合仿真计算。图8为本发明实施例中的气动力插值模型和子母弹高速分离多体动力学模型间实时耦合计算过程图。
子母弹高速分离流固耦合仿真计算结果将获得子母弹高速分离过程的安全包络，对子母弹高速分离安全性进行动态评估，最终获得最优的子母弹高速分离初始边界条件。评估内容包含分离过程中是否会发生卡滞、子弹与母弹外包络是否发生碰撞；子弹分离后的姿态角和姿态角速度是否满足控制系统要求等。图9为本发明实施例中的多工况子母弹高速分离过程子弹俯仰角速度、偏航角速度和滚转角速度曲线，可根据该仿真结果判断子弹的分离角速度是否满足控制系统的起控要求。
本领域技术人员知道，除了以纯计算机可读程序代码方式实现本发明提供的系统、装置及其各个模块以外，完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统、装置及其各个模块以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同程序。所以，本发明提供的系统、装置及其各个模块可以被认为是一种硬件部件，而对其内包括的用于实现各种程序的模块也可以视为硬件部件内的结构；也可以将用于实现各种功能的模块视为既可以是实现方法的软件程序又可以是硬件部件内的结构。
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本申请的实施例和实施例中的特征可以任意相互组合。

Claims (10)

1.一种子母弹分离安全性评价方法，其特征在于，包括：
步骤S1：计算子母弹分离各个特征点的气动力系数；
步骤S2：根据气动力系数构建子母弹分离变构型气动力响应面模型；
步骤S3：构建子母弹分离过程中变构型变质量组合体多体动力学模型；
步骤S4：通过多体动力学模型与气动力响应面模型耦合模拟，构建子母弹分离流固耦合动力学模型；
步骤S5：对子母弹分离流固耦合动力学模型不同的初始边界条件进行仿真计算，获得子母弹分离过程的安全包络，对子母弹分离安全性进行动态评估，获得最优的子母弹分离初始边界条件；
所述步骤S1采用：所述的子母弹分离过程中各个特征点气动力系数应用计算流体力学软件Fluent计算，揭示变构型子母弹组合体的非定常气动特征；
子母弹分离时，其周围的流场为三维粘性非定常的湍流流场，湍流流场采用分离涡模拟方法模拟；
所述的分离涡模拟方法，在薄边界层中分离涡模拟使用RANS湍流模型，RANS为雷诺时均模拟；在非薄边界层区域使用LES方法，LES为大涡模拟方法；
所述步骤S2包括：根据计算获得的气动力系数应用三次样条插值方法构建子母弹分离变构型气动力响应面模型；
所述步骤S3包括：考虑非线性气动载荷和非线性非连续接触碰撞的因素，采用牛顿欧拉矢量力学的方法构建分离过程中变构型、变质量组合体多体动力学模型；
所述步骤S4采用：所述的子母弹分离流固耦合动力学模型的构建基本思路是：利用接口程序连接子母弹分离变构型气动力响应面模型求解器和子母弹分离变构型变质量多体动力学模型求解器，气动力响应面模型求解器实时提供子母弹在当前姿态下的气动载荷，多体动力学模型求解器实时提供子母弹的分离姿态数据，构建多体非连续接触动态分离安全性评估的时空数值仿真平台；
所述步骤S5采用：采用蒙特卡洛方法赋予子母弹分离过程流固耦合动力学模型不同的初始边界条件，利用构建的时空数值仿真平台中的接口程序进行仿真计算；接口程序将当前气动载荷传递给考虑非线性非连续接触的多体动力学模型求解器，多体动力学模型求解器计算出当前气动载荷作用下子母弹的姿态数据；接口程序再将子母弹姿态数据传递给气动力响应面模型求解器用于子母弹气动力计算；如此反复迭代完成子母弹分离过程的流固耦合仿真计算。
2.根据权利要求1所述的子母弹分离安全性评价方法，其特征在于：所述的子母弹分离时，其周围的流场为三维粘性非定常的湍流流场，湍流流场采用分离涡模拟方法模拟，其控制方程的运输方程形式为：
Figure FDA0003590407010000021
其中，ρ为空气密度，
Figure FDA0003590407010000022
为流场通量，t为时间，div()为散度，u为速度矢量，div(grad())表示梯度的散度，Γ为扩散系数，S为源项；
3.根据权利要求1所述的子母弹分离安全性评价方法，其特征在于，所述三次样条插值方法具体步骤如下：
已知气动力响应面模型S(t)满足插值条件yi＝f(xi,ti)，i＝1,2,3...,n,并且有S(ti)＝yi，i＝1,2,...,n及其相应的边界条件为：
(1)第一类边界条件：S'(t1)＝f1'，...，S'(tn)＝fn'；
(2)第二类边界条件：S”(t1)＝f1”，...，S”(tn)＝f”n；
yi表示第i个插值节点处的插值函数序列值；
xi表示第i个插值节点处的自变量序列值；
ti表示第i个插值节点；
n表示插值节点i的最大值；
S(ti)表示第i个插值节点处的气动力序列值；
S'(t1)表示第1个插值节点处的一阶导数；
f1'表示插值函数第1个插值节点处的一阶导数；
S'(tn)表示第n个插值节点处的一阶导数；
fn'表示插值函数第n个插值节点处的一阶导数；
S”(t1)表示第1个插值节点处的二阶导数；
f1”表示插值函数第1个插值节点处的二阶导数；
S”(tn)表示第n个插值节点处的二阶导数；
f”n表示插值函数第n个插值节点处的二阶导数；
利用三次样条函数的二阶导数值S”(ti)＝mi，i＝1,2,3...,n表达样条函数S(t)，由于在区间[ti,ti+1]上是次数不超过三次的多项式，因此在区间上是线性的，具体表示为
Figure FDA0003590407010000031
其中，ti为插值节点，hi＝ti+1-ti；
S”(ti)表示第i个插值节点处的二阶导数；
mi表示插值节点t处的二阶导数；
S”(t)表示插值节点t处的二阶导数；
对上式进行两次积分并利用插值条件，得积分常数，于是得到三次样条函数为：
Figure FDA0003590407010000032
采用第一类边界条件和第二类边界条件求得mi，i＝1,2,3...,n，由于在区间[ti-1,ti]和[ti,ti+1]上的表达式不同，有要保证在节点ti处的光滑性，即ti处左右导数相等，S'i(ti-0)＝S'i+1(ti+0)，S'i(ti-0)表示插值节点i处的一阶左导数，S'i+1(ti+0)表示插值节点i处的一阶右导数；由上述条件得到关于mi的三对角方程组，由于方程组的系数矩阵是对角占优矩阵，所以所求方程有唯一解，用追赶法求出mi，代入上式即得到所求的三次样条插值函数S(t)，即所需的气动力响应面模型。
4.根据权利要求1所述的子母弹分离安全性评价方法，其特征在于，所述的牛顿欧拉矢量力学的方法构建分离过程中变构型、变质量组合体多体动力学模型为：以变构型变质量子母弹分离机构中每一个物体为单元，建立固结在刚体上的坐标系，刚体的位置相对于一个公共参考基进行定义，其位置坐标统一为刚体坐标系基点的笛卡尔坐标与坐标系的方位坐标，方位坐标选用欧拉角或欧拉参数；单个物体位置坐标在三维系统中为6个；对于由N个刚体组成的变构型变质量子母弹分离机构，位置坐标阵q中的坐标个数为6N，由于铰约束的存在，这些位置坐标不独立；
完整约束条件下变构型变质量子母弹分离的多体动力学微分-代数方程表示为：
Figure FDA0003590407010000033
其中，q、
Figure FDA0003590407010000034
分别为多体系统的广义坐标、广义坐标对时间的一阶倒数、广义坐标对时间的二阶导数，M(q,t)为多体系统的广义质量矩阵，Φ(q,t)和Φq(q,t)为多体系统的约束函数向量及其对广义坐标q的雅克比矩阵，λ为约束的拉格朗日乘子，
Figure FDA0003590407010000041
为广义外力向量，包含燃气弹射力Fg、法向接触力FN、切向接触力FT和气动力
Figure FDA0003590407010000042
的广义力。
5.根据权利要求1所述的子母弹分离安全性评价方法，其特征在于，所述步骤S4包括：建立含非连续接触碰撞的多体动力学模型与气动力响应面模型耦合模拟方法，构建多体非连续接触动态分离安全性评估的时空数值仿真平台。
6.根据权利要求1所述的子母弹分离安全性评价方法，其特征在于，所述步骤S5包括：评估内容包含分离过程中是否会发生卡滞、子弹与母弹外包络是否发生碰撞、子弹分离后的姿态角和姿态角速度是否满足控制系统要求。
7.一种子母弹分离安全性评价系统，其特征在于，包括：
模块M1：计算子母弹分离各个特征点的气动力系数；
模块M2：根据气动力系数构建子母弹分离变构型气动力响应面模型；
模块M3：构建子母弹分离过程中变构型变质量组合体多体动力学模型；
模块M4：通过多体动力学模型与气动力响应面模型耦合模拟，构建子母弹分离流固耦合动力学模型；
模块M5：对子母弹分离流固耦合动力学模型不同的初始边界条件进行仿真计算，获得子母弹分离过程的安全包络，对子母弹分离安全性进行动态评估，获得最优的子母弹分离初始边界条件；
所述步骤S1采用：所述的子母弹分离过程中各个特征点气动力系数应用计算流体力学软件Fluent计算，揭示变构型子母弹组合体的非定常气动特征；
子母弹分离时，其周围的流场为三维粘性非定常的湍流流场，湍流流场采用分离涡模拟方法模拟；
所述的分离涡模拟方法，在薄边界层中分离涡模拟使用RANS湍流模型，RANS为雷诺时均模拟；在非薄边界层区域使用LES方法，LES为大涡模拟方法；
所述步骤S2包括：根据计算获得的气动力系数应用三次样条插值方法构建子母弹分离变构型气动力响应面模型；
所述步骤S3包括：考虑非线性气动载荷和非线性非连续接触碰撞的因素，采用牛顿欧拉矢量力学的方法构建分离过程中变构型、变质量组合体多体动力学模型；
所述步骤S4采用：所述的子母弹分离流固耦合动力学模型的构建基本思路是：利用接口程序连接子母弹分离变构型气动力响应面模型求解器和子母弹分离变构型变质量多体动力学模型求解器，气动力响应面模型求解器实时提供子母弹在当前姿态下的气动载荷，多体动力学模型求解器实时提供子母弹的分离姿态数据，构建多体非连续接触动态分离安全性评估的时空数值仿真平台；
所述步骤S5采用：采用蒙特卡洛方法赋予子母弹分离过程流固耦合动力学模型不同的初始边界条件，利用构建的时空数值仿真平台中的接口程序进行仿真计算；接口程序将当前气动载荷传递给考虑非线性非连续接触的多体动力学模型求解器，多体动力学模型求解器计算出当前气动载荷作用下子母弹的姿态数据；接口程序再将子母弹姿态数据传递给气动力响应面模型求解器用于子母弹气动力计算；如此反复迭代完成子母弹分离过程的流固耦合仿真计算。
8.根据权利要求7所述的子母弹分离安全性评价系统，其特征在于：
所述的子母弹分离时，其周围的流场为三维粘性非定常的湍流流场，湍流流场采用分离涡模拟方法模拟，其控制方程的运输方程形式为：
Figure FDA0003590407010000051
其中，ρ为空气密度，
Figure FDA0003590407010000052
为流场通量，t为时间，div()为散度，u为速度矢量，div(grad())表示梯度的散度，Γ为扩散系数，S为源项；
所述的分离涡模拟方法，在薄边界层中分离涡模拟使用RANS湍流模型，RANS为雷诺时均模拟；在非薄边界层区域使用LES方法，LES为大涡模拟方法。
9.根据权利要求7所述的子母弹分离安全性评价系统，其特征在于：所述三次样条插值方法具体步骤如下：
已知气动力响应面模型S(t)满足插值条件yi＝f(xi,ti)，i＝1,2,3...,n,并且有S(ti)＝yi，i＝1,2,...,n及其相应的边界条件为：
(1)第一类边界条件：S'(t1)＝f1'，S'(tn)＝f'n；
(2)第二类边界条件：S”(t1)＝f1”，S”(tn)＝f”n；
yi表示第i个插值节点处的插值函数序列值；
xi表示第i个插值节点处的自变量序列值；
ti表示第i个插值节点；
n表示插值节点i的最大值；
S(ti)表示第i个插值节点处的气动力序列值；
S'(t1)表示第1个插值节点处的一阶导数；
f1'表示插值函数第1个插值节点处的一阶导数；
S'(tn)表示第n个插值节点处的一阶导数；
fn'表示插值函数第n个插值节点处的一阶导数；
S”(t1)表示第1个插值节点处的二阶导数；
f1”表示插值函数第1个插值节点处的二阶导数；
S”(tn)表示第n个插值节点处的二阶导数；
f”n表示插值函数第n个插值节点处的二阶导数；
利用三次样条函数的二阶导数值S”(ti)＝mi，i＝1,2,3...,n表达样条函数S(t)，由于在区间[ti,ti+1]上是次数不超过三次的多项式，因此在区间上是线性的，具体表示为
Figure FDA0003590407010000061
其中，ti为插值节点，hi＝ti+1-ti；
S”(ti)表示第i个插值节点处的二阶导数；
mi表示插值节点t处的二阶导数；
S”(t)表示插值节点t处的二阶导数；
对上式进行两次积分并利用插值条件，得积分常数，于是得到三次样条函数为：
Figure FDA0003590407010000062
采用第一类边界条件和第二类边界条件求得mi，i＝1,2,3...,n，由于在区间[ti-1,ti]和[ti,ti+1]上的表达式不同，有要保证在节点ti处的光滑性，即ti处左右导数相等，S'i(ti-0)＝S'i+1(ti+0)，S'i(ti-0)表示插值节点i处的一阶左导数，S'i+1(ti+0)表示插值节点i处的一阶右导数；由上述条件得到关于mi的三对角方程组，由于方程组的系数矩阵是对角占优矩阵，所以所求方程有唯一解，用追赶法求出mi，代入上式即得到所求的三次样条插值函数S(t)，即所需的气动力响应面模型；
所述的牛顿欧拉矢量力学的方法构建分离过程中变构型、变质量组合体多体动力学模型为：以变构型变质量子母弹分离机构中每一个物体为单元，建立固结在刚体上的坐标系，刚体的位置相对于一个公共参考基进行定义，其位置坐标统一为刚体坐标系基点的笛卡尔坐标与坐标系的方位坐标，方位坐标选用欧拉角或欧拉参数；单个物体位置坐标在三维系统中为6个；对于由N个刚体组成的变构型变质量子母弹分离机构，位置坐标阵q中的坐标个数为6N，由于铰约束的存在，这些位置坐标不独立；
完整约束条件下变构型变质量子母弹分离的多体动力学微分-代数方程表示为：
Figure FDA0003590407010000071
其中，q、
Figure FDA0003590407010000072
分别为多体系统的广义坐标、广义坐标对时间的一阶倒数、广义坐标对时间的二阶导数，M(q,t)为多体系统的广义质量矩阵，Φ(q,t)和Φq(q,t)为多体系统的约束函数向量及其对广义坐标q的雅克比矩阵，λ为约束的拉格朗日乘子，
Figure FDA0003590407010000073
为广义外力向量，包含燃气弹射力Fg、法向接触力FN、切向接触力FT和气动力
Figure FDA0003590407010000074
的广义力；
在所述模块M4中：建立含非连续接触碰撞的多体动力学模型与气动力响应面模型耦合模拟方法，构建多体非连续接触动态分离安全性评估的时空数值仿真平台；
在所述模块M5中：评估内容包含分离过程中是否会发生卡滞、子弹与母弹外包络是否发生碰撞、子弹分离后的姿态角和姿态角速度是否满足控制系统要求。
10.一种子母弹分离仿真平台，其特征在于，采用权利要求1至6中任一项所述的子母弹分离安全性评价方法，或者权利要求7至9中任一项所述的子母弹分离安全性评价系统，对子母弹高速分离进行仿真，得到安全性的仿真结果并展示给用户。
CN202110950777.4A 2021-08-18 2021-08-18 子母弹分离安全性评价方法及系统及子母弹分离仿真平台 Active CN113761653B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110950777.4A CN113761653B (zh) 2021-08-18 2021-08-18 子母弹分离安全性评价方法及系统及子母弹分离仿真平台

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110950777.4A CN113761653B (zh) 2021-08-18 2021-08-18 子母弹分离安全性评价方法及系统及子母弹分离仿真平台

Publications (2)

Publication Number Publication Date
CN113761653A CN113761653A (zh) 2021-12-07
CN113761653B true CN113761653B (zh) 2022-06-07

Family

ID=78790361

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110950777.4A Active CN113761653B (zh) 2021-08-18 2021-08-18 子母弹分离安全性评价方法及系统及子母弹分离仿真平台

Country Status (1)

Country Link
CN (1) CN113761653B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111482533A (zh) * 2020-03-10 2020-08-04 西安近代化学研究所 一种二次起爆装置用分体式电缆分离装置
CN112964450A (zh) * 2021-02-07 2021-06-15 中国空气动力研究与发展中心超高速空气动力研究所 一种用于风洞多体分离试验模型轨迹预测方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003097453A2 (en) * 2002-05-21 2003-11-27 Nir Padan System and method for enhancing the payload capacity, carriage efficiency, and adaptive flexibility of external stores mounted on an aerial vehicle
CN111630953B (zh) * 2012-05-29 2015-03-11 湖北航天飞行器研究所 一种子弹超声速抛撒安全性分析方法
CN107991053B (zh) * 2017-11-06 2020-04-10 中国航天空气动力技术研究院 高速风洞投放模型试验飞行器分离轨迹预测方法及系统
CN112067248B (zh) * 2020-07-27 2022-09-23 中国航天空气动力技术研究院 一种两级运动的九自由度捕获轨迹试验装置及方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111482533A (zh) * 2020-03-10 2020-08-04 西安近代化学研究所 一种二次起爆装置用分体式电缆分离装置
CN112964450A (zh) * 2021-02-07 2021-06-15 中国空气动力研究与发展中心超高速空气动力研究所 一种用于风洞多体分离试验模型轨迹预测方法

Also Published As

Publication number Publication date
CN113761653A (zh) 2021-12-07

Similar Documents

Publication Publication Date Title
Panagiotopoulos et al. CFD transonic store separation trajectory predictions with comparison to wind tunnel investigations
Lijewski et al. Time-accurate computational fluid dynamics approach to transonic store separation trajectory prediction
Zhang et al. Applications of dynamic hybrid grid method for three-dimensional moving/deforming boundary problems
Xue et al. Derivation and validation of wind tunnel free-flight similarity law for store separation from aircraft
Hua et al. Effect of elastic deformation on the trajectory of aerial separation
Tian et al. A numerical method for multi-body separation with collisions
CN106570242A (zh) 低空大动压整体式整流罩高速分离流固耦合仿真方法
CN113761653B (zh) 子母弹分离安全性评价方法及系统及子母弹分离仿真平台
Rafeie et al. The Aerodynamic and Dynamic Analysis ofThree Common 4.5 mm Caliber Pellets in a Transonic Flow
Brock et al. Dynamic CFD simulations of the supersonic inflatable aerodynamic decelerator (SAID) ballistic range tests
Welterlen et al. Application of viscous, Cartesian CFD to aircraft stores carriage and separation simulation
Kanyshev et al. Methods of computing thrust vector coordinates for aircrafts equipped with thrust vector control
Murman et al. Cartesian-grid simulations of a canard-controlled missile with a spinning tail
Aliabadi et al. Parallel finite element computation of the dynamics of large ram air parachutes
Baayen Vortexje-an open-source panel method for co-simulation
CN102890734A (zh) 飞行模拟器中建立飞行结冰降阶模型的方法
Demir Computational Fluid Dynamics Modelling of Store Separation for Transonic Generic Store
Lasek et al. Modelling of external store separation
Coleman et al. Numerical simulation of a store separation event from an F-15E aircraft
Yang et al. Turbulent study of roll pitch yaw phenomena for a finned store released from a wing regarding different values of angle of attack
Lijewski et al. Comparison of transonic store separation trajectory predictions using the PEGASUS/DXEAGLE and BEGGAR codes
Potter et al. CFD Based Reduced Order Modeling for Hypersonic Vehicles Using CREATE (TM)-AV Kestrel
CN103902753A (zh) 飞行模拟器中建立飞行结冰降阶模型的方法
Ivanov et al. Numerical investigation of the EXPERT reentry vehicle aerothermodynamics along the descent trajectory
Sahu et al. Numerical computations of unsteady aerodynamics of maneuvering projectiles

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Wang Bolan

Inventor after: Li Weiwei

Inventor after: Liu Guang

Inventor after: Gong Yuzhe

Inventor after: Xu Bin

Inventor after: Hu Shan

Inventor after: Ge Zhishan

Inventor after: Yu Liujian

Inventor after: Liu Luguang

Inventor after: Zhang Baogang

Inventor before: Wang Bolan

Inventor before: Liu Guang

Inventor before: Xu Bin

Inventor before: Hu Shan

Inventor before: Ge Zhishan

Inventor before: Yu Liujian

Inventor before: Liu Luguang

Inventor before: Zhang Baogang

Inventor before: Li Weiwei