CN110210093B - 一种撞击载荷工程模型可靠性的验证方法及装置 - Google Patents
一种撞击载荷工程模型可靠性的验证方法及装置 Download PDFInfo
- Publication number
- CN110210093B CN110210093B CN201910433630.0A CN201910433630A CN110210093B CN 110210093 B CN110210093 B CN 110210093B CN 201910433630 A CN201910433630 A CN 201910433630A CN 110210093 B CN110210093 B CN 110210093B
- Authority
- CN
- China
- Prior art keywords
- model
- impact load
- course curve
- time
- load
- 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
Links
Classifications
-
- 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
-
- 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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明实施例涉及一种撞击载荷工程模型可靠性的验证方法及装置。所述方法包括:获取预先建立的大型商用飞机的几何模型;获取预先建立的大型商用飞机和刚性墙的有限元计算模型;利用所述有限元计算模型,计算所述大型商用飞机的几何模型以预设撞击速度撞击所述刚性墙的第一数值;利用经过修正的Riera模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第二数值;利用撞击载荷工程模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第三数值;将所述第三数值与第一数值、第二数值分别进行比较,以验证撞击载荷工程模型的可靠性。
Description
技术领域
本发明实施例涉及计算机技术领域,尤其涉及一种撞击载荷工程模型可靠性的验证方法及装置。
背景技术
目前大型商用飞机撞击核电站,导致的核电站整体结构响应分析方法主要为力时程分析方法以及飞射物目标相互作用分析方法,且上述方法由美国核管理委员会(NRC)认可,可以用于大型商用飞机撞击核电站的相关评估。
1968年Riera提出著名的冲击方程(也称为Riera模型),为研究大型商用飞机撞击核电站问题提出了理论指导,冲击方程如下:
F(t)=Pc[x(t)]+μ[x(t)]V(t)2
1977年Hornyik通过实验结果和上述冲击方程计算结果进行对比,对Riera模型进行修正得到如下公式:
F(t)=Pc[x(t)]+αμ[x(t)]V(t)2
然而经过修正的Riera模型在实际应用过程中面临如下问题:大型商用飞机结构复杂,体积较大,通过理论计算不能获得准确的PC值,而通过开展1:1的原型实验获取PC值代价太高。
发明内容
鉴于此,为解决上述技术问题或部分技术问题,本发明实施例提供一种撞击载荷工程模型可靠性的验证方法及装置。
第一方面,本发明实施例提供一种撞击载荷工程模型可靠性的验证方法,所述方法包括:
获取预先建立的大型商用飞机的几何模型;
获取预先建立的大型商用飞机和刚性墙的有限元计算模型;
利用所述有限元计算模型,计算所述大型商用飞机的几何模型以预设撞击速度撞击所述刚性墙的第一数值;
利用经过修正的Riera模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第二数值;
利用撞击载荷工程模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第三数值;
将所述第三数值与第一数值、第二数值分别进行比较,以验证撞击载荷工程模型的可靠性;
其中,所述撞击载荷工程模型可以通过以下方式得到:
确定压损载荷总冲量IP(t)与冲击载荷总冲量IF(t);
计算压损载荷总冲量IP(t)与冲击载荷总冲量IF(t)之比γ=IP/IF;
根据压损载荷总冲量IP(t)与冲击载荷总冲量IF(t)之比,将经过修正的Riera模型转换为撞击载荷工程模型。
在一个可能的实施方式中,所述利用所述有限元计算模型,计算所述大型商用飞机的几何模型以预设撞击速度撞击所述刚性墙的第一数值,包括:
对大型商用飞机的几何模型所采用的网格大小和网格加密模型的计算结果进行对比,研究有限元计算模型的网格尺寸效应;
在有限元计算模型的网格尺寸效应满足预设的精度要求之后,利用所述有限元计算模型,计算所述大型商用飞机的几何模型以预设撞击速度撞击所述刚性墙的第一数值。
在一个可能的实施方式中,所述利用撞击载荷工程模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第三数值,包括:
确定撞击载荷工程模型中参数;
在确定撞击载荷工程模型中参数之后,利用撞击载荷工程模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第三数值。
在一个可能的实施方式中,所述确定撞击载荷工程模型中参数,包括:
提取所述第一数值,获得刚性墙受到的撞击力载荷时程曲线,和机身平均速度时程曲线;
对机身平均速度时程曲线进行积分获得位移时程曲线,从而获得机身质量线密度时程曲线;
对机身质量线密度时程曲线进行积分获得已碰撞质量位移曲线,从而获得已碰撞质量时程曲线;
对机身平均速度求导获得机身平均加速度时程曲线,则任一时刻的压损载荷:
Pc(t)=a(t)[M总-M(t)]
分别对撞击力载荷时程曲线、压损载荷时程曲线、动载荷项时程曲线进行积分,获得撞击力载荷时程曲线对应的总冲量、压损载荷时程曲线对应的总冲量、动载荷项时程曲线对应的总冲量;
根据所述撞击力载荷时程曲线对应的总冲量、所述压损载荷时程曲线对应的总冲量、所述动载荷项时程曲线对应的总冲量,以及经过修正的Riera模型,确定撞击载荷工程模型中参数。
在一个可能的实施方式中,所述撞击载荷工程模型为:
第二方面,本发明实施例提供一种撞击载荷工程模型可靠性的验证装置,所述装置包括:
第一获取模块,用于获取预先建立的大型商用飞机的几何模型;
第二获取模块,用于获取预先建立的大型商用飞机和刚性墙的有限元计算模型;
第一计算模块,用于利用所述有限元计算模型,计算所述大型商用飞机的几何模型以预设撞击速度撞击所述刚性墙的第一数值;
第二计算模块,用于利用经过修正的Riera模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第二数值;
第三计算模块,用于利用撞击载荷工程模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第三数值;
验证模块,用于将所述第三数值与第一数值、第二数值分别进行比较,以验证撞击载荷工程模型的可靠性;
其中,所述撞击载荷工程模型可以通过以下方式得到:
确定压损载荷总冲量IP(t)与冲击载荷总冲量IF(t);
计算压损载荷总冲量IP(t)与冲击载荷总冲量IF(t)之比γ=IP/IF;
根据压损载荷总冲量IP(t)与冲击载荷总冲量IF(t)之比,将经过修正的Riera模型转换为撞击载荷工程模型。
在一个可能的实施方式中,所述第一计算模块具体用于:
对大型商用飞机的几何模型所采用的网格大小和网格加密模型的计算结果进行对比,研究有限元计算模型的网格尺寸效应;
在有限元计算模型的网格尺寸效应满足预设的精度要求之后,利用所述有限元计算模型,计算所述大型商用飞机的几何模型以预设撞击速度撞击所述刚性墙的第一数值。
在一个可能的实施方式中,所述第三数值计算模块包括:
参数确定子模块,用于确定撞击载荷工程模型中参数;
第三数值计算子模块,用于在确定撞击载荷工程模型中参数之后,利用撞击载荷工程模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第三数值。
在一个可能的实施方式中,所述参数确定子模块具体用于:
提取所述第一数值,获得刚性墙受到的撞击力载荷时程曲线,和机身平均速度时程曲线;
对机身平均速度时程曲线进行积分获得位移时程曲线,从而获得机身质量线密度时程曲线;
对机身质量线密度时程曲线进行积分获得已碰撞质量位移曲线,从而获得已碰撞质量时程曲线;
对机身平均速度求导获得机身平均加速度时程曲线,则任一时刻的压损载荷:
Pc(t)=a(t)[M总-M(t)]
分别对撞击力载荷时程曲线、压损载荷时程曲线、动载荷项时程曲线进行积分,获得撞击力载荷时程曲线对应的总冲量、压损载荷时程曲线对应的总冲量、动载荷项时程曲线对应的总冲量;
根据所述撞击力载荷时程曲线对应的总冲量、所述压损载荷时程曲线对应的总冲量、所述动载荷项时程曲线对应的总冲量,以及经过修正的Riera模型,确定撞击载荷工程模型中参数。
在一个可能的实施方式中,所述撞击载荷工程模型为:
本发明实施例提供的技术方案,通过撞击载荷工程模型,解决了不能获得准确的PC值的问题,可以在不知飞机压损载荷的情况下计算大型商用飞机以任意速度撞击核电站时,核电站受到的冲击载荷随时间的变化,并且根据有限元计算模型以及经过修正的Riera模型对撞击载荷工程模型的可靠性进行验证。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明实施例中记载的一些实施例,对于本领域普通技术人员来讲,还可以根据这些附图获得其他的附图
图1是本发明实施例的撞击载荷工程模型可靠性的验证方法的实施流程示意图;
图2是本发明实施例的有限元计算模型示意图;
图3是本发明实施例的参数α和γ与飞机初始速度的关系式的计算示意图;
图4是本发明实施例的飞机以300m/s速度撞击刚性墙,对比三种模型获得的冲击载荷时程曲线;
图5是本发明实施例的撞击载荷工程模型可靠性的验证装置的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为便于对本发明实施例的理解,下面将结合附图以具体实施例做进一步的解释说明,实施例并不构成对本发明实施例的限定。
在本发明实施例中提成一种撞击载荷工程模型,其中撞击载荷工程模型为:
所述撞击载荷工程模型可以通过以下方式得到:确定压损载荷总冲量IP(t)与冲击载荷总冲量IF(t);计算压损载荷总冲量IP(t)与冲击载荷总冲量IF(t)之比γ=IP/IF;根据压损载荷总冲量IP(t)与冲击载荷总冲量IF(t)之比,将经过修正的Riera模型转换为撞击载荷工程模型。由于压损载荷和冲击载荷作用的时间一致,假设撞击过程中压损载荷总冲量IP(t)与冲击载荷总冲量IF(t)之比γ=IP/IF,那么经过修正的Riera模型可以转化为撞击载荷工程模型。
对于上述撞击载荷工程模型,本发明实施例中可以通过如下方法验证其可靠性:
如图1所示,为本发明实施例所提供的一种撞击载荷工程模型可靠性的验证方法的实施流程示意图,该方法具体可以包括以下步骤:
S101,获取预先建立的大型商用飞机的几何模型;
在本发明实施例中,可以预先建立大型商用飞机的几何模型,其外形与实际相符,考虑了飞机的主要结构,其结构尺寸按照真实情况或类似型号选取。除了机身蒙皮,还考虑了内部主要结构如机身框、机身桁条、地板梁、机翼翼梁、机翼翼肋、机翼桁条、尾翼翼梁、尾翼翼肋、尾翼桁条以及发动机。
例如预先使用Solidworks建立波音737-800商用飞机的几何模型,获取与先建立的波音737-800商用飞机的几何模型。
S102,获取预先建立的大型商用飞机和刚性墙的有限元计算模型;
预先建立大型商用飞机和刚性墙的有限元计算模型,例如使用Hypermesh建立商用飞机和刚性墙的有限元计算模型,如图2所示,飞机有限元模型仅包含壳单元和梁单元,壳单元采用*MAT_JOHNSON_COOK(MAT015)作为材料本构模型并且采用Belytschko-Tsay单元算法,使用沙漏控制可进行单点积分;梁单元采用*MAT_PLASTIC_KINEMATIC(MAT003)作为材料本构模型并且采用Hughes-Liu截面积分单元算法,壳单元和梁单元采用共节点建模,发动机使用钢材,其余部位使用铝合金材料,材料参数如表一、表二、表三所示:
表1钢材及铝合金的材料属性
A | B | C | m | n | D1 | D2 | D3 | D4 | D5 | |
铝 | 369 | 684 | 0.0083 | 1.7 | 0.73 | 0.13 | 0.13 | -1.5 | 0.011 | 0 |
钢 | 350 | 275 | 0.022 | 1.0 | 0.36 | 0.05 | 3.44 | -2.12 | 0.002 | 0.61 |
表2 MAT015参数
σ0/MPa | Et/MPa | β | C | p | Fs | |
铝 | 503 | 5000 | 0 | 40 | 5 | 0.05 |
表3 MAT003参数
S103,利用所述有限元计算模型,计算所述大型商用飞机的几何模型以预设撞击速度撞击所述刚性墙的第一数值;
利用所述有限元计算模型,计算所述大型商用飞机的几何模型以预设撞击速度撞击所述刚性墙的第一数值,其中在计算第一数值之前,对大型商用飞机的几何模型所采用的网格大小和网格加密模型的计算结果进行对比,研究有限元计算模型的网格尺寸效应,在有限元计算模型的网格尺寸效应满足预设的精度要求之后,利用所述有限元计算模型,计算所述大型商用飞机的几何模型以预设撞击速度撞击所述刚性墙的第一数值。
例如,对飞机模型所采用的网格大小和网格加密50%模型的计算结果进行对比研究数值模型的网格尺寸效应,在满足精度要求之后开展不同速度下飞机撞击工况的数值模拟计算,飞机分别选取200m/s,250m/s,300m/s三个速度撞击刚性墙。
S104,利用经过修正的Riera模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第二数值;
利用经过修正的Riera模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第二数值,例如开展不同速度下飞机撞击工况的Riera模型数值计算,飞机分别选取200m/s,250m/s,300m/s三个速度撞击刚性墙,计算三个Riera模型数值。
S105,利用撞击载荷工程模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第三数值;
结合第一数值,即数值模拟结果,提出撞击载荷工程模型相关参数的确定方法:提取所述第一数值,获得刚性墙受到的撞击力载荷时程曲线,和机身平均速度时程曲线;对机身平均速度时程曲线进行积分获得位移时程曲线,从而获得机身质量线密度时程曲线;对机身质量线密度时程曲线进行积分获得已碰撞质量位移曲线,从而获得已碰撞质量时程曲线;对机身平均速度求导获得机身平均加速度时程曲线,则任一时刻的压损载荷:
Pc(t)=a(t)[M总-M(t)]
分别对撞击力载荷时程曲线、压损载荷时程曲线、动载荷项时程曲线进行积分,获得撞击力载荷时程曲线对应的总冲量、压损载荷时程曲线对应的总冲量、动载荷项时程曲线对应的总冲量;根据所述撞击力载荷时程曲线对应的总冲量、所述压损载荷时程曲线对应的总冲量、所述动载荷项时程曲线对应的总冲量,以及经过修正的Riera模型,确定撞击载荷工程模型中参数。
例如,提取数值计算结果获得刚性墙受到的撞击力载荷时程曲线F-t和机身平均速度时程曲线v-t,对机身平均速度时程曲线积分获得位移时程曲线从而获得机身质量线密度时程曲线μ-t(已知线密度位移曲线),对线密度位移曲线积分获得已碰撞质量位移曲线从而获得已碰撞质量时程曲线M-t,对机身平均速度求导获得机身平均加速度时程曲线a-t,则每一时刻的压损载荷:
Pc(t)=a(t)[M总-M(t)]
分别对撞击力载荷时程曲线F-t,静载荷项(压损载荷)时程曲线Pc-t,动载荷项时程曲线μv2-t进行积分,获得三者的总冲量IF,IP,IP',代入修正后的Riera理论模型公式积分形式求得参数α值,另一个参数γ=IP/IF。
在确定撞击载荷工程模型中参数之后,利用撞击载荷工程模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第三数值,例如开展不同速度下飞机撞击工况的撞击载荷工程模型数值计算,飞机分别选取200m/s,250m/s,300m/s三个速度撞击刚性墙,计算三个撞击载荷工程模型数值。
S106,将所述第三数值与第一数值、第二数值分别进行比较,以验证撞击载荷工程模型的可靠性。
将所述第三数值与第一数值、第二数值分别进行对比分析,验证载荷工程模型的可靠性。
例如飞机以275m/s速度撞击刚性墙时,通过相关参数的确定方法确定得到两个参数α和γ的值分别为0.87和0.106,而将V=275m/s代入图3所示参数α和γ与飞机初始速度的关系式计算得到α和γ的值分别为0.868和0.101,两个参数的误差分别为0.2%和4.7%;数值模拟和撞击载荷工程模型计算得到的撞击力峰值分别为359MN和375MN,误差为4.5%,两者的撞击力冲量分别为17.00MN·s和17.14MN·s,误差为0.8%。
又例如,例如飞机以150m/s速度撞击刚性墙时,通过相关参数的确定方法确定得到两个参数α和γ的值分别为0.6和0.357,而将V=150m/s代入图3所示参数α和γ与飞机初始速度的关系式计算得到α和γ的值分别为0.597和0.343,两个参数的误差分别为0.5%和3.9%;数值模拟和工程模型计算得到的撞击力峰值分别为76.3MN和81.7MN,误差为7.1%,两者的撞击力冲量分别为7.44MN·s和7.41MN·s,误差为-0.4%。
又例如,如图4所示,通过撞击载荷工程模型计算获得的冲击载荷时程曲线与数值模拟/Riera理论模型计算获得的冲击载荷时程曲线相比,无论包络线还是载荷峰值的偏差都比较小,表明我们提出的撞击载荷工程模型相对可靠,可以用于实际工程计算中。
通过上述对本发明实施例提供的技术方案的描述,通过撞击载荷工程模型,解决了不能获得准确的PC值的问题,可以在不知飞机压损载荷的情况下计算大型商用飞机以任意速度撞击核电站时,核电站受到的冲击载荷随时间的变化,并且根据有限元计算模型以及经过修正的Riera模型对撞击载荷工程模型的可靠性进行验证。
相对于上述方法实施例,本发明实施例还提供了一种撞击载荷工程模型可靠性的验证装置,如图5所示,可以包括:第一获取模块510、第二获取模块520、第一计算模块530、第二计算模块540、第三计算模块550、验证模块560。
第一获取模块510,用于获取预先建立的大型商用飞机的几何模型;
第二获取模块520,用于获取预先建立的大型商用飞机和刚性墙的有限元计算模型;
第一计算模块530,用于利用所述有限元计算模型,计算所述大型商用飞机的几何模型以预设撞击速度撞击所述刚性墙的第一数值;
第二计算模块540,用于利用经过修正的Riera模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第二数值;
第三计算模块550,用于利用撞击载荷工程模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第三数值;
验证模块560,用于将所述第三数值与第一数值、第二数值分别进行比较,以验证撞击载荷工程模型的可靠性;
其中,所述撞击载荷工程模型可以通过以下方式得到:
确定压损载荷总冲量IP(t)与冲击载荷总冲量IF(t);
计算压损载荷总冲量IP(t)与冲击载荷总冲量IF(t)之比γ=IP/IF;
根据压损载荷总冲量IP(t)与冲击载荷总冲量IF(t)之比,将经过修正的Riera模型转换为撞击载荷工程模型。
根据本申请提供的一种具体实施方式,所述第一计算模块530具体用于:
对大型商用飞机的几何模型所采用的网格大小和网格加密模型的计算结果进行对比,研究有限元计算模型的网格尺寸效应;
在有限元计算模型的网格尺寸效应满足预设的精度要求之后,利用所述有限元计算模型,计算所述大型商用飞机的几何模型以预设撞击速度撞击所述刚性墙的第一数值。
根据本申请提供的一种具体实施方式,所述第三数值计算模块550包括:
参数确定子模块551,用于确定撞击载荷工程模型中参数;
第三数值计算子模块552,用于在确定撞击载荷工程模型中参数之后,利用撞击载荷工程模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第三数值。
根据本申请提供的一种具体实施方式,所述参数确定子模块551具体用于:
提取所述第一数值,获得刚性墙受到的撞击力载荷时程曲线,和机身平均速度时程曲线;
对机身平均速度时程曲线进行积分获得位移时程曲线,从而获得机身质量线密度时程曲线;
对机身质量线密度时程曲线进行积分获得已碰撞质量位移曲线,从而获得已碰撞质量时程曲线;
对机身平均速度求导获得机身平均加速度时程曲线,则任一时刻的压损载荷:
Pc(t)=a(t)[M总-M(t)]
分别对撞击力载荷时程曲线、压损载荷时程曲线、动载荷项时程曲线进行积分,获得撞击力载荷时程曲线对应的总冲量、压损载荷时程曲线对应的总冲量、动载荷项时程曲线对应的总冲量;
根据所述撞击力载荷时程曲线对应的总冲量、所述压损载荷时程曲线对应的总冲量、所述动载荷项时程曲线对应的总冲量,以及经过修正的Riera模型,确定撞击载荷工程模型中参数。
根据本申请提供的一种具体实施方式,所述撞击载荷工程模型为:
上述装置中各个模块的功能和作用的实现过程具体详见上述方法中对应步骤的实现过程,在此不再赘述。
通过上述对本发明实施例提供的技术方案的描述,通过撞击载荷工程模型,解决了不能获得准确的PC值的问题,可以在不知飞机压损载荷的情况下计算大型商用飞机以任意速度撞击核电站时,核电站受到的冲击载荷随时间的变化,并且根据有限元计算模型以及经过修正的Riera模型对撞击载荷工程模型的可靠性进行验证。
专业人员应该还可以进一步意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、计算机软件或者二者的结合来实现,为了清楚地说明硬件和软件的可互换性,在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
结合本文中所公开的实施例描述的方法或算法的步骤可以用硬件、处理器执行的软件模块,或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种撞击载荷工程模型可靠性的验证方法,其特征在于,所述方法包括:
获取预先建立的大型商用飞机的几何模型;
获取预先建立的大型商用飞机和刚性墙的有限元计算模型;
利用所述有限元计算模型,计算所述大型商用飞机的几何模型以预设撞击速度撞击所述刚性墙的第一数值;
利用经过修正的Riera模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第二数值;
利用撞击载荷工程模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第三数值;
将所述第三数值与第一数值、第二数值分别进行比较,以验证撞击载荷工程模型的可靠性;
其中,所述撞击载荷工程模型可以通过以下方式得到:
确定压损载荷总冲量IP(t)与冲击载荷总冲量IF(t);
计算压损载荷总冲量IP(t)与冲击载荷总冲量IF(t)之比γ=IP/IF;
根据压损载荷总冲量IP(t)与冲击载荷总冲量IF(t)之比,将经过修正的Riera模型转换为撞击载荷工程模型;
所述利用所述有限元计算模型,计算所述大型商用飞机的几何模型以预设撞击速度撞击所述刚性墙的第一数值,包括:
对大型商用飞机的几何模型所采用的网格大小和网格加密模型的计算结果进行对比,研究有限元计算模型的网格尺寸效应;
在有限元计算模型的网格尺寸效应满足预设的精度要求之后,利用所述有限元计算模型,计算所述大型商用飞机的几何模型以预设撞击速度撞击所述刚性墙的第一数值;
所述利用撞击载荷工程模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第三数值,包括:
确定撞击载荷工程模型中参数;
在确定撞击载荷工程模型中参数之后,利用撞击载荷工程模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第三数值;
所述确定撞击载荷工程模型中参数,包括:
提取所述第一数值,获得刚性墙受到的撞击力载荷时程曲线,和机身平均速度时程曲线;
对机身平均速度时程曲线进行积分获得位移时程曲线,从而获得机身质量线密度时程曲线;
对机身质量线密度时程曲线进行积分获得已碰撞质量位移曲线,从而获得已碰撞质量时程曲线;
对机身平均速度求导获得机身平均加速度时程曲线,则任一时刻的压损载荷:
Pc(t)=a(t)[M总-M(t)]
分别对撞击力载荷时程曲线、压损载荷时程曲线、动载荷项时程曲线进行积分,获得撞击力载荷时程曲线对应的总冲量、压损载荷时程曲线对应的总冲量、动载荷项时程曲线对应的总冲量;
根据所述撞击力载荷时程曲线对应的总冲量、所述压损载荷时程曲线对应的总冲量、所述动载荷项时程曲线对应的总冲量,以及经过修正的Riera模型,确定撞击载荷工程模型中参数;
所述撞击载荷工程模型为:
2.一种撞击载荷工程模型可靠性的验证装置,其特征在于,所述装置包括:
第一获取模块,用于获取预先建立的大型商用飞机的几何模型;
第二获取模块,用于获取预先建立的大型商用飞机和刚性墙的有限元计算模型;
第一计算模块,用于利用所述有限元计算模型,计算所述大型商用飞机的几何模型以预设撞击速度撞击所述刚性墙的第一数值;
第二计算模块,用于利用经过修正的Riera模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第二数值;
第三计算模块,用于利用撞击载荷工程模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第三数值;
验证模块,用于将所述第三数值与第一数值、第二数值分别进行比较,以验证撞击载荷工程模型的可靠性;
其中,所述撞击载荷工程模型可以通过以下方式得到:
确定压损载荷总冲量IP(t)与冲击载荷总冲量IF(t);
计算压损载荷总冲量IP(t)与冲击载荷总冲量IF(t)之比γ=IP/IF;
根据压损载荷总冲量IP(t)与冲击载荷总冲量IF(t)之比,将经过修正的Riera模型转换为撞击载荷工程模型;
所述第一计算模块具体用于:
对大型商用飞机的几何模型所采用的网格大小和网格加密模型的计算结果进行对比,研究有限元计算模型的网格尺寸效应;
在有限元计算模型的网格尺寸效应满足预设的精度要求之后,利用所述有限元计算模型,计算所述大型商用飞机的几何模型以预设撞击速度撞击所述刚性墙的第一数值;
所述第三数值计算模块包括:
参数确定子模块,用于确定撞击载荷工程模型中参数;
第三数值计算子模块,用于在确定撞击载荷工程模型中参数之后,利用撞击载荷工程模型,计算所述大型商用飞机的几何模型以所述预设撞击速度撞击所述刚性墙的第三数值;
所述参数确定子模块具体用于:
提取所述第一数值,获得刚性墙受到的撞击力载荷时程曲线,和机身平均速度时程曲线;
对机身平均速度时程曲线进行积分获得位移时程曲线,从而获得机身质量线密度时程曲线;
对机身质量线密度时程曲线进行积分获得已碰撞质量位移曲线,从而获得已碰撞质量时程曲线;
对机身平均速度求导获得机身平均加速度时程曲线,则任一时刻的压损载荷:
Pc(t)=a(t)[M总-M(t)]
分别对撞击力载荷时程曲线、压损载荷时程曲线、动载荷项时程曲线进行积分,获得撞击力载荷时程曲线对应的总冲量、压损载荷时程曲线对应的总冲量、动载荷项时程曲线对应的总冲量;
根据所述撞击力载荷时程曲线对应的总冲量、所述压损载荷时程曲线对应的总冲量、所述动载荷项时程曲线对应的总冲量,以及经过修正的Riera模型,确定撞击载荷工程模型中参数;
所述撞击载荷工程模型为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910433630.0A CN110210093B (zh) | 2019-05-23 | 2019-05-23 | 一种撞击载荷工程模型可靠性的验证方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910433630.0A CN110210093B (zh) | 2019-05-23 | 2019-05-23 | 一种撞击载荷工程模型可靠性的验证方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110210093A CN110210093A (zh) | 2019-09-06 |
CN110210093B true CN110210093B (zh) | 2023-09-26 |
Family
ID=67788353
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910433630.0A Active CN110210093B (zh) | 2019-05-23 | 2019-05-23 | 一种撞击载荷工程模型可靠性的验证方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110210093B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110610048B (zh) * | 2019-09-18 | 2023-02-28 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种针对偏心撞击试验的动载荷因数计算方法 |
CN114004032B (zh) * | 2021-11-11 | 2022-07-22 | 南京理工大学 | 一种抗冲击载荷的板梁组合结构优化设计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2372598C1 (ru) * | 2008-05-13 | 2009-11-10 | Открытое акционерное общество "Информационные спутниковые системы" имени академика М.Ф. Решетнева | Способ испытаний космического аппарата на ударные воздействия |
CN103995917A (zh) * | 2014-04-17 | 2014-08-20 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种全尺寸飞机着陆撞击的起落架载荷仿真方法 |
CN107742025A (zh) * | 2017-10-16 | 2018-02-27 | 北京航空航天大学 | 一种预测飞机金属蒙皮冲击后疲劳寿命的方法 |
-
2019
- 2019-05-23 CN CN201910433630.0A patent/CN110210093B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2372598C1 (ru) * | 2008-05-13 | 2009-11-10 | Открытое акционерное общество "Информационные спутниковые системы" имени академика М.Ф. Решетнева | Способ испытаний космического аппарата на ударные воздействия |
CN103995917A (zh) * | 2014-04-17 | 2014-08-20 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种全尺寸飞机着陆撞击的起落架载荷仿真方法 |
CN107742025A (zh) * | 2017-10-16 | 2018-02-27 | 北京航空航天大学 | 一种预测飞机金属蒙皮冲击后疲劳寿命的方法 |
Non-Patent Citations (2)
Title |
---|
飞机撞击下核电厂安全壳载荷曲线相关问题研究;温丽晶等;《安全与环境学报》;20151225(第06期);全文 * |
飞机模型高速撞击钢筋混凝土载荷特性实验研究;温丽晶等;《爆炸与冲击》;20180704(第04期);第811-818页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110210093A (zh) | 2019-09-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ghoreyshi et al. | Reduced order unsteady aerodynamic modeling for stability and control analysis using computational fluid dynamics | |
Reuther et al. | Aerodynamic shape optimization of wing and wing-body configurations using control theory | |
Murman | Reduced-frequency approach for calculating dynamic derivatives | |
Alkishriwi et al. | A large-eddy simulation method for low Mach number flows using preconditioning and multigrid | |
CN110210093B (zh) | 一种撞击载荷工程模型可靠性的验证方法及装置 | |
Currao et al. | Hypersonic fluid-structure interaction on a cantilevered plate | |
Housman et al. | Efficient near-field to mid-field sonic boom propagation using a high-order space marching method | |
US11188692B2 (en) | Turbulent boundary layer modeling via incorporation of pressure gradient directional effect | |
Larsen | Impact loads on circular cylinders | |
Vallabh et al. | Investigation of nozzle contours in the CSIR supersonic wind tunnel | |
Min et al. | Hybrid Navier-Stokes/free-wake method for modeling blade-vortex interactions | |
US20200285709A1 (en) | Turbulent Boundary Layer Modeling via Incorporation of Pressure Gradient Directional Effect | |
Chaudhry et al. | Numerical evaluation of the hydrodynamic impact characteristics of the air launched AUV upon water entry | |
Murman et al. | Dynamic analysis of atmospheric-entry probes and capsules | |
US11835054B2 (en) | Method for automatic detection of axial cooling fan rotation direction | |
Maxwell | Efficient Design of Viscous Waveriders with CFD Verification and Off-Design Performance Analysis | |
Roy et al. | Estimation of numerical error for 3D inviscid flows on Cartesian grids | |
Dadone et al. | Fast convergence of viscous airfoil design problems | |
Kim et al. | Nonlinearc aeroelastic analysis of an airfoil using CFD-based indicial approach | |
Mohseni et al. | Improvement of compressible vorticity confinement method by combining it with vortex feature detection methods | |
Bartels et al. | Computation of shock buffet onset for a conventional and supercritical airfoil | |
Farsadi et al. | Expansion of indicial function approximations for 2-D subsonic compressible aerodynamic loads | |
Anderson et al. | Accurate solutions, parameter studies and comparisons for the euler and potential flow equations | |
Otero et al. | Direct numerical simulations for adjoint-based optimal flow and noise control of a backward-facing step | |
Abdullah et al. | Measurement Correction of Two-Dimensional Wind Tunnel with Porous Walls Using Singularity Method |
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 |