结构损伤一致多尺度分析的关键是在不同尺度的模型衔接起来进行计算，为此首先需要将前述跨尺度界面单元衔接方法引入ABAQUS，按软件指定的格式引入跨尺度界面上须满足的位移约束方程。但是由于K2本身是奇异的，同时计算机有效位数是有限的，α过大会导致系统方程病态而使计算失效。根据所建立的结构一致多尺度模型并通过宏细观变量的跨尺度关联，便可实现结构损伤的并发多尺度计算。......
2023-08-26
结构损伤多尺度模拟：裂纹扩展准则3.4.2
【摘要】：图3.15疲劳裂纹扩展速率分区图3.15中I区的特性是存在一个“疲劳门槛”，即循环应力强度因子幅值ΔKth，在该区内应力强度因子幅值低于门槛值，疲劳裂纹基本不扩展。在高速率区，疲劳裂纹k扩展速率较快，出于简化计算和突出对II区的研究，认为裂纹扩展速率恒为，进入高速率区的判断准则为该裂纹的应力强度因子幅值ΔKth，J＞0.8 KC，J。
主要考虑两种工程中常见荷载工况下的裂纹扩展准则，一是准静态加载，另一是疲劳荷载。准静态加载导致的破坏是由材料达到强度后不能抵挡所施加的静态荷载而导致的；而疲劳是在低于或者远低于强度的荷载循环作用下，材料中的疲劳裂纹萌生和不断扩展导致的破坏，两种情况具有不同的损伤与失效机理。
1.准静态裂纹扩展准则
在准静态荷载下，裂纹扩展的驱动力为应力强度因子。根据线弹性断裂力学理论，当裂纹尖端的应力强度因子达到断裂韧度时发生裂纹扩展，即裂纹扩展的准则是：
这里KI为I型裂纹的应力强度因子，KI，C为I型断裂韧度。
基于上述准则，任意扩展步后计算得到所有裂纹尖端的应力强度因子，其中第J条裂纹的I型（拉伸）和II型（剪切）应力强度因子为
。那么第J条裂纹的等效应力强度因子可以计算为：
其裂纹扩展准则为：
这里KC，J为第J条裂纹尖端的断裂韧度（toughness），根据下面不同情况选取不同值：
（1）在界面上：KC，J＝KC，interface
（2）在砂浆中：KC，J＝KC，mortar
在每一步分析中，首先施加测试载荷F0或d0（分别对应于应力控制和位移控制的边界条件）。求解该问题可以获得所有裂纹尖端的应力强度因子KJ，然后确定“危险”（即处于扩展临界状态）的裂纹尖端，β为接近断裂韧度的裂纹，即
这里Ntips为所有裂纹尖端总数。
具有“危险”的裂纹J的扩展方向由以下方法确定：
（1）若尖端位于界面上，裂纹将沿着界面扩展，迅速地扩展到长轴的两端。在后面的计算步中，裂纹将被视为一般裂纹，即断裂韧度取为KC，mortar。
（2）若尖端位于砂浆，裂纹将往指定方向扩展Δl长度，扩展方向由最大周向拉应力原则确定，即
2.疲劳裂纹扩展准则
疲劳荷载下的裂纹扩展规律由Paris提出，建立了应力强度因子和裂纹扩展速率之间的关系，是当今工程中应用广泛的疲劳裂纹扩展理论：
其中a为裂纹长度，da／dN为裂纹扩展速率，C、m为材料常数，可由实验数据拟合得到，ΔK为应力强度因子幅值。
疲劳裂纹扩展的三阶段理论引入了疲劳扩展的应力强度因子门槛值ΔKth，认为低于此门槛时裂纹不扩展。疲劳裂纹扩展的三个阶段（区）如图3.15所示。
图3.15 疲劳裂纹扩展速率分区
图3.15中I区的特性是存在一个“疲劳门槛”，即循环应力强度因子幅值ΔKth，在该区内应力强度因子幅值低于门槛值，疲劳裂纹基本不扩展。II区为应力强度因子幅值大于ΔKth的疲劳裂纹扩展特性，在该区内的裂纹扩展速率da／dN与应力强度因子幅值ΔK的关系服从Paris公式。在III区，ΔKmax已经接近材料的断裂韧度KI，C，裂纹扩展速率将呈失稳状态急剧增快直至迅速断裂。
实验研究［17］表明，混凝土的裂纹扩展与加载的频率相关，当加载频率较低时，混凝土构件裂纹长度也会像金属材料那样随着疲劳荷载次数的增加而缓慢增加，并且扩展速率也服从Paris公式。根据上述理论，对于疲劳荷载，可引入以下的裂纹扩展率计算方式：
（1）在低速率区（I区），疲劳裂纹扩展速率很低，认为裂纹不发生扩展，该裂纹的扩展门槛值ΔKth，J＝0.2 KC，J。其中J表示不同的裂纹，KC，J为该条裂纹的断裂韧度。
（2）在中速率区（II区），认为裂纹扩展速率符合满足Paris公式：
其中混凝土材料参数C、m的选取参考文献［17］的实验数据，m取值一般为1.0～2.5之间，裂纹扩展速率一般在10－5 mm／c与10－2 mm／c之间。
（3）在高速率区（III区），疲劳裂纹k扩展速率较快，出于简化计算和突出对II区的研究，认为裂纹扩展速率恒为
，进入高速率区的判断准则为该裂纹的应力强度因子幅值ΔKth，J＞0.8 KC，J。
需要说明的是，实际计算过程中，一个计算步不等于一个荷载循环步。疲劳模拟中不会把每一次荷载循环都逐个计算，这会导致计算效率很低，实际上，荷载只有循环到一定次数之后，裂纹才出现明显的扩展。因此，计算中实行的是控制Δa策略，即
在每一个“计算步”（记为i上标，以区别于表示裂纹的下标J）中，控制裂纹扩展速率最大的那个裂纹（记为β）所扩展的长度为Δai，然后计算出ΔNi，然后遍历所有的裂纹，计算裂纹扩展量
。最后对ΔNi求和得到总寿命：
有关结构损伤多尺度模拟与分析的文章
- 详细阅读
-
结构损伤多尺度模拟：含裂纹的脆性固体模型详细阅读
为了简化起见，假设混凝土中的微裂纹是三维钱币状的。在建立内含大量分布微观裂纹固体的细观力学模型以前，需要首先讨论孤立裂纹的情况，然后再拓展到含大量分布微裂纹的无限大固体。......
2023-08-26
-
结构损伤多尺度模拟-跨尺度演化规律详细阅读
目前对在疲劳损伤过程中微裂纹成核与扩展问题的研究，大都是在单一尺度下进行。疲劳损伤累积过程中的微裂纹成核与扩展阶段的疲劳损伤过程是一个跨尺度演化行为，在细观尺度下，表现为微裂纹的成核与扩展造成的微观缺陷的增长，在宏观尺度下，表现为微裂纹扩展的群体效应造成宏观疲劳损伤累积和寿命消耗。......
2023-08-26
-
结构损伤多尺度模拟：分形概念详细阅读
分形维数可对此类自相似现象的整体特性提供数学描述。分形研究的对象是具有不规则性和自相似性的无序系统，定量描述这种自相似性的参数是分形维数，分形维数的变化是连续的。可以用分形维数来表征广泛存在于自然界的一类无序、复杂、奇异客体。......
2023-08-26
-
结构损伤多尺度模拟与分析详细阅读
对这些易损部位就需要考虑材料与结构在细观尺度下存在的细观缺陷及其非连续、非均匀特征进行损伤分析，在这些易损局部的复杂几何构造细节和在内部客观存在的细观缺陷，导致局部热点应力升高、触发损伤演化。分析其组成材料的性能、含量、界面、缺陷等，基于对这些细观构造与缺陷的细观分析结果来确定这些易损部位在宏观尺度下结构模型中的近似于均匀材料的等效性能。......
2023-08-26
-
结构损伤多尺度模拟案例-数值模拟成果详细阅读
本节将展示一个典型的混凝土试样由微裂纹扩展导致破坏的数值模拟案例，以说明所提出的混凝土细观多裂纹模型可以有效地模拟混凝土中裂纹扩展导致构件破坏的全过程。图3.20混凝土单轴拉伸试样对任何一种混凝土配比，将随机产生8～10个数值分析的试样。图3.21微裂纹扩展导致试样破坏的全过程模拟如图3.21是试样初始状态下的裂纹分布状。......
2023-08-26
-
结构损伤多尺度模拟与分析验证详细阅读
在考虑材料损伤演化与失效过程的有限元分析中，损伤演化方程与结构平衡方程、几何方程、物理方程一起，构成基本方程组。目前耦合损伤的有限元计算方法主要有：全解耦方法。本章在上一章对两种多尺度方法基本理论介绍的基础上，分别阐述了两种多尺度分析方法在有限元计算软件ABAQUS中实施的算法流程与具体步骤，并通过算例分析以验证两种损伤多尺度分析方法各自的有效性和可行性。......
2023-08-26
-
结构损伤的多尺度模拟与分析详细阅读
因此，结构损伤的量化其实已转化为与特定结构构成和内力分配关系相关的求解过程。对于不同类型的结构，有必要针对其特定的构造形式专门研究其结构损伤的量化方法。......
2023-08-26
相关推荐