CN113704888B - 一种单元应力筛选方法 - Google Patents

一种单元应力筛选方法 Download PDF

Info

Publication number
CN113704888B
CN113704888B CN202110970081.8A CN202110970081A CN113704888B CN 113704888 B CN113704888 B CN 113704888B CN 202110970081 A CN202110970081 A CN 202110970081A CN 113704888 B CN113704888 B CN 113704888B
Authority
CN
China
Prior art keywords
stress
value
unit stress
unit
area
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
CN202110970081.8A
Other languages
English (en)
Other versions
CN113704888A (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 CN202110970081.8A priority Critical patent/CN113704888B/zh
Publication of CN113704888A publication Critical patent/CN113704888A/zh
Application granted granted Critical
Publication of CN113704888B publication Critical patent/CN113704888B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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/14—Force analysis or force optimisation, e.g. static or dynamic forces

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)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本申请属于飞机强度设计领域，特别涉及一种单元应力筛选方法。包括：步骤一、获取预定区域内的单元应力数据；步骤二、对当前区域内的所有单元按照应力升序进行排序；步骤三、计算当前区域内的单元应力均值；步骤四、计算当前区域内的单元应力中位值；步骤五、计算当前区域内的所述单元应力均值与所述单元应力中位值的差值；步骤六、判断所述单元应力均值与所述单元应力中位值的差值是否小于阈值；若是，则取最大单元应力值为当前区域的应力；若否，则去掉最大单元应力值，单元总数减1，返回步骤三。本申请根据同一区域内的单元应力中位值和单元应力平均值的差值决定区域应力水平，能够快速给出该区域单元应力的合理水平值。

Description

一种单元应力筛选方法
技术领域
本申请属于飞机强度设计领域，特别涉及一种单元应力筛选方法。
背景技术
随着航空结构设计与制造水平的不断提高，飞行器设计性能指标也越来越高，结构重量系数作为影响飞行器性能的主要因素之一备受关注，而强度设计指标选取的是否合适直接决定了飞行器结构重量系数。结构单元的应力水平是结构设计最基础的强度指标，在强度校核和结构优化过程中广泛使用。同一个区域的单元应力值会分布在一定范围内，但为了控制实际参与优化的单元应力约束的数量，往往需要采用一种方法在同一个区域内选择有限个单元应力来代表该区域的整体应力水平，如果取最大应力值，导致设计太保守，而且偶尔出现的异常单元应力会使结果严重失真；如果取平均值，有可能漏掉某些集中的高应力区域，导致设计偏危险。
目前工程应用时大多直接取最大值以确保设计安全，但这种方式会直接导致结构超重，从而影响未来飞行器整体性能的提升。
因此，希望有一种技术方案来克服或至少减轻现有技术的至少一个上述缺陷。
发明内容
本申请的目的是提供了一种单元应力筛选方法，以解决现有技术存在的至少一个问题。
本申请的技术方案是：
一种单元应力筛选方法，包括：
步骤一、获取预定区域内的单元应力数据；
步骤二、对当前区域内的所有单元按照应力升序进行排序；
步骤三、计算当前区域内的单元应力均值；
步骤四、计算当前区域内的单元应力中位值；
步骤五、计算当前区域内的所述单元应力均值与所述单元应力中位值的差值；
步骤六、判断所述单元应力均值与所述单元应力中位值的差值是否小于阈值；
若是，则取最大单元应力值为当前区域的应力；
若否，则去掉最大单元应力值，单元总数减1，返回步骤三。
在本申请的一个实施方式中，步骤三中，所述计算当前区域内的单元应力均值包括：
其中，为单元应力均值，σ(i)为对应单元应力值，n为单元总数。
在本申请的一个实施方式中，步骤四中，所述计算当前区域内的单元应力中位值包括：
其中，σM为单元应力中位值，n为单元总数。
在本申请的一个实施方式中，步骤五中，所述计算当前区域内的所述单元应力均值与所述单元应力中位值的差值包括：
其中，Δσ为单元应力均值与单元应力中位值的差值。
在本申请的一个实施方式中，步骤六中，所述阈值为所述单元应力均值的10％。
发明至少存在以下有益技术效果：
本申请的单元应力筛选方法，根据同一区域内的单元应力中位值和单元应力平均值的差值决定区域应力水平，能够快速给出该区域单元应力的合理水平值，为强度校核和优化设计中的应力筛选提供支持。
附图说明
图1是本申请一个实施方式的单元应力筛选方法流程图。
具体实施方式
为使本申请实施的目的、技术方案和优点更加清楚，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行更加详细的描述。在附图中，自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本申请一部分实施例，而不是全部的实施例。下面通过参考附图描述的实施例是示例性的，旨在用于解释本申请，而不能理解为对本申请的限制。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。下面结合附图对本申请的实施例进行详细说明。
在本申请的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本申请和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本申请保护范围的限制。
下面结合附图1对本申请做进一步详细说明。
本申请提供了一种单元应力筛选方法，包括以下步骤：
S001、获取预定区域内的单元应力数据；
S002、对当前区域内的所有单元按照应力升序进行排序；
S003、计算当前区域内的单元应力均值；
S004、计算当前区域内的单元应力中位值；
S005、计算当前区域内的单元应力均值与单元应力中位值的差值；
S006、判断单元应力均值与单元应力中位值的差值是否小于阈值；
若是，则取最大单元应力值为当前区域的应力；
若否，则去掉最大单元应力值，单元总数减1，返回S003。直至筛选出当前区域的单元应力，结束当前区域单元应力筛选。
本申请的单元应力筛选方法，对于一个包含n个单元的区域，步骤S003中，计算当前区域内的单元应力均值包括：
其中，为单元应力均值，σ(i)为对应单元应力值，n为单元总数。
本申请的单元应力筛选方法，步骤S004中，计算当前区域内的单元应力中位值包括：
其中，σM为单元应力中位值，n为单元总数。
本申请的单元应力筛选方法，步骤S005中，计算当前区域内的单元应力均值与单元应力中位值的差值包括：
其中，Δσ为单元应力均值与单元应力中位值的差值。
在本申请的优选实施例中，步骤S006中，阈值设置为单元应力均值的10％。
在本申请的一个实施方式中，实现对两个区域内的单元应力筛选，每个区域内均包含6个有限元单元，单元中间标识了该单元的编号和该单元对应的应力值。区域1内的单元应力数据参见表1，区域2内的单元应力数据参见表2。
表1
区表域21
区本域1实施例中，对区域1内的单元区域应2力筛选过程如下：
对区域1内的6个单元按照应力值从小到大排序，结果为：3-123MPa、6-155MPa、2-168MPa、5-178MPa、1-232MPa、4-256MPa，则按照公式(1)计算得到区域1内的6个单元的单元应力均值为185.3MPa，按照公式(2)计算得到区域1内的6个单元的单元应力中位值σM为173MPa，按照公式(3)计算得到区域1内单元应力均值与单元应力中位值的差值Δσ为12.3MPa，判断该差值小于10％的单元应力平均值(18.53MPa)，则直接选择当前最大单元应力值256MPa作为区域1的应力，用于后续的优化设计和强度校核。
本实施例中，对区域2内的单元应力筛选过程如下：
对区域2内的6个单元按照应力值从小到大排序，结果为：10-131MPa、11-145MPa、7-157MPa、8-176MPa、9-234MPa、12-359MPa，则按照公式(1)计算得到区域2内的6个单元的单元应力均值σ为200.3MPa，按照公式(2)计算得到区域2内的6个单元的单元应力中位值σM为166.5MPa，按照公式(3)计算得到区域2内单元应力均值与单元应力中位值的差值Δσ为33.8MPa，判断该差值大于10％的单元应力平均值(20.03MPa)，则去掉最大单元应力值359MPa，单元总数n＝n-1，返回S003，按照公式(1)计算得到区域2内的5个单元的单元应力均值为168.6MPa，按照公式(2)计算得到区域2内的5个单元的单元应力中位值σM为157MPa，按照公式(3)计算得到区域2内单元应力均值与单元应力中位值的差值Δσ为11.8MPa，判断该差值小于10％的单元应力平均值(16.86MPa)，选择当前最大单元应力值234MPa作为区域2的应力，用于后续的优化设计和强度校核。
本申请的单元应力筛选方法，根据同一区域的单元应力中位值和单元应力平均值的差值，判断是否剔除最高单元应力，逐步缩小该区域单元应力中位值和单元应力平均值的差值，直至达到设计要求，进而选取此时剩余单元应力的最大值作为该区域的应力水平值，即减小了直接取最大应力值带来的重量冗余，也避免了直接取均值引入的潜在欠强度风险，能够快速给出该区域单元应力的合理水平值，为强度校核和优化设计中的应力筛选提供了支持。
以上所述，仅为本申请的具体实施方式，但本申请的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本申请揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本申请的保护范围之内。因此，本申请的保护范围应以所述权利要求的保护范围为准。

Claims (3)

1.一种单元应力筛选方法，其特征在于，包括：
步骤一、获取预定区域内的单元应力数据；
步骤二、对当前区域内的所有单元按照应力升序进行排序；
步骤三、计算当前区域内的单元应力均值；
步骤四、计算当前区域内的单元应力中位值；
步骤五、计算当前区域内的所述单元应力均值与所述单元应力中位值的差值；
步骤六、判断所述单元应力均值与所述单元应力中位值的差值是否小于阈值；
若是，则取最大单元应力值为当前区域的应力；
若否，则去掉最大单元应力值，单元总数减1，返回步骤三；
步骤三中，所述计算当前区域内的单元应力均值包括：
其中，为单元应力均值，σ(i)为对应单元应力值，n为单元总数；
步骤四中，所述计算当前区域内的单元应力中位值包括：
其中，σM为单元应力中位值，n为单元总数。
2.根据权利要求1所述的单元应力筛选方法，其特征在于，步骤五中，所述计算当前区域内的所述单元应力均值与所述单元应力中位值的差值包括：
其中，△σ为单元应力均值与单元应力中位值的差值。
3.根据权利要求1所述的单元应力筛选方法，其特征在于，步骤六中，所述阈值为所述单元应力均值的10％。
CN202110970081.8A 2021-08-23 2021-08-23 一种单元应力筛选方法 Active CN113704888B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110970081.8A CN113704888B (zh) 2021-08-23 2021-08-23 一种单元应力筛选方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110970081.8A CN113704888B (zh) 2021-08-23 2021-08-23 一种单元应力筛选方法

Publications (2)

Publication Number Publication Date
CN113704888A CN113704888A (zh) 2021-11-26
CN113704888B true CN113704888B (zh) 2024-02-23

Family

ID=78654196

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110970081.8A Active CN113704888B (zh) 2021-08-23 2021-08-23 一种单元应力筛选方法

Country Status (1)

Country Link
CN (1) CN113704888B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104517012A (zh) * 2014-12-25 2015-04-15 芜湖市汽车产业技术研究院有限公司 一种计算非线性结构的应变时间历程的方法和装置
CN107391867A (zh) * 2017-07-31 2017-11-24 吴锦 一种冲压零件的回弹补偿方法和装置
CN108009370A (zh) * 2017-12-13 2018-05-08 中国飞机强度研究所 一种结构应力敏度求解方法
CN108038335A (zh) * 2017-12-28 2018-05-15 北京数设科技有限公司 一种确定飞机蒙皮单元应力载荷的方法与装置
CN109359360A (zh) * 2018-09-30 2019-02-19 国家超级计算天津中心 一种基于局部特征的结构应力处理方法
CN112100774A (zh) * 2020-09-16 2020-12-18 哈尔滨理工大学 一种基于变密度法的应力和应变能双约束的拓扑优化方法
WO2021056977A1 (zh) * 2019-09-24 2021-04-01 重庆美的通用制冷设备有限公司 管路振动控制方法、计算机设备、存储介质及管路系统
CN113239584A (zh) * 2021-04-26 2021-08-10 云南大学 一种优化增材制造方法及系统

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104517012A (zh) * 2014-12-25 2015-04-15 芜湖市汽车产业技术研究院有限公司 一种计算非线性结构的应变时间历程的方法和装置
CN107391867A (zh) * 2017-07-31 2017-11-24 吴锦 一种冲压零件的回弹补偿方法和装置
CN108009370A (zh) * 2017-12-13 2018-05-08 中国飞机强度研究所 一种结构应力敏度求解方法
CN108038335A (zh) * 2017-12-28 2018-05-15 北京数设科技有限公司 一种确定飞机蒙皮单元应力载荷的方法与装置
CN109359360A (zh) * 2018-09-30 2019-02-19 国家超级计算天津中心 一种基于局部特征的结构应力处理方法
WO2021056977A1 (zh) * 2019-09-24 2021-04-01 重庆美的通用制冷设备有限公司 管路振动控制方法、计算机设备、存储介质及管路系统
CN112100774A (zh) * 2020-09-16 2020-12-18 哈尔滨理工大学 一种基于变密度法的应力和应变能双约束的拓扑优化方法
CN113239584A (zh) * 2021-04-26 2021-08-10 云南大学 一种优化增材制造方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《Research on the Highly Accelerated Stress Screening (HASS) Test Process of Civil Aircraft Airborne Equipment》;Peng Lu 等;《IOP Conf. Series: Materials Science and Engineering》;1043 *
《基于单元应力级数展开的下限原理有限元法》;李春光 等;《工程力学》;第32卷(第10期);38-43,59 *

Also Published As

Publication number Publication date
CN113704888A (zh) 2021-11-26

Similar Documents

Publication Publication Date Title
CN112100774A (zh) 一种基于变密度法的应力和应变能双约束的拓扑优化方法
CN105760619B (zh) 一种考虑疲劳可靠性的复合材料铺层优化设计方法
CN113704888B (zh) 一种单元应力筛选方法
CN106410073A (zh) 一种铝合金电池盖的制备方法、铝合金电池盖及电子设备
CN108564270A (zh) 一种贮存失效风险下的伽玛型单元备件需求量计算方法
CN108520145A (zh) 一种贮存失效风险下的威布尔型单元备件需求量计算方法
CN107336433A (zh) 一种3d打印自适应厚度分层方法
CN106011537A (zh) 一种细晶高强韧β钛合金及其制作方法
CN101512671B (zh) 水反应器燃料包壳管及其制造方法
CN114741915A (zh) 针对防隔热一体化结构热防护方案的优化方法及系统
CN117455222A (zh) 基于分布式异构流水车间组调度问题的求解方法
CN103870614B (zh) 一种结构概率优化设计方法
CN111274624B (zh) 一种基于rbf代理模型的多工况异形节点拓扑优化设计方法
CN104097098A (zh) 一种大型曲面薄壁件的多点定位方法
CN108090252A (zh) 一种数控铣床横梁的优化设计方法及优化设计的横梁
CN113887098B (zh) 一种基于疲劳应力鲁棒性和可靠性的金属结构拓扑优化方法
CN112749884A (zh) 一种飞行器电子元器件国产化替代风险决策方法
CN104897009A (zh) 一种耐烧蚀夹层复结构尾翼
CN109635315A (zh) 面向胚胎硬件细胞重用策略的可靠性建模和设计指导方法
CN112173168B (zh) 一种滑橇式起落架检查周期计算方法
CN110990939B (zh) 防冰腔结构可靠性设计方法
CN111021438A (zh) 一种低桩承台基础竖向承载力验算方法
Khurd et al. Analysis of two wheeler suspension spring by using fea for different materials
CN115203829A (zh) 一种基于截面系数的纵梁断面优化方法和纵梁
CN110175775A (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