CN106055748B - 固体捆绑火箭捆绑位置的参数化判定方法 - Google Patents
固体捆绑火箭捆绑位置的参数化判定方法 Download PDFInfo
- Publication number
- CN106055748B CN106055748B CN201610343825.2A CN201610343825A CN106055748B CN 106055748 B CN106055748 B CN 106055748B CN 201610343825 A CN201610343825 A CN 201610343825A CN 106055748 B CN106055748 B CN 106055748B
- Authority
- CN
- China
- Prior art keywords
- rocket
- binding
- model
- finite element
- binding position
- 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.)
- Expired - Fee Related
Links
- 230000027455 binding Effects 0.000 title claims abstract description 68
- 238000009739 binding Methods 0.000 title claims abstract description 68
- 239000007787 solid Substances 0.000 title claims abstract description 20
- 238000000034 method Methods 0.000 title claims abstract description 19
- 230000008859 change Effects 0.000 claims abstract description 6
- 238000012805 post-processing Methods 0.000 claims description 8
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 238000003786 synthesis reaction Methods 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 5
- 239000000284 extract Substances 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000005096 rolling process Methods 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 14
- 230000008569 process Effects 0.000 abstract description 5
- 238000005457 optimization Methods 0.000 abstract description 2
- 238000010348 incorporation Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 6
- 201000002197 heart lymphoma Diseases 0.000 description 6
- 201000002113 hereditary lymphedema I Diseases 0.000 description 6
- 208000031223 plasma cell leukemia Diseases 0.000 description 6
- 238000004088 simulation Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 239000003380 propellant Substances 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Classifications
-
- 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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及固体捆绑火箭捆绑位置的参数化判定方法，属于航空航天领域。本发明的判定方法，首先采用梁‑集中质量模型对火箭模型进行简化，选取第一个捆绑组合位置为原始位置，利用有限元法对整箭进行有限元建模并解算；然后改变捆绑位置，并通过第一个捆绑组合位置的模型得到得到PCL语言描述的参数化形式模型，将新的捆绑位置替换原模型的捆绑位置后再次建模并解算，再采用相同方法解算剩余捆绑位置的模型。本发明的判断方法避免了固体捆绑火箭改变捆绑位置的有限元建模过程中给设计人员带来的重复劳动，提高了对火箭力学特性分析的设计效率，可以结合工程领域已有技术对固体捆绑火箭实现优化设计。
Description
技术领域
本发明涉及固体捆绑火箭捆绑位置的参数化判定方法，属于航空航天领域。
背景技术
运载火箭是将地面有效载荷送入空间轨道的有效技术手段，我国当前正在进行的探月工程、空间站工程和火星计划都与运载火箭技术息息相关。其中，火箭的力学特性分析是运载火箭设计阶段的关键部分。针对助推器的捆绑位置的不同，为了得到最佳的捆绑位置，必须对运载火箭进行参数化的建模与分析。
运载火箭力学特性包括静力学特性与动力学特性。质量特性在运载火箭力学特性分析中占有重要的地位，而推进剂的质量在火箭总质量的占绝大部分。针对运载火箭推进剂质量建模方式的不同，至今已形成三类建模分析方法：梁-集中质量模型、壳-集中质量模型和全箭实体模型。其中梁-集中质量模型在模型经济性和计算精度上有一定的优势。另外，由于改变运载火箭运载火箭助推器的捆绑位置和捆绑形式，每次都需要重新建立有限元模型，所以在初始设计阶段，捆绑位置的变化会带来有限元建模上的重复操作，严重影响设计效率。因此建立参数化的分析模型，使得重复性的建模过程简化并易于实施，是十分有必要的。
火箭结构示意图如图1所示，并按图2建立火箭坐标系，采用“I-II-III-IV”火箭坐标系统，即右手法则的笛卡尔坐标系，坐标系y轴正向与火箭坐标系III象限线重合，坐标系z轴正向与火箭坐标系II象限线重合，坐标系x轴正向遵守右手法则指向火箭箭轴负方向(即整流罩指向火箭尾段的方向)，坐标系原点位于整流罩von Karman曲线理论尖点。
针对固体捆绑火箭的捆绑形式，分析模型采用的三杆连接形式如图3所示，此捆绑位置离坐标系原点较近，为前捆绑。球铰连接形式如图4所示，由4个捆绑连杆和2个主捆绑接头组成，主捆绑结构球头位于距芯级150mm(水平)，距助推器150mm(水平)的位置，此捆绑位置离坐标系原点较远，为后捆绑。各取火箭助推器前后捆绑连接点5个位置，并对其进行排列组合，得到25个捆绑位置组合形式(如表2所示)，按表中的顺序对25个捆绑位置的固体捆绑火箭建立参数化的模型，进行仿真和结果的后处理与分析。其中，计算所需的箭体结构有限元模型主要由梁单元模拟，结构质量由集中质量单元模拟。捆绑连接部分：捆绑连杆采用杆单元模拟；主捆绑连接接头采用梁模型模拟，球头部分采用自由度释放(DOF Release)的方法释放梁单元的3个转动自由度，以实现球铰的功能。各发动机机架均采用梁模型模拟。
Patran和Nastran软件是国际公认的大型通用有限元分析软件。在众多的有限元分析软件中，Patran和Nastran软件提供的二次开发接口具有一定优势，其中脚本接口(Patran Scripting Interface)是在PCL语言的基础上进行的。本发明利用PCL语言驱动Patran和Nastran软件来进行固体捆绑火箭的建模、计算和后处理。
发明内容
本发明的目的是为了解决现有技术在固体捆绑运载火箭力学特性分析中的设计效率问题，提出一种针对固体捆绑运载火箭捆绑位置的参数化判定方法。
本发明的目的是通过以下技术方案实现的：
本发明的固体捆绑火箭捆绑位置的参数化判定方法，具体判定步骤如下：
步骤1：采用梁-集中质量模型对火箭模型进行简化，选取任意一个捆绑组合位置为原始位置，利用有限元法对整箭进行有限元建模，建模过程可以在MSC.Patran中实现。
步骤2：将步骤1中得到的有限元模型提交给有限元软件进行计算。计算能通过商业有限元软件Abaqus、MSC.Nastran、Ansys、HyperWorks中任意一款实现。
步骤3：根据步骤2得到的计算结果，从中提取站点截面的气动力、操纵力和晃动力，根据力的合成法则式1合成球头载荷与连杆轴力载荷；提取整箭的动特性的特征根，并按从小到大排列，按式2计算得到整箭的共振频率。将计算得到的数据利用MATLAB软件写入文本文件并画图。
其中，Ncombine是合成之后的力，Naero为气动力，Nmanoeuver为操纵力，Nslosh为晃动力；f为频率，w为圆频率，即得到的整箭的动特性特征根。
步骤4：改变捆绑位置，得到新的捆绑位置。利用步骤1得到的原始绑位有限元模型，得到PCL语言描述的参数化形式模型，该过程可以通过有限元软件MSC.Patran中实现。将新的捆绑位置替换原模型的捆绑位置，替换方式以文本读写方式进行，替换位置可以在MSC帮助文档中得到关键字sgm_const_curve_2d_circle_v2，改变关键字函数中的捆绑位置坐标即可完成PCL的二次开发。将改变捆绑位置后得到的参数化有限元模型提交给有限元软件进行计算，重复步骤2与步骤3，完成有限元计算和结果的后处理。
步骤5：重复该步骤4，直至得到所有捆绑位置的合成球头载荷与连杆轴力载荷，以及整箭的共振频率，以三连杆轴力合成载荷、球头合成载荷最小准则以及一阶模态频率最低准则对结果进行评价，从而选取最佳捆绑组合形式。
有益效果
本发明实现了固体捆绑火箭捆绑位置的参数化模型的自动化建立，避免了固体捆绑火箭改变捆绑位置的有限元建模过程中给设计人员带来的重复劳动。将建模与仿真和结果后处理模块集成为统一的参数化分析方法，解决了现有技术在固体捆绑改变捆绑位置中的火箭力学特性分析的设计效率问题。参数化分析的写入参数与执行程序分离，可以结合工程领域已有技术对固体捆绑火箭实现优化设计。
附图说明
图1为固体捆绑火箭结构示意图；
图2为火箭模型坐标系示意图；
图3为三连杆捆绑连接示意图；
图4为球头捆绑连接示意图；
1-有效载荷；2-整流罩；3-芯二级；4-助推器；5-芯一级。
具体实施方式
下面结合实施例和附图对本发明内容作进一步说明。
实施例
以判定某型号火箭的捆绑位置为例，该火箭的结构示意图如图1所示，按图2建立火箭坐标系，采用“I-II-III-IV”火箭坐标系统，即右手法则的笛卡尔坐标系，坐标系y轴正向与火箭坐标系III象限线重合，坐标系z轴正向与火箭坐标系II象限线重合，坐标系x轴正向遵守右手法则指向火箭箭轴负方向(即整流罩指向火箭尾段的方向)，坐标系原点位于整流罩von Karman曲线理论尖点。
针对固体捆绑火箭的捆绑形式，分析模型采用的三杆连接形式，如图3所示，此捆绑位置离坐标系原点较近，为前捆绑。球铰连接形式如图4所示，由4个捆绑连杆和2个主捆绑接头组成，主捆绑结构球头位于距芯级150mm(水平)，距助推器150mm(水平)的位置，此捆绑位置离坐标系原点较远，为后捆绑。各取火箭助推器前后捆绑连接点5个位置，并对其进行排列组合，得到25个捆绑位置组合形式(如表2所示)，按表中的顺序对25个捆绑位置的固体捆绑火箭建立参数化的模型，进行仿真和结果的后处理与分析。其中，计算所需的箭体结构有限元模型主要由梁单元模拟，结构质量由集中质量单元模拟。捆绑连接部分：捆绑连杆采用杆单元模拟；主捆绑连接接头采用梁模型模拟，球头部分采用自由度释放(DOFRelease)的方法释放梁单元的3个转动自由度，以实现球铰的功能。各发动机机架均采用梁模型模拟。
对上述固体捆绑火箭捆绑位置的参数化判定方法，具体判定步骤如下：
步骤1：采用梁-集中质量模型对火箭模型进行简化，以表2中的第一个捆绑组合位置为原始位置，利用有限元法对整箭进行有限元建模，建模过程可以在MSC.Patran中实现。
步骤2：将步骤1中得到的有限元模型提交给有限元软件进行计算。计算能通过商业有限元软件Abaqus、MSC.Nastran、Ansys、HyperWorks中任意一款实现。
步骤3：参数化结果后处理，根据步骤2得到的计算结果，从中提取站点截面的气动力、操纵力和晃动力，根据力的合成法则(式1)合成球头载荷与连杆轴力载荷；提取整箭的动特性的特征根，并按从小到大排列，按式2计算得到整箭的共振频率。将计算得到的数据利用MATLAB软件写入文本文件并画图。
其中，Ncombine是合成之后的力，Naero为气动力，Nmanoeuver为操纵力，Nslosh为晃动力；f为频率，w为圆频率，即得到的整箭的动特性特征根。
步骤4：按照表1顺序改变捆绑位置，得到新的捆绑位置。利用步骤1得到的原始绑位有限元模型，得到PCL语言描述的参数化形式模型，该过程可以通过有限元软件MSC.Patran中实现。将新的捆绑位置替换原模型的捆绑位置，替换方式以文本读写方式进行，替换位置可以在MSC帮助文档中得到关键字sgm_const_curve_2d_circle_v2，改变关键字函数中的捆绑位置坐标即可完成PCL的二次开发。将改变捆绑位置后得到的参数化模型提交给有限元软件进行计算，重复步骤2与步骤3，完成有限元计算和结果的后处理。
步骤5：重复该步骤4，得到25个捆绑位置的合成球头载荷与连杆轴力载荷，以及整箭的共振频率。通过对比分析，如表2所示，以三连杆轴力合成载荷、球头合成载荷最小准则以及一阶模态频率最低准则对结果进行了评价，选取了前捆绑点为2，后捆绑点为3时的捆绑组合形式为最佳捆绑形式。
表1为捆绑组合形式表格
表2为各捆绑组合计算结果
Claims (1)
1.固体捆绑火箭捆绑位置的参数化判定方法，其特征是具体判定步骤如下：
步骤1：采用梁-集中质量模型对火箭模型进行简化，选取任意一个捆绑组合位置为原始位置，利用有限元法对整箭进行有限元建模；
步骤2：将步骤1中得到的有限元模型提交给有限元软件进行计算；
步骤3：根据步骤2得到的计算结果，从中提取站点截面的气动力、操纵力和晃动力，根据力的合成法则式1合成球头载荷与连杆轴力载荷；提取整箭的动特性的特征根，并按从小到大排列，按式2计算得到整箭的共振频率；将计算得到的数据利用MATLAB软件写入文本文件并画图；
其中，Ncombine是合成之后的力，Naero为气动力，Nmanoeuver为操纵力，Nslosh为晃动力；f为频率，w为圆频率，即得到的整箭的动特性特征根；
步骤4：改变捆绑位置，得到新的捆绑位置；利用步骤1得到的原始绑位有限元模型，得到PCL语言描述的参数化形式模型，该过程可以通过有限元软件MSC.Patran中实现；将新的捆绑位置替换原模型的捆绑位置，替换方式以文本读写方式进行，替换位置可以在MSC帮助文档中得到关键字sgm_const_curve_2d_circle_v2，改变关键字函数中的捆绑位置坐标即可完成PCL的二次开发；将改变捆绑位置后得到的参数化有限元模型提交给有限元软件进行计算，重复步骤2与步骤3，完成有限元计算和结果的后处理；
步骤5：重复该步骤4，直至得到所有捆绑位置的合成球头载荷与连杆轴力载荷，以及整箭的共振频率，以三连杆轴力合成载荷、球头合成载荷最小准则以及一阶模态频率最低准则对结果进行评价，从而选取最佳捆绑组合形式。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610343825.2A CN106055748B (zh) | 2016-05-23 | 2016-05-23 | 固体捆绑火箭捆绑位置的参数化判定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610343825.2A CN106055748B (zh) | 2016-05-23 | 2016-05-23 | 固体捆绑火箭捆绑位置的参数化判定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106055748A CN106055748A (zh) | 2016-10-26 |
| CN106055748B true CN106055748B (zh) | 2019-02-05 |
Family
ID=57177485
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610343825.2A Expired - Fee Related CN106055748B (zh) | 2016-05-23 | 2016-05-23 | 固体捆绑火箭捆绑位置的参数化判定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106055748B (zh) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2001064509A2 (en) * | 2000-02-29 | 2001-09-07 | Lockheed Martin Corporation | Aerodynamic fairing |
| CN103455646A (zh) * | 2012-05-31 | 2013-12-18 | 北京宇航系统工程研究所 | 一种主捆绑连接结构动力学模拟方法 |
-
2016
- 2016-05-23 CN CN201610343825.2A patent/CN106055748B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2001064509A2 (en) * | 2000-02-29 | 2001-09-07 | Lockheed Martin Corporation | Aerodynamic fairing |
| CN103455646A (zh) * | 2012-05-31 | 2013-12-18 | 北京宇航系统工程研究所 | 一种主捆绑连接结构动力学模拟方法 |
Non-Patent Citations (1)
| Title |
|---|
| 捆绑火箭助推器与芯级间捆绑联接技术应用进展;冯韶伟 等;《导弹与航天运载技术》;20120630(第6期);第20-23页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106055748A (zh) | 2016-10-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109325271B (zh) | 一种运载火箭结构有限元模型参数化建模方法 | |
| CN107491616A (zh) | 一种适用于格栅构型舵面的结构有限元参数化建模方法 | |
| CN103455645B (zh) | 一种全箭模态提取方法 | |
| CN110889166B (zh) | 基于热力耦合约束的航空用轴承支架轻量化设计方法 | |
| CN105005642A (zh) | 一种三维模型批量格式转换及轻量化方法 | |
| CN102201017B (zh) | 汽车车身的优化设计方法 | |
| CN112364288B (zh) | 一种反应堆多物理场耦合计算系统及方法 | |
| CN104123400A (zh) | 基于力法的Global-Local细节有限元分析方法 | |
| Krist et al. | Kestrel Results at Liftoff Conditions for a Space Launch System Configuration in Proximity to the Launch Tower | |
| CN111898295A (zh) | 一种变刚度复合材料层合板的有限元建模方法 | |
| CN106055748B (zh) | 固体捆绑火箭捆绑位置的参数化判定方法 | |
| CN102830978A (zh) | 一种可视化3d编程系统及应用 | |
| Ferreres et al. | A skew mu toolbox (SMT) for robustness analysis | |
| CN107742042A (zh) | 一种cae仿真前处理系统及方法 | |
| CN109002597A (zh) | 一种基于数字化双胞胎技术的虚实结合运维仿真调试平台 | |
| CN110874517A (zh) | 一种快速验证fpga互连线连接正确的方法 | |
| CN106777640A (zh) | 一种基于pdms的模型阵列复制的方法 | |
| CN104133923A (zh) | 一种车身快速改型设计方法 | |
| CN106407600B (zh) | 一种基于流-固耦合及代理模型的入水仿真分析方法 | |
| CN108921950A (zh) | 三维裂缝模拟的方法以及相关装置 | |
| CN112182873A (zh) | 一种基于硬件系统对航空发动机模型加速的方法 | |
| CN105786786A (zh) | 一种web报表数据异步填充的实现方法 | |
| Péron et al. | A mixed overset grid/immersed boundary approach for CFD simulations of complex geometries | |
| CN110992475A (zh) | 水利水电工程大体积混凝土快速计算任意部位工程量的方法、系统及介质 | |
| CN105783931A (zh) | 一种电子地图编译方法及装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190205 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |