在ABAQUS中对应力的部分理解
在ABAQUS中对应力的部分理解
关于abaqus中 mises， s11 s22 s33 ，s12，tresca pressure， max principal，mid principal，min principal。简单地理解，
在ABAQUS中，一般是把X轴当做1轴，Y轴当做2轴，Z轴当做3轴；那么：
S11就是X轴向的应力，正值为拉应力，负值为压应力;
S22就是Y轴向的应力，正值为拉应力，负值为压应力;
S33就是Z轴向的应力，正值为拉应力，负值为压应力;
S12就是在YZ平面上，沿Y向的剪力;
S13就是在YZ平面上，沿Z向的剪力;
S23就是在XZ平面上，沿Z向的剪力;
由于剪力的对称性：S12=S21，S13=S31，S23=S32
Mises应力是即第四强度理论，根据能量守恒原理，用于判断材料是否屈服的应力准则，即Mises准则，一般使用于判断延性比较好的材料，对于脆性材料，一般采用第一强度理论。
使用ABAQUS计算应力强度因子
------------------------------------------------------------------------------------------------------- 如何使用ABAQUS计算应力强度因子 Simwefanhj（fanhjhj@https://www.360docs.net/doc/9b1465030.html,） 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所示。
Abaqus 中显示动力学分析步骤
准静态分析——ABAQUS/Explicit 准静态过程（guasi-static process） 在过程进行的每一瞬间，系统都接近于平衡状态，以致在任意选取的短时间dt内，状态参量在整个系统的各部分都有确定的值，整个过程可以看成是由一系列极接近平衡的状态所构成，这种过程称为准静态过程。无限缓慢地压缩和无限缓慢地膨胀过程可近似看作为准静态过程。准静态过程是一种理想过程，实际上是办不到的。 准静态原为一个热力学概念，在这里引用主要是指模型在加载的过程中任意时刻所经历的中间状态都可近似地视为静力状态，因此当加载过程进行得无限缓慢时，在各个时刻模型所处的状态就可近似地看作是静态，该过程便是准静态过程。准静态啮合过程仿真主要考虑的是弧齿锥齿轮副在加载时的接触状态，以及齿面和齿根的应力变化规律，其前提是不考虑齿轮副惯性的影响。 ABAQUS/Explicit准静态分析 显式求解方法是一种真正的动态求解过程，它的最初发展是为了模拟高速冲击问题，在这类问题的求解中惯性发挥了主导性作用。当求解动力平衡的状态时，非平衡力以应力波的形式在相邻的单元之间传播。由于最小稳定时间增量一般地是非常小的值，所以大多少问题需要大量的时间增量步。 在求解准静态问题上，显式求解方法已经证明是有价值的，另外ABAQUS/Explicit在求解某些类型的静态问题方面比ABAQUS/Standard更容易。在求解复杂的接触问题时，显式过程相对于隐式过程的一个优势是更加容易。此外，当模型很大时，显式过程比隐式过程需要较少的系统资源。 将显式动态过程应用于准静态问题需要一些特殊的考虑。根据定义，由于一个静态求解是一个长时间的求解过程，所以在其固有的时间尺度上分析模拟常常在计算上是不切合实际的，它将需要大量的小的时间增量。因此，为了获得较经济的解答，必须采取一些方式来加速问题的模拟。但是带来的问题是随着问题的加速，静态平衡的状态卷入了动态平衡的状态，在这里惯性力成为更加起主导作用的力。目标是在保持惯性力的影响不显著的前提下用最短的时间进行模拟。 准静态（Quasi-static）分析也可以在ABAQUS/Standard中进行。当惯性力可以忽略时，在ABAQUS/Standard中的准静态应力分析用来模拟含时间相关材料响应（蠕变、膨胀、粘弹性和双层粘塑性）的线性或非线性问题。关于在ABAQUS/Standard中准静态分析的更多信息，请参阅ABAQUS分析用户手册（ABAQUS Analysis User’s Manual）的第6.2.5节“Quasi-static analysis”。 1. 显式动态问题类比 假设两个载满了乘客的电梯。在缓慢的情况下，门打开后你步入电梯。为了腾出空间，邻近门口的人慢慢地推他身边的人，这些被推的人再去推他身边的人，如此继续下去。这种扰动在电梯中传播，直到靠近墙边的人表示他们无法移动为止。一系列的波在电梯中传播，直到每个人都到达了一个新的平衡位置。如果你稍稍加快速度，你会比前面更用力地推动你身边的人，但是最终每个人都会停留在与缓慢的情况下相同的位置。 在快速情况下，门打开后你以很高的速度冲入电梯，电梯里的人没有时间挪动位置来重新安排他们自己以便容纳你。你将会直接地撞伤在门口的两个人，而其他人则没有受到影响。
abaqus计算应力强度因子
重庆大学 课题：Abaqus计算裂纹应力强度因子 学院： 专业： 学号： 姓名：
一、计算裂纹应力强度因子
问题描述：以无限大平板含有一单边裂纹为例，裂纹长度为a=10mm，平板宽度h=30，弹性模量E=210000Pa，泊松比v=0.33，在远场受双向均布拉应力。 使用Abaqus计算该问题： 1、进入part模块 建立平板part，平板的尺寸相对于裂纹足够大，本例尺寸为50x30 （mm）；使用Partation Face:sketch工具，将part分隔成如图1形式 图1 2、进入property模块 建立弹性材料；截面选择平面问题的solid，homogeneous；赋予截 面。 3、进入Assembly模块 实体的类型（instance type）选择independent。 4、进入mesh模块 划分单元格如图2所示。
图2 5、进入interaction模块 指定裂纹special/creak/assign seam；生成裂纹crack 1， special/crack/create；special/crack/edit，对两个裂纹进行应力奇异的 设置。 6、进入step模块 在initial步之后建立static，general步；在 output/history output requests/create/中创建输出变量。 7、进入load模块 定义位移和荷载边界，如图3所示。
图3 8、进入job模块，提交计算 Mises应力分布见图4，在.dat文件中（图5）查看应力强度因子。 图4
图5 计算解析解： 由公式F=1.12?0.23(a/h)+10.6(a/h)2?21.71(a/h)3+30.38(a/h)4 计算得解析解为k=1001 应力强度因子误差为0.09% 二、误差分析 改变板的长度，其他条件不变 1.当长度L=100时 误差为0.5% 2.当板长L=30
Abaqus中应力应变的理解
在ABAQUS 中对应力的部分理解 1、三维空间中任一点应力有6个分量yz xz xy z y ,,,σσσσσσ，，x ，在ABAQUS 中分别对应S11，S22，S33，S12，S13，S23。 2、一般情况下，通过该点的任意截面上有正应力及其剪应力作用。但有一些特殊截面，在这些截面上仅有正应力作用，而无剪应力作用。称这些无剪应力作用的面为主截面，其上的正应力为主应力，主截面的法线叫主轴，主截面为互相正交。主应力分别以321,,σσσ表示，按代数值排列（有正负号）为321σσσ≥≥。其中321,,σσσ在ABAQUS 中分别对应Max. Principal 、Mid. Principal 、Min. Principal ，这三个量在任何坐标系统下都是不变量。 可利用最大主应力判断一些情况：比如混凝土的开裂，若最大主应力（拉应力）大于混凝土的抗拉强度，则认为混凝土开裂，同时通过显示最大主应力的法线方向，可以大致表示出裂缝的开裂方向等。 利用最小主应力，可以查看实体中残余压应力的大小等。 3、弹塑性材料的屈服准则 3.1、Mises 屈服准则 22 13232 2 212)()()(S σ σσσσ σσ=-+-+- 其中s σ为材料的初始屈服应力。 在三维空间中屈服面为椭圆柱面；在二维空间中屈服面为椭圆。 Mises 等效应力的定义为：(牵扯到张量知识) 其中 S 为偏应力张量，其表达式为 其中为应力， I 为单位矩阵，p 为等效压应力（定义如下）： , 也就是我们常见的 )(31z y x p σσ σ ++= 。 还可以具体表达为： 其中 , , 为偏应力张量（反应塑 性变形形状的变化）。 q 在ABAQUS 中对应 Mises ，它有6个分量（随坐标定义的不同而变化）S11，S22，S33，S12，S13，S23 3.2、Trasca 屈服准则 主应力间的最大差值=2k
ABAQUS定义真实应力和真实应变
ABAQUS 中定义真实应力和真实应变 在ABAQUS 中必须用真实应力和真实应变定义塑性.ABAQUS 需要这些值并对应地在输入文件中解释这些数据。 然而，大多数实验数据常常是用名义应力和名义应变值给出的。这时，必须应用公式将塑性材料的名义应力（变）转为真实应力（变）。 考虑塑性变形的不可压缩性，真实应力与名义应力间的关系为： 00l A lA =， 当前面积与原始面积的关系为： 00l A A l = 将A 的定义代入到真实应力的定义式中，得到： 00 ()nom F F l l A A l l σσ=== 其中0 l l 也可以写为1nom ε+。 这样就给出了真实应力和名义应力、名义应变之间的关系： (1)nom nom σσε=+ 真实应变和名义应变间的关系很少用到，名义应变推导如下： 0001nom l l l l l ε-= =- 上式各加1，然后求自然对数，就得到了二者的关系： ln(1)nom εε=+ ABAQUS 中的*PLASTIC 选项定义了大部分金属的后屈服特性。ABAQUS 用连接给定数据点的一系列直线来逼近材料光滑的应力－应变曲线。可以用任意多的数据点来逼近实际的材料性质;所以，有可能非常逼真地模拟材料的真实性质。在*PLASTIC 选项中的数据将材料的真实屈服应力定义为真实塑性应变的函数。选项的第一个数据定义材料的初始屈服应力，因此，塑性应变值应该为零。 在用来定义塑性性能的材料实验数据中，提供的应变不仅包含材料的塑性应变，而是包括材料的总体应变。所以必须将总体应变分解为弹性和塑性应变分量。弹性应变等于真实应力与杨氏模量的比值，从总体应变中减去弹性应变，就得到了塑性应变，其关系为： /pl t el t E ε εεεσ=-=- 其中pl ε是真实塑性应变，t ε是总体真实应变，el ε是真实弹性应变。
abaqus有限元分析
Abaqus分析报告 （齿轮轴） 名称：Abaqus齿轮轴 姓名： 班级： 学号： 指导教师：
一、简介 所分析齿轮轴来自一种齿轮泵，通过用abaqus软件对齿轮轴进行有限元分析和优化。齿轮轴装配结构图如图1，分析图1中较长的齿轮轴。 图1.齿轮轴装配结构图 二、模型建立与分析 通过part、property、Assembly、step、Load、Mesh、Job等步骤建立齿轮轴模型，并对其进行分析。 1.part 针对该齿轮轴，拟定使用可变型的3D实体单元，挤压成型方式。 2.材料属性 材料为钢材，弹性模量210Gpa，泊松比0.3。
3.截面属性 截面类型定义为solid，homogeneous。 4.组装 组装时选择dependent方式。 5.建立分析步 本例用通用分析中的静态通用分析（Static，General）。 6.施加边界条件与载荷 对于齿轮轴，因为采用静力学分析，考虑到前端盖、轴套约束，而且根据理论，对受力部分和轴径突变的部分进行重点分析。 边界条件：分别在三个轴径突变处采用固定约束，如图2。 载荷：在Abaqus中约束类型为pressure，载荷类型为均布载荷,分别施加到齿轮接触面和键槽面，根据实际平衡情况，两力所产生的绕轴线的力矩方向相反，大小按比例分配。 均布载荷比计算： 矩形键槽数据： 长度：8mm、宽度：5mm、高度：3mm、键槽所在轴半径：7mm 键槽压力面积：S1 = 8x3=24mm2 平均受力半径：R1=6.5mm 齿轮数据：= 齿轮分度圆半径：R2 =14.7mm、压力角：20°、 单个齿轮受力面积：S2 ≈72mm2 通过理论计算分析，S1xR1xP1=S2xR2xP2，其中，P1为键槽均布载荷
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积分则可以完全避免这种复杂的处理过程。 为了计算二维情况下的J积分，ABAQUS定义了围绕裂纹尖端由单元组成的环形的积分域，如下图所示。 图2 ABAQUS中围线的定义 ABAQUS在计算围线积分时，采用的是先计算出围线上面所取的若干个离散点处J积分值，然后乘以每个点对应的加权值后，所有点相加来近似地求解出围线积分，即J积分的值和，进而得到复合裂纹的应力强度因子和。 2两条共线裂纹应力强度因子的算例分析 2.1共线双裂纹在压缩荷载作用下应力强度因子的解析解 有许多学者对含有裂纹的无限大板，裂纹尖端的应力强度因子进行了研究。Zhu Z M[5] 等从理论和实验两个方面都做了详细的研究与探讨。基于前人的研究结果，Zhu Z M 给出了共线裂纹的应力函数及其应力强度因子的基本公式，并就共线双裂纹问题进行了研究，给出了裂纹应力强度因子精确的解析解。 图3压缩载荷作用下的含有共线双裂纹的无限大板 2.2 ABAQUS计算共线裂纹应力强度因子
abaqus后处理中各应力解释个人收集修订版
a b a q u s后处理中各应力 解释个人收集修订版 IBMT standardization office【IBMT5AB-IBMT08-IBMT2C-ZZT18】
ABAQUS中的壳单元S33代表的是壳单元法线方向应力，S11 S22 代表壳单元面内的应力。因为壳单元的使用范围是“沿厚度方向应力为0”，也即沿着法相方向应力为0，且满足几何条件才能使用壳单元，所以所有壳单元的仿真结果应力查看到的S33应力均为0。 S11 S22 S33 实体单元是代表X Y Z三个方向应力，但壳单元不是，另外壳单元只有S12，没有S13，S23。 LE----真应变（或对数应变） LEij---真应变 ... 应变分量； PE---塑性应变分量； PEEQ---等效塑性应变 ABAQUS Field Output Stresses S stress components and invariants 应力分量和变量 SVAVG volume-averaged stress components and invariants (Eulerian only) MISESMAX 最大 Mises 应力 TSHR transverse shear stress(for thick shells)横向剪切应力 CTSHR transverse shear stress in stacked continuum shells 连续堆垛壳横向剪切应力 TRIAX stress triaxiality 应力三轴度 VS stress in the elastic-viscous network 弹粘性网格应力 PS stress in the plastic-viscous
abaqus6.11一个从初学到精通粘弹性的分析的经验积累
问题积累（待续） 1.abaqus如何调整图例的大小，就是云图左上角那个图框，字太小了看不清！！ 直接设置图例的字体大小就可以：工具栏viewport>viewport annotation options>legend(选项卡)>text(选项)>set font(按钮)>size，修改size选项中的数字，就可以修改图例大小了。 2.cohesive element ＡＢＡＱＵＳ 在6.11使用cohesive element，定义cohesive材料属性的时候主要步骤： 1.定义一个材料的名字，比如cohesive，不要去定义任何属性（弹性，弹塑性等等）。 2.打开工具栏model--edit keywords，在inp中手动添加材料的各种属性。 PS: 定义section的时候选cohesive，element control选sweep，element type选cohesive，这些是使用cohesive element的基本步骤。 zero thickness的cohesive section设定abaqus所谓的 zero-thickness，其实就是定义cohesive section的initial thickness=1.0。你可以在定义section的时候定义（specify），也可以用系统默认的thickness（也是1.0），这样有关cohesive element 的计算当中，就有displacement（位移）=strain（应变）*thickness ( 1.0 )=strain的数值。我们知道从1914年Ingless和1921年Griffith提出断裂力学开始，一直到60年代都停留在线弹性断裂力
最新Abaqus中应力应变的理解
在ABAQUS 中对应力的部分理解 1、三维空间中任一点应力有6个分量y z xz xy z y ,,,σσσσσσ，，x ，在ABAQUS 中分别对应S11，S22，S33，S12，S13，S23。 2、一般情况下，通过该点的任意截面上有正应力及其剪应力作用。但有一些特殊截面，在这些截面上仅有正应力作用，而无剪应力作用。称这些无剪应力作用的面为主截面，其上的正应力为主应力，主截面的法线叫主轴，主截面为互相正交。主应力分别以321,,σσσ表示，按代数值排列（有正负号）为321σσσ≥≥。其中321,,σσσ在ABAQUS 中分别对应Max. Principal 、Mid. Principal 、Min. Principal ，这三个量在任何坐标系统下都是不变量。 可利用最大主应力判断一些情况：比如混凝土的开裂，若最大主应力（拉应力）大于混凝土的抗拉强度，则认为混凝土开裂，同时通过显示最大主应力的法线方向，可以大致表示出裂缝的开裂方向等。 利用最小主应力，可以查看实体中残余压应力的大小等。 3、弹塑性材料的屈服准则 3.1、Mises 屈服准则 22132322212)()()(S σσσσσσσ=-+-+- 其中s σ为材料的初始屈服应力。 在三维空间中屈服面为椭圆柱面；在二维空间中屈服面为椭圆。 Mises 等效应力的定义为：(牵扯到张量知识) 其中 S 为偏应力张量，其表达式为 其中为应力， I 为单位矩阵，p 为等效压应力（定义如下）： , 也就是我们常见的)(3 1z y x p σσσ++=。 还可以具体表达为： 其中 , , 为偏应力张量（反应塑性 变形形状的变化）。 q 在ABAQUS 中对应 Mises ，它有6个分量（随坐标定义的不同而变化）S11，S22，S33，S12，S13，S23 3.2、Trasca 屈服准则 主应力间的最大差值=2k
abaqus裂纹模拟心得
abaqus裂纹模拟心得 baqus裂纹模拟心得（Contour Integral不是XFEM） 最近由于项目需要，做了一些裂纹相关的模拟，在此把一些心得体会贴到论坛上与大家分享，如有不当之处，欢迎大家指正！ 本帖主要侧重于介绍裂纹定义过程中各个选项的意义，具体的操作过程论坛里已经有高手做了很好的教程，至于断裂力学理论推荐大家看一下沈成康写的《断裂力学》一书。裂纹的定义和输出需要用到interaction模块和step模块： 一、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（定义裂纹扩展方向）：这里定义的其实是一个虚拟的裂纹扩展方向，定义了这个参考方向后，我们才能通过输出的角度判断裂纹扩展方向，可以通过两种方法： o q vector：输入一个方向，用来作为计算裂纹的扩展方向的参考方向； o normal to crack plane：crack plane表示裂纹的对称面（当裂纹在一个平面内时，可能需要分开定义多个裂纹），这种方法下我们只需定义裂纹面的法线方向，通过（t表示裂纹尖端的切线）, 会在每个节点得出一个q方向（如下图）； o 注意：q的方向对输出的应力强度因子，J积分等都会有影响，一般情况下，q最好在裂纹平面内，且垂直于裂尖线的切线，否则算出的应力强度因子，J积分值等等在不同围线积分中会差别较大。 二、step模块 定义好了裂纹相关参数后，我们需要返回step模块定义输出变量： 步骤：菜单/output/history output requests/create，domain：crack，可以输出的值包括：J-integral，Ct-integral，stress intensity factor，T-stress —J-integral ：用于应变率无关材料的准静态分析过程，包括线弹性，非线性弹性，弹塑性材料（单调加载工况）的静态分析。J-integral的优点是和积分路径无关，从而可以避开尖端塑性区的
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积分值等等在不同围线积分中会差别较大。 二、step模块 定义好了裂纹相关参数后，我们需要返回step模块定义输出变量： 步骤：菜单/output/history output requests/create，domain：crack，可以输出的值包括：J-integral，Ct-integral，stress intensity factor，T-stress —J-integral ：用于应变率无关材料的准静态分析过程，包括线弹性，非线性弹性，弹塑性材料（单调加载工况）的静态分析。J-integral的优点是和积分路径无关，从而可以避开尖端塑性区的影响。 —Ct-integral：用于蠕变分析（一般较少用到） —应力强度因子： （1）只能用于分析线弹性材料，表示裂纹尖端的应力场强度； （2）有三个应力强度因子K1，K2，K3，分别对应于张开型，滑开型和撕开型裂纹的应力强度因子 （3）在输出应力强度因子时也会输出一个J-integral值，因为算法不同，这个值
基于ABAQUS的热应力分析
1.1基于ABAQUS的热应力分析 1.1.1 温度场数据处理 (1)打开INP_Generator.exe，出现如下软件界面： 图1.数据处理软件 (2)点击“浏览”按钮，选择由FLUENT导出的inp文件所在路径，如下图 所示： 图2.路径选择 (3)点击“生成”按钮，则在inp文件所在路径下自动生成包含多个温度场的 ABAQUS输入文件ABAQUSinputfile.inp。 图3.生成包含连续温度场INP文件
1.1.2 复材工装模板热应力分析 (1)打开ABAQUS，导入inp文件后，打开Tools菜单下“Set - Manager”， 如下图所示。检查是否有名为“PID6”的set，若没有则创建一个名为 “PID*”的set，set为模板整体。（“*”为任意数字或字母） 图4.创建SET (2)打开Plug-ins菜单下“CAC Project - Composite Analyse”，弹出如下界面。 在Step1标签中输入用到的材料名称并选择工作路径；在Step2中定义铺 层信息，可通过右键删除或添加行；按照Step3和Step4的提示，使用 ABAQUS/CAE自身功能完成剩余分析工作。 (a)
(b) (c) 图5.定义材料及铺层 (3)进入Load模块，定义垂直于模板表面平面部分的局部坐标系。选择“Tools” 菜单下“Datum”，Type选择“CSYS”Method选择“3Points”，然后默认点击“Continue”按钮。依次在模板表面选择坐标原点、X轴上点和XY面上的点，生成局部坐标。 图6.定义模板局部坐标系 (4)点击“Create Boundary Condition”按钮，弹出边界条件定义对话框。
abaqus6_11常见问题及解决办法
1.abaqus如何调整图例的大小，就是云图左上角那个图框，字太小了 看不清！！ 直接设置图例的字体大小就可以： 工具栏viewport>viewport annotation options>legend(选项卡)>text(选项)>set font(按钮)>size，修改size选项中的数字，就可以修改图例大小了。 2.cohesive element ＡＢＡＱＵＳ 在6.11使用cohesive element，定义cohesive材料属性的时候主要步骤： 1.定义一个材料的名字，比如cohesive，不要去定义任何属性（弹性，弹塑性等等）。 2.打开工具栏model--edit keywords，在inp中手动添加材料的各种属性。 PS: 定义section的时候选cohesive，element control选sweep，element type选cohesive，这些是使用cohesive element的基本步骤。 zero thickness的cohesive section设定abaqus所谓的 zero-thickness，其实就是定义cohesive section的initial thickness=1.0。你可以在定义section的时候定义（specify），也可以用系统默认的thickness（也是1.0），这样有关cohesive element 的计算当中，就有displacement（位移）=strain（应变）*thickness ( 1.0 )=strain的数值。我们知道从1914年Ingless和1921年Griffith提出断裂力学开始，一直到60年代都停留在线弹性断裂力学(LEFM)的层次。后来由於发现在裂纹尖端进入塑性区后用LEFM仍
Abaqus裂纹模拟心得(Contour Integral不是XFEM)
Abaqus裂纹模拟心得（Contour Integral不是XFEM） 最近由于项目需要，做了一些裂纹相关的模拟，在此把一些心得体会贴到论坛上与大家分享，如有不当之处，欢迎大家指正！ 本帖主要侧重于介绍裂纹定义过程中各个选项的意义，具体的操作过程论坛里已经有高手做了很好的教程，至于断裂力学理论推荐大家看一下沈成康写的《断裂力学》一书。 裂纹的定义和输出需要用到interaction模块和step模块： 一、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（定义裂纹扩展方向）：这里定义的其实是一个虚拟的裂纹扩展方向，定义了这个参考方向后，我们才能通过输出的角度判断裂纹扩展方向，可以通过两种方法： o q vector：输入一个方向，用来作为计算裂纹的扩展方向的参考方向； o normal to crack plane：crack plane表示裂纹的对称面（当裂纹在一个平面内时，可能需要分开定义多个裂纹），这种方法下我们只需定义裂纹面的法线方向，通过（t表示裂纹尖端的切线）, 会在每个节点得出一个q方向（如下图）； o 注意：q的方向对输出的应力强度因子，J积分等都会有影响，一般情况下，q最好在裂纹平面内，且垂直于裂尖线的切线，否则算出的应力强度因子，J积分值等等在不同围线积分中会差别较大。
Abaqus中应力应变的理解
在ABAQUS 中对应力得部分理解 1、三维空间中任一点应力有6个分量y z xz xy z y ,,,σσσσσσ，，x ，在ABAQUS 中分别对应S11，S22，S33，S12，S13，S23。 2、一般情况下，通过该点得任意截面上有正应力及其剪应力作用。但有一些特殊截面，在这些截面上仅有正应力作用，而无剪应力作用。称这些无剪应力作用得面为主截面，其上得正应力为主应力，主截面得法线叫主轴，主截面为互相正交。主应力分别以321,,σσσ表示，按代数值排列（有正负号）为321σσσ≥≥。其中321,,σσσ在ABAQUS 中分别对应Max 、 Principal 、Mid 、 Principal 、Min 、 Principal ，这三个量在任何坐标系统下都就是不变量。 可利用最大主应力判断一些情况：比如混凝土得开裂，若最大主应力（拉应力）大于混凝土得抗拉强度，则认为混凝土开裂，同时通过显示最大主应力得法线方向，可以大致表示出裂缝得开裂方向等。 利用最小主应力，可以查瞧实体中残余压应力得大小等。 3、弹塑性材料得屈服准则 3、1、Mises 屈服准则 22132322212)()()(S σσσσσσσ=-+-+- 其中s σ为材料得初始屈服应力。 在三维空间中屈服面为椭圆柱面；在二维空间中屈服面为椭圆。 Mises 等效应力得定义为：(牵扯到张量知识) 其中 S 为偏应力张量，其表达式为 其中为应力， I 为单位矩阵，p 为等效压应力（定义如下）： , 也就就是我们常见得)(31z y x p σσσ++=。 还可以具体表达为： 其中 , , 为偏应力张量（反应塑 性变形形状得变化）。 q 在ABAQUS 中对应 Mises ，它有6个分量（随坐标定义得不同而变化）S11，S22，S33，S12，S13，S23 3、2、Trasca 屈服准则 主应力间得最大差值=2k 若明确了321σσσ≥≥，则有k =-)(2 131σσ，若不明确就需要分别两两求差值，瞧哪个最大。 ABAQUS 中得Trasca 等效应力就就是“主应力间得最大差值” 3、3 ABAQUS 中得Pressure----等效压应力 即为上面提到得p ：, 也就就是我们常见得)(3 1z y x p σσσ++=。
含圆孔和裂纹板应力强度因子分析
《断裂力学》 大作业 题目：含圆孔和裂纹板应力强度因子分析. 姓名： 学号： 专业： 授课教师： ^
一、问题描述 含多裂纹矩形板受垂直方向拉伸载荷作用，如图 1 所示，计算中心裂纹尖端的应力强度因子KⅠ和KⅡ，并讨论其随即和参数L、h、a、D、 等的变化规律，写一篇分析报告。 图1. 含三条裂纹矩形板受垂直拉伸载荷作用 要求 （1）报告中计算所用到的分析方法和模型应阐述清楚，并写出必要的计算公式。 （2）绘制应力强度因子随几何参数的变化曲线。 （3）列出必要的参考文献 二、理论分析
— 在线弹性断裂力学中，I型裂纹尖端的应力场为： (1sin sin) 222 (1sin sin) 222 cos cos 222 3 3 3 x y xy σ σ τ θθθ θθθ θθθ ? =- ? ? ? =+ ? ? ? = ? ? I型裂纹尖端的位移场为： 1)cos(1cos) 22 1)sin sin 22 3 3 u v κ κ θθ θθ ? =-- ? ? ? ?=+ ?? 其中： 34 3 1 ν κν ν - ? ? =?- ?+ ? 平面应变 平面应力 同理，对II型裂纹尖端的应力场： (2cos cos) 222 cos sin cos 222 (1sin sin) 222 3 3 3 x y xy σ σ τ θθθ θθθ θθθ ? =+ ? ? ? = ? ? ? =- ? ? 显然，位移场和应力场均可以表示成应力强度因子的形式。通过对裂纹尖端的应力应变场分析来求解对应的应力强度因子，便是传统有限元求解应力强度因子的原理。而对于I、II复合型裂纹尖端的应力强度因子，可通过它们的叠加获得。 确定应力强度因子的方法有3大类：解析法、数值解法和实验方法。解析法只能计算简单问题，大多数问题需要采用数值解法，当前工程中广泛采用的数值解法是有限单元法。随着有限元法的发展，有限元
ABAQUS中的断裂力学及裂纹分析总结
ABAQUS中的断裂力学及裂纹分析总结（转自simwe） （1） 做裂纹ABAQUS有几种常见方法。最简单的是用debond命令, 定义 *FRACTURE CRITERION, TYPE=XXX, 参数。。。 ** *DEBOND, SLAVE=XXX, MASTER=XXX, time increment=XX 0,1, …… ...... time,0 要想看到开裂特别注意需要在指定的开裂路径上定义一个*Nset，然后在 *INITIAL CONDITIONS, TYPE=CONTACT中定义 master, slave, 及指定的Nset 这种方法用途其实较为有限。 （2） 另一种方法，在interaction模块，special, 定义crack seam, 网格最好细化，用collapse element模拟singularity. 这种方法可以计算J积分，应力强度因子等常用的断裂力学参数. 裂尖及奇异性定义: 在interaction-special,先定义crack, 定义好裂尖及方向, 然后在singularity选择：midside node parameter: 输入0.25, 然后选Collapsed element side, duplicate nodes，8节点单元对应(1/r)+(1/r^1/2)奇异性。 这里midside node parameter选0.25对应裂尖collapse成1/4节点单元。如果midside nodes 不移动到1/4处, 则对应(1/r)奇异性, 适合perfect plasticity的情况. 网格划分: 裂尖网格划分有一些技巧需要注意，partition后先处理最外面的正方形，先在对角线和边上
采用abaqus的cae进行力学问题的分析
采用abaqus的cae进行力学问题的分析，其对模型的处理存在很多的技巧，对abaqus的一些分析技巧进行一些概述，希望对大家有所帮助。 abaqus的多图层绘图 abaqus的cae默认一个视区仅仅绘出一个图形，譬如contor图，变形图，x-y曲线图等，其实在abaqus里面存在一个类似于origin里面的图层的概念，对于每个当前视区里面的图形都可以建立一个图层，并且可以将多个图层合并在一个图形里面，称之为Overlay Plot，譬如你可以在同一副图中，左边绘出contor图，右边绘出x-y图等等，并且在abaqus里面的操作也是很简单的。 1.首先进入可视化模块，当然要先打开你的模型数据文件(.odb) 2.第一步要先创建好你的图形，譬如变形图等等 3.进入view里面的overlay plot，点击creat，创建一个图层，现在在viewport layer里出现了你创建的图层了 4.注意你创建的图层，可以看到在visible 下面有个选择的标记，表示在视区里面你的图层是否可见，和autocad里面是一样，取消则不可见current表示是否是当前图层，有些操作只能对当前图层操作有效，同cad name是你建立图层的名称，其他的属性值和你的模型数据库及图形的类型有关，一般不能改动的。 5.重复2-4步就可以创建多个图层了 6.创建好之后就可以选择plot/apply，则在视区显示出所有的可见的图层 子结构的概述 1.什么是子结构 子结构也叫超单元的(两者还是有点区别的，文后会谈到)，子结构并不是abaqus里面的新东东，而是有限元里面的一个概念，所谓子结构就是将一组单元组合为一个单元(称为超单元)，注意是一个单元，这个单元和你用的其他任何一种类型的单元一样使用。 2.为什么要用子结构 使用子结构并不是为了好玩，凡是建过大型有限元模型的兄弟们都可能碰到过计算一个问题要花几个小时，一两天甚至由于单元太多无法求解的情况，子结构正是针对这类问题的一种解决方法，所以子结构肯定是对一个大型的有限元模型的，譬如在求解非线性问题的时候，因为对于一个非线性问题，系统往往经过多次迭代，每次这个系统的刚度矩阵都会被重新计算，而一般来说一个大型问题往往有很大一部分的变形是很小的，把这部分作为一个子结构，其刚度矩阵仅要计算一次，大大节约了计算时间。 3.哪些情况可以使用子结构 前面提到的非线性问题，包括了很小变形的或者线弹性部分可以使用子结构，特别是当模型中有很多相同的部分时，提到的最多的一个例子就是桌子的四条腿，四条腿作为子结构(因为基本时