CN104217094A - 用于计算结构的疲劳以及疲劳破坏的方法 - Google Patents
用于计算结构的疲劳以及疲劳破坏的方法 Download PDFInfo
- Publication number
- CN104217094A CN104217094A CN201410201290.6A CN201410201290A CN104217094A CN 104217094 A CN104217094 A CN 104217094A CN 201410201290 A CN201410201290 A CN 201410201290A CN 104217094 A CN104217094 A CN 104217094A
- Authority
- CN
- China
- Prior art keywords
- change
- stress
- damage
- strain
- load
- 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 64
- 230000008569 process Effects 0.000 title abstract description 7
- 230000006378 damage Effects 0.000 claims abstract description 88
- 230000008859 change Effects 0.000 claims abstract description 74
- 238000004364 calculation method Methods 0.000 claims abstract description 39
- 239000000463 material Substances 0.000 claims description 47
- 230000003111 delayed effect Effects 0.000 claims description 18
- 150000001875 compounds Chemical class 0.000 claims description 9
- 239000000126 substance Substances 0.000 claims description 9
- 238000004458 analytical method Methods 0.000 claims description 7
- 238000004519 manufacturing process Methods 0.000 claims description 5
- 230000036962 time dependent Effects 0.000 claims description 5
- 239000013536 elastomeric material Substances 0.000 claims description 3
- 238000012986 modification Methods 0.000 claims description 3
- 230000004048 modification Effects 0.000 claims description 3
- 238000004590 computer program Methods 0.000 claims 2
- 238000005457 optimization Methods 0.000 abstract description 3
- 230000035882 stress Effects 0.000 description 66
- 239000002131 composite material Substances 0.000 description 35
- 230000004087 circulation Effects 0.000 description 13
- 238000009825 accumulation Methods 0.000 description 11
- 230000006399 behavior Effects 0.000 description 9
- 238000010276 construction Methods 0.000 description 7
- 238000013461 design Methods 0.000 description 7
- 239000000047 product Substances 0.000 description 7
- 238000004088 simulation Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 238000006073 displacement reaction Methods 0.000 description 5
- 238000012360 testing method Methods 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 238000005452 bending Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000005484 gravity Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000012938 design process Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000003446 memory effect Effects 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 239000004033 plastic Substances 0.000 description 2
- 229920003023 plastic Polymers 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000005050 thermomechanical fatigue Methods 0.000 description 2
- 229920002430 Fibre-reinforced plastic Polymers 0.000 description 1
- 241000276425 Xiphophorus maculatus Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 239000000109 continuous material Substances 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006355 external stress Effects 0.000 description 1
- 239000011151 fibre-reinforced plastic Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000021715 photosynthesis, light harvesting Effects 0.000 description 1
- 238000013031 physical testing Methods 0.000 description 1
- 238000002910 structure generation Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M5/00—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L5/00—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本申请提供一种用于计算结构的疲劳以及疲劳破坏的方法。本发明涉及用于虚拟地预测结构的耐久性性能从而能够优化该耐久性性能的方法。在第一步骤，通过一系列计算点对该结构进行建模。然后针对每个点确定由载荷循环所引起的应力及应变。然后，预测并存储由于载荷循环而引起的累积损伤。为了预测，首先使用滞后算子作为载荷随时间变化的函数计算应力沿着滞后分支的一部分的变化，其次使用应力的变化以及所存储的累积损伤，计算损伤的变化。然后基于这些属性的变化来计算应力及应变的进一步的变化，以确定新的适合的滞后分支。然后，作为载荷随时间进一步变化的函数计算应力沿着适合的滞后分支的另外的部分的进一步的变化。最后相应地制造该结构。
Description
技术领域
本发明涉及预测结构的耐久性性能以及疲劳特性的领域。
背景技术
由于轻质材料允许使用较少的生态足迹来提供具有优越机械属性的产品，交通运输、风能及机械行业的制造商越来越多地使用轻质材料。尤其是在交通运输(汽车、航空……)中，大量采用轻质设计材料(例如复合材料)将会是满足日益严格的二氧化碳排放标准的唯一可行的路径。然而目前，在工业设计及开发过程中大规模使用轻质材料受限于缺乏预测性的建模工具来预测轻质材料结构的宏观特性。如果可以基于虚拟模型做出较好的机械预测，则可以基于虚拟模拟(而不是基于最新的且昂贵的物理测试)来优化产品。
复合材料通常可以承受多个载荷循环(即，具有好的疲劳特性)。但是，甚至改变部件的刚度的第一破坏性事件发生地相当早并且发生在相当小的载荷循环下。对于产品设计工程师来说，这意味着他们不能设计出具有“完全不会发生疲劳破坏”的要求的产品，因为这将导致典型的结构的超裕度设计。反之，工程师希望能够预测整个生命周期的疲劳特性-包括刚度的变化，以及不同的破坏机制的交叉影响。通过在产品设计确定方法中使用对疲劳特性的更好的理解，将可以实现更好、更轻且更经济的设计。
例如，在汽车应用中，对疲劳载荷的分解和合成在过去的几十年中已大大地改进。针对当前由金属制造的汽车，正在使用成熟的虚拟设计优化方法，这使得汽车制造商能够按需定制他们的汽车以足以满足实际使用，并且能够量身打造他们的汽车以满足不同的市场的要求，同时确保避免超裕度设计。这涉及针对不同的道路施加完全可变的载荷的复杂的载荷安排。反之，对于复合材料来说，还没有成熟的工具来实际预测产品性能。这意味着在设计过程的较后阶段需要昂贵的物理测试方法来检查/确认复合材料设计是否满足要求。
要解决的问题是适当引进新的虚拟方法以预测复合材料结构的耐久性并且将该新的方法包括在可以用于预测并优化机械结构的耐久性性能的数值设备或工具集中。创新方法以及相关工具集的实用性将使企业能够改进他们用于其轻质产品的设计方法，以较低的成本实现较好的产品质量。此成功的方法的必备要素是：
·在模型中记入累进损伤的能力，即，预测复合材料中的损伤随时间的演化并且在模拟模型中记录实际损伤状态的能力。
·考虑到实际损伤时，预测每个时刻的损伤累积(必须预测朝向下一个步骤的损伤演化)的能力。这需要预测模型中的每个点处的损伤的增加(dD/dN)的所谓的“滞后算子”方法。
·记入可变的循环载荷的影响的能力。这样的时变载荷是在运行(例如，在汽车或航天结构或者风力涡轮机中)经历中的结构，所以必须将时变载荷设计成能承受此运行经历。可变幅值载荷经常由作为块载荷的简化的信号来代替，这是由于简化的信号简化了所建立的模拟和测试。但是如果在设计过程中只使用块载荷(即，随着时间的推移恒定的载荷)，则这不足以确保结构在其生命周期期间可以承受实际的负载。此外，在这些情况中，将与轻质设计的范例相矛盾的大的安全因素增加至设计中。因此，疲劳的模拟不得不也记入可变循环载荷的影响。
有限元(FE)分析(FEA)是众所周知的用于对机械结构进行建模的数值技术。模型包括多个小且简单的“元”。每个这样的元定义相当简单的机械问题，例如其中以许多集成或计算点的形式计算机械方程的板状结构。在建立数值FEA模型之后，可以建立包括彼此相关的大量的简单的“元”问题的完整结构的矩阵方程。解此方程使得能够基于许多小的元问题的解决方案来近似结构的机械特性。要说明的是，本技术的状态是基于对机械FE模型的疲劳预测。FEA的典型结果是基于结构的FE模型机械预测结构中的应力及应变。
已知的一个技术是使用SN曲线以及线性损伤累积用于金属结构的疲劳预测。SN曲线将在发生疲劳破坏之前(即材料损坏之前)施加至材料的应力或载荷循环的幅值表示成循环量的函数。填入SN曲线的数据通常通过测试获得，即，从材料的简单的标准化样本中获得(即，所谓的“样本测试”)。针对所选择的载荷方向(例如，纵向载荷或横向载荷)得出SN曲线。使用雨流计数方法，可以将复杂的载荷循环分成一组具有不同幅值的简单的载荷循环。使用这些幅值，可以然后根据SN曲线得出破坏之前的循环数目。
因此，使用SN曲线预测破坏可以处理可变幅值载荷，但是不能将多轴的情况考虑其中。另外的缺点是，不可以考虑累进损伤的影响，以及因此不可以考虑结构中的刚度的减小和应力的重新分布。这意味着基于雨流的方法只可以解决其中在应力与应变之间存在固定关系的情况下的应力-应变特性及损伤累积，并且相同的应力-应变总是产生相同的损伤。
因此，针对复合材料结构使用SN曲线非常受限。对于可变幅值载荷来说非常典型的是，完成最大载荷循环(其对破坏作出最大贡献)耗费相当长的时间，这是由于其内部存在许多嵌套循环。在这种情况下，当载荷循环(如在雨流方法中)那样被完成时，只考虑循环的方法可能不再是合理的。
传统的疲劳损伤累积使用雨流计数以及线性迈因纳-帕姆格林(Miner-Palmgren)损伤累积。在复合材料的情况下，由于损伤影响被连接至应力-应变滞后回线或分支，所以疲劳特性由于若干不同的破坏机制而随时间变化。这些回线通常彼此嵌套，以使得在一个大的回线闭合之前可以存在在其中打开并闭合的许多小的回线。所以非常典型的是，完成最大载荷循环(其对破坏作出最大贡献)耗费相当长的时间，这是由于其内部存在许多嵌套循环。在这种情况下，当载荷循环如在雨流方法中那样被完成时，只考虑循环的方法可能不再是合理的。
在“Brokate,M；Dressler,K；Krejci,P:Rainflow counting and energydissipation in elastoplasticity,Eur.J.Mech.A/Solids15,.705-737(1996)”中提出了用于损伤累积的方法，该方法可以考虑由于外部变化(例如，温度)所造成的疲劳及材料特性的变化。该方法使用两个滞后算子—一个用于应力-应变特性并且一个用于损伤特性。使用这些具有固定的应力-应变特性的算子以及基于SN曲线的损伤演化，将导致如之前提及的雨流方法及基于SN曲线的方法那样的相同的特性。
在“Nagode,M.&Hack,M.,The damage operator approach,creepfatigue and visco-plastic modelling in thermo-mechanical fatigue,SAEInternational Journal of Materials&Manufacturing,4(1),632-637.doi:10.4271/2011-01-0485,(2011)”中，此方法被扩展至其中应力-应变特性及损伤特性取决于改变应力-应变特性和损伤特性二者的外部作用(温度)。即使此现有技术能够记入外部参数，但是其仍然不能记入如用于复合材料的特性(其中，应力-应变特性及损伤也取决于累积至某点的损伤)。
在相关技术应用“Van Paepegem,W.；Development and finite elementimplementation of a damage model for fatigue of fiber reinforced polymers,Ph.D.thesis,Department of Material Science and Engineering,Ghentuniversity,2002.”中，公开了用于使用N-Jump算法模拟累进损伤以及复合材料疲劳模拟的方法。以刻意选择的间隔使用完整的有限元分析来计算损伤累积及应力重新分布。计算并进一步推算相应的循环中的材料的刚度的降低。此方法的缺点是其只可应用于大载荷并且因此不可以用于可变幅值载荷。
本发明的目的是克服上面所提及的相关技术中的缺点。
由于本发明也可适用于由其他材料制成的结构，所以本申请不意图将本发明的范围唯一地限制成复合材料结构。因此，本发明也适用于具有下述属性的弹性体材料或者任何其他材料：在结构由于疲劳破坏而损坏之前，损伤会由于载荷循环在材料中累积。
发明内容
因此，本申请在第一方面涉及根据权利要求1的方法。
优点是，由于使用滞后算子并且考虑到在载荷循环的完整的幅值之间的时间步的事实，所以可以考虑嵌套回线循环以及多个具有可变幅值的载荷循环。
优点是，由于考虑到材料属性的变化的事实，所以可以预测针对长的载荷循环的疲劳破坏。例如，这些长的载荷循环在虚拟耐久性测试追踪模拟中不表示短的距离，而可以在虚拟测试追踪模拟中表示数千公里。因此，预测将会更精确并且与现实世界中的情况相对应。
另外的优点是，由于考虑到所累积的损伤以及因此材料特性的变化的事实，所以精确地模拟应力的局部重新分布，并且因此，疲劳破坏的预测更精确。
由于对疲劳的更精确的模拟，本发明的另外的优点是需要复合材料结构的较少的超裕度设计。这导致了针对结构的耐久性和性能要求的给定的集合的较轻的结构。
根据第一方面的实施方式，方法还包括根据属性的变化将复合材料结构的模型更新至更新后的模型，并且，针对每个计算点确定由载荷循环引起的新的应力及应变，所述新的应力及应变限定滞后回线。
由于复合材料结构的模型被更新并且计算点中的应力及应变被重新计算的事实，所以本发明的优点是，考虑了应力的全局重新分布以及贯穿复合材料结构的刚度的减小或增加。
本发明的第二方面由根据权利要求8的系统来实现。
在第三方面，本发明由根据权利要求11的软件程序产品来实现。
附图说明
根据下面的描述以及附图将会更好的理解本发明，在附图中：
图1示出了根据本发明的第一方面的方法的优选实施方式；
图2示出了附接至墙壁的梁作为承受外力或载荷的复合材料结构的示例；
图3示出了作为时间的函数的载荷的变化；
图4示出了具有嵌套回线循环的滞后回线；
图5示出了作为相对载荷循环的函数的复合材料的相对刚度；
图6示出了滞后曲线的变化；
图7示出了根据本发明的第一方面的方法的另外的实施方式；以及
图8示出了根据本发明的第二方面的系统的实施方式。
具体实施方式
图1示出了根据本发明的第一方面的优选实施方式的方法的若干步骤。通过此方法，可以预测复合材料结构的耐久性性能。换言之，可以预测复合材料结构在其由于疲劳破坏而损坏之前将能够承受多少载荷循环。
在第一步(1)中，获得其耐久性性能将被预测的复合材料结构的模型。在此第一步中，通过一系列包括一个或更多个计算点的元对该结构完成建模。这可以是有限元(FE)模型。如在技术中众所周知的那样，这样的模型将结构划分成一系列简单的元并且通过简单的方程描述结构的每个元之间的交互，反之，当要考虑作为整体的结构时，必需要对一组复杂的偏微分方程进行求解。此模型通常包括结构的几何属性以及将结构分成一系列这些简单的元的方法。紧接着，例如，模型还包括外力怎样对结构起作用，并且进一步包括结构的材料属性，例如结构的刚度。在图1中，此模型被表示为标称FE模型以指示此模型是在方法的开始处所使用的初始模型。
图2中示出了此结构的说明性示例，其中，梁结构(21)被附接至墙壁(23)。除了保持将梁附接至墙壁的力，梁还承受由重力引起的力Fg(24)及外力F(22)。外力(22)可以例如是由通过附接至梁的另一结构产生的重力或扭力引起的。作为模型，可以与所有这些力所施加的方向的指示一起提供梁的CAD图。基于此CAD图，然后将梁划分成一系列包括计算点的元，并且建立定义元之间的交互的矩阵方程。
当结构承受一系列载荷循环(即，幅值随时间变化的一个或更多个外力)时，将发生复合材料结构中的疲劳破坏。为了预测耐久性性能，方法还需要这些载荷循环作为输入。在图1中，这被表示为载荷历程(2)。图3示出了此载荷(y轴)可以怎样随时间(x轴)以不同的载荷循环(31)进行变化。载荷(通常以牛顿(逐次地，牛顿米)或等同单位来指定的力或力矩)的变化可以不具有重复性也不具有固定的幅值。返回至图2的示例，图3的载荷循环将然后表示外力(22)怎样对梁结构起作用。结构可以承受具有不同的载荷循环的不同的载荷。载荷循环通常将与其中通过模拟或者测量可以获得或使用结构的现实世界情况相对应。例如，当结构是车辆的车轴时，载荷循环将与由于道路的颠簸、拐角、刹车等所引起的力相对应。然后可以通过在此情况期间模拟此情况或通过实际测量对车轴作用的外力来建立载荷历程。
在方法的接下来的步骤(3)中，确定由载荷循环(2)引起的应力及应变(4)。由于应力由载荷循环产生并且包括在结构的模型中，因此应力将沿着结构建立。众所周知，应力用来表达相邻元或彼此施加的连续材料的粒子的内力。通常以兆帕(MPa)或牛顿每平方毫米(N/mm2)来表示应力。可以通过使用以标称FE模型(1)作为输入进行有限元分析(FEA)的有限元(FE)解算器(3)来完成这些应力及应变的计算。FE解算器因此还将也由FE模型(1)定义的外力怎样对结构起作用以及材料属性作为输入。因此，FE解算器将计算所谓的节点处的位移。根据这些计算每个元中的相关应力(用σ表示)及应变(用ε表示)的值。可以针对元中的许多计算点计算这些应力。对于复合材料结构来说，这通常是在不同的板层完成的。在其中计算应力及应变的所有的计算点中可以进行疲劳分析。应变指示由应力所引起的结构的元的位移。因为位移是相对量度，所以其由载荷位移与全局位移的比值来定义。
在图2的示例中，在力F和力Fg的作用下，将会沿着梁建立应力并且引起梁向下弯曲。使用FE解算器计算应力及应变，可以知道梁的每个部分将会怎样弯曲以及其怎样与梁中建立的内部应力相关联。所选择的模型的元越小，应力及应变的计算将会越精确，但是进行计算也将会耗费越多的时间。
方法然后进行至步骤(5)和步骤(6)，其中，预测并存储(8)由于载荷循环引起的累积损伤。现在将更详细地说明这些步骤。
步骤(3)将应力及应变与外部应力或载荷相关联，但是只针对静态情况，即，只针对不变化的外力(22)。然而，为了预测复合材料结构的疲劳，方法还需要考虑由载荷历程(2)所定义的载荷随时间的复杂变化。在结构承受此载荷循环时，从FE解算器获得的结构的每个计算点中的应力与应变的关系将展示出记忆效应。换言之，应变不只取决于材料，还取决于应力的历程，即，如何建立到应力的路径。图4中示出了复合材料结构的计算点中应力及应变的变化以及应力及应变从初始时刻t0直到时刻t6如何变化。图3的时间戳t0-t6与图4的时间戳t0-t6相对应。图4中的路径t0-t6因此示出了复合材料结构的计算点中的应力σ及应变ε在变化载荷即图3的载荷循环的影响下将怎样变化。图3的载荷循环将遵循图4所表示的滞后曲线、回线或分支(41)(42)。由于载荷随时间的周期变化是复杂的，所以滞后曲线或回线是嵌套的。换言之，新的回线(42)例如t3-t4-t5将在另一回线(41)例如t2-t6结束之前开始。因此，由载荷循环所引起的应力及应变限定了滞后回线(42)(41)。
在步骤(5)中，使用滞后算子(5)，作为载荷随时间变化的函数现在计算应力(7)沿着滞后回线的一部分的变化。基于载荷历程(2)以及FE解算器(3)的输出，滞后算子将计算与下一个时间步t1(见图3和图4)相关的应力σ1(7)(见图4)。可以在“Nagode,M.&Hack,M.,Thedamage operator approach,creep fatigue and visco-plastic modelling inthermo-mechanical fatigue,SAE International Journal of Materials&Manufacturing,4(1),632-637.doi:10.4271/2011-01-0485,(2011)”中发现滞后算子的可能的实现方式。因此，计算应力沿着滞后回线的一部分的变化作为载荷随时间变化的函数。由于复合材料结构中的滞后影响，滞后算子不只考虑内部为单个回线的情形，还要考虑回线是否已被嵌套至另一回线中的情形。为了实现此目的，滞后算子(5)追踪材料的记忆。在步骤(5)中，内部状态变量(8)还将追踪方法在方法的步骤(5)的每个执行中恰好沿着滞后回线。状态信息因此追踪材料记忆。换言之，滞后算子(5)包括存储滞后回线(41)中的位置信息作为状态信息。
返回至图2的梁的示例，滞后算子(5)将针对某个值的所施加的应力F(22)计算某时刻t1处的梁的每个计算点的应力。当施加至梁的力以零开始(未施加力)时，力然后在时刻t3时变成最大并且在时刻t6时返回零(再次未施加力)，残余的应力将不返回至零，并且因此，即使未施加力时梁也将保持弯曲至某程度。这示出了梁的材料的记忆效应。
在另外的步骤(6)中，使用根据步骤(5)的应力(7)的变化以及所存储的累积损伤，计算损伤的变化(6)，并且因此，还计算每个计算点中的包括刚度的结构的属性(14)的变化。由于方法是重复的过程，因此当第一次执行步骤(6)时，累积损伤的值将具有初始预定值，例如零。在复合材料结构中，每个载荷循环(31)将在结构内部产生损伤的变化。在一定量的循环之后，损伤的所有的变化的总和(即累积损伤)将超过某一阈值并且将发生疲劳破坏。在复合材料结构中，可以如下面的等式(eq)那样来表征累进损伤或损伤的变化：
换言之，由载荷循环的某载荷循环(N)引起的材料的某方向l中的损伤的变化取决于已经由载荷循环引起的总的累积损伤(dl)及应力σ。应力反过来还取决于已经累积的损伤dl。图5示出了针对复合材料的预定的应力值以相对于载荷循环的相对数目(N/Nf)绘出相对刚度的变化(E/E0)的情况。如附图中明显的示出的那样，刚度的变化取决于载荷循环。刚度反过来可以通过下面的等式(eq)进一步与累积损伤相关联：
其中
El是材料在方向l中的刚度
El 0是在施加载荷之前的原始刚度
dl是总共的累积损伤
下面的等式(eq)与图5的曲线相关，并且损伤与刚度的上面的关系然后产生如在步骤6中所计算的累进损伤：
其中
是相对应力值，以及
X1是指示结构在静态应力值下将会被永久地破坏的静态强度常数，
c1，c2，c3，c4及c5是与通过材料样板的应力测量获得的与材料相关的常数，以及
dl是结构在方向l中的总共的累积损伤。
上面的方程只是关于如何导出结构的属性以及损伤的变化的示例，但本发明不限于此。通常，计算损伤的变化可以直接使用如图5所示出的特征数据，该特征数据将属性和损伤的变化与累积损伤以及应力的变化相关联。或者，如【eq3】所示，计算损伤的变化还可以包括使用用于获得损伤的变化的解析式。解析式因而是特征数据的近似。
根据刚刚获得的累进损伤以及累积损伤的先前值导出累积损伤的新的值。此新的值然后被存储作为状态信息(8)以作为材料属性的一部分直到过程进行下一个重复为止。相似地，计算作为材料属性的刚度El的新的值。
参考图2的梁的示例，在方法的步骤(6)中，可以针对梁的每个计算点采用公式计算累进损伤并且记入获得的应力(7)的值、根据状态信息的已经累积的损伤以及材料常数并获得损伤的变化的值。根据损伤的变化，然后将状态信息(8)中的累积损伤更新至新的值。
由于材料刚度因为材料的损伤累积已经变化，所以如图4中所描绘的滞后曲线(41)不再更新。复合材料结构的每个计算点中的应力与应变之间的关系已经变化并且需要重新计算。因此，在步骤(15)中，基于属性(14)的变化计算(15)应力及应变(16)的进一步的变化(15)以确定(7)新的适合的滞后分支。应当注意的是，由于应力-应变特性的连续变化，回线不再需要闭合，因此回线也可以被注释为分支以指示它们未闭合的事实。换言之，因为刚度El的变化，所以应力及应变在结构的模型的计算点中局部地重新分布。就物理意义而言，刚度值El是滞后曲线的斜率的量度。由于在应力的重新计算中未考虑复合材料结构的元之间的交互，所以应力的重新分布被认为是局部的。如图6中所示，应力及应变(16)的进一步的变化将限定适合的滞后分支(61)。
参考图2的梁的示例，更新应力及应变(16)将意味着梁将或多或少弯曲。这甚至在载荷历程(2)的不同的循环中施加同一力时也会是此情况。例如，当梁中的元的刚度减小时，梁将会在下一个载荷循环中对于同一施加的力而弯曲的更多。
在步骤(15)之后，方法返回至步骤(5)，其中，计算应力(7)沿着所述适合的滞后分支的另外的部分的进一步的变化作为载荷随时间的进一步的变化的函数。参考图3和图4，这意味着计算复合材料结构的计算点中的应力在下一个时间步t2处的新值。图6示出了新的滞后曲线(61)导致时间步t2处的新的应力值σ2’。如果不使用新的材料刚度值更新曲线，则将会获得另一错误的应力值σ2。为了能够沿着滞后曲线进行，在计算应力(7)沿着适合的滞后分支的一部分的变化时需要在每个步骤(5)中使用所存储的关于包括滞后曲线中的当前位置的材料记忆的信息。
根据时间步t2处的新的应力值σ2’，方法再次进行至步骤(6)，其中，如之前所描述的那样获得累进损伤以及刚度的新的值。然后，继续沿着步骤(5)、(6)及(15)的循环，同时连续更新状态信息(8)，即，材料记忆(例如，滞后分支位置)以及材料的属性的变化(例如，刚度及累积损伤)。在某时刻tn，在载荷循环(2)中，复合材料结构的计算点中的一个计算点将超过某损伤累积值(以及从而相关刚度值)，并且因此，将发生疲劳破坏。在那时，因为已知疲劳破坏将发生在载荷历程(2)中的哪个时刻，所以复合材料结构的耐久性性能是已知的。在那时，所获得的复合材料结构的疲劳特性可以是符合要求的并且该结构可以通过使用如由标称FE模型(1)所指定的材料属性及尺寸来制造。如果不符合复合材料结构的要求，则可以调整模型并且重复预测耐久性性能。复合材料结构的耐久性性能的方法优化可以被获得。
参考梁的示例，当进行重复的步骤(5)、(6)及(15)时，在某时刻tn处，在梁的某计算点中，损伤累积将会使刚度过小。换言之，梁将会根据需求或者其已被折断而弯曲的更厉害。
图7示出了改进根据本发明的第一方面的方法的精度的另外的步骤。如之前所指出的那样，损伤的累积导致复合材料结构的元中的刚度的局部下降。这还会导致复合材料结构中的应力的全局重新分布。这只可以通过新的有限元计算来说明。因此，当进行应力的全局重新分布时，在步骤(12)中，使用属性的变化将结构的模型更新为更新后的模型(13)，并且针对每个计算点确定由载荷循环(2)引起的新的应力及应变(4)，所述新的应力及应变(4)限定滞后分支。在进行步骤(12)之前，确定(9)是更新结构的模型(12)还是计算应力及应变(16)的另外的变化。
当确定重新进行FEA分析时，使用状态信息(8)保存系统的状态：
保存所有的局部损伤累积状态；
保存滞后算子的所有的内部状态(材料记忆)；以及
保存载荷历程(2)中的时间。
然后，将材料属性(例如刚度)的变化应用至有限元模型(12)，并且因此，获得损伤的FE模型以及更新后的模型(13)。然后，通过FE解算器(3)计算更新的应力及应变(4)并且滞后算子在其停止处重新开始。
根据本发明的第二方面，本发明还涉及用于进行根据本发明的第一方面的方法的系统。图8中示出了此系统的实施方式。
系统包括用于通过一系列包括一个或更多个计算点的简单元来对复合材料结构进行建模的装置(81)。这可以由其中使用计算机的输入装置(鼠标、键盘……)绘制复合材料结构的模型的计算机来实现。使用此计算机，通过定义元及FE模型、外力以及材料属性来进一步完成模型。此外，装置(81)还可以包括将结构的物理样本数字化的浏览装置以及在其上进一步完成FE建模的计算机。
系统还包括用于确定载荷历程(2)及载荷循环的装置(82)。前面提及的装置可以包括另一计算机或者与包括在装置(81)中的那样的同一计算机，以通过模拟确定载荷循环。装置(82)还可以包括测试装置以通过实际物理实验测量载荷循环。
系统还可以包括用于计算根据本发明的第一方面的方法的步骤(3)、(5)、(6)、(9)、(12)及(15)的装置(83)。可以以程序语言写入这些步骤并且编译这些步骤以在包括在计算装置(83)中的处理器上运行这些步骤。装置(83)还可以包括用于追踪状态信息(8)的存储器(86)。此存储器可以是计算装置(83)的RAM存储器或硬盘驱动器。
系统还包括用于相应地制造复合材料结构的装置(84)。
Claims (14)
1.一种用于虚拟地预测结构的耐久性性能从而能够优化所述耐久性性能的方法，其中，通过一系列计算点来对所述结构进行建模，其中，
针对每个计算点，
-确定(3)由载荷循环(2)所引起的应力及应变(4)，所述应力及应变(4)限定滞后分支，以及
-预测(5)(6)并存储(8)由于所述载荷循环而引起的累积损伤，以及
-使用滞后算子(5)，作为所述载荷随时间变化的函数计算所述应力(7)沿着滞后分支的一部分的变化，以及
-使用所述应力(7)的所述变化以及所存储的所述累积损伤，计算(6)所述损伤的变化，并且因此，还计算所述结构的属性(14)在所述计算点处的变化，所述属性(14)包括刚度，以及
-基于所述属性(14)的所述变化来计算(15)所述应力及应变(16)的进一步的变化，以确定新的适合的滞后分支，以及
-作为所述载荷随时间进一步变化的函数，计算(5)所述应力(7)沿着所述适合的滞后分支的另外的部分的进一步的变化，以及
-在方法的最后，相应地制造所述结构。
其中，所述结构由下述材料制成：当经受所述载荷循环时，由于疲劳破坏，在损坏之前累积损伤的材料。
2.根据权利要求1所述的方法，其中，确定由载荷循环引起的所述应力及应变还包括对所述结构的有限元模型(1)进行有限元分析，所述有限元分析基于由所述载荷循环施加在所述结构上的外力来进行。
3.根据权利要求1或2所述的方法，还包括使用所述属性的变化将所述结构的所述模型更新(12)为更新后的模型(13)，并且，针对每个计算点，确定由所述载荷循环(2)所引起的新的应力及应变(4)，所述新的应力及应变(4)限定滞后分支。
4.根据权利要求3所述的方法，还包括：在所述更新(12)之前，确定(9)是进行所述更新(12)还是进行(15)对所述应力及应变(16)的进一步的变化的所述计算。
5.根据权利要求1至4所述的方法，其中，使用所述滞后算子(5)还包括将所述滞后分支的位置的信息存储为状态信息(8)，并且其中，计算所述应力(7)沿着所述适合的滞后分支的一部分的变化使用存储的信息。
6.根据权利要求1至5所述的方法，其中，计算所述损伤的变化还包括使用特征数据，所述特征数据将损伤和属性的所述变化与所述累积损伤以及所述应力的所述变化相关联。
7.根据权利要求6所述的方法，其中，计算所述损伤的所述变化还包括使用用于获得所述损伤的所述变化的解析式，所述解析式是所述特征数据的近似。
8.一种用于虚拟预测结构的耐久性性能从而能够优化所述耐久性性能的系统，所述系统包括：
-用于通过一系列计算点来对所述结构进行建模的装置(81)，以及
-用于针对每个计算点，确定由载荷循环所引起的应力及应变的装置(83)，所述载荷循环限定滞后分支，以及
-用于确定所述载荷循环的装置(82)，以及
-用于针对每个计算点，预测由于所述载荷循环而引起的累积损伤的装置(83)，以及
-用于针对每个计算点，存储所述累积损伤的装置(86)，以及
-用于使用滞后算子，作为所述载荷随时间变化的函数，针对每个计算点计算所述应力(7)沿着所述滞后分支的一部分的变化的装置(85)，以及
-用于使用所述应力(7)的所述变化以及所存储的累积损伤，针对每个计算点计算所述损伤的变化的装置(85)，以及因此，用于针对每个计算点计算所述结构的属性在所述计算点处的变化的装置，所述属性包括刚度，以及
-用于基于所述属性(14)的所述变化，针对每个计算点计算所述应力及应变的进一步的变化，以确定新的适合的滞后分支(7)的装置(85)，以及
-用于针对每个计算点，作为所述载荷的进一步的变化的函数，计算所述应力(7)沿着所述适合的滞后分支的一部分的变化的装置(85)，以及
-用于相应地制造所述结构的装置(84)。
9.根据权利要求8所述的系统，其中，用于确定由载荷循环引起的所述应力及应变的装置还包括用于对所述结构的有限元模型(1)进行有限元分析的装置，所述有限元分析基于由所述载荷循环施加在所述结构上的外力来进行。
10.根据权利要求8或9所述的系统，还包括：用于使用所述属性的变化对所述结构的所述模型进行更新(12)的装置，以及用于针对每个计算点，确定由所述载荷循环(2)所引起的新的应力及应变(4)的装置，所述新的应力及应变(4)限定滞后分支。
11.一种用于执行权利要求1的所述方法的软件步骤(3)、(5)、(6)及(15)的计算机程序产品。
12.根据权利要求11所述的计算机程序产品，进一步执行权利要求3的所述方法的软件步骤(12)。
13.根据权利要求1至7中任一项所述的方法，其中，所述材料是由复合材料或弹性体材料构成的组中的任何一种材料。
14.根据权利要求8至10中任一项所述的系统，其中，所述材料是由复合材料或弹性体材料构成的组中的任何一种材料。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| EP13167410.3 | 2013-05-13 | ||
| EP13167410.3A EP2803968B1 (en) | 2013-05-13 | 2013-05-13 | A process for calculating fatigue and fatigue failure of structures |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104217094A true CN104217094A (zh) | 2014-12-17 |
| CN104217094B CN104217094B (zh) | 2019-02-19 |
Family
ID=48470735
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410201290.6A Active CN104217094B (zh) | 2013-05-13 | 2014-05-13 | 用于计算结构的疲劳以及疲劳破坏的方法 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US20140336954A1 (zh) |
| EP (1) | EP2803968B1 (zh) |
| CN (1) | CN104217094B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106383052A (zh) * | 2016-08-19 | 2017-02-08 | 北京工业大学 | 一种考虑夹杂的金属疲劳载荷下薄弱区域确定方法 |
| CN110573402A (zh) * | 2017-04-28 | 2019-12-13 | 通用电气全球采购有限责任公司 | 检测推进子系统退化的监控系统 |
| CN111091310A (zh) * | 2020-01-17 | 2020-05-01 | 中联重科股份有限公司 | 挖掘设备健康监测系统及方法 |
| CN111122358A (zh) * | 2020-01-13 | 2020-05-08 | 上海工程技术大学 | 一种考虑滞弹性能的镁合金疲劳寿命的确定方法 |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150274314A1 (en) * | 2014-03-28 | 2015-10-01 | Bell Helicopter Textron Inc. | Aircraft Prognostics Health System |
| CN104573172B (zh) * | 2014-11-19 | 2017-10-31 | 新疆金风科技股份有限公司 | 风力发电机组中的结构件的疲劳分析方法和疲劳分析装置 |
| PT3237873T (pt) * | 2014-12-23 | 2019-04-01 | Ore Catapult Development Services Ltd | Ensaio à fadiga |
| US10378426B2 (en) | 2016-12-16 | 2019-08-13 | Caterpillar Inc. | Useful life prediction system for cooling components |
| US10691847B2 (en) * | 2017-01-13 | 2020-06-23 | Sap Se | Real-time damage determination of an asset |
| CN107817081B (zh) * | 2017-08-31 | 2019-12-20 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种座舱盖静力试验台刚度分布模拟方法 |
| CN109885874B (zh) * | 2019-01-11 | 2022-12-23 | 华东理工大学 | 一种基于abaqus的多轴蠕变疲劳预测方法 |
| CN110414161B (zh) * | 2019-08-01 | 2023-04-07 | 广州建筑股份有限公司 | 一种梭形滞回曲线的参数确定方法及系统 |
| CN112560162B (zh) * | 2019-09-24 | 2024-05-10 | 上海汽车集团股份有限公司 | 一种动力总成悬置路谱载荷的缩减方法及装置 |
| CN113396437B (zh) * | 2019-09-30 | 2024-04-09 | 日立建机株式会社 | 疲劳管理系统 |
| CN111474044B (zh) * | 2020-04-24 | 2021-03-19 | 南京航空航天大学 | 一种考虑多级疲劳损伤的编织陶瓷基复合材料疲劳迟滞回线的预测方法 |
| CN113378327B (zh) * | 2021-07-02 | 2022-09-02 | 吉林重通成飞新材料股份公司 | 风电叶片后缘抗开裂设计方法、装置、设备及存储介质 |
| CN113792414B (zh) * | 2021-08-20 | 2024-03-29 | 北京工业大学 | 一种变幅加载下陶瓷基复合材料疲劳寿命的预测方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2012118059A (ja) * | 2010-11-09 | 2012-06-21 | Nippon Steel Corp | 部品の疲労破壊評価装置、部品の疲労破壊評価方法、及びコンピュータプログラム |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5419200A (en) * | 1994-06-14 | 1995-05-30 | The United States Of America As Represented By The Secretary Of The Army | Method for assessing the effects of loading forces on a composite material structure |
| WO2010080196A2 (en) * | 2009-01-08 | 2010-07-15 | Battelle Memorial Institute | Path-dependent cycle counting and multi-axial fatigue evaluation of engineering structures |
| US20140156229A1 (en) * | 2012-12-03 | 2014-06-05 | Caterpillar Inc. | Fatigue-based topology optimization method and tool |
-
2013
- 2013-05-13 EP EP13167410.3A patent/EP2803968B1/en active Active
-
2014
- 2014-05-09 US US14/273,781 patent/US20140336954A1/en not_active Abandoned
- 2014-05-13 CN CN201410201290.6A patent/CN104217094B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2012118059A (ja) * | 2010-11-09 | 2012-06-21 | Nippon Steel Corp | 部品の疲労破壊評価装置、部品の疲労破壊評価方法、及びコンピュータプログラム |
Non-Patent Citations (5)
| Title |
|---|
| MARKO NAGODE等: ""The Damage Operator Approach: Fatigue, Creep and Viscoplasticity Modeling in Thermo-Mechanical Fatigue"", 《SAE INTERNATIONAL JOURNAL OF MATERIALS AND MANUFACTURING》 * |
| S.SUBRAMANIAN等: ""A cumulative damage model to predict the fatigue life of composite laminates including the effect of a fibre-matrix interphase"", 《INTERNATIONAL JOURNAL OF FATIGUE》 * |
| 刘灵灵等: ""多轴加载下缺口件弹塑性有限元分析及疲劳寿命预测"", 《机械设计》 * |
| 胡青春等: ""基于ANSYS Workbench的注胚模腔疲劳寿命研究"", 《塑料工业》 * |
| 苏承东等: ""砂岩巴西劈裂疲劳破坏过程中变形与强度特征的试验研究"", 《岩石力学与工程学报》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106383052A (zh) * | 2016-08-19 | 2017-02-08 | 北京工业大学 | 一种考虑夹杂的金属疲劳载荷下薄弱区域确定方法 |
| CN110573402A (zh) * | 2017-04-28 | 2019-12-13 | 通用电气全球采购有限责任公司 | 检测推进子系统退化的监控系统 |
| CN111122358A (zh) * | 2020-01-13 | 2020-05-08 | 上海工程技术大学 | 一种考虑滞弹性能的镁合金疲劳寿命的确定方法 |
| CN111122358B (zh) * | 2020-01-13 | 2022-05-31 | 上海工程技术大学 | 一种考虑滞弹性能的镁合金疲劳寿命的确定方法 |
| CN111091310A (zh) * | 2020-01-17 | 2020-05-01 | 中联重科股份有限公司 | 挖掘设备健康监测系统及方法 |
| CN111091310B (zh) * | 2020-01-17 | 2020-12-08 | 中联重科股份有限公司 | 挖掘设备健康监测系统及方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| EP2803968B1 (en) | 2019-05-08 |
| EP2803968A1 (en) | 2014-11-19 |
| US20140336954A1 (en) | 2014-11-13 |
| CN104217094B (zh) | 2019-02-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104217094A (zh) | 用于计算结构的疲劳以及疲劳破坏的方法 | |
| Lin et al. | Damage detection in the cable structures of a bridge using the virtual distortion method | |
| Kong et al. | Fatigue life prediction of parabolic leaf spring under various road conditions | |
| Kong et al. | Mission profiling of road data measurement for coil spring fatigue life | |
| Ambassa et al. | Fatigue life prediction of an asphalt pavement subjected to multiple axle loadings with viscoelastic FEM | |
| Li et al. | A unified frequency domain fatigue damage modeling approach for random-on-random spectrum | |
| Wu et al. | Development of sensor placement optimization tool and application to large-span cable-stayed bridge | |
| Wang et al. | Life prediction for aircraft structure based on Bayesian inference: towards a digital twin ecosystem | |
| Su | A road load data processing technique for durability optimization of automotive products | |
| CN110955933B (zh) | 一种基于响应面法的机械结构模糊疲劳可靠度计算方法 | |
| CN114491825B (zh) | 一种汽车翼子板振动强度分析方法 | |
| Brudnak et al. | Durability test time reduction methods | |
| Raphael et al. | Sensor data interpretation in bridge monitoring—A case study | |
| Du et al. | Fatigue life prediction of the gear box in tracked vehicles based on running simulation tests | |
| Kulisiewicz | Modeling and identification of nonlinear mechanical systems under dynamic complex loads | |
| Rios et al. | Synthetic data generation for the creation of bridge digital twins what-if scenarios | |
| Nabiyan et al. | Structural identification for dynamic strain estimation in wind turbine towers | |
| Jakobsson et al. | Data driven modeling and estimation of accumulated damage in mining vehicles using on-board sensors | |
| Lin | A New Method of d'Alembert's Principle Finite Element Based Fatigue Calculation with Input of Loads and Accelerations | |
| Aflatooni et al. | Condition monitoring and rating of bridge components in a rail or road network by using SHM systems within SRP | |
| Stephens et al. | Computational fatigue life prediction of welded and non-welded ground vehicle components | |
| Saleem et al. | Health monitoring of ultra high fiber performance reinforced concrete communication tower using machine learning algorithms | |
| Fedorchenko et al. | Cycle Counting Methods of the Aircraft Engine. | |
| Papadioti et al. | Fatigue monitoring in metallic structures using vibration measurements | |
| Ceravolo et al. | Strategy for the maintenance and monitoring of electric road infrastructures based on recursive lifetime prediction |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C53 | Correction of patent of invention or patent application | ||
| CB02 | Change of applicant information |
Address after: Holland Brabant Applicant after: SIEMENS INDUSTRY SOFTWARE LTD. Address before: Holland Brabant Applicant before: LMS INTERNATIONAL NV |
|
| COR | Change of bibliographic data |
Free format text: CORRECT: APPLICANT; FROM: LMS INTERNATIONAL N. V. TO: SIEMENS INDUSTRY SOFTWARE NV Free format text: CORRECT: ADDRESS; FROM: |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |