灵活裂纹裁减模拟 (活动裂缝)

本文目录导航：

• 灵活裂纹裁减模拟
• 各类钢材断裂要素汇总

灵活裂纹裁减模拟

一、数值流形方法的能源学求解格局

前面所探讨的裂纹裁减疑问，都只限于静态加载或缓慢加载下的断裂疑问和裂纹的稳态裁减疑问，它们都属于断裂静力学的领域。

但是在实践工程中，会经常遇到含裂纹的结构在高加载速率下的断裂疑问或裂纹在失稳后的极速裁减疑问，而这两类疑问都属于断裂能源学识题。

断裂能源学是钻研惯性效应不能疏忽的断裂力学识题（范天佑，1990），断裂能源学识题的求解方法清楚不同于断裂静力学识题，在灵活加载时，试件除发生弹塑性变形外，外部各质点的自在振动取得必定的减速度，从而发生惯性力，即所谓的灵活加载时的惯性效应（Chen E P et al.，1976）。

数值流形方法在解决能源学识题时，在以后步的计算中，各单元承袭了前一期间步的速度，而不是像解决静力学识题一样置以后期间步的速度为零。

这种解决方法思考了单元的惯性效应，但没有思考单元的阻尼，对普通能源学识题精度不够，本节在数值流形方法的能源学求解格局剖析的基础上，启动改良，以更好地模拟普通能源学识题。

在数值流形方法的求解公式中，能源疑问与静力疑问的关键区别在于品质矩阵，在能源学识题中，这是一个很关键的矩阵，当期间步比拟小时，惯性力矩阵将管理着整个资料体一切各点的静止和稳固。

思考以后期间步时，（u（x，y，t）v（x，y，t））示意为单元e的任一点（x，y）与期间相关的位移，假定M示意单位面积的品质，e是q个笼罩（U，U，…，U）的交加。单位面积上的惯性力为

岩石断裂与挫伤

式中：

［T（x，y）］=（TTT…T）

｛D（t）｝=｛D（t）D（t）D（t）…D（t）｝

在单元e中，由惯性力惹起的势能为

岩石断裂与挫伤

设{D}={0}是在期间步起始时的单元位移，{D（Δ）}={D}是在期间步结束时的位移，Δ是期间步长。则

岩石断裂与挫伤

所以

岩石断裂与挫伤

式中：，是期间步开局时的单元速度。该期间步结束时的速度{V（Δ）}是

岩石断裂与挫伤

则势能变为

岩石断裂与挫伤

上式可合成为两局部，第一局部对应于刚度矩阵，为

岩石断裂与挫伤

第二局部对应于荷载向量，为

岩石断裂与挫伤

其中：r，s=1，…，12。

由上可知，数值流形方法对能源学识题的求解是显式算法，即应用以后期间步的位移和速度求下一时步的位移和速度，在每一期间开局时，承袭了前一时步的速度和位移，且位移是用泰勒级数倒退的，其截断误差是三阶的，当期间步取得较大时，误差也较大。另外该算法也是有条件的，在计算环节中期间步长的选用必定要小于某一个值，即（王勖成等，1995）：，其中：T为求解系统中最小尺寸单元的最小固有振动周期，因此求解系统中最小单元的尺寸将选择其期间步长的选用。

这种算法比拟适宜于动摇疑问的求解，一方面是由于这种方法的求解特点正好和波的流传特点相分歧；另一方面，钻研波的流传环节要求驳回小的期间步长，也合乎该算法的要求。

普通的结构能源学识题结构的能源相应中低频成分占主导，从计算精度上思考，准许驳回较大的期间步长，通常驳回Newmark无条件稳固隐式算法，它是借用能源有限元法的求解思维，驳回能源有限元中的Newmark方法来求解能源学识题。

在有限元方法中，对能源学识题的模拟驳回的方程为结构能源学方程（张国新等，2002；张子明等，2001）：

岩石断裂与挫伤

其中：M为品质矩阵；C为阻尼矩阵；d为位移增量；，区分是位移速度和减速度。K=K+K+K+K，K是刚度矩阵，K、K区分为块体及不延续面之间的接触矩阵，K是解放矩阵；F为总载荷向量，F=F+F+F-F+F+F+F，F是外载荷向量，F是体积力向量，F是已知解放位移惹起的等效载荷向量，F是初应力向量，F、F区分为法向和切向接触惹起的等效载荷向量，F为接触面之间的摩擦力惹起的等效载荷向量。

对普通的结构能源学识题，能源方程式（12-26）通常驳回Newmark解法。Newmark积分法实质上是线性减速度法的一种推行，它驳回如下假定（王勖成等，1995）：

岩石断裂与挫伤

其中：α、β是按积分精度和稳固性要求而选择的两个参数。当β≥0.5，α≥0.25（0.5+β）时，Newmark法是一种无条件稳固算法，而且求解时可以驳回比普通差分法大得多的期间步长，它比拟适宜于求解结构能源学识题。在这种方法中t+Δt时辰的位移解答d是经过满足t+Δt时辰的静止方程：而失掉的。

应用Newmark法求解静止方程的算法步骤如下：

1.初始计算

（1）构成刚度矩阵K，品质矩阵M和阻尼矩阵C。

（2）给定d、和。

（3）选用期间步长Δt、参数α和β，并计算积分常数：

岩石断裂与挫伤

（4）构成有效刚度矩阵：

岩石断裂与挫伤

（5）三角合成：

岩石断裂与挫伤

2.关于每一期间步长

（1）计算期间t+Δt的有效载荷：

岩石断裂与挫伤

（2）求解期间t+Δt的位移：

岩石断裂与挫伤

（3）计算期间t+Δt的速度和减速度：

岩石断裂与挫伤

综上可见，入选用适宜的参数后，Newmark法是一种无条件稳固的隐式算法，即期间步长Δt的大小不影响解的稳固性，此时Δt的选用关键依据解的精度要求而选用。

与数值流形方法中最后的能源学求解方法相比，这种算法是以对求逆失掉较大的期间步长（彭自强，2003）。

Newmark算法思考了能源求解疑问中的阻尼效应。

本章中能源学识题求解驳回Newmark算法。

二、灵活裂纹计算模型

与静态加载时一样，当灵活载荷作用于含裂纹的物体时，在裂纹尖端也将发生应力集中现象，这种裂纹尖端的应力集中现象用灵活应力强度因子来形容，而灵活应力强度因子岂但是载荷、裂纹长度和带裂纹体几何尺寸的函数，而且是期间的函数，有文献提出了灵活应力强度因子K（t）与相应的静态应力强度因子K（0）的相关式：

岩石断裂与挫伤

式中k（V）是裂纹裁减速度的函数，而与几何状态有关。其中k（V）可近似表白为

岩石断裂与挫伤

式中：C为Rayleigh外表波速；h是一个弹性波速的函数，有

岩石断裂与挫伤

其中：c、c区分为资料的收缩波速和剪切波速，有

岩石断裂与挫伤

其中：ρ、G、ν区分为资料的密度、剪切模量和泊松比。

依据以上的探讨，可以由资料的性质求出裂纹的灵活断裂应力强度因子，并依据相应的准绳来对灵活裂纹的裁减启动判别。灵活裂纹裁减的判据很多，最便捷适用的是应力强度因子判据，对Ⅰ型裂纹：K（t）≥K（V），K（V）是纯Ⅰ型裂纹的灵活断裂韧性，它是一个和裂纹裁减速度有关的函数，但其变动不是很大，因此在实践运行中，通常以为它是一个不变的资料常数（刘再华等，1996）。

关于平面Ⅰ、Ⅱ型混合灵活裂纹疑问来说，裂纹尖端的应力场和位移场的方式与式（12-3）齐全一样，所不同的就是把式（12-3）中的静态应力强度因子换为灵活应力强度因子即可（刘再华等，1996）。

虽然方式一样，但是由于灵活应力强度因子是期间的函数，所以灵活疑问的应力场和位移场也都是期间的函数，因此与静态疑问有着实质的不同。

上方就依据以上灵活裂纹裁减的计算原理对灵活裂纹裁减疑问启动模拟。

三、算例剖析

冲击载荷是工程通常中经常出现的一种灵活载荷，也是惹起大少数工程资料破坏的关键外载方式。

当冲击载荷作用于含裂纹等毛病的结构物时，这些结构物的灵活断裂行为与静态状况下有很大的差异。

冲击载荷作用于物体会发生应力波在物体中流传，当应力波作用于裂纹时，裂纹尖端的应力强度因子将因资料结构和裂纹形式的变动而具备不同的灵活照应，并具备不同的起裂和裁减行为。

同时，裂纹的灵活行为也对应力波的流传具备不同的散射作用，因此存在着各种波与裂纹间的相互作用相关。

所以对含裂纹毛病的结构在冲击载荷作用下灵活行为的钻研，不只有钻研资料的灵活力学功能，而且要对工程结构的安保评价等启动钻研。

上方应用数值流形方法对冲击载荷作用下裂纹的灵活裁减行为启动模拟。

为了与相关实验结果启动对比剖析，驳回姚学锋等（1996）中的实验模型启动模拟。计算模型为如图12-7所示的含偏置裂纹的简支梁，为了能够更好地观察到裂纹的曲裂现象，简支梁的宽度取得较大，详细尺寸如图所示，所受的载荷为一阶跃的冲击载荷，其大小为10N。梁的弹性常数为：弹性模量E=100GPa，泊松比ν=0.2。关于不同的偏离距离a启动模拟，模拟结果如图12-8所示。

图12-7 含偏置裂纹的简支梁（单位：cm）

由图12-8的模拟结果可以看出，在不同偏置距离的状况下，裂纹的裁减门路表现出不同的曲裂现象，其现象和要素剖析如下：

（1）裂纹开局起裂时裁减方向沿着初始裂纹方向，裁减到达必定长度后，开局偏离原来的方向出现曲裂，其方向是朝着冲击载荷作用点的方向。

（2）初始裂纹的偏置距离越大，曲裂出现的期间越早，沿初始裂纹方向裁减的长度也越小，曲裂现象越重大。

（3）裂纹向上裁减到必定长度后，又突然偏离了原来的方向而出现第二次曲裂，最终与上边界贯串。

图12-8 不同偏置距离时的裂纹裁减门路

发生曲裂现象的关键要素是：当冲击载荷发生的笔挺波作用于裂纹时，裂纹尖端作用的力是Ⅰ型和Ⅱ型复合载荷，裂纹的裁减属于Ⅰ-Ⅱ型复合裂纹。

偏置距离越小，Ⅰ型应力重量所占有的比重越大，所以刚开局起裂时，裂纹沿着原来的方向裁减。

随着偏置距离的参与，裂纹的曲裂水平也有所参与。

当裂纹裁减到梁的上半部时，裂纹在惯性作用下继续作曲裂静止，但遭到压应力的作用，使得裂纹的裁减速度有所减慢。

由于裂纹的裁减造成梁内的应力从新散布，在裂纹尖端的下部出现了一个新的核心点，并在该点的上部和下局部别构成了相应的拉应力重量和压应力重量，这些重量构成了新的弯矩使裂纹发生了第二次曲裂，并最终造成裂纹抵达上边界而发生结构的全体破坏。

图12-9为不同偏置距离状况下的裂纹尖端两类应力强度因子随期间的变动相关。从图中可以看到：

（1）K均大于零，而K有正有负。正K和负K的意义表现裂纹外表相对位移的方向雷同。

（2）应力强度因子在开局的时刻比拟小，随着裂纹裁减长度的参与，裂纹尖端的应力强度因子也迅速参与，最后造成裂纹贯串，在行将贯串的最后时辰，应力强度因子出现了动摇。但在整个环节中，K一直是比拟小的，这说明在曲裂环节中Ⅰ型裂纹裁减占主导作用，而Ⅱ型裂纹仅占无所谓位置。

图12-9 不同偏置距离时的应力强度因子随期间的变动曲线图（期间单位：s；应力强度因子单位：10N/m）

（3）从裂纹裁减期间过去看，当偏置距离较小时，裂纹裁减到梁顶部即贯串的期间较短，而随着偏置距离的参与，裂纹裁减的门路也随着延长，裂纹的贯串期间也增大。

（4）从裂纹裁减环节中应力强度因子的大小来看，在同一时辰，随着裂纹偏置距离的参与，裂纹尖端的应力强度因子是逐渐减小的，与通常剖析结果分歧，由于随着偏置距离的参与，在雷同的外力作用下，裂纹受力减小，因此裂纹尖端的应力强度因子也相应地有所减小。

各类钢材断裂要素汇总

钢材种类单一，每种钢材功能、化学成分或合金种类和含量的差异，选择了其商品称号。

但是，断裂韧性值在钢材选用时极为关键，却难以适用于一切钢材。

断裂要素关键分为两大类：一是钢材冶炼环节中的合金元素参与与热解决造成的显微组织变动，扭转了原有的功能；二是炼钢和浇注环节中发生的毛病，特意是集中毛病（如气孔、夹杂等）在轧制时尤为敏感。

在同一化学成分的钢中，这些毛病在不同炉次和钢坯部位的变动，影响着钢材品质。

钢材韧性关键取决于显微结构和毛病的分散水平，而非化学成分。

热解决后韧性会清楚变动。

深化摸索钢材功能及其断裂要素，需把握物理冶金学与显微组织与钢材韧性的相关。

铁素体-珠光体钢，约占钢总产量的绝大少数，通常含碳量在0.05%～0.20%。

它们经过参放大批合金元素以提高屈服强度和韧性。

铁素体-珠光体的显微组织由铁素体、珠光体、可溶合金和Fe3C组成。

在低碳钢中，渗碳体颗粒在铁素体晶粒边界和晶粒中，而当碳含量高于0.02%时，大局部Fe3C构成片状结构称为珠光体。

含碳量在0.10%～0.20%时，珠光体含量占10%～25%。

珠光体使钢材韧性参与，而其含量与散布影响钢材功能。

解决工艺影响钢材的冲击功能。

水淬火钢的冲击功能优于退火或正火钢，快冷阻止了渗碳体在晶界构成，促使铁素体晶粒变细。

热轧条件对冲击功能有关键影响，较低的终轧温度会降落冲击转变温度，增大冷却速度，促使铁素体晶粒变细，提高钢材韧性。

厚板因冷却速度较慢，铁素体晶粒比薄板粗大，因此韧性较差。

正火解决可以改善钢板功能。

热轧后，钢材或者发生各向同性，构成珠光体带、夹杂晶界与轧制方向分歧的定向韧性钢。

铁素体-可溶合金元素的参与，旨在消费在某些环境温度下的固溶体软化钢，提高晶格摩擦应力。

锰的参与可降落晶粒边界渗碳体薄膜构成的趋向，减小铁素体晶粒尺寸，发生少量而粗大的珠光体颗粒。

在低碳钢中，锰含量的参与降落约50%的转变温度，但高锰含量对正火钢的冲击功能影响清楚。

镍可改善铁-碳合金韧性，参与量到达2%时，防止热轧态和正火钢晶界渗碳体的构成，实质降落开局转变温度，升高夏比冲击曲线峰值。

磷是铁素体的稳固剂，参与钢中参与δi值和铁素体晶粒尺寸，提高韧性。

硅作为脱氧剂，同时有益于提高冲击功能，参与量为0.1%，可使50%转变温度升高约44℃。

铝与氮联合生成AlN，细化了铁素体晶粒，降落转变温度约40℃。

氧在晶界发生偏析造成铁合金晶连续裂，经过参与碳、锰、硅、铝和锆去除游离氮，减轻氧脆化疑问。

亚共析钢的含碳量在0.3%～0.8%，先共析铁素体是延续相，珠光体在奥氏体晶粒内构成，占显微组织的35%～100%。

珠光体强度比先共析铁素体高，限度了铁素体的流动，提高钢的屈服强度和应变软化率。

珠光体含量参与，限度造用增强，形变环节中会在低平和/或高应变率时构成微型解理裂纹。

在铁素体片之间、相邻汇集组织中的铁素体晶粒内有择优取向。

贝氏体钢在含碳量为0.10%的低碳钢中参与0.05%钼和硼，可优化奥氏体-铁素体转变。

上贝氏体和下贝氏体构成，由针状铁素体和分散的碳化物组成。

转变温度的降落细化了贝氏体铁素体片，提高了抗拉强度。

下贝氏体针状铁素体晶粒尺寸仅为上贝氏体奥氏体晶粒尺寸的1/2，断裂特色取决于针状铁素体晶粒尺寸，回火下贝氏体的夏比冲击功能优于未回火上贝氏体。

马氏体钢经过提前奥氏体转变成铁素体或珠光体，极速冷却使奥氏体转变成马氏体，构成现实的断裂特色。

马氏体具备粗大的四面体铁素体或针状马氏体，碳原子在剪切环节中来不迭分散，造成马氏体晶粒拉长，晶格收缩。

马氏体转变后的结构为马氏体和剩余奥氏体的混合结构。

回火是稳固钢功能的关键步骤，经过转变剩余奥氏体和参与弥散性碳化物含量，降落屈服强度，同时提高冲击韧性。

中强度钢的回火包含转变剩余奥氏体和参与弥散性碳化物含量，同时降落屈服强度。

在特定温度范围内回火，可防止构成脆性珠光体，提高抗拉塑性。

回火脆性是可逆的，经过再次加热至临界范围，可以恢复资料的功能。

高强度钢经过淬火和回火、淬火和回火前奥氏体变形、退火和时效消费积淀软化钢，以及应变和再回火或回火期应变，进一步提高钢的强度。

不锈钢关键由铁-铬、铁-铬-镍合金和其余改善力学功能与抗蚀才干的元素组成。

不锈钢的防蚀性依赖于在金属外表生成的可防止进一步氧化的铬氧化物层。

不锈钢在氧化气氛中防蚀，但在恢复气氛中铬氧化层遭到侵害。

铬和镍含量的参与提高了抗蚀性，镍片面优化铁的钝化性。

参与碳以改善力学功能和保障奥氏体不锈钢功能的稳固性。

不锈钢依据显微组织分类为马氏体不锈钢、铁素体不锈钢和奥氏体不锈钢。

马氏体不锈钢属于铁-铬合金，经过奥氏体化和后续热解决生成马氏体。

铁素体不锈钢含铬约14%～18%，镍含量为8%，铬含量为18%（300型）能使奥氏体相十分稳固，奥氏体不锈钢相似于铁素体不锈钢，不能经过马氏体转变而软化。

奥氏体不锈钢在冷冻温度下都无法能解理断裂，适用于导弹系统和液氧贮存箱，但不稳固，或者出现塑性变形造成脆性断裂。

稳固化解决经过参与钛或铌与碳构成合金碳化物，防止铬耗尽和应力侵蚀裂纹。

奥氏体不锈钢在高温下用于压力容器，防止侵蚀和抗蠕变。

在焊接和高温环境下，铌或钛碳化物会在晶粒内和晶界积淀，造成裂纹发生，影响经常使用寿命，必定给予高度注重。