CN112100563B - 多轴载荷等效处理方法、装置、计算机设备和介质 - Google Patents
多轴载荷等效处理方法、装置、计算机设备和介质 Download PDFInfo
- Publication number
- CN112100563B CN112100563B CN202010953720.5A CN202010953720A CN112100563B CN 112100563 B CN112100563 B CN 112100563B CN 202010953720 A CN202010953720 A CN 202010953720A CN 112100563 B CN112100563 B CN 112100563B
- Authority
- CN
- China
- Prior art keywords
- load
- equivalent
- loads
- decomposed
- damage
- 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
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 26
- 230000006378 damage Effects 0.000 claims abstract description 166
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 59
- 238000012360 testing method Methods 0.000 claims abstract description 52
- 238000000034 method Methods 0.000 claims abstract description 45
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 238000012545 processing Methods 0.000 claims abstract description 18
- 239000000463 material Substances 0.000 claims description 19
- 238000004590 computer program Methods 0.000 claims description 15
- 238000009825 accumulation Methods 0.000 claims description 11
- 238000007781 pre-processing Methods 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 3
- 230000006735 deficit Effects 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 description 15
- 208000027418 Wounds and injury Diseases 0.000 description 11
- 208000014674 injury Diseases 0.000 description 11
- 230000008569 process Effects 0.000 description 11
- 239000013598 vector Substances 0.000 description 7
- 230000000694 effects Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012942 design verification Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- 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)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种多轴载荷等效处理方法、装置、计算机设备和介质，以提高试验台架加载的等效载荷的等效性、准确性。方法部分包括：将多轴载荷投影到载荷空间中，以得到若干个方向上的分解载荷；分别对若干个方向上的分解载荷进行伪损伤计算，以获取分解载荷的伪损伤；根据分解载荷的伪损伤和多轴载荷循环次数确定分解载荷的总伪损伤；初始化等效载荷循环次数，并构建与若干个方向的分解载荷对应的预设信号形式的等效载荷式；根据等效载荷循环次数和等效载荷式确定等效载荷对应的总伪损伤式；令若干个方向的分解载荷的总伪损伤分别等于等效载荷对应的总伪损伤式，以获取载荷等效关系式并进行迭代求解以获取等效载荷式的参数；根据参数得到等效载荷式对应的等效载荷。
Description
技术领域
本发明涉及汽车耐久试验技术领域，尤其涉及一种多轴载荷等效处理方法、装置、计算机设备和介质。
背景技术
底盘等零部件耐久性能开发一般基于试验场耐久试验获取，而零部件所受的载荷谱是一种随机的，为了加快零件设计验证试验。需要将随机载荷谱转化为由若干等幅正弦波载荷谱组成的块谱进行加载，因为块谱有加载简单、试验时间短等特点。这种通过某种方法将随机谱转化为块谱的过程称为载荷等效。值得注意的是，载荷等效的等效性主要体现在两点，一是随机载荷与块谱载荷加载下零部件的疲劳寿命要相等；二是两者载荷加载下零部件的失效形式相同，即疲劳失效的位置和机理相同。
一般来说，载荷等效方法是基于单轴载荷进行的。然而，底盘的零部件都工作在多轴载荷状态下。因此，现有技术中常用单轴载荷分析方法将随机载荷转化为块谱时，会将多轴载荷分成多个单轴载荷进行单轴载荷等效，从而得到各单轴载荷转化若干正弦波载荷谱。
但是上述现有技术中，是将各轴载荷作为单轴载荷看作是相互独立的，其没有考虑各轴载荷之间的相位关系，对于底盘零件来说，大部分零部件都无法对各方向载荷解耦，所以对载荷之间的相位关系比较敏感，例如，转向节、轴节等。可见，单轴载荷等效方法只能处理单轴问题。例如，对于多轴问题，该方法只能将纵向、侧向、垂向3各方向的载荷分别看作是三个独立的单轴问题进行载荷等效。但实际上，这3个方向的载荷并不是完全独立的，其所受的多轴载荷是存在一定的相位关系，这种利用单轴等效方法解决多轴问题的解决方式，不能有效地解决受力复杂的底盘零件的载荷等效问题。因此，现有技术中，运用单轴载荷等效分析方法最终会影响试验台架加载的等效载荷的等效性。
发明内容
本发明提供一种多轴载荷等效处理方法、装置、计算机设备和介质，以解决现有技术中，会影响试验台架加载的等效载荷的等效性的问题。
一种多轴载荷等效处理方法，包括：
采集目标零部件在试验路面实际所受的多轴载荷，并确定多轴载荷循环次数；
将所述多轴载荷投影到载荷空间中，以得到若干个方向上的分解载荷；
分别对所述若干个方向上的分解载荷进行伪损伤计算，以获取所述分解载荷的伪损伤；
根据所述分解载荷的伪损伤和多轴载荷循环次数确定所述分解载荷的总伪损伤；
初始化等效载荷循环次数，并构建与所述若干个方向的分解载荷对应的预设信号形式的等效载荷式；
根据所述等效载荷循环次数和所述等效载荷式确定所述等效载荷对应的总伪损伤式；
令所述若干个方向的分解载荷的总伪损伤分别等于所述等效载荷对应的总伪损伤式，以获取载荷等效关系式；
对所述载荷等效关系式进行迭代求解，以获取所述等效载荷式的参数；
根据所述参数得到所述等效载荷式对应的等效载荷。
进一步地，所述将所述多轴载荷投影到载荷空间中，以得到若干个方向上的分解载荷，包括：
确定所述多轴载荷中不共线载荷的数量；
根据所述多轴载荷中不共线载荷的数量确定投影方向数量和投影方向；
按照所述投影方向数量和投影方向将所述多轴载荷投影到载荷空间中，以得到投影方向数量对应的若干个所述投影方向上的分解载荷。
进一步地，所述分别对所述若干个方向上的分解载荷进行伪损伤计算，以获取所述分解载荷的伪损伤，包括：
利用雨流计数法对所述分解载荷进行雨流计数处理，以获取所述分解载荷对应的多个应力循环；
根据所述目标零部件材料对应的S-N曲线计算各个所述应力循环对应的伪损伤值；
利用损伤累积理论对各个所述应力循环对应的伪损伤值进行累积处理，以获取所述分解载荷的伪损伤。
进一步地，所述各个应力循环对应的伪损伤值通过下式计算：
表示所述应力循环，/>表示伪损伤值，N表示目标零部件的疲劳寿命，b为对应目标零部件对应材料的S-N曲线的斜率，S0为目标零部件对应材料的S-N曲线的斜截距；
所述分解载荷的伪损伤通过下式计算：
表示所述分解载荷的伪损伤。
进一步地，所述构建与所述若干个方向的分解载荷对应的预设信号形式的等效载荷式，包括：
构建与所述若干个方向的分解载荷对应的正弦信号形式的等效载荷式。
进一步地，所述对所述载荷等效关系式进行迭代求解，以获取所述等效载荷式的参数，包括：
通过最小二乘法对所述载荷等效关系式进行迭代求解，以获取所述等效载荷式的参数，其中，各所述方向上等效载荷的总伪损伤与对应的分解载荷的总伪损伤之差的平方和取最小值为所述迭代求解的目标。
进一步地，所述将所述多轴载荷投影到载荷空间中，以得到若干个方向上的分解载荷之前，所述方法还包括：
对所述多轴载荷进行预处理，所述预处理包括去除毛刺、均值修正和滤波处理；
所述将所述多轴载荷投影到载荷空间中，以得到若干个方向上的分解载荷，包括：
将经过所述预处理的所述多轴载荷投影到所述载荷空间中，以得到所述若干个方向上的分解载荷。
一种多轴载荷等效处理装置，包括：
采集模块，用于采集目标零部件在试验路面实际所受的多轴载荷，并确定多轴载荷循环次数；
投影模块，用于将所述多轴载荷投影到载荷空间中，以得到若干个方向上的分解载荷；
计算模块，用于分别对所述若干个方向上的分解载荷进行伪损伤计算，以获取所述分解载荷的伪损伤；
第一确定模块，用于根据所述分解载荷的伪损伤和多轴载荷循环次数确定所述分解载荷的总伪损伤；
构建模块，用于初始化等效载荷循环次数，并构建与所述若干个方向的分解载荷对应的预设信号形式的等效载荷式；
第二确定模块，用于根据所述等效载荷循环次数和所述等效载荷式确定所述等效载荷对应的总伪损伤式；
第一获取模块，用于令所述若干个方向的分解载荷的总伪损伤分别等于所述等效载荷对应的总伪损伤式，以获取载荷等效关系式；
第二获取模块，用于对所述载荷等效关系式进行迭代求解，以获取所述等效载荷式的参数；
第三获取模块，用于根据所述参数得到所述等效载荷式对应的等效载荷。
一种计算机设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现前述多轴载荷等效处理方法的步骤。
一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现如前述多轴载荷等效处理方法的步骤。
本发明所提供的其中一个方案中，提供了一种多轴载荷等效处理方法，通过投影法对多轴载荷进行拆解成多个分解载荷，将多轴问题转化为单个载荷问题，并分别对若干个方向上的分解载荷进行伪损伤计算，并获取所述分解载荷的伪损伤；最后依据各个分解载荷的伪损伤构建相应载荷等效关系式，并对载荷等效关系式进行迭代求解，以获取所述等效载荷式的参数；最后根据所述参数得到所述等效载荷式对应的等效载荷，从而把复杂的多轴随机的多轴载荷等效地转化为等效载荷进行台架试验，由于保留了分解载荷是由多轴载荷分解得到，且通过各个分解载荷计算出的等效载荷之间保留多轴载荷的相位关系，从而提高了等效载荷的等效性，也就提高了后续零部件耐久试验的等效载荷的准确性。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对本发明实施例的描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例中多轴载荷等效处理方法的一个流程示意图；
图2是本发明实施例中多轴载荷的一个投影示意图；
图3是本发明实施例中多轴载荷的另一个投影示意图；
图4是本发明实施例中多轴载荷等效处理装置的一个结构示意图；
图5是本发明实施例中计算机设备的一个结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
在描述本发明所提供的多轴载荷等效处理方法之前，为了便于理解和描述本发明实施例，首先对本发明所涉及的概念、术语相关以及相关关系先进行介绍。
1、关于伪损伤
在计算汽车的某个零部件的疲劳损伤，以用于确定该零部件的耐久性能时，一般会利用有限元法得到零部件所受载荷下的应力或应变时间历程，再根据零部件材料的应力-寿命曲线(S-N曲线)或应变-寿命曲线(E-N曲线)，用损伤累积理论计算得到零部件的损伤值，此方法得出的损伤值通常称为真实损伤。也就是说真实损伤是载荷在零部件的作用效果，与载荷和零部件相关。其中，S-N曲线是以材料标准试件疲劳强度为纵坐标，以疲劳寿命N为横坐标，表示一定循环特征下标准试件的疲劳强度与疲劳寿命之间关系的曲线；E-N曲线是以材料标准试件应变幅度为纵坐标，以疲劳寿命N为横坐标，表示一定循环特征下标准试件的应变幅度与疲劳寿命之间关系的曲线。
与上述真实损伤不同，伪损伤不考虑具体的零部件，而是直接把任意施加于零部件的载荷信号，例如力、力矩、位移、加速度等看作是广义应力σG。并将该广义应力σG拆解为若干个应力循环每一个应力循环/>都根据对应的S-N曲线计算对应的损伤值/>再利用损伤累积理论计算得到该广义应力σG对应的总损伤值DG。这种对广义应力进行损伤累积得到的损伤值称为伪损伤。
其中，示例性的，单个应力循环产生的伪损伤/>可通过如下方式计算：
N表示疲劳寿命，b为对应材料的S-N曲线的斜率，S0为对应材料的S-N曲线的斜截距。
依据上式(1)，对广义应力σG拆解出的若干个应力循环对应的伪损伤/>进行叠加，便可得广义应力σG对应的总伪损伤为DG为：
2、关于伪损伤数
可以理解，上式中S0是与S-N曲线中的参数，为进一步使伪损伤与材料参数无关，本发明中定义dG为伪损伤数：
则依据上式，便可以得到广义应力σG对应的伪损伤与伪损伤数的关系为：
可以理解，由于伪损伤的绝对数值无实质性意义，作用是在于比较两组不同的载荷对零部件产生的疲劳损伤的程度，因此，通常以两组载荷对应的伪损伤之比，即伪损伤比来表述这种程度关系。例如，当伪损伤比接近1时，则可以认为两组载荷的强度相近，对零部件将产生相近的疲劳损伤。
假设有存在两组广义应力循环和/>可得两者损伤比为：
观察上式便可发现，S-N曲线的截距S0与损伤比无关，伪损伤之比即为伪损伤数之比。而斜率b则直接影响损伤比。
需要说明的是，令一般地，β取值在1至10之间。根据经验，研究零部件对象为焊接零件或者裂纹扩展现象时，β＝3；研究零部件对象为典型的汽车零件时(如零件表面相对粗糙的零件)，β＝5；研究零部件对象为具有光滑表面的零件时，β＝7。上述β取值在此仅为示例性说明，并不造成限定。
3、关于等效载荷和实测载荷
如背景所述，在进行零部件的耐久试验时，需用等效载荷代替随机交变载荷(也即实测得到的多轴载荷)用于加载试验，需让等效载荷产生的伪损伤Deq要等于多轴载荷的伪损伤DTarget。其中，多轴载荷为整车耐久性试验采集而来的载荷谱所获取，简而言之，多轴载荷是零部件实际受到的随机交变的多个载荷；等效载荷是用于某种固定形式，用于替代多轴载荷的载荷，等效载荷用于后续耐性试验时加载使用。值得说明的是，一个完整的耐久性试验包含多种路面，对应多种路况，不同路况下所受载荷情况有所不同，每种路面对应的载荷的循环次数也不尽相同。
为了便于理解后续的多轴载荷等效处理过程，这里先以单轴载荷为例说明下载荷循环次数、总伪损伤数、对应等效载荷的关系，单轴载荷的总伪损伤还与单轴载荷的循环次数关联。
设单轴载荷的总伪损伤数为：
nk为第k路面对应的循环次数；
dk为第k路面的伪损伤数。
该单轴载荷的等效载荷的伪损伤数为：
等效载荷的伪损伤要等于该单轴载荷的等效载荷的总伪损伤：
所以等效台架试验的等效载荷的应该为：
上述对本发明涉及到的伪损伤和伪损伤数等术语或概念进行了简单的介绍，便于说明后续的多轴载荷等效关系式以及对应的多轴载荷等效处理方法，在下文提到的多轴载荷等效处理方法将利用到上述所提及的术语或概念。在一实施例中，如图2所示，提供了一种多轴载荷等效处理方法，以完整的实施例进行介绍，包括如下步骤：
S10：采集目标零部件在试验路面实际所受的多轴载荷，并确定多轴载荷循环次数。
需要说明的是，汽车的诸如底盘等零部件所受的载荷是一种随机变化的载荷，并且并不是单轴、单一载荷，许多零部件的受载荷状态是处于多轴载荷力状态下，也就是收到的是多个方向的多轴载荷，因此，针对多轴载荷的零部件，本发明中称为目标零部件，例如底盘。在本发明中，可以获取利用试验场进行测试得到的目标零部件的载荷谱进行采集，从而得到汽车的目标零部件在试验路面实际所受的多轴载荷，并多轴载荷循环次数。
其中，试验路面通常为多个，不同试验路面的轴载荷循环次数可能会不一样。因此，为了后续的等效处理，本发明需获取不同试验路面的多轴载荷以及对应试验路面上的多轴载荷循环次数，便于确定后续不同试验路面上的多轴载荷的总伪损伤。
S20：将所述多轴载荷投影到载荷空间中，以得到若干个方向上的分解载荷。
汽车的底盘等目标零部件在试验路面上经常会承受多个方向的多轴载荷，本发明中，在对多轴载荷进行等效时，可利用投影法将单轴载荷的情况拓展到多轴载荷的情况，即将利用投影法将多轴载荷等效分解为多个单轴载荷进行等效，具体地，在获取到利用试验场试验并采集得到的多轴载荷之后，采样投影法将多轴载荷投影到载荷空间中，以得到若干个方向上的分解载荷。下面对将所述多轴载荷投影到载荷空间中，以得到若干个方向上的分解载荷的过程做进一步地解释和说明。
其中，在一实施方式中，将所述多轴载荷投影到载荷空间中，以得到若干个方向上的分解载荷，具体包括如下步骤：
S21：确定所述多轴载荷中不共线载荷的数量。
S22：根据所述多轴载荷中不共线载荷的数量确定投影方向数量和投影方向。
在采集目标零部件实际所受的多轴载荷的载荷图谱之后，可确定所述多轴载荷中不共线载荷的数量，根据所述多轴载荷中不共线载荷的数量确定投影方向数量和投影方向，其中，投影方向数量大于不共线载荷的数量。
S23：按照所述投影方向数量和投影方向将所述多轴载荷投影到载荷空间中，以得到投影方向数量对应的若干个所述投影方向上的分解载荷。
在得到投影方向数量和投影方向之后，按照投影方向数量和投影方向将多轴载荷投影到载荷空间中，以得到投影方向数量对应的若干个方向上的分解载荷。
需要说明的是，示例性的，根据经验获得如下投影方式较佳：如图2，如果多轴载荷是两个不共线的载荷，如图2中的不共线载荷L1(t)和L2(t)，则可将载荷空间分解为12个投影方向数量，并且在具体实施中，可在同个半圆面上均匀取点作为投影方向从而确定12个方向，随后依据分解出的12个方向进行投影；如果多轴载荷是三个不共线的载荷，如图3中所示的不共线载荷L1(t)、L2(t)和L3(t)，则可将载荷空间分解为25个投影方向数量，对应的25个投影方向设在同个半球面，并且在具体实施中，可在同个半球面上按照规范取25个点，从而确定出25个方向，最后依据分解的25个方向进行投影。需要说明的是，上述例子仅为示例性说明，本发明不做限制。
可以理解，在单轴载荷情况下，对于零部件上的某点x，在某时刻t，在单轴载荷L(t)的作用下，该点x始终处于单轴应力状态，且有σij(x,t)＝cij(x)·L(t)，其中，σij(x,t)表示零部件点x在单轴载荷L(t)下的应力循环，cij(x)表示该单轴载荷L(t)对应的载荷系数。
在多轴载荷情况下，某点x的应力循环σij(x,t)是多轴载荷Lk(t)的线性组合，也就是多轴载荷可拆解出的多个分解载荷的线性组合，也即存在如下关系式：
其中：
为应力循环σij(x,t)的第k个载荷系数；
L1(t)、2(t)、...、Ln(t)为多轴载荷分解出的时刻t的n个分解载荷；
cij(x)为系数向量，是多轴载荷L(t)到点x应力循环σij(x,t)的映射；
L(t)为时刻t的分解载荷向量。
可见，是随点x的位置变化而变化，与单轴载荷的情况类似，在实际开发过程中，希望用等效载荷来模拟耐久性试验的效果，因此，当零部件承受多个方向的多轴载荷时，等效载荷也需要从多个方向加载到零部件上，从而得到等效的目的。
由于希望等效载荷与试验场得到的多轴载荷对于零件上的任意一点的应力循环是等效的，也就是希望等效载荷和实测的多轴载荷都有相等的强度。本发明中，可利用伪损伤来描述这种强度，不同的是多轴载荷等效方法的这种等效性要求对任意的一组线性组合均成立。
假设多轴载荷Lk(t)的任意线性组合Lc(t)，和等效载荷的任意线性组合分别如下所示：
Lc(t)＝c1·L1(t)+c2·L2(t)+…+cn·Ln(t)＝c·L(t)；
其中，上两式中的ck为任意实数，ck的任意性使得Lc(t)和包含了零部件上任意一点分别在Lk(t)和/>载荷作用下的某种应力值。这些应力值包含任意一点上的应力循环的各个应力循环分量σij(x,t)和能够通过各应力循环分量的线性组合表示的应力，例如临界平面法中各个临界平面的法向和切向应力。
在本发明中，设fd(L(t))代表对多轴载荷L(t)对应的伪损伤数，基于求取等效载荷的需求，对于任意实数的ck，Lc(t)和应当产生相同的伪损伤数，也即存在：
值得强调的是，对于任意的一组实数ck，实数ck代表了n维空间内的一个向量c＝{c1,c2,…,cn}T，在进行等效载荷求取时，可不需要研究空间内所有的向量，见如下说明：
对于两组线性相关的向量c1和c2，由于是线性相关关系，设c2＝a·c1，如果Lc1(t)＝c1·L(t)和两组信号有相同的伪损伤数，那么有如下关系：
因此，Lc2(t)＝c2·L(t)和两组信号的损伤数也相等：
假设对于某个应力循环σ(t)的雨流计数结果为σi(t)，则a·σ(t)的雨流计数结果为a·σi(t)，因此还存在如下关系：
所以当则有如下关系：
可以看出，对于一系列的线性相关的向量中，只需要考虑其中的一个向量即可，因此，可归一化这组向量使其满足这样就可以把一个多轴载荷的无限种可能的空间离散为若干个有限的组合，这也是在步骤S20中可以将多轴载荷拆解成多个不同方向的分解载荷的依据。
在一实施方式中，将所述多轴载荷投影到载荷空间中，以得到若干个方向上的分解载荷之前，还需对实测的多轴载荷进行如下预处理：去除毛刺、均值修正和雨流滤波处理，将经过所述预处理的所述多轴载荷投影到所述载荷空间中，以得到所述若干个方向上的分解载荷。
可以理解，由于实测的测试条件等限制，采集的载荷谱难免存在干扰信息，故本发明在进行投影分解处理之前，需对实测得到的载荷谱进行滤波、去毛刺等预处理，以从经过预处理后的载荷谱中得到符合规范的多轴载荷，有利于后续的等效载荷处理计算，降低后续计算的难度和提高计算准确性。
从上述描述可见，在本发明实施例中，可将多轴载荷投影进行分解，一得到若干个方向上的分解载荷，从而可以将多轴字载荷的无限种可能的空间离散为若干个有限的组合，也可以极大地减少了后续等效处理的计算量，也提高了可实施性。
S30：分别对所述若干个方向上的分解载荷进行伪损伤计算，以获取所述分解载荷的伪损伤。
在将多轴载荷投影到载荷空间中，以得到若干个方向上的分解载荷之后，便可分别对每个方向上的分解载荷进行伪损伤计算，以获取不同方向上的分解载荷的伪损伤。
其中，步骤S30中，也即分别对所述若干个方向上的分解载荷进行伪损伤计算，以获取所述分解载荷的伪损伤，具体包括如下步骤：
S31：利用雨流计数法对所述分解载荷进行雨流计数处理，以获取所述分解载荷对应的多个应力循环。
S32：根据所述目标零部件材料对应的S-N曲线计算各个所述应力循环对应的伪损伤值。
S33：利用损伤累积理论对各个所述应力循环对应的损伤值进行累积处理，以获取所述分解载荷的伪损伤。
对于步骤S31-S32，由于在描述本发明实施例之前，已经对伪损伤等相关术语进行描述，在此不重复描述，其中，分解载荷的伪损伤的计算方式可如下所示：
设分解载荷用σG表示，首先利用雨流统计法将该将分解载荷σG拆解为若干个应力循环每一个应力循环/>都根据目标零部件材料对应的S-N曲线计算每个应力循环对应的损伤值/>
其中，示例性的，单个应力循环产生的伪损伤值/>可通过如下方式计算：
表示所述应力循环，/>表示伪损伤值，N表示目标零部件的疲劳寿命，b为对应目标零部件对应材料的S-N曲线的斜率，S0为目标零部件对应材料的S-N曲线的斜截距。
可以理解，利用上述计算方式便可知道分解载荷的每一个应力循环对应的伪损伤值/>
在得到分解载荷σG的每一个应力循环对应的伪损伤值/>之后，针对某个分解载荷，利用损伤累积理论，例如，利用米勒线性损伤累积理论对该分解载荷σG拆解出的若干个应力循环/>对应的伪损伤值/>进行叠加，便可得该分解载荷σG对应的伪损伤为DG为：
S40：根据所述分解载荷的伪损伤和多轴载荷循环次数确定所述分解载荷的总伪损伤。
可以理解，在经过步骤得到每个方向上的分解载荷的伪损伤之后，便可根据所述分解载荷的伪损伤和多轴载荷循环次数确定所述分解载荷的总伪损伤。其中，此过程可依据多轴载荷的试验路面数量和试验路面对应的多轴载荷循环次数，确定出述分解载荷在不同试验路面的伪损伤进行叠加，便可回到分解载荷的总伪损伤。
需要说明的是，在前述中已经知道伪损伤之比即为伪损伤数之比，故在根据所述分解载荷的伪损伤和多轴载荷循环次数确定所述分解载荷的总伪损伤时，可对应确定出分解载荷的总伪损伤数便可。此时，对应的，每个分解载荷的总伪损伤数可表示为：
此时，dTarget表示每个分解载荷的总伪损伤数，nk为第k路面对应的循环次数dk为第k路面时对应的伪损伤数，伪损伤数的计算可参阅前述描述。
S50：初始化等效载荷循环次数，并构建与所述若干个方向的分解载荷对应的预设信号形式的等效载荷式。
在发明中，为了得到与实测得到的多轴载荷等效的等效载荷，需预先构建某种预设信号形式的等效载荷式，并初始该等效载荷式的循环次数，本发明中为了区别多轴载荷循环次数，将初始该等效载荷式的循环次数称为等效载荷循环次数。需要说明的是，该预设信号形式可以是指正弦波形式，或者其他信号形式的等效载荷，本发明不做限制。
值得说明的是，为了缩短后续载荷加载试验时间和简化载荷迭代求解的过程，本发明实施例中采用正弦波载荷作为等效载荷代替实测的多轴载荷，因此，本发明可构建如下形式的等效载荷式：
其中，分别表示多轴载荷的若干方向上的分解载荷对应的等效载荷式。
另外，等效载荷循环次数一般可事先取定的，一般底盘零件问题取10万次到30万次左右，值得注意的是，在后续构建载荷等效关系式时，还会涉及等效载荷的初始值选，由于初始值对最后迭代解决结果影响不大，这里不展开详细说明，本发明中一般取10kN左右即可。
S60：根据所述等效载荷循环次数和所述等效载荷式确定所述等效载荷对应的总伪损伤式。
等效载荷的总伪损伤与等效载荷循环次数相关，此过程可对应参阅前述计算分解载荷的总伪损伤的过程，在此不做描述。
S70：令所述若干个方向的分解载荷的总伪损伤分别等于所述等效载荷对应的总伪损伤式，以获取载荷等效关系式。
依据前述提及的若载荷等效，则相应的伪损伤和伪损伤数为1的原理，本发明可令所述若干个方向的分解载荷的总伪损伤等于对应的等效载荷的总伪损伤式，以获取载荷等效关系式。
值得说明的是，经过前述说明，若伪损伤为1，则伪损伤数也应为1，故，在步骤S60-S70中，可根据所述等效载荷循环次数和所述等效载荷式确定所述等效载荷对应的总伪损伤数式，并令所述若干个方向的分解载荷的总伪损伤数分别等于所述等效载荷对应的总伪损伤数式，也可获取相应的载荷等效关系式，本发明将以采用伪损伤数构建载荷等效关系式为例进行说明。
需要说明的是，等效载荷的构造方法是多样的，能够使等效载荷和多轴载荷在各方向上的分解载荷的总伪损伤都尽量接近是构造的等效载荷的目标，因此，令所述若干个方向的分解载荷的总伪损伤数等于对应的等效载荷的总伪损伤数式之后，可得如下载荷等效关系式：
其中，表示构建的等效载荷式对应的总伪损伤数式，fd(c1·L(t))表示分解载荷L(t)对应的总伪损伤数。
S80：对所述载荷等效关系式进行迭代求解，以获取所述等效载荷式的参数。
S90：根据所述参数得到所述等效载荷式对应的等效载荷。
在得到载荷等效关系式之后，便可对所述载荷等效关系式进行迭代求解，以获取所述等效载荷式的参数，也即以等幅值的正弦波等效载荷为例，最后可求得频率、幅值和相位参数以及循环次数等参数。需要说明的是，在一实施方式中，可通过最小二乘和最优解方法求解，其中，迭代求解目标为各所述方向上等效载荷的总伪损伤与对应的分解载荷的总伪损伤之差的平方和取最小值为所述迭代求解的目标；或者各方向上分解载荷的总伪损伤数与对应等效载荷的总伪损伤数之差的平方和取最小值为目标或者，也即迭代求解目标如下式所示：
因此，利用前述载荷等效关系式进行迭代求解便可求解得到相应的等效载荷和载荷循环次数，具体地迭代求解过程在此不展开描述。
本发明提供了一种多轴载荷等效处理方法，通过投影法对多轴载荷进行拆解成多个分解载荷，将多轴问题转化为单个载荷问题，并分别对若干个方向上的分解载荷进行伪损伤计算，并获取所述分解载荷的伪损伤；最后依据各个分解载荷的伪损伤构建相应载荷等效关系式，并对载荷等效关系式进行迭代求解，以获取所述等效载荷式的参数；最后根据所述参数得到所述等效载荷式对应的等效载荷，从而把复杂的多轴随机的多轴载荷等效地转化为等效载荷进行台架试验，由于保留了分解载荷是由多轴载荷分解得到，通过各个分解载荷计算出的等效载荷之间保留多轴载荷的相位关系，从而提高了等效载荷的等效性，也就提高了后续零部件耐久试验的准确性。
值得说明的是，本方法可建出多轴载荷对应的等效载荷，理论上无论多少轴的多轴载荷均可都适合，适应性较强，并且通过雨流投影计数的方式，将单轴载荷等效方法拓展到多轴载荷问题上，具体地将单轴载荷拆解分多个分解载荷进行后续等效，从而可以将载荷空间内的无限个方向的问题离散为有限个方向的问题，便于构建载荷等效关系式和迭代求解；另外，在迭代求解过程中，可运用最小二乘的思想求解近似解，而不是精确解，提高了方案的可实施性；另外，将多轴随机的载荷转化为正弦波载荷，简化了后续试验等效载荷的加载难度、大幅提升了试验效率的同时还保留了多轴载荷的相位关系，提升了后续试验的准确性。
应理解，上述实施例中各步骤的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本发明实施例的实施过程构成任何限定。
在一实施例中，提供一种多轴载荷等效处理装置，该多轴载荷等效处理装置与上述实施例中多轴载荷等效处理方法中一一对应。如图4所示，该胎压展示装置包括采集模块101、投影模块102、计算模块103、第一确定模块104、构建模块105、第二确定模块106、第一获取模块107、第二获取模块108、和第三获取模块109。各功能模块详细说明如下：
采集模块101，用于采集目标零部件在试验路面实际所受的多轴载荷，并确定多轴载荷循环次数；
投影模块102，用于将所述多轴载荷投影到载荷空间中，以得到若干个方向上的分解载荷；
计算模块103，用于分别对所述若干个方向上的分解载荷进行伪损伤计算，以获取所述分解载荷的伪损伤；
第一确定模块104，用于根据所述分解载荷的伪损伤和多轴载荷循环次数确定所述分解载荷的总伪损伤；
构建模块105，用于初始化等效载荷循环次数，并构建与所述若干个方向的分解载荷对应的预设信号形式的等效载荷式；
第二确定模块106，用于根据所述等效载荷循环次数和所述等效载荷式确定所述等效载荷对应的总伪损伤式；
第一获取模块107，用于令所述若干个方向的分解载荷的总伪损伤分别等于所述等效载荷对应的总伪损伤式，以获取载荷等效关系式；
第二获取模块108，用于对所述载荷等效关系式进行迭代求解，以获取所述等效载荷式的参数；
第三获取模块109，用于根据所述参数得到所述等效载荷式对应的等效载荷。
在一实施例中，投影模块具体用于：
确定所述多轴载荷中不共线载荷的数量；
根据所述多轴载荷中不共线载荷的数量确定投影方向数量和投影方向；
按照所述投影方向数量和投影方向将所述多轴载荷投影到载荷空间中，以得到投影方向数量对应的若干个所述投影方向上的分解载荷。
在一实施例中，计算模块，具体用于：
利用雨流计数法对所述分解载荷进行雨流计数处理，以获取所述分解载荷对应的多个应力循环；
根据所述目标零部件材料对应的S-N曲线计算各个所述应力循环对应的伪损伤值；
利用损伤累积理论对各个所述应力循环对应的损伤值进行累积处理，以获取所述分解载荷的伪损伤。
在一实施例中，所述各个应力循环对应的伪损伤值通过下式计算：
表示所述应力循环，/>表示伪损伤值，N表示目标零部件的疲劳寿命，b为对应目标零部件对应材料的S-N曲线的斜率，S0为目标零部件对应材料的S-N曲线的斜截距；
所述分解载荷的伪损伤通过下式计算：
表示所述分解载荷的伪损伤。
在一实施例中，构建模块具体用于：
构建与所述若干个方向的分解载荷对应的正弦信号形式的等效载荷式。
在一实施例中，第二获取模块具体用于：
通过最小二乘法对所述载荷等效关系式进行迭代求解，以获取所述等效载荷式的参数，其中，各所述方向上等效载荷的总伪损伤与对应的分解载荷的总伪损伤之差的平方和取最小值为所述迭代求解的目标。
在一实施例中，所述装置还包括预处理模块；
所述预处理模块，用于对所述多轴载荷进行预处理，所述预处理包括去除毛刺、均值修正和滤波处理；
所述投影模块，用于将经过所述预处理的所述多轴载荷投影到所述载荷空间中，以得到所述若干个方向上的分解载荷。
本发明提供了一种多轴载荷等效处理装置，通过投影法对多轴载荷进行拆解成多个分解载荷，将多轴问题转化为单个载荷问题，并分别对若干个方向上的分解载荷进行伪损伤计算，并获取所述分解载荷的伪损伤；最后依据各个分解载荷的伪损伤构建相应载荷等效关系式，并对载荷等效关系式进行迭代求解，以获取所述等效载荷式的参数；最后根据所述参数得到所述等效载荷式对应的等效载荷，从而把复杂的多轴随机的多轴载荷等效地转化为等效载荷进行台架试验，由于保留了分解载荷是由多轴载荷分解得到，通过各个分解载荷计算出的等效载荷之间保留多轴载荷的相位关系，从而提高了等效载荷的等效性，也就提高了后续零部件耐久试验的准确性。
关于多轴载荷等效处理装置的具体限定可以参见上文中对于多轴载荷等效处理方法的限定，在此不再赘述。上述多轴载荷等效处理装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中，提供了一种计算机设备，其内部结构图可以如图5所示。该计算机设备包括通过***总线连接的处理器、存储器和网络接口。该控制器可通过网络接口获取所需信息，其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括存储介质。该存储器存储有计算机程序。该存储器为计算机程序的运行提供环境。该计算机程序被处理器执行时以实现一种多轴载荷等效处理方法。
在一个实施例中，提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现以下步骤：
采集目标零部件在试验路面实际所受的多轴载荷，并确定多轴载荷循环次数；
将所述多轴载荷投影到载荷空间中，以得到若干个方向上的分解载荷；
分别对所述若干个方向上的分解载荷进行伪损伤计算，以获取所述分解载荷的伪损伤；
根据所述分解载荷的伪损伤和多轴载荷循环次数确定所述分解载荷的总伪损伤；
初始化等效载荷循环次数，并构建与所述若干个方向的分解载荷对应的预设信号形式的等效载荷式；
根据所述等效载荷循环次数和所述等效载荷式确定所述等效载荷对应的总伪损伤式；
令所述若干个方向的分解载荷的总伪损伤分别等于所述等效载荷对应的总伪损伤式，以获取载荷等效关系式；
对所述载荷等效关系式进行迭代求解，以获取所述等效载荷式的参数；
根据所述参数得到所述等效载荷式对应的等效载荷。
在一个实施例中，提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如下步骤：
采集目标零部件在试验路面实际所受的多轴载荷，并确定多轴载荷循环次数；
将所述多轴载荷投影到载荷空间中，以得到若干个方向上的分解载荷；
分别对所述若干个方向上的分解载荷进行伪损伤计算，以获取所述分解载荷的伪损伤；
根据所述分解载荷的伪损伤和多轴载荷循环次数确定所述分解载荷的总伪损伤；
初始化等效载荷循环次数，并构建与所述若干个方向的分解载荷对应的预设信号形式的等效载荷式；
根据所述等效载荷循环次数和所述等效载荷式确定所述等效载荷对应的总伪损伤式；
令所述若干个方向的分解载荷的总伪损伤分别等于所述等效载荷对应的总伪损伤式，以获取载荷等效关系式；
对所述载荷等效关系式进行迭代求解，以获取所述等效载荷式的参数；
根据所述参数得到所述等效载荷式对应的等效载荷。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，仅以上述各功能单元、模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能单元、模块完成，即将所述装置的内部结构划分成不同的功能单元或模块，以完成以上描述的全部或者部分功能。
以上所述实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围，均应包含在本发明的保护范围之内。
Claims (10)
1.一种多轴载荷等效处理方法，其特征在于，包括：
采集目标零部件在试验路面实际所受的多轴载荷，并确定多轴载荷循环次数；
将所述多轴载荷投影到载荷空间中，以得到若干个方向上的分解载荷；
分别对所述若干个方向上的分解载荷进行伪损伤计算，以获取所述分解载荷的伪损伤；
根据所述分解载荷的伪损伤和多轴载荷循环次数确定所述分解载荷的总伪损伤；
初始化等效载荷循环次数，并构建与所述若干个方向的分解载荷对应的预设信号形式的等效载荷式；
根据所述等效载荷循环次数和所述等效载荷式确定所述等效载荷对应的总伪损伤式；
令所述若干个方向的分解载荷的总伪损伤分别等于所述等效载荷对应的总伪损伤式，以获取载荷等效关系式；
对所述载荷等效关系式进行迭代求解，以获取所述等效载荷式的参数；
根据所述参数得到所述等效载荷式对应的等效载荷。
2.如权利要求1所述的多轴载荷等效处理方法，其特征在于，所述将所述多轴载荷投影到载荷空间中，以得到若干个方向上的分解载荷，包括：
确定所述多轴载荷中不共线载荷的数量；
根据所述多轴载荷中不共线载荷的数量确定投影方向数量和投影方向；
按照所述投影方向数量和投影方向将所述多轴载荷投影到载荷空间中，以得到投影方向数量对应的若干个所述投影方向上的分解载荷。
3.如权利要求1所述的多轴载荷等效处理方法，其特征在于，所述分别对所述若干个方向上的分解载荷进行伪损伤计算，以获取所述分解载荷的伪损伤，包括：
利用雨流计数法对所述分解载荷进行雨流计数处理，以获取所述分解载荷对应的多个应力循环；
根据所述目标零部件材料对应的S-N曲线计算各个所述应力循环对应的伪损伤值；
利用损伤累积理论对各个所述应力循环对应的伪损伤值进行累积处理，以获取所述分解载荷的伪损伤。
4.如权利要求3所述的多轴载荷等效处理方法，其特征在于，所述各个应力循环对应的伪损伤值通过下式计算：
表示所述应力循环，/>表示伪损伤值，N表示目标零部件的疲劳寿命，b为对应目标零部件对应材料的S-N曲线的斜率，S0为目标零部件对应材料的S-N曲线的斜截距；
所述分解载荷的伪损伤通过下式计算：
表示所述分解载荷的伪损伤。
5.如权利要求1-4任一项所述的多轴载荷等效处理方法，其特征在于，所述构建与所述若干个方向的分解载荷对应的预设信号形式的等效载荷式，包括：
构建与所述若干个方向的分解载荷对应的正弦信号形式的等效载荷式。
6.如权利要求5所述的多轴载荷等效处理方法，其特征在于，所述对所述载荷等效关系式进行迭代求解，以获取所述等效载荷式的参数，包括：
通过最小二乘法对所述载荷等效关系式进行迭代求解，以获取所述等效载荷式的参数，其中，各所述方向上等效载荷的总伪损伤与对应的分解载荷的总伪损伤之差的平方和取最小值为所述迭代求解的目标。
7.如权利要求1-4任一项所述的多轴载荷等效处理方法，其特征在于，所述将所述多轴载荷投影到载荷空间中，以得到若干个方向上的分解载荷之前，所述方法还包括：
对所述多轴载荷进行预处理，所述预处理包括去除毛刺、均值修正和滤波处理；
所述将所述多轴载荷投影到载荷空间中，以得到若干个方向上的分解载荷，包括：
将经过所述预处理的所述多轴载荷投影到所述载荷空间中，以得到所述若干个方向上的分解载荷。
8.一种多轴载荷等效处理装置，其特征在于，包括：
采集模块，用于采集目标零部件在试验路面实际所受的多轴载荷，并确定多轴载荷循环次数；
投影模块，用于将所述多轴载荷投影到载荷空间中，以得到若干个方向上的分解载荷；
计算模块，用于分别对所述若干个方向上的分解载荷进行伪损伤计算，以获取所述分解载荷的伪损伤；
第一确定模块，用于根据所述分解载荷的伪损伤和多轴载荷循环次数确定所述分解载荷的总伪损伤；
构建模块，用于初始化等效载荷循环次数，并构建与所述若干个方向的分解载荷对应的预设信号形式的等效载荷式；
第二确定模块，用于根据所述等效载荷循环次数和所述等效载荷式确定所述等效载荷对应的总伪损伤式；
第一获取模块，用于令所述若干个方向的分解载荷的总伪损伤分别等于所述等效载荷对应的总伪损伤式，以获取载荷等效关系式；
第二获取模块，用于对所述载荷等效关系式进行迭代求解，以获取所述等效载荷式的参数；
第三获取模块，用于根据所述参数得到所述等效载荷式对应的等效载荷。
9.一种计算机设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1至7任一项多轴载荷等效处理方法的步骤。
10.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至7任一项多轴载荷等效处理方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010953720.5A CN112100563B (zh) | 2020-09-11 | 2020-09-11 | 多轴载荷等效处理方法、装置、计算机设备和介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010953720.5A CN112100563B (zh) | 2020-09-11 | 2020-09-11 | 多轴载荷等效处理方法、装置、计算机设备和介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112100563A CN112100563A (zh) | 2020-12-18 |
| CN112100563B true CN112100563B (zh) | 2024-04-19 |
Family
ID=73751425
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010953720.5A Active CN112100563B (zh) | 2020-09-11 | 2020-09-11 | 多轴载荷等效处理方法、装置、计算机设备和介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112100563B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114371018B (zh) * | 2021-11-30 | 2024-02-27 | 江铃汽车股份有限公司 | 一种汽车动力总成悬置支架的台架耐久试验方法及*** |
Citations (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7677094B1 (en) * | 2008-09-30 | 2010-03-16 | Nissan Technical Center North America, Inc. | Method and system for testing road load durability of truck rear bed |
| JP2013044667A (ja) * | 2011-08-25 | 2013-03-04 | Ihi Corp | 多軸疲労寿命評価方法 |
| CN105718633A (zh) * | 2016-01-15 | 2016-06-29 | 重庆长安汽车股份有限公司 | 一种底盘件的载荷分析方法 |
| CN105784339A (zh) * | 2016-02-24 | 2016-07-20 | 安徽工业大学 | 非线性结构零件损伤循环计数方法及其疲劳寿命分析方法 |
| KR101706508B1 (ko) * | 2016-07-12 | 2017-02-16 | 한국에너지기술연구원 | 풍력 발전기의 피로 해석 및 등가하중 해석 시스템 |
| CN106739918A (zh) * | 2016-12-13 | 2017-05-31 | 广州汽车集团股份有限公司 | 车身姿态主动控制***及主动控制方法 |
| CN108733879A (zh) * | 2018-04-12 | 2018-11-02 | 南京航空航天大学 | 一种多轴非比例载荷下剪应力幅值计算方法 |
| CN109241577A (zh) * | 2018-08-14 | 2019-01-18 | 南京航空航天大学 | 一种预测多轴疲劳寿命的修正fs法 |
| DE102017217561A1 (de) * | 2017-10-04 | 2019-04-04 | Bayerische Motoren Werke Aktiengesellschaft | Verfahren zur Betriebsfestigkeitsanalyse eines Bauteils |
| CN110082120A (zh) * | 2019-04-24 | 2019-08-02 | 一汽-大众汽车有限公司 | 加速耐久性试验的路线规划方法和装置 |
| JP2019203765A (ja) * | 2018-05-23 | 2019-11-28 | 株式会社豊田中央研究所 | 大変形域又は多軸応力下の応力ひずみ特性計測装置、及び大変形域又は多軸応力下の応力ひずみ特性解析装置 |
| CN110646216A (zh) * | 2019-09-11 | 2020-01-03 | 湖南猎豹汽车股份有限公司 | 一种由路谱载荷生成台架试验载荷的方法 |
| CN110955950A (zh) * | 2018-09-25 | 2020-04-03 | 湖南工业大学 | 一种基于耗散能预测任意应力比下焊缝疲劳寿命的方法 |
| CN111159871A (zh) * | 2019-12-23 | 2020-05-15 | 北京工业大学 | 一种基于路径曲线积分的随机多轴循环计数方法 |
| CN111222199A (zh) * | 2019-11-13 | 2020-06-02 | 中国汽车工程研究院股份有限公司 | 用户与试验场关联时关键指标选择与等效计算方法 |
| CN111289246A (zh) * | 2018-12-06 | 2020-06-16 | 广州汽车集团股份有限公司 | 一种双离合自动变速器驻车机构耐久试验方法 |
| CN111581893A (zh) * | 2020-04-03 | 2020-08-25 | 上海理工大学 | 一种电驱动总成机械***可靠性试验载荷谱的编制方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8649985B2 (en) * | 2009-01-08 | 2014-02-11 | Battelle Memorial Institute | Path-dependent cycle counting and multi-axial fatigue evaluation of engineering structures |
| NO3074325T3 (zh) * | 2013-12-23 | 2018-02-24 | ||
| US11471982B2 (en) * | 2017-08-18 | 2022-10-18 | The Regents Of The University Of Michigan | Unified fatigue life evaluation method for welded structures |
-
2020
- 2020-09-11 CN CN202010953720.5A patent/CN112100563B/zh active Active
Patent Citations (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7677094B1 (en) * | 2008-09-30 | 2010-03-16 | Nissan Technical Center North America, Inc. | Method and system for testing road load durability of truck rear bed |
| JP2013044667A (ja) * | 2011-08-25 | 2013-03-04 | Ihi Corp | 多軸疲労寿命評価方法 |
| CN105718633A (zh) * | 2016-01-15 | 2016-06-29 | 重庆长安汽车股份有限公司 | 一种底盘件的载荷分析方法 |
| CN105784339A (zh) * | 2016-02-24 | 2016-07-20 | 安徽工业大学 | 非线性结构零件损伤循环计数方法及其疲劳寿命分析方法 |
| KR101706508B1 (ko) * | 2016-07-12 | 2017-02-16 | 한국에너지기술연구원 | 풍력 발전기의 피로 해석 및 등가하중 해석 시스템 |
| CN106739918A (zh) * | 2016-12-13 | 2017-05-31 | 广州汽车集团股份有限公司 | 车身姿态主动控制***及主动控制方法 |
| DE102017217561A1 (de) * | 2017-10-04 | 2019-04-04 | Bayerische Motoren Werke Aktiengesellschaft | Verfahren zur Betriebsfestigkeitsanalyse eines Bauteils |
| CN108733879A (zh) * | 2018-04-12 | 2018-11-02 | 南京航空航天大学 | 一种多轴非比例载荷下剪应力幅值计算方法 |
| JP2019203765A (ja) * | 2018-05-23 | 2019-11-28 | 株式会社豊田中央研究所 | 大変形域又は多軸応力下の応力ひずみ特性計測装置、及び大変形域又は多軸応力下の応力ひずみ特性解析装置 |
| CN109241577A (zh) * | 2018-08-14 | 2019-01-18 | 南京航空航天大学 | 一种预测多轴疲劳寿命的修正fs法 |
| CN110955950A (zh) * | 2018-09-25 | 2020-04-03 | 湖南工业大学 | 一种基于耗散能预测任意应力比下焊缝疲劳寿命的方法 |
| CN111289246A (zh) * | 2018-12-06 | 2020-06-16 | 广州汽车集团股份有限公司 | 一种双离合自动变速器驻车机构耐久试验方法 |
| CN110082120A (zh) * | 2019-04-24 | 2019-08-02 | 一汽-大众汽车有限公司 | 加速耐久性试验的路线规划方法和装置 |
| CN110646216A (zh) * | 2019-09-11 | 2020-01-03 | 湖南猎豹汽车股份有限公司 | 一种由路谱载荷生成台架试验载荷的方法 |
| CN111222199A (zh) * | 2019-11-13 | 2020-06-02 | 中国汽车工程研究院股份有限公司 | 用户与试验场关联时关键指标选择与等效计算方法 |
| CN111159871A (zh) * | 2019-12-23 | 2020-05-15 | 北京工业大学 | 一种基于路径曲线积分的随机多轴循环计数方法 |
| CN111581893A (zh) * | 2020-04-03 | 2020-08-25 | 上海理工大学 | 一种电驱动总成机械***可靠性试验载荷谱的编制方法 |
Non-Patent Citations (10)
| Title |
|---|
| A path-dependent cycle counting method for variable-amplitude multi-axial loading;DONG P;《International Journal of Fatigue》;第32卷(第4期);720-734 * |
| Life prediction techniques for variable amplitude multiaxial fatigue—part 1：theories;WANG C H;《Journal of Engineering Materials and Technology》;第118卷(第3期);367-370 * |
| 基于CAE仿真和矩阵标定的汽车零部件载荷测试方法;程凯;计算机辅助工程;第25卷(第02期);43-47 * |
| 基于损伤等效的多轴疲劳试验谱编制研究;赵勇铭;航空动力学报;第24卷(第09期);2026-2032 * |
| 基于时域损伤编辑法的整车道路模拟试验谱编制;花菲菲;郑松林;冯金芝;于佳伟;;机械强度;第37卷(第05期);964-971 * |
| 基于等损伤的试验场关联及其在Excel中的实现;钟志宏;工业技术创新;第04卷(第05期);95-98 * |
| 基于等效应力的试验载荷谱开发与研究;黄庆;汽车技术(第02期);53-57 * |
| 基于雨流法的汽车底盘件疲劳耐久试验方法研究;武文超;汽车技术(第12期);46-49、60 * |
| 多轴随机载荷下疲劳寿命预测方法的研究;宋宇;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;第2015年卷(第05期);C029-28 * |
| 目标用户道路谱与试验场道路谱的载荷当量等效模拟研究;吴建国;汽车技术(第07期);21-24 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112100563A (zh) | 2020-12-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Cumbo et al. | Kalman-based load identification and full-field estimation analysis on industrial test case | |
| Döhler et al. | Subspace-based damage detection under changes in the ambient excitation statistics | |
| WO2019035098A1 (en) | METHOD FOR EVALUATING UNIFORM FATIGUE RESISTANCE FOR WELDED STRUCTURES | |
| Braccesi et al. | The frequency domain approach in virtual fatigue estimation of non-linear systems: The problem of non-Gaussian states of stress | |
| CN112100563B (zh) | 多轴载荷等效处理方法、装置、计算机设备和介质 | |
| CN112781616B (zh) | 星敏感器在轨测量低频误差分析方法、装置和存储介质 | |
| Xu et al. | Continuous test suite augmentation in software product lines | |
| JP2020109651A (ja) | 材料モデルの安定性に基づく制約付き数値較正のためのシステム及び方法 | |
| CN106997410B (zh) | 一种基于模态应变能的损伤发生的判断方法 | |
| Cannizzaro et al. | Explicit dynamic response of damaged beams with application to uncertain and identification problems | |
| CN105512372B (zh) | 模型化的星载数据处理仿真测试方法 | |
| CN114065373A (zh) | 一种汽车控制臂台架耐久试验方法、装置以及设备 | |
| Forbes et al. | Generation of numerical artefacts for geometric form and tolerance assessment | |
| Friswell | Damage identification using inverse methods | |
| CN115204020B (zh) | 电驱桥***强度分析方法、***、存储介质及测试设备 | |
| Bel Knani et al. | Fatigue damage assessment of a car body-in-white using a frequency-domain approach | |
| Gordis et al. | Minimum condition number by orthogonal projection row selection of artificial boundary conditions for finite element model update and damage detection | |
| CN113704085A (zh) | 用于检查技术***的方法和设备 | |
| Yan | A bayesian approach for identification of structural crack using strain measurements | |
| CN114741887B (zh) | 保险杠耐久性分析方法、装置、计算机设备、存储介质 | |
| Luo et al. | A rapid method to predict biaxial fatigue life of automotive wheels using proper orthogonal decomposition and radial basis function algorithm | |
| Liu et al. | Stochastic analysis for in-plane dynamic responses of low-speed uniformity of tires due to geometric defects | |
| Ding et al. | A Simplified Small Scale Yielding Estimate Method for Determination of Second Elastic-Plastic Fracture Mechanics Parameter | |
| Mihalache et al. | Crack reconstruction in ferromagnetic materials using nonlinear FEM–BEM scheme and neural networks | |
| US12032880B2 (en) | Method and device for simulating a technical system |
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 |