abaqus应力强度因子计算的例子
ABAQUS热应力分析实例详解
热应力分析实例详解学习要点通过实例分析，学习如何进行热应力分析，并掌握ABAQUS/CAE 的以下功能：1）在Material 功能模块中，定义线胀系数；2）在Load 功能模块中，使用预定义场（predefined field）来定义温度场；实例1：带孔平板的热应力分析定义材料属性——Property Property——Material——Edit——steelMechanical——Elastic, 输入弹性模量和泊松比定义材料属性——Property Property——Material——Edit——steelMechanical——Expansion, 输入线胀系数定义边界条件——Load定义边界条件——Load定义边界条件——Load固支边界条件使用预定义场定义初始温度Load——PredefinedField Manager使用预定义场使模型温度升高至120℃网格划分——Mesh结果分析——Visualization小结在ABAQUS中进行热应力分析的基本步骤：⏹定义线胀系数⏹定义初始温度场⏹定义分析步中的温度场实例2：法兰盘感应淬火的残余应力场模拟问题描述：◆表面感应淬火是一种工程中常用的热处理工艺，其原理是使用感应器来对工件的局部进行加热，然后迅速冷却，从而使工件表面产生残余压应力，抵消工作载荷所产生的一部分拉应力。
◆表面感应淬火可显著提高工件弯曲疲劳抗力和扭转疲劳抗力，工件表面产生的马氏体具有良好的耐磨性。
实例2：法兰盘感应淬火的残余应力场模拟 本例中的法兰盘经淬火后，由试验测得法拉盘的内圆角表面残余压应力约为-420MPa。
法拉盘的一端固定，另一端的整个端面受向下的面载荷p=100MPa，法拉盘内孔直径为24mm，材料的弹性模量为210000MPa，泊松比为0.3，线胀系数为1.35e-5/ ℃。
要求：模拟分析感应淬火所产生的残余应力场，并分析此残余应力场在缓和应力集中方面所起的作用。
使用ABAQUS计算应力强度因子
------------------------------------------------------------------------------------------------------- 如何使用ABAQUS计算应力强度因子Simwefanhj（fanhjhj@）2011.9.9------------------------------------------------------------------------------------------------------- 问题描述：以无限大平板含有一贯穿裂纹为例，裂纹长度为10mm（2a），在远场受双向均布拉应力σ＝100N/mm2。
按解析解，此I型裂纹计算出的应力＝396.23（N.mm-3/2）强度因子πσaK=I以下为使用ABAQUS6.10的计算该问题的过程。
第一步：进入part模块①建立平板part（2D Planar;Deformation;shell），平板的尺寸相对于裂纹足够大，本例的尺寸为100×50（mm）。
②使用Partation Face:sketch工具，将part分隔成如图1形式。
图1第二步：进入property模块①建立弹性材料；②截面选择平面问题的solid，homogeneous；③赋予截面。
第三步：进入Assembly模块不详述。
需注意的是：实体的类型（instance type）选择independent。
第四步：进入mesh模块除小圈内使用CPS6单元外，其它位置使用CPS8单元离散（图2）。
裂纹尖端的奇异在interaction模块中（图4）考虑。
图2第五步：进入interaction模块①指定裂纹special/creak/assign seam，选中示意图3中的黄色线，done！②生成裂纹crack 1，special/crack/create,name：crack 1，type: contour integral.当提示选择裂纹前端时，选则示意图的红圈区域，当提示裂纹尖端区域时选择红圈的圆心，用向量q表示裂纹扩展方向（示意图3绿色箭头）。
第二章 应力强度因子的计算
第二章 应力强度因子的计算K --应力、位移场的度量⇒K 的计算很重要,计算K 值的几种方法: 1.数学分析法:复变函数法、积分变换； 2.近似计算法：边界配置法、有限元法； 3.实验标定法：柔度标定法； 4.实验应力分析法：光弹性法.§2-1 三种基本裂纹应力强度因子的计算一、无限大板Ⅰ型裂纹应力强度因子的计算K Z ξ→=→ⅠⅠ计算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,cos1222bbbπππξξξ==⇒sin()sincos 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 值之比。
使用ABAQUS计算应力强度因子
使用ABAQUS计算应力强度因子应力强度因子（Stress Intensity Factor，简称SIF）是应力场的一种特征参数，用于描述应力状态下混合模式断裂的倾向性。
它在断裂力学和疲劳断裂力学中起着非常重要的作用。
在ABAQUS软件中，可以通过线性弹性断裂力学方法来计算应力强度因子。
ABAQUS中计算SIF的方法通常分为两步：1.求解应力场2.计算SIF在求解应力场时，可以采用以下几种途径：1.固定边界条件：如果边界条件已知并且不会发生变化，则可以直接固定边界条件来求解应力场。
这种方法适用于简单的几何形状和加载情况。
2.施加约束：对于复杂几何形状和加载情况，可以施加约束来求解应力场。
例如，可以在加载边界上施加位移或力，并在其他边界上施加自由边界条件。
ABAQUS软件将通过求解线性弹性方程来获得应力场。
3.等效边界法：对于无法通过上述两种方法求解应力场的情况，可以采用等效边界法。
该方法将复杂几何体简化为等效的几何体，通过在等效边界上施加约束来求解应力场。
然后，可以使用所得的应力场计算SIF。
在计算SIF时，可以采用两种方法：1.J积分方法：这是一种基于应变能的方法，通过计算闭合路径上的应力和应变来计算SIF。
ABAQUS提供了J积分的计算方法，可以直接计算SIF。
2.基于位移法：这是一种基于位移的方法，通过计算表面位移场的奇异性来计算SIF。
ABAQUS也提供了这种方法的计算选项。
计算SIF的步骤一般如下：1.定义几何模型和输入材料参数。
2.设置边界条件和加载条件。
3.运行ABAQUS求解应力场。
4.运行相应的计算器（如J计算器或位移计算器）以计算SIF。
5.根据得到的SIF结果进行进一步的断裂力学分析。
需要注意的是，计算SIF是一个相对复杂的过程，需要对模型几何形状、边界条件、加载条件和材料参数等进行仔细考虑和设置。
此外，模型的网格划分和数值求解的精度也会对计算结果产生影响，因此需要进行适当的验证和后处理分析。
基于ABAQUS的货叉三维裂纹应力强度因子有限元分析
基于ABAQUS的货叉三维裂纹应力强度因子有限元分析货叉是一种常用于起重机械的重要零件，承受着大量的动态和静态荷载。
在使用过程中，货叉可能会受到裂纹的影响，从而降低其强度和安全性。
因此，对货叉的裂纹应力强度因子进行分析是非常必要的。
裂纹应力强度因子是评估裂纹尖端应力场的参数，它可以用来判断裂纹的扩展情况以及材料的断裂行为。
基于ABAQUS的有限元分析可以用来计算货叉在裂纹尖端处的应力强度因子。
该分析要求以下几个步骤：1. 建立货叉的三维有限元模型：模型要包括真实的几何形状和材料性质。
可以使用ABAQUS提供的建模工具，如Part模块和Assembly模块，来构建模型。
此外，还需考虑货叉的边界条件和加载方式。
2.设置裂纹：在模型中引入裂纹，它可以是表面裂纹或体内裂纹。
可以使用ABAQUS提供的功能来创建裂纹和裂纹前沿。
3.划分网格：为了计算裂纹应力强度因子，需要划分网格并分配单元类型和单元属性。
合理的网格划分可以提高计算精度和效率。
4.应用荷载：根据实际情况，在模型中施加与实际工作状况相对应的荷载。
荷载类型可以包括静态荷载、动态荷载或者其他较为复杂的荷载。
5.运行分析：设置好所有必要的计算参数后，可以运行分析并计算货叉的裂纹应力强度因子。
6.结果分析：根据计算结果，可以评估货叉中裂纹的状态和扩展情况。
一般来说，如果裂纹应力强度因子超过了材料的断裂韧性，则裂纹有可能扩展，从而降低货叉的强度和安全性。
在进行有限元分析时，需要注意模型的合理性和准确性。
同时，还应考虑到材料的非线性特性和可能的影响因素，以获得较为准确的分析结果。
总之，基于ABAQUS的货叉三维裂纹应力强度因子有限元分析可以用来评估货叉中裂纹的状态和扩展情况，为提高货叉的安全性和可靠性提供科学依据。
含圆孔和裂纹板应力强度因子分析
《断裂力学》大作业题目：含圆孔和裂纹板应力强度因子分析姓名：学号：专业：授课教师：一、问题描述含多裂纹矩形板受垂直方向拉伸载荷作用，如图 1 所示，计算中心裂纹尖端的应力强度因子KⅠ和KⅡ，并讨论其随即和参数L、h、a、D、 等的变化规律，写一篇分析报告。
图1. 含三条裂纹矩形板受垂直拉伸载荷作用要求（1）报告中计算所用到的分析方法和模型应阐述清楚，并写出必要的计算公式。
（2）绘制应力强度因子随几何参数的变化曲线。
（3）列出必要的参考文献二、理论分析在线弹性断裂力学中，I型裂纹尖端的应力场为：(1sin sin)222(1sin sin)222cos cos222333xyxyσστθθθθθθθθθ⎧=-⎪⎪⎪=+⎨⎪⎪=⎪⎩I型裂纹尖端的位移场为：1)cos(1cos)221)sin sin2233uvκκθθθθ⎧=--⎪⎪⎨⎪=+⎪⎩其中：3431νκνν-⎧⎪=⎨-⎪+⎩平面应变平面应力同理，对II型裂纹尖端的应力场：(2cos cos)222cos sin cos222(1sin sin)222333xyxyσστθθθθθθθθθ⎧=+⎪⎪⎪=⎨⎪⎪=-⎪⎩显然，位移场和应力场均可以表示成应力强度因子的形式。
通过对裂纹尖端的应力应变场分析来求解对应的应力强度因子，便是传统有限元求解应力强度因子的原理。
而对于I、II复合型裂纹尖端的应力强度因子，可通过它们的叠加获得。
确定应力强度因子的方法有3大类：解析法、数值解法和实验方法。
解析法只能计算简单问题，大多数问题需要采用数值解法，当前工程中广泛采用的数值解法是有限单元法。
随着有限元法的发展，有限元在断裂力学中的应用越来越普及。
近年来，计算机技术得到了迅猛发展，许多功能强大的有限元软件也相继问世，大型通用有限元程序abaqus就是当前工程中应用最广泛的有限元软件之一。
采用abaqus软件计算裂纹尖端的应力强度因子，通过阅读abaqus 的帮助文件，得到abaqus基于有限元方法在线弹性范围内计算应力强度因子的原理。
Abaqus 齿轮力和应力计算
Ken Youssefi
Mechanical Engineering Dept.
14
Surface Strength Analysis
The basic surface deterioration Scoring
If the surface asperity welding and tearing cause a transfer of metal from one surface to the other, the resulting surface damage is called scoring. If the local welding of asperities becomes so extensive that the surfaces no longer slide on each other, the resulting failure is called seizure. Initial scoring on 4340 steel helical gear Moderate scoring on a 3310 steel spur gear.
Wt
F
Substituting for x and introducing p (circular pitch),
The form factor y is called Lewis form factor. Substituting P = / p and Y = y
Lewis’ equation, where
Ken Youssefi Mechanical Engineering Dept.
6
Modification of Lewis’ Equation
Assumptions made in deriving Lewis’ equation
abaqus截面拉应力积分
Abaqus截面拉应力积分1. 简介Abaqus是一款强大的有限元分析软件，可用于模拟和分析各种结构和材料的力学行为。
在工程实践中，经常需要对截面进行应力分析，以评估结构的性能和安全性。
Abaqus提供了一种方便的方法来计算截面上的拉应力积分，从而得到准确的结果。
本文将介绍如何使用Abaqus进行截面拉应力积分，并提供详细步骤和示例。
2. 模型准备在进行截面拉应力积分之前，首先需要准备好模型。
模型可以通过Abaqus提供的建模工具创建，也可以导入其他CAD软件中创建好的模型。
下面以一个简单的梁模型为例来说明：```这是一个简单的矩形梁模型，我们将对其截面进行拉应力积分。
3. 定义切割区域在Abaqus中，我们可以通过定义切割区域来指定要进行拉应力积分的区域。
切割区域定义了一个封闭曲线或多边形，它将被用作拉应力积分的边界。
在本例中，我们将定义一个矩形切割区域，以便对梁的截面进行拉应力积分。
以下是定义切割区域的步骤：1.打开Abaqus软件并加载模型。
2.在”Part”模块中选择要进行拉应力积分的截面。
3.进入”Part”模块的编辑模式。
4.在工具栏上选择”Create Sketch”工具，并在截面上创建一个矩形。
5.使用”Constraints”工具将矩形的边界与截面对齐。
6.完成矩形的创建后，选择它并使用右键菜单中的”Properties…“选项打开属性窗口。
7.在属性窗口中，将切割区域类型设置为”Integral Measure”.8.确定好切割区域后，退出编辑模式并保存模型。
4. 设置拉应力积分在定义好切割区域后，接下来需要设置拉应力积分。
以下是设置拉应力积分的步骤：1.在Abaqus主界面上选择”Job Manager”.2.右键单击要进行拉应力积分的作业，并选择”Edit Job…“选项。
3.在作业编辑器中，选择要进行拉应力积分的截面部分，并使用右键菜单中的”Create Field Output…“选项。
基于ABAQUS的两种应力强度因子计算方法对比
21
50.363 2
5
282.285 291 2
25
1 411.426 456
22
49.705
5.25
285.476 037 7
27.562 5
1 498.749 198
23
49.082 3
5.5
288.533 458 1
30.25
1 586.934 019
24
48.524 8
5.75
291.667 209 7
5.062 5
740.943 363 2
11
74.684 7
2.5
295.999 907 7
6.25
739.999 769 2
12
65.193 5
2.75
270.994 561 2
7.562 5
745.235 043 2
13
59.114 6
3
256.652 390 1
9
769.957 170 4
14
，则其截距
2 ABAQUS中两种SIF计算方法的实 现[5-=]
2.1 ABAQUS内置方法计算SIF
ABAQUS 中建立半宽,=100 ??，半高-=200 mm，
厚度1 mm的平板模型（ 2），该平板中心预制有一
条裂纹，其半裂纹长度 20 mm; 弹性模量 2!
105MPa，泊松 0.25的“弹性各向异性”固体材料；
31
K FACTOR ESTIMATES
CRACK CRACKFRONT
CONTOURS
(平均:75%)
NAME
NODE SET
+1.054e+03
abaqus算例
ABAQUS实例操作吴敏2013202100035一、型钢梁建模分析1.1 问题描述一型钢梁，尺寸如图所示，利用软件分析其内力。
材料特性：弹性模量E=2.1e11N/m2,泊松比μ=0.3，屈服强度ƒy=3.45e8N/m2。
1.2创建部件点击创建部件按钮，在对话框中设置参量如右图：模型空间设置为三维的，类型为可变性的，基本特征为实体，可拉伸，比例设为1.1.2生成三维模型首先，在二维的环境下，输入横截面的各点坐标，然后再输入深度6m，便可生成如下图型：1.3创建材料和截面属性1.3.1创建材料先输入弹性模量，泊松比，以及屈服应力，塑性应变，点击确认即可。
1.3.2 创建截面属性并赋予给部件名字命名为section-beam，种类为实体，类型为均质，其他值保持默认，点击确认，接着选择整个部件，将截面性质赋予之。
1.4 定义装配件点击装配功能模块，选择部件为非独立实体其他保持为默认值点解确认即可。
1.5 设置分析步选择分析步模块，点击create instance 在对话框里面，输入名字为step,procedure type设置为general ，在下拉菜单中选择static general项，保持其他参数不变，点击确认。
1.6 定义荷载和边界条件选择荷载模块：①施加荷载在create load对话框中，名字设置为load,step项中选择为step,将荷载设置为pressure，其他值保持不变，点击继续，在荷载的大小后面输入3.5e5，其他参数不变，完成荷载的定义。
②定义便捷条件在对话框中将step 设置为initial,将施加边界条件的方式设置为位移/转角，保持其余参数不变，点击确认。
在弹出的对话框中选择U1=U2=UR2=UR3=0，即对选中面施加铰接约束，点击ok。
同样的方式在另一边同样设置。
1.7 划分网格在列表中选择功能模块，对模型进行网格划分，将环境栏中的object项设为part，即为部件划分网格。
ABAQUS计算J积分细节
ABAQUS计算J积分细节Abaqus计算J积分，主要是指派裂纹及定义裂纹的方向，同时在step中的历程模型：10×50裂纹：定义一个尖端；另一方面指明裂纹的扩展方向（1.0,0.0），说白了就是X积分数值,图中输入10，则输出10个J积分值计算包含裂纹尖端的包络区域的面积即为J积分避开裂纹尖端塑形区域的不可计算的特性。
同时J积分的计算数值与积分路径无从以上图例可以看出J积分数值区域稳定。
疑问：为何计算多个积分点，是否最后的稳定数值就是需要计算的J积分数值？J积分应该是数值，而不是多个不同的数值。
我个人觉得最后的稳定数值应该是需要计算的积分数值。
从dat文件到inp文件，找到积分区域。
pickseted12以及pickseted13都是节点4，坐标如图所示，在cad模型中的位置如箭头指向，即裂纹尖端。
详细的需要看一下abaqus帮助文档，关于J积分的计算细节。
积分点的个数的意义我还没有搞清楚。
对于应力强度因子K，表征裂纹尖端受力的一个参量，在裂纹尖端的应力场的一定范围内，不同的节点计算数值大体是相同的。
计算应力强度因子：可以利用abaqus直接输出，也可以利用公式计算应力强度因子，以下为利用有限元法计算应力强度因子：计算应力强度因子：从上图可以看出，计算应力强度应力的点与计算J积分的点是一致的。
Abaqus计算的应力强度因子为裂尖处的应力强度因子。
下面我们利用有限元法计算y=0处的应力强度因子，最后外推到裂尖处的应力强度因子。
选取不同的半Abaqus计算J积分注意事项：(1)一、Interaction模块1.1 预制裂纹（步骤：菜单/special/crack/assign seam）注意：并不是作裂纹分析都要定义seam，如果你的裂纹不是一条缝，而是一个缺口，则不需要assign seam，直接走下一步（定义裂纹）就行。
1.2 创建裂纹（步骤：菜单/special/crack/create，type：contour integral）—crack front：crack front是用来定义第一围线积分的区域，2D下我们可以选择包围裂尖点的面，3D则选择包围裂尖线的面；另外还有一种定义crack front的方法，就是直接选择裂尖点（2D）或裂尖线3D），用这个方法定义crack front不需要再定义下一步的crack tip/line，比较简便，两种方法算出的结果没有明显的差别，其实只是影响积分路线的问题，但是J积分值是路径无关的，看个人喜好吧—crack tip/line：这个比较好理解就是裂尖点（2D）或线（3D），如果我们在上一步中用方法二定义crack front，这一步就直接跳过了—crack extension direction（定义裂纹扩展方向）：这里定义的其实是一个虚拟的裂纹扩展方向，定义了这个参考方向后，我们才能通过输出的角度判断裂纹扩展方向，可以通过两种方法：（1）q vector：输入一个方向，用来作为计算裂纹的扩展方向的参考方向；（2）normal to crack plane：crack plane表示裂纹的对称面（当裂纹在一个平面内时，可能需要分开定义多个裂纹），这种方法下我们只需定义裂纹面的法线方向，通过（t表示裂纹尖端的切线）, 会在每个节点得出一个q方向；（3）注意：q的方向对输出的应力强度因子，J积分等都会有影响，一般情况下，q最好在裂纹平面内，且垂直于裂尖线的切线，否则算出的应力强度因子，J积分值等等在不同围线积分中会差别较大。
ABAQUS计算J积分细节
ABAQUS计算J积分细节Abaqus计算J积分，主要是指派裂纹及定义裂纹的方向，同时在step中的历程输出变量中定义J积分。
如下图所示：模型：10×50定义裂纹：裂纹：定义一个尖端；另一方面指明裂纹的扩展方向（1.0,0.0），说白了就是X轴的方向。
Step中输出裂纹参量：积分数值,图中输入10，则输出10个J积分值J积分的理解：计算包含裂纹尖端的包络区域的面积即为J积分尖端塑形区域的不可计算的特性。
同时J积分的计算数值与积分路径无关。
以下为六个计算从以上图例可以看出J积分数值区域稳定。
疑问：为何计算多个积分点，是否最后的稳定数值就是需要计算的J积分数值？J积分应该是数值，而不是多个不同的数值。
我个人觉得最后的稳定数值应该是需要计算的积分数值。
从dat文件到inp文件，找到积分区域。
pickseted12以及pickseted13都是节点4，坐标如图所示，在cad模型中的位置如箭头指向，即裂纹尖端。
详细的需要看一下abaqus帮助文档，关于J积分的计算细节。
积分点的个数的意义我还没有搞清楚。
对于应力强度因子K，表征裂纹尖端受力的一个参量，在裂纹尖端的应力场的一定范围内，不同的节点计算数值大体是相同的。
计算应力强度因子：可以利用abaqus直接输出，也可以利用公式计算应力强度因子，以下为利用有限元法计算应力强度因子：计算应力强度因子：从上图可以看出，计算应力强度应力的点与计算J积分的点是一致的。
Abaqus计算的应力强度因子为裂尖处的应力强度因子。
下面我们利用有限元法计算y=0处的应力强度因子，最后外推到裂尖处的应力强度因子。
选取不同的半径r值计算得到的应力Abaqus计算J积分注意事项：(1)一、Interaction模块(2) 1.1 预制裂纹（步骤：菜单/special/crack/assign seam）(3)注意：并不是作裂纹分析都要定义seam，如果你的裂纹不是一条缝，而是一个缺口，则不需要assign seam，直接走下一步（定义裂纹）就行。
abaqus强度折减法算例
abaqus强度折减法算例ABAQUS STRENGTH REDUCTION FACTOR EXAMPLE应力折减系数（strength reduction factor,SRF），简称力学折减系数，是用来精确表示在较低的应力水平下材料强度降低的一个量化概念。
它是应力折减技术的重要参数，常用来预测材料和结构在损伤情况下的抗压强度及机构的稳定性。
它实质上是材料应力变形和损伤特性曲线上的斜率，用于评估在特定应力水平下材料的强度折减程度。
Strength reduction factor (SRF), also known as mechanical reduction factor, is a quantitative concept used to accurately represent the reduction in material strength at lower stress levels. It is an important parameter of stress reduction technology, which is often used to predict the compressive strength of the material and the stability of the mechanism under damage conditions. It is essentially the slope of the stress-strain and damage characteristics curve used to evaluate the degree of strength reduction of the material ata specific stress level.1. ABAQUS 力学折减系数例1. ABAQUS Mechanical Reduction Factor Example假设材料的杨氏模量为E，极限应力为σy，ABAQUS中可以使用力学折减系数来评估在低应力下的材料抗压强度。
基于ABAQUS的复杂裂纹应力强度因子的研究
基于ABAQUS的复杂裂纹应力强度因子的研究摘要：在实际工程领域，裂纹问题一直都是一个影响结构安全性的难题。
本文从现实工程角度出发，针对复杂裂纹的断裂破坏机理，采用基于ABAQUS的数值模拟技术，进行了深入的研究。
为了证明采用数值模拟技术研究复杂裂纹的可靠性，本文首先利用ABAQUS对其他学者做出的理论值进行了验证，结果表明理论和数值具有高度的吻合性。
然后进行了复杂裂纹在不同情况下的裂尖应力强度因子进行了研究。
从而为现实工程建设中一些不便于实验研究的复杂情况，提出了研究的可行性方法，具有较强的实践意义。
关键词：裂纹；应力强度因子；ABAQUS；数值模拟0 引言断裂力学是研究带裂纹体的强度和裂纹扩展规律的一门学科。
断裂力学的最早理论可以追溯到1920年，为了研究玻璃、陶瓷等脆性材料的实际强度比理论强度低的原因，Griffith[1-2]提出了在材料中存在裂纹的设想，而Irwin[3-4]在1957年提出了应力强度因子以及其后形成的断裂韧度的概念后，断裂力学理论出现了重大的突破，奠定了线弹性断裂力学的基础。
当前国内外许多学者[5-8]已将多裂纹问题作为主攻方向，但是受到诸多条件的制约，目前并未取得明显的成果。
1复杂裂纹研究的现实依据图（a）是位于四川的几大主要断裂带的分布图，在宏观上形成了V型的断层结构。
汶川Ms8.0级地震，就是龙门山断裂带在构造应力场长期作用下，积累的应变能突然释放的结果。
图（b）是2002年11月3日，阿拉斯加中部迪纳利断层发生了一次M w7.9级地震。
从图（b）可以看出三条主断裂层形成了典型的的分支结构，究断层的性质对于预知、了解地质作用，估计可能发生的破坏具有重要参考价值。
图1 现实V型、Y型裂纹存在模型2 Y型裂纹的应力强度因子研究2.1有限元建模中间带有裂纹的双轴压缩试件，如图2所示，在相关断裂力学的试验分析中被广泛应用。
故此，本文以此典型试件作为研究对象，进行数值模拟计算。
应力强度因子的计算
应力强度因子的计算应力强度因子（Stress Intensity Factor）是应用于裂纹尖端的一个参数，用于描述裂纹尖端应力场的强度和分布情况，是计算裂纹扩展速率和破裂韧性的重要参数。
本文将详细介绍应力强度因子的计算方法。
一、引言在构件中存在裂纹时，应力场的分布将发生变化，通常存在一个应力集中区域，即裂纹尖端。
在裂纹尖端附近，裂纹两侧的应力强度具有很大的梯度，因此需要引入应力强度因子来准确描述和分析裂纹尖端的应力状态。
二、应力强度因子的定义应力强度因子可以描述裂纹尖端应力场的强度和分布情况。
对于模式I或拉应力模式下的裂纹，应力强度因子K是一个标量，具有长度的物理意义。
对于一种给定的应力场，应力强度因子K与应力强度因子K对应的应力场是相似的。
此外，由于应力强度因子K的引入，裂纹尖端附近的应力场能够用一个等效应力来代替，从而使裂纹尖端的破坏准则能够使用等效应力来描述。
三、常用的计算方法1.解析方法解析方法是通过对裂纹尖端附近应力场的数学分析，推导出裂纹尖端的应力强度因子。
常用的方法有：格里菲斯公式、韦尔奇定理、赵万江公式等。
这些方法通常需要对裂纹尖端应力场进行严格的数学推导和分析，适用于简单几何形状的裂纹。
2.应力分析方法应力分析方法是通过有限元分析、边界元分析等数值方法，对裂纹附近的应力场进行数值模拟，进而计算应力强度因子。
通过数值模拟可以得到更为复杂的几何形状下的应力强度因子。
通常需要使用计算机软件进行模拟和计算。
3.基于实验的方法基于实验的方法是通过实验测定裂纹尖端的应力强度因子，从而得到一种实验估算的方法。
常用的实验方法有高约束比压缩试验法、断口法、几何函数法等。
与解析方法和数值方法相比，实验方法具有直接、可靠、全面的优点，但通常对实验设备和技术要求较高。
四、应力强度因子的应用应力强度因子的计算在材料科学、工程结构分析和破坏力学等领域具有广泛的应用价值。
它可用于计算裂纹扩展速率、破断韧性、疲劳寿命等。
ABAQUS计算裂纹尖端应力强度因子有效性的算例研究
ABAQUS计算裂纹尖端应力强度因子有效性的算例研究发表时间：2018-09-11T11:34:12.223Z 来源：《新材料.新装饰》2018年3月下作者：汪波[导读] 在实际工程领域中，相当部分的脆性材料总是不可避免的存在着裂纹或是缺陷。
在实际环境中材料的受力往往是相当复杂的。
基于ABAQUS平台的裂纹仿真软件，它具有简单易用的特点。
（成都理工大学工程技术学院，四川乐山 614000）摘要：在实际工程领域中，相当部分的脆性材料总是不可避免的存在着裂纹或是缺陷。
在实际环境中材料的受力往往是相当复杂的。
基于ABAQUS平台的裂纹仿真软件，它具有简单易用的特点。
通过算例分析验证表明，该软件的计算结果具有较高的精度，完全可以用于实际工程问题的计算，通过分析验证表明该软件的设计是成功的。
此外，今后可以在它的基础上进行更多功能扩展，从而使它拥有分析更为复杂问题的能力。
关键词：裂纹；应力强度因子；断裂力学；ABAQUS引言材料在成型和加工过程中在其内部造成了很多缺陷，而其破坏正好均源于构件内部的微小裂纹，所以研究带裂纹的物体力学性能具有十分重要的意义。
图1存在于岩石和混凝土地面中的裂缝1920年， Griffith[1-2]提出了在材料中存在裂纹的设想，而从Irwin[]3-4]在1957年提出了应力强度因子以及其后形成的断裂韧度的概念后，断裂力学理论出现了重大的突破，奠定了线弹性断裂力学的基础。
1基本原理近年来以数值分析为基础的手段来解决断裂力学相关问题的技术得到了广泛的发展应用，并且不断的调整完善。
该技术在一定程度上较好的克服了实验条件下的不足。
对于线弹性断裂力学而言，裂尖区域的位移场、应力、应变场由应力强度因子决定，故而通过有限元计算的结果来得到具体的应力强度因子的值是线弹性断裂力学中用有限元法的基本要求。
1.1 ABAQUS求解裂纹尖端的应力强度因子传统的有限元在计算裂纹尖端的应力强度因子的时候，无可避免地遇到裂尖复杂应力场和位移场的计算，J积分则可以完全避免这种复杂的处理过程。
有限元Abaqus实例讲解——压力容器不连续区域应力分析
Abaqus实例讲解——压力容器不连续区域应力分析问题描述：某高压容器设计压力为P=16MPa，材料的弹性模量为2e5MPa,泊松比为0.3。
筒体内径为R1=775mm，壁厚t=100mm，封头内径R2=800mm，厚度t2=48mm，筒体削边长度L=95mm。
对该高压容器筒体与封头连接区域进行应力分析。
注：本实例所使用的软件为Abaqus6.10-1实例来源：中国水利水电出版社，赵腾伦，《ABAQUS6.6在机械工程中的应用》运行ABAQUS/CAE，修改模型名称为Model-1为pressure-vesel，保存模型为pressure-vessel.cae。
1.双击模型树中模型pressure-vessel下的，弹出Create Part对话框，在Modeling Space 栏中选择Axisymmetric，Type栏中选择Deformable，Base Feature栏中选择Shell，Approximate size中输入4，单击Continue按钮，进入草图环境。
2.单击工具箱中的，在提示区中依次输入下列坐标值：（0.775,1.2）、（0.775,0）、（0.875,1.2）、（0.875,1.2），每输入一个点的坐标以后回车，单击鼠标中毽。
3.单击工具箱中的，在提示区一次输入下列坐标值（输入完一个坐标值回车一次）：（0.8,1.2）、（0.848,1.2）、（0.0,2.0）、（0.0,2.048）、（0.0,1.2），创建5个点。
4.单击工具箱中的，选择点（0.0,1.2）为圆心，以点（0.8,1.2）为起点，点（0.0,2.0）为终点作一条圆弧。
5.同样的方法，以点（0.0,1.2）为圆心，以点（0.848,1.2）为起点，点（0.0,2.048）为终点作另一条圆弧。
6.单击工具箱中的，单击第4步中创建的5个点（多选按Shift键），单击鼠标中键，删除这5个点。
7.单击工具箱中的，连接两个圆弧左上部的端点。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
------------------------------------------------------------------------------------------------------- 如何使用ABAQUS计算应力强度因子
Simwefanhj（fanhjhj@）
2011.9.9
------------------------------------------------------------------------------------------------------- 问题描述：以无限大平板含有一贯穿裂纹为例，裂纹长度为10mm（2a），在远场受双向均布拉应力σ＝100N/mm2。
按解析解，此I型裂纹计算出的应力＝396.23（N.mm-3/2）
强度因子π
σa
K=
I
以下为使用ABAQUS6.10的计算该问题的过程。
第一步：进入part模块
①建立平板part（2D Planar;Deformation;shell），平板的尺寸相对于裂纹足够大，本例的尺寸为100×50（mm）。
②使用Partation Face:sketch工具，将part分隔成如图1形式。
图1
第二步：进入property模块
①建立弹性材料；
②截面选择平面问题的solid，homogeneous；
③赋予截面。
第三步：进入Assembly模块
不详述。
需注意的是：实体的类型（instance type）选择independent。
第四步：进入mesh模块
除小圈内使用CPS6单元外，其它位置使用CPS8单元离散（图2）。
裂纹尖端的奇异在interaction模块中（图4）考虑。
图2
第五步：进入interaction模块
①指定裂纹special/creak/assign seam，选中示意图3中的黄色线，done！
②生成裂纹crack 1，special/crack/create,name：crack 1，type: contour integral.
当提示选择裂纹前端时，选则示意图的红圈区域，当提示裂纹尖端区域时选择红圈的圆心，用向量q表示裂纹扩展方向（示意图3绿色箭头）。
用同样的方法建立crack 2（示意图3中的蓝色区域）。
special/crack/edit，对两个裂纹进行应力奇异的设置，如图4所示。
图3
第六步：进入step模块
①在initial步之后建立static，general步。
②在output/history output requests/create/中创建2个输出变量，如图5所示。
（分别定义crack 1和crack2 ）。
其中Number of contours的意义是：SIF是通过J 积分算出的，J积分本来是线积分，在有限元计算中不方便，后来有人研究可以用环形区域积分代替，Abaqus里J积分的计算就是这样做的，积分围数1就是指计算J积分时取裂尖处奇异单元外第一圈单元作为环形区域算J积分，不同积分围数算出的应力强度因子应该相差不多，如果相差太多，就说明算的有问题。
图5
第七步：进入load模块
定义位移和荷载边界，如图6所示。
图6
第八步：进入job模块
提交计算。
Mises应力分布见图7，在.dat文件中（图8）查看应力强度因子。
图7
图8
可见与I K 解析解的误差仅有2.5％。