CN107145676A - 一种复杂载荷条件下的结构层级式有限元仿真分析方法 - Google Patents
一种复杂载荷条件下的结构层级式有限元仿真分析方法 Download PDFInfo
- Publication number
- CN107145676A CN107145676A CN201710357712.2A CN201710357712A CN107145676A CN 107145676 A CN107145676 A CN 107145676A CN 201710357712 A CN201710357712 A CN 201710357712A CN 107145676 A CN107145676 A CN 107145676A
- Authority
- CN
- China
- Prior art keywords
- finite element
- submodel
- model
- main structure
- conditions
- 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
Links
Classifications
-
- 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
- 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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/20—Configuration CAD, e.g. designing by assembling or positioning modules selected from libraries of predesigned modules
-
- 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)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开的一种复杂载荷条件下的结构层级式有限元仿真分析方法,包括以下步骤:S1:建立舱体主结构全尺寸有限元模型作为1级母模型;S2:针对主结构失效方式进行有限元计算,计算方法涵盖结构稳定性分析和静力强度分析;S3:对主结构失稳载荷低于1.5倍设计载荷的局部失稳位置进行危险性评估,选定细化分析区域;S4:对主结构强度剩余系数低于1.3的应力危险区域进行危险性评估,选定细化分析区域;S5:以S3、S4中的细化分析区域建立有限元1级子模型;S6:建立子模型与母模型的有限元连接模型;S7:将含有过渡模型的子模型提交有限元计算;S8:重复S1~S7步骤,直至模型足够精细。
Description
技术领域
本发明涉及运载器、飞行器强度设计领域,具体涉及舱体主结构强度的数值分析方法。
背景技术
当前的航天发展形势给仿真行业提出了空前的挑战。一方面,轻量化设计使得结构设计裕度不断降低,这给仿真预试的准确性提出了严格的要求。另一方面,密集发射、密集设计又使得工作效率需要不断提升,落实在仿真行业即为不断地提升仿真建模和分析效率。对于复杂载荷、复杂布局的舱体结构,仿真方法准确性和计算效率的矛盾更为突出,因此必须研究一种可以兼顾准确性和计算效率的有限元分析方法。
传统的有限元分析常采用多模型计算解决以上矛盾。建立整体模型用于分析结构整体稳定性和主传力部件的强度失效问题。建立局部模型用于分析局部区域在承受复杂载荷下可能发生的失效问题。这种分析思想隐藏了三个隐患。首先,整体模型的简化有可能忽略结构传力路径上的失效行为。例如,部件之间的连接若采用刚性约束代替,则无法分析出连接失效形式。其次,局部模型往往与整体模型割裂。局部模型的分析区域往往由设计师凭借工程经验确定,可能与整体模型的危险区域并不一致。而且,局部模型的边界常为固支或简支边界,与整体模型并无关系。这种割裂使得局部模型分析的覆盖度和精度有一定程度的降低。
本发明公开了一种复杂载荷条件下的结构层级式有限元仿真分析方法,该方法能够准确、有效地计算复杂舱体的失效位置与承载能力,可以广泛应用于航天运载器、飞行器的舱体结构设计。
发明内容
本发明提供一种复杂载荷条件下的结构层级式有限元仿真分析方法,吸取了传统有限元的多模型分析思想,但通过识别破坏行为、引入多模型之间的边界关系,建立多模型之间的联系,避免传统多模型分析覆盖性和精度降低的隐患。这种分析方法的计算效率和计算精度已得到多次的静力试验验证。
本发明通过以下技术方案实现:
一种复杂载荷条件下的结构层级式有限元仿真分析方法,包括以下步骤:
S1:建立舱体主结构全尺寸有限元模型作为1级母模型;
S2:针对主结构失效方式进行有限元计算,计算载荷为设计载荷,计算方法涵盖结构稳定性分析和静力强度分析;
S3:进行主结构稳定性分析后处理,分析失稳载荷低于1.5倍设计载荷的局部失稳位置,进行危险性评估后选定细化分析区域;
S4:进行主结构静力强度分析后处理,分析强度剩余系数低于1.3的应力危险区域,进行危险性评估后选定细化分析区域;
S5:以S3、S4中的细化分析区域建立有限元1级子模型,视模型进行单元简化;
S6:建立子模型与母模型的有限元过渡模型,该模型可以真实地传递母模型内部边界区域的位移场和应力场;
S7:将含有过渡模型的子模型提交有限元计算;
S8:如果1级子模型仍不能详细反映结构破坏方式,可以1级子模型作为2级母模型,重复S1~S6步骤,直至模型足够精细。
较佳的,步骤S1中蒙皮、桁条、框、加强件均简化为板单元,铆钉或螺栓简化为梁单元。梁单元与板单元之间采用共节点或局部多点约束连接。
较佳的,步骤S5中,蒙皮仍简化为板单元,桁条、框、加强件简化为体单元或板单元,铆钉或螺栓简化为体单元或梁单元。各部件之间采用最为真实的连接模拟方式。
较佳的,步骤S5中的连接模拟方式为三维有限元模型的接触关系,包括法向与切向的接触关系。
较佳的,步骤S6中的连接模型细化为:单纯的子模型,其外部边界的场变量由母模型内部边界提供;子模型、母模型混合模型,子模型外部边界与母模型内部边界通过共节点、多点约束、壳与实体耦合约束等方式连接。
附图说明
图1为本发明一舱体结构示意图;
图2为本发明一舱体结构有限元分析1级母模型示意图;
图3为本发明一舱体结构有限元分析1级子模型示意图;
图4为本发明一舱体结构有限元分析1级子模型失效机制图;
图5为本发明提供一种复杂载荷条件下的结构层级式有限元仿真分析方法流程图。
具体实施方式
下面结合实施例对本发明作详细说明,本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式,但本发明的保护范围不限于下述的实施例。
如图1所示,为本发明实施例的一个舱体主结构示意图,结构在接头处承受集中载荷。
本发明提供的复杂载荷条件下的结构层级式有限元仿真分析方法,较为详细的,包括以下步骤:
步骤1:建立舱体结构全尺寸有限元模型作为母模型。如图2右侧,蒙皮简化为弱刚性板单元,桁条、框、接头简化为板单元;如图2左侧,铆钉或螺栓简化为梁单元。梁单元与板单元之间采用共节点约束;
步骤2:针对舱体结构全尺寸模型进行有限元计算,选用软件为ABAQUS,选用计算步为Static,General/Static,Riks,前者用于计算结构强度,后者用于计算结构稳定性,设计载荷为3000kN;
步骤3:根据稳定性分析结果,失稳载荷超过1.5倍设计载荷。根据强度分析结果,应力危险区域在接头附近,表现为接头和接头下方的连接件破坏。综合分析结果,选取接头附近作为细化分析区域,如图3所示;
步骤4:将母模型的接头附近区域删除,替换为子模型,即接头、框、桁、力扩散板、螺栓、铆钉均简化为体单元,蒙皮简化为板单元。各体单元之间均采用接触连接,蒙皮与其他部件之间采用刚性约束连接;
步骤5:将子模型、母模型之间的边界均用刚性约束连接,如图3所示;
步骤6:将新的修正模型提交静力有限元计算,仍选用Static,General计算步,接头处仍施加3000kN设计载荷;
步骤7:将修正模型进行有限元后处理,可以辨认出接头、连接件各剖面的应力云图,如图4所示。
该方法通过母模型分析精确识别稳定性破坏模式和主结构强度破坏模式,并大致判断出连接件的破坏模式,进一步通过子模型精确识别部件和连接件的破坏行为。通过子模型、母模型的边界耦合,使得分析规模大幅度降低,并保证了子模型的边界准确性。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
Claims (5)
1.一种复杂载荷条件下的结构层级式有限元仿真分析方法,其特征在于,包括以下步骤:
S1:建立舱体主结构全尺寸有限元模型作为1级母模型;
S2:针对主结构失效方式进行有限元计算,计算载荷为设计载荷,计算方法涵盖结构稳定性分析和静力强度分析;
S3:进行主结构稳定性分析后处理,分析失稳载荷低于1.5倍设计载荷的局部失稳位置,进行危险性评估后选定细化分析区域;
S4:进行主结构静力强度分析后处理,分析强度剩余系数低于1.3的应力危险区域,进行危险性评估后选定细化分析区域;
S5:以S3、S4中的细化分析区域建立有限元1级子模型,视模型进行单元简化;
S6:建立子模型与母模型的有限元连接模型,该模型可以真实地传递母模型内部边界区域的位移场和应力场;
S7:将含有过渡模型的子模型提交有限元计算;
S8:如果1级子模型仍不能详细反映结构破坏方式,可以1级子模型作为2级母模型,重复S1~S7步骤,直至模型足够精细。
2.如权利要求1所述的复杂载荷条件下的结构层级式有限元仿真分析方法,其特征在于,步骤S1中蒙皮、桁条、框、加强件均简化为板单元,铆钉或螺栓简化为梁单元;梁单元与板单元之间采用共节点或局部多点约束连接。
3.如权利要求1所述的复杂载荷条件下的结构层级式有限元仿真分析方法,其特征在于,步骤S5中蒙皮简化为板单元,桁条、框、加强件简化为体单元或板单元,铆钉或螺栓简化为体单元或梁单元;各部件之间采用最为真实的连接模拟方式。
4.如权利要求3所述的复杂载荷条件下的结构层级式有限元仿真分析方法,其特征在于,步骤S5中所述的连接模拟方式为三维有限元模型的接触关系,包括法向与切向的接触关系。
5.如权利要求1所述的复杂载荷条件下的结构层级式有限元仿真分析方法,其特征在于,步骤S6中所述的连接模型细化为:单纯的子模型,其外部边界的场变量由母模型内部边界提供;子模型、母模型混合模型,子模型外部边界与母模型内部边界通过共节点、多点约束、壳与实体耦合约束等方式连接。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710357712.2A CN107145676B (zh) | 2017-05-19 | 2017-05-19 | 一种复杂载荷条件下的结构层级式有限元仿真分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710357712.2A CN107145676B (zh) | 2017-05-19 | 2017-05-19 | 一种复杂载荷条件下的结构层级式有限元仿真分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107145676A true CN107145676A (zh) | 2017-09-08 |
CN107145676B CN107145676B (zh) | 2020-07-14 |
Family
ID=59777053
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710357712.2A Active CN107145676B (zh) | 2017-05-19 | 2017-05-19 | 一种复杂载荷条件下的结构层级式有限元仿真分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107145676B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108763603A (zh) * | 2018-01-12 | 2018-11-06 | 温州燃机发电有限公司 | 一种采用非线性有限元整体—局部联合建模进行汽门紧固螺栓耦合应力分析的方法 |
CN109684675A (zh) * | 2018-12-04 | 2019-04-26 | 中国航空工业集团公司西安飞机设计研究所 | 一种基于abaqus板体耦合的航空薄壁支架强度分析方法 |
CN109710982A (zh) * | 2018-12-04 | 2019-05-03 | 中国特种飞行器研究所 | 基于有限元模拟分析的部件级模拟试验件设计验证方法 |
CN109753687A (zh) * | 2018-12-04 | 2019-05-14 | 中国航空工业集团公司西安飞机设计研究所 | 一种复合式气动弹性动力学建模方法 |
CN110737997A (zh) * | 2019-08-05 | 2020-01-31 | 中国船舶科学研究中心(中国船舶重工集团公司第七0二研究所) | 复杂载荷作用下超大型浮体结构极限强度的分析方法 |
CN111651881A (zh) * | 2020-06-01 | 2020-09-11 | 中国第一汽车股份有限公司 | 一种锁铆仿真失效参数的简化方法 |
CN112560311A (zh) * | 2020-12-17 | 2021-03-26 | 中航沈飞民用飞机有限责任公司 | 一种快速自动提取结构信息并对结构进行分析和优化的方法 |
CN112818479A (zh) * | 2021-01-15 | 2021-05-18 | 北京理工大学 | 一种有限元子模型边界条件加载方法 |
CN113239592A (zh) * | 2021-05-19 | 2021-08-10 | 北京航空航天大学 | 一种面向复杂多层次结构产品仿真应用的层级式载荷分解方法 |
CN114626235A (zh) * | 2022-03-21 | 2022-06-14 | 北京航空航天大学 | 一种面向复杂多层级结构产品的载荷分解改进方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7430499B2 (en) * | 2005-03-23 | 2008-09-30 | The Boeing Company | Methods and systems for reducing finite element simulation time for acoustic response analysis |
CN103452495A (zh) * | 2013-09-11 | 2013-12-18 | 上海大学 | 一种确定复杂工况条件下钻具接头极限工作扭矩的方法 |
CN104408238A (zh) * | 2014-11-11 | 2015-03-11 | 上海宇航系统工程研究所 | 一种桁式半硬壳结构轴压稳定性的数值分析方法 |
CN106055782A (zh) * | 2016-05-27 | 2016-10-26 | 中铁第四勘察设计院集团有限公司 | 一种隧道施工救生舱舱体强度有限元模拟分析方法 |
CN106156449A (zh) * | 2016-08-31 | 2016-11-23 | 中航沈飞民用飞机有限责任公司 | 一种复合材料机翼壁板优化设计方法 |
-
2017
- 2017-05-19 CN CN201710357712.2A patent/CN107145676B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7430499B2 (en) * | 2005-03-23 | 2008-09-30 | The Boeing Company | Methods and systems for reducing finite element simulation time for acoustic response analysis |
CN103452495A (zh) * | 2013-09-11 | 2013-12-18 | 上海大学 | 一种确定复杂工况条件下钻具接头极限工作扭矩的方法 |
CN104408238A (zh) * | 2014-11-11 | 2015-03-11 | 上海宇航系统工程研究所 | 一种桁式半硬壳结构轴压稳定性的数值分析方法 |
CN106055782A (zh) * | 2016-05-27 | 2016-10-26 | 中铁第四勘察设计院集团有限公司 | 一种隧道施工救生舱舱体强度有限元模拟分析方法 |
CN106156449A (zh) * | 2016-08-31 | 2016-11-23 | 中航沈飞民用飞机有限责任公司 | 一种复合材料机翼壁板优化设计方法 |
Non-Patent Citations (3)
Title |
---|
刘松林: "基于Gurson模型的压力连接过程仿真及连接点强度预测", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑 (月刊 )》 * |
曾清香: "飞行器折叠翼设计及其运动动态性能分析与试验研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑 (月刊 )》 * |
王云开 等: "基于有限元仿真计算的低压环境舱强度分析", 《科技资讯》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108763603A (zh) * | 2018-01-12 | 2018-11-06 | 温州燃机发电有限公司 | 一种采用非线性有限元整体—局部联合建模进行汽门紧固螺栓耦合应力分析的方法 |
CN108763603B (zh) * | 2018-01-12 | 2022-06-14 | 温州燃机发电有限公司 | 一种汽门紧固螺栓耦合应力分析的方法 |
CN109684675A (zh) * | 2018-12-04 | 2019-04-26 | 中国航空工业集团公司西安飞机设计研究所 | 一种基于abaqus板体耦合的航空薄壁支架强度分析方法 |
CN109710982A (zh) * | 2018-12-04 | 2019-05-03 | 中国特种飞行器研究所 | 基于有限元模拟分析的部件级模拟试验件设计验证方法 |
CN109753687A (zh) * | 2018-12-04 | 2019-05-14 | 中国航空工业集团公司西安飞机设计研究所 | 一种复合式气动弹性动力学建模方法 |
CN109710982B (zh) * | 2018-12-04 | 2023-05-26 | 中国特种飞行器研究所 | 基于有限元模拟分析的部件级模拟试验件设计验证方法 |
CN109684675B (zh) * | 2018-12-04 | 2023-04-18 | 中国航空工业集团公司西安飞机设计研究所 | 一种基于abaqus板体耦合的航空薄壁支架强度分析方法 |
CN110737997A (zh) * | 2019-08-05 | 2020-01-31 | 中国船舶科学研究中心(中国船舶重工集团公司第七0二研究所) | 复杂载荷作用下超大型浮体结构极限强度的分析方法 |
CN110737997B (zh) * | 2019-08-05 | 2024-01-26 | 中国船舶科学研究中心(中国船舶重工集团公司第七0二研究所) | 复杂载荷作用下超大型浮体结构极限强度的分析方法 |
CN111651881A (zh) * | 2020-06-01 | 2020-09-11 | 中国第一汽车股份有限公司 | 一种锁铆仿真失效参数的简化方法 |
CN112560311A (zh) * | 2020-12-17 | 2021-03-26 | 中航沈飞民用飞机有限责任公司 | 一种快速自动提取结构信息并对结构进行分析和优化的方法 |
CN112818479B (zh) * | 2021-01-15 | 2022-09-02 | 北京理工大学 | 一种有限元子模型边界条件加载方法 |
CN112818479A (zh) * | 2021-01-15 | 2021-05-18 | 北京理工大学 | 一种有限元子模型边界条件加载方法 |
CN113239592B (zh) * | 2021-05-19 | 2022-07-01 | 北京航空航天大学 | 一种面向复杂多层次结构产品仿真应用的层级式载荷分解方法 |
CN113239592A (zh) * | 2021-05-19 | 2021-08-10 | 北京航空航天大学 | 一种面向复杂多层次结构产品仿真应用的层级式载荷分解方法 |
CN114626235A (zh) * | 2022-03-21 | 2022-06-14 | 北京航空航天大学 | 一种面向复杂多层级结构产品的载荷分解改进方法 |
CN114626235B (zh) * | 2022-03-21 | 2024-06-07 | 北京航空航天大学 | 一种面向复杂多层级结构产品的载荷分解改进方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107145676B (zh) | 2020-07-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107145676A (zh) | 一种复杂载荷条件下的结构层级式有限元仿真分析方法 | |
CN102456082B (zh) | 二维-三维相结合的卫星总装工艺规划与仿真方法 | |
Meneghetti et al. | The peak stress method to estimate the mode I notch stress intensity factor in welded joints using three-dimensional finite element models | |
CN106123771B (zh) | 基于激光位移传感器的轮对尺寸测量结果可信度的区间模糊评判方法 | |
CN103063401B (zh) | 一种钢桁架拱桥风洞试验模型的制备方法 | |
Lin et al. | Automating closed-loop structural safety management for bridge construction through multisource data integration | |
CN109710982A (zh) | 基于有限元模拟分析的部件级模拟试验件设计验证方法 | |
CN105045989A (zh) | 一种基于虚拟维修的飞行器零部件维修性分析方法 | |
Adibfar et al. | Creation of a mock-up bridge digital twin by fusing intelligent transportation systems (ITS) Data into Bridge Information Model (BrIM) | |
Basri et al. | Conceptual design and simulation validation based finite element optimisation for tubercle leading edge composite wing of an unmanned aerial vehicle | |
KR102316603B1 (ko) | 건축물 내진성능 평가를 위한 초고속 데이터생성 엔진 및 그 방법 | |
Pervan et al. | Development of parametric CAD model and structural analysis of the car jack | |
CN108984909A (zh) | 一种基于Mar-Lin模型的含大损伤飞机复合材料结构剩余强度分析方法 | |
CN109726435B (zh) | 一种接头钉群载荷计算的方法 | |
CN104699873B (zh) | 一种飞机壁板开口结构数值分析方法 | |
CN111272764B (zh) | 大型智能临时平台非接触图像识别移动管控系统及方法 | |
Shukla et al. | Automated design optimization for the P2 and P8 hypersonic inlets | |
CN106844915B (zh) | 基于ansys节点耦合的盾构机箱涵吊机梁结构设计方法 | |
Coiro et al. | Design of low-speed aircraft by numerical and experimental techniques developed at DPA | |
CN105260502B (zh) | 确定各种因素对飞机刹车系统平均修复时间影响力的方法 | |
CN106709183A (zh) | 采用金属板替换的复合材料混杂接头的细节应力仿真方法 | |
Xiao et al. | Architecture-level particular risk modeling and analysis for a cyber-physical system with AADL | |
Margheri et al. | An uncertainty quantification analysis in a simplified problem of urban pollutant dispersion by means of ANOVA-POD/Kriging-based response surfaces | |
CN105868511B (zh) | 一种确定飞机加筋板单元及载荷的方法与装置 | |
Telenta et al. | Detached Eddy Simulation of the flow around a simplified vehicle sheltered by wind barrier in transient yaw crosswind |
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 |