CN106557599A - 一种空间力系载荷谱编制方法 - Google Patents

一种空间力系载荷谱编制方法 Download PDF

Info

Publication number
CN106557599A
CN106557599A CN201510633822.8A CN201510633822A CN106557599A CN 106557599 A CN106557599 A CN 106557599A CN 201510633822 A CN201510633822 A CN 201510633822A CN 106557599 A CN106557599 A CN 106557599A
Authority
CN
China
Prior art keywords
load
testpieces
coordinate
sigma
pressure heart
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
CN201510633822.8A
Other languages
English (en)
Other versions
CN106557599B (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.)
AVIC Aircraft Strength Research Institute
Original Assignee
AVIC Aircraft Strength Research Institute
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 AVIC Aircraft Strength Research Institute filed Critical AVIC Aircraft Strength Research Institute
Priority to CN201510633822.8A priority Critical patent/CN106557599B/zh
Publication of CN106557599A publication Critical patent/CN106557599A/zh
Application granted granted Critical
Publication of CN106557599B publication Critical patent/CN106557599B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design 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)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明提供了一种空间力系载荷谱编制方法,其特征在于,包括以下步骤:步骤一、载荷分区简化令Pix、Piy和Piz依次为节点i在X向、Y向和Z向上的气动力,根据提供的分布节点载荷分别计算出n种载荷工况下X向载荷、Y向载荷和Z向载荷的总合力及相应压心坐标,对于n种载荷工况,m个气动载荷节点在X向、Y向和Z向的合力为Fix、Fiy和Fiz

Description

一种空间力系载荷谱编制方法
技术领域
本发明属于飞行器结构强度技术领域,涉及一种空间力系载荷谱编制方法。
背景技术
任何一种疲劳强度试验载荷谱的编制都是保证试验成功的关键环节。对于平行力系,尤其是飞机翼面载荷谱的编制,工程上已经给出了几种成熟的分区方法和明确的理论描述,但是对于对称试验件的空间力系,如何快速其疲劳载荷谱目前还没有跟好的方法。
本发明提出的空间力系载荷谱编制方法,由于是首次使用,所以进行了数值仿真,数值仿真结果表明,这种方法是可行的,根据此方法编制的疲劳载荷谱已经成功完成了试验。
发明内容
本发明的目的是:一种空间力系载荷谱编制方法。
1.一种空间力系载荷谱编制方法,其特征在于,包括以下步骤:
步骤一、载荷分区简化
令Pix、Piy和Piz依次为节点i在X向、Y向和Z向上的气动力,根据提供的分布节点载荷分别计算出n种载荷工况下X向载荷、Y向载荷和Z向载荷的总合力及相应压心坐标,对于n种载荷工况,m个气动载荷节点在X向、Y向和Z向的合力为Fix、Fiy和Fiz,其中
Fix的压心坐标(Yi,Zi):
Fiy的压心坐标(Xi,Zi):
Fiz的压心坐标(Xi,Yi):
坐标(Yi,Zi)为Y-Z平面X向载荷总压心,过(Yi,Zi)做X轴的平行线,将试验件分成两个区,每个区上载荷分别为Fxi1(y,z)和Fxi2(y,z),通过调整载荷Fxi1(y,z)和Fxi2(y,z)可以保证X向的总合力通过X向载荷总压心;
坐标(Xi,Zi)为X-Z平面Y向载荷总压心,过(Xi,Zi)做Z轴的平行线,将试验件分成两个区,每个区上载荷分别为Fyi1(x,z)和Fyi2(x,z),通过调整载荷Fyi1(x,z)和Fyi2(x,z)可以保证Y向的总合力通过Y向载荷总压心;
对Z向载荷的处理是在X-Y投影面进行分区,分区方法与Y向载荷不同,不同之处在于分区线只通过压心(Xi,Yi)但方向不受限制,每个区上载荷分别为Fzi1(x,y)和Fzi2(x,y),可以通过调整分区线的方向改变Fzi1(x,y)和Fzi2(x,y)的大小,将Z向载荷作用点的X坐标和Y坐标调整到压心误差范围内;
对Fxi1(y,z)、Fxi2(y,z)、Fyi1(x,z)、Fyi2(x,z)、Fzi1(x,y)、Fzi2(x,y)所在区域进行进一步分区,载荷分区依次进行,最终得到胶布带的粘贴位置和各个胶布带的载荷值,同时,根据载荷的大小和作用点位置可以计算出加载杠杆的长度和比例;
步骤二、在有限元模型上模拟胶布带粘贴
用计算机先在试验件的有限元模型上进行了胶布带粘贴,计算机粘贴胶布带的结果如果发生重叠则对载荷重新分区简化,直到不发生重叠为止;
步骤三、计算机模拟加载
用计算机在试验件的有限元模型上沿每级载荷的作用线进行了钢索布置,以验证试验安装中连接件在空间是否会发生干涉,并依此确定载荷是否有必要重新进行分区简化;
步骤四、计算验证
计算验证可以给出试验件上的最大法向应力,在计算中根据法向应力与试验件面板的剥离强度和胶布带承载能力判断疲劳试验过程中试验件面板是否会发生局部剥离、胶布带是否会被拉脱。
附图说明
图1是Y向载荷压心调整原理图;
图2是Z向载荷压心调整原理图;
图3是载荷计算的程序流程图;
图4是胶布带粘贴示意图;
图5是计算模型上布置加载钢索示意图;
图6是结构有限元模型;
图7是胶布带区域加密结构有限元模型;
图8是试验载荷应变云图。
具体实施方式
发明原理:
本发明应用计算机软件对试验件承受的分布气动载荷在三个投影平面(X-Z、Y-Z、X-Y)上分别对对应的轴向载荷用交互分区法进行载荷简化,再按照载荷简化得到的胶布带位置用计算机在试验件的有限元模型上进行粘贴和计算机模拟加载,使得载荷的简化最终达到既能满足载荷精度要求、又能满足三个载荷方向的胶布带不重叠、加载连接件在空间不相互干涉的要求为止。按照胶布带粘贴位置施加局部压力载荷可以确定试验件的危险部位,为试验中应变片的布置提供依据,同时计算机模拟加载还可以给出试验件上的最大法向应力,在计算中能够根据法向应力与试验件面板的剥离强度和胶布带承载能力判断疲劳试验过程中试验件面板是否会发生局部剥离、胶布带是否会被拉脱。下面结合附图和实施例对本发明做进一步说明:
步骤一、载荷分区简化
令Pix、Piy和Piz依次为节点i在X向、Y向和Z向上的气动力,根据提供的分布节点载荷分别计算出n种载荷工况下X向载荷、Y向载荷和Z向载荷的总合力及相应压心坐标,对于n种载荷工况,m个气动载荷节点在X向、Y向和Z向的合力为Fix、Fiy和Fiz,其中
Fix的压心坐标(Yi,Zi):
Fiy的压心坐标(Xi,Zi):
Fiz的压心坐标(Xi,Yi):
见图1,坐标(Yi,Zi)为Y-Z平面X向载荷总压心,过(Yi,Zi)做X轴的平行线,将试验件分成两个区,每个区上载荷分别为Fxi1(y,z)和Fxi2(y,z),通过调整载荷Fxi1(y,z)和Fxi2(y,z)可以保证X向的总合力通过X向载荷总压心;
坐标(Xi,Zi)为X-Z平面Y向载荷总压心,过(Xi,Zi)做Z轴的平行线,将试验件分成两个区,每个区上载荷分别为Fyi1(x,z)和Fyi2(x,z),通过调整载荷Fyi1(x,z)和Fyi2(x,z)可以保证Y向的总合力通过Y向载荷总压心;
见图2,对Z向载荷的处理是在X-Y投影面进行分区,分区方法与Y向载荷不同,不同之处在于分区线只通过压心(Xi,Yi)但方向不受限制,每个区上载荷分别为Fzi1(x,y)和Fzi2(x,y),可以通过调整分区线的方向改变Fzi1(x,y)和Fzi2(x,y)的大小,将Z向载荷作用点的X坐标和Y坐标调整到压心误差范围内;
对Fxi1(y,z)、Fxi2(y,z)、Fyi1(x,z)、Fyi2(x,z)、Fzi1(x,y)、Fzi2(x,y)所在区域进行进一步分区,载荷分区依次进行,最终得到胶布带的粘贴位置和各个胶布带的载荷值,同时,根据载荷的大小和作用点位置可以计算出加载杠杆的长度和比例(图3)。
步骤二、在有限元模型上模拟胶布带粘贴
由于载荷的简化是分别在三个投影面上进行的,载荷简化计算给出的是三个方向上的胶布带粘贴坐标。考虑到胶布带实际粘贴,不同方向载荷的胶布带位置可能会局部重叠,用计算机先在试验件的有限元模型上进行了胶布带粘贴(图4)。计算机粘贴胶布带的结果如果发生重叠则对载荷重新分区简化,直到不发生重叠为止。
步骤三、计算机模拟加载
试验载荷分三个主轴方向施加,施加不同方向载荷的连接件(加载钢索、杠杆)要在试验件壳体内部布置,载荷简化完成后,用计算机在试验件的有限元模型上沿每级载荷的作用线进行了钢索布置,以验证试验安装中连接件在空间是否会发生干涉,依此确定载荷是否有必要重新进行分区简化、是否需要按照钢索布置情况设计专用的加载夹具(图5)。
步骤四、计算验证
计算验证的目的是验证试验载荷的正确性、确定试验件的危险部位,为试验中应变片的布置提供依据;同时计算验证可以给出试验件上的最大法向应力,在计算中根据法向应力与试验件面板的剥离强度和胶布带承载能力判断疲劳试验过程中试验件面板是否会发生局部剥离、胶布带是否会被拉脱。

Claims (1)

1.一种空间力系载荷谱编制方法,其特征在于,包括以下步骤:
步骤一、载荷分区简化
令Pix、Piy和Piz依次为节点i在X向、Y向和Z向上的气动力,根据提供的分布节点载荷分别计算出n种载荷工况下X向载荷、Y向载荷和Z向载荷的总合力及相应压心坐标,对于n种载荷工况,m个气动载荷节点在X向、Y向和Z向的合力为Fix、Fiy和Fiz,其中
F i x = Σ j = 1 m P i x , ( i = 1 , n )
F i y = Σ j = 1 m P i y , ( i = 1 , n )
F i z = Σ j = 1 m P i z , ( i = 1 , n )
Fix的压心坐标(Yi,Zi):
Y i = Σ j = 1 m ( P i j x Y j ) F i x , Z i = Σ j = 1 m ( P i j x Z j ) F i x , ( i = 1 , n )
Fiy的压心坐标(Xi,Zi):
X i = Σ j = 1 m ( P i j y X j ) F i y , Z i = Σ j = 1 m ( P i j y Z j ) F i y , ( i = 1 , n )
Fiz的压心坐标(Xi,Yi):
X i = Σ j = 1 m ( P i j z Y j ) F i z , Y i = Σ j = 1 m ( P i j z Y j ) F i z , ( i = 1 , n )
坐标(Yi,Zi)为Y-Z平面X向载荷总压心,过(Yi,Zi)做X轴的平行线,将试验件分成两个区,每个区上载荷分别为Fxi1(y,z)和Fxi2(y,z),通过调整载荷Fxi1(y,z)和Fxi2(y,z)可以保证X向的总合力通过X向载荷总压心;
坐标(Xi,Zi)为X-Z平面Y向载荷总压心,过(Xi,Zi)做Z轴的平行线,将试验件分成两个区,每个区上载荷分别为Fyi1(x,z)和Fyi2(x,z),通过调整载荷Fyi1(x,z)和Fyi2(x,z)可以保证Y向的总合力通过Y向载荷总压心;
对Z向载荷的处理是在X-Y投影面进行分区,分区方法与Y向载荷不同,不同之处在于分区线只通过压心(Xi,Yi)但方向不受限制,每个区上载荷分别为Fzi1(x,y)和Fzi2(x,y),可以通过调整分区线的方向改变Fzi1(x,y)和Fzi2(x,y)的大小,将Z向载荷作用点的X坐标和Y坐标调整到压心误差范围内;
对Fxi1(y,z)、Fxi2(y,z)、Fyi1(x,z)、Fyi2(x,z)、Fzi1(x,y)、Fzi2(x,y)所在区域进行进一步分区,载荷分区依次进行,最终得到胶布带的粘贴位置和各个胶布带的载荷值,同时,根据载荷的大小和作用点位置可以计算出加载杠杆的长度和比例;
步骤二、在有限元模型上模拟胶布带粘贴
用计算机先在试验件的有限元模型上进行了胶布带粘贴,计算机粘贴胶布带的结果如果发生重叠则对载荷重新分区简化,直到不发生重叠为止;
步骤三、计算机模拟加载
用计算机在试验件的有限元模型上沿每级载荷的作用线进行了钢索布置,以验证试验安装中连接件在空间是否会发生干涉,并依此确定载荷是否有必要重新进行分区简化;
步骤四、计算验证
计算验证可以给出试验件上的最大法向应力,在计算中根据法向应力与试验件面板的剥离强度和胶布带承载能力判断疲劳试验过程中试验件面板是否会发生局部剥离、胶布带是否会被拉脱。
CN201510633822.8A 2015-09-29 2015-09-29 一种空间力系载荷谱编制方法 Active CN106557599B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510633822.8A CN106557599B (zh) 2015-09-29 2015-09-29 一种空间力系载荷谱编制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510633822.8A CN106557599B (zh) 2015-09-29 2015-09-29 一种空间力系载荷谱编制方法

Publications (2)

Publication Number Publication Date
CN106557599A true CN106557599A (zh) 2017-04-05
CN106557599B CN106557599B (zh) 2019-11-15

Family

ID=58414681

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510633822.8A Active CN106557599B (zh) 2015-09-29 2015-09-29 一种空间力系载荷谱编制方法

Country Status (1)

Country Link
CN (1) CN106557599B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107813961A (zh) * 2017-10-25 2018-03-20 北京富力通达科技有限公司 多通道协调加载系统飞续飞谱实验方法
CN109214131A (zh) * 2018-10-30 2019-01-15 中国运载火箭技术研究院 一种误差优化的静力试验载荷设计方法及系统
CN111539068A (zh) * 2020-04-22 2020-08-14 中国飞机强度研究所 一种飞机结构静强度试验数据对比显示方法
CN112763304A (zh) * 2020-12-23 2021-05-07 北京机电工程研究所 一种疲劳试验加载谱形成方法和装置及疲劳性能测试方法
CN116698471A (zh) * 2023-08-07 2023-09-05 四川腾盾科技有限公司 一种飞行器舵面静强度试验方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5930155A (en) * 1997-03-25 1999-07-27 Hitachi Metals, Ltd. Method of evaluating endurance of vehicle wheel by computer simulation
CN104075868A (zh) * 2014-05-30 2014-10-01 西北工业大学 用于飞机襟、缝翼系统可靠性试验的气动载荷加载方法
CN104374570A (zh) * 2014-11-24 2015-02-25 中国航空动力机械研究所 直升机传动系统部件寿命的获取方法
CN104778372A (zh) * 2015-04-27 2015-07-15 中国直升机设计研究所 一种直升机蜂窝夹层许用缺陷确定方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5930155A (en) * 1997-03-25 1999-07-27 Hitachi Metals, Ltd. Method of evaluating endurance of vehicle wheel by computer simulation
CN104075868A (zh) * 2014-05-30 2014-10-01 西北工业大学 用于飞机襟、缝翼系统可靠性试验的气动载荷加载方法
CN104374570A (zh) * 2014-11-24 2015-02-25 中国航空动力机械研究所 直升机传动系统部件寿命的获取方法
CN104778372A (zh) * 2015-04-27 2015-07-15 中国直升机设计研究所 一种直升机蜂窝夹层许用缺陷确定方法

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107813961A (zh) * 2017-10-25 2018-03-20 北京富力通达科技有限公司 多通道协调加载系统飞续飞谱实验方法
CN107813961B (zh) * 2017-10-25 2020-09-01 北京富力通达科技有限公司 多通道协调加载系统飞续飞谱实验方法
CN109214131A (zh) * 2018-10-30 2019-01-15 中国运载火箭技术研究院 一种误差优化的静力试验载荷设计方法及系统
CN109214131B (zh) * 2018-10-30 2022-12-27 中国运载火箭技术研究院 一种误差优化的静力试验载荷设计方法及系统
CN111539068A (zh) * 2020-04-22 2020-08-14 中国飞机强度研究所 一种飞机结构静强度试验数据对比显示方法
CN112763304A (zh) * 2020-12-23 2021-05-07 北京机电工程研究所 一种疲劳试验加载谱形成方法和装置及疲劳性能测试方法
CN112763304B (zh) * 2020-12-23 2023-11-07 北京机电工程研究所 一种疲劳试验加载谱形成方法和装置及疲劳性能测试方法
CN116698471A (zh) * 2023-08-07 2023-09-05 四川腾盾科技有限公司 一种飞行器舵面静强度试验方法
CN116698471B (zh) * 2023-08-07 2023-11-07 四川腾盾科技有限公司 一种飞行器舵面静强度试验方法

Also Published As

Publication number Publication date
CN106557599B (zh) 2019-11-15

Similar Documents

Publication Publication Date Title
CN106557599A (zh) 一种空间力系载荷谱编制方法
EP3318993A1 (en) Truss antenna reflector deployment dynamics modelling method based on multi-body analysis test
EP2191928A1 (en) Method of, computer readable medium and system for determination of spot weld failure using finite element analysis with stress and location sensitivity scale factors
US9770873B2 (en) System and method for optimizing composite laminate structures
CN104075868A (zh) 用于飞机襟、缝翼系统可靠性试验的气动载荷加载方法
CN109573098B (zh) 一种全尺寸疲劳试验高垂尾垂向载荷加载设计方法
CN105808829A (zh) 一种基于cpu+gpu异构并行计算的透平机械叶片固有频率特性分析方法
CN111209699A (zh) 一种基于车辆悬架静载荷的自动化提取方法及系统
CN109299579B (zh) 大展弦比飞机风洞测力试验数据修正方法
CN106547963B (zh) 一种风机偏航制动盘极限强度的校核方法及装置
CN106599405B (zh) 一种计算主减速器与机身任意连接点位置实际载荷的方法
CN104029826B (zh) 一种飞机结构强度试验扣重计算方法
US9026406B2 (en) Method, apparatus and computer program product for determining the strain induced at a selected point in a stiffened panel structure in response to a load, taking into account one or more out of plane (OOP) effects
CN109490114A (zh) 一种全尺寸疲劳试验襟翼载荷加载方法
CN109726435B (zh) 一种接头钉群载荷计算的方法
CN111044302B (zh) 基于振动试验耦合系统的夹具有效性验证优化方法
CN114818410A (zh) 一种对接实体焊缝疲劳强度的仿真方法
CN116698471B (zh) 一种飞行器舵面静强度试验方法
CA2829683C (en) Method for high fidelity modeling of an aircraft electrical power system
CN110321571A (zh) 一种蜂窝板壳结构的力学参数数值提取方法
CN109592074B (zh) 一种飞机结构疲劳试验加载频率的确定方法
CN109902372B (zh) 一种基于有限元分析的电池卷芯模拟方法
CN114770977B (zh) 一种自动铺丝工装的设计方法、装置、设备及存储介质
CN103366094A (zh) 一种基于力载荷和力边界条件的确定复合材料结构强度的方法
CN115408909A (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