CN110309619B - 一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法 - Google Patents
一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法 Download PDFInfo
- Publication number
- CN110309619B CN110309619B CN201910618452.9A CN201910618452A CN110309619B CN 110309619 B CN110309619 B CN 110309619B CN 201910618452 A CN201910618452 A CN 201910618452A CN 110309619 B CN110309619 B CN 110309619B
- Authority
- CN
- China
- Prior art keywords
- fiber
- stress
- ceramic matrix
- matrix composite
- matrix
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- 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
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法,在陶瓷基复合材料中,由于纤维/基体间界面的载荷传递作用,断裂纤维还具有一定的承载能力。本发明考虑了断裂纤维的继续承载,给出了判断纤维断裂位置的方法,计算了断裂纤维在基体裂纹平面处的应力,最终得到了陶瓷基复合材料的强度,可以大大提高预测结果的精度。
Description
技术领域
本发明属于复合材料力学性能预测领域,具体涉及一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法。
背景技术
纤维增强陶瓷基复合材料由于具有低密度和高耐温性,在航空发动机热端部件的应用上已经显示出巨大的发展潜力。强度是陶瓷基复合材料最重要的性能之一,进行陶瓷基复合材料结构设计要求设计人员能准确预测材料的强度。现有的陶瓷基复合材料强度预测方法往往是根据纤维强度预测陶瓷基复合材料的强度,然而这些方法,如文献[李龙彪,宋迎东,孙志刚.单向纤维增强陶瓷基复合材料单轴拉伸行为[J].复合材料学报,2008,25(04):154-160]中的方法,都认为纤维在某处断裂后,整根纤维都失去承载能力,他们将纤维断裂等效为纤维体积分数的降低,都没有考虑断裂纤维的承载能力。事实上,由于纤维/基体间界面的载荷传递作用,断裂纤维还具有一定的承载的能力。因此,若准确预测陶瓷基复合材料的强度,必须要考虑断裂纤维的承载能力。
发明内容
本发明针对现有技术中的不足,提供一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法。
为实现上述目的,本发明采用以下技术方案:
一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法,其特征在于,包括如下步骤:
步骤1:将陶瓷基复合材料中每根纤维等分为多个纤维微元;
步骤3:给定陶瓷基复合材料的外载σ;
步骤4:以任一基体裂纹平面为基准面,划定影响区,只分析影响区域内的纤维微元;
为优化上述技术方案,采取的具体措施还包括:
进一步地,步骤5中,计算第i根纤维在基体裂纹平面处的应力的方法为:
以该基体裂纹平面为基准面,由近到远在影响区域内依次遍历第i根纤维上的微元,对于第j个微元,选取纤维应力分布模型,计算得到该微元的应力判断该微元是否断裂,若则该微元断裂并停止遍历;若影响区域内所有微元都未断裂,则该纤维在基体裂纹平面处的应力为否则计算该纤维在基体裂纹平面处的应力的方法为式中rf是纤维半径,vf是纤维体积分数,τ是纤维/基体间界面剪应力,Lbreak为断裂微元到基体裂纹平面的距离。
进一步地,步骤5中,选取Weibull模型来确定基体裂纹平面的位置,Weibull模型中基体裂纹等间距分布,相邻两个基体裂纹平面的间距L由下式确定:
式中Lsat为饱和基体裂纹密度,σ0和m为经验参数;根据相邻两个基体裂纹平面的间距L来计算该微元离最近的基体裂纹平面的距离Lflaw。
本发明的有益效果是:本发明提供了一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法,在陶瓷基复合材料中,由于纤维/基体间界面的载荷传递作用,断裂纤维还具有一定的承载能力。本发明考虑了断裂纤维的继续承载,给出了判断纤维断裂位置的方法,计算了断裂纤维在基体裂纹平面处的应力,最终得到了陶瓷基复合材料的强度,可以大大提高预测结果的精度。
附图说明
图1是本发明预测陶瓷基复合材料强度的流程图。
图2是本发明影响区范围示意图。
图3是本发明计算每根纤维在基体裂纹平面处应力的流程图。
图4是本发明遍历纤维上微元的顺序示意图。
图中:1-纤维,2-基体,3-影响区,4-基体裂纹,5-纤维微元。
具体实施方式
现在结合附图对本发明作进一步详细的说明。
本发明提供一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法,如图1所示,该方法包括以下步骤:
步骤1:将陶瓷基复合材料中每根纤维等分为很多个纤维微元。
步骤3:给定陶瓷基复合材料的初始外载σ。
步骤5:如图3所示,计算每根纤维在基体裂纹平面处的应力。以基体裂纹平面为基准面,由近到远在影响区域内依次遍历该纤维上的微元,遍历顺序如图4所示。对于第i根纤维上的第j个微元,选取纤维应力分布模型,计算得到该微元的应力判断该微元是否断裂,若则该微元断裂并停止遍历。计算该纤维在基体裂纹平面处的应力,若影响区域内所有微元都未断裂,则该纤维在基体裂纹平面处的应力为否则计算该纤维在基体裂纹平面处的应力的方法为
实施例:
为了使本技术领域的人员更好地理解本申请中的技术方案和有益效果,本申请在实施例中针对SiC/SiC陶瓷基复合材料进行强度预测。将SiC/SiC陶瓷基复合材料中每根纤维等分为很多个纤维微元,每个微元长度为Δl=0.02mm。本实施例中SiC纤维强度分布如下:
表1 SiC纤维各级缺陷参数
表中Sj表示第j级缺陷的强度,P(Aj)表示长度为Lf=25mm的纤维上出现第j级缺陷的概率。由表1中的P(Aj)计算出各级缺陷出现在纤维微元上的概率用计算机产生一个随机数ξ(0<ξ<1),若ξ<Pj,则认为该微元上存在第j类缺陷,微元的强度等于最低级缺陷对应的强度。
选取剪滞模型来计算纤维微元应力,第i根纤维上第j个微元的应力为:
选取Weibull模型来确定基体裂纹平面的位置。Weibull模型中基体裂纹等间距分布,相邻两个基体裂纹平面的间距L由下式确定:
式中Lsat为饱和基体裂纹密度,σ0和m为经验参数,计算得当前外载下的基体裂纹间距为L=0.41mm。根据相邻两个基体裂纹平面的间距L即可求得该微元离最近的基体裂纹平面的距离Lflaw
根据纤维在基体裂纹平面处的应力,计算纤维断裂后SiC/SiC陶瓷基复合材料的平均应力为由于SiC/SiC陶瓷基复合材料未断裂。增大SiC/SiC陶瓷基复合材料的外载并重新计算纤维断裂后SiC/SiC陶瓷基复合材料的平均应力。取断裂前的最大平均应力为其强度,最终计算得本实施例中SiC/SiC陶瓷基复合材料的强度为581.6MPa,试验结果是570.7MPa,误差为1.9%。
本实施例还计算了当不考虑断裂现在承载时的SiC/SiC陶瓷基复合材料的强度,为207.3MPa,与试验结果误差高达63.7%。因此本发明提供的一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法能显著提高预测精度。
需要注意的是,发明中所引用的如“上”、“下”、“左”、“右”、“前”、“后”等的用语,亦仅为便于叙述的明了,而非用以限定本发明可实施的范围,其相对关系的改变或调整,在无实质变更技术内容下,当亦视为本发明可实施的范畴。
以上仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,应视为本发明的保护范围。
Claims (4)
1.一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法,其特征在于,包括如下步骤:
步骤1:将陶瓷基复合材料中每根纤维等分为多个纤维微元;
步骤3:给定陶瓷基复合材料的外载σ;
步骤4:以任一基体裂纹平面为基准面,划定影响区,只分析影响区域内的纤维微元;
步骤5中,计算第i根纤维在基体裂纹平面处的应力的方法为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910618452.9A CN110309619B (zh) | 2019-07-09 | 2019-07-09 | 一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910618452.9A CN110309619B (zh) | 2019-07-09 | 2019-07-09 | 一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110309619A CN110309619A (zh) | 2019-10-08 |
CN110309619B true CN110309619B (zh) | 2020-05-26 |
Family
ID=68079999
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910618452.9A Active CN110309619B (zh) | 2019-07-09 | 2019-07-09 | 一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110309619B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111474054B (zh) * | 2020-04-21 | 2021-06-29 | 南京航空航天大学 | 一种无机非金属纤维强度分布测算方法 |
CN111751260B (zh) * | 2020-06-30 | 2021-12-21 | 南京航空航天大学 | 陶瓷基纤维束复合材料横截面积和孔隙率测量装置及方法 |
CN113139300B (zh) * | 2021-05-14 | 2023-09-12 | 南京工程学院 | 陶瓷微孔致裂强度定量预测方法、装置、计算机设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105631148A (zh) * | 2015-12-31 | 2016-06-01 | 南京航空航天大学 | 应力氧化环境下单向陶瓷基复合材料力学性能分析方法 |
CN109583052A (zh) * | 2018-11-15 | 2019-04-05 | 西南交通大学 | 纤维增强陶瓷基复合材料微结构精确建模及模量计算方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9562844B2 (en) * | 2014-06-30 | 2017-02-07 | Baker Hughes Incorporated | Systems and devices for sensing corrosion and deposition for oil and gas applications |
CN105912827B (zh) * | 2016-07-06 | 2018-10-16 | 北京航空航天大学 | 一种预报复合材料纤维拉伸失效的能量判据 |
CN106777562B (zh) * | 2016-11-29 | 2018-05-04 | 南京航空航天大学 | 一种确定陶瓷基复合材料强度的方法 |
CN109781546B (zh) * | 2018-12-29 | 2020-09-22 | 南京航空航天大学 | 一种编织陶瓷基复合材料拉伸强度的预测方法 |
CN109684598B (zh) * | 2018-12-29 | 2020-08-21 | 南京航空航天大学 | 一种通过迟滞耗散能预测编织陶瓷基复合材料高温疲劳寿命的方法 |
CN109711100B (zh) * | 2019-01-23 | 2020-12-08 | 南京航空航天大学 | 一种编织陶瓷基复合材料疲劳寿命的预测方法 |
-
2019
- 2019-07-09 CN CN201910618452.9A patent/CN110309619B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105631148A (zh) * | 2015-12-31 | 2016-06-01 | 南京航空航天大学 | 应力氧化环境下单向陶瓷基复合材料力学性能分析方法 |
CN109583052A (zh) * | 2018-11-15 | 2019-04-05 | 西南交通大学 | 纤维增强陶瓷基复合材料微结构精确建模及模量计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110309619A (zh) | 2019-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110309619B (zh) | 一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法 | |
Ferdous et al. | Testing and modelling the fatigue behaviour of GFRP composites–Effect of stress level, stress concentration and frequency | |
Jollivet et al. | Damage of composite materials | |
Wang et al. | A stochastic model for the growth of matrix cracks in composite laminates | |
Karakuzu et al. | Impact characterization of glass/epoxy composite plates: An experimental and numerical study | |
Zweben | Tensile failure of fiber composites. | |
Koricho et al. | Innovative tailored fiber placement technique for enhanced damage resistance in notched composite laminate | |
CN106599507B (zh) | 一种改进b-k准则用于含纤维桥接影响复合材料多向层板分层预测的方法 | |
Okabe et al. | Numerical modeling of progressive damage in fiber reinforced plastic cross-ply laminates | |
Hocheng et al. | Prediction of critical thrust force for tubular composite in drilling-induced delamination by numerical and experimental analysis | |
Carraro et al. | Fatigue damage and stiffness evolution in composite laminates: a damage-based framework | |
Watson et al. | A critical flaw size approach for predicting the strength of bolted glass connections | |
Torabi et al. | Extension of the virtual isotropic material concept to mixed mode I/II loading for predicting the last-ply-failure of U-notched glass/epoxy laminated composite specimens | |
Torabi et al. | On the ability of the notch fracture mechanics in predicting the last-ply-failure of blunt V-notched laminated composite specimens: A hard problem can be easily solved by conventional methods | |
Cao et al. | Analysis of ductile–brittle competitive fracture criteria for tension process of 7050 aluminum alloy based on elastic strain energy density | |
Feng et al. | Investigation on the buckling and postbuckling performance of aero stiffened composite panels under axial compression | |
Jollivet et al. | Fractography, a powerful tool for identifying and understanding fatigue in composite materials | |
Nairn | Applications of finite fracture mechanics for predicting fracture events in composites | |
Smith et al. | A Comparison of Transverse Cracking Phenomena in (0/90) s and (90/0) s CFRP Laminates | |
Li et al. | FE parametric study on the longitudinal tensile strength and damage mechanism of Z‐pinned laminates | |
Guo et al. | Fatigue life prediction considering strength contribution of fibre layers with different orientations for CFRP laminates at high temperature | |
Tang et al. | Predictions of thermo-mechanical behavior of a Nicalon-CAS 0°-90° ceramic matrix composite from constituent materials properties | |
Asp et al. | Effects of pores and voids on the interlaminar delamination toughness of a carbon/epoxy composite | |
Wang et al. | Flexural tests and analysis of notched specimens of glass fiber reinforced composite | |
Abbasi et al. | Shear strengthening of RC beams with FRP composites: database of FE simulations and analysis of studied parameters |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |