CN103294856B - 一种复合材料层压板低速冲击损伤的数值模拟方法 - Google Patents
一种复合材料层压板低速冲击损伤的数值模拟方法 Download PDFInfo
- Publication number
- CN103294856B CN103294856B CN201310195634.2A CN201310195634A CN103294856B CN 103294856 B CN103294856 B CN 103294856B CN 201310195634 A CN201310195634 A CN 201310195634A CN 103294856 B CN103294856 B CN 103294856B
- Authority
- CN
- China
- Prior art keywords
- model
- laminate
- damage
- delamination
- interlayer
- 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
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明属于航空复合材料静损伤领域,特别是涉及到一种复合材料层压板低速冲击损伤的数值模拟方法,包括建立层压板单层模型的步骤、建立层压板的层间模型的步骤、建立层压板的树脂裂纹模型的步骤、建立冲头模型的步骤、组合成层压板冲击模型的步骤以及进行求解的步骤。本发明采用准静态载荷代替复杂的动态冲击载荷,减小了计算规模,提高了计算效率;采用Abaqus/Standard求解器,缩短计算时间;更好的模拟出在低速冲击载荷下较早发生的树脂裂纹,较好的模拟出分层损伤的面积和形状;很好的模拟了树脂裂纹与层间分层的相互作用,得到了准确的分层损伤的面积和形状;更好的模拟出分层损伤的面积和形状。
Description
技术领域
本发明属于航空复合材料静损伤领域,特别是涉及到一种复合材料层压板低速冲击损伤的数值模拟方法。
背景技术
复合材料层压板受到外来物的低速冲击损伤可能发生如分层、基体裂纹、纤维断裂和表面破坏等一系列目视勉强可见的损伤。存在此类损伤的层压板压缩强度可降低到无损板的50%左右。
国内外对复合材料低速冲击损伤的数值计算主要有三种方法:以强度破坏为准则的计算模型、断裂力学模型和界面单元模型。强度破坏模型不能很好的模拟分层和树脂裂纹的逐渐扩展过程,但是由于冲击损伤持续时间很短,这种模拟方法也被广泛应用。
断裂力学模型对网格构型和网格密度很敏感且需要在模型中布置初始裂纹,同时,这种方法还需要自适应网格技术来跟踪裂纹尖端的变化。
界面单元模型可以方便的模拟分层和树脂裂纹的起始和扩展,但是无法模拟界面中由冲击引起的压缩载荷对分层破坏的限制作用。
对复合材料层压板低速冲击损伤的数值模拟多采用动态载荷加载模型,耗时耗力,计算效率较低。
发明内容
本发明的目的是:提供一种能够减小计算规模、提高计算效率的一种复合材料层压板低速冲击损伤的数值模拟方法。
本发明的技术方案是:一种复合材料层压板低速冲击损伤的数值模拟方法,利用到Abaqus商业软件,包括以下步骤:
步骤一、采用连续壳单元SC8R,建立层压板单层模型,单元尺寸为H;
步骤二、采用界面单元COH3D8,建立层压板的层间模型,单元尺寸为H/5;
步骤三、采用界面单元COH3D8,建立层压板的树脂裂纹模型,单元尺寸为H/5;
步骤四、采用刚体单元,建立冲头模型;
步骤五、将单层模型、层间模型、树脂裂纹模型及冲头模型组合成层压板冲击模型;其中,层间模型引入层间摩擦力F,模拟层间压缩力对分层损伤的限制作用,摩擦系数为u;
步骤六、将层压板冲击模式提交到Abaqus/Standard求解器,进行求解;
步骤七、将求解结果与试验结果进行对比,若模拟结果与试验结果一致,则模拟完成,否则返回步骤五,修改摩擦系数u,继续进行求解,直至模拟结果与试验结果一致。
本发明的优点是:本发明采用采用准静态载荷代替复杂的动态冲击载荷,减小了计算规模,提高了计算效率;采用Abaqus/Standard求解器,较Abaqus/Explicit求解器较大的提高了计算效率,缩短计算时间;在单层面内纤维方向垂直布置界面单元模拟树脂裂纹,更好的模拟出在低速冲击载荷下较早发生的树脂裂纹;在层间布置界面单元模拟由树脂裂纹引起的层间分层,较好的模拟出分层损伤的面积和形状;很好的模拟了树脂裂纹与层间分层的相互作用,得到了准确的分层损伤的面积和形状;在受压载的层间引入层间摩擦载荷模拟压载对分层损伤的限制作用,更好的模拟出分层损伤的面积和形状。
附图说明
图1是层压板单层模型示意图;
图2是层间模型示意图;
图3是树脂裂纹模型示意图;
图4是冲头模型示意图;
图5是层压板冲击模型示意图。
图6是模拟结果示意图;
图7是试验结果示意图。
具体实施方式
下面结合附图对本发明做进一步详细描述,请参阅图1至图6。
一种复合材料层压板低速冲击损伤的数值模拟方法,利用到Abaqus商业软件,包括以下步骤:
步骤一、采用连续壳单元SC8R,建立层压板单层模型,单元尺寸为H。层压板包含多个单层模型,单层模型厚度为具有相同铺层角度的多个预浸料单向带或编织布的总厚度;如图1所示;
步骤二、采用界面单元COH3D8,建立层压板的层间模型,单元尺寸为H/5。层间模型厚度取值为层压板模型厚度的千分之一,并将其置于两个单层模型之间来模拟层间分层损伤;如图2所示;
步骤三、采用界面单元COH3D8,建立层压板的树脂裂纹模型,单元尺寸为H/5。树脂裂纹模型厚度取值为层压板模型厚度的千分之一,并将其垂直置于单层模型的中面位置,放置方向为单层模型的纤维方向;如图3所示;
步骤四、采用刚体单元,建立冲头模型;如图4所示;
步骤五、将单层模型、层间模型、树脂裂纹模型及冲头模型组合成层压板冲击模型,由对称关系,采用1/4模型进行计算分析;如图5所示;其中,相邻单层模型之间引入层间摩擦力F,模拟层间压缩载荷对分层损伤的限制作用,摩擦系数为u;
步骤六、将层压板冲击模式提交到Abaqus/Standard求解器,进行求解;
步骤七、将求解结果与试验结果进行对比,若模拟结果与试验结果一致,则模拟完成,否则返回步骤五,修改摩擦系数u,继续进行求解,直至模拟结果与试验结果一致。如图6及图7所示。
Claims (1)
1.一种复合材料层压板低速冲击损伤的数值模拟方法,利用到Abaqus商业软件,其特征在于,包括以下步骤:
步骤一、采用连续壳单元SC8R,建立层压板单层模型,单元尺寸为H;
步骤二、采用界面单元COH3D8,建立层压板的层间模型,单元尺寸为H/5;
步骤三、采用界面单元COH3D8,建立层压板的树脂裂纹模型,单元尺寸为H/5;
步骤四、采用刚体单元,建立冲头的模型;
步骤五、将单层模型、层间模型、树脂裂纹模型及冲头模型组合成层压板冲击模型;其中,层间模型引入层间摩擦力F,模拟层间压缩力对分层损伤的限制作用,摩擦系数为u;
步骤六、将层压板冲击模式提交到Abaqus/Standard求解器,进行求解;
步骤七、将求解结果与试验结果进行对比,若模拟结果与试验结果一致,则模拟完成,否则返回步骤五,修改摩擦系数u,继续进行求解,直至模拟结果与试验结果一致。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310195634.2A CN103294856B (zh) | 2013-01-05 | 2013-05-23 | 一种复合材料层压板低速冲击损伤的数值模拟方法 |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310002904.3 | 2013-01-05 | ||
CN201310002904 | 2013-01-05 | ||
CN2013100029043 | 2013-01-05 | ||
CN201310195634.2A CN103294856B (zh) | 2013-01-05 | 2013-05-23 | 一种复合材料层压板低速冲击损伤的数值模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103294856A CN103294856A (zh) | 2013-09-11 |
CN103294856B true CN103294856B (zh) | 2016-04-20 |
Family
ID=49095713
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310195634.2A Active CN103294856B (zh) | 2013-01-05 | 2013-05-23 | 一种复合材料层压板低速冲击损伤的数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103294856B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117610344A (zh) * | 2023-11-10 | 2024-02-27 | 中山大学 | 一种飞行汽车海上迫降结构损伤评价方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101002203A (zh) * | 2004-07-02 | 2007-07-18 | 恩格纽伊蒂有限公司 | 压碎模拟 |
CN102607935A (zh) * | 2011-12-27 | 2012-07-25 | 中国飞机强度研究所 | 一种含冲击损伤复合材料层压板剩余压缩强度的测算方法 |
-
2013
- 2013-05-23 CN CN201310195634.2A patent/CN103294856B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101002203A (zh) * | 2004-07-02 | 2007-07-18 | 恩格纽伊蒂有限公司 | 压碎模拟 |
CN102607935A (zh) * | 2011-12-27 | 2012-07-25 | 中国飞机强度研究所 | 一种含冲击损伤复合材料层压板剩余压缩强度的测算方法 |
Non-Patent Citations (2)
Title |
---|
复合材料低速冲击损伤分析方法;刘德博,关志华,陈建华,黎增山;《北京航空航天大学学报》;20120331;第38卷(第3期);参见424页,右栏第1.2.3节 * |
复合材料铺层板低速冲击作用下损伤的有限元分析;张彦,来新民,朱平,梁新华;《上海交通大学学报》;20060831(第8期);第1-4节,图3、5 * |
Also Published As
Publication number | Publication date |
---|---|
CN103294856A (zh) | 2013-09-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Abir et al. | On the relationship between failure mechanism and compression after impact (CAI) strength in composites | |
Li et al. | Repeated low-velocity impact response and damage mechanism of glass fiber aluminium laminates | |
CN103914593B (zh) | 层状复合材料力学行为的虚拟预测方法 | |
CN106777769B (zh) | 预测低速冲击下复合材料多层厚板渐进失效的有限元方法 | |
CN103559390B (zh) | 一种基于平均失效指数的复合材料π形胶接连接结构拉伸强度预测方法 | |
Fedele et al. | Three-dimensional effects induced by FRP-from-masonry delamination | |
CN103984802B (zh) | 一种蜂窝夹层结构的有限元建模简化方法 | |
Guida et al. | Certification by birdstrike analysis on C27J fullscale ribless composite leading edge | |
CN106066913A (zh) | 复杂复合材料结构等效材料性能多尺度计算方法 | |
CN107451307A (zh) | 一种多尺度计算复杂复合材料结构等效刚度矩阵的方法 | |
Xiaohui et al. | A new zig-zag theory and C 0 plate bending element for composite and sandwich plates | |
CN104133930A (zh) | 一种复合材料层合板损伤过程模拟方法 | |
CN103729498A (zh) | 一种大型复合材料加筋壁板结构设计方法 | |
CN106202597A (zh) | 一种复合材料加筋壁板结构优化分析方法 | |
Littell et al. | Effect of microscopic damage events on static and ballistic impact strength of triaxial braid composites | |
CN103294856B (zh) | 一种复合材料层压板低速冲击损伤的数值模拟方法 | |
CN106354965A (zh) | 一种含分层复合材料层压板剩余强度的有限元分析方法 | |
CN115391947A (zh) | 轨道交通车辆复合材料结构虚拟分析方法及系统 | |
Tong et al. | Material parameter modeling and solution technique using birth–death element for notched metallic panel repaired with bonded composite patch | |
Wang | Study on the application of carbon fiber composite materials in high-speed trains | |
Zhang et al. | Research on the limit operating pressure of a bedded salt gas storage cavern: A case study | |
Zhang et al. | Failure patterns of shear keys and seismic resistance of a gravity dam with longitudinal joints | |
Zang et al. | Simulation on high velocity impact process of windshield by SPH/FEM coupling method | |
CN106897483A (zh) | 一种基于任意两单元间插入内聚力建立冰雹模型的方法 | |
Yi | The progressive failure analysis of uni-directional fibre reinforced composite laminates |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |