CN110309619A - 一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法 - Google Patents

一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法 Download PDF

Info

Publication number
CN110309619A
CN110309619A CN201910618452.9A CN201910618452A CN110309619A CN 110309619 A CN110309619 A CN 110309619A CN 201910618452 A CN201910618452 A CN 201910618452A CN 110309619 A CN110309619 A CN 110309619A
Authority
CN
China
Prior art keywords
fiber
ceramic matric
matric composite
stress
infinitesimal
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.)
Granted
Application number
CN201910618452.9A
Other languages
English (en)
Other versions
CN110309619B (zh
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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201910618452.9A priority Critical patent/CN110309619B/zh
Publication of CN110309619A publication Critical patent/CN110309619A/zh
Application granted granted Critical
Publication of CN110309619B publication Critical patent/CN110309619B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/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：将陶瓷基复合材料中每根纤维等分为多个纤维微元；
步骤2：根据纤维强度分布和纤维强度模型，确定各纤维微元的强度其中i、j表示第i根纤维上的第j个微元；
步骤3：给定陶瓷基复合材料的外载σ；
步骤4：以任一基体裂纹平面为基准面，划定影响区，只分析影响区域内的纤维微元；
步骤5：计算每根纤维在该基体裂纹平面处的应力
步骤6：根据纤维在该基体裂纹平面处的应力计算纤维断裂后的陶瓷基复合材料平均应力
步骤7：若则陶瓷基复合材料未断裂，增大陶瓷基复合材料的外载为σ+Δσ并返回步骤4；否则陶瓷基复合材料断裂，取陶瓷基复合材料断裂前的最大平均应力为其强度S，即
为优化上述技术方案，采取的具体措施还包括：
进一步地，步骤4中，影响区长度为式中rf是纤维半径，vf是纤维体积分数，τ是纤维/基体间界面剪应力。
进一步地，步骤5中，计算第i根纤维在基体裂纹平面处的应力的方法为：
以该基体裂纹平面为基准面，由近到远在影响区域内依次遍历第i根纤维上的微元，对于第j个微元，选取纤维应力分布模型，计算得到该微元的应力判断该微元是否断裂，若则该微元断裂并停止遍历；若影响区域内所有微元都未断裂，则该纤维在基体裂纹平面处的应力为否则计算该纤维在基体裂纹平面处的应力的方法为式中rf是纤维半径，vf是纤维体积分数，τ是纤维/基体间界面剪应力，Lbreak为断裂微元到基体裂纹平面的距离。
进一步地，步骤5中，选取剪滞模型来计算纤维微元应力第i根纤维上第j个微元的应力为：
式中Lflaw为该微元离最近的基体裂纹平面的距离，d是脱粘区长度，是粘结区纤维应力。
进一步地，步骤5中，选取Weibull模型来确定基体裂纹平面的位置，Weibull模型中基体裂纹等间距分布，相邻两个基体裂纹平面的间距L由下式确定：
式中Lsat为饱和基体裂纹密度，σ0和m为经验参数；根据相邻两个基体裂纹平面的间距L来计算该微元离最近的基体裂纹平面的距离Lflaw。
进一步地，步骤6中，计算纤维断裂后陶瓷基复合材料平均应力的方法为：其中Nf是纤维总数，vf是纤维体积分数。
本发明的有益效果是：本发明提供了一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法，在陶瓷基复合材料中，由于纤维/基体间界面的载荷传递作用，断裂纤维还具有一定的承载能力。本发明考虑了断裂纤维的继续承载，给出了判断纤维断裂位置的方法，计算了断裂纤维在基体裂纹平面处的应力，最终得到了陶瓷基复合材料的强度，可以大大提高预测结果的精度。
附图说明
图1是本发明预测陶瓷基复合材料强度的流程图。
图2是本发明影响区范围示意图。
图3是本发明计算每根纤维在基体裂纹平面处应力的流程图。
图4是本发明遍历纤维上微元的顺序示意图。
图中：1-纤维，2-基体，3-影响区，4-基体裂纹，5-纤维微元。
具体实施方式
现在结合附图对本发明作进一步详细的说明。
本发明提供一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法，如图1所示，该方法包括以下步骤：
步骤1：将陶瓷基复合材料中每根纤维等分为很多个纤维微元。
步骤2：根据纤维强度分布确定各纤维微元的强度
步骤3：给定陶瓷基复合材料的初始外载σ。
步骤4：如图2所示，以任一基体裂纹平面为基准面，划定影响区，影响区长度为
步骤5：如图3所示，计算每根纤维在基体裂纹平面处的应力。以基体裂纹平面为基准面，由近到远在影响区域内依次遍历该纤维上的微元，遍历顺序如图4所示。对于第i根纤维上的第j个微元，选取纤维应力分布模型，计算得到该微元的应力判断该微元是否断裂，若则该微元断裂并停止遍历。计算该纤维在基体裂纹平面处的应力，若影响区域内所有微元都未断裂，则该纤维在基体裂纹平面处的应力为否则计算该纤维在基体裂纹平面处的应力的方法为
步骤6：计算纤维断裂后的陶瓷基复合材料平均应力
步骤7：若则陶瓷基复合材料未断裂，增大陶瓷基复合材料的外载为σ+Δσ并返回步骤4；否则陶瓷基复合材料断裂，取陶瓷基复合材料断裂前的最大平均应力为其强度，即
实施例：
为了使本技术领域的人员更好地理解本申请中的技术方案和有益效果，本申请在实施例中针对SiC/SiC陶瓷基复合材料进行强度预测。将SiC/SiC陶瓷基复合材料中每根纤维等分为很多个纤维微元，每个微元长度为Δl＝0.02mm。本实施例中SiC纤维强度分布如下：
表1 SiC纤维各级缺陷参数
表中Sj表示第j级缺陷的强度，P(Aj)表示长度为Lf＝25mm的纤维上出现第j级缺陷的概率。由表1中的P(Aj)计算出各级缺陷出现在纤维微元上的概率用计算机产生一个随机数ξ(0＜ξ＜1)，若ξ＜Pj，则认为该微元上存在第j类缺陷，微元的强度等于最低级缺陷对应的强度。
本实施例中给定SiC/SiC陶瓷基复合材料的初始外载为σ＝400MPa，以任一基体裂纹平面为基准面，计算得当前外载下影响区长度为以该基体裂纹平面为基准面，由近到远在影响区域内依次遍历该纤维上的微元。
选取剪滞模型来计算纤维微元应力，第i根纤维上第j个微元的应力为：
式中Lflaw为该微元离最近的基体裂纹平面的距离，d是脱粘区长度，是粘结区纤维应力。
选取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 (6)

1.一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法，其特征在于，包括如下步骤：
步骤1：将陶瓷基复合材料中每根纤维等分为多个纤维微元；
步骤2：根据纤维强度分布和纤维强度模型，确定各纤维微元的强度其中i、j表示第i根纤维上的第j个微元；
步骤3：给定陶瓷基复合材料的外载σ；
步骤4：以任一基体裂纹平面为基准面，划定影响区，只分析影响区域内的纤维微元；
步骤5：计算每根纤维在该基体裂纹平面处的应力
步骤6：根据纤维在该基体裂纹平面处的应力计算纤维断裂后的陶瓷基复合材料平均应力
步骤7：若则陶瓷基复合材料未断裂，增大陶瓷基复合材料的外载为σ+Δσ并返回步骤4；否则陶瓷基复合材料断裂，取陶瓷基复合材料断裂前的最大平均应力为其强度S，即
2.如权利要求1所述的一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法，其特征在于：步骤4中，影响区长度为式中rf是纤维半径，vf是纤维体积分数，τ是纤维/基体间界面剪应力。
3.如权利要求1所述的一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法，其特征在于：步骤5中，计算第i根纤维在基体裂纹平面处的应力的方法为：
以该基体裂纹平面为基准面，由近到远在影响区域内依次遍历第i根纤维上的微元，对于第j个微元，选取纤维应力分布模型，计算得到该微元的应力判断该微元是否断裂，若则该微元断裂并停止遍历；若影响区域内所有微元都未断裂，则该纤维在基体裂纹平面处的应力为否则计算该纤维在基体裂纹平面处的应力的方法为式中rf是纤维半径，vf是纤维体积分数，τ是纤维/基体间界面剪应力，Lbreak为断裂微元到基体裂纹平面的距离。
4.如权利要求3所述的一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法，其特征在于：步骤5中，选取剪滞模型来计算纤维微元应力第i根纤维上第j个微元的应力为：
式中Lflaw为该微元离最近的基体裂纹平面的距离，d是脱粘区长度，是粘结区纤维应力。
5.如权利要求4所述的一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法，其特征在于：步骤5中，选取Weibull模型来确定基体裂纹平面的位置，Weibull模型中基体裂纹等间距分布，相邻两个基体裂纹平面的间距L由下式确定：
式中Lsat为饱和基体裂纹密度，σ0和m为经验参数；根据相邻两个基体裂纹平面的间距L来计算该微元离最近的基体裂纹平面的距离Lflaw。
6.如权利要求1所述的一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法，其特征在于：步骤6中，计算纤维断裂后陶瓷基复合材料平均应力的方法为：其中Nf是纤维总数，vf是纤维体积分数。
CN201910618452.9A 2019-07-09 2019-07-09 一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法 Active CN110309619B (zh)

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 true CN110309619A (zh) 2019-10-08
CN110309619B 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)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111474054A (zh) * 2020-04-21 2020-07-31 南京航空航天大学 一种无机非金属纤维强度分布测算方法
CN111751260A (zh) * 2020-06-30 2020-10-09 南京航空航天大学 陶瓷基纤维束复合材料横截面积和孔隙率测量装置及方法
CN113139300A (zh) * 2021-05-14 2021-07-20 南京工程学院 陶瓷微孔致裂强度定量预测方法、装置、计算机设备及存储介质

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150377765A1 (en) * 2014-06-30 2015-12-31 Baker Hughes Incorporated Systems and devices for sensing corrosion and deposition for oil and gas applications
CN105631148A (zh) * 2015-12-31 2016-06-01 南京航空航天大学 应力氧化环境下单向陶瓷基复合材料力学性能分析方法
CN105912827A (zh) * 2016-07-06 2016-08-31 北京航空航天大学 一种预报复合材料纤维拉伸失效的能量判据
CN106777562A (zh) * 2016-11-29 2017-05-31 南京航空航天大学 一种确定陶瓷基复合材料强度的方法
CN109583052A (zh) * 2018-11-15 2019-04-05 西南交通大学 纤维增强陶瓷基复合材料微结构精确建模及模量计算方法
CN109684598A (zh) * 2018-12-29 2019-04-26 南京航空航天大学 一种通过迟滞耗散能预测编织陶瓷基复合材料高温疲劳寿命的方法
CN109711100A (zh) * 2019-01-23 2019-05-03 南京航空航天大学 一种编织陶瓷基复合材料疲劳寿命的预测方法
CN109781546A (zh) * 2018-12-29 2019-05-21 南京航空航天大学 一种编织陶瓷基复合材料拉伸强度的预测方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150377765A1 (en) * 2014-06-30 2015-12-31 Baker Hughes Incorporated Systems and devices for sensing corrosion and deposition for oil and gas applications
CN105631148A (zh) * 2015-12-31 2016-06-01 南京航空航天大学 应力氧化环境下单向陶瓷基复合材料力学性能分析方法
CN105912827A (zh) * 2016-07-06 2016-08-31 北京航空航天大学 一种预报复合材料纤维拉伸失效的能量判据
CN106777562A (zh) * 2016-11-29 2017-05-31 南京航空航天大学 一种确定陶瓷基复合材料强度的方法
CN109583052A (zh) * 2018-11-15 2019-04-05 西南交通大学 纤维增强陶瓷基复合材料微结构精确建模及模量计算方法
CN109684598A (zh) * 2018-12-29 2019-04-26 南京航空航天大学 一种通过迟滞耗散能预测编织陶瓷基复合材料高温疲劳寿命的方法
CN109781546A (zh) * 2018-12-29 2019-05-21 南京航空航天大学 一种编织陶瓷基复合材料拉伸强度的预测方法
CN109711100A (zh) * 2019-01-23 2019-05-03 南京航空航天大学 一种编织陶瓷基复合材料疲劳寿命的预测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李龙彪等: ""单向纤维增强陶瓷基复合材料单轴拉伸行为"", 《复合材料学报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111474054A (zh) * 2020-04-21 2020-07-31 南京航空航天大学 一种无机非金属纤维强度分布测算方法
CN111474054B (zh) * 2020-04-21 2021-06-29 南京航空航天大学 一种无机非金属纤维强度分布测算方法
CN111751260A (zh) * 2020-06-30 2020-10-09 南京航空航天大学 陶瓷基纤维束复合材料横截面积和孔隙率测量装置及方法
CN111751260B (zh) * 2020-06-30 2021-12-21 南京航空航天大学 陶瓷基纤维束复合材料横截面积和孔隙率测量装置及方法
CN113139300A (zh) * 2021-05-14 2021-07-20 南京工程学院 陶瓷微孔致裂强度定量预测方法、装置、计算机设备及存储介质
CN113139300B (zh) * 2021-05-14 2023-09-12 南京工程学院 陶瓷微孔致裂强度定量预测方法、装置、计算机设备及存储介质

Also Published As

Publication number Publication date
CN110309619B (zh) 2020-05-26

Similar Documents

Publication Publication Date Title
Jollivet et al. Damage of composite materials
Wang et al. Fatigue behavior of cracked steel plates strengthened with different CFRP systems and configurations
CN110309619A (zh) 一种考虑断裂纤维承载的陶瓷基复合材料强度预测方法
Quaresimin et al. Damage evolution under cyclic multiaxial stress state: A comparative analysis between glass/epoxy laminates and tubes
Askarinejad et al. Fracture and mixed-mode resistance curve behavior of bamboo
Srinivasa et al. Fracture morphology of carbon fiber reinforced plastic composite laminates
Xie et al. Size effect research of tensile strength of bamboo scrimber based on boundary effect model
Cho et al. Distribution of matrix cracks in a uniaxial ceramic composite
Carraro et al. Fatigue damage and stiffness evolution in composite laminates: a damage-based framework
Yao et al. A validation of a modified Paris relation for fatigue delamination growth in unidirectional composite laminates
Yao et al. Thickness effects on fibre-bridged fatigue delamination growth in composites
Sabaghi et al. Critical strain energy release rate of woven carbon/epoxy composites subjected to thermal cyclic loading
Shao et al. Scaling analysis of the tensile strength of bamboo fibers using Weibull statistics
Larco et al. The effects of fibre volume fraction on a glass-epoxy composite material
Talreja A mechanisms‐based reliability model for fatigue of composite laminates
Qiao et al. Mixed-mode fracture of hybrid material bonded interfaces under four-point bending
Zhuk et al. Compressive behaviour of thin-skin stiffened composite panels with a stress raiser
CN110322937A (zh) 一种考虑制备工艺的陶瓷基复合材料强度预测方法
Abbasi et al. Shear strengthening of RC beams with FRP composites: database of FE simulations and analysis of studied parameters
Kaveh Determining the crack propagation in DCB composite specimens by fiber bridging modeling
Ke et al. Fatigue performance of carbon fiber‐reinforced polymer strengthened single‐edge cracked steel plates: Effects of adhesive property and CFRP arrangement
Kulkarni et al. Effect of processing parameters on Tensile Strength of GFRP with Liquid Silicon Rubber and reinforced with fine Aluminium Powder and Silica Powder
Liu et al. An IIPG-based finite element framework in MOOSE for modeling fiber reinforced composite failure governed by extrinsic cohesive laws
Meirinhos et al. Tapered laminates under static and fatigue tension loading
Jones et al. Fracture toughness prediction of composite materials

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