CN106777617B - 一种二级箱间段结构优化设计方法 - Google Patents
一种二级箱间段结构优化设计方法 Download PDFInfo
- Publication number
- CN106777617B CN106777617B CN201611103436.9A CN201611103436A CN106777617B CN 106777617 B CN106777617 B CN 106777617B CN 201611103436 A CN201611103436 A CN 201611103436A CN 106777617 B CN106777617 B CN 106777617B
- Authority
- CN
- China
- Prior art keywords
- top connection
- tank section
- design
- optimization
- fork unit
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; 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 OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; 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)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Filling Or Discharging Of Gas Storage Vessels (AREA)
Abstract
本发明属于航天舱段结构设计领域,涉及一种二级箱间段结构优化设计方法,本发明主要包括二级箱间段结构的叉子单元上接头拓扑优化和设计优化两部分,具体为:根据结构载荷形式,对箱间段叉子单元进行拓扑优化,得到粗略结构形式;根据具体的结构载荷及总体结构设计要求,在拓扑优化得到的结构形式基础上进行设计优化,得到箱间段叉子单元接头的具体参数;对箱间段叉子单元及相邻的火箭壳体进行联合建模,校核二者互为真实边界后的结构强度;对上接头结构的具体几何参数进行优化,获得详细的结构参数。本发明能够提高杆系部段叉子单元的接头集中载荷的扩散效率,降低接头应力水平,解决箱间段叉子单元接头的局部强度问题。
Description
技术领域
本发明属于航天舱段结构设计领域,涉及一种二级箱间段结构优化设计方法。
背景技术
结构优化、特别是结构拓扑优化技术在各类结构设计中具有重要意义。基于先进的有限元数值分析技术发展起来的现代结构拓扑优化技术已经在航空航天、交通运输等部门得到了成功的应用,如在飞机等关键部件设计中获得了可观的优化收益。
芯二级氢箱和二级氧箱之间通过杆系形式的二级箱间段连接,在试验或者飞行中主要受轴拉/轴压载荷作用。对传统设计给出的结构形式经过有限元分析表明:箱间段杆系上接头在轴拉工况下发生三点弯曲,存在大面积的高应力区,不满足强度要求。为能得到满足设计要求的上接头结构形式,必须通过拓扑优化的手段重新设计上接头结构形式,使得该区域应力水平下降,达到预期性能指标。
发明内容
针对现有技术的不足,本发明提供一种二级箱间段结构优化设计方法,该方法能够提高杆系部段叉子单元的接头集中载荷的扩散效率,降低接头应力水平,解决箱间段叉子单元接头的局部强度问题。本发明适合于运载火箭杆系部段叉子单元的接头设计,有利于降低接头应力水平,提高结构承载能力。
为了达到上述目的,本发明的技术方案为:
一种二级箱间段结构优化设计方法,主要包括二级箱间段结构的叉子单元上接头拓扑优化和设计优化两部分,具体步骤如下:
第一步,对二级箱间段结构的叉子单元上接头进行拓扑优化
传统设计给出的叉子单元上接头结构,在部段轴拉载荷作用下,结构局部存在很高的集中应力,不满足强度要求。为能得到满足设计要求的上接头结构形式,本发明通过拓扑优化的手段重新设计上接头结构形式,使得该区域应力水平下降。
1.1)在HyperMesh中建立原二级箱间段结构的上接头模型,用一块实体的设计域包围原有上接头,建立上接头拓扑优化模型;并考虑螺栓的安装需求,在设计域中预留出螺帽所占用的空间。
1.2)简化上接头拓扑优化模型,将螺孔约束,并在上接头的轴向加载集中拉力,近似等效轴拉工况,完成上接头拓扑优化的模型建立工作;
1.3)因为拓扑优化只为寻找最佳的传力路径,据此给出概念设计,并不需要进行十分精确的非线性分析,考虑到接头的实际约束和受力情况,优化模型予以简化,直接将螺孔约束,并在上接头的轴向加载集中拉力来近似等效轴拉工况。
以上接头拓扑优化模型刚度最大(最小柔顺性)为设计目标,给定结构重量上限,对优化设计域内的材料密度分布进行设计,在OptiStruct中定义求解二级箱间段结构的最小柔顺性优化问题,得到上接头拓扑优化设计结果,即二级箱间段接头的具体参数。
依据OptiStruct软件中拓扑优化问题的求解算法(SIMP),构造最小柔顺性优化问题:
(a)目标:二级箱间段结构的最小柔顺性设计使结构整体刚度最大。
(b)约束:确定体积用量,作为约束上限。所述的体积用量为整体二级箱间段结构的10%。
(c)设计变量:优化设计域内单元密度。
1.4)根据步骤1.3)得到二级箱间段上接头的具体参数及其他部件原始参数,在有限元分析软件Abaqus中建立箱间段叉子单元及相邻火箭壳体的联合模型,进行有限元分析,得到有利于箱间段上接头集中应力的扩散的上接头结构形式,该上接头结构形式能够避免局部三点弯曲效应,降低接头应力水平,结构满足强度要求。
第二步,对二级箱间段结构的叉子单元上接头进行设计优化
在第一步1.4)得到的上接头结构形式的基础上,考虑上接头自身强度及总体要求、工艺等,对上接头结构的具体几何参数进行优化,获得详细的结构参数。经过设计优化后,叉子单元接头材料为钛合金ZTC4。该结构形式满足总体结构要求,构型比较简单,易于加工制造,且有利于装配以及集中应力扩散,经有限元分析和地面静力试验验证,该结构满足强度要求。
本发明的有益效果为:在叉子单元上接头优化前,该结构在轴拉载荷作用下存在三点弯曲问题,大面积材料应力、应变水平过高,不满足强度要求。采用本发明的“三爪”结构形式的上接头后,结构一次性通过地面静力试验考核,应力水平得到极大的降低,结构承载力明显提高。
附图说明
图1为拓扑优化模型;
图2为为上接头拓扑优化历程曲线;
图3为设计优化后叉子单元上接头结构形式。
具体实施方式
以下结合具体实施方式对本发明做进一步说明。
一种二级箱间段结构优化设计方法,主要包括二级箱间段结构的叉子单元上接头拓扑优化和设计优化两部分,整个过程如图2所示,步骤如下:
第一步,根据结构载荷形式,对二级箱间段结构的叉子单元进行拓扑优化,得到粗略结构形式,通过拓扑优化的手段重新设计上接头结构形式,使得该区域应力水平下降。
1.1)在HyperMesh中建立原二级箱间段结构的上接头模型,用一块实体的设计域包围原有上接头,建立上接头拓扑优化模型;并考虑螺栓的安装需求,在设计域中预留出螺帽所占用的空间。
1.2)简化上接头拓扑优化模型,将螺孔约束,并在上接头的轴向加载集中拉力,近似等效轴拉工况,完成上接头拓扑优化的模型建立工作;
1.3)以上接头拓扑优化模型刚度最大(最小柔顺性)为设计目标,给定结构重量上限,对优化设计域内的材料密度分布进行设计,在OptiStruct中定义求解二级箱间段结构的最小柔顺性优化问题,得到上接头拓扑优化设计结果,即二级箱间段接头的具体参数。
依据OptiStruct软件中拓扑优化问题的求解算法(SIMP),构造最小柔顺性优化问题:
(a)目标:二级箱间段结构的最小柔顺性设计使结构整体刚度最大。
(b)约束:确定体积用量,作为约束上限。所述的体积用量为整体二级箱间段结构的10%。
(c)设计变量:优化设计域内单元密度。
1.4)根据步骤1.3)得到二级箱间段上接头的具体参数及其他部件原始参数,在有限元分析软件Abaqus中建立箱间段叉子单元及相邻火箭壳体的联合模型,进行有限元分析,得到有利于箱间段上接头集中应力的扩散的上接头结构形式。
第二步,对二级箱间段结构的叉子单元上接头进行设计优化
在第一步1.4)得到的上接头结构形式的基础上,考虑上接头自身强度及总体要求、工艺等,对上接头结构的具体几何参数进行优化,获得详细的结构参数。结果表明上接头“三爪”结构的结构形式(如图3所示)有利于箱间段接头集中应力的扩散,避免了局部三点弯曲效应,降低了接头应力水平,结构满足强度要求。
Claims (1)
1.一种二级箱间段结构优化设计方法,其特征在于以下步骤:
第一步,对二级箱间段结构的叉子单元上接头进行拓扑优化,通过拓扑优化的手段重新设计上接头结构形式,使得上接头结构应力水平下降;
1.1)建立原二级箱间段结构的上接头模型,用一块实体的设计域包围原有上接头,建立上接头拓扑优化模型;在设计域中预留螺帽的空间,用于安装螺栓;
1.2)简化上接头拓扑优化模型,将螺孔约束,并在上接头的轴向加载集中拉力,近似等效轴拉工况,完成上接头拓扑优化的模型建立工作;
1.3)以上接头拓扑优化模型刚度最大、最小柔顺性为设计目标,给定二级箱间段结构重量上限,对优化设计域内的材料密度分布进行设计,求解二级箱间段结构的最小柔顺性优化问题,得到上接头拓扑优化设计结果,即二级箱间段接头的具体参数;
1.4)根据步骤1.3)得到二级箱间段上接头的具体参数及其他部件原始参数,在有限元分析软件中建立箱间段叉子单元及相邻火箭壳体的联合模型,进行有限元分析,得到有利于箱间段上接头集中应力的扩散的上接头结构形式;
第二步,对二级箱间段结构的叉子单元上接头进行设计优化;
在第一步1.4)得到的上接头结构形式的基础上,考虑上接头自身强度及总体要求、工艺,对上接头结构的几何参数进行优化,获得结构参数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611103436.9A CN106777617B (zh) | 2016-12-05 | 2016-12-05 | 一种二级箱间段结构优化设计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611103436.9A CN106777617B (zh) | 2016-12-05 | 2016-12-05 | 一种二级箱间段结构优化设计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106777617A CN106777617A (zh) | 2017-05-31 |
| CN106777617B true CN106777617B (zh) | 2019-06-14 |
Family
ID=58883965
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611103436.9A Expired - Fee Related CN106777617B (zh) | 2016-12-05 | 2016-12-05 | 一种二级箱间段结构优化设计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106777617B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109325256B (zh) * | 2018-08-16 | 2023-01-17 | 上海宇航系统工程研究所 | 偏置载荷助推受力均匀的一二级连接部段结构设计方法 |
| CN110705084B (zh) * | 2019-09-26 | 2023-05-05 | 内蒙动力机械研究所 | 一种复合材料壳体的快速设计软件系统 |
| CN114492130B (zh) * | 2022-01-26 | 2024-10-11 | 上海交通大学 | 一种基于杆系结构模型的载荷扩散结构优化方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009004364A2 (en) * | 2007-06-29 | 2009-01-08 | Airbus Uk Limited | Improvements in elongate composite structural member |
| CN104359694A (zh) * | 2014-11-06 | 2015-02-18 | 北京强度环境研究所 | 液体火箭系统级试车工作模态及工作变形试验方法 |
| CN105468822A (zh) * | 2015-11-16 | 2016-04-06 | 北京宇航系统工程研究所 | 一种火箭助推分离冲击载荷确定方法 |
-
2016
- 2016-12-05 CN CN201611103436.9A patent/CN106777617B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009004364A2 (en) * | 2007-06-29 | 2009-01-08 | Airbus Uk Limited | Improvements in elongate composite structural member |
| CN104359694A (zh) * | 2014-11-06 | 2015-02-18 | 北京强度环境研究所 | 液体火箭系统级试车工作模态及工作变形试验方法 |
| CN105468822A (zh) * | 2015-11-16 | 2016-04-06 | 北京宇航系统工程研究所 | 一种火箭助推分离冲击载荷确定方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106777617A (zh) | 2017-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109583058B (zh) | 基于有限元方法的汽车电池箱结构设计方法 | |
| CN109726484B (zh) | 基于独立连续映射法的多材料连续体结构拓扑优化设计方法 | |
| CN106777617B (zh) | 一种二级箱间段结构优化设计方法 | |
| CN113609584B (zh) | 重复使用火箭发动机推力室可靠性优化设计方法 | |
| CN111191401B (zh) | 高超声速飞行器气动/控制/结构多学科协同优化方法 | |
| CN110175430B (zh) | 一种叶片极限状态下多目标自动优化方法 | |
| Lujano-Rojas et al. | Contract design of direct-load control programs and their optimal management by genetic algorithm | |
| CN107169209B (zh) | 一种偏置集中力作用下的结构设计及优化方法 | |
| CN106650045B (zh) | 一种二级箱间段结构优化设计与强度分析方法 | |
| CN108491576B (zh) | 一种复合材料机翼开口补强的优化设计方法 | |
| CN111737908B (zh) | 一种基于动载荷静力等效的蒙皮桁条结构快速动态优化设计方法 | |
| CN117669016A (zh) | 全局/局部耦合的宽速域乘波-机翼布局的优化设计方法 | |
| CN106202693A (zh) | 一种基于参数化建模的加筋壁板结构抗振动疲劳优化方法 | |
| CN114282310B (zh) | 一种基于自适应加点代理模型的气弹结构耦合寻优方法 | |
| Bortolotti et al. | An efficient approach to explore the solution space of a wind turbine rotor design process | |
| CN109766638B (zh) | 一种燃油喷嘴及其轻量化设计方法 | |
| CN105335549A (zh) | 一种上下翼面变形叠加方法 | |
| CN111859545A (zh) | 一种考虑升力匹配的宽速域高升阻比机翼优化设计方法 | |
| CN115859717A (zh) | 一种网壳结构装配式节点的拓扑优化设计方法 | |
| CN115663818A (zh) | 交直流配电网分布式凸优化方法及装置 | |
| CN110826182B (zh) | 基于顶点法和序贯优化策略的飞行器结构气动弹性设计法 | |
| CN108248825A (zh) | 一种飞行器翼面结构 | |
| CN113111562B (zh) | 基于变形技术的风电机组主机架优化设计方法 | |
| CN114462260B (zh) | 基于超级计算机的磁约束聚变融合数值处理系统 | |
| CN115470591A (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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190614 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |