CN117540494B - 一种航空舱段结构降阶分析方法、系统及设备 - Google Patents
一种航空舱段结构降阶分析方法、系统及设备 Download PDFInfo
- Publication number
- CN117540494B CN117540494B CN202410032291.6A CN202410032291A CN117540494B CN 117540494 B CN117540494 B CN 117540494B CN 202410032291 A CN202410032291 A CN 202410032291A CN 117540494 B CN117540494 B CN 117540494B
- Authority
- CN
- China
- Prior art keywords
- unit cell
- working condition
- generalized
- freedom
- bending
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 30
- 238000004088 simulation Methods 0.000 claims abstract description 117
- 239000011159 matrix material Substances 0.000 claims abstract description 93
- 230000000737 periodic effect Effects 0.000 claims abstract description 76
- 238000000034 method Methods 0.000 claims abstract description 50
- 230000004044 response Effects 0.000 claims abstract description 48
- 239000013598 vector Substances 0.000 claims abstract description 33
- 210000004027 cell Anatomy 0.000 claims description 256
- 238000005452 bending Methods 0.000 claims description 93
- 238000010008 shearing Methods 0.000 claims description 77
- 238000006073 displacement reaction Methods 0.000 claims description 60
- 230000009471 action Effects 0.000 claims description 18
- 230000003993 interaction Effects 0.000 claims description 17
- 230000008878 coupling Effects 0.000 claims description 16
- 238000010168 coupling process Methods 0.000 claims description 16
- 238000005859 coupling reaction Methods 0.000 claims description 16
- 238000013519 translation Methods 0.000 claims description 14
- 230000001808 coupling effect Effects 0.000 claims description 13
- 230000000694 effects Effects 0.000 claims description 13
- 238000004590 computer program Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 6
- 230000007935 neutral effect Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 238000005457 optimization Methods 0.000 claims description 4
- 238000003860 storage Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract description 6
- 238000004364 calculation method Methods 0.000 description 14
- 238000000265 homogenisation Methods 0.000 description 9
- 239000000463 material Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 230000006399 behavior Effects 0.000 description 3
- 239000011365 complex material Substances 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000012916 structural analysis Methods 0.000 description 1
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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
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
本发明公开了一种航空舱段结构降阶分析方法、系统及设备，涉及航空板壳结构设计技术领域。方法包括：构建单胞模型并仿真，进而确定当前单胞结构的等效刚度矩阵；基于所述等效刚度矩阵和所述周期性排列规则，利用有限元建模软件构建航空舱段结构的整体简化模型；以实际工况下的宏观载荷为条件，对整体简化模型进行有限元仿真，得到航空舱段结构中关注细节区域处的广义应力向量以确定航空舱段结构的细节力学响应，进而优化所述航空舱段结构。本发明通过构建航空舱段结构的降阶模型并进行多尺度分析，在维持仿真精度的基础上降低计算量，进而提高航空舱段结构的仿真分析效率。
Description
技术领域
本发明涉及航空板壳结构设计技术领域，特别是涉及一种航空舱段结构降阶分析方法、系统及设备。
背景技术
航空舱段多为板壳结构，板壳结构是航空工程应用中常用的轻量化结构，随着制造、加工工艺的进步，现代航空板壳结构通常包括多种材料和复杂拓扑结构。在结构设计特别是初始设计阶段，建立包含详细材料的微观结构的全尺寸单元进行仿真往往是不经济甚至不可行的。这使得结构细节的变形响应情况难以准确获得，预测承受复杂耦合载荷作用下的结构响应也变得困难。因此，应用有显著更低自由度的均质壳模型以预测结构的整体力学响应往往是必要的。经典的板壳理论，如Kirchhoff理论和Reissner理论以及现有有限元软件中相应的壳单元已被广泛用于轻量化结构的分析和设计。值得注意的是，对于有简单力学行为的结构（例如层合板和功能梯度板）可以进行解析或半解析的板壳模型降阶，但对于具有复杂材料或拓扑结构的结构（如由加筋层合板或夹层板组成的飞机机身截面）仍然需要发展均匀化方法来计算壳体的等效截面性能。
近年来，一些研究人员在板壳均质化领域做出了显著的贡献。Hodges等人将变分渐近梁截面分析方法扩展到板壳结构的均匀化，并在此基础上进一步发展了结构基因力学。此外，在具有三向周期性的材料渐近均匀化理论的基础上，有学者发明了相应的板壳渐进均匀化方法。然而，当考虑剪切变形和力时，具有周期单元格的板壳结构不表现出周期力学行为，并且不能直接使用周期边界条件和Hill条件进行均匀化分析。Hodges等人发明的变分渐近方法仍然主要应用于梁的截面分析。在变分渐近方法和随后发展的结构基因组力学中，也没有给出与周期性边界条件不同的保持物理一致性的边界条件。至于壳的渐近均匀化方法，徐亮等人认识到周期边界条件的不适用性，提出了一种假设应变场为线性变化的壳的渐近均匀化方法。然而，周期性结构单元的非周期力学响应是由面外剪力而不是剪切应变引起的。因此，其基于应变的均匀化方法无法处理弯曲与剪切存在耦合或两个方向上的剪切存在耦合作用的结构。
综上，目前需要发展一套方法，可以在保持与原始三维结构实际力学行为一致的同时，实现准确而有效的板壳结构降阶和多尺度分析，进而完成优化航空舱段结构的设计流程。
发明内容
本发明的目的是提供一种航空舱段结构降阶分析方法、系统及设备，能够构建航空舱段结构的降阶模型并进行多尺度分析，在维持仿真精度的基础上降低计算量，进而提高航空舱段结构的优化效率。
为实现上述目的，本发明提供了如下方案：
一种航空舱段结构降阶分析方法，包括：
确定航空舱段结构中任一单胞结构为当前单胞结构；所述航空舱段结构为筒状结构；所述航空舱段结构是将多个所述单胞结构按照周期性排列规则进行排列的得到的；所述周期性排列规则为所述单胞结构沿所述航空舱段的端面圆周呈周期性排列；所述单胞结构沿所述航空舱段的径向呈周期性排列；
利用有限元建模软件对当前单胞结构进行有限元建模处理，得到单胞模型；
基于工况种类调整单胞模型的周期性边界条件对单胞模型进行仿真，得到单胞仿真结果；所述单胞仿真结果包括在不同虚拟节点自由度上施加单位载荷时单胞结构的位移场、应变场和应力场；
基于所述单胞仿真结果，确定当前单胞结构的等效刚度矩阵；
基于所述等效刚度矩阵和所述周期性排列规则，利用有限元建模软件构建航空舱段结构的整体简化模型；
获取航空舱段结构中的关注细节区域；所述关注细节区域的面积小于所述单胞结构的面积；
以实际工况下的宏观载荷为条件，对整体简化模型进行有限元仿真，得到航空舱段结构中关注细节区域处的广义应力向量；
根据所述单胞仿真结果和所述广义应力向量，确定航空舱段结构中关注细节区域处的力学响应；所述力学响应包括位移场、应变场和应力场；
基于关注细节区域处的力学响应，优化所述航空舱段结构。
可选的，所述航空舱段结构具有周期性特点。
可选的，利用有限元建模软件对当前单胞结构进行有限元建模处理，得到单胞模型，包括：
利用有限元建模软件对当前单胞结构进行网格划分，得到多个网格节点；
以当前单胞结构中心点为原点O，以第一周期性方向为方向，以第二周期性方向为/>方向，以垂直于/>平面的方向为/>方向，构建三维坐标系；
确定第一正端面和第一负端面；所述第一正端面和所述第一负端面在所述第一周期性方向上平行设置；
将所述第一正端面和所述第一负端面上方向坐标相等且/>方向坐标相等的网格节点进行配对处理；
确定第二正端面和第二负端面；所述第二正端面和所述第二负端面在所述第二周期性方向上平行设置；
将所述第二正端面和所述第二负端面上方向坐标相等且/>方向坐标相等的网格节点进行配对处理；
基于板壳广义应变类型设置单胞模型中的虚拟节点自由度；所述虚拟节点自由度与板壳广义应变类型一一对应；所述广义应变包括方向拉伸、/>方向拉伸、面内剪切、方向弯曲、/>方向弯曲、扭转、/>方向面外剪切和/>方向面外剪切；
获取参考点和参考直线；所述参考直线为当前单胞结构中过坐标原点且垂直于板壳中性面的任一直线；
基于单胞刚体位移约束条件，设置参考点和参考直线处的参数，约束当前单胞结构的刚体位移；所述单胞刚体位移约束条件为设置参考点处方向平动、/>方向平动、以方向平动、/>方向转动、/>方向转动和以/>方向转动均为0，并向参考直线施加分布耦合约束。
可选的，基于工况种类调整单胞模型的周期性边界条件对单胞模型进行仿真，得到单胞仿真结果，包括：
基于周期性边界条件，分别在第一工况下和第二工况下对单胞模型进行仿真，得到广义应力作用下的单胞仿真结果；
基于周期性边界条件，分别在第三工况下和第四工况下对单胞模型进行仿真，得到广义弯矩作用下的单胞仿真结果；
根据广义弯矩作用下的单胞仿真结果确定扰动面力场施加量；
基于扰动面力场施加量和修正周期性边界条件，分别在第五工况下、第六工况下、第七工况下、第八工况下和第九工况下对单胞模型进行仿真，得到广义应力耦合作用下的单胞仿真结果。
可选的，所述第一工况为在方向拉伸、/>方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲和扭转对应的虚拟节点自由度上均同时施加单位载荷；
所述第二工况在所有虚拟节点自由度上分别施加单位载荷；
所述第三工况为在方向弯曲对应的虚拟节点自由度上施加单位载荷，同时在方向拉伸、/>方向弯曲、扭转和/>方向面外剪切对应的虚拟节点自由度上均施加0载荷；
所述第四工况为在方向弯曲对应的虚拟节点自由度上施加单位载荷，同时在方向拉伸、/>方向弯曲、扭转和/>方向面外剪切对应的虚拟节点自由度上均施加0载荷；
所述第五工况为在方向面外剪切对应的虚拟节点自由度上施加单位载荷；
所述第六工况为在方向面外剪切对应的虚拟节点自由度上施加单位载荷；
所述第七工况为在方向面外剪切和/>方向面外剪切对应的虚拟节点自由度上同时施加单位载荷；
所述第八工况为在方向面外剪切对应的虚拟节点自由度上施加单位载荷，依次分别在/>方向拉伸、/>方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲和扭转对应的虚拟节点自由度上施加单位载荷；
所述第九工况为在方向面外剪切对应的虚拟节点自由度上施加单位载荷，依次分别在/>方向拉伸、/>方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲和扭转对应的虚拟节点自由度上施加单位载荷。
可选的，所述广义应力作用下的单胞仿真结果包括：第一工况下6个虚拟节点自由度的位移，以及第二工况下8个虚拟节点自由度的力学响应；
所述广义弯矩作用下的单胞仿真结果包括：第三工况下单胞扰动位移场、单胞面力和方向弯曲对应的虚拟节点自由度上的位移，以及第四工况下单胞扰动位移场、单胞面力和/>方向弯曲对应的虚拟节点自由度上的位移；
所述广义应力耦合作用下的单胞仿真结果包括第五工况下方向面外剪切对应的虚拟节点自由度的应变能、第六工况下/>方向面外剪切对应的虚拟节点自由度的应变能、第七工况下的/>方向面外剪切和/>方向面外剪切相互作用时的单胞应变能、第八工况下/>方向面外剪切与/>方向拉伸、/>方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲和扭转中任一广义应变类型相互作用时的单胞应变能和第九工况下/>方向面外剪切与/>方向拉伸、/>方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲和扭转中任一广义应变类型相互作用时的单胞应变能。
可选的，基于所述单胞仿真结果，确定当前单胞结构的等效刚度矩阵，包括：
根据广义应变-广义应力关系，构建初始柔度矩阵；所述初始柔度矩阵为8阶矩阵；
确定所述初始柔度矩阵中的第一参数、第二参数、第三参数和第四参数；所述第一参数为初始柔度矩阵中的元素；所述第二参数为初始柔度矩阵中的元素/>；所述第三参数为初始柔度矩阵中的元素；所述第四参数为初始柔度矩阵中除所述第一参数、所述第二参数和所述第三参数之外的参数；
基于所述广义应力作用下的单胞仿真结果，确定所述初始柔度矩阵中的第一参数；
基于所述广义弯矩作用下的单胞仿真结果、所述第五工况下方向面外剪切对应的虚拟节点自由度的应变能和所述第六工况下/>方向面外剪切对应的虚拟节点自由度的应变能，确定所述初始柔度矩阵中的第二参数；
基于所述第七工况下的方向面外剪切和/>方向面外剪切相互作用时的单胞应变能，确定所述初始柔度矩阵中的第三参数；
基于第八工况下方向面外剪切与/>方向拉伸、/>方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲和扭转中任一广义应变类型相互作用时的单胞应变能，以及第九工况下方向面外剪切与/>方向拉伸、/>方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲和扭转中任一广义应变类型相互作用时的单胞应变能，确定所述初始柔度矩阵中的第四参数，得到柔度矩阵；
确定所述柔度矩阵的逆矩阵为当前单胞结构的等效刚度矩阵。
可选的，根据所述单胞仿真结果和所述广义应力向量，确定航空舱段结构中关注细节区域处的力学响应，包括：
确定任一广义应力类型为当前广义应力类型；
确定当前广义应力类型对应的广义应力向量分量，与确定第二工况下当前广义应力类型对应虚拟节点自由度的力学响应之积，为航空舱段结构中关注细节区域处在当前广义应力类型对应虚拟节点自由度上的力学响应；所述广义应力向量包括多个广义应力向量分量；所述广义应力向量分量与所述广义应力类型一一对应。
一种航空舱段结构降阶分析系统，包括：
当前单胞结构确定模块，用于确定航空舱段结构中任一单胞结构为当前单胞结构；所述航空舱段结构为筒状结构；所述航空舱段结构是将多个所述单胞结构按照周期性排列规则进行排列的得到的；所述周期性排列规则为所述单胞结构沿所述航空舱段的端面圆周呈周期性排列；所述单胞结构沿所述航空舱段的径向呈周期性排列；
单胞模型构建模块，用于利用有限元建模软件对当前单胞结构进行有限元建模处理，得到单胞模型；
单胞仿真模块，用于基于工况种类调整单胞模型的周期性边界条件对单胞模型进行仿真，得到单胞仿真结果；所述单胞仿真结果包括在不同虚拟节点自由度上施加单位载荷时单胞结构的位移场、应变场和应力场；
等效刚度矩阵确定模块，用于基于所述单胞仿真结果，确定当前单胞结构的等效刚度矩阵；
整体简化模型构建模块，用于基于所述等效刚度矩阵和所述周期性排列规则，利用有限元建模软件构建航空舱段结构的整体简化模型；
关注细节区域获取模块，用于获取航空舱段结构中的关注细节区域；所述关注细节区域的面积小于所述单胞结构的面积；
广义应力向量确定模块，用于以实际工况下的宏观载荷为条件，对整体简化模型进行有限元仿真，得到航空舱段结构中关注细节区域处的广义应力向量；
细节力学响应确定模块，用于根据所述单胞仿真结果和所述广义应力向量，确定航空舱段结构中关注细节区域处的力学响应；所述力学响应包括位移场、应变场和应力场；
一种电子设备，包括存储器及处理器，所述存储器用于存储计算机程序，所述处理器运行所述计算机程序以使所述电子设备执行所述的一种航空舱段结构降阶分析方法；所述存储器为可读存储介质。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明提供的一种航空舱段结构降阶分析方法、系统及设备，通过仿真得到航空板壳结构单胞在不同载荷工况下的变形，进一步计算等效降阶板壳模型的本构矩阵（等效刚度矩阵），进而为飞机宏观板壳模型有限元分析提供降阶模型单元，将复杂拓扑、复杂材料组分的三维结构进行简化建模，处理成仅包含少量板壳单元的仿真模型，从而高效地对整机模型进行仿真分析，快速准确地预测各种复杂板壳类结构的力学响应，显著提升结构分析与仿真的效率。本发明具有良好的通用性，适用于超材料、复合材料、夹芯、加筋等多种形式的航空板壳结构，且随着进一步发展，相关方法在屈曲、振动、多场分析等方面具有很大的应用潜力。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例1中航空舱段结构降阶分析方法流程图；
图2为本发明实施例2中航空舱段结构降阶分析方法流程图；
图3为本发明实施例2中航空舱段结构三维单胞模型示意图；
图4为本发明实施例2中整体舱段结构全尺寸模型示意图；
图5为本发明实施例2中舱段结构受载情况示意图；
图6为本发明实施例2中板壳简化模型扰度仿真结果示意图；
图7为本发明实施例2中结构简化模型与全尺寸模型分析结果对比图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的目的是提供一种航空舱段结构降阶分析方法、系统及设备，能够构建航空舱段结构的降阶模型并进行多尺度分析，在维持仿真精度的基础上降低计算量，进而提高航空舱段结构的优化效率。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例1。
如图1所示，本实施例提供了一种航空舱段结构降阶分析方法，包括：
步骤101：确定航空舱段结构中任一单胞结构为当前单胞结构。航空舱段结构为筒状结构。航空舱段结构是将多个单胞结构按照周期性排列规则进行排列的得到的。周期性排列规则为单胞结构沿航空舱段的端面圆周呈周期性排列。单胞结构沿航空舱段的径向呈周期性排列。航空舱段结构具有周期性特点。
步骤102：利用有限元建模软件对当前单胞结构进行有限元建模处理，得到单胞模型。
步骤103：基于工况种类调整单胞模型的周期性边界条件对单胞模型进行仿真，得到单胞仿真结果。单胞仿真结果包括在不同虚拟节点自由度上施加单位载荷时单胞结构的位移场、应变场和应力场。
步骤104：基于单胞仿真结果，确定当前单胞结构的等效刚度矩阵。
步骤105：基于等效刚度矩阵和周期性排列规则，利用有限元建模软件构建航空舱段结构的整体简化模型。
步骤107：以实际工况下的宏观载荷为条件，对整体简化模型进行有限元仿真，得到航空舱段结构中关注细节区域处的广义应力向量。
步骤108：根据单胞仿真结果和广义应力向量，确定航空舱段结构中关注细节区域处的力学响应。力学响应包括位移场、应变场和应力场。
步骤102，包括：
步骤102-1：利用有限元建模软件对当前单胞结构进行网格划分，得到多个网格节点。
步骤102-2：以当前单胞结构中心点为原点O，以第一周期性方向为方向，以第二周期性方向为/>方向，以垂直于/>平面的方向为/>方向，构建三维坐标系。
步骤102-3：确定第一正端面和第一负端面。第一正端面和第一负端面在第一周期性方向上平行设置。
步骤102-4：将第一正端面和第一负端面上方向坐标相等且/>方向坐标相等的网格节点进行配对处理。
步骤102-5：确定第二正端面和第二负端面。第二正端面和第二负端面在第二周期性方向上平行设置。
步骤102-6：将第二正端面和第二负端面上方向坐标相等且/>方向坐标相等的网格节点进行配对处理。
步骤102-7：基于板壳广义应变类型设置单胞模型中的虚拟节点自由度。虚拟节点自由度与板壳广义应变类型一一对应。广义应变包括方向拉伸、/>方向拉伸、面内剪切、方向弯曲、/>方向弯曲、扭转、/>方向面外剪切和/>方向面外剪切。
步骤102-8：获取参考点和参考直线。参考直线为当前单胞结构中过坐标原点且垂直于板壳中性面的任一直线。
步骤102-9：基于单胞刚体位移约束条件，设置参考点和参考直线处的参数，约束当前单胞结构的刚体位移。单胞刚体位移约束条件为设置参考点处方向平动、/>方向平动、以/>方向平动、/>方向转动、/>方向转动和以/>方向转动均为0，并向参考直线施加分布耦合约束。
步骤103，包括：
步骤103-1：基于周期性边界条件，分别在第一工况下和第二工况下对单胞模型进行仿真，得到广义应力作用下的单胞仿真结果。
步骤103-2：基于周期性边界条件，分别在第三工况下和第四工况下对单胞模型进行仿真，得到广义弯矩作用下的单胞仿真结果。
步骤103-3：根据广义弯矩作用下的单胞仿真结果确定扰动面力场施加量。
步骤103-4：基于扰动面力场施加量和修正周期性边界条件，分别在第五工况下、第六工况下、第七工况下、第八工况下和第九工况下对单胞模型进行仿真，得到广义应力耦合作用下的单胞仿真结果。
其中，第一工况为在方向拉伸、/>方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲和扭转对应的虚拟节点自由度上均同时施加单位载荷。
第二工况在所有虚拟节点自由度上分别施加单位载荷；
第三工况为在方向弯曲对应的虚拟节点自由度上施加单位载荷，同时在/>方向拉伸、/>方向弯曲、扭转和/>方向面外剪切对应的虚拟节点自由度上均施加0载荷。
第四工况为在方向弯曲对应的虚拟节点自由度上施加单位载荷，同时在/>方向拉伸、/>方向弯曲、扭转和/>方向面外剪切对应的虚拟节点自由度上均施加0载荷。
第五工况为在方向面外剪切对应的虚拟节点自由度上施加单位载荷。
第六工况为在方向面外剪切对应的虚拟节点自由度上施加单位载荷。
第七工况为在方向面外剪切和/>方向面外剪切对应的虚拟节点自由度上同时施加单位载荷。
第八工况为在方向面外剪切对应的虚拟节点自由度上施加单位载荷，依次分别在/>方向拉伸、/>方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲和扭转对应的虚拟节点自由度上施加单位载荷。
第九工况为在方向面外剪切对应的虚拟节点自由度上施加单位载荷，依次分别在/>方向拉伸、/>方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲和扭转对应的虚拟节点自由度上施加单位载荷。
具体的，广义应力作用下的单胞仿真结果包括：第一工况下6个虚拟节点自由度的位移，以及第二工况下8个虚拟节点自由度的力学响应。
广义弯矩作用下的单胞仿真结果包括：第三工况下单胞扰动位移场、单胞面力和方向弯曲对应的虚拟节点自由度上的位移，以及第四工况下单胞扰动位移场、单胞面力和/>方向弯曲对应的虚拟节点自由度上的位移。
广义应力耦合作用下的单胞仿真结果包括第五工况下方向面外剪切对应的虚拟节点自由度的应变能、第六工况下/>方向面外剪切对应的虚拟节点自由度的应变能、第七工况下的/>方向面外剪切和/>方向面外剪切相互作用时的单胞应变能、第八工况下/>方向面外剪切与/>方向拉伸、/>方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲和扭转中任一广义应变类型相互作用时的单胞应变能和第九工况下/>方向面外剪切与/>方向拉伸、/>方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲和扭转中任一广义应变类型相互作用时的单胞应变能。
步骤104，包括：
步骤104-1：根据广义应变-广义应力关系，构建初始柔度矩阵。初始柔度矩阵为8阶矩阵。
步骤104-2：确定初始柔度矩阵中的第一参数、第二参数、第三参数和第四参数。第一参数为初始柔度矩阵中的元素；第二参数为初始柔度矩阵中的元素/>；第三参数为初始柔度矩阵中的元素；第四参数为初始柔度矩阵中除第一参数、第二参数和第三参数之外的参数。
步骤104-3：基于广义应力作用下的单胞仿真结果，确定初始柔度矩阵中的第一参数。
步骤104-4：基于广义弯矩作用下的单胞仿真结果、第五工况下方向面外剪切对应的虚拟节点自由度的应变能和第六工况下/>方向面外剪切对应的虚拟节点自由度的应变能，确定初始柔度矩阵中的第二参数。
步骤104-5：基于第七工况下的方向面外剪切和/>方向面外剪切相互作用时的单胞应变能，确定初始柔度矩阵中的第三参数。
步骤104-6：基于第八工况下方向面外剪切与/>方向拉伸、/>方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲和扭转中任一广义应变类型相互作用时的单胞应变能，以及第九工况下/>方向面外剪切与/>方向拉伸、/>方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲和扭转中任一广义应变类型相互作用时的单胞应变能，确定初始柔度矩阵中的第四参数，得到柔度矩阵。
步骤104-7：确定柔度矩阵的逆矩阵为当前单胞结构的等效刚度矩阵。
步骤108，包括：
步骤108-1：确定任一广义应力类型为当前广义应力类型。
步骤108-2：确定当前广义应力类型对应的广义应力向量分量，与确定第二工况下当前广义应力类型对应虚拟节点自由度的力学响应之积，为航空舱段结构中关注细节区域处在当前广义应力类型对应虚拟节点自由度上的力学响应。广义应力向量包括多个广义应力向量分量。广义应力向量分量与广义应力类型一一对应。
实施例2。
如图2，所示，本实施例提供了一种基于三维单胞的航空舱段结构优化方法如下。
步骤1.航空三维板壳结构单胞建模。
如图3，进行分析对象的三维结构单胞建模，具体步骤如下。
步骤1.1根据板壳的内部结构、材料布局和周期性特点，选择适当单胞进行建模，划分有限元网格，航空舱段结构如图4。
步骤1.2对单胞模型的两个周期性方向上的端面网格节点进行节点配对。
步骤1.3设置与板壳广义应变相对应的虚拟节点自由度。
步骤1.4约束单胞的刚体位移。约束单胞刚体位移不限于将参考点与过坐标原点的线耦合，通过与其他点耦合，并令参考点位移等于单胞此处刚体平动及转动大小亦可实现刚体约束。
步骤2.单胞多工况变形仿真。
对多种典型工况下的结构单胞响应进行仿真计算，为等效板壳属性计算和结构细节应力恢复奠定基础，具体步骤如下。
步骤2.1仿真计算包含拉、弯、扭、面内剪切在内的6个单位广义应力作用工况下的板壳单胞的应力、应变响应。
步骤2.1.1基于方程约束施加周期性边界条件。
步骤2.1.2施加各方向单位拉、弯、扭、面内剪切载荷，进行有限元计算。
步骤2.2仿真计算用于计算扰动场的包含额外约束的广义弯矩作用下的2个工况下的板壳单胞的应力、应变响应。
步骤2.2.1基于方程约束施加周期性边界条件。
步骤2.2.2通过约束虚拟节点自由度的方式施加额外广义应变约束。
步骤2.2.3施加各方向单位弯曲载荷，进行有限元计算。
步骤2.3仿真计算2个面外剪切以及13个面外剪切与其他广义应力耦合作用工况下的板壳单胞的应力、应变响应。
步骤2.3.1基于方程约束施加包含扰动场修正的周期性边界条件。
步骤2.3.2在周期性边界面上施加面力扰动场。
步骤2.3.3施加各方向单位剪切载荷，进行有限元计算。
步骤2.3.4施加两个方向剪切耦合作用的载荷，进行有限元计算。
步骤2.3.5施加剪切与其他广义应力耦合作用的载荷，进行有限元计算。
步骤3.等效板壳属性计算。
计算结构等效板壳刚度矩阵，具体步骤如下。
步骤3.1基于步骤2.1.2所得有限元仿真结果，计算柔度矩阵左上角的6*6阶薄板柔度矩阵系数。柔度矩阵左上角的6*6阶矩阵柔度系数的求解方式不限于基于应力-应变关系式计算，通过类似步骤3.2的应变能计算方式亦可实现求解。
步骤3.2基于步骤2.3.3所得有限元仿真结果，计算柔度矩阵对角线上的两个剪切柔度系数。
步骤3.3基于步骤2.3.4所得有限元仿真结果，计算柔度矩阵两个方向剪切耦合项系数。
步骤3.4基于步骤2.3.5所得有限元仿真结果，计算柔度矩阵中的其余耦合项系数。
步骤3.5对所得柔度矩阵求逆，得到相应等效板壳刚度矩阵。
步骤4.板壳简化模型仿真。
利用所得到的等效板壳刚度矩阵，应用有限元仿真软件中的板壳单元建立原三维结构的低自由度简化板壳模型，施加待分析的结构服役过程中实际工况下的宏观载荷，进行有限元仿真。
步骤5.板壳结构单胞细节应力恢复。
利用上述所得有限元仿真结果，计算并输出板壳结构单胞的细节应力场。
进一步地，所述步骤1.2中的节点配对步骤特征在于：板壳结构单胞的有限元网格模型，在两个周期性方向上形成两对相平行的端面，以单胞中心点为坐标原点。在网格为对称离散的情况下，方向正端面的网格节点与/>方向负端面上有相同/>坐标的节点配对，相应地，/>方向正端面的网格节点与/>方向负端面上有相同/>坐标的节点配对。其中，/>代表三个坐标方向，/>分别为网格节点在/>方向的坐标值。
进一步地，所述步骤1.3中的设置虚拟节点自由度的特征在于：设置8个虚拟节点自由度，定义第个虚拟节点自由度的位移大小为/>，其数值大小分别对应板壳的8个广义应变：/>，分别代表/>方向拉伸，面内剪切，方向弯曲，扭转，/>方向面外剪切。
进一步地，所述步骤1.4中的约束单胞刚体位移的特征在于：实现方法为设置参考点，将参考点与板壳单胞中过坐标原点的垂直于板壳中性面的一条线施加分布耦合约束，同时令该参考点的三个方向平动及三个转动为0。
进一步地，所述步骤2.1.1、步骤2.2.1中的周期性边界条件，其方程约束形式为公式（1）。
（1）。
其中，分别表示/>方向正负端面上的点在/>方向上的位移，/>分别表示单胞在/>方向上的长度。
进一步地，所述步骤2.1.2、步骤2.2.3、步骤2.3.3中施加载荷的实现方式为在步骤1.3所设置虚拟节点自由度上施加载荷，在其8个分量上所施加的载荷大小即分别等于对板壳施加的8个广义应力：，其分别是/>方向拉力，面内剪力，/>方向弯矩，扭矩，/>方向面外剪力。
进一步地，在所述步骤2.1.2中施加的单位载荷是分别依次施加在所对应的节点自由度上，在所述步骤2.2.3中施加的单位载荷是分别依次施加在所对应的节点自由度上，在所述步骤2.3.3中施加的单位载荷是分别依次施加在/>所对应的节点自由度上。
进一步地，所述步骤2.2.2中通过约束虚拟节点自由度的方式施加额外广义应变约束：在对应虚拟节点自由度施加单位载荷的工况中，施加额外约束为令广义应变对应的虚拟节点自由度为0。在/>对应虚拟节点自由度施加单位载荷的工况中，施加额外约束为令广义应变/>对应的虚拟节点自由度为0。/>
进一步地，所述步骤2.3.1中包含扰动场修正的周期性边界条件，其约束方程为公式（2）。
（2）。
其中，为扰动位移场，分别等于由步骤2.2.3仿真得到的单胞在/>对应工况下/>方向正端面上的节点在/>方向上的位移，/>分别为该工况中在/>所对应的节点自由度上施加的载荷大小。
进一步地，所述步骤2.3.2中在周期性边界面上施加面力扰动场的实现方式为在单胞在方向正、负端面节点的/>方向上施加大小为/>的扰动面力场，其中/>具体数值等于步骤2.2.3仿真得到的单胞在/>对应工况下/>方向正端面上的节点在/>方向上的面力。
进一步地，所述步骤2.3.4中施加的耦合作用载荷的实现方式为同时在所对应的节点自由度上施加单位载荷。
进一步地，所述步骤2.3.5中施加的耦合作用载荷的实现方式为在所对应的节点自由度上施加单位载荷的同时，依次分别在/>所对应的节点自由度上施加单位载荷，以及在/>所对应的节点自由度上施加单位载荷的同时，依次分别在所对应的节点自由度上施加单位载荷。
进一步地，所述步骤3中计算结构的等效板壳刚度矩阵：该刚度矩阵的求解通过对柔度矩阵/>求逆实现，柔度矩阵描述板壳结构8个广义应变与8个广义应力之间的关系，见公式（3）。/>
（3）。
其中，为柔度矩阵的待求柔度系数。
进一步地，所述步骤3.1中柔度矩阵左上角的6*6阶矩阵柔度系数计算，实现方式的根据在于如式（3）所示的广义应力应变关系式，具体地，待求系数的数值分别等于步骤2.1.2的有限元仿真中在第/>个虚拟节点自由度上施加单位载荷时第/>个虚拟节点自由度的位移计算结果。
进一步地，所述步骤3.2中柔度矩阵对角线上的剪切柔度系数计算，实现方式的根据在于宏细观应变能守恒，具体地，记步骤2.2.3中在第4、5个虚拟节点自由度上施加单位载荷时有限元仿真所得第4、5个虚拟节点自由度的位移计算结果为，记步骤2.3.3中在第7、8个虚拟节点自由度上施加单位载荷时有限元仿真所得结构单胞应变能为，则对应待求系数/>的值分别为公式（4）所示。
（4）。
进一步地，所述步骤3.3中柔度矩阵剪切耦合项柔度系数计算的特征在于：实现方式的根据在于宏细观应变能守恒，具体地，记步骤2.3.4有限元仿真所得结构单胞应变能为，则对应待求系数/>的值为公式（5）所示。
（5）。
进一步地，所述步骤3.4中柔度矩阵其余系数计算的特征在于：实现方式的根据在于宏细观应变能守恒，具体地，记步骤2.3.5中，第个虚拟节点自由度与第个虚拟节点自由度同时施加单位载荷时有限元仿真所得结构单胞应变能分别为/>，则对应待求系数/>的值为公式（6）所示。/>
（6）。
进一步地，所述步骤5板壳结构三维单胞细节应力恢复的特征在于：根据场的线性叠加方法，记步骤2.1.2、步骤2.3.3中在第个虚拟节点自由度上施加单位载荷时有限元仿真所得结构单胞位移场、应变场、应力场分别为/>，记步骤4中简化板壳模型在宏观载荷下仿真得到的所关注细节处广义应力向量为/>，则可计算得到该结构细节处的三维单胞实际位移场、应变场、应力场，见公式（7）。
（7）。
其中为广义应力向量/>的8个广义应力分量。即实现了板壳结构的细节处力学响应的快速求解。
下面，以由蒙皮、框、长桁和接头组成的舱段（具体材料参数见表1）为例，对本发明进行具体说明。
。
步骤1.航空板壳结构三维单胞建模。
首先，根据板壳的内部结构、材料布局和周期性特点，选择适当单胞进行有限元建模，划分有限元网格。对单胞模型的两个周期性方向上的端面网格节点进行节点配对，如图2所示。以单胞中心点为坐标原点构建坐标系，将方向正端面的网格节点与/>方向负端面上有相同/>坐标的节点配对，相应地，/>方向正端面的网格节点与/>方向负端面上有相同/>坐标的节点配对。
其次，设置8个虚拟节点自由度，定义第个虚拟节点自由度的位移大小为/>，其数值大小分别对应板壳的8个广义应变：/>，分别代表/>方向拉伸，面内剪切，/>方向弯曲，扭转，/>方向面外剪切。
再次，设置额外参考点进行刚体位移约束，将参考点与板壳单胞中过坐标原点的垂直于板壳中性面的一条线施加分布耦合约束，同时令该参考点的三个方向平动及三个转动为0。
步骤2.单胞多工况变形仿真。
对多种典型工况下的结构单胞响应进行仿真计算。
首先，通过如公式（1）所示的方程约束在配对的节点上施加周期性边界条件。同时在所对应的节点自由度上施加单位载荷，并进行有限元仿真，计算包含拉、弯、扭、面内剪切在内的6个单位广义应力作用工况下的板壳单胞的应力、应变响应。读取仿真结果：得到在第/>个虚拟节点自由度上施加单位载荷时，第/>个虚拟节点自由度的位移，记为/>。得到在第/>个虚拟节点自由度上施加单位载荷时结构单胞的位移场、应变场、应力场，分别记为/>。
其次，在对应虚拟节点自由度上施加单位载荷，同时令广义应变对应的虚拟节点自由度为0，并进行有限元仿真。在/>对应虚拟节点自由度上施加单位载荷，同时令广义应变/>对应的虚拟节点自由度为0，并进行有限元仿真。计算扰动场的包含额外约束的广义弯矩作用下的2个工况下的板壳单胞的应力、应变响应。读取仿真结果：得到单胞在/>对应工况下/>方向正端面上的节点在方向上的位移，记为扰动位移场/>。得到的单胞在/>对应工况下/>方向正端面上的节点在/>方向上的面力，记为/>。得到单胞在第4、5个虚拟节点自由度上施加单位载荷时，第4、5个虚拟节点自由度的位移，记为/>。
再次，通过如公式（2）所示的方程约束在配对的节点上施加包含扰动场修正的周期性边界条件。同时，在单胞方向正、负端面节点的/>方向上施加大小为/>的扰动面力场。
在所对应的节点自由度上分别单独施加单位载荷，并进行有限元仿真，计算2个面外剪切单独作用工况下的单胞的力学响应。读取仿真结果：得到单胞在第7、8个虚拟节点自由度上施加单位载荷时的应变能，记为/>。
其次，在所对应的节点自由度上同时施加单位载荷，并进行有限元仿真，计算2个面外剪切耦合作用工况下的单胞的力学响应。读取仿真结果，得到单胞应变能，记为。/>
再次，在所对应的节点自由度上施加单位载荷的同时，依次分别在所对应的节点自由度上施加单位载荷，并进行有限元仿真，计算/>方向面外剪切与其他广义应力耦合作用工况下的板壳单胞的力学响应。再次，在/>所对应的节点自由度上施加单位载荷的同时，依次分别在/>所对应的节点自由度上施加单位载荷，并进行有限元仿真，计算/>方向面外剪切与其他广义应力耦合作用工况下的板壳单胞的力学响应。读取仿真结果，分别得到单胞在第/>个虚拟节点自由度与第个虚拟节点自由度同时施加单位载荷时有限元仿真所得结构单胞应变能，记为/>。
步骤3.等效板壳属性计算。
计算结构等效板壳刚度矩阵。
根据在于如式（3）所示的广义应力应变关系式，得到柔度矩阵左上角的6×6阶薄板柔度矩阵系数为：。
根据宏细观应变能守恒，减去弯曲引起的应变能部分，根据公式（4）计算得到柔度矩阵中对角线上的两个剪切柔度系数。根据公式（5）计算得到两个方向剪切耦合项系数。根据公式（6）计算得到其余耦合项系数，即得到全部柔度系数。
最后对柔度矩阵求逆，得到板壳结构的等效刚度矩阵/>。
步骤4.板壳简化模型仿真。
基于所得到的等效板壳刚度矩阵，应用有限元仿真软件中的板壳单元建立原三维结构的低自由度简化板壳模型，施加待分析的结构服役过程中实际工况下的宏观载荷，进行有限元仿真，得到整体结构在所关注细节处的广义应力向量，记为。
步骤5.板壳结构单胞细节应力恢复。
基于场的线性叠加原理，利用上述步骤得到有限元仿真结果，根据公式（7）计算得到结构细节处的三维单胞实际位移场、应变场、应力场，实现大型复杂板壳结构的细节处力学响应的快速求解。计算细节应力场的方法包括但不限于所述的线性叠加方法，亦可将广义单元应力以边界条件形式施加于单胞重新进行有限元计算的求解。
至此基于三维单胞的航空板壳结构降阶分析方法应用完毕。
建立舱段结构的单胞有限元模型，按照上述步骤，求解舱段结构的等效刚度，并应用具有该等效板壳属性的由48个壳单元组成的降阶模型来模拟该整体结构在飞机上的实际宏观受载情况。该舱段在重力、内压、升力以及来自机翼和尾翼的弯矩作用下处于平衡状态，如图5所示。通过有限元仿真计算机舱简化模型的挠度，结果如图6所示。同时为验证方法的正确性，采用1790352个单元的全尺寸有限元模型进行相同工况下的仿真。舱段挠度计算结果对比如图7所示，最大相对误差仅为不到5%。
实施例3。
为了执行上述实施例1对应的方法，以实现相应的功能和技术效果，下面提供了一种航空舱段结构降阶分析系统，包括：
当前单胞结构确定模块，用于确定航空舱段结构中任一单胞结构为当前单胞结构。航空舱段结构为筒状结构。航空舱段结构是将多个单胞结构按照周期性排列规则进行排列的得到的。周期性排列规则为单胞结构沿航空舱段的端面圆周呈周期性排列。单胞结构沿航空舱段的径向呈周期性排列。
单胞模型构建模块，用于利用有限元建模软件对当前单胞结构进行有限元建模处理，得到单胞模型。
单胞仿真模块，用于基于工况种类调整单胞模型的周期性边界条件对单胞模型进行仿真，得到单胞仿真结果。单胞仿真结果包括在不同虚拟节点自由度上施加单位载荷时单胞结构的位移场、应变场和应力场。
等效刚度矩阵确定模块，用于基于单胞仿真结果，确定当前单胞结构的等效刚度矩阵。
整体简化模型构建模块，用于基于等效刚度矩阵和周期性排列规则，利用有限元建模软件构建航空舱段结构的整体简化模型。
关注细节区域获取模块，用于获取航空舱段结构中的关注细节区域。关注细节区域的面积小于单胞结构的面积。
广义应力向量确定模块，用于以实际工况下的宏观载荷为条件，对整体简化模型进行有限元仿真，得到航空舱段结构中关注细节区域处的广义应力向量。
细节力学响应确定模块，用于根据单胞仿真结果和广义应力向量，确定航空舱段结构中关注细节区域处的力学响应。力学响应包括位移场、应变场和应力场。
实施例4。
本实施例提供了一种电子设备，包括存储器及处理器，存储器用于存储计算机程序，处理器运行计算机程序以使电子设备执行实施例1或实施例2所述的一种航空舱段结构降阶分析方法。存储器为可读存储介质。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。
Claims (8)
1.一种航空舱段结构降阶分析方法，其特征在于，包括：
确定航空舱段结构中任一单胞结构为当前单胞结构；所述航空舱段结构为筒状结构；所述航空舱段结构是将多个所述单胞结构按照周期性排列规则进行排列的得到的；所述周期性排列规则为所述单胞结构沿所述航空舱段的端面圆周呈周期性排列；所述单胞结构沿所述航空舱段的径向呈周期性排列；
利用有限元建模软件对当前单胞结构进行有限元建模处理，得到单胞模型；
基于工况种类调整单胞模型的周期性边界条件对单胞模型进行仿真，得到单胞仿真结果；所述单胞仿真结果包括在不同虚拟节点自由度上施加单位载荷时单胞结构的位移场、应变场和应力场；
基于所述单胞仿真结果，确定当前单胞结构的等效刚度矩阵；
基于所述等效刚度矩阵和所述周期性排列规则，利用有限元建模软件构建航空舱段结构的整体简化模型；
获取航空舱段结构中的关注细节区域；所述关注细节区域的面积小于所述单胞结构的面积；
以实际工况下的宏观载荷为条件，对整体简化模型进行有限元仿真，得到航空舱段结构中关注细节区域处的广义应力向量；
根据所述单胞仿真结果和所述广义应力向量，确定航空舱段结构中关注细节区域处的力学响应；所述力学响应包括位移场、应变场和应力场；
基于关注细节区域处的力学响应，优化所述航空舱段结构；
利用有限元建模软件对当前单胞结构进行有限元建模处理，得到单胞模型，包括：
利用有限元建模软件对当前单胞结构进行网格划分，得到多个网格节点；
以当前单胞结构中心点为原点O，以第一周期性方向为方向，以第二周期性方向为/>方向，以垂直于/>平面的方向为/>方向，构建三维坐标系；
确定第一正端面和第一负端面；所述第一正端面和所述第一负端面在所述第一周期性方向上平行设置；
将所述第一正端面和所述第一负端面上方向坐标相等且/>方向坐标相等的网格节点进行配对处理；
确定第二正端面和第二负端面；所述第二正端面和所述第二负端面在所述第二周期性方向上平行设置；
将所述第二正端面和所述第二负端面上方向坐标相等且/>方向坐标相等的网格节点进行配对处理；
基于板壳广义应变类型设置单胞模型中的虚拟节点自由度；所述虚拟节点自由度与板壳广义应变类型一一对应；所述广义应变包括方向拉伸、/>方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲、扭转、/>方向面外剪切和/>方向面外剪切；
获取参考点和参考直线；所述参考直线为当前单胞结构中过坐标原点且垂直于板壳中性面的任一直线；
基于单胞刚体位移约束条件，设置参考点和参考直线处的参数，约束当前单胞结构的刚体位移；所述单胞刚体位移约束条件为设置参考点处方向平动、/>方向平动、以/>方向平动、/>方向转动、/>方向转动和以/>方向转动均为0，并向参考直线施加分布耦合约束；
基于工况种类调整单胞模型的周期性边界条件对单胞模型进行仿真，得到单胞仿真结果，包括：
基于周期性边界条件，分别在第一工况下和第二工况下对单胞模型进行仿真，得到广义应力作用下的单胞仿真结果；
基于周期性边界条件，分别在第三工况下和第四工况下对单胞模型进行仿真，得到广义弯矩作用下的单胞仿真结果；
根据广义弯矩作用下的单胞仿真结果确定扰动面力场施加量；
基于扰动面力场施加量和修正周期性边界条件，分别在第五工况下、第六工况下、第七工况下、第八工况下和第九工况下对单胞模型进行仿真，得到广义应力耦合作用下的单胞仿真结果。
2.根据权利要求1所述的一种航空舱段结构降阶分析方法，其特征在于，所述航空舱段结构具有周期性特点。
3.根据权利要求1所述的一种航空舱段结构降阶分析方法，其特征在于，
所述第一工况为在方向拉伸、方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲和扭转对应的虚拟节点自由度上均同时施加单位载荷；
所述第二工况在所有虚拟节点自由度上分别施加单位载荷；
所述第三工况为在方向弯曲对应的虚拟节点自由度上施加单位载荷，同时在/>方向拉伸、/>方向弯曲、扭转和/>方向面外剪切对应的虚拟节点自由度上均施加0载荷；
所述第四工况为在方向弯曲对应的虚拟节点自由度上施加单位载荷，同时在/>方向拉伸、/>方向弯曲、扭转和方向面外剪切对应的虚拟节点自由度上均施加0载荷；
所述第五工况为在方向面外剪切对应的虚拟节点自由度上施加单位载荷；
所述第六工况为在方向面外剪切对应的虚拟节点自由度上施加单位载荷；
所述第七工况为在方向面外剪切和/>方向面外剪切对应的虚拟节点自由度上同时施加单位载荷；
所述第八工况为在方向面外剪切对应的虚拟节点自由度上施加单位载荷，依次分别在/>方向拉伸、/>方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲和扭转对应的虚拟节点自由度上施加单位载荷；
所述第九工况为在方向面外剪切对应的虚拟节点自由度上施加单位载荷，依次分别在/>方向拉伸、/>方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲和扭转对应的虚拟节点自由度上施加单位载荷。
4.根据权利要求1所述的一种航空舱段结构降阶分析方法，其特征在于，
所述广义应力作用下的单胞仿真结果包括：第一工况下6个虚拟节点自由度的位移，以及第二工况下8个虚拟节点自由度的力学响应；
所述广义弯矩作用下的单胞仿真结果包括：第三工况下单胞扰动位移场、单胞面力和方向弯曲对应的虚拟节点自由度上的位移，以及第四工况下单胞扰动位移场、单胞面力和/>方向弯曲对应的虚拟节点自由度上的位移；
所述广义应力耦合作用下的单胞仿真结果包括第五工况下方向面外剪切对应的虚拟节点自由度的应变能、第六工况下/>方向面外剪切对应的虚拟节点自由度的应变能、第七工况下的/>方向面外剪切和/>方向面外剪切相互作用时的单胞应变能、第八工况下/>方向面外剪切与/>方向拉伸、/>方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲和扭转中任一广义应变类型相互作用时的单胞应变能和第九工况下/>方向面外剪切与/>方向拉伸、方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲和扭转中任一广义应变类型相互作用时的单胞应变能。
5.根据权利要求4所述的一种航空舱段结构降阶分析方法，其特征在于，基于所述单胞仿真结果，确定当前单胞结构的等效刚度矩阵，包括：
根据广义应变-广义应力关系，构建初始柔度矩阵；所述初始柔度矩阵为8阶矩阵；
确定所述初始柔度矩阵中的第一参数、第二参数、第三参数和第四参数；所述第一参数为初始柔度矩阵中的元素，其中/>；所述第二参数为初始柔度矩阵中的元素/>，其中/>；所述第三参数为初始柔度矩阵中的元素/>，其中；所述第四参数为初始柔度矩阵中除所述第一参数、所述第二参数和所述第三参数之外的参数；
基于所述广义应力作用下的单胞仿真结果，确定所述初始柔度矩阵中的第一参数；
基于所述广义弯矩作用下的单胞仿真结果、所述第五工况下方向面外剪切对应的虚拟节点自由度的应变能和所述第六工况下/>方向面外剪切对应的虚拟节点自由度的应变能，确定所述初始柔度矩阵中的第二参数；
基于所述第七工况下的方向面外剪切和/>方向面外剪切相互作用时的单胞应变能，确定所述初始柔度矩阵中的第三参数；
基于第八工况下方向面外剪切与/>方向拉伸、/>方向拉伸、面内剪切、/>方向弯曲、方向弯曲和扭转中任一广义应变类型相互作用时的单胞应变能，以及第九工况下/>方向面外剪切与/>方向拉伸、/>方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲和扭转中任一广义应变类型相互作用时的单胞应变能，确定所述初始柔度矩阵中的第四参数，得到柔度矩阵；
确定所述柔度矩阵的逆矩阵为当前单胞结构的等效刚度矩阵。
6.根据权利要求5所述的一种航空舱段结构降阶分析方法，其特征在于，根据所述单胞仿真结果和所述广义应力向量，确定航空舱段结构中关注细节区域处的力学响应，包括：
确定任一广义应力类型为当前广义应力类型；
确定当前广义应力类型对应的广义应力向量分量，与确定第二工况下当前广义应力类型对应虚拟节点自由度的力学响应之积，为航空舱段结构中关注细节区域处在当前广义应力类型对应虚拟节点自由度上的力学响应；所述广义应力向量包括多个广义应力向量分量；所述广义应力向量分量与所述广义应力类型一一对应。
7.一种航空舱段结构降阶分析系统，其特征在于，包括：
当前单胞结构确定模块，用于确定航空舱段结构中任一单胞结构为当前单胞结构；所述航空舱段结构为筒状结构；所述航空舱段结构是将多个所述单胞结构按照周期性排列规则进行排列的得到的；所述周期性排列规则为所述单胞结构沿所述航空舱段的端面圆周呈周期性排列；所述单胞结构沿所述航空舱段的径向呈周期性排列；
单胞模型构建模块，用于利用有限元建模软件对当前单胞结构进行有限元建模处理，得到单胞模型；利用有限元建模软件对当前单胞结构进行网格划分，得到多个网格节点；以当前单胞结构中心点为原点O，以第一周期性方向为方向，以第二周期性方向为/>方向，以垂直于/>平面的方向为/>方向，构建三维坐标系；确定第一正端面和第一负端面；所述第一正端面和所述第一负端面在所述第一周期性方向上平行设置；将所述第一正端面和所述第一负端面上/>方向坐标相等且/>方向坐标相等的网格节点进行配对处理；确定第二正端面和第二负端面；所述第二正端面和所述第二负端面在所述第二周期性方向上平行设置；将所述第二正端面和所述第二负端面上/>方向坐标相等且/>方向坐标相等的网格节点进行配对处理；基于板壳广义应变类型设置单胞模型中的虚拟节点自由度；所述虚拟节点自由度与板壳广义应变类型一一对应；所述广义应变包括/>方向拉伸、/>方向拉伸、面内剪切、/>方向弯曲、/>方向弯曲、扭转、/>方向面外剪切和/>方向面外剪切；获取参考点和参考直线；所述参考直线为当前单胞结构中过坐标原点且垂直于板壳中性面的任一直线；基于单胞刚体位移约束条件，设置参考点和参考直线处的参数，约束当前单胞结构的刚体位移；所述单胞刚体位移约束条件为设置参考点处/>方向平动、/>方向平动、以/>方向平动、/>方向转动、/>方向转动和以/>方向转动均为0，并向参考直线施加分布耦合约束；
单胞仿真模块，用于基于工况种类调整单胞模型的周期性边界条件对单胞模型进行仿真，得到单胞仿真结果；所述单胞仿真结果包括在不同虚拟节点自由度上施加单位载荷时单胞结构的位移场、应变场和应力场；基于周期性边界条件，分别在第一工况下和第二工况下对单胞模型进行仿真，得到广义应力作用下的单胞仿真结果；基于周期性边界条件，分别在第三工况下和第四工况下对单胞模型进行仿真，得到广义弯矩作用下的单胞仿真结果；根据广义弯矩作用下的单胞仿真结果确定扰动面力场施加量；基于扰动面力场施加量和修正周期性边界条件，分别在第五工况下、第六工况下、第七工况下、第八工况下和第九工况下对单胞模型进行仿真，得到广义应力耦合作用下的单胞仿真结果；
等效刚度矩阵确定模块，用于基于所述单胞仿真结果，确定当前单胞结构的等效刚度矩阵；
整体简化模型构建模块，用于基于所述等效刚度矩阵和所述周期性排列规则，利用有限元建模软件构建航空舱段结构的整体简化模型；
关注细节区域获取模块，用于获取航空舱段结构中的关注细节区域；所述关注细节区域的面积小于所述单胞结构的面积；
广义应力向量确定模块，用于以实际工况下的宏观载荷为条件，对整体简化模型进行有限元仿真，得到航空舱段结构中关注细节区域处的广义应力向量；
细节力学响应确定模块，用于根据所述单胞仿真结果和所述广义应力向量，确定航空舱段结构中关注细节区域处的力学响应；所述力学响应包括位移场、应变场和应力场；
航空舱段结构优化模块，用于基于关注细节区域处的力学响应，优化所述航空舱段结构。
8.一种电子设备，其特征在于，包括存储器及处理器，所述存储器用于存储计算机程序，所述处理器运行所述计算机程序以使所述电子设备执行权利要求1至6中任一项所述的一种航空舱段结构降阶分析方法；所述存储器为可读存储介质。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410032291.6A CN117540494B (zh) | 2024-01-10 | 2024-01-10 | 一种航空舱段结构降阶分析方法、系统及设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410032291.6A CN117540494B (zh) | 2024-01-10 | 2024-01-10 | 一种航空舱段结构降阶分析方法、系统及设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117540494A CN117540494A (zh) | 2024-02-09 |
| CN117540494B true CN117540494B (zh) | 2024-04-05 |
Family
ID=89790414
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410032291.6A Active CN117540494B (zh) | 2024-01-10 | 2024-01-10 | 一种航空舱段结构降阶分析方法、系统及设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117540494B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112560167A (zh) * | 2020-11-10 | 2021-03-26 | 北京航空航天大学 | 机翼结构力学高保真降阶仿真方法、电子设备、存储介质 |
| CN116305589A (zh) * | 2023-05-22 | 2023-06-23 | 北京航空航天大学 | 一种直升机桨叶结构降阶分析方法、系统、设备及介质 |
| WO2023142333A1 (zh) * | 2022-01-25 | 2023-08-03 | 大连理工大学 | 一种薄壁筒壳模型修正方法及系统 |
| WO2023151233A1 (zh) * | 2022-02-09 | 2023-08-17 | 南京航空航天大学 | 基于虚拟应变能的金属材料多轴疲劳寿命预测方法与系统 |
| CN116738707A (zh) * | 2023-06-08 | 2023-09-12 | 西安交通大学 | 部分周期性换热器通道的等效力学性能预测方法及系统 |
| CN116822167A (zh) * | 2023-06-08 | 2023-09-29 | 西安交通大学 | 换热器热力耦合性能多尺度分析方法、系统、介质及设备 |
-
2024
- 2024-01-10 CN CN202410032291.6A patent/CN117540494B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112560167A (zh) * | 2020-11-10 | 2021-03-26 | 北京航空航天大学 | 机翼结构力学高保真降阶仿真方法、电子设备、存储介质 |
| WO2023142333A1 (zh) * | 2022-01-25 | 2023-08-03 | 大连理工大学 | 一种薄壁筒壳模型修正方法及系统 |
| WO2023151233A1 (zh) * | 2022-02-09 | 2023-08-17 | 南京航空航天大学 | 基于虚拟应变能的金属材料多轴疲劳寿命预测方法与系统 |
| CN116305589A (zh) * | 2023-05-22 | 2023-06-23 | 北京航空航天大学 | 一种直升机桨叶结构降阶分析方法、系统、设备及介质 |
| CN116738707A (zh) * | 2023-06-08 | 2023-09-12 | 西安交通大学 | 部分周期性换热器通道的等效力学性能预测方法及系统 |
| CN116822167A (zh) * | 2023-06-08 | 2023-09-29 | 西安交通大学 | 换热器热力耦合性能多尺度分析方法、系统、介质及设备 |
Non-Patent Citations (2)
| Title |
|---|
| 杨留义 ; 谭惠丰 ; 曹宗胜 ; .基于单胞有限元模型的织物复合材料弯曲刚度预报.复合材料学报.2017,(05),全文. * |
| 陈利 ; 许正亚 ; 马振杰 ; 李嘉禄 ; .基于单胞数字化模型的复合材料本构关系三维数值模拟.复合材料学报.2007,(03),全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117540494A (zh) | 2024-02-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Li et al. | An element-free study of variable stiffness composite plates with cutouts for enhanced buckling and post-buckling performance | |
| Guruswamy | A review of numerical fluids/structures interface methods for computations using high-fidelity equations | |
| CN110287637B (zh) | 一种弹塑性屈曲承载力计算方法 | |
| CN107220461A (zh) | 一种变刚度复合材料板壳结构高效优化方法 | |
| CN112560167A (zh) | 机翼结构力学高保真降阶仿真方法、电子设备、存储介质 | |
| Zhuang et al. | Modal and aeroelastic analysis of trapezoidal corrugated-core sandwich panels in supersonic flow | |
| CN106096257A (zh) | 一种非线性索单元分析方法及系统 | |
| CN112580241A (zh) | 一种基于结构降阶模型的非线性气动弹性动响应分析方法 | |
| Li et al. | A peridynamic model for the nonlinear static analysis of truss and tensegrity structures | |
| CN110955941A (zh) | 基于向量场的复合材料结构优化设计方法及设备 | |
| Wang et al. | Voronoi polygonal hybrid finite elements with boundary integrals for plane isotropic elastic problems | |
| Zhang et al. | A morphing wing with cellular structure of non-uniform density | |
| Li et al. | An efficient implementation of aeroelastic tailoring based on efficient computational fluid dynamics-based reduced order model | |
| Islam et al. | An efficient model for laminated composite thin-walled beams of open or closed cross-section and with or without in-filled materials | |
| Zhang et al. | Optimization and analysis of composite sandwich box beam for solar drones | |
| Zhang et al. | Mathematical modeling and dynamic characteristic analysis of a novel parallel tracking mechanism for inter-satellite link antenna | |
| Ren et al. | An adaptive triangular element of absolute nodal coordinate formulation for thin plates and membranes | |
| CN117171864B (zh) | 一种梁结构线性振动预测方法 | |
| CN117540494B (zh) | 一种航空舱段结构降阶分析方法、系统及设备 | |
| Ge et al. | Static analysis of defective sandwich beam by Chebyshev quadrature element method | |
| Liu et al. | A new layout optimization method for stiffened panels based on ground stiffener structure (GSS) and thickness penalty | |
| Liu et al. | Buckling and post-buckling analysis of butterfly-shaped auxetic core sandwich plates based on variational asymptotic method | |
| Wang et al. | Deformation performance analysis of thin plates based on a deformation decomposition method | |
| CN115906333A (zh) | 一种桁架结构的几何非线性等效板动力学建模与响应分析方法 | |
| Werthen et al. | Comparison of different cross-sectional approaches for the structural design and optimization of composite wind turbine blades based on beam models |
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 |