CN106897483A - 一种基于任意两单元间插入内聚力建立冰雹模型的方法 - Google Patents
一种基于任意两单元间插入内聚力建立冰雹模型的方法 Download PDFInfo
- Publication number
- CN106897483A CN106897483A CN201710011726.9A CN201710011726A CN106897483A CN 106897483 A CN106897483 A CN 106897483A CN 201710011726 A CN201710011726 A CN 201710011726A CN 106897483 A CN106897483 A CN 106897483A
- Authority
- CN
- China
- Prior art keywords
- unit
- hail
- cohesive force
- model
- abaqus
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Stored Programmes (AREA)
- Electrophonic Musical Instruments (AREA)
Abstract
本发明涉及有限元建模领域,旨在提供一种基于任意两单元间插入内聚力建立冰雹模型的方法。包括下述过程:建立含所需冰雹部件的有限元模型;基于ABAQUS的PYTHON脚本语言实现任意两单元间插入内聚力的建模,形成带有内聚力单元的冰雹模型;使用带有内聚力单元的冰雹网格模型,结合复合材料层合板模型,实现冰雹冲击复合材料层合板的有限元仿真计算。本发明利用ABAQUS‑PYTHON脚本语言提供一种基于任意两单元间插入内聚力建立冰雹模型的方法,通过该方法能实现ABAQUS现有模块无法实现的任意两单元间插入内聚力的功能,建立带有内聚力单元的冰雹模型,且基于脚本语言的高效建模与ABAQUS完全兼容。
Description
技术领域
本发明是关于有限元建模领域,特别涉及一种基于任意两单元间插入内聚力建立冰雹模型的方法。
背景技术
当前,内聚力模型是弹塑性断裂力学中一个被广泛应用和研究的计算模型。内聚力模型避免了线弹性断裂力学中的裂纹尖端应力奇异性,计算了开裂过程中应力以及断裂能,使其在裂纹扩展,复合材料分层和脱粘,器件胶接等领域得到广泛重视和运用。
目前,基于内聚力模型的内聚力界面单元已经内嵌在ABAQUS商业有限元软件,基于ABAQUS有限元软件CAE建模功能可以在模型的内部面与面之间插入内聚力界面单元,这种插入内聚力界面单元的方式可以满足复合材料层合板分层和脱粘,器件胶接等有限元分析,但面-面间插入内聚力界面单元的方式因其局限性并不能满足所有的有限元分析,冰雹冲击复合材料层合板需要在冰雹模型的任意两单元间插入内聚力单元,现有的ABAQUS商业软件模块尚无法做到。因此,提供一种基于任意两单元间插入内聚力建立冰雹模型的方法尤为重要。
发明内容
本发明要解决的技术问题是,克服现有技术中的不足,提供一种基于任意两单元间插入内聚力建立冰雹模型的方法。
为解决上述技术问题,本发明的解决方案是:
提供一种基于任意两单元间插入内聚力建立冰雹模型的方法,包括下述过程:
一、建立含所需冰雹部件的有限元模型;
在ABAQUS-PART模块建立所需部件并命名,假定命名为“A”,在ABAQUS-PROPERTY模块分别为材料单元和内聚力界面单元建立两种材料属性,在ABABQUS-MESH模块为部件划分网格,此时“A”带有网格属性;
二、基于ABAQUS的PYTHON脚本语言实现任意两单元间插入内聚力的建模,形成带有内聚力单元的冰雹模型;
(1)读取ABAQUS冰雹部件的网格模型建立孤立网格模型,并单独命名为“A-1”,将“A-1”所有的节点信息和单元信息全部删除;
(2)复制“A”所有的节点并将节点组建成单元赋予给“A-1”,此部分的单元即为“A-1”的材料单元;
(3)将“A-1”公用面上的节点进行复制,形成节点对,组建零厚度内聚力单元;
(4)将材料单元和内聚力单元分别建立集合并赋予材料属性,建成带有内聚力单元的冰雹网格模型“A-1”;
三、使用带有内聚力单元的冰雹网格模型,结合复合材料层合板模型,实现冰雹冲击复合材料层合板的有限元仿真计算。
与现有技术相比,本发明的有益效果是:
本发明利用ABAQUS-PYTHON脚本语言提供一种基于任意两单元间插入内聚力建立冰雹模型的方法,通过该方法能实现ABAQUS现有模块无法实现的任意两单元间插入内聚力的功能,建立带有内聚力单元的冰雹模型,且基于脚本语言的高效建模与ABAQUS完全兼容。
附图说明
图1为本发明对所提出的基于任意两单元间插入内聚力建立冰雹模型的PYTHON语言实现流程图。
图2为本发明实施例中冰雹模型中内聚力单元集合的示意图。
图3为本发明实施例中生成的带有内聚力单元的冰雹网格模型。
具体实施方式
首先需要说明的是,本发明是计算机技术在有限元建模领域的一种应用。在本发明的实现过程中,会涉及到多个软件功能模块的应用。申请人认为,如在仔细阅读申请文件、准确理解本发明的实现原理和发明目的以后,在结合现有公知技术的情况下,本领域技术人员完全可以运用其掌握的软件编程技能实现本发明。凡本发明申请文件提及的均属此范畴,申请人不再一一列举。
一种基于任意两单元间插入内聚力建立冰雹模型的方法,包括下述过程:
一、建立含所需冰雹部件的有限元模型;
在ABAQUS-PART模块建立所需部件并命名,假定命名为“A”,在ABAQUS-PROPERTY模块分别为材料单元和内聚力界面单元建立两种材料属性,在ABABQUS-MESH模块为部件划分网格,此时“A”带有网格属性;
二、基于ABAQUS的PYTHON脚本语言实现任意两单元间插入内聚力的建模,形成带有内聚力单元的冰雹模型;
(1)读取ABAQUS冰雹部件的网格模型建立孤立网格模型,并单独命名为“A-1”,将“A-1”所有的节点信息和单元信息全部删除;
(2)复制“A”所有的节点并将节点组建成单元赋予给“A-1”,此部分的单元即为“A-1”的材料单元;
(3)将“A-1”公用面上的节点进行复制,形成节点对,组建零厚度内聚力单元;
(4)将材料单元和内聚力单元分别建立集合并赋予材料属性,建成带有内聚力单元的冰雹网格模型“A-1”。
三、使用带有内聚力单元的冰雹网格模型,结合复合材料层合板模型,实现冰雹冲击复合材料层合板的有限元仿真计算。
以下用基于上述过程应用的具体实施例子来说明本发明的实现效果:
首先在ABAQUS-PART模块建立球形冰雹部件,并在MESH模块划分网格,在PROPERTY模块设置材料单元和内聚力单元的材料属性。根据流程图如图1使用ABAQUS-PYTHON脚本语言,复制部件网格模型为孤立网格模型,删除孤立网格模型上的节点和单元信息,再复制部件网格模型的节点建立孤立网格模型材料单元,再将孤立网格模型共同面的节点复制形成节点对,组建零厚度内聚力单元;最后将材料单元和内聚力单元分别建立集合赋予材料属性。形成的任意两单元间的内聚力单元集合如图2,形成的带有内聚力单元的冰雹网格模型如图3。
本发明利用ABAQUS-PYTHON脚本语言一种基于任意两单元间插入内聚力建立冰雹模型的方法,通过该方法实现ABAQUS现有模块无法实现的任意两单元间插入内聚力的功能,建立带有内聚力的冰雹模型,且基于脚本语言的高效建模与ABAQUS完全兼容。
最后,需要注意的是,以上列举的仅是本发明的具体实施例。显然,本发明不限于以上实施例,还可以有很多变形。本领域的普通技术人员能从本发明公开的内容中直接导出或联想到的所有变形,均应认为是本发明的保护范围。
Claims (1)
1.一种基于任意两单元间插入内聚力建立冰雹模型的方法,其特征在于,包括下述过程:
一、建立含所需冰雹部件的有限元模型;
在ABAQUS-PART模块建立所需部件并命名,假定命名为“A”,在ABAQUS-PROPERTY模块分别为材料单元和内聚力界面单元建立两种材料属性,在ABABQUS-MESH模块为部件划分网格,此时“A”带有网格属性;
二、基于ABAQUS的PYTHON脚本语言实现任意两单元间插入内聚力的建模,形成带有内聚力单元的冰雹模型;
(1)读取ABAQUS冰雹部件的网格模型建立孤立网格模型,并单独命名为“A-1”,将“A-1”所有的节点信息和单元信息全部删除;
(2)复制“A”所有的节点并将节点组建成单元赋予给“A-1”,此部分的单元即为“A-1”的材料单元;
(3)将“A-1”公用面上的节点进行复制,形成节点对,组建零厚度内聚力单元;
(4)将材料单元和内聚力单元分别建立集合并赋予材料属性,建成带有内聚力单元的冰雹网格模型“A-1”;
三、使用带有内聚力单元的冰雹网格模型,结合复合材料层合板模型,实现冰雹冲击复合材料层合板的有限元仿真计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710011726.9A CN106897483B (zh) | 2017-01-08 | 2017-01-08 | 一种基于任意两单元间插入内聚力建立冰雹模型的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710011726.9A CN106897483B (zh) | 2017-01-08 | 2017-01-08 | 一种基于任意两单元间插入内聚力建立冰雹模型的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106897483A true CN106897483A (zh) | 2017-06-27 |
CN106897483B CN106897483B (zh) | 2019-05-17 |
Family
ID=59198461
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710011726.9A Active CN106897483B (zh) | 2017-01-08 | 2017-01-08 | 一种基于任意两单元间插入内聚力建立冰雹模型的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106897483B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113297670A (zh) * | 2021-05-19 | 2021-08-24 | 西北工业大学 | 基于近场动力学对冰雹冲击飞机复合材料层合板的建模方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103605865A (zh) * | 2013-11-28 | 2014-02-26 | 北京航空航天大学 | 一种基于内聚力模型的复合材料多向铺层板分层扩展行为模拟方法 |
CN105158256A (zh) * | 2015-09-30 | 2015-12-16 | 山东大学 | 一种复合材料在线健康监测系统和监测方法 |
CN105760623A (zh) * | 2016-03-16 | 2016-07-13 | 中国直升机设计研究所 | 一种直升机复合材料主桨叶许用缺陷确定方法 |
EP3098263A1 (en) * | 2015-05-29 | 2016-11-30 | ATLANTA - Componentes Para Calçado, Lda. | Composite material, obtaining method and use thereof |
-
2017
- 2017-01-08 CN CN201710011726.9A patent/CN106897483B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103605865A (zh) * | 2013-11-28 | 2014-02-26 | 北京航空航天大学 | 一种基于内聚力模型的复合材料多向铺层板分层扩展行为模拟方法 |
EP3098263A1 (en) * | 2015-05-29 | 2016-11-30 | ATLANTA - Componentes Para Calçado, Lda. | Composite material, obtaining method and use thereof |
CN105158256A (zh) * | 2015-09-30 | 2015-12-16 | 山东大学 | 一种复合材料在线健康监测系统和监测方法 |
CN105760623A (zh) * | 2016-03-16 | 2016-07-13 | 中国直升机设计研究所 | 一种直升机复合材料主桨叶许用缺陷确定方法 |
Non-Patent Citations (3)
Title |
---|
SHO FUNAI: "Hail Ice Impact on Composite Structures at Glancing Angles", 《UNIVERSITY OF CALIFORNIA》 * |
林荣会: "夹层玻璃受冰雹冲击的有限元模拟和冲击点的动态响应", 《硅酸盐通报》 * |
黄兴: "航空飞行器冰雹冲击研究进展", 《佳木斯大学学报(自然科学版)》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113297670A (zh) * | 2021-05-19 | 2021-08-24 | 西北工业大学 | 基于近场动力学对冰雹冲击飞机复合材料层合板的建模方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106897483B (zh) | 2019-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hofmann et al. | A grain based modeling study of mineralogical factors affecting strength, elastic behavior and micro fracture development during compression tests in granites | |
CN108804787A (zh) | 基于批量插入内聚力单元模拟岩桥贯通的方法 | |
Viola et al. | Soft core plane state structures under static loads using GDQFEM and cell method | |
CN108804735A (zh) | 低速冲击下复合材料层合板渐进失效的多尺度预测方法 | |
CN103914593A (zh) | 层状复合材料力学行为的虚拟预测方法 | |
McElroy et al. | Simulation of delamination–migration and core crushing in a CFRP sandwich structure | |
Zheng et al. | Arch-dam crack deformation monitoring hybrid model based on XFEM | |
CN106874547A (zh) | 一种预测纤维全缠绕复合材料气瓶爆破压力的方法 | |
CN106897483A (zh) | 一种基于任意两单元间插入内聚力建立冰雹模型的方法 | |
CN107193559A (zh) | 一种基于元模型的复杂系统虚拟样机设计方法 | |
Paśnik et al. | Application of the CZM Technique to Delamination Analysis of Coupled Laminate Beams | |
Tong et al. | Material parameter modeling and solution technique using birth–death element for notched metallic panel repaired with bonded composite patch | |
CN110472307B (zh) | 一种显示系统光学胶剥离仿真方法 | |
Akpinar et al. | Effect of protrusion at the ends of bondline in single lap joints under tension and bending | |
Dixit et al. | Finite element analysis of quasi‐static indentation of woven fabric textile composites using different nose shape indenters: Finite‐Elemente‐Analyse von quasi statischen Eindrücken mit verschiedenen Indentergeometrien bei textilen Kompositen | |
CN106156441A (zh) | 宏纤维压电复合材料压电片动力学的作动力识别方法 | |
CN103294856B (zh) | 一种复合材料层压板低速冲击损伤的数值模拟方法 | |
Mi et al. | A research on loss factor of embedded and co-cured composite damping structure | |
Thomson et al. | Numerical prediction of the ballistic performance of hygrothermally aged CFRP laminates using a multi-scale modelling approach | |
Zheng et al. | Effect of Notch Depth on the Butt-Joint under Izod Impact Test with Material Properties and Mechanics Analysis | |
Nocera et al. | Correlation studies for the in-plane analysis of masonry walls based on macroscopic FE models with damage | |
Benedetti et al. | Nonlinear free vibrations analysis of cracked composite stiffened plates via X-Ritz approach | |
Ghelli et al. | Comparison between low velocity impact and quasi-static indentation tests on CFRP composite laminates | |
Chen et al. | Experimental Study and Finite Analysis on Flexural Capacity of Laminated Glass | |
Shor et al. | Through-Thickness Element Splitting for Simulation of Delamination in Laminated Composite Materials |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |