CN105740566B - 一种预测层状复合材料层内损伤和层间分层的有限元方法 - Google Patents
一种预测层状复合材料层内损伤和层间分层的有限元方法 Download PDFInfo
- Publication number
- CN105740566B CN105740566B CN201610090758.8A CN201610090758A CN105740566B CN 105740566 B CN105740566 B CN 105740566B CN 201610090758 A CN201610090758 A CN 201610090758A CN 105740566 B CN105740566 B CN 105740566B
- Authority
- CN
- China
- Prior art keywords
- damage
- abaqus
- model
- composite
- interface
- 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
- 239000002131 composite material Substances 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000032798 delamination Effects 0.000 title claims abstract description 14
- 238000006073 displacement reaction Methods 0.000 claims abstract description 18
- 239000007787 solid Substances 0.000 claims abstract description 12
- 239000010410 layer Substances 0.000 claims description 17
- 239000000463 material Substances 0.000 claims description 12
- 239000000835 fiber Substances 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 7
- 239000002356 single layer Substances 0.000 claims description 3
- 238000005336 cracking Methods 0.000 claims description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 6
- 229910052799 carbon Inorganic materials 0.000 description 6
- 239000011347 resin Substances 0.000 description 6
- 229920005989 resin Polymers 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 150000001875 compounds Chemical class 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000011230 binding agent Substances 0.000 description 1
- 230000005494 condensation Effects 0.000 description 1
- 238000009833 condensation Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000002023 wood Substances 0.000 description 1
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)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及复合材料损伤领域,旨在提供一种预测层状复合材料层内损伤和层间分层的有限元方法。该种预测层状复合材料层内损伤和层间分层的有限元方法,包括过程:建立复合材料损伤模型、基于ABAQUS软件实现提出的损伤模型。本发明基于有限元方法,使用ABAQUS用户子程序模块,对层内实体单元和层间界面单元进行实时组装,统一求解节点位移和节点力,同时预测层状复合材料层内损伤和层间分层,能全面真实地预测层状复合材料的损伤演化过程。
Description
技术领域
本发明是关于复合材料损伤领域,特别涉及一种预测层状复合材料层内损伤和层间分层的有限元方法。
背景技术
当前,复合材料正广泛应用于航空航天、风力发电、压力容器、汽车等高新技术领域。层状复合材料结构是实现高刚度、高强度轻量化设计的重要途径之一。然而,目前复合材料轻量化设计遇到的挑战问题在于其复杂的损伤失效机理,特别是层内损伤和层间分层之间的相互作用。
国内外对层状复合材料的失效机理和强度进行了大量的理论和数值分析。对于复合材料层内损伤,涉及到失效准则、损伤本构模型、损伤演化模型和有限元算法。对于复合材料层间分层,主要采用内聚力模型预测。近年来,有限元法已应用于复合材料结构损伤研究中。然而,国内外同时预测层状复合材料层内损伤和层间分层的有限元方法很少。
ABAQUS软件可以对复合材料的纤维和基体以及接触界面进行渐进损伤和失效分析。ABAQUS中对复合材料分层破坏的模拟有两种方式,即虚拟裂纹闭合技术和内聚力模型。然而虚拟裂纹闭合技术不能模拟裂纹萌生,ABAQUS自带的内聚力模型为双线性,没有其他的内聚力形状,不能模拟零厚度内聚力界面单元。
发明内容
本发明的主要目的在于克服现有技术中的不足,提供一种结合ABAQUS-UMAT和UEL子程序模块发展数值技术,预测复合材料损伤演化特性的有限元方法。为解决上述技术问题,本发明的解决方案是:
提供一种预测层状复合材料层内损伤和层间分层的有限元方法,包括下述过程:
一、建立复合材料损伤模型;
二、基于ABAQUS软件实现提出的损伤模型;
所述过程一具体包括下述步骤:
步骤(1):建立层内损伤模型:
各向异性层内损伤应力-应变本构方程为:σ=Cd:ε:
其中,σ为应力,ε为应变,Cd(1-di(i=1,2,3),C)是指各向异性损伤四阶弹性张量,所述C为未损伤的四阶弹性张量,di(i=1,2,3)分别代表纤维断裂、基体开裂和纤维/基体界面分离损伤变量;
步骤(2):建立层间分层模型:
对于单一模式分层,内聚力损伤模型描述为界面牵引力Ti与位移跳[[ui]]之间的本构关系:
其中,Ti为界面牵引力,e为自然对数的底数,Ti c(i=1,2,3)为最大牵引力,[[ui]]为位移跳,为损伤变量,[[ui]]c(i=1,2,3)为临界位移跳,[[ui]]f(i=1,2,3)为失效时的位移跳;
对于混合模式分层,内聚力损伤模型描述为界面牵引力T与位移跳[[u]]之间的本构关系:
其中,T为界面牵引力,e为自然对数的底数,[[u]]为位移跳,[[u]]c为临界位移跳,tc为最大牵引力,ds是对应于混合模式分层的损伤变量;
步骤(3):求解整体刚度方程:
其中,A代表实体单元和界面单元刚度组装,nbulk和ncoh分别为实体单元和内聚力单元的数量,KL为实体单元刚度,Kc为界面单元刚度,Δu为节点位移增量,Fext为外部节点力;
所述过程二具体包括下述步骤:
步骤(4):使用FORTRAN语言编写ABAQUS用户子程序数值代码,定义材料参数和内聚界面单元(当ABAQUS没有提供我们需要的材料和单元类型时,可以运用ABAQUS用户子程序来自己定义);
使用ABAQUS-UMAT用户材料子程序,定义过程一建立的复合材料损伤模型中的变量(应力、应变、损伤变量等),使用ABAQUS-UEL用户单元子程序定义零厚度内聚界面单元;
步骤(5):在ABAQUS中建立复合材料层合板的有限元模型,单层复合材料层合板采用实体单元,层与层之间的界面采用内聚界面单元;
步骤(6):根据过程一中建立的复合材料损伤模型,基于ABAQUS-UMAT用户材料子程序模块预测层内损伤演化过程,基于ABAQUS-UEL用户单元子程序模块预测分层过程。
与现有技术相比,本发明的有益效果是:
本发明基于有限元方法,使用ABAQUS用户子程序模块,对层内实体单元和层间界面单元进行实时组装,统一求解节点位移和节点力,同时预测层状复合材料层内损伤和层间分层,能全面真实地预测层状复合材料的损伤演化过程。
附图说明
图1为本发明预测层状复合材料损伤特性的流程图。
图2为实施例含中心孔的碳纤维树脂基层状复合材料的有限元模型示意图。
图3为实施例含中心孔的碳纤维树脂基层状复合材料基体损伤过程示意图。
图4为实施例含中心孔的碳纤维树脂基层状复合材料分层示意图。
图5为实施例含中心孔的碳纤维树脂基层状复合材料数值模拟结果与实验结果的对比示意图。
具体实施方式
下面结合附图与具体实施方式对本发明作进一步详细描述:
首先需要说明的是,本发明是计算机技术在复合材料损伤预测领域的一种应用。在本发明的实现过程中,会涉及到多个软件功能模块的应用。申请人认为,如在仔细阅读申请文件、准确理解本发明的实现原理和发明目的以后,在结合现有公知技术的情况下,本领域技术人员完全可以运用其掌握的软件编程技能实现本发明。前述软件功能模块包括但不限于:ABAQUS用户子程序数值代码等,凡本发明申请文件提及的均属此范畴,申请人不再一一列举。
如图1所示,本发明预测层状复合材料损伤演化特性的流程如下:
(1)在ABAQUS中建立复合材料层合板的有限元模型,单层复合材料层合板采用实体单元,层与层之间的界面采用内聚界面单元;
(2)给定第一次迭代时的初始节点位移;
(3)使用ABAQUS-UMAT用户材料子程序,计算层内实体单元应力、单元应变、损伤演化特性和单元刚度,使用ABAQUS-UEL用户单元子程序计算内聚力界面单元牵引力、损伤演化特性和单元刚度;
(4)在ABAQUS主文件里,为实体单元和内聚力界面单元组装单元刚度和节点力,计算节点位移;
(5)判断是否产生新的失效纤维和基体单元?若是,对实体单元更新损伤变量和折减刚度,提高外部荷载增量作为下次迭代,转到(3)继续计算,若否,则进入下一步;
(6)判断层状复合材料是否垮塌?若否,在下个载荷步中提高外部荷载增量,转到(3)继续计算,若垮塌,则计算结束。
本发明的一个具体实施例是模拟含中心孔的碳纤维树脂基层状复合材料的拉伸失效过程,材料为T700/8911,孔直径为10mm。由于结构的对称性,取四分之一建立有限元模型,如图2所示。基体的损伤演化过程如图3所示,其中(a)、(b)、(c)分别对应应变为0.672%、1.001%、1.622%三个阶段。图4为分层失效过程,其中(a)、(b)、(c)分别对应应力为σmax=τmax=5MPa、σmax=τmax=10MPa、σmax=τmax=20MPa。从图5可以看出,本发明的模拟结果与实验结果吻合较好。
本发明结合层内损伤和层间分层模型,通过有限元数值模拟,能够准确地预测带有中心孔的碳纤维树脂基层状复合材料层状复合材料的损伤演化过程,为深入阐明复合材料结构的损伤失效特性、提升轻量化强度设计水平提供了技术支撑。
最后,需要注意的是,以上列举的仅是本发明的具体实施例。显然,本发明不限于以上实施例,还可以有很多变形。本领域的普通技术人员能从本发明公开的内容中直接导出或联想到的所有变形,均应认为是本发明的保护范围。
Claims (1)
1.一种预测层状复合材料层内损伤和层间分层的有限元方法,其特征在于,包括下述过程:
一、建立复合材料损伤模型;
二、基于ABAQUS软件实现提出的损伤模型;
所述过程一具体包括下述步骤:
步骤(1):建立层内损伤模型:
各向异性层内损伤应力-应变本构方程为:σ=Cd:ε;
其中,σ为应力,ε为应变,Cd(1-di(i=1,2,3),C)是指各向异性损伤四阶弹性张量,所述C为未损伤的四阶弹性张量,di(i=1,2,3)分别代表纤维断裂、基体开裂和纤维/基体界面分离损伤变量;
步骤(2):建立层间分层模型:
对于单一模式分层,内聚力损伤模型描述为界面牵引力Ti与位移跳[[ui]]之间的本构关系:
其中,Ti为界面牵引力,e为自然对数的底数,Ti c(i=1,2,3)为最大牵引力,[[ui]]为位移跳,为损伤变量,[[ui]]c(i=1,2,3)为临界位移跳;
对于混合模式分层,内聚力损伤模型描述为界面牵引力T与位移跳[[u]]之间的本构关系:
其中,T为界面牵引力,e为自然对数的底数,[[u]]为位移跳,[[u]]c为临界位移跳,tc为最大牵引力,ds是对应于混合模式分层的损伤变量;
步骤(3):求解整体刚度方程:
其中,A代表实体单元和界面单元刚度组装,nbulk和ncoh分别为实体单元和内聚力单元的数量,KL为实体单元刚度,Kc为界面单元刚度,Δu为节点位移增量,Fext为外部节点力;
所述过程二具体包括下述步骤:
步骤(4):使用FORTRAN语言编写ABAQUS用户子程序数值代码,定义材料参数和内聚界面单元;
使用ABAQUS-UMAT用户材料子程序,定义过程一建立的复合材料损伤模型中的变量,使用ABAQUS-UEL用户单元子程序定义零厚度内聚界面单元;
步骤(5):在ABAQUS中建立复合材料层合板的有限元模型,单层复合材料层合板采用实体单元,层与层之间的界面采用内聚界面单元;
步骤(6):根据过程一中建立的复合材料损伤模型,基于ABAQUS-UMAT用户材料子程序模块预测层内损伤演化过程,基于ABAQUS-UEL用户单元子程序模块预测分层过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610090758.8A CN105740566B (zh) | 2016-02-18 | 2016-02-18 | 一种预测层状复合材料层内损伤和层间分层的有限元方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610090758.8A CN105740566B (zh) | 2016-02-18 | 2016-02-18 | 一种预测层状复合材料层内损伤和层间分层的有限元方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105740566A CN105740566A (zh) | 2016-07-06 |
CN105740566B true CN105740566B (zh) | 2018-06-22 |
Family
ID=56245285
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610090758.8A Expired - Fee Related CN105740566B (zh) | 2016-02-18 | 2016-02-18 | 一种预测层状复合材料层内损伤和层间分层的有限元方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105740566B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106503292B (zh) * | 2016-09-20 | 2018-04-24 | 浙江大学 | 预测低速冲击下复合材料层合板渐进失效的有限元方法 |
CN108021721A (zh) * | 2016-10-31 | 2018-05-11 | 大连比亚迪汽车有限公司 | 客车薄壁复合材料间再夹芯结构的有限元建模方法、装置 |
CN106777769B (zh) * | 2017-01-08 | 2018-04-24 | 浙江大学 | 预测低速冲击下复合材料多层厚板渐进失效的有限元方法 |
CN107463720A (zh) * | 2017-03-17 | 2017-12-12 | 浙江工业大学 | 一种评估多层薄膜膜基结合强度的方法 |
CN107832492B (zh) * | 2017-10-11 | 2021-07-06 | 大连理工大学 | 一种基于内聚力模型的钢结构腐蚀疲劳损伤计算方法 |
CN108595754A (zh) * | 2018-03-20 | 2018-09-28 | 南京航空航天大学 | 层间增韧复合材料层合板的仿真方法 |
CN108256281B (zh) * | 2018-03-26 | 2020-03-10 | 中国矿业大学 | 一种考虑搭接界面形貌及搭接物梯度性质的强度预测方法 |
CN109115629A (zh) * | 2018-07-18 | 2019-01-01 | 西安近代化学研究所 | 一种层间增韧复合材料损伤与演化失效的计算方法 |
CN111368389B (zh) * | 2019-10-11 | 2023-11-07 | 暨南大学 | 一种预测复合材料层合板失效强度的方法 |
CN110941922A (zh) * | 2019-11-21 | 2020-03-31 | 北京航空航天大学 | 一种预测复合材料层压板层间和纤维带间损伤的仿真方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005078556A (ja) * | 2003-09-03 | 2005-03-24 | Yokohama Rubber Co Ltd:The | 構造体の有限要素モデル作成方法および構造体のシミュレーション方法 |
CN103970969A (zh) * | 2014-05-30 | 2014-08-06 | 哈尔滨工业大学 | 一种使用有限元方法模拟编织复合材料双轴试验确定材料性能参数的方法 |
CN105205223A (zh) * | 2015-08-27 | 2015-12-30 | 湘潭大学 | 一种建立含缺陷的材料模型的有限元建模方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009190427A (ja) * | 2008-02-12 | 2009-08-27 | Yokohama Rubber Co Ltd:The | タイヤのシミュレーション方法 |
-
2016
- 2016-02-18 CN CN201610090758.8A patent/CN105740566B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005078556A (ja) * | 2003-09-03 | 2005-03-24 | Yokohama Rubber Co Ltd:The | 構造体の有限要素モデル作成方法および構造体のシミュレーション方法 |
CN103970969A (zh) * | 2014-05-30 | 2014-08-06 | 哈尔滨工业大学 | 一种使用有限元方法模拟编织复合材料双轴试验确定材料性能参数的方法 |
CN105205223A (zh) * | 2015-08-27 | 2015-12-30 | 湘潭大学 | 一种建立含缺陷的材料模型的有限元建模方法 |
Non-Patent Citations (1)
Title |
---|
模糊随机损伤力学及模糊随机损伤有限元在岩土工程中的应用;王亚军;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20091015(第10期);第C038-57页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105740566A (zh) | 2016-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105740566B (zh) | 一种预测层状复合材料层内损伤和层间分层的有限元方法 | |
He | A review of finite element analysis of adhesively bonded joints | |
Camanho et al. | Hybrid titanium–CFRP laminates for high-performance bolted joints | |
Liu et al. | Recent developments on damage modeling and finite element analysis for composite laminates: A review | |
Huang et al. | Characterization and analysis of carbon fibre-reinforced polymer composite laminates with embedded circular vasculature | |
Ribeiro et al. | Damage model and progressive failure analyses for filament wound composite laminates | |
Wang et al. | Stepped flush repairs for primary composite structures | |
Purimpat et al. | Effect of fiber angle orientation on a laminated composite single-lap adhesive joint | |
Fu et al. | Strength prediction of composite π joint under bending load and study of geometric and material variations effects | |
McElroy | Use of an enriched shell finite element to simulate delamination-migration in a composite laminate | |
Yousefsani et al. | Analytical solution of stress field in adhesively bonded composite single-lap joints under mechanical loadings | |
Vyas et al. | Computational implementation of a novel constitutive model for multidirectional composites | |
Yao et al. | A novel total fatigue life model for delamination growth in composite laminates under generic loading | |
Sgambitterra et al. | A robust three-node shell element for laminated composites with matrix damage | |
Zhang et al. | A novel energy‐based framework for characterizing the strain‐softening behavior of CFRP composites using cyclic loading | |
Zheng et al. | Improving damping performance of fiberglass-reinforced resin matrix composites by designing viscoelastic sandwich layer | |
Tong et al. | Material parameter modeling and solution technique using birth–death element for notched metallic panel repaired with bonded composite patch | |
Masmanidis et al. | Modeling the progressive failure of scarf joints for wind turbine blade repair | |
Wang et al. | Development of a finite element model for progressive damage analysis of composite laminates subjected to low velocity impact | |
Versino | Refined theories and discontinuous Galerkin methods for the analysis of multilayered composite structures | |
Allegri et al. | A simplified approach to the damage tolerance design of asymmetric tapered laminates. Part II: Methodology validation | |
Harvey et al. | Numerical and analytical study of delamination in composite laminates | |
Aldoumani et al. | The robustness of carbon fibre members bonded to aluminium connectors in aerial delivery systems | |
Ogasawara et al. | Linear and nonlinear torsional behavior of unidirectional CFRP and GFRP | |
Han et al. | Multiscale analysis of the correlation of processing parameters on viscidity of composites fabricated by automated fiber placement |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into 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: 20180622 |
|
CF01 | Termination of patent right due to non-payment of annual fee |