基于ANSYS的三维贯穿裂纹的断裂参数计算
ANSYSworkbench裂纹分析
基于ANSYS Workbench的表面裂纹计算By Yan Fei本教程使用ANSYS Workbench17.0 进行试件表面裂纹的分析，求应力强度因子。
需要提前说明的是，本案例没有工程背景，仅为说明裂纹相的计算方法，因此参数取值比较随意，大量设置都采用了默认值。
1.背景知识传统的强度设计思想把材料视为无缺陷的均匀连续体，而实际工程构件中存在多种缺陷，断裂力学是从20实际50年代末期发展起来的一门弥补了传统强度设计思想严重不足的新的学科，是专门研究含缺陷或裂纹的物体在外界条件作用下构件的强度、裂纹扩展趋势以及疲劳寿命的科学。
断裂力学是从构件内部具有初始缺陷这一实际情况出发，研究在外部荷载下的裂纹扩展规律，从而提出带裂纹构件的安全设计准则。
a 张开型裂纹b 滑开型裂纹c 撕开型裂纹图 1 裂纹的分类使用弹性力学方法可以求得，在裂纹尖端处的应力的解析解为无穷大，此时应力值已经失去意义，一般采用应力强度因子作为判断结构是否安全的指标。
目前的断裂力学研究主要集中在I型裂纹的开裂，数值计算工具也多集中在I型裂纹的计算上，因此以I型裂纹为例。
图2 裂纹尖端坐标系含有裂纹的无限大平板的I 型裂纹尖端附近的应力为：)(23cos 2sin 223sin 2sin 12cos223sin 2sin 12cos20ⅠⅠⅠr O r K r K r K xyy x +=+=-=其中，K Ⅰ叫Ⅰ型裂纹的应力强度因子。
2.ANSYS Workbench 裂纹分析2.1.分析模型的建立1 建立一个静力分析步，材料使用默认，需要说明的是，现有计算技术下，断裂力学计算一般都采用线弹性材料，考虑到断裂中塑性区一般都不大，线弹性的假设还是可以接受的。
图3 分析步设置2 建立几何模型，本案例使用spaceclaim 建立几何模型。
图4 试件平面图图5 试件立体图3 分网格，必须采用四面体网格。
本文划分单元特征尺寸1mm。
图 6 网格设置图7 分网效果4 划分网格完成以后，首先进行一次静力计算，确保所有设置正确，对ANSYS Workbench比较熟悉的同学可以省略这一步，静力计算时，试件的两个端面一个约束位移，另一个加1000N的力，方向沿试件轴向，使试件受拉。
基于ANSYS的断裂参数的计算
基于ANSYS的断裂参数的计算1 引言断裂事故在重型机械中是比较常见的，我国每年因断裂造成的损失十分巨大。
一方面，由于传统的设计是以完整构件的静强度和疲劳强度为依据，并给以较大的安全系数，但是含裂纹在役设备还是常有断裂事故发生。
另一方面，对于一些关键设备，缺乏对不完整构件剩余强度的估算，让其提前退役，从而造成了不必要的浪费。
因此，有必要对含裂纹构件的断裂参量进行评定，如应力强度因了和J积分。
确定应力强度因了的方法较多，典型的有解析法、边界配位法、有限单元法等。
对于工程上常见的受复杂载荷并包含不规则裂纹的构件，数值模拟分析是解决这些复杂问题的最有效方法。
本文以某一锻件中取出的一维断裂试样为计算模型，介绍了利用有限元软件ANSYS计算应力强度因子。
2 断裂参量数值模拟的理论基础对于线弹性材料裂纹尖端的应力场和应变场可以表述为：(1)其中K是应力强度因子，r和θ是极坐标参量，可参见图1，(1)式可以应用到三个断裂模型的任意一种。
图1 裂纹尖端的极坐标系(2)应力强度因子和能量释放率的关系：G=K/E"(3)其中：G为能量释放率。
平面应变：E"=E/(1-v2)平面应力：E=E"3 求解断裂力学问题断裂分析包括应力分析和计算断裂力学的参数。
应力分析是标准的ANSYS 线弹性或非线性弹性问题分析。
因为在裂纹尖端存在高的应力梯度，所以包含裂纹的有限元模型要特别注意存在裂纹的区域。
如图2所示，图中给出了二维和三维裂纹的术语和表示方法。
图2 二维和三维裂纹的结构示意图3.1 裂纹尖端区域的建模裂纹尖端的应力和变形场通常具有很高的梯度值。
场值得精确度取决于材料，几何和其他因素。
为了捕获到迅速变化的应力和变形场，在裂纹尖端区域需要网格细化。
对于线弹性问题，裂纹尖端附近的位移场与成正比，其中r是到裂纹尖端的距离。
在裂纹尖端应力和应变是奇异的，并且随1/变化而变化。
为了产生裂纹尖端应力和应变的奇异性，裂纹尖端的划分网格应该具有以下特征：·裂纹面一定要是一致的。
基于ANSYS有限元软件裂纹扩展模拟
基于ANS Y S有限元软件裂纹扩展模拟刘 莎3 张 芳(武汉铁路局武昌东站技术科)(十堰东风商用车研发中心) 摘 要 从能量释放率准则出发,用AN SYS软件作为平台,进行二次开发来模拟二维复合加载下裂纹的扩展。
裂纹扩展路径的模拟是模拟裂纹扩展中的难点。
重点描述了模拟裂纹扩展路径。
关键词 裂纹 能量释放率 裂纹扩展 Paris公式0 前言 裂纹扩展有限元模拟研究涉及三个问题:理论基础、扩展控制参量及模拟方法。
理论基础直接影响有限元方程构成和具体实施的难易程度,应用全量理论还是增量理论;采用非线性弹性假设还是考虑扩展过程中能量耗散的真实弹塑性本构关系、屈服条件、小变形、有限变形或大变形理论等等。
就目前看来,研究主要以非线性弹性及小变形理论为主,且大多数采用V on M ises屈服准则,包括能量耗散在内的真实弹塑性及大变形理论的有限元研究者也有,但研究不很系统。
裂纹扩展控制参数与断裂理论发展及裂纹扩展测量技术有关。
扩展控制参数研究是目前弹塑性断裂问题有限元热门课题之一。
扩展模拟控制参数主要有下述几种:J R-∃a控制,J R -CTOA联合控制,载荷控制P-∃a及载荷线位移控制LLD-∃a,能量释放准则控制G-∃a。
在J R-∃a控制的实施过程中,J R阻力曲线作为材料特性,并假设与样本几何性及加载过程无关。
如上所述,此假设的合理性尚存在某些疑问,已有证据表明,当裂纹扩展量增大时, J R阻力曲线的样本几何依赖性明显增大。
尽管如此,在目前裂纹扩展有限元研究中,J R-∃a 控制仍是应用最广泛的方法之一,包括在裂纹扩展量大的情况下,其中原因之一是J R阻力曲线属于远场J,而远场的有限元实施具有相当的数值稳定性。
裂纹小量扩展后,CTOA基本保持常数的特性使J R-CTOA联合作为扩展准则具有潜在的发展前景,因为当J R-CTOA 联合使用时,可避开临界CTOA测量这个难点,即有限元实施时的小量扩展阶段应用J R-∃a控制,同时连续计算CTOA,当CTOA为常数开始点时,也几乎是J控制失效点,随后的裂纹扩展则用CTOA代替J R作为控制参数。
基于ansys的胶粘结构界面开裂有限元计算
基于ansys的胶粘结构界面开裂有限元计算胶粘结构在机械工程、航天航空和汽车工业等领域中起着关键作用。
然而，界面的开裂问题一直是一个令人头痛的难题。
在本文中，我们将探讨如何使用Ansys有限元软件来进行基于胶粘结构界面开裂的计算。
1. 胶粘结构界面开裂的挑战胶粘结构界面开裂问题是由于应力集中引起的，在许多应用中都是影响结构完整性的重要因素。
胶粘接过程中的温度梯度和载荷变化会导致界面开裂。
因此，准确地计算胶粘结构界面的应力分布和开裂扩展行为对于保证结构的可靠性至关重要。
2. 使用Ansys进行胶粘结构界面开裂的有限元计算Ansys是一种广泛应用于工程领域的有限元软件，它提供了多种功能和模块，可以用于各种结构分析。
在胶粘结构界面的开裂问题中，我们可以利用Ansys的有限元模块进行模拟和计算。
首先，我们需要建立一个精确模型来描述胶粘结构界面的几何形状和材料性质。
Ansys提供了几何建模工具，可以根据实际情况创建并优化模型。
在模型建立完成后，我们需要为界面设置适当的材料属性和边界条件，以反映真实的工程情况。
接下来，我们可以利用Ansys的有限元求解器来进行界面开裂的模拟计算。
有限元方法基于将结构分割成许多小单元，利用数学方法求解每个单元上的力学方程，并找到整个结构的应力分布。
在胶粘结构界面开裂的有限元计算中，我们可以使用线性弹性有限元或非线性有限元方法来模拟接触压力和开裂行为。
为了更准确地模拟胶粘结构界面开裂的过程，我们还可以考虑温度和湿度等因素的影响。
Ansys提供了耦合热和湿度分析的功能模块，可以帮助我们更好地理解胶粘结构界面的开裂机制。
3. 结果分析和验证在进行胶粘结构界面开裂的有限元计算后，我们可以通过Ansys提供的后处理工具来分析模拟结果。
这些工具可以帮助我们可视化应力分布、应变分布和开裂扩展情况，并对计算结果进行验证。
验证是一个重要的步骤，可以与实验结果比较，并评估模拟计算的准确性和可靠性。
断裂力学参量[整理版]
![断裂力学参量[整理版]](https://img.taocdn.com/s3/m/b8a2a748814d2b160b4e767f5acfa1c7aa0082f3.png)
ANSYS求解断裂力学参量的理论方法工程上,线弹性断裂力学中常用应力强度因子K、J积分、G能量释放率这三个参量来描述裂纹场。
ANSYS软件能较好地计算裂纹周围区域的应力分布，并能计算裂纹的应力强度因子K、J积分以及能量释放率G等，其特点是简单、经济、精度高。
下面主要介绍在ANSYS中如何求解应力强度因子K和J积分。
（1）求解应力强度因子ANSYS软件中提供了所谓的“位移外推”法(displacement extrapolation) 来计算应力强度因子[5]。
在线弹性范围内,对于三维裂纹,裂纹尖端的局部位移场与应力强度因子的关系为[6]：)2)22IIIIIIKu kGKv kGKwG⎧=+⎪⎪⎪⎪=+⎨⎪⎪⎪=⎪⎩式中: u、v、w—如图2.5所示裂纹尖端局部直角坐标系下裂纹前端位移；r—如图2.5所示裂纹尖端局部柱坐标系下坐标；G—材料剪切模量；K I、K II、K III—应力强度因子；v—为泊松比；34()3()1vk vv-⎧⎪=⎨-⎪+⎩平面应变或轴对称平面应力当利用裂纹尖端节点的位移进行计算时，应力强度因子和裂纹面节点的位移差存在下列关系：IIIIIIKKK⎧=⎪⎪⎪⎪=⎨⎪⎪⎪=⎪⎩三维裂纹的局部坐标在使用有限元法进行应力强度因子计算时,由于常规单元在裂纹尖端存在奇异性,为使计算准确,必须在裂纹尖端使用细小的单元;如果使用奇异元,即使用二次三角(或五面体)单元,并将靠近裂纹尖端的中间节点置于1/4处,则位于沿裂纹尖端的单元边上的应力和应变与1/消除了奇异性，也就是说，可以用相对比较稀疏的单元得到精度较高的结果。
（2）求解J积分J积分定义为一个围绕裂尖的线积分(二维) 或一个围绕裂纹前沿的面积分。
它用计算裂纹尖端的奇异应力和应变,与积分路径无关。
为了避开裂纹尖点的奇异性,取得较好的精度,积分路径一般取得离裂纹尖点较远。
J积分形式如图2.6所示，其表达式如下：()yxx yuuJ Wdy t t dsx yΓΓ∂∂=-+∂∂⎰⎰式中：W—应变能密度(单位体积应变能)；Г—围绕裂纹尖点任意路径；xt—X 方向的作用向量,x x xy yt nσσ=+;yt—Y方向的作用向量,y y xy xt nσσ=+;n—积分路径的外法向向量；s —积分路径距离；围绕裂纹尖端的任意一条J 积分路径在ANSYS 中，为了计算位移向量的偏导数x u x ∂∂与y u y ∂∂，将积分路径向x 正负方向分别移动Δx/2，并求出路径Γ+Δx/2上各点的位移u x1和u y 1以及路径Γ-Δx/2上各点的u x 1和u y 1，则：2121()()x x x y y y u x u u xu y u u y∂∂=-∆⎧⎪⎨∂∂=-∆⎪⎩ ANSYS 具有强大的后处理功能,利用此功能,在求解后可以通过ANSYS 通用后处理器中的单元列表功能,很方便地把各变量映射到自定义的路径中去。
采用ANSYS仿真模拟软件建立三维混凝土试件实体裂纹扩展的模拟
MSHKEY,0
!*
CM,_Y,VOLU
VSEL, , , , 3
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
!施加下端约束
FLST,2,1,5,ORDE,1
FITEM,2,18
!通水管道施压
FLST,2,2,5,ORDE,2
FITEM,2,12
FITEM,2,-13
/GO
!*
SFA,P51X,1,PRES,-2000000
!裂纹面施压
FLST,2,4,5,ORDE,2
FITEM,2,14
FITEM,2,-17
/GO
!*
SFA,P51X,1,PRES,2000000
Fini !退出前处理器
wpro,,90.000000, !旋转工作平面
CSWPLA,100,1,1,1, !在工作平面位置建立局部坐标100，类型为柱坐标
FLST,3,1,6,ORDE,1
FITEM,3,1
VGEN, ,P51X, , , ,45, , , ,1 !旋转长方体
wpro,,-90.000000, !旋转回原工作平面
!(4)定义材料参数
MP,EX,1,1.668E10 !弹性模量
MP,PRXY,1,0.3 !泊松比
!(5)建立剖面几何模型
BLOCK,-0.015,0.015,-0.025,0.025,-0.0005,0.0005, !建立一个长方体
WPSTYLE,,,,,,,,1
Ansys 断裂力学理论
第四章断裂力学文献来源：/document/200707/article796_2.htm4.1 断裂力学的定义在许多结构和零部件中存在的裂纹和缺陷，有时会导致灾难性的后果。
断裂力学在工程领域的应用就是要解决裂纹和缺陷的扩展问题。
断裂力学是研究载荷作用下结构中的裂纹是怎样扩展的，并对有关的裂纹扩展和断裂失效用实验的结果进行预测。
它是通过计算裂纹区域和破坏结构的断裂参数来预测的，如应力强度因子，它能估算裂纹扩展速率。
一般情况下，裂纹的扩展是随着作用在构件上的循环载荷次数而增加的。
如飞机机舱中的裂纹扩展，它与机舱加压及减压有关。
此外，环境条件，如温度、或大范围的辐射都能影响材料的断裂特性。
典型的断裂参数有：与三种基本断裂模型相关的应力强度因子(K I,K II,K III)(见图4-1)；J积分，它定义为与积分路径无关的线积分，用于度量裂纹尖端附近奇异应力与应变的强度；能量释放率(G)，它反映裂纹张开或闭合时功的大小；注意--在本节大部分的图形中裂纹的宽度被放大了许多倍。
图4-1 裂缝的三种基本模型4.2 断裂力学的求解求解断裂力学问题的步骤为：先进行线弹性分析或弹塑性静力分析，然后用特殊的后处理命令、或宏命令计算所需的断裂参数。
本章我们集中讨论下列两个主要的处理过程。
裂纹区域的模拟；计算断裂参数。
4.2.1 裂纹区域的模拟在断裂模型中最重要的区域，是围绕裂纹边缘的部位。
裂纹的边缘，在2D模型中称为裂纹尖端，在3D模型中称为裂纹前缘。
如图4-2所示。
图4-2 裂纹尖端和裂纹前缘在线弹性问题中，在裂纹尖端附近(或裂纹前缘)某点的位移随而变化，γ是裂纹尖端到该点的距离，裂纹尖端处的应力与应变是奇异的，随1/变化。
为选取应变奇异点，相应的裂纹面需与它一致，围绕裂纹顶点的有限元单元应该是二次奇异单元，其中节点放到1/4边处。
图4-3表示2-D和3-D模型的奇异单元。
图4-3 2-D和3-D模型的奇异单元4.2.1.1 2-D断裂模型对2D断裂模型推荐采用PLANE2单元，其为六节点三角形单元。
在ANSYS中计算裂缝应力强度因子的技巧
在ANSYS中计算裂缝应力强度因子的技巧在ANSYS中计算裂缝应力强度因子的技巧裂缝应力强度因子用ANSYS中怎么求呀。
另外，建模时，裂纹应该怎么处理呀，难道只有画出一条线吗？首先说一下裂纹怎么画，其实裂纹很简单啊。
只要画出裂纹的上下表面（线）就可以了，即使是两个面（线）重合也一定要是两个面（线）；如果考虑道对称模型就更好办了，裂纹尖点左面用一个面（线），右边用另外一个面（线），加上对称边界约束。
再说一下裂尖点附近网格的划分。
ansys提供了一个kscon的命令，主要是使得crack tip的第一层单元变成奇异单元，用来模拟断裂奇异性(singularity)。
当然这个步骤不是必须的，有的人说起用ansys算强度因子的时候就一定要用奇异单元，其实是误区（原因下面解释）好了，回到强度因子的计算。
其实只要学过一些断裂力学都知道，K的求法很多。
就拿Mode I的KI来说吧，Ansys自己提供了一个办法（displacement extrapolation），中文可能翻译作“位移外推”法，其实就是根据解析解的位移公式来对计算数据进行fitting的。
分3步走，如果你已经算完了：第一步，先定义一个crack-tip的局部坐标系，这是ansys帮助文件中说的，其实如果你的裂纹尖端就是整体坐标原点的话，而且你的x-axis就顺着裂纹，就没有什么必要了。
第二步，定义一个始于crack-tip的path,什么什么？path怎么定义？？看看帮助吧，在索引里面查找fracture mechanics，找到怎么计算断裂强度因子。
（my god,我这3步全是在copy 帮助中的东东啊）。
第三步，Nodal Calcs>Stress Int Factr ，别忘了，这是在后处理postproc中啊。
办法是好，可是对于裂纹尖端的单元网格依赖性很大，所以用kscon制造尖端奇异单元很重要。
curtain的经验是path路径取的越靠近cracktip得到的强度因子就越大，所以单元最好是越fine越好啊。
用ANSYS作裂纹走向预测的计算技巧
用ANSYS作裂纹走向预测的计算技巧ANSYS是一个广泛应用于工程领域的有限元分析软件，用于模拟和解决各种工程问题。
在裂纹走向预测方面，ANSYS提供了多种功能和技巧。
本篇文章将介绍ANSYS在裂纹走向预测方面的计算技巧，并提供一些实用的方法和建议。
以下是一些值得关注的关键步骤和技巧：1.建立准确的模型：在进行裂纹走向预测之前，需要建立一个符合实际情况的准确模型。
模型的准确性对于预测结果的准确性至关重要。
在建模过程中，需要考虑材料的性质、裂纹的大小和方向以及与裂纹配合的部件的几何形状。
2.材料参数的输入：ANSYS提供了材料数据库，可以选择标准材料参数。
然而，在一些情况下，需要自定义材料参数。
这涉及到材料的宏观和微观力学性质。
这些材料参数包括弹性模量、屈服强度、破坏韧性等。
正确输入材料参数对于准确预测裂纹走向至关重要。
3.边界条件的设置：边界条件对于裂纹行为的模拟非常关键。
在模型中正确设置边界条件将能够准确预测裂纹的行为。
对于裂纹走向预测，需要考虑材料的加载状态和应力分布。
要模拟真实情况下材料的力学行为，可以设置边界条件来模拟真实的受力情况。
4.裂纹尺寸的输入：在模拟裂纹行为时，需要定义裂纹的尺寸。
ANSYS提供了多种定义裂纹尺寸的方法，包括手动输入和自动生成。
在裂纹走向预测中，可以通过输入不同的裂纹尺寸来模拟不同的裂纹形态，然后预测不同的裂纹走向。
5.工程应力的加载：工程应力加载是模拟实际工程问题的关键步骤之一、通过在模型中应用工程应力，可以模拟裂纹行为的响应。
可以在ANSYS中使用加载边界条件来模拟不同的加载条件，例如拉伸、压缩或弯曲。
6. 材料损伤准则的选择：裂纹走向预测中，需要选择适当的材料损伤准则。
材料损伤准则用于预测裂纹扩展方向和速率。
ANSYS提供了多种材料损伤准则，如J-Integral、CTOD等。
选择适当的材料损伤准则可以提高预测结果的准确性。
7.结果分析和后处理：在模拟完成后，需要对结果进行分析和后处理。
基于ANSYS的穿透裂纹三维板数值分析
（ ３ ）
…
板中心裂纹
式中： 厂—— 围绕裂纹尖 端 的一条任 意 反时 针 回路 ， 起端始 于裂 卜 弹塑性条件下 ， 回路 ，上任 意点 （ Ｙ ，）的应变 能密
图 １ 函数 ｔ） ０ 取值 图 表 １ 中 心裂 纹 计 算 结 果
・
５ ６・
第３ ８卷 第 ２ ９期
２ ２ 月 １ 年１ ０ ０
山 西 建 筑
裂纹 的力学特性 。板表面均匀压力荷载及约束如图 ８ 图 ９所示 。 阶梯荷载分别计算 。模拟结果与解 析解对 比， ， 得到三 种不 同位 置
３ 计 算方 法及 理论依 据
３ １ Ｊ积 分 ．
３ １０ ． １ ４６ ５ ．６ ６２ ０ ．２ ７ ７５ ．７ ９ ３０ ．３
误 差／ ％
１８ ．３ １８ ．１ １８ ．２ １８ ．２ １８ ．２
Ｋ ＡＪ Ｃｌ Ｃ
３１９ ．２ ４ ６４ ．９ ６２８ ．５ ７８３ ．２ ９３８ ． ８
误差／ ％
图０三 不 位 裂 情 的（）值 线 １ 种 同置 纹 况 ， 取 曲。 为 ÷
［ ６
篓 研 院 力 度 子 册Ｍ 北 ：学 版 究 ・ 强 因 手 ［ ・京科 出 应 ］
Ｏｎ ｎ ｍ ｅ ｉ ｎ ｌ ｓｓ ｏ ｅ ｅ ｒ ｔ ｄ ｃ ａ ｋ ｔ ｒ ｅ ｄ ｍ ｅ ｓｏ ｌ ｔ ａ ｅ ｎ ＡＮＳ ｕ ｒｃ ａ ａ ｙ ｉ ｆｐ ｎ ｔ ａ ｅ ｒ ｃ ｈ ｅ ． ｉ ｎ ｉ ｎ ｐ ａ ｅ ｂ ｓ ｄ ｏ ＹＳ
， ，
１ ｉａ ．２ｒ ２
爹哼． 卧 参考文献 ： ｘ ． ［ ］ 康颖安． １ 断裂力学的发展 与研 究现 状 ［ ］ 湖 南工程 学 院学 Ｊ．
基于ANSYS的三维应力强度因子的计算
基于ANSYS的三维应力强度因子的计算摘要:本文在总结断裂力学各种行为研究方法的基础上，采用有限单元法思想，利用ANSYS软件建立裂纹体有限元模型。
通过计算得出I型裂纹尖端的应力强度因子，其计算结果与理论值吻合良好。
表明模型的选取和网格的划分是合理的，具有可靠的精度，其结果完全可以指导工程设计。
关键词：断裂力学；有限单元法；ANSYS；应力强度因子1 前言随着现代生产技术的高速发展，新材料、新工艺在航空、航天、压力容器、核反应堆、机械、土木等领域得到了广泛应用，结构在高速、高温、高压等环境中使用时，按照传统强度理论设计的结构在应用中却出现大量的断裂事故。
目前关于断裂力学的研究，理论上只能求解较简单的模型或做出较强的假设条件；通过实验探求其规律性的成本较高、周期较长。
因此对于断裂力学问题，尤其是三维裂纹问题目前大多借助数值方法进行研究。
裂纹问题严格来说都是三维问题，并且工程中最后发生事故的裂纹问题的物体(如机械构件，土建结构)总是有限弹性的。
因此对有限弹性裂纹问题进行三维分析在实际工程上有重要的现实意义[1]。
迄今在平面断裂力学中已形成了一整套相当成熟的计算方法，但在三维断裂问题中，尤其在表面裂纹方面，还有很多问题有待进一步探讨，本文正是在这方面进行了探索和研究。
通常板表面裂纹应力强度因子可以统一表示为：(1)只是不同的解给修正因子F赋予不同的表达式。
由Newman的研究成果可以看出，他主要考虑了拉伸或弯曲载荷下半椭圆表面裂纹的各种裂纹形状（a/c）和板宽、板厚对裂纹前缘应力强度因子的影响，而没有考虑板长的影响，也就是没有考虑平行边界对应力强度因子的影响。
2 裂纹类型与有限元法2.1 裂纹类型在断裂力学中，按裂纹受力情况和裂纹面的相对位移方向将裂纹分为三种基本类型[2]：即张开型（I型），滑移型（II型），撕裂型（III型），如图1。
两种或三种基本类型的组合称之为复合型裂纹问题。
图1中所示三种基本类型的裂纹模型的受力特点如下：I型裂纹受垂直于裂纹面的拉应力作用；II型裂纹受平行于裂纹面而垂直于裂纹前缘的剪应力作用，又称为面内剪切型；III型裂纹受既平行于裂纹面又平行于裂纹前缘的剪应力作用，对应于反平面剪切，又称为面外剪切或纵向剪切。
断裂参量在ANSYS中的计算分析
$! 结论
"# $ ! 通过以上的计算分析， 说明断裂参量在 %&’(’ 中的计算是可行的， 为复杂裂纹在复杂载荷下的断裂判 据的计算提供了有效的方法 "# % ! 宏观裂纹的倾角对于应力强度因子和 # 积分值都有一定程度的影响， 但是对应力强度因子的影响大一 点， 对于 # 积分的影响不太明显。 "# & ! 宏观裂纹的长度在构件的失稳破坏过程中起着主导作用， 对于应力强度因子和 # 积分的影响都比较显 著% "# " ! 应力强度因子和 # 积分是相关联的， 它们之间有着一定的内在关系。 参考文献：
$# 断裂参量数值计算的理论基础
56 5 # 应力强度因子数值计算的理论基础 在线弹性断裂力学中， 由于裂纹尖端的应力场的强弱程度主要由应力强度因子 " 这个参量来描述， 故 通过它可以建立 " 9 " #$ 的断裂准则， 来解决工程实际的断裂问题。 " 的大小与外载的性质、 裂纹及裂纹弹性体几何形状等因素有关的一个量， 写成通式是 " % & ! %（ ’ "， #）% & !， "， # 式中 ! 是与裂纹有关的几何形状因子， " 是裂纹承受的应力， # 是裂纹长度。 56 7 # ! 积分数值计算的理论基础 ! 积分是断裂力学中的一个重要概念， 其数值是一个与积分回路无关的常数， 即具有守恒性， 它也反映了 裂纹尖端的某种力学特性或应力应变场强度， 同时在分析中可能避开裂纹尖端这个难以直接严密分析的区 域。这里利用 ! 积分的定义来求解其值。
/0 1234 5 "64， 7%8 9:4; （ <:=3>?@:4? AB +:CD342C3E F4;24::>24; ， GD243 H42I:>J2?K AB +2424; 34L M:CD4AEA;K， N6ODA6 1234;J6 ..)PPQ ， GD243） <51 =.402： B>3C?6>: =3>3@:?:>J； J?>:JJ 24?:4J2?K B3C?A>； & 24?:;>3E； %&’(’ 8>2,4(*,： G3EC6E3?2A4 @:?DAL AB B>3C?6>: =3>3@:?:>J 24 %&’(’ 2J 24?>AL6C:LR MD: >:J6E?J JDAS ?D3? ?D: & 5 24?:;>3E 34L J?>:JJ 24?:4J2?K B3C?A> 3>: I3>23TE: S2?D ?D: C>3CU A>2:4?3?2A4 34L C>3CU E:4;?D >:J=:C?2I:EKR MD:J: =>AI2L: I3E2L @:?DAL BA> ?D: C3EC6E3?2A4 AB CA@=E:V C>3CUJ 64L:> CA@=E:V EA3LJR
基于ANSYS的三维疲劳裂纹模拟
疲劳裂纹扩展问题在物理上是几何非连续现象，数学上一般不具有简单理论解，多数情况下只能依赖数值法，如边界元法，有限元法，扩展有限元法，无网格法等等。
截至目前，有限元法仍被认为最具通用性。
相当多商用有限元软件已集成裂纹参量计算的部分有关功能。
虽然利用有限元法研究疲劳裂纹问题的研究已十分丰富，但绝大多数研究往往集中于特定的三维裂纹（如纯I型）或是二维复合型裂纹，任意三维疲劳裂纹扩展模拟仍比较少。
这是因为一方面，材料试验数据有限，复合三维裂纹扩展理论还有待于完善，还没有一种裂纹扩展理论可以适用于所有不同材料不同载荷形式的裂纹扩展情况，在这方面目前存在的问题远多于已取得的共识；另一个重要原因是，复合裂纹引起的三维裂纹扩展路径为非平面裂纹面，形状复杂，在商用有限元软件中往往很难用简单图元表达该空间曲面，造成后续的有限元建模困难。
因此，用有限元法模拟任意三维裂纹扩展的关键技术之一在于建立非平片裂纹有限元模型，并保证裂纹前缘附近具有良好的网格质量。
为解决这一问题，近年出现了一些专用的裂纹前后处理软件（主要与ABQUS衔接），如ADAPCRACK3D、ZENCRACK、CRACK3D及康奈尔大学断裂力学小组开发的Franc3D （本身是边界元法前处理程序，但可作为有限元模型的前处理程序）等。
这些程序在建模思路、计算精度与计算效率方面，各有千秋。
对ANSYS来说，其自身带可编程的APDL语言，实现三维裂纹非平面扩展并非不可能。
当然，能够实现和计算结果准确是两回事。
这里涉及诸多方面，既有理论的约束，又有模型的简化。
要实现裂纹前缘的表达，具有一定的图形学基础是必须的。
一般采用样条曲线拟合法表示裂纹前缘。
为什么不采用差值法呢，这里涉及到计算误差及建模问题。
如果裂纹前缘局部呈凹性，采用较大的步长时，插值法所得曲线可能出现自相交，如图所示。
引起裂纹前缘局部凹陷主要由初始裂纹形状呈凹性或数值模拟精度等原因引起。
例如裂纹扩展计算中采用了较大的步长而发散或应力强度因子计算结果出现波动，拟合法可减少曲线的拐点数，较插值法稳定。
ansys断裂力学技巧
Ansys断裂力学裂纹和瑕疵在很多结构和零部件中会出现，有时会导致严重的后果。
断裂力学就是研究裂纹扩散问题的学科。
12.1 断裂力学的理解断裂力学就是解决结构在外载荷作用下，裂纹和瑕疵如何扩散的问题。
它包含裂纹扩散相应的解析预报和实验结果验证。
解析预报是通过断裂参数的计算得出的，如裂纹区域的应力强度因子，它可以用来评估裂纹的生长率。
最具典型的是，裂纹的长度随着一些循环载荷的每一次作用而增长，如飞机上机舱的增压-减压。
另外，环境的情况，如温度或光线的照射等，都会影响某些材料的断裂性能。
在研究中，断裂问题需重点研究的典型参数如下：●应力强度因子（K I, K II和K III），是断裂的三个基本形式。
●J-积分，是一种不受线路影响的线积分，用来测量裂纹端点的奇异应力和应变。
●能量释放率（G），它代表裂纹开始和终止处的能量的大小。
12.2 求解断裂力学问题求解断裂力学问题包括执行线弹性或弹塑性静态分析，以及使用专用的后处理命令或宏来计算需要的断裂参数。
此处分成两个部分来介绍：●裂纹区域的建模●计算断裂参数12.2.1裂纹区域的建模断裂模型中最重要的部分就是裂纹边界的部分。
在ansys中，在二维模型和三位模型中，分别将裂纹的边界看成是裂纹端点和裂纹前端。
如图12.1所示。
r是距离裂纹端点的长度。
裂裂纹面应该是重合纹端点处的应力和应变是奇异的，的，裂纹端点（或裂纹前端）附近的单元应该是二次的，即角点之间有中间节点。
这种单元被称为奇异单元。
12.2.1.1 二维断裂模型二维断裂模型的推荐单元类型是PLANE2，6节点的三角实体单元。
裂纹端点附近的单元的第一行是奇异的，如图12.2（a）所示。
前处理模块PREP7的命令(Main Menu> Preprocessor> Meshing> Size Cntrls> Concentrat KPs> Create)可以定义某关键点附近的单元划分的大小，在断裂模型中特别有用。
ansys 裂纹 坐标系
ansys 裂纹坐标系ANSYS裂纹坐标系引言：当研究材料的裂纹行为时，了解裂纹坐标系的概念是非常重要的。
裂纹坐标系是一种由裂纹定义的坐标系，用于描述裂纹的位置和方向。
在使用ANSYS软件进行裂纹有限元分析时，裂纹坐标系的正确使用是确保分析结果准确性和可靠性的关键。
本文将逐步介绍ANSYS裂纹坐标系的概念和使用方法。
第一部分：裂纹坐标系的概念和定义（300字）裂纹坐标系是一种通过选定的坐标轴和角度来描述裂纹位置和方向的数学框架。
裂纹坐标系的原点通常位于裂纹的起点，而坐标轴则沿着裂纹的主轴方向定向。
在裂纹坐标系中，通常使用极坐标或笛卡尔坐标来描述裂纹的位置和方向。
第二部分：裂纹坐标系的应用（500字）在裂纹有限元分析中，使用裂纹坐标系可以帮助我们确定载荷作用下裂纹尖端应力场的分布。
通过定义裂纹坐标系，可以将裂纹尖端的应力场分解为正应力、应变和剪应力分量，这有助于我们更好地理解裂纹的行为。
裂纹坐标系的应用主要包括：1. 确定裂纹尖端应力场的形状和大小。
2. 计算裂纹尖端的应力强度因子和应力强度因子范围。
3. 分析裂纹扩展速率和裂纹扩展路径。
4. 预测裂纹的疲劳寿命和断裂特性。
在ANSYS软件中，我们可以通过定义一个裂纹坐标系来进行裂纹有限元分析。
首先，需要在裂纹起点附近选择一个点作为裂纹坐标系的原点。
然后，可以沿着裂纹的主轴方向定义一个坐标轴，通常选择裂纹的延伸方向或裂纹尖端切向作为主轴。
在使用ANSYS时，可以使用命令行输入APDL 命令或者使用图形界面来定义裂纹坐标系。
第三部分：在ANSYS中定义裂纹坐标系（700字）在ANSYS中，我们可以使用APDL命令行或者图形界面来定义裂纹坐标系。
如果选择使用APDL命令行，可以使用命令'LCS'来定义裂纹坐标系。
该命令有多个选项可以用来定义裂纹坐标系的原点和坐标轴。
首先，使用'LCS'命令定义裂纹坐标系的原点。
例如，可以使用'LCS,0,0,0'命令将坐标系原点设置为全局坐标系中的点(0,0,0)。
利用ANSYS进行断裂分析
利用ANSYS进行断裂分析初次试做断裂分析，希望有这方面经验的高手能发表些经验之谈！这个模型由两种材料组成：表面镀层为铝,基底为钢。
目的是对表面镀层的剥离过程进行分析。
目前这个模型是个假想的简化模型，初步目标是实现剥离过程的模拟。
裂纹扩展是通过接触单元生死功能实现的。
基层和镀层由接触单元连接，然后通过断裂判断准则确定要杀死的失效的接触单元。
第一版(没有加断裂判断准则，强行逐个杀死界面接触单元)：fini/clear/filn,crack1/PREP7!*ET,1,PLANE182!*KEYOPT,1,1,2KEYOPT,1,3,1KEYOPT,1,4,0KEYOPT,1,6,0KEYOPT,1,10,0!*rect,0,100,0,100rect,0,100,100,110lesi,1,,,10lesi,2,,,10esha,2!*MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,210e3MPDATA,PRXY,1,,0.3MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,2,,70MPDATA,PRXY,2,,0.33amesh,1lesi,5,,,10lesi,6,,,2mat,2amesh,2lsel,s,,,3nsll,s,1cm,c1,nodelsel,s,,,5nsll,s,1cm,t1,nodensel,s,loc,xd,all,uxnsel,s,loc,yd,all,uyd,all,uxmp,mu,3,0/COM, CONTACT PAIR CREATION - START CM,_NODECM,NODECM,_ELEMCM,ELEMCM,_LINECM,LINECM,_AREACM,AREA/GSA V,cwz,gsav,,tempMP,MU,3,0MA T,3R,3REAL,3ET,2,169ET,3,172R,3,,,100,0.1,0,RMORE,,,1.0E20,0.0,1.0,RMORE,0.0,0,1.0,,1.0,0.5RMORE,0,0.5,1.0,0.0,KEYOPT,3,2,0KEYOPT,3,3,0KEYOPT,3,4,0KEYOPT,3,5,0KEYOPT,3,7,0KEYOPT,3,8,0KEYOPT,3,9,0KEYOPT,3,10,0KEYOPT,3,11,0KEYOPT,3,12,5! Generate the target surfaceNSEL,S,,,T1CM,_TARGET,NODETYPE,2ESLN,S,0ESURF,ALLCMSEL,S,_ELEMCM! Generate the contact surfaceNSEL,S,,,C1CM,_CONTACT,NODETYPE,3ESLN,S,0ESURF,ALLALLSELESEL,ALLESEL,S,TYPE,,2ESEL,A,TYPE,,3ESEL,R,REAL,,3/PSYMB,ESYS,1/PNUM,TYPE,1/NUM,1EPLOTESEL,ALLESEL,S,TYPE,,2ESEL,A,TYPE,,3ESEL,R,REAL,,3CMSEL,A,_NODECMCMDEL,_NODECMCMSEL,A,_ELEMCMCMDEL,_ELEMCMCMSEL,S,_LINECMCMDEL,_LINECMCMSEL,S,_AREACMCMDEL,_AREACM/GRES,cwz,gsavCMDEL,_TARGETCMDEL,_CONTACT/COM, CONTACT PAIR CREATION - END lsel,s,,,7nsll,s,1cm,s1,node!Gradient surface loadSFGRAD,PRES,0,X,0,-0.1,sf,all,pres,-0.1nsel,allesel,all!save/solutime,1deltim,1,1,1solve/post1plns,s,1anty,,resttime,1.1ekill,140solve/post1plns,s,1/soluanty,,resttime,1.2ekill,140ekill,139solve/post1plns,s,1/soluanty,,resttime,1.3ekill,140ekill,139ekill,138solve/post1plns,s,1/soluanty,,resttime,1.4ekill,140ekill,139ekill,138ekill,137solve/post1plns,s,1第二版(加了断裂自动判断准则)。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
基于ANSYS的三维贯穿裂纹的断裂参数计算
据一些工业化国家统计，因材料和结构的破坏所造成的损失占国民经济生产总值的8%-12%多。
破坏事故所造成的人员伤亡的损失更不可估量。
我国作为一个发展中国家，在这方面的情况比西方发达国家更严重。
因此无论是为了减少破坏事故的损失还是研发满足现代工业所需要的新材料，都要求对材料的破断过程有科学的、全面的、定量化的认识。
三维裂纹作为工程中常见的裂纹形式，早在六十年代初就有不少研究者开始研究，到现在已有大量的文献资料论及这一问题，出现了一些有特点的分析方法。
工程上常见的表面裂纹的断裂分析，由于其实质是三维问题，也几乎同时开始被人们所关注。
三维裂纹问题的危害极大，断裂造成了大量的灾难性事故发生，这使得断裂力学在机械工程、海洋工程、核工程，特别是今天的航空航天工程中受到更广泛的重视和深入研究。
因此对含三维裂纹结构断裂特性尤其对三维裂纹体的应力强度因子的研究有重要的现实意义。
本文使用ANSYS成功的计算了三维贯穿裂纹的应力强度因子，为计算三维裂纹提供了一种便捷方式。
1.模型的建立
图1 三维贯穿裂纹模型
本文三维裂纹模型长度为L，高度为H，宽度为W，裂纹半长为a，裂纹位于模型的中心部位。
几何参数见表1。
模型的为线弹性材料，其弹性模量为2.1E11Pa，泊松比为0.3。
模型的边界条件为：底端固定，顶端承受拉应力σ为2E6Pa。
表1 模型的几何参数
本文采用二维奇异单元PLANE183建立二维的裂纹模型，然后通过拉伸并使用三维奇异单元SOLID186来建立三维贯穿裂纹模型。
图2-图5给出了二维裂纹模型和三维裂纹模型。
在13.0中对应力强度因子的计算增加了一种计算方法即互动积分法（Interaction Integrals ），这种方法与计算J积分的主域积分法类似。
在二维问题进行面积分，在三维问题中进行体积分来获得应力强度因子。
这种方法与传统的位移扩展法相比精度高，需要的单元数少。
图2 二维裂纹模型图3 二维裂纹模型
裂纹尖端网格
图4 三维贯穿裂纹有限元模型图5 三维贯穿裂纹局
部图
2.结果与讨论
图6到图7给出了三维贯穿裂纹的应力等效云图，通过图7可知，SOLID186单元可以很好的模拟出裂纹尖端的应力场。
对于应力强度因子，本文计算图5中的三个位置即位置1（0.159，0.25 ，0.00），位置2（0.159，0.25 ，0.04）和位置3（0.159，0.25 ，0.1），并取8条积分线，计算结果列于表2中，通过表2可知，位置1和位置3的各条应力强度因子积分都相同，而位置2只有后三条的应力强度因子积分值与位置1和位置3相同。
出现这种差异的原因，图8给出了解释，主要是2位置的等效应力与位置1和位置3存在差异，最终造成断裂参数的计算差异。
这也说明了对于该模型，不能采用对称的计算。
由表2还可以推出，位置2的I型应力强度因子要大于位置1和位置3，当外载开始增加是，裂纹将从模型的内部开始扩展，这对于实际工程有这指导意义。
图6 三维贯穿裂纹的等效应力云图图7 三维贯穿裂纹的局部放
大图
表2 应力强度因子
图8 裂纹前缘的等效应力云图与位置关系图
3.结论
通过以上计算可以得出一下结论：
（1）互动积分法可以方便的计算三维裂纹问题；
（2）对于本文模型，内部裂纹的I型应力强度因子大于外表面的I型应力强度因子，这对于实际工程有这指导意义。