CN106202691A - 一种预报复合材料基体开裂的失效判据 - Google Patents

一种预报复合材料基体开裂的失效判据 Download PDF

Info

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
Application number
CN201610527455.8A
Other languages
English (en)
Inventor
熊峻江
杜屹森
朱云涛
云新尧
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN201610527455.8A priority Critical patent/CN106202691A/zh
Publication of CN106202691A publication Critical patent/CN106202691A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical 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)下,开裂前基体储存的弹性应变能密度分别为
u 22 = σ 22 2 2 E 22 u 12 = τ 12 2 2 G 12 u 23 = τ 23 2 2 G 23 - - - ( 1 )
其中,下标1、2分别代表复合材料的纵向和横向,3代表复合材料的厚度方向。u22代表基体开裂前的横向拉压弹性应变能密度,u12代表基体开裂前的面内剪切弹性应变能密度,u23代表基体开裂前的面外剪切弹性应变能密度。E22代表基体的横向拉伸和压缩模量,G12代表基体的面内剪切模量,G23代表基体的面外剪切模量。σ22代表基体的横向正应力,τ12和τ23分别代表基体的面内和面外切应力。
步骤三、计算基体开裂时释放的弹性应变能
在复杂应力状态(即σ22、τ12和τ23)下,基体发生开裂失效时释放的弹性应变能密度分别为
u 22 , c r = Y 2 2 E 22 u 12 , c r = S 12 2 2 G 12 u 23 , c r = S 23 2 2 G 23 - - - ( 2 )
其中,u22,cr代表基体开裂时释放的临界横向拉压弹性应变能密度,u12,cr代表基体开裂时释放的临界面内剪切弹性应变能密度,u23,cr代表基体开裂时释放的临界面外剪切弹性应变能密度。Y代表基体的横向拉伸和压缩强度,S12和S23代表基体的面内和面外剪切强度。
步骤四、建立基体开裂的失效判据
定义f22、f12和f23分别为基体横向拉/压失效、面内剪切失效以及面外剪切失效的能量失效指数,其表达式可写为
f 22 = u 22 / u 22 , c r = ( σ 22 Y ) 2 f 2 = u 12 / u 12 , c r = ( τ 12 S 12 ) 2 f 23 = u 23 / u 23 , c r = ( τ 23 S 23 ) 2 - - - ( 3 )
将横向拉压失效、面内剪切失效以及面外剪切失效的能量失效指数相加,得到基体开裂总失效指数,即
F = f 22 + f 12 + f 23 = ( σ 22 Y ) 2 + ( τ 12 S 12 ) 2 + ( τ 23 S 23 ) 2 - - - ( 4 )
基体开裂的失效判据可写为F≥1。
优点及功效:本发明提供的失效判据则可避免裂纹尖端的应力场奇异问题,更好地描述复合材料基体开裂部位受力的强弱程度,具有使用简单、方便,效率高的优点,能为复合材料结构的损伤容限设计提供技术支持。
附图说明
图1为基体的可能破坏模式示意图。
图2为基体开裂失效判据发明方法流程图。
图中符号说明如下:
图1中的σ22为基体的横向正应力分量,τ12和τ23分别为基体的面内和面外切应力分量。
具体实施方式
图2为本发明所述方法的流程框图,本发明分四步实现,具体为:
步骤一、确定复杂应力状态下复合材料基体主要失效模式
复杂应力状态下,复合材料基体的主要失效模式为相邻纤维间的基体开裂。在横向拉/压力以及剪切力作用下,单层板相邻纤维之间的树脂基体可能发生的失效模式包括横向拉伸/压缩开裂和面内/面外剪切开裂(如图1所示)。
步骤二、计算基体开裂前储存的弹性应变能
在复杂应力状态(即σ22、τ12和τ23)下,开裂前基体储存的弹性应变能密度分别为
u 22 = σ 22 2 2 E 22 u 12 = τ 12 2 2 G 12 u 23 = τ 23 2 2 G 23 - - - ( 1 )
其中,下标1、2分别代表复合材料的纵向和横向,3代表复合材料的厚度方向。u22代表基体开裂前的横向拉压弹性应变能密度,u12代表基体开裂前的面内剪切弹性应变能密度,u23代表基体开裂前的面外剪切弹性应变能密度。E22代表基体的横向拉伸和压缩模量,G12代表基体的面内剪切模量,G23代表基体的面外剪切模量。σ22代表基体的横向正应力,τ12和τ23分别代表基体的面内和面外切应力。
步骤三、计算基体开裂时释放的弹性应变能
在复杂应力状态(即σ22、τ12和τ23)下,基体发生开裂失效时释放的弹性应变能密度分别为
u 22 , c r = Y 2 2 E 22 u 12 , c r = S 12 2 2 G 12 u 23 , c r = S 23 2 2 G 23 - - - ( 2 )
其中,u22,cr代表基体开裂时释放的临界横向拉压弹性应变能密度,u12,cr代表基体开裂时释放的临界面内剪切弹性应变能密度,u23,cr代表基体开裂时释放的临界面外剪切弹性应变能密度。Y代表基体的横向拉伸和压缩强度,S12和S23代表基体的面内和面外剪切强度。
步骤四、建立基体开裂的失效判据
定义f22、f12和f23分别为基体横向拉/压失效、面内剪切失效以及面外剪切失效的能量失效指数,其表达式可写为
f 22 = u 22 / u 22 , c r = ( σ 22 Y ) 2 f 2 = u 12 / u 12 , c r = ( τ 12 S 12 ) 2 f 23 = u 23 / u 23 , c r = ( τ 23 S 23 ) 2 - - - ( 3 )
将横向拉压失效、面内剪切失效以及面外剪切失效的能量失效指数相加,得到基体开裂总失效指数,即
F = f 22 + f 12 + f 23 = ( σ 22 Y ) 2 + ( τ 12 S 12 ) 2 + ( τ 23 S 23 ) 2 - - - ( 4 )
基体开裂的失效判据可写为F≥1。

Claims (1)

1.一种预报复合材料基体开裂的失效判据,该判据的具体步骤如下:
步骤一、确定复杂应力状态下复合材料基体主要失效模式
复杂应力状态下,复合材料基体的主要失效模式为相邻纤维间的基体开裂。在横向拉/压力以及剪切力作用下,单层板相邻纤维之间的树脂基体可能发生的失效模式包括横向拉伸/压缩开裂和面内/面外剪切开裂。
步骤二、计算基体开裂前储存的弹性应变能
在复杂应力状态(即σ22、τ12和τ23)下,开裂前基体储存的弹性应变能密度分别为
u 22 = σ 22 2 2 E 22 u 12 = τ 12 2 2 G 12 u 23 = τ 23 2 2 G 23 - - - ( 1 )
其中,下标1、2分别代表复合材料的纵向和横向,3代表复合材料的厚度方向。u22代表基体开裂前的横向拉压弹性应变能密度,u12代表基体开裂前的面内剪切弹性应变能密度,u23代表基体开裂前的面外剪切弹性应变能密度。E22代表基体的横向拉伸和压缩模量,G12代表基体的面内剪切模量,G23代表基体的面外剪切模量。σ22代表基体的横向正应力,τ12和τ23分别代表基体的面内和面外切应力。
步骤三、计算基体开裂时释放的弹性应变能
在复杂应力状态(即σ22、τ12和τ23)下,基体发生开裂失效时释放的弹性应变能密度分别为
u 22 , c r = Y 2 2 E 22 u 12 , c r = S 12 2 2 G 12 u 23 , c r = S 23 2 2 G 23 - - - ( 2 )
其中,u22,cr代表基体开裂时释放的临界横向拉压弹性应变能密度,u12,cr代表基体开裂时释放的临界面内剪切弹性应变能密度,u23,cr代表基体开裂时释放的临界面外剪切弹性应变能密度。Y代表基体的横向拉伸和压缩强度,S12和S23代表基体的面内和面外剪切强度。
步骤四、建立基体开裂的失效判据
定义f22、f12和f23分别为基体横向拉/压失效、面内剪切失效以及面外剪切失效的能量失效指数,其表达式可写为
f 22 = u 22 / u 22 , c r = ( σ 22 Y ) 2 f 12 = u 12 / u 12 , c r = ( τ 12 S 12 ) 2 f 23 = u 23 / u 23 , c r = ( τ 23 S 23 ) 2 - - - ( 3 )
将横向拉压失效、面内剪切失效以及面外剪切失效的能量失效指数相加,得到基体开裂总失效指数,即
F = f 22 + f 12 + f 23 = ( σ 22 Y ) 2 + ( τ 12 S 12 ) 2 + ( τ 23 S 23 ) 2 - - - ( 4 )
基体开裂的失效判据可写为F≥1。
功能特征:避免复合材料基体裂纹尖端的应力场奇异问题,能更好地描述复合材料基体开裂部位受力的强弱程度,具有使用简单、方便,效率高的优点。
CN201610527455.8A 2016-07-06 2016-07-06 一种预报复合材料基体开裂的失效判据 Pending CN106202691A (zh)

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)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108549743A (zh) * 2018-03-13 2018-09-18 东南大学 一种纤维增强复合材料动态拉伸失效评估方法

Citations (5)

* Cited by examiner, † Cited by third party
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 东南大学 一种梁柱式木结构框架整体抗火性能的分析系统及分析方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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