CN110119560B - 一种基于曲率半径校核通信卫星软波导强度的方法 - Google Patents

一种基于曲率半径校核通信卫星软波导强度的方法 Download PDF

Info

Publication number
CN110119560B
CN110119560B CN201910364758.6A CN201910364758A CN110119560B CN 110119560 B CN110119560 B CN 110119560B CN 201910364758 A CN201910364758 A CN 201910364758A CN 110119560 B CN110119560 B CN 110119560B
Authority
CN
China
Prior art keywords
waveguide
curvature radius
soft
checking
soft waveguide
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
CN201910364758.6A
Other languages
English (en)
Other versions
CN110119560A (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.)
China Academy of Space Technology CAST
Original Assignee
China Academy of Space Technology CAST
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 China Academy of Space Technology CAST filed Critical China Academy of Space Technology CAST
Priority to CN201910364758.6A priority Critical patent/CN110119560B/zh
Publication of CN110119560A publication Critical patent/CN110119560A/zh
Application granted granted Critical
Publication of CN110119560B publication Critical patent/CN110119560B/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/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Waveguides (AREA)
  • Details Of Aerials (AREA)

Abstract

本发明公开了一种基于曲率半径校核通信卫星软波导强度的方法，(1)建立软波导有限元模型；(2)选取软波导节点；(3)在整星有限元模型上施加激励载荷，并计算所选取的软波导节点在每一个频点或每一个时间步的位移响应；(4)在所选取节点的初始坐标上叠加步骤(3)计算得到的位移响应，得到承受激励时所选择节点的坐标位置；(5)计算曲率半径；(6)根据产品使用场景校核软波导的强度是否符合安全裕度要求。本发明可以充分利用位移计算精度高的优点，有效规避了位移一次求导后计算得到的应力精度低、应力量级偏大的缺陷，可有效提高校核的准确性。

Description

一种基于曲率半径校核通信卫星软波导强度的方法
技术领域
本发明涉及一种基于曲率半径校核通信卫星软波导强度的方法，属于机械技术领域。
背景技术
软波导是一种常用的通信卫星星上波导形式，其结构构型为由一组皱褶环依序连接构成的管状结构(如图1所示，主体构型类似于波纹管)，通常用来实现跨舱段的波导连接。为保证卫星功能性能的完好性，在整星波导布局和力学试验前，需要对软波导强度进行校核，避免软波导因承受力学环境载荷而开裂损坏，进而影响整星功能和性能。
当前通信卫星对于软波导的强度校核主要是基于类比方法和强度分析方法。类比方式主要是对以往星上软波导的环境试验数据进行梳理，并与拟使用的软波导进行对比分析，评估其安全性。但是由于软波导的截面尺寸、长度、所处舱段位置等因素均存在较大的差别，这一方法的准确度较差。强度分析方式主要是建立软波导的准确有限元模型，将其带入整星模型进行分析，计算软波导上的应力并与其许用应力进行比较，校核其强度，当应力计算较为准确时，这一方法的精度远高于前述类比方法。但是根据力学常识，应力为位移的一阶导数与弹性模量的乘积，由于软波导本身为皱褶环依序连接构成的管状结构，当网格稀疏时，虽然根据当前有限元软件普遍采用的最小位移变分原理，可以计算求解得到较为准确的位移值，但是其导数极有可能不连续，导致对其取一阶导数计算得到的应力畸变严重，计算值将严重偏大。若要精确计算其应力水平，则需要极其密集的有限元网格。由于计算量通常与网格密度的平方呈正比，计算量将极其庞大。综上所述，当前对软波导尽享强度校核的方法存在着不够准确和/或计算量难以承受的缺陷。
根据材料力学常识，结构内部某区域的应力与其局部的曲率半径有关，但是由于工程结构的复杂性，在一般的工程结构设计实践及各种相关标准中，没有以曲率半径形式给出的强度校核准则。一个原因是对于复杂工程结构，计算结构的曲率半径的复杂性远远超出了计算其应力的复杂性；另一个原因是许用强度仅取决于材料本身的特性，而许用曲率半径不仅取决于材料本身的特性，还与结构形式息息相关，不仅同样材料构成的不同结构的许用曲率半径不同，而且同一个结构内部的局部结构形式不同，也会导致不同局部结构的许用曲率半径不同。
发明内容
本发明的技术解决问题：为克服现有技术的不足，提供一种基于曲率半径校核通信卫星软波导强度的方法，有效规避了位移一次求导后计算应力精度低、应力量级偏大的缺陷，可有效提高校核的准确性。
本发明的技术解决方案：
一种基于曲率半径校核通信卫星软波导强度的方法，具体步骤为：
(1)建立软波导有限元模型：有限元模型中节点的排列方式与软波导母线方向平行，然后将软波导有限元模型导入已建立好的整星有限元模型中；
(2)选取软波导节点：在每一个皱褶环上各选取一个节点，形成一个节点组，并记录节点组中所有节点的初始坐标位置；
(3)在整星有限元模型上施加激励载荷，并计算所选取的软波导节点在每一个频点或每一个时间步的位移响应；
(4)在所选取节点的初始坐标上叠加步骤(3)计算得到的位移响应，得到承受激励时所选择节点的坐标位置；
(5)根据数学常识，任意3个点可以在三维空间唯一确定一个圆，则在第(2)步所选取的节点组中，基于每三个相邻的软波导节点，根据其承受激励时的坐标位置，即可计算得到其曲率半径；
(6)按照下式计算公式，根据产品使用场景校核软波导的强度是否符合安全裕度要求，
Figure BDA0002047828340000031
上式中，MS为安全裕度；rs为最小曲率半径，为步骤(5)中计算得到的所有曲率半径中的最小值；rf为许用曲率半径，由软波导自身特性决定；fs为安全系数，根据产品使用场景确定。
步骤(1)中模型上的每一个皱褶环沿波导轴线方向至少需要用4个单元表征其力学特性。
步骤(1)中模型上的每一个皱褶环用8个单元表征其力学特性。
步骤(2)中该节点组中所有节点均应位于一条与母线方向平行的直线上。
本发明与现有技术相比具有如下有益效果：
本发明提出了一种首先利用软波导承受载荷激励后在星上位移的计算其曲率半径，再用曲率半径校核软波导强度的方法。一方面，由于软波导结构形式相对固定，为由一组皱褶环依序连接构成的类似于波纹管的管状结构，且在卫星工程实践中通常直接选用型谱化产品，可以通过力学试验或分析，明确确定其许用的曲率半径，替代直接校核其材料的许用强度的方法，作为校核软波导强度的准则，存在校核技术路线的可行性；另一方面，这一方法可以充分利用位移计算精度高的优点，有效规避了位移一次求导后计算应力精度低、应力量级偏大的缺陷，可有效提高校核的准确性。
附图说明
图1为本发明适用的软波导构型示意图；
图2为软波导截面定义示意图；
图3为本发明的基于曲率半径校核软波导强度的流程图；
图4为基于相邻的三个点求解曲率半径的示意图；
图5为软波导有限元模型的示意图；
图6为选取软波导有限元模型上一组与轴线方向平行的节点的示意图；
图7某型号通信卫星软波导安装位置示意图；
具体实施方式
下面结合附图对本发明作进一步详细的描述。
本发明的技术解决方案主要包括：选取节点、计算相对位移、计算承受载荷后的坐标位置、计算曲率、校核曲率，具体实现步骤如下：
一种基于曲率半径校核通信卫星软波导强度的方法，具体步骤为：
(1)建立软波导有限元模型，有限元模型中节点的排列方式需与软波导母线方向平行，如图5和图6所示，然后将软波导有限元模型代入已建立好的整星有限元模型中；
(2)选取软波导节点：在每一个皱褶环上各选取一个节点，形成一个节点组，如图5所示，并记录节点组中所有节点的初始坐标位置；
(3)在整星有限元模型上施加激励载荷，并计算所选取的软波导节点在每一个频点或每一个时间步的位移响应；
(4)在所选取节点的初始坐标上叠加步骤(3)计算得到的位移响应，得到承受激励时所选择节点的坐标位置；
(5)根据数学常识，任意3个点可以在三维空间唯一确定一个圆，则按照图4所示的方法，在第(2)步所选取的节点组中，基于每三个所选取的相邻的软波导节点，根据其承受激励时的坐标位置，即可计算得到其曲率半径；
(6)按照下式所示的计算公式，根据产品使用场景校核软波导的强度是否符合安全裕度要求。
Figure BDA0002047828340000041
上式中，MS为安全裕度；rs为最小曲率半径，为第(5)步骤中计算得到的所有曲率半径中的最小值；rf为许用曲率半径，由软波导自身特性决定；fs为安全系数，根据产品使用场景确定。
步骤(1)中模型上的每一个皱褶环沿波导轴线方向至少需要用4个单元表征其力学特性，优选8个单元表征其力学特性。
步骤(2)中该节点组中所有节点均应位于一条与母线方向平行的直线上。
实施例
对于某型号通信卫星，卫星上安装了四根软波导，如图7所示。其具体实施方式：
1)首先逐一建立软波导的有限元模型，然后将软波导模型代入整星模型；
2)在每根软波导上选取一组与母线方向平行的节点，并记录其初始坐标位置；
3)在整星有限元模型上分别沿X、Y、Z方向施加表1所示的激励载荷，并计算所选取的软波导节点在每一个频点的位移响应；
表1激励载荷
Figure BDA0002047828340000051
4)将计算得到的位移响应叠加到每一个节点的初始坐标上，获取其在激励作用下的坐标位置；
5)计算承受激励时所选取的每一个节点处的曲率半径；
6)按照本专利方法校核曲率半径安全裕度如表2所示(表中E截面和H截面定义如图2所示)，校核时的许用曲率半径由软波导承研厂家给出，安全系数取1.5，要求安全裕度大于0。由表可见，最小安全裕度为0.72，满足大于0的使用判据要求，校核结果为软波导强度可以满足要求。采用同样的有限元模型，计算得到的最大应力及与许用强度强度的校核结果如表3所示，校核时的许用应力由软波导承研厂家给出，安全系数取1.5，要求安全裕度大于0。由表可见，最小安全裕度为-0.89，不满足大于0的使用判据要求，校核结果为软波导强度不满足要求。
表2基于曲率半径的软波导强度校核表
Figure BDA0002047828340000061
表3基于应力的软波导强度校核表
Figure BDA0002047828340000062
Figure BDA0002047828340000071
7)上述软波导随整星参与了力学试验考核，并顺利通过，验证了本专利提出方法的有效性。
本发明提出了一种首先利用软波导承受载荷激励后在星上位移的计算其曲率半径，再用曲率半径校核软波导强度的方法。一方面，由于软波导结构形式相对固定，为由一组皱褶环依序连接构成的类似于波纹管的管状结构，且在卫星工程实践中通常直接选用型谱化产品，可以通过力学试验或分析，明确确定其许用的曲率半径，替代直接校核其材料的许用强度的方法，作为校核软波导强度的准则，存在校核技术路线的可行性；另一方面，这一方法可以充分利用位移计算精度高的优点，有效规避了位移一次求导后计算应力精度低、应力量级偏大的缺陷，可有效提高校核的准确性。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。

Claims (4)

1.一种基于曲率半径校核通信卫星软波导强度的方法，其特征在于，具体步骤为：
(1)建立软波导有限元模型：有限元模型中节点的排列方式与软波导母线方向平行，然后将软波导有限元模型导入已建立好的整星有限元模型中；
(2)选取软波导节点：在每一个皱褶环上各选取一个节点，形成一个节点组，并记录节点组中所有节点的初始坐标位置；
(3)在整星有限元模型上施加激励载荷，并计算所选取的软波导节点在每一个频点或每一个时间步的位移响应；
(4)在所选取节点的初始坐标上叠加步骤(3)计算得到的位移响应，得到承受激励时所选择节点的坐标位置；
(5)根据数学常识，任意3个点可以在三维空间唯一确定一个圆，则在第(2)步所选取的节点组中，基于每三个相邻的软波导节点，根据其承受激励时的坐标位置，即可计算得到其曲率半径；
(6)按照下式计算公式，根据产品使用场景校核软波导的强度是否符合安全裕度要求，
Figure FDA0002047828330000011
上式中，MS为安全裕度；rs为最小曲率半径，为步骤(5)中计算得到的所有曲率半径中的最小值；rf为许用曲率半径，由软波导自身特性决定；fs为安全系数，根据产品使用场景确定。
2.如权利要求1所述的一种基于曲率半径校核通信卫星软波导强度的方法，其特征在于，步骤(1)中模型上的每一个皱褶环沿波导轴线方向至少需要用4个单元表征其力学特性。
3.如权利要求2所述的一种基于曲率半径校核通信卫星软波导强度的方法，其特征在于，步骤(1)中模型上的每一个皱褶环用8个单元表征其力学特性。
4.如权利要求1所述的一种基于曲率半径校核通信卫星软波导强度的方法，其特征在于，步骤(2)中该节点组中所有节点均应位于一条与母线方向平行的直线上。
CN201910364758.6A 2019-04-30 2019-04-30 一种基于曲率半径校核通信卫星软波导强度的方法 Active CN110119560B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910364758.6A CN110119560B (zh) 2019-04-30 2019-04-30 一种基于曲率半径校核通信卫星软波导强度的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910364758.6A CN110119560B (zh) 2019-04-30 2019-04-30 一种基于曲率半径校核通信卫星软波导强度的方法

Publications (2)

Publication Number Publication Date
CN110119560A CN110119560A (zh) 2019-08-13
CN110119560B true CN110119560B (zh) 2022-10-11

Family

ID=67520347

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910364758.6A Active CN110119560B (zh) 2019-04-30 2019-04-30 一种基于曲率半径校核通信卫星软波导强度的方法

Country Status (1)

Country Link
CN (1) CN110119560B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101834331A (zh) * 2009-03-10 2010-09-15 中国电子科技集团公司第二十三研究所 3公分压窄型截面软波导及其制造方法
CN107085631A (zh) * 2017-03-30 2017-08-22 中国航空工业集团公司西安飞机设计研究所 一种基于舱门细节模型的强度校核方法
WO2018028284A1 (zh) * 2016-08-09 2018-02-15 苏州数设科技有限公司 飞机结构强度模型的创建方法及装置
CN109684693A (zh) * 2018-12-12 2019-04-26 中国航空工业集团公司西安飞机设计研究所 一种基于有限元分析预计加筋壁板后屈曲的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101834331A (zh) * 2009-03-10 2010-09-15 中国电子科技集团公司第二十三研究所 3公分压窄型截面软波导及其制造方法
WO2018028284A1 (zh) * 2016-08-09 2018-02-15 苏州数设科技有限公司 飞机结构强度模型的创建方法及装置
CN107085631A (zh) * 2017-03-30 2017-08-22 中国航空工业集团公司西安飞机设计研究所 一种基于舱门细节模型的强度校核方法
CN109684693A (zh) * 2018-12-12 2019-04-26 中国航空工业集团公司西安飞机设计研究所 一种基于有限元分析预计加筋壁板后屈曲的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
折叠展开波导在外力作用下横截面的可能性形变力学分析;陈霄等;《机械》;20080831;第35卷(第08期);第41-43页 *

Also Published As

Publication number Publication date
CN110119560A (zh) 2019-08-13

Similar Documents

Publication Publication Date Title
CN109858112B (zh) 基于结构应力监测结果的数值反演分析方法
CN108121865B (zh) 一种基于伴随变量法的含有多不确定参数的结构区间响应传播分析方法
CN112784356A (zh) 一种飞机结构静强度的试验方法
CN109359336A (zh) 一种基于多目标优化的绑扎桥相似畸变模型构建方法
CN115455793A (zh) 基于多尺度模型修正的高层结构复杂构件受力分析方法
CN110110406B (zh) 一种基于Excel计算平台实现LS-SVM模型的边坡稳定性预测方法
CN104424380A (zh) 机械变形测量仪的模拟
CN117421701B (zh) 一种自升式平台桩腿三维空间姿态分布式监测方法
CN110119560B (zh) 一种基于曲率半径校核通信卫星软波导强度的方法
CN107609234B (zh) 基于不确定性的稳健性分析方法及系统
CN110580391B (zh) 一种柔性结构的基频模态测量方法
CN111159936B (zh) 一种基于广义时域有限差分的电缆接头热场计算方法
Cary et al. Current practice unstructured grid CFD results for 3rd AIAA High Lift Prediction Workshop
Chwalowski et al. FUN3D analyses in support of the First Aeroelastic Prediction Workshop
CN116227045B (zh) 一种结构试件的局部应力应变场构造方法及系统
CN115357976B (zh) 基于结构性能的大跨度钢结构施工方案的迭代计算方法
CN114254533B (zh) 考核疲劳振动对产品组部件固定角度影响和预测的方法
CN115374667A (zh) 基于钢结构施工变形对结构设计性能影响的逆分析方法
CN110733665B (zh) 一种水面飞行器船底水载荷的三维预置加载方法
CN113139316A (zh) 一种利用频响函数指标辨识二维梁结构刚度薄弱环节的方法
CN112347686A (zh) 一种仪器安装结构的螺栓校核方法
Wieseman Structural dynamics modeling of HIRENASD in support of the aeroelastic prediction workshop
CN112182739A (zh) 一种飞行器结构非概率可信可靠性拓扑优化设计方法
CN114818445B (zh) 一种飞机结构静力叠加振动试验中静力加载点确定方法
CN110686615B (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