CN105760605A - 复杂编织结构陶瓷基复合材料疲劳寿命预测方法 - Google Patents
复杂编织结构陶瓷基复合材料疲劳寿命预测方法 Download PDFInfo
- Publication number
- CN105760605A CN105760605A CN201610097455.9A CN201610097455A CN105760605A CN 105760605 A CN105760605 A CN 105760605A CN 201610097455 A CN201610097455 A CN 201610097455A CN 105760605 A CN105760605 A CN 105760605A
- Authority
- CN
- China
- Prior art keywords
- fiber
- strain
- interface
- tau
- unloading
- 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/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/12—Cloth
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)
Abstract
本发明公开了一种复杂编织结构陶瓷基复合材料疲劳寿命预测方法,首先计算该循环数下疲劳性能,并计算该循环数下纤维失效百分数,然后确定纤维失效百分数与纤维失效临界值的关系,计算单胞尺度疲劳性能,得到该循环下的最大应变再确定与最大失效应变εmax的关系,最后得到材料的疲劳寿命曲线。本发明提出了考虑纤维/基体/孔隙的微观尺度模型和考虑经纱/纬纱/孔洞的单胞多尺度预测模型,该模型克服了细观力学法难以直接应用于结构复杂的编织材料和宏观唯象法依赖大量试验且只能预测特定材料疲劳寿命的缺陷,将宏、细观相结合,给出了复杂编织结构的细观应力应变场,在精确地预测出该材料疲劳寿命曲线的同时扩展了材料的应用范围。
Description
技术领域
本发明涉及一种基于多尺度模型的复杂编织结构陶瓷基复合材料疲劳寿命预测方法。
背景技术
复杂编织陶瓷基复合材料具有高比强、高比模、耐高温、耐腐蚀和低密度等优良性能,在空天飞行器高温防护***具有广泛的需求。材料在其使用过程中,由于受载荷和环境等因素的影响,会逐渐产生构件损伤以至于破坏,其主要破坏形式之一是疲劳损伤。这种疲劳破坏对飞行器的危害巨大,是飞行器防护***结构中常见的失效形式。因此,准确地预测复杂编织结构陶瓷基复合材料疲劳寿命,是排除事故隐患、制定健康修复计划和延长使用寿命的保障。
由于单向陶瓷基复合材料存在非纤维方向力学性能弱等缺点,其应用范围受到了限制。二维,2.5维和三维复杂编织结构陶瓷基复合材料的出现,克服了单向复合材料的缺点,同时在厚度方向上纤维束整体化更高,增加了材料层间剪切强度,减少了分层现象,并提高了复合材料抗冲击性能和弯曲疲劳性能,因此大大扩展了陶瓷基复合材料的应用范围。
然而由于复杂编织结构陶瓷基复合材料是一种新型结构材料,国内外还没有高效的方法预测其疲劳寿命,也未见公开的发明专利。孙志刚(孙志刚,许仁红,宋迎东.陶瓷基复合材料低循环拉—拉疲劳寿命预测[J].机械工程学报,2012,12:31-36.)采用细观力学法对正交铺设陶瓷基复合材料低循环拉—拉疲劳寿命进行了预测,得到了材料S-N曲线。该论文计算结果与试验结果较为一致,但是并没有给出复杂编织结构复合材料的疲劳寿命回线。李龙彪(李龙彪.长纤维增强陶瓷基复合材料疲劳损伤模型与寿命预测[D].南京航空航天大学,2010.)采用刚度平均法将二维编织结构简化为层合板,再利用细观力学法对材料疲劳寿命进行了预测。然而该方法并没有实质性的创新且组分材料应力场不能直接从解析表达式中得到,因此这种方法不能够精确地预测出复杂编织结构的疲劳寿命曲线。当前,如何准确的预测复杂编织结构陶瓷基复合材料的疲劳寿命曲线是本技术领域重要而难以解决的问题。
发明内容
针对上述现有技术的不足,本发明的目的是提供一种复杂编织结构陶瓷基复合材料疲劳寿命预测方法,以解决现有技术存在的不能够精确地预测出复杂编织结构的疲劳寿命曲线的问题。
为实现上述目的,本发明采用以下技术方案:
一种复杂编织结构陶瓷基复合材料疲劳寿命预测方法,包括如下步骤:
(1)建立复杂编织结构陶瓷基复合材料的单胞模型;
(2)初始加载至疲劳峰值载荷,认定纬纱单元都已经失效;
(3)进入卸载-加载过程,假设当前循环数为Cycle;
(4)计算该循环数下微观尺度疲劳性能;
(5)计算该循环下纤维失效百分数P(T);
(6)判断纤维失效百分数P(T)与纤维失效临界值q*的关系;若P(T)>q*,进入步骤(9);否则,进入步骤(7);
(7)计算单胞尺度疲劳性能,得到Cycle循环下的最大应变
(8)判断与最大失效应变εmax的关系;若材料发生疲劳破坏,进入步骤(8);否则,Cycle=Cycle+1,同时进入步骤(3);
(9)输出循环次数Cycle,疲劳失效,程序终止。
所述步骤(1)中,通过ANSYS建模仿真软件建立单胞模型。
所述步骤(2)中,失效的纬纱单元还具备承载能力,即单元常数不为零,只是弹性常数发生缩减,缩减后的弹性常数取初始值的1%。
所述步骤(3)中,假定当前循环数Cycle=1。
所述步骤(4)中,当复合材料出现损伤时,假设复合材料应变等于未损伤的纤维应变,即
式中,εc表示复合材料应变,Ef表示纤维弹性模量,L表示基体裂纹间距,σf(x)表示纤维轴向应力,αc表示复合材料热膨胀系数,αf表示纤维热膨胀系数,ΔT表示复合材料制备温度与工作温度之差;
分别考虑基于卸载/重新加载时纤维相对基体在界面脱粘区内滑移和反向滑移的损伤机制,分析以下四种情况下应力应变关系:
(a)界面部分脱粘,卸载/重新加载时纤维相对基体在界面脱粘区完全滑移;
(b)界面部分脱粘,卸载/重新加载时纤维相对基体在界面脱粘区部分滑移;
(c)界面完全脱粘,卸载/重新加载时纤维相对基体在界面脱粘区部分滑移;
(d)界面完全脱粘,卸载/重新加载时纤维相对基体在界面脱粘区完全滑移;
当界面部分脱粘时,卸载/重新加载界面滑移时卸载应变为:
当界面部分脱粘时,卸载/重新加载界面滑移时重新加载应变为:
当界面完全脱粘时,卸载/重新加载界面滑移时卸载应变为:
当界面完全脱粘时,卸载/重新加载界面滑移时重新加载应变为:
上述四式中,εc_pu表示界面部分脱粘时的卸载应变,εc_pr界面部分脱粘时的重新加载应变,εc_fu表示界面完全脱粘时的卸载应变,εc_fr表示界面完全脱粘时的重新加载应变,σ表示复合材料轴向应力,Vf表示纤维体积含量,τi表示界面剪应力,rf表示纤维半径,y表示卸载界面反向滑移长度,z表示重新加载新界面滑移长度,Ld表示界面脱粘长度。
所述步骤(5)中,纤维失效时,假设完好纤维和断裂纤维之间载荷分配服从总体承担准则:
上式中,T为完好纤维在基体裂纹平面承担的应力,<Tb>为断裂纤维承担的应力,P(T)为纤维失效概率:
上式中,σc为纤维在特征长度δc内特征强度,σo(N)为第N次循环时纤维参考强度;τi(N)为第N次循环下纤维/基体界面剪应力;
上式中,lo为参考长度,σo为纤维在参考长度lo内参考强度,mf为纤维威布尔模量。
所述步骤(6)中,纤维失效体积分数临界值为:
上式中,mf为纤维威布尔模量。
所述步骤(7)中,先给定计算载荷步数N,设定当前载荷步Loop=1,再提取N次循环下单胞的最大应变。
所述步骤(8)中,最大应变失效准则认为当材料应变达到最大值时,复合材料发生失效,即:
其中,为允许的最大应变;
当上式成立时,材料发生疲劳失效。
所述步骤(9)中,如此反复循环,最终得到材料的疲劳寿命曲线。
本发明的有益效果是:本发明的方法考虑了基体开裂、纤维断裂、界面滑移和界面磨损等失效机理。将界面剪应力衰退模型、纤维强度衰退模型、纤维随机失效模型和最大应变失效模型相结合,提出了考虑纤维/基体/孔隙的微观尺度模型和考虑经纱/纬纱/孔洞的单胞尺度预测模型。本发明提出的多尺度模型可以给出复杂编织结构的细观应力应变场,因此能够精确地预测出复杂编织结构陶瓷基复合材料的疲劳寿命。
附图说明
图1是2.5维单胞模型结构尺寸参数;
图2是浅交弯联2.5维C/SiC编织复合材料完整的单胞模型和网格划分结果;
图3是加卸载时纤维相对基体滑移示意图;
图4是疲劳寿命预测流程图;
图5是室温环境下2.5D-C/SiC复合材料在拉-拉疲劳载荷作用下S-N理论预测值与试验值对比曲线;
图6是800℃空气环境下2.5D-C/SiC复合材料在拉-拉疲劳载荷作用下S-N理论预测值与试验值对比曲线。
具体实施方式
下面结合附图对本发明作更进一步的说明。
如图4所示,本发明的基于多尺度的复杂编织结构陶瓷基复合材料疲劳寿命预测方法,包括如下步骤:
(1)用ANSYS建模仿真软件建立复杂编织结构陶瓷基复合材料的单胞模型;
(2)初始加载至疲劳峰值载荷,认定纬纱单元都已经失效,纬纱单元失效并不意味着该处单元成为一个“孔洞”,它还具备一定承载能力,即单元常数不为零,只是弹性常数发生缩减,缩减后的弹性常数取初始值的1%;
(3)进入卸载-加载过程,假设当前循环数为Cycle,假定当前循环数Cycle=1;
(4)计算该循环数下微观尺度疲劳性能;
当复合材料出现损伤时,假设复合材料应变等于未损伤的纤维应变,即
式中,εc表示复合材料应变,Ef表示纤维弹性模量,L表示基体裂纹间距,σf(x)表示纤维轴向应力,αc表示复合材料热膨胀系数,αf表示纤维热膨胀系数,ΔT表示复合材料制备温度与工作温度之差;
分别考虑基于卸载/重新加载时纤维相对基体在界面脱粘区内滑移和反向滑移的损伤机制,分析(1)界面部分脱粘,卸载/重新加载时纤维相对基体在界面脱粘区完全滑移;(2)界面部分脱粘,卸载/重新加载时纤维相对基体在界面脱粘区部分滑移;(3)界面完全脱粘,卸载/重新加载时纤维相对基体在界面脱粘区部分滑移;(4)界面完全脱粘,卸载/重新加载时纤维相对基体在界面脱粘区完全滑移四种情况下应力应变关系;
(a)当界面部分脱粘时,卸载/重新加载界面滑移时卸载应变为:
(b)当界面部分脱粘时,卸载/重新加载界面滑移时重新加载应变为:
(c)当界面完全脱粘时,卸载/重新加载界面滑移时卸载应变为:
(d)当界面完全脱粘时,卸载/重新加载界面滑移时重新加载应变为:
上述四式中,εc_pu表示界面部分脱粘时的卸载应变,εc_pr界面部分脱粘时的重新加载应变,εc_fu表示界面完全脱粘时的卸载应变,εc_fr表示界面完全脱粘时的重新加载应变,σ表示复合材料轴向应力,Vf表示纤维体积含量,τi表示界面剪应力,rf表示纤维半径,y表示卸载界面反向滑移长度,z表示重新加载新界面滑移长度,Ld表示界面脱粘长度。
(5)计算该循环下纤维失效百分数P(T);
纤维失效时,假设完好纤维和断裂纤维之间载荷分配服从总体承担准则:
上式中,T为完好纤维在基体裂纹平面承担的应力,<Tb>为断裂纤维承担的应力,P(T)为纤维失效概率:
上式中,σc为纤维在特征长度δc内特征强度,σo(N)为第N次循环时纤维参考强度;τi(N)为第N次循环下纤维/基体界面剪应力。
上式中,lo为参考长度,σo为纤维在参考长度lo内参考强度,mf为纤维威布尔模量。
(6)判断纤维失效百分数P(T)与纤维失效临界值q*的关系;若P(T)>q*,进入步骤(9);否则,进入步骤(7);
纤维失效体积分数临界值为:
上式中,mf为纤维威布尔模量,本发明中取mf=5;
(7)计算单胞尺度疲劳性能,得到Cycle循环下的最大应变
先给定计算载荷步数N,设定当前载荷步Loop=1,再提取N次循环下单胞的最大应变;
(8)判断与最大失效应变εmax的关系,若材料发生疲劳破坏,进入步骤(8);否则,Cycle=Cycle+1,同时进入步骤(3);
最大应变失效准则认为当材料应变达到最大值时,复合材料发生失效,即:
其中,为允许的最大应变;
当上式成立时,材料发生疲劳失效;
(9)输出循环次数Cycle,疲劳失效,程序终止;如此反复循环,最终得到材料的疲劳寿命曲线。
下面结合附图和具体实施例对本发明进行详细说明。
实施例:
本发明以2.5维C/SiC复合材料为例,对室温环境下2.5D-C/SiC复合材料在拉-拉疲劳载荷作用下S-N曲线进行了预测。表1给出了室温纤维束基本材料属性。
表1室温下C/SiC纤维束材料基本属性
(1)基于有限元法,建立2.5维单胞模型。根据周期性特点,选取单胞模型结构尺寸如图1所示。该单胞模型共5个参数。Js为经纱跨度;Jh为经纱高度;Jw为经纱宽度;Ww为纬纱宽度;Wh为纬纱高度,尺寸参数通过显微照片测得。通过ANSYS商用软件,得到了2.5维编织C/SiC复合材料细观结构模型,编织结构几何尺寸见表2。在进行单胞模型有限元分析过程中,本发明选用的是ANSYS软件包中的Solid45三维实体单元对单胞模型进行网格划分,该单元有8个节点,且每个节点具有三个方向自由度:即沿X、Y和Z三个方向平动自由度。图2给出了浅交弯联2.5维C/SiC编织复合材料完整的单胞模型和网格划分结果,该单胞共有1408个单元,2005个节点。
表2单胞模型几何尺寸
(2)初始加载至疲劳峰值载荷,认定纬纱单元都已经失效。纬纱单元失效并不意味着该处单元成为一个“孔洞”,它还具备一定承载能力,即单元常数不为零,只是弹性常数发生缩减,缩减后的弹性常数取初始值的1%。
(3)进入卸载-加载过程,假设当前循环数为Cycle。假定当前循环数Cycle=1。
(4)计算该循环数下微观尺度疲劳性能。当复合材料出现损伤时,假设复合材料应变等于未损伤的纤维应变。即
式中,εc表示复合材料应变,Ef表示纤维弹性模量,L表示基体裂纹间距,σf(x)表示纤维轴向应力,αc表示复合材料热膨胀系数,αf表示纤维热膨胀系数,ΔT表示复合材料制备温度与工作温度之差。
分别考虑了基于卸载/重新加载时纤维相对基体在界面脱粘区内滑移和反向滑移的损伤机制,分析了(1)界面部分脱粘,卸载/重新加载时纤维相对基体在界面脱粘区完全滑移;(2)界面部分脱粘,卸载/重新加载时纤维相对基体在界面脱粘区部分滑移;(3)界面完全脱粘,卸载/重新加载时纤维相对基体在界面脱粘区部分滑移;(4)界面完全脱粘,卸载/重新加载时纤维相对基体在界面脱粘区完全滑移四种情况下应力应变关系。图3给出了纤维相对基体滑移的过程。
(a)当界面部分脱粘时,卸载/重新加载界面滑移时卸载应变为:
(b)当界面部分脱粘时,卸载/重新加载界面滑移时重新加载应变为:
(c)当界面完全脱粘时,卸载/重新加载界面滑移时卸载应变为:
(d)当界面完全脱粘时,卸载/重新加载界面滑移时重新加载应变为:
上述四式中,εc_pu表示界面部分脱粘时的卸载应变,εc_pr界面部分脱粘时的重新加载应变,εc_fu表示界面完全脱粘时的卸载应变,εc_fr表示界面完全脱粘时的重新加载应变,σ表示复合材料轴向应力,Vf表示纤维体积含量,τi表示界面剪应力,rf表示纤维半径,y表示卸载界面反向滑移长度,z表示重新加载新界面滑移长度,Ld表示界面脱粘长度。
(5)计算该循环下纤维失效百分数P(T)。纤维失效时,假设完好纤维和断裂纤维之间载荷分配服从总体承担准则:
上式中,T为完好纤维在基体裂纹平面承担的应力,<Tb>为断裂纤维承担的应力,P(T)为纤维失效概率:
上式中,σc为纤维在特征长度δc内特征强度,σo(N)为第N次循环时纤维参考强度;τi(N)为第N次循环下纤维/基体界面剪应力。
上式中,lo为参考长度,σo为纤维在参考长度lo内参考强度,mf为纤维威布尔模量。
当纤维断裂时,断裂处纤维不在承载。远离断裂点处,通过纤维/基体间界面剪应力传递载荷,使得纤维承担的载荷逐渐增加。
定义断裂纤维承担应力恢复到其断裂前应力所需要的界面滑移长度为:
在距离基体裂纹平面±lf范围内,纤维断裂的概率密度函数f(x)为:
式中:σo(N)为第N次循环时纤维参考强度;τi(N)为第N次循环下纤维/基体界面剪应力。
联立以上各式,得到断裂纤维承担的平均应力〈Tb〉为:
联立以上各式,得:
(6)判断纤维失效百分数P(T)与纤维失效临界值q*的关系。若P(T)>q*,进入步骤(9);否则,进入步骤(7)。纤维失效体积分数临界值为:
上式中,mf为纤维威布尔模量,本文取mf=5。
(7)计算单胞尺度疲劳性能,得到Cycle循环下的最大应变先给定计算载荷步数N,设定当前载荷步Loop=1,再提取N次循环下单胞的最大应变,具体步骤见图4。
(8)判断与最大失效应变εmax的关系。若材料发生疲劳破坏,进入步骤(8);否则,Cycle=Cycle+1,同时进入步骤(3)。最大应变失效准则认为当材料应变达到某一值时,复合材料发生失效。即:
当上式成立时,材料发生疲劳失效。
(9)输出循环次数Cycle,疲劳失效,程序终止。如此反复循环,最终得到材料的疲劳寿命曲线。图5给出了室温环境下2.5D-C/SiC复合材料在拉-拉疲劳载荷作用下S-N理论预测值与试验值对比曲线。图6给出了800℃空气环境下2.5D-C/SiC复合材料在拉-拉疲劳载荷作用下S-N理论预测值与试验值对比曲线。
通过对比可见,本发明的方法能够有效预测复杂编织结构陶瓷基复合材料疲劳寿命。
虽然本发明已以较佳实施例公开如上,但他们并不是用来限定本发明的,任何熟习此技艺者,在不脱离本发明之精神和范围内,自当可做各种变化或润饰,因此本发明的保护范围应当已本申请的专利保护范围所界定的为准。
Claims (10)
1.一种复杂编织结构陶瓷基复合材料疲劳寿命预测方法,其特征在于:包括如下步骤:
(1)建立复杂编织结构陶瓷基复合材料的单胞模型;
(2)初始加载至疲劳峰值载荷,认定纬纱单元都已经失效;
(3)进入卸载-加载过程,假设当前循环数为Cycle;
(4)计算该循环数下微观尺度疲劳性能;
(5)计算该循环下纤维失效百分数P(T);
(6)判断纤维失效百分数P(T)与纤维失效临界值q*的关系;若P(T)>q*,进入步骤(9);否则,进入步骤(7);
(7)计算单胞尺度疲劳性能,得到Cycle循环下的最大应变
(8)判断与最大失效应变εmax的关系;若材料发生疲劳破坏,进入步骤(8);否则,Cycle=Cycle+1,同时进入步骤(3);
(9)输出循环次数Cycle,疲劳失效,程序终止。
2.根据权利要求1所述的复杂编织结构陶瓷基复合材料疲劳寿命预测方法,其特征在于:所述步骤(1)中,通过ANSYS建模仿真软件建立单胞模型。
3.根据权利要求1所述的复杂编织结构陶瓷基复合材料疲劳寿命预测方法,其特征在于:所述步骤(2)中,失效的纬纱单元还具备承载能力,即单元常数不为零,只是弹性常数发生缩减,缩减后的弹性常数取初始值的1%。
4.根据权利要求1所述的复杂编织结构陶瓷基复合材料疲劳寿命预测方法,其特征在于:所述步骤(3)中,假定当前循环数Cycle=1。
5.根据权利要求1所述的复杂编织结构陶瓷基复合材料疲劳寿命预测方法,其特征在于:所述步骤(4)中,当复合材料出现损伤时,假设复合材料应变等于未损伤的纤维应变,即
式中,εc表示复合材料应变,Ef表示纤维弹性模量,L表示基体裂纹间距,σf(x)表示纤维轴向应力,αc表示复合材料热膨胀系数,αf表示纤维热膨胀系数,△T表示复合材料制备温度与工作温度之差;
分别考虑基于卸载/重新加载时纤维相对基体在界面脱粘区内滑移和反向滑移的损伤机制,分析以下四种情况下应力应变关系:
(a)界面部分脱粘,卸载/重新加载时纤维相对基体在界面脱粘区完全滑移;
(b)界面部分脱粘,卸载/重新加载时纤维相对基体在界面脱粘区部分滑移;
(c)界面完全脱粘,卸载/重新加载时纤维相对基体在界面脱粘区部分滑移;
(d)界面完全脱粘,卸载/重新加载时纤维相对基体在界面脱粘区完全滑移;
当界面部分脱粘时,卸载/重新加载界面滑移时卸载应变为:
当界面部分脱粘时,卸载/重新加载界面滑移时重新加载应变为:
当界面完全脱粘时,卸载/重新加载界面滑移时卸载应变为:
当界面完全脱粘时,卸载/重新加载界面滑移时重新加载应变为:
上述四式中,εc_pu表示界面部分脱粘时的卸载应变,εc_pr界面部分脱粘时的重新加载应变,εc_fu表示界面完全脱粘时的卸载应变,εc_fr表示界面完全脱粘时的重新加载应变,σ表示复合材料轴向应力,Vf表示纤维体积含量,τi表示界面剪应力,rf表示纤维半径,y表示卸载界面反向滑移长度,z表示重新加载新界面滑移长度,Ld表示界面脱粘长度。
6.根据权利要求1所述的复杂编织结构陶瓷基复合材料疲劳寿命预测方法,其特征在于:所述步骤(5)中,纤维失效时,假设完好纤维和断裂纤维之间载荷分配服从总体承担准则:
上式中,T为完好纤维在基体裂纹平面承担的应力,<Tb>为断裂纤维承担的应力,P(T)为纤维失效概率:
上式中,σc为纤维在特征长度δc内特征强度,σo(N)为第N次循环时纤维参考强度;τi(N)为第N次循环下纤维/基体界面剪应力;
上式中,lo为参考长度,σo为纤维在参考长度lo内参考强度,mf为纤维威布尔模量。
7.根据权利要求1所述的复杂编织结构陶瓷基复合材料疲劳寿命预测方法,其特征在于:
所述步骤(6)中,纤维失效体积分数临界值为:
上式中,mf为纤维威布尔模量。
8.根据权利要求1所述的复杂编织结构陶瓷基复合材料疲劳寿命预测方法,其特征在于:所述步骤(7)中,先给定计算载荷步数N,设定当前载荷步Loop=1,再提取N次循环下单胞的最大应变。
9.根据权利要求1所述的复杂编织结构陶瓷基复合材料疲劳寿命预测方法,其特征在于:所述步骤(8)中,最大应变失效准则认为当材料应变达到最大值时,复合材料发生失效,即:
其中,为允许的最大应变;
当上式成立时,材料发生疲劳失效。
10.根据权利要求1所述的复杂编织结构陶瓷基复合材料疲劳寿命预测方法,其特征在于:所述步骤(9)中,如此反复循环,最终得到材料的疲劳寿命曲线。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2015109566207 | 2015-12-17 | ||
CN201510956620 | 2015-12-17 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105760605A true CN105760605A (zh) | 2016-07-13 |
Family
ID=56331133
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610097455.9A Pending CN105760605A (zh) | 2015-12-17 | 2016-02-23 | 复杂编织结构陶瓷基复合材料疲劳寿命预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105760605A (zh) |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106294967A (zh) * | 2016-08-02 | 2017-01-04 | 浙江大学 | 一种考虑荷载频率的水泥基材料疲劳失效概率和疲劳应变概率模型的建立方法 |
CN106570288A (zh) * | 2016-11-09 | 2017-04-19 | 中国科学院深圳先进技术研究院 | 编织结构的建模方法及装置 |
CN107301285A (zh) * | 2017-06-16 | 2017-10-27 | 南京航空航天大学 | 基于剩余寿命预测的非电子产品序贯验证试验方法 |
CN107747919A (zh) * | 2017-10-13 | 2018-03-02 | 广东工业大学 | 一种纤维横截面积测量方法及*** |
CN108204925A (zh) * | 2016-12-16 | 2018-06-26 | 海口未来技术研究院 | 复合材料的疲劳寿命预测方法及预测*** |
CN108427850A (zh) * | 2018-03-20 | 2018-08-21 | 东南大学 | 一种复合材料剪切强度包络线的获取方法 |
CN108984981A (zh) * | 2018-09-04 | 2018-12-11 | 北京航空航天大学 | 一种模拟陶瓷基复合材料失效机理的统一代表性体积元 |
CN109002649A (zh) * | 2018-08-20 | 2018-12-14 | 东汉新能源汽车技术有限公司 | 碳纤维增强复合材料构件的疲劳寿命预测方法 |
CN109583052A (zh) * | 2018-11-15 | 2019-04-05 | 西南交通大学 | 纤维增强陶瓷基复合材料微结构精确建模及模量计算方法 |
CN109598098A (zh) * | 2019-01-23 | 2019-04-09 | 南京航空航天大学 | 一种考虑环境影响的编织陶瓷基复合材料稳态基体开裂应力的预测方法 |
CN109614755A (zh) * | 2018-12-29 | 2019-04-12 | 南京航空航天大学 | 一种通过迟滞耗散能预测编织陶瓷基复合材料高温疲劳纤维/基体界面剪应力的方法 |
CN109632887A (zh) * | 2018-12-29 | 2019-04-16 | 南京航空航天大学 | 一种通过迟滞耗散能监测编织陶瓷基复合材料高温疲劳损伤累积的方法 |
CN109632530A (zh) * | 2018-12-29 | 2019-04-16 | 南京航空航天大学 | 一种编织陶瓷基复合材料热机械疲劳迟滞回线的预测方法 |
CN109657412A (zh) * | 2019-01-22 | 2019-04-19 | 北京航空航天大学 | 一种陶瓷基复合材料的疲劳寿命的多尺度预测方法 |
CN109670272A (zh) * | 2019-01-23 | 2019-04-23 | 南京航空航天大学 | 一种编织陶瓷基复合材料过载疲劳迟滞回线的预测方法 |
CN109684598A (zh) * | 2018-12-29 | 2019-04-26 | 南京航空航天大学 | 一种通过迟滞耗散能预测编织陶瓷基复合材料高温疲劳寿命的方法 |
CN109840392A (zh) * | 2019-03-26 | 2019-06-04 | 南京航空航天大学 | 一种陶瓷基复合材料疲劳极限的快速预测方法 |
CN109920495A (zh) * | 2019-03-28 | 2019-06-21 | 南京航空航天大学 | 一种编织陶瓷基复合材料强度的多尺度预测方法 |
CN109992850A (zh) * | 2019-03-15 | 2019-07-09 | 南京航空航天大学 | 一种陶瓷基复合材料应力氧化环境剩余拉伸强度预测方法 |
CN110348058A (zh) * | 2019-06-12 | 2019-10-18 | 南京航空航天大学 | 一种陶瓷基复合材料在高温应力环境下的剩余强度计算方法 |
CN110633556A (zh) * | 2019-10-24 | 2019-12-31 | 南京航空航天大学 | 一种陶瓷基复合材料流固耦合响应计算方法 |
CN110688790A (zh) * | 2019-08-30 | 2020-01-14 | 南京航空航天大学 | 一种复杂编织结构陶瓷基复合材料拉伸失效模拟方法 |
CN111198141A (zh) * | 2020-03-11 | 2020-05-26 | 烟台宏远氧业股份有限公司 | 基于细观力学的聚合物基复合材料疲劳寿命预测方法 |
CN111400922A (zh) * | 2020-03-23 | 2020-07-10 | 南京航空航天大学 | 单向陶瓷基复合材料任意应变加卸载应力应变行为计算方法 |
CN113139300A (zh) * | 2021-05-14 | 2021-07-20 | 南京工程学院 | 陶瓷微孔致裂强度定量预测方法、装置、计算机设备及存储介质 |
CN114692468A (zh) * | 2022-04-24 | 2022-07-01 | 西南交通大学 | 一种连续纤维增强陶瓷基复合材料跨尺度损伤预测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6361888B1 (en) * | 1999-01-19 | 2002-03-26 | The Board Of Trustees Of The University Of Illinois | Toughening of ceramic composites by transformation weakening of interphases |
CN103871059A (zh) * | 2014-03-13 | 2014-06-18 | 南京航空航天大学 | 纤维增强复合材料等效弹性参数的计算方法 |
CN103913378A (zh) * | 2014-03-21 | 2014-07-09 | 南京航空航天大学 | 陶瓷基复合材料拉伸应力应变曲线的测试方法 |
CN104866690A (zh) * | 2015-06-16 | 2015-08-26 | 南京航空航天大学 | 单向陶瓷基复合材料任意加卸载应力应变行为预测方法 |
-
2016
- 2016-02-23 CN CN201610097455.9A patent/CN105760605A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6361888B1 (en) * | 1999-01-19 | 2002-03-26 | The Board Of Trustees Of The University Of Illinois | Toughening of ceramic composites by transformation weakening of interphases |
CN103871059A (zh) * | 2014-03-13 | 2014-06-18 | 南京航空航天大学 | 纤维增强复合材料等效弹性参数的计算方法 |
CN103913378A (zh) * | 2014-03-21 | 2014-07-09 | 南京航空航天大学 | 陶瓷基复合材料拉伸应力应变曲线的测试方法 |
CN104866690A (zh) * | 2015-06-16 | 2015-08-26 | 南京航空航天大学 | 单向陶瓷基复合材料任意加卸载应力应变行为预测方法 |
Non-Patent Citations (2)
Title |
---|
C.D.M.LILJEDAHL 等: "Residual stresses in structures reinforced with adhesively bonded straps designed to retard fatigue crack growth", 《COMPOSITE STRUCTURES》 * |
杨福树: "2.5维编织陶瓷基复合材料疲劳行为研究", 《中国优秀硕士学位论文全文数据库工程科技I辑》 * |
Cited By (40)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106294967B (zh) * | 2016-08-02 | 2019-10-29 | 浙江大学 | 一种考虑荷载频率的水泥基材料疲劳失效概率和疲劳应变概率模型的建立方法 |
CN106294967A (zh) * | 2016-08-02 | 2017-01-04 | 浙江大学 | 一种考虑荷载频率的水泥基材料疲劳失效概率和疲劳应变概率模型的建立方法 |
CN106570288A (zh) * | 2016-11-09 | 2017-04-19 | 中国科学院深圳先进技术研究院 | 编织结构的建模方法及装置 |
CN106570288B (zh) * | 2016-11-09 | 2020-02-21 | 中国科学院深圳先进技术研究院 | 编织结构的建模方法及装置 |
CN108204925A (zh) * | 2016-12-16 | 2018-06-26 | 海口未来技术研究院 | 复合材料的疲劳寿命预测方法及预测*** |
CN107301285B (zh) * | 2017-06-16 | 2020-11-20 | 南京航空航天大学 | 基于剩余寿命预测的非电子产品序贯验证试验方法 |
CN107301285A (zh) * | 2017-06-16 | 2017-10-27 | 南京航空航天大学 | 基于剩余寿命预测的非电子产品序贯验证试验方法 |
CN107747919A (zh) * | 2017-10-13 | 2018-03-02 | 广东工业大学 | 一种纤维横截面积测量方法及*** |
CN108427850A (zh) * | 2018-03-20 | 2018-08-21 | 东南大学 | 一种复合材料剪切强度包络线的获取方法 |
CN109002649A (zh) * | 2018-08-20 | 2018-12-14 | 东汉新能源汽车技术有限公司 | 碳纤维增强复合材料构件的疲劳寿命预测方法 |
CN108984981A (zh) * | 2018-09-04 | 2018-12-11 | 北京航空航天大学 | 一种模拟陶瓷基复合材料失效机理的统一代表性体积元 |
CN108984981B (zh) * | 2018-09-04 | 2022-10-04 | 北京航空航天大学 | 一种模拟陶瓷基复合材料失效机理的统一代表性体积元 |
CN109583052B (zh) * | 2018-11-15 | 2022-04-05 | 西南交通大学 | 纤维增强陶瓷基复合材料微结构精确建模及模量计算方法 |
CN109583052A (zh) * | 2018-11-15 | 2019-04-05 | 西南交通大学 | 纤维增强陶瓷基复合材料微结构精确建模及模量计算方法 |
CN109632530B (zh) * | 2018-12-29 | 2020-09-22 | 南京航空航天大学 | 一种编织陶瓷基复合材料热机械疲劳迟滞回线的预测方法 |
CN109614755A (zh) * | 2018-12-29 | 2019-04-12 | 南京航空航天大学 | 一种通过迟滞耗散能预测编织陶瓷基复合材料高温疲劳纤维/基体界面剪应力的方法 |
CN109684598A (zh) * | 2018-12-29 | 2019-04-26 | 南京航空航天大学 | 一种通过迟滞耗散能预测编织陶瓷基复合材料高温疲劳寿命的方法 |
CN109632530A (zh) * | 2018-12-29 | 2019-04-16 | 南京航空航天大学 | 一种编织陶瓷基复合材料热机械疲劳迟滞回线的预测方法 |
CN109632887A (zh) * | 2018-12-29 | 2019-04-16 | 南京航空航天大学 | 一种通过迟滞耗散能监测编织陶瓷基复合材料高温疲劳损伤累积的方法 |
CN109632887B (zh) * | 2018-12-29 | 2020-08-14 | 南京航空航天大学 | 一种监测编织陶瓷基复合材料高温疲劳损伤累积的方法 |
CN109657412A (zh) * | 2019-01-22 | 2019-04-19 | 北京航空航天大学 | 一种陶瓷基复合材料的疲劳寿命的多尺度预测方法 |
CN109657412B (zh) * | 2019-01-22 | 2022-05-03 | 北京航空航天大学 | 一种陶瓷基复合材料的疲劳寿命的多尺度预测方法 |
CN109670272A (zh) * | 2019-01-23 | 2019-04-23 | 南京航空航天大学 | 一种编织陶瓷基复合材料过载疲劳迟滞回线的预测方法 |
CN109598098A (zh) * | 2019-01-23 | 2019-04-09 | 南京航空航天大学 | 一种考虑环境影响的编织陶瓷基复合材料稳态基体开裂应力的预测方法 |
CN109670272B (zh) * | 2019-01-23 | 2020-08-14 | 南京航空航天大学 | 一种编织陶瓷基复合材料过载疲劳迟滞回线的预测方法 |
CN109992850A (zh) * | 2019-03-15 | 2019-07-09 | 南京航空航天大学 | 一种陶瓷基复合材料应力氧化环境剩余拉伸强度预测方法 |
CN109840392B (zh) * | 2019-03-26 | 2019-10-11 | 南京航空航天大学 | 一种陶瓷基复合材料疲劳极限的快速预测方法 |
CN109840392A (zh) * | 2019-03-26 | 2019-06-04 | 南京航空航天大学 | 一种陶瓷基复合材料疲劳极限的快速预测方法 |
CN109920495A (zh) * | 2019-03-28 | 2019-06-21 | 南京航空航天大学 | 一种编织陶瓷基复合材料强度的多尺度预测方法 |
CN110348058A (zh) * | 2019-06-12 | 2019-10-18 | 南京航空航天大学 | 一种陶瓷基复合材料在高温应力环境下的剩余强度计算方法 |
CN110688790B (zh) * | 2019-08-30 | 2021-06-22 | 南京航空航天大学 | 一种复杂编织结构陶瓷基复合材料拉伸失效模拟方法 |
CN110688790A (zh) * | 2019-08-30 | 2020-01-14 | 南京航空航天大学 | 一种复杂编织结构陶瓷基复合材料拉伸失效模拟方法 |
CN110633556B (zh) * | 2019-10-24 | 2020-05-26 | 南京航空航天大学 | 一种陶瓷基复合材料流固耦合响应计算方法 |
CN110633556A (zh) * | 2019-10-24 | 2019-12-31 | 南京航空航天大学 | 一种陶瓷基复合材料流固耦合响应计算方法 |
CN111198141A (zh) * | 2020-03-11 | 2020-05-26 | 烟台宏远氧业股份有限公司 | 基于细观力学的聚合物基复合材料疲劳寿命预测方法 |
CN111400922A (zh) * | 2020-03-23 | 2020-07-10 | 南京航空航天大学 | 单向陶瓷基复合材料任意应变加卸载应力应变行为计算方法 |
CN113139300A (zh) * | 2021-05-14 | 2021-07-20 | 南京工程学院 | 陶瓷微孔致裂强度定量预测方法、装置、计算机设备及存储介质 |
CN113139300B (zh) * | 2021-05-14 | 2023-09-12 | 南京工程学院 | 陶瓷微孔致裂强度定量预测方法、装置、计算机设备及存储介质 |
CN114692468A (zh) * | 2022-04-24 | 2022-07-01 | 西南交通大学 | 一种连续纤维增强陶瓷基复合材料跨尺度损伤预测方法 |
CN114692468B (zh) * | 2022-04-24 | 2023-09-19 | 西南交通大学 | 一种连续纤维增强陶瓷基复合材料跨尺度损伤预测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105760605A (zh) | 复杂编织结构陶瓷基复合材料疲劳寿命预测方法 | |
CN105701312A (zh) | 复杂编织结构陶瓷基复合材料疲劳迟滞行为预测方法 | |
Huang et al. | Review of the mechanical properties of a 3D woven composite and its applications | |
Sevenois et al. | Fatigue damage modeling techniques for textile composites: review and comparison with unidirectional composite modeling techniques | |
WO2018107730A1 (zh) | 复合材料的疲劳寿命预测方法及预测*** | |
Tan et al. | Finite element model for compression after impact behaviour of stitched composites | |
CN111024486B (zh) | 一种单向陶瓷基复合材料蠕变行为预测方法 | |
CN104866690A (zh) | 单向陶瓷基复合材料任意加卸载应力应变行为预测方法 | |
Wang et al. | Modelling of damage evolution in braided composites: recent developments | |
Capela et al. | Fatigue strength of tubular carbon fibre composites under bending/torsion loading | |
Hong et al. | Safety evaluation of 70 MPa-capacity type III hydrogen pressure vessel considering material degradation of composites due to temperature rise | |
Yang et al. | Multiscale modeling and failure analysis of an 8-harness satin woven composite | |
Chen et al. | Numerical analysis and experimental investigation of wind turbine blades with innovative features: Structural response and characteristics | |
Dewangan et al. | Multiscale modelling of 3D orthogonal woven composite under ballistic impact using FEM | |
Dewangan et al. | Multi‐scale modeling of three‐dimensional angle interlock woven composite subjected to ballistic impact using FEM | |
Li et al. | Micromechanical finite element analyses of fire-retarded woven fabric composites at elevated temperatures using unit cells at multiple length scales | |
Pankow et al. | Experimental investigation on the deformation response of hybrid 3D woven composites | |
Tomkova et al. | Evaluation of effective thermal conductivities of porous textile composites | |
Jin et al. | Structural effects of three-dimensional angle-interlock woven composite undergoing bending cyclic loading | |
Yang et al. | A progressive damage model for three-dimensional woven composites under ballistic impact loading | |
Ohyama et al. | Estimation of mechanical behavior of braided composites based on mesh superposition method | |
Rokkam et al. | A peridynamics-FEM approach for crack path prediction in fiber-reinforced composites | |
Bhuiyan et al. | A multiscale approach for progressive fatigue failure modeling of a woven composite RVE | |
Sorini et al. | Effects of adiabatic heating on the high strain rate deformation response of triaxially braided polymer matrix composites | |
Tong et al. | Finite element analysis of woven fabric laminates structural strength |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160713 |
|
RJ01 | Rejection of invention patent application after publication |