材料强度的统计计算.pptVIP

计算力学课程报告作者：张鹏第1页，共26页，星期日，2025年，2月5日强度统计计算与多尺度力学计算力学课程报告作者：张鹏第2页，共26页，星期日，2025年，2月5日算法介绍1脆性材料的强度对裂纹分布非常敏感材料强度微裂纹平均密度微裂纹平均长度长度的涨落密度的涨落微裂纹串接过程取决于裂纹间的强相互作用影响原因第3页，共26页，星期日，2025年，2月5日算法介绍2为揭示强度对微裂纹分布的敏感性，需要采用统计计算方法材料破坏概率（与统计变量的关系）裂纹扩展后的期望长度拟和破坏曲线……微裂纹平均密度微裂纹平均长度长度的涨落密度的涨落统计计算第4页，共26页，星期日，2025年，2月5日问题描述1无限大平板，包含N个共线裂纹，无穷远处作用有均匀拉应力σ∞第5页，共26页，星期日，2025年，2月5日问题描述2a：半裂纹长c：裂纹间距a和c都是统计变量，它们的统计分布用函数f(a)、p(c)来表示，c-,c+,a-,a+是c和a的下届和上届。第6页，共26页，星期日，2025年，2月5日问题描述3目标求出σ∞作用下的材料破坏概率用直接数值模拟进行校验用适当的分布函数拟和破坏概率曲线第7页，共26页，星期日，2025年，2月5日简化问题主要考虑相邻的两个微裂纹之间的强相互作用第8页，共26页，星期日，2025年，2月5日裂尖应力强度因子A点的SIF（stressintensifyfactor）：其中F是无量纲函数第9页，共26页，星期日，2025年，2月5日一个特例（a’=a’’=a0）N个长度相同的裂纹，间距不同（即a’=a’’=a0）则f(a)是Diracdelta函数δ(a-a0)当a’=a’’=a0时，c的分布函数p(c)是一个正态分布，c的取值范围为（c-,c+），平均间距记为c0第10页，共26页，星期日，2025年，2月5日不同c0/a0比率下的KA－c/a0图当a’=a’’=a0时，KA是一个c/a0单调递减函数第11页，共26页，星期日，2025年，2月5日σ∞的临界值σth∞σth∞：满足使裂纹扩展的σ∞最小值KIC是基体断裂刚度第12页，共26页，星期日，2025年，2月5日对应给定σ∞的临界裂纹间距ccr1ccr1：σ∞作用下不至于连通的最小裂纹间距第13页，共26页，星期日，2025年，2月5日不同σ∞对应的不同情况如果σ∞σth∞，应力小于使基体断裂的最小应力，裂纹不会扩展。如果σ∞≥σth∞，则间距小于ccr1的相邻裂纹将连通。第14页，共26页，星期日，2025年，2月5日裂纹长度与间距分布函数的变化裂纹连通后，裂纹长度和间距的分布函数p(c)和f(a)将改变Heavisidestep函数连通概率没有连通的裂纹连通的裂纹第15页，共26页，星期日，2025年，2月5日扩展后的微裂纹长度及间距的期望值第16页，共26页，星期日，2025年，2月5日计算流程图a0+c0ccr1f1(a)+p1(c)+ccr2KAKIC?σ∞+f(a)+p(c)KAKIC?Repeat第17页，共26页，星期日，2025年，2月5日重复n次后裂纹长度第k次扩展后的裂纹间距期望值第k次扩展的临界裂纹长度总循环次数M第18页，共26页，星期日，2025年，2月5日破坏概率根据WLT以及有关的统计学知识，我们可以得到材料的破坏概率Pfail为1（M＝1）（M1）其中第19页，共26页，星期日，2025年，2月5日计算力学课程报告作者：张鹏