CN114117826B - 一种针对浮动摩擦片跨尺度疲劳寿命预测方法及装置 - Google Patents

一种针对浮动摩擦片跨尺度疲劳寿命预测方法及装置 Download PDF

Info

Publication number
CN114117826B
CN114117826B CN202210068797.3A CN202210068797A CN114117826B CN 114117826 B CN114117826 B CN 114117826B CN 202210068797 A CN202210068797 A CN 202210068797A CN 114117826 B CN114117826 B CN 114117826B
Authority
CN
China
Prior art keywords
fatigue life
stress
equivalent
sample
floating friction
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
CN202210068797.3A
Other languages
English (en)
Other versions
CN114117826A (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.)
Institute of Mechanics of CAS
Original Assignee
Institute of Mechanics of CAS
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 Institute of Mechanics of CAS filed Critical Institute of Mechanics of CAS
Priority to CN202210068797.3A priority Critical patent/CN114117826B/zh
Publication of CN114117826A publication Critical patent/CN114117826A/zh
Application granted granted Critical
Publication of CN114117826B publication Critical patent/CN114117826B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/04Ageing analysis or optimisation against ageing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Abstract

本发明公开了一种针对浮动摩擦片跨尺度疲劳寿命预测方法及装置,包括以下步骤:步骤S1、设置齿状等效试样,将浮动摩擦片服役过程中遭受到的载荷谱加载到齿状等效试样上,通过DIC分析技术分析得到齿状等效试样的齿根应力状态;步骤S2、将所述齿根应力状态进行分析,并利用参数公式将齿根应力状态转换为应力集中等效参数和最大主应力;步骤S3、设置标准试样,设计多组疲劳寿命实验,并基于疲劳寿命实验的实验结果得到疲劳寿命预测模型。本发明通过齿状等效试样的实验,保证了预测结果的准确性,通过标准试样的实验,保证了实验方法的简便性,做到了准确性与简便性的统一。

Description

一种针对浮动摩擦片跨尺度疲劳寿命预测方法及装置
技术领域
本发明涉及浮动摩擦片寿命预测技术领域,具体涉及一种针对浮动摩擦片跨尺度疲劳寿命预测方法及装置。
背景技术
对于浮动摩擦片的疲劳寿命预测是相关工程产品在服役中保持结构完整性、可靠性、稳定性关键技术之一。经典断裂疲劳理论和大量疲劳实验表明,浮动摩擦片的齿根位置由于应力集中,是裂纹萌生的主要区域。结合理论分析与浮动摩擦片的实际应用,要想对浮动摩擦片的疲劳寿命进行精确的预测,至少需要考虑以下几个因素:应力集中等效参数、最大主应力、残余应力、粗糙度、应力比、齿根的应力状态等。其中特别需要指出的是齿根的应力状态,经典断裂疲劳理论和大量疲劳实验表明,浮动摩擦片的齿根位置由于应力集中,是裂纹萌生的主要区域,因而,准确量化这一因素的影响对浮动摩擦片的疲劳寿命的精准预测至关重要。
当前并不具备方便的进行浮动摩擦片的疲劳寿命实验的设备,若考虑前文提到的多种因素,强行对其进行直接的多种类型的疲劳寿命实验,还会耗费大量的人力、物力和财力。而直接利用标准试样进行材料级别的实验,又无法考虑齿根的应力状态,导致预测结果不准。
发明内容
本发明的目的在于提供一种针对浮动摩擦片跨尺度疲劳寿命预测方法及装置,以解决现有技术中强行对其进行直接的多种类型的疲劳寿命实验,还会耗费大量的人力、物力和财力,而直接利用标准试样进行材料级别的实验,又无法考虑齿根的应力状态,导致预测结果不准的技术问题。
为解决上述技术问题,本发明具体提供下述技术方案:
一种针对浮动摩擦片跨尺度疲劳寿命预测方法,包括以下步骤:
步骤S1、设置齿状等效试样,将浮动摩擦片服役过程中遭受到的载荷谱加载到齿状等效试样上,通过DIC分析技术分析得到齿状等效试样的齿根应力状态,所述载荷谱由所述浮动摩擦片在服役过程中采集得到,所述齿状等效试样表征与浮动摩擦片具有等效齿状结构的试样,并用于等效获得浮动摩擦片的齿根应力状态;
步骤S2、将所述齿根应力状态进行分析,并利用参数公式将齿根应力状态转换为应力集中等效参数和最大主应力;
步骤S3、设置标准试样,基于所述应力集中等效参数和最大主应力,以及所述标准试样的残余应力、粗糙度、应力比为标准试样设计多组疲劳寿命实验,并基于疲劳寿命实验的实验结果得到疲劳寿命预测模型,所述标准试样用于等效获得应力集中等效参数、最大主应力、残余应力和粗糙度对浮动摩擦片疲劳特性的影响。
作为本发明的一种优选方案,所述将浮动摩擦片服役过程中遭受到的载荷谱加载到齿状等效试样上,通过DIC分析技术分析得到齿状等效试样的齿根应力状态,包括:
将所述载荷谱作为输入载荷,对齿状等效试样开展疲劳寿命实验,依托摄像机对疲劳寿命实验全过程进行监测;
通过DIC分析技术分析疲劳寿命实验的实验结果得到所述齿状等效试样的齿根应力状态。
作为本发明的一种优选方案,所述将所述齿根应力状态进行分析,并利用参数公式将齿根应力状态转换为应力集中等效参数和最大主应力,包括:
考量残余应力的影响,将所述齿根应力状态依托第一参数公式得到齿根处的等效应力水平、应力集中等效参数和最大主应力,所述第一参数公式为:
Figure 100002_DEST_PATH_IMAGE002
式中,σapp1表征最大主应力,Kt为应力集中等效参数,σeff表征等效应力水平,r0表征残余应力影响区深度,a表征齿根近似的椭圆形结构的长轴长,b表征齿根近似的椭圆形结构的短轴长,σ0表征载荷谱中施加在齿根处的外载荷, r表征残余应力影响区深度的积分变量,σres表征残余应力;
在整个疲劳寿命实验的时间段内,依托基于韦伯概率分布的第二参数公式得到基于等效应力水平的最大估计应力,所述第二参数公式为:
Figure 100002_DEST_PATH_IMAGE004
式中,σm表征最大估计应力,P表征概率密度;
考虑表面粗糙度影响,通过第三参数公式基于所述最大估计应力得到用于预测试样的最终应力水平,所述第三参数公式为:
Figure 100002_DEST_PATH_IMAGE006
式中,σ预测试样表征预测试样的最终应力水平,
Figure 100002_DEST_PATH_IMAGE008
表征第一修正系数,Ra表征粗糙度,所述预测试样表征待预测疲劳寿命的试样。
作为本发明的一种优选方案,所述基于疲劳寿命实验的实验结果得到疲劳寿命预测模型,包括:
通过改变所述应力集中等效参数和最大主应力,以及集合所述残余应力、粗糙度、应力比的值对标准试样进行多组的疲劳寿命实验;
通过分析多组的疲劳寿命实验的实验结果对所述应力集中等效参数和最大主应力,以及所述残余应力、粗糙度、应力比的参数影响进行评估,基于评估结果反推算出浮动摩擦片的疲劳寿命预测模型,所述疲劳寿命预测模型的函数表达式为:
Figure 100002_DEST_PATH_IMAGE010
式中,N标准试样表征标准试样的疲劳寿命,N预测试样表征预测试样的疲劳寿命,σ预测试样表征标准试样的最终应力水平,σ标准试样表征预测试样的最终应力水平,
Figure 100002_DEST_PATH_IMAGE012
表征第二修正参数。
作为本发明的一种优选方案,在利用所述疲劳寿命预测模型对遭受固定载荷谱的多个浮动摩擦片进行疲劳寿命的预测时,仅利用与对应的浮动摩擦片具有等效齿状结构的齿状等效试样进行疲劳寿命实验即可,所述多个浮动摩擦片具有固定形状且具有不同的残余应力和不同的粗糙度,包括:
选取多组具有固定形状且具有不同的残余应力和不同的粗糙度的浮动摩擦片,并为多组所述浮动摩擦片对应设置多组齿状等效试样;
分别对多组齿状等效试样进行疲劳寿命实验分析得到实验结果,以及对选定的标准试样进行疲劳寿命实验分析得到实验结果,并基于粗糙度和残余应力构建出多组浮动摩擦片中表征N标准试样和σ标准试样的非线性映射关系的多组疲劳寿命曲线,所述疲劳寿命曲线的函数表达式标记为
Figure 100002_DEST_PATH_IMAGE014
,式中,
Figure 100002_DEST_PATH_IMAGE016
表征第i组浮动摩擦片对应的标准试样的最终应力水平,
Figure 100002_DEST_PATH_IMAGE018
表征第i组浮动摩擦片对应的标准试样的疲劳寿命,
Figure 100002_DEST_PATH_IMAGE020
表征所述非线性映射关系的函数体,
Figure 100002_DEST_PATH_IMAGE022
表征第i组浮动摩擦片的粗糙度,
Figure 100002_DEST_PATH_IMAGE024
表征第i组浮动摩擦片的残余应力,i为计量常数;
基于疲劳寿命预测模型在疲劳寿命曲线中求得多组齿状等效试样的疲劳寿命,所述多组齿状等效试样的疲劳寿命的计算公式为:
Figure 100002_DEST_PATH_IMAGE026
式中,
Figure DEST_PATH_IMAGE027
表征第i组浮动摩擦片对应的齿状等效试样的疲劳寿命,
Figure 507695DEST_PATH_IMAGE016
表征第i组浮动摩擦片对应的齿状等效试样的最终应力水平;
将多组齿状等效试样的疲劳寿命等价作为多组浮动摩擦片的疲劳寿命,以实现对浮动摩擦片的疲劳寿命预测。
作为本发明的一种优选方案,所述多组齿状等效试样的实验结果和标准试样的实验结果在进行疲劳寿命预测模型运算时需要进行归一化处理。
作为本发明的一种优选方案,所述对所述浮动摩擦片进行疲劳寿命实验也可实现对所述疲劳寿命预测模型进行验证。
作为本发明的一种优选方案,本发明提供了一种根据所述的针对浮动摩擦片跨尺度疲劳寿命预测方法的预测装置,包括:
齿状等效试样实验模块,用于对齿状等效试样进行疲劳寿命实验及分析;
标准试样实验模块,用于对标准试样进行疲劳寿命实验及分析;
模型构建及验证模块,用于构建疲劳寿命预测模型及验证。
作为本发明的一种优选方案,所述齿状等效试样实验模块、标准试样实验模块和模型构建及验证模块中均集成有通讯模块和数据存储模块,所述通讯模块用于在所述齿状等效试样实验模块、标准试样实验模块和模型构建及验证模块中进行数据交互,所述数据存储模块用于在所述齿状等效试样实验模块、标准试样实验模块和模型构建及验证模块中进行数据存储。
作为本发明的一种优选方案,所述齿状等效试样与浮动摩擦片采用同种制作材料。
本发明与现有技术相比较具有如下有益效果:
本发明通过齿状等效试样的实验获得应力集中等效参数、最大主应力、残余应力、粗糙度,保证了预测结果的准确性,通过标准试样的实验获得了齿根的应力状态,保证了实验方法的简便性,做到了准确性与简便性的统一,最终形成了一套浮动摩擦片跨尺度疲劳寿命预测方法及模型,模型在应用上还具有很好的广泛适应性。
附图说明
为了更清楚地说明本发明的实施方式或现有技术中的技术方案,下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍。显而易见地,下面描述中的附图仅仅是示例性的,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图引伸获得其它的实施附图。
图1为本发明实施例提供的针对浮动摩擦片跨尺度疲劳寿命预测方法流程图;
图2为本发明实施例提供的等效应力水平示意图;
图3(1)为本发明实施例提供的基于粗糙度的标准试样疲劳寿命曲线示意图;
图3(2)为本发明实施例提供的基于残余应力的标准试样疲劳寿命曲线示意图
图4为本发明实施例提供的预测装置结构框图;
图5为本发明实施例提供的实验件模型结构示意图。
图中的标号分别表示如下:
1-齿状等效试样实验模块;2-标准试样实验模块;3-模型构建及验证模块。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,浮动摩擦片的寿命预测需要考虑多种因素,其中必不可少的包含上文中提到的应力集中等效参数、最大主应力、残余应力、粗糙度、应力比、齿根的应力状态等。在本发明提出的方法中,将参数主要分为三类。第一类参数为结构相关参数,如齿根的应力状态,第二类参数为材料相关参数,如应力集中等效参数、最大主应力、残余应力、粗糙度,第三类参数为可通过公式转化的参数如应力比等,若考虑前文提到的多种因素,强行对其进行直接的多种类型的疲劳寿命实验,还会耗费大量的人力、物力和财力。而直接利用标准试样进行材料级别的实验,又无法考虑齿根的应力状态,导致预测结果不准,因此本发明提供了一种针对浮动摩擦片跨尺度疲劳寿命预测方法,通过齿状等效试样的实验获得第二类参数,保证了预测结果的准确性,通过标准试样的实验获得了第一类参数,保证了实验方法的简便性,做到了准确性与简便性的统一。
一种针对浮动摩擦片跨尺度疲劳寿命预测方法,包括以下步骤:
步骤S1、设置齿状等效试样,将浮动摩擦片服役过程中遭受到的载荷谱加载到齿状等效试样上,通过DIC分析技术分析得到齿状等效试样的齿根应力状态,DIC分析技术即数字图像相关分析技术,所述载荷谱由所述浮动摩擦片在服役过程中采集得到,所述齿状等效试样表征与浮动摩擦片具有等效齿状结构的试样,并用于等效获得浮动摩擦片的齿根应力状态;
所述将浮动摩擦片服役过程中遭受到的载荷谱加载到齿状等效试样上,通过DIC分析技术分析得到齿状等效试样的齿根应力状态,包括:
将所述载荷谱作为输入载荷,对齿状等效试样开展疲劳寿命实验,依托摄像机对疲劳寿命实验全过程进行监测;
通过DIC分析技术分析疲劳寿命实验的实验结果得到所述齿状等效试样的齿根应力状态。
步骤S1的操作主要对象为一个齿状等效试样,这也是本实施例有别于传统方法的一个重要特点。传统的方法主要为用材料的疲劳力学特性预测结构的疲劳性能,有一些还会结合上有限元的方法。本方法考虑了传统方法中忽视的齿根的应力状态这一重要因素。第一步中的输入为浮动摩擦片遭受的载荷谱,载荷谱是在浮动摩擦片服役过程中采集到的。将该载荷谱作为输入载荷,对齿状等效试样开展疲劳寿命实验,依托摄像机对实验全过程进行监测,通过对DIC结果的分析,可以得到齿状等效试样的齿根应力状态,即第一类参数。
步骤S2、将所述齿根应力状态进行分析,并利用参数公式将齿根应力状态转换为应力集中等效参数和最大主应力,即为第二类参数;
步骤S2的操作为一个分析过程,主要的目的为将所有的参数统一到第二类参数进行表征。具体包括两个方面:第一个方面为将齿根的应力状态依托第一参数公式基于图2在考虑残余应力的影响下得到齿根处的等效应力水平,其中r0是残余应力影响区深度,齿根近似为椭圆形结构(a为长轴长,b为短轴长),Kt为应力集中系数,σ0为施加的外载荷。
所述将所述齿根应力状态进行分析,并利用参数公式将齿根应力状态转换为应力集中等效参数和最大主应力,包括:
考量残余应力的影响,将所述齿根应力状态依托第一参数公式得到齿根处的等效应力水平、应力集中等效参数和最大主应力,所述第一参数公式为:
Figure DEST_PATH_IMAGE029
式中,σapp1表征最大主应力,Kt为应力集中等效参数,σeff表征等效应力水平,r0表征残余应力影响区深度,a表征齿根近似的椭圆形结构的长轴长,b表征齿根近似的椭圆形结构的短轴长,σ0表征载荷谱中施加在齿根处的外载荷, r表征残余应力影响区深度的积分变量,σres表征残余应力;
在整个疲劳寿命实验的时间段内,依托基于韦伯概率分布的第二参数公式得到基于等效应力水平的最大估计应力,所述第二参数公式为:
Figure DEST_PATH_IMAGE031
式中,σm表征最大估计应力,P表征概率密度;
考虑表面粗糙度影响,通过第三参数公式基于所述最大估计应力得到用于预测试样的最终应力水平,所述第三参数公式为:
Figure DEST_PATH_IMAGE032
式中,σ预测试样表征预测试样的最终应力水平,
Figure 715954DEST_PATH_IMAGE008
表征第一修正系数,Ra表征粗糙度,所述预测试样表征待预测疲劳寿命的试样。
步骤S3、设置标准试样,标准试样选取标准无具体要求,一般只要满足疲劳实验的实验要求即可,如实验机无法给出具体尺寸要求,推荐采用GB/T 15248-2008中6.1.3 板材试样,基于所述应力集中等效参数和最大主应力,以及所述标准试样的残余应力、粗糙度、应力比为标准试样设计多组疲劳寿命实验,并基于疲劳寿命实验的实验结果得到疲劳寿命预测模型,所述标准试样用于等效获得应力集中等效参数、最大主应力、残余应力和粗糙度对浮动摩擦片疲劳特性的影响。
步骤S3的操作主要对象为标准试样。通过上述步骤S2的操作,将第一类参数和第三类参数都转化到了第二类参数中进行表征。通过改变这些第二类参数的值并对标准试样进行一系列的疲劳实验,通过分析得到的一系列数据,进而可以评估上述不同参数的影响,进而可以反推出浮动摩擦片的疲劳寿命预测模型。利用本实施例提到的步骤S1至S3中的实验和分析,结合疲劳寿命预测模型,可以对遭受固定载荷谱的固定形状的各种不同的残余应力和不同的粗糙度的浮动摩擦片的疲劳寿命进行预测。通过对一个完整的浮动摩擦片进行实验也可以对本疲劳寿命预测模型进行验证。
所述基于疲劳寿命实验的实验结果得到疲劳寿命预测模型,包括:
通过改变所述应力集中等效参数和最大主应力,以及结合所述残余应力、粗糙度、应力比的值对标准试样进行多组的疲劳寿命实验;
通过分析多组的疲劳寿命实验的实验结果对所述应力集中等效参数和最大主应力,以及所述残余应力、粗糙度、应力比的参数影响进行评估,基于评估结果反推算出浮动摩擦片的疲劳寿命预测模型,所述疲劳寿命预测模型的函数表达式为:
Figure DEST_PATH_IMAGE033
式中,N标准试样表征标准试样的疲劳寿命,N预测试样表征预测试样的疲劳寿命,σ预测试样表征标准试样的最终应力水平,σ标准试样表征预测试样的最终应力水平,
Figure 700965DEST_PATH_IMAGE012
表征第二修正参数。
本实施例提到的模型在应用上还具有很好的广泛适应性。在本模型的应用中,基于已有数据,如想对材质相同的其他载荷谱或其他形状的另一种浮动摩擦片进行疲劳寿命预测,只需利用该种浮动摩擦片的等效齿状结构进行实验即可,而无需再进行标准试样的实验即可得到对应的疲劳寿命预测模型。
在利用所述疲劳寿命预测模型对遭受固定载荷谱的多个浮动摩擦片进行疲劳寿命的预测时,仅利用与对应的浮动摩擦片的等效齿状结构的齿状等效试样进行疲劳寿命实验即可,所述多个浮动摩擦片具有固定形状且具有不同的残余应力和不同的粗糙度,包括:
选取多组具有固定形状且具有不同的残余应力和不同的粗糙度的浮动摩擦片,并为多组所述浮动摩擦片对应设置多组齿状等效试样;
如图3(1)和3(2)所示,分别对多组齿状等效试样进行疲劳寿命实验分析得到实验结果,以及对选定的标准试样进行疲劳寿命实验分析得到的实验结果,并构建出多组浮动摩擦片中基于粗糙度和残余应力对应表征N标准试样和σ标准试样非的线性映射关系的多组疲劳寿命曲线,所述疲劳寿命曲线的函数表达式标记为
Figure 726690DEST_PATH_IMAGE014
,式中,
Figure 134669DEST_PATH_IMAGE016
表征第i组浮动摩擦片对应的标准试样的最终应力水平,
Figure DEST_PATH_IMAGE034
表征第i组浮动摩擦片对应的标准试样的疲劳寿命,
Figure DEST_PATH_IMAGE035
表征所述非线性映射关系的函数体,
Figure 710400DEST_PATH_IMAGE022
表征第i组浮动摩擦片的粗糙度,
Figure DEST_PATH_IMAGE036
表征第i组浮动摩擦片的残余应力,i为计量常数;
基于疲劳寿命预测模型在疲劳寿命曲线中求得多组齿状等效试样的疲劳寿命,所述多组齿状等效试样的疲劳寿命的计算公式为:
Figure DEST_PATH_IMAGE037
式中,
Figure 51383DEST_PATH_IMAGE018
表征第i组浮动摩擦片对应的齿状等效试样的疲劳寿命,
Figure DEST_PATH_IMAGE038
表征第i组浮动摩擦片对应的齿状等效试样的最终应力水平;
将多组齿状等效试样的疲劳寿命等价作为多组浮动摩擦片的疲劳寿命,以实现对浮动摩擦片的疲劳寿命预测。
所述多组齿状等效试样的实验结果和标准试样的实验结果在进行疲劳寿命预测模型运算时需要进行归一化处理。
所述对所述浮动摩擦片进行疲劳寿命实验也可实现对所述疲劳寿命预测模型进行验证。
如图4所示,基于上述针对浮动摩擦片跨尺度疲劳寿命预测方法,本发明提供了一种装置,包括:
齿状等效试样实验模块1,用于对齿状等效试样进行疲劳寿命实验及分析;
标准试样实验模块2,用于对标准试样进行疲劳寿命实验及分析;
模型构建及验证模块3,用于构建疲劳寿命预测模型及验证。
所述齿状等效试样实验模块、标准试样实验模块和模型构建及验证模块中均集成有通讯模块和数据存储模块,所述通讯模块用于在所述齿状等效试样实验模块、标准试样实验模块和模型构建及验证模块中进行数据交互,所述数据存储模块用于在所述齿状等效试样实验模块、标准试样实验模块和模型构建及验证模块中进行数据存储。
所述齿状等效试样与浮动摩擦片采用同种制作材料。实验件模型如图5所示:
(1)齿状等效试样示意图:示意图为基本形状,对该件进行实验并进行分析,可获得齿根处的应力状态,获得第一类参数。
(2)标准试样形状示意图:示意图为基本形状,实际应用过程中对其施加人工缺口等各种特征,利用此类型试件可获得第二类参数对疲劳特性的影响。
(3)浮动摩擦片示意图:示意图为某类型浮动摩擦片基本形状,利用本专利提出的方法,可获得浮动摩擦片的疲劳寿命预测模型,对该实际浮动摩擦片进行疲劳实验,可验证本方法和模型的可靠性。
预测试样为一个完整的浮动摩擦片,属于一个结构,是最终要预测的试样,但直接对其进行疲劳寿命实验较为复杂;标准试样是一个常规的试样,该试样可以方便的进行疲劳寿命试验,本发明目的为利用标准试验对预测试样的疲劳寿命进行预测。但是由于预测试样本身具有一定的结构特征,采用标准试样对预测试样进行寿命预测需要了解预测试样的受力状况。直接采用现有的标准试样进行实验,难点也在于无法确定预测试样的受力状况,进而无法得到标准试样的实验条件;齿状等效试样是预测试样和标准试验之间的过度,对齿状等效试样模拟预测式样进行实验,输入为预测试样的实验条件,输出为近似的预测式样的受力状况,对该受力状况进行分析处理可以得到标准试样的实验条件。总结,通过齿状等效试样的过度,实现了利用标准试样对预测试样的疲劳寿命的预测。本发明通过齿状等效试样的实验获得应力集中等效参数、最大主应力、残余应力、粗糙度,保证了预测结果的准确性,通过标准试样的实验获得了齿根的应力状态,保证了实验方法的简便性,做到了准确性与简便性的统一,最终形成了一套浮动摩擦片跨尺度疲劳寿命预测方法及模型,模型在应用上还具有很好的广泛适应性。
同样上述过程可实现对所述浮动摩擦片进行疲劳寿命实验也可实现对所述疲劳寿命预测模型进行验证。
以上实施例仅为本申请的示例性实施例,不用于限制本申请,本申请的保护范围由权利要求书限定。本领域技术人员可以在本申请的实质和保护范围内,对本申请做出各种修改或等同替换,这种修改或等同替换也应视为落在本申请的保护范围内。

Claims (8)

1.一种针对浮动摩擦片跨尺度疲劳寿命预测方法,其特征在于,包括以下步骤:
步骤S1、设置齿状等效试样,将浮动摩擦片服役过程中遭受到的载荷谱加载到齿状等效试样上,通过DIC分析技术分析得到齿状等效试样的齿根应力状态,所述载荷谱由所述浮动摩擦片在服役过程中采集得到,所述齿状等效试样表征与浮动摩擦片具有等效齿状结构的试样,并用于等效获得浮动摩擦片的齿根应力状态;
所述将浮动摩擦片服役过程中遭受到的载荷谱加载到齿状等效试样上,通过DIC分析技术分析得到齿状等效试样的齿根应力状态,包括:
将所述载荷谱作为输入载荷,对齿状等效试样开展疲劳寿命实验,依托摄像机对疲劳寿命实验全过程进行监测;
通过DIC分析技术分析疲劳寿命实验的实验结果得到所述齿状等效试样的齿根应力状态;
步骤S2、将所述齿根应力状态进行分析,并利用参数公式将齿根应力状态转换为应力集中等效参数和最大主应力;
所述将所述齿根应力状态进行分析,并利用参数公式将齿根应力状态转换为应力集中等效参数和最大主应力,包括:
考量残余应力的影响,将所述齿根应力状态依托第一参数公式得到齿根处的等效应力水平、应力集中等效参数和最大主应力,所述第一参数公式为:
Figure DEST_PATH_IMAGE002
式中,σapp1表征最大主应力,Kt为应力集中等效参数,σeff表征等效应力水平,r0表征残余应力影响区深度,a表征齿根近似的椭圆形结构的长轴长,b表征齿根近似的椭圆形结构的短轴长,σ0表征载荷谱中施加在齿根处的外载荷, r表征残余应力影响区深度的积分变量,σres表征残余应力;
在整个疲劳寿命实验的时间段内,依托基于韦伯概率分布的第二参数公式得到基于等效应力水平的最大估计应力,所述第二参数公式为:
Figure DEST_PATH_IMAGE004
式中,σm表征最大估计应力,P表征概率密度;
考虑表面粗糙度影响,通过第三参数公式基于所述最大估计应力得到预测试样的最终应力水平,所述第三参数公式为:
Figure DEST_PATH_IMAGE006
式中,σ预测试样表征预测试样的最终应力水平,
Figure DEST_PATH_IMAGE008
表征第一修正系数,Ra表征粗糙度,所述预测试样表征待预测疲劳寿命的试样;
步骤S3、设置标准试样,基于所述应力集中等效参数和最大主应力,以及所述标准试样的残余应力、粗糙度、应力比为标准试样设计多组疲劳寿命实验,并基于疲劳寿命实验的实验结果得到疲劳寿命预测模型,所述标准试样用于等效获得应力集中等效参数、最大主应力、残余应力和粗糙度对浮动摩擦片疲劳特性的影响。
2.根据权利要求1所述的一种针对浮动摩擦片跨尺度疲劳寿命预测方法,其特征在于:所述基于疲劳寿命实验的实验结果得到疲劳寿命预测模型,包括:
通过改变所述应力集中等效参数和最大主应力,以及结合所述残余应力、粗糙度、应力比的值对标准试样进行多组的疲劳寿命实验;
通过分析多组的疲劳寿命实验的实验结果对所述应力集中等效参数和最大主应力,以及所述残余应力、粗糙度、应力比的参数影响进行评估,基于评估结果反推算出浮动摩擦片的疲劳寿命预测模型,所述疲劳寿命预测模型的函数表达式为:
Figure DEST_PATH_IMAGE010
式中,N标准试样表征标准试样的疲劳寿命,N预测试样表征预测试样的疲劳寿命,σ预测试样表征标准试样的最终应力水平,σ标准试样表征预测试样的最终应力水平,
Figure DEST_PATH_IMAGE012
表征第二修正参数。
3.根据权利要求2所述的一种针对浮动摩擦片跨尺度疲劳寿命预测方法,其特征在于:在利用所述疲劳寿命预测模型对遭受固定载荷谱的多个浮动摩擦片进行疲劳寿命的预测时,仅利用与对应的浮动摩擦片具有等效齿状结构的齿状等效试样进行疲劳寿命实验即可,所述多个浮动摩擦片具有固定形状且具有不同的残余应力和不同的粗糙度,包括:
选取多组具有固定形状且具有不同的残余应力和不同的粗糙度的浮动摩擦片,并为多组所述浮动摩擦片对应设置多组齿状等效试样;
分别对多组齿状等效试样进行疲劳寿命实验分析得到实验结果,以及对选定的标准试样进行疲劳寿命实验分析得到实验结果,并基于粗糙度和残余应力构建出多组浮动摩擦片中表征N标准试样和σ标准试样的非线性映射关系的多组疲劳寿命曲线,所述疲劳寿命曲线的函数表达式标记为
Figure DEST_PATH_IMAGE014
,式中,
Figure DEST_PATH_IMAGE016
表征第i组浮动摩擦片对应的标准试样的最终应力水平,
Figure DEST_PATH_IMAGE018
表征第i组浮动摩擦片对应的标准试样的疲劳寿命,
Figure DEST_PATH_IMAGE020
表征所述非线性映射关系的函数体,
Figure DEST_PATH_IMAGE022
表征第i组浮动摩擦片的粗糙度,
Figure DEST_PATH_IMAGE024
表征第i组浮动摩擦片的残余应力,i为计量常数;
基于疲劳寿命预测模型在疲劳寿命曲线中求得多组齿状等效试样的疲劳寿命,所述多组齿状等效试样的疲劳寿命的计算公式为:
Figure DEST_PATH_IMAGE026
式中,
Figure 641323DEST_PATH_IMAGE018
表征第i组浮动摩擦片对应的齿状等效试样的疲劳寿命,
Figure 89622DEST_PATH_IMAGE016
表征第i组浮动摩擦片对应的齿状等效试样的最终应力水平;
将多组齿状等效试样的疲劳寿命等价作为多组浮动摩擦片的疲劳寿命,以实现对浮动摩擦片的疲劳寿命预测。
4.根据权利要求3所述的一种针对浮动摩擦片跨尺度疲劳寿命预测方法,其特征在于,所述多组齿状等效试样的实验结果和标准试样的实验结果在进行疲劳寿命预测模型运算时需要进行归一化处理。
5.根据权利要求4所述的一种针对浮动摩擦片跨尺度疲劳寿命预测方法,其特征在于:所述对所述浮动摩擦片进行疲劳寿命实验也可实现对所述疲劳寿命预测模型进行验证。
6.一种根据权利要求1-5任一项所述的针对浮动摩擦片跨尺度疲劳寿命预测方法的预测装置,其特征在于,包括:
齿状等效试样实验模块(1),用于对齿状等效试样进行疲劳寿命实验及分析;
标准试样实验模块(2),用于对标准试样进行疲劳寿命实验及分析;
模型构建及验证模块(3),用于构建疲劳寿命预测模型及验证。
7.根据权利要求6所述的预测装置,其特征在于,所述齿状等效试样实验模块、标准试样实验模块和模型构建及验证模块中均集成有通讯模块和数据存储模块,所述通讯模块用于在所述齿状等效试样实验模块、标准试样实验模块和模型构建及验证模块中进行数据交互,所述数据存储模块用于在所述齿状等效试样实验模块、标准试样实验模块和模型构建及验证模块中进行数据存储。
8.根据权利要求7所述的预测装置,其特征在于,所述齿状等效试样与浮动摩擦片采用同种制作材料。
CN202210068797.3A 2022-01-21 2022-01-21 一种针对浮动摩擦片跨尺度疲劳寿命预测方法及装置 Active CN114117826B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210068797.3A CN114117826B (zh) 2022-01-21 2022-01-21 一种针对浮动摩擦片跨尺度疲劳寿命预测方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210068797.3A CN114117826B (zh) 2022-01-21 2022-01-21 一种针对浮动摩擦片跨尺度疲劳寿命预测方法及装置

Publications (2)

Publication Number Publication Date
CN114117826A CN114117826A (zh) 2022-03-01
CN114117826B true CN114117826B (zh) 2022-04-29

Family

ID=80360936

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210068797.3A Active CN114117826B (zh) 2022-01-21 2022-01-21 一种针对浮动摩擦片跨尺度疲劳寿命预测方法及装置

Country Status (1)

Country Link
CN (1) CN114117826B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104634676A (zh) * 2015-02-13 2015-05-20 中国北方车辆研究所 一种控制摩擦片动态试验载荷的方法及实现该方法的系统
CN111289555A (zh) * 2020-03-27 2020-06-16 中国北方车辆研究所 一种摩擦片热疲劳载荷全场测试方法
CN111931393A (zh) * 2019-04-28 2020-11-13 沈阳工业大学 一种废旧机床关键零件剩余寿命预测评估方法
CN112504594A (zh) * 2020-12-29 2021-03-16 中国北方车辆研究所 浮动摩擦片冲击和振动疲劳耦合加载装置
CN113268939A (zh) * 2021-06-07 2021-08-17 四川大学 基于瞬态混合润滑理论的谐波减速接触疲劳寿命预测方法
WO2021227925A1 (zh) * 2020-05-09 2021-11-18 清华大学 基于主载荷模式的缺口件疲劳寿命预测方法及预测装置

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8297130B2 (en) * 2008-11-12 2012-10-30 The United States Of America As Represented By The Secretary Of The Air Force Microtesting rig with variable compliance loading fibers for measuring mechanical properties of small specimens
CN109522569B (zh) * 2017-09-19 2022-10-04 浙江大学 基于Weibull方程和残余变形的混凝土疲劳寿命预测方法和装置
CN112926163B (zh) * 2021-04-07 2022-04-29 唐山学院 基于Workbench和Matlab构建螺栓疲劳寿命预测模型的方法及装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104634676A (zh) * 2015-02-13 2015-05-20 中国北方车辆研究所 一种控制摩擦片动态试验载荷的方法及实现该方法的系统
CN111931393A (zh) * 2019-04-28 2020-11-13 沈阳工业大学 一种废旧机床关键零件剩余寿命预测评估方法
CN111289555A (zh) * 2020-03-27 2020-06-16 中国北方车辆研究所 一种摩擦片热疲劳载荷全场测试方法
WO2021227925A1 (zh) * 2020-05-09 2021-11-18 清华大学 基于主载荷模式的缺口件疲劳寿命预测方法及预测装置
CN112504594A (zh) * 2020-12-29 2021-03-16 中国北方车辆研究所 浮动摩擦片冲击和振动疲劳耦合加载装置
CN113268939A (zh) * 2021-06-07 2021-08-17 四川大学 基于瞬态混合润滑理论的谐波减速接触疲劳寿命预测方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Establishment and Verification of Multiaxis Fatigue Life Prediction Model;Zhuo Fu 等;《Hindawi》;20211231;第1-6页 *
Impacting load control of floating supported friction plate and its experimental verification;Keyan Ning 等;《Journal of Physics Conference Series》;20171231;第1-8页 *
摩擦片非线性损伤累积计算与寿命预测(英文);王玉 等;《机床与液压》;20170928;第45卷(第18期);第23-26页 *
磨损加疲劳载荷下的协同疲劳行为;梁益龙 等;《摩擦学学报》;20150115;第35卷(第1期);第59-65页 *
高能摩擦副齿部啮合非线性冲击损伤新算法及疲劳寿命预测;李杰 等;《中国公路学报》;20200408;第33卷(第8期);第110-118页 *

Also Published As

Publication number Publication date
CN114117826A (zh) 2022-03-01

Similar Documents

Publication Publication Date Title
Gu et al. A new method of accelerated life testing based on the Grey System Theory for a model-based lithium-ion battery life evaluation system
Danciu et al. The 2014 Earthquake Model of the Middle East: ground motion model and uncertainties
CN109685289B (zh) 高炉炉况顺行预测方法、装置及系统
Goulet et al. Data‐driven post‐earthquake rapid structural safety assessment
CN107436963B (zh) 一种基于Copula函数多元退化的O型橡胶密封圈寿命预测方法
Tovo On the fatigue reliability evaluation of structural components under service loading
Chen et al. A data heterogeneity modeling and quantification approach for field pre-assessment of chloride-induced corrosion in aging infrastructures
Geyskens et al. Bayesian prediction of elastic modulus of concrete
Lipke Project duration forecasting
Zhang A fuzzy residual strength based fatigue life prediction method
Kassem et al. Quantification of the inherent uncertainty in the relaxation modulus and creep compliance of asphalt mixes
CN114117826B (zh) 一种针对浮动摩擦片跨尺度疲劳寿命预测方法及装置
Messner et al. Initial Development of an Improved Creep-Fatigue Design Method That Avoids the Separate Evaluation of Creep and Fatigue Damage and Eliminates the Requirement for Stress Classification
Hinrichs The impact of measurement uncertainty on the producer’s and user’s risks, on classification and conformity assessment: an example based on tests on some construction products
Aflatooni et al. Synthetic rating system for railway bridge management
Dominguez et al. A statistical model for fatigue crack growth under random loads including retardation effects
US20090089005A1 (en) Generation of test specifications based on measured data points
CN110895624B (zh) 基于最大熵谱估计的加速贮存与自然贮存退化数据一致性检验法
CN114741284B (zh) 任务可靠性评估方法、装置、计算机设备和存储介质
CN115640542A (zh) 基于贝叶斯理论的智能电能表状态评估方法及评估装置
CN111881259B (zh) 一种基于文本挖掘的设备故障概率评估方法及系统
CN110210711B (zh) 一种电动汽车电能计量准确度评估方法
CN113742248A (zh) 一种基于项目测量数据进行组织过程预测的方法及系统
Ye et al. Conformity control of chloride diffusivity for reinforced concrete structures of Hong Kong–Zhuhai–Macau sea-link project in marine environment
CN105528657A (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