CN110686915B - 多应力加速试验剖面确定方法、系统、介质及设备 - Google Patents

多应力加速试验剖面确定方法、系统、介质及设备 Download PDF

Info

Publication number
CN110686915B
CN110686915B CN201911018342.5A CN201911018342A CN110686915B CN 110686915 B CN110686915 B CN 110686915B CN 201911018342 A CN201911018342 A CN 201911018342A CN 110686915 B CN110686915 B CN 110686915B
Authority
CN
China
Prior art keywords
stress
test
acceleration
temperature
information
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
CN201911018342.5A
Other languages
English (en)
Other versions
CN110686915A (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.)
Shanghai Space Precision Machinery Research Institute
Original Assignee
Shanghai Space Precision Machinery Research Institute
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 Shanghai Space Precision Machinery Research Institute filed Critical Shanghai Space Precision Machinery Research Institute
Priority to CN201911018342.5A priority Critical patent/CN110686915B/zh
Publication of CN110686915A publication Critical patent/CN110686915A/zh
Application granted granted Critical
Publication of CN110686915B publication Critical patent/CN110686915B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M99/00Subject matter not provided for in other groups of this subclass
    • G01M99/005Testing of complete machines, e.g. washing-machines or mobile phones

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明提供了一种多应力加速试验剖面确定方法、系统、介质及设备,包括:步骤1:确定试验应力类型;步骤2:根据试验应力类型信息,建立温湿电多应力加速模型;步骤3:建立试验对象有限元分析模型;步骤4:确定初始加速试验剖面,选取初始加速试验载荷量级参数、初始加速试验载荷周期参数;步骤5:创建温湿电应力综合加速因子模型;步骤6:计算试验对象的疲劳寿命,再将疲劳寿命除以一个周期时长即可得到破坏循环数;步骤7:计算试验对象在不同应力量级下的破坏循环数;步骤8:计算试验对象的加速比;步骤9:判断试验对象加速比是否满足预定试验要求;步骤10:获取多应力加速试验指导信息;本发明能够降低试验样件的大量消耗。

Description

多应力加速试验剖面确定方法、系统、介质及设备
技术领域
本发明涉及试验剖面确定方法领域,具体地,涉及一种多应力加速试验剖面确定方法、系统、介质及设备,尤其涉及一种多应力条件下加速试验载荷量级与周期确定方法。
背景技术
由于机电装备加速寿命试验有助于建立产品和工艺设计的质量而导致更好的产品设计和生产过程,进一步提高平均无故障时间,大大节约成本。传统加速寿命试验方法往往是针对单一应力条件的。但从机电类产品或试验件的失效模式分析中可以发现,温度应力、湿度应力和电应力的综合作用是导致其破坏失效的主要因素。因此,机电类产品或试验件在多应力加速试验中需要考虑多种环境和电应力载荷,即需要分析多应力条件下的加速寿命试验技术。由于试验剖面中的载荷量级以及作用周期的综合作用是复杂的,仿真分析可以辅助技术员快速获得试验剖面设计所需的各项关键参数。
专利文献CN107390668B公开了一种确定防滑刹车控制装置加严可靠性试验剖面的方法,用相同数字电路防滑刹车控制装置高加速寿命试验的数据确定试验应力,形成加严的可靠性试验剖面。该发明中的加严试验剖面中的试验应力不仅包含环境试验的极值条件,还包含大于环境极值条件的情况,按照该发明的加严试验剖面进行可靠性试验,确定和纠正该防滑刹车控制装置的故障隐患。该专利并不能很好地适用于多应力加速试验剖面确定中。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种多应力加速试验剖面确定方法、系统、介质及设备。
根据本发明提供的一种多应力加速试验剖面确定方法,包括:步骤1:确定试验应力类型,获取试验应力类型信息;对于机电类试验对象的分析需包括温度应力、湿度应力和电应力;步骤2:根据试验应力类型信息,建立温湿电多应力加速模型,以破坏循环数表征,获取试验对象产品特征寿命信息;步骤3:建立试验对象有限元分析模型,得到不同应力类型的不同载荷量级下的应力应变响应结果,确定试验对象的工作应力极限,获取试验对象工作应力极限信息;步骤4:确定初始加速试验剖面,选取初始加速试验载荷量级参数、初始加速试验载荷周期参数,获取初始加速试验载荷量级信息、初始加速试验载荷周期参数信息;初始加速试验载荷量级的选取为试验对象的工作应力极限,初始加速试验载荷周期参数依据行业标准的剖面设计一般方法给出。步骤5:创建温湿电应力综合加速因子模型,获取温湿电应力综合加速因子信息;步骤6:利用修正的Coffin-Manson方程计算试验对象的疲劳寿命,再将疲劳寿命除以一个周期时长即可得到破坏循环数,获取破坏循环数信息;步骤7:根据破坏循环数信息,计算试验对象在不同应力量级下的破坏循环数,通过最小二乘法确定温湿电应力综合加速因子模型中待定系数a、b、c、d、e的值,得到基于寿命数据针对试验对象的温湿电多应力条件下的综合加速模型,获取湿电多应力条件下综合加速模型信息;步骤8:根据湿电多应力条件下综合加速模型信息,计算试验对象的加速比,获取试验对象加速比信息;步骤9:根据试验对象加速比信息,判断试验对象加速比是否满足预定试验要求,若不满足调整试验对象加速试验剖面参数,继续执行步骤8再次计算,若满足则确定最终的多应力加速试验剖面,获取多应力加速试验剖面确定信息;步骤10:根据多应力加速试验剖面确定信息,获取多应力加速试验指导信息;试验对象加速试验剖面参数包括加速试验载荷量级参数、加速试验周期参数。
优选地,所述步骤2包括:步骤2.1:采用温湿电多应力加速模型计算试验对象的产品特征寿命,其中,试验对象的产品特征寿命计算公式如下:
Figure BDA0002246404550000021
其中,η为试验对象的产品特征寿命(以破坏循环数表征),I为电流(单位:A),H为相对湿度(单位:%),f为循环频率(单位:Hz),Tmax为最高温度(单位:℃),Tmin为最低温度(单位:℃),Ea为激活能取1eV,k为玻尔兹曼常数(值为8.617×10-5 eV/K),a、b、c、d、e为待定系数。
优选地,所述步骤5包括:步骤5.1:计算温湿电应力综合加速因子,其中,温湿电应力综合加速因子的计算公式如下:
Figure BDA0002246404550000022
其中,AF为温湿电应力综合加速因子,Is、Hs、fs、Tsmax、Tsmin为加速应力条件参数,I0、H0、f0、T0max、T0min为正常工作应力条件,η0为试验对象的产品特征寿命,I为电流(单位:A),H为相对湿度(单位:%),f为循环频率(单位:Hz),Tmax为最高温度(单位:℃),Tmin为最低温度(单位:℃),Ea为激活能(单位:eV),k为玻尔兹曼常数(值为8.617×10-5eV/K),a、b、c、d、e为待定系数。
优选地,步骤6包括:步骤6.1:利用修正的Coffin-Manson方程计算试验对象的疲劳寿命,所述修正的Coffin-Manson方程如下:
Figure BDA0002246404550000031
Figure BDA0002246404550000032
β=-0.442-6×10-4tm+1.74×10-2ln(1+nf)
其中,Nf为材料的疲劳寿命,Δγp为等效塑性剪切应变范围,εpmax为等效塑性应变最大值,εpmin为等效塑性应变最小值,εf为疲劳韧性系数,β为疲劳韧性指数,tm为平均温度,nf为材料的循环周期数。
根据本发明提供的一种多应力加速试验剖面确定系统,包括:模块1:确定试验应力类型,获取试验应力类型信息;对于机电类试验对象的分析需包括温度应力、湿度应力和电应力;模块2:根据试验应力类型信息,建立温湿电多应力加速模型,以破坏循环数表征,获取试验对象产品特征寿命信息;模块3:建立试验对象有限元分析模型,得到不同应力类型的不同载荷量级下的应力应变响应结果,确定试验对象的工作应力极限,获取试验对象工作应力极限信息;模块4:确定初始加速试验剖面,选取初始加速试验载荷量级参数、初始加速试验载荷周期参数,获取初始加速试验载荷量级信息、初始加速试验载荷周期参数信息;初始加速试验载荷量级的选取为试验对象的工作应力极限,初始加速试验载荷周期参数依据行业标准的剖面设计一般方法给出。模块5:创建温湿电应力综合加速因子模型,获取温湿电应力综合加速因子信息;模块6:利用修正的Coffin-Manson方程计算试验对象的疲劳寿命,再将疲劳寿命除以一个周期时长即可得到破坏循环数,获取破坏循环数信息;模块7:根据破坏循环数信息,计算试验对象在不同应力量级下的破坏循环数,通过最小二乘法确定温湿电应力综合加速因子模型中待定系数a、b、c、d、e的值。得到基于寿命数据针对试验对象的温湿电多应力条件下的综合加速模型,获取湿电多应力条件下综合加速模型信息;模块8:根据湿电多应力条件下综合加速模型信息,计算试验对象的加速比,获取试验对象加速比信息;模块9:根据试验对象加速比信息,判断试验对象加速比是否满足预定试验要求,若不满足调整试验对象加速试验剖面参数,继续执行模块8再次计算,若满足则确定最终的多应力加速试验剖面,获取多应力加速试验剖面确定信息;模块10:根据多应力加速试验剖面确定信息,获取多应力加速试验指导信息;试验对象加速试验剖面参数包括加速试验载荷量级参数、加速试验周期参数。
优选地,所述模块2包括:模块2.1:采用温湿电多应力加速模型计算试验对象的产品特征寿命,其中,试验对象的产品特征寿命计算公式如下:
Figure BDA0002246404550000041
其中,η为试验对象的产品特征寿命(以破坏循环数表征),I为电流(单位:A),H为相对湿度(单位:%),f为循环频率(单位:Hz),Tmax为最高温度(单位:℃),Tmin为最低温度(单位:℃),Ea为激活能取1eV,k为玻尔兹曼常数(值为8.617×10-5 eV/K),a、b、c、d、e为待定系数。
优选地,所述模块5包括:模块5.1:计算温湿电应力综合加速因子,其中,温湿电应力综合加速因子的计算公式如下:
Figure BDA0002246404550000042
其中,AF为温湿电应力综合加速因子,Is、Hs、fs、Tsmax、Tsmin为加速应力条件参数,I0、H0、f0、T0max、T0min为正常工作应力条件,η0为试验对象的产品特征寿命,I为电流(单位:A),H为相对湿度(单位:%),f为循环频率(单位:Hz),Tmax为最高温度(单位:℃),Tmin为最低温度(单位:℃),Ea为激活能(单位:eV),k为玻尔兹曼常数(值为8.617×10-5eV/K),a、b、c、d、e为待定系数。
优选地,模块6包括:模块6.1:利用修正的Coffin-Manson方程计算试验对象的疲劳寿命,所述修正的Coffin-Manson方程如下:
Figure BDA0002246404550000043
Figure BDA0002246404550000044
β=-0.442-6×10-4tm+1.74×10-2ln(1+nf)
其中,Nf为材料的疲劳寿命,Δγp为等效塑性剪切应变范围,εpmax为等效塑性应变最大值,εpmin为等效塑性应变最小值,εf为疲劳韧性系数,β为疲劳韧性指数,tm为平均温度,nf为材料的循环周期数。
根据本发明提供的一种存储有计算机程序的计算机可读存储介质,所述计算机程序被处理器执行时实现多应力加速试验剖面确定方法的步骤。
根据本发明提供的一种多应力加速试验剖面确定设备,包括:控制器;
所述控制器包括所述的存储有计算机程序的计算机可读存储介质,所述计算机程序被处理器执行时实现多应力加速试验剖面确定方法的步骤;或者,所述控制器包括多应力加速试验剖面确定系统。
与现有技术相比,本发明具有如下的有益效果:
1、本发明方法可以改变以多次摸底试验经逐步试探建立加速试验剖面的工作模式,降低试验样件的大量消耗;
2、本发明方法运用于多个运载与武器型号的加速试验剖面设计中,大大缩短试验周期;
3、本发明可以科学定量的进行加速试验剖面设计,大大缩短试验周期,有效减少终的多应力加速试验剖面,为后续多应力加速试验的开展提供指导试验对象消耗。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为本发明实施例中的多应力条件下加速试验载荷量级与周期确定方法流程示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。
如图1所示,根据本发明提供的一种多应力加速试验剖面确定方法,包括:步骤1:确定试验应力类型,获取试验应力类型信息;对于机电类试验对象的分析需包括温度应力、湿度应力和电应力;步骤2:根据试验应力类型信息,建立温湿电多应力加速模型,以破坏循环数表征,获取试验对象产品特征寿命信息;步骤3:建立试验对象有限元分析模型,得到不同应力类型的不同载荷量级下的应力应变响应结果,确定试验对象的工作应力极限,获取试验对象工作应力极限信息;步骤4:确定初始加速试验剖面,选取初始加速试验载荷量级参数、初始加速试验载荷周期参数,获取初始加速试验载荷量级信息、初始加速试验载荷周期参数信息;初始加速试验载荷量级的选取为试验对象的工作应力极限,初始加速试验载荷周期参数依据行业标准的剖面设计一般方法给出。步骤5:创建温湿电应力综合加速因子模型,获取温湿电应力综合加速因子信息;步骤6:利用修正的Coffin-Manson方程计算试验对象的疲劳寿命,再将疲劳寿命除以一个周期时长即可得到破坏循环数,获取破坏循环数信息;步骤7:根据破坏循环数信息,计算试验对象在不同应力量级下的破坏循环数,通过最小二乘法确定温湿电应力综合加速因子模型中待定系数a、b、c、d、e的值,得到基于寿命数据针对试验对象的温湿电多应力条件下的综合加速模型,获取湿电多应力条件下综合加速模型信息;步骤8:根据湿电多应力条件下综合加速模型信息,计算试验对象的加速比,获取试验对象加速比信息;步骤9:根据试验对象加速比信息,判断试验对象加速比是否满足预定试验要求,若不满足调整试验对象加速试验剖面参数,继续执行步骤8再次计算,若满足则确定最终的多应力加速试验剖面,获取多应力加速试验剖面确定信息;步骤10:根据多应力加速试验剖面确定信息,获取多应力加速试验指导信息;试验对象加速试验剖面参数包括加速试验载荷量级参数、加速试验周期参数。
优选地,所述步骤2包括:步骤2.1:采用温湿电多应力加速模型计算试验对象的产品特征寿命,其中,试验对象的产品特征寿命计算公式如下:
Figure BDA0002246404550000061
其中,η为试验对象的产品特征寿命(以破坏循环数表征),I为电流(单位:A),H为相对湿度(单位:%),f为循环频率(单位:Hz),Tmax为最高温度(单位:℃),Tmin为最低温度(单位:℃),Ea为激活能取1eV,k为玻尔兹曼常数(值为8.617×10-5 eV/K),a、b、c、d、e为待定系数。
优选地,所述步骤5包括:步骤5.1:计算温湿电应力综合加速因子,其中,温湿电应力综合加速因子的计算公式如下:
Figure BDA0002246404550000062
其中,AF为温湿电应力综合加速因子,Is、Hs、fs、Tsmax、Tsmin为加速应力条件参数,I0、H0、f0、T0max、T0min为正常工作应力条件,η0为试验对象的产品特征寿命,I为电流(单位:A),H为相对湿度(单位:%),f为循环频率(单位:Hz),Tmax为最高温度(单位:℃),Tmin为最低温度(单位:℃),Ea为激活能(单位:eV),k为玻尔兹曼常数(值为8.617×10-5eV/K),a、b、c、d、e为待定系数。
优选地,步骤6包括:步骤6.1:利用修正的Coffin-Manson方程计算试验对象的疲劳寿命,所述修正的Coffin-Manson方程如下:
Figure BDA0002246404550000071
Figure BDA0002246404550000072
β=-0.442-6×10-4tm+1.74×10-2ln(1+nf)
其中,Nf为材料的疲劳寿命,Δγp为等效塑性剪切应变范围,εpmax为等效塑性应变最大值,εpmin为等效塑性应变最小值,εf为疲劳韧性系数,β为疲劳韧性指数,tm为平均温度,nf为材料的循环周期数。
具体地,在一个实施例中,一种多应力条件下加速试验载荷量级与周期确定方法,包括以下步骤:
步骤1:确定试验应力类型,对于机电类试验对象的分析需包括温度应力、湿度应力和电应力。
步骤2:建立温-湿-电多应力加速模型:
Figure BDA0002246404550000073
其中,η为试验对象的产品特征寿命(以破坏循环数表征),I为电流(单位:A),H为相对湿度(单位:%),f为循环频率(单位:Hz),Tmax为最高温度(单位:℃),Tmin为最低温度(单位:℃),Ea为激活能(单位:eV),总结国内外最新研究成果,结合项目仿真分析,此处取1,k为玻尔兹曼常数(值为8.617×10-5eV/K),a、b、c、d、e为待定系数。
步骤3:建立试验对象有限元分析模型,基于仿真分析得到不同应力类型的不同载荷量级下的应力应变响应结果,进一步确定试验对象的工作应力极限。
步骤4:确定初始加速试验剖面。加速试验剖面是由在两个极限应力之间的温度、湿度和电应力等环境应力综合作用的循环周期构成的,其中温度循环的参数包括温度极限、温度变化率和循环数,湿度参数包括应力量级和驻留时间,电应力参数包括应力量级和驻留时间。初始加速试验载荷量级的选取为试验对象的工作应力极限,初始加速试验载荷周期参数依据行业标准的剖面设计一般方法给出。
步骤5:温-湿-电应力综合加速因子模型:
Figure BDA0002246404550000081
其中,Is、Hs、fs、Tsmax、Tsmin为加速应力条件,I0、H0、f0、T0max、T0min为正常工作应力条件,η0为试验对象的产品特征寿命,I为电流(单位:A),H为相对湿度(单位:%),f为循环频率(单位:Hz),Tmax为最高温度(单位:℃),Tmin为最低温度(单位:℃),Ea为激活能(单位:eV),k为玻尔兹曼常数(值为8.617×10-5eV/K),a、b、c、d、e为待定系数。
步骤6:基于步骤3的温-湿-电应力仿真分析结果,利用修正的Coffin-Manson方程计算试验对象的疲劳寿命,再将疲劳寿命除以一个周期时长即可得到破坏循环数。修正的Coffin-Manson方程如下:
Figure BDA0002246404550000082
Figure BDA0002246404550000083
β=-0.442-6×10-4tm+1.74×10-2ln(1+nf)
其中,Nf为材料的疲劳寿命,Δγp为等效塑性剪切应变范围,εpmax为等效塑性应变最大值,εpmin为等效塑性应变最小值,εf为疲劳韧性系数,β为疲劳韧性指数,tm为平均温度,nf为材料的循环周期数。
步骤7:利用步骤6)计算的试验对象在不同应力量级下的破坏循环数,通过最小二乘法确定温-湿-电应力综合加速因子模型中待定系数a、b、c、d、e的值。得到基于寿命数据针对试验对象的温-湿-电多应力条件下的综合加速模型。
步骤8:利用步骤7)中针对试验对象的温-湿-电多应力条件下的综合加速模型,计算试验对象的加速比。
步骤9:判断加速比是否满足试验要求,若不满足调整试验对象加速试验剖面参数,包括加速试验载荷量级与周期参数,返回步骤8)再次计算,若满足则确定最终的多应力加速试验剖面。
本领域技术人员可以将本发明提供的多应力加速试验剖面确定方法,理解为本发明提供的多应力加速试验剖面确定系统的一个实施例。即,所述多应力加速试验剖面确定系统可以通过执行所述多应力加速试验剖面确定方法的步骤流程实现。
根据本发明提供的一种多应力加速试验剖面确定系统,包括:模块1:确定试验应力类型,获取试验应力类型信息;对于机电类试验对象的分析需包括温度应力、湿度应力和电应力;模块2:根据试验应力类型信息,建立温湿电多应力加速模型,以破坏循环数表征,获取试验对象产品特征寿命信息;模块3:建立试验对象有限元分析模型,得到不同应力类型的不同载荷量级下的应力应变响应结果,确定试验对象的工作应力极限,获取试验对象工作应力极限信息;模块4:确定初始加速试验剖面,选取初始加速试验载荷量级参数、初始加速试验载荷周期参数,获取初始加速试验载荷量级信息、初始加速试验载荷周期参数信息;初始加速试验载荷量级的选取为试验对象的工作应力极限,初始加速试验载荷周期参数依据行业标准的剖面设计一般方法给出。模块5:创建温湿电应力综合加速因子模型,获取温湿电应力综合加速因子信息;模块6:利用修正的Coffin-Manson方程计算试验对象的疲劳寿命,再将疲劳寿命除以一个周期时长即可得到破坏循环数,获取破坏循环数信息;模块7:根据破坏循环数信息,计算试验对象在不同应力量级下的破坏循环数,通过最小二乘法确定温湿电应力综合加速因子模型中待定系数a、b、c、d、e的值,得到基于寿命数据针对试验对象的温湿电多应力条件下的综合加速模型,获取湿电多应力条件下综合加速模型信息;模块8:根据湿电多应力条件下综合加速模型信息,计算试验对象的加速比,获取试验对象加速比信息;模块9:根据试验对象加速比信息,判断试验对象加速比是否满足预定试验要求,若不满足调整试验对象加速试验剖面参数,继续执行模块8再次计算,若满足则确定最终的多应力加速试验剖面,获取多应力加速试验剖面确定信息;模块10:根据多应力加速试验剖面确定信息,获取多应力加速试验指导信息;试验对象加速试验剖面参数包括加速试验载荷量级参数、加速试验周期参数。
优选地,所述模块2包括:模块2.1:采用温湿电多应力加速模型计算试验对象的产品特征寿命,其中,试验对象的产品特征寿命计算公式如下:
Figure BDA0002246404550000091
其中,η为试验对象的产品特征寿命(以破坏循环数表征),I为电流(单位:A),H为相对湿度(单位:%),f为循环频率(单位:Hz),Tmax为最高温度(单位:℃),Tmin为最低温度(单位:℃),Ea为激活能取1eV,k为玻尔兹曼常数(值为8.617×10-5 eV/K),a、b、c、d、e为待定系数。
优选地,所述模块5包括:模块5.1:计算温湿电应力综合加速因子,其中,温湿电应力综合加速因子的计算公式如下:
Figure BDA0002246404550000092
其中,AF为温湿电应力综合加速因子,Is、Hs、fs、Tsmax、Tsmin为加速应力条件参数,I0、H0、f0、T0max、T0min为正常工作应力条件,η0为试验对象的产品特征寿命,I为电流(单位:A),H为相对湿度(单位:%),f为循环频率(单位:Hz),Tmax为最高温度(单位:℃),Tmin为最低温度(单位:℃),Ea为激活能(单位:eV),k为玻尔兹曼常数(值为8.617×10-5eV/K),a、b、c、d、e为待定系数。
优选地,模块6包括:模块6.1:利用修正的Coffin-Manson方程计算试验对象的疲劳寿命,所述修正的Coffin-Manson方程如下:
Figure BDA0002246404550000101
Figure BDA0002246404550000102
β=-0.442-6×10-4tm+1.74×10-2ln(1+nf)
其中,Nf为材料的疲劳寿命,Δγp为等效塑性剪切应变范围,εpmax为等效塑性应变最大值,εpmin为等效塑性应变最小值,εf为疲劳韧性系数,β为疲劳韧性指数,tm为平均温度,nf为材料的循环周期数。
根据本发明提供的一种存储有计算机程序的计算机可读存储介质,所述计算机程序被处理器执行时实现多应力加速试验剖面确定方法的步骤。
根据本发明提供的一种多应力加速试验剖面确定设备,包括:控制器;
所述控制器包括所述的存储有计算机程序的计算机可读存储介质,所述计算机程序被处理器执行时实现多应力加速试验剖面确定方法的步骤;或者,所述控制器包括多应力加速试验剖面确定系统。
本发明方法可以改变以多次摸底试验经逐步试探建立加速试验剖面的工作模式,降低试验样件的大量消耗;本发明方法运用于多个运载与武器型号的加速试验剖面设计中,大大缩短试验周期;本发明可以科学定量的进行加速试验剖面设计,大大缩短试验周期,有效减少终的多应力加速试验剖面,为后续多应力加速试验的开展提供指导试验对象消耗。
本领域技术人员知道,除了以纯计算机可读程序代码方式实现本发明提供的系统及其各个装置、单元、单元以外,完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统及其各个装置、单元、单元以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同功能。所以,本发明提供的系统及其各项装置、单元、单元可以被认为是一种硬件部件,而对其内包括的用于实现各种功能的装置、单元、单元也可以视为硬件部件内的结构;也可以将用于实现各种功能的装置、单元、单元视为既可以是实现方法的软件单元又可以是硬件部件内的结构。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。

Claims (6)

1.一种多应力加速试验剖面确定方法,其特征在于,包括:
步骤1:确定试验应力类型,获取试验应力类型信息;
步骤2:根据试验应力类型信息,建立温湿电多应力加速模型,以破坏循环数表征,获取试验对象产品特征寿命信息;
步骤3:建立试验对象有限元分析模型,得到不同应力类型的不同载荷量级下的应力应变响应结果,确定试验对象的工作应力极限,获取试验对象工作应力极限信息;
步骤4:确定初始加速试验剖面,选取初始加速试验载荷量级参数、初始加速试验载荷周期参数,获取初始加速试验载荷量级信息、初始加速试验载荷周期参数信息;
步骤5:创建温湿电应力综合加速因子模型,获取温湿电应力综合加速因子信息;
步骤6:计算试验对象的疲劳寿命,再将疲劳寿命除以一个周期时长即可得到破坏循环数,获取破坏循环数信息;
步骤7:根据破坏循环数信息,计算试验对象在不同应力量级下的破坏循环数,确定温湿电应力综合加速因子模型中待定系数值,得到基于寿命数据针对试验对象的温湿电多应力条件下的综合加速模型,获取温湿电多应力条件下综合加速模型信息;
步骤8:根据温湿电多应力条件下综合加速模型信息,计算试验对象的加速比,获取试验对象加速比信息;
步骤9:根据试验对象加速比信息,判断试验对象加速比是否满足预定试验要求,若不满足调整试验对象加速试验剖面参数,继续执行步骤8再次计算,若满足则确定最终的多应力加速试验剖面,获取多应力加速试验剖面确定信息;
步骤10:根据多应力加速试验剖面确定信息,获取多应力加速试验指导信息;
试验对象加速试验剖面参数包括加速试验载荷量级参数、加速试验周期参数;
所述步骤2包括:
步骤2.1:采用温湿电多应力加速模型计算试验对象的产品特征寿命,其中,试验对象的产品特征寿命计算公式如下
Figure FDA0003004258030000011
其中,η为试验对象的产品特征寿命,I为电流,H为相对湿度,f为循环频率,Tmax为最高温度,Tmin为最低温度,Ea为激活能量,Ea的单位取1eV,k为玻尔兹曼常数,a、b、c、d、e为待定系数;
所述步骤5包括:
步骤5.1:计算温湿电应力综合加速因子,其中,温湿电应力综合加速因子的计算公式如下:
Figure FDA0003004258030000021
其中,AF为温湿电应力综合加速因子,Is、Hs、fs、Tsmax、Tsmin为加速应力条件参数,I0、H0、f0、T0max、T0min为正常工作应力条件,η0为试验对象的产品特征寿命,I为电流,H为相对湿度,f为循环频率,Tmax为最高温度,Tmin为最低温度,Ea为激活能,k为玻尔兹曼常数,a、b、c、d、e为待定系数。
2.根据权利要求1所述的多应力加速试验剖面确定方法,其特征在于,步骤6包括:
步骤6.1:利用修正的Coffin-Manson方程计算试验对象的疲劳寿命,所述修正的Coffin-Manson方程如下:
Figure FDA0003004258030000022
Figure FDA0003004258030000023
β=-0.442-6×10-4tm+1.74×10-2ln(1+nf)
其中,Nf为材料的疲劳寿命,Δγp为等效塑性剪切应变范围,εpmax为等效塑性应变最大值,εpmin为等效塑性应变最小值,εf为疲劳韧性系数,β为疲劳韧性指数,tm为平均温度,nf为材料的循环周期数。
3.一种多应力加速试验剖面确定系统,其特征在于,包括:
模块1:确定试验应力类型,获取试验应力类型信息;
模块2:根据试验应力类型信息,建立温湿电多应力加速模型,以破坏循环数表征,获取试验对象产品特征寿命信息;
模块3:建立试验对象有限元分析模型,得到不同应力类型的不同载荷量级下的应力应变响应结果,确定试验对象的工作应力极限,获取试验对象工作应力极限信息;
模块4:确定初始加速试验剖面,选取初始加速试验载荷量级参数、初始加速试验载荷周期参数,获取初始加速试验载荷量级信息、初始加速试验载荷周期参数信息;
模块5:创建温湿电应力综合加速因子模型,获取温湿电应力综合加速因子信息;
模块6:计算试验对象的疲劳寿命,再将疲劳寿命除以一个周期时长即可得到破坏循环数,获取破坏循环数信息;
模块7:根据破坏循环数信息,计算试验对象在不同应力量级下的破坏循环数,确定温湿电应力综合加速因子模型中待定系数值,得到基于寿命数据针对试验对象的温湿电多应力条件下的综合加速模型,获取温湿电多应力条件下综合加速模型信息;
模块8:根据温湿电多应力条件下综合加速模型信息,计算试验对象的加速比,获取试验对象加速比信息;
模块9:根据试验对象加速比信息,判断试验对象加速比是否满足预定试验要求,若不满足调整试验对象加速试验剖面参数,继续执行模块8再次计算,若满足则确定最终的多应力加速试验剖面,获取多应力加速试验剖面确定信息;
模块10:根据多应力加速试验剖面确定信息,获取多应力加速试验指导信息;
试验对象加速试验剖面参数包括加速试验载荷量级参数、加速试验周期参数;
所述模块2包括:
模块2.1:采用温湿电多应力加速模型计算试验对象的产品特征寿命,其中,试验对象的产品特征寿命计算公式如下
Figure FDA0003004258030000031
其中,η为试验对象的产品特征寿命,I为电流,H为相对湿度,f为循环频率,Tmax为最高温度,Tmin为最低温度,Ea为激活能量,Ea的单位取1eV,k为玻尔兹曼常数,a、b、c、d、e为待定系数;
所述模块5包括:
模块5.1:计算温湿电应力综合加速因子,其中,温湿电应力综合加速因子的计算公式如下:
Figure FDA0003004258030000032
其中,AF为温湿电应力综合加速因子,Is、Hs、fs、Tsmax、Tsmin为加速应力条件参数,I0、H0、f0、T0max、T0min为正常工作应力条件,η0为试验对象的产品特征寿命,I为电流,H为相对湿度,f为循环频率,Tmax为最高温度,Tmin为最低温度,Ea为激活能,k为玻尔兹曼常数,a、b、c、d、e为待定系数。
4.根据权利要求3所述的多应力加速试验剖面确定系统,其特征在于,模块6包括:
模块6.1:利用修正的Coffin-Manson方程计算试验对象的疲劳寿命,所述修正的Coffin-Manson方程如下:
Figure FDA0003004258030000041
Figure FDA0003004258030000042
β=-0.442-6×10-4tm+1.74×10-2ln(1+nf)
其中,Nf为材料的疲劳寿命,Δγp为等效塑性剪切应变范围,εpmax为等效塑性应变最大值,εpmin为等效塑性应变最小值,εf为疲劳韧性系数,β为疲劳韧性指数,tm为平均温度,nf为材料的循环周期数。
5.一种存储有计算机程序的计算机可读存储介质,其特征在于,所述计算机程序被处理器执行时实现权利要求1至2中任一项所述的多应力加速试验剖面确定方法的步骤。
6.一种多应力加速试验剖面确定设备,其特征在于,包括:控制器;
所述控制器包括权利要求5所述的存储有计算机程序的计算机可读存储介质,所述计算机程序被处理器执行时实现权利要求1至2中任一项所述的多应力加速试验剖面确定方法的步骤;或者,所述控制器包括权利要求3至4中任一项所述的多应力加速试验剖面确定系统。
CN201911018342.5A 2019-10-24 2019-10-24 多应力加速试验剖面确定方法、系统、介质及设备 Active CN110686915B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911018342.5A CN110686915B (zh) 2019-10-24 2019-10-24 多应力加速试验剖面确定方法、系统、介质及设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911018342.5A CN110686915B (zh) 2019-10-24 2019-10-24 多应力加速试验剖面确定方法、系统、介质及设备

Publications (2)

Publication Number Publication Date
CN110686915A CN110686915A (zh) 2020-01-14
CN110686915B true CN110686915B (zh) 2021-05-25

Family

ID=69114018

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911018342.5A Active CN110686915B (zh) 2019-10-24 2019-10-24 多应力加速试验剖面确定方法、系统、介质及设备

Country Status (1)

Country Link
CN (1) CN110686915B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111062183A (zh) * 2020-03-06 2020-04-24 北京强度环境研究所 多应力强化试验剖面编制方法、装置及计算机存储介质
CN111208447A (zh) * 2020-03-06 2020-05-29 上海度普新能源科技有限公司 功率模块测试系统
CN111579185B (zh) * 2020-04-08 2022-05-03 中国电子产品可靠性与环境试验研究所((工业和信息化部电子第五研究所)(中国赛宝实验室)) 潜艇用电子设备可靠性加速试验方法和装置
CN112034345B (zh) * 2020-09-28 2023-05-16 上海电气集团股份有限公司 一种车用电机高温耐久测试方法
CN114323427A (zh) * 2021-12-24 2022-04-12 中国电子科技集团公司第四十九研究所 一种基于压力传感器的多应力加速贮存试验方法
CN114414463B (zh) * 2021-12-28 2023-11-14 北京遥感设备研究所 一种铝基复合材料光学系统长期贮存稳定性验证方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102445338B (zh) * 2011-11-24 2013-08-14 北京航空航天大学 一种航天驱动组件综合应力加速寿命试验方法
CN103364703B (zh) * 2012-11-27 2015-07-29 深圳清华大学研究院 Led多应力下可靠性快速评测方法
CN103246821B (zh) * 2013-05-21 2016-08-03 北京航空航天大学 一种基于仿真的多应力小样本加速寿命试验方案设计优化方法
CN104251784B (zh) * 2013-06-27 2017-09-01 上海机电工程研究所 机电产品整机综合应力可靠性加速试验方法
CN103616326B (zh) * 2013-12-10 2017-06-20 中国人民解放军军械工程学院 通过温湿度、电应力加速退化试验获得雷达寿命的方法
CN104156599B (zh) * 2014-08-13 2017-03-08 北京航空航天大学 一种考虑多机理的加速寿命试验载荷谱确定方法
CN104298853B (zh) * 2014-08-26 2017-04-05 北京航空航天大学 一种多应力多机理作用下柱塞泵产品综合加速因子优化方法
CN104156621B (zh) * 2014-08-26 2017-04-26 北京航空航天大学 一种基于理论模型的柱塞泵加速因子确定方法
CN106021928B (zh) * 2016-05-20 2019-02-22 中国人民解放军总参谋部第六十研究所 一种综合应力加速试验方法
US10739397B2 (en) * 2017-05-10 2020-08-11 International Business Machines Corporation Accelerated wafer testing using non-destructive and localized stress
CN107390668B (zh) * 2017-06-26 2019-07-16 西安航空制动科技有限公司 确定防滑刹车控制装置加严可靠性试验剖面的方法
CN110362957B (zh) * 2019-08-01 2022-05-03 西北工业大学 一种航空发动机关键件寿命管理方法及装置

Also Published As

Publication number Publication date
CN110686915A (zh) 2020-01-14

Similar Documents

Publication Publication Date Title
CN110686915B (zh) 多应力加速试验剖面确定方法、系统、介质及设备
CN108204925B (zh) 复合材料的疲劳寿命预测方法及预测系统
Agius et al. Sensitivity and optimisation of the Chaboche plasticity model parameters in strain-life fatigue predictions
US11115295B2 (en) Methods and systems for online monitoring using a variable data
CN103310051B (zh) 一种全寿命周期内的车载信息终端故障率预测方法
CN107784178B (zh) 一种基于多故障机理耦合的机械结构可靠性分析方法
CN101620045A (zh) 基于时间序列的步进应力加速退化试验可靠性评估方法
CN111382546B (zh) 一种基于支持向量机建模对发电机绝缘系统进行寿命预测的方法
KR101526313B1 (ko) 피로수명 예측방법
WO2008074526A3 (en) Method, system and computer program for testing software applications based on multiple data sources
Liu et al. Steady ratcheting strains accumulation in varying temperature fatigue tests of PMMA
CN109213098B (zh) 调整运行参数的方法、装置、电子设备和计算机可读介质
CN104063569A (zh) 一种基于emd去噪和渐消记忆的设备剩余寿命预测方法
CN108228979B (zh) 一种用于发电机结构的低周疲劳寿命分析方法
EP2988186B1 (en) Detection of consistent steady-state behavior in process plants
Lin et al. A novel product remaining useful life prediction approach considering fault effects
CN111475925A (zh) 一种电力设备的状态评估方法、装置及存储介质
Bogdanov et al. Probabilistic analysis of the fatigue crack growth based on the application of the monte-carlo method to unigrow model
Sunder Spectrum load fatigue—underlying mechanisms and their significance in testing and analysis
CN103488826A (zh) 基于经验加速模型的退化量分布参数建模外推方法
KR101658748B1 (ko) 신뢰성 설계 지원 장치 및 신뢰성 설계 지원 방법
CN1249421C (zh) 基于灰色模型的动调陀螺仪寿命预测方法
KR101293982B1 (ko) 탄성중합체의 시뮬레이션 방법
CN107290603B (zh) 一种产品可靠性评价方法及装置
Lapira et al. A systematic approach to intelligent maintenance of production systems with a framework for embedded implementation

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