第六章-连续损伤力学.pptVIP

2021/2/7 * 式中,g和k是材料参数。积分此式,利用初始条件(N=0时,ω=0)和破坏条件( N= Nf时,ω=1),很容易导出损伤度ω随循环周次N变化的关系: 将(6)式代入(5)中得: 积分此式,在连续载荷情况下,并考虑初始条件和破坏条件（ N= N*时,ω=1 ）,则有: 即: 可见疲劳损伤演变方程式（5）也是与线性累积损伤理论等价的。 （6） （7） 2021/2/7 * 但不少实验证明,许多材料的参数k与循环应力幅σa相关,即k = k (σa) 。相应地,式（6）应写为: 可见,当参数k与σa无关,损伤演变方程式（8）符合线性累积损伤理论(图a);当参数k与σa相关(图b) ,则容易由图推出: （8） 说明属于非线性累积。 2021/2/7 * 对两级循环加载的疲劳问题,即k=1,2,且σa1 ＞σa2 ,其ω~N/Nf关系如下图所示。若先加σa1 ,循环周次ΔN1;再加σa2 ,循环周次ΔN2时断裂。若已知σa1 和σa2下的疲劳寿命分别为N1f和N2f,则可得到: 若改变加载次序:先加σa2 ,相应周次ΔN2 ;再加σa1 ,相应周次ΔN1 ,而断裂得到的累积结果为: 说明若参数k与σa相关,则累积的结果与加载次序有关。 2021/2/7 * 2021/2/7 * 2）高周疲劳损伤模型 （1）基于微塑性分析的模型 对于高周疲劳,通常循环应力的最大值低于屈服极限（但高于疲劳极限）,材料仍会发生疲劳破坏。这是由于材料某些局部的细观组织发生了塑性变形,即所谓微塑性的不可逆变形。循环应力造成的循环微塑性应变是高周疲劳的主要微观机制。 假定细观结构的应力-应变关系与宏观的关系一致,设晶粒屈服后,其微塑性应力-应变关系仍用幂律表示如下: 式中是 微塑性应变;ka和na是材料参数。 （9） 2021/2/7 * Lemaitre认为微塑性损伤演变方程可以有类似于宏塑性损伤演变方程的形式,参照弹塑性损伤公式并用累积微塑性应变率 代替累积宏塑性应变率 ,有: 这里δ0和α0是材料参数。在多维情况下,可设: 这里k、?也是材料参数。将此式代入（10）,得到微塑性损伤演变方程的具体形式: 式中,Tr是三轴性因子; , （10） 或 （11） （12） 2021/2/7 * 在单轴应力情况下, , , 以及Tr =1,则由式（12）可得单轴应力下微塑性损伤演变方程: 若循环应力是比例加载,并且 并考虑到一次加载中,损伤变化很小而加以忽略。对式（12）积分得: 从而得到疲劳损伤演化方程: （13） （14） 或 2021/2/7 * 对于稳定的高周疲劳问题,积分上式,并利用初始条件(N=0时,ψ=1)和破坏条件( N=Nf时,ψ =0),可得: 式中 在1维情况下,对拉伸脉动循环, , ,疲劳损伤演变方程简化为: 积分此式,利用初始条件与破坏条件,可得式（15）,但疲劳寿命表达式应是: （15） （16） （17） 或 2021/2/7 * 代入式（17）,单轴疲劳损伤演变方程改为: 对连续变化的载荷,积分上式,取破坏条件N=N*,ω=1时,有: 得: 疲劳损伤演变方程与线性累积疲劳损伤理论是一致的。 （18） （19） 2021/2/7 * (2) Choboche模型 Chaboche研究疲劳损伤演变时,考虑到平均应力 的影响,疲劳损伤的演变方程表示为: 式中,f是损伤度的函数;β、b是材料参数,σ0和σ1分别为脉动循环应力和对称循环应力下的疲劳极限,σa是材料的拉伸强度极限。且b与平均应力有关: （a）对 ,这里β1也是材料参数,于是疲劳损伤演变方程写为 （20） （21） 2021/2/7 * 积分上式,利用初始条件,得到: （b