CN113239556B - 任意加速度功率谱密度的疲劳损伤速率倍数估算方法 - Google Patents

任意加速度功率谱密度的疲劳损伤速率倍数估算方法 Download PDF

Info

Publication number
CN113239556B
CN113239556B CN202110558067.7A CN202110558067A CN113239556B CN 113239556 B CN113239556 B CN 113239556B CN 202110558067 A CN202110558067 A CN 202110558067A CN 113239556 B CN113239556 B CN 113239556B
Authority
CN
China
Prior art keywords
spectral density
power spectral
natural frequency
damping ratio
acceleration power
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
CN202110558067.7A
Other languages
English (en)
Other versions
CN113239556A (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.)
General Engineering Research Institute China Academy of Engineering Physics
Original Assignee
General Engineering Research Institute China Academy of Engineering Physics
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 General Engineering Research Institute China Academy of Engineering Physics filed Critical General Engineering Research Institute China Academy of Engineering Physics
Priority to CN202110558067.7A priority Critical patent/CN113239556B/zh
Publication of CN113239556A publication Critical patent/CN113239556A/zh
Application granted granted Critical
Publication of CN113239556B publication Critical patent/CN113239556B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/08—Probabilistic or stochastic CAD
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/04—Ageing analysis or optimisation against ageing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明公开了任意加速度功率谱密度的疲劳损伤速率倍数估算方法，包括利用采集的振动加速度数据通过自谱分析得到加速度功率谱密度A1(f)和A2(f)；根据***的固有频率和阻尼比将加速度功率谱密度A1(f)和A2(f)转换为应力功率谱密度G1(f)和G2(f)；选取能够覆盖对象实际阻尼比和固有频率的区间，将阻尼比区间和固有频率区间进行等分，形成离散的阻尼比ξ数列和固有频率fn数列，再采用正交法将两个数列进行组合，计算出不同fn和β组合下的一系列
Figure DDA0003078056020000011
和
Figure DDA0003078056020000012
根据得到的一系列
Figure DDA0003078056020000013
和
Figure DDA0003078056020000014
分别计算出一系列的
Figure DDA0003078056020000015
和
Figure DDA0003078056020000016
的比值，取最小值作为疲劳损伤速率倍数的最终估计值。通过本发明可以实现对任意加速度功率谱密度曲线之间的损伤速率倍数的估算，进而提高加速试验中试验时间设计的准确度。

Description

任意加速度功率谱密度的疲劳损伤速率倍数估算方法
技术领域
本发明涉及，具体是任意加速度功率谱密度的疲劳损伤速率倍数估算方法。
背景技术
在基于随机振动环境的加速试验中，需要计算不同加速度功率谱密度条件下疲劳损伤速率倍数，以获得加速后的试验时间。
传统方法只能计算相同谱形、不同量值加速度功率谱密度条件下的疲劳损伤速率倍数。在线弹性范围内，在不改变加速度功率谱形状、仅提高谱密度量值时，应力循环的分布不变，应力值表现为整体增大，均方根应力增大的比例与加速度均方根增大的比例相同。功率谱密度比值是均方根比值的二次方。所以，在加速度功率谱的谱形相同的情况下，很容易得到不同谱密度量值下的应力比例关系，已知SN曲线斜率就可以计算出疲劳损伤速率倍数。
对于谱形不同的情况则需要掌握对象的详细结构模型参数，通过有限元法或简化动力学模型等转换为应力功率谱密度或者应力时域历程才能进行计算。应力功率谱密度已知时可以采用各种频域疲劳损伤计算方法，如窄带法、DK法等。应力时域历程已知时可以采用雨流计数法等。分别计算得到不同条件下的疲劳损伤速率，从而得到损伤速率倍数，但传统方法只能计算相同谱形、不同量值加速度功率谱密度条件下的疲劳损伤速率倍数。而对于谱形不同的情况则需要转换为应力功率谱密度才能进行计算，一般需要掌握对象的结构模型参数。在结构模型参数未知的情况下是无法进行疲劳损伤速率的计算。
发明内容
本发明的目的在于克服现有技术的不足，提供任意加速度功率谱密度的疲劳损伤速率倍数估算方法，包括如下步骤：
步骤一，利用采集的振动加速度数据通过自谱分析得到加速度功率谱密度A1(f)和A2(f)；
步骤二，根据***的固有频率和阻尼比将加速度功率谱密度A1(f)和A2(f)转换为应力功率谱密度G1(f)和G2(f)；选取覆盖对象实际阻尼比和固有频率的区间，将阻尼比区间和固有频率区间进行等分，形成离散的阻尼比ξ数列和固有频率fn数列，再采用正交法将两个数列进行组合，计算损伤损率时对全部组合进行遍历；
步骤三，利用频域疲劳寿命预测法，计算出应力的概率密度分布p(S)和期望峰值率Vp，并计算出单位时间内的累积损伤和计算出不同fn和β组合下的一系列
Figure BDA0003078056000000011
和/>
Figure BDA0003078056000000012
步骤四，根据得到的一系列
Figure BDA0003078056000000013
和/>
Figure BDA0003078056000000014
分别计算出一系列的/>
Figure BDA0003078056000000015
和/>
Figure BDA0003078056000000016
的比值，取最小值作为疲劳损伤速率倍数的最终估计值。
进一步的，所述的根据***的固有频率和阻尼比将加速度功率谱密度A1(f)和A2(f)转换为应力功率谱密度G1(f)和G2(f)，所采用的传递函数形式为一阶***基础激励下的响应函数：
G(f)＝H(f)·A(f)
Figure BDA0003078056000000021
β＝f/fn
式中：H为加速度功率谱密度到应力功率谱密度的传递函数；β为频率比；fn为对象的固有频率；ξ为阻尼比，取值范围为[0.01，0.1]；
所述的计算出不同fn和β组合下的一系列
Figure BDA0003078056000000022
和/>
Figure BDA0003078056000000023
采用如下计算公式：
Figure BDA0003078056000000024
式中：S为应力幅值；k为材料S-N曲线斜率；N为循环次数；C为材料常数。
进一步的，所述的计算出一系列的
Figure BDA0003078056000000026
和/>
Figure BDA0003078056000000027
的比值，采用如下公式：
Figure BDA0003078056000000025
本发明的有益效果是：不需要知道对象结构模型参数，可以实现对任意加速度功率谱密度曲线之间的损伤速率倍数的估算，进而提高加速试验中试验时间设计的准确度。
附图说明
图1为任意加速度功率谱密度的疲劳损伤速率倍数估算方法的流程示意图。
具体实施方式
下面结合附图进一步详细描述本发明的技术方案，但本发明的保护范围不局限于以下所述。
如图1所示，任意加速度功率谱密度的疲劳损伤速率倍数估算方法，包括如下步骤：
步骤一，利用采集的振动加速度数据通过自谱分析得到加速度功率谱密度A1(f)和A2(f)；
步骤二，根据***的固有频率和阻尼比将加速度功率谱密度A1(f)和A2(f)转换为应力功率谱密度G1(f)和G2(f)；选取覆盖对象实际阻尼比和固有频率的区间，将阻尼比区间和固有频率区间进行等分，形成离散的阻尼比ξ数列和固有频率fn数列，再采用正交法将两个数列进行组合，计算损伤损率时对全部组合进行遍历；
步骤三，利用频域疲劳寿命预测法，计算出应力的概率密度分布p(S)和期望峰值率Vp，并计算出单位时间内的累积损伤和计算出不同fn和β组合下的一系列
Figure BDA0003078056000000031
和/>
Figure BDA0003078056000000032
步骤四，根据得到的一系列
Figure BDA0003078056000000033
和/>
Figure BDA0003078056000000034
分别计算出一系列的/>
Figure BDA0003078056000000035
和/>
Figure BDA0003078056000000036
的比值，取最小值作为疲劳损伤速率倍数的最终估计值。
所述的根据***的固有频率和阻尼比将加速度功率谱密度A1(f)和A2(f)转换为应力功率谱密度G1(f)和G2(f)，所采用的传递函数形式为一阶***基础激励下的响应函数：
G(f)＝H(f)·A(f)
Figure BDA0003078056000000037
β＝f/fn
式中：H为加速度功率谱密度到应力功率谱密度的传递函数；β为频率比；fn为对象的固有频率；ξ为阻尼比，取值范围为[0.01，0.1]；
所述的计算出不同fn和β组合下的一系列
Figure BDA0003078056000000038
和/>
Figure BDA0003078056000000039
采用如下计算公式：
Figure BDA00030780560000000310
式中：S为应力幅值；k为材料S-N曲线斜率；N为循环次数；C为材料常数。
所述的计算出一系列的
Figure BDA00030780560000000311
和/>
Figure BDA00030780560000000312
的比值，采用如下公式：
Figure BDA00030780560000000313
具体的，利用实测振动加速度数据通过自谱分析得到加速度功率谱密度A1(f)和A2(f)，f为频率、单位为Hz。
根据***的固有频率和阻尼比将加速度功率谱密度A1(f)和A2(f)转换为应力功率谱密度G1(f)和G2(f)，采用的传递函数形式为一阶***基础激励下的响应函数。
Figure BDA0003078056000000041
式中：H-加速度功率谱密度到应力功率谱密度的传递函数；β-频率比；fn-对象的固有频率；ξ-阻尼比，一般取值范围为[0.01，0.1]。
对象的阻尼比和固有频率准确值未知时，根据相似产品的已有数据(可以是实测值、也可以是仿真计算结果)，已知对象阻尼比和固有频率的大致范围。选取一个能够覆盖对象实际阻尼比和固有频率的区间。将阻尼比区间和固有频率区间进行等分，形成离散的阻尼比数列和固有频率数列，再采用正交法将两个数列进行组合，计算损伤损率时对全部组合进行遍历。
例如，阻尼比大约在0.01～0.05之间，可以划分为[0.01，0.02，0.03，0.04，0.05]的阻尼比数列；固有频率大约在80Hz～160Hz之间，可以划分为[80，90，100，110，120，130，140，150，160]的固有频率数列。正交考虑即有5×9＝45种组合。
利用频域疲劳寿命预测法，如窄带分布法、Wirsching-Light法(WL法)、Tovo-Benasciutti(TB法)法和Dirlik(DK法)等，计算出应力的概率密度分布p(S)和期望峰值率Vp。按照Miner线性累积损伤理论，计算出单位时间内的累积损伤，这里也称损伤速率。计算出不同fn和β组合下的一系列
Figure BDA0003078056000000042
和/>
Figure BDA0003078056000000043
Figure BDA0003078056000000044
设fn有i个数，β有j个数
分别计算G1(f)i，j和G2(f)i，j，也即第i个fn和第j个β对应的应力功率谱密度，计算p1(S)i，j和
Figure BDA0003078056000000045
得到/>
Figure BDA0003078056000000046
计算/>
Figure BDA0003078056000000047
和/>
Figure BDA0003078056000000048
得到/>
Figure BDA0003078056000000049
式中：S-应力幅值；k-材料S-N曲线斜率；N-循环次数；C-材料常数计算损伤速率倍数时会消掉，可用单位1代替。
计算出一系列的
Figure BDA00030780560000000410
和/>
Figure BDA00030780560000000411
的比值，取最小值作为疲劳损伤速率倍数的最终估计值。
Figure BDA00030780560000000412
以上所述仅是本发明的优选实施方式，应当理解本发明并非局限于本文所披露的形式，不应看作是对其他实施例的排除，而可用于各种其他组合、修改和环境，并能够在本文所述构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围，则都应在本发明所附权利要求的保护范围内。

Claims (3)

1.任意加速度功率谱密度的疲劳损伤速率倍数估算方法，其特征在于，包括如下步骤：
步骤一，利用采集的振动加速度数据通过自谱分析得到加速度功率谱密度A1(f)和A2(f)；
步骤二，根据***的固有频率和阻尼比将加速度功率谱密度A1(f)和A2(f)转换为应力功率谱密度G1(f)和G2(f)；选取覆盖对象实际阻尼比和固有频率的区间，将阻尼比区间和固有频率区间进行等分，形成离散的阻尼比ξ数列和固有频率fn数列，再采用正交法将两个数列进行组合，计算损伤损率时对全部组合进行遍历；
步骤三，利用频域疲劳寿命预测法，计算出应力的概率密度分布p(S)和期望峰值率Vp，并计算出单位时间内的累积损伤和计算出不同fn和β组合下的一系列
Figure FDA0004160017360000011
和/>
Figure FDA0004160017360000012
步骤四，根据得到的一系列
Figure FDA0004160017360000013
和/>
Figure FDA0004160017360000014
分别计算出一系列的/>
Figure FDA0004160017360000015
和/>
Figure FDA0004160017360000016
的比值，取最小值作为疲劳损伤速率倍数的最终估计值。
2.权利要求1所述的任意加速度功率谱密度的疲劳损伤速率倍数估算方法，其特征在于，所述的根据***的固有频率和阻尼比将加速度功率谱密度A1(f)和A2(f)转换为应力功率谱密度G1(f)和G2(f)，所采用的传递函数形式为一阶***基础激励下的响应函数：
G(f)＝H(f)·A(f)
Figure FDA0004160017360000017
β＝f/fn
式中：H为加速度功率谱密度到应力功率谱密度的传递函数；β为频率比；fn为对象的固有频率；ξ为阻尼比，取值范围为[0.01,0.1]；
所述的计算出不同fn和β组合下的一系列
Figure FDA0004160017360000018
和/>
Figure FDA0004160017360000019
采用如下计算公式：
Figure FDA00041600173600000110
式中：S为应力幅值；k为材料S-N曲线斜率；C为材料常数。
3.权利要求2所述的任意加速度功率谱密度的疲劳损伤速率倍数估算方法，其特征在于，所述的计算出一系列的
Figure FDA00041600173600000111
和/>
Figure FDA00041600173600000112
的比值，采用如下公式：
Figure FDA00041600173600000113
CN202110558067.7A 2021-05-21 2021-05-21 任意加速度功率谱密度的疲劳损伤速率倍数估算方法 Active CN113239556B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110558067.7A CN113239556B (zh) 2021-05-21 2021-05-21 任意加速度功率谱密度的疲劳损伤速率倍数估算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110558067.7A CN113239556B (zh) 2021-05-21 2021-05-21 任意加速度功率谱密度的疲劳损伤速率倍数估算方法

Publications (2)

Publication Number Publication Date
CN113239556A CN113239556A (zh) 2021-08-10
CN113239556B true CN113239556B (zh) 2023-07-14

Family

ID=77138031

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110558067.7A Active CN113239556B (zh) 2021-05-21 2021-05-21 任意加速度功率谱密度的疲劳损伤速率倍数估算方法

Country Status (1)

Country Link
CN (1) CN113239556B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117686935B (zh) * 2024-01-30 2024-05-14 湖北工业大学 一种基于电压概率密度的电池rul预测方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106482827A (zh) * 2016-11-11 2017-03-08 北京航空航天大学 基于互功率谱函数模态参数识别的电子产品振动dlp方法
CN107103162A (zh) * 2017-05-26 2017-08-29 中国人民解放军国防科学技术大学 一种基于疲劳损伤累积理论的振动加速试验方法及***

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NO3074325T3 (zh) * 2013-12-23 2018-02-24
CN109857977B (zh) * 2019-03-08 2022-10-18 北京工业大学 一种变温振动下基于频域的疲劳寿命计算方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106482827A (zh) * 2016-11-11 2017-03-08 北京航空航天大学 基于互功率谱函数模态参数识别的电子产品振动dlp方法
CN107103162A (zh) * 2017-05-26 2017-08-29 中国人民解放军国防科学技术大学 一种基于疲劳损伤累积理论的振动加速试验方法及***

Also Published As

Publication number Publication date
CN113239556A (zh) 2021-08-10

Similar Documents

Publication Publication Date Title
Yang et al. The role of the turbulent scales in the settling velocity of heavy particles in homogeneous isotropic turbulence
CN105008887B (zh) 使用非接触测量和动态响应重构技术的涡轮机叶片疲劳寿命分析
CN110595765A (zh) 基于vmd和fa_pnn风电机组齿轮箱故障诊断方法
CN109212966B (zh) 一种多工况动态基准化的机械设备剩余寿命预测方法
CN107885928B (zh) 考虑测量误差的步进应力加速性能退化可靠性分析方法
CN107436983B (zh) 一种基于多元样本差异的o型橡胶密封圈寿命预测方法
CN113239556B (zh) 任意加速度功率谱密度的疲劳损伤速率倍数估算方法
CN107084709B (zh) 一种多弹性对径流变化驱动因素的定量分割方法
KR101452171B1 (ko) 타워형 구조물의 변위 추정 방법
CN108629092B (zh) 一种基于尺寸效应修正的涡轮盘分区可靠性分析方法
CN112697691A (zh) 一种航空合金钢服役环境下腐蚀损伤分布规律的研究方法
CN107045659A (zh) 预测光伏发电功率的方法及装置
CN107764257A (zh) 一种惯性器件数值模拟方法
CN106326530A (zh) 一种基于右矩阵分式模型的时变结构模态参数辨识方法
Schwartz Reply to comments by G. Foster et al., R. Knutti et al., and N. Scafetta on “Heat capacity, time constant, and sensitivity of Earth’s climate system”
CN116307039A (zh) 一种计及气象差异性的光伏出力智能预测方法
CN106990011A (zh) 一种加速寿命试验的温度加速应力谱的确定方法
CN114046903A (zh) 一种高压电缆缆芯的实时温度预测方法及***
CN111144230B (zh) 一种基于vmd的时域载荷信号的去噪方法
CN113486295A (zh) 基于傅里叶级数的臭氧总量变化预测方法
Bogdanoff et al. On nonstationary cumulative damage models
CN105987821B (zh) 低气压环境应力下单机产品的加速贮存试验方法
Abbasi et al. On enhanced sensitivity of nonparametric EWMA control charts for process monitoring
KR102230397B1 (ko) 변형률 및 가속도값에 기반한 구조물의 변위 추정방법
CN114564487A (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