裂纹尖端应力场,应力强度因子
Ⅰ型裂纹尖端圆弧对应力强度因子影响的数值研究
第 ４期
２１ ００年 ２月
科
学
技
术
与
工
程
Ｖｏ ． ０ Ｎｏ ４ Ｆ ｂ． ０１ １１ ． ｅ ２ ０ ⑥ ２ ０ Ｓ ｉ Ｔ ｃ Ｅ ｇ ｇ ０１ ｃ ． ｅ ｈ． ｎ ｎ ．
１ ７ —１ １ （ ０ ０ ４ ０ ６ －５ ６ １ ８５ ２ １ ） － １ ９ ０
Ｓｉｎｅ ｅｈｏｏｙａｄＥ ￣ｎｅ ｎ ｃ ｃ ｃｎｌｇ ｎ ｎ ｅｆ ｇ ｅ Ｔ ｉ
算 和工 程 力学 。
９２ ６
科
学
技
术
与
工
程
１ Ｏ卷
０７ ， ． ％ 分析结果 准确 可信 。该 曲线 为带 圆弧 裂纹 的 应力 强度 因子测试 和带 圆弧裂 纹构 件强度 计算 提供
了依据 。
ｔｓ导ｉ －ｎｓ ｉ ｎ
Ｋ ／ ￣ １＝ ， ｆ ２ １
￣ ０ ． ３０ ・ ５
摘
要
应 力 强度 因子 表 征 了裂 纹 尖 端 奇 异 应 力 场 的 强度 ， 它是 研 究裂 纹 扩 展 规 律 和 带 裂 纹 构 件 强度 的 基 础 。 采 用 有 限 元
法 ， 受 均布 荷 载 作用 存 在 边 缘 Ｉ 裂 纹 的平 面 板 进行 了数 值 分 析 。 研 究 了裂 纹 尖端 圆弧 对 应 力 强度 因 子 的 影 响 ， 别 计 算 对 型 分 了具 有不 同裂 尖 圆 弧 ｒ Ｉ ｎ的 型裂 纹 的 应 力 强度 因子 。 采 用 应 力 法计 算 不 同 半 径 处 的表 观 应 力 强度 因子 ， 值 到 裂 尖 圆 弧 而 插 得 。 根据 计 算 结 果 ， 制 Ｋ ｒ曲线 ， 用 最 小 二 乘 法拟 合 至 尖 裂纹 （ｏ ０ 即得 理 想 尖裂 纹 的应 力 强度 因 子 ， 解 析 解 相 差仅 绘 ． ｎ 利 ｒ＝ ） 与 ０ ７ 。该 曲线 为 带 圆弧 裂 纹 的应 力 强 度 因 子测 试 和 带 圆 弧裂 纹 构 件 强度 计 算 提 供 了依 据 。 ．％ 关键 词 裂纹尖端圆弧 应 力 强度 因子 有 限元 法
matlab裂纹扩展计算公式
MATLAB是一种用于算法开发、数据分析、数据可视化和数值计算的高级技术计算语言和交互式环境。
在工程和科学领域中，MATLAB经常被用来进行裂纹扩展计算。
裂纹扩展计算是研究材料断裂行为的重要课题，也是工程设计和材料研发中不可或缺的一环。
在MATLAB中进行裂纹扩展计算，需要使用一些特定的公式和算法。
下面将介绍一些常用的裂纹扩展计算公式及在MATLAB中的实现方法。
1. 裂纹长度计算公式裂纹长度是裂纹扩展计算中的重要参数，通常使用Griffith裂纹力学理论进行计算。
Griffith裂纹力学理论认为，裂纹的扩展速度与应力强度因子成正比。
根据Griffith裂纹力学理论，裂纹长度计算公式如下：\[ a = \left(\frac{2\gamma E^*}{\pi\sigma^2}\right)^{\frac{1}{2}} \]其中，a为裂纹长度，γ为材料的表面能，E*为有效断裂韧度，σ为应力。
在MATLAB中，可以使用如下代码实现裂纹长度的计算：```matlabfunction a = crackLength( gamma, E_star, sigma )a = sqrt(2 * gamma * E_star / (pi * sigma^2));end```2. 裂纹扩展速率计算公式裂纹扩展速率是裂纹扩展过程中的另一个关键参数。
根据线弹性断裂力学理论，裂纹扩展速率与应力强度因子的变化率成正比。
裂纹扩展速率计算公式如下：\[ \frac{da}{dt} = C(\Delta K)^n \]其中，\(\frac{da}{dt}\)为裂纹扩展速率，C为材料常数，\(\Delta K\)为应力强度因子的变化量，n为指数。
在MATLAB中，可以使用如下代码实现裂纹扩展速率的计算：```matlabfunction da_dt = crackGrowthRate( C, delta_K, n )da_dt = C * delta_K^n;end```3. 应力强度因子计算公式应力强度因子是裂纹扩展计算中的关键参数，描述了裂纹尖端应力场的分布。
应力场强度因子
应力场强度因子
应力场强度因子是研究材料断裂行为的重要参数之一。
它是描述材料在受到外力作用下，裂纹尖端应力场的强度和分布情况的物理量。
应力场强度因子的大小和方向对材料的断裂行为有着重要的影响。
应力场强度因子的计算方法有多种，其中最常用的是Williams和Landel的方法。
该方法基于线弹性力学理论，通过对裂纹尖端应力场的分析，得出了应力场强度因子的计算公式。
该公式中包含了裂纹尖端应力场的强度和分布情况，因此可以用来预测材料的断裂行为。
应力场强度因子的大小和方向对材料的断裂行为有着重要的影响。
当应力场强度因子达到一定的临界值时，裂纹尖端的应力场会达到材料的断裂强度，从而导致材料的断裂。
因此，应力场强度因子可以用来预测材料的断裂强度和断裂模式。
除了预测材料的断裂行为外，应力场强度因子还可以用来优化材料的设计和制备。
通过对应力场强度因子的分析，可以确定材料的最大承载能力和断裂模式，从而优化材料的设计和制备过程。
应力场强度因子是研究材料断裂行为的重要参数之一。
它可以用来预测材料的断裂行为、优化材料的设计和制备过程，对于提高材料的性能和可靠性具有重要的意义。
裂纹扩展k判据名词解释
裂纹扩展k判据名词解释
裂纹扩展K判据是用于研究裂纹在材料中扩展的一种方法。
K判据中的K值代表裂纹尖端的应力强度因子，它是衡量裂纹扩展能力的一个重要参数。
K值的大小取决于裂纹尖端的应力场和材料的弹性性质。
裂纹扩展K判据主要分为三类：基于线性弹性理论的判据、基于弹塑性理论的判据和基于断裂力学理论的判据。
其中，基于线性弹性理论的判据主要适用于裂纹扩展速度较慢、材料较脆性的情况下；基于弹塑性理论的判据适用于裂纹扩展速度较快、材料较韧性的情况下；基于断裂力学理论的判据则适用于复杂的载荷条件和裂纹扩展的非线性情况。
裂纹扩展K判据的应用十分广泛，可以用于分析材料的疲劳寿命、评估结构的可靠性、预测裂纹扩展的趋势等。
准确地计算裂纹扩展K值和选择合适的K判据对于材料的性能分析和结构设计具有重要意义。
T型焊接件的表面裂纹应力强度因子的计算
中图分类号：U661
文献标识码：A
文章编号：1006-7973（2021）02-0057-02
引言 T 型焊接接头是工程结构中常采用的典型结构，其疲劳 断裂性能严重影响着焊接结构的使用安全[1]。应力强度因子 是表征裂纹尖端应力场的重要参量[2]，研究含裂纹结构强度 时必须计算裂纹处应力强度因子大小，以判断结构工作时的 安全状态。随着有限元理论的发展和有限元软件功能的完善， 现广泛借助有限元软件求解应力强度因子[3]。 在实际工程结构中，由于疲劳、缺陷、应力集中等会出 现各种形式的裂纹，表面裂纹是其中常见的裂纹形式。本文 以 T 型焊接件为研究对象进行焊缝表面裂纹的热-结构耦合 应力强度因子的计算，利用稳态热分析模块和静力学分析模 块对焊接件进行稳态热应力分析，可为同类型应力强度因子 的计算提供借鉴，对工程的实际应用有重要的应用价值。 二、T 型焊接件有限元模型建立 1．T 型焊接件结构分析 本文中 T 型焊接件为结构钢直接焊接而成，T 型焊接件 所用材料为工程常用材料 Q235A 钢，其材料参数见表 1 所 示，结构模型如图 1 所示。
表 1 Q235A 钢材料参数
参数
值
参数
值
弹性模量/MPa
2.1
屈服强度/MPa
≥235
密度/（kg/m³）
7,900
泊松比
0.277
2．表面裂纹的插入 将用 Solidworks 软件建立的三维立体模型导入 ANSYS 软件中，工程上 3D 表面裂纹的形状一般采用半椭圆来描述。 在一侧的焊缝表面建立局部坐标系，此局部坐标系用来给裂 纹定位，插入表面半椭圆裂纹（工程中最常见的裂纹类型）， 裂纹的长半轴（裂纹长度）为 0.008m，裂纹的短半轴（裂 纹深度）为 0.003m，影响区域为 0.008m，裂纹前缘的径 向划分为 20 份，周向划分为 16 份（需为 8 的倍数），最大 绕线半径为 0.0035m。 3．网格划分 利用 Ansys 分析求解获得足够精度的应力强度因子的 关键是单元的选取和网格的划分。进行网格划分，由于焊接 件结构简单，直接用 workbench 进行四面体网格主导方法 进行划分。 利用有限元法求解应力强度因子，裂纹前沿周围网格的 划分是关键[3]。为了兼顾模拟结果的精度，裂纹所在面需要 更精密的网格划分，基体的网格划分太粗糙，所以用局部控 制 Sizing 设置该焊接面网格尺寸 0.001m，重新细化该面网 格，提高计算精度，T 型焊接件模型网格划分结果如图 2 所 示。
(整理)应力强度因子的计算
第二章 应力强度因子的计算K --应力、位移场的度量⇒K 的计算很重要,计算K 值的几种方法: 1.数学分析法:复变函数法、积分变换； 2.近似计算法：边界配置法、有限元法； 3.实验标定法：柔度标定法； 4.实验应力分析法：光弹性法.§2-1 三种基本裂纹应力强度因子的计算一、无限大板Ⅰ型裂纹应力强度因子的计算K Z ⅠⅠx ®=计算K 的基本公式，适用于Ⅱ、Ⅲ型裂纹.1.在“无限大”平板中具有长度为2a 的穿透板厚的裂纹表面上，距离x b =±处各作用一对集中力p .Re Im x Z y Z σ'=-ⅠⅠRe Im y Z y Z σ'=+ⅠⅠRe xy y Z τ'=-Ⅰ选取复变解析函数：222()Z z b π=- 边界条件：a.,0x y xy z σστ→∞===.b.,z a <出去z b =±处裂纹为自由表面上0,0y xy στ==。
c.如切出xy 坐标系内的第一象限的薄平板，在x 轴所在截面上内力总和为p 。
y '以新坐标表示：Z=⇒lim()K Zξξ→==Ⅰ2.在无限大平板中,具有长度为2a的穿透板厚的裂纹表面上，在距离1x a=±的范围内受均布载荷q作用.利用叠加原理:微段→集中力qdx→dK=Ⅰ⇒K=⎰Ⅰ令cos cosx a aθθ==,cosdx a dθθ=⇒111sin()1cos22()cosaa aaaK daθθθ--==Ⅰ当整个表面受均布载荷时,1a a→.⇒12(aaK-==Ⅰ3.受二向均布拉力作用的无限大平板,在x轴上有一系列长度为2a,间距为2b 的裂纹.边界条件是周期的: a. ,y x z σσσ→∞==.b.在所有裂纹内部应力为零.0,,22y a x a a b x a b =-<<-±<<±在区间内0,0y x y στ==c.所有裂纹前端y σσ> 单个裂纹时Z =又Z 应为2b 的周期函数⇒sinzZ πσ=采用新坐标:z a ξ=-⇒sin()a Z πσξ+=当0ξ→时,sin,cos 1222bbbπππξξξ==⇒sin()sin cos cos sin22222a a a bbbbbπππππξξξ+=+σcossin222a a bbbπππξ=+2222[sin()]()cos 2cos sin(sin)2222222a a a a a bbbbbb bπππππππξξξ+=++22[sin()](sin )2cos sin22222a a a a bbbbbπππππξξ⇒+-=sinaZ ξπσ→⇒=sinlim aK ξπσ→⇒===Ⅰ=取w M =修正系数,大于1,表示其他裂纹存在对K Ⅰ的影响. 若裂纹间距离比裂纹本身尺寸大很多(2125a b ≤)可不考虑相互作用,按单个裂纹计算.二、无限大平板Ⅱ、Ⅲ型裂纹问题应力强度因子的计算 1.Ⅱ型裂纹应力强度因子的普遍表达形式(无限大板):lim (K Z ξξ→=Ⅱ2.无限大平板中的周期性的裂纹,且在无限远的边界上处于平板面内的纯剪切力作用.τsin()zZ z πτ=sin()()a Z πτξξ+=lim ()K ξξ→⇒==Ⅱ3.Ⅲ型裂纹应力强度因子的普遍表达形式(无限大板):lim ()K ξξ→=Ⅲ4.周期性裂纹:K =§2-2 深埋裂纹的应力强度因子的计算1950年，格林和斯内登分析了弹性物体的深埋的椭圆形裂纹邻域内的应力和应变，得到椭圆表面上任意点,沿y 方向的张开位移为:1222022(1)x z y y a c=--其中:202(1)ay E μσ-=Γ.Γ为第二类椭圆积分.有φϕ= (于仁东书) 1222220[sin ()cos ]a d cπϕϕϕ=+⎰(王铎书)1962年,Irwin 利用上述结果计算在这种情况下的应力强度因子σ原裂纹面11cos ,sin z x ρϕρϕ==又222222221111221x z c x a z a c a c+=⇒+= ⇒ρ=假设:椭圆形裂纹扩展时,其失径ρ的增值r 与ρ成正比.r f ρ= (f 远小于1)r f ρ⇒==边缘上任一点(,)p x z ''',有:1()sin (1)sin (1)x r f f x ρϕρϕ'=+=+=+1()cos (1)z r f z ρϕ'=+=+11(,),(,)p x z p x z '''⇒均在0y =的平面内. 222242222(1)c x a z f a c a c ''''''⇒+=+=⇒新的裂纹面仍为椭圆.长轴(1)c f c '=+,短轴(1)a f a '=+. ⇒y 向位移22002(1)2(1)(1)(1)a f a y f y E E μσμσϕϕ'--+'===+原有裂纹面:222220()1x z ya c y ++=扩展后裂纹面:222220()1x z y a c y '''++='''以1x x '=,1z z '=,代入⇒原有裂纹面的边缘y 向位移y ',有2222211112222222011(1)(1)x z x z y y a c f a f c'=-+=--'''++2222221111112222221(12)(12)12()x z x z x z f f f a c a c a c----=--++2f =2222200022(1)2y fy f f y fy ''⇒==+又f =⇒2y '=设各边缘的法向平面为平面应变,有:31)sin sin ]22v k θθ=+- 其中34k μ=-当θπ=时24(1)v K E μ-=222216(1)2I r K E μπ-⇒=22021E ()41I K y acπμ⇒=-又202(1)ay E μσϕ-=14122222()(sin cos )I a K c a cϕϕφ⇒=+在椭圆的短轴方向上，即2πϕ=，有I ImaxK K φ== 危险部位 →椭圆片状深埋裂纹的应力强度因子当a c =时→圆片状裂纹，2πφ=2I K π⇒=§2-3 半椭圆表面裂纹的应力强度因子计算一、表面浅裂纹的应力强度因子当a B （板厚）→线裂纹⇒可以忽略后自由表面对A 点应力强度的影响 欧文假设：半椭圆片状表面线裂纹I K 与深埋椭圆裂纹的I K 之比等于边裂纹平板与中心裂纹平板的I K 值之比。
应力强度因子的数值计算方法
应力强度因子的数值计算方法应力强度因子是用来描述裂纹尖端应力场的重要参数，它在研究裂纹扩展、断裂行为等问题中具有重要的应用价值。
本文将介绍应力强度因子的数值计算方法，包括解析方法和数值方法。
一、解析方法解析方法是指通过求解弹性力学方程，得到应力场的解析表达式，进而计算应力强度因子。
常见的解析方法有：1. 爱尔兰函数法：该方法适用于轴对称问题，通过引入爱尔兰函数，将弹性力学方程转化为常微分方程，进而得到应力强度因子的解析表达式。
2. 奇异积分法：该方法适用于不规则裂纹形状或复杂载荷情况。
通过奇异积分的性质，将应力场分解为奇异和非奇异两部分，进而得到应力强度因子的解析表达式。
3. 线性弹性断裂力学方法：该方法通过建立合适的应力强度因子与裂纹尺寸之间的关系，利用裂纹尖端应力场的奇异性，通过分析弹性力学方程的边界条件，得到应力强度因子的解析表达式。
二、数值方法数值方法是指通过数值计算的方式，求解弹性力学方程，得到应力场的数值解，从而计算应力强度因子。
常见的数值方法有：1. 有限元法：有限元法是一种广泛应用的数值方法，通过将结构离散为有限个单元，建立节点间的关系，利用数值方法求解离散方程组，得到应力场的数值解，进而计算应力强度因子。
2. 边界元法：边界元法是一种基于边界积分方程的数值方法，通过将边界上的应力场表示为边界积分方程的形式，利用数值方法对积分方程进行离散求解，得到应力场的数值解，进而计算应力强度因子。
3. 区域积分法：区域积分法是一种基于区域积分方程的数值方法，通过将应力场表示为积分方程的形式，利用数值方法对积分方程进行离散求解，得到应力场的数值解，进而计算应力强度因子。
以上介绍了应力强度因子的数值计算方法，包括解析方法和数值方法。
解析方法适用于问题简单、载荷条件规则的情况，可以得到解析表达式并具有较高的精度；数值方法适用于问题复杂、载荷条件不规则的情况，通过数值计算可以得到应力场的数值解，并利用数值解计算应力强度因子。
断裂力学3裂纹尖端应力场和位移场计算
25
线弹性裂尖场特点
❖ 三种情况下的K场有相似的形式，分别由应力强度 因子决定着其场的强度。SIF取决于外加载荷，而 且与构件几何、裂纹尺寸有关，但是与( )坐标 无关。在K场范围内，应力和应变均正比于SIF，所 以SIF是裂纹尖端附近应力、应变场强度的表征， 是描述裂尖场强度的参数。
应力强度因子
——名义应力，即裂纹位置上按无裂纹计算的应力 ——裂纹尺寸，即裂纹长或深
——形状系数，与裂纹大小、位置有关 应力强度因子单位：N.m-3/2
28
应力强度因子
鉴于应力强度因子的重要性，在断裂力学这门科学近半个 世纪的快速发展中，应力强度因子的分析计算一直是一个 经久不衰的研究课题，这可从这方面的专著（如二十世纪 七十年代Sih的专著和近期的专著）和专门的应力强度因子 手册可见一斑。从研究方法上，从解析的Westergaard stress function、 Muskhelishvili stress function 到解 析的或半解析的Green Function、Singular Integral Equation、Conforming Mapping（保形映射）, 及数值 方法如Boundary Collocation Method， Finite Element Method (有限元法)和Boundary Element Method (边界元法)。
23
通过前面的推导，各种类型裂尖应力和来自移场可表示为若上标写成II、III，代表II型或III型裂纹。 裂纹尖端应力场是渐进解，仅仅适合于裂纹尖端附近
24
线弹性裂尖场特点
❖ 三种变形情况下裂纹尖端应力场和应变场都具有奇异 性，即在裂纹尖端处，应力和应变为无穷大，这种不 真实的性质是由于所采用的本构关系所决定的，即认 为材料能承受无限大的应力，且应变与应力呈线性关 系。另外，在上述的分析中，裂纹假设成理想的尖裂 纹，即裂纹尖端曲率为无穷大。实际上，裂纹尖端不 可避免地会出现塑性区，并且裂纹尖端地曲率是有限 的，但是在塑性区很小的情况下，在围绕裂尖的一个 环状区域内K场是适用的。
ABAQUS计算裂纹尖端应力强度因子有效性的算例研究
ABAQUS计算裂纹尖端应力强度因子有效性的算例研究摘要：在实际工程领域中，相当部分的脆性材料总是不可避免的存在着裂纹或是缺陷。
在实际环境中材料的受力往往是相当复杂的。
基于ABAQUS平台的裂纹仿真软件，它具有简单易用的特点。
通过算例分析验证表明，该软件的计算结果具有较高的精度，完全可以用于实际工程问题的计算，通过分析验证表明该软件的设计是成功的。
此外，今后可以在它的基础上进行更多功能扩展，从而使它拥有分析更为复杂问题的能力。
关键词：裂纹；应力强度因子；断裂力学；ABAQUS引言材料在成型和加工过程中在其内部造成了很多缺陷，而其破坏正好均源于构件内部的微小裂纹，所以研究带裂纹的物体力学性能具有十分重要的意义。
图1存在于岩石和混凝土地面中的裂缝1920年， Griffith[1-2]提出了在材料中存在裂纹的设想，而从Irwin[]3-4]在1957年提出了应力强度因子以及其后形成的断裂韧度的概念后，断裂力学理论出现了重大的突破，奠定了线弹性断裂力学的基础。
1基本原理近年来以数值分析为基础的手段来解决断裂力学相关问题的技术得到了广泛的发展应用，并且不断的调整完善。
该技术在一定程度上较好的克服了实验条件下的不足。
对于线弹性断裂力学而言，裂尖区域的位移场、应力、应变场由应力强度因子决定，故而通过有限元计算的结果来得到具体的应力强度因子的值是线弹性断裂力学中用有限元法的基本要求。
1.1 ABAQUS求解裂纹尖端的应力强度因子传统的有限元在计算裂纹尖端的应力强度因子的时候，无可避免地遇到裂尖复杂应力场和位移场的计算，J积分则可以完全避免这种复杂的处理过程。
为了计算二维情况下的J积分，ABAQUS定义了围绕裂纹尖端由单元组成的环形的积分域，如下图所示。
图2 ABAQUS中围线的定义ABAQUS在计算围线积分时，采用的是先计算出围线上面所取的若干个离散点处J积分值，然后乘以每个点对应的加权值后，所有点相加来近似地求解出围线积分，即J积分的值和，进而得到复合裂纹的应力强度因子和。
裂纹尖端应力场,应力强度因子
GI
=
s2a E
GI
=
s2 (1 2 )a E
GI
=
K
2 I
E
GI
=
K2I (1 2 ) E
Sliding mode analysis or Mode 2
For problems with crack tip under shear loading, Airy’s stress
function is taken as
f
ij
and
KI Ys
a
where Y = configuration factor
From Hooke’s law, displacement field can be obtained as
u
2(1 E
)
KI
r 2
cos
2
2
1
sin2
2
v
2(1 E
)
KI
r 2
sin
2
compatibility Conditions
0 4
Then
sx
2 y2
,
sy
2 x2
,
txy
2 xy
For problems with crack tip Westergaard introduced Airy’s stress
function as
Re[Z] y Im[Z]
Where Z is an analytic complex function
sz = ν (sx+sy) where ν is Poisson's ratio.
For plane problems, the equilibrium conditions are:
对应力强度因子Ｋ物理含义的一点理解
对应⼒强度因⼦Ｋ物理含义的⼀点理解 线弹性断裂⼒学中的的应⼒强度因⼦Ｋ是⼀个容易让初学者不甚明了的量。
即使已经了解了Irwin理论的推导过程和Ｋ因⼦的得出⽅法，对于究竟Ｋ是个什么东东？它的物理意义究竟是什么还是难以说得清楚。
在此想谈谈个⼈的⼀点理解。
从何处谈起呢？还是从欧⽂的“裂纹长度ａ修正”谈起吧。
欧⽂为了将其的线弹性断裂⼒学理论，既Ｋ因⼦的理论推⼴到裂尖出现了微⼩塑性区的情形，提出了修正裂纹长度的补救⽅法。
认为，裂纹尖端出现塑性区会导致整个受⼒构件的柔度提⾼，降低了结构的承载能⼒，这种情况就等价于⼀个长度更长的裂纹存在于结构中。
因此，在原来的裂纹长度ａ的基础上将裂纹长度向更长的⽅向修正为ａ+ｒp，仍然使⽤线弹性⽅法来评估结构的断裂问题。
表⾯看来这个修正法不⽆道理，但我必须说这个修正法是修正的⽅向发⽣了错误了！裂纹长度ａ往较短的⽅向作出修正才是物理意义上合理的。
加长的修正法只能是为了使⽤欧⽂的线弹性理论的不得不为之的权宜之计。
事实说明，这种修正法在实际中确实并未明显增加其理论的适⽤范围，很快就只能让位于更合理些的弹塑性断裂⼒学理论了。
因此，表⾯看起来是往长或短的⽅向修正的⼆选⼀，其实带有物理意义上的原则性差别。
裂纹变长和裂尖出现塑性区虽然就宏观效果来看都是导致“结构的柔度提⾼”，但从对Ｋ值的实际影响却是刚好相反的。
较长裂纹会在外载荷不变的情况下增⼤Ｋ值，⽽形成塑性区则会导致裂尖钝化，结果却是降低应⼒集中程度，从⽽势必降低Ｋ值。
另外，从对促进裂纹开裂的动⼒⾓度看，形成塑性区会耗散掉⼀部分弹性变形能，因此是⼀个阻裂因素。
⽽裂纹加长则是⼀个促裂因素。
从这个⾓度⾔，出现塑性区后裂纹长度应该向缩短的⽅向修正才是合理的。
但这样⼀来，缩短裂纹，结构的柔度⼜应该降低，加之向缩短裂纹的⽅向修正意味着进⼊到了⾮线性弹性区甚⾄裂纹空腔区，这样线弹性⽅程成⽴的基础也不存在了。
由此可见欧⽂修正法遇到了⽆法两全的⽭盾！究竟应该如何修正才是合理的？究竟是应该以结构柔度变化为依据还是应该以促裂或阻裂为依据呢？要回答这个问题就必须要真正理解应⼒强度因⼦Ｋ究竟表达了什么样的物理含义这个极具关键性的问题。
金属学及热处理课后习题答案解析第六章
⾦属学及热处理课后习题答案解析第六章第六章⾦属及合⾦的塑性变形和断裂2）求出屈服载荷下的取向因⼦，作出取向因⼦和屈服应⼒的关系曲线，说明取向因⼦对屈服应⼒的影响。
答：1）需临界临界分切应⼒的计算公式：τk=σs cosφcosλ，σs为屈服强度=屈服载荷/截⾯积需要注意的是：在拉伸试验时，滑移⾯受⼤⼩相等，⽅向相反的⼀对轴向⼒的作⽤。
当载荷与法线夹⾓φ为钝⾓时，则按φ的补⾓做余弦计算。
2）c osφcosλ称作取向因⼦，由表中σs和cosφcosλ的数值可以看出，随着取向因⼦的增⼤，屈服应⼒逐渐减⼩。
cosφcosλ的最⼤值是φ、λ均为45度时，数值为0.5，此时σs为最⼩值，⾦属最易发⽣滑移，这种取向称为软取向。
当外⼒与滑移⾯平⾏（φ=90°）或垂直（λ=90°）时，cosφcosλ为0，则⽆论τk数值如何，σs均为⽆穷⼤，表⽰晶体在此情况下根本⽆法滑移，这种取向称为硬取向。
6-2 画出铜晶体的⼀个晶胞，在晶胞上指出：1）发⽣滑移的⼀个滑移⾯2）在这⼀晶⾯上发⽣滑移的⼀个⽅向3）滑移⾯上的原⼦密度与{001}等其他晶⾯相⽐有何差别4）沿滑移⽅向的原⼦间距与其他⽅向有何差别。
答：解答此题⾸先要知道铜在室温时的晶体结构是⾯⼼⽴⽅。
1）发⽣滑移的滑移⾯通常是晶体的密排⾯，也就是原⼦密度最⼤的晶⾯。
在⾯⼼⽴⽅晶格中的密排⾯是{111}晶⾯。
2）发⽣滑移的滑移⽅向通常是晶体的密排⽅向，也就是原⼦密度最⼤的晶向，在{111}晶⾯中的密排⽅向<110>晶向。
3）{111}晶⾯的原⼦密度为原⼦密度最⼤的晶⾯，其值为2.3/a2，{001}晶⾯的原⼦密度为1.5/a24）滑移⽅向通常是晶体的密排⽅向，也就是原⼦密度⾼于其他晶向，原⼦排列紧密，原⼦间距⼩于其他晶向，其值为1.414/a。
6-3 假定有⼀铜单晶体，其表⾯恰好平⾏于晶体的（001）晶⾯，若在[001]晶向施加应⼒，使该晶体在所有可能的滑移⾯上滑移，并在上述晶⾯上产⽣相应的滑移线，试预计在表⾯上可能看到的滑移线形貌。
基于ANSYS的裂纹尖端应力强度因子研究
’。 －，、 ＾ ， －。 、 ，、 ＾－ ， － 、 ，、 ＾ ， ， ’ ’，、 一 ＿
Ｇ＝ ４ＧＰａ ４ 。
裂 纹是低 应力 断裂 的主要 原 因 ，
采用 Ｐ Ａ Ｅ Ｌ Ｎ ２单元 划 分 网 格 。焊 缝 区域 网格 划 分 单
。 ， 。 是 最 危 险 的 ，也 是 多 年来 实 验 和 元 设 置 为 ２ 热影 响 区设 置 为 ４ 其 他 区域设 置 为 ６ 裂纹 。图 ４ 研 究 的 主体 ，当 实 际 裂 纹 是 复 合 区域网格尽量细分。划分网格后的整体模型如图 ３ 型裂 纹 时 ， 往 作 为 Ｉ 处 理 ， 往 型 这
载 荷 的作 用 下 也会 继 续 扩 展 。 型 Ｉ
种 材 料 属 性 ，分 别
对应 焊 接接 头 的每 个 区域 。本文 通 过 改变 屈 服强 度
搭 配 组 合 来 模 表 如
ｌ 。其 它 材 料 属 性 设 置 ： Ｅ＝１ ８ Ｐ ， ＝ ．４， ８ Ｇ ａ １ Ｇ ａ ／ ０ ３ Ｇ＝ ０ Ｐ ， ｘ
值， 代入 式 （ ） 外 推 得 到 ３并
Ｋ ｌ ２ ｕ／ Ｉ ｍ Ｇ 、兰 ＝ ｉ
Ｖ ｒ
（） ４
［ ］ 范引鹤 ， ３ 阮米庆 ． 焊接接头力 学不均 匀性的研究 ［ ］ Ｊ． 航空 学报 ，
１９ （ ）６ ８ ６ ０ ９ ６ ５ ：３ — ４ ．
３ 用 ＡＮＳ ＹＳ模 拟 焊接 接 头求 解应 力 强度 因子
４ 结果 与讨 论
）２ｉ２ ］ ＋ｓ ｎ０
二
（） １ （） ２
由于 焊 接接 头 每个 区域 的力 学 性 能定 义 不 同 ，计 算 出来 的应 力 强度 因子 值 就反 映 出材 料 力 学性 能 对应 力 强 度 因子 的 影 响规 律 。Ａ Ｓ Ｓ 算 结果 如 表 ２ ＮＹ 计 。
应力强度因子的数值计算方法
应力强度因子的数值计算方法引言一、理论计算方法1.弹性理论解法弹性理论解法是应力强度因子计算中最常用的一种方法。
它假设材料是弹性线性的，并忽略了材料的塑性变形。
常用的解法有Westergaard解和Westergaard-Hankel解。
2.能量解法能量解法是一种基于弹性力学的解法，通过计算裂纹尖端处的应力场能量和应变能量来计算应力强度因子。
常用的解法有Line-spring法和Irwin法。
3.有限元法有限元法是一种数值计算方法，通过将复杂的问题离散化为多个小区域，并在每个小区域上建立适当的数学模型进行计算。
通过求解离散化的方程组，可以得到裂纹尖端处的应力强度因子。
有限元法可以处理各种复杂的边界条件和几何形状的问题，并且可以考虑非线性和塑性变形。
这使得它成为计算应力强度因子的一种重要方法。
二、实验计算方法实验计算方法主要是通过设计和进行试验来测量裂纹尖端区域的应力和应变场，然后根据测量数据计算应力强度因子。
常用的方法有：1.发光全场法发光全场法是一种全场应变测量技术，通过在被测结构表面涂覆一层发光材料，然后利用高速摄像机记录结构在加载过程中的应变分布。
通过分析图像数据，可以得到裂纹尖端区域的应力和应变场，进而计算应力强度因子。
2.特征裂纹法特征裂纹法是一种利用疲劳试验得到应力强度因子的方法。
通过在试样上开几何形状确定的裂纹，然后在加载过程中观察裂纹的扩展行为，通过测量裂纹长度和加载荷载的关系，可以计算应力强度因子。
3.数值模拟法数值模拟法是一种将实验和数值计算相结合的方法。
通过建立几何和材料特性相似的数值模型，并在模型中模拟加载过程，可以得到裂纹尖端区域的应力和应变场，进而计算应力强度因子。
三、应力强度因子的应用1.疲劳断裂评估基于应力强度因子的计算结果，可以对工程结构在疲劳载荷下的断裂寿命进行评估和预测。
这对于提高结构的可靠性和安全性具有重要意义。
2.材料断裂韧性评定3.裂纹扩展行为研究通过分析应力强度因子的变化规律，可以研究裂纹在不同加载条件下的扩展行为，揭示断裂的机理和规律。
断裂力学的关键参数-应力强度因子
小刘-LZP08-07原文材料或构件中存在宏观裂纹，这些裂纹产生的原因一般为如下几个方面：应力强度因子是表征材料断裂的重要参量，是表征外力作用下弹性物体裂纹尖端附近应力场强度的一个参量。
1957年, 欧文(Irwin,G.R.)建立了以应力强度因子为参量的裂纹扩展准则——应力强度因子准则,从而成功地解释了低应力脆断事故。
应力强度因子的概念：应力强度因子是断裂力学在研究应力作用下考虑应力和裂纹尺寸这两个因素对裂纹扩展影响而引入的新参数，记为K，它反映了裂纹顶端附近应力强弱程度。
对于普通的构件，一般形状的裂纹应力强度因子属于KⅠ型。
应力强度因子与作用在构件裂纹顶端处的名义应力σ及裂纹尺寸α之间存在如下的普遍关系。
上式中的Y为表征含裂纹构件几何形状的一个无因次系数。
应力强度因子的分类：对于不同的裂纹扩展类型有不同的应力强度因子。
可以用下图表示：K1,K2,K3,分别对应于张开型,滑开型和撕开型裂纹的应力强度因子。
张开型(Ⅰ型)裂纹应力强度因子KⅠ是线弹性断裂力学中一个重要断裂参量.设外载和结构均以裂纹2a为对称。
工程上Ⅰ型裂纹出现的最多，最危险，研究最深入。
是低应力脆断的主要原因。
应力强度因子的应用：由张开型的应力强度因子表达式可以看出，KⅠ仅由裂纹长度和名义应力确定。
若已知裂纹长度和名义应力，则KⅠ为定值，并确定了裂纹能否扩展。
由此，我们可以用KⅠ来建立某个条件并判断构件的裂纹是否扩展。
比如，某一有一个2α长度的穿透裂纹的平板，在均匀拉应力作用下，KⅠ值随外应力增大。
当外应力σ增大到一定程度时，裂纹达到失稳状态，此时，即使外力不再增加，裂纹也会迅速扩展，直到断裂。
这说明此时材料已达到KⅠ的极值。
这个极值称为材料的断裂韧性，记为KⅠc。
可见，KⅠc 表示的是材料的一种力学性能，它与试件的几何形状、受力情况、试验环境以及加载方式等有关，其值可以用试验测定。
显而易见，带裂纹的零部件产生脆断的临界条件为：上式称为脆性断裂判断式，即说明当带张开型裂纹的机械零件的应力强度，因子KⅠ达到断裂韧性KⅠc时，零件即断裂。
i型应力强度因子为负值
i型应力强度因子为负值
I型应力强度因子为负值可能表示应力场与裂纹面不垂直，或者裂纹面与受力方向不垂直。
在断裂力学中，I型应力强度因子是描述裂纹尖端应力场强度的参数，其值的大小和符号可以反映裂纹尖端的应力状态和裂纹扩展的趋势。
如果I型应力强度因子为负值，则表示裂纹尖端的应力场与裂纹面不垂直，即裂纹面与受力方向不垂直。
这种情况下，裂纹尖端的应力场可能存在一定的剪切分量，这可能会导致裂纹发生剪切断裂。
另外，如果I型应力强度因子为负值，也可能是由于计算过程中使用的参数不正确或者模型简化不当导致的。
因此，在计算过程中需要仔细检查和校验所使用的参数和模型。
需要注意的是，对于具体的工程应用，需要综合考虑不同类型和方向的裂纹对结构性能的影响，进行详细的断裂分析和评估，以确保结构的安全性和可靠性。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
In stress analysis each point, x,y,z, of a stressed solid undergoes the stresses; sx sy, sz, txy, txz,tyz. With reference to figure 2.3, when a body is loaded and these loads are within the same plane, say the x-y plane, two different loading
s
contraction of lateral surfaces
X
occurs, and, a
2. plane strain (PSN), when the
Crack Plane
sz sz sz sz
specimen is thick enough to avoid contraction in the
断裂力学 第三讲
裂纹尖端应力场，应力强度因子
问题
• 有哪几种典型的裂纹扩展形式？ • I、II、III型裂纹的裂纹尖端应力应变计位移的分布形式？ • 为什么用裂纹强度因子可以表征线弹性材料的断裂过程？ • 裂纹强度因子与能量释放率之间的关系
弹性理论简述
• 应力；
• Airy应力函数；
• 变形；
• Airy函数的复变函数
Mode I - Opening mode: where the crack surfaces separate symmetrically with respect to the plane occupied by the crack prior to the deformation (results from normal stresses perpendicular to the crack plane); Mode II - Sliding mode: where the crack surfaces glide over one another in opposite directions but in the same plane (results from in-plane shear); and Mode III - Tearing mode: where the crack surfaces are displaced in the crack plane and parallel to the crack front (results from out-of-plane shear).
crack length, applied stress s and shape factor Y( which accounts for finite size
of the component and local geometric features).
The Airy stress function.
conditions are possible:
s
1. plane stress (PSS), when the
thickness of the body is
comparable to the size of the
y syy
Thickness
B
s
Thickness B
s
plastic zone and a free
sz = ν (sx+sy) where ν is Poisson's ratio.
• 应变；
表示形式；
• 应力应变关系；
• Westergaard应力函数；
• •
平衡方程；
•
平面应力和平面应变问题；
三种裂纹尖端的线性 弹性应力场解；
• 二维问题弹性解；
LINEAR ELASTIC FRACTURE MECHANICS (LEFM)
For LEFM the structure obeys Hooke’s law and global behavior is linear and if any local small scale crack tip plasticity is ignored The fundamental principle of fracture mechanics is that the stress field around a crack tip being characterized by stress intensity factor K which is related to both the stress and the size of the flaw. The analytic development of the stress intensity factor is described for a number of common specimen and crack geometries below. The three modes of fracture
s
s
thickness z-direction.
s
Plane Stress
Plane Strain
In the former case, the overall stress state is reduced to the three
components; sx, sy, txy, since; sz, txz, tyz= 0, while, in the latter case, a normal stress, sz, is induced which prevents the z displacement, ez = w = 0. Hence, from Hooke's law:
LINEAR ELASTIC FRACTURE MECHANICS (Contd.)
In the 1950s Irwin [7] and coworkers introduced the concept of stress intensity
factor, which defines the stress field around the crack tip, taking into account