CN113434959B - 一种三维空间可行域内的管路弯折设计方法 - Google Patents

一种三维空间可行域内的管路弯折设计方法 Download PDF

Info

Publication number
CN113434959B
CN113434959B CN202110717390.4A CN202110717390A CN113434959B CN 113434959 B CN113434959 B CN 113434959B CN 202110717390 A CN202110717390 A CN 202110717390A CN 113434959 B CN113434959 B CN 113434959B
Authority
CN
China
Prior art keywords
pipeline
layout
bending
feasible region
section
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
CN202110717390.4A
Other languages
English (en)
Other versions
CN113434959A (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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN202110717390.4A priority Critical patent/CN113434959B/zh
Publication of CN113434959A publication Critical patent/CN113434959A/zh
Application granted granted Critical
Publication of CN113434959B publication Critical patent/CN113434959B/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/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft 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
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/14Pipes
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Pipeline Systems (AREA)

Abstract

本发明一种三维空间可行域内的管路弯折设计方法,属于航空航天飞行器的管路设计领域;首先管路布局的三维空间可行域,然后在布局可行域内确定管路弯折数量k,最后以管路的动力学指标为目标实施优化设计。本发明对于飞行器管路系统的任意两个出入口端口之间的管路布局设计提供了方法,发明了一段式~多段式弯折方案。能够对不同管路布局进行分析和选择,便于工程中多种布局管路结构的快速建模和频率特性对比。本发明以避开发动机风扇转子转动频率、泵设备运行频率、降低振动应力等动力学物理量为目标,能够达到错开飞行器的危险共振频率、降低管路结构的振动应力的目的,从而提高飞行器管路结构的振动可靠性。

Description

一种三维空间可行域内的管路弯折设计方法
技术领域
本发明属于航空航天飞行器的管路设计领域,具体涉及一种三维空间可行域内的管路弯折设计方法,即在已知管路两端接口(即入口、出口)位置的情况下,在两端 口之间的三维空间可行域内进行管形布局设计的方法。
背景技术
航空航天各类飞行器中,管路系统经常用来连接飞行器上两个或者多个附件设备, 是各种流体介质(液压油、燃油、滑油、空气等)相互传输的通道。管路的敷设经常 是在设计好的飞行器结构上,根据附件设备接口位置来实施管路布局设计。
各种管路结构一直是航空航天飞行器结构失效的薄弱环节。一方面,航空航天管路系统经常处于恶劣的振动冲击环境中,发动机旋转部件、泵设备运行、燃气喷射等 引起的振动激励均会传递给管路结构。另一方面,管路的动力学特性与管形布局密切 相关,在较宽的发动机工作频段内,管路极易引起共振或者应力过大,严重影响飞行 器的寿命及其安全性,通过合理的管形设计是提高管路振动可靠性的必由之路。
目前,飞行器管路的动力学全流程正向设计研究还处于探索阶段,其中基本管形(例如L形、U形或Z形)布局设计较为成熟。但是对于更加复杂的管路,其布局设 计主要借助于商用软件(例如UG、CATIA等)建模,再校核管路动响应的方法。该方 法存在以下三个问题:(1)管路设计手册标准仅是一般工程规则,具体到某些管路布 局设计时,仍然有多个布局参数(如弯折尺寸和位置)需要设计者来确定,此时就会 引入较大的经验性。(2)实际工程中,飞行器中各种管路的数量多达上千根,各个管 路之间形式差异很大,传统方法没有体现出不同管路的内在结构几何特征与动力学特 性的关系。若要逐一研究每根不规则管路的动态特性,分析工作量巨大。(3)在传统 的管路故障溯源和改进设计中,难以确定布局尺寸对管路结构的动态影响规律,无法 快速地甄别关键结构参数,排故中经常采用试凑调整的方法,盲目性较大,且效率低 下。而造成上述问题的根本原因是,在三维空间内设计未知的管形布局结构缺乏方法。 显然,必须将管路布局参数与其动力学影响联系起来,才能真正打通管路的动力学正 向设计。因此,研究基于动力学的飞行器管路三维空间布局设计方法具有重要的意义。
目前管路系统三维布局优化设计方法多以最短路径、路径寻优、避障和防止干涉为目标,主要以改善优化算法的效率、搜索策略和收敛性等为突破点[1-5],包括遗传 算法、NSGA-Ⅱ、粒子群算法等,这类研究代表性的文献有:
[1]熊勇,张加,余嘉俊,张本任,梁萱卓,朱奇舸.船舶三维管路智能布局优化算法[J].计算机应用,2020,40(07):2164-2170.
[2]柳强,焦国帅.基于改进NSGA-Ⅱ的航空发动机管路多目标布局优化[J].计算机集成制造系统,2018,24(05):1217-1227.
[3]王长涛,王哲,高治军,孙亮亮,朱毅.一种三维建筑单管路自动布局方法设计[J].沈阳建筑大学学报(自然科学版),2018,34(01):115-121.
[4]吴宏超.复杂产品中的管路自动布局与优化技术研究[D].北京理工大 学,2015.
[5]白晓兰.航空发动机管路布局智能设计方法研究[D].东北大学,2013.
由于飞行器、发动机上的管路系统处于恶劣的振动环境,上述管路布局设计方法未考虑管路系统在振动环境下的模态特性和共振失效问题。文献[6]虽然以控制动力学 指标为目标,研究了发动机弧面管路的布局变量优化设计方法,但是其管形布局形式 是给定的,即已有管形布局尺寸优化。但是工程中还存在一类管路未知的情况,即需 要在一个三维空间中进行合理的管路布局设计。
[6]刘伟,朱宏彦,赵玉杰,岳珠峰.发动机机匣弧面不规则管路的布局参数多目标优化设计[J].航空动力学报,2021,36(01):148-156.
专利号是CN110654333A的一种专用车液压管路布置方法,公开的管路布局方式是将液压钢管由传统的车体两侧布置改变为在副车架内部布置,避免了较多的液压管路 裸露在车体两侧容易与环境发生干涉而破损的问题。一体化管路布置立式大型发酵罐 及其管路布置方法,专利号CN108004120A,公开了立式罐体管路的布置方式,将立式 罐体具有的各种接口连接的相应各路工艺管路固定在立式管架上,各路工艺管路上安 装有相应控制阀门,主要目的是确保夹套(列管)内的冷却水排尽,以提高立式大型 发酵罐的生产效率。
上述专利均未考虑管路系统在振动环境下的模态特性和共振失效问题,目前在关于管路动力学布局方面的专利处于空白。
发明内容
要解决的技术问题:
为了避免现有技术的不足之处,本发明提出一种三维空间可行域内的管路弯折设计方法,通过参数化的布局参数,实现的布局设计方法。在一个三维空间内,根据管 路的两端接口位置、接口引入方向直线段要求、管径、弯曲半径等要求,提出了布局 可行域的概念,实现三维空间可行域内管路布局参数化建模,并以管路的动力学指标 为目标实施优化设计。本发明以避开发动机风扇转子转动频率、泵设备运行频率、降 低振动应力等动力学物理量为目标,能够达到错开飞行器的危险共振频率、降低管路 结构的振动应力的目的,从而提高飞行器管路结构的振动可靠性。
本发明的技术方案是:一种三维空间可行域内的管路弯折设计方法,其特征在于具体步骤如下:
步骤一:确定布局可行域;
首先,通过确定管路所要连接的入口A0、出口B0的坐标位置和入口方向nA,出口 方向nB;然后,依据航空航天管路设计标准,确定入口A0、出口B0处需要预留的接口 过渡段,接口过渡段的总长度为Ls+Rb,其中,Ls为管路接口过渡段的直线段,Rb为 管路接口过渡段的圆弧段的弯曲半径;
最后,确定出管路布局的三维空间可行域,其中,布局可行域是以A1B1为对角线 的长方体空间,A1和B1的坐标为:
Figure BDA0003135406250000031
管路布局可行域的长、宽、高分别为:
Figure BDA0003135406250000032
步骤二:在布局可行域内确定管路弯折数量k,具体根据A1B1的长度进行确定:
Figure BDA0003135406250000041
Figure BDA0003135406250000042
其中,k=0表示可行域空间不足,无法在可行域中布局;当2Rb≤A1B1<4Rb,采 用一段弯折,管路为直接连接布局可行域的对角线;当4Rb≤A1B1<6Rb,采用一段、两 段弯折管形,管路在布局可行域内确定一个弯折点;当A1B1≥6Rb,采用一段、两段、 三段弯折管形,管路在布局在可行域内确定两个弯折点;
步骤三:首先,将步骤二确定的弯折点坐标参数化,并赋予初始值,建立上述不 同布局的管路模型后,进一步计算分析管路结构的模态特性;然后,根据实际飞行器 的振动频率环境,管路通过对弯折点的优化设计来错开飞行器的危险共振频率;最后, 确保管路能够获得符合动力学特性的管形布局设计。
本发明的进一步技术方案是:所述步骤一中,为了确定布局可行域的界限,需要确定A1和B1点的坐标,并考虑该接口过渡段的方向,其中,入口A0的直线引入段方向 为nA,出口B0直线引出段的方向为nB,该方向用归一化的列矢量形式表示为:
Figure BDA0003135406250000043
式中,i,j,k分别表示直角坐标系x,y,z三个轴正向单位矢,(aA,bA,cA)、(aB,bB,cB)分别表示入口和出口直线段的归一化方向,即
Figure BDA0003135406250000044
则A1和B1点的坐标采用向量表示如下:
rA1=rA0+(Ls+Rb)nA
rB1=rB0-(Ls+Rb)nB
得到坐标表达式为:
Figure BDA0003135406250000045
本发明的进一步技术方案是:所述步骤二中,布局可行域内的一段式弯折布局方式为,将布局可行域的对角线A1 B1采用直线连接。
本发明的进一步技术方案是:所述步骤二中,布局可行域内的两段式弯折布局方式:首先确定折弯点C1坐标为(xC1,yC1,zC1),然后将A1、B1分别与C1点直线连接,即 得到两段式弯曲的管路布局形式。
本发明的进一步技术方案是:所述步骤二中,折弯点C1的位置确定方法:采用六种特殊位置情况,对一阶固有频率做评估;所述六种特殊位置是指折弯点C1分别位于 长方体布局可行域的六个顶点处,所得到的管路布局形式。
本发明的进一步技术方案是:所述步骤二中,布局可行域内的三段式弯折布局方式:首先确定两个折弯点的坐标,折弯点C1坐标为(xC1,yC1,zC1),折弯点C2坐标为 (xC2,yC2,zC2);然后将A1点和C1点、C1点和C2点、C2点和B1点分别采用直线连接,即得 到三段式弯曲的管路布局形式。
本发明的进一步技术方案是:所述步骤二中,折弯点C1的坐标应满足:
Figure BDA0003135406250000051
折弯点C2的坐标应满足:
Figure BDA0003135406250000052
本发明的进一步技术方案是:所述步骤三中,建立管路模型后计算管路的固有频率;然后判断初始管路是否存在共振风险,如不存在风险则完成设计;如存在风险则 进行优化设计;
所述优化设计具体为:首先根据实际飞行器的振动频率环境,设计管路的固有频率目标,即错开飞行器的危险共振频率;然后对确定的两个折弯点进行多目标遗传优 化算法,获得优化后的两个折弯点的坐标。
有益效果
本发明的有益效果在于:本发明对于飞行器管路系统的任意两个出入口端口之间的管路布局设计提供了方法,发明了一段式~多段式弯折方案。能够对不同管路布局进 行分析和选择,便于工程中多种布局管路结构的快速建模和频率特性对比。具体效果 如下:
(1)发明中给出了管路入口、出口接头的过渡段的尺寸预留方法。同时考虑了该过渡段的方向(采用归一化的矢量形式),给出了过渡段终点的坐标表达式。
Figure BDA0003135406250000061
上式是考虑了出入口坐标、直线段要求、弯曲半径和各自引出引入方向四个因素的综合表达式,只要输入上述四个信息,就能直接确定布局可行域的对角线坐标。
(2)发明中将管路的布局空间进行了长方体量化,明确出布局可行域的概念,将管形布局优化的设计空间限定在布局可行域内,便于量化实施设计的开展。并且在布 局可行域内提出明确的弯折数量k与对角线A1 B1的长度相关,即
Figure BDA0003135406250000062
在布局可行域内,当2Rb≤A1B1<4Rb,只能采用一段弯折,为直接连接布局可行域的对角线。当4Rb≤A1B1<6Rb,可采用一段、两段弯折管形。当A1B1≥6Rb,可采用一段、 两段、三段弯折管形。弯折点约多,对于管路布局可行域的空间要求越大,管路的可 设计性越强,但是计算工作量也随之增大。设计人员可以根据该式,简便快速地确定 在可行域内的管路的可弯折数量。
(3)发明中针对两段弯折管形,根据布局可行域的边界提出了六种特殊的弯折点时布局形式,通过计算六种特殊的两段弯折管形布局的固有频率来评估两段弯折管形 的动力学特性,设计人员可以根据这六种特殊情况的固有频率判断两段弯折管形的频 率区间,合理选择弯折点的位置来满足不同的工况情形。
(4)三段弯折管形具有更好的动力学设计性,发明给出了三段弯折避开多个共振频率的弯折点位置优化方法。并且针对三段弯折管形中弯折点C1和C2计算公式复杂, 且直接优化时寻点效率低下、甚至不收敛等问题,根据各个弯折段的尺寸要求,发明 中给出了通过预先对C1点和C2进行划定域,即:
|xC1|≤l1-4Rb |xC2|≥l1-2Rb
|yC1|≤l2-4Rb,|yC2|≥l2-2Rb
|zC1|≤l3-4Rb |zC2|≥l3-2Rb
该方法实现了两个弯折点C1和C2点的快速寻优,大大提高了优化效率。并且能够获得多目标避开共振的管路动力学设计目标。
附图说明
图1入口、出口处需要预留的接口过渡段总长度示意图;
图2空间管路入口、出口和布局可行域示意图;
图3可行域内一段式弯曲的管路布局形式示意图;
图4可行域内两段式弯折的管路布局形式示意图;
图5六种特殊的C1点的位置时的两段弯曲管形布局示意图;
图6可行域内三段式弯折的管路布局形式示意图;
图7三段弯折管形中两个弯折点坐标位置的优化算法步骤图;
图8本发明实施例中空间管路入口、出口和布局可行域案例示意图;
图9本发明实施例中一段式弯折管路布局示意图;
图10本发明实施例中两段式弯折管路布局示意图;
图11本发明实施例中六种特殊的C1点的位置时的两段弯折管形布局示意图;
图12本发明实施例中二段式弯折管路布局示意图。
附图标记说明:A0为管路所要连接的入口的坐标位置,B0为管路所要连接的出口的坐标位置,A1为入口坐标位置A0经过过渡段后到达的坐标点,B1为出口坐标位置B0经 过过渡段后到达的坐标点,nA为管路所要连接的入口方向,nB为管路所要连接的出口 方向,Ls为管路接口过渡段的直线段,Rb为管路接口过渡段的圆弧段的弯曲半径,D0为管路外径。
具体实施方式
下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
本发明一种三维空间可行域内的管路弯折设计方法,在一个任意的三维直角坐标空间内,要实现合理的管路布局设计,具体步骤如下:
(1)确定布局可行域。通过管路所要连接的入口A0、出口B0的坐标位置和入口 方向nA,出口方向nB,再根据直线段要求Ls、弯曲半径Rb、管路外径D0等信息,确 定出管路布局的三维空间可行域。其中,布局可行域是以A1B1为对角线的长方体空间, A1和B1的坐标为:
Figure BDA0003135406250000081
管路布局可行域的长、宽、高分别为:
l1=|xA1-xB1|
l2=|yA1-yB1|
l3=|zA1-zB1|
(2)在布局可行域内确定在可行域内的弯折数量k,根据实际飞行器管路结构形式,主要分为一段弯折、两段弯折、三段弯折形式或者更多数量的弯折。确定的方法 需要根据A1B1的长度进行如下判别:
Figure BDA0003135406250000082
(3)根据布局可行域对角线的长度,选取不同的弯折数量,注意弯折点的确定应考虑弯曲角度和弯曲半径的限制,即可获得对应的布局管形。其中:当2Rb≤A1B1<4Rb, 只能采用一段弯折,为直接连接布局可行域的对角线。当4Rb≤A1B1<6Rb,可采用一段、 两段弯折管形,在可行域内确定一个弯折点。当A1B1≥6Rb,可采用一段、两段、三段 弯折管形,在可行域内确定两个弯折点。
(4)优化设计中首先将弯折点坐标参数化,并赋予初始值,建立上述不同布局的管路模型后,进一步计算分析管路结构的模态特性。根据实际飞行器的振动频率环境, 管路通过对弯折点的优化设计来错开飞行器的危险共振频率。确保管路能够获得符合 动力学特性的管形布局设计。
对于上述步骤得到的管路排布方式进行优化,得到最终优化排布结果。
下面依据上述四个步骤,详细说明设计方法:
(1)管路布局可行域的确定
航空航天管路设计标准规定,管路从端头沿着引出方向(介质传输方向)应具有一截接口过渡段,如图1所示,直线过渡段的总长度应包括两个部分,即直线段长度Ls和管路的弯曲半径Rb。所以接口过渡段的总长度要求为Ls+Rb
其中,直线段的长度Ls由各行业的设计标准而定,例如航空常标准规定为
Figure BDA0003135406250000091
式中,D0为所设计管道的外径。
同时,管路弯折时的弯曲半径Rb也有相关规定,一般而言,Rb取为四倍的管道外径。即
Rb=4D0 (2)
因此,在管路接口引出方向上应考虑接口过渡段所占据的长度,管路的布局设计应在除去接口过渡段之外的空间进行。如图2所示,假设两个设备的两个接口A0和B0的之间布置一段管路,其中,入口A0和出口B0坐标用向量分别表示为:
Figure BDA0003135406250000092
设计的管路从入口A0引入、从出口B0引出,为了避免管路的突然弯折,设计要求 出入口需要有一个过渡段。在两个接口各自的引出方向上,入口A0和出口B0经过过渡 段后分别到达A1和B1点,所以根据式(1),管路沿着接口A0和B0引出方向的长度,应 至少预留的接口过渡段长度A0 A1,B0 B1的最小长度应为:
|A0A1|=|B0B1|=Ls+Rb (4)
管路的实际布局空间应扣除各自的接口过渡段后,即在A1B1的范围内进行布局设计,以A1B1为对角线的长方体区域称为布局可行域。即图中的虚线长方体区域。若建 立图示以A1为坐标原点的直角坐标系,管路布局可行域的长、宽、高分别为:
Figure BDA0003135406250000101
为了确定布局可行域的界限,需要确定A1和B1点的坐标,此时需要考虑该过渡段的方向,这里令入口A0的直线引入段方向为nA,出口B0直线引出段的方向为nB。该方 向用归一化的列矢量形式表示为:
Figure BDA0003135406250000102
式中,i,j,k分别表示直角坐标系x,y,z三个轴正向单位矢。(aA,bA,cA)、(aB,bB,cB)分别表示入口和出口直线段的归一化方向。即
Figure BDA0003135406250000103
则A1和B1点的坐标采用向量表示如下:
rA1=rA0+(Ls+Rb)nA (8)
rB1=rB0-(Ls+Rb)nB (9)
具体的坐标表达式为:
Figure BDA0003135406250000104
(2)在可行域内设计弯折数量
一般而言,管路弯折数量与对角线A1 B1的长度相关,假设在布局可行域内确定在可行域内的弯折数量k,根据实际飞行器管路结构形式,主要分为一段弯折(k=1)、 两段弯折(k=2)、三段弯折形式(k=3),更多数量的弯折仅在复杂长管路中使用, 且需要配合布置卡箍。本发明仅讨论三段弯折以内的情形。确定的弯折数量的方法, 需要根据A1B1的长度进行如下判别:
Figure BDA0003135406250000105
其中,k=0表示可行域空间不足,无法在可行域中布局,需要将可行域进行拓展后再进行布局设计,此种情形将另做发明讨论,本发明仅讨论布局可行域的情形,从 式(11)可以看出,若越长,则可选择的弯折数量越多。下面具体对三种情况进行说 明如下:
(a)一段式弯折的管路布局形式(k=1)
一般而言,两点间的管路布局优先服从路径最短、弯折最少的要求。此时仅需将A1和B1点直接连接,形成直接一段式弯折的管路布局形式,如图3所示。这种布局方 式是接口端头A0和B0最直接的一种布局方式,路径最短且弯折最少。从动力学角度来 讲,也是刚度最大,质量最小,因此这种管形布局的基频最高。但是一段式弯折的管 路布局方式为固定(无法修改),因此无法使用该方案进行管路调频设计。
可以看到,对于图3中可行域内一段式弯折布局方式,它是将布局域的对角线A1 B1直接连接起来。由于弯折点A1和B1处存在弯曲半径Rb,因此,A1 B1的长度应大于2 倍的弯曲半径,即
A1B1≥2Rb (12) 。
(b)两段式弯折的管路布局形式(k=2)
布局可行域内两段式弯曲的管路布局方法,实际上是在一段式弯折的基础上,在长方体可行布局域内寻找一点C1,假设其坐标为(xC1,yC1,zC1),并将A1和B1均与C1点连接,这样就形成了两段式弯曲的管路布局形式,如图4所示。
两段式弯折的管路布局形式对布局可行域的尺寸有进一步的要求,即应满足A1C1和C1B1两段管路的弯曲半径的需要,即长方体布局可行域的长、宽、高均应大于2倍 的弯曲半径。由于弯折点A1、B1存在弯曲半径Rb,且C1处有两个弯曲半径。因此,A1 B1的长度至少应大于4倍的弯曲半径,才能符合在布局可行域实施两折线弯曲布局的条 件。即
A1B1≥4Rb (13) 。
另外,C1点的位置确定方法,可以依据管路系统避开共振频率的要求进行优选,通过优化设计计算获得。为了快速确定不同位置C1点对管路系统动态特性的影响,可 以先采用六种特殊位置情况,对一阶固有频率做一个评估。这六种特殊位置即是C1点 分别在长方体布局可行域的三个顶点处,如图5所示,其中(a)(d)为弯折最少,(c) (f)为弯曲角度最小(所有弯折均为90°)。显然,C1点的位置若在可行域之内,其 一阶固有频率必然小于上述六种情形,且大于一段弯折的情形。
(c)三段式弯折的管路布局形式(k=3)
如图8所示,所谓三段式弯折的管路布局,即在可行域内有三段管路。实际上是 在两段式弯折的基础上再增加一个点C2,假设其坐标是(xC2,yC2,zC2)。由于弯折点 A1、B1和C1处存在弯曲半径Rb,因此,各个弯折段的长度均应大于2倍的弯曲半径, 即A1 B1的长度应大于等于6倍的弯曲半径,才能符合布局可行域条件。
A1B1≥6Rb (14)
对比式(12)、式(13)和式(15)可以看到,弯折数量越多,对可行域的尺寸要 求越大,因此,实际管路布局设计中应根据可行域的空间尺寸,确定合理的弯折数量。
C1和C2点应根据管路的引入方向和引出方向,确在一定范围内,这样能够保证飞行器管路设计中弯折角度应≥90度。在图6所示的三折弯曲走向布局形式中,为了避 免管形干涉,C1和C2点应满足以下条件:
Figure BDA0003135406250000121
式(15)用坐标形式表示为:
Figure BDA0003135406250000131
Figure BDA0003135406250000132
Figure BDA0003135406250000133
Figure BDA0003135406250000134
Figure BDA0003135406250000135
式(16)计算公式复杂,直接采用优化确定C1和C2点的方法,容易造成寻点效率 低下、甚至不收敛等问题,较为有效的方法是首先给C1点划定域,由于要给另外两个 弯折段留有足够的布局空间,C1点的坐标布置范围应该在区间内:
Figure BDA0003135406250000136
同理,C2点的坐标应满足:
Figure BDA0003135406250000137
综上,弯折数量与可行域对角线的长度相关。当2Rb≤A1B1<4Rb,只能采用一段弯折,为直接连接布局可行域的对角线。当4Rb≤A1B1<6Rb,可采用一段、两段弯折管形, 在可行域内确定一个弯折点。当A1B1≥6Rb,可采用一段、两段、三段弯折管形,在可 行域内确定两个弯折点,该情形可以实施弯折点的优化设计。
综合所需的优化目标(动力学设计一般要求避开系统的工况运行频率),再约束条件式(17)、(18),采用多目标遗传算法,可以通过优化计算获得两个弯折点的最优位 置。具体方法如图7所示。
实施例:
本发明一种三维空间可行域内的管路弯折设计方法的具体实施例,如图8所示,若要在入口A0和出口B0之间设计管路,分别采用本发明提出的一段式弯折、两段式弯 折和三段式弯折方案进行设计。这里假设入口A0的坐标为(0,0,0),引入过渡段的 方向为
Figure BDA0003135406250000138
出口B0的坐标为(500,400,300),引入过渡段的方向为
Figure BDA0003135406250000139
坐标单位为mm。需设计的管路直径D为12mm,弯曲半径Rb为4倍的管径,即48mm, 直线段Ls的长度要求为16mm。
根据本发明,首先根据式(10)确定出A1和B1,如图6所示,以A1和B1为对角线 的长方体区域即为管路的布局可行域。
Figure BDA0003135406250000141
然后根据式(11)确定可选的弯折数量,对比可知:
Figure BDA0003135406250000142
因此,在该布局可行域内,可以实施一段弯折~三段弯折的布局方式,通过本发明提出的管形布局设计方法,可以设置不同弯折数量后,即可得到对应的管形布局。下 面列出了一段弯折~三段弯折时,管形布局结果。
(a)一段式弯折布局
根据图3所示的原理,将本例中布局可行域的对角线A1B1直接连接,并设置给定 的弯曲半径Rb,即可得到一段式弯折布局管路,如图9所示。一段弯折管形无设计点, 结构形式固定不可调,对应的动力学特性也是固定的,计算图9所示的一段弯折管路 结构的固有频率如表1所示。
(b)两段式弯折布局
根据图4所示的原理和式(9)中对C1点的位置要求,对于本例中,C1点的选取在 布局可行域范围内即可,即:
64≤xC1≤500
0≤yC1≤336
-300≤zC1≤0
在该范围内选择任意一点,即可形成一段式弯折管路布局,如图10所示
按照图5所示的六种特殊的C1点的位置时的两段弯折管形布局设计方法,获得对应于本例所适用的六种两段式弯折布局如图9所示。这六种管形布局对应的固有频率 特性如表1.设计人员可以根据这六种特殊情况判断两段弯折管形的频率区间。并选择 合理的管形布局来满足不同的情形。例如要避开液压泵的脉动频率[550,600]区间,可 以选用(a)、(c)、(d),若要避开发动机转子叶片的工作频率约186Hz,可以选用管 形(f)。
表1上述六种特殊的C1点的位置时的两段弯折管形布局的固有频率
Figure BDA0003135406250000151
(c)三段式弯折布局
根据图5所示的原理和式(10)、式(11)中对C1点和C2点的位置要求,对于本 例中,若在前述C1点的基础上进一步选取C2点,即可形成一段式弯折管路布局,如图 12所示。
(d)不同管形布局下结构的固有频率特性对比分析
将上述三种管形的入口A0和出口B0固定,计算管路的固有频率,表2列出了上述 三种管形的各点坐标和固有频率。能够看到,采用上述不同管形布局时,对应的管路 结构的固有频率也有所不同,例如若要提高管路的一阶基频,可以使用一段式弯折布 局。若要避开400Hz~460Hz的频带,则应选择三段式弯折管路布局,可见,本发明提 出的三维空间可行域内的管路布局设计方法为管路的避开共振的设计提供了参考。能 够在两个出入端口坐标位置不变的情形下,通过管形设计来实现错开共振频率并降低 结构振动应力,提高了管路结构的振动可靠性。
表2三种管形的各点坐标和固有频率
Figure BDA0003135406250000152
Figure BDA0003135406250000161
(e)三段弯折管形中两个弯折点的优化设计
若表2中三段弯折管形给出的两个初始弯折点C1(300,136,100)和C2 (360,256,200)的坐标无法满足动力学设计要求,例如对于发动机管路要求管路的固 有频率要避开典型的外部激振频率,假设管路要求避开的两个外部激励频率范围:一 是风扇转动激励频率区间[100,220]Hz,二是液压泵的脉动工作频率[540,620]Hz。 可以采用本发明提出的三段弯折管形中两个弯折点坐标位置的优化算法进行优化设 计,将避开上述两个频率区间作为设计目标,即
f2≤100
f3≥220
f5≤520
f6≥620
优化计算后得到的两个弯折点C1和C2的坐标分别为:
xC1=149.49xC2=323.6
yC1=6.48,yC2=377.79
zC1=22.68zC2=281.6
将上述最优坐标带入建模对比计算,可以看到三段弯折管路对弯折点优化前后的效果和管形对比,如表3所示。可以看到,优化后的管形各阶固有频率均错开了两个 激振频段,完成了最优设计管形结构设计。
表3三段弯折管形动力学特性优化前后对比
Figure BDA0003135406250000162
Figure BDA0003135406250000171
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和 宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (7)

1.一种三维空间可行域内的管路弯折设计方法,其特征在于具体步骤如下:
步骤一:确定布局可行域;
首先,通过确定管路所要连接的入口A0、出口B0的坐标位置和入口方向nA,出口方向nB;然后,依据航空航天管路设计标准,确定入口A0、出口B0处需要预留的接口过渡段,接口过渡段的总长度为Ls+Rb,其中,Ls为管路接口过渡段的直线段,Rb为管路接口过渡段的圆弧段的弯曲半径;
最后,确定出管路布局的三维空间可行域,其中,布局可行域是以A1B1为对角线的长方体空间,A1和B1的坐标为:
Figure FDA0003928622460000011
管路布局可行域的长、宽、高分别为:
Figure FDA0003928622460000012
步骤二:在布局可行域内确定管路弯折数量k,具体根据A1B1的长度进行确定:
Figure FDA0003928622460000013
Figure FDA0003928622460000014
其中,k=0表示可行域空间不足,无法在可行域中布局;当2Rb≤A1B1<4Rb,采用一段弯折,管路为直接连接布局可行域的对角线;当4Rb≤A1B1<6Rb,采用一段、两段弯折管形,管路在布局可行域内确定一个弯折点;当A1B1≥6Rb,采用一段、两段、三段弯折管形,管路在布局在可行域内确定两个弯折点;
步骤三:首先,将步骤二确定的弯折点坐标参数化,并赋予初始值,建立不同布局的管路模型后,进一步计算分析管路结构的模态特性;然后,根据实际飞行器的振动频率环境,管路通过对弯折点的优化设计来错开飞行器的危险共振频率;最后,确保管路能够获得符合动力学特性的管形布局设计;
所述步骤一中,为了确定布局可行域的界限,需要确定A1和B1点的坐标,并考虑该接口过渡段的方向,其中,入口A0的直线引入段方向为nA,出口B0直线引出段的方向为nB,该方向用归一化的列矢量形式表示为:
Figure FDA0003928622460000021
式中,i,j,k分别表示直角坐标系x,y,z三个轴正向单位矢,(aA,bA,cA)、(aB,bB,cB)分别表示入口和出口直线段的归一化方向,即
Figure FDA0003928622460000022
则A1和B1点的坐标采用向量表示如下:
rA1=rA0+(Ls+Rb)nA
rB1=rB0-(Ls+Rb)nB
得到坐标表达式为:
Figure FDA0003928622460000023
2.根据权利要求1所述三维空间可行域内的管路弯折设计方法,其特征在于:所述步骤二中,布局可行域内的一段式弯折布局方式为,将布局可行域的对角线A1 B1采用直线连接。
3.根据权利要求1所述三维空间可行域内的管路弯折设计方法,其特征在于:所述步骤二中,布局可行域内的两段式弯折布局方式:首先确定折弯点C1坐标为(xC1,yC1,zC1),然后将A1、B1分别与C1点直线连接,即得到两段式弯曲的管路布局形式。
4.根据权利要求3所述三维空间可行域内的管路弯折设计方法,其特征在于:所述步骤二中,折弯点C1的位置确定方法:采用六种特殊位置情况,对一阶固有频率做评估;所述六种特殊位置是指折弯点C1分别位于长方体布局可行域的六个顶点处,所得到的管路布局形式。
5.根据权利要求1所述三维空间可行域内的管路弯折设计方法,其特征在于:所述步骤二中,布局可行域内的三段式弯折布局方式:首先确定两个折弯点的坐标,折弯点C1坐标为(xC1,yC1,zC1),折弯点C2坐标为(xC2,yC2,zC2);然后将A1点和C1点、C1点和C2点、C2点和B1点分别采用直线连接,即得到三段式弯曲的管路布局形式。
6.根据权利要求5所述三维空间可行域内的管路弯折设计方法,其特征在于:所述步骤二中,折弯点C1的坐标应满足:
Figure FDA0003928622460000031
折弯点C2的坐标应满足:
Figure FDA0003928622460000032
7.根据权利要求1所述三维空间可行域内的管路弯折设计方法,其特征在于:所述步骤三中,建立管路模型后计算管路的固有频率;然后判断初始管路是否存在共振风险,如不存在风险则完成设计;如存在风险则进行优化设计;
所述优化设计具体为:首先根据实际飞行器的振动频率环境,设计管路的固有频率目标,即错开飞行器的危险共振频率;然后对确定的两个折弯点进行多目标遗传优化算法,获得优化后的两个折弯点的坐标。
CN202110717390.4A 2021-06-28 2021-06-28 一种三维空间可行域内的管路弯折设计方法 Active CN113434959B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110717390.4A CN113434959B (zh) 2021-06-28 2021-06-28 一种三维空间可行域内的管路弯折设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110717390.4A CN113434959B (zh) 2021-06-28 2021-06-28 一种三维空间可行域内的管路弯折设计方法

Publications (2)

Publication Number Publication Date
CN113434959A CN113434959A (zh) 2021-09-24
CN113434959B true CN113434959B (zh) 2023-03-03

Family

ID=77754873

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110717390.4A Active CN113434959B (zh) 2021-06-28 2021-06-28 一种三维空间可行域内的管路弯折设计方法

Country Status (1)

Country Link
CN (1) CN113434959B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3900990A1 (de) * 1989-01-12 1990-08-02 Bundy Gmbh Verfahren zum rechnergestuetzten auslegen des verlaufs, zum herstellen und einrichtung zum herstellen von raeumlich winklig gebogenen, kleinkalibrigen rohrleitungen, insbes. von kraftstoff- oder bremsleitungen im kraftfahrzeugbau
CN102262704A (zh) * 2011-08-19 2011-11-30 泰戈特(北京)工程技术有限公司 管路系统布置图的生成方法和装置
CN109740256A (zh) * 2019-01-02 2019-05-10 北京航空航天大学 一种融合装备维修性的线缆布局设计优化方法
CN110320931A (zh) * 2019-06-20 2019-10-11 西安爱生技术集团公司 基于航向控制律的无人机避障航路规划方法
CN110442940A (zh) * 2019-07-25 2019-11-12 天津大学青岛海洋技术研究院 一种高效的船舶多管路并行智能布局设计方法
CN111597648A (zh) * 2020-04-23 2020-08-28 西北工业大学 一种航空发动机机匣弧面管路的参数化建模方法
CN111666637A (zh) * 2020-06-29 2020-09-15 江苏科技大学 一种基于知识的船舶管路布局优化设计方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2829993B1 (en) * 2013-07-25 2020-09-30 Dassault Systèmes Design of a path connecting a first point to a second point in a three-dimensional scene
CN107985626B (zh) * 2017-10-27 2021-02-19 南京航空航天大学 一种基于可变构型空天飞行器的气动布局设计方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3900990A1 (de) * 1989-01-12 1990-08-02 Bundy Gmbh Verfahren zum rechnergestuetzten auslegen des verlaufs, zum herstellen und einrichtung zum herstellen von raeumlich winklig gebogenen, kleinkalibrigen rohrleitungen, insbes. von kraftstoff- oder bremsleitungen im kraftfahrzeugbau
CN102262704A (zh) * 2011-08-19 2011-11-30 泰戈特(北京)工程技术有限公司 管路系统布置图的生成方法和装置
CN109740256A (zh) * 2019-01-02 2019-05-10 北京航空航天大学 一种融合装备维修性的线缆布局设计优化方法
CN110320931A (zh) * 2019-06-20 2019-10-11 西安爱生技术集团公司 基于航向控制律的无人机避障航路规划方法
CN110442940A (zh) * 2019-07-25 2019-11-12 天津大学青岛海洋技术研究院 一种高效的船舶多管路并行智能布局设计方法
CN111597648A (zh) * 2020-04-23 2020-08-28 西北工业大学 一种航空发动机机匣弧面管路的参数化建模方法
CN111666637A (zh) * 2020-06-29 2020-09-15 江苏科技大学 一种基于知识的船舶管路布局优化设计方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
"?"弯曲管形在飞行器高压高速流体管道系统中应用的力学机理探讨;张雅玫等;《应用力学学报》;20200815(第04期);第48-53页 *
FSI research in pipeline systems – A review of the literature;Li, SJ;Karney, BW;Liu, GM;《Journal of Fluids & Structures》;20150725;第57卷;第277-297页 *
三维管线布局优化;王志红等;《计算机辅助设计与制造》;19971117(第11期);第22-25页 *
基于NSGA-Ⅱ的三维空间直角管路布局优化;焦国帅等;《控制工程》;20181120(第11期);第108-113页 *

Also Published As

Publication number Publication date
CN113434959A (zh) 2021-09-24

Similar Documents

Publication Publication Date Title
Muralikrishna et al. Heat exchanger design targets for minimum area and cost
CN111666637A (zh) 一种基于知识的船舶管路布局优化设计方法
CN109977475B (zh) 一种用于三维复杂管路流固耦合计算的动网格更新方法
CN101567026B (zh) 一种输液管道导流结构减振弯头优化设计方法
CN113434959B (zh) 一种三维空间可行域内的管路弯折设计方法
CN112613134B (zh) 一种基于涡流分布的阀体结构优化方法
TW201809574A (zh) 對於流體分佈網路具有減少之壓力損失的u形彎管(u-bends)
CN108062444B (zh) 弯管参数确定方法
CN111597648B (zh) 一种航空发动机机匣弧面管路的参数化建模方法
CN116380472B (zh) 一种大涵道比发动机核心机试验中的进气装置
CN109614747B (zh) 大展弦比机翼中管道布局设计方法
CN113221354A (zh) 一种管道弯曲变形拟合算法
Gupta et al. Symmetrical laminar channel flow with wall suction
EP1563216A2 (en) A flexible, tubular device e.g. a bellows
CN116541949A (zh) 一种在拓展域内管路布局设计方法
CN110005895B (zh) 变径弯管及利用变径弯管的管道布置连接方法
Kyaw Oo D’Amore et al. Numerical study on modelling perforated elements using porous baffle interface and porous region
CN108870086B (zh) 一种长直集输管道热位移综合控制系统及方法
CN112347583A (zh) 一种增压站含双缺陷弯管极限内压的计算方法
CN112228314B (zh) 排气管配管、衍生结构、客车空调压缩机及仿真方法
CN1545616A (zh) 基于闭式波纹管刚度确定的模型的方法及装置
Peiris et al. Prediction of Flow Path Pressure Drops in Curved Galleries for Additively Manufactured Hydraulic Manifolds
Singh et al. Design the Shell and Tube Heat Exchanger with the help of Programming using MATLAB software
CN219912193U (zh) 一种用于trt或bprt的进出口管道结构
CN201122816Y (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
GR01 Patent grant