CN110705137B - 一种应力幅值和均值的确定方法以及装置 - Google Patents

一种应力幅值和均值的确定方法以及装置 Download PDF

Info

Publication number
CN110705137B
CN110705137B CN201910780182.1A CN201910780182A CN110705137B CN 110705137 B CN110705137 B CN 110705137B CN 201910780182 A CN201910780182 A CN 201910780182A CN 110705137 B CN110705137 B CN 110705137B
Authority
CN
China
Prior art keywords
stress
maximum
tensor
liangcha
zhang
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
CN201910780182.1A
Other languages
English (en)
Other versions
CN110705137A (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.)
CRRC Qingdao Sifang Co Ltd
Original Assignee
CRRC Qingdao Sifang Co Ltd
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 CRRC Qingdao Sifang Co Ltd filed Critical CRRC Qingdao Sifang Co Ltd
Priority to CN201910780182.1A priority Critical patent/CN110705137B/zh
Publication of CN110705137A publication Critical patent/CN110705137A/zh
Application granted granted Critical
Publication of CN110705137B publication Critical patent/CN110705137B/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
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本申请公开了一种应力幅值和均值的确定方法以及装置,在获取到待评价对象在各个工况下的应力张量后,获取每两个不同工况下的应力张量差,并从各个应力张量差中确定最大应力张量差以及最大应力张量差对应的两种工况,并根据最大应力张量差以及最大应力张量差对应的两种工况,确定待评价对象的应力幅值及应力均值。由于每个应力张量差能够准确地表征从一种工况过渡到另一种工况时产生的应力变化,使得最大应力张量差能够准确地表征待评价对象的最大应力变量,即应力幅值,因而,根据最大应力张量差及对应工况确定的应力幅值和应力均值能够准确地表征待评价对象的最大应力幅值及均值,提高应力幅值和均值的准确性,并提高结构疲劳强度的评估准确性。

Description

一种应力幅值和均值的确定方法以及装置
技术领域
本申请涉及疲劳强度评价技术领域,尤其涉及一种应力幅值和均值的确定方法以及装置。
背景技术
目前,通常使用Goodman方法评价疲劳强度;而且,在利用Goodman方法进行疲劳强度评价之前需要获取待评价对象的应力幅值和均值。然而,对于多轴疲劳评价来说,因多个工况之间的主应力方向不同而导致应力幅值和均值的获取是十分困难的。
在现有技术中,通常使用最大主应力法和最大主应力投影法来确定应力幅值和均值。其中,最大主应力法和最大主应力投影法均基于“最大幅值方向就是最大主应力方向”的假设来确定应力幅值和均值,但是在实际应用中最大幅值方向与最大主应力方向重合的概率极低(例如,当压应力绝对值远大于拉应力时,最大幅值方向与最大主应力方向是不重合的),如此导致利用最大主应力法或最大主应力投影法所确定的应力幅值和均值误差较大。
发明内容
为了解决现有技术中存在的以上技术问题,本申请提供一种应力幅值和均值的确定方法以及装置,能够准确地确定出应力幅值和均值,从而能够提高对结构疲劳强度的评估准确性。
为了实现上述目的,本申请提供的技术方案如下:
本申请实施例提供一种应力幅值和均值的确定方法,包括:
获取待评价对象在各个工况下的应力张量,并将每两个不同工况下的应力张量作差,得到各个应力张量差;
从所述各个应力张量差中确定最大应力张量差;其中,所述最大应力张量差由在第一工况下的应力张量和在第二工况下的应力张量进行作差得到;
根据所述最大应力张量差确定待评价对象的应力幅值,并根据所述在第一工况下的应力张量和所述在第二工况下的应力张量,确定待评价对象的应力均值。
可选的,所述从所述各个应力张量差中确定最大应力张量差,具体包括:
根据所述各个应力张量差,获取各个应力张量差的最大主应力;
根据所述各个应力张量差的最大主应力中的最大值,确定最大应力张量差。
可选的,所述根据所述最大应力张量差确定待评价对象的应力幅值,具体包括:
将所述最大应力张量差的最大主应力的二分之一作为待评价对象的应力幅值。
可选的,所述根据所述在第一工况下的应力张量和所述在第二工况下的应力张量,确定待评价对象的应力均值,具体包括:
根据所述在第一工况下的应力张量和所述在第二工况下的应力张量,得到目标张量和;
根据所述目标张量和,确定待评价对象的应力均值。
可选的,所述根据所述目标张量和,确定待评价对象的应力均值,具体包括:
将所述目标张量和在应力幅值方向上应力分量的二分之一作为待评价对象的应力均值;其中,所述应力幅值方向为最大应力张量差的最大主应力方向。
可选的,所述目标张量和在应力幅值方向上应力分量,具体为:
Figure GDA0002268004390000021
式中,
Figure GDA0002268004390000022
为目标张量和在应力幅值方向上分量;
Figure GDA0002268004390000023
为目标张量和在x轴方向的应力分量;
Figure GDA0002268004390000024
为目标张量和在y轴方向的应力分量;
Figure GDA0002268004390000025
为目标张量和在z轴方向的应力分量;
Figure GDA0002268004390000026
为目标张量和在x=y方向上的应力分量;
Figure GDA0002268004390000027
为目标张量和在y=z方向上的应力分量;
Figure GDA0002268004390000028
为目标张量和在z=x方向上的应力分量;l为最大应力张量差的最大主应力方向与x轴的夹角余弦值;m为最大应力张量差的最大主应力方向与y轴的夹角余弦值;n为最大应力张量差的最大主应力方向与z轴的夹角余弦值。
可选的,所述根据所述在第一工况下的应力张量和所述在第二工况下的应力张量,得到目标张量和,具体为:
σA+B=σAB
式中,σA+B为目标张量和;σA为在第一工况下的应力张量;σB为在第二工况下的应力张量。
可选的,所述获取待评价对象在各个工况下的应力张量,具体包括:
利用有限元计算方法,获取待评价对象在各个工况下的应力张量。
可选的,还包括:
根据所述待评价对象的应力幅值和应力平均值,绘制古德曼图,以便根据所述古德曼图进行疲劳强度的评价。
本申请实施例还提供了一种应力幅值和均值的确定装置,包括:
张量获取单元,用于获取待评价对象在各个工况下的应力张量,并将每两个不同工况下的应力张量作差,得到各个应力张量差;
第一确定单元,用于从所述各个应力张量差中确定最大应力张量差;其中,所述最大应力张量差由在第一工况下的应力张量和在第二工况下的应力张量进行作差得到;
第二确定单元,用于根据所述最大应力张量差确定待评价对象的应力幅值,并根据所述在第一工况下的应力张量和所述在第二工况下的应力张量,确定待评价对象的应力均值。
与现有技术相比,本申请实施例至少具有以下优点:
本申请实施例提供的应力幅值和均值的确定方法,在获取到待评价对象在各个工况下的应力张量后,先将每两个不同工况下的应力张量作差,得到各个应力张量差,再从各个应力张量差中确定最大应力张量差以及该最大应力张量差所对应的第一工况和第二工况,以便根据该最大应力张量差确定待评价对象的应力幅值,并根据在第一工况下的应力张量和在第二工况下的应力张量,确定待评价对象的应力均值。其中,由于每个应力张量差能够准确地表征在从一种工况过渡到另一种工况时所产生的应力变化,因而,当从所有应力张量差中筛选出最大应力张量差时,该最大应力张量差能够准确地表征待评价对象的最大应力变量,即应力幅值,。如此,根据最大应力张量差所确定的应力幅值能够准确地表征待评价对象的最大应力幅值,且根据该最大应力张量差对应的在第一工况下的应力张量和所述在第二工况下的应力张量所确定的应力均值能够准确地表征待评价对象的最大应力均值,使得待评价对象的应力幅值和应力均值更接近于实际应力幅值和实际应力均值,提高应力幅值和均值的准确性,从而能够提高对结构疲劳强度的评估准确性。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1为本申请实施例提供的应力幅值和均值的确定方法的流程图;
图2为本申请实施例提供的一个等效单元的疲劳强度评价对比图;
图3为本申请实施例提供的应力幅值和均值的确定装置的结构示意图。
具体实施方式
目前,通常使用最大主应力法和最大主应力投影法来确定应力幅值和均值。下面分别介绍这两种方法。
最大主应力法确定应力幅值和应力均值的具体过程为:获取待评价对象在所有工况下的应力最大值σmax与应力最小值σmin,以便根据公式(1)确定待评价对象的应力幅值,以及根据公式(2)确定待评价对象的应力均值。
Figure GDA0002268004390000041
Figure GDA0002268004390000042
式中,σa为待评价对象的应力幅值;σm为待评价对象的应力均值;σmax为待评价对象在所有工况下的应力最大值;σmin为待评价对象在所有工况下的应力最小值。
最大主应力投影法确定应力幅值和应力均值的具体过程为:首先,获取待评价对象在所有工况下的应力最大值σmax,并将σmax的应力方向作为应力幅值方向;其次,将待评价对象在所有工况下的应力向该应力幅值方向投影,并选择最小的投影值作为应力最小值σmin;最后,根据公式(1)确定待评价对象的应力幅值,以及根据公式(2)确定待评价对象的应力均值。
基于上述提供的两种确定应力幅值和均值的方法,发明人经过研究发现,上述两种确定应力幅值和均值的方法存在以下问题:
由于最大主应力法是将应力最大值σmax与应力最小值σmin作为标量进行加减的,但是在实际中应力最大值σmax的方向与应力最小值σmin的方向是不相同的,使得应力最大值σmax的方向与应力最小值σmin的方向之间存在夹角,因而,根据公式(1)和(2)所确定的应力幅值和应力均值是存在误差的,而且,该误差会随着应力最大值σmax的方向与应力最小值σmin的方向之间夹角的变化而发生改变,且当应力最大值σmax的方向与应力最小值σmin的方向之间夹角为90°时,该误差达到最大。例如,对于一个简支梁来说,载荷的变化使得梁的下表面的拉应力在100MPa和90MPa之间徘徊,由此可知该两个工况的应力幅值是5MPa;但是在利用最大主应力法确定应力幅值时,由于应力最大值为100MPa且应力最小值为0MPa,因而,利用最大主应力法确定的应力幅值为50MPa,如此使得利用最大主应力法确定的应力幅值高于实际应力幅值,此时,在根据基于最大主应力法确定的应力幅值进行疲劳强度评价时,所得的评价结果比较保守。
另外,由于最大主应力法和最大主应力投影法均是基于“最大幅值方向就是最大主应力方向”的假设来确定应力幅值和均值的,但是在实际应用中最大幅值方向与最大主应力方向重合的概率极低(例如,当压应力绝对值远大于拉应力时,最大幅值方向与最大主应力方向是不重合的),从而导致利用最大主应力法或最大主应力投影法所确定的应力幅值和均值误差较大。
此外,当利用投影法确定应力幅值和均值时,因在投影法中通常会忽略待评价对象主要应力特征,从而导致利用投影法所确定的应力幅值和均值误差较大,此时,在根据基于最大主应力投影法确定的应力幅值进行疲劳强度评价时,所得的评价结果比较冒险。
为了解决上述两种确定应力幅值和均值的方法所存在的技术问题,本申请实施例基于应力张量提供了一种的应力幅值和均值的确定方法,在该方法中,在获取到待评价对象在各个工况下的应力张量后,先将每两个不同工况下的应力张量作差,得到各个应力张量差,再从各个应力张量差中确定最大应力张量差以及该最大应力张量差所对应的第一工况和第二工况,以便根据该最大应力张量差确定待评价对象的应力幅值,并根据在第一工况下的应力张量和在第二工况下的应力张量,确定待评价对象的应力均值。其中,由于每个应力张量差能够准确地表征在从一种工况过渡到另一种工况时所产生的应力变化,因而,当从所有应力张量差中筛选出最大应力张量差时,该最大应力张量差能够准确地表征待评价对象的最大应力变量,即应力幅值,。如此,根据最大应力张量差所确定的应力幅值能够准确地表征待评价对象的最大应力幅值,且根据该最大应力张量差对应的在第一工况下的应力张量和所述在第二工况下的应力张量所确定的应力均值能够准确地表征待评价对象的最大应力均值,使得待评价对象的应力幅值和应力均值更接近于实际应力幅值和实际应力均值,提高应力幅值和均值的准确性,从而能够提高对结构疲劳强度的评估准确性。
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
方法实施例
参见图1,该图为本申请实施例提供的应力幅值和均值的确定方法的流程图。
本申请实施例提供的应力幅值和均值的确定方法,包括步骤S1-S5:
S1:获取待评价对象在各个工况下的应力张量。
待评价对象用于表示需要进行疲劳强度评价的对象;而且,本申请实施例不限定待评价对象的具体实施方式,例如,待评价对象可以包括车辆,也可以包括车辆中的至少一个部件,也可以包括在对结构进行有限元计算时所获得的至少一个等效单元。
需要说明的是,对于每个部件(或每个等效单元)来说,均可以采用相同的方法进行应力幅值和均值的确定。因而,为了便于理解和解释本申请实施例提供的应力幅值和均值的确定方法,后续步骤中均将以待评价对象包括车辆中的一个部件(或一个等效单元)为例进行说明。
工况是指车辆行驶状况以及其他类型状况,本申请对此不做具体限定。例如,工况可以包括:刹车、左转、直行等。
应力张量用于记录待评价对象在一种工况下的多轴应力状态数据或单轴应力状态数据。
作为一种实施方式,为了能够提高应力幅值和均值的准确性,本申请实施例还提供了步骤S1的一种具体实施方式,在该实施方式中,S1具体包括:利用有限元计算方法,获取待评价对象在各个工况下的应力张量。
本申请实施例不限定有限元计算方法的具体实施方式,可以采用任一种现有或未来出现的有限元计算方法来获取待评价对象在各个工况下的应力张量。
在本申请实施例中,先需要将结构(例如,车辆或部件等的结构)服役过程中的所有可能的工况,抽象成N种载荷和约束的组合;再基于该N种载荷和约束的组合对该结构进行有限元计算,以便得到待评价对象(例如,该待评价对象可以包括至少一个等效单元)在各个工况下的应力张量。其中,N为正整数,且N≥1。
以上为步骤S1的具体实施方式。
S2:将每两个不同工况下的应力张量作差,得到各个应力张量差。
应力张量差用于表示在从一种工况过渡到另一种工况时所产生的应力变化。
作为一种实施方式,当存在M种工况时,则S2具体为:利用公式(3),将每个工况下的应力张量与其他工况下的应力张量作差,得到M*(M-1)个应力张量差。
σi-j=σij   (3)
式中,σi-j表示在第i个工况下的应力张量与在第j个工况下的应力张量之间的应力张量差;σi为第i个工况下的应力张量;σj为在第j个工况下的应力张量。
在本申请实施例中,通常需要将每个工况下的应力张量与其他工况下的应力张量进行作差,以便获得从每一工况过渡到其他工况时的应力变化,以便后续能够根据每一工况过渡到其他工况时的应力变化,准确地确定出待评价对象的最大应力变量,即应力幅值,。
S3:从各个应力张量差中确定最大应力张量差。
其中,最大应力张量差由在第一工况下的应力张量和在第二工况下的应力张量进行作差得到。
在本申请实施例中,可以通过比较每个应力张量差的最大主应力来确定最大应力张量差,因而,本申请实施例提供了步骤S3的一种实施方式,在该实施方式中,S3具体包括步骤S31-S32:
S31:根据各个应力张量差,获取各个应力张量差的最大主应力。
本申请实施例不限定根据应力张量差确定该应力张量差的最大主应力的计算过程,可以采用现有或未来出现的任一种能够根据应力张量差确定该应力张量差的最大主应力的计算方法来执行步骤S31。
S32:根据各个应力张量差的最大主应力中的最大值,确定最大应力张量差。
在本申请实施例中,先从各个应力张量差的最大主应力中筛选出最大值,并将该最大值所对应的应力张量差作为最大应力张量差。
作为示例,当在步骤S31中得到第1个应力张量差的最大主应力至第T个应力张量差的最大主应力,且第Y个应力张量差的最大主应力大于除了第Y个应力张量差的最大主应力以外的其他任一应力张量差的最大主应力时,则步骤S32具体可以为:从第1个应力张量差的最大主应力、第2个应力张量差的最大主应力、……、和第T个应力张量差的最大主应力中进行最大值筛选,得到第Y个应力张量差的最大主应力,并将第Y个应力张量差确定为最大应力张量差。
以上为步骤S3的具体实施方式。
S4:根据最大应力张量差确定待评价对象的应力幅值。
在本申请实施例中,由于最大应力张量差能够准确地表征待评价对象的最大应力变量,即应力幅值,因而,根据最大应力张量差所确定的应力幅值能够准确地表征待评价对象的最大应力幅值,使得根据最大应力张量差所确定的应力幅值更接近于实际应力幅值,如此提高了应力幅值的准确性。
作为一种实施方式,步骤S4具体可以为:将最大应力张量差的最大主应力的二分之一作为待评价对象的应力幅值。
在本申请实施例中,应力幅值方向就是最大应力张量差的最大主应力的方向。
S5:根据在第一工况下的应力张量和在第二工况下的应力张量,确定待评价对象的应力均值。
在本申请实施例中,由于最大应力张量差是由在第一工况下的应力张量和在第二工况下的应力张量进行作差得到的,因而,根据在第一工况下的应力张量和所述在第二工况下的应力张量所确定的应力均值能够准确地表征待评价对象的最大应力均值,使得根据在第一工况下的应力张量和所述在第二工况下的应力张量所确定的应力均值更接近于实际应力均值,如此提高了应力均值的准确性。
作为一种实施方式,S5具体可以包括步骤S51-S52:
S51:根据在第一工况下的应力张量和所述在第二工况下的应力张量,得到目标张量和。
作为一种实施方式,S51具体可以为:根据在第一工况下的应力张量和所述在第二工况下的应力张量,利用公式(4),得到目标张量和。
σA+B=σAB   (4)
式中,σA+B为目标张量和;σA为在第一工况下的应力张量;σB为在第二工况下的应力张量。
S52:根据目标张量和,确定待评价对象的应力均值。
作为一种实施方式,可以根据目标张量和在应力幅值方向的应力分量确定待评价对象的应力均值,因而,本申请实施例提供了步骤S52的一种实施方式,在该实施方式中,步骤S52具体可以为:将目标张量和在应力幅值方向上应力分量的二分之一作为待评价对象的应力均值。
其中,应力幅值方向为最大应力张量差的最大主应力方向。
另外,目标张量和在应力幅值方向上应力分量,可以利用公式(5)确定,且公式(5)具体为:
Figure GDA0002268004390000091
式中,
Figure GDA0002268004390000101
为目标张量和在应力幅值方向上分量;
Figure GDA0002268004390000102
为目标张量和在x轴方向的应力分量;
Figure GDA0002268004390000103
为目标张量和在y轴方向的应力分量;
Figure GDA0002268004390000104
为目标张量和在z轴方向的应力分量;
Figure GDA0002268004390000105
为目标张量和在x=y方向上的应力分量;
Figure GDA0002268004390000106
为目标张量和在y=z方向上的应力分量;
Figure GDA0002268004390000107
为目标张量和在z=x方向上的应力分量;l为最大应力张量差的最大主应力方向(也就是,应力幅值方向)与x轴的夹角余弦值;m为最大应力张量差的最大主应力方向(也就是,应力幅值方向)与y轴的夹角余弦值;n为最大应力张量差的最大主应力方向(也就是,应力幅值方向)与z轴的夹角余弦值。
在本申请实施例中,目标张量和在应力幅值方向上应力分量
Figure GDA0002268004390000108
的二分之一能够准确地表征待评价对象的最大应力均值,使得确定的应力均值更接近于实际应力均值,如此提高了应力均值的准确性。
以上为步骤S5的具体实施方式。
需要说明的是,本申请实施例不限定步骤S4和S5的执行顺序,可以依次执行步骤S4和S5,也可以依次执行步骤S5和S4,还可以同时执行步骤S4和S5。
另外,在本申请实施例中,在确定了待评价对象的应力幅值和应力平均值之后,还可以根据待评价对象的应力幅值和应力平均值进行疲劳强度的评价。因而,本申请实施例还提供了应力幅值和均值的确定方法的另一种实施方式,在该实施方式中,应力幅值和均值的确定方法除了包括步骤S1-S5以外,还包括步骤S6:
S6:根据待评价对象的应力幅值和应力平均值,绘制古德曼图,以便根据所述古德曼图进行疲劳强度的评价。
古德曼图就是Goodman图,而且,古德曼图用于进行疲劳强度的评价。
为了便于理解和解释S6,下面结合示例进行说明。
作为示例,假设待评价对象包括通过对结构进行有限元计算所获得的至少一个等效单元时,则S6具体可以为:将至少一个等效单元的应力幅值和应力平均值绘制于Goodman图中,以便根据Goodman图进行疲劳强度的评价。
此外,为了能够证明本申请实施例提供的应力幅值和均值的确定方法的准确性,如图2所示,本申请实施例还提供了一个等效单元的疲劳强度评价对比图,该图能够有效地证明:最大主应力法确定的应力幅值和应力均值过于保守,而且,最大主应力投影法(最小主应力投影)以及最大主应力投影法(全应力投影)确定的应力幅值和应力均值均过于冒险,只有本申请实施例提供的应力幅值和均值的确定方法)所确定的应力幅值和均值更接近于实际的应力幅值和均值。由此可知,本申请实施例提供的应力幅值和均值的确定方法能够准确地确定应力幅值和应力均值,以便后续在根据该应力幅值和应力均值进行疲劳强度评价时能够准确地获得待评价对象的疲劳强度结果,提高了结构疲劳强度的评估准确性。
以上为方法实施例提供的应力幅值和均值的确定方法的具体实施方式,在该实施方式中,在获取到待评价对象在各个工况下的应力张量后,先将每两个不同工况下的应力张量作差,得到各个应力张量差,再从各个应力张量差中确定最大应力张量差以及该最大应力张量差所对应的第一工况和第二工况,以便根据该最大应力张量差确定待评价对象的应力幅值,并根据在第一工况下的应力张量和在第二工况下的应力张量,确定待评价对象的应力均值。其中,由于每个应力张量差能够准确地表征在从一种工况过渡到另一种工况时所产生的应力变化,因而,当从所有应力张量差中筛选出最大应力张量差时,该最大应力张量差能够准确地表征待评价对象的最大应力变量,即应力幅值,。如此,根据最大应力张量差所确定的应力幅值能够准确地表征待评价对象的最大应力幅值,且根据该最大应力张量差对应的在第一工况下的应力张量和所述在第二工况下的应力张量所确定的应力均值能够准确地表征待评价对象的最大应力均值,使得能够准确地确定应力幅值和均值,从而能够提高结构疲劳强度的评估准确性。
基于上述方法实施例提供的应力幅值和均值的确定方法,本申请实施例还提供了一种应力幅值和均值的确定装置,下面结合附图进行解释和说明。
装置实施例
装置实施例提供的应力幅值和均值的确定装置的技术详情请参照上述方法实施例。
参见图3,该图为本申请实施例提供的应力幅值和均值的确定装置的结构示意图。
本申请实施例提供的应力幅值和均值的确定装置30,包括:
张量获取单元31,用于获取待评价对象在各个工况下的应力张量,并将每两个不同工况下的应力张量作差,得到各个应力张量差;
第一确定单元32,用于从所述各个应力张量差中确定最大应力张量差;其中,所述最大应力张量差由在第一工况下的应力张量和在第二工况下的应力张量进行作差得到;
第二确定单元33,用于根据所述最大应力张量差确定待评价对象的应力幅值,并根据所述在第一工况下的应力张量和所述在第二工况下的应力张量,确定待评价对象的应力均值。
作为一种实施方式,为了提高应力幅值和均值的准确性,所述张量获取单元31,具体包括:
第一获取子单元,用于根据所述各个应力张量差,获取各个应力张量差的最大主应力;
第一确定子单元,用于根据所述各个应力张量差的最大主应力中的最大值,确定最大应力张量差。
作为一种实施方式,为了提高应力幅值和均值的准确性,所述第二确定单元33,具体包括:
第二确定子单元,用于将所述最大应力张量差的最大主应力的二分之一作为待评价对象的应力幅值。
作为一种实施方式,为了提高应力幅值和均值的准确性,所述第二确定单元33,具体包括:
第三确定子单元,用于根据所述在第一工况下的应力张量和所述在第二工况下的应力张量,得到目标张量和;
第四确定子单元,用于根据所述目标张量和,确定待评价对象的应力均值。
作为一种实施方式,为了提高应力幅值和均值的准确性,所述第四确定子单元,具体用于:
将所述目标张量和在应力幅值方向上应力分量的二分之一作为待评价对象的应力均值;其中,所述应力幅值方向为最大应力张量差的最大主应力方向。
作为一种实施方式,为了提高应力幅值和均值的准确性,所述目标张量和在应力幅值方向上应力分量,具体为:
Figure GDA0002268004390000131
式中,
Figure GDA0002268004390000132
为目标张量和在应力幅值方向上分量;
Figure GDA0002268004390000133
为目标张量和在x轴方向的应力分量;
Figure GDA0002268004390000134
为目标张量和在y轴方向的应力分量;
Figure GDA0002268004390000135
为目标张量和在z轴方向的应力分量;
Figure GDA0002268004390000136
为目标张量和在x=y方向上的应力分量;
Figure GDA0002268004390000137
为目标张量和在y=z方向上的应力分量;
Figure GDA0002268004390000138
为目标张量和在z=x方向上的应力分量;l为最大应力张量差的最大主应力方向与x轴的夹角余弦值;m为最大应力张量差的最大主应力方向与y轴的夹角余弦值;n为最大应力张量差的最大主应力方向与z轴的夹角余弦值。
作为一种实施方式,为了提高应力幅值和均值的准确性,所述第三确定子单元,具体为:
σA+B=σAB
式中,σA+B为目标张量和;σA为在第一工况下的应力张量;σB为在第二工况下的应力张量。
作为一种实施方式,为了提高应力幅值和均值的准确性,所述张量获取单元31,具体用于:
利用有限元计算方法,获取待评价对象在各个工况下的应力张量。
作为一种实施方式,为了提高应力幅值和均值的准确性,所述装置30还包括:
绘制单元,用于根据所述待评价对象的应力幅值和应力平均值,绘制古德曼图,以便根据所述古德曼图进行疲劳强度的评价。
以上为本申请实施例提供的应力幅值和均值的确定装置30的具体实施方式,在该实施方式中,在获取到待评价对象在各个工况下的应力张量后,先将每两个不同工况下的应力张量作差,得到各个应力张量差,再从各个应力张量差中确定最大应力张量差以及该最大应力张量差所对应的第一工况和第二工况,以便根据该最大应力张量差确定待评价对象的应力幅值,并根据在第一工况下的应力张量和在第二工况下的应力张量,确定待评价对象的应力均值。其中,由于每个应力张量差能够准确地表征在从一种工况过渡到另一种工况时所产生的应力变化,因而,当从所有应力张量差中筛选出最大应力张量差时,该最大应力张量差能够准确地表征待评价对象的最大应力变量,即应力幅值,。如此,根据最大应力张量差所确定的应力幅值能够准确地表征待评价对象的最大应力幅值,且根据该最大应力张量差对应的在第一工况下的应力张量和所述在第二工况下的应力张量所确定的应力均值能够准确地表征待评价对象的最大应力均值,使得待评价对象的应力幅值和应力均值更接近于实际应力幅值和实际应力均值,提高应力幅值和均值的准确性,从而能够提高结构疲劳强度的评估准确性。
应当理解,在本申请中,“至少一个(项)”是指一个或者多个,“多个”是指两个或两个以上。“和/或”,用于描述关联对象的关联关系,表示可以存在三种关系,例如,“A和/或B”可以表示:只存在A,只存在B以及同时存在A和B三种情况,其中A,B可以是单数或者复数。字符“/”一般表示前后关联对象是一种“或”的关系。“以下至少一项(个)”或其类似表达,是指这些项中的任意组合,包括单项(个)或复数项(个)的任意组合。例如,a,b或c中的至少一项(个),可以表示:a,b,c,“a和b”,“a和c”,“b和c”,或“a和b和c”,其中a,b,c可以是单个,也可以是多个。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制。虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明。任何熟悉本领域的技术人员,在不脱离本发明技术方案范围情况下,都可利用上述揭示的方法和技术内容对本发明技术方案做出许多可能的变动和修饰,或修改为等同变化的等效实施例。因此,凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所做的任何简单修改、等同变化及修饰,均仍属于本发明技术方案保护的范围内。

Claims (5)

1.一种应力幅值和均值的确定方法,其特征在于,包括:
获取待评价对象在各个工况下的应力张量,并将每两个不同工况下的应力张量作差,得到各个应力张量差;
从所述各个应力张量差中确定最大应力张量差;其中,所述最大应力张量差由在第一工况下的应力张量和在第二工况下的应力张量进行作差得到;所述最大应力张量差是根据所述各个应力张量差的最大主应力中的最大值所确定的;
根据所述最大应力张量差确定待评价对象的应力幅值,并根据所述在第一工况下的应力张量和所述在第二工况下的应力张量,确定待评价对象的应力均值;所述应力幅值为所述最大应力张量差的最大主应力的二分之一;所述应力均值为目标张量和在应力幅值方向上应力分量的二分之一;所述目标张量和是根据所述在第一工况下的应力张量和所述在第二工况下的应力张量所确定的;所述应力幅值方向为最大应力张量差的最大主应力方向;
根据所述待评价对象的应力幅值和应力平均值,绘制古德曼图,以便根据所述古德曼图进行疲劳强度的评价。
2.根据权利要求1所述的方法,其特征在于,所述目标张量和在应力幅值方向上应力分量,具体为:
Figure FDA0004128749020000011
式中,
Figure FDA0004128749020000012
为目标张量和在应力幅值方向上分量;
Figure FDA0004128749020000013
为目标张量和在x轴方向的应力分量;
Figure FDA0004128749020000014
为目标张量和在y轴方向的应力分量;
Figure FDA0004128749020000015
为目标张量和在z轴方向的应力分量;
Figure FDA0004128749020000016
为目标张量和在x=y方向上的应力分量;
Figure FDA0004128749020000017
为目标张量和在y=z方向上的应力分量;
Figure FDA0004128749020000018
为目标张量和在z=x方向上的应力分量;l为最大应力张量差的最大主应力方向与x轴的夹角余弦值;m为最大应力张量差的最大主应力方向与y轴的夹角余弦值;n为最大应力张量差的最大主应力方向与z轴的夹角余弦值。
3.根据权利要求1所述的方法,其特征在于,所述目标张量和的确定过程,具体为:
σA+B=σAB
式中,σA+B为目标张量和;σA为在第一工况下的应力张量;σB为在第二工况下的应力张量。
4.根据权利要求1所述的方法,其特征在于,所述获取待评价对象在各个工况下的应力张量,具体包括:
利用有限元计算方法,获取待评价对象在各个工况下的应力张量。
5.一种应力幅值和均值的确定装置,其特征在于,包括:
张量获取单元,用于获取待评价对象在各个工况下的应力张量,并将每两个不同工况下的应力张量作差,得到各个应力张量差;
第一确定单元,用于从所述各个应力张量差中确定最大应力张量差;其中,所述最大应力张量差由在第一工况下的应力张量和在第二工况下的应力张量进行作差得到;所述最大应力张量差是根据所述各个应力张量差的最大主应力中的最大值所确定的;
第二确定单元,用于根据所述最大应力张量差确定待评价对象的应力幅值,并根据所述在第一工况下的应力张量和所述在第二工况下的应力张量,确定待评价对象的应力均值;所述应力幅值为所述最大应力张量差的最大主应力的二分之一;所述应力均值为目标张量和在应力幅值方向上应力分量的二分之一;所述目标张量和是根据所述在第一工况下的应力张量和所述在第二工况下的应力张量所确定的;所述应力幅值方向为最大应力张量差的最大主应力方向;
绘制单元,用于根据所述待评价对象的应力幅值和应力平均值,绘制古德曼图,以便根据所述古德曼图进行疲劳强度的评价。
CN201910780182.1A 2019-08-22 2019-08-22 一种应力幅值和均值的确定方法以及装置 Active CN110705137B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910780182.1A CN110705137B (zh) 2019-08-22 2019-08-22 一种应力幅值和均值的确定方法以及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910780182.1A CN110705137B (zh) 2019-08-22 2019-08-22 一种应力幅值和均值的确定方法以及装置

Publications (2)

Publication Number Publication Date
CN110705137A CN110705137A (zh) 2020-01-17
CN110705137B true CN110705137B (zh) 2023-04-21

Family

ID=69193725

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910780182.1A Active CN110705137B (zh) 2019-08-22 2019-08-22 一种应力幅值和均值的确定方法以及装置

Country Status (1)

Country Link
CN (1) CN110705137B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110727989B (zh) * 2018-06-28 2023-08-04 中车大同电力机车有限公司 结构疲劳强度分析方法、装置以及计算机可读存储介质
CN111695214B (zh) * 2020-05-26 2023-09-08 上海拜安传感技术有限公司 基于统计模型的岸桥起重机疲劳损伤的确定方法
CN111950101A (zh) * 2020-08-07 2020-11-17 中铁轨道交通装备有限公司 一种基于张量的多轴疲劳分析方法

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20030059425A (ko) * 2001-12-29 2003-07-10 삼성테크윈 주식회사 고주기 피로 수명 해석 방법
JP2011108032A (ja) * 2009-11-18 2011-06-02 Ihi Corp 多結晶金属材の疲労強度評価装置と方法
JP2013044667A (ja) * 2011-08-25 2013-03-04 Ihi Corp 多軸疲労寿命評価方法
JP2013044666A (ja) * 2011-08-25 2013-03-04 Ihi Corp 多軸疲労寿命評価方法
CN103604688A (zh) * 2013-12-01 2014-02-26 北京航空航天大学 一种基于临界面法的塑性金属材料多轴高周疲劳寿命预测方法
CN104392130A (zh) * 2014-11-21 2015-03-04 南京衍达软件科技有限公司 多轴疲劳的最损伤载荷方向的确定方法及其应用
CN106021688A (zh) * 2016-05-13 2016-10-12 中车青岛四方机车车辆股份有限公司 一种轨道交通机车的车轮强度的评估方法以及装置
WO2017017343A1 (fr) * 2015-07-27 2017-02-02 Safran Aircraft Engines Procédé de détermination d'une contrainte alternée maximale admissible pour une pièce soumise a des sollicitations cycliques; unité de détermination d'une telle contrainte
CN106951582A (zh) * 2017-01-24 2017-07-14 浙江工业大学 一种c型独立液货罐y型接头疲劳评估筛分准则的确定方法
CN108733879A (zh) * 2018-04-12 2018-11-02 南京航空航天大学 一种多轴非比例载荷下剪应力幅值计算方法
CN108959692A (zh) * 2018-04-28 2018-12-07 南京航空航天大学 一种多轴非比例载荷下平均剪应力计算方法
CN109190177A (zh) * 2018-08-06 2019-01-11 同济大学 一种轨道车辆车体焊缝疲劳可靠性评估方法及系统
JP2019008630A (ja) * 2017-06-27 2019-01-17 国立大学法人大阪大学 破壊現象の解析装置、破壊現象の解析方法、及びプログラム
CN109543355A (zh) * 2019-01-03 2019-03-29 付远 一种基于有限元法的应力-循环次数(s-n)的疲劳失效算法

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20030059425A (ko) * 2001-12-29 2003-07-10 삼성테크윈 주식회사 고주기 피로 수명 해석 방법
JP2011108032A (ja) * 2009-11-18 2011-06-02 Ihi Corp 多結晶金属材の疲労強度評価装置と方法
JP2013044667A (ja) * 2011-08-25 2013-03-04 Ihi Corp 多軸疲労寿命評価方法
JP2013044666A (ja) * 2011-08-25 2013-03-04 Ihi Corp 多軸疲労寿命評価方法
CN103604688A (zh) * 2013-12-01 2014-02-26 北京航空航天大学 一种基于临界面法的塑性金属材料多轴高周疲劳寿命预测方法
CN104392130A (zh) * 2014-11-21 2015-03-04 南京衍达软件科技有限公司 多轴疲劳的最损伤载荷方向的确定方法及其应用
WO2017017343A1 (fr) * 2015-07-27 2017-02-02 Safran Aircraft Engines Procédé de détermination d'une contrainte alternée maximale admissible pour une pièce soumise a des sollicitations cycliques; unité de détermination d'une telle contrainte
CN106021688A (zh) * 2016-05-13 2016-10-12 中车青岛四方机车车辆股份有限公司 一种轨道交通机车的车轮强度的评估方法以及装置
CN106951582A (zh) * 2017-01-24 2017-07-14 浙江工业大学 一种c型独立液货罐y型接头疲劳评估筛分准则的确定方法
JP2019008630A (ja) * 2017-06-27 2019-01-17 国立大学法人大阪大学 破壊現象の解析装置、破壊現象の解析方法、及びプログラム
CN108733879A (zh) * 2018-04-12 2018-11-02 南京航空航天大学 一种多轴非比例载荷下剪应力幅值计算方法
CN108959692A (zh) * 2018-04-28 2018-12-07 南京航空航天大学 一种多轴非比例载荷下平均剪应力计算方法
CN109190177A (zh) * 2018-08-06 2019-01-11 同济大学 一种轨道车辆车体焊缝疲劳可靠性评估方法及系统
CN109543355A (zh) * 2019-01-03 2019-03-29 付远 一种基于有限元法的应力-循环次数(s-n)的疲劳失效算法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Multiaxial fatigue criterion for a high-density polyethylene thermoplastic;Berrehili, A等;《FATIGUE & FRACTURE OF ENGINEERING MATERIALS & STRUCTURES》;20100630;第33卷(第6期);第345-357页 *
Validating Fatigue Safety Factor Calculation Methods for Cardiovascular Stents;Marrey, Ramesh等;《JOURNAL OF BIOMECHANICAL ENGINEERING-TRANSACTIONS OF THE ASME》;20180630;第140卷(第6期);第1-32页 *
多轴随机载荷下疲劳寿命预测方法的研究;宋宇;《中国优秀硕士学位论文全文数据库多轴随机载荷下疲劳寿命预测方法的研究》;20150531;第C029-28页 *
考虑材料非线性和辐辋过盈装配的车轮径向疲劳特性研究;王海霞等;《汽车工程》;20131231;第35卷(第9期);第822-826页 *
随机载荷下金属结构多轴疲劳寿命预测方法的研究;葛敬冉;《中国博士学位论文全文数据库基础科学辑》;20170228;第A004-20页 *

Also Published As

Publication number Publication date
CN110705137A (zh) 2020-01-17

Similar Documents

Publication Publication Date Title
CN110705137B (zh) 一种应力幅值和均值的确定方法以及装置
CN110442973B (zh) 一种车辆关键零部件的耐久度测试方法、系统及存储介质
JP5576567B2 (ja) 異常の発生を検知する方法、装置及びコンピュータプログラム
CN117076935B (zh) 数字孪生辅助的机械故障数据轻量级生成方法及系统
CN113554624B (zh) 异常检测方法、装置及计算机存储介质
CN117091754B (zh) 基于数据分析的大型设备故障检测方法及系统
CN112528388B (zh) 一种悬置支架强度分析方法、装置、存储介质及终端设备
CN112229654B (zh) 板簧试验方法、装置、电子设备和计算机可读介质
CN115183963A (zh) 桥梁异常状态检测方法、系统、设备及介质
CN114759904A (zh) 数据处理方法、装置、设备、可读存储介质及程序产品
US11416371B2 (en) Method and apparatus for evaluating and selecting signal comparison metrics
Li et al. Multichannel intelligent fault diagnosis of hoisting system using differential search algorithm‐variational mode decomposition and improved deep convolutional neural network
CN111639664A (zh) 一种基于多策略模式的线条痕迹批量比对系统
CN104200002A (zh) 一种粘性阻尼振动信号中的模态参数提取方法
US11733302B2 (en) Complicated system fault diagnosis method and system based on multi-stage model
CN117629636B (zh) 一种燃机滚动轴承的健康评估、故障诊断方法及系统
JP7272143B2 (ja) 検査機管理システム及び検査機管理方法
CN117526318B (zh) 配电网状态评估方法、装置、电子设备、存储介质
CN115905372B (zh) 一种用于核电配电柜的抗震结构分析方法及系统
CN114297782A (zh) 载荷谱处理方法、装置、设备及存储介质
CN110598768B (zh) 一种齿轮故障分类方法、分类装置及可读存储介质
Zhuang et al. Improved Wavelet Threshold Denoising and Multi-scale Principal Component Analysis Method with Application to TE Process Monitoring
KR20240057269A (ko) 불균형 데이터에 대한 전처리 방법 및 이를 이용한 최적 학습 모델 결정 시스템
CN116541768A (zh) 故障诊断方法、装置、设备及计算机可读存储介质
CN113901646A (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