CN109470549B - 增材制造材料p-s-n曲线表征方法及其应用 - Google Patents

增材制造材料p-s-n曲线表征方法及其应用 Download PDF

Info

Publication number
CN109470549B
CN109470549B CN201811042140.XA CN201811042140A CN109470549B CN 109470549 B CN109470549 B CN 109470549B CN 201811042140 A CN201811042140 A CN 201811042140A CN 109470549 B CN109470549 B CN 109470549B
Authority
CN
China
Prior art keywords
curve
fatigue
distribution
life
additive manufacturing
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
Application number
CN201811042140.XA
Other languages
English (en)
Other versions
CN109470549A (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 CN201811042140.XA priority Critical patent/CN109470549B/zh
Publication of CN109470549A publication Critical patent/CN109470549A/zh
Application granted granted Critical
Publication of CN109470549B publication Critical patent/CN109470549B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00Investigating strength properties of solid materials by application of mechanical stress
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0058Kind of property studied
    • G01N2203/0069Fatigue, creep, strain-stress relations or elastic constants
    • G01N2203/0073Fatigue
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02Details not specific for a particular testing method
    • G01N2203/0202Control of the test
    • G01N2203/021Treatment of the signal; Calibration
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02Details not specific for a particular testing method
    • G01N2203/06Indicating or recording means; Sensing means
    • G01N2203/067Parameter measured for estimating the property

Abstract

本发明涉及材料疲劳性能表征领域,特别涉及增材制造材料P‑S‑N曲线表征方法及其应用。采用双峰对数正态分布描述给定应力水平下的疲劳寿命分布特性的模型,并建立了分布参数估计方法;在此基础上,建立了疲劳P‑S‑N曲线参数估计方法。获取给定应力水平下的对疲劳试验结过,通过双峰对数正态分布的数据处理,并建立分布参数估计方法,形成P‑S‑N的曲线表征描述材料的疲劳寿命分布特征,比单峰的对数正态分布模型更为合理,能更好的反映由于工艺特性等因素导致的疲劳寿命分散性问题,在有限试件试验的基础上获得的结果更合理。

Description

增材制造材料P-S-N曲线表征方法及其应用
技术领域
本发明涉及材料疲劳性能表征领域,特别涉及增材制造材料 P-S-N曲线表征方法及其应用。
背景技术
制造是国内外关注的焦点,有明显优势:相对传统工艺,在前期 研发和样品试制阶段,可以实现低成本、高效率、便于修改;在复杂 件生产中如模具制造、产品直接成型,可优化设计,降低成本,缩短 周期,实现原来无法实现的产品结构,亦可减少零件数量、提高系统 可靠性;对破损的部件可以实现修复,不需要重复生产,利用率大大 提高。因此,制造成为了国内外研究热点,增材制造金属材料/结构 在航空航天上的应用越来越广泛。
结构强度性能是增材制造材料/结构能否用于结构的主要考核指 标,目前,随着制造工艺(含后续热处理)的改进,增材制造材料/结 构的静力学性能有了很大提高,不低于甚至超过原材料板材、锻件的 性能,为在主承力结构上的应用奠定了良好的基础;
现代飞机结构按耐久性/损伤容限思想设计。在使用过程中,结 构要承受大量的交变载荷作用,交变载荷作用下的疲劳/断裂是结构 最主要的失效形式之一,因此,需要评估增材制造金属材料/结构的 疲劳性能。其中,疲劳性能测试是评估材料性能的重要手段。按相关 标准要求,需要测试得到材料的P-S-N曲线,即给定应力比R下应力峰 值σmax与对应一定可靠度P的可靠寿命NP的关系曲线。按目前规范和 标准要求,通常要进行成组疲劳试验,并认为给定应力水平下的寿命 服从对数正态分布,进行参数估计,确定P-S-N曲线参数。但增材制 造金属材料由于其材料和成型工艺特性,疲劳性能存在如下特点:
1)材料内部普遍存在初始缺陷,如气泡、孔洞等,由于工艺特 点,上述缺陷目前无法避免,而传统的无损检测方法往往无法发现上 述缺陷。上述缺陷在材料内部随机分布,导致其疲劳失效模式与常规 金属材料不同,常规金属材料由于工艺相对成熟与稳定,同一批次试 件的疲劳失效模式基本相同。但对增材制造金属材料,即使是同一批 次加工的试件,其疲劳失效模式也不同,可分为两类:一类为材料内 部微小缺陷导致的失效,通常从缺陷处萌生裂纹导致破坏;一类是与 常规材料一致的破坏形式,即内部没有明显可见缺陷。由于疲劳失效 模式不同,在进行疲劳性能表征时,采用传统的方法(认为给定应力 水平下的寿命值来源于同一个母体,采用对数正态分布描述其分布) 往往是不合适的。
2)增材制造工艺分散性大。增材制造金属材料的性能受原材料、 设备、工艺流程、环境条件等的影响很大,不同批次加工得到的试件 的疲劳性能存在明显差异,即使严格按工艺标准进行加工,也会出现 明显高于常规材料的工艺分散性。
由于上述特点,在采用已有的P-S-N曲线描述方法表征增材制造 材料的疲劳P-S-N曲线时,会出现很多问题,如:
1)对同一批次的材料试件,由于缺陷的随机分布,部分试件可 能存在缺陷,部分试件可能碰巧没有缺陷,导致疲劳寿命不同,不属 于同一个母体。采用传统的单峰分布函数(如单峰对数正态分布)描述 给定应力水平下的疲劳寿命分布不适用。
2)由于给定应力水平下的疲劳寿命分散性大,采用传统的基于 单峰分布随机变量模型的P-S-N曲线描述增材制造金属材料的疲劳性 能时,高可靠度下的可靠寿命会很短,难以满足结构设计要求,限制 了增材制造结构的应用。
3)由于分散性大,在进行疲劳性能测试时,试件数往往比常规 材料疲劳试验试件数要大得多,但这样必然造成时间和经济成本增 加。
4)由于客观存在的工艺分散性,不同批次材料的疲劳性能分散 性会明显高于常规材料,采用传统的疲劳P-S-N曲线也难以描述这种 特点。
目前,在采用传统的疲劳P-S-N曲线描述增材制造金属材料的疲 劳性能时,由于不符合其工艺和疲劳失效特点,造成疲劳性能分散性 大,可靠寿命低,大大限制了增材制造金属材料/结构的应用。从调 查到的资料来看,目前增材制造结构主要用于不承力或次承力结构, 对于主承力结构,迫切需要解决增材制造金属材料的疲劳性能P-S-N 曲线表征方法。
发明内容
(一)要解决的技术问题
本发明的目的是提供增材制造材料P-S-N曲线表征方法及其应 用,解决单峰表征结果不理想的问题。
(二)技术方案
为了解决上述技术问题,本发明提供一种增材制造材料P-S-N曲 线表征方法,其包括:
在预设应力比R下对n个试件进行多个应力作用下的疲劳试验, 获得与不同应力对应的疲劳寿命,n>1;
采用极大似然估计法建立似然函数
Figure BDA0001792367690000031
对每一个应 力下的疲劳寿命样本进行处理,迭代求解得到参数α,μ1,σ1,μ2,σ2, 其中,α为权重,0≤α≤1,μ1、μ2分别为2个分布的数学期望,σ1、σ2分别为双峰分布的对数寿命标准差;
利用
Figure BDA0001792367690000041
P(N≥NP)=P,xP=lgNP, 采用数值解法,及参数α,μ1,σ1,μ2,σ2求解可靠寿命NP,P为可靠 度,
Figure BDA0001792367690000042
采用双峰对数正态分布描述给定应力水平下对应可靠度P的可 靠寿命NP随应力峰值σmax变化的P-S-N曲线。
在一些实施例中,优选为,P-S-N曲线描述中,采用幂函数式进 行描述,所述幂函数为:
Figure BDA0001792367690000043
其中P为可靠度;σmax为应力峰值;mP为幂;NP为 可靠寿命;CP为曲线参数。
在一些实施例中,优选为,mP的求解公式为:
Figure BDA0001792367690000044
在一些实施例中,优选为,曲线参数CP的求解公式为:
Figure BDA0001792367690000045
在一些实施例中,优选为,所述双峰对数正态分布为:
概率密度函数f(x)=αf1(x)+(1-α)f2(x)
分布函数
Figure BDA0001792367690000046
设给定应力水平下的疲劳寿命为N,用随机变量Y表示,取其 对数为x=lgY。
在一些实施例中,优选为,在疲劳试验中,试件的破坏形式包括: 表面及亚表面起裂;内部起裂。
本发明还提供了一种上述P-S-N曲线表征方法的应用,其应用的 材料包括:铸件、C/C、C/SiC、颗粒增强金属材料。
(三)有益效果
本发明提供的技术方案中获取给定应力水平下的对疲劳试验结 过,通过双峰对数正态分布的数据处理,并建立分布参数估计方法, 形成P-S-N的曲线表征描述材料的疲劳寿命分布特征,比单峰的对数 正态分布模型更为合理,能更好的反映由于工艺特性等因素导致的疲 劳寿命分散性问题,在有限试件试验的基础上获得的结果更合理。
附图说明
图1为本发明增材制造材料疲劳P-S-N曲线表征方法示意图;
图2为本发明疲劳试验采用的试件图;
图3为本发明疲劳试验结果示意图;
图4a为本发明表面起裂类试件的破坏形式示意图;
图4b为本发明内部起裂类试件的破坏形式示意图;
图5为720MPa应力水平下单峰模型的疲劳寿命疲劳分布直方图 (FDH)和概率密度曲线;
图6为760MPa应力水平下单峰模型的疲劳寿命疲劳分布直方图 (FDH)和概率密度曲线;
图7为800MPa应力水平下单峰模型的疲劳寿命疲劳分布直方图 (FDH)和概率密度曲线;
图8为720MPa双峰模型的疲劳寿命疲劳分布直方图(FDH)和概率 密度曲线;
图9为760MPa双峰模型的疲劳寿命疲劳分布直方图(FDH)和概率 密度曲线;
图10为800MPa双峰模型的疲劳寿命疲劳分布直方图(FDH)和概 率密度曲线;
图11为720MPa疲劳寿命疲劳分布直方图(FDH)和概率密度曲线 对比图;
图12为760MPa疲劳寿命疲劳分布直方图(FDH)和概率密度曲线 对比图;
图13为800MPa疲劳寿命疲劳分布直方图(FDH)和概率密度曲线 对比图;
图14为单峰模型和双峰模型的S-N曲线拟合图;
图15为90%可靠度的P-S-N曲线;
图16为95%可靠度的P-S-N曲线;
图17为99%可靠度下的P-S-N曲线;
图18为99.9%可靠度下的P-S-N曲线。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细 描述。以下实例用于说明本发明,但不用来限制本发明的范围。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定, 术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接, 也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电 连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个 元件内部的连通。“第一”“第二”“第三”“第四”不代表任何的序列关系, 仅是为了方便描述进行的区分。对于本领域的普通技术人员而言,可 以具体情况理解上述术语在本发明中的具体含义。“当前”在执行某动 作之时的时刻,文中出现多个当前,均为随试件流逝中实时记录。
基于现有单峰表征结果不理想的问题,本发明给出了P-S-N曲线 表征方法及其应用。
下面将通过基础设计、扩展设计及替换设计对产品、方法等进行 详细描述。
本发明提供一种P-S-N曲线表征方法,结合具体试验进行说明, 如图1所示,其包括:
步骤110,准备疲劳测试试件;
由于本表征方法可适用于铸件、C/C、C/SiC、颗粒增强金属材料 等,因此,本例中以TA15钛合金为例,采用激光选区熔化成形工 艺制造,再机械加工成图2所示标准圆棒试件。试件表面粗糙度 Ra=0.8,同轴度为0.03,垂直度为0.04。
步骤120,在预设应力比R下对n个试件进行多个应力作用下对 试件进行疲劳试验,获得与不同应力对应的疲劳寿命,n>1;
在该疲劳试验中应力分别取720MPa、760MPa、800MPa,应力 比取R=0.1,获得成组疲劳试验结果见图3。其中试件为多组,对应 不同的应力。在其他的实施例中,可以将增加应力的数目,获取更多 应力下的疲劳寿命结果。
试件的破坏形式主要分为两种:一种为表面及亚表面起裂,见 4a,与常规材料圆棒试件破坏形式一致;一种为内部起裂,见图4b, 从试件内部的气泡处萌生裂纹,裂纹扩展,最后断裂。
步骤130,构建疲劳寿命分布的双峰模型
同一批次加工的增材制造金属材料存在不同的疲劳失效形式,成 组疲劳寿命试验结果包含了多种失效特点,考虑到疲劳寿命的随机特 性,可采用双峰对数正态分布描述疲劳寿命分布。
利用如下方法对三种应力水平下的疲劳试验得到的疲劳寿命数 据进行双峰对数正态分布模型的参数估计。
设给定应力水平下的疲劳寿命为N,用随机变量Y表示,取其对 数为X=lgY,采用如下由2个对数正态分布函数线性加权组成的双峰 分布函数描述X的分布
概率密度函数f(x)=αf1(x)+(1-α)f2(x)
分布函数
Figure BDA0001792367690000081
式中,α为权重,0≤α≤1,一般情况下α与应力水平无关。f1(x)、 f2(x)为两个概率密度函数,通常可取为对数正态分布,见式(2),对应 的分布函数见式(3)。其中,μ1、μ2分别为2个分布的数学期望;σ1、σ2分别为2个分布的对数寿命标准差。f1(x)、f2(x)均为单峰对数正态分布。 当α=0或α=1时式(2)和式(3)退化成如式(1)所示单峰对数正态分布。
Figure BDA0001792367690000082
Figure BDA0001792367690000083
对给定应力水平下的疲劳寿命样本,采用极大似然估计法(MLE) 估计参数α,μ1,σ1,μ2,σ2
根据MLE原理建立似然函数
Figure BDA0001792367690000084
与单峰对数正态分布函数相比,式(7)由2个对数正态分布加权得 到,概率密度函数有2个峰,形成双峰分布函数。
步骤140,双峰模型的参数估计
设某应力水平下进行了n个试件的疲劳试验,获得了n个疲劳寿命 数据。
(1)分布参数估计
对某给定应力水平下的疲劳寿命样本,采用极大似然估计法(MLE)估计参数α,μ1,σ1,μ2,σ2
根据MLE原理建立似然函数
Figure BDA0001792367690000091
方程两侧取对数,得到对数似然函数:
Figure BDA0001792367690000092
将(6)分别对α,μ1,σ1,μ2,σ2求导并令其为0得到似然方程:
Figure BDA0001792367690000093
整理得:
Figure BDA0001792367690000094
式(8)为非线性方程组,采用Newton-Raphson迭代法求解上述非 线性方程组。记FL(θ)=(fL,1(θ),fL,2(θ),fL,3(θ),fL,4(θ),fL,5(θ))T=0,其中 θ=(α,μ1122)T。假设已迭代到第k次,在点θ(k)=(α(k)1 (k)1 (k)2 (k)2 (k))T对 fL,1,fL,2,fL,3,fL,4,fL,5进行Taylor展开,并忽略高阶项,得到式(9)。令 LT(θ)=(lT,1(θ),lT,2(θ),lT,3(θ),lT,4(θ),lT,5(θ))T,取LT(θ)作为FT(θ)的近似,则LT(θ)=0的 根就是FT(θ)=0的根的近似。
Figure BDA0001792367690000101
令LT(k+1))=0,则,
Figure BDA0001792367690000102
其中Jacobi矩阵J为FT(θ)在θ(k)的导数:
Figure BDA0001792367690000103
则得到如下迭代格式
Figure BDA0001792367690000111
即,θ(k+1)=θ(k)-J-1FL(k)),k=0,1,2…。
在实际计算中,上述迭代格式需要求解Jacobi矩阵的逆矩阵,计 算中采用如式(13)中格式。
Figure BDA0001792367690000112
求解后即得到α,μ1,σ1,μ2,σ2
(2)可靠寿命估计
取可靠度为P,即可靠寿命为NP,则NP满足:
P(N≥NP)=P
若记xP=lgNP,由式(1)则有
Figure BDA0001792367690000113
采用数值解法即可由估计的分布参数α,μ1,σ1,μ2,σ2及给定P值 计算对应的xp,则
Figure BDA0001792367690000115
步骤150,描述疲劳P-S-N曲线
采用幂函数式描述给定应力水平下对应可靠度P的可靠寿命NP随应力峰值σmax变化的P-S-N曲线
Figure BDA0001792367690000114
式中,P为可靠度;σmax为应力峰值;mP为幂;NP为可靠寿命;CP为曲线参数。
当取P=50%时,为中值P-S-N曲线,简称为S-N曲线。
步骤160,P-S-N曲线参数估计
基于某应力比R下三种或三种以上应力水平σmax,j(j=1,…,k,k≥3) 下的成组疲劳试验,按2估计每种应力水平下的分布参数 α,μ1,σ1,μ2,σ2和给定可靠度要求P下的可靠寿命NP,j(j=1,…,k), 得到的各个应力水平下的可靠寿命(NP)和应力峰值(σmax)数据对(NP, σmax)j(j=1,…,k)。
将式(15)两边取对数,由(NP,σmax)j(j=1,…,k)数据采用线性回 归方法估计得到相关参数,估计式如下:
Figure BDA0001792367690000121
将求取的参数带入步骤150,获得疲劳P-S-N曲线。
下面,依然基于上述的疲劳样本测试,进行单峰模型和双峰模型 的对比,以验证双峰模型增材制造材料P-S-N曲线表征中的准确性。
图2为应力比R=0.1,应力分别为720MPa,760MPa,800MPa 的疲劳试验结果。
对该试验结果进行处理分析:
分析一:单峰对数正态分布模型表征S-N曲线
假定疲劳寿命样本来源于同一个母体,采用式(17)所示对数正态 分布描述给定应力水平下的疲劳寿命分布,设给定应力水平下的疲劳 寿命为N,用随机变量Y表示,取X=lgY。有:
Figure BDA0001792367690000122
式中,μ、σ为随机变量X的数学期望和对数寿命标准差。
采用单峰对数正态分布描述疲劳寿命分布,按式(18)、式(19)进 行参数估计和统计分析。
设给定应力水平下的疲劳寿命样本为yi(i=1,…,n),按极大似然 法(MLE)估计得到分布参数,结果见式(18)。
Figure BDA0001792367690000131
1)中值寿命和对数寿命标准差
中值寿命N50和可靠寿命NP按式(19)估计
Figure BDA0001792367690000132
式中,P为可靠度;up为标准正态分布分位点。
3个应力水平下的中值寿命和对数寿命标准差估计值见表1。
表1中值寿命和对数寿命标准差
应力峰值σ<sub>max</sub>/MPa 对数中值寿命 中值寿命/cycle 对数寿命标准差 试件数
720 5.66 457088 0.316 15
760 5.33 213796 0.348 17
800 4.98 95499 0.370 22
从表中可以看出,虽然每个应力水平下的有效试件数均远大于相 关标准要求,但3种应力水平下的对数寿命标准差均在0.3以上,为 《材料性能手册》中TA15钛合金板材和锻件圆棒试件的对数寿命标 准差的1.5倍以上。
2)单峰对数正态分布模型的概率密度曲线图
将三种应力水平下的疲劳寿命疲劳分布直方图(FDH)和概率密度 曲线绘制入图5-7。
3)分析
从5-7中可以看出,三种应力水平下,对数正态分布的概率密度 曲线与FDH的分布情况均有较大的差别。三种应力水平下的FDH均 呈现除了较为明显的双峰现象,因此采用传统单峰对数正态分布对增 材制造钛合金的疲劳寿命进行描述的无法取得良好的拟合效果。
4)基于单峰分布的S-N曲线
将表1中三个峰值应力和对应的中值寿命(σmax,N50)j(j=1,…,3) 代入式(16),估计得到基于单峰分布的幂函数式S-N曲线:
Figure BDA0001792367690000141
单双峰分布的S-N曲线拟合图件图14。
分析二,双峰对数正态分布模型
1)采用步骤160的方法对三种应力水平下了疲劳试验得到的疲 劳寿命数据进行了双峰对数正态分布模型的参数估计,结果见表2。 三种应力水平下的概率密度曲线见图8-10。
表2双峰对数正态分布模型分布参数
应力水平/MPa α μ<sub>1</sub> σ<sub>1</sub> μ<sub>2</sub> σ<sub>2</sub>
720 0.508 5.43 0.235 5.90 0.144
760 0.296 4.87 0.143 5.53 0.158
800 0.528 4.69 0.138 5.31 0.242
2)基于双峰分布的S-N曲线
取P=50%,计算对应的N50,将三个峰值应力和对应的中值寿命 (σmax,N50)j(j=1,…,3)代入式(16),估计得到基于双峰分布的幂函数式 S-N曲线:
Figure BDA0001792367690000142
单双峰分布的S-N曲线拟合图件图14。
对比分析:单峰模型和双峰模型对比分析
将三种应力水平下,两种分布模型的概率密度曲线和疲劳寿命频 率分布直方图一并绘制入图11-图13,从图中可以看出,双峰对数正 态分布更能够反应疲劳寿命频率分布直方图的双峰特点,与单峰对数 正态分布相比,双峰对数正态分布能够更加真实的反应AM钛合金的 疲劳寿命分布特点。图中点虚线表示概率(0.13%,99.87%)分布带, 对应于对数正态分布±3σ分布带,从中可以看出采用双峰对数正态分 布的分散带远远窄于单峰对数正态分布。
对比分析二单双峰模型的P-S-N曲线表征
取常用的寿命可靠度要求,P=90%、95%、99%、99.9%,分别 计算基于单峰分布和双峰分布的NP,将(σmax,NP)j(j=1,2,3)数据列入 表9和表10,并分别估计对应的P-S-N曲线,见图15-图18。
表3单峰分布模型下应力峰值,可靠性与寿命NP的关系
Figure BDA0001792367690000151
表4双峰分布模型下应力峰值,可靠性与寿命NP的关系
Figure BDA0001792367690000152
由于双峰对数正态分布函数更为合理,得到的P-S-N曲线也更为 合理。从表3-表4和图15-图18中可以看出,采用双峰对数正态分布 描述疲劳寿命分布,高可靠度要求下P-S-N曲线明显高于基于单峰对 数正态分布的P-S-N曲线,即同样的应力水平下,当用更为合理的双 峰对数正态分布描述疲劳寿命分布时,其可靠寿命NP明显高于服从 单峰对数正态分布的NP值,为增材制造金属材料在主承力结构上的 应用奠定了良好的基础。
本发明为进行增材制造金属材料的疲劳性能表征,考虑增材制造 工艺的特点,提出了一种采用双峰对数正态分布描述给定应力水平下 的疲劳寿命分布特性的模型,并建立了分布参数估计方法;在此基础 上,建立了疲劳P-S-N曲线参数估计方法。
优点:
(1)采用双峰对数正态分布描述增材制造金属材料的疲劳寿命 分布特征,比单峰对数正态分布模型更为合理,能更好的反映由于工 艺特性导致的疲劳寿命分散性;
(2)基于单峰对数正态分布模型的P-S-N曲线过于保守,在考虑 可靠度要求时的疲劳寿命过低,造成结构寿命过低的情况,为了保证 结构的使用安全,则必须降低应力水平,从而导致结构重量增加。采 用本文方法得到的P-S-N曲线在高可靠度要求下明显高于基于单峰 对数正态模型的P-S-N曲线,结果更为科学合理;
(3)采用传统的单峰对数正态分布描述寿命分布时,由于分散性 过大,为了能够获得有效的试验结果,必须增加试件数。而采用本文 的方法,能用较少的试件数获得符合要求的试验结果;
(4)本方法可推广应用到铸件、C/C、C/SiC、颗粒增强金属材料 及其他在材料制备时容易产生初始缺陷的材料疲劳P-S-N曲线表征。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明, 凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等, 均应包含在本发明的保护范围之内。

Claims (7)

1.一种增材制造材料P-S-N曲线表征方法,其特征在于,包括:
准备试件;
在预设应力比R下对n个试件进行多个应力作用下的疲劳试验,获得与不同应力对应的疲劳寿命,n>1;
采用极大似然估计法建立似然函数
Figure FDA0002408700910000011
对每一个应力下的疲劳寿命样本进行处理,迭代求解得到参数α,μ1,σ1,μ2,σ2,其中,α为权重,0≤α≤1,μ1、μ2分别为2个分布的数学期望,σ1、σ2分别为双峰分布的对数寿命标准差;
利用
Figure FDA0002408700910000012
P(N≥NP)=P,xP=lgNP,采用数值解法,及参数α,μ1,σ1,μ2,σ2求解可靠寿命NP,P为可靠度,
Figure FDA0002408700910000013
采用双峰对数正态分布描述给定应力水平下对应可靠度P的可靠寿命NP随应力峰值σmax变化的P-S-N曲线。
2.如权利要求1所述的增材制造材料P-S-N曲线表征方法,其特征在于,P-S-N曲线描述中,采用幂函数式进行描述,所述幂函数为:
Figure FDA0002408700910000014
其中P为可靠度;σmax为应力峰值;mP为幂;NP为可靠寿命;CP为曲线参数。
3.如权利要求2所述的增材制造材料P-S-N曲线表征方法,其特征在于,mP的求解公式为:
Figure FDA0002408700910000021
NP,j(j=1,…,k)为给定可靠度要求P下的可靠寿命,σmax,j(j=1,…,k)为基于应力比R下应力水平,k为迭代次数。
4.如权利要求2所述的增材制造材料P-S-N曲线表征方法,其特征在于,曲线参数CP的求解公式为:
Figure FDA0002408700910000022
NP,j(j=1,…,k)为给定可靠度要求P下的可靠寿命,σmax,j(j=1,…,k)为基于应力比R下应力水平,k为迭代次数。
5.如权利要求1-4任一项所述的增材制造材料P-S-N曲线表征方法,其特征在于,所述双峰对数正态分布为:
概率密度函数f(x)=αf1(x)+(1-α)f2(x)
分布函数
Figure FDA0002408700910000023
设给定应力水平下的疲劳寿命为N,用随机变量Y表示,取其对数为x=lgY。
6.如权利要求1-4任一项所述的增材制造材料P-S-N曲线表征方法,其特征在于,在疲劳试验中,试件的破坏形式包括:表面及亚表面起裂;内部起裂。
7.一种权利要求1-4任一项所述增材制造材料P-S-N曲线表征方法的应用,其特征在于,应用材料包括:铸件、C/C、C/SiC、颗粒增强金属材料。
CN201811042140.XA 2018-09-07 2018-09-07 增材制造材料p-s-n曲线表征方法及其应用 Active CN109470549B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811042140.XA CN109470549B (zh) 2018-09-07 2018-09-07 增材制造材料p-s-n曲线表征方法及其应用

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811042140.XA CN109470549B (zh) 2018-09-07 2018-09-07 增材制造材料p-s-n曲线表征方法及其应用

Publications (2)

Publication Number Publication Date
CN109470549A CN109470549A (zh) 2019-03-15
CN109470549B true CN109470549B (zh) 2020-07-28

Family

ID=65663074

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811042140.XA Active CN109470549B (zh) 2018-09-07 2018-09-07 增材制造材料p-s-n曲线表征方法及其应用

Country Status (1)

Country Link
CN (1) CN109470549B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112464490B (zh) * 2020-12-10 2022-11-25 北京航空航天大学 增材钛合金的dfr确定方法
CN112949020B (zh) * 2020-12-10 2022-08-12 北京航空航天大学 增材钛合金的疲劳寿命的确定方法、系统、设备与介质
CN113742948B (zh) * 2021-08-23 2022-07-22 西安石油大学 一种超高强度抽油杆p-s-n曲线拟合新模型及方法
CN113536490B (zh) * 2021-08-27 2022-05-17 西安石油大学 一种抽油杆疲劳寿命预测方法及系统
CN115628998B (zh) * 2022-10-27 2023-09-19 深海技术科学太湖实验室 节点疲劳试验的可靠性确定方法、装置及存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101122560A (zh) * 2007-08-23 2008-02-13 南京航空航天大学 机械结构的裂纹扩展率和裂纹扩展寿命预测方法
CN102509023A (zh) * 2011-11-24 2012-06-20 北京航空航天大学 航天驱动组件综合应力加速寿命试验损伤累积模型的建模方法
CN104036120A (zh) * 2014-05-20 2014-09-10 北京航空航天大学 一种测定材料与构件疲劳s-n曲线性能的单点法
CN105716935A (zh) * 2016-01-28 2016-06-29 中国科学院力学研究所 预测试样尺寸对疲劳寿命影响的方法
CN105740625A (zh) * 2016-01-31 2016-07-06 太原科技大学 一种齿轮的实时剩余寿命预测方法
CN106383959A (zh) * 2016-09-23 2017-02-08 南京航空航天大学 一种基于最大熵模型的材料疲劳寿命的预测方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100324900B1 (ko) * 1999-05-25 2002-02-28 조양호 피로수명 예측 방법
JP4488957B2 (ja) * 2005-05-26 2010-06-23 財団法人鉄道総合技術研究所 疲労状態解析装置及び疲労状態解析プログラム

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101122560A (zh) * 2007-08-23 2008-02-13 南京航空航天大学 机械结构的裂纹扩展率和裂纹扩展寿命预测方法
CN102509023A (zh) * 2011-11-24 2012-06-20 北京航空航天大学 航天驱动组件综合应力加速寿命试验损伤累积模型的建模方法
CN104036120A (zh) * 2014-05-20 2014-09-10 北京航空航天大学 一种测定材料与构件疲劳s-n曲线性能的单点法
CN105716935A (zh) * 2016-01-28 2016-06-29 中国科学院力学研究所 预测试样尺寸对疲劳寿命影响的方法
CN105740625A (zh) * 2016-01-31 2016-07-06 太原科技大学 一种齿轮的实时剩余寿命预测方法
CN106383959A (zh) * 2016-09-23 2017-02-08 南京航空航天大学 一种基于最大熵模型的材料疲劳寿命的预测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
一种P-S-N曲线的快速测定方法;陈一龙 等;《汽车工艺与材料》;20041231(第7期);101-104页 *
估计三种常用疲劳应力-寿命模型P-S-N曲线的统一经典极大似然法;赵永翔 等;《应用力学学报》;20010331;第18卷(第1期);83-90页 *

Also Published As

Publication number Publication date
CN109470549A (zh) 2019-03-15

Similar Documents

Publication Publication Date Title
CN109470549B (zh) 增材制造材料p-s-n曲线表征方法及其应用
Molent A review of equivalent pre‐crack sizes in aluminium alloy 7050‐T7451
Markov et al. Development of a new process for forging plates using intensive plastic deformation
CN108613889B (zh) 基于循环寿命的钛合金钝缺口疲劳强度损失系数评估方法
CN104729962A (zh) Gh4169合金锻件晶粒度分析及预测方法
CN109142528B (zh) 一种高强度钛合金高温超高周疲劳寿命预测方法
CN113987681A (zh) 缺口-尺寸效应下耦合应变能梯度的结构疲劳寿命评估方法
CN104732003A (zh) 一种基于可靠性的磨削工艺评估方法
Fan et al. Competitive effect of residual stress and surface roughness on the fatigue life of shot peened S42200 steel at room and elevated temperature
CN109213963B (zh) 一种激光沉积成型钛合金材料疲劳寿命统计分析方法
CN113604649B (zh) 42CrMo合金结构钢空心轴类锻件热处理方法及其制备的空心轴类锻件
Zamani et al. On the influence of riveting process parameters on fatigue life of riveted lap joint
Liu et al. Efficient methods for time-dependent fatigue reliability analysis
CN112100806B (zh) 一种表征复合材料本构曲线的获取方法
Ball et al. A detailed evaluation of the effects of bulk residual stress on fatigue in aluminum
CN111175127A (zh) 一种基于零件变形程度控制的tc4合金机匣锻件验收方法
CN116384809A (zh) 一种工程实体质量分析评价方法
Christiansen et al. Multi-objective optimization of die geometry in ingot forging
Main et al. Outcomes of research into small fatigue crack nucleation and growth in AA7085-T7452
CN109870257B (zh) 一种板材厚度方向淬火残余应力分布预测方法
Corran et al. The development of anomaly distributions for machined holes in aircraft engine rotors
CN108664716B (zh) 一种基于超高周疲劳性能的高强钛合金热处理工艺评价模型
Shu et al. A comparative study on constitutive equations and artificial neural network model to predict high-temperature deformation behavior in Nitinol 60 shape memory alloy
Shao et al. Research on springback compensation for multi-point forming of corrosion-resistant aluminum alloy hyperbolic component
CN110565135B (zh) 薄壁零件的电镀方法

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