CN114065672A - 一种液体运载火箭增压输送系统管路变形补偿量计算方法 - Google Patents

一种液体运载火箭增压输送系统管路变形补偿量计算方法 Download PDF

Info

Publication number
CN114065672A
CN114065672A CN202111265640.1A CN202111265640A CN114065672A CN 114065672 A CN114065672 A CN 114065672A CN 202111265640 A CN202111265640 A CN 202111265640A CN 114065672 A CN114065672 A CN 114065672A
Authority
CN
China
Prior art keywords
pipeline
working condition
deformation
under
axial
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.)
Granted
Application number
CN202111265640.1A
Other languages
English (en)
Other versions
CN114065672B (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.)
CASIC Rocket Technology Co
Original Assignee
CASIC Rocket Technology Co
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 CASIC Rocket Technology Co filed Critical CASIC Rocket Technology Co
Priority to CN202111265640.1A priority Critical patent/CN114065672B/zh
Priority claimed from CN202111265640.1A external-priority patent/CN114065672B/zh
Publication of CN114065672A publication Critical patent/CN114065672A/zh
Application granted granted Critical
Publication of CN114065672B publication Critical patent/CN114065672B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明涉及一种液体运载火箭增压输送系统管路变形补偿量的计算方法,包括如下步骤:按照火箭从总装测试到最终发射的先后顺序,将火箭的工作工况分为气密增压工况、推进剂加注工况、射前增压工况和飞行工况;梳理每种工况下影响管路变形的载荷类型;分析每种工况下各种载荷与管路变形量的定量关系,得出计算公式,根据公式计算得出各载荷作用下的管路端点位移及管路热变形量,合并计算得出每种工况下管路端点位移总量及管路热变形量总量;按照定量关系,将每种工况下管路端点位移总量和管路热变形量总量折算成管路要求的补偿量;汇总各工况下的波纹管的轴向和横向补偿量,得出最终要求的波纹管补偿量范围,为波纹管设计提供参考和依据。

Description

一种液体运载火箭增压输送系统管路变形补偿量计算方法
技术领域
本发明属于运载火箭增压输送系统领域,具体涉及一种液体运载火箭增压输送系统管路变形补偿量的计算方法。
背景技术
液体火箭增压输送系统管路在火箭加注、飞行过程中需要承受温度、内压以及其他外部载荷的作用,在这些载荷的耦合作用下,管路会发生热胀冷缩变形以及其他相对位移,为了保证管路的稳定和安全,减少由于管变形产生的内应力,针对直径较大的管路,常常采用波纹管补偿器进行变形补偿(吸收)。内压载荷和位移载荷是波纹管补偿器设计时重要输入参数指标。内压一般由动力系统直接给出,位移载荷即波纹管需要吸收的管路变形补偿量,由管路设计根据火箭在不同的工况下的载荷情况计算得出。大直径管路一端往往与贮箱(或发动机)连接,由于火箭工况比较复杂,贮箱在温度、内压、轴压等不同载荷作用下发生变形,导致与之相连的管路端点发生位移,这部分的影响在管路的补偿设计中同样不可忽视,现有技术鲜有全面考虑各种工况的影响,在进行管路的补偿计算时更多的侧重于热补偿的计算,即管路在高温或低温下的热胀冷缩变形量计算,依据热变形量进行的波纹管的设计往往过于保守、粗放。因此,如何总结出一个全面系统的管路变形补偿量计算方法对于波纹管的精准、合理设计至关重要。
发明内容
本发明提出了一种液体运载火箭增压输送系统管路变形补偿量的计算方法,该计算方法综合考虑了火箭从测试到发射过程中的各个工况下管路的变形,为管路波纹管补偿器的设计提供参考和依据。
本发明提出一种液体运载火箭增压输送系统管路变形补偿量的计算方法,具体技术方案包括如下步骤:
S1、按照火箭从总装测试到最终发射的先后顺序,将火箭的工作工况分为气密增压工况、推进剂加注工况、射前增压工况和飞行工况;
S2、系统梳理每种工况下影响管路变形的主要载荷类型;气密增压工况载荷类型主要为内压,推进剂加注工况载荷类型主要为温度,射前增压工况主要载荷类型为温度和内压,飞行工况载荷类型主要为温度、内压、轴压和发动机推力;
S3、分析每种工况下各种载荷与管路变形量的定量关系,得出计算公式,根据公式计算得出各载荷作用下的管路端点位移及管路热变形量,合并计算得出每种工况下管路端点位移总量及管路热变形量总量;
S4、分析管路变形量与波纹管补偿量之间的定量关系,按照定量关系,将每种工况下管路端点位移总量和管路热变形量总量,折算成管路要求的补偿量,即得到每种中工况下管路中波纹管的轴向和横向补偿量;
S5、汇总各工况下的波纹管的轴向和横向补偿量,得出最终要求的波纹管补偿量范围,为波纹管设计提供参考和依据。
进一步地,步骤S3中,内压载荷引起的贮箱变形量计算公式:
Figure BDA0003326890600000021
Figure BDA0003326890600000022
Figure BDA0003326890600000023
式中:
δ1—贮箱柱段轴向伸长量,mm;
δ2—柱段半径伸长量mm;
P—内压力,MPa;
H—贮箱柱段长度,mm;
E—贮箱材料弹性模量,MPa;
d—贮箱壁厚,mm;
r—贮箱半径,mm;
μ—贮箱材料泊松比。
进一步地,步骤S3中,温度载荷引起的贮箱或管路变形量计算公式:
低温推进剂加注后,假设贮箱柱段和前后底均由初始温度T0变化到T1,任意两点距离为L,则低温变形量为:
ΔL=αL(T1-T0)
式中:
ΔL-变形量,mm;
L-任意两点间距离,mm。
α-材料的热膨胀系数,mm/(mm·K);
T0-初始温度,取室温,K;
T1-介质温度,取推进剂的温度,K;
进一步地,步骤S3中,轴压载荷F作用下贮箱柱段轴向变形量δ3、径向变形量δ4计算公式如下
Figure BDA0003326890600000031
Figure BDA0003326890600000032
式中:
δ3-贮箱柱段轴向变形量,mm;
δ4-贮箱柱段径向变形量,mm;
F-轴压力,N;
H-贮箱柱段长度,mm;
E-贮箱材料弹性模量,MPa;
d-贮箱壁厚,mm;
r-贮箱半径,mm;
μ-贮箱材料泊松比。
进一步地,步骤S3中,发动机推力引起的变形量发动机专业直接提供。
进一步地,所述步骤S3中,合并计算具体为:将各个工况下的各载荷作用产生的管路端点轴向位移和径向位移分别相加;将各个工况下的各载荷作用产生管路轴向热变形量和径向热变形量分别相加
本发明与现有技术相比的有益效果是:
本发明提出了一种液体运载火箭增压输送系统管路变形补偿量的计算方法,该计算方法全面考虑了火箭从测试到发射过程中的各个工况下不同载荷对管路变形的影响,分别计算每个载荷作用下管路的变形量,最后综合计算得出管路总的变形补偿量,采用该方法计算得出数据能够更加准确的用于指导波纹管的精准、合理设计,且具有一定的通用性。
附图说明
附图1为管路变形补偿量计算方法流程图;
附图2为低温燃料贮箱预冷回流管布局图。
附图说明:1-端点1,2-端点2,3-端点3,4-管段1,5-管段2,6-波纹管1,7-波纹管2,8-波纹管3,9-低温燃料贮箱。
具体实施方式
为了使本发明的结构的功能更加清晰,下面结合附图和具体实施方式对本发明的技术方案作进一步说明。
以低温燃料贮箱上的预冷回流管为例,结合附图2对本发明提出的管路变形量计算方法进行详细说明。附图2中,端点1为贮箱柱段法兰与管路连接位置,端点2为设置在贮箱后端壳上的止推支架,为固定点(位移量为0),端点3为管路与发动机上的法兰连接位置;管段1为端点1与端点2之间的直管段,管段2为端点2与端点3之间的Z形管段;管段1上设置一处不带钢丝网套的波纹管1,管段2上设置两个钢丝网套波纹管。预冷回流管从发动机出口(端点3)引出,从舱壁穿出后沿贮箱侧壁分布,引至贮箱柱段与贮箱上的法兰口(端点1)相连,其中H为柱段总长,H1是实际有效的柱段长度。
一种液体运载火箭增压输送系统管路变形补偿量的计算方法,具体步骤包括:
S1、按照火箭从总装测试到最终发射的先后顺序,将火箭的工作工况分为气密增压工况、推进剂加注工况、射前增压工况和飞行工况。
S2、梳理每种工况下影响管路变形的主要载荷类型,分别为气密增压工况载荷类型为内压,推进剂加注工况载荷类型为温度,射前增压工况载荷类型为温度和内压,飞行工况载荷类型为温度、内压、轴压和发动机推力。
S3、计算得出每种工况下各载荷作用下的管路端点位移及管路热变形量,合并计算得出每种工况下管路端点位移总量及管路热变形量总量。
(1)气密增压工况:载荷类型为内压,贮箱在内压下膨胀,柱段发生轴向和径向伸长,端点1产生轴向和径向位移,计算得出:
端点1的轴向位移量为:
Figure BDA0003326890600000051
端点1的径向位移量为:
Figure BDA0003326890600000052
(2)推进剂加注工况:载荷类型为温度,贮箱在低温下收缩,柱段H1轴向和径向缩短导致端点1产生轴向和径向位移,计算得出:
端点1的轴向位移量为:
ΔL1=αH1(T1-T0)
端点1的径向位移量为:
ΔL2=αr(T1-T0)
管段1轴向收缩,管段2的收缩分解为轴向和径向两个方向的分量。计算得出:
管段1的轴向变形量为:
ΔL1x=αL1(T1-T0)
管段2的轴向变形量为:
ΔL2x=αL2x(T1-T0)
管段2的径向变形量为:
ΔL2y=αL2y(T1-T0)
(3)射前增压工况:载荷类型为内压+温度,各端点位移及管路变形量为加注工况与气密增压工况的叠加;
端点1的轴向位移量为:δ1+ΔL1
端点1的径向位移量为:δ2+ΔL2
管段1的轴向变形量为:ΔL1x
管段2的轴向变形量为:ΔL2x
管段2的径向变形量为:ΔL2y
(4)飞行工况:载荷类型为温度+内压+轴压+发动机推力,在射前增压基础上,考虑在轴压及发动机推力作用下各端点位移量。各变形量计算时在射前增压工况基础上进行叠加。
在轴压作用下,贮箱柱段发生轴向缩短和径向伸长,导致端点1产生轴向和径向位移,计算得出:
端点1的轴向位移量为:
Figure BDA0003326890600000061
端点1的径向位移量为;
Figure BDA0003326890600000062
发动机推力引起发动机机架轴向变形量Δy,径向变形量Δx,由发动机专业直接提供,导致端点3产生轴向和径向位移:
端点3的轴向位移量为:Δy
端点3的径向位移量为:Δx
合并计算飞行状况下各载荷导致的端点位移量和管路变形量为:
端点1的轴向位移量为:δ1+ΔL13
端点1的径向位移量为:δ2+ΔL24
端点3的轴向位移量为:Δy
端点3的径向位移量为:Δx
管段1的轴向变形量为:ΔL1x
管段2的轴向变形量为:ΔL2x
管段2的径向变形量为:ΔL2y
上述式中各物理量含义为:
δ1-贮箱柱段轴向伸长量,mm;
δ2-柱段半径伸长量mm;
P-内压力,MPa;
H1-实际有效贮箱柱段长度,mm;
E-贮箱材料弹性模量,MPa;
d-贮箱壁厚,mm;
r-贮箱半径,mm;
μ-贮箱材料泊松比。
ΔL-变形量,mm;
L-任意两点间距离,mm。
α-材料的热膨胀系数,mm/(mm·K);
T0-初始温度,取室温,K;
T1-介质温度,取推进剂的温度,K;
δ3-贮箱柱段轴向变形量,mm;
δ4-贮箱柱段径向变形量,mm;
F-轴压力,N;
将每种工况下位移总量及管路热变形量总量进行汇总如表1和表2所示。
表1
Figure BDA0003326890600000071
表2
Figure BDA0003326890600000072
S4、分析管路变形量(端点位移和自身热变形)与波纹管补偿量之间的定量关系,得到结论如下:
A.管路端点位移与波纹管补偿量之间的关系:与贮箱相连的管路端点的轴向或径向位移会直接引起与之相关的管路中的波纹管发生轴向或横向变形。如图2所示,端点1轴向(或径向)位移会引起波纹管1的轴向(或横向)变形;端点3的轴向(或径向)位移会引起波纹管2(或波纹管3)的横向变形等(此处假设端点3的轴向位移全部被波纹管2吸收,径向位移全部被波纹管3吸收)。
B.管路热变形量与波纹管量之间的关系:管路自身的伸长或缩短会直接引起与之相连的波纹管发生轴向或横向变形。如图2所示,管段1的轴向缩短(或伸长)会引起波纹管1的轴向伸长(或压缩)变形,管段3的轴向缩短(或伸长)会引起波纹管2的横向变形,管段3的径向缩短(或伸长)会引起波纹管3的横向变形(此处假设管段3的轴向变形全部被波纹管2横向吸收,管段3的径向变形全部被波纹管3横向吸收)。
按照以上定量关系,将每种工况下管路端点位移总量和管路热变形量总量,折算成管路要求的补偿量,即得到每种工况下管路波纹管的轴向和横向补偿量,如表3所示。
表3
Figure BDA0003326890600000081
S5、对比分析表3中每一个波纹管在不同工况下的轴向和横向补偿量,得出最终要求的波纹管补偿量范围,为波纹管设计提供参考和依据。

Claims (6)

1.一种液体运载火箭增压输送系统管路变形补偿量的计算方法,其特征在于,包括如下步骤:
S1、按照火箭从总装测试到最终发射的先后顺序,将火箭的工作工况分为气密增压工况、推进剂加注工况、射前增压工况和飞行工况;
S2、系统梳理每种工况下影响管路变形的主要载荷类型;
气密增压工况载荷类型主要为内压,推进剂加注工况载荷类型主要为温度,射前增压工况主要载荷类型为温度和内压,飞行工况载荷类型主要为温度、内压、轴压和发动机推力;
S3、分析每种工况下各种载荷与管路变形量的定量关系,得出计算公式,根据公式计算得出各载荷作用下的管路端点位移及管路热变形量,合并计算得出每种工况下管路端点位移总量及管路热变形量总量;
S4、分析管路变形量与波纹管补偿量之间的定量关系,按照定量关系,将每种工况下管路端点位移总量和管路热变形量总量,折算成管路要求的补偿量,即得到每种中工况下管路中波纹管的轴向和横向补偿量;
S5、汇总各工况下的波纹管的轴向和横向补偿量,得出最终要求的波纹管补偿量范围,为波纹管设计提供参考和依据。
2.根据权利要求1所述的一种液体运载火箭增压输送系统管路变形补偿量的计算方法,其特征在于,步骤S3中,内压载荷引起的贮箱发生膨胀变形量计算公式:
Figure FDA0003326890590000011
Figure FDA0003326890590000012
式中:
δ1—贮箱柱段轴向伸长量,mm;
δ2—柱段半径伸长量mm;
P—内压力,MPa;
H—贮箱柱段长度,mm;
E—贮箱材料弹性模量,MPa;
d—贮箱壁厚,mm;
r—贮箱半径,mm;
μ—贮箱材料泊松比。
3.根据权利要求1所述的一种液体运载火箭增压输送系统管路变形补偿量的计算方法,其特征在于,步骤S3中,温度载荷引起的贮箱或管路变形量计算公式:
低温推进剂加注后,假设贮箱柱段和前后底均由初始温度T0变化到T1,任意两点距离为L,则低温变形量为:
ΔL=αL(T1-T0)
式中:
ΔL—变形量,mm;
L—任意两点间距离,mm。
α—材料的热膨胀系数,mm/(mm·K);
T0—初始温度,取室温,K;
T1—介质温度,取推进剂的温度,K。
4.根据权利要求1所述的一种液体运载火箭增压输送系统管路变形补偿量的计算方法,其特征在于,步骤S3中,轴压载荷F作用下贮箱柱段轴向变形量δ3、径向变形量δ4计算公式如下:
Figure FDA0003326890590000021
Figure FDA0003326890590000022
式中:
δ3—贮箱柱段轴向变形量,mm;
δ4—贮箱柱段径向变形量,mm;
F—轴压力,N;
H—贮箱柱段长度,mm;
E—贮箱材料弹性模量,MPa;
d—贮箱壁厚,mm;
r—贮箱半径,mm;
μ—贮箱材料泊松比。
5.根据权利要求1所述的一种液体运载火箭增压输送系统管路变形补偿量的计算方法,其特征在于,步骤S3中,发动机推力引起的变形量发动机专业直接提供。
6.根据权利要求1所述的一种液体运载火箭增压输送系统管路变形补偿量的计算方法,其特征在于,所述步骤S3中,合并计算具体为:将各个工况下的各载荷作用产生的管路端点轴向位移和径向位移分别相加;将各个工况下的各载荷作用产生管路轴向热变形量和径向热变形量分别相加。
CN202111265640.1A 2021-10-28 一种液体运载火箭增压输送系统管路变形补偿量计算方法 Active CN114065672B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111265640.1A CN114065672B (zh) 2021-10-28 一种液体运载火箭增压输送系统管路变形补偿量计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111265640.1A CN114065672B (zh) 2021-10-28 一种液体运载火箭增压输送系统管路变形补偿量计算方法

Publications (2)

Publication Number Publication Date
CN114065672A true CN114065672A (zh) 2022-02-18
CN114065672B CN114065672B (zh) 2024-06-25

Family

ID=

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114595618A (zh) * 2022-05-10 2022-06-07 中建安装集团有限公司 一种竖井管道固定支架荷载计算系统、方法及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112128498A (zh) * 2020-09-02 2020-12-25 航天科工火箭技术有限公司 一种火箭侧壁管路的补偿装置及方法
WO2021097738A1 (zh) * 2019-11-21 2021-05-27 大连理工大学 基于改进型Smith预估器的航空发动机H∞控制方法
CN112989513A (zh) * 2021-03-04 2021-06-18 河南科技大学 一种试验与数值计算相结合获取轴承工作游隙的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021097738A1 (zh) * 2019-11-21 2021-05-27 大连理工大学 基于改进型Smith预估器的航空发动机H∞控制方法
CN112128498A (zh) * 2020-09-02 2020-12-25 航天科工火箭技术有限公司 一种火箭侧壁管路的补偿装置及方法
CN112989513A (zh) * 2021-03-04 2021-06-18 河南科技大学 一种试验与数值计算相结合获取轴承工作游隙的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
华寅淞;丁关锦;张游;堵同亮;铁文军;: "运载火箭贮箱排气管路中波纹管力学环境地面模拟验证分析", 飞控与探测, no. 01, 25 January 2019 (2019-01-25), pages 71 - 76 *
吴兵;杨蓉;周江帆;徐国梁;: "CZ-5大型结构及低温贮箱静强度试验技术研究", 强度与环境, no. 01, 15 February 2017 (2017-02-15), pages 1 - 10 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114595618A (zh) * 2022-05-10 2022-06-07 中建安装集团有限公司 一种竖井管道固定支架荷载计算系统、方法及存储介质
CN114595618B (zh) * 2022-05-10 2022-07-22 中建安装集团有限公司 一种竖井管道固定支架荷载计算系统、方法及存储介质

Similar Documents

Publication Publication Date Title
CN114065672A (zh) 一种液体运载火箭增压输送系统管路变形补偿量计算方法
CN114065672B (zh) 一种液体运载火箭增压输送系统管路变形补偿量计算方法
Sorour et al. Limit load analysis of thin-walled as-fabricated pipe bends with low ovality under in-plane moment loading and internal pressure
US9970357B2 (en) Superelastic fluid conduit for a gas turbine engine
CN101907212B (zh) 一种外压内置万能平衡型膨胀节
CN104976278A (zh) 一种液压式扭矩补偿装置
Lak et al. Advancement in cryogenic propulsion system performance through propellant densification
CN213576156U (zh) 一种火箭侧壁管路的补偿装置
CN112128498B (zh) 一种火箭侧壁管路的补偿装置及方法
CN103672280B (zh) 长距离宽温区导管补偿装置
Nagai et al. Status of H-II rocket first stage propulsion system
CN112128498A (zh) 一种火箭侧壁管路的补偿装置及方法
CN114483377B (zh) 一种倾斜导杆式推力架泵前管路力平衡补偿结构
CN110489825A (zh) 一种大通径类空气管路的补偿设计方法
CN212778046U (zh) 一种下进下出式过冷器
CN103678879A (zh) 一种输送管载荷分析方法
CN203703479U (zh) 一种压力平衡型管路补偿器
Yoo et al. Static and dynamic structural analyses for a 750 kN class liquid rocket engine with TVC actuation
Korte et al. Multidisciplinary approach to linear aerospike nozzle optimization
Praveen et al. Thermo-structural Design of an AA2014 mounting between a Structure and Liquid Oxygen Propellant Tank of a Launch Vehicle
RU2399563C2 (ru) Двигательная установка ракетного блока
Kurosu et al. LE-X-Japanese next liquid booster engine
US3237975A (en) Balanced fluid coupling joint
Swanson et al. Design analysis of the Ares I Pogo accumulator
CN210800399U (zh) 一种万向抗震波纹管膨胀节

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