基于ANSYS的三维疲劳裂纹模拟

coolhand

疲劳裂纹扩展问题在物理上是几何非连续现象，数学上一般不具有简单理论解，多数情况下只能依赖数值法，如边界元法，有限元法
，扩展有限元法，无网格法等等。截至目前，有限元法仍被认为最具通用性。相当多商用有限元软件已集成裂纹参量计算的部分有关功能。虽然利用有限元法研究疲劳裂纹问题的研究已十分丰富，但绝大多数研究往往集中于特定的三维裂纹（如纯I型）或是二维复合型裂纹，任意三维疲劳裂纹扩展模拟仍比较少。这是因为一方面，材料试验数据有限，复合三维裂纹扩展理论还有待于完善，还没有一种裂纹扩展理论可以适用于所有不同材料不同载荷形式的裂纹扩展情况，在这方面目前存在的问题远多于已取得的共识；另一个重要原因是，复合裂纹引起的三维裂纹扩展路径为非平面裂纹面，形状复杂，在商用有限元软件中往往很难用简单图元表达该空间曲面，造成后续的有限元建模困难。因此，用有限元法模拟任意三维裂纹扩展的关键技术之一在于建立非平片裂纹有限元模型，并保证裂纹前缘附近具有良好的网格质量。为解决这一问题，近年出现了一些专用的裂纹前后处理软件（主要与ABQUS衔接），如ADAPCRACK3D、ZENCRACK、CRACK3D及康奈尔大学断裂力学小组开发的Franc3D（本身是边界元法前处理程序，但可作为有限元模型的前处理程序）等。这些程序在建模思路、计算精度与计算效率方面，各有千秋。
   对ANSYS来说，其自身带可编程的APDL语言，实现三维裂纹非平面扩展并非不可能。当然，能够实现和计算结果准确是两回事。这里涉及诸多方面，既有理论的约束，又有模型的简化。

coolhand

要实现裂纹前缘的表达，具有一定的图形学基础是必须的。一般采用样条曲线拟合法表示裂纹前缘。为什么不采用差值法呢，这里涉及到计算误差及建模问题。如果裂纹前缘局部呈凹性，采用较大的步长时，插值法所得曲线可能出现自相交，如图所示。引起裂纹前缘局部凹陷主要由初始裂纹形状呈凹性或数值模拟精度等原因引起。例如裂纹扩展计算中采用了较大的步长而发散或应力强度因子计算结果出现波动，拟合法可减少曲线的拐点数，较插值法稳定。

coolhand

下面给出一种三维裂纹非平面扩展的建模方法。整个裂纹面用一组Coons面进行逼近，每个Coons面均由两条相邻的扩展轨迹线形成，如图2所示。将该裂纹面沿着法向与反法向进行增厚，形成两组实体。为构造几何不连续性，两组体中间并不共面。为保证裂纹前缘附近具有良好的网格质量，在裂纹面的自然延伸方向，再形成两个实体，如图3所示。为简洁，这两组实体称为裂纹体。图4显示了该网格划分的思想，即裂纹体被结构包含，裂纹面又包含在裂纹体中。如果一个含裂纹结构按上述方法进行实体剖分，由于裂纹体是比较规则的实体，因而裂纹前缘附近的网格比较容易控制。图5(a) 是整个结构的网格划分示意图；图5(b)为裂纹体网格，图5(c)为裂尖附近典型截面网格（使用了1/4边中节点奇异元）。

coolhand

用ANSYS实现的一个简单算例（算例来自 Buchholz F.-G., Chergui A. Richard H.A. Fracture analyses and experimental results of crack growth under general mixed mode loading conditions[J]. Engng. Fracture Mech, 2007?)
图6为几何参数 图7为有限元模型 图8为扩展后形状  图9为文献中实验结果
欢迎交流，不怕拍砖

coolhand

从试验结果看，其裂纹面起始扩展路径十分粗糙，裂纹面由非共面的多重裂纹源组成，反映出复合型三维裂纹扩展的复杂性。显然，目前通过有限元分析还难以模拟出这种裂纹面前缘的多源性与扩展方向存在一定的随机性，不能反映其细节特征。

jazzydzj

有米有个简单通用的例子把裂纹扩展的模拟下，比如说一块板或者什么的

郑俊雅策划

一步步超越自己！

llijianguo

看看，学习一下！谢啦，楼主！

yunyun5258

经典呢，多谢楼主，辛苦了

wangjibo1208

额，对俺来说还是有点高深~~~

qilin307609

由于毕业在即，本人寻求 朋友们谁能做螺旋压缩弹簧的疲劳裂纹萌生和扩展的计算机模拟
（注：用ANSYS二次开发做） 谁做出来了，我给5000元人民币,联系QQ524125655

coolhand

11# qilin307609

没有涉及过弹簧的裂纹扩展计算。
实际的结构其失效机理为疲劳失效？起裂位置？裂纹为张开型还是扭转撕裂型？
如果模型不是很复杂，可传个简化模型上来帮你免费一试，也帮我考核一下我的程序能不能计算。

hfwh888

用数值做这个会有说服力么？

coolhand

用数值做这个会有说服力么？
hfwh888 发表于 2011-3-2 16:33

若力学模型不合理、基本材料数据不充分，难言可信度。
断裂力学很多方面理论还不够完善，数值成功的例子比较少，大多还以实验为主。
但作为研究一种辅助手段，数值分析还是很有必要。实验存在代价昂贵、周期长等不足。

CHENLONG

楼主的这个例子很经典，学习了！

lcm315

向楼上的大家学习啊，我也想请教一个问题，在模拟结构损伤时，1、利用降低刚度的方法；2、模拟实际裂纹，这两者的差距有多大呢

dreamfly9999

楼主强大，谢谢

dreamfly9999

楼主，留个联系方式，咨询点问题，谢了。516721764 2# coolhand

歌菲

楼主很厉害呀，我用重合网格法（s-version FEM）有模拟过你这个例子。

云弓

受教了！刚接触裂纹扩展问题，看楼主的文章还是有些太深奥！！！