CN111159871A - 一种基于路径曲线积分的随机多轴循环计数方法 - Google Patents
一种基于路径曲线积分的随机多轴循环计数方法 Download PDFInfo
- Publication number
- CN111159871A CN111159871A CN201911338574.9A CN201911338574A CN111159871A CN 111159871 A CN111159871 A CN 111159871A CN 201911338574 A CN201911338574 A CN 201911338574A CN 111159871 A CN111159871 A CN 111159871A
- Authority
- CN
- China
- Prior art keywords
- point
- time
- counting method
- axis
- strain
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000010354 integration Effects 0.000 title claims abstract description 11
- 230000006378 damage Effects 0.000 claims abstract description 13
- DMSMPAJRVJJAGA-UHFFFAOYSA-N benzo[d]isothiazol-3-one Chemical compound C1=CC=C2C(=O)NSC2=C1 DMSMPAJRVJJAGA-UHFFFAOYSA-N 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 3
- 208000027418 Wounds and injury Diseases 0.000 claims 1
- 208000014674 injury Diseases 0.000 claims 1
- 239000000463 material Substances 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 2
- 238000009661 fatigue test Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 229910000851 Alloy steel Inorganic materials 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
Images
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
Description
技术领域
本发明涉及疲劳强度领域,特指一种基于路径曲线积分的随机多轴循环计数方法。
背景技术
很多机器在服役时承受复杂的载荷状态,如航空航天飞行器、高速列车、大型搬运机械,其零部件关键部位危险点处于多轴应力状态,应力分量通常随时间任意变化并且幅值很大,疲劳裂纹容易萌生和扩展,导致服役寿命急剧降低,所以精准预测变幅随机多轴疲劳寿命意义重大。
过去几十年来,单轴疲劳强度理论发展逐渐完善,但多轴疲劳强度理论还需要补充和优化。预测随机多轴疲劳寿命可以简单分为三个步骤:1)通过计数方法,将随机变幅载荷历程分解成多个恒幅载荷历程,即多个半循环(反复)。2)利用合适的损伤模型,计算每个反复的损伤参量,进而得出损伤。3)借助损伤累积准则,将每个反复的损伤累加为总损伤,估算疲劳寿命。计数方法是变幅随机多轴寿命预测的前处理流程,其合理性直接影响疲劳寿命评估精度。现存的计数方法有很多缺点,如多轴雨流计数方法可能会丢失辅助计数通道的峰谷值、基于等效应变的雨流计数方法会遗漏谷值符号、依赖损伤模型的计数方法不具有通用性。因此,提出不依赖损伤模型且综合考虑分量通道信息的计数方法具有重要意义。
发明内容
本发明针对变幅多轴疲劳寿命预测的发展,考虑了加载路径对于疲劳寿命的影响,克服了计数方法对疲劳损伤模型的依赖性和遗漏计数通道谷值符合的缺点,提出了一种基于路径曲线积分的随机多轴疲劳循环计数方法。
本发明所提供的一种基于路径曲线积分的随机多轴疲劳循环计数方法,其步骤为:
其中是起始时刻1的投影点和时刻t的投影点之间的加载路径,是起始时刻start的投影点指向时刻t投影点的向量,是时刻t-1的投影点指向时刻t投影点的向量,ds是空间中加载路径的弧长微元,sign(x)是符号函数,表达式如下
其中x为任意实数。
步骤3):计算所有时刻点的积分值Yt,找出对应Yt的绝对值abs(Yt)最大的投影点Am,将Am点前面和后面的历程调换顺序,将Am置为起始点Astart。
步骤4):计算新序列所有时刻点的Yt,找出对应Yt最大的投影点Am,将起始点Astart和点Am之间的历程记为一个半循环(反复)。
步骤5):对于点Am和终止点Aend之间的历程,首先将Am置为起始点,然后按照步骤4)处理。
步骤6):若在步骤4)中出现Yt值非单调增加的子历程,将该子历程按照步骤4)、5)进行处理,直至计出所有的半循环(反复)。
与现有技术相比,本发明具有如下有益效果:
本发明提出一种基于路径曲线积分的随机多轴循环计数方法,该方法通过计算应变空间中各点相对于参考起始点的曲线积分值Yt,找出积分值最大的投影点Am,将A0与Am之间的历程计为一个反复,对未处理的历程和积分值Yt出现转折的历程按照上述步骤进行递归处理,直至计出所有反复。计数过程中不需要考虑材料常数,对损伤模型无依赖性,解决了基于等效应变的计数方法丢失载荷谷值符号的问题,在单轴状态下可退化到雨流计数法,易于程序实现,便于工程计算。
附图说明
图1为基于路径曲线积分的随机多轴计数方法流程图。
图2为待计数的随机多轴疲劳试验的拉-扭载荷时间历程。
图4为计数结果及每个反复相对计出的反复起点的Yt-t图。
具体实施方式
结合附图说明本发明的具体实施方式。
本发明通过拉-扭变幅随机多轴疲劳试验做进一步说明,试件材料为En15R合金钢,采用控制应变加载,加载波形为斜波,轴向应变幅为0.80%,扭向应变幅为1.37%,试验采集了341个数据点。
本发明所提供的一种基于路径曲线积分的随机多轴疲劳循环计数方法,其步骤为:
其中是起始时刻1的投影点和时刻t的投影点之间的加载路径,是起始时刻1的投影点指向时刻t的投影点的向量,是时刻t-1的投影点指向时刻t的投影点的向量,ds是空间中加载路径的弧长微元,sign(x)是符号函数,表达式如下
其中x为任意实数。
步骤3):按照时间先后顺序将投影点编号为A1,A1,…,A340,A341,计算所有时刻点的Yt值,找出对应Yt的绝对值abs(Yt)最大的时刻300,将载荷重排列为A300,…A341,A1,…,A299,将A300置为A1,得到新序列A1,A1,…,A340,A341。
步骤4):对于新序列,计算每个时刻点的积分值Yt,找出对应Yt最大的投影点Am(第一次执行本步骤时m=68),将A1和Am之间的历程计为一个半循环。这里m的值会随着递归次数更新。
步骤5):将Am和Aend之间的历程重新编号为A1,…,Aend(第一次执行本步骤时end=341-68=274),按照步骤3)处理。这里end的值会随着递归次数更新。
Claims (2)
1.一种基于路径曲线积分的随机多轴循环计数方法,其特征在于:该方法的实现步骤如下,
其中是起始时刻1的投影点和时刻t的投影点之间的加载路径,是起始时刻start的投影点指向时刻t投影点的向量,是时刻t-1的投影点指向时刻t投影点的向量,ds是空间中加载路径的弧长微元,sign(x)是符号函数,表达式如下
其中x为任意实数;
步骤3):计算所有时刻点的积分值Yt,找出对应Yt的绝对值abs(Yt)最大的投影点Am,将Am点前面和后面的历程调换顺序,将Am置为起始点Astart;
步骤4):计算新序列所有时刻点的Yt,找出对应Yt最大的投影点Am,将起始点Astart和点Am之间的历程记为一个半循环;
步骤5):对于点Am和终止点Aend之间的历程,首先将Am置为起始点,然后按照步骤4)处理;
步骤6):若在步骤4)中出现Yt值非单调增加的子历程,将该子历程按照步骤4)、5)进行处理,直至计出所有的半循环。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911338574.9A CN111159871B (zh) | 2019-12-23 | 2019-12-23 | 一种基于路径曲线积分的随机多轴循环计数方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911338574.9A CN111159871B (zh) | 2019-12-23 | 2019-12-23 | 一种基于路径曲线积分的随机多轴循环计数方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111159871A true CN111159871A (zh) | 2020-05-15 |
CN111159871B CN111159871B (zh) | 2024-03-26 |
Family
ID=70557797
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911338574.9A Active CN111159871B (zh) | 2019-12-23 | 2019-12-23 | 一种基于路径曲线积分的随机多轴循环计数方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111159871B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112100563A (zh) * | 2020-09-11 | 2020-12-18 | 广州汽车集团股份有限公司 | 多轴载荷等效处理方法、装置、计算机设备和介质 |
CN112749683A (zh) * | 2021-01-27 | 2021-05-04 | 吉林大学 | 一种保留载荷时序的雨流计数方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110288790A1 (en) * | 2009-01-08 | 2011-11-24 | Battelle Memorial Institute | Path-Dependent Cycle Counting And Multi-Axial Fatigue Evaluation Of Engineering Structures |
JP2013044667A (ja) * | 2011-08-25 | 2013-03-04 | Ihi Corp | 多軸疲労寿命評価方法 |
CN106600066A (zh) * | 2016-12-19 | 2017-04-26 | 南京理工大学 | 一种基于scada数据的风力发电机齿轮箱疲劳寿命估计方法 |
CN109883709A (zh) * | 2019-03-08 | 2019-06-14 | 北京工业大学 | 一种基于相对等效应变的随机多轴热机计数方法 |
-
2019
- 2019-12-23 CN CN201911338574.9A patent/CN111159871B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110288790A1 (en) * | 2009-01-08 | 2011-11-24 | Battelle Memorial Institute | Path-Dependent Cycle Counting And Multi-Axial Fatigue Evaluation Of Engineering Structures |
JP2013044667A (ja) * | 2011-08-25 | 2013-03-04 | Ihi Corp | 多軸疲労寿命評価方法 |
CN106600066A (zh) * | 2016-12-19 | 2017-04-26 | 南京理工大学 | 一种基于scada数据的风力发电机齿轮箱疲劳寿命估计方法 |
CN109883709A (zh) * | 2019-03-08 | 2019-06-14 | 北京工业大学 | 一种基于相对等效应变的随机多轴热机计数方法 |
Non-Patent Citations (3)
Title |
---|
乐京霞;王思宇;刘玉亮;陈鹏飞;姚国全;: "FPSO单点系泊***导缆器疲劳寿命预报方法", 中国舰船研究, no. 02 * |
于强;: "多轴非比例加载下疲劳短裂纹扩展速率的探讨", 机械强度, no. 01 * |
孙斌斌;王英玉;: "一个新的多轴疲劳循环计数法", 江苏航空, no. 1 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112100563A (zh) * | 2020-09-11 | 2020-12-18 | 广州汽车集团股份有限公司 | 多轴载荷等效处理方法、装置、计算机设备和介质 |
CN112100563B (zh) * | 2020-09-11 | 2024-04-19 | 广州汽车集团股份有限公司 | 多轴载荷等效处理方法、装置、计算机设备和介质 |
CN112749683A (zh) * | 2021-01-27 | 2021-05-04 | 吉林大学 | 一种保留载荷时序的雨流计数方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111159871B (zh) | 2024-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106840877B (zh) | 一种基于应力的多轴小裂纹全寿命预测方法 | |
CN106096073B (zh) | 一种基于损伤力学非概率区间分析模型的金属疲劳裂纹全寿命预估方法 | |
Carpenter et al. | Accuracy of shock capturing in two spatial dimensions | |
CN112487683B (zh) | 一种考虑残余应力影响的结构件高周疲劳寿命预测方法 | |
CN111159871A (zh) | 一种基于路径曲线积分的随机多轴循环计数方法 | |
Colledani et al. | Design and evaluation of in-line product repair strategies for defect reduction in the production of electric drives | |
KR102030213B1 (ko) | 성형 강판 패널들의 스냅-스루 좌굴 예측 시스템 및 방법 | |
CN109657412B (zh) | 一种陶瓷基复合材料的疲劳寿命的多尺度预测方法 | |
CN112180065A (zh) | 一种金属材料的j-c本构模型建立方法 | |
CN116629509A (zh) | 一种基于关键机器的流水车间成组调度快速评价方法 | |
CN113948163B (zh) | 一种用于重复使用火箭发动机涡轮泵的高低周复合疲劳寿命预测方法 | |
US11554433B2 (en) | Methods and systems for wire electric discharge machining and validation of parts | |
CN111859729A (zh) | 考虑多弹丸随机分布的喷丸模型对轮盘寿命的计算方法 | |
Colledani et al. | Analysis of In-line Quality-oriented Assembly Strategies in the Production of Electric Drives | |
US11288424B2 (en) | Adaptive time scaling for accelerating explicit finite element analysis | |
CN110287619B (zh) | 一种多轴随机载荷疲劳寿命预测方法 | |
Kurek et al. | Including of ratio of fatigue limits from bending and torsion for estimation fatigue life under cyclic loading | |
Maywald et al. | Prediction of Geometrically Induced Localization Effects Using a Subset of Nominal System Modes | |
Kirthan et al. | Computational analysis of fatigue crack growth based on stress intensity factor approach in axial flow compressor blades | |
Arévalo et al. | XFEM-based fracture mechanics analysis of aeronautical structures affected by residual stresses | |
Viharos et al. | Support Vector Machine (SVM) based general model building algorithm for production control | |
Wołowiec et al. | Simulation and control of tool steel quenching process | |
Al-Asady et al. | FEA based durability using strain-life models for different medium carbon steel as fabrication materials for an automotive component | |
Ricardo | Influence of node release in crack propagation simulation under variable amplitude loading | |
CN112180886B (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 |