CN104020254B - 一种测定复合材料剩余强度与剩余寿命的应变控制方法 - Google Patents

一种测定复合材料剩余强度与剩余寿命的应变控制方法 Download PDF

Info

Publication number
CN104020254B
CN104020254B CN201410213306.5A CN201410213306A CN104020254B CN 104020254 B CN104020254 B CN 104020254B CN 201410213306 A CN201410213306 A CN 201410213306A CN 104020254 B CN104020254 B CN 104020254B
Authority
CN
China
Prior art keywords
formula
strain
fatigue
value
residual
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.)
Expired - Fee Related
Application number
CN201410213306.5A
Other languages
English (en)
Other versions
CN104020254A (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.)
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 CN201410213306.5A priority Critical patent/CN104020254B/zh
Publication of CN104020254A publication Critical patent/CN104020254A/zh
Application granted granted Critical
Publication of CN104020254B publication Critical patent/CN104020254B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

一种测定复合材料剩余强度与剩余寿命的应变控制方法,该方法有三大步骤:步骤一、应变控制剩余强度模型;步骤二、应变控制剩余强度的随机模型;步骤三、模型参数估计。本发明简单实用、操作方便、计算精度高。本发明在测试技术领域里具有较好的实用价值和广阔地应用前景。

Description

一种测定复合材料剩余强度与剩余寿命的应变控制方法
技术领域
本发明提供一种测定复合材料剩余强度与剩余寿命的应变控制方法,属于试验测试技术领域。
背景技术
复合材料剩余强度与剩余寿命测定方法是其疲劳寿命评估的重要前提,由于复合材料疲劳损伤的复杂性,难以采用单一方式定义复合材料损伤,因此,人们先后提出了各种基于刚度降、裂纹密度、裂纹长度等概念的疲劳损伤模型;然而,这些模型难以通过试验方法方便地测定,为此,本发明提出一种测定复合材料剩余强度与剩余寿命的应变控制方法,该方法简单实用、操作方便、计算精度高,能充分而合理地表征复合材料疲劳损伤物理特性与唯象的试验数据规律,具有重要的学术意义和工程应用价值。
发明内容
1、目的:本发明的目的是提供了一种测定复合材料剩余强度与剩余寿命的应变控制方法,该方法具有简单实用、操作方便、计算精度高,并能合理表征其损伤规律等优点。
2、技术方案:本发明一种测定复合材料剩余强度与剩余寿命的应变控制方法,该方法具体步骤如下:
步骤一、应变控制剩余强度模型
疲劳损伤导致强度下降,随时间变化的复合材料有效模量降可表示为
dR ( n ) dn = - f ( r , s , ω ) R b - 1 ( n ) - - - ( 1 )
式中,f(r,s,ω)为最大疲劳应力s、加载频率ω和应力比r的函数。在不考虑加载顺序效应及不改变应力水平的情况下,对上式积分,得到
n=f(r,s,ω)[R0-R(n)]b(2)
式中,R0为拟合强度极限。对于给定的加载频率ω和应力比r,f(r,s,ω)=f(s),则式(2)为
n=f(s)[R0-R(n)]b(3)
式(3)即为剩余强度R-疲劳应力s-疲劳应力循环次数n的关系曲面。根据S-N曲线规律,S-N曲线常采用幂函数式表示:
N=C(S-S0)m(4)
式中,C和m为材料常数,S为疲劳强度,S0为拟合疲劳极限。由式(4)可得
f(s)=C(s-S0)m(5)
将式(5)代入式(3),可获得应力控制剩余强度的方程
n=C(s-S0)m[R0-R(n)]b(6)
静强度R0和剩余强度R(n)可分别由下式求得:
R0=E0εf(7)
R(n)=E(n)εf(8)
式中,εf为断裂应变,E0为初始模量,E(n)为剩余模量。
将式(7)和式(8)代入式(6),得到
n=C0(s-S0)m[E0-E(n)]b(9)
式中,在指定疲劳应力s的条件下,剩余模量R(n)与疲劳应变ε(n)之间存在如下关系:
E ( n ) = s ϵ ( n ) - - - ( 10 )
将式(10)代入式(9),可获得应变控制疲劳剩余强度模型
n = C 0 ( s - S 0 ) m [ E 0 - s ϵ ( n ) ] b - - - ( 11 )
步骤二、应变控制剩余强度的随机模型
将式(11)随机化,即得到应变控制剩余强度的随机模型
n p = C 0 ( s - S 0 ) m [ E 0 - s ϵ ( n ) ] b · exp [ u p kσ ^ ] - - - ( 12 )
n pγ = C 0 ( s - S 0 ) m [ E 0 - s ϵ ( n ) ] b · exp { σ · [ k ^ u p + t γ 1 n + u p 2 ( k ^ 2 - 1 ) ] } - - - ( 13 )
对式(11)随机化,并取对数,得到
Y=a0+a1x1+a2x2+U(14)
式中,Y=lgn,a0=lgC0,a1=m,a2=b,x1=lg(s-S0),U=lgX(n),且U为正态随机变量N[0,σ2]。由式(14)可知,Y为正态随机变量N[a0+a1x1+a2x22],则根据极大似然法,得到
a 0 = y ‾ - a 1 x ‾ 1 + a 2 x ‾ 2 - - - ( 15 )
a 1 = L 12 L 20 - L 22 L 10 L 12 L 21 - L 11 L 22 - - - ( 16 )
a 2 = L 21 L 10 - L 11 L 20 L 12 L 21 - L 11 L 22 - - - ( 17 )
σ = Σ i = 1 l ( y i - a 0 - a 1 x 1 i + a 2 x 2 i ) 2 l - - - ( 18 )
式中
y ‾ = 1 l Σ i = 1 l y i - - - ( 19 )
x ‾ = 1 l Σ i = 1 l x 1 i - - - ( 20 )
x ‾ 2 = 1 l Σ i = 1 l x 2 i - - - ( 21 )
L 11 = Σ i = 1 l ( x 1 i - x ‾ 1 ) 2 - - - ( 22 )
L 22 = Σ i = 1 l ( x 2 i - x ‾ 2 ) 2 - - - ( 23 )
L 12 = Σ i = 1 l ( x 1 i - x ‾ 1 ) ( x 2 i - x ‾ 2 ) - - - ( 24 )
L21=L12(25)
L 10 = Σ i = 1 l ( x 1 i - x ‾ 1 ) ( y i - y ‾ ) - - - ( 26 )
L 20 = Σ i = 1 l ( x 2 i - x ‾ 2 ) ( y i - y ‾ ) - - - ( 27 )
步骤三、模型参数估计
式(15)至式(17)是待定常数E0和S0的二元函数,因此,需要先求出的E0和S0值,再由式(15)至式(18)获得a0、a1、a2和σ。具体的求解步骤如下:
(1)首先,令残差平方和函数
Q ( E 0 , S 0 ) = Σ i = 1 l ( y i - a 0 - a 1 x 1 i - a 2 x 2 i ) 2 - - - ( 28 )
(2)确定E0和S0的取值范围
E0∈(Emax,Emax+△]
S0∈[0,S0min)
式中,Emax=max{E1,E2,…,El},其中Ei(i=1,2,…,l)为剩余模量试验数据;△为一有限值;S0min=min{s1,s2,…,sl},其中si(i=1,2,…,l)为试验疲劳应力取值。
(3)给定一组E0和S0的初始值并分别给定E0和S0的取值步长△1和△2,按式(28)计算Q(E0,S0)的值,寻找Q(E0,S0)的最小值点对应的E0和S0值。
(4)再由上面求解的E0和S0值,按式(15)至式(18)得到a0、a1、a2和σ,最终获得
C 0 = 1 0 y ‾ - a 1 x ‾ 1 - a 2 x ‾ 2 ϵ f b - - - ( 29 )
m = L 12 L 20 - L 22 L 10 L 12 L 21 - L 11 L 22 - - - ( 30 )
b = L 21 L 10 - L 11 L 20 L 12 L 21 - L 11 L 22 - - - ( 31 )
将式(29)至式(31)代入式(12)和式(13)即可。
3、优点及功效:本发明一种测定复合材料剩余强度与剩余寿命的应力控制方法,其特点是简单实用、操作方便、计算精度高。
附图说明
图1为是本发明所述方法的流程框图。
图中符号说明如下:Q为残差平方和函数,E0、S0、C0、m和b均为待定常数。
具体实施方式
图1为本发明所述方法的流程框图,本发明分三大步骤实现,具体为:
步骤一、应变控制剩余强度模型
疲劳损伤导致强度下降,随时间变化的复合材料有效模量降可表示为
dR ( n ) dn = - f ( r , s , ω ) R b - 1 ( n ) - - - ( 32 )
式中,f(r,s,ω)为最大疲劳应力s、加载频率ω和应力比r的函数。在不考虑加载顺序效应及不改变应力水平的情况下,对上式积分,得到
n=f(r,s,ω)[R0-R(n)]b(33)
式中,R0为拟合强度极限。对于给定的加载频率ω和应力比r,f(r,s,ω)=f(s),则式(33)为
n=f(s)[R0-R(n)]b(34)
式(34)即为剩余强度R-疲劳应力s-疲劳应力循环次数n的关系曲面。根据S-N曲线规律,S-N曲线常采用幂函数式表示:
N=C(S-S0)m(35)
式中,C和m为材料常数,S为疲劳强度,S0为拟合疲劳极限。由式(35)可得
f(s)=C(s-S0)m(36)
将式(36)代入式(34),可获得应力控制剩余强度的方程
n=C(s-S0)m[R0-R(n)]b(37)
静强度R0和剩余强度R(n)可分别由下式求得:
R0=E0εf(38)
R(n)=E(n)εf(39)
式中,εf为断裂应变,E0为初始模量,E(n)为剩余模量。
将式(38)和式(39)代入式(37),得到
n=C0(s-S0)m[E0-E(n)]b(40)
式中,在指定疲劳应力s的条件下,剩余模量R(n)与疲劳应变ε(n)之间存在如下关系:
E ( n ) = s ϵ ( n ) - - - ( 41 )
将式(41)代入式(40),可获得应变控制疲劳剩余强度模型
n = C 0 ( s - S 0 ) m [ E 0 - s ϵ ( n ) ] b - - - ( 42 )
步骤二、应变控制剩余强度的随机模型
将式(42)随机化,即得到应变控制剩余强度的随机模型
n p = C 0 ( s - S 0 ) m [ E 0 - s ϵ ( n ) ] b · exp [ u p kσ ^ ] - - - ( 43 )
n pγ = C 0 ( s - S 0 ) m [ E 0 - s ϵ ( n ) ] b · exp { σ · [ k ^ u p + t γ 1 n + u p 2 ( k ^ 2 - 1 ) ] } - - - ( 13 )
对式(42)随机化,并取对数,得到
Y=a0+a1x1+a2x2+U(45)
式中,Y=lgn,a0=lgC0,a1=m,a2=b,x1=lg(s-S0),U=lgX(n),且U为正态随机变量N[0,σ2]。由式(45)可知,Y为正态随机变量N[a0+a1x1+a2x22],则根据极大似然法,得到
a 0 = y ‾ - a 1 x ‾ 1 + a 2 x ‾ 2 - - - ( 46 )
a 1 = L 12 L 20 - L 22 L 10 L 12 L 21 - L 11 L 22 - - - ( 47 )
a 2 = L 21 L 10 - L 11 L 20 L 12 L 21 - L 11 L 22 - - - ( 48 )
σ = Σ i = 1 l ( y i - a 0 - a 1 x 1 i + a 2 x 2 i ) 2 l - - - ( 49 )
式中
y ‾ = 1 l Σ i = 1 l y i - - - ( 50 )
x ‾ = 1 l Σ i = 1 l x 1 i - - - ( 51 )
x ‾ 2 = 1 l Σ i = 1 l x 2 i - - - ( 52 )
L 11 = Σ i = 1 l ( x 1 i - x ‾ 1 ) 2 - - - ( 53 )
L 22 = Σ i = 1 l ( x 2 i - x ‾ 2 ) 2 - - - ( 54 )
L 12 = Σ i = 1 l ( x 1 i - x ‾ 1 ) ( x 2 i - x ‾ 2 ) - - - ( 55 )
L21=L12(56)
L 10 = Σ i = 1 l ( x 1 i - x ‾ 1 ) ( y i - y ‾ ) - - - ( 57 )
L 20 = Σ i = 1 l ( x 2 i - x ‾ 2 ) ( y i - y ‾ ) - - - ( 58 )
步骤三、模型参数估计
式(46)至式(48)是待定常数E0和S0的二元函数,因此,需要先求出的E0和S0值,再由式(46)至式(49)获得a0、a1、a2和σ。具体的求解步骤如下:
(1)首先,令残差平方和函数
Q ( E 0 , S 0 ) = Σ i = 1 l ( y i - a 0 - a 1 x 1 i - a 2 x 2 i ) 2 - - - ( 59 )
(2)确定E0和S0的取值范围
E0∈(Emax,Emax+△]
S0∈[0,S0min)
式中,Emax=max{E1,E2,…,El},其中Ei(i=1,2,…,l)为剩余模量试验数据;△为一有限值;S0min=min{s1,s2,…,sl},其中si(i=1,2,…,l)为试验疲劳应力取值。
(3)给定一组E0和S0的初始值并分别给定E0和S0的取值步长△1和△2,按式(59)计算Q(E0,S0)的值,寻找Q(E0,S0)的最小值点对应的E0和S0值。
(4)再由上面求解的E0和S0值,按式(46)至式(49)得到a0、a1、a2和σ,最终获得
C 0 = 1 0 y ‾ - a 1 x ‾ 1 - a 2 x ‾ 2 ϵ f b - - - ( 60 )
m = L 12 L 20 - L 22 L 10 L 12 L 21 - L 11 L 22 - - - ( 61 )
b = L 21 L 10 - L 11 L 20 L 12 L 21 - L 11 L 22 - - - ( 62 )
将式(60)至式(62)代入式(43)和式(44)即可。

Claims (1)

1.一种测定复合材料剩余强度与剩余寿命的应变控制方法,其特征在于:该方法具体步骤如下:
步骤一、应变控制剩余强度模型
疲劳损伤导致强度下降,随时间变化的复合材料有效模量降表示为
式中,f(r,s,ω)为最大疲劳应力s、加载频率ω和应力比r的函数;在不考虑加载顺序效应及不改变应力水平的情况下,对上式积分,得到
n=f(r,s,ω)[R0-R(n)]b(2)
式中,R0为拟合强度极限,对于给定的加载频率ω和应力比r,f(r,s,ω)=f(s),则式(2)为
n=f(s)[R0-R(n)]b(3)
式(3)即为剩余强度R-疲劳应力s-疲劳应力循环次数n的关系曲面;根据S-N曲线规律,S-N曲线常采用幂函数式表示:
N=C(S-S0)m(4)
式中,C和m为材料常数,S为疲劳强度,S0为拟合疲劳极限;由式(4)得
f(s)=C(s-S0)m(5)
将式(5)代入式(3),获得应力控制剩余强度的方程
n=C(s-S0)m[R0-R(n)]b(6)
R0和剩余强度R(n)分别由下式求得:
R0=E0εf(7)
R(n)=E(n)εf(8)
式中,εf为断裂应变,E0为初始模量,E(n)为剩余模量;
将式(7)和式(8)代入式(6),得到
n=C0(s-S0)m[E0-E(n)]b(9)
式中,在指定疲劳应力s的条件下,剩余模量E(n)与疲劳应变ε(n)之间存在如下关系:
将式(10)代入式(9),获得应变控制疲劳剩余强度模型
步骤二、应变控制剩余强度的随机模型
将式(11)随机化,即得到应变控制剩余强度的随机模型
对式(11)随机化,并取对数,得到
Y=a0+a1x1+a2x2+U(14)
式中,Y=lgn,a0=lgC0,a1=m,a2=b,x1=lg(s-S0),U=lgX(n),且U为正态随机变量N[0,σ2];由式(14)可知,Y为正态随机变量N[a0+a1x1+a2x22],则根据极大似然法,得到
式中
L21=L12(25)
步骤三、模型参数估计
式(15)至式(17)是待定常数E0和S0的二元函数,因此,需要先求出E0和S0值,再由式(15)至式(18)获得a0、a1、a2和σ;具体的求解步骤如下:
(1)首先,令残差平方和函数
(2)确定E0和S0的取值范围
E0∈(Emax,Emax+Δ]
S0∈[0,S0min)
式中,Emax=max{E1,E2,…,El},其中Ei(i=1,2,…,l)为剩余模量试验数据;Δ为一有限值;S0min=min{s1,s2,…,sl},其中si(i=1,2,…,l)为试验疲劳应力取值;
(3)给定一组E0和S0的初始值并分别给定E0和S0的取值步长Δ1和Δ2,按式(28)计算Q(E0,S0)的值,寻找Q(E0,S0)的最小值点对应的E0和S0值;
(4)再由上面求解的E0和S0值,按式(15)至式(18)得到a0、a1、a2和σ,最终获得
将式(29)至式(31)代入式(12)和式(13)即可。
CN201410213306.5A 2014-05-20 2014-05-20 一种测定复合材料剩余强度与剩余寿命的应变控制方法 Expired - Fee Related CN104020254B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410213306.5A CN104020254B (zh) 2014-05-20 2014-05-20 一种测定复合材料剩余强度与剩余寿命的应变控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410213306.5A CN104020254B (zh) 2014-05-20 2014-05-20 一种测定复合材料剩余强度与剩余寿命的应变控制方法

Publications (2)

Publication Number Publication Date
CN104020254A CN104020254A (zh) 2014-09-03
CN104020254B true CN104020254B (zh) 2016-01-13

Family

ID=51437125

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410213306.5A Expired - Fee Related CN104020254B (zh) 2014-05-20 2014-05-20 一种测定复合材料剩余强度与剩余寿命的应变控制方法

Country Status (1)

Country Link
CN (1) CN104020254B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105044146B (zh) * 2015-07-02 2017-12-22 西安电子科技大学 一种两相复合材料的随机热均化分析方法
CN105956236B (zh) * 2016-04-22 2019-03-12 西安交通大学 一种双重更新的四因素随机退化模型齿轮寿命预测方法
CN106202906B (zh) * 2016-07-06 2018-10-09 北京航空航天大学 一种腐蚀疲劳性能表征和寿命估算方法
CN106383962B (zh) * 2016-09-30 2019-06-28 北京航空航天大学 一种热压成型平面编织复合材料的残余热应力估算方法
CN108204925B (zh) * 2016-12-16 2020-03-20 海口未来技术研究院 复合材料的疲劳寿命预测方法及预测系统
CN108681641B (zh) * 2018-05-18 2021-01-08 北京航空航天大学 一种考虑缺陷尺寸效应的复合材料缺陷容限性能表征方法
CN109446594A (zh) * 2018-10-10 2019-03-08 北京航空航天大学 一种测定复合材料疲劳s-n曲线的剩余强度控制方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6778914B1 (en) * 2000-03-17 2004-08-17 University Of Delaware Dynamic interphase-loading apparatus and method of using the same
CN101231222A (zh) * 2008-02-15 2008-07-30 上海理工大学 一种剩余强度和剩余寿命的快速、无损预测方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6778914B1 (en) * 2000-03-17 2004-08-17 University Of Delaware Dynamic interphase-loading apparatus and method of using the same
CN101231222A (zh) * 2008-02-15 2008-07-30 上海理工大学 一种剩余强度和剩余寿命的快速、无损预测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
A strain-based residual strength model of carbon fibre/epoxy composites based on CAI and fatigue residual strength concepts;Jun-Jiang Xiong etal;《Composite Structures》;20071017;第85卷(第1期);第30页"术语表",第37-38页 4. Residual fatigue strength surface model *

Also Published As

Publication number Publication date
CN104020254A (zh) 2014-09-03

Similar Documents

Publication Publication Date Title
CN104020254B (zh) 一种测定复合材料剩余强度与剩余寿命的应变控制方法
CN104035455B (zh) 一种测定复合材料剩余强度与剩余寿命的应力控制方法
CN105808865A (zh) 一种低温疲劳性能表征与寿命估算的方法
CN105466772B (zh) 一种基于临界面法的多轴短裂纹扩展寿命预测方法
CN105067457B (zh) 一种腐蚀裂纹扩展性能表征与寿命估算的方法
Yin et al. Fractional description of mechanical property evolution of soft soils during creep
Zhang et al. Identifying accurate crack initiation and propagation thresholds in siliceous siltstone and limestone
CN104809311B (zh) 一种基于多因素融合修正的结构件剩余寿命预测方法
CN102129512A (zh) 基于Paris公式的疲劳寿命分析方法
CN107436983A (zh) 一种基于多元样本差异的o型橡胶密封圈寿命预测方法
Zolotukhin et al. Machine learning in reservoir permeability prediction and modelling of fluid flow in porous media
CN104977226B (zh) 岩石密度测量方法与岩石密度测量装置
Hossain et al. MHD free convection and mass transfer flow through a vertical oscillatory porous plate with Hall, ion-slip currents and heat source in a rotating system
CN104655505A (zh) 一种基于仪器化球压入技术的残余应力检测方法
CN111680440A (zh) 不同温度下数字岩心声学特性分析方法
Zhang et al. Influences of test method and loading history on permeability of tight reservoir rocks
CN102928171B (zh) 航天器总漏率测试结果的不确定度确定方法
Pourahmadian et al. On the elastic-wave imaging and characterization of fractures with specific stiffness
CN105550457A (zh) 动态演化模型校正方法及系统
CN103528844B (zh) 基于经验模态分解的结构损伤预警方法
CN106202906B (zh) 一种腐蚀疲劳性能表征和寿命估算方法
CN104568640A (zh) 一种水泥浆稳定性评价方法
Konovalov et al. Numerical estimation of effective elastic properties of elastomer composites under finite strains using spectral element method with CAE Fidesys
CN106885891B (zh) 由直剪和单轴压缩试验确定砂土邓肯‑张模型参数的方法
Wang et al. Reference strain γr in hyperbolic modeling of dynamic shear modulus of soils

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160113

Termination date: 20200520