CN106202691A - 一种预报复合材料基体开裂的失效判据 - Google Patents
一种预报复合材料基体开裂的失效判据 Download PDFInfo
- Publication number
- CN106202691A CN106202691A CN201610527455.8A CN201610527455A CN106202691A CN 106202691 A CN106202691 A CN 106202691A CN 201610527455 A CN201610527455 A CN 201610527455A CN 106202691 A CN106202691 A CN 106202691A
- Authority
- CN
- China
- Prior art keywords
- matrix
- cracking
- represent
- shear
- strain energy
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)
- Casting Or Compression Moulding Of Plastics Or The Like (AREA)
Abstract
一种预报复合材料基体开裂的失效判据,该方法有四大步骤:步骤一、确定复杂应力状态下复合材料基体主要失效模式;步骤二、计算基体开裂前储存的弹性应变能;步骤三、计算基体开裂时释放的弹性应变能;步骤四、建立基体开裂的失效判据;本发明特点是提出了一种预报复合材料基体开裂的失效判据,可避免裂纹尖端应力场奇异问题,更好地描述复合材料基体开裂部位受力的强弱程度,具有使用简单、方便,效率高的优点。
Description
技术领域
本发明给出了一种预报复合材料基体开裂的失效判据,属于复合材料设计领域。
背景技术
在纤维增强复合材料层合板的使用过程中,基体开裂是起始最早、最常发生的一类失效模式。尤其是对于具有90°叠层的层板复合材料,90°层中基体是主要的载荷承担者,在外部载荷远小于其极限载荷时,由于基体断裂应变小于纤维的断裂应变,基体横向层将首先垂直于拉伸载荷方向发生开裂。由于复合材料基体开裂尖端存在应力场奇异问题,传统的局部应力/应变方法难以表征复合材料基体开裂部位受力的大小,而能量判据则可避免裂纹尖端的应力场奇异问题,更好地描述复合材料基体开裂部位受力的强弱程度。因此,本发明提出了一种预报复合材料基体开裂的失效判据,为复合材料结构的损伤容限设计提供技术支持。
发明内容
本发明提出了一种预报复合材料基体开裂的判据,为复合材料的损伤容限设计提供技术支持。本发明的技术方案如下:
步骤一、确定复杂应力状态下复合材料基体主要失效模式
复杂应力状态下,复合材料基体的主要失效模式为相邻纤维间的基体开裂。在横向拉/压力以及剪切力作用下,单层板相邻纤维之间的树脂基体可能发生的失效模式包括横向拉伸/压缩开裂和面内/面外剪切开裂(如图1所示)。
步骤二、计算基体开裂前储存的弹性应变能
在复杂应力状态(即σ22、τ12和τ23)下,开裂前基体储存的弹性应变能密度分别为
其中,下标1、2分别代表复合材料的纵向和横向,3代表复合材料的厚度方向。u22代表基体开裂前的横向拉压弹性应变能密度,u12代表基体开裂前的面内剪切弹性应变能密度,u23代表基体开裂前的面外剪切弹性应变能密度。E22代表基体的横向拉伸和压缩模量,G12代表基体的面内剪切模量,G23代表基体的面外剪切模量。σ22代表基体的横向正应力,τ12和τ23分别代表基体的面内和面外切应力。
步骤三、计算基体开裂时释放的弹性应变能
在复杂应力状态(即σ22、τ12和τ23)下,基体发生开裂失效时释放的弹性应变能密度分别为
其中,u22,cr代表基体开裂时释放的临界横向拉压弹性应变能密度,u12,cr代表基体开裂时释放的临界面内剪切弹性应变能密度,u23,cr代表基体开裂时释放的临界面外剪切弹性应变能密度。Y代表基体的横向拉伸和压缩强度,S12和S23代表基体的面内和面外剪切强度。
步骤四、建立基体开裂的失效判据
定义f22、f12和f23分别为基体横向拉/压失效、面内剪切失效以及面外剪切失效的能量失效指数,其表达式可写为
将横向拉压失效、面内剪切失效以及面外剪切失效的能量失效指数相加,得到基体开裂总失效指数,即
基体开裂的失效判据可写为F≥1。
优点及功效:本发明提供的失效判据则可避免裂纹尖端的应力场奇异问题,更好地描述复合材料基体开裂部位受力的强弱程度,具有使用简单、方便,效率高的优点,能为复合材料结构的损伤容限设计提供技术支持。
附图说明
图1为基体的可能破坏模式示意图。
图2为基体开裂失效判据发明方法流程图。
图中符号说明如下:
图1中的σ22为基体的横向正应力分量,τ12和τ23分别为基体的面内和面外切应力分量。
具体实施方式
图2为本发明所述方法的流程框图,本发明分四步实现,具体为:
步骤一、确定复杂应力状态下复合材料基体主要失效模式
复杂应力状态下,复合材料基体的主要失效模式为相邻纤维间的基体开裂。在横向拉/压力以及剪切力作用下,单层板相邻纤维之间的树脂基体可能发生的失效模式包括横向拉伸/压缩开裂和面内/面外剪切开裂(如图1所示)。
步骤二、计算基体开裂前储存的弹性应变能
在复杂应力状态(即σ22、τ12和τ23)下,开裂前基体储存的弹性应变能密度分别为
其中,下标1、2分别代表复合材料的纵向和横向,3代表复合材料的厚度方向。u22代表基体开裂前的横向拉压弹性应变能密度,u12代表基体开裂前的面内剪切弹性应变能密度,u23代表基体开裂前的面外剪切弹性应变能密度。E22代表基体的横向拉伸和压缩模量,G12代表基体的面内剪切模量,G23代表基体的面外剪切模量。σ22代表基体的横向正应力,τ12和τ23分别代表基体的面内和面外切应力。
步骤三、计算基体开裂时释放的弹性应变能
在复杂应力状态(即σ22、τ12和τ23)下,基体发生开裂失效时释放的弹性应变能密度分别为
其中,u22,cr代表基体开裂时释放的临界横向拉压弹性应变能密度,u12,cr代表基体开裂时释放的临界面内剪切弹性应变能密度,u23,cr代表基体开裂时释放的临界面外剪切弹性应变能密度。Y代表基体的横向拉伸和压缩强度,S12和S23代表基体的面内和面外剪切强度。
步骤四、建立基体开裂的失效判据
定义f22、f12和f23分别为基体横向拉/压失效、面内剪切失效以及面外剪切失效的能量失效指数,其表达式可写为
将横向拉压失效、面内剪切失效以及面外剪切失效的能量失效指数相加,得到基体开裂总失效指数,即
基体开裂的失效判据可写为F≥1。
Claims (1)
1.一种预报复合材料基体开裂的失效判据,该判据的具体步骤如下:
步骤一、确定复杂应力状态下复合材料基体主要失效模式
复杂应力状态下,复合材料基体的主要失效模式为相邻纤维间的基体开裂。在横向拉/压力以及剪切力作用下,单层板相邻纤维之间的树脂基体可能发生的失效模式包括横向拉伸/压缩开裂和面内/面外剪切开裂。
步骤二、计算基体开裂前储存的弹性应变能
在复杂应力状态(即σ22、τ12和τ23)下,开裂前基体储存的弹性应变能密度分别为
其中,下标1、2分别代表复合材料的纵向和横向,3代表复合材料的厚度方向。u22代表基体开裂前的横向拉压弹性应变能密度,u12代表基体开裂前的面内剪切弹性应变能密度,u23代表基体开裂前的面外剪切弹性应变能密度。E22代表基体的横向拉伸和压缩模量,G12代表基体的面内剪切模量,G23代表基体的面外剪切模量。σ22代表基体的横向正应力,τ12和τ23分别代表基体的面内和面外切应力。
步骤三、计算基体开裂时释放的弹性应变能
在复杂应力状态(即σ22、τ12和τ23)下,基体发生开裂失效时释放的弹性应变能密度分别为
其中,u22,cr代表基体开裂时释放的临界横向拉压弹性应变能密度,u12,cr代表基体开裂时释放的临界面内剪切弹性应变能密度,u23,cr代表基体开裂时释放的临界面外剪切弹性应变能密度。Y代表基体的横向拉伸和压缩强度,S12和S23代表基体的面内和面外剪切强度。
步骤四、建立基体开裂的失效判据
定义f22、f12和f23分别为基体横向拉/压失效、面内剪切失效以及面外剪切失效的能量失效指数,其表达式可写为
将横向拉压失效、面内剪切失效以及面外剪切失效的能量失效指数相加,得到基体开裂总失效指数,即
基体开裂的失效判据可写为F≥1。
功能特征:避免复合材料基体裂纹尖端的应力场奇异问题,能更好地描述复合材料基体开裂部位受力的强弱程度,具有使用简单、方便,效率高的优点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610527455.8A CN106202691A (zh) | 2016-07-06 | 2016-07-06 | 一种预报复合材料基体开裂的失效判据 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610527455.8A CN106202691A (zh) | 2016-07-06 | 2016-07-06 | 一种预报复合材料基体开裂的失效判据 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106202691A true CN106202691A (zh) | 2016-12-07 |
Family
ID=57465939
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610527455.8A Pending CN106202691A (zh) | 2016-07-06 | 2016-07-06 | 一种预报复合材料基体开裂的失效判据 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106202691A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108549743A (zh) * | 2018-03-13 | 2018-09-18 | 东南大学 | 一种纤维增强复合材料动态拉伸失效评估方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102364489A (zh) * | 2011-10-25 | 2012-02-29 | 陈志勇 | 木材复杂各向异性本构关系模型的数值模拟方法 |
CN102507744A (zh) * | 2011-10-16 | 2012-06-20 | 浙江大学 | 一种检测碳纤维复合材料破坏失效的声发射装置及方法 |
CN103592175A (zh) * | 2013-11-13 | 2014-02-19 | 北京航空航天大学 | 一种基于连续损伤力学退化模型的复合材料结构失效分析方法 |
US20140067285A1 (en) * | 2011-02-25 | 2014-03-06 | Board Of Supervisors Of Louisiana State University And Agricultural And Mechanical College | Fatigue Monitoring for Composite Materials |
CN105512394A (zh) * | 2015-12-08 | 2016-04-20 | 东南大学 | 一种梁柱式木结构框架整体抗火性能的分析系统及分析方法 |
-
2016
- 2016-07-06 CN CN201610527455.8A patent/CN106202691A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140067285A1 (en) * | 2011-02-25 | 2014-03-06 | Board Of Supervisors Of Louisiana State University And Agricultural And Mechanical College | Fatigue Monitoring for Composite Materials |
CN102507744A (zh) * | 2011-10-16 | 2012-06-20 | 浙江大学 | 一种检测碳纤维复合材料破坏失效的声发射装置及方法 |
CN102364489A (zh) * | 2011-10-25 | 2012-02-29 | 陈志勇 | 木材复杂各向异性本构关系模型的数值模拟方法 |
CN103592175A (zh) * | 2013-11-13 | 2014-02-19 | 北京航空航天大学 | 一种基于连续损伤力学退化模型的复合材料结构失效分析方法 |
CN105512394A (zh) * | 2015-12-08 | 2016-04-20 | 东南大学 | 一种梁柱式木结构框架整体抗火性能的分析系统及分析方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108549743A (zh) * | 2018-03-13 | 2018-09-18 | 东南大学 | 一种纤维增强复合材料动态拉伸失效评估方法 |
CN108549743B (zh) * | 2018-03-13 | 2019-03-12 | 东南大学 | 一种纤维增强复合材料动态拉伸失效评估方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Rong et al. | Effects of geometric configurations of corrugated cores on the local impact and planar compression of sandwich panels | |
Zhu et al. | Low-velocity impact response of multilayer foam core sandwich panels with composite face sheets | |
CN108549743B (zh) | 一种纤维增强复合材料动态拉伸失效评估方法 | |
CN106202598B (zh) | 一种复合材料冲击损伤后剩余压缩强度的分析方法 | |
Li et al. | Low-velocity impact and compression-after-impact behaviour of flax fibre-reinforced composites | |
Fajrin et al. | Flexural behaviour of hybrid sandwich panel with natural fiber composites as the intermediate layer | |
Kormaníková et al. | Finite element analysis of damage modeling of fiber reinforced laminate plate | |
Daniel | Influence of core properties on the failure of composite sandwich beams | |
CN106202691A (zh) | 一种预报复合材料基体开裂的失效判据 | |
Zhang et al. | Damage behaviors of foam sandwiched composite materials under quasi-static three-point bending | |
Zhifei et al. | Interfacial debonding of coated-fiber-reinforced composites under tension-tension cyclic loading | |
Joseph et al. | The est model for predicting progressive damage and failure of open hole bending specimens | |
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 | |
Harvey et al. | Modeling of delamination propagation in composite laminated beam structures | |
Tran et al. | Blast resistance of hybrid elastomeric composite panels | |
Kim et al. | Analysis on the Impact of Composite by Using FEM | |
Watanabe et al. | New Solution for Material Damage Characterization of CFRP Laminate with Filament Winding Structure Using a Hexagonal-Shaped Mandrel | |
Upadhyay et al. | Modeling snow slab failure in propagation saw test using Drucker-Prager model | |
Bang et al. | Simulation analysis on the compression property of sandwich composite | |
Mei et al. | The mechanical performance analysis of grid structures with foam under the tests of quasi-static and dynamic loadings | |
Found et al. | Modelling the impact behaviour of thin CERP panels | |
Reis et al. | Mechanical behaviour of sandwich beams manufactured with glass or jute fiber in facings and cork agglomerates as core | |
MASSARWA et al. | Multiscale progressive damage analyses for fiber reinforced composites subjected to biaxial loading | |
Wang et al. | Buckling response of metallic foam composite sandwich panels subjected to edgewise compression | |
Psarras et al. | Performance of composite plates after multi-site impacts |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20161207 |
|
WD01 | Invention patent application deemed withdrawn after publication |